基于PFC3D的鱼鳅坡滑坡运动过程分析

以贵州省开阳县鱼鳅坡滑坡为研究对象,采用颗粒流离散元（PFC3D）对其破坏运动过程进行数值模拟。采用Ball-Wall建模方法建立滑坡模型,对滑坡不同关键部位颗粒进行位移、速度监测,阐明其破坏运动特征。结果表明,降雨为鱼鳅坡滑坡的直接诱发因素。该滑坡在破坏初始阶段以蠕滑变形为主,随着变形量的增加,滑坡体不断挤压坡脚,滑坡岩土体到达应力平衡极限,坡脚产生剪切破坏,并向上牵引发展,滑坡发生整体滑动,斜坡变形破坏模式为蠕滑-拉裂,按照力学条件为牵引式破坏。滑坡滑动最高时速12.4 m/s,最大滑移80 m,滑动阶段持续50 s。研究成果可为对该类滑坡影响范围预测,以及工程措施的制定具有一定的参考意义。     

基于颗粒流方法的堆积层滑坡运动过程模拟

以颗粒流离散元为研究方法对勉县杨家湾十组堆积层滑坡破坏方式与运动过程进行数值模拟研究。通过PFC^2D双轴模拟试验所标定的岩土体宏观模拟参数与室内试验所获取的宏观实测参数进行对比,确定堆积层滑坡所需的颗粒细观参数,然后将标定的细观参数代入堆积层滑坡模型,对滑坡破坏方式及运动过程进行模拟研究。结果表明：滑坡破坏在初始阶段蠕滑变形累积,滑坡体挤压坡脚,直至坡脚产生剪切破坏,并向上牵引发展,使得滑坡整体顺接触面破坏下滑并堆积于坡脚,表现为典型的牵引式渐进破坏,结果与实际情况基本吻合。研究认为,采用颗粒流方法对堆积层滑坡破坏与运动过程的模拟研究具有较高的适用性,对该类滑坡防治具有一定的参考意义。     

PFC滑坡模拟二、三维建模方法研究

滑坡运动过程模拟避免了试验尺度与监测手段的限制,可以详细观察破坏过程,是定量评估滑坡灾变风险的重要研究手段。开展滑坡运动过程模拟的首要工作是建立滑坡模型。颗粒流程序（PFC) 虽然是滑坡运动过程模拟应用最广泛的程序之一,但在建立滑坡模型的前处理方面较弱,使得其在滑坡运动过程模拟中的推广应用受到限制。对此,本文指出了Ball-Ball和Ball-Wall两种建模方法的适用性、优缺点及滑体滑床边界确定方法,并以2014年地震触发的红石岩滑坡为例,以Brick填充法为基础,从获取地形数据、确定滑体和滑床区域、建立滑体和滑床几何模型、生成颗粒模型4个步骤出发,提出了基于数字等高线地形图建立PFC二、三维复杂滑坡模型的前处理方法,弥补了PFC软件前处理的不足,从而为今后基于PFC的滑坡运动过程模拟提供有益的帮助。     

颗粒流离散元在滑坡运动过程模拟中的应用

介绍了颗粒流离散元的分析原理,以某岩质边坡为实例,构建了计算模型。根据实际资料,在计算模型中加了节理,节理的力学参数弱于岩体。在下方的抗滑挡土墙完全失去作用的情况下,用此模型对该滑坡的运动过程进行了计算预测,分析了整体的变形、不同部位的位移和应力等运动特点,上部粉质粘土区系拉裂引起的下滑,下部碎石土区因挡墙损坏而产生沿节理面的滑动,但拉应力变化不大。指出了在该滑坡存在2个滑动面,并针对2个滑动面提供了相应的治理建议。     

金沙江白格滑坡运动过程特征数值模拟与危险性预测研究

