Comparative study of time and frequency domain BEM approaches in frictional contact problem for antiplane crack under harmonic loading

Two different boundary element methods (BEM) for crack analysis in two dimensional (2-D) antiplane, homogeneous, isotropic and linear elastic solids by considering frictional contact of the crack edges are presented. Hypersingular boundary integral equations (BIE) in time-domain (TD) and frequency domain (FD), with corresponding elastodynamic fundamental solutions are applied for this purpose. For evaluation of the hypersingular integrals involved in BIEs a special regularization process that converts the hypersingular integrals to regular integrals is applied. Simple regular formulas for their calculation are presented. For the problems solution while considering frictional contact of the crack edges a special iterative algorithm of Udzava's type is elaborated and used. Numerical results for crack opening, frictional contact forces and dynamic stress intensity factors (SIFs) are presented and discussed for a finite III-mode crack in an infinite domain subjected to a harmonic crack-face loading and considering crack edges frictional contact interaction using the TD and FD approaches.     

The Time-Averaged Elastodynamic J-Integral

The time-averaged path independent J-integral for a stationary crack subjected to time-harmonic elastic waves is introduced. It can be determined from remote fields providing an alternative approach to compute the stress intensity factors. The J-integral is evaluated for a semi-infinite crack impinged by a plane sheer wave at an oblique angle.     

Scattering of SH-waves by an interface cavity

Summary. The scattering of the SH-wave and dynamic stress concentrations near an arbitrary cavity situated at the planar interface separating two different elastic media are investigated. The total wave field can be obtained by superposition of the free field and the scattered field. The free field is composed of the incident, reflected and refracted waves. The scattered wave fields in adjacent media are expressed respectively, and the method of wave functions expansion is applied to obtain the solutions for these fields. The scattered wave functions can be expanded into Hankel-Fourier series with unknown coefficients. In solving for the unknown coefficients according to the boundary conditions for the total wave field at the interface and at the cavity wall, the non-orthogonality makes the system of equations for the unknown coefficients infinite and coupling each other. Another key point is to extend each scattered wave field from its own half-plane domain into the full plane domain by a certain way keeping the total wave field unchanged for the non-orthogonal Fourier integrals around the cavity. Finally, the scattering of the SH wave by an interface ellipse with different ratios between long and short axis is considered, and the distributions of dynamic stress concentration factors at the cavity wall are presented.     

Interaction integrals for fracture analysis of functionally graded magnetoelectroelastic materials

This paper presents the domain form of interaction integrals based on three independent formulations for computation of stress intensity factors, electric displacement intensity factors and magnetic induction intensity factors for cracks in functionally graded magnetoelectroelastic materials. Conservation integrals of J-type are derived based on the governing equations for magnetoelectroelastic media and the crack tip asymptotic fields of homogeneous magnetoelectroelastic medium as auxiliary fields. Each of the formulations differs in the way auxiliary fields are imposed in the evaluation of interaction integrals and each of them results in a consistent form of the interaction integral in the sense that extra terms naturally appear in their derivation to compensate for the difference in the chosen crack tip asymptotic fields of homogeneous and functionally graded magnetoelectroelastic medium. The additional terms play an important role of ensuring domain independence of the presented interaction integrals. Comparison of numerically evaluated intensity factors through the three consistent formulations with those obtained using displacement extrapolation method is presented by means of two examples.     

An Elliptical Crack in an Anisotropic Elasic Medium Subjected to a Constant External Field

An isolated elliptical crack in an infinite orthotropic elastic medium is considered. An efficient numerical algorithm of the solution of the problem for a crack subjected to a constant external field is proposed. The calculation of the crack opening vector and the stress intensity factors on the crack edge is reduced to regular 2D-integrals. These integrals may be simply calculated numerically for an arbitrary orientation of the crack plane with respect to the principal axes of the anisotropy of the medium. Examples of the calculation of the crack opening vector and stress intensity factors are presented.     

Transient analysis of stress waves around two rectangular cracks under impact load

