岩石破裂与渗流耦合过程细观力学模型

从岩石细观非均匀性的特点出发,应用弹性损伤力学、Biot渗流力学理论,对岩石破裂过程渗流-应力-损伤耦合模型进行了描述.算例表明,该模型能较好的模拟出岩石材料在水压和载荷作用下裂纹扩展路径不规则发展及稳态、瞬态渗流过程,裂纹扩展对渗流路径和流动过程起到明显的控制作用.     

岩石破裂过程THMD耦合数值模型研究

从岩石的细观非均匀性特点出发,应用损伤力学、热力学和渗流力学理论,建立了岩体热(温度)-水(渗流)-岩(应力)-损伤耦合数值模型(THMD model),把岩石(体)THM耦合问题的研究从应力状态分析深入到损伤、破坏过程分析之中。探讨了THM耦合作用下岩石材料的细观结构损伤及其诱发的材料力学性能演化机制,并运用所提出方法计算温度-渗流-应力耦合作用下井筒近场围岩的稳定性,模拟得到的岩体破坏过程、应力分布、AE特性及渗流特性变化与现场标定结果有着一致的规律性,初步证明了该数值模型的合理性和有效性。THMD模型以简单的数值模型表征了岩石(体)中热、水、岩及损伤之间复杂的作用关系,为从细观损伤演化揭示宏观岩体温度-渗流-应力耦合破坏机制提供了一种新的数值分析方法。     

岩体中饱和渗流应力耦合模型研究进展

在大坝渗流与控制、各种水工隧洞和交通隧道以及地下各种洞室开挖、采矿工程与油藏开采中的水力劈裂、 岩坡和坝基的稳定性研究等领域中经常会遇到岩体的饱和渗流应力耦合问题.本文对这些领域中的渗流应力耦合问题分为6种力学模型: 等效连续介质模型、 裂隙网络模型、 双重介质模型、 断裂力学模型、 连续损伤力学模型和统计模型.前3种主要从经典的黏弹塑性本构着手考虑, 着重对渗流场的处理, 后3种则是从考虑岩体在耦合作用下发生损伤破裂行为方面切入, 着重于岩体结构内部发生质的改变带来的更加复杂的耦合效应.对这些模型进行了详细的介绍, 指出了这些模型的优缺点与适用范围, 给出了一些有代表性的研究成果, 并对该研究领域中未来的研究方向进行了展望.     

含损伤演化的TM耦合数值模型及其应用研究

基于生物力学中的Wolff法则，发展了一种连续体拓扑优化的新方法. 有待优化的结构 被看作是一块遵从Wolff法则生长的骨骼，把寻找其最优拓扑的过程比拟为骨骼的重建/生 长过程. 采用骨骼的重建/生长规律作为准则更新材料分布，直至达到一个平衡状态并由此 获得结构的最优拓扑. 算例表明了所提出方法的有效性.     

考虑软硬物质双变形特征的脆性岩石损伤本构模型研究

基于脆性岩石变形力学特征,将脆性岩石视为由岩石颗粒骨架等硬物质以及裂隙和孔隙等软物质组成;考虑软硬物质变形特点,分别建立了软物质与硬物质的本构关系以及软硬物质与RVE之间的关系。然后,基于几何损伤理论,建立了考虑软硬物质双变形特征的脆性岩石变形破坏模拟方法,根据软硬物质及其与脆性岩石宏观变形关系以及峰值点法,分别提出了力学参数和统计参数的合理确定方法。最后,通过实例分析对本文模型和方法进行验证,结果表明本文模型能够对脆性岩石变形破坏过程进行模拟,脆性岩石宏观变形被视为由软物质变形和硬物质变形两部分组成具有可行性与合理性。     

有限应变下蠕变—弹塑性—损伤耦合本构模型的数值模拟

本文提出了一个模拟有限元应变下各向同性蠕变－弹塑性－损伤耦合本构行为的数值模型。模型基于不可逆过程的热动力学和内状态变量理论、以及应力就变本构关系的超弹性形式描述和变形梯度的弹性及非弹性部分的乘式分解。     

单轴拉伸条件下脆性岩石微裂纹损伤模型研究

利用断裂力学、损伤力学和均匀化原理,对脆性岩石单轴拉伸条件下的力学特性进行分析,建立了脆性岩石的微裂纹损伤本构模型.首先对岩石内部微裂纹的统计分布规律进行分析,给出了理论分析过程中微裂纹分布的假设条件,在此基础上,参考已有研究成果,得到含细长微裂纹脆性岩石有效弹性参数的计算公式.然后,对岩石内部单一微裂纹进行断裂力学和损伤力学分析,得到了扩展裂纹尖端的应力强度因子计算公式,在一定微裂纹断裂扩展准则和断裂扩展速率的假设基础上,利用积分原理,得到了岩石整体的损伤变量和损伤演化方程,由此建立单轴拉伸条件下脆性岩石的微裂纹损伤本构模型.最后,通过一花岗岩的单轴拉伸试验结果对微裂纹损伤本构模型进行了验证.     

岩体水力劈裂的应力-渗流-损伤耦合模型研究

利用无单元法追踪裂纹扩展的优势研究模拟岩体水力劈裂的数值分析模型。定义损伤变量以描述岩体的损伤程度对岩体渗透性和强度的影响,建立了岩体水力劈裂的应力-渗流-损伤耦合分析模型,并编制了基于无单元法的计算程序。综合考虑应力、渗流及损伤之间的相互作用影响,分析了具有初始裂纹的岩体平面应力模型的裂纹扩展过程,指出考虑渗流作用时的裂纹扩展角大于不考虑渗流作用时的裂纹扩展角。同时,裂纹内水压力和渗流的作用对于裂纹的扩展方向和扩展长度具有较大的影响。     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。     

动力损伤后的脆性岩石静力蠕变断裂模型研究

应力波动力扰动下脆性岩石的静力蠕变特性,对深部地下工程围岩变形的评价有重要的实践意义.动力载荷作用导致的局部细观裂纹损伤严重影响脆性岩石蠕变力学行为.基于细观裂纹扩展与应力关系模型、动力扰动损伤演化函数、静动力载荷演化路径函数与黏弹性本构模型,提出一种应力波动力扰动下脆性岩石蠕变断裂特性的宏细观力学模型.其中动力损伤通过控制岩石内部细观裂纹数量变化实现.模型描述了应力波动力扰动下岩石的应变时间演化曲线,解释了岩石动力扰动下蠕变失效特性.研究了不同应力波幅值及周期影响下的脆性岩石应变-时间关系曲线,并通过试验结果验证了模型的合理性.讨论了动力损伤变化形式,突变发生时刻,突变量的大小对岩石蠕变失效特性的影响.分析了应力波幅值、周期对岩石动态动力损伤效应以及蠕变失效特性的影响.主要研究结果:动力损伤的变化值越大,岩石蠕变失效发生时间越短.冲击载荷扰动期间,动力损伤发生的时刻及增加的形式,对动力扰动后的岩石应变及蠕变破坏时间影响很小.动力损伤变化量随应力波幅值增加、周期减小而加速增大.应力波幅值越大、周期越小,岩石发生蠕变失效时间越短.     

水岩耦合变形破坏过程及机理研究进展

坝址水力劈裂、煤矿和隧道突水以及水力压裂工程都是岩体损伤演化积累诱发破坏 (灾变) 的过程,然而能否正确认识和评价这些过程取决于人们对水岩耦合变形破坏过程和机理的研究程度. 阐述了岩石渗流 --变形 --破坏特性及机理研究的进程, 侧重于物理实验成果. 综述了水岩耦合变形破坏过程分析模型和数值计算, 特别是有限单元法处理裂纹的思路, 并介绍了作者课题组近年来所取得的一些代表性成果. 给出了当前理论 (数值) 研究中所面临的几个具体问题和挑战, 并对今后研究工作进行了力所能及的展望.      

