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.     

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.     

The size effect on void growth in ductile materials

We have extended the Rice-Tracey model (J. Mech. Phys. Solids 17 (1969) 201) of void growth to account for the void size effect based on the Taylor dislocation model, and have found that small voids tend to grow slower than large voids. For a perfectly plastic solid, the void size effect comes into play through the ratio εl/R₀, where l is the intrinsic material length on the order of microns, ε the remote effective strain, and R₀ the void size. For micron-sized voids and small remote effective strain such that εl/R₀?0.02, the void size influences the void growth rate only at high stress triaxialities. However, for sub-micron-sized voids and relatively large effective strain such that εl/R₀>0.2, the void size has a significant effect on the void growth rate at all levels of stress triaxiality. We have also obtained the asymptotic solutions of void growth rate at high stress triaxialities accounting for the void size effect. For εl/R₀>0.2, the void growth rate scales with the square of mean stress, rather than the exponential function in the Rice-Tracey model (1969). The void size effect in a power-law hardening solid has also been studied.     

Influence of stress triaxiality and temperature on void growth and ductile/brittle transition behavior of 40 Cr steel

Examined experimentally are the influence of stress triaxiality and temperature on the growth of microvoids and the ductile/brittle transition (DBT) macrobehavior of 40 Cr steel subjected to two different heat treatments. This is accomplished by testing more than 300 smooth and notched specimens over a temperature range of 20°C to −196°C. Changes in the microstructure morphology are examined by scanning electron microscopy (SEM) and identified with fracture data on a surface constructed from the uniaxial strain ε_c at fracture, the stress triaxiality R_σ and the temperature T. While stress triaxiality has a significant influence on the DBT temperature T_c, it does not affect the ratio of the average radius of voids R_o to that of inclusions R_i. The ratio R_o/R_i is found to increase with temperature and remains constant in specimens with different notch radii regardless of the temperature. Empirical relations between T_c and R_σ⁻ and R_o/R_i and T are proposed to better understand how macrofracture parameters are influenced by microstructure entities.     

Influence of stress triaxiality on damage and crack tip opening displacement parameters for steels

The dependence of the void growth parameter on the local stress triaxiality and local effective plastic strain near the crack tip of ductile materials provides the motivation to seek for parameters that could rank the ductility of steels. Experimental data for AS 1405-180, AS 1204-350, HY-80 and C---Mn steels show that the crack tip opening displacement (CTOD) at initiation δ_c decreases with increasing crack tip stress triaxiality. This trend is confirmed by analysis. As the critical local effective plastic strain ε_ec also decreases with increasing local stress triaxiality, the ratio δ_c/ε_ec is found to remain nearly constant or independent of the local constraint, i.e., the stress triaxiality. These parameters are given for a class of steels in this paper. Their association with the resistance to ductile fracture remains to be investigated.     

Influences of particle size and interface energy on the stress concentration induced by the oblate spheroidal particle and the void nucleation mechanism

Separation of the particle–matrix interface and breakage of the second-phase particle are two main void nucleation mechanisms, which are directly associated with the stress concentration factors (SCFs) at the interface and within the particle, respectively. This work investigates the coupled effects of particle size and particle shape on these stress concentrations by solving an infinite solid containing an oblate spheroidal particle under remote stress boundary condition. The phenomenological strain plasticity theory by Fleck–Hutchinson [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] is adopted to capture the size effect, various particle aspect ratios are considered to depict the particle shape effect and an interfacial energy concept is introduced to settle the double-traction equilibrium problem at the matrix–particle interface. By using a Ritz procedure, solutions about the stress concentrations are numerically achieved and three main results are found. First, the interfacial normal stress near the particle pole, the interfacial shear stress and the particle opening stress are dramatically elevated and their distributions are significantly modified by decrease in the particle size. Second, this particle size effect is influenced by the remote effective strain, remote stress triaxiality and the interfacial energy to different extent. Finally, the particle shape effect is coupled with this particle size effect, and the more oblate the particle is, the more significant the size effect on SCF elevation is. These findings are helpful for us to understand deeply the void nucleation mechanism at the micron scale.     

Micromechanics of coalescence in ductile fracture

