FLAC数值模拟上机题
计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为： 密度ρ＝2500 kg/m3, 弹性模量E＝1×108 Pa，泊松比μ＝0.3，
抗拉强度σt＝0.8×106 Pa，内聚力C＝4.2×104 Pa，摩擦角φ＝17°，膨胀角Δ＝20°
试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。 附换算公式：
1 kN/m3= 100 kg/m3
剪切弹性模量：G?3D
E
2(1??)E3(1?2?)
体积弹性模量：K?1.坡角为30°时
图1 倾角为30°的边坡(单位：m)
输入源代码： set log on
setlogfile stability30.log genzon brick &
p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10 genzon brick &
p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 74.64 0 60 p4 100 2 40 & p5 74.64 2 60 p6 100 0 60 p7 100 2 60 size 30 1 10
modelmohr
Plot Show
Add surface blue Add axes yellow
prop dens 2500 bulk 9.26e7 shear 3.79e7 coh 4.2e4 ten 8e5 friction 17 dilation 20 fix x y z range z -0.1 0.1 fix x range x -0.1 0.1 fix x range x 99.9 100.1 fix y
set gravity 0 0 -10
defcalfos minf=0 maxf=2 loop while maxf-minf>0.01 fs=(maxf+minf)/2.0 refric=atan(0.31/fs)*180/3.14 recoh=42000/fs command inisxx 0.0 syy 0.0 szz 0.0 sxy 0.0 sxz 0.0 szz 0.0 inixvel 0.0 yvel 0.0 zvel 0.0 inixdis 0.0 ydis 0.0 zdis 0.0 profricrefriccohrecoh setmech ratio 1e-5 solve step 5000 printfs end_command aa=mech_ratio ifaa<1e-5 then minf=fs else maxf=fs end_if end_loop end calfos save 2.sav
输入指令显示图形
可以得到fs=1.476，边坡是稳定的。 2.坡角为45°时 修改对应的尺寸即可 genzon brick &
p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10 genzon brick &
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 2 60 size 30 1 10