斜入射地震作用下坡面形態(tài)對碎石土邊坡穩(wěn)定性影響研究1

邱 敏 陳 行 田明杰 龍雪梅 鄒祖銀

?

斜入射地震作用下坡面形態(tài)對碎石土邊坡穩(wěn)定性影響研究1

邱 敏1）陳 行1）田明杰1）龍雪梅1）鄒祖銀1，2）

1）四川農(nóng)業(yè)大學(xué)土木工程學(xué)院，都江堰611830?2）西南交通大學(xué)土木工程學(xué)院，成都610031

本文研究了在震源近場時，地震波斜入射條件下，不同坡面形態(tài)對碎石土邊坡穩(wěn)定性的影響。利用FLAC3D有限差分軟件，建立了多種坡面形態(tài)的邊坡模型。采用對比方法，分析了邊坡在地震波作用下的穩(wěn)定性。結(jié)果表明：臨空坡面對輸入地震波有放大效應(yīng)；坡面的水平位移與加速度峰值隨高度增加而增大，隨坡度增加而增大；隨著坡面凹度的增大，剪應(yīng)力增量向坡肩處集中，坡肩愈加不穩(wěn)定；隨著坡面凸度的增大，剪應(yīng)力增量向坡腰處集中，坡腰愈加不穩(wěn)定；臺階對坡面的加速度放大效應(yīng)有削弱作用，臺階數(shù)越多越穩(wěn)定。同時給出了邊坡防治的幾點建議：減小坡面的坡度；對凹型坡加強坡肩位置處的防治，對凸型坡加強坡腰位置處的防治；對高邊坡可設(shè)置多級臺階增加穩(wěn)定性。

坡面形態(tài) 地震 斜入射 碎石土 穩(wěn)定性

引言

近年來我國的地震頻發(fā)，在地震力的作用下極易觸發(fā)邊坡滑坡。在地震中由于滑坡造成的次生災(zāi)害十分常見，尤其是大型和巨型滑坡造成的災(zāi)害非常嚴(yán)重。在發(fā)生過地震的災(zāi)區(qū)存在著數(shù)量眾多的碎石土堆積體邊坡，其坡面形態(tài)多種多樣，多數(shù)處于亞穩(wěn)定狀態(tài)，當(dāng)余震發(fā)生時極可能再次滑坡（黃潤秋，2011）。因此，對碎石土邊坡穩(wěn)定性的研究具有重要意義。

地震波輸入的入射角對邊坡的動力響應(yīng)影響較大（鄧龍勝等，2012）。一般在計算邊坡的地震反應(yīng)時，地震波輸入通常采用垂直向上入射的剪切波方式，這對于震源距邊坡較遠(yuǎn)時是合理的；然而當(dāng)震源距邊坡較近時，地震波并非垂直入射，而是以一定的傾斜角度向上傳播，這時地震動呈現(xiàn)出空間變化特性（苑舉衛(wèi)等，2010），其中一個有效的證據(jù)是利用強震記錄可統(tǒng)計出低頻段S波的入射角均值為56o左右（尤紅兵等，2009）。

目前，很多學(xué)者對地震作用下邊坡的穩(wěn)定性進(jìn)行了相關(guān)研究。如：黃博等（2013）基于半無限彈性空間地震波的傳播理論，對斜入射地震波在土體中的動應(yīng)力路徑進(jìn)行了分析；李山有等（2005）應(yīng)用顯示有限元法，對傾斜斷層臺階在地震體波斜入射下的地震反應(yīng)進(jìn)行了數(shù)值模擬；何劉等（2014）采用自制的單向振動試驗臺，對不同坡面形態(tài)的土質(zhì)邊坡失穩(wěn)破壞過程進(jìn)行了模型試驗研究；盧坤林等（2014）從失穩(wěn)邊坡資料統(tǒng)計分析、模型試驗和理論計算三個方面，研究了坡面形態(tài)對邊坡穩(wěn)定的影響。但是，對斜入射地震波作用下的碎石土邊坡以及坡面形態(tài)對邊坡穩(wěn)定性影響的研究還較少。

本文在前人的大量研究基礎(chǔ)之上，基于有限元差分軟件FLAC3D展開了多種坡面形態(tài)的對比分析，其中包括不同坡度坡面的對比，不同凹凸度坡面的對比和不同臺階數(shù)坡面的對比。本文得出的坡面形態(tài)對邊坡穩(wěn)定性的影響規(guī)律及關(guān)系，可為碎石土邊坡在斜入射地震波作用下的穩(wěn)定性研究提供基礎(chǔ)數(shù)據(jù)，也為類似邊坡的防治提供了相關(guān)的技術(shù)依據(jù)。

