平面桁架计算程序.pdf

! PTA 说明 : 1,数据输入文件为 DATAIN.TXT 中,数据依次为 (NE,NJ,NR,NP,NP); (X,Y)*NJ ; ! (IJ1,IJ2,A,EI)*NE;(JR1,JR2,JR3,JR4)*NR;(PJ1,PJ2,PJ3)*NP;(PF1,PF2,PF3,PF4)*NF ! 2,NE 单元总数； NJ 节点总数； NR 约束总数； NP节点荷载总数； NF非节点荷载数； X,Y 节点坐标； ! IJ(NE,1),IJ(NE,2) 单元期终点； A 单元面积； ZI 截面惯性矩； JR(NR,1) 约束 结点号 ; ! JR(NE,2:4) 横向竖向转动约束 (1); PJ(NP,1:3) 分别为结点号、荷载类型、荷载 值； ! PF(NF,1:4) 分别为单元号、荷载类型，荷载值、荷载据起点距离 !主程序 PROGRAM Plane_Truss_Analysis implicit none integer NE,NJ,NP,NF,NR,N real E integer,allocatable,DIMENSION(:,:)::IJ,JR real,allocatable,DIMENSION(:)::A,P,X,Y !动态数组定义 real,allocatable,DIMENSION(:,:)::PJ,PF,TK OPEN(1,FILE=datain.TXT,STATUS=OLD) !打开文件大 datain.txt 并存到标号 1 中 OPEN(2,FILE=dataout.TXT,STATUS=NEW) !建立文件 dataout存储计算结果 READ(1,*)NE,NJ,NR,NP,NF,E N=NJ*2 !位移总数 allocate(X(1:NJ),Y(1:NJ),IJ(1:NE,2),JR(1:NR,4),A(1:NE),PJ(1:NP,3),PF(1:NF,4),TK(1:N,1:N),P(1:N) ) WRITE(2,10)NE,NJ,NR,NP,NF,E ! 打印表头 10 FORMAT(/1X, *********** 平面桁架力计算 PTA***********//4X, 单元数 NE=,I2,12X, 结点数 NJ=,I2,14X, 支座数 NR=,I2,/4X, 结点荷载数 NP=,I2,8X, 非节点荷载数 NF=,I2,8X, 弹性模量 E=,E12.4) CALL INPUT(NE,NJ,NR,NP,NF,X,Y,IJ,A,JR,PJ,PF) !数据输入 CALL TSM(NE,NJ,E,X,Y,IJ,A,TK,N) ! 形成结构原始刚度矩阵 CALL JLP(NE,NJ,NP,NF,X,Y,IJ,PJ,PF,P,N) !形成结构综合节点荷载阵列 CALL ISC(NR,JR,TK,P,N) ! 引入结构约束条件 CALL GAUSS(TK,P,N) ! 高斯消去法计算结构坐标下的节点位移 CALL MVN(NE,NJ,NF,E,X,Y,IJ,A,PF,P,N) ! 计算单元杆端内力 CLOSE(1) CLOSE(2) deallocate(X,Y,IJ,JR,A,P,PJ,PF,TK) END PROGRAM Plane_Truss_Analysis !原始数据输入 SUBROUTINE INPUT(NE,NJ,NR,NP,NF,X,Y,IJ,A,JR,PJ,PF) DIMENSION X(NJ),Y(NJ),IJ(NE,2),A(NE),JR(NR,4),PJ(NP,3),PF(NF,4) READ(1,*) (X(I),Y(I),I=1,NJ) !从 1 中顺序读取节点坐标 READ(1,*