岩土边坡稳定性分析新方法与工程应用课件

1、岩土边坡稳定性分析新方法与工程应用1岩土边坡稳定性分析新方法岩土边坡稳定性分析新方法与工程应用与工程应用解放军理工大学工程兵工程学院解放军理工大学工程兵工程学院中中 铁铁 西西 北北 科科 学学 研研 究究 院院岩土边坡稳定性分析新方法与工程应用2第一章第一章 绪论绪论 第二章第二章 边坡稳定性极限平衡法显式解答边坡稳定性极限平衡法显式解答 第三章第三章 边坡稳定性极限平衡法统一计算框架边坡稳定性极限平衡法统一计算框架 第四章第四章 锚固边坡稳定性计算分析方法锚固边坡稳定性计算分析方法 第五章第五章 morgenstern-pricemorgenstern-price法的改进法的改进 第六章第

2、六章 严格严格janbujanbu法的改进法的改进 第七章第七章 边坡临界滑动场的提出与数值模拟方法边坡临界滑动场的提出与数值模拟方法第八章第八章 基于严格条分法的边坡临界滑动场基于严格条分法的边坡临界滑动场 第九章第九章 边坡临界滑动场法的工程应用边坡临界滑动场法的工程应用 第十章第十章 应用临界滑动场法计算土压力应用临界滑动场法计算土压力第十一章第十一章 应用临界滑动场法计算地基承载力应用临界滑动场法计算地基承载力 第十二章第十二章 总结总结岩土边坡稳定性分析新方法与工程应用3边坡稳定性极限平衡法显式解答边坡稳定性极限平衡法显式解答 岩土边坡稳定性分析新方法与工程应用4传统方法传统方法：通

3、过对条块间作用力方式进行假设，使：通过对条块间作用力方式进行假设，使滑体满足部分或全部平衡条件。滑体满足部分或全部平衡条件。 除瑞典法外，安全系数是隐含于平衡方程或方除瑞典法外，安全系数是隐含于平衡方程或方程组，需迭代求解。如需严格满足平衡条件，安全程组，需迭代求解。如需严格满足平衡条件，安全系数求解过程非常复杂且可能出现不收敛。系数求解过程非常复杂且可能出现不收敛。显式方法显式方法：通过对滑面正应力分布形状进行假设，：通过对滑面正应力分布形状进行假设，使滑体满足全部平衡条件。使滑体满足全部平衡条件。 安全系数为显式表达式，求解过程简单，不用安全系数为显式表达式，求解过程简单，不用划分条块，不

4、存在不收敛问题。划分条块，不存在不收敛问题。岩土边坡稳定性分析新方法与工程应用521baab (x)u(x) e(x)y=z(x) t(x) wkcw (xc,yc)s(x)y=g(x)yxo(x)aba1a2(a)滑体受力情况及滑面正应力分布假设滑体受力情况及滑面正应力分布假设abaa311abaa322岩土边坡稳定性分析新方法与工程应用6滑面正应力分布假设滑面正应力分布假设假设滑面两端正应力由简单的微分条块平衡条件确定，滑面中间两正应力值为两个待定参数，滑面正应力分布为3次样条函数。thrust forceu(x)(x)(x)(x)(x)kw(x)w(x)1 xxxfxxxxxcxxufx

5、xkxxwxcsintan1cossintan1cossincos00滑面两端正应力计算滑面两端正应力计算滑面应力分布滑面应力分布 xxxx32211 211121aabaaaaxbxaxx 122212aabaaaaxbxaxx 212121213abababaxaxaxaaaabaaxaxbxxba岩土边坡稳定性分析新方法与工程应用7极限平衡方程组极限平衡方程组 babacwkdxxdxxsx (7a) wdxxsxdxxbaba (7b) 0dxxxxsxxxsyxxsxbacc (7c) xcxxuxfxstan1平衡条件：水平力平衡、竖直力平衡、力矩平衡平衡条件：水平力平衡、竖直力平

