Crack tip fields in ductile crystals

Results on the asymptotic analysis of crack tip fields in elastic-plastic single crystals are presented and some preliminary results of finite element solutions for cracked solids of this type are summarized. In the cases studied, involving plane strain tensile and anti-plane shear cracks in ideally plastic f c c and b c c crystals, analyzed within conventional small displacement gradient assumptions, the asymptotic analyses reveal striking discontinuous fields at the crack tip.For the stationary crack the stress state is found to be locally uniform in each of a family of angular sectors at the crack tip, but to jump discontinuously at sector boundaries, which are also the surfaces of shear discontinuities in the displacement field. For the quasi-statically growing crack the stress state is fully continuous from one near-tip angular sector to the next, but now some of the sectors involve elastic unloading from, and reloading to, a yielded state, and shear discontinuities of the velocity field develop at sector boundaries. In an anti-plane case studied, inclusion of inertial terms for (dynamically) growing cracks restores a discontinuous stress field at the tip which moves through the material as an elastic-plastic shock wave. For high symmetry crack orientations relative to the crystal, the discontinuity surfaces are sometimes coincident with the active crystal slip planes, but as often lie perpendicular to the family of active slip planes so that the discontinuities correspond to a kinking mode of shear.The finite element studies so far attempted, simulating the ideally plastic material model in a small displacement gradient type program, appear to be consistent with the asymptotic analyses. Small scale yielding solutions confirm the expected discontinuities, within limits of mesh resolution, of displacement for a stationary crack and of velocity for quasi-static growth. Further, the discontinuities apparently extend well into the near-tip plastic zone. A finite element formulation suitable for arbitrary deformation has been used to solve for the plane strain tension of a Taylor-hardening crystal panel containing, a center crack with an initially rounded tip. This shows effects due to lattice rotation, which distinguishes the regular versus kinking shear modes of crack tip relaxation. and holds promise for exploring the mechanics of crack opening at the tip.     

Numerical analysis of blunting of a crack tip in a ductile material under small-scale yielding and mixed mode loading

The blunting of the tip of a crack in a ductile material is analysed under the conditions of plane strain, small-scale yielding, and mixed mode loading of Modes I and II. The material is assumed to be an elastic-perfectly plastic solid with Poisson's ratio being 1/2. The stress and strain fields for a sharp crack under mixed mode loading are first determined by means of elastic-plastic finite element analysis. It is shown that only one elastic sector exists around the crack tip, in contrast with the possibility of existence of two elastic sectors as discussed by Gao. The results obtained for a sharp crack are used as the boundary conditions for the subsequent numerical analysis of crack tip blunting under mixed mode loading, based on slip line theory. The characteristic shapes of the blunted crack tip are obtained for a wide range of Mode I and Mode II combinations, and found to resemble the tip of Japanese sword. Also the stress field around the blunted crack tip is determined.     

Influence of isotropic strain hardening on plane strain dynamic growth in elastic-plastic materials

In this paper, dynamic crack growth in an elastic-plastic material is analysed under mode I, plane strain, small-scale yielding conditions using a finite element procedure. The material is assumed to obey J₂ incremental theory of plasticity with isotropic strain hardening which is of the power-law type under uniaxial tension. The influence of material inertia and strain hardening on the stress and deformation fields near the crack tip is investigated. The results demonstrate that strain hardening tends to oppose the role of inertia in decreasing plastic strains and stresses near the crack tip. The length scale near the crack tip over which inertia effects are dominant also diminishes with increase in strain hardening. A ductile crack growth criterion based on the attainment of a critical crack tip opening displacement is used to obtain the dependence of the theoretical dynamic fracture toughness on crack speed. It is found that the resistance offered by the elastic-plastic material to high speed crack propagation may be considerably reduced when it possesses some strain hardening.     

The dynamic near crack-line fields for mode II crack growth in an elastic perfectly-plastic solid

The dynamic near crack-line fields for mode II crack growth in an elastic perfectly-plastic solid are investigated under plane strain and plane stress conditions. In each case, by expanding the plastic fields and the governing equations in the coordinate y, the problem is reduced to solving a system of nonlinear ordinary differential equations which is similar to that of mode III derived by Achenbach and Z.L.Li. An approximate solution for small values of x is obtained and matched with the elastic field of a blunt crack at the elastic-plastic boundary. The crack growth criterion of critical strain is employed to determine the value of K _II of the far-field that would be required for a steadily growing crack.     

Orowan mechanism for abrupt localized flow and fracture initiation in metals ahead of an yielding notch

