第三讲：滑坡稳定性计算的几个问题

1、中国地质大学中国地质大学( (武汉武汉) )20132013年年7 7月月2222日日滑坡稳定性计算方法研究滑坡稳定性计算方法研究苏爱军内容内容一、滑坡稳定计算方法分类二、稳定性计算的几个基本概念二、稳定性计算的几个基本概念三、稳定性计算方法的改进三、稳定性计算方法的改进四、稳定性计算参数的确定四、稳定性计算参数的确定五、孔隙水压力计算五、孔隙水压力计算滑坡稳定性计算的几个问题滑坡稳定性计算的几个问题类型亚类备注确定性方法极限平衡法现行主要适应土体的条分法：如瑞典条分法，Bishop法，Janbu法，Spencer法，Morgenstern-Spencer法等；对岩土体均适应的条分法：如传递系

2、数法，改进的传递系数法，Hoek楔体分析法，Sarma法，实用递推法，精确递推法等弹塑性理论法塑性极限平衡法：适于土质斜坡，假定土体为理想刚塑性体，按Mohr-Coulomb屈服准则确定稳定系数数值分析法：适于岩土质斜坡，用弹塑(粘)有限元等数值法，计算斜坡应力分布状况，按Mohr-Coulomb破坏准则计算出破坏点和塑性区分布状况，据此确定稳定系数非确定性方法破坏概率法通过引入稳定性计算参数的分布概率（如C、概率分布等），计算滑坡稳定系数fos值的理论分布可靠度分析法稳定性计算时，通过引入各种荷载分项系数求得稳定系数fos值的可靠度（一）滑坡稳定计算方法分类 非确定性方法是以极限平衡法为基础

3、，通过引入稳定性计算参数的分布概率求得稳定系数的概率，从而判断滑坡的稳定性。与此同时，伴随近代土力学的发展，在确定边坡不同工作条件下的强度指标、孔隙水压力等方面也形成了一些行之有效的办法。这一完备的理论体系也随着非确定性分析的耦合而更加完善。近年来，三维力学模型的改善、可靠度理论的引入，使得边坡稳定性分析服务于工程的能力得到较大改善。 极限平衡分析法（包括现行的各种条分法）研究历史可追溯到上世纪20年代或更早。自上世纪三十年代Fellenius（1927，1936）提出土坝圆弧滑动稳定分析之后，经Bishop（1952），Morgenstern和Price（1965），Spencer（1966

4、），Janbu（1973），Sarma（1979），Chen（1983）等的改进，迄今已经形成能够满足静力、静力和力矩平衡的几十种滑坡稳定性评价方法，其中包括适用于各种形状滑动面的通用条分法。 一直以来，在工程界进行滑坡稳定性评价及滑坡治理的滑坡推力求解时，基本上都是采用基于极限平衡原理的条分法的各种方法。（二）极限平衡分析法 Petterson K.E首先提出了瑞典圆弧法，Fellenius W.和Taylor D.W.对其进行了改进。该方法最初在瑞典被采用，故此称为瑞典圆弧法（或瑞典条分法）。瑞典圆弧法的两个最重要的假定是：（1）滑动面为圆弧，稳定系数fos是围绕圆心的抗滑力矩与滑动力矩之

5、比值；（2）计算时不考虑土条之间的相互作用力，每个土条底部的反作用力可以直接由该土条上的荷载算出。 式中： 为滑动面上可能提供的抗滑力矩； 滑动力矩。0MMfosrrM0M 1、堤防工程设计规范（GB50286-98）8.2.4规定土堤抗滑稳定计算可采用瑞典圆弧滑动法（主要适用于在较为均质的土质边坡中进行最危险滑动面搜索并进行抗滑稳定计算，同时用于已知为圆弧型滑动面的滑坡的稳定性评价，适用该类滑坡的评价方法还有Bishop法、Spencer法等）。 2、建筑地基基础设计规范（PGB500072002）和岩土工程勘察规范（GB50021-2001）规定采用传递系数法。类似的评价方法有美国陆军师团

6、法、罗厄法等，主要是针对已知滑动面的滑坡而言的，可用于滑动面形状各异的土质、岩土质滑坡。Janbu法、普通条分法、Morgenstern & Spencer法（Morgenstern & Spencer，1965）等同样适用该类滑坡。Sarma法(1979)建立在斜分条基础上，是专门针对通常受软弱结构面切割的岩质边坡产生滑坡而开发的。对于受结构面控制的规模不大岩质边坡具有独特的优越性。方法应用条件及特征稳定图法（Janbu，1968；Duncan等，1987）对大多数分析来说，精度足够，比详细的计算机分析速度快瑞典条分法（Fellenius，1927，1936）圆弧滑面，定转动中心，各块间作用

7、合力平行滑面，满足力矩平衡，不满足水平、垂直力平衡改进Bishop法（Bishop，1955）圆弧滑面，拟合圆弧与转心，各块间作用力水平，条间切向力X为零，满足力矩平衡，满足垂直力平衡，不满足水平力平衡Janbu普通条分法（Janbu，1968）适用于任何形状滑动面，满足力矩平衡，满足垂直、水平力的平衡，允许侧向力方向变化，常遇到数值问题。简化计算条间切向力X0，再对稳定系数作修正Hoek楔体分析法(1974)楔形滑面，各滑面均为平面，以各滑面总抗滑力和楔体总下滑力确定稳定系数传递系数法(1977)圆弧或非圆弧滑面。条块间合力方向与上一条块滑面平行(Xi/Ei=tgi-1)Sarma法(197

8、9)非圆弧滑面或楔形滑面等复杂滑面。认为除平面和圆弧面外，滑体必先破裂成相互滑动的块体才能滑动，该方法以保证块体处于极限平衡状态为准，确定稳定系数Morgenstern & Spencer法（Morgenstern & Spencer，1965）适用于任何形状滑动面，满足力矩平衡，满足垂直、水平力的平衡，垂直与水平作用力之比（X/E）存在与水平方向坐标的函数关系(X/E=f(x)允许侧向力方向变化，常遇到数值问题Spencer法（Spencer，1965）圆弧滑面，或拟合中心圆弧。垂直与水平作用力之比（X/E）为一给定常值 主要条分法的特点类型亚类及分类原则 具体方法与假设条件 适用条件1 、

9、 精 确条 分 法（条块和整体满足静力与力矩平衡）（1）条块间水平作用力与垂直作用力之比满足一定的函数关系斯宾塞法（Spencer）、摩根斯坦普莱斯法和摩根斯坦陈祖煜法、简化潘一孙法等（假定条块间力的作用方向）主要适用土质滑坡，弧型滑动面。对非弧型滑动面存在收敛性问题（2）条块间作用力合力作用点的位置满足一定的函数关系简布法（Janbu法）(条块间作用力合力作用点位于条块高度的约三分之一处)适用土质、岩土质和岩质滑坡，任何形状滑动面，常遇到收敛性问题条分法的分类及其适用条件2 、半精确条分法（ 1 ）静力平衡（条块间水平作用力与垂直作用力之比满足一定的函数关系）不平衡推力法、传递系数法、改进的

10、传递系数法（假定条块间合力的作用方向平行于上一条块的滑动面）适用土质、岩土质和岩质滑坡，任何形状的滑动面Sarma法（假定条块间水平作用力与垂直作用力之比满足摩尔-库仑定理，垂直条分上限解法）简化的简布法（垂直条分下限解法，假定条块间垂直作用力等于零）美国陆军师团法（条块间合力倾角坡度和滑面平均倾角相等）罗厄法（条块间合力倾角与上一条块坡面和滑面倾角的相等）（ 2 ）力矩平衡或力矩平衡及部分静力平衡瑞典圆弧法（满足围绕滑动面圆心的力矩平衡，不考虑静力平衡和条块的力矩平衡）适用土质滑坡，圆弧型滑动面毕肖普法1（满足围绕滑动面圆心的力矩平衡，条块间垂直作用力平衡）适用土质滑坡，圆弧型滑动面，个别情

11、况不收敛毕肖普法2（条块和整体满足力矩平衡，条块间垂直作用力平衡）适用土质、岩土质和岩质滑坡，任何形状的滑动面，个别情况不收敛 表 2 条分法各种方法比较表 所满足的平衡条件 计算方法 整体力矩 土条力矩 垂直力 水平力 滑面形式 瑞典圆弧法 圆弧 简化毕肖普法 圆弧 简布法 任意 斯宾塞法 任意 摩根斯坦 普赖斯法 任意 沙尔玛法 任意 不平衡推力法 任意 内容内容一、一、滑坡稳定计算方法分类滑坡稳定计算方法分类二、稳定性计算的几个基本概念三、稳定性计算方法的改进三、稳定性计算方法的改进四、稳定性计算参数的确定四、稳定性计算参数的确定五、孔隙水压力计算五、孔隙水压力计算滑坡稳定性计算方法研究

12、滑坡稳定性计算方法研究（一）稳定系数与安全系数有关滑坡或边坡稳定性评价及其防治方面的文献，对于稳定系数与安全系数的概念应用比较混乱，许多文献将安全系数与稳定系数混为一谈。早期很多学者将稳定系数称之为稳定安全系数，如毕肖普等将土坡沿整个滑裂面的抗剪强度f与实际产生的剪应力之比定义为稳定安全系数，而后来的学者，如陈祖煜2、钱家欢3等直接称其为安全系数。因此，有必要厘清稳定系数与安全系数的概念，以免给初学者的学习带来困难。1、稳定系数 就单纯满足静力平衡条件而言，稳定系数可定义为沿可能滑动面上的抗剪强度 累积值与剪应力累积值之比值。在既满足静力平衡同时又满足力矩平衡的条件下，还应同时满足在此稳定系数

13、条件下可能滑动岩土体对其重心的力矩平衡。fcf 从理论上讲，滑坡体或边坡在特定荷载下的应力、变形关系可以采用合适的数学模型如实表达（例如应力应变间的非线性性质，卸荷增荷性质，屈伏破坏性质等等）。通过合适的数学模型精确合理地给出各点应力、变形值，从而得到滑坡或边坡稳定问题的精确解。 就滑坡这种抗剪问题而言，采用合适的数学模型可以得到每一点上的剪应力，并求得每一点上的抗剪强度f ，同时可算出每一点上的局部稳定系数。如果每一点上的稳定系数fos均大于1，整个滑坡从抗剪切方面来说是安全的。如果有个别点的稳定系数fos小于1，则这些部位的将自动等于f，应力向临近点转移，这些部位进入屈伏（甚至破坏）状态。

14、由于这些屈伏区是孤立的，局限于小范围内，没有形成连贯的破坏面，因此，在特定荷载下该滑坡仍是稳定的。相反，如果屈伏的部位已经连成一个连贯性的破坏面，甚至得不到一个满足平衡要求的解答，说明该滑坡在特定荷载下已不能维持稳定。 基于刚体极限平衡理论，采用条分法求解稳定系数可通过三个途径实现： 强度储备法； 超载法； 准超载法。 （1）强度储备法 强度储备法是将岩土体的抗剪强度参数（、c）除以fos，并将fos逐渐增大或缩小。每取一个fos值就进行一次稳定性分析，直到滑坡或边坡沿某一连贯性可能滑动面达到所设定的平衡条件为止。这个fos值就是沿某一连贯性可能滑动面产生滑动的稳定系数。所求解稳定系数的方法称

15、为强度储备法。 （2） 超载法 超载法是将荷载乘以fos，并将fos逐渐增大或缩小，每取一个fos值就进行一次稳定性分析，直到滑坡或边坡沿某一连贯性可能滑动面达到所设定的平衡条件为止。这个fos值就是沿某一连贯性可能滑动面产生滑动的稳定系数。所求解稳定系数的方法称为超载法。 强度储备法的稳定系数fos1与超载法的稳地系数fos2不相等，与准超载法稳定系数fos2相等： 强度储备法的稳地系数fos1： 1=（Wcostg/fos1 +c L /fos1）/ Wsin 超载法的稳定系数fos2： 1=（fos2Wcostg+c L ）/ fos2 Wsin 其中：Fos1、Fos2分别为强度储备法

16、、超载法的滑坡稳定系数； W为滑体重量（kN）； 为滑动面倾角（）； 为滑动面内摩擦角（） C为滑动面内聚力（kPa）； L为滑动面长度(m)。WCLtgWfosCLfosCLtgWCLtgWfosfosfoscos1coscos22212CLtgWfosfosCLCLtgWfosCLtgWfosfoscos11coscos22221 则由上式可知，当滑带土的内摩擦角和内聚力均不等于零时，强度储备法和超载法的滑坡稳定系数是不相等的。对于滑带土为非粘土（C=0 ）的滑坡，fos2=fos1 成立的充分必要条件是滑坡处于临界状态（fos2=1，即= ），否则 。因此超载法也不适合计算滑带土为非粘土

17、（ C=0）的滑坡的稳定系数。只有滑带土为纯粘土（ =0）时，强度储备法和超载法的滑坡稳定系数才相等（ ，fos2=fos1 ）。综上所述，除了极个别条件下（ =0），强度储备法和超载法的滑坡稳定系数相等外，用超载法求得的滑坡稳定系数与用强度储备法求得的滑坡稳定系数既不相等，二者之比也不是一个确定的常数，超载法的概念不太明确，不确定性较大。 （3） 准超载法 为克服准超载法的上述不足，本人以强度储备法为基础，综合储备法稳定性计算与超载法剩余推力计算的习惯与优点，提出了准超载法的新概念。 滑坡稳定性计算中的准超载法是沿某一连贯性可能滑动面各点将切向荷载乘以fos，并将fos逐渐增大或缩小，每取一

18、个fos值就进行一次稳定性分析，直到滑坡或边坡该可能滑动面达到所设定的平衡条件为止。这个fos值就是沿某一连贯性可能滑动面产生滑动的稳定系数。所求解稳定系数的方法称为准超载法。 强度储备法的稳定系数fos1与超载法的稳地系数fos2不相等，与准超载法稳定系数fos3相等： 强度储备法的稳地系数fos1： 1=（Wcostg/fos1 +c L /fos1）/ Wsin 超载法的稳定系数fos2： 1=（fos2Wcostg+c L ）/ fos2 Wsin 准超载法的稳定系数fos3： 1=（Wcostg+c L ）/ fos3 Wsin 其中：Fos1、Fos2、 Fos3分别为强度储备法、

19、超载法和准超载法的滑坡稳定系数； W为滑体重量（kN）； 为滑动面倾角（）； 为滑动面内摩擦角（） C为滑动面内聚力（kPa）； L为滑动面长度(m)。W 准超载法的稳定系数与强度储备法的稳地系数相等，但准超载法的剩余推力P2是强度储备法剩余推力P1 的fos倍： 1=（Wtg+c L）/(fosP1) =（Wtg/fos+c L /fos）/ P1 剩余推力：P1WP2= fosP1WP2= fosP12 、安全系数 由于滑坡体或边坡岩土体的的不均匀性、物理力学性质的各向异性以及作用荷载不确定性等原因，实际上我们难以用一个或几个简单的数学模型精确合理地给出各点应力、变形值，从而得到滑坡或边坡

20、稳定问题的精确解。某些情况下当稳定系数fos大于1时，滑坡在特定荷载下可能仍然是不稳定的。因此，我们必须根据计算过程中的荷载组合、所假定的边界条件和不同的计算方法以及保护对象的安全等级等，给出相对应的、大于1的经验值K（安全系数）。当稳定系数fosK时，我们认为滑坡体或边坡不会沿整个滑裂面产生整体滑动，滑坡体或边坡有足够的安全储备，是安全的；当稳定系数fosK时，不能保证滑坡体或边坡不会沿整个滑裂面产生整体滑动，此时定义滑坡体或边坡是不安全的。此K值为经验值，而非计算值。我们将其定义为安全系数。 1、准超载法的特点 由准超载法定义及表达式可知，强度储备法和准超载法的稳定系数值相同，但在极限平衡

21、条件下计算所得到的沿滑动面的滑动力和抗滑力不一样。准超载法沿滑动面的滑动力和抗滑力强度储备法的fos倍。同理，采用垂直条分法进行滑坡推力计算时，准超载法计算所得到的滑坡推力是强度储备法的k倍（k为给定安全系数）。（二）定义准超载法的意义？ 以往滑坡防治工程设计所采用的滑坡设计阻滑力是由强度储备法计算所得到的滑坡推力乘以荷载系数来确定的。荷载系数与所防治的滑坡安全等级有关，但荷载系数确定的的相关条件不够明确，因此，设计人员的经验不同，同等情况下所给出的荷载系数可以不同，一般情况下荷载系数取 1.2。准超载法将与滑坡安全等级密切相关的荷载系数与安全系数等同起来，使各种计算方法之间具有可比性，在很大

22、程度上克服了荷载系数使用时的不确定性，其考虑问题的出发点与现行岩土工程勘察规范（GB50021-2001）相关规定相吻合。2、准超载法的应用价值？（三）（三）现行条分法的现行条分法的3个共性问题个共性问题 众多的垂直条分法通过引入条块间作用力假定的方式来求解本质上静不定问题，其合理性一直受到人们普遍关注。Morgenstern& Price最早提出了解的合理性限制，认为所获得的解必须满足以下条件： （1）条块间不产生拉力； （2）作用于土条界面上的剪力不超过按摩尔库仑法则提供的抗剪强度。 同时，Morgenstern & Price还提出了这样的推论：不同的关于土条作用力的假定，只要满足上述两

23、个合理性条件限制，相应的稳定系数彼此相差不大。 然而，大量的工程实践表明，滑动面形状不规则（滑动面起伏较大）的岩土质和岩质滑坡，由于人为认定的最危险滑动面的不准确性，或条块间作用力相关假定的不合理性等原因，采用现有条分法进行稳定性评价时，常常会出现以下三个方面的问题： 1、滑体条块间出现拉应力（Morgenstern&Price提出的合理性限制条件1得不到满足），计算所得到的滑体整体稳定系数偏大 滑体条块间出现拉应力的原因可以是多方面的：其一，勘察所认定的滑动面起伏较大，并不是统一的最危险滑动面，滑坡体实际上由两大部份或几大部份组成。各部分可能滑动面形态不同，稳定程度各异，因而导致所计算的一些

24、条块间出现拉应力。其二，由于外部荷载施加的部位和方式不一，本应为整体性的滑坡演变为前缘局部牵引式滑坡。滑坡在某种状态下达到分段平衡，靠近剪出口段稳定系数小，远离剪出口段稳定系数大，所计算的两段结合部位一些条块间必然出现拉应力。其三，所假定的条块间作用力方向与实际作用力方向有明显偏差。 2、当滑动面出现反倾时，一些方法可能出现后一滑块超覆于前一滑块的滑动 当滑动面出现反倾时，无论是不平衡推力法还是是传递系数法，均会出现后一滑块超覆于前一滑块的滑动。 罗厄法和简化的简布法在某些情况下也存在同样的问题。从物理意义上说，这是不合理的。满足静力与力矩平衡条件的现行条分法的各种方法或不收敛，或得到的解不合

25、理。滑动面 3、同时满足静力与力矩平衡条件的现行条分法的各种方法对于滑动面形状不规则（滑动面起伏较大）的岩土质、岩质滑坡，尤其是大型岩质滑坡均不同程度的存在收敛性问题（不能保证在合理条件下的无条件收敛） 上述三个问题在很大程度上是由于现行条分法的各种方法计算公式常常难以满足上述两个合理性限制条件所决定的。要解决好三个问题，必须对条块间作用力的合力方向随横坐标的变化进行深入研究，提出比现行各种条分法更合理的条块间作用力合力方向随横坐标或滑动面倾角变化的函数关系。 对于第一问题，通过增加相应的必要条件，在提高现行满足静力平衡条件的条分法的各种计算方法的精度的同时，使问题得到解决。 对于第一和第二个

26、问题同时存在的情况，通过构筑特定的滑体条块间剩余下滑力的合力倾角的分布函数，建立了应用范围更广泛、精度更高的实用递推法，使问题得到圆满解决。 通过增加相应的必要条件并构筑合适的滑体条块间剩余下滑力的合力倾角的分布函数后，可以解决第三个问题。内容内容一、一、滑坡稳定计算方法分类滑坡稳定计算方法分类二、稳定性计算的几个基本概念二、稳定性计算的几个基本概念三、稳定性计算方法的改进四、稳定性计算参数的确定四、稳定性计算参数的确定五、孔隙水压力计算五、孔隙水压力计算滑坡稳定性计算方法研究滑坡稳定性计算方法研究（一）从传递系数法说起 传递系数法是当前我国规程规范应用最广的方法。为什么没有得到国际公认！i-

27、1iii+1iiiii-1ENQTREi-1ii+1nNi=QicosTi=Qisin视滑坡后缘至剪出口n-1条块滑动面处于极限状态（稳定系数为1），并达到静力平衡，将n-1个条块的剩余推力全部投影到第n条块滑动面上，然后将第n条块滑动面上投影的抗滑力除以滑动力，由此求得滑坡稳定系数。传递系数法视滑坡后缘至剪出口n-1条块滑动面处于极限状态（稳定系数为1），并达到静力平衡。考虑第n-1条块对第n条块的作用力和第n条块自身的重力、渗透压力等，将第n条块滑动面上抗滑力除以滑动力，由此求得滑坡稳定系数Fos 。其结果是，滑坡后缘至剪出口n-1条块的稳定系数为1，不同于第n条块的稳定系数Fos ！（1

28、）静力平衡的原则不一样！i-1iii+1iiiii-1ENQTREi-1ii+1nNi=QicosTi=Qisin其它条分法则不同，将n个条块滑动面上的剪切强度参数tg、C全部除稳定系数Fos ，并使n个条块沿滑动面上全部处于静力平衡，此时得Fos 即为滑坡稳定系数。其结果是n个条块滑动面的稳定系数都是相等的。 （2）传递系数法展示了滑坡极限平衡的不完整过程，稳定系数的大小与最后一个条块的划分密切相关，具有较大的随意性。其它条分法反映的是滑坡整体平衡状态，最后一个条块的划分方式影响不大。（3）传递系数法在科学意义方面的不足 从某种程度上说，传递系数法考虑滑坡渐进破坏的过程，更贴近实际一些。 但

29、是，传递系数法没有严格地遵从国际通行的基于刚体极限平衡理论的稳定系数定义。从科学研究的角度来看，其具有横向不可对比性和相对不确定性。 从技术角度，传递系数法的上述差异只是导致计算的稳定系数偏大，通过安全系数的调整，基本不影响在工程方面的应用！ （二）改进的传递系数法1、改进原则：静力平衡的原则与国际接轨，即基于刚体极限平衡理论的稳定系数定义原则与国际接轨滑坡整体稳定系数与各条块稳定系数相等；进行滑坡稳定系数和滑坡推力计算时，计算机能自动识别条块间出现拉应力的部位，并自动改正；引入准超载法，统一荷载系数和安全系数。i-1iii+1iiiii-1ENQTREi-1ii+1nNi=QicosTi=Q

30、isin0111111111111aiiiakikjjiiakikjjiinainnijjinainnijiiERfosTRTfosETTRRfos2、改进的传递系数法的基本公式j=cos(i-i+1) - sin(i-i+1)tani+1/fos改进的传递系数法的传递系数：j=cos(i-i+1) - sin(i-i+1)tani+1/fos传递系数法的传递系数：j=cos(i-i+1) - sin(i-i+1)tani+1 tgsin - )cos( =j1+j1+jj1+jj111)(111)(ninijTnjiTninijRnjiRfos传递系数法改进的传递系数法的滑坡推力公式En=

31、kw sin+ k En-1 cos(n-1 n )-sin (n-1 n )tg / k -(wcostg+cl) = kwsin+ En-1 cos(n-1 n )-sin (n-1 n )tg /k -(wcostg+cl) 传递系数法的滑坡推力公式 En= kw sin+ En-1 cos(n-1 n)-sin (n-1 n)tg - (wcostg+cl) 二者的差别在于用前者用sin (n-1 n )tg /k代替了后者的sin (n-1 n )tg ，后者是前者的简化计算方法。由此导致了实际应用时二者4方面差异： （1）公式本身的差别 当条块间作用力出现负值时，改进的传递系数法能

32、通过计算机自动识别并改正，传递系数法做不到这一点； 当条块间作用力均大于零时，改进的传递系数法与传递系数法稳定性计算式形式相同，但传递系数不同：改进传递系数法的传递系数：j=cos(i-i+1) - sin(i-i+1)tani+1/fos传递系数法的传递系数：j=cos(i-i+1) - sin(i-i+1)tani+1 （2）给定剪切强度参数C 、 ，当稳定系数小于1时，传递系数法计算所得到的稳定系数小于改进的传递系数法特解计算所得到的稳定系数；当稳定系数大于1时，传递系数法计算所得到的稳定系数大于改进的传递系数法特解计算所得到的稳定系数。稳定系数与1差别越大，改进的传递系数法特解与传递系

33、数法计算所得到的稳定系数差别也越大。一般情况下，当稳定系数小于1.5时，二者相对差不超过2%。 （3）给定稳定系数进行反演分析、求解滑动面剪切强度参数时，给定稳定系数等于1，二者求得的剪切强度参数完全一致；对于正在蠕变的滑坡，给定稳定系数小于1，传递系数法求得的剪切强度参数比改进的传递系数法特解求得的剪切强度参数大；对于处于稳定的滑坡，给定稳定系数大于1，传递系数法求得的剪切强度参数比改进的传递系数法特解求得的剪切强度参数小。 （4）给定安全系数，传递系数法计算的滑坡推力小于改进的传递系数法特解计算的滑坡推力。给定的安全系数越大，传递系数法计算的滑坡推力越小。一般情况下，当安全系数小于1.25

34、时（安全等级为级的滑坡基本工况所对应的安全系数），二者相对差不超过15%。滑坡推力越大，二者相对差越小。当滑坡推力大于1000kN时，传递系数法计算所得到的滑坡推力一般大于改进的传递系数法特解90， 当滑坡推力很小时（如单宽滑坡推力小于250kN/m），二者相对差将超过15%。 滑 方 项 坡 法 目石佛寺千将坪磨子坪寨坝山水美地改 进 传 系 数法稳定系数1.1331.0661.2991.1961.109滑 坡 推 力（kPa/m）13364.9828565.78138.411126.58传系数法稳定系数1.1351.0701.3121.2111.113滑 坡 推 力（kPa/m）12237

35、.1927794.9490.831067.32相对差稳定系数241%1%4滑坡推力8%3%34%5%典型滑坡稳定系数及最后条块的滑坡推力（精度1，安全系数1.25）Table 1 List of result about stability coefficient and thrusting force of typical landslides 滑体条块间不出现拉应力的情况滑体条块间出现拉应力的情况110150190230110150190230高程(m)0100200300400(m)F26F25F24F23F22F21F20F19J3P长江学校1234567891011121314151

36、718192021地下水位重庆市丰都南竹园滑坡计算剖面示意图Fig. 2 Sketch map of calculated section 50100150200250300350400450-10000100020003000400050006000P=0 未考虑负的剩余下滑力(安 全系数取1.2) 考虑了负的剩余下滑力(安 全系数取1.2)剩余推力(KPa)滑坡水平距离(m)改进传递系数法和传递系数法滑坡推力对比图（安全系数1.2）Fig. 3 Contrast chart of thrusting force of improved transfer coefficient method

37、 and transfer coefficient method计算方法传递系数法改进传递系数法稳定系数1.15011.0792最后条块滑坡推力(kP)918.91839.0表2 稳定系数及滑坡推力 (安全系数取1.20)Table 2 List of result about stability coefficient and thrusting force从计算结果看，由于条块出现了负的剩余下滑力，传递系数法计算的稳定系数比改进传递系数法的稳定系大0.071。取安全系数1.2时，传递系数法计算得到的各条块的滑坡推力均比改进传递系数法小，最后一条块的滑坡推力小50。WiNiTiEiPiXiX

38、i-1Ei-1Pi-1iQiqiPtiPdiTi-NiiLibi i（一）半精确条分法的通用公式Uti-1UtiUiuti-1utiui地 下 水 面 线 Hwi+1Hwi滑坡条块划分及受力分析滑坡条块划分及受力分析滑坡条块划分及地下水压力分析图滑坡条块划分及地下水压力分析图其中：稳定系数计算递推公式（苏爱军，2008年） ： 累计求和公式表达式：式中： 为第 i条块水平推力（kN）； 为第 i条块垂直推力（kN）； 为第i条块滑体的重力（kN）； 为第i条块地表垂直荷载（kN）； 为第i条块所受的水平地震荷载（kN）； 为第i条块所受的人工抗滑力（kN）； 为第i条块滑面的倾角（）； 为第i

39、条块所受的人工抗滑力的倾角（）； 为条块间第 i条块作用力合力的倾角（）； 为第i条块靠剪出口侧侧面水压力（kN）； 为第i条块滑面水的扬压力（kN）； 为第i条块滑带土的内摩擦角（）； 为第i条块滑带土的粘聚力(kPa)； 为第i条块滑动面长度(m)。式中：fos为稳定系数滑面组合实用递推法改进的传递系数法传递系数法陆军师团法罗厄法简化的毕肖普法简化Janbu下限解上限解勘察的滑面工况11.050 1.064 1.067 1.054 1.064 1.031 1.030 1.062工况20.972 0.985 0.986 0.976 0.986 0.954 0.953 0.986花坪滑坡花坪滑

40、坡C各种工况下的稳定性计算结果各种工况下的稳定性计算结果部分半精确条分法和精确条分法稳定系数计算结果比较Table 7 Stability coefficients calculated by some semi-accurate slice methods and accurate slice methods 计算模型半精确条分法精确条分法采用数值模拟值计算SU法实 用递 推法传递系数法改进的传递系数法罗厄法美国陆军工程师团法简 化的 简布法斯 宾赛法摩 根斯坦-陈法实际*1.4921.479 1.490 1.5721.5251.5061.5021.384 1.4851.463*1.2691

41、.265 1.274 1.3291.3061.2891.2851.183 1.2711.252简化*1.4241.416 1.422 1.5021.4591.4451.4341.316 1.4181.396*1.2111.211 1.216 1.2681.2481.2361.2261.125 1.2131.193 计算结果表明：（1）半精确条分法中，自传递系数法、改进的传递系数法1、罗厄法、美国陆军工程师团法、实用递推法、SU法至简化的简布法，稳定系数逐渐减小；计算结果稳定系数最大的传递系数法比最小的简化的简布法大12.3%13.6%。（2）上述方法中，仅SU法所计算的稳定系数均位于精确条分法

42、中的斯宾赛法和摩根斯坦-陈法所计算的稳定系数之间。以斯宾赛法和摩根斯坦-陈法所计算的稳定系数平均值为标准，SU法所计算的稳定系数，实际模型仅大0.28%0.34%，简化模型仅大0.64%0.674%。因此，SU法计算精度很高，完全满足工程精度要求。（3）SU法和实用递推法所计算的稳定系数与将数值模拟值计算得到的条块间作用力倾角代入力的平衡方程计算所得到的稳定系数最接近。其中，SU法与后者的稳定系数相差0.08%0.75%；实用递推法与后者的稳定系数相差0.07%0.42%。因此，计算滑坡稳定系数时，实用递推法的计算精度也满足工程精度要求。（三）精确条分法的通用公式WiTiEiPiXiXi-1E

43、i-1Pi-1iQiqiPtiPdiNiiLibi iHiHpiHi-1HdiUigiei滑坡条块受力分析图滑坡条块受力分析图Uti-1UtiUiuti-1utiui地下水面线Hti-1HtiHwiHwi-1滑坡条块地下水压力分析图滑坡条块地下水压力分析图精确条分法递推公式（苏爱军，2008年） （3-35）累计求和公式表达式：iHiHtiHdigie式中：-第i条块左侧面水压力作用点到滑面的高度（m）；-第i条块左侧面条块间作用力作用点到滑面的高度（m）；-第i条块水平地震力作用点到滑面的高度（m）；-第i条块滑动面水压力作用点到滑面中点的力臂长度（m）-第i条块滑体及其上部重力作用点到滑面

44、中点的力臂长度（m）。实用递推法通用公式其中：内容内容一、一、滑坡稳定计算方法分类滑坡稳定计算方法分类二、稳定性计算的几个基本概念二、稳定性计算的几个基本概念三、稳定性计算方法的改进三、稳定性计算方法的改进四、稳定性计算参数的确定五、孔隙水压力计算五、孔隙水压力计算滑坡稳定性计算方法研究滑坡稳定性计算方法研究 滑带剪切强度参数的确定方法有三种滑带剪切强度参数的确定方法有三种-试验法、试验法、反演分析法和经验法反演分析法和经验法 在条件许可时，首先考虑反演分析在条件许可时，首先考虑反演分析方法，方法，c值按值按滑坡厚度确定，滑坡厚度确定，值用反演分析结果；其次，值用反演分析结果；其次，试样具试样

45、具有代表性，试验结果可靠时，可以将饱和固结快剪有代表性，试验结果可靠时，可以将饱和固结快剪值乘值乘0.85-0.9，c值乘值乘0.30-0.4作为设计剪切强度指标作为设计剪切强度指标；经验值仅作为参考。上述剪切强度均为有效应力状；经验值仅作为参考。上述剪切强度均为有效应力状态下的取值。态下的取值。 1、试验方法、试验方法 目前试验确定c、值的方法很多，既可通过原位试验测定，也可以在试验室内通过剪切试验测定，原位试验的特点是简捷、快速，并大体上在原位应力和边界条件下进行。缺点是边界条件无法精确确定和控制，这给试验结果的解释带来了困难，常用的原位试验中，除十字板剪切试验可直接测定饱和软粘土的不排水

46、强度外，其它如标准贯入试验，静力触探试验和旁压试验等均属间接测试方法，常带有地区性或经验性，相反，室内试验的主要优点是边界条件比较明确，且容易控制，饱和试样中孔隙水压力的测试成功，有效应力原理的证实，大大推进和提高了对土强度特性的试验研究和认识。但不论哪种试验，都要求被测土的抗动程度最小，室内试验，除某些探索强度基本规律的试验外，总要采取土样，取样将不可避免地引起应力释放和土的结构扰动。因此，试验结果的可靠性将首先取决于取样的质量。如果土样业已剧烈扰动，那么精细的试验将失去意义，因为再精细的试验也不能补偿土样扰动的影响。 该方法最大的不足是由于滑带土的不均匀性导致试验结果离散性太大而难以应用。

47、尽管目前对采取部位进行了限定，如分别对滑坡前、中、后三部分的滑带土进行采取，试验尽可能模拟各部位滑坡的受力特征等，仍难以达到满意的效果，从三峡库区滑坡勘察成果分析，值的变异系数0.30-0.5，c值的变异系数0.35-0.6。值增大，c值减小，而c值增大值减小的现象。对和平广场8个土层采用竖井取大样进行室内饱和固结快剪试验表明，同一个大样做的试验结果，尽管具有较好的线性相关性，但仍然存在因物质成份的差异、导致试验结果出现偏差的现象，由此可见在大范围内取样试验结果的代表性是较差的。对库水位突降引起的滑坡，可采用固结不排水剪指标，利用总应力法进行稳定性及推力计算；对于形成稳定渗流场的滑坡，模拟稳定

48、渗流期，可采用固结排水剪强度指标，采用有效应力法进行稳定性及推力计算。一般性情下，考虑地下水浮力及渗透压力的情况下，均可采用固结排水剪试验的剪切强度指标，利用有效应力法进行稳定性及推力计算。QHJ2-1 深度3.0m QHJ2-1 深度3.0mQJ4-1深度4.5-4.7m QJ4-1深度4.5-4.7m QHJ1-2深度7.5m QHJ1-2深度7.5m QHJ3-2深度4.7-4.9m QHJ3-2深度4.7-4.9m QHJ3-1深度3.0-3.2m QHJ3-1深度3.0-3.2m QHJ2-2深度0.5m QHJ2-2深度0.5m QHJ1-4 深度9.2-9.4m QHJ1-4 深

49、度9.2-9.4m 排除试验样代表性所带来的误差我们可以得到以下规律： 1)就三峡库区而言，对同一试验方法（如饱和固结快剪），在有效主压应力150-700kPa时（相当于饱水时土层厚度12-60m），主压应力（1）越大，值越小、c值越大，即摩尔包线是曲线而非直线的特征不明显，摩尔包线基本上为直线。 2)对于侏罗系紫红色碎屑岩中的滑坡，当主压应力1超过700Kpa后，抗剪强度增加不明显既c值明显增大，而值明显减小。 3)饱和固结快剪的值与排水反复剪的残余强度值相差较小，但内聚力c值前者比后者大2-3倍，有的甚至大4倍。 4)饱和固结快剪值比常规三轴排水剪值大2-5。根据上述特点，若采用试验方法确