6、衡、力矩平衡破坏准则：破坏准则：滑面应力分布滑面应力分布 xxxx32211岩土边坡稳定性分析新方法与工程应用8 dxxxfxsdxxxuxcfwkdxxxfxsdxxxfxsbasbascbasbas32211tan1tan1 tan1tan1 (12a) basbasbasbasdxxxfxsdxxxuxcxsfwdxxxfxsdxxxfxs32211tan1tan1tan1tan1 (12b) dxxxxsyxsxdxxxxsxsyxcxxuxxxfbaccbaccstan32211 (12c)含含3个末知量（个末知量（1 , 2 and fs）的）的3个平衡方程：个平衡方程：岩土边坡稳

7、定性分析新方法与工程应用9安全系数求解安全系数求解33222111111afaafaafasss (13a)33222111111bfbbfbbfbsss (13b)3221132211eeedddfs (13c)简化的平衡方程组：简化的平衡方程组：2210221011111gfgfgtftftssss2210221021111gfgfgsfsfsssss232221213121103020123222121312110302011111gesetefgesetefgesetegdsdtdfgdsdtdfgdsdtdfsssss岩土边坡稳定性分析新方法与工程应用10001223tftftfss

8、s033232qtfptfss33233223223223pqqpqqtfs根据卡尔丹公式，上式根有4种组合情况：（1）1个实根，两个共轭复根；（2）3个重零根；（3）3个实根中，有两个相等；（4）3个不相等的实根。对实际边坡稳定性问题，不存在2个或3个不等的安全系数同时使边坡达到极限平衡状态，因此式p1-22只有1个实根，其余为无意义的复根，该实根为安全系数显解：安全系数显解：安全系数安全系数3次代数方程：次代数方程：岩土边坡稳定性分析新方法与工程应用11算例算例1：与理论解比较：与理论解比较the theory of plasticity : (x)the present method :

9、 (x)line of thrustkpa5040302010 045111.4 kpa10 m无粘性土边坡，坡顶受均布荷载作用，其极限荷载有理论解，sokolovskii解。c=10 kpa and =30理论极限荷载111.4 kpa理论安全系数fs=1.0实际计算安全系数fs=1.032结论：结论：尽管本法对应的滑面正应力分布与理论解有一定的尽管本法对应的滑面正应力分布与理论解有一定的差别，但计算的安全系数或坡面极限荷载与理论解非常接差别，但计算的安全系数或坡面极限荷载与理论解非常接近，据此可认为，建议假设的滑面正应力分布形状在边坡近，据此可认为，建议假设的滑面正应力分布形状在边坡工程应

10、用中是可以接受的。工程应用中是可以接受的。岩土边坡稳定性分析新方法与工程应用12算例算例2：与：与spencer法比较法比较surface 4surface 3surface 2surface 1-6.0 m-9.0 m 9.0 m 3.0 m(4)(2)(1)0.015.0 m12(3)table 2. comparison of factors of safety computedwater pressureconditionslip surfacepresentspencermorgensternand pricesimplifiedbishopordinary11.5601.5591.

11、5591.5341.49621.5841.6161.62831.1661.2111.1971.0790.922wet slope41.1091.1501.14112.0352.0352.0352.0111.93522.0492.0872.10431.7441.8361.8231.4291.229dry slope41.7091.7721.765note: the interslice function used in morgenstern-price method is half sine.layerc(1)18.8 kn/m320.0 kpa18.0(2)18.5 kn/m340.0 kp

12、a22.0(3)18.4 kn/m325.0 kpa26.0(4)18.0 kn/m310.0 kpa12.0结论结论：本例中，本文方：本例中，本文方法与法与spencer法计算安法计算安全系数最大误差不超过全系数最大误差不超过5%；spencer法不光滑法不光滑的滑面正应力分布可用的滑面正应力分布可用建议的光滑分布形式代建议的光滑分布形式代替；内力分布的尚在合替；内力分布的尚在合理范围内。理范围内。岩土边坡稳定性分析新方法与工程应用13 (x)400300200100 0150100 50 0t(x)e(x)present methodspencer methodpresent method

13、spencer methodsimplified bishopmethodu(x)present methodspencer method(a) lines of thrust forces(b) distribution of total normal stresses along the slip surface(c) magnitude of internal forceskn/mkpa岩土边坡稳定性分析新方法与工程应用14400300200100 0150100 50 0t(x)e(x)present methodspencer methodpresent methodspencer

14、methodu(x) (x)kn/mkpaspencer methodpresent method(a) lines of thrust forces(b) distribution of total normal stresses along the slip surface(c) magnitude of internal forces岩土边坡稳定性分析新方法与工程应用15t(x)e(x)u(x) (x)200016001200 800 400 0300200 100 0kn/mpresent methodspencer method(a) lines of thrust forces(b

15、) distribution of total normal stresses along the slip surface(c) magnitude of internal forcespresent methodspencer methodpresent methodspencer methodsimplified bishopmethodkpa岩土边坡稳定性分析新方法与工程应用16u(x) (x)2200200016001200 800 400 0400300 200 100 0t(x)e(x)present methodspencer methodkn/mkpapresent meth

16、odspencer method(a) lines of thrust forces(b) distribution of total normal stresses along the slip surface(c) magnitude of internal forcespresent methodspencer methodsimplified bishopmethod岩土边坡稳定性分析新方法与工程应用17边坡稳定性极限平衡法统一计算框架边坡稳定性极限平衡法统一计算框架岩土边坡稳定性分析新方法与工程应用18理论与工程背景理论与工程背景：现今条分法有十几种之多，每种方法都曾：现今条分法有十

17、几种之多，每种方法都曾得到不同程度地应用，已积累了大量的使用经验。许多国家得到不同程度地应用，已积累了大量的使用经验。许多国家的规范规定，对同一边坡特别是土石坝问题，宜采用多种方的规范规定，对同一边坡特别是土石坝问题，宜采用多种方法同时计算，比较结果，根据经验判断其合理性。然后各种法同时计算，比较结果，根据经验判断其合理性。然后各种方法是不同的提出者根据当时的理解形成的计算格式，后人方法是不同的提出者根据当时的理解形成的计算格式，后人大都沿用原始格式计算。而这些原始计算格式，各种算法的大都沿用原始格式计算。而这些原始计算格式，各种算法的不统一，也不便于人们理解各自方法的优缺点。用上述显示不统一

18、，也不便于人们理解各自方法的优缺点。用上述显示解格式将所有条分法统一起来，使计算原理更为清晰、计算解格式将所有条分法统一起来，使计算原理更为清晰、计算过程更为便捷、一些常见的数值困难得到克服。过程更为便捷、一些常见的数值困难得到克服。统一格式要点统一格式要点：将现有：将现有1212种条分法分成种条分法分成4 4种平衡条件组合，每种平衡条件组合，每种组合推导出安全系数显示表达式（基于滑面正应力修正模种组合推导出安全系数显示表达式（基于滑面正应力修正模式），式），1212种条分法对应滑面正应力可用通式表达，通过自动种条分法对应滑面正应力可用通式表达，通过自动迭代求解与传统方法意义一致的安全系数。迭

19、代求解与传统方法意义一致的安全系数。岩土边坡稳定性分析新方法与工程应用19u(x)x/cosx/cose-e/2e+e/2t+t/2t-t/2hthkcwxwxqxxqyxxba0(x)(x)(x)u(x)0(x)(x)g.w.lyt(x)s(x)g(x)taeatbqx(x)qy(x)(xc , yc) x yoeb(a)(b)(c)baabxceedxqwksbaabyttdxqws tacacatbcbcbbacxcccycyyexatyyexbtdxgyqgsywkxxqwssys5 . 05 . 0平衡方程：平衡方程：cufs1破坏准则：破坏准则：简化平衡方程：简化平衡方程：baba

20、sxsdxcuffdxfs11babasysdxcusffdxfs111bacbabasdxrmdxrcudxrf岩土边坡稳定性分析新方法与工程应用20现有条分法考虑现有条分法考虑4 4种平衡条件组合种平衡条件组合，即：，即：（1 1）考虑所有平衡条件即：水平、垂直力及力矩平衡）考虑所有平衡条件即：水平、垂直力及力矩平衡 ( (简称简称hvmhvm组合组合) )。spencerspencer法法、morgernstern-pricemorgernstern-price法法、sarma2sarma2法法、sarma3sarma3法法、correiacorreia法法。（2 2）考虑垂直方向力的平

