Analysis of Seismic Response of Bridge across Earthquake Fault with Different Input Modes of Seismic Action
-
摘要: 大地震在近断层场地产生强烈地震动的同时,还会由于断层错动直接导致基岩甚至上覆土层破裂,在断层两侧产生显著差异性永久位移,造成位于断层附近或跨越断层的工程结构破坏。因此,跨断层桥梁面对的地震作用是断层两侧桥墩处场地的不同地震动,包括存在永久性位移的地震动。本文以垂直跨越走滑断层的多跨简支梁桥为例,基于OpenSees有限元模拟平台建立了桥梁结构的三维计算模型,计算分析了不同地震作用输入模式下桥梁结构的地震反应及其差异。考虑的地震作用模式包括:(1)断层两侧场地的地震作用视为相同的无永久位移的地震动,即无永久位移的一致地震动作用模式;(2)断层主动盘一侧场地的地震作用具有永久位移地震动,被动盘一侧采用无永久位移地震动,即具有永久位移的非一致地震动作用模式;(3)在断层主动盘一侧场地以静力方式施加断层错动位移,而被动盘一侧场地固定不动,即断层错动位移静力作用模式。计算结果分析表明,不考虑永久位移的一致地震动作用模式的地震动输入会导致严重低估桥梁反应计算结果,这也说明地震动的断层两侧永久性位移差异会显著增大桥梁结构反应;而一致地震动作用叠加断层错动永久位移静力作用的结果与非一致地震动作用模式的结果非常接近。为此,在某种程度上说,跨断层桥梁结构地震反应可采用一致地震动作用叠加断层错动位移静力作用的桥梁结构反应来近似模拟。Abstract: A large earthquake can produce strong ground motion at the near-fault site, and at the same time, it may directly break up the bedrock and even the overlying soil layers due to fault rupture, which leads to significant differential permanent displacement on both sides of the fault, and results in severe damage of the structure located near or crossing fault. Therefore, the seismic action on the bridge across fault is different on both sides of the fault, including the ground motion with permanent displacement. A simply supported girder bridge vertically across strike-slip fault is taken as an example, and a three-dimensional fnite-element model is developed using the earthquake engineering simulation software framework OpenSees (Open System for Earthquake Engineering Simulation). The seismic response of the bridge structure under different modes of seismic action is analyzed. The modes of seismic action considered include: (1) Seismic action of sites on both sides of the fault is regarded as the same ground motion without permanent displacement, that is, a consistent ground motion mode of seismic action without permanent displacement; (2) Seismic action of sites is regarded as a ground motion with permanent displacement on the active side of the fault, and a ground motion without permanent displacement on the passive side of the fault, that is, a non-consistent ground motion mode of seismic action with permanent displacement; (3) Fault dislocation displacement is applied to the site on the active side of fault, while the site on the passive side is fixed, that is, a static force mode of seismic action with fault dislocation displacement. The analysis results show that a consistent ground motion mode of seismic action without permanent displacement leads to a significant underestimation of the bridge structure responses, which indicates that the difference of permanent displacements of ground motions on both sides of the fault significantly increases the bridge structure response; The result from the superposition of the consistent ground motion mode and static force mode of fault dislocation permanent displacement is very close to the result from the non-consistent ground motion mode. Therefore, to a certain extent, the seismic response of cross-fault bridge structures can be approximately simulated by combining a consistent ground motion mode and a static force mode of fault dislocation displacement.
-
引言
隧道在穿越地震活动断层带时可能会遭受严重破坏,然而隧道选线常不可避免地穿越地震活动断裂带,如国外穿越东非大裂谷的肯尼亚内罗毕—马拉巴标准轨铁路(内马铁路)隧道和我国穿越龙门山断裂带的广元—甘肃高速公路隧道及穿越了多条地震活动断裂带的高黎贡山铁路隧道(何川等,2014;Tsinidis等,2020;文鑫涛等,2021)。解决好隧道穿越地震活动断裂带问题对于我国川藏铁路建设、推进“西部大开发”政策及保障国家和人民生命财产安全具有重要意义。
关于地震和断层错动对隧道结构的影响,相关学者进行了广泛研究(耿萍等,2012;刘学增等,2013;Baziar等,2014;张景等,2017;Yan等,2020)。已有研究表明,相较于受地震动或断层错动单一因素作用,修筑于地震活动断裂带的隧道在受到地震动和断层错动联合作用时,隧道整体结构将发生更剧烈的破坏(Fan等,2020)。Fan等(2020)通过振动台试验研究了三维正断层错动与地震动联合作用下跨断层隧道的地震响应,并得出了断层滑移导致地震动作用下隧道整体刚度下降的结论。Anastasopoulos等(2008)研究了正断层破裂变形与地震动联合作用下深埋沉管隧道特性,认为在正断层错动与强地震动联合作用下隧道结构安全将受到极大威胁。Shen等(2020)以汶川地震中龙溪隧道为背景,通过振动台试验分析了穿越断层隧道震害,提出隧道结构会经历断层运动和地震运动阶段的地震破坏,且断层错动对隧道结构的破坏较地震作用严重。Zhen等(2022)采用数值方法研究了断层错动和后续地震激励对隧道的影响,考虑隧道岩石-隧道界面的弱化,分析了断层破裂变形与地震动联合作用下隧道的变形响应和破坏机理。闫高明等(2019)通过试验研究了地震动和断层错动联合作用下柔性接头设计对隧道结构减震抗断错的效果。然而目前针对地震动和断层错动联合作用的研究较少,且普遍通过模型试验和数值模拟分析地震动与断层错动联合作用下隧道的变形响应和破坏机理,未提出可用于工程设计的理论分析方法。
一般采用拟静力法进行地下结构抗震理论分析,主要包括自由场变形法和土体-结构相互作用法(St John等,1987;Hashash等,2001),将地震动作用等效为静荷载,用静力弹性计算模型分析隧道结构变形及内力变化,通过静态荷载模拟地下结构的动力响应。王明年等(2011)基于弹性地基梁理论,提出了地震波沿隧道减震结构轴线方向作用的理论解。Yu等(2018)考虑了地层变化及刚度突变对隧道结构纵向响应的影响,提出了穿越土岩变化地层的隧道结构纵向地震响应解析方法。刘国钊等(2020)基于不同的分析方法提出了地下管道在断层错动作用下的纵向响应。
受上述工作的启发,本文基于拟静力方法,建立地震动和断层错动联合作用的隧道纵向响应理论模型。基于Winkler弹性地基梁的计算假定,沿隧道纵向施加静力荷载,以模拟地震动作用,在断层一侧施加相应断层错动位移,以实现断层的位错作用,通过隧道结构静力平衡及连续边界条件,求解沿隧道纵向挠曲线方程。借助MATLAB软件编制相应的计算程序,分析在地震动和断层错动联合作用下隧道结构纵向变形及内力分布特征,讨论断层错动位移、两侧围岩地基系数与断层破碎带地基系数比及围岩与隧道结构刚度比等参数对隧道结构变形及内力的影响,为穿越地震活动断裂带的隧道结构设计提供参考。
1. 理论模型与解析解
1.1 理论模型
弹性地基梁理论为岩土工程中常用的计算方法,其核心思想是将地基离散为刚性支座上一系列独立弹簧,且地基表面任意点的位移与该点单位面积上所受的压力成正比(Liang等,2020)。本文基于弹性地基梁理论,考虑围岩和结构之间的相互作用,将地震动作用下的自由场位移yu简化为理想的正弦波形(St John等,1987),并沿隧道纵向以静荷载形式施加该地震动,同时在理论模型上盘施加相应错动位移(图1),通过静力弹性地基梁理论进行求解,以研究地震动和断层错动联合作用下隧道的纵向响应。
计算假定如下:
(1)隧道结构简化为二维平面内实截面梁,且隧道结构在断层错动作用下是xoy平面内的变形;
(2)隧道结构与围岩的相互作用通过地基弹簧的形式体现;
(3)忽略自重应力及构造应力等初始应力;
(4)忽略断层错动的时间动力效应。
基于上述假设,隧道在断层错动作用下,假设隧道与围岩始终紧密贴合,而在地基弹簧作用下,隧道结构内部会产生相应的剪力、弯矩等内力。
由地震动作用引起的土体自由场位移yu可简化为理想的正弦波形,地震波如图2所示,自由场位移函数可表示为:
$$ {y_{\rm{u}}} = {y_{\max }}\sin \left( {\frac{{2{\text{π}} }}{L}x + {\gamma _0}} \right) $$ (1) 式中,ymax为地震自由场峰值位移;L为地震波波长;
$ {\gamma _0} $ 为相位角,可通过改变相位角实现地震剪切波在隧道结构中的纵向移动;x为距原点的距离。根据弹性地基梁理论,隧道结构弯曲控制方程如下:
$$ EI\frac{{{{\rm{d}}^4}{y_{\rm{t}}}}}{{{\rm{d}}{x^4}}} + P = 0 $$ (2) 式中,EI为隧道结构抗弯刚度;P为隧道结构与围岩之间的相互作用力。
地基与结构紧密接触,地基对结构的作用可简化为连续独立弹簧的作用:
$$ P = K\left( {{y_{\rm{t}}} - {y_{\rm{u}}} \pm {y_{\rm{f}}}} \right) $$ (3) 式中,yt为隧道结构绝对位移;yf为断层错动位移,对于倾滑正断层,断层上盘相对下降,取正值,对于倾滑逆断层,断层上盘相对上升,取负值;K为地基反力系数,其为关于地震波波长L的函数(Vesić,1961),可通过式(4)确定。
$$ K = \frac{{8{\text{π}} {E_{\text{s}}}(1 - {v_{\text{s}}})}}{{(3 - 4{v_{\text{s}}})(1 + {v_{\text{s}}})}}\frac{b}{L} $$ (4) 式中,b为梁的有效计算宽度,在此为隧道直径;Es为围岩弹性模量;vs为围岩泊松比。
将式(1)和式(3)代入式(2)可得到隧道基本控制方程为:
$$ EI\frac{{{{\rm{d}}^4}{y_{\rm{t}}}}}{{{\rm{d}}{x^4}}} + K{y_{\rm{t}}} = K{y_{\max }}\sin \left( {\frac{{2{\text{π}} }}{L}x + {\gamma _0}} \right) \pm K{y_{\rm{f}}} $$ (5) 1.2 解析解
根据式(5),得到隧道全段控制方程为:
$$ \left\{ \begin{array}{l} EI\dfrac{{{{\rm{d}}^4}{y_{\rm{t}}}}}{{{\rm{d}}{x^4}}} + K{y_{\rm{t}}} = K{y_{\max }}\sin \left( {wx + {\gamma_0}} \right),\quad{x < - {l_{\rm{f}}}} \\ EI\dfrac{{{{\rm{d}}^4}{y_{\rm{t}}}}}{{{\rm{d}}{x^4}}} + {K_{\rm{f}}}{y_{\rm{t}}} = {K_{\rm{f}}}{y_{\max {\rm{f}}}}\sin \left( {wx + {\gamma_0}} \right),\quad{ - {l_{\rm{f}}} \leqslant x < 0} \\ EI\dfrac{{{{\rm{d}}^4}{y_{\rm{t}}}}}{{{\rm{d}}{x^4}}} + {K_{\rm{f}}}{y_{\rm{t}}} = {K_{\rm{f}}}{y_{\max {\rm{f}}}}\sin \left( {wx + {\gamma_0}} \right) \pm {K_{\rm{f}}}{y_{\rm{f}}},\quad {0 \leqslant x < {l_{\rm{f}}}} \\ EI\dfrac{{{{\rm{d}}^4}{y_{\rm{t}}}}}{{{\rm{d}}{x^4}}} + K{y_{\rm{t}}} = K{y_{\max }}\sin \left( {wx + {\gamma_0}} \right) \pm K{y_{\rm{f}}},\quad { x\geqslant {l_{\rm{f}}}} \\ \end{array} \right. $$ (6) 式中,Kf为断层破碎带地基系数;ymaxf为断层破碎带中地震自由场峰值位移。
若x< lf,隧道结构处于断层下盘,求解控制方程得到隧道结构挠曲线方程为:
$$ {y_{\rm{t}}} = {{\rm{e}}^{\alpha x}}\left( {{A_1}\cos \alpha x + {B_1}\sin \alpha x} \right) + {{\rm{e}}^{ - \alpha x}}\left( {{C_1}\cos \alpha x + {D_1}\sin \alpha x} \right) + \frac{{a{y_{\max }}\sin \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}} $$ (7) 式中,
$\alpha = \sqrt [\uproot{4}{{4}}]{{K/4 EI}}$ ;A1、B1、C1、D1为待定积分常数;$ a = K/EI $ ;$w = 2{\text{π}} /L$ 。将隧道结构假设为岩土体中两端自由的半无限长梁,当
$ x \to - \infty $ 时,应有:$$ {y_{\rm{t}}} = \frac{{a{u_{\max }}\sin \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}} $$ (8) 结合式(7)和式(8)可知,C1=0、D1 =0,式(7)可简化为:
$$ {y_{\rm{t}}} = {{\rm{e}}^{\alpha x}}\left( {{A_1}\cos \alpha x + {B_1}\sin \alpha x} \right) + \frac{{a{y_{\max }}\sin \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}} $$ (9) 同样的,若
$ x \geqslant {l_{\text{f}}} $ ,当$ x \to + \infty $ 时,有A4=0、B4=0,则隧道结构全段挠曲线方程为:$$ \left\{ \begin{gathered} {y_{\rm{t}}} = {{\rm{e}}^{\alpha x}}\left( {{A_1}\cos \alpha x + {B_1}\sin \alpha x} \right) + \frac{{a{y_{\max }}\sin \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}},\quad {x < - {l_{\rm{f}}}} \\ {y_{\rm{t}}} = {{\rm{e}}^{{\alpha _{\rm{f}}}x}}\left( {{A_2}\cos {\alpha _{\rm{f}}}x + {B_2}\sin {\alpha _{\rm{f}}}x} \right) + {{\rm{e}}^{ - {\alpha _{\rm{f}}}x}}\left( {{C_2}\cos {\alpha _{\rm{f}}}x + {D_2}\sin {\alpha _{\rm{f}}}x} \right) + \frac{{{a_{\rm{f}}}{y_{\max {\rm{f}}}}\sin \left( {wx + {\gamma _0}} \right)}}{{{a_{\rm{f}}} + {w^4}}},\quad { - {l_{\rm{f}}} \leqslant x < 0} \\ {y_{\rm{t}}} = {{\rm{e}}^{{\alpha _{\rm{f}}}x}}\left( {{A_3}\cos {\alpha _{\rm{f}}}x + {B_3}\sin {\alpha _{\rm{f}}}x} \right) + {{\rm{e}}^{ - {\alpha _{\rm{f}}}x}}\left( {{C_3}\cos {\alpha _{\rm{f}}}x + {D_3}\sin {\alpha _{\rm{f}}}x} \right) + \frac{{{a_{\rm{f}}}{y_{\max {\rm{f}}}}\sin \left( {wx + {\gamma _0}} \right)}}{{{a_{\rm{f}}} + {w^4}}} \pm {y_{\rm{f}}},\quad {0 \leqslant x < {l_{\rm{f}}}} \\ {y_{\rm{t}}} = {{\rm{e}}^{ - \alpha x}}\left( {{C_4}\cos \alpha x + {D_4}\sin \alpha x} \right) + \frac{{a{y_{\max }}\sin \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}} \pm {y_{\rm{f}}},\quad { x\leqslant {l_{\rm{f}}}} \\ \end{gathered} \right. $$ (10) 式中,
${\alpha _{\rm{f}}} = \sqrt[{\uproot{4}{{4}}}]{{{K_{\rm{f}}}/4 EI}}$ ;af =Kf /EI;Ai、Bi、Ci、Di (i=1,2,3,4)为待定系数。对于隧道结构任意截面的弯矩、剪力,可由隧道结构挠度二阶导数和三阶导数乘以−EI求得,则隧道结构全段弯矩
$ {M_{\rm{t}}} $ 、剪力$ {Q_{\rm{t}}} $ 可(12)表示为:$$ \left\{ \begin{array}{l} {M_{\rm{t}}} = - EI\left( {2{\alpha ^2}{{\rm{e}}^{\alpha x}}\left( { - {A_1}\sin \alpha x + {B_1}\cos \alpha x} \right) - {w^2}\dfrac{{a{y_{\max }}\sin \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}}} \right),\quad {x < - {l_{\rm{f}}}} \\ {{M_{\rm{t}}} = - EI\left( \begin{gathered} 2\alpha _{\rm{f}}^2{{\rm{e}}^{{\alpha _{\rm{f}}}x}}\left( {{B_2}\cos {\alpha _{\rm{f}}}x - {A_2}\sin {\alpha _{\rm{f}}}x} \right) + 2\alpha _{\rm{f}}^2{{\rm{e}}^{ - {\alpha _{\rm{f}}}x}}\left( { - {D_2}\cos {\alpha _{\rm{f}}}x + {C_2}\sin {\alpha _{\rm{f}}}x} \right) \\ - {w^2}\frac{{{a_{\rm{f}}}{y_{\max {\rm{f}}}}\sin \left( {wx + {\gamma _0}} \right)}}{{{a_{\rm{f}}} + {w^4}}} \\ \end{gathered} \right)},\quad { - {l_{\rm{f}}} \leqslant x < 0} \\ {{M_{\rm{t}}} = - EI\left( \begin{gathered} 2\alpha _{\rm{f}}^2{{\rm{e}}^{{\alpha _{\rm{f}}}x}}\left( {{B_3}\cos {\alpha _{\rm{f}}}x - {A_3}\sin {\alpha _{\rm{f}}}x} \right) + 2\alpha _{\rm{f}}^2{{\rm{e}}^{ - {\alpha _{\rm{f}}}x}}\left( { - {D_3}\cos {\alpha _{\rm{f}}}x + {C_3}\sin {\alpha _{\rm{f}}}x} \right) \\ - {w^2}\frac{{{a_{\rm{f}}}{y_{\max {\rm{f}}}}\sin \left( {wx + {\gamma _0}} \right)}}{{{a_{\rm{f}}} + {w^4}}} \\ \end{gathered} \right)},\quad {0 \leqslant x < {l_{\rm{f}}}} \\ {M_{\rm{t}}} = - EI\left( {2{\alpha ^2}{{\rm{e}}^{ - \alpha x}}\left( {{C_4}\sin \alpha x - {D_4}\cos \alpha x} \right) - \dfrac{{{w^2}a{y_{\max }}\sin \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}}} \right),\quad { x\geqslant {l_{\rm{f}}} } \\ \end{array} \right. $$ (11) $$ \left\{ \begin{array}{l} {Q_{\rm{t}}} = - EI\left( {2{\alpha ^3}{{\rm{e}}^{\alpha x}}\left( { - {A_1}(\cos \alpha x + \sin \alpha x) + {B_1}(\cos \alpha x - \sin \alpha x)} \right) - \dfrac{{{w^3}a{y_{\max }}\cos \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}}} \right),\quad {x < - {l_{\rm{f}}}} \\ {{Q_{\rm{t}}} = - EI\left( \begin{gathered} 2\alpha _{\rm{f}}^3{{\rm{e}}^{{\alpha _{\rm{f}}}x}}\left( { - {A_2}(\cos {\alpha _{\rm{f}}}x + \sin {\alpha _{\rm{f}}}) + {B_2}(\cos {\alpha _{\rm{f}}}x - \sin {\alpha _{\rm{f}}})} \right) \\ + 2{\alpha _{\rm{f}}}^3{{\rm{e}}^{ - {\alpha _{\rm{f}}}x}}\left( {{C_2}\left( {\cos {\alpha _{\rm{f}}}x - \sin {\alpha _{\rm{f}}}x} \right) + {D_2}\left( {\cos {\alpha _{\rm{f}}}x + \sin {\alpha _{\rm{f}}}x} \right)} \right) - \frac{{{w^3}{a_{\rm{f}}}{y_{\max {\rm{f}}}}\cos \left( {wx + {\gamma _0}} \right)}}{{{a_{\rm{f}}} + {w^4}}} \\ \end{gathered} \right)},\quad { - {l_{\rm{f}}} \leqslant x < 0} \\ {{Q_{\rm{t}}} = - EI\left( \begin{gathered} 2\alpha _{\rm{f}}^3{{\rm{e}}^{{\alpha _{\rm{f}}}x}}\left( { - {A_3}(\cos {\alpha _{\rm{f}}}x + \sin {\alpha _{\rm{f}}}) + {B_3}(\cos {\alpha _{\rm{f}}}x - \sin {\alpha _{\rm{f}}})} \right) \\ + 2\alpha _{\rm{f}}^33{{\rm{e}}^{ - {\alpha _{\rm{f}}}x}}\left( {{C_3}\left( {\cos {\alpha _{\rm{f}}}x - \sin {\alpha _{\rm{f}}}x} \right) + {D_3}\left( {\cos {\alpha _{\rm{f}}}x + \sin {\alpha _{\rm{f}}}x} \right)} \right) - \frac{{{w^3}{a_{\rm{f}}}{y_{\max {\rm{f}}}}\cos \left( {wx + {\gamma _0}} \right)}}{{{a_{\rm{f}}} + {w^4}}} \\ \end{gathered} \right)},\quad {0 \leqslant x < {l_{\rm{f}}}} \\ {Q_{\rm{t}}} = - EI\left( {2{\alpha ^3}{{\rm{e}}^{ - \alpha x}}\left( {{C_4}(\cos \alpha x - \sin \alpha x) + {D_4}(\cos \alpha x + \sin \alpha x)} \right) - \dfrac{{{w^3}a{y_{\max }}\cos \left( {wx + {\gamma _0}} \right)}}{{a + {w^4}}}} \right) ,\quad {x \geqslant {l_{\rm{f}}}} \\ \end{array} \right. $$ (12) 同时在
$ x = - {l_{\rm{f}}} $ 、$ x = 0 $ 、$ x = {l_{\rm{f}}} $ 处隧道结构挠度、转角、弯矩、剪力应满足式(13)至式(15)的连续性条件:$$ \underset{x\to -{l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}=\underset{x\to -{l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}},\;\; \underset{x\to -{l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime }=\underset{x\to -{l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime },\;\; EI\underset{x\to -{l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime }=EI\underset{x\to -{l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime },\;\; EI\underset{x\to -{l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime \text{}\prime }=EI\underset{x\to -{l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime \text{}\prime } $$ (13) $$ \underset{x\to {0}^-}{\mathrm{lim}}{y}_{\rm{t}}=\underset{x\to {0}^+}{\mathrm{lim}}{y}_{\rm{t}},\;\; \underset{x\to {0}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime }=\underset{x\to {0}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime },\;\; EI\underset{x\to {0}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime }=EI\underset{x\to {0}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime },\;\; EI\underset{x\to {0}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime \text{}\prime }=EI\underset{x\to {0}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime \text{}\prime }$$ (14) $$ \underset{x\to {l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}=\underset{x\to {l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}},\;\; \underset{x\to {l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime }=\underset{x\to {l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime },\;\; EI\underset{x\to {l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime }=EI\underset{x\to {l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime },\;\; EI\underset{x\to {l}_{\rm{f}}{}^-}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime \text{}\prime }=EI\underset{x\to {l}_{\rm{f}}{}^+}{\mathrm{lim}}{y}_{\rm{t}}{}^{\prime \text{}\prime \text{}\prime } $$ (15) 式中,
$ {y_{\rm{t}}}^\prime $ 、$ {y_{\rm{t}}}^{\prime \prime } $ 、$ {y_{\rm{t}}}^{\prime \prime \prime } $ 分别为隧道结构挠度的对x一阶、二阶、三阶导数。2. 模型验证
内马铁路一期项目三标段隧道位于东非大裂谷东翼区(甘星球等,2021),受东非大裂谷活跃的地质作用影响显著,线路经过的活动断层地段岩体较破碎、节理裂隙发育,如图3所示。本文以该隧道工程为背景,根据前文建立的理论模型,采用MATLAB软件编程计算解析解,并利用ABAQUS有限元分析软件建模计算数值解进行对比,验证本文建立的理论模型及解析解的有效性。
使用Midas SoilWorks软件进行一维自由场分析,分别选择正常围岩剖面和断层破碎带剖面进行计算(图3中Fb9-2断层),以弱风化粗面岩作为一维场地分析的基岩,地震波选取峰值地面加速度为0.1 g的El Centro波(图4),波长L=280 m,地震波由基岩处输入。通过计算得出地震自由场峰值位移ymax=0.027 m,断层破碎带中地震自由场峰值位移ymaxf =0.052 m。利用ABAQUS有限元分析软件建立数值模型,隧道结构采用线性两结点梁单元模型进行模拟,选择欧拉-伯努利梁,仅考虑梁的弯曲变形,忽略其剪切变形,纵向计算长度为560 m。将有限元模型梁单元截面简化为矩形,梁截面惯性矩与隧道结构相等,即I=173.63 m4 ,梁截面等效宽度b=6.2 m,基于惯性矩和等效宽度得出梁单元截面高度h=6.96 m,将有限元梁模型划分为5 600个单元,在ABAQUS软件Interaction模块中选择接地弹簧,并在每个单元结点上均连接上接地弹簧,通过接地弹簧模拟隧道与围岩的相互作用,将地震波等效荷载施加于梁上,断层错动位移从地基弹簧远端施加,计算得到地基系数K=1 083.4 MPa,Kf =332.7 MPa。
算例中隧道结构截面为马蹄形,隧道穿过倾角为90°的倾向逆断层(图3中Fb9-2断层),断层错动位移设为0.05 m,进行模型验证时认为隧道结构穿过均质岩体,即K=Kf =1 083.4 MPa,ymax=ymaxf=0.027 m,且忽略隧道截面形状对隧道纵向受力变形的影响,相关计算参数如表1所示。将MATLAB理论模型解析解与有限元模型数值解进行对比分析,结果如图5所示。
表 1 计算参数Table 1. Calculation parameter table材料 弹性模量/GPa 围岩泊松比 围岩密度/(kg·m−3) 地基系数/MPa 隧道截面惯性矩/m4 隧道宽度/m 强风化粗面岩 6.5 0.32 2 400 1 083.4 — — 断层破碎带 2.0 0.30 2 200 332.7 — — 隧道衬砌 35.0 0.20 — — 173.63 6.2 由图5(a)可知,在
$ {\gamma _0} = 0 $ 时刻,地震动作用在隧道结构上,导致隧道结构产生相应的弯曲变形,因为断层错动,隧道结构右半部分随断层上盘整体向上抬升,隧道结构挠度整体呈S形分布,同时隧道结构变形引起隧道内力变化,由图5(b)、图5(c)可知,隧道结构弯矩响应明显大于剪力响应,同时在断层错动界面,即x=0点处,由断层错动引起的内力变化明显大于由地震动引起的内力变化,说明在断层错动区域,隧道结构变形及内力变化主要由断层错动引起。改变相位角使地震剪切波在隧道结构中移动,在$ {\gamma _0} ={\text{π}} /2 $ 时刻,隧道变形及内力均发生变化。由图5(a)可知,隧道结构挠度在断层交界面处变化更迅速,峰值挠度略小于$ {\gamma _0} = 0 $ 时。由图5(b)可知,地震剪切波移动导致隧道结构弯矩由正弯矩变为负弯矩,由顶部受拉变为底部受拉。同时断层交界面处峰值负弯矩明显减小,由${{ - 7}} \times {\text{1}}{{\text{0}}^{\text{5}}}$ kN·m变为${{ - 6}} \times {\text{1}}{{\text{0}}^{\text{5}}}$ kN·m,这是因为在$ {\gamma _0} = {\text{π}} /2 $ 时刻,地震剪切波在断层交界面使隧道产生正弯矩,抵消了一部分断层错动引起的负弯矩。对于隧道结构剪力,由图5(c)可知,地震剪切波移动对隧道结构峰值剪力的影响较小,断层交界面两侧剪力不再对称,断层交界面左侧剪力略小于右侧剪力。总体而言,理论解析解与有限元模型数值解吻合较好,隧道结构截面挠度误差较小,对于隧道截面内力,断层错动产生的内力峰值相对误差均小于5%,满足工程精度要求,验证了本文计算方法的有效性。
3. 参数敏感性分析
断层破碎带中的岩层成因复杂,与两侧岩体在力学特性上存在明显差异,断层破碎带等岩层性质突变位置也是隧道震害最集中的位置。断层错动位移是断层错动最直观的参数(孙飞等,2019;Zhong等,2020),而围岩与隧道结构刚度比对拟静力方法的影响较大(Wang,1993)。为更直观地体现参数变化带来的影响,基于前文计算参数,采用单一变量原则,进一步分析各参数对地震动和断层错动联合作用下隧道结构纵向响应的影响规律。
3.1 断层错动位移的影响
为探究断层错动位移对隧道结构变形及内力的影响,设断层破裂带宽度为10 m(−lf<x< lf),其余参数同前文,断层错动位移取0、3、5、7、10 cm进行计算分析,结果如图6所示。由图6可知,当断层错动位移为零时,隧道仅受地震动作用,在断层破碎带中隧道挠度变化不明显,然而对于隧道结构内力,由于断层破碎带岩层性质与左右侧岩层性质相差较大,在地震动作用下岩层性质的突变使断层破碎带交界面岩层存在位移差,使隧道结构弯矩在破碎带交界面上产生负弯矩。当断层错动位移不为零时,断层破碎带内隧道截面挠度随着断层错动位移的增大而增大。隧道截面弯矩、剪力在断层错动交界面附近50 m的范围内随断层错动位移的增大而增大,且剪力反弯点与弯矩反弯点一致。另外,截面弯矩在断层错动界面左右侧呈反对称分布,最大值出现在错动界面两侧;截面剪力基本呈对称分布,剪力最大值出现在断层交界面,且由于存在断层破裂带,断层破碎带位置处剪力先骤减后突增,这是因为剪力的变化主要由断层错动作用下的围岩抗力导致,而在断层破碎带处围岩软硬突变的位置中,隧道结构产生的位移不相等,出现隧道结构被剪切的效果。断层错动位移增加使隧道弯矩和剪力接近线性增加,而不改变弯矩和剪力沿隧道纵向分布的影响范围。
3.2 两侧围岩地基系数与断层破碎带地基系数比的影响
断层破碎带是隧道施工中常见的不良地质,为探究断层破碎带围岩特性对隧道结构变形及内力的影响,引入参数λ=K/Kf表示两侧围岩地基系数与断层破碎带地基系数的比值,保持K不变,改变断层破碎带地基系数Kf,分别取参数λ、0.1λ、0.2λ、5λ、10λ进行分析,结果如图7所示。由图7可知,随着λ的减小,隧道结构弯矩变化平缓,隧道截面弯矩峰值增大,而随着Kf的减小,弯矩产生明显突变。隧道结构剪力受λ的影响最大,随着λ的增大,Kf减小,隧道剪力由1.0×105 kN增至2.5×105 kN,在断层破碎带交界面剪力骤减至−1.0×105 kN,这是因为断层破碎带围岩性质差,与两侧围岩相差过大,导致隧道结构剪力急剧变化,如图7(c)所示,这种情况下隧道结构易发生破坏,因此当隧道穿越断层破碎带时,应关注两侧围岩地基系数与断层破碎带地基系数比,有需要时对断层破碎带围岩进行补强加固。
3.3 围岩与隧道结构刚度比的影响
围岩与隧道结构刚度比对拟静力方法的影响较大,围岩与隧道结构刚度比也是工程实际中关注的重点。以围岩与隧道结构刚度比a=K/EI为基准,其他参数保持不变,通过改变隧道结构刚度EI,分别取0.01a、0.1a、a、10a、100a进行分析,结果如图8所示。
由图8(a)可知,随着隧道结构刚度EI的不断减小,围岩与隧道结构刚度比a增大,土质越硬,对隧道结构的约束作用增加,由于断层错动作用,隧道结构挠度变化更迅速,挠度变化范围减小,且在断层破碎带与两侧围岩分界面上出现挠度突变,这是因为隧道结构相对更“柔”,地震动作用在断层破碎带类软弱交界面对隧道结构的影响更大。柔性结构使隧道结果在地震动和断层错动联合作用下发生较大的变形,隧道结构内力响应减小,如图8(b)、图8(c)所示。
随着隧道结构刚度EI的增大,围岩与隧道结构刚度比a减小,隧道结构越来越“刚”,由图8(a)可知,隧道结构由地震动作用引起的位移明显减小,同时在断层破碎带内,断层引起的隧道结构挠度变化更平缓,挠度变化范围增大。然而,图8(b)、图8(c)揭示了由于隧道结构刚度增大,在地震动和断层错动联合作用下隧道结构即便产生微小形变,隧道内力(弯矩、剪力)也发生较大变化。同时,在断层破碎带两侧交界面上,由于地震动通过围岩软硬交界面,隧道结构产生了相对位移差,同样导致了隧道截面剪力急剧变化,不利于隧道结构的正常使用。因此,隧道结构穿越活动断裂带时,不能一味增大隧道刚度,应适当考虑采取柔性减震抗断错措施。
4. 结论
通过本文研究得出以下结论:
(1)建立了地震动和断层错动联合作用下隧道结构纵向响应理论模型,提出了更接近工程实际的理论分析方法,并基于实际工程背景,与数值算例进行对比,验证了该方法的有效性。
(2)断层错动位移增加使隧道结构截面弯矩、剪力峰值接近线性增加,不改变截面弯矩、剪力沿隧道纵向分布的影响范围,同时断层破碎带界面出现了截面剪力突变。
(3)两侧围岩地基系数与断层破碎带地基系数比对隧道结构截面剪力的影响较大,随着断层破碎带地基系数Kf的减小,断层破碎带围岩性质较差,导致在地震动和断层错动联合作用下隧道结构剪力在断层破碎带界面急剧减小,隧道结构易发生破坏。
(4)围岩与隧道结构刚度比a对隧道结构的影响显著。随着隧道结构刚度EI的增大,围岩与隧道结构刚度比a减小,地震动引起的隧道位移明显减小,断层错动作用引起的隧道挠度变化范围增大,挠度变化更平缓。然而隧道结构内力响应对于围岩与隧道结构刚度比的变化较敏感,隧道结构穿越活动断裂带时,不能一味增大隧道刚度,应适当考虑采取柔性减震抗断错措施。
-
表 1 原始地震动时程信息
Table 1. Original ground motion parameters
方向 PGA /g PGV/(cm·s−1) PGD/cm 横桥向 0.73 133.33 113.87 顺桥向 0.79 28.09 25.52 表 2 合成的地震动时程信息
Table 2. Synthetic ground motion parameters
类别 PGA /g PGV /(cm·s−1) PGD /cm 横桥向不考虑永久位移 0.60 111.97 105.39 横桥向考虑永久位移 0.60 117.79 138.88 顺桥向不考虑永久位移 0.60 118.89 110.84 表 3 横桥向P3和P4墩顶相对位移最大值和残余值
Table 3. Maximum and residual relative displacement of the pier top at P3 and P4 in the transverse direction
工况 最大值/cm
(残余值/cm)P2 P3 P4 P5 P6 1 9.0580
(−0.1654)4.6155
(−0.0405)3.8598
(−0.0254)6.4269
(−0.1254)5.5834
(−0.1095)2 9.8470
(0.2830)4.4210
(−1.2270)4.1927
(1.1377)6.3951
(−0.0789)5.5772
(−0.0757)3 9.2286
(0.1558)4.7349
(−1.4285)4.7175
(1.2923)6.4440
(−0.1007)5.5611
(−0.1431)表 4 横桥向P3和P4墩底剪力最大值和残余值
Table 4. The maximum and residual shear force at the bottom of pier P3 and P4 in the transverse direction
工况 剪力最大值/kN 残余剪力/kN P3 P4 P3 P4 1 7.9340$ \times {10}^{3} $ 7.6760$ \times {10}^{3} $ 0.0044$ \times {10}^{3} $ −0.0084$ \times {10}^{3} $ 2 7.6930$ \times {10}^{3} $ 8.1211$ \times {10}^{3} $ 1.7479$ \times {10}^{3} $ −1.7280$ \times {10}^{3} $ 3 8.5876$ \times {10}^{3} $ 8.9535$ \times {10}^{3} $ 1.9306$ \times {10}^{3} $ −1.9704$ \times {10}^{3} $ 表 5 横桥向P3和P4墩底弯矩最大值和残余值
Table 5. The maximum and residual bending moment at the bottom of pier P3 and P4 in the transverse direction
工况 弯矩最大值 /(kN∙m) 残余弯矩 /(kN∙m) P3 P4 P3 P4 1 1.7790$ \times {10}^{5} $ 1.6317$ \times {10}^{5} $ 0.0016$ \times {10}^{4} $ 0.0328$ \times {10}^{4} $ 2 1.6921$ \times {10}^{5} $ 1.7596$ \times {10}^{5} $ 4.7504$ \times {10}^{4} $ 4.6370$ \times {10}^{4} $ 3 1.7830$ \times {10}^{5} $ 1.9733$ \times {10}^{5} $ 5.2868$ \times {10}^{4} $ 5.2781$ \times {10}^{4} $ 表 6 P3和P4墩底扭矩最大值和残余值
Table 6. The maximum and residual torque at the bottom of pier P3 and P4 in the transverse direction
工况 扭矩最大值 /(kN∙m) 残余扭矩 /(kN∙m) P3 P4 P3 P4 1 5.4362$ \times {10}^{4} $ 0.6780$ \times {10}^{4} $ 0.0225$ \times {10}^{4} $ 0.0054$ \times {10}^{4} $ 2 5.3635$ \times {10}^{4} $ 2.2612$ \times {10}^{4} $ 2.0436$ \times {10}^{4} $ 2.0624$ \times {10}^{4} $ 3 5.2122$ \times {10}^{4} $ 1.8589$ \times {10}^{4} $ 1.7523$ \times {10}^{4} $ 1.7337$ \times {10}^{4} $ -
陈令坤, 张楠, 夏禾, 2016. 方向脉冲及竖向效应对高铁桥梁地震响应影响. 振动工程学报, 29(4): 704—713Chen L. K. , Zhang N. , Xia H. , 2016. Effects of directivity pulse and vertical earthquake on seismic response of high-speed railway bridge. Journal of Vibration Engineering, 29(4): 704—713. (in Chinese) 杜修力, 韩强, 李忠献等, 2008.5. 12汶川地震中山区公路桥梁震害及启示. 北京工业大学学报, 34(12): 1270—1279Du X. L. , Han Q. , Li Z. X. , et al. , 2008. The seismic damage of bridges in the 2008 Wenchuan earthquake and lessons from its damage. Journal of Beijing University of Technology, 34(12): 1270—1279. (in Chinese) 管仲国, 游瀚, 郭河, 2016. 近断层斜拉桥弹塑性索与阻尼器组合横向减震. 同济大学学报(自然科学版), 44(11): 1653—1659Guan Z. G. , You H. , Guo H. , 2016. Responses of cable-stayed bridge transversely isolated with elasto-plastic cable pairs and fluid viscous damper subjected to near fault ground motions. Journal of Tongji University (Natural Science), 44(11): 1653—1659. (in Chinese) 惠迎新, 王克海, 李冲, 2014. 跨断层地表破裂带桥梁震害研究及抗震概念设计. 公路交通科技, 31(10): 51—57 doi: 10.3969/j.issn.1002-0268.2014.10.009Hui Y. X. , Wang K. H. , Li C. , 2014. Study of seismic damage and seismic conceptual design of bridges across fault surface rupture zones. Journal of Highway and Transportation Research and Development, 31(10): 51—57. (in Chinese) doi: 10.3969/j.issn.1002-0268.2014.10.009 惠迎新, 王克海, 2015 a. 基于多点激励位移输入模型的跨断层桥梁地震动输入方法. 东南大学学报(自然科学版), 45(3): 557—562Hui Y. X. , Wang K. H. , 2015 a. Earthquake motion input method for bridges crossing fault based on multi-support excitation displacement input model. Journal of Southeast University (Natural Science Edition), 45(3): 557—562. (in Chinese) 惠迎新, 王克海, 2015 b. 跨断层桥梁地震响应特性研究. 桥梁建设, 45(3): 70—75Hui Y. X. , Wang K. H. , 2015 b. Study of seismic response features of bridges crossing faults. Bridge Construction, 45(3): 70—75. (in Chinese) 惠迎新, 王克海, 吴刚等, 2015 c. 跨断层桥梁地震响应分析及合理跨越角度研究. 振动与冲击, 34(13): 6—11, 17Hui Y. X. , Wang K. H. , Wu G. , et al. , 2015 c. Seismic responses of bridges crossing faults and their best crossing angles. Journal of Vibration and Shock, 34(13): 6—11, 17. (in Chinese) 惠迎新, 台玉吉, 王克海等, 2017. 跨断层桥梁抗震若干问题探讨. 地震工程学报, 39(5): 870—875, 889Hui Y. X. , Tai Y. J. , Wang K. H. , et al. , 2017. Discussion about earthquake resistance of the bridges crossing the active fault. China Earthquake Engineering Journal, 39(5): 870—875, 889. (in Chinese) 江辉, 朱晞, 倪永军, 2011. 基于谱分析与RC梁桥检验的脉冲型近场地震动强度度量指标. 铁道学报, 33(3): 91—99 doi: 10.3969/j.issn.1001-8360.2011.03.016Jiang H. , Zhu X. , Ni Y. J. , 2011. Intensity indice of impulsive near-fault earthquake ground motion based on spectra analysis and RC bridge check. Journal of the China Railway Society, 33(3): 91—99. (in Chinese) doi: 10.3969/j.issn.1001-8360.2011.03.016 李鸿晶, 陆鸣, 温增平等, 2009. 汶川地震桥梁震害的特征. 南京工业大学学报(自然科学版), 31(1): 24—29Li H. J. , Lu M. , Wen Z. P. , et al. , 2009. Characteristics of bridge damages in Wenchuan earthquake. Journal of Nanjing University of Technology (Natural Science Edition), 31(1): 24—29. (in Chinese) 李建中, 管仲国, 2011. 基于性能桥梁抗震设计理论发展. 工程力学, 28(S2): 24—30, 53Li J. Z. , Guan Z. G. , 2011. Performance-based seismic design for bridges. Engineering Mechanics, 28(S2): 24—30, 53. (in Chinese) 李帅, 张凡, 颜晓伟等, 2017. 近断层地震动合成方法及其对超大跨斜拉桥地震响应影响. 中国公路学报, 30(2): 86—97, 106Li S. , Zhang F. , Yan X. W. , et al. , 2017. Synthetic method for near-fault ground motions and its influence on seismic response of super-span cable-stayed bridge. China Journal of Highway and Transport, 30(2): 86—97, 106. (in Chinese) 李小军, 1992. 场地土层对地震地面运动影响的分析方法. 世界地震工程, 8(2): 49—60. 李小军, 彭青, 2001. 不同类别场地地震动参数的计算分析. 地震工程与工程振动, 21(1): 29—36Li X. J. , Peng Q. , 2001. Calculation and analysis of earthquake ground motion parameters for different site categories. Earthquake Engineering and Engineering Vibration, 21(1): 29—36. (in Chinese) 李小军, 贺秋梅, 亓兴军, 2012. 地震动速度脉冲对大跨斜拉桥减震控制的影响. 应用基础与工程科学学报, 20(2): 272—285Li X. J. , He Q. M. , Qi X. J. , 2012. Seismic mitigation control effects of long-span cable-stayed bridges to ground motions with velocity pulse. Journal of Basic Science and Engineering, 20(2): 272—285. (in Chinese) 刘爱文, 夏珊, 徐超, 2008. 汶川地震交通系统震害及震后抢修. 震灾防御技术, 3(3): 243—450.Liu A. W., Xia S., Xu C., 2008. Damage and emergency recovery of the transportation systems after Wenchuan earthquake. Technology for Earthquake Disaster Prevention, 3(3): 243–250. (in Chinese) 刘启方, 袁一凡, 金星等, 2006. 近断层地震动的基本特征. 地震工程与工程振动, 26(1): 1—10Liu Q. F. , Yuan Y. F. , Jin X. , et al. , 2006. Basic characteristics of near-fault ground motion. Earthquake Engineering and Engineering Vibration, 26(1): 1–10. (in Chinese) 石岩, 王东升, 孙治国, 2014. 近断层地震动下减隔震桥梁地震反应分析. 桥梁建设, 44(3): 19—24Shi Y. , Wang D. S. , Sun Z. G. , 2014. Analysis of seismic response of seismically mitigated and isolated bridge subjected to near-fault ground motion. Bridge Construction, 44(3): 19—24. (in Chinese) 王东升, 杨海红, 王国新, 2005. 考虑邻梁碰撞的多跨长简支梁桥落梁震害分析. 中国公路学报, 18(3): 54—59Wang D. S. , Yang H. H. , Wang G. X. , 2005. Seismic analysis of girders falling down in multi-span long simply supported bridges with adjacent pounding effects. China Journal of Highway and Transport, 18(3): 54—59. (in Chinese) 王东升, 郭迅, 孙治国等, 2009. 汶川大地震公路桥梁震害初步调查. 地震工程与工程振动, 29(3): 84—94Wang D. S. , Guo X. , Sun Z. G. , et al. , 2009. Damage to highway bridges during Wenchuan earthquake. Journal of Earthquake Engineering and Engineering Vibration, 29(3): 84—94. (in Chinese) 王东升, 孙治国, 李晓莉等, 2010. 汶川大地震曲线梁桥震害及破坏机理分析. 防灾减灾工程学报, 30(5): 572—579Wang D. S. , Sun Z. G. , Li X. L. , et al. , 2010. Seismic damage investigation and failure mechanism analysis of curved bridges in Wenchuan earthquake. Journal of Disaster Prevention and Mitigation Engineering, 30(5): 572—579. ( in Chinese 杨怀宇, 李建中, 2015. 断层地震动对隔震桥梁地震响应的影响. 同济大学学报(自然科学版), 43(8): 1144—1152Yang H. Y. , Li J. Z. , 2015. Response analysis of seismic isolated bridge under influence of fault-crossing groundmotions. Journal of Tongji University (Natural Science), 43(8): 1144—1152. (in Chinese) 袁一凡, 田启文, 2012. 工程地震学. 北京: 地震出版社. 中华人民共和国交通运输部, 2014. JTG B02—2013 公路工程抗震规范. 北京: 人民交通出版社.Ministry of Transport of the People's Republic of China, 2014. JTG B02—2013 Specification of seismic design for highway engineering. Beijing: China Communications Press. (in Chinese) 中华人民共和国交通运输部, 2020. JTG/T 2231-01—2020 公路桥梁抗震设计规范. 北京: 人民交通出版社.Ministry of Transport of the People's Republic of China, 2020. JTG/T2231-01—2020 Specifications for seismic design of highway bridges. Beijing: China Communications Press. (in Chinese) 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 2010. GB 50011—2010 建筑抗震设计规范. 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People’s Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, 2010. GB 50011—2010 Code for seismic design of buildings. Beijing: China Architecture & Building Press. (in Chinese) 周正华, 张艳梅, 孙平善等, 2003. 断层场地震害研究综述. 地震工程与工程振动, 23(5): 38—41 doi: 10.3969/j.issn.1000-1301.2003.05.006Zhou Z. H. , Zhang Y. M. , Sun P. S. , et al. , 2003. A summary of damage study on fault site. Earthquake Engineering and Engineering Vibration, 23(5): 38—41. (in Chinese) doi: 10.3969/j.issn.1000-1301.2003.05.006 Byers W. G. , Edwards C. , Tang A. , et al. , 2000. Performance of transportation systems after the 1999 Kocaeli earthquake. Earthquake Spectra, 16(S1): 403—410. Caltrans, 2012. Bridge design practice . Sacramento: California Department of Transportation. Caltrans, 2019. Seismic design criteria: version 2.0. Sacramento: California Department of Transportation. CEN, 2004. EN 1998-5—2004 Eurocode 8: design of structures for earthquake resistance Part 5: foundation, retaining structures and geotechnical aspects. Brussels: European Committee for Standardization, 14—20. Chang K. C. , Chang D. W. , Tsai M. H. , et al. , 2000. Seismic performance of highway bridges. Earthquake Engineering and Engineering Seismology, 2(1): 55—77. Elgamal A., Lu J., Mackie K., 2014. MSBrindge: OpenSees pushover and earthquake analysis of multi-span bridges - User Manual. San Diego: Department of Structural Engineering, University of California. Ghasemi H., Cooper J. D., Imbsen R., et al., 2000. The November 1999 Duzce earthquake: post-earthquake investigation of the structures on the TEM. Washington: Federal Highway Administration. Goel R. , Qu B. , Tures J. , et al. , 2014. Validation of fault rupture-response spectrum analysis method for curved bridges crossing strike-slip fault rupture zones. Journal of Bridge Engineering, 19(5): 06014002. doi: 10.1061/(ASCE)BE.1943-5592.0000602 Goel R. K., Chopra A. K., 2008. Analysis of ordinary bridges crossing fault-rupture zones. In: The 14 th World Conference on Earthquake Engineering. Beijing. Goel R. K. , Chopra A. K. , 2009 a. Linear analysis of ordinary bridges crossing fault-rupture zones. Journal of Bridge Engineering, 14(3): 203—215. doi: 10.1061/(ASCE)1084-0702(2009)14:3(203) Goel R. K. , Chopra A. K. , 2009 b. Nonlinear analysis of ordinary bridges crossing fault-rupture zones. Journal of Bridge Engineering, 14(3): 216—224. doi: 10.1061/(ASCE)1084-0702(2009)14:3(216) Güney D. , Acar M. , Özlüdemir M. T. , et al. , 2010. Investigation of post-earthquake displacements in viaducts using Geodetic and Finite Element Methods. Natural Hazards and Earth System Sciences, 10(12): 2579—2587. doi: 10.5194/nhess-10-2579-2010 He W. L. , Agrawal A. K. , 2008. Analytical model of ground motion pulses for the design and assessment of seismic protective systems. Journal of Structural Engineering, 134(7): 1177—1188. doi: 10.1061/(ASCE)0733-9445(2008)134:7(1177) Hsu Y. T. , Fu C. C. , 2004. Seismic effect on highway bridges in Chi Chi earthquake. Journal of Performance of Constructed Facilities, 18(1): 47—53. doi: 10.1061/(ASCE)0887-3828(2004)18:1(47) Imbsen R. A. , Roblee C. J. , Yashinsky M. , et al. , 2000. Impact on highway structures. Earthquake Spectra, 16(S1): 411—435. Kitagawa Y. , Hiraishi H. , 2004. Overview of the 1995 Hyogo-Ken Nanbu earthquake and proposals for earthquake mitigation measures. Journal of Japan Association for Earthquake Engineering, 4 (3): 1—29. doi: 10.5610/jaee.4.3_1 Li S. , Zhang F. , Wang J. Q. , et al. , 2017. Effects of near-fault motions and artificial pulse-type ground motions on super-span cable-stayed bridge systems. Journal of Bridge Engineering, 22(3): 04016128. doi: 10.1061/(ASCE)BE.1943-5592.0001008 Mackie K. R. , Lu J. C. , Elgamal A. , 2012. Performance-based earthquake assessment of bridge systems including ground-foundation interaction. Soil Dynamics and Earthquake Engineering, 42: 184—196. doi: 10.1016/j.soildyn.2012.05.023 Makris N. , Chang S. P. , 2000. Effect of viscous, viscoplastic and friction damping on the response of seismic isolated structures. Earthquake Engineering & Structural Dynamics, 29(1): 85—107. Mavroeidis G. P. , Papageorgiou A. S. , 2003. A mathematical representation of near-fault ground motions. Bulletin of the Seismological Society of America, 93(3): 1099—1131. doi: 10.1785/0120020100 Park S. W. , Ghasemi H. , Shen J. , et al. , 2004. Simulation of the seismic performance of the Bolu Viaduct subjected to near‐fault ground motions. Earthquake Engineering & Structural Dynamics, 33(13): 1249—1270. Roussis P. C. , Constantinou M. C. , Erdik M. , et al. , 2003. Assessment of performance of seismic isolation system of bolu viaduct. Journal of Bridge Engineering, 8(4): 182—190. doi: 10.1061/(ASCE)1084-0702(2003)8:4(182) Sehhati R. , Rodriguez-Marek A. , ElGawady M. , et al. , 2011. Effects of near-fault ground motions and equivalent pulses on multi-story structures. Engineering Structures, 33(3): 767—779. doi: 10.1016/j.engstruct.2010.11.032 Sengupta A. , Quadery L. , Sarkar S. , et al. , 2016. Influence of bidirectional near-fault excitations on RC bridge piers. Journal of Bridge Engineering, 21(7): 04016034. doi: 10.1061/(ASCE)BE.1943-5592.0000836 Somerville P. G., 2003. Magnitude scaling of the near fault rupture directivity pulse. Physics of the Earth and Planetary Interiors, 137(1—4): 201—212. Wang N. , Elgamal A. , Shantz T. , 2017. Recorded seismic response of the Samoa Channel Bridge-foundation system and adjacent downhole array. Soil Dynamics and Earthquake Engineering, 92: 358—376. doi: 10.1016/j.soildyn.2016.09.034 Yan G. Y. , Chen F. Q. , 2015. Seismic performance of midstory isolated structures under near-field pulse-like ground motion and limiting deformation of isolation layers. Shock and Vibration, 2015: 730612. Yang S. , Mavroeidis G. P. , 2018. Bridges crossing fault rupture zones: a review. Soil Dynamics and Earthquake Engineering, 113: 545—571. doi: 10.1016/j.soildyn.2018.03.027 -