An analysis of dynamic,ductile crack growth in a double edge cracked specimen

Dynamic crack growth is analysed numerically for a plane strain double edge cracked specimen subject to symmetric impulsive tensile loading at the two ends. The material behavior is described in terms of an elastic-viscoplastic constitutive model that accounts for ductile fracture by the nucleation and subsequent growth of voids to coalescence. Two populations of second phase particles are represented, including large inclusions or inclusion colonies with low strength, which result in large voids near the crack tip at an early stage, and small second phase particles, which require large strains before cavities nucleate. The crack growth velocities determined here are entirely based on the ductile failure predictions of the material model, and thus the present study is free from ad hoc assumptions regarding appropriate dynamic crack growth criteria. Adiabatic heating due to plastic dissipation and the resulting thermal softening are accounted for in the analyses. Different prescribed impact velocities, inclusion spacings and values of the inclusion nucleation stress are considered. Predictions for the dynamic crack growth behavior and for the time variation of crack tip characterizing parameters are obtained for each case analyzed.     

A numerical investigation of ductile fracture initiation in a high-strength low-alloy steel

In this work, static and drop-weight impact experiments, which have been conducted using three-point bend fracture specimens of a high-strength low-alloy steel, are analysed by performing finite-element simulations. The Gurson constitutive model that accounts for the ductile failure mechanisms of microvoid nucleation, growth and coalescence is employed within the framework of a finite deformation plasticity theory. Two populations of second-phase particles are considered, including large inclusions which initiate voids at an early stage and small particles which require large strains to nucleate voids. The most important objective of the work is to assess quantitatively the effects of material inertia, strain rate sensitivity and local adiabatic temperature rise (due to conversion of plastic work into heat) on dynamic ductile crack initiation. This is accomplished by comparing the evolution histories of void volume fraction near the notch tip in the static analysis with the dynamic analyses. The results indicate that increased strain hardening caused by strain rate sensitivity, which becomes important under dynamic loading, plays a benign role in considerably slowing down the void growth rate near the notch tip. This is partially opposed by thermal softening caused by adiabatic heating near the notch tip.     

A numerical analysis of failure characteristics of ductile layers in laminated composites

In this study, the failure of the ductile layers from collinear, multiple and delaminating cracks that occur in laminated composite systems was studied using a constitutive relationship that accounts for strength degradation resulting from the nucleation and growth of voids. The results indicate that, in laminated composites, void nucleation and growth ahead of the cracks occur at a much faster rate because of evolution of much higher stress values in the interface region. Except for short crack extensions, collinear and multiple cracks develop crack resistance curves similar to that seen for a crack in the ductile layer material as a homogenous isotropic cases. For delaminating crack cases, the fracture behaviour is strongly influenced by the delamination length. The resistance of the ductile layers to crack extension can be significantly reduced by short delamination lengths; however, for large delamination lengths the resistance to crack extension becomes greater than that seen for the ductile material. The results also show that, if the crack tip is at the interface, similar maximum stress values develop in the ductile layers as in the fracture test of the same ductile material, suggesting that ductile–brittle fracture transition behaviour of the ductile layers is dependent upon the extent of the cracks in the brittle layers and fracture characteristics of the brittle layers.     

A simple damage-accumulation model for constraint effects on ductile fracture

A simple model is presented to account for the effects of void-type damage on crack initiation and propagation in ductile steels under plane strain conditions by virtue of elementary fracture mechanics solutions. Multiple primary voids from large inclusions are uniformly distributed ahead of the crack tip. The growth of these primary voids is followed by nucleation of a large population of secondary voids from second-phase particles. A critical accumulative damage based on the length ratio of the damage zone to the spacing of primary voids, is employed as a failure criterion, including contributions from two populations of voids. Damage accumulation depends much on the strain and stress states such as stress triaxilities, which are extracted from existing results instead of detailed computation. Results show the dependence of fracture toughness on the size of damage zones associated with constraints. Initiation of crack growth is insensitive to the constraints since nucleation of fine voids is determined by local deformation. The model captures the transition in mechanisms from void-by-void growth to multiple void interactions in terms of a decreasing trend in the slopes of fracture resistance curves. At high constraints and large damage zone, a steady-state crack advance is identified with constant toughness. Damage accumulation from the growth of primary voids determines subsequent crack growth resistance and the study demonstrates its dependences on the crack-tip constraints.     

Local approach to ductile fracture and dynamic strain aging

Experiments on smooth and notched round specimens on a C–Mn steel used in nuclear industry are performed at different temperatures under quasi-static loadings, revealing dynamic strain aging (DSA). The behavior is highly dependent on temperature and strain rate, and a drop in fracture strain is observed. Fracture surface observation on notched tensile specimens shows classical ductile fracture mechanisms with growth and coalescence of voids. The apparent strain hardening behavior at each temperature and strain rate is taken into account to compute the void growth with the Rice and Tracey model and with a damage law developed from unit cell computations. It is shown that the apparent strain hardening at large strains is of major importance to correctly predict fracture with the Rice and Tracey model, but its influence on the void growth law is of minor importance. In particular, the stress triaxiality ratio within the notch is increased due to the negative strain rate sensitivity. The ductility drop observed in DSA domain is then partly explained, but void nucleation and void growth in presence of strain bands should be included in the fracture modeling of such materials.     

An FEM study on crack tip blunting in ductile fracture initiation

Ductile fracture is initiated by void nucleation at a characteristic distance (I_c) from the crack tip and propagated by void growth followed by coalescence with the tip. The earlier concepts expressed I_c in terms of grain size or inter-particle distance because grain and particle boundaries form potential sites for void nucleation. However, Srinivas et al. (1994) observed nucleation of such voids even inside the crack tip grains in a nominally particle free Armco iron. In an attempt to achieve a unified understanding of these observations, typical crack-tip blunting prior to ductile fracture in a standard C(T) specimen (Mode I) was studied using a finite element method (FEM) supporting large elasto-plastic deformation and material rotation. Using a set of experimental data on Armco iron specimens of different grain sizes, it is shown that none of the locations of the maxima of the parameters stress, strain and strain energy density correspond to I_c. Nevertheless, the size of the zone of intense plastic deformation, as calculated from the strain energy density distribution ahead of the crack tip in the crack plane, compares well with the experimentally measured I_c. The integral of the strain energy density variation from the crack tip to the location of void nucleation is found to be linearly proportional to J_IC. Using this result, an expression is arrived at relating I_c to J_IC and further extended to CTOD_c.     

Effect of pre-strain on fracture toughness of ductile structural steels under static and dynamic loading

The ductile fracture process consists of void nucleation, growth and coalescence. The whole ductile process can be divided into two successive steps: (I) the initial state to void nucleation, followed by (II) void growth up to void coalescence. Based on this suggestion, resistance to ductile fracture could be divided into the resistance to stage I and stage II, and accordingly the whole fracture toughness could be regarded to be due to contributions from stages I and II. The fracture toughness contributed from the two steps is, respectively, denoted as void nucleation-contributed fracture toughness and void growth-contributed fracture toughness. The effect of plastic pre-strain on the fracture toughness of ductile structural steels under static and dynamic loading (4.9 m/s) within the ductile fracture range was evaluated by summing contributions due to void nucleation-contributed and void growth-contributed fracture toughness. The effect of strain rate on fracture toughness was also investigated by the same means. The results show that both plastic pre-strain and high-speed loading decrease the void nucleation-contributed fracture toughness while their effects on the void growth-contributed fracture toughness depend on the variations in strength and ductility. Moreover, fracture toughness of structural steels generally decreases with increasing strain rate.     

Dynamic ductile crack growth and transition to cleavage – a cell model approach

The fracture behavior of ferritic steel in the transition regime is controlled by the competition between ductile tearing and cleavage. Many test specimens that failed by catastrophic cleavage showed significant amounts of ductile tearing prior to cleavage fracture. The transition from ductile tearing to cleavage has been attributed to the increase in constraint and sampling volume associated with ductile crack growth. This work examines the role of dynamic ductile crack growth on the fracture mode transition by way of a cell model of the material. The cell model incorporates the effects of stress triaxiality and strain rate on material failure characteristics of hole growth and coalescence. Loading rate and microstructure effects on the stress fields that evolve with rapid (ductile) crack growth are systematically studied. The stress fields are employed to compute the Weibull stress which provides probability estimates for the susceptibility to cleavage fracture. A center-cracked panel subjected to remote tension is the model problem under study. The computational model uses an elastic-viscoplastic constitutive relation which incorporates enhanced strain rate hardening at high strain rates. Adiabatic heating due to plastic dissipation and the resulting thermal softening are also accounted for. Under dynamically high loading rate, our model shows the crack speed achieves its peak value soon after crack initiation and quickly falls off to slower speeds with further crack growth. Remarkably, the Weibull stress follows a similar pattern which suggests that the transition to the cleavage fracture is most likely to occur, if at all, at the peak speed of ductile crack growth. Key words: Dynamic fracture, ductile tearing, crack growth, transition regime, cleavage fracture, cell model, finite element.     

