Damage of ductile materials deforming under multiple plastic or viscoplastic mechanisms

This work presents a model to represent ductile failure (i.e. failure controlled by nucleation, growth and coalescence) of materials whose irreversible deformation is controlled by several plastic or viscoplastic deformation mechanisms. In addition work hardening may result from both isotropic and kinematic hardening. Damage is represented by a single variable representing void volume fraction. The model uses an additive decomposition of the plastic strain rate tensor. The model is developed based on the definition of damage dependant effective scalar stresses. The model is first developed within the generalized standard material framework and expressions for Helmholtz free energy, yield potential and dissipation potential are proposed. In absence of void nucleation, the evolution of the void volume fraction is governed by mass conservation and damage does not need to be represented by state variables. The model is extended to account for void nucleation. It is implemented in a finite element software to perform structural computations. The model is applied to three case studies: (i) failure by void growth and coalescence by internal necking (pipeline steel) where plastic flow is either governed by the Gurson–Tvergaard–Needleman model or the Thomason model, (ii) creep failure (Grade 91 creep resistant steel) where viscoplastic flow is controlled by dislocation creep or diffusional creep and (iii) ductile rupture after pre-compression (aluminum alloy) where kinematic hardening plays an important role.     

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.     

Predicting failure modes and ductility of dual phase steels using plastic strain localization

Ductile failure of metals is often treated as the result of void nucleation, growth and coalescence. Various criteria have been proposed to capture this failure mechanism for various materials. In this study, ductile failure of dual phase steels is predicted in the form of plastic strain localization resulting from the incompatible deformation between the harder martensite phase and the softer ferrite matrix. Microstructure-level inhomogeneity serves as the initial imperfection triggering the instability in the form of plastic strain localization during the deformation process. Failure modes and ultimate ductility of two dual phase steels are analyzed using finite element analyses based on the actual steel microstructures. The plastic work hardening properties for the constituent phases are determined by the in-situ synchrotron-based high-energy X-ray diffraction technique. Under different loading conditions, different failure modes and ultimate ductility are predicted in the form of plastic strain localization. It is found that the local failure mode and ultimate ductility of dual phase steels are closely related to the stress state in the material. Under plane stress condition with free lateral boundary, one dominant shear band develops and leads to final failure of the material. However, if the lateral boundary is constrained, splitting failure perpendicular to the loading direction is predicted with much reduced ductility. On the other hand, under plane strain loading condition, commonly observed necking phenomenon is predicted which leads to the final failure of the material. These predictions are in reasonably good agreement with experimental observations.     

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

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

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.     

A micromechanics based damage model for composite materials

The predictive capacity of ductile fracture models when applied to composite and multiphase materials is related to the accuracy of the estimated stress/strain level in the second phases or reinforcements, which defines the condition for damage nucleation. Second phase particles contribute to the overall hardening of the composite before void nucleation, as well as to its softening after their fracture or decohesion. If the volume fraction of reinforcement is larger than a couple of percents, this softening can significantly affect the resistance to plastic localization and cannot be neglected. In order to explicitly account for the effect of second phase particles on the ductile fracture process, this study integrates a damage model based on the Gologanu–Leblond–Devaux constitutive behavior with a mean-field homogenization scheme. Even though the model is more general, the present study focuses on elastic particles dispersed in an elasto-plastic matrix. After assessing the mean-field homogenization scheme through comparison with two-dimensional axisymmetric finite element calculations, an extensive parametric study is performed using the integrated homogenization-damage model. The predictions of the integrated homogenization-damage model are also compared with experimental results on cast aluminum alloys, in terms of both the fracture strain and overall stress–strain curves. The study demonstrates the complex couplings among the load transfer to second phase particles, their resistance to fracture, the void nucleation mode, and the overall ductility.     

Modelling of plastic flow localisation and damage development in friction stir welded 6005A aluminium alloy using physics based strain hardening law

Plastic flow localisation and ductile failure during tensile testing of friction stir welded aluminium specimens are investigated with a specific focus on modelling the local, finite strain, hardening response. In the experimental part, friction stir welds in a 6005A-T6 aluminium alloy were prepared and analysed using digital image correlation (DIC) during tensile testing as well as scanning electron microscopy (SEM) on polished samples and on fracture surfaces. The locations of the various regions of the weld were determined based on hardness measurements, while the flow behaviour of these zones was extracted from micro-tensile specimens cut parallel to the welding direction. The measured material properties and weld topology were introduced into a 3D finite element model, fully coupled with the damage model. A Voce law hardening model involving a constant stage IV is used within an enhanced Gurson type micro-mechanical damage model, accounting for void nucleation, growth and coalescence, as well as void shape evolution. The stage IV hardening, observed in Simar et al. (2010), was found to increase the stiffness during plastic flow localisation as well as to postpone the onset of fracture as determined by the void coalescence criterion. Furthermore, the presence of a second population of voids was concluded to strongly affect the fracture strain of the high strength regions of the welds. This modelling effort links the microstructure and process parameters to macroscopic parameters relevant to the optimisation of the welds.     

Modeling stress state dependent damage evolution in a cast Al–Si–Mg aluminum alloy

