結(jié)力電算程序連續(xù)梁計(jì)算應(yīng)用報(bào)告

NNYNYNRW=(1)求截面內(nèi)力(2)求影響線(3)種：12341（均布荷載左端距本單元左端點(diǎn)距離，對(duì)于其他荷載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) 荷載大小 非結(jié)點(diǎn)荷載:PF(NE,1)=荷載大小PF(NE,2)=距離1PF(NE,3)=距離2 PF(NE,4)=單元號(hào)PF(NE,5)=荷載類型號(hào)FORMAT(1X,'桿長,慣性矩FORMAT(1X,'節(jié)點(diǎn)荷載大小,對(duì)應(yīng)未知數(shù)序號(hào)PJ(I,1),PJ(I,2)')ENDIF ENDIF 計(jì)算第NHF個(gè)非結(jié)點(diǎn)荷截引起的等效結(jié)點(diǎn)荷截列陣F0 DIMENSIONPF(NF,5),F0(2),GC(NE) GOTO200 GOTO200 GOTO200 計(jì)算第NE0個(gè)單元的單元?jiǎng)偠染仃嘍K(2,2) SUBROUTINEDG(NE0,NE,E0,GC,GX,DK)DIMENSIONGC(NE),GX(NE),DK(2,2)DO15DO10J=1,2 集成總體剛度矩陣,ZKl,ZK2分別存主對(duì)角元素和付對(duì)角元素 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計(jì)算單元桿端內(nèi)力（桿端彎矩，桿端剪力，支座反力）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 彎矩包絡(luò)圖值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 五、算例結(jié)利用程序計(jì)算下面例題（E0=5,GX=2計(jì)算截面內(nèi)計(jì)算第二跨跨中彎矩剪421250031300632單元數(shù) 支承類型 1非節(jié)點(diǎn)荷載個(gè)數(shù) 彈性模量 5.0000截面數(shù) 桿長,慣性矩節(jié)點(diǎn)荷載大小,對(duì)應(yīng)未知數(shù)序號(hào) 位結(jié)點(diǎn)號(hào)=結(jié)點(diǎn)號(hào)=結(jié)點(diǎn)號(hào)=結(jié)點(diǎn)號(hào)=結(jié)點(diǎn)號(hào)=單元號(hào)= 右端彎矩單元號(hào)=單元號(hào)=右端彎矩單元號(hào)=右端彎矩截面編號(hào)，截面所在單元編號(hào)，截面距離，截面彎矩，截面 計(jì)算影計(jì)算第二跨跨中彎矩剪力影響線以及第三個(gè)支座影響42626262630303033 桿長,慣性矩GC(NE),GX(NE) 彎