梁静力试验报告.docx

梁静力试验报告 　　连续梁静力分析程序设计　　班级土木二班学号姓名卢显滨　　一、连续梁静力分析程序总体设计　　1.程序的功能与算法设计　　连续梁计算程序的功能设计　　计算程序的功能要求：　　可计算均质同一材料的多跨连续梁可进行5种荷载作用下的受力分析可计算任意指定截面的内力　　可求任意指定截面的弯矩、剪力影响线可求任意指定支座的反力影响线　　可计算均布活荷载作用下的弯矩和剪力包络图尽量减少初始数据的输入尽可能提高程序质量尽量节省内存单元连续梁计算程序的算法设计　　为实现上述功能和要求，本程序采用如下算法：　　采用矩阵位移法分析，以便能计算大多数多跨连续梁总刚度矩阵采用变半带宽一维存储方法，以便节省内存用置0置1法处理已知支座位移情况　　结点信息、单元信息、荷载信息、支承情况分类输入，相同信息可自动生成，减少数据输入量　　求解位移时，考虑系数矩阵是三对角线阵，用高斯消元法较简单、便捷　　求内力包络图时，先将单元按奇偶数标记，欲求某跨跨中最大正弯矩时，应在该跨布置活荷载，然后向两侧隔跨布置；欲求某支座截面最大负弯矩时，应在该支座相邻两跨布置活荷载，然后向两侧隔跨布置。2.主程序框图设计　　由矩阵位移法可知，连续梁内力计算的矩阵位移法一般步骤为：离散化，划分单元，编结点码、位移码，为结构分析准备数据，即输入或生成计算用初始数据。　　单元分析，计算单元刚度矩阵，按单元定位向量组装成总体刚度矩阵当有已知支座位移时进行支承条件处理求解以结点位移为未知数的线性方程组计算各单元的杆端力计算指定截面内力　　求指定截面内力影响线求指定支座反力影响线求内力包络图　　综合程序功能要求和矩阵位移法的计算步骤，连续梁静力分析主程序流程图如图1所示。　　图1程序流程图　　由流程图可知，主程序包含13个子程序：输入数据、单刚、集装整体刚度矩阵、等效结点荷载、生成荷载列阵、支承条件处理、方程求解、输出位移、计算杆端力、计算截面内力、求影响线、求包络图和求包络图时的荷载信息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)=，ZK2(1)=，P(1)=；若NZ(2)=0，则ZK1(NJ)=，ZK2(NJ-1)=，P(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种，第一种为奇数跨布置活荷载，第二种为偶数跨布置活荷载，其它情况是中间某支座相邻单元布置活荷载，然后隔跨布置。　　比较各种荷载情况下截面内力的大小，求出最大值和最小值。8.输出格式子程序　　求影响线的输出子程序SCGS，按任务分为弯矩、剪力影响线和支反力影响线，输出时避免输出相同截面、支座的影响线数据，尽可能节省内存。求包络图的输出子程序LXLB，将最大值小于0和最小值大于0的截面内力值，都作0处理。　　各子程序的关系见图2。　　图2子程序调用关系　　二、程序的验算与应用　　1.计算指定截面内力　　算例1：如图3，连续梁受集中力、集中力偶、