江婷+沈振中+聶柏松+劉沖+談家誠(chéng)

摘要：基于離散元法，結(jié)合某水電站廠房后緣順層巖質(zhì)邊坡工程，研究緩傾角多結(jié)構(gòu)面巖質(zhì)邊坡的地震響應(yīng)規(guī)律及錨桿的加固效果，提出了該邊坡抗震加固優(yōu)化方案。研究表明：未加固順層塊狀巖質(zhì)邊坡在地震過(guò)程中位移呈現(xiàn)階梯式增長(zhǎng)，地震結(jié)束時(shí)刻邊坡巖塊松散崩落；加固后邊坡位移顯著減小，且能在地震作用下保持穩(wěn)定；錨桿的剪力響應(yīng)與其附近結(jié)構(gòu)面的相對(duì)位移相關(guān)性強(qiáng)，表明錨桿附近巖層面間的相對(duì)位移是激發(fā)錨桿加固作用的主要因素，在坡面出露且出露點(diǎn)低的巖層面在地震作用下相對(duì)位移較大，對(duì)應(yīng)的錨桿承受剪力也較大，據(jù)此可更合理地布置錨桿。

關(guān)鍵詞：巖質(zhì)邊坡；結(jié)構(gòu)面；離散元法；地震響應(yīng)；錨桿；加固

中圖分類號(hào)：TU 457 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-1683（2017）04-0136-06

Abstract：The seismic responses and anchor reinforcement effects of a bedding stereoplasm rock slope with multiple structural planes behind the powerhouse of a hydro-power plant by use of the DEM，and a optimal layout plan of the anchors is proposed.The research shows that displacements of the unreinforced bedding stereoplasm slope under seismic action is present the stepped increases，and the rock blocks fall at the end of earthquake.Displacements of reinforced slope diminish significantly，and the slope can remain stable in the earthquake.The shear force response of the anchor elements shows good agreement with the relative displacement response of the lay of terrain around the anchor element，which means that the relative displacement of the lay of terrain around the anchor is the main stimulation factor of the anchor reinforcement effects in the earthquake and the relative displacement of the lay of terrain emerging on slope surface and at a lower elevation is larger so that the shear force of the corresponding anchor is larger.Hereby this paper put forward the reasonable arrangement schemes of the anchor.

Key words：rock slope；structural plane；DEM；seismic response；anchor；strengthening

節(jié)理化緩傾角巖質(zhì)邊坡在水電工程中大量存在，地震誘發(fā)其失穩(wěn)是影響水電工程建設(shè)順利進(jìn)行的關(guān)鍵問(wèn)題[1-4]。加強(qiáng)對(duì)地震作用下邊坡的動(dòng)力響應(yīng)研究，為邊坡穩(wěn)定進(jìn)行預(yù)報(bào)預(yù)警提供力學(xué)理論依據(jù)，極具工程意義[5-8]。

巖質(zhì)邊坡位移的動(dòng)力響應(yīng)較為復(fù)雜，目前主要通過(guò)現(xiàn)場(chǎng)監(jiān)測(cè)法、模型試驗(yàn)法和數(shù)值模擬法得到邊坡位移動(dòng)力響應(yīng)?，F(xiàn)場(chǎng)監(jiān)測(cè)方法是最直接的方法，可信度高，但由于地震的不可預(yù)測(cè)性和監(jiān)測(cè)周期長(zhǎng)的原因，現(xiàn)場(chǎng)監(jiān)測(cè)方法人力和物力成本較高。室內(nèi)模型試驗(yàn)不必進(jìn)行復(fù)雜全尺寸試驗(yàn)來(lái)獲得真實(shí)試驗(yàn)成果，比較經(jīng)濟(jì)，但目前大多進(jìn)行的是簡(jiǎn)單概化的模型試驗(yàn)研究，地質(zhì)條件與加載條件較單一，無(wú)法精確模擬實(shí)際工程。數(shù)值模擬方法可以模擬分析復(fù)雜巖土體內(nèi)部的響應(yīng)現(xiàn)象和破壞機(jī)理，并把結(jié)果展示出來(lái)，因而許多工程中應(yīng)用較多的是基于數(shù)值模擬方法來(lái)進(jìn)行邊坡位移動(dòng)力響應(yīng)的。數(shù)值模擬法主要有有限元法和離散元法等。例如張江偉[9]利用有限元法，基于強(qiáng)度折減法探討了邊坡在地震作用下的穩(wěn)定性判定方法；彭寧波[10]等采用FLAC3D分析了波型以及振幅、及頻率等地震參數(shù)對(duì)順層巖質(zhì)邊坡的動(dòng)力響應(yīng)的影響；廖少波[11]等采用3DEC軟件，基于離散元方法分析地震波在節(jié)理巖體中的傳播特性，揭示地震作用對(duì)塊狀巖質(zhì)邊坡穩(wěn)定性的影響機(jī)制；鄭允[12]等采用UDEC軟件提出一種計(jì)算爆破開(kāi)挖時(shí)巖質(zhì)邊坡安全系數(shù)的節(jié)點(diǎn)擬靜力法，提出利用最大節(jié)點(diǎn)位移時(shí)程曲線來(lái)判斷邊坡臨界狀態(tài)。

目前對(duì)巖質(zhì)邊坡的失穩(wěn)判據(jù)已有了較深入的研究，一旦發(fā)現(xiàn)邊坡存在失穩(wěn)可能便要采取加固措施，而利用錨桿進(jìn)行邊坡治理是最常用的方法[13-15]，但如何合理有效地布置錨桿是加固設(shè)計(jì)的難題。本文基于離散元法，結(jié)合某水電站工程，應(yīng)用通用離散元程序?qū)Φ卣鹱饔孟聨r質(zhì)邊坡加固前后的動(dòng)力響應(yīng)進(jìn)行研究，分析地震作用下加固與未加固邊坡的位移狀況以及錨桿的受剪情況，提出錨桿優(yōu)化布置的建議。

[BT2+*8]1 離散元法動(dòng)力計(jì)算模型

離散元法在連續(xù)介質(zhì)力學(xué)的基礎(chǔ)上，加入接觸的概念，并將其作為連續(xù)體之間的邊界，并采用連續(xù)力學(xué)方法描述連續(xù)體和接觸的力學(xué)行為[16-17]。在UDEC模型中，把巖體定義為相互接觸的塊體的組合，塊體代表巖體物理力學(xué)性質(zhì)，接觸代表結(jié)構(gòu)面兩側(cè)塊體間的相互作用。結(jié)構(gòu)面的切向剛度較低，可不斷地調(diào)整結(jié)構(gòu)面的物理力學(xué)參數(shù)（如切向剛度、法向剛度、內(nèi)摩擦角、黏聚力等）從而不斷地改變結(jié)構(gòu)面宏觀物理力學(xué)性能，從而使模型不斷地接近實(shí)際構(gòu)造。

UDEC軟件將巖體模型剖分成塊體，將塊體定義為可變形的，再講塊體劃分為由有限差分的三角形單元組成的網(wǎng)格，在計(jì)算過(guò)程中單元各節(jié)點(diǎn)都是可以運(yùn)動(dòng)的。當(dāng)研究區(qū)域的邊界條件及作用力系改變時(shí)，一些塊體將產(chǎn)生加速度并發(fā)生位移，塊體的位置將改變，移動(dòng)后的塊體將與其周圍的塊體發(fā)生“疊合”，出現(xiàn)新的作用力系，隨著力的傳遞，更多的塊體發(fā)生位移。

1.1 差分網(wǎng)格構(gòu)建

1.2 錨桿結(jié)構(gòu)模擬

UDEC軟件中采用cable單元[19]來(lái)模擬錨桿，即通過(guò)對(duì)錨桿、漿體和巖體單元邊界的節(jié)點(diǎn)的相對(duì)位移，表示錨桿、灌漿體以及巖土體之間的相對(duì)滑動(dòng)。在建模計(jì)算過(guò)程中，錨桿被分為許多雙彈簧單元，彈簧1表述錨桿單元之間的拉撥力，彈簧2表述錨桿單元與砂漿間的黏結(jié)力，見(jiàn)圖1。通過(guò)這些單元物理力學(xué)性質(zhì)和位移可求出錨桿的變形情況和應(yīng)力狀態(tài)。

1.3 邊界條件確定與地震波輸入

固定邊界會(huì)反彈地震波，而黏滯邊界[20]是在模型邊界施加法向阻尼器和切向阻尼器，阻尼器產(chǎn)生切向力和法向黏性力，將反射波引起的應(yīng)力抵消，能夠達(dá)到吸收入射波能量的目的。因此地震響應(yīng)計(jì)算時(shí)須采用黏滯邊界。

在地震響應(yīng)計(jì)算前須對(duì)地震波進(jìn)行基線校正和濾波處理。基線校正是在原始地震波時(shí)程上加入一個(gè)低頻的多項(xiàng)式函數(shù)波型，使動(dòng)力計(jì)算結(jié)束后模型速度和加速度為零。濾波處理是濾掉地震波中的高頻部分，因?yàn)橥ㄟ^(guò)功率譜分析發(fā)現(xiàn)輸入地震波的低頻部分的能量占了地震波總能量的絕大部分，而且離散化的介質(zhì)存在能量傳播的上限頻率，只有當(dāng)輸入的荷載頻率小于這個(gè)上限頻率時(shí)計(jì)算結(jié)果才有意義。

2 工程實(shí)例

2.1 地質(zhì)工程概況

某水電站主廠房后緣邊坡，開(kāi)挖后形成高約30 m的順層塊狀巖質(zhì)邊坡。根據(jù)工程地質(zhì)勘察報(bào)告可知，一級(jí)、二級(jí)和三級(jí)邊坡的傾角分別是72°、53°和72°，由傾角為56°的連續(xù)等距巖層面組和傾角為146°的不連續(xù)等距節(jié)理組切割形成塊狀邊坡。該邊坡具備有該類邊坡的一般特點(diǎn)，即巖性單一，邊坡結(jié)構(gòu)面主要是層面和節(jié)理，節(jié)理的間距與層厚成正比，層面結(jié)合程度較好，在地震作用下可能發(fā)生局部崩塌性破壞，危及廠房安全。根據(jù)邊坡工程地質(zhì)情況和可能的破壞型式以及錨桿布置情況，在對(duì)其合理簡(jiǎn)化的基礎(chǔ)上，建立如圖2所示的剖面。巖體和結(jié)構(gòu)面物理力學(xué)參數(shù)見(jiàn)表1和表2。工程中采用全長(zhǎng)粘結(jié)式砂漿錨桿，其布置間距為2.5 m，錨桿布設(shè)角度和長(zhǎng)度見(jiàn)表3，錨桿物理力學(xué)參數(shù)見(jiàn)表4。巖體材料采用理想彈塑性模型，破壞準(zhǔn)則為M-C準(zhǔn)則，采用面接觸的Coulomb滑動(dòng)模型模擬邊坡的結(jié)構(gòu)面和節(jié)理面。

在地震過(guò)程中，未加固順層塊狀巖質(zhì)邊坡的監(jiān)測(cè)點(diǎn)位移在3～5 s、9～10 s、14～16三個(gè)時(shí)間段突增，其他時(shí)間段保持平穩(wěn)，呈現(xiàn)階梯式增長(zhǎng)，這說(shuō)明位移具有累積效應(yīng)，邊坡巖塊崩落發(fā)生在巖塊位移累積到一定值時(shí)，而不是地震加速度最大時(shí)。在地震結(jié)束時(shí)刻，一級(jí)邊坡坡面巖塊松散滑移崩落，部分已滑落到坡腳，邊坡處于不穩(wěn)定狀態(tài)。一級(jí)邊坡上的巖塊易滑移崩落，是因?yàn)樵诘卣疬^(guò)程中，巖塊所處位置越低處，所處位置為潛在滑動(dòng)面時(shí)滑動(dòng)體積越大，越易滑動(dòng)。由分析可知，該邊坡需進(jìn)行加固。

2.3 加固后邊坡的動(dòng)力響應(yīng)特性

2.3.1 邊坡位移響應(yīng)分析

UDEC能夠記錄邊坡在地震過(guò)程中任意時(shí)刻下位移最大的節(jié)點(diǎn)的位移值，加固前后的邊坡最大位移時(shí)程曲線見(jiàn)圖6，加固后地震結(jié)束時(shí)刻邊坡單元節(jié)點(diǎn)位移情況見(jiàn)圖7。

全長(zhǎng)粘結(jié)式錨桿與周圍巖體有很大的粘結(jié)剛度，能承受很大的軸向拉撥力，使錨桿在地震過(guò)程中不發(fā)生滑脫。在地震結(jié)束時(shí)刻，加固前邊坡發(fā)生破壞，其最大位移值達(dá)到11.97 m，此時(shí)坡體的最大位[CM（22]移代表坡體巖塊的最大滑落距離；加固后邊坡坡面完好、未發(fā)生巖塊滑落，最大位移為0.20 m，比未加固邊坡減少了近98%，且最大位移曲線趨于平穩(wěn)沒(méi)有上升趨勢(shì)。這表明錨桿加固效果顯著，加固后邊坡在地震作用下未發(fā)生破壞。

2.3.2 錨桿剪力分析

對(duì)于順層邊坡，全長(zhǎng)粘結(jié)式錨桿的主要作用是阻止邊坡巖體沿著層面滑動(dòng)破壞，即抗剪作用大于抗拉作用，故本文重點(diǎn)分析錨桿的剪力。典型時(shí)刻各錨桿最大剪力時(shí)程曲線見(jiàn)圖8，地震結(jié)束時(shí)刻錨桿剪力分布情況見(jiàn)圖9，錨桿截面上的剪力對(duì)錨桿上任意一點(diǎn)矩為順時(shí)針轉(zhuǎn)向時(shí)為正，反之為負(fù)。

地震作用下，巖層面錯(cuò)動(dòng)較大因此巖層面附近錨桿剪力較大。當(dāng)錨桿穿過(guò)多個(gè)巖層面時(shí)，錨桿與在坡面出露的巖層面相交的部位剪力較大，且出露[CM（22]點(diǎn)越低剪力越大，因?yàn)槌雎饵c(diǎn)越低的巖層面為潛在滑動(dòng)面時(shí)滑動(dòng)體積越大，結(jié)構(gòu)面相對(duì)位移越大。結(jié)合邊坡的具體情況分析可知，第一級(jí)邊坡的錨桿剪力最大，第三級(jí)邊坡上錨桿剪力次之，第二級(jí)邊坡上的錨桿剪力最小，這和計(jì)算結(jié)果相吻合。

2.4 加固錨桿布置優(yōu)化

加固后邊坡在給定的地震荷載下能保持穩(wěn)定，但從巖層面對(duì)巖質(zhì)邊坡動(dòng)力響應(yīng)影響及加固后邊坡的錨桿剪力分析來(lái)看，工程中所采用的加固措施可以適當(dāng)優(yōu)化，具體如下：

（1）第二級(jí)邊坡坡度較緩，且坡面并無(wú)巖層面出露，錨桿的加固作用并沒(méi)有充分發(fā)揮，地震結(jié)束時(shí)刻各錨桿最大剪力值也驗(yàn)證這個(gè)推論，第二級(jí)邊坡的錨桿最大剪力僅為34 kN。因此，可考慮不在該級(jí)邊坡布置錨桿；

（2）第三級(jí)邊坡最下部3支錨桿所穿過(guò)的巖層面在整個(gè)邊坡坡面均未出露，這三支錨桿的抗震加固作用也沒(méi)有充分發(fā)揮，且地震結(jié)束時(shí)刻這3支錨桿的最大剪力僅為15.7 kN。因此，可考慮不布置編號(hào)為5、6、7的3支錨桿；

（3）對(duì)于第三級(jí)邊坡最上部編號(hào)為1、2、3和4的4支錨桿，從地震結(jié)束時(shí)刻錨桿剪力值上看，最大剪力只有56，相對(duì)較小，但這4支錨桿均穿過(guò)了在第三級(jí)邊坡坡面出露的巖層面，取消這4支錨桿可能會(huì)導(dǎo)致坡體沿巖層面滑動(dòng)，所以不能取消布置這4支錨桿。

綜上，取消布置編號(hào)為5～14的錨桿。

2.5 優(yōu)化加固后邊坡的動(dòng)力響應(yīng)特性

優(yōu)化加固后地震結(jié)束時(shí)刻邊坡單元節(jié)點(diǎn)位移情況見(jiàn)圖6，優(yōu)化加固前后最大位移時(shí)程曲線見(jiàn)圖7。

在地震過(guò)程中，優(yōu)化加固邊坡未發(fā)生破壞，其位移和原加固邊坡位移在方向和大小上基本保持一致。地震結(jié)束時(shí)刻優(yōu)化加固后的邊坡保持穩(wěn)定，其最大位移時(shí)程曲線趨于穩(wěn)定值，大小為0.21 m，僅比原加固邊坡多了0.01 m。可判斷優(yōu)化加固方案與原加固方案效果相差不大，均能使邊坡在給定地震荷載下保持穩(wěn)定。

3 結(jié)論

對(duì)邊坡進(jìn)行地震響應(yīng)過(guò)程模擬，有助于人們及時(shí)采取適當(dāng)?shù)募庸檀胧苊鉃?zāi)難發(fā)生。本文采用UDEC，針對(duì)某水電站廠房后緣順層塊狀巖質(zhì)邊坡，計(jì)算分析了緩傾角多結(jié)構(gòu)面復(fù)雜巖質(zhì)邊坡的動(dòng)力響應(yīng)規(guī)律及錨桿的加固效果。主要結(jié)論如下：

（1）離散元法能夠很好地模擬含有巖層面、節(jié)理面等非連續(xù)結(jié)構(gòu)面的巖質(zhì)邊坡的動(dòng)力響應(yīng)，即可以得到任何時(shí)刻下坡體與錨桿單元各節(jié)點(diǎn)的變形狀態(tài)與應(yīng)力狀態(tài)，模擬出坡體漸進(jìn)破壞的過(guò)程，其非連續(xù)性的破壞特征是有限元法等連續(xù)性方法所難以模擬的。

（2）基于該方法研究發(fā)現(xiàn)：未加固順層塊狀巖質(zhì)邊坡在地震作用下，其破壞模式以一級(jí)邊坡上局部巖塊掉落為主，因?yàn)橛袔r層面出露且出露點(diǎn)低的坡面部位更可能發(fā)生巖塊滑落。而加固后邊坡能夠在地震作用下保持穩(wěn)定。穿過(guò)在坡面出露的巖層面的錨桿才能真正發(fā)揮作用，據(jù)此可優(yōu)化錨桿的布置方案。

在其他工程中，巖質(zhì)邊坡的結(jié)構(gòu)面分布形式多樣，其破壞形式也很多。本文只是討論了順層塊狀巖質(zhì)邊坡加固前后的地震響應(yīng)情況，而在其他結(jié)構(gòu)面分布形式下的邊坡穩(wěn)定性則有待進(jìn)一步研究。

參考文獻(xiàn)（References）：

[1] 言志信，張森，張學(xué)東.地震邊坡失穩(wěn)機(jī)理及穩(wěn)定性分析探討[J].工程地質(zhì)學(xué)報(bào)，2010，18（6）：844-850.（YAN Zhi-xin，ZHANG Sen，ZHANG Xue-dong .Failure mechanism and stability analysis of slope under earthquake[J].Journal of Engineering Geology，2010，18（6）：844-850.（in Chinese））

[2] 石崇，王盛年，劉琳.地震作用下陡巖崩塌顆粒離散元數(shù)值模擬研究[J].巖石力學(xué)與工程學(xué)報(bào)，2013，32（增1）：2798-2805.（SHI Chong，WANG Sheng-nian，LIU Lin.Research of avalanche disaster numerical simulation based on granular discrete element method of high-steep slope under seismic loads[J].Chinese Journal of Rock Mechanics and Engineering，2013，32（s1）：2798-2805.（in Chinese））

[3] 賀續(xù)文，劉忠，廖彪，等.基于離散元法的節(jié)理巖體邊坡穩(wěn)定性分析[J].巖土力學(xué)，2011，32（7）：2199-2204.（HE Xu-wen，LIU Zhong，LIAO Biao.Stability analysis of jointed rock slopes based on discrete element method[J].Rock and Soil Mechanics，2011，32（7）：2199-2204.（in Chinese））

[4] 寧宇，徐衛(wèi)亞，鄭文棠.應(yīng)用離散元強(qiáng)度折減對(duì)復(fù)雜邊坡進(jìn)行穩(wěn)定性分析[J].巖土力學(xué)，2007，28（增）：569-574.（NING Yu，XU Wei-ya，ZHENG Wen-tang.Application in complicated high slope with strength reduction method based on discrete element method[J].Rock and Soil Mechanics，2007，28（s）：569-574.（in Chinese））

[5] 劉天蘋，李世海，劉曉宇.節(jié)理化巖質(zhì)邊坡隨機(jī)結(jié)構(gòu)面有限元和離散元耦合計(jì)算方法研究[J].巖石力學(xué)與工程學(xué)報(bào)，2014，33（增1）：3114-3122.（LIU Tian-ping，LI Shi-hai，LIU Xiao-yu.Study of dem and fem combination method for stochastic structural plane in jointed rock slope[J].Chinese Journal of Rock Mechanics and Engineering，2014，33（s1）：3114-3122.（in Chinese））

[6] He J M，Li S D，Yin Y P.Study of seismic response of colluvium accumulation slope by particle flow code[J].Granular Matter，2010，12：483-490 （10.1007/s10035-010-0213-8）.

[7] 楊國(guó)香，伍法權(quán)，董金玉.地震作用下巖質(zhì)邊坡動(dòng)力響應(yīng)特性及變[JP+2]形破壞機(jī)制研究[J].巖石力學(xué)與工程學(xué)報(bào)，2012，31（4）：696-702.（YANG Guo-xinag，WU Fa-quan，DONG Jin-yu.Study of dynamic response characters and failure mechanism of rock slope under earthquake[J].Chinese Journal of Rock Mechanics and Engineering，2012，31（4）：696-702.（in Chinese））

[8] 鄭穎人，葉海林，黃潤(rùn)秋，等.地震邊坡穩(wěn)定性分析探討[J].地震工程與工程振動(dòng)，2010，30（2）：66-73.（ZHENG Ying-ren，YE Hai-lin，HUANG Run-qiu，et al.Study on the stability analysis of earthquake slope[J].Journal of Earthquake Engineering and Engineering Vibration，2010，30（2）：66-73.（in Chinese））

[9] 張江偉.土質(zhì)邊坡地震動(dòng)力響應(yīng)特牲分析及穩(wěn)定性評(píng)價(jià)[D].北京：中國(guó)地震局地球物理研究所，2015（ZHANG Jiang-wei，Seismic dynamic response analysis and stability evaluation of soil slope[D].Beijing：Insititute of Geophysics，China Earthquake Administration，2015.（in Chinese））

