土坡稳定性分析.docx

第七章土坡稳定性分析第一节概述土坡就是由土体构成、具有倾斜坡面的土体，它坡眉坡皿的简单外形如图7-1所示。一般而言，土坡有两种类坡'型。由自然地质作用所形成的土坡称为天然土坡，如顶/誓山坡、江河岸坡等；由人工开挖或回填而形成的土坡直底一圣L称为人工土 （边）坡，如基坑、土坝、路堤等的边坡。土坡在各种内力和外力的共同作用下，有可能产生剪图7-1 土坡各部位名称切破坏和土体的移动。如果靠坡面处剪切破坏的面积很大，则将产生一部分土体相对于另一部分土体滑动的现象，称为滑坡。土体的滑动一般 系指土坡在一定范围内整体地沿某一滑动面向下和向外移动而丧失其稳定性。除设计或施 工不当可能导致土坡的失稳外，外界的不利因素影响也触发和加剧了土坡的失稳，一般有 以下几种原因：1. 土坡所受的作用力发生变化：例如，由于在土坡顶部堆放材料或建造建筑物而使坡 顶受荷。或由于打桩振动，车辆行驶、爆破、地震等引起的振动而改变了土坡原来的平衡 状态；2. 土体抗剪强度的降低：例如，土体中含水量或超静水压力的增加；3. 静水压力的作用：例如，雨水或地面水流入土坡中的竖向裂缝，对土坡产生侧向压 力，从而促进土坡产生滑动。因此，粘性土坡发生裂缝常常是土坡稳定性的不利因素，也 是滑坡的预兆之一。在土木工程建筑中，如果土坡失去稳定造成塌方，不仅影响工程进度，有时还会危及 人的生命安全，造成工程失事和巨大的经济损失。因此，土坡稳定问题在工程设计和施工 中应引起足够的重视。天然的斜坡、填筑的堤坝以及基坑放坡开挖等问题，都要演算斜坡的稳定性，亦既比 较可能滑动面上的剪应力与抗剪强度。这种工作称为稳定性分析。土坡稳定性分析是土力 学中重要的稳定分析问题。土坡失稳的类型比较复杂，大多是土体的塑性破坏。而土体塑 性破坏的分析方法有极限平衡法、极限分析法和有限元法等。在边坡稳定性分析中，极限 分析法和有限元法都还不够成熟。因此，目前工程实践中基本上都是采用极限平衡法。极 限平衡方法分析的一般步骤是：假定斜坡破坏是沿着土体内某一确定的滑裂面滑动，根据 滑裂土体的静力平衡条件和莫尔一库伦强度理论，可以计算出沿该滑裂面滑动的可能性， 即土坡稳定安全系数的大小或破坏概率的高低，然后，再系统地选取许多个可能的滑动面， 用同样的方法计算其稳定安全系数或破坏概率。稳定安全系数最低或者破坏概率最高的滑 动面就是可能性最大的滑动面。本章主要讨论极限平衡方法在斜坡稳定性分析中的应用，并简要介绍有限元法的概念。第二节无粘性土坡稳定性分析无粘性土坡即是由粗颗粒土所堆筑的土坡。相对而言，无粘性土坡的稳定性分析比较 简单，可以分为下面二种情况进行讨论。一、均质的干坡和水下坡均质的干坡系指由一种土组成，完全在水位以上的无粘性土坡。水下土坡亦是由一种 土组成，但完全在水位以下，没有渗透水流作用的无粘性土坡。在上述二种情况下，只要 土坡坡面上的土颗粒在重力作用下能够保持稳定，那么，整个土坡就是稳定的。在无粘性土坡表面取一小块土体来进行分析(图7-2)，设该小块土体的重量为W，其法 向分力N = Wcosa，切向分力T = Wsina。法向分力产生摩擦阻力，阻止土体下滑，称为 抗滑力，其值为R = Ntg9 =Wcosa tg中。切向分力T是促使小土体下滑的滑动力。则土体的稳定安全系数F为：S抗滑力 R W cos a tg。tg。(7-1)F =s 滑动力 T W sin atga式中：。土的内摩擦角(° );a 土坡坡角(° )。图7-2无粘性土坡由上式可见，当a =。时，Fs=l，即其抗滑力等于滑动力，土坡处于极限平衡状态， 此时的a就称为天然休止角。当a。时，土坡就是稳定的。为了使土坡具有足够的安全 储备，一般取Fs=1.11.5。二、有渗透水流的均质土坡图7-3渗透水流逸出的土坡当边坡的内、外出现水位差时，例如基坑排水、坡外水位下降时，在挡水土堤内形成渗流场，如果浸润线在下游坡面逸出(图7-3)，这时，在浸润线以下，下游坡内的土体除了受到重力作用外，还受到由于水的 渗流而产生的渗透力作用，因而使下游边坡的稳定性降低。渗流力可用绘流网的方法求得。作法是先绘制流网，求滑弧范围内每一流网网格的平 均水力梯度，从而求得作用在网格上的渗透(流)力：(7-2)J =y iA式中：Y*水的重度;Ai网格的面积。求出每一个网格上的渗透力Ji后，便可求得滑弧范围内渗透力的合力“。将此力作为 滑弧范围内的外力(滑动力)进行计算，在滑动力矩中增加一项：AM，= Tjlj(7-3)式中：'4距圆心的距离。如果水流方向与水平面呈夹角。，则沿水流方向的渗透力j = y i。在坡面上取土体V 中的土骨架为隔离体，其有效的重量为y'V。分析这块土骨架的稳定性，作用在土骨架上 的渗透力为J = jV =y wiV。因此，沿坡面的全部滑动力，包括重力和渗透力为T =y'Vsina+y iVcos(a0)(7-4)坡面的正压力为*N =y' V cos a -y iV sin(a - 0)(7-5)*则土体沿坡面滑动的稳定安全系数：F = N tg = y' V cos a - y *iV sin(a - 0 )tg任s T y'Vsina+y iVcos(a -0)*式中：i渗透坡降；y '土的浮重度；y *水的重度；"土的内摩擦角。若水流在逸出段顺着坡面流动，即。=a。这时，流经路途ds的水头损失为曲，所 以，有(7-7). dh .-I = = sin a ds将其代入式(76)，得：F = y E (7-8) 'sat '土坡稳定安全系数降低y'/ysat。因此，要保持同图7-4渗透水流未逸出的土坡由此可见，当逸出段为顺坡渗流时， 样的安全度，有渗流逸出时的坡角比没有 渗流逸出时要平缓得多。为了使土坡的设 计既经济又合理，在实际工程中，一般要 在下游坝址处设置排水棱体，使渗透水流 不直接从下游坡面逸出(图7-4)。这时的 下游坡面虽然没有浸润线逸出，但是，在下游坡内，浸润线以下的土体仍然受到渗透力的作用。这种渗透力是一种滑动力，它将降 低从浸润线以下通过的滑动面的稳定性。这时深层滑动面如图7-4中虚线表示)的稳定性可 能比下游坡面的稳定性差，即危险的滑动面向深层发展。这种情况下，除了要按前述方法 验算坡面的稳定性外，还应该用圆弧滑动法验算深层滑动的可能性。第三节 粘性土坡的稳定性分析一般而言，粘性土坡由于剪切而破坏的滑动 面大多数为一曲面，一般在破坏前坡顶先有张裂 缝发生，继而沿某一曲线产生整体滑动。图7-5 中的实线表示一粘性土坡滑动面的曲面，在理论 分析时可以近似地将其假设为圆弧，如图中虚线 表示。为了简化计算，在粘性土坡的稳定性分析 中，常假设滑动面为圆弧面。建立在这一假定上 的稳定性分析方法称为圆弧滑动法。这是极限平 衡方法的一种常用分析方法。一、整体圆弧滑动法瑞典的彼得森(K.E.Petterson)于1915年采用圆弧滑 动法分析了边坡的稳定性。此后，该法在世界各国的土 木工程界得到了广泛的应用。所以，整体圆弧滑动法也 被称为瑞典圆弧法。如图7-6,表示一个均质的粘性土坡，它可能沿圆 弧面AC滑动。土坡失去稳定就是滑动土体绕圆心0发 生转动。这里把滑动土体当成一个刚体，滑动土体的重图7-5粘性土坡的滑动面量W为滑动力，将使土体绕圆心O旋转，滑动力矩图7-6整体圆弧滑动受力示意图Ms=W d (d为通过滑动土体重心的竖直线与圆心。的水平距离)。抗滑力矩Mr由两部分组成：滑动面AC上粘聚力产生的抗滑力矩，值为 caCR；滑动土体的重量W在滑动面上的反力所产生的抗滑力矩。反力的大小和方向与土的内摩擦角中值有关。当中=0时，滑动面是一个光滑曲面，反力的方向必定垂直于 滑动面，即通过圆心。，它不产生力矩，所以，抗滑力矩只有前一项caACR。这时， 可定义粘性土坡的稳定安全系数为:F =抗滑力矩=m r = c AC R(7-9)s 一滑动力矩mR Wd此式即为整体圆弧滑动法计算边坡稳定安全系数的公式。注意，它只适用于。=0的情况。若。力0，则抗滑力与滑动面上的法向力有关，其求解可参阅下面的条分法。二、瑞典条分法所谓瑞典条分法，就是将滑动土体竖直分成若干个土条，把土条看成是刚体，分别求 出作用于各个土条上的力对圆心的滑动力矩和抗滑力矩，然后按公式(7-9)求土坡的稳定安 全系数。把滑动土体分成若干个土条后，土条的两个侧面分别存在着条块间的作用力(图7-7)。作用在条块i上的力，除了重力叫外，条块侧面ac和bd上作用有法向力P、Pi+i，切向 力H.、H.+，法向力的作用点至滑动弧面的距离为h.、h.+。滑弧段cd的长度;其上作 用着法向力质和切向力T.，T.包括粘聚阻力c.l.和摩擦阻力Ntg ©。考虑到条块的宽度111111i不大，Wj和N.可以看成是作用于cd弧段的中点。在所有的作用力中，P、H1在分析前一 土条时已经出现，可视为已知量，因此，待定的未知量P1+i. H1+i. h"、N1和T15个。 每个土条可以建立三个静力平衡方程，即，Fx1=o, Z F0和? M1=和一个极限平如果把滑动土体分成n个条块，则n个条块之间的分界面就有(n-1)个。分界面上的未 知量为3(n-1)，滑动面上的未知量为2n个，还有待求的安全系数匚，未知量总个数为(5n-2)， 可以建立的静力平衡方程和极限平衡方程为4n个。待求未知量与方程数之差为(n-2)。而一 般条分法中的n在10以上。因此，这是一个高次的超静定问题。为使问题求解，必须进行 简化计算。瑞典条分法假定滑动面是一个圆弧面，并认为条块间的作用力对土坡的整体稳定 性影响不大，故而忽略不计。或者说，假定条块两侧的作用力大小相等，方向相反且作用 于同一直线上。图7-8中取条块i进行分析，由于不考虑条块间的作用力，根据径向力的静 力平衡条件，有：N=Wcos 0 .(7-10)11 i根据滑动弧面上的极限平衡条件，有：T = TJF =(c. l.+N.tg 甲.)/ F(7-11)1 f1 s 111 Q 1 s式中：Tfi条块i在滑动面上的抗剪强度；气一滑动圆弧的稳定安全系数。另外，按照滑动土体的整体力矩平衡条件，外力对圆心力矩之和为零。在条块的三个(7-12)(7-13)作用力中，法向力N.通过圆心不产生力矩。重力叫产生的滑动力矩为: 1EWi d = EWi R ；in们滑动面上抗滑力产生的抗滑力第、：1ZTR :L " + N皿i . R iFs 滑动土体的整体力矩平衡，即 SM= 0 ,故有：EW.*d = ET.*R(7-14)将式(7-12)和式(7-13)代入式(7-14)，并进行简化，得：Z (c l + W cos 0 tg© )F =2 Wsin 0 (7-15)i i式(7-15)是最简单的条分法计算公式，因为它是由瑞典人费伦纽斯(W.Fellenius)等首 先提出的，所以称为瑞典条分法，又称为费伦纽斯条分法。从分析过程可以看出，瑞典条分法是忽略了土条块之间力的相互影响的一种简化计算 方法，它只满足于滑动土体整体的力矩平衡条件，却不满足土条块之间的静力平衡条件。 这是它区别于后面将要讲述的其它条分法的主要特点。由于该方法应用的时间很长，积累 了丰富的工程经验，一般得到的安全系数偏低，即误差偏于安全，所以目前仍然是工程上 常用的方法。三、毕肖甫条分法毕肖甫(A.N.Bishop)于1955年提出一个考虑条块间侧面力的土坡稳定性分析方法，称 为毕肖甫条分法。此法仍然是圆弧滑动条分法。在图7-9中，从圆弧滑动体内取出土条.进行分析。作用在条块.上的力，除了重力w外，滑动面上有切向力t和法向力N.，条块的侧面分别有法向力p广p.+1和切向力h、 H+1。假设土条处于静力平衡状态，根据竖向力的平衡条件，应有：'1+12 F广0(7-16)(7-11)W. +AH = N. cos 0. + T sin 0.N. cos 0. = W. +AH - T sin 0.根据满足土坡稳定安全系数Fs的极限平衡条件有：T =(弓.1+地既©.” Fs式中:图7-9毕肖甫法条块作用力分析将式(7-11)代入式(7-16)，整理后得：, c lW. +AHsin 0.N =7s = .cos 0. + 衅tm0s% = cos 0 + M0. fsrW +AH -c l 、sin 0Fs(7-17)(7-18)考虑整个滑动土体的整体力矩平衡条件，各个土条的作用力对圆心的力矩之和为零。这时条块之间的力P.和耳成对出现，大小相等，方向相反，相互抵消，对圆心不产生力矩。 滑动面上的正压力N.通过圆心，也不产生力矩。因此，只有重力W,和滑动面上的切向力 T.对圆心产生力矩。按式(7-14 )'SW d 二ZTRi ii将式(7-11)代入上式，得S WR sin 0. =S2(c./. + N.tg©.) R将式(7-17)的值代入上式，s简化后得Scb. + (W. +AH )tg©.m y (7-19)乙 W sin 0.这就是毕肖甫条分法计算土坡稳定安全系数Fs的一般公式。式中的乙凡=巧+1-巧，仍 然是未知量。如果不引进其它的简化假定，式(7-19)仍然不能求解。毕肖甫进一步假定A * =0，实际上也就是认为条块间只有水平作用力P.，而不存在切向作用力H。于是式(7-19) 进一步简化为：W sin 0(7-20)此式称为简化的毕肖甫公式。式中的参数m&包含有稳定安全系数F,。因此，不能直接 求出土坡的稳定安全系数Fs，而需要采用试算的办法，迭代求算Fs值。为了便于迭代计算， 已编制成mQ-0关系曲线，如图7-10。'图7-10 mQ曲线图试算时，可以先假定F, =1.0，由图7-10查出各个0 .所相应的m0值，并将其代入式（7-20） 中，求得边坡的稳定安全系数F。若F'与Fs之差大于规定的误差，用F'查m0，再次计 算出稳定安全系数F"，此如这样反复迭代计算，直至前后两次计算的稳定安全系数非常接s近，满足规定精度的要求为止。通常迭代总是收敛的，一般只要试算34次，就可以满足 迭代精度的要求。与瑞典条分法相比，简化的毕肖甫法是在不考虑条块间切向力的前提下，满足力的多 边形闭合条件，也就是说，隐含着条块间有水平力的作用，虽然在公式中水平作用力并未 出现。所以它的特点是：（1）满足整体力矩平衡条件；（2）满足各个条块力的多边形闭合条件， 但不满足条块的力矩平衡条件；（3）假设条块间作用力只有法向力没有切向力；（4）满足极限 平衡条件。由于考虑了条块间水平力的作用，得到的稳定安全系数较瑞典条分法略高一些。 很多工程计算表明，毕肖甫法与严格的极限平衡分析法，即满足全部静力平衡条件的方法（如 下述的简布法）相比，结果甚为接近。由于计算过程不很复杂，精度也比较高，所以，该方 法是目前工程中很常用的一种方法。四、普遍条分法（简布法，N.Janbu）普遍条分法的特点是假定条块间水平作用力的位置。在这一假定前提下，每个土条块 都满足全部的静力平衡条件和极限平衡条件，滑动土体的整体力矩平衡条件也自然得到满 足。而且，它适用于任何滑动面，而不必规定滑动面是一个圆弧面，所以称为普遍条分法。 它是由简布提出的，又常称为简布法。图7-11简布法条块作用力分析从图7-11 (a)滑动土体ABC中取任意条块i进行静力分析。作用在条块上的力及其作 用点见图7-11(b)所示。按照静力平衡条件I - 0,得：W. +AH = N. cos 0 iN cos0 = W +AHIF =0，得：X+ T sin 0.-T sin0：(7-16)AP = T cos 0 一 N sin 0(7-21)将式(7-16 )代入式(7-21)整理后得：'八 sin2 0 '八AP = T cos0 +- (W +AH.)tg0.''I 'cos七 J '''(7-22)根据极限平衡条件，考虑土坡稳定安全系数f sT. =+ N .tg.)s由式(7-16)得：(7-11)N=i1cos 0i(W. + AH. - T sin0.)(7-23)1FT =ii(7-24)代入式(7-11),整理后得;1c.i. + co； (W. +AH. tg.)i:tg6tg1 +1.Fs将式(7-24)代入式(7-22),得:1 sec20(7-25)P = 1 0 t、c.l. cos0. + (W, + AH,)tg0. - (W, + AH,)tg0.S 1 +.i-FS图7-12条块侧面法向力图7-12表示作用在土条条块侧面的法向力 P,显然有 P1=A P1, P2=P1+A P=L P1+ P2,依此类推，有：P = AP(7-26)ijj=1若全部土条条块的总数为，则有：P =£ AP = 0i=1(7-27)将式(7-25)代入式(7-27)，得:Z7Sez2c.l. cos0. + (W + AH.)tg。. £(W + AH.)tg0. =0Ftg0 - tgQ 1111 ii1 iiS + FS整理后得:Vsec20乙cl cos0 + (W +AH )tg© _- i i i i i i 1 + tg0 tg / Fs£ (W + AH )tg0i_i'ii i£ cb + (W +AH 硕i i ii i m0- (W +AH )sin 0iii(7-28)比较毕肖甫公式(7-19)和简布公式(7-28)，可以看出两者很相似，但分母有差别，毕肖 甫公式是根据滑动面为圆弧面，滑动土体满足整体力矩平衡条件推导出的。简布公式则是利用力的多边形闭合和极限平衡条件，最后从£AP = 0得出。显然这些条件适用于任何形ii=1式的滑动面而不仅仅局限于圆弧面，在式(7-28)中， H.仍然是待定的未知量。毕肖甫没有 解出 H，而让 H=Q，从而成为简化的毕肖甫公式而简布法则是利用条块的力矩平衡 条件，因而整个滑动土体的整体力矩平衡也自然得到满足。将作用在条块上的力对条块滑弧段中点O.取矩(图7-11(b)，并让EM o,= 0。重力W.和滑弧段上的力N.和T.均通过o,不产生力矩。条块间力的作用点位置已确定，故有:-（p +ap）h + Ah -1 1 1 1r 1+ P h - - AX jtgq = 011 1 2117略去高阶微量整理后得：H AX - P Ah -Aph = 0 i i i i i iH = P -Ah + AP(7-29)AX,i AX,AH, = H,+- - %(7-30)式(7-29)表示土条间切向力与法向力之间的关系。式中符号见图7-11。由公式(7-25)、(7-26)、(7-27)、(7-28)、(7-29)和(7-30)，利用迭代法可以求得普遍条分法的边坡稳定安全系数丹。 其步骤如下：'(1) 假定AH, = 0，利用式(7-28)，迭代求第一次近似 的边坡稳定安全系数Fs1。(2) 将Fs1和AH, = 0代入式(7-25)，求相应的AP (对每 一条块，从1到n)。(3) 用式(7-26) P = AP,求条块间的法向力(对每j=1一条块，从1到n)。(4) 将P.和乙P.代入式(7-29)和(7-30)，求条块间的切向作用力日,(对每一条块，从1到 研呸H,。(5) 将乙H.重新代入式(7-28)，迭代求新的稳定安全系数 Fs2。'如果Fs2"Fs1>A (为规定的计算精度)，重新按上述 步骤(2)(5)进行第二轮计算。如此反复进行，直至F v -Fs (k) s令 AH = 0F 式（7-28）AP 式（7-25） 11 FP：式(7-26)H 式（7-29）AH,式（7-30）F,式（7-28）End图7-13简布法计算程序流程(k-1)M为止。Fs(k)就是该假定滑动面的稳定安全系数。边坡真正的稳定安全系数还要计算很多滑动面，进行比较，找出最危险的滑动面，其边坡稳定安全系数才是真正的安全系数。这种计算工作量相当浩繁， 一般要在计算机上计算。用普遍条分法计算一个滑动面稳定安全系数的流程如图7-13。【例题7-1】 一简单的粘性土坡，高25m，坡比1：2，辗压土的重度y =20kN/m3，内摩 擦角4 =26.6° (相当于tg4 =0.5)，粘结力c =10kN/m2，滑动圆心0点如图7-14所示，试分别用瑞典条分法和简化毕肖甫法求该滑动圆弧的稳定安全系数，并对结果进行比较。27S图7-14 例题7-1图Z W sin9 = 3584kNii边坡稳定安全系数Fs分别计算各条块的重量叫，按照瑞典条分法和简化毕肖Z cl = 650kNi i解：为了使例题计算简单，只将滑动土体分成6个土条， 滑动面长度I，滑动面中心与过圆心铅垂线的圆心角e，然后 甫法进行稳定分析计算。1 .瑞典条分法瑞典条分法分项计算结果见例表7-1。ZW cos9 tg© = 4228kNi i i£ (Wcos 9 .tg。. + cL) 4228 + 650 广=Vil= 1.36乙 W sin93584ii2.简化毕肖甫法根据瑞典条分法得到计算结果Fs=1.36，由于毕肖甫法的稳定安全系数稍高于瑞典条分 法。设Fs1=1.55,按简化的毕肖甫条分法列表分项计算，结果如例表7-2。S£ Cg。. = 5417kNm9i条块编号9 i(° )W.1(kN)sin9icos9iW sin9 ii(kN)W cos9.(kN)'W cos 9 .tg。. '(kN)'l. 1 (m)c.l.(kN)-1-9.93412.5-0.1720.985-71.0406.32038.08000160001.00160080010.0100113.2923750.2300.9735462311115610.5105227.3726250.4600.88812072331116611.5115343.6021500.6900.7241484155777914.0140459.55487.50.8620.50742024712411.0110例题7-1瑞典条分法计算成果例表7-1例表7-2例题7-1毕肖普法分项计算成果编号cos9isin9isin 9 tg。iisin9 tg。 fiFsM9iW sin 9 iicb.W.tg。.1Icb + W tg。.m-10.985-0.172-0.086-0.0550.93-7180206.3307.801.000001.00010080090010.9730.2300.1150.0741.0475461001188123020.8880.4600.2300.1481.03612071001313136430.7240.6900.3450.2230.94714841001075124140.5070.8620.4310.2780.78542050243.8374.3安全系数£5.)5417£W sin9= 3 =1,51i i毕肖甫法稳定安全系数公式中的滑动力EW.sin9f与瑞典条分法相同。Fs1-Fs2 = 0.04,误差较大。按F = 1.51,s2£ 呐+ 忡.=5404.8进行第二次迭代计算，结果列于例表7-3中。 n "5.)=业= 1.507稳定安全系数乙 W sin93586i i气一。=。.。3,十分接近，因此，可以认为F = 1.51。例表7-3例题7-1毕肖普法第二次迭代计算成果编号cos0isin 0.sin 0 tg。.sin 0 tg。 iFF%W sin0.平W tg。IIc b + W tg。.m0-10.985-0.172-0.086-0.0570.928-7180206.3308.501.000.0001.00010080090010.9730.2300.1150.0761.04554610011881232.520.8880.4600.2300.1521.040120710013131358.630.7240.6900.3450.2280.952148410010751234.240.5070.8620.4310.2850.79242050243.8371计算结果表明，简化毕肖甫条分法的稳定安全系数较瑞典条分法高，约大0.15，与一 般结论相同。五、有限元法从瑞典条分法到普遍条分法的基本思路都是把滑动土体分成有限宽度的土条，把土条 当成刚体，根据滑动土体的静力平衡条件和极限平衡条件，求得滑动面上力的分布，从而 可以计算出边坡稳定安全系数。但是，因为土体是变形体，而并非是刚体，所以，引用 分析刚体的办法来分析变形体，并不满足变形协调条件，因而计算出的滑动面上的应力状 态不可能是真实的。有限元法就是把土坡当成变形体，按照土的变形特性，计算出土坡内 的应力分布，然后，再把圆弧滑动面的概念引入其中，验算滑动土体的整体抗滑稳定性。将土坡划分成许多单元体如图7-15所示。用有限元法可以计算出每个土单元的应力、 应变和每个结点的结点力和位移。这种计算目前已经成为土石坝应力变形分析的常用方法， 有各种现成的程序可供应用。图7-16表示的是一座土坝采用有限元法分析得到的竣工时坝 体的剪应变分布图，可以清楚看出坝坡在重力的作用下剪切变形的轨迹类似于滑弧面。图7-15 土坝的有限元网格和滑孤面图7-16某坝竣工后的剪应变分布（有限元法分析）土坡的应力计算出来以后，再引入圆弧滑动面的概念。图7-16中表示一个可能的圆弧 滑动面。把可能的圆弧滑动面划分成若干个小弧段，小弧段 I.上的应力用弧段中点的 应力代表，其值可以按照有限元法应力分析的结果，根据弧段中点所在的单元的应力确定， 表示为b，b，t 。如果小弧段 l.与水平线的倾角为。，则作用在弧段上的法向应力 和剪应力分别为"b . = 2（b . +b .） -2（b . b .）cos20. +t . sin20.（7-31）(7-32)t . = t . cos 20. - ； 9 . -。 )sin 20.根据莫尔一库伦强度理论，该点土的抗剪强度为T f. W +b n tg©.将滑动面上所有小弧段的剪应力和抗剪强度分别求出来以后，再累加求得沿着滑动面 总的剪切力E T /和抗剪力E T 。因此，边坡稳定安全系数为"'Fs =乙 <c. + a 试 tg©.) Mj 乙t . AZ.(7-33)i=1i=1很显然，有限元分析方法的优点是把边坡稳定分析与坝体的应力和变形分析结合起来。 这时，滑动土体自然满足静力平衡条件而不必如条分法那样引入人为的假定。但是，当边 坡接近失稳时，滑裂面通过的大部分土单元处于临近破坏状态，这时，用有限元法分析边 坡内的应力和变形所需要的土的基本特性，如变形特性，强度特性等均变得十分复杂，因 此，要提出一种能反映土体实际受力状况的计算模型是很不容易的。如果说在边坡稳定性 分析中极限平衡分析法是当前工程上主要应用的方法，那么，有限元方法则是一种潜在的 具有很大发展前景的方法。六、最危险滑裂面的确定方法和容许安全系数(一)最危险滑裂面的位置以上介绍的是计算某个位置已经确定的滑动面稳定安全系数的几种方法。这一稳定安 全系数并不代表边坡的真正稳定性，因为边坡的滑动面是任意选取的。假设边坡的一个滑 动面，就可计算其相应的安全系数。真正代表边坡稳定程度的稳定安全系数应该是稳定安 全系数中的最小值。相应于边坡最小的稳定安全系数的滑动面称为最危险滑动面，它才是 土坡真正的滑动面。确定土坡最危险滑动面圆心的位置和半径大小是稳定分析中最繁琐、工作量最大的工 作。需要通过多次的计算才能完成。这方面费伦纽斯(W.Fellenius)提出的经验方法，对于较 快地确定土坡最危险的滑动面很有帮助。费伦纽斯认为，对于均匀粘性土坡，其最危险的滑动面一般通过坡趾。在©=0法的边 坡稳定分析中，最危险滑弧圆心的位置可以由图7-17( a)中pi和p 2夹角的交点确定。、 P2的值与坡角a大小的关系，可由表7-1查用。11对于©0的土坡，最危险滑动面的圆心位置如图7-17( b)所示。首先按图7-17( b)中 所示的方法确定DE线。自£点向DE延线上取圆心O'。之，通过坡趾A分别作圆弧， AC】、AC?、，并求出相应的边坡稳定安全系数Fs1、Fs2。表7-1各种坡角的八、八值坡角a坡度1 : m 260°1 : 0.5829°40°45°1 : 1.028°37°33° 41，1 : 1.526°35°26° 34，1 : 2.025°35°18° 26，1 : 3.026°35°14° 02，1 : 4.025°36°11° 19，1 : 5.025°39°然后，再用适当的比例尺标在相应的圆心点上，并且连接成安全系数Fs随圆心位置的变化 曲线。曲线的最低点即为圆心在DE线上时安全系数的最小值。但是真正的最危险滑弧圆心 并不一定在DE线上。通过这个最低点，引DE的垂直线FG。在FG线上，在DE延线的最 小值前后再定几个圆心O" o,2，用类似步骤确定FG线上对应于最小安全系数的圆心， 这个圆心。才被认为是通过坡趾滑出时的最危险滑动圆弧的中心。当地基土层性质比填土软弱，或者坝坡不是单一的土坡，或者坝体填土种类不同、强 度互异时，最危险的滑动面就不一定从坡趾滑出。这时寻找最危险滑动面位置就更为繁琐。 实际上，对于非均质的、边界条件较为复杂的土坡，用上述方法寻找最危险滑动面的位置将是十分困难的。随着计算机技术的发展和普及，目前可以采用最优化方法，通过随机搜 索，寻找最危险的滑动面的位置。国内已有这方面的程序可供使用。表7-2碾压式土石坝坝坡容许稳定安全系数运用条件工程等级IIIIIIIV、V正常运用条件1.301.251.201.15非常运用条件I1.201.151.101.05非常运用条件II1.101.051.051.00注：正常运用条件系指:(1) 水库水位处于正常高水位(或设计洪水位)与死水位之间的各种水位下的稳定渗流期；(2) 水库水位在上述范围内的经常性正常降落；(3) 抽水蓄能电站的水库水位的经常性变化和降落。非常运用条件I系指：(1) 施工期；(2) 校核洪水位下有可能形成稳定渗流的情况；(3) 水库水位的非常降落，如自校核洪水位降落、降落至死水位以下、大流量快速泄空等；(4) 正常运用条件遭遇地震。非常运用条件II系指以上非常运用条件(1)(3)再遭遇地震的情况。(二)边坡容许安全系数在土坡稳定的分析中，从土体材料的强度指标到计算方法，很多因素都无法准确确定。 因此，如果计算得到的土坡稳定安全系数等于1或稍大于1,并不表示边坡的稳定性能得到 可靠的保证。安全系数必须满足一个最起码的要求，称为容许安全系数。容许安全系数值 是以过去的工程经验为依据并以各种规范的形式确定。因此采用不同的抗剪强度试验方法 和不同的稳定分析方法所得到的安全系数差别甚大，所以在应用规范所给定的土坡稳定容 许安全系数时，一定要注意它所规定的试验方法和计算方法。表7-2为1984年水电部颁布的碾压土石坝设计规范(SDJ21884)中的边坡容许安 全系数表。表中除注明者外，均适用于瑞典圆弧法。对I、I级的中、高土石坝以及一些 复杂的情况，应同时采用毕肖甫法或其它更严格的方法(如普遍条分法等)进行计算。此时安 全系数的容许值，应比表中所规定的值略微提高10%左右。对于I级土石坝在正常使用条 件下，安全系数不得小于1.5。第四节天然土坡的稳定问题天然土体由于形成的自然环境、沉积时间以及应力历史等因素不同，性质比人工填土 要复杂得多，边坡稳定分析仍然可按上述方法进行，但在强度指标的选择上要更为慎重。一、裂隙硬粘土的边坡稳定性硬粘土通常为超固结土，其应力一应变关 i °系曲线属应变软化型曲线，如图7-18所示。这/X类土如果也按一般的天然土坡稳定分析办法，/ 认为剪切过程中密度不变，故宜采用