馬斌　李昕堯　梁超

摘要：運(yùn)用有限差分軟件FLAC，建立了呷爬滑坡體三維模型，分別施加1倍和3倍El-Centro地震波。計(jì)算結(jié)果表明，1倍地震波作用下會發(fā)生一定的永久位移，但不會滑動失穩(wěn)；3倍地震波作用下將發(fā)生滑動失穩(wěn)。進(jìn)一步發(fā)現(xiàn)邊坡在1倍和3倍地震波激勵下的變形滑移特征不同，第一種情況屬于累積效應(yīng)階段，變形主要集中于滑體O，第二種情況屬于觸發(fā)效應(yīng)階段，滑體N的滑動變形同樣劇烈。文中從滑坡體的組成物質(zhì)和材料參數(shù)、滑動區(qū)域的地形因素和滑動區(qū)域的動力響應(yīng)三方面進(jìn)行了分析。運(yùn)用Newmark有限滑動位移法，基于滑體N和滑體O的加速度響應(yīng)時(shí)程，計(jì)算了兩者的永久位移，結(jié)果與數(shù)值計(jì)算吻合良好。表明滑體O處的動力響應(yīng)會產(chǎn)生更大的永久位移，從而驗(yàn)證了數(shù)值計(jì)算的結(jié)果。同時(shí)，基于結(jié)構(gòu)動力響應(yīng)而計(jì)算得到的永久位移也可以作為衡量地震作用對該處結(jié)構(gòu)變形影響的指標(biāo)，從而估計(jì)地震對局部結(jié)構(gòu)的破壞作用。

關(guān)鍵詞：呷爬滑坡；變形特征；Newmark 方法；永久位移；加速度響應(yīng)

中圖分類號：TU443 文獻(xiàn)標(biāo)識碼：A 文章編號：1672-1683（2017）05-0110-06

對于人口分布稀疏的山區(qū)，地震觸發(fā)的滑坡、山體崩塌以及泥石流等災(zāi)害所造成的損失往往比地震造成的直接損失更大。對于20世紀(jì)末到21世紀(jì)初發(fā)生的臺灣集集地震、汶川地震、雅安蘆山地震，均伴隨多次大規(guī)?；潞推渌紊鸀?zāi)害，其人員傷亡和財(cái)產(chǎn)損失不亞于一場地震。因此，地表運(yùn)動激勵下的滑坡體穩(wěn)定性分析逐漸被巖土工程界和地震工程界所重視。近年來，專家學(xué)者對于地震作用下邊坡的失穩(wěn)機(jī)制和動力特性進(jìn)行了大量研究。何蘊(yùn)龍等通過有限元對滑坡體進(jìn)行了動力分析，可以近似計(jì)算滑坡所受到的地震作用；祁升文等對各種情況下滑坡的動力響應(yīng)進(jìn)行了數(shù)值模擬與分析，總結(jié)出了高、低邊坡的動力反應(yīng)形式；言志信等對于順層巖質(zhì)滑坡體地震激勵下的加速度響應(yīng)規(guī)律、頻譜特性、破裂面發(fā)生位置以及失穩(wěn)判定條件進(jìn)行了研究，并對耦合地震波作用、多級平臺和黃土地質(zhì)等情況進(jìn)行了較為全面的研究；徐光興等結(jié)合Newmark方法，提出了根據(jù)加速度預(yù)測邊坡永久位移的數(shù)學(xué)模型；張雪東等運(yùn)用FLAG3D有限差分商業(yè)軟件分析了錦屏庫區(qū)呷爬滑坡在靜力作用下的變形特征和穩(wěn)定性；殷躍平等結(jié)合大量的現(xiàn)場資料，分析了大光包巨型滑坡的失穩(wěn)機(jī)制和動力響應(yīng)。本文根據(jù)呷爬滑坡體變形監(jiān)測圖及地形圖，建立了滑坡體模型，運(yùn)用FLAC對其動力作用下的變形特征進(jìn)行了分析，并結(jié)合New-mark數(shù)值計(jì)算法利用邊坡地震響應(yīng)計(jì)算了邊坡的永久位移。

1模型建立

1.1滑坡體地質(zhì)概況

呷爬滑坡位于錦屏水電站庫區(qū)上游右岸，距離壩址約11.5 km。雅礱江流經(jīng)該滑體處的河道較順直，流向約為N12°E。該處河谷寬闊，左岸山體較為平緩，并有河漫灘、階地發(fā)育；右岸山體呈反向坡，岸坡的坡度約為35°。

雅礱江岸坡巖層產(chǎn)狀在呷爬滑坡處呈N 10～30°E/SE ∠85°，巖層主要由三疊系雜谷腦組深-灰黑色粉砂質(zhì)板巖夾變質(zhì)細(xì)砂巖組成。

呷爬滑坡體的平面投影呈長條狀，具有完整的圈谷地貌?；w靠近雅礱江一側(cè)高程較低，為l 655 m，與枯水期雅礱江水位相同，靠近山體一側(cè)高程較高，為2 120m，滑體前后緣高差為465 m?；麦w的縱長約為880m，寬度在260到300 m之間，表面面積約為Q 28 km²，滑坡殘?bào)w體積約為1 300萬m³。

滑坡表面可見3級平臺，3級平臺高程分別為1775～1840 m、1880～1925m、2010～2100 m，坡度為17～18°，原為耕地。在第一和第二級平臺中發(fā)育有兩條縱向沖溝。其中，一條位于滑坡軸線的北側(cè)（1號溝），另一條位于滑坡軸線的南側(cè)（2號溝），南側(cè)沖溝比北側(cè)沖溝更為深長。

在滑坡體臨江側(cè)到中部，表層巖體主要由塊碎石土組成，塊碎石土的成分為粉砂質(zhì)板巖、泥巖和礫巖；在滑坡中部到靠近山體一側(cè)，表層巖體主要由孤石，塊碎石土，解體不完全的粉砂質(zhì)板巖、變質(zhì)細(xì)砂巖和礫巖等組成；滑體物質(zhì)分層明顯，據(jù)鉆孔揭露，孔深27.97～39.57 m以上，為褐黃、灰黃色塊碎石土，以下至77.08～81.80 m為破碎的巖石碎塊夾黏土，局部巖體保留有原始層面。

在滑坡臨江側(cè)和兩條沖溝內(nèi)，可見明顯的滑帶出露。出露的滑帶在剪出口處反翹，根據(jù)地址勘察結(jié)果，其產(chǎn)狀為N15°E/SE∠22°?；瑤У闹饕M成物質(zhì)為灰黑色泥夾碎石并帶有大量石英顆粒，其中碎石成分為炭質(zhì)板巖、變質(zhì)細(xì)砂巖等，滑帶在滑體范圍內(nèi)平均厚度為1～3 m，局部可達(dá)8 m，上述滑帶組成物質(zhì)在干燥時(shí)剛度較大，含有一定水分后剛度顯著減小。

1.2建立數(shù)值模型

模型建立是FLAC^3D的短板，本文采用先在CAD中生成面域，導(dǎo)入AN SYS劃分網(wǎng)格，再利用清華大學(xué)鄭文棠博士開發(fā)的ANSYS2FLAC3D程序?qū)?fù)雜的地質(zhì)體模型數(shù)據(jù)變換為FLAC³⁰內(nèi)部代碼，從而實(shí)現(xiàn)模型導(dǎo)入。

為了模擬地震波在滑坡體中的傳播，避免在不同介質(zhì)的分界面上產(chǎn)生與實(shí)際情況不符的反射和折射，必須設(shè)置合理的邊界條件。本文在有限元模型四周網(wǎng)格設(shè)置自由場邊界條件，它的原理是在邊界節(jié)點(diǎn)上添加水平和垂直阻尼器來達(dá)到邊界條件與主體網(wǎng)格的耦合，并將自由場網(wǎng)格所受的不平衡力傳遞給主體網(wǎng)格。自由場邊界條件對入射角大于30°的波具有良好的衰減效果。對于呷爬滑坡體的基巖，由于其變形模量較大，假定為剛性地基，因此底部可不設(shè)置自由場邊界條件，直接輸入地震加速度時(shí)程。邊界條件見圖2。

本文采用局部阻尼形式，局部阻尼通過增加或減少節(jié)點(diǎn)的質(zhì)量，達(dá)到衰減振動的目的。阻尼比取為0.157。表1給出了滑坡體數(shù)值模型建立需要的材料。

根據(jù)文獻(xiàn)[15-16]的研究，網(wǎng)格的最大尺寸應(yīng)不大于最小波長的1/10到1/8，兼顧模型的計(jì)算速度，網(wǎng)格最大尺寸取為15 m，局部網(wǎng)格更小。endprint

2邊坡的變形特征分析

2.1地震波輸入

根據(jù)呷爬滑坡現(xiàn)場覆蓋層剪切波速大約在460～490 m/s，地表覆蓋層平均厚度約為90 m。根據(jù)文獻(xiàn)[17]規(guī)定，該場地為二類場地，因此選擇適用于二類場地的El-Centro南北向地震波經(jīng)濾波和基線調(diào)整后，作為輸入地震波，取其持續(xù)時(shí)間為28 s，峰值加速度為3.417 m/s²，加速度時(shí)程曲線見圖3。

邊坡在動力作用下產(chǎn)生的位移場直接反應(yīng)了其變形特征和破壞程度。加速度場作為動力響應(yīng)的最基本資料，具有重要意義。因此本文主要從位移場分析滑坡體的變形特征和穩(wěn)定狀態(tài)，并基于New-mark-β法和模型加速度響應(yīng)計(jì)算了模型不同點(diǎn)處的永久位移。

2.2不同狀態(tài)下的邊坡變形特征

圖4給出了滑坡體位移監(jiān)測點(diǎn)布置隋況。為了分析不同情況下的滑坡體變形特征，本文分別將El-Centro加速度時(shí)程乘以1和3的系數(shù)輸入，其中1倍和3倍加速度時(shí)程峰值為3.417 m/s²和10.251m/s²，分別代表一般地震作用和較極端的地震作用，從而綜合分析兩種情況下滑坡的變形特征和失穩(wěn)破壞情況，以揭示同一滑坡體在不同峰值特性動力輸入情況下的變形特征的不同。對于單個(gè)單元，如果應(yīng)力等于屈服應(yīng)力，將產(chǎn)生較大的塑性變形。但對于整個(gè)滑坡體模型，只有當(dāng)滑體或滑動面上所有單元均處于塑眭變形狀態(tài)且應(yīng)變大于一定值時(shí)，滑動才會發(fā)生?？梢灶A(yù)見，當(dāng)?shù)卣饎舆_(dá)到峰值時(shí)，監(jiān)測點(diǎn)對應(yīng)單元容易進(jìn)入塑性狀態(tài)，但當(dāng)?shù)卣鸷笃诘卣饎于呌谛》椒€(wěn)振動時(shí)，塑性單元將受到周圍處于彈性變形的單元或其它邊界約束條件的影響，使其變形受到限制。這時(shí)，若監(jiān)測點(diǎn)位移仍具有變形趨勢（即殘余速度），則表明監(jiān)測點(diǎn)附近單元或滑面對應(yīng)單元均進(jìn)入塑性變形階段，位移持續(xù)增大，滑坡體破壞。因此，本文將監(jiān)測點(diǎn)產(chǎn)生一定的永久變形，并且在地震動結(jié)束時(shí)仍然具有殘余速度作為滑坡體滑動失穩(wěn)的判據(jù)。圖5給出了滑坡體位移云圖和關(guān)鍵點(diǎn)監(jiān)測位移時(shí)程，施加1倍地震波后，所有測點(diǎn)先急劇變形，而后在15 s左右變形停止，基本處于穩(wěn)定狀態(tài)，直至地震結(jié)束，不符合滑動破壞判據(jù)，因此認(rèn)為滑坡體產(chǎn)生了一定的永久位移，但并未滑動失穩(wěn)。3倍地震波輸入后，滑體N和滑體O持續(xù)大幅變形，地震結(jié)束時(shí)變形呈發(fā)散狀態(tài)，監(jiān)測點(diǎn)仍具有殘余速度（E點(diǎn)速度為0.17 m/s），符合滑動判據(jù)，因此認(rèn)為滑體N和滑體O處于滑動失穩(wěn)狀態(tài)；小滑坡體測點(diǎn)變化趨勢與1倍地震輸入相同，因此處于穩(wěn)定狀態(tài)。

值得注意的是兩種情況下滑坡體的變形滑移特征不同，施加1倍El-Centro地震波時(shí)，位移集中分布于滑體O的凸起處，最大產(chǎn)生3.3 m左右的永久位移。施加3倍El-Centro地震波情況下，滑體N的滑動位移顯著增加，與滑體O相近。由圖5（b）可見，位移較大區(qū)域分為兩塊，分別位于滑體N和滑體O。根據(jù)實(shí)際算例表明，隨著地震動峰值的繼續(xù)增大，滑體N的運(yùn)動速度更快，位移更大，滑動更加劇烈。

張倬元等對于巖體邊坡的滑動機(jī)制理論可以較好地解釋本文的計(jì)算結(jié)果，該理論認(rèn)為地震對巖體的作用主要表現(xiàn)在累積效應(yīng)和觸發(fā)效應(yīng)兩個(gè)方面，前者主要引起邊坡巖土體結(jié)構(gòu)的松動，軟弱結(jié)構(gòu)面錯(cuò)位等失穩(wěn)前兆；后者主要表現(xiàn)在地震作用造成的邊坡中軟弱夾層流變和液化以及邊坡瞬間失穩(wěn)等。該理論與Newmark有限滑動位移理論一樣，都認(rèn)為地震可能造成滑坡體產(chǎn)生一定的永久位移，但邊坡是否滑動取決于位移的大小以及巖土體內(nèi)部的應(yīng)力應(yīng)變情況，滑動后的情況則與地形，滑體物質(zhì)組成，滑體的啟動條件等諸多復(fù)雜因素相關(guān)。

因此，可以認(rèn)為施加1倍地震波時(shí)，滑體O的累積效應(yīng)非常明顯，產(chǎn)生較大的永久位移。一旦滑動區(qū)由累積效應(yīng)過渡到觸發(fā)效應(yīng)后，滑體N的滑動位移顯著增加，與滑體O相近。主要原因歸納為三點(diǎn)：

（1）整個(gè)滑坡體前緣至中部，表層組成物質(zhì)主要為塊碎石土；中部及后緣，表層組成物質(zhì)主要為含孤石及塊石較多的塊碎石土、未完全解體的變質(zhì)粉砂質(zhì)板巖及變質(zhì)細(xì)砂巖、礫巖?；wO巖體組織較為松散，更易變形，從數(shù)值模型材料參數(shù)上來看其模量較??；

（2）滑體O位于坡面局部凸起處，其下滑坡度較緩，滑體N所處位置坡面較陡，一旦進(jìn)入觸發(fā)效應(yīng)階段，相同條件下其速度和位移勢必更加劇烈；

（3）施加1倍El-Centro波時(shí)，滑體本身的動力反應(yīng)導(dǎo)致滑體O的永久位移要大于滑體N。分析過程如第3節(jié)所示。

從數(shù)值分析結(jié)果看，小滑坡的穩(wěn)定狀態(tài)好于主滑坡。

3利用Newmark方法估計(jì)邊坡動力響應(yīng)對變形的影響

3.1屈服加速度的確定

根據(jù)Newmark方法，假設(shè)滑動面為圓弧形，見圖6，由動力穩(wěn)定系數(shù)定義，有下面兩式：

本文運(yùn)用自編強(qiáng)度折減法程序，對組成模型的所有材料的內(nèi)聚力和內(nèi)摩擦角進(jìn)行折減，得到滑體N為最危險(xiǎn)滑動體，其安全系數(shù)為1.4。然后保持滑體N的材料參數(shù)不變，對其他材料進(jìn)行參數(shù)折減，得到最危險(xiǎn)滑動體為滑體O，其安全系數(shù)為1.5。見圖7，與動力計(jì)算得出的兩個(gè)滑動危險(xiǎn)區(qū)域相對應(yīng)。

本文將兩個(gè)滑體的滑動面概化為圓弧，見圖8，由式（7）可分別求出其屈服加速度。

3.2滑體N和滑體O的變形情況分析

選擇滑動處某點(diǎn)的加速度響應(yīng)時(shí)程，根據(jù)Ne-wmark有限滑動位移原理，方向指向坡外，絕對值大于屈服加速度的加速度時(shí)程部分，才能產(chǎn)生永久位移，稱為有效加速度。

圖9中分別為滑體N和滑體O上某點(diǎn)的有效加速度時(shí)程以及由有效加速度時(shí)程二次積分求得的永久位移。計(jì)算得滑體N永久位移為2.185 m，滑體O永久位移為3.102 m。與數(shù)值模型結(jié)果符合較好。

4結(jié)語

（1）根據(jù)呷爬滑坡體的數(shù)值分析結(jié)果，在地震達(dá)到一定程度時(shí)，滑體N和滑體O將首先發(fā)生滑動，并且在1倍和3倍El-Centro波作用下具有不同的變形滑動機(jī)制：在1倍地震波作用下，滑坡體處于累積效應(yīng)階段，滑體O的永久位移較大；在3倍地震波作用下，轉(zhuǎn)入觸發(fā)效應(yīng)階段，滑體N的滑動速度及位移顯著增加，與滑體O相近。本文分別從滑體物質(zhì)組成和材料參數(shù)、滑動地形和動力響應(yīng)三個(gè)方面進(jìn)行了分析。從數(shù)值計(jì)算結(jié)果看，小滑坡的穩(wěn)定性要好于主滑坡。

（2）本文根據(jù)Newmark有限滑動位移理論，基于滑體N和滑體O處某點(diǎn)的加速度響應(yīng)時(shí)程，估計(jì)了滑坡體的永久位移，與數(shù)值計(jì)算結(jié)果對比符合較好。它還可以作為邊坡動力響應(yīng)對滑坡體影響的一個(gè)指標(biāo)，可以估計(jì)地震作用對局部區(qū)域的變形影響。但必須指出，由于動力響應(yīng)的非線性和共振現(xiàn)象，同一區(qū)域內(nèi)不同點(diǎn)的加速度時(shí)程所產(chǎn)生的永久位移必然存在差異，并不是所有加速度時(shí)程都可以代表整個(gè)滑坡體進(jìn)行永久位移的計(jì)算，因而，在加速度時(shí)程的選取方面還有待進(jìn)一步改進(jìn)。endprint