越顶破坏模式下沉埋桩受荷段及沉埋段推力算法

闫玉平

doi:10.16030/j.cnki.issn.1000-3665.202211007

越顶破坏模式下沉埋桩受荷段及沉埋段推力算法

闫玉平

中铁第一勘察设计院集团有限公司，陕西西安　710043

基金项目: 国家自然科学基金项目（41831290）

详细信息

作者简介: 闫玉平（1989-），男，博士，工程师，主要从事铁路设计、边坡稳定性分析、边坡加固等方面的研究工作。E-mail：2420907827@qq.com

中图分类号: P642.23

收稿日期: 2022-11-02

修回日期: 2022-12-03

刊出日期: 2023-05-15

Calculation method of thrust force of the embedded stabilizing piles under the overtop-sliding failure mode

YAN Yuping

China Railway First Survey and Design Institute Group Co. Ltd., Xi’an, Shaanxi　710043, China

Received Date: 02 November 2022

Revised Date: 03 December 2022

Publish Date: 15 May 2023

摘要

摘要:
滑坡推力的确定对于抗滑桩设计极其重要，沉埋桩作为对传统抗滑桩的优化，其受荷段推力的研究目前主要借助于模型试验和数值模拟，缺乏深入的理论分析。为了建立沉埋桩后侧受荷段及桩顶沉埋段滑坡推力计算方法，针对沉埋桩加固的基岩-覆盖层式滑坡，基于潜在越顶破坏模式，由桩顶位置将越顶滑面分为顶部、底部两段，其中，顶部滑面上水平方向合力即为沉埋段推力，可由积分求得，底部滑面上各个方向的力求解方法与此类似；在此基础上，利用刚体极限平衡理论对底部滑面与桩受荷段所围滑体进行受力分析，进而可得受荷段推力计算公式。实例分析表明：理论算法所得沉埋段与受荷段推力值与FLAC3D结果非常接近，其中，受荷段推力、沉埋段推力、设桩位置处总推力随沉埋比的增大而分别非线性减小、增大、减小；沉埋比位于0~0.67范围内时，沉埋段与受荷段推力之比由0缓慢增大到0.30~0.50，随着沉埋比增大到0.8，该比值急剧增大到1.47~2.12；一般沉埋深度下，沉埋段推力小于受荷段推力。沉埋桩推力的理论研究对于桩体内力优化、沉埋深度确定具有重要的现实意义，将有助于该桩型的进一步推广应用。
- 沉埋桩 /
- 越顶破坏 /
- 沉埋段推力 /
- 受荷段推力 /
- 极限平衡理论
Abstract:
Embedded piles act as an optimization structure compared with the traditional stabilizing pile. The determination of the thrust on the loading section is based mainly on the model test and numerical simulation, and there is a lack of in-depth theoretical analysis. For the bedrock-talus landslide reinforced by embedded piles, according to the potential overtop-sliding failure mode, the slide surface can be divided into top and bottom sections by the position of pile top, and the horizontal resultant force of the top section can be obtained by integration, which is the so-called thrust of the embedded section. Similarly, the force on the bottom section of the overtop-sliding surface can also be obtained. Based on the limit equilibrium theory, the force analysis of the sliding mass enclosed by the bottom sliding surface and the load section of the pile can be carried out, and the thrust on the loaded section can also be obtained. Example analyses show that the thrust of the embedded section and the loaded section obtained by the theoretical method are very consistent with the results of FLAC3D, the resultant force of the loading section decreases nonlinearly with the increase of the ratio, while the resultant force of the embedded section presents an opposite trend. With the increase of the embedded ratio from 0 to 0.67, the thrust ratio of embedded section and loading section increases slowly from 0 to 0.3−0.5. With the increase of embedded ratio from 0.67 to 0.8, the ratio increases sharply to 1.47−2.12. Generally, the thrust of embedded section is less than that of loaded section. The theoretical research of the thrust of the embedded pile is of great practical significance for the optimization of the pile internal force and the determination of the pile embedded depth, which will promote the further application of this structure.
- embedded stabilizing piles /
- overtop-sliding failure /
- thrust of the embedded section /
- thrust of the loaded section /
- the limit equilibrium theory

HTML全文

抗滑桩至今已有70多年的实践历史，其中运用最多的结构形式是单排桩^[1-3]。随着滑（边）坡治理规模的需要及设计理论的完善，抗滑桩逐渐朝着大直径^[4-5]、多排桩^[6-7]的方向发展，治理方案趋于综合化^[8]，结构型式趋于多样化^[9]，抗滑桩受力趋于更加合理^[10-12]。

抗滑桩加固坡体设计主要包括桩体受力分析与加固后坡体稳定性评价^[13-16]。传统抗滑桩，桩顶一般位于坡面，所以桩体需要承担后侧滑体传来的全部滑坡推力，这样桩体前侧滑体抗力无法有效发挥且桩身受力偏大^[6，17-18]。沉埋桩是在桩位不变的情况下，将桩顶埋入坡面一定深度，已有研究表明，该桩结构可以有效克服全长桩的不足^[17-18]。目前针对沉埋桩桩体受力的研究主要如下：熊治文^[19]通过模型试验与数值模拟，对沉埋桩受力分布规律、沉埋深度做了探讨，得出沉埋桩较全长桩受力更小的结论；雷文杰等^[20-21]利用模型试验与数值模拟分析了设桩位置不变时，不同桩顶埋深所对应设桩位置的滑坡总推力和桩身内力，发现沉埋深度增加可以有效降低后侧推力与桩身内力，且沉埋深度对设桩截面处总推力影响不大；雷文杰等^[22]采用有限元法对沉埋桩加固后的实际滑坡做了分析，得到沉埋段土体推力随其深度增加而增大，桩顶土体自身可以承受一定的推力，沉埋桩桩身受力较全长桩明显降低，桩身内力和桩后推力分布随沉埋深度增加变得更为合理等结果；宋雅坤等^[23]结合大型物理模型试验和数值方法分析了沉埋段和受荷段滑坡推力，初步得出加固后坡体稳定系数随沉埋段的增加而减小，受荷段与沉埋段推力随埋深变化趋势与雷文杰等^[22]的结论一致。但是针对沉埋桩受力的理论研究，目前鲜有文献报道。

在加固后坡体稳定性评价方面，Xiao^[17]在假定滑体为无黏性土的基础上，利用极限平衡理论和传递系数法得出了后排桩沉埋深度的计算方法，该方法可将桩顶埋深与坡体稳定系数建立关系，同样得出加固后坡体稳定系数随着沉埋深度的增加而降低，同时桩身内力也在减小的结论。Yan等^[18]基于沉埋桩加固滑坡可能发生的越顶破坏模式，结合传递系数法与极限平衡理论框架下的变分法，推导了对黏性土和无黏性土都适用的后排桩顶沉埋深度的解析解，得出与Xiao^[17]一致的结论。

综上，目前针对加固后坡体稳定性分析已经有了较为系统全面的认识，而桩体受力特别是受荷段推力的研究主要借助于模型试验和数值模拟，缺乏深入的理论研究。鉴于此，本文在Yan等^[18]所提越顶破坏模式的基础上，由桩顶位置将越顶滑面分为顶部、底部两段，通过积分求得两段滑面上合力，其中，顶部滑面上水平方向合力即为沉埋段推力；然后借助刚体极限平衡理论对底部段滑面与桩体受荷段所围滑体进行受力分析，求得受荷段推力。

1. 计算方法

桩体后缘传来的滑坡推力，一部分作用在桩上，不妨称此部分桩体为受荷段（图1）；另一部分越过桩顶滑体继续向前传递，通过前侧滑体自身稳定性提供抗力^[17]，桩顶距坡面深度部分不妨称之为沉埋段。

图 1. 沉埋桩加固滑坡横截面示意图

Figure 1. Sketch map of the cross section of a landslide reinforced with embedded pile

下载: 全尺寸图片幻灯片

Yan等^[18]提出的越顶破坏模式如图2所示，其中，极坐标系下越顶滑面BDC方程ρ(θ)及其正应力函数σ(θ)分别如式（1）（2）所示^[18]。

图 2. 沉埋桩越顶破坏模式

Figure 2. Overtop-sliding failure mode of the embedded piles

下载: 全尺寸图片幻灯片

$\rho \left( \theta \right) = A{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}$

(1)

$\begin{split} \sigma \left( \theta \right) &= \frac{{A\gamma }}{{1 + 9{{\left( {\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}} \right)}^2}}}{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}\left( {\sin \theta - 3\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\cos \theta } \right) \\ &+B{{\rm{e}}^{2\tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }} - \dfrac{c}{{\tan \varphi }} \end{split}$

(2)

式中：θ——极角/（°）；

F_s——与越顶破坏模式对应的稳定系数；

$\gamma$ ——滑体重度/（kN·m⁻³）；

c、φ——黏聚力/kPa、内摩擦角/（°）；

A、B——积分常数。

越顶滑面段BD、DC，受荷段BDH、HDE块的受力分析如图3所示。做如下假定：

图 3. 各块体受力分析模型

Figure 3. Mechanical analysis model of each typical rigid block

下载: 全尺寸图片幻灯片

（1）越顶滑面与桩体受荷段之间的初始滑面BHE处于极限状态且满足库仑强度准则。

（2）初始滑面BHE上只有一个折点H，无折点情况为其特解，多折点求解方法类似一个折点。

（3）BDH和HDE块都为刚体，除越顶滑面段BD外其他边界上力作用点都位于相应作用面中点。

下面对图3所示的越顶滑面BDC、刚性块体BDH和HDE分别进行受力分析。

1）越顶滑面BD段上合力求解

图3（a）中越顶滑面BDC上各点处的正应力函数如式（2）所示，根据摩尔-库仑强度理论，切向应力函数τ(θ)可表示为式（3）。滑面BDC上长度增量记为dl，任意一点切向力和法向力增量分别记为dT和dN，则dl、dT和dN可分别表示为式（4）、式（5）、式（6）。

$\begin{split} &\tau \left( \theta \right) = \sigma \left( \theta \right)\tan \varphi + c = \\ &\frac{{A\gamma \dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}}}{{1 + 9{{\left( {\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}} \right)}^2}}}{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}\left( {\sin \theta - 3\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\cos \theta } \right) + B\frac{{\tan \varphi }}{{{F_{\rm{s}}}}}{{\rm{e}}^{2\tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }} \end{split}$

(3)

${\rm{d}}l = \frac{{\rho \left( \theta \right)}}{{\cos \varphi }}{\rm{d}}\theta = \frac{{A{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}}}{{\cos \varphi }}{\rm{d}}\theta$

(4)

