淤泥区狭长基坑抗隆起稳定性分析

2023-01-26卢瀚

城市道桥与防洪 2022年12期

卢瀚

[上海市政工程设计研究总院（集团）有限公司，上海市200092]

0 引言

珠三角地区分布的深厚流塑状淤泥具有剪切强度低、含水量高、压缩性高、渗透性低的特点，是全国最软的软土之一[1]。软土的力学特性决定了基坑设计和施工难度大、造价高，且易出现工程事故。在软土地区建设综合管廊类狭长工程时，基坑的坑底抗隆起稳定成了开挖建设期间安全性的控制因素之一[2，3]。

传统的规范解析方法[4]采用考虑土体剪切强度指标、值的圆弧滑动模式来计算坑底抗隆起稳定性。但该方法不考虑狭长基坑两侧支护与土体的相互作用，计算得出的圆弧滑动面往往超过了基坑宽度，与实际情况不符。近年来，国内外学者针对不同因素对坑底抗隆起稳定性的影响进行了研究。张飞[2]通过离心机试验和数值分析得出基坑抗隆起稳定性应特别重视坑内外水位和墙体插入深度。俞建霖[5]通过对比现有解析法结合实际案例得出软弱土中规范要求的安全系数偏高并提出了考虑基坑宽度影响的抗隆起稳定性算法。王洪新[6]认为小窄基坑具有更好的稳定性，提出了一种考虑尺寸效应的基坑稳定性安全系数解析算法。郑刚[7]基于不连续布局优化法对土体强度在深度上的不均匀性对坑底抗隆起稳定性进行了研究，并指出应考虑窄基坑的空间效应对安全系数的提高作用。陈孝湘[8]基于强度折减法对狭长基坑的坑底隆起破坏模式进行了研究，得出内支撑的部位和基坑宽度对稳定性均有影响。

现基于强度折减法，采用土体硬化（HS）模型对珠三角深厚淤泥地区的某综合管廊基坑工程破坏模式进行了研究，分析了基坑宽度、深度、坑底加固等对坑底抗隆起稳定性的影响。

1 数值分析方法简介

数值分析方法[9]可根据土层分布的性质、基坑宽度对土体变形的影响、土体开挖和支撑架设的施工过程等综合分析，合理地描述基坑破坏时的性状。数值分析选取的本构模型和参数应能反应土体塑性性质，以准确模拟主动和被动土压力，从而合理地模拟基坑变形[10]。

1.1 土体硬化（HS）模型

HS模型为Plaxis软件中的一种本构模型，可描述多种土体的破坏和变形特征[9，10]。该模型可同时考虑土体的剪切硬化和压缩硬化，可适用于分析土体开挖卸荷类问题，已成为基坑数值分析中应用最多的模型之一。HS模型共有11个参数，包括刚度参数：三轴固结排水试验的参考割线模量Eref50、固结试验的参考切线模量Erefocd、三轴固结排水卸载再加载试验的参考卸载再加载模量Erefur、刚度应力水平参数m；强度参数：有效粘聚力c、有效内摩擦角φ、剪胀角ψ；高级参数：卸载再加载泊松比νur、参考应力pref、破坏比Rf、正常固结下侧压力系数K0。

本构模型中原状土剪切强度c和φ可通过地勘报告获得。根据Plaxis材料手册[11]：K0=1-sin（φ），pref=100 kPa，νur=0.2，黏性土的剪胀角ψ一般取为0，软土m取0.8，Rf=0.9。

对于刚度参数通常取Eref50=Erefocd，Erefur=3Erefocd，Erefocd≈Es1-2，其中Es1-2为压力区间100 kPa～200 kPa下的土体压缩模量。但对于软土地区，该刚度关系并非完全适用甚至会得到不合理的结果。王卫东等[12]对上海地区软土层HS模型刚度参数进行了试验研究和工程验证，并收集了多个国内外研究得出的参数；刘关虎等[13]对珠海深厚淤泥区的HS模型刚度参数进行了均匀试验和反分析，并汇总了国内外各地区刚度参数间及与压缩模量的关系。根据上述文献的研究成果和数据收集成果，初步得出软土层的硬化土刚度关系为Eref50=（0.5～1.8）Erefocd，Erefur=（2.9～13）Erefocd，Erefocd≈Es1-2。

坑底土体进行加固后性质会发生较大的变化，黄绍铭和高大钊[14]指出水泥土的参考割线模量Eref50与其无侧限抗压强度qu之间的大致关系为：Eref50=126qu。

1.2 强度折减法

有限元强度折减法最早应用于分析边坡稳定问题[15]，已有大量算例和分析研究表明该方法使用简单、准确，近年来在基坑工程数值分析中也得到了广泛应用。强度折减法的基本原理为逐步降低土层的剪切强度c和tan φ，代入有限元模型中进行计算，直到计算不收敛时，位移-强度折减系数曲线的转折点即为安全系数RF，其计算表达式（1）如下，式中：φreduced和creduced为折减后的内摩擦角及粘聚力。

2 工程实例分析

2.1 工程概况

珠三角某市地下综合管廊新建工程，拟建场地属三角洲堆积（河、海）地貌。管廊总长度约7 km，基坑剖面如图1所示：开挖深度约为6.0 m，宽度约10.85 m，采用PU400×170拉森钢板桩支护，竖向设两道ϕ609×16钢管支撑，坑内采用三轴搅拌桩进行满堂加固。

图1 综合管廊基坑剖面图

2.2 有限元模型的建立

该综合管廊基坑工程是一种狭长形基坑，采用Plaxis 2D建立二维有限元模型，围护桩采用板单元模拟，支撑采用锚杆单元模拟，土体采用高精度三角形平面应变单元模拟，有限元模型如图2所示。

图2 综合管廊基坑有限元模型

基于强度折减法的安全系数获得方法为：土层强度折减系数RF由初始值1.0以0.01的步长逐步增加，直至计算无法收敛，选取坑底变形量最大的节点建立隆起量与强度折减系数的关系曲线，曲线的转折点对应的强度折减系数即为安全系数。基坑开挖步骤模拟如表1所列。

表1 基坑开挖步骤模拟一览表

2.3 算例验证

以该工程深厚淤泥区的试验段基坑为算例，根据勘察报告，②层淤泥的基本地质参数见表2所列，坑底加固土的有效内摩擦角和有效粘聚力参数参考文献[16]，取值分别为50 kPa和40°。同时建立有限元模型与现场监测结果相互验证，以获得适合该工程的其他物理力学参数。

表2 土层基本物理力学参数表

根据本文1.1节得出的软土HS模型各刚度参数间的近似关系，结合试验段围护桩水平位移监测数据进行反分析，得出②层淤泥刚度参数关系为：坑底加固土的无侧限抗压强度为qu=0.5 MPa，各刚度参数关系按常规取值。该算例HS模型参数见表3所列。基坑开挖至坑底后的支护桩水平位移与监测数据对比结果如图3所示。

图3 数值模拟曲线与监测结果对比图

表3 HS模型参数表

从图3可以看出数值模拟的曲线能够与监测数据较好地吻合，可以较为真实地反映该基坑围护桩的变形趋势，最大深层水平位移发生在地面下4.5 m处，有限元计算值为14.38 mm，实测值为15.1 mm，说明HS本构模型及相应的物理力学参数选取有一定的合理性，符合工程计算要求。

2.4 基坑破坏形态

数值分析模拟的基坑破坏时形态如图4（a）所示：深厚淤泥下狭长基坑破坏表现为坑底土体发生较大的隆起变形，围护桩底端向坑内移动，同时坑底土体沿围护桩底面滑动时不再表现为完整的圆弧，而是呈现出“ω”形的破裂面。强度折减法计算得出的基坑安全系数为2.24，较规范圆弧滑动法得出的抗隆起安全系数1.98提高了13%，同时试验段的基坑监测数据表明坑底最大隆起量仅为5.48 mm，远未达到隆起破坏的临界点，说明对侧围护桩及土体对狭长基坑的破坏形态和安全系数有较显著的影响，有限元计算方法更能反映该基坑工程的真实情况。

图4 基坑破坏形态图示

3 坑底抗隆起稳定性的影响因素分析

3.1 基坑深度和宽度对坑底抗隆起稳定性的影响

现对各深度、宽度下的基坑坑底抗隆起稳定性进行了比较分析：基坑深度h由初始值6 m以1 m步长增加为9 m，各深度下基坑宽度w为6～24 m，其他参数保持不变，得出安全系数随基坑深度和宽度的关系曲线如图5所示。以基坑深度h=6 m为例，不同基坑宽度下的坑底隆起破坏形态如图6所示。

图5 不同开挖深度下安全系数与基坑宽度的关系曲线

图6 不同宽度基坑的坑底隆起破坏形态图示

图5分析结果表明：同样深度下的基坑，其坑底抗隆起稳定性安全系数随着基坑宽度的增加而逐渐降低；基坑越狭长（宽深比越小），基坑稳定性越高，同时坑底抗隆起安全系数随深度的增加而降低得越快。从图6有限元分析的基坑破坏形态可以看出，随着基坑宽度的增加，坑底滑动面叠加的部分逐渐减少，直至相互分离而呈现为两侧相互独立圆弧滑动面，坑底抗隆起安全系数也基本不再随基坑宽度的增加而变化。

3.2 坑底加固对坑底抗隆起稳定性的影响

深厚淤泥类软土地区基坑设计时，常采用坑内土体加固的措施，以提高基坑开挖期间被动区土体的抗力。为研究坑底加固土体厚度对基坑稳定性的作用，现保持基坑宽度为6 m，将基坑开挖深度由初始值h=6 m以1 m步长增加为9 m，各开挖深度下坑底加固厚度由2 m以1m步长增加至9 m，得出安全系数随基坑深度和坑底加固厚度的关系曲线如图7所示。