GeoStudio软件在土坡饱和—非饱和渗流分析中的应用

在对饱和-非饱和渗流数学模型简要分析的基础上,对影响渗流有限元计算分析结果的一些参数选取方法进行讨论,如非饱和土的水土特征曲线、渗透系数及毛细水头等。利用GeoStud io软件对降雨条件下的土质边坡内部渗流场进行模拟,分析地下水位线、坡体饱和度、压力水头随降雨时间的变化规律。模拟结果与实际情况比较吻合,可为土质边坡稳定性的分析提供参考依据。     

基于饱和-非饱和渗流理论的基坑防渗措施分析

应用饱和-非饱和渗流理论,对某临河基坑进行了渗流计算,评价了基坑的渗透稳定性,并确定了基坑的最佳防渗措施及相应的参数。计算结果表明：无防渗措施时,基坑坡面土体在渗透压力作用下是不稳定的,必须采取防渗措施。从便于施工的角度考虑,采用钢板桩防渗,其位置定在距基坑右边坡6.0 m处。钢板桩长度对渗流场影响显著,钢板桩应穿过粉砂层,这样才能起到较好的防渗效果,有利于基坑边坡的稳定;同时也不宜太长,太长不仅浪费材料,效果也不会更好,最佳长度应为15.4 m。     

降雨条件下岩土饱和-非饱和渗流分析

基于Buckley-Leverett两相渗流方程,提出了新的岩土饱和-非饱和渗流数学模型,利用有限差分方法给出了隐式压力显示饱和度数值求解方法,编制了饱和-非饱和渗流计算程序.结合工程实例,模拟了降雨入渗条件下边坡岩体渗流场孔隙压力变化和含水饱和度变化,模拟结果验证了所提出的模型对饱和-非饱和渗流的有效性.     

降雨条件下土坡饱和-非饱和渗流分析

分析了降雨入渗的过程,提出了降雨概念模型:讨论了如何用饱和-非饱和渗流数值方法分析边坡在不同雨型下的渗流场;分析结果表明:雨水入渗过程受初始渗流场影响很大,初始含水量的大小影响湿润锋的推进。暴雨情形下,雨水入渗缓慢,浅部土体保持一定吸力,雨停后渗流场仍有较大变化。     

考虑饱和-非饱和渗流的土坡极限分析

给出了均质土坝的下游坝坡的安全系数的计算方法。均质土坝的下游坝坡有可能会通过滑裂面发生破坏。土坝中的土坡通常处于非饱和状态。非饱和土坡的安全系数计算需要考虑吸力对抗剪强度的贡献以及土坡中的非饱和渗流。给出处于饱和-非饱和渗流状态下的土坡的安全系数,有助于评价均质土坝的安全系数。下游坝坡的安全系数计算方法有：极限平衡法、上限解法和下限解法,该算法适用于非饱和土坡,且是在饱和土坡安全系数的计算方法上修正得到的。算例中非饱和土坡安全系数的计算考虑了吸力对抗剪强度的贡献。考虑非饱和渗流理论的土坡安全系数计算方法通常更加接近现场实际情况,并且对于同一坝坡,考虑非饱和渗流计算出的土坡安全系数要比饱和渗流理论计算出的安全系数大。     

考虑饱和-非饱和渗流作用的土质边坡稳定性分析

饱和-非饱和土的研究一直是岩土工程界的一个热点问题。土体中孔隙水压力的存在使得对岩土体的力学性状的分析变得更加复杂化。基于土体渗透系数与基质吸力之间的相互联系,对理想的各向同性土的水库堤坝进行了非稳定渗流场的数值模拟,得出不同水位、同一水位因不同骤升速率而产生的渗流场的变化规律,并对土体固-液耦合研究进行了初步探讨。把堤坝渗流场的动态变化在时间和空间上进行离散,分别引入到坝坡稳定性的分析中,以极限平衡原理和Morgenstern-Price条分法为基础对坝坡产生滑移的最小安全系数变化规律进行了分析,从而得出渗流-滑坡的内在联系。     

裂隙性膨胀土饱和-非饱和渗流分析

