鄧東平，李 亮

（中南大學(xué) 土木工程學(xué)院，長沙 410075）

1 引 言

在邊坡穩(wěn)定性分析方法中，極限平衡法作為一種非常有效的分析方法被廣泛應(yīng)用[1]。目前，極限平衡法多為豎直條分法[2－10]，且在安全系數(shù)求解時對條間力所建立的假設(shè)條件也是針對豎直條分法[11－14]。然而，豎直條分法盡管有其優(yōu)勢，但當(dāng)邊坡成層巖土體時，豎向條分體內(nèi)部受巖土體性質(zhì)變化影響較大，這樣會給計算帶來極大的不便[15]，因此，這時水平條分法更具有適用性。

文獻(xiàn)[16]在一般條分時通過借用豎直條分法的條間力假設(shè)條件進(jìn)行了安全系數(shù)的推導(dǎo)，發(fā)現(xiàn)當(dāng)條分不再為豎直條分時，計算得的安全系數(shù)與豎直條分法相差較大。同時，在針對水平條分法研究中，國內(nèi)外大多數(shù)文獻(xiàn)關(guān)注于有加固措施的邊坡穩(wěn)定性，如加筋土[17－19]、錨桿（索）[20]等。文獻(xiàn)[1, 15]指出，這些研究沒有針對普通邊坡的穩(wěn)定性分析，而且要么有些計算假定過于復(fù)雜，不利于工程應(yīng)用，要么有些過于簡單，使計算結(jié)果誤差較大。另外，目前針對一般邊坡穩(wěn)定性的水平條分法，所涉及到的算例分析一般是認(rèn)為滑動面比較適合水平條分劃分，而對于一般型式下的滑動面，如存在一段弓形體時，水平條分則表現(xiàn)出條分界限不明，弓形體以下水平條分出現(xiàn)兩個滑動底面的缺點(diǎn)。

綜上所述可知，在邊坡穩(wěn)定性水平條分法極限平衡分析中，存在如下不足：①過少對水平條分的條間力進(jìn)行分析；②沒有將豎直條分法與水平條分法的條間力假設(shè)條件相對應(yīng)；③沒有在滑動面為一般型式的條件下，對水平條分法如何進(jìn)行穩(wěn)定性分析采取相應(yīng)處理方法。

本文通過對水平條分和豎直條分中條間力的關(guān)系進(jìn)行分析，建立起水平條分法與豎直條分法相對應(yīng)的條間力假設(shè)條件，并對滑動面存在弓形體的一般情況，將水平條分法采用斜條分與水平條分相結(jié)合的模式進(jìn)行處理。在此基礎(chǔ)上，推導(dǎo)出了一些方法在水平條分時的安全系數(shù)計算公式。當(dāng)使用圓弧滑動面和任意曲線滑動面兩種型式時，通過算例分析，并與豎直條分法和以往研究成果相對比，驗證了本文針對水平條分研究的可行性。

2 水平條分存在的問題及處理方法

在水平條分法中，以往研究成果對如下問題研究不足：①條間力的分析，如圖 1(a)所示；②滑動面型式?jīng)]有針對更為普遍的情況，如圖1(b)、(c)、(d)所示滑動面存在一段弓形體。

圖1 水平條分條間力和存在弓形體的滑動面Fig.1 Inter-slice forces of horizontal slice and sliding surface with presence of toxoplasma

2.1 水平條分條間力分析

當(dāng)水平條分法直接使用豎直條分法對條間力的假定時需研究其可行性。如圖2所示，abcd為水平土條為豎直土條分別為水平土條 i上下底面的條間法向力和剪切力，分別為豎直土條i左右側(cè)面的條間法向力和剪切力，分別為水平土條i和豎直土條 i共有滑動底面的法向力和剪切力，分別為水平土條i和豎直土條i的重力。

圖2 水平條分條間力與豎直條分條間力關(guān)系Fig.2 Relation of inter-slice forces between horizontal slice and vertical slice

圖3 水平條分條間力假設(shè)條件分析Fig.3 Assumption analysis of inter-slice forces of horizontal slice

（2）對于簡化Bishop法和簡化Janbu法，豎直條分：

（3）對于M-P法，

2.2 一般滑動面情況下的處理方法

從圖1（b）、（c）、（d）中的3種滑動面情況可知：當(dāng)滑動面存在一段弓形體時，若使用水平條分，則會出現(xiàn)兩個滑動底面，因而對該處土條的受力求解帶來麻煩，并造成水平條分界限不清的缺點(diǎn)。為此，本文采用如圖4（a）所示的條分法，在下滑點(diǎn)A及其水平延長線以下的弓形體采用水平間距為且過A點(diǎn)的斜條分（條分?jǐn)?shù)為n2），下滑點(diǎn)A及其水平延長線以上仍然采用水平條分（條分?jǐn)?shù)為n1）。

圖4 一般滑動面情況下水平條分法的處理Fig.4 Treatment of horizontal slice method undergeneral situation of sliding surface

按照2.1節(jié)所述，在斜條分中，為了建立起與豎直條分相等效的合理條間力假設(shè)條件，本文以式（2）為基礎(chǔ)，得到了如下相應(yīng)計算方法的斜條條間力假設(shè)條件。

（1）對于瑞典法：

（2）對于簡化Bishop法和簡化Janbu法：

（3）對于M-P法：

3 安全系數(shù)計算方法

如圖5所示，水平條分的土條i上各力對滑動底面法向、切向的力平衡和對滑動底面中心點(diǎn)的力矩平衡如式（3），斜條分的土條j上各力對滑動底面法向、切向的力平衡和對滑動底面中心點(diǎn)的力矩平衡如式（4）。

圖5 水平條分和斜條分受力分析Fig.5 Stress analysis of horizontal slice and oblique slice

（1）瑞典法

由單個土條體的x軸和y軸平衡，以及滑動土體對滑動圓弧圓心的力矩和滿足或土條底面下滑力和抗滑力的總和滿足分別得到安全系數(shù)的計算式（5）、（6）。

（2）簡化Bishop法

由單個土條體的y軸平衡，以及滑動土體對滑動圓弧圓心的力矩和滿足得到安全系數(shù)的計算公式如式（7）。

（3）簡化Janbu法

由單個土條體的x軸和y軸平衡，得到安全系數(shù)的計算公式如式（9）。

同時，有些文獻(xiàn)[21－22]對豎直條分法中簡化Janbu法進(jìn)行了改進(jìn)，即將式（9）分子中的直接乘入到分母，本文也對式（9）進(jìn)行同樣改進(jìn)，得到安全系數(shù)的計算公式如式（10）。

（4）M-P法

由單個土條體的滑動底面法向和切向力平衡，以及對土條滑動底面中心點(diǎn)的力矩平衡，得到安全系數(shù)的計算公式如式（11），及式中水平條分和斜條分時的各參數(shù)量的計算如式（12）～（21）。

水平條分時：

斜條分時：

水平條分時

斜條分時

4 計算對比分析

4.1 均質(zhì)邊坡

[23]，選擇 3個算例進(jìn)行對比分析，如圖6～8所示。其中，對于任意曲線滑動面，采取文獻(xiàn)[21]方法。算例參見表1。

算例1～3計算得的最小安全系數(shù)見表2，得到的臨界滑動面如圖6～8所示。

從表2可知：本文各方法在水平條分時，計算得的安全系數(shù)與豎直條分法和文獻(xiàn)[23]方法較接近，說明本文對水平條分不足所采取的處理方法可行；水平條分法時，兩種瑞典法計算的安全系數(shù)基本一致；水平和豎直條分法下的一般簡化Janbu法（即簡化Janbu法式（9））計算得的安全系數(shù)偏小，改進(jìn)的簡化Janbu法（即簡化Janbu法式（10））計算得的安全系數(shù)與簡化Bishop法和M-P法頗為接近，說明改進(jìn)簡化Janbu法對邊坡的穩(wěn)定性評判更準(zhǔn)確些。

從圖6～8可知：水平條分法和豎直條分法下，本文各方法得到的臨界任意曲線滑動面和臨界圓弧滑動面與文獻(xiàn)[23]得到的臨界滑動面較接近；臨界任意曲線滑動面近似于圓弧滑動面，但表現(xiàn)為非圓弧型；水平條分下的兩種瑞典法得到的臨界圓弧滑動面基本相同；一般簡化Janbu法得到的臨界滑動面的范圍較其他方法要大。

表1 算例1～3邊坡參數(shù)Table 1 Slope parameters of examples 1-3

表2 均質(zhì)邊坡算例中最小安全安全系數(shù)對比Table 2 Contrast of the minimum FOS in homogeneous slope examples

圖6 算例1中的臨界滑動面Fig.6 Critical sliding surfaces in example 1

圖7 算例2中的臨界滑動面Fig.7 Critical sliding surfaces in example 2

圖8 算例3中的臨界滑動面Fig.8 Critical sliding surfaces in example 3

4.2 分層邊坡

分別參考文獻(xiàn)[1, 24－25]，選擇3個算例進(jìn)行對比分析，如圖9～11所示。其中，對于任意曲線滑動面，采取文獻(xiàn)[21]方法。

算例4：本邊坡由3層水平土層組成，其下層為基巖。計算參數(shù)見表3。在文獻(xiàn)[1]中，假定滑動面為圓弧滑動面時，瑞典法計算得的最小安全系數(shù)為1.321，Bishop法為1.355，Chen-M法為1.357，文獻(xiàn)[1]法為1.367。

算例5：本邊坡由2層土層組成，計算參數(shù)見表3。在文獻(xiàn)[24]中，假定滑動面為螺旋線時，由剛體單元上限法得出的最小安全系數(shù)為1.01。

算例6：該邊坡由2層水平土層組成，其計算參數(shù)見表3。在文獻(xiàn)[25]中，假定滑動面為分段螺旋線時，通過極限上限分析法計算得的最小安全系數(shù)為1.000。

算例4～6計算得的最小安全系數(shù)見表4，得到的臨界滑動面如圖9～11所示。

表3 算例4～6邊坡參數(shù)Table 3 Slope parameters of examples 4-6

表4 分層邊坡算例中最小安全系數(shù)對比Table 4 Contrast of the minimum FOS in layered slope examples

從表4可知：在分層土坡中，本文水平條分法下的各方法計算得的安全系數(shù)與豎直條分法和文獻(xiàn)方法相也非常接近，兩種瑞典法計算得的安全系數(shù)也基本一致；改進(jìn)的簡化 Janbu法（即簡化 Janbu法式（10））計算得的安全系數(shù)仍較一般簡化Janbu法（即簡化Janbu法式（9））優(yōu)。

從圖9～11可知：本文水平條分法和豎直條分法下各方法得到的臨界滑動面（包括任意曲線滑動面和圓弧滑動面）與文獻(xiàn)得到的臨界滑動面相接近，臨界任意曲線滑動面相似于圓弧滑動面，但表現(xiàn)為非圓弧型；水平條分法下兩種瑞典法得到的臨界圓弧滑動面仍基本相同。

圖9 算例4中的臨界滑動面Fig.9 Critical sliding surfaces in example 4

圖10 算例5中的臨界滑動面Fig.10 Critical sliding surfaces in example 5

圖11 算例6中的臨界滑動面Fig.11 Critical sliding surfaces in example 6

5 結(jié) 論

（1）本文對水平條分和豎直條分中條間力的關(guān)系進(jìn)行了分析，建立起水平條分法與豎直條分法相對應(yīng)的條間力假設(shè)條件，并對水平條分遇到滑動面存在一段弓形體時，采用取斜條分與水平條分相結(jié)合模式，通過算例對比分析可知，文中各方法在水平條分時計算得的安全系數(shù)與豎直條分法和以往研究成果相接近，得到的臨界滑動面也頗為相似，說明本文針對水平條分不足所采取的處理方法可行；

（2）任意曲線滑動面方法與圓弧滑動面方法計算得的安全系數(shù)接近，但較圓弧滑動面方法整體上偏小，且得到的臨界滑動面與圓弧滑動面相似，但表現(xiàn)為非圓弧型；

（3）水平條分法中兩種瑞典法計算得的安全系數(shù)和臨界滑動面基本一致，豎直和水平條分法中一般Janbu法計算得的安全系數(shù)偏小，改進(jìn)簡化Janbu法計算得的安全系數(shù)與簡化Bishop法和M-P法頗為接近，且改進(jìn)簡化Janbu法較一般Janbu法得到的臨界滑動面與其他方法更接近一些。

參 考 文 獻(xiàn)

[1]陳昌富, 楊宇. 邊坡穩(wěn)定性分析水平條分法及其進(jìn)化計算[J]. 湖南大學(xué)學(xué)報(自然科學(xué)版), 2004, 31(3): 72－75.CHEN Chang-fu, YANG Yu. Horizontal slice method of slope stability analysis based on evolutionary algorithm[J].Journal of Hunan University(Natural Sciences), 2004,31(3): 72－75.

[2]朱大勇, 鄧建輝, 臺佳佳. 簡化 Bishop法嚴(yán)格性的論證[J]. 巖石力學(xué)與工程學(xué)報, 2007, 26(3): 455－458.ZHU Da-yong, DENG Jian-hui, TAI Jia-jia. Theoretical verification of rigorous nature of simplified Bishop method[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(3): 455－458.

[3]朱大勇, 李焯芬, 黃茂松, 等. 對 3種著名邊坡穩(wěn)定性計算方法的改進(jìn)[J]. 巖石力學(xué)與工程學(xué)報, 2005, 24(2):183－194.ZHU Da-yong, LEE C F, HUANG Mao-song, et al.Modifications to three well-known methods of slope stability analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(2): 183－194.

[4]鄒廣電. 邊坡穩(wěn)定分析條分法的一個全局優(yōu)化算法[J].巖土工程學(xué)報, 2002, 24(3): 309－312.ZOU Guang-dian. A global optimization method of the slice method for slope stability analysis[J]. Chinese Journal of Geotechnical Engineering, 2002, 24(3): 309－312.

[5]萬文, 曹平, 馮濤, 等. 基于加速混合遺傳算法搜索復(fù)雜邊坡的最危險滑動面[J]. 巖土工程學(xué)報, 2006, 28(4):475－479.WAN Wen, CAO Ping, FENG Tao, et al. Searching for the most dangerous failure surface of complex slope based on accelerating hybrid genetic algorithm[J].Chinese Journal of Geotechnical Engineering, 2006,28(4): 475－479.

[6]陳祖煜. 土坡穩(wěn)定分析通用條分法及其改進(jìn)[J]. 巖土工程學(xué)報, 1983, 5(4): 11－27.CHEN Zu-yu. The generalized method of slices for slope stability analysis and its modifications[J]. Chinese Journal of Geotechnical Engineering, 1983, 5(4): 11－27.

[7]方玉樹. 邊坡穩(wěn)定性分析條分法最小解研究[J]. 巖土工程學(xué)報, 2008, 30(3): 331－335.FANG Yu-shu. The lowest solution of slice method for slope stability analysis[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(3): 331－335.

[8]吳振君, 葛修潤. 求解邊坡矢量和安全系數(shù)的條分法[J]. 巖土力學(xué), 2009, 30(8): 2337－2342.WU Zhen-jun, GE Xiu-run. Solving vector sum factor of safety of slope by method of slices[J]. Rock and Soil Mechanics, 2009, 30(8): 2337－2342.

[9]童志怡, 陳從新, 徐健, 等. 邊坡穩(wěn)定性分析的條塊穩(wěn)定系數(shù)法[J]. 巖土力學(xué), 2009, 30(5): 1393－1398.TONG Zhi-yi, CHEN Cong-xin, XU Jian, et al. A slice-stability method for stability analysis of slopes[J].Rock and Soil Mechanics, 2009, 30(5): 1393－1398.

[10]高玉峰, 范昭平. 多點(diǎn)、多向地震作用下非圓弧滑面邊坡穩(wěn)定分析通用條分法[J]. 巖土力學(xué), 2010, 31(12):3816－3822.GAO Yu-feng, FAN Zhao-ping. Generalized slice method for stability analysis of slope with non-circular slide face under multipoint and multidirection seismic ground motions[J]. Rock and Soil Mechanics, 2010, 31(12):3816－3822.

[11]殷宗澤. 土工原理[M]. 北京: 中國水利水電出版社,2007.

[12]BISHOP A W. The use of slip circle for the stability analysis of slopes[J]. Geotechnique, 1955, 5(1): 7－17.

[13]MORGENSTERN N R, PRICE V E. The analysis of the stability of general slip surfaces[J]. Geotechnique, 1965,15(l): 79－93.

[14]SPENCER E. A method of analysis of the stability of embankments assuming parallel inter-slice forces[J].Geotechnique, 1967, 17(1): 11－26.

[15]陳昌富, 楊宇, 龔曉南. 基于遺傳算法地震荷載作用下邊坡穩(wěn)定性分析水平條分法[J]. 巖石力學(xué)與工程學(xué)報,2003, 22(11): 1919－1923.CHEN Chang-fu, YANG Yu, GONG Xiao-nan.Horizontal slice method of slope stability analysis under seismic load based on genetic algorithm[J]. Chinese Journal of Rock Mechanics and Engineering, 2003,22(11): 1919－1923.

[16]楊明成. 基于力平衡求解安全系數(shù)的一般條分法[J].巖石力學(xué)與工程學(xué)報, 2005, 24(7): 1216－1221.YANG Ming-cheng. A general slice method for determination of factors of safety based on force equilibrium[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(7): 1216－1221.

[17]SHAHGHOLI M, FAKHER A, JONES C J F P.Horizontal slice method of analysis[J]. Géotechnique,2001, 51(10): 881－885.

[18]NOURI H, FAKHER A, JONES C J F P. Development of horizontal slice method for seismic stability analysis of reinforced slopes and walls[J]. Geotextiles and Geomembranes, 2006, 24(3): 175－187.

[19]NOURI H, FAKHER A, JONES C J F P. Evaluating the effects of the magnitude and amplification of pseudo-static acceleration on reinforced soil slopes and walls using the limit equilibrium horizontal slices method[J]. Geotextiles and Geomembranes, 2008, 26(3):263－278.

[20]魏作安, 朱彬, 萬玲, 等. 邊坡錨固工程分析中的水平條分法[J]. 土木建筑與環(huán)境工程, 2010, 32(3): 53－56.WEI Zuo-an, ZHU Bin, WAN Lin, et al. Horizontal slice method for the stability analysis of reinforced slopes[J].Journal of Civil Architectural and Environmental Engineering, 2010, 32(3): 53－56.

[21]鄧東平, 李亮, 趙煉恒. 基于 Janbu法的邊坡整體穩(wěn)定性滑動面搜索新方法[J]. 巖土力學(xué), 2011, 32(3): 891－898.DENG Dong-ping, LI Liang, ZHAO Lian-heng. A new method of sliding surface searching for general stability of slope based on Janbu method[J]. Rock and Soil Mechanics, 2011, 32(3): 891－898.

[22]鄧東平, 李亮, 趙煉恒. 穩(wěn)定滲流條件下土坡穩(wěn)定性分析的一種新方法[J]. 工程地質(zhì)學(xué)報, 2011, 19(1): 29－36.DENG Dong-ping, LI Liang, ZHAO Lian-heng. A new method of stability analysis of slope under the condition of steady seepage[J]. Journal of Engineering Geology,2011, 19(1): 29－36.

[23]Rocscience Inc. Slide verification manual[R]. Toronto:Rocscience Inc., 2003.

[24]王根龍, 伍法權(quán), 張軍慧. 非均質(zhì)土坡穩(wěn)定性分析評價的剛體單元上限法[J]. 巖石力學(xué)與工程學(xué)報, 2008,28(增刊 2): 3425－3430.WANG Gen-long, WU Fa-quan, ZHANG Jun-hui. Upper bound approach of rigid elements for inhomogeneous soil slope stability analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 28(Supp. 2): 3425－3430.

[25]KUMAR J, SAMUI P. Stability determination for layered soil slopes using the upper bound limit analysis[J].Geotechnical and Geological Engineering, 2006, 24(6):1803－1819.