FLAC3D数值模拟上机报告计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为：容重r＝2500kg/m3,弹性模量E＝1×108Pa，泊松比μ＝0.3，抗拉强度σt＝0.8×106Pa，内聚力C＝4.2×104Pa，摩擦角φ＝17°试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。附换算公式：1kN/m3=100kg/m3剪切弹性模量：881100.38510()2(1)2(10.3)EGPa体积弹性模量：881100.83310()3(12)3(120.3)EKPa一、坡度为30°的情况4025.36604010030°图1倾角为30°的边坡(单位：m)算例分析：命令流：new;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;建立网格模型genzonebrickp0000p110000p2020p30040size50110genzonebrickp040040p1100040p240240p374.64060p4100240&p574.64260p6100060p7100260size30110;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;定义本构模型modelmohrpropdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;设置边界条件fixxyzrangez-0.10.1fixxrangex99.9100.1fixxrangex-0.10.1fixy;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;设置重力加速度setgravity=10.0;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;设定初始条件inisxx0.0syy0.0szz0.0sxy0.0sxz0.0szz0.0inixvel0.0yvel0.0zvel0.0inixdis0.0ydis0.0zdis0.0;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;安全系数求解defcalfosminf=0.1maxf=2.2loopwhilemaxf-minf>0.01fs=(maxf+minf)/2.0refric=atan(0.30/fs)*180/3.14recoh=42000/fscommandprofricrefriccohrecohsetmechratio1e-5solvestep5000printfsend_commandaa=mech_ratioifaa<1e-5thenminf=fselsemaxf=fsend_ifend_loopendcalfos图1网格剖分图图2速度矢量图图3速度等值线图图4位移等值线图图5剪应变增量云图图6最小主应力云图图7最大主应力云图安全系数：因此，最终计算边坡稳定性系数为1.453516二、坡度为45°的情况1004060404045°图2倾角为45°的边坡(单位：m)算例分析：命令流：new;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;建立网格模型genzonebrickp0000p110000p2020p30040size50110genzonebrickp040040p1100040p240240p360060p4100240&p560260p6100060p7100260size30110;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;定义本构模型modelmohrpropdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;设置边界条件fixxyzrangez-0.10.1fixxrangex99.9100.1fixxrangex-0.10.1fixy;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;设置重力加速度setgravity=10.0;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;设定初始条件inisxx0.0syy0.0szz0.0sxy0.0sxz0.0szz0.0inixvel0.0yvel0.0zvel0.0inixdis0.0ydis0.0zdis0.0;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;安全系数求解defcalfosminf=0.1maxf=2.2loopwhilemaxf-minf>0.01fs=(maxf+minf)/2.0refric=atan(0.30/fs)*180/3.14recoh=42000/fscommandprofricrefriccohrecohsetmechratio1e-5solvestep5000printfsend_commandaa=mech_ratioifaa<1e-5thenminf=fselsemaxf=fsend_ifend_loopendcalfos图1网格剖分图图2速度矢量图图3速度等值线图图4位移等值线图图5剪应变增量云图图6最小主应力云图图7最大主应力云图安全系数：因此，最终边坡的稳定性系数为1.10894三、坡度为60°的情况100406048.454060°图3倾角为60°的边坡(单位：m)算例分析：命令流：new;＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞;建立网格模型genzonebrickp000...