1 方案分析

1.1 數(shù)值模擬理論基礎(chǔ)

FLAC3D采用完全非線性分析方法并基于顯式差分法，使用由周圍區(qū)域真實密度得出的網(wǎng)格節(jié)點集中質(zhì)量，求解全部運動方程。采用摩爾-庫倫模型，程序可以自動計算永久變形，并同時模擬壓縮波和剪切波的傳播及兩者耦合作用時對材料的影響（陳育民等，2009）。因此，F(xiàn)LAC3D可以模擬邊坡在地震動作用下的完全非線性響應(yīng)。

1.2 斜入射地震作用的實現(xiàn)

理論上，任何一種地震波的入射方式都可以分解為垂直和水平兩個方向上的震動（陳曉利等，2011）。首先根據(jù)波動理論，將靜態(tài)邊界上兩個方向的波動場分別分解為自由波場（已分解得到的兩個方向的震動波）和散射波場（由局部場地效應(yīng)引起）兩部分（廖振鵬，2002）。然后，將地震波動的問題轉(zhuǎn)化為波源的問題，即將地震波轉(zhuǎn)化為直接作用在靜態(tài)邊界上法向和切向的等效荷載。

基于人工邊界等效荷載法（張如林等，2010），總波場u可分解為自由波場u和散射波場u：

由位移連續(xù)條件和力學(xué)平衡條件，靜態(tài)人工邊界上任意一點B的運動方程可寫為：

（2）

式中，BT、BN為質(zhì)量矩陣；BT、BN分別為靜態(tài)人工邊界節(jié)點上平面內(nèi)切向和法向的阻尼器系數(shù)；BT、BN為剛度矩陣；、分別為模擬自由波場需施加在靜態(tài)人工邊界節(jié)點上的切向和法向荷載；、分別為模擬散射波場需施加在靜態(tài)人工邊界點上的切向和法向荷載，可由靜態(tài)人工邊界提供，故無需在邊界節(jié)點上額外施加。

自由波場的模擬，需要首先設(shè)定已知的原自由波場為0(，，)和0(，，)。其中，在靜態(tài)人工邊界上B點產(chǎn)生的切向位移和法向位移分別為0(B，B，)和0(B，B，)；切向應(yīng)力和法向應(yīng)力分別為0(B，B，)和0(B，B，)，B和B為靜態(tài)人工邊界上節(jié)點B的坐標(biāo)。

為準(zhǔn)確實現(xiàn)波動輸入，還需要在靜態(tài)人工邊界上施加連續(xù)分布的等效荷載，以使靜態(tài)人工邊界上產(chǎn)生的位移和應(yīng)力與原自由波場中的位移與應(yīng)力相同，即滿足以下條件：

（4）

式中，(B，B，)和(B，B，)、(B，B，)和(B，B，)分別為施加的等效荷載在靜態(tài)人工邊界上B點產(chǎn)生的位移和應(yīng)力。

如圖1所示，根據(jù)一般力學(xué)分析中的脫離體概念，可將靜態(tài)人工邊界與其附加的阻尼元件脫離，則靜態(tài)人工邊界上B點的應(yīng)力為：

（6）

將式（3）代入到式（6）可得：

（8）

聯(lián)立式（4）、式（7）和式（8），根據(jù)已知的自由波場求出人工邊界上的速度和應(yīng)力，進(jìn)而分別計算出靜態(tài)人工邊界上B點需要施加的切向和法向等效荷載：

（10）

這樣便實現(xiàn)了地震波轉(zhuǎn)化為直接作用在靜態(tài)邊界上法向和切向的等效荷載，就可以在FLAC3D里實現(xiàn)地震波斜入射條件下碎石土邊坡動力響應(yīng)的數(shù)值模擬。

1.3 坡面形態(tài)的選擇

碎石土邊坡的坡面形態(tài)復(fù)雜多樣，自然風(fēng)化形成的邊坡和地震波作用后形成的邊坡均可分為四種基本的類型：凸面、凹面、順直面和臺階面，以及它們的相互組合形成的組合形式。為了便于研究不同坡面形態(tài)對斜入射地震波作用下碎石土邊坡穩(wěn)定性的影響，筆者選擇了以下簡化對比分析：

（1）順直型坡面的坡度對碎石土邊坡穩(wěn)定性的影響；

（2）順直型（即不凹不凸）、凸面型和凹面型坡面的凹凸度，對碎石土邊坡穩(wěn)定性的影響；

（3）臺階型坡面的臺階數(shù)對碎石土邊坡穩(wěn)定性的影響。

2 算例驗證及分析

2.1 算例模型

為驗證前面介紹的在FLAC3D中實現(xiàn)斜入射地震波作用的數(shù)值模擬方法，本文的算例模型采用二維計算模型，如圖2所示。算例中模型的各種介質(zhì)參數(shù)為：剪切模量為16.7MPa；彈性模量為40MPa；泊松比為0.2；容重為2000kg/m3。

在地震反應(yīng)分析時，輸入的地震波通常采用的是SV波，這里僅進(jìn)行SV波輸入下的驗證。入射波采用正弦波位移，位移波時程如圖3所示，幅值為1m，持續(xù)時間為2s，波速s為1400m/s，分別進(jìn)行0°垂直向上入射和30°斜入射。在計算模型中選取兩個監(jiān)測點A和B，其中A為自由表面節(jié)點，B為模型底部節(jié)點，如圖2所示。

圖2 計算模型及監(jiān)測點圖

圖3 輸入位移圖

2.2 算例結(jié)果及分析

在入射角分別為0°和30°的SV波作用下，通過FLAC3D計算可得到位移數(shù)值解，根據(jù)地震波反射原理推導(dǎo)出位移理論解，如圖4—圖7所示。

圖4 SV波0°垂直入射下A點位移時程

圖5 SV波0°垂直入射下B點位移時程

圖6 SV波30°垂直入射下A點位移時程

圖7 SV波30°垂直入射下B點位移時程

從圖4可以看出，SV波在0°垂直入射時，由于斜坡自由表面的反射作用，入射波和反射波在自由表面處幅值疊加，自由表面節(jié)點A的位移幅值為入射位移幅值的1.5倍，數(shù)值解和理論解符合很好。從圖5可以看出，在底部節(jié)點B的位移時程中可看到兩個波，其中第一個周期的正弦波為入射波，第二個周期的正弦波為從自由表面反射回來的波，它們在到達(dá)底部之后的時間內(nèi)，其位移值幾乎維持在0附近，即在底面靜態(tài)人工邊界上能夠全部透射出去。這表明本文采用的靜態(tài)人工邊界和等效荷載輸入方法具有較高的計算精度。從圖6和圖7可以看出，SV波在30°斜入射下恰好垂直于斜坡，入射波和反射波在自由表面處幅值疊加，自由表面節(jié)點A的位移幅值為入射位移幅值的1.73倍，底部節(jié)點B的位移幅值為入射位移幅值的0.87倍。無論是0°垂直入射，還是30°斜入射，表面A點和底部B點的數(shù)值計算結(jié)果與理論解都非常符合，這表明采用人工邊界和等效荷載波動輸入的方式，可以進(jìn)行碎石土邊坡在斜入射地震波作用下的分析。

3 數(shù)值模型的建立

3.1 模型尺寸及監(jiān)測點布置

模型基本尺寸為：長250m；寬60m；高130m。為便于對比分析，可將各種坡面形態(tài)的邊坡模型進(jìn)行以下分組。

（1）坡度影響組

建立3個不同坡面坡度的模型，其坡度比分別為：1:1.73、1:1.3和1:1。在3個不同坡度模型的坡面上分別布置坡肩A點、坡腰B點和坡腳C點，如圖8所示。

圖8 坡度影響組模型尺寸及檢測點

（2）凹凸型影響組

建立5個不同坡面凹凸度的模型，其分別為：凹型2個、不凹不凸型（即順直型）1個和凸型2個。在5個不同凹凸度模型的坡面上分別布置坡肩A點、坡腰B點和坡腳C點，如圖9所示。

圖9 凹凸度影響組模型尺寸及檢測點

（3）臺階級數(shù)影響組

建立3個不同坡面臺階數(shù)的模型，其分別為：無臺階（即順直型）、一級臺階和兩級臺階。在無臺階模型的坡面上分別布置坡肩A1點、坡腰B1點和坡腳C1點；在一級臺階模型的坡面上分別布置坡肩A2點、坡腰B2點和坡腳C2點；在兩級臺階模型的坡面上分別布置坡肩A3點、坡腰B3點和坡腳C3點，如圖10所示。

