UDEC流固耦合-边坡渗流实例解析

合集下载

UDEC滑坡实例步骤

1、加载UDEC进入DOC环境后输入giic或者gui命令，然后进入主菜单2、Model option 选择合适条件通常情况下，你可以使用默认域联系（domain-logic）检测模式。

如果你想监测任何块体的位移，这些块体可能从隧道顶部分离或掉落，你应该使用“cell-space detection”模式跟踪位移和下落块体的潜在接触。

3、命名并且保存文件4、New block 建模，根据需要设置模型的长30 宽15415 410 405 400 395 390 385 380 375 370415 410 405 400 395 390 385 380 375 370 420 425 430 435 440 445 450 455 m5、Bound 调节边界，与实际相符6、Crack 添加节理，（层状岩体是否按节理处理？）岩层20°∠34°，J1产状60°∠15°J2产状为35°∠47°，J3产状为95°∠89°（怎么将不同产状节理进行转换？）路线设计好，为后来开挖做好准备。

7、execute 执行文件8、zone 执行长度为0.5的最大区域边界，划分网格9、Zone material 创建一个或者几个块体材料属性，选择一种本构模型，本次选择的是Mohr-Coulomb模型prop mat 1 den--2143 bu=30e9 sh=18e9 c=1.2e5f=21 t--2e5prop mat 2 den=2260 bu=40e9 sh=24e9 c=1.5e5f=28 t--2.5e5prop mat 3 den--2300 bu=50e9 sh=28e9 c=3.5e5f=32 t--3.5e510、JointMat 节理属性，选择;liexi 1prop jmat 1 jkn=4e9 jks=2.Oe9 j卜18 jc=0 jt=0jieli2prop jmat 2 jkn=4e9 jks=2.Oe9 jf=18 jc=le4 jt=0;huadai;bianxingtidijiemian3prop jmat 3 jkn=4e9 jks=2e9 jf--10jc=2e3 jt=0;duanceng4prop jmat 4 jkn=2e9 jks=0.7e9 jf=13 jc=le3 jt=011、instiu -bound 边界条件左右边xvelocity为0，下边界yvelocity为0保存slope112、utility 访问不同变量，可以监测将要下滑的块体垂直和水平位移13、setting 设置重力9.8114、run -solve 自动检测平衡，在分析的不同阶段保存项目为slope2、3、415、build -cut or fill 如果删除某一块后则可立即保存，评价岩体的稳定性用solve工具实现，计算稳定状态得到开挖后岩体周边的位移通过y位移等值线图来说明，点击“Contour-Motion/ydisplace”输出，保存slope516、接下来可以评价有弱节理的岩体中因开挖引起节理岩体的反应，首先回到slope1，A分支包含关于强节理岩体开挖的保存文件，B分支是一个新的分支，你可以用来执行弱节理的模拟。

基于流-固耦合理论的水体下盾构隧道施工渗流规律

基于流-固耦合理论的水体下盾构隧道施工渗流规律随着水下隧道工程建设技术日趋完善，水下盾构隧道已经成为跨越水域的重要交通工程之一。

在隧道工程施工过程中，由于存在着稳定的补给水源，渗流场和应力场相互作用、相互联系，使得隧道围岩孔隙压力发生变化，隧道围岩极易发生失稳，造成严重的工程事故。

近年来，国内外许多学者对隧道施工过程中流-固耦合问题进行了研究。

YI et al［1］通过数值模拟研究分析了盾构隧道开挖后孔隙水压力和地层沉降的分布规律；GALLI et al［2］对二维模型和三维模型计算结果进行了对比，认为三维模型更加符合工程实际；纪佑军等［3］采用数值模拟，认为注浆可以有效控制隧道围岩变形和地下水涌入；原华等［4］通过对上海越江隧道东线中段的数值模拟分析，得出管片在渗流作用下的受力更加趋于“平均”，地下水有利于隧道衬砌管片的受力；吉小明等［5］通过分析渗流作用下隧道开挖过程发现，隧道周边最大剪应力和位移都有所增大；马龙祥等［6］假设当软土盾构隧道达到稳定渗流状态时，作用于衬砌外壁的总水头为一个未知常量，计算出稳定渗流时隧道的闭合解析解；易小明等［7］通过对比分析流-固耦合作用下海底隧道开挖数值模拟的结果与现场监测数据，认为两者具有一致性；贾瑞华等［8］通过数值模拟分析得出了隧道开挖面处的渗流场、位移场以及孔隙水压力的分布规律。

然而，国内外对高水头压力作用时盾构隧道施工扰动下地表的沉降量、围岩渗透性、管片内力分布等方面的研究并不多。

近几年来，穿越江、河的盾构隧道工程越来越多。

因此，研究高水头压力作用下盾构隧道的地表沉降、孔隙水压力分布、施工扰动下围岩的渗透性以及隧道衬砌管片的受力特征，对于盾构隧道施工及运营期间的安全稳定具有重要意义。

本文选取太原地铁2号线双塔西街—大南门站下穿迎泽湖区段为工程背景，采用有限差分软件FLAC3D内嵌的FISH语言二次开发建立三维数值模型，得出不同水头压力作用下的盾构隧道在流-固耦合效应下的渗流场、位移场以及管片衬砌的孔隙水压力的分布规律，以期为水下隧道施工工程风险控制提供理论依据和技术支持。

渗流-应力耦合作用下边坡支护稳定性分析

渗流-应力耦合作用下边坡支护稳定性分析渗流-应力耦合作用下边坡支护稳定性分析在岩土工程领域中，边坡稳定性一直是一个重要的研究方向。

边坡工程中的稳定性问题往往涉及到多种因素，如地下水渗流、土体的力学性质等。

其中，渗流-应力耦合作用是边坡稳定性分析中的重要考虑因素之一。

本文将重点分析渗流-应力耦合作用对边坡支护稳定性的影响并进行详细讨论。

首先，我们需要了解渗流-应力耦合作用的基本概念。

渗流-应力耦合作用是指地下水渗流与土体应力状态相互影响的过程。

在边坡工程中，当地下水渗流进入土体中时，水力压力会改变土体的力学性质，从而对边坡的稳定性产生影响。

同时，土体的应力状态也会影响地下水的渗流路径和速度。

因此，边坡稳定性的分析必须考虑渗流-应力耦合作用。

接下来，我们将详细分析渗流-应力耦合作用对边坡支护稳定性的影响。

