Stress intensity factor histories of a steadily propagating crack under 3-D combined mode traction

Elastodynamic stress intensity factor histories of an unbounded solid containing a semi-infinite plane crack that propagates at a constant velocity under 3-D time-independent combined mode loading are considered. The fundamental solution, which is the response of point loading, is obtained. Then, stress intensity factor histories of a general loading system are written out in terms of superposition integrals. The methods used here are the Laplace transform methods in conjunction with the Wiener-Hopf technique.

A contact crack problem in an infinite plate

The problem of an infinite plate with crack of length 2a loaded by the remote tensile stress P and a pair of concentrated forces Q is discussed. The value of the force Q for the initial contact of crack face is investigated and the contact length elevated, while the Q force increases. The problem is solved assuming that the stress intensity factor vanishes at the end point of the contact portion. By the Fredholm integral equation for the multiple cracks, the reduction of stress intensity factor due to Q is found. (C) 1999 Elsevier Science Ltd. All rights reserved.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90degrees and 180degrees is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress sigma(infinity) and electric field E-infinity or (sigma(infinity), E-infinity) is replaced by strain e and electric displacement D-infinity or (epsilon(infinity), D-infinity). Mixed conditions (sigma(infinity),D-infinity) and (epsilon(infinity),E-infinity) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E-y(infinity),sigma(y)(infinity)) for negative E-y(infinity) and (D-y(infinity), epsilon(y)(infinity)) for positive D-y(infinity) while the mechanical conditions sigma(y)(infinity) or epsilon(y)infinity are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth. (C) 2002 Published by Elsevier Science Ltd.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

Micromechanics analysis on evolution of crack in viscoelastic materials

A preliminary analysis on crack evolution in viscoelastic materials was presented. Based on the equivalent inclusion concept of micro-mechanics theory, the explicit expressions of crack opening displacement delta and energy release rate G were derived, indicating that both delta and G are increasing with time. The equivalent modulus of the viscoelastic solid comprising cracks was evaluated. It is proved that the decrease of the modulus comes from two mechanisms: one is the viscoelasticity of the material; the other is the crack opening which is getting larger with time.

Singularity of dynamic stress fields around an interface crack between viscoelastic bodies

The singular nature of the dynamic stress fields around an interface crack located between two dissimilar isotropic linearly viscoelastic bodies is studied. A harmonic load is imposed on the surfaces of the interface crack. The dynamic stress fields around the crack are obtained by solving a set of simultaneous singular integral equations in terms of the normal and tangent crack dislocation densities. The singularity of the dynamic stress fields near the crack tips is embodied in the fundamental solutions of the singular integral equations. The investigation of the fundamental solutions indicates that the singularity and oscillation indices of the stress fields are both dependent upon the material constants and the frequency of the harmonic load. This observation is different from the well-known -1/2 oscillating singularity for elastic bi-materials. The explanation for the differences between viscoelastic and elastic bi-materials can be given by the additional viscosity mismatch in the case of viscoelastic bi-materials. As an example, the standard linear solid model of a viscoelastic material is used. The effects of the frequency and the material constants (short-term modulus, long-term modulus and relaxation time) on the singularity and the oscillation indices are studied numerically.

Crack propagation in structures subjected to periodic excitation

In the present paper, a simple mechanical model is developed to predict the dynamic response of a cracked structure subjected to periodic excitation, which has been used to identify the physical mechanisms in leading the growth or arrest of cracking. The structure under consideration consists of a beam with a crack along the axis, and thus, the crack may open in Mode I and in the axial direction propagate when the beam vibrates. In this paper, the system is modeled as a cantilever beam lying on a partial elastic foundation, where the portion of the beam on the foundation represents the intact portion of the beam. Modal analysis is employed to obtain a closed form solution for the structural response. Crack propagation is studied by allowing the elastic foundation to shorten (mimicking crack growth) if a displacement criterion, based on the material toughness, is met. As the crack propagates, the structural model is updated using the new foundation length and the response continues. From this work, two mechanisms for crack arrest are identified. It is also shown that the crack propagation is strongly influenced by the transient response of the structure.

Effect of a laser pre-quenched steel substrate surface on the crack driving force in a coating-steel substrate system

A mechanical model of a coating/laser pre-quenched steel substrate specimen with a crack oriented perpendicular to the interface between the coating and the hardened layer is developed to quantify the effects of the residual stress and hardness gradient on the crack driving force in terms of the J-integral. It is assumed that the crack tip is in the middle of the hardened layer of the pre-quenched steel substrate. Using a composite double cantilever beam model, analytical solutions can be derived, and these can be used to quantify the effects of the residual stress and the hardness gradient resulting from the pre-quenched steel substrate surface on the crack driving force. A numerical example is presented to investigate how the residual compressive stress, the coefficient linking microhardness and yield strength and the Young's modulus ratio of the hardened layer to the coating influence the crack driving force for a given crack length. (C) 2007 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Atomistic simulations of crack nucleation and intergranular fracture in bulk nanocrystalline nickel

We report large scale molecular dynamics simulations of dynamic cyclic uniaxial tensile deformation of pure, fully dense nanocrystalline Ni, to reveal the crack initiation, and consequently intergranular fracture is the result of coalescence of nanovoids by breaking atomic bonds at grain boundaries and triple junctions. The results indicate that the brittle fracture behavior accounts for the transition from plastic deformation governed by dislocation to one that is grain-boundary dominant when the grain size reduces to the nanoscale. The grain-boundary mediated plasticity is also manifested by the new grain formation and growth induced by stress-assisted grain-boundary diffusion observed in this work. This work illustrates that grain-boundary decohesion is one of the fundamental deformation mechanisms in nanocrystalline Ni.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

On crack path stability in a layered material

Crack paths in an elastic layer on top of a substrate are considered. Crack growth is initiated from an edge crack in the layer. The plane of the initially straight crack forms an angle to the free surface. The load consists of a pair of forces applied at the crack mouth and parallel to the interface. Crack paths are calculated using a boundary element method. Crack growth is assumed to proceed along a path for which the mode II stress intensity factor vanishes. The inclination and the length of the initial crack are varied. The effect of two different substrates on the crack path evolution is demonstrated. A crack path initially leading perpendicularly to the interface is shown to be directionally unstable for a rigid substrate. Irrespective of its initial angle, the crack does not reach the interface, but reaches the free surface if the layer is infinitely long. At finite layer length the crack reaches the upper free surface if the initial crack inclination to the surface is small enough. For an inextendable flexible substrate, on the other hand, the crack reaches the interface if its initial inclination is large enough. For the flexible substrate an unstable path parallel with the sides of an infinitely long layer is identified. The results are compared with experimental results and discussed in view of characterisation of directionally unstable crack paths. The energy release rate for an inclined edge crack is determined analytically.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

Lattice instability at a fast moving crack tip

A molecular dynamics method is used to analyze the dynamic propagation of an atomistic crack tip. The simulation shows that the crack propagates at a relatively constant global velocity which is well below the Rayleigh wave velocity. However the local propagation velocity oscillates violently, and it is limited by the longitudinal wave velocity. The crack velocity oscillation is caused by a repeated process of crack tip blunting and sharpening. When the crack tip opening displacement exceeds a certain critical value, a lattice instability takes place and results in dislocation emissions from the crack tip. Based on this concept, a criterion for dislocation emission from a moving crack tip is proposed. The simulation also identifies the emitted dislocation as a source for microcrack nucleation. A simple method is used to examine this nucleation process. (C) 1996 American Institute of Physics.