岩石圆盘径向压缩破坏的Voronoi子块体单元DDA方法模拟

非连续变形分析(DDA)方法是计算离散可变形块体系统力学行为的数值计算方法，可通过子块体单元DDA方法模拟岩石的开裂破坏。考虑到Voronoi多边形颗粒与细观尺度下岩石矿物晶粒形态的相似性，提出一种基于随机圆的Voronoi颗粒单元模型生成方法；并通过完整及带预制裂纹岩石圆盘径向压缩破坏的模拟，验证岩石破裂问题Voronoi子块体单元DDA模拟方法的适用性。结果表明，当子块体单元数较小时，圆盘表现出更高的整体强度；随着子块体单元数的增大，起裂处位置更接近真实，开裂破坏路径更清晰；子块体单元数较大时不同倾角预制裂纹圆盘破坏的模拟结果与实验结果高度吻合，并能有效反映圆盘中心加工小孔对开裂破坏路径的影响。使用Voronoi子块体单元DDA方法能够有效模拟岩石的开裂破坏过程，为进一步开展基于Voronoi颗粒单元模型的岩石开裂破坏模拟创造了条件。     

岩石力学中的Fuzzy数学方法

鉴于岩石力学问题中的许多因素具有强烈的“模糊性”,本文利用Fuzzy数学来研究工程中的岩体力学问题;并针对露天矿山岩体变形破坏问题,应用Fuzzy概率测度理论导出了岩体变形破坏的Fuzzy概率公式,且给出了数值方法,实例计算表明,本方法适用于分析和研究岩石力学尤其是矿山岩体力学问题。     

岩石破裂滑动过程及其声发射数值模拟

运用岩石破裂与失稳过程分析RFPA2D系统, 对滑动形成过程进行了数值模拟研究。滑动岩体似为非均匀弹 -脆岩石材料, 模拟结果再现了滑动形成从变形到破坏直至失稳的全过程及其声发射规律。     

采动诱发岩体移动破坏过程数值模拟研究

本文基于岩石介质的宏观非线性主要是由非均质性和各向异性造成的, 应用新的数值计算软件RFPA(2D), 对采动引起岩体失稳破坏的全过程进了数值模拟研究。     

考虑岩体破坏分区的岩石爆破爆炸荷载历程研究

针对岩石爆破爆炸荷载历程中未联合考虑岩石爆破动态过程和炮孔周围岩体破坏分区的不足,开展了考虑岩体破坏分区的岩石爆破爆炸荷载历程及其适用性研究。联合岩石爆破动态过程和岩体破坏分区的理论解,推导了考虑岩体破坏分区的岩石爆破爆炸荷载理论公式,比较了考虑岩体破坏分区的岩石爆破爆炸荷载历程与实测炮孔爆炸压力曲线,开展了单孔爆破现场试验和相应条件下3种爆炸荷载工况的数值模拟,并对爆破振动现场实测和数值模拟结果进行了对比。研究结果表明：考虑岩体破坏分区的爆炸荷载历程包括上升段和衰减段Ⅰ、Ⅱ、Ⅲ,上升段持续时间极短,衰减段持续时间较长且主要由填塞情况控制;考虑岩体破坏分区的爆炸荷载历程理论计算结果与实测炮孔爆炸压力曲线的变化趋势一致,验证了考虑岩体破坏分区的岩石爆破爆炸荷载理论公式的可靠性;考虑岩体破坏分区的爆炸荷载工况下,单孔爆破振动波形的数值模拟结果与现场实测结果的主要特征一致,该荷载工况下质点峰值振速计算结果与现场实测值偏差率最小,绝大部分在7%以内,表明了其应用于数值模拟的优越性;考虑岩体破坏分区的爆炸荷载可随岩石爆破系统条件的变化而动态调整,其可靠性好、适应性强、应用效果佳。     

EVALUATION ON STRESS INTENSITY FACTOR OF CRACK UNDER DYNAMIC LOAD USING NUMERICAL MANIFOLD METHOD

相较于传统有限元,数值流形方法（numerical manifold method, NMM） 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子（dynamic stress intensity factor,DSIF） 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性.     

动载下裂纹应力强度因子计算的数值流形元法

相较于传统有限元,数值流形方法（numerical manifold method, NMM） 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子（dynamic stress intensity factor,DSIF） 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性.      

高边坡岩体渐进性破坏粘弹塑性有限元数值模拟

根据边坡实际地质模型 ,基于弹塑性与粘弹 -粘塑性理论的本构方程 ,通过有限元模拟分析 ,定量地揭示和模拟再现了高边坡岩体破裂、变形、破坏及失稳前后锁固段岩体渐进性破坏的机制和过程 ,探讨了高压水流作用下滑坡启程剧动的破坏机理。     

高陡岩质边坡稳定性三维离散元分析

某高陡岩质边坡地质条件复杂、软弱结构面发育、开挖高度大、坡度陡、临空面多,为边坡变形提供了有利的空间,边坡多处出现失稳破坏迹象。通过对边坡工程地质条件调查,岩体结构特征和边坡开挖等影响因素的分析,认为边坡变形主要发生在强风化强卸荷岩体内,受软弱结构面的控制比较明显,表现为结构面组合控制的块体变形失稳破坏模式。采用3DEC数值模拟软件,模拟了边坡开挖后坡体变形特征,数值模拟结果表明,边坡浅表层块体以及控制性块体稳定性差,可能导致边坡产生整体失稳。     

澜沧江某水电站右坝肩工程边坡倾倒变形问题的数值 模拟研究

澜沧江某水电站处于青藏高原东部边缘地带,属于高山峡谷地形地貌,高地应力环境,岩体卸荷裂隙很发育,使得倾倒变形和岩体质量、断裂活动及地震构造一样成为影响工程边坡岩体稳定的主要因素。针对工程边坡的大变形问题可采用离散元的数值模拟分析方法。通过建立理论开挖和工程边坡开挖离散元模型,可分别得出倾倒变形破坏机理发展过程为初期弱倾倒变形岩体的层内剪切错动、强倾倒变形岩体的层内拉张变形、强倾倒变形岩体的切层张-剪破裂及极强倾倒破裂岩体的折断张裂(坠覆)破裂和工程边坡的变形范围、确定开挖面及加固方式等。通过工程边坡模型的计算结果和现场地质调查成果的比较表明,计算结果和实际情况基本吻合。     

非连续变形分析方法的改进及应用1）

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现. 其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征. 研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳. 改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

应用真三维岩体建模仿真技术推进岩石力学教学改革

随着计算机技术的发展和研究项目最新成果的取得, 把岩体真三维建模仿真技术引入到岩石 力学教学改革与实践中, 挖掘教学资源潜力; 开展岩体结构、岩石变形与破坏过程的数值试 验, 再现许多物理实验所不能观察到的力学现象, 例如岩体的空间真三维结构、岩石破坏过 程块体滑落等, 从而使学生对岩石的变形与破坏过程有更加清晰的认识, 提高岩石力学实验 的教学质量.     

非连续变形分析方法的改进及应用

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现.其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征.研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳.改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

基于持续同调的边坡土颗粒应力链失稳分析

采用颗粒离散元方法和持续同调理论,研究了内排土场堆叠至不同高度时的边坡稳定性。为便于研究,现采用一水平金属板向下施加压力,代替不同厚度土层的重力荷载,对边坡在竖向荷载作用下的失稳破坏过程进行了颗粒离散元模拟。研究了二维边坡土颗粒速度总矢量、边坡失稳破坏时滑移开裂面的角度以及边坡坡顶y方向的平均速度等宏观响应过程,并构建了自然堆积下边坡堆积体颗粒的法向力链无向网络模型。最后,用持续同调方法对边坡坡顶颗粒接触力链网络的拓扑特征进行分析,获得条码图,建立了岩体结构持续同调特征与失稳演化的关系。本文为研究边坡失稳拓扑识别提供了一种新方法,从而可以有效预测边坡的失稳破坏。     

工程岩体地质力学作用过程的数值模拟

工程岩体地质力学作用过程的数值模拟, 是研究岩石工程设计、施工问题的有效途径。本文提出了工程岩体计算数据智能化准备与确定子系统、地质力学作用过程的仿真数值模拟及其快速图形处理子系统, 并将其集成为一体。系统采用人机交互界面及多系统的自动转换, 使整个分析系统能充分发挥其模拟分析及图形处理功能, 方便实用。     

刚性元方法和块状岩体稳定性分析

本文提出的刚性元方法是一种数值分析方法,它将变形体离散为一些刚性块体(剐性元)和块体之间的可变形薄层(节理元)的组合体,在各种载荷条件下求出块体的运动以及薄层的变形和应力分布。本文联合岩石系统的失稳准则来使用刚性元方法,成功地研究了块状岩体的稳定性问题。计算实例表明,本文提出的方法可以正确地反映块状岩体的失稳机制,是一种合理和可行的岩体稳定性分析手段。     

不连续岩体变形破坏过程中AE的空间演化与意义

通过一些典型节理化岩体模型变形破坏过程中的AE事件的定位分析和破裂面形态的研究,探讨了不连续岩体变形破坏过程的空间演化特征。结果表明：单节理岩体变形破坏的空间动态特征受节理角控制,随着节理角从小变大,岩体的变形破坏方式从AE在局部集中的突发式失稳向AE弥漫于整个节理层面的类似渐进式的破坏形式发展；失稳破裂面从简单的整体状岩-膏界面变成形态复杂的层内复合破裂面；复杂节理岩体的变形破坏受一条（组）主导构造控制,强AE事件集中发生在不连续构造交汇部位,在包含多条主构造的岩体模型中最终失稳破裂面只与最后阶段AE的位置一致。     

差异流动变形破坏研究

工程岩体出现许多流变现象, 它最突出的特点是差异流动变形破坏。本文从流变现象及变形量测结果论述差异流动变形破坏特征, 并且根据蠕变流动变形破坏与时间的关系进行了差异流动变形分类, 这些类型同时存在于一个岩体或岩块内。     

岩石的记忆功能及对岩体变形的影响

岩石的记忆功能包括残余应力对材料力学性质的影响, 应力的“冻结”特性和对固结应力的记忆。本文叙述这三类应力的理论基础、力学特性和计算方法, 通过一些工程算例指明其对岩体变形及岩体稳定性分析的重要性。     

岩石劈裂破坏过程的计算模拟

岩体的破坏主要是岩体中裂纹的扩展与贯通的结果.受压岩体的裂纹尖端,可能处于压剪状态,也可能处于拉剪状态.从断裂力学的角度,研究了地下工程围岩中的裂纹在各种受力条件下(包括拉剪、压剪两种情况)的起裂、扩展及贯通以至最终形成劈裂裂缝的特征及各个阶段的判据.针对压剪断裂,裂纹之间存在有粘结力与摩擦强度的特征,提出了岩石压剪断裂性判据.最后以受压岩体中的斜裂纹扩展过程分析为例,模拟了岩体中裂纹逐渐发展成劈裂破坏的过程,并说明了各个阶段的计算方法.