宋民崇，余云燕

(1.中國市政工程華北設計研究總院有限公司，天津 300074;2.蘭州交通大學 土木工程學院,甘肅 蘭州 730070)

?

三維數值模擬在某高速公路邊坡穩(wěn)定性分析中的應用

宋民崇1，余云燕2

(1.中國市政工程華北設計研究總院有限公司，天津 300074;2.蘭州交通大學 土木工程學院,甘肅 蘭州 730070)

摘要：近年來三維數值模擬在邊坡穩(wěn)定性分析中被廣泛應用，文中以某滑坡為工程背景，運用FLAC3D軟件分析滑坡的穩(wěn)定狀態(tài)，對滑坡體的位移、應力、應變、最大不平衡力和特定記錄點的位移與時步關系進行分析，模擬出滑動面的位置，綜合評價滑坡的穩(wěn)定狀態(tài)，從而為其它同類工程提供有意義的經驗。

關鍵詞：滑坡；穩(wěn)定性分析；FlAC3D；數值模擬

隨著國家經濟的發(fā)展，對基礎設施的投資日益加大，特別是西部大開發(fā)的開展，鐵路、公路建設量逐年增長，而且不斷向山區(qū)延伸，因地形、地貌等條件的約束，不可避免地遇到大量的邊坡工程，有些地區(qū)甚至出現高陡邊坡，不少公路、鐵路線段常有滑坡等地質災害發(fā)生，造成運輸中斷，大量人員傷亡和嚴重的經濟損失。所以進行邊坡穩(wěn)定性研究非常重要。

1工程概況

1.1地形地貌

坡體位于西南山區(qū)，地貌形態(tài)屬于侵蝕構造中山地貌，山體主要由板巖、千枚巖等組成。境內山巒重疊，山坡多在35°以上，懸崖峭壁隨處可見，山脈為石質山地，北側為下古生界碧口群石英巖夾變質砂巖組成的巖質岸坡，上部覆蓋第四系薄層坡積碎石土，多為巖石風化殘積物。

1.2地層巖性

滑坡區(qū)地層巖性主要為第四系坡積碎石、滑坡堆積塊石、碎石和土、下元古界碧口群，巖性主要有板巖、砂質板巖、變質砂巖、石英片巖。從新到老依次如下：

第四系坡積碎石層：分布于滑坡北側平緩山坡上?；尹S色，主要由碎石和土組成，碎石成分為變砂巖、板巖等，分級差，成棱角狀，結構松散，局部具有架空現象。

滑坡堆積塊石、碎石層：分布于南側，主要由塊、碎石和土組成，地表植被發(fā)育，富含腐殖質。前緣巖體破碎，多彎曲變形。

下元古界碧口群板巖：巖石的力學性能較好，具有抗風化能力強、巖體致密而堅硬的特點。變質砂巖夾石英片巖，新鮮面呈青灰色，微具片理，變余細粒結構，厚層狀結構，單層厚0.6～0.8 m，巖性致密堅硬。

1.3水文條件

所處區(qū)域降水量較大，對坡體的影響主要為兩種途徑：①基巖裂隙水?；鶐r裂隙水主要受大氣降水補給，沿層理、節(jié)理、裂隙下滲向河谷排泄。②孔隙潛水。賦存于第四系松散覆蓋層中，受大氣降水補給，排泄以地下徑流和蒸發(fā)為主。

地下水水質及水化學類型：經地表水取樣并進行水質分析，不含侵蝕性CO2、水化學類型基本為HCO3-Ca-Mg、SO4-Ga型水?；聟^(qū)地下水對混凝土無腐蝕性，對鋼筋混凝土結構中的鋼筋無腐蝕性。

滑坡體目前主要受自然條件的影響，包括地形地貌、地層巖性、水文地質條件等，沒有人類活動如施工擾動的影響。某滑坡工程地質斷面圖如圖1所示。

圖1　某滑坡工程地質斷面圖

2數值計算

本文采用FLAC3D對滑坡的穩(wěn)定安全系數進行計算，基于強度折減法的數值計算相對于傳統(tǒng)的極限平衡法所得到的安全系數有很多優(yōu)點：可以模擬出復雜的邊界條件，地質結構特征，各種作用荷載和巖土體的各種力學特性，并且可以考慮巖土體的應力與應變的關系?；贔LAC3D的強度折減法不需要預先假定滑面，在求解安全系數時就能夠自動搜索出滑移面的位置和形狀，可以對滑坡的漸近破壞過程進行一定的分析。

強度折減方法應用到有限差分或有限元分析中可以表述為:保持巖土體的重力加速度為常數，通過逐步減小抗剪強度指標，將c，φ值同時除以折減系數Fs，得到一組新的強度指標ci，φi，然后進行有限元分析，反復計算直至邊坡達到臨界破壞狀態(tài)，此時采用的強度指標與巖土體原具有的強度指標之比即為該邊坡的安全系數Fs。公式如下:

2.1計算模型的建立

