Transient response of a Griffith crack between dissimilar piezoelectric layers under anti-plane mechanical and in-plane electrical impacts

The problem of an interface crack between dissimilar piezoelectric layers under mechanical and electrical impacts is formulated by using integral transform and Cauchy singular integral equation methods. The dynamic stress intensity factor and dynamic energy release rate (DERR) are determined through use of the obtained solutions and the effects of the loading ratio, the geometry of crack configuration and the combination of material parameters on the above two quantities are discussed. The numerical calculations indicate that the electrical load can promote or retard the crack growth depending on its magnitude, direction and the existence of the mechanical load and that with the increase of the value of ratio of two material parameters, some material parameters will inhibit the crack growth. On the other hand, some material parameters play the contrary roles. In addition, the geometry of the crack configuration has the significant effects on the DERR. Finally the results are compared with those obtained in a previous investigation.     

Dynamic analysis of interfacial crack problems in anisotropic bi-materials by a time-domain BEM

A time-domain boundary element method (BEM) together with the sub-domain technique is applied to study dynamic interfacial crack problems in two-dimensional (2D), piecewise homogeneous, anisotropic and linear elastic bi-materials. The bi-material system is divided into two homogeneous sub-domains along the interface and the traditional displacement boundary integral equations (BIEs) are applied on the boundary of each sub-domain. The present time-domain BEM uses a quadrature formula for the temporal discretization to approximate the convolution integrals and a collocation method for the spatial discretization. Quadratic quarter-point elements are implemented at the tips of the interface cracks. A displacement extrapolation technique is used to determine the complex dynamic stress intensity factors (SIFs). Numerical examples for computing the complex dynamic SIFs are presented and discussed to demonstrate the accuracy and the efficiency of the present time-domain BEM.     

Dynamic crack propagation in matrix involving inclusions by a time-domain BEM

The dynamic crack propagation in composites with inclusions is analyzed by a time-domain boundary element method (BEM) in conjunction with the sub-region technique. The crack-growth direction and the crack-tip instantaneous velocity are determined by the maximum circumferential stress criterion. The instantaneous velocity is well smoothed by a bisection technique. New crack-tip elements of inconstant length are added to the active crack-tip to simulate the fast crack growth. Square root shape functions are adopted as to describe the proper asymptotic behavior in the vicinity of the crack-tips. The computation time for the dynamic fracture problems in composites with multiple inclusions is reduced by a numerical method. The influences of the inclusions on the shielding ratio, the crack growing path and the crack-tip instantaneous speed are well investigated.     

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.     

Material properties of piezoelectric composites by BEM and homogenization method

Applications of boundary element method (BEM) to piezoelectric composites in conjunction with homogenization approach for determining their effective material properties are discussed in this paper. The composites considered here consist of inclusion and matrix phases. The homogenization model for composites with inhomogeneities is developed and introduced into a BE formulation to provide an effective means for estimating overall material constants of two-phase composites. In this model, a representative volume element (RVE) is used whose volume average stress and strain are calculated by the boundary tractions and displacements of the RVE. Thus BEM is suitable for performing calculations on average stress and strain fields of the composites. Numerical results for a piezoelectric plate with circular inclusions are presented to illustrate the application of the proposed micromechanics––BE formulation.     

Transient dynamic crack analysis in piecewise homogeneous, anisotropic and linear elastic composites by a spatial symmetric Galerkin-BEM

Transient elastodynamic analysis of two-dimensional, piecewise homogeneous, anisotropic and linear elastic solids containing interior and interface cracks is presented in this paper. To solve the initial boundary value problem, a spatial symmetric time-domain boundary element method is developed. Stationary cracks subjected to impact loading conditions are considered. Elastodynamic fundamental solutions for homogenous, anisotropic and linear elastic solids are implemented. The piecewise homogeneous, anisotropic and linear elastic solids are modeled by the multi-domain technique. The spatial discretization is performed by a symmetric Galerkin-method, while a collocation method is utilized for the temporal discretization. An explicit time-stepping scheme is obtained for computing the unknown boundary data. Numerical examples are presented and discussed to show the effects of the interface cracks, the material anisotropy, the material combination and the dynamic loading on the dynamic stress intensity factors.     

A two-dimensional time-domain boundary element method for dynamic crack problems in anisotropic solids

A time-domain boundary element method (BEM) for transient dynamic crack analysis in two-dimensional, homogeneous, anisotropic and linear elastic solids is presented in this paper. Strongly singular displacement boundary integral equations (DBIEs) are applied on the external boundary of the cracked body while hypersingular traction boundary integral equations (TBIEs) are used on the crack-faces. The present time-domain method uses the quadrature formula of Lubich for approximating the convolution integrals and a collocation method for the spatial discretization of the time-domain boundary integral equations. Strongly singular and hypersingular integrals are dealt with by a regularization technique based on a suitable variable change. Discontinuous quadratic quarter-point elements are implemented at the crack-tips to capture the local square-root-behavior of the crack-opening-displacements properly. Numerical examples for computing the dynamic stress intensity factors are presented and discussed to demonstrate the accuracy and the efficiency of the present method.     

