基于纤维模型的高墩连续刚构桥地震行波效应分析

Size: px

Start display at page:

Download "基于纤维模型的高墩连续刚构桥地震行波效应分析"

侧蜀文
5 years ago
Views:

1 第 4 卷第 3 期 2009 年 9 月震灾防御技术 Technology for Earthquake Disaster Prevention Vol. 4, No. 3 Sep., 2009 张磊, 韩强, 杜修力,2009. 基于纤维模型的高墩连续刚构桥地震行波效应分析. 4(3): 基于纤维模型的高墩连续刚构桥 1 地震行波效应分析张磊韩强杜修力 ( 北京工业大学城市与工程安全减灾省部共建教育部重点实验室, 北京 ) 摘要本文基于在 MSC.Marc 中加入纤维模型编制的 THUFIBER 程序的基础上, 分别采用纤维模型和三维实体单元模型模拟桥墩和上部结构, 建立高墩大跨连续刚构桥有限元模型同时, 采用 2008 年汶川地震波, 分析了行波效应对四川某高墩大跨连续刚构桥非线性地震反应的影响结果表明 : 行波效应对高墩桥的墩顶位移和墩底影响较大, 墩顶位移减小了 20% 50%; 而墩底增大了 20% 30% 关键词 : 大跨度桥梁纤维模型地震响应行波效应高墩引言我国西部地区多为高地震烈度的山岭重丘区或黄土深壑, 地形地貌和地质条件复杂, 山区桥梁结构通常采用多联连续梁或连续刚构, 下部一般为高墩, 且墩高相差悬殊, 属于典型的非规则桥梁高墩桥梁结构复杂, 多采用薄壁空心墩, 长细比较大, 地震反应较中低桥墩复杂此外, 高墩大跨桥梁空间尺度大, 自振周期长, 在强震作用下由于变形过大而容易发生破坏, 通常表现为桥梁支座系统或局部构件的连接对大变形的不适应性, 缺乏空间变形能力, 特别是单柱式高桥墩情况更为突出但是, 我国现行公路工程抗震设计规范 (JTJ ) ( 中华人民共和国交通部部颁标准,1990) 仅适用于跨径不超过 150m 的钢筋混凝土和预应力钢筋混凝土梁桥圬工或钢筋混凝土拱桥的抗震设计在桥墩墩高方面, 现行规范只对墩高不超过 30m 的桥梁抗震计算和设计作了相关的规定在现行抗震规范中, 除欧洲规范考虑了地震动的空间变化外, 其余规范都采用一致地震动输入地震现场的实测资料表明, 即使在 50m 的范围内, 地基各点振动的幅值和相位也有很大差别 ( 梁爱虎等,1996) 行波法是对大跨度桥梁采取非一致输入地震反应的有效方法, 假定地基条件一致, 地震波沿地表面以一定速度传播, 结构各支承点处地震波波形不变, 只是存在时间滞后和振幅衰减尽管行波法因上述假设存在诸多局限性, 但它仍在一定程度上反映了地震波传播的基本特征 1 基金项目国家自然科学基金项目 ( , ) 资助 [ 收稿日期 ] [ 作者简介 ] 张磊, 男, 生于 1983 年硕士研究生主要从事桥梁抗震方面的研究 curious123@163.com

将会明显影响桥墩的承载力和延性大量已有研究表明, 基于材料本构模型的钢筋混凝土杆件纤维模型, 是分析此类结构非线性问题的较好方法 ( 江见鲸等,2005) 因此, 本文基于在通用有限元软件 MSC.