21、衡和对选定的求矩中心的力矩平衡）考虑垂直方向力的平衡和对选定的求矩中心的力矩平衡（简称（简称vmvm组合）。组合）。简化简化bishopbishop法法。（3 3）考虑水平、垂直力的平衡（简称）考虑水平、垂直力的平衡（简称hvhv组合）。组合）。简化简化janbujanbu法法、美国陆军工程师团法美国陆军工程师团法、lowe-karafiathlowe-karafiath法法、sarma1sarma1法法。严格严格janbujanbu法法自动考虑了力矩平衡，求解过程中只用了两个力的平自动考虑了力矩平衡，求解过程中只用了两个力的平衡，因此归属衡，因此归属hvhv组合。组合。（4 4）仅对选定的求

22、矩中心的力矩平衡（简称）仅对选定的求矩中心的力矩平衡（简称m m组合），此法为组合），此法为瑞典法瑞典法。岩土边坡稳定性分析新方法与工程应用21滑面正应力修正滑面正应力修正 xxx0 babxx1 abaxx2 xxx2211211112vm、hv组合：m组合：岩土边坡稳定性分析新方法与工程应用22安全系数求解安全系数求解33222111111afaafaafasss (14a) 33222111111bfbbfbbfbsss (14b) 3221132211eeedddfs (14c)babasxbassdxcuffdxfsdxfs111202101 (13a) babasybassdxcu

