1. 引言
软土地区超深地连墙开挖容易引发较大的地层变形，并且与地层侧向土压力的变化特征有较大关系，基坑开挖后会产生较大的卸荷效应，导致地应力重分布。当软土地区基坑工程开挖深度大时，开挖后还会产生较大的地层位移，这又会反作用于土压力。因此，要更好地解决超深地连墙施工中的变形和破坏问题就需要弄清超深地连墙施工全过程的土体变形特征。本研究依托工程的土体具有“含水量高、孔隙比大、强度低及压缩性高”的特点，土体本身也会随着时间产生蠕变，这对土压力的大小也会产生影响，具有明显的时间效应 [1] [2] [3] [4] [5]。
目前分析地下连续墙施工引起的地层响应的计算方法主要有两种：极限平衡方法与数值方法。尽管极限平衡分析方法简便，机理清晰，但其较难反映周围土体的应力变化和变形性状，因此不能确定地层侧向位移和地面沉降大小 [6] [7] [8]。而数值方法却可以弥补极限平衡方法的不足，基于连续体介质力学数值方法能够模拟软粘土的应力与变形性态，刻画地下连续墙施工引起的土体变形和受力的动态变化特征。另外，分析地下连续墙施工引起的地层响应评价方法研究主要集中在对以往工作从折减系数(朗肯理论)、破坏形式和搜索最危险滑动面方法(库伦理论)这几个方面的改进和修正 [9] [10] [11]。基于库伦理论的评价方法在工程中应用较广泛，但以往的研究都是假定成槽周边土体为均质，对层状地层结构中地下连续墙成槽开挖稳定性评价方法的研究较少 [12] [13] [14] [15]。
综上可知，超深地下连续墙成槽施工过程对周围土体原有的应力和变形会产生更大影响，对土体扰动程度以及土体的动态变化特征进行系统的研究就显得很有意义。如果土体扰动较小，超深地下连续墙成槽施工过程效应将并不明显，后续的设计中也不需专门考虑其作用；反之，如果影响作用较大，且考虑影响下的变形计算结果较之原设计计算影响作用要大，则更应引起设计和施工单位的高度重视。本文将依托上海市域线机场联络线梅富路工作井基坑工程采用三维数值模拟方法开展超深基坑施工全过程的变形和受力分析，研究成果可为工程实践提供参考。
2. 工程概况
本研究依托工程为上海市域线机场联络线(西段) JCXSG-3标梅富路工作井，该工作井西临梅富路，南侧为春申塘，坑底深度约为31米。工作井基坑工程范围现状场地位于梅南路上；东侧为梅强路，南侧临近220 kV高压电线塔；西侧为梅富路市政泵站，1~2层混凝土结构，距离基坑最近处约49 m，北侧为空地。梅富路工作井采用1200 mm厚度的地下连续墙围护加内支撑，接头为钢性接头，地墙两侧加固采用
850三轴水泥土搅拌桩加固。工作井位置及周边环境平面图如图1所示。梅富路工作井所处场地位于上海市中部，且地势较为平坦，钻探标高一般在4.48~4.91 m之间。超深地连墙穿越的土层主要为填土、淤泥质粘土、粉质粘土等软土地层。
(a) 梅富路工作井所处位置 
(b) 梅富路工作井周边环境平面图
Figure 1. Location of Meifu Road working shaft
图1. 梅富路工作井地理位置
本工程基坑开挖深度达到36.5 m，地连墙插入深度达到66 m，施工过程中存在许多安全隐患，所以本研究将利用ABAQUS有限元软件对梅富路工作井的开挖施工全过程进行三维数值模拟分析，模拟中考虑了坑外土体性质、坑内支撑体系以及开挖施工步骤等，通过模拟结果探究了超深基坑开挖导致的土体及围护结构地连墙的变形情况，可为工程提供参考。
3. 三维数值计算模型和计算参数
3.1. 计算模型
本模拟选取整个梅富路工作井及周围土体作为研究对象，三维数值模型如图2所示。考虑到边界效应的影响，地层模型尺寸确定为400 m × 250 m × 100 m，并将其沿深度方向根据地勘报告划分为8个土层，采用Mohr-Coulomb本构模型；按照对应土层参数，采用三维实体单元C3D8进行土体建模。为计算方便，地下连续墙建模将各幅地连墙建立为一个整体，且由于地连墙沿幅宽方向尺寸远大于厚度方向尺寸，所以选用壳单元S4R进行模拟，平面尺寸为两个工作井25 m × 25 m，连廊15 m × 90 m，地连墙高度66 m，厚度1.2 m。支撑体系中包括支撑横梁、连系梁和格构柱，均采用梁单元B31进行模拟，其中横梁截面尺寸为第一、三、五道钢筋混凝土支撑1000 m × 1000 m，第七道钢筋混凝土支撑1200 m × 1000 m，第九道钢筋混凝土支撑1400 m × 1000 m，其余均为钢管支撑
800 mm，t = 200 mm，连系梁和格构柱截面尺寸均为600 m × 600 m。模型中相关约束为墙体和土体绑定，支撑横梁和墙体绑定，连系梁、格构柱和支撑横梁之间绑定，格构柱埋入土体部分与土体施加嵌入式约束。边界条件为模型四周设置法向约束，底部设置固定约束，地表自由。整个模型一共305,720个网格结点，288,685个网格单元。
(a) 地层模型 
(b) 地连墙模型 
(c) 支撑体系模型
Figure 2. 3D numerical simulation model of Meifu road working shaft
图2. 梅富路工作井三维数值模型
3.2. 计算参数
三维数值模拟分析中采用的土体本构模型是Mohr-Coulomb模型，依据地勘报告提供的参数值，综合分析后将土体分成8层进行建模，各土层及对应的建模参数如表1所示。其中，建模时土体的弹性模量按照地勘给出的压缩模量的2~5倍进行取值。
Table 1. Mechanical parameters of ground strata
表1. 地层力学参数
4. 数值模拟步骤
本模拟根据工作井基坑实际开挖步骤进行数值模拟，计算步骤如下：
1) 第一步为地层的地应力平衡；
2) 第二步为超深地连墙及格构柱的施工，在本步骤重激活事先杀死的地连墙和格构柱模型及其约束来模拟地连墙和格构柱的施工；
3) 后续开始基坑开挖以及支撑体系的施工模拟按照工序分十步杀死基坑内土体，同时在每一步重激活该层支撑体系完成基坑开挖和支护。分析步骤如表2。
Table 2. Procedures of 3D numerical simulation
表2. 三维数值模拟步骤
5. 结果分析
5.1. 连廊方向地连墙侧向位移
通过三维数值分析可得到基坑开挖完成后工作井连廊方向地连墙侧向位移随开挖施工的变化典型云图如图3所示。由图3可以看出，超深地连墙沿深度方向变形呈两边小中间大的鼓形，并且中间最大位移发生在开挖面处。随着开挖的进行，最大位移位置在逐渐下移，模拟的变形规律与现场监测结果吻合。总体上，变形分布在沿连廊方向呈连廊中间大，两边工作井处小的鼓形。
Figure 3. Displacement contour of diaphragm wall along connection corridor (unit: m)
图3. 连廊方向地连墙侧向位移云图(单位：m)
连廊中间部分侧向位移变化与现场监测数据对比如图4所示。可以看出地连墙侧向的位移变化规律与实际监测结果基本一致，但也存在一定的差异，主要体现在地连墙下半部分发生朝向基坑外侧的反弯变形，这是由于模拟中选取的土体本构为Mohr-Coulomb模型所致，该本构模型在岩土开挖模拟中会出现较明显的地表膨胀隆起的现象，导致基坑底部边缘土体发生隆起；此外，模拟结果中地连墙侧向位移最大值有一定偏差，这可能与计算参数和实际存在差异所致，但整体变形规律一致。
5.2. 超深地连墙周围土体变形特征
随着深基坑的开挖，基坑外土体由于卸载作用产生应力释放发生变形。在刚开始开挖基坑时，由于地应力的释放，基坑两侧浅层土体发生向基坑内侧的变形，最大影响范围达到50 m。由于此时支撑体系还未完工，所以地连墙会随着土体变形发生向基坑内侧的倾覆，导致地连墙墙趾有向基坑外侧移动的趋势，引起周围土体发生与地表土体相反的变形。但由于支撑体系的刚度大于基坑内土体的刚度，所以随着开挖的进行以及支撑体系的作用，基坑开挖对地表土体的影响范围也越来越小，基坑外土体的变形主要集中在开挖面两侧，并且在开挖面与基坑壁交角处发生明显的地表隆起情况，要特别注意超深基坑开挖过程中地连墙的反弯变形和反弯点的动态变化，这与水土压力随着深度的变化关系较大。整个开挖过程中土体全过程变形如图5所示。
Figure 4. Comparison analysis of lateral displacement of diaphragm wall along connection corridor
图4. 连廊中部地连墙侧向位移对比分析
6. 结论
本文以上海市域线机场联络线梅富路工作井超深地连墙为研究对象，通过三维数值模拟和现场实测研究，获得了软土地区超深地连墙施工全过程的土体变形特征，探究了超深地连墙侧向变形模式，并将计算结果与数值模拟结果进行对比分析，得到如下结论：
1) 通过ABAQUS有限元软件对深基坑开挖过程进行了三维数值模拟研究，获得了超深基坑开挖全过程对超深地连墙变形的主要影响因素，其最大影响范围是50米左右；
2) 数值分析中全面考虑工作井施工、支撑体系以及周围土体性质，模拟了超深地连墙施工全过程的动态演化特征以及基坑外土体的变形过程，并与现场监测数据进行对比分析，发现软土地层中超深基坑中间最大位移是发生在开挖面处，且随着开挖深度不断下移；在开挖面与基坑壁交角处会发生明显的地表隆起情况，因此超深地连墙施工中要重点重视开挖面的土体变形状态。
致谢
本项目得到了上海市科委重大科技攻关项目(19DZ1201004)的资助，在此一并致谢。
NOTES
*通讯作者。