${\rm{d}}T = \tau \left( \theta \right){\rm{d}}l = \left[ {\frac{{A\gamma \dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}}}{{1 + 9{{\left( {\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}} \right)}^2}}}{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}\left( {\sin \theta - 3\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\cos \theta } \right) + B\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}{{\rm{e}}^{2\tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}} \right]\dfrac{{A{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}}}{{\cos \varphi }}{\rm{d}}\theta$

(5)

${\rm{d}}N = \sigma \left( \theta \right){\rm{d}}l = \left[ {\frac{{A\gamma }}{{1 + 9{{\left( {\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}} \right)}^2}}}{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}\left( {\sin \theta - 3\dfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\cos \theta } \right) + B{{\rm{e}}^{2\tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }} - \dfrac{c}{{\tan \varphi }}} \right]\dfrac{{A{{\rm{e}}^{ - \tfrac{{\tan \varphi }}{{{F_{\rm{s}}}}}\theta }}}}{{\cos \varphi }}{\rm{d}}\theta$

(6)

令滑面BDC上x和y方向上力的增量分别为dF_x和dF_y，则有：

$\left\{ \begin{gathered} {\rm{d}}{F_x} = {\rm{d}}N\cos \left( {\theta - \varphi } \right) + {\rm{d}}T\sin \left( {\theta - \varphi } \right) \\ {\rm{d}}{F_y} = {\rm{d}}N\sin \left( {\theta - \varphi } \right) - {\rm{d}}T\cos \left( {\theta - \varphi } \right) \\ \end{gathered} \right.$

(7)

以O’点为矩心，可求得滑面BD段上的水平向合力F_xBD、竖向合力F_yBD及对O’点的力矩M_BD，即：

$\left\{ \begin{gathered} {F_{xBD}} = \int\limits_{{\theta _B}}^{{\theta _D}} {{\rm{d}}{F_x}} = \int\limits_{{\theta _B}}^{{\theta _D}} {\frac{{\left[ {\sigma \left( \theta \right)\cos \left( {\theta - \varphi } \right) + \tau \left( \theta \right)\sin \left( {\theta - \varphi } \right)} \right]\rho \left( \theta \right)}}{{\cos \varphi }}{\rm{d}}\theta } \\ {F_{yBD}} = \int\limits_{{\theta _B}}^{{\theta _D}} {{\rm{d}}{F_y}} = \int\limits_{{\theta _B}}^{{\theta _D}} {\frac{{\left[ {\sigma \left( \theta \right)\sin \left( {\theta - \varphi } \right) - \tau \left( \theta \right)\cos \left( {\theta - \varphi } \right)} \right]\rho \left( \theta \right)}}{{\cos \varphi }}{\rm{d}}\theta } \\ {M_{BD}} = c{A^2}\int\limits_{{\theta _B}}^{{\theta _D}} {{{\rm{e}}^{ - 2\tfrac{{\tan \varphi }}{{{F_s}}}\theta }}{\rm{d}}\theta } \\ \end{gathered} \right.$

(8)

2）块体BDH受力分析

BDH块受力分析模型如图3（b）所示，B、D、H各点坐标分别记为（x_B，y_B）、（x_D，y_D）、（x_H，y_H），同时令BH中点坐标为（x₁，y₁），HD中点坐标为（x₂，y₂），BH和HD长度记为l_BH和b。由静力平衡条件（以O’为矩心）及假定（1）可得：

$\left\{ \begin{gathered} \sum {{F_x} = 0} \\ \sum {{F_y} = 0} \\ \sum {{M_{O'}}} = 0 \\ {T_1} = {N_1}\frac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}} + \frac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}{l_{BH}} \\ \end{gathered} \right.$

(9)

式中：c_m、ϕ_m——原滑面BHE上的黏聚力、内摩擦角；

T₁、N₁——BH面上的切向力、法向力/kN。

为简化计算，曲线BD按直线边计算，则其重力（G₁）可表示为如下形式：

${G_1} = \frac{1}{2}\gamma {l_{BH}}b\sin \left( {{\text{π}} - \eta - \omega } \right)$

(10)

令BDH重心到O’点的水平距离为d₁，则可求得如下BH线方程以及经过BH中点的垂线方程：

$\left\{ \begin{gathered} \frac{{{y_H} - {y_B}}}{{{x_H} - {x_B}}}x - y + \frac{{{x_H}{y_B} - {x_B}{y_H}}}{{{x_H} - {x_B}}} = 0 \\ \frac{{{x_H} - {x_B}}}{{{y_H} - {y_B}}}x + y - {y_1} - \frac{{{x_H} - {x_B}}}{{{y_H} - {y_B}}}{x_1} = 0 \\ \end{gathered} \right.$

(11)

O’到线段BH及其过中点垂线的距离分别记为l₁和l₂，其表达式如式（12）所示：

$\left\{ \begin{gathered} {l_1} = \frac{{\left| {\dfrac{{{y_H} - {y_B}}}{{{x_H} - {x_B}}}{x_0} - {y_0} + \dfrac{{{x_H}{y_B} - {x_B}{y_H}}}{{{x_H} - {x_B}}}} \right|}}{{\sqrt {{{\left( {\dfrac{{{y_H} - {y_B}}}{{{x_H} - {x_B}}}} \right)}^2} + 1} }} \\ {l_2} = \frac{{\left| {\dfrac{{{x_H} - {x_B}}}{{{y_H} - {y_B}}}\left( {{x_0} - {x_1}} \right) + {y_0} - {y_1}} \right|}}{{\sqrt {{{\left( {\dfrac{{{x_H} - {x_B}}}{{{y_H} - {y_B}}}} \right)}^2} + 1} }} \\ \end{gathered} \right.$

(12)

同理，可得DH线方程及经过DH中点的垂线方程（为了节省篇幅，在此不再给出），记O’到线段DH及其过中点垂线的距离分别记为l₃和l₄，将式（10）—（12）代入式（9），则有：

$\left\{ \begin{gathered} {F_{xBD}} + {N_1}\sin \eta - {T_1}\cos \eta - {T_2}\cos \omega - {N_2}\sin \omega = 0 \\ {F_{yBD}} + {G_1} - {N_1}\cos \eta - {T_1}\sin \eta + {T_2}\sin \omega - {N_2}\cos \omega = 0 \\ {M_{BD}} - {G_1}{d_1} + {N_1}{l_2} - {T_1}{l_1} - {T_2}{l_2} - {N_2}{l_4} = 0 \\ {T_1} = {N_1}\frac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}} + \frac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}{l_{BH}} \\ \end{gathered} \right.$

(13)

式中：η、ω——BH边、DH边与水平面的夹角/（°）；

T₂、N₂——DH面上的切向力、法向力/kN。

求解方程组（13），可解得各个边界上的力如式（14）所示：

$\left\{ \begin{gathered} {N_1} ={ \dfrac{{\left[ {{F_{xBD}}\sin \omega + \left( {{F_{yBD}} + {G_1}} \right)\cos \omega } \right]\left( {{l_3}\cos \omega + {l_4}\sin \omega } \right) + \left[ {{l_3}\sin \eta + {l_1}\sin \omega - \left( {{l_3}\cos \omega + {l_4}\sin \omega } \right)\sin \left( {\omega + \eta } \right)} \right]\dfrac{{{c_m}}}{{{F_{\rm{s}}}}}{l_{BH}} - \left( {{M_{BD}} - {G_1}{d_1}} \right)\sin \omega + \left( {{F_{yBD}} + {G_1}} \right){l_3}}}{{{l_2}\sin \omega - {l_3}\cos \eta + \left( {{l_3}\cos \omega + {l_4}\sin \omega } \right)\left[ {\cos \left( {\omega + \eta } \right) + \sin \left( {\omega + \eta } \right)\dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]}}} \\ {T_1} = {N_1}\dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}} + \dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}l{}_{BH} \\ {N_2} = {F_{xBD}}\sin \omega + \left( {{F_{yBD}} + {G_1}} \right)\cos \omega - \dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}l{}_{BH}\sin \left( {\omega + \eta } \right) - \left[ {\cos \left( {\omega + \eta } \right) + \sin \left( {\omega + \eta } \right)\dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]{N_1} \\ {T_2} = {F_{xBD}}\cos \omega - \left( {{F_{yBD}} + {G_1}} \right)\sin \omega - \dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}l{}_{BH}\cos \left( {\omega + \eta } \right) + \left[ {\sin \left( {\omega + \eta } \right) - \cos \left( {\omega + \eta } \right)\dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]{N_1} \\ \end{gathered} \right.$

(14)

3）块体HDE受力分析

HDE块受力分析模型如图3（c）所示，记HE和DE长度分别为j和k。由静力平衡条件（以DE中点为矩心）及假定（1）可得式（15）。求解式（15），则可解得边界HE和DE上的力如式（16）（17）所示。

$\left\{ \begin{gathered} {N_2}\cos \xi + {T_2}\sin \xi + {N_3}\sin \delta - {T_3}\cos \delta - {N_p} = 0 \\ {N_2}\sin \xi - {T_2}\cos \xi - {N_3}\cos \delta - {T_3}\sin \delta - {T_p} + {G_2} = 0 \\ \begin{split} &\left( {{N_3}\sin \delta - {T_3}\cos \delta } \right)\frac{b}{2}\cos \xi - \left( {{N_3}\cos \delta + {T_3}\sin \delta } \right)\frac{b}{2}\sin \xi + \left( {{N_2}\sin \xi - {T_2}\cos \xi } \right)\frac{j}{2}\cos \delta - \\ &\left( {{N_2}\cos \xi + {T_2}\sin \xi } \right)\frac{j}{2}\sin \delta + \frac{j}{3}{G_2}\sin \left( {\frac{\text{π} }{2} - \delta } \right) = 0 \end{split} \\ {T_3} = {N_3}\frac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}} + \frac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}j \\ \end{gathered} \right.$

(15)

$\left\{ \begin{gathered} {N_3} = \dfrac{{s\dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}j + {N_2}j\tan \left( {\delta - \xi } \right) + {T_2}j - \dfrac{{2{G_2}j\sin \left( {\dfrac{\text{π} }{2} - \delta } \right)}}{{3\cos \left( {\delta - \xi } \right)}}}}{{b\left[ {\tan \left( {\delta - \xi } \right) - \dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]}} \\ {T_3} = \dfrac{{\left[ {{N_2}j\tan \left( {\delta - \xi } \right) + {T_2}j - \dfrac{{2{G_2}j\sin \left( {\dfrac{\text{π} }{2} - \delta } \right)}}{{3\cos \left( {\delta - \xi } \right)}}} \right]\dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}} + b\dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}j\tan \left( {\delta - \omega } \right)}}{{b\left[ {\tan \left( {\delta - \omega } \right) - \dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]}} \\ \end{gathered} \right.$

(16)

$\left\{ \begin{gathered} \begin{split} {N_{\rm{p}}} &= {N_2}\cos \xi + {T_2}\sin \xi + \dfrac{{b\dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}} + {N_2}\tan \left( {\delta - \xi } \right) + {T_2} - \dfrac{{2{G_2}\sin \left( {\dfrac{\text{π} }{2} - \delta } \right)}}{{3\cos \left( {\delta - \xi } \right)}}}}{{b\left[ {\tan \left( {\delta - \xi } \right) - \dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]}}j\sin \delta \\ &-\dfrac{{\left[ {{N_2}\tan \left( {\delta - \xi } \right) + {T_2} - \dfrac{{2{G_2}\sin \left( {\dfrac{\text{π} }{2} - \delta } \right)}}{{3\cos \left( {\delta - \xi } \right)}}} \right]\dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}} + b\dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}\tan \left( {\delta - \xi } \right)}}{{b\left[ {\tan \left( {\delta - \xi } \right) - \dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]}}j\cos \delta \end{split} \\ \begin{split} {T_{\rm{p}}} &= {N_2}\sin \xi - {T_2}\cos \xi - \dfrac{{b\dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}} + {N_2}\tan \left( {\delta - \xi } \right) + {T_2} - \dfrac{{2{G_2}\sin \left( {\dfrac{\text{π} }{2} - \delta } \right)}}{{3\cos \left( {\delta - \xi } \right)}}}}{{b\left[ {\tan \left( {\delta - \xi } \right) - \dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]}}j\cos \delta \\ &- \dfrac{{\left[ {{N_2}\tan \left( {\delta - \xi } \right) + {T_2} - \dfrac{{2{G_2}\sin \left( {\dfrac{\text{π} }{2} - \delta } \right)}}{{3\cos \left( {\delta - \xi } \right)}}} \right]\dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}} + b\dfrac{{{c_{\rm{m}}}}}{{{F_{\rm{s}}}}}\tan \left( {\delta - \xi } \right)}}{{b\left[ {\tan \left( {\delta - \xi } \right) - \dfrac{{\tan {\varphi _{\rm{m}}}}}{{{F_{\rm{s}}}}}} \right]}}j\sin \delta + {G_2} \end{split} \\ \end{gathered} \right.$

