岩石三维破裂过程的数值模拟研究

采用细观弹性损伤模型和有限元计算方法实现岩石三维破裂过程的数值模拟。考虑到岩石非均匀性的本质特征，通过引入简单直观的单元本构模型，采用细观单元材料性质退化的办法，利用位移加载来实现岩石逐渐破裂过干旱；模拟单轴压缩、单轴拉伸和剪切破裂3种基本试验，得到岩石非线性应力-应变曲线和不同载荷阶段三维损伤破裂演化系列图像；分析细观非均匀性对岩石宏观破裂力学行为的影响。试验研究表明，三维破裂比二维破裂更为复杂，RFPA^3D可以有效地模拟脆性材料的三维破裂。     

基于细观结构表征的岩石破裂热–力耦合模型及应用

 岩石热破裂的研究只有考虑各种矿物组分造成的岩石的非均匀性,才能更客观地反映岩石热破裂的本质。利用数字图像处理技术数字化表征岩石内部矿物颗粒的几何形态,充分考虑岩石真实的细观结构,结合细观损伤力学和热弹性理论,建立能更客观的分析岩石热–力耦合作用下破裂过程的数值模型。以花岗岩为例,运用数值模型研究花岗岩在温度和压缩荷载共同作用下的力学行为和破裂过程。研究结果表明,温度对岩石的力学性质和破裂演化过程影响显著,热破裂裂纹多发生在矿物颗粒边界处,并沿颗粒边界扩展,局部会形成闭合多边形,其热破裂演化过程与试验结果基本相符,从而验证了数值模型的合理性和有效性,该数值模型为细观尺度定量研究岩石热破裂提供一种新的方法。     

岩石破裂过程分析系统并行计算方法研究

岩石工程灾害与岩石破裂过程失稳密切相关。大型岩石工程破裂过程数值分析需要高效、准确、强大的计算能力支持。一般传统串行计算方法难以满足要求，大规模并行计算是解决这一难题的有效途径。岩石破裂过程分析系统是研究岩石破裂过程的一个重要数值分析工具。代有限元方法和数值计算方法，在消息传递并行环境下，在岩石破裂过程分析系统串行单机版的基础上，结合现利用区域分解和主从编程模式，采用分布存储稀疏线性迭代并行求解方法，在Linux机群上实现应力分析模块中有限元计算的并行处理。通过Windows和Linux协调处理策略，有效地把原有的前后处理功能和机群系统强大的计算能力结合起来，建立岩石破裂过程分析RFPA^3D．Parallel并行分析系统。算例结果表明，并行程序具有很高的加速比和并行效率，能够快速完成三维条件下300万单元的大规模岩石破裂过程分析。应用RFPA^3D-Parallel并行分析系统模拟地壳介质中广泛存在的龟裂现象，再现非均匀介质破坏和裂纹演化过程，从而显示该系统广泛的应用前景。     

热应力作用下的岩石破裂过程分析

热应力引起的岩石破裂称为岩石的热破裂，它是热和力之间相互耦合作用的结果。岩石热破裂研究的工程意义重大。根据岩体介质变形及其热力学的理论基础，充分考虑岩石的非均匀性和热固耦合作用，在原有的岩石破裂过程分析系统的基础上，建立了具有热固耦合作用的岩石热破裂分析模型。数值模型再现岩石的热破裂过程，并反映岩石热破裂的规律。运用数值模型，对含有单个内嵌颗粒的岩石试件在温度变化过程中的热开裂进行了数值模拟。研究结果表明：在温度升高过程中，如果内嵌颗粒的热膨胀系数大于基质的热膨胀系数，在基质内产生径向裂纹；如果内嵌颗粒的热膨胀系数小于基质热膨胀系数，便在基质内产生环向裂纹。数值模拟结果与试验结果有较好的一致性。RFPA^2D-thermal模型为从细观力学角度上分析岩石的热破裂过程和机制提供了一种新的方法。     

大理岩破坏过程的三维细观弹塑性损伤模拟研究

为了考虑岩石细观局部塑性变形,基于应变空间理论导出细观弹塑性损伤模型,采用有限元计算方法实现岩石三维破裂过程的数值模拟。提出细观破坏单元网格消去法,利用位移加载来实现岩石逐渐破裂过程;模拟二维、三维大理岩三点弯曲梁弹塑性损伤破坏试验,得到岩石非线性应力–应变曲线和不同载荷阶段弹塑性损伤破裂演化系列图像;分析细观非均匀性对岩石宏观破裂力学行为的影响。试验研究表明,三维破裂比二维破裂更为复杂,本文所建模型可以有效地模拟岩体的三维破坏过程。     

水力压裂时岩石破裂压力数值计算

 以多孔介质流体渗流与岩石应力–应变耦合理论为基础,推导水力压裂时岩石破裂压力数值计算方程,提出一种全新的岩石破裂压力计算方法。该方法考虑岩石变形与流体渗流的全耦合作用,采用有限元法数值计算技术,能模拟计算流体向地层渗滤情况下井眼周围地层有效应力的时空分布,因此,该方法不但能求得岩石的破裂压力,还能同时求得岩石的破裂时间。建立的破裂压力数值计算理论克服传统破裂压力解析计算方法的许多不足,如无法精确计算井眼周围孔隙压力的升高导致应力集中加剧对破裂压力的影响,无法计算地层破裂的确切时间等问题。该计算理论的建立,实现岩石破裂压力数值计算的突破和计算精度的提高,为水力压裂时岩石破裂压力的计算找到了新的理论和方法。     

关于计算岩石力学发展的几点思考

关于计算岩石力学发展的几点思考王泳嘉冯夏庭（东北大学资源与土木工程学院，沈阳，１１０００６）计算岩石力学是指用计算机对复杂的岩石力学问题进行数值计算或模拟试验，也就是说在计算机上模拟岩石力学的物理过程，研究其规律性，核实试验数据的可靠性，检验理论预...     

岩石剪切破坏过程的RFPA2D数值模拟

运用岩石破裂与失稳过程分析RFPA^2D系统，对岩石剪切破裂过程进行了数值模拟研究。将剪切岩石视为非均匀弹－脆岩石材料，模拟结果再现了岩石剪切滑动形成从变形到破坏直至失稳的全过程，剪断面是滑动形成的主要形式，并首先在试样一端出现，然后再形成由一端及里的剪断面扩展直至另一端最后产生剪断面破坏贯通，形成统一的滑动面，岩石剪切破裂面分形维数还对应着剪断面的粗糙程度和力学行为。     

岩石破裂过程分析(RFPA2D)系统的细观单元本构关系及验证

从岩石的细观非均匀性特点出发，应用弹性损伤力学，对岩石破裂过程分析数值（RFPA^2D)模拟系统的本构关系进行了描述。用RFPA^2D系统模拟了岩石试样在单轴或双轴载荷作用下的强度特征和破坏形态，通过与实验结果对比验证了该本构关系和数值模拟系统的合理性和有效性。     

非均质岩石材料宏细观力学参数的关系研究

利用基于岩石矿物种类及其含量提出的岩石矿物细胞元随机性参数赋值方法,针对两矿物细胞元混合的非均质岩石模型进行参数赋值,通过数值试验研究了非均质岩石材料矿物细胞元含量及其力学参数对岩石宏观力学参数的影响;根据数值试验结果,探讨了岩石材料矿物细胞元弹性模量与宏观模量之间的函数关系。研究结果表明,岩石宏观弹性模量与矿物细胞元弹性模量、矿物细胞元含量都具有十分显著的线性相关关系,即基于岩石矿物种类及其含量的矿物细胞元随机性参数赋值方法建立的岩石有限元模型宏细观力学参数之间呈显著的线性关系;宏观力学参数可用细观力学参数的数学期望值即加权平均值来表示。