资源描述：
《空间桁架结构程序设计(fortran)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、空间桁架静力分析空间桁架静力分析程序及算例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、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)支
3、座反力数组FL(NF*I-2)存放受约束的I节点X方向的反力FL(NF*I-1)存放受约束的I节点Y方向的反力FL(NF*I)存放受约束的I节点Z方向的反力中间变量NPF=NF*NP二维总刚度矩阵的最大行数-16-空间桁架静力分析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结构一维总
4、刚度矩阵的总容量CKK(NWK)结构一维总刚度矩阵NN结构矩阵方程的方程总数（去掉约束）NNMNNM=NN+1V(NN)已知节点荷载列阵数组，回代完成后为存放结构位移PP(NPF)所有节点荷载列阵数组2、空间桁架结构有限元分析程序源代码!主程序（读入文件，调用总计算程序，输出结果）CHARACTERIDFUT*20,OUTFUT*20WRITE(*,*)'InputDataFilename:'READ(*,*)IDFUTOPEN(11,FILE=IDFUT,STATUS='OLD')WRITE(*,*)'OutputFilename:'READ(*,*)O
5、UTFUTOPEN(12,FILE=OUTFUT,STATUS='UNKNOWN')WRITE(12,*)'*****************************************'WRITE(12,*)'*ProgramforAnalysisofSpaceTrusses*'WRITE(12,*)'*SchoolofCivilEngineeringCSU*'WRITE(12,*)'*2012.6.25DesignedByMuZhaoxiang*'WRITE(12,*)'***************************************
6、**'WRITE(12,*)''WRITE(12,*)'*****************************************'WRITE(12,*)'*************TheInputData****************'WRITE(12,*)'*****************************************'WRITE(12,100)READ(11,*)NF,NP,NE,NM,NR,NCF,NDWRITE(12,110)NF,NP,NE,NM,NR,NCF,ND100FORMAT(6X,'TheGeneralI
7、nformation'/2X,'NF',5X,'NP',5X,'NE',5X,'NM',5X,'NR',&5X,'NCF',5X,'ND')110FORMAT(2X,I2,6I7)NPF=NF*NPNDF=ND*NFCALLANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF)END!********************************************************************!总计算程序-16-空间桁架静力分析SUBROUTINEANALYSE(NF,NP,NE,NM,NR,NCF,ND,N
8、PF,NDF)DIMENSIONX(NP),Y(NP),Z(NP)