A hypersingular time‐domain BEM for 2D dynamic crack analysis in anisotropic solids

A hypersingular time‐domain boundary element method (BEM) for transient elastodynamic crack analysis in two‐dimensional (2D), homogeneous, anisotropic, and linear elastic solids is presented in this paper. Stationary cracks in both infinite and finite anisotropic solids under impact loading are investigated. On the external boundary of the cracked solid the classical displacement boundary integral equations (BIEs) are used, while the hypersingular traction BIEs are applied to the crack‐faces. The temporal discretization is performed by a collocation method, while a Galerkin method is implemented for the spatial discretization. Both temporal and spatial integrations are carried out analytically. Special analytical techniques are developed to directly compute strongly singular and hypersingular integrals. Only the line integrals over an unit circle arising in the elastodynamic fundamental solutions need to be computed numerically by standard Gaussian quadrature. An explicit time‐stepping scheme is obtained to compute the unknown boundary data including the crack‐opening‐displacements (CODs). Special crack‐tip elements are adopted to ensure a direct and an accurate computation of the elastodynamic stress intensity factors from the CODs. Several numerical examples are given to show the accuracy and the efficiency of the present hypersingular time‐domain BEM. Copyright © 2008 John Wiley & Sons, Ltd.     

Electroelastic analysis of an interface anti-plane shear crack in a layered piezoelectric plate

The problem of an anti-plane interface crack in a layered piezoelectric plate composed of two bonded dissimilar piezoelectric ceramic layers subjected to applied voltage is considered. It is assumed that the crack is either impermeable or permeable. An integral transform technique is employed to reduce the problem considered to dual integral equations, then to a Fredholm integral equation by introducing an auxiliary function. Field intensity factors and energy release rate are obtained in explicit form in terms of the auxiliary function. In particular, by solving analytically a resulting singular integral equation, they are determined explicitly in terms of given electromechanical loadings for the case of two bonded layers of equal thickness. Some numerical results are presented graphically to show the influence of the geometric parameters on the field intensity factors and the energy release rate.     

Transient analysis of a piezoelectric strip with a permeable crack under anti-plane impact loads

The problem of a through permeable crack situated in the mid-plane of a piezoelectric strip is considered under anti-plane impact loads for two cases. The first is that the strip boundaries are free of stresses and of electric displacements, and the second is that the strip boundaries are clamped rigid electrodes. The method adopted is to reduce the mixed initial-boundary value problem, by using integral transform techniques, to dual integral equations, which are further transformed into a Fredholm integral equation of the second kind by introducing an auxiliary function. The dynamic stress intensity factor and energy release rate in the Laplace transform domain are obtained in explicit form in terms of the auxiliary function. Some numerical results for the dynamic stress intensity factor are presented graphically in the physical space by using numerical techniques for solving the resulting Fredholm integral equation and inverting Laplace transform.     

Stress intensity factor analysis of a three-dimensional interface crack between dissimilar anisotropic materials

A new numerical method to calculate the stress intensity factors (SIFs) of a three-dimensional interface crack between dissimilar anisotropic materials was developed. In this study, the M-integral method was employed for mode separation of the SIFs. The moving least-square method was utilized to calculate the M-integral. Using the M-integral with the moving least-square method, SIFs can be automatically calculated with only the nodal displacements from the finite element method (FEM). Here, SIFs analyses of some typical three-dimensional problems are demonstrated. Excellent agreement was achieved between the numerical results obtained by the present method and the corresponding results proposed by other researchers. In addition, the SIFs of a single-edge crack, a through crack, and a semi-circular crack between two anisotropic solids in three-dimensional structures were analyzed.     

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.     

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.     

Dynamic fracture behavior of an internal interfacial crack between two dissimilar magneto-electro-elastic plates

