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.     

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.     

Effects of non-singular stresses on crack-tip fields for pressure-sensitive materials,Part 1: Plane strain case

Mode I near-tip stress fields for elastic perfectly plastic pressure-sensitive materials under plane strain and small-scale yielding conditions are presented. A Coulomb-type yield criterion described by a linear combination of the effective stress and the hydrostatic stress is adopted in the analysis. The finite element computational results sampled at the distance of a few crack opening displacements from the tip show that, as the pressure sensitivity increases, the magnitudes of the normalized radial and hoop stress ahead of the tip decrease, the total angular span of the singular plastic sectors decreases, and the angular span of the elastic sectors bordering the crack surfaces increases. When non-singular T stresses are considered along the boundary layer of the small-scale yielding model, the near-tip stresses decrease as the T stress decreases. The plastic zone shifts toward the crack surfaces as the T stress increases. When the discontinuities of the radial stress and the out-of-plane normal stress along the border between the plastic sector and the elastic sector are allowed, the angular variations of the asymptotic crack-tip fields agree well with those of the finite element computations. Variation of the Q stresses for pressure-sensitive materials can be found from the asymptotic solutions when the plastic zone size ahead of the tip is relatively larger than the crack opening displacement. In addition the T stress is shown to have strong effects on the plastic zone sizes and shapes which could affect the toughening of pressure-sensitive materials.     

Plane-strain mixed-mode near-tip fields in elastic perfectly plastic solids under small-scale yielding conditions

Within the context of the small-strain approach, plane-strain mixed-mode near-tip fields of a stationary crack in an elastic perfectly plastic Mises solid under small-scale yielding conditions are examined by finite element methods. Steady-state stress fields in the immediate vicinity of the crack tip are obtained as the remote loading of the elastic K-field increases. Asymptotic crack-tip solutions consisting of constant stress sectors, centered fan sectors, and an elastic sector are then constructed accordingly. The asymptotic crack-tip stress solutions agree well with the numerical results for a whole spectrum of mixed-mode loadings. Our mixed-mode near-tip solution with an elastic sector differs from that of Saka et al. by one (plastic) constant stress sector situated between the elastic sector and the neighbouring fan sector. The effect of the existence of the elastic sector on the near-tip fields is discussed in the light of the computational results. The plastic mixity factor of the near-tip field is given as a function of the elastic mixity factor of the prescribed K-field. This function is well bounded by that of the perfectly plastic limit of the corresponding solutions for power-law hardening materials given by Shih. Some issues pertaining to the numerical procedures such as the implementation of the small-scale yielding assumption are also addressed.     

Constraint effects on plastic crack-tip fields for plane strain mode-I, II and III cracks in non-hardening materials

To explore constraint effects on fully plastic crakc-tip fields, analytical solutions are examined for mode-I, II and III loading in non-hardening materials under plane strain conditions. The results reveal that under mode-II and III loading the crack-tip stress fields are unique, and thus can be characterized by a `single parameter'. Under mode-I loading, however, the crack-tip stress field is non-unique but can be characterized by two sets of solutions or `two parameters'. One set of the solutions is the well-known Prandtl field and the other is a plastic T-stress field. This conclusion corroborates the observation of McClintock (1971) that the slip-line field is non-unique for plane strain tensile cracks. A two-term plastic solution which combines the Prandtl field and the plastic T-stress field with two parameters B ₁ and B ₂ can then characterize the crack-tip stress field of plane strain mode-I crack over the plastic region and quantify the magnitude of crack-tip constraints. These characters are similar to those for hardening materials. Analyses and examples show that the two-term plastic solution can match well with the slip-line field or finite element results over plastic region. Thus the parameters B ₁ and B ₂ can be used to characterize the constraint level for mode-I finite-sized crack specimens in non-hardening materials under plane strain conditions.     

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.     

Elastoplastic crack analysis for pressure-sensitive dilatant materials-Part II: Interface cracks

The problem of a plane strain crack lying along an interface between a rigid substrate and an elastic-plastic material has been studied. The elastic-plastic material exhibits pressure-sensitive yielding and plastic volumetric deformation. Two-term expansions of the asymptotic solutions for both closed frictionless and open crack-tip models have been obtained. The Mises effective stress in the interfacial crack-tip fields is a decreasing function of the pressure-sensitivity in both open and closed-crack tip models. The variable-separable solution exists for most pressure-sensitive materials and the limit values for existence of the variable-separable solution vary with the strain-hardening exponents. The governing equations become singular as the pressure-sensitivity limit is approached. Strength of the leading stress singularity is equal 1/(n+1) for both crack-tip models, regardless of the pressure-sensitivity. The second-order fields have been solved as an additional eigenvalue problem and the elasticity terms do not enter the second-order solutions as n2. The second exponents for the closed crack model are negative for the weak strain hardening, whereas the negative second-order eigenvalue in the open crack model slightly grows with the pressure-sensitivity. The second-order solutions are of significance in characterising the crack-tip fields. The leading-order solution contains the dominant mode I components for both open and closed crack-tip models when the materials do not have substantial strain hardening. The second-order solutions are generally mode-mixed and depend significantly on the pressure-sensitivity.     

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.     

