《连续梁计算程序》由会员分享，可在线阅读，更多相关《连续梁计算程序（7页珍藏版）》请在金锄头文库上搜索。
1、连续梁计算程序连续梁计算程序! 连续梁静力计算程序PROGRAM LXLDIMENSION GC(20),GX(20),PJ(20,2),PF(10,4),DK(2,2),P(45),F0(2),WY(2),F(2),ZK1(45),ZK2(45)WRITE(*,11)11 FORMAT(1X,单元数,支承类型,节点荷载,非节点荷载,弹性模量)READ(*,*)NE,NZ,NP,NF,E0NJ=NE+1! 输入初始数据 CALL SRSJ(NE,NP,NF,GC,GX,PJ,PF)! 形成总荷载矩阵 PCALL XCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)! 集成整体刚度矩阵
2、ZKl，ZK2CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)! 支承条件处理CALL ZCCL(NZ,NJ,ZK1,ZK2,P)! 方程求解-求位移 PCALL FCQJ(NJ,ZK1,ZK2,P)! 输出位移CALL SCWY(NJ,P)! 计算杆端内力CALL GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX)! 计算结束END ! *! 输入初始数据 SRSJ 子程序! *SUBROUTINE SRSJ(NE,NP,NF,GC,GX,PJ,PF)DIMENSiON GC(NE),GX(NE),PJ(NP,2),PF(NF,4)! GC(
3、NE)：杆长 GX(NE)：惯性矩! 结点荷载:PJ(NP,1)=荷载大小 PJ(NP,2)=对应未知数序号 ! 非结点荷载 : PF(NE,1)=荷载大小 PF(NE,2)=距离 PF(NE,3)=单元号 PF(NE,4)=荷载类型号WRITE(*,11)11 FORMAT(1X,杆长，惯性矩 GC(NE),GX(NE)READ(*,*)(GC(I),GX(I),I=1,NE)IF(NP.GT.0) THENWRITE(*,22)22 FORMAT(1X,荷载大小，对应未知数序号 PJ(I,1)，PJ(I,2)READ(*,*)(PJ(I,1),PJ(I,2),I=1,NP)END IFIF
4、(NF.GT.0) THENWRITE(*,33) 33 FORMAT(1X,非结点荷载值，距离，单元号+荷载类型号)READ(*,*)(PF(I,J),J=1,4),I=1,NF)END IFEND ! *#*! 计算第 NHF 个非结点荷截引起的等效结点荷截列阵 F0! *SUBROUTINE DJH(NHF,NE,NF,PF,F0,GC)DIMENSION PF(NF,4),F0(2),GC(NE)! G荷载值，C作用距离，NT-单元号，ID荷载类型号G=PF(NHF,1)C=PF(NHF,2)NT=INT(PF(NHF,3)+0.1)ID=INT(PF(NHF,4)+0.1)BL=GC
5、(NT)D=BL-CC1=C/BLC2=C1*C1GOTO(10,20), ID! 均布荷载作用下的等效荷载列阵 F0(2)10 F0(1)=G*C*C*(6.-8.*C1+3.*C2)/12.0F0(2)=-G*C*C*C*(4.-3.*C1)/12.0/BLGOTO 200! 集中力作用下的等效荷截列阵 P0(2)20 F0(1)=G*C*D*D/BL/BLF0(2)=-G*D*C2200 RETURNEND! *! 计算第 NE0 个单元的单元刚度矩阵 DK(2，2)! *SUBROUTINE DG(NE0,NE,E0,GC,GX,DK)DIMENSION GC(NE),GX(NE),D
6、K(2,2)DO 15 I=1,2DO 10 J=1,2DK(I,J)=1.010 CONTINUE15 CONTINUE! DL杆长，DI-惯性矩，S-线刚度DL=GC(NE0)DI=GX(NE0)S=E0*DI/DLDK(1,1)=4.0*SDK(1,2)=2.0*SDK(2,1)=2.0*SDK(2,2)=4.0*SEND! *! 集成总体刚度矩阵，ZKl，ZK2 分别存主对角元素和付对角元素! *SUBROUTINE JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)DIMENSION ZK1(NJ),ZK2(NJ),DK(2,2),GC(NE),GX(NE)DO 10 I
7、=1,NJZK1(I)=0.0ZK2(I)=0.010 CONTINUEDO 100 I=1,NECALL DG(I,NE,E0,GC,GX,DK)ZK1(I)=ZK1(I)+DK(1,1)ZK2(I)=ZK2(I)+DK(1,2)ZK1(I+1)=ZK1(I+1)+DK(2,2)100 CONTINUEEND! *! 形成总荷裁矩阵! *SUBROUTINE XCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)DIMENSION P(NJ),PJ(NP,2),PF(NF,4),F0(2),GC(NE)DO 10 I=1,NJP(I)=0.010 CONTINUEIF(NP.GT.0)
8、 THENDO 20 I=1,NPJ=INT(PJ(I,2)+0.1)P(J)=PJ(I,1)20 CONTINUEEND IFIF(NF.GT.0) THENDO 100 NF0=1,NFCALL DJH(NF0,NE,NF,PF,F0,GC)ND=INT(PF(NF0,3)+0.1)P(ND)=P(ND)+F0(1)P(ND+1)=P(ND+1)+F0(2)100 CONTINUEEND IFEND! *! 进行支承条件处理! *SUBROUTINE ZCCL(NZ,NJ,ZK1,ZK2,P)DIMENSION ZK1(NJ),ZK2(NJ),P(NJ)GOTO(10,20,30,20),
9、NZ10 GOTO 10020 ZK1(1)=1.0P(1)=0.0ZK2(1)=0.0IF(NZ.EQ.4) GOTO 30GOTO 10030 ZK1(NJ)=1.0ZK2(NJ-1)=0.0P(NJ)=0.0100 CONTINUEEND! *! 解方程求节点位移 P! *SUBROUTINE FCQJ(NJ,ZK1,ZK2,P)DIMENSION ZK1(NJ),ZK2(NJ),P(NJ)DO 10 I=1,NJ-1ZK1(I+1)=ZK1(I+1)-ZK2(I)*ZK2(I)/ZK1(I)P(I+1)=P(I+1)-ZK2(I)*P(I)/ZK1(I)10 CONTINUEP(NJ)
10、=P(NJ)/ZK1(NJ)DO 20 I=1,NJ-1P(NJ-I)=(P(NJ-I)-ZK2(NJ-I)*P(NJ-I+1)/ZK1(NJ-I)20 CONTINUEEND! *! 输出位移! *SUBROUTINE SCWY(NJ,P)DIMENSION P(NJ)WRITE(*,10)10 FORMAT(1X,：位 移 ：；：)DO 100 I=1,NJWRITE(*,20)I,P(I)20 FORMAT(1X,结点号=,I2,5X,F10.4)100 CONTINUEEND! *! 计算单元杆端内力! *SUBROUTINE GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,
11、PF,GC,GX)DIMENSION P(NJ),WY(2),F(2),F0(2),DK(2,2),PF(NF,4)DIMENSION GC(NE),GX(NE)WRITE(*,2)2 FORMAT(1X,.各单元杆端内力.)DO 200 NE0=1,NECALL DG (NE0,NE,E0,GC,GX,DK)WY(1)=P(NE0)WY(2)=P(NE0+1)DO 35 I=1,2F(I)=0.0DO 30 J=1,2F(I)=F(I)+DK(I,J)*WY(J)30 CONTINUE35 CONTINUEIF(NF.GT.0) THENDO 50 I=1,NFIF(INT(PF(I,3)+0.1).EQ.NE0) THENCALL DJH(I,NE,NF,PF,F0,GC)DO 40 J=1,2F(J)=F(J)-F0(J)40 CONTINUEEND IF50 CONTINUEEND IFWRITE(*,150) NE0,F(1),F(2)150 FORMAT(1X,单元号=,I2,5X,左端弯矩=,F9.3,2X,右端弯矩,F9.3)200 CONTINUEWRITE(*,11)11 FORMAT(1X,= 计算结束 =)END
