孫朝燚　陳從新　鄭允　張偉　馬力　張海娜　張亞鵬

摘???要：針對巖質(zhì)反傾邊坡滑動傾倒復(fù)合破壞的研究不足，首先建立了邊坡滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞的地質(zhì)模型，具體可分為滑動區(qū)、塊狀傾倒區(qū)和彎曲傾倒區(qū);然后根據(jù)三個分區(qū)巖塊的破壞機(jī)制，提出了各個分區(qū)的力學(xué)模型;并基于極限平衡理論和懸臂梁模型，提出了邊坡滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞的逐步分析方法;最后通過一個工程實例驗證了所提地質(zhì)模型和分析方法的正確性.?研究結(jié)果表1明：巖質(zhì)反傾邊坡滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞的穩(wěn)定性由滑動-塊狀傾倒復(fù)合破壞區(qū)域控制;塊狀傾倒區(qū)域?qū)儆谥鲃悠茐膮^(qū)域，滑動區(qū)屬于被動破壞區(qū)域，治理加固時應(yīng)重點加固塊狀傾倒破壞區(qū).

關(guān)鍵詞：巖質(zhì)反傾邊坡;橫向節(jié)理;邊坡穩(wěn)定性;數(shù)值模擬;分析方法

中圖分類號：P642????????????????????????????????文獻(xiàn)標(biāo)志碼：A

Analysis?of?Anti-dip?Rock?Slopes?Against?Composite?Toppling?Failure

SUN?Chaoyi1，2，CHEN?Congxin1，ZHENG?Yun1？，ZHANG?Wei1，

MA?Li1，ZHANG?Haina1，2，ZHANG?Yapeng1，2

（1.?State?Key?Laboratory?of?Geomechanics?and?Geotechnical?Engineering，Institute

of?Rock?and?Soil?Mechanics，Chinese?Academy?of?Sciences，Wuhan?430071，China;

2.?University?of?Chinese?Academy?of?Sciences，Beijing?100049，China）

Abstract：Aiming?at?the?insufficient?research?on?the?complex?sliding-toppling?failure?of?anti-dip?rock?slopes，the?geological?model?of?complex?sliding-block-flexural?toppling?failure?was?constructed，firstly，which?contains?the?sliding，block?toppling?and?flexural?toppling?zones.?Next，based?on?the?failure?mechanisms?of?these?three?zones，their?mechanical?models?were?established，respectively.?Then，the?stability?analysis?approach?of?the?complex?sliding-block-flexural?toppling?failure?was?proposed?according?to?the?limit?equilibrium?theory?and?the?cantilever?beam?model.?Finally，one?case?study?was?performed?for?practical?verification?on?the?proposed?geomechanics?models?and?analysis?approach.?The?results?show?that?the?stability?of?the?complex?sliding-block-flexural?toppling?zones?is?controlled?by?the?block-flexural?toppling?zones.?Further，the?block?and?flexural?toppling?zones?are?active?and?passive?failure?zones，respectively.?Consequently，the?block?toppling?zone?plays?a?vital?role?in?slope?design?treatment.

Key?words：anti-dip?rock?slopes;transverse?joints;slope?stability;numerical?simulation;analysis?approach

孫朝燚等：巖質(zhì)反傾邊坡復(fù)合傾倒破壞分析

傾倒破壞是巖質(zhì)反傾邊坡的主要失穩(wěn)模式之一，常見于露天采礦、水利水電、公路鐵路等各類工程邊坡中，這些邊坡的失穩(wěn)給工程建設(shè)和人民生活造成了嚴(yán)重的危害[1-6].?Goodman和Bray[7]從單個巖層破壞的力學(xué)機(jī)制出發(fā)，將傾倒破壞分為彎曲傾倒、塊狀傾倒和塊狀-彎曲復(fù)合傾倒三種類型.?若巖質(zhì)邊坡中存在一組與邊坡走向近似、插入坡內(nèi)的主控裂隙組如層理、頁理等，邊坡類似疊合懸臂梁，由于重力作用，巖塊承受著拉伸和壓縮彎曲應(yīng)力，當(dāng)每個巖塊中的拉應(yīng)力超過巖塊的抗拉強(qiáng)度時，邊坡就會發(fā)生傾倒失穩(wěn)破壞，這種破壞模式稱為彎曲傾倒破壞.?若上述邊坡中還存在一組與主控裂隙組成大角度相交的橫向裂隙組，此時巖塊不抗拉，由于重力作用，巖塊可能繞底面發(fā)生轉(zhuǎn)動，這種破壞模式稱為塊狀傾倒破壞.?塊狀-彎曲復(fù)合傾倒破壞是兩種破壞模式的組合，指邊坡中既有塊狀傾倒的巖塊，也有發(fā)生彎曲傾倒破壞的巖塊.

