Simple formula for dynamic stress intensity factor of pre-cracked Charpy specimen

It has been reported that consideration of inertia effect is important to estimate the stress intensity factor for impact loading. In the present paper a simple expression for the dynamic stress intensity factor of a pre-cracked Charpy specimen is derived by making use of an approximate solution for one-dimensional beam equations. The time variation of the dynamic stress intensity factor obtained retaining only the fundamental mode of bending vibration is in good agreement with the two-dimensional finite element solutions for step-loading and also for final-peak-sawtooth-pulse loading.     

Stress intensity factor analyses of three-dimensional interface cracks using tetrahedral finite elements

A numerical method is proposed for evaluating the stress intensity factors of a three-dimensional bimaterial interfacial crack using tetrahedral finite elements. This technique is based on the M₁-integral method and employs the moving least-squares approximation. Stress or strain in the M₁-integral equation is automatically approximated from the nodal displacements obtained by the finite element analysis using the moving least-squares method. Therefore, the presented method needs no elemental information from the finite element analysis. In this study, stress intensity factor analyses of some typical three-dimensional interface crack problems using the tetrahedral finite elements are demonstrated.     

A Method to Extract Interface Stress Intensity Factors for Bimaterial Problems using Interlayers

A numerical method is presented here to determine stress intensity factors for interface cracks in plane, isotropic, elastic bimaterial fracture problems. The method relies on considering a companion problem wherein a very thin elastic interlayer with a crack, is artificially inserted between the two material regions of the original bimaterial problem. Modes I and II stress intensity factors are obtained for the companion problem using the modified virtual crack closure method. These stress intensity factors for the companion problem are then converted to the stress intensity factors for the original interface crack problem with the help of a universal relation. This universal relation between the stress intensity factors of the two problems is established by considering an asymptotic problem where the thickness of the interlayer is small compared with all other length scales. Two benchmark problems are considered to demonstrate the effectiveness of the interlayer approach for determining interface stress intensity factors.     

A complex variable boundary element method for solving interface crack problems

The problem of finite bimaterial plates with an edge crack along the interface is studied. A complex variable boundary element method is presented and applied to determine the stress intensity factor for finite bimaterial plates. Using the pseudo-orthogonal characteristic of the eigenfunction expansion forms and the well-known Bueckner work conjugate integral and taking the different complex potentials as auxiliary fields, the interfacial stress intensity factors associated with the physical stress-displacement fields are evaluated. The effects of material properties and crack geometry on stress intensity factors are investigated. The numerical examples for three typical specimens with six different combinations of the bimaterial are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

Diffraction of horizontal shear waves by a moving interface crack

Summary An analysis of the diffraction of horizontally polarized shear waves by a finite crack moving on a bimaterial interface is carried out. Fourier transform method is used to reduce the mixed boundary value problem to the solutions of two pairs of dual integral equations. These equations are further reduced to a pair of Fredholm integral equations of the second kind. The dynamic stress intensity factors are obtained for several values of wave number, incident angle, crack velocity, and material constants.With 7 Figures     

Initiation, propagation and arrest of an interface crack subjected to controlled stress wave loading

An experimental study has been conducted to investigate the initiation, propagation, and arrest of bimaterial interface cracks subjected to controlled stress wave loading in the form of a tensile dilatational stress wave pulse. The tensile pulse is generated by detonating lead azide explosive in a specially designed specimen. Dynamic loading of the bimaterial interface results in crack initiation, propagation, and arrest, all in the same experiment. This failure event is observed using photoelasticity in conjunction with high speed photography. Full field data from the experimentally obtained isochromatic fringe patterns is analyzed to determine time histories of various fracture parameters such as the crack tip speed, the dynamic complex stress intensity factor, the energy release rate, and the mixity. The experimental data is also used to quantify the values of the dynamic initiation and arrest toughness and to evaluate a recently proposed dynamic interface fracture criterion. This revised version was published online in August 2006 with corrections to the Cover Date.     

Analysis of surface cracks using the line-spring boundary element method and the virtual crack extension technique

The authors have developed a new line-spring boundary element method which couples the line-spring model with the boundary element method to deal with the problem of a surface cracked plate. However, the drawback of the line-spring model is that a reliable stress intensity factor could not be directly obtained near the free surface intersection. Therefore, the virtual crack extension technique is employed in a post-processor of the line-spring boundary element method to obtain the stress intensity factor at the crack front-free surface intersection. Theoretical analysis is described. Stress intensity factors for surface cracks are calculated to verify the proposed method. The interaction of two surface cracks is also investigated. The solutions obtained by the line-spring boundary element method show that the method proposed is efficient and reasonably accurate.     