(17)

式中：T₃、N₃——HE面上切向力、法向力/kN；

T_p、N_p——DE面上切向力、受荷段推力/kN；

δ——HE边与水平面的夹角/（°）；

ξ——HD边与竖向的夹角/（°）；

G₂——块体HDE的重力/kN。

4）沉埋段推力求解

图3（a）中越顶滑面DC段水平方向合力即为沉埋段推力，记为F_xDC，求解方法同滑面BD段，表达式如下：

$F_{x D C}=\int_{\theta_{D}}^{\theta_C} \mathrm{d} F_{x}=\int_{\theta_{D}}^{\theta_{C}} \frac{[\sigma(\theta) \cos (\theta-\varphi)+\tau(\theta) \sin (\theta-\varphi)] \rho(\theta)}{\cos \varphi} \mathrm{d} \theta$

(18)

上述求解过程给出了越顶破坏模式下，沉埋段与受荷段推力的严格解答，该求解过程可以借助MATLAB完成^[24]。

2. 实例分析

四川境内宝成铁路沿线某中风化大理岩上覆碎石土的基覆式滑坡主横断面如图4所示，通过现场地质勘查与室内试验，坡体主要物理力学参数如表1所示。利用传递系数法^[25-26]求得自然坡体稳定系数为1.06。设计安全系数为1.20，拟采用沉埋桩对该坡体加固，桩体截面尺寸为2 m×3 m，平面外桩间距为5 m，设桩位置见图4，该位置滑体厚度为15.1 m。