Plane strain elastic-plastic stresses are determined in Mises yielding solid at the root of an yielding crack like notch. This external edge notch is infinitely deep, and has a small finite (fixed) flank angle with a small tip root blunting radius. A boundary value type approach has been followed throughout, to solve this famous Orowan-lrwin problem. Firstly, a fictitious elastic stress field is calculated, considering a misfit in the bulk volume loading; these elastic stress expressions are valid when the notch is fully loaded. Secondly, the plastic stresses are determined in the compressibility gradient, maintaining the continuity of stresses and their derivatives at the yielded-unyielded interface. Our calculations reveal that: Orowan mechanism is fairly dominant below the notch root, as well as on ± 45° planes. It is concluded that the flow-localization in the Mises solid is due to a reverse slip, caused by the sudden release of a favourable critical mismatch stress concentration. Some elastic strain energy density is seen to be getting released from the bulk volume, while unloading the misfit load. The mismatch has been created entirely due to the compressibility-incompressibility difference, as suggested by Orowan.

Following Orowan, it is shown here that, before the onset of a stable crack extension, the increase in stress concentration at the notch tip root, is directly proportional to the strength of mismatch strain-localization below the notch, and inversely proportional to the plane strain plastic zone size on the crack extension plane. For a large scale yielding situation, compressive stresses and pure distortion regions are seen to occur at a far field within the plastic enclave.     

Determination of K Id at or below NDTT using instrumented drop-weight testing

Two kinds of surface cracked specimens are numerically analyzed by the three dimensional elastic-plastic FEM. Near tip regions are divided into fine elements, and the stress and displacement fields at the crack tip are compared with HRR solutions.At first, surface cracked specimens subjected to bending are analyzed by changing the aspect ratio and depth/thickness ratio. The effect of the loading condition, crack shape and the crack depth on the stress and displacement fields are discussed. Then the pipe with surface crack subjected to bending is analyzed and the availability of the J-integral concept to the LBB analysis is discussed.In every specimen, it is shown that in the regions very near to the crack tip, the displacement field is similar to HRR solutions of plane strain. In the outer regions, however, the stress and displacement fields depend strongly on the shape, thickness, and loading conditions.     

Tensile crack tip fields in elastic-ideally plastic crystals

Crack tip stress and deformation fields are analyzed for tensile-loaded ideally plastic crystals. The specific cases of (0 1 0) cracks growing in the [1 0 1] direction, and (1 0 1) cracks in the [0, 1, 0] direction, are considered for both fcc and bcc crystals which flow according to the critical resolved shear stress criterion. Stationary and quasistatically growing crack fields are considered. The analysis is asymptotic in character; complete elastic-plastic solutions have not been determined. The near-tip stress state is shown to be locally constant within angular sectors that are stressed to yield levels at a stationary crack tip, and to change discontinuously from sector to sector. Near tip deformations are not uniquely determined but fields involving shear displacement discontinuities at sector boundaries are required by the derived stress state. For the growing crack both stress and displacement must be fully continuous near the tip. An asymptotic solution is given that involves angular sectors at the tip that elastically unload from, and then reload to, a plastic state. The associated near-tip velocity field then has discontinuities of slip type at borders of the elastic sectors. The rays, emanating from the crack tip, on which discontinuities occur in the two types of solutions are found to lie either parallel or perpendicular to the family of slip plane traces that are stressed to yield levels by the local stresses. In the latter case the mode of concentrated shear along a ray of discontinuity is of kink type. Some consequences of this are discussed in terms of the dislocation generation and motion necessary to allow the flow predicted macroscopically.     

An elastic-plastic finite element analysis of crack tip fields under biaxial loading conditions

A square plate containing a central crack and subjected to biaxial stresses has been studied by a finite element analysis. An elastic analysis shows that the crack opening displacement and stress of separation ahead of the crack tip are not affected by the mode of biaxial loading and therefore the stress intensity factor adequately describes the crack tip states in an elastic continuum.

An elastic-plastic analysis involving more than localized yielding at the crack tip provides different solutions of crack tip stress fields and crack face displacements for the different modes of biaxial loading.

The equi-biaxial loading mode causes the greatest separation stress but the smallest plastic shear ear and crack displacement. The shear loading system induces the maximum size of shear ear and crack displacement but the smallest value of crack tip separation stress.

     

Higher order asymptotic solutions of V-notch tip fields for damaged nonlinear materials under antiplane shear loading

