劉超群，傅 偉，郭 威

(1.江西省高速公路投資集團有限責任公司，江西 南昌 330025;2.江西省高速公路投資集團有限責任公司撫州管理中心，江西撫州 330086;3.江西省交通科學研究院，江西 南昌 330200)

0 引言

對于高填方路堤邊坡等由于地下水位變化、降雨等造成路基土體的干濕交替頻繁，路基土在干濕循環(huán)作用下強度減弱，路堤邊坡安全系數(shù)降低，導致滑坡的現(xiàn)象發(fā)生，干濕循環(huán)是誘發(fā)失穩(wěn)邊坡的重要原因?；跇O限平衡理論分析邊坡穩(wěn)定的方法，沒有考慮邊坡內(nèi)部應(yīng)力～應(yīng)變的關(guān)系，計算得到的只是假定滑動面上平均安全系數(shù)，無法求出滑體內(nèi)部以及滑動面上的真實反力與內(nèi)力。有限元理論的不斷發(fā)展與完善以及計算機技術(shù)的日益成熟，逐漸成為工程研究和分析領(lǐng)域的重要方法。促使相應(yīng)的分析軟件也大量涌現(xiàn)，本文采用ABAQUS 軟件對邊坡穩(wěn)定性進行計算。將連續(xù)求解域離散為一組，使其按某種方式相互聯(lián)結(jié)形成有限個單元集合體，用近似函數(shù)和導數(shù)來表示每個單元體內(nèi)各結(jié)點數(shù)值和插值函數(shù)，結(jié)點未知量若被求出，利用相應(yīng)插值函數(shù)可求出單元內(nèi)場函數(shù)近似值。喻波等［1］用ANSYS 建立邊坡模型，結(jié)合ABAQUS 軟件中的Mohr－Coulomb 準則，通過邊坡穩(wěn)定性分析預(yù)測邊坡潛在滑動面，得出邊坡穩(wěn)定性系數(shù)。Cai、Ugai(2004)［2］通過滲流和變形非藕合方法，采用有限元強度折減法分析了降雨對邊坡穩(wěn)定性的影響規(guī)律。劉成［3］通過有限元強度折減法，使用ABAQUS 軟件，變動計算參數(shù)分析路堤邊坡的土性參數(shù)、邊界條件、相對滲透值及浸潤線分布對邊坡穩(wěn)定性的影響。

1 強度折減基本原理

邊坡在自重或外部荷載作用下發(fā)生破壞時，塑性區(qū)貫穿于整個邊坡形成滑裂帶，邊坡整體將沿滑裂帶滑移。導致滑裂帶上部的邊坡體形成機動結(jié)構(gòu)，而其有限元的計算結(jié)果不收斂，基于此提出了強度折減系數(shù)法。

強度折減法是通過不斷折減巖土體的剪切強度，使邊坡從穩(wěn)定到破壞的變化過程，從中找出邊坡巖土體的最薄弱部分。通過將邊坡土體的強度參數(shù)粘聚力c 以及內(nèi)摩擦角φ 的正切值tanφ 同除一個相同的折減系數(shù)F，可得到新的一組ce、φe值，將其作為新的計算參數(shù)輸入，再次試算收斂與否。不間斷的調(diào)整折減系數(shù)F 進行試算，直至相應(yīng)的折減系數(shù)F 使計算結(jié)果剛好收斂(F 若微小增加10－3，計算結(jié)果不收斂)。這時所對應(yīng)的Fs 稱為邊坡的最小穩(wěn)定安全系數(shù)，邊坡處在極限狀態(tài)，剛剛發(fā)生剪切破壞。

強度參數(shù)的折減按以下公式進行:

2 紅砂巖填料干濕循環(huán)試驗

采用ABAQUS 軟件進行路堤邊坡穩(wěn)定性分析需要填料的強度參數(shù)(粘聚力和內(nèi)摩擦角)，本文通過室內(nèi)模擬紅砂巖填料的干濕循環(huán)進行大直剪試驗，得到紅砂巖填料在干濕循環(huán)作用下的強度參數(shù)，為路堤邊坡的穩(wěn)定性分析提供數(shù)據(jù)。

2.1 試驗儀器