MICROMECHANICAL MODELLING OF DAMAGE AND FRACTURE OF DUCTILE MATERIALS

A micromechanical model of ductile damage by void nucleation, growth and coalescence is widely and successfully applied to describe phenomena of ductile tearing. The model's fundamental principles, and especially the constitutive equations of Gurson, Tvegaard and Needleman (GTN-model), are briefly described. Some of the material parameters of the GTN-model are calibrated by performing cell model calculations, which is a method of determining the structural behaviour of a single void in a plastic material. The approach is used to study the dependence of material strength and toughness on microstructural features of nodular cast iron. The numerical simulations were realized within the FE-program ABAQUS by a user-supplied material model.     

Numerical simulations of fracture initiation in ductile solids under mode I, dynamic loading

In this paper, an overview of some recent computational studies by the authors on ductile crack initiation under mode I, dynamic loading is presented. In these studies, a large deformation finite element procedure is employed along with the viscoplastic version of the Gurson constitutive model that accounts for the micro-mechanical processes of void nucleation, growth and coalescence. A three-point bend fracture specimen subjected to impact, and a single edge notched specimen loaded by a tensile stress pulse are analysed. Several loading rates are simulated by varying the impact speed or the rise time and magnitude of the stress pulse. A simple model involving a semi-circular notch with a pre-nucleated circular hole situated ahead of it is considered. The growth of the hole and its interaction with the notch tip, which leads to plastic strain and porosity localization in the ligament connecting them, is simulated. The role of strain-rate dependence on ductile crack initiation at high loading rates, and the specimen geometry effect on the variation of dynamic fracture toughness with loading rate are investigated.     

Ductile to brittle transition of an A508 steel characterized by Charpy impact test: Part I: experimental results

This study is devoted to the ductile-brittle transition behavior of a French A508 Cl3 (16MND5) steel. Due to its importance for the safety assessment of PWR vessels, a full characterization of this steel with Charpy V-notch test in this range of temperature was undertaken. The aim of this study is to provide a wide experimental database and microstructural observations to supply, calibrate and validate models used in a local approach methodology. Mechanical and fracture properties of the steel have been investigated over a wide range of temperatures and strain-rates. Effects of impact velocity on ductile-brittle transition curve, on ductile tearing and on notch temperature rise are presented and discussed. A detailed study of ductile crack initiation and growth in Charpy specimens is also carried out. From fractographic investigations of the microvoids nucleation around carbide second phase particles, a plastic strain threshold for nucleation is determined for this material. A508 Cl3 steels undergo a transition in fracture toughness properties with temperature, due to a change in fracture mode from microvoids coalescence to cleavage fracture. A systematic investigation on the nature and the position of cleavage triggering sites and on any change in the ductile to brittle transition (DBT) range has been carried out. This leads to the conclusion that manganese sulfide inclusions do not play an increasing role with increasing test temperature as recently mentioned in other studies on A508 Cl3 steel with a higher sulfur content. In a companion paper [Tanguy et al., Engng. Fract. Mech., in press], the numerical simulation of the Charpy test in the ductile-brittle transition range using fully coupled local approach to fracture is presented.     

A continuum damage mechanics model for void growth and micro crack initiation

A damage mechanics model is proposed to study the void growth and crack initiation. J₂ incremental flow theory along with a damage variable is used to model the material behaviour in elasto-plastic regime. Large deformation (large rotation and finite strain) finite element analysis is carried out for five different cases. In all the cases it is observed that the triaxiality and the plastic strain play an important role in void growth and crack initiation in ductile material. A failure curve is obtained for the material AISI-1090 spheroidised steel. Finally, it is concluded that the critical value of the damage variable can be taken as a crack initiation parameter.     

Temperature dependence of fracture toughness J IC and ductility for BCC materials in the transition region

The temperature dependence of ductility, strength and fracture toughness for a BCC material undergoing predominantly linear elastic behavior at low temperatures and elastic-plastic behavior at higher temperatures is examined. A model, based on ductile fracture mechanisms involving void nucleation followed by cavity growth and void coalescence, is developed to relate the fracture toughness parameter J _IC with temperature. Two general equations for linear elastic and elastic plastic regimes of J _IC versus temperature T, are obtained. Applications of this model to experimental data obtained on a carbon steel show that J IC varies with T ² at low temperatures and with T at higher temperatures, thus defining a transition temperature.     

Crack straining-based spall model

