On strain and damage interactions during tearing: 3D in situ measurements and simulations for a ductile alloy (AA2139-T3)

Strain and damage interactions during tearing of a ductile Al-alloy with high work hardening are assessed in situ and in 3D combining two recently developed experimental techniques, namely, synchrotron laminography and digital volume correlation. Digital volume correlation consists of registering 3D laminography images. Via simultaneous assessments of 3D strain and damage at a distance of 1-mm ahead of a notch root of a thin Compact Tension-like specimen, it is found that parallel crossing slant strained bands are active from the beginning of loading in a region where the crack will be slanted. These bands have an intermittent activity but are stable in space. Even at late stages of deformation strained bands can stop their activity highlighting the importance of plasticity on the failure process rather than damage softening. One void is followed over the loading history and seen to grow and orient along the slant strained band at very late stages of deformation. Void growth and strain are quantified. Gurson–Tvergaard–Needleman-type simulations using damage nucleation for shear, which is based on the Lode parameter, are performed and capture slant fracture but not the initial strain fields and in particular the experimentally found slant bands. The band formation and strain distribution inside and outside the bands are discussed further using plane strain simulations accounting for plastic material heterogeneity in soft zones.     

Vapor pressure and void size effects on failure of a constrained ductile film

To achieve certain properties, semiconductor adhesives and molding compounds are made by blending filler particles with polymer matrix. Moisture collects at filler particle/polymer matrix interfaces and within voids of the composite. At reflow temperatures, the moisture vaporizes. The rapidly expanding vapor creates high internal pressure on pre-existing voids and particle/matrix interfaces. The simultaneous action of thermal stresses and internal vapor pressure drives both pre-existing and newly nucleated voids to grow and coalesce causing material failure. Particularly susceptible are polymeric films and adhesives joining elastic substrates, e.g. Ag filled epoxy. Several competing failure mechanisms are studied including: near-tip void growth and coalescence with the crack; extensive void growth and formation of an extended damaged zone emanating from the crack; and rapid void growth at highly stressed sites at large distances ahead of the crack, leading to multiple damaged zones. This competition is driven by the interplay between stress elevation induced by constrained plastic flow and stress relaxation due to vapor pressure assisted void growth.A model problem of a ductile film bonded between two elastic substrates, with a centerline crack, is studied. The computational study employs a Gurson porous material model incorporating vapor pressure effects. The formation of multiple damaged zones is favored when the film contains small voids or dilute second-phase particle distribution. The presence of large voids or high vapor pressure favor the growth of a self-similar damage zone emanating from the crack. High vapor pressure accelerates film cracking that can cause device failures.     

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.     

A constitutive model for elastoplastic solids containing primary and secondary voids

In many ductile metallic alloys, the damage process controlled by the growth and coalescence of primary voids nucleated on particles with a size varying typically between 1 and 100 μm, is affected by the growth of much smaller secondary voids nucleated on inclusions with a size varying typically between 0.1 and 3 μm. The goal of this work is first to quantify the potential effect of the growth of these secondary voids on the coalescence of primary voids using finite element (FE) unit cell calculations and second to formulate a new constitutive model incorporating this effect. The nucleation and growth of secondary voids do essentially not affect the growth of the primary voids but mainly accelerate the void coalescence process. The drop of the ductility caused by the presence of secondary voids increases if the nucleation strain decreases and/or if their volume fraction increases and/or if the primary voids are flat. A strong coupling is indeed observed between the shape of the primary voids and the growth of the second population enhancing the anisotropy of the ductility induced by void shape effects. The new micromechanics-based coalescence condition for internal necking introduces the softening induced by secondary voids growing in the ligament between two primary voids. The FE cell calculations were used to guide and assess the development of this model. The use of the coalescence condition relies on a closed-form model for estimating the evolution of the secondary voids in the vicinity of a primary cavity. This coalescence criterion is connected to an extended Gurson model for the first population including the effect of the void aspect ratio. With respect to classical models for single void population, this new constitutive model improves the predictive potential of damage constitutive models devoted to ductile metal while requiring only two new parameters, i.e. the initial porosity of second population and a void nucleation stress, without any additional adjustment.     

