黏性土坡地震失穩(wěn)機制及其在黃土斜坡地震穩(wěn)定性分析中的應(yīng)用①

趙文琛， 吳志堅,2， 陳豫津

(1.中國地震局蘭州地震研究所 黃土地震工程重點實驗室，甘肅 蘭州 730000；2.甘肅省巖土防災(zāi)工程技術(shù)研究中心，甘肅 蘭州 730000)

黏性土坡地震失穩(wěn)機制及其在黃土斜坡地震穩(wěn)定性分析中的應(yīng)用①

趙文琛1， 吳志堅1,2， 陳豫津1

(1.中國地震局蘭州地震研究所 黃土地震工程重點實驗室，甘肅 蘭州 730000；2.甘肅省巖土防災(zāi)工程技術(shù)研究中心，甘肅 蘭州 730000)

摘要：基于地震作用下黏性土坡失穩(wěn)滑動特點，以土體應(yīng)力狀態(tài)及其變化分析邊坡失穩(wěn)過程。通過分析地震作用下邊坡不同部位土體應(yīng)力狀態(tài)和剪應(yīng)力變化，結(jié)合實際地震邊坡失穩(wěn)破壞特征，提出黏性土坡地震三段式滑動失穩(wěn)機制。在分析該滑動失穩(wěn)機制與有限元強度折減法之間應(yīng)力關(guān)聯(lián)的基礎(chǔ)上，將兩者結(jié)合應(yīng)用于實際黃土地震滑坡動力穩(wěn)定性分析。依據(jù)此考慮得到的動力安全系數(shù)相比較其他方法，與極限平衡法得到的結(jié)果更為接近。

關(guān)鍵詞：黏性土坡滑動失穩(wěn)機制； 邊坡地震穩(wěn)定性分析； 有限元強度折減法

0引言

地震邊坡失穩(wěn)機理和邊坡地震穩(wěn)定性分析是構(gòu)成邊坡地震穩(wěn)定性研究的兩大主要內(nèi)容。隨著室內(nèi)外土體測試技術(shù)與土工數(shù)值分析方法的進步，兩者都得到了長足的發(fā)展并形成了各自的理論體系[1-8]。然而，目前對于邊坡地震失穩(wěn)機理與穩(wěn)定性分析的研究少有交叉和融合，將二者綜合考慮用于邊坡穩(wěn)定性研究的理論和方法為數(shù)更少。

地震邊坡失穩(wěn)機理往往偏重于以地震慣性力和孔隙水壓力為根本誘因進行定性分析[1-4]，對坡體應(yīng)力分布及變化特征關(guān)注較少，導(dǎo)致大多數(shù)的地震邊坡失穩(wěn)機理很難與目前邊坡地震穩(wěn)定性分析法中較為成熟的有限元方法相結(jié)合。

以極限平衡法和有限元法為代表的邊坡地震穩(wěn)定性分析方法大都專注于給出邊坡動力安全系數(shù)和動力響應(yīng)特征，而忽視對地震邊坡失穩(wěn)機理的研究[5-8]。因此，探尋合適的地震邊坡失穩(wěn)機理，發(fā)掘其與邊坡地震穩(wěn)定性分析方法間的內(nèi)在聯(lián)系，提出考慮地震邊坡失穩(wěn)機理的邊坡地震穩(wěn)定性分析方法，具有重要的理論意義和應(yīng)用價值。

1反映邊坡土體應(yīng)力變化特征的地震邊坡失穩(wěn)機制

邊坡土體應(yīng)力是地震邊坡失穩(wěn)機理與邊坡地震穩(wěn)定性分析方法的重要聯(lián)系，要將兩者結(jié)合應(yīng)用于分析邊坡地震穩(wěn)定性，則需分析邊坡失穩(wěn)過程中土體應(yīng)力變化特征。

1.1地震作用下黏性土坡的受力特點及滑動模式

