Tensile cracks in polycrystalline ice under transient creep: Part II — Numerical simulations

This paper discusses predictions of a numerical model presented in the companion paper (Nanthikesan and Shyam Sunder, 1995) to analyze tensile cracks in polycrystalline ice undergoing transient creep. The numerical model is based on the internal state variable constitutive theory of transient creep in ice developed by Shyam Sunder and Wu (1989a,b, 1990). The finite element model uses the boundary layer approach of Rice (1968), in conjunction with a mid-point crack-tip element and reduced integration, to simulate the asymptotic stress and deformation fields in the vicinity of the crack tip, including incompressible creep deformations.

The problem of a stationary, traction-free, tensile (mode I) crack is analyzed to predict the size, shape and time evolution of the creep-dominated fracture process zone surrounding the crack-tip. The numerical simulations quantify the effects of transient creep, material strain hardening, fabric anisotropy, loading rate, temperature, and finite fracture test-specimen boundary on the development of the creep zone. A range of stress-intensity rates from 1 to 100 kPa s⁻¹ and temperatures from −5° to −25°C is considered in the simulations.

The results from a comprehensive numerical simulation study show that: (i) transient creep increases the creep zone size by more than an order of magnitude over that for a power-law creeping material, i.e., about 40 times for the isotropic, equiaxed granular ice tested by Jacka (1984); (ii) material strain hardening significantly affects the creep zone size, i.e., the creep zone for the transversely-isotropic columnar-grained ice tested by Sinha (1978), with the crack loaded in the plane of isotropy, is about 4 times smaller than that for the granular isotropic ice; (iii) fabric anisotropy increases the size of the creep zone by a factor of at least two for cracks in the transversely-isotropic, columnar-grained ice loaded in the plane of isotropy; (iv) the Riedel and Rice (1980) equation, which was derived for an isotropic power-law creeping material subjected to a suddenly applied constant stress-intensity, overestimates the creep zone size by a factor of 4.2 for a constant stress-intensity rate loading; and (v) as the crack size increases, linear elastic fracture mechanics becomes increasingly applicable at lower loading rates and higher temperatures.     

Crack-tip fields in elastic-plastic material under plane stress mode I loading

New results on the crack-tip fields in an elastic power-law hardening material under plane stress mode I loading are presented. Using a generalized asymptotic expansion of the stress function, higher-order terms are found which have newly-discovered characteristics. A series solution is obtained for the elastic-plastic crack-tip fields. The expansion of stress fields contains both the and terms where t_i is real and t_k is complex; the terms σ⁽ⁱ⁾ _pq(θt_i) and σ^(k) _rsθt_k) are real and complex functions of θ respectively. Comparing the results with that for the plane strain mode I loading shows that: (1) the effect of higher-order solutions on the crack-tip fields is much smaller; and (2) the path-independent integral J also controls the second-order or third-order term in the asymptotic solutions of the crack-tip fields for most of the engineering materials (1 < n < 11) in plane stress, while the J-integral does not control the second and the third-order terms for the plane strain mode I case for n > 3. These theoretical results imply that the crack-tip fields can be well characterized by the J-integral, and can be used as a criterion for fracture initiation under plane stress mode I loading. This is in agreement with existing full-field solutions and experimental data that J at crack growth initiation is essentially independent of in-plane specimen geometry. The comparison confirms the theoretical asymptotic solutions developed in this study. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

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.     

Strain-energy density and fracture-process zone. Report 1. Theoretical premises

A representation is proposed for the total strain-energy density (SED) in dimensionless form. The -distribution of the elastic and plastic parts of the total SED, and also the strain-energy-density factor S for the full range of mixed fracture modes under plane strain and a plane stress state are presented. Differences, as defined by the strain-energy-density factor, are established between fracture conditions under plane strain and a plane stress state. A general relationship between strain-energy density and a material parameter (the size of the fracture-process zone) is brought to light in dimensionless form. The relationship obtained predicts the pattern of continuous variation in the size of the fracture-process zone from static to cyclic material fracture. A material constant can be considered one of the consequences of fracture. Equations are proposed for calculation of the size of the fracture-process zone under a plane stress state and plane strain over the full range of mixed deformation modes with respect to standard mechanical properties of the materials.Translated from Problemy Prochnosti, No. 10, pp. 3–17, October, 1995.     

Complete theoretical analysis for higher order asymptotic terms and the HRR zone at a crack tip for Mode I and Mode II loading of a hardening material

Summary A complete development for the first two terms of the crack tip fields for both Mode I and Mode II loading of a hardening material in either plane stress or plane strain is performed, including the elastic deformation in the analysis. It is shown that the determination of the order of the second term depends on bothn and whether plane stress or plane strain is considered. In addition, regions of HRR dominance at a crack tip for the field variables are estimated. Comparison of the analytic predictions with finite element results indicates that the analytic results for the zone of HRR dominance are in agreement with numerical predictions.     

Tensile cracks in polycrystalline ice under transient creep: Part I — Numerical modeling

The interaction between creep deformations and a stationary or growing crack is a fundamental problem in ice mechanics. Knowledge concerning the physical mechanisms governing this interaction is necessary: (1) to establish the conditions under which linear elastic fracture mechanics can be applied in problems ranging from ice-structure interaction to fracture toughness testing; and (2) to predict the ductile-to-brittle transition in the mechanical behavior of ice and, especially, the stability and growth of cracks subjected to crack-tip blunting by creep deformations. This requires a quantitative estimate of the creep zone surrounding a crack-tip, i.e., the zone within which creep strains are greater than the elastic strains.

The prediction of the creep zone in previous ice mechanics studies is based on the theory developed by Riedel and Rice (1980) for tensile cracks in creeping solids. This theory is valid for a stationary crack embedded in an isotropic material obeying an elastic, power-law creep model of deformation and for a suddenly applied uniform far-field tension load that is held constant with time. The deformation of ice at strain-rates ahead of a crack (i.e., 10⁻⁶ to 10⁻² s⁻¹) is dominated, however, by transient (not steady power-law) creep and the loading, in general, is not instantaneous and constant.

A numerical model is developed in this paper to investigate the role of transient creep and related physical mechanisms in predicting the size, shape and time evolution of the creep zone surrounding the tip of a static crack in polycrystalline ice. The model is based on the fully consistent tangent formulation derived in closed form (Shyam Sunder et al., 1993) and used in the solution of the physically-based constitutive theory developed by Shyam Sunder and Wu (1989a, b) for the multiaxial behavior of ice undergoing transient creep. The boundary value problem involving incompressible deformations ahead of a stationary, traction-free mode I crack in a semi-infinite medium is modeled and solved by a finite element analysis using the boundary layer approach of Rice (1968). This model is verified by comparing its predictions with (i) the known theoretical solutions for the elastic and HRR asymptotic stress and strain fields in an elastic-plastic material of the Ramberg-Osgood type, and (ii) the creep zone size for an isotropic material obeying the elastic power-law creep model of deformation.     

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.     

Transient crack-tip fields for mixed-mode power law creep

Mixed-mode stationary crack-tip fields are obtained for an elastic-nonlinear viscous power law creeping solid under conditions of plane strain and small-scale creep. Power law exponents of 2 and 5 are considered which are representative of the creep response of a wide range of ceramics and metals. The imposed far-field mixity ranges from pure mode I to pure mode II. Crack tip fields are calculated during the transient regime using a detailed finite element analysis and are shown to be governed by a Hutchinson-Rice-Rosengren type singularity over the inner one fifth of the creep zone. Dominance of universal mixed-mode near-tip fields within the inner creep zone is found for several mixtures of far-field mode I and mode II. The pronounced effects of the amount of mixity on the size and shape of the creep zone as well as on the time required to reach extensive creep conditions are determined. For a creep exponent of 5, it is estimated that the creep zone grows about seven times faster in mode II than in mode I, with a corresponding decrease in the transition time from small-scale to extensive creep. For a creep exponent of 2, the creep zone grows about six times faster in mode II than in mode I. Finally, the mixed-mode creep fields are used to assess possible beneficial effects of crack deflection or branching in metals and ceramics at elevated temperatures.     

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.     

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.     

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 tensile crack in creeping solids with large damage near the crack tip

