A modified Gurson model and its application to punch-out experiments

Recent experimental evidence has reiterated that ductile fracture is a strong function of stress triaxiality. Under high stress triaxiality loading, failure occurs as a result of void growth and subsequent necking of inter-void ligaments while under low stress triaxiality failure is driven by shear localization of plastic strain in these ligaments due to void rotation and distortion. The original Gurson model is unable to capture localization and fracture for low triaxiality, shear-dominated deformations unless void nucleation is invoked. A phenomenological modification to the Gurson model that incorporates damage accumulation under shearing has been proposed. Here we further extend the model and develop the corresponding numerical implementation method. Several benchmark tests are performed in order to verify the code. Finally, the model is utilized to model quasi-static punch-out experiments on DH36 steel. It is shown that the proposed modified Gurson model, in contrast to the original model, is able to capture the through-thickness development of cracks as well as the punch response. Thus, the computational fracture approaches based on the modified Gurson model may be applied to shear-dominated failures.     

On the numerical implementation of a shear modified GTN damage model and its application to small punch test

Gurson-type models have been widely used to predict failure during sheet metal forming process. However, a significant limitation of the original GTN model is that it is unable to capture fracture under relatively low stress triaxiality. This paper focused on the fracture prediction under this circumstance, which means shear-dominated stress state. Recently, a phenomenological modification to the Gurson model that incorporates damage accumulation under shearing has been proposed by Nahshon and Hutchinson. We further calibrated new parameters based on this model in 22MnB5 tensile process and developed the corresponding numerical implementation method. Lower stress triaxiality were realized by new-designed specimens. Subsequently, the related shear parameters were calibrated by means of reverse finite element method and the influences of new introduced parameters were also discussed. Finally, this shear modified model was utilized to model the small punch test (SPT) on 22MnB5 high strength steel. It is shown that the shear modification of GTN model is able to predict failure of sheet metal forming under wide range of stress state.     

Use of a modified Gurson model approach for the simulation of ductile fracture by growth and coalescence of microvoids under low, medium and high stress triaxiality loadings

In the paper the modified Gurson model is developed for the simulation of damage growth and ductile fracture under low, medium and high stress triaxiality loadings. A new coalescence criterion is introduced based on a simple assumption that singular value of the effective stress triggers the coalescence of microvoids in materials. According to the introduced approach the void coalescence described by means of the modified Gurson model is not only determined by the so-called critical, constant void volume fraction but also by the stress triaxiality ratio. Computational simulations have been carried out for Al 2024–T351 aluminum specimens. In order to find some improvements of micromechanical damage models, two different approaches have been compared for modeling the shear driven microvoid coalescence under low stress triaxiality loadings.     

3D modelling of plug failure in resistance spot welded shear-lab specimens (DP600-steel)

Ductile plug failure of resistance spot welded shear-lab specimens is studied by full 3D finite element analysis, using an elastic-viscoplastic constitutive relation that accounts for nucleation and growth of microvoids to coalescence (The Gurson model). Tensile properties and damage parameters are based on uni-axial tensile testing of the basis material, while the modelled tensile response of the shear-lab specimens is compared to experimental results for the case of a ductile failure near the heat affected zone (HAZ). A parametric study for a range of weld diameters is carried out, which makes it possible to numerically relate the weld diameter to the tensile shear force (TSF) and the associated displacement, u _TSF , respectively. Main focus in the paper is on modelling the localization of plastic flow and the corresponding damage development in the vicinity of the spot weld, near the HAZ. For decreasing weld diameter, localization of plastic flow may be observed to occur in the weld nugget, introducing significant shearing. Due to these competing mechanisms a critical transition radius of the weld may be found. However, due to the limitation of the Gurson model in describing ductile failure at very low stress triaxiality, further analysis of the shear failure is omitted.     

Use of a modified Gurson model for the failure behaviour of the clinched joint on Al6061 sheet

This paper discusses the ductile fracture behaviour of the clinched joint on Alloy 6061 sheets. Failure behaviour of the clinched joint is associated with the nucleation, growth and coalescence of voids within the microstructure. Various corrections to the original Gurson model are proposed to allow for instability and final fracture of the material. This paper is concerned with the application of the modified Gurson–Tvergaard–Needleman damage developed by Ken Nahshon and Zhenyu Xue. The results of the tension tests are compared with those of numerical analysis to obtain the initial void volume fraction f₀ and the shear damage coefficients k_w. The modified Gurson–Tvergaard–Needleman model is used to describe the failure behaviour and the shear strength of the clinched joint. Good agreement is shown between the experiment results and the numerical results on the location of fracture and the maximum failure load.     

Numerical modeling of elasto-plastic porous materials with void shape effects at finite deformations

A new constitutive model for elasto-plastic (rate-independent) porous materials subjected to general three-dimensional finite deformations is presented. The new model results from simple modifications of an earlier model of Kailasam and Ponte Castañeda (1997, 1998) [40], [41] so that it reproduces the exact spherical and cylindrical shell solution (composite sphere and composite cylinder assemblage) under purely hydrostatic loadings, while predicting (by calibration) accurately the void shape evolution according to the recent “second-order” model of Danas and Ponte Castañeda [17]. Furthermore, the present model is based on a rigorous homogenization method which is capable of predicting both the constitutive behavior and the microstructure evolution of porous materials. The microstructure is described by voids of arbitrary ellipsoidal shapes and orientations and as a result the material exhibits deformation-induced (or morphological) anisotropy at finite deformations. This is in contrast with the well-known Gurson [32] model which assumes that the voids remain spherical during the deformation process and thus the material remains always isotropic. The present model is implemented numerically in a finite element program where a three-dimensional thin-sheet (butterfly) specimen is subjected to a combination of shear and traction loading conditions in order to examine the effect of stress triaxiality and shearing upon material failure. The ability of the present model to take into account the nontrivial evolution of the microstructure and especially void shape effects leads to the prediction of material failure even at low stress triaxialities and small porosities without the use of additional phenomenological damage criteria. At high stress triaxialities, the present model gives similar predictions as the Gurson model.     

Influence of Triaxial Stress State on Ductile Fracture Strength of Polycrystalline Nickel

In this contribution, the effect associated with stress triaxiality on ductile damage evolution in high purity nickel has been investigated from both experimental and theoretical points of view. Tensile tests on smooth and notched round bar specimens were performed to calibrate the fracture strain in a wide range of stress triaxiality. The capability of the Gurson model to reproduce and predict physical failure behaviour was examined. It was shown that stress triaxiality played a major role on damage evolution as demonstrated by the progressive reduction of material ductility under increasing triaxial states of stress.     

Continuum damage mechanics modelling incorporating stress triaxiality effect on ductile damage initiation

Micromechanical modelling of void nucleation in ductile metals indicates that strain required for damage initiation reduces exponentially with increasing stress triaxiality. This feature has been incorporated in a continuum damage mechanics (CDM) model, providing a phenomenological relationship for the damage threshold strain dependence on the stress triaxiality. The main consequences of this model modification are that the failure locus is predicted to change as function of stress triaxiality sensitivity of the material damage threshold strain and that high triaxial fracture strain is expected to be even lower than the threshold strain at which the damage processes initiate at triaxiality as low as 1/3. The proposed damage model formulation has been used to predict ductile fracture in unnotched and notched bars in tension for two commercially pure α‐iron grades (Swedish and ARMCO iron). Finally, the model has been validated, predicting spall fracture in a plate‐impact experiment and confirming the capability to capture the effect of the stress state on material fracture ductility at very high stress triaxiality.     

Fracture Criteria for Automobile Crashworthiness Simulation of Wrought Aluminium Alloy Components

