基于Mindlin解的搅拌桩施工挤土效应预测分析——上海某近地铁区间隧道工程案例分析
Squeezing Effect Prediction and Analysis of Tri-Axial Cement-Mixing Pile (TCMP) Based on Mindlin Equation—A Case Study of a Project near the Subway Tunnel in Shanghai

作者: 龚灵迪 , 王自民 , 孙久平 :上海两港市政工程有限公司,上海;

关键词: 三轴水泥土搅拌桩挤土效应Mindlin理论解土体位移Tri-Axial Cement-Mixing Pile (TCMP) Squeezing Effect of Pile-Driving Mindlin Equation Soil Movement

摘要: 水泥土搅拌桩是地基处理的常用方法。在大直径钻孔灌注桩成孔施工时,为了防止坍孔,以确保邻近地层稳定,经常采用三轴水泥土搅拌桩对钻孔桩周围地层进行加固处理。但三轴水泥土搅拌桩沉桩施工因为其高压注浆和快速钻进原因,会对邻近土层产生挤压作用,这种挤压作用将导致邻近土体的位移。为此,需要根据邻近地铁区间隧道的运营维护要求,适当控制搅拌桩的注浆压力和沉桩速度,将搅拌桩的挤土不利影响降至最低程度。本文基于Mindlin理论,采用理论解析的方法,就三轴水泥土搅拌桩施工对邻近土体及邻近隧道的影响进行了公式推导和计算分析,并将计算结果和试验监测结果进行了对比。本案例表明,采用理论解析方法来预测搅拌桩的挤土效应具有快速方便的特点,在搅拌桩施工参数选择决策中,不失为一种有益的方法。

Abstract: Tri-axial cement-mixing pile (TCMP) is a very common option for soil improvement. During pile driving of bored pile with large diameter, TCMP is often adopted in order to avoid caving of bored hole, and accordingly to ensure the safety of the neighbouring subway tunnel. But unfortunately, TCMP can also lead to rather large movement of surrounding soil due to the high grouting pressure and rapid boring speed when TCMP driving, which is called pile-driving’s squeezing effect. To meet the requirement of the displacement of subway tunnel, it is necessary to control the grouting pressure and boring speed of TCMP in order to reduce the adverse effect to the minimum. Based on Mindlin Equation, this paper has derived the analytical formula of the displacement induced to TCMP driving and carried out relevant calculation, and finally, this paper has compared this calculation with the site measurement. This paper has shown that this analytical method, being fast and convenient, is a recommendable one in predicting the squeezing effect of TCMP driving, and is very helpful in selecting appropriate grouting pressure and boring speed of TCMP.

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] 得到了均匀各向同性体三维弹性方程的解。

Figure 1. Sketch of Mindlin model

图1. Mindlin解模型示意图

σ x = P x 8 π ( 1 ν ) { 1 2 ν R 1 3 + 3 x 2 R 1 5 ( 1 2 ν ) ( 5 4 ν ) R 2 3 + 3 ( 3 4 ν ) x 2 R 2 5 + 4 ( 1 ν ) ( 1 2 ν ) R 2 ( R 2 + z + c ) 2 [ 3 x 2 ( 3 R 2 + z + c ) R 2 2 ( R 2 + z + c ) ] 6 c R 2 5 [ 3 c ( 3 2 ν ) ( z + c ) + 5 x 2 z R 2 2 ] } (1)

σ y = P x 8 π ( 1 ν ) { 1 2 ν R 1 3 + 3 y 2 R 1 5 ( 1 2 ν ) ( 3 4 ν ) R 2 3 + 3 ( 3 4 ν ) y 2 R 2 5 + 4 ( 1 ν ) ( 1 2 ν ) R 2 ( R 2 + z + c ) 2 [ 1 y 2 ( 3 R 2 + z + c ) R 2 2 ( R 2 + z + c ) ] 6 c R 2 5 [ c ( 1 2 ν ) ( z + c ) + 5 y z 2 R 2 2 ] } (2)