對于上述三種傾倒破壞模式，國內(nèi)外很多學(xué)者已做了大量工作并取得了顯著的成果[8-17].?然而由于巖石的脆性和節(jié)理的不規(guī)則性，巖質(zhì)反傾邊坡還可能發(fā)生滑動-傾倒復(fù)合破壞.?Alejano等[1]通過西班牙的一個失穩(wěn)礦山邊坡，詳細(xì)分析了邊坡上部塊狀傾倒-下部圓弧滑動的破壞機(jī)制，并提出了相應(yīng)的穩(wěn)定性分析方法.?Amini和Mohtarami等[18-19]通過室內(nèi)模型試驗，研究了邊坡上部圓弧滑動-下部彎曲傾倒破壞的復(fù)合破壞機(jī)制，并推導(dǎo)了安全系數(shù)計算公式.?當(dāng)橫向節(jié)理或裂隙并未由坡腳貫通到坡頂時，坡腳巖層由于巖塊高寬比較小，可能沿橫向節(jié)理發(fā)生平面滑動破壞，中部巖層可能發(fā)生塊狀傾倒破壞，而上部巖層由于未被橫向節(jié)理切割可能發(fā)生彎曲傾倒破壞，這類破壞可稱之為滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞.

本文首先建立邊坡滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞的地質(zhì)模型和力學(xué)模型，然后基于極限平衡理論和懸臂梁模型提出邊坡復(fù)合破壞的穩(wěn)定性分析方法，最后通過山西渾源某花崗巖礦山邊坡來驗證所提地質(zhì)力學(xué)模型和分析方法的正確性.

1???復(fù)合傾倒破壞分析

1.1???復(fù)合傾倒破壞過程

底部橫向節(jié)理未完全貫通的傾倒破壞過程如圖1所示.?橫向節(jié)理將巖層切割成離散的巖塊，巖塊區(qū)域的上部巖塊由于高寬比大，在自重和層間推力作用下發(fā)生塊狀傾倒變形，推擠下部巖塊，發(fā)生滑動-塊狀傾倒復(fù)合破壞;未被橫向節(jié)理貫通的頂部巖層失去前緣巖塊的支撐抗力之后形成懸臂段，在自重作用下發(fā)生懸臂彎曲，形成彎曲傾倒破壞;邊坡完全傾倒破壞后，形成“L”型的破壞形態(tài).

復(fù)合傾倒破壞主要為邊坡前緣和中部巖塊發(fā)生滑動-塊狀傾倒復(fù)合破壞，而后伴隨后緣巖層發(fā)生彎曲傾倒破壞.?因此，可將此類復(fù)合破壞分為3個分區(qū)：滑動區(qū)、塊狀傾倒區(qū)和彎曲傾倒區(qū)（如圖1所示）.

1.2???地質(zhì)幾何模型

巖質(zhì)反傾邊坡復(fù)合傾倒破壞的分析模型如圖2所示，圖中h為邊坡高度，hi為巖塊i的高度，t為巖塊厚度，α為邊坡坡度，θ為裂隙帶傾角，β為反傾結(jié)構(gòu)面法線傾角，θr為裂隙帶與反傾結(jié)構(gòu)面法線夾角，θr?=?θ?-?β，η為反傾結(jié)構(gòu)面傾角，η?=?90°?-?β.

破壞面以上任意巖塊的自重為：

wi?=?γhi?t.??????????????????????????????（1）

式中：γ為巖體重度;i為自坡腳向上的四邊形巖塊編號，1≤i≤itotal（1為橫向節(jié)理貫通的坡腳第一個巖塊，itotal為橫向節(jié)理貫通的最后一個巖塊）.

1.3???力學(xué)分析方程

復(fù)合傾倒破壞是一個十分復(fù)雜的過程，要想構(gòu)建完全精確的力學(xué)分析方程是不太可能的，但是可以通過合理的假設(shè)和簡化得到滿足工程精度要求的分析方程.?本文沿用了以下幾條基本假設(shè)來簡化復(fù)合傾倒破壞的分析過程[7-15].

1）底部橫向節(jié)理起始于坡腳位置，且滑動區(qū)和塊狀傾倒區(qū)巖塊沿橫向節(jié)理發(fā)生破壞.

2）巖塊發(fā)生滑動破壞時，假設(shè)巖層間無相互錯動，只考慮層面法向力Pi，底滑面滿足極限摩擦平衡條件Si?=?Ni?tan?φ1，Si、Ni和φ1分別為滑動破壞面的切向力、法向力和摩擦角.

3）巖塊發(fā)生塊狀傾倒破壞時，層間相互作用力簡化成集中力，作用點位于巖塊界面最上端，層面滿足極限摩擦平衡條件Qi?=?Ti?tan?φ2，Qi、Ti和φ2分別為傾倒層面的切向力、法向力和摩擦角.

4）所有潛在破壞巖層具有相同的安全系數(shù)，且

都等于邊坡的整體安全系數(shù)Fs?.

根據(jù)反傾邊坡的復(fù)合傾倒過程，從坡腳向坡頂?shù)姆较蛑鸩綄r塊進(jìn)行穩(wěn)定性分析，首先對坡腳第1個巖塊進(jìn)行力學(xué)分析，如圖3所示.

對于坡腳第1個巖塊，按照剪切滑移進(jìn)行分析，沿破壞面和垂直破壞面建立平衡方程如下：

w1sin?θ?+?P1cos?θr?=?N1tan?φ1?/Fs，???（2）

N1?=?w1cos?θ?-?P1sin?θr?.????（3）

整理上式可得發(fā)生剪切滑移破壞時所需第2個巖塊作用的最小推力為：

P1?=??.????（4）

