复合地层盾构开挖面极限支护力上限分析

2021-03-29代仲海胡再强

工程科学与技术 2021年2期

代仲海，胡再强

(西安理工大学岩土工程研究所，陕西西安710048)

随着城市化建设高速发展，开发地下空间资源已成为构建资源节约型城市的重要手段之一[1]，如依靠地铁建设缓解地上交通拥挤的情况。盾构法作为隧道施工的工法之一，以其施工速度快、周围环境不受干扰等优点而被广泛采用[2]。盾构施工常常遇见复合地层的情况，盾构在复合地层掘进过程中若土仓压力控制不当，可能导致开挖面失稳，严重破坏周边建筑物，并严重威胁施工人员生命安全。由于地下岩土体特性及受力情况的复杂，关于盾构开挖面稳定性的研究主要集中于数值模拟与理论解析研究方面。数值模拟的优势在于能够模拟盾构施工中复杂的地层情况及各类荷载，可为实际工程评估提供借鉴。Senent等[3]利用FLAC3D对上软下硬复合地层的极限支护压力与开挖面的软硬比进行研究，得出上部地层较软时，支护压力不足会导致开挖面局部破坏的结论。何小辉等[4]利用ABAQUS对上软下硬复合地层进行模拟，研究了隧道的变形特征。Zhang等[5]基于网格自适应有限元法，研究了上软下硬地层中矩形隧道开挖面的稳定性。理论解析法大致可分为极限平衡法和极限分析法。已有研究多基于极限平衡法。陈强[6]推导出上软下硬复合地层中极限支护压力的计算公式。李康[7]引入了地层复合角的概念，研究了上软下硬复合地层中盾构隧道施工引起的地表沉降及开挖面支护力。同时，极限分析因其严密的理论框架也被部分学者采用[8-9]。吕玺琳等[10]采用极限分析法推导了2维村山氏模型盾构开挖面极限支护压力最小值的计算公式。冯利坡等[11]采用极限分析上限法建立盾构隧道开挖面3维对数螺旋破面模型，并推导出极限支护压力计算公式。在极限分析的框架内，部分学者对盾构隧道上覆土体成层情况及土体黏聚力非均质的情况进行了研究[12-14]。

可见，对于盾构穿越复合地层开挖面稳定性解析方法的研究仍然较少。解析方法可为参数研究提供极大便利，因此，有必要对盾构穿越复合地层开挖面稳定性理论开展进一步研究。

支护压力计算方面，俞超杰等[15]运用FLAC3D软件对不同泥浆液位下盾构隧道开挖面主被动极限支护力进行研究，得出盾构推进过程中支护力非均匀分布与均匀分布的计算结果十分接近，故本文将盾构开挖面上的支护力简化为均布荷载。同时，学者们给出不同的破坏模式。Jancsecz等[16]采用楔形体极限平衡模型，计算最小极限支护力；Soubra[17]采用多块体破坏模型并运用极限分析得到了极限支护力上限解；Takano等[18]通过X射线断层扫描仪得出地层竖直破坏面类似于对数螺旋，故本文采用对数螺线破坏面作为复合地层盾构隧道开挖面破坏模式。

作者基于极限分析上限定理，构建了复合地层中盾构隧道开挖面对数螺旋破坏模式；将盾构开挖面上的支护力简化为均布荷载，计算了外力功率与内能耗损率，并利用非线性优化方法求得极限支护力的最优解，且与现有文献进行对比验证。着重讨论了地层参数、土体重度、隧道直径及深度系数对极限支护力及破坏模式的影响，并给出临界滑动面示意图。本文推导出的复合地层开挖面极限支护力解答简便而高效，为实际工程提供了一定参考。

1 计算原理和计算模型

1.1 计算原理

极限分析上限定理可表述为对于任意一个运动学许可的破坏机构，令其外力所做功率与内能耗散率相等，便能够得到破坏荷载或极限荷载的一个不安全的上限，所确定的荷载不小于真实极限荷载[19]，即：

1.2 计算模型

如图1所示，根据相关联流动法则的要求[20]，挖面前方塌落块体边界线由3条对数螺旋线AF、FE和EB组成，塌落块体AFEB以角速度 ω绕点O作刚体旋转破坏，点O为极坐标原点。隧道埋深为C，隧道直径为D；盾构穿越上部地层高度为nD，下部地层为（1-n）D，其中n为深度系数。上下部分地层抗剪强度参数分别为 c1、 φ1和 c2、 φ2。

图1 盾构开挖面破坏模型Fig.1 Failure mechanism of the shield excavation face

由此可得，对数螺旋线AF、FE和EB的表达式为：

2 外功率与内能耗散率

2.1 外功率

构建速度相容的破坏机构后，利用极限分析上限定理求解极限支护力的关键在于外功率和能耗的计算。本文所讨论的外力有土体重力和开挖面上的支护力。首先，计算土体重力功率，如图1所示，重力功率可由OEF、OEB区域功率之和减去OAF、OAB区域功率得到。考虑区域OEF，取一微元如图2所示。

图2 重力功率计算示意图Fig.2 Schematic diagram for calculating power of the gravity

盾构推进过程中支护压力简化为均匀分布，如图3所示。

图3 极限支护力功率计算示意图Fig.3 Schematic diagram of limit supporting forcepower calculation

图3中，P为开挖面上一点，相应极角为θ，由几何关系得：

2.2 内能耗散率

由于塌落块体作刚体旋转破坏，土体内部的能量耗散率忽略不计。将各速度间断面上的微分面积与速度间断线切线方向速度及土体黏聚力相乘，并沿整个间断面积分，可得AF、FE、BE上的内能耗损率：

则总的内能耗散率为：

2.3 极限支护力的求解

3 验证

表1 极限支护力对比结果Tab.1 Comparison resultsof ultimate support force

由表1可知，利用本文解答优化得到的结果非常接近文献[21]中的结果，最大误差不超过5%，由此验证了本文方法的可行性。

4 参数分析

为了解各因素对于盾构穿越复合地层稳定性的影响，计算不同参数下的开挖面极限支护力，并进行了相应的探讨。参数取值如下：φ1=φ2=20°，c1=c2=10 kPa ，D=8 m， C/D=1.5，γ =20 kN/m3。

4.1 地层强度参数

图4为不同深度系数n下，极限支护力随地层内摩擦角的变化。

由图4可知：当上部地层内摩擦角一定，随着下部地层内摩擦角的增大，极限支护力显著降低，这是由于土体自稳能力得到了改善；且深度系数越小，这一趋势越明显。深度系数为0.5时，上部地层内摩擦角为10°，下部地层内摩擦角由10°变化到30°，极限支护力减小56%，这比传统内摩擦角取均值时的计算结果减小了33%。同样地，当下部地层内摩擦角不变，随着上部地层内摩擦角增大，极限支护力逐渐降低。