王偉興, 何文勇, 高　楊

(1.貴州大學(xué) 資源與環(huán)境工程學(xué)院， 貴州 貴陽　550025;　2.貴州省交通規(guī)劃勘察設(shè)計研究院股份有限公司, 貴州 貴陽　550081)

?

基于Geo-Studio和Flac 2D的邊坡流固耦合作用分析

王偉興1,2, 何文勇2, 高楊1,2

(1.貴州大學(xué) 資源與環(huán)境工程學(xué)院， 貴州 貴陽550025;2.貴州省交通規(guī)劃勘察設(shè)計研究院股份有限公司, 貴州 貴陽550081)

[摘要]以貴州省凱里至羊甲高速公路小苗嶺巖堆邊坡RK9+130～RK9+410右側(cè)滑坡項目為依托，通過邊坡應(yīng)力場與滲流場的耦合作用分析該邊坡因地下水位抬升發(fā)生失穩(wěn)的原因。采用Geo-Studio軟件中的滲流分析模塊Seep/W和邊坡穩(wěn)定性分析模塊Slope/W分析不同水力梯度下邊坡的安全系數(shù)，采用Flac 2D分析邊坡隨著地下水位的抬升塑性區(qū)的變化。分析和論證了由于降雨入滲，使坡體內(nèi)部產(chǎn)生滲流場，在流固耦合作用下使邊坡發(fā)生失穩(wěn)甚至破壞。

[關(guān)鍵詞]流固耦合； 邊坡穩(wěn)定性； Geo-Studio； Flac 2D； 塑性區(qū)

0前言

工程巖土體中存在應(yīng)力場、滲流場等多場并存且相互影響、相互作用，構(gòu)成一種耦合關(guān)系。滲流場產(chǎn)生的滲流體積力[1]作用于巖土介質(zhì)，會使巖土介質(zhì)應(yīng)力場和位移場發(fā)生變化；巖土介質(zhì)應(yīng)力狀態(tài)的變化將改變其滲流性質(zhì)，即應(yīng)力場、位移場的改變使得巖土介質(zhì)的孔隙比、孔隙率發(fā)生變化；同時由于多孔介質(zhì)的滲透系數(shù)[2]與其孔隙的分布情況關(guān)系很大，孔隙比、孔隙率的變化必然引起介質(zhì)滲透性能即滲透系數(shù)的改變，巖土介質(zhì)的滲流場受到影響也會發(fā)生改變。所以說巖土滲流性質(zhì)的變化又將改變滲流場的分布規(guī)律，同時也改變了巖土介質(zhì)的滲透力；應(yīng)力場對滲流場影響的實(shí)質(zhì)是應(yīng)力場改變了巖土介質(zhì)中孔隙的分布狀況，從而改變了巖土介質(zhì)的滲透特性。當(dāng)滲流場和應(yīng)力場發(fā)生改變時，常常伴有一系列的工程問題產(chǎn)生。本文以小苗嶺邊坡為研究對象，由于邊坡開挖過程中出現(xiàn)持續(xù)降雨導(dǎo)致地下水位抬升，使得坡體內(nèi)部的滲流場與應(yīng)力場發(fā)生改變且相互作用即流固耦合作用，通過Geo-Studio和Flac 2D分析流固耦合作用對邊坡穩(wěn)定性的影響。

1工程概況

1.1施工設(shè)計

貴州省凱里至羊甲高速公路路基挖方路段邊坡RK9+130～RK9+410，設(shè)計為六級邊坡，第一級坡比1∶1.25，第二～六級坡比1∶1.5。該邊坡于2013年9月開始開挖，2013年12月底開挖至第一級平臺后，坡體后側(cè)出現(xiàn)零星裂縫，2014年3月初，于填方段坡體后側(cè)發(fā)現(xiàn)開裂，且逐漸發(fā)展現(xiàn)已形成較大錯臺。2014年3月7日～3月28日期間對滑坡進(jìn)行了深層位移監(jiān)測，監(jiān)測結(jié)果顯示，滑坡最大變形深度達(dá)35.5 m，變形速率為9.0～13.0 mm/d，處于孺滑狀態(tài)且該期間場區(qū)范圍內(nèi)出現(xiàn)持續(xù)降雨。目前裂縫已基本貫通，統(tǒng)一滑面基本形成，最大裂縫已擴(kuò)展至路基右側(cè)260 m處，裂縫最大錯距1.7 m。

1.2地層巖性

