On the crack-tip stress singularity of interfacial cracks in transversely isotropic piezoelectric bimaterials

In this paper, characteristics of the interface crack-tip stress and electric displacement fields in transversely isotropic piezoelectric bimaterials are studied. The authors have proven, within the framework of the generalized Stroh formalism for piezoelectric bimaterials, that there is no coexistence of the parameters (oscillating) and κ (non-oscillating) in the interface crack-tip generalized stress field for all transversely isotropic piezoelectric bimaterials. This leads to the classification of piezoelectric bimaterials into one group that exhibits the oscillating property in the interface crack-tip generalized stress field and the other that does not. Fifteen (15) pair-combinations of six (6) piezoelectric materials PZT-4, PZT-5H, PZT-6B, PZT-7A, P-7, and BaTiO₃, which are commonly used in practice, are numerically analyzed in this study, and the results backup the above theoretical conclusions. Moreover, the associated eigenvectors for such material systems (with either =0 or κ=0) are also obtained numerically, and the result show that there still exist four linear independent associate eigenvectors for each bimaterial.     

ANALYSIS OF CRACK-TIP SINGULARITIES FOR AN INTERFACIAL PERMEABLE CRACK IN METAL/PIEZOELECTRIC BIMATERIALS

By modeling metal as a special piezoelectric material with extremely small piezoelec- tricity and extremely large permittivity,we have obtained the analytical solutions for an interfacial permeable crack in metal/piezoelectric bimaterials by means of the generalized Stroh formalism. The analysis shows that the stress fields near a permeable interfacial crack tip are usually with three types of singularities:r~(-1/2 iε)and r~(-1/2).Further numerical calculation on the oscillatory indexεare given for 28 types of metal/piezoelectric bimaterials combined by seven commercial piezoelectric materials: PZT-4,BaTiO_3,PZT-5H,PZT-6B,PZT-7A,P-7 and PZT-PIC 151 and four metals:copper,silver,lead and aluminum,respectively.The explicit expressions of the crack tip energy release rate(ERR)and the crack tip generalized stress intensity factors(GSIF)are obtained.It is found that both the ERR and GSIF are independent of the electric displacement loading,although they seriously depends on the mechanical loadings.     

Arbitrarily oriented crack near interface in piezoelectric bimaterials

Arbitrarily oriented crack near interface in piezoelectric bimaterials is considered. After deriving the fundamental solution for an edge dislocation near the interface, the present problem can be expressed as a system of singular integral equations by modeling the crack as continuously distributed edge dislocations. In the paper, the dislocations are described by a density function defined on the crack line. By solving the singular integral equations numerically, the dislocation density function is determined. Then, the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) at the crack tips are evaluated. Subsequently, the influences of the interface on crack tip SIFs, EDIF, and the mechanical strain energy release rate (MSERR) are investigated. The J-integral analysis in piezoelectric bimaterals is also performed. It is found that the path-independent of J₁-integral and the path-dependent of J₂-integral found in no-piezoelectric bimaterials are still valid in piezoelectric bimaterials.     

On the dynamic behaviour of a piezoelectric laminate with multiple interfacial collinear cracks

This paper investigates the dynamic behaviour of a piezoelectric laminate containing multiple interfacial collinear cracks subjected to steady-state electro-mechanical loads. Both the permeable and impermeable boundary conditions are examined and discussed. Based on the use of integral transform techniques, the problem is reduced to a set of singular integral equations, which can be solved using Chebyshev polynomial expansions. Numerical results are provided to show the effect of the geometry of interacting collinear cracks, the applied electric fields, the electric boundary conditions along the crack faces and the loading frequency on the resulting dynamic stress intensity and electric displacement intensity factors.     

Subinterface crack in an anisotropic piezoelectric bimaterial

