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

1、静力分析程序计算一、 斜拉桥模型基本信息：计算模型1、 结构体系：半漂浮体系；2、 主梁：单元号1-23（1） 跨径：中跨：100m；边跨：50m；（2） 截面参数：换算截面：A=58.8m2,惯性矩：I=2.20m4,每延米重：BRO=-180.00KN/m3、 塔墩：双塔，单元号24-37（1）索塔高30m，塔跨比为1/3.33，墩高20m（2）截面参数：换算截面：A=12.00m2,惯性矩：I=4.58m4,每延米重：BRO=-335.46KN/m4、 斜拉索：单元号38-57（1）竖琴形单索面形式，索距：主梁8m，索塔6m；（2）截面参数：换算截面：A=0.200m2,每延米重：BRO

2、=-10.00KN/m二、 程序修改部分：（1）结构参数NJ=40,NE=57,N=111,UBW1=12,NB1=9,NJX=120,M=6DATA 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,

3、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,5

4、0,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/（3） 修改对应的截面参数，并前后统一；（4） 子程序1) 平面梁单元和铰接杆单元刚度矩阵BEMK、ROADK：REAL X(40) Y(40)2) 自重荷载作用下的等价集中力BEMLOAD：SUBROUTINE BEMLOD(BRO,L,XM,XN,PF1,PF2,IAA) 其中梁单元IAA=1.0，铰接杆单元IAA=0.0;（5）主程序1）整体刚度矩阵和荷载等价集中

5、力列阵:塔墩：IZ=4.58AX=12.00BRO=-335.46CALL 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.20BRO=-10.0DO 170 I=1,6DO 170 J=1,6KS(I,J)=0.0 170 CONTINUE CALL RODK(IO,JO,K

6、E,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.00IZ=4.58BRO=-335.46CALL 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

7、,PF1,PF2,1)CALL ELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJX)斜拉索：AX=0.200BRO=-10.00 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 ELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJX)三、计算主程序：PLANE ANALYSIS OF Cable-Stayed BRIDGEPROGRAM Cable_Stayed BRI

8、DGEinteger ubw,f,ubw1real l,iz,lx,lyPARAMETER (NJ=40,NE=57,N=111,UBW1=12,NB1=9,NJX=120,M=6)integer io(NE),jo(NE),rl1(NB1),RLS(NJX),CRL(NJX),NIJ(M)REAL KO(N,UBW1),R(M,M),RT(M,M),KE(M,M),KR(M,M),KS(M,M)REAL A(NE,7),XU(N),p(N),xv(NJX),pij(M),pij2(M),fij(M)REAL PF1(M),PF2(M),X(NJ),Y(NJ)CREADING ARRAYSDA

9、TA 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,

10、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,

11、0,0,24,0,0,18,0,0,12,0,0,6,0,0,0,-10,-20/CBEGIN THE CALCULATION DO 10 I=1,NB1J=RL1(I)RLS(J)=1 10CONTINUEWRITE(*,20) 20FORMAT(1X,18H开始计算约束表rls=)WRITE(*,40) RLS 40FORMAT(1X,/24I3)do 60 I=1,njxK=0DO 50 J=1,IK=K+RLS(J) 50CONTINUECRL(I)=K 60CONTINUEWRITE(*,70) 70FORMAT(1X,4HCRL=)WRITE(*,40) CRLWRITE(*,10

12、0)100FORMAT(1X,12H开始计算总刚)F=0UBW=0WW=0YE=2.5E+07110 F=F+1IF(F.GT.23) GOTO 150IZ=2.20AX=63.1BRO=-180.00115 CALL BEMK(io,jo,ke,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)DO 120 I=1,6DO 120 J=1,6KS(I,J)=0.0KS(I,J)=KE(I,J) 120CONTINUECALL BEMLOD(BRO,L,XM,XN,PF1,PF2,1)DO 140 I=1,6PF2(I)=PF1(I) 140CONTINUEGOTO

13、270 150 IF(F.GT.37) GOTO 160IZ=4.58AX=12.00BRO=-335.46CALL 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)GOTO 270160 AX=0.20BRO=-10.0DO 170 I=1,6DO 170 J=1,6KS(I,J)=0.0 1

14、70 CONTINUE 180CALL 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) 270 DO 280 J=1,3NIJ(J)=I1+J-1NIJ(J+3)=J1+J-1 280CONTINUECALL COMBK(NJ,N,UBW,NIJ,CRL,RLS,KO,KS,P,PF2,M,NJX,u

15、bw1)WRITE(*,285) F,UBW 285FORMAT(1X,2HF=,I3,4HUBW=,I3) WW=WW+L*BROIF (F.LT.NE) GOTO 110CALL CHOLS(N,UBW,KO) CALL CHOLEJ(N,UBW,P,XU,KO)CDEVELOPMENT OF UDO 320 I=1,3*NJIF (RLS(I).NE.1) GOTO 310XV(I)=0.0GOTO 320 310XV(I)=XU(I-CRL(I) 320CONTINUEWRITE(*,325) 325FORMAT(25X,14H节点位移xv(i)=)DO 330 I=1,NJWRITE

16、(*,335) I,(XV(3*(I-1)+J),J=1,3) 330CONTINUE 335FORMAT(1X,I3,3(8X,F11.7)CSULUTION OF ELEMENT FORCEWRITE(*,395) 395FORMAT(1X,14H开始计算单元力)f=0 400f=f+1IF(F.GT.23) GOTO 425IZ=2.2AX=58.8BRO=-180.00 402 CALL BEMK(IO,JO,KE,I1,J1,YE,AX,IZ,NE,F,L,LX,LY,X,Y,XM,XN)DO 405 J=1,3PIJ(J)=XV(I1+J-1)PIJ(J+3)=XV(J1+J-1)

17、 405CONTINUEA(F,1)=FCALL BEMLOD(BRO,L,XM,XN,PF1,PF2,1)CALL AAXA(KE,PIJ,FIJ,M)DO 407 I=1,6A(F,I+1)=FIJ(I)-PF1(I) 407 CONTINUEGOTO 445 425IF (F.GT.37) GOTO 437AX=12.00IZ=4.58BRO=-335.46CALL 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)CA

18、LL ELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJX) 420 CONTINUEGOTO 445 437AX=0.200BRO=-10.00 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 ELMFOC(F,NE,R,XV,A,PF1,KE,I1,J1,M,NJX) 445IF (F.LT.NE) GOTO 400WRITE(*,450) 450FORMAT(25X,13H单元力a(I

19、,J)=)DO 460 I=1,NEWRITE(*,465) (A(I,J),J=1,7) 460CONTINUE 465FORMAT(1X,f4.0,3(8X,F15.5)/12X,F15.5,2(8X,F15.5)C OPEN (6,FILE='OUT.txt') WRITE(*,466) WW 466FORMAT(1X,F18.4)END四、运行结果：开始计算约束表rls= 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

20、0 0 0 0 0 0 0 1 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 CRL= 1 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 4 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6

21、 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 7 8 9 开始计算总刚 F= 1UBW= 8 F= 2UBW= 10 F= 3UBW= 10 F= 4UBW= 10 F= 5UBW= 10 F= 6UBW= 10 F= 7UBW= 12 F= 8UBW= 12 F= 9UBW= 12 F= 10UBW= 12 F= 11UBW= 12 F= 12UBW= 12 F= 13UBW= 12 F= 14UBW= 12 F= 15UBW= 12 F= 16UBW= 12 F= 17UBW=

22、12 F= 18UBW= 12 F= 19UBW= 12 F= 20UBW= 12 F= 21UBW= 12 F= 22UBW= 12 F= 23UBW= 12 F= 24UBW= 12 F= 25UBW= 12 F= 26UBW= 12 F= 27UBW= 12 F= 28UBW= 12 F= 29UBW= 12 F= 30UBW= 12 F= 31UBW= 12 F= 32UBW= 12 F= 33UBW= 12 F= 34UBW= 12 F= 35UBW= 12 F= 36UBW= 12 F= 37UBW= 12 F= 38UBW= 12 F= 39UBW= 12 F= 40UBW= 1

23、2 F= 41UBW= 12 F= 42UBW= 12 F= 43UBW= 12 F= 44UBW= 12 F= 45UBW= 12 F= 46UBW= 12 F= 47UBW= 12 F= 48UBW= 12 F= 49UBW= 12 F= 50UBW= 12 F= 51UBW= 12 F= 52UBW= 12 F= 53UBW= 12 F= 54UBW= 12 F= 55UBW= 12 F= 56UBW= 12 F= 57UBW= 12节点位移xv(i)= 1 0.0000000 0.0000000 0.0000000 2 0.0025195 -0.0011915 -0.0005210 3

24、 0.0107661 -0.0022712 -0.0011454 4 -0.0001902 -0.0148611 -0.0028911 5 -0.0001244 0.0065998 -0.0022557 6 0.0188943 -0.0028653 -0.0015702 7 -0.0002560 -0.0384054 -0.0030558 8 -0.0000779 0.0198689 -0.0011090 9 0.0292416 -0.0033239 -0.0018110 10 -0.0003097 -0.0642556 -0.0034072 11 -0.0000411 0.0249372 -

25、0.0001844 12 0.0400289 -0.0036709 -0.0017523 13 -0.0003502 -0.0923223 -0.0035329 14 -0.0000144 0.0230835 0.0006411 15 0.0501347 -0.0039001 -0.0016258 16 -0.0003964 -0.1340584 -0.0022193 17 0.0000002 0.0150259 0.0013366 18 0.0596950 -0.0040129 -0.0015770 19 -0.0003999 -0.1383182 -0.0020371 20 0.00000

26、00 0.0000000 0.0015173 21 -0.0008017 0.0000000 -0.0015156 22 -0.0004010 -0.1382917 0.0020403 23 -0.0605021 -0.0040129 0.0015790 24 -0.0008017 0.0150059 -0.0013340 25 -0.0004044 -0.1340251 0.0022228 26 -0.0509291 -0.0039002 0.0016282 27 -0.0007870 0.0230359 -0.0006368 28 -0.0004506 -0.0922293 0.00353

27、73 29 -0.0408036 -0.0036709 0.0017567 30 -0.0007601 0.0248511 0.0001894 31 -0.0004911 -0.0641349 0.0034088 32 -0.0299759 -0.0033239 0.0018209 33 -0.0007231 0.0197509 0.0011111 34 -0.0005447 -0.0383049 0.0030478 35 -0.0195393 -0.0028653 0.0015913 36 -0.0006764 0.0064991 0.0022479 37 -0.0006106 -0.014

28、8617 0.0028762 38 -0.0112488 -0.0022712 0.0011768 39 -0.0026823 -0.0011915 0.0005494 40 0.0000000 0.0000000 0.0000000 开始计算单元力单元力a(I,J)= 1. -28.93215 400.72119 -0.00171 28.93215 1399.27881 -4986.50977 2. 2877.67090 530.67322 4986.50098 -2877.67090 909.32678 -6495.67188 3. 5270.57715 685.35114 6495.66

29、846 -5270.57715 754.64886 -6771.60107 4. 7256.05859 584.46387 6771.59961 -7256.05859 855.53613 -7859.44580 5. 9164.03125 475.45593 7859.45020 -9164.03125 964.54407 -9825.05273 6. 12978.41113 1846.24072 9825.06543 -12978.41113 -406.24072 -829.28906 7. 12978.40723 406.23758 829.32111 -12978.40723 1033

30、.76245 -3355.76172 8. 10576.90625 717.37811 3355.75830 -10576.90625 722.62189 -3394.49365 9. 7993.31641 1115.07861 3394.44458 -7993.31641 324.92145 -252.87531 10. 5208.30615 1613.81909 252.84082 -5208.30615 906.18091 4691.56934 11. 2689.44531 461.28223 -4691.27490 -2689.44531 -101.28223 5206.42334 1

31、2. 88.86657 1800.45007 -5206.22754 -88.86657 1799.54993 5215.23730 13. 2687.74927 -99.80148 -5214.02832 -2687.74927 459.80148 4700.86865 14. 5205.50391 905.30182 -4701.48584 -5205.50391 1614.69824 -255.23047 15. 7985.83154 320.54352 255.21814 -7985.83154 1119.45654 -3431.77295 16. 10565.60352 715.36

32、249 3431.68530 -10565.60352 724.63751 -3451.03296 17. 12992.79297 1045.78076 3451.05029 -12992.79297 394.21918 -828.50745 18. 12992.79980 -394.21631 828.49707 -12992.79980 1834.21631 -9728.13477 19. 9204.16016 957.28094 9728.12402 -9204.16016 482.71906 -7820.63965 20. 7292.67969 850.87030 7820.63770

33、 -7292.67969 589.12970 -6770.09229 21. 5303.12891 753.01569 6770.08789 -5303.12891 686.98431 -6507.19824 22. 2907.93335 909.41510 6507.19873 -2907.93335 530.58490 -4997.29639 23. -0.00676 1400.35559 4997.29395 0.00676 399.64438 0.00293 24. 37422.12109 -117.78848 5374.62207 -34067.51953 117.78848 -65

34、54.89160 25. 34067.53125 -117.78899 6554.89600 -30712.93164 117.78899 -7740.41943 26. 30712.92578 -117.79594 7740.39893 -28700.16406 117.79594 -8467.93848 27. 23938.27734 1295.09399 8467.90527 -21925.51562 -1295.09399 -723.23242 28. 18356.83203 619.46649 723.27936 -16344.07129 -619.46649 2966.28955

35、29. 12466.17480 -180.00641 -2966.17725 -10453.41504 180.00641 1860.29736 30. 6643.14990 -306.03296 -1860.48303 -4630.39014 306.03296 -0.20698 31. 37422.08203 88.82239 -5843.95215 -34067.48047 -88.82239 6734.69434 32. 34067.47266 88.82223 -6734.69385 -30712.87305 -88.82223 7630.82080 33. 30712.86914

36、88.81074 -7630.82520 -28700.10742 -88.81074 8184.88086 34. 23938.22852 -1272.62048 -8184.82910 -21925.46680 1272.62048 575.22009 35. 18357.04102 -604.34418 -575.13354 -16344.28027 604.34418 -3023.67529 36. 12466.90137 186.46631 3023.63916 -10454.14160 -186.46631 -1878.98108 37. 6642.93604 309.08929

37、1878.95837 -4630.17627 -309.08929 0.08079 38. -3483.25537 200.00000 0.00000 3783.25537 200.00000 0.00000 39. -2871.13330 160.00000 0.00000 3111.13330 160.00000 0.00000 40. -2391.85376 120.00000 0.00000 2571.85376 120.00000 0.00000 41. -2324.97363 80.00000 0.00000 2444.97363 80.00000 0.00000 42. -4737.97363 40.00000 0.00000 4797.97363 40.00000 0.00000 43. -3100.73706 -200.00000 0.00000 3400.73706 -200.00000 0.00000 44. -2859.24854 -190.00000 0.00000 3099.24854 -190.00000 0.00000 45. -3391.25928 -120.00000 0.00000 3571.25928 -120.00000 0.00000 46. -3169.49341 -80.00000 0.00000 3289.4934