Loading [MathJax]/jax/output/SVG/jax.js
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

地应力对岩体爆破块度分布特征的影响研究

曹泽铭, 杨建华, 叶志伟, 冷振东, 姚池, 张小波

曹泽铭, 杨建华, 叶志伟, 冷振东, 姚池, 张小波. 地应力对岩体爆破块度分布特征的影响研究[J]. 岩土工程学报, 2024, 46(10): 2202-2211. DOI: 10.11779/CJGE20230714
引用本文: 曹泽铭, 杨建华, 叶志伟, 冷振东, 姚池, 张小波. 地应力对岩体爆破块度分布特征的影响研究[J]. 岩土工程学报, 2024, 46(10): 2202-2211. DOI: 10.11779/CJGE20230714
CAO Zeming, YANG Jianhua, YE Zhiwei, LENG Zhendong, YAO Chi, ZHANG Xiaobo. Influences of in-situ stress on distribution characteristics of rock blasting fragmentation[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(10): 2202-2211. DOI: 10.11779/CJGE20230714
Citation: CAO Zeming, YANG Jianhua, YE Zhiwei, LENG Zhendong, YAO Chi, ZHANG Xiaobo. Influences of in-situ stress on distribution characteristics of rock blasting fragmentation[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(10): 2202-2211. DOI: 10.11779/CJGE20230714

地应力对岩体爆破块度分布特征的影响研究  English Version

基金项目: 

国家自然科学基金项目 52179102

国家自然科学基金项目 51969015

江西省自然科学基金项目 20224BAB214079

江西省自然科学基金项目 20232BAB204090

详细信息
    作者简介:

    曹泽铭(1999—),男,硕士研究生,主要从事岩石动力学及工程爆破数值模拟方面的研究工作。E-mail: 406000210058@email.ncu.edu.cn

    通讯作者:

    杨建华, E-mail: yangjianhua86@ncu.edu.cn

  • 中图分类号: TU432

Influences of in-situ stress on distribution characteristics of rock blasting fragmentation

  • 摘要: 深地工程高地应力岩体爆破开挖时常出现岩体破碎难、大块率高的问题。为深入理解这一问题,采用动力有限元方法模拟不同地应力条件下的岩体爆破开裂过程,引入爆生裂纹分形维数和岩体爆破块度图像识别方法,研究地应力对岩体爆破爆生裂纹分布和块度分布特征的影响,并通过不同埋深和地应力水平下的岩体爆破工程实例进行验证。研究结果表明:地应力对炮孔远区爆生裂纹扩展的抑制作用以及非静水地应力对爆生裂纹扩展的导向作用,导致了岩体爆破大块的产生;随着地应力水平的提高,岩体爆破平均块度尺寸、最大块度尺寸、块度不均匀系数和大块率均显著增加;炮孔近区的爆生裂纹扩展几乎不受地应力的影响,从而使得地应力对岩体爆破小块度的尺寸影响不大;垂直于炮孔轴向的两个主地应力差别较大时,不利于爆破破岩,侧压力系数约为0.75时,岩体破碎效果较好。
    Abstract: During blasting excavation of deep rock masses under high in-situ stress, the rock masses are often difficult to be fully broken, frequently resulting in large fragments. To understand this problem more deeply, the blast-induced rock cracks under various in-situ stress conditions are simulated by using the dynamic finite element method. The fractal dimension theory and the image recognition method are introduced to investigates the influence of in-situ stress on the distribution characteristics of rock blasting fragmentation. The project cases blasted at different depths and in-situ stress levels are also employed to study this problem. The results show that the generation of large fragments during blasting in deep rock masses owes to the inhibiting effect of in-situ stress on the growth of blast-induced cracks in the far field of blastholes, and the orienting effects of the non-hydrostatic in-situ stress on the propagation of blast-induced cracks. As the in-situ stress level increases, average size, the maximum size, nonuniform coefficient and large block rate of rock fragmentation increase significantly. The size of the smaller fragments is little affected by the in-situ stress as these fragments are generated in the vicinity of blastholes and the propagation of blasting cracks in this zone is almost unaffected by the in-situ stress. When the maximum and minimum principal in-situ stresses perpendicular to the blasthole axis differ greatly, it is adverse to rock fragmentation by blasting. In comparison, a better rock fragmentation is achieved when the lateral coefficient of the principal in-situ stresses is about 0.75.
  • 扩体预制桩由内芯预制桩和外围扩体材料固结体组成,一般可通过就地搅拌水泥土或机械扩孔并向孔内灌注扩体材料浆体,再植入内芯同轴预制桩的工艺成桩[1]。近年来,随着厂拌高质量扩体浆材的利用和成桩工艺的创新[2],扩体预制桩的桩身质量和承载性能大幅提高,开始被广泛应用于高层建筑、桩板路桥等承载力和沉降控制要求较高的建构筑物桩基础。

    与单一材料桩相比,扩体预制桩竖向承载时涉及到内芯预制桩、外围扩体和周围土3种不同材质的相互作用,荷载传递机制十分复杂[3-4]。然而,目前有关扩体预制桩承载机制的理论分析研究极为有限。考虑芯桩-扩体-周围土协同工作对承载性能的影响,Wang等[5]在假定芯桩与外围水泥土界面固定黏结的基础上,基于荷载传递法建立了劲性搅拌桩荷载传递简化计算模型;俞建霖等[6]针对刚性基础下砼芯水泥土桩复合地基,给出了砼芯、外围水泥土、周围土三者的受力平衡微分方程;结合引扩孔灌浆扩体预制桩工艺特点,张浩等[7]在考虑芯桩−扩体材料−周围土径向相互作用的基础上,提出了扩体预制桩整体复合模量表达式,并给出了全截面竖向受压模式下扩体桩的承载变形解析计算方法;考虑芯桩-水泥土-周围土界面的非线性力学特征,Zhu等[8]通过三折线模型建立了刚性基础下劲性搅拌桩的荷载传递理论分析模型,极大完善了扩体预制桩荷载传递机制的分析理论。然而,这些研究多集中于桩顶全截面受压工况,即桩顶处芯桩与扩体近似等沉变形受荷,较适用于建筑筏板、桥台等刚性基础下的桩基工况。

    近年来,随着扩体预制桩在桩板式无土路基结构、高桩承台结构中的广泛应用,芯桩单独承载的扩体预制桩荷载传递机制研究不容忽视。既有试验研究表明[9],芯桩单独承压下扩体预制桩中芯桩-扩体-土的荷载传递机制与全截面受压下显著不同。顾士坦等[10]基于Mindlin弹性解和剪滞原理,仅考虑芯桩-外围水泥土的相互作用,给出了芯桩单独受力时劲性搅拌桩芯桩荷载传递规律的理论分析方法;任连伟等[11]、Zhou等[12]基于荷载传递法,分别给出了高喷插芯组合桩和静钻根植桩芯桩受力时荷载传递简化计算方法,但芯桩-水泥土界面完全线弹性剪切模型的假定,无法全面反映此类扩体预制桩的工作状态。

    据此,本文针对芯桩单独承载工况下全长扩体预制桩中芯桩-扩体-周围土的荷载传递特点,在充分考虑芯桩-扩体、扩体-土界面非线性相互作用特性以及两界面剪切性能发挥交互影响的基础上,构建了全长扩体预制桩芯桩承载荷载传递计算分析模型,并参考工程实例,对全长扩体预制桩芯桩承载荷载传递机制及其影响因素进行了分析。

    图 1为典型全长扩体预制桩芯桩单独承载工况下桩土作用示意图。荷载P直接作用于内芯预制桩顶部,芯桩首先向下发生相对外围扩体材料的变形,并在芯桩-扩体界面产生剪切作用;外围扩体材料发生剪切变形,进一步将荷载传递至周围土。同时随着荷载向下传递,桩端处发生向下刺入变形趋势,并受到地基土的反力作用。因此芯桩承载工况下上部荷载是经由芯桩-扩体-周围土相互作用传递至周围土的。芯桩承载扩体预制桩的承载变形机制涉及到芯桩-扩体、扩体-周围土内外两个相互作用界面的协同工作。扩体预制桩整体受到周围土侧阻力τs、芯桩端阻力Rp和扩体端阻力Rc,内部则存在芯桩与扩体材料的剪切作用力τc

    图  1  承载变形机制示意图
    Figure  1.  Diagrammatic sketch of bearing mechanism

    图 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)为芯桩与土的相对位移。

    既有研究表明[8, 13],理想弹塑性模型能较好的反映水泥土、水泥砂浆等水泥基胶凝材料与土之间的荷载传递关系,因此将扩体-周围土界面模型简化为理想弹塑性模型,如图 2(a)所示,即

    τs={ksscs (scs<sfs)τfs(scssfs) (2)
    图  2  界面模型示意图
    Figure  2.  Diagrammatic sketch of interface model

    式中:τsscs为扩体-土界面侧阻力和相对位移;ks为扩体-土剪切刚度;τfssfs界面极限侧阻力和极限相对位移。

    芯桩与水泥土等水泥基胶凝材料的剪切试验[14]表明,在较小相对位移下界面剪切应力即表现出较高水平,随着相对位移的增大,界面应力达到极限值后会出现迅速降低,界面呈脆性破坏模式,如图 2(b)所示。因此,芯桩-扩体界面不宜采用完全线弹性剪切模型,可简化为图 2(b)理想弹性-破坏模型,即

    τc={kcspc(spc<sfc)τfc (spcsfc) (3)

    式中:τcspc为芯桩-扩体界面侧阻力和相对位移;kc为芯桩-扩体界面剪切刚度,由试验得τc-spc曲线的直线段斜率确定;τfcsfc为界面极限侧阻力和相对位移;τfc为界面残余侧阻力,与扩体材料水泥掺量、龄期等有关,约为界面极限侧阻力的(0.057~0.153)倍[14]

    据此,分析时作如下基本假定:①内芯预制桩和外围扩体材料均为线弹性体;②针对目前以水泥土、水泥砂浆为扩体材料的扩体预制桩,扩体剪切刚度远大于周围土,且厚度相对有限[1-2],受荷时扩体剪切变形较小,忽略不计;③假定芯桩-扩体荷载传递为理想弹性-破坏模型,扩体-土荷载传递为理想弹塑性模型。

    图 3所示,取深度z处“扩体预制桩-土”单元进行受力分析,根据单元体静力平衡条件可得

    图  3  单元体受力变形示意图
    Figure  3.  Diagrammatic sketch of forces and deformations on elements
    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处芯桩-扩体、扩体-土界面的侧阻力;upuc为芯桩、扩体外侧周长。

    若受力变形后芯桩单元体和扩体单元体的压缩量分别为ΔpΔc,根据图 3(b)变形模式可得

    dspc(z)=ΔcΔpdscs(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)

    式中:EpEc为芯桩、扩体的弹性模量;ApAc为芯桩、扩体的横截面积。

    芯桩承载扩体预制桩与土相互作用模式见图 4

    图  4  扩体预制桩-土相互作用模式
    Figure  4.  Reamed precast pile-soil interaction model

    随着桩顶荷载增加,扩体-土界面和芯桩-扩体界面将由完全弹性状态、从浅而深逐渐进入塑性或破坏状态。设桩长L,扩体-土界面塑性区深度L1,芯桩-扩体界面破坏区深度L2。一般情况下,由于芯桩-扩体界面剪切强度比扩体-土界面要大,工程实践中为充分发挥芯桩的承载性能,往往会通过调配扩体材料强度和厚度来保证芯桩-扩体共同作用,避免芯桩-扩体界面提前破坏[1],此时L1往往大于L2,如图 4(a)

    当外围扩体强度较小或厚度较薄等特殊情况下,芯桩单独承载桩顶附近芯桩和扩体应力差异显著,芯桩-扩体界面可能先出现破坏,而后扩体-土界面逐渐进入塑性状态,桩-土相互作用为图 4(b)所示模式。

    (1)L2L1工况

    根据芯桩-扩体、扩体-土两界面相互作用状态的不同,将扩体预制桩划分为3个区域,见图 4(a)

    a)当L1 < zL时,芯桩-扩体界面、扩体-土界面均处于弹性阶段,由式(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)

    式中:C1C2C3C4为方程待定系数;

    a1=[(δ+α)+(δ+α)24(δαγβ)]/2
    a2=[(δ+α)+(δ+α)24(δαγβ)]/2
    a3=[(δ+α)(δ+α)24(δαγβ)]/2
    a4=[(δ+α)(δ+α)24(δαγβ)]/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)/Gsrc为扩体外壁半径;Gs为桩周土剪切模量;rm为影响半径[15]

    b)当L2 < zL1时,扩体-周围土界面处于塑性状态,芯桩-扩体界面处于弹性阶段,

    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)

    式中:D1D2D3D4为方程待定系数。

    进而,芯桩轴力Pp2(z)和扩体轴力Pc2(z)

    Pp2(z)=EpAp[γ+αα(D1eαzD2eα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(rmrc)/τsfrcln(rmrc)GsGs

    c)当0≤zL2时,扩体-周围土界面处于塑性状态,芯桩-扩体界面处于破坏阶段,

    d2spc(z)dz2=g, d2scs(z)dz2=f} (15)

    式中:g=(upEcAc+upEpAp)τfcucτfsEcAcf=ucτfsEcAcupτfcEcAc

    求解方程,可得z处芯桩-扩体相对位移spc3(z)、扩体-周围土相对位移scs3(z)

    spc3(z)=12gz2+F1z+F2scs3(z)=12fz2+F3z+F4} (16)

    式中:F1F2F3F4为方程待定系数。

    进而,芯桩轴力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=PPc(z)|z=0=0} (19)

    桩端边界条件为

    Pp(z)|z=L=kLApsp(z)|z=LPc(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)

    综上,针对L2L1工况,根据边界和连续条件(19)~(22),可得关于C1~C4D1~D4F1~F4L1、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 < LL2=0),上述计算退化为关于C1~C4D1~D4L1共9个未知量的求解问题,此时将式(8),(9),(12),(13)代入边界和连续条件(19)~(21)即可得出9个求解方程,进而可得该工况下桩身受力响应;③随着桩顶荷载的增加,芯桩-扩体界面也开始出现脆性破坏(0 < L1 < L,0 < L2 < L1),即进入式(7)~(18)计算工况;④当扩体-土界面全部进入塑性区而芯桩-扩体界面局部出现破坏时(L1=L,0 < L2 < L1),上述计算退化为关于D1~D4F1~F4L2共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 < zL时,芯桩-扩体界面、扩体-土界面均处于弹性阶段;②L1 < zL2时,芯桩-扩体界面处于破坏状态,而外部扩体-土界面仍处于弹性状态;③0 < zL1时,扩体-土界面处于塑性状态,芯桩-扩体界面处于破坏阶段。

    其中,在L2 < zL和0 < zL1区域,桩-土相互作用状态与前述2.3节中扩体桩上、下两段区域桩-土相互作用状态相同,芯桩-扩体、扩体-土两界面的微分控制方程与式(7),(15)相同,求解方法也相同,在此不再赘述。

    L1 < zL2时,由式(2),(3),(6)可得此时芯桩-扩体和扩体-土界面微分方程:

    d2spc(z)/dz2=ξd2scs(z)/dz2=δscs(z)ζ} (23)

    式中:ξ=upτfc/EpApζ=upτfc/EcAc

    求解方程,可得该区域z处芯桩位移spc2(z)、扩体-周围土相对位移scs2(z)

    spc2(z)=12ξz2+J1z+J2zscs2(z)=J3eζz+J4eζz+ζ/ξ} (24)

    式中:J1J2J3J4为方程待定系数。

    进而,芯桩轴力Pp2(z)、扩体轴力Pc2(z)和芯桩竖向变形s2(z)

    Pp2(z)=EpAp[ξz+J1]Pc2(z)=EcAc[J3eζz+J4eζz]} (25)
    s2(z)=spc2(z)+scs2(z)+ss2(z) (26)

    式中:s2(z)=τs(z)rcln(rm/rc)/Gs

    据此,可对L2 > L1工况下扩体预制桩荷载传递进行计算分析。其中,桩顶和桩端边界条件同式(19),(20);桩身受力变形连续条件为

    spc2(z)|z=L2=spc1(z)|z=L2+scs2(z)|z=L2=scs1(z)|z=L2+Pp2(z)|z=L2=Pp1(z)|z=L2+Pc2(z)|z=L2=Pc1(z)|z=L1+spc2(z)|z=L2=sfc} (27)
    spc3(z)|z=L1=spc2(z)|z=L1+scs3(z)|z=L1=scs2(z)|z=L1+Pp3(z)|z=L1=Pp2(z)|z=L1+Pc3(z)|z=L1=Pc3(z)|z=L1+scs3(z)|z=L1=sfs} (28)

    式中:spc1(z)scs2(z)Pp1Pc1(z)分别为L2 < zL芯桩和扩体的变形和轴力;spc2(z)scs2(z)Pp2(z)Pc2(z)分别为L1 < zL2芯桩和扩体的变形和轴力;spc3(z)scs3(z)Pp3(z)Pc3(z)分别为0 < zL1芯桩和扩体的变形和轴力。

    考虑不同计算条件,采用与2.3节相同的求解方法,可得L2 > L1工况下扩体桩桩身轴力与位移。

    水泥土扩体预制桩现场试验参见文献[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。

    采用本文方法对现场试桩进行计算,如图 56。各级荷载作用下,桩顶沉降计算值大于实测值,芯桩轴力计算值整体上与实测值相比偏小,这是由于实际芯桩后植入过程中会对水泥土和周围土产生挤密作用,使界面黏结强度增大[2],桩体整体承载性能提升,但计算分析中界面参数未考虑该部分的影响。

    图  5  实测与计算沉降-荷载曲线对比
    Figure  5.  Comparison between measured and calculated curves of displacement-load
    图  6  芯桩轴力与桩身压缩量的实测与计算结果对比
    Figure  6.  Comparison of measured and calculated axial force of core pile and settlement compression of pile

    从芯桩轴力分布规律来看,芯桩轴力沿深度逐渐减小,且随着上部荷载增加,该变化趋势越来越明显。各级荷载下桩端轴力约占桩顶总荷载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
    下载: 导出CSV 
    | 显示表格

    德上高速合肥-枞阳段地处安徽中南部,部分路段采用了“桩柱一体式”桩板路基结构。考虑沿线基岩较浅且上覆沉积物工程性能较好,传统打入式预制桩沉桩困难,因此采用了引扩孔灌浆植桩工艺。芯桩采用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/
    MPa
    fsk/
    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
    下载: 导出CSV 
    | 显示表格

    为了获取芯桩-砂浆、砂浆-土界面荷载传递参数,采用自研界面剪切装置(图 7)开展界面剪切试验,试验方法同直接剪切试验,每组试验3个试样,分别考虑100,200,300 kPa 3种不同的法向应力,参考《土工试验方法标准》(GB/T50123—2019),每个法向力采用一次轻轻施加。受制于剪切试样尺寸(ϕ61.8 mm×20 mm)较小,芯桩混凝土采用同强度C80灌浆料替代,所用土样和水泥砂浆皆取自工程现场。

    图  7  界面剪切试验示意图
    Figure  7.  Diagrammatic sketch of interface shear tests

    制作水泥砂浆-土试样时,为保证试样土重度与工程实际的一致性,首先参考地勘报告调整土样含水率,并以勘察报告土体重度为基准,根据下剪切盒体积,计算每个试样所需土的重量,将其填入下剪切盒至剪切界面高度,然后在其上剪切盒内浇入水泥砂浆。制作芯桩-水泥砂浆试样时,首先将替代灌浆料浇入模具中,待标准养护28 d后再在其上浇入水泥砂浆。制样完成后,所有试样标准养护14 d后开展界面剪切试验,部分成果如图 8所示。其中,散点曲线为试验测试结果,黑色实折线为试验数据拟合模型曲线。

    图  8  界面剪切试验结果
    Figure  8.  Results of interface shear tests

    依据界面剪切试验,取芯桩-砂浆界面极限侧阻力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),桩顶沉降变化量降低。

    图  9  实测值与计算值对比
    Figure  9.  Comparison between measured and calculated values

    采用本文方法计算的荷载-位移曲线与实测结果具有较好的一致性,同级荷载作用下计算桩顶沉降相较实测值整体偏大。当扩径为700 mm时,各级荷载水平下桩顶沉降平均误差率6.7%;随着扩径直径增大,计算值误差率随之增大。当扩径直径由700 mm增加至800,1000,1200 mm时,桩顶沉降平均误差率分别增加至24.8%,34.6%,31.1%。分析其原因,上述计算分析未考虑芯桩后植入施工的扩挤影响,尤其是后植入芯桩对孔内砂浆挤密及砂浆渗滤增强桩-土界面剪切性能的作用效应,同时简化理论模型中忽略外围扩体剪切变形的基本假定也对计算结果有影响。考虑相关技术规程和工程实践中[1]为了充分发挥芯桩高承载力的特性和合理控制孔内灌浆后植入桩的挤土效应,推荐扩孔直径宜为芯桩直径的1.5倍,即扩体厚度多控制在100~200 mm,因此本文计算方法适用于常规扩体预制桩的工程计算分析。

    基于上述实例,进一步计算分析复合桩身荷载传递机制。以扩径直径700 mm(内芯预制桩500 mm)水泥砂浆扩体预制桩为例,图 10所示为计算的桩身轴力分布曲线。不同荷载水平下芯桩轴力沿桩身呈逐渐减小的变化趋势,随着桩顶荷载水平的增大,在上部0~3 m内芯桩轴力减小幅度相对较大,3~8 m内芯桩轴力减小幅度变小且相对稳定,而桩端附近(8~10 m)芯桩轴力减小幅度再次变大。

    图  10  扩体预制桩桩身轴力
    Figure  10.  Axial forces of reamed precast pile

    与之对应,外围水泥砂浆截面轴力沿桩身呈逐渐增大变化趋势,在0~3 m内外围砂浆轴力快速增加,说明芯桩单独承载工况下,桩顶附近芯桩和外围砂浆的相互作用较强;随着深度增加,外围砂浆轴力增加幅度降低且相对稳定,与芯桩趋于共同作用;桩端(8~10 m)附近,由于芯桩与外围砂浆刚度的显著差异,以及桩端角砾岩较高桩端反力作用,芯桩相对外围砂浆“下刺入”过程中,两者相互作用增强,外围砂浆轴力呈明显增大的变化趋势,如图 10(b)所示。

    图 11所示为水泥砂浆扩体预制桩内、外两界面阻力分布曲线。内侧芯桩-水泥砂浆界面阻力沿桩身呈两端大、中间小的分布规律。这主要是由于在桩顶、桩端附近芯桩和外围砂浆的刚度差异导致两者界面相对位移趋势明显,引起较强烈的界面相互作用。其中桩顶荷载4000 kN时,芯桩-水泥砂浆界面在桩顶0.2 m范围内出现剪切破坏,界面阻力降为残余侧阻力。

    图  11  扩体预制桩界面阻力
    Figure  11.  Interfacial resistances of reamed precast pile

    对于外侧水泥砂浆-土界面,界面阻力沿桩身整体呈逐渐减小分布。当桩顶荷载较大时,水泥砂浆-土界面由上而下逐渐进入塑性状态,界面阻力达到极限侧阻力。随着桩顶荷载增加,界面塑性区长度增大,桩顶荷载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%。外围水泥砂浆对芯桩具有“扩径”作用。

    图  12  不同扩径情况下桩身轴力
    Figure  12.  Axial forces under different reaming diameters

    当扩体直径相对较大(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%,说明外围扩体材料刚度的改善,可提升内芯预制桩向周围传递荷载的效率。

    图  13  不同扩体材料情况下桩身轴力
    Figure  13.  Axial forces under different peripheral materials

    由扩体截面轴力分布曲线可以看出,当扩体采用刚度相对较低(0.8 GPa)的水泥土时,芯桩与扩体的显著相互作用仅发生在桩顶(1 m)和桩端(2 m)附近有限区域内;随着扩体刚度增加,芯桩与扩体的显著相互作用区域逐渐增大,两者共同作用效应增强。在工程实践中可通过改善扩体材料性能,增强内芯预制桩与外围扩体的共同作用。

    (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   深埋隧洞开挖爆破后的大块岩石

    Figure  1.   Large fragments created by blasting in a deep tunnel

    图  2   高地应力岩体爆破计算模型

    Figure  2.   Computational model for rock blasting under high in-situ stress

    图  3   花岗岩单孔爆破裂纹分布试验与数值模拟结果对比

    Figure  3.   Comparison of blast-induced cracks in granite between experimental and numerical results

    图  4   爆生裂纹计盒维数计算流程图

    Figure  4.   Flow chart for calculating box-counting fractal dimension of blast-induced cracks

    图  5   基于ImageJ图像处理的岩体爆破块度统计流程图

    Figure  5.   Flow chart for analyzing rock blasting fragmentation based on ImageJ

    图  6   不同地应力水平下的岩体爆生裂纹分布

    Figure  6.   Distributions of blast-induced cracks in rock under different in-situ stress levels

    图  7   不同地应力水平下爆生裂纹分形维数和裂纹面积比

    Figure  7.   Fractal dimensions and scaled areas of blast-induced cracks under different in-situ stress levels

    图  8   不同地应力水平下的岩体爆破块度分布曲线

    Figure  8.   Distribution curves of rock blasting fragmentation under different in-situ stress levels

    图  9   不同地应力水平下岩体爆破块度的特征尺寸

    Figure  9.   Characteristic sizes of rock blasting fragmentation under different in-situ stress levels

    图  10   不同地应力水平下爆破块度不均匀系数和大块率

    Figure  10.   Nonuniform coefficients and large block rates of rock blasting fragmentation under different in-situ stress levels

    图  11   不同侧压力系数下的岩体爆生裂纹分布(σv=30 MPa)

    Figure  11.   Distributions of blast-induced cracks in rock under different lateral coefficients (σv=30 MPa)

    图  12   不同侧压力系数下爆生裂纹分形维数和裂纹面积比(σv=30 MPa)

    Figure  12.   Fractal dimensions and scaled areas of blast-induced cracks under different lateral coefficients (σv=30 MPa)

    图  13   不同侧压力系数下的岩体爆破块度分布曲线(σv=30 MPa)

    Figure  13.   Distribution curves of rock blasting fragmentation under different lateral coefficients (σv=30 MPa)

    图  14   不同侧压力系数下的岩体爆破块度特征尺寸

    Figure  14.   Characteristic sizes of rock blasting fragmentation under different lateral coefficients

    图  15   不同侧压力系数下爆破块度不均匀系数和大块率

    Figure  15.   Nonuniform coefficients and large block rates of rock blasting fragmentation under different lateral coefficients

    图  16   斜井不同埋深岩体爆破后的爆堆

    Figure  16.   Piles of rock blasting at different depths in inclined shaft

    图  17   不同埋深岩体爆破块度分布曲线

    Figure  17.   Distribution curves of rock blasting fragmentation at different depths

    图  18   不同埋深岩体爆破块度特征尺寸对比

    Figure  18.   Comparison of characteristic sizes of rock blasting fragmentation at different depths

    表  1   炸药材料参数

    Table  1   Parameters of explosive

    ρ/(kg⋅m-3) D/(m⋅s-1) C1/GPa C2/GPa R1 R2 ω E0/GPa
    880 5170 348.6 11.29 7.0 2.0 0.24 5.025
    下载: 导出CSV

    表  2   花岗岩RHT模型参数

    Table  2   RHT model parameters of granite

    参数 数值 参数 数值
    密度ρ0/(kg⋅m-3) 2660 初始孔隙度α0 0
    压碎压力pel/MPa 125 孔隙度指数N 3.0
    压实压力pcomp/GPa 6.0 参考压缩应变率EOC/s-1 3×10-5
    Hugoniot多项式参数A1/GPa 25.7 参考拉伸应变率EOT/s-1 3×10-6
    Hugoniot多项式参数A2/GPa 37.84 失效压缩应变率EC/s-1 3×1025
    Hugoniot多项式参数A3/GPa 21.29 失效拉伸应变率ET/s-1 3×1025
    侵蚀塑性应变EPSF 2.0 最小失效应变EPM 0.015
    状态方程参数B0 1.22 压缩应变率指数βc 0.026
    状态方程参数B1 1.22 拉伸应变率指数βt 0.007
    状态方程参数T1/GPa 25.7 压缩屈服面参数gc 0.53
    状态方程参数T2 0.0 拉伸屈服面参数gt 0.7
    剪切模量G/GPa 21.9 失效面参数A 2.44
    单轴抗压强度fc/MPa 167.8 剪切模量缩减系数ξ 0.5
    拉压强度比ft 0.04 损伤参数D1 0.04
    剪压强度比fs 0.21 损伤参数D2 1.0
    失效面指数n 0.76 残余应力强度参数Af 0.25
    拉压子午比参数Q0 0.68 残余应力强度指数nf 0.62
    罗德角相关参数B 0.05 Grunesien参数Γ 0.0
    下载: 导出CSV
  • [1] 谢和平, 李存宝, 高明忠, 等. 深部原位岩石力学构想与初步探索[J]. 岩石力学与工程学报, 2021, 40(2): 217-232.

    XIE Heping, LI Cunbao, GAO Mingzhong, et al. Conceptualization and preliminary research on deep in situ rock mechanics[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(2): 217-232. (in Chinese)

    [2]

    HE C L, YANG J, YU Q. Laboratory study on the dynamic response of rock under blast loading with active confining pressure[J]. International Journal of Rock Mechanics and Mining Sciences, 2018, 102: 101-108. doi: 10.1016/j.ijrmms.2018.01.011

    [3]

    MITCHELL T R, WANG Z, ARAOS M, et al. Experimental and numerical investigation into the fracture patterns induced by blast-loading under unconfined and confined conditions[J]. Rock Mechanics and Rock Engineering, 2023, 56(4): 2433-2455. doi: 10.1007/s00603-022-03195-x

    [4]

    YANG L Y, YANG A Y, CHEN S Y, et al. Model experimental study on the effects of in situ stresses on pre-splitting blasting damage and strain development[J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 138: 104587. doi: 10.1016/j.ijrmms.2020.104587

    [5] 张凤鹏, 闫广亮, 郝琪琪, 等. 单向压应力对砂岩漏斗爆破过程的影响[J]. 东北大学学报(自然科学版), 2021, 42(2): 220-225.

    ZHANG Fengpeng, YAN Guangliang, HAO Qiqi, et al. Influence of unidirectional compressive stress on the crater blasting of sandstone[J]. Journal of Northeastern University (Natural Science), 2021, 42(2): 220-225. (in Chinese)

    [6] 徐颖, 顾柯柯, 葛进进, 等. 装药不耦合系数对初始地应力下岩石爆破裂纹扩展影响的试验研究[J]. 爆破, 2022, 39(4): 1-9.

    XU Ying, GU Keke, GE Jinjin, et al. Experimental study on effect of charge decoupling coefficient on crack propagation in rock by blasting under initial in-situ stress[J]. Blasting, 2022, 39(4): 1-9. (in Chinese)

    [7] 杨建华, 孙文彬, 姚池, 等. 高地应力岩体多孔爆破破岩机制[J]. 爆炸与冲击, 2020, 40(7): 118-127.

    YANG Jianhua, SUN Wenbin, YAO Chi, et al. Mechanism of rock fragmentation by multi-hole blasting in highly-stressed rock masses[J]. Explosion and Shock Waves, 2020, 40(7): 118-127. (in Chinese)

    [8]

    HAN H Y, FUKUDA D, XIE J B, et al. Rock dynamic fracture by destress blasting and application in controlling rockbursts in deep underground[J]. Computers and Geotechnics, 2023, 155: 105228. doi: 10.1016/j.compgeo.2022.105228

    [9]

    CÓRDOVA E, GOTTREUX I, ANANI A, et al. Blasting and preconditioning modelling in underground cave mines under high stress conditions[J]. Journal of the Southern African Institute of Mining and Metallurgy, 2021, 121(2): 71-80.

    [10] 高启栋, 卢文波, 冷振东, 等. 岩石爆破中孔内起爆位置对爆炸能量传输的调控作用研究[J]. 岩土工程学报, 2020, 42(11): 2050-2058. doi: 10.11779/CJGE202011010

    GAO Qidong, LU Wenbo, LENG Zhendong, et al. Regulating effect of detonator location in blast-holes on transmission of explosion energy in rock blasting[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(11): 2050-2058. (in Chinese) doi: 10.11779/CJGE202011010

    [11] 陈榕, 武智勇, 郝冬雪, 等. 高应力下石英砂三轴剪切颗粒破碎演化规律及影响[J]. 岩土工程学报, 2023, 45(8): 1713-1722.

    CHEN Rong, WU Zhiyong, HAO Dongxue, et al. Evolution rules and effects of particle breakage for quartz sand in triaxial shear tests under high pressures[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(8): 1713-1722. (in Chinese)

    [12]

    XIE L X, LU W B, ZHANG Q B, et al. Analysis of damage mechanisms and optimization of cut blasting design under high in situ stresses[J]. Tunnelling and Underground Space Technology, 2017, 66: 19-33. doi: 10.1016/j.tust.2017.03.009

    [13]

    LI X H, ZHU Z M, WANG M, et al. Numerical study on the behavior of blasting in deep rock masses[J]. Tunnelling and Underground Space Technology, 2021, 113: 103968. doi: 10.1016/j.tust.2021.103968

    [14]

    DEHGHAN BANADAKI M M, MOHANTY B. Numerical simulation of stress wave induced fractures in rock[J]. International Journal of Impact Engineering, 2012, 40: 16-25.

    [15] 杨仁树, 肖成龙, 陈程, 等. 基于分形理论不同装药量的爆破动焦散线实验研究[J]. 振动与冲击, 2020, 39(14): 80-86, 93.

    YANG Renshu, XIAO Chenglong, CHEN Cheng, et al. Experimental study on the blasting dynamic caustics under different charge weight based on the fractal theory[J]. Journal of Vibration and Shock, 2020, 39(14): 80-86, 93. (in Chinese)

    [16]

    LI X D, LIU K W, YANG J C, et al. Numerical study on blast-induced fragmentation in deep rock mass[J]. International Journal of Impact Engineering, 2022, 170: 104367. doi: 10.1016/j.ijimpeng.2022.104367

  • 期刊类型引用(1)

    1. 朱晟,陆雪妮,汤洪洁,孟涛,路德任. 大粒径砂砾料的缩尺试验方法及应用. 水利规划与设计. 2022(12): 5-10 . 百度学术

    其他类型引用(1)

图(18)  /  表(2)
计量
  • 文章访问数:  0
  • HTML全文浏览量:  0
  • PDF下载量:  0
  • 被引次数: 2
出版历程
  • 收稿日期:  2023-07-26
  • 网络出版日期:  2024-03-24
  • 刊出日期:  2024-09-30

目录

/

返回文章
返回