Quasi-statically growing crack-tip fields in elastic perfectly plastic pressure-sensitive materials under plane strain conditions

Quasi-statically growing crack-tip fields in elastic perfectly plastic pressure-sensitive materials under plane strain conditions are investigated in this paper. The materials are assumed to follow the Drucker-Prager yield criterion and the normality flow rule. The asymptotic mode I crack-tip fields are assumed to follow the five-sector assembly of Drugan et al. (1982) for Mises materials. The crack-tip sectors, in turns, from the front of the crack tip are a constant stress sector, a centered fan sector, a non-singular plastic sector, an elastic sector and finally a trailing non-singular plastic sector bordering the crack face. The results of the asymptotic analysis show that as the pressure sensitivity increases, the plastic deformation shifts to the front of the tip, the angular span of the elastic unloading sector increases, and the angular span of the trailing non-singular plastic sector bordering the crack surface decreases. As the pressure sensitivity increases to about 0.6, the angular span of the trailing non-singular plastic sector almost vanishes. The effects of the border conditions between the centered fan sector and the first non-singular plastic sector on the solutions of the crack-tip fields for both Mises and pressure-sensitive materials are investigated in details. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

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.     

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.     

Asymptotic fields near the crack tip in elastic-perfectly plastic crystals

The basic equations of plane strain problem for the elastic-perfectly plastic crystals with double slip systems have been presented in the basis of three dimensional flow theory of crystal plasticity. Using these equations the stationary crack tip stress and deformation fields are analysed for tensile load. The fields involve an elastic angular sector and are fully continuous. An asymptotic solution is also obtained for the steadily growing crack that consists of five angular sectors: two plastic angular sectors in the front of the crack tip connected with the boundary on which the associated velocity field has discontinuities; a secondary plastic angular sector near the crack face; two elastically unload angular sectors connected with the boundary on which the discontinuity of the associated velocity field occurs. The asymptotoic solution is not unique. A family of solutions is obtained. Finally, the application of these solutions on both FCC and BCC crystals is discussed.     

Plane-stress crack-tip fields for perfectly plastic orthotropic materials

Plane stress mode I crack-tip fields for perfectly plastic orthotropic materials are studied. Plastic orthotropy is described by Hill's quadratic yield function. The construction of crack-tip fields is based on the general crack-tip field analysis for elastic perfectly plastic materials given by Rice [1] and guided by the corresponding low-hardening power-law solutions. Two very different types of plane-stress crack-tip fields emerge as plastic orthotropy is varied. The first one consists of a centered fan sector in front of the crack tip and two neighboring constant stress sectors. The second one consists of a constant stress sector in front of the crack tip, a constant stress sector bordering the crack face, and a centered fan sector between the two constant stress sectors. All the perfectly plastic crack-tip solutons are verified by the corresponding low-hardening power-law solutions. General trends of crack-tip stress solutions as functions of plastic orthotropy and implications of these solutions to the design of ductile composite materials are discussed.

---

Résumé On étudie les champs de contraintes planes de mode I à l'extrémité d'une fissure, dans les matériaux orthotropiques parfaitement plastique. L'orthotropie plastique est décrite par la fonction quadratique de plasticité de Hill. On base les constructions des champs de constraintes sur l'analyse générale des constraintes à l'extrémité d'une fissure fournie par Rice pour les matériaux élastiques parfaitement plastiques, que l'on règle par les lois paraboliques caractérisant un faible écrouissage. Lorsque l'on modifie l'orthotropie plastique, il apparaît deux types de champs de contraintes à l'extrémité de la fissure très différents. Le premier comporte un secteur en éventail centré sur le front de fissure, et deux secteurs voisins à contraintes constantes. Le second consiste en une secteur à contrainte au bord de la surface de la fissure, et un secteur en éventail centré sur les deux secteurs à contraintes constantes. Toutes les solutions relatives à une extrémité de fissures parfaitement plastique sont vérifiées par les fonctions paraboliques d'écrouissage faible correspondantes. On discute des tendances générales que suivent les solutions pour les contraintes en extrémité de fissure selon l'orthrotropie plastique, et des implications que comportent ces solutions dans la conception de matériaux composites ductiles.

     

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.     

