Development of a multiaxial fatigue theory by considering constraint effects on small mixed-mode cracks

The effects of the transverse strain (the normal strain in the crack-line direction) on the near-tip fields of small shallow surface cracks (Case A cracks) in power-law hardening materials are investigated by finite element analyses. The small Case A cracks are under plane stress, general yielding, and mixed mode I and II conditions. Constant effective stress contours representing the intense straining zones near the tip, deformed crack-tip profiles and near-tip mode mixity factors are presented for different transverse strains in the crack-line direction. Based on the concept of characterization of fatigue crack growth by the cyclic J-integral, the effects of the transverse strain on J are investigated. The results suggest that the fatigue life prediction based on multiaxial fatigue theories and the critical plane approach should include the constraint effects due to the transverse strain. Consequently, the concept of constant fatigue life contour on the Γ-plane in multiaxial fatigue theories is generalized to the constant fatigue life surface in the Γ-space where the shear strain and the two normal strains are the three axes. Finally, a damage parameter as a function of the shear strain and the two normal strains is proposed for evaluation of fatigue damage under multiaxial loading conditions.     

Finite deformation analysis of crack tip fields in plastically compressible hardening–softening–hardening solids

Crack tip fields are calculated under plane strain small scale yielding conditions. The material is characterized by a finite strain elastic–viscoplastic constitutive relation with various hardening–softening–hardening hardness functions. Both plastically compressible and plastically incompressible solids are considered. Displacements corresponding to the isotropic linear elastic mode I crack field are prescribed on a remote boundary. The initial crack is taken to be a semi-circular notch and symmetry about the crack plane is imposed. Plastic compressibility is found to give an increased crack opening displacement for a given value of the applied loading. The plastic zone size and shape are found to depend on the plastic compressibility, but not much on whether material softening occurs near the crack tip.On the other hand, the near crack tip stress and deformation fields depend sensitively on whether or not material softening occurs. The combination of plastic compressibility and softening(or softening–hardening) has a particularly strong effect on the near crack tip stress and deformation fields.     

Near-tip temperature distribution for mode III stationary crack in power-hardening materials

In this paper the near-tip temperature distribution for Mode III stationary crack in power-hardening materials is studied within the framework of thermodynamics. The hardening of materials is characterized as internal state variables to investigate the plastic dissipation. Numerical results indicate that the local temperature rise near crack tip by taking account of hardening is appreciably lower than those in perfectly-plastic case, and thus, are in better accord with the experimental results.The project is supported by National Natural Science Foundation of China     

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.     

ASYMPTOTIC ANALYSIS OF PLANE-STRAIN MODE I STEADY-STATE CRACK GROWTH IN TRANSFORMATION TOUGHENINGCERAMICS(Ⅱ)

I.IntroductionThemechanismoffracturetoughnessenhancementinceramicshasbeenwidelystudiedsincetheearlyl98Os.ThepioneeringconstitutivemodeldevelopedbyBudianskyetal.ll1includestheeffectofplasticdiIation,butneglectsthetransformation-inducedshearstrain.Andtheirc…     

Asymptotic fields for dynamic crack growth in non-associative pressure sensitive materials

Asymptotic near-tip fields are analyzed for a plane strain Mode I crack propagating dynamically in non-associative elastic–plastic solids of the Drucker–Prager type with an isotropic linear strain hardening response. Eigen solutions are obtained over a range of material parameters and crack speeds, based on the assumption that asymptotic solutions are variable-separable and fully continuous. A limiting speed, beyond which a tendency to slope discontinuity in angular distributions of stresses and velocities is detected, is found to deviate from the associative models. At low strain-hardening rates, the onset of the plastic potential corner zone ahead of the crack-tip imposes another limit to the crack speed. Correspondingly, those limits imply the limits to the degree of non-associativity at a given crack speed. In addition, a tendency to slope discontinuity in the angular radial stress distribution sets another limit on the non-associativity at vanishing hardening rates.     

Finite deformation analysis of crack-tip opening in elastic-plastic materials and implications for fracture

Analyses of the stress and strain fields around smoothly-blunting crack tips in both non-hardening and hardening elastic-plastic materials, under contained plane-strain yielding and subject to mode I opening loads, have been carried out by use of a finite element method suitably formulated to admit large geometry changes. The results include the crack-tip shape and near-tip deformation field, and the crack-tip opening displacement has been related to a parameter of the applied load, the J-integral. The hydrostatic stresses near the crack tip are limited due to the lack of constraint on the blunted tip, limiting achievable stress levels except in a very small region around the crack tip in power-law hardening materials. The J-integral is found to be path-independent except very close to the crack tip in the region affected by the blunted tip. Models for fracture are discussed in the light of these results including one based on the growth of voids. The rate of void-growth near the tip in hardening materials seems to be little different from the rate in non-hardening ones when measured in terms of crack-tip opening displacement, which leads to a prediction of higher toughness in hardening materials. It is suggested that improvement of this model would follow from better understanding of void-void and void-crack coalescence and void nucleation, and some criteria and models for these effects are discussed. The implications of the finite element results for fracture criteria based on critical stress or strain, or both, is discussed with respect to transition of fracture mode and the angle of initial crack-growth. Localization of flow is discussed as a possible fracture model and as a model for void-crack coalescence.     

ASYMPTOTIC ANALYSIS OF PLANE-STRAIN MODE I STEADY-STATE CRACK GROWTH IN TRANSFORMATION TOUGHENING CERAMICS(Ⅱ)

Based on a constitutive law which includes the shear components of transformation plasticity, the asymptotic solutions to near-tip fields of plane-strain mode I steadity propagating cracks in transformed ceramics are obtained for the case of linear isotropic hardening. The stress singularity, the distributions of stresses and velocities at the crack tip are determined for various material parameters. The factors influencing the near-tip fields are discussed in detail.Project supported by the National Natural Science Foundation of China     

On neck propagation in amorphous glassy polymers under plane strain tension

Unlike metals, necking in polymers under tension does not lead to further localization of deformation, but to propagation of the neck along the specimen. Finite element analysis is used to numerically study necking and neck propagation in amorphous glassy polymers under plane strain tension during large strain plastic flow. The constitutive model used in the analyses features strain-rate, pressure, and temperature dependent yield, softening immediately after yield and subsequent orientational hardening with further plastic deformation. The latter is associated with distortion of the underlying molecular network structure of the material, and is modelled here by adopting a recently proposed network theory developed for rubber elasticity. Previous studies of necking instabilities have almost invariably employed idealized prismatic specimens; here, we explicitly account for the unavoidable grip sections of test specimens. The effects of initial imperfections, strain softening, orientation hardening, strain-rate as well as of specimen geometry and boundary conditions are discussed. The physical mechanisms for necking and neck propagation, in terms of our constitutive model, are discussed on the basis of a detailed parameter study.     

Phenomenological nonlinear modelling of glassy polymers

Compared with the numerous works into the constitutive equations for the mechanical behaviour of metals, very little attention has been devoted to those of polymers. However, a model is required to describe both the complex shape of the stress–strain curves and strain rate sensitivity of glassy polymers. In this Note, a unified viscoelastic-viscoplastic model is presented in which the nonlinear pre-yield behaviour, the strain softening and the strain hardening are described by internal variables, in analogy with the models developed for metals. In order to check the predictive capability of the model, the numerical results are compared with the experimental data (monotone, creep and relaxation tests) of a typical amorphous glassy polymer. To cite this article: F. Zaïri et al., C. R. Mecanique 333 (2005).     

Void interaction and coalescence in polymeric materials

Beyond pressure-sensitivity, plastic deformation of glassy polymers exhibits intrinsic softening followed by progressive rehardening at large strains. This highly nonlinear stress–strain behavior is captured by a constitutive model introduced in this work. In the first part of the paper, we focus on void growth and coalescence in an axisymmetric representative material volume consisting of a single large void and a population of discrete microvoids. Our study shows that microvoid cavitation, enhanced by strain softening, accelerates the process of void coalescence resulting in brittle-like failure at lowered stresses and strains. Pressure-sensitivity also reduces stress-carrying capacity as well as influences the strain for void coalescence; plastic dilatancy effects are relatively milder. In the second part of the paper, we introduce a population of discrete spherical voids within a three-dimensional computational model to study void growth and damage ahead of a crack front. Our studies reveal a distinctive change in the deformed void shape from oblate to prolate when strain softening is followed by high rehardening at large plastic strains. By contrast, an extended strain softening regime promotes oblacity and facilitates multiple void interaction and their cooperative growth over large distances ahead of the crack front. This multi-void failure mechanism is exacerbated by pressure-sensitivity.     

