Numerical simulations of hole growth and ductile fracture initiation under mixed-mode loading

Numerical simulations of ductile fracture initiation caused by the interaction between a notch tip and a nearby hole under mixed-mode loading involving modes I and II are performed. Attention is restricted to plane strain, small-scale yielding conditions. The Gurson constitutive model that accounts for the ductile failure mechanisms of micro-void nucleation, growth and coalescence is employed within the framework of a finite deformation plasticity theory. The failure of the ligament connecting the notch tip and the hole by either microvoid coalescence or by intense plastic strain localization is modelled. The effect of mode-mixity on the notch tip deformation, hole growth and the critical value of J at fracture initiation is examined. The dominant failure mechanism which is operative near the notch tip for various extents of mixity of modes I and II is identified.     

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.     

The fracture behaviour of notched specimens of polymethylmethacrylate

The fracture behaviour of notched specimens of polymethylmethacrylate has been examined for a wide range of geometries in Charpy impact tests, and in tensile and slow bend fracture tests. It was found that the failure of the very sharply notched specimens was consistent with linear elastic fracture mechanics and defined a constant fracture toughness K _IC for a constant notch tip radius, whereas the blunt notched specimens failed at a constant critical stress at the root of the notch.     

Effect of lattice orientation on crack tip constraint in ductile single crystals

In this work, the effect of lattice orientation on the fields prevailing near a notch tip is investigated pertaining to various constraint levels in FCC single crystals. A modified boundary layer formulation is employed and numerical solutions under mode I, plane strain conditions are generated by assuming an elastic–perfectly plastic FCC single crystal. The analysis is carried out corresponding to different lattice orientations with respect to the notch line. It is found that the near‐tip deformation field, especially the development of kink or slip shear bands is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the notch tip are also strongly influenced by the level of T ‐stress. The present results clearly establish that ductile single crystal fracture geometries would progressively lose crack tip constraint as the T ‐stress becomes more negative irrespective of lattice orientation. Also, the near‐tip field for a range of constraint levels can be characterized by two‐parameters such as K–T or J–Q as in isotropic plastic solids.     

Influence of notch root radius on ductile rupture fracture toughness evaluation with Charpy-V type specimens

The blunt notch fracture toughness of four types of carbon-manganese steel (ASTM A516 grade 70) has been determined by J-integral tests on Charpy-V type samples with different values of notch root radius, ρ. J-ρ plots, determined using specimens with a notch depth to width ratio, a/w, equal to 0.5, have shown the existence of a limiting ρ value (ρ_eff) below which applied J-intergral values at fracture initiation are constant. These ρ_eff values have been seen to depend only on second-phase particle distribution and not on their volume fraction or on the steel ferritic grain size. The procedure for deriving J-integral values at the onset of stable crack growth from J resistance curves in the case of notches has also been discussed. Experiments with Charpy specimens with a/w = 0.2 do not allow the derivation of meaningful J-ρ plots. In all cases, a ductile fracture criterion based on the constancy of the notch tip strain at rupture initiation has been proved when ρ >ρ_eff.     

Mixed-mode fracture toughness testing of concrete beams in three-point bending

Results obtained for mixed-mode fracture toughness parameters K _c , G _c , J _c , G _F (plane strain mixed-mode stress-intensity factor, energy release rate, J-integral and fracture energy, respectively) for small notched concrete beams in bending indicate that all these parameters decrease with x/S (x is the distance from support, S is the span) in general to values near midspan consistent with Mode I results. All the parameters except J _c vary with notch depth in a similar manner for each notch location.     

Fracture toughness and hydrogen embrittlement of pipes with notches

A technique for experimental determination of fracture toughness and hydrogen embrittlement of pipes made of API 5L X52 steel is described. The tests were performed using arc-shaped specimens with a notch cut out from pipes under the conditions of a three-point bend. The fracture toughness was determined in terms of the J-integral and the stress intensity factor at the notch tip. The value of K _ρ,c was established using the volumetric method based on the experimentally measured critical load and the results of the FEM calculation of the distribution of elastic-plastic stresses ahead of the notch tip, and J _ρ,c was determined using the method of separable functions. The effect of hydrogen embrittlement was studied using electrolytically prehydrogenated specimens.     

