1. 引言
随着中国城市化程度不断加快,城市基础设施建设远远落后于城市化进程的矛盾日益显露,越来越多的城市需要进行地下空间的建设,因此基坑和隧道工程大量涌现,并且开挖深度逐渐增大。由于深基坑开挖施工的地质条件和环境条件日益复杂,呈现出很强的复杂性和独特性,工程事故率和损失也越来越大。
由于台州长期受海侵和海退的影响,沿海地区广泛分布有深厚的海陆交互相沉积软黏土。该类软黏土抗剪强度低、承载力低、灵敏度大、空隙比大、极易被扰动,且被扰动后土体性质急剧降低难以恢复,易产生触变、流变,易引起周边土地不均匀变形沉降,这样的地质条件给工程设计和施工提出了严峻的挑战 [1]。付宪章 [2]、古存军 [3] 等以深厚软土地区施工深基坑作为具体案例,分析深基坑支护设计要点和施工过程中的关键控制点,控制地表不均匀沉降;吴波 [4] 等通过现场监测数据,研究具体基坑在施工过程中,长期的地表变形规律;奚家米 [5] 等以上海虹桥SOHO深基坑为研究对象,运用数值计算软件对比研究了不同支护形式对基坑变形的影响;马永峰 [6] 等以杭州某软土深基坑工程为例,分析基坑开挖引起土体变形及对临近地铁隧道的影响。
本文结合台州市中心站S2线车站东端头井深开挖施工技术,利用Midas GTS NX有限元分析软件对车站端头井深开挖过程进行数值模拟计算与分析,为工程建设提供技术依据和保障,达到保证安全、优化设计、指导施工的目的,同时为台州以后的城市地下空间基坑和隧道设计与施工提供理论和技术经验。
2. 工程背景
2.1. 工程概况
台州中心站为台州市域铁路S2线与S1线的换乘站,位于椒江区葭芷街道下北山村市府大道和东官河南侧,采用明挖法施工。
S2线台州中心站为地下三层车站,车站主体深度约31.41 m,支护形式为地下连续墙+内支撑,地下连续墙墙厚1.2 m。长度约为32.8 m,宽度约为31.9 m,处围护结构地连墙最深达71 m,端头井设计尺寸大、开挖深度深。工程范围内为冲海积平原区大范围分布的软土层,主要为淤泥,其具有易触变性、高压缩性、强度低等特性,基坑施工变形控制难度大,安全风险高。
2.2. 场地工程地质条件
根据地质勘察报告,本工程场地属冲海积平原区,地形平坦开阔,深厚层软土为本工程的主要特殊岩土,软土地基的强度低、稳定性差和不均匀沉降及变形大等,是主要工程地质问题。勘探深度范围内的土层之上而下分为10个主要土层,特殊岩土主要包括填土、淤泥和黏土。土体主要物理力学参数如表1所示。
2.3. 端头井施工方案
针对本工程地质软土的特点,根据“时空效应”理论,对基坑端头井围护和开挖过程中的时空效应进行认真分析,严格控制基坑变形。
S2线东端头井围护结构采用地下连续墙 + 内支撑的支护型式。端头井地连墙围护结构封闭后,开始施作冠梁及第一道混凝土支撑,当达到设计强度后,开始基坑开挖。基坑开挖严格遵守“纵向分段,水平分层,随挖随撑,严禁超挖”的原则进行施工。每个工作面内均采用流水作业方式,做到开挖、支撑、监控量测等相关工序紧密配合。
S2线车站为地下三层结构,设六道支撑,其中端头井第一、二、四道为砼支撑,其余均为钢支撑。以各道支撑为界,竖向相应分为六层放坡开挖。基坑中部设有格构柱,格构柱两侧各设置1台挖掘机,两侧土体同时开挖,高差控制在2 m以内。开挖时,每一个放坡段长度一般不超过12 m,开挖时间和钢支撑的安装时间控制在16小时和8小时以内。开挖中及时测定支撑安装点,监控支撑端部中心位置控制在容许限值内。
3. 端头井施工数值模拟
3.1. 有限元模拟及模拟工况
基坑工程在开挖施工过程中所处的状态比较复杂,面临的问题及困难也很多。为保证工程本身及周围环境的安全,需要充分了解施工过程中结构及土体的受力变形情况。
结合台州中心站S2东端头井深基坑开挖实际情况,建立计算模型,如图1所示。为了减小有限元模型中边界约束条件对计算结果产生的不利影响,计算模型的边界范围在各个方向上均大于三倍的端头井深度。具体计算时,计算区域在水平向上的长度为210 m,宽度为147 m,在竖直方向上,取竖直长度为120 m。整个计算模型划分75,467个单元。
计算模型中,根据场地地层情况,将土体划分为10层,土体采用三维实体单元模拟,定义土体为弹塑性体,采用修正Mohr-Colomb本构,相关参数取值情况如表1所示;地连墙采用三维结构单元模拟,厚度为1.2 m,采用C35混凝土材料,桩基和支撑采用一维结构单元模拟,均为混凝土材料,地连墙、桩基和支撑结构均定义为弹性体,模型顶面为自由面,无约束;底面采用固定端约束;4个侧面均只约束法向,其余方向自由无约束。
在模型仿真计算中,计算模型的网格划分和计算结果的收敛性直接相关,在对基坑内部土体和坑周土体网格划分时,对坑内土体进行加密划分,控制网格尺寸,保证计算精度和模型收敛性,坑外至模型边界处尺寸较大,以提升计算效率。
3.2. 工况模拟
1) 初始位移清零。在进行模型仿真计算时,需要考虑土体自身的固结沉降以及周边建构筑物引起的地表沉降,因此,需要在基坑施工的步骤前进行模型的初始位移清零。
2) 地下车站深基坑端头井施加第1道支撑后,开挖至第2道支撑底;
3) 地下车站深基坑端头井施加第2道支撑后,开挖至第3道支撑底;
4) 地下车站深基坑端头井施加第3道支撑后,开挖至第4道支撑底;
5) 地下车站深基坑端头井施加第4道支撑后,开挖至第5道支撑底;
6) 地下车站深基坑端头井施加第5道支撑后,开挖至第6道支撑底;
7) 地下车站深基坑端头井施加第6道支撑后,开挖至端头井坑底。
4. 计算结果分析
本文重点研究开挖过程对端头井周边土体沉降、地连墙侧向位移和内支撑轴力分布的影响。
4.1. 实测计算结果对比
在端头井开挖前,于北侧地连墙距东墙面5 m左右处设置测斜管,由于现场施工条件制约,实测采集到开挖第一层和开挖第五层的地连墙实际测斜数据,和模型计算数据对比,整理如图2和图3所示侧向位移图。
![](//html.hanspub.org/file/3-2751305x8_hanspub.png?20220217080704770)
Figure 2. Lateral displacement of diaphragm wall in the first step of excavation
图2. 第1步开挖地连墙侧向位移
![](//html.hanspub.org/file/3-2751305x9_hanspub.png?20220217080704770)
Figure 3. Lateral displacement of diaphragm wall in the fifth step of excavation
图3. 第5步开挖地连墙侧向位移
由图2和图3可知,在开挖第一层时,实际测量值与理论计算值的变化趋势一致,即侧向位移都是随着深度的增加,先增大后减小,侧向位移最大值都出现在地下9 m左右,该处实测最大位移值为1.83 mm,模型计算最大位移值为3.07 mm;在开挖第五层时,实际测量值与理论计算值的变化趋势同样一致,也为侧向位移都是随着深度的增加,先增大后减小,侧向位移最大值都出现在地下21.3 m左右,该处实测最大位移值为12.53 mm,计算最大位移值为9.75 mm。根据测量值和理论值变化趋势和转折点位置,考虑有可能出现的测量误差,计算模型与工程实际结果匹配,计算结果参考性强。
4.2. 地表沉降结果分析
提取图1中A-A’断面的计算结果,绘制地表沉降图,如图4所示。由计算结果可知,周边土体沉降随着开挖层数的增加不断增大;根据图中每个开挖步骤对应的最大沉降量可知,增速随着开挖深度的增加会适当放缓。且开挖对地表沉降的影响范围约60 m左右,约等于两倍的端头井宽度(31.9 m),当距离端头井边缘超过60 m处,地表土地沉降不足1 mm,可以忽略不计;由图4可知,周边土体在开挖的过程中,地表土体沉降呈现的沉降槽形状和位置不变,即变形量大小关系随着距离端头井边缘的远近变化规律一致,最大的位置几乎不变,均为距离基坑边15 m左右,在第六步开挖施工完成后,该区域沉降值最大。
在端头井开挖过程中,土层的应力会发生重分布,会造成土体的变形沉降。地表沉降为槽型,原因是距离端头井较近位置处有围护结构加固,地表沉降变形受基坑开挖影响较小,距离端头井边缘稍远处,土体挤压作用明显,地表沉降变形大,远离基坑边缘的位置,开挖对土层的应力分布逐渐变小,因此地表沉降逐渐减小。
4.3. 地连墙水平位移结果分析
根据计算模型结果整理得到,基坑开挖至坑底时,地连墙的侧向变形云图如图5所示。地连墙各开挖步骤对应的最大侧向位移如表2所示。
![](Images/Table_Tmp.jpg)
Table 2. Maximum lateral displacement and maximum displacement position of each excavation step
表2. 各开挖步侧向位移最大值及最大位移位置表
由图5和表2分析可知:随着开挖的不断深入,地连墙最大侧向位移量也在不断变大,地连墙最大位移发生点也在不断下移,最大位移点分布在坑底附近,为2.12 cm。具体来说,开挖初期最大位移点出现在开挖面以下,随着开挖的进行,最大位移点最终处于开挖面以上位置,均在开挖面附近。
结合各个开挖步骤,在开挖及布置内支撑以后,发生最大水平方向位移处在不断下移,说明开挖产生应力重分布,会伴随着对地连墙体产生力的作用,并且在开挖面处作用效果最明显。在内支撑施工完成以后,该位置墙体位移变小,内支撑对地连墙荷载的传递作用明显。
4.4. 地连墙水平位移结果分析
根据计算模型整理得到基坑分步开挖过程中,各道支撑最大轴力随开挖进行计算值如表3所示。
![](Images/Table_Tmp.jpg)
Table 3. Calculation value of support axial force of each layer
表3. 各层支撑轴力计算值
由表3可知,第1道支撑轴力在开挖到第二层时达到最大,随着开挖的进行,支撑轴力逐渐减小;第2道、第3道、第4道、第5道和第6道支撑轴力均随着开挖的进行,轴力值逐渐增大,开挖到坑底时达到最大值;其中开挖至坑底时,第4道支撑轴力大于其他支撑结构的支撑轴力,最大值达到6792 kN。
通过表3可以看到,支撑轴力会随着开挖进行,迅速对地连墙等围护结构产生力的作用,以控制变形。端头井在开挖过程中,开挖不断向下,内支撑受到的压力会快速增大。但是随着后一道支撑的施加,轴力得到分担,轴力的变化趋势产生拐点,随后随着支撑开挖,变化趋于平缓。
5. 结论
本文采用摩尔–库伦模型对深大基坑端头井过程进行模拟计算,研究深厚软土地区基坑开挖对地表沉降的影响、开挖过程中地连墙侧向变形情况以及基坑横撑轴力变化情况。
数值计算结果表明,在深厚软土地层基坑开挖过程中,基坑周围土层会出现沉降现象,并且随着开挖的进行,沉降量会不断变大。沉降最大值出现在距离基坑边缘15 m处,即0.5个基坑宽度处,沉降的影响范围约60 m左右,为两倍的基坑宽度。
在基坑开挖的过程中,作为围护结构的地连墙会发生水平位移,位移量随着开挖的进行不断变大,在开挖至坑底时达到最大值,位移最大值出现的位置和开挖面位置相近。
基坑开挖过程中,内支撑受到压力迅速增加,随着后一道支撑完成施工,上一道支撑轴力会得到分担;开挖到第二层时,第一道支撑轴力达到最大值,开挖到坑底时,其他支撑轴力达到最大值。