路段區(qū)覆蓋層為殘、坡積層粉質(zhì)粘(Qel+dl)，坡積層角礫土(Qdl)及崩積層(Qc)碎石土。下伏基巖為二疊系上統(tǒng)吳家坪組(P2w)砂巖夾炭質(zhì)泥巖及志留系中統(tǒng)翁項組(S2w)砂巖。

1.3水文地質(zhì)

① 地下水。場區(qū)地下水為基巖裂隙水及第四系松散層類孔隙水，地下水主要靠大氣降水垂直補(bǔ)給。

② 水位的監(jiān)測。在RK9+260斷面(見圖1)共布置深部位移監(jiān)測孔3個即JCK03、JCK04、JCK05(見圖2)。

圖1　RK9+260斷面圖Figure 1　RK9+260 section

圖2　JCK03、JCK04、JCK05測斜管內(nèi)水位變化圖Figure 2　JCK03, JCK05, JCK04 water level variation in the inclined tube

由圖2可知： 在一個月內(nèi)的水位監(jiān)測中，JCK03、JCK04、JCK05內(nèi)水位均有不同程度抬升，尤其JCK04、JCK05內(nèi)水位抬升幅度較大約20 m，JCK03內(nèi)水位基本穩(wěn)定，2014-03-18到03-30期間邊坡區(qū)域范圍內(nèi)出現(xiàn)連續(xù)降雨，從圖中可知從2014-03-18開始水位變化幅度增大，在2014-03-30前后監(jiān)測孔相繼被剪斷，此時邊坡處于孺滑階段[3]，可以認(rèn)為2014-03-30所測得的水位為該監(jiān)測期間的最高水位，而自2014-03-02開始水位變化幅度較小且水位較為穩(wěn)定，此時的水位為該監(jiān)測期間的最低水位。由于滑動面的范圍遠(yuǎn)超過測斜管JCK05的位置，為此在計算當(dāng)中應(yīng)用的浸潤線位置應(yīng)該延伸到邊坡后緣，根據(jù)所測得的水位資料繪出5條不同水位情況見圖3。

圖3　不同水位條件下的浸潤線位置Figure 3　The position of the line at different water level

由圖3中可知: 邊坡后緣最高水位高程為930 m，最低水位高程為910 m，變化幅度為20 m，與監(jiān)測的水位資料情況相符，邊坡前緣水位基本穩(wěn)定水位高程為855 m，因此需要分析水位高程分別為910、915、920、925、930 m(即后緣水頭高度為70、75、80、85、90 m，前緣水頭高度為15 m)時坡體內(nèi)部流固耦合作用對邊坡穩(wěn)定性的影響。

1.4計算參數(shù)

碎石土的彈性模量為30 MPa，泊松比0.25，干重度20 kN/m3，飽和重度22 kN/m3，內(nèi)摩擦角為25°，粘聚力為10 kPa，滲透系數(shù)3.3e-4 m/sec。

2Geo-Studio的流固耦合作用分析

2.1Geo-Studio軟件耦合計算原理

在GeoStudio系統(tǒng)軟件中Seep/W模塊可以完成孔隙水壓力的計算并可以得到不同水力梯度下的浸潤線，然后將計算所得的孔隙水壓力的結(jié)果導(dǎo)入Slope/W模塊中得到在不同水力梯度下邊坡的安全系數(shù)即流固耦合作用[4，5]下的安全系數(shù)。

2.2Geo-Studio模型建立

由于通過Seep/W模塊和Slope/W模塊進(jìn)行耦合計算時，為了能更加真實(shí)的反應(yīng)實(shí)際情況，因此滑動面設(shè)定為“指定滑動面位置”的方式進(jìn)行耦合計算，根據(jù)現(xiàn)場深部位移監(jiān)測曲線(見圖1)據(jù)測斜管被剪斷掉的位置來確定滑動面的位置見圖1。因此以斷面RK9+260為依托建立計算模型見圖4。

圖4　Geo-Studio計算模型Figure 4　Geo-Studio calculation model

2.3基于Geo-Studio的不同水力梯度下的邊坡安全系數(shù)(見圖5)

不同水力梯度對邊坡穩(wěn)定性的影響(見表1)。

