劉 茂 楊紅娟 錢江澎

( ①四川省地質(zhì)工程勘察院 成都 610072)

( ②成都理工大學(xué)核技術(shù)與自動化工程學(xué)院 成都 610059)

0 引 言

英國著名土力學(xué)家畢肖普( Bishop A W) 于1955 年發(fā)表了圓弧條分法在土質(zhì)邊坡穩(wěn)定性中的應(yīng)用一文，Bishop 法由此正式誕生。Bishop( 1955)考慮條塊間力的作用，假設(shè)條間力為水平力，條塊土體界面上的剪應(yīng)力不能超過該界面的最大抗剪強(qiáng)度值，該法根據(jù)每個(gè)條塊的豎向力平衡條件和整體力矩平衡條件( 所有條塊所受作用力的力矩之和為0)推導(dǎo)出來的，由于條塊間切向力為未知量，Bishop A W 進(jìn)一步假設(shè)ΔHi= Hi+1－ Hi= 0，得到了簡化Bishop 法，簡化Bishop 法( John Krahn，2004; 李廣信等，2013) 滿足力多邊形閉合條件。

國內(nèi)外很多學(xué)者( Duncan，1996; 陳祖煜，2003; 鄭穎人等，2007; 蘇愛軍，2008) 對邊坡穩(wěn)定分析方法的計(jì)算精度進(jìn)行了討論和分析，簡化Bishop 法計(jì)算出的圓弧滑動面安全系數(shù)與嚴(yán)格條分法( 如Morgenstern-Price 法和Spencer 法) 的安全系數(shù)十分接近，一般誤差只有3%～4%，甚至在1%以內(nèi)。當(dāng)時(shí)，簡化Bishop法也可以看成嚴(yán)格條分法。Duncan( 1996) 認(rèn)為簡化Bishop 法在所有情況下都是精確的( 除了遇到數(shù)值分析困難情況外) ，其局限性表現(xiàn)在僅適用于圓弧滑裂面以及有時(shí)會遇到數(shù)值分析問題( 當(dāng)mai=cosαi+sinαitanφi/Fs＜0.2，可以認(rèn)為簡化Bishop 法中存在數(shù)值分析問題) 。

朱大勇等( 2007) 對簡化Bishop 法的嚴(yán)格性進(jìn)行了論證分析，簡化Bishop 法滿足所有平衡條件;劉茂( 2011) 對各種穩(wěn)定性計(jì)算方法進(jìn)行了總結(jié)，張年學(xué)等( 2018) 認(rèn)為條間剪力對其沒有明顯影響，鄧濤等( 2019) 基于通用條分原理，采用微分形式的土條受力平衡推導(dǎo)出了錨框支護(hù)邊坡地震動力穩(wěn)定分析方法。簡化Bishop 法被公認(rèn)為是一個(gè)計(jì)算圓弧形滑面穩(wěn)定性的簡便、精確的方法，該穩(wěn)定性計(jì)算方法已被廣泛應(yīng)用于圓弧形滑面的穩(wěn)定性分析中，但卻沒有比較簡易的剩余下滑推力解析方法，本文從簡化Bishop 法基本原理出發(fā)，利用力多邊形法則，推導(dǎo)出基于簡化Bishop 法的剩余下滑推力解析計(jì)算方法。

1 條塊i 的下滑推力模型

根據(jù)簡化Bishop 法的假設(shè)條件，條塊i 受力模型如下:

圖1 條塊i 的受力模型Fig. 1 The force model acting on slice i

根據(jù)條塊所受外力方向、條塊傾角、條塊所受力的固有特征，利用力多邊形法則和單個(gè)條塊處于靜力平衡條件可知，單個(gè)條塊力多邊形模型可分為兩種模型( Ⅰ型和Ⅱ型) 。

Ⅰ型為條塊所受抗滑力方向與坡向相反，再根據(jù)條塊i 的力多邊形模型，Ⅱ型可具體分為Ⅱ1型( 多邊形模型水平力之和大于等于0( 圖2 b、圖2d) ) 和Ⅱ2型( 多邊形模型水平力之和小于0( 圖2a、圖2c) ) ，典型模型見圖2。

具體模型表達(dá)式如下:

Ⅱ型為條塊所受抗滑力方向與坡向相同，再根據(jù)條塊i 的力多邊形模型，Ⅰ型可具體分為Ⅰ1型( 多邊形模型水平力之和大于等于0( 圖3 b、圖3d) ) 和Ⅰ2型( 多邊形模型水平力之和小于0( 圖3a、圖3c) ) ，典型模型見圖3。

具體模型表達(dá)式如下:

圖2 典型受力模型( Ⅰ型)Fig. 2 The typical force model of slice i( model Ⅰ)

圖3 典型受力模型( Ⅱ型)Fig. 3 The typical force model of slice i( model Ⅱ)

Ⅱ1型:

－Qisinθi+Ticosαi+( Ni+Ui) sinαi+KsWi≥0

Ⅱ2型:

－Qisinθi+Ticosαi+( Ni+Ui) sinαi+KsWi＜0

上述模型表達(dá)式中，Qi為外力( kN·m－1) ( 順時(shí)針為正，逆時(shí)針為負(fù)) ，Wi為條塊 i 的重 力( kN·m－1) ，Ti為條塊i 的抗滑力( kN·m－1) ( 與坡向相同為負(fù)，與坡向相反為正) ，Ni為條塊i 的垂直壓力( kN·m－1) ，Ui為條塊i 的孔隙水壓力( kN·m－1) ，Ks為水平加速度系數(shù)( g) ，αi為條塊i 的滑動面傾角( 與坡向相同為正，與坡向相反為負(fù)) ，θi為條塊i所受外力的傾角( 順時(shí)針為負(fù)，逆時(shí)針為正) 。

2 簡化Bishop 法剩余下滑推力解析方法

根據(jù)上述條塊i 作用力的力多邊形模型可知，對于條塊i，根據(jù)力多邊形法則，考慮垂直方向力的平衡時(shí)有:

根據(jù)滿足安全系數(shù)為Fs時(shí)的極限平衡條件可得:

將式( 2) 代入式( 1) ，整理后可得:

計(jì)算剩余下滑推力時(shí)，還必須滿足一定規(guī)定，規(guī)定如下:

(1) 土體抗拉強(qiáng)度很弱，可忽略不計(jì)，因此，剩余下滑推力不能為負(fù)值，若為負(fù)值時(shí)取為0。

(2) 上一條塊的剩余下滑推力直接作用于相鄰下一條塊上，作用方向?yàn)樗健?/p>

根據(jù)以上規(guī)定和兩種模型可知，條塊i 的剩余下滑推力Fi的解析表達(dá)式為式( 7) :

當(dāng)Fi小于等于0 時(shí)，條塊i 的剩余下滑推力為0，但剩余下滑力計(jì)算過程中不取0，按照實(shí)際計(jì)算值進(jìn)行計(jì)算。

當(dāng)Ti指向坡向方向時(shí)取負(fù)值，否則為正，Ti的判據(jù)表達(dá)式見式( 8) 。

其余公式和符號同上。

3 簡化Bishop 法剩余下滑推力與傳遞系數(shù)法隱式解剩余下滑推力比較與討論

3.1 圓弧形滑面模型

如圖4 所示的某圓弧滑動面的滑坡，該圓弧形滑動面滑坡模型數(shù)據(jù)見表1。

圖4 圓弧滑動面滑坡模型Fig. 4 Landslide model of the slip circle

表1 圓弧滑動面滑坡模型數(shù)據(jù)Table 1 The data of landslide model of the slip circle

3.2 剩余下滑推力對比分析

兩種方法的剩余下滑推力計(jì)算，模型采用上述圓弧形滑動面滑坡模型，滑體天然重度為19.6 kN·m－3，飽和重度為20.1 kN·m－3，內(nèi)黏聚力為18 kPa，內(nèi)摩察角為15°，分別得到各個(gè)安全系數(shù)ks=1.00 ～1.35 安全儲備下的剩余下滑推力，剩余下滑推力曲線如圖5、圖6。

圖5 簡化Bishop 法剩余下滑推力曲線Fig. 5 The curve map of residual sliding force based on simplified Bishop method

圖6 傳遞系數(shù)法剩余下滑推力曲線Fig. 6 The curve map of residual sliding force based on transfer coefficient method