Calculation of the energy J-integral for bodies with notches and cracks

The approximate solutions for calculation of the energy J-integral of a body both with a notch and with a crack under elastic-plastic loading have been obtained. The crack is considered as the limit case of a sharp notch. The method is based on stress concentration analysis near a notch/crack tip and the modified Neuber's approach. The HRR-model and the method based on an equation of equilibrium were also employed to calculate the J-integral. The influence of the strain hardening exponent on the J-integral is discussed. New aspects of the two-parameter J ^* _c-fracture criterion for a body with a short crack are studied. A theoretical investigation of the effect of the applied critical stress (or the crack length) on the strain fields ahead of the crack tip has been carried out.     

Finite element simulations of notch tip fields in magnesium single crystals

Recent experiments using three point bend specimens of Mg single crystals have revealed that tensile twins of \(\{10\bar{1}2\}\) -type form profusely near a notch tip and enhance the fracture toughness through large plastic dissipation. In this work, 3D finite element simulations of these experiments are carried out using a crystal plasticity framework which includes slip and twinning to gain insights on the mechanics of fracture. The predicted load–displacement curves, slip and tensile twinning activities from finite element analysis corroborate well with the experimental observations. The numerical results are used to explore the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning is rationalized from the variation of normal stress ahead of the notch tip. Further, deflection of the crack path at twin–twin intersections observed in the experiments is examined from an energy standpoint by modeling discrete twins close to the notch root.     

Effect of notch root radius on fracture toughness of Ti–18Al–8Nb alloy

Abstract

The effect of notch root radius on the mode I fracture toughness of Ti–18Al–8Nb alloy in beta solution treated and water quenched condition was investigated. The apparent fracture toughness K _IA was found to be independent of the notch root radius below a critical notch root radius ρ ₀ and subsequently increase linearly with the square root of notch root radius ρ^1/2 beyond ρ ₀. The critical notch root radius in this alloy was found to be ～50 μm. The results were explained on the basis of strain controlled fracture model.     

Singular and non-singular approaches for predicting fatigue crack growth behavior

In this work, three classes of mechanisms that can cause load sequence effects on fatigue crack growth are discussed: mechanisms acting before, at or after the crack tip. After reviewing the crack closure idea, which is based on what happens behind the crack tip, quantitative models are proposed to predict the effects at the crack tip due to crack bifurcation. To predict the behavior ahead of the crack tip, a damage accumulation model is proposed. In this model, fatigue cracking is assumed caused by the sequential failure of volume elements or tiny εN specimens in front of the crack tip, calculated by damage accumulation concepts. The crack is treated as a sharp notch with a small, but not zero radius, avoiding the physically unrealistic singularity at its tip. The crack stress concentration factor and a strain concentration rule are used to calculate the notch root strain and to shift the origin of a modified HRR field, resulting in a non-singular model of the strain distribution ahead of the crack tip. In this way, the damage caused by each load cycle, including the effects of residual stresses, can be calculated at each element ahead of the crack tip using the correct hysteresis loops caused by the loading. The proposed approach is experimentally validated and extended to predict fatigue crack growth under variable amplitude loading, assuming that the width of the volume element broken at each cycle is equal to the region ahead of the crack tip that suffers damage beyond its critical value. The reasonable predictions of the measured fatigue crack growth behavior in steel specimens under service loads corroborate this simple and clear way to correlate da/dN and εN properties.     

THE ENGINEERING FLAW ASSESSMENT METHOD (EFAM)

This communication gives an introduction to the Engineering Flaw Assessment Method (EFAM) presently being developed at GKSS. The EFAM consists of various documents describing experimental procedures for determining fracture properties under various conditions, as well as analytical procedures for estimating the crack tip opening displacement in terms of δ₅ , the rate of δ₅ , (dδ₅ /dt), and the J-integral as driving force parameters.     

Surface effects on mode-I crack tip fields: A numerical study

Surface energy often significantly influences the deformation and failure behavior of materials and devices at the nanoscale. However, how it alters the local deformation around a crack tip remains unclear. In the present paper, we investigate the surface effects on the near-tip fields of a mode-I blunt crack (or notch). The theory of surface elasticity is incorporated into the finite element method. It is found that when the curvature radius of the crack root shrinks to nanometers, surface effects considerably affect the local stress distributions near the crack tip. We also calculate the J-integral, which is almost independent of surface effects except when the integral path approaches the crack tip. This demonstrates that surface effects are localized in a small zone around the crack tip, where the classical fracture mechanics solutions neglecting surface effects should be modified.     

The influence of weld strength mismatch on crack-tip constraint in single edge notch bend specimens

Previous work by Dodds and Anderson provides a framework to quantify finite size and crack depth effects on cleavage fracture toughness when failure occurs at deformation levels where J no longer uniquely describes the state of stresses and strains in the vicinity of the crack tip. Size effects on cleavage fracture are quantified by defining a value termed J _SSY: the J to which an infinite body must be loaded to achieve the same likelihood of cleavage fracture as in a finite body. In weld metal fracture toughness testing, mismatch between weld metal and baseplate strength can alter deformation patterns, which complicate size and crack depth effects on cleavage fracture toughness. This study demonstrates that there is virtually no effect of ±20 percent mismatch on J _SSYif the distance from the crack tip to the weld/plate interface (L _min) exceeds 5 mm. At higher levels of overmatch (50 to 100%), it is no longer possible to parameterize the departure of J _SSYfor a weldment from that for a homogeneous SE(B) based on L _min alone. Weld geometry significantly influences the accuracy with which J _SSYfor a welded SE(B) can be approximated by J _SSYfor a homogeneous specimen at these extreme overmatch levels.     

DAMAGE MECHANICS APPROACH TO REMOVE THE CONSTRAINT DEPENDENCE OF ELASTIC–PLASTIC FRACTURE TOUGHNESS

It is now generally agreed that the applicability of a one-parameter J-based ductile fracture approach is limited to so-called high constraint crack geometries, and that the elastic-plastic fracture toughness J_1c, is not a material constant but strongly specimen geometry constraint-dependent. In this paper, the constraint effect on elastic-plastic fracture toughness is investigated by use of a continuum damage mechanics approach. Based on a new local damage theory for ductile fracture(proposed by the author) which has a clear physical meaning and can describe both deformation and constraint effects on ductile fracture, a relationship is described between the conventional elastic-plastic fracture toughness, J_1c, and crack tip constraint, characterized by crack tip stress triaxiality T. Then, a new parameter J_dc (and associated criterion, J_d=J_dc) for ductile fracture is proposed. Experiments show that toughness variation with specimen geometry constraint changes can effectively be removed by use of the constraint correction procedure proposed in this paper, and that the new parameter J_dc is a material constant independent of specimen geometry (constraint). This parameter can serve as a new parameter to differentiate the elastic-plastic fracture toughness of engineering materials, which provides a new approach for fracture assessments of structures. It is not necessary to determine which laboratory specimen matches the structural constraint; rather, any specimen geometry can be tested to measure the size-independent fracture toughness J_dc. The potential advantage is clear and the results are very encouraging.     

An analytical expression for the J2‐integral of an interfacial crack in orthotropic bimaterials

This paper presents a new analytical expression relating the J₂‐integral and stress intensity factors (SIF) in an in‐plane traction‐free crack between two orthotropic elastic solids using the complex function method. The singular oscillatory near tip field of a bimaterial interfacial crack is usually characterized by a pair of SIFs. In linear elastic interfacial fracture mechanics, the majority of numerical and experimental methods rely on the analytical equations relating J_k‐integrals and SIFs. Although an analytical equation relating J₁‐integral or strain energy release rate and SIFs is available, a similar relation for J₂‐integral in debonded anisotropic solids is non‐existent. Using this new analytical expression, in conjunction with the values of J_k, the SIFs can be computed without the need for an auxiliary relation. An example with known analytical solutions for SIFs is presented to show the variation of the J₂‐integral near the crack tip of a bimaterial orthotropic plate. Different bimaterial combinations are considered, and the effect of material mismatch on J_k is demonstrated.     