50、定滑带土的剪切强度指标，首先应保证取样的代表性并尽可能少扰动，应尽可能采用竖井勘探方法取大样，对同一组样最好同时进行快剪、饱和快剪、饱和固结快剪及三轴排水剪与不排水剪试验，以便进行不同试验方法的结果比较分析。2、反演分析对于处于蠕滑阶段的滑坡，一般假定稳定系数等于0.98（特殊的如基岩滑坡考虑滑坡局部应力集中、产生渐进破坏，可假定稳定系数等于0.9），利用刚体极限平衡理论方法进行c、值参数的反演。一般情况下，当滑体厚度小于5m时，c值的影响比值大，而厚度大于5m的滑坡值影响比c值大，且厚度越大，值影响越大。在进行反演分析时一般假定c值求值。山田刚二等在滑坡和斜坡崩坍及其防治一书中指出，由土体无

51、侧限抗压试验测得qu，然后通过qu求得c，c值与滑体厚度呈正比关系。同时他还指出当试验值与上述经验值相差太大时，最好不采用试验值。岩土工程手册（中国建筑工程出版社）亦持相同观点。我委在三峡库区滑坡防治工程设计时基本上是采用上述方法求得滑带土的c、值，再参照试验经验值及试验指标。在进行反演分析时应特别注意以下几点：1)应尽可能地模拟滑坡蠕滑时的边界条件，尤其是地下水水位，如果难以做到，则可取勘探时雨季最高地下水位，这时的计算结果相对偏安全的。2)选择的反演分析剖面与主滑剖面一致，否则成果过于偏安全。3)用作反演分析的刚体极限平衡理论方法，应与设计用的稳定性及推力计算方法一致，一般采用设计规范所推