首先，渗流-应力耦合作用会导致边坡内的水力压力变化。

当地下水渗流进入边坡内部时，水位的上升会增加土体中的水力压力，从而增加边坡的自重。

这会使边坡受到更大的重力作用，从而增加了边坡的倾覆风险。

此外，水力压力的增加还可能导致边坡土体的饱和度增加，从而引起土体的液化现象，使边坡更容易失稳。

其次，渗流-应力耦合作用还会改变边坡的有效应力分布。

当地下水渗流进入边坡内部时，水力压力的增加会降低土体的有效应力，在边坡支护结构处产生更大的水平承载力和竖向支撑力。

这会增加了边坡支护结构的稳定性，防止边坡的下滑和滑动。

然而，如果渗流-应力耦合作用导致边坡土体饱和度增加，土体的强度将大大降低，对支护结构的稳定性构成威胁。

最后，我们需要采取合适的工程措施来应对渗流-应力耦合作用对边坡支护稳定性的影响。

一方面，可以通过合理设计排水系统，控制地下水的渗流路径和速度，减小渗流对边坡稳定性的影响。

另一方面，可以选择适当的支护结构，提高边坡的抗滑能力和抗液化能力，确保边坡在渗流-应力耦合作用下的稳定性。

综上所述，渗流-应力耦合作用是影响边坡支护稳定性的重要因素。

(完整word版)UDEC模拟实例与解析

UDEC 实例翻译与命令解析中铁隧道集团科研所——珠穆浪玛UDEC 实例翻译与命令解析翻译：珠穆朗玛1 地震诱发地层坍塌 Seismic-Induced Groundfall1.1 问题描述本例展示使用 UDEC 模拟分析地震诱发地层坍塌的一类的问题,模型见图 1.1,该模型基于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿 34-1-554 切割断面的一个剖面图的结构和尺寸. 用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高 5m,宽 10m.假定两个连续节理交叉平面分析:一个角度为 45 度,另一个为-9 度,两者节理间距均为 5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳定性。

围岩参数来自试验室平均测试数值，假定岩石块体参数如下：假定块体仅具有弹性行为，节理假定符合库伦滑动准则，选择典型的教课书数值作为节理参数，如下：初始应力状态按各向同性估计为24Mpa（假定垂直荷载由覆盖深度大约800m 的岩层产生）。

1.2 UDEC 分析UDEC 模拟顺序分三个阶段,首先,模型在初始应力状态下进行无超采固结.其次,进行开挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图1.2.超采正上方和下方的块体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题域的外周边界引入粘滞边界用以消除波的反射.从而模拟有限的岩体,地震事件用施加到模型顶部y 方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中,施加1.25Mpa 的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部, 施加的有效影响应力应该是1.25 MPa/2, or 0.625 MPa.0.02 秒后的开挖面拱顶的应力分布见图1.3,两点的位移被监测,1 点位于开挖面的左角,点2 位于拱顶块体的右角, 图1.4 的位移时间曲线显示两点本质上是弹性反应.本例关心的问题是在模型顶部施加的速度和计算速度的对比,下面的公式可以用以估计施加的波速.使用这个方程,施加的最大波速大概是0.04m/sec,图1.5 显示的峰值波速小于0.06m/sec. 估计的波速和监测波速的不同在于使用的围岩模量.而是没有考虑节理变形的相等变形模量.在第二个案例中,施加应力波峰值12.5 Mpa(有效应力6.25Mpa).0.02 秒后的开挖拱顶应力分布见图1.6.该图显示出拱顶岩体不受力,表面该块体已经松散并正在下落.对于关心的问题,后来三个时间的几何体和应力分布见图1.8 至图1.10.在问题的顶部预测的波速(从上面的方程)是0.4m/sec.从模型中计算的波速见图1.11,再次,由于使用的是原岩弹性模量而不是岩体的变形模量导致预测和监测的波速之间的差异.1.3 节包含了该模型的数据列表,该列表包含了一个FISH 函数(show)被用来创建坍塌的动画文件,每隔0.02 秒俘获一个显示的图片.通过改变FISH 参数time_int 可以改变动画帧的间隔.视图的总数也可以通过改变snap_shot 的数值进行改变.为了显示80 帧的显示图片而创建的该电影文件需要大概13MB 的硬盘空间.1.3 数据文件列表Example 1.1 SEISMIC.DATtitleSEISMIC INDUCED ROOF COLLAPSE 地震诱发拱顶坍塌;round 0.01; define original boundary of modeled region 定义模型区域的原始边界block -25,-20 -25,20 25,20 25,-20; generate joint pattern over entire original region 在整个原始区域生成节理形态jregion id 1 -25,-25 -25,25 25,25 25,-25jset 45,0 200,0 0,0 5.0,0 (0,0) range jreg 1jset -9,0 200,0 0,0 5.0,0 (0,0) range jreg 1; put in joints needed for the later excavation 为了后面开挖而设置的节理crack -5.01,-2.51 5.01,-2.51crack -5.01, 2.51 5.01, 2.51crack -5,-2.5 -5,2.5crack 5,-2.5 5,2.5crack 2.25,2.5 1.93,5.0; generate fdef zones and assign joint properties (mat=1 & jmat=1;default) 生成单元和设置节理参数generate edge 9.0 range -30,30 -30,30prop mat=1 d=0.00300 k=39060 g=31780prop jmat=1 jkn=20000 jks=20000prop jmat=1 jf=30.0; apply boundary conditions and initial conditions to 在地应力下施加边界条件和初始条件; consolidate model under field stressesbound stress=-24.0, 0.0, -24.0 ygrad=-0.3 0 -0.3insitu stress=-24.0, 0.0, -24.0 ygrad=-0.3 0 -0.3bound yvel 0.0 range -26,26 -21,-19grav 0.0 -10.0; track the x-displacement, and y-displacement over time 追踪位移hist solvehist xdis=0,7 ydis=0,7 type 1solve rat 1e-5; save consolidated statesave seismic1.sav; make excavationdelete range -5,5 -2.5,2.5solve rat 1e-5; save excavated statesave seismic2.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.04 m/sec);; set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -1.25 yhist=cos(100.0,0.0195) range -26 26 19 21;reset time hist disp rothist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 mass; 0.02 sec.cyc time 0.02save seismic3.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.4 m/sec); set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -12.5 yhist=cos(100.0,0.0195) range -26 26 19 21reset time hist disphist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 masssave seismov.sav;; 0.02 sec. —————————————————————————————————————UDEC 实例翻译与命令解析中铁隧道集团科研所——珠穆浪玛cyc time 0.02save seismic4.sav; 0.25 sec.cyc time 0.23save seismic5.sav; 0.50 sec.cyc time 0.25save seismic6.sav; 0.75 seccyc time 0.25save seismic7.sav;rest seismov.sav; make a movie of the groundfall;wind -12 12 -12 12set ovtol 0.05plot block vel max 2.0 blue stress max 50movie onmovie file = seismic.dcxmovie step 1000step 400003 隧道支护荷载Tunnel Support Loading3.1 问题陈述本例模拟展示了UDEC 在检查衬砌隧道方面的应用,着重强调了荷载在混凝土衬砌中的发展,本例也解释了模拟连续建造操作中独立阶段的模拟程序.隧道系统的理想几何体见图3.1.系统包含在海床下大约70m(中线)深度,中线间距12m 的两个隧道, 初始水位在隧道中线上方110m 处.服务隧道直径5.24m,衬砌厚度37cm.主隧道直径8.22m,衬砌厚度46cm.服务隧道先于主隧道开挖和衬砌.随后设置主隧道衬砌,水位上升增加到100m.—————————————————————————————————————UDEC 实例翻译与命令解析中铁隧道集团科研所——珠穆浪玛施工顺序是:(1)开挖服务隧道excavation of the service tunnel;(2)衬砌服务隧道lining of the service tunnel; (3)开挖主隧道excavation of the main tunnel; (4)衬砌主隧道lining of the main tunnel; and (5)升高水位raising of the water level.分析的目的是评价每个施工阶段服务隧道和主隧道支护状况.本例的材料参数见下:岩体——开挖隧道的围岩参数为：弹性模量elastic modulus 0.89 GPa泊松比Poisson’s ratio 0.35单轴抗压强度uniaxial compressive strength 3.5 MPa粘聚力cohesion 1 MPa密度density 1340 kg/m3混凝土衬砌——弹性模量为24 GPa ，泊松比为0.19. 假定衬砌为线弹性材料。

