摘要：FLAC 3D是當(dāng)前的在巖土工程領(lǐng)域應(yīng)用較普遍的數(shù)值分析軟件，在邊坡穩(wěn)定性分析時(shí)具有強(qiáng)大的計(jì)算功能，在建立模型通過不斷單一的改變不同的巖土體的參數(shù)值，得出不同巖土體參數(shù)對(duì)滑坡面的影響程度。

關(guān)鍵詞：FLAC3D;安全系數(shù);滑坡面;數(shù)值模擬

1 概述

由于巖土體介質(zhì)的復(fù)雜性，分析邊界和工況的多變，用經(jīng)典力學(xué)求解復(fù)雜的巖體工程問題的解析解已經(jīng)很難，甚至不可能。隨著電子計(jì)算機(jī)的廣泛使用，數(shù)值分析方法在巖土工程、工程地質(zhì)領(lǐng)域獲得了迅速的發(fā)展，為研究巖土工程和工程地質(zhì)問題提供了強(qiáng)有力的手段。數(shù)值模擬理論和方法發(fā)展迅速，工程地質(zhì)數(shù)值方法的不斷成熟和完善，使得解決工程地質(zhì)問題更加廣泛，研究的課題更加深入。一方面，飛速發(fā)展的工程地質(zhì)學(xué)不斷地提出新的難題，用現(xiàn)成的數(shù)學(xué)、力學(xué)理論對(duì)其無法作出確切的描述，工程地質(zhì)數(shù)值方法為解決這類問題提供了可能的手段;另一方面，各種數(shù)值方法的不斷成功應(yīng)用，深化了人們對(duì)許多工程地質(zhì)現(xiàn)象的理解，并有力地推動(dòng)了工程地質(zhì)學(xué)科的定量化進(jìn)程。

FLAC（Fast Lagrangian Analysis of Continua）數(shù)值分析軟件是為滿足巖土工程、隧道工程、采礦工程、道路與鐵道工程等學(xué)科領(lǐng)域研究的有限差分計(jì)算軟件，目前，該軟件已得到全世界 70 多個(gè)國家廣泛應(yīng)用，它有兩個(gè)計(jì)算程序版本即二維和三維兩個(gè)計(jì)算程序。FLAC3D（Three Dimensional Fast Lagrangian Analysis of Continua）是由美國明尼蘇達(dá)Itasca Consulting Goup Inc. 軟件公司開發(fā)的通用程序，它是 FLAC2D 在三維空間上的拓展，可以用于模擬巖石、土質(zhì)以及其它材料的三維結(jié)構(gòu)受力特性，尤其用于分析達(dá)到屈服極限強(qiáng)度時(shí)的塑性流動(dòng)，廣泛應(yīng)用于支護(hù)設(shè)計(jì)和評(píng)價(jià)，邊坡穩(wěn)定性分析、施工設(shè)計(jì)（模擬開挖和填筑）、拱壩穩(wěn)定性分析、河谷演化進(jìn)程再現(xiàn)等。本文通過建立數(shù)值模型作為分析平臺(tái)，在未支護(hù)狀態(tài)下，采用強(qiáng)度折減法分析邊坡在達(dá)到臨界破壞時(shí)，所產(chǎn)生的塑性區(qū)、位移等值圖、剪應(yīng)變邊分布圖等來確定滑動(dòng)面的分布情況。通過不斷的改變巖、土體的參數(shù)值（粘聚力、內(nèi)摩擦力、彈性模量等），找出不同參數(shù)對(duì)于滑動(dòng)的面的影響，以及敏感程度。

2 計(jì)算模型

采用如下模型分析：邊坡高15米，坡度1：1，坡腳到左側(cè)的距離為20米，坡頂?shù)接覀?cè)邊界的的距離為40米，坡腳向下方延伸15米深。

3計(jì)算結(jié)果處理：

通過不斷的改變單一的參數(shù)值，得到不同的滑動(dòng)面，借用經(jīng)濟(jì)上的敏感性因素分析法，來確定不同參數(shù)對(duì)滑動(dòng)面的影響情況。

當(dāng)改變粘聚力時(shí)，得到粘聚力和安全系數(shù)有如下關(guān)系見表（一）

表（一）

粘聚力與安全系數(shù)的的關(guān)系

粘聚力

安全系數(shù)

粘聚力改變比

安全系數(shù)變化比

