Geostudio及FLAC3D

1、高程（同1450, " ? »- 七 .1350 二力二：；：；%_.,壁.；* ,； ；:一 Vk. ZK3+'二；“、4旌产4T1 - I T"W4x L *- -r_噫犍赤k里端蜀学a.4*军宝灌冢.1350130c1250Ejl四系全新统；四系全新统|坡堆积层S第四系冲洪 积层推测潜在滑面ZK41286.01. 13007.8(12582。西河.12501200碎块石土I层界线目剖面方向1200地质工程数值法GeoStudio上机实习作业内 容：用GeoStudio的Seep/w模块模拟滑坡不同工况下的地下水渗 流曲线。要求：提供模拟结果文件，现场

2、验收操作相关资料：?滑坡工程地质剖面图;?岩土参数;?渗流模拟工况1计算模型的建立根据2#滑坡的工程地质条件，选择该滑坡钻探剖面（图 1）作为模拟剖面。其方位 NW113 ,滑坡前缘高程1190m,后缘高程1398m,前缘和后部较薄，中部较厚，根据钻孔 揭露，滑体物质主要是灰黑色碎块石土，滑带土主要为灰黑色含角砾粉质粘土、角砾质土或含角砾泥炭质土，滑带厚度0.61.5m,滑床为炭质页岩和炭质灰岩的碎块石土。基于以上条件，对该堆积层滑坡的坡体地质模型进行网格剖分，共剖分为2829个网格单元，如图 2所示。1450高程) 113 一比例尺1:1000.14001150050100150200250

3、 300350400450 500 550600650-1150滑动面及潸在滑动面钻孔编号滑动方向及滑动方向I孔口高程图1 2#滑坡工程地质剖面图操作步骤：1.打开 Geostudio ：2.操作界面尺寸设置:欢迎下载Grid fpaang 任rg LnX: 1I： 1Eng. Uhits; HnetersnrrmX； 0,077637Y； 0.Q773S7J Esplay GridZ 访蜃区网怎Dose3.导入图形欢迎下载4.检查图形5.设置材料属性欢迎下载Kcn Hydraulic Ccndurtivrt Funrti&nsaddDrlfrUPCCfNanwr11TpprCur/e

4、 Fit:即 pTDTFxart必SEIbts：13 44Z2Z?IS 怎找网c 3u店Lio口 QJil咕E Scpe of Functor0 gdit Data PointeAdd Parasklc必以由世J _oe Aug £zkK Axis Range;to4naly&e.Remiens., L'neu.PqinUIMateridlsn.Boundary Conditions.Hydraulic FunctionsAir FunctionsSpatial FunctionsInKial Water Table .Flux Sections.Generate P

5、lan View.6.将材料属性赋予给图形区域欢迎下载7.划分网格欢迎下载8.设置边界条件，修改颜色9.赋予边界条件欢迎下载.dny如 f pE-ktii 眄区"10.检查模型是否正确国 Venf/ZOp-timize DataInfcir 四占 ticin:'.nh'irg 7】三t 三a二门 hoce 区 a：Uzhet to an & arrer：.Verifyiig infinite dements.Verifying secondary nodes for sach element.,.Verifying that eadn element has

6、a positwe area.,.Verifyinfl that not 司1 elements are NULL dements.,.Verifying ma terigl props比sVerify infl boundary 8nditions Verifyino flux sections, 一Cheddno settings.Verifying integraticn ElesVerificatiori completed - error（sl。丸,日mingG）,11.保存，并solve欢迎下载欢迎下载12.查看结果与预期结果相一致，上层碎石土的渗流力明显大于滑床的。2计算参数与工

7、况参数的取值主要由勘查资料中的岩土体物理力学试验、水文地质试验、工程地质类比以及经验值等综合确定。具体为：模型材料的天然、 饱和容重等由渗透试验确定；饱和体积含水量和饱和渗透系数，类比附近滑坡或者同一类型滑坡的物质组成及坡体结构等，确定出参数取值范围，然后根据野外钻探的实际资料，通过参数调整，反演出2#滑坡的实际地下水位，获取后续模拟的参数，渗流场模拟过程中所涉及到的各层岩土体参数见表1。采用GEO-Studio软件中的SEEP/W模块进行滑坡渗流场模拟时， 视滑坡岩土体材料渗透 系数和体积含水量为坡体孔隙水压力的函数。当岩土体处于饱和状态时，渗透系数和体积含水量为定值；而当岩土体处于非饱和状

8、态时其渗透系数和体积含水量应该通过试验来确定。但鉴于试验的难度，一般取经验值。本节进行2#滑坡的渗流场*II拟时采用SEEP/W模块自带的Van Genuchten经验曲线和饱和状态时的参数来确定当坡体材料处于非饱和状态时的渗透 系数和体积含水量与孔隙水压力的函数关系，分别如图3图8所示。渗流模拟工况：一、自重 +地下水；二、自重+50年一遇15日连续强降雨工况（50年一 遇降雨强度为24小时最大值130.0mm）。表1 2#滑坡渗流场模拟各层岩土体参数选取饱和饱和体积含水量（）饱和渗透岩性重度(kN/m3)重度(kN/m 3)系数(x 17m/s)天然C(KPa)d()饱和C(KPa)M)滑

9、坡堆积层19.720.522.86.938333227滑带19.720.523.26.922.52219.120.7滑床20.421.422.40.86042523610欢迎下载操作步骤:1.输入模型2.输入材料属性欢迎下载113.输入材料属性（因为模型是直接从seep/w里面导入的，没有取消seep里面使用的划分网格）欢迎下载12SIS eprE/啊>1QLriL I'B E*W #.-MArl im|h| 中 lh* .hww£c lbw施 ZPJSMM 1>*5.draw slip surface and exit range£rf+ iw

10、63;»« vyln >hf Shck .u点S ，6卜 Hj'ipwR" * I口占昌吟 -<BIBI M 可Driw Up uM-heh Ef*t ard 金at 虹prMrii 。访jih Pone>Ei >口£4<(-KXtinT*费15即*M1= 毋# gfLlT*E白做，yi 口即一*Q，* M/审_"尸曲题1 0 111 ”片jgw- I6.查看区域属性欢迎下载13把:工遥View Object Information,bility- ower：-CDkjr-Model- Uni: Weigh

11、t-CQhestjnPtii-Phi-EF4>hr-Cojlcnb 21.4kN/in3 5； kPa 武> 00由Pg V.ater Pre* Peiris As1旧畤JrrP7.画孔隙水压力云图IIl 1 % 之，二三/广义的球El3中 Km Slnis T" rikiE q*.L(W>t *f r c b r , ：r v r 1.即£/ S44tk4»» fSIMWBPqNERif E+f电电*win LVgaa卜 SPT匕 4 PCu * u占小、七<< H)Bi 册-FOkt:EMrn£«r*

12、u8.查看数据欢迎下载14If二一isFf7-sr-at下MStop159.Solvelb chengheng.gz： - SLOPE/W SOL», o_ =File Help回 SLOPE/WOO7二口口如 File; dnenqshsng.gsrAnalysis: Slope Staisili+yMinimum Factor of Safetyr*tonent ForceOrdinary:1.367=Bishop;L395-Jaritxi:-1366M -P;L3941.394Slip StJrfece 新125 of 125Searching foCritical Slip

13、SurfaceGte叶欢迎下载*46C.&S232.可以看到安全系数为1.394Slip Surfaces欢迎下载1610.极限平衡法对应的每个条块受力图国 Free Body Diagram & Force PolygonFactor of S/ety1 35424Phi Angle36nC (Strength)52 kPaEC iTorce)12B13kNPg< Water-9B1 取 kPa1-=i|Pore Water Force-241.92 kNPore Jr PressLweOkPaPore Ar ForceOkNPhi g Angle0sSlice Wid

