實(shí)驗(yàn)三平面剛架程序設(shè)計(jì)

1、實(shí)驗(yàn)三 平面剛架程序設(shè)計(jì)開 始標(biāo)題及數(shù)組說明（讀入題目序號NO）NO是否為零?形成整體剛度矩陣計(jì)算并打印各桿軸力解方程并打印桿端位移形成結(jié)點(diǎn)載荷結(jié) 束子程序READ子程序MKE子程序MAKE子程序MR子程序MF是否讀入數(shù)據(jù)并打印子程序MULV6子程序CALM子程序MK子程序MULV子程序TARN子程序SOLV子程序MADE子程序PE一 平面剛架內(nèi)力和位移計(jì)算總框圖二 平面桁架靜力分析源程序（PFSAP.FOR）C ANALYSIS PROGRAM FOR PLANE FRAME REAL K(200,200),KE(6,6),AKE(6,6),X(100),Y(100),AL(100), &a

2、mp; 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,*)N

3、O 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,N5 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

4、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,NPJ30 P(JPJ(I)=P(JPJ(I)+PJ(I) CALL SLOV(K,P,D,N) WRITE(8,'(/2(26(1H*),A)') 'RESULTS OF CALC

5、ULATION' 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(/'N

6、O.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)

7、 THEN CALL PE(FE,IP,JPF,PF,AL) DO 90 I=1,690 FF(I)=FF(I)-FE(I) ENDIF100 CONTINUE WRITE(8,110) IE,(FF(I),I=1,6)110 FORMAT(I5,2X,6G12.5)130 CONTINUE GOTO 1END SUBROUTINE READ(NJ,N,NE,NM,NPJ,NPF,JN,X,Y,JE,JEAI,EAI, & JPJ,PJ,JPF,PF) REAL X(100),Y(100),EAI(3,100),PJ(100),PF(2,100) INTEGER JE(2,100),J

8、N(3,100),JPJ(100),JPF(2,100),JEAI(100), & TITLE(20) READ(6,'(20A4)') (TITLE(I),I=1,20) WRITE(8,'(/7X,20A4)')TITLE READ (6,*) NJ,N,NE,NM,NPJ,NPF WRITE(8,'(/3(5X,A4,1H:,I2)')'NJ=',NJ, & 'N=',N,'NE=',NE,'NM=',NM,'NPJ=',NPJ,'NPF

9、=',NPF WRITE (8,5)5 FORMAT(/4X,'NO. (1) (2) (3)',10X,'X ',8X,'Y') READ (6,10)(JN(J,I),J=1,3),X(I),Y(I),I=1,NJ)10 FORMAT(2(3I5,2G16.4) DO 20 I=1,NJ20 WRITE (8,'(2X,1H(,I2,1H),3I6,4X,2F10.3)') I,JN(1,I),JN(2,I), & JN(3,I),X(I),Y(I) WRITE (8,30)30 FORMAT(/10X,

10、9;ELEMENT NO. NODE-1 NODE-2 MATERIALS') READ (6,40)(JE(1,I),JE(2,I),JEAI(I),I=1,NE)40 FORMAT(5(3I5) DO 50 I=1,NE50 WRITE (8,'(14X,I2,3(7X,I3)') I,JE(1,I),JE(2,I),JEAI(I) READ(6,*)(EAI(I,J),I=1,3),J=1,NM) WRITE(8,60)(J,(EAI(I,J),I=1,3),J=1,NM)60 FORMAT(/3X,'NO.MAT',6X,'ELASTIK

11、 MODULUS',8X, & 'AREA',5X,'MOMENT OF INERTIA'/(I6,9X,3G16.6) IF(NPJ.EQ.0) GOTO 90 WRITE(8,'(/20X,16H NODEL LOADS )') WRITE(8,'(16XA)')' NO.DISP. VALUE' READ (6,70) (JPJ(I),PJ(I),I=1,NPJ)70 FORMAT (5(I5,G16.4) DO 80 I=1,NPJ80 WRITE(8,'(14X,I7,F16.3)

12、') JPJ(I),PJ(I)90 CONTINUE IF(NPF.EQ.0) GOTO 130 WRITE(8,'(/20X,16HNON-NODEL LOADS )') WRITE(8,'(11X,A,8X,A,9X,A)')'NO.E NO.LOAD.MODEL','A','C' READ (6,100) (JPF(1,I),JPF(2,I),PF(1,I),PF(2,I),I=1,NPF)100 FORMAT (2(2I5,2G16.4) DO 110 I=1,NPF110 WRITE(8,120)

13、 (JPF(J,I),J=1,2),PF(1,I),PF(2,I)120 FORMAT(6X,2I8,10X,2F10.3)130 CONTINUE RETURN END SUBROUTINE MKE(KE,IE,JE,JEAI,EAI,X,Y,AL) REAL KE(6,6),X(100),Y(100),EAI(3,100),AL(100),L INTEGER JE(2,100),JEAI(100) II=JE(1,IE) JJ=JE(2,IE) MT=JEAI(IE) L=SQRT(X(JJ)-X(II)*2+(Y(JJ)-Y(II)*2) AL(IE)=L A1=EAI(1,MT)*EA

14、I(2,MT)/L A2=EAI(1,MT)*EAI(3,MT)/L*3 A3=EAI(1,MT)*EAI(3,MT)/L*2 A4=EAI(1,MT)*EAI(3,MT)/L KE(1,1)=A1 KE(1,4)=-A1 KE(2,2)=12*A2 KE(2,3)=6*A3 KE(2,5)=-12*A2 KE(2,6)=6*A3 KE(3,3)=4*A4 KE(3,5)=-6*A3 KE(3,6)=2*A4 KE(4,4)=A1 KE(5,5)=12*A2 KE(5,6)=-6*A3 KE(6,6)=4*A4 DO 10 I=1 ,6 DO 10 K=I ,610 KE(K,I)=KE(I,

15、K) RETURN END SUBROUTINE MR(R,IE,JE,X,Y) REAL R(6,6),X(100),Y(100),L,CX,CY INTEGER JE(2,100) I=JE(1,IE) J=JE(2,IE) L=SQRT(X(J)-X(I)*2+(Y(J)-Y(I)*2) CX=(X(J)-X(I)/L CY=(Y(J)-Y(I)/L DO 10 J=1,6 DO 10 I=1,610 R(I,J)=0. DO 20 I=1,4,3 R(I,I)=CX R(I,I+1)=CY R(I+1,I)=-CY R(I+1,I+1)=CX20 R(I+2,I+2)=1. RETUR

16、N END SUBROUTINE MAKE(KE,R,AKE) REAL KE(6,6),R(6,6),RT(6,6),TMP(6,6),AKE(6,6) CALL TRAN(R,RT) CALL MULV(RT,KE,TMP) CALL MULV(TMP,R,AKE) RETURN END SUBROUTINE CALM(M,IE,JN,JE) INTEGER M(6),JN(3,100),JE(2,100),IE DO 10 I=1,3 M(I)=JN(I,JE(1,IE)10 M(I+3)=JN(I,JE(2,IE) RETURN END SUBROUTINE MK(K,AKE,M) R

17、EAL K(200,200),AKE(6,6) INTEGER M(6) DO 10 I=1,6 DO 10 J=1,6 IF(M(I).NE.0.AND.M(J).NE.0) & K(M(I),M(J)=K(M(I),M(J)+AKE(I,J)10CONTINUE RETURN END SUBROUTINE PE(FE,IP,JPF,PF,AL) REAL FE(6),PF(2,100),AL(100),L INTEGER JPF(2,100) A=PF(1,IP) C=PF(2,IP) L=AL(JPF(1,IP) IND=JPF(2,IP) DO 5 I=1,65 FE(I)=0

18、. GOTO(10,20,30,40,50,60),IND10 FE(2)=(7*A/20+3*C/20)*L FE(3)=(A/20+C/30)*L*2 FE(5)=(3*A/20+7*C/20)*L FE(6)=-(A/30+C/20)*L*2 RETURN20 FE(5)=A*C*3*(2*L-C)/2/L*3 FE(2)=A*C-FE(5) FE(3)=A*C*2*(6*L*L-8*C*L+3*C*C)/12/L/L FE(6)=-A*C*3*(4*L-3*C)/12/L/L RETURN30 FE(2)=A*(L-C)*2*(L+2*C)/L*3 FE(3)=A*C*(L-C)*2/

19、L*2 FE(5)=A-FE(2) FE(6)=-A*C*2*(L-C)/L*2 RETURN40 FE(2)=-6*A*C*(L-C)/L*3 FE(3)=A*(L-C)*(L-3*C)/L*2 FE(5)=-FE(2) FE(6)=A*C*(3*C-2*L)/L*2 RETURN50 FE(1)=A*(1-C/L) FE(4)=A*C/L RETURN60 FE(1)=C*L/2. FE(4)=FE(1) RETURN END SUBROUTINE MULV6(A,B,C) REAL C(6),A(6,6),B(6) DO 10 I=1,6 C(I)=0. DO 10 J=1,610 C(

20、I)=C(I)+A(I,J)*B(J) RETURN END SUBROUTINE MF(P,AFE,M) REAL P(100),AFE(6) INTEGER M(6) DO 10 I=1,6 IF(M(I).NE.0)P(M(I)=AFE(I)+P(M(I)10 CONTINUE RETURN END SUBROUTINE SLOV(AK,P,D,N) REAL AK(200,200),P(100),D(100) DO 5 I=1,1005 D(I)=P(I) DO 10 K=1,N-1 DO 10 I=K+1,N C=-AK(K,I)/AK(K,K) DO 20 J=I,N20 AK(I

21、,J)=AK(I,J)+C*AK(K,J)10 D(I)=D(I)+C*D(K) D(N)=D(N)/AK(N,N) DO 40 I=N-1,1,-1 DO 30 J=I+1,N30 D(I)=D(I)-AK(I,J)*D(J)40 D(I)=D(I)/AK(I,I)RETURN END SUBROUTINE MADE(IE,JN,JE,D,ADE) REAL ADE(6),D(100) INTEGER IE,JN(3,100),JE(2,100) DO 3 I=1,63 ADE(I)=0 DO 10 I=1,3 IF (JN(I,JE(1,IE).NE.0) ADE(I)=D(JN(I,JE

22、(1,IE) IF (JN(I,JE(2,IE).NE.0) ADE(I+3)=D(JN(I,JE(2,IE)10 CONTINUE RETURN END SUBROUTINE TRAN(R,RT) REAL R(6,6),RT(6,6) DO 10 I=1,6 DO 10 J=1,610 RT(I,J)=R(J,I) RETURN END SUBROUTINE MULV(A,B,C) REAL A(6,6),B(6,6),C(6,6) DO 10 I=1 , 6 DO 10 J=1 , 6 C(I,J)=0. DO 10 K=1 , 610 C(I,J)=C(I,J)+A(I,K)*B(K,

23、J) RETURN END三 課后習(xí)題（1） 準(zhǔn)備原始數(shù)據(jù)確定結(jié)點(diǎn)、劃分單元，建立整體坐標(biāo)系與局部坐標(biāo)系如下圖所示。自由結(jié)點(diǎn)位移編碼如下圖所示。q1(0,0,0)2(1,2,3)4(0,0,0)/xyLLL3(0,0,0)設(shè)本題標(biāo)題為EXAMPLE(5.1)，整理并計(jì)算本題所需的數(shù)據(jù)填入下表：題號1標(biāo)題EXAMPLE(5.1)基本數(shù)據(jù)結(jié)點(diǎn)總數(shù)自由度單元總數(shù)單元類型數(shù)結(jié)點(diǎn)載荷數(shù)非結(jié)點(diǎn)載荷數(shù)433101結(jié)點(diǎn)數(shù)據(jù)起止號結(jié)點(diǎn)位移號結(jié)點(diǎn)坐標(biāo)結(jié)點(diǎn)位移號結(jié)點(diǎn)坐標(biāo)uvxyuvxy120000.0.1235.0.3400010.0.0005.5.單元數(shù)據(jù)起止單元1端節(jié)點(diǎn)號2端節(jié)點(diǎn)號單元類型號1端節(jié)點(diǎn)號2端節(jié)點(diǎn)號

24、單元類型號1端節(jié)點(diǎn)號2端節(jié)點(diǎn)號單元類型號1端結(jié)點(diǎn)號2端節(jié)點(diǎn)號單元類型號1端節(jié)點(diǎn)號2端節(jié)點(diǎn)號單元類型號121231241單元類型數(shù)據(jù)類型號彈性模量橫截面積慣性矩13.E80.50.04167結(jié)點(diǎn)載荷數(shù)據(jù)起止序號位移號數(shù)值位移號數(shù)值位移號數(shù)值位移號數(shù)值位移號數(shù)值15非結(jié)點(diǎn)載荷起止號單元號類型參數(shù)a參數(shù)c單元號類型參數(shù)a參數(shù)c12114.84.834（2） 建立數(shù)據(jù)文件PFSAP.IN并輸入以下數(shù)據(jù)：1 EXAMPLE(5.1)4,3,3,1,0,10,0,0,0.,0.,1,2,3,5.,0.0,0,0,10.,0.,0,0,0,5.,5.1,2,1,2,3,1,2,4,13.E8,0.5,0.041671,1,4.8,4.80（3