各向异性岩体在静态与动态联合荷载作用下的隧道失稳机制

《Geohazard Mechanics》：Tunnel instability mechanisms in anisotropic rock masses under combined static and dynamic loading

　　

在本研究中，研究人员采用数值模拟方法，探讨了由不连续面存在及其产状引起的岩体各向异性对隧道在静态与动态联合荷载作用下行为的影响。通过节理岩体模型(JRM)模拟不同的节理倾角和地应力比，分析了隧道衬砌内力(TLIF)和位移响应。动态荷载以水平传播的剪切波加速度时程施加，峰值加速度为0.3g，主频代表中等地震激励。结果表明，各向异性显著影响应力重分布和变形模式，尤其是在低应力比(K=0.5)和节理倾角约为60°的情况下，弯矩和剪力相对于无节理的各向同性参考模型可增加高达250%。相反，轴向力和位移可能根据节理构型而减少。这些发现为节理岩体中地震诱发的不稳定机制提供了见解，并有助于提高对地质灾害相关动态荷载下隧道性能的理解。

**论文解读：各向异性岩体在静态与动态联合荷载作用下的隧道失稳机制**

**1. 研究背景与问题**

在自然界中，由于岩土颗粒的层状排列以及岩体中结构面（如裂隙、节理和断层）的存在，大多数岩体表现出各向异性特性。这种各向异性导致地下洞室周围应力状态和位移分布不均匀，产生不同的破坏深度，并可能在相对较低的应力水平和意外应力条件下引发应力诱发破坏（如近壁退化）。在地震活跃区域，这些机制可触发显著的不稳定性，包括局部坍塌、过度变形和衬砌超应力，构成强烈地面运动期间的主要地质灾害风险。

以往研究主要集中于静态荷载条件或假设围岩为各向同性，未能充分捕捉节理诱发各向异性或地震荷载对隧道响应的影响。例如，Wang等人提出了考虑全滑移界面的任意形状衬砌隧道解析解，但仅限于静水应力场；Xu等人通过现场调查和数值分析研究了层状岩层中二次衬砌的损伤，但未系统考虑动态荷载；Karimi-Khajelangi等人结合离散裂隙网络(DFN)和离散元法(DEM)分析各向异性对孔周应力分布的影响，但未涉及隧道地震响应。此外，现有研究很少系统地涵盖节理倾角和地应力比在地震荷载下的参数分析，限制了对隧道内力和塑性区演化综合影响的理解。因此，本研究旨在通过数值模拟，系统探究节理诱发各向异性与动态荷载的耦合效应，填补这一研究空白，为节理岩体中隧道的抗震设计提供新见解。

**2. 研究内容与意义**

研究人员采用节理岩体模型(JRM)模拟各向异性介质，评估方向性刚度、节理倾角和地应力比如何共同影响隧道变形和衬砌内力，为地震活跃、层状地质构造中的隧道设计提供新见解。研究结果表明，特定节理取向（尤其是低应力比K=0.5下的60°倾角）可导致剪力和弯矩显著放大，增幅超过250%。各向异性条件，特别是具有两组相交节理的介质，显著影响应力传递机制，导致隧道衬砌与围岩之间复杂且高度方向性的相互作用。敏感性分析确认，刚度比(En/Et)和节理强度参数（内聚力、摩擦角）在应力传递机制和变形行为中起关键作用。这些发现系统耦合了节理控制各向异性与动态荷载，为理解方向性刚度、节理激活和应力波传播如何控制隧道变形和破坏机制提供了新视角，对节理岩体中的抗震设计和灾害评估具有实际价值。

**3. 关键技术方法**

本研究主要采用数值模拟方法，使用节理岩体模型(JRM)模拟各向异性岩体。模型基于弹性-理想塑性本构关系，考虑正交各向异性刚度（法向En和平行Et）和节理面的Mohr-Coulomb强度准则。通过在不同节理倾角（0°、30°、45°、60°、90°）和地应力比（K=0.5、1、1.5）下，模拟隧道开挖和地震荷载施加，分析衬砌内力和位移响应。动态荷载以水平传播的剪切波加速度时程施加，峰值加速度0.3g，持续30秒。边界条件采用自由场边界和柔性地基，以吸收反射波。所有模型基于平面应变假设，材料参数保持一致性以确保结果可比性。

