A finite element analysis of dynamic fracture initiation by ductile failure mechanisms in a 4340 steel

In some recent dropweight impact experiments [5] with pre-notched bend specimens of 4340 steel, it was observed that considerable crack tunneling occurred in the interior of the specimen prior to gross fracture initiation on the free surfaces. The final failure of the side ligaments happened because of shear lip formation. The tunneled region is characterized by a flat, fibrous fracture surface. In this paper, the experiments of [5] (corresponding to 5 m/s impact speed) are analyzed using a plane strain, dynamic finite element procedure. The Gurson constitutive model that accounts for the ductile failure mechanisms of micro-void nucleation, growth and coalescence is employed. The time at which incipient failure was observed near the notch tip in this computation, and the value of the dynamic J-integral, J d, at this time, compare reasonably well with experiments. This investigation shows that J-controlled stress and deformation fields are established near the notch tip whenever J d , increases with time. Also, it is found that the evolution of micro-mechanical quantities near the notch root can be correlated with the time variation of J d .The strain rate and the adiabatic temperature rise experienced at the notch root are examined. Finally, spatial variations of stresses and deformations are analyzed in detail.

Dynamic growth of tensile cracks by ductile and brittle fracture mechanisms in a viscoplastic material

In this paper, a finite element analysis of steady-state dynamic crack growth under Mode I, plane strain, small-scale yielding conditions is performed in a rate dependent plastic material characterized by the over-stress model. The main objective of the paper is to obtain theoretically the dependence of dynamic fracture toughness on crack speed. Crack propagation due to a ductile (micro-void) mechanism or a brittle (cleavage) mechanism, as well as transition from one mode to another are considered. The conversion from ductile to brittle has been observed experimentally but has received very little attention using analytical methods. Local fracture criteria based on strains and stresses are used to describe ductile and brittle fracture mechanisms. The results obtained in this paper are in general agreement with micro-structural observations of mode conversion during fracture initiation. Finally, the particular roles played by material rate sensitivity and inertia are examined in some detail.

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.

Influence of anvil supports and overhangs in dynamic analysis of three point bend testing

Optimum design of dynamic fracture test rigs demands a thorough appreciation of beam vibration under impact. Analyses invariably presume rigid anvils, and neglect overhang effects. The beam response predicted analytically and numerically in this paper highlights the significant role of anvil rigidity and beam overhangs on the impact dynamics of three point bend (3PB) specimens.

Asymptotic and finite element analyses of mode III dynamic crack growth at a ductile-brittle interface

In this work, dynamic crack growth along a ductile-brittle interface under anti-plane strain conditions is studied. The ductile solid is taken to obey the J(2) flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behavior. Firstly, the asymptotic near-tip stress and velocity fields are derived. These fields are assumed to be variable-separable with a power singularity in the radial coordinate centered at the crack tip. The effects of crack speed, strain hardening of the ductile phase and mismatch in elastic moduli of the two phases on the singularity exponent and the angular functions are studied. Secondly, full-field finite element analyses of the problem under small-scale yielding conditions are performed. The validity of the asymptotic fields and their range of dominance are determined by comparing them with the results of the full-field finite element analyses. Finally, theoretical predictions are made of the variations of the dynamic fracture toughness with crack velocity. The influence of the bi-material parameters on the above variation is investigated.

Constraint loss under dynamic loading in rate independent plastic solids

The objectives of this paper are to examine the loss of crack tip constraint in dynamically loaded fracture specimens and to assess whether it can lead to enhancement in the fracture toughness at high loading rates which has been observed in several experimental studies. To this end, 2-D plane strain finite element analyses of single edge notched (tension) specimen and three point bend specimen subjected to time varying loads are performed. The material is assumed to obey the small strain J(2) flow theory of plasticity with rate independent behaviour. The results demonstrate that a valid J-Q field exists under dynamic loading irrespective of the crack length and specimen geometry. Further, the constraint parameter Q becomes strongly negative at high loading rates, particularly in deeply cracked specimens. The variation of dynamic fracture toughness K-dc with stress intensity rate K for cleavage cracking is predicted using a simple critical stress criterion. It is found that inertia-driven constraint loss can substantially enhance K-dc for (K) over dot > 10(5) MPa rootm/s.

