High order asymptotic field for plane stress crack problem

This paper presents an exact analysis for high order asymptotic field of the plane stress crack problem. It has been shown that the second order asymptotic field is not an independent eigen field and should be matched with the elastic strain term of the first order asymptotic field. The second order stress field ahead of the crack tip is quite small compared with the first order stress field. The stress field ahead of crack tip is characterized by the HRR field. Hence the J integral can be used as a criterion for crack initiation.

Collapsed isoparametric element as a singular element for a crack normal to the bi-material interface

It is shown that the variable power singularity of the strain field at the crack tip can be obtained by the simple technique of collapsing quadrilateral isoparametric elements into triangular elements around the crack tip and adequately shifting the side-nodes adjacent to this crack tip. The collapsed isoparametric elements have the desired singularity at crack tip along any ray. The strain expressions for a single element have been derived and in addition to the desired power singularity, additional singularities are revealed. Numerical examples have shown that triangular elements formed by collapsing one side lead to excellent results.

Crack growth by grain boundary cavitation in the transient and extensive creep regimes

Crack growth due to cavity growth and coalescence along grain boundaries is analyzed under transient and extensive creep conditions in a compact tension specimen. Account is taken of the finite geometry changes accompanying crack tip blunting. The material is characterized as an elastic-power law creeping solid with an additional contribution to the creep rate arising from a given density of cavitating grain boundary facets. All voids are assumed present from the outset and distributed on a given density of cavitating grain boundary facets. The evolution of the stress fields with crack growth under three load histories is described in some detail for a relatively ductile material. The full-field plane strain finite element calculations show the competing effects of stress relaxation due to constrained creep, diffusion and crack tip blunting. and of stress increase due to the instantaneous elastic response to crack growth. At very high crack growth rates the Hui-Riedel fields dominate the crack tip region. However. the high growth rates are not sustained for any length of time in the compact tension geometry analyzed. The region of dominance of the Hui-Riedel field shrinks rapidly so that the near-tip fields are controlled by the HRR-type field shortly after the onset of crack growth. Crack growth rates under various conditions of loading and spanning the range of times from small scale creep to extensive creep are obtained. We show that there is a strong similarity between crack growth history and the behaviour of the C(t) and C(t) parameters. so that crack growth rates correlate rather well with C(t) and C(t). A relatively brittle material is also considered that has a very different near-tip stress field and crack growth history.

Fatigue crack-growth rate in ferrite martensite dual-phase steel

An empirical study is made on the fatigue crack growth rate in ferrite-martensite dual-phase (FMDP) steel. Particular attention is given to the effect of ferrite content in the range of 24.2% to 41.5% where good fatigue resistance was found at 33.8%. Variations in ferrite content did not affect the crack growth rate when plotted against the effective stress intensity factor range  which was assumed to follow a linear relation with the crack tip stress intensity factor range ΔK. A high  corresponds to uniformly distributed small size ferrite and martensite. No other appreciable correlation could be ralated to the microstructure morphology of the FMDP steel. The closure stress intensity factor , however, is affected by the ferrite content with  reaching a maximum value of 0.7. In general, crack growth followed the interphase between the martensite and ferrite.

Dividing the fatigue crack growth process into Stage I and II where the former would be highly sensitive to changes in ΔK and the latter would increase with ΔK depending on the  ratio. The same data when correlated with the strain energy density factor range ΔS showed negligible dependence on mean stress or R ratio for Stage I crack growth. A parameter α involving the ratio of ultimate stress to yield stress, percent reduction of area and R is introduced for Stage II crack growth so that the  data for different R would collapse onto a single curve with a narrow scatter band when plotted against αΔS.

Fracture analysis for 3-dimensional bodies with a surface crack

This paper deals with fracture analyses in 3-dimensional bodies containing a surface crack. A general solution of stress-strain fields at crack tip is proposed. Based on the stress-strain fields obtained, a high-order 3-dimensional special element is established to calculate the stress intensity factors in a plate with a surface crack. The variation of stress intensity factors with geometric parameters is investigated.

High-order asymptotic field of tensile plane-strain nonlinear crack problems

In this paper, the governing equations and the analytical method of the secondorder asymptotic field for the plane-straln crack problems of mode I have been presented. The numerical calculation has been carried out. The amplitude coefficients k2 of the second term of the asymptotic field have been determined after comparing the present results with some fine results of the finite element calculation. The variation of coefficients k2 with changes of specimen geometry and developments of plastic zone have been analysed. It is shown that the second term of the asymptotic field has significant influence on the near-crack-tip field. Therefore, we may reasonably argue that both the J-integral and the coefficient k2 can beeome two characterizing parameters for crack initiation.

A Finite Element Analysis of a Crack Penetrating or Deflecting into an Interface in a Thin Laminate

The crack tip driving force of a crack growing from a pre-crack that is perpendicular to and terminating at an interface between two materials is investigated using a linear fracture mechanics theory. The analysis is performed both for a crack penetrating the interface, growing straight ahead, and for a crack deflecting into the interface. The results from finite element calculations are compared with asymptotic solutions for infinitesimally small crack extensions. The solution is found to be accurate even for fairly large amounts of crack growth. Further, by comparing the crack tip driving force of the deflected crack with that of the penetrating crack, it is shown how to control the path of the crack by choosing the adhesion of the interface relative to the material toughness.