In automobile crashworthiness simulation, the prediction of plastic deformation and fracture of each significant, single component is critical to correctly represent the transient energy absorption through the car structure. There is currently a need, in the commercial FEM community, for validated material fracture models which adequately represent this phenomenon. The aim of this paper is to compare and to validate existing numerical approaches to predict failure with test data. All studies presented in this paper were carried out on aluminium wrought alloys: AlMgSi1.F31 and AlMgSiCu‐T6. A viscoplastic material law, whose parameters are derived from uniaxial tensile and compression tests at various strain rates, is developed and presented herein. Fundamental ductile fracture mechanisms such as void nucleation, void growth, and void coalescence as well as shear band fracture are present in the tested samples and taken into consideration in the development of the fracture model. Two approaches to the prediction of fracture initiation are compared. The first is based on failure curves expressed by instantaneous macroscopic stresses and strains (i. e. maximum equivalent plastic strain vs. stress triaxiality). The second approach is based on the modified Gurson model and uses state variables at the mesoscopic scale (i. e. critical void volume fraction). Notched tensile specimens with varying notch radii and axisymmetric shear specimens were used to produce ductile fractures and shear band fractures at different stress states. The critical macroscopic and mesoscopic damage values at the fracture initiation locations were evaluated using FEM simulations of the different specimens. The derived fracture criteria (macroscopic and mesoscopic) were applied to crashworthiness experiments with real components. The quality of the prediction on component level is discussed for both types of criteria.     

Stress triaxiality effect on void nucleation in ductile metals

The stress triaxiality effect on the strain required for void nucleation by particle‐matrix debonding has been investigated by means of micromechanical modelling. A unit‐cell model considering an elastic spherical particle embedded in an elastic‐plastic matrix was developed to the purpose. Particle‐matrix decohesion was simulated through the progressive failure of a cohesive interface. It has been shown that the parameters of matrix‐particle cohesive interface are correlated with macroscopic material properties. Here, a simple relationship for the maximum cohesive opening at interface failure as a function of material fracture toughness and yield stress has been derived. Results seem to confirm that, increasing stress triaxiality, the strain at which void nucleation is predicted to occur decreases exponentially in a similar way as for fracture strain. This result has substantial implications in modelling of ductile damage because it indicates that if the stress triaxiality is high enough, ductile fracture can occur at plastic strain lower than that necessary to nucleate damage for moderate or low stress triaxiality regime.     

APPLICATION OF LOCAL DAMAGE MODELS TO THE NUMERICAL ANALYSIS OF DUCTILE RUPTURE

Abstract— Two damage models were implemented into the finite element program ADINA to study the correlation between microscopical damage and macroscopical material failure. In the first model, based on the Gurson yield function the nucleation, growth and the coalescence of voids were incorporated into the constitutive relations. In the second model the void growth was determined according to the Rice and Tracey model using the von Mises yield function, and material failure was simulated by eliminating the elements where the critical void growth ratio was exceeded. The numerical results for the local and global behaviour of the specimens were compared with experiments. The generality of the damage parameters was checked by investigating several specimen geometries. Both damage models deliver qualitatively consistent results with regard to the influence of the stress triaxiality on the void growth and on the beginning of the material failure. However, the Gurson model gives a more accurate numerical simulation because the damage development and the stress drop continue after the onset of void coalescence while the critical void growth model causes less convergence problems in the simulation of large crack extension. The J _n-curve was estimated on the basis of both models.     

Collapse and coalescence of spherical voids subject to intense shearing: studied in full 3D

