平面杆结构代码.docxVIP

图1.计算简图图1所示为本算例的计算简图，已标注了单元号，结点号，以及荷载的数值，此次编程采用的是fortran语言，源代码见下面。FORTRAN源程序代码说明如下：!LOC(I,1),LOC(I,2),AREA(I),FINT(I)!第I个单元的节点号I,J，截面面积A，截面惯性矩I。!CX(I),CY(J)!第J个节点的X,Y坐标!IFIX(K)!第K个给定约束的自由度号!NN：节点个数 NE：单元个数 ND：总自由度个数 NFIX:被约束的自由度个数!NP,NVD!集中荷载作用点个数和横向均布荷载力荷载组数!如果NP≠0，顺序输入NP组数! I F(3*I-2) F(3*I-1) F(3*I)! 节点号 X方向集中力 Y方向集中力弯矩!如果NVD≠0，顺序输入NVD组数! IE Q! 单元号均布荷载集度DIMENSION LOC(9,2),IFIX(6),AREA(9),FINT(9),CX(8),CY(8),F(24),FR(24),FE(9,6) COMMON NN,NE,ND,NFIX,EOPEN(1,file=F:\FORTRAN\1.txt,status=old)OPEN(2,file=F:\FORTRAN\2.txt) !请修改文件路径READ(1,*)NN,NE,ND,NFIX,E WRITE(2,105)NN,NE,ND,NFIX,E 105 FORMAT(4X,NN NE ND NFIX E/4I5,E11.4) READ(1,*)(LOC(I,1),LOC(I,2),AREA(I),FINT(I),I=1,NE) WRITE(2,106) 106 FORMAT(/1X,ELEMENT NODE1 NODE2 AREA FINT) WRITE(2,107)(I,LOC(I,1),LOC(I,2),AREA(I),FINT(I),I=1,NE) 107 FORMAT(1X,3I6,2E12.4) READ(1,*)(CX(I),CY(I),I=1,NN) WRITE(2,108) 108 FORMAT(/4X,NODE X-COORD Y-COORD) WRITE(2,109)(I,CX(I),CY(I),I=1,NN) 109 FORMAT(4X,I3,3X,2E12.4) READ(1,*)(IFIX(K),K=1,NFIX) WRITE(2,110) (IFIX(K),K=1,NFIX) 110 FORMAT(/3X,IFIX=,6I4) CALL CST(LOC,IFIX,AREA,FINT,CX,CY,F,FR,FE) WRITE(2,111) 111 FORMAT(/4X,NODE,5X,X-DISP,8X,Y-DISP,9X,THITA) WRITE(2,112) (I,F(3*I-2),F(3*I-1),F(3*I),I=1,NN) 112 FORMAT(2X,I5,3E15.6) !WRITE(2,113) !113 FORMAT(/4X,REACTION NODAL FORCES FROM THE EQUATION,/4X,NODE,5X,X-LOAD,8X,Y-LOAD,9X,MOMENT) !WRITE(2,114)(I,FR(3*I-2),FR(3*I-1),FR(3*I),I=1,NN) !114 FORMAT(2X,I5,3E15.6) WRITE(2,115) 115 FORMAT(/4X,ELEMENT,4X,AI-F,7X,SHEAR-Q,7X,MOMENT-M) DO 10 I=1,NE 10 WRITE(2,116)I, FE(I,1),FE(I,2),FE(I,3),FE(I,4),FE(I,5),FE(I,6) 116 FORMAT(2X,I4,3E15.6/6X,3E15.6) STOP END 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),F(ND), FR(ND),F1(ND), FE(NE,6),P1(NE,6),P2(NE,6), GK(ND,ND),GK1(ND,ND),AL(NE),TT(NE,6,6),BKL(NE,6,6), T(6,6),CSS(NE),SNN(NE),EK(6,6),EKL