Plane Stress Crack Tip Field Around a Rapidly Growing Ductile/Rigid Interface Crack

Plane stress dynamic crack growth along a ductile/rigid interface is investigated. The ductile material is taken to be ideally plastic and obey the J₂ flow theory of plasticity. Under steady-state conditions, the asymptotic structure of the crack-tip stress, velocity and strain fields has been obtained. The study reveals that two types of crack-tip sectors exist, namely uniform and nonuniform plastic sectors and that the stress, strain and velocity fields are bounded (nonsingular) in all sectors. In a uniform sector, the rectangular Cartesian components of the stress, strain and velocity fields are constant, and there is no plastic strain accumulation. In a nonuniform sector, the stress, strain and velocity components at a point depend on the angular position of the point in the crack-tip polar coordinate system and are governed by a system of simultaneous ordinary differential equations. This is a sector plastic strains can accumulate. A general crack-tip sector assembly is obtained for a practical range of crack growth speeds. Several nontrivial families of admissible solutions of the crack-tip fields based on this general assembly of uniform and nonuniform crack-tip sectors are presented and discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

The mode III full-field solution in elastic materials with strain gradient effects

The near-tip asymptotic field and full-field solution are obtained for a mode III crack in an elastic material with strain gradient effects. The asymptotic analysis shows that, even though the near-tip field is governed by a single parameter B (similar to the mode III stress intensity factor), the near-tip field is very different from the classical K_III field; stresses have r ^-3/2 singularity near the crack tip, and are significantly larger than the classical K _III field within a zone of size l to the crack tip, where l is an intrinsic material length, depending on microstructures in the material. This high-order stress singularity, however, does not violate the boundness of strain energy around a crack tip. The parameter B of the near-tip asymptotic field has been determined for two anti-plane shear loadings: the remotely imposed classical K _III field, and the arbitrary shear stress tractions on crack faces. The mode III full-field solution is obtained analytically for an elastic material with strain gradient effects subjected to remotely imposed classical K _III field. It shows that the near-tip asymptotic field dominates within a zone of size 0.5 l to the crack tip, while strain gradient effects are clearly observed within 5l. It is also shown that the conventional way to evaluate the crack tip energy release rate would lead to an incorrect, infinite value. A new evaluation gives a finite crack tip energy release rate, and is identical to the J-integral. This revised version was published online in August 2006 with corrections to the Cover Date.     

Analytical forms of higher-order asymptotic elastic-plastic crack-tip fields in a linear hardening material under antiplane shear

An analytical study of the higher-order asymptotic solutions of the stress and strain fields near the traction-free crack tip under antiplane shear in a linear hardening material is investigated. The results show that every term of the asymptotic fields is controlled by both elasticity and plasticity and all the higher-order asymptotic fields are governed by linear nonhomogeneous equations. The first four term solutions are presented analytically and the first four terms are described by two independent parameters J and K ₂. The amplitude of the second order term solution is only dependent on the material properties, but independent of loading and geometry. This paper focuses on the case with traction-free crack surface boundary conditions. The effects of different crack surface boundary conditions, such as clamped and mixed surfaces, on the crack-tip fields are also presented. Comparison of multi-term solution with leading term solution, and finite element solution in an infinite strip with semi-infinite crack under constant displacements along the edges is provided.     

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

Dynamic fields near a crack tip growing in an elastic-perfectly-plastic solid

A complete asymptotic solution is given for the fields in the neighborhood of the tip of a steadily advancing crack in an incompressible elastic-perfectly-plastic solid.For Mode I crack growth in the plane strain condition, the following noteworthy results are revealed: (1) The entire crack tip in steady crack growth is surrounded by a plastic region, and no elastic unloading is predicted by the complete dynamic asymptotic solution. Thus, the elastic unloading region predicted by the result of neglecting the important influence of the inertia terms in the equations of motion. (2) Unlike the quasi-static solution, the dynamic solution yields strain fields with a logarithmic singularity everywhere near the crack tip. (3) The stress field varies throughout the entire crack tip neighborhood, but does display behavior which can be approximated by a constant field followed by an essentially centered-fan field and then by another constant field, especially for small crack growth speeds. Indeed, the stress field reduces to that for the stationary crack, as the crack tip velocity - measured by the Mach number, M - reduces to zero; the strain field, however, does not reduce to that for the static solution, as M vanishes. (4) There are two shock fronts emanating from the crack tip across which certain stress and strain components undergo jump discontinuities. The location of the shock fronts and the magnitude of the jumps depend on the crack growth speed. The stress jump vanishes while the strain jump becomes unbounded, as the crack tip speed goes to zero.Finally, the Mode III steady-state crack growth is reviewed and, on the basis of Mode I and Mode III results, it is concluded that ductile fracture criteria for nonstationary cracks must be based on solutions which include the inertia effects, and that for this purpose, quasi-static solutions may be inadequate. Then, a possible ductile fracture criterion is suggested and discussed.One interesting feature of the complete dynamic asymptotic solution is that, unlike the quasi-static solution, it yields the same strain singularity for all three fracture modes.     