對于坡腳第1個巖塊，按照塊狀傾倒破壞進(jìn)行分析，以巖塊破壞面下端點為轉(zhuǎn)折點，建立力矩平衡方程如下：

T1·h1+w1sin?β?·0.5h1-w1cos?β?·0.5t=T1tan?φ2?/Fs·t.

（5）

整理上式可得發(fā)生塊狀傾倒破壞時所需第2個巖塊作用的最小推力為：

T1?=?0.5w1cos?β?.??????（6）

P1、T1的相對大小決定坡腳第1個巖塊的潛在破壞模式.?若P1?≤?T1，則坡腳第1個巖塊潛在的破壞模式為剪切滑移破壞;若P1?>?T1，則坡腳第1個巖塊潛在的破壞模式為塊狀傾倒破壞.

同理，巖塊i按照剪切滑移進(jìn)行分析時，力學(xué)模型如圖4（a）所示，沿破壞面和垂直破壞面建立平衡方程如下：

wi?sin?θ+Pi?cos?θr?=Ni?tan?φ1?/Fs?+min（Ti-1，Pi-1）cos?θr，

（7）

Ni?=wi?cos?θ-Pi?sin?θr?+min（Ti-1，Pi-1）sin?θr?.??????（8）

整理上式可得發(fā)生剪切滑移破壞時所需上部巖塊作用的最小推力為：

Pi?=?min（Ti-1，Pi-1）+.

（9）

巖塊i按照塊狀傾倒破壞進(jìn)行分析時，力學(xué)模型如圖4（b）所示，以巖塊破壞面下端點為轉(zhuǎn)折點，建立力矩平衡方程如下：

Ti·hi?-?min（Ti-1，Pi-1）·hi-1?+?wi?（sin?β?·0.5hi?-

cos?β?·0.5t）=Ti?tan?φ2?/Fs·t.?????（10）

整理上式可得發(fā)生塊狀傾倒破壞時所需上部巖塊作用的最小推力為：

Ti??=?.

（11）

若塊狀傾倒區(qū)中先有部分巖塊和滑動區(qū)巖塊發(fā)生滑動-塊狀傾倒復(fù)合破壞，塊狀傾倒區(qū)后部存在殘留巖塊，巖塊處于臨空狀態(tài)，層間作用力變?yōu)榱?，力學(xué)模型如圖5所示，以巖塊破壞面下端點為轉(zhuǎn)折點，建立力矩平衡方程如下：

wi?sin?β（0.5hi+0.5t·tan?θr）=wi?cos?β·0.5t.?????（12）

整理上式可得處于臨空狀態(tài)的巖塊發(fā)生塊狀傾倒破壞時所需滿足的最小高寬比χmin為：

χmin?=??=??-?tan?θr?.???????（13）

滑動區(qū)和塊狀傾倒區(qū)的巖塊全部發(fā)生破壞后，對彎曲傾倒區(qū)形成反坡臨空面，其后緣巖塊處于懸臂狀態(tài).?因為下部巖塊破壞后與上部巖塊分離，層間作用力變?yōu)榱悖詮澢鷥A倒區(qū)的巖塊穩(wěn)定性問題可轉(zhuǎn)化為“斜置懸臂梁”問題，左保成[4]、盧海峰等[20]、陳從新等[21]基于圖6所示的斜置懸臂梁模型，推導(dǎo)了單巖層折斷所需滿足的最小高度hfmin?為：

hf?min?=?.???（14）

1.4???破壞類型和分析方法

本文所提反傾邊坡滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞，主要為邊坡前緣和中部巖塊先發(fā)生滑動-塊狀傾倒復(fù)合破壞，而后伴隨后緣巖層發(fā)生彎曲傾倒破壞.?從復(fù)合傾倒破壞過程來看，塊狀傾倒區(qū)的巖塊可能全部和滑動區(qū)巖塊先發(fā)生滑動-塊狀傾倒復(fù)合破壞，也可能只有部分和滑動區(qū)巖塊先發(fā)生滑動-塊狀傾倒復(fù)合破壞.

1.4.1???塊狀傾倒區(qū)巖塊全部發(fā)生復(fù)合破壞

以式（9）和（11）作為迭代方程，取Fs?=?1，從坡腳巖塊向橫向節(jié)理貫通的最后一個巖塊計算，將每一步迭代得到的min（Ti-1，Pi-1）值作為下一步迭代計算中的Ti和Pi，取Fi?=?min（Ti，Pi）（1≤i≤itotal），即邊坡失穩(wěn)破壞時所需的附加外力.?若橫向節(jié)理貫通的最后一個巖塊計算所得附加外力F

itotal?

此外滑動區(qū)與塊狀傾倒區(qū)的邊界可通過式（15）獲得[21].若巖塊ist及ist+1破壞時需要上部巖塊作用的推力滿足式（15），則[1，ist]屬于潛在滑動區(qū)，（ist，itotal]屬于潛在塊狀傾倒區(qū).

P

ist≤T

ist，

P

ist+1>T

ist+1.???????（15）

式中：P

ist和P

ist+1分別為式（9）求得巖塊ist及ist+1破壞時需要上部巖塊作用的推力;T

ist和T

ist+1分別為式（11）求得巖塊ist及ist+1破壞時需要上部巖塊作用的