Prediction of ductile fracture in tension by bifurcation,localization, and imperfection analyses

In this investigation, it is shown that the onset of ductile fracture in tension can be interpreted as the result of a supercritical bifurcation of homogeneous deformation and that this fact can be applied to predict ductile fracture initiation of materials with general imperfections or flaws. We focus on one dimensional quasi-static simple tension for rate-independent isotropic plastic materials. For deformation beyond the bifurcation point, multiple equilibrium paths appear. The homogeneous deformation, as one of the equilibrium paths, loses stability while the inhomogeneous paths are stable, thus indicating the occurrence of strain localization. This investigation also provides a physical example for the application of the Lambert W function in material localization analyses. Material instability is treated as the instability of a static system with dynamic perturbation. We also address the presence of microvoids in a power law plastic material as an unfolding of the supercritical pitchfork bifurcation. The imperfect system, idealized as spherical voids within the plastic matrix, is analyzed using the familiar Gurson model which is based on the presumption of a randomly voided material and characterized by the volume fraction of voids. If, in addition, the sizes of the microvoids are known, this then provides a length scale for the imperfection zone. In this manner, relevance to the sample size effects of strain-to-failure for ductile fracture initiation is addressed by considering separate zones with variations in void volume fractions. Fracture initiation predictions are presented and compare very well to existing experimental results.     

Microstructure and fracture characteristic of Mg–Al–Zn–Si3N4 composites

Magnesium matrix composites reinforced by three dimensional (3-D) network structure were fabricated by pressure-less infiltration technology. The 3-D network structure reinforcement and its composites exhibited special topology structure and different fracture characteristic. Metal matrix fractured in a ductile mode is manifested by small dimples and craters on the fracture surface. When the volume fraction of reinforcement is not in excess of 6%, the composite had an improved fracture toughness. This is because of the relatively homogeneous Si₃N₄ particles distributed in the metal matrix and the occurrence of interface reaction product such as MgAlO₂ spinel phase. With the increases of volume fraction of reinforcement (>6%), the fracture toughness decreases slowly at the initial stages and then decreases rapidly towards the end. Therefore, the main fracture failure mechanisms consist of crack nucleation, growth, coalescence and crack propagation.     

An interactive micro-void shear localization mechanism in high strength steels

Ductility of high strength steels is often restricted by the onset of a void-sheet mechanism in which failure occurs by a micro-void shear localization process. For the first time, the micro-void shear instability mechanism is identified here by examining the interactions occurring within a system of multiple embedded secondary particles (carbides ∼10-100 nm), through a finite element based computational cell modeling technique (in two and three dimensions). Shear deformation leads to the nucleation of micro-voids as the secondary particles debond from the surrounding alloy matrix. The nucleated micro-voids grow into elongated void tails along the principal shear plane and coalesce with the micro-voids nucleated at neighboring particles. At higher strains, the neighboring particles are driven towards each other, further escalating the severity of the shear coalescence effect. This shear driven nucleation, growth and coalescence mechanism leads to a decrease in the load-bearing surface in the shear plane and a terminal shear instability occurs. The mechanism is incorporated mathematically into a hierarchical steel model. The simulated response corresponds to experimentally observed behavior only when the micro-void shear localization mechanism is considered.     

A multiresolution continuum simulation of the ductile fracture process

With the advancement in computational science that is stepping into the Exascale era and experimental techniques that enable rapid reconstruction of the 3D microstructure, quantitative microstructure simulations at an unprecedented fidelity level are giving rise to new possibilities for linking microstructure to property. This paper presents recent advances in 3D computational modeling of ductile fracture in high toughness steels. Ductile fracture involves several concurrent and mutually interactive mechanisms at multiple length scales of microstructure. With serial sectioning tomographic techniques, a digital dataset of microstructure features associated with the fracture process has been experimentally reconstructed. In this study, primary particles are accurately and explicitly modeled while the secondary particles are modeled by a two scale multiresolution continuum model. The present numerical simulation captures detailed characteristics of the fracture process, such as zigzag crack morphology, critical void growth ratios, local stress triaxiality variation, and intervoid ligament structure. For the first time, fracture toughness is linked to multiscale microstructures in a realistic large 3D model.     

