3D Anisotropic Elastoplastic-Damage Model and its Application in Simulating the Behavior of Rock Materials

A 3D anisotropic elastoplastic-damage model was presented based on continuum damage mechanics theory. In this model, the tensor decomposition technique is employed. Combined with the plastic yield rule and damage evolution, the stress tensor in incremental format is obtained. The derivate eigenmodes in the proposed model are assumed to be related with the uniaxial behavior of the rock material. Each eigenmode has a corresponding damage variable due to the fact that damage is a function of the magnitude of the eigenstrain. Within an eigenmodes, different damage evolution can be used for tensile and compressive loadings. This model was also developed into finite element code in explicit format, and the code was integrated into the well-known computational environment ABAQUS using the ABAQUS/Explicit Solver. Numerical simulation of an uniaxial compressive test for a rock sample is used to examine the performance of the proposed model, and the progressive failure process of the rock sample is unveiled.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Damage evaluation and damage localization of rock

Knowledge of damage accumulation and corresponding failure evolution are prerequisite for effective maintenance of civil engineering so as to avoid disaster. Based on statistical mesoscopic damage mechanics, it was revealed that there are three stages in the process of deformation, damage and failure of multiscale heterogeneous elastic-brittle medium. These are uniformly distributed damage, localized damage and catastrophic failure. In order to identify the transitions from scattering damage to macroscopically localized one, a condition for damage localization was given. The experiments of rock under uniaxial compression with the aid of observations of acoustic emission and speckle correlation do support the concept of localization. This provides a potential approach to properly evaluate damage accumulation in practice. In addition, it is found in the experiments that catastrophic failure displays critical sensitivity. This gives a helpful clue to the prediction of catastrophic failure. (C) 2004 Elsevier Ltd. All rights reserved.

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Evolution of Localized Damage Zone in Heterogeneous Media

Evolution of localized damage zone is a key to catastrophic rupture in heterogeneous materials. In the present article, the evolutions of strain fields of rock specimens are investigated experimentally. The observed evolution of fluctuations and autocorrelations of strain fields under uniaxial compression demonstrates that the localization of deformation always appears ahead of catastrophic rupture. In particular, the localization evolves pronouncedly with increasing deformation in the rock experiments. By means of the definition of the zone with high strain rate and likely damage localization, it is found that the size of the localized zone decreases from the sample size at peak load to an eventual value. Actually, the deformation field beyond peak load is bound to suffer bifurcation, namely an elastic unloading part and a continuing but localized damage part will co-exist in series in a specimen. To describe this continuous bifurcation and localization process observed in experiments, a model on continuum mechanics is developed. The model can explain why the decreasing width of localized zone can lead stable deformation to unstable, but it still has not provided the complete equations governing the evolution of the localized zone.

Relationship between load/unload response ratio and damage variable and its application

The physics-based parameter: load/unload response ratio (LURR) was proposed to measure the proximity of a strong earthquake, which achieved good results in earthquake prediction. As LURR can be used to describe the damage degree of the focal media qualitatively, there must be a relationship between LURR and damage variable (D) which describes damaged materials quantitatively in damage mechanics. Hence, based on damage mechanics and LURR theory, taking Weibull distribution as the probability distribution function, the relationship between LURR and D is set up and analyzed. This relationship directs LURR applied in damage analysis of materials quantitatively from being qualitative earlier, which not only provides the LURR method with a more solid basis in physics, but may also give a new approach to the damage evaluation of big scale structures and prediction of engineering catastrophic failure. Copyright (c) 2009 John Wiley & Sons, Ltd.

Relationship between load/unload response ratio and damage variable and its application

The physics-based parameter: load/unload response ratio (LURR) was proposed to measure the proximity of a strong earthquake, which achieved good results in earthquake prediction. As LURR can be used to describe the damage degree of the focal media qualitatively, there must be a relationship between LURR and damage variable (D) which describes damaged materials quantitatively in damage mechanics. Hence, based on damage mechanics and LURR theory, taking Weibull distribution as the probability distribution function, the relationship between LURR and D is set up and analyzed. This relationship directs LURR applied in damage analysis of materials quantitatively from being qualitative earlier, which not only provides the LURR method with a more solid basis in physics, but may also give a new approach to the damage evaluation of big scale structures and prediction of engineering catastrophic failure. Copyright (c) 2009 John Wiley & Sons, Ltd.

Multiscale Coupling in Complex Mechanical Systems

Multiscale coupling attracts broad interests from mechanics, physics and chemistry to biology. The diversity and coupling of physics at different scales are two essential features of multiscale problems in far-from-equilibrium systems. The two features present fundamental difficulties and are great challenges to multiscale modeling and simulation. The theory of dynamical system and statistical mechanics provide fundamental tools for the multiscale coupling problems. The paper presents some closed multiscale formulations, e.g., the mapping closure approximation, multiscale large-eddy simulation and statistical mesoscopic damage mechanics, for two typical multiscale coupling problems in mechanics, that is, turbulence in fluids and failure in solids. It is pointed that developing a tractable, closed nonequilibrium statistical theory may be an effective approach to deal with the multiscale coupling problems. Some common characteristics of the statistical theory are discussed.

Validation of component assembly model and extension to plasticity

Material potential energy is well approximated by '' pair-functional '' potentials. During calculating potential energy, the orientational and volumetric components have been derived from pair potentials and embedding energy, respectively. Slip results in plastic deformation, and slip component has been proposed accordingly. Material is treated as a component assembly, and its elastic, plastic and damage properties are reflected by different components respectively. Material constitutive relations are formed by means of assembling these three kinds of components. Anisotropy has been incorporated intrinsically via the concept of component. Theoretical and numerical results indicate that this method has the capacity of reproducing some results satisfactorily, with the advantages of physical explicitness, etc. (c) 2007 Elsevier Ltd. All rights reserved.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

三峡库区巴东组泥灰岩溶蚀作用与岩体力学性质弱化研究

Abstract The karsrt erosion engineering geology became a highlight problem in recent years, in particularly, the karst erosion of marlite of Badong formation made the rock mechanics weaken in Three Gorges Reservoir area, which reduces the safety of slope. During the immigrant construction, many high slopes have been formed, whose instabilities problems pose serious threats to the safety of the people and properties. The accidents of the slope failure take place now and then. By testing, it has been found that the karst erosion pattern and dissolution rate of marlite are not weaker than that of the pure limestone. Furthermore, owning to the weathering and unloading, the karst erosion of the marlite will reach certain depth of the slope, which is named infiltrated karst erosion. The karst erosion made the rock mass quality of slope or foundation worse in a large scale. The karst erosion geological disasters, taken place or not, has become the main restrictive factors to the social stability and economic development. Thus the karst erosion process and mechanism of marlite of Badong formation are studied as the main content of this dissertation. The weakening characteristic of rock mass mechanics parameters are studied along with the rock mass structure deformation and failure processes in the course of the karst erosion. At first, the conditions and influencing factors of the karst erosion are analyzed in the investigative region, on the basis of different karst erosion phenomenon of the marlite and different failure modes of slope. Then via indoor the karst erosion tests, it is analyzed that the karst erosion will change the rock mass composition and its structure. Through test, the different karst erosion phenomena between micro and macro have been observed, and the karst erosion mechanism of the marlite has been summarized. Damage theory is introduced to explain the feature of dissolution pore and the law of crack propagation in the marlite. By microscope and the references data, it can be concluded that the karst erosion process can be divided into rock minerals damage and rock structural damage. And the percent of karst erosion volume is named damage factor, which can be used to describe the quantify karst erosion degree of marlite. Through test, the rock mechanical properties in the different period of karst erosion are studied. Based on the damage mechanics theory and the test result, the relation between the karst erosion degree of marlite and weakening degree of mechanical properties is summarized. By numerical simulations, the karst erosive rock mass mechanics is verified. The conclusion is drawn as below: to the rock mass of marlite, the karst erosion damage made mechanics parameters variation, the deformation modulus, cohesion, and inter friction angle reduce as the negative exponent with the increasing of the karst erosion volume, however, the Poisson ratio increases as the positive exponent with the karst erosion volume increasing. It should be noticed that the deduced formulations are limited to the test data and certain conditions. It is suitable to the rock mass parametric weakening process after the karst erosion of marlite in Three Gorges Reservoir area. Based on the failure types of marlite slope in the field, the karst erosion and weathering process of rock mass are analyzed. And the evolution law of deformation and failure of the marlite mass is studied. The main failure feature of the marlite slope is the karst erosive structure subsidence mode in Three Gorges Reservoir area. The karst erosive structure subsidence mode is explained as follows: the rock mass undergoes the synthetic influence, such as weathering, unloading, corrosion, and so on, many pores and cavities have been formed in the rock mass interior, the rock mass quality is worsen and the rock mass structure is changed, and then the inherent structure of rock mass is collapsed under its gravity, therefore, the failure mode of compaction and subsidence take place. Finally, two examples are used to verify the rock mass parameters in Three Gorges Reservoir area, and the relationship between the marlite slope stability and the time of karst erosion is proposed.

