Dynamics of void collapse in shocked energetic materials: physics of void–void interactions

This work presents the response of a porous energetic material subjected to severe transient loading conditions. The porosities, represented by voids, entirely change the response of an otherwise homogeneous material. The variations in terms of energy distribution and maximum temperature reached in the material in the presence of heterogeneities (voids) but in the absence of chemical reactions are studied. This study also accounts for void–void interactions to enhance the understanding of the localization of energy in the material. It is observed that relative position of voids can have important consequence on energy distribution as well as rise in temperature of the energetic material. The relative position of voids further influences the interaction of secondary shock waves generated during the collapse of one void with the downstream voids. This interaction can either enhance or diminish the strength of the shock depending on the location of downstream voids. This work also reveals that the findings from mutual void–void interactions can be used to study systems with multiple voids. This is shown by analyzing systems with 10–25 % void volume fraction. The effect of void–void interactions are connected to the overall response of a chemically inert porous material to imposed transient loads.     

Void growth to coalescence in a non-local material

The size-effect in metals containing distributed spherical voids is analyzed numerically using a finite strain generalization of a length scale dependent plasticity theory. Results are obtained for stress-triaxialities relevant in front of a crack tip in an elastic-plastic metal. The influence of different material length parameters in a multi-parameter theory is studied, and it is shown that the important length parameter is the same as under purely hydrostatic loading. It is quantified how micron scale voids grow less rapidly than larger voids, and the implications of this in the overall strength of the material is emphasized. The size effect on the onset of coalescence is studied, and results for the void volume fraction and the strain at the onset of coalescence are presented. It is concluded that for cracked specimens not only the void volume fraction, but also the typical void size is of importance to the fracture strength of ductile materials.     

Nonlocal plasticity effects on interaction of different size voids

A nonlocal elastic–plastic material model is used to show that the rate of void growth is significantly reduced when the voids are small enough to be comparable with a characteristic material length. For a very small void in the material between much larger voids the competition between an increased growth rate due to the stress concentrations around the larger voids and a reduced growth rate due to the nonlocal effects is studied. The analyses are based on an axisymmetric unit cell model with special boundary conditions, which allow for a relatively simple investigation of a full three dimensional array of spherical voids. It is shown that the high growth rate of very small voids predicted by conventional plasticity theory is not realistic when the effect of a characteristic length, dependent on the dislocation structure, is accounted for.     

Modeling of ductile fracture: Significance of void coalescence

In this paper void coalescence is regarded as the result of localization of plastic flow between enlarged voids. We obtain the failure criterion for a representative material volume (RMV) in terms of the macroscopic equivalent strain (E_c) as a function of the stress triaxiality parameter (T) and the Lode angle (θ) by conducting systematic finite element analyses of the void-containing RMV subjected to different macroscopic stress states. A series of parameter studies are conducted to examine the effects of the initial shape and volume fraction of the primary void and nucleation, growth, and coalescence of secondary voids on the predicted failure surface E_c(T, θ). As an application, a numerical approach is proposed to predict ductile crack growth in thin panels of a 2024-T3 aluminum alloy, where a porous plasticity model is used to describe the void growth process and the expression for E_c is calibrated using experimental data. The calibrated computational model is applied to predict crack extension in fracture specimens having various initial crack configurations and the numerical predictions agree very well with experimental measurements.     

Stress Magnification due to Stretching and Bending of Thin Ligaments between Voids

Stress magnification in thin ligaments between small and large cylindrical voids is obtained by matching the inner field approximation by beam theory to the outer rigid-body field in the bulk of the material. A void between two larger voids is modeled as a large hole within a strip of straight edges (boundaries of the holes with infinite radii of curvature). Both stretching and bending types of loading are applied to the strip. Comparison of different orders of stress magnification for different geometries and loading conditions is made. It is shown that the order of stress magnification in thin ligaments is (R/δ) ⁿ , where n=1/2 in the ligament between one small and one large void, n=1 in the ligament between one small void and two large voids, or between two small and two large voids, and n=2 in the ligament between a large void and a small void coalescing with another large void. The relevance of these results for the study of material failure by void growth and coalescence is discussed.     

The modified Gurson model accounting for the void size effect

The Gurson model [J. Engrg. Mater. Technol. 99 (1977) 2] has been widely used to study the deformation and failure of metallic materials containing microvoids. The void volume fraction is the only parameter representing voids since the void size does not come into play in the Gurson model. Based on the Taylor dislocation model [Proc. R. Soc. (Lond.) A145 (1934) 362; J. Int. Metals 62 (1938) 307], we extend the Gurson model to account for the void size effect. It is shown that the yield surfaces for micron- and submicron-sized voids are significantly larger than that given by the Gurson model. For a voided, dilating material subject to uniaxial tension, the void size has essentially no effect on the stress–strain curve at small initial void volume fraction. However, as the initial void volume fraction increases, the void size effect may become significant.     

Generalized variational principles of the viscoelastic body with voids and their applications

From the Boltzmann‘ s constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given. By using the variational integral method, the convolution-type functional was given and the corresponding generalized variational principles and potential energy principle of viscoelastic solids with voids were presented. It can be shown that the variational principles correspond to the differential equations and theinitial and boundary conditions of viscoelastic body with voids. As an application, a generalized variational principle of viscoelastic Timoshenko beams with damage was obtained which corresponds to the differential equations of generalized motion and the initial and boundary conditions of beams. The variational principles provide a way for solving problems of viscoelastic solids with voids.     

Modeling of crack growth in ductile solids: a three-dimensional analysis

A population of several spherical voids is included in a three-dimensional, small scale yielding model. Two distinct void growth mechanisms, put forth by [Int. J. Solids Struct. 39 (2002) 3581] for the case of a two-dimensional model containing cylindrical voids, are well contained in the model developed in this study for spherical voids. A material failure criterion, based on the occurrence of void coalescence in the unit cell model, is established. The critical ligament reduction ratio, which varies with stress triaxiality and initial porosity, is used to determine ligament failure between the crack tip and the nearest void. A comparison of crack initiation toughness of the model containing cylindrical voids with the model containing spherical voids reveals that the material having a sizeable fraction of spherical voids is tougher than the material having cylindrical voids. The proposed material failure determination method is then used to establish the fracture resistance curve (J–R curve) of the material. For a ductile material containing a small volume fraction of microscopic voids initially, the void by void growth mechanism prevails, which results in a J–R curve having steep slope. On the other hand, for a ductile material containing a large volume fraction of initial voids, the multiple voids interaction mechanism prevails, which results in a flat J–R curve. Next, the effect of T-stress on fracture resistance is examined. Finally, nucleation and growth of secondary microvoids and their effects on void coalescence are briefly discussed.     

A phase field model for failure in interconnect lines due to coupled diffusion mechanisms

We describe a diffuse interface, or phase field model for simulating electromigration and stress-induced void evolution and growth in interconnect lines. Microstructural evolution is tracked by defining an order parameter, which takes on distinct uniform values within solid material and voids, and varying rapidly from one to the other over narrow interfacial layers associated with the void surfaces. The order parameter is governed by a form of the Cahn-Hilliard equation. An asymptotic analysis demonstrates that the zero contour of order parameter tracks the motion of a void evolving by coupled surface and lattice diffusion, driven by stress, electron wind and vacancy concentration gradients. Efficient finite element schemes are described to solve the modified Cahn-Hilliard equation, as well as the equations associated with the accompanying mechanical, electrical and bulk diffusion problems. The accuracy and convergence of the numerical scheme is investigated by comparing results to known analytical solutions. The method is applied to solve various problems involving void growth and evolution in representative interconnect geometries.     

EVOLUTION OF VOIDS IN DUCTILE POROUS MATERIALS AT HIGH STRAIN RATE

In the present work,a dynamic damage model in ductile materials underthe application of dynamic general stresses loading is presented.The evolution equationof ductile voids has the closed form,in which work-hardening,the change of surfaceenergy of voids,rate-dependent,inertial effects are taken into account.Theexpressions of critical stresses for the growth and compaction of voids are directlyobtained from the evolution equations of voids.From the expressions,the resultobtained by Carroll and Holt,as a special example,is given.Numerical analysis ofthe model indicates that the growth of voids is sensitive to the strain rates.The voidsgrow quickly as the increase of strain rates.It is also shown that the influence of theinertial effects on the void growth is great at high loading rates.It appears to resist thegrowth of voids.In addition,a dynamic collapse model of ductile voids is alsoproposed,which can be applied to study the problems of compaction in powder andother materials.     

Calibration,validation, and verification including uncertainty of a physically motivated internal state variable plasticity and damage model

In this paper, we illustrate a formal calibration, validation, and verification process that includes uncertainty in an internal state variable plasticity-damage model that is implemented in a finite element code. The physically motivated continuum model characterizes damage evolution by incorporating material uncertainty due to microstructural spatial clustering. The uncertainty analysis is performed by introducing material variation through model validation and verification. The effect of variability in microstructural clustering and boundary conditions on the sensitivities and uncertainty of the plasticity-damage evolution for the 7075 aluminum alloy is characterized. The results show the potential of this methodology in the evaluation of material response uncertainty due to microstructure spatial clustering and its effect on damage evolution. For damage evolution, we have shown that the initial isotropic damage evolved into an anisotropic form as the deformation increased which is consistent with experimentally observed behavior for 7075 aluminum alloy in literature. Also, the sensitivities were found to be consistent with the physics of damage progression for this particular type of material. Through the sensitivity analysis, the initial defect size and number density of cracked particles are important at the beginning of deformation. As damage evolves, more voids are nucleated and grow and the sensitivity analysis illustrates this as well. Then, voids combine with each other and coalescence becomes the main driver, which is also confirmed by the sensitivity analysis. This work also shows that the microstructurally based damage evolution equations provide an accurate representation of the damage progression due to large intermetallic particles. Finally, we illustrate that the initial variation in the microstructure clustering can lead to about ±7.0%, ±8.1%, and ±9.75% variation in the elongation to failure strain for torsion, tensile, and compressive loading, respectively.     

Bounds and estimates for the effective yield surface of porous media with a uniform or a nonuniform distribution of voids

This paper aims at studying the effects of a nonuniform distribution of voids on the macroscopic yield response of porous media with a rigid-perfectly plastic matrix. For this purpose, a semi-analytical model, recently proposed by Bilger et al. [Bilger, N., Auslender, F., Bornert, M., Masson, R., 2002. New bounds and estimates for porous media with rigid perfectly plastic matrix. C. R. Mecanique 330, 127–132], is extended to more general situations where the local porosity can fluctuate. The microstructure is described by a generalized Hashin-type assemblage of hollow spheres and the distribution of the local porosity is obtained from a three-dimensional simulated microstructure. The matrix layer around the voids is discretized into concentric sub-layers so as to take better into account the plasticity gradient along the radial direction. Classical homogenization techniques then provide new self-consistent estimates and upper bounds for the macroscopic yield surface. These results are compared first to the predictions of the Gurson model and its extensions and then to numerical results derived from three-dimensional Fast Fourier Transform (FFT) calculations carried out with the same material porosity distribution. A good agreement is obtained with the three-dimensional FFT calculations and with Gurson–Tvergaard's predictions even for high triaxiality and without fitting any parameter. Nevertheless, when the heterogeneous distribution of voids tends to form clusters, the proposed model fails to capture the properties of the macroscopic yield surface for large triaxiality factors.     

Ductile fracture by cavity nucleation between larger voids

A mechanism of ductile fracture involving the interaction of relatively large voids with small-scale voids is studied by a computational model. The larger voids are described as circular cylindrical holes arranged in a doubly periodic array in the initial state. In the matrix material between these voids the nucleation and growth of much smaller voids is accounted for by using approximate constitutive equations for a ductile, porous medium. The computations show bands of highly localized straining and void growth, initiating at the surfaces of larger voids and growing into the matrix material, until the bands connect two neighbouring voids. The materials are analysed both under plane strain conditions and under conditions approximating those in a round tensile bar. The failure strains obtained under different principal stress ratios show rather good agreement when plotted against a measure of the stress-triaxiality.     

Relations between a micro-mechanical model and a damage model for ductile failure in shear

Gurson type constitutive models that account for void growth to coalescence are not able to describe ductile fracture in simple shear, where there is no hydrostatic tension in the material. But recent micro-mechanical studies have shown that in shear the voids are flattened out to micro-cracks, which rotate and elongate until interaction with neighbouring micro-cracks gives coalescence. Thus, the failure mechanism is very different from that under tensile loading. Also, the Gurson model has recently been extended to describe failure in shear, by adding a damage term to the expression for the growth of the void volume fraction, and it has been shown that this extended model can represent experimental observations. Here, numerical studies are carried out to compare predictions of the shear-extended Gurson model with the shear failures predicted by the micro-mechanical cell model. Both models show a strong dependence on the level of hydrostatic tension. Even though the reason for this pressure dependence is different in the two models, as the shear-extended Gurson model does not describe voids flattening out and the associated failure mechanism by micro-cracks interacting with neighbouring micro-cracks, it is shown that the trends of the predictions are in good agreement.     

The mechanical behaviour of ductile materials damaged by two generations of voids

Large strain finite element method is employed to investigate the damaging effect of two generations of voids in ductile materials. An axisymmetric cylinder embedding an initially spherical void is chosen as the model cell. Secondary voids will initiate around the initial void when the local stress/strain in the matrix increases to certain critical conditions. This event is numerically simulated through an empty element technique. The interaction between these two generations of voids has been proved to be favourable to the voiding condition, thus accelerating the material damage, characterized by the value of the overall elastic modulus which may undergo drastic drop when nearing final fracture.     