An asymptotic analysis of stationary mode I crack in creeping solids with large damage near crack tip is conducted. To consider the damage effect, Kachanov damage evolution law is utilized and incorporated into the power-law creep constitutive equation. With the compatibility equation, a nonlinear eigenvalue problem which can be solved by numerical approaches is established. From this result, the distribution of stress and strain rate are obtained with the coupling effect of damage and creep under plane stress condition. Also the influence of material parameters on the stress is examined. According to the result, it is shown that the creep exponent n and damage parameter (=/(1+k)) have a significant effect on determining the eigenvalue s and angular distribution of stress and strain rate near the crack tip. The creep exponent n plays the role to soften and damage parameter plays the role to harden the material near the crack tip. The stress and strain rate show quite different behavior from those of HRR problem.     

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.     

A Three-dimensional Numerical Study of Mixed Mode (I and II) Crack Tip Fields in Elastic–plastic Solids

The stress fields near a crack front in a ductile solid are essentially three-dimensional (3D) in nature. The objective of this paper is to investigate the structure of these fields and to establish the validity of two-dimensional (2D) plane stress and plane strain approximations near the crack front under mixed mode (combined modes I and II) loading. To this end, detailed 3D and 2D small strain, elastic–plastic finite element simulations are carried out using a boundary layer (small scale yielding) formulation. The plastic zones and radial, angular and thickness variations of the stresses are studied corresponding to different levels of remote elastic mode mixity and applied load, as measured by the plastic zone size with respect to the plate thickness. The 3D results are compared with those obtained from 2D simulations and asymptotic solutions. It is found that, in general, plane stress conditions prevail at a distance from the crack front exceeding half the plate thickness, although it could be slightly smaller for mode II predominant loading. The implications of the 3D stress distribution on micro-void growth near the crack front are briefly discussed.     

Creep crack growth in an elastic-creeping material Part II: mode I

The quasi-static growth of a crack in an elastic-creeping material under mode I loading is investigated. The creep strain rate of the material is assumed to be governed by a power law involving the stress and creep strain. The major emphasis of this investigation is on elastic-primary creep response. The asymptotic crack tip fields for a quasi-statically extending crack under conditions of plane strain and plane stress are developed. The asymptotic fields are unambiguously determined in terms of the instantaneous crack speed and material parameters and are independent of the prior crack history, specimen geometry, and loading. A plane strain finite element analysis is performed to determine the complete stress and strain fields. These fields are compared with the asymptotic ones to establish the zone of dominance of the crack tip fields. The zone of dominance can be a very small fraction of the size of the creep zone attending the crack up.

---

Résumé On étudie la croissance quasi-statique d'une fissure dans un matériau sujet à fluage élastique sous une sollicitation de Mode I. On suppose que la vitesse de déformation en fluage du matériau est régie par une loi parabolique comportant la contrainte et la déformation. L'accent est surtout mis sur la réponse au fluage élastique primaire.On établit les configurations asymptotiques règnant à l'extrémité de la fissure, dans le cas d'une fissure quasi-statique en extension sour état plan de dilatation ou sous état plan de tension. Ces champs sont déterminés sans ambiguïté par la vitesse instantanée de la fissuration et par les paramétres du matériau; ils sont indépendants de l'histoire primitive de la fissure, de la géométrie de l'éprouvette et de la sollicitation. On effectue une analyse par éléments finis en déformations planes pour déterminer complètement les champs de contrainte et de dilatation, que l'on compare aux champs asymptotiques aux fins d'établir la zone de prédominance des champs à l'extrémité de la fissure. Cette zone peut être une fraction très petite de la taille de la zone où s'effectue un fluage au voisinage de l'extrémité de la fissure.

     

A criterion for plane strain,fully plastic,quasi-steady crack growth