A finite element analysis of plane strain dynamic crack growth at a ductile-brittle interface

In this work, steady, dynamic crack growth under plane strain, small-scale yielding conditions along a ductile-brittle interface is analysed using a finite element procedure. The ductile solid is taken to obey the J ₂ flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behaviour. The objectives of this work are to establish the validity of an asymptotic solution for this problem which has been derived recently [12], and to examine the effect of changing the remote (elastic) mode-mixity on the near-tip fields. Also, the influence of crack speed on the stress fields and crack opening profiles near the propagating interface crack tip is assessed for various bi-material combinations. Finally, theoretical predictions are made for the variation of the dynamic fracture toughness with crack speed for crack growth under a predominantly tensile mode along ductile-brittle interfaces. Attention is focused on the effect of mismatch in stiffness and density of the constituent phases on the above aspects.     

Analytic and numerical studies on mode I and mode II fracture in elastic-plastic materials with strain gradient effects

This paper presents a study on fracture of materials at microscale (∼1 μm) by the strain gradient theory (Fleck and Hutchinson, 1993; Fleck et al., 1994). For remotely imposed classical K fields, the full-field solutions are obtained analytically or numerically for elastic and elastic-plastic materials with strain gradient effects. The analytical elastic full-field solution shows that stresses ahead of a crack tip are significantly higher than their counterparts in the classical K fields. The sizes of dominance zones for mode I and mode II near-tip asymptotic fields are 0.3l and 0.5l,while strain gradient effects are observed within land 2l to the crack tip, respectively, where l is the intrinsic material length in strain gradient theory and is on the order of microns in strain gradient plasticity (Fleck et al., 1994; Nix and Gao, 1998; Stolken and Evans, 1997). The Dugdale–Barenblatt type plasticity model is obtained to provide an estimation of plastic zone size for mode II fracture in materials with strain grain effects. The finite element method is used to investigate the small-scale-yielding solution for an elastic-power law hardening solid. It is found that the size of the dominance zone for the near-tip asymptotic field is the intrinsic material lengthl. For mode II fracture under the small-scale-yielding condition, transition from the remote classical K _IIfield to the near-tip asymptotic field in strain gradient plasticity goes through the HRR field only when K _IIis relatively large such that the plastic zone size is much larger than the intrinsic material length l. For mode I fracture under small-scale-yielding condition, however, transition from the remote classical K _I field to the near-tip asymptotic field in strain gradient plasticity does not go through the HRR field, but via a plastic zone. This revised version was published online in July 2006 with corrections to the Cover Date.     

Crack-tip fields for matrix cracks between dissimilar elastic and creeping materials

The stress fields near the tip of a matrix crack terminating at and perpendicular to a planar interface under symmetric in-plane loading in plane strain are investigated. The bimaterial interface is formed by a linearly elastic material and an elastic power-law creeping material in which the crack is located. Using generalized expansions at the crack tip in each region and matching the stresses and displacements across the interface in an asymptotic sense, a series asymptotic solution is constructed for the stresses and strain rates near the crack tip. It is found that the stress singularities, to the leading order, are the same in each material; the stress exponent is real. The oscillatory higher-order terms exist in both regions and stress higher-order term with the order of O(r°) appears in the elastic material. The stress exponents and the angular distributions for singular terms and higher order terms are obtained for different creep exponents and material properties in each region. A full agreement between asymptotic solutions and the full-field finite element results for a set of test examples with different times has been obtained.     

Singular fields of a mode III interface crack in a power law hardening bimaterial

