Macroscopic yield criteria for plastic anisotropic materials containing spheroidal voids

The combined effects of void shape and matrix anisotropy on the macroscopic response of ductile porous solids is investigated. The Gologanu–Leblond–Devaux’s (GLD) analysis of an rigid-ideal plastic (von Mises) spheroidal volume containing a confocal spheroidal cavity loaded axisymmetrically is extended to the case when the matrix is anisotropic (obeying Hill’s [Hill, R., 1948. A theory of yielding and plastic flow of anisotropic solids. Proc. Roy. Soc. London A 193, 281–297] anisotropic yield criterion) and the representative volume element is subjected to arbitrary deformation. To derive the overall anisotropic yield criterion, a limit analysis approach is used. Conditions of homogeneous boundary strain rate are imposed on every ellipsoidal confocal with the cavity. A two-field trial velocity satisfying these boundary conditions are considered. It is shown that for cylindrical and spherical void geometries, the proposed criterion reduces to existing anisotropic Gurson-like yield criteria. Furthermore, it is shown that for the case when the matrix is considered isotropic, the new results provide a rigorous generalization to the GLD model. Finally, the accuracy of the proposed approximate yield criterion for plastic anisotropic media containing non-spherical voids is assessed through comparison with numerical results.     

Void growth and coalescence in anisotropic plastic solids

Large strain finite element calculations of unit cells subjected to triaxial axisymmetric loadings are presented for plastically orthotropic materials containing a periodic distribution of aligned spheroidal voids. The spatial distribution of voids and the plastic flow properties of the matrix are assumed to respect transverse isotropy about the axis of symmetry of the imposed loading so that a two-dimensional axisymmetric analysis is adequate. The parameters varied pertain to load triaxiality, matrix anisotropy, initial porosity and initial void shape so as to include the limiting case of penny-shaped cracks. Attention is focussed on comparing the individual and coupled effects of void shape and material anisotropy on the effective stress–strain response and on the evolution of microstructural variables. In addition, the effect of matrix anisotropy on the mode of plastic flow localization is discussed. From the results, two distinct regimes of behavior are identified: (i) at high triaxialities, the effect of material anisotropy is found to be persistent, unlike that of initial void shape and (ii) at moderate triaxialities the influence of void shape is found to depend strongly on matrix anisotropy. The findings are interpreted in light of recent, microscopically informed models of porous metal plasticity. Conversely, observations are made in relation to the relevance of these results in the development and calibration of a broader set of continuum damage mechanics models.     

On ductile fracture initiation toughness: Effects of void volume fraction,void shape and void distribution

This paper studies the effects of the initial relative void spacing, void pattern, void shape and void volume fraction on ductile fracture toughness using three-dimensional, small scale yielding models, where voids are assumed to pre-exist in the material and are explicitly modeled using refined finite elements. Results of this study can be used to explain the observed fracture toughness anisotropy in industrial alloys. Our analyses suggest that simplified models containing a single row of voids ahead of the crack tip is sufficient when the initial void volume fraction remains small. When the initial void volume fraction becomes large, these simplified models can predict the fracture initiation toughness (J_Ic) with adequate accuracy but cannot predict the correct J–R curve because they over-predict the interaction among growing voids on the plane of crack propagation. Consequently, finite element models containing multiple rows of voids should be used when the material has large initial void volume fraction.     

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.     

RVE-based studies on the coupled effects of void size and void shape on yield behavior and void growth at micron scales