边坡坡面流-渗流耦合冲蚀的模型实验

ｐｅｒｉｍｅｎｔａｌａｐｐａｒａｔｕｓｎｄａｍｅｔｈｏｄｓｆｏｒｅｘｐｌｏｒｉｇｎｔｈｅｓｃｏｕｒｉｇｎｐｒｏｃｅｓｓｏｆｓｌｏｐｅｓｏｉｌｄｅａｃｔｈｍｅｎｔ（ｓｈｅｅｔｅｒｏｓｉｏｎ）ｂｙｓｈａｌｌｏｗｆｌｏｗｓ，ａｔｔｅｍｐｔｉｎｇｔｏｒｅｖｅａｌｔｈｅｃｏｕｐｌｉｇｎｍｅｃｈａｎｉｓｍｏｆｓｅｅｐａｇｅｏｎｓｌｏｐｅｅｒｏｓｉｏｎ，ａｎｄｔｏｐｒｅｄｉｃｔｔｏｏｇｐｒａｐｉｃｈａｌｃｈａｎｇｅｓｏｆｓｌｏｐｅｍｏｒｅａｃｃｕｒａｔｅｌｙ，ｗｈｉｃｈｃａｎｂｅ
法的可行性。关键词：岩土工程；坡面侵蚀；渗流；坡面流；实验；ＡＢＡＱｕｓ中图分类号：ＴＶ６９７文献标志码：Ａ文章编号：２０９５ —２７８３（２０１３）０５ —０３６３— ０４
Ｅｘｐｅｒｉｍｅｎｔａｌｒｅｓｅａｒｃｈｏｎｓｌｏｐｅｓｏｉｌｄｅｔａｃｈｍｅｎｔｂｙｓｈａｌｌｏｗｌｏｆｗｓａｎｄｕｐｗａｒｄｓｅｅｐａｇｅｓ
ｉｎｔｒｏｄｕｃｄｅｉｎｔｏｔｈｅｔｈｅｏｉｒｓｅｏｆｓｌｏｅｐｅｒｏｓｉｏｎｉｎｒｅｌａｔｉｏｎｔＯｓｅｅｐａｇｅｅｆｆｅｃｔｎｄａｉｎｔｏｄｓｉｅｇｎｏｆｐｒｏｔｅｃｔｉｏｎａｇａｉｎｓｔｅｏｓｒｉｏ．Ｂａｎｓｄｅｏｎｔｈｅｎｕｍｅｒ —

水中填筑围堰边坡稳定的流—固耦合分析

水中填筑围堰边坡稳定的流—固耦合分析通过对流-固耦合的特征和渗透边界的分析，研究了在水中填筑围堰边坡时会受到的影响和应该注意的问题，重点研究了在使用流-固耦合时的特征时间和边界确定，这两个方面的精确度直接影响了最终的计算结果。

标签：围堰；流-固耦合分析；边坡稳定渗流是广泛存在于各大水利工程中的常见问题，渗漏会造成很多严重的后果，甚至是威胁到整个堤坝的安全和稳定，现在大多数的研究都将目光集中在被水浸泡的部位，这些部位在长期的浸泡中变成软土，不具有强大的支撑力。

这些研究虽然取得了不少的成果，但是在实际应用中仍旧没有很好的效果，主要原因是是因为忽视了水渗透的特性和边界，导致整个堤坝的稳定性都受到了影响。

一、水对边坡的作用想要真正的处理好渗透问题就必须要正确的认识到水对边坡的影响和作用。

水是水利工程中最主要的主体，它的特性影响了整个工程的建设，水作为一种流动的液体，具有很多固体所没有的特点，也对堤坝岩土主体造成了不小的影响，长时间的积累会让岩土体的工作效果大大折扣。

其中，水对边坡的作用主要表现在以下几点上：（一）减少岩石的有效重量岩石的有效重量是指岩石在建设中的受到地心引力对水产生的压力作用，起到了稳固堤坝的作用。

而水作为一种流动液体，可以变成细小的分子渗透在岩石的各个孔隙之中，在渗透的过程中会产生浮力效应，将岩石受到的重力减小，这样就让边坡整体的抗滑力下降，造成整个边坡失稳。