Asymptotic analysis of growing crack stress/deformation fields in porous ductile metals and implications for stable crack growth

Asymptotic stress and deformation fields near a quasi-statically growing plane strain tensile crack tip in porous elastic-ideally plastic material, characterized by the Gurson-Tvergaard yield condition and associated flow rule, are derived for small uniform porosity levels throughout the range 0 to 4.54 percent. The solution configuration resembles that for crack growth in fully dense, elastically compressible, elastic-ideally plastic Huber-Mises material for this porosity range, except that the angular extents and border locations of near-tip solution sectors vary with porosity level, as do the stress and deformation fields within sectors. Increasing porosity is found to result in a dramatic reduction in maximum hydrostatic stress level, greater than that for a stationary crack; it also causes a significant angular redistribution of stresses, particularly for a range of angles ahead of the crack and adjacent to the crack flank. The near-tip deformation fields derived are employed to generalize a previously-developed, successful ductile crack growth criterion. Our model predicts that for materials having the same initial slopes of their crack growth resistance curves, but different levels of uniform porosity, higher porosity results in a substantially greater propensity for stable crack growth.     

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

The asymptotic structures of crack-tip stress and deformation fields are investigated numerically for quasi-static and dynamic crack growth in isotropic linear hardening elastic-plastic solids under mode I, plane stress, and small-scale yielding conditions. An Eulerian type finite element scheme is employed. The materials are assumed to obey the von Mises yield criterion and the associated flow rule. The ratio of the crack-tip plastic zone size to that of the element nearest to the crack tip is of the order of 1.6 × 10⁴. The results of this study strongly suggest the existence of crack-tip stress and strain singularities of the type r ^s (s < 0) at r=0, where r is the distance to the crack tip, which confirms the asymptotic solutions of variable-separable type by Amazigo and Hutchinson [1] and Ponte Castañeda [2] for quasi-static crack growth, and by Achenbach, Kanninen and Popelar [3] for dynamic crack propagation. Both the values of the parameter s and the angular stress and velocity field variations from the present full-field finite element analysis agree very well with those from the analytical solutions. It is found that the dominance zone of the r ^s-singularity is quite large compared to the size of the crack-tip active plastic zone. The effects of hardening and inertia on the crack-tip fields as well as on the shape and size of the crack-tip active plastic zone are also studied in detail. It is discovered that as the level of hardening decreases and the crack propagation speed increases, a secondary yield zone emerges along the crack flank, and kinks in stress and velocity angular variations begin to develop. This dynamic phenomenon observable only for rapid crack growth and for low hardening materials may explain the numerical difficulties, in obtaining solutions for such cases, encountered by Achenbach et al. who, in their asymptotic analysis, neglected the existence of reverse yielding zones along the crack surfaces.     

An Experimental Study on the Elastic–Plastic Fracture in a Ductile Material under Mixed-Mode Loading

Abstract:  Stereo vision was used to measure the crack-tip parameters, such as J integral, plastic mixity and elastic mixity of mixed-mode fracture specimens, and to study the applicability of the Shih's plane strain solution to the mixed-mode crack-tip fields. The fracture specimen used in this study was a compact tension shear (CTS) specimen made of 2024-O aluminum. The in-plane strain and stress fields near the mixed-mode crack tip of the CTS specimen were determined using the deformation field measured by the stereo vision. It is observed that the J integral values computed along rectangular contours surrounding the mixed-mode crack-tip approach constant values after r / h  > 0.5. The in-plane strains determined experimentally at several points near the crack tip and at several radial lines emerging from the crack tip are compared with the values calculated using Shih's plane-strain solution and the HRR slope, named after the investigations of Hutchinson, Rice and Rosengren respectively. It is found that the measured values follow the trends of the Shih's plane-strain solution. The elastic mixity evaluated using the measured crack-tip stress fields is close to that obtained from analytical solution. However, the evaluated plastic mixity deviates from the analytical solution.     

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.     