图 4. 沉埋桩加固宝成铁路沿线某滑坡横断面图

Figure 4. Cross section of a practical landslide reinforced with embedded stabilizing piles

下载: 全尺寸图片幻灯片

表 1. 抗滑桩及坡体主要物理力学参数

Table 1. Main physical and mechanical parameters of the landslide and piles

材料类型	土体重度 /（kN·m^–3）	黏聚力 /kPa	内摩擦角 /（°）	弹性模量 /MPa	泊松比
碎石土	22	11	20.6	50	0.33
大理岩	23	450	37.0	600	0.25
抗滑桩	25			30 000	0.22
注：抗滑桩视为弹性体，故不考虑其黏聚力和内摩擦角。

| Show Table

DownLoad: CSV

为了便于表述，引入Yan等^[18]对沉埋比（η）的定义，即：η=桩顶埋深/设桩点滑体厚度。通过Yan等^[18]的方法可得，设计安全系数为1.20时，η=0.364，即桩顶最大沉埋深度为5.5 m。利用本文方法可得此沉埋深度下沉埋段推力及受荷段推力分别为439.9 kN/m和2949.6 kN/m，其比值为0.149。

为了进一步验证本文方法的合理性，利用FLAC3D进行数值模拟。数值模型含47436个8节点六面体单元（图5），坡体采用服从摩尔-库仑屈服准则和关联流动法则的理想弹塑性本构模型，桩体视为弹性材料，利用结构单元模拟。模型前后左右4个边界采用水平位移约束，底面采用水平和竖向位移约束。采用强度折减法进行数值模拟，得到沉埋深度为5.5 m时，该滑坡的稳定系数为1.21，此时临界滑动面如图6所示，该结果与Yan等^[18]的计算结果非常相近，在此不再赘述。

图 5. 沉埋桩加固滑坡数值模型

Figure 5. Numerical model of the practical landslide reinforced with embedded pile

下载: 全尺寸图片幻灯片

图 6. 数值模拟得到的临界滑动面