Micro-mechanical 2D cell model studies have revealed ductile failure during intense shearing to be governed by the interaction of neighbouring voids, which collapse to micro-cracks and continuously rotate and elongate until coalescence occurs. For a three-dimensional void structure, this implies significant straining of the matrix material located on the axis of rotation. In particular, the void surface material is severely deformed during shearing and void surface contact is established early in the deformation process. This 3D effect intensifies with decreasing stress triaxiality and complicates the numerical analysis, which is also reflected in published literature. Rather than moving towards very low triaxiality shearing, work has focused on extracting wide-ranging results for moderate stress triaxiality (T ~ 1), in order to achieve sufficient understanding of the influence of initial porosity, void shape, void orientation etc. The objective of this work is to expand the range of stress triaxiality usually faced in 3D cell model studies, such that intense shearing is covered, and to bring forward details on the porosity and void shape evolution. The overall material response is presented for a range of initial material configurations and loading conditions. In addition, a direct comparison to corresponding 2D cell model predictions for circular cylindrical voids under plane strain shearing is presented. A quantitatively good agreement of the two model configurations (2D vs. 3D) is obtained and similar trends are predicted. However, the additional layer of matrix material, connecting voids in the transverse direction, is concluded to significantly influence the void shape evolution and to give rise to higher overall ductility. This 3D effect is demonstrated for various periodic distributions of voids.     

Effects of the stress state on plastic deformation and ductile failure: Experiment and numerical simulation using a newly designed tension‐shear specimen

The stress state is one of the most notable factors that dominates the initiation of ductile fracture. To examine the effects of the stress state on plasticity and ductile failure, a new tension‐shear specimen that can cover a wide range of stress triaxialities was designed. A fracture locus was constructed in the space of ductility and stress triaxiality for two typical steels based on a series of tests. It is observed that the equivalent plastic strain at failure exhibits a nonmonotonic variation with increasing the value of stress triaxiality. A simple damage model based on the ductility exhaustion concept was used to simulate the failure behaviour, and a good agreement is achieved between simulation results and experimental data. It is further shown that consideration of fracture locus covering a wide range of stress triaxialities is a key to an accurate prediction.     

Calibration procedures for a computational model of ductile fracture

A recent extension of the Gurson constitutive model of damage and failure of ductile structural alloys accounts for localization and crack formation under shearing as well as tension. When properly calibrated against a basic set of experiments, this model has the potential to predict the emergence and propagation of cracks over a wide range of stress states. This paper addresses procedures for calibrating the damage parameters of the extended constitutive model. The procedures are demonstrated for DH36 steel using data from three tests: (i) tension of a round bar, (ii) mode I cracking in a compact tension specimen, and (iii) shear localization and mode II cracking in a shear-off specimen. The computational model is then used to study the emergence of the cup-cone fracture mode in the neck of a round tensile bar. Ductility of a notched round bar provides additional validation.     

Ductile failure behavior of polycrystalline Al 6061-T6 under shear dominant loading

Ductile failure in polycrystalline aluminum alloys under pure shear as well as with superposed tension and compression loading is explored through the modified Arcan shear experiments. Specimens obtained through tests interrupted at various stages of deformation and failure evolution are examined through quantitative microscopy to discern the mechanisms of failure and to evaluate the local strain evolution quantitatively. Fractographic observations are used to identify the onset and evolution of damage processes during deformation and failure of these aluminum alloys. Local strain levels are estimated from measurements of the change in grain size with deformation and used to indicate that the local values of failure strains are likely to be much larger than that estimated from strains averaged over characteristic specimen dimensions such as the gage length or the specimen diameter. Lower bound estimates of the failure strain in low triaxiality conditions are obtained from the experiments. It is shown that strain-to-failure decreases monotonically with stress triaxiality in stark contrast with recent works where a reverse behavior in low stress triaxiality levels was reported. Eventual failure that occurs through void growth and coalescence is shown to be restricted to a very small region within the localized deformation band.     

Dynamic vs. quasi-static shear failure of high strength metallic alloys: Experimental issues