Influence of martensite content and morphology on the toughness and fatigue behavior of high-martensite dual-phase steels

A series of high-martensite dual-phase (HMDP) steels exhibiting a 0.3 to 0.8 volume fraction of martensite (V m ), produced by intermediate quenching (IQ) of a vanadium and boron-containing microalloyed steel, have been studied for toughness and fatigue behavior to supplement the contents of a recent report by the present authors on the unusual tensile behavior of these steels. The studies included assessment of the quasi-static and dynamic fracture toughness and fatigue-crack growth (FCG) behavior of the developed steels. The experimental results show that the quasi-static fracturetoughness (K ICV ) increases with increasing V m in the range between V m =0.3 and 0.6 and then decreases, whereas the dynamic fracture-toughness parameters (K ID , K D , and J ID ) exhibit a significant increase in their magnitudes for steels containing 0.45 to 0.60 V m before achieving a saturation plateau. Both the quasi-static and dynamic fracture-toughness values exhibit the best range of toughnesses for specimens containing approximately equal amounts of precipitate-free ferrite and martensite in a refined microstructural state. The magnitudes of the fatigue threshold in HMDP steels, for V m between 0.55 and 0.60, appear to be superior to those of structural steels of a similar strength level. The Paris-law exponents (m) for the developed HMDP steels increase with increasing V m , with an attendant decrease in the pre-exponential factor (C).

A comparative study of dynamic,ductile fracture initiation in two specimen configurations

In this work a single edge notched plate (SEN(T)) subjected to a tensile stress pulse is analysed, using a 2D plane strain dynamic finite element procedure. The interaction of the notch with a pre-nucleated hole ahead of it is examined. The background material is modelled by the Gurson constitutive law and ductile failure by microvoid coalescence in the ligament connecting the notch and the hole is simulated. Both rate independent and rate dependent material behaviour is considered. The notch tip region is subjected to a range of loading rates j by varying the peak value and the rise time of the applied stress pulse. The results obtained from these simulations are compared with a three point bend (TPB) specimen subjected to impact loading analysed in an earlier work [3] The variation of J at fracture initiation, J(c), with average loading rate j is obtained from the finite element simulations. It is found that the functional relationship between J(c) and j is fairly independent of the specimen geometry and is only dependent on material behaviour.

Dynamic initiation toughness measu~ment using double Grooved discs

