Discrete dislocation simulations of the development of a continuum plastic zone ahead of a stationary Mode III crack

We present results from discrete dislocation simulations showing the development of the plastic zone in front of a Mode III crack under constant load. We find that the equilibrated zone is circular, in agreement with continuum mechanics predictions of the elastic-perfectly plastic Mode III crack. The size of the equilibrated zone scales as the square of the applied load (K_III), also in agreement with the continuum results. The zone approaches saturation exponentially, with a time that scales as K_III²/σ_p³, where σ_p is the Peierls stress. These results delineate conditions under which the classical, continuum predictions of elastic-plastic fracture mechanics are applicable.     

An experimental investigation of crack-path directional stability

Using the criterion that a crack will extend perpendicular to the maximum circumferential stress,σ _θ, we show that the directional stability of crack growth is governed by the location of microcrack initiation ahead of the crack tip. At distances greater than a geometrical radiusr _o, the maximum value ofσ _θ deviates from the position of symmetry. Thus, if we assume that the physical processes involved in fracture lead to crack initiation at a distancer _c ahead of the crack tip, the criterion for directional stability isr _o>r _c. Experimental and theoretical values ofr _o verify that, asr _o becomes small, the crack's directional stability deteriorates. Observing that a lengthwise compressive stress increasesr _o, a center-cracked specimen was developed which allows the application of controlled lengthwise compression independently of the opening-mode load. A detailed photoelastic analysis of the specimen has provided the value ofr _o as a function of the crack length. The value ofr _o is then compared with the expected microcrack initiation distances in ductile fracture. By applying sufficient lengthwise compression, we are able to make the crack grow straight and obtain numerous data points from this specimen which would otherwise be directionally unstable. The results indicate that, as the total lengthwise tensile stress at the crack tip increases, the fracture toughness also increases. Using this information we can then adjustK _Ic for zero lengthwise loading and obtain a geometry independent fracture toughness.     

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.     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Tensile-shear transition in mixed mode I/III fracture

The propensity of the transition of fracture type in either brittle or ductile cracked solid under mixed-mode I and III loading conditions is investigated. A fracture criterion based on the competition of the maximum normal stress and maximum shear stress is utilized. The prediction of the fracture type is determined by comparing τ_max/σ_max at a critical distance from the crack tip to the material strength ratio τ_C/σ_C, i.e., (τ_max/σ_max)<(τ_C/σ_C) for tensile fracture and (τ_max/σ_max)>(τ_C/σ_C) for shear fracture, where σ_C (τ_C) is the fracture strength of materials in tension (shear). Mixed mode I/III fracture tests were performed using circumferentially notched cylindrical bars made of PMMA and 7050 aluminum alloy. Fracture surface morphology of the specimens reveals that: (1) for the brittle material, PMMA, only tensile type of fracture occurs, and (2) for the ductile material, 7050 aluminum alloy, either tensile or shear type of fracture occurs depending on the mode mixity. The transition (in ductile material) or non-transition (in brittle material) of the fracture type and the fracture path observed in experiments were properly predicted by the theory. Additional test data from open literature are also included to validate the proposed theory.     

A new finite element method for strain gradient theories and applications to fracture analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J₂ deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.     

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

Fracture of bicrystal metal/ceramic interfaces: A study via the mechanism-based strain gradient crystal plasticity theory

Two continuum mechanical models of crystal plasticity theory namely, conventional crystal plasticity theory and mechanism-based crystal plasticity theory, are used to perform a comparative study of stresses that are reached at and ahead of the crack tip of a bicrystal niobium/alumina specimen. Finite element analyses are done for a stationary crack tip and growing cracks using a cohesive modelling approach. Using mechanism-based strain gradient crystal plasticity theory the stresses reached ahead of the crack tip are found to be two times larger than the stresses obtained from conventional crystal plasticity theory. Results also show that strain gradient effects strongly depend on the intrinsic material length to the size of plastic zone ratio (l/R₀). It is found that the larger the (l/R₀) ratio, the higher the stresses reached using mechanism-based strain gradient crystal plasticity theory. An insight into the role of cohesive strength and work of adhesion in macroscopic fracture is also presented which can be used by experimentalists to design better bimaterials by varying cohesive strength and work of adhesion.     

