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.     

Finite element solutions for plane strain mode I crack with strain gradient effects

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom ω_i are introduced in addition to the conventional three translational degrees of freedom u_i. ω_i is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l₁. Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale.     

Crack tip singular fields in ductile crystals with taylor power-law hardening. I: Anti-plane shear

Asymptotic singular solutions of the HRR type are presented for anti-plane shear cracks in ductile crystals. These are assumed to undergo Taylor hardening with a power-law relation between stress and strain at sufficiently large strain. Results are given for several crack orientations in fcc and bcc crystals. The neartip region divides into angular sectors which are the maps of successive flat segments and vertices on the yield locus. Analysis is simplified by use of new general integrals of crack tip singular fields of the HRR type. It is conjectured that the single crystal HRR fields are dominant only over part of the plastic region immediately adjacent to the crack tip, even at small scale yielding, and that their domain of validity vanishes as the perfectly plastic limit is approached. This follows from the fact that while in the perfectly plastic limit the HRR stress states approach the correct discontinuous distributions of the complete elasticideally plastic solutions for crystals (Rice and Nikolic, J. Mech. Phys. Solids33, 595 (1985)), the HRR displacement fields in that limit remain continuous. Instead, the complete elastic-ideally plastic solutions have discontinuous displacements along planar plastic regions emanating from the tip in otherwise elastically stressed material. The approach of the HRR stress fields to their discontinuous limiting distributions is illustrated in graphical plots of results. A case examined here of a fcc crystal with a crack along a slip plane is shown to lead to a discontinuous near-tip stress state even in the hardening regime.Through another limiting process, the asymptotic solution for the near-tip field for an isotropic material is also derived from the present single crystal framework.     

Role of elasticity in elastic–plastic fracture: Analytical bi-linear crack-tip fields and finite element analysis

A solution for Model-I plane strain crack tip fields in a bi-linear elastic–plastic material is presented. The elastic–plastic Poisson's ratio is introduced to characterize the influence of elastic deformation on the near tip constraint. Attention is focused on the distribution of elastic/plastic strain energy in the sensitive region of the forward sector ahead of a crack tip. The present study shows that the elastic strain energy can be higher than the plastic strain energy in this sensitive sector while large amount of the plastic strain energy develops outside this sector around the crack tip. The effect of elastic deformation in this sensitive region on the structure of crack-tip fields is considerable and the assumption in some important solutions for crack-tip fields reported in literature that the elastic deformation is small and can be ignored is therefore not physically reasonable. Besides, finite element analysis is carried out to validate the analytical solution and good agreement between them is found. It is seen that the present solution with T-stress can properly describe the crack-tip fields under various constraints for different specimens and an analytical relation is established between the critical value of J-integral, J_c, and T-stress for elastic–plastic fracture.     

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.     

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 → ∞.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

Field structure at mode III dynamically propagating crack tip in elastic-viscoplastic materials

An elastic-viscoplastic mechanics model is used to investigate asymptotically the mode Ⅲ dynamically propagating crack tip field in elastic-viscoplastic materials. The stress and strain fields at the crack tip possess the same power-law singularity under a linear-hardening condition. The singularity exponent is uniquely determined by the viscosity coefficient of the material. Numerical results indicate that the motion parameter of the crack propagating speed has little effect on the zone structure at the crack tip. The hardening coefficient dominates the structure of the crack-tip field. However, the secondary plastic zone has little influence on the field. The viscosity of the material dominates the strength of stress and strain fields at the crack tip while it does have certain influence on the crack-tip field structure. The dynamic crack-tip field degenerates into the relevant quasi-static solution when the crack moving speed is zero. The corresponding perfectly-plastic solution is recovered from the linear-hardening solution when the hardening coefficient becomes zero.     

Intergranular strain evolution near fatigue crack tips in polycrystalline metals