Conducting cracks in dissimilar piezoelectric media

Complete stress and electric fields near the tip of a conducting crack between two dissimilar anisotropic piezoelectric media, are obtained in terms of two generalized bimaterial matrices proposed in this paper. It is shown that the general interfacial crack-tip field consists of two pairs of oscillatory singularities. New definitions of real-valued stress and electric field intensity factors are proposed. Exact solutions of the stress and electric fields for basic interface crack problems are obtained. An alternate form of the J integral is derived, and the mutual integral associated with the J integral is proposed. Closed form solutions of the stress and electric field intensity factors due to electromechanical loading and the singularities for a semi-infinite crack as well as for a finite crack at the interface between two dissimilar piezoelectric media, are also obtained by using the mutual integral.     

Finite element calculation of stress intensity factors for interfacial crack using virtual crack closure integral

This paper presents a successful implementation of the virtual crack closure integral method to calculate the stress intensity factors of an interfacial crack. The present method would compute the mixed-mode stress intensity factors from the mixed-mode energy release rates of the interfacial crack, which are easily obtained from the crack opening displacements and the nodal forces at and ahead of the crack tip, in a finite element model. The simple formulae which relate the stress intensity factors to the energy release rates are given in three separate categories: an isotropic bimaterial continuum, an orthotropic bimaterial continuum, and an anisotropic bimaterial continuum. In the example of a central crack in a bimaterial block under the plane strain condition, comparisons are made with the exact solution to determine the accuracy and efficiency of the numerical method. It was found that the virtual crack closure integral method does lead to very accurate results with a relatively coarse finite element mesh. It has also been shown that for an anisotropic interfacial crack under the generalized plane strain condition, the computed stress intensity factors using the virtual crack closure method compared favorably with the results using the J integral method applied to two interacting crack tip solutions. In order for the stress intensity factors to be used as physical variables, the characteristic length for the stress intensity factors must be properly defined. A study was carried out to determine the effects of the characteristic length on the fracture criterion based the mixed-mode stress intensity factors. It was found that the fracture criterion based on the quadratic mixture of the normalized stress intensity factors is less sensitive to the changes in characteristic length than the fracture criterion based on the total energy release rate along with the phase angle.This work has been supported by ONR, with Dr. Y. Rajapakse as the program official.     

Line-spring model for surface cracks in a reissner plate

In this paper the line-spring model developed by Rice and Levy for a surface crack in elastic plates is reconsidered. The problem is formulated by using Reissner's plate bending theory. For the plane strain problem of a strip containing an edge crack and subjected to tension and bending new expressions for stress intensity factors are used which are valid up to a depth-to-thickness ratio of 0.8. The stress intensity factors for a semi-elliptic and a rectangular crack are calculated. Considering the simplicity of the technique and the severity of the underlying assumptions, the results compare rather well with the existing finite element solutions.     

A definition and evaluation procedure of generalized stress intensity factors at cracks and multi-material wedges

This paper proposes a definition of generalized stress intensity factors that includes classical definitions for crack problems as special cases. Based on the semi-analytical solution obtained from the scaled boundary finite-element method, the singular stress field is expressed as a matrix power function with its dimension equal to the number of singular terms. Not only real and complex power singularities but also power-logarithmic singularities are represented in a unified expression without explicitly determining the type of singularity. The generalized stress intensity factors are evaluated directly from the scaled boundary finite-element solution for the singular stress field by following standard stress recovery procedures in the finite element method. The definition and evaluation procedure are valid to multi-material wedges composed of any number of isotropic and anisotropic materials. Numerical examples, including a cracked homogeneous plate, a bimaterial plate with an interfacial crack, a V-notched bimaterial plate and a crack terminating at a material interface, are analyzed. Features of this unified definition are discussed.     

Elastodynamics of a crack on the bimaterial interface

The paper is an application of boundary integral equations to the problem of a crack located on the bimaterial interface under time-harmonic loading. A system of linear algebraic equations is derived for solving the problem numerically. The distributions of the displacements and tractions at the bimaterial interface are obtained and analysed for the case of a penny-shaped crack under normal tension-compression wave. The dynamic stress intensity factors (normal and shear modes) are also computed. The results are compared with those obtained for the static case.     

Partition of unity enrichment for bimaterial interface cracks