[10] 彭寧波.錨固巖質(zhì)邊坡地震動(dòng)力響應(yīng)及錨固機(jī)理研究[D].蘭州：蘭州大學(xué)，2014.（Peng Ningbo，Study on seismic responses and anchoring mechanism of anchored rock slope[D]..Lanzhou：Lanzhou University，2014.（in Chinese）

[11] 廖少波.強(qiáng)震作用下塊狀巖體邊坡穩(wěn)定性研究[D].武漢：中國(guó)地質(zhì)大學(xué)，2013.（LIAO Shao-bo.Study on stability of block rock slope under strong earthquake[D].Wuhan：China University of Geosciences，2013.（in Chinese）

[12] 鄭 允 ，陳從新，朱璽璽，等.基于UDEC 的巖質(zhì)邊坡開(kāi)挖爆破節(jié)點(diǎn)擬靜力穩(wěn)定性計(jì)算方法[J].巖石力學(xué)與工程學(xué)報(bào)，2014，33（增3）：3932-3940.（ZHENG Yun，CHEN Cong-xin，ZHU Xi-xi，et al.Node quasi-static stability analysis of rock slope underExcavation blasting based on udec[J].Chinese Journal of Rock Mechanics and Engineering，2014，33（s3）：3932-3940.（in Chinese））

[13] 葉海林，鄭穎人，黃潤(rùn)秋，等.巖質(zhì)邊坡描桿支護(hù)參數(shù)地震敏感性分析[J].巖土工程學(xué)報(bào)，2011，28（9）：1374-1379.（YE Hai-lin，ZHENG Ying-ren，HUANG Run-qiu，et al.Sensitivity analysis of parameters for bolts in rock slopes under earthquakes[J].Chinese Journal of Geotechnical Engineering，2011，28（9）：1374-1379.（in Chinese））

[14] 吳月秀，劉泉聲，盧海峰，等.UDEC 中Saeb-Amadei 節(jié)理本構(gòu)模型的二次開(kāi)發(fā)[J].巖石力學(xué)與工程學(xué)報(bào)，2012，31（增1）：3022-3030.（WU Yue-xiu，LIU Quan-sheng，LU Hai-feng，et al.Implementation of saeb-amadei model into user-definedjoint constitutive model in udec，[J].Chinese Journal of Rock Mechanics and Engineering，2012，31（s1）：3022-3030.（in Chinese））

[15] 高保彬，高佳佳，袁東升.基于UDEC 的大采高覆巖破裂的模擬與分析[J].湖南科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，28（2）：1-6.（GAO Bao-bin，GAO Jia-jia，YUAN Dong-sheng.Simulation and analysis of overlying strata fracture in large mining height top-coal caving face of coal seams based on UDEC[J].Journal of Hunan University of Science & Technology：Natural Science Edition，2013，28（2）：1-6.（in Chinese））

[16] Eberhardt E，Stead D，Coggan J S.Numerical analysis of initiation and progressive failure in natural rock slopes-the 1991 Randa rockslide[J].International Journal of Rock Mechanics and Mining Sciences，2004，41：69-87（10.1016/S1365-1609（03）00076-5）.

[17] Bonilla-Sierra V，Scholtes L，Donze FV.DEM analysis of rock bridges and the contribution to rock slope stability in the case of translational sliding failures[J].International Journal Of Rock Mechanics And Mining Sciences，2015，80：67-78.

[18] Kuhlmeyer R L，J.Lysmer.Finite Element Method Accuracy for Wave Propagation Problems[J].soil Mech & Foundations Div，ASCE，1973，99：421-427（10.1007/s11803-013-0196-7）.

[19] Itasca Consulting Group Inc.Theory and background[Z].Minnesota：Itasca Consulting Group，2002.

[20] Cundall P.A.UDEC-A generalized District Element Program for Modeling Jointed Roc[Z].European Research Office，U.S.Army.1980.

[21] SL203-97.水工建筑物抗震設(shè)計(jì)規(guī)范[S].（SL203-97.Code for seismic design of hydraulic structures.（in Chinese））