A high order of asymptotic solution of the singular fields near the tip of a mode III interface crack for pure power law hardening bimaterials is obtained by using the hodograph transformation. It is found that the zero order of the asymptotic solution corresponds to the assumption of a rigid substrate at the interface, and the first order of it is deduced in order to satisfy completely two continuity conditions of the stress and displacement across the interface in the asymptotic sense. The singularities of stress and strain of the zero order asymptotic solutions are –1/(n ₁+1) and –n/(n ₁+1) respectively (n=n ₁, n ₂ is the hardening exponent of the bimaterials). The applicability conditions of the asymptotic solutions are determined for both zero and first orders. It is proved that the Guo-Keer solution [23] is limited in some conditions. The angular functions of the singular fields for this interface crack problem are first expressed by closed form.     

Crack growth by grain boundary cavitation in the transient and extensive creep regimes

Crack growth due to cavity growth and coalescence along grain boundaries is analyzed under transient and extensive creep conditions in a compact tension specimen. Account is taken of the finite geometry changes accompanying crack tip blunting. The material is characterized as an elastic-power law creeping solid with an additional contribution to the creep rate arising from a given density of cavitating grain boundary facets. All voids are assumed present from the outset and distributed on a given density of cavitating grain boundary facets. The evolution of the stress fields with crack growth under three load histories is described in some detail for a relatively ductile material. The full-field plane strain finite element calculations show the competing effects of stress relaxation due to constrained creep, diffusion and crack tip blunting, and of stress increase due to the instantaneous elastic response to crack growth. At very high crack growth rates the Hui-Riedel fields dominate the crack tip region. However, the high growth rates are not sustained for any length of time in the compact tension geometry analyzed. The region of dominance of the Hui-Riedel field shrinks rapidly so that the near-tip fields are controlled by the HRR-type field shortly after the onset of crack growth. Crack growth rates under various conditions of loading and spanning the range of times from small scale creep to extensive creep are obtained. We show that there is a strong similarity between crack growth history and the behaviour of the C(t) and C _t parameters, so that crack growth rates correlate rather well with C(t) and C _t .A relatively brittle material is also considered that has a very different near-tip stress field and crack growth history.Visiting Professor, Brown University, August 1988 through December 1989.     

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.     

Dynamic crack growth along an elastoplastic bimaterial interface

Summary The near-tip stress and deformation rate fields of a crack dynamically propagating along an interface between dissimilar elastic-plastic bimaterials are presented in this paper. The elastic-plastic materials are characterised by theJ ₂-flow theory with linear plastic hardening. The solutions are assumed to be of variable-separable form with a power-law singularity in the radial direction. Two distinct solutions corresponding to the tensile and shear solutions exist with slightly different singularity strengths and very different mixities at the crack tip. The phenomenon of discrete and determinate mixities at the interfacial crack tip is confirmed in dynamic crack growth. This is not an artifact of the variable-separable solution assumption, arising from the linear-hardening material model. The dynamic crack analysis shows that the mixity of the near-tip field is mainly determined by the given material parameters and affected slightly by the crack propagation velocity. A significant variation of the mixity is observed near to the coalescing point of the tensile and shear solutions. The strength of the singularity is almost determined by the smaller strain-hardening alone, and dynamic inertia decreases the stress intensity. The asymptotic solutions reveal that the crack propagation velocity changes only the stress field of the tensile mode significantly. With increasing the crack propagation velocity, the stress singularity of the tensile solutions decreases obviously and the stress triaxiality at the tip (=0) falls considerably at the unity effective stress. These observations imply that the fracture toughness of the interface crack under tensile mode may be significantly higher than that under quasi-static conditions.     

Analysis of elastoplastic sharp notches

In this paper the elastoplastic solutions with higher-order terms for apex V-notches in power-law hardening materials have been discussed. Two-term expansions of the plane strain and the plane stress solutions have been obtained. It has been shown that the leading-order singularity approaches the value for a crack when the notch angle is not too large. In plane strain cases the elasticity does not enter the second-order solutions when the notch opening angle is too small. For a large notch angle, the two-term expansions of the plane strain near-tip fields are described by a single amplitude parameter. The plane stress solutions generally contain the elasticity terms. The boundary layer formulations based on the small-strain plasticity theory confirm that a dominance zone exists ahead of the notch tip. Finite element results give good agreement to the asymptotic solutions under both plane strain and plane stress conditions. The second-order terms cannot improve the predictions significantly. The near-tip fields are dominated by a single parameter. Finite element calculations under the finite strain J ₂-flow plasticity theory revealed that the finite strains can only affect local characterization of the asymptotic solution. The asymptotic solution has a large dominance zone around the notch tip. For an apex notch bounded to a rigid substrate the leading-order singularity falls with the notch angle significantly more slowly than in the homogeneous material. It vanishes at the notch angle about 135° for all power-hardening exponents. The elasticity effects enter the second-order solutions when the notch angle becomes large enough. The tip fields are characterized by the hydrostatic stress and the shear stress ahead of the notch.     

