《Computers and Geotechnics》：Evaluation of tunnel face stability using slip curves based on lower bound theory
编辑推荐：
本文针对隧道开挖中掌子面稳定性评估的关键技术难题，提出了一种基于下限理论的创新分析方法。传统方法多假设破坏区内滑移线为直线，并需预设难以确定的侧向土压力系数，导致评估不够精确。本研究通过将破坏区内的滑移线建模为对数螺旋曲线，对梯形非对称微元体建立力与力矩平衡方程，并推导出微分方程进行求解，从而建立了无需侧向土压力系数的新评估体系。计算验证表明，该方法能够准确复现先前离心模型试验得到的掌子面稳定性指标，为盾构法和NATM法等隧道工程的设计与施工提供了更可靠的理论工具。
当我们乘坐地铁或驾车穿行隧道时，很少会想到支撑在我们头顶的岩土经历了怎样复杂的力学平衡。在隧道工程，无论是采用盾构法还是新奥法（NATM），开挖过程中掌子面的稳定性都是关乎工程安全与效率的核心。如果掌子面支撑压力不足，前方土体可能发生失稳坍塌；而压力过高，则可能造成不必要的工程浪费甚至引发地表隆起。因此，如何精确评估和设定平衡掌子面土体所需的支撑压力，是隧道工程师们孜孜以求的目标。
长久以来，评估这一压力的主流理论方法，大多借鉴了筒仓效应的思想，即假设破坏区由两条垂直的滑移线围成，形成一个类似“筒仓”的区域。然而，无论是先前的隧道模型试验还是离心机试验都表明，在松散地层中，实际的滑移破坏线往往是曲线而非直线。强行将其简化为直线，无疑会引入误差。更棘手的是，在传统的地压理论中，侧向土压力系数是一个关键的、但通常难以准确确定的参数。这两个问题——不真实的直线滑移线假设和不确定的侧向土压力系数——共同限制了现有理论评估的精度。此外，基于上限定理的数值分析方法虽然功能强大，但需要构建复杂的静力许可速度场并进行大规模矩阵运算，在实际工程设计中应用不便。面对这些挑战，亟需一种既能更真实反映物理现象，又便于工程应用的新评估方法。
为了应对上述挑战，日本埼玉大学的研究团队在《Computers and Geotechnics》上发表了一项研究，提出了一种基于下限定理、利用曲线滑移线评估掌子面稳定性的新方法。研究团队从基本的力学平衡原理出发，构建了一个由一条垂直滑移线（g）和一条等角对数螺旋曲线滑移线（h）围成的二维破坏区模型。在这个破坏区内，他们提取了一个梯形非对称微元体，对其建立竖向力、水平力以及力矩的平衡方程。通过对微元体上应力的张量方向进行分析，并结合滑移线上的应力状态关系（基于莫尔-库仑破坏准则），他们成功地推导出了控制该力学系统的微分方程。这个方法的巧妙之处在于，它无需预设侧向土压力系数，也无需像传统上限分析法那样进行复杂的空间离散化。通过数值求解这个微分方程，可以得到掌子面土压力的下限解，从而评估其稳定性。研究还采用了对数螺旋曲线来模拟外滑移线，这与土力学中考虑三轴剪切角的常规做法一致，使得模型更能贴合实际观测到的曲线破坏形态。
为了验证所提方法的合理性与适用性，研究人员进行了多组参数计算，并与先前的离心模型试验结果进行了对比。
方法的理论基础与模型建立
研究首先详细阐述了基于曲线滑移线的掌子面压力估算理论。该方法本质上是一种极限平衡法，其核心创新是用曲线来刻画与开挖相关的塌陷区。如图2所示，破坏区由掌子面上方的垂直滑移线g和一条等角对数螺旋线h围成。通过建立并求解微元体（如图3所示）的平衡方程组，可以导出描述现象的微分方程，从而构成极限平衡法的下限解。对数螺旋线的形状参数（m, z）由覆土厚度C、隧道直径D以及地表破坏区宽度L等几何条件决定，确保了曲线能真实地连接地表与掌子面下部边缘。
微元体的力学平衡分析
在确定的滑移线框架下，研究对一个位于破坏区内部、以滑移线g为上底、曲线h为下底的梯形非对称微元体进行了细致的力学分析。如图3所示，分别考虑了微元体上竖向力、水平力的平衡以及绕其几何中心的力矩平衡，并推导出相应的平衡方程（文中式6, 7, 8）。通过数学变换，这些方程最终被整理为关于角度θ的微分方程形式（文中式9, 10, 11），为后续数值求解奠定了基础。
滑移线上的应力状态关系
为了封闭方程组，研究引入了滑移线上的应力条件。基于莫尔-库仑破坏准则，作用在滑移线g和h上的剪应力与法向应力满足关系（文中式12）。更重要的是，通过分析滑移线g上的应力状态（如图4, 5, 6所示），并结合覆土压力为γC sinθ这一条件，可以推导出该滑移线上法向应力σ_g与剪应力τ_g之间的一个重要线性关系（文中式14），其中系数ζ_g = 2 cot(π/2 - φ)。这个关系式的引入，巧妙地避免了对难以确定的侧向土压力系数的依赖，这是本方法的关键优势之一。
计算验证与结果分析
通过数值求解上述微分方程组，研究对提出的方法进行了系统验证。计算结果表明，所提出的解能够正确反映摩擦材料由于覆土和摩擦阻力而产生的行为特性。尤为重要的是，计算得到的掌子面稳定性指标与先前一项离心模型试验所获得的结果能够充分且准确地吻合。这强有力地证明了，采用曲线滑移线模型能够更真实地模拟实际破坏模式，从而提高掌子面土压力评估的精度。
研究结论与重要意义
本研究成功开发并验证了一种用于评估隧道掌子面稳定性的新理论方法。该方法的核心贡献与重要意义体现在以下几个方面：
- 1.
模型的物理真实性提升：摒弃了传统理论中将所有滑移线视为垂直直线的简化假设，转而采用更符合实验观测（特别是在松散土层中）的曲线滑移线，尤其是对数螺旋线，来刻画破坏区边界，使理论模型更贴近实际物理现象。
- 2.
理论框架的简化与实用化：该方法基于下限定理，通过求解微元体平衡导出的微分方程来获得解，无需进行复杂的空间离散化或大规模矩阵运算，避免了传统上限分析法计算上的复杂性，使其更易于在工程设计和实践中应用。
- 3.
关键参数依赖的消除：通过引入滑移线上的应力状态关系，本方法在推导过程中完全不需要预先设定或估算侧向土压力系数这一在传统地压理论中的未确定参数，从而减少了评估过程中的主观性和不确定性。
- 4.
良好的验证效果：研究通过参数计算证实了该方法的合理行为，并且其计算结果与先前的离心机模型试验数据吻合良好，证明了该方法在准确预测掌子面稳定性指标方面的有效性。
综上所述，这项研究为隧道掌子面稳定性分析提供了一种新颖、严谨且更可靠的极限平衡方法。它不仅增强了我们对掌子面破坏机制的理论理解，而且为隧道工程，特别是盾构掘进中的泥水压力或土压平衡设定，提供了更具指导意义的理论工具，有助于实现更安全、更经济的设计与施工。未来，研究者计划将这一方法从目前的二维问题推广到更复杂的三维情况，以拓展其应用范围。