由于黏聚力的存在，黏性土坡不會像無黏性土坡一樣沿坡面滑動[9]，在地震作用下其滑動面多表現(xiàn)為上陡下緩，破壞和滑動通常呈現(xiàn)出滑裂面高于坡腳、滑裂面穿過坡腳和滑裂面低于坡腳三種類型(圖1)。

圖1　地震作用下黏性土坡滑動示意圖Fig.1　Sliding diagram of cohesive soil slope under earthquake

與靜力作用下不同，處在地震作用下的土體，其應(yīng)力狀態(tài)是由地震形成的動應(yīng)力和自重形成的靜應(yīng)力的疊加，且動應(yīng)力的大小和方向時刻都在變化。地震作用下的垂直動應(yīng)力σn和水平動應(yīng)力σs定義為[10]：

(1)

(2)

(3)

(4)

式中:ρ為土體密度；vn、vs分別為土體中質(zhì)點的豎向及水平向振動速度；cn為縱波波速；cs為剪切波速；K為體積變形模量；G為剪切模量。

取動應(yīng)力最大值和最不利于穩(wěn)定的方向，則邊坡一點土體應(yīng)力狀態(tài)如圖2所示[10]。一般認為地震以橫向作用為主，故取σs>σn[8]。

圖2　地震作用下土體應(yīng)力狀態(tài)Fig.2　Stress state of soil under earthquake

任意的黏性土坡都能簡化為圖3所示的平面斜坡形式。借鑒王恭先[11]提出的滑坡三段式滑動模式，以斜坡不同位置土體應(yīng)力狀態(tài)及變化來描述斜坡失穩(wěn)狀態(tài)，將地震作用以上述動應(yīng)力的形式引入，按幾何形態(tài)特征將黏性土邊坡劃分為坡底、坡面、坡頂三段，分別分析地震作用下黏性土坡受力特點。

圖3　黏性土坡分段受力示意圖Fig.3　Distribution of the stress in different segments of cohesive soil

如圖3所示，對于三種類型的黏性土滑坡在失穩(wěn)過程中受滑裂面形成及滑動的影響，坡頂段土體的豎直向大主應(yīng)力基本保持不變，水平向小主應(yīng)力減小。在坡頂段土體的擠壓和坡底段土體的阻礙作用下，坡面段土體大致與坡面平行的大主應(yīng)力和與其垂直的小主應(yīng)力均增大。對于前兩種類型的滑坡，滑裂面在坡面出露，坡底段土體只會受到微弱的擠壓作用。而對于第三種類型的滑坡，滑裂面在坡底段表面出露，坡底段土體由于受到坡面段土體和坡頂段土體的推力，致使其分別平行和垂直于滑動面的大主應(yīng)力和小主應(yīng)力均增大。受地震橫豎向拉張作用(σs和σn)影響，坡頂段土體大主應(yīng)力和小主應(yīng)力均被削弱，坡面段土體大主應(yīng)力被削弱，小主應(yīng)力轉(zhuǎn)變?yōu)槔瓚?yīng)力后得到加強，坡底段土體大主應(yīng)力和小主應(yīng)力也都遭到削弱。

1.2地震作用下黏性土坡失穩(wěn)土體應(yīng)力變化特征

為半定量化分析地震作用下黏性土坡失穩(wěn)滑動過程中不同受力段土體的應(yīng)力狀態(tài)變化，以圖3所示的一般平面土坡為分析對象，分別選取位于坡底段、坡面段及坡頂段內(nèi)且埋深均為h，鉛垂向應(yīng)力σz相等的B、M、T三點處的土體(圖4)，分析計算不同受力段土體在地震作用下邊坡失穩(wěn)前后的應(yīng)力狀態(tài)，揭示其剪應(yīng)力變化規(guī)律。

圖4　平面黏性土坡土體受力分析示意圖Fig.4　Schematic diagram of soil stress in plane cohesive soil slope

圖5為邊坡穩(wěn)定時(σ0)和地震作用下開始失穩(wěn)時(σe)不同受力段的應(yīng)力狀態(tài)。受滑裂面形成及滑動的影響，坡頂段T點土體小主應(yīng)力σ3t逐漸減小；受坡頂段的擠壓和坡底段的阻礙作用，坡面段M點土體大主應(yīng)力σ1m增大，小主應(yīng)力σ3m則由壓應(yīng)力逐漸減小并轉(zhuǎn)變?yōu)槔瓚?yīng)力；在坡面段蠕動土體的擠壓下，坡底段B點土體大主應(yīng)力σ1b和小主應(yīng)力σ3b都逐漸增大。受地震橫豎向拉張作用(σs和σn)影響，坡頂段T點土體大主應(yīng)力σ1t和小主應(yīng)力σ3t均被削弱，坡面段M點土體大主應(yīng)力σ1m被削弱，小主應(yīng)力σ3m轉(zhuǎn)變?yōu)槔瓚?yīng)力后得到加強，坡底段B點土體大主應(yīng)力σ1b和小主應(yīng)力σ3b也都遭到削弱。

圖5　地震作用下邊坡失穩(wěn)前后各段土體應(yīng)力變化Fig.5　Variation of soil stress in different segments before and after the slope instability under earthquake

考慮到邊坡土體應(yīng)力狀態(tài)的復(fù)雜性，為方便定量計算，假定土體均質(zhì)且滿足摩爾-庫倫強度準則，將位于坡面附近土體水平向應(yīng)力σx定義為與土體埋深z有關(guān)的函數(shù)[12]:

(5)

式中:K0為側(cè)應(yīng)力系數(shù)；γ為土體重度；β為坡角。

邊坡坡頂段(T點)剪應(yīng)力計算結(jié)果為：

(6)

(7)

式中:t=σ3te/σ3t0，0

邊坡坡面段(M點)剪應(yīng)力計算結(jié)果為

(8)

(9)

式中:m1=σ1me/σ1m0,m1>1；m3=σ3me/σ3m0,0

邊坡坡底段(B點)剪應(yīng)力計算結(jié)果為

(10)

(11)

式中：b1=σ1be/σ1b0,b1>1；b3=σ3be/σ3b0,b3>1；θ為坡底段土體大主應(yīng)力方向與豎直方向的夾角。

根據(jù)以上剪應(yīng)力計算結(jié)果，分別作出三種狀態(tài)時不同受力段土體的摩爾圓(圖6)。從圖中可以看出，與無地震作用下邊坡失穩(wěn)時(σ)相比，地震作用下(σe)不同受力部位土體的剪應(yīng)力受到的影響不同：坡頂段得到較大幅度增長，坡面段呈現(xiàn)小幅增長，坡底段則表現(xiàn)為小幅降低?？傮w來看，地震作用下邊坡失穩(wěn)過程中剪應(yīng)力大小及變化趨勢表現(xiàn)為：坡頂段初始剪應(yīng)力較小但增長迅速；坡面段初始剪應(yīng)力較大且增長較快；坡底段初始剪應(yīng)力較小且增長較慢。

圖6　邊坡失穩(wěn)前后不同滑動段土體剪應(yīng)力變化Fig.6　Shear stress variation of soil in different segments before and after the slope instability

在地震荷載長持時作用下，邊坡土體更容易發(fā)生拉破壞[10,13-14]。上述分析中，坡頂段土體受地震作用影響發(fā)生因小主應(yīng)力由壓應(yīng)力轉(zhuǎn)變?yōu)槔瓚?yīng)力的拉張型剪切破壞，從而較好地解釋了地震邊坡上部往往呈現(xiàn)出拉破壞特征且存在拉張型裂縫的現(xiàn)象。

1.3黏性土坡三段式滑動地震失穩(wěn)機制

邊坡失穩(wěn)滑動是由于其內(nèi)部分土體發(fā)生剪切破壞所致，而這種剪切破壞的發(fā)生與否則取決于土體的剪應(yīng)力大小和變化趨勢。位于邊坡體內(nèi)不同部位的土體，其受力特點和應(yīng)力水平本就不同。在邊坡失穩(wěn)前其剪應(yīng)力水平不同，在邊坡失穩(wěn)過程中，隨著應(yīng)力狀態(tài)變化的差異，不同位置土體剪應(yīng)力的變化趨勢更是不盡相同。但就土體破壞發(fā)展程度而言，剪應(yīng)力在失穩(wěn)前原本就大，且在失穩(wěn)過程中增長快的土體接近破壞的速度更快。

圖6所示的地震作用下不同受力段土體剪應(yīng)力初始水平和變化趨勢表明，坡頂段可能同時甚至先于坡面段發(fā)生破壞。一些震后頂部存在大量深長拉張型裂縫但整體卻未發(fā)生失穩(wěn)的黏性土坡的存在[10,13-14]，也能說明地震作用下此類邊坡的失穩(wěn)可能始于坡頂段的破壞。

綜合考慮地震作用下黏性土坡不同受力段應(yīng)力變化特征和失穩(wěn)破壞特征，將地震作用以應(yīng)力的形式表征，以斜坡不同位置土體應(yīng)力狀態(tài)及變化來描述斜坡地震失穩(wěn)的特征，研究提出黏性土坡三段式滑動地震失穩(wěn)機制：地震作用下，坡頂段受地震拉張作用的影響，側(cè)向支撐力迅速減小從而發(fā)生主動土壓破壞，形成后緣張裂縫;與此同時，坡面段因無法承受坡頂段下滑土體逐漸增大的推力、自身上覆巖土產(chǎn)生的下滑力和地震慣性力而沿原本就已存在的強度較低的相對軟弱帶(面)發(fā)生蠕動;在坡面段和坡頂段剩余下滑力的推擠下，若坡底段土體不發(fā)生破壞，則滑動面在坡面出露，滑體在坡面剪出，若坡底段土體產(chǎn)生被動土壓破壞并形成破壞面，則滑動面空間得以完全貫通三段土體，滑坡開始加速整體滑移。

2考慮黏性土坡三段式滑動地震失穩(wěn)機制的有限元強度折減法

有限元強度折減法是目前較為成熟的邊坡地震穩(wěn)定性分析方法，若要將其與黏性土坡地震三段式滑動失穩(wěn)機制相結(jié)合，則應(yīng)找尋該方法中與剪應(yīng)力有關(guān)的定義并加以分析和利用。

在有限元強度折減法中，邊坡失穩(wěn)是通過土體強度的折減變低導(dǎo)致土體破壞來實現(xiàn)的。折減系數(shù)增幅的大小與土體接近破壞速度的快慢有關(guān)。

既然實際邊坡失穩(wěn)時土體應(yīng)力變化下的剪應(yīng)力初始大小和變化趨勢決定了土體接近破壞的速度，而在有限元強度折減法中土體接近破壞的速度又取決于折減系數(shù)的增長速度，因此可以通過對不同受力段折減系數(shù)進行調(diào)整來實現(xiàn)黏性土坡三段式滑動地震失穩(wěn)機制和有限元強度折減法的結(jié)合。

在運用有限元強度折減法進行邊坡地震穩(wěn)定性分析時，由于不同受力段都在邊坡受力達到極限平衡時停止折減，即折減時間相同，所以要滿足不同受力段土體接近破壞時速度不同的特征，三段土體的折減系數(shù)初值應(yīng)相等，終值應(yīng)保證在使邊坡受力達到極限平衡狀態(tài)的前提下，滿足坡頂段和坡面段較大、坡底段較小的條件。

3數(shù)值模擬驗證

運用上述考慮黏性土坡三段式滑動地震失穩(wěn)機制的有限元強度折減法進行邊坡地震穩(wěn)定性分析，選用定義更為明確合理的斜坡土體位移突變點法，并結(jié)合破裂面貫通及計算收斂兩種失穩(wěn)判據(jù)來確定邊坡動力安全系數(shù)[15-16]。

選用具有集合定義和場變量功能的ABAQUS有限元軟件，利用其集合定義功能完成對邊坡不同受力段的劃分，將坡頂段、坡面段和坡底段折減系數(shù)ωt、ωm和ωb定義為同一場變量?？紤]到坡頂段與坡面段破壞先后順序難以確定，在保證邊坡能夠達到極限平衡狀態(tài)的條件下，分為ωt>ωm>ωb、ωt=ωm>ωb及ωm>ωt>ωb三種情況進行試算。

選取2013年岷縣漳縣MS6.6地震誘發(fā)的永光村西側(cè)滑坡，以驗證拓展的有限元強度折減法的可靠性。該滑坡為地震誘發(fā)的典型黃土滑坡[17-18]，坡高270m，平均坡度為28.7°，斜坡簡化模型如圖7所示。根據(jù)已有的現(xiàn)場測試和室內(nèi)土體力學(xué)實驗結(jié)果，滑坡基巖上覆土體分為三層，各層土體物理力學(xué)參數(shù)如表1所列。土體計算采用理想彈塑性本構(gòu)模型及莫爾-庫侖屈服準則。邊坡的有限元分析模型和網(wǎng)格劃分如圖7所示。模型土體均采用四邊形四節(jié)點平面應(yīng)變單元，網(wǎng)格劃分為763個單元、837個節(jié)點。模型左右邊界采用黏彈性邊界條件，底部為人工邊界。

表 1　各層材料參數(shù)

在動力計算過程中，水平向地震荷載從模型底部由左到右輸入。輸入地震波采用岷縣漳縣6.6級地震中距震中18km的岷縣臺記錄的主震加速度時程(圖8)。其地震峰值加速度為220gal，卓越頻譜為4.5～5.5Hz，輸入地震動作用時間為40s。

圖7　邊坡簡化模型Fig.7　Simplified model of the slope

圖8　地震波加速度時程曲線Fig.8　Acceleration time-history curves of input seismic wave

有限元強度折減法折減系數(shù)和水平向位移計算結(jié)果如圖9所示。折減系數(shù)的取值范圍為0.5～1.2，當折減系數(shù)取1.00時，計算不收斂，邊坡達到極限平衡狀態(tài)，即該邊坡以折減系數(shù)定義的動力安全系數(shù)ωc為1.00。

圖9　有限元強度折減法計算結(jié)果Fig.9　Results of strength reduction FEM

該拓展方法首先應(yīng)保證極限平衡狀態(tài)時各段折減系數(shù)的取值范圍能夠達到原方法所得的折減系數(shù)ωc；其次，由于邊坡滑動時土體中同時存在殘余強度和峰值強度，且前者約為后者的70%～90%[19]，取80%作為極限平衡狀態(tài)時折減系數(shù)變化幅度的上限，即折減系數(shù)最大取值不超過1.25ωc。選取符合上述條件，且滿足ωt> ωm>ωb、ωt=ωm>ωb及ωm>ωt>ωb三種不同條件的三組折減系數(shù)(表2)進行試算。圖10為拓展方法折減系數(shù)和水平向位移的計算結(jié)果，分別畫出兩種方法得到的折減系數(shù)-位移曲線(圖11)。采用以邊坡土體位移突變點法為主的綜合判據(jù)法確定邊坡動力安全系數(shù)，并對兩種方法得到的動力安全系數(shù)加以分析比較。通過對比不同組計算結(jié)果(如表2所列)可以看出：當ωt>ωm>ωb時拓展方法所得動力安全系數(shù)較大；當ωt=ωm>ωb時兩種方法所得動力安全系數(shù)接近；當ωm>ωt>ωb時拓展方法所得動力安全系數(shù)較小。

圖10　拓展方法計算結(jié)果Fig.10　Calculation results using the expanding method

圖11　折減系數(shù)-邊坡頂點水平向位移曲線Fig.11　Relationship between horizontal displacement at the top of slope and reduction factor

該滑坡采用極限平衡法中簡化Bishop法計算得到的動力安全系數(shù)為0.851，采用瑞典條分法及Janbu法得到的動力安全系數(shù)分別為0.891和0.895，拓展后折減系數(shù)滿足ωt>ωm>ωb條件的有限元強度折減法得到的結(jié)果為0.880。該結(jié)果顯然與極限平衡法得出的結(jié)果更為接近，從數(shù)值模擬層面驗證了黏性土坡三段式滑動地震失穩(wěn)機制的合理性。

這說明考慮黏性土坡三段式滑動地震失穩(wěn)機制的動力有限元強度折減法應(yīng)用于黃土斜坡地震穩(wěn)定性分析是合理可行的，以其為依據(jù)確定的動力安全系數(shù)是安全可靠的。

表 2　折減系數(shù)及邊坡安全系數(shù)

4結(jié)論

(1) 邊坡不同受力段土體剪應(yīng)力受到地震作用的影響不同：坡頂段得到較大幅度增長，坡面段呈現(xiàn)出小幅增長，坡底段則表現(xiàn)為小幅降低?？傮w來看，地震作用下邊坡失穩(wěn)過程中剪應(yīng)力大小及變化趨勢表現(xiàn)為：坡頂段初始剪應(yīng)力較小但增長迅速，坡面段初始剪應(yīng)力較大且增長較快，坡底段初始剪應(yīng)力較小且增長較慢。

(2) 結(jié)合黏性土坡滑動模式和地震失穩(wěn)破壞特征，提出了黏性土坡三段式滑動地震失穩(wěn)機制。地震作用下，坡頂段受地震拉張作用的影響，側(cè)向支撐力迅速減小從而發(fā)生主動土壓破壞，形成后緣張裂縫。與此同時，坡面段因無法承受坡頂段下滑土體逐漸增大的推力、自身上覆巖土產(chǎn)生的下滑力和地震慣性力而沿原本就已存在的強度較低的相對軟弱帶(面)發(fā)生蠕動。在坡面段和坡頂段剩余下滑力的推擠下，若坡底段土體不發(fā)生破壞，則滑動面在坡面出露，滑體在坡面剪出，若坡底段土體產(chǎn)生被動土壓破壞并形成破壞面，則滑動面空間得以完全貫通三段土體，滑坡開始加速整體滑移。

(3) 邊坡土體剪應(yīng)力水平及其變化趨勢是地震邊坡失穩(wěn)機理與邊坡地震穩(wěn)定性分析方法的重要聯(lián)系，通過對有限元強度折減法中折減系數(shù)進行調(diào)整，可以充分體現(xiàn)黏性土坡三段式滑動地震失穩(wěn)機制。

(4) 考慮黏性土坡三段式滑動地震失穩(wěn)機制的有限元強度折減法所得的動力安全系數(shù)與極限平衡法所得結(jié)果十分接近。此方法的提出，實現(xiàn)了地震邊坡失穩(wěn)機理與邊坡地震穩(wěn)定性分析方法的初步結(jié)合，為實際工程的黏性土坡地震穩(wěn)定性分析提供了較為可靠的數(shù)值模擬依據(jù)。

參考文獻(References)

[1]KramerSL.EffectsofEarthquakesonDamsandEmbankmentsGeotechnicalEarthquakeEngineering[M].USA,NewJersey:PrenticeHall,1995.

[2]李天池.地震與滑坡的關(guān)系及地震滑坡預(yù)測的探討(節(jié)錄)[C]//滑坡文集(第二集).北京:中國鐵道出版社,1979：127-132.

LITian-chi.DiscussionoftheRelationshipbetweenEarthquakeandtheLandslideandEarthquakeLandslidePrediction(Excerpt)[C]//LandslideAnthology(SecondVolume).Beijing:ChinaRailwayPublishingHouse,1979:127-132.(inChinese)

[3]祁生文,伍法權(quán), 劉春玲,等.地震邊坡穩(wěn)定性的工程地質(zhì)分析[J].巖石力學(xué)與工程學(xué)報,2004,23(16):2792-2797.

QISheng-wen，WUFa-quan,LIUChun-ling,etal.EngineeringGeologyAnalysisonStabilityofSlopeunderEarthquake[J].ChineseJournalofRockMechanicsandEngineering,2004,23(16):2792-2797.(inChinese)

[4]張振中.黃土地震災(zāi)害預(yù)測[M].北京:地震出版社, 1999.

ZHANGZhen-zhong.TheEarthquakeDisasterPredictioninLoessAreas[M].Beijing:SeismologicalPress,1999.(inChinese)

[5]TerzaghiK.MechanismsofLandslides,EngineeringinGeology(Berdey)Volume[M].US:GeologicalSocietyofAmerica,1950.

[6]NewmarkNM.EffectsofEarthquakesonDamsandEmbankments[J].Geotechnique,1965,15(2):139-160.

[7]劉紅帥,薄景山,劉德東.巖土邊坡地震穩(wěn)定性分析研究評述[J].地震工程與工程振動,2005,25(1):164-171.

LIUHong-shuai,BOJing-shan,LIUDe-dong.ReviewonStudyofSeismicStabilityAnalysisofRock-soilSlopes[J].EarthquakeEngineeringandEngineeringVibration,2005,25(1):164-171.(inChinese)

[8]葉海林,黃潤秋,鄭穎人.地震作用下邊坡穩(wěn)定性安全評價的研究[J].地下空間與工程學(xué)報,2009,5(6):1248-1257.

YEHai-lin,HUANGRun-qiu,ZHENGYing-ren.TheStudyofSafetyEvaluationonSlopeSeismicStability[J].ChineseJournalofUndergroundSpaceandEngineering,2009,5(6):1248-1257.(inChinese)

[9]陳仲頤,周景星,王洪瑾.土力學(xué)[M].北京:清華大學(xué)出版社,2004.

CHENZhong-yi,ZHOUJing-xing,WANGHong-jin.SoilMechanics[M].Beijing:TsinghuaUniversityPress,1994.(inChinese)

[10]言志信,張森,張學(xué)東,等.地震邊坡失穩(wěn)機理及穩(wěn)定性分析[J].工程地質(zhì)學(xué)報,2010,18(6):844-850.

YANZhi-xin,ZHANGSen,ZHANGXue-dong,etal.FailureMechanismandStabilityAnalysisofSlopeunderEarthquake[J].JournalofEngineeringGeology,2010,18(6):884-850.(inChinese)

[11]王恭先.滑坡學(xué)與滑坡防治技術(shù)[M].北京:中國鐵道出版社,2004.

WANGGong-xian.LandslideandLandslideControlTechnology[M].Beijing:ChinaRailwayPublishingHouse,2004.(inChinese)

[12]王國體.以土體應(yīng)力狀態(tài)計算邊坡安全系數(shù)的方法[J].中國工程科學(xué),2006,8(12):80-84.

WANGGuo-ti.TheMethodandFormulaforStabilityFactorofSlopewithStateofSoilOriginalStress[J].EngineeringScience,2006,8(12):80-84.(inChinese)

[13]鄭穎人,葉海林,黃潤秋.地震邊坡破壞機制及其破裂面的分析探討[J].巖石力學(xué)與工程學(xué)報,2009,28(8):1714-1723.

ZHENGYing-ren,YEHai-lin,HUANGRun-qiu.AnalysisandDiscussionofFailureMechanismandFractureSurfaceofSlopeunderEarthquake[J].ChineseJournalofRockMechanicsandEngineering,2009,28(8):1714-1723.(inChinese)

[14]鐘秀梅,王謙,張國信,等.天水市水眼寨滑坡特征及穩(wěn)定性分析[J].地震工程學(xué)報,2014,36(4):887-891.

ZHONGXiu-mei,WANGQian,ZHANGGuo-xin,etal.CharacteristicsandStabilityAnalysisofShuiyanzhaiLandslideinTianshuiCity[J].ChinaEarthquakeEngineeringJournal,2014,36(4):887-891.(inChinese)

[15]鄭穎人,趙尚毅.有限元強度折減法在土坡與巖坡中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報,2004,23(19):3381-3388.

ZHENGYing-ren,ZHAOShang-yi.ApplicationofStrengthReductionFEMinSoilandRockSlope[J].ChineseJournalofRockMechanicsandEngineering,2004,23(19):3381-3388.(inChinese)