2 3 期张磊等 : 基于纤维模型的高墩连续刚构桥地震行波效应分析 313 因此, 基于行波法的多点输入在一定程度上反映了地震波的传播效应对大跨度桥梁的地震反应特性的影响本文在通用有限元程序 MSC.Marc 中加入基于纤维模型编制的 THUFIBER 程序的基础上, 采用纤维模型模拟钢筋混凝土高墩, 上部结构则采用实体单元建模, 分析行波效应对高墩连续刚构桥的地震反应影响 1 钢筋混凝土桥墩的纤维模型现阶段, 我国高墩桥梁多为钢筋混凝土桥墩在动力荷载作用下, 桥墩的钢筋可能屈服, 混凝土会开裂甚至会破碎, 且桥墩受到的轴力将会发生变化, 轴力和间复杂的耦合关系, 将会明显影响桥墩的承载力和延性大量已有研究表明, 基于材料本构模型的钢筋混凝土杆件纤维模型, 是分析此类结构非线性问题的较好方法 ( 江见鲸等,2005) 因此, 本文基于在通用有限元软件 MSC.Marc 基础上开发的纤维模型有限元程序 THUFIBER( 叶列平等,2006), 建立桥墩的有限元模型 THUFIBER 程序将构件截面划分成若干混凝土纤维或钢纤维 ( 图 1), 根据截面曲率和应变计算得到每根纤维的应变, 再由应变根据钢筋和混凝土的滞回本构模型 ( 图 2) 得到纤维应力, 最后积分得到截面的和轴力由于采用了比较准确的材料本构模型, 因而可以较好地模拟钢筋混凝土构件在复杂轴力下的强度延性和滞回行为 ( 陆新征等, 2007) 图 1 THUFIBER 程序中截面纤维划分 Fig.1 Fiber partition in THUFIBER program 图 2 THUFIBER 程序中混凝土及钢筋纤维的本构模型 Fig.2 Concrete and reinforced fiber constitutive model in THUFIBER program 2 多点激励的运动方程对于大跨度桥梁结构, 在多点非一致地震激励下, 其运动方程与一致地震激励下的运动方程不同, 其表达式为 : Mss Msb X&& s Css Csb X& s Kss Ksb Xs = M bs Mbb X&& C b bs Cbb X& K b bs Kbb X b Fb (1) 式中,{X b } { X & b } { X && b } 和 {X s } { X & s } { X && s } 分别表示在地震激励下桥梁结构支座节点和非支座节点的绝对位移速度和加速度向量 ;{F b } 表示地震动作用于桥梁结构支座节点的力

3 314 震灾防御技术 4 卷展开式 (1) 的第一行可得 : [M ss ]{ X && s } +[C ss ]{ X & s } +[K ss ]{X s }= [M sb ]{ X && b } [C sb ]{ X & b } [K sb ]{X b } (2) 为求解方程 (1), 通常把结构的绝对位移分解成拟静力位移和动相对位移, 其中动相对位移由结构的动力响应引起, 而拟静力位移由结构支承点处的位移引起, 即 : Xs Xm X = + (3) Xb Xb 0 式中,{X m } 表示因支承运动 {X b } 所引起的自由节点处的拟静力位移 ;{X} 为自由节点相对于地面的动力位移由于 {K sb } 表示因支承单位位移在自由节点上产生的力, 而 {K ss } 表示自由节点单位位移所产生的力, 因此 {X m } 和 {X b } 满足条件由此得到 {X m } 和 {X b } 的关系为 : [K ss ]{X m }+[K sb ]{X b }={0} (4) {X m }= [K ss ] 1 [K sb ]{X b } (5) 将式 (3) 第一个方程代入式 (2), 则 [M ss ]{X}+[C ss ]{ X & } +[K ss ]{X}={P(t)} (6) 其中,{P(t)} 表示等效地震作用向量, 它的表达式为 : {P(t)}= ([M ss ]{ X && m} +[M sb ]{ X && b }) ([C ss ]{ X & m} +[C sb ]{ X & b} )-([K ss ]{X m }+[K sb ](X b )) (7) 式 (7) 右端第二项表示结构与支座的阻尼耦联, 由于比较小, 通常可以忽略同时据式 (4) 和式 (5), 则式 (7) 可简化为 : {P(t)}=([M ss ][K ss ] 1 [K sb ] [M sb ]){ X && b } (8) 3 实例分析本算例以位于北京至昆明高速公路四川境内雅安至泸沽段栗子树大桥中的一联连续刚构为原型, 适当修改作为计算模型 ( 梁智垚,2007) 全桥孔跨径布置为 50m+2 80m+50m, 其计算简图见图 3; 主梁为单箱单室变截面预应力混凝土箱梁, 其截面形式见图 4; 桥墩为薄壁空心方形墩, 其截面形式见图 5; 不同墩高的连续刚构墩柱基本参数见表 1 本文在进行多点激励与行波效应分析时, 不考虑基础土 - 结构相互作用, 采用大质量法来模拟行波效应, 即在桥梁结构的支座节点上, 建立一个单元点质量单元, 给这一节点单元设定很大的质量 ( 通常为结构总质量的 100 万倍以上 ) 对结构进行动力分析时, 释放支座节点相应激励方向的自由度约束, 然后在支座节点激励方向施加很大的力, 这个力为质量单元的质量与地面加速度的乘积通过这种方法, 可间接地把地面加速度施加到支座节点上在通用有限元软件 MSC.Marc 进行数值模拟时, 高墩采用纤维模型来模拟, 主梁采用实体单元来模拟, 主梁与桥墩之间的连接采用软件中的 LINK100 单元 ( 可以约束节点的所有自由度 ), 而主梁与桥墩的连接和主梁模型如图 6 所示全桥的有限元模型如图 7 所示

