Transient responses of a multilayered spherically isotropic piezoelectric hollow sphere

The dynamic solution of a multilayered spherically isotropic piezoelectric hollow sphere subjected to radial dynamic loads is obtained. By the method of superposition, the solution is divided into two parts: one is quasi-static and the other is dynamic. The quasi-static part is derived by the state-space method, and the dynamic part is obtained by the method of separation of variables coupled with the initial parameter method as well as the orthogonal expansion technique. By using the quasi-static and dynamic parts, the electric boundary conditions as well as the electric continuity conditions, a Volterra integral equation of the second kind with respect to a function of time is derived, which can be solved successfully by means of the interpolation method. The displacements, stresses and electric potentials are finally obtained. The present method is suitable for a multilayered spherically isotropic piezoelectric hollow sphere consisting of arbitrary layers and subjected to arbitrary spherically symmetric dynamic loads. Finally, numerical results are presented and discussed.     

Dynamic solution of a multilayered orthotropic piezoelectric hollow cylinder for axisymmetric plane strain problems

The dynamic solution of a multilayered orthotropic piezoelectric infinite hollow cylinder in the state of axisymmetric plane strain is obtained. By the method of superposition, the solution is divided into two parts: one is quasi-static and the other is dynamic. The quasi-static part is derived by the state space method, and the dynamic part is obtained by the separation of variables method coupled with the initial parameter method as well as the orthogonal expansion technique. By using the obtained quasi-static and dynamic parts and the electric boundary conditions as well as the electric continuity conditions, a Volterra integral equation of the second kind with respect to a function of time is derived, which can be solved successfully by means of the interpolation method. The displacements, stresses and electric potentials are finally obtained. The present method is suitable for a multilayered orthotropic piezoelectric infinite hollow cylinder consisting of arbitrary layers and subjected to arbitrary axisymmetric dynamic loads. Numerical results are finally presented and discussed.     

层合球面各向同性热释电空心球的瞬态响应

运用叠加原理，将层合球面各向同性热释电空心球的球对称动力学问题的解分成准静态和动 态两部分，准静态部分首先运用状态空间法给出了显式表达式，然后运用分离变量法、初参 数法和特征函数展开技术，给出了动态部分的表示式，再结合内外表面上的电学边界条件和 界面上的电学连续条件，导出一个关于时间函数的第二类Volterra积分方程，运用插值法 可成功地给出此积分方程的高精度数值解，最终可求得原问题的位移、应力、电位移以及电 势的响应. 此方法适用任意层数且各层是任意厚度的层合热释电空心球作用随时间以任意形 式变化的球对称温度场. 文中还给出了数值结果.     

Isolated crack in three-dimensional piezoelectric solid. Part II: Stress intensity factors for circular crack

This is part II of the work concerned with finding the stress intensity factors for a circular crack in a solid with piezoelectric behavior. The method of solution involves reducing the problem to a system of hypersingular integral equations by application of the unit concentrated displacement discontinuity and the unit concentrated electric potential discontinuity derived in part I [1]. The near crack border elastic displacement, electric potential, stress and electric displacement are obtained. Stress and electric displacement intensity factors can be expressed in terms of the displacement and the potential discontinuity on the crack surface. Analogy is established between the boundary integral equations for arbitrary shaped cracks in a piezoelectric and elastic medium such that once the stress intensity factors in the piezoelectric medium can be determined directly from that of the elastic medium. Results for the penny-shaped crack are obtained as an example.     

Analytical solution for the electroelastic dynamics of a nonhomogeneous spherically isotropic piezoelectric hollow sphere

Summary By introduction of a special dependent variable and separation of variables technique, the electroelastic dynamic problem of a nonhomogeneous, spherically isotropic hollow sphere is transformed to a Volterra integral equation of the second kind about a function of time. The equation can be solved by means of the interpolation method, and the solutions for displacements, stresses, electric displacements and electric potential are obtained. The present method is suitable for a piezoelectric hollow sphere with an arbitrary thickness subjected to arbitrary mechanical and electrical loads. Numerical results are presented at the end.The work was supported by the National Natural Science Foundation of China (No. 10172075 and No. 10002016).     

New numerical method for volterra integral equation of the second kind in piezoelastic dynamic problems

The elastodynamic problems of piezoelectric hollow cylinders and spheres under radial deformation can be transformed into a second kind Volterra integral equation about a function with respect to time, which greatly simplifies the solving procedure for such elastodynamic problems. Meanwhile, it becomes very important to find a way to solve the second kind Volterra integral equation effectively and quickly. By using an interpolation function to approximate the unknown function, two new recursive formulae were derived, based on which numerical solution can be obtained step by step. The present method can provide accurate numerical results efficiently. It is also very stable for long time calculating.     

