基于Flac3D的滑坡安全性分析

杜軒

摘要： 在我國西南部地區，滑坡已經成為高頻率發生、難以預測、危害程度大的一種地質災害。如何有效預防和及時發現滑坡已經成為道路和防災部門的首要任務。本文以昆民吳家營某處邊坡為例，應用Flac3D技術對邊坡的剪切變形趨勢和安全系數進行模擬和計算，對其穩定性和預防措施提出一些建議。

Abstract： In southwest China， the landslide has become a kind of geological disaster， which is a kind of high frequency occurrence， difficult to predict， and a great harm. How to effectively prevent and timely discovery of the landslide has become the primary task of the Road and Disaster Prevention Department. In this paper， taking a slope in Wujiaying in Kunming as an example， the Flac3D technology is used to simulate and calculate the shear deformation trends and the safe coefficient of the slope， and some suggestions on its stability and the preventive measures are put forward.

關鍵詞： 滑坡穩定性；Flac3D；模型；防治

Key words： slope stability；flac3D；models；prevention

中圖分類號：P642.22 文獻標識碼：A 文章編號：1006-4311（2016）07-0111-03

0 引言

我國尤其是西南山區存在非常多的自然邊坡，每年都有大量滑坡災害發生。影響邊坡穩定性的內部因素有邊坡巖土體類型、邊坡形態、邊坡地址構造和地形地貌、地下水等；外部因素有地震、風化作用、植被、氣候條件和人工活動等。隨著城市擴張和各地道路的興修，一些人工活動變得頻繁，如切坡、削坡、坡頂加載及其他破壞其邊坡兩面性的活動，最后導致邊坡失穩，形成滑坡。去年深圳12.26滑坡事故就是由于非法人為坡頂加載造成的。表1是由中國地質環境信息網所給出的《全國地質災害通報》，可了解我國近幾年地質災害情況：

從表1可知滑坡災害發生數量比其他災害之和還要多，因此如何有效預防和判斷滑坡成為當務之急。Flac3D軟件是一款以有限差分法分析連續介質的工程軟件，處理能力快捷且強大，廣泛應用于邊坡穩定性分析、隧道、礦山、基坑開挖等工程領域。本文中，應用Flac3D邊坡安全系數的強度折減法，求解其剪切變形趨勢，判斷到其穩定性，從而采取對應的工程措施進行預防和控制。

1 工程概況

邊坡（見圖1）位于云南昆明呈貢區吳家營村。邊坡坡體為一級臺階，發育完整，分布連續，為典型的山地丘陵地貌。長約120-150m，坡高10-15m，坡角為30°-45°左右，坡面向陽，表面有植被覆蓋，多為低矮灌木。坡體土質為粉質粘土和砂土組成，坡角和坡頂無人工擾動痕跡。由于其坡體坡度較大，土質為軟弱土體，因而選取此邊坡分析更具有典型性意義。

2 Flac3D建模與求解

2.1 Flac3D簡介

Flac3D軟件是美國Itasca咨詢有限責任公司開發研制的一種巖土工程專業分析軟件。由于其強大而快速的計算功能和廣泛的模擬功能，尤其是在各種變形問題上具有獨特的分析能力優勢，使的Flac3D在國際巖土工程界非常流行。其一般的建模求解過程如圖2所示。

2.2 初始模型建立

圖3所示為模擬的邊坡體，選取部分坡體，基體長20m，寬10m，高13m，滑面高10m，長10m，邊坡角度為45°（模擬最大坡角）。假定邊坡為連續均勻介質的巖土體，經室內試驗測得其平均密度ρ為1.97×103kg/m3，變形模量E為80MPa，泊松比ν為0.37，體積模量K=E/3（1-2ν）=103MPa，剪切模量G=E/2（1+ν）=29.2MPa，抗拉強度為12.38kPa，內摩擦角？準=20°剪脹角Ψ=18°。模型每隔一米設置一等分，共計5513個單元格。其中X、Y、Z軸方向如圖3所示。

2.3 命令流及求解

強度折減法中邊坡穩定的安全系數定義為：使邊坡剛好達到臨界破壞狀態時，對巖、土 體的抗剪強度進行折減的程度，即定義安全系數為巖土體的實際抗剪強度與臨界破壞時的折減后剪切強度的比值。通過不斷地增加折減系數，反復計算，直至邊坡達到臨界破壞，此時得到的折減系數即為安全系數 F。在不考慮地表徑流和地下水的情況下，根據邊坡的各種物理力學參數編寫其命令如下：

