1. 引言
在邻近地铁区间隧道区域进行大直径钻孔灌注桩施工时,常常采用三轴深层水泥土搅拌桩预先加固钻孔灌注桩桩周地层,以防止孔壁坍塌。但是,三轴水泥土搅拌桩施工时喷浆压力大、下沉速度大,挤土效应显著,往往导致邻近地下管线和构筑物的损坏。因此,针对三轴搅拌桩施工,需要首先对其挤土效应进行预测,为选择合适的喷浆压力及沉桩速度提供依据。
对于桩基施工及其它岩土施工活动如盾构(顶管)隧道掘进等的挤土效应,预测分析的方法有很多,其中基于Mindlin基本解 [1] 的理论解析方法具有快速易学的优点,是最为常用的理论解析方法之一。至今,采用此法的工程实践积累了较为丰富的研究成果。
魏刚,徐日庆(2005)基于Mindlin基本解推导了盾构隧道掌子面荷载作用下土体的应力和位移计算公式,提出了软土隧道盾构法施工引起的纵向地面变形预测的方法 [2];许有俊等人(2019)针对土压平衡矩形顶管正面附加推力作用,采用基于Mindlin基本解导出的土体位移公式,对地表隆起变形影响进行了分析研究,并提出了预测顶管施工过程地表隆起变形的方法 [3];秦世伟等人(2013)从Mindlin基本解导出了静压桩沉桩施工时小孔扩张应力下半无限体内任一点处的应力和位移计算公式,并据此就静压桩沉桩施工对临近隧道的影响进行了分析研究,提出了保护隧道的措施 [4];谷冠思(2019)针对钻孔灌注桩基桩施工全过程对临近地铁隧道土体扰动行为,采用Mindlin基本解进行了研究,在此基础上提出了控制孔壁坍塌的建议 [5]。上述学者的研究方法为三轴深层水泥土搅拌桩施工的挤土效应预测分析提供了有益的借鉴。
本文依托上海浦东创新中路附近的某地铁二号线跨线桥工程,采用基于Mindlin基本解的理论解析方法就三轴水泥土搅拌桩施工的挤土效应进行研究,推导了三轴搅拌桩施工时邻近土体的位移计算公式,并以本工程作为案例进行计算分析。
2. 三轴搅拌桩施工挤土效应的Mindlin理论解析
2.1. Mindlin弹性基本解公式
事实上,所谓Mindlin弹性解就是在Boussinesq解答的基础上推导出竖向或水平向荷载作用在半无限弹性体内部时,半无限体内任一点的应力与位移。图1为水平向荷载作用在半无限弹性体内部时Mindlin弹性解示意图。1936年,Mindlin [1] 得到了均匀各向同性体三维弹性方程的解。
![](//html.hanspub.org/file/10-2751239x9_hanspub.png?20210826084223994)
Figure 1. Sketch of Mindlin model
图1. Mindlin解模型示意图
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
式中,
,
;
c——水平力作用点与土体表面的垂直距离(m);
z——计算点与土体表面间的垂直距离(m);
x——水平力作用点与计算点之间在X轴上的投影距离(m);
P——水平集中力的大小(N);
G——土体剪切模量(MPa),可通过弹性模量E折算
;
——土体泊松比。
2.2. 三轴水泥土搅拌桩挤土效应的Mindlin理论解析
2.2.1. 简化假定
三轴水泥土搅拌桩施工过程中,喷浆压力和沉桩压力均会对桩侧土体产生压力。搅拌桩的挤土效应实质上就是在半无限土体中作用了水平荷载,因此,可以采用 Mindlin理论求解。图2为三轴搅拌桩Mindlin理论解析的模型示意图。
为了分析验算方便,特做如下简化假定:
1) 所研究的土体为各向同性线弹性体;
2) 喷浆压力和沉桩钻进压力合并成单一荷载考虑;
3) 三轴搅拌桩横截面等效成外包圆截面,每次加载段长度固定,为d1;
4) 假定搅拌桩施工过程荷载分段施加,位移场分段累加:前一阶段为后一阶段的初始位移场。
因此,求解某一段搅拌桩施工荷载对计算点的位移时,需要先按圆周对Mindlin理论解进行积分,然后对各段结果进行位移场叠加。
![](//html.hanspub.org/file/10-2751239x23_hanspub.png?20210826084223994)
Figure 2. Mindlin model sketch of TMCP
图2. 三轴搅拌桩 Mindlin解模型示意图
2.2.2. 单段加载时的位移公式
先求算分段情况下的位移场。图3和图4为圆形桩与计算点的位置关系计算简图(俯视图)。综合图3和图4,r的计算在转角
取任意值都可按式(10)计算。
(10)
将式(7)、式(9)中的
替换为
,令:
(11)
(12)
设搅拌桩各加载段长度为
,加载压力为p,则
(13)
将式(13)替换式(7)、式(9)中的P,并按圆周积分,可得某段加载对计算点的水平和竖向位移为:
![](//html.hanspub.org/file/10-2751239x32_hanspub.png?20210826084223994)
Figure 3. Calculation sketch at the interval
图3.
时计算简图
![](//html.hanspub.org/file/10-2751239x35_hanspub.png?20210826084223994)
Figure 4. Calculation sketch at the interval
图4.
时计算简图
(14)
(15)
式中,
——搅拌桩的桩头加载段长度,本工程取0.5 m;
p——搅拌桩的桩头荷载集度(Pa)。
从式(14)和式(15)可以看出,邻近土体的水平向位移
和竖向位移
与等效圆圆周上的压力p密切相关,也与施工时加载段长度
有关。而加载段长度
是由搅拌桩机械设备决定的,为一定值,即一旦施工设备确定,则加载段长度
也即确定。所以,控制邻近土体的水平向位移
和竖向位移
只要控制压力p就行了。
所以Mindlin理论解可以十分方便地匡算出搅拌桩挤土效应(即邻近土体的水平向位移
和竖向位移
),以此作为依据,调整相应的注浆压力和沉桩速度,以达到确保周边环境安全之目的。
2.2.3. 整桩施工时的位移公式
将搅拌全过程切分为不同个小阶段,每个阶段产生的位移场成为下一阶段的初始位移场。整个过程在土体内确定一点产生的位移可以表示为:
(16)
式中,
为每个小阶段喷浆压力对确定一点土体造成的位移;
n——划分的小阶段数目。
3. 案例分析
3.1. 工程概况
计算案例依托于上海市浦东新区创新中路某跨线桥工程 [6] [7] [8]。该桥跨越正在运营的地铁2号线区间隧道,跨线桥的桥墩基桩采用钻孔灌注桩。为防止钻孔灌注桩成孔施工时发生坍孔,先采用三轴水泥土搅拌桩对钻孔灌注桩区域的地层进行地基加固,三轴搅拌桩直径850 mm,加固深度22 m。搅拌桩和区间隧道最小距离为2.83米。桩基与盾构隧道位置关系如图5所示。
![](//html.hanspub.org/file/10-2751239x52_hanspub.png?20210826084223994)
Figure 5. Subway tunnel and pile foundation of bridge
图5. 桥梁桩基和地铁隧道的位置关系
3.2. 地质参数
地基土土层自地面而下分布如下:① 填土,② 粉质黏土,③ 淤泥质粉质黏土,③ 夹层粉质黏土,④ 淤泥质粘土,⑤ 1粘土,⑤ 3-1粉质黏土,⑤ 3-2粘质粉土,⑤ 3-3粉质黏土,⑤ 3-4粉质黏土夹粉性土,⑦ 2粉砂,⑧ 1粉质黏土,⑧ 2-1粉质黏土粉砂互层。
各层土的计算参数如表1所示。
![](Images/Table_Tmp.jpg)
Table 1. Properties of soils for Mindlin equation
表1. 土层的主要计算参数一览表
3.3. 现场量测
为了了解三轴搅拌桩施工对邻近土体的影响,施工期间在邻近隧道边缘布置了测斜管,以便获得深层土体的水平位移。测点布置如图6所示。测斜管深度为22米,距离隧道外缘0.5 m。
![](//html.hanspub.org/file/10-2751239x53_hanspub.png?20210826084223994)
Figure 6. Layout of inclinometer casing
图6. 测斜管布置示意图
3.4. 计算结果及其分析
根据荷载情况、土层情况,按式(13)和(14)采用Mathcad软件或Matlab软件编程计算,可得测斜管水平位置沿深度的水平位移和竖向位移。计算时土层弹性模量及泊松比取值如表1所示。
3.4.1. 搅拌桩施工导致邻近土体的水平位移
计算得到搅拌桩施工对邻近土体的水平位移位移和现场实测水平位移,如图7所示。结果表明,计算值在4 m埋深处达到最大值3 mm,实测值也在大致深度处达到最大值4.7 mm。位移变化规律是:理论分析与现场实验相比,位移沿深度变化趋势相同,位移峰值点出现在较浅层土体(0~5 m)处。理论解得出的峰值位移偏小,可能是理论计算时选取的喷浆压力偏小所致。
![](//html.hanspub.org/file/10-2751239x54_hanspub.png?20210826084223994)
Figure 7. Horizontal displacement of surrounding soil induced to TMCP driving
图7. 搅拌桩施工导致邻近土体的水平位移
3.4.2. 搅拌桩施工导致邻近土体的竖向位移
计算得到搅拌桩施工对邻近土体的竖向位移,如图8所示。结果表明,在浅层土体表现为土体隆起,随深度增加,位移减小,直至反向达到峰值,随后沉降又逐渐减弱。
![](//html.hanspub.org/file/10-2751239x55_hanspub.png?20210826084223994)
Figure 8. Vertical displacement of surrounding soil induced to TMCP driving
图8. 搅拌桩施工对土体竖向位移影响
3.4.3. 计算结果讨论
Mindlin基本解是在假定土体为各向同性的线弹性体基础上获得的,其解答适合于小变形情况。实际上土体是各向异性的弹塑性体,因此采用Mindlin公式得到的计算结果是近似的。
由于三轴搅拌桩施工对土体的扰动是在深层土体中发生的,总体上符合小变形假定。对于变形相对微小(即实际变形远小于结构如隧道的尺寸)情况,只要土层参数准确,计算结果的误差在工程允许范围内,其结论对于施工参数的初步选择仍具有指导意义,有助于快速了解注浆压力及沉桩速度对邻近土体位移的影响程度。
4. 结论
1) 搅拌桩施工时的挤土效应本质上就是施工过程中在大地半无限体内某一区域作用了分布荷载,因此,可以采用Mindlin解来预测邻近土体的位移;
2) 三轴搅拌桩的荷载效应计算时假定在其外包等效圆圆周上作用了径向分布力p,然后依此对Mindlin公式进行积分得到搅拌桩施工全过程中邻近土体的位移;
3) 尽管Mindlin解是在各向同性线弹性体假定基础上得到的,它仍可以真实反映搅拌桩施工时注浆压力及沉桩速度与挤土效应(及邻近土体的位移)的关系,可以帮助现场工程师快速获得施工参数控制的关键信息,这是采用复杂的三维有限元数值计算软件难以做到的。
致谢
本文公式推导得到同济大学赖允瑾和谷冠思两位同志的帮助,特此致谢。
基金项目
2019年上海浦东新区城建系统科学技术研究项目,项目编号:2019-015。