空间桁架结构程序设计

空间桁架静力分析程序及算例1、变量及数组说明

2、空间桁架结构有限元分析程序源代码

!主程序（读入文件，调用总计算程序，输出结果）

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,*)'* 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,6I7)

NPF=NF*NP

NDF=ND*NF

CALL ANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF)

END

!********************************************************************

!总计算程序

SUBROUTINE ANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF)

DIMENSION X(NP),Y(NP),Z(NP),MM(NE),ME(ND,NE),IT(NF,NP),RR(ND,NR), NAE(NE),&

AE(1,2),PF(4,NCF),LMT(NDF,NE),MAXA(NPF),CKK(1000),V(NPF),DIST(NPF),&

PP(NPF),FF(NPF),SG(NE),SM(NE)

READ(11,*)(X(I),Y(I),Z(I),I=1,NP)

READ(11,*)(MM(I),ME(1,I),ME(2,I),NAE(I),I=1,NE)

READ(11,*)(RR(1,J),RR(2,J),J=1,NR)

READ(11,*)(AE(1,J),J=1,2)

WRITE(12,120)

WRITE(12,121)(I,X(I),Y(I),Z(I),I=1,NP)

WRITE(12,130)

WRITE(12,131)(MM(I),ME(1,I),ME(2,I),NAE(I),I=1,NE)

WRITE(12,140)

WRITE(12,141)(INT(RR(1,J)),RR(2,J),J=1,NR)

WRITE(12,150)

WRITE(12,151)(AE(1,J),J=1,2)

IF(NCF/=0)THEN

READ(11,*)((PF(I,J),I=1,4),J=1,NCF)

WRITE(12,160)

WRITE(12,161)(INT(PF(1,J)),PF(2,J),PF(3,J),PF(4,J),J=1,NCF)

ENDIF

120 FORMAT(/6X,'The Information of Joints'/2x,'Joint',5X,'X',5X,'Y',5X,'Z')

121 FORMAT(1X,I4,

130 FORMAT(/6X,'The Information of Members'/2x,'Member',2X,'START',4X,'END',6X,'NAE')

131 FORMAT(1X,I4,3I8)

140 FORMAT(/6X,'The Information of SUPPORTS'/2x,'Joint',5X,'S')

141 FORMAT(1X,I4,

150 FORMAT(/6X,'The Information of Sections'/4x,'E0',8X,'A0')

151 FORMAT(1X,,

160 FORMAT(/6X,'The Loading at Joints'/2x,'Joint',5X,'FX',5X,'FY',7X,'FZ')

161 FORMAT(1X,I4,

CALL FLMT(NP,NE,NN,NNM,NR,RR,ND,NF,NDF,ME,IT,LMT)

CALL FMAXA(NNM,NE,LMT,MAXA,NWK,NPF,NDF)

CALL LP(V,PP,IT,PF,NN,NCF,NF,NP,NPF)

CALL CONKB(NP,NE,NM,NWK,ME,X,Y,Z,AE,NAE,LMT,MAXA,CKK,NNM)

ISH=1

CALL LDLT(CKK,MAXA,NN,ISH,IOUT,NWK,NNM)

CALL REBACK(CKK,V,MAXA,NN,NWK,NNM)

CALL DISPLS(NP,NE,NPF,NM,NN,IT,V,DIST,AE,NAE,X,Y,Z,PP,FF,SG,SM,ME,NR,RR,NF)

END

!********************************************************************

!矩阵转置子程序

SUBROUTINE MAT(M,N,A,B)

DIMENSION A(M,N),B(N,M)

DO I=1,M

DO J=1,N