Dynamic response of elastic layer on stiff foundation under time harmonic surface vertical concentrated load

In this paper, the line-load integral equation method proposed in reference [1] is first used for solving the elastodynamic problems. A set of one-dimensional regular integral equation is derived for calculating the dynamic response of elastic layer on stiff foundation under time harmonic surface vertical concentrated load. And the numerical solution of the integral equation is obtained.     

SCATTERING OF THE HARMONIC STRESS WAVE BY CRACKS IN FUNCTIONALLY GRADED PIEZOELECTRIC MATERIALS

The present paper considers the scattering of the time harmonic stress wave by a single crack and two collinear cracks in functionally graded piezoelectric material （FGPM）. It is assumed that the properties of the FGPM vary continuously as an exponential function. By using the Fourier transform and defining the jumps of displacements and electric potential components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement and electric potential components across the crack surface are expanded in a series of Jacobi polynomials. Numerical examples are provided to show the influences of material properties on the dynamic stress and the electric displacement intensity factors.     

Boundary element analysis of interaction between an elastic rectangular inclusion and a crack

The interaction between an elastic rectangular inclusion and a kinked crack inan infinite elastic body was considered by using boundary element method. The new complexboundary integral equations were derived. By introducing a complex unknown function H(t)related to the interface displacement density and traction and applying integration by parts,the traction continuous condition was satisfied automatically. Only one complex boundaryintegral equation was obtained on interface and involves only singularity of order l/ r. Toverify the validity and effectiveness of the present boundary element method, some typicalexamples were calculated. The obtained results show that the crack stress intensity factorsdecrease as the shear modulus of inclusion increases. Thus, the crack propagation is easiernear a softer inclusion and the harder inclusion is helpful for crack arrest.     

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.     

Finite-part integral and boundary element method to solve three-dimensional crack problems in piezoelectric materials

Using the hypersingular integral equation method based on body force method, a planar crack in a three-dimensional transversely isotropic piezoelectric solid under mechanical and electrical loads is analyzed. This crack problem is reduced to solve a set of hypersingular integral equations. Compare with the crack problems in elastic isotropic materials, it is shown that for the impermeable crack, the intensity factors for piezoelectric materials can be obtained from those for elastic isotropic materials. Based on the exact analytical solution of the singular stresses and electrical displacements near the crack front, the numerical method of the hypersingular integral equation is proposed by the finite-part integral method and boundary element method, which the square root models of the displacement and electric potential discontinuities in elements near the crack front are applied. Finally, the numerical solutions of the stress and electric field intensity factors of some examples are given.     

Penny shaped crack in a piezoelectric cylinder surrounded by an elastic medium subjected to combined in-plane mechanical and electrical loads

The fracture problem of a penny shaped crack in a piezoelectric ceramic cylinder surrounded by an infinite elastic medium under in-plane normal mechanical and electrical loads is considered with the electric continuous boundary conditions on the crack surface. By using the potential theory and Hankel transform, a system of dual integral equations is obtained, and expressed to a Fredholm integral equation of the second kind. The mechanical and electrical field equations and all sorts of field intensity factors of mode I are obtained, and the numerical values of various field intensity factors for PZT-6B piezoelectric ceramic surrounded by several different elastic media are graphically shown for a uniform load and a ring-shaped load, respectively. And the effects of the size of the piezoelectric cylinder and the elastic material properties on various field intensity factors are obtained.     

DYNAMIC BEHAVIOR OF TWO PARALLEL SYMMETRY CRACKS IN MAGNETO-ELECTRO-ELASTIC COMPOSITES UNDER HARMONIC ANTI-PLANE WAVES

The dynamic behavior of two parallel symmetry cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied by Schmidt method. By using the Fourier transform, the problem can be solved with a pair of dual integral equations in which the unknown variable is the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surface were expanded in a series of Jacobi polynomials. The relations among the electric filed, the magnetic flux and the stress field were obtained. From the results, it can be obtained that the singular stresses in piezoelectric/piezomagnetic materials carry the same forms as those in a general elastic material for the dynamic anti-plane shear fracture problem. The shielding effect of two parallel cracks was also discussed.     

A crack in a piezoelectric layer bonded to a dissimilar elastic layer under transient load

