貫通型單臺階巖體結(jié)構(gòu)面剪切性質(zhì)及應(yīng)用

黃達(dá)，黃潤秋，雷鵬

?

貫通型單臺階巖體結(jié)構(gòu)面剪切性質(zhì)及應(yīng)用

黃達(dá)1, 2，黃潤秋3，雷鵬1

(1. 重慶大學(xué)土木工程學(xué)院，重慶，400045；2. 重慶大學(xué)煤礦災(zāi)害動力學(xué)與控制國家重點實驗室，重慶，400044；3. 成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國家重點實驗室，四川成都，610059)

利用二維顆粒流(PFC2D)方法模擬研究臺階高長比(/)對貫通型單臺階巖體結(jié)構(gòu)面剪切變形及強(qiáng)度影響。研究結(jié)果表明：臺階損傷破壞模式可分為壓切和張剪破壞2種基本類型，隨著/增大，壓切模式逐漸轉(zhuǎn)化為張剪模式，但其與法向應(yīng)力基本無關(guān)。當(dāng)臺階張剪損傷破壞時(＞0.45)，受彎矩及剪切力偶合作用，臺階根部易張裂；當(dāng)臺階壓切損傷破壞時(≤0.45)，主要受斜向上的剪切力作用而使其尖角剪斷，沿剪切面呈現(xiàn)明顯的爬坡效應(yīng)，具剪脹特性。針對實際巖體結(jié)構(gòu)面臺階高長比/常較小(壓切模式)，建立同時考慮結(jié)構(gòu)面爬坡與啃斷效應(yīng)的剪切強(qiáng)度經(jīng)驗公式，并運用于邊坡穩(wěn)定性計算的參數(shù)取值。

巖體力學(xué)；結(jié)構(gòu)面；剪切變形；剪切強(qiáng)度；顆粒流

巖體結(jié)構(gòu)面地質(zhì)力學(xué)性質(zhì)對巖體穩(wěn)定性及加固設(shè)計均有非常重要的影響。巖體結(jié)構(gòu)面力學(xué)性質(zhì)與壁巖特性(壁巖成分、強(qiáng)度、結(jié)構(gòu)構(gòu)造等)、結(jié)構(gòu)面結(jié)合狀態(tài)(充填度、充填物成分、充填膠結(jié)強(qiáng)度)和結(jié)構(gòu)面起伏形態(tài)有關(guān)[1?6]。尤其是無充填型結(jié)構(gòu)面表面起伏形態(tài)對其剪切變形、損傷破壞及強(qiáng)度等性質(zhì)影響顯著[5?6]。結(jié)構(gòu)面粗糙起伏一般可分為平直結(jié)構(gòu)面、鋸齒狀結(jié)構(gòu)面、波狀起伏結(jié)構(gòu)面、臺階狀結(jié)構(gòu)面這4種基本類 型[6]。目前主要集中于諸如規(guī)則鋸齒狀及波狀起伏型結(jié)構(gòu)面的破損及強(qiáng)度等力學(xué)效應(yīng)研究。研究發(fā)現(xiàn)，當(dāng)法向應(yīng)力相對較低時，結(jié)構(gòu)面在剪切時沿突齒斜面滑移爬坡，出現(xiàn)剪脹現(xiàn)象；而當(dāng)法向應(yīng)力超過某一臨界值時，傾斜粗糙面上的滑動被抑制，將剪切啃斷鋸 齒[7]。HOMAND等[8]基于低法向應(yīng)力(小于5 MPa)、小截面尺寸(長150 mm×寬150 mm)和較小齒狀起伏的花崗巖結(jié)構(gòu)面循環(huán)直剪試驗，分析了剪切前后鋸齒結(jié)構(gòu)面表面面積的衰減變形規(guī)律，定量描述了低法向應(yīng)力循環(huán)剪切條件下鋸齒狀結(jié)構(gòu)面磨損隨法向應(yīng)力增強(qiáng)的特性。周輝等[9]研究了水泥砂漿鋸齒狀結(jié)構(gòu)面在不同起伏高度、不同剪切速率和不同法向壓力條件下聲發(fā)射參數(shù)的變化規(guī)律和發(fā)生機(jī)制。李海波等[10?12]采用混凝土試件開展了鋸齒型結(jié)構(gòu)面剪切強(qiáng)度及變形等力學(xué)特性直剪試驗研究。目前，針對具凸起臺階狀結(jié)構(gòu)面的剪切性質(zhì)方面的研究相對較少。郭志[13]采用石膏試件開展了臺階高為1.0，0.5，1.0和1.5 cm且高長比 (/)為0.06，0.16，0.30和0.50臺階狀結(jié)構(gòu)面的直剪試驗，發(fā)現(xiàn)峰值強(qiáng)度與殘余強(qiáng)度之差隨著高長比/的增大而增大。但試驗試件太小，尺寸效應(yīng)明顯，且由于試驗設(shè)備等限制而沒有揭示結(jié)構(gòu)面損傷破裂演化。利用顆粒流程序(particle flow code，PFC)開展巖石力學(xué)試驗其主要優(yōu)點在于避免了設(shè)定材料宏觀本構(gòu)的經(jīng)驗主觀行為，而是通過顆粒接觸黏結(jié)機(jī)理從細(xì)觀尺度仿真巖石礦物顆粒及其間的黏結(jié)與摩擦來實現(xiàn)宏觀材料的組構(gòu)，反映了顆粒膠結(jié)型巖石類材料的本質(zhì)細(xì)觀結(jié)構(gòu)屬性[14]。POTYONDY等[15?17]采用PFC模擬了多種巖石宏細(xì)觀力學(xué)特性，較充分地闡述或論證了PFC模擬在研究巖石裂隙擴(kuò)展、細(xì)觀損傷演化及力學(xué)性質(zhì)等方面的科學(xué)優(yōu)勢。為了豐富對臺階狀結(jié)構(gòu)面剪切特性的研究，本文作者首先從理論上探討臺階狀結(jié)構(gòu)面剪切特性，進(jìn)而利用顆粒流程序開展較大長度且具有不同高長比(/)單臺階巖體結(jié)構(gòu)面試件直剪數(shù)值試驗。從顆粒細(xì)觀位移、裂紋演化及剪切應(yīng)力?位移曲線等方面研究臺階狀結(jié)構(gòu)面剪切變形的細(xì)宏觀多尺度演化過程；建立綜合考慮結(jié)構(gòu)面基本摩擦角、巖石黏聚力b及臺階高長比/的臺階狀結(jié)構(gòu)面剪切強(qiáng)度經(jīng)驗公式。