4 3 期张磊等 : 基于纤维模型的高墩连续刚构桥地震行波效应分析 315 图 3 桥梁计算图 ( 单位 :m) Fig.3 Sketch of bridge for calculation (unit: m) 图 4 主梁横截面图 ( 单位 :cm) Fig.4 Cross-section of girder (unit: cm) 图 5 桥墩横截面图 Fig.5 Cross-section of pier 表 1 墩柱基本参数 Table 1 Basic parameters of piers 30m 墩 60m 墩 90m 墩边长 a(m) 边长 b(m) 厚度 d(%) 配筋率 (%) 全桥结构自振频率和振型全桥自振频率和振型如图 8 所示, 括号中数值为该阶振型对应的自振频率值由图 8 可知, 本桥的第 1 阶振型为体系纵向飘移, 基频较小, 为 0.46 Hz, 是第 2 阶振型的 63%, 第 4 阶振型的 27%, 该振型在所有振型中占有比较大的优势 ; 该桥第 2 到第 4 阶振型均为横弯, 横弯振型对高墩内力反应起比较大的作用 3.2 地震波选取与工况考虑到行波效应对高墩桥梁的地震反应影响, 以及桥梁所处地质地形条件等, 选用了 2008 年汶川地震波 ( 地点 :51JYC), 图 9 为其加速度时程和 Fourier 谱曲线 2008 年汶川地震波记录时间为 250s, 为了简化和减少计算时间, 同时也为了便于对比分析, 取中间的第 40s 至第 140s, 共 100s 进行计算, 时间间隔 0.005s, 分以下 4 种工况对行波效应进行分析

5 316 震灾防御技术 4 卷图 6 主梁与桥墩的连接及主梁模型 Fig.6 Link between main beam and pier, and main beam model 图 7 全桥的有限元模型 Fig.7 Finite element model of the bridge (a) 第 1 阶振型 (0.461Hz) (b) 第 2 阶振型 (0.735Hz) (c) 第 3 阶振型 (1.328Hz) (d) 第 4 阶振型 (1.734Hz) Fig.8 图 8 前四阶振型 The first 4 modal at 0.461Hz, 0.735Hz, 1.328Hz and 1.734Hz

