空间桁架结构程序设计(Fortran).docVIP

空间桁架静力分析程序及算例 1、变量及数组说明 输入数据 控制数据 NF 单个节点的自由度数 NP 结构离散节点的总数 NE 结构离散单元的总数 NM 结构中单元不同的特征数类的总数 NR 结构受约束节点的总数 NCF 结构受外荷载作用的节点总数 ND 一个单元的节点总数 几何数据 X(NP) 节点X坐标数组 Y(NP) 节点Y坐标数组 Z(NP) 节点Z坐标数组 ME(ND,NE) 单元节点信息存储矩阵 ME(1,NE)存储杆件始端节点号 ME(2,NE)储存杆件末端节点号 RR(2,NR) 结构约束信息矩阵 RR(1,NR)存放受有约束的节点号 RR(2,NR)存放节点位移约束情况 单元特征数据 AE(2,IN) 单元特征数类数组 AE(1,IN)单元的弹性模量 AE(2,IN)单元的横截面面积 NAE(NE) 单元特征类信息存储数组 荷载数据 PF(4,NCF) 外荷载信息数组 PF(1,NCF)存放外荷载作用的节点号 PF(2,NCF)存放X方向的外荷载 PF(3,NCF)存放Y方向的外荷载 PF(4,NCF)存放Z方向的外荷载 位移 DIST(NPF) 节点位移数组 DIST(NF*I-2)存放I节点X方向的位移 DIST(NF*I-1)存放I节点Y方向的位移 DIST(NF*I) 存放I节点Z方向的位移 力 SG(NE) 单元内力数组 SM(NE) 单元截面应力数组 FL(NF*NR) 支座反力数组 FL(NF*I-2)存放受约束的I节点X方向的反力 FL(NF*I-1)存放受约束的I节点Y方向的反力 FL(NF*I)存放受约束的I节点Z方向的反力 NPF=NF*NP 二维总刚度矩阵的最大行数 NDF=ND*NF 一个单元的自由度总数（2*3=6） IN 单元特征类总数 AKE(2,2) 单元在局部坐标系中的刚度局矩阵 BL 杆件单元长度 T(2,6) 坐标转换矩阵 TAK(6,6) 单元在总体坐标系中的刚度矩阵 IT(NF,NP) 节点联系数组 LMT(NDF,NE) 单元联系数组 MAXA(NPF) 结构二维总刚度矩阵主对角元地址数组 NWK 结构一维总刚度矩阵的总容量 CKK(NWK) 结构一维总刚度矩阵 NN 结构矩阵方程的方程总数（去掉约束） NNM NNM=NN+1 V(NN) 已知节点荷载列阵数组，回代完成后为存放结构位移 PP(NPF) 所有节点荷载列阵数组 CHARACTER IDFUT*20,OUTFUT*20 WRITE(*,*) Input Data File name: READ (*,*)IDFUT OPEN (11,FILE=IDFUT,STATUS=OLD) WRITE(*,*) Output File name: READ (*,*)OUTFUT OPEN(12,FILE=OUTFUT,STATUS=UNKNOWN) WRITE(12,*)***************************************** WRITE(12,*)* Program for Analysis of Space Trusses * WRITE(12,*)* School of Civil Engineering CSU * WRITE(12,*)* 2012.6.25 Designed By MuZhaoxiang * WRITE(12,*)***************************************** WRITE(12,*) WRITE(12,*)***************************************** WRITE(12,*)*************The Input Data**************** WRITE(12,*)***************************************** WRITE(12,100) READ(11,*)NF,NP,NE,NM,NR,NCF,ND WRITE(12,110)NF,NP,NE,NM,NR,NCF,ND 100 FORMAT(6X,The General Information/2X,NF,5X,NP,5X,NE,5X,NM,5X,NR, 5X,NCF,5X,ND) 110 FORMAT(2X,I2,6