1 結(jié)構(gòu)面凸起臺階壓剪破裂理論 分析

具有臺階狀起伏形態(tài)的外傾結(jié)構(gòu)面在斜坡(坡體常為層狀或塊裂結(jié)構(gòu))的淺表生改造帶常見，特別是在反傾層狀巖質(zhì)斜坡的淺表層，順層斜坡也偶有出現(xiàn)。斜/邊坡穩(wěn)定狀態(tài)明顯受外傾貫通型結(jié)構(gòu)面的凸起巖石臺階所控制性，即受巖石臺階規(guī)模及強(qiáng)度等綜合控制。根據(jù)一般的固體力學(xué)理論，可得到2點單臺階結(jié)構(gòu)面剪切性質(zhì)的定性認(rèn)知：1) 當(dāng)臺階高遠(yuǎn)大于臺階長時，臺階如同懸臂梁一樣，臺階根部既有剪力作用又有彎矩作用，在這種復(fù)合應(yīng)力的作用下，臺階根部巖石更易于張剪復(fù)合破壞甚至彎曲折斷；2) 當(dāng)臺階高遠(yuǎn)小于臺階長時，臺階巖石主要受剪切力作用而彎矩很小，結(jié)構(gòu)面易剪斷臺階尖角滑移。在絕大多數(shù)情況下，在凸起臺階前方(剪切方向)巖體結(jié)構(gòu)面的延伸長度要比凸起臺階的高度大得多，故真實巖體的臺階型結(jié)構(gòu)面常滿足遠(yuǎn)小于。

若臺階長度足夠長且假設(shè)剪裂面通過臺階腳點，則建立如圖1所示的結(jié)構(gòu)面凸起臺階剪切破裂模型。圖1中：虛斜線為假設(shè)的剪切破裂面；為臺階高；為臺階長；n和分別為結(jié)構(gòu)面上法向應(yīng)力和剪切應(yīng)力；′n和′分別為作用在斜面上的正應(yīng)力和切應(yīng)力；為臺階剪切破裂面傾角。破裂面上正應(yīng)力′n和剪應(yīng)力′分別為：

圖1 凸起巖石臺階剪切破裂模型