（二）降低岩石的强度水中含有很多的化学物质和矿物质，这些物质在水的渗透过程中会和岩石内部的物质发生反应，改变岩石的内部结构，降低岩石整体上的强度。

岩土主体的强度下降会对整体的稳定性造成一定程度的影响。

（三）对岩石具有较大的冲击力水在运动的过程中会对岩石产生较大的冲击力，日积月累之后会造成岩石的变形或蚀变，造成岩石的不稳定，当遇到比较的降雨或者冲刷时会在巨大压力的作用下发生崩塌，造成大型的滑坡事故。

二、流-固耦合分析研究随着科学技术和计算机的发展，对边坡稳定性的计算也越来越精确和实用，尤其是在上世界70年代出现的有限元法等理论，给各种复杂渗透问题的计算提供了便利。

基于流固耦合的边坡稳定性研究

基于流固耦合的边坡稳定性研究基于流固耦合的边坡稳定性研究摘要：随着工程建设的不断发展，边坡工程在土木工程中的重要性日益凸显。

然而，由于气候变化和自然灾害等因素的影响，边坡的稳定性成为一个严重的问题。

为了有效地评估和预测边坡的稳定性，研究者们采用了基于流固耦合的方法，通过模拟边坡受到降雨等外力作用时发生的水流和土壤运动相互耦合的过程，从而更真实地模拟边坡的稳定性。

关键词：边坡；流固耦合；稳定性；模拟1.引言边坡是指山体或土体斜坡面与地面之间的交界处，其稳定性是土木工程中的一个关键问题。

然而，在实际工程中，边坡往往受到多种外力的作用，如降雨引起的水流，地震引起的地表振动等。

这些外力，尤其是降雨引起的水流对边坡的稳定性产生了重要影响。

因此，研究者们开始采用流固耦合的方法来研究边坡的稳定性。

2.基于流固耦合的边坡模拟流固耦合是指边坡受到外力作用时，水流与土壤的相互作用和相互影响的过程。

基于此原理，流固耦合模拟将水流和土壤运动联系在一起，以更真实地模拟边坡的稳定性。

2.1 水流模拟在流固耦合的边坡模拟中，首先需要对水流进行模拟。

水流模拟一般采用流体力学的数值计算方法，如有限元法或有限差分法。

通过数值模拟，可以计算边坡上任意位置的水流速度、水压和水位等参数。

2.2 土壤运动模拟在水流模拟完成后，需要模拟土壤的运动过程。

土壤运动模拟一般采用岩土力学中的塑性本构模型，如Mohr-Coulomb准则。

通过建立合适的数学模型和边界条件，可以计算边坡上土壤的运动过程，包括土体的位移、形变等。

2.3 流固耦合在完成水流模拟和土壤运动模拟后，需要将两者耦合在一起。

通过将土壤运动过程中产生的位移和形变输送到水流模拟中，同时将水流对土壤的作用考虑到土壤运动模拟中，使得水流与土壤之间相互影响和相互耦合。

3.边坡稳定性评估与预测基于流固耦合的边坡模拟可以有效地评估和预测边坡的稳定性。

通过模拟边坡在不同外力作用下的变形和破坏过程，可以得到边坡的位移、形变和应力等参数。

隧道工程中的流固耦合问题研究

隧道工程中的流固耦合问题研究一、引言地下隧道建设是人类工程史上的伟大壮举，既改善了城市交通环境，又方便了人们的出行。

然而，隧道围岩的流固耦合问题一直是影响隧道施工和使用的关键因素。

本文将探讨隧道工程中的流固耦合问题，并提出相应的解决方法。

二、隧道工程中的流固耦合问题隧道建设中的流固耦合问题指隧道中负责流动的液体与固体隧道围岩之间相互作用的问题。

1. 套筒压力在隧道施工的过程中，为加固隧道围岩结构，常常会使用套筒增加隧道的支撑能力。

但是，在套筒与岩体之间的空间中，常常会产生水压。

当水压达到一定程度时，套筒就会发生变形，使得隧道支撑能力下降。

2. 稳定性问题在一些地形较为崎岖的地区，隧道的稳定性会受到地形的影响。

在这种情况下，隧道施工过程中，液体的流动会对隧道的安全性产生影响，隧道围岩的稳定性也会受到相应的影响。

3. 水压问题在某些条件下，隧道中液体的流动会导致水压的问题。

水压会导致隧道中的固体结构发生移位，对隧道的结构稳定性产生影响。

三、解决方法隧道工程中的流固耦合问题在一定程度上会影响隧道建设的进度和质量。

下面列举一些解决方法：1. 施工中的控制在隧道施工的过程中，水压的控制是保证隧道施工进度和质量的关键因素之一。

当施工人员发现水压过高时，应及时采取措施以控制水压。

2. 模拟实验通过对隧道施工中不同情况下流体在固体结构中的流动状况进行模拟实验，可以找到更加准确的解决方法。

模拟实验过程中，可以对不同变量进行控制，从而找到最合适的施工方法。

3. 选择合适的隧道设计方案在进行隧道设计时，考虑流固耦合问题对隧道施工和使用的影响，选择合适的隧道设计方案。

通过优化隧道设计方案，可以减轻隧道中的水压问题，提高隧道结构稳定性。

四、相关技术和设备在解决隧道建设中的流固耦合问题过程中，相关技术和设备的运用是不可或缺的。

下面介绍一些主要的技术和设备：1. 隧道测量技术隧道测量技术是对隧道施工过程中液体和固体结构的流动状况进行监测和记录的技术。

基于UDEC的渗流诱发岩质边坡失稳分析

