先进位移-约束法（ADCM）：软土地层闭合盾构隧道短期力学行为的半解析框架

《TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY》：The advanced displacement-confinement method for design of mechanized shield tunnels in soft ground

　　

研究人员提出了一种名为先进位移-约束法（Advanced Displacement-Confinement Method, ADCM）的半解析框架，用于评估软土地层中闭合式机械化圆形盾构隧道的短期土体与衬砌相互作用。该方法整合了多个独立分析算法，可直接估算盾构机（TBM）压力分布、地层变形及衬砌受力。其本构模型结合了弹塑性剪切硬化行为，通过SHIELD模型引入各向异性约束比λTBM(θ)，并利用RINGS模型计算考虑盾构机护盾相互作用的地层位移与应力场。ADCM避免了复杂的各向异性应力场和迭代本构模型，仅需标准勘察数据即可完成分析，计算效率显著优于三维数值模拟。通过五个实际工程案例验证，ADCM预测的沉降槽与衬砌受力与现场监测结果吻合良好，平均均方根误差约为1 mm。该方法特别适用于设计初期参数研究和实时施工监控，能够将地表沉降预警阈值直接与TBM操作参数相关联。

研究背景与意义

在城市地下空间开发中，软土地层盾构隧道施工引起的地层变形与控制是岩土工程的核心挑战。传统设计依赖经验公式或计算密集型的三维数值模拟，前者难以精确量化TBM操作参数（如支护压力、盾构间隙）对地层响应的影响，后者则耗时过长且参数校准复杂，无法满足施工实时决策需求。现有本构模型（如Hardening Soil模型）虽能模拟土体硬化行为，但需大量专项试验参数，且对OCR（超固结比）等初始条件敏感，限制了其在常规项目中的应用。为此，研究人员开发了ADCM框架，旨在建立从TBM操作参数到地表沉降的可逆分析链，实现高效、透明的短期力学行为预测。

关键技术方法

研究人员采用五个已竣工的软土盾构隧道案例进行验证，涵盖泥水盾构（Slurry TBM）和土压平衡盾构（EPBM），直径6.27~12.6 m，埋深16.9~33 m，地质条件包括黏土、粉土、砂土及含砾地层。核心技术包括：1）MCSH（Modified Cam-Clay with Shear Hardening）本构模型，在标准Mohr-Coulomb模型中引入HS模型的剪切硬化定律，仅需常规勘察参数（摩擦角φ、黏聚力c'、卸载模量E_ur）；2）SHIELD模型，基于各向异性约束比λ_TBM(θ)计算TBM支护压力分布；3）RINGS算法，通过径向增量法求解弹性-塑性分区应力与位移场，并引入护盾间隙限制最大变形；4）PCTL（预制混凝土管片）结构模型，将地层荷载转化为等效弹簧刚度，计算衬砌内力与收敛变形；5）体积损失解析算法，综合地层应变与衬砌收敛推导地表沉降曲线。

研究结果

2.2.2 弹性-理想塑性应力场

基于Lee (2018)的各向异性解析解近似，研究人员推导了弹性与塑性区的径向和切向应力表达式。弹性区采用各向同性形式（K₀=1），但引入SHIELD模型计算的各向异性卸荷比λ_TBM(θ)。塑性区通过Mohr-Coulomb强度准则与平衡微分方程求解，得到塑性半径r_p(θ)随卸荷比变化的显式公式。结果表明，当卸荷比超过临界值λ_TBM^e时，塑性区从隧道开挖面开始扩展。

2.2.3 硬化定律

研究人员在MCSH模型中嵌入HS模型的剪切硬化屈服面，通过归一化硬化函数r_q控制偏应力q的增长。硬化参数由参考压力p_ref（通常取1 atm）、塑性硬化参数γ_p和破坏比R_f（通常0.9）确定。与标准MC模型相比，MCSH在达到屈服前即产生塑性剪切应变，降低了峰值偏应力并增加了径向位移，同时避免了HS模型对OCR的敏感性。

2.2.4 流动法则

采用非关联流动法则计算塑性应变增量，通过径向积分累加弹性与塑性应变。剪切模量G由卸载模量E_ur和泊松比ν确定，径向位移u_r通过切向应变ε_θ与半径r的乘积获得。计算表明，选择远场半径r_∞≥20r₀和径向步长Δr≤0.25 m可保证精度，初始应力需设为竖向有效应力σ'_v以避免应变累积误差。

2.3 与MC和HS模型的比较

对比显示，MCSH模型介于MC与HS模型之间：在屈服前遵循MC模型的弹性应力路径，屈服后产生塑性应变，总变形量接近HS模型OCR=1~2的中间值。HS模型对初始OCR敏感，正常固结条件下（OCR=1）塑性区无限扩展，而MCSH无需考虑OCR即可获得稳定响应，更适合常规设计。

2.4 盾构相互作用

通过缩放Vlachopoulos-Diederichs纵向位移曲线（LDP）至盾构间隙Δ_shield，研究人员量化了护盾对变形的限制作用。当塑性半径R*超过临界值7.32时，地层全塑性失稳，该指标可作为掌子面稳定性预警阈值。算法在达到最大位移后停止应变累积，符合实际护盾-土体接触的非线性行为。

3. PCTL结构模型

基于RINGS输出的径向应力σ'_r,RINGS(θ)和TBM支护压力P'_SHIELD(θ)，研究人员计算了衬砌有效应力σ'_lining(θ)。平均竖向与水平荷载分别由拱顶与拱底应力均值及侧向土压力系数K_TBM确定，衬砌刚度采用RINGS输出的平均卸载模量E_s。通过Duddeck-Erdmann解析方程计算弯矩、轴力和收敛变形，验证了浅埋与深埋隧道的适用性。

4. 地层损失与地表移动

体积损失V_L由平均切向应变ε?_θ与衬砌净收敛u?/R之和解析获得，无需经验假设。结合O’Reilly-New半经验公式计算地表沉降槽，其中沉降槽宽度参数K_i与Rankine主动土压力系数K_a相关（0.3~0.5）。

4.1 位移-约束曲线（DCC）

研究人员绘制了地表最大沉降随TBM卸荷比（1-λ_TBM）的变化曲线。最优卸荷比范围为0.2~0.8，对应可控的地表沉降；低于0.2时塑性半径超过隧道埋深导致失稳，高于0.8时对沉降影响微弱但增加刀具磨损。DCC可直接将沉降预警与TBM操作参数关联。

5. 验证与案例研究

五个案例的实测与计算峰值沉降均方根误差为0.5~1.5 mm。反演得到的λ_TBM值在0.35~0.93之间，多数处于0.6~0.8的最优区间。案例4（λ_TBM=0.93）因欠压导致较大变形，但受小盾构间隙（15 mm）限制未引发失稳；案例5（λ_TBM=0.85）因大间隙（220 mm）产生近200 mm隧道收敛，但高注浆压力抑制了地表沉降。衬砌弯矩利用率与体积损失正相关，大直径隧道因严格控制沉降而受力更低。

结论与讨论

ADCM首次实现了TBM操作-地层响应-衬砌受力-地表沉降的可逆分析链，其半解析特性使计算时间缩短至秒级，显著优于三维数值模拟。研究证实，软土盾构短期沉降主要由掌子面约束控制，孔隙水压力消散与长期固结对绿field沉降贡献极小。该方法为设计阶段参数优化和施工期实时风险管控提供了量化工具，未来可结合监测系统实现动态反演与注浆 mitigation 策略优化。论文发表于《Tunnelling and Underground Space Technology》。