23、sffdxfsdxfs11111202101 (13b) babacbababasdxrdxrmdxrcudxrdxrf202101202101 (13c)岩土边坡稳定性分析新方法与工程应用233323322)hvm(3223223pqqpqqtfs1211)hm(22qppfs2222)hv(22qppfs321321)m(eeedddfs安全系数表达式：安全系数表达式：岩土边坡稳定性分析新方法与工程应用24各种条分法滑面正应力初始值通式各种条分法滑面正应力初始值通式 xxexxt21条分法条间力假设通式条分法条间力假设通式dxdecuqwksxc0000dxdtscuqwsy0000101

24、0210010001ssecuwkqscuqwcxy微分条块力平衡条件微分条块力平衡条件滑面正应力初始值通式滑面正应力初始值通式考虑过破坏准则考虑过破坏准则cufs001岩土边坡稳定性分析新方法与工程应用25条分法总结与有关参数计算 表 1 平衡条件编号方法假设1(x)2(x)1(x)2(x)fs垂直力水平力力矩1瑞典法 eq.(32) eq.(33a)eq.(33b)eq.(27)c2简化 bishop 法 eq.(34) eq.(35a)eq.(35b)eq.(23)sc3简化 janbu 法 eq.(34) eq.(35a)eq.(35b)eq.(26)ss4工程师团法 eq.(36)

25、eq.(37a)eq.(37b)eq.(26)ss5lowe & karafiath eq.(38) eq.(39a)eq.(39b)eq.(26)ss6sarma (i) eq.(40) eq.(41a)eq.(41b)eq.(26)ss7spencer eq.(42) eq.(43a)eq.(43b)eq.(20)sss8morgenstern & price eq.(44) eq.(45a)eq.(45b)eq.(20)sss9sarma (ii) eq.(46) eq.(47a)eq.(47b)eq.(20)sss10sarma (iii) eq.(48) eq.(49

26、a)eq.(49b)eq.(20)sss11correria eq.(50) eq.(51a)eq.(51b)eq.(20)sss12严格 janbu 法 eq.(52) eq.(53a)eq.(53b)eq.(26)ssas注: c = 考虑, s =满足, as = 自动满足；表中所列公式在 part 2 部分。岩土边坡稳定性分析新方法与工程应用26(1) ordinary method (fellenius, 1936) xexsxt and 0dxdt (32) xsx1, 02 x (33a) 01 x , 02 x (33b) (2) simplified bishop metho

27、d (bishop, 1955) 0 xt (34) 01 x , 02 x (35a) 01 x , 02 x (35b) (3) simplified janbu method (janbu, et.al, 1956) 0 xt (34) 01 x , 02 x (35a) 01 x , 02 x (35b)岩土边坡稳定性分析新方法与工程应用27 (4) corps of engineers method (u.s. army corps of engineers, 1970) xexgxt (36) xgx1 , 02 x (37a) xxx11 , 02 x (37b)(5) low-

28、karafiath method (low & karafiath, 1960) xexsxgxt21 (38) xsxgx211, 02 x (39a) xxx11 , 02 x (39b)(6) sarma method (i)(sarma, 1979) xhxcxxpxextwv0v0 (40) xxv01 , xhxcxxpxwv0v02 (41a) xxx11 , xxx22 (41b)岩土边坡稳定性分析新方法与工程应用28 (7) spencer method (spencer, 1967; 1973) xext (42) x1 , 02 x (43a) 01 x , 02

29、 x (43b) (8) morgenstern-price method (morgenstern & spencer, 1965) xexfxt (44) xfx1 , 02 x (45a) xxx11 , 02 x (45b)(9) sarma method (ii) (sarma, 1973) xhxcxxpxextwv0v0 (46) xxv01 , xhxcxxpxwv0v02 (47a) xxx11 , xxx22 (47b)岩土边坡稳定性分析新方法与工程应用29 (10) sarma method (iii) (sarma, 1973) xhxcxxpxexfxtwv0

30、v0 (48) xxfxv01 , xhxcxxpxfxwv0v02 (49a) xxx11 , xxx22 (49b) (11) correia method (correia, 1989) xfxt (50) 01 x ; xfx2 (51a) 01 x ; xxx22 (51b)(12) rigorous janbu method (janbu, 1954; 1973) xhxwkxhqxehyxextcxtt5 . 0 (52) tyx1 xhxwkxhqxehxcxt5 . 02 (53a) tyx 1 , xxx22 (53b)岩土边坡稳定性分析新方法与工程应用30严格严格janb

31、u法等数值导数平滑处理法等数值导数平滑处理1=-22y2y0y1 x1 x0 x2/2/22=-11y2y0y1y2y0y1x1 x0 x2=bx1=a x0 x2120yyxy10/ab 5 . 015 . 02 202111020yyyyxy x0 位于区间中部 x0 接近区间端部岩土边坡稳定性分析新方法与工程应用31内力计算与检验内力计算与检验 xaxcadqwksexe xayadqwstxt xaaxcxaataaxcytedqwksxatyedgqgswkxqwsssxy5 . 0 xtxhxcxxpxexfwsvvv条间水平力：条间水平力：条间竖向力：条间竖向力：推力线位置：推力

32、线位置：条间局部安全系数：条间局部安全系数：岩土边坡稳定性分析新方法与工程应用32严格条分法待定参数严格条分法待定参数 计计算算mmxext for method 7 (57a) mmmxfxext for method 8 (57b) mmmmwmsmxhxcxxpxefxtvv for method 9 (57c) mmmmmwmsmxfxhxcxxpxefxtvv for method 10 (57d)mmxfxt for method 11 (57e)in which 2baxm.对于不严格条分法，可将其对应的滑面正应力分布对于不严格条分法，可将其对应的滑面正应力分布进行一次修正，得到

33、严格意义上的安全系数进行一次修正，得到严格意义上的安全系数岩土边坡稳定性分析新方法与工程应用33 是 否 fs0 , 0 0(x) a1, a1, a2, a2, a3, a3; b1, b1, b2, b2, b3, b3; d1, d2, d3 ; e1, e2, e3 t0, t1, t2 ; s0, s1, s2 ; g0, g1, g2 t0, t1, t2 p1, q1 p2, q2 p, q fs(hvm) fs(vm) fs(hv) fs(m) e(x), t(x) fs ? fs(hvm) 开始 结束 岩土边坡稳定性分析新方法与工程应用34算例与比较算例与比较circular

34、 slip surfacegeneral slip surface 1 230.00.0layerc18.2 kn/m320.0 kpa32.018.0 kn/m325.0 kpa30.018.5 kn/m340.0 kpa18.018.8 kn/m340.0 kpa28.0岩土边坡稳定性分析新方法与工程应用35figure 3. comparison of factors of safety (circular slip surface)11.11.21.31.4ordinarysimplified bishopsimplified janbucorps of engineerslowe &

35、amp; karafiathsarma ispencermorgenstern & pricesarma iisarma iiicorreriarigorous janbufactor of safetyconventionalrigorous岩土边坡稳定性分析新方法与工程应用36figure 4. comparison of factors of safety( general slip surface)0. 911. 11. 21. 3ordinarysimplified bishopsimplified janbucorps of engineerslowe & kara

36、fiathsarma i spencermorgenstern & pricesarma iisarma iiicorreriarigorous janbufactor of safetyconventionalrigorous岩土边坡稳定性分析新方法与工程应用37vteev(r)pwt(r)e(r)v(r)pwt(r)e(r)u(r)u(r)v(r)pwt(r)e(r)u(r) kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200k

37、n/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0v(r)pwt(r)e(r)u(r) yt(r) yt(r) yt(r) yt(r)method 3:simplified janbumethod 4:corps of engineersmet

38、hod 2:simplified bishopmethod 1:ordinarytevvte(r) kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0teevpwvpwtuuv(r)pwt(

39、r)e(r)u(r) kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0v(r)pwt(r)eu(r) yt yt yt(r) yt(r)method 7:spencermethod 8:morgenstern & pricemethod 6:sarma (i)method 5:lowe & karafiath岩土边坡稳定性分析新方法与工程应用38t kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0

40、 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 ytvteev(r)pwt(r)e(r)vpwtu(r)uvpwteu kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0vpweu yt(r) yt yt ytm

41、ethod 11:correriamethod 12:rigorous janbumethod 10:sarma (iii)method 9:sarma (ii)vteev(r)pwt(r)e(r)v(r)pwt(r)e(r)u(r)u(r)v(r)pwt(r)e(r)u(r) kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 80

42、0 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0v(r)pwt(r)e(r)u(r) yt(r) yt(r) yt(r) yt(r)method 3:simplified janbumethod 4:corps of engineersmethod 2:simplified bishopmethod 1:ordinary岩土边坡稳定性分析新方法与工程应

43、用39t kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 ytvteev(r)pwt(r)e(r)vpwtu(r)uvpwteu kpa 800 600 400 200 0 -200kn

44、/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0vpweu yt(r) yt yt ytmethod 11:correriamethod 12:rigorous janbumethod 10:sarma (iii)method 9:sarma (ii)tevvte(r) kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.

45、0 0 -1.0 -2.0 kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0teevpwvpwtuuv(r)pwt(r)e(r)u(r) kpa 800 600 400 200 0 -200kn/m 8000 6000 4000 20000-2000 2.0 1.0 0 -1.0 -2.0v(r)pwt(r)eu(r) yt yt yt(r) yt(r)method 7:spencermethod 8:morgenstern & pricemethod 6:sarma (i)meth

46、od 5:lowe & karafiath岩土边坡稳定性分析新方法与工程应用40锚固边坡稳定性计算分析方法锚固边坡稳定性计算分析方法岩土边坡稳定性分析新方法与工程应用41传统方法传统方法：将锚固力作为边坡面外力直接计入所处条块，：将锚固力作为边坡面外力直接计入所处条块，再用传统条分法直接计算安全系数。再用传统条分法直接计算安全系数。缺点是引起的滑面正应力分布在个别点有突变，与实际严缺点是引起的滑面正应力分布在个别点有突变，与实际严重不符，内力也更加不合理，有时安全系数计算出现不收重不符，内力也更加不合理，有时安全系数计算出现不收敛。敛。本方法本方法：利用无限楔的弹性应力解答计算锚固力引

47、起滑面：利用无限楔的弹性应力解答计算锚固力引起滑面正应力，再与无锚固力作用下滑面正应力叠加，经过调整正应力，再与无锚固力作用下滑面正应力叠加，经过调整修正后使边坡满足极限平衡条件，直接得到已知锚固力的修正后使边坡满足极限平衡条件，直接得到已知锚固力的安全系数，或给定的安全系数所需的锚固力。安全系数，或给定的安全系数所需的锚固力。岩土边坡稳定性分析新方法与工程应用42te (xp , yp)r(b)(a)ppr 1 2ippp0kcbapp1pp2ox(xc , yc)y岩土边坡稳定性分析新方法与工程应用432121212121212121sin2sin2sinsin2cos2cos2iirpr

48、2sinrprpbu锚固力引起的滑面正应力分布锚固力引起的滑面正应力分布径向应力：径向应力：滑面正应力：滑面正应力：孔隙水压力：孔隙水压力： ppxx02211111pssbfcuf破坏准则：破坏准则： abaxxbabxx21;应力修正：应力修正：岩土边坡稳定性分析新方法与工程应用44简化式：简化式： baxpcpdxwks (13a) bayppdxws (13b) cpypcpxpbaccccxxpyypdxgsywkxxwssys5 . 05 . 0 (13c)baspbabasxpxsdxbfdxcufpfdxfs111 (15a)babaspbasypysdxbsfdxcusfpf

49、dxfs1111 (15b)baspbasppcbasdxbrfdxcurfmmdxfrr111 (15c)极限平衡方程极限平衡方程岩土边坡稳定性分析新方法与工程应用45安全系数解答安全系数解答pp01bcupc令0根据无锚作用下边坡稳定计算结果确定，建议用spencer法bacsxpxbabassdxfpfdxfsdxfs1 112211 (18a) bacsypybabassdxsfpfdxfsdxfs1 11112211 (18b)babappcbabacbasmmdxrdxrdxrdxrdxrf22112211 (18c)33233223223223pqqpqqtfs岩土边坡稳定性分析

50、新方法与工程应用46所需锚固力解答所需锚固力解答 sxfs1 ; syfs11 ; srfrr1 (27a)cufsu1 ; 11bfsb (27b)令bauxbapxbbappxbappxdxfpdxdxdx111022011 (29a)bauybapybbappybappydxsfpdxsdxdx111022011 (29b) babababprprpbabaucbarrpdxrmdxdxdxrmdxdx2211202101 (29c)33233223223223pqqpqqtp岩土边坡稳定性分析新方法与工程应用47layer =18 kn/m3c=15. kpa=28layer =18

51、kn/m3c=15. kpa=42-113 kn/m307 kn/mt(p)t(c)e(0)t(0)e(c)e(p)101 kpa390 kpa(0)(c)(p)lt(c)lt(0)lt(p)(b)p300.0 m15.0 m(a)(c)p=0. fs=0.998p=300 kn/m fs=1.287(fs=1.357, conventional)fs=1.5 p=485 kn/m算例算例 1岩土边坡稳定性分析新方法与工程应用48算例算例 2table 1. values of factors of safety and load factorsb=0.0b=0.25b=0.5b=0.75b=

52、1.0no. p=0 p=1.fs=1.2 p=1.fs=1.2 p=1.fs=1.2 p=1.fs=1.2 p=1.fs=1.2fs0fs pfs pfs pfs pfs p11.1731.3370.1761.3070.2151.2770.2751.2470.3831.2170.63121.1471.3140.3361.2840.4081.2530.5191.2230.7151.1931.14631.0931.2650.6441.2340.7771.2030.9791.1711.3231.1412.04241.0811.2520.7151.2210.8611.1901.0811.1601.4

53、541.1292.22351.0711.2460.7551.2150.9061.1831.1331.1521.5121.1212.27661.0581.2410.7951.2080.9541.1751.1931.1421.5941.1102.40171.0261.2330.8571.1991.0081.1651.2231.1311.5561.0972.13981.0311.2640.7481.2290.8671.1941.0311.1591.2711.1251.65791.0591.3240.5631.2890.6431.2530.7511.2170.9021.1821.130101.0551

54、.3740.4911.3310.5621.2880.6751.2450.7911.2010.994110.9961.6260.3511.4230.5061.2230.9071.0294.3070.844-1.241 p1 = p2 = p3 = 1000 kn/m 1 2 3 4 5 6 7 8 9 10 11 layer =18 kn/m3 c=10. kpa =33 pp2 pp1 hard stratum 45 pp3 30 40.0 m 20.0 m 0.0 m layer =19 kn/m3 c=0. kpa =20 60 岩土边坡稳定性分析新方法与工程应用49morgenstern

55、-pricemorgenstern-price法的改进法的改进岩土边坡稳定性分析新方法与工程应用50morgenstern-pricemorgenstern-price法是国际公认的性能最好的能满足所有法是国际公认的性能最好的能满足所有平衡条件的严格条分法，是国际通用首选的计算任意滑面安平衡条件的严格条分法，是国际通用首选的计算任意滑面安全系数的极限平衡法。全系数的极限平衡法。morgenstern-pricemorgenstern-price法基本点是，假设条间正切力比法基本点是，假设条间正切力比x x/ /e e在滑在滑体内分布用一函数体内分布用一函数 f f( (x x) ) 表示（表示

56、（ f f( (x x) )为条间力函数），为条间力函数）， 为为比例系数比例系数。应用两个方向的力平衡条件及力矩平衡条件解。应用两个方向的力平衡条件及力矩平衡条件解出出 及安全系数及安全系数f fs s。由于计算过程相当繁琐，在我国始终未得到广泛应用，不能由于计算过程相当繁琐，在我国始终未得到广泛应用，不能与国际学术界接轨。我国规范规定用剩余推力法，实际上是与国际学术界接轨。我国规范规定用剩余推力法，实际上是morgenstern-pricemorgenstern-price法的特例，同时也严格法的特例，同时也严格janbujanbu法的特例。法的特例。岩土边坡稳定性分析新方法与工程应用51

57、 slip surface phreatic surface line of thrust (n) (i) (1) (2) surface loads (a) pi fiei ei kcwi (b) hi hi/2 zi zi-1 ui =uibiseci i sk= (nitani+cbiseci )/fs ni wi qi bi fi-1ei-1 i ei-1 岩土边坡稳定性分析新方法与工程应用52iissiiiiiiiisiiiiiiiirtffffefffesincostancossinsincostancossin111力平衡方程力平衡方程（考虑两个方向力平衡条件及mohr-coul

58、omb准则）iiiiiiiiiiciiibcuqwkwrsectancossincosiiiiiciiiqwkwtsincossin条块抗滑力（除条间力的贡献）：条块抗滑力（除条间力的贡献）：条块下滑力（除条间力的贡献）：条块下滑力（除条间力的贡献）：岩土边坡稳定性分析新方法与工程应用53iisiiiiirtfee111siiiiiiiifffsincostancossin条间力递推方程：条间力递推方程：1111sincostancossinisiiiiiiiifffsiiiiiiiifff11111111sincostancossin0;00nee端部条件：端部条件：111111ninnij

59、jininnijjisttrrf安全计算公式：安全计算公式：岩土边坡稳定性分析新方法与工程应用54iiiiiciiiiiiiiiiiiihqhwkefefbbzebzesin22tan2tan21111力矩平衡方程：力矩平衡方程：111;iiiiiizemzem条间力矩定义：条间力矩定义：iiiiiciiiiiiiiiiihqhwkeebefefbmmsin2tan221111条间力矩递推方程：条间力矩递推方程：0;00nmm端部条件：端部条件：niiiiiiniiiiiiciiiiefefbhqhwkeeb11111sin2tan 计算公式：计算公式：岩土边坡稳定性分析新方法与工程应用55 no choose f(x) yes division of slices ri , ti assume fs0, 0 i , i-1 fs ei start fs 1 and 2 ? i , i-1 fs end 岩土边坡稳定性分析新方法与工程应用56 non-circular slip