．
ＹｕＹｅｑｉｎｇＨｕＪｉｎｇｑｉａｎｇＭｕＣｈａｎｇｐｉｎｇ
（ＳｉｃｈｕａｎＭｅｔａｌｌｒｕｇｉｃａｌＤｅｓｉｇｎ＆ＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ）
余业清胡敬强母昌平
（四川省冶金设计研究院）
摘
要以某矿为工程背景，采用理论分析和数值模拟计算相结合的方法，对受渗流载荷作用
下的露天边坡稳定性进行了研究。结果表明：液体沿节理面以一定的水力坡度发生渗流并产生孔
ｂｉｌｉｔｙｗｉｌｌｈａｐｐｅｎ．
ＫｅｙｗｏｒｄｓＲｏｃｋｓｌｏｐｅ，Ｓｅｅｐａｇｅ，ＵＤＥＣ，Ｓｔａｂｉｌｉｔｙ
自然条件下，岩体经过长期地质构造作用，内部产生复杂、多样的结构面。岩体的结构特性决定了岩体的力学性质，也是引起岩质边坡失稳破坏的重要因素。影响边坡稳定性的因素可分为内因和外因，内因主要是岩体的岩性和结构，外因包括地下水、人工开挖及地震作用等。边坡的稳定性分析研究始于２０世纪２０年代，其分析方法可概括为确定性方法和非确定性方法两类］确定性方法以极限平衡分析法和数值计算为代表；非确定陛方法包括分析几何法、模糊数学法、非线性理论和人工智能法等。目前工程实践中主要应用确定性方法。随着计算机技术的快速发展，离散单元法引、有限单元法、有限差分法及边界单元法等数值计算方法在边坡稳定性分析中的运用日趋广泛。

《2024年流固耦合渗流规律研究》范文

《流固耦合渗流规律研究》篇一一、引言流固耦合渗流规律研究是近年来多学科交叉领域的重要研究方向，涉及到岩土工程、水利工程、环境工程等多个领域。

随着人类对地下资源的开发利用和环境保护意识的提高，对流固耦合渗流规律的研究显得尤为重要。

本文旨在探讨流固耦合渗流的基本原理、研究方法及实际应用，以期为相关领域的工程实践提供理论支持。

二、流固耦合渗流的基本原理流固耦合渗流是指流体在多孔介质中流动时，与固体骨架产生相互作用，导致流体和固体骨架的变形和应力分布发生变化的过程。

这种相互作用涉及到流体流动、固体变形、渗流力学、岩土力学等多个学科领域。

流固耦合渗流的基本原理包括渗流力学原理和岩土力学原理。

渗流力学原理主要研究流体在多孔介质中的流动规律，包括达西定律、渗流速度场、压力场等；岩土力学原理则主要研究岩土体的应力应变关系、变形特性等。

在流固耦合作用下，流体和固体骨架的相互作用会导致渗流场和应力场的相互影响，形成复杂的流固耦合系统。

三、研究方法流固耦合渗流规律的研究方法主要包括理论分析、数值模拟和物理模拟等。

理论分析主要是通过建立数学模型，对流固耦合渗流过程进行理论推导和分析。

数值模拟则是利用计算机软件，对流固耦合渗流过程进行数值计算和模拟分析。

物理模拟则是通过实验手段，对实际工程中的流固耦合渗流过程进行模拟和观测。

在实际研究中，通常将这三种方法相结合，互相验证和补充。

例如，可以通过理论分析和数值模拟，对渗流场和应力场的相互影响进行定量分析；通过物理模拟，对理论分析和数值模拟的结果进行验证和比对。

四、实际应用流固耦合渗流规律的研究在多个领域具有广泛的应用。

在岩土工程领域，可以应用于地下水渗透、堤坝渗透、地下洞室稳定性等领域；在水利工程领域，可以应用于水库蓄水、河岸防护、水利工程安全监测等领域；在环境工程领域，可以应用于地下水污染防治、土壤固废处理等领域。

以地下水渗透为例，通过对地下水位、水头分布、渗流速度等参数的研究，可以掌握地下水的流动规律和渗透能力，为地下水资源开发和利用提供科学依据。

《流固耦合渗流规律研究》范文

《流固耦合渗流规律研究》篇一一、引言流固耦合现象广泛存在于自然界和工程领域中，如地下水渗流、岩土工程、多孔介质流动等。

这些现象涉及到流体与固体之间的相互作用，对渗流规律的研究具有重要的理论意义和实际应用价值。

本文旨在研究流固耦合渗流规律，通过对相关理论的分析和实验研究，揭示流固耦合渗流的基本规律和特点。

二、流固耦合渗流理论基础流固耦合渗流是指在多孔介质中，流体与固体骨架之间的相互作用导致流体在介质中的流动规律发生变化。

其理论基础主要包括多孔介质渗流理论、弹性力学理论、流固耦合理论等。

多孔介质渗流理论是研究流体在多孔介质中流动的规律和特性的理论。

该理论认为多孔介质中的流体流动受到介质孔隙结构、流体性质、外部条件等多种因素的影响。

弹性力学理论则是研究固体在外力作用下的变形和应力分布的理论。

在流固耦合过程中，固体骨架的变形会影响流体的流动规律，而流体的流动又会对固体骨架产生作用力，二者相互影响、相互制约。

流固耦合理论则是将多孔介质渗流理论和弹性力学理论相结合，研究流体与固体之间的相互作用及其对渗流规律的影响。

该理论认为，在流固耦合过程中，流体与固体之间的相互作用会导致介质孔隙结构的变化，从而影响流体的流动规律。

三、实验研究方法为了深入研究流固耦合渗流规律，本文采用实验研究方法。

具体步骤如下：1. 制备多孔介质样品：选用不同类型的多孔介质样品，如砂土、黏土等，进行制备和加工。

2. 设计实验装置：设计一套流固耦合渗流实验装置，包括供水系统、压力传感器、数据采集系统等。

3. 进行实验：将多孔介质样品置于实验装置中，通过改变流体压力、流量等参数，观察流体在介质中的流动规律和特点。

4. 数据处理与分析：对实验数据进行处理和分析，得出流固耦合渗流的基本规律和特点。

四、实验结果与分析通过实验研究，我们得出以下结论：1. 流固耦合渗流过程中，流体与固体之间的相互作用会导致介质孔隙结构的变化，从而影响流体的流动规律。

2. 流体压力、流量等参数对流固耦合渗流规律具有重要影响。

《流固耦合渗流规律研究》范文

《流固耦合渗流规律研究》篇一一、引言流固耦合渗流规律研究是当前多学科交叉领域的一个重要研究方向，涉及力学、物理学、地质学等多个学科。

随着科学技术的不断发展，流固耦合渗流规律在工程实践中应用越来越广泛，如地下水动力学、油气开采、地质灾害防治等领域。

本文旨在通过分析流固耦合渗流的基本原理和规律，探讨其在实际应用中的一些关键问题。

二、流固耦合渗流的基本原理流固耦合渗流是指在多孔介质中，流体和固体骨架之间相互作用而形成的流动过程。

这种流动过程具有复杂的物理机制，包括多孔介质的渗透性、流体与固体骨架的相互作用力等。

