ISSN 1000-3665 CN 11-2202/P
  • 中文核心期刊
  • GeoRef收录期刊
  • Scopus 收录期刊
  • 中国科技核心期刊
  • DOAJ 收录期刊
  • CSCD(核心库)来源期刊
  • 《WJCI 报告》收录期刊
欢迎扫码关注“i环境微平台”

基于能量法的轴横向荷载作用下单桩受力变形分析

张玲, 陈金海, 欧强

张玲, 陈金海, 欧强. 基于能量法的轴横向荷载作用下单桩受力变形分析[J]. 水文地质工程地质, 2020, 47(5): 81-91. DOI: 10.16030/j.cnki.issn.1000-3665.201912065
引用本文: 张玲, 陈金海, 欧强. 基于能量法的轴横向荷载作用下单桩受力变形分析[J]. 水文地质工程地质, 2020, 47(5): 81-91. DOI: 10.16030/j.cnki.issn.1000-3665.201912065
ZHANG Ling, CHEN Jinhai, OU Qiang. Deformation analysis of pile under combined axial and lateral loads by using the energy method[J]. Hydrogeology & Engineering Geology, 2020, 47(5): 81-91. DOI: 10.16030/j.cnki.issn.1000-3665.201912065
Citation: ZHANG Ling, CHEN Jinhai, OU Qiang. Deformation analysis of pile under combined axial and lateral loads by using the energy method[J]. Hydrogeology & Engineering Geology, 2020, 47(5): 81-91. DOI: 10.16030/j.cnki.issn.1000-3665.201912065

基于能量法的轴横向荷载作用下单桩受力变形分析

基金项目: 

国家自然科学基金项目 51678231

详细信息
    作者简介:

    张玲(1982-), 女, 博士, 副教授, 主要从事桩基础及软土地基处理方面的研究工作。E-mail: zhanglhd@163.com

  • 中图分类号: TU473.1+1