Nonlinear response of anti-plane shear crack

Investigated in this communication is the nonlinear response of an anti-plane shear crack. A nonlinear hardening/softening relation is assumed between the resultant shear stress and shear strain for determining the state of affairs near the crack tip. Amplification of the local stress fields is obtained in terms of the material parameters by assuming that nonlinearity is confined to the immediate vicinity of the crack tip region.     

Steady crack growth in elastic–plastic fluid-saturated porous media

An asymptotic solution is obtained for stress and pore pressure fields near the tip of a crack steadily propagating in an elastic–plastic fluid-saturated porous material displaying linear isotropic hardening. Quasi-static crack growth is considered under plane strain and Mode I loading conditions. In particular, the effective stress is assumed to obey the Drucker–Prager yield condition with associative or non-associative flow-rule and linear isotropic hardening is adopted. Both permeable and impermeable crack faces are considered. As for the problem of crack propagation in poroelastic media, the behavior is asymptotically drained at the crack-tip. Plastic dilatancy is observed to have a strong effect on the distribution and intensity of pore water pressure and to increase its flux towards the crack-tip.     

A comparison of cohesive zone modeling and classical fracture mechanics based on near tip stress field

A mode III crack with a cohesive zone in a power-law hardening material is studied under small scale yielding conditions. The cohesive law follows a softening path with the peak traction at the start of separation process. The stress and strain fields in the plastic zone, and the cohesive traction and separation displacement in the cohesive zone are obtained. The results show that for a modest hardening material (with a hardening exponent N = 0.3), the stress distribution in a large portion of the plastic zone is significantly altered with the introduction of the cohesive zone if the peak cohesive traction is less than two times yield stress, which implies the disparity in terms of the fracture prediction between the classical approach of elastic–plastic fracture mechanics and the cohesive zone approach. The stress distributions with and without the cohesive zone converge when the peak cohesive traction becomes infinitely large. A qualitative study on the equivalency between the cohesive zone approach and the classical linear elastic fracture mechanics indicates that smaller cracks require a higher peak cohesive traction than that for longer cracks if similar fracture initiations are to be predicted by the two approaches.     

On the importance of thermo-elastic cooling in the fracture of glassy polymers at high rates

In a previous thermo-mechanical analysis [Estevez, R., Basu, S., van der Giessen, E., 2005. Analysis of temperature effects near mode I cracks in glassy polymers. Int. J. Fract. 132, 249–273] in which shear yielding of the bulk and failure by crazing were accounted for, we examined which of these two viscoplastic processes contributed to heat in mode I fracture. The present study completes this work by investigating the conditions for thermo-elastic cooling prior to crack propagation as reported experimentally by Rittel [Rittel, D., 1998. Experimental investigation of transient thermo-elastic effects in dynamic fracture. Int. J. Solids Struct. 35, 2959–2973] and Bougaut and Rittel [Bougaut, O., Rittel, D., 2001. On crack tip cooling during dynamic crack propagation. Int. J. Solids Struct. 38, 2517–2532] on high strain rate loading of PMMA. To this end, coupled thermo-mechanical finite element simulations are carried out by accounting for the thermo-elastic source, in addition to the heat sources related to shear yielding and crazing. The bulk as well as cohesive zone parameters for crazing realistically describe PMMA as they are obtained from detailed calibration experiments. Our results show that if significant thermo-elastic cooling has to be observed in the vicinity of the crack tip of a polymeric material, suppression of shear yielding as well as suppression of crazing is necessary. It seems that at these high strain rates a brittle fracture mechanism activated at very high stresses takes over from crazing, or at least that craze initiation occurs for stress levels very different to those for quasi-static conditions.     

Finite element analysis for steady-state hydride-induced fracture in metals by composite model