Internal state variable rate equations are cast in a continuum framework to model void nucleation, growth, and coalescence in a cast Al–Si–Mg aluminum alloy. The kinematics and constitutive relations for damage resulting from void nucleation, growth, and coalescence are discussed. Because damage evolution is intimately coupled with the stress state, internal state variable hardening rate equations are developed to distinguish between compression, tension, and torsion straining conditions. The scalar isotropic hardening equation and second rank tensorial kinematic hardening equation from the Bammann–Chiesa–Johnson (BCJ) Plasticity model are modified to account for hardening rate differences under tension, compression, and torsion. A method for determining the material constants for the plasticity and damage equations is presented. Parameter determination for the proposed phenomenological nucleation rate equation, motivated from fracture mechanics and microscale physical observations, involves counting nucleation sites as a function of strain from optical micrographs. Although different void growth models can be included, the McClintock void growth model is used in this study. A coalescence model is also introduced. The damage framework is then evaluated with respect to experimental tensile data of notched Al–Si–Mg cast aluminum alloy specimens. Finite element results employing the damage framework are shown to illustrate its usefulness.     

Effect of large elastic strains on cavitation instability predictions for elastic–plastic solids

For an infinite solid containing a void, the cavitation instability limit is defined as the remote stress–and strain state, at which the void grows without bound, driven by the elastic energy stored in the surrounding material. Such cavitation limits have been analysed by a number of authors for metal plasticity as well as for nonlinear elastic solids. The analyses for elastic–plastic solids are here extended to consider the effect of a large initial yield strain, and it is shown how the critical stress value decays for increasing value of the yield strain. Analyses are carried out for remote hydrostatic tension as well as for more general axisymmetric remote stress field, with an initially spherical void. Different levels of strain hardening are considered.     

An analysis of ductile rupture in notched bars

Ductile fracture in axisymmetric and plane strain notched tensile specimens is analyzed numerically, based on a set of elastic-plastic constitutive relations that account for the nucleation and growth of microvoids. Final material failure by void coalescence is incorporated into the constitutive model via the dependence of the yield function on the void volume fraction. In the analyses the material has no voids initially; but as the voids nucleate and grow, the resultant dilatancy and pressure sensitivity of the macroscopic plastic flow influence the solution significantly. Considering both a blunt notch geometry and a sharp notch geometry in the computations permits a study of the relative roles of high strain and high triaxiality on failure. Comparison is made with published experimental results for notched tensile specimens of high-strength steels. All axisymmetric specimens analyzed fail at the center of the notched section, whereas failure initiation at the surface is found in plane strain specimens with sharp notches, in agreement with the experiments. The results for different specimens are used to investigate the circumstances under which fracture initiation can be represented by a single failure locus in a plot of stress triaxiality vs effective plastic strain.     

Void growth and coalescence in ductile solids with stage III and stage IV strain hardening

State of the art ductile fracture models often rely on simple power laws to describe the strain hardening of the matrix material. Power laws do not distinguish between the two main stages of hardening observed in polycrystals, referred to as stage III and stage IV hardening, and which emerge from the evolution of the dislocation substructure. The aim of this study is to couple a physics based strain hardening law including these two stages to a micromechanics based ductile damage model. One of the main motivations is that, the stage IV constant hardening rate stage, occurring only at large strain, will be attained in most ductile failure problems if not at the overall level of deformation, at least locally around the growing voids. Furthermore, proper modelling of the stage III involving dislocation storage and recovery terms and the transition to stage IV provides a link with the underlying physical mechanisms of deformation and with the microstructure. First, in order to evaluate the effects of the stage III and stage IV hardening on void growth and coalescence, an extensive parametric study is performed on two-dimensional (2D) axisymmetric finite element (FE) unit cell calculations, using a Kocks-Mecking type hardening law. The cell calculations demonstrate that accounting for the stage IV hardening can have a profound effect on delaying void coalescence and increasing the ductility. The magnitude of the recovery term during stage III has also a significant effect on the void growth rate. Then, the Kocks-Mecking law is incorporated into the Gologanu-Leblond-Devaux (GLD) porous plasticity model supplemented by two different versions of the Thomason void coalescence criterion. The predictions of the damage model are in good agreement with the results of the FE calculations in terms of the stress-strain curves, the evolution of void shape and porosity, as well as the strain value at the onset of void coalescence.     

Determination of material intrinsic length and strain gradient hardening in microbending process

Microbending experiments of pure aluminum show that the springback angles increase with the decrease of foil thickness, which indicates obvious size effects and attributes to plastic strain gradient hardening. Then a constitutive model, taking into accounts both plastic strain and plastic strain gradient hardening, is proposed to analyze the microbending process of thin foil. The model is based on the relationship between shear yield stress and dislocation density, in which the material intrinsic length is related to material properties and average grain numbers along the characteristic scale direction of part. It is adopted in analytical model to calculate the non-dimensional bending moment and predict the springback angle after microbending. It is confirmed that the predictions by the proposed hardening model agree well with the experimental data, while those predicted by the classical plasticity model cannot capture such size effects. The contribution of plastic strain gradient increases with the decrease of foil thickness and is independent on the bending angle.     