Among the multitude of test specimen geometries used for dynamic fiacture toughness evaluation, the most widely uscd specimen is lhc Chavpy specimen due its simple geomclry and availability of testing machines. The standard Chatpy specimen dimensions may llOl always give plane st~ain condilions and hence, it may be necessary Io coilduct lcs/s using specimens of dillEvcnt thicknesses to establish the plane strain K~a. An axisymmct/ic specimen, on the otlaev hand would always give flow constraints l~n a nominal specimen thickness i~rcspcctive of the test matctial. The notched disk specimen pVOl)oscd by Bcrn:ud ctal. [1] for static and dynamic initiation toughness measurement although p~ovicles plain-strain conditions, the crack plopagatcs at an angle to the direction of applied load. This makes inteq~retation of the test results difficult us it ~Ccluivcs ~actial slices to be cut fiom the fractured specimen to ascertain the angle o1 crack growth and a linite element model l~)r tl);t{ pa~ticulat ctack o~icntalion.

Dynamic strain ageing in Ni-base superalloy 720Li

An experimental investigation into the dynamic strain ageing (DSA) of a wrought Ni-base superalloy 720Li was conducted. Characteristics of jerky, flow have been studied at intermediate temperatures of 350, 400 and 450 degrees C at strain-rates between 10(-3) and 10(-5) s(-1). Serrations of Type C are predominant within the temperature/strain-rate range explored. The major characteristics of the serrations-i.e. (a) critical plastic strain for onset of serrations, epsilon(c); (b) average stress decrement, Delta sigma(avg); and (c) strain increment between serrations. Delta epsilon(BS)-have been examined at selected temperatures and strain-rates. Negative strain-rate sensitivity was observed in the DSA regime. However. temperature did not influence tensile properties such as yield strength, ultimate strength. elongation, reduction in area, and work hardening rate or fracture features in DSA regime. Analysis of the results Suggests that locking of the mobile dislocations by substitutional alloying elements is responsible for the DSA in alloy 720Li.

Modeling of dynamic material behavior in hot deformation: Forging of Ti-6242

A new method of modeling material behavior which accounts for the dynamic metallurgical processes occurring during hot deformation is presented. The approach in this method is to consider the workpiece as a dissipator of power in the total processing system and to evaluate the dissipated power co-contentJ = ∫o σ ε ⋅dσ from the constitutive equation relating the strain rate (ε) to the flow stress (σ). The optimum processing conditions of temperature and strain rate are those corresponding to the maximum or peak inJ. It is shown thatJ is related to the strain-rate sensitivity (m) of the material and reaches a maximum value(J max) whenm = 1. The efficiency of the power dissipation(J/J max) through metallurgical processes is shown to be an index of the dynamic behavior of the material and is useful in obtaining a unique combination of temperature and strain rate for processing and also in delineating the regions of internal fracture. In this method of modeling, noa priori knowledge or evaluation of the atomistic mechanisms is required, and the method is effective even when more than one dissipation process occurs, which is particularly advantageous in the hot processing of commercial alloys having complex microstructures. This method has been applied to modeling of the behavior of Ti-6242 during hot forging. The behavior of α+ β andβ preform microstructures has been exam-ined, and the results show that the optimum condition for hot forging of these preforms is obtained at 927 °C (1200 K) and a strain rate of 1CT•3 s•1. Variations in the efficiency of dissipation with temperature and strain rate are correlated with the dynamic microstructural changes occurring in the material.

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 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.

Dynamic analysis of three point bend specimens under impact

A simple one dimensional inertial model is presented for transient response analysis of notched beams under impact, and extracting dynamic initiation toughness values. The model includes the effects of striker mass interactions, and contact deformations of the beam. Displacement time history of the striker mass is applied to the model as forcing function. The model is validated by comparison with the experimental investigation on ductile aluminium 6061 alloy and brittle polymer, PMMA.

Flow Instabilities and fracture in Ti-6Al-4V deformed in compression at 298 K to 673 K

Uniaxial compression tests were conducted on Ti-6Al-4V specimens in the strain-rate range df 0.001 to 1 s(-1) and temperature range of 298 to 673 K. The stress-strain curves exhibited a peak flow stress followed by flow softening. Up to 523 K, the specimens cracked catastrophically after the flow softening started. Adiabatic shear banding was observed in this regime. The fracture surface exhibited both mode I and II fracture features. The state of stress existing in a compression test specimen when bulging occurs is responsible for this fracture. The instabilities observed in the present tests are classified as ''geometric'' in nature and are state-of-stress dependant, unlike the ''intrinsic'' instabilities, which are dependant on the dynamic constitutive behavior of the material.

The effect of T-stress on plane strain dynamic crack growth in elastic-plastic materials

In this paper, the effects of T -stress on steady, dynamic crack growth in an elastic-plastic material are examined using a modified boundary layer formulation. The analyses are carried out under mode I, plane strain conditions by employing a special finite element procedure based on moving crack tip coordinates. The material is assumed to obey the J (2) flow theory of plasticity with isotropic power law hardening. The results show that the crack opening profile as well as the opening stress at a finite distance from the tip are strongly affected by the magnitude and sign of the T -stress at any given crack speed. Further, it is found that the fracture toughness predicted by the analyses enhances significantly with negative T -stress for both ductile and cleavage mode of crack growth.