降雨入渗条件下的裂隙性膨胀土边坡稳定性与其渗流场分布密切相关,对渗流分析中的积水深度进行了数值模拟,结果显示,不同的积水深度对渗流场的影响很小,因此,在进行裂隙性膨胀土渗流分析时,积水深度为0的假定是可以接受的。基于膨胀土开裂裂隙的规模及渗透性的不同,提出了考虑裂隙系统的膨胀土边坡渗流分析方法,计算结果表明,土体开裂显著改变了膨胀土内部的渗流场分布。该方法可以更好地模拟坡渗流场随时间的变化规律,以及更合理地解释了降雨入渗引起的膨胀土边坡浅层滑动机制。     

考虑饱和-非饱和渗流、温度和应力耦合的三维有限元程序研制

在探索多场耦合理论基础上,研制出了一个考虑饱和-非饱和渗流场和温度场对应力场作用的三维弹塑性有限元程序,该程序采用了以下简化假设：不考虑力学过程和渗流过程对温度场产生的影响,而只考虑变温对力学过程产生的影响;不考虑力学过程和温度场对渗流场产生的影响,而只考虑渗流场对力学过程产生的影响。经与弹性力学解析解对比,验证了该程序计算结果是合理的。     

岩体裂隙饱和/非饱和渗流机理初步研究

基于裂隙局部开度为平行板模型的假定,按一定分布抽样生成裂隙面上开度随机数,并分别对单裂隙饱和及非饱和渗流进行了数值模拟.裂隙饱和及非饱和渗流在渗流机理上都是由于水势差引起的,只不过在非饱和渗流情形下,湿润相水压力小于非湿润相空气的压力.饱和渗流时,水流主要集中在大开度,即“沟槽”中;而非饱和渗流时,水流运动多发生在小开度,即“狭缝”中.     

三维饱和-非饱和渗流场求解及其逸出面边界条件处理

对三维饱和-非饱和渗流场的有限单元法求解进行了详细论述,在边界条件中明确列出了非饱和渗流逸出面边界条件的数学表达式。为了理论上能够严密地给出非饱和逸出面的边界条件和数值处理,根据渗流场固有的物理特性,定义了极限含水率和非饱和渗流逸出面蒸逸系数两个渗透基本特性参数。针对渗流场饱和与非饱和逸出面分界线(逸出线)位置事先不知道的问题,提出在这种边界非线性渗流场迭代求解过程中,应根据中间解通过边界条件类型的转化来准确地甄别这一分界线,完整地获得精细地处理饱和及非饱和渗流逸出面边界条件的思路与方法。     

饱和-非饱和非稳定渗流的数值模拟

介绍了饱和-非饱和渗流的计算程序数值模拟方法,该方法模拟了降雨过程中边坡内孔隙水压力的变化,克服了传统降雨边界处理方法的缺陷,考虑了降雨过程中地表边界的条件转化,为以后的径流渗流耦合的模拟提供了良好的基础,为降雨诱发滑坡研究提供了定量分析手段。     

TOUGHREACT搭接FLAC3D算法

多孔介质THM多场耦合分析在诸如核废料地质储存、干热岩和CO2地质埋藏等地质科学领域广泛应用,鉴于其计算的复杂性往往需要借助数值手段。基于角点-中心点网格转换以及三维三点拉格朗日插值方法,通过编制程序将TOUGHREACT和FLAC3D搭接,实现了多相多组分流体-温度-岩石力学的耦合计算分析。为了验证程序的效率及适用性,构建了Salah场地深部咸水层CO2储存数值模型,模拟CO2连续注入3a的地下流体运移和储、盖层的岩石力学变化,并将地表位移与inSAR实测值进行了对比分析。以此为基础进一步模拟了中国西部某盆地CO2地质储存三维筒注模型。结果表明:搭接程序能平滑连接TOUGHREACT和FLAC3D,其考虑岩石力学影响下的注入压力和地表位移等成果符合理论及实际监测规律,同时得到岩石孔隙率、渗透率随CO2注入时间增加而变大的结论。     

边界面模型在FLAC3D中的开发及验证

利用FLAC3D程序提供的二次开发平台UDM在VC++编译环境下实现了一基于修正剑桥模型的简化边界面模型的开发。介绍了边界面模型的基本原理,并给出了开发关键技术和具体实施方法。通过三轴加卸载试验、不排水静三轴试验、不排水动三轴试验等几种数值试验与修正剑桥模型进行了对比研究。计算结果验证了所开发的边界面模型的正确性及相对修正剑桥模型的优越性。     

基于FLAC3D和DEM数据的缓倾斜煤层开采沉陷分析

