文章快速检索    
  同济大学学报(自然科学版)  2018, Vol. 46 Issue (10): 1384-1391.  DOI: 10.11908/j.issn.0253-374x.2018.10.009
0

引用本文  

宫全美, 张润来, 余杰, 周顺华. 桩底盾构施工引起的桩基承载力损失计算[J]. 同济大学学报(自然科学版), 2018, 46(10): 1384-1391. DOI: 10.11908/j.issn.0253-374x.2018.10.009.
GONG Quanmei, ZHANG Runlai, YU Jie, ZHOU Shunhua. Calculating Method of Pile Bearing Capacity Loss Induced by the Tunnel Excavation Underneath the Pile[J]. Journal of Tongji University (Natural Science), 2018, 46(10): 1384-1391. DOI: 10.11908/j.issn.0253-374x.2018.10.009

基金项目

国家自然科学基金(41472247)

第一作者

宫全美(1967—),女,教授,博士生导师,工学博士,主要研究方向为铁道工程、线路动力学.E-mail: gongqm@tongji.edu.cn

通信作者

张润来(1991—),男,博士生,主要研究方向为盾构隧道.E-mail: 1510061@tongji.edu.cn

文章历史

收稿日期:2017-09-29
桩底盾构施工引起的桩基承载力损失计算
宫全美 , 张润来 , 余杰 , 周顺华     
同济大学 道路与交通工程教育部重点实验室,上海 201804
摘要:提出一种桩底盾构施工引起的桩基承载能力损失计算方法.先实现考虑卸载过程的双曲线荷载传递函数编译,然后利用该桩土接触模型得到隧道开挖前和开挖后的Q~s曲线;取s=50 mm时对应的荷载为桩基极限承载能力,采用桩基极限承载力损失百分比作为隧道开挖对桩基承载性能影响的评价指标.结合杭州地铁1号线某工点实例,分析了桩基承载曲线变化特征、桩体内力变化规律和承载力损失影响因素.案例中,盾构施工体积损失率控制在0.5%时,桩基承载力损失值为22%.隧道开挖后,桩体中存在一点侧摩阻力不变,在该点上部桩体侧摩阻力增大,在该点下部桩体侧摩阻力有所减小.承载力损失值随着体积损失率的增大而增大;桩基的初始荷载水平越大,承载力损失值越大;桩底与隧道顶部的距离越大,桩基承载力损失值越小.
关键词桩基    盾构隧道    承载力    
Calculating Method of Pile Bearing Capacity Loss Induced by the Tunnel Excavation Underneath the Pile
GONG Quanmei , ZHANG Runlai , YU Jie , ZHOU Shunhua     
Key Laboratory of Road and Traffic Engineering of the Ministry of Education, Tongii University, Shanghai 201804, China
Abstract: A method was proposed to calculate the pile bearing capacity loss induced by the tunnel excavation underneath the pile. The hyperbolic load-transfer function considering the effect of tunnel excavation at the pile-soil interface was compiled and based on it, the Q~s curves of the pile and that after tunnel excavation were figured out. Ultimate bearing capacity of pile was confirmed where pile settlement was fifty millimeter, and it was used as the indicator to measure the influence caused by the excavation. The variation of bearing curve and pile forces as well as influencing factors to the bearing capacity were analyzed based on one work point of Hangzhou Metro Line 1. The bearing capacity loss of pile is 22% when volume loss rate is controlled at 0.5%. After tunnel excavation, there is a position in the pile where the lateral friction is unchanged, while the lateral friction increases upper the position and decreases bellow the position. The bearing capacity loss has a positive correlation with volume loss rate and Initial load of pile, however, has a negative correlation with the distance between the bottom of the pile and the top of the tunnel.
Key words: pile foundation    shield tunnel    bearing capacity    

地铁盾构隧道穿越建筑物不可避免地会下穿建筑物桩基,从而造成桩底、桩侧土应力场发生改变,会对既有建筑物桩基的变形和受力情况产生不利的影响,从而影响上部结构安全.