σ z = P x 8 π ( 1 ν ) { 1 2 ν R 1 3 + 3 ( z c ) 2 R 1 5 1 2 ν R 2 3 + 3 ( 3 4 ν ) ( z + c ) 2 R 2 5 6 c R 2 5 [ c + ( 1 2 ν ) ( z + c ) + 5 z ( z + c ) 2 R 2 2 ] } (3)

τ x y = P y 8 π ( 1 ν ) { 1 2 ν R 1 3 + 3 x 2 R 1 5 1 2 ν R 2 3 + 3 ( 3 4 ν ) x 2 R 2 5 + 4 ( 1 ν ) ( 1 2 ν ) R 2 ( R 2 + z + c ) 2 + [ 1 x 2 ( 3 R 2 + z + c ) R 2 2 ( R 2 + z + c ) ] 6 c z R 2 5 [ 1 5 x 2 R 2 2 ] } (4)

τ y z = P x y 8 π ( 1 ν ) { 3 ( z c ) R 1 5 + 3 ( 3 4 ν ) ( z + c ) R 2 5 6 c R 2 5 [ 1 2 ν + 5 z ( z + c ) R 2 2 ] } (5)

τ z x = P 8 π ( 1 ν ) { ( 1 2 ν ) ( z c ) R 1 3 + 3 x 3 ( z c ) R 1 5 ( 1 2 ν ) ( z c ) R 2 3 + 3 ( 3 ν ) x 2 ( z + c ) R 2 5 + 6 c R 2 5 [ z ( z + c ) ( 1 2 ν ) z 2 5 x z ( z + c ) R 2 2 ] } (6)

u x = P 16 π G ( 1 ν ) { 3 4 ν R 1 + x 2 R 1 3 + 1 R 2 + ( 3 4 ν ) x 2 R 2 3 + 2 c z R 2 3 ( 1 3 x 2 R 2 2 ) + 4 ( 1 ν ) ( 1 2 ν ) R 2 + z + c [ 1 x 2 R 2 ( R 2 + z + c ) ] } (7)

u y = P x y 16 π G ( 1 ν ) { 1 R 1 3 + 3 4 ν R 2 3 6 c z R 2 3 4 ( 1 ν ) ( 1 2 ν ) R 2 ( R 2 + z + c ) 2 } (8)

u z = P x 16 π G ( 1 ν ) { z c R 1 3 + ( 3 4 ν ) ( z c ) R 2 3 6 c z ( z + c ) R 2 5 + 4 ( 1 ν ) ( 1 2 ν ) R 2 ( R 2 + z + c ) } (9)

式中, R 1 = [ r 2 + ( z c ) 2 ] 1 2 R 2 = [ r 2 + ( z + c ) 2 ] 1 2

c——水平力作用点与土体表面的垂直距离(m);

z——计算点与土体表面间的垂直距离(m);

x——水平力作用点与计算点之间在X轴上的投影距离(m);

P——水平集中力的大小(N);

G——土体剪切模量(MPa),可通过弹性模量E折算 G = E 2 ( 1 + ν )

ν ——土体泊松比。

2.2. 三轴水泥土搅拌桩挤土效应的Mindlin理论解析

2.2.1. 简化假定

三轴水泥土搅拌桩施工过程中,喷浆压力和沉桩压力均会对桩侧土体产生压力。搅拌桩的挤土效应实质上就是在半无限土体中作用了水平荷载,因此,可以采用 Mindlin理论求解。图2为三轴搅拌桩Mindlin理论解析的模型示意图。

为了分析验算方便,特做如下简化假定:

1) 所研究的土体为各向同性线弹性体;

2) 喷浆压力和沉桩钻进压力合并成单一荷载考虑;

3) 三轴搅拌桩横截面等效成外包圆截面,每次加载段长度固定,为d1