Mode I plane stress crack-tip constraint for elastic-plastic materials with different strain hardening

Finite element analyses and simulations have been undertaken to investigate the triaxial constraint in the crack-tip regions of a stationary crack and a steady-state growing crack under mode I plane stress for elastic-plastic materials with different strain hardening. The results show that the triaxial constraint in the crack-tip region is independent of specimen geometry, and material strain hardening, both for a stationary and an extending crack quasi-statically. The triaxial constraints for the various configurations examined are in better accordance with those required by the HRR solution for a stationary crack, which defines the low and similar constraints in crack-tip regions for different material strain hardening in the plane stress case. Along the entire ligament ahead of a crack tip, the constraint level transites gradually from that defined by the HRR solution within the near tip zone to that characterized by the stress intensity factor K _I in the far field.     

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.     

Characteristic crack-tip distances in fracture criteria: Is crack propagation discontinuous?

The present paper describes a possible mechanism for discontinuous crack advance in which surface separation occurs initially not at the crack-tip itself but within the crack-tip plastic zone of size r_p, at the mid-point of the crack-tip characteristic distance d (identified here with the finite growth step Δa), i.e., at the region of maximum opening tensile stress, spreading towards (and also away from) the crack-tip. The crack extension occurs when the crack-tip is reached and full opening over the distance d is completed.Finite element analyses show that this mechanism causes the formation of a rippled crack face surface in elastic-plastic materials in which irreversible plastic deformations take place during each growth step, in sharp contrast with the smooth surface created in ideal elastic materials in which all deformations are fully reversible. Some pictorial evidence of void formation ahead of the crack tip and of ripples during propagation, found in the literature, is presented.Although the present analysis is from a continuum standpoint it is acknowledged that micro structural features and mechanisms can condition the fracture events taking place in the process zone.The implication to the brittle-ductile transition of the dependence of the energy release rate, G^ΔΞ, on the ratio q (=Δa/r_p) is also discussed.     

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.     

Steady crack growth in a thin, ductile plate under small-scale yielding conditions: Three-dimensional modeling

This work employs high resolution, finite element computations to investigate key features of the elastic–plastic fields near a steadily advancing crack at quasi-static rates under three-dimensional, small-scale yielding conditions. The model represents a structurally thin component constructed of a material (e.g., Al and Ti alloys) with flow stress and fracture toughness properties that together limit the size of the in-plane plastic zone during steady-growth to no more than several multiples of the plate thickness. The computational approach generalizes the streamline integration procedure used previously for two-dimensional studies into three dimensions to represent steady-state growth on a fixed mesh in a boundary-layer framework. The plate thickness provides the only geometrical length scale. Crack extension occurs at the remotely applied, fixed loading without the need for a local growth criterion. In the first computations of this type, the present work considers a straight crack front advancing under local and global mode I loading with zero T-stress in a moderately hardening material. Applied remote loads at steady growth generate plastic zone sizes ahead of the advancing crack front ranging from 0.25 to 6.4 times the thickness. Key results include: (1) the crack-front fields exhibit a self-similar scaling characterized by a non-dimensional loading parameter; (2) three-dimensional effects extend to distances of approximately 1.5–2.5 times the thickness ahead of the advancing crack front for key values of this loading parameter, beyond which the fields (elastic–plastic then linear-elastic at greater distances) become uniform over the thickness; and (3) crack opening profiles on the outside surface reveal a “wedge-like”, opening shape which simplifies the definition of a crack-tip opening angle.     

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.     

Element free galerkin study of steady quasi-static crack growth in plane strain tension in elastic-plastic materials

Employing an extension of the convective mesh technique, an Element Free Galerkin (EFG) based formulation for steady quasi-static crack growth in elastic-plastic materials undergoing small scale yielding is established. In accordance with the steady state condition, a parallel path constitutive law integration scheme is implemented into the formulation and a body force type term is introduced for the contribution of crack tip plasticity. A detailed numerical analysis is performed for steady quasi-static growth of a mode I crack under plane strain conditions in elastic-perfectly plastic materials. The numerical solutions to this problem confirm the existence of an elastic unloading wedge. The computed near-tip stress distributions and the crack opening displacements for this problem provide good agreements with the corresponding asymptotic solutions and demonstrate the validity of the EFG based formulation. A rough estimate of the range of validity of the asymptotic solutions is also made based on the numerical solutions. In addition, the effect of strain hardening on the steady quasi-static crack growth is investigated. Both power hardening and linear hardening models are considered.     