在流固耦合渗流中，流体在多孔介质中流动时，会对固体骨架产生作用力，而固体骨架的变形也会对流体流动产生影响。

这种相互作用力的大小和方向取决于多孔介质的性质、流体和固体骨架的物理性质以及流体流动的物理条件等因素。

因此，在研究流固耦合渗流规律时，需要综合考虑这些因素。

三、流固耦合渗流的规律分析1. 多孔介质中的流体流动规律在多孔介质中，流体受到的阻力主要是由多孔介质的孔隙结构和流体的性质所决定的。

流体在多孔介质中流动时，会受到粘滞力、惯性力等多种力的作用，这些力的作用大小和方向取决于多孔介质的渗透性、孔隙大小分布等因素。

因此，在研究多孔介质中的流体流动规律时，需要综合考虑这些因素。

2. 固体骨架的变形规律在流固耦合渗流中，固体骨架的变形也是非常重要的一个方面。

当流体在多孔介质中流动时，会对固体骨架产生作用力，导致固体骨架发生变形。

这种变形的大小和方向取决于流体和固体骨架的相互作用力以及固体骨架的物理性质等因素。

因此，在研究固体骨架的变形规律时，也需要考虑这些因素。

3. 流固耦合渗流的总体规律流固耦合渗流的总体规律是指在多孔介质中，流体和固体骨架之间的相互作用力所形成的流动过程。

这种流动过程具有复杂的物理机制和数学模型。

在实际应用中，需要根据具体的工程问题建立相应的数学模型和求解方法，从而得出可靠的结论和解决方案。

《流固耦合渗流规律研究》

《流固耦合渗流规律研究》篇一一、引言流固耦合问题作为现代科学研究中具有广泛性的领域，是多种复杂自然现象以及工程现象的基本反映。

特别地，流固耦合渗流规律的研究，对于理解流体在多孔介质中的运动、传输和变形过程具有重要意义。

本文旨在探讨流固耦合渗流规律的研究现状、方法及进展，为相关领域的研究提供参考。

二、流固耦合渗流的基本概念流固耦合渗流是指流体在多孔介质中流动时，由于流体与固体骨架的相互作用，导致固体骨架发生变形，进而影响流体流动的过程。

在这个过程中，流体与固体骨架相互依赖、相互影响，共同决定着渗流的运动规律。

三、研究现状目前，流固耦合渗流规律的研究主要集中于岩土工程、石油工程、环境工程等领域。

这些领域中的研究主要集中在多孔介质的力学性质、流体的流动特性以及流固耦合的相互作用机制等方面。

随着计算机技术的发展，数值模拟方法在流固耦合渗流规律的研究中得到了广泛应用。

四、研究方法（一）理论分析理论分析是研究流固耦合渗流规律的基础。

通过建立数学模型，描述流体在多孔介质中的流动过程以及固体骨架的变形过程，进而分析两者之间的相互作用机制。

（二）实验研究实验研究是验证理论分析的重要手段。

通过设计实验装置，模拟实际工程中的流固耦合渗流过程，观察并记录实验现象，为理论分析提供依据。

（三）数值模拟数值模拟是研究流固耦合渗流规律的重要手段。

通过建立数值模型，利用计算机技术对流体在多孔介质中的流动过程进行模拟，从而揭示流固耦合的渗流规律。

五、研究进展近年来，流固耦合渗流规律的研究取得了重要进展。

一方面，理论分析方面取得了突破性进展，建立了更加完善的数学模型，为深入研究提供了理论基础。

另一方面，实验研究和数值模拟方面的技术手段不断更新，提高了研究的准确性和可靠性。

此外，多学科交叉融合的研究方法也为流固耦合渗流规律的研究提供了新的思路和方法。

六、结论与展望通过对流固耦合渗流规律的研究，我们深入理解了流体在多孔介质中的运动、传输和变形过程。

《2024年流固耦合渗流规律研究》范文

《流固耦合渗流规律研究》篇一一、引言流固耦合渗流规律研究是现代工程领域中一个重要的研究方向，主要涉及到流体在多孔介质中与固体骨架的相互作用过程。

随着工业、环境、地质等领域的快速发展，对多孔介质中流体流动的规律及与固体骨架的耦合作用的研究越来越受到重视。

本文旨在通过对流固耦合渗流规律的研究，为相关领域的工程实践提供理论依据和指导。

二、流固耦合渗流基本概念流固耦合渗流是指流体在多孔介质中流动时，与固体骨架发生相互作用，导致流体流动特性和固体骨架的变形或移动过程相互影响的一种现象。

多孔介质中流体流动受到诸多因素的影响，如介质的物理性质、流体的性质以及边界条件等。

而流固耦合作用则涉及到流体对固体骨架的应力作用以及固体骨架变形对流体流动的影响。

三、流固耦合渗流规律研究方法针对流固耦合渗流规律的研究，目前主要采用实验、理论分析和数值模拟等方法。

1. 实验方法：通过设计实验装置，模拟多孔介质中流体的流动过程，观察并记录流体流动特性和固体骨架的变形情况，从而得出流固耦合渗流的规律。

2. 理论分析：基于多孔介质力学、渗流力学等理论，建立流固耦合渗流的数学模型，通过解析或数值方法求解，得出流体流动特性和固体骨架变形的规律。

3. 数值模拟：利用计算机技术，建立多孔介质和流体流动的数值模型，通过模拟流体在多孔介质中的流动过程，得出流固耦合渗流的规律。

四、流固耦合渗流规律研究进展近年来，随着计算机技术的快速发展，数值模拟方法在流固耦合渗流规律研究中得到了广泛应用。

研究者们通过建立更加复杂和精细的数值模型，考虑更多的影响因素，如介质的非均质性、流体的非达西流动等，使得研究结果更加符合实际情况。

同时，实验方法和理论分析也在不断发展和完善，为流固耦合渗流规律的研究提供了更加丰富的手段和思路。

五、研究结论及展望通过对流固耦合渗流规律的研究，我们得出了以下结论：1. 流固耦合渗流过程中，流体流动特性和固体骨架的变形相互影响，使得渗流规律变得更加复杂。

《流固耦合作用下粗糙裂隙岩体渗流及滑移失稳机理研究》

《流固耦合作用下粗糙裂隙岩体渗流及滑移失稳机理研究》一、引言随着工程建设的不断深入，岩体工程稳定性的研究逐渐成为重要的研究方向。