In this paper, the problem of a subinterface crack in an anisotropic piezoelectric bimaterial is analyzed. A system of singular integral equations is formulated for general anisotropic piezoelectric bimaterial with kernel functions expressed in complex form. For commonly used transversely isotropic piezoelectric materials, the kernel functions are given in real forms. By considering special properties of one of the bimaterial, various real kernel functions for half-plane problems with mechanical traction-free or displacement-fixed boundary conditions combined with different electric boundary conditions are obtained. Investigations of half-plane piezoelectric solids show that, particularly for the mechanical traction-free problem, the evaluations of the mechanical stress intensity factors (electric displacement intensity factor) under mechanical loadings (electric displacement loading) for coupled mechanical and electric problems may be evaluated directly by considering the corresponding decoupled elastic (electric) problem irrespective of what electric boundary condition is applied on the boundary. However, for the piezoelectric bimaterial problem, purely elastic bimaterial analysis or purely electric bimaterial analysis is inadequate for the determination of the generalized stress intensity factors. Instead, both elastic and electric properties of the bimaterial’s constants should be simultaneously taken into account for better accuracy of the generalized stress intensity factors.     

Transient response of a cracked piezoelectric strip under arbitrary electro-mechanical impact

By using the well-developed integral transform methodology, the dynamic response of stress and electric displacement around a finite crack in an infinite piezoelectric strip are investigated under arbitrary dynamic anti-plane loads. The dynamic stress intensity factors and electric displacement are obtained analytically. It is shown that the dynamic crack-tip stress and electric field still have a square-root singularity. Numerical computations for the dynamic stress intensity factor show that the electric load has a significant influence on the dynamic response of stress field. The higher the ratio of the crack length to the width of the strip, the higher the peak value of the dynamic stress intensity factor is. On the other hand, the dynamic response of the electric field is determined solely by the applied electric load. The electric field will promote or retard the propagation of the crack depending on the time elapse since the application of the external electro-mechanical loads. The project supported by the National Natural Science Foundation of China and the Post-Doctor Science Foundation of China     

压电体中孔边Ⅲ型界面裂纹的动应力强度因子

采用Green函数法研究界面上含圆孔边界径向有限长度裂纹的两半无限压电材料对SH波的散射和裂纹尖端动应力强度因子问题.首先构造出具有半圆型凹陷半空间的位移Green函数和电场Green函数,然后采用裂纹"切割"方法构造孔边裂纹,并根据契合思想和界面上的连接条件建立起求解问题的定解积分方程.最后作为算例,给出了孔边界面裂纹尖端动应力强度因子的计算结果图并进行了讨论.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

An interface crack with partially electrically conductive crack faces under antiplane mechanical and in-plane electric loadings

An interface crack in a bimaterial piezoelectric space under the action of antiplane mechanical and in-plane electric loadings is analyzed. One zone of the crack faces is electrically conductive while the other part is electrically permeable. All electro-mechanical values are presented using sectionally-analytic vector-functions and a combined Dirichlet-Riemann boundary value problem is formulated. An exact analytical solution of this problem is obtained. Simple analytical expressions for the shear stress, electric field and also for mechanical displacement jump of the crack faces are derived. These values are also presented graphically along the corresponding parts of the material interface. Singular points of the shear stress, electric field and electric displacement jump are found. Their intensity factors are determined as well. Intensity factors variations with respect to the external electric field and different ratios between the electrically conductive and electrically permeable crack face zones are also demonstrated.     

Thermoelectroelastic Green's function and its application for bimaterial of piezoelectric materials

Summary For a two-dimensional piezoelectric plate, the thermoelectroelastic Green's functions for bimaterials subjected to a temperature discontinuity are presented by way of Stroh formalism. The study shows that the thermoelectroelastic Green's functions for bimaterials are composed of a particular solution and a corrective solution. All the solutions have their singularities, located at the point applied by the dislocation, as well as some image singularities, located at both the lower and the upper half-plane. Using the proposed thermoelectroelastic Green's functions, the problem of a crack of arbitrary orientation near a bimaterial interface between dissimilar thermopiezoelectric material is analysed, and a system of singular integral equations for the unknown temperature discontinuity, defined on the crack faces, is obtained. The stress and electric displacement (SED) intensity factors and strain energy density factor can be, then, evaluated by a numerical solution at the singular integral equations. As a consequence, the direction of crack growth can be estimated by way of strain energy density theory. Numerical results for the fracture angle are obtained to illustrate the application of the proposed formulation. Received 10 November 1997; accepted for publication 3 February 1998     