由圖5 和圖6 可知，簡化Bishop 法剩余下滑推力最大值部位幾乎不隨安全系數(shù)變化，幾乎位于坡體的中部部位條塊，傳遞系數(shù)法隱式解剩余下滑推力最大值隨著安全系數(shù)增大從坡體中上部向坡體中部遷移，但出現(xiàn)最大值的條塊范圍分布范圍較小;整個(gè)坡體條塊的剩余下滑推力趨勢與傳遞系數(shù)法隱式解基本一致。

圖7 剩余下滑推力差值曲線Fig. 7 The difference curve map of residual sliding force between the simplified Bishop method and transport coefficient method

簡化Bishop 法剩余下滑推力比傳遞系數(shù)法隱式解剩余下滑推力大，兩者的剩余下滑推力差值在坡體的中上部( 約坡體的1/3) 隨著水平距離增大而增大，隨著安全系數(shù)增大，兩者的差值在減小; 安全系數(shù)為ks=1.0～1.15 時(shí)，坡體中部和中下部差值隨著安全系數(shù)增大差值呈增大的趨勢，只是在最后1～2 個(gè)條塊趨于減少趨勢，安全系數(shù)為ks=1.15 ～1.35時(shí)，坡體中部和中下部差值隨著安全系數(shù)增大差值呈減小的趨勢。

由圖8、圖9 可知，簡化Bishop 法相鄰條塊剩余下滑推力增量趨勢是一致的，從上往下總體呈現(xiàn)上升－下降－平緩－下降的趨勢，但是，簡化Bishop 法相鄰條塊剩余下滑推力增量在坡體上部( 1/3 部分) 隨著安全系數(shù)升高而增大，坡體其余部位幾乎不隨安全系數(shù)變化而變化，而傳遞系數(shù)法相鄰條塊剩余下滑推力增量各個(gè)條塊間隨著安全儲備的增大而增大，坡體中部基本上呈現(xiàn)線形增加。

由圖10 可知，簡化Bishop 法剩余下滑推力比傳遞系數(shù)法隱式解剩余下滑推力大，整個(gè)坡體條塊剩余下滑推力的平均值在坡體的中上部( 1/3 部位) ，簡化Bishop 法較傳遞系數(shù)法隱式解整個(gè)坡體條塊剩余下滑推力的增加百分比平均值一般可增加15%～35%，最大可增加45%，且隨著安全儲備增大，整個(gè)坡體條塊剩余下滑推力的增加百分比平均值在減小。

當(dāng)ks＜1.25 時(shí)，增加百分比從坡頂至坡腳呈增加－減少－平緩－快速增加的趨勢，且增加或減少的幅度隨著安全系數(shù)增加而減小; 當(dāng)ks＞1.25 時(shí)，增加百分比從坡頂至坡腳呈快速減少－減小的趨勢，且減少的幅度隨著安全系數(shù)增加而減小。

圖8 簡化Bishop 法相鄰條塊剩余下滑推力增量曲線Fig. 8 The residual slope thrust increment curve of adjacent bars by simplified Bishop method

圖9 傳遞系數(shù)法相鄰條塊剩余下滑推力增量曲線Fig. 9 The residual slope thrust increment curve of adjacent bars by transfer coefficient method

圖10 剩余下滑推力增減百分比曲線Fig. 10 The percentage increase curve map of residual sliding force when compared to transport coefficient method

表2 剩余下滑推力增減百分比Table 2 The percentage increase of residual sliding force when compared to transport coefficient method

當(dāng)安全系數(shù)較小( ks＜1.1) 時(shí)，坡體最后1 ～4 個(gè)條塊推力增幅較大，一般可增幅70%～100%，最大可增幅約10 倍甚至更大，該方法對于安全系數(shù)較小的邊坡最后幾個(gè)條塊的剩余下滑推力有待進(jìn)一步分析研究。

3.3 剩余下滑推力對比結(jié)果分析探討

針對以上剩余下滑推力的對比分析，總體上呈現(xiàn)出簡化Bishop 法剩余下滑推力高于傳遞系數(shù)法隱式解的剩余下滑推力，增大的幅度也不同，分析探討如下:

( 1) 對于同一個(gè)滑坡( 邊坡) ，兩種計(jì)算穩(wěn)定性的方法不一致，由于兩種穩(wěn)定性計(jì)算方法都有其局限性，簡化Bishop 法主要適用于圓弧形滑面滑坡，且只考慮條間力方向?yàn)樗椒较?，傳遞系數(shù)法主要受相鄰條塊傾角差和突變的影響，且滑面傾角不宜過大，評價(jià)出的滑坡( 邊坡) 穩(wěn)定性不同，安全系數(shù)和穩(wěn)定系數(shù)的差值，即安全儲備大小影響剩余下滑推力值。

(2) 簡化Bishop 法剩余下滑推力比傳遞系數(shù)法隱式解剩余下滑推力大，究其原因在于簡化Bishop法本身忽略了條塊間的豎向摩擦力，而傳遞系數(shù)法考慮了條塊間的豎向摩擦力，且考慮了相鄰條塊滑面傾角的影響。

( 3) 隨著安全儲備的增大，兩者的剩余下滑推力差值減小，原因在于這兩種計(jì)算方法都為隱式解法，同時(shí)折減c、tanΦ，安全系數(shù)變大，條塊的垂直壓力減小，下滑力增大可減少條塊間的相互影響，可減小兩者剩余下滑力的差值。

4 算 例

本次選取一個(gè)經(jīng)典算例進(jìn)行簡化Bishop 法和傳遞系數(shù)法剩余下滑推力比較，邊坡剖面和土層參數(shù)( 圖11) ，該邊坡假設(shè)為永久邊坡，邊坡的安全等級為一級邊坡，取穩(wěn)定安全系數(shù)為ks=1.35，計(jì)算得到簡化Bishop 法和傳遞系數(shù)法剩余下滑推力值( 圖12) 和兩者的剩余下滑推力增減百分比曲線( 圖13) 。

圖11 典型算例的地質(zhì)剖面圖Fig. 11 The geological profile of a typical example

5 結(jié) 論

( 1) 評價(jià)圓弧形滑動面滑坡時(shí)，簡化Bishop 法是一種公認(rèn)精確、簡便的“嚴(yán)格”方法，該方法在我國被普遍應(yīng)用，但卻沒有與之相適應(yīng)的剩余下滑推力解析計(jì)算方法。

( 2) 根據(jù)條塊的受力條件、力的多邊形法則和Ti的方向( 指向坡向方向時(shí)取負(fù)值，否則為正) ，把單個(gè)條塊力多邊形模型分為兩類:Ⅰ型( Ti≥0) 和Ⅱ型( Ti＜0) ，依據(jù)上述兩類模型，利用條塊垂直力平衡和極限平衡條件以及條塊水平力平衡推導(dǎo)出了剩余下滑推力的解析計(jì)算公式。

圖12 簡化Bishop 法和傳遞系數(shù)法剩余下滑推力曲線Fig. 12 The curve map of residual sliding force based on simplified Bishop method and transfer coefficient method

圖13 典型算例剩余下滑推力增減百分比曲線Fig. 13 The difference curve map of residual sliding force between the simplified Bishop method and transport coefficient method about typical example

(3) 簡化Bishop 法和傳遞系數(shù)法隱式解剩余下滑推力比較可知，簡化Bishop 法剩余下滑推力比傳遞系數(shù)法隱式解剩余下滑推力大，簡化Bishop 法整個(gè)坡體條塊剩余下滑推力的增加百分比平均值一般可增加15%～35%，最大可增加45%，且隨著安全儲備增大，整個(gè)坡體條塊剩余下滑推力的增加百分比平均值在減小。當(dāng)安全系數(shù)較小( ks＜1.1) 時(shí)，坡體最后1～4 個(gè)條塊推力增幅較大，一般可增幅70%～100%，最大可增幅約10 倍甚至更大，該方法對于安全系數(shù)較小的滑坡( 邊坡) 最后幾個(gè)條塊的剩余下滑推力有待進(jìn)一步分析研究。

( 4) 該剩余下滑推力解析法為隱式解法，且隱式解法屬于強(qiáng)度儲備系數(shù)法，該方法比較符合滑坡( 邊坡) 下滑破壞的規(guī)律，為圓弧形滑動面滑坡( 邊坡) 的剩余下滑推力提供了一定的理論依據(jù)。