Non-singular stress and velocity fields for crack-tip plastic zones and conditions for uniqueness

Non-singular plastic stress and velocity fields, for the tip of a crack of finite thickness and root radius, are developed as an elastic-plastic crack model that is likely to be more physically realistic than the classical infinitesimal crack with a plastic crack-tip singularity. With a non-singular plastic zone the velocity-field equations are not uniquely determined by the boundary conditions, under large geometrical changes, and they must therefore have the form of a wide set of kinematically-admissible velocity fields. These virtual velocity fields are used to establish the critical work-hardening rate to give a sufficient condition for uniqueness of the crack-tip velocity field in elastic-plastic fracture; it is shown that proof of uniqueness of the velocity field is likely to be an essential requirement for the valid application of elastic-plastic fracture mechanics.The elastic infinitesimal-crack model is shown to give an inadequate representation of the circumferential T-stress distribution at the surface of a crack of finite root radius, and this requires the adoption of a finite-thickness elliptical crack model to give approximate consistency between the elastic stress field and the non-singular plastic stress field at the crack tip.     

Measuring overload effects during fatigue crack growth in bainitic steel by synchrotron X-ray diffraction

In this work we present the results of in situ synchrotron X-ray diffraction measurements of fatigue crack-tip strain fields following a 100% overload (OL) under plane strain conditions. The study is made on a bainitic steel with a high toughness and fine microstructure. This allowed a very high (60 μm) spatial resolution to be achieved so that fine-scale changes occurring around the crack-tip were captured along the crack plane at the mid-thickness of the specimen. We have followed the crack as it grew through the plastic/residually stressed zone associated with the OL crack location. We observed two effects; one when the enhanced plastic zone is ahead of the crack and one after it has been passed. Regarding the former it was found that the compressive stress at the crack-tip initially falls sharply, presumably due to the increased plastic stretch caused by the OL. This is associated with a concomitant fall in peak tensile stress at K_max, the elastic excursion between K_min and K_max remaining essentially unchanged from before OL. Subsequently discontinuous closure as seen previously for plane stress caused by crack face contact at the OL location limits the elastic strain range experienced by the crack tip and thereby retards crack growth.     

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.     

Effects of plastic anisotropy on crack-tip behaviour

For a crack in a homogeneous material the effect of plastic anisotropy on crack-tip blunting and on the near-tip stress and strain fields is analyzed numerically. The full finite strain analyses are carried out for plane strain under small scale yielding conditions, with purely symmetric mode I loading remote from the crack-tip. In cases where the principal axes of the anisotropy are inclined to the plane of the crack it is found that the plastic zones as well as the stress and strain fields just around the blunted tip of the crack become non-symmetric. In these cases the peak strain on the blunted tip occurs off the center line of the crack, thus indicating that the crack may want to grow in a different direction. When the anisotropic axes are parallel to the crack symmetry is retained, but the plastic zones and the near-tip fields still differ from those predicted by standard isotropic plasticity.     

Fully plastic crack-tip fields for plane strain shallow-cracked specimens under pure bending

Detailed finite element (PE) analyses are performed to study the effect of crack depth on crack-tip constraint at full yielding for pure bending of plane strain single-edge-cracked specimens. Analyses are based on small-strain formulations and perfect plasticity. The crack depth a/W ranges from 0.1 to 0.7, and the deformation is applied up to the limiting state of full plasticity where crack-tip stresses reach steady-state limiting values.At load levels smaller than the limit load (contained yielding), the crack-tip constraint (stress triaxiality) gradually decreases as a/W decreases, but, at load levels close to the limit load (or at the limit load), it decreases very sharply. In terms of a/W, tractable closed-form approximations for fully plastic crack-tip stress and strain fields are proposed, and fully plastic values of crack-tip stresses are re-phrased in terms of the Q-parameter [1, 2]. The role of crack-tip strains on fracture of shallow-cracked bending specimens is briefly discussed.     

