Dynamic damage and failure in viscoplastic materials

In this paper, a dynamic damage model in ductile solids under the application of a dynamic mean tensile stress is developed. The proposed model considers void nucleation and growth as parts of the damage process under intense dynamic loading (strain rates epsilon greater than or equal to 10(3) s(-1)). The evolution equation of the ductile void has the closed form, in which work-hardening behavior, rate-dependent contribution and inertial effects are taken into account. Meanwhile, a plate impact test is performed for simulating the dynamic fracture process in LY12 aluminum alloy. The damage model is incorporated in a hydrodynamic computer code, to simulate the first few stress reverberations in the target as it spalls and postimpact porosity in the specimen. Fair agreement between computed and experimental results is obtained. Numerical analysis shows that the influence of inertial resistance on the initial void growth in the case of high loading rate can not be neglected. It is also indicated that the dynamic growth of voids is highly sensitive to the strain rates.

Embryo-damage induced nucleation of microcracks in an aluminum-alloy under impact loading

The nucleation of microdamage under dynamic loading was investigated through planar impact experiments accomplished with a light gas gun. The microscopic observation of recovered and sectioned specimens showed that microcracks were nucleated only by cracking of brittle particles inside material. However, for comparison the in situ static tensile tests on the same material conducted with a scanning electron microscope showed that the microcracks were nucleated by many forms those were fracture of ductile matrix, debonding particles from matrix and cracking of brittle particles. The quantitative metallographic observations of the specimens subjected to impact loading showed that most of the cracked particles were situated on grain boundaries of the aluminium matrix. These facts suggested the concept of critical size and incubation time of submicroscopic cavities in the dynamic case and the mechanism of embryo-damage induced nucleation by fracture of brittle particles in the aluminium alloy under impact loading was proposed.

Spall damage in aluminum-alloy

A void growth relations for ductile porous materials under intense dynamic general loading condition is presented. The mathematical model includes the influence of inertial effects, material rate sensitivity, as well as the contribution of void surface energy and material work-hardening. Numerical analysis shows that inertia appears to resist the growth of voids. The inertial effects increase quickly with the loading rates. The theoretical analysis suggests that the inertial effects cannot be neglected at high loading rates. Plate-impact tests of aluminum alloy are performed with light gas gun. The processes of dynamic damage in aluminum alloy are successfully simulated with a finite-difference dynamic code in which the theoretical model presented in this paper is incorporated.

Statistical description of pattern evolution in damage-fracture

A model of dynamical process and stochastic jump has been put forward to study the pattern evolution in damage-fracture. According to the final states of evolution processes, the evolution modes can be classified as globally stable modes (GS modes) and evolution induced catastrophic modes (ElC modes); the latter are responsible for fracture. A statistical description is introduced to clarify the pattern evolution in this paper. It is indicated that the appearance of fracture in disordered materials should be depicted by probability distribution function.

Influence of secondary void damage in the matrix material around voids

The mechanism of ductile damage caused by secondary void damage in the matrix around primary voids is studied by large strain, finite element analysis. A cylinder embedding an initially spherical void, a plane stress cell with a circular void and plane strain cell with a cylindrical or a flat void are analysed under different loading conditions. Secondary voids of smaller scale size nucleate in the strain hardening matrix, according to the requirements of some stress/strain criteria. Their growth and coalescence, handled by the empty element technique, demonstrate distinct mechanisms of damage as circumstances change. The macroscopic stress-strain curves are decomposed and illustrated in the form of the deviatoric and the volumetric parts. Concerning the stress response and the void growth prediction, comparisons are made between the present numerical results and those of previous authors. It is shown that loading condition, void growth history and void shape effect incorporated with the interaction between two generations of voids should be accounted for besides the void volume fraction.

Statistical modeling of damage evolution in spallation

In order to understand the mechanism of the incipient spallation in rolled metals, a one dimensional statistical mode1 on evolution of microcracks in spallation was proposed. The crack length appears to be the fundamental variable in the statistical description. Two dynamic processes, crack nucleation and growth, were involved in the model of damage evolution. A simplified case was examined and preliminary correlation to experimental observations of spallation was made.

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.

“Deborah Numbers”, Coupling Multiple Space and Time Scales and Governing Damage Evolution to Failure

Two different spatial levels are involved concerning damage accumulation to eventual failure. nucleation and growth rates of microdamage nN* and V*. It is found that the trans-scale length ratio c*/L does not directly affect the process. Instead, two independent dimensionless numbers: the trans-scale one * * ( V*)including the * **5 * N c V including mesoscopic parameters only, play the key role in the process of damage accumulation to failure. The above implies that there are three time scales involved in the process: the macroscopic imposed time scale tim = /a and two meso-scopic time scales, nucleation and growth of damage, (* *4) N N t =1 n c and tV=c*/V*. Clearly, the dimensionless number De*=tV/tim refers to the ratio of microdamage growth time scale over the macroscopically imposed time scale. So, analogous to the definition of Deborah number as the ratio of relaxation time over external one in rheology. Let De be the imposed Deborah number while De represents the competition and coupling between the microdamage growth and the macroscopically imposed wave loading. In stress-wave induced tensile failure (spallation) De* < 1, this means that microdamage has enough time to grow during the macroscopic wave loading. Thus, the microdamage growth appears to be the predominate mechanism governing the failure. Moreover, the dimensionless number D* = tV/tN characterizes the ratio of two intrinsic mesoscopic time scales: growth over nucleation. Similarly let D be the “intrinsic Deborah number”. Both time scales are relevant to intrinsic relaxation rather than imposed one. Furthermore, the intrinsic Deborah number D* implies a certain characteristic damage. In particular, it is derived that D* is a proper indicator of macroscopic critical damage to damage localization, like D* ∼ (10–3~10–2) in spallation. More importantly, we found that this small intrinsic Deborah number D* indicates the energy partition of microdamage dissipation over bulk plastic work. This explains why spallation can not be formulated by macroscopic energy criterion and must be treated by multi-scale analysis.

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.

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.

固体力学进展及应用: 庆贺李敏华院士90华诞文集

《固体力学进展及应用：庆贺李敏华院士90华诞文集》收录了近代固体力学基础理论及其应用领域的重要科技成果和最新进展。作者是在同体力学领域工作多年的资深研究员，他们来自各行各业，有丰富的科研与丁作经验。他们提供的论文在相当程度上反映当前同体力学的发展现状与成就，并能看出发展趋势，对未来研究的课题选择有参考价值。《固体力学进展及应用：庆贺李敏华院士90华诞文集》还收集了李敏华院士的珍贵照片和纪念李敏华院士90华诞的庆贺和回忆文章，具有重要的史料价值。

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.