FLAC3D数值模拟上机报告.doc

中国地质大学（武汉）China University of GeosciencesFLAC3D数值模拟上机报告计算模型分别如图1、2、3所示，边坡倾角分别为30、45、60，岩土体参数为：密度2500 kg/m3, 弹性模量E1108 Pa，泊松比0.3，抗拉强度t0.8106 Pa，内聚力C4.2104 Pa，摩擦角17，膨胀角20试用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 10gen 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;=;设置边界条件fix x y z range z -0.1 0.1fix x range x 99.9 100.1fix x range x -0.1 0.1fix y;=;初始地应力的生成model elasprop density 2000 bulk 3e9 shear 1e9set gravity 0 0 -10solveini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0;=;安全系数求解model mohrprop density 2500.0 bulk 8.3E7 shear 3.8E7 coh 42000.0 tens 0.8E6 friction 17 dilation 20solve fos file slope3dfos1.sav associated安全系数：最终计算边坡稳定性系数为1.45310数值模拟上机报告图1 网格剖分图 图2 速度矢量图图3 速度等值线图 图4 位移等值线图图5 剪应变增量云图二 坡度为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;=;设置边界条件fix x y z range z -0.1 0.1fix x range x 99.9 100.1fix x range x -0.1 0.1fix y;=;初始地应力的生成model elasprop density 2000 bulk 3e9 shear 1e9set gravity 0 0 -10solveini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0;=;安全系数求解model mohrprop density 2500.0 bulk 8.3E7 shear 3.8E7 coh 42000.0 tens 0.8E6 friction 17 dilation 20solve fos file slope3dfos1.sav associated安全系数：最终边坡的稳定性系数为1.14图1 网格剖分图 图2 速度矢量图图3 速度等值线图 图4 位移等值线图图5 剪应变增量云图三 坡度为60的情况图3 倾角为60的边坡(单位：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 51.55 0 60 p4 100 2 40 & p5 51.55 2 60 p6 100 0 60 p7 100 2 60 size 30 1 10;=;设置边界条件fix x y z range z -0.1 0.1fix x range x 99.9 100.1fix x range x -0.1 0.1fix y;=;初始地应力的生成model elasprop density 2000 bulk 3e9 shear 1e9set gravity 0 0 -10solveini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0;=;安全系数求解model mohrprop density 2500.0 bulk 8.3E7 shear 3.8E7 coh 42000.0 tens 0.8E6 friction 17 dilation 20solve fos file slope3dfos1.sav associated安全系数：最终边坡的稳定性系数为0.928 图1 网格剖分图 图2 速度矢量图图3 速度等值线图 图4 位移等值线图 图5 剪应变增量云图四 坡度为60的边坡开挖情况算例分析：命令流： new;= 建立网格模型 gen zone brick p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10gen zone brick p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 50 0 50 p4 100 2 40 p5 50 2 50 p6 100 &0 50 p7 100 2 50 size 30 1 10gen zone brick p0 53 0 50 p1 100 0 50 p2 53 2 50 p3 63 0 60 p4 100 2 50 p5 63 2 60 p6 100& 0 60 p7 100 2 60 size 15 1 10gen zone brick p0 45.77 0 50 p1 53 0 50 p2 45.77 2 50 p3 51.55 0 60 p4 53 2 50 p5 51.55 2 &60 p6 63 0 60 p7 63 2 60 size 15 1 10 group exc1gen zone wedge p0 40 0 40 p1 50 0 50 p2 40 2 40 p3 45.77 0 50 p4 50 2 50 p5 45.77 2 50 &size 30 1 10 group exc2group section1 range y 0 2 group exc1group section2 range y 0 2 group exc2attach face;=;定义本构模型mod elaspro density 2500 bulk 0.83e8 she 0.38e8 ;=;设置边界条件fix x y z range z -.1 .1fix x range x 99.9 100.1fix x range x -0.1 0.1fix y;=;设置重力加速度set gravity 0 0 -10.0;= ;设定初始条件 ini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0 ;= ;开挖 mod mohrpro density 2500 bulk 0.83e8 she 0.38e8 fric 17 coh 4.2e4 ten 0.8e