多求解器FEA三维建模方法用于隧道稳定性评估

《International Journal of Geo-Engineering》：Multi solver FEA 3D modelling approach for tunnel stability assessment

　　

数值模拟广泛应用于地下工程中，以评估开挖引起的变形以及隧道与支护系统之间的相互作用。在本研究中，研究人员采用三维有限元方法分析了隧道稳定性，该方法结合了地应力初始化（geostatic initialization）与隐式-动态准静态步骤（implicit–dynamic quasi-static step）。该技术增强了屈服后收敛性，并有效捕捉了开挖过程中的瞬态力学响应。该框架应用于一个砂岩隧道断面，预测位移与现场监测数据进行了验证，在三个拱顶监测点显示出2–12%的相对误差一致性。网格敏感性分析表明，减少网格细化仅导致精度轻微下降，证实了建模框架的稳健性，关键输出差异仅为2.5%，从而确认了网格收敛性和稳定性。还使用隧道拱顶的节点应变能（strain-energy）输出进行了基于能量的评估。将释放的岩体能量与支护系统的工作能力进行比较表明，支护在其能量吸收极限内运行良好。此外，观察到的瞬态能量峰的逐步减小反映了以支护屈服行为为主的力学稳定化。总体而言，研究结果证明了多求解器准静态框架用于隧道稳定性评估的有效性，并为岩石-支护相互作用机制提供了有价值的见解，这些见解与工程设计及未来方法论发展相关。

**论文解读：多求解器FEA三维建模方法用于隧道稳定性评估**

**研究背景与问题**
在地下工程中，隧道稳定性评估是核心挑战，尤其在岩体具有非均匀强度、变形和能量耗散特性的开挖场景中。传统数值模拟方法（如静态或准静态求解）虽然能估算应力重分布和开挖后变形，但难以充分捕捉瞬态现象（如能量快速释放、局部屈服或动态稳定化），这些因素可能显著影响失稳的触发。近年来，求解器选择和本构模型的进步使得更全面的隧道-支护相互作用分析成为可能，但现有研究多关注单一求解器，缺乏对地应力初始化（geostatic initialization）与隐式-动态准静态（implicit–dynamic quasi-static）步骤结构化交替使用的方法学评估。此外，在数据有限的大型交通基础设施项目中，本构校准、求解器收敛性和稳定性指标解读存在不确定性。本研究旨在通过一种实用的多求解器准静态方法，结合能量分析，为隧道稳定性评估提供可复现的工程框架，填补方法学空白。

**研究内容与结论**
研究人员针对Texanna双隧道（位于阿尔及利亚东西高速公路，总长2 km，埋深130 m，穿过Albo–Aptian复理石单元，主要为砂岩与泥岩互层），采用三维有限元模型（FEA）在Abaqus中模拟隧道开挖。模型使用Mohr–Coulomb本构定律（Mohr–Coulomb (MC) model）代表砂岩的剪切主导行为，支持系统由HEB-180钢拱和0.25 m厚C30喷射混凝土组成。研究通过交替使用地应力初始化步骤和隐式-动态准静态步骤，模拟了分阶段开挖，并将预测位移与三个拱顶监测点的现场数据对比，相对误差在2–12%之间，一致性达98.6%（隐式-动态）和87.23%（地应力）。网格敏感性分析显示，减少约1000个元素（从5156到4156）后，位移输出仅变化2.5%，证实模型稳健性。能量分析基于隧道拱顶节点535的应变能输出，支持系统理论工作能力约1,002,016.6 J，而节点应变能峰值（约-42,000 J和21,000 J）远低于该值，表明支护系统在其能量吸收极限内运行良好。瞬态能量峰逐步减小，反映了屈服行为主导的力学稳定化。该研究发表在《International Journal of Geo-Engineering》。

**关键技术方法**
（1）三维有限元模型（3D FEA model）在Abaqus中建立，尺寸200 m×200 m×20 m，采用自由网格划分（medial axis算法），共5,156个元素，支持系统简化为板单元。（2）地应力初始化（geostatic initialization）步骤用于建立初始应力场（K?条件），确保静力平衡。（3）隐式-动态准静态步骤（implicit–dynamic quasi-static step）用于模拟开挖，控制加载速率使惯性效应可忽略。（4）能量评估基于节点应变能（nodal strain-energy）输出，比较释放岩体能量与支护系统工作能力。（5）网格敏感性分析通过减少约1,000个元素验证收敛性。（6）参数敏感性分析对杨氏模量（Young's modulus, E）、黏聚力（cohesion, c）和内摩擦角（friction angle, φ）独立降低15%，评估模型稳健性。现场监测数据来自三个拱顶监测点，样本为砂岩隧道段。

**研究结果**
**数值输出（Numerical outputs）**
通过隐式-动态和地应力求解器获得的垂直位移最大值分别为4.183 cm和4.536 cm。与三个监测点实测位移（4.9 cm、4.1 cm、4.6 cm、4.0 cm等）对比，一致性范围为85.36%–98.6%。网格减少1,000个元素后，位移变化极小（4.032 cm和4.18 cm），证实网格收敛性。参数敏感性分析表明，杨氏模量对变形影响最大，黏聚力和内摩擦角主要影响塑性变形范围，但参数变化15%后数值响应保持稳定，确认方法稳健性。

**能量分析（Energy analysis）**
节点535的应变能输出显示两个明显峰值（约-31,000 J和-42,000 J，对应正峰值约20,000 J和21,000 J），能量历史整体呈稳态演化后逐渐稳定，峰值均远低于支护系统工作能力（1,002,016.6 J）。瞬态能量峰逐步减小反映支护屈服行为导致的力学稳定化，而非惯性效应。

**结果讨论（Results discussions）**
数值预测与现场监测的良好一致性（98.6%和87.23%）表明多求解器交替工作流程可提供可靠的变形预测。网格敏感性分析中仅2.5%的差异证明模型在运输尺度评估中足够稳健。能量分析表明支护系统在能量吸收极限内运行，瞬态能量衰减反映稳定化机制。该框架在数据有限条件下仍能有效评估隧道稳定性，尤其适用于缺乏详细实验室测试的工程项目。

**总结与结论**
本研究通过结构化多求解器三维有限元框架评估了隧道稳定性。数值预测与现场监测的位移一致性高达约98%，求解器交替展示了准静态开挖条件下一致的变形趋势。网格敏感性分析表明，即使元素数量显著变化，数值输出保持稳定，证明采用的离散化方案对于运输尺度工程评估足够稳健。基于局部节点应变能演化与支护系统理论工作能力比较的能量解读表明，安装的支护系统在其预定能量吸收能力内运行。本研究仍受限于数据有限的工程环境，缺乏详细的岩体各向异性实验室表征和先进本构校准。因此，该框架应被视为一种实用的工程方法，用于初步稳定性评估和求解器验证，而非替代全面的现场特定岩土力学表征。总体而言，所提出的建模策略代表了一种实用且与现场对齐的框架，适用于大型交通基础设施项目中常见的数据有限条件下的隧道稳定性评估。其可复现性、低数据需求以及与现场测量结果的强验证使其成为设计验证和风险导向决策的合适工具。未来扩展可结合各向异性或应变软化本构模型，以获得更多现场特定数据。