The deformation field near a steady fatigue crack includes a plastic zone in front of the crack tip and a plastic wake behind it, and the magnitude, distribution, and history of the residual strain along the crack path depend on the stress multiaxiality, material properties, and history of stress intensity factor and crack growth rate. An in situ, full-field, non-destructive measurement of lattice strain (which relies on the intergranular interactions of the inhomogeneous deformation fields in neighboring grains) by neutron diffraction techniques has been performed for the fatigue test of a Ni-based superalloy compact tension specimen. These microscopic grain level measurements provided unprecedented information on the fatigue growth mechanisms. A two-scale model is developed to predict the lattice strain evolution near fatigue crack tips in polycrystalline materials. An irreversible, hysteretic cohesive interface model is adopted to simulate a steady fatigue crack, which allows us to generate the stress/strain distribution and history near the fatigue crack tip. The continuum deformation history is used as inputs for the micromechanical analysis of lattice strain evolution using the slip-based crystal plasticity model, thus making a mechanistic connection between macro- and micro-strains. Predictions from perfect grain-boundary simulations exhibit the same lattice strain distributions as in neutron diffraction measurements, except for discrepancies near the crack tip within about one-tenth of the plastic zone size. By considering the intergranular damage, which leads to vanishing intergranular strains as damage proceeds, we find a significantly improved agreement between predicted and measured lattice strains inside the fatigue process zone. Consequently, the intergranular damage near fatigue crack tip is concluded to be responsible for fatigue crack growth.     

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.     

A reexamination of plasticity-induced crack closure in fatigue crack propagation

The crack closure concept is often used to consider the R-ratio and overload effects on fatigue crack growth. The presumption is that when the crack is closed, the external load produces negligible fatigue damage in the cracked component. The current investigation provides a reassessment of the frequently used concept with an emphasis on the plasticity-induced crack closure. A center cracked specimen made of 1070 steel was investigated. The specimen was subjected to plane-stress mode I loading. An elastic–plastic stress analysis was conducted for the cracked specimens using the finite element method. By applying the commonly used one-node-per-cycle debonding scheme for the crack closure simulations, it was shown that the predicted crack opening load did not stabilize when the extended crack was less than four times of the plastic zone size. The predicted opening load was strongly influenced by the plasticity model used. When the elastic–perfectly plastic (EPP) stress–strain relationship was used together with the kinematic hardening plasticity theory, the predicted crack opening load was found to be critically dependent on the element size of the finite element mesh model. For R = 0, the predicted crack opening load was greatly reduced when the finite element size became very fine. The kinematic hardening rule with the bilinear (BL) stress–strain relationship predicted crack closure with less dependence on the element size. When a recently developed cyclic plasticity model was used, the element size effect on the predicted crack opening level was insignificant. While crack closure may occur, it was demonstrated that cyclic plasticity persisted in the material near the crack tip. The cyclic plasticity was reduced but not negligible when the crack was closed. The traditional approaches may have overestimated the effect of crack closure in fatigue crack growth predictions.     

A modified upper bound approach to limit analysis for plane strain deeply cracked specimens

The classical upper bound approach of limit analysis is based on assumption of rigid blocks of deformation that move between lines of tangential displacement discontinuity. This assumption leads to considerable simplification but often at cost of higher estimate of the actual load. Moreover, in many cases, it does not give a correct shape of the plastic field. In order to overcome these limitations a modified upper bound approach is proposed in this article. The proposed approach is basically an energetic approach but unlike the classical upper bound approach it is capable of including presence of statically governed stress field. As an application, of proposed approach, theoretical plane strain solutions are presented for deeply cracked fracture mechanics specimens (single edge cracked specimen in pure bending – SE (PB), single edge cracked specimen in three-point bending – SE (B), and compact tension – C (T) specimens). Plane strain plasticity problem in rigid elastic–plastic mono-material (homogeneous) was solved to evaluate useful parameters like limit load, plastic eta function (η_p) and plastic rotation factor (r_p) and in bi-material (mismatch welds) to evaluate mismatch limit load, for deeply cracked specimens. New kinematically admissible velocity fields are proposed for SE (B) and C (T) specimens. Proposed theoretical solutions were confirmed by classical slip-line field solutions, wherever available, and by detailed elastic–plastic finite element analysis with Von-Mises yield criterion. Good agreement was found between proposed solutions and results obtained from the classical slip-line field theory and finite element analysis.     

Effect of crack-tip shape on the near-tip field in glassy polymer

