Geostudio及FLAC3D.docx
《Geostudio及FLAC3D.docx》由会员分享，可在线阅读，更多相关《Geostudio及FLAC3D.docx（22页珍藏版）》请在冰豆网上搜索。
Geostudio及FLAC3D
Geostudio及FLAC3D
《地质工程数值法》
GeoStudio上机实习作业
内容：
用GeoStudio的Seep/w模块模拟滑坡不同工况下的地下水渗流曲线。
要求：
提供模拟结果文件，现场验收操作。
相关资料：
滑坡工程地质剖面图；
岩土参数；
渗流模拟工况。
1计算模型的建立
根据2#滑坡的工程地质条件，选择该滑坡钻探剖面（图1）作为模拟剖面。
其方位NW113°，滑坡前缘高程1190m，后缘高程1398m，前缘和后部较薄，中部较厚，根据钻孔揭露，滑体物质主要是灰黑色碎块石土，滑带土主要为灰黒色含角砾粉质粘土、角砾质土或含角砾泥炭质土，滑带厚度0.6～1.5m，滑床为炭质页岩和炭质灰岩的碎块石土。
基于以上条件，对该堆积层滑坡的坡体地质模型进行网格剖分，共剖分为2829个网格单元，如图2所示。
图12#滑坡工程地质剖面图
操作步骤：
1.打开Geostudio：
2.操作界面尺寸设置：
3.导入图形
4.检查图形
5.设置材料属性
6.将材料属性赋予给图形区域
7.划分网格
8.设置边界条件，修改颜色
9.赋予边界条件
10.检查模型是否正确
11.保存，并solve
12.查看结果
与预期结果相一致，上层碎石土的渗流力明显大于滑床的。
2计算参数与工况
参数的取值主要由勘查资料中的岩土体物理力学试验、水文地质试验、工程地质类比以及经验值等综合确定。
具体为：
模型材料的天然、饱和容重等由渗透试验确定；饱和体积含水量和饱和渗透系数，类比附近滑坡或者同一类型滑坡的物质组成及坡体结构等，确定出参数取值范围，然后根据野外钻探的实际资料，通过参数调整，反演出2#滑坡的实际地下水位，获取后续模拟的参数，渗流场模拟过程中所涉及到的各层岩土体参数见表1。
采用GEO-Studio软件中的SEEP/W模块进行滑坡渗流场模拟时，视滑坡岩土体材料渗透系数和体积含水量为坡体孔隙水压力的函数。
当岩土体处于饱和状态时，渗透系数和体积含水量为定值；而当岩土体处于非饱和状态时其渗透系数和体积含水量应该通过试验来确定。
但鉴于试验的难度，一般取经验值。
本节进行2#滑坡的渗流场模拟时采用SEEP/W模块自带的VanGenuchten经验曲线和饱和状态时的参数来确定当坡体材料处于非饱和状态时的渗透系数和体积含水量与孔隙水压力的函数关系，分别如图3~图8所示。
渗流模拟工况：
一、自重+地下水；二、自重+50年一遇15日连续强降雨工况（50年一遇降雨强度为24小时最大值130.0mm）。
表12#滑坡渗流场模拟各层岩土体参数选取
岩性
天然
重度
（kN/m3）
饱和
重度
（kN/m3）
饱和体积
含水量（%）
饱和渗透
系数
（×10-7m/s）
抗剪强度指标
天然
饱和
C（KPa）
φ（°）
C（KPa）
φ（°）
滑坡堆积层
19.7
20.5
22.8
6.9
38
33
32
27
滑带
19.7
20.5
23.2
6.9
22.5
22
19.1
20.7
滑床
20.4
21.4
22.4
0.8
60
42
52
36
操作步骤：
1.输入模型
2.输入材料属性
3.输入材料属性（因为模型是直接从seep/w里面导入的，没有取消seep里面使用的划分网格）
4.画地下水位
5.drawslipsurfaceandexitrange
6.查看区域属性
7.画孔隙水压力云图
8.查看数据
9.Solve
可以看到安全系数为1.394
10.极限平衡法对应的每个条块受力图
11.查看不同的方法对应的不同安全系数
12.画孔隙水压力图
13.输出滑坡稳定性报告
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.坡角为30°时
图1倾角为30°的边坡（单位：
m）
输入源代码：
setlogon
setlogfilestability30.log
genzonbrick&
p0000p110000p2020p30040size50110
genzonbrick&
p040040p1100040p240240p374.64060p4100240&
p574.64260p6100060p7100260size30110
modelmohr
Plot
Show
Addsurfaceblue
Addaxesyellow
propdens2500bulk9.26e7shear3.79e7coh4.2e4ten8e5friction17dilation20
fixxyzrangez-0.10.1
fixxrangex-0.10.1
fixxrangex99.9100.1
fixy
setgravity00-10
defcalfos
minf=0
maxf=2
loopwhilemaxf-minf>0.01
fs=（maxf+minf）/2.0
refric=atan（0.31/fs）*180/3.14
recoh=42000/fs
command
inisxx0.0syy0.0szz0.0sxy0.0sxz0.0szz0.0
inixvel0.0yvel0.0zvel0.0
inixdis0.0ydis0.0zdis0.0
profricrefriccohrecoh
setmechratio1e-5
solvestep5000
printfs
end_command
aa=mech_ratio
ifaa<1e-5then
minf=fs
else
maxf=fs
end_if
end_loop
end
calfos
save2.sav
输入指令显示图形
可以得到fs=1.476，边坡是稳定的。
2.坡角为45°时
修改对应的尺寸即可
genzonbrick&
p0000p110000p2020p30040size50110
genzonbrick&
p040040p1100040p240240p360060p4100240&
p560260p6100060p7100260size30110
升级会员 