gen zone brick p0 0 0 0 p1 2 0 0 p2 0 10 0 p3 0 0 3 size 3 1 3

gen zone brick p0 2 0 0 p1 20 0 0 p2 2 10 0 p3 2 0 3 size 17 1 3 ratio 1.03 1 1

gen zone brick p0 2 0 3 p1 20 0 3 p2 2 10 3 p3 12 0 13 p4 20 10 3 p5 12& 10 13

P6 20 0 13 P7 20 10 13 SIZE 17 1 17 ratio 1.03 1 1 // 建立如圖3所示模型體

fix x y z range z -0.1 0.1 // 設置邊界條件，即模型體只能在X和Z軸

fix x range x 19.9 20.1 方向上發生形變，在Y軸無形變

fix x range x -0.1 0.1

fix y

model elas //設置彈性本構模型

prop density 2000 bulk 3e9 shear 1e9 //設置力學參數和重力加速度

set gravity 0 0 -10.0

solve

ini xdisp 0 ydisp 0 zdisp 0 //位移場和速度場清零

ini xvel 0 yvel 0 zvel 0

model mohr //設置庫倫-摩爾模型

prop density 1970 bulk 1.03e8 shear 2.92e7 & //設置物理力學參數

coh 12380.0 tens 1.0e6 fric 20 dila 18

solve fos file slope3dfos1.sav associated //求解安全系數

plot fos cont ssi outline on vel red //繪制安全系數、剪切應變增量圖

得到的模型安全系數、剪切應變增量圖見圖4。

從圖中可以明顯看到塑性貫通區域即潛在滑動面，因滑動面外側區域各網格點速度明顯大于其它區域，尤其在坡腳處更加明顯。而且其邊坡安全系數（FoS value）為1.04，稍稍大于臨界破壞值1，說明這一區域已經出現明顯滑動并即將發生破壞。

3 治理工作

針對上述滑坡存在安全隱患和潛在破壞性，可以考慮以下幾個方面的治理：

①傳統的治理方式如在坡腳設置擋土墻，坡體打錨桿、錨索、鋼絲繩網和抗滑樁等，這些方案是在滑坡治理中采用最多、見效最快的。一般為道路邊坡防護所采用（圖5為錨索擋土墻）。

②高度重視排水工作，坡體邊界外設置截水溝，在滑坡坡體上設置樹枝狀或網狀排水系統，將坡面徑流匯集、旁引于滑坡體外排除，采用各種形式的盲溝、仰斜排水孔等措施，將滑坡體內地下水導出，增強坡體的自然穩定性。

③對于滑坡坡度較大者，可以在坡頂清方卸載回填于坡腳，以減少其下滑力。在坡體種植根系發達的喬木、灌木，增強坡體的土體粘聚力。采用固結灌漿或電化學加固法加強邊坡巖土體強度。

④針對其有潛在的滑坡風險，可以設置滑坡位移監測儀，并定期觀察坡體地面是否有開裂、脫落或鼓脹等情況，及時掌握其動態，做好預防工作。

4 結束語

針對滑坡災害的頻發性和危害程度，如何有效預防和應對滑坡的發生已經成了重中之重。Flac3D軟件通過有效的模型建立和求解，對滑坡的剪切變形趨勢和安全系數進行了具體化的圖形數據顯示，有效的體現了滑坡的危險性和潛在性，使人們能夠有效的預防和控制滑坡災害的發生。隨著國家工程建設持續發展，現代化城市不斷地擴張，解決邊坡失穩問題及滑坡防治技術的研究，對我國經濟建設和發展有著特殊的社會和科技價值。

參考文獻：

[1]高剛.降雨入滲條件下邊坡穩定性數值分析[D].蘭州大學，2013.

[2]中國地質環境信息網 全國地質災害通報06-13年度 http：//www.cigem.gov.cn.

[3]宋民崇，余云燕，趙德安.FLAC軟件在邊坡穩定性分析中的應用[J].黑龍江工程學院學報：自然科學版，2009，12（4）.

[4]胡輝.基于Flac3D強度折減理論的邊坡穩定性分析[J].公路與汽運，2012（4）.

[5]李懷玉.錨固巖土邊坡支護效果的評價方法研究[D].長沙理工大學，2010.

[6]賈春玲.邊坡穩定性研究方法[J].中國高新技術企業，2010（10）.