在岩体工程中，流固耦合作用下的渗流及滑移失稳现象是常见的地质灾害，如岩体滑坡、渗水等。

因此，对流固耦合作用下粗糙裂隙岩体渗流及滑移失稳机理的研究具有重要的理论意义和实际应用价值。

本文旨在探讨流固耦合作用下粗糙裂隙岩体的渗流特性及其滑移失稳机理，为岩体工程稳定性分析和灾害预防提供理论支持。

二、研究背景及意义流固耦合是指流体与固体相互作用的过程，其中流体对岩体的渗透作用和岩体变形对流体流动的影响相互影响。

在岩体工程中，粗糙裂隙是常见的地质结构，其渗流特性和滑移失稳机理对岩体的稳定性具有重要影响。

因此，研究流固耦合作用下粗糙裂隙岩体的渗流及滑移失稳机理，有助于深入了解岩体的渗流特性和变形特性，提高岩体工程的稳定性和安全性。

三、研究内容和方法（一）研究内容本文研究内容主要包括：1. 分析流固耦合作用下粗糙裂隙岩体的渗流特性；2. 研究粗糙裂隙岩体的滑移失稳机理；3. 探讨流固耦合作用下粗糙裂隙岩体的稳定性分析方法。

（二）研究方法本文采用理论分析、数值模拟和现场试验相结合的方法进行研究。

具体包括：1. 理论分析：通过建立数学模型，分析流固耦合作用下粗糙裂隙岩体的渗流特性和滑移失稳机理；2. 数值模拟：采用有限元法、离散元法等数值模拟方法，对粗糙裂隙岩体的渗流特性和滑移失稳过程进行模拟分析；3. 现场试验：结合实际工程案例，进行现场试验和监测，验证理论分析和数值模拟结果的正确性。

四、实验结果和讨论（一）渗流特性分析通过理论分析和数值模拟，发现流固耦合作用下粗糙裂隙岩体的渗流特性受到多种因素的影响，如裂隙形态、孔隙率、流体性质等。

其中，裂隙形态对渗流特性的影响最为显著。

不同形态的裂隙对流体的渗透性、流动速度和压力分布等具有不同的影响。

此外，孔隙率和流体性质也会对渗流特性产生影响。

（二）滑移失稳机理研究通过对粗糙裂隙岩体的滑移失稳过程进行数值模拟和现场试验，发现滑移失稳的机理主要包括两个方面：一是由于流体渗透作用导致岩体内部应力分布发生变化，从而引起岩体的变形和破坏；二是由于外部因素（如地震、降雨等）的作用，导致岩体发生滑移失稳。

基于UDEC平台的裂隙岩体边坡渗流应力耦合分析

基于UDEC平台的裂隙岩体边坡渗流应力耦合分析关千军【摘要】根据裂隙岩体边坡渗流的基本特征和基本理论,阐述了离散元软件UDEC 渗流—应力耦合的基本计算理论.利用UDEC作为分析平台建立数学模型,在考虑到渗流场与应力场耦合的作用下对裂隙岩体边坡的稳定性进行了全面分析.对边坡的应力、位移、渗流、裂隙分布等方面进行了研究,重点分析了在渗流与应力耦合的作用下边坡的变形破坏规律及自身稳定性.结果表明利用离散元软件UDEC分析流固耦合是一种有效、可行的方法,同时也为裂隙岩体边坡流固耦合方面的研究提供较高的工程应用价值.【期刊名称】《山西建筑》【年(卷),期】2017(043)036【总页数】3页(P51-53)【关键词】渗流应力耦合;UDEC;裂隙岩体;边坡稳定【作者】关千军【作者单位】太原市市政工程设计研究院,山西太原 030000【正文语种】中文【中图分类】TU4521 概述近年来，由边坡渗流作用引起的工程事故已屡见不鲜。

边坡内的地下水在流动过程中会产生压力并直接作用于岩体上，控制着边坡的应力场，从而间接地影响边坡内节理裂隙的发育情况，对节理裂隙的渗透性起着至关重要的作用；节理裂隙的渗透性反过来又影响着边坡的渗流场，这种应力与渗流相互制约的作用称作二者的耦合作用[1]。

边坡岩体内发育的裂隙一方面影响边坡的表观特征及地质构造作用，另一方面也控制着边坡体内渗流场与应力场的耦合作用。

因此，对裂隙岩体边坡的渗流—应力耦合分析有重要的意义，对其流固耦合的数值模拟研究也变得极为迫切。

UDEC分析平台是一种基于离散单元法，侧重于研究大变形块体模拟的软件。

现已被广泛地应用在岩土工程和岩土力学领域，在节理裂隙边坡应力与渗流耦合方面标新立异并取得了优异的成果及广泛的认可。

Cundall[2]于1971年首次提出了离散单元法，将块体在受力后变形以及根据破坏准则允许断裂的离散单元法称为UDEC。

在国内于1986年的第一届全国岩石力学数值计算及模型试验讨论会上，王泳嘉[3]和剑万禧[4]首次介绍了离散单元法的基本原理及几个实际工程应用的案例。

渗流的流固耦合问题及应用

渗流的流固耦合问题及应用第18卷第5期1999年10月岩石力学与工程学报ChineseJournalofRockMechanicsandEngineering18(5):497～502Oct.,1999渗流的流固耦合问题及应用徐曾和(东北大学岩石破裂与失稳研究中心沈阳110006)博士学位论文摘要渗流过程中的孔隙改变,既影响流体质量,又会引起介质渗透率的变化,导致非线性流固耦合作用。

因此研究应力与孔隙压力共同作用下孔隙改变的机制是重要的。

首先考察了经典渗流力学的基本假定,说明忽略总应力对孔隙改变的作用是经典渗流力学及其他非耦合理论不能研究流固耦合作用的基本原因。

饱和多孔介质是不溶混的混合物,内部孔隙结构对多孔介质整体和孔隙改变的响应方程有影响。

导出了多孔介质的响应方程及有效应力系数Α1、孔隙改变的响应方程和有效应力系数Α2以及孔隙压缩模量Kp。

证明有效应力系数Α1和Α2小于1,Kp小于固体基质的相应模量。

它们都取决于组成多孔介质的固体基质的力学性质和孔隙结构。

还证明了对于孔隙的变形,总应力比孔隙压力更重要。

