《桥梁结构振动与抗震》动力电算分析.doc

静力分析程序计算 斜拉桥模型基本信息： 计算模型 结构体系：半漂浮体系； 主梁：单元号1-23 跨径：中跨：100m；边跨：50m； 截面参数： 换算截面：A=58.8m2, 惯性矩：I=2.20m4, 每延米重：BRO=-180.00KN/m 塔墩：双塔，单元号24-37 （1）索塔高30m，塔跨比为1/3.33，墩高20m （2）截面参数： 换算截面：A=12.00m2, 惯性矩：I=4.58m4, 每延米重：BRO=-335.46KN/m 斜拉索：单元号38-57 （1）竖琴形单索面形式，索距：主梁8m，索塔6m； （2）截面参数： 换算截面：A=0.200m2, 每延米重：BRO=-10.00KN/m 程序修改部分： （1）结构参数 NJ=40,NE=57,N=111,UBW1=12,NB1=9,NJX=120,M=6 DATA IO/20,17,14,11,8,5,4,7,10,13,16,19,22,25,28,31,34,37,36,33, $ 30,27,24,1,2,3,6,9,12,15,40,39,38,35,32,29,26, $ 17,14,11,8,5,19,16,13,10,7,24,27,30,33,36,22,25,28,31,34/ DATA JO/17,14,11,8,5,4,7,10,13,16,19,22,25,28,31,34,37,36,33,30, $ 27,24,21,2,3,6,9,12,15,18,39,38,35,32,29,26,23, $ 18,15,12,9,6,18,15,12,9,6,23,26,29,32,35,23,26,29,32,35/ DATA RL1/1,2,3,58,59,62,118,119,120/ DATA X/-50,-50,-50,-50,-58,-50,-42,-66,-50,-34,-74,-50,-26, $ -82,-50,-18,-90,-50,-10,-100,100,10,50,90,18,50,82, $ 26,50,74,34,50,66,42,50,58,50,50,50,50/ DATA Y/-20,-10,0,0,0,6,0,0,12,0,0,18,0,0,24,0,0,30,0,0,0, $ 0,30,0,0,24,0,0,18,0,0,12,0,0,6,0,0,0,-10,-20/ 修改对应的截面参数，并前后统一； 子程序 平面梁单元和铰接杆单元刚度矩阵BEMK、ROADK：REAL X(40) Y(40) 自重荷载作用下的等价集中力BEMLOAD： SUBROUTINE BEMLOD(BRO,L,XM,XN,PF1,PF2,IAA) 其中梁单元IAA=1.0，铰接杆单元IAA=0.0; （5）主程序 1）整体刚度矩阵和荷载等价集中力列阵: 塔墩：IZ=4.58 AX=12.00 BRO=-335.46 CALL BEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN) CALL BEMR(LX,LY,L,R,RT) CALL BEMLOD(BRO,L,XM,XN,PF1,PF2,1) CALL AAXA(RT,PF1,PF2,M) CALL AAXAA(KE,R,KR,M) CALL AAXAA(RT,KR,KS,M) 斜拉索：AX=0.20 BRO=-10.0 DO 170 I=1,6 DO 170 J=1,6 KS(I,J)=0.0 170 CONTINUE CALL RODK(IO,JO,KE,I1,J1,YE,AX,NE,F,L,LX,LY,X,Y,XM,XN) CALL BEMR(LX,LY,L,R,RT) CALL BEMLOD(BRO,L,XM,XN,PF1,PF2,0) CALL AAXA(RT,PF1,PF2,M) CALL AAXAA(KE,R,KR,M) CALL AAXAA(RT,KR,KS,M) 2）计算单元力 塔墩：AX=12.00 IZ=4.58 BRO=-335.46 CALL BEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN) CALL BEMR(LX,LY,L,R,RT) CALL BEMLOD(BRO,L,XM,XN,PF1,PF2,1) CALL ELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJