结构力学电算实习程序.docVIP

PROGRAM LXL DIMENSION GC(20),GX(20),PJ(20,2),PF(10,4),DK(2,2),P(45) DIMENSION WY(2),ZK1(45),ZK2(45) DIMENSION F(2),F0(2),F1(20),F2(20) DIMENSION FF(2),FF0(2),FF1(20),FF2(20) character*12 dat0 write(*,*)输入:初始数据文件名 read(*,10)dat0 10 format(A12) open(8,file=dat0,status=old) open(9,file=连续梁内力计算结果.txt,status=new) write(9,11) 11 format(8x,*************连续梁内力计算*****************) READ(8,*)NE,NZ,NP,NF,E0 WRITE(9,12)NE,NZ,NP,NF,E0 12 FORMAT(1X,单元数=,6x,I3,5x,支承类型=,I3,5x,节点荷载个数=I3,/,1x,非节点荷载个数=,5x,I3,弹性模量=,F12.4) NJ=NE+1 CALL SRSJ(NE,NP,NF,GC,GX,PJ,PF) PRINT*,计算连续梁内力请输入1 PRINT*,计算指定截面弯矩和剪力请输入2 PRINT*,计算所有截面弯矩和剪力请输入3 PRINT*,计算指定截面弯矩和剪力影响线请输入4 PRINT*,计算指定支座反力影响线请输入5 PRINT*,计算均布活荷载作用下的包络图请输入6 READ*,I !********************************** ! 计算连续梁内力 !********************************** IF(I==1) THEN CALL XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,FF0) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) CALL ZCCL(NZ,NJ,ZK1,ZK2,P) CALL FCQJ(NJ,ZK1,ZK2,P) CALL GDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF0,FF1,FF2) CALL SC(NE,F1,F2,FF1,FF2,NJ,P) !********************************* ! 计算指定截面弯矩和剪力 !********************************* ELSE IF(I==2)THEN CALL XCP(NJ,NP,NF,NE,P,PJ,PF,GC,F0,FF0) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) CALL ZCCL(NZ,NJ,ZK1,ZK2,P) CALL FCQJ(NJ,ZK1,ZK2,P) CALL GDL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF0,FF1,FF2) PRINT*,请输入要计算内力截面位置(离连续梁最左端的距离) READ*,WZ MT=0.0 DO I=1,NE ML=MT MT=MT+GC(I) IF((WZ=ML).AND.(WZ=MT))THEN WZ=WZ-ML NTWZ=I ENDIF ENDDO CALL WJJL(NE,F1,FF1,NF,PF,SHEAR,TORTION,WZ,NTWZ) PRINT*, print*,弯矩=,TORTION,剪力=,SHEAR PRINT*, call SC(NE,F1,F2,FF1,FF2,NJ,P) !********************************** ! 计算所有截面弯矩和剪力 　　 　　 !********************************** ELSE IF(I==3) THEN CALL NL(NE,NJ,NF,E0,P,WY,DK,PF,GC,GX,F,F0,F1,F2,FF,FF0,FF1,FF2,SHEAR,TORTION,ZK1,ZK2,NZ,PJ,NP) call SC(NE,F1,F2,FF1,FF2,NJ,P) !***********