Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

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 omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(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(1). 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. (C) 2002 Elsevier Science Ltd. All rights reserved.

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 length l, 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.

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.

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

Asymptotic analysis of plane-strain mode I steady-state crack growth in transformation toughening ceramics (II)

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.

Plastic zone in front of mode I crack in phenolphthalein polyether ketone

The plastic zone size and crack opening displacement of phenolphthalein polyether ketone (PEK-C) at various conditions were investigated. Both of them increase with increasing temperature (decreasing strain rate), i.e. yield stress steadily falls. Thus, the mechanism increasing the yield stress leads to increased constraint in the crack tip and a corresponding reduction in the crack opening displacement and the plastic deformation zone. The effect of the plastic deformation on the fracture toughness is also discussed.

DETERMINATION OF INITIATION VALUE OF MODE-I INTERLAMINAR FRACTURE-TOUGHNESS OF UNIDIRECTIONAL POLYMER COMPOSITES

The use of interlaminar fracture tests to measure the delamination resistance of unidirectional composite laminates is now widespread. However, because of the frequent occurrence of fiber bridging and multiple cracking during the tests, it leads to artificially high values of delamination resistance, which will not represent the behavior of the laminates. Initiation fracture from the crack starter, on the other hand, does not involve bridging, and should be more representative of the delamination resistance of the composite laminates. Since there is some uncertainty involved in determining the initiation value of delamination resistance in mode I tests in the literature, a power law of the form G(IC) = A.DELTA alpha(b) (where G(IC) is mode I interlaminar fracture toughness and DELTA alpha is delamination growth) is presented in this paper to determine initiation value of mode I interlaminar fracture toughness. It is found that initiation values of the mode I interlaminar fracture toughness, G(IC)(ini), can be defined as the G(IC) value at which 1 mm of delamination from the crack starter has occurred. Examples of initiation values determined by this method are given for both carbon fiber reinforced thermoplastic and thermosetting polymers.

A STUDY OF MODE-I DELAMINATION RESISTANCE OF A THERMOPLASTIC COMPOSITE

This paper describes the mode I delamination behaviour of a unidirectional carbon-fibre/poly(phenylene ether ketone)(PEK-C) composite. Tests have been performed on double cantilever beam (DCB) specimens. Several data reduction schemes are used to obtain the critical strain energy release rate, G(IC), and the results are compared. It is shown that when using a DCB test to determine the fracture toughness, corrections must be employed. The experimental methods have been described for ascertaining the correction terms, and the results are consistent after modification. Some of the authors' results are different from those of other authors, particularly the negative correction term for crack length, the larger exponent (n > 3) in the relationship C = Ra(n), and decrements of flexural modulus with the crack growth when using the simple beam theory to predict the bending behaviour of DCB specimens. The possible reasons are discussed.

Investigation Of Mode Ii Crack Growth Following A Very High Speed Impact

A recoverable plate impact testing technology has been developed for studying fracture mechanisms of mode II crack. With this technology, a single duration stress pulse with submicrosecond duration and high loading rates, up to 10(8) MPam(1/2)s(-1), can be produced. Dynamic failure tests of Hard-C 60# steel were carried out under asymmetrical impacting conditions with short stress-pulse loading. Experimental results show that the nucleation and growth of several microcracks ahead of the crack tip, and the interactions between them, induce unsteady crack growth. Failure mode transitions during crack growth, both from mode I crack to mode II and from brittle to ductile fracture, were observed. Based on experimental observations, a discontinuous crack growth model was established. Analysis of the crack growth mechanisms using our model shows that the shear crack extension is unsteady when the extending speed is between the Rayleigh wave speed c(R) and the shear wave speed c(S). However, when the crack advancing speed is beyond c(S), the crack grows at a steady intersonic speed approaching root 2c(S). It also shows that the transient mechanisms, such as nucleation, growth, interaction and coalescence among microcracks, make the main crack speed jump from subsonic to intersonic and the steady growth of all the subcracks causes the main crack to grow at a stable intersonic speed.

Mixed-mode fracture mechanisms near the fatigue threshold of aisi-316 stainless-steel

Near threshold, mixed mode (I and II), fatigue crack growth occurs mainly by two mechanisms, coplanar (or shear) mode and branch (or tensile) mode. For a constant ratio of ΔK_I/ΔK_II the shear mode growth shows a self-arrest character and it would only start again when ΔK_I and ΔK_II are increased. Both shear crack growth and the early stages of tensile crack growth, are of a crystallographic nature; the fatigue crack proceeds along slip planes or grain boundaries. The appearance of the fracture surfaces suggest that the mechanism of crack extension is by developing slip band microcracks which join up to form a macrocrack. This process is thought to be assisted by the nature of the plastic deformation within the reversed plastic zone where high back stresses are set up by dislocation pile-ups against grain boundaries. The interaction of the crack tip stress field with that of the dislocation pile-ups leads to the formation of slip band microcracks and subsequent crack extension. The change from shear mode to tensile mode growth probably occurs when the maximum tensile stress and the microcrack density in the maximum tensile plane direction attain critical values.

The Flow Theory of Mechanism-Based Strain Gradient Plasticity

The flow theory of mechanism-based strain gradient (MSG) plasticity is established in this paper following the same multiscale, hierarchical framework for the deformation theory of MSG plasticity in order to connect with the Taylor model in dislocation mechanics. We have used the flow theory of MSG plasticity to study micro-indentation hardness experiments. The difference between deformation and flow theories is vanishingly small, and both agree well with experimental hardness data. We have also used the flow theory of MSG plasticity to investigate stress fields around a stationary mode-I crack tip as well as around a steady state, quasi-statically growing crack tip. At a distance to crack tip much larger than dislocation spacings such that continuum plasticity still applies, the stress level around a stationary crack tip in MSG plasticity is significantly higher than that in classical plasticity. The same conclusion is also established for a steady state, quasi-statically growing crack tip, though only the flow theory can be used because of unloading during crack propagation. This significant stress increase due to strain gradient effect provides a means to explain the experimentally observed cleavage fracture in ductile materials [J. Mater. Res. 9 (1994) 1734, Scripta Metall. Mater. 31 (1994) 1037; Interface Sci. 3(1996) 169].

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 J2 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.

Dislocations emission and crack extension at the atomistic crack tip in body-centered-cubic metal Mo

The behaviors of a crack in body-centered-cubic metal Mo under different loading modes were studied using the molecular dynamics method. Dislocation emission was observed near the crack tip in response to mode II loading with theta = 0 degrees in which theta is the inclination angle of the slip plane with respect to the crack plane, and two full dislocations were observed at the stress level of K-II = 1.17 MPa m(1/2) without any evidence of crack extension. Within the range of 0 degrees less than or equal to theta less than or equal to 45 degrees, crack extension was observed in response to mode I loading, and the effect of crystal orientation on the crack propagation was studied, The crack propagated along the [111] slip direction without any evidence of dislocations emission.

Crack Tip Field and J-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 conventional J_2 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-I K-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, the J-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases.