Processing math: 39%

竖向力-温度荷载-扭矩加载路径下单桩承载变形研究

江杰, 陈秋怡, 黄中正, 陈朝棋, 欧孝夺

江杰, 陈秋怡, 黄中正, 陈朝棋, 欧孝夺. 竖向力-温度荷载-扭矩加载路径下单桩承载变形研究[J]. 工程力学, 2023, 40(11): 130-139. DOI: 10.6052/j.issn.1000-4750.2022.01.0125
引用本文: 江杰, 陈秋怡, 黄中正, 陈朝棋, 欧孝夺. 竖向力-温度荷载-扭矩加载路径下单桩承载变形研究[J]. 工程力学, 2023, 40(11): 130-139. DOI: 10.6052/j.issn.1000-4750.2022.01.0125
JIANG Jie, CHEN Qiu-yi, HUANG Zhong-zheng, CHEN Chao-qi, OU Xiao-duo. RESEARCH ON BEARING DEFORMATION OF SINGLE PILE UNDER VERTICAL FORCE-THERMAL LOAD-TORQUE LOADING PATH[J]. Engineering Mechanics, 2023, 40(11): 130-139. DOI: 10.6052/j.issn.1000-4750.2022.01.0125
Citation: JIANG Jie, CHEN Qiu-yi, HUANG Zhong-zheng, CHEN Chao-qi, OU Xiao-duo. RESEARCH ON BEARING DEFORMATION OF SINGLE PILE UNDER VERTICAL FORCE-THERMAL LOAD-TORQUE LOADING PATH[J]. Engineering Mechanics, 2023, 40(11): 130-139. DOI: 10.6052/j.issn.1000-4750.2022.01.0125

竖向力-温度荷载-扭矩加载路径下单桩承载变形研究

基金项目: 国家自然科学基金项目(52068004);广西大学学科交叉科研项目(2022JCB012);广西重点研发计划项目(AB19245018)
详细信息
    作者简介:

    陈秋怡(1998−),女,广西人,硕士生,主要从事桩基工程理论研究(E-mail: chenqyqy@foxmail.com)

    黄中正(1996−),男,广西人,硕士生,主要从事桩基工程理论研究(E-mail: 1536215838@qq.com)

    陈朝棋(1997−),男,四川人,硕士生,主要从事桩基工程理论研究(E-mail: 1769794511@qq.com)

    欧孝夺(1970−),男(壮族),广西人,教授,博士,博导,主要从事特殊土的特性研究(E-mail: ouxiaoduo@163.com)

    通讯作者:

    江 杰(1979−),男,湖北人,研究员,博士,博导,主要从事桩基工程理论与应用研究(E-mail: jie_jiang001@126.com)

  • 中图分类号: TU473.1

RESEARCH ON BEARING DEFORMATION OF SINGLE PILE UNDER VERTICAL FORCE-THERMAL LOAD-TORQUE LOADING PATH

  • 摘要: 正常服役过程能量桩受复杂的力学行为,目前针对多向荷载作用下能量桩的承载变形特性研究相对较少。为探讨能量桩在竖向力、温度荷载和扭矩共同作用下的承载特性,依次施加桩顶竖向力、温度荷载和桩顶扭矩,通过考虑温度荷载对桩侧摩阻力及边界条件的影响,基于荷载传递法及边界元法构建了桩身位移控制方程,提出竖向力→温度荷载→扭矩加载路径下的能量桩承载变形特性的分析方法,通过与已有试验和ABAQUS有限元结果进行对比,均具有比较好的吻合度。研究表明,温度荷载会改变单桩的荷载传递特征,影响桩身轴力和桩侧摩阻力分布。竖向力→温度荷载→扭矩加载路径下,温度变化引起的附加荷载会导致单桩抗扭承载力降低。进一步的参数分析表明,增大竖向荷载会使桩侧极限环向摩阻力减小,导致能量桩单桩抗扭能力降低26.2%(75%PuPu为竖向极限荷载);随着长径比的增加,桩身承载力逐渐增大,可选取合适的长径比抵消温度荷载带来的影响;随着温度增量的增大,桩身变形量逐渐增大,其中桩身0.6L以上部分变形较大,因此作为能量桩使用的工程桩基需要对地基上部进行加固。
    Abstract: The energy pile exhibits complex mechanical behavior, and there is relatively little research on the bearing deformation characteristics of the energy pile under multi-direction loads. To investigate the bearing characteristics of energy piles under vertical load, thermal load, and torque, the thermal load and torque were applied after preloading the vertical load to the pile top. By considering the effect of thermal load on the shaft resistance and boundary conditions of the pile, the pile displacement control equation is derived based on the load transfer method and boundary element method, and the analysis method of the bearing deformation characteristics of energy pile under vertical force → thermal load → torque loading path is proposed. The proposed method is in good agreement with the existing experiments and ABAQUS finite element results. The results show that thermal loading changes the load transfer characteristics of the monopile and affects the axial force and lateral frictional resistance distribution of the pile. Under the vertical force → thermal load → torque loading path, the additional load caused by the temperature change leads to the reduction of the monopile torsional resistance. Further parametric analysis shows that increasing the vertical load decreases the ultimate circumferential frictional resistance on the pile side, resulting in a 26.2% reduction (75%Pu, where Pu is the ultimate load) in the torsional resistance of the energy pile monopile. With the increase in the length-diameter ratio, the bearing capacity gradually increases, and thus the effect of thermal load can be offset by selecting a suitable length-diameter ratio. With the increase of temperature increment, the deformation increases gradually, and the deformation of the upper 0.6L is more significant, so the engineering pile foundation used as an energy pile needs to be reinforced on the upper part.
  • 近年来,基于传统地源热泵技术发展起来的能量桩逐渐受到国内外学者的广泛关注[1-3]。对于高层建筑等大型结构物使用的能量桩而言,除了受到竖向荷载、温度荷载作用外,可能还会承受水平荷载、扭矩荷载和动荷载。鉴于问题的复杂性,本文先讨论竖向力、温度和扭矩作用下能量桩的承载变形研究,以期更好地推广能量桩的应用。

    对于竖向荷载下能量桩的受力变形特性,国内外学者进行了相关的试验。LALOUI等[4]和BOURNE-WEBB等[5]开展了能量桩现场原位实验,研究温度变化对桩身力学行为的影响;路宏伟等[6]和蒋刚等[7]开展竖向力和温度荷载作用下的现场试验,揭示了桩体温度变化和桩顶竖向荷载会引起能量桩荷载传递机制和桩身内力的变化;方鹏飞等[8] 结合现场试验探究升温工况下能量桩的工作机制,重点研究了能量桩桩身内力的变化规律;王成龙等[9]认为桩顶和桩端约束受温度荷载的影响,并通过开展相应的模型试验,分析不同约束条件下桩身位移和应力变化规律,并进一步研究了温度零点随约束不同的变化规律。鉴于试验费用高昂,试验条件难以控制,许多学者开始用理论方法研究能量桩的受力变形特性。KNELLWOLF等[10]将荷载传递法应用于能量桩,采用弹簧模拟温度荷载下上部结构与桩的相互作用,构造了竖向力和温度荷载作用下的荷载传递方程;PASTEN等[11]进一步将温度变形的影响考虑在桩段压缩量中,分析了温度荷载作用下能量桩位移变化规律;徐新丽等[12]基于弹性有限单元分析模型,通过对桩身和桩周土体进行受力分析,将温度荷载简化为单元内力中,最后结合静力平衡,建立了竖向力和温度荷载作用下能量桩桩身平衡方程,并结合London[4]试验和OUYANG等[13]提出的混合传递法进行了对比验证。费康等[14]引入双曲线荷载传递函数,将温度荷载作用考虑到单元压缩量中,结合单元增量平衡方程,揭示了任意竖向力和温度荷载作用下的桩身变形、桩身内力分布。

    综上,针对竖向力和温度荷载作用下能量桩承载变形特性,在试验和理论方面研究取得较多成果。但目前针对多向荷载作用下能量桩的研究仍较少,缺乏足够的理论研究,限制其进一步的推广应用。因此,本文对竖向力→温度荷载→扭矩加载路径下能量桩的承载变形展开研究。根据能量桩的工作条件,依次施加桩顶竖向力、温度荷载和桩顶扭矩。基于荷载传递法和边界单元法提出竖向力→温度荷载→扭矩加载路径下的能量桩单桩计算方法,并通过与已有试验和ABAQUS有限元结果进行验证,最后分析竖向荷载、长径比和温度对能量桩承载变形的影响。

    本文分析基于以下假定:

    1) 竖向力和温度荷载作用下桩的径向变形较小[15],故计算分析中只考虑桩的轴向变形。

    2) 桩为线弹性等截面桩,桩-土界面处的荷载传递函数均符合双曲线形式,如图1所示,桩侧和桩端荷载传递函数分别表示为:

    τ = Δa+bΔ (1)
    τ = Δba+bΔb (2)

    式中:ττ分别为桩侧剪应力和桩端阻力;ΔΔb分别为桩-土相对位移和桩端位移;ab、aʹbʹ均为模型参数。对桩侧,根据CAO等[16]建议确定,a=Δu/τub=Rsf/τu,其中:Δu为桩体达极限摩阻力对应的极限位移,FLEMING等[17]建议可取0.5%~2.0%DD为桩的直径;Rsf为桩侧破坏比,根据费康等[11]建议取为1.0;τu为桩-土界面的破坏剪应力,τu=K0σtanφ,其中K0为静止土压力系数,K0=1sinφφ为桩-土界面摩擦角; σ 为竖向应力。对桩端,根据RANDOLPH[18]建议确定,a=πD(1ν)/(8G)b=Rbf/τu,其中:Gs为桩端土的剪切模量;ν为桩端土的泊松比;Rbf为桩端破坏比,取为1.0;τu为桩端界面阻力破坏值,根据极限平衡理论确定,τu=ζcNc+ζqγhNq,其中ζcζq为桩端为方形、圆形时的形状系数,NcNq为承载力系数。

    3) 桩周土体性质、桩侧和桩端荷载传递函数不受温度影响。

    4) 取温度增加、桩体位移向下和桩侧阻力向上为正。

    图  1  能量桩桩-土界面模型
    Figure  1.  Pile-soil interface model of energy pile

    将桩长为L的能量桩划分为N个单元,对单元和节点自上向下编号,建立竖向力和温度荷载作用下单桩承载变形特性的简化分析方法,如图2所示。

    图  2  能量桩沿桩长离散示意图
    Figure  2.  Discrete diagram of energy pile along pile length

    取任意桩体单元分析,由增量平衡条件,有:

    ΔPi=ΔPi+1+ΔQi (3)
    ΔQi=ksiπDΔLiΔXi+ΔXi+12 (4)

    式中:ΔPi为单元i上截面的轴力增量;ΔPi+1为单元i下截面的轴力增量;ΔQi为单元i桩侧阻力的增量;ksi为单元i中点处桩侧抗剪切刚度,即荷载传递函数的切线斜率;ΔLi为单元i长度;D为桩体直径;ΔXi为节点截面i的位移增量;ΔXi+1为节点截面i+1的位移增量。

    ΔXi=ΔXi+1+Δsi (5)

    式中,Δsi为单元的压缩量。考虑竖向力和温度荷载的相互作用,有:

    Δsi=ΔPi+ΔPi+12EAΔLiΔT (6)

    式中:E为桩身弹性模量;A为单元截面面积;ΔT为自由膨胀量,根据线热膨胀理论,ΔT=αTΔTΔLiαT为桩身混凝土线膨胀系数,ΔT为温度增量。

    联立式(3)~式(6),分离变量,有:

    ΔPi=(EAΔLi+πDksiΔLi4)ΔXi(EAΔLiπDksiΔLi4)ΔXi+1+EAαTΔT (7)
    ΔPi+1=(EAΔLiπDksiΔLi4)ΔXi(EAΔLi+πDksiΔLi4)ΔXi+1+EAαTΔT (8)

    写成矩阵形式有:

    \left\{ \begin{matrix} \Delta {P_i} \\ \Delta {P_{i + 1}} \end{matrix} \right\} = \left[ \begin{matrix} k^i_{11} & - k^i_{12} \\ k^i_{21} & - k^i_{22} \end{matrix} \right]\left\{ \begin{matrix} \Delta {X_i} \\ \Delta {X_{i + 1}} \end{matrix} \right\} + {\boldsymbol{C}} (9)

    式中:{{\boldsymbol{k}}^i}为桩身单元i刚度矩阵。其中:

    k_{11}^i = k_{22}^i = \frac{{EA}}{{\Delta {L_i}}} + \frac{{\pi D{k_{si}} \cdot \Delta {L_i}}}{4}\;,
    k_{12}^i = k_{21}^i = \frac{{EA}}{{\Delta {L_i}}} - \frac{{\pi D{k_{si}} \cdot \Delta {L_i}}}{4}\;,
    C = EA{\alpha _T}\Delta T 。

    结合单元上、下截面连续性,将每个桩身单元的刚度矩阵进行组装,得到竖向力和温度荷载作用下桩身整体矩阵通式为:

    {{\boldsymbol{K}}^i} = \left[ \begin{matrix} - (k_{11}^i + k_{22}^{i - 1}) & k_{12}^i \\ k_{21}^i & - (k_{22}^i + k_{11}^{i + 1}) \end{matrix}\right] (10)

    接着考虑桩顶边界第一个单元的轴力增量及桩端边界第N+1节点反力得到式(11)、式(12)。

    桩顶节点1:

    k_{11}^1 \cdot \Delta {X_1} - k_{12}^1 \cdot \Delta {X_2} = \Delta {F_{}} - C (11)

    桩端节点N+1:

    \Delta {P_{N + 1}} = {k_{\text{b}}}A \cdot \Delta {X_{N + 1}} (12)

    式中, {k_{\text{b}}} 为桩端土体的抗压刚度系数。将式(12)代入式(8),有:

    k_{21}^N \cdot \Delta {X_N} - (k_{22}^N + {k_{\text{b}}}{A_{}})\Delta {X_{N + 1}} = - C (13)

    综上,可得到竖向力和温度荷载作用下能量桩的控制方程:

    \Delta {\boldsymbol{F}} = {\boldsymbol{K}} {\Delta {\boldsymbol{X}}} (14)

    其中,桩的总刚度矩阵 {\boldsymbol{K}} 为:

    {\boldsymbol{K}} = {\left[ {\begin{array}{*{20}{l}} {k_{11}^1}\qquad\;\;\;\;\;\;\;{-k_{12}^1}\;\qquad{}\qquad{}\qquad{}\\ {k_{21}^1}\qquad{ - \left( {k_{22}^1 + k_{11}^2} \right)}\qquad{k_{12}^2}\qquad{}\qquad{}\\ \quad \ddots \qquad\qquad \ddots \qquad\qquad \ddots \qquad{}\qquad{}\\ {}\qquad{k_{21}^i}\qquad{ - \left( {k_{22}^i + k_{11}^{i + 1}} \right)}\qquad{k_{12}^{i + 1}}\qquad{}\\ {}\qquad\qquad \ddots\qquad \qquad \ddots \qquad\qquad \ddots \qquad{}\\ {}\qquad\qquad{}\qquad\qquad{}\quad\qquad{k_{21}^N}\qquad{ - \left( {k_{22}^N + {k_{\rm b}}A} \right)} \end{array}} \right]_{\left( {N + 1} \right) \times \left( {N + 1} \right)}}

    对方程组求解可得到桩体节点位移增量 \Delta X 。求解过程如下:① 将荷载分为多级荷载增量,从第一个荷载增量 {\Delta {\boldsymbol{F}} _1} ,先施加荷载增量的1/2,即 {\Delta {\boldsymbol{F}} _{1/2}} ,根据桩侧、桩端初始刚度建立桩身刚度矩阵,分别求出 {\Delta {\boldsymbol{F}} _{1/2}} 的节点位移;② 根据上一增量结束阶段的节点位移、应力建立刚度矩阵,求出增量中点位移;③ 用新求出的桩身节点节点位移求出增量 {\Delta {\boldsymbol{F}} _1} 的中点刚度矩阵,接着由式(14)求出第一级增量下的位移、应力;④ 重复步骤②~步骤③,直到最后一个增量,求出桩体节点位移。

    根据能量桩的实际运行情况,依次进行竖向力加载和温度加载。竖向力加载中桩顶视为自由,温度荷载加载过程中上部结构对桩顶的约束作用,采用刚度系数为 {k_{\text{t}}} 对节点1的刚度矩阵进行修正,得到式(15):

    (k_{11}^1 + {k_{\text{t}}}{A_{}})\Delta {X_1} - k_{12}^1\Delta {X_2} = - C (15)

    为求解竖向力→温度荷载→扭矩加载路径能量桩承载变形特性,根据江杰等 [19]和邹新军等[20],通过边界元法,单元 i 表面作用有竖向力、温度产生的桩侧竖向摩阻力 {\tau _{\text{v}}}(i) 和扭矩作用下的环向摩阻力 {\tau _{\text{t}}}(i) ,如图3所示,假设桩身摩阻力满足式(16),有:

    图  3  桩周土体剪切作用
    Figure  3.  Shear action of soil around piles
    \tau (i) = \sqrt {[ {{{[{\tau _{\text{t}}}(i)]}^2} + {{[{\tau _{\text{v}}}(i)]}^2}} ]} \leqslant {\tau _{\text{f}}}(i) (16)

    竖向力→温度加载结束后,扭矩作用下的环向极限摩阻力为:

    {\tau _{{\text{tf}}}}(i) = \sqrt {[ {{{[{\tau _{\text{f}}}(i)]}^2} - {{[{\tau _{\text{v}}}(i)]}^2}} ]} (17)

    式中:{\tau _{\text{f}}}(i)为单元 i 的极限摩阻力; {\tau _{\text{v}}}(i) 为单元 i 的竖向摩阻力;{\tau _{{\text{tf}}}}(i)为单元 i 的环向极限摩阻力。

    扭矩作用下,桩单元的静力平衡条件为:

    \frac{{{\rm d}T({\textit{z}})}}{{{\rm d}{\textit{z}}}} = 0.5\pi {\tau _{\text{t}}}({\textit{z}}){D^2} (18)

    式中:T(\textit{z})\textit{z}处扭矩; {\tau _{\text{t}}}(i) 为深度\textit{z}处桩侧环向摩阻力。

    扭矩作用下桩的变形条件为:

    \frac{{{\rm d}\theta ({\textit{z}})}}{{{\rm d}{\textit{z}}}} = - \frac{{T({\textit{z}})}}{{{J_{\text{p}}}{G_{\text{p}}}}} (19)

    式中:\theta (\textit{z})为扭转角; {G_{\text{p}}} 为桩单元的剪切模量; {J_{\text{p}}} 为桩单元的极惯性矩。

    根据江杰等[21],扭转角与扭转位移的关系可以表示为:

    \theta ({\textit{z}}) = \frac{{2{s_{\text{t}}}({\textit{z}})}}{D} (20)

    式中,{s_{{\rm t}}}({\textit z})为深度{\textit z}处的扭转位移。

    \frac{{{\rm d}{\theta ^2}({\textit{z}})}}{{{\rm d}{{\textit{z}}^2}}} = {\alpha ^2}\theta ({\textit{z}}) (21)

    式中,\alpha {\text{ = }}\sqrt {\dfrac{{{k_\theta }\pi {D^3}}}{{4{G_{\text{p}}}{J_{\text{p}}}}}}。其中, {k_\theta } 可以表示为:

    {k_\theta }({\textit{z}}) = \frac{{{\tau _{\text{t}}}({\textit{z}})}}{{{s_{\text{t}}}({\textit{z}})}} (22)

    桩顶、桩端边界条件分别如下:

    桩顶:

    T(0) = {T_1} (23)

    桩端边界条件引入POULOS[22]的建议,有:

    {\theta _n} = \frac{{{\theta _{n + 1}} - {\theta _{n - 1}}}}{2} = \frac{{3{T_{\text{b}}}}}{{16{{G_{\text{s}}'}}r_0^3}} (24)

    式中: {T_1} 为桩顶扭矩; {\theta _i} 为第 i 单元的桩身扭转角( i = n - 1,{\text{ }}n,{\text{ }}n + 1 );{T_{\text{b}}} = {T_1} - \displaystyle\sum\limits_{i = 0}^{n - 1} {0.5\pi {\tau _{\text{t}}}} (i){D^2}h为桩端扭矩。

    \left[ \begin{matrix} {{B_0}}&2&{}&{}&{}\\ 1&{{B_1}}&1&{}&{}\\ \ddots & \ddots & \ddots &{}&{}\\ {}&1&{{B_i}}&1&{}\\ {}& \ddots & \ddots & \ddots &{}\\ {}&{}&1&{{B_{n - 1}}}&1\\ {}&{}&{}&2&{{B_n}} \end{matrix} \right]\left[ \begin{matrix} {{\theta _0}}\\ {{\theta _1}}\\ \vdots \\ {{\theta _i}}\\ \vdots \\ {{\theta _{n - 1}}}\\ {{\theta _n}} \end{matrix} \right] = \left[ \begin{matrix} {\dfrac{{ - 2{T_1}h}}{{{G_{\rm p}}{J_{\rm p}}}}}\\ 0\\ \vdots \\ 0\\ \vdots \\ 0\\ {\dfrac{{ - 3{T_{\rm b}}h}}{{8G'r_0^3}}} \end{matrix} \right] (25)

    式中,{B_i} = - (\alpha _i^2{h^2} + 2)

    采用有限差分进行求解:① 假设桩身扭转角\theta 为任意非零矩阵,根据式(17)、式(20)和式(22)求得桩身扭转矩阵;② 由式(25)求得到桩身扭转角{\theta ^k};③ 用{\theta ^k}求出新的桩身扭转刚度矩阵,由式(25)求得到新的{\theta ^{k + 1}}。当|{\theta ^k} - {\theta ^k}^{ + 1}|小于限定值时,可得到桩身扭转角\theta ,若|{\theta ^k} - {\theta ^k}^{ + 1}|不满足要求则重复步骤②~步骤③直至满足限定值。

    本文计算依次进行竖向力加载、温度荷载加载和扭矩加载,有:

    1) 竖向力加载

    \Delta T = 0,将竖向力划分为若干增量步,根据土层参数和桩身参数得到初始刚度矩阵,并采用中点增量法对式(14)进行求解,得到桩体节点位移增量。

    2) 温度荷载加载

    以竖向力加载结束的桩身节点位移作为温度加载的初始加载状态,将温度荷载划分为若干增量步,求解过程与竖向力加载大致相同,最后得到竖向力→温度加载下的桩基响应。

    3) 扭矩加载

    计算竖向力加载和温度荷载加载结束后的桩侧竖向摩阻力,通过式(17)计算桩侧环向极限摩阻力,取深度z处的微段桩段{\rm d}{\textit z}进行分析,建立桩身扭矩角整体方程,利用有限差分法,对式 (25)进行求解,通过控制|{\theta ^k} - {\theta ^{k+1}}{ }|迭代误差得到桩身节点扭转角,最后输出竖向力→温度荷载→扭矩加载路径下能量桩的桩基响应。

    1) 模型简介

    本文以文献[23]中的能量桩为例,桩径D=0.88 m,桩长L=19.6 m。通过ABAQUS有限元软件进行模拟,在桩-土界面设置接触对模拟实际桩-土界面的相互作用。桩-土界面的摩擦角根据试验取为25°,桩身和桩周砂土分别采用弹性本构模型和Mohr-Coulomb理想弹塑性模型。桩-土其他材料参数见表1

    表  1  桩-土材料参数
    Table  1.  Pile and soil parameters
    参数砂土混凝土
    密度ρ/(kg/m3)15712640
    弹性模量E/MPa 1127.8×103
    泊松比ν0.20.25
    热膨胀系数α/(με/℃)102.2
    下载: 导出CSV 
    | 显示表格

    2) 边界条件及网格划分

    根据工程实际,土体侧面约束径向位移,底部和顶部视为固定边界和自由面边界条件。

    通过自动平衡法来实现地应力平衡,在地应力平衡分析步,仅对完整的土体模型进行地应力平衡;在静力分析步,将桩单元激活,从而实现地应力平衡。在静力分析步施加竖向力,竖向力加载结束后,导入温度场模拟能量桩受热状态,最后达到指定温度后在桩顶施加扭矩。桩体和土体均采用C3D8网格类型,桩体模型通过结构化网格技术进行划分,总共5760个桩体网格单元;土体模型首先对土体定义Partition,最后结合扫掠技术划分,总共45 000个土体网格单元,如图4所示。

    图  4  三维数值模型
    Figure  4.  3D numerical model

    3) 桩身轴力分布

    纯升温工况下,桩身的轴力在数值上沿深度呈现先增大、后减小的规律,大致在桩身0.6L处达到最大值。图5为不同温度增量下沿桩长方向的桩身轴力图,随着温度增量的增大,桩身轴力逐渐增大,其中桩顶和桩端的轴力较小,桩身中下部轴力较大。可以看出,本文方法计算结果、有限元结果与NG等 [23]的试验数据都具有较好的一致性。

    图  5  温度荷载下桩身轴力分布
    Figure  5.  Axial force distribution of pile under thermal load

    4) 桩侧竖向摩阻力分布

    纯升温工况下,桩身体积变化导致与桩周土体将产生相对位移,从而桩侧竖向摩阻力发生改变。图6为不同温度增量作用下沿桩长方向的桩身竖向摩阻力分布图。

    图6可知,计算值与模拟值有很好的一致性,桩侧摩阻力沿深度的变化规律与龚建清等[24]的模拟趋势一致。中性点(桩侧摩阻力为0的点)位置大约在桩长0.6L处,与BOURNE-WEBB等[5]现场实测的中性点位置相近,桩侧竖向摩阻力数值上沿桩长方向呈现上部分先增大后减小至0,下部分从0开始反向增大的分布,这与LAOUI等[4]研究成果一致。本文算例工况采用NG等[23]的进行计算,中密砂土层,弹性模量较小,对桩端的约束较小,故中性点位置在0.6L处,而LAOUI等[4]的桩端持力层刚度大,可提供更大的桩端约束,因此其中性点更接近桩端。此外,可以看出,纯温度荷载作用下(\Delta T > 0),导致桩身上半部分的摩阻力为负值,\Delta T = 30 ℃相对于\Delta T = 15 ℃,能量桩最大侧摩阻力(绝对值)有所增大。

    图  6  温度荷载下桩侧摩阻力分布
    Figure  6.  Distribution of pile shaft resistance under thermal load

    将纯温度加载工况下的计算值与相关试验数据和有限元对比,证明了方法的准确性。

    本文的计算理论不仅可以分析竖向力→温度荷载→扭矩加载路径下的能量桩承载变形特性,也可以用来分析竖向力→温度荷载加载路径下能量桩计算。根据蒋刚等[7]试验参数,重新建立数值模型,将竖向力→温度荷载加载路径下单桩计算值、模拟值与蒋刚等[7]试验值进行对比。图7为桩身的轴力随深度变化曲线,对比单一竖向力加载、竖向力→温度荷载工况可知,温度荷载会改变桩身轴力分布。由图7可知,本文计算值、模拟值与试验值接近,因此可证明本文计算方法的正确性。

    图  7  竖向力→温度荷载加载路径桩身轴力分布曲线
    Figure  7.  Axial force distribution curve under vertical force → thermal loading path

    本文基于NG等[23]的试验,对温度增量为30 ℃的单桩施加扭矩,通过与有限元进行对比,验证竖向力→温度荷载→扭矩加载路径下单桩承载变形计算方法的正确性。

    为研究竖向力→温度荷载→扭矩加载路径下,桩顶扭矩对能量桩环向摩阻力的影响。对竖向力→温度荷载结束后施加扭矩后,为了验证本文理论方法的正确性,通过有限元分别对桩顶施加扭矩100 kN·m、200 kN·m、300 kN·m、400 kN·m,提取扭矩作用后的桩身环向摩阻力后与计算模拟结果对比如图8所示,4组计算值与数值模拟结果均比较接近,证明了本文方法在计算竖向力→温度荷载→扭矩加载路径下桩侧环向摩阻力的正确性。

    图8可知,当桩顶扭矩分别为100 kN·m、200 kN·m、300 kN·m、400 kN·m时,桩身环向摩阻力随深度呈现线性增大的趋势,总体上,竖向力→温度荷载→扭矩加载路径下桩身环向摩阻力会随着桩顶扭矩的增大而增大,当桩顶扭矩增大到一定范围后,环向摩阻力的增长较之前缓慢。原因是,扭矩较小时,桩-土界面进入弹性阶段,温度荷载引起土体的抵抗作用发挥明显,从而导致环向摩阻力随深度增大;当扭矩增大到一定数值时,环向摩阻力趋势趋于一致,表明桩顶扭矩起主导作用。

    图  8  桩侧环向摩阻力分布
    Figure  8.  Distribution of circumferential friction resistance on pile side

    对于竖向力→温度荷载→扭矩加载路径下单桩计算,通过确定桩侧极限环向摩阻力,并通过式(24)考虑桩端边界条件,最后得到T-θ曲线。取竖向荷载P=1000 kN时,研究扭转角随不同桩顶扭矩作用下的变化。整体上来看,本文计算得到的曲线和有限元结果较为吻合,证明了竖向力→温度荷载→扭矩加载路径下单桩承载变形计算方法具有较好的准确性。从图9可看出,ΔT=30 ℃工况下单桩的极限扭矩[25]为645 kN·m,对比ΔT=0 ℃工况可知,能量桩的极限扭矩降低了20.37%,说明温度荷载对能量桩的极限扭矩有影响。

    图  9  T{{ - }}\theta曲线对比
    Figure  9.  Comparison of T{{ - }}\theta curves

    为研究竖向力→温度荷载→扭矩加载路径下单桩的影响因素,分别对竖向力、长径比和温度增量进行分析。

    为研究不同竖向荷载(根据NG等[23]的试验,取温度荷载增量为30 ℃)对能量桩荷载传递特性的影响,分别施加25% Pu、50% Pu、75% Pu的竖向力(Pu为竖向极限荷载)。图10给出3种荷载工况下,桩顶扭矩角随扭矩的发展情况。可见,温度增量和桩顶扭矩为定值,随着竖向力的增加,25% Pu、50% Pu、75% Pu对应的极限扭矩分别为680 kN·m、590 kN·m和515 kN·m,极限扭矩降低幅度为26.2%,说明了竖向力→温度荷载→扭矩加载路径下,在达到竖向极限荷载的过程中,由竖向力产生的桩侧竖向摩阻力会不断逼近极限值,进而影响桩侧极限环向摩阻力,最终导致单桩极限扭矩受到影响,此结果验证了式(17)的正确性。

    图  10  三种竖向力工况下T{{ - }}\theta曲线对比
    Figure  10.  Comparison of T{{ - }}\theta curves under three vertical load conditions

    图11为25%Pu、50%Pu、75%Pu(温度增量30 ℃,扭矩350 kN·m)三种工况下的桩身扭转角随深度的分布情况。由图11可知,三种竖向力工况下桩身扭转角存在较大的差异。25%Pu、50%Pu和75%Pu作用下桩顶桩身扭转角分别为0.020、0.023和0.027,而桩顶扭转角最大增幅约35%。因此,对于竖向力→温度荷载→扭矩加载路径下,随着竖向力的增大,桩身变形会受到较大的影响。因此工程桩基中的能量桩需适当控制桩顶竖向力水平。

    为探究桩身长径比对承载力的影响规律,选取3组不同长径比的桩进行分析,获得相应的桩身承载力(桩顶所能承受竖向力与扭矩极值)包络图如图12所示。

    图  11  三种竖向力工况下桩身扭转角分布情况
    Figure  11.  Distribution of pile torque angle under three vertical load conditions
    图  12  不同长径比L/D对桩顶承载力包络线的影响
    Figure  12.  Effect of different ratio L/D on the envelopes of pile top bearing capacity

    由图12可以看出,相同温度增量下,当长径比L/D=10时(保持桩径不变,改变桩身长度),桩身承载能力较小;当长径比L/D=30时,对应的扭矩极限承载力相对于L/D=20增加约 1.1倍,表明随着桩身长径比的增加,桩身承载力得到提高,承载力包络线有所外扩。

    为研究温度增量对能量桩承载变形特性的影响。图13给出了温度增量分别为10 ℃、20 ℃、30 ℃、40 ℃工况下,桩顶扭矩-扭矩角分布情况。可以看出,增大温度增量,单桩抗扭能力会受到影响。

    图  13  不同温度增量下T{\text{ - }}\theta 曲线对比
    Figure  13.  Comparison of T{\text{ - }}\theta curves under different temperature increments

    图14为不同温度增量(扭矩为500 kN·m)工况下的桩身扭转角分布情况。由图14可知,四种温度工况下的沿深度分布的桩身扭转角存在较大的差异。在竖向力→温度荷载→扭矩加载路径下,随着温度增量的增大,桩身变形会受到影响。对比温度增量10 ℃、20 ℃、30 ℃、40 ℃的桩身扭转角沿深度的分布规律,可知桩端扭转角受温度影响较大,说明升温引起的桩身体积变化,会导致桩端阻力受影响。此外还可以看出,大约0.6L以上扭转角较大,0.6L以下扭转角较小,原因是随着温度荷载的增大,桩端约束有所提高,故靠近桩端附近的变形较小。因此,竖向力→温度荷载→扭矩加载路径下的单桩变形主要发生在桩身0.6L以上,故在能量桩施工过程中要注意对浅层地基的保护。

    图  14  不同温度增量下桩身扭转角分布情况
    Figure  14.  Distribution of torsion angle of pile under different temperature increments

    本文基于荷载传递法,将温度荷载考虑到单元压缩量中,建立了竖向力→温度荷载→扭矩加载路径下能量桩承载特性的分析方法,并对桩身内力及变形进行了分析,主要结论如下:

    (1) 温度荷载作用下(\Delta T > 0),中性点大致位于桩身0.6L处,桩端轴力数值较小,桩身中部轴力数值较大。

    (2) 竖向力→温度荷载→扭矩加载路径下,能量桩相对于传统桩的极限扭矩降低20.37%,说明温度荷载对单桩极限扭矩产生影响。工程设计可适当提高配筋率,以弱化温度效应的影响。

    (3) 对竖向力的参数分析表明,改变竖向力会对能量桩的抗扭承载力和扭转角产生较大影响。由竖向力产生的桩侧竖向摩阻力会影响桩侧环向极限摩阻力,导致单桩抗扭能力降低26.2% (75%Pu),桩身变形量增大35% (75%Pu)。

    (4) 对桩身长径比参数分析表明,竖向力→温度荷载→扭矩加载路径下单桩扭矩极限承载力随着桩身长径比(保持桩径不变)增加有所提高,因此可通过适当提高长径比来提高桩-土极限承载力。

    (5) 对温度增量的参数分析表明,改变温度增量会对能量桩的承载力和变形产生影响,且单桩变形主要发生在桩身0.6L以上,桩端附近的变形较小。因此在实际能量桩工程中,需要考虑浅层地基的保护。

  • 图  1   能量桩桩-土界面模型

    Figure  1.   Pile-soil interface model of energy pile

    图  2   能量桩沿桩长离散示意图

    Figure  2.   Discrete diagram of energy pile along pile length

    图  3   桩周土体剪切作用

    Figure  3.   Shear action of soil around piles

    图  4   三维数值模型

    Figure  4.   3D numerical model

    图  5   温度荷载下桩身轴力分布

    Figure  5.   Axial force distribution of pile under thermal load

    图  6   温度荷载下桩侧摩阻力分布

    Figure  6.   Distribution of pile shaft resistance under thermal load

    图  7   竖向力→温度荷载加载路径桩身轴力分布曲线

    Figure  7.   Axial force distribution curve under vertical force → thermal loading path

    图  8   桩侧环向摩阻力分布

    Figure  8.   Distribution of circumferential friction resistance on pile side

    图  9   T{{ - }}\theta曲线对比

    Figure  9.   Comparison of T{{ - }}\theta curves

    图  10   三种竖向力工况下T{{ - }}\theta曲线对比

    Figure  10.   Comparison of T{{ - }}\theta curves under three vertical load conditions

    图  11   三种竖向力工况下桩身扭转角分布情况

    Figure  11.   Distribution of pile torque angle under three vertical load conditions

    图  12   不同长径比L/D对桩顶承载力包络线的影响

    Figure  12.   Effect of different ratio L/D on the envelopes of pile top bearing capacity

    图  13   不同温度增量下T{\text{ - }}\theta 曲线对比

    Figure  13.   Comparison of T{\text{ - }}\theta curves under different temperature increments

    图  14   不同温度增量下桩身扭转角分布情况

    Figure  14.   Distribution of torsion angle of pile under different temperature increments

    表  1   桩-土材料参数

    Table  1   Pile and soil parameters

    参数砂土混凝土
    密度ρ/(kg/m3)15712640
    弹性模量E/MPa 1127.8×103
    泊松比ν0.20.25
    热膨胀系数α/(με/℃)102.2
    下载: 导出CSV
  • [1]

    GUO Y, ZHANG G Z, LIU S Y. Investigation on the thermal response of full-scale PHC energy pile and ground temperature in multi-layer strata [J]. Applied Thermal Engineering, 2018, 143: 836 − 848.

    [2]

    LORIA A F R, VADROT A, LALOUI L. Analysis of the vertical displacement of energy pile groups [J]. Geomechanics for Energy and the Environment, 2018, 16: 1 − 14. doi: 10.1016/j.gete.2018.04.001

    [3] 陆浩杰, 吴迪, 孔纲强, 等. 循环温度作用下饱和黏土中摩擦型桩变形特性研究[J]. 工程力学, 2020, 37(5): 156 − 165. doi: 10.6052/j.issn.1000-4750.2019.07.0353

    LU Haojie, WU Di, KONG Gangqiang, et al. Displacement characteristics of friction piles embedded in saturated clay subjected to thermal cycles [J]. Engineering Mechanics, 2020, 37(5): 156 − 165. (in Chinese) doi: 10.6052/j.issn.1000-4750.2019.07.0353

    [4]

    LALOUI L, NUTH M, VULLIET L. Experimental and numerical investigations of the behaviour of a heat exchanger pile [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2006, 30(8): 763 − 781.

    [5]

    BOURNE-WEBB P J, AMATYA B, SOGA K, et al. Energy pile test at Lambeth College, London: Geotechnical and thermodynamic aspects of pile response to heat cycles [J]. Géotechnique, 2009, 59(3): 237 − 248.

    [6] 路宏伟, 蒋刚, 王昊, 等. 摩擦型能源桩荷载–温度现场联合测试与承载性状分析[J]. 岩土工程学报, 2017, 39(2): 334 − 342.

    LU Hongwei, JIANG Gang, WANG Hao, et al. In-situ tests and thermo-mechanical bearing characteristics of friction geothermal energy piles [J]. Chinese Journal of Geotechnical Engineering, 2017, 39(2): 334 − 342. (in Chinese)

    [7] 蒋刚, 李仁飞, 王昊, 等. 摩擦型能源桩热–力耦合全过程承载性能分析[J]. 岩石力学与工程学报, 2019, 38(12): 2525 − 2534.

    JIANG Gang, LI Renfei, WANG Hao, et al. Numerical analysis of the bearing capacity of floating energy piles during the full process of thermal-mechanical coupling [J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(12): 2525 − 2534. (in Chinese)

    [8] 方鹏飞, 高翔, 娄扬, 等. 夏季工况下正常服役地热能源桩承载性能原位试验研究[J]. 岩石力学与工程学报, 2021, 40(5): 1032 − 1042.

    FANG Pengfei, GAO Xiang, LOU Yang, et al. Field test on the bearing behaviors of geothermal energy piles in natural service under the summer condition [J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(5): 1032 − 1042. (in Chinese)

    [9] 王成龙, 刘汉龙, 孔纲强, 等. 工作荷载下温度循环对桩基变形与应力的影响分析[J]. 岩土力学, 2016, 37(S1): 317 − 322. doi: 10.16285/j.rsm.2016.S1.042

    WANG Chenglong, LIU Hanlong, KONG Gangqiang, et al. Influence of circular temperature on the strain and stress of energy piles under a working load [J]. Rock and Soil Mechanics, 2016, 37(S1): 317 − 322. (in Chinese) doi: 10.16285/j.rsm.2016.S1.042

    [10]

    KNELLWOLF C, PERON H, LALOUI L. Geotechnical analysis of heat exchanger piles [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2011, 137(10): 890 − 902. doi: 10.1061/(ASCE)GT.1943-5606.0000513

    [11]

    PASTEN C, SANTAMARINA J C. Thermally induced long-term displacement of thermoactive piles [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2014, 140(5): 70 − 75.

    [12] 徐新丽, 蒋刚, 路宏伟, 等. 能源桩热-力半耦合弹性理论分析方法[J]. 南京工业大学学报(自然科学版), 2019, 41(1): 121 − 128.

    XU Xinli, JIANG Gang, LU Hongwei, et al. Elasticity theory of energy pile under thermal-mechanical semi-coupling [J]. Journal of Nanjing University of Technology (Natural Science Edition), 2019, 41(1): 121 − 128. (in Chinese)

    [13]

    OUYANG Y, SOGA K, LEUNG Y F. Numerical back-analysis of energy pile test at Lambeth College, London [M]// Geo-frontiers Congress 2011: Advances in Geotechnical Engineering, 2011: 440 − 449.

    [14] 费康, 戴迪, 洪伟. 能量桩单桩工作特性简化分析方法[J]. 岩土力学, 2019, 40(1): 70 − 80, 90. doi: 10.16285/j.rsm.2017.2266

    FEI Kang, DAI Di, HONG Wei. A simplified method for working performance analysis of single energy piles [J]. Rock and Soil Mechanics, 2019, 40(1): 70 − 80, 90. (in Chinese) doi: 10.16285/j.rsm.2017.2266

    [15]

    FAIZAL M, BOUAZZA A, HABERFIELD C, et al. Axial and radial thermal responses of a field-scale energy pile under monotonic and cyclic temperature changes [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2018, 144(10): 04018072. doi: 10.1061/(ASCE)GT.1943-5606.0001952

    [16]

    CAO W P, CHEN Y M, WOLFE W E. New load transfer hyperbolic model for pile-soil interface and negative skin friction on single piles embedded in soft soils [J]. International Journal of Geomechanics, 2014, 14(1): 92 − 100. doi: 10.1061/(ASCE)GM.1943-5622.0000289

    [17]

    FLEMING W, WELTMAN A J, RANDOLPH M F, et al. Piling engineering [M]. New York: CRC Press, 2008.

    [18]

    RANDOLPH M F, WROTH C P. An analysis of the vertical deformation of pile groups [J]. Géotechnique, 1979, 29(4): 423 − 439.

    [19] 江杰, 王顺苇, 欧孝夺, 等. 黏土地基中桩顶扭矩-竖向荷载加载路径下单桩承载特性分析[J]. 岩土力学, 2020, 41(11): 3573 − 3582.

    JIANG Jie, WANG Shunwei, OU Xiaoduo, et al. Analysis of the bearing characteristics of single pile under the T→V loading path in clay ground [J]. Rock and Soil Mechanics, 2020, 41(11): 3573 − 3582. (in Chinese)

    [20] 邹新军, 王亚雄, 徐洞斌. Gibson地基中V-T联合受荷桩承载力分析[J]. 工程力学, 2020, 41(11): 3573 − 3582. doi: 10.6052/j.issn.1000-4750.2014.01.0073

    ZOU Xinjun, WANG Yaxiong, XU Dongbin. Bearing capacity analysis of piles under V-T combined loading in Gibson subsoil [J]. Engineering Mechanics, 2020, 41(11): 3573 − 3582. (in Chinese) doi: 10.6052/j.issn.1000-4750.2014.01.0073

    [21] 江杰, 王顺苇, 欧孝夺, 等. 膨胀土地基中单桩受扭非线性分析[J]. 工程力学, 2020, 37(11): 219 − 227. doi: 10.6052/j.issn.1000-4750.2020.04.0243

    JIANG Jie, WANG Shunwei, OU Xiaoduo, et al. Nonlinear torsional analysis of a single pile in an expansive soil foundation [J]. Engineering Mechanics, 2020, 37(11): 219 − 227. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.04.0243

    [22]

    POULOS H G. Torsional response of piles [J]. Journal of the Geotechnical Engineering Division, 1975, 101(10): 1019 − 1035. doi: 10.1061/AJGEB6.0000203

    [23]

    NG C W W, SHI C, GUNAWAN A, et al. Centrifuge modelling of heating effects on energy pile performance in saturated sand [J]. Canadian Geotechnical Journal, 2015, 52(8): 1045-1057.

    [24] 龚建清, 彭文哲. 倾斜荷载下能量桩受力变形特性三维有限元分析[J]. 岩土工程学报, 2021, 43(11): 2105 − 2111.

    GONG Jianqing, PENG Wenzhe. Three-dimensional finite element analysis of stress and deformation characteristics of energy piles under inclined loads [J]. Chinese Journal of Geotechnical Engineering, 2021, 43(11): 2105 − 2111. (in Chinese)

    [25] GB 50007−2011, 建筑地基基础设计规范 [S]. 北京: 中国建筑工业出版社, 2012.

    GB 50007−2011, Code foe design of building for foundation [S]. Beijing: China Construction Industry Publishing Houses, 2012. (in Chinese)

  • 期刊类型引用(0)

    其他类型引用(3)

图(14)  /  表(1)
计量
  • 文章访问数:  210
  • HTML全文浏览量:  61
  • PDF下载量:  42
  • 被引次数: 3
出版历程
  • 收稿日期:  2022-01-27
  • 修回日期:  2022-05-23
  • 网络出版日期:  2022-06-24
  • 刊出日期:  2023-11-24

目录

/

返回文章
返回