Geostudio及FLAC3D.doc

地质工程数值法GeoStudio上机实习作业内 容：用GeoStudio的Seep/w模块模拟滑坡不同工况下的地下水渗流曲线。要 求：提供模拟结果文件，现场验收操作。相关资料： 滑坡工程地质剖面图； 岩土参数； 渗流模拟工况。1 计算模型的建立根据2#滑坡的工程地质条件，选择该滑坡钻探剖面（图1）作为模拟剖面。其方位NW113，滑坡前缘高程1190m，后缘高程1398m，前缘和后部较薄，中部较厚，根据钻孔揭露，滑体物质主要是灰黑色碎块石土，滑带土主要为灰黒色含角砾粉质粘土、角砾质土或含角砾泥炭质土，滑带厚度0.61.5m，滑床为炭质页岩和炭质灰岩的碎块石土。基于以上条件，对该堆积层滑坡的坡体地质模型进行网格剖分，共剖分为2829个网格单元，如图2所示。图1 2#滑坡工程地质剖面图操作步骤：1. 打开Geostudio：2. 操作界面尺寸设置：3. 导入图形4. 检查图形5. 设置材料属性6. 将材料属性赋予给图形区域7. 划分网格8. 设置边界条件，修改颜色9. 赋予边界条件10. 检查模型是否正确11. 保存，并solve12. 查看结果与预期结果相一致，上层碎石土的渗流力明显大于滑床的。2 计算参数与工况参数的取值主要由勘查资料中的岩土体物理力学试验、水文地质试验、工程地质类比以及经验值等综合确定。具体为：模型材料的天然、饱和容重等由渗透试验确定；饱和体积含水量和饱和渗透系数，类比附近滑坡或者同一类型滑坡的物质组成及坡体结构等，确定出参数取值范围，然后根据野外钻探的实际资料，通过参数调整，反演出2#滑坡的实际地下水位，获取后续模拟的参数，渗流场模拟过程中所涉及到的各层岩土体参数见表1。采用GEO-Studio软件中的SEEP/W模块进行滑坡渗流场模拟时，视滑坡岩土体材料渗透系数和体积含水量为坡体孔隙水压力的函数。当岩土体处于饱和状态时，渗透系数和体积含水量为定值；而当岩土体处于非饱和状态时其渗透系数和体积含水量应该通过试验来确定。但鉴于试验的难度，一般取经验值。本节进行2#滑坡的渗流场模拟时采用SEEP/W模块自带的Van Genuchten经验曲线和饱和状态时的参数来确定当坡体材料处于非饱和状态时的渗透系数和体积含水量与孔隙水压力的函数关系，分别如图3图8所示。渗流模拟工况：一、自重+地下水；二、自重+50年一遇15日连续强降雨工况（50年一遇降雨强度为24小时最大值130.0mm）。表1 2#滑坡渗流场模拟各层岩土体参数选取岩性天然重度(kN/m3)饱和重度(kN/m3)饱和体积含水量(%)饱和渗透系数（10-7m/s）抗剪强度指标天然饱和C(KPa)()C(KPa)()滑坡堆积层19.720.522.86.938333227滑 带19.720.523.26.922.52219.120.7滑床20.421.422.40.860425236操作步骤：1. 输入模型2.输入材料属性3.输入材料属性（因为模型是直接从seep/w里面导入的，没有取消seep里面使用的划分网格）4.画地下水位5.draw slip surface and exit range6.查看区域属性7. 画孔隙水压力云图8. 查看数据9. Solve可以看到安全系数为1.39410.极限平衡法对应的每个条块受力图11. 查看不同的方法对应的不同安全系数12. 画孔隙水压力图13. 输出滑坡稳定性报告FLAC3D数值模拟上机题计算模型分别如图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剪切弹性模量：体积弹性模量：1.坡角为30时图1 倾角为30的边坡(单位：m)输入源代码：set log onset logfile stability30.loggen zon brick &p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10gen zon 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 10model mohrPlotShowAdd surface blueAdd axes yellowprop dens 2500 bulk 9.26e7 shear 3.79e7 coh 4.2e4 ten 8e5 friction 17 dilation 20fix x y z range z -0.1 0.1fix x range x -0.1 0.1 fix x range x 99.9 100.1fix yset gravity 0 0 -10def calfosminf=0maxf=2loop while maxf-minf0.01 fs=(maxf+minf)/2.0 refric=atan(0.31/fs)*180/3.14 recoh=42000/fs commandini sxx 0.0 syy 0.0 szz 0.0 sxy 0.0 sxz 0.0 szz 0.0ini xvel 0.0 yvel 0.0 zvel 0.0ini xdis 0.0 ydis 0.0 zdis 0.0pro fric refric coh recoh set mech ratio 1e-5solve step 5000print fs end_commandaa=mech_ratioif aa1e-5 thenminf=fselsemaxf=fsend_ifend_loopendcalfossave 2.sav输入指令显示图形可以得到fs=1.476，边坡是稳定的。2. 坡角为45时修改对应的尺寸即可gen zon brick &p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0