14、th017B3CmMi 小 Height0.056527 mBase Length口 2修9 mRa4A 4run|pdi 5 nQapy DiagramCloseTiffie. f11.查看不同的方法对应的不同安全系数欢迎下载17林应一 M iifiiji th , 唐 0 二 C3 P | K 3 学杂洌击，*| 0 &E4的口山1度/加湎时1.394*1.395 12.画孔隙水压力图欢迎下载18 Draw GraphDeleteNamepori 占ter pressureName:pore-'Aaber pressureata fron:SJioe DataPwe'

15、; Wster Pressure| Su rm (Y) vc. Avpfflw 0 VaIlmsUndo .PrefeercejaoEty_|%EiibchMateralru怔FaintsFolnftLpkj/jMfltpid 田undarU-ground Surface UmXOPEW OoticraSlBW.P.Cmditnm_rc 卜二5a:"nd. =J| Acndhg Arrani甘耻 Cdgr; | |s£t ” MashAnaly«/ ter tcrxrfrent lLaad3 SorrLMds d| SurdiagB LoadsSurdiarge

16、Lfiad QbadrgFT ISFont”. A'id Regular办自Point L.abe&!O*St?m. LJdds:BReqia力 Labels:12#xis LabiBEEIQNoda Ltb«t:8-conwjr nabtte-ioEM 5:始欢迎下载叵 Auto Refresh叵 Use Auto Sc出RefreshX Axis Range:toCloseSflitCTZ 的口 OtjfrM/ Ak>V Aclure二小比MTdh/ CcmtMr Lm»回 T«m 0ridtLin«*1 TErskin Gac