In this paper, the anti-plane problem for an interfacial crack between two dissimilar magneto-electro-elastic plates subjected to anti-plane mechanical and in-plane magneto-electrical impact loadings is investigated. Four kinds of crack surface conditions are adopted: magneto-electrically impermeable (Case 1), magnetically impermeable and electrically permeable (Case 2), magnetically permeable and electrically impermeable (Case 3), and magneto-electrically permeable (Case 4). The position of the interfacial crack is arbitrary. The Laplace transform and finite Fourier transform techniques are employed to reduce the mixed boundary-value problem to triple trigonometric series equations in the Laplace transform domain. Then the dislocation density functions and proper replacements of the variables are introduced to reduce the series equations to a standard Cauchy singular integral equation of the first kind. The resulting integral equation together with the corresponding single-valued condition is approximated as a system of linear algebra equations, which can easily be solved. Field intensity factors and energy release rates are determined and discussed. The effects of loading combination parameters on dynamic energy release rate are plotted for Cases 1-3. On the other hand, since the magneto-electrically permeable condition is perhaps more physically reasonable for type III crack, the effect of the crack configuration on the dynamic fracture behavior of the crack tips is studied in detail for Case 4. The results could be useful for the design of multilayered magneto-electro-elastic structures and devices.     

Transient response of a cracked magnetoelectric material under the action of in-plane sudden impacts

Dynamic analysis of a crack embedded in a magnetoelectric material is made when subjected to in-plane mechanical, electric and magnetic impacts. The Laplace and Fourier transforms are applied to reduce the associated initial- and mixed-boundary value problem to dual integral equations, and then to singular integral equations with Cauchy kernel. By numerically solving the resulting equation, the dynamic field intensity factors as well as CODs, and energy release rates near the crack tip are evaluated and presented graphically. The effects of applied magnetic and electric impacts on crack growth are discussed. Obtained results show that, different from the static results, applied magnetic and electric impacts can strongly affect dynamic stress intensity factors.     

Edge cracked piezoelectric ceramic block under electromechanical impact loading

The dynamic field intensity factors and energy release rates in a piezoelectric ceramic block containing an edge crack with the condition of continuous electric crack faces under electromechanical impact loading are obtained. Integral transform method is used to reduce the problem to two pairs of dual integral equations, which are then expressed to an Fredholm integral equation of the second kind. Numerical values on the dynamic stress intensity factor and dynamic energy release rate are obtained to show the influence of the geometry and electric field.     

Transient response of an insulating crack near to the interface between two piezoelectric half-planes under electromechanical impacts by BEM

A time-domain boundary element method (BEM) together with the sub-domain technique is applied to study transient response of an insulating crack near to the interface between two anisotropic piezoelectric half-planes under electromechanical impacts. The present time-domain BEM uses a quadrature formula for the temporal discretization to approximate the convolution integrals and a collocation method for the spatial discretization. Quadratic quarter-point elements are implemented at the crack tip. A displacement extrapolation technique is used to determine the dynamic stress intensity factors (DSIFs) and the dynamic electrical displacement intensity factor. Numerical examples are presented to show the effects of load combination, geometric configuration and material combination on dynamic intensity factors and dynamic energy release rate.     

The shielding effect of the imperfect interface on a mode III permeable crack in a layered piezoelectric sensor

The mechanical model is established for a piezoelectric sensor with a mode III permeable crack parallel to the imperfect interface. Fracture analysis is performed by the standard methods of Fourier transform and singular integral equation. Three conclusions are drawn: (a) the imperfect interface has a shielding effect on the crack parallel and very near to it; (b) the shielding effect depends on the structural stiffness and the distance between the crack and interface; (c) for the electrically permeable crack, mechanical imperfection has more remarkable shielding effect than dielectric imperfection does.     

Transient thermoelectroelastic response of a functionally graded piezoelectric strip with a penny-shaped crack

Transient response of a penny-shaped crack in a plate of a functionally graded piezoelectric material (FGPM) is studied under thermal shock loading conditions. It is assumed that the thermoelectroelastic properties of the strip vary continuously along the thickness of the strip, and that the crack faces are completely insulated. By using both the Laplace and Hankel transforms, the thermal and electromechanical problems are reduced to a singular integral equation and a system of singular integral equations which are solved numerically. The intensity factors vs. time for various crack size, crack position and material nonhomogeneity are obtained.     

Transient response of a piezoelectric material with a semi-infinite mode-III crack under impact loads

The problem of a semi-infinite impermeable mode-III crack in a piezoelectric material is considered under the action of impact loads. For the case when a pair of concentrated anti-plane impact loads and electric displacements are exerted symmetrically on the upper and lower surfaces of the crack, the asymptotic electroelastic field ahead of the crack tip is determined in explicit form. The dynamic intensity factors of electroelastic field and dynamic mechanical strain energy release rate are obtained. The obtained results can be taken as fundamental solutions, from which general results may directly be evaluated by integration. The method adopted is to reduce the mixed initial-boundary value problem, by using the Laplace and Fourier transforms, into two simultaneous dual integral equations. One may be converted into an Abel's integral equation and the other into a singular integral equation with Cauchy kernel, and the solutions of both equations can be determined in closed-form, respectively. For some particular cases, the present results reduce to the previous results.