Figure 6. Critical slip surface obtained with the shear strength reduction method

下载: 全尺寸图片幻灯片

图7所示为数值模拟得到的越顶破坏模式下滑坡的水平应力分布云图，由此易得沉埋段与受荷段推力分别为428.9 kN/m和2 788.9 kN/m，其比值为0.154。可见，数值结果与理论计算结果非常吻合，这说明了本文方法的合理性。

图 7. 沉埋桩加固滑坡的水平应力云图

Figure 7. Contour of horizontal stress of the reinforced landslide using the embedded pile

下载: 全尺寸图片幻灯片

3. 讨论

下面对受荷段推力、沉埋段推力、设桩位置处总推力与沉埋比的关系逐一讨论。不同沉埋比对应不同的越顶滑面与稳定系数，为了论述的完整性，首先利用文献[18]与FLAC3D分别获得加固后坡体稳定系数（F_s）随沉埋深度的变化情况，如图8所示，可见沉埋比（η）越大，坡体稳定系数越小。

图 8. 加固后坡体稳定系数随η的变化情况

Figure 8. Variation of stability coefficient of reinforced landslide with the embedded ratio

下载: 全尺寸图片幻灯片

图9为理论方法与FLAC3D得出的受荷段推力及沉埋段推力随η的变化情况。可见，FLAC3D结果与理论计算结果非常接近，随着η由0增大到1，沉埋段推力由0非线性增大到2000 kN/m左右，受荷段推力则由4300 kN/m非线性减小到0，同时，η位于0.7～0.8区间时，受荷段与沉埋段推力分别急剧减小与增加。

图 9. 受荷段及沉埋段推力随η的变化情况

Figure 9. Variation of the thrust behind the pile and the thrust in the submerged section with the embedded ratio

下载: 全尺寸图片幻灯片

将沉埋段推力与受荷段推力叠加可得不同沉埋比时设桩位置处总推力的变化曲线（图10），可见2种方法求得的设桩位置处总推力非常接近，且变化趋势一致。η为0时，桩位处总滑坡推力明显大于η为1时的情况，前者为4350 kN/m，后者为2040 kN/m，两者比值为2.13，这可由图8解释，全长桩时，加固后坡体稳定系数最大，而无桩情况下，坡体稳定系数最小。

图 10. 设桩处总滑坡推力随η的变化情况

Figure 10. Variation of the total thrust at the section of pile position with the embedded ratio

下载: 全尺寸图片幻灯片

进一步将2种方法得到的沉埋段与受荷段推力之比随η的变化情况绘于图11，可见2种方法求得的比值随η非线性增大，且变化趋势非常一致。当η位于0～0.67时，沉埋段推力与桩后受荷段推力之比由0缓慢增大到0.30～0.50，当η在0.67～0.80范围内时，该比值急剧增大到1.47～2.12，此时抗滑桩作用效果已经非常小，说明了该沉埋深度在工程上需要避免。此外，图11中也标出了沉埋段推力与桩后受荷段推力之比为1时，η为0.72～0.74。综合分析可得，沉埋深度合理时，沉埋段推力小于受荷段推力。

图 11. 沉埋段、受荷段推力之比与η的关系曲线

Figure 11. Variation curve of thrust ratio of the submerged section and the loading section with the embedded ratio

下载: 全尺寸图片幻灯片

4. 结论

（1）针对沉埋桩加固基覆式滑坡，基于潜在越顶破坏模式，由桩顶位置将越顶滑面分为顶部、底部两段。对顶部滑面上的力进行积分可得沉埋段推力，利用刚体极限平衡理论对底部滑面与受荷段所围滑体进行受力分析可得受荷段推力，进而可得受荷段与沉埋段推力随桩顶埋深的变化规律、不同沉埋深度下设桩位置处滑坡推力、受荷段与沉埋段推力比值等。

（2）理论计算所得沉埋段与受荷段推力值与FLAC3D结果非常接近，其中，受荷段推力、沉埋段推力、设桩位置处总推力随沉埋比的增大而分别非线性减小、增大、减小。

（3）沉埋比为0和1时，桩位处总设计滑坡推力分别为4350 kN/m与2040 kN/m，前者与后者比值为2.13；沉埋比为0.7～0.8时，受荷段与沉埋段推力分别急剧减小与增加。沉埋比为0～0.67时，沉埋段与受荷段推力之比由0缓慢增大到约0.3～0.5，沉埋比由0.67增大到0.80时，该比值急剧增大到1.47～2.12。

（4）为了充分发挥沉埋桩加固效果，一般沉埋深度下，沉埋段推力应小于受荷段推力。

图 1 沉埋桩加固滑坡横截面示意图

Figure 1.

下载: 全尺寸图片幻灯片

图 2 沉埋桩越顶破坏模式

Figure 2.

下载: 全尺寸图片幻灯片

图 3 各块体受力分析模型

Figure 3.

下载: 全尺寸图片幻灯片

图 4 沉埋桩加固宝成铁路沿线某滑坡横断面图

Figure 4.

下载: 全尺寸图片幻灯片

图 5 沉埋桩加固滑坡数值模型

Figure 5.

下载: 全尺寸图片幻灯片

图 6 数值模拟得到的临界滑动面

Figure 6.

下载: 全尺寸图片幻灯片

图 7 沉埋桩加固滑坡的水平应力云图

Figure 7.

下载: 全尺寸图片幻灯片

图 8 加固后坡体稳定系数随η的变化情况

Figure 8.

