空间刚架计算程序及说明.docVIP

设计意图及总体框架： 大多数空间刚架都属于超静定结构，因此难于用一般的结构力学方法进行计算，如用矩阵位移法则由于阶数太大难于手算，但如果结合有限元方法编写程序进行计算，则会有很大的普遍性，且方便快捷。 本程序能计算只承受任意的非节点集中力、均布荷载和集中力偶矩的所有刚架（包括平面刚架、平面板架）的节点内力和位移，但不适用有支座沉降和温度变化情况，计算出节点力之后，可以很方便地画出结构内力图，从而整个问题都得到求解。 程序设计的整个过程可以参阅相关书籍，空间刚架的程序设计和平面刚架的程序设计的最大区别在于空间刚架计算的坐标变换矩阵较复杂和由于维数和阶数的增加所带来的子程序之间衔接的复杂性，除此之外，一些平面刚架程序的子函数可以直接用到空间刚架程序中来。 本程序一共包括10子函数：局部坐标下单元刚度计算子程序KEP()、坐标变换矩阵计算子程序CR()、固端反力计算子程序FIXF()、单元定位向量计算子程序（）、指示数组计算子程序LD()、约束处理子程序FCC()、解线性方程组子程序DECOM()、矩阵赋零子程序ZERO()、矩阵转置子程序TRAN()和矩阵相乘子程序DOT()。具体算法见源程序清单，输入输出都用文件进行。 注：1,为了方便求坐标变换矩阵，可以在建立局部坐标系时使y’//XOY,而这一点我们总是可以做到的 2,在上面的坐标系下，可以将刚架单元分为竖直单元和一般单元，分别求坐标转换矩阵，具体做法见程序。 3,单元刚度矩阵可以用平面刚架和平面板架的结果叠加。 变量及符号： NN——节点总数， NM——材料类型数 NA——截面类型数 NE——单元总数， NC——受约束自由度数 NF——节点自由度数 NP2——非节点荷载数 ND——单元节点数 NFD——单元自由度数 N——结构自由度总数 NT——一维数组A的容量 X(NN)——X坐标数组， Y(NN)——Y坐标数组 Z(NN)——Z坐标数组 MEA(NE,4)——单元信息数组，其中: MEA(I,1)——I单元i节点编号 MEA(I,2)——I单元j节点编号 MEA(I,3)——I单元截面编号 MEA(I,4)——I单元材料编号 AEU(NM,3)——材料信息数组，其中: AEU(NM1,1)——弹性模量EO(E) AEU(NM1,2)——剪切模量GO(G) AEU(NM1,3)——泊松比UO(u) AAI(NA,4)——截面信息数组，其中： AAI(NA1,1)——NA1截面面积A AAI(NA1,2)——绕Z轴惯性矩IZ AAI(NA1,3)——绕Y轴惯性矩IY AAI(NA1,4)——绕X轴惯性矩IP P2(NP2,4)——非节点荷载数组，其中： P2(I,1)——作用节点号， P2(I,2)——荷载类型码（见程序） P2(I,3)——作用位置参数（见程序） P2(I,4)——荷载大小 JC(NC,2)——约束信息数组，其中： JC(I,1)——约束节点号 JC(I,2)——约束代码取1，2，3，4，5，6 代码 1:U=0, 2:V=0, 3:W=0, 4:SITA-X=0, 5:SITA-Y=0, 6:SITA-Z=0 ANG(NE)——单元角数组，即Y与Y’所成的角度，用以求坐标变换矩阵 AKEP(12,12)——局部坐标系下单元刚度矩阵 AKE(12,12)——整体坐标系下单元刚度矩阵 R(12,12)——坐标变换矩阵 RT(12,12)——R(12,12)的转置矩阵 POP(12)——局部坐标系下固端力向量 PO(12)——整体坐标系下固端力向量 A(NT)——一维变带宽存储总刚度矩阵一维数组 P(N)——节点荷载向量 UPE(12)——局部坐标系下节点位移 IS(12)——定位向量数组 LD(600)——指示数组 源程序清单： C 3D-FRAME STRUCTURE ANALYAIS PROGRAM MAIN /*主程序开始