Mean void growth-based spall models that avoid the complications of nucleation have been successfully applied to the problem of ductile spallation. However, similar models based on mean crack growth, applicable to brittle spallation, are not promising. This is because it remains to be demonstrated how an appropriate mean crack size is chosen to identify the brittle spall strength as the threshold pressure for crack growth. In the authors' view, the solidity evolution is not merely a consequence of the nucleation and growth of cracks but also determined by the crack straining caused by the relaxed tensile pressure. This paper presents a crack straining-based spall model which assumes that the inelastic volumetric strain caused by the relaxed tensile pressure at a critical fragment volume is the main factor governing the solidity evolution during the process of coalescence and fragmentation. Such an approach makes the modelling of brittle spall compatible with that of ductile. The developed model requires only two additional parameters to those required by the chosen constitutive model and equation of state. We present the comparisons between experimental and computational simulation of the free surface velocity history of the target in a plane impact plate experiment. These show very good agreement.     

Simulation of Stable Crack Growth for Welded Joints Including Strength Mismatching

1.IotroductionStrengthmismatchingeffectsonfracturetough-nessofferritesteelsremainakeyissueforthesafetyassessmentofstructures.Structuralandpressureves-selsteelsgenerallyexhibitincreaseinfracturetough-nessoverthefirstfewdistanceofstablecrackgrowth.Laboratorytestingoffracturespecimenstomeasureresistancecurves(R-curves)consistentlyrevealsamarkedeffectofstrengthmismatchingonR-curves.Fortheweldedjoint,whoseweldmetalstrengthishigherthanthatofthebasemetal(overmatching),yieldahigherRcurvecomparingwit…     

Nanoscale void nucleation and growth and crack tip stress evolution ahead of a growing crack in a single crystal

A constrained three-dimensional atomistic model of a cracked aluminum single crystal has been employed to investigate the growth behavior of a nanoscale crack in a single crystal using molecular dynamics simulations with the EAM potential. This study is focused on the stress field around the crack tip and its evolution during fast crack growth. Simulation results of the observed nanoscale fracture behavior are presented in terms of atomistic stresses. Major findings from the simulation results are the following: (a) crack growth is in the form of void nucleation, growth and coalescence ahead of the crack tip, thus resembling that of ductile fracture at the continuum scale; (b) void nucleation occurs at a certain distance ahead of the current crack tip or the forward edge of the leading void ahead of the crack tip; (c) just before void nucleation the mean atomic stress (or equivalently its ratio to the von Mises effective stress, which is called the stress constraint or triaxiality) has a high concentration at the site of void nucleation; and (d) the stress field ahead of the current crack tip or the forward edge of the leading void is more or less self-similar (so that the forward edge of the leading void can be viewed as the effective crack tip).     

A finite element analysis of fracture initiation in ductile/brittle periodically layered composites

A near-tip plane strain finite element analysis of a crack terminating at and normal to the interface in a laminate consisting of alternate brittle and ductile layers is conducted under mode-I loading. The studies are carried out for a system representing steel/alumina composite laminate. The Gurson constitutive model, which accounts for the ductile failure mechanisms of microvoid nucleation, growth and coalescence, is employed within the framework of small deformation plasticity theory. Evolution of plastic zone and damage in the ductile layer is monitored with increasing load. High plastic strain localization and microvoid damage accumulation are found to occur along the brittle/ductile interface at the crack-tip. Fracture initiation in the ductile phase is predicted and the conditions for crack renucleation in the brittle layer ahead of the crack are established for the system under consideration. Ductile fracture initiation has been found to occur before plasticity spreads in multiple ductile layers. Effects of material mismatch and yield strength on the plastic zone evolution are briefly discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Effect of residual stresses on ductile crack growth resistance

It is well known that residual stresses influence the ductile fracture behaviour. In this paper, a numerical study was performed to assess the effect of residual stresses on ductile crack growth resistance of a typical pipeline steel. A modified boundary layer model was employed for the analysis under plane strain, Mode I loading condition. The residual stress fields were introduced into the finite element model by the eigenstrain method. A sharp crack was embedded in the center of the weld region. The complete Gurson model has been applied to simulate the ductile fracture by microvoid nucleation, growth and coalescence. Results show that tensile residual stresses can significantly reduce the crack growth resistance when the crack growth is small compared with the length scale of the tensile residual stress field. With the crack growth, the effect of residual stresses on the crack growth resistance tends to diminish. The effect of residual stress on ductile crack growth resistance seems independent of the size of geometrically similar welds. When normalized by the weld zone size, the ductile crack growth resistance collapses into one curve, which can be used to assess the structural integrity and evaluate the effect of residual stresses. It has also been found that the effect of residual stresses on crack growth resistance depends on the initial void volume fraction f₀, hardening exponent n and T-stress.