選用自主研發(fā)的“一種干濕循環(huán)大型應(yīng)變控制式直剪裝置”對紅砂巖填料進行干濕循環(huán)試驗研究。該儀器主要由垂直加載系統(tǒng)、水平加載系統(tǒng)、量測系統(tǒng)、反力框架、干濕循環(huán)模擬系統(tǒng)以及剪切盒等6 大部分組成。一共配備2 個剪切盒，其尺寸為0.5 m×0.5 m×0.3 m(長寬高)，儀器實物圖如圖1 所示。

圖1 試驗儀器實物圖

2.2 試驗用料

試驗用料取自江西某高速紅砂巖料場，最大粒徑60 mm(表1)，巖性主要是石英砂巖，部分夾有泥質(zhì)粉砂巖。由于試件制備困難，干濕循環(huán)模擬周期較長，未開展平行試驗。

表1 紅砂巖填料的粒組含量

2.3 試驗步驟

1)制備試樣。取足量紅砂巖填料按照施工所用最佳含水率拌勻并悶料，然后將其分層搗實填筑于剪切盒內(nèi)，每層高度保持在10 cm 上下?？刂谱钌蠈油翗拥谋砻婢嚯x上剪切盒面下3 ～5 cm，整平最上層土體表面，使得傳壓板與試樣表面接觸充分。將每層土體擊實至規(guī)定要求的層高，將表面刨毛后再填筑上一層，重復上述步驟直至最上層土體填筑完成。

2)干濕循環(huán)模擬。完成紅砂巖土樣制備之后，將水注入儲水環(huán)中，水慢慢滲入試樣中。使試樣浸泡在水中24 h，儲水環(huán)中一直保持足夠的水以保證浸泡過程的完整，模擬干濕循環(huán)中的濕潤狀態(tài)。完成浸泡后，將加熱裝置開啟，烘干試樣24 h，模擬干濕循環(huán)中的干燥狀態(tài)。每一次的“烘干～浸泡”的模擬一次干濕循環(huán)。重復“烘干～浸泡”次數(shù)，使其達到試驗所設(shè)定的干濕循環(huán)次數(shù)。

3)測量設(shè)備安裝，記錄垂直、水平方向的千斤頂和百分表的讀數(shù)。當水平方向荷載讀數(shù)不再增加或剪切變形急驟增長時，可認為試樣已剪損，清理剪切盒。若無上述兩種情況出現(xiàn)，應(yīng)控制剪切變形達試樣直徑的1/15 ～1/10，方可停止試驗。應(yīng)控制試樣在5 ～10 min 內(nèi)達到剪切破壞。

圖2 為剪切試驗圖。

圖2 剪切試驗圖

2.4 試驗結(jié)果分析

歷經(jīng)0、1、2、3、4 次干濕循環(huán)的紅砂巖填料干濕循環(huán)直剪試驗結(jié)果如表2 所示。

表2 紅砂巖填料直剪試驗結(jié)果

根據(jù)庫侖定理，擬合抗剪強度與垂直壓力的關(guān)系曲線，可以得到紅砂巖填料的抗剪強度參數(shù)如表3 所示。

由表3 中的強度參數(shù)，可直觀看出，隨干濕循環(huán)次數(shù)的增多干濕循環(huán)次數(shù)紅砂巖填料粘聚力和內(nèi)摩擦角都呈減小的態(tài)勢，第1 次和第2 次干濕循環(huán)作用下衰減強烈，第3 次和第4 次干濕循環(huán)作用時強度參數(shù)衰減作用減弱。

表3 紅砂巖填料強度參數(shù)

3 有限元模型建立

采用ABAQUS 分析路堤邊坡穩(wěn)定性首先得建立相關(guān)的模型，確定適宜模型尺寸、邊界條件，選定正確的本構(gòu)模型以及參數(shù)才能合理地對路堤邊坡進行穩(wěn)定性分析。

3.1 模型尺寸

模型邊界的大小直接影響到應(yīng)力應(yīng)變的分布以及計算的精確度，鄭穎人［4］對路基邊坡穩(wěn)定性有限元計算的研究表明，路堤邊坡坡角到左端邊界的距離為坡高的1.5 倍，且上下邊界總高不低于2 倍坡高時，計算精度最為理想。如圖3。