Delayed hydride cracking, which is observed in hydride-forming metals, due to the precipitation of hydrides near the crack tip, is investigated under conditions of constant temperature and crack velocity, plane strain and small-scale hydride-precipitation. The coupling of the operating physical processes of hydrogen-diffusion, hydride precipitation and material deformation is taken into account. The material is assumed to be an elastic composite made of hydrides and solid solution, with properties depending locally on the volume fraction of the hydrides. In the present analysis, the composite elastic properties have been derived by a generalized self consistent model for particulate composites. With respect to hydride-precipitation, two cases have been considered: (i) precipitation in a homogeneous medium with elastic properties, equal to the effective properties of the composite and (ii) precipitation in an inhomogeneous medium, where the expanding hydride has different elastic properties than those of the surrounding solid solution. The differences between the near-tip field distributions, produced by the two precipitation models, are relatively small. The effect of the hydrogen concentration far from the crack tip, on the near-tip field is also studied. It is shown that for small crack growth velocities, near the threshold stress intensity factor, the remote hydrogen concentration weakly affects the normalized stress distribution in the hydride-precipitation zone, which is controlled by the thermodynamically required hydrostatic stress, under hydrogen chemical equilibrium. However, for values of the applied stress intensity factor and the crack tip velocity, away from the threshold stress intensity factor and crack arrest, the effect of remote hydrogen concentration on the normalized near-tip stress field is strong. Reduction of the remote hydrogen concentration generally leads to reduction of the hydride-precipitation zone and increase of the near-tip stresses. Also reduction of the remote hydrogen concentration leads to distributions closer to those under hydrogen chemical equilibrium.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

A weighted residual method for elastic-plastic analysis near a crack tip and the calculation of the plastic stress intensity factors

In this paper,a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample.The elastic-plastic solutions of the crack tip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field,to calculate the plastic stress intensity factors,are also developed.Therefore,a complete analysis based on the calculation both for the crack tip field and for the whole crack body field is provided.     

A viscoplastic study of crack-tip deformation and crack growth in a nickel-based superalloy at elevated temperature

Viscoplastic crack-tip deformation behaviour in a nickel-based superalloy at elevated temperature has been studied for both stationary and growing cracks in a compact tension (CT) specimen using the finite element method. The material behaviour was described by a unified viscoplastic constitutive model with non-linear kinematic and isotropic hardening rules, and implemented in the finite element software ABAQUS via a user-defined material subroutine (UMAT). Finite element analyses for stationary cracks showed distinctive strain ratchetting behaviour near the crack tip at selected load ratios, leading to progressive accumulation of tensile strain normal to the crack-growth plane. Results also showed that low frequencies and superimposed hold periods at peak loads significantly enhanced strain accumulation at crack tip. Finite element simulation of crack growth was carried out under a constant ΔK-controlled loading condition, again ratchetting was observed ahead of the crack tip, similar to that for stationary cracks.A crack-growth criterion based on strain accumulation is proposed where a crack is assumed to grow when the accumulated strain ahead of the crack tip reaches a critical value over a characteristic distance. The criterion has been utilized in the prediction of crack-growth rates in a CT specimen at selected loading ranges, frequencies and dwell periods, and the predictions were compared with the experimental results.     

Asymptotic crack-tip fields for dynamic crack growth in compressive power-law hardening materials

The effect of material compressibility on the stress and strain fields for a mode-I crack propagating steadily in a power-law hardening material is investigated under plane strain conditions. The plastic deformation of materials is characterized by the J₂ flow theory within the framework of isotropic hardening and infinitesimal displacement gradient. The asymptotic solutions developed by the present authors [Zhu, X.K., Hwang K.C., 2002. Dynamic crack-tip field for tensile cracks propagating in power-law hardening materials. International Journal of Fracture 115, 323–342] for incompressible hardening materials are extended in this work to the compressible hardening materials. The results show that all stresses, strains, and particle velocities in the asymptotic fields are fully continuous and bounded without elastic unloading near the dynamic crack tip. The stress field contains two free parameters σ_eq0 and s₃₃₀ that cannot be determined in the asymptotic analysis, and can be determined from the full-field solutions. For the given values of σ_eq0 and s₃₃₀, all field quantities around the crack tip are determined through numerical integration, and then the effects of the hardening exponent n, the Poisson ratio ν, and the crack growth speed M on the asymptotic fields are studied. The approximate behaviors of the proposed solutions are discussed in the limit of ν → 0.5 or n → ∞.