2018年10月、11月于金沙江川藏交界处江达县波罗乡白格村先后发生两次体积约2400×104 m3和850×104 m3的滑坡,两次滑坡平均运动距离1400 m,堵塞金沙江形成堰塞湖。首次形成的堵江滑坡坝天然溃决,未造成人员伤亡;然而第2次滑坡堵塞第1次滑坡自然溃口,导致堰塞湖库容迅速增加到3.85×108 m3。政府部门立即开展抢险工作,通过人工修建溢洪道的方法成功泄洪,极大程度上降低洪水风险。本文利用PFC3D颗粒流软件模拟两次滑坡的发生、运动、堆积过程,并在反演结果的基础上对白格滑坡滑源区残留潜在不稳定部分未来失稳的运动路径和堆积范围进行预测,对其危险性进行科学评价。结果表明:（1）滑坡在重力作用下失稳,除了受初始势能的影响外,微地貌也是决定滑坡运动路径与距离的关键因素之一;（2）PFC3D颗粒流数值模拟方法适用于类似于白格滑坡这类碎屑流类型的滑坡,两次滑坡反演得到的堆积厚度、堆积范围均与真实结果相近;（3）利用两次事件反演所得参数,可以预测若滑源区潜在不稳定部分同时失稳,则形成约70 m高的滑坡坝,可能再次堵塞金沙江。     

颗粒流方法及PFC2D程序

归纳总结了颗粒流模拟方程产生的背景,比较了与其他模拟方法的异同之处,介绍了颗粒流模拟方法解题的步骤及其在岩土工程领域的应用实例。     

基于颗粒流模拟的堆积层滑坡过程分析

基于颗粒流离散元程序模拟了堆积层滑坡变形破坏及滑动的全过程。结合现场地质调查资料,总结了拜殿乡Ⅰ号滑坡的成因机制及影响因素,分析表明该滑坡是在内外因共同作用下的结果,其主要原因是滑体物质松散,雨水入渗形成软弱滑动带,导致坡体的变形滑动。采用数值试验确定土体细观参数与宏观力学性质的关系,建立拜殿乡Ⅰ号滑坡数值模型,考虑降雨对滑带土强度的软作用化,对滑坡变形破坏及滑动全过程进行了模拟。模拟结果表明：滑带土强度降低后坡脚出现剪切破坏集中区,后缘处出现拉张裂破坏集中区,随着坡体缓慢的蠕动变形,坡脚的裂隙范围逐渐扩大,坡体中后缘出现多条拉张裂隙;根据滑体的运动速度将滑动过程分为蠕动变形、滑动加速、滑动减速和缓慢堆积4个阶段;滑坡滑动情况和最终堆积状态表明拜殿乡Ⅰ号滑坡严重威胁龙潭村小学及乡村公路的安全,模拟结果对该滑坡的工程治理具有重要的指导价值。     

基于PFC2D的贵阳西二环滑坡破坏过程模拟

近年来,随着计算机技术的发展,数值模拟方法在滑坡领域得到广泛应用.目前的数值模拟方法中,主要包括确定性分析方法和非确定性分析方法两类,而确定性分析方法又可分为连续介质分析与非连续介质分析方法.其中,连续介质数值分析方法有有限单元法、边界元法、有限差分法等,非连续介质分析方法有块体离散元法、颗粒离散元法、关键块体理论、不...     

大光包滑坡运动特征及其过程分析

大光包滑坡是512汶川地震触发的规模最大的滑坡。滑坡的形成机制和运动学特征引起国内外学者的广泛关注。本文基于滑坡体运动-堆积特征的现场调查与分析,对滑坡的运动学特征进行了系统的阐述。为了再现滑坡的破坏和运动过程,应用物理模拟试验方法进行了对比研究。(1)通过对滑坡擦痕、标志性地物、滑坡堆积体植被、坡表块石倾向、滑坡裂缝、以及其他滑坡典型运动特征调查与分析研究,确定了滑坡的滑动方向、滑动距离、运动速度、运动特征值等滑坡运动特征参数,并系统地论述了滑坡扬尘、滑坡气浪、滑坡舌前缘岩性混杂堆积带以及滑坡碎屑流等典型运动-堆积现象的形成机理; (2)滑坡运动过程物理模拟试验表明,大光包滑坡破坏运动过程中存在前缘锁固段岩体剪断迸射、滑坡气浪、滑坡扬尘和急刹车效应等地质现象,且滑坡运动过程具有显著的阶段性; (3)大光包滑坡的运动过程可以概括为以下4个主要阶段,即:快速启动高速滑动急刹车制动拆离滑动。     

基于Flac-3D的盐源县玻璃村滑坡变形破坏机理数值模拟

