馬斌　李昕堯　梁超

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

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

中圖分類號：TU443 文獻標識碼：A 文章編號：1672-1683（2017）05-0110-06

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

1模型建立

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

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

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

呷爬滑坡體的平面投影呈長條狀，具有完整的圈谷地貌。滑體靠近雅礱江一側(cè)高程較低，為l 655 m，與枯水期雅礱江水位相同，靠近山體一側(cè)高程較高，為2 120m，滑體前后緣高差為465 m。滑坡體的縱長約為880m，寬度在260到300 m之間，表面面積約為Q 28 km²，滑坡殘體體積約為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ì)細砂巖和礫巖等組成；滑體物質(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ì)細砂巖等，滑帶在滑體范圍內(nèi)平均厚度為1～3 m，局部可達8 m，上述滑帶組成物質(zhì)在干燥時剛度較大，含有一定水分后剛度顯著減小。

1.2建立數(shù)值模型

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

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

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

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

2邊坡的變形特征分析

2.1地震波輸入

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

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

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

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

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

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

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

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

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

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

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

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

3.1屈服加速度的確定

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

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

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

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

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

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

4結(jié)語

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

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