推力.

邊坡安全系數(shù)Fs主要由滑動-塊狀傾倒復(fù)合破壞區(qū)域的巖塊穩(wěn)定性決定.?先假定一個初始安全系數(shù)?F（0）

s??，以式（9）和（11）作為迭代方程，從坡腳巖塊開始向該復(fù)合破壞區(qū)域的最后一個巖塊計算，求出附加外力F，直到假定的Fs滿足F為零的要求.

1.4.2???塊狀傾倒區(qū)巖塊部分發(fā)生復(fù)合破壞

以式（9）和（11）作為迭代方程，取Fs?=?1，從坡腳巖塊向橫向節(jié)理貫通的最后一個巖塊計算，將每一步迭代得到的min（Ti-1，Pi-1）值作為下一步迭代計算中的Ti和Pi，取Fi?=?min（Ti，Pi）（1≤i≤itotal），即邊坡失穩(wěn)破壞時所需的附加外力.?若橫向節(jié)理貫通的最后一個巖塊計算所得附加外力F

itotal?≥?0，前面存在附加外力F

itotal??0，則巖塊處于穩(wěn)定狀態(tài).?若巖塊計算所得的附加外力Fi?=?0，則巖塊處于臨界平衡狀態(tài).

塊狀傾倒區(qū)中，[1，im]區(qū)域巖塊發(fā)生復(fù)合破壞形成反坡臨空面，其后緣（im，itotal]區(qū)域的巖塊失去層間作用力，其穩(wěn)定性通過式（13）計算分析.?若邊坡后緣?（im，itotal]區(qū)域的塊狀傾倒區(qū)巖塊全部發(fā)生塊狀傾倒破壞，則需通過式（14）計算彎曲傾倒區(qū)巖塊穩(wěn)定性.

邊坡安全系數(shù)Fs也由滑動-塊狀傾倒復(fù)合區(qū)域的巖塊穩(wěn)定性決定.?通過先假定一個初始安全系數(shù)F（0）

s??，迭代求解坡腳巖塊至該復(fù)合破壞區(qū)域的最后一個巖塊，直到假定的Fs滿足巖塊附加外力F為零的要求.

2???工程實例數(shù)值模擬

山西省渾源縣某花崗巖礦于2017年8月23日發(fā)生復(fù)合傾倒破壞.?傾倒破壞主體位于1?950?m高程平臺下部5～15?m范圍內(nèi)，傾倒破壞區(qū)內(nèi)巖塊沿裂隙帶發(fā)生滑動-塊狀傾倒復(fù)合破壞，往正面臨空方向傾倒，后緣巖塊因正面臨空發(fā)生彎曲傾倒變形，形成“L”形的凹腔，邊坡傾倒破壞前后的概貌如圖7所示.

通過傾倒破壞區(qū)域的地質(zhì)調(diào)查可知，反傾邊坡中存在三組優(yōu)勢結(jié)構(gòu)面和底部裂隙帶：第一組結(jié)構(gòu)面J1產(chǎn)狀為67°∠73°，傾向北東，節(jié)理密度約為3.63條/m，裂隙間距的平均值為27.6?cm;與坡面形態(tài)近似垂直，解除了傾倒巖體的側(cè)向約束（因圖7為反傾邊坡復(fù)合傾倒破壞的側(cè)視圖，結(jié)構(gòu)面J1垂直于反傾邊坡坡面，故未展示J1的分布）.?第二組結(jié)構(gòu)面J2產(chǎn)狀為136°∠81°，傾向北西，節(jié)理密度約為2條/m，裂隙間距的平均值為50?cm;結(jié)構(gòu)面呈張性，延伸性較好，對傾倒破壞起主控作用.?第三組結(jié)構(gòu)面J3產(chǎn)狀為172°∠25°～45°，傾向南，節(jié)理密度約為1.5條/m，裂隙間距的平均值為65?cm;延伸性較差，對邊坡穩(wěn)定性影響較小.?反傾邊坡底部發(fā)育一條裂隙帶，順坡向傾向外，傾角30°，厚0.2?m，延伸約13?m.?結(jié)合邊坡地質(zhì)結(jié)構(gòu)，考慮陡傾結(jié)構(gòu)面J2、緩傾結(jié)構(gòu)面J3和底部裂隙帶，得到反傾邊坡地質(zhì)剖面圖（如圖8所示）.

2.1???數(shù)值模型和材料參數(shù)

反傾邊坡傾倒破壞的誘發(fā)因素是持續(xù)降雨，但是期間持續(xù)降雨的累積降雨量較小，坡體中未見明顯水流，故本文不考慮裂隙水在巖體裂隙間的動靜水壓力[22-23]，只考慮降雨對結(jié)構(gòu)面力學(xué)參數(shù)的弱化作用.?因此，數(shù)值模擬主要包括天然工況和降雨工況，采用UDEC的摩爾庫侖模型，計算所采用的材料參數(shù)通過室內(nèi)試驗、工程類比和參數(shù)反演綜合得出[24]（見表11），天然工況采用天然材料參數(shù)，降雨工況采用飽和材料參數(shù).?基于圖8所示的反傾邊坡地質(zhì)剖面圖，建立了反傾邊坡的數(shù)值模型（如圖9所示）.