Crack tip field andJ-integral with strain gradient effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventionalJ ₂ deformation theory. No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory. Two typical crack problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-IK-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory. The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, theJ-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases. The project supported by the National Natural Science Foundation of China (19704100 and 10202023) and the Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20)     

Effect of temperature on the fracture modes of BHW-35 steel under mode-II loading

Finite element method (FEM) has been used to analyze the stress and strain fields and the stress tri-axial levels around the tip of the crack under mode- II loading. The results show that: under mode- II loading, the direction of the maximum tensile stress and that of the maximum tri-axial levels (R _o ) exist at an angle of –75. 3° from the original crack plane; the maximum shear stress andR _o = 0 exist along the original crack plane.Mode- II loading experiment using BHW-35 steel at different temperatures show that there are two kinds of fracture mode, opening mode (or tensile mode) and sliding mode (or shear mode). A decrease in temperature causes the fracture mode to change from shear mode to tensile mode. For BHW-35 steel, this critical temperature is about –90 C. Actually, under any kind of loading mode (mode I . mode II , mode III or mixed mode), there always exist several kinds of potenital fracture modes (for example, opening mode, sliding mode, tearing mode or mixed mode). The effect of temperature under mode- II loading is actually related to the change of the elastic-plastic properties of the material.     

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.     

On the effects of characteristic lengths in bending and torsion on Mode III crack in couple stress elasticity

The problem of a stationary semi-infinite crack in an elastic solid with microstructures subject to remote classical K_III field is investigated in the present work. The material behavior is described by the indeterminate theory of couple stress elasticity developed by Koiter. This constitutive model includes the characteristic lengths in bending and torsion and thus it is able to account for the underlying microstructure of the material as well as for the strong size effects arising at small scales. The stress and displacement fields turn out to be strongly influenced by the ratio between the characteristic lengths. Moreover, the symmetric stress field turns out to be finite at the crack tip, whereas the skew-symmetric stress field displays a strong singularity. Ahead of the crack tip within a zone smaller than the characteristic length in torsion, the total shear stress and reduced tractions occur with the opposite sign with respect to the classical LEFM solution, due to the relative rotation of the microstructural particles currently at the crack tip. The asymptotic fields dominate within this zone, which however has limited physical relevance and becomes vanishing small for a characteristic length in torsion of zero. In this limiting case the full-field solution recovers the classical K_III field with square-root stress singularity. Outside the zone where the total shear stress is negative, the full-field solution exhibits a bounded maximum for the total shear stress ahead of the crack tip, whose magnitude can be adopted as a measure of the critical stress level for crack advancing. The corresponding fracture criterion defines a critical stress intensity factor, which increases with the characteristic length in torsion. Moreover, the occurrence of a sharp crack profile denotes that the crack becomes stiffer with respect to the classical elastic response, thus revealing that the presence of microstructures may shield the crack tip from fracture.     

Characterization of crack-tip strain singularity in brittle-microcracking materials by means of the photoelastic-coating method

The strain fields ahead of crack tips in rock, mortar, and graphite were measured using a photoelastic coating method. A transparent ferroelectric ceramic was used as a coating material to observe the photoelastic fringe pattern. A coating plate of 110–150 μm thick was placed on single-edge-notch specimens, and photoelasticity experiments were conducted in three-point bending under a polarizing microscope. The results show that the principal-strain difference ahead of the crack tip is given by $$\Delta \in = \Delta \in _o [(J/\sigma _{ult} )/r]^m $$ whereσ _ult is the ultimate tensile strength,r is the distance from the crack tip, and9? _o andm are material constants. Based on this observation, the use of theJ _Ic concept in determining the fracture toughness of brittle-microcracking materials is discussed.     

Plane-strain crack problems in microstructured solids governed by dipolar gradient elasticity