Deformation analysis of pile under combined axial and lateral loads by using the energy method

  • 摘要:

    轴向和横向荷载作用下单桩的受力变形分析是桩基研究的重点内容之一。单桩在水平荷载作用下会产生一定的水平位移与弯矩,而此时作用轴向荷载会使得桩体出现一定的压曲与附加弯矩,以致轴横向荷载作用下的单桩受力变形与单独作用水平荷载或轴向荷载的单桩存在较大的区别。故本文基于能量法,首先分别建立轴横向荷载作用下单桩的受力变形能量方程以及桩周土体能量方程,然后考虑桩土变形协调与一定的桩土相互作用,基于最小势能原理得到单桩变形控制微分方程,并采用幂级数法进行求解,最终得到轴横向荷载作用下单桩受力变形分析的幂级数解答。通过编程计算,将本文方法计算结果与试验结果、数值分析结果、规范法计算结果进行对比分析,验证了本文方法的合理性和可行性。在此基础上,基于本文解答进行了影响参数分析,结果表明:桩体长径比、桩土弹性模量比、桩周土模量深度变化系数均对轴横向受荷单桩的桩身水平位移与最大弯矩值有一定的影响,其中桩周土模量深度变化系数以不小于0.6为宜。

    Abstract:

    An analysis of deformation and inner force of pile under combined axial and lateral loads is one of the most important aspects for pile foundation. Under the action of lateral load, the single pile will produce certain lateral displacement and bending moment, at this time, the axial load will make the pile appear certain buckling and additional bending moment. The deformation of a single pile under the action of combined axial and lateral loads is different from that under the action of the lateral load or axial load. By using the energy method, this study firstly deduces the energy equations for the pile under combined axial and lateral loads and the surrounding soils, respectively. Combining with deformation compatibility and certain interaction for the pile-soil, the governing differential equation of pile deformation is obtained based on the principle of the minimum potential energy. A power series method is employed to solve the equation, and the power series solutions for the deformation and inner force of pile under the axial and lateral loads is finally obtained. To validate the presented solutions in this study, the results from current solutions are compared with the test results, numerical results, and the calculation results from the code method, showing a good agreement with the existed methods. The influence parameters are discussed. The results indicate that the ratio of length to pile diameter, the ratio of elastic modulus of pile to soil and the coefficient of modulus of soil around pile all have certain influence on the lateral displacement and the maximum bending moment of pile under the axial and lateral loads. It is found that the coefficient of modulus of soil around pile should be no less than 0.6.

  • 在桥梁、海洋等桩基的施工与服役过程中, 桩基除承受竖向荷载外, 往往还可能同时受到风荷载、汽车制动力、波浪冲刷、地震作用等水平荷载的作用。因此轴向、横向荷载作用下单桩的受力变形分析非常重要。目前高层建筑、桥梁及海洋工程的桩基设计中, 常规的做法是将轴向荷载与横向荷载分开考虑, 分别计算轴向与横向荷载作用下桩身的内力与位移, 而后基于小变形假定进行叠加, 该种叠加的方法虽计算相对简单, 但未考虑轴向荷载与横向荷载作用引起的耦合效应, 不能充分反映轴横向荷载作用下单桩的受力变形实际。

    已有的室内试验及数值分析结果均表明[1-5]: 轴向荷载对水平受荷桩的承载性状有较大影响, 如“P-△效应”等, 且轴横向荷载存在较为复杂的相互作用。日本学者横山幸满曾指出[6]: 应力叠加的原理并不适用桩体同时承受轴向与横向荷载的情况。

    为此, 轴横向荷载共同作用下单桩的受力变形应考虑轴横向应力的耦合作用, 目前相关的理论研究成果主要有: 赵明华等[7-8]基于“m法”, 求得了考虑轴横向荷载共同作用下的柔性桩解析解, 并对基桩的特性进行了深入的研究。李微哲等[9]通过假定地基系数满足(mz+C)的线性增长, 探讨分析了轴横向荷载偏心情况下的单桩内力与位移; 梁发云等[10]基于小变形假定对轴横向荷载共同作用下的单桩进行了研究, 并针对其主要影响因素进行了分析。张磊等[11-12]利用地基反力法求得了轴横向受荷桩的半解析解, 研究了自由段桩长及地基反力系数的分布情况等因素对单桩内力与变形的影响; Liang等[13]基于Mindlin方程提出了轴横向受荷桩的有限元解。张玲等[14]基于Winkler弹性地基梁理论, 利用相邻单元的连续条件与矩阵传递法研究了地基成层性对轴横向受荷桩的影响。周立朵等[15]基于p-y曲线法研究了轴横向荷载作用下的单桩桩土相互作用, 得到了荷载-位移曲线、弯矩分布等规律。Zhu等[16]通过矩阵传递法研究了轴向力及其沿桩身的分布对水平受荷响应的影响。

    但上述理论分析方法中, m值的选取与p-y曲线形式的确定不能直接从土工试验中获得, 而能量法的优点在于所需的土体参数较易获得, 且可在一定程度上考虑桩土复杂的相互作用。因此有学者尝试采用能量法研究荷载作用下单桩的受力变形[17-18]。而对于水平荷载作用下的单桩, Sun等[19]建立了土体在柱坐标下的位移方程, 并采用能量法去考虑水平受荷桩的受力与变形, Zhang等[20]、徐洞斌等[21]分别通过能量原理建立了水平荷载作用下桩土体的能量方程, 并通过有限差分求得了单桩在水平荷载作用下的受力响应。但上述研究仅针对水平受荷桩而言, 未考虑轴向荷载的影响。为此, 本文拟在前人的基础上, 基于最小势能原理, 并结合幂级数法得到轴横向荷载作用下单桩桩身水平位移及内力的解答, 以期为工程实际中轴横向荷载下桩基的计算分析提供参考。

    单桩在轴横向荷载作用时, 其桩身因水平荷载产生了较大的水平位移与弯矩, 而轴向力作用使得桩身会产生附加弯矩, 出现一定的压曲, 这个附加弯矩的存在又会使得桩身的水平位移增大, 而轴向力引起的桩身水平位移是非线性的, 故建立如图 1所示模型。

    图  1  单桩受力计算模型示意图
    Figure  1.  Schematic diagram of the calculation model for pile

    图 1所示的桩-土体系中: 桩长为L(m); 桩身半径为R(m); 桩身抗弯刚度为EPIP; 桩顶作用水平力H0, 轴向荷载P0, 弯矩M0。桩周土体弹性模量为Es, 泊松比为νs; 考虑桩周土体模量沿深度的变化, 其变化形式参考文献[20], 桩周土体在地面处的弹性模量为Es1, 在桩底处弹性模量为Es2, 桩身范围内距桩顶z深度处土体的弹性模量为:

    Es=Es2[αs+(1αs)z/L] (1)

    式中: αs--桩周土弹性模量沿深度的变化系数, αs= Es1/ Es2

    本文采用Sun等[19]提出的土体在柱坐标下的位移模式:

    {u(r,θ,z)=u(z)φ(r)cosθv(r,θ,z)=u(z)φ(r)sinθw=0 (2)

    式中: u, w, v--沿径向、轴向、转角方向的位移, 如图 2所示;

    图  2  坐标系与位移分量图
    Figure  2.  The coordinate system and displacement components

    r, θ, z--柱坐标径向、轴向与转角方向;

    u(z)--桩身随深度z变化的位移函数;

    φ(r)--r方向土体的无量纲位移函数。

    根据弹性力学几何关系与本构关系可得:

    [εrrεθθεzzγrθγθzγzr]=[u( dφ dr)cosθ00u( dφ dr)sinθ(du dz)φsinθ(du dz)φcosθ] (3)
    [σrrσθθσzzτrθτθzτzr]=A0[u dφ drcosθu dφ drsinθdu dzφsinθ du dzφcosθ] (4)

    式中: εrr, εθθ, εzz--r, θ, z方向的正应变分量;

    γ, γθz, γzr--, θz, zr方向的剪应变分量;

    σrr, σθθ, σzz--r, θ, z方向的正应力分量;

    τ, τθz, τzr--, θz, zr方向的剪应力分量。

    A0=[Es(1νs)(1+νs)(12νs)000Esνs(1+νs)(12νs)000Esνs(1+νs)(12νs)0000Es2(1+νs)0000Es2(1+νs)0000Es2(1+νs)]

    则桩土体系中各部分能量分别可表示为[22]:

    Πpile =12EpIp0L( d2u dz2)2 dz (5)
    Πsoil =1202π00σijεijr dr dθ dz (6)
    Πf=cτ20lz( du dz)2 dz (7)
    Πp=P020L( du dz)2 dz (8)
    ΠH=H0u|z=0 (9)
    ΠM=M0( du dz)|z=0 (10)

    式中: c--桩截面周长;

    τ--桩侧摩阻力, 根据文献[23]可采用静载试验获得或根据地区经验值取得;

    Πpile, Πsoil, Πf, Πp, ΠH, ΠM--桩身应变能、土体应变能、桩侧土体摩阻力势能、桩顶竖向荷载势能、桩顶水平荷载势能及桩顶弯矩势能。

    桩土体系总能量方程为:

    Π=Πpile +Πsoil +Πr+Πp+ΠH+ΠM=1202π00σijεijr dr dθ dz+12EpIp0L( d2u dz2)2 dz+cτ20Lz( du dz)2 dz+P020L( du dz)2 dzH0u|z=0M0( du dz)|z=0 (11)

    将式(3)、式(4)代入式(11), 可得桩土体系总能量方程为:

    Π=12EpIp0L( d2u dz2)2 dz+π2R0L[Es(34νs)2(1+νs)(12νs)u2( dφ dr)2+Es1+νs( du dz)2φ2]r dr dz+π2RL[Es(34νs)2(1+νs)(12νs)u2( dφ dr)2+Es1+νs( du dz)2φ2]r dr dz+cτ20Lz( du dz)2 dz+P020L( du dz)2 dzH0u|z=0M0( du dz)|z=0 (12)

    由最小势能原理可知, δП=0, 由式(12)求变分可得式(13):

    {{EpIpd4udz4πEs21+vsddz{[αs+(1αs)zL]dudz}Rφ2rdr+π(34vs)Es2[αs+(1αs)z/L]2(1+vs)(12vs)R(dφdr)2rdrucτ(dudz+zd2udz2)P0 d2u dz2}δudz+{EpIpd3udz3P0dudzcτzdudzπEs2αs1+vsR(dφdr)2rdrdudzH0}δu|z=0+{EpIpd3udz3+P0dudz+cτzdudz+πEs2[αs+(1αs)z/L]1+vsR(dφdr)2rdrdudz}δu|z=L+(EpIpd2udz2M0)δdudz|z=L+(EpIpd2udz2)δdudz|z=L}=0 (13)

    故由式(13)可推得桩身位移控制方程(14)及沿桩身任意深度处弯矩M及剪力H的表达式(17):

    EpIp d4u dz4πEs21+νs d dz{[αs+(1αs)zL]du dz}Rφ2r dr+π(34νv)Es2[αs+(1αs)z/L]2(1+νs)(12νs)R(dφ dr)2r drucτ( du dz+z d2u dz2)P0d2udz2=0(0zL) (14)

    其中令:

    n1=1R2Rφ2r dr,n2=R(dφ dr)2r dr (15)

    即桩身位移控制方程可化为:

    d4u dz4πn1R2Es2[αs+(1αs)z/L](1+νs)EpIp d2u dz2(1+νs)(cτz+P0)(1+νs)EpIp d2u dz2[πn1R2Es2(1αs)(1+νs)LEpIp+cτEpIp]du dz+π(34νs)n2Es2[αs+(1αs)z/L]2(1+νs)(12νs)EpIpu=0(0zL) (16)

    桩身任意深度处弯矩M及剪力H为:

    {M=EpIp d2u dz2H=EpIp d3u dz3πn1Es2αsR2+(1+νs)(P0+cτz)(1+νs)du dz (17)

    同理, 式(12)对φ取变分可得土体位移控制方程:

    r d2φ dr2+dφ dr(ψR)2rφ=0 (18)

    式中: ψ--无量纲系数, 决定了水平位移随径向r方向减小的速度。

    ψ可以表示为:

    ψ2=X1+X2X3+X4 (19)

    其中Χ1, Χ2, Χ3, Χ4分别为:

    X1=2n2R20L[αs+(1αs)zL](du dz)2 dzX2=n2R(34νs)n2(1+2n1)(12νs)u2(L)X3=(34νs)n212νs0L[αs+(1αs)zL](u)2 dzX4=12R(34νs)(1+2n1)n2(12νs)u2(L)

    式(18)满足边界条件: 当r=R时, φ(r)=1, 当r→∞时, φ(r→∞)=0, 则有:

    n1=12[K0(γ)]2{[K1(γ)]2[K0(γ)]2} (20)
    n2=12[K0(γ)]2{[K1(γ)+γK0(γ)]2(γ2+1)[K1(γ)]2} (21)

    式中: K0(·), K1(·)--修正后的第二类贝塞尔函数。

    式(14)可采用幂级数法求解, 具体求解过程参考文献[7], 故可得其解为:

    u(z)=j=03ajξj (22)

    其中: ξj=0nCj,nzn

    Cj,n+4=(A+Q)(n+1)(n+2)Cj,n+2(n+1)(n+2)(n+3)(n+4)+(B+T)(n+1)2Cj,n+4FCj,nGCj,n1(n+1)(n+2)(n+3)(n+4)(n1)Cj,4=2(A+Q)Cj,2+(B+C)Cj,1FCj,04×3×2×1Cj,k={0j=k;j,k=0,1,2,31j=k;j,k=0,1,2,3

    式中: a0, a1, a2, a3--待定系数。

    且:

    {A=πn1R2Es2αs(1+νs)EpIpB=πn1R2Es2(1αs)(1+νs)LEpIpT=cτEpIpQ=P0EpIpF=π(34νs)n2Es2αs2(1+νs)(12νs)EpIpG=π(34νs)n2Es2αs2(1+νs)(12νs)LEpIp (23)

    由式(22)可知桩身位移方程的幂级数函数均可由系数a0, a1, a2, a3表示, 而a0, a1, a2, a3可通过联立桩顶与桩底的边界条件进行求解。

    其中边界条件可由式(13)推得, 表达式为:

    {EpIpd3u dz3πn1E2αsR2+(1+νv)P0(1+νs)du dzH0=0(z=0)EpIpd2u dz2+M0=0(z=0)(桩顶自由)du dz=0(z=0)(桩顶嵌固) (24)
    {EpIp d2u dz2=0(z=L) (桩底自由) [EpIp d3u dz3πn1Es2αsR2+(1+νs)(P0+cτz)(1+νs)du dzπEs2R2(1+νs)(34νs)(1+2n1)n2(12νs)=0](z=L) (桩底自由) du dz=0(z=L) (桩底嵌固) u=0(z=L) (桩底嵌固)  (25)

    假定桩顶位移u0、转角φ0、弯矩M0以及剪力H0已知, 结合式(22)、(24)、(25)可得:

    {u|z=0=a0=u0 du dz|z=0=a1=φ0 d2u dz2|z=0=2a2=M0EpIp d3u dz3|z=0=6a3=H0EpIp+(A+Q)φ0 (26)

    由上述边界条件可求得a0, a1, a2, a3的值。再选定合适的项数, 即可求得水平位移的函数表达式, 其具体流程如图 3所示。当桩底嵌入岩层或者桩底土层较坚硬的时候, 认为桩端为固定边界。为方便求解, 下文将采用桩顶自由、桩底固定的边界条件进行求解验证。

    图  3  桩身水平位移求解程序图
    Figure  3.  Program diagram for solving horizontal displacement of pile

    Kerisel等[24]对水平受荷单桩进行了室内模型试验。模型桩桩长L=4.65 m, 桩径D=0.36 m, 桩顶自由, 桩底嵌固。受到的水平荷载H0=60 kN, 弯矩M0=69 kN·m, 桩身的弹性模量Ep=20 GPa, 土体的弹性模量Es=9.233 MPa。土体的泊松比νs=0.3, 桩周土模量沿深度的变化系数αs=1。Filho等[25]通过有限元与边界元软件验证了其试验结果, 图 4给出了本文方法计算结果与Kerisel等模型试验结果、Filho等有限元结果、《公路桥涵地基与基础设计规范》(以下简称规范法)推荐的简化方法计算结果的对比分析。

    图  4  水平荷载作用下的桩身水平位移
    Figure  4.  Lateral displacement of pile due to the horizontal load

    图 4可知, 本文方法计算所得桩身水平位移结果与上述几种方法所得结果都比较接近, 而且相比较于Filho等有限元方法, 本文计算方法所得的桩身水位位移更接近于Kerisel模型试验结果。就桩顶(z=0 m)水平位移而言, Kerisel模型试验结果为9.7 mm; Filho等有限元方法结果为10.1 mm, 相对比试验结果误差为4.12%;本文方法计算所得桩顶位移为10.08 mm, 相对比试验结果误差为3.9%。由此可见本文方法可用于分析水平受荷桩的受力变形分析。

    肖世伟等[26]利用有限元方法对马桥河桥钻孔桩基足尺试验所得结果进行了模拟分析, 并将实测数据进行反算, 给出了地基系数m=1.398×104 kN/m4与桩身弹性模量Eh=3.188×107 kN/m2, 试验的桩体参数分别为: 桩长L=17.0 m, 桩径D=0.85 m, 土层泊松比νs=0.3, 桩周土模量沿深度的变化系数αs=1, 则公式(1)即简化为Es=Es2, 土体弹性模量Es2=6 MPa, 本文方法计算所得桩顶位移与文献[26]结果如表 1所示。此外, 为比较本文计算桩身位移沿桩身的分布情况, 以文中参数建立弹性有限元模型, 取桩顶轴向荷载为313.60 kN的情况, 本文方法计算的桩身位移与有限元模拟的桩身位移对比结果如图 5所示。

    表  1  本文方法与文献[26]结果对比
    Table  1.  Comparison of the result of the method in this paper with that of the reference[26]
    下载: 导出CSV 
    | 显示表格
    图  5  桩身水平位移图
    Figure  5.  Lateral displacement of pile

    表 1可知, 本文方法计算所得结果与文献[26]所得结果基本吻合, 且随着轴向荷载的增大, 对桩顶水平位移的影响也随之增大, 且由图 5可知, 本文方法计算所得的桩身位移与弹性有限元分析所得结果拟合较好, 可应用于工程实际。

    Cho等[27]对水平受荷的嵌岩桩基进行了多组现场横向荷载试验, 该工程位于美国北卡罗来纳州的一个桥梁替换项目, 其中一组其桩侧有2层地基土, 第一层为深度L0=1 m的砂土层, 第二层为深度L1=3.3 m的泥板岩; 桩长为L=4.3 m, 桩体半径R=0.38 m, 桩体弹性模量Ep=24.8 GPa。

    Yang等[28]采用地基反力法与变分原理对该工程进行了验证分析, 并根据Cho等试验所得结果的p-y曲线进行反算, 假设桩侧路基反力模量Kh0在地面处为0, 沿深度线性增加, 求得首层土其增加系数为nhs, 0=19.5 MN/m3。泥板岩内的路基反力模量Kh1=180.6 MPa, nhs, 1=393.3 MN/m3, 路基反力模量的经验公式如式(27)所示, 其计算模型如图 6所示。

    { 岩体:K=EM(D/Dr)e2νs(EpIpEMD4)0.284土体:K=EsD(1νs2)Dref[EsD4EpIp]1/12Dr=0.305 m,Dref =1.0 m (27)
    图  6  Yang等计算模型图
    Figure  6.  Schematic diagram of the Yang calculation model

    式中: K--路基反力模量;

    EM--岩体模量。

    利用本文能量法对该算例进行验证分析, 桩体与土体参数如上所述, 本文所用的土体弹性模量与路基反力模量的换算公式为[29]:

    Kh=1.3Es1νs2(16EsR4EpIp)112 (28)

    式中: Kh--路基反力模量。

    桩体在泥板岩内基本不发生水平位移, 因此为简便计算, 本文利用砂土层的参数进行验证, 由式(28)可得砂土层底土体弹性模量Es2=17.73 MPa, 取桩周土模量沿深度的变化系数αs=0.1, 土体泊松比取为νs=0.3, 计算结果如图 7所示, 其中Yang所用p-y曲线法为利用Cho试验结果反算求得。

    图  7  水平荷载作用下的桩身水平位移
    Figure  7.  Lateral displacement of pile due to the horizontal load

    图 7可知, 本文方法所得的桩顶位移较文献[28]更接近实测值。

    取桩顶自由, 桩底嵌固, 桩的长径比为L/D=15, 取桩周土模量沿深度的变化系数αs=1, 土的泊松比取为νs=0.3, 桩土的弹性模量比值Ep/Es2分别为100, 500, 1 000, 2 000与5 000。桩顶水平荷载与轴向荷载分为取为H0=120 kN, P0=60 kN。其计算结果如图 8所示。

    图  8  桩土弹性模量比变化下的桩身水平位移和桩身弯矩
    Figure  8.  Horizontal displacement and bending moment of pile under the variation of elastic modulus ratio of pile to soil

    图 8可知, 桩身水平位移随着桩土模量比的增大而逐渐增大, 桩顶处水平位移在桩土模量比Ep/Es2=5 000处急剧增大, 在轴向压力的作用下, 桩身最大弯矩点不断向下移动。而当轴向荷载由0变为60 kN时, 其桩顶水平位移减小, 桩身最大弯矩点下移, 轴向荷载的存在使得桩身出现一定的压曲现象, 桩身的水平位移与弯矩均出现了增大的情况。

    取桩顶自由, 桩底嵌固, 桩土弹性模量比Ep/Es2=1 000, 取桩周土模量沿深度的变化系数αs=1, 土体的泊松比为νs=0.3, 桩体长径比L/D分别为5, 10, 15, 20, 25。桩顶水平荷载与轴向荷载分为取为H0=120 kN, P0=60 kN。为比较方便, 将桩长进行归一化可得图 9

    图  9  长径比变化下的桩身水平位移和桩身弯矩
    Figure  9.  Horizontal displacement and bending moment of pile under the variation of length-diameter ratio

    图 9可见, 当桩顶轴向荷载为60 kN时, 桩的工作性状会随桩长径比的变化而发生改变。当桩的长径比L/D大于10时, 桩身出现反弯。其原因是随桩的长径比增大, 桩体表现出柔性桩性质, 近地面处的土体发生屈服, 导致桩顶水平位移增大; 而当长径比L/D小于10时, 桩体受力近似于刚性桩, 桩身位移曲线近似为直线, 且桩身不出现反弯点。此外, 桩的长径比对桩身弯矩也有较大影响。随着桩的长径比L/D的增大, 桩身中上部的弯矩变化明显, 桩身最大弯矩减小, 且最大弯矩点逐渐上移。当长径比变化时, 轴向荷载对桩身位移与弯矩的影响与桩土模量比变化时的影响类似。

    取桩顶自由, 桩底嵌固, 桩的长径比L/D=15, 取桩周土模量沿深度的变化系数αs=1, 土的泊松比取为νs=0.3, 桩土弹性模量比值Ep/Es2=1 000。桩顶水平荷载取H0=120 kN, 轴向荷载P0分别为0, 600, 1 200, 2 400, 4 800 kN。计算结果如图 10所示。

    图  10  轴向荷载变化下的桩身水平位移和桩身弯矩
    Figure  10.  Horizontal displacement and bending moment of pile under the variation of axial load

    图 10可知, 轴向荷载对桩身水平位移的影响主要位于桩身上部2 m, 约为10倍桩径范围内。随着轴向荷载的增加, 桩身最大弯矩点不断下移, 其影响范围约位于地面以下4 m, 约为20倍桩径范围内。与文献[10]所得结论基本一致。而且, 随着P0/H0比值的增大, 轴向荷载对桩身水平位移及弯矩的影响明显增大, 其主要原因在于由轴向荷载引起的桩身挠曲变形是非线性的。

    桩周土体的弹性模量对桩身内力及变形有较大影响, 由式(1)可知, 系数αs决定了桩周土体弹性模量沿深度的变化趋势, 故对该系数做参数分析。

    取桩顶自由, 桩底嵌固, 桩的长径比为L/D=15, 土的泊松比取为νs=0.3, 桩土的弹性模量比值Ep/Es2=1 000, αs分别取为1, 0.8, 0.6, 0.4, 0.2。桩顶水平荷载与轴向荷载分为取为H0=120 kN, P0=60 kN, 计算结果如图 11所示。

    图  11  αs变化下的桩身水平位移和桩身弯矩
    Figure  11.  Horizontal displacement and bending moment of pile under the variation of αs

    图 11可知, 当桩顶轴向荷载由0增加到60 kN时, 桩顶水平位移从3.17 mm增大为3.86 mm, 增大比例为21.7%, 桩身最大弯矩由22.2 kN·m增大到31.36 kN·m, 增大比例为41.3%。而当桩顶轴向荷载不变时, 随桩周土模量沿深度的变化系数αs从1减小到0.2时, 桩顶水平位移的增大比例分别为19.17%、47.67%、94.3%、190.2%, 桩顶弯矩的增大比例分别为10.9%、25.7%、50.7%、95.7%。可知当桩周土模量沿深度的变化系数过小时, 即上部土层弹性模量过小, 会导致桩顶位移急剧增大, 桩身弯矩急剧增大, 最大弯矩点逐渐下移, 故在实际工程中可采取对浅部土层进行换填等方法去改善土质, 从而控制桩顶的水平位移与桩身弯矩。

    (1) 本文基于最小势能原理建立了轴横向荷载作用下单桩的能量方程, 一定程度上考虑了轴横向荷载的耦合作用, 将桩顶的轴横向荷载、桩侧摩阻力及桩土相互作用视为一个整体进行分析, 并给出了幂级数解答。

    (2) 参数分析结果表明: ①桩土弹性模量比对单桩内力与变形有较大影响。②桩身长径比小于10时, 桩体受力近似刚性桩容易引起桩身失稳。③轴向荷载对桩身水平位移与弯矩的影响主要集中与桩身上部20倍桩径范围内, 而对深处影响较小。④桩周土模量沿深度的变化系数αs小于0.6时, 单桩桩身水平位移与弯矩受其影响很大。

    (3) 本文方法可进一步推广到多层地基以及桩身变截面等情况。此外, 为求解方便, 本文仅以土体弹性模量沿深度线性变化的角度考虑土体的非线性, 未考虑土体的弹塑性变形, 这将在后续的工作中进一步深入研究。

  • 图  1   单桩受力计算模型示意图

    Figure  1.   Schematic diagram of the calculation model for pile

    图  2   坐标系与位移分量图

    Figure  2.   The coordinate system and displacement components

    图  3   桩身水平位移求解程序图

    Figure  3.   Program diagram for solving horizontal displacement of pile

    图  4   水平荷载作用下的桩身水平位移

    Figure  4.   Lateral displacement of pile due to the horizontal load

    图  5   桩身水平位移图

    Figure  5.   Lateral displacement of pile

    图  6   Yang等计算模型图

    Figure  6.   Schematic diagram of the Yang calculation model

    图  7   水平荷载作用下的桩身水平位移

    Figure  7.   Lateral displacement of pile due to the horizontal load

    图  8   桩土弹性模量比变化下的桩身水平位移和桩身弯矩

    Figure  8.   Horizontal displacement and bending moment of pile under the variation of elastic modulus ratio of pile to soil

    图  9   长径比变化下的桩身水平位移和桩身弯矩

    Figure  9.   Horizontal displacement and bending moment of pile under the variation of length-diameter ratio

    图  10   轴向荷载变化下的桩身水平位移和桩身弯矩

    Figure  10.   Horizontal displacement and bending moment of pile under the variation of axial load

    图  11   αs变化下的桩身水平位移和桩身弯矩

    Figure  11.   Horizontal displacement and bending moment of pile under the variation of αs

    表  1   本文方法与文献[26]结果对比

    Table  1   Comparison of the result of the method in this paper with that of the reference[26]

    下载: 导出CSV
  • [1] Meyerhof G G. Rupture surfaces in sand under oblique loads of the Soil Mechanics and Foundations Division[J]. ASCE, 1956(discussion). Journal, 82(SM3): 152-158.
    [2]

    MEYERHOF G G, HANNA A M. Ultimate bearing capacity of foundations on layered soils under inclined load[J]. Canadian Geotechnical Journal, 1978, 15(4): 565-572. DOI: 10.1139/t78-060

    [3] 赵春风, 王卫中, 赵程, 等. 砂土中竖向和弯矩荷载下单桩水平承载特性试验研究[J]. 岩石力学与工程学报, 2013, 32(1): 184-190. DOI: 10.3969/j.issn.1000-6915.2013.01.025

    ZHAO C F, WANG W Z, ZHAO C, et al. Lateral bearing capacity of single piles under vertical and moment load in sand[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(1): 184-190. (in Chinese) DOI: 10.3969/j.issn.1000-6915.2013.01.025

    [4] 赵春风, 刘丰铭, 邱志雄, 等. 砂土中竖向和水平荷载共同作用下的单桩承载特性研究[J]. 岩土工程学报, 2015, 37(1): 183-190.

    ZHAO C F, LIU F M, QIU Z X, et al. Study on bearing behavior of a single pile under combined vertical and lateral loads in sand[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(1): 183-190. (in Chinese)

    [5] 郑刚, 王丽. 成层土中倾斜荷载作用下桩承载力有限元分析[J]. 岩土力学, 2009, 30(3): 680-687. DOI: 10.3969/j.issn.1000-7598.2009.03.019

    ZHENG G, WANG L. Finite element analysis of bearing capacity of pile under inclined load in layered soil[J]. Rock and Soil Mechanics, 2009, 30(3): 680-687. (in Chinese) DOI: 10.3969/j.issn.1000-7598.2009.03.019

    [6] 横山幸满.桩结构物的计算方法和计算实例[M].唐业清,吴庆荪,译.北京:中国铁道出版社,1984:35—44.

    横山幸満. Calculation method and example of pile structure[M]. TANG Y Q, WU Q S, trans. Beijing: China Railway Publishing House, 1984: 35-44. (in Chinese)

    [7] 赵明华. 轴向和橫向荷载同时作用下的桩基计算[J]. 湖南大学学报(自然科学版), 1987, 14(2): 68-81.

    ZHAO M H. The calculation of piles under simultaneous axial and lateral loading[J]. Journal of Hunan University(Natural Science), 1987, 14(2): 68-81. (in Chinese)

    [8] 赵明华. 倾斜荷载下基桩的受力研究[D]. 长沙: 湖南大学, 2001.

    ZHAO M H. Study on the behavior of piles under inclined loads[D]. Changsha: Hunan University, 2001. (in Chinese)

    [9] 李微哲, 赵明华, 单远铭, 等. 倾斜偏心荷载下基桩内力位移分析[J]. 中南公路工程, 2005, 30(3): 53-57. DOI: 10.3969/j.issn.1674-0610.2005.03.015

    LI W Z, ZHAO M H, SHAN Y M, et al. Analysis of single pile under eccentric and inclined loading[J]. Journal of Central South Highway Engineering, 2005, 30(3): 53-57. (in Chinese) DOI: 10.3969/j.issn.1674-0610.2005.03.015

    [10] 梁发云, 陈海兵, 宋著. 轴向荷载对横向受荷桩承载性状影响的参数分析[J]. 土木建筑与环境工程, 2011, 33(6): 25-30.

    CHEN H B, LIANG F Y, SONG Z. Parameter analysis of effects of axial load on the lateral response of piles[J]. Journal of Civil, Architectural & Environment Engineering, 2011, 33(6): 25-30. (in Chinese)

    [11] 张磊, 龚晓南, 俞建霖. 考虑土体屈服的纵横荷载单桩变形内力分析[J]. 岩土力学, 2011, 32(8): 2441-2445. DOI: 10.3969/j.issn.1000-7598.2011.08.032

    ZHANG L, GONG X N, YU J L. Analysis of deformation and internal force for single pile under combined vertical and lateral loads considering yielding of soil[J]. Rock and Soil Mechanics, 2011, 32(8): 2441-2445. (in Chinese) DOI: 10.3969/j.issn.1000-7598.2011.08.032

    [12] 张磊, 龚晓南, 俞建霖. 纵横荷载下单桩地基反力法的半解析解[J]. 哈尔滨工业大学学报, 2011, 43(6): 96-100. DOI: 10.3969/j.issn.1009-1971.2011.06.017

    ZHANG L, GONG X N, YU J L. Semi-analytical solutions for single piles subjected to combined vertical and lateral loads by subgrade reaction method[J]. Journal of Harbin Institute of Technology, 2011, 43(6): 96-100. (in Chinese) DOI: 10.3969/j.issn.1009-1971.2011.06.017

    [13]

    LIANG F Y, CHEN H B, CHEN S L. Influences of axial load on the lateral response of single pile with integral equation method[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2012, 36(16): 1831-1845. http://www.researchgate.net/profile/Fayun_Liang/publication/263300513_Influences_of_axial_load_on_the_lateral_response_of_single_pile_with_integral_equation_method/links/561b32e108ae6d1730899326.pdf

    [14] 张玲, 赵明华, 赵衡. 倾斜荷载下桩柱式桥墩受力变形分析传递矩阵法[J]. 中国公路学报, 2015, 28(2): 69-76. DOI: 10.3969/j.issn.1001-7372.2015.02.009

    ZHANG L, ZHAO M H, ZHAO H. Transfer matrix method for deformation of pile type bridge pier under axial transverse load[J]. China Journal of Highway and Transport, 2015, 28(2): 69-76. (in Chinese) DOI: 10.3969/j.issn.1001-7372.2015.02.009

    [15] 周立朵, 孔纲强, 彭怀风, 等. 倾斜荷载下群桩承载特性理论分析[J]. 岩土力学, 2017, 38(9): 2647-2654.

    ZHOU L D, KONG G Q, PENG H F, et al. Theoretical analysis of bearing capacities of pile group under oblique load[J]. Rock and Soil Mechanics, 2017, 38(9): 2647-2654. (in Chinese)

    [16]

    ZHU M X, ZHANG Y B, GONG W M, et al. Generalized solutions for axially and laterally loaded piles in multilayered soil deposits with transfer matrix method[J]. International Journal of Geomechanics, 2017, 17(4): 04016104.

    [17] 张小敏, 祝启坤. 能量法在搅拌桩复合地基稳定性分析中的应用[J]. 水文地质工程地质, 2003, 30(4): 79-82.

    ZHANG X M, ZHU Q K. Application of energy analysis method in the stability study of deep mixing composite foundation[J]. Hydrogeology & Engineering Geology, 2003, 30(4): 79-82. (in Chinese)

    [18] 邹新军, 王英宇, 赵明华, 等. 基于能量法的非均质地基中单桩受扭弹塑性分析[J]. 水文地质工程地质, 2016, 43(2): 54-61. https://www.swdzgcdz.com/article/id/20160210

    ZOU X J, WANG Y Y, ZHAO M H, et al. Elastic-plastic analysis of a torsional single pile in non-homogeneous subsoil based on energy principle[J]. Hydrogeology and Engineering Geology, 2016, 43(2): 54-61. (in Chinese) https://www.swdzgcdz.com/article/id/20160210

    [19]

    SUN K M. Laterally loaded piles in elastic media[J]. Journal of Geotechnical Engineering, 1994, 120(8): 1324-1344. http://www.onacademic.com/detail/journal_1000035951786010_ac30.html

    [20]

    ZHANG L Y, ERNST H, EINSTEIN H H. Nonlinear analysis of laterally LoadedRock-socketed shafts[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2000, 126(11): 955-968. http://www.onacademic.com/detail/journal_1000037704851710_7b7e.html

    [21] 徐洞斌. H-T联合作用下基桩内力位移理论分析及试验研究[D]. 长沙: 湖南大学, 2015.

    XU D B. Theoretical analysis and experimental study of internal forces and deformation of single pile under H-T combined loading[D]. Changsha: Hunan University, 2015. (in Chinese)

    [22] 赵明华, 刘伟浩, 尹平保, 等. 基于能量法的陡坡段桥梁基桩屈曲稳定性分析[J]. 中南大学学报(自然科学版), 2016, 47(2): 586-592.

    ZHAO M H, LIU W H, YIN P B, et al. Buckling analysis of bridge piles in steep slopes based on energy method[J]. Journal of Central South University(Science and Technology), 2016, 47(2): 586-592. (in Chinese)

    [23] 建筑桩基技术规范: GJ94-2008[S]. 北京: 中国建筑工业出版社, 2008.

    Technical Code for Building Pile Foundations: GJ94-2008[S]. Beijing: China Architecture & Building Press, 2008. (in Chinese)

    [24]

    KE′RISEL J, ADAM M. Calcul des forces horizontales applicables aux fondations profondes dans les argiles el limons[J]. Annales d'Institut Technique du Batiment e des Travaux Publics, Paris 1967, 239: 1653-1694.

    [25]

    MATOS FILHO R, MENDONCA A V, DEPAIVA J B. Static boundary element analysis of piles submitted to horizontal and vertical loads[J]. Engineering Analysis With Boundary Elements, 2005, 29(3): 195-203. http://www.sciencedirect.com/science?_ob=ShoppingCartURL&_method=add&_eid=1-s2.0-S0955799704001614&originContentFamily=serial&_origin=article&_ts=1468008048&md5=271e8d3faeac079ac9282910998c9970

    [26] 肖世伟. 轴向力对横向受力桩的影响[J]. 西南交通大学学报, 1990, 25(4): 58-62.

    XIAO S W. Effect of axial force on laterally-loaded piles[J]. Journal of Southwest Jiaotong University, 1990, 25(4): 58-62. (in Chinese)

    [27]

    HWAN CHO K, CLARK S C, KEANEY B D, et al. Laterally loaded drilled shafts embedded in soft rock[J]. Transportation Research Record: Journal of the Transportation Research Board, 2001, 1772(1): 3-11. http://www.onacademic.com/detail/journal_1000039817278910_fdc2.html

    [28]

    YANG K, LIANG R. Numerical solution for laterally loaded piles in a two-layer soil profile[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2006, 132(11): 1436-1443.

    [29]

    Bowles J E. Foundation analysis and design[M]. 5th Ed. New York: McGraw Hill, 1996.

  • 期刊类型引用(6)

    1. 丁晓勇,许能权,邢皓枫. 砂土场地型钢桩现场载荷试验及其承载特性分析. 水文地质工程地质. 2024(03): 69-79 . 本站查看
    2. 尹平保,伍晗曦,贺炜,谭李俊英,刘江波,方敏. 推扭(H-T)荷载作用下刚性矩形基础承载特性研究. 水文地质工程地质. 2024(04): 125-134 . 本站查看
    3. 阳吉宝. 基于尖点突变模型确定静载试验的基桩极限承载力. 水文地质工程地质. 2024(05): 114-123 . 本站查看
    4. 艾启胜,王瑞杰,陈文德,任旭,彭亮,张超超. 基于改进深度神经网络的桩身内力测试数据自动定时采集算法. 土木工程与管理学报. 2024(05): 51-56 . 百度学术
    5. 吴琼. 基于动荷载作用的建筑地基施工二次变形预警技术研究. 建筑机械. 2022(02): 37-41 . 百度学术
    6. 张文扬. 地面堆载作用下邻近桩基侧向变形分析. 建筑机械. 2022(10): 60-64 . 百度学术

    其他类型引用(6)

图(11)  /  表(1)
计量
  • 文章访问数:  380
  • HTML全文浏览量:  179
  • PDF下载量:  703
  • 被引次数: 12
出版历程
  • 收稿日期:  2019-12-29
  • 修回日期:  2020-03-12
  • 网络出版日期:  2020-12-09
  • 刊出日期:  2020-08-31

目录

/

返回文章
返回