4) 假定搅拌桩施工过程荷载分段施加,位移场分段累加:前一阶段为后一阶段的初始位移场。

因此,求解某一段搅拌桩施工荷载对计算点的位移时,需要先按圆周对Mindlin理论解进行积分,然后对各段结果进行位移场叠加。

Figure 2. Mindlin model sketch of TMCP

图2. 三轴搅拌桩 Mindlin解模型示意图

2.2.2. 单段加载时的位移公式

先求算分段情况下的位移场。图3图4为圆形桩与计算点的位置关系计算简图(俯视图)。综合图3图4,r的计算在转角 θ 取任意值都可按式(10)计算。

r = ( r 0 2 + ( r 0 + r 1 ) 2 2 ( r 0 + r 1 ) r 0 cos θ ) 1 2 (10)

将式(7)、式(9)中的 x 2 替换为 [ ( r 0 + r 1 ) cos θ r 0 ] 2 ,令:

R 1 = [ r 0 2 + ( r 0 + r 1 ) 2 2 ( r 0 + r 1 ) r 0 cos θ + ( z c ) 2 ] 1 2 (11)

R 2 = [ r 0 2 + ( r 0 + r 1 ) 2 2 ( r 0 + r 1 ) r 0 cos θ + ( z + c ) 2 ] 1 2 (12)

设搅拌桩各加载段长度为 d 1 ,加载压力为p,则

d P = p d 1 r d θ (13)

将式(13)替换式(7)、式(9)中的P,并按圆周积分,可得某段加载对计算点的水平和竖向位移为:

Figure 3. Calculation sketch at the interval θ [ 90 , 90 ]

图3. θ [ 90 , 90 ] 时计算简图

Figure 4. Calculation sketch at the interval θ [ 90 , 270 ]

图4. θ [ 90 , 270 ] 时计算简图

U x ( z ) = p d 1 r 16 π G ( 1 ν ) 0 2 π { 3 4 ν R 1 + [ ( r 0 + r 1 ) cos θ r 0 ] 2 R 1 3 + 1 R 2 + ( 3 4 ν ) [ ( r 0 + r 1 ) cos θ r 0 ] 2 R 2 3 + 2 c z R 2 3 ( 1 3 [ ( r 0 + r 1 ) cos θ r 0 ] 2 R 2 2 ) + 4 ( 1 ν ) ( 1 2 ν ) R 2 + z + c [ 1 [ ( r 0 + r 1 ) cos θ r 0 ] 2 R 2 ( R 2 + z + c ) ] } d θ (14)

U z ( z ) = p d 1 r 16 π G ( 1 ν ) 0 2 π [ ( r 0 + r 1 ) cos θ r 0 ] { z c R 1 3 + ( 3 4 ν ) ( z c ) R 2 3 6 c z ( z + c ) R 2 5 + 4 ( 1 ν ) ( 1 2 ν ) R 2 ( R 2 + z + c ) } d θ (15)

式中, d 1 ——搅拌桩的桩头加载段长度,本工程取0.5 m;

p——搅拌桩的桩头荷载集度(Pa)。

从式(14)和式(15)可以看出,邻近土体的水平向位移 U x ( z ) 和竖向位移 U z ( z ) 与等效圆圆周上的压力p密切相关,也与施工时加载段长度 d 1 有关。而加载段长度 d 1 是由搅拌桩机械设备决定的,为一定值,即一旦施工设备确定,则加载段长度 d 1 也即确定。所以,控制邻近土体的水平向位移 U x ( z ) 和竖向位移 U z ( z ) 只要控制压力p就行了。

所以Mindlin理论解可以十分方便地匡算出搅拌桩挤土效应(即邻近土体的水平向位移 U x ( z ) 和竖向位移 U z ( z ) ),以此作为依据,调整相应的注浆压力和沉桩速度,以达到确保周边环境安全之目的。

2.2.3. 整桩施工时的位移公式