任何排除总应力影响的理论,其出发点均不尽合理。

研究了流固耦合机制。

在渗流过程中,介质整体变形和孔隙变化是应力和孔隙压力相互作用的结果;孔隙改变会影响两相物质之间的扩散力和流体的质量守恒方程;扩散力和孔隙压力对两相物质的动量守恒有影响。

在混合物理论的基本框架内建立了渗流耦合问题的基本方程,、非线性耦合、非牛顿流体渗流的耦合作用、饱和多孔介质的动力响应等;,应用混合物理;,Biot的三维固结理论只考虑了多孔介质整体压缩(膨胀),、线性的相互作用。

,并得到了平面应变和广义平面应力状态下,。

通过这些解答可以看出,对于稳定渗流,线性耦合理论(Biot理论);,线性耦合理论与经典渗流力学有明显的量的差别,但没有性质的变化。

与弹性:若不考虑渗透率的变化,可压缩流体渗流引起的非线性相互作用以一个高阶小量体现出来,因此与近似的非耦合分析的结果相差不远。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

岩体节理边坡的稳定性受边坡后面水位所影响。

本例中，水位阶段升高直到边坡失稳，当节理内的流体压力增加（节理内的有效正应力降低）以至于边坡面节理的抗剪强度达到极限时，边坡产生破坏。

包含规则节理边坡的几何模型见图8.1。

水位在坡址按四个阶段逐渐升高，四个阶段的水位高程分别是6m、8m、9m和10m。

在每一阶段均进行一个稳定流体的
分析。

图8.1节理边坡渗流问题的几何体
假定下面的材料参数为节理边坡的岩体参数。

8.2 UDEC分析UDEC分析模型见图8.2, 本问题通过设置SET flow steady命令进行稳定流模拟,水位的上升采用BOUND pp pygrad命令改变每个阶段的流
体压力梯度来实现.数据文件见8.3节.
首先,边坡在重力荷载下达到平衡,然后在右侧的水位上升到坡址上高度6m处;左侧的水位保持在坡址的水平上,在这个流体压力条件下,边坡是稳定的.此条件下
的稳定流态见图8.3.
接下来,右侧的水位上升到8m,8m水位高度条件下稳定流态见图8.4.系统再次达到平衡.然后水位上升至9m.流体状态见图9.5,边坡仍旧保持稳定.
最后,水位达到边坡的顶部,此条件下的流态见图8.6.在10m水位时,边坡滑动,见
图8.7的岩石楔形体.
注意.本问题采用了命令SET delc off进行设置,它将阻止接触后面破坏的岩石楔形体被删除,否则,与这些接触相关的区域将变成外部区域,并且流体压力将会消失. 当楔形体产生大位移时.这将很难确定实际的流体压力.保守的估计是假定流
体压力没有消.
Example 8.1 SLOPEFLO.DAT
; ================================================
;
; --- fluid flow test run ---
; --- slope : 10 m high ---
;
; --- 2 joint sets : 20 and 80 deg.
; friction = 30 deg.
;
; --- r.h.s. water level : 6 m --- no failure ---
;
8m
;
9m
;
10 m --- failure ---
;
; ================================================ ;
round 0.05
;设置倒圆角数值
set delc off
;设置破坏体不被删除
;
block 0,-5 0,0 5,0 11,10 22,10 22,-5
;定义问题几何体的外轮廓
;节理切割块体
jset 20,0 100,0 0,0 2,0 (5,1)
;倾角20
jset 80,0 100,0 0,0 3,0 (5,0)
;倾角80
;
del area 0.1
;删除面积小于0.1的块体
;;plot block
;显示块体
gen auto 10
;生成变长10的三角形单元
;;plot block zones
;显示单元划分网格
;参数设置
prop mat=1 dens=0.0025 k=16667 g=10000
prop mat=1 jkn=10000 jks=10000 jfric=45.0
prop mat=1 jperm=1.0e8 azero=0.0005 ares=0.0002 ;
change jmat=1 jcons=2
;初始应力条件
insitu str -0.125 0 -0.25 ygrad 0.0125 0 0.025
;边界条件
bound xvel=0 range -1,1 -6 1
bound yvel=0 range -1,31 -6 -4
bound xvel=0 range 21,23 -6 11
;设置重力加速度
grav 0 -10
;记录时步和点的位移、不平衡力hist n=100 xdis 11,10 ydis 11,10
hist unbal
;求解和保存
solve rat 1e-5
save slfl1.sav
;查看计算结果
;plot block dis yel
;plot hist 1
;plot hist 2
; --------------------------------------
;
; flow --- r.h.s. water at y=6m
;
; --- no failure ---
;
; --------------------------------------
fluid dens 0.001
;流体密度
bound imperm range -1,31 -6 -4
;不可渗透边界
bound pygrad -0.01 range -1,1 -6 1
;左侧水位0,x=0,y=-6~1
bound pp=0.06 pygrad -0.01 range 21,23 -6 6 ; 右侧升高水位到6m,x=22,y=-6~6
;
set flow steady
;设置稳定流
;
hist xdis 5.92 1.54 ydis 5.92 1.54
;记录点位移
reset displ jdis
prop mat=1 jfric=25
cyc 500
save slfl2.sav
;计算结果显示
;plot bl vf yel
显示流体矢量-黄色
;plot bl dis yel
显示块体位移
;plot bl pp yel
显示节理孔隙水压力
; --------------------------------------
;
; flow --- r.h.s. water raised to y=8m
;
; --- some slip, but no block failure ---
;
; --------------------------------------
bound pp=0.08 pygrad -0.01 range 21,23 -6 8 cyc 500
save slfl3.sav
;
;plot bl vf yel
流体矢量-黄色
;plot bl pp red dis yel
同时显示孔隙水压力(红色)和位移(黄色)
;plot bl sh yel
显示节理剪切位移数量
; --------------------------------------
;
; flow --- r.h.s. water raised to y=9m
;
; --- no failure ---
;
; --------------------------------------
bound pp=0.09 pygrad -0.01 range 21,23 -6 9
cyc 1000 save slfl4.sav
;
; --------------------------------------
;
; flow --- r.h.s. water raised to y=10m
;
; --- failure ---
;
; --------------------------------------
bound pp=0.10 pygrad -0.01 range 21,23 -6 11
cyc 15000
save slfl5.sav
return。