Fig. 1 Shear fracture model for a raised rock step

由庫侖準(zhǔn)則可知臺階剪裂面上抗剪應(yīng)力′r為

式中：b和分別為完整巖石的黏聚力和摩擦角。

(4)

此時，臺階的剪切應(yīng)力為

式(5)對求導(dǎo)，則式(5)在極值條件下可求得斜面傾角=45°?b/2，再代入式(5)，得臺階的臨界剪切應(yīng)力為

由于凸起臺階最不利剪切破裂面與水平面的夾角=45°?b/2，而此時臺階高和長需滿足＜tan(45°?b/2)。

2 臺階狀結(jié)構(gòu)面直剪試驗PFC模擬

2.1 數(shù)值模型及參數(shù)標(biāo)定

圖2所示為臺階狀結(jié)構(gòu)面PFC2D數(shù)值模型，其中，正方形模型邊長為0.4 m，貫通閉合型單臺階巖體結(jié)構(gòu)面兩側(cè)顆粒用深色示意。定義6道墻作為邊界條件(側(cè)墻5，6和2，3均以結(jié)構(gòu)面為分界點)，對構(gòu)建下部剪切盒的邊界墻(如圖2中1，2和6墻)進(jìn)行法向變形約束，對上部剪切盒兩側(cè)墻(如圖2中3和5墻)賦予一定的同向等量移動速度(剪切速率)即可對試樣施加水平剪切推力(圖2中箭頭表示墻運動方向)，恒定的法向荷載可通過編制伺服系統(tǒng)程序不斷地調(diào)整上部墻(如圖2中4墻)的位移速度實現(xiàn)。

圖2 臺階狀結(jié)構(gòu)面PFC2D數(shù)值模型

模型臺階長均設(shè)定為0.2 m，臺階高長比/設(shè)定為0.10，0.20，0.30，0.40，0.45，0.50和0.60，法向荷載n設(shè)為3.5，5.0，6.5，8.0和12.0 MPa，剪切速率為0.005 mm/s。選用較符合巖石細(xì)觀顆粒間組構(gòu)的平行黏結(jié)接觸模型(parallel-bonded model)，其不但可模擬巖石細(xì)觀顆粒間的法向和切向受力方式，而且可模擬顆粒間因膠結(jié)形成的力矩[15]。

通過單軸壓縮、直剪數(shù)值試驗測試巖體宏觀參數(shù)與室內(nèi)試驗結(jié)果相匹配來實現(xiàn)細(xì)觀參數(shù)的標(biāo)定。較簡單有效的標(biāo)定步驟為：1) 將平行黏結(jié)法向強(qiáng)度b,m和切向強(qiáng)度b,m設(shè)為較大的值，通過改變顆粒接觸模量c和平行黏結(jié)模量來匹配巖石宏觀彈性模量；2) 改變顆粒剛度比n/s和平行黏結(jié)剛度比來匹配巖石泊松比；3) 改變b,m和b,m來匹配單軸壓縮強(qiáng)度；4) 改變顆粒摩擦因數(shù)等細(xì)觀參數(shù)來匹配巖石內(nèi)摩擦角和黏聚力。選用重慶地區(qū)常發(fā)育臺階狀裂隙的三疊系下統(tǒng)嘉陵江組白云巖作為原型，經(jīng)過顆粒間細(xì)觀逼近宏觀力學(xué)參數(shù)的目標(biāo)調(diào)試，最終確定一組較可靠的平行黏結(jié)模型細(xì)觀參數(shù)，如表1所示。表2所示為PFC數(shù)值試驗宏觀參數(shù)與室內(nèi)試驗結(jié)果的比較，圖3所示為在5 MPa法向應(yīng)力作用下白云巖直剪試驗與PFC直剪模擬剪應(yīng)力?剪切位移曲線圖。由表2和圖3可知：表1中的巖石的細(xì)觀參數(shù)較真實地反映了試驗白云巖的細(xì)觀組構(gòu)及力學(xué)性質(zhì)。平直近光滑的貫通型白云巖結(jié)構(gòu)面的基本摩擦角r為35.5°。

表1 PFC模型細(xì)觀參數(shù)

Table 1 Meso-mechanical parameters of PFC model

表2 室內(nèi)試驗測得試樣宏觀參數(shù)與PFC模擬結(jié)果比較

Table 2 Comparison of macro properties for the specimens tested in lab and by PFC simulation

1—直剪試驗；2—PFC模擬。

