边坡稳定性分析的数值模拟

1.FLAC3D数值模拟上机题计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为：密度ρ＝2500kg/m3,弹性模量E＝1×108Pa，泊松比μ＝0.3，抗拉强度σt＝0.8×106Pa，内聚力C＝4.2×104Pa，摩擦角φ＝17°，膨胀角Δ＝20°。试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。附换算公式：1kN/m3=100kg/m3剪切弹性模量：)1(2EG；体积弹性模量：)21(3EK4025.36604010030°1004060404045°图1倾角为30°的边坡(单位：m)图2倾角为45°的边坡(单位：m)100406048.454060°图3倾角为60°的边坡(单位：m)实例分析：1）坡角为30°时的边坡情况：计算代码（模式）：new；开始一个新的分析genzonebrickp0000p110000p2020p30040&size50110；生成下面的矩形，沿x、y、z三房向分为50，1，10分genzonebrick&p040040p1100040p240240p374.64060&p4100240p574.64260p6100060p7100260&size30110；生成上面的梯形，沿x、y、z三房向分为30，1，10分fixzrangez-0.10.1；固定模型底面fixxrangex-0.10.1；固定模型左面fixxrangex99.9100.1；固定模型右面fixyrangey-0.10.1；固定模型前面fixyrangey1.92.1；固定模型后面modelmohr；库伦摩尔模型propcoh=4.2e4ten=8e5fric=17；力学参数赋值inidens=2500setgra=0,0,-9.8；重力设置propbulk8.3e7shear3.85e7inizvel0；z方向初始速度为0inixdisp0ydisp0zdisp0；xyz方向初始位移为0plotcreateslope；创建一个斜坡plotaddaxes；添加坐标轴plotaddblockplotshowsolvefosfileslope3dfos.savassociated强度折减法求解FLAC3D3.00ItascaConsultingGroup,Inc.Minneapolis,MNUSASettings:ModelPerspective16:50:15SatJun072008Center:X:5.000e+001Y:1.000e+000Z:3.000e+001Rotation:X:0.000Y:0.000Z:0.000Dist:2.775e+002Mag.:1Ang.:22.500SurfaceMagfac=0.000e+000FLAC3D3.00ItascaConsultingGroup,Inc.Minneapolis,MNUSAStep25701ModelPerspective22:14:18SatJun072008Center:X:5.000e+001Y:1.000e+000Z:3.000e+001Rotation:X:0.000Y:0.000Z:0.000Dist:2.775e+002Mag.:1Ang.:22.500SurfaceMagfac=0.000e+000VelocityMaximum=4.905e-007Linestyle图4网格剖分图图5速度矢量图FLAC3D3.00ItascaConsultingGroup,Inc.Minneapolis,MNUSAStep25701ModelPerspective22:17:17SatJun072008Center:X:5.000e+001Y:1.000e+000Z:3.000e+001Rotation:X:0.000Y:0.000Z:0.000Dist:2.775e+002Mag.:1Ang.:22.500SurfaceMagfac=0.000e+000ContourofVelocityMag.Magfac=0.000e+0000.0000e+000to5.0000e-0085.0000e-008to1.0000e-0071.0000e-007to1.5000e-0071.5000e-007to2.0000e-0072.0000e-007to2.5000e-0072.5000e-007to3.0000e-0073.0000e-007to3.5000e-0073.5000e-007to4.0000e-0074.0000e-007to4.5000e-0074.5000e-007to4.9050e-007Interval=5.0e-008FLAC3D3.00ItascaConsultingGroup,Inc.Minneapolis,MNUSAStep25701ModelPerspective22:20:36SatJun072008Center:X:5.000e+001Y:1.000e+000Z:3.000e+001Rotation:X:0.000Y:0.000Z:0.000Dist:2.775e+002Mag.:1Ang.:22.500SurfaceMagfac=0.000e+000ContourofDisplacementMag.Magfac=0.000e+0000.0000e+000to5.0000e-0025.0000e-002to1.0000e-0011.0000e-001to1.5000e-0011.5000e-001to2.0000e-0012.0000e-001to2.5000e-0012.5000e-001to3.0000e-0013.0000e-001to3.2479e-001Interval=5.0e-002图6速度等值线图图7位移等值线图最终计算边坡的稳定性系数为：Fs＝1.472）坡角为45°时的边坡情况：代码：newgenzonebrickp0000p110000p2020p30040size50110genzonebrick&p040040p1100040p240240p360060&p4100240p560260p6100060p7100260&；建立模型size30110fixzrangez-0.10.1；固定底面fixxrangex-0.10.1；固定左面fixxrangex99.9100.1；固定右面fixyrangey-0.10.1；固定前面fixyrangey1.92.1；固定后面modelmohr；摩尔库伦模型propcoh=4.2e4ten=8e5fric=17；参数赋值inidens=2500setgra=0,0,-9.8propbulk8.3e7shear3.85e7inizvel0rangez060y02x0100；初始速度为0plotcreatexxx；创建一个名为xxx的新视图plotaddaxes；添加坐标轴plotaddblock；根据不同的模型变量用不同的颜色绘出单元体面plotshow；屏幕上显示当前视图solvefosassociated自动查找安全因子，实施关联流动规则即膨胀角等于摩擦角solvefosfileslope3dfos.sav；前solvefos为自动查找安全因子，后半为把最后不平衡力写进指定的文件名中这最后两句可以一次写完：solvefosfileslope3dfos.savassociated图8网格剖分图图9速度矢量图图10速度等值线图图11位移等值线图最终计算边坡的稳定性系数为：Fs＝1.133）坡角为60°时的边坡：代码：newgenzonebrickp0000p110000p2020p30040size50110genzonebrick&p040040p1100040p240240p351.55060&p4100240p551.55260p6100060p7100260&size30110；创建模型fixzrangez-0.10.1；固定底面fixxrangex-0.10.1；固定左面fixxrangex99.9100.1；固定右面fixyrangey-0.10.1；固定前面fixyrangey1.92.1；固定后面modelmohr；摩尔库伦模型propcoh=4.2e4ten=8e5fric=17；参数赋值inidens=2500setgra=0,0,-9.8propbulk8.3e7shear3.85e7inizvel0rangez060y02x0100；初始速度为0plotcreateslope；创建一个斜坡plotaddaxes；添加坐标轴plotaddblockplotshowsolvefosfileslope3dfos.savassociated；强度折减法求解图12网格剖分图图13速度矢量图FLAC3D3.00ItascaConsultingGroup,Inc.Minneapolis,MNUSAStep29011ModelPerspective09:59:41SunJun082008Center:X:5.000e+001Y:1.000e+000Z:3.000e+001Rotation:X:0.000Y:0.000Z:0.000Dist:2.775e+002Mag.:1Ang.:22.500SurfaceMagfac=0.000e+000ContourofDisplacementMag.Magfac=0.000e+0000.0000e+000to1.0000e-0011.0000e-001to2.0000e-0012.0000e-001to3.0000e-0013.0000e-001to4.0000e-0014.0000e-001to5.0000e-0015.0000e-001to6.0000e-0016.0000e-001to7.0000e-0017.0000e-001to8.0000e-0018.0000e-001to9.0000e-0019.0000e-001to1.0000e+0001.0000e+000to1.0562e+000Interval=1.0e-001图14速度等值线图图15位移等值线图最终计算边坡的稳定性系数为：Fs＝0.94。4）模拟开挖过程：选取坡角为60°情况下，模拟边坡开挖过程中边坡稳定性系数的变化。网格剖分图如下图16，拟要开挖有二个部分，每次开挖后计算的边坡稳定性系数Fs=1.24。//建立网格模型//newgenzonewedgep051.55260p160260p251.55060p340240p460060p540040size8120group1；开挖块体一genzonebrickp050050p153050p250250p360060p453250p560260p663060p763260size3110group2；开挖块体二genzonebrickp053050p1100050p253250p363060p4100250p563260p6100060p7100260size47110genzonebrickp040040p1100040p240240p350050p4100240p550250p6100050p7100250size50110genzonebrickp0-0.800p110000p2-0.820p3-0.8040p410020p5-0.8240p6100040p7100240size84120；创建模型fixzrangez-0.10.1；固定底面fixxrangex-0.9-0.8；固定左面fixxrangex99.9100.1；固定右面fixy；固定y方向modelelas；弹塑性模型propcoh=4.2e10ten=8e10fric=17；参数赋值propdens=2500modelmohr；库伦摩尔模型propcoh