The present study aims at determining the elastic stress and displacement fields around the tips of a finite-length crack in a microstructured solid under remotely applied plane-strain loading (mode I and II cases). The material microstructure is modeled through the Toupin-Mindlin generalized continuum theory of dipolar gradient elasticity. According to this theory, the strain-energy density assumes the form of a positive-definite function of the strain tensor (as in classical elasticity) and the gradient of the strain tensor (additional term). A simple but yet rigorous version of the theory is employed here by considering an isotropic linear expression of the elastic strain-energy density that involves only three material constants (the two Lamé constants and the so-called gradient coefficient). First, a near-tip asymptotic solution is obtained by the Knein-Williams technique. Then, we attack the complete boundary value problem in an effort to obtain a full-field solution. Hypersingular integral equations with a cubic singularity are formulated with the aid of the Fourier transform. These equations are solved by analytical considerations on Hadamard finite-part integrals and a numerical treatment. The results show significant departure from the predictions of standard fracture mechanics. In view of these results, it seems that the classical theory of elasticity is inadequate to analyze crack problems in microstructured materials. Indeed, the present results indicate that the stress distribution ahead of the crack tip exhibits a local maximum that is bounded. Therefore, this maximum value may serve as a measure of the critical stress level at which further advancement of the crack may occur. Also, in the vicinity of the crack tip, the crack-face displacement closes more smoothly as compared to the standard result and the strain field is bounded. Finally, the J-integral (energy release rate) in gradient elasticity was evaluated. A decrease of its value is noticed in comparison with the classical theory. This shows that the gradient theory predicts a strengthening effect since a reduction of crack driving force takes place as the material microstructure becomes more pronounced.     

Direct Evaluation of Cohesive Law in Mode I of Pinus pinaster by Digital Image Correlation

The direct identification of the cohesive law in pure mode I of Pinus pinaster is addressed in this work. The approach couples the double cantilever beam (DCB) test with digital image correlation (DIC). Wooden beam specimens loaded in the radial-longitudinal (RL) fracture propagation system are used. The strain energy release rate in mode I (G _I) is uniquely determined from the load–displacement curve by means of the compliance-based beam method (CBBM). This method relies on the concept of equivalent elastic crack length (a _eq) and therefore does not require the monitoring of crack propagation during test. DIC measurements are processed with two different purposes. Firstly, the physical evidence of a _eq is discussed with regard to actual estimation of the crack length based on post-processing full-field displacement measurements. Secondly, the crack tip opening displacement in mode I (w _I) is determined from the displacements near the initial crack tip. The cohesive law in mode I (σ _I???w _I) is then identified by numerical differentiation of the G _I???w _I relationship. The methodology and accuracy on this reconstruction are addressed. Moreover, the proposed procedure is validated by finite element analyses including cohesive zone modelling. It is concluded that the proposed data reduction scheme is adequate for assessing the cohesive law in pure mode I of P. pinaster.     

An Experimental Study of Mixed Mode Crack Initiation and Growth in Functionally Graded Materials

Quasi-static mixed mode crack initiation and growth in functionally graded materials (FGMs) was studied through fracture experiments on polymer-based FGMs manufactured by selective ultraviolet irradiation poly(ethylene carbon monoxide)—a photo-sensitive copolymer that becomes more brittle and stiffer under ultraviolet irradiation. The objective of the study was to determine whether crack kinking criteria for homogeneous materials, e.g., maximum hoop stress criterion, also hold for FGMs. Single edge notched tension specimens with different spatial variations of Young's modulus, failure stress and failure strain, were tested. Near tip mode mixity was introduced either by inclining the crack to the remote loading direction, as in the case of homogeneous materials, or to the direction of material gradient, or both. A full-field digital image correlation technique was used to measure in real-time the displacement field around the crack tip while it propagated through the graded material, and to extract the fracture parameters of stress intensity factor K _I and K _II , and the T-stress. It was found that the nonsingular T-stress term in the asymptotic expansion for stresses plays a very important role in accurately measuring fracture parameters. It was also found that the maximum tangential stress criterion can be applied to the case of FGMs to predict crack kinking provided that the effect of the T-stress is accounted for and the process zone size is small compared to the intrinsic material gradient length scale. However, for accurate crack path prediction at a length scale comparable to the material gradient, detailed material property information is required. In general, the crack will propagate towards a region that exhibits less fracture toughness, but, unlike the case of homogeneous materials, along a path where K _II is not necessarily equal to zero.     

Finite-deformation analysis of the crack-tip fields under cyclic loading

