一平面桁架静力分析源程序PTSAP.FOR.doc

一、平面桁架静力分析源程序——PTSAP．FOR C C C PLANE TRUSS STRUCTURAL ANALYSIS PROGRAM C C MAIN PROGRAM DIMENSION JE(2,100),JS(3,50),X(100),Y(100),E(100),A(100), PJ(3,50),P(200),F(100) REAL KB(200,50) OPEN(5,FILE=PTR.DAT,STATUS=OLD) OPEN(6,FILE=PTW.DAT,STATUS=NEW) C INPUT CONTROL PARAMETERS OF STRUCTURE READ(5,*) NJ,NE,NS,NPJ CALL DAT(NJ,NE,NS,NPJ,NW,NEQ,X,Y,JE,JS,PJ,A,E) CALL SSM(NEQ,NE,NJ,NW,JE,X,Y,E,A,KB) CALL FPJ(NEQ,NPJ,PJ,P) CALL CSE(NEQ,NS,NW,JS,KB,P) CALL BGS(NEQ,NJ,NW,KB,P) CALL TRS(NEQ,NE,NJ,NS,NPJ,JE,JS,X,Y,E,A,PJ,P,F) CLOSE(5) CLOSE(6) END C C SUBPROGRAM- DATA INPUT AND PRINT SUBROUTINE DAT(NJ,NE,NS,NPJ,NW,NEQ,X,Y,JE,JS,PJ,A,E) DIMENSION JE(2,NE),JS(3,NS),X(NJ),Y(NJ),E(NE),A(NE),PJ(3,NPJ) READ(5,*)(X(I),Y(I),I=1,NJ) READ(5,*)((JE(I,J),I=1,2),J=1,NE) READ(5,*) ((JS(I,J),I=1,3),J=1,NS) READ(5,*) ((PJ(I,J),I=1,3),J=1,NPJ) READ(5,*) (A(I),E(I),I=1,NE) NW=0 DO 70 I=1,NE IW=IABS(JE(1,I)-JE(2,I)) IF(IW.GT.NW) NW=IW 70 CONTINUE NW=2*(NW+1) NEQ=2*NJ WRITE(6,75) 75 FORMAT(1X,PLANE TRUSS STRUCTURE ANALYSIS/) C OUTPUT CONTROL PARAMETERS OF STRUCTURE WRITE(6,130) WRITE(6,140) NE,NJ,NS,NPJ,NEQ,NW C OUTPUT DATA OF ELEMENTS AND JOINTS WRITE(6,150) WRITE(6,78) 78 FORMAT(/1X,NODAL COODINATES) WRITE(6,80) (I,X(I),Y(I),I=1,NJ) WRITE(6,90) (J,(JE(I,J),I=1,2),A(J),E(J),J=1,NE) WRITE(6,100) ((JS(I,J),I=1,3),J=1,NS) WRITE(6,110) ((PJ(I,J),I=1,3),J=1,NPJ) 80 FORMAT(7X,NODE,7X,X,14X,Y,/(1X,I10,2E15.6)) 90 FORMAT(/4X,ELEMENT,4X,NODE-I,4X,NODE-J,7X,A,14X,E /(1X,3I10,2E15.6)) 100 FORMAT(/1X,CONSTRAINTED NODES/7X,NODE, 8X,XX,8X,YY/(1X,3I10)) 110 FORMAT(/1X,JOINT LOADS/7X,NODE,8X,PX,11X,PY /(1X,F10.0,2E15.6)) 130 FORMAT(1X,STRUCTURAL PARAMETERS//1X,TOTAL,1X,NO. OF) 140 FORMAT(1X,ELEMENTS=,I5/1X,JOINTS=,