A hybrid finite element analysis of interface cracks

A versatile hybrid finite element scheme consisting of special crack-tip elements and crack face contact elements is developed to analyse a partially closed interface crack between two dissimilar anisotropic elastic materials. The crack-tip element incorporates higher-order asymptotic solutions for an interfacial crack tip. These solutions are obtained from complex variable methods in Stroh formalism. For a closed interfacial crack tip, a generalized contact model in which the crack-tip oscillation is eliminated is adopted in the calculation. The hybrid finite element modelling allows the stress singularity at an open and closed crack tip to be accurately treated. The accuracy and convergence of the developed scheme are tested with respect to the known interface crack solutions. Utilizing this numerical scheme, the stress intensity factors and contact zone are calculated for a finite interface crack between a laminated composite material.     

Continuous near-tip fields for a dynamic crack propagating in a power-law elastic-plastic material

By use of the J ₂ flow theory and the rectangular components of field quantities, the near-tip asymptotic fields are studied for a dynamic mode-I crack propagating in an incompressible power-law elastic-plastic material under the plan strain conditions. Through assuming that the stress and strain components near a dynamic crack tip are of the same singularity, the present paper constructs with success the fully continuous dominant stress and strain fields. The angular variations of these fields are identical with those corresponding to the dynamic crack propagation in an elastic-perfectly plastic material (Leighton et al., 1987). The dynamic asymptotic field does not reduce to the quasi-static asymptotic field in the limit as the crack speed goes to zero. This revised version was published online in July 2006 with corrections to the Cover Date.     

Three-dimensional numerical simulations of dynamic fracture in silicon carbide reinforced aluminum

The three-point bending test by Kolsky-bar apparatus is a convenient technique to test the dynamic fracture properties of materials. This paper presents detailed three-dimensional finite element simulations of a silicon particle reinforced aluminum (SiC_p/Al) experiment (Li et al., [Proceedings of the US Army Symposium on Solid Mechanics]. In the simulations, the interaction between the input bar and the specimen is modeled by coupled boundary conditions. The material model includes large plastic deformations, strain-hardening and strain-rate hardening mechanisms. Furthermore, crack initiation and propagation processes are simulated by a cohesive element model. The simulation results quantitatively agree with the experimental measurements on three fronts: (1) the structural response of the specimen, (2) the time of unstable crack propagation, and (3) the local deformations at the crack-tip zone. The simulations reveal crack propagation characteristics, including crack-tip plastic deformation, crack front curving, and crack velocity profile. The effectiveness of Kolsky-bar type fracture tests is verified. It is shown that a rate-independent cohesive model can describe the complicated dynamic elastic–plastic fracture process in the SiC_p/Al material.     

Analytical solutions of fully plastic crack-tip higher order fields under antiplane shear

Analytical solutions of higher order fields in a fully plastic power-law hardening material are presented. By the use of hodograph transformation and asymptotic analysis the stress and strain exponents, angular distributions of shear stresses and strains are analytically determined. Special cases, such as linearly elastic, perfectly plastic materials are discussed. Similar characteristics between mode III and mode I plane strain, and mode II plane stress are examined. Comparison of four-term asymptotic solutions with exact and leading term solutions in an infinite strip with a semi-infinite crack under constant displacements along its edges is provided.     

Damage Field Ahead of a Tensile Crack in an Elastic-Plastic and Viscoplastic Material

This paper investigates the interaction between a macroscopic crack and microscopic damage in an elastic-plastic and viscoplastic material subjected to tensile in-plane loading. The aim is to predict the fracture conditions by accounting for void accumulation in the vicinity of the crack-tip. A power law relates the stress to the strain of the material. The damage, which invokes the growth and coalescence of microvoids, is confined to a small circular zone surrounding the crack-tip. At the onset of crack extension, the applied stress for small-scale and large-scale yielding solutions is found to be proportional to a₀ ^-1/(n+1), where 2a₀is the initial crack length and n is the strain hardening exponent of the material. For small-scale yielding, the conditions required for fatigue crack growth and steady-state creep are determined. In particular, the variations of the normalized crack length with the number of loading cycles and the time required for failure are shown for various strain hardening exponents, applied loading, and material damage parameters. This revised version was published online in July 2006 with corrections to the Cover Date.