岩石声发射Kaiser效应的数值模拟试验研究

运用岩石破裂过程分析RFPA^2D系统,对岩石破裂过程中声发射Kaiser效应进行了数值模拟试验研究,得到了与物理试验现象十分吻合的模拟结果。     

模拟岩石材料脆性破裂过程的三维离散元模型

发展一种能够模拟岩石材料脆性破裂过程的三维不规则、可变形块体离散元模型. 一方面,在裂纹扩展过程中动态地将潜在破坏的连续块体沿潜在破坏方向细化为若干子块体,并在子块体之间的界面上设置连接型弹簧;另一方面,连接型弹簧在满足张拉-剪切复合破坏准则时发生脆性破裂,转变为接触型弹簧,实现材料由连续到非连续的破裂. 借助动态松弛技术完成求解,通过计算实例验证该方法的适用性.     

黄土损伤与流变耦合模型及参数研究

以试验为基础,对不同含水量Q3黄土的蠕变特征进行了分析研究,结合黄土结构损伤特性及其对黄土力学特征影响,通过考虑结构损伤对变形模量的劣化来反映黄土结构性对蠕变规律的影响,提出了黄土损伤与流变耦合作用的本构模型。该模型不仅具有参数少的优点,而且能够准确描述黄土的突发性破坏特征。利用结构性软土流变的研究成果,给出了黄土流变模型参数的确定方法,并对模型参数随非饱和黄土含水量变化的规律进行了研究。数值模拟结果表明,该模型能够很好地描述黄土的加速蠕变阶段的变形特征。     

混凝土三点弯曲试件破坏过程的数值模拟

根据混凝土试件三点弯曲试验的物理模型,用MFPA数值模拟软件分别对砂浆和混凝土的三点弯曲试件的破坏全过程进行模拟．给出两个试件的载荷与加载点位移关系曲线,并进行了对比分析．     

岩样单轴压缩变形破坏与能量特征研究

基于在伺服试验机上得到的同直径不同长度大理岩样的单轴压缩试验结果,研究了尺寸对岩石变形破坏与能量特征的影响规律.结果表明,岩石长度对峰值应力前的轴向变形特性没有显著影响,但明显改变峰后的轴向变形特性;岩样长度对峰后塑性变形特性影响不大;直径一定时随着长度的增大,岩样破裂形式由竖向劈裂变为剪切破坏.随着长径比的增加,岩石破坏应变能呈幂律衰减趋势,且当长径比增加至2以后,便基本趋于恒定,不再明显减小;而岩石破坏能呈先减小后增大的非线性变化趋势.     

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

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

岩石爆破分形损伤模型研究

在Ｋｕｓ和Ｔｈｏｒｎｅ等人研究基础上,利用分形理论及应用成果提出了一种新的岩石爆破损伤模型。新模型将岩体中各种结构弱面视为初始损伤,并将其分形维数作为岩石性质的主要参量,用损伤能量耗散率及与分形维数的关系表征爆破过程中的岩石损伤演化规律,从而弥补了现有模型未考虑初始损伤影响和将损伤演化归结为体积应变函数的不足。还通过模拟试验等手段验证了新模型的如上观点。     

三维非均匀脆性材料破坏过程的数值模拟

采用有限元方法模拟了三维均匀固体材料在宏观等效力学性质和破坏过程。首先采用格形（lattice）方法把试件离散成三维均匀网格,在每个单元格中将材料按照均匀处理,根据给定的统计规则来确定不同单元格中的材料常数以反映材料的非均匀性。然后对非均匀脆性材料选用简单的本构关系与断裂准则,采用自适应选取载荷步长对试件进行加载,通过非平衡迭代技术对刚度矩阵进行不断修正,实现了非均匀脆性材料的弹性行为及破坏过程的数值模拟。在此基础上,通过数值计算研究了材料的非均匀分布对宏观等效力学性质和破坏过程的影响,给出了破坏全过程的非线性载荷－位移曲线以及不同载荷阶段的三维损伤破坏的演化图。