Modification of the Gurson Model for shear failure

Recent experimental evidence points to limitations in characterizing the critical strain in ductile fracture solely on the basis of stress triaxiality. A second measure of stress state, such as the Lode parameter, is required to discriminate between axisymmetric and shear-dominated stress states. This is brought into the sharpest relief by the fact that many structural metals have a fracture strain in shear, at zero stress triaxiality, that can be well below fracture strains under axisymmetric stressing at significantly higher triaxiality. Moreover, recent theoretical studies of void growth reveal that triaxiality alone is insufficient to characterize important growth and coalescence features. As currently formulated, the Gurson Model of metal plasticity predicts no damage change with strain under zero mean stress, except when voids are nucleated. Consequently, the model excludes shear softening due to void distortion and inter-void linking. As it stands, the model effectively excludes the possibility of shear localization and fracture under conditions of low triaxiality if void nucleation is not invoked. In this paper, an extension of the Gurson model is proposed that incorporates damage growth under low triaxiality straining for shear-dominated states. The extension retains the isotropy of the original Gurson Model by making use of the third invariant of stress to distinguish shear dominated states. The importance of the extension is illustrated by a study of shear localization over the complete range of applied stress states, clarifying recently reported experimental trends. The extension opens the possibility for computational fracture approaches based on the Gurson Model to be extended to shear-dominated failures such as projectile penetration and shear-off phenomena under impulsive loadings.     

On the mechanism of void growth and the effect of straining mode in ductile materials

Finite element analysis is employed to investigate void growth embedded in elastic–plastic matrix material. Axisymmetric and plane stress conditions are considered. The simulation of void growth in a unit cell model is carried out over a wide range of triaxial tensile stressing or large plastic straining for various strain hardening materials to study the mechanism of void growth in ductile materials. Triaxial tension and large plastic strain encircling around the void are found to be of most importance for driving void growth. The straining mode of incremental loading which favors the necessary strain concentration around void for its growth can be characterized by the vanishing condition of a parameter called “the third invariant of generalized strain rate”. Under this condition, it accentuates the internal strain concentration and the strain energy stored/dissipated within the material layer surrounding the void. Experimental results are cited to justify the effect of this loading parameter.     

Coupling effects of void size and void shape on the growth of prolate ellipsoidal microvoid

The combined effects of void size and void shape on the void growth are studied by using the classical spectrum method. An infinite solid containing an isolated prolate spheroidal void is considered to depict the void shape effect and the Fleck-Hutchinson phenomenological strain gradient plasticity theory is employed to capture the size effects. It is found that the combined effects of void size and void shape are mainly controlled by the remote stress triaxiality. Based on this, a new size-dependent void growth model similar to the Rice-Tracey model is proposed and an important conclusion about the size-dependent void growth is drawn: the growth rate of the void with radius smaller than a critical radius r_c may be ignored. It is interesting that r_c is a material constant independent of the initial void shape and the remote stress triaxiality.The project supported by the National Natural Science Foundation of China (A10102006) and the New Century Excellent Talents in Universities of China. The English text was polished by Keren Wang.     

The evolution of voids in the plasticity strain hardening gradient materials

The main aim of this paper is to opens out the meso-mechanism of void growth and coalescence in the matrix materials with graded strain-hardening exponent distribution. For this end, detailed finite element computations of a representative cylindrical cell containing a spherical void have been carried out. According to the FE analyses, significant effects of the strain-hardening exponent gradient (SEG) in the matrix on the void growth and coalescence are revealed: (1) In the homogeneous materials, the void growth and coalescence are slightly dependent on the strain-hardening exponent, however, the SEG distribution in the matrix can increase remarkably the void growth rate and decrease seriously the void coalescence strain. (2) The critical void shapes in the homogeneous materials are mainly governed by the macroscopic stress triaxiality, but due to earlier plastic flow localization in the softer matrix layer, the SEG distribution in the matrix has very significant effects on the deformed void shapes, especially when the stress triaxiality is lower. (3) When the triaxial stress levels are lower, in the homogeneous materials, the shape change mode of the void evolution is dominate so the void growth rate is very low; however, the SEG distribution in the matrix can bring the volume change mode out, as a result of increasing the void growth rate. (4) Comparisons of the numerical results with the existing damage model indicate that the classic damage model cannot give satisfying prediction to the void growth in both the homogeneous strain-hardening matrix and the SEG materials. On the basis of large numbers of numerical computations, a new damage model, which can uniformly describe the void growing in the homogeneous and plasticity gradient materials, is suggested. A mass of element computations have validated that the new damage model can give satisfying agreement with the FE results of cell model.     

