资源描述：
《平面桁架杆单元fortran程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、PROGRAMBAR3DDIMENSIONMEA(100,4),JZ(50,2),AAI(10),AEU(10,2),R(2,6),RT(6,2)DIMENSIONAKEP(2,2),AKE(6,6),P(300),X(100),Y(100),Z(100)DIMENSIONLD(300),IS(6),A(50000),UPE(2),PPE(2),UE(6),RA(6,2)OPEN(1,FILE='FNAME1_t.txt')OPEN(2,FILE='FNAME2_t.txt')READ(1,*)NN,NE,NM,NA,NCwrite
2、(*,*)NN,NE,NM,NA,NCNF=3ND=2NFD=NF*NDN=NN*NFDO50I=1,NN50READ(1,*)K,X(I),Y(I),Z(I),(P(NF*(I-1)+J),J=1,NF)READ(1,*)((JZ(I,J),J=1,2),I=1,NC)DO100I=1,NE100READ(1,*)IE,(MEA(I,J),J=1,4)READ(1,*)(AAI(I),I=1,NA)READ(1,*)((AEU(I,J),J=1,2),I=1,NM)CLOSE(1)WRITE(2,460)NN,NE,NM,NA,N
3、CWRITE(2,465)(I,X(I),Y(I),Z(I),P(3*I-2),P(3*I-1),P(3*I),I=1,NN)WRITE(2,470)((JZ(I,J),J=1,2),I=1,NC)write(2,475)(I,(Mea(i,j),j=1,4),i=1,Ne)WRITE(2,480)(I,AAI(I),I=1,NA)WRITE(2,485)(I,(AEU(I,J),J=1,2),I=1,NM)400FORMAT(5I5)405FORMAT(I5,6F10.0)410FORMAT(2I5)415FORMAT(5I5)4
4、20FORMAT(F10.0)425FORMAT(2F10.0)460FORMAT(/5X,'THEINPUTOFNN,NE,NM,NA,NC'//(5X,5I5))465FORMAT(/5X,'THEINPUTOFX,Y,Z,P'//(5X,I5,6F10.2))470FORMAT(/5X,'THEINPUTOFJZ'//(5X,I5,5X,I5))475FORMAT(/5X,'THEINPUTOFMEA'//(5X,I5,5X,4I5))480FORMAT(/5X,'THEINPUTOFAAI'//(5X,I5,5X,F20.6
5、))485FORMAT(/5X,'THEINPUTOFAEU'//(5X,I5,5X,2E15.6))CALLFLD(NN,NE,MEA,NF,ND,N,NT,LD)DO500I=1,NT500A(I)=0.0DO600IE=1,NECALLKEP(NN,NA,IE,MEA,AEU,AAI,X,Y,Z,AKEP)!CALLCR(NN,IE,MEA,X,Y,Z,R)CALLTRAN(2,6,R,RT)CALLDOT(6,2,2,RT,AKEP,RA)CALLDOT(6,2,6,RA,R,AKE)CALLFIS(IE,MEA,NF,ND
6、,NFD,IS)DO560I=1,NFDDO560J=1,NFDIF(IS(I)-IS(J))560,520,520520NI=IS(I)IJ=LD(NI)-NI+IS(J)A(ij)=A(ij)+AKE(I,J)560CONTINUE600CONTINUECALLFCC(NC,N,NT,NF,JZ,LD,A)CALLDECOM(N,NT,A,P,LD)WRITE(2,850)(I,P(3*I-2),P(3*I-1),P(3*I),I=1,NN)WRITE(2,860)DO800IE=1,NECALLKEP(NN,NA,IE,MEA
7、,AEU,AAI,X,Y,Z,AKEP)CALLCR(NN,IE,MEA,X,Y,Z,R)CALLFIS(IE,MEA,NF,ND,NFD,IS)DO750I=1,NFDNI=IS(I)UE(I)=P(NI)750CONTINUECALLDOT(2,NFD,1,R,UE,UPE)CALLDOT(2,2,1,AKEP,UPE,PPE)NA1=MEA(IE,3)AO=AAI(NA1)STRESS=PPE(2)/AOWRITE(2,880)IE,PPE(2),STRESS800CONTINUE850FORMAT(//25X,'NODALD
8、ISPLACEMENT'//5X,'NODENo.',12X,'U',12X,'V',12X,'W'/(7X,I5,3E15.6))860FORMAT(//15X,'AXIALFORCEANDSTRESSOFELEMENT'/11X,