2.2 臺階損傷破壞模式及演化

作為一種離散元方法，顆粒流程序能從細(xì)宏觀多尺度模擬結(jié)構(gòu)面的剪切變形和破壞。基于裂紋發(fā)展、分布及位移矢量的分析表明，直立臺階的損傷破壞模式主要受其高長比的影響，而法向應(yīng)力的影響相對很少。按臺階高長比的不同，變形破壞模式可分為2類：當(dāng)臺階高長比＞0.45時，凸起臺階將張剪破壞；當(dāng)臺階高長比≤0.45時，凸起臺階將壓切破壞。與第2節(jié)理論分析凸起臺階沿最不利剪切破裂面壓切破壞時需滿足＜tan(45°?b/2)=0.479基本吻合，顯示了顆粒流方法的可靠性。

2.2.1 張剪損傷破壞

圖4所示為結(jié)構(gòu)面臺階高長比=0.5、法向應(yīng)力n=5 MPa時剪切過程中剪應(yīng)力及累積微裂紋數(shù)隨剪切變形的演化曲線。根據(jù)剪應(yīng)力?位移及微裂紋?位移曲線，將剪切過程分為3個階段：彈性階段(圖4中點之前)，峰前累進(jìn)損傷階段(圖4中段)，峰后破壞階段(圖4中點之后)。峰前累積損傷段剪應(yīng)力?位移曲線呈小幅多頻性波動，微裂紋穩(wěn)步增加；在峰值點后應(yīng)力跌落幅度較大，跌落過程中微裂紋數(shù)顯著增加。

(a) 剪應(yīng)力?位移曲線；(b) 微裂紋?位移曲線▲—剪切狀態(tài)1；●—剪切狀態(tài)2；■—剪切狀態(tài)3

圖5所示為圖4中3個不同剪切狀態(tài)時刻結(jié)構(gòu)面凸起臺階裂紋擴(kuò)展及位移矢量。臺階上角點(上尖角)類似于點荷載頂壓及剪切推力的綜合作用下，致使上角點左上方(后上方)巖體首先出現(xiàn)張剪復(fù)合裂紋(圖5(a))。由于臺階高長比相對較大，隨著剪切力的增加，凸起臺階處的位移矢量表明臺階既受彎矩作用又受剪力作用，故臺階根部易于在彎矩作用下張裂而形成張拉裂紋(圖5(b)~(c))。兩角點附近的宏觀張性裂紋傾伏向下且傾伏角小于25°(圖5(c))，隨高長比減小而逐漸減小，當(dāng)=0.45時幾乎為0°(近水平)。

(a) 剪切狀態(tài)1；(b) 剪切狀態(tài)2；(c) 剪切狀態(tài)3

2.2.2 壓切損傷破壞

圖6所示為結(jié)構(gòu)面臺階高長比=0.4、法向應(yīng)力n=5 MPa時剪切過程中剪應(yīng)力及累積微裂紋數(shù)隨剪切變形的演化曲線。相對于圖4中的張剪模式，壓切模式存在2點明顯區(qū)別：峰后應(yīng)力跌落幅度(即峰值與殘余應(yīng)力差)相對較小；累進(jìn)損傷及峰后破壞過程中產(chǎn)生相對更多的微裂紋。

(a) 剪應(yīng)力?位移曲線；(b) 微裂紋?位移曲線

圖7所示為圖6中3個不同剪切狀態(tài)時刻結(jié)構(gòu)面凸起臺階裂紋擴(kuò)展及位移矢量。凸起臺階處的位移矢量表明臺階主要受斜向上的剪切力作用，這就是水平剪切推力作用下使相對較小的臺階尖角剪斷的變形響應(yīng)(沿如圖7(c)所示的剪裂面剪斷)，致使臺階根部及上尖角附近應(yīng)力集中，臺階根部及上尖角附近巖體相繼壓剪損傷，造成此兩處微裂紋發(fā)育且宏觀上呈壓碎狀(如圖7(b))。壓切破壞模式沿剪切面呈現(xiàn)明顯的爬坡效應(yīng)，結(jié)構(gòu)面出現(xiàn)張裂架空(見圖7(c))，具剪脹性質(zhì)。

(a) 剪切狀態(tài)1；(b) 剪切狀態(tài)2；(c) 剪切狀態(tài)3

3 結(jié)構(gòu)面剪切強(qiáng)度

結(jié)構(gòu)面抗剪強(qiáng)度一般通過Mhor?Coulomb屈服準(zhǔn)則確定。由前述的討論可知，當(dāng)臺階高長比> 0.45時，臺階是在彎矩和剪力聯(lián)合作用下張剪破壞，顯然這時利用Mhor?Coulomb準(zhǔn)則不恰當(dāng)。而實際結(jié)構(gòu)面臺階的高長比均比較小，故這里只探討當(dāng)＜0.45的情況。

為了使臺階狀巖體結(jié)構(gòu)面強(qiáng)度公式的物理意義更明確，構(gòu)建其剪切強(qiáng)度公式為

4 邊坡穩(wěn)定性算例

假設(shè)順向貫通型單臺階結(jié)構(gòu)面控制的某白云巖邊坡如圖9所示，坡高為16 m，坡角為65°，離坡頂左側(cè)10 m。固定結(jié)構(gòu)面點的位置不變，調(diào)整點位置，設(shè)置臺階高長比/分別為0.1，0.2，0.3和0.4。結(jié)構(gòu)面基本摩擦角為35.5°，巖石黏聚力b為 9.05 MPa，摩擦角為38.8°，其他參數(shù)見表2。下面采用極限平衡法(平面滑移)和強(qiáng)度折減法分別計算這4種/工況時邊坡的穩(wěn)定性。

圖9 單臺階結(jié)構(gòu)面邊坡模型

4.1 極限平衡計算

臺階狀潛在滑動面和圖6所示的微裂紋分布圖表明臺階狀滑動帶主要沿臺階起伏區(qū)域損傷貫通。故假設(shè)如圖9所示臺階起伏的虛線區(qū)域為滑動帶，且臺階起伏的中線視為宏觀滑動面。按照圖9計算模型，則穩(wěn)定性系數(shù)，即

式中：R為抗滑力；為宏觀滑裂面的傾角；為滑體的質(zhì)量。

臺階狀結(jié)構(gòu)面抗剪強(qiáng)度采用經(jīng)驗式(8)確定，則抗滑力R可表示為

4.2 強(qiáng)度折減法計算

強(qiáng)度折減法的基本原理是將巖體抗剪強(qiáng)度參數(shù)黏聚力和摩擦角同時乘以一折減系數(shù)s，得到一組新的′和′，然后作為新的材料參數(shù)進(jìn)行計算，當(dāng)邊坡失穩(wěn)時，對應(yīng)的s稱為邊坡的安全系數(shù)。

采用FLAC?2D進(jìn)行強(qiáng)度折減法安全系數(shù)計算，計算模型與圖9所示的一致。臺階狀結(jié)構(gòu)面采用接觸面單元interface模擬。邊坡左邊界、右邊界及下邊界均采用法向位移約束。利用剪應(yīng)變增量、塑性區(qū)范圍及連通狀態(tài)和迭代求解的不收斂性綜合判定臨界穩(wěn)定狀態(tài)。同時對結(jié)構(gòu)面及圖9中虛線區(qū)域內(nèi)潛在滑動帶巖體的抗剪強(qiáng)度參數(shù)按相同的折減系數(shù)s折減，從而求得邊坡安全系數(shù)。

圖10和圖11所示分別為單臺階(高長比=0.1)結(jié)構(gòu)面邊坡臨界滑動時塑性區(qū)與剪切應(yīng)變增量分布圖。此時塑性區(qū)基本貫通，剪切滑動帶基本形成，且大致沿圖9中所假設(shè)的損傷帶滑動。臺階屈服且其頂部尖角剪應(yīng)變增量明顯集中。

圖10 邊坡塑性區(qū)

圖11 邊坡剪切應(yīng)變增量

4.3 2種方法比較分析

表3所示為2種方法穩(wěn)定性計算結(jié)果。雖然極限平衡計算的穩(wěn)定性系數(shù)與強(qiáng)度折減法計算的安全系數(shù)s存在一定的差別(這與其基本概念的定義不同相關(guān))，其中強(qiáng)度折減法計算結(jié)果相對偏大。但這2種方法計算結(jié)果隨臺階高長比的變化規(guī)律基本一致，也可表明本文建立的強(qiáng)度估算式(7)及式(8)基本合理，同時也反映了顆粒流程序計算剪切強(qiáng)度的可靠性。

表3 極限平衡法和強(qiáng)度折減法穩(wěn)定性計算結(jié)果

Table 3 Results of stability calculation using limit equilibrium methods and strength reduction methods

5 結(jié)論

1)巖體結(jié)構(gòu)面臺階可分為壓切損傷破壞(臺階高長比≤0.45)和張剪損傷破壞(＞0.45)模式。隨著增大，彎矩作用逐漸增強(qiáng)，由壓切逐漸轉(zhuǎn)化為張剪模式。但其損傷破壞模式與法向應(yīng)力基本無關(guān)。

2) 對于張剪損傷破壞模式，凸起臺階既受彎矩作用又受剪力作用，臺階根部易于在彎矩作用下張裂而形成張拉裂紋。臺階兩角點附近的宏觀張性裂紋傾伏向下且傾伏角隨減小而減小，當(dāng)=0.45時幾乎為0°。

3) 對于壓切損傷破壞模式，臺階主要受斜向上的剪切力作用而使其尖角剪斷，沿剪切面呈現(xiàn)明顯的爬坡效應(yīng)，具剪脹性質(zhì)。

4) 實際巖體結(jié)構(gòu)面臺階高長比常較小，針對壓切損傷破壞模式，建立了同時考慮結(jié)構(gòu)面爬坡效應(yīng)和啃斷效應(yīng)的剪切強(qiáng)度經(jīng)驗公式，通過邊坡穩(wěn)定性算例分析，驗證了此公式的合理性。

[1] 劉博, 李海波, 劉亞群. 循環(huán)剪切荷載作用下巖石節(jié)理變形特性試驗研究[J]. 巖土力學(xué), 2013, 34(9): 2475?2481. LIU Bo, LI Haibo, LIU Yaqun. Experimental study of deformation behavior of rock joints under cyclic shear loading[J]. Rock and Soil Mechanics, 2013, 34(9): 2475?2481.

[2] 王姣, 李長冬, 雍睿, 等. 三峽庫區(qū)侏羅系不同類型結(jié)構(gòu)面抗剪強(qiáng)度對比[J]. 煤田地質(zhì)與勘探, 2014, 42(5): 61?66. WANG Jiao, LI Changdong, YONG Rui, et al. Comparison study on shear strength parameters of Jurassic different rock discontinuities in Three Gorges Reservoir region[J]. Coal Geology & Exploration, 2014, 42(5): 61?66.

[3] INDRARATNA B, JAYANATHAN M, BROWN T. Shear strength model for overconsolidated clay-infilled idealised rock joints[J]. Géotechnique, 2008, 58(1): 55?65.