Significant progress has been recently made in modelling the onset of void coalescence by internal necking in ductile materials. The aim of this paper is to develop a micro-mechanical framework for the whole coalescence regime, suitable for finite-element implementation. The model is defined by a set of constitutive equations including a closed form of the yield surface along with appropriate evolution laws for void shape and ligament size. Normality is still obeyed during coalescence. The derivation of the evolution laws is carefully guided by coalescence phenomenology inferred from micromechanical unit-cell calculations. The major implication of the model is that the stress carrying capacity of the elementary volume vanishes as a natural outcome of ligament size reduction. Moreover, the drop in the macroscopic stress accompanying coalescence can be quantified for many initial microstructures provided that the microstructure state is known at incipient coalescence. The second part of the paper addresses a more practical issue, that is the prediction of the acceleration rate δ in the Tvergaard-Needleman phenomenological approach to coalescence. For that purpose, a Gurson-like model including void shape effects is used. Results are presented and discussed in the limiting case of a non-hardening material for different initial microstructures and various stress states. Predicted values of δ are extremely sensitive to stress triaxiality and initial spacing ratio. The effect of initial porosity is significant at low triaxiality whereas the effect of initial void shape is emphasized at high triaxiality.     

Ductile fracture of materials with high void volumefraction

In the present paper, axisymmetric cell models containing one or two voids and athree-dimensional cell model containing two voids have been used to investigate void size andspacing effect on the ductile fracture in materials with high initial void volume fraction. They areperformed for round smooth and round notched specimens under uniaxial tension. The examplematerial used for comparison is a nodular cast iron material GGG-40 with initial void volumefraction of 7.7%. The parameters considered in this paper are void size and shape foraxisymmetric cell models containing a single void, and void distribution pattern foraxisymmetric and 3D cell models containing two voids of different sizes. The results obtainedfrom these cell models by using FEM calculations are compared with the Gurson model, theGurson–Tvergaard–Needleman model, the Rice–Tracey model and the modified Rice–Traceymodel. It can be stated that the influence of void size and void spacing on the growth in volumeof voids is very large, and it is dependent on the distribution of voids. Using non-uniform voiddistribution, the results of axisymmetric cell models can explain how a void can grow in anunstable state under very low stress triaxiality at very small strain as observed in experiments.Calculations using cell models containing two voids give very different results about the stableand unstable growth of voids which are strongly dependent on the configuration of cell model.     

The criterion of critical void growth ratio and its applications

In the present paper, the new material parameter of ductile fracturev _gc with an obvious microscopic physical meaning is emphasized, which can be determined by macroscopic method. The theoretical and experimental background for the criterion of critical void growth ratio is expounded from both macroscopic and microscopic points of view. The relationship between the macroform parameterv _gc and microform oneR _c /R _go of critical void growth ratio is investigated. The engineering experimental procedure for determining the value ofv _gc of materials has been recommended. The values ofv _gc for some materials with different heat treatments have been estimated. It is found that the parameterv _gc can be effectively used for selecting materials and their heat treatments. A modified ASPEF criterion accounting for stress states has been proposed. Related expression of modified ASPEF andv _gc has been derived.     

The effect of straining mode on void growth

Large strain finite element method is employed to investigate the effect of straining mode on void growth. Axisymmetric cell model embedded with spherical void is controlled by constant triaxiality loading, while plane-stress model containing a circular void is loaded by constant ratio of straining. Elastic-plastic material is used for the matrix in both cases. It is concluded that, besides the known effect of triaxiality, the straining mode which intensifies the plastic concentration around the void is also a void growth stimulator. Experimental results are cited to justify the computation results. This paper is jointly supported by the National Natural Science Foundation of China (19872064), the Chinese Academy of Sciences (KJ951-1-201) and the Laboratory for Nonlinear mechanics of Continuous Media of the Institute of Mechanics     

Atomistic insights into dislocation-based mechanisms of void growth and coalescence

One of the low-temperature failure mechanisms in ductile metallic alloys is the growth of voids and their coalescence. In the present work we attempt to obtain atomistic insights into the mechanisms underpinning cavitation in a representative metal, namely Aluminum. Often the pre-existing voids in metallic alloys such as Al have complex shapes (e.g. corrosion pits) and the defromation/damage mechanisms exhibit a rich size-dependent behavior across various material length scales. We focus on these two issues in this paper through large-scale calculations on specimens of sizes ranging from 18 thousand to 1.08 million atoms. In addition to the elucidation of the dislocation propagation based void growth mechanism we highlight the observed length scale effect reflected in the effective stress-strain response, stress triaxiality and void fraction evolution. Furthermore, as expected, the conventionally used Gurson's model fails to capture the observed size-effects calling for a mechanistic modification that incorporates the mechanisms observed in our (and other researchers') simulation. Finally, in our multi-void simulations, we find that, the splitting of a big void into a distribution of small ones increases the load-carrying capacity of specimens. However, no obvious dependence of the void fraction evolution on void coalescence is observed.     

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

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

Scale-dependent plasticity potential of porous materials and void growth