Ductile fracture of metals by void nucleation, growth and coalescence under positive stress triaxiality is well admitted. This is not the case when metals are submitted to negative stress triaxiality. The present work aims at contributing to a better understanding of the competition between micro-mechanisms at the origin of failure of metals when submitted to shear-pressure loading at low and high strain rates. With this aim in view, experiments were carried out on Ti–6Al–4V shear-compression samples involving a stress triaxiality range comprised between −0.2 and −0.5. Results show that the material failure is the consequence of a void growth induced process. At high strain rate, due to the localization of the deformation within adiabatic shear bands, the failure of the material occurs earlier, leading to maximum shear strain smaller at high strain rate than at low strain rate. Impact tests were also carried out on Kalthoff and Winkler type double notched plates. They showed that the interaction between tension and shear waves leads to a complex Mode I–Mode II crack propagation.     

基于微孔贯通细观损伤模型的金属韧性断裂分析

韧性材料断裂过程通常可看作是材料内部微孔洞的形核、扩展及相互贯通的积累。经典的Gurson- Tvergaard (GT)模型能够很好地模拟具有变形均匀、各向同性的孔洞的萌生及扩展过程;但无法模拟由孔洞贯通而引起的局部变形过程,因此需要对其修正,引入相应的孔洞贯通准则。该文采用两种贯通准则对经典GT模型进行修正,即Thomason的塑性极限载荷准则和临界等效塑性应变准则。借助用户自定义程序UMAT将采用这两种贯通准则修正的GT本构关系嵌入至商用有限元软件ABAQUS中,从而可通过对金属材料应力状态和断裂机理的分析控制孔洞的贯通。以一组含有不同缺口根半径的圆棒拉伸试验件为例,分析了该类金属构件自孔洞萌生至最终断裂的整个损伤演化过程,并与试验数据进行了对比,验证了该模型的有效性。该文还讨论了金属断裂过程中应力三轴度对微裂纹萌生与扩展的影响。     

Modeling of failure mode of laser welds in lap-shear specimens of HSLA steel sheets

Failure mode of laser welds in lap-shear specimens of high strength low alloy (HSLA) steel sheets is investigated in this paper. The experiments for laser welds in lap-shear specimens under quasi-static loading conditions are briefly reviewed first. The experimental results showed that the laser welds failed in a ductile necking/shear failure mode and the ductile failure was initiated at a distance away from the crack tip near the boundary of the base metal and heat affected zone. In order to understand the failure mode of these welds, finite element analyses under plane strain conditions were conducted to identify the effects of the different plastic behaviors of the base metal, heat affected zone, and weld zone as well as the weld geometry on the ductile failure. The results of the reference finite element analysis based on the homogenous material model show that the failure mode is most likely to be a middle surface shear failure mode in the weld. The results of the finite element analysis based on the multi-zone non-homogeneous material models show that the higher effective stress–plastic strain curves of the weld and heat affected zones and the geometry of the weld protrusion result in the necking/shear failure mode in the load carrying sheet. The results of another finite element analysis based on the non-homogeneous material model and the Gurson yield function for porous materials indicate that the consideration of void nucleation and growth is necessary to identify the ductile failure initiation site that matches well with the experimental observations. Finally, the results of this investigation indicate that the failure mode of the welds should be examined carefully and the necking/shear failure mode needs to be considered for development of failure or separation criteria for welds under more complex loading conditions.     

Three-dimensional cell model analyses of void growth in ductile materials

Three-dimensional micromechanical models were developed to study the damage by void growth in ductile materials. Special emphasis is given to the influence of the spatial arrangement of the voids. Therefore, periodical void arrays of cubic primitive, body centered cubic and hexagonal structure are investigated by analyzing representative unit cells. The isotropic behaviour of the matrix material is modelled using either v. Mises plasticity or the modified Gurson-Tvergaard constitutive law. The cell models are analyzed by the large strain finite element method under monotonic loading while keeping the stress triaxiality constant. The obtained mesoscopic deformation response and the void growth of the unit cells show a high dependence on the value of triaxiality. The spatial arrangement has only a weak influence on the deformation behaviour, whereas the type and onset of the plastic collapse behaviour are strongly affected. The parameters of the Gurson-Tvergaard model can be calibrated to the cell model results even for large porosity, emphasizing its usefulness and justifying its broad applicability.