1、书书书Journal of Engineering Geology工程地质学报10049665/2022/30(6)-1867-07胡新丽，徐楚，张杰豪 2022 多层滑带滑坡抗滑桩桩后推力分布计算方法研究J 工程地质学报，30(6):18671873 doi:1013544/jcnkijeg20220708Hu Xinli，Xu Chu，Zhang Jiehao 2022 Study on distribution and calculation of thrust on stabilizing pile of the multi-sliding zones landslideJ Journ
2、alof Engineering Geology，30(6):18671873 doi:1013544/jcnkijeg20220708多层滑带滑坡抗滑桩桩后推力分布计算方法研究*胡新丽徐楚张杰豪(中国地质大学(武汉)工程学院，武汉 430074，中国)摘要多层滑带滑坡演化过程具有多级滑体相对运动、应力非连续性特点，由于不同层滑体间的相对运动，多层滑带滑坡内的推力分布形式往往比较复杂，亟需提出针对多层滑带滑坡的抗滑桩桩后推力分布计算方法。本文采用水平极限平衡法对多层滑带滑坡进行分析，根据水平微段极限平衡状态下的受力分析，可计算抗滑桩桩后每层滑体内的滑坡推力。在分析竖直方向的受力平衡时，需要考虑
3、不同层滑体向下的受力传递，从而导致滑坡推力分布函数在滑带处存在突变。据此推导出了多层滑带滑坡的抗滑桩桩后推力分布函数，并通过数值模拟进行了工程案例验算分析，验算结果表明了该方法的可靠性。本文提出的多层滑带滑坡抗滑桩桩后推力分布计算方法，为多层滑带滑坡的防治理论提供了一定的参考。关键词多层滑带;滑坡推力;抗滑桩;水平极限平衡法;受力分析中图分类号:P642.22文献标识码:Adoi:1013544/jcnkijeg20220708*收稿日期:20221005;修回日期:20221205基金项目:国家自然科学基金重点国际(地区)合作研究项目(资助号:42020104006)This researc
4、h is supported by the NSFC Key Projects of International Cooperation and Exchanges(Grant No 42020104006)第一作者简介:胡新丽(1968)，女，博士，教授，博士生导师，主要从事地质灾害防治及岩土体稳定性评价等研究工作 E-mail:huxinlicugeducn通讯作者简介:徐楚(1994)，男，博士生，主要从事滑坡地质灾害防治研究工作 E-mail:xuchu cugeducnSTUDY ON DISTIBUTION AND CALCULATION OF THUST ON STABI-LIZ
5、ING PILE OF THE MULTISLIDING ZONES LANDSLIDEHU XinliXU ChuZHANG Jiehao(Faculty of Engineering，China University of Geosciences，Wuhan 430074，China)AbstractThe evolution process of landslides with multi-sliding zones is characterized by the relative motion ofmulti-layer and stress discontinuities Due t
6、o the mutual dislocation and pile-soil interaction between sliding bodiesin different layers，the thrust distribution of the multi-sliding zones landslide is often complex It is urgent to put for-ward the calculation method for the thrust distribution on stabilizing pile of multi-sliding zones landsl
7、ide In this pa-per，the horizontal limit equilibrium method is used to analyze the multi-sliding zones landslide According to thestress analysis under the limit equilibrium state of horizontal micro segment，the landslide thrust on stabilizing pilein each layer is obtained When analyzing the force bal
8、ance in the vertical direction，it is necessary to consider thedownward force transmission of sliding mass in different layers，resulting in the sudden change of landslide thrustdistribution function at the sliding zone The thrust distribution function on stabilizing pile of multi-sliding zoneslandsli
9、de is derived and verified by numerical simulation，which means that the function is reliable The thrust cal-culation method of multi-sliding zones landslide is proposed in this paper，which provides a certain idea for thetreatment of multi-sliding zones landslideKey wordsMulti-sliding zones;Landslide
10、 thrust;Stabilizing pile;Horizontal limit equilibrium method;Force a-nalysis0引言多层滑带滑坡是一种危害非常严重并且分布广泛的自然地质灾害，严重影响着我国人民的生命财产安全(李守定等，2004;黄润秋，2007;任伟中等，2009;代贞伟等，2016;向家松等，2017;Wanget al，2018;Zhang et al，2018;Tang et al，2019;郭长宝等，2022;Xu et al，2022;朱鸿鹄等，2022)。抗滑桩是目前开展滑坡治理使用最为广泛的防治结构之一，在我国滑坡地质灾害治理中广泛应用(
11、胡新丽等，2005;李长冬，2009;李新哲等，2021)。在进行抗滑桩设计时，需要确定滑坡推力的大小和分布(Ke et al，2006;Powrie et al，2007;戴自航等，2007;Zhang et al，2015;Ke et al，2019;Xiong etal，2019;Yan et al，2020)。多层滑带滑坡演化过程具有多级滑体相对运动、应力非连续性特点，滑坡推力分布和单层滑带滑坡不同，亟需对多层滑带滑坡推力计算方法进行深入研究。目前国内外研究学者主要通过理论计算和模型试验等方法确定滑坡的推力分布形式，已经取得了较多研究成果。戴自航(2002)通过理论推导和试验分析将滑坡
12、推力分布模式分为三角形、梯形和矩形 3 种;杨涛等(2006)通过假设圆弧组合面，使用Janbu 法计算滑坡推力;肖世国(2010)对似土质滑坡进行了水平微段极限平衡法的分析，推导出了滑坡推 力 的 计 算 公 式;郑 宏(2007)和 凌 道 盛 等(2013)采用迭代的方法对三维滑坡稳定性和滑坡推力进行了研究;Sharafi et al(2016)通过模型试验和数值模拟研究了坡顶加载条件下，不同因素对抗滑桩桩后推力分布的影响;陈鑫等(2019)基于土拱效应解释了滑坡推力呈“中间大、两头小”抛物线分布规律的原因;蔡强等(2016)通过模型试验表明钢管抗滑短桩桩后滑坡推力呈现“S”型分布;Zh
13、anget al(2018)基于有限差分方法并考虑不同力学参数的影响，研究了含多滑带的滑坡推力分布形式。这些研究成果在一定程度上解决了实际工程中滑坡推力计算的问题，但滑坡推力分布模式只有三角形、矩形、梯形或抛物线形等形式，而且不适用于多层滑带滑坡。综上所述，国内外研究学者对抗滑桩桩后推力进行了深入研究，但对多层滑带滑坡推力分布的相关研究还刚刚起步。单层滑带滑坡和多层滑带滑坡的滑坡推力在水平方向上的分布特征相似，在竖直方向上分布特征有所区别。相关学者一般认为竖直方向上的滑坡推力的分布图式为三角形、矩形、梯形或者抛物线形。由于不同层滑体间的相对运动，多层滑带滑坡抗滑桩所受滑坡推力分布规律和单层滑坡
14、抗滑桩所受滑坡推力分布有所区别，目前缺少对于多层滑带滑坡抗滑桩桩后推力分布的研究。本文通过改进水平极限平衡法，考虑了滑带处的受力传递，提出了多层滑带滑坡抗滑桩桩后推力分布计算方法，并通过数值模拟进行了工程案例验算分析，为多层滑带滑坡的防治理论提供了一定的参考。1理论模型建立如图 1 所示的理论分析模型，抗滑桩后侧滑体沿滑面滑移达到极限平衡状态，z 为沿桩侧的计算深度，方向竖直向下。自桩顶起向下沿着 z 轴方向取微段 dz，进行隔离体受力分析。为了便于分析问题，该计算模型采用如下假定:(1)在同一层滑体内，该微段的上下表面的剪应力大小 相等，根据肖世国(2010)的验证，可按库仑摩擦定律近似取为
15、 tan bz，一般取 0.20.4 范围内;(2)抗滑桩在滑坡演化过程中始终保持竖直方向，即桩后侧正应力 q 始终保持水平方向;(3)计算深度 z 采用等效重力法(选用底边为 y，高为 z 的矩形面积代替微段上面实际滑体重量图形面积);(4)设桩处滑带倾角不为 0。如图 1 所示的分析模型，以抗滑桩顶点为坐标轴原点，抗滑桩桩后侧水平方向为 y 轴正方向，竖直向下方向为 z 轴正方向。对于桩侧任一计算点，到滑带的水平距离为 y，计算深度 z 采用等效重力法(选用底边为 y，高为 z 的矩形面积代替微段上面实际滑体重量图形面积);W 为微段上部滑体的等效重力(在深层滑体中时，需要考虑浅层滑带对深
16、层滑体的作用力 1，1)，c0、0为桩和滑体界面的黏聚力和内摩擦角，q，q和，分别为微段对应的桩土界面和滑带上的正应力和剪应力，q和 取值按8681Journal of Engineering Geology工程地质学报2022图 1浅层滑体内滑坡推力分析模型Fig 1Thrust analysis model of landslide in shallow sliding massa 理论分析模型;b 微段受力模型式(1)。q=c0+qtan 0=c+tan(1)2滑坡推力公式推导由于不同层滑体间的相对运动，多层滑带启动顺序不同时，其滑坡推力分布存在一定差异。当深层滑带先启动时，浅层滑体内滑坡推力计算如下:根据图 1 所示的分析模型，对水平微段隔离体进行受力分析，当桩后滑体处于极限平衡状态时，根据水平和竖直两个方向的静力平衡条件，可以解出桩后侧面的正应力:Wy+qdz+ydw+dz=cot 1dz+(W+dW)(y dzcot 1)(2)y+dz+cot 1dz=qdz+(y dzcot 1)(3)式中:W=z。联立式(2)和式(3)可求得式(4)浅层滑体推力q1:q1=(Wcot