2018年7月19日凌晨5点40分左右,四川省凉山州盐源县桃子乡玻璃村发生一起特大型滑坡地质灾害,为了查明其变形破坏机理,通过现场调查、室内外试验及数值模拟等方法,在查明滑坡地质背景基础上,对其发育特征进行分析,并基于Flac-3D有限元软件强度折减法计算滑坡天然及饱和工况下的稳定性系数.通过对比分析两种工况下塑性区及...     

延安新区马兰黄土湿陷特性的PFC2D模拟

为了深入了解延安新区马兰黄土的粒间作用,揭示其湿陷机理。首先利用扫描电子显微镜观察了黄土的微观结构特征,并且利用IPP图像处理软件对黄土颗粒定量分析,论证了PFC2D软件模拟黄土粒间作用的可行性。然后,基于离散元理论,利用颗粒流PFC2D程序模拟了延安新区不同埋深原状黄土颗粒的分布情况及颗粒间相互作用,以一种直观的方法再现了天然黄土的微观结构特征。对比分析了不同深度黄土颗粒分布及接触力分布与湿陷性的关系。结果表明:骨架颗粒间的作用决定了土体的结构强度,尤其以团粒结构的接触力为主。随着深度的增加,黄土骨架颗粒由点接触方式逐渐过渡到团粒重叠接触,粒间接触力明显增强。PFC2D模拟方法为黄土湿陷性成因机理研究提供了新的思路。      

基于FLAC3D的饱和-非饱和渗流分析

当使用基于有限差分法的FLAC3D计算渗流场时,自由面位置过高,非饱和区渗透系数被定义成常数,这与实际情况不符。运用变单元渗透系数法,对其进行二次开发,并提出了修正的变单元渗透系数法,对非饱和区的渗透系数进行较为严格的理论修正,得出准确的连续光滑的自由面,求解出饱和-非饱和渗流场。利用有限元法与该方法进行对比,论证了该方法的可行性和有效性。提出了附着接触面的方法,解决了三维复杂边界中寻找溢出点的困难;并以实际的三维复杂工程为例,计算了某边坡的渗流场分布情况,经验证,结果准确。研究成果对研究饱和-非饱和渗流以及该软件中与渗流有关的计算具有有指导意义。     

复杂开挖过程FLAC3D力学仿真代码生成系统研究

针对直接采用FLAC3D编制复杂开挖过程力学仿真程序所遇到的困难,开发了FLAC3D力学仿真程序代码生成系统(generating code system,简称GCS)。介绍了如何利用矿用地质软件Surpac,按空间、次序、结构功能及工程地质体属性对复杂开挖过程进行剖分,形成可编程的组件;着重解决了模型单元质心定位、单元捕捉、重组和应用程序变量到FLAC3D变量、函数、语句、数据结构之间映射等问题;使用VC++开发并给出了一个特定回采开挖过程的FLAC3D力学仿真代码生成系统实例,该应用程序以模型单元数据库为数据源,通过一系列与用户交互的页面,由用户对回采开挖过程所需步骤和参数进行选择和确认,并将这些步骤和参数传递给仿真代码生成系统,由系统自动完成FLAC3D力学仿真代码的生成。FLAC3D力学仿真程序代码生成系统架起了地质三维几何造型软件与三维计算软件之间的桥梁,可大大提高编制FLAC3D三维力学仿真代码的效率,降低程序编制的出错率,有利于应用FLAC3D软件更为深广地解决复杂岩土开挖过程中的力学问题。     

TOUGHREACT搭接FLAC3D算法

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

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

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

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

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

基于扩展有限元法对土体滑坡破坏过程的模拟

边坡失稳一般经历滑裂面萌生、扩展和贯通的过程。近年发展起来的扩展有限元法（XFEM）可以合理地模拟这一过程,但在实际应用中仍未克服需要预先设定滑裂面的起始位置,以及滑裂面前端扩展方向的判断精度较低等缺陷。首先介绍滑裂面萌生和扩展的判断和程序实现机制。然后提出了确定滑裂面起始位置的方法,重点描述了其中的自动判断裂缝萌生位置的方法,即根据单元内的应力状态求得拉应力水平和剪应力水平,根据其相对关系并结合应力历史进行判断。接着介绍了联合运用扇形控制域和圆形控制域确定滑裂面前端扩展方向的新方法。该方法可提高XFEM对滑坡破坏过程的模拟精度。最后通过两个边坡失稳算例的计算分析验证所提方法的可靠性。     

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

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