The present paper extends the Gurson and GLD models [Gurson, A.L., 1977. Continuum theory of ductile rupture by void nucleation and growth, Part I—yield criteria and flow rules for porous ductile media. J. Mech. Phys. Solids 99, 2–15; Gologanu, M., Leblond, J.B., Devaux, J., 1993. Approximate models for ductile metals containing non-spherical voids—case of axisymmetric prolate ellipsoidal cavities. J. Mech. Phys. Solids 41, 1723–1754; Gologanu, M., Leblond, J.B., Devaux, J., 1994. Approximate models for ductile metals containing non-spherical voids—case of axisymmetric oblate ellipsoidal cavities. J. Eng. Mater. Technol. 116, 290–297] to involve the coupled effects of void size and void shape on the macroscopic yield behavior of non-linear porous materials and on the void growth. A spheroidal representative volume element (RVE) under a remote axisymmetric homogenous strain boundary condition is carefully analyzed. A wide range of void aspect ratios covering the oblate spheroidal, spherical and prolate spheroidal void are taken into account to reflect the shape effect. The size effect is captured by the Fleck–Hutchinson phenomenological strain gradient plasticity theory [Fleck, N.A., Hutchinson, J.W., 1997. Strain gradient plasticity. In: Hutchinson, J.W., Wu, T.Y. (Eds.), Advance in Applied Mechanics, vol. 33, Academic Press, New York, pp. 295–361]. A new size-dependent damage model like the Gurson and GLD models is developed based on the traditional minimum plasticity potential principle. Consequently, the coupled effects of void size and void shape on yield behavior of porous materials and void growth are discussed in detail. The results indicate that the void shape effect on the yield behavior of porous materials and on the void growth can be modified dramatically by the void size effect and vice versa. The applied stress triaxiality plays an important role in these coupled effects. Moreover, there exists a cut-off void radius r_c, which depends only on the intrinsic length l₁ associated with the stretch strain gradient. Voids of effective radius smaller than the critical radius r_c are less susceptible to grow. These findings are helpful to our further understanding to some impenetrable micrographs of the ductile fracture surfaces.     

Verification of the transferability of micromechanical parameters by cell model calculations with visco-plastic materials

Finite element (FE) calculations of a cylindrical cell containing a spherical hole have been performed under large strain conditions for varying triaxiality with three different constitutive models for the matrix material, i.e. rate independent plastic material with isotropic hardening, visco-plastic material under both isothermal and adiabatic conditions, and porous plastic material with a second population of voids nucleating strain controlled. The “mesoscopic” stress-strain and void growth responses of the cell are compared with predictions of the modified Gurson model in order to study the effects of varying triaxiality and strain rate on the critical void volume fraction. The interaction of two different sizes of voids was modelled by changing the strain level for nucleation and the stress triaxiality. The study confirms that the void volume fraction at void coalescence does not depend significantly on the triaxiality if the initial volume fraction of the primary voids is small and if there are no secondary voids. The strain rate does not affect f_c either. The results also indicate that a single internal variable, f, is not sufficient to characterize the fracture processes in materials containing two different size-scales of void nucleating particles.     

Plastic potential and yield function of porous materials with aligned and randomly oriented spheroidal voids

Based on an energy approach, the plastic potential and yield function of a porous material containing either aligned or randomly oriented spheroidal voids are developed at a given porosity and pore shape. The theory is applicable to both elastically compressible and incompressible matrix and, it is proved that, in the incompressible case, the theory with spherical and aligned spheroidal voids also coincides with Ponte Castaneda's bounds of the Hashin-Shtrikman and Willis types, respectively. Comparison is also made between the present theory and those of Gurson and Tvergaard, with a result giving strong overall support of this new development. For the influence of pore shape, the yield function and therefore the stress-strain curve of the isotropic porous material are found to be stiffest when the voids are spherical, and those associated with other pore shapes all fall below these values, the weakest one being caused by the disc-shaped voids. The transversely isotropic nature of the yield function and stress-strain curves of a porous material containing aligned pores are also demonstrated as a function of porosity and pore shape, and it is further substantiated with a comparison with an exact, local analysis when the void shape becomes cylindrical.     

Modeling of ductile fracture: Application of the mechanism-based concepts

This paper summarizes our recent studies on modeling ductile fracture in structural materials using the mechanism-based concepts. We describe two numerical approaches to model the material failure process by void growth and coalescence. In the first approach, voids are considered explicitly and modeled using refined finite elements. In order to predict crack initiation and propagation, a void coalescence criterion is established by conducting a series of systematic finite element analyses of the void-containing, representative material volume (RMV) subjected to different macroscopic stress states and expressed as a function of the stress triaxiality ratio and the Lode angle. The discrete void approach provides a straightforward way for studying the effects of microstructure on fracture toughness. In the second approach, the void-containing material is considered as a homogenized continuum governed by porous plasticity models. This makes it possible to simulate large amount of crack extension because only one element is needed for a representative material volume. As an example, a numerical approach is proposed to predict ductile crack growth in thin panels of a 2024-T3 aluminum alloy, where a modified Gologanu–Leblond–Devaux model [Gologanu, M., Leblond, J.B., Devaux, J., 1993. Approximate models for ductile metals containing nonspherical voids – Case of axisymmetric prolate ellipsoidal cavities. J. Mech. Phys. Solids 41, 1723–1754; Gologanu, M., Leblond, J.B., Devaux, J., 1994. Approximate models for ductile metals containing nonspherical voids – Case of axisymmetric oblate ellipsoidal cavities. J. Eng. Mater. Tech. 116, 290–297; Gologanu, M., Leblond, J.B., Perrin, G., Devaux, J., 1995. Recent extensions of Gurson’s model for porous ductile metals. In: Suquet, P. (Ed.) Continuum Micromechanics. Springer-Verlag, pp. 61–130] is used to describe the evolution of void shape and void volume fraction and the associated material softening, and the material failure criterion is calibrated using experimental data. The calibrated computational model successfully predicts crack extension in various fracture specimens, including the compact tension specimen, middle crack tension specimens, multi-site damage specimens and the pressurized cylindrical shell specimen.     

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.     

Void growth and coalescence in single crystals

Void growth and coalescence in single crystals are investigated using crystal plasticity based 3D finite element calculations. A unit cell involving a single spherical void and fully periodic boundary conditions is deformed under constant macroscopic stress triaxiality. Simulations are performed for different values of the stress triaxiality, for different crystal orientations, and for low and high work-hardening capacity. Under low stress triaxiality, the void shape evolution, void growth, and strain at the onset of coalescence are strongly dependent on the crystal orientation, while under high stress triaxiality, only the void growth rate is affected by the crystal orientation. These effects lead to significant variations in the ductility defined as the strain at the onset of coalescence. An attempt is made to predict the onset of coalescence using two different versions of the Thomason void coalescence criterion, initially developed in the framework of isotropic perfect plasticity. The first version is based on a mean effective yield stress of the matrix and involves a fitting parameter to properly take into account material strain hardening. The second version of the Thomason criterion is based on a local value of the effective yield stress in the ligament between the voids, with no fitting parameter. The first version is accurate to within 20% relative error for most cases, and often more accurate. The second version provides the same level of accuracy except for one crystal orientation. Such a predictive coalescence criterion constitutes an important ingredient towards the development of a full constitutive model for porous single crystals.     

不同应力三维度条件下孔洞的演变及修正的Gurson模型

本文对含不同形状孔洞的幂硬化材料的圆柱体胞模型，运用控制宏观应力三维工的方法进行了有限元分析。计算结果表明：１．孔洞初始形状，应力三维度对孔洞的长大有重要影响；２．Ｇｕｓｏｎ模型对孔洞长大规律的描述是不准确的，不准确度与孔洞初始形状，应力三维度有关，修正后的Ｇｕｒｓｏｎ模型与有限元结果吻合较好；３．在低应力三维度区，孔洞以及形状改变为主，在高应力三维度区，孔洞以扩张为主；     

考虑三轴约束时孔洞的聚合机理及有效能量准则

通过体胞分析方法,对不同状孔洞在从光滑试样到裂纹试样的三轴应力场中的聚合机理进行了较精解的有限元分析,计算结果表明：（１）孔洞的相互靠近和横向扩展是导致相邻孔洞发生内颈缩聚合的两种基本机制,在应力三维度Ｒσ等于１．２５附近,这两种机制发生较明显的变化。（２）单纯以孔洞体积分数ｆＣ概念为基础的材料破坏参数一般敏感于应力三维度,不能很好地预报不同三轴应力场中材料的破坏,在此基础上,提出了描述孔洞聚合的     

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.     

A criterion for the onset of void coalescence under combined tension and shear

Depending on the relative positions of voids and on the loading conditions, shear loading components can play an important role in the void coalescence process leading to ductile fracture. Yet, most void coalescence criteria including the original criterion of Thomason, and its various extensions/improvements, take only normal loads into account and neglect the contribution from shear loads to coalescence. Shear can affect both the stress/strain at the onset of coalescence and the direction of deformation localization. In this paper, first, the predictive capabilities of different coalescence criteria without shear effect are critically assessed and the expressions involved in the original Thomason criterion are fine-tuned by comparing with 3D finite element calculations performed on a unit cell containing a spheroidal void. Then, the improved Thomason criterion is theoretically extended—by using limit load analysis—to incorporate the effect of shear. The predictions of this new coalescence criterion are in good agreement with the results produced by 3D finite element calculations, for both loadings involving or not a shear component.     

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.