The three-dimensional response of two rectangular cracks in an infinite elastic medium to impact load is investigated in this paper. Fourier and Laplace transforms are applied and the problem is reduced to that of solving dual integral equations in the Laplace transform domain. To solve these equations, the crack surface displacement is expanded in a double series of functions which are zero outside of the cracks. The unknown coefficients accompanied in that series are solved with the aid of the Schmidt method. The dynamic stress intensity factors are computed numerically.     

The gauss-hermite numerical integration method for the solution of the plane elastic problem of semi-infinite periodic cracks

The problem of parallel semi-infinite periodic cracks subjected to a transversely directed load in an infinite isotropic and elastic medium under conditions of plane stress or plane strain can be reduced to the solution of a Cauchy-type singular integral equation along one of the cracks. This equation can be transformed into a system of linear equations by means of an approximation of the integrals through the Gauss-Hermite procedure and application of the equation to distinct points along the faces of the crack. Stress intensity factors thus determined for the crack tips under constant load along the cracks are in satisfactory agreement with corresponding values derived previously.     

Effect of stress wave propagation phenomenon on the determination of strain energy density theory parameters and dynamic J-integral

Dynamic loading for stationary cracks leads to results that are many times greater in magnitude than their static counterparts. If the dynamic loading is in the form of impact type, stress wave propagation effects become dominant. FRAC3D program comprises enriched element formulation which doesn't require excessive mesh refinement around crack tip for accuracy. Strain energy density (SED) theory parameters and dynamic J-integral are sought in this study to simulate and understand wave propagation phenomenon in detail. Structures under the effect of wave propagations yield more reliable J-integral values by taking the average of the results from multiple domain sizes. Governed by stress waves, space-time variations of minimum energy density locations strongly influence fracture characterization for straight and curved crack fronts. Details given in numerical examples section of this paper make a great contribution to understanding of the response for cracked structures subjected to sudden loading.     

Trifurcation of a crack due to plane SH-waves in an infinite elastic medium

The dynamic anti-plane problem of trifurcation of a semi-infinite crack due to incidence of two linearly varying plane SH-waves with non-parallel wave fronts in an infinite elastic medium has been considered. The semi-infinite crack is assumed to trifurcate when the plane waves intersect the crack tip. The problem has been solved using the self-similar technique, which is based on the observation that certain field variables show dynamic similarities. The results include the expressions for shear stress in the planes of the cracks and the stress intensity factors at the crack tips. Numerical calculations have been carried out to show the variations of stress intensity factors at the crack tips with the angle of skew for different values of the crack tip velocity and angle of incidence.     

Transient response of a surface crack subjected to dynamic anti-plane concentrated loadings

In this study, the transient response of a surface crack in an elastic solid subjected to dynamic anti-plane concentrated loadings is investigated. The angles of the surface crack and the half-plane are 60° and 90°. In analyzing this problem, an infinite number of diffracted and reflected waves generated by the crack tip and edge boundaries must be taken into account and it will make the analysis extremely difficult. The solutions are determined by superposition of the proposed fundamental solution in the Laplace transform domain and by using the method of image. The fundamental solution to be used is the problem for applying exponentially distributed traction on the crack faces. The exact transient solutions of dynamic stress intensity factor are obtained and expressed in formulations of series form. The solutions are valid for an infinite length of time and have accounted for the contribution of an infinite number of diffracted waves. The explicit value of the dynamic overshot for the perpendicular surface crack is obtained from the analysis. Numerical results are evaluated which indicate that the dynamic stress intensity factors will oscillate near the correspondent static values after the first three or six waves have passed the crack tip.     

Extension of a crack due to plane SH wave in a pre-stressed infinite elastic medium

The problem of extension of an infinitesimal flaw into a plane crack due to two linearly varying plane SH-waves with non-parallel wave fronts in an infinite elastic medium which is initially in a state of uniform anti-plane shear, has been considered. Fracture is assumed to initiate at a point a finite time after the waves intersect there. The crack is assumed to extend non-symmetrically along the trace of the wave intersection. The method of analysis is based on the observation that certain field quantities show dynamic similarity. The results include expressions for the stress intensity factors at the crack tips and the rate of energy flux into the crack edges. Numerical calculations are carried out to obtain stress intensity factors and the rate of energy flux into the crack tips for different values of the parameters.     

A comparative study of three BEM for transient dynamic crack analysis of 2-D anisotropic solids

Three different boundary element methods (BEM) for transient dynamic crack analysis in two-dimensional (2-D), homogeneous, anisotropic and linear elastic solids are presented. Hypersingular traction boundary integral equations (BIEs) in frequency- domain, Laplace-domain and time-domain with the corresponding elastodynamic fundamental solutions are applied for this purpose. In the frequency-domain and the Laplace-domain BEM, numerical solutions are first obtained in the transformed domain for discrete frequency or Laplace-transform parameters. Time-dependent results are subsequently obtained by means of the inverse Fourier-transform and the inverse Laplace-transform algorithm of Stehfest. In the time-domain BEM, the quadrature formula of Lubich is adopted to approximate the arising convolution integrals in the time-domain BIEs. Hypersingular integrals involved in the traction BIEs are computed through a regularization process that converts the hypersingular integrals to regular integrals, which can be computed numerically, and singular integrals which can be integrated analytically. Numerical results for the dynamic stress intensity factors are presented and discussed for a finite crack in an infinite domain subjected to an impact crack-face loading.     

Transient analysis of dynamic crack propagation in piezoelectric materials

Abstract

In this paper, the transient analysis of semi‐infinite propagating cracks in piezoelectric materials subjected to dynamic anti‐plane concentrated body force is investigated. The crack surface is assumed to be covered with an infinitesimally thin, perfectly conducting electrode that is grounded. In analyzing this problem, it has characteristic lengths and a direct attempt towards solving this problem by transform and Wiener‐Hopf techniques (Noble, 1958) is not applicable. In order to solve this problem, a new fundamental solution for propagating cracks in piezoelectric materials is first established and the transient response of the propagating crack is obtained by superposition of the fundamental solution in the Laplace transform domain. The fundamental solution to be used is the responses of applying exponentially distributed traction in the Laplace transform domain on the propagating crack surface. Taking into account the quasi‐static approximation, exact analytical transient solutions for the dynamic stress intensity factor and the dynamic electric displacement intensity factor are obtained by using the Cagniard‐de Hoop method (Cagnard, 1939; de Hoop, 1960) of Laplace inversion and are expressed in explicit forms. Numerical calculations of dynamic intensity factors are evaluated and the results are discussed in detail. The transient solutions for stationary cracks have been shown to approach the corresponding static values after the shear wave of the piezoelectric material has passed the crack tip.     

An analytical expression for the J2‐integral of an interfacial crack in orthotropic bimaterials

This paper presents a new analytical expression relating the J₂‐integral and stress intensity factors (SIF) in an in‐plane traction‐free crack between two orthotropic elastic solids using the complex function method. The singular oscillatory near tip field of a bimaterial interfacial crack is usually characterized by a pair of SIFs. In linear elastic interfacial fracture mechanics, the majority of numerical and experimental methods rely on the analytical equations relating J_k‐integrals and SIFs. Although an analytical equation relating J₁‐integral or strain energy release rate and SIFs is available, a similar relation for J₂‐integral in debonded anisotropic solids is non‐existent. Using this new analytical expression, in conjunction with the values of J_k, the SIFs can be computed without the need for an auxiliary relation. An example with known analytical solutions for SIFs is presented to show the variation of the J₂‐integral near the crack tip of a bimaterial orthotropic plate. Different bimaterial combinations are considered, and the effect of material mismatch on J_k is demonstrated.     

Effect of Couple-Stresses on the Transient Dynamic Stress Intensity Factors for a Crack in an Infinite Elastic Medium Under an Impact Stress Wave

This study applies linearized couple-stress theory to evaluate the dynamic stresses around a crack in an infinite elastic medium that is subjected to an incoming shock stress wave impinging normal to the crack. The boundary conditions with respect to the crack are reduced to dual integral equations using a Fourier transform in the Laplace domain. To solve these equations, the differences in the displacement and rotation at the crack are expanded by a series of functions that are zero-valued outside the crack in the Laplace domain.     

The interaction of two inclined cracks with dynamic stress wave loadings

To gain insight into the phenomenon of the interaction of stress waves with material defects and the linkage of two cracks, the transient response of two semi-infinite inclined cracks subjected to dynamic loading is examined. The solutions are obtained by the linear superposition of fundamental solutions in the Laplace transform domain. The fundamental solution is the exponentially distributed traction on crack faces proposed by Tsai and Ma [1]. The exact closed form solutions of stress intensity factor histories for these two inclined cracks subjected to incident plane waves and diffracted waves are obtained explicitly. These solutions are valid for the time interval from initial loading until the first wave scattered at one crack tip returns to the same crack tip after being diffracted by another crack tip. The result shows that the contribution of diffracted waves to stress intensity factors is much less than the incident waves. The probable crack propagation direction is predicted from the fracture criterion of maximum circumferential tensile stress. The linkage of these two cracks is also investigated in detail.     

Theoretical simulations of a propagating crack subjected to in-plane stress wave loading by caustic method

The optical method of caustics for measuring the dynamic stress intensity factor in a transient process is investigated in this study. The transient full-field solutions of a propagating crack contained in an infinite medium subjected to step-stress wave and ramp-stress wave loadings are used to establish the exact equations of the initial and caustic curves. The results of the stress intensity factor obtained from the caustic method are compared with theoretical predictions and some experiments. The results demonstrate that a significant deviation can occur in the determination of the dynamic stress intensity factor from shadow spot measurements. The factors, such as screen distance, magnitude of loading, crack speed and rising time which can influence the accuracy of the experimental measurements are discussed in detail. In addition, the valid region of the dynamic stress singular field for the propagating crack is discussed in detail and it gives a better understanding of the appropriate region of measurements for investigators. This revised version was published online in July 2006 with corrections to the Cover Date.     

Transient response of a finite crack subjected to dynamic anti-plane loading

In this study, the transient response of a finite crack in an elastic solid subjected to dynamic antiplane loading is investigated. Two specific loading situations, a body force near the finite crack and a concentrated point loading applied on the crack face, are analyzed in detail. In analyzing this problem, an infinite number of diffracted waves generated by two crack tips must be taken into account which will make the analysis extremely difficult. The solutions are determined by superposition of proposed fundamental solutions in the Laplace transform domain.The fundamental solutions to be used are the problems for applying exponentially distributed traction and screw dislocation to the crack faces and along the crack-tip line respectively. Exact transient closed-form solutions for the dynamic stress intensity factor are obtained and expressed in very simple and compact formulations. The solutions are valid for an infinite length of time and have accounted for the contributions of an infinite number of diffracted waves. Numerical calculations for the two problems are evaluated and results indicate that the dynamic stress intensity factors will oscillate near the corresponding static values after the first three waves have passed through the specified crack tip.     

Perfectly matched layers for transient elastodynamics of unbounded domains

One approach to the numerical solution of a wave equation on an unbounded domain uses a bounded domain surrounded by an absorbing boundary or layer that absorbs waves propagating outward from the bounded domain. A perfectly matched layer (PML) is an unphysical absorbing layer model for linear wave equations that absorbs, almost perfectly, outgoing waves of all non‐tangential angles‐of‐incidence and of all non‐zero frequencies. In a recent work [Computer Methods in Applied Mechanics and Engineering 2003; 192: 1337–1375], the authors presented, inter alia, time‐harmonic governing equations of PMLs for anti‐plane and for plane‐strain motion of (visco‐) elastic media. This paper presents (a) corresponding time‐domain, displacement‐based governing equations of these PMLs and (b) displacement‐based finite element implementations of these equations, suitable for direct transient analysis. The finite element implementation of the anti‐plane PML is found to be symmetric, whereas that of the plane‐strain PML is not. Numerical results are presented for the anti‐plane motion of a semi‐infinite layer on a rigid base, and for the classical soil–structure interaction problems of a rigid strip‐footing on (i) a half‐plane, (ii) a layer on a half‐plane, and (iii) a layer on a rigid base. These results demonstrate the high accuracy achievable by PML models even with small bounded domains. Copyright © 2004 John Wiley & Sons, Ltd.