Microcrack interaction brittle rock subjected to uniaxial tensile loads

A micromechanics-based model is proposed to describe unstable damage evolution in microcrack-weakened brittle rock material. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress–strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to be consistent with the experimental results.     

Constitutive relationship of brittle rock subjected to dynamic uniaxial tensile loads with microcrack interaction effects

A micromechanical model is proposed to describe both stable and unstable damage evolution in microcrack-weakened brittle rock material subjected to dynamic uniaxial tensile loads. The basic idea of the present model is to classify the constitution relationship of rock material subjected to dynamic uniaxial tensile loads into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress-strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to consistent with the experimental results.     

Zonal fracturing mechanism in deep crack-weakened rock masses

The mechanical behaviors of deep crack-weakened rock masses are different from those of shallow crack-weakened rock masses. The surrounding rock in shallow crack-weakened rock mass engineering is classified into loose zone, plastic zone and elastic zone, while the surrounding rock in deep crack-weakened rock mass engineering is classified into fractured zone and non-fractured zone, which occur alternatively. It is assumed that the deep rock masses contain one joint set, in which the probability density function describing the distribution of sizes is assumed to follow the Rayleigh distribution, and the probability density function describing the distribution of spacing is assumed to follow the Weibull distribution. On the basis of strength criterion of deep rock mass, the near-field stress redistribution around circular opening induced by excavation is determined. The strong interaction among cracks is investigated by using the dislocation model. The nucleation, growth, interaction and coalescence of cracks were analyzed based on the strain energy density factor theory. When cracks coalesce, failure of deep crack-weakened rock masses occurs, fractured zone is formed. Then, size and quantity of fractured zone and non-fractured zone are given out. The size and quantity of fractured zone increase with decreasing strength of rock mass. The size and quantity of fractured zone increase with increasing in situ stress. Zonal fracturing phenomenon occurs once value of in situ stress is larger than the unaxial compressive strength of rock masses. The size and quantity of fractured zone decrease with increasing λ when p₂ > p₁. The size and quantity of fractured zone increase with increasing λ when p₂ < p₁.     

岩石单轴压缩作用下变形局部化的梯度塑性解

采用梯度塑性理论研究单轴压缩作用下岩石变形局部化,得到了单轴压缩作用下岩石变形局部化带宽度的一维、二维解析解,为实验测定内部材料长度参数提供了理论依据.     

寒区冻融环境条件下岩石损伤扩展研究探讨

本文借助于岩石材料冻融循环CT扫描实验，对岩石材料在冻融循环条件下损伤扩展特性进行研究，着重探讨了冻融循环次数对岩石材料损伤扩展的影响，并对冻融循环条件下岩石损伤扩展本构关系进行了探讨。     

考虑裂纹闭合效应的岩石损伤本构关系

岩石中的预存裂纹只有在一定的法向压应力即裂纹闭合应力的作用下才可能闭合，其闭合过程与其方位和外加应力场有关，并且，即使对于裂纹已经完全闭合的岩石，如果裂纹闭合应力不同，则岩石的应力应交关系也不相同。本文建立了考虑裂纹闭合效应的岩石细观损伤力学模型，分析了裂纹闭合应力对岩石损伤演化过程和应力应变关系的影响。数值结果表明裂纹闭合应力显著地改变岩石的应力应变关系，表现为随裂纹闭合应力的增加，岩石的轴向应变变化较小，侧向应变和体积应变则大为增加。     

Triaxial compressive behavior of rock with mesoscopic heterogenous behavior: Strain energy density factor approach

The strain energy density factor approach is used in conjunction with a micromechanics model to investigate the condition and direction of shear failure for brittle rock subjected to triaxial compression. Moderate confinement in addition to localized deformation and damage are considered. Quantified are the effects of the various geometric and load parameters that involve the interaction of microcrack, friction and the confining pressure such that the path of the wing crack is taken into account. The influence of all microcracks with different orientations are introduced into the constitutive relation. The closed-form solution for the complete stress–strain relation of rock containing microcracks is obtained. It is shown that the complete stress–strain relationship includes linear, nonlinear hardening, rapid stress drop and strain softening effects. The theoretical results show that deviation of the direction of wing cracks from the line of the pre-existing crack decreases with increasing confinement pressure and friction coefficient. Theoretical predictions and experimental results show good agreement.     

不同应变率下煤岩破坏特征及其本构模型

利用直径50 mm的分离式霍普金森压杆,对煤岩展开20～100 s?1动态应变率下的单轴冲击压缩试验,结合高速摄影分析其变形破坏特征,并建立基于Weibull统计分布和Drucker-Prager破坏准则的煤岩动态强度型统计损伤本构模型。试验结果表明:（1）煤岩动态应力-应变曲线存在明显的非线性特征,随应变率升高,动态抗压强度与弹性模量均呈线性增长且增幅显著,破坏形态由低应变率下的轴向劈裂破坏向高应变率下的压碎破坏过渡;（2）在动态应变率20～100 s?1下,煤岩破坏后碎块具有明显的分形特性,破碎块度分维值为1.9～2.2,且随着应变率的升高,煤岩破碎程度增大,碎块块度减小;（3）基于Weibull分布参数F₀、m和应变率的关系,修正煤岩的本构模型,并与试验结果进行对比,验证该模型的合理性。     

Computer simulation of effect of grain size distribution on Weibull parameters

Computer simulation of crack nucleation associated with unstable equilibrium in polycrystal consisting of grains of random sizes and orientations are presented. Micro-stress and micro-strain fluctuations were taken into account. Distribution functions for probability of fracture for different parameters of grain structure were obtained. It is shown that these functions may be approximated by three-parameter Weibull distribution with high accuracy. Ascertained are the dependencies of Weibull parameters on both grain sizes and their variance.     

基于0.2定围压比状态下的约束混凝土性能研究

通过对混凝土施加0.2定围压比,使单轴受压混凝土裂缝发展受到约束,研究了混凝土在接近无裂纹发展状态下的受压力学性能。本文以在MTS815.02岩石电液伺服系统上对不同强度等级的混凝土试件在定围压比为0.2时的单调、循环加卸载受压常规三轴试验为基础,研究分析了0.2定围压比对混凝土受压应力-应变规律的影响和围压对混凝土裂缝发展的限制效果,并根据0.2定围压比下混凝土材料的特殊性能,提出了无裂纹发展状态下混凝土上升段表达式及参数。分析结果表明,0.2定围压比下,各强度混凝土本构关系均表现出了类似金属材料的本构关系,且低强混凝土试件应力-应变的提高幅度优于中强混凝土试件。综上所述,0.2定围压比改善了混凝土在受压时无征兆失效(脆性破坏)的力学性能。     

隧洞围岩力学性态预报的随机介质模型

本文基于系统学观点,利用随机过程以及概率与数理统计方法,视岩体为一随机不确定性系统,建立了描述变形概率密度函数演化所服从的方程,详细讨论了变形发展的概率预报,最后给出了一些有意义的结论和意见。     

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

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

岩石化学腐蚀-围压细观损伤本构模型与离散元模拟研究

基于Lemaitre应变等价性假设理论,假定受水化学-力耦合损伤的岩石微元强度服从Weibull分布,考虑化学腐蚀与围压耦合作用对岩石力学参数的影响,通过核磁共振技术与损伤力学理论,引入细观化学损伤变量与力损伤变量,并认为微元破坏符合SMP准则,建立岩石化学腐蚀-力耦合损伤本构模型,并采用理论推导的方法得出所需的模型参数。同时基于颗粒离散元方法,引入参数半径乘子来改变颗粒间的黏结接触尺寸,从而模拟水化学损伤,采用平直节理模型对水化学作用后的岩石进行三轴压缩模拟,得到了水化学作用和不同围压下的岩石三轴应力-应变模拟曲线。通过对比所构建的岩石化学腐蚀-力耦合损伤本构模型理论曲线、离散元模拟曲线和试验曲线,结果表明三者吻合度较好,能够很好地反映岩石在化学腐蚀和围压耦合作用下的力学特性与破坏特征,并通过离散元方法得到了岩石在三轴压缩过程中裂纹的产生与分布情况。     

非Gauss分布随机波浪力对海洋平台的非线性作用

为避免求解决定Maikov过程转移概率密度的Fokker—Planck方程，基于尺度分离的假设导出了一组描述非线性海洋平台受非Gauss分布随机波浪载荷作用所产生响应的矩量的常微分方程组。矩量方程清楚地反映出分别对应随机载荷和结构响应的两种不同统计特性的相互关系。由于矩量方程不依赖载荷的概率分布的具体细节，以它来模拟随机激励作用下的非线性系统将免于Monte Carlo方法所面临的正确模拟载荷概率分布的困难任务。将摄动法用于矩量方程可使线性化不再需要，这样就不会因为线性化而产生不可预料的误差。     

DOUBLE-MEDIUM CONSTITUTIVE MODEL OF GEOLOGICAL MATERIAL IN UNIAXIAL TENSION AND COMPRESSION

Based on elasto-plasticity and damage mechanics, a double-medium constitutive model of geological material under uniaxial tension and compression was presented, on the assumption that rock and soil materials are the pore-fracture double-medium, and porous medium has no damage occurring, while fracture medium has damage occurring with load. To the implicit equation of the model, iterative method was adopted to obtain the complete stress-strain curve of the material. The result shows that many different distributions (uniform distribution, concentrated distribution and random distribution) of fractures in rock and soil material are the essential reasons of the daedal constitutive relations. By the reason that the double-medium constitutive model separates the material to be porous medium part, which is the main body of elasticity, and fracture medium part, which is the main body of damage, it is of important practical values and theoretical meanings to the study on failure of rock and soil or materials containing damage.     

考虑微裂纹相互作用的的岩石微观力学弹塑性损伤模型研究

本文建立基于微裂纹扩展的岩石弹塑性损伤微观力学模型。用自洽方法考虑裂隙间相互影响,压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的微观损伤,基于应变能密度准则用Newton迭代法求复合型断裂的翼裂纹扩展长度,并采用微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映塑性内变量对硬化函数的影响;建立岩石的弹塑性损伤本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性损伤模型的程序。从围压和微裂隙长度等因素分析弹塑性损伤模型的岩石的损伤和宏观塑性特性。     

Non-linear constitutive relations for magnetostrictive materials

In this paper, non-linear deformation behavior of magnetostrictive materials is studied and three magnetoelastic coupling constitutive models are developed. The standard square (SS) constitutive model is developed by means of truncating the polynomial expansion of the Gibbs free energy. The hyperbolic tangent (HT) constitutive equations, which involve a hyperbolic tangent magnetic-field dependence, are proposed to model the magnetic-field-induced strain saturation of magnetostrictive materials in the region of intense magnetic fields. A new model based on density of domain switching (DDS) is established in terms of the basic truth that magnetic domain switching underlies magnetostrictive deformation. In this model, it is assumed that the relation between density of domain switching, defined by the quantity of magnetic domains switched by per unit magnetic field and magnetic field can be described by a density function with normal distribution. The moduli in these constitutive models can be determined by a material function that is proposed to describe the dependence of the peak piezomagnetic coefficient on the compressive pre-stress for one-dimensional cases based on the experimental results published. The accuracy of the non-linear constitutive relations is evaluated by comparing the theoretical values with experimental results of a Terfenol-D rod operated under both compressive pre-stress and bias magnetic field. Results indicate that the SS constitutive equations can accurately predict the experimental results under a low or moderate magnetic field while the HT model can, to some extent, reflect the trend of saturation of magnetostrictive strain under a high magnetic field. The model based on DDS, which is more effective in simulating the experimental curves, can capture the main characteristics of the mechanism of magnetoelastic coupling deformation of a Terfenol-D rod, such as the notable dependence of magnetoelastic response on external stress and the saturation of magnetostrictive strain under intense magnetic fields. In addition, the SS constitutive relation for a general three-dimensional problem is discussed and an approach to characterize the modulus tensors is proposed.     

Creep rupture of specimens with random material properties

Possible sources of the observed scatter in creep deformation and rupture time are discussed. The scatter due to random material properties in creep is considered.A constitutive equation describing these random material properties is formulated. The rupture time of a specimen under constant load is calculated using Hoff s theory of ductile rupture. Rupture is shown to occur at finite elongation, in contrast to Hoff's original analysis.The probability distribution of the rupture time is determined with extremum value analysis and compared with experimentally observed ones.     

Upper and lower bounds for constitutive relation of crack-weakened rock masses under dynamic compressive loads

A micro-mechanics-based model is proposed to investigate the rate-dependent constitutive relation for crack-weakened rock masses subjected to dynamic compressive loads. The present micro-mechanical model reveals that the nucleation, growth and coalescence of sliding cracks dominate the failure and macroscopic properties of crack-weakened rock masses subjected to dynamic compressive loads. The interactions among multiple parallel sliding cracks in crack-weakened rock masses subjected to dynamic compressive loads are examined asymptotically in an explicit and quantitative manner in order to reveal fully their so-called shielding and magnification effects on the stress–strain relation. Based on the micro-mechanical framework and the asymptotic analysis, analytical upper and lower bounds are proposed for the rate-relation for rock masses containing multiple rows of echelon cracks subjected to dynamic compressive loads. The factors that affect the rate-dependent properties of crack-weakened rock masses have been analyzed. The strain energy density factor approach, which is related to crack growth velocity and dynamic fracture toughness of rock material, is employed in the analysis. The rate-dependent constitutive relation of crack-weakened rock masses is derived from micro-mechanical framework and the asymptotic analysis. The closed-form explicit expression for the rate-dependent constitutive relation of rock masses containing echelon cracks subjected to dynamic compressive loads is obtained. Finally, the present model is used to analyze the complete stress–strain relation and strength for jointed rock masses at shiplock slope of the Three Gorges Dam.     

Micromechanical definition of an entropy for quasi-static deformation of granular materials

A micromechanical theory is formulated for quasi-static deformation of granular materials, which is based on information theory. A reasoning is presented that leads to the definition of an information entropy that is appropriate for quasi-static deformation of granular materials. This definition is based on the hypothesis that relative displacements at contacts with similar orientations are independent realisations of a random variable. This hypothesis is made plausible based on the results of Discrete Element simulations. The developed theory is then used to predict the elastic behaviour of granular materials in terms of micromechanical quantities. The case considered is that of two-dimensional assemblies consisting of non-rotating particles with an elastic contact constitutive relation. Applications of this case are the initial elastic (small-strain) deformation of granular materials. Theoretical results for the elastic moduli, relative displacements, energy distribution and probability density functions are compared with results obtained from the Discrete Element simulations for isotropic assemblies with various average numbers of contacts per particle and various ratios of tangential to normal contact stiffness. This comparison shows that the developed information theory is valid for loose systems, while a theory based on the uniform-strain assumption is appropriate for dense systems.