中国地质大学数值模拟上机报告教程.doc

FLAC3D数值模拟上机报告 计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为： 容重r＝2500 kg/m3, 弹性模量E＝1×108 Pa，泊松比μ＝0.3， 抗拉强度σt＝0.8×106 Pa，内聚力C＝4.2×104 Pa，摩擦角φ＝17° 试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。 附 换算公式： 1 kN/m3= 100 kg/m3 剪切弹性模量： 体积弹性模量： 一、 坡度为30°的情况 图1 倾角为30°的边坡(单位：m) 算例分析： 命令流： new ;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞ ;建立网格模型 gen zone brick p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10 gen zone 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 ; ＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞ ;定义本构模型 model mohr prop density 2500.0 bulk 8.3E7 shear 3.8E7 coh 42000.0 tens 0.8E6 friction 17 ; ＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞ ;设置边界条件 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 ; ＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞ ;设定初始条件 ini sxx 0.0 syy 0.0 szz 0.0 sxy 0.0 sxz 0.0 szz 0.0 ini xvel 0.0 yvel 0.0 zvel 0.0 ini xdis 0.0 ydis 0.0 zdis 0.0 ; ＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞ ;安全系数求解 def calfos minf=0.1 maxf=2.2 loop while maxf-minf0.01 fs=(maxf+minf)/2.0 refric=atan(0.30/fs)*180/3.14 recoh=42000/fs command pro fric refric coh recoh set mech ratio 1e-5 solve step 5000 print fs end_command aa=mech_ratio if aa1e-5 then minf=fs else maxf=fs end_if end_loop end calfos  图1 网格剖分图 图2 速度矢量图 图3 速度等值线图 图4 位移等值线图 图5 剪应变增量云图 图6 最小主应力云图 图7 最大主应力云图 安全系数： 因此，最终计算边坡稳定性系数为1.453516 二、 坡度为45°的情况 图2 倾角为45°的边坡(单位：m) 算例分析： 命令流： new ; ＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞ ;建立网格模型 gen zone brick p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10 gen zone 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 ; ＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞ ;定义本构模型 model mohr prop density 2500.0 bul