根據地質資料建立三維計算模型，FLAC3D建模工程浩大、繁瑣，本文采用ANSYS建立模型，然后將模型導入FLAC3D進行數值計算，邊坡前緣的腳點為坐標原點(0,0,0) ，X方向沿著坡腳方向向右為正，Y方向以滑動向前為負，Z向上為正，向下為負，高度為150 m，邊坡長150 m，寬150 m，初始地應力為自重地應力場。邊界條件為：X,Y,Z方向為固定約束，滑坡體主要由碎石土組成，滑床主要由板巖構成。計算參數選取根據勘測資料，按照表1取值，計算工況為天然狀態(tài)。

表1　巖土體力學參數

在計算巖土體的剪切模量和體積模量時，應根據以下公式計算：

式中：E為巖土體彈性模量，B為巖土體體積模量，S為巖土體剪切模量，μ為材料泊松比。

2.2數值計算結果

應用FLAC3D軟件進行數值模擬，計算得到滑坡安全系數為0.97，如圖2所示，并對位移、應力、應變、最大不平衡力和特定記錄點的位移與時步關系進行分析。圖3為ANSYS網格劃分示意圖，圖4為滑坡體示意圖。

圖2　安全系數計算結果

圖3　ANSYS網格劃分示意圖

圖4　滑坡體示意圖

由于滑坡體主要成分為碎石土，力學性質較低，滑床為板巖，較為堅固，滑坡內部應力場由于外界條件的影響產生變化，本文應用TECPLOT圖形處理軟件出圖，下面對計算結果進行分析。

圖5為滑坡體總位移云圖，圖6為滑坡體總體位移中線位置剖面云圖，在總體位移云圖中可以看到，滑坡體最大位移為15 cm，滑坡體位移從后緣中心位置向外逐漸減小，最大位移產生在滑坡體的后緣中心位置處，從中線位置剖面云圖中可以看出最大位移是從滑坡體后緣中心位置向下逐漸減小。

圖5滑坡體總位移云圖

圖6　滑坡體總位移中線位置剖面云圖

圖7為滑坡體X方向位移云圖，圖8為滑坡體X方向中線位置位移剖面云圖。由圖7、圖8可以看出，滑坡體在X方向的最大位移為3 cm，在X方向上產生位移很小，產生的位移是由于滑坡體一側略微高于另一側，在總位移云圖中可以看出藍色區(qū)域位移為0，其它區(qū)域位移為毫米級。

圖7　滑坡體X方向位移云圖

圖8　滑坡體X方向中線位置位移剖面云圖

圖9為滑坡體Y方向位移云圖，圖10為滑坡體Y方向中線位置位移剖面云圖。由圖9、圖10可以看出，滑坡體在Y方向產生位移，即為滑坡滑動方向，最大位移為11 cm，在Y方向產生位移主要是滑坡向下滑動導致的，在中線剖面云圖中也看出位移是由上到下逐漸減小。

圖9　滑坡體Y方向位移云圖

圖10　滑坡體Y方向中線位置位移剖面云圖

圖11為滑坡體Z方向位移云圖，圖12為滑坡體Z方向中線位置位移剖面云圖，從圖11可以看出，滑坡體最大位移為11 cm，即在上下方向上產生豎向位移，從圖12可以看出，Z方向位移也是由滑坡體后緣中心位置由上到下逐漸減小的。

圖11　滑坡體Z方向位移云圖

圖12　滑坡體Z方向中線位置位移剖面云圖

圖13為主剖面最大主應力云圖，圖14為主剖面最小主應力云圖，由圖13、圖14可以看出，滑坡體最大主應力在分布上基本平行于斜坡臨空面，最大主應力約為0.2～2.8 MPa，最小主應力約為0.05～0.55 MPa，都為壓應力。

圖13　主剖面最大主應力云圖

圖14　主剖面最小主應力云圖

圖15為滑坡體X方向應力云圖，圖16為滑坡體Y方向應力云圖，圖17為滑坡體Z方向應力云圖，由圖15～17可以看出，X方向應力較小，這與X方向上產生位移較小相符合，Y方向的應力最大值出現在坡腳，Z方向應力較大，即豎直方向上產生的應力值較大，并且隨著深度的增加而增加，這符合地應力賦存規(guī)律。

圖15　滑坡體X方向應力云圖

圖16　滑坡體Y方向應力云圖

圖17　滑坡體Z方向應力云圖

圖18、圖19分別為滑坡體剪應變增量云圖和剪應變增變率云圖，由圖18、圖19可以看出，最大值出現在滑坡體的后緣，因此，滑坡體的破壞應該是后緣先開裂。在剪應變增量云圖中同時可以看出滑坡體滑動面的位置。

圖18　滑坡體剪應變增量云圖

圖19　滑坡體剪應變增變率云圖

圖20　滑坡體內點(75,95,120)X方向位移變化歷史圖