目前,国内外学者在隧道开挖对建筑物桩基影响研究上主要采用的方法有现场实测法、理论解析法、模型试验法和数值法.数值法由于其适应性强,分析结果详细,数据查看方便,被更广泛应用.Mroueh等[1]、Lee等[2]在有限元中分别考虑了土体应力释放因素与隧道桩基的相互作用.朱逢斌等[3]通过数值模拟和离心试验结果的对比分析,证明了采用数值模拟的方法来进行盾构隧道开挖对邻近桩基影响的研究是可行的.霍军帅等[4]用有限元法分析盾构隧道施工对城际铁路桩基的影响并对加固方案进行分析.王卫东等[5]对有限元中桩土接触单元所采用的接触模型进了修正,将桩土界面的双曲线模型嵌套到ABAQUS软件中.袁海平等[6]用有限差分法分析盾构推进过程中不同工况下桥桩结构受力、水平变形、地层沉降的变化规律.陈明等[7]等对桩土界面的荷载传递模型进行了改进,考虑了基坑开挖引起的卸荷效应.李新星等[8]采用有限元法模拟隧道近距离穿越高架桩基的施工过程并提出保护措施.

综合目前数值分析法的研究结果,首先,主要集中在桩基变形、受力分析和盾构施工措施方面,较少涉及隧道下穿施工引起的桩基承载能力损失研究,尚缺乏较为直观的承载力损失分析方法与评判指标.其次,目前盾构侧穿桩基研究较多,然而盾构桩底下穿施工对桩基影响较大,但桩底下穿施工研究较少.最后,现在盾构施工对桩基影响的数值模拟研究大都直接采用商业数值软件中内置的桩土接触模型,即理想弹塑性或者刚塑性模型,在卸载阶段剪应力与相对位移之间的关系与加载段相同,这与隧道开挖过程中桩土接触卸载特征不同,应采用能够考虑卸载过程的桩土接触模型来模拟这一过程.

因此,针对盾构桩底下穿施工情况展开研究,分析其对桩基承载性能的影响.桩土接触模型采用能够考虑卸载过程的双曲线荷载传递函数模型,在ABAQUS环境中进行编译实现,并验证其可靠性.提出明确的桩基承载力计算方法;分析隧道开挖引起桩基位移和受力变化;对盾构施工控制水平、桩基初始承载水平和桩底与隧道顶部距离3个影响因素进行分析.

1 考虑卸载过程双曲线荷载传递函数 1.1 考虑卸载过程的双曲线荷载传递函数介绍

Zhang等[9]在土-结构动力相互作用(SSI)循环剪切试验结果[10-11]提出的双曲线解析计算模型基础上,提出了考虑受隧道开挖影响以及桩土界面滑移特性的荷载传递函数.如图 1所示为循环SSI试验中剪应力τ与相对位移u的关系,如图 2所示为传递函数中剪应力和相对位移之间的双曲线关系,OAAB为加载段,AC为卸载段,CD为反向加载段,k0为界面的初始剪切刚度,δir为残余界面滑移值,τult(z)为极限侧摩阻力.

图 1 循环SSI试验中剪应力与相对位移的关系 Fig.1 Relationship between the shear stress and the relative displacement in the SSI experiments
图 2 考虑卸载过程的双曲线荷载传递函数 Fig.2 Hyperbolic load-transfer function considering unloading process

式(1)为双曲线关系解析计算式,用于计算侧摩阻力;式(2)计算界面的初始剪切刚度;式(3)计算极限侧摩阻力.

$ \tau \left( z \right) = \left\{ {\begin{array}{*{20}{l}} {\frac{\delta }{{1/{k_0} + \delta /{\tau _{{\rm{ult}}}}\left( z \right)}},}&{{\rm{加载}}OA{\rm{和}}AB}\\ {\tau \left( z \right) = {k_0}\left( {\delta - {\delta _{{\rm{ir}}}}} \right),}&{{\rm{卸载}}AC}\\ {\frac{{{\delta _{{\rm{ir}}}} - \delta }}{{1/{k_0} + \left( {{\delta _{{\rm{ir}}}} - \delta } \right)/{\tau _{{\rm{ult}}}}\left( z \right)}},}&{{\rm{反向加载}}CD} \end{array}} \right. $ (1)
$ {k_0} = \frac{{\eta {\tau _{{\rm{ult}}}}\left( z \right)}}{{\left( {1 - \eta } \right){\delta _2}}} $ (2)
$ {\tau _{{\rm{ult}}}}\left( z \right) = {c_{\rm{a}}} + {\delta _{\rm{m}}}\tan {\varphi _{\rm{a}}} $ (3)

式中:τ(z)为侧摩阻力;k0为界面的初始剪切刚度;τult(z)为极限侧摩阻力;δ为桩土相对位移;δir为残余界面滑移值;ca为界面黏聚力;σm为侧向压力;$φ$a为界面摩擦角;η为效率系数,是某一时刻桩土界面的剪切力与界面极限剪切力的比值;δ2表示桩土界面的剪切力与界面极限剪切力比值为η时所对应的桩土相对位移值.

1.2 考虑卸载过程的双曲线接触程序编译实现

费康等[12]开始采用古德曼(Goodman)接触模型来模拟接触面上的错动滑裂或者开裂.Goodman认为,接触面在受力之前完全吻合,接触面上剪切力和相对位移之间模拟为非线性弹性关系,接触面的本构关系可用式(4)表示.本文采用Goodman单元非线性特征来实现双曲线的荷载传递函数.将双曲线荷载传递函数关系式(1)对δ求一阶导数,如式(5)所示,可得图 2中曲线上任意一点的斜率,然后代入式(4)中.

$ \left\{ {\begin{array}{*{20}{c}} {\Delta {\tau _1}}\\ {\Delta {\tau _2}} \end{array}} \right\} = \left[ {\begin{array}{*{20}{c}} {{k_{{\rm{s1}}}}}&0\\ 0&{{k_{{\rm{s2}}}}} \end{array}} \right]\left\{ {\begin{array}{*{20}{c}} {\Delta {\gamma _1}}\\ {\Delta {\gamma _2}} \end{array}} \right\} $ (4)
$ \begin{array}{l} \frac{{{\rm{d}}{\tau _{n1}}\left( z \right)}}{{{\rm{d}}{\delta _{n1}}}} = {k_{sn1}} = \\ \left\{ {\begin{array}{*{20}{l}} {\frac{1}{{{k_0}{{\left( {\frac{1}{{{k_0}}} + \frac{{{\delta _{n1}}}}{{{\tau _{{\rm{ult}}}}\left( z \right)}}} \right)}^2}}}}&{OA{\rm{和}}AB}\\ {{k_0}}&{AC}\\ { - \frac{1}{{{k_0}{{\left( {\frac{1}{{{k_0}}} + \frac{{{{\left( {{\delta _{{\rm{ir}}}} - {\delta _{n1}}} \right)}^2}}}{{{\tau _{{\rm{ult}}}}\left( z \right)}}} \right)}^2}}}}&{CD} \end{array}} \right. \end{array} $ (5)

式中:Δτ1、Δτ2分别为方向1、方向2的剪应力;ks1ks1分别为方向1、方向2的刚度;Δγ1、Δγ2分别为方向1、方向2的相对位移;δn1(n=1, 2)分别为方向1、方向2的相对位移;τn1(z)(n=1, 2)分别为方向1、方向2的剪切力.

对于上述双曲线接触模型计算每一时步都需要进行受力状态判定.为了便于判断考虑卸载过程的双曲线接触函数上一点的受力状态,将式(5)中的全部δ用摩阻力τ的函数式进行表达,得到根据曲线上任意一点的τ值来确定该点斜率的表达式(6).

$ \left\{ \begin{array}{l} \left| {{f_1}\left( \tau \right)} \right| = \frac{{{k_0}{{\left( {\left( {{c_{\rm{a}}} + {\sigma _{\rm{m}}}\tan {\varphi _{\rm{a}}}} \right) - \tau \left( z \right)} \right)}^2}}}{{{{\left( {{c_{\rm{a}}} + {\sigma _{\rm{m}}}\tan {\varphi _{\rm{a}}}} \right)}^2}}}\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;OA{\rm{和}}AB{\rm{和}}CD\\ {f_2}\left( \tau \right) = {k_{\rm{0}}}\;\;\;AC \end{array} \right. $ (6)

式中:f1(τ)、f2(τ)为不同加卸载段的刚度.

图 3所示,点A(δ0τ0)处于加载阶段时,根据双曲线切线斜率计算所得的增加后的剪应力τ1如式(7),而根据双曲线的初始切线斜率计算所得的增加后的剪应力τ1′如式(8).为了保证加载阶段剪应力与相对位移关系为双曲线,则增加后的剪应力值取τ1τ1′中绝对值中较小值,即相对位移增加Δδ后的最终剪应力值τ3取值可按照式(9),对于卸载时同理.其中:τ0δ0为点A的剪切力和位移,Δδ为相对位移增加量,τ1δ1为点A加载后根据双曲线得到的剪切力和位移,τ1δ1为点A加载后根据初始刚度k0得到的剪切力和位移,τ3为加载后受力状态判定的最终剪应力取值,τ2δ2为点A卸载后根据双曲线得到的剪切力和位移,τ2δ2为点A卸载后根据初始刚度k0得到的剪切力和位移,kc为双曲线的斜率.

图 3 计算点受力状态判定示意 Fig.3 Schematic diagram about the stress determination of the calculation point
$ {\tau _1} = {\tau _0} + {k_{\rm{c}}}\Delta \delta $ (7)
$ {{\tau '}_1} = {\tau _0} + {k_0}\Delta \delta $ (8)
$ {\tau _3} = \left\{ {\begin{array}{*{20}{l}} {{\tau _1},}&{\left| {{\tau _1}} \right| < \left| {{{\tau '}_1}} \right|}\\ {{{\tau '}_1},}&{\left| {{\tau _1}} \right| \ge \left| {{{\tau '}_1}} \right|} \end{array}} \right. $ (9)

在应力积分算法选择时有多种,如基本增量法、四阶四段龙格库塔增量法.本文采用计算精度更高的四阶四段龙格库塔法.采用Fortran语言完成考虑卸载过程的双曲线接触程序编译.

2 双曲线接触程序可靠性验证 2.1 滑块验证

图 4所示,将一尺寸为1.0 m×1.0 m×1.0 m的三维弹性立方体(滑块A)放置在一尺寸为3.0 m×3.0 m×3.0 m的刚体上,刚体完全固定.刚体的弹性模量为1 GPa,泊松比取零,滑块弹性模量10 MPa,泊松比0.4;两滑块间的摩擦系数取0.5,效率系数η取0.8,接触采用考虑卸载过程的双曲线接触程序.整个过程分为4个阶段,第1阶段在滑块A的正上方施加压力P,取值为200 kPa,建立2个部件之间的接触;第2阶段使滑块A沿y轴正向水平移动0.01 m;第3阶段使滑块A沿y轴负方向运动0.01 m回到原点;第4阶段让滑块A从原点开始沿y轴负方向运动0.01 m.得到摩擦力与相对位移的关系如图 5所示.

图 4 两滑块的相对位置关系 Fig.4 Position of two sliding block
图 5 摩擦力与相对位移的关系 Fig.5 Relationship between the friction and relative displacement

图 5可知,在HE段(0 s~1 s)以及FG段(2 s~3 s)滑块受到的摩擦力以及滑块与固定刚体之间的相对位移呈双曲线关系;在EF段(1 s~2 s)滑块受到的摩擦力以及滑块与固定刚体之间的相对位移呈直线关系.

将根据理论公式计算所得的摩擦力、相对位移的结果也绘制在图 5中.看出采用编译的接触模型计算所得摩擦力与相对位移的关系曲线图与理论计算所得的曲线基本重合.由此可知,有限元计算结果与理论公式计算结果一致,说明考虑卸载的双曲线接触模型正确.

2.2 模型试验比较验证

为了验证程序的可靠性,采用了Ng等[13]接近桩底正下方单线隧道开挖对既有桩基影响的三维离心试验数据进行比较验证.隧道开挖通过对开挖部位的土体施加对应体积损失率的位移边界来实现.开挖引起的土体位移模式见图 6.土体损失率η与等效土体损失参数g之间的关系如式(10)所示.

$ \eta = \frac{{{S_{{\rm{loss}}}}}}{{\pi {R^2}}} = 2\int_{R - \frac{g}{\pi }\theta }^R {r{\rm{d}}r} \int_0^\pi {{\rm{d}}\theta } = \frac{g}{R} - \frac{{{g^2}}}{{3{R^2}}} $ (10)
图 6 隧道开挖引起的土体位移模式 Fig.6 Displacement mode of the soil caused by tunnel excavation

式中:η为土体损失率;g为等效土体损失参数;Sloss为土体损失体积;R为隧道半径.

接触模型验证中参数如表 1.土体采用弹塑性模型、摩尔库伦破坏准则;桩体采用线弹性本构模型.盾构开挖引起的体积损失率为1%,体积损失率1%对应的g值为0.031 m.

下载CSV 表 1 有限元计算中主要参数 Tab.1 Main parameters in finite element model

图 7为调用程序进行有限元计算得出的隧道开挖前后桩基的轴力分布曲线与离心试验得到的桩基轴力分布曲线对比图,其中D为桩直径.从图中可以看出,隧道开挖前后通过有限元分析所得到的桩基轴力沿深度的变化规律与离心试验得到的轴力分布规律一致且数值接近,说明程序可靠.

图 7 有限元计算轴力与离心试验轴力对比 Fig.7 Comparison of axial force between finite element method and centrifugal test
3 桩基承载力损失计算 3.1 桩基承载力损失表示方法

桩基的承载性能是在保证桩身材料强度的前提下用沉降来描述,由《桩基手册》可知,在单桩静载试验的Q~s曲线中,一般取沉降值40 mm~60 mm对应的荷载值为单桩的极限承载力.我国桥梁和地基基础设计规范等也根据上部结构受力要求,规定了桥墩和房屋桩基允许的沉降值,如《建筑地基基础设计规范》规定,均匀沉降量不超过50 mm.故取沉降等于50 mm对应的荷载Q为单桩的极限承载力.令Q1为隧道开挖前桩基的极限承载力,Q2为隧道开挖后桩基的极限承载力,隧道开挖具有卸载作用能够引起桩基沉降,那么会导致Q2小于Q1,以式(11)表示极限承载力损失值e.

$ e = \frac{{{Q_1} - {Q_2}}}{{{Q_1}}} \times 100\% $ (11)
3.2 隧道开挖引起的桩基承载曲线变化特征

以杭州地铁1号线某工点作为分析实例.在该工点中,地铁盾构从一办公楼正下方穿过,办公楼桩基采用钻孔灌注桩,共3排,横向间距2 m,纵向间距3 m,前后排3根,中间2根;桩长为8.5 m,桩直径0.5 m,盾构外径为6.2 m.盾构隧道与建筑物桩基的相对位置关系见图 8.为消除边界效应的影响,模型的长宽深为60 m×60 m×60 m.桩的直径为0.5 m,桩长8.5 m.土体采用摩尔库伦理想弹塑性模型,土体参数见表 2;桩体采用线弹性本构模型,由于在桩长范围内一共有2层土,故桩与土体的接触分两段施加,接触面参数见表 3.

图 8 隧道与桩基的横剖面(单位:m) Fig.8 Transverse section between the tunnel and pile(unit: m)
下载CSV 表 2 土层主要力学性质参数 Tab.2 Main mechanical parameters of soil layers
下载CSV 表 3 桩土接触面的参数设置 Tab.3 Parameters of pile-soil interface

通过数值模拟得到隧道开挖前和开挖后的桩基Q~s曲线.首先,在没有隧道施工时,进行静载试验,桩顶荷载分8级施加,每一级荷载大小为100 kN,总共加载到800 kN,提取数据得到隧道开挖前的桩基Q~s曲线.然后,在有隧道施工时,先在桩顶施加254 kN的初始工作荷载Q0,模拟桩基已有的受荷状态.第二步进行隧道开挖卸载,通过对开挖部位土体施加对应体积损失率的位移边界来实现,开挖引起的体积损失率取软土地区常用的最大值0.5%.第三步逐级施加桩顶荷载,每一级荷载大小为100 kN,总共加载到800 kN,提取数据得到开挖后的桩基Q~s曲线.

数值模拟结果见图 9OAB为隧道开挖前的桩基Q~s曲线,OACD为隧道开挖后的桩基Q~s曲线.可以看出在A点为在初始工作荷载作用发生的桩基沉降,A点所对应的荷载为桩基初始荷载,OA段为初始荷载导致的桩基沉降.从A点到C点桩基承受荷载没有发生变化,但桩基发生沉降,这是由隧道开挖卸载所引起的,因此以AC段表达隧道开挖过程.图中虚线表示沉降为50 mm,其与OAB曲线的交点所对应的荷载为开挖前的极限承载能力Q1,其与OACD曲线的交点所对应的荷载为开挖后的极限承载能力Q2.因此,图 9能充分反应隧道开挖前后桩基承载性能变化情况,可被用来计算隧道开挖引起的桩基承载力损失计算.

图 9 桩基的Q~s曲线 Fig.9 Q~s curve of the pile

图 9可知,隧道开挖前桩基的极限承载力Q1为508 kN.盾构施工体积损失率控制在0.5%时,隧道开挖后的桩基的极限承载力Q2为400.3 kN.由式(11)计算,得到因隧道开挖导致的桩基承载力损失值为22%.初始工作荷载Q0为254.0 kN,小于开挖后的桩基极限承载力Q2,因此,盾构施工保持该控制水平时,降低了部分桩基承载能力,但不影响桩基础正常使用.

3.3 隧道开挖引起桩基位移和受力变化

针对上述实例,对隧道开挖前后桩基位移和受力变化进行以下分析.以Sp代表桩基的沉降,Ss代表桩周土的沉降,图 10所示为开挖前后桩基与桩周土的沉降变化.隧道开挖前后桩基的总沉降分别为8.5 mm、27.8 mm,且桩的总沉降值均大于桩周土的总沉降.

图 10 桩基与桩周土体的变形 Fig.10 Deformation of pile and soil around pile

图 11所示为隧道开挖前以及开挖引起的桩土相对位移和相对位移增量沿深度分布的情况.可知在隧道开挖前以及开挖后各阶段中,相对位移δ(δ=Sp-Ss)均为正值.这意味着在整个桩长范围内,桩的沉降值均比桩周土的沉降值大,若规定侧摩阻力向上为正,则桩受到的均为正摩阻力的作用.但相对位移增量在0~1.05D范围内为正值,则开挖导致该范围侧摩阻力增大;相对位移增量在1.05D~1.40D范围内为负值,该范围内侧摩阻力有所减小.在1.05D处侧摩阻力大小不发生改变.

图 11 桩土相对位移和桩土相对位移增量分布 Fig.11 Distribution of the relative displacement and the relative displacement increment between pile and soil

图 12a中所示,侧摩阻力变化量与上述观点一致.此外,图 12a还表明侧摩阻力分布沿深度先增大后减小.桩基侧摩阻力与桩侧土压力和桩土相对位移的发挥程度均相关.在深度小于1.10D的范围内,桩侧土压力线性增大,桩土相对位移逐渐减小,但桩侧土压力增大对侧摩阻力的影响比较显著,故侧摩阻力表现为逐渐增大;在1.10D~1.40D深度范围内,由于桩体的压缩导致桩土相对位移继续减小,虽然在此深度范围内桩侧土压力也线性增大,但其对侧摩阻力的影响程度小于桩土相对位移减小所产生的影响程度,故侧摩阻力在该范围内逐渐减小.图 12b所示为开挖前后桩基轴力沿深度分布图,可以看出隧道开挖后桩基轴力都有所减小.

图 12 桩基侧摩阻力和轴力分布 Fig.12 Distribution of lateral frictional resistance fore on pile and axial force in pile
3.4 桩基承载力损失的影响因素分析

盾构隧道桩底施工时造成桩基承载力损失的因素主要有土层力学性质、盾构施工控制水平、桩基初始荷载水平和桩底与隧道顶部距离.土层力学性质直接影响土层提供侧摩阻力和桩端反力的能力,反应到图 9桩基Q~s曲线上为直接影响OABOACD曲线特征.盾构施工控制水平是对整个盾构施工过程的笼统表述,学术界用地层体积损失率来整体表达盾构施工控制水平,在图 9桩基Q~s曲线上为AC段特征,当体积损失率越大时,AC段造成的沉降就越大.桩基初始荷载水平为桩基承载能力已经发挥的程度,在图 9桩基Q~s曲线上为反应A点在OB段中所处的位置.桩底与隧道顶部距离越近,隧道开挖卸载造成的桩基沉降越大,在图 9桩基Q~s曲线上主要反应在AC段.

由于土层性质过于复杂,桩基所处的持力层、侧摩阻力层和盾构穿越断面地层等不同,对承载力影响的结果都有很大不同,因此对于地层性质的影响常作具体问题具体分析.工程中盾构施工控制水平是重中之重;桩基的初始荷载水平除了影响承载力损失以外,也用于衡量是否具有足够的储备承载力来承受盾构穿越造成的承载力损失;桩底与隧道顶部距离是工程设计时需要解决的关键问题.因此,在上述实例的基础上,对盾构施工控制水平、桩基初始荷载水平和桩底与隧道顶部距离3个因素做规律性分析.

(1) 盾构施工控制水平.用体积损失率来表示盾构施工控制水平,通过不同的体积损失率来计算盾构施工控制水平对单桩承载性能的影响.计算结果如图 13所示,盾构施工控制水平对桩基的承载性能影响显著,承载力损失值随着体积损失率的增大而增大.当体积损失率增长到1.5%时,桩基的承载力损失值增加到40%.因此,在盾构下穿桩基时应努力提高盾构控制水平.

图 13 不同体积损失率下的承载力损失值 Fig.13 Loss of pile capacity at different loss rate

(2) 初始荷载水平.为了表示桩基在初始情况下的受荷状态,用式(12)来表达桩基的初始荷载水平ζQ0为初始桩顶荷载,Q1为隧道开挖前桩基的极限承载力.分别进行ζ为10%、20%、30%、40%、50%的情况下隧道开挖时桩基承载性能影响计算.计算结果如图 14所示,桩基的初始荷载水平越大,开挖引起的桩基承载力损失值也越大,即桩基在初始情况下的承载力发挥得就越多,隧道开挖对桩基的承载力影响也更为显著.但是,桩基初始荷载从10%提高到50%时,承载力损失量只是提高了5%.因此,对于初始荷载水平的敏感性要低于盾构施工控制水平,降低体积损失率能够更加有效地降低承载力损失.

$ \zeta = \frac{{{Q_0}}}{{{Q_1}}} \times 100\% $ (12)
图 14 不同初始承载水平下的承载力损失值 Fig.14 Loss of pile capacity at different initial load level

(3) 桩底与隧道顶部距离.h为桩底与隧道顶部距离,进行h为1D~6D(增长梯度为0.5D)的隧道开挖计算.除h发生改变以外,其他参数保持不变,体积损失率为0.5%,桩及土体的尺寸、桩土的物理力学参数设置同上.图 15为桩基承载力损失值与不同h/D的关系图.由图可知,承载力损失值随着h/D的增大而曲线减小.当h/D增加到6时,隧道开挖引起的桩基承载力损失值减小到近10%.据此可认为,当h/D≥6以后,隧道开挖对桩基的影响已较小.

图 15 不同桩底与隧道顶部距离下的承载力损失值 Fig.15 Loss of pile capacity at different distance between pile bottom and the roof of the tunnel
4 结论

(1) 提出一种盾构桩底施工桩基承载能力损失计算方法.先实现考虑卸载过程的双曲线荷载传递函数编译,然后利用该桩土接触模型得到隧道开挖前和开挖后的Q~s曲线;取s=50 mm时对应的荷载为桩基极限承载能力,采用桩基极限承载力损失百分比作为隧道开挖对桩基承载性能影响的评价指标.并结合实例,分析了桩基承载曲线变化特征、桩体内力变化规律和承载力损失影响因素.

(2) 在实例中,隧道开挖前桩基的极限承载力Q1为508.0 kN.盾构施工体积损失率控制在0.5%时,隧道开挖后的桩基的极限承载力Q2为400.3 kN.隧道开挖导致的桩基承载力损失值为22%.初始工作荷载Q0为254.0 kN,因此,盾构施工保持该控制水平时降低了部分桩基承载能力,但不影响桩基础正常使用.隧道开挖后,桩体中存在一点侧摩阻力不变,在该点上部桩体侧摩阻力增大,在该点下部桩体侧摩阻力有所减小.

(3) 盾构隧道桩底施工时桩基承载力损失的影响因素主要有土层力学性质、盾构施工控制水平、桩基初始荷载水平和桩底与隧道顶部距离.承载力损失值随着体积损失率的增大而增大;桩基的初始荷载水平越大,承载力损失值也越大;桩底与隧道顶部的距离越大,桩基承载力损失值越小;承载力损失对盾构施工控制水平更加敏感,降低体积损失率能够更加有效地降低承载力损失.

参考文献
[1]
MROUEH H, SHAHROUR I. Three-dimensional finite element analysis of the interaction between tunneling and pile foundations[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2007, 26(3): 217
[2]
LEE G T K, NG C W W. Effects of advancing open face tunneling on an existing loaded pile[J]. Journal of Geotechnical & Geoenvironmental Engineering, 2005, 131(2): 193
[3]
朱逢斌, 杨平, ONG C W. 盾构隧道开挖对邻近桩基影响数值分析[J]. 岩土工程学报, 2008, 30(2): 298
ZHU Fengbi, YANG Ping, ONG C W. Numerical analysis on influence of shield tunnel excavation to neighboring pile[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(2): 298 DOI:10.3321/j.issn:1000-4548.2008.02.025
[4]
霍军帅, 王炳龙, 周顺华. 地铁盾构隧道下穿城际铁路地基加固方案安全性分析[J]. 中国铁道科学, 2011, 32(5): 71
HUO Junshuai, WANG Binglong, ZHOU Shunhua. Safety analysis of foundation reinforcement scheme for shield tunnel under-passing intercity railway[J]. Chinese Railway Science, 2011, 32(5): 71
[5]
王卫东, 李永辉, 吴江斌. 超长灌注桩桩土界面剪切模型及其有限元模拟[J]. 岩土力学, 2012, 33(12): 3818
WANG Weidong, LI Yonghui, WU Jiangbin. Pile-soil interface shear model of super long bored pile and its FEM simulation[J]. Rock and Soil Mechanics, 2012, 33(12): 3818
[6]
袁海平, 王斌, 朱大勇, 等. 盾构近距侧穿高架桥桩的施工力学行为研究[J]. 岩石力学与工程学报, 2014, 33(7): 1457
YUAN Haiping, WANG Bin, ZHU Dayong, et al. Mechanical behaviours of a shield tunnel adjacent to existing viaduct pile foundations[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(7): 1
[7]
陈明, 李镜培. 开挖卸荷桩土界面荷载传递模型的修正与验证[J]. 哈尔滨工业大学学报, 2014, 46(10): 75
CHEN Ming, LI Jingpei. Improved hyperbolic model of load-transfer for pile-soil interface and its verification considering deep excavation[J]. Journal of Harbin Institute of Technology, 2014, 46(10): 75 DOI:10.11918/j.issn.0367-6234.2014.10.013
[8]
李新星, 杨志豪. 盾构近距离穿越高架桩基的施工影响与保护措施[J]. 岩土力学, 2015, 36(S1): 537
LI Xinxing, YANG Zhihao. Influences of construction of side-crossing shield tunnel on adjacent pile foundation and reinforcement effect of protection measures[J]. Rock and Soil Mechanics, 2015, 36(S1): 537
[9]
ZHANG R J, ZHENG J J, YU S. Responses of piles subjected to excavation-induced vertical soil movement considering unloading effect and interfacial slip characteristics[J]. Tunnelling and Underground Space Technology, 2013, 36(2): 66
[10]
CLOUGH G W, DUNCAN J M. Finite element analyses of retaining wall behavior[J]. Asce Soil Mechanics Foundation Division Journal, 1971, 97(12): 1657
[11]
KISHIDA H, UESUGI M. Tests of the interface between sand and steel in the simple shear apparatus[J]. Journal of the Physical Society of Japan, 1987, 7(1): 102
[12]
费康, 张建伟. ABAQUS在岩土工程中的应用[M]. 北京: 中国水利水电出版社, 2013
FEI Kang, ZHANG Jianwei. Application of ABAQUS in geotechnical engineering[M]. Beijing: China Water and Power Press, 2013
[13]
NG C W W, LU H, PENG S Y. Three- dimensional centrifuge modeling of the effects of twin tunneling on an existing pile[J]. Tunneling and Underground Space Technology Research, 2013, 35(4): 189