Discontinuous growth mechanisms of mode II crack under high-speed impact conditions

A recoverable plate impact testing technology has been used for studying the growth mechanisms of mode II crack. The results show that interactions of microcracks ahead of a crack tip cause the crack growth unsteadily. Failure mode transitions of materials were observed. Based on the observations, a discontinuous crack growth model was established. Analysis shows that the shear crack grows unsteady as the growth speed is between the Rayleigh wave speed c(R) and the shear wave speed c(s); however, when the growth speed approaches root 2c(s), the crack grows steadily. The transient microcrack growth makes the main crack speed to jump from subsonic to intersonic and the steady growth of all the sub-cracks leads the main crack to grow stably at an intersonic speed.

Subsonic interface crack with crack face contact

A steady-state subsonic interface crack propagating between an elastic solid and a rigid substrate with crack face contact is studied. Two cases with respective to the contact length are considered, i.e., semi-infinite and finite crack face contact. Different from a stationary or an open subsonic interface crack, stress singularity at the crack tip in the present paper is found to be non-oscillatory. Furthermore, in the semi-infinite contact case, the singularity of the stress field near the crack tip is less than 1/2. In the finite contact case, no singularity exists near the crack tip, but less than 1/2 singularity does at the end of the contact zone. In both cases, the singularity depends on the linear contact coefficient and the crack speed. Asymptotic solutions near the crack tip are given and analyzed. In order to satisfy the contact conditions, reasonable region of the linear contact coefficient is found. In addition, the solution predicts a non-zero-energy dissipation rate due to crack face contact.

A discrete dislocation analysis of fatigue crack growth

Cyclic loading of a plane strain mode I crack under small scale yielding is analyzed using discrete dislocation dynamics. The dislocations are all of edge character, and are modeled as line singularities in an elastic solid. At each stage of loading, superposition is used to represent the solution in terms of solutions for edge dislocations in a half-space and a non-singular complementary solution that enforces the boundary conditions, which is obtained from a linear elastic, finite element solution. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated into the formulation through a set of constitutive rules. An irreversible relation between the opening traction and the displacement jump across a cohesive surface ahead of the initial crack tip is also specified, which permits crack growth to emerge naturally. It is found that crack growth can occur under cyclic loading conditions even when the peak stress intensity factor is smaller than the stress intensity required for crack growth under monotonic loading conditions; however below a certain threshold value of ΔKI no crack growth was seen.

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113-2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent -1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.

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.

Steady-state motion of a Mode-III crack on imperfect interfaces

Mishuris, Gennady; Movchan, N.V.; Movchan, A.B., (2006) 'Steady-state motion of a Mode-III crack on imperfect interfaces', Quarterly Journal of Mechanics and Applied Mathematics 59(4) pp.487-516 RAE2008

A multi-scale atomistic-continuum modelling of crack propagation in a two-dimensional macroscopic plate

A novel multi-scale seamless model of brittle-crack propagation is proposed and applied to the simulation of fracture growth in a two-dimensional Ag plate with macroscopic dimensions. The model represents the crack propagation at the macroscopic scale as the drift-diffusion motion of the crack tip alone. The diffusive motion is associated with the crack-tip coordinates in the position space, and reflects the oscillations observed in the crack velocity following its critical value. The model couples the crack dynamics at the macroscales and nanoscales via an intermediate mesoscale continuum. The finite-element method is employed to make the transition from the macroscale to the nanoscale by computing the continuum-based displacements of the atoms at the boundary of an atomic lattice embedded within the plate and surrounding the tip. Molecular dynamics (MD) simulation then drives the crack tip forward, producing the tip critical velocity and its diffusion constant. These are then used in the Ito stochastic calculus to make the reverse transition from the nanoscale back to the macroscale. The MD-level modelling is based on the use of a many-body potential. The model successfully reproduces the crack-velocity oscillations, roughening transitions of the crack surfaces, as well as the macroscopic crack trajectory. The implications for a 3-D modelling are discussed.

Multiscale numerical modelling of crack propagation in two-dimensional metal plate

A novel multiscale model of brittle crack propagation in an Ag plate with macroscopic dimensions has been developed. The model represents crack propagation as stochastic drift-diffusion motion of the crack tip atom through the material, and couples the dynamics across three different length scales. It integrates the nanomechanics of bond rupture at the crack tip with the displacement and stress field equations of continuum based fracture theories. The finite element method is employed to obtain the continuum based displacement and stress fields over the macroscopic plate, and these are then used to drive the crack tip forward at the atomic level using the molecular dynamics simulation method based on many-body interatomic potentials. The linkage from the nanoscopic scale back to the macroscopic scale is established via the Ito stochastic calculus, the stochastic differential equation of which advances the tip to a new position on the macroscopic scale using the crack velocity and diffusion constant obtained on the nanoscale. Well known crack characteristics, such as the roughening transitions of the crack surfaces, crack velocity oscillations, as well as the macroscopic crack trajectories, are obtained.