3.2 模型材料

本文的重點是探究斜入射地震波作用下坡面形態(tài)對碎石土邊坡變形的影響，因此，假定邊坡為均質(zhì)碎石土邊坡，采用吳周明等（2012）給出的物理力學(xué)參數(shù)取值，如表1所示。

表1 碎石土的物理力學(xué)參數(shù)表

3.3 模型邊界條件與阻尼設(shè)置

在模型四周設(shè)置自由場邊界，對生成的自由場網(wǎng)格進(jìn)行g(shù)roup賦值，保證在人工邊界處的傳播特性與在原連續(xù)介質(zhì)中保持一致。這樣可較好地模擬無限域的效果，使向上的面波在邊界上不會發(fā)生扭曲。在模型底部設(shè)置靜態(tài)邊界，但輸入地震波作用時需將地震波轉(zhuǎn)換為人工邊界上的等效荷載輸入。

在FLAC3D中有瑞利阻尼、局部阻尼和滯后阻尼。考慮到瑞利阻尼的理論與常規(guī)動力分析方法類似，而且實踐證明瑞利阻尼計算得到的加速度響應(yīng)規(guī)律比較符合實際。因此，本文采用瑞利阻尼，最小臨界阻尼比min取5%，最小中心頻率min取16.8。

3.4 地震波的選擇與轉(zhuǎn)換

本文選取汶川地震波的SV波，圖11為汶川地震波加速度時程，其加速度峰值為0.3g，持時為20s，入射角為56o。按照上述方案分析中的斜入射實現(xiàn)方法，利用數(shù)學(xué)軟件MATLAB計算出汶川地震波的等效荷載，然后編成table文件，再輸入到FLAC3D中，從而實現(xiàn)地震波的斜入射。

4 各組模型的計算及結(jié)果分析

4.1 坡度影響組

坡度影響組分別計算和記錄了3種不同坡度的碎石土邊坡沿坡面上各監(jiān)測點的加速度、速度和位移。為了描述地震動力作用下碎石土邊坡坡面加速度的反映規(guī)律，筆者設(shè)定坡面各監(jiān)測點的加速度峰值PGA與模型底部輸入地震波的加速度峰值的比值為PGA放大系數(shù)（徐光興等，2008）。根據(jù)各監(jiān)測點觀測到的數(shù)據(jù)，處理后可得到不同坡度的碎石土邊坡各監(jiān)測點水平方向PGA放大系數(shù)關(guān)系圖（見圖12）和各監(jiān)測點水平位移圖（見圖13）。

圖12 各監(jiān)測點水平方向PGA放大系數(shù)

圖13 各監(jiān)測點水平位移

由圖12可知，通過對比坡面各高度監(jiān)測點水平方向的PGA放大系數(shù)后發(fā)現(xiàn)，水平方向PGA放大系數(shù)呈現(xiàn)出沿坡面隨高度增加而增大的規(guī)律；通過對比不同坡度的碎石土邊坡坡面的各監(jiān)測點水平方向PGA放大系數(shù)后發(fā)現(xiàn)，PGA放大系數(shù)呈現(xiàn)出隨坡度增加而增大的規(guī)律。

由圖13可知，3個碎石土邊坡坡面各監(jiān)測點均有不同大小的水平位移，坡度比為1:1.73 的碎石土邊坡總體水平位移在0.25—0.78m范圍內(nèi)；坡度比為1:1.3的碎石土邊坡總體水平位移在0.5—1.75m范圍內(nèi)；坡度比為1:1的碎石土邊坡總體水平位移在0.82—2.89m范圍內(nèi)。這表明在斜入射地震波的作用下碎石土邊坡均不穩(wěn)定，但不同坡度的碎石土邊坡的不穩(wěn)定程度是有差異的，它們分別表現(xiàn)為：坡度比為1:1.73的碎石土邊坡相對較為穩(wěn)定；坡度比為1:1.3的碎石土邊坡局部滑動，但整體還未破壞；坡度比為1:1的碎石土邊坡大范圍滑動，邊坡整體破壞。從碎石土邊坡各監(jiān)測點水平位移的大小可以得出，坡面的水平位移隨高度增加而增大，隨坡度增加而增大，這與坡面的PGA放大系數(shù)是相吻合的。

4.2 凹凸度影響組

凹凸度影響組分別計算和記錄了5種不同凹凸度的碎石土邊坡沿坡面上各監(jiān)測點的加速度、速度和位移，并繪制出了各邊坡的剪應(yīng)變增量云圖，如圖14和圖15所示，其中，凹型2>凹型1；凸型2>凸型1。

由圖14可知，剪應(yīng)變增量隨邊坡凹度增加，從順直型的近似均布逐步沿坡面向上轉(zhuǎn)移到坡肩集中分布，可以推斷凹型坡的豎向位移主要集中在坡肩處。由圖15可知，剪應(yīng)變增量隨邊坡凸度增加，從順直型的近似均布逐步沿坡面向中下轉(zhuǎn)移到坡腰集中分布，可以推斷凸型坡的水平位移主要集中在坡腰處。因此，凹型邊坡和凸型邊坡的破壞形式是不同的。

單獨對比凹型1和凹型2邊坡坡肩的豎向位移后發(fā)現(xiàn)，凹型2的終時豎向位移為-2.65m，大于凹型1的終時豎向位移-1.20m，如圖16所示。單獨對比凸型1和凸型2邊坡坡腰的水平位移后發(fā)現(xiàn)，凸型2的終時水平位移為2.35m，大于凸型1的終時水平位移1.20m，如圖17所示。由此可以得出，凹型邊坡的坡肩是潛在滑坡面，凸型邊坡的坡腰下方是潛在滑坡面，并且凹度越大，坡肩的豎向位移越大，邊坡越不穩(wěn)定；凸度越大，坡腰的水平位移越大，邊坡越不穩(wěn)定。

圖16 凹型坡面坡肩A點位移

圖17 凸型坡面坡腰B點位移

4.3 臺階級數(shù)影響組

臺階級數(shù)影響組分別計算和記錄了3種不同臺階級數(shù)的邊坡沿坡面上各監(jiān)測點的加速度、速度和位移，并繪制出了各監(jiān)測點位移-時間關(guān)系圖，如圖18所示。

由圖18可知，被臺階分割形成的各小坡段內(nèi)，坡面位移和整坡的坡面位移變化規(guī)律大致相同，沿坡面隨高度增加而增大；坡度和總高度相同但臺階數(shù)不同的3個邊坡，它們各自的總體豎向位移有所差異，其中無臺階邊坡的總體豎向位移最大，在-01.25—0.12m范圍內(nèi)；一級臺階邊坡的總體豎向位移次之，在-0.96—0.07m范圍內(nèi)；兩級臺階邊坡的總體豎向位移最小，在-00.65—0.08m范圍內(nèi)。這表明坡度和高度相同的碎石土邊坡，在斜入射地震波的作用下坡面的臺階級數(shù)越多邊坡的穩(wěn)定性越好。

綜上筆者認(rèn)為，加速度在沿坡面?zhèn)鬟f時，臺階對邊坡坡面的加速度峰值放大效應(yīng)有減弱的作用，從而減小了位移量。同時，臺階將一個長又高的坡面分割成了幾個小坡面，降低了坡面的高度，從而減少了坡面加速度峰值的放大效應(yīng)，繼而減小了位移量。

5 結(jié)論

本文利用FLAC3D軟件進(jìn)行了各種坡面形態(tài)的碎石土邊坡在斜入射地震波作用下的動力響應(yīng)模擬實驗，通過大量的數(shù)值計算和不同模型的對比，對斜入射地震波作用下的碎石土邊坡得到了如下認(rèn)識：

（1）各種坡面形態(tài)的碎石土邊坡均有較大的位移，表明碎石土邊坡在斜入射地震波作用下大多是不穩(wěn)定的，存在滑坡的安全隱患，相關(guān)部門需重視類似碎石土邊坡的防治工作。

（2）在斜入射的情況下，坡面凹凸度對碎石土邊坡潛在滑坡面的位置有明顯影響，且與一般垂直入射的影響不同：隨著坡面凹度的增大，剪應(yīng)力增量向坡肩處集中，坡肩愈加不穩(wěn)定；隨著坡面凸度的增大，剪應(yīng)力增量向坡腰處集中，坡腰愈加不穩(wěn)定。臺階對坡面的加速度放大效應(yīng)有削弱作用，臺階數(shù)越多越穩(wěn)定。

