《桥梁结构振动与抗震》静力电算分析报告

《桥梁结构振动与抗震》课程电算程序结果报告书（一）指导教师姓名学号系、年级一、斜拉桥模型基本信息：1、结构体系：漂浮体系；2、主梁：（1）跨径：中跨：100m；边跨：50m；（2）截面参数：换算截面：A=6.31m2,惯性矩：I=2.20m4,每延米重：BRO=-180.00KN/m3、塔墩：双塔（1）索塔高40m，塔跨比为2/5，墩高16m（2）截面参数：换算截面：A=12.00m2,惯性矩：I=4.58m4,每延米重：BRO=-335.46KN/m4、斜拉索（1）竖琴形单索面形式，索距：主梁10m，索塔8m；（2）截面参数：换算截面：A=0.200m2,每延米重：BRO=-10.00KN/m二、程序修改部分：（1）结构参数NJ=34,NE=47,N=94,UBW1=12,NB1=8,NJX=102,M=6DATAIO/17,15,12,9,6,4,7,10,13,16,19,22,25,28,31,29,26,23,20,1,2,3,5,8,11,34,33,32,30,27,24,6,9,12,15,5,8,11,14,28,25,22,19,30,27,24,21/DATAJO/15,12,9,6,4,7,10,13,16,19,22,25,28,31,29,26,23,20,18,2,3,5,8,11,14,33,32,30,27,24,21,5,8,11,14,7,10,13,16,30,27,24,21,29,26,23,20/DATARL1/1,2,3,50,53,100,101,102/DATAX/-50,-50,-50,-50,-50,-60,-40,-50,-70,-30,-50,-80,-20,-50,-90,-10,-100,100,10,90,50,20,80,50,30,70,50,40,60,50,50,50,50,50/DATAY/-16,-8,0,0,10,0,0,20,0,0,30,0,0,40,0,0,0,0,0,0,40,0,0,30,0,0,20,0,0,10,0,0,-8,-16/（2）单元编号主梁：○1~○19；塔墩：○20~○31；斜拉索：○32~○47；如下图所示。（3）修改对应的截面参数，并前后统一；（4）子程序1)平面梁单元和铰接杆单元刚度矩阵BEMK、ROADK：REALX(34)Y(34)2)自重荷载作用下的等价集中力BEMLOAD：SUBROUTINEBEMLOD(BRO,L,XM,XN,PF1,PF2,IAA)其中梁单元IAA=1.0，铰接杆单元IAA=0.0;（5）主程序1）整体刚度矩阵和荷载等价集中力列阵:塔墩：IZ=4.58AX=12.00BRO=-335.46CALLBEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,1)CALLAAXA(RT,PF1,PF2,M)CALLAAXAA(KE,R,KR,M)CALLAAXAA(RT,KR,KS,M)斜拉索：AX=0.20BRO=-10.0DO170I=1,6DO170J=1,6KS(I,J)=0.0170CONTINUECALLRODK(IO,JO,KE,I1,J1,YE,AX,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,0)CALLAAXA(RT,PF1,PF2,M)CALLAAXAA(KE,R,KR,M)CALLAAXAA(RT,KR,KS,M)2）计算单元力塔墩：IZ=4.58AX=12.00BRO=-335.46CALLBEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,1)CALLELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJX)斜拉索：AX=0.200BRO=-10.00CALLRODK(IO,JO,KE,I1,J1,YE,AX,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,0)CALLELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJX)三、计算主程序：PLANEANALYSISOFCable-StayedBRIDGEcPLANEANALYSISOFCable-StayedBRIDGEPROGRAMCable_StayedBRIDGEintegerubw,f,ubw1reall,iz,lx,lyPARAMETER(NJ=34,NE=47,N=94,UBW1=12,NB1=8,NJX=102,M=6)integerio(NE),jo(NE),rl1(NB1),RLS(NJX),CRL(NJX),NIJ(M)REALKO(N,UBW1),R(M,M),RT(M,M),KE(M,M),KR(M,M),KS(M,M)REALA(NE,7),XU(N),p(N),xv(NJX),pij(M),pij2(M),fij(M)REALPF1(M),PF2(M),X(NJ),Y(NJ)CREADINGARRAYSDATAIO/17,15,12,9,6,4,7,10,13,16,19,22,25,28,31,29,26,23,20,$1,2,3,5,8,11,34,33,32,30,27,24,$6,9,12,15,5,8,11,14,28,25,22,19,30,27,24,21/DATAJO/15,12,9,6,4,7,10,13,16,19,22,25,28,31,29,26,23,20,18,$2,3,5,8,11,14,33,32,30,27,24,21,$5,8,11,14,7,10,13,16,30,27,24,21,29,26,23,20/DATARL1/1,2,3,50,53,100,101,102/DATAX/-50,-50,-50,-50,-50,-60,-40,-50,-70,-30,-50,-80,-20,$-50,-90,-10,-100,100,10,90,50,20,80,50,30,70,50,40,60,50,$50,50,50,50/DATAY/-16,-8,0,0,10,0,0,20,0,0,30,0,0,40,0,0,$0,0,0,0,40,0,0,30,0,0,20,0,0,10,0,0,-8,-16/CBEGINTHECALCULATIONopen(5,file='uu.dat',status='new')DO10I=1,NB1J=RL1(I)RLS(J)=110CONTINUEWRITE(5,20)20FORMAT(1X,18H开始计算约束表rls=)WRITE(5,40)RLS40FORMAT(1X,/24I3)do60I=1,njxK=0DO50J=1,IK=K+RLS(J)50CONTINUECRL(I)=K60CONTINUEWRITE(5,70)70FORMAT(1X,4HCRL=)WRITE(5,40)CRLWRITE(5,100)100FORMAT(1X,12H开始计算总刚)F=0UBW=0WW=0YE=2.5E+07110F=F+1IF(F.GT.19)GOTO150IZ=2.20AX=63.1BRO=-180.00115CALLBEMK(io,jo,ke,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)DO120I=1,6DO120J=1,6KS(I,J)=0.0KS(I,J)=KE(I,J)120CONTINUECALLBEMLOD(BRO,L,XM,XN,PF1,PF2,1)DO140I=1,6PF2(I)=PF1(I)140CONTINUEGOTO270150IF(F.GT.31)GOTO160IZ=4.58AX=12.00BRO=-335.46CALLBEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,1)CALLAAXA(RT,PF1,PF2,M)CALLAAXAA(KE,R,KR,M)CALLAAXAA(RT,KR,KS,M)GOTO270160AX=0.20BRO=-10.0DO170I=1,6DO170J=1,6KS(I,J)=0.0170CONTINUE180CALLRODK(IO,JO,KE,I1,J1,YE,AX,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,0)CALLAAXA(RT,PF1,PF2,M)CALLAAXAA(KE,R,KR,M)CALLAAXAA(RT,KR,KS,M)270DO280J=1,3NIJ(J)=I1+J-1NIJ(J+3)=J1+J-1280CONTINUECALLCOMBK(NJ,N,UBW,NIJ,CRL,RLS,KO,KS,P,PF2,M,NJX,ubw1)WRITE(5,285)F,UBW285FORMAT(1X,2HF=,I3,4HUBW=,I3)WW=WW+L*BROIF(F.LT.NE)GOTO110CALLCHOLS(N,UBW,KO)CALLCHOLEJ(N,UBW,P,XU,KO)CDEVELOPMENTOFUDO320I=1,3*NJIF(RLS(I).NE.1)GOTO310XV(I)=0.0GOTO320310XV(I)=XU(I-CRL(I))320CONTINUEWRITE(5,325)325FORMAT(25X,14H节点位移xv(i)=)DO330I=1,NJWRITE(5,335)I,(XV(3*(I-1)+J),J=1,3)330CONTINUE335FORMAT(1X,I3,3(8X,F11.7))CSULUTIONOFELEMENTFORCEWRITE(5,395)395FORMAT(1X,14H开始计算单元力)f=0400f=f+1IF(F.GT.19)GOTO425IZ=2.2AX=6.31BRO=-180.00402CALLBEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)DO405J=1,3PIJ(J)=XV(I1+J-1)PIJ(J+3)=XV(J1+J-1)405CONTINUEA(F,1)=FCALLBEMLOD(BRO,L,XM,XN,PF1,PF2,1)CALLAAXA(KE,PIJ,FIJ,M)DO407I=1,6A(F,I+1)=FIJ(I)-PF1(I)407CONTINUEGOTO445425IF(F.GT.31)GOTO437AX=12.00IZ=4.58BRO=-335.46CALLBEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,1)CALLELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJX)420CONTINUEGOTO445437AX=0.200BRO=-10.00CALLRODK(IO,JO,KE,I1,J1,YE,AX,NE,F,L,LX,LY,X,Y,XM,XN)CALLBEMR(LX,LY,L,R,RT)CALLBEMLOD(BRO,L,XM,XN,PF1,PF2,0)CALLELMFOC(F,NE,R,XV