The plane-strain crack subjected to mode I cyclic loading under small scale yielding was analysed. The influence of the load range, load ratio and overload on the near-tip deformation-, stress- and strain-fields was studied. Although the near-tip zones of appreciable cyclic plastic flow for all loading regimes matched closely one another, when scaled with (ΔK/σ_Y)², the activities of plastic flow within them manifested dependence on K_max and K_min, as well as on overload. Cyclic trajectories of the crack-tip opening displacement (CTOD) converged to stable self-similar loops of the sizes proportional to ΔK² and positions in CTOD-K plane dependent on the maximum K along the whole loading route, including an overload. Computed near-tip deformation evidenced plastic crack advance, this way visualising of the Laird–Smith concept of fatigue cracking. This crack growth by blunting-resharpening accelerated with rising ΔK and was halted by an overload. Crack closure upon unloading had no place. The affinities were revealed between computed near-tip stress–strain variables and the experimental trends of the fatigue crack growth rate, such as its dependence on K_max and K_min (or ΔK and K_max), and retardation by overload. Thus, the effects of loading parameters on fatigue cracking, hitherto associated with crack closure, are attributable to the stress–strain fields in front of it as the direct drives of the key fatigue constituents – damage accumulation and bond breaking.     

Stress distribution about a slowly growing crack determined by photoelastic-coating method

An experimental stress-analysis technique using a birefringent coating is reported for determining the stress distribution about a slowly growing crack. The maximum error of the test method for a large strain gradient is found to be less than 10 percent. For a plate with an internal crack, the experimentally determined stress distribution compares favorably with two numerical solutions. Comparison of stresses about an internal or double-edge crack to those about a single-edge crack indicates that the isochromatics bend over to about 45 deg with the plane of the crack in the former and are inclined at about 60 deg in the latter. Also, the stresses for a single-crack tip vary as the inverse square root of the radius, while the stresses for a double-crack tip follow anr ^?1/4 law more closely.     

Effects of characteristic material lengths on mode III crack propagation in couple stress elastic–plastic materials

The asymptotic fields near the tip of a crack steadily propagating in a ductile material under Mode III loading conditions are investigated by adopting an incremental version of the indeterminate theory of couple stress plasticity displaying linear and isotropic strain hardening. The adopted constitutive model is able to account for the microstructure of the material by incorporating two distinct material characteristic lengths. It can also capture the strong size effects arising at small scales, which results from the underlying microstructures. According to the asymptotic crack tip fields for a stationary crack provided by the indeterminate theory of couple stress elasticity, the effects of microstructure mainly consist in a switch in the sign of tractions and displacement and in a substantial increase in the singularity of tractions ahead of the crack-tip, with respect to the classical solution of LEFM and EPFM. The increase in the stress singularity also occurs for small values of the strain hardening coefficient and is essentially due to the skew-symmetric stress field, since the symmetric stress field turns out to be non-singular. Moreover, the obtained results show that the ratio η introduced by Koiter has a limited effect on the strength of the stress singularity. However, it displays a strong influence on the angular distribution of the asymptotic crack tip fields.     

Mode I crack tip field with strain gradient effects

A plane strain mode I crack tip field with strain gradient effects is investigated. A new strain gradient theory is used. An elastic-power law hardening strain gradient material is considered and two hardening laws, i. e. a separation law and an integration law are used respectively. As for the material with the separation law hardening, the angular distributions of stresses are consistent with the HRR field, which differs from the stress results^[19]; the angular distributions of couple stresses are the same as the couple stress results^[19]. For the material with the integration law hardening, the stress field and the couple stress field can not exist simultaneously, which is the same as the conclusion^[19], but for the stress dominated field, the angular distributions of stresses are consistent with the HRR field; for the couple stress dominated field, the angular distributions of couple stresses are consistent with those in Ref. [19]. However, the increase in stresses is not observed in strain gradient plasticity because the present theory is based on the rotation gradient of the deformation only, while the crack tip field of mode I is dominated by the tension gradient, which will be shown in another paper. Supported by the National Science Foundation of China (No. 19704100), Science Foundation of Chinese Academy of Sciences (Project KJ951-1-20), CAS K. C. Wong Post-doctoral Research Award Fund and the Post Doctoral Science Fund of China.