Lateral bearing characteristics of inclined pile-cap system installed in sandy ground
-
摘要: 通过开展室内模型试验,研究了砂土中斜桩-承台系统的水平承载特性,重点分析了桩身倾角及桩身抗弯刚度两个因素对桩顶水平位移、桩身弯矩分布的影响。利用有限元方法开展了与模型试验相对应的三维有限元分析,两者所得结果的对比表明所采用的有限元分析方法能较好地反映斜桩在水平荷载下的受力和变形特性。为了更好地与工程实际相对应,基于上述数值分析方法开展了原型尺寸的数值参数分析,考虑了桩体的弹塑性力学特性,探究土体弹性模量、桩身刚度以及桩身长度对斜桩-承台系统水平承载能力的影响。计算结果表明,相较于改变桩身的抗弯刚度,土体弹性模量的变化对斜桩-承台系统水平承载能力的影响更为显著;当桩长小于25 m(25倍桩径)时,桩长的增加能较明显降低桩顶水平位移响应,当桩长超过该值时,斜桩水平承载能力的提升效率下降。此外,基于试验数据开展了多元回归分析,获得了斜桩-承台系统桩顶水平位移的半经验预测公式,该公式的预测值与试验实测值及数值参数分析计算值的吻合度均较高,具有良好的适用性,有望能为砂土中斜桩工程的设计和分析提供一定的参考依据。Abstract: The lateral bearing characteristics of inclined pile-cap systems installed in sandy soil are comprehensibly studied via a series of laboratory model tests, and the influences of inclination angle and flexural rigidity of pile on the lateral displacement at pile top and the distribution of pile bending moment are analyzed. Based on the model tests, a suite of three-dimensional finite element analyses are performed, and the favorable comparison of the test and numerical results convinces the capability of the adopted finite element analysis procedure in replicating the stress and deformation characteristics of the inclined piles in the model tests. Subsequently, a number of numerical parametric analyses based on the prototype models are implemented, accounting for the elastic-plastic mechanical properties of the pile, as well as the variations in elastic modulus of soil, flexural rigidity and length of the pile. The analysis results show that the variation in the elastic modulus of soil has a more significant effect on the lateral bearing capacity of the inclined pile-cap system than the change in pile flexural rigidity. When the pile length is less than 25 times the pile diameter(D), the increase of pile length can effectively reduce the lateral displacement response at the pile top. Otherwise, increasing the pile length has negligible effects on enhancing the lateral performance of the inclined pile-cap system. Furthermore, the multivariate regression analysis is performed based on the test data of lateral displacement at pile top, from which a semi-empirical equation can be derived. The performance of this semi-empirical equation is favorably examined by both the relevant experimental data and the numerical parametric analysis results, which likely provides a useful reference for the design and analysis of the inclined pile-cap system in sandy ground.
-
0. 引言
扩体预制桩由内芯预制桩和外围扩体材料固结体组成,一般可通过就地搅拌水泥土或机械扩孔并向孔内灌注扩体材料浆体,再植入内芯同轴预制桩的工艺成桩[1]。近年来,随着厂拌高质量扩体浆材的利用和成桩工艺的创新[2],扩体预制桩的桩身质量和承载性能大幅提高,开始被广泛应用于高层建筑、桩板路桥等承载力和沉降控制要求较高的建构筑物桩基础。
与单一材料桩相比,扩体预制桩竖向承载时涉及到内芯预制桩、外围扩体和周围土3种不同材质的相互作用,荷载传递机制十分复杂[3-4]。然而,目前有关扩体预制桩承载机制的理论分析研究极为有限。考虑芯桩-扩体-周围土协同工作对承载性能的影响,Wang等[5]在假定芯桩与外围水泥土界面固定黏结的基础上,基于荷载传递法建立了劲性搅拌桩荷载传递简化计算模型;俞建霖等[6]针对刚性基础下砼芯水泥土桩复合地基,给出了砼芯、外围水泥土、周围土三者的受力平衡微分方程;结合引扩孔灌浆扩体预制桩工艺特点,张浩等[7]在考虑芯桩−扩体材料−周围土径向相互作用的基础上,提出了扩体预制桩整体复合模量表达式,并给出了全截面竖向受压模式下扩体桩的承载变形解析计算方法;考虑芯桩-水泥土-周围土界面的非线性力学特征,Zhu等[8]通过三折线模型建立了刚性基础下劲性搅拌桩的荷载传递理论分析模型,极大完善了扩体预制桩荷载传递机制的分析理论。然而,这些研究多集中于桩顶全截面受压工况,即桩顶处芯桩与扩体近似等沉变形受荷,较适用于建筑筏板、桥台等刚性基础下的桩基工况。
近年来,随着扩体预制桩在桩板式无土路基结构、高桩承台结构中的广泛应用,芯桩单独承载的扩体预制桩荷载传递机制研究不容忽视。既有试验研究表明[9],芯桩单独承压下扩体预制桩中芯桩-扩体-土的荷载传递机制与全截面受压下显著不同。顾士坦等[10]基于Mindlin弹性解和剪滞原理,仅考虑芯桩-外围水泥土的相互作用,给出了芯桩单独受力时劲性搅拌桩芯桩荷载传递规律的理论分析方法;任连伟等[11]、Zhou等[12]基于荷载传递法,分别给出了高喷插芯组合桩和静钻根植桩芯桩受力时荷载传递简化计算方法,但芯桩-水泥土界面完全线弹性剪切模型的假定,无法全面反映此类扩体预制桩的工作状态。
据此,本文针对芯桩单独承载工况下全长扩体预制桩中芯桩-扩体-周围土的荷载传递特点,在充分考虑芯桩-扩体、扩体-土界面非线性相互作用特性以及两界面剪切性能发挥交互影响的基础上,构建了全长扩体预制桩芯桩承载荷载传递计算分析模型,并参考工程实例,对全长扩体预制桩芯桩承载荷载传递机制及其影响因素进行了分析。
1. 荷载传递机制与模型假定
1.1 工作机理
图 1为典型全长扩体预制桩芯桩单独承载工况下桩土作用示意图。荷载P直接作用于内芯预制桩顶部,芯桩首先向下发生相对外围扩体材料的变形,并在芯桩-扩体界面产生剪切作用;外围扩体材料发生剪切变形,进一步将荷载传递至周围土。同时随着荷载向下传递,桩端处发生向下刺入变形趋势,并受到地基土的反力作用。因此芯桩承载工况下上部荷载是经由芯桩-扩体-周围土相互作用传递至周围土的。芯桩承载扩体预制桩的承载变形机制涉及到芯桩-扩体、扩体-周围土内外两个相互作用界面的协同工作。扩体预制桩整体受到周围土侧阻力τs、芯桩端阻力Rp和扩体端阻力Rc,内部则存在芯桩与扩体材料的剪切作用力τc。
由图 1扩体预制桩的承载变形特征,可得芯桩竖向变形s(z)的表达式:
s(z)=spc(z)+sc(z)+scs(z)+ss(z)=sp(z)+ss(z)。 (1) 式中:spc(z),scs(z)为深度z处芯桩-扩体界面和扩体-周围土界面的相对位移;sc(z),ss(z)为z处扩体材料和周围土的剪切变形;sp(z)为芯桩与土的相对位移。
1.2 模型简化及基本假定
既有研究表明[8, 13],理想弹塑性模型能较好的反映水泥土、水泥砂浆等水泥基胶凝材料与土之间的荷载传递关系,因此将扩体-周围土界面模型简化为理想弹塑性模型,如图 2(a)所示,即
τs={ksscs (scs<sfs)τfs(scs⩾sfs)。 (2) 式中:τs,scs为扩体-土界面侧阻力和相对位移;ks为扩体-土剪切刚度;τfs,sfs界面极限侧阻力和极限相对位移。
芯桩与水泥土等水泥基胶凝材料的剪切试验[14]表明,在较小相对位移下界面剪切应力即表现出较高水平,随着相对位移的增大,界面应力达到极限值后会出现迅速降低,界面呈脆性破坏模式,如图 2(b)所示。因此,芯桩-扩体界面不宜采用完全线弹性剪切模型,可简化为图 2(b)理想弹性-破坏模型,即
τc={kcspc(spc<sfc)τ′fc (spc⩾sfc)。 (3) 式中:τc,spc为芯桩-扩体界面侧阻力和相对位移;kc为芯桩-扩体界面剪切刚度,由试验得τc-spc曲线的直线段斜率确定;τfc,sfc为界面极限侧阻力和相对位移;τ′fc为界面残余侧阻力,与扩体材料水泥掺量、龄期等有关,约为界面极限侧阻力的(0.057~0.153)倍[14]。
据此,分析时作如下基本假定:①内芯预制桩和外围扩体材料均为线弹性体;②针对目前以水泥土、水泥砂浆为扩体材料的扩体预制桩,扩体剪切刚度远大于周围土,且厚度相对有限[1-2],受荷时扩体剪切变形较小,忽略不计;③假定芯桩-扩体荷载传递为理想弹性-破坏模型,扩体-土荷载传递为理想弹塑性模型。
2. 荷载传递方程与计算
2.1 计算单元分析
如图 3所示,取深度z处“扩体预制桩-土”单元进行受力分析,根据单元体静力平衡条件可得
dPp(z)=−τc(z)up dz, dPc(z)=τc(z)up dz−τs(z)uc dz∘} (4) 式中:Pp(z),Pc(z)为z处芯桩和扩体的截面轴力;τc(z),τs(z)为z处芯桩-扩体、扩体-土界面的侧阻力;up,uc为芯桩、扩体外侧周长。
若受力变形后芯桩单元体和扩体单元体的压缩量分别为Δp和Δc,根据图 3(b)变形模式可得
dspc(z)=Δc−Δp,dscs(z)=−Δc。} (5) 进而,基于胡克定律由式(4),(5)联合,可得芯桩-扩体、扩体-周围土的界面微分控制方程:
d2spc(z)dz2=(upEcAc+upEpAp)τc(z)−ucEcAcτs(z),d2scs(z)dz2=ucEcAcτs(z)−upEcAcτc(z)} (6) 式中:Ep,Ec为芯桩、扩体的弹性模量;Ap,Ac为芯桩、扩体的横截面积。
2.2 荷载传递模式
芯桩承载扩体预制桩与土相互作用模式见图 4。
随着桩顶荷载增加,扩体-土界面和芯桩-扩体界面将由完全弹性状态、从浅而深逐渐进入塑性或破坏状态。设桩长L,扩体-土界面塑性区深度L1,芯桩-扩体界面破坏区深度L2。一般情况下,由于芯桩-扩体界面剪切强度比扩体-土界面要大,工程实践中为充分发挥芯桩的承载性能,往往会通过调配扩体材料强度和厚度来保证芯桩-扩体共同作用,避免芯桩-扩体界面提前破坏[1],此时L1往往大于L2,如图 4(a)。
当外围扩体强度较小或厚度较薄等特殊情况下,芯桩单独承载桩顶附近芯桩和扩体应力差异显著,芯桩-扩体界面可能先出现破坏,而后扩体-土界面逐渐进入塑性状态,桩-土相互作用为图 4(b)所示模式。
2.3 方程建立与计算
(1)L2≤L1工况
根据芯桩-扩体、扩体-土两界面相互作用状态的不同,将扩体预制桩划分为3个区域,见图 4(a)。
a)当L1 < z≤L时,芯桩-扩体界面、扩体-土界面均处于弹性阶段,由式(2),(3),(6)可得
d2spc(z)dz2=αspc(z)+βscs(z),d2scs(z)dz2=γspc(z)+δscs(z)。} (7) 式中:α=upkc/EcAc+upkc/EpAp;β=−ucks/EcAc;γ=−upkc/EcAc;δ=ucks/EcAc。
求解方程,可得芯桩-扩体的相对位移spc1(z)、扩体-周围土的相对位移scs1(z):
spcl(z)=1γ(C1a21ea1z+C2a22ea2z+C3a23ea3z+C4a24ea4z)−δγ(C1ea1z+C2ea2z+C3ea3z+C4ea4z),scs1(z)=C1ea1z+C2ea2z+C3ea3z+C4ea4z} (8) 式中:C1,C2,C3,C4为方程待定系数;
a1=√[(δ+α)+√(δ+α)2−4(δα−γβ)]/2; a2=−√[(δ+α)+√(δ+α)2−4(δα−γβ)]/2; a3=√[(δ+α)−√(δ+α)2−4(δα−γβ)]/2; a4=−√[(δ+α)−√(δ+α)2−4(δα−γβ)]/2。 由式(5),(8)可得芯桩轴力和扩体轴力:
Pp1(z)=EpAp(δ−γγ)(C1a1ea1z+C2a2ea2z+C3a3ea3z+C4a4ea4z)−EpApγ(C1a13ea1z+C2a23ea2z+C3a33ea3z+C4a43ea4z),Pc1(z)=−EcAc(C1a1ea1z+C2a2ea2z+C3a3ea3z+C4a4ea4z)。} (9) 此时,由式(1)可得芯桩竖向变形s1(z):
s1(z)=sp1(z)+ss1(z) = spc1(z)+scs1(z)+ss1(z)。 (10) 式中:ss1(z)=τs(z)rcln(rm/rc)/Gs;rc为扩体外壁半径;Gs为桩周土剪切模量;rm为影响半径[15]。
b)当L2 < z≤L1时,扩体-周围土界面处于塑性状态,芯桩-扩体界面处于弹性阶段,
d2spc(z)dz2=αspc(z)+κ,d2scs(z)dz2=γspc(z)−κ。} (11) 式中:κ=−ucτfs/EcAc。
求解方程,可得z处芯桩-扩体相对位移spc2(z)、扩体-周围土相对位移scs2(z):
spc2(z)=D1e√αz+D2e−√αz−κα,scs2(z)=γD1αe√αz+γD2αe−√αz−κ2(γα+1)z2+D3z+D4。} (12) 式中:D1,D2,D3,D4为方程待定系数。
进而,芯桩轴力Pp2(z)和扩体轴力Pc2(z)为
Pp2(z)=−EpAp[γ+α√α(D1e√αz−D2e−√αz)−κ(γα+1)z+D3],Pc2(z)=−EcAc[γD1√αe√αz−γD2√αe−√αz−κ(γα+1)z+D3]。} (13) 此时,芯桩竖向变形s2(z)为
s2(z)=sp2(z)+ss2(z) = spc2(z)+scs2(z)+ss2(z)。 (14) 式中:ss2(z)=τsfrcln(rm−rc)/τsfrcln(rm−rc)GsGs。
c)当0≤z≤L2时,扩体-周围土界面处于塑性状态,芯桩-扩体界面处于破坏阶段,
d2spc(z)dz2=g, d2scs(z)dz2=f∘} (15) 式中:g=(upEcAc+upEpAp)τ′fc−ucτfsEcAc;f=ucτfsEcAc−upτ′fcEcAc。
求解方程,可得z处芯桩-扩体相对位移spc3(z)、扩体-周围土相对位移scs3(z):
spc3(z)=12gz2+F1z+F2,scs3(z)=12fz2+F3z+F4。} (16) 式中:F1,F2,F3,F4为方程待定系数。
进而,芯桩轴力Pp3(z)和扩体轴力Pc3(z)为
Pp3(z)=−EpAp(gz+fz+F1+F3),Pc3(z)=−EcAc(fz+F3)。} (17) 此时,芯桩竖向变形s3(z):
s3(z)=spc3(z)+scs3(z)+ss3(z)。 (18) 式中:ss2(z)=τsfrcln(rm/rc)/Gs。
d)边界与连续条件
根据扩体桩受力特点,桩顶边界条件为
Pp(z)|z=0=P,Pc(z)|z=0=0。} (19) 桩端边界条件为
Pp(z)|z=L=kLApsp(z)|z=L,Pc(z)|z=L=kLAcscs(z)|z=L。} (20) 式中:kL为桩端土刚度系数,基于Randolph等[16]刚性体压入弹性半空间解答,kL=4GL/4GL[πrc(1−ν)][πrc(1−ν)],GL为桩端土体剪切模量;ν为桩端土泊松比。
桩身受力变形连续条件:
spc2(z)|z=L1−=spc1(z)|z=L1+,scs2(z)|z=L1−=scs1(z)|z=L1+,Pp2(z)|z=L1−=Pp1(z)|z=L1+,Pc2(z)|z=L1−=Pc1(z)|z=L1+,scs2(z)|z=L1=γD1αe√αL1+γD2αe−√αL1− κ2(γα+1)L12+D3L1+D4=sfs。} (21) spc3(z)|z=L2−=spc2(z)|z=L2+,scs3(z)|z=L2−=scs2(z)|z=L2+,Pp3(z)|z=L2−=Pp2(z)|z=L2+,Pc3(z)|z=L2−=Pc2(z)|z=L2+,spc3(z)|z=L2=12gL22+F1L2+F2=sfc。} (22) 综上,针对L2≤L1工况,根据边界和连续条件(19)~(22),可得关于C1~C4、D1~D4、F1~F4、L1、L2的14个方程,求解未知量并回代入式(8)~(10)、(12)~(14)、(16)~(18),即可得扩体桩桩身轴力(Pp(z),Pc(z))与位移(spc(z),scs(z),s(z))。
考虑不同计算条件影响:①当桩顶荷载相对较小时,全桩长范围芯桩-扩体、扩体-土界面均处于弹性状态时(L1=0,L2=0),上述计算可退化为关于C1~C4共4个未知量的求解问题,仅需将式(8)、(9)代入边界条件(19)、(20)即可得出4个求解方程,进而可得该工况下桩身受力响应;②若随着荷载增加,扩体-土界面开始部分出现塑性区、但芯桩-扩体界面全桩长处于弹性阶段时(0 < L1 < L,L2=0),上述计算退化为关于C1~C4、D1~D4、L1共9个未知量的求解问题,此时将式(8),(9),(12),(13)代入边界和连续条件(19)~(21)即可得出9个求解方程,进而可得该工况下桩身受力响应;③随着桩顶荷载的增加,芯桩-扩体界面也开始出现脆性破坏(0 < L1 < L,0 < L2 < L1),即进入式(7)~(18)计算工况;④当扩体-土界面全部进入塑性区而芯桩-扩体界面局部出现破坏时(L1=L,0 < L2 < L1),上述计算退化为关于D1~D4、F1~F4、L2共9个未知量的求解问题,此时将式(12),(13),(16),(17)代入边界和连续条件(19)~(21)即可得出9个求解方程;⑤当扩体-土界面全部进入塑性区且芯桩-扩体界面全部破坏时(L1= L2=L),上述计算退化为关于F1~F4共4个未知量的求解问题,将式(16),(17)代入边界条件(19),(20)即可得出4个求解方程;⑥若实际工况是非均质层状土,可通过桩长范围内各土层参数的加权平均进行简化计算;当然也可通过桩身离散的矩阵系数传递法,采用界面相对位移迭代计算的方法近似求解。
(2)L2 > L1工况
如图 4(b)所示,当芯桩-扩体破坏区深度L2大于外部扩体-土塑性区深度L1时,沿桩长同样可以将扩体预制桩划分为3个区域:①L2 < z≤L时,芯桩-扩体界面、扩体-土界面均处于弹性阶段;②L1 < z≤L2时,芯桩-扩体界面处于破坏状态,而外部扩体-土界面仍处于弹性状态;③0 < z≤L1时,扩体-土界面处于塑性状态,芯桩-扩体界面处于破坏阶段。
其中,在L2 < z≤L和0 < z≤L1区域,桩-土相互作用状态与前述2.3节中扩体桩上、下两段区域桩-土相互作用状态相同,芯桩-扩体、扩体-土两界面的微分控制方程与式(7),(15)相同,求解方法也相同,在此不再赘述。
当L1 < z≤L2时,由式(2),(3),(6)可得此时芯桩-扩体和扩体-土界面微分方程:
d2spc(z)/dz2=ξ,d2scs(z)/dz2=δscs(z)−ζ。} (23) 式中:ξ=upτ′fc/EpAp;ζ=upτ′fc/EcAc。
求解方程,可得该区域z处芯桩位移s′pc2(z)、扩体-周围土相对位移s′cs2(z)。
s′pc2(z)=12ξz2+J1z+J2z,s′cs2(z)=J3e√ζz+J4e−√ζz+ζ/ξ。} (24) 式中:J1,J2,J3,J4为方程待定系数。
进而,芯桩轴力P′p2(z)、扩体轴力P′c2(z)和芯桩竖向变形s′2(z):
P′p2(z)=−EpAp[ξz+J1],P′c2(z)=−EcAc[J3e√ζz+J4e−√ζz],} (25) s′2(z)=s′pc2(z)+s′cs2(z)+s′s2(z)。 (26) 式中:s′2(z)=τs(z)rcln(rm/rc)/Gs。
据此,可对L2 > L1工况下扩体预制桩荷载传递进行计算分析。其中,桩顶和桩端边界条件同式(19),(20);桩身受力变形连续条件为
s′pc2(z)|z=L2−=s′pc1(z)|z=L2+,s′cs2(z)|z=L2−=s′cs1(z)|z=L2+,P′p2(z)|z=L2−=P′p1(z)|z=L2+,P′c2(z)|z=L2−=P′c1(z)|z=L1+,s′pc2(z)|z=L2=sfc,} (27) s′pc3(z)|z=L1−=s′pc2(z)|z=L1+,s′cs3(z)|z=L1−=s′cs2(z)|z=L1+,P′p3(z)|z=L1−=P′p2(z)|z=L1+,P′c3(z)|z=L1−=P′c3(z)|z=L1+,s′cs3(z)|z=L1=sfs。} (28) 式中:s′pc1(z),s′cs2(z),P′p1,P′c1(z)分别为L2 < z≤L芯桩和扩体的变形和轴力;s′pc2(z),s′cs2(z),P′p2(z),P′c2(z)分别为L1 < z≤L2芯桩和扩体的变形和轴力;s′pc3(z),s′cs3(z),P′p3(z),P′c3(z)分别为0 < z≤L1芯桩和扩体的变形和轴力。
考虑不同计算条件,采用与2.3节相同的求解方法,可得L2 > L1工况下扩体桩桩身轴力与位移。
3. 实例验证与荷载传递规律分析
3.1 水泥土扩体预制桩实例
水泥土扩体预制桩现场试验参见文献[17]。试桩桩长10 m,芯桩为PHC-400AB(95)管桩,外围为厚300 mm水泥土桩,水泥掺量18%,水泥土90 d抗压强度≥1.6 MPa。计算中,C80管桩弹性模量38 GPa,水泥土弹性模量800 MPa。针对界面荷载传递模型参数取值,根据文献[14]混凝土-水泥土界面极限侧阻力约抗压强度的0.064~0.259,临界相对位移约1.23~2.40 mm,界面残余侧阻力约为极限侧阻力的0.057~0.153,取芯桩-水泥土界面极限侧阻力为250 kPa,临界相对位移sfc=2 mm,残余侧阻力τ′fc=25 kPa。基于剪切位移法[15],以芯桩-水泥土极限剪应力估算外围水泥土剪切变形约为1.4×10-4 m。可见,虽然外围水泥土包裹厚度达300 mm,其剪切变形仍较小,可忽略不计。水泥土-周围土界面荷载传递参数取值,以文献[13]水泥土桩-土试验结果作为参考,取极限侧阻力τfs=65 kPa,ks=20 kPa/mm。依据桩底土体的性质,结合文献[17]现场载荷试验桩端阻力和桩端沉降量,反演计算取桩端土刚度系数kL为1000 kPa/mm。
采用本文方法对现场试桩进行计算,如图 5,6。各级荷载作用下,桩顶沉降计算值大于实测值,芯桩轴力计算值整体上与实测值相比偏小,这是由于实际芯桩后植入过程中会对水泥土和周围土产生挤密作用,使界面黏结强度增大[2],桩体整体承载性能提升,但计算分析中界面参数未考虑该部分的影响。
从芯桩轴力分布规律来看,芯桩轴力沿深度逐渐减小,且随着上部荷载增加,该变化趋势越来越明显。各级荷载下桩端轴力约占桩顶总荷载18.4%~21.3%,说明水泥土扩体预制桩整体上表现为端承摩擦桩的承载特性。进一步对比芯桩轴力计算结果与现场实测结果,两者在桩顶、桩端附近差异性较小,这主要因为理论计算中桩端刚度由现场实测反演所得,且桩顶荷载边界与现场试桩相同。在桩身中部范围芯桩轴力存在差异,且实测桩身轴力沿桩身分布有波动性,这是由于实际桩周水泥土就地搅拌时均匀性不能得到很好保证,从而影响内部芯桩受力。参考文献[17]对试桩沉降压缩量进行计算分析,如图 6(b)和表 1所示,可以看出本文方法也能较好反映桩身压缩变形特征。
表 1 试桩沉降压缩量对比分析Table 1. Comparative analysis of settlement compression of pile计算分析 桩顶荷
载/kN桩顶
沉降/mm桩身压
缩/mm桩端沉
降/mm桩身压缩
沉降比/%试桩[17] 2200 10.23 2.82 7.41 27.57 本文方法 2200 9.83 2.73 7.10 27.76 3.2 水泥砂浆扩体预制桩实例
德上高速合肥-枞阳段地处安徽中南部,部分路段采用了“桩柱一体式”桩板路基结构。考虑沿线基岩较浅且上覆沉积物工程性能较好,传统打入式预制桩沉桩困难,因此采用了引扩孔灌浆植桩工艺。芯桩采用PHC-500AB(100)管桩,外围采用厂拌充盈水泥砂浆,28 d抗压强度 > 15 MPa,桩长10 m。为揭示水泥砂浆扩体预制桩的工程特性,开展扩孔直径分别为700,800,1000,1200 mm(对应扩体厚度100,150,250,350 mm)的静载试验,试验区土层性质见表 2。
表 2 土层主要物理力学指标Table 2. Physico-mechanical indexes of soil layers土层 h/
mγ/
(kN·m-3)c/
kPaϕ/
(°)Es1-2/
MPafsk/
kPa②粉土 6.1 19.3 20.6 10.2 7.50 140 ②1粉质黏土 2.4 19.1 34.6 12.2 8.59 180 ③粗砂 2.8 — — — — 250 ⑦1全风化角砾岩 — — 51.9 12.5 — 350 为了获取芯桩-砂浆、砂浆-土界面荷载传递参数,采用自研界面剪切装置(图 7)开展界面剪切试验,试验方法同直接剪切试验,每组试验3个试样,分别考虑100,200,300 kPa 3种不同的法向应力,参考《土工试验方法标准》(GB/T50123—2019),每个法向力采用一次轻轻施加。受制于剪切试样尺寸(ϕ61.8 mm×20 mm)较小,芯桩混凝土采用同强度C80灌浆料替代,所用土样和水泥砂浆皆取自工程现场。
制作水泥砂浆-土试样时,为保证试样土重度与工程实际的一致性,首先参考地勘报告调整土样含水率,并以勘察报告土体重度为基准,根据下剪切盒体积,计算每个试样所需土的重量,将其填入下剪切盒至剪切界面高度,然后在其上剪切盒内浇入水泥砂浆。制作芯桩-水泥砂浆试样时,首先将替代灌浆料浇入模具中,待标准养护28 d后再在其上浇入水泥砂浆。制样完成后,所有试样标准养护14 d后开展界面剪切试验,部分成果如图 8所示。其中,散点曲线为试验测试结果,黑色实折线为试验数据拟合模型曲线。
依据界面剪切试验,取芯桩-砂浆界面极限侧阻力368 kPa,临界相对位移2.3 mm,残余侧阻力74 kPa;取砂浆-土界面极限摩阻力90 kPa,临界相对位移3 mm。计算中,内芯预制桩弹性模量38 GPa,扩体砂浆弹性模量7 GPa,参考《公路桥涵地基与基础设计规范》(JTG3363—2019)桩端岩石地基竖向抗力系数取值,桩端全风化角砾岩刚度系数kL取2500 kPa/mm。
图 9为不同扩径工况扩体预制桩荷载-位移曲线。随着扩径直径增加,扩体预制桩的承载性能大体呈增大变化趋势,当扩径直径由700 mm增加到1000 mm时,同级荷载桩顶沉降降低约44.3%~59.5%;而扩径的进一步增加(1200 mm),桩顶沉降变化量降低。
采用本文方法计算的荷载-位移曲线与实测结果具有较好的一致性,同级荷载作用下计算桩顶沉降相较实测值整体偏大。当扩径为700 mm时,各级荷载水平下桩顶沉降平均误差率6.7%;随着扩径直径增大,计算值误差率随之增大。当扩径直径由700 mm增加至800,1000,1200 mm时,桩顶沉降平均误差率分别增加至24.8%,34.6%,31.1%。分析其原因,上述计算分析未考虑芯桩后植入施工的扩挤影响,尤其是后植入芯桩对孔内砂浆挤密及砂浆渗滤增强桩-土界面剪切性能的作用效应,同时简化理论模型中忽略外围扩体剪切变形的基本假定也对计算结果有影响。考虑相关技术规程和工程实践中[1]为了充分发挥芯桩高承载力的特性和合理控制孔内灌浆后植入桩的挤土效应,推荐扩孔直径宜为芯桩直径的1.5倍,即扩体厚度多控制在100~200 mm,因此本文计算方法适用于常规扩体预制桩的工程计算分析。
3.3 荷载传递规律分析
基于上述实例,进一步计算分析复合桩身荷载传递机制。以扩径直径700 mm(内芯预制桩500 mm)水泥砂浆扩体预制桩为例,图 10所示为计算的桩身轴力分布曲线。不同荷载水平下芯桩轴力沿桩身呈逐渐减小的变化趋势,随着桩顶荷载水平的增大,在上部0~3 m内芯桩轴力减小幅度相对较大,3~8 m内芯桩轴力减小幅度变小且相对稳定,而桩端附近(8~10 m)芯桩轴力减小幅度再次变大。
与之对应,外围水泥砂浆截面轴力沿桩身呈逐渐增大变化趋势,在0~3 m内外围砂浆轴力快速增加,说明芯桩单独承载工况下,桩顶附近芯桩和外围砂浆的相互作用较强;随着深度增加,外围砂浆轴力增加幅度降低且相对稳定,与芯桩趋于共同作用;桩端(8~10 m)附近,由于芯桩与外围砂浆刚度的显著差异,以及桩端角砾岩较高桩端反力作用,芯桩相对外围砂浆“下刺入”过程中,两者相互作用增强,外围砂浆轴力呈明显增大的变化趋势,如图 10(b)所示。
图 11所示为水泥砂浆扩体预制桩内、外两界面阻力分布曲线。内侧芯桩-水泥砂浆界面阻力沿桩身呈两端大、中间小的分布规律。这主要是由于在桩顶、桩端附近芯桩和外围砂浆的刚度差异导致两者界面相对位移趋势明显,引起较强烈的界面相互作用。其中桩顶荷载4000 kN时,芯桩-水泥砂浆界面在桩顶0.2 m范围内出现剪切破坏,界面阻力降为残余侧阻力。
对于外侧水泥砂浆-土界面,界面阻力沿桩身整体呈逐渐减小分布。当桩顶荷载较大时,水泥砂浆-土界面由上而下逐渐进入塑性状态,界面阻力达到极限侧阻力。随着桩顶荷载增加,界面塑性区长度增大,桩顶荷载3000,3500,4000 kN时,对应水泥砂浆-土界面塑性区长度分别为3.8,6.0,7.5 m。
进一步分析不同扩孔直径对扩体桩承载特性的影响,以桩顶荷载3000 kN工况为例,如图 12所示。随深度增加芯桩轴力逐渐减小,桩端处芯桩轴力占总截面轴力约56.1%,45.9%,33.3%,26.5%,扩体直径越大,桩端处芯桩轴力占比越小;砂浆截面轴力随深度增加而逐渐增大,占总截面轴力比例也逐渐增大,桩端处砂浆截面轴力占比约43.9%,54.1%,66.7%,73.5%。外围水泥砂浆对芯桩具有“扩径”作用。
当扩体直径相对较大(1200 mm)时,在桩顶(1 m)附近芯桩-砂浆的荷载传递作用出现明显降低现象,这主要是因为,随着扩体截面增大,其单位面积应力水平降低,单独承载芯桩与外围砂浆的相对位移势能增强,两者共同作用逐渐减低。由此可见,虽然扩体直径增加可以提升复合桩身的整体承载性能,但较厚外围扩体也可能制约内部芯桩承载能力的发挥。工程实践中宜注意芯桩和外围扩体的工作协同性。针对水泥砂浆扩体预制桩,以工程现场试验和本计算算例分析来看,对于直径500 mm的芯桩,外围砂浆扩体直径不宜大于1000 mm。
考虑外围扩体材料性能的影响,计算对比扩体弹性模量为0.8,4,7,10 GPa时扩体预制桩桩身内力,如图 13。随着扩体材料模量增大,芯桩轴力逐渐减小而外围扩体截面轴力逐渐增大。当扩体模量由较低的0.8 GPa(如水泥土),提升至7~10 GPa(如水泥砂浆)时,芯桩轴力由总截面轴力的93%,降低至77%~80%;而外围扩体截面轴力由总截面轴力的6%,提升至19%~23%,说明外围扩体材料刚度的改善,可提升内芯预制桩向周围传递荷载的效率。
由扩体截面轴力分布曲线可以看出,当扩体采用刚度相对较低(0.8 GPa)的水泥土时,芯桩与扩体的显著相互作用仅发生在桩顶(1 m)和桩端(2 m)附近有限区域内;随着扩体刚度增加,芯桩与扩体的显著相互作用区域逐渐增大,两者共同作用效应增强。在工程实践中可通过改善扩体材料性能,增强内芯预制桩与外围扩体的共同作用。
4. 结论
(1)根据芯桩承载全长扩体预制桩的受力变形特点,综合考虑芯桩-扩体、扩体-周围土内外两界面剪切特性的发挥与耦合作用,推导建立了芯桩承载扩体预制桩的荷载传递理论计算方法,工程案例分析验证了该方法的可靠性。
(2)外围扩体可改善桩身整体承载性能,当扩体直径(芯桩直径500 mm)由700 mm增加至1000 mm时,桩顶沉降降低约44.3%~59.5%;外围扩体可将芯桩荷载有效传递至周围土体;竖向荷载作用下,随着深度增加,芯桩轴力逐渐减小,而外围扩体轴力逐渐增大,对芯桩具有“扩径”作用。
(3)扩体刚度的增加可提升芯桩荷载传递效率,当扩体弹性模量由0.8 GPa增加至7~10 GPa时,芯桩轴力占总截面轴力比例由93%降低至77%~80%,而扩体截面轴力占比由6%增加至19%~23%,芯桩与扩体共同作用效应增强。
-
表 1 模型桩的尺寸、力学参数表
Table 1 Dimensions and mechanical parameters of model piles
名称 弹性模量/GPa 密度/(kg·m-3) 外径/mm 内径/mm 有机玻璃 3.2 180 20 — 铝 71.0 2700 20 — 钢 210.0 7900 20 16 表 2 试验方案列表
Table 2 Test schemes
组号 材料 编号 桩身倾角/(°) 外径/mm 内径/mm M1 有机玻璃 P11 0 20 — P12 10 20 — P13 20 20 — P14 30 20 — P15 45 20 — M2 铝 P21 0 20 — P22 10 20 — P23 20 20 — P24 30 20 — P25 45 20 — M3 钢 P31 0 20 16 P32 10 20 16 P33 20 20 16 P34 30 20 16 P35 45 20 16 表 3 桩基-承台系统水平极限承载力
Table 3 Lateral ultimate bearing capacities of pile-cap system
材料 桩身倾角/(°) 水平极限承载力/N 建议水平极限承载力/N 位移梯度法 位移限制法 有机玻璃 0 45.5 96.5 45.5 10 71.0 122.0 71.0 20 96.5 173.0 96.5 30 122.0 — 122.0 45 173.0 — 173.0 铝 0 71 224 71 10 122 275 122 20 173 — 173 30 224 — 224 45 275 — 275 钢 0 173 275 173 10 224 377 224 20 377 479 377 30 428 530 428 45 479 — 479 表 4 桩身倾斜角度对斜桩桩身最大弯矩的影响
Table 4 Influence of inclination angle of pile on maximum bending moment
编号 倾角/(°) 最大弯矩/(N·m) 相对直桩的减幅/% P11 0 15.29 — P12 10 14.17 7.3 P13 20 12.29 19.6 P14 30 9.99 34.7 P15 45 8.43 44.9 P21 0 22.70 — P22 10 19.80 12.8 P23 20 17.35 23.6 P24 30 15.67 30.9 P25 45 13.57 40.2 P31 0 37.94 — P32 10 35.84 5.5 P33 20 33.87 10.7 P34 30 32.59 14.1 P35 45 29.33 22.8 表 5 承台和桩的计算参数表
Table 5 Computational parameters of caps and piles
名称 弹性模量/GPa 密度/(kg·m-3) 泊松比 有机玻璃桩 3.2 180 0.3 铝桩及承台 71 2700 0.2 钢桩及承台 124 7900 0.2 表 6 地基土计算参数表
Table 6 Computational parameters of foundation soil
重度/(kN·m-3) 内摩擦角/(°) 弹性模量/MPa 泊松比 17.6 40.9 63 0.3 -
[1] 郑刚, 何晓佩, 周海祚, 等. 基坑斜-直交替支护桩工作机理分析[J]. 岩土工程学报, 2019, 41(增刊1): 97-100. doi: 10.11779/CJGE2019S1025 ZHENG Gang, HE Xiaopei, ZHOU Haizuo, et al. Working mechanism of inclined-vertical retaining piles in excavations[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(S1): 97-100. (in Chinese) doi: 10.11779/CJGE2019S1025
[2] GUO L, ZHANG Z Q, WANG X B, et al. Stability analysis of transmission tower foundations in permafrost equipped with thermosiphons and vegetation cover on the Qinghai-Tibet Plateau[J]. International Journal of Heat and Mass Transfer, 2018, 121: 367-376. doi: 10.1016/j.ijheatmasstransfer.2018.01.009
[3] ZHENG G, YANG X Y, ZHOU H Z, et al. A simplified prediction method for evaluating tunnel displacement induced by laterally adjacent excavations[J]. Computers and Geotechnics, 2018, 95: 119-128. doi: 10.1016/j.compgeo.2017.10.006
[4] 吕凡任. 倾斜荷载作用下斜桩基础工作性状研究[D]. 杭州: 浙江大学, 2004. LÜ Fanren. Study on Behavior of Battered Piles under Inclined Load[D]. Hangzhou: Zhejiang University, 2004. (in Chinese)
[5] 王新泉, 陈永辉, 安永福, 等. 塑料套管现浇混凝土桩倾斜对承载性能影响的模型试验研究[J]. 岩石力学与工程学报, 2011, 30(4): 834-842. WANG Xinquan, CHEN Yonghui, AN Yongfu, et al. Model test study of effect of inclination on bearing behaviors of plastic tube cast-in-place concrete pile[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(4): 834-842. (in Chinese)
[6] SABBAGH T T, AL-SALIH O, AL-ABBOODI I. Experimental investigation of batter pile groups behaviour subjected to lateral soil movement in sand[J]. International Journal of Geotechnical Engineering, 2020, 14(7): 705-716. doi: 10.1080/19386362.2019.1585596
[7] MEYERHOF G G, YALCIN A S. Behaviour of flexible batter piles under inclined loads in layered soil[J]. Canadian Geotechnical Journal, 1993, 30(2): 247-256. doi: 10.1139/t93-021
[8] ZHANG L M, MCVAY M C, LAI P W. Centrifuge modelling of laterally loaded single battered piles in sands[J]. Canadian Geotechnical Journal, 1999, 36(6): 1074-1084. doi: 10.1139/t99-072
[9] HAZZAR L, HUSSIEN M N, KARRAY M. Numerical investigation of the lateral response of battered pile foundations[J]. International Journal of Geotechnical Engineering, 2017, 11(4): 376-392. doi: 10.1080/19386362.2016.1224030
[10] BAJAJ P, YADU L, CHOUKSEY S K. Study on vertical and batter piles subjected to lateral loads in different non-cohesive sub-soil conditions[J]. International Journal of Geotechnical Engineering, 2020, 14(6): 603-613. doi: 10.1080/19386362.2018.1564181
[11] KONG D S, DENG M X, LI Y Z. Experimental study on mechanical deformation characteristics of inclined and straight alternating pile groups[J]. Advances in Civil Engineering, 2020, 2020: 1-11.
[12] RAJASHREE S S, SITHARAM T G. Nonlinear finite-element modeling of batter piles under lateral load[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2001, 127(7): 604-612. doi: 10.1061/(ASCE)1090-0241(2001)127:7(604)
[13] 黄剑慧, 陈红兵, 刘杰伟, 等. 群桩中的斜桩受力分析与计算[J]. 广东土木与建筑, 2007, 14(8): 21-23. HUANG Jianhui, CHEN Hongbing, LIU Jiewei, et al. Force analysis and calculation of inclined piles in pile groups[J]. Guangdong Architecture Civil Engineering, 2007, 14(8): 21-23. (in Chinese)
[14] LING D, REN T, WANG Y. A p-y curve method for horizontal bearing characteristics of single batter pile in sands[J]. Rock Soil Mech, 2013, 34(1): 155-162.
[15] REESE L C, COX W R, KOOP F D. Analysis of laterally loaded piles in sand[C]//Proceedings of the Sixth Annual Offshore Technology Conference. Texas, 1974.
[16] 曹卫平, 吴奇兴, 夏冰, 等. 考虑上拔荷载影响的水平受荷斜桩p-y曲线[J]. 岩土工程学报, 2020, 42(7): 1189-1198. doi: 10.11779/CJGE202007002 CAO Weiping, WU Qixing, XIA Bing, et al. p-y curve of laterally loaded batter piles considering effect of uplift loads[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(7): 1189-1198. (in Chinese) doi: 10.11779/CJGE202007002
[17] 樊文甫. 水平受荷斜桩承载变形性状数值分析[D]. 西安: 西安建筑科技大学, 2016. FAN Wenfu. Numerical Analyses on the Behavior of Horizontally Loaded Batter Piles[D]. Xi'an: Xi'an University of Architecture and Technology, 2016. (in Chinese)
[18] 曹卫平, 葛欣. 水平受荷斜桩承载变形性状及荷载传递机理分析[J]. 西安建筑科技大学学报(自然科学版), 2017, 49(5): 624-629. CAO Weiping, GE Xin. Analyses of bearing capacity and load transfer mechanism of horizontally loaded batter piles[J]. Journal of Xi'an University of Architecture & Technology (Natural Science Edition), 2017, 49(5): 624-629. (in Chinese)
[19] 王恩钰, 周海祚, 郑刚, 等. 基坑倾斜桩支护的变形数值分析[J]. 岩土工程学报, 2019, 41(增刊1): 73-76. doi: 10.11779/CJGE2019S1019 WANG Enyu, ZHOU Haizuo, ZHENG Gang, et al. Numerical analyses of deformation of inclined pile-retained excavations[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(S1): 73-76. (in Chinese) doi: 10.11779/CJGE2019S1019
[20] 赵爽, 吴君涛, 张日红, 等. 砂土中斜桩单桩水平承载与变形特性数值分析[J]. 中南大学学报(自然科学版), 2022, 53(2): 579-588. ZHAO Shuang, WU Juntao, ZHANG Rihong, et al. Numerical analysis of horizontal bearing capacity and deflection behaviors of single batter pile in sand[J]. Journal of Central South University (Science and Technology), 2022, 53(2): 579-588. (in Chinese)
[21] 土工试验方法标准: GB/T 50123—2019[S]. 北京: 中国计划出版社, 2019. Standard for Soil Test Method: GB/T 50123—2019[S]. Beijing: China Planning Press, 2019. (in Chinese)
[22] 建筑基桩检测技术规范: JGJ 106—2014[S]. 北京: 中国建筑工业出版社, 2014. Technical Code for Testing of Building Foundation Piles: JGJ 106—2014[S]. Beijing: China Architecture & Building Press, 2014. (in Chinese)
[23] 袁志林, 段梦兰, 陈祥余, 等. 水平荷载下导管架平台桩基础的非线性有限元分析[J]. 岩土力学, 2012, 33(8): 2551-2560. YUAN Zhilin, DUAN Menglan, CHEN Xiangyu, et al. Nonlinear finite element analysis of jacket platform pile foundations under lateral loads[J]. Rock and Soil Mechanics, 2012, 33(8): 2551-2560. (in Chinese)
[24] 刘晋超, 熊根, 朱斌, 等. 砂土海床中大直径单桩水平承载与变形特性[J]. 岩土力学, 2015, 36(2): 591-599. LIU Jinchao, XIONG Gen, ZHU Bin, et al. Bearing capacity and deflection behaviors of large diameter monopile foundations in sand seabed[J]. Rock and Soil Mechanics, 2015, 36(2): 591-599. (in Chinese)
[25] 张磊, 周腾, 芮瑞, 等. 软黏土场地桩基桥墩系统地震损伤特性的三维有限元动力分析[J]. 振动与冲击, 2021, 40(23): 109-119. ZHANG Lei, ZHOU Teng, RUI Rui, et al. Three-dimensional finite element analyses of clay-pile-pier systems subjected to seismic motions[J]. Journal of Vibration and Shock, 2021, 40(23): 109-119. (in Chinese)
[26] 混凝土结构设计规范: GB 50010—2010[S]. 北京: 中国建筑工业出版社, 2011. Code for Design of Concrete Structures: GB 50010—2010[S]. Beijing: China Architecture & Building Press, 2011. (in Chinese)
[27] 郑刚, 王玉萍, 程雪松, 等. 基坑倾斜桩支护性能及机理大型模型试验研究[J]. 岩土工程学报, 2021, 43(9): 1581-1591. doi: 10.11779/CJGE202109002 ZHENG Gang, WANG Yuping, CHENG Xuesong, et al. Large-scale model tests on performance and mechanism of inclined retaining structures of excavations[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(9): 1581-1591. (in Chinese) doi: 10.11779/CJGE202109002
[28] 郑刚, 王丽. 竖向及水平荷载加载水平、顺序对单桩承载力的影响[J]. 岩土工程学报, 2008, 30(12): 1796-1804. http://www.cgejournal.com/cn/article/id/13062 ZHENG Gang, WANG Li. Effect of loading level and sequence of vertical and lateral load on bearing capacity of single pile[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(12): 1796-1804. (in Chinese) http://www.cgejournal.com/cn/article/id/13062
-
期刊类型引用(1)
1. 朱晟,陆雪妮,汤洪洁,孟涛,路德任. 大粒径砂砾料的缩尺试验方法及应用. 水利规划与设计. 2022(12): 5-10 . 百度学术
其他类型引用(1)
-
其他相关附件