Spallation analysis with a closed trans-scale formulation damage evolution

A closed, trans-scale formulation of damage evolution based on the statistical microdamage mechanics is summarized in this paper. The dynamic function of damage bridges the mesoscopic and macroscopic evolution of damage. The spallation in an aluminium plate is studied with this formulation. It is found that the damage evolution is governed by several dimensionless parameters, i.e., imposed Deborah numbers De* and De, Mach number M and damage number S. In particular, the most critical mode of the macroscopic damage evolution, i.e., the damage localization, is deter-mined by Deborah number De+. Deborah number De* reflects the coupling and competition between the macroscopic loading and the microdamage growth. Therefore, our results reveal the multi-scale nature of spallation. In fact, the damage localization results from the nonlinearity of the microdamage growth. In addition, the dependence of the damage rate on imposed Deborah numbers De* and De, Mach number M and damage number S is discussed.

Damage evolution, localization and failure of solids subjected to impact loading

In order to reveal the underlying mesoscopic mechanism governing the experimentally observed failure in solids subjected to impact loading, this paper presents a model of statistical microdamage evolution to macroscopic failure, in particular to spallation. Based on statistical microdamage mechanics and experimental measurement of nucleation and growth of microcracks in an Al alloy subjected to plate impact loading, the evolution law of damage and the dynamical function of damage are obtained. Then, a lower bound to damage localization can be derived. It is found that the damage evolution beyond the threshold of damage localization is extremely fast. So, damage localization can serve as a precursor to failure. This is supported by experimental observations. On the other hand, the prediction of failure becomes more accurate, when the dynamic function of damage is fitted with longer experimental observations. We also looked at the failure in creep with the same idea. Still, damage localization is a nice precursor to failure in creep rupture.

A micromechanical damage theory for brittle materials with small cracks

For brittle solids containing numerous small cracks, a micromechanical damage theory is presented which accounts for the interactions between different small cracks and the effect of the boundary of a finite solid, and includes growth of the pre-existing small cracks. The analysis is based on a superposition scheme and series expansions of the complex potentials. The small crack evolution process is simulated through the use of fracture mechanics incorporating appropriate failure criteria. The stress-strain relations are obtained from the micromechanics analysis. Typical examples are given to illustrate the potential capability of the proposed theory. These results show that the present method provides a direct and efficient approach to deal with brittle finite solids containing multiple small cracks. The stress-strain relation curves are evaluated for a rectangular plate containing small cracks.