The physical nature of a crack tip is not absolutely sharp but blunt with finite curvature. In this paper, the effects of crack-tip shape on the stress and deformation fields ahead of blunted cracks in glassy polymers are numerically investigated under Mode I loading and small scale yielding conditions. An elastic–viscoplastic constitutive model accounting for the strain softening upon yield and then the subsequently strain hardening is adopted and two typical glassy polymers, one with strain hardening and the other with strain softening–rehardening are considered in analysis. It is shown that the profile of crack tip has obvious effect on the near-tip plastic field. The size of near-tip plastic zone reduces with the increase of curvature radius of crack tip, while the plastic strain rate and the stresses near crack tip enhance obviously for two typical polymers. Also, the plastic energy dissipation behavior near cracks with different curvatures is discussed for both materials.     

Asymptotic fields in steady crack growth with linear strain-hardening

The asymptotic stress and velocity fields of a crack propagating steadily and quasi-statically into an elastic-plastic material are presented. The material is characterized by J₂-flow theory with linear strain- hardening. The possibility of reloading on the crack flanks is taken into account. The cases of anti-plane strain (mode III), plane strain (modes I and II), and plane stress (modes I and II) are considered. Numerical results are given for the strength of the singularity and for the distribution of the stress and velocity fields in the plastic loading, elastic unloading and plastic reloading regions, as functions of the strain-hardening parameter. An attempt is made to make a connection with the perfectly-plastic solutions in the limit of vanishing strain-hardening.     

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.     

Energy balance and identification of hardening moduli in plastic deformation processes

The hardening moduli H_r and H_d of plastic deformation associated with the free energy and dissipation function in a representative material element are defined analytically and specified experimentally for three materials. Besides the stress–strain curve and work expended during the deformation process, variation of the hardening moduli with plastic deformation is also determined for austenitic steel, austenitic-ferritic steel and Fe–Si alloy.     

A generalized formulation of the concept of nonhardening region

A stable hysteresis loop is realized under a cyclic loading between fixed strain limits for metals. To describe this phenomenon the idea of memory hypersphere was proposed by Chabocheet al. [1979] and it has been extended by Ohno [1982] as the concept of nonhardening (strain) region, which assumes that an isotropic hardening does not proceed when a plastic strain exists inside a hypersphere in the plastic strain space. A generalized formulation for this concept is presented in this article, which concerns isotropic-kinematic hardening materials, not limited to metals, undergoing a large deformation with a material rotation.     

平面幂律塑性Ⅰ型裂纹尖端应力场全解

在航空航天、船舶、石油管道和核电等领域,服役结构或部件在长期极端条件下运行,不可避免地会产生裂纹,因此,为研究含裂纹结构的准静态断裂行为,必须了解裂纹尖端附近区域的应力应变场特点.对于幂律材料裂纹构元,研究平面应变和平面应力条件下Ⅰ型裂纹尖端应力场的解析分布.基于能量密度等效和量纲分析,推导了能量密度中值点代表性体积单元(representative volume element, RVE)的等效应力解析方程,并定义其为应力因子,进而针对有限平面应变和平面应力紧凑拉伸(compact tension, CT)试样和单边裂纹弯曲(single edge bend, SEB)试样,以应力因子作为应力特征量,并构造用于表征裂尖等效应力等值线的蝶翅轮廓式和扇贝轮廓式三角特殊函数,提出描述幂律塑性条件下平面I型裂纹尖端应力场的半解析模型.该半解析模型形式简单,对CT和SEB试样的裂尖应力场的预测结果与有限元分析的结果比较表明,两者之间均密切吻合,模型公式可直接用于预测Ⅰ型裂纹尖端应力分布,方便于断裂安全评价和理论发展.     

Singular plastic fields in wedge indentation of pressure sensitive solids

The paper examines singular plastic fields induced near the tip of a wedge indentating a pressure sensitive solid. Plane strain conditions are assumed and material response is modelled by the small strain Drucker–Prager rigid/plastic constitutive law. A standard separation of variables solution is numerically generated for pure power-law hardening. Three possible measures of wall roughness are studied with an attempt to expose the coupling between wall friction and material pressure sensitivity. Sample calculations illustrate that stress singularity decreases with increasing friction, wedge angle and hardening exponent, but increases with pressure sensitivity. At large values of the hardening exponent, when the material is nearly perfectly plastic, effective stress contours approach the slip line limit. The concept of indentation index is introduced as a possible estimate for average indentation pressure.