Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 53, 3, pp. 673-682, Warsaw 2015 DOI: 10.15632/jtam-pl.53.3.673 WAVE PROPAGATION IN PIEZOELECTRIC RINGS WITH RECTANGULAR CROSS-SECTIONS Xiaoming Zhang School of Mechanical and Power Engineering, Henan Polytechnic University, Jiaozuo, P.R. China e-mail: zxmworld11@163.com Youchao Wang School of Materials Science and Engineering, Henan Polytechnic University, Jiaozuo, P.R. China e-mail: wangyc@hpu.edu.cn Huitao Chen School of Mechanical and Power Engineering, Henan Polytechnic University, Jiaozuo, P.R. China e-mail: huitaochen@hpu.edu.cn The ring ultrasonic transducers are widely used in the ocean engineering andmedical fields. This paper employs an extended orthogonal polynomial approach to solve the guided wave propagation in two-dimensional structures, i.e. piezoelectric rings with rectangular cross- -sections. The extended polynomial approach can overcome the drawbacks of the conven- tional orthogonal polynomial approach which can be used to solve wave propagation in one-dimensional structures. Through numerical comparison with the available results for a rectangular aluminumbar, the validity of the present approach is illustrated.The dispersion curves and displacement and electric potential distributions of various rectangular piezo- electric rings are calculated, and the effects of different radius to thickness ratios, width to height ratios and polarizing directions on the dispersion curves are illustrated. Keywords: piezoelectric rings, orthogonal polynomial, guided wave, dispersion curves 1. Introduction With the development ofmaterials and advances inmanufacturing technology, piezoelectricma- terials having electromechanical coupling effects have foundextensive applications inmanysmart devices, and the behavior of the selected wave mode can directly affect the performance of the devices. Thus, it is very important to study the wave characteristics in piezoelectric structures. Much effort has beenmade both theoretically and experimentally to studywave propagation in piezoelectric structures by scientists and engineers,which is essential for the application of piezo- electricmaterials.Wave propagation along piezoelectric cylindrical rods of hexagonal 6, 622, and 6mm crystal symmetry was discussed byWilson andMorrison (1977). The propagating nature of the elastic and electric wave in bone andporousPZTwas investigated byChakraborty (2009). Using the extended Durbin method, Ing et al. (2013) investigated the transient elastic waves propagating in a two-layered piezoelectric medium The wave propagation behavior in layered piezoelectric structures and functionally graded piezoelectric material structures has also been studied bymany researches with different methods, such as the transfer matrix method (Cai et al., 2001), the layer elementmethod (Han et al., 2004), the orthogonal polynomial seriesmethod (Yu andMa, 2008; Yu et al., 2013; Singh and Rokne, 2013), and so on. As a common structure, hollow cylinder has been paid considerable attention on the wave propagation. For piezoelectric media, Paul and Venkatesan (1987) and Shul’ga (2002) studied three dimensional electroelastic waves and the axisymmetric waves in a hollow piezoelectric 674 X. Zhang et al. ceramic cylinder. Puzyrev andStorozhev (2011) studied the problemof electroelastic waves pro- pagating in piezoelectric hollow cylinders of sector cross section and analyzed mode asymptotic behavior and amplitude distributions of wave characteristics. Zenkour (2012) developed an ana- lytical solution to the axisymmetric problemof a radially polarized piezoelectric hollow cylinders subjected to electric, thermal and mechanical load. Using the Legendre orthogonal polynomial approach developed by Lefebvre et al. (1999) to solve thewaves inmultilayered plates, Yu et al. (2009) investigated thewave characteristics in functionally graded piezoelectric hollow cylinders. These investigations focused on one-dimensional structures, i.e. the hollow cylinders are infinite in axial direction andhave a finite dimension in only one direction. But in practical applications, many piezoelectric elements have very finite dimensions in two directions. One-dimensionalmo- dels are not suitable for these structures. Thus, it is significant to study thewave characteristics in 2-D piezoelectric structures, such as a piezoelectric ring with rectangular cross-section. The ring ultrasonic transducer has been widely used in ocean engineering and medical fields. But few investigations on the wave propagation in ring transducers have been reported. In this paper, we present an extended orthogonal polynomial series approach to solve the wave propagation in a 2-D structure, i.e. a piezoelectric ring with rectangular cross-section. The present approach can overcome the drawbacks of the conventional orthogonal polynomial approachwhich can onlydealwith the one-dimensional structures that have afinite dimension in only one direction, such as the axially infinite hollow cylinder, horizontally infinite flat plate.The dispersion curves and displacement and electric potential distributions of various piezoelectric rings with rectangular cross sections are shown. The effects of different width to height ratios, radius to thickness ratios andpolarizingdirections onthewave characteristics arealso illustrated. The investigating results can be used to direct the design and optimization of the piezoelectric ring transducers. 2. Mathematics and formulation of the problem In this seciton, we derive the analytical formulation of the problem in cylindrical coordinate (r,θ,z) with the z-axis coincidingwith the axis of the ring.Considering an orthotropic ringwith rectangular cross-section, h is height in z direction and a, b denote the inner and outer radius respectively, as shown in Fig. 1. The radius to thickness ratio is defined as η = b/(b−a) and the width to height ratio is d/h. The polarizing direction is in r or z direction. In this paper, traction free and open circuit boundary conditions are assumed. Fig. 1. Schematic of a piezoelectric ring with rectangular cross-section In cylindrical coordinate, the dynamic equation for the piezoelectric ring is governed by ∂Trr ∂r + 1 r ∂Trθ ∂θ + ∂Trz ∂z + Trr−Tθθ r = ρ ∂2ur ∂t2 ∂Trθ ∂r + 1 r ∂Tθθ ∂θ + ∂Tθz ∂z + 2Trθ r = ρ ∂2uθ ∂t2 (2.1) Wave propagation in piezoelectric rings... 675 ∂Trz ∂r + 1 r ∂Tθz ∂θ + ∂Tzz ∂z + Trz r = ρ ∂2uz ∂t2 ∂Dr ∂r + 1 r ∂Dθ ∂θ + ∂Dz ∂z + Dr r =0 where ui, Tij andDi denotemechanical displacement, the stress and electric displacement com- ponents respectively; ρ is the density of the material. The strain-displacement relations are εrr = ∂ur ∂r εθθ = 1 r ∂uθ ∂θ + ur r εzz = ∂uz ∂z εθz = 1 2 (∂uθ ∂z + ∂uz r∂θ ) εrz = 1 2 (∂ur ∂z + ∂uz ∂r ) εrθ = 1 2 (1 r ∂ur ∂θ + ∂uθ ∂r − uθ r ) (2.2) where εij denotes the strain. The constitutive equations for the piezoelectric ring with radial polarizing direction can be written in the following form Tθθ =C11εθθ +C12εzz +C13εrr +e31 ∂Φ ∂r Tzz = ( C12εθθ+C22εzz +C23εrr +e32 ∂Φ ∂r ) I(r,z) Trr = ( C13εθθ+C23εzz +C33εrr+e33 ∂Φ ∂r ) I(r,z) Trz = ( 2C44εrz +e24 ∂Φ ∂z ) I(r,z) Trθ = ( 2C55εrθ +e15 1 r ∂Φ ∂θ ) I(r,z) Tθz =2C66εθzI(r,z) (2.3) and Dθ =2e15εrθ− ǫ11 1 r ∂Φ ∂θ Dz = ( 2e24εrz − ǫ22 ∂Φ ∂z ) I(r,z) Dr = ( e31εθθ+e32εzz +e33εrr− ǫ33 ∂Φ ∂r ) I(r,z) (2.4) whereΦ denotes the electric potential.Cij, eij and ǫij are the elastic, piezoelectric and dielectric coefficients given in the crystallographic axes, respectively. I(r,z) is the rectangular window function, introduced so as tomeet the stress-free boundary conditions (Trr =Trθ =Trz =Tθz = Tzz =Dr =Dz =0 at the four boundaries), defined as I(y,z) = { 1 0¬ y¬ d and 0¬ z¬h 0 elsewhere (2.5) For a free harmonic plane wave propagating in the circumferential direction in a ring, we assume the displacement components, to be of the form ur(r,θ,z,t) = exp(ikbθ− iωt)U(r,z) uθ(r,θ,z,t) = exp(ikbθ− iωt)V (r,z) uz(r,θ,z,t) = exp(ikbθ− iωt)W(r,z) ϕ(r,θ,z,t) = exp(ikbθ− iωt)X(r,z) (2.6) whereU(r,z), V (r,z) andW(r,z) denote themechanical displacement amplitudes in the radial, circumferential and axial directions respectively, andX(r,z) represents the amplitude of electric potential. ω is the angular frequency, and k is the magnitude of the wave vector. 676 X. Zhang et al. Substituting Eqs. (2.2)-(2.6) into Eq. (2.1), the governing differential equations in terms of mechanical displacement and electric potential components, gives [C33(r 2U,rr+rU,r )−C11U −k 2b2C55U +C44r 2U,zz−ikb(C11+C55)V +(C23−C12)rW,z +ikb(C13+C55)rV,r+(C23+C44)r 2W,rz+e33(r 2X,rr+rX,r )−e31rX,r−k 2b2e15X +e24r 2X,zz ]I(r,z)+ [C33r 2U,r+C13r(ikbV +U)+C23r 2W,z+e33r 2X,r ]I(r,z),r +[C44r 2(U,z+W,r )+e24r 2X,z ]I(r,z),z=−ρr 2ω2U [C55(r 2V,rr+rV,r )− (C55+(kb) 2C11)V +ikb(C13+C55)rU,r+C66r 2V,zz +ikb(C11+C66)U +ikb(C12+C66)rW,z+(e31+e15)rX,r+2e15X]I(r,z) +C66(r 2V,z+ikbrW)I(r,z),z+[C55(r 2V,r+rV +ikbrU)+e15rX]I(r,z),r=−ρr 2ω2V [C44(r 2W,rr+rW,r )+C22r 2W,zz−(kb) 2C66W +(C12+C44)rU,z+(C23+C44)r 2U,rz +ikb(C12+C66)rV,z+e24rX,z+(e24+e32)r 2X,rz ]I(r,z)+ [C12r(ikbV +U) +C23r 2U,r+C22r 2W,z+e32r 2X,r ]I(r,z),z+[C44r 2(W,r+U,z )+e24r 2X,z ]I(r,z),r =−ρr2ω2V [e33(r 2U,rr+rU,r )−k 2b2e15U+e31rU,r+e24r 2U,zz+(e31+e15)rV,r−e15V +(e24+e32)r 2W,rz+e24rW,z−ǫ33(r 2X,rr+rX,r )− ǫ22r 2X,zz+(kb) 2ǫ11X]I(r,z) + [e24r 2U,z+e24r 2W,r−ǫ22r 2X,z ]I(r,z),z +[e31rU+e33r 2U,r+e31rV +e32r 2W,z−ǫ33r 2X,r ]I(r,z),r=0 (2.7) where subscript comma indicates partial derivative. To solve the coupled wave equation, we expand U(r,z), V (r,z), W(r,z) and X(r,z) into products of two Legendre orthogonal polynomial series U(r,z) = ∞ ∑ m,j=0 p1m,jQm(r)Qj(z) V (r,z) = ∞ ∑ m,j=0 p2m,jQm(r)Qj(z) W(r,z) = ∞ ∑ m,j=0 p3m,jQm(r)Qj(z) X(r,z) = ∞ ∑ m,j=0 p4m,jQm(r)Qj(z) (2.8) where pim,j (i=1,2,3,4) is the expansion coefficients and Qm(r)= √ 2m+1 b−a Pm (2r− b−a b−a ) Qn(z) = √ 2n+1 h Pn (2z−h h ) (2.9) with Pm and Pn representing the mth and the nth Legendre polynomial. The summation over the polynomials can be halted at some finite valuem=M andn=N, when higher order terms become essentially negligible. Multiplying each equation by Qj(r)Ql(z)e −jωt with j and l running respectively from zero toM and zero toN, and integrating over z fromzero tohand r froma to b and taking advantage of the orthonormality of the polynomials Qm(r) and Qn(z), Eqs. (2.7) can be reorganized into a form of the system problem A jlmn 11 p 1 m,n+A jlmn 12 p 2 m,n+A jlmn 13 p 3 m,n+A jlmn 14 p 4 m,n =−ω 2Mjlmnp 1 m,n A jlmn 21 p 1 m,n+A jlmn 22 p 2 m,n+A jlmn 23 p 3 m,n+A jlmn 24 p 4 m,n =−ω 2Mjlmnp 2 m,n A jlmn 31 p 1 m,n+A jlmn 32 p 2 m,n+A jlmn 33 p 3 m,n+A jlmn 34 p 4 m,n =−ω 2Mjlmnp 3 m,n A jlmn 41 p 1 m,n+A jlmn 42 p 2 m,n+A jlmn 43 p 3 m,n+A jlmn 44 p 4 m,n =0 (2.10) Wave propagation in piezoelectric rings... 677 whereA jlmn αβ (α,β =1,2,3,4) andMjlmn are the elements of a non-symmetric matrix. Equations (2.10)4 can be written as p4m,n =− ( A jlmn 44 ) −1( A jlmn 41 p 1 m,n+A jlmn 42 p 2 m,n+A jlmn 43 p 3 m,n ) (2.11) Substituting Equation (2.11) into equations (2.10)1, (2.10)2 and (2.10)3, gives [ A jlmn 11 −A jlmn 14 ( A jlmn 44 ) −1 A jlmn 41 ] p1m,n+ [ A jlmn 12 −A jlmn 14 ( A jlmn 44 ) −1 A jlmn 42 ] p2m,n + [ A jlmn 13 −A jlmn 14 ( A jlmn 44 ) −1 A jlmn 43 ] p3m,n =−ω 2Mjlmnp 1 m,n [ A jlmn 21 −A jlmn 24 ( A n,m 44 ) −1 A jlmn 41 ] p1m,n+ [ A jlmn 22 −A jlmn 24 ( A jlmn 44 ) −1 A jlmn 42 ] p2m,n + [ A jlmn 23 −A jlmn 24 ( A jlmn 44 ) −1 A jlmn 43 ] p3m,n =−ω 2Mjlmnp 1 m,n [ A jlmn 31 −A jlmn 34 ( A jlmn 44 ) −1 A jlmn 41 ] p1m,n+ [ A jlmn 32 −A jlmn 34 ( A jlmn 44 ) −1 A jlmn 42 ] p2m,n + [ A jlmn 33 −A jlmn 34 ( A jlmn 44 ) −1 A jlmn 43 ] p3m,n =−ω 2Mjlmnp 1 m,n (2.12) Then, Eqs. (2.12) can be recognized into     A jlmn 11 A jlmn 12 A jlmn 13 A jlmn 21 A jlmn 22 A jlmn 23 A jlmn 31 A jlmn 32 A jlmn 33          p1m,n p2m,n p3m,n      =−ω2    Mjlmn 0 0 0 Mjlmn 0 0 0 Mjlmn         p1m,n p2m,n p3m,n      (2.13) So,Eq. (2.13) forms the eigenvalue problemtobe solved.Theeigenvectors pim,n (i=1,2,3) allow the components of the displacement and p4m,n determines the electric potential distribution.The eigenvalue ω2 gives angular frequency. 3. Numerical results The computer programs in terms of the extended orthogonal polynomial approach have been written using Mathematica to calculate the dispersion curves and displacement and electric potential distributions for the piezoelectric rings. The physical properties of the piezoelectric material, PZT-4, are listed in Table 1. Here, the elastic constants of the radial polarizing ring and axial polarizing rings are the same to have a clear comparison. Table 1.Material parameters of the piezoelectric materials Property C11 C12 C13 C22 C23 C33 C44 C55 C66 13.9 7.8 7.4 13.9 7.4 11.5 2.56 2.56 3.05 r-polarization e15 e24 e31 e32 e33 ǫ11 ǫ22 ǫ33 ρ 12.7 12.7 −5.2 −5.2 15.1 650 650 560 7.5 a-polarization e34 e16 e23 e21 e22 ǫ11 ǫ22 ǫ33 ρ 12.7 12.7 −5.2 −5.2 15.1 650 560 650 7.5 Units:Cij [10 10N/m2], ǫij [10 −11F/m], eij [C/m 2], ρ [103kg/m3] 678 X. Zhang et al. 3.1. Approach validation To the authors’ knowledge, there are not published results on the wave propagation for a piezoelectric ring with rectangular cross-section so far. In order to check the effectiveness of the present approach and validate the computer program, we calculate a 16mmby 5mm rectangu- lar aluminum bar and make a comparison with previous results. The material parameters are ρ = 2.7 · 103kg/m3, C11 = C22 = C33 = 10.78 · 10 10Pa, C12 = C13 = C23 = 5.494 · 10 10Pa, C44 = C55 = C66 = 2.645 · 10 10Pa. Figure 2 is the corresponding dispersion curves, of which lines are fromLoveday (2006), and dotted lines are obtained from the present approach. As can be seen, the agreement between the present approach and the previous results is quite good. Fig. 2. Dispersion of propagating waves in a rectangular waveguide; lines: PhilipW. Loveday’s results, dotted lines: the authors’ results 3.2. Guided waves in piezoelectric rings with rectangular cross-sections Figure 3 shows the dispersion curves of the first four ordermodes for the PZT-4 ring with a square cross section and for the corresponding non-piezoelectric one with h=1mm, a=9mm, b = 10mm and η = 10. It can be seen that piezoelectricity has a significant effect on the dispersion curves. For any one specific mode, the phase velocities of non-piezoelectric ring are smaller than those of the corresponding piezoelectric one, and the piezoelectric effect is very little on the low order modes at low frequency and becomes stronger as the wave number and mode order increase. The wave number is usually very big and the operating frequency is very high in micro-scale SAWdevices. So, the piezoelectric effects will be prominent. Fig. 3. Dispersion curves of the first four order modes for a square ring with η=10: (a) phase velocity spectra, (b) frequency spectra, solid line, piezoelectric; dotted line, non-piezoelectric Figures 4 and 5 show the dispersion curves of the PZT-4 rings with different radius to thickness ratios (η = 10, η = 2) and different width to height ratios (d/h = 1/2, d/h = 1/4, d/h = 1/10), respectively. Figure 6 illustrates the dispersion curves of the PZT-4 ring with Wave propagation in piezoelectric rings... 679 Fig. 4. Phase velocity dispersion curves for piezoelectric square rings: (a) η=10, (b) η=2 Fig. 5. Phase velocity dispersion curves for piezoelectric rectangular rings with different width to thickness ratios: (a) d/h=1/2, (b) d/h=1/4, (c) d/h=1/10 Fig. 6. Phase velocity dispersion curves for a piezoelectric rectangular ring with axial polarization 680 X. Zhang et al. Fig. 7. Mechanical displacement and electric potential profiles of the first mode for a piezoelectric square ring with η=2 at kd=4.1 Fig. 8. Mechanical displacement and electric potential profiles of the first mode for a piezoelectric square ring with η=2 at kd=4.1 Wave propagation in piezoelectric rings... 681 axial polarization. From the curves, we find that different radius to thickness ratios and width to height ratios and polarizing direction all have significant influence on the dispersion curves. The first two modes have no cut-off frequencies, which is different from that for an infinite hollow cylinder in which only the first mode has no cut-off frequencies. In an infinite hollow cylinder, only the thickness direction is a finite dimension, but there are two finite dimensions in a rectangular ring.With the width to height ratio increasing, the difference between the first mode dispersion curve and the second one becomes small, and the cut-off frequencies become small. Thedisplacement and electric potential profiles of thefirstmode for a square ringwithη=10 at kd=4.1 and kd=40.1 are respectively shown in Figs. 7 and 8.We notice that displacement u and v and electric potential distributions are symmetry and displacementw is antisymmetry in axial direction, This is because the geometry and material propreties are symmetric in axial direction. The displacement and electric potential profiles distribute mainly near the outside edge at small wavenumber case and distribute around the four boundaries at big wavenumber case. 4. Conclusions and prospects The formulation to analyze the guidedwave in piezoelectric ringswith rectangular cross-sections using the extended orthogonal polynomial approach has beenpresented in this paper.According to the numerical results, we can draw the following conclusions: (a) The effects of the piezoelectricity on dispersion curves become stronger with the wave number andmode order increasing. (b) The width to height ratio, radius to thickness ratio and polarization all can significantly influence the guided wave characteristics in piezoelectric rings. (c) The displacement and electric potential distributions are symmetry in axial direction and distribute mainly near the outside edge. So, through changing the width to height ratio and the radius to thickness ratio of the piezoelectric ring, we can obtain the ring transducers with the dispersion features and field distributions that we want. We consider that the present approach could be of interest in non-destructive testing evalu- ation, and can deal with 2D structureswithmore complex cross sections andmulti-field coupled 2D structures. Acknowledgments The work is supported by the National Natural Science Foundation of China (No. 11272115) and Doctoral Fund of HenanPolytechnicUniversity (No. B2009-81) andFoundation forDistinguishedYoung Scholars of Henan Polytechnic University (No. J2013-08) and by the high-performance grid computing platform of Henan Polytechnic University. References 1. CaiC., LiuG.R., LamK.Y., 2001,A technique formodellingmultiple piezoelectric layers,Smart Materials And Structures, 10, 689 2. Chakraborty A., 2009, Wave Propagation in Porous Piezoelectric Media, Computer Modeling in Engineering and Sciences, 40, 105 3. Han X., Liu G.R., Ohyoshi T., 2004, Dispersion and characteristic surfaces of waves in hybrid multilayered piezoelectric circular cylinders,Computational Mechanics, 33, 334-344 682 X. Zhang et al. 4. Ing Y.S., Liao H.F., Huang K.S., 2013, The extended Durbin method and its application for piezoelectric wave propagation problems, International Journal of Solids and Structures, 50, 4000-4009 5. Lefebvre J.E., Zhang V., Gazalet J., Gryba T., 1999, Legendre polynomial approach for modeling free-ultrasonic waves in multilayered plates, Journal of Applied Physics, 85, 3419-3427 6. Loveday P.W., 2006, Numerical comparison of patch and sandwich piezoelectric transducers for transmitting ultrasonic waves, Smart Structures and Materials. International Society for Optics and Photonics, 6166, 1-8 7. Paul H.S., Venkatesan M., 1987, Vibrations of a hollow circular cylinder of piezoelectric cera- mics,The Journal of the Acoustical Society of America, 82, 952-956 8. Puzyrev V., Storozhev V., 2011, Wave propagation in axially polarized piezoelectric hollow cylinders of sector cross section, Journal of Sound and Vibration, 330, 4508-4518 9. Shul’gaH.A., 2002, Propagation of harmonic waves in anisotropic piezoelectric cylinders. Homo- geneous piezoceramic waveguides, International Applied Mechanics, 38, 933-953 10. Singh B.M., Rokne J., 2013, Propagation of SH waves in layered functionally gradient piezoelectric-piezomagnetic structures,Philosophical Magazine, 93, 1690-1700 11. WilsonL.O.,Morrison J.A., 1977,Wave propagation in piezoelectric rods of hexagonal crystal symmetry, Journal of Mechanics and Applied Mathematics, 30, 387-395 12. Yu J.G., Lefebvre J.E., GuoY.Q., 2013,Wave propagation inmultilayered piezoelectric sphe- rical plates,Acta Mechanica, 224, 1335-1349 13. Yu J.G., Ma Q.J., 2008, Circumferential wave in functionally graded piezoelectric cylindrical curved plates,Acta Mechanica, 198, 171-190 14. Yu J.G., Wu B., Chen G.Q., 2009, Wave characteristics in functionally graded piezoelectric hollow cylinders,Archive of Applied Mechanics, 79, 807-824 15. Zenkour A.M., 2012, Piezoelectric behavior of an inhomogeneous hollow cylinder with thermal gradient, International Journal of Thermophysics, 33, 1288-1301 Manuscript received May 6, 2014; accepted for print February 17, 2015