**4. 研究结果**

**4.1 单一节理组的影响**

通过模拟不同节理倾角模型，发现地应力比(K)对隧道响应有显著影响。在低K值（如K=0.5）下，节理对内力传递和减少的影响更强，变化幅度更大；而在高K值下，影响趋于均匀。节理倾角方面，剪切力和弯矩对倾角高度敏感，变化范围从-50%到+250%。最强效应出现在倾角60°和K=0.5时，轴向力减少约20%，而剪力和弯矩增加2.5倍。这表明节理优先耗散轴向应力，同时集中剪切和弯曲需求于衬砌。变形方面，倾角小于45°时，岩体刚度方向与水平运动对齐，减少变形；倾角大于45°时，弱方向对齐，增加变形。总体响应呈非线性非单调，源于弹性变形和沿节理塑性滑动的组合。

**4.2 敏感性分析**

针对各向异性参数（En/Et、内聚力、摩擦角）的敏感性分析，以K=0.5和倾角30°为代表案例。结果显示，随着En/Et增大，最大轴向力减小，弯矩增大；最大剪力出现在平行刚度是垂直刚度两倍时。内聚力增加减少了剪力和弯矩，但略增轴向力，且超过100 kPa后影响减弱。摩擦角增大（>15°）显著减少TLIF，但无剪切阻力（C=0, φ=0°）时，轴向力增加而剪力和弯矩分别减少50%和25%。这些规律表明，节理强度参数调节动态荷载下的滑动激活，高强度限制滑移，低强度则吸收地震能量但产生不对称变形。

**4.3 两组节理的影响**

模拟两组节理（不同倾角组合）显示，其影响可为正或负，取决于节理与地震波的相互作用。当两组节理倾角为45°和135°时，剪力、弯矩和位移增幅最大，但轴向力减少；最大轴向力增幅出现在15°和105°组合时，同时位移减少最小。当一组水平、另一组倾角小于60°时，TLIFs减少，位移轻微增加（<3%）。这表明节理方向性激活导致多方向应力重分布，特定组合可显著放大或减少内力。

**4.4 塑性区分析**

观察峰值地震荷载下K=0.5时的塑性区发展，发现各向异性显著影响塑性区形态。各向同性模型形成V形塑性区，延伸至表面；但水平节理衰减地震波，防止塑性区形成；45°节理产生极小塑性区；垂直节理使塑性区从隧道上部延伸至表面。两组节理配置中，当一组垂直且另一组倾斜时，塑性区扩大，促进连续剪切带形成；而两组近水平节理时，无塑性区发展。塑性区扩展方向与最大动剪切应变方向一致，而非节理平面。

**5. 讨论与结论**

讨论部分指出，隧道响应复杂性源于地应力比控制的应力重分布和节理诱导方向性刚度的交互作用。地震荷载作为水平剪切波施加，进一步增加动态交互。工程启示方面，忽略各向异性可能导致低估隧道变形和内力，特别是弱刚度方向与地震波传播方向对齐时。陡倾或垂直节理需特别关注，因其促进连续剪切带和扩大塑性区，威胁衬砌完整性。低K值下，各向异性岩石对地震响应更敏感，需保守评估不利节理取向和低刚度比。

研究结论部分翻译：本研究通过JRM数值模型，系统解决了节理诱导各向异性如何影响隧道内力、变形和塑性区发展的问题。结果表明，特定节理取向（如60°倾角）在低应力比下可显著放大剪力和弯矩，增幅超250%；两组节理介质的各向异性条件产生复杂方向性相互作用。敏感性分析确认刚度比和节理强度参数的关键作用。这些发现为抗震设计和灾害评估提供了新见解，但需注意研究基于平面应变假设，未来应纳入三维效应、节理连续性和多向地震输入，并扩展至深埋隧道或强震条件。