SINGULAR BEHAVIORS OF INTERFACIAL CRACKS IN 2D MAGNETOELECTROELASTIC BIMATERIALS

The singular characteristics of stress, electric displacement and magnetic induction fields near the tip of impermeable interracial cracks in two-dimensional magnetoelectroelastic bimaterials are studied using the generalized Stroh formalism. Two types of singularities are obtained： one is the oscillating singularity 1/2±iε, the other is the non-oscillating singularity 1/2±κ. It is found that the non-zero parameters ε and κ cannot coexist for one transversely isotropic MEE bimaterial, a similar result is obtained for transversely isotropic piezoelectric bimaterials.     

Impact response of an interface crack in a hybrid piezoelectric laminate

Summary  In a hybrid laminate containing an interfacial crack between piezoelectric and orthotropic layers, the dynamic field intensity factors and energy release rates are obtained for electro-mechanical impact loading. The analysis is performed within the framework of linear piezoelectricity. By using integral transform techniques, the problem is reduced to the solution of a Fredholm integral equation of the second kind, which is obtained from one pair of dual integral equations. Numerical results for the dynamic stress intensity factor show the influence of the geometry and electric field. Received 29 June 2001; accepted for publication 3 December 2001     

Dynamic response of a cracked piezoelectric ceramic under arbitrary electro-mechanical impact

The dynamic response of a cracked piezoelectric ceramic under in-plane electric and anti-plane mechanical impact is investigated by the integral transform method. The electric and mechanical loads are assumed to be arbitrary functions of time. It is shown that the dynamic crack-tip stress and electric displacement fields still have a square-root singularity. Numerical computations for the dynamic stress intensity factor show that the electric load has a significant influence on the dynamic response of stress field. On the other hand, the dynamic response of the electric field is determined solely by the applied electric field. The electric field will promote or retard the propagation of the crack depending on the time elapsed since the application of the external electro-mechanical loads.     

反平面扩展裂纹的J积分与COD

给出一个以任意速率扩展的反平面裂纹与路径无关的Ｊ积分，证明Ｊ积分扩展裂纹尖端的张开位移（动态ＣＯＤ）之间有的简单的关系，Ｊ积分与能量释放率，动应力强度因子之间也有简单关系，利用这些关系，给出了动态ＣＯＤ与动应力强度因子之间的关系式。     

平行于异材界面的三维椭圆平片裂纹分析

讨论了拉伸载荷作用下平行于两相材料界面的椭圆平片裂纹问题．首先,使用有限部积分概念和两相材料界面完全接合时的点力基本解导出了一组以裂纹表面位移差为未知函数的超奇异积分方程组．该组方程表明,此时三种裂纹模型同时存在;其次,在数值求解该组方程的过程中,未知函数裂纹表面位移差被近似为位移差的基本密度函数与多项式之积．基本密度函数反映了裂纹前沿应力奇性性态;最后,以拉伸载荷为例,讨论了椭圆平片裂纹与界面的距离、裂纹形状比和不同材料组合对应力强度因子的影响,并以图表形式给出。     

The scattering of P-waves by a piezoelectric particle with FGPM interfacial layers in a polymer matrix

Propagation of P-wave in an unbounded elastic polymer medium which contains a set of nested concentric spherical piezoelectric inhomogeneities is formulated. The polymer matrix is made of Epoxy and is isotropic; each phase of the inhomogeneity is made of a different piezoelectric material and is radially polarized and has spherical isotropy. Note that the individual phases are homogeneous, and all interfaces are perfectly bonded. The scattered displacement and electric potentials in the matrix are expressed in terms of spherical wave vector functions and Legendre functions, respectively. The transmitted displacement and electric potentials within each phase of the piezoelectric particle are expressed in terms of Legendre functions. The equations of motion and electrostatics in each phase of the piezoelectric inhomogeneity lead to a system of coupled second order differential equations, which is solved using the generalized Frobenius series. The present theory is extended to the case where the core of the inhomogeneity is made of PZT-4 and its coating is made of functionally graded piezoelectric material (FGPM) whose microstructural composition varies smoothly from PZT-4 at the core–coating interface to Epoxy at the coating–matrix interface. The effects of different types of variation in the electro-mechanical properties of FGPM on scattering cross-section and other electro-mechanical fields are addressed. The present theory is valid for arbitrary coating thickness, and arbitrary frequencies.     