[16]雷天,吳志堅,陳拓.天水市某典型黃土邊坡塊石支護穩(wěn)定性效果分析[J].地震工程工程學(xué)報,2014,36(4):968-974.

LEITian,WUZhi-jian,CHENTuo.StabilityAnalysisofaTypicalLoessSlopeSupportedbyBlockStoneRetainingWallinTianshui[J].ChinaEarthquakeEngineeringJournal,2014,36(4):968-974.(inChinese)

[17]WUZhi-jian,SUNJun-jie,CHENYu-jin,etal.AnalysisofDisaster-causingMechanismofLoessLandslidesInducedbytheMinxian-ZhangxianMS6.6Earthquake,China[C]//Proceedingof6thChina-JapanGeotechnicalSymposium,Hokkaido,2015.

[18]車愛蘭,吳志堅,彭冬,等.黃土斜坡震害面波勘探調(diào)查及其動力穩(wěn)定性分析——以甘肅岷縣、漳縣MS6.6地震為例[J].地震工程學(xué)報,2013,35(4):724-729.

CHEAi-lan,WUZhi-jian,PENGDong,etal.SurfaceWaveInvestigationandDynamicStabilityAnalysisforEarthquake-inducedLoessLandslides[J].ChinaEarthquakeEngineeringJournal,2013,35(4):724-729.(inChinese)

[19]譚文輝,任奮華,苗勝軍.峰值強度與殘余強度對邊坡加固的影響研究[J].巖土力學(xué),2007,28(增刊):616-618.

TANWen-hui,RENFen-hua,MIAOSheng-jun,etal.InfluenceofParametersofPeakStrengthandResidualStrengthontheReinforcementofSlopes[J].RockandSoilMechanics,2007,28(Supp):616-618.(inChinese)

EarthquakeInstabilityMechanismofCohesiveSoilSlopeandItsApplicationinSeismicStabilityAnalysisofLoessSlope

ZHAOWen-chen1,WUZhi-jian1,2,CHENYu-jin1

(1.Key Laboratory of Loess Earthquake Engineering, Lanzhou Institute of Seismology, CEA, Lanzhou 730000, Gansu, China;2. Geotechnical Disaster Prevention Engineering Technology Research Center of Gansu Province, Lanzhou 730000, Gansu, China)

Abstract：Based on the earthquake instability characteristics of cohesive soil slopes, in this study, we used the stress analysis method to analyze the process of slope instability. By analyzing the stress state and shear stress variation features of soil in different parts of the slope during earthquakes, combined with the failure characteristics of actual slopes during earthquakes, we developed our proposed three-segment sliding instability mechanism of cohesive soil slopes. From our analysis results of the stress correlation between the above-mentioned mechanism combined with the strength reduction finite element method (FEM), we carried out a dynamic stability analysis of an actual loess seismic landslide. A comparison of our results with those of other methods reveals that the dynamic safety factor obtained by the strength reduction FEM, which considers the three-segment sliding instability mechanism of cohesive soil slopes, was very close to that determined by the limit equilibrium method.

Key words：sliding instability mechanism of cohesive soil slope; seismic stability analysis of slope; strength reduction FEM

收稿日期：①2015-12-31

基金項目：國家自然科學(xué)基金資助項目(41472297，51478444)；青海省交通建設(shè)科技項目“鐵路建設(shè)對雅丹地貌的環(huán)境影響與減輕技術(shù)研究”

作者簡介：趙文琛(1991-)，男，碩士，主要從事巖土地震工程方面的研究。E-mail：wenchenzhao1219@126.com。 通信作者：吳志堅(1974-)，男，安徽安慶人，博士，研究員，主要從事巖土地震工程與凍土工程領(lǐng)域的研究。E-mail：zhijianlz@163.com。

中圖分類號:TD853.34

文獻標志碼：A

文章編號:1000-0844(2016)03-0452-08

DOI:10.3969/j.issn.1000-0844.2016.03.0452