New Trends in Rock Mechanics

Some new results in rock mechanics and some applications are presented. First, it is shown how an elastic/viscoplastic nonassociated constitutive equation can be derived for rock-like and particulate materials. This constitutive equation is based on the fundamental concept of compressibility/dilatancy. Damage and microcracking energy are discussed. Examples are given to demonstrate that it is possible to predict where the rock near underground openings is dilatable and where is compressible, where and when short-term or long-term failure is to be expected, the amount of rock involved in evolutive damage, etc.     

A procedure to determine nonassociated constitutive equations for geomaterials

A procedure for determining a phenomenological elastic/viscoplastic nonassociated constitutive equation for geomaterials is presented. For this purpose, triaxial test data obtained with either a “true” or a classical triaxial device are necessary. The constitutive equation is aimed at describing such geomaterial properties as creep, irreversible compressibility or dilatancy, work-hardening, damage, and failure. Long-term failure can also be described with this model. According to the procedure, first the elastic parameters are determined from unloading tests (which follow short-term creep tests), then the yield function is determined, and finally the viscoplastic potential. No a priori assumption is made concerning the form of the yield function or of the viscoplastic potential; their expressions are obtained from the data by using the procedure suggested here. Examples for sand and rock salt are given. Comparisons of the model predictions with the experimental data are discussed.     

不同应力差下非线性蠕变研究

深部工程处于复杂的环境场(应力场、温度场和渗流场)之中,引起的地下工程围岩流变破坏事故屡见不鲜。本文通过三轴蠕变实验,根据岩体流变力学和蠕变原理对影响深部岩体蠕变的因素进行了研究。用实验数据进行数值计算,得出了深部软岩巷道围岩在高地应力作用下的非线性蠕变模型,建立了应力差-时间-蠕变的定量关系,利用MATLAB软件得到了不同情况下的蠕变-时间曲线;应力差-时间-蠕变曲面。从而验证了非线性流变理论,找到了引起蠕变的主要因素应力差。并利用应力差理论分析了深部巷道围岩的蠕变破坏机理。     

基于改进西原模型巷道围岩流变参数反演分析

当应力达到一定水平时,岩石经过衰减和稳定蠕变过程之后会发生加速蠕变破坏,但是传统西原体模型无法更有效地反映该加速蠕变阶段。为了解决这一问题,引入一个非线性流变元件,并将其串联在西原体模型上,组成一个新的六元件模型,该模型可充分反映岩石蠕变的三个阶段。基于改进的西原体模型本构方程,应用Laplace变换,得到了盾构圆形巷道围岩的径向位移变化规律;根据巷道围岩收敛变形现场实测数据,反演得到了岩石黏弹塑性蠕变参数;由巷道围岩径向位移计算结果和实测数据对比分析表明:两者吻合较好;当t ≤25 d时,巷道围岩处于加速变形阶段,其径向位移变化速率随时间的增长而不断增加,当t>25 d时,巷道围岩进入稳定变形阶段,随时间的推移,其径向围岩位移变化速率逐渐趋于稳定。     

Viscoplastic creep of rocks around a lined tunnel

A rock-support analysis is done assuming that the rock behaves elastic/viscoplastically while the support can be any kind of nonlinear support. The case of circular tunnels is considered and the formulation of the mathematical problem is discussed. Furthermore, the creep of rock around a tunnel with nonlinear supports (called “yieldable” or “self-adapting” supports) to limit the closure of the tunnel is studied. It is shown that the ultimate ground reaction is not unique since it depends on the loading history such as excavation layout and support design. For relatively small pressures exerted by the support, it is shown that the onset of failure by dilatancy can be determined. Two kinds of solutions are given: a simplified approximate one, which is easy to apply and with which the qualitative discussion of the solution is simple and revealing, and a more rigorous general elaborate numerical solution obtained by using computer programs.     

Effect of time dependent compressibility on non-linear viscoelastic wave propagation

The work presented consists essentially of two parts: the first deals with the development of a non-linear constitutive equation for a three-dimensional viscoelastic material with instantaneous and time dependent compressibility; the second deals with the solution of some specific wave propagation problems for three simple three-dimensional geometries. The constitutive equation is based on the existence of elastic and creep potentials and is expressed in terms of single memory integrals with non-linear kernels. The wave propagation problems are solved by numerical integration along the characteristics of the governing equations. The primary conclusion drawn deals with the effect of time dependent compressibility on the dynamic stress, strain and velocity fields. Results indicate that the dynamic response of even slightly time dependent compressible materials varies dramatically from those assumed to have only an instantaneous elastic compressibility.     

2021-07期目录

在陈宗基院士关于地下硐室长期稳定性力学问题研究的基础上,采用Sadovsky院士关于复杂地质岩体的等级构造学说,围绕深部岩体非均匀构造与封闭应力固有的统计力学属性,研究了岩体非均匀变形与封闭应力特性,以及深部硐室围岩的长期稳定性等两个力学问题。给出了岩体非均匀构造与封闭应力的数学表征;根据质量守恒定律得到了计算深部硐室围岩长期变形的一般公式;得到了围岩变形中劈裂扩容变形占主要部分的结论,并且阐明了深部围岩卸荷时更易出现劈裂破坏的原因。给出了劈裂破坏形态的演进序列与扩容位移的计算方法。将围岩松动圈范围、破裂区位置和边壁位移的计算结果与锦屏一级电站厂房现有监测数据进行了对比,两者相当吻合。     

Displacement of surrounding rock in a deep circular hole considering double moduli and strength-stiffness degradation

The problem of cavity stability widely exists in deep underground engineering and energy exploitation. First, the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak strength drop model considering the rock's characteristics of constant modulus and double moduli. Then, the orthogonal non-associative flow rule is used to establish the displacement of the surrounding rock under constant modulus and double moduli, respectively, considering the stiffness degradation and dilatancy effects in the plastic region and assuming that the elastic strain in the plastic region satisfies the elastic constitutive relationship. Finally, the evolution of the displacement in the surrounding rock is analyzed under the effects of the double modulus characteristics, the strength drop, the stiffness degradation, and the dilatancy. The results show that the displacement solutions of the surrounding rock under constant modulus and double moduli have a unified expression. The coefficients of the expression are related to the stress field of the original rock, the elastic constant of the surrounding rock, the strength parameters, and the dilatancy angle. The strength drop, the stiffness degradation, and the dilatancy effects all have effects on the displacement. The effects can be characterized by quantitative relationships.     

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

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

一种新的岩石非线性黏弹塑性蠕变模型研究

为了克服传统元件组合模型不能描述岩石蠕变过程中非线性特征的缺陷,首先根据加速蠕变阶段的应变和应变率随蠕变时间急剧增大的特点,建立黏塑性应变与蠕变时间的指数函数关系并提出非线性黏塑性体．将该非线性黏塑性体与广义Burgers蠕变模型串联,建立可以描述岩石全蠕变过程的非线性黏弹塑性蠕变模型,根据叠加原理得到一维应力状态下的轴向蠕变方程．然后基于塑性力学理论指出岩石三维蠕变本构方程建立过程中的不足之处,并给出非线性黏弹塑性蠕变模型合理的三维蠕变方程．最后采用不同应力水平下砂岩轴向蠕变试验对模型合理性进行验证,结果表明：拟合曲线与试验曲线吻合度较高,所建蠕变模型能够很好地描述砂岩在不同应力水平下的蠕变变形规律,尤其对加速蠕变阶段的非线性特征描述效果很好,验证了模型的合理性．     

Effect of the silicon-carbide micro- and nanoparticle size on the thermo-elastic and time-dependent creep response of a rotating Al–SiC composite cylinder

The history of stresses and creep strains of a rotating composite cylinder made of an aluminum matrix reinforced by silicon carbide particles is investigated. The effect of uniformly distributed SiC micro- and nanoparticles on the initial thermo-elastic and time-dependent creep deformation is studied. The material creep behavior is described by Sherby’s constitutive model where the creep parameters are functions of temperature and the particle sizes vary from 50 nm to 45.9 μm. Loading is composed of a temperature field due to outward steady-state heat conduction and an inertia body force due to cylinder rotation. Based on the equilibrium equation and also stress-strain and strain-displacement relations, a constitutive second-order differential equation for displacements with variable and time-dependent coefficients is obtained. By solving this differential equation together with the Prandtl–Reuss relation and the material creep constitutive model, the history of stresses and creep strains is obtained. It is found that the minimum effective stresses are reached in a material reinforced by uniformly distributed SiC particles with the volume fraction of 20% and particle size of 50 nm. It is also found that the effective and tangential stresses increase with time at the inner surface of the composite cylinder; however, their variation at the outer surface is insignificant.     

岩石含水蠕变损伤模型的开发与应用

摘要：为研究岩石含水率对蠕变的影响,将含水开关与蠕变损伤阀值引入蠕变模型,实现了含水劣化与蠕变损伤在本构关系上的耦合,建立了岩石含水蠕变损伤模型。模型与某软岩的三轴压缩蠕变试验结果取得了较好的拟合效果,并由麦夸特法+通用全局优化（LM-UGO）算法反演出了模型参数随岩石含水率的变化规律。在C++语言和FLAC3D中Fish语言的环境下,通过推导模型本构方程的三维差分形式,并利用FLAC3D所提供的接口实现了模型的二次开发。通过对开发的岩石含水蠕变损伤模型进行算例验证和分析,验证了模型的可靠性,也得出了软岩巷道蠕变变形随围岩含水率增大而增加的结论,进一步反映出蠕变模拟计算中考虑岩石含水率的必要性。     

Effect of dilatancy on the strain localization of water-saturated elasto-viscoplastic soil

It is well known that geomaterials such as soils exhibit an increase in volume during shearing deformation, referred to as dilatancy. Dilatancy is a typical property of such granular materials as soils and is closely related to changes in the microstructure. Normally consolidated clay exhibits negative dilatancy or contractancy, namely, a decrease in volume during shearing. On the other hand, overconsolidated clay shows positive dilatancy, namely, an increase in volume during shearing. The aim of the present paper is to study the effects of the microstructure, such as dilatancy and permeability, on the strain localization of water-saturated clay using an elasto-viscoplastic constitutive model. Based on the non-linear kinematic hardening theory and a Chaboche type of viscoplasticity model, an elasto-viscoplastic model for both normally consolidated and overconsolidated clays is proposed; the model can address both negative and positive dilatancies. Firstly, the instability of the model under undrained creep conditions is analyzed in terms of the accelerating creep failure. The analysis shows that clay with positive dilatancy is more unstable than clay with negative dilatancy. Secondly, a finite element analysis of the deformation of water-saturated clay is presented with focus on the numerical results under plane strain conditions. From the present numerical analysis, it is found that both dilatancy and permeability prominently affect shear strain localization behavior.     

聚合物注射成型流动残余应力的数值分析

建立了可压缩黏弹性聚合物熔体在薄壁型腔中充模/保压过程中非等温、非稳态流动 的数学模型,用数值方法实现了注射成型过程中流动应力和取向建立及松弛过程的模拟,研 究了熔体温度、模具温度和注射速率等工艺条件对分子冻结取向的影响,取得了与实验相符 的结果.     

改进的非定常Bingham岩石蠕变模型及参数辨识

传统Bingham模型适用于描述蠕变的前两个阶段,为了更加准确地描述岩体蠕变的三个阶段,本文提出了一种改进的非定常Bingham蠕变模型,并推导出其蠕变本构方程．该模型在传统Bingham模型的基础上串联一个非定常黏壶元件,并考虑岩体弹性模量及原有黏壶元件蠕变参数的非定常化．基于该模型对泥岩的蠕变试验曲线进行拟合及参数辨识,结果表明该模型具有较高的拟合精度,可以准确地描述泥岩蠕变的三个阶段,并实现了在数值软件中的应用,为实际工程提供了一种简单有效的蠕变本构模型．     

海相沉积软土蠕变BP神经网络本构模型

海相沉积软土具有很强的蠕变特性,传统分级构建蠕变本构模型方法不实用,且很难真正反映岩土流变的非线性特性。为此,本文引入具有超强非线性映射和容错能力的BP神经网络模型,通过改进BP算法,根据江门软土的室内直剪蠕变试验结果,建立了海相沉积软土BP神经网络蠕变本构模型,避免了传统方法为满足试验曲线变化规律和蠕变特性而需要建立复杂的本构数学表达式。最后,利用上海地区软土蠕变实验结果对本文提出的方法进行了验证,并对BP神经网络蠕变模型在描述软土流变方面的特点进行了讨论。结果表明,本文建模方法简单,并能很好地描述软土的非线性蠕变问题。     

Experimental stress analysis by photo-rheologic method

In the photo-rheologic method used in this paper, the time-dependent creep strain rate and the instantaneous elastic-plastic strain rate are considered. Unlike ordinary photoplasticity, the strain history and the time dependence are accounted for in this stress analysis by characteristic values independent of strain rate. This method is also suitable for analyzing finite deformations. Using this method, the time-dependent deformations and stresses in a body can be analyzed with the same accuracy as photoelasticity from the instant of loading to the static-equilibrium state during plastic flow or creep. Fundamental concepts and procedures are described. Two examples are used to demonstrate possible applications of this method.     

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

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