In the present paper, a boundary value problem about the macroscopic response and its microscopic mechanism of a representative spherical cell with a spherical microvoid under axisymmetric triaxial tension has been theoretically investigated. To capture the size effects of local plastic deformation in the matrix, the strain gradient constitutive theory including the rotation and the stretch gradients developed by Fleck and Hutchinson [Strain gradient plasticity, in: J.W. Hutchinson, T.Y. Wu (Eds.), Advance in Applied Mechanics, vol. 33, Academic Press, New York, 1997, p. 295] is adopted. By means of the principle of minimum plasticity potential and the Lagrange multipliers method, the self-contained displacement field within the matrix has been computationally determined. Based on these, a size-dependent constitutive potential theory for porous material has been developed. The results indicate clearly that the microvoid evolution predicted by the present constitutive model displays very significant dependences on the void size especially when the radius a of microvoids is comparable with the intrinsic characteristic length l of the matrix. And when the void radius a is much lager than l, the present model can retrogress automatically to the Gurson model improved by Wang and Qin [Acta Mech. Solid. Sin. 10 (1989) 127, in Chinese].     

Lattice orientation effects on void growth and coalescence in fcc single crystals

Void growth and coalescence in fcc single crystals were studied using crystal plasticity under uniaxial and biaxial loading conditions and various orientations of the crystalline lattice. A 2D plane strain unit cell with one and two cylindrical voids was employed using three-dimensional 12 potentially active slip systems. The results were compared to five representative orientations of the tensile axis on the stereographic triangle. For uniaxial tension conditions, the void volume fraction increase under the applied load is strongly dependent on the crystallographic orientation with respect to the tensile axis. For some orientations of the tensile axis, such as [1 0 0] or [1 1 0], the voids exhibited a growth rate twice as fast compared with other orientations ([1 0 0], [2 1 1]). Void growth and coalescence simulations under uniaxial loading indicated that during deformation along some orientations with asymmetry of the slip systems, the voids experienced rotation and shape distortion, due mainly to lattice reorientation. Coalescence effects are shown to diminish the influence of lattice orientation on the void volume fraction increase, but noteworthy differences are still present. Under biaxial loading conditions, practically all differences in the void volume fraction for different orientations of the tensile axes during void growth vanish. These results lead to the conclusion that at microstructural length scales in regions under intense biaxiality/triaxiality conditions, such as crack tip or notched regions, the plastic anisotropy due to the initial lattice orientation has only a minor role in influencing the void growth rate. In such situations, void growth and coalescence are mainly determined by the stress triaxiality, the magnitude of accumulated strain, and the spatial localization of such plastic strains.     

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

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

A microscopic damage model considering the change of void shape and application in the void closing

A microscopic damage model of ellipsoidal body containing ellipsoidal void for nonlinear matrix materials is developed under a particular coordinate. The change of void shape is considered in this model. The viscous restrained equation obtained from the model is affected by stress ∑_(ij), void volume fraction f, material strain rate exponent m as well as the void shape. Gurson's equation is modified from the numerical solution. The modified equation is suitable for the case of nonlinear matrix materials and changeable voids. Lastly, the model is used to analyze the closing process of voids.     

On the micromechanics of void growth by prismatic-dislocation loop emission

Experimental evidence and recent molecular dynamics simulations of void growth indicate that prismatic dislocation loop emission by externally applied stresses is a viable mechanism of void growth under shock loading conditions when diffusive processes are given no time to operate. In this paper, the process of growth by loop emission is studied in a model system comprised of a void in an infinite linearly elastic and isotropic solid loaded axisymmetrically by remote applied stresses. First, the interaction between applied stresses, the stress field of a single dislocation loop or a pile-up of loops next to the void, the surface energy expenditure on void surface change, and the lattice resistance to the motion of loops is reviewed. The necessary condition for interstitial loop emission is used to determine the equilibrium positions of the loops as well as the maximum number of loops in a pile-up under given applied stresses. For the parameters of the model-material with purely hydrostatic loading, the numerical results yield a volume change for the void, which when normalized by the initial undeformed volume, exhibits a strong dependence on the size of the void for radii less than ∼400 times the lattice Burgers vector. For larger voids, the normalized volume change was found to be independent of the void radius.     

用三维电镜图像数据测试空穴形状

用碳化硅强化金属基复合材料制作的不同切口的拉伸试样进行了拉伸试验并使用三维电子扫描显微镜对拉伸延性断面进行了微小空穴的三维形状测试，分析了在不同应力三维度下空穴聚合时空穴几何形状的变化，为准确模拟金属基复合材料在多向外载荷下的损伤过程，判断延性损伤机理提供了科学依据。