Analytical solutions for crack-tip sectors in perfectly plastic Mises materials under mixed in-plane and out-of-plane shear loading conditions

In this paper, analytical solutions for asymptotic crack-tip plastic sectors in perfectly plastic Mises materials are derived under mixed in-plane and out-of-plane shear loading conditions. Plastic strains in crack-tip plastic sectors are considered to be singular and non-singular. Sectors with singular plastic strains have the solution of centered fan type, and sectors with non-singular plastic strains have the solution of either centered fan or constant stress type. The requirement of stress continuity along the border between a constant stress and a centered fan sectors is then discussed. Discontinuities of the normal and out-of-plane shear stresses in the radial direction between two constant stress sectors are assumed in assembling the crack-tip fields under mixed mode II/III and I/III conditions. Crack-tip fields under mixed mode II/III and I/III conditions with small contributions of mode III are then presented to show the existence of asymptotic crack-tip fields for perfectly plastic materials under mixed in-plane and out-of-plane shear loading conditions. The trends of the angular variations of the mode III stresses under the mixed mode II/III and I/III conditions are generally in agreement with those of the available asymptotic and finite element analyses for low strain hardening materials.     

Effect of residual stresses on interface crack growth by void expansion mechanism

Crack growth along an interface between two adjacent elastic–plastic materials in a layered solid is analysed, using special interface elements to represent the fracture process ahead of the crack-tip. These interface elements account for ductile failure by the nucleation and growth of voids to coalescence. In these elements the stress components normal to the interface and the shear stresses are given by equilibrium with the surrounding material, and the stress component tangential to the interface is determined by the requirement of compatibility with the surrounding material in the tangential direction. It is assumed that the layers are sufficiently thick, so that the plastic regions around the crack-tip are much smaller than the thickness of the nearest layers. The analyses focus on the effect of initial residual stresses in the layered material, or on T-stress components induced during loading. The results show that the value of the T-stress component in the softer material adjacent to the interface crack plays the dominant role, such that a negative value of this stress component gives a significant increase of the interface fracture toughness.     

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.     

Modeling the brittle–ductile transition in ferritic steels: dislocation simulations

We present a model for the brittle–ductile transition in ferritic steels based on two dimensional discrete dislocation simulations of crack-tip plasticity. The sum of elastic fields of the crack and the emitted dislocations defines an elasto–plastic crack field. Effects of crack-tip blunting of the macrocrack are included in the simulations. The plastic zone characteristics are found to be in agreement with continuum models, with the added advantage that the hardening behavior comes out naturally in our model. The present model is composed of a macrocrack with microcracks ahead of it in its crack-plane. These microcracks represent potential fracture sites at internal inhomogeneities, such as brittle precipitates. Dislocations that are emitted from the crack-tip account for plasticity. When the tensile stress along the crack plane attains a critical value σ _F over a distance fracture is assumed to take place. The brittle–ductile transition curve is obtained by determining the fracture toughness at various temperatures. Factors that contribute to the sharp upturn in fracture toughness with increasing temperature are found to be: the increase in dislocations mobility, and the decrease in tensile stress ahead of the macrocrack tip due to increase in blunting, and the slight increase in fracture stress of microcracks due to increase in plasticity at the microcrack. The model not only predicts the sharp increase in fracture toughness near the brittle–ductile transition temperature but also predicts the limiting temperature above which valid fracture toughness values cannot be estimated; which should correspond to the ductile regime. The obtained results are in reasonable agreement when compared with the existing experimental data.