SURFACE EFFECT ON NANOSIZED VOID GROWTH IN A RIGID-PERFECTLY PLASTIC MATERIAL

The influence of the surface effect on the nanosized spherical void growth in a rigidperfectly plastic material is analyzed and the mechanism of the nanosized void growth with high triaxiality is given. Based on the Rice and Tracey model for a macro void growth, the present model is proposed to account for the nanosized void growth under a uniform remote strain rate field with consideration on the surface effect. It is concluded that the surface effect yields an evident resistant influence on the nanosized void growth. That is, this influence decays as the void radius increases. With high triaxiality, the nanosized void growth is divided into two stages： the initial stage and the mature stage. At the first stage, the void grows slowly and the influence of surface effect is relatively weak, whereas at the second stage, the influence is significant and the void grows drastically.     

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

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

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.     

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

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

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.     

三轴应力场中不同形状孔洞的长大及其新模型

对不同形状孔洞在从光滑试样到裂纹试样这样广泛三轴应力场中的长大规律,本文通过控制体胞宏观应力三维度的方法进行了精确的有限元分析,计算结果表明：（1）孔洞的体积改变和形状变化是孔洞演化的两种基本机制,在不同的三轴应力场中,这两种机制的作用不同;（2）现有模型对孔洞长大规律的描述是不准确的,由它们得到的临界孔洞扩张比参数HGC与临界孔洞体积分数fc不具备一一对应关系,因此不以很好地反映也洞的实际扩张。在此基础上,提出了一个描述孔洞长大的新模型,与四种常用的现有模型相比,该模型不仅能更好地描述不同三轴应力场中孔洞的长大,而且能反映不同应力三维度水平下材料破坏模式的变化。     

Anisotropic finite-strain models for porous viscoplastic materials with microstructure evolution

In this work, we propose a constitutive model for the finite-strain, macroscopic response of porous viscoplastic solids, accounting for deformation-induced changes in the size, shape and distribution of the voids. The model makes use of consistent homogenization estimates obtained by the “iterated variational linear comparison” procedure of Agoras and Ponte Castañeda (2013) to characterize both the instantaneous effective response of the porous material and the evolution of the underlying microstructure. The proposed model applies for general, three-dimensional loading conditions and can be implemented numerically for use in standard FEM codes. We also investigate the interplay between the evolution of the microstructure and the macroscopic stress–strain response, in the context of displacement-controlled, plane strain loading (bi-axial straining) of initially isotropic, porous, rigid-plastic materials with power-law hardening. We focus on the effect of strain triaxiality, and consider both extensional and contractile loading conditions leading to porosity growth and collapse, respectively. For both types of loadings, it is found that the macroscopic behavior of the material exhibits an initial hardening regime followed by a softening regime at sufficiently large strains. Consistent with earlier models and experimental results, the softening regime for extensional loadings is a consequence of porosity growth. On the other hand, the softening behavior predicted for contractile loadings is found to be a consequence of void collapse, i.e., of rapid changes in the average shape of the pores leading to crack-like shapes. For both types of loadings, the transition from hardening to softening in the macroscopic response can be identified with the onset of macroscopic strain localization. The associated critical conditions at the onset of localization are determined as a function of the strain triaxiality. The type of localization band ranges from dilatational to compaction bands as the bi-axial straining varies from uniaxial extension to uniaxial contraction.     

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.     

Influences of initial porosity,stress triaxiality and Lode parameter on plastic deformation and ductile fracture

Local mechanical properties in aluminum cast components are inhomogeneous as a consequence of spatial distribution of microstructure,e.g.,porosity,inclusions,grain size and arm spacing of secondary dendrites.In this work,the effect of porosity is investigated.Cast components contain voids with different sizes,forms,orientations and distributions.This is approximated by a porosity distribution in the following.The aim of this paper is to investigate the influence of initial porosity,stress triaxiality and Lode parameter on plastic deformation and ductile fracture.A micromechanical model with a spherical void located at the center of the matrix material,called the representative volume element(RVE),is developed.Fully periodic boundary conditions are applied to the RVE and the values of stress triaxiality and Lode parameter are kept constant during the entire course of loading.For this purpose,a multi-point constraint(MPC)user subroutine is developed to prescribe the loading.The results of the RVE model are used to establish the constitutive equations and to further investigate the influences of initial porosity,stress triaxiality and Lode parameter on elastic constant,plastic deformation and ductile fracture of an aluminum die casting alloy.     