The higher order solutions of stress and deformation fields near the tip of a sharp V-notch in a power-law hardening material with continuous damage formation are analytically investigated under antiplane shear loading condition. The interaction between a macroscopic sharp notch and distributed microscopic damage is considered by describing the effect of damage in terms of a damage variable in the framework of damage mechanics. A deformation plasticity theory coupled with damage and a damage evolution law are formulated. A hodograph transformation is employed to determine the solution of damaged nonlinear notch problem in the stress plane. Then, inversion of the stress plane solution to the physical plane is performed. Consequently, higher order terms in the asymptotic solutions of the notch tip fields are obtained. Analytical expressions of the dominant and second order singularity exponents and associated angular distribution functions of notch tip stress and strain are presented. Effects of damage and strain hardening exponents and notch angle on the singular behavior of the notch tip quantities are discussed detailly. It is found that damage can lead to a weaker singularity of the dominant term of stress on one hand, but to stronger singularities of the second order term of stress and the dominant and second order terms of strain compared to that for undamaged case on the other. Also, both hardening exponent and notch angle have important effects on the notch tip quantities. Moreover, reduction of the notch tip solutions to a damaged nonlinear crack problem is carried out, and higher order solutions of the crack tip fields are obtained. Effects of damage and hardening exponents on the dominant and second order terms in the crack tip solutions are detailly discussed. Discussions on some other special cases are also presented, which shows that if damage exponent equals to zero, then the present solutions can be easily reduced to the solutions for undamaged cases. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

An interface crack in a coating-substrate composite with mixed-mode Dugdale corrections

Considering both plane stress and plane strain conditions, the plastic zone size and the crack tip opening displacement of an interface crack between a coating and a semi-infinite substrate under a normal load on the crack surfaces are investigated by the mixed-mode Dugdale model. In the model, stresses applied in the plastic zones satisfy the Von Mises yield criterion. The plastic zone size can be calculated by satisfying the condition that the complex stress intensity factors vanish. After the plastic zone size is solved, the crack tip opening displacement can be obtained by dislocation theories. In numerical examples, a uniform load is considered, and the effects of the normalized elastic modulus (the ratio of the elastic modulus of the coating to the elastic modulus of the substrate) and the normalized crack depth (the ratio of the coating thickness to the interface crack length) on the normalized plastic zone size and the normalized crack tip opening displacement are examined. Numerical examples show in the case of thin coatings, the value of the normalized plastic zone size decreases with increasing the normalized elastic modulus.     

Energy density approach for calculation of three-dimensional inelastic strain and stress at the crack tip in compact tension specimens of polycarbonate

An energy-based method is utilized for calculating elastic-plastic strains and stresses near fatigue crack tip in specimens of Merlon polycarbonate. The stress redistribution caused by the plastic yielding around the crack tip is taken into account so that theoretical crack tip strain is improved. The estimated values of crack tip strain based on an energy density approach are compared with experimental results obtained from an embedded grid moire technique and embedded strain gages. Large-scale yielding seems to dominate near the crack tip. In fact, the measured strain is in agreement with the elastic solution, which means, in reality, only small-scale yielding takes place near the crack tip. The strain in the mid-plane (plane strain) is found to be higher than in the surface plane (plane stress). The experimental and theoretical results are in good agreement.     

Fatigue crack closure study based on whole-field displacements

This study is concerned with crack tip strain field fluctuations at loads below the point of crack closure in fatigue cycling. Moiré interferometry was used to investigate crack tip fields in compact tension specimens, cracked under constant stress intensity range and fixed R-ratio conditions. An elastic-plastic finite element model of simulated closure was developed to provide a theoretical cross-reference for the moiré studies. The ‘stretched zone’, which is believed to be the most significant source of closure effects, was simulated by inserting a constant thickness strip of elements into the crack before unloading from the maximum load point. Analysis of the crack tip fields in the experimental and theoretical cases was made in terms of crack face opening profiles, compliance changes and elastic stress intensity parameters. The latter were inferred through stress and displacement measurements made along circular and radial paths relative to the crack tip. Closure on the stretched zone was found to generate non-proportional loading in the crack tip field, so that the resulting stress changes were not well characterized by the asymptotic elastic equations. It is concluded firstly, that significant strain fluctuations occur below the point of closure load and that these should not be ignored in crack propagation studies. Secondly, the effective stress intensity range in fatigue cycling is not simply related to the open-crack stress intensity range and the need therefore remains for R-ratio and geometry effects to be treated as variables in crack propagation data collection programmes.     

Evaluation of crack tip stress and strain fields under nonproportional loading in a viscoelastic material

The crack tip strain and stress fields in a viscoelastic material under nonproportional loading conditions are evaluated. In order to evaluate the strain field, the crack tip displacement field is measured by using the digital image correlation (DIC) technique. This displacement field is then approximated by using the theoretically obtained crack tip displacement field in viscoelastic materials. The result shows that the approximation method can smoothly reconstruct the experimentally obtained displacement field. From the approximated displacement field, the crack tip strain field can be precisely obtained by using the differential form of the theoretical displacement. On the other hand, the crack tip stress field is analyzed by using the stress function. This suggests that the strain and stress fields can be independently evaluated. In addition, different time dependencies between stress and strain fields near the crack tip are observed. Based on this experiment, we can discuss the several criteria for the crack propagation directions in viscoelastic materials.     

A finite element analysis of plane strain dynamic crack growth in materials displaying the Bauschinger effect