通過表1可知: 隨著水位的不斷升高，邊坡安全系數(shù)逐漸下降，當(dāng)水位上升至90 m即監(jiān)測期間的最高水位時，邊坡安全系數(shù)約為F=1.00，可以認(rèn)為邊坡處于極限平衡狀態(tài)[6，7]，通過現(xiàn)場測斜資料可知JCK03、JCK04、JCK05在此水位條件下相繼被剪斷，即邊坡開始出現(xiàn)滑動，與計算結(jié)果基本相符，如果水位繼續(xù)上升，邊坡將會失穩(wěn)而演變成滑坡。

圖5　不同水位時的邊坡穩(wěn)定性分析Figure 5　Analysis of slope stability under different water levels

表1 不同水位下安全系數(shù)Table1 Safetyfactorunderdifferentwaterlevels水頭/m方法OrdinaryBishopJanbuM-P701.2091.2611.1891.207751.1751.2201.1561.173801.1361.1741.1191.135851.0891.1211.0731.088901.0361.0641.0211.035

3Flac 2D的流固耦合作用分析

3.1Flac 2D模型建立

以RK9+260斷面(見圖1)為依托建立模型如圖6所示。

圖6　Flac 2D計算模型Figure 6　Flac 2D computing model

3.2基于Flac 2D的不同水力梯度下邊坡的塑性區(qū)變化

分析當(dāng)后緣水頭高度分別為70、75、80、85、90 m，前緣水頭高度為15 m時，邊坡塑性區(qū)的變化(見圖7)。

圖7　不同水位時的塑性區(qū)Figure 7　the plastic zone under different water level

結(jié)果分析：

通過不同水力梯度下的塑性區(qū)對比分析可以看出：由于邊坡前后緣存在一定的水力梯度[8，9]，在此水力梯度的作用下坡體內(nèi)部產(chǎn)生滲流場，隨著水位的抬升，在坡腳處產(chǎn)生塑性區(qū)，并且該塑性區(qū)有逐漸增大的趨勢，當(dāng)水位上升至90 m時，塑性區(qū)貫通，從圖15中可以看出: 坡體表面出現(xiàn)拉應(yīng)力區(qū)，若水位繼續(xù)上升，該邊坡將誘發(fā)成滑坡，從圖7可以看出: 坡體后緣出現(xiàn)拉應(yīng)力區(qū)，而通過現(xiàn)場可知邊坡后緣存在拉裂縫[10，11]，這與實(shí)際情況相符。

4監(jiān)測對比分析

在該邊坡坡面設(shè)置地表位移監(jiān)測點(diǎn)JCD18見圖1所示，邊坡從2013年9月份開始開挖到2014年8月份JCD18的位移曲線如圖8所示。

圖8　JCD18(ΔX,ΔY,ΔZ)位移曲線圖Figujre 8　JCD18 (ΔX,ΔY,ΔZ) displacement curve

從圖8中可以看出: 自2013年12月起邊坡開始發(fā)生變形，此時邊坡開挖至一級平臺，坡體后側(cè)出現(xiàn)零星裂縫，邊坡繼續(xù)開挖且場區(qū)出現(xiàn)持續(xù)降雨，邊坡變形逐漸增大，通過Geo-Studio和Flac 2D分析中可以看出邊坡地下水位逐漸抬升，邊坡安全系數(shù)[12，13]逐漸下降，在坡腳處邊坡塑性區(qū)[14，15]逐漸增大和邊坡位移量逐漸增大，這與圖8所示的位移變化情況基本相符。

5結(jié)論

① 利用GeoStudio軟件中Seep/W模塊和Slope/W模塊分析邊坡內(nèi)部流固耦合作用是一種基于剛體極限平衡理論的分析方法，得到不同水力梯度下的安全系數(shù)是流固耦合作用對邊坡穩(wěn)定性影響的一個定量的評價。

② 利用Flac 2D分析邊坡流固耦合作用，是基

于彈塑性本構(gòu)模型的，對于流固耦合作用的影響是一個定性評價，主要揭示該邊坡演變成滑坡的過程。

③ 通過2種不同的分析方法可知，隨著地下水位的抬升，邊坡安全系數(shù)逐漸下降直至出現(xiàn)孺滑狀態(tài)，同樣邊坡塑性區(qū)逐漸從坡腳向后緣延伸，最終達(dá)到塑性區(qū)貫通，即邊坡發(fā)生失穩(wěn)乃至破壞。

④ 由于地下水位的抬升使得坡體內(nèi)部出現(xiàn)水力梯度差，產(chǎn)生滲流場并與坡體內(nèi)部的應(yīng)力場相互作用，降低邊坡的安全儲備，對于工程造成損失甚至產(chǎn)生危險，因此在工程施工前應(yīng)注重排水的重要性。