State‐of‐the‐art review on extended stress intensity factor concepts

The stress intensity factor concept for describing the stress field at pointed crack or slit tips is well known from fracture mechanics. It has been substantially extended since Williams' basic contribution (1952) on stress fields at angular corners. One extension refers to pointed V‐notches with stress intensities depending on the notch opening angle. The loading‐mode‐related simple notch stress intensity factors K₁, K₂ and K₃ are introduced. Another extension refers to rounded notches with crack shape or V‐notch shape in two variants: parabolic, elliptic or hyperbolic notches (‘blunt notches’) on the one hand and root hole notches (‘keyholes’ when considering crack shapes) on the other hand. Here, the loading‐mode‐related generalised notch stress intensity factors K_1ρ, K_2ρ and K_3ρ are defined. The concepts of elastic stress intensity factor, notch stress intensity factor and generalised notch stress intensity factor are extended into the range of elastic–plastic (work‐hardening) or perfectly plastic notch tip or notch root behaviour. Here, the plastic notch stress intensity factors K_1p, K_2p and K_3p are of relevance. The elastic notch stress intensity factors are used to describe the fatigue strength of fillet‐welded attachment joints. The fracture toughness of brittle materials may also be evaluated on this basis. The plastic notch stress intensity factors characterise the stress and strain field at pointed V‐notch tips. A new version of the Neuber rule accounting for the influence of the notch opening angle is presented.     

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.     

Strain rate concentration and dynamic stress concentration for double‐edge‐notched specimens subjected to high‐speed tensile loads

Engineering plastics provide superior performance to ordinary plastics for wide range of the use. For polymer materials, dynamic stress and strain rate may be major factors to be considered when the strength is evaluated. Recently, high‐speed tensile test is being recognized as a standard testing method to confirm the strength under dynamic loads. In this study, therefore, high‐speed tensile test is analysed by the finite element method; then, the maximum dynamic stress and strain rate are discussed with varying the tensile speed and maximum forced displacement. The maximum strain rate increases with increasing the tensile speed u/t, but the strain rate concentration factor is found to be constant independent of tensile speed, which is defined as the maximum strain rate appearing at the notch root over the average nominal strain rate at the minimum section . It is found that the strain rate at the notch root depends on the dynamic stress rate at the notch root and independent of the notch root radius ρ. It is found that the difference between the static and dynamic maximum stress concentration (σ_yA,max ? σ_yA,st) at the notch root is proportional to the tensile speed when u/t = 5000 mm/s. Strain rate concentration factors are also discussed with varying the notch depth and specimen length. Based on the elastic strain rate concentration factor, the master curve is obtained useful for understanding the impact fracture of polycarbonate for the wide range of temperature and impact speed.     

Tensile tests and analyses of notched specimens fabricated from high strength steels using a generalized yield model

This article summarises the uniaxial tension tests for 20 notched bars fabricated from high strength steel Q345 specified in Chinese National Standards. The effects of the notch radius, r, and that of the notch depth ratio, d/D , on the ductility and fracture resistance of this high strength steel are examined. The experimental data are further analyzed using a generalized yield model together with an elliptical fracture stress envelope originally proposed by the first author. The experimental results demonstrate that cracks initiate at the notched section, with the fracture surface filled with many dimples and shearing marks. Specimens with a sharper notch radius (a smaller r) and a larger notch depth (a smaller d/D ratio) show poor ductility, but high fracture strength. The stress field computed from the numerical procedure including the generalized yield model indicates that the crack initiation occurs at the centre of the notched section which experiences the highest stress triaxiality ratio (σ_m/σ_seq) . As the stresses at the notched section reach the limiting values determined from the elliptical fracture criterion, macroscopic fracture failure in the notched bar occurs.