Plane strain fracture by hole growth in ordinary-sized parts of low-to-medium strength steels is essentially rigid-plastic, and may be approximated as non-hardening. Quasi-steady crack growth for such materials is predicted for crack-tip fields approximated by a pair of slip lines, such as unequally grooved specimens in tension and deep singly-face-cracked specimens under combined bending and tension. The crack growth increment a is given in terms of material parameters, far-field geometry, and loadings and their increments.For the rigid-plastic, non-hardening approximation, stress and strain increment fields for growing cracks are identical to those for stationary cracks. For fields with a pair of symmetric slip-lines, the flanks of the decohering zone turn out to be rigid, and the decohering zone does not affect the crack-tip opening angle (CTOA), which then depends only on the micromechanisms of hole nucleation, growth and linkage by flow localization or fine cracking. These mechanisms are in turn approximately controlled by the near-field plasticity parameters: the angle of the slip plane _s, and the normal stress and displacement increment across the slip plane _s and u_s. Note the three-parameter characterization of the near-tip fields, in contrast to the one- or two-parameter characterization in elastic or nonlinear elastic fracture mechanics.A sliding off and shear-cracking model for a growing crack, based on a hole growth equation, gives an approximate CTOA in terms of _s, _s, and material parameters. When hole nucleation strain is negligible, the estimated CTOA exhibits an inverse exponential dependence on _s and a higher order parabolic dependence on _s. For a given material, a series of fully plastic crack growth experiments is suggested to determine the approximate material parameters needed to characterize the dependence of CTOA on _s and _s, or from kinematics, of the shear strain behind the slip plane, _f, on _s.     

Effect of Plate Thickness on Crack-Tip Plasticity

This paper presents an analytical method for determining the three-dimensional stress fields in plates with a through-the-thickness crack, especially under elastic-plastic conditions. Using the generalised plane strain theory in conjunction with the deformation theories of plasticity, exact solutions are obtained for the effects of plate thickness on the crack-tip plastic zone size and a plastic constraint factor, which is shown to correlate well with published finite element solutions.     

Creep crack growth in an elastic-creeping material Part I: mode III

The quasi-static growth of a crack in an elastic-creeping material under antiplane shear or mode III loading is investigated. The creep response of the material is assumed to be governed by a power-law between the creep strain rate, creep strain, and stress. While this law is capable of describing elastic-primary, secondary, or tertiary creep, the major emphasis of this paper is on crack growth in the elastic-primary creep regime. The asymptotic crack-tip stress and strain fields for a quasi-statically extending crack in an elastic-primary creeping material are developed. This is followed by a finite element analysis to determine the complete stress and strain fields within the confines of small scale yielding. These fields are then compared with the asymptotic ones to establish the size of the zone of dominance of the crack-tip fields.

---

Résumé On étudie la croissance quasi-statique d'une fissure dans un matériau en phase de fluage élastique sous des cisaillements anti-planaires ou sous une sollicitation de mode III.La réponse du matériau sur fluage est supposée gouvernée par une loi parabolique qui lie la vitesse de déformation par fluage, la déformation par fluage et la contrainte. Si cette loi est à même de décrire les fluages élastiques-primaire, secondaire et tertiaire-l'accent est ici principalement placé sur la croissance d'une fissure sous un fluage élastique-primaire. On établit les champs de contrainte et de déformation asymptotiques à l'extrémité de la fissure dans le cas d'une fissure en croissance quasi-statique, dans un matériau soumis à fluage élastique primaire.On procède ensuite à une analyse par éléments finis afin de déterminer complètement les champs de contraintes et de déformation aux confins d'une zone de petite taille en déformation plastique. On compare ces champs aux champs asymptotiques, en vue d'établir la taille de la zone où les champs à l'extrémité de la fissure s'exercent de manière dominante.

     

Creep fracture in Al-SiC metal-matrix composites

Creep fracture behaviour of pure aluminium-matrix composites with 10–30 vol% SiC particulates at 623 K is reported. A comparison of tensile and compression creep data shows the existence of a transition stress. Above this transition stress no steady state creep is observed in tension. This transition stress is related to a transition from intergranular to transgranular fracture. The origin of transition stress is perhaps associated with the diffusional relaxation of stress concentration at the matrix/particle interface by lattice diffusion. The intergranular creep fracture of composites appears to be similar to that of unreinforced aluminium and it is power-law creep controlled. The transgranular creep fracture occurs by void nucleation and growth. The nucleation strain for voids is quite small and hence the tertiary stage starts before the end of the primary stage. The ductile fracture models overestimate the strain to fracture and do not predict the observed stress dependence of strain to fracture.