邊坡模型高25?m，長35?m;坡高15?m，坡度65°;內(nèi)部發(fā)育兩組結(jié)構(gòu)面：陡傾結(jié)構(gòu)面J2傾角81°、間距1?m，緩傾結(jié)構(gòu)面J3傾角25°、間距0.8m;底部發(fā)育一條裂隙帶，傾角30°，延伸約13?m.

2.2???天然工況模擬

天然工況下數(shù)值計算收斂，總共計算229?867步，最終邊坡位移矢量圖和坡面位移量如圖10所示.?邊坡整體位移變形較小，說明該邊坡在天然狀態(tài)下是穩(wěn)定的.?從圖10中具體來看：邊坡坡腳巖塊位移很小，底部裂隙面閉合且無滑動;坡頂巖塊在自重作用下，有水平向外的傾倒趨勢，但傾倒位移量較小，坡面位移量隨坡高的增加近似呈線性增長，最大位移量4.4×10-3?m，位于靠近坡肩的位置，邊坡整體呈穩(wěn)定狀態(tài).

2.3???降雨工況模擬

雨水的弱化作用導(dǎo)致邊坡巖體和結(jié)構(gòu)面的力學(xué)參數(shù)降低，在飽和力學(xué)參數(shù)條件下，數(shù)值計算不收斂，說明該邊坡在降雨狀態(tài)下是失穩(wěn)的.

圖11所示為降雨工況下，邊坡計算50萬步時的位移矢量圖和坡面位移量.

從圖11中可以看出，降雨工況下邊坡整體位移變形較大，坡面位移量隨坡高不再呈線性增長，而是呈現(xiàn)明顯的分區(qū)特征，具體可分為滑動區(qū)、塊狀傾倒區(qū)和彎曲傾倒區(qū)，巖塊最大位移7.5×10-2?m，位于塊狀傾倒區(qū)后緣頂部的位置.?根據(jù)邊坡的位移矢量特征（箭頭方向表1示位移方向，箭頭長度表1示位移大?。┛蛇M(jìn)一步得出三個分區(qū)的巖塊位移特征：滑動區(qū)集中在坡腳位置，主要為前三個巖塊，位移方向為平行于底部裂隙面向下，巖塊各高度位移大小近似相等，呈整體向下剪切滑動的特征.?塊狀傾倒區(qū)主要位于滑動區(qū)上方，位移方向近似水平向外，巖塊頂部位移大，底部位移小，位移矢量近似垂直于反傾結(jié)構(gòu)面，呈現(xiàn)沿底部裂隙面下端點向外轉(zhuǎn)動傾倒的特征.?彎曲傾倒區(qū)巖塊底部沒有貫通的裂隙帶，發(fā)生彎曲傾倒變形，處于彎曲變形狀態(tài).?因此復(fù)合傾倒主要破壞為塊狀傾倒區(qū)巖塊先發(fā)生傾倒破壞，擠壓推動前部巖塊，觸發(fā)前緣巖塊發(fā)生滑動破壞.

圖12所示為降雨工況下，邊坡計算50萬步時的塑性區(qū)和局部破壞圖.?彎曲傾倒區(qū)位于塊狀傾倒區(qū)的上方，巖塊底部后側(cè)出現(xiàn)局部拉伸屈服破壞，前側(cè)出現(xiàn)局部剪切破壞（圖12中小圓圈表1示拉伸破壞塑性區(qū)，星號表1示剪切破壞塑性區(qū)），呈現(xiàn)向外彎曲傾倒的特征.?因彎曲傾倒區(qū)的巖塊彎曲拉應(yīng)力未達(dá)到其抗拉強(qiáng)度，故彎曲傾倒區(qū)的巖塊未出現(xiàn)傾倒破壞，處于彎曲變形狀態(tài).?因此，降雨工況下邊坡先發(fā)生前緣和中部滑動-塊狀傾倒破壞，并伴隨后緣巖塊彎曲傾倒變形.

此外，從滑動區(qū)最后一個巖塊失穩(wěn)的局部放大圖（A1）可看出滑動區(qū)巖塊的破壞特征：該巖塊底部與裂隙帶（破壞面）壓密閉合，相互剪切錯動，發(fā)生壓剪滑移破壞;巖塊后部出現(xiàn)圖中所示的明顯張裂隙，張裂隙呈現(xiàn)頂部窄底部寬的形態(tài)，這主要是由于后部巖塊沿底部下端點發(fā)生轉(zhuǎn)動傾倒，上部變形大于下部，上部擠推前部巖塊，致使滑動區(qū)后緣巖塊的下部形成“張裂隙”，此張裂隙可作為滑動區(qū)和塊狀傾倒區(qū)的分界面.?從塊狀傾倒區(qū)的最后一個巖塊失穩(wěn)的局部放大圖（A2）可看出：塊狀傾倒區(qū)的巖塊在底部和后部存在兩條張裂隙，底部張裂隙沿破壞面呈現(xiàn)前窄后寬的形態(tài)，后部張裂隙沿巖塊高度呈現(xiàn)頂寬底窄的形態(tài)，表1現(xiàn)出明顯的塊狀傾倒特征.?因此，從巖塊失穩(wěn)的局部放大圖進(jìn)一步肯定了該邊坡的滑動-塊狀傾倒復(fù)合破壞特征.