在FLAC3D計算時，可在滑坡體內設置一點(75,95,120)，記錄這一點隨時步的變化，同時也可以記錄滑坡體最大不平衡力的變化歷史。圖20～22分別為滑坡體內某一點X方向、Y方向和Z方向的位移變化歷史，圖23為最大不平衡力變化歷史。由圖20～22可以看出，X方向產生微小的位移，到計算結束時，記錄點有2.2 cm的位移，隨著時步的增加，位移是從0逐漸增加的，記錄點Y方向產生8.1 cm的位移，即指向臨空面方向，也是滑坡滑動方向，記錄點在Z方向上產生7 cm的位移，即向下有7 cm的位移，此記錄點在滑坡體內中心位置。對記錄點的選取有一定的代表性，記錄點要能反映滑坡體的穩(wěn)定狀態(tài)，通過位移的變化來說明滑坡體所處狀態(tài)。由圖23可以看出，最大不平衡力逐漸趨于0，計算是收斂的。

圖21　滑坡體內點(75,95,120)Y方向位移變化歷史

圖22　滑坡體內點(75,95,120)Z方向位移變化歷史

圖23　最大不平衡力變化歷史

圖24、圖25分別為滑坡體位移矢量圖和速度矢量圖。由圖24、圖25可以看出，這兩個矢量云圖的矢量方向基本是重合的，顯示出了滑坡體的運動趨勢，即向著主滑方向滑動，主要是Z方向和Y方向，這兩張云圖同時也顯示了滑坡體的范圍。

圖24　滑坡體位移矢量圖

圖25　滑坡體速度矢量圖

3結論

本文主要是對某高速公路滑坡進行穩(wěn)定性分析，應用FLAC3D軟件對其進行了定量的計算，所做的主要工作和得到的相關結論可以概括為以下幾點：

1) 根據地質資料運用ANSYS建立三維計算模型，然后將其導入FLAC3D軟件進行計算，計算得到滑坡體安全系數為0.97，并對位移、應力、應變、最大不平衡力和特定記錄點的位移與時步關系進行定量的分析。

2)通過三維數值模擬，對滑坡體運動趨勢進行分析，在其主滑方向即Y方向最大位移是11 cm，豎直向下即Z方向的位移為11 cm，X方向上位移較小，最大位移僅為3 cm。

3)對滑坡體應力的分析得出Y方向的應力最大值出現在坡腳，Z方向的應力值較大，并且隨著深度的增加而增加，符合地應力賦存規(guī)律，對剪應變增量的分析可以得出滑坡的破壞將是后緣先開裂，同時也可以判定出滑坡體滑動面的位置。

滑坡體位于高速公路沿線，目前處于不穩(wěn)定狀態(tài)，所處區(qū)域多雨且地震活動頻繁，為了保證日后公路運營的安全，建議采取邊坡防護、治理措施。

參考文獻

[1]徐邦棟.滑坡分析與防治[M].北京:中國鐵道出版社,2001.

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

[3]劉波,韓彥輝.FLAC原理、實例與應用指南[M].北京:人民交通出版社,2005.

[4]王恭先,徐峻齡,劉光代,等.滑坡學與滑坡防治技術[M].北京:中國鐵道出版社，2007.

[5]宋民崇，余云燕，趙德安.FLAC軟件在邊坡穩(wěn)定性分析中的應用[J].黑龍江工程學院學報(自然科學版)，2009，23(4)：24-27.

[6]陳育民,徐鼎平.FLAC/FLAC3D基礎與工程實例[M].北京:中國水利水電出版社,2009.

[7]許強,黃潤秋,程謙恭,等.三峽庫區(qū)泄灘滑坡滑帶土特征研究[J].工程地質學報，2003，11(4)：354-359.

[8]張永安.滇中紅層路塹邊坡與路堤穩(wěn)定性研究[D].昆明：昆明理工大學，2010.

[9]莫素華.邊坡穩(wěn)定性研究方法[J].中國水運(學術版)，2007，7(2):100-102.

[10] 俸錦福;郝勇;顏丙山.邊坡穩(wěn)定性分析理論及方法綜述[J].中國水運(下半月),2009,9(7):183-184.

[責任編輯：劉文霞]

Application of three-dimensional numerical simulation tothe stability analysis of a highway slope

SONG Minchong1，YU Yunyan2

(1.North China Municipal Engineering Design and Research Institute Co.,Ltd，Tianjin 300074，China；2.College of Civil Engineering，Lanzhou Jiaotong University，Lanzhou 730070，China)

Abstract：In recent years, three-dimensional numerical simulation has been widely used in landslide stability analysis.Taking a landslids as the engineering background,this paper makes use of FLAC3Dsoftware to analyze the stability of this landslide, such as landslide displacement, stress, strain, the maximum unbalanced force and the relationship between the displacement of a specific record point and the time step, and simulates the position of the sliding surface, and the comprehensive evaluation of the stability of the landslide,in order to provide a meaningful experience for othersimilar projects.

Key words：landslide; stability analysis; FlAC3D; numerical simulation

中圖分類號：TU457

文獻標識碼：A

文章編號：1671-4679(2016)01-0004-06

作者簡介：宋民崇(1983-)，男，工程師，研究方向：道路與橋梁.

收稿日期：2015-09-23