有限元钢架源程序

C-FRAME.FOR-(FOR FRAME PROBLEMS)- DIMENSION LOC(3,2),IFIX(6),AREA(3),FINT(3),CX(4),CY(4), 1 F(12),FR(12),FE(3,6) COMMON NN,NE,ND,NFIX,EC- OPEN (5,FILE=F1.TXT,STATUS=OLD) OPEN (6,FILE=OUTF1.) READ(5,*) NN,NE,ND,NFIX,E WRITE(6,105) NN,NE,ND,NFIX,E105 FORMAT(4X,NN NE ND NFIX E/4I5,E11.4) READ(5,*) (LOC(I,1),LOC(I,2),AREA(I),FINT(I),I=1,NE) WRITE(6,108)108 FORMAT(/1X,ELEMENT NODE1 NODE2 AREA FINT) WRITE(6,110) (I,LOC(I,1),LOC(I,2),AREA(I),FINT(I),I=1,NE)110 FORMAT(1X,3I6,2E12.4) READ(5,*) (CX(J),CY(J),J=1,NN) WRITE(6,112)112 FORMAT(/4X,NODE X-COORD Y-COORD) WRITE(6,115) (J,CX(J),CY(J),J=1,NN)115 FORMAT(4X,I3,3X,2E12.4) READ(5,*) (IFIX(K),K=1,NFIX) WRITE(6,118) (IFIX(K),K=1,NFIX)118 FORMAT(/3X,IFIX=,6I4)C- CALL CST(LOC,IFIX,AREA,FINT,CX,CY,F,FR,FE)C- WRITE(6,130)130 FORMAT(/4X,NODE,5X,X-DISP,8X,Y-DISP,9X,THITA) WRITE(6,135) (I,F(3*I-2),F(3*I-1),F(3*I),I=1,NN)135 FORMAT(2X,I5,3E15.6) WRITE(6,140)140 FORMAT(/4X,reaction nodal forces from the equations, 1 /4X,NODE,5X,X-LOAD,8X,Y-LOAD,9X,MOMENT) WRITE(6,145) (I,FR(3*I-2),FR(3*I-1),FR(3*I),I=1,NN)145 FORMAT(2X,I5,3E15.6) WRITE(6,150)150 FORMAT(/4X,ELEMENT,4X,AXI-F,7X,SHEAR-Q,7X,MOMENT-M) DO 10 I=1,NE10 WRITE(6,155) I,FE(I,1),FE(I,2),FE(I,3),FE(I,4),FE(I,5),FE(I,6)155 FORMAT(2X,I4,3E15.6/6X,3E15.6) STOP ENDC- SUBROUTINE CST(LOC,IFIX,AREA,FINT,CX,CY,F,FR,FE) DIMENSION LOC(NE,2),IFIX(NFIX),AREA(NE),FINT(NE),CX(NN),CY(NN), 1 F(ND),FR(ND),FE(NE,6), 2 P1(3,6),P2(3,6),GK(12,12),GK1(12,12),AL(3),TT(3,6,6),BKL(3,6,6), 3 T(6,6),CSS(3),SNN(3),EK(6,6),EKL(6,6),EKK(3,6,6),XX(6),BA(6,6) COMMON NN,NE,ND,NFIX,EC- DO 5 I=1,NE I1=LOC(I,1) I2=LOC(I,2) X12=CX(I2)-CX(I1) Y12=CY(I2)-CY(I1) AL(I)=SQRT(X12*2+Y12*2) CSS(I)=X12/AL(I) SNN(I)=Y12/AL(I)5 CONTINUEC FOR CONCENTRAL FORCE, SET TO A NODE READ (5,*) NP,NVD IF (NP.EQ.0) GO TO 15 DO 10 II=1,NP READ (5,*) I,F(3*I-2),F(3*I-1),F(3*I)10 CONTINUE15 CONTINUEC FOR VERTICAL DISTRIBUTED LOADING OF AN ELEMENT IF (NVD.EQ.0) GO TO 25 DO 20 II=1,NVD READ(5,*) IE,Q I1=LOC(IE,1) I2=LOC(IE,2) P1(IE,3*I1-1)= Q*AL(IE)/2. P1(IE,3*I1) = Q*AL(IE)*2/12. P1(IE,3*I2-1)= Q*AL(IE)/2. P1(IE,3*I2) =-Q*AL(IE)*2/12.20 CONTINUE25 CONTINUEC- DO 30 I=1,ND DO 30 J=1,ND30 GK(I,J)=0.0C- DO 100 I=1,NE DO 35 J=1,6 DO 35 K=1,6 EKL(J,K)=0.0 EK(J,K)=0.035 T(J,K)=0.0C- EAL=E*AREA(I)/AL(I) EIL1=E*FINT(I)/AL(I) EIL2=E*FINT(I)/AL(I)*2 EIL3=E*FINT(I)/AL(I)*3 EKL(1,1)=EAL EKL(2,2)=12.*EIL3 EKL(3,3)= 4.*EIL1 EKL(4,4)=EAL EKL(5,5)=12.*EIL3 EKL(6,6)= 4.*EIL1 EKL(3,2)= 6.*EIL2 EKL(4,1)=-EAL EKL(5,2)=-12.*EIL3 EKL(5,3)= -6.*EIL2 EKL(6,2)= 6.*EIL2 EKL(6,3)= 2.*EIL1 EKL(6,5)= -6.*EIL2 DO 50 II=1,5 DO 50 JJ=II+1,650 EKL(II,JJ)=EKL(JJ,II) DO 55 K=1,6 DO 55 L=1,6 EKK(I,K,L)=EKL(K,L)C WRITE(6,140) I,K,L,EKL(K,L),EKK(I,K,L)140 FORMAT(1X,IKL,EK,3I4,2E12.5)55 CONTINUEC- T(1,1)= CSS(I) T(1,2)=-SNN(I) T(2,1)= SNN(I) T(2,2)= CSS(I) T(3,3)= 1. T(4,4)= CSS(I) T(4,5)=-SNN(I) T(5,4)= SNN(I) T(5,5)= CSS(I) T(6,6)= 1.C- DO 60 J=1,6 DO 60 K=1,6 TT(I,J,K)=T(J,K)60 P2(I,J)=P2(I,J)+T(J,K)*P1(I,K)C- DO 65 II=1,6 DO 65 J=1,6 BA(II,J)=0.0 DO 65 K=1,6 BA(II,J)=BA(II,J)+TT(I,II,K)*EKL(K,J)65 CONTINUE DO 70 II=1,6 DO 70 J=1,6 EK(II,J)=0.0 DO 70 K=1,6 EK(II,J)=EK(II,J)+BA(II,K)*TT(I,J,K)70 CONTINUEC WRITE(6,120) (I,II,J,EK(II,J),TT(I,II,J),II=1,6),J=1,6)120 FORMAT(1X,I,II,J,EK,TT=,3I4,2E12.5)C- DO 85 INODE=1,2 NODEI=LOC(I,INODE) DO 85 IDOFN=1,3 NROWS=(NODEI-1)*3+IDOFN NROWE=(INODE-1)*3+IDOFN F(NROWS)=F(NROWS)+P2(I,NROWE) DO 85 JNODE=1,2 NODEJ=LOC(I,JNODE) DO 85 JDOFN=1,3 NCOLS=(NODEJ-1)*3+JDOFN NCOLE=(JNODE-1)*3+JDOFN85 GK(NROWS,NCOLS)=GK(NROWS,NCOLS)+EK(NROWE,NCOLE)100 CONTINUE DO 88 I=1,ND DO 88 J=1,ND88 GK1(I,J)=GK(I,J) WRITE(6,160)160 FORMAT(/4X,nodal forces from applied loads, 1 /4X,NODE,5X,X-LOAD,8X,Y-LOAD,9X,MOMENT) WRITE(6,165) (I,F(3*I-2),F(3*I-1),F(3*I),I=1,NN)165 FORMAT(2X,I5,3E15.6)C WRITE(6,170) (I,J,GK(I,J),J=1,ND),I=1,ND)170 FORMAT(1X,I,J,GK1,2I4,E12.5,2X,2I4,E12.5)C- DO 90 I=1,NFIX IX=IFIX(I)90 GK(IX,IX)=GK(IX,IX)*1.0E15 CALL GAUSS(GK,F,ND)C- DO 92 I=1,ND FR(I)=0.0 DO 94 J=1,ND94 FR(I)=FR(I)+GK1(I,J)*F(J)92 FR(I)=FR(I)-F(I)C- DO 95 I=1,NE I1=LOC(I,1) I2=LOC(I,2) XX(1)=F(3*I1-2) XX(2)=F(3*I1-1) XX(3)=F(3*I1) XX(4)=F(3*I2-2) XX(5)=F(3*I2-1) XX(6)=F(3*I2) DO 96 J=1,6 DO 96 K=1,6 BKL(I,J,K)=0.0 DO 96 L=1,6 BKL(I,J,K)=BKL(I,J,K)+EKK(I,J,L)*TT(I,K,L)96 CONTINUE DO 97 J=1,6 FE(I,J)=0.0 DO 97 K=1,6 FE(I,J)=FE(I,J)+BKL(I,J,K)*XX(K)97 CONTINUE DO 98 J=1,698 FE(I,J)=FE(I,J)-P1(I,J)95 CONTINUE RETURN ENDC-GAUSSS METHOD- SUBROUTINE GAUSS(A,B,N) DIMENSION A(N,N),B(N) DO 1 I=1,N I1=