基于分段模拟的岩石蠕变模型

为了更好地实现对岩石蠕变特性的描述,基于岩石蠕变的阶段特征,提出分段模拟岩石蠕变的建模思路。应用变分数阶微积分改进的软体元件与弹簧串联得到改进的Maxwell模型实现对岩石稳定蠕变阶段的模拟;同时基于应变触发方式改进蠕变体元件实现对岩石非稳定蠕变阶段的模拟。综合改进的Maxwell模型与改进的蠕变体元件构建一个新的三元件蠕变模型。应用岩盐典型蠕变曲线数据对三元件蠕变模型进行了验证。结果显示构建的三元件蠕变模型可以很好地模拟岩盐的三阶段蠕变,分段模拟思路是可靠的。     

单轴压缩下绿砂岩长期强度的尺寸效应研究

岩石的蠕变特性是影响岩体工程稳定性的重要因素,而岩石的长期强度是确定岩体工程长期稳定的一个重要指标。由于岩石材料的非均质性,其长期强度具有明显的尺寸效应。为了研究岩石长期强度的尺寸效应,首先,在幂函数模型基础上,基于损伤力学理论,建立了能够描述岩石蠕变全过程的非线性蠕变损伤模型;然后,把运用该模型计算得到的单轴压缩蠕变数值模拟结果与室内单轴压缩蠕变试验结果进行对比,验证了模型的正确性;最后,采用所提出的模型对7个不同尺寸的岩样进行了单轴压缩蠕变数值模拟,并对岩石长期强度尺寸效应进行了分析。数值模拟结果表明：随着试样尺寸的逐渐增大,岩石长期强度值逐渐减小,当试样尺寸增大到一定程度时,岩石长期强度稳定在一个特定值附近。     

基于西原模型的变参数蠕变方程

岩石材料的蠕变具有明显的非线性特征,在蠕变模型中考虑蠕变参数随应力、时间等因素而变化是解决非线性蠕变问题的方法之一。试验已经证实,应力、时间及渗流等外部条件会引起岩石某些力学参数的改变,因此,将蠕变参数视为变量符合岩石材料的本质属性。以西原蠕变模型为基础,研究了各蠕变参数与应力及时间的关系。根据试验结果,通过非线性最小二乘法从蠕变试验数据中求出不同应力作用下不同时刻的各个蠕变参数,再根据参数随应力、时间的分布规律经非线性拟合得到各蠕变参数与应力及时间的表达式,进而可得到变参数的蠕变方程。蠕变参数的变化是岩石非线性特征的充分体现,从损伤力学的角度来讲,蠕变参数的变化反映了岩石内部损伤不断加剧、累积,材料特性不断劣化的演变过程。     

润扬大桥北锚碇岩石流变特性的试验研究

为了解润扬长江大桥北锚基础岩石的流变力学特性,针对北锚基础区域的岩样,通过模拟其受力的实际情况,采用岩石双轴流变试验机进行了单轴压缩蠕变试验,并选用广义的Kelvin模型进行了参数拟合分析,掌握了岩石蠕变的基本规律,并获得了粘滞系数等相关的流变力学参数。     

软弱破碎围岩隧道炭质页岩蠕变特性试验研究

采集贵广铁路天平山隧道典型软弱破碎围岩岩样,在室内开展三轴常规压缩和三轴压缩蠕变试验。基于岩石试验成果,分析炭质页岩的蠕变力学特性,低应力水平条件下岩石仅出现减速蠕变和等速蠕变阶段,在更高的应力水平条件下,岩石并未出现加速蠕变的现象,但围压在一定程度上提高岩石的屈服强度,弹性模量亦有随着围压增加而增大的趋势。据此,建立依托工程炭质页岩蠕变全过程的黏弹塑性应变软化蠕变力学模型,推导该蠕变力学模型在三维应力状态下的本构方程,并确定模型的参数,运用最小二乘法对模型的参数进行辨识。通过与试验曲线对比显示,所建蠕变模型能够描述依托工程炭质页岩的蠕变力学特征,可以用来研究该隧道软弱破碎围岩的蠕变性质和稳定性。     

岩石蠕变声发射特性研究现状与展望

蠕变是岩石一种重要的力学特性,与工程的长期安全稳定和安全密切相关。声发射技术作为一种无损检测手段,被广泛应用于岩石变形破坏研究领域。本文阐述了近年来基于声发射技术开展岩石蠕变特性方面的研究进展,主要包括不同岩性、加载方式、应力水平下岩石蠕变声发射特性,岩石蠕变破坏声发射前兆特征。在此基础上,提出了几点今后需要进一步深入研究的问题。     

甘肃甘谷裂隙泥岩剪切蠕变行为及其修正模型研究

岩石在长期荷载作用下表现出典型的蠕变特性,裂隙泥岩在长期荷载作用下的力学性质不同于完整岩石。对甘肃甘谷赵家窑滑坡钻孔取芯的裂隙泥岩开展剪切蠕变试验,探讨轴向压应力、预制裂隙长度等因素对试样剪切蠕变特征的影响。试验结果表明:轴向荷载和裂隙长度对试样剪切蠕变特征有明显的影响。轴向压应力越大,试样剪切蠕变变形越不明显;裂隙越长,试样剪切蠕变变形越显著。应力较低时,如Cvisc模型等蠕变模型能够较好反映泥岩的蠕变力学特性。而在加速蠕变阶段,大多数模型却不能准确体现泥岩蠕变曲线的非线性变化。改进的西原模型克服了传统模型不能对非线性蠕变特征进行描述的缺点,且模型参数较少,能较好地适用工程软弱岩体。不同载荷作用下泥质软岩试样剪切蠕变过程中,裂隙泥岩剪切破坏带形成表现为渐进与突发相结合。泥岩的剪切蠕变力学模型可解释泥岩滑坡的多种现象,边坡潜在滑移带内泥质软岩发生的剪切蠕变则是泥岩滑坡的内在诱因。     

单轴压缩下岩石蠕变失稳破坏过程数值模拟

在岩石破裂过程分析（RFPA2D）系统的基础上,考虑岩石损伤过程的时间因素影响,引入岩石细观单元蠕变本构方程,建立了考虑流变效应的岩石破裂过程RFPA2D数值模型。应用该模型模拟了恒定荷载作用下岩石的蠕变破坏过程,得到了岩石蠕变破裂的3个典型阶段：初始蠕变阶段、稳定蠕变阶段和加速蠕变阶段,模拟结果同实验室试验所观察到的现象十分吻合。这表明考虑流变效应的RFPA2D数值模型适用于模拟岩石的蠕变破坏这一复杂的、非线性演化问题。此外,数值模拟还揭示了岩石的宏观蠕变破坏实质上是细观层次上单元损伤累计的结果,这些结论对岩石工程的长期稳定性研究具有重要的理论指导和实践意义。     

三轴压缩下粉砂质泥岩蠕变本构模型研究

依据粉砂质泥岩三轴压缩蠕变试验结果,分析了岩石的蠕变力学特性。选取线性黏弹性Burgers模型来描述岩石的衰减蠕变和稳定蠕变特性。引入非线性黏塑性元件,将其与Burgers模型串联起来,建立一个新的六元件非线性黏弹塑性Burgers蠕变本构模型。采用非线性Burgers模型来描述岩石的加速蠕变特性。应用Levenberg-Marquardt算法对蠕变试验曲线进行辨识,得出各级应力水平下岩石轴向与径向的蠕变模型参数。对比试验曲线和模型拟合曲线,二者吻合较好,表明线性Burgers模型与非线性Burgers模型可以准确描述岩石3个蠕变阶段的力学特性。模型参数辨识结果表明,在同一级应力水平下岩石轴向与径向的蠕变模型参数各自独立并不统一,模型参数取值时应考虑岩石蠕变的各向异性以及轴向与径向蠕变量、蠕变速率大小的差异。在前7级应力水平下,模型参数应取以轴向应变辨识得到的Burgers模型参数;在第8级应力水平下,模型参数应取以径向应变辨识得到的Burgers模型参数;在破裂应力水平下,模型参数应取以径向应变辨识得到的非线性Burgers模型参数。     

岩石蠕变的应力-应变比分析

对岩石蠕变的应力-应变比分析方法进行了探索。蠕变是岩石重要的力学性质之一。岩石在各个应力水平下的蠕变曲线具有几何相似性。根据岩石蠕变的应力-应变比曲线可以近一步分析得到等时模量曲线。在岩石蠕变本构描述十分复杂、至今没有成熟的理论的情况下,应力-应变比分析方法具有十分重要的现实意义。     

基于FLAC3D的改进Burgers蠕变损伤模型的二次开发研究

岩石的流变变形对于岩质边坡和地下结构产生重要影响,其时效特性越来越引起人们的关注。然而通用分析软件所提供的岩石流变本构模型均有其特定的使用范围,往往不能满足工程数值分析的需要。考虑到岩体中参数劣化造成的影响,采用了一个变参数的改进Burgers蠕变损伤模型,该模型认为岩体的蠕变参数是随时间因子逐渐弱化的,这弥补了FLAC3D自身流变模型的不足,并利用FLAC3D（Version 2.1）提供的二次开发程序接口,分析了蠕变损伤模型二次开发程序运行的基本原理,给出了该模型具体的程序框图和代码编写中应该注意的关键问题。介绍了该模型的几种退化形式,通过算例分别验证了该模型黏弹性、塑性和损伤的力学性质,分析了损伤参数对蠕变曲线的影响,并由此介绍了验证流变模型二次开发程序正确性的一般步骤。由于采用了面向对象的编程技术,FLAC3D二次开发接口简单实用,研究思路可为其他流变本构模型的二次开发提供参考。     

Effect of simulated sampling disturbance on creep behaviour of rock salt

Summary This article presents the results of an experimental study of creep behaviour of a rock salt under uniaxial compression as a function of prestrain, simulating sampling disturbance. The prestrain was produced by radial compressive loading of the specimens prior to creep testing. The tests were conducted on an artifical salt to avoid excessive scattering of the results.The results obtained from several series of single-stage creep tests show that, at short-term, the creep response of salt is strongly affected by the preloading history of samples. The nature of this effect depends upon the intensity of radial compressive preloading, and its magnitude is a function of the creep stress level. The effect, however, decreases with increasing plastic deformation, indicating that large creep strains may eventually lead to a complete loss of preloading memory.     

考虑初始损伤和蠕变损伤的岩石蠕变全过程本构模型

考虑到天然岩石存在不同程度初始损伤以及蠕变过程中岩石受载后裂隙扩展而导致的新损伤,对具有初始损伤的岩石蠕变特性进行全面描述。根据不闭合结构面应力与法向变形之间的关系,提出裂隙岩石塑性变形体元件,描述岩石蠕变过程中的瞬时塑性变形。引入初始损伤影响因子,建立具有初始损伤的岩石损伤变量演化方程,构建模拟岩石加速蠕变的蠕变损伤体元件。将裂隙岩石塑性变形体和蠕变损伤体与描述瞬时弹性变形和黏弹性变形的广义开尔文模型进行串联组合,形成能够反映具有初始损伤的岩石瞬时弹-塑性变形、稳定蠕变和加速蠕变的蠕变全过程本构模型,提出了进行少量蠕变试验既能解析模型参数的方法,在不同应力水平下模型理论曲线与蠕变试验曲线吻合。     

一种模拟岩石蠕变的数值流形方法

为了保障岩土工程结构能长期正常使用,需要对其蠕变变形进行分析。“时步-初应变”法是一种常用的计算岩石蠕变的方法。数值流形方法是一种新兴的数值计算方法,常用于计算节理岩体的变形,但尚未被试用于计算蠕变变形。在原数值流形方法的程序中增加了基于“时步-初应变”法的计算模块,通过对广义开尔文模型进行的模拟,显示新程序可以正确反映岩石的黏弹性蠕变趋势,并能够计算包含节理的岩体的蠕变变形,改进后的数值流形方法不但能够模拟岩石的线弹性变形,而且可以模拟岩石的黏弹性蠕变,比原流形方法更能全面地模拟岩石的变形,扩展了数值流形方法在岩土工程中的使用范围。     

基于分数阶导数的非定常蠕变本构模型研究

针对锦屏二级水电站深埋长大隧洞围岩在开挖过程中呈现明显的流变劣化特征,基于分数阶微积分理论,用Abel黏壶替换西原模型中的牛顿黏壶,推导出基于分数阶导数的蠕变本构方程,同时考虑岩石在蠕变过程中,尤其是所受应力超过其长期强度时,蠕变参数是非定常的特征,把蠕变参数的非定常性引入本构方程中,建立了基于分数阶导数的非定常蠕变本构模型。结合锦屏二级水电站引水隧洞主要围岩之一的大理岩的室内剪切蠕变试验的结果,对模型中的参数进行反演,结果表明,基于分数阶导数的非定常蠕变本构模型能够在蠕变初期和拐点处与试验结果很好地吻合,并且可以克服西原模型不能反映蠕变第3阶段的不足。通过对模型中参数的敏感性进行分析,得出分数阶阶数、非定常参数对岩石蠕变的影响规律。研究结果表明,提出的基于分数阶导数的非定常蠕变模型能够很好地反映大理岩的蠕变全过程。     

Enhanced double-mechanism creep laws for salt rocks

The double-mechanism creep law (DM model) is an empirical creep constitutive model widely employed in Brazilian salt rock mechanics. This model often presents good performance in steady-state creep prediction. However, transient creep is not accounted for, and whenever early creep estimates are important, the contribution of this phase might be meaningful. This work adds value by presenting two alternatives to account for transient creep in the DM model. The first alternative couples the transient function from Sandia’s multi-mechanism deformation model to the DM model steady-state creep rate (EDMT model). The second alternative couples the DM model response to Norton’s power law when the strain rate given by the latter remains lower than the one from the former (EDMP model). These models can be implemented in numerical simulators at small code extensions of the DM model implementations. Applications from previous works by the authors are revisited to validate the formulations based on experimental data. EDMT and EDMP models differ in the formulation of transient creep and, consequently, in the time of transition between the transient and the steady-state phases. Both methods were successful in treating transient creep and in simulating experimental results.     

锦屏大理岩蠕变损伤演化细观力学特征的数值模拟研究

岩石的蠕变损伤和断裂是岩石流变效应的重要表现形式,但其损伤演化过程往往难以直观观测,为此,采用二维颗粒流数值模拟方法(PFC2D)对岩石的蠕变损伤和断裂的细观力学机制进行了分析。在锦屏大理岩室内试验基础上,利用颗粒流应力腐蚀模型(PSC),建立了能反映其短期和长期强度特征的柱状岩样数值模型,并开展了大量数值试验。结果表明,蠕变损伤的演化过程与暂态的损伤演化过程具有明显的差异。在岩样蠕变损伤过程中,其内部微裂纹多沿加载方向开裂且分布均匀,先快速增加再稳定扩展,最后则发生快速断裂。当荷载较小时,岩样宏观上呈现劈裂破坏特征,当荷载较大时,岩样呈现剪切破坏特征。在岩样蠕变损伤初始和稳定演化阶段的前期,荷载大小对岩样的损伤演化过程影响不大;在稳定演化阶段的后期至断裂过程中,低荷载下岩样的损伤增速比高荷载下快。