In this paper dynamic crack growth in an elastic-plastic material is analyzed under mode I plane strain small-scale yielding conditions using a finite element procedure. The main objective of this paper is to investigate the influence of anisotropic strain hardening on the material resistance to rapid crack growth. To this end, materials that obey an incremental plasticity theory with linear isotropic or kinematic hardening are considered. A detailed study of the near-tip stress and deformation fields is conducted for various crack speeds. The results demonstrate that kinematic hardening does not oppose the role of inertia in decreasing the plastic strains and stresses near the crack tip with increase in crack speed to the same extent as isotropic strain hardening. A ductile crack growth criterion based on the attainment of a critical crack opening displacement at a small micro-structural distance behind the tip is used to obtain the dependence of the theoretical dynamic fracture toughness with crack speed. It is found that for any given level of strain hardening, the dynamic fracture toughness displays a much more steep increase with crack speed over the quasi-static toughness for the kinematic hardening material as compared to the isotropic hardening case.     

The range of applicability of the small-scale yielding approach in theoretical investigations of yield and fracture at a crack tip

Theoretical analysis of the processes of plastic deformation and fracture in the vicinity of a sharp crack tip becomes easier if plasticity is confined to the immediate vicinity of the crack tip, for linear elastic stress intensity factors can be used to give the surrounding elastic field, and it is possible to employ a boundary layer formulation to determine the elastic-plastic field in the immediate vicinity of the crack tip. Consideration is given to the range of applicability of the small-scale yielding approach using a simple model in which slip is confined to the crack plane. A range of stress-displacement laws is considered, the object being to assess the effect of a material's flow characteristics on the range of applicability. The conclusions are discussed with respect to plane strain crack propagation in mild steel.     

Interfacial crack-tip constraints and <Emphasis Type="Italic">J</Emphasis>-integral for bi-materials with plastic hardening mismatch

This paper investigates interfacial crack tip stress fields and the J-integral for bi-materials with plastic hardening mismatch via detailed elastic-plastic finite element analyses. For small scale yielding, the modified boundary layer formulation with the elastic T-stress is employed. For fully plastic yielding, plane strain single-edge- cracked specimens under pure bending are considered. Interfacial crack tip stress fields are explained by modified Prandtl slip-line fields. It is found that, for bi-materials consisting of two elastic-plastic materials, increasing plastic hardening mismatch increases both crack-tip stress constraint in the lower hardening material and the J-contribution there. The implication of asymmetric J-integral in bi-materials is also discussed.     

The strain-hardening effect in HRR plane fields according to T-criterion

The strain-hardening effect on fracture is investigated with the aid of the T-criterion using HRR stress fields [1–3] around a crack tip in a power hardening material. Using the appropriate components of strain energy density for the elastic-plastic as well as a nearly elastic expression of the T-criterion, we find the fracture angles, as well as fracture stresses in materials possessing an elastic extended up to a perfectly plastic behavior, by considering plane mixed-mode deformation at the crack tip.Significant influence of the strain hardening coefficient, n on the fracture stress, as well as the hardening parameter mainly appeared in plane strain conditions. This phenomenon was observed almost independently of the solution applied, which provides a nearly elastic or an elastic-plastic expression of the T-criterion describing the fracture conditions.     

Stress analysis of sharply notched plates and measurement of notch tip blunting

An experimental stress analysis technique for the determination of stress intensity acting at a crack tip has been developed and applied to several problems. The approach uses precise, continuous mechanical measurement of the small changes in separation of the crack faces near the crack tip during loading. From these results and the opening mode linear elastic displacement field equation the stress intensities are calculated.

The technique has been used to measure the change in stress intensity acting on a crack progressing along a line of rivet holes and to determine the interaction effects in a parallel crack array. Agreement with the available theoretical solutions is good.

The technique has also been used to explore the onset of elastic-plastic behavior at the crack tip. It has been shown that the deviation from elastic behaviour can be analyzed in terms of equivalent elastic response to determine the crack tip dislocation as a function of load. The results are in good agreement with the relatively crude theoretical models now available.     

Elastic-plastic stress-strain calculation in notched bodies subjected to non-proportional loading

An analytical method for calculating notch tip stresses and strains in elastic-plastic isotropic bodies subjected to non-proportional loading sequences is presented. The key elements of the two proposed models are generalized relationships between elastic and elastic-plastic strain energy densities, and the material constitutive relations. These two models form the lower and the upper limits of the actual energy densities at the notch tip. Each method consists of a set of seven linear algebraic relations that can easily be solved for elastic-plastic strain and stress increments, knowing the hypothetical notch tip elastic stress history and the material stress-strain curve. Results of the validation show that the proposed methods compare well with finite element data and each solution set forms the limits of a band within which actual notch tip strains fall.