A finite element investigation of quasi-static and dynamic asymptotic crack-tip fields in hardening elastic-plastic solids under plane stress

This is the second half of a two-part finite element investigation of quasi-static and dynamic crack growth in hardening elastic-plastic solids under mode I plane stress, steady state, and small-scale yielding conditions. The hardening materials are assumed to obey the von Mises yield criterion and the associated flow rule, and are characterized by a Ramberg-Osgood type power-law effective stress-strain curve. The asymptotic feature of the crack-tip stress and deformation fields, and the influence of hardening and crack propagation speed on these fields as well as on the size and shape of the crack-tip active plastic zone, are addressed in detail. The results of this study strongly suggest the existence of stress and strain singularities of the type [ln(R _o/r )]^s (s>0) at r=0, where r is the distance to the crack tip and R ₀ is a length scaling parameter, which is consistent with the predictions of asymptotic analyses of variable-separable type by Gao et al. [1–4]. Difficulties in estimating the values of R ₀ and s by fitting the results of the present full-field study to the type of singularities shown above are analyzed, and quantititive differences between the results of this study and those of the asymptotic analyses are discussed. As expected, findings presented here share many similarities with those reported in the first part of this study [5] for crack growth in linear hardening solids. A prominent common feature of crack growth in these two types of hardening materials is that as the level of hardening decreases and the crack propagation speed increases, a secondary yield zone emerges along the crack surface, and kinks in the angular variations of the stress and velocity fields begin to develop near where elastic unloading is taking place.     

Elastic-plastic analysis of an arbitrarily-oriented mode III crack touching an interface

In this paper we investigate a semi-infinite crack terminating at an arbitrarily oriented interface between two elastic-plastic materials under an anti-plane shear loading. An analytical solution is first developed for general power-law hardening materials under a mode III loading. If both materials have the same hardening exponent, the formulation results in a nonlinear eigenequation which can be solved numerically. The stress singularities are determined as a function of two material constants: the hardening exponent n and parameter G which represents the relative resistance of the two materials. In addition to the power of the singularity, the stress, strain and displacement asymptotic fields are also determined. If the hardening exponents are not the same, the leading order terms of an expansion model ensure the stress continuity across the interface. The results show that the stress singularity mainly depends upon the material having the larger hardening exponent, with the highest stresses in the material having the smaller hardening exponent. By taking the hardening exponent n , the perfectly plastic bimaterial problem is studied. It has been found that if the crack lies in the less stiff material, the entirely plastic asymptotic fields around the crack tip can be determined. On the other hand, if the crack lies in the stiffer material, the crack-tip fields are partially elastic and partially plastic. For both cases, unique asymptotic fields can be determined explicitly. For those cases when the materials present a strain hardening property, different mathematical models are established.     

Effect of loading rate on dynamic fracture initiation toughness of brittle materials

Numerical simulation is carried out to investigate the effect of loading rate on dynamic fracture initiation toughness including the crack-tip constraint. Finite element analyses are performed for a single edge cracked plate whose crack surface is subjected to uniform pressure with various loading rate. The first three terms in the Williams’ asymptotic series solution is utilized to characterize the crack-tip stress field under dynamic loads. The coefficient of the third term in Williams’ solution, A ₃, was utilized as a crack tip constraint parameter. Numerical results demonstrate that (a) the dynamic crack tip opening stress field is well represented by the three term solution at various loading rate, (b) the loading rate can be reflected by the constraint, and (c) the constraint A ₃ decreases with increasing loading rate. To predict the dynamic fracture initiation toughness, a failure criterion based on the attainment of a critical opening stress at a critical distance ahead of the crack tip is assumed. Using this failure criterion with the constraint parameter, A ₃, fracture initiation toughness is determined and in agreement with available experimental data for Homalite-100 material at various loading rate.     