17、k Shadingqp Wrteec回 VWfl£of5Kid sip surges Alrrf5-Jr5Fof £ S*gnant Figurei! 1-回mw 也jtp看 41Pbirfaces： 国 ;团Sip9urfeoeShading: |.J团'totfSlras im&ThdaiME： jihrn Jxs191、工玄 _<1»-凰Hl囹hi同 丁 y归事.1了 F E中产产二- H早（ d.F4 JI M |»IM T| VirtJE , 113.输出滑坡稳定性报告?i >"* igsrinvua Pf

18、r, *>«" f 坤 “a 力-m 审 Lif > 孙| 申瑜II*Lidir InEn'man*ra）MS Eritrw|i-畤1 k Eiling.SIom- Sl¥b<f>% JEHfr ,串 JEih * LMipC-iirxEnrlrK ,Lnaa Mjoflwt g k 配片* -Ir* q/ hid- SlmcT itoSudK.Ale- IniformationkAnrtKH Nmlm J>ftri inr. "L"叫中箕m上m但用5（也11小& l ”EnckoM 一 I. ,

19、. ij. ；8-4修* ailLlI/l%>1Pr川get SMllngi， 的lELlgU： irlr w EHrjucXrM iiM-s*1 rdffMiFiurri&'h taaiM4p. 4M= i:a Hntifth Uwts： e* 收11内峰中” #r wi41： *.门） 'I VifiniL20欢迎下载flaCd数值模拟上机题计算模型分别如图1、2、3所示，边坡倾角分别为30。、45。、60。，岩土体参数为:密度p= 2500 kg/m3,弹性模量E= 1X108 Pa,泊松比科=0.3,抗拉强度b t=0.8X106 Pa,内聚力C=4.2X

20、104 Pa,摩擦角e = 17° ,膨胀角A= 20°试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。附换算公式：1 kN/m 3= 100 kg/m 3剪切弹性模量：G E一2(1)体积弹性模量：K一E一3(1 2 )输入源代码：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 b

21、rick &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 10 model mohr21欢迎下载ShowAdd 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

22、.1 0.1fix x range x 99.9 100.1fix yset gravity 0 0 -10def calfosminf=0maxf=2loop while maxf-minf>0.01fs=(maxf+minf)/ 2.0refric=atan(0.31/fs)*180/ 3.14recoh=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.0 pro fric refric coh recoh set mech ratio 1e-5solve step 5000 print fsend_commandaa=mech_ratioif aa<1e-5 then minf=fselsemaxf=fsend_ifend_loopendcalfossave 2.sav22欢迎下载甘 H-ACJD 匕臼事 DKpLty Option e. Plot Windw Ht* CcHTTurund WT nd ow ft Jllijr >11 I. vr mup