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

?

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

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

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

本文研究了在震源近場(chǎng)時(shí)，地震波斜入射條件下，不同坡面形態(tài)對(duì)碎石土邊坡穩(wěn)定性的影響。利用FLAC3D有限差分軟件，建立了多種坡面形態(tài)的邊坡模型。采用對(duì)比方法，分析了邊坡在地震波作用下的穩(wěn)定性。結(jié)果表明：臨空坡面對(duì)輸入地震波有放大效應(yīng)；坡面的水平位移與加速度峰值隨高度增加而增大，隨坡度增加而增大；隨著坡面凹度的增大，剪應(yīng)力增量向坡肩處集中，坡肩愈加不穩(wěn)定；隨著坡面凸度的增大，剪應(yīng)力增量向坡腰處集中，坡腰愈加不穩(wěn)定；臺(tái)階對(duì)坡面的加速度放大效應(yīng)有削弱作用，臺(tái)階數(shù)越多越穩(wěn)定。同時(shí)給出了邊坡防治的幾點(diǎn)建議：減小坡面的坡度；對(duì)凹型坡加強(qiáng)坡肩位置處的防治，對(duì)凸型坡加強(qiáng)坡腰位置處的防治；對(duì)高邊坡可設(shè)置多級(jí)臺(tái)階增加穩(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ā)生時(shí)極可能再次滑坡（黃潤(rùn)秋，2011）。因此，對(duì)碎石土邊坡穩(wěn)定性的研究具有重要意義。

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

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

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

1 方案分析

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

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

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

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

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

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

（2）

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

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

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

（4）

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

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

（6）

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

（8）

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

（10）

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

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

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

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

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

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

2 算例驗(yàn)證及分析

2.1 算例模型

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

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

圖2 計(jì)算模型及監(jiān)測(cè)點(diǎn)圖

圖3 輸入位移圖

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

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

圖4 SV波0°垂直入射下A點(diǎn)位移時(shí)程

圖5 SV波0°垂直入射下B點(diǎn)位移時(shí)程

圖6 SV波30°垂直入射下A點(diǎn)位移時(shí)程

圖7 SV波30°垂直入射下B點(diǎn)位移時(shí)程

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

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

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

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

（1）坡度影響組

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

圖8 坡度影響組模型尺寸及檢測(cè)點(diǎn)

（2）凹凸型影響組

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

圖9 凹凸度影響組模型尺寸及檢測(cè)點(diǎn)

（3）臺(tái)階級(jí)數(shù)影響組

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

3.2 模型材料

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

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

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

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

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

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

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

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

4.1 坡度影響組

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

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

圖13 各監(jiān)測(cè)點(diǎn)水平位移

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

由圖13可知，3個(gè)碎石土邊坡坡面各監(jiān)測(cè)點(diǎ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的碎石土邊坡相對(duì)較為穩(wěn)定；坡度比為1:1.3的碎石土邊坡局部滑動(dòng)，但整體還未破壞；坡度比為1:1的碎石土邊坡大范圍滑動(dòng)，邊坡整體破壞。從碎石土邊坡各監(jiān)測(cè)點(diǎn)水平位移的大小可以得出，坡面的水平位移隨高度增加而增大，隨坡度增加而增大，這與坡面的PGA放大系數(shù)是相吻合的。

4.2 凹凸度影響組

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

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

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

圖16 凹型坡面坡肩A點(diǎn)位移

圖17 凸型坡面坡腰B點(diǎn)位移

4.3 臺(tái)階級(jí)數(shù)影響組

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

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

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

5 結(jié)論

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

張如林，樓夢(mèng)麟，2010．基于FLAC3D的斜入射地震波作用的數(shù)值模擬方法研究．土木工程學(xué)報(bào)，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)練項(xiàng)目（編號(hào)：040606693）

2014-11-21

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

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