将搅拌全过程切分为不同个小阶段,每个阶段产生的位移场成为下一阶段的初始位移场。整个过程在土体内确定一点产生的位移可以表示为:

u t o t a l = u 1 + u 2 + u 3 + + u i + + u n (16)

式中, u i 为每个小阶段喷浆压力对确定一点土体造成的位移;

n——划分的小阶段数目。

3. 案例分析

3.1. 工程概况

计算案例依托于上海市浦东新区创新中路某跨线桥工程 [6] [7] [8]。该桥跨越正在运营的地铁2号线区间隧道,跨线桥的桥墩基桩采用钻孔灌注桩。为防止钻孔灌注桩成孔施工时发生坍孔,先采用三轴水泥土搅拌桩对钻孔灌注桩区域的地层进行地基加固,三轴搅拌桩直径850 mm,加固深度22 m。搅拌桩和区间隧道最小距离为2.83米。桩基与盾构隧道位置关系如图5所示。

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所示。

Table 1. Properties of soils for Mindlin equation

表1. 土层的主要计算参数一览表

3.3. 现场量测

为了了解三轴搅拌桩施工对邻近土体的影响,施工期间在邻近隧道边缘布置了测斜管,以便获得深层土体的水平位移。测点布置如图6所示。测斜管深度为22米,距离隧道外缘0.5 m。

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)处。理论解得出的峰值位移偏小,可能是理论计算时选取的喷浆压力偏小所致。

Figure 7. Horizontal displacement of surrounding soil induced to TMCP driving

图7. 搅拌桩施工导致邻近土体的水平位移

3.4.2. 搅拌桩施工导致邻近土体的竖向位移

计算得到搅拌桩施工对邻近土体的竖向位移,如图8所示。结果表明,在浅层土体表现为土体隆起,随深度增加,位移减小,直至反向达到峰值,随后沉降又逐渐减弱。

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。

文章引用: 龚灵迪 , 王自民 , 孙久平 (2021) 基于Mindlin解的搅拌桩施工挤土效应预测分析——上海某近地铁区间隧道工程案例分析。 土木工程, 10, 787-796. doi: 10.12677/HJCE.2021.108089

参考文献

[1] Mindlin, R.D. (1936) Force at a Point in the Interior of a Semi Infinite Solid. Physics, 7, 195-202.
https://doi.org/10.1063/1.1745385

[2] 魏纲, 徐日庆. 软土隧道盾构法施工引起的纵向地面变形预测[J]. 岩土工程学报, 2005, 27(9): 1077-1081.

[3] 许有俊, 张治华, 史明, 等. 土压平衡矩形顶管正面附加推力对地表隆起变形影响研究[J]. 西安建筑科技大学学报(自然科学版), 2019, 51(6): 814-824.

[4] 秦世伟, 周艳坤, 莫泷. 静压桩沉桩施工对临近隧道的影响[J]. 上海大学学报(自然科学版), 2013, 19(5): 527-533.

[5] 谷冠思. 桩基施工全过程对临近地铁隧道土体扰动行为研究[D]: [硕士学位论文], 上海: 同济大学, 2019.

[6] 龚灵迪, 周萌, 赖允瑾, 谷冠思. 三轴搅拌桩施工挤土效应的试验研究[J]. 土木工程, 2020, 9(9): 864-874.
https://doi.org/10.12677/HJCE.2020.99091

[7] 汤新春, 张荣, 赖允瑾, 谷冠思. 搅拌桩和钻孔桩先后施工对邻近地铁隧道的影响分析[J]. 土木工程, 2020, 9(9): 853-863.
https://doi.org/10.12677/HJCE.2020.99090

[8] 黄芹, 陈耀华, 赖允瑾. 邻近地铁隧道钻孔灌注桩套打施工的防坍孔力学机理及工程应用[J], 现代隧道技术, 2018, 55(S2): 419-425.

分享
Top