Grain boundary versus transgranular ductile failure

The competition between intergranular and intragranular fracture is investigated using a bilayer damage model, which incorporates the relevant microstructural features of aluminium alloys with precipitate free zones (PFZ) nearby the grain boundary. One layer represents the grain behaviour: due to precipitation, it presents a high yield stress and low hardening exponent. The other layer represents the PFZ which has the behaviour of a solid solution: it is much softer but with a much higher strain hardening capacity. In both layers, void growth and coalescence is modelled using an enhanced Gurson-type model incorporating the effects of the void aspect ratio and of the relative void spacing. The effects on the ductility (i) of the flow properties of each zone, (ii) of the relative thickness of the PFZ, and (iii) of the particles spacing and volume fraction in the PFZ are elucidated. Comparisons are made with experimental data. Based on the previous analysis, qualitative understanding of trends in the fracture toughness of aluminium alloys can be gained in order to provide a link with the thermal treatment process.     

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.     

Mode II fracture behavior of a Zr-based bulk metallic glass

Shear band formation and fracture are characterized during mode II loading of a Zr-based bulk metallic glass. The measured mode II fracture toughness, K_IIc=75±4 MPa√m, exceeds the reported mode I fracture toughness by ∼4 times, suggesting that normal or mean stresses play a significant role in the deformation process at the crack tip. This effect is explained in light of a mean stress modified free volume model for shear localization in metallic glasses. Thermal imaging of deformation at the mode II crack tip further reveals that shear bands initiate, arrest, and reactivate along the same path, indicating that flow in the shear band leads to permanent changes in the glass structure that retain a memory of the shear band path. The measured temperature increase within the shear band is a fraction of a degree. However, heat dissipation models indicate that the temperature could have exceeded the glass transition temperature for less than 1 ms immediately after the shear band formed. It is shown that this time scale is sufficient for mechanical relaxation slightly above the glass transition temperature.     

An indentation technique for estimating the energy density as fracture toughness with Berkovich indenter for ductile bulk materials

A technique is proposed to estimate the energy density as fracture toughness for ductile bulk materials with an indentation system equipped with a Berkovich indenter based on the theory of plastic deformation energy transforming into the indentation energy of fracture. With progressive increase of penetration loads, the material damage is exhibited on the effective elastic modulus. A quadratic polynomial relationship between the plastic penetration depth and penetration load, and an approximate linear relationship between logarithmic plastic penetration depth and logarithmic effective elastic modulus are exhibited by indentation investigation with Berkovich indenter. The parameter of damage variable is proposed to determine the critical effective elastic modulus at the fracture point. And the strain energy density factor is calculated according to the equations of penetration load, plastic penetration depth and effective elastic modulus. The fracture toughness of aluminum alloy and stainless steel are evaluated by both indentation tests and K_IC fracture toughness tests. The predicted S_cr values of indentation tests are in good agreement with experimental results of CT tests.     

Constraint effects in adhesive joint fracture

Constraint effects in adhesive joint fracture are investigated by modelling the adherents as well as a finite thickness adhesive layer in which a single row of cohesive zone elements representing the fracture process is embedded. Both the adhesive and the adherents are elastic-plastic with strain hardening. The bond toughness Γ (work per unit area) is equal to Γ₀+Γ_p, where Γ₀ is the intrinsic work of fracture associated with the embedded cohesive zone response and Γ_p is the extra contribution to the bond toughness arising from plastic dissipation and stored elastic energy within the adhesive layer. The parameters of the model are identified from experiments on two different adhesives exhibiting very different fracture properties. Most of the tests were performed using the wedge-peel test method for a variety of adhesives, adherents and wedge thicknesses. The model captures the constraint effects resulting from the change in Γ_p: (i) the plastic dissipation increases with increasing bond line thickness in the fully plastic regime and then decreases to reach a constant value for very thick adhesive layers; (ii) the plastic dissipation in the fully plastic regime increases drastically as the thickness of the adherent decreases. Finally, this model is used to assess a simpler approach which consists of simulating the full adhesive layer as a single row of cohesive elements.     

