Analytical Solution for the Longitudinal Response of Tunnels under Combined Seismic-fault Misalignment
-
摘要: 隧道穿越地震活动断层带时可能遭受严重破坏,以内马铁路一期三标段工程为依托,对地震动和断层错动联合作用下隧道结构纵向响应进行研究。针对穿越断层破碎带的隧道结构,基于地下结构抗震拟静力法,将其简化为弹性地基梁,并将地震动和断层错动位移简化为作用在隧道上的静荷载,建立地震动和断层错动联合作用下隧道纵向响应理论模型并进行求解,利用有限元数值模拟方法验证解析解的正确性。通过解析解进行参数敏感性分析,揭示断层错动位移、两侧围岩地基系数与断层破碎带地基系数比及围岩与隧道结构刚度比对隧道结构纵向响应的影响规律。研究结果表明,断层错动位移增加使隧道结构截面弯矩、剪力峰值接近线性增加,隧道内力沿隧道纵向分布的影响范围不变,且断层破碎带界面出现了截面剪力突变;两侧围岩地基系数与断层破碎带地基系数比对隧道结构截面剪力的影响较大,在地震动和断层错动联合作用下隧道结构剪力在断层破碎带界面急剧减小;随着围岩与隧道结构刚度比的减小,地震动引起的隧道挠度减小,断层错动作用引起的隧道挠度变化范围增大,同时隧道结构内力响应明显增大。Abstract: China is one of the countries that has the most severe seismic hazards in the world, and tunnels may suffer severe damage when crossing seismically active fault zones. In this paper, the longitudinal response of the tunnel structure under the joint action of earthquake and fault misalignment is investigated based on the third standard section of the first phase of the Nairobi-Malaba railway as a project. For the tunnel structure crossing the fault fragmentation zone, based on the idea of seismic anthropomorphism of underground structure, it is simplified into elastic foundation beam, and the ground vibration and fault misalignment displacement are simplified into static load acting on the tunnel, the theoretical model of tunnel longitudinal response under the joint action of earthquake - fault misalignment is established and solved, and the correctness of the analytical solution is verified by using the finite element numerical simulation method. Parametric sensitivity analysis is carried out to reveal the influence of the fault misalignment displacement, the ratio of the foundation coefficient of the surrounding rock to the foundation coefficient of the fault fragmentation zone, and the ratio of the surrounding rock to the tunnel structure stiffness on the longitudinal response of the tunnel structure. The results show that the increase of fault misalignment displacement makes the tunnel structure cross-sectional moment and tunnel structure cross-sectional shear peak close to linear increase, the influence range of tunnel internal force distribution along the tunnel longitudinal direction remains unchanged, and there is a sudden change of cross-sectional shear at the interface of fault fragmentation zone; the ratio of surrounding rock foundation coefficient and fault fragmentation zone foundation coefficient on both sides has a greater influence on the tunnel structure cross-sectional shear, under the combined effect of earthquake - fault misalignment tunnel The shear force of the tunnel structure decreases sharply at the interface of the fault zone; the ratio of surrounding rock to tunnel structure stiffness decreases, the tunnel deflection caused by the earthquake decreases, and the range of deflection caused by fault misalignment increases, while the internal force response of the tunnel structure increases significantly.
-
Key words:
- Tunnel /
- Seismic shaking /
- Fault dislocation /
- Analytical solution /
- Foundation beam /
- Quasi-static method
-
引言
颗粒阻尼器作为被动阻尼器,具有耐久性好、可靠度高、操作方便、安装简单等优点(鲁正等,2013)。颗粒阻尼技术的耗能机理为颗粒间碰撞耗能与主体结构冲击耗能(鲁正等,2014)。目前已有多位学者对颗粒阻尼器进行试验研究、理论分析与数值模拟。鲁正等(2012)通过3层钢框架附加颗粒阻尼器的振动台试验表明附加较小质量比(2.25%)的颗粒阻尼器可减小主体结构的响应。Fang等(2018)引入致密颗粒流理论,建立了NOPD(非阻塞型颗粒阻尼器)能量耗散定量模型理论,结果表明,NOPD能量耗散率随振动强度的增大而增大,随粒径的增大而减小。Hu等(2008)通过使用离散元法模拟发现质量比和材料密度均对阻尼性能有一定影响。Luo等(2019)提出多层颗粒阻尼器应用于缩尺比例为1:20的单塔自锚式悬索桥,发现多层颗粒阻尼器可在不同烈度地震激励下实现显著减震效果。在Papalou等(2015)开展的试验中,以颗粒阻尼器替换帕特农神庙中破损的石鼓,研究了质量比(颗粒质量与古柱之比)、阻尼器位置与颗粒尺寸等参数对颗粒阻尼器有效性的影响。
由于颗粒间碰撞属于高度非线性行为,目前关于颗粒阻尼器的数值模拟分2种:一种是使用离散元素法单独研究颗粒的力学行为;另一种是(施卫星等,2013)输入颗粒阻尼器的减震参数,进行有限元分析时,将多颗粒简化为单个颗粒,忽略了颗粒间摩擦耗能。本文在Papalou等(2015)试验的基础上,应用PFC3D、FLAC3D软件实现离散-有限的耦合。PFC3D是基于离散单元法模拟和分析颗粒体的软件,计算原理主要是力-位移定律和牛顿第二运动定律,采用显式有限差分法进行循环迭代求解,通过调幅可考虑不同地震动强度及不同地震动强度下结构反应标准差的变化规律(贾晗曦等,2019)。FLAC3D为有限元分析软件,采用显式拉格朗日算法和混合-离散分区技术,可准确发现模拟材料的塑性破坏和流动,由于无须形成刚度矩阵,因此,采用较小的计算资源可求解大范围的三维工程问题。本文首次通过PFC3D与FLAC3D软件耦合模拟附有多颗粒阻尼器的帕特农神庙多鼓石柱,充分考虑颗粒间摩擦与碰撞耗能,并探究地震强度和频率对阻尼器抗震性能的影响及阻尼器位置对结构体系抗震的影响,研究结果对于我国类似碑型古建筑抗震保护具有借鉴意义。
1. 耦合方法
PFC3D与FLAC3D软件各自以插件方式运行于彼此的环境中,采用界面耦合的方法模拟,界面耦合的基本原理是在可能发生耦合作用的空间建立1个单元wall,wall可作为2个空间信息传递的媒介,实现连续-非连续耦合作用。附着于FLAC3D模型元素上的1个wall单元如图1所示,图中GH表示位于PFC3D模型组件上的接触点,同时定义G为位于wall上对应于GH的接触点,由于接触部位存在变形,同时存在拉伸、剪切及扭转作用,因此GH与G点空间位置可能存在不一致的情况。
$ {b}_{i} $ 表示wall角点及其附着部位FLAC3D的zone结点或结构单元结点的坐标,$ {a}_{i} $ 为3个表征三角形的面积。定义
$ {r}_{i} $ 为GH与G的距离($ {r}_{i}=GH-{b}_{i} $ ),接触传递至FLAC3D模型元素的总接触力与总弯矩(PFC3D模型中的接触可承受弯矩)分别为F和M,其中M按下式计算:$$ M={M}_{\mathrm{b}}+(G-GH)\times F $$ (1) 式中,
$ {M}_{b} $ 、$ (G-GH)\times F $ 分别为接触中的弯矩及接触位置错位产生的弯矩。在迭代过程中,F、M均为已知量。由此,FLAC3D模型元素在wall附着部位即角点
$ {b}_{i} $ 处,承受的来自于接触耦合作用GH的集中力$ {F}_{i} $ 和弯矩$ {M}_{i} $ 计算如下:$$ {F}_{i}={a}_{i}/\sum _{i=1}^{3}{a}_{i}F $$ (2) $$ {M}_{i}={r}_{i}F $$ (3) 2. 试验结果验证
2.1 计算模型与验证方案
帕特农神庙位于希腊雅典卫城的最高处石灰岩的山岗上,是卫城最重要的主体建筑,建成于公元前432年,神庙平面呈长方形,由46根高34 英尺(9.88 m)大理石柱构成,现仅留有1座石柱林立的外壳,如图2(a)所示。首先使用FLAC3D以zone单元对神庙古柱结构进行建模,然后导入PFC3D使用ball单元建立颗粒阻尼器模型,其次在古柱结构与颗粒阻尼器界面处以wall耦合交互信息,最后在FLAC3D中分析结构模型抗震性能,并与Papalou等(2015)试验结果进行对比。在此基础上,分析地震强度、频率及阻尼器位置对结构抗震性能的影响。
Papalou等(2015)试验中将颗粒填入石制空鼓,替换已破损的石鼓,分析替换后颗粒阻尼器对帕特农神庙古柱减震的影响,研究表明,选择适当的参数组合可降低30%以上的古柱动态响应。与额外增加同等质量配重相比较,增加颗粒可以更加有效地控制共振响应(张奎等,2017)。为便于比较分析,本文数值模型与Papalou等(2015)试验保持一致,即柱高为2992 mm,自重为1707 kg,由11个直径不同的空鼓组成,顶鼓直径为445 mm,底鼓直径为584 mm,鼓高均为272 mm,柱身材料为大理石,柱体底座为
$ 900\;\mathrm{m}\mathrm{m}\times 700\;\mathrm{m}\mathrm{m}\times 140\;\mathrm{m}\mathrm{m} $ 大理石板,如图2(b)所示。古柱PFC3D与FLAC3D耦合仿真模型如图2(c)所示。2.2 验证结果
依据Papalou等(2015)试验模型,输入Kalamata地震波,如图3所示,古柱顶鼓中心点为荷载作用下的监测位点。将古柱顶鼓替换为包含颗粒阻尼器的空鼓,阻尼器颗粒粒径均为50 mm,颗粒个数为32个,颗粒为钢球。数值计算与试验得到的位移时程曲线如图4(a)所示,由图4(a)可知,数值计算结果与试验结果基本一致。数值计算与试验得到的加速度功率谱如图4(b)所示,由图4(b)可知,数值计算结果与试验结果基本一致,最大误差为17.5%。综上所述,PFC3D与FLAC3D耦合方法可模拟颗粒阻尼器与结构模型,实现连续-非连续分析,可充分考虑颗粒间的相互作用,真实有效地分析颗粒阻尼器抗震性能,可用于分析颗粒阻尼器参数对结构减震性能的影响。
3. 柱抗震性能影响因素分析
3.1 地震强度的影响
基于前述模型及参数,将Kalamata地震波加速度幅值调整为0.05 g、0.1 g、0.2 g、0.245 g,作为结构底部振动输入。将古柱顶鼓替换为包含颗粒阻尼器的空鼓,颗粒粒径为50 mm,颗粒个数为32个,颗粒为钢球。不同地震强度激励下顶鼓滞回曲线(柱顶点抗力-位移曲线)如图5所示,由图5可知,随着地震强度的加强,滞回曲线包围的面积越大,说明颗粒阻尼器消耗结构能量越大;不同加速度下滞回曲线整体呈梭形且形状相似,当激励幅值较小时,古柱变形能力较好,表现出较强的抗震性能。
取每个加速度幅值下20个振动循环,求得加速度幅值为0.05 g时的等效黏滞阻尼系数为13.03%~15.27%,加速度幅值为0.1 g时的等效黏滞阻尼系数为16.50%~19.88%,加速度幅值为0.2 g时的等效黏滞阻尼系数为20.55%~23.34%,加速度幅值为0.245 g时的等效黏滞阻尼系数为21.17%~24.51%,可知随着加速度幅值的增加,等效黏滞阻尼系数增加,结构耗能能力增强。
不同地震强度下结构顶鼓加速度和位移响应如表1所示,由表1可知,附加颗粒阻尼器古柱加速度均值响应和位移均值响应均小于未附加颗粒阻尼器古柱;加速度均方根减震率为(未附加颗粒阻尼器古柱响应−附加颗粒阻尼器古柱响应)/未附加颗粒阻尼器古柱响应,最大值为39.6%,最小值为13.9%,说明颗粒阻尼器在不同地震强度下均表现出良好的耗能能力;地震强度对结构加速度均方根减震率的影响较大,对位移均方根减震率的影响较小。
表 1 不同地震强度下结构顶鼓加速度和位移响应Table 1. Acceleration and displacement response of the structure with different earthquake intensity地震波
类型加速度
幅值/g阻尼器
位置加速度
峰值/$ \mathit{g} $加速度
均方根/$ \mathit{g} $位移峰值/
mm位移均
方根/mm无阻尼器
加速度均方根/$ \mathit{g} $无阻尼器位移
均方根/mm加速度均方根
减震率/%位移均方根
减震率/%Kalamata1986 0.245 顶鼓 1.044 49 0.212 76 13.526 0 1.926 0.264 42 2.889 0 13.9 33.3 0.200 顶鼓 0.694 80 0.100 40 7.717 2 1.718 0.166 30 2.583 9 39.6 31.0 0.100 顶鼓 0.458 66 0.084 55 5.885 0 0.942 0.114 54 1.318 8 26.2 28.5 0.050 顶鼓 0.250 91 0.047 46 3.451 0 0.093 0.064 07 0.126 2 25.9 25.8 3.2 输入地震波的影响
分别将宁河波、Kobe波、Kalamata波加速度幅值调整为0.2 g,地震波傅里叶幅值谱如图6所示,作为结构底部振动输入。经分析,宁河波主频为0.732 Hz。Kobe波主频为1.4 Hz,Kalamata波主频为2.614 Hz,古柱的前4阶自振频率分别为1.2、1.44、1.86、2 Hz。将顶鼓替换为包含颗粒阻尼器的空鼓,颗粒参数同上节。
损耗能量(输入总能量-弹性振动能)、加速度(由宁河波、Kobe波、Kalamata波作用下结构顶部19 s内加速度时程响应)与频率三维图如图7所示。由图7可知,当频率低且加速度较小时,结构损耗能量低,当频率低且加速度较大时,结构损耗能量高;低频区结构损耗能量低,高频区结构损耗能量高;当入射波主频频率增大时,三维图像趋于平缓,加速度对阻尼器的影响减小,阻尼器对结构减震作用增强。
不同激励频率下结构顶鼓加速度和位移响应如表2所示,由表2可知,不同地震波作用下附加颗粒阻尼器古柱加速度均值响应和位移均值响应均小于未附加颗粒阻尼器古柱;不同地震波作用下,颗粒阻尼器对古柱的减震效果不同,Kalamata波作用下古柱加速度均方根减震率和位移均方根减震率最大,宁河波作用下古柱加速度均方根减震率和位移均方根减震率最小。综上所述,颗粒阻尼器对入射波频率敏感,对结构减震效果的影响显著。
表 2 不同激励频率下结构顶鼓加速度和位移响应Table 2. Acceleration and displacement responses of the structure with different excitation frequencies地震波
类型加速度
幅值/g阻尼器
位置加速度
峰值/$ \mathit{g} $加速度
均方根/$ \mathit{g} $位移峰值/
mm位移均
方根/mm无阻尼器加速度
均方根/$ \mathit{g} $无阻尼器位移
均方根/mm加速度均方根
减震率/%位移均方根
减震率/%宁河波 0.2 顶鼓 0.589 5 0.089 9 7.100 8 1.621 0 0.119 3 2.133 0 24.6 24.0 Kobe波 0.2 顶鼓 0.621 3 0.094 6 7.487 8 1.513 7 0.146 8 2.300 8 35.5 34.2 Kalamata波 0.2 顶鼓 0.694 8 0.100 4 7.717 2 1.718 2 0.166 3 2.583 9 39.6 31.0 3.3 颗粒阻尼器的影响
将颗粒阻尼器分别置于1~11层鼓处,分析阻尼器位置对古柱减震效果的影响,入射波为加速度幅值0.2 g的Kalamata波,颗粒参数同上节。加速度响应如图8所示,位移响应如图9所示。由图8、9可知,阻尼器位置对减震效果的影响显著,在11层鼓处布置颗粒阻尼器对结构峰值加速度及均值加速度的控制效果最优;在7~11层鼓处布置颗粒阻尼器对结构峰值位移、均值位移的控制效果优于在1~6层鼓处布置阻尼器,这与Papalou等(2015)的试验结果一致;另外,在1、2层布置颗粒阻尼器时,结构均值位移与峰值位移明显大于阻尼器布置在3~11层时,可能因为在地震波激励作用下大理石柱具有一定滤波效果,古柱下部频率高于上部,因此阻尼器布置在1、2层鼓处时增大了结构响应,因此产生较大位移。应在1~6层鼓处减少颗粒阻尼器数量,以防形成鞭梢效应,造成控制效果较差。
A、B古柱形成双柱结构体系(图10),入射波为加速度幅值0.2 g的Kalamata波,颗粒阻尼器布置于A柱11鼓处,颗粒参数同上节。A柱峰值加速度及位移如图11所示,B柱峰值加速度及位移如图12所示。由图11、12可知,考虑双柱结构体系,A、B柱峰值加速度及位移均有所减低,且A柱减震效果优于单柱,可知替换破损的空鼓(由前文可知在7~11鼓处布置阻尼器优于在1~6鼓处布置阻尼器)时,采用结构整体分析更符合实际情况。这与(Veeramuthuvel等,2016)实验在飞行器给定位置设置颗粒阻尼器,在较远点位仍具有减震效果结论一致。
4. 结论
本文应用PFC3D与FLAC3D耦合技术,研究颗粒阻尼器对帕特农神庙多鼓古柱抗震性能的影响,得出以下结论。
(1)数值分析结果与试验结果吻合度较高,表明PFC3D与FLAC3D耦合方法可用于分析颗粒阻尼器对结构减震性能的影响。
(2)地震强度较小时,颗粒阻尼器仍可产生较好的减震效果。不同地震波激励对颗粒阻尼器性能的影响显著。
(3)设计颗粒阻尼器时,采用结构整体分析更符合实际情况。颗粒阻尼器布置在结构中上部抗震效果更好。
-
表 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 -
甘星球, 徐锋, 王晓伟等, 2021. 断层错动隧道地震响应规律及减震模拟研究. 现代隧道技术, 58(3): 100—106Gan X. Q. , Xu F. , Wang X. W. , et al. , 2021. Simulation study on seismic response laws and seismic mitigation measures of tunnels under fault dislocation. Modern Tunnelling Technology, 58(3): 100—106. (in Chinese) 耿萍, 吴川, 唐金良等, 2012. 穿越断层破碎带隧道动力响应特性分析. 岩石力学与工程学报, 31(7): 1406—1413Geng P. , Wu C. , Tang J. L. , et al. , 2012. Analysis of dynamic response properties for tunnel through fault fracture zone. Chinese Journal of Rock Mechanics and Engineering, 31(7): 1406—1413. (in Chinese) 何川, 李林, 张景等, 2014. 隧道穿越断层破碎带震害机理研究. 岩土工程学报, 36(3): 427—434He C. , Li L. , Zhang J. , et al. , 2014. Seismic damage mechanism of tunnels through fault zones. Chinese Journal of Geotechnical Engineering, 36(3): 427—434. (in Chinese) 刘国钊, 乔亚飞, 何满潮等, 2020. 活动性断裂带错动下隧道纵向响应的解析解. 岩土力学, 41(3): 923—932Liu G. Z. , Qiao Y. F. , He M. C. , et al. , 2020. An analytical solution of longitudinal response of tunnels under dislocation of active fault. Rock and Soil Mechanics, 41(3): 923—932. (in Chinese) 刘学增, 王煦霖, 林亮伦, 2013.75°倾角正断层黏滑错动对公路隧道影响的模型试验研究. 岩石力学与工程学报, 32(8): 1714—1720Liu X. Z. , Wang X. L. , Lin L. L. , 2013. Model experiment on effect of normal fault with 75° dip angle stick-slip dislocation on highway tunnel. Chinese Journal of Rock Mechanics and Engineering, 32(8): 1714—1720. (in Chinese) 孙飞, 张志强, 易志伟, 2019. 正断层黏滑错动对地铁隧道结构影响的模型试验研究. 岩土力学, 40(8): 3037—3044, 3053Sun F. , Zhang Z. Q. , Yi Z. W. , 2019. Model experimental study of the influence of normal fault with stick-slip dislocation on subway tunnel structure. Rock and Soil Mechanics, 40(8): 3037—3044, 3053. (in Chinese) 王明年, 崔光耀, 2011. 高烈度地震区隧道设置减震层的减震原理研究. 土木工程学报, 44(8): 126—131Wang M. N. , Cui G. Y. , 2011. Study of the mechanism of shock absorption layer in the supporting system of tunnels in highly seismic areas. China Civil Engineering Journal, 44(8): 126—131. (in Chinese) 文鑫涛, 李华玥, 段乙好等, 2021. 2020年中国大陆地震灾害损失述评. 震灾防御技术, 16(4): 651—656.Wen X. T., Li H. Y., Duan Y. H., et al., 2021. Earthquake disasters loss on Chinese mainland in 2020. Technology for Earthquake Disaster Prevention, 16(4): 651−656. (in Chinese) 闫高明, 申玉生, 高波等, 2019. 穿越黏滑断层分段接头隧道模型试验研究. 岩土力学, 40(11): 4450—4458Yan G. M. , Shen Y. S. , Gao B. , et al. , 2019. Experimental study of stick-slip fault crossing segmental tunnels with joints. Rock and Soil Mechanics, 40(11): 4450—4458. (in Chinese) 张景, 何川, 耿萍等, 2017. 穿越软硬突变地层盾构隧道纵向地震响应振动台试验研究. 岩石力学与工程学报, 36(1): 68—77 doi: 10.13722/j.cnki.jrme.2016.0103Zhang J. , He C. , Geng P. , et al. , 2017. Shaking table tests on longitudinal seismic response of shield tunnel through soft-hard stratum junction. Chinese Journal of Rock Mechanics and Engineering, 36(1): 68—77. (in Chinese) doi: 10.13722/j.cnki.jrme.2016.0103 Anastasopoulos I. , Gerolymos N. , Drosos V. , et al. , 2008. Behaviour of deep immersed tunnel under combined normal fault rupture deformation and subsequent seismic shaking. Bulletin of Earthquake Engineering, 6(2): 213—239. doi: 10.1007/s10518-007-9055-0 Baziar M. H. , Nabizadeh A. , Lee C. J. , et al. , 2014. Centrifuge modeling of interaction between reverse faulting and tunnel. Soil Dynamics and Earthquake Engineering, 65: 151—164. doi: 10.1016/j.soildyn.2014.04.008 Fan L. , Chen J. L. , Peng S. Q. , et al. , 2020. Seismic response of tunnel under normal fault slips by shaking table test technique. Journal of Central South University, 27(4): 1306—1319. doi: 10.1007/s11771-020-4368-0 Hashash Y. M. A. , Hook J. J. , Schmidt B. , et al. , 2001. Seismic design and analysis of underground structures. Tunnelling and Underground Space Technology, 16(4): 247—293. doi: 10.1016/S0886-7798(01)00051-7 Liang L. J. , Xu C. J. , Zhu B. T. , et al. , 2020. Theoretical method for an elastic infinite beam resting on a deformable foundation with a local subsidence. Computers and Geotechnics, 127: 103740. doi: 10.1016/j.compgeo.2020.103740 Shen Y. S. , Wang Z. Z. , Yu J. , et al. , 2020. Shaking table test on flexible joints of mountain tunnels passing through normal fault. Tunnelling and Underground Space Technology, 98: 103299. doi: 10.1016/j.tust.2020.103299 St John C. M. , Zahrah T. F. , 1987. Aseismic design of underground structures. Tunnelling and Underground Space Technology, 2(2): 165—197. doi: 10.1016/0886-7798(87)90011-3 Tsinidis G. , De Silva F. , Anastasopoulos I. , et al. , 2020. Seismic behaviour of tunnels: from experiments to analysis. Tunnelling and Underground Space Technology, 99: 103334. doi: 10.1016/j.tust.2020.103334 Vesić A. B. , 1961. Bending of beams resting on isotropic elastic solid. Journal of the Engineering Mechanics Division, 87(2): 35—53. doi: 10.1061/JMCEA3.0000212 Wang J. N., 1993. Seismic design of tunnels: a simple state-of-the-art design approach. New York: Parsons Brinckerhoff Quade & Douglas. Inc. Yan G. M. , Shen Y. S. , Gao B. , et al. , 2020. Damage evolution of tunnel lining with steel reinforced rubber joints under normal faulting: an experimental and numerical investigation. Tunnelling and Underground Space Technology, 97: 103223. doi: 10.1016/j.tust.2019.103223 Yu H. T. , Zhang Z. W. , Chen J. T. , et al. , 2018. Analytical solution for longitudinal seismic response of tunnel liners with sharp stiffness transition. Tunnelling and Underground Space Technology, 77: 103—114. doi: 10.1016/j.tust.2018.04.001 Zhen C. , Qian S. , Gui-Min Z. , et al. , 2022. Response and mechanism of a tunnel subjected to combined fault rupture deformation and subsequent seismic excitation. Transportation Geotechnics, 34: 100749. doi: 10.1016/j.trgeo.2022.100749 Zhong Z. L. , Wang Z. , Zhao M. , et al. , 2020. Structural damage assessment of mountain tunnels in fault fracture zone subjected to multiple strike-slip fault movement. Tunnelling and Underground Space Technology, 104: 103527. doi: 10.1016/j.tust.2020.103527 -