下载: 全尺寸图片幻灯片

图 9 受荷段及沉埋段推力随η的变化情况

Figure 9.

下载: 全尺寸图片幻灯片

图 10 设桩处总滑坡推力随η的变化情况

Figure 10.

下载: 全尺寸图片幻灯片

图 11 沉埋段、受荷段推力之比与η的关系曲线

Figure 11.

下载: 全尺寸图片幻灯片

表 1 抗滑桩及坡体主要物理力学参数

Table 1. Main physical and mechanical parameters of the landslide and piles

材料类型	土体重度 /（kN·m^–3）	黏聚力 /kPa	内摩擦角 /（°）	弹性模量 /MPa	泊松比
碎石土	22	11	20.6	50	0.33
大理岩	23	450	37.0	600	0.25
抗滑桩	25			30 000	0.22
注：抗滑桩视为弹性体，故不考虑其黏聚力和内摩擦角。

下载: 导出CSV

参考文献(26)

[1]	ITO T,MATSUI T,HONG W P. Design method for stabilizing piles against landslide: One row of piles[J]. Soils and Foundations,1981,21(1):21 − 37. doi: 10.3208/sandf1972.21.21
[2]	RANDOLPH M F. The response of flexible piles to lateral loading[J]. Géotechnique,1981,31(2):247 − 259.
[3]	常晓军,葛伟亚,于洋,等. 福建省永泰县东门旗山滑坡诱发机理与防治措施[J]. 吉林大学学报(地球科学版),2019,49(4):1063 − 1072. [CHANG Xiaojun,GE Weiya,YU Yang,et al. Mechanism and mitigation measures of Qishan landslide of Yongtai in Fujian Province[J]. Journal of Jilin University (Earth Science Edition),2019,49(4):1063 − 1072. (in Chinese with English abstract)
[4]	CHOW Y K. Analysis of piles used for slope stabilization[J]. International Journal for Numerical and Analytical Methods in Geomechanics,1996,20(9):635 − 646. doi: 10.1002/(SICI)1096-9853(199609)20:9<635::AID-NAG839>3.0.CO;2-X
[5]	HASSIOTIS S,CHAMEAU J L,GUNARATNE M. Design method for stabilization of slopes with piles[J]. Journal of Geotechnical and Geoenvironmental Engineering,1997,123(4):314 − 323. doi: 10.1061/(ASCE)1090-0241(1997)123:4(314)
[6]	郑颖人, 雷文杰, 赵尚毅, 等. 抗滑桩设计中的两个问题[J]. 公路交通科技, 2005, 22（增刊1）: 45 − 51 ZHENG Yingren, LEI Wenjie, ZHAO Shangyi, et al. Two problems in design methods of anti-slide piles[J]. Journal of Highway and Transportation Research and Development, 2005, 22（Sup 1）: 45 − 51. （in Chinese with English abstract）
[7]	闫玉平,肖世国. 双排抗滑桩加固滑坡的前桩后侧推力算法[J]. 水文地质工程地质,2021,48(4):55 − 63. [YAN Yuping,XIAO Shiguo. A calculation method for thrust on the fore piles of double-row stabilizing piles used to reinforce landslides[J]. Hydrogeology & Engineering Geology,2021,48(4):55 − 63. (in Chinese with English abstract)
[8]	王文沛,殷跃平,王立朝,等. 排水抗滑桩技术研究现状及展望[J]. 水文地质工程地质,2022,49(2):73 − 83. [WANG Wenpei,YIN Yueping,WANG Lichao,et al. Status and prospects of studies on anti-slide shaft technology[J]. Hydrogeology & Engineering Geology,2022,49(2):73 − 83. (in Chinese with English abstract)
[9]	向俐蓉,陈伟志,郭在旭,等. 大型滑坡抗滑桩-桩板结构受力变形研究[J]. 水文地质工程地质,2021,48(2):125 − 131. [XIANG Lirong,CHEN Weizhi,GUO Zaixu,et al. A study of the deformation of anti-slide pile and pile-plate structure in large landslide[J]. Hydrogeology & Engineering Geology,2021,48(2):125 − 131. (in Chinese with English abstract)
[10]	邓时容,肖世国. 嵌固段顶部拓宽型抗滑桩计算方法[J]. 中国地质灾害与防治学报,2022,33(4):84 − 91. [DENG Shirong,XIAO Shiguo. Calculation method of stabilizing piles with broadened top at the built-in section[J]. The Chinese Journal of Geological Hazard and Control,2022,33(4):84 − 91. (in Chinese with English abstract)
[11]	陶波,佴磊,伍法权,等. 论静止侧压力对抗滑桩的作用[J]. 吉林大学学报(地球科学版),2006,36(5):837 − 840. [TAO Bo,NAI Lei,WU Faquan,et al. The effect imposed on anti-slide pile by the static lateral pressure[J]. Journal of Jilin University (Earth Science Edition),2006,36(5):837 − 840. (in Chinese with English abstract)
[12]	肖世国. 考虑桩间土体抗滑作用的单排抗滑桩受力计算方法[J]. 中国地质灾害与防治学报,2020,31(1):89 − 94. [XIAO Shiguo. A calculation method for one row stabilizing piles considering resistance of soil between adjacent piles[J]. The Chinese Journal of Geological Hazard and Control,2020,31(1):89 − 94. (in Chinese with English abstract)
[13]	POULOS H G. Design of reinforcing piles to increase slope stability[J]. Canadian Geotechnical Journal,1995,32(5):808 − 818. doi: 10.1139/t95-078
[14]	钟昌茂,邱恩喜,魏永幸,等. 基于桩体抗剪强度的复合地基路堤稳定性分析方法[J]. 水文地质工程地质,2020,47(5):100 − 107. [ZHONG Changmao,QIU Enxi,WEI Yongxing,et al. Stability analysis method for the composite foundation embankment based on pile shear strength[J]. Hydrogeology & Engineering Geology,2020,47(5):100 − 107. (in Chinese with English abstract)
[15]	任青阳,赵梦园,谢忠伟,等. 抗滑桩应变特征与内力非线性研究[J]. 水文地质工程地质,2021,48(2):114 − 124. [REN Qingyang,ZHAO Mengyuan,XIE Zhongwei,et al. A study of the strain characteristics and internal force nonlinearity of anti-slide pile[J]. Hydrogeology & Engineering Geology,2021,48(2):114 − 124. (in Chinese with English abstract)
[16]	李宁,刘冠麟,许建聪,等. 降雨条件下抗滑桩边坡三维稳定性分析[J]. 中国地质灾害与防治学报,2018,29(3):38 − 46. [LI Ning,LIU Guanlin,XU Jiancong,et al. 3D numerical analysis of the stability of a slope reinforced with piles under rainfall conditions[J]. The Chinese Journal of Geological Hazard and Control,2018,29(3):38 − 46. (in Chinese with English abstract)
[17]	XIAO Shiguo. A simplified approach for stability analysis of slopes reinforced with one row of embedded stabilizing piles[J]. Bulletin of Engineering Geology and the Environment,2017,76(4):1371 − 1382. doi: 10.1007/s10064-016-0934-y
[18]	YAN Yuping,XIAO Shiguo. A calculation method for the embedded depth of stabilizing piles in reinforced slopes[J]. International Journal for Numerical and Analytical Methods in Geomechanics,2020,44(7):1077 − 1092. doi: 10.1002/nag.3052
[19]	熊治文. 深埋式抗滑桩的受力分布规律[J]. 中国铁道科学,2000,21(1):48 − 56. [XIONG Zhiwen. Force distribution rule of deeply buried anti slide pile[J]. China Railway Science,2000,21(1):48 − 56. (in Chinese with English abstract)
[20]	雷文杰, 郑颖人, 冯夏庭. 沉埋桩的有限元设计方法探讨[J]. 岩石力学与工程学报, 2006, 25（增刊1）: 2924 − 2929 LEI Wenjie, ZHENG Yingren, FENG Xiating. Study on finite element design methods of slope stabilized by deeply buried anti-slide piles[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25（Sup 1）: 2924 − 2929. （in Chinese with English abstract）
[21]	雷文杰,郑颖人,王恭先,等. 沉埋桩加固滑坡体模型试验的机制分析[J]. 岩石力学与工程学报,2007,26(7):1347 − 1355. [LEI Wenjie,ZHENG Yingren,WANG Gongxian,et al. Mechanism analysis of slope reinforcement with deeply buried piles with model test[J]. Chinese Journal of Rock Mechanics and Engineering,2007,26(7):1347 − 1355. (in Chinese with English abstract)
[22]	雷文杰,郑颖人,冯夏庭. 滑坡加固系统中沉埋桩的有限元极限分析研究[J]. 岩石力学与工程学报,2006,25(1):27 − 33. [LEI Wenjie,ZHENG Yingren,FENG Xiating. Limit analysis of slope stabilized by deeply buried piles with finite element method[J]. Chinese Journal of Rock Mechanics and Engineering,2006,25(1):27 − 33. (in Chinese with English abstract)
[23]	宋雅坤, 郑颖人, 雷文杰. 沉埋式抗滑桩机制模型试验数值分析研究[J]. 岩土力学, 2007, 28（增刊1）: 63 − 68 SONG Yakun, ZHENG Yingren, LEI Wenjie. Analysis of the experimental mechanism of deeply buried anti-slide pile by FEM[J]. Rock and Soil Mechanics, 2007, 28（Sup1）: 63 − 68. （in Chinese with English abstract）
[24]	The MathWorks Inc. Global optimization toolbox: User’s guide （R2018b）[M/OL]. （2018-09-01）[2022-11-10]. http://www.mathworks.com/help/pdf doc/gads/gads tb.pdf.
[25]	徐邦栋. 滑坡分析与防治[M]. 北京: 中国铁道出版社, 2001 XU Bangdong. Landslide analysis and prevention[M]. Beijing: China Railway Publishing House, 2001. （in Chinese）
[26]	铁道部第二勘察设计院. 抗滑桩设计与计算[M]. 北京: 中国铁道出版社, 1983 The Second Survey and Design Institute of Ministry of Railways. Design and calculation of anti-slide pile[M]. Beijing: China Railway Publishing House, 1983. （in Chinese）

施引文献

期刊类型引用(3)

1.	汤科，赵华，唐雪峰，程儒凇. 双排桩最优排距确定及桩间土压力计算方法. 水文地质工程地质. 2025(02): 104-113 . 本站查看
2.	尹平保，伍晗曦，贺炜，谭李俊英，刘江波，方敏. 推扭（H-T）荷载作用下刚性矩形基础承载特性研究. 水文地质工程地质. 2024(04): 125-134 . 本站查看
3.	马露，王坤，陶思源，程希莹，常山，朱学敏. 基于有效应力法的单桩负摩阻力解析解. 水文地质工程地质. 2024(05): 79-86 . 本站查看