52、荐的不平衡推力法。按滑坡厚度选取C值滑坡厚度(m)C 值(kPa)55101015152020 注：引自岩土工程手册3、经验方法岩土工程手册推荐的经验公式及数据如下：1)我国铁路沿线粘性土滑带，经多次直剪试验得出，残余强度与塑性指数Ip，液性指数IL的经验关系式如下：2=2.4278-1.2229Ip-0.1172IL相关系数为0.86。2)我国宝成线上大型堆积层滑坡，其黑色炭质页岩风化的粘土滑带c=20Kpa，=712；紫红色泥岩风化之粘土滑带c=10Kpa，=57；黄土质粉质粘土滑带c=20Kpa，=1316。3)陕西南部裂隙粘土质滑带c=68Kpa，=8。岩土工程勘察设计手册引用日本专家

53、系统数据认为滑坡斜面坡度与滑带土c、值相关（详见表1）。从三峡库区滑坡勘察防治成果看，侏罗系中下统红层中的崩坡积层滑坡，从三峡库区滑坡勘察防治成果看，侏罗系中下统红层中的崩坡积层滑坡，滑带土的内摩擦角一般滑带土的内摩擦角一般1720；三叠系上统沙镇溪组及侏罗系上统蓬来；三叠系上统沙镇溪组及侏罗系上统蓬来镇组崩坡积及残坡积土层滑坡滑带土的内摩擦角一般镇组崩坡积及残坡积土层滑坡滑带土的内摩擦角一般713；三叠系巴东；三叠系巴东组中滑坡滑带土内摩擦角组中滑坡滑带土内摩擦角1619；二叠系嘉陵江组及更古老的碳酸盐岩；二叠系嘉陵江组及更古老的碳酸盐岩中的滑带土的内摩擦角一般中的滑带土的内摩擦角一般182