⑤ 通過現(xiàn)場地表位移監(jiān)測資料的對比，進(jìn)一步驗(yàn)證了基于GeoStudio和Flac 2D的流固耦合作用分析的正確性和與現(xiàn)場情況的一致性，同時也證明了滲流場存在的真實(shí)性。

[參考文獻(xiàn)]

[1]仵彥卿,張悼元.巖體水力學(xué)導(dǎo)論[M].成都:西南交通大學(xué)出版社,1995.

[2]《工程地質(zhì)手冊》編寫委員會.工程地質(zhì)手冊[M].北京:中國建筑工業(yè)出版社,1992.

[3]JTGD30-2004,公路路基設(shè)計規(guī)范[S].

[4]解宏偉,揚(yáng)龍,蘇曉波.二維滲流場基本物理量的有限元計算法[J].青海大學(xué)學(xué)報:自然科學(xué)版,1999,17(4):21-23.

[5]姜立春,吳愛祥,陳嘉生.高陡邊坡地下水滲流場三維有限元分析及其實(shí)例研究[J].湘潭礦業(yè)學(xué)院學(xué)報,2003,18(3):24-28.

[6]李艷春．土質(zhì)學(xué)與土力學(xué)[M]．北京：中國建材工業(yè)出版社，2005．

[7]陳祖煜. 土質(zhì)邊坡穩(wěn)定分析：原理-方法-程序[M].北京：中國水利水電出版社,2003.

[8]柳群義，朱自強(qiáng)，何現(xiàn)啟，等.水位漲落對庫岸滑坡孔隙水壓力影響的非飽和滲流分析[J].巖土力學(xué),2008,29(S1):85-89.

[9]謝定義,姚仰平,黨發(fā)寧.高等土力學(xué)[M].北京: 高等教育出版社,2008:223-224.

[10]祁生文,伍法權(quán),蘭恒星,等.錦屏一級水電站普斯羅溝左岸深部裂縫成因的工程地質(zhì)分析[J].巖土工程學(xué)報，2002,24(5):596-599.

[11]張倬元,王士天,王蘭生,等.工程地質(zhì)分析原理[M].北京:地質(zhì)出版社,1994.

[12]GB50330-2002,建筑邊坡工程技術(shù)規(guī)范[S].

[13]趙尚毅,鄭穎人,時衛(wèi)民,等. 用有限元強(qiáng)度折減法求邊坡穩(wěn)定安全系數(shù)[J]. 巖土工程學(xué)報，2002，24(3)：343-346.

[14]袁文伯,陳進(jìn).軟化巖層中巷道的塑性區(qū)與破碎區(qū)[J].煤炭學(xué)報,1986,11(3):77-85.

[15]吳家龍.彈性力學(xué)[M].上海: 同濟(jì)大學(xué)出版社，1987.

Analysis of Fluid Solid Coupling of Slope Based on Geo-Studio and Flac 2D

WANG Weixin1,2, HE Wenyong2, GAO Yang1,2

(1.Resource & Environment Engineering College of GuiZhou university, Guiyang， Guizhou 550025，China;2.Guizhou Transportation Planning Survey & Design Academe CO.,LTD Guiyang， Guizhou 550081,China)

[Abstract]depend on the project that a highway of Xiao Miaoling from Kaili to Yangjia rock pile slope RK9+130 to RK9+410 right landslide in Guizhou province, the analysis of the reason that the slope is unstability due to underground water level uplift by the stress field and seepage field. calculate the safety factor of slope by Slope/W and Seep/W of Geo-Studio software, Flac 2D was used to analyze the change of the plastic zone of the slope. Analysis and demonstration of the rainfall infiltration, so that the field was produced in the slope, under the flow and solid coupling to make the slope instability and even damage.

[Key words]fluid structure coupling; slope stability; geo-studio; flac 2D; plastic zone

[中圖分類號]U 416.1+4

[文獻(xiàn)標(biāo)識碼]A

[文章編號]1674—0610(2016)02—0178—04

[作者簡介]王偉興(1990—)，男，吉林長春人，碩士研究生，從事邊坡穩(wěn)定性分析和治理方面研究，研究方向：巖土體工程性質(zhì)及其應(yīng)用。

[基金項目]貴州省交通運(yùn)輸廳科技項目(2012-022-132)

[收稿日期]2015—09—22