圖11和圖12的UDEC模擬結(jié)果均一致反映了底部橫向節(jié)理未貫通的反傾邊坡的復(fù)合破壞機(jī)制，表1現(xiàn)為前緣和中部滑動-塊狀傾倒破壞，而后伴隨后緣巖塊彎曲傾倒變形.

為便于理論公式的計算分析，提出了反傾邊坡簡化地質(zhì)模型，即不考慮緩傾結(jié)構(gòu)面J3的影響.?對該簡化地質(zhì)模型采用相同的飽和材料參數(shù)，計算50萬步可得降雨工況下邊坡位移矢量和局部破壞如圖13所示.

從圖13中簡化地質(zhì)模型邊坡的位移矢量來看，依然呈現(xiàn)滑動區(qū)、塊狀傾倒區(qū)和彎曲傾倒區(qū)三個分區(qū)特征，滑動區(qū)和塊狀傾倒區(qū)的位移矢量特征和原始地質(zhì)模型完全一致，彎曲傾倒區(qū)由于不考慮緩傾結(jié)構(gòu)面J3，巖塊的變形比原始地質(zhì)模型稍微偏小.?總體來看，簡化地質(zhì)模型的邊坡位移矢量特征和塑性區(qū)分布也表1現(xiàn)出滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞機(jī)制，這和圖11與圖12原始邊坡地質(zhì)模型所得結(jié)果一致.?此外，對比分析滑動區(qū)巖塊的局部失穩(wěn)圖A1′和A1，后部張裂隙形狀均呈現(xiàn)壓剪滑移破壞特征;塊狀傾倒區(qū)的局部失穩(wěn)圖A2′和A2，底部和后部兩條張裂隙形狀均一致呈現(xiàn)出塊狀傾倒特征.?因此從邊坡位移矢量和局部破壞特征可知該反傾邊坡中的緩傾結(jié)構(gòu)面J3對邊坡的破壞無顯著影響，采用簡化地質(zhì)模型研究邊坡破壞機(jī)制和分析方法是合理的.

3???對比分析

文中實例邊坡的復(fù)合傾倒破壞機(jī)制，通過數(shù)值模擬具體表1現(xiàn)為邊坡前緣和中部巖塊發(fā)生滑動-塊狀傾倒破壞，后緣巖塊發(fā)生彎曲傾倒變形.?在表11所示的材料參數(shù)條件下采用本文所提復(fù)合傾倒破壞分析方法對該反傾邊坡簡化地質(zhì)模型進(jìn)行計算分析，計算模型如圖14所示.

按照1.4節(jié)的分析步驟，首先取Fs?=?1，以式（9）和（11）作為迭代方程，從坡腳巖塊（i=1）向橫向節(jié)理貫通的最后一個巖塊（i=8）計算，將每一步迭代得到的min（Ti-1，Pi-1）值作為下一步迭代計算中的Ti和Pi，取Fi?=?min（Ti，Pi）（1≤i≤8），即邊坡失穩(wěn)破壞時所需的附加外力.?在飽和材料參數(shù)條件下，滑動-塊狀傾倒破壞區(qū)域的各巖塊發(fā)生滑動破壞或塊狀傾倒破壞所需附加外力Pi和Ti，見表12.?從表12中各巖塊滑動破壞或塊狀傾倒破壞時所需附加外力來看，依據(jù)較小附加外力所對應(yīng)的類型即為巖塊的潛在破壞類型，可清楚地判定各巖塊的破壞類型，如表12所示.?此外，第8號坡頂巖塊所需附加外力小于零，處于失穩(wěn)狀態(tài)，頂部巖塊傾倒加壓于下部塊體是導(dǎo)致邊坡發(fā)生滑動-塊狀傾倒復(fù)合破壞的根本原因，故對該類邊坡治理加固時應(yīng)及早采取措施，重點加固塊狀傾倒破壞區(qū).

進(jìn)一步對理論方法和數(shù)值模擬的計算結(jié)果作對比分析，如表13所示.?從表13中可看出，利用本文所提的理論方法和數(shù)值模擬得到的結(jié)果近似一致，兩種方法計算的邊坡穩(wěn)定性結(jié)果均處于失穩(wěn)狀態(tài)，兩種方法確定的滑動區(qū)和塊狀傾倒區(qū)的范圍基本相同，均與實際相符.?此外，對于彎曲傾倒區(qū)的巖塊，通過式（14）計算得巖塊最小折斷高度為24?m，遠(yuǎn)大于該邊坡彎曲傾倒區(qū)的巖塊懸臂長度10?m，因此巖塊僅處于彎曲變形狀態(tài)，并未發(fā)生彎曲傾倒破壞，這也與數(shù)值模擬的結(jié)果一致.

邊坡安全系數(shù)Fs以式（9）和（11）作為迭代方程，從坡腳巖塊開始向該復(fù)合破壞區(qū)域的最后一個巖塊計算，求出附加外力F，直到假定的Fs滿足F為零的要求.?對該實例邊坡采用天然材料參數(shù)計算得安全系數(shù)Fs為1.35，F(xiàn)s?>1說明邊坡處于穩(wěn)定狀態(tài);采用飽和材料參數(shù)計算得安全系數(shù)Fs為0.96，F(xiàn)s?<1說明邊坡處于失穩(wěn)狀態(tài).?從邊坡安全系數(shù)角度考慮，理論方法和數(shù)值模擬的結(jié)果亦一致.