Experimental investigation and modelling of compressibility induced by damage in carbon black-reinforced natural rubber

While natural rubber is commonly considered as an incompressible material, this study shows how carbon black-reinforced natural rubber (NR-CB), when subjected to various mechanical loading conditions (uniaxial, hydrostatic, monotonic, cyclic), is affected by volume change. Experiments show a volume variation even for low straining values and a significant volume change for large elongations. Moreover, volume change can be either reversible or not, depending on the loading conditions. It is related to a competition between void growth, chain orientation, and stress softening. At a microscopic scale, in situ Scanning Electron Microscopy (SEM) examinations and image analysis allow one to record damage and microscopic volume change as a function of elongation. Therefore the volume change measured at the microscopic scale is equal to the macroscopic one. Based on the experimental results, this paper shows that the hypothesis of incompressibility is worth being revisited. Thus, a nearly compressible approach was considered, where the strain energy is assumed to be the sum of spherical and deviatoric parts that are both affected by damage. The model was then implemented in a finite-element code. Good agreement was obtained between experimental results and model predictions for low triaxiality test conditions.     

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.     

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.     

Rupture mechanisms in combined tension and shear—Micromechanics

A micromechanics model based on the theoretical framework of plastic localization into a band introduced by Rice is developed. The model consists of a planar band with a square array of equally sized cells, with a spherical void located in the centre of each cell. The periodic arrangement of the cells allows the study of a single unit cell for which fully periodic boundary conditions are applied. The micromechanics model is applied to analyze failure by ductile rupture in experiments on double notched tube specimens subjected to combined tension and torsion carried out by the present authors. The stress state is characterized in terms of the stress triaxiality and the Lode parameter. Two rupture mechanisms can be identified, void coalescence by internal necking at high triaxiality and void coalescence by internal shearing at low triaxiality. For the internal necking mechanism, failure is assumed to occur when the deformation localizes into a planar band and is closely associated with extensive void growth until impingement of voids. For the internal shearing mechanism, a simple criterion based on the attainment of a critical value of shear deformation is utilized. The two failure criteria capture the transition between the two rupture mechanisms successfully and are in good agreement with the experimental result.     

Evaluation of Ductile Damage Parameters under High Stress Triaxiality

The GTN model proposed by Gurson, Tvergaard and Needleman has been widely applied to predict ductile fracture. However, the evaluation of the GTN model under high stress triaxiality has only been reported in a few studies. In this paper, a series of tensile tests on round notched specimens were performed to evaluate the applicability of the GTN model parameters under high stress triaxiality. The evaluation was carried out by comparing the predicted load-displacement curves with experimental results. It was observed the GTN model parameters only depend on the material except the critical void volume fraction. The influence of stress triaxiality on the critical void volume fraction was discussed. A further discussion about the construction of a new void coalescence criterion for the GTN model was also presented in this paper.     

On the extension of the Gurson-type porous plasticity models for prediction of ductile fracture under shear-dominated conditions

One of the major drawbacks of the Gurson-type of porous plasticity models is the inability of these models to predict material failure under low stress triaxiality, shear dominated conditions. This study addresses this issue by combining the damage mechanics concept with the porous plasticity model that accounts for void nucleation, growth and coalescence. In particular, the widely adopted Gurson–Tvergaard–Needleman (GTN) model is extended by coupling two damage parameters, representing the volumetric damage (void volume fraction) and the shear damage, respectively, into the yield function and flow potential. The effectiveness of the new model is illustrated through a series of numerical tests comparing its performance with existing models. The current model not only is capable of predicting damage and fracture under low (even negative) triaxiality conditions but also suppresses spurious damage that has been shown to develop in earlier modifications of the GTN model for moderate to high triaxiality regimes. Finally the modified GTN model is applied to predict the ductile fracture behavior of a beta-treated Zircaloy-4 by coupling the proposed damage modeling framework with a recently developed J₂–J₃ plasticity model for the matrix material. Model parameters are calibrated using experimental data, and the calibrated model predicts failure initiation and propagation in various specimens experiencing a wide range of triaxiality and Lode parameter combinations.