6 3 期张磊等 : 基于纤维模型的高墩连续刚构桥地震行波效应分析 317 图 9 汶川波的加速度时程曲线和 Fourier 谱 Fig.9 Time-history acceleration curve and Fourier spectrum of Wenchuan seismic waves 工况一 :P1 墩滞后 0.1s,P2 墩滞后 0.26s,P3 墩滞后 0.42s,P4 墩滞后 0.52s; 工况二 :P1 墩滞后 0.06s,P2 墩滞后 0.17s,P3 墩滞后 0.28s,P4 墩滞后 0.34s; 工况三 :P1 墩滞后 0.05s,P2 墩滞后 0.13s,P3 墩滞后 0.21s,P4 墩滞后 0.26s; 工况四 : 各墩一致输入 3.3 考虑行波效应影响的地震反应姜海鹏等 (2001) 郑史雄等(1997) 和马坤全 (1997) 的研究表明, 高墩刚构桥在地震反应中上部结构受力较小, 结构内力反应最大值基本都出现在墩底与墩顶处, 桥墩的动力响应将对整个桥梁工程抗震起到关键的作用因此本文对桥墩墩顶相对位移以及桥墩顶与墩底内力进行了分析 (1) 墩顶纵向位移分别计算了一致输入下和考虑行波效应下墩顶的纵向位移时程, 中间桥墩的位移时程见图 10 计算结果表明: 一致激励 ( 工况四 ) 时, 各桥墩墩顶位移大致相等, 这与主梁的横向刚度较大有关 ; 无论一致激励还是考虑行波效应, 可以很明显的看出纵向位移的衰减 ; 当水平视波速为 1000m/s 时 ( 工况三 ), 各墩墩顶的纵向位移时程与一致激励 ( 工况四 ) 下大致相等 ; 但当水平视波速为 500m/s( 工况一 ) 和 750m/s( 工况二 ) 时, 各墩墩顶的纵向位移与一致激励相比有较大差别其中工况一各墩墩顶的最大纵向位移分别为工况四的 60% 58% 56% 52% 49%; 工况二各墩墩顶的最大纵向位移分别为工况四的 79% 80% 83% 82% 81% 可见, 当地震波水平视波速较小时, 行波效应对于桥梁桥墩的纵向位移的影响是有利的 (2) 桥墩内力分别计算了一致输入下和考虑行波效应时桥墩的受力情况, 由于高墩桥梁多为弯曲破坏和剪切破坏, 因此, 抗震设计时更加关注桥墩所受和的最大值, 表 2 和表 3 分别列出了各个桥墩墩顶与墩底所受和最大值 ( 绝对值 ) 通过以上分析结果可以看出, 本桥在综合考虑一致激励和行波效应后, 主墩最大值和最大值都发生在 P2 墩墩底, 也即最高墩的墩底, 其中最大为 kn m, 最大值为 kn; 同时桥墩墩底所受的和都大于墩顶的相关值, 可见桥墩的墩底是地震分析时的危险截面, 也即控制截面 ;P0 和 P4 墩顶和随着水平视波速的增加而逐渐减小, 其余各墩墩顶和随着水平视波速的增加反而逐渐增大 ; 而墩底的除最高墩 P2 墩以外, 都随着水平视波速的增加而增大 ; 各墩墩底的随着水平视波速的增加却逐渐减小

result of the top-section of piers P0 墩顶 P1 墩顶 P2 墩顶 P3 墩顶 P4 墩顶工况工况一 32428.80 3247.

50 49631.10 4932.20 81968.10 8139.19 42873.70 4289.34 30632.30 3011.75 工况三 30370.20 3087.

7 318 震灾防御技术 4 卷图 10 P2 墩墩顶的纵向位移时程 Fig.10 Horizontal displacement of the top of P2 pier Table 2 表 2 主墩墩顶截面地震反应结果 Seismic response result of the top-section of piers P0 墩顶 P1 墩顶 P2 墩顶 P3 墩顶 P4 墩顶工况工况一工况二工况三工况四

8 3 期张磊等 : 基于纤维模型的高墩连续刚构桥地震行波效应分析 319 Table 3 表 3 主墩墩底截面地震反应结果 Seismic response result of the bottom-section of piers 工况 P0 墩底 P1 墩底 P2 墩底 P3 墩底 P4 墩底工况一工况二工况三工况四结论本文基于纤维模型编制的 THUFIBER 程序基础上, 在通用有限元程序 MSC.Marc 中加入采用纤维模型和三维实体单元模型, 分别模拟钢筋混凝土桥墩和桥梁上部结构, 建立了高墩大跨连续刚构桥的有限元模型并采用汶川地震记录, 针对四川境内某一实际桥梁, 分析了行波效应对其地震反应的影响, 获得了以下结论 : (1) 地震动行波效应对算例桥梁的位移影响较大水平视波速为 500m/s 时, 各墩墩顶的最大纵向位移分别为一致激励时的 60% 58% 56% 52% 49%; 水平视波速为 750m/s 时, 各墩墩顶的最大纵向位移分别为工况四时的 79% 80% 83% 82% 81% (2) 在算例中, 最高墩 P2 墩底的和各墩墩底的随着水平视波速的增加逐渐减小由于目前工程实际中使用一致激励输入, 所以通常的设计方法是偏于不安全的为此, 在进行高墩桥梁抗震设计时, 建议考虑行波效应的影响 (3) 考虑地震动的行波效应时, 与一致激励比较, 结构的动力反应在某些区段某些截面上可能加大, 也可能降低, 这种情况在本算例主墩的和中都有体现, 说明地震动的行波效应对结构反应较为复杂, 其影响与结构的特性计算峰值和截面位置等因素有关参考文献江见鲸, 陆新征, 叶列平,2005. 混凝土结构有限元分析. 北京 : 清华大学出版社. 姜海鹏, 王航, 杨庆山,2001. 空间变异地震作用下大型连续刚构桥的延性性能研究. 铁道工程学报,6(2): 梁爱虎, 陈厚群, 侯顺载,1996. 随机地震动场激励下拱坝多点输入的抗震可靠度分析. 地震工程与工程振动,16(1): 梁智垚,2007. 非规则高墩桥梁抗震设计理论研究. 上海 : 同济大学. 陆新征, 叶列平等,2007. 考虑地震行波效应大型高架桥梁破坏模拟. 工程抗震与加固改造,29(3):1 5. 马坤全,1997. 连续刚架桥抗震延性分析. 上海铁道大学学报 ( 自然科学版 ),18(4):6 16. 叶列平,Qazi A.U., 马千里, 陆新征,2006. 高强钢筋对框架结构抗震破坏机制和性能控制的研究. 工程抗震与加固改造,28(1): 郑史雄, 奚绍中, 杨建忠,1997. 大跨度刚构桥的地震反应分析. 西南交通大学学报,32(6): 中华人民共和国交通部部颁标准,1990. 公路工程抗震设计规范 (JTJ004-89). 北京 : 人民交通出版社.

9 320 震灾防御技术 4 卷 Seismic Traveling Wave Effect of Continuous Rigid-frame Bridge with High Piers Base on Fiber Model Zhang Lei, Han Qiang and Du Xiuli (The Key Lab. of Urban Security and Disaster Engineering of Ministry of Education, Beijing University of Technology, Beijing , China) Abstract In this paper, with application of THUFIBER program, fiber model is used to simulate bridge piers and solid element is used to simulate the superstructure in program MSC.Marc. The effects of traveling wave on the nonlinear seismic response of long-span continuous rigid-frame bridge with high piers are analyzed under 2008 Wenchuan earthquake wave excitation. Out results show that traveling wave has a great effect on displacement and shear force of high piers of the bridge. The displacement of top of pier decreases 20% to 50% and shear force at base of pier increases 20% to 30%. Key words: Long-span bridge; Fiber model; Seismic response; Traveling wave effect; High pier

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva 39 6 2011 12 Journal of Fuzhou University Natural Science Edition Vol 39 No 6 Dec 2011 DOI CNKI 35-1117 /N 20111220 0901 002 1000-2243 2011 06-0923 - 07 350108 105 m 14 69% TU311 3 A Seismic analysis of