Partition of unity enrichment techniques are developed for bimaterial interface cracks. A discontinuous function and the two‐dimensional near‐tip asymptotic displacement functions are added to the finite element approximation using the framework of partition of unity. This enables the domain to be modelled by finite elements without explicitly meshing the crack surfaces. The crack‐tip enrichment functions are chosen as those that span the asymptotic displacement fields for an interfacial crack. The concept of partition of unity facilitates the incorporation of the oscillatory nature of the singularity within a conforming finite element approximation. The mixed‐mode (complex) stress intensity factors for bimaterial interfacial cracks are numerically evaluated using the domain form of the interaction integral. Good agreement between the numerical results and the reference solutions for benchmark interfacial crack problems is realized. Copyright © 2004 John Wiley & Sons, Ltd.     

Determination of thermal stress intensity factors for an interface crack under vertical uniform heat flow

In the case where an interface crack exists in an infinite two-dimensional elastic bimaterial, the crack surface is insulated under traction-free conditions and the uniform heat flow vertical to the crack from an infinite boundary is given, temperature and stress potentials are obtained by using the complex variable approach to solve Hubert problems, and the results are used to obtain thermal stress intensity factors. The mode II thermal stress intensity factor only occurs if both the shear moduli, as well as the Poisson's ratios in the upper and lower material, are the same. Otherwise, mode I and II thermal stress intensity factors exist but the value of the mode I thermal stress intensity factor is much smaller than that of mode II.     

Kinking of an interface crack in an orthotropic bimaterial

We investigated the asymptotic problem of a kinked interface crack in an orthotropic bimaterial under in‐plane loading conditions. The stress intensity factors at the tip of the kinked interface crack are described in terms of the stress intensity factors of the interface crack prior to the kink combined with a dimensionless matrix function. Using a modified Stroh formalism and an orthotropy rescaling technique, the matrix function was obtained from the solutions of the corresponding problem in transformed bimaterial. The effects of orthotropic and bimaterial parameters on the matrix function were examined. A reduction in the number of dependent material parameters on the matrix function was made using the modified Stroh formalism. Moreover, the explicit dependence of one orthotropic parameter on the matrix function was determined using an orthotropic rescaling technique. The effects of the other material parameters on the matrix function were numerically examined. The energy release rate was obtained for a kinked interface crack in an orthotropic bimaterial.     

A methodology for measuring interface fracture properties of composite materials

A methodology is presented for measuring interface fracture properties of composite materials. A bimaterial Brazilian disk specimen with a crack along the interface is employed. The specimen is analyzed by means of the finite element method and a conservative integral to determine stress intensity factors as a function of loading angle and crack length. A weight function is employed to determine the effect of residual curing stresses on the stress intensity factors. These are combined to determine the critical interface energy release rate _ic as a function of stress intensity phase angle for tests carried out on a glass/epoxy material pair.     

Weight function calculation for surface cracks using the line-spring model

A numerical method for calculating weight functions for surface cracks in plates and shells is proposed. Thick-shell finite elements are used to create the discrete model of a body with a through-wall flaw. Line-spring elements transform the through-wall flaw into a surface crack. A quadratic line-spring element is presented. Weight functions for some semielliptical surface cracks in a plate have been calculated. The weight functions obtained may be used for computing stress intensity factors related to two-dimensional stress fields at the crack surface.     

Fracture mechanics parameters for cracks on a slightly undulating interface

Typical bimaterial interfaces are non-planar due to surface facets or roughness. Crack-tip stress fields of an interface crack must be influenced by non-planarity of the interface. Consequently, interface toughness is affected. In this paper, the crack-tip fields of a finite crack on an elastic/rigid interface with periodic undulation are studied. Particular emphasis is given to the fracture mechanics parameters, such as the stress intensity factors, crack-tip energy release rate, and crack-tip mode mixity. When the amplitude of interface undulation is very small relative to the crack length (which is the case for rough interfaces), asymptotic analysis is used to convert the non-planarity effects into distributed dislocations located on the planar interface. Then, the resulting stress fields near the crack tip are obtained by using the Fourier integral transform method. It is found that the stress fields at the crack tip are strongly influenced by non-planarity of the interface. Generally speaking, non-planarity of the interface tends to shield the crack tip by reducing the crack-tip stress concentration.     

An interaction integral method for computing mixed mode stress intensity factors for curved bimaterial interface cracks in non-uniform temperature fields

In finite element analysis the interaction integral has been a useful tool for computing the stress intensity factors for fracture analysis. This work extends the interaction integral to account for non-uniform temperatures in the calculation of stress intensity factors for three dimensional curvilinear cracks either in a homogeneous body or on a bimaterial interface. First, the derivation of the computational algorithm, which includes the additional terms developed by the non-zero gradient of the temperature field, is presented in detail. The algorithm is then implemented in conjunction with commercial finite element software to calculate the stress intensity factors of a crack undergoing non-uniform temperatures on both a homogeneous and a bimaterial interface. The numerical results displayed path independence and showed excellent agreement with available analytical solutions.