[4] 孫輔庭, 佘成學(xué), 萬利臺. 充填水泥漿巖石節(jié)理峰值剪切強(qiáng)度模型[J]. 巖石力學(xué)與工程學(xué)報, 2014, 33(12): 2481?2489. SUN Futing, SHE Chengxue, WAN Litai. A peak shear strength model for cement filled rock joints[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(12): 2481?2489.

[5] BARTON N. Shear strength criteria for rock,rock joints,rockfill and rock masses: problems and some solutions[J]. Journal of Rock Mechanics & Geotechnical Engineering, 2013, 5(4): 249?261.

[6] 周瑞光. 巖體結(jié)構(gòu)面形態(tài)的力學(xué)效應(yīng)[J].工程勘察, 1981, 9(1): 46?51. ZHOU Ruiguang. Mechanics effect of discontinuity surface morphology in rock mass[J]. Geotechnical Investigation & Surveying, 1981, 9(1): 46?51.

[7] EINSTEIN H H, VENEZIANO D, BAECHER G B, et al. The effect of discontinuity persistence of rock slope stability[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 1983, 20(5): 227?236.

[8] HOMAND F, BELEM T, SOULEY M. Friction and degradation of rock structural plane surfaces under shear loads[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2001, 25(10): 973?999.

[9] 周輝, 孟凡震, 張傳慶, 等. 結(jié)構(gòu)面剪切過程中聲發(fā)射特性的試驗研究[J]. 巖石力學(xué)與工程學(xué)報, 2015, 34(S1): 2827?2836. ZHOU Hui, MENG Fanzhen, ZHANG Chuanqing, et al. Experimental study of acoustic emission characteristic of discontinuity under shearing condition[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(S1): 2827?2836.

[10] 李海波, 劉博, 馮海鵬, 等. 模擬巖石結(jié)構(gòu)面試樣剪切變形特征和破壞機(jī)制研究[J]. 巖土力學(xué), 2008, 29(7): 1741?1746. LI Haibo, LIU Bo, FENG Haipeng, et al. Study of deformability behavior and failure mechanism by simulating rock joints sample under different loading conditions[J]. Rock and Soil Mechanics, 2008, 29(7): 1741?1746.

[11] 沈明榮, 張清照. 規(guī)則齒型結(jié)構(gòu)面剪切特性的模型試驗研究[J]. 巖石力學(xué)與工程學(xué)報, 2010, 29(4): 713?719. SHEN Mingrong, ZHANG Qingzhao. Experimental study of shear deformation characteristics of rock mass discontinuities with regular surface roughness[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(4): 713?719.

[12] 張清照, 沈明榮, 丁其文. 結(jié)構(gòu)面在剪切狀態(tài)下的力學(xué)特性研究[J]. 水文地質(zhì)工程地質(zhì), 2012, 39(2): 37?41. ZHANG Qingzhao, SHEN Mingrong, DING Qiwen. Study on the mechanical properties of rock mass discontinuity under shear condition[J]. Hydrogeology and Engineering Geology, 2012, 39(2): 37?41.

[13] 郭志. 實用巖體力學(xué)[M]. 北京: 地震出版社, 1996: 30?68. GUO Zhi. Practical rock mass mechanics[M]. Beijing: Seismological Press, 1996: 30?68.

[14] Itasca Consulting Group Inc. Manual of particle flow code in 2-dimension: theory and background[M]. 3.1st ed. Minneapolis: Itasca Consulting Group Inc, 2004: 1?4.

[15] POTYONDY D O, CUNDALL P A. A bonded-particle model for rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(8): 1329?1364.

[16] POTYONDY D O. The bonded-particle model as a tool for rock mechanics research and application: current trends and future directions[J]. Geosystem Engineering, 2014, 17(6): 342?369.

[17] BEWICK R P, KAISER P K, BAWDEN W F. Shear rupture under constant normal stiffness boundary conditions[J]. Tectonophysics, 2014, 634(11): 76?90.

(編輯 劉錦偉)

Shear properties of through-going discontinuity with a rock step and its application

HUANG Da1, 2, HUANG Runqiu3, LEI Peng1

(1. School of Civil Engineering, Chongqing University, Chongqing 400045, China;2. State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, China;3. State Key Laboratory of Geohazard Prevention and Geoenvironment Protection,Chengdu University of Technology, Chengdu 610059, China)

The influence of the ratio () of height () to length () of a rock step on shear deformation and strength for through-going discontinuity was investigated by 2D particle flow code (PFC2D) simulation method. The results show that the damage-failure modes of rock step can be concluded into two basic types: compression-shear failure and tension-shear failure. With the increase of the/of the rock step, the damage-failure mode gradually evolves from compression-shear to tension-shear failure. And the damage-failure modes do not pertain to the normal stress. For the tension-shear failure mode (＞0.45), the tensile crack is easily formed by rifting in the root of the rock step as the raised step is subjected to bending moment and shear force. For the compression-shear failure mode (≤0.45), the shear dilatation characteristics can be shown mainly due to the action of the up-sloping shear force which can cause that the sharp corners of the rock step sheared off and generate climbing effect along shear surface. For the/of a rock step in the actual slopes is often small relatively, and an empirical shear strength equation is established in which both the

rock mass mechanics; discontinuity; shear deformation; shear strength; particle flow code (PFC) effects of climbing along shear surface and the rock step gnawed are considered. The empirical equation can be applied to the estimation of shear strength parameter of the discontinuity with a rock step to calculate the slope stability.

10.11817/j.issn.1672-7207.2016.06.027

TD313

A

1672?7207(2016)06?2015?08

2015?06?25；

2015?08?09

國家自然科學(xué)基金資助項目(41472245，41130745)；重慶市國土房管科技計劃項目(CQGT-KJ-2014049)；中央高校基本科研業(yè)務(wù)費重大項目(106112016CDJZR208804)(Projects(41172243, 41130745) supported by the National Natural Science Foundation of China; Project(CQGT-KJ-2014049) supported by the Chongqing Administration of Land, Resources and Housing; Project (106112016CDJZR208804) supported by the Fundamental Research Funds for the Central Universities)

黃達(dá)，博士，教授，博士生導(dǎo)師，從事巖土工程和工程地質(zhì)研究；E-mail：dahuang@cqu.edu.cn