宁夏石嘴山矿区位于西部黄河流域,其煤矿采空区沉陷导致地表生态和环境问题频发,对其采煤沉陷分析将对西部黄河流域煤矿区的环境修复有一定的积极作用。为研究缓倾斜煤层采空区围岩应力与位移场演化特征,以宁夏石嘴山矿区为对象,基于FLAC3D数值模拟软件,建立缓斜煤层开采三维数值模型,计算分析采空区围岩应力、塑性区及位移变化规律,并基于两时相DEM叠加统计分析地表位移变化,与数值模拟结果进行相互验证。结果表明:地下开采引起应力重分布,采空区顶板及煤柱出现明显的应力集中现象,最大主应力呈现从煤层顶板向地表递减的变化趋势;越靠近采空区顶部的岩层垂直位移越大,随着远离采空区逐渐减少,开采完成后地表垂直位移最大值约12 m;随着采空区面积的不断增大,采空区四周及角隅处塑性区逐步延伸扩大,且以剪切破坏为主;地面沉陷盆地不对称,2个沉降中心均发生在沉陷盆地中部且偏下山方向,下山方向比上山方向影响范围更大;数值模拟计算的沉降量与两时相DEM叠加统计分析的变化量结果及趋势基本一致,研究成果可为煤炭安全开采提供参考依据,为地表沉降监测提供新方法。      

无压渗流分析的复合单元法

针对无压渗流问题,在传导矩阵调整法的基础上,对穿越自由面的单元用复合单元来处理,每个复合单元具有两套结点水头,一套用于饱和区的水头插值;另一套则用于无水区。通过变分原理,推出了控制方程并在程序中实现。应用复合单元法可以较好地解决自由面的插值问题,消除了穿过自由面单元的水上部分对自由面附近结点水头的贡献,同时对逸出面边界条件的处理也更加方便和准确。算例表明了该方法的有效性。     

基于FLAC3D含SMP强度准则黏弹塑性模型的二次开发

在分析岩土工程时通常会用到强度准则,在众多强度准则中,SMP准则的物理意义最为明确,其屈服函数表达式是非线性形式,对其进行二次开发研究具有较为重要的意义。利用FLAC3D的二次开发平台,使用VC++环境在FLAC3D软件中实现了含SMP强度准则的黏弹塑性模型的开发,其中的SMP强度准则推广到了黏性土,通过一个简单的算例验证了程序编制的正确性和可靠性,所采用的相关二次开发的思路可以为其他强度准则及流变本构模型的二次开发提供参考。     

基于FLAC3D的煤矿采空区变形数值模拟分析

本文选取拟建国道丹阿线新华至英安段煤矿采空区作为工程背景,在有限差分软件FLAC^3D基础上,建立采空区模型。研究了在自重应力作用下,采空区上覆岩层及地表变形规律,并且得出了开采过程中的应力-应变规律。分析结果表明：由于地下煤层的开采,原有的地质状态发生了变化;最大主应力发生在采空区的边界以及顶板位置,出现了明显的应力集中现象,在顶板的上部,其最大主应力呈现出从上至下逐渐增大的变化趋势,主应力主要表现为压应力;在采空区的顶板位置最容易发生沉降变形破坏,越靠近采空区的顶部位置,产生的垂向位移则越大,岩层随着远离采空区,产生的垂向位移则逐渐减少;塌陷盆地不对称,地表沉降位移的最大值发生在塌陷盆地中部、偏于下山方向;多层采空区开采结束后,沉降位移达到最大值,最大位移达到了791 mm。通过概率积分法的理论计算,得到研究区域最大沉降值为830 mm,得出概率积分法和FLAC^3D数值模拟方法两者的结果相互吻合。     

尾砂胶结充填体蠕变模型及在FLAC3D 二次开发中的实验研究

关于岩石蠕变特性的研究已有诸多成果,但对胶结充填体的蠕变特性还缺乏系统地研究。在室内充填体单轴蠕变试验的基础上,利用Hoek-Kelvin模型表征充填体的蠕变特性,运用粒子群优化算法,对选用的蠕变模型的参数进行辨识,并研究各蠕变参数对应力水平的敏感程度。研究结果表明,在不同分级应力水平加载下,弹性模量EH和黏滞系数η变化较小,而参数EK对应力水平变化较敏感。利用FLAC3D软件二次开发所建立的蠕变模型,用开发的模型进行分级加载蠕变的数值模拟,计算结果与试验数据吻合。