Charged dislocations in piezoelectric bimaterials

In some piezoelectric semiconductors and ceramic materials, dislocations can be electrically active and could be even highly charged. However, the impact of dislocation charges on the strain and electric fields in piezoelectric and layered structures has not been presently understood. Thus, in this paper, we develop, for the first time, a charged three-dimensional dislocation loop model in an anisotropic piezoelectric bimaterial space to study the physical and mechanical characteristics which are essential to the design of novel layered structures. We first develop the analytical model based on which a line-integral solution can be derived for the coupled elastic and electric fields induced by an arbitrarily shaped and charged three-dimensional dislocation loop. As numerical examples, we apply our solutions to the typical piezoelectric AlGaN/GaN bimaterial to analyze the fields induced by charged square and elliptic dislocation loops. Our numerical results show that, except for the induced elastic (mechanical) displacement, charges along the dislocation loop could substantially perturb other induced fields. In other words, charges on the dislocation loop could significantly affect the traditional dislocation-induced stress/strain, electric displacement, and polarization fields in piezoelectric bimaterials.     

Effects of electric field on crack growth for a penny-shaped dielectric crack in a piezoelectric layer

The assumptions of impermeable and permeable cracks give rise to significant errors in determining electro-elastic behavior of a cracked piezoelectric material. The former simply imposes that the permittivity or electric displacement of the crack interior vanishes, and the latter neglects also the effects of the dielectric of an opening crack interior. Considering the presence of the dielectric of an opening crack interior and the permeability of the crack surfaces for electric field, this paper analyzes electro-elastic behavior induced by a penny-shaped dielectric crack in a piezoelectric ceramic layer. In the cases of prescribed displacement or prescribed stress at the layer surfaces, the Hankel transform technique is employed to reduce the problem to Fredholm integral equations with a parameter dependent nonlinearly on the unknown functions. For an infinite piezoelectric space, a closed-form solution can be derived explicitly, while for a piezoelectric layer, an iterative technique is suggested to solve the resulting nonlinear equations. Field intensity factors are obtained in terms of the solution of the equations. Numerical results of the crack opening displacement intensity factors are presented for a cracked PZT-5H layer and the effect of applied electric field on crack growth are examined for both cases. The results indicate that the fracture toughness of a piezoelectric ceramic is affected by the direction of applied electric fields, dependent on the elastic boundary conditions.     

Closed-form solution for two collinear cracks in a piezoelectric strip

Under the permeable electric boundary condition, the problem of two collinear anti-plane shear cracks lying at the mid-plane of a piezoelectric strip is investigated. By using the Fourier transform, the associated problem is reduced to a singular integral equation. Solving the resulting equation analytically, the electro-elastic field intensity factors and energy release rates at the inner and outer crack tips can be determined in explicit form. Numerical results for PZT-5H piezoelectric ceramic are also presented graphically. The results reveal that the effect of electric field on crack growth in piezoelectric materials is dependent on applied elastic displacement.     

Interface crack-tip generalized stress field and stress intensity factors in transversely isotropic piezoelectric bimaterials

Transversely isotropic piezoelectric (TIP) bimaterials with an impermeable interface crack have been classified [Int. J. Frac. 119 (2003) L41] into two classes corresponding to the vanishing of the two singularity parameters or κ. It is shown in the present paper that the related eigenvalue problems for either =0 or κ=0 are not degenerate. The crack-tip generalized stress fields are obtained subsequently. A new definition of crack-tip intensity factors is presented for interface cracks in practical TIP bimaterial of practical interest. Such defined intensity factors are real numbers, which dominate the maximum crack-tip stress singularity and do not generate any phase angle change under any dimension system transformation for physical quantities.