实验三 平面刚架程的序设计.doc

实验三 平面刚架程序设计 平面刚架内力和位移计算总框图 平面桁架静力分析源程序（PFSAP.FOR） C ANALYSIS PROGRAM FOR PLANE FRAME REAL K(200,200),KE(6,6),AKE(6,6),X(100),Y(100),AL(100), EAI(3,100),PJ(100),PF(2,100),R(6,6),P(100),FF(6), FE(6),D(100),ADE(6),DE(6),RT(6,6),AFE(6),F(3) INTEGER JE(2,100),JN(3,100),JPJ(100),JPF(2,100),M(6), JEAI(100),NO OPEN (6,FILE=PFSAP.IN,STATUS=OLD) OPEN (8,FILE=PFSAP.OUT,STATUS=NEW) 1 READ (6,*)NO IF(NO.EQ.0)STOP WRITE (8,(/A5,I3,A1))(NO.=,NO,) CALL READ(NJ,N,NE,NM,NPJ,NPF,JN,X,Y,JE,JEAI,EAI,JPJ,PJ,JPF,PF) DO 5 I=1,N P(I)=0. DO 5 J=1,N 5 K(I,J)=0 DO 10 IE=1,NE CALL MKE(KE,IE,JE,JEAI,EAI,X,Y,AL) CALL MR(R,IE,JE,X,Y) CALL MAKE(KE,R,AKE) CALL CALM(M,IE,JN,JE) CALL MK(K,AKE,M) 10 CONTINUE DO 20 IP=1,NPF CALL MR(R,JPF(1,IP),JE,X,Y) CALL TRAN(R,RT) CALL PE(FE,IP,JPF,PF,AL) CALL MULV6 (RT,FE,AFE) CALL CALM(M,JPF(1,IP),JN,JE) CALL MF(P,AFE,M) 20 CONTINUE DO 30 I=1,NPJ 30 P(JPJ(I))=P(JPJ(I))+PJ(I) CALL SLOV(K,P,D,N) WRITE(8,(/2(26(1H*),A))) RESULTS OF CALCULATION WRITE(8,40) 40 FORMAT(/5X,NO.N,4X,X-DISPLACEMENT,2X, Y-DISPLACEMENT,3X,ANG.ROT.(RAD)) DO 60 KK=1,NJ DO 50 II=1,3 F(II)=0. I1=JN(II,KK) 50 IF(I1.GT.0)F(II)=D(I1) 60 WRITE(8,70)KK,F(1),F(2),F(3) 70 FORMAT(I8,2X,3G16.5) WRITE(8,80) 80 FORMAT(/NO.E,5X,N(1),8X,Q(1),8X,M(1), 8X,N(2),8X,Q(2),8X,M(2)) DO 130 IE=1,NE CALL MADE(IE,JN,JE,D,ADE) CALL MKE(KE,IE,JE,JEAI,EAI,X,Y,AL) CALL MR(R,IE,JE,X,Y) CALL MULV6(R,ADE,DE) CALL MULV6(KE,DE,FF) DO 100 IP=1,NPF IF (JPF(1,IP).EQ.IE) THEN CALL PE(FE,IP,JPF,PF,AL)