The effect of porosity distribution on ductile failure

The effect of a nonuniform distribution of porosity on flow localization and failure in a porous material is analyzed numerically. The void density distribution and properties used to characterize the material behavior were obtained from measurements on partially consolidated and sintered iron powder. The calculations were carried out using an elastic viscoplastic constitutive relation for porous plastic solids. Local material failure is incorporated into the model through the dependence of the flow potential on void volume fraction. The region modelled is a small portion of a larger body, subject to various triaxial stress conditions. Both plane strain and axisymmetric deformations are considered with imposed periodic boundary conditions. Interactions between regions with higher void fractions promote plastic flow localization into a band. Local failure occurs by void growth and coalescence within the band. The results suggest a failure criterion based on a critical void volume fraction that is only weakly dependent on stress history. The critical void fraction does. however, depend on the initial void distribution and material hardening characteristics.     

一般加载规律的弹塑性本构关系

将有关文献给出一般加载规律一维全量理论的简单模型推广到一般加载规律的一维增量理论,进而推广到一般加载规律的多维增量理论,在此基础上,建立了推导一般加载规律的多维增量理论的本构关系的一种途径。应用这种途径,从应力空间的加载函数和应变空间的加载函数出发,推导了等向强化材料和被加热的等向强化材料的一般加载规律的弹塑性本构关系的两种表示形式。理论和实例均表明,这种途径对等向强化材料、随动强化材料和理想弹塑性材料均适用。     

Strain gradient crystal plasticity analysis of a single crystal containing a cylindrical void

The effects of void size and hardening in a hexagonal close-packed single crystal containing a cylindrical void loaded by a far-field equibiaxial tensile stress under plane strain conditions are studied. The crystal has three in-plane slip systems oriented at the angle 60° with respect to one another. Finite element simulations are performed using a strain gradient crystal plasticity formulation with an intrinsic length scale parameter in a non-local strain gradient constitutive framework. For a vanishing length scale parameter the non-local formulation reduces to a local crystal plasticity formulation. The stress and deformation fields obtained with a local non-hardening constitutive formulation are compared to those obtained from a local hardening formulation and to those from a non-local formulation. Compared to the case of the non-hardening local constitutive formulation, it is shown that a local theory with hardening has only minor effects on the deformation field around the void, whereas a significant difference is obtained with the non-local constitutive relation. Finally, it is shown that the applied stress state required to activate plastic deformation at the void is up to three times higher for smaller void sizes than for larger void sizes in the non-local material.     

Theoretical study of void closure in nonlinear plastic materials

Void closing from a spherical shape to a crack is investigated quantitatively in the present study. The constitutive relation of the Void-free matrix is assumed to obey the Norton power law. A representative volume element （RVE） which includes matrix and void is employed and a Rayleigh-Ritz procedure is developed to study the deformation-rates of a spherical void and a penny-shaped crack. Based on an approximate interpolation scheme, an analytical model for void closure in nonlinear plastic materials is established. It is found that the local plastic flows of the matrix material are the main mechanism of void deformation. It is also shown that the relative void volume during the deformation depends on the Norton exponent, on the far-field stress triaxiality, as well as on the far-field effective strain. The predictions of void closure using the present model are compared with the corresponding results in the literature, showing good agreement. The model for void closure provides a novel way for process design and optimization in terms of elimination of voids in billets because the model for void closure can easily be applied in the CAE analysis.     

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.     

Influence of void nucleation on ductile shear fracture at a free surface

An approximate continuum model of a ductile, porous material is used to study the influence of the nucleation and growth of micro-voids on the formation of shear bands and the occurrence of surface shear fracture in a solid subject to plane strain tension. Bifurcation into diffuse modes is analysed for a plane strain tensile specimen described by these constitutive relations, which account for a considerable plastic dilatancy due to void growth and for the possibility of non-normality of the plastic flow law. In particular, bifurcation into surface wave modes and the possible influence of such modes triggering shear bands is investigated. For solids with initial imperfactions such as a surface undulation, a local material inhomogeneity on an inclusion colony, the inception and growth of plastic flow localization is analysed numerically. Both the formation of void-sheets and the final growth of cracks in the shear bands is described numerically. Some special features of shear band development in the solid obeying non-normality are studied by a simple model problem.     

Post-yield material nonlinearity: Optimal homogeneous shear-frame sections and hysteretic behavior

A nonlinear hardening rule that describes a general yield surface change for frame members is proposed. The rule is valid for homogeneous materials and is defined by a nonlinear constitutive model. The model is given as the stress–strain relationship through the section thickness. The post-yield behavior of shear frame members is assessed using a hardening index parameter, a plastic strain coefficient, and a smooth and optimal mathematical function to model the geometry of structural wide-flange and tube sections. This enables a continuous distribution of post-yield curvatures to exist along the member length and results in accurate finite-element member displacement predictions. Steady-state hysteresis responses are determined for frame sections using the model’s nonlinearly degrading stiffness. The hardening rule is used to model the detailed material behavior of shear frame members that are subjected to lateral loads. The results of the model are validated through experimental published literature.