FLAC3D数值模拟上机报告材料.docx
《FLAC3D数值模拟上机报告材料.docx》由会员分享,可在线阅读,更多相关《FLAC3D数值模拟上机报告材料.docx(15页珍藏版)》请在冰豆网上搜索。
FLAC3D数值模拟上机报告材料
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
剪切弹性模量:
体积弹性模量:
一坡度为30°的情况
图1倾角为30°的边坡(单位:
m)
算例分析:
命令流:
new
;==========================================================
建立网格模型
genzonebrickp0000p110000p2020p30040size50110
genzonebrickp040040p1100040p240240p374.64060p4100240&
p574.64260p6100060p7100260size30110
;=========================================================
;设置边界条件
fixxyzrangez-0.10.1
fixxrangex99.9100.1
fixxrangex-0.10.1
fixy
;======================
;初始地应力的生成
modelelas
propdensity2000bulk3e9shear1e9
setgravity00-10
solve
inixdisp0ydisp0zdisp0
inixvel0yvel0zvel0
;===================================
;安全系数求解
modelmohr
propdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17dilation20
solvefosfileslope3dfos1.savassociated
安全系数:
最终计算边坡稳定性系数为1.453
图1网格剖分图图2速度矢量图
图3速度等值线图图4位移等值线图
图5剪应变增量云图
二坡度为45°的情况
图2倾角为45°的边坡(单位:
m)
算例分析:
命令流:
new
;==========================================================
建立网格模型
genzonebrickp0000p110000p2020p30040size50110
genzonebrickp040040p1100040p240240p360060p4100240&
p560260p6100060p7100260size30110
;=========================================================
;设置边界条件
fixxyzrangez-0.10.1
fixxrangex99.9100.1
fixxrangex-0.10.1
fixy
;======================
;初始地应力的生成
modelelas
propdensity2000bulk3e9shear1e9
setgravity00-10
solve
inixdisp0ydisp0zdisp0
inixvel0yvel0zvel0
;===================================
;安全系数求解
modelmohr
propdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17dilation20
solvefosfileslope3dfos1.savassociated
安全系数:
最终边坡的稳定性系数为1.14
图1网格剖分图图2速度矢量图
图3速度等值线图图4位移等值线图
图5剪应变增量云图
三坡度为60°的情况
图3倾角为60°的边坡(单位:
m)
算例分析:
命令流:
new
;==========================================================
建立网格模型
genzonebrickp0000p110000p2020p30040size50110
genzonebrickp040040p1100040p240240p351.55060p4100240&
p551.55260p6100060p7100260size30110
;=========================================================
;设置边界条件
fixxyzrangez-0.10.1
fixxrangex99.9100.1
fixxrangex-0.10.1
fixy
;======================
;初始地应力的生成
modelelas
propdensity2000bulk3e9shear1e9
setgravity00-10
solve
inixdisp0ydisp0zdisp0
inixvel0yvel0zvel0
;===================================
;安全系数求解
modelmohr
propdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17dilation20
solvefosfileslope3dfos1.savassociated
安全系数:
最终边坡的稳定性系数为0.928
图1网格剖分图图2速度矢量图
图3速度等值线图图4位移等值线图
图5剪应变增量云图
四坡度为60°的边坡开挖情况
算例分析:
命令流:
new
;==========================================================
建立网格模型
genzonebrickp0000p110000p2020p30040size50110
genzonebrickp040040p1100040p240240p350050p4100240p550250p6100&050p7100250size30110
genzonebrickp053050p1100050p253250p363060p4100250p563260p6100&060p7100260size15110
genzonebrickp045.77050p153050p245.77250p351.55060p453250p551.552&60p663060p763260size15110groupexc1
genzonewedgep040040p150050p240240p345.77050p450250p545.77250&size30110groupexc2
groupsection1rangey02groupexc1
groupsection2rangey02groupexc2
attachface
;=========================================================
;定义本构模型
modelas
prodensity2500bulk0.83e8she0.38e8
;================================================================
;设置边界条件
fixxyzrangez-.1.1
fixxrangex99.9100.1
fixxrangex-0.10.1
fixy
;======================
;设置重力加速度
setgravity00-10.0
;===================================
;设定初始条件
inixdisp0ydisp0zdisp0
inixvel0yvel0zvel0
;===================================
;开挖
modmohr
prodensity2500bulk0.83e8she0.38e8fric17coh4.2e4ten0.8e6dila20
modnullrangegroupsection1
modnullrangegroupsection2
solvefosfileslope3dfos1.savassociated
安全系数:
最终边坡的稳定性系数为1.36
图1网格剖分图图2速度矢量图
图3速度等值线图图4位移等值线图
图5剪应变增量云图