有限元实验.docVIP

实验二 平面桁架程序设计 平面桁架内力和位移计算的总框图 平面桁架静力分析源程序（PTSAP.FOR） C ANALYSIS PROGRAM FOR PLANE TRUSS DIMENSION JE(2,100),JEA(100),JN(2,100),M(4),EA(2,25), JPJ(50),PJ(50),P(200),FE(4),T(4,4),TL(20),AKE(4,4),F(4), AK(200,100),X(100),Y(100),FF(100) OPEN(6,FILE=PTSAP.IN,STATUS=OLD) OPEN(8,FILE=PTSAP.OUT,STATUS=NEW) 1 READ(6,10)TL 10 FORMAT(20A4) READ(6,11)NJ,N,NNE,NMT,NPJ 11 FORMAT(5I5) IF(NJ.EQ.0)STOP WRITE(8,10)TL WRITE(8,12)NJ,N,NNE,NMT,NPJ 12 FORMAT(10X,NUMBER OF NODES =,I5, /10X,NUMBER OF DEGREE OF FREEDOM =,I5, /10X,NUMBER OF ELEMENTS =,I5, /10X,NUMBER OF MATERIALS =,I5, /10X,NUMBER OF NODAL LOADS =,I5) READ(6,96) ((JN(J,I),J=1,2),X(I),Y(I),I=1,NJ) WRITE(8,97)(I,(JN(J,I),J=1,2),X(I),Y(I),I=1,NJ) 96 FORMAT(2(2I5,2F10.0)) 97 FORMAT(3X,NO.N,2X,NO.DISP.(X,Y.),2X, X-COORDINATE Y-COORDINATE/(3I6,9X,2G14.6)) READ(6,98)((JE(J,I),J=1,2),JEA(I),I=1,NNE) 98 FORMAT(15I5) WRITE(8,103)(I,(JE(J,I),J=1,2),JEA(I),I=1,NNE) 103 FORMAT(8X,NO.E,5X,1(NODE),3X,2(NODE),3X ,NO.MAT/4(I10)) READ(6,105)((EA(I,J),I=1,2),J=1,NMT) 105 FORMAT(6f10.0) WRITE(8,106)(J,(EA(I,J),I=1,2),J=1,NMT) 106 FORMAT(3X,NO.MAT,6X,ELASTIK MODULUS,8X, AREA/(I6,9X,2G16.6)) READ(6,107) (JPJ(I),PJ(I),I=1,NPJ) 107 FORMAT(5(I5,F10.0)) WRITE(8,108)(JPJ(I),PJ(I),I=1,NPJ) 108 FORMAT(20X,NODAL LOAD(NO.DISP.,VALUE)/4(I6,G12.5)) ND=0 DO 113 IE=1,NNE CALL CALM(M,IE,JE,JN) MX=0 MI=N DO 112 I=1,4 L=M(I) IF(L.EQ.0)GOTO 112 IF(L.GT.MX) MX=L IF(L.LT.MI) MI=L 112 CONTINUE 113 IF(ND.LT.MX-MI)ND=MX-MI ND=ND+1 WRITE(8,114)ND 114 FORMAT(/10X,5(* *), RESULTS OF CALCULATION, 5(* *)/15X,SEMI-BAND WIDTH=,I3) IF(ND.GT.100)STOP DO 120 J=1,ND DO 120 I=1,N 120 AK(I,J)=0. DO 125 I