计算命令流如下： new

gen zone bri p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10

gen zone bri p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 60 0 60 p4 100 2 40 p5 60 2 60 p6 100 0 60& p7 100 1 60 size 30 1 10 model mohr

prop density 2500.0 bulk 8.3e7 shear 3.8e7 c 42000.0 tens 0.8E6 & friction 17 dilation 20 fix x y z range z -0.1 0.1 fix x range x 99.9 100.1 fix x range x -0.1 0.1 fix y

set gravity = 10.0 plot add axe red plot con dis

;定义循环终止条件 def calfos ait1=0.02 k11=0.6 k12=2.0 loop while (k12-k11)>ait1 fs=(k12+k11)/2.0

refric=(atan((tan(17*pi/180))/fs))*180/pi recoh=42000/fs ;折减实现过程 command ini sxx 0.0 syy 0.0 szz 0.0 sxy 0.0 sxz 0.0 syz 0.0 ini xvel 0.0 yvel 0.0 zvel 0.0 ini xdis 0.0 ydis 0.0 zdis 0.0 pro fric refric coh recoh set mech ratio 1e-5 solve step 5000 print fs end_command aa=mech_ratio if aa<1e-5 then k11=fs else k12=fs end_if end_loop end calfos

SAVE results.sav

6

solve fos file slop45.sav associated 计算结果如下：

图2-a，网格剖分图 图2-b，速度矢量图

图2-c，速度等值线图 图2-d，位移等值线图

最终计算边坡的稳定性系数为：Fs＝1.15625

分析：

3D

45°边坡稳定性系数采用的是FLAC自编的强度折减法求解，稳定性系数 1.16>1，从稳定性系数系数可以判断该边坡处于安全状态。坡面最大速度为0.64mm/s，随着深度的增加，竖向应力逐渐增大。坡肩处出现下沉，最大值达到5,34m。

48.45404010060

7

60°

图3 倾角为60°的边坡(单位：m)

计算命令流如下： new

gen zone brick p0 0 0 0 p1 100 0 0 p2 0 1 0 p3 0 0 40 size 50 1 10

gen zone brick p0 40 0 40 p1 100 0 40 p2 40 1 40 p3 51.55 0 60 p4 100 1 40 p5 51.55 1 60& p6 100 0 60 p7 100 1 60 size 30 1 10 attach face range z 39.9 40.1 fix x range x -0.1 0.1 fix x range x 99.9 100.1 fix y

fix z range z 0.1 -0.1 model elas

prop density 2500 bulk 3e9 shear 1e9 set gravity 0 0 -10 solve

ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0 model mohr

prop density 2500 bulk 8.33e7 shear 3.846e7 c 42000 fric 17 ten 800000 solve fos file slope60.sav associated 计算结果如下：

图3-a，网格剖分图 图3-b，速度矢量图

图3-c，速度等值线图 图3-d，位移等值线图

最终计算边坡的稳定性系数为：Fs＝0.94

8

分析：

3D

60°边坡稳定性系数采用的是FLAC内置的强度折减法求解，稳定性系数0.94<1，从稳定性系数系数可以判断该边坡处于不安全状态。坡面最大速度为0.576mm/s，随着深度的增加，竖向应力逐渐增大。坡肩处出现下沉，最大值达到4.14m，故应采取措施以保证边坡安全。

48.4540开挖后坡面原始坡面34010060 图4 边坡开挖算例分析 计算命令流如下： new

gen zone brick p0 0 0 0 p1 100 0 0 p2 0 1 0 p3 0 0 40 size 100 1 10

gen zone brick p0 40 0 40 p1 100 0 40 p2 40 1 40 p3 50 0 50 p4 100 1 40 p5 50 1 50 p6 100 0 50 p7 100 1 50 size 60 1 5

gen zone brick p0 53 0 50 p1 100 0 50 p2 53 1 50 p3 63 0 60 p4 100 1 50 p5 63 1 60 p6 100 0 60 p7 100 1 60 size 47 1 5

attach face range z 39.9 40.1 attach face range z 49.9 50.1 fix x range x -0.1 0.1 fix x range x 99.9 100.1 fix y

fix z range z 0.1 -0.1 model elas

prop density 2500 bulk 3e9 shear 1e9 set gravity 0 0 -10 solve

ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0 model mohr

prop density 2500 bulk 8.33e7 shear 3.846e7 c 42000 fric 17 ten 800000

9

60°45°45°

solve fos file slope60kaiwa.sav associated 计算结果如下：

图4-a，网格剖分图 图4-b，速度矢量图

图4-c，速度等值线图 图4-d，位移等值线图

最终计算边坡的稳定性系数为：Fs＝1.26

分析：

3D

60°开挖后边坡稳定性系数采用的是FLAC内置的强度折减法求解，稳定性系数1.26>1，从稳定性系数系数可以判断该边坡处于安全状态。坡面最大速度为0.51mm/s，随着深度的增加，竖向应力逐渐增大。坡肩处出现下沉，最大值达到3.24m。

10