54、2。c值与滑坡厚度有关。值与滑坡厚度有关。滑坡斜面与滑带土C、值的关系滑坡斜面坡度()实例数()C(kPa)相关系数010469.05.800.944101513214.83.400.951152012720.70.400.93520259523.64.400.94125304027.94.000.930302730.07.300.944 注：引自岩土工程勘察设计手册内容内容一、一、滑坡稳定计算方法分类滑坡稳定计算方法分类二、稳定性计算的几个基本概念二、稳定性计算的几个基本概念三、稳定性计算方法的改进三、稳定性计算方法的改进四、稳定性计算参数的确定四、稳定性计算参数的确定五、孔隙水压力计算滑坡

55、稳定性计算方法研究滑坡稳定性计算方法研究（一）滑体及滑带中的地下水位及孔隙（一）滑体及滑带中的地下水位及孔隙水压力水压力“治滑先治水”，滑体与滑带土中地下水储存状态及孔隙水压力一直是滑坡研究中的一个重要研究课题。目前比较统一的观点是滑带粘土颗粒间的孔隙非常小，而且土体中孔隙所占的容积比在砂中的要大得多，在这种粘土层中的上面或下面若存在很大水压的含水层时，这种压力是否能提高孔隙水压力尚有不少疑问（因为实际上水并不是完全的流体，而具有某种程度的粘性）。但是用含水层的水位来代替滑带土中的水位，在某种程度上与实际情况还是较吻合的。为了探讨上述观点的准确性，在进行三峡库区和平广场滑坡勘察过程中，本人专门

56、布置了三个地下水位长观孔，对滑坡堆积体中的地下水进行分层观测（和平广场滑坡堆积体的渗透系数为i10-6i10-5cm/sec）。具体步骤是这样的，堵塞滑面以下的孔移段，对滑坡堆积体中的地下水安装双层观测孔，观测滑带及其以上3m的地下水位与潜水面地水位的差值及随时间的变化。通过一个多月的观测表明，潜水面比滑带附近的水位高出0.53m，二水位随时间变化具有同步性，土层渗透系数越低水位差越大。图2-3 分层观测管结构图3铁纱铁纱砂砾石卡盘粘土3（孔径110mm）2外层观测管1外层观测管（直径91mm）粘土砂砾石（直径32mm）30cm1-5m1230m滑坡堆积体基岩图 万 州和平广场滑坡群地下水位长观孔观测曲线( a ) AZ3地下水位历时曲线0.000.501.001.502.002.506月2 1日7月6 日7月1 2日7月1 5日7月1 8日7月2 9日8月3 日观测日期地下水位埋深（m ）( b ) AZ23 地下水位历时曲线0.001.002.003.004.005.006.007.008.007月3 日7月6 日7月1 2日7月1 5日7月1 8日7月2 9日8月3 日观测日