0

0.4721

10

0.7328

100

55.22

20

0.9941

100

110.57

30

1.3119

100

177.89

40

1.6223

100

243.63

經(jīng)線性擬合，如圖（2）所示。

如圖（2）

改變內(nèi)摩擦角時(shí)，得到內(nèi)摩擦角和安全系數(shù)有如下關(guān)系見表（二）

表（二）

內(nèi)摩擦角與安全系數(shù)的的關(guān)系

內(nèi)摩擦角

安全系數(shù)

內(nèi)摩擦角增加率

安全系數(shù)變化率

10

0.7832

15

0.9222

50

17.75

20

1.1961

100

52.72

25

1.3371

150

70.72

30

1.6127

200

105.91

經(jīng)線性擬合，如圖（3）所示。

如圖（3）

3）、改變彈性模量時(shí)，得到彈性模量和安全系數(shù)有如下關(guān)系見表（三）

表（三）

彈性模量與安全系數(shù)的的關(guān)系

彈性模量

安全系數(shù)

彈性模量改變率

安全系數(shù)變化率

0

1.2654

10

1.2655

100

0.01

20

1.2659

200

0.04

30

1.2652

300

-0.02

40

1.2655

400

0.01

從表（三）和圖（三）都可以清楚的知道當(dāng)彈性模量的變化達(dá)到百分之百是安全系數(shù)的變化率還不足百分之零點(diǎn)一。

4）、改變抗拉強(qiáng)度時(shí)，得到抗拉強(qiáng)度和安全系數(shù)有如下關(guān)系見表（四）

表（四）

抗拉強(qiáng)度與安全系數(shù)的的關(guān)系

抗拉強(qiáng)度

安全系數(shù)

抗拉強(qiáng)度增加率

安全系數(shù)變化率

50

1.2123

100

1.2097

100

-0.21

150

1.2158

200

0.29

200

1.212

300

-0.02

250

1.2178

400

0.45

繪抗拉強(qiáng)度增量的改變率和安全系數(shù)變化率的關(guān)系如圖（三）所示，

從表（四）和圖（四）都可以清楚的知道當(dāng)彈性模量的變化達(dá)到百分之百是安全系數(shù)的變化率還不足百分之一。

4 計(jì)算結(jié)果分析：

粘聚力增加率和安全系數(shù)變化率的關(guān)系的線性擬合方程為 ，其中R值為0.9991。內(nèi)摩擦力的增加率和安全系數(shù)變化率的線性擬合方程為 ，其中R值為0.9928。粘聚力和內(nèi)摩擦角的增加率與安全系數(shù)的變化率均符合線性的變化，且粘聚力的斜率要比內(nèi)摩擦角的斜率大?？估瓘?qiáng)度和彈性模量的增加率與安全系數(shù)的變化則沒有明顯的變化，均在百分之一以內(nèi)。從影響因素的敏感性角度來講粘聚力的變化對(duì)安全系數(shù)的影響最大，內(nèi)摩擦角次之，抗拉強(qiáng)度和彈性模量基本不產(chǎn)生影響。

5結(jié)語

FLAC 3D是當(dāng)前的在巖土工程領(lǐng)域應(yīng)用較普遍的數(shù)值分析軟件，本文以該軟件為基礎(chǔ)模擬均質(zhì)土體邊坡滑坡面的安全系數(shù)，通過改變不同的巖土體的參數(shù)最終的出不同參數(shù)對(duì)滑坡面的影響程度，為應(yīng)用該軟件從事巖土設(shè)計(jì)提供參考。

參考文獻(xiàn)：

[1] 陳育民，劉漢龍. FLAC/FLAC3D基礎(chǔ)與工程實(shí)例[M].北京：中國水利水電出版社，2013.

[2] 張書偉，林" 杭，任連偉.FLAC3D在巖土工程中的應(yīng)用[M].北京：中國水利水電出版社，2011.

[3] 張玉燈.FLAC3D在巖質(zhì)邊坡穩(wěn)定性分析中的應(yīng)用[J].路基工程，2008（6）：164—165.

[4] 王仲原.FLAC3D在邊坡中的應(yīng)用分析[J].科技資訊，2009（16）：75—76.

[5] 陳勇河，周德泉.用FLAC-3D分析邊坡的穩(wěn)定性[J].山西建筑，2008 34（22）：18—19.