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

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

脆性岩石蠕变裂纹成核宏细观力学机理研究

脆性岩石内部细观裂纹扩展、贯通及成核影响下的脆性蠕变行为, 对深部地下工程围岩微震及岩爆事件评价有着重要意义. 然而, 目前能够解释裂纹成核损伤突变影响下, 脆性岩石完整蠕变宏细观力学机理模型的研究很少. 本文基于脆性岩石亚临界裂纹扩展模型、裂纹-应变-声发射事件相关的损伤模型及裂纹成核损伤时间演化路径函数, 提出了一种脆性岩石裂纹成核损伤突变影响下的蠕变宏细观力学模型. 裂纹成核损伤时间演化路径函数通过岩石内部裂纹成核损伤突变大小$\Delta D_{CN}$及相邻裂纹成核损伤时间差$\Delta t$进行定义, 该函数可以结合岩石声发射监测试验数据定义的岩石损伤数据确定. 通过与试验结果对比分析验证模型的合理性. 并讨论了裂纹成核损伤大小、相邻裂纹成核损伤时间间隔、及裂纹成核数量对脆性岩蠕变裂纹长度、裂纹速率、轴向应变及应变率的影响. 该结果对于更加合理、经济、高效的深部地下工程施工及设计提供了一定的理论支持.      

A micro–macro method for predicting the shear strength of brittle rock under compressive loading

Microcracks have great significance for shear strength of brittle rock in compression. A major challenge of this area is to establish the correlation of microcracks and macroscopic shear strength. A new micro–macro method is presented to predict the shear strength of brittle rock in compression. This method incorporates the microcrack model suggested by Ashby, Mohr–Coulomb failure criterion and a crack-strain relation. This crack–strain relation is presented to link the crack growth and axial strain by combining the micro and macro definitions from rock damage. The shear strength and stress–strain relationship of Jinping marble are theoretically investigated in detail. The rationality of this suggested method is verified by using the experimental results founded on Jinping marble. Effects of the initial microcrack size, friction coefficient and confining pressure on internal friction angle, cohesion, and shear strength are also discussed.     

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.     

Damage accumulation for growth of anti-plane crack in work-hardening material

Elastic–plastic solutions of an anti-plane crack in an infinite body are used in conjunction with a continuum damage model to describe the conditions necessary for the onset of crack instability, fatigue crack propagation due to cyclic loading, and rates of crack growth due to time dependent events. A power law relates the stress to the strain of the material. The damage, which invokes nucleation, growth and coalescence of microvoids due to elevated strain, is confined to the plastic zone surrounding the crack tip. For applied loading below the yield stress, the small-scale and large-scale yielding solutions are used to determine the influence of strain hardening on crack instability and failure. Crack growth due to cyclic loading and time-dependent deformations are studied using the small-scale yielding solution of the deformation theory of plasticity.     

Localization of deformation and stress–strain relation for mesoscopic heterogeneous brittle rock materials under unloading

Stress redistribution induced by excavation of underground engineering and slope engineering results in the unloading zone in parts of surrounding rock masses. The mechanical behaviors of crack-weakened rock masses under unloading are different from those of crack-weakened rock masses under loading. A micromechanics-based model has been proposed for brittle rock material undergoing irreversible changes of their microscopic structures due to microcrack growth when axial stress is held constant while lateral confinement is reduced. The basic idea of the present model is to classify the constitution relation of rock material 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 closed-form explicit expression for the complete stress–strain relation of rock materials containing cracks under unloading is obtained. The results show that the complete stress–strain relation and the strength of rock materials under unloading depend on the crack spacing, the fracture toughness of rock materials, orientation of the cracks, the crack half-length and the crack density parameter.     

Discrete element modeling on the crack evolution behavior of brittle sandstone containing three fissures under uniaxial compression

Based on experimental results of brittle, intact sandstone under uniaxial compression, the micro-parameters were firstly confirmed by adopting particle flow code(PFC2D). Then, the validation of the simulated models were cross checked with the experimental results of brittle sandstone containing three parallel fissures under uniaxial compression. The simulated results agreed very well with the experimental results, including the peak strength, peak axial strain, and ultimate failure mode. Using the same microparameters, the numerical models containing a new geometry of three fissures are constructed to investigate the fissure angle on the fracture mechanical behavior of brittle sandstone under uniaxial compression. The strength and deformation parameters of brittle sandstone containing new three fissures are dependent to the fissure angle. With the increase of the fissure angle, the elastic modulus, the crack damage threshold,and the peak strength of brittle sandstone containing three fissures firstly increase and secondly decrease. But the peak axial strain is nonlinearly related to the fissure angle. In the entire process of deformation, the crack initiation and propagation behavior of brittle sandstone containing three fissures under uniaxial compression are investigated with respect tothe fissure angle. Six different crack coalescence modes are identified for brittle sandstone containing three fissures under uniaxial compression. The influence of the fissure angle on the length of crack propagation and crack coalescence stress is evaluated. These investigated conclusions are very important for ensuring the stability and safety of rock engineering with intermittent structures.     

岩爆岩石断裂的微观结构形貌分析及岩爆机理

利用扫描电镜(SEM)对岩爆岩石断口微观形貌特征进行研究分析,从微观角度探索岩爆产生的机理。通过对平顶山十二矿岩爆现场取样对其断口形貌特征与地应力和岩石成分之间关系进行研究。巷道围岩劈裂岩块断口形貌多呈台阶状,劈裂面与地应力最大主应力方向平行,岩石断口属拉张断裂,劈裂纹的产生主要是脆性断裂;岩爆抛射出的岩块断口形貌非常复杂,裂面与切应力(最大主应力)方向平行或相交,不同平面内的微裂纹通过与岩爆裂纹间的微裂纹或受撕裂作用形成台阶,表面不平整,属于拉张或剪切型断裂。岩石细观成分对岩爆的影响也较大,结晶程度高、结构致密的硬脆岩石更易发生岩爆。     

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.     

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.     

An experimental study on fracture mechanical behavior of rock-like materials containing two unparallel fissures under uniaxial compression

Strength and deformability characteristics of rock with pre-existing fissures are governed by cracking behavior. To further research the effects of pre-existing fissures on the mechanical properties and crack coalescence process, a series of uniaxial compression tests were carried out for rock-like material with two unparallel fissures.In the present study, cement, quartz sand, and water were used to fabricate a kind of brittle rock-like material cylindrical model specimen. The mechanical properties of rock-like material specimen used in this research were all in good agreement with the brittle rock materials. Two unparallel fissures(a horizontal fissure and an inclined fissure) were created by inserting steel during molding the model specimen.Then all the pre-fissured rock-like specimens were tested under uniaxial compression by a rock mechanics servocontrolled testing system. The peak strength and Young's modulus of pre-fissured specimen all first decreased and then increased when the fissure angle increased from 0?to 75?.In order to investigate the crack initiation, propagation and coalescence process, photographic monitoring was adopted to capture images during the entire deformation process.Moreover, acoustic emission(AE) monitoring technique was also used to obtain the AE evolution characteristic of prefissured specimen. The relationship between axial stress, AE events, and the crack coalescence process was set up: when a new crack was initiated or a crack coalescence occurred, thecorresponding axial stress dropped in the axial stress–time curve and a big AE event could be observed simultaneously.Finally, the mechanism of crack propagation under microscopic observation was discussed. These experimental results are expected to increase the understanding of the strength failure behavior and the cracking mechanism of rock containing unparallel fissures.     

A physically motivated anisotropic tensorial representation of damage with separate functions for void nucleation,growth, and coalescence

A phenomenological anisotropic damage progression formulation for porous ductile metals with second phases is described through mechanisms of void nucleation, growth and coalescence. The model is motivated from fracture mechanisms and microscale physical observations. To describe the creation of new pores, the decohesion at the particle–matrix interface and the fragmentation of second phase particles, the void-crack nucleation equation is related to several microstructural parameters (fracture toughness, length scale parameter, particle size, volume and fraction of second phase), the plastic strain level, and the stress state. Nucleation is represented by a general symmetric second rank tensor, and its components are proportional to the absolute value of the plastic strain rate components. Based on the Rice and Tracey model, void growth is a scalar function of the trace of damage tensor and the positive triaxiality. Like nucleation, coalescence is a second rank tensor governed by the plastic strain rate tensor and the stress state. The coalescence threshold is related to the void length scale for void impingement and void sheet mechanisms. The coupling of damage with the Bammann–Chiesa–Johnson (BCJ) plasticity model is written in the thermodynamic framework and derives from the concept of effective stress assuming the hypothesis of energy equivalence. A full-implicit algorithm is used for the stress integration and the determination of the consistent tangent operator. Finally, macroscale correlations to cast A356 AL alloy and wrought 6061-T6 AL alloy experimental data are completed with predictive void-crack evolution to illustrate the applicability of the anisotropic damage model.     

A fracture damage model for Jointed rock masses and its application to the stability analysis of dam abutments

Based on continuum damage mechanics, for jointed rock masses, a fracture damage model is presented in this paper. First, the damage tensors are defined through the elastic-flexibility of intact rock and the equivalent elastic-damage flexibility for rock mass. Then, by the self-consistent principle of solid mechanics, the equivalent elastic-damage flexibility tensors involving the interaction between multicracks are deduced. The damage evolution law is proposed involving the mechanism of crack propagation process: frictional sliding, crack kinking, growing of branched tension cracks, interlinking of the microcracks near branched crack tips leading to the breakthrough of macro-cracks and finally the failure of rock mass. Thus the evolution of damage variables reasonably unified with the process of crack propagation is given. Finally, a plastic-brittle damage constitutive relation including brittle coupled strain rate, developed and applied to the stability analysis of complicated rock foundation of a dam in China, is described in this paper.     

Constitutive and damage model for 63Sn37Pb solder under uniaxial and torsional cyclic loading

A constitutive model with Ohno–Wang kinematic hardening rule is developed and employed to simulate the isothermal cyclic behavior of Sn–Pb solder under uniaxial and torsional loading. An implicit constitutive integration scheme is presented for inelastic flow of solder. Then a modified low cycle fatigue life prediction model is put forward in which the sum of maximum shear strain range and normal strain range based on the critical plane concept is adopted to replace the uniaxial strain range used by Stolkarts et al. [Stolkarts, V., Keer, L.M., Fine, M.E., 1999. Damage evolution governed by microcrack nucleation with application to the fatigue of 63Sn–37Pb solder. J. Mech. Phys. Solids 47, 2451–2468]. Comparison of the experimental results and simulation verifies that the stress strain hysteresis loops and peak stress decline curve of solder can be reasonably modeled over a wide range of loading conditions with implement of damage coupled constitutive model, and the lifetime estimations of 63Sn37Pb solder based on the assumption of microcrack nucleation governed damage is effective to provide a conservative prediction.     

Fractal analysis of fracture in concrete

Experimental results indicate that propagation paths of cracks in concrete are often irregular, producing rough fracture surfaces which are fractal. Based on dynamic analysis of microcrack coalescence, this paper presents a statistical fractal model to describe the damage evolution of concrete. The model demonstrates that the mechanism of fracture surfaces formed in concrete is closely related to the dynamic processes of the cascade coalescence of microcracks. A unimodal relation between the fractal dimension and the coalescence threshold can qualitatively explain the relation between fractal dimension and fracture energy.     

DAMAGE MODEL FOR HARDENING AND SOFTENING MATERIAL IN GENERAL STATE OF STRESS

The present paper is aimed to simulate progression of damage,hardeningand softening response in brittle materials such as concrete or rock in general state ofstress.Similar shape of surfaces for yield,failure and damage progressing areavailable,and softening strain is treated as plasticity.Then,the proposed model isapplied to solving several boundary value problems.     

Damage and self-similarity in fracture

Consider applications of damage mechanics to material failure. The damage variable introduced in damage mechanics quantifies the deviation of a brittle solid from linear elasticity. An analogy between the metastable behavior of a stressed brittle solid and the metastable behavior of a superheated liquid is established. The nucleation of microcracks is analogous to the nucleation of bubbles in the superheated liquid. In this paper we have applied damage mechanics to four problems. The first is the instantaneous application of a constant stress to a brittle solid. The results are verified by applying them to studies of the rupture of chipboard and fiberglass panels. We then obtain a solution for the evolution of damage after the instantaneous application of a constant strain. It is shown that the subsequent stress relaxation can reproduce the modified Omori’s law for the temporal decay of aftershocks following an earthquake. Obtained also are the solutions for application of constant rates of stress and strain. A fundamental question is the cause of the time delay associated with damage and microcracks. It is argued that the microcracks themselves cause random fluctuations similar to the thermal fluctuations associated with phase changes.     

A micromechanics-based strain gradient damage model for fracture prediction of brittle materials – Part II: Damage modeling and numerical simulations

In this paper, we established a strain-gradient damage model based on microcrack analysis for brittle materials. In order to construct a damage-evolution law including the strain-gradient effect, we proposed a resistance curve for microcrack growth before damage localization. By introducing this resistance curve into the strain-gradient constitutive law established in the first part of this work (Li, 2011), we obtained an energy potential that is capable to describe the evolution of damage during the loading. This damage model was furthermore implemented into a finite element code. By using this numerical tool, we carried out detailed numerical simulations on different specimens in order to assess the fracture process in brittle materials. The numerical results were compared with previous experimental results. From these studies, we can conclude that the strain gradient plays an important role in predicting fractures due to singular or non-singular stress concentrations and in assessing the size effect observed in experimental studies. Moreover, the self-regularization characteristic of the present damage model makes the numerical simulations insensitive to finite-element meshing. We believe that it can be utilized in fracture predictions for brittle or quasi-brittle materials in engineering applications.     

Experimental study and numerical modeling of brittle fracture of carbonate rock under uniaxial compression

The brittle carbonate rock taken from the Tarim Oilfield is tested in laboratory under uniaxial compression. The acoustic emission (AE) is used to monitor the microcracking activity in rock during the experiment. Moreover, the 3D tomograms of carbonate rock after uniaxial compression are obtained by using CT imaging technology, which indicates that microcracks mutually interconnect and eventually form macroscopic fractures after failure. The PFC2D is used to model the behavior of brittle rock including microcracks propagation. The stress–strain curve and cracks distribution in rock model are obtained from the PFC simulation. The numerical results agree with the experimental test well.     

岩石紧凑拉伸断裂过程及其尺寸效应的三维数值模拟

采用最近开发的三维岩石破裂过程分析软件RFPA3D模拟单边裂纹紧凑拉伸断裂过程。试验中五个不同尺寸的岩样具有相同的力学性质参数分布,模拟结果得到了裂纹扩展中的应力场、位移场和声发射的空间分布以及单边裂纹扩展贯通的过程。单边裂纹拉伸断裂的路径是一个复杂的空间三维曲面,三维裂纹比二维裂纹更为复杂。分析了岩石试样的峰值强度和试样尺寸之间的关系。随着岩样尺寸的增加,峰值强度逐渐减小,并且延性破坏特征更加明显,模拟结果满足岩石的尺寸效应规律。最后模拟了三组不同均匀性的试样拉伸破坏过程,结果表明细观上的非均匀性对岩石尺寸效应有很大影响,随着非均匀性的增加,岩石宏观强度随之提高,即使在均匀材料中一样存在尺寸效应。