考虑时空效应的杭州软黏土超深基坑地表沉降分析方法
Analysis Method for Ground Settlement Induced by Ultra-Deep Excavation in Hangzhou Soft Clay Considering Time-Space Effect
通讯作者: 应宏伟,教授,博士生导师;E-mail:ice898@zju.edu.cn.
责任编辑: 王一凡
收稿日期: 2023-05-18 修回日期: 2023-08-24 接受日期: 2024-05-21
基金资助: |
|
Received: 2023-05-18 Revised: 2023-08-24 Accepted: 2024-05-21
作者简介 About authors
熊一帆(1999—),硕士生,从事地下工程方面的研究.
利用PLAXIS 3D软件和软土蠕变模型,建立杭州中心项目超深基坑群B2基坑的三维数值模型,深入分析了时空因素对坑外地表沉降的影响.引入互补误差函数和三折线模型,利用数值结果对既有互补误差函数进行修正,提出了考虑时空效应的地表沉降快速分析方法.结果表明:软黏土蠕变既诱发围护墙的附加侧移从而引起土体沉降,又诱发坑外土体产生不依赖于围护墙侧移的沉降;忽略软土蠕变对预测超深基坑坑外地表沉降的影响不弱于其对墙体侧移的影响;基坑深度和开挖速率相同时,开挖面积直接决定施工时间的长短,从而影响软黏土蠕变诱发的地表沉降.
关键词:
The PLAXIS 3D software and the soft soil creep model were used to establish a 3D numerical model of the B2 excavation within an ultra-deep excavation group, allowing for analysis of the effects of time-space factors on ground settlements. Then, the complementary error function and trilinear model were introduced. After that, the numerical results were used to correct the complementary error function. Finally, a rapid approach was proposed for predicting ground settlements while accounting for the temporal and spatial influences. The results indicate that soft soil creep not only causes additional wall deflections resulting in soil settlements but also induces soil settlements independent of these deflections. The impact of ignoring soft soil creep on prediction of ground settlement is not weaker than that of wall deflections. The excavation area determines the construction duration when the excavation depth and rate are identical, impacting ground settlements induced by soft soil creep.
Keywords:
本文引用格式
熊一帆, 应宏伟, 张金红, 程康, 李冰河.
XIONG Yifan, YING Hongwei, ZHANG Jinhong, CHENG Kang, LI Binghe.
当前,研究基坑开挖诱发坑外地表沉降的方法大多分为以下几种.一是现场实测法.王卫东等[4]整理统计了35个工程的实测数据,给出了上海软土地区深基坑地表最大沉降δvm与围护墙最大侧移δhm比值的平均值为0.84,该数值与江晓峰等[5]归纳总结的平均值0.83极为接近.应宏伟等[6]基于杭州某深厚软黏土基坑的实测数据进行分析,发现时空特性是影响坑外地表沉降的重要因素,且δvm≈0.6δhm.二是理论解析法.González等[7]利用虚拟镜像技术和弹性力学理论推导了圆孔收缩诱发土体位移的解析公式;Zheng等[8]利用González等[7]所推公式和面积等效原理,进一步推导了围护墙变位诱发坑外地表沉降的公式.Qian等[9]基于平面应变下的位移-位移弹性边值问题,推导了刚性围护墙水平变位诱发坑外地表沉降的显式解;以此为基础,许多学者进一步推导了柔性围护墙变位下坑外地表沉降的理论解答[10⇓-12].当围护墙的变位模式为抛物线时,所有理论解答给出的关系均为δvm/δhm≈0.4.三是数值模拟法.Mu等[13]利用实测数据对某基坑进行反分析建模,总结数值模型中土体位移的衰减规律,提出了考虑小应变特性的土体三维位移场经验公式.Goh等[14]通过大量有限元分析确定了影响地表沉降的关键参数,并开发了人工神经网络模型来估计地表最大沉降.Ying等[15]利用既有工程资料建立了杭州地铁基坑的标准化模型,基于标准化模型进行三维数值模拟,提出了一种针对杭州地区地铁基坑坑外地表沉降的估算方法.
基坑开挖下坑外软土的变形和强度会随时间发生明显变化,实测研究成果[4⇓-6]虽能在一定程度上反映软土的蠕变特性,却无法详尽考虑不同时效或非时效因素对坑外地表沉降的影响;解析方面,既有理论方法没有考虑软土的蠕变效应,且大多基于平面应变假定,无法考虑坑外地表的空间效应;利用数值模拟方法的研究成果[13⇓-15]中,HS-Small模型[13,15]由于能够考虑土体的小应变特性常为学者们所使用,然而其无法考虑与时间相关的软黏土蠕变.另一方面,坑外地表沉降的性状受到地层条件、基坑尺寸和围护结构型式等不同因素的影响,一般很难对坑外地表沉降进行准确的评估.然而,周边环境复杂的软土顺作法超深基坑一般在围护结构体系设计、土方开挖及施工上存在很多相似性,对某地区特定地层条件下典型的超深基坑案例进行建模,以期给出该地区超深基坑开挖引发地表沉降的可行分析方法.
综上所述,本文利用PLAXIS 3D软件和软土蠕变模型,建立杭州某超深基坑的三维数值模型,与实测数据的对比验证了模型的合理性;基于参考算例展开参数讨论,深入分析了时空因素对围护墙侧移、地表沉降及二者关系的影响;引入互补误差函数和三折线模型,利用数值分析结果对既有互补误差函数进行修正,提出了考虑时空效应的地表沉降快速分析方法,并应用于其他超深基坑工程验证该方法的适用性.
1 模型建立与验证
1.1 工程概况
杭州中心项目超深顺作基坑采用分坑施工方法,Cheng等[16]分析了B2分坑开挖的变形性状,平剖面如图1所示.图中:z为纵坐标;q为地连墙测斜孔;d为地表沉降观测点.如图1(a)所示,主要地层自上而下分别为:①填土、②软黏土、③粉质黏土、④圆砾、⑤强风化、⑥中风化砂岩;基坑开挖深度为30.2 m,平面尺寸约98 m×92 m;围护结构采用厚1.2 m、深50 m的地连墙结合6道混凝土满堂对撑和被动区裙边加抽条加固.如图1(b)所示为地连墙测斜以及坑外地表沉降监测点,自2019年4月1日第1道支撑形成开始,至2019年11月16日地下室基础底板形成强度,基坑开挖施工可分为如表1所示的6个主要工况.
图1
表1 基坑施工进度表
Tab.1
工况 | 施工内容 | 时长/d |
---|---|---|
1 | 开挖至-8.3 m,第2道支撑形成强度 | 30 |
2 | 开挖至-13.8 m,第3道支撑形成强度 | 37 |
3 | 开挖至-17.9 m,第4道支撑形成强度 | 43 |
4 | 开挖至-22.1 m,第5道支撑形成强度 | 32 |
5 | 开挖至-26.3 m,第6道支撑形成强度 | 31 |
6 | 开挖至-30.2 m,底板形成强度 | 57 |
1.2 三维数值计算
采用PLAXIS 3D有限元软件,考虑对称性,将前述B2分坑的1/4作为三维数值分析对象.为降低模型边界对计算精确性的影响,地连墙至水平边界的距离取5倍开挖深度,竖向边界取至地连墙墙底以下约20 m,模型尺寸为200 m×197 m×70 m.至于计算模型的边界条件则为:下边界位移约束,上边界位移自由,水平边界水平方向位移约束,竖直方向位移自由;设置每个开挖步中临空面为孔压为0的边界面.
计算分析时采用的土体模型及其参数如表2所示.表中:l为层厚;γ为土体重度;c'为土体有效黏聚力;φ'为土体有效内摩擦角;λ*为修正压缩指数;κ*为修正膨胀指数;μ*为修正蠕变指数;
表2 土体计算参数
Tab.2
本构模型 | 土层 | l/m | γ/(kN·m-3) | c'/kPa | φ'/(°) | λ* | κ* | μ* | |||
---|---|---|---|---|---|---|---|---|---|---|---|
SSC | ②软黏土 | 20.4 | 17.8 | 5 | 25 | 0.06 | 0.006 | 0.002 5 | |||
γ0.7 | m | ||||||||||
HS-Small | ①填土 | 6.5 | 17.0 | 5 | 28 | 4.5 | 3.8 | 22.7 | 79.4 | 0.000 2 | 0.80 |
③粉质黏土 | 13.5 | 19.3 | 25 | 30 | 6.5 | 5.4 | 32.4 | 113.4 | 0.000 2 | 0.80 | |
加固土 | 23.7 | 19.3 | 60 | 25 | 15.0 | 15.0 | 85.0 | 170.0 | 0.002 0 | 0.12 | |
E/MPa | υ | ||||||||||
MC | ④圆砾 | 2.2 | 19.0 | 1 | 18 | 36 | 0.30 | ||||
⑤强风化砂岩 | 6.9 | 21.0 | 30 | 20 | 90 | 0.25 | |||||
⑥中风化砂岩 | 20.5 | 22.0 | 50 | 30 | 3 600 | 0.25 |
图2
1.3 模型验证
图3
图3
基坑变形计算值与实测值对比
Fig.3
Comparison of excavation deformations between calculated and measured value
2 参数分析
将所建数值模型作为参考算例,以工况6作为对比工况,进一步探讨不同时空因素对基坑变形,尤其是坑外地表沉降的影响.值得说明的是,墙体、支撑刚度对超深基坑来说一般都足够大,且杭州主城区典型软黏土的刚度也基本相同,因此本文不对以上参数展开分析,仅聚焦影响软黏土超深基坑坑外地表沉降的5个关键因素.本文利用修正系数对参考算例中δvm和δvm/δhm进行修正,使方法具有普适性,而非仅给出定性的结论.
2.1 开挖速率
软黏土基坑的蠕变变形受开挖时间t的影响,通过改变表1中不同工况的时长,图4给出了基坑开挖速率v=Vs/t对基坑变形的影响,Vs为基坑需开挖的土方体积.图中:δm为基坑最大变形量;vm=δm/t,为最大变形的平均变化速率;vhm和vvm分别为围护墙最大侧移和地表最大沉降的平均变化速率.由图4(a)可知,地表沉降受开挖速率的影响与围护墙侧移相当,考虑到地表最大沉降增量一般小于围护墙最大侧移增量,可以说明软土蠕变诱发坑外土体产生了不依赖于围护墙侧移的沉降;由图4(b)可知,随着开挖速率的提高,基坑变形减小,而最大变形的平均变化速率却有所提高.工程实践中,过大的变形和变形速率都是不被允许的,因此设计施工人员应选择合理的开挖速率;由图4(c)可知,随着开挖速率的提高,δvm/δhm逐渐下降但下降速率趋缓,这是因为开挖速率的提高意味着时效变形的减小,而非时效地表沉降仅与非时效墙体侧移的形态和大小有关,所以δvm/δhm存在极限小值.
图4
从图4可得土方开挖速率对参考算例(v=1 184 m3/d)的δvm和δvm/δhm的修正系数:
式中:ξv、ηv分别为土方开挖速率对参考算例的δvm和δvm/δhm的修正系数;v0为单位开挖速率;R为Pearson相关系数.
2.2 软土层厚度
图5
图5
软土层厚度对基坑变形的影响
Fig.5
Influence of soft soil thickness on excavation deformations
可得软土层厚度对参考算例(ls=20.4 m)的δvm和δvm/δhm的修正系数:
式中:
2.3 软土修正蠕变指数
PLAXIS 3D软件中软土蠕变模型的修正蠕变指数的计算公式[20]如下:
式中:Cα为蠕变指数,可由室内试验得到,若无试验数据时可利用经验公式Cc/Cα≈15~25来确定[20],Cc为压缩指数;e0为土体的初始孔隙比.
图6
图6
软土修正蠕变指数对基坑变形的影响
Fig.6
Influence of modified creep index of soft soil on excavation deformations
从图6可得到软土修正蠕变指数对参考算例(μ*=0.002 5)的δvm和δvm/δhm的修正系数:
式中:
2.4 开挖面积
图7
由图7可得开挖面积对参考算例(S=9 016 m2)的δvm和δvm/δhm的修正系数:
式中:ξS、ηS分别为基坑开挖面积对参考算例的δvm和δvm/δhm的修正系数;S0为单位开挖面积.
2.5 基坑长宽比
在不改变基坑开挖面积的条件下,图8给出了不同长宽比L1/L2对基坑变形的影响,δhm、δvm取自边长为L1的围护墙及其坑外地表.这里仅考虑 1/3<L1/L2<3的情况,是因为实际工程中一般会利用基坑空间效应分块或分仓开挖,避免基坑施工接近平面应变状态.可以看出,随着L1/L2的提高,基坑变形基本呈线性提高,显然这是因为随着L1/L2逐渐增大,矩形基坑边长为L1中部的空间效应减弱,逐渐接近平面应变状态;同时也可发现 δvm/δhm 随着L1/L2的提高逐渐增长但增长速率趋缓.
图8
图8
长宽比对基坑变形的影响
Fig.8
Influence of length-width ratio on excavation deformations
由图8可得基坑长宽比对参考算例(L1/L2=1.065)的δvm和δvm/δhm的修正系数:
式中:
2.6 δvm和δvm/δhm拟合公式的建立
基于数值模型的参考算例中不同开挖深度H下的坑外地表沉降最大值,图9给出了基坑开挖深度与坑外地表沉降最大值的关系.
图9
拟合公式如下:
式中:H0为单位开挖深度.结合前文考虑时空因素的各项修正系数即可得到杭州软黏土地区任意超深基坑地表沉降的最大值:
如果已有确定墙体最大水平位移δhm的经验关系,也可按下式确定坑外地表最大沉降δvm:
3 地表沉降的快速分析方法
3.1 平行于围护墙方向
式中:A为原点至拐点的距离,拐点指平行于围护墙剖面上地表沉降为沉降最大值1/2的点.
图10
图10
平行于围护墙方向的地表沉降示意图
Fig.10
Schematic diagram of ground settlements parallel to retaining wall
然而,参数A为由实测数据拟合得到的经验值,并不适用于时间效应显著的杭州软黏土地区,需要利用本文数值结果对A重新拟合,如图11所示.拟合公式如下:
图11
利用式(15)、(17),结合y=49.0 m的实测沉降,图12给出了平行于围护墙方向的计算值与本文参考算例的实测数据对比.可见计算值与实测值较为吻合,验证了拟合参数的合理性.
图12
图12
地表沉降计算值与实测值对比(平行于围护墙方向)
Fig.12
Comparison of ground settlements between calculated and measured value (parallel to retaining wall)
3.2 垂直于围护墙方向
引入《基坑工程技术规范》[23]中的三折线模型对垂直于围护墙方向的地表沉降进行描述,三折线模型定义了坑外4点的沉降与沉降最大值的关系,4点分别为(0, 0.2δvm)、(H/2, δvm)、(2H, 0.1δvm)、(4H, 0).
利用式(12)、(15)和(17),图13给出了垂直于围护墙方向计算值与参考算例的实测值对比.可见,除距围护墙32.0 m处受道路交通荷载扰动的监测点以外,其他各点计算值与实测值吻合较好.
图13
图13
地表沉降计算值与实测值对比(垂直于围护墙方向)
Fig.13
Comparison of ground settlements between calculated and measured value (perpendicular to retaining wall)
4 应用与验证
上述超深基坑地表沉降快速分析方法是基于杭州中心项目超深基坑群中最早施工的B2基坑实测数据和三维数值分析建立的,为检验其适用性,将该方法应用于该超深基坑群中第2期施工的B1基坑.
4.1 工程概况
杭州中心B1基坑深度为30.2 m,开挖面积为 5 219 m2,平面非严格矩形,尺寸约70.9 m×64.0 m,分析对象位于临近市政道路的短边,L1/L2=0.903;如图14所示,坑外埋设了3×4个地表沉降监测点,垂直于围护墙方向的监测点由近及远分别距围护墙2.0、7.0、17.0、45.5 m,平行于围护墙方向上,d1~d3沉降观测断面分别距较近的角部17.3、29.0、11.2 m;自2020年4月3日第1道支撑形成开始,至2020年8月17日地下室基础底板形成强度,基坑开挖施工总计约137 d,开挖速率约 1 150 m3/d;在土层条件和围护结构设计参数等方面,B1与B2基坑一致.从平面上看,B1基坑规模明显比B2小,同样开挖深度下坑外地表沉降受开挖尺寸的影响更大;从时间因素考虑,二者开挖速率大抵相当,故B1基坑开挖施工的历时明显小于B2基坑,受软黏土蠕变的影响相对更小.
图14
4.2 计算流程
(1) 将v、ls、μ*、S、L1/L2代入式(1)、(3)、(6)、(8)、(10),求得ξv=1.01、
(2) 将ξv、
(3) 将H和L1=64.0 m代入式(17)求得H=22.1,26.3,30.2 m时的A=9.1,7.9,7.0 m.
(4) 将δvm、L1和A代入式(15)可求得不同开挖深度下y=29.0,17.3,11.2 m剖面的δvm.
(5) 将不同深度下y=29.0,17.3,11.2 m剖面的δvm代入三折线模型即可求得不同深度下相应剖面任意点处的δv.
4.3 结果对比
图15给出了本文方法与实测数据的地表沉降对比.可见,计算值与实测值较为相符,可以认为本文方法可以较为准确地预测杭州软黏土地区超深基坑坑外地表沉降.同时还可发现,B1基坑的坑外地表沉降计算值和实测值均明显小于B2基坑,与B1基坑平面尺寸更小、施工历时更短、时空效应更为有利的客观分析是吻合的.
图15
图15
B1基坑地表沉降计算值与实测值对比
Fig.15
Comparison of ground settlements of excavation B1 between calculated and measured value
5 结论
采用软土蠕变模型,建立杭州中心项目超深基坑群B2基坑的三维数值模型;基于参考算例展开参数讨论,深入分析了坑外地表沉降的时空效应,获得了考虑时空因素的坑外地表最大沉降的修正系数;引入互补误差函数和三折线模型,利用数值结果对既有互补误差函数进行修正,提出了考虑时空效应的地表沉降快速分析方法.结果表明:
(1) 软黏土蠕变既诱发围护墙的附加侧移从而引起土体沉降,又诱发坑外土体产生不依赖于围护墙侧移的沉降;忽略软土蠕变对预测超深基坑坑外地表沉降的影响不弱于其对墙体侧移的影响.
(2) 基坑深度和开挖速率相同时,开挖面积直接决定施工时间的长短,从而影响软黏土蠕变诱发的地表沉降.
(3) 本文方法可以较为准确地预测杭州深厚软黏土超深基坑的坑外地表沉降,也为建立其他软黏土地区超深基坑开挖环境影响的快速评估方法提供了思路.
需要说明的是,本文研究的是城市中心区软黏土超深基坑工程的环境影响预测方法,由于此类工程的重要性,其围护墙刚度、插入深度以及支撑刚度一般都足够大,所以提出的预测方法未考虑这些常规因素的影响.
参考文献
紧邻基坑同步施工下坑间隧道的变形特性
[J].
DOI:10.16183/j.cnki.jsjtu.2013.10.009
[本文引用: 1]
采用非线性平面有限元方法研究紧邻基坑同步施工下坑间土体与隧道的变形特性.将该方法运用到上海某基坑工程中的隧道位移预测,实测数据与有限元结果比较吻合,验证了其可行性.通过参数分析,讨论了基坑之间的距离、开挖宽度对坑间土体竖向位移的影响,并由此推算出了使得坑间地表土体平均竖向位移接近零的最优基坑间距与开挖宽度的关系.同时讨论了盾构隧道、明挖暗埋隧道与基坑的距离对于隧道竖向位移、支护结构侧移的影响.研究表明,开挖完成后,坑间盾构隧道的沉降比周围土体的沉降略小,而有围护结构的明挖暗埋隧道则表现为上抬,距离基坑越远,上抬量越小.
Movement of tunnel between two adjacent concurrent excavations
[J].
基坑开挖对临近明挖暗埋隧道竖向变形的影响机理
[J].
DOI:10.16183/j.cnki.jsjtu.2018.11.004
[本文引用: 1]
结合上海地区某邻近明挖暗埋隧道基坑工程案例,采用数值方法模拟了深基坑开挖的过程,并通过与实测数据进行对比验证了模型的合理性.针对隧道与基坑之间的连接墙进行对比分析,揭示了连接墙所起的作用.通过分析相邻地连墙竖向位移以及其两侧土体位移,得出地连墙与隧道之间以及墙土之间的相互影响规律.根据地连墙墙体轴力分布,研究墙体在开挖过程中的受力形态,揭示临近基坑开挖对地下结构(明挖暗埋隧道)的影响机理.结果表明,对于明挖隧道与基坑共墙的情况,隧道的隆起主要由地连墙以及连接墙的共同影响所致.明挖隧道受到基坑地连墙变形的影响明显大于受周围土体位移的影响.地连墙发生竖向隆起的原因为坑内土体提供的摩擦力大于坑外土体提供的摩擦力.
Mechanism analysis of influence of deep excavation on deformation of nearby cut-and-cover tunnel
[J].
地铁车站深基坑开挖变形智能多步预测方法
[J].
DOI:10.16183/j.cnki.jsjtu.2022.419
[本文引用: 1]
为更好预测深基坑开挖期间地下连续墙的侧向位移变形,基于长短期记忆神经网络(LSTM)智能算法理论构建了LSTM多步预测模型.首先对多步预测模型的多输出策略进行论述,其次详细介绍了LSTM多步预测模型的构建方法,并对模型输入集空间维度和时间维度两项超参数进行探究,以提高模型的预测准确度.最后依托某富水砂土深基坑工程实例,分析了模型预测值与实际监测值的差异.3个典型监测点的数据分析结果表明LSTM多步预测模型具有较强的泛化能力,相关算法对深基坑开挖变形预测方法的改进优化具有参考价值.
Multi-step prediction of excavation deformation of subway station based on intelligent algorithm
[J].
上海地区深基坑周边地表变形性状实测统计分析
[J].
Statistical analysis of characteristics of ground surface settlement caused by deep excavations in Shanghai soft soils
[J].
基于实测数据的上海地区超深基坑变形特性研究
[J].
Deformation characteristics of ultra-deep foundation pit in Shanghai based on measured data
[J].
杭州深厚软黏土中某深大基坑的性状研究
[J].
Characteristics of a large and deep soft clay excavation in Hangzhou
[J].
Patterns of soil deformations around tunnels. Application to the extension of Madrid Metro
[J].
Analytical prediction for ground movements due to deep excavations in soils
[J].
A displacement controlled method for evaluating ground settlement induced by excavation in clay
[J].
挡墙水平变位诱发地表沉降的显式解析解
[J].
Analytical formula for ground settlement induced by horizontal movement of retaining wall
[J].
Closed-form solution for excavation-induced ground settlement profile in clay
[J].
A theoretical study on ground surface settlement induced by a braced deep excavation
[J].
Small strain based method for predicting three-dimensional soil displacements induced by braced excavation
[J].
Numerical study of the effects of groundwater drawdown on ground settlement for excavation in residual soils
[J].
An efficient method for evaluating the ground surface settlement of Hangzhou metro deep basement considering the excavation process
[J].
Observed performance of a 30.2 m deep-large basement excavation in Hangzhou soft clay
[J].
上海土体小应变硬化模型整套参数取值方法及工程验证
[J].
On HSS model parameters for Shanghai soils with engineering verification
[J].
复杂环境下某深厚软土基坑的实测性状研究
[J].
Measured characteristics of a deep soft soil excavation in complex environment
[J].
软土基坑逆作开挖的时效分析
[J].
Time effect analysis of top-down excavation of foundation pits in soft clay
[J].
Distributions of ground movements parallel to deep excavations in clay
[J].
/
〈 |
|
〉 |
