微裂纹细观损伤理论

材料的宏细观破坏理论是当前固体力学和材料科学研究的一个重要课题。本文在对连续损伤理论和细观损伤理论进行评述的基础上，着重研究了脆性材料中微裂纹细观损伤问题。本文建立了一套完整的细观损伤理论来分析二维多裂纹体问题。该理论的基本方法是基本解叠加法，此方法直接考虑了微裂纹之间的相互作用以及有限边界的影响。通过叠加原理，使在裂纹面和外边界满足边界条件，用边界配置法化控制方程组为线性方程组，进行数值求解。本文以裂纹密度为参量，针对微裂纹随机分布和平行分布两种情况，计算了无限大体中代表性体元（VRVE）和多裂纹有限体的有效弹性模量。数值计算结果表明，本文所用方法具有统一与直能的优点，采用此法所得模量与试验结果吻合，在处理多裂纹体问题时计算效率高、精度好，对求解多裂纹问题非常有效。此外，通过建立微裂纹晶内扩展准则和穿晶扩展准则，分析了微裂纹扩展连接直至裂纹形成、扩展这一全过程的细观力学行为，对微裂纹的损伤演化过程进行了直接模拟，计算了含微裂纹矩形板的宏观应国变关系曲线。本文进一步提出了三维微裂纹相互作用的数学分析方法 — 扁球坐标和位移函数法，并采用边界配置法或裂纹面面力平均化方法进行求解。数值结果表明，扁球坐标和位移函数法分析三维微裂纹的相互作用问题是有效可行的。最后，本文提出了埋入基体的镶嵌体胞模型，建立了计算非均质体有效弹性模量的解析表达式。该式从理论上讲是严格的，且具有形式简单、内涵丰富及有效弹性模量能显式表达等优点。针对球体含球形夹杂、裂纹及旋转扁球体含球形夹杂、裂纹等不同体胞结构计算了其有效弹性模量，并与其他细观力学方法所得结果进行了比较。本文还将埋入基体的镶嵌体胞模型进行了发展，研究了二相颗粒复合材料的弹塑性本构关系（基体为弹性而颗粒为塑性材料），计算了球体含球形颗粒用旋转扁球体含扁球状颗粒两种体胞结构的宏观应力 - 应变曲线。

Moment Tensor Analysis of the Acoustic Emission Source in the Rock Damage Process

To further investigate the mechanism of acoustic emission (AE) in the rock fracture experiment, moment tensor analysis was carried out. The AE sources characterized by crack sizes, orientations and fracture modes, are represented by a time-dependent momen

Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

A thermodynamic framework for the evolution of damage in concrete under fatigue

A closed-form expression for the dual of dissipation potential is derived within the framework of irreversible thermodynamics using the principles of dimensional analysis and self-similarity. Through this potential, a damage evolution law is proposed for concrete under fatigue loading using the concepts of damage mechanics in conjunction with fracture mechanics. The proposed law is used to compute damage in a volume element when a member is subjected to fatigue loading. The evolution of damage from microcracking to macrocracking of the entire member is captured through a series of volume elements failing one after the other. The number of loading cycles to failure of the member is obtained as the summation of number of cycles to failure for each individual volume element. A parametric study is conducted to determine the effect of the size of the volume element on the model's prediction of fatigue life. A global damage index is also defined, and the residual moment carrying capacity of damaged beams is evaluated. Through a deterministic sensitivity analysis, it is found that the load range and maximum aggregate size are the most influencing parameters on the fatigue life of a plain concrete beam.

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.

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.

Microdamage Evolution, Energy Dissipation and Their Trans-Scale Effects on Macroscopic Failure

The process of damage evolution concerns various scales, from micro- to macroscopic. How to characterize the trans-scale nature of the process is on the challenging frontiers of solid mechanics. In this paper, a closed trans-scale formulation of damage evolution based on statistical microdamage mechanics is presented. As a case study, the damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that the following dimensionless numbers: reduced Mach number M, damage number S, stress wave Fourier number P, intrinsic Deborah number D*, and the imposed Deborah number De*, govern the whole process of deformation and damage evolution. The evaluation of P and the estimation of temperature increase show that the energy equation can be ignored as the first approximation in the case of spallation. Hence, apart from the two conventional macroscopic parameters: the reduced Mach number M and damage number S, the damage evolution in spallation is mainly governed by two microdamage-relevant parameters: the Deborah numbers D* and De*. Higher nucleation and growth rates of microdamage accelerate damage evolution, and result in higher damage in the target plate. In addition, the mere variation in nucleation rate does not change the spatial distribution of damage or form localized rupture, while the increase of microdamage growth rate localizes the damage distribution in the target plate, which can be characterized by the imposed Deborah number De*.

Statistical and stochastic description of microvoid evolution observed in dynamic failure of ductile materials

A brief review is presented of statistical approaches on microdamage evolution. An experimental study of statistical microdamage evolution in two ductile materials under dynamic loading is carried out. The observation indicates that there are large differences in size and distribution of microvoids between these two materials. With this phenomenon in mind, kinetic equations governing the nucleation and growth of microvoids in nonlinear rate-dependent materials are combined with the balance law of void number to establish statistical differential equations that describe the evolution of microvoids' number density. The theoretical solution provides a reasonable explanation of the experimentally observed phenomenon. The effects of stochastic fluctuation which is influenced by the inhomogeneous microscopic structure of materials are subsequently examined (i.e. stochastic growth model). Based on the stochastic differential equation, a Fokker-Planck equation which governs the evolution of the transition probability is derived. The analytical solution for the transition probability is then obtained and the effects of stochastic fluctuation is discussed. The statistical and stochastic analyses may provide effective approaches to reveal the physics of damage evolution and dynamic failure process in ductile materials.