多孔Mooney-Rivlin材料矩形板的单向拉伸

本文利用不同可压超弹性材料大变形的Mooney-Rivlin应变能函数研究了含有多个微孔的矩形板在单向拉伸作用下的有限变形和受力分析。首先利用不可压条件得到了文中所给的含有某种对称性分布的多个微孔的矩形板的变形模式函数，其中所含的一个参数可由远离微孔的无穷远处的变形状态确定，另一个参可用最小势能原理导出变分近似解。文中详细分析了板中微孔（一个，三个和五个）随载荷作用的增长情况和微孔边缘应力的分布情况，并进行了比较。讨论了微孔的个数和排列方式，微孔的孔间距离等因素对微孔增长和应力分布的影响。     

Prediction of ductile fracture in tension by bifurcation,localization, and imperfection analyses

In this investigation, it is shown that the onset of ductile fracture in tension can be interpreted as the result of a supercritical bifurcation of homogeneous deformation and that this fact can be applied to predict ductile fracture initiation of materials with general imperfections or flaws. We focus on one dimensional quasi-static simple tension for rate-independent isotropic plastic materials. For deformation beyond the bifurcation point, multiple equilibrium paths appear. The homogeneous deformation, as one of the equilibrium paths, loses stability while the inhomogeneous paths are stable, thus indicating the occurrence of strain localization. This investigation also provides a physical example for the application of the Lambert W function in material localization analyses. Material instability is treated as the instability of a static system with dynamic perturbation. We also address the presence of microvoids in a power law plastic material as an unfolding of the supercritical pitchfork bifurcation. The imperfect system, idealized as spherical voids within the plastic matrix, is analyzed using the familiar Gurson model which is based on the presumption of a randomly voided material and characterized by the volume fraction of voids. If, in addition, the sizes of the microvoids are known, this then provides a length scale for the imperfection zone. In this manner, relevance to the sample size effects of strain-to-failure for ductile fracture initiation is addressed by considering separate zones with variations in void volume fractions. Fracture initiation predictions are presented and compare very well to existing experimental results.     

受有两级空洞损伤时韧性材料的力学行为

本文利用大应变有限元方法研究了两级空洞对韧性材料的损伤作用.模型是以轴对称圆柱基体作为胞元,内含一初始的球型空洞.基体内的应力/应变随胞元外载的增大而达到临界状态,从而在围绕初级空洞的基体内将萌生次级空洞.后者是由空单元实现的.两级空洞的交互作用被证明将促进材料中的空洞化现象从而加速损伤并导至材料的总体弹性模量值在临近破断时急剧下降.     

A physical model for nucleation and early growth of voids in ductile materials under dynamic loading

Spall fracture and other rapid tensile failures in ductile materials are often dominated by the rapid growth of voids. Recent research on the mechanics of void growth clearly shows that void nucleation may be represented as a bifurcation phenomenon, wherein a void forms spontaneously followed by highly localized plastic flow around the new void. Although thermal, viscoplastic, and work hardening effects all play an essential role in the earliest stages of nucleation and growth, the flow becomes dominated by spherical radial inertia, which soon causes all voids to grow asymptotically at the same rate, regardless of differences in initial conditions or constitutive details, provided only that there is the same density of matrix material and the same excess loading history beyond the cavitation stress.These two facts, initiation by bifurcation at a cavitation stress, at which a void first appears, and rapid domination by inertia, are used to postulate a simple, but physically realistic, model for nucleation and early growth of voids in a ductile material under rapid tensile loading. A reasonable statistical distribution for the cavitation stress at various nucleation sites and a simple similarity solution for inertially dominated void growth permit a simple calculation of the initiation and early growth of porosity in the material.Parametric analyses are presented to show the effect that loading rate, peak loading stress, density of nucleation sites, physical properties of the material, etc. have on the applied pressure and distribution of void sizes when a critical porosity is reached.     

Propagation of wave at the boundary surface of transversely isotropic thermoelastic material with voids and isotropic elastic half-space

The purpose of this research is to study the effect of voids on the surface wave propagation in a layer of a transversely isotropic thermoelastic material with voids lying over an isotropic elastic half-space. The frequency equation is derived after developing a mathematical model for welded and smooth contact boundary conditions. The dispersion curves giving the phase velocity and attenuation coefficient via wave number are plotted graphically to depict the effects of voids and anisotropy for welded contact boundary conditions. The specific loss and amplitudes of the volume fraction field, the normal stress, and the temperature change for welded contact are obtained and shown graphically for a particular model to depict the voids and anisotropy effects. Some special cases are also deduced from the present investigation.