（3）邊坡防治的幾點建議：減小坡面的坡度；對凹型坡加強坡肩位置處的防治，對凸型坡加強坡腰位置處的防治；對高邊坡可設(shè)置多級臺階增加穩(wěn)定性。

陳曉利，李楊，洪啟宇等，2011．地震作用下邊坡動力響應(yīng)的數(shù)值模擬研究．巖石學(xué)報，27（6）：1899—1908．

陳育民，徐鼎平，2009．FLAC/FLAC3D基礎(chǔ)與工程實例．北京：中國水利水電出版社，181—184.

鄧龍勝，范文，2012．黃土邊坡動力響應(yīng)的影響效應(yīng)研究．工程地質(zhì)學(xué)報，20（4）：483—490．

何劉，吳光，趙志明，2014．坡面形態(tài)對邊坡動力變形破壞影響的模型試驗研究．巖土力學(xué)，33（1）：111—117．

黃博，凌道盛，丁浩，2013．斜入射地震波在土體中產(chǎn)生的動應(yīng)力路徑及試驗?zāi)M．巖土工程學(xué)報，35（2）：276—283．

黃潤秋，2011．汶川地震地質(zhì)災(zāi)害后效應(yīng)分析．工程地質(zhì)學(xué)報，19（2）：145—151．

李山有，馬強，韋慶海，2005. 地震體波斜入射下的斷層臺階地震反應(yīng)分析．地震研究，28（3）：277—281．

廖振鵬，2002．工程波動理論導(dǎo)引．北京：科學(xué)出版社，81—92．

盧坤林，朱大勇，2014．坡面形態(tài)對邊坡穩(wěn)定性影響的理論與試驗研究．巖石力學(xué)與工程學(xué)報，33（1）：35—42．

吳周明，傅軍鍵，2012．碎石土邊坡穩(wěn)定性及影響因素的數(shù)值模擬研究．江西理工大學(xué)學(xué)報，33（1）：30—34．

徐光興，姚令侃，李朝紅等，2008．邊坡地震動力響應(yīng)規(guī)律及地震動參數(shù)影響研究．巖土工程學(xué)報，30（6）：918—923．

尤紅兵，趙鳳新，榮棉水，2009．地震波斜入射時水平層狀場地的非線性地震反應(yīng)．巖土工程學(xué)報，31（2）：234—240．

苑舉衛(wèi)，杜成斌，劉志明，2010．基于設(shè)計地震動的地震波斜入射波動輸入研究．四川大學(xué)學(xué)報（工程科學(xué)版），42（5）：250—255．

張如林，樓夢麟，2010．基于FLAC3D的斜入射地震波作用的數(shù)值模擬方法研究．土木工程學(xué)報，43（S1）： 22—27．

Effect of Slope Shape on Gravel-Soil Slope Stability in Obliquely Incident Seismic Wave

Qiu Min1）, Chen Hang1）, Tian Mingjie1）, Long Xuemei1）and Zou Zuyin1,2）

1）College of Civil Engineering, Sichuan Agricultural University, Dujiangyan 611830, China?2）College of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China

In this paper we study the different effect of slope shape on gravel-soil slope stability under condition of obliquely incident seismic wave at the source near field. We use the FLAC3D to set up slope model with a variety of slope shape. We also use the comparison method to analyze the stability of slope under seismic action. Our results showed that faced slope has amplification effect on entered earthquake waves; the level displacement and acceleration peak of slope surface increase with the slope height and gradient increase. With increase of slope surface concave degree, shear stress increment concentrates on slope shoulder, slope shoulder become not stable; with increase of slope surface convex degree, shear stress increment concentrates on slope waist, slope waist become not stable. Slope steps can weaken the acceleration amplification effect of slope surface, the more the step number the slope has, the more stable it is. Finally, we suggest the following countermeasures to increase the stability of slope in the future, i.e. reducing the slope gradient; strengthening the concave slope; and setting more steps.

Slope shape; Earthquake; Obliquely incident; Gravel-soil; Stability

四川農(nóng)業(yè)大學(xué)創(chuàng)新訓(xùn)練項目（編號：040606693）

2014-11-21

邱敏，男，生于1992年。四川農(nóng)業(yè)大學(xué)土木工程學(xué)院本科生，土木工程專業(yè)。E-mail：qiuminqm@foxmail.com

鄒祖銀，男，生于1974年。副教授。主要從事道路與鐵路工程的科研和教學(xué)工作。E-mail：zouzyin@gmail.com