Failure process analysis of cement under explosive loading with a generalized driven nonlinear threshold model

The microstructural heterogeneity and stress fluctuation play important roles in the failure process of brittle materials. In this paper, a generalized driven nonlinear threshold model with stress fluctuation is presented to study the effects of microstructural heterogeneity on continuum damage evolution. As an illustration, the failure process of cement material under explosive loading is analyzed using the model. The result agrees well with the experimental one, which proves the efficiency of the model.

Elastic and statistically brittle (ESB) model, damage localization and catastrophic rupture

In this paper, an elastic and statistically brittle (ESB) model is applied to the process of damage evolution induced catastrophic rupture and the influence of localization and softening on catastrophic rupture is discussed. According to the analysis, the uncertainty of catastrophic rupture should be attributed to the unknown scale of localized zone. Based on the elastic and statistically brittle model but local mean field approximation, the relation between the scale of localized zone and catastrophic rupture is obtained and then justified with experiments. These results can not only give a deeper understanding of the mechanism governing catastrophic rupture, but also provide a possible tool to foresee the occurrence of catastrophic rupture.

MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion

In this work, we introduce the Object Kinetic Monte Carlo (OKMC) simulator MMonCa and simulate the defect evolution in three different materials. We start by explaining the theory of OKMC and showing some details of how such theory is implemented by creating generic structures and algorithms in the objects that we want to simulate. Then we successfully reproduce simulated results for defect evolution in iron, silicon and tungsten using our simulator and compare with available experimental data and similar simulations. The comparisons validate MMonCa showing that it is powerful and flexible enough to be customized and used to study the damage evolution of defects in a wide range of solid materials.

Ultrasonic and X-ray computed tomography characterization of progressive fracture damage in low-porous carbonate rocks

This paper studies the fracturing process in low-porous rocks during uniaxial compressive tests considering the original defects and the new mechanical cracks in the material. For this purpose, five different kinds of rocks have been chosen with carbonate mineralogy and low porosity (lower than 2%). The characterization of the fracture damage is carried out using three different techniques: ultrasounds, mercury porosimetry and X-ray computed tomography. The proposed methodology allows quantifying the evolution of the porous system as well as studying the location of new cracks in the rock samples. Intercrystalline porosity (the smallest pores with pore radius < 1 μm) shows a limited development during loading, disappearing rapidly from the porosimetry curves and it is directly related to the initial plastic behaviour in the stress–strain patterns. However, the biggest pores (corresponding to the cracks) suffer a continuous enlargement until the unstable propagation of fractures. The measured crack initiation stress varies between 0.25 σp and 0.50 σp for marbles and between 0.50 σp and 0.85 σp for micrite limestone. The unstable propagation of cracks is assumed to occur very close to the peak strength. Crack propagation through the sample is completely independent of pre-existing defects (porous bands, stylolites, fractures and veins). The ultrasonic response in the time-domain is less sensitive to the fracture damage than the frequency-domain. P-wave velocity increases during loading test until the beginning of the unstable crack propagation. This increase is higher for marbles (between 15% and 30% from initial vp values) and lower for micrite limestones (between 5% and 10%). When the mechanical cracks propagate unstably, the velocity stops to increase and decreases only when rock damage is very high. Frequency analysis of the ultrasonic signals shows clear changes during the loading process. The spectrum of treated waveforms shows two main frequency peaks centred at low (~ 20 kHz) and high (~ 35 kHz) values. When new fractures appear and grow the amplitude of the high-frequency peak decreases, while that of the low-frequency peak increases. Besides, a slight frequency shift is observed towards higher frequencies.

Investigation of RuO2/Ta2O5 thin film evolution by thermogravimetry combined with mass spectrometry

The thermal evolution process of RuO2–Ta2O5/Ti coatings with varying noble metal content has been investigated under in situ conditions by thermogravimetry combined with mass spectrometry. The gel-like films prepared from alcoholic solutions of the precursor salts (RuCl3·3H2O, TaCl5) onto titanium metal support were heated in an atmosphere containing 20% O2 and 80% Ar up to 600 °C. The evolution of the mixed oxide coatings was followed by the mass spectrometric ion intensity curves. The cracking of retained solvent and the combustion of organic surface species formed were also followed by the mass spectrometric curves. The formation of carbonyl- and carboxylate-type surface species connected to the noble metal was identified by Fourier transform infrared emission spectroscopy. These secondary processes–catalyzed by the noble metal–may play an important role in the development of surface morphology and electrochemical properties. The evolution of the two oxide phases does not take place independently, and the effect of the noble metal as a combustion catalyst was proved.

Response and damage evaluation of reinforced concrete frames subjected to blast loading

Bomb attacks carried out by terrorists, targeting high occupancy buildings, have become increasingly common in recent times. Large numbers of casualties and property damage result from overpressure of the blast followed by failing of structural elements. Understanding the blast response of multi-storey buildings and evaluating their remaining life have therefore become important. Response and damage analysis of single structural components, such as columns or slabs, to explosive loads have been examined in the literature, but the studies on blast response and damage analysis of structural frames in multi-storey buildings is limited and this is necessary for assessing the vulnerability of them. This paper investigates the blast response and damage evaluation of reinforced concrete (RC) frames, designed for normal gravity loads, in order to evaluate their remaining life. Numerical modelling and analysis were carried out using the explicit finite element software, LS DYNA. The modelling and analysis takes into consideration reinforcement details together and material performance under higher strain rates. Damage indices for columns are calculated based on their residual and original capacities. Numerical results generated in the can be used to identify relationships between the blast load parameters and the column damage. Damage index curve will provide a simple means for assessing the damage to a typical multi-storey building RC frame under an external bomb circumstance.