圖3 路基邊坡有限元模型邊界范圍

3.2 本構(gòu)模型

參考邊坡土體材料的實際特性，本文選定基于摩爾－庫倫(Mohr－Coulomb)準則的理想彈塑性的本構(gòu)模型。該模型假定:某處的剪應(yīng)力等于抗剪強度時，該處產(chǎn)生應(yīng)力破壞，該處作用的剪切強度和正應(yīng)力呈線性相關(guān)?；诓牧掀茐膽?yīng)力狀態(tài)的摩爾圓提出了摩爾－庫倫塑性模型，破壞線是與摩爾圓相切的直線，其屈服準則如下:

τ=c－σ

圖4 為ABAQUS 中的Mohr－Coulomb 破壞面。

當摩擦角在0≤φ ＜90°時，摩爾－庫倫模型退化為與圍壓無關(guān)的Tresca 模型，π 平面上的屈服面為正六角形，當φ=90°時，摩爾－庫倫模型將演化為與Rankine 模型，π 平面上的屈服面為正三角形。

圖4 ABAQUS 中的Mohr－Coulomb 破壞面

3.3 參數(shù)取值

為了正確的分析路基邊坡穩(wěn)定性的情況，要合理的確定參數(shù)的取值。采用摩爾－庫侖準則的強度折減方法分析邊坡的穩(wěn)定性時，影響邊坡安全系數(shù)的土體參數(shù)為土體的容重、彈性模量、泊松比、粘聚力和內(nèi)摩擦角。研究表明彈性模量與泊松比對邊坡穩(wěn)定安全系數(shù)的計算影響極小。如鄭穎人等［4］的計算表明泊松比對安全系數(shù)計算結(jié)果幾乎沒有影響，只影響塑性區(qū)的大小;然而彈性模量對邊坡的變形和位移大小有影響，但對安全系數(shù)的計算也沒有影響，因而在考慮邊坡的干濕循環(huán)次數(shù)中可取剛度系數(shù)為一個定值。類似極限平衡法，有限元強度折減法中影響邊坡安全系數(shù)的土體參數(shù)主要是強度指標粘聚力和內(nèi)摩擦角。地基參數(shù)可以通過地勘和設(shè)計文件獲得，而路堤填料的抗剪強度參數(shù)選自上述室內(nèi)紅砂巖填料的模擬干濕循環(huán)試驗。

4 有限元模擬工程實例

4.1 模型建立

選取江西某高速一高填路段位置斷面，填方高達14.2 m。紅砂巖填料路堤，地基分為粘土和頁巖兩層。幾何尺寸如圖5 所示。

圖5 紅砂巖填料路堤邊坡斷面圖(單位:m)

邊界約束:邊坡底面采用完全約束條件，左右采用水平約束條件，上部為自由邊界。土體采用摩爾－庫侖破壞準則與非關(guān)聯(lián)流動法則的理想彈塑性本構(gòu)模型。單元為CPE8R(8 節(jié)點平面應(yīng)變四邊形減縮積分單元)，平坦地基上紅砂巖填料路堤邊坡有限元計算中網(wǎng)格劃分如圖6 所示。

圖6 紅砂巖填料路堤邊坡網(wǎng)格劃分圖

4.2 參數(shù)選取

模型分紅砂巖填料路堤、地基淺層粘土和地基深層頁巖3 層，每層都當成均質(zhì)材料考慮。紅砂巖填料的粘聚力和內(nèi)摩擦角采用室內(nèi)試驗數(shù)據(jù)，彈性模量取42 MPa，泊松比取0.33，重度為19.6 kN/m3。查閱地勘資料，得到地基粘土的土性參數(shù)(不考慮地基粘土和頁巖的干濕循環(huán)效應(yīng))，匯總該處邊坡紅砂巖填料和地基粘土、頁巖的力學參數(shù)如表4 所示。

表4 邊坡斷面各層巖土力學參數(shù)

4.3 計算結(jié)果分析

當紅砂巖填料路堤不經(jīng)歷干濕循環(huán)時，不同的折減系數(shù)(Fs)條件下紅砂巖填料路堤邊坡的等效塑性應(yīng)變分布情況如圖7 所示。可以看到，塑性區(qū)首先在路堤邊坡的坡腳產(chǎn)生，隨著折減系數(shù)的增大，塑性區(qū)逐步發(fā)展擴大，最后直至趨于貫通。當折減系數(shù)為Fs=2.128 時，坡腳塑性區(qū)基本貫通坡面，計算也不再收斂。因此，判斷紅砂巖填料路堤邊坡的穩(wěn)定系數(shù)為2.128。

圖7 不同折減系數(shù)下紅砂巖填料路堤邊坡等效塑性應(yīng)變云圖

依照上述模型，取歷經(jīng)1、2、3、4 次干濕循環(huán)下紅砂巖填料強度參數(shù)計算出該邊坡的穩(wěn)定系數(shù)如表5 所示，繪出紅砂巖填料路堤邊坡穩(wěn)定系數(shù)隨干濕循環(huán)次數(shù)的變化曲線如圖8。

表5 干濕循環(huán)下紅砂巖填料路堤邊坡穩(wěn)定系數(shù)

圖8 干濕循環(huán)下紅砂巖填料路堤邊坡穩(wěn)定系數(shù)曲線

分析圖8、表5 可知，該處紅砂巖填料邊坡的穩(wěn)定性隨干濕循環(huán)次數(shù)的增多而減小，在干濕循環(huán)作用下邊坡的穩(wěn)定性系數(shù)由2.128 降低到1.672，降幅約為21%。穩(wěn)定系數(shù)前3 次干濕循環(huán)作用時下衰減較快，第4 次干濕循環(huán)作用下穩(wěn)定系數(shù)衰減較小，穩(wěn)定系數(shù)趨于穩(wěn)定。

5 結(jié)論

1)隨干濕循環(huán)次數(shù)的增多干濕循環(huán)次數(shù)紅砂巖填料粘聚力和內(nèi)摩擦角都呈減小的態(tài)勢，第1 次和第2 次干濕循環(huán)作用下衰減強烈，第3 次和第4次干濕循環(huán)作用時強度參數(shù)衰減作用減弱

2)紅砂巖填料邊坡的穩(wěn)定性隨干濕循環(huán)次數(shù)的增多而減小，在干濕循環(huán)作用下邊坡的穩(wěn)定性系數(shù)由2.128 降低到1.672，降幅約為21%。

［1］喻 波，黃政寧.基于ANSYS 和ABAQUS 強度折減邊坡穩(wěn)定性分析［J］.公路工程，2008，33(2):47－53.

［2］CAI F，UGAI K.Numerieal analysis of rainfall effects on slope stability［J］.International Journal of Geomeehanies，ASCE，2004，4(2):69－78.

［3］劉 成.路(河)堤邊坡穩(wěn)定性數(shù)值分析及圖表法研究［D］.大連:大連理工大學，2010.

［4］鄭穎人，趙尚毅.有限元強度折減法，在土坡與巖坡中的應(yīng)用［J］.巖石力學與工程學報，2004，23(19):3381－3388.

［5］崔 溦，張志耕，閆澍旺.膨脹土的干濕循環(huán)性狀及其在邊坡穩(wěn)定性分析中的應(yīng)用［J］.水利與建筑工程學報，2010，8(5).

［6］趙艷林，曾召田，呂海波.干濕循環(huán)對膨脹土變形指標的影響［J］.桂林工學院學報，2009，29(4).

［7］王建華，高玉琴.干濕循環(huán)過程導致水泥改良土強度衰減機理研究［J］.中國鐵道科學，2006，27(5):23－27.

［8］萬少石.涉水邊坡穩(wěn)定性的三維強度折減有限元分析［D］.大連:大連理工大學，2009.

［9］李 萍，付 強，陳小念.基于是強度折減法分析邊坡穩(wěn)定［J］.甘肅科學學報，2007，19(4):63－66.

［10］賈蒼琴，黃茂松，王貴和.非飽和非穩(wěn)定滲流作用下土坡穩(wěn)定分析的強度折減有限元方法［J］.巖石力學與工程學報，2007，26(6):1290－1296.

［11］陳建斌，郝紅玉，盧方偉.非飽和膨脹土強度特性試驗研究［J］.土工基礎(chǔ)，2008，22(6):76－80.