Asymptotic crack-tip fields in anisotropic power law material under antiplane shear

A hodograph transformation in conjunction with an appropriate affine transformation are both used to investigate the strain and stress fields near the crack tip in an anisotropic power law material under antiplane shear. Stress and strain exponents as well as angular distributions for the asymptotic stress and strain fields are obtained analytically. All the stress strain exponents are independent of material anisotropy, and the effect of material anisotropy on the asymptotic stress and strain field is discussed including higher order terms.     

J-Dominance in Mixed Mode Ductile Fracture Specimens

The asymptotic mixed mode crack tip fields in elastic-plastic solids are scaled by the J-integral and parameterized by a near-tip mixity parameter, M _{_p} . In this paper, the validity and range of dominance of these fields are investigated. To this end, small strain elastic-plastic finite element analyses of mixed mode fracture are first performed using a modified boundary layer formulation. Here, a two term expansion of the elastic crack tip field involving the stress intensity factor |K| the elastic mixity parameter M _{_e} as well as the T-stress is prescribed as remote boundary conditions. The analyses are conducted for different values of M _{_e} and the T-stress. Next, several commonly used mixed mode fracture specimens such as Compact Tension Shear (CTS), Four Point Bend (4PB), and modified Compact Tension specimen are considered. Here, the complete range of loading from contained yielding to large scale yielding is analyzed. Further, different crack to width ratios and strain hardening exponents are considered. The results obtained establish that the mixed mode asymptotic fields dominate over physically relevant length scales in the above geometries, except for predominantly mode I loading and under large scale yielding conditions. This revised version was published online in August 2006 with corrections to the Cover Date.     

Stress, deformation and damage fields near the tip of a crack in a damaged nonlinear material

The stress, strain, displacement and damage fields near the tip of a crack in a power-law hardening material with continuous damage formation under antiplane longitudinal shear loading are investigated analytically. The interaction between a major crack and distributed microscopic damage is considered by describing the effect of damage in terms of a damage variable D. A deformation plasticity theory coupled with damage and a damage evolution law are formulated. A hodograph transformation is employed to determine the singularity and angular distribution of the crack-tip quantities. Consequently, analytical solutions for the antiplane shear crack-tip fields are obtained. Effects of the hardening exponent n and the damage exponent m on the crack-tip fields are discussed. It is found that the present crack-tip stress and strain solutions for damaged nonlinear material are similar to the well-known HRR fields for virgin materials. However, damage leads to a weaker singularity of stress, and to a stronger singularity of strain compared to that for virgin materials, respectively. The stress associated with damage always falls below the HRR field for virgin material; but the distribution of strain associated with damage lies slightly above the HRR field for r/(J/0) > 1.5 while the difference becomes negligible when r/(J/0) > 2. The limiting distributions of stress and strain may indeed be given by the HRR field.     

A perturbation analysis of combined mode I and III dynamic crack propagation

Summary In the present paper the asymptotic stress and deformation fields of dynamic crack extension in materials with linear plastic hardening under combined mode I (plane strain and plane stress) and anti-plane shear loading conditions (mode III) are investigated. The governing equations of the asymptotic crack-tip fields are formulated from two groups of angular functions, one for the in-plane mode and the other for the anti-plane shear mode. It was assumed that all stresses and deformations are of separable functional forms ofr and , which represent the polar coordinates centered at the actual crack tip. Perturbation solutions of the governing equations were obtained. The singularity behavior and the angular functions of the crack-tip in-plane and the anti-plane stresses obtained from the perturbation analysis show that, regardless of the mixity of the crack-tip field and the strain-hardening, the in-plane stresses under the combined mode I and mode III conditions have stronger singularity in the whole mixed mode steady-state crack growth than that of the anti-plane shear stresses. The anti-plane shear stresses perturbed from the plane strain mode I solutions lose their singularity for small strain hardening, whereas the angular stress functions perturbed from the plane stress mode I have a nearly analogous uniform distribution feature compared to pure mode III cases. An obvious deviation from the unperturbed solution is generally to be observed under combined plane strain mode I and anti-plane mode III conditions, especially for a large Mach number in a material with small strain-hardening; but not under plane stress and mode III conditions. The crack propagation velocity decreases the singularities of both pure mode and perturbed crack-tip fields.