连续梁静力分析程序设计.docxVIP

连续梁静力分析程序设计 班级 土木二班 学号 1083310424 姓名 卢显滨 一、连续梁静力分析程序总体设计 1.程序的功能与算法设计 1.1连续梁计算程序的功能设计 计算程序的功能要求： （1）可计算均质同一材料的多跨连续梁 （2）可进行5种荷载作用下的受力分析 （3）可计算任意指定截面的内力 （4）可求任意指定截面的弯矩、剪力影响线 （5）可求任意指定支座的反力影响线 （6）可计算均布活荷载作用下的弯矩和剪力包络图 （7）尽量减少初始数据的输入 （8）尽可能提高程序质量 （9）尽量节省内存单元 1.2连续梁计算程序的算法设计 为实现上述功能和要求，本程序采用如下算法： （1）采用矩阵位移法分析，以便能计算大多数多跨连续梁 （2）总刚度矩阵采用变半带宽一维存储方法，以便节省内存 （3）用置0置1法处理已知支座位???情况 （4）结点信息、单元信息、荷载信息、支承情况分类输入，相同信息可自动生成，减少数据输入量 （5）求解位移时，考虑系数矩阵是三对角线阵，用高斯消元法较简单、便捷 （6）求内力包络图时，先将单元按奇偶数标记，欲求某跨跨中最大正弯矩时，应在该跨布置活荷载，然后向两侧隔跨布置；欲求某支座截面最大负弯矩时，应在该支座相邻两跨布置活荷载，然后向两侧隔跨布置。 2.主程序框图设计 由矩阵位移法可知，连续梁内力计算的矩阵位移法一般步骤为： （1）离散化（数据化），划分单元，编结点码、位移码，为结构分析准备数据，即输入或生成计算用初始数据。 （2）单元分析，计算单元刚度矩阵，按单元定位向量组装成总体刚度矩阵 （3）当有已知支座位移时进行支承条件处理 （4）求解以结点位移为未知数的线性方程组 （5）计算各单元的杆端力 （6）计算指定截面内力 （7）求指定截面内力影响线 （8）求指定支座反力影响线 （9）求内力包络图 综合程序功能要求和矩阵位移法的计算步骤，连续梁静力分析主程序流程图如图1所示。 说明部分 输入数据 集装整体刚度矩阵 开始 结束 生成综合荷载列阵azhe 支承条件处理 解线性方程组 输出结点位移 计算并输出杆端力 计算指定截面内力 求影响线 求包络图 图1 程序流程图 由流程图可知，主程序包含13个子程序：输入数据（SRSJ）、单刚（DG）、集装整体刚度矩阵（JCZK）、等效结点荷载（DJH）、生成荷载列阵（XCP）、支承条件处理（ZCCL）、方程求解（FCQJ）、输出位移（SCWY）、计算杆端力（GDL）、计算截面内力（JMNI）、求影响线（SCGS）、求包络图（LXLB）和求包络图时的荷载信息（HZXX） 3.子程序功能设计 1.初始数据输入子程序SRSJ 根据输入的单元数、结点数、结点荷载数、非结点荷载数，分别输入各单元的杆长、惯性矩、结点荷载、非结点荷载，存贮在GC、GX、PJ、PF数组中，为结构分析准备数据。求截面内力、影响线、包络图时还需要截面信息，存在JM数组中。 2.集装整体刚度矩阵子程序JCZK 根据各单元信息，调用计算单元刚度矩阵的子程序DG，按单元定位向量组装，用ZK1存储整体刚度阵的主对角线元素，副对角线元素存在ZK2中。 3.生成总荷载列阵子程序XCP 作用在连续梁上的荷载分为直接结点荷载PJ和非节点荷载数PF，若有非结点荷载，调用DJH子程序生成等效结点荷载F0，根据定位向量生成综合荷载列阵P。 4.支承条件处理子程序ZCCL 用数组NZ存储连续梁两端支承情况，有零位移输0,否则输1。用置0置1法处理刚度矩阵K和荷载列阵P。若NZ(1)=0，则ZK1(1)=1.0，ZK2(1)=0.0，P(1)=0.0；若 NZ(2)=0，则ZK1(NJ)=1.0，ZK2(NJ-1)=0.0，P(NJ)=0.0（NJ为结点数）。 5.方程求解子程序FCQJ 因刚度方程中的系数矩阵是三对角线的稀疏矩阵，用高斯消元法较简捷。即通过方程的相加减，让系数矩阵的下三角元素化为零，可得最末一个位移X(NJ)，再回代到上一个方程中求出X(NJ-1)，一直迭代到求出所有位移为止。 6.计算并输出杆端力子程序GDL 由矩阵位移法原理，计算杆端力时要用到结点位移、单元刚度矩阵和等效结点荷载，需调用DG、DJH子程序。 7.计算截面内力子程序JMNI 单元杆端力和荷载作用在简支梁上，叠加得到截面内力，内力值存在LJM中。 内力应分截面左值和右值，LJM(1,1) 、LJM(1,2) 、LJM(1,3) 、LJM(1,4)分别为左弯矩、左剪力、右弯矩、右剪力。在求支座反力影响线时，用支座右端剪力减去左端剪力，并特殊考虑荷载作用在所求截面位置时左右截面内力值的变化。 8.求包络图时的荷载信息子程序HZXX 活荷载分布情况一共有NJ种（结点数），第一种为奇数跨布置活荷载，第二种为偶数跨