Potential flow model of cavitation-induced interfacial fracture in a confined ductile layer

Fracture of a thin ductile layer sandwiched between stiff substrates often results from growth and coalescence of microscopic cavities ahead of an extending crack. Cavitation induced by plastic flow in a confined, ductile layer is analyzed here to evaluate the interfacial fracture toughness of such sandwich structures. For rigid-plastic materials, a new method is proposed in which the potential flow field of a fluid is used to approximate the plastic deformation. The principle of virtual work rate is applied to determine the equivalent traction-separation law. The method is demonstrated and validated for spherically symmetric cavity growth, for which an exact solution exists. We then study in detail the growth of an initially spherical cavity in a cylindrical bar of finite length subject to uniform traction at its ends. The results show that the stress-separation curves depend strongly on initial cavity size and the strain-hardening exponent, and weakly on the nominal strain. The method has clear advantages over numerical methods, such as finite-element analysis, for parametric study of cavity growth with large plastic deformation.     

The dynamic growth of a single void in a viscoplastic material under transient hydrostatic loading

We have examined the problem of the dynamic growth of a single spherical void in an elastic-viscoplastic medium, with a view towards addressing a number of problems that arise during the dynamic failure of metals. Particular attention is paid to inertial, thermal and rate-dependent effects, which have not previously been thoroughly studied in a combined setting. It is shown that the critical stress for unstable growth of the void in the quasistatic case is strongly affected by the thermal softening of the material (in adiabatic calculations). Thermal softening has the effect of lowering the critical stress, and has a stronger influence at high strain hardening exponents. It is shown that the thermally diffusive case for quasistatic void growth in rate-dependent materials is strongly affected by the initial void size, because of the length scale introduced by the thermal diffusion. The effects of inertia are quantified, and it is demonstrated that inertial effects are small in the early stages of void growth and are strongly dependent on the initial size of the void and the rate of loading. Under supercritical loading for the inertial problem, voids of all sizes achieve a constant absolute void growth rate in the long term. Inertia first impedes, but finally promotes dynamic void growth under a subcritical loading. For dynamic void growth, the effect of rate-hardening is to reduce the rate of void growth in comparison to the rate-independent case, and to reduce the final relative void growth achieved.     

The effects of void cluster size on ductile fracture

The effects of void clustering on ductile fracture are studied by modeling a discrete set of randomly distributed clusters. Each cluster consists of four, equally-spaced, cylindrical voids. The spacing between the clusters is held constant while the spacing between the voids is varied. A Eulerian finite element program is used to numerically solve the boundary value problems. A salient feature of the previous investigations is that both the ultimate stress and the fracture strain are functions of the void distribution. In contrast, the ultimate stress remains constant while the fracture strain changes with the void cluster diameter in the current investigation.     

Energy dissipation mechanisms in ductile fracture

The objective is to investigate energy dissipation mechanisms that operate at different length scales during fracture in ductile materials. A dimensional analysis is performed to identify the sets of dimensionless parameters which contribute to energy dissipation via dislocation-mediated plastic deformation at a crack tip. However, rather than using phenomenological variables such as yield stress and hardening modulus in the analysis, physical variables such as dislocation density, Burgers vector and Peierls stress are used. It is then shown via elementary arguments that the resulting dimensionless parameters can be interpreted in terms of competitions between various energy dissipation mechanisms at different length scales from the crack tip; the energy dissipations mechanisms are cleavage, crack tip dislocation nucleation and also dislocation nucleation from a Frank-Read source. Therefore, the material behavior is classified into three groups. The first two groups are the well-known intrinsic brittle and intrinsic ductile behavior. The third group is designated to be extrinsic ductile behavior for which Frank-Read dislocation nucleation is the initial energy dissipation mechanism. It is shown that a material is predicted to exhibit extrinsic ductility if the dimensionless parameter bρ_disl^1/2 (b is Burgers vector, ρ_disl is dislocation density) is within a certain range defined by other dimensionless parameters, irrespective of the competition between cleavage and crack tip dislocation nucleation. The predictions compare favorably to the documented behavior of a number of different classes of materials.     

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.