4???結(jié)???論

1）基于底部橫向裂隙未完全貫通的反傾邊坡，建立了包含滑動區(qū)、塊狀傾倒區(qū)和彎曲傾倒區(qū)的復(fù)合傾倒破壞地質(zhì)模型.

2）利用三個分區(qū)的巖塊破壞特點，基于極限平衡理論和懸臂梁模型，提出了巖質(zhì)反傾邊坡復(fù)合傾倒破壞自下而上逐步分析的理論方法.

3）巖質(zhì)反傾邊坡滑動-塊狀傾倒-彎曲傾倒復(fù)合破壞的穩(wěn)定性由滑動-塊狀傾倒復(fù)合破壞區(qū)域控制.

4）滑動-塊狀傾倒復(fù)合破壞區(qū)中塊狀傾倒區(qū)域

屬于主動破壞區(qū)域，滑動區(qū)屬于被動破壞區(qū)域，治理加固時應(yīng)重點加固塊狀傾倒破壞區(qū).

參考文獻(xiàn)

[1]????ALEJANO?L?R，G？MEZ-MáRQUEZ?I，MARTíNEZ-ALEGRíA?R.?Analysis?of?a?complex?toppling-circular?slope?failure?[J].?Engineering?Geology，2010，114（1/2）：93—104.

[2]????邱俊，任光明，王云南.?層狀反傾-順傾邊坡傾倒變形形成條件及發(fā)育規(guī)模特征[J].?巖土力學(xué)，2016，37（S2）：513—524.

QIU?J，REN?G?M，WANG?Y?N.?Characteristics?of?forming?conditions?and?development?scale?of?toppling?in?anti-dip?and?dip?stratified?slopes?[J].?Rock?and?Soil?Mechanics，2016，37（S2）：513—524.?（In?Chinese）

[3]????劉順昌.?如美水電站巖質(zhì)邊坡傾倒破壞機(jī)理研究[D].武漢：中國地質(zhì)大學(xué)工程學(xué)院，2013：7—14.

LIU?S?C.?Study?on?toppling?failure?mechanism?of?rock?slope?in?Rumei?hydropower?station?[D].?Wuhan：Faculty?of?Engineering，China?University?of?Geosciences，2013：7—14.?（In?Chinese）

[4]????左保成.?反傾巖質(zhì)邊坡破壞機(jī)理研究[D].?武漢：中國科學(xué)院武漢巖土力學(xué)研究所，2004：7—13.

ZUO?B?C.?Study?on?losing?stability?mechanics?of?counter-tilt?rock?slopes?[D].?Wuhan：Institute?of?Rock?and?Soil?Mechanics，Chinese?Academy?of?Sciences，2004：7—13.?（In?Chinese）

[5]????盧海峰，劉泉聲，陳從新.?反傾巖質(zhì)邊坡懸臂梁極限平衡模型的改進(jìn)[J].?巖土力學(xué)，2012，33（2）：577—584.

LU?H?F，LIU?Q?S，CHEN?C?X.?Improvement?of?cantilever?beam?limit?equilibrium?model?of?counter-tilt?rock?slopes?[J].?Rock?and?Soil?Mechanics，2012，33（2）：577—584.?（In?Chinese）

[6]????譚強(qiáng)，高明忠，魏進(jìn)兵，等.?巖質(zhì)邊坡淺表1層改造蠕動破壞機(jī)理探索與穩(wěn)定性分析[J].?湖南大學(xué)學(xué)報（自然科學(xué)版），2014，41（11）：103—109.

TAN?Q，GAO?M?Z，WEI?J?B，et?al.?A?case?study?on?creep?destruction?resulted?from?exogenetic?reconstruction?and?stability?evaluation?of?rock?slope?[J].?Journal?of?Hunan?University?（Natural?Sciences），2014，41（11）：103—109.?（In?Chinese）

[7]????GOODMAN?R?E，BRAY?J?W.?Toppling?of?rock?slopes[C]//Proceedings?of?ASCE?Specialty?Conference，Rock?Engineering?for?Foundations?and?Slopes：Vol.2.?Colorado：Boulder，1976：201—234.

[8]????ADHIKARY?D?P，DYSKIN?A?V，JEWELL?R?J，et?al.?A?study?of?the?mechanism?of?flexural?toppling?failure?of?rock?slopes?[J].?Rock?Mechanics?&?Rock?Engineering，1997，30（2）：75—93.

[9]????ZHENG?Y，CHEN?C?X，LIU?T?T，et?al.?Stability?analysis?of?rock?slopes?against?sliding?or?flexural-toppling?failure?[J].?Bulletin?of?Engineering?Geology?&?the?Environment，2017（1）：1—21.

[10]??BOBET?A.?Analytical?solutions?for?toppling?failure?[J].?International?Journal?of?Rock?Mechanics?and?Mining?Sciences，1999，36（7）：971—980.

[11]??LIU?C?H，JAKSA?M?B，MEYERS?A?G.?Toppling?mechanisms?of?rock?slopes?considering?stabilization?from?the?underlying?rock?mass?[J].?International?Journal?of?Rock?Mechanics?and?Mining?Sciences，2008，45（8）：348—354.

[12]??LIU?C?H，JAKSA?M?B，MEYERS?A?G.?Improved?analytical?solution?for?toppling?stability?analysis?of?rock?slopes?[J].?International?Journal?of?Rock?Mechanics?and?Mining?Sciences，2008，45（8）：1361—1372.

[13]??LIU?C?H，JAKSA?M?B，MEYERS?A?G.?A?transfer?coefficient?method?for?rock?slope?toppling?[J].?Canadian?Geotechnical?Journal，2009，46（1）：1—9.

[14]??劉才華，陳從新.?地震作用下巖質(zhì)邊坡塊體傾倒破壞分析[J].?巖石力學(xué)與工程學(xué)報，2010，29（S1）：3193—3198.

LIU?C?H，CHEN?C?X.?Analysis?of?toppling?failure?of?rock?slopes?due?to?earthquakes?[J].?Chinese?Journal?of?Rock?Mechanics?and?Engineering，2010，29（S1）：3193—3198.?（In?Chinese）

[15]??鄭允，陳從新，劉婷婷，等.?坡頂荷載作用下巖質(zhì)邊坡傾倒破壞分析[J].?巖土力學(xué)，2015，36（9）：2639—2647.

ZHENG?Y，CHEN?C?X，LIU?T?T，et?al.?Analysis?of?toppling?failure?of?rock?slopes?under?the?loads?applied?on?the?top?[J].?Rock?and?Soil?Mechanics，2015，36（9）：2639—2647.?（In?Chinese）

[16]??AMINI?M，VESHADI?M?A.?Stability?analysis?of?rock?slopes?against?block-flexure?toppling?failure?[J].?Rock?Mechanics?&?Rock?Engineering，2012，45（4）：519—532.

[17]??AMINI?M，GHOLAMZADEH?M，KHOSRAVI?M?H.?Physical?and?theoretical?modeling?of?rock?slopes?against?block-flexure?toppling?failure?[J].?International?Journal?of?Mining?and?Geo-Engineering，2015，49（2）：155—171.

[18]??MOHTARAMI?E，JAFARI?A，AMINI?M.?Stability?analysis?of?slopes?against?combined?circular-toppling?failure?[J].?International?Journal?of?Rock?Mechanics?&?Mining?Sciences，2014，67（2）：43—56.

[19]??AMINI?M，ARDESTANI?A，KHOSRAVI?M?H.?Stability?analysis?of?slide-toe-toppling?failure?[J].?Engineering?Geology，2017，228—96.

[20]??盧海峰，劉泉聲，陳從新.?反傾巖質(zhì)邊坡懸臂梁極限平衡模型的改進(jìn)[J].?巖土力學(xué)，2012，33（2）：577—584.

LU?H?F，LIU?Q?S，CHEN?C?X.?Improvement?of?cantilever?beam?limit?equilibrium?model?of?counter-tilt?rock?slopes?[J].?Rock?and?Soil?Mechanics，2012，33（2）：577—584.?（In?Chinese）

[21]??陳從新，鄭允，孫朝燚.?巖質(zhì)反傾邊坡彎曲傾倒破壞分析方法研究[J].?巖石力學(xué)與工程學(xué)報，2016，35（11）：2174—2187.

CHEN?C?X，ZHENG?Y，SUN?C?Y.?Analytical?approach?on?flexural?toppling?failure?of?counter-tilt?slopes?of?layered?rock?[J].?Chinese?Journal?of?Rock?Mechanics?and?Engineering，2016，35（11）：2174—2187.?（In?Chinese）

[22]??舒繼森，唐震，才慶祥.?水力學(xué)作用下順層巖質(zhì)邊坡穩(wěn)定性研究[J].?中國礦業(yè)大學(xué)學(xué)報，2012，41（4）：521—525.

SHU?J?S，TANG?Z，CAI?Q?X.?Research?on?stability?of?bedding?rock?slopes?under?hydraulic?pressure?[J].?Journal?of?China?University?of?Ming?&?Technology，2012，41（4）：521—525.?（In?Chinese）

[23]??鄭俊杰，魚歡，郭震山，等.?降雨條件下拓寬路堤穩(wěn)定性數(shù)值分析[J].?湖南大學(xué)學(xué)報（自然科學(xué)版），2017，44（7）：150—155.

ZHENG?J?J，YU?H，GUO?Z?S，et?al.?Numerical?analysis?of?the?stability?of?widened?embankment?under?rainfall?conditions?[J].?Journal?of?Hunan?University?（Natural?sciences），2017，44（7）：150—155.?（In?Chinese）

[24]??陳從新，肖國峰，魯祖德，等.?新利花崗巖礦邊坡穩(wěn)定性研究[R].?武漢：中國科學(xué)院武漢巖土力學(xué)研究所，2015：41—92.

CHEN?C?X，XIAO?G?F，LU?Z?D，et?al.??Report?on?the?stability?evaluation?of?Xinli?granite?mine?slope[R].?Wuhan：Institute?of?Rock?and?Soil?Mechanics，Chinese?Academy?of?Sciences，2015：41—92.?（In?Chinese）