多层滑带滑坡抗滑桩桩后推力分布计算方法研究_胡新丽.pdf

书书书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 Journalof Engineering Geology，30(6):18671873 doi:1013544/jcnkijeg20220708多层滑带滑坡抗滑桩桩后推力分布计算方法研究*胡新丽徐楚张杰豪(中国地质大学(武汉)工程学院，武汉 430074，中国)摘要多层滑带滑坡演化过程具有多级滑体相对运动、应力非连续性特点，由于不同层滑体间的相对运动，多层滑带滑坡内的推力分布形式往往比较复杂，亟需提出针对多层滑带滑坡的抗滑桩桩后推力分布计算方法。本文采用水平极限平衡法对多层滑带滑坡进行分析，根据水平微段极限平衡状态下的受力分析，可计算抗滑桩桩后每层滑体内的滑坡推力。在分析竖直方向的受力平衡时，需要考虑不同层滑体向下的受力传递，从而导致滑坡推力分布函数在滑带处存在突变。据此推导出了多层滑带滑坡的抗滑桩桩后推力分布函数，并通过数值模拟进行了工程案例验算分析，验算结果表明了该方法的可靠性。本文提出的多层滑带滑坡抗滑桩桩后推力分布计算方法，为多层滑带滑坡的防治理论提供了一定的参考。关键词多层滑带;滑坡推力;抗滑桩;水平极限平衡法;受力分析中图分类号:P642.22文献标识码:Adoi:1013544/jcnkijeg20220708*收稿日期:20221005;修回日期:20221205基金项目:国家自然科学基金重点国际(地区)合作研究项目(资助号:42020104006)This research 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-LIZING 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 to 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 landslide 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 balance 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 zoneslandslide 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 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)。抗滑桩是目前开展滑坡治理使用最为广泛的防治结构之一，在我国滑坡地质灾害治理中广泛应用(胡新丽等，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)通过理论推导和试验分析将滑坡推力分布模式分为三角形、梯形和矩形 3 种;杨涛等(2006)通过假设圆弧组合面，使用Janbu 法计算滑坡推力;肖世国(2010)对似土质滑坡进行了水平微段极限平衡法的分析，推导出了滑坡推 力 的 计 算 公 式;郑 宏(2007)和 凌 道 盛 等(2013)采用迭代的方法对三维滑坡稳定性和滑坡推力进行了研究;Sharafi et al(2016)通过模型试验和数值模拟研究了坡顶加载条件下，不同因素对抗滑桩桩后推力分布的影响;陈鑫等(2019)基于土拱效应解释了滑坡推力呈“中间大、两头小”抛物线分布规律的原因;蔡强等(2016)通过模型试验表明钢管抗滑短桩桩后滑坡推力呈现“S”型分布;Zhanget al(2018)基于有限差分方法并考虑不同力学参数的影响，研究了含多滑带的滑坡推力分布形式。这些研究成果在一定程度上解决了实际工程中滑坡推力计算的问题，但滑坡推力分布模式只有三角形、矩形、梯形或抛物线形等形式，而且不适用于多层滑带滑坡。综上所述，国内外研究学者对抗滑桩桩后推力进行了深入研究，但对多层滑带滑坡推力分布的相关研究还刚刚起步。单层滑带滑坡和多层滑带滑坡的滑坡推力在水平方向上的分布特征相似，在竖直方向上分布特征有所区别。相关学者一般认为竖直方向上的滑坡推力的分布图式为三角形、矩形、梯形或者抛物线形。由于不同层滑体间的相对运动，多层滑带滑坡抗滑桩所受滑坡推力分布规律和单层滑坡抗滑桩所受滑坡推力分布有所区别，目前缺少对于多层滑带滑坡抗滑桩桩后推力分布的研究。本文通过改进水平极限平衡法，考虑了滑带处的受力传递，提出了多层滑带滑坡抗滑桩桩后推力分布计算方法，并通过数值模拟进行了工程案例验算分析，为多层滑带滑坡的防治理论提供了一定的参考。1理论模型建立如图 1 所示的理论分析模型，抗滑桩后侧滑体沿滑面滑移达到极限平衡状态，z 为沿桩侧的计算深度，方向竖直向下。自桩顶起向下沿着 z 轴方向取微段 dz，进行隔离体受力分析。为了便于分析问题，该计算模型采用如下假定:(1)在同一层滑体内，该微段的上下表面的剪应力大小 相等，根据肖世国(2010)的验证，可按库仑摩擦定律近似取为 tan bz，一般取 0.20.4 范围内;(2)抗滑桩在滑坡演化过程中始终保持竖直方向，即桩后侧正应力 q 始终保持水平方向;(3)计算深度 z 采用等效重力法(选用底边为 y，高为 z 的矩形面积代替微段上面实际滑体重量图形面积);(4)设桩处滑带倾角不为 0。如图 1 所示的分析模型，以抗滑桩顶点为坐标轴原点，抗滑桩桩后侧水平方向为 y 轴正方向，竖直向下方向为 z 轴正方向。对于桩侧任一计算点，到滑带的水平距离为 y，计算深度 z 采用等效重力法(选用底边为 y，高为 z 的矩形面积代替微段上面实际滑体重量图形面积);W 为微段上部滑体的等效重力(在深层滑体中时，需要考虑浅层滑带对深层滑体的作用力 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 1+c1 c0)(1+tan 1cot 1)cot 1 tan 1+(1+tan 1cot 1)tan 0+(c1)cot 1(cot 1 tan 1)cot 1 tan 1+(1+tan 1cot 1)tan 0(4)当深层滑带先启动时，深层滑体内滑坡推力计算如下:根据图 2 所示的分析模型，对水平微段隔离体进行受力分析，当桩后滑体处于极限平衡状态时，根据水平和竖直两个方向的静力平衡条件，可以解出桩后侧面的正应力:Wy+qdz+ydw图 2深层滑体内滑坡推力分析模型Fig 2Thrust analysis model of landslide in deep sliding massa 理论分析模型;b 微段受力模型=cot 1dz+(W+dW)(y dzcot 1)+dz(5)y+dz=qdz+(y dzcot 1)+cot 1dz(6)联立式(5)和式(6)可求得式(7)深层滑体推力q2:q2=(Wcot 2 c2 c0)(1 tan 2cot 2)tan 2+cot 2+(1 tan 2cot 2)tan 0+(c2)cot 2(cot 2+tan 2)tan 2+cot 2+(1 tan 2cot 2)tan 0(7)对于深层滑体来说，上部受到的重力可以等效为浅层滑带对下部作用的正应力和剪应力竖直向下的分量，以此做等效代换，如图 3。图 3W 等效分析图Fig 3W equivalent analysis diagram可求得等效 W 为式(8):W=(z2 z1)+cos 1sin 1(8)联立式(6)和式(1)可解出:=q1cot 1 ccot 11+tan cot 1(9)=c+q1cot 1 ccot 11+tan cot 1tan(10)将式(9)和式(10)代入式(8)可求得:W=(z2 z1)c1sin 1)(1+tan 1cot 1)1+tan 1cot 1tan bcot 1(cos 1 tan 1sin 1)+968130(6)胡新丽等:多层滑带滑坡抗滑桩桩后推力分布计算方法研究(q1 c1cot 1)(cos 1 tan 1sin 1)1+tan 1cot 1tan bcot 1(cos 1 tan 1sin 1)(11)式中:c0，0为桩和滑体界面的黏聚力和内摩擦角;c1，1，1和 c2，2，2分别为浅层滑带和深层滑带的黏聚力、内摩擦角和计算点水平向对应滑面上点的切线倾角;b1，b2分别为浅层滑体和深层滑体的内摩擦角。当浅层滑带先启动时，浅层滑体内滑坡推力计算如下:图 4浅层滑体内滑坡推力分析模型Fig 4Thrust analysis model of landslide in shallow sliding massa 理论分析模型;b 微段受力模型根据如图 4 所示的分析模型，对水平微段隔离体进行受力分析，当桩后滑体处于极限平衡状态时，根据水平和竖直两个方向的静力平衡条件，可以解出桩后侧面的正应力:Wy+qdz+ydw=cot 1dz+(W+dW)(y dzcot 1)+dz(12)y+dz=qdz+(y dzcot 1)+cot 1dz(13)式中:W=z。联立式(12)和式(13)可求得式(14)浅层滑体推力 q1:q1=(Wcot 1 c1 c0)(1 tan 1cot 1)tan 1+cot 1+(1 tan 1cot 1)tan 0+(c1)cot 1(cot 1+tan1)tan 1+cot 1+(1 tan 1cot 1)tan 0(14)当浅层滑带先启动时，深层滑体内滑坡推力计算如下:根据如图 5 所示的分析模型，对水平微段隔离体进行受力分析，当桩后滑体处于极限平衡状态时，根据水平和竖直两个方向的静力平衡条件，可以解图 5深层滑体内滑坡推力分析模型Fig 5Thrust analysis model of landslide in deep sliding massa 理论分析模型;b 微段受力模型出桩后侧面的正应力:Wy+qdz+ydw=cot 1dz+(W+dW)(y dzcot 1)+dz(15)y+dz=qdz+(y dzcot 1)+cot 1dz(16)联立式(15)和式(16)可求得式(17)深层滑体推力 q2:q2=(Wcot 2 c2 c0)(1 tan 2cot 2)tan 2+cot 2+(1 tan 2cot 2)tan 0+(c2)cot 2(cot 2+tan 2)tan 2+cot 2+(1 tan 2cot 2)tan 0(17)对于深层滑体来说，上部受到的重力可以等效为浅层滑带对下部作用的正应力和剪应力竖直向下的分量，以此做等效代换，如图 6 所示。图 6W 等效分析图Fig 6W equivalent analysis figure可求得等效 W 为式(18):W=(z2 z1)+cos 1+sin 1(18)联立式(13)和式(1)可解出:=q cot 1+ccot 11+tan cot 1(19)=c+q cot 1+ccot 11+tan cot 1tan(20)将式(19)和式(20)代入式(18)可求得:0781Journal of Engineering Geology工程地质学报2022W=(z2 z1)+c1sin 1)(1+tan 1cot 1)1+tan 1cot 1+tan bcot 1(cos 1+tan 1sin 1)+(q1+c1cot 1)(cos 1+tan 1sin 1)1+tan 1cot 1+tan bcot 1(cos 1+tan 1sin 1)(21)基于该方法，将两层滑带滑坡推广到 n 层滑带滑坡，可得到如下结论:当多层滑带滑坡的 n 层滑带由深到浅依次启动时，第 i 层滑体内滑坡推力公式。当 1 i n 时:qi=(Wcot i+ci c0)(1+tan icot i)cot i tani+(1+tan icot i)tan 0+(ci)cot i(cot i tan i)cot i tani+(1+tan icot i)tan 0(22)当 i=n 时:qi=(Wcot i ci c0)(1 tan icot i)tan i+cot i+(1 tan icot i)tan 0+(ci)cot i(cot i+tan i)tan i+cot i+(1 tan icot i)tan 0(23)当多层滑带滑坡的 n 层滑带由浅到深依次启动时，第 i 层滑体内滑坡推力为:qi=(Wcot i ci c0)(1 tan icot i)tan i+cot i+(1 tan icot i)tan 0+(ci)cot i(cot i+tan i)tan i+cot i+(1 tan icot i)tan 0(24)3工程案例验算为了验证前文多层滑带滑坡推力分布计算方法的可靠性，选取马家沟多层滑带滑坡为研究对象，将理论计算与数值模拟计算结果进行对比分析。本文多层滑带滑坡数值模拟的主要研究目的是得到滑坡推力分布和抗滑桩内力分布特征。根据现场监测数据(张玉明，2018)，使用 FLAC3D软件采用局部三维模型计算，建立计算模型如图 7 所示，滑坡模型长625 m，宽 22 m，高 216 m，共划分 71 195 个节点和87 600 个网格单元。滑坡重点研究区域为多层滑带和设桩位置附近，在此区域附近加密网格重点分析。模型的边界条件设定为速度为 0 的固定边界，四周及底面均设置为法向约束，坡表不设约束。根据现场布设情况，在模型中部建立 3 根矩形抗滑桩，桩长 40 m，截面尺寸为 2.0 m3.0 m，具体尺寸和位置如图 7 所示。抗滑桩位置浅层滑带、中层滑带及深层滑带深度分别为 12 m、18 m 和 26 m。计算模图 7马家沟多层滑带滑坡局部三维数值模型Fig 7Local three dimensional numerical model of Majiagoumulti-sliding sliding zone landslide-pile system型材料参数如表 1 所示，取值参考了土工试验数据和前人研究成果(张玉明，2018)。考虑数值模拟计算的复杂性，本文未能考虑地下水作用的影响。表 1计算模型材料物理力学参数取值表Table 1Values of physical and mechanical parametersof calculation model materials岩土材料类型密度/kg m3黏聚力/kPa内摩擦角/()体积模量/Pa剪切模量/Pa孔隙率滑体211035279310646106040浅层滑带211010156710628106035中层滑带245015182810715107025深层滑带245010142810715107025基岩25002000507510840108020抗滑桩25001610912109表 2计算模型桩土接触面参数取值表Table 2Parameter values of pile-soil interfaceof calculation model接触面类型法向刚度/kPa m1切向刚度/kPa m1黏聚力/kPa内摩擦角/()滑体-抗滑桩151071510728216浅层滑带-抗滑桩26106261066412中层滑带-抗滑桩851068510612144深层滑带-抗滑桩85106851068112滑床-抗滑桩5010850108160040采用实体单元建立抗滑桩，在桩土接触界面采用接触面单元进行模拟，抗滑桩和岩土体接触面参数如表 2 所示。通过数值模拟计算，得到了马家沟多层滑带滑坡推力分布图式，如图 8 所示。根据数值模拟计算178130(6)胡新丽等:多层滑带滑坡抗滑桩桩后推力分布计算方法研究图 8多层滑带滑坡推力计算结果对比图Fig 8Comparison of thrust calculation results ofmulti-layer sliding zone landslide结果，通过监测滑带处位移可知滑坡在演化过程中滑带整体由深到浅依次启动，同时滑坡位移由前缘逐步向后缘发展。因此水平微段极限平衡法理论计算中设定 3 层滑带由深到浅依次启动;通过多次试算，当理论计算取=0.35 时，理论计算结果和数值模拟计算结果基本一致，计算结果如图 8 所示。数值模拟计算结果显示多层滑带滑坡推力随深度增加而不断增大，在浅层滑体内符合三角形分布，在中层滑体和深层滑体内滑坡推力数值虽有波动，但基本符合梯形分布，波动原因可能是数值模型网格划分不够精细。由图 8 可知，在浅层滑带处滑坡推力突然增大，在中层滑带处滑坡推力突然减小。其原因是中层滑体相对于浅层滑体向下滑动，浅层滑体下方受到顺坡向的摩擦力，因此浅层滑带处推力突然增大。由图 7 可知，马家沟滑坡中层滑带为中前缘局部滑带，在桩前与深层滑带形成了整体滑动，同时滑坡位移前缘逐步向后缘发展，桩后区域中层滑体相对于深层滑体顺坡向滑动，中层滑带处受到逆坡向的摩擦力，因此中层滑带处滑坡推力突然减小。由数值模拟和水平微段极限平衡法理论计算两种结果对比可知，两者抗滑桩桩后推力分布模式与大小基本相似，说明了本文方法的可靠性。4结论(1)在前人研究的基础上，分析得出了多层滑带滑坡推力分布函数存在分段的特征;根据水平微段极限平衡状态下的受力分析，通过竖直方向和水平方向的受力平衡求解出了每层滑体内的滑坡推力。(2)分析竖直方向的受力平衡时，需要考虑不同层滑体向下的受力传递，导致滑坡推力分布函数在滑带处存在突变;通过理论推导将两层滑带滑坡推力公式推广到 n 层滑带滑坡中，得到了多层滑带滑坡抗滑桩桩后推力计算分布函数。(3)选取马家沟多层滑带滑坡为工程案例，将水平微段极限平衡法理论与数值模拟计算结果进行了对比分析，对比结果验证了本文多层滑带滑坡抗滑桩桩后推力分布计算方法的可靠性，为多层滑带滑坡的防治理论提供了一定的参考。参考文献Cai Q，Li Q K，Shi S W，et al 2016 Study of mechanical characteristicsof short anti-sliding steel pipe pile by physical model testJ ockand Soil Mechanics，37(S2):679684Chen X，Xiang X C，Liu K，et al 2019 Thrust distribution law of anti-slide pile under small pile spacingJ Geological Science andTechnology Information，38(6):157164Dai Z W，Yin Y P，Wei Y J，et al 2016 Deformation and failure mecha-nism of outang landslide in Three Gorges eservoir areaJ Journalof Engineering Geology，24(1):4455Dai Z H，Chen L J 2007 Two numerical solutions of laterally loaded pilesinstalled in multi-layered soils by m methodJ Chinese Journal ofGeotechnical Engineering，29(5):690696Dai Z H 2002 Study on the distribution law of landslide-thrust and re-sistance of sliding mass acting on antislide pilesJ Chinese Journalof ock Mechanics and Engineering，21(4):517521Guo C B，Zhang Y S，Liu D T，et al 2022 Centrifuge model test of reacti-vationmechanismofJiangdingyaancientlandslideinGansuProvinceJ Journal of Engineering Geology，30(1):164176Hu X L，Wang L Q，Tang H M，et al 2005 Several key problems for opti-mum design on anti-slide pile of bank landslide in Three Gorgeseservoir regionJ Geological Science and Technology Informa-tion，24(S1):121125Huang Q 2007 Large-scale landslides and their sliding mechanisms inChina since the 20thcenturyJ Chinese Journal of ock Mechanicsand Engineering，26(3):433454Ke W，Liu Q J，Zhang C 2019 Kinematic bending of single piles in lay-ered soilJ Acta Geotechnica，14(1):101110Ke Y，Liang 2006 Numerical solution for laterally loaded piles in atwo-layer soil profileJ Journal of Geotechnical Geoenvironmen-tal Engineering，132(11):14361443Li C D 2009 Study on interaction mechanism between anti-slide pile andlandslide mass and pile optimizationD Wuhan:China Universityof GeosciencesLi S D，Li X，Zhang N X，et al 2004 Sedimentation characteristics of theJurassic sliding-prone stratum in the Three Gorges eservoir area andtheir influence on physical and mechanical properties of rockJ2781Journal of Engineering Geology工程地质学报2022Journal of Engineering Geology，12(4):385389Li X Z，Li T，Pan D，et al 2021 Optimal reinforcement position of anti-slide piles based on stress level of sliding zone J Journal of Engi-neering Geology，29(3):640646Ling D S，Qi S C，Chen F，et al 2013 A limit equilibrium method basedon Morgenstern-price method for 3D slope stability analysisJ Chi-nese Journal of ock Mechanics and Engineering，32(1):107116Powrie W，Smethurst J 2007 Monitoring and analysis of the bending be-haviour of discrete piles used to stabilise a railway embankmentJGotechnique，57(8):663677en W Z，Bai S W，Tang X J，et al 2009 esearch on deformation andfailure mechanisms of landslide with double sliding plane and rock-layer reverseJ Chinese Journal of ock Mechanics and Engineer-ing，28(3):552559Sharafi H，Sojoudi Y 2016 Experimental and numerical study of pile-sta-bilized slopes under surface load conditions J International Jour-nal of Civil Engineering，14(4):221232Tang H M，Wasowski J，Juang C H 2019 Geohazards in the Three GorgeseservoirArea，China-LessonslearnedfromdecadesofresearchJ Engineering Geology，261:105267Wang J，Su A，Liu Q，et al 2018 Three-dimensional analyses of the slid-ing surface distribution in the Huangtupo No 1 riverside sliding massin the Three Gorges eservoir area of China J Landslides，15(7):14251435Xiang J S，Wen B P，Chen M，et al 2017 Activity response of a landslidewith complex structure to fluctuation of reservoir water level:A casestudy of the Chaiwan Landslide in the Three Gorges eservoirJHydrogeology Engineering Geology，44(4):7177Xiao S G 2010 Approximate theoretical solution of distribution modes oflandslidethrustonanti-slidingpilesinsoil-likeslopesorlandslidesJ Chinese Journal of Geotechnical Engineering，32(1):120123Xiong X，Shi Z，Xiong Y，et al 2019 Unsaturated slope stability aroundthe Three Gorges eservoir under various combinations of rainfalland water level fluctuation J Engineering Geology，261:105231Xu C，Hu X，Niu L，et al 2022 Physical model test of the deformationbehaviorandevolutionaryprocessofthemulti-slidingzonelandslideJ Bulletin of Engineering Geology and the Environment，81(10):115Yan Y P，Xiao S G 2020 A calculation method for the embedded depthof stabilizing piles in reinforced slopes J International Journal forNumerical and Analytical Methods in Gemechanics，44(7):10771092Yang T，Zhou D P，Zhang J Y，et al 2006 Distribution of landslide thruston anti-slide pilesJ Chinese Journal of Geotechnical Engineering，28(3):322326Zhang L，Zhao M，Zou X 2015 Behavior of laterally loaded piles in mul-tilayered soilsJ International Journal of Geomechanics，15(2):06014017Zhang Y M 2018 Multi-field characteristics and evolution mechanism ofMajiagoulandslide-stablizingpilessystemunderreservoiroperationsD Wuhan:China University of GeosciencesZhang Y，Hu X，Tannant D D，et al 2018 Field monitoring and deforma-tion characteristics of a landslide with piles in the Three Gorges es-ervoir areaJ Landslides，15(3):581592ZhengH2007Arigorousthree-dimensionallimitequilibriummethodJ Chinese Journal of ock