结力电算程序连续梁计算应用报告

NNYNYNRW=(1)求截面内力(2)求影响线(3)种：12341（均布荷载左端距本单元左端点距离，对于其他荷载PROGRAMDIMENSIONCALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,NJ,ZK1,ZK2,P)CALLFCQJ(NJ,ZK1,ZK2,P)CALLSCWY(NJ,P)CALLSCGDL(NE,GDM)CALLSCJMNL(FJM,FJMNL,NJM)END END SRSJ子程序 SUBROUTINESRSJ(NE,NP,NF,GC,GX,PJ,PF,FJM,NJM) 荷载大小 非结点荷载:PF(NE,1)=荷载大小PF(NE,2)=距离1PF(NE,3)=距离2 PF(NE,4)=单元号PF(NE,5)=荷载类型号FORMAT(1X,'杆长,惯性矩FORMAT(1X,'节点荷载大小,对应未知数序号PJ(I,1),PJ(I,2)')ENDIF ENDIF 计算第NHF个非结点荷截引起的等效结点荷截列阵F0 DIMENSIONPF(NF,5),F0(2),GC(NE) GOTO200 GOTO200 GOTO200 计算第NE0个单元的单元刚度矩阵DK(2,2) SUBROUTINEDG(NE0,NE,E0,GC,GX,DK)DIMENSIONGC(NE),GX(NE),DK(2,2)DO15DO10J=1,2 集成总体刚度矩阵,ZKl,ZK2分别存主对角元素和付对角元素 SUBROUTINEJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)DIMENSIONZK1(NJ),ZK2(NJ),DK(2,2),GC(NE),GX(NE)DO10I=1,NJ DO CALLDG(I,NE,E0,GC,GX,DK) SUBROUTINEXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)DO10I=1,NJ DO ENDIFIF(NF.GT.0)THENDO100NF0=1,NF ENDIF SUBROUTINEZCCL(NZ,NJ,ZK1,ZK2,P) P SUBROUTINEFCQJ(NJ,ZK1,ZK2,P)DIMENSIONZK1(NJ),ZK2(NJ),P(NJ)DO10I=1,NJ-1 DO20I=1,NJ-1 SUBROUTINESCWY(NJ,P)DIMENSIONP(NJ) MX, 移 c************************************************c计算单元杆端内力（杆端弯矩，杆端剪力，支座反力）c************************************************DIMENSIONP(NJ),WY(2),F(2),F0(2),DK(2,2),PF(NF,5),GDM(NE,2)DIMENSIONGC(NE),GX(NE),FZZ(NJ),FY(NE,2)DO200CALLDG(NE0,NE,E0,GC,GX,DK)DO35I=1,2DO30 DO50I=1,NFCALLDJH(I,NE,NF,PF,F0,GC)DO40J=1,2 ENDIF ENDIFDO60ELSEENDIFELSEELSEIF(ID.EQ.4)THENELSEENDIFEND ENDIF DO65J=1,NFENDEND DO70DO80J=1,NFENDIFEND SUBROUTINESCGDL(NE,GDM)DIMENSIONGDM(NE,2) DO10 SUBROUTINEJMNLA(FJM,FJMNL,NJM,FY,GDM,GC,PF,NE,NF)DO100J=1,NJMDO20K=1,NFELSEELSEIF(R.GT.C)THENENDIFELSEELSEIF(R.GT.C)THENENDIFENDELSEIF(ID==3)THENELSEIF(R.GT.C)THENENDELSEIF(ID==4)THENELSEIF(R.GT.C)THENENDELSEIF(ID==5)THENELSEENDIFENDIFENDIF ENDIF100CONTINUE DIMENSIONFJM(NJM,2),FJMNL(NJM,2)1力 DIMENSIONGC(NE),GX(NE),FJM(1,2),PF(1,5),NYXX(NE),FJMNL(1,2), FORMAT(1X,'杆长,惯性矩 FORMAT(1X,'DO400DO300I=0,NYXX(J)DO19K=1,J-119ENDCALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,NJ,ZK1,ZK2,P)CALLCALLJMNLA(FJM,FJMNL,NJM,FY,GDM,GC,PF,NE,NF) 300 弯矩包络图值DIMENSION FORMAT(1X,'杆长,惯性矩 DO6K=1,NEEND CALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,NJ,ZK1,ZK2,P)CALLJMNLA(FJM,FJMNL,NJM,FY,GDM,GC,PF,NE,NF)DO11 DO8K=1,NEEND CALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,NJ,ZK1,ZK2,P)CALLFCQJ(NJ,ZK1,ZK2,P)CALLJMNLA(FJM,FJMNL,NJM,FY,GDM,GC,PF,NE,NF)DO21 DO30J=1,NJM DO100DO40K=1,NE DO41I=NT+2,NE,2 DO42I=NT-1,1,2 CALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,NJ,ZK1,ZK2,P)CALLJMNLA(FJM,FJMNL,NJM,FY,GDM,GC,PF,NE,NF)DO44K=1,NE DO45I=NT+1,NE,2 DO46I=NT-2,1,2 CALLXCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC)CALLJCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX)CALLZCCL(NZ,NJ,ZK1,ZK2,P)CALLJMNLA(FJM,FJMNL,NJM,FY,GDM,GC,PF,NE,NF) DO50K=1,NJM 五、算例结利用程序计算下面例题（E0=5,GX=2计算截面内计算第二跨跨中弯矩剪421250031300632单元数 支承类型 1非节点荷载个数 弹性模量 5.0000截面数 杆长,惯性矩节点荷载大小,对应未知数序号 位结点号=结点号=结点号=结点号=结点号=单元号= 右端弯矩单元号=单元号=右端弯矩单元号=右端弯矩截面编号，截面所在单元编号，截面距离，截面弯矩，截面 计算影计算第二跨跨中弯矩剪力影响线以及第三个支座影响42626262630303033 杆长,惯性矩GC(NE),GX(NE) 弯