Summary  A piezoelectric layer bonded to the surface of an elastic structure is considered. The piezoelectric and the elastic layers are infinite along the x-axis and have finite thickness in the y-direction. The polarization direction of the piezoelectric material is along the y-axis. By means of the method of singular integral equations, the solution in a Laplace transform plane is demonstrated. Laplace inversion yields the results in the time domain. Numerical values of the crack tip fields under in-plane transient electromechanical loading are obtained. The influence of layers thickness on stress and electric displacement intensity factors is investigated. Received 16 March 2000; accepted for publication 16 August 2000     

PMN-PT层/弹性基底结构中声表面波特性分析

研究了PMN-PT 压电层/弹性(金刚石) 基底结构中表面波的传播特性,压电层表面是机械自由的,电学边界条件分为电学开路和电学短路,压电层与基底之间采用理想连接. 得到了满足控制方程和边界条件的电弹场以及弹性波在结构中传播时的频散方程,通过数值算例分析了压电材料PMN-PT 的极化方向对弹性波频散曲线和机电耦合系数的影响,以及不同极化方向时弹性位移和电势随结构深度方向的变化,结果可为PMN-PT 压电材料在高频声表面波器件中的应用提供有价值的理论参考.      

Mechanics of indentation for piezoelectric thin films on elastic substrate

Frictionless normal indentation problem of rigid flat-ended cylindrical, conical and spherical indenters on piezoelectric film, which is either in frictionless contact with or perfectly bonded to an elastic half-space (substrate), is investigated. Both conducting and insulating indenters are considered. With Hankel transform, the general solutions of the homogeneous governing equations for the piezoelectric layer and the elastic half-space are presented. Using the boundary conditions for a vertical point force or a point electric charge, and the boundary conditions on the film/substrate interface, the Green’s functions can be obtained by solving sets of simultaneous linear algebraic equations. The solution of the indentation problem is obtained by integrating these Green’s functions over the contact area with unknown surface tractions or electric charge distribution, which will be determined from the boundary conditions on the contact surface between the indenter and the film. The solution is expressed in terms of dual integral equations that are converted to a Fredholm integral equation of the second kind and solved numerically. Numerical examples are also presented. The comparison between two film/substrate bonding conditions is made. It shows that the indentation rigidity of the film/substrate system is lower when the film is in frictionless contact with the substrate. The effects of the Young’s modulus and Poisson’s ratio of the elastic substrate, indenter electrical condition and indenter prescribed electric potential on the indentation responses are presented.     

Thermopiezoelectric interaction of macro- and micro-cracks in piezoelectric medium

Considered is the interaction of macro-and micro-cracks in an anisotropic piezoelectric solid. The Green’s function and principle of superposition are used to formulate a system of singular integral equations for solving the unknown temperature discontinuity and elastic displacement-electric potential. The residual heat flux, stress and electric displacement on the microcrack are evaluated directly from the near-tip field of main crack. Numerical results for stress and electric displacement intensity factors in a three-crack system are obtained to illustrate the application of the method.     

Dynamic behavior of piezoelectric/magnetostrictive composite hollow cylinder

The dynamic behavior of a multilayered, perfectly bonded piezoelectric/magnetostrictive composite hollow cylinder under radial deformation is investigated. The superposition method, the state space method as well as the separation of variables method are elegantly integrated in the solution approach. The governing equations are finally transformed into two Volterra integral equations of the second kind with respect to two functions of time. The elastic, electric and magnetic fields are finally obtained according to solving the integral equations. Free vibrations and transient responses are demonstrated by numerical experiments.     

SH波在压电材料条中垂直界面裂纹处的散射

研究了SH波在压电材料条中裂纹处的散射.压电材料条两侧涂有相同梯度参数的两个半无限大功能梯度材料,裂纹垂直于界面.通过Fourier变换,利用边界条件把问题转化为柯西核奇异积分方程,然后利用Chebyshev多项式对奇异积分方程进行数值求解.通过数值计算,分析讨论了压电条的几何参数和SH波频率对标准动应力强度因子的影响.     

Numerical analysis for a crack in piezoelectric material under impact

In this paper, a numerical analysis of impact interfacial fracture for a piezoelectric bimaterial is provided. Starting from the basic equilibrium equation, a dynamic electro-mechanical FEM formulation is briefly presented. Then, the path-independent separated dynamic J integral is extended to piezoelectric bimaterials. Based on the relationship of the path-independent dynamic J integral and the stress and electric displacement intensity factors, the component separation method is used to calculate the stress and electric displacement intensity factors for piezoelectric bimaterials in this finite-element analysis. The response curves of the dynamic J integral, the stress and electric displacement intensity factors are obtained for both homogeneous material (PZT-4 and CdSe) and CdSe/PZT-4 bimaterial. The influences of the piezoelectricity and the electro-mechanical coupling factor on these responses are discussed. The effects of an applied electric field are also discussed.