- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《桥梁结构振动与抗震》动力电算分析
静力分析程序计算
斜拉桥模型基本信息:
计算模型
结构体系:半漂浮体系;
主梁:单元号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,
文档评论(0)