This paper presents the analysis and numerical results of rotation, propagation angle, and attenuation angle upon the waves propagating in the piezoelectric body. Via considering the centripetal and Coriolis accelerations in the piezoelectric equations with respect to a rotating frame of reference, wave velocities and attenuations are derived and plotted graphically. It is demonstrated that rotation speed vector can affect wave velocities and make the piezoelectric body behaves as if it was damping. Besides, the effects of propagation angle and attenuation angle are presented. Critical point is found when rotation speed is equal to wave frequency, around which wave characteristics change drastically.
1. Introduction
The gyroscopes of rotating motion sensors have important applications in automobiles, video cameras, smart weapon systems, machine control, robotics, and navigation. Recently, using piezoelectric materials to make vibratory gyroscopes has been of increasing interest. Piezoelectric gyroscopes can make use of rotation-induced frequency shifts in surface acoustic wave (SAW) or bulk acoustic wave (BAW) piezoelectric resonators to measure angular rates. The basic behaviors of a piezoelectric gyroscope are governed by the equations of a rotating piezoelectric body, which consist of the equations of linear piezoelectricity with rotation-related Coriolis and centrifugal accelerations.
As far as we know, the research of rotation-affected vibration or wave was started by Huston [1] who investigated the effect of “rigid-body” rotation on wave propagation velocities in elastic media. Further, rotation effect was studied in the “in-plane” vibration of rotating circular disks [2]. It was found that the inclusion of Coriolis and centripetal accelerations leads to the result that the medium behaves dispersive and anisotropic [3]. In another contribution [4], the reflection and refraction of plane waves were considered. The effect of rotation on surface acoustic waves was derived theoretically in a perturbation treatment of the Coriolis force for an isotropic medium [5]. The problem of wave propagation in a rotating random infinite magnetothermoviscoelastic medium was studied and a coupled dispersion relation for longitudinal and transverse waves was deduced to determine the effect of viscoelasticity, relaxation times, and rotation on the phase velocity of the coupled waves [6]. In the work of Wauer [7], the propagation of waves in a conducting piezoelectric solid was studied for the case when the entire medium rotates with a uniform angular velocity. Destrade and Saccomandi [8] raised and addressed two questions related to elastic motions and found some finite amplitude transverse waves in rotating incompressible elastic solids with general shear response. Auriault [9] revealed that free wave propagation in non-Galilean rotating media gives rise to two dispersive waves which are coupled dilatational-shear waves. The propagation of surface (Rayleigh) waves over a rotating orthorhombic crystal was studied [10], in which the secular equation for the surface wave speed was found explicitly. In the work of Ting [11], the Stroh formalism for surface waves in an anisotropic elastic half-space was extended to the case when the half-space rotates about an axis with a constant rotation rate. Auriault [12] investigated wave propagation in elastic porous media which are saturated by incompressible viscous Newtonian fluids when the porous media are in rotation with respect to a Galilean frame. Yang [13] presented a review of analyses on the vibrations of rotating piezoelectric structures for applications in piezoelectric angular rate sensors. Propagation of plane waves in a micropolar porous elastic solid rotating with a uniform angular velocity was investigated [14]. The paper [15] dealt with the propagation of body waves in a rotating, generalized thermoelastic solid by using Cardano’s and perturbation methods. A two-dimensional problem in electromagnetic micropolar generalized thermoelastic medium subjected to mechanical force or thermal source was investigated [16]. Biryukov et al. [17] investigated the gyroscopic effect in arbitrary crystals by taking into account the medium rotation. Recently, the paper [18] considered the propagation of body waves in a homogenous isotropic, rotating, generalized thermoelastic solid with voids. Wegert et al. [19] analysed theoretical upper bounds for the size of the gyroscopic effect on the frequency of guided acoustic waves in (piezo)elastic media, which are valid in the regime of small rotation rates as compared to the frequency of the guided acoustic wave. The contribution [20] was aimed at the effects of rotation on the propagation of harmonic plane waves under two-temperature thermoelasticity theory. Kothari and Mukhopadhyay [21] analyzed the effects of rotation on the propagation of harmonic plane waves in an unbounded thermoelastic media rotating with a uniform angular velocity. The investigation [22] was performed with the effect of rotation on an infinite circular cylinder subjected to certain boundary conditions.
As stated above, it is seen that many achievements have been done about the rotation effects on waves. This paper prefers to investigate the rotation on the inhomogeneous waves in piezoelectric body which have been researched by Yuan and his colleagues [23–26]. To our knowledge, no systematic empirical research exists addressing the question of inhomogeneous wave propagation in a rotating piezoelectric body. Our work here is to present the analysis and result for this problem in the framework of inhomogeneous wave theory.
The paper is organized in the following manner. In the next section, the basic equations for motion in a rotating piezoelectric solid and their wave dispersion equations to harmonic waves are given. Next, using the inhomogeneous wave theory, we recast the dispersion equations in a general complex form which separable real solutions to define the phase velocity and attenuation are admitted. Thus, we can discuss the wave phase velocities, attenuations with three independent parameters: propagation angle, attenuation angle, and rotation speed. Finally, in Sections 3 and 4, numerical results are presented and conclusions are inferred, respectively.
2. Basic Governing Equations
We consider a linear homogeneous piezoelectric body shown in Figure 1, and M is the material point rotating with the speed vector Ω(=Ω1e1+Ω2e2+Ω3e3). It should be mentioned that throughout this paper, all equations are expressed in the inertial frame ox1x2x39, in which there are base vector e1, e2, and e3 along three axes, respectively.
The rotating piezoelectric body.
Thus, the momentum balance in a piezoelectric body can be written as
(1)ρ[∂2u∂t2+Ω×(Ω×u)+2Ω×∂u∂t]=∇·σ,
and equivalently, in component form:
(2)ρ[∂2uj∂t2+εjikεkmnΩiΩmun+2εjikΩi×∂uk∂t]=σij,i.
In the above equation, ρ is the mass density, t is the time variable, u is the displacement vector, σ is the Cauchy stress tensor, and εjik is the permutation tensor. The subscripts range from 1 to 3. On account of rotation, the term Ω×(Ω×u) denotes the centripetal acceleration, and due to the time-varying motion, 2Ω×(∂u/∂t) corresponds to the Coriolis acceleration [3].
Further, the electric field can be described by the electrostatic equation
(3)Di,i=0,
where Di is the electric displacement vector, and with material equations
(4)σij=Cijklεkl-ekijEk,Di=ϵijEj+eiklεkl,
where σij are the strain tensor and Ek the electric field vector while Cijkl, ekij, and ϵij are the elasticity, piezoelectricity, and permittivity tensors of the material. The Einstein summation is implied in the above equations over the repeated subscripts.
The electric field vector can be derived from an electric potential, that is,
(5)Ek=-φ,k,
where φ is the electric potential. The geometric relationship between the strain and the displacement tensors is defined as
(6)εkl=12(uk,l+ul,k).
Eliminating εkl and Ek from (4), (5), and (6) yields
(7)σij=Cijkluk,l+ekijφ,k,Di=-ϵijφ,j+eikluk,l.
Then, instituting (7) into (2) and (3), one obtains a set of wave equations for rotating piezoelectric body
(8)ρ[∂2uj∂t2+εjikεkmnΩiΩmun+2εjikΩi∂uk∂t]-Cijkluk,li-ekijφ,ki=0-ϵijφ,ji+eikluk,li=0
which contains four equations (three elastic wave equations and one electric wave equation that is associated with the elastic waves by the material relationship equation (4) and three independent unknown variables: u1, u2, u3, and φ).
3. The Inhomogeneous Wave Solutions
Generally, the wave equations of (8) can be solved by introducing complex monochromatic plane wave functions, such as,
(9)ui=UieI(kjxj-ωt),φ=ΨeI(kjxj-ωt),
where kj is the complex wave vector, ω is the wave circular frequency, I is the imaginary unit (I=-1), t is the time variable, and (Ui,Ψ) are the complex amplitudes of displacements and electric potential, respectively. Inserting (9) into (8) gives
(10)ρω2[-Uj+εjikεkmnΩiωΩmωUn-2IεjikΩiωUk]+CijklUkklki+ekijkkkiΨ=0,ϵijkjkiΨ-eiklUkklki=0.
A nontrivial solution of these four linear homogeneous equations for U1, U2, U3, and Ψ exists only if the determinant of the coefficients vanishes, which yields the governing dispersion relation
(11)detG=0,
in which the elements gij(i,j=1,2,3,4) of the matrix G are
(12)[g11g12g13g14g21g22g23g24g31g32g33g34g41g42g43g44]{U1U2U3Ψ}={0000},
where
(13)g1s=ρω2[-δ1s+ε1ikεkmsΩiωΩmω-2Iε1isΩiω]+Ci1slklki,s=1,2,3,g14=eki1kkkiΨ,g2s=ρω2[-δ2s+ε2ikεkmsΩiωΩmω-2Iε2isΩiω]+Ci2slklki,s=1,2,3,g24=eki2kkkiΨ,g3s=ρω2[-δ3s+ε3ikεkmsΩiωΩmω-2Iε3isΩiω]+Ci3slklki,s=1,2,3,g34=eki3kkkiΨ,g41=-ei1lklki,g42=-ei2lklki,g43=-ei3lklki,g44=ϵijkjki.
Further, with the help of inhomogeneous wave theory [23, 25], assume that the complex wave vector can be decomposed in terms of wave propagation direction as
(14)kj=Pj+iAj=Pnj+iAmj,
where Pj is the propagation vector with its magnitude of P=PjPj, Aj is the attenuation vector with its magnitude of A=AjAj, and (nj,mj) are the unit vectors along the propagation direction (normal to the equiphase plane) and the perpendicular to the plane of constant amplitude (normal to the equiamplitude plane), respectively. Generally, nj≠mj represents an inhomogeneous wave problem while nj=mj represents a special case of a homogeneous wave problem.
Further, the unit vectors (nj,mj) can be further expressed in terms of the angle θ between nj and x3, the angle γ between nj and mj as shown in Figure 2. Via (14), we obtain
(15){n1,n2}={sinθ,cosθ}T,{m1,m2}={sin(θ+γ),cos(θ+γ)}T,njmj=cosγ.
Correspondingly, the wave vector ki can be expressed in terms of one complex number, the propagation angle θ, and the attenuation angle γ, such that,
(16)k1=Psinθ+iAsin(θ+γ),k2=Pcosθ+iAcos(θ+γ).
Inserting (16) into the dispersion equation (11) and then decomposing it into the real and imaginary parts leads to solvable equations in terms of P and A for the given attenuation angle γ, propagation angle θ, and rotation speed Ω(17)DR(P,A)=0,DI(P,A)=0,A∈0∪R+,
where DR and DI are the operators on P andA, which are nonlinear and coupled algebraic equations in terms of (P,A). According to the definitions of P and A in (14), the right solution of P and A should be real-valued. Therefore, the domain of attenuation angle γ is determined by the condition that P and A are nonnegative real numbers (here only the positive direction of wave propagation is considered). Thus, the wave propagates with the phase velocity cp and the nonnegative attenuation A, which agrees with the Sommerfeld radiation condition, that is, vanishing at infinity. Generally, there are three roots of (P,A) that are related to three elastic wave modes: one quasilongitudinal (QL) and two quasitransverse (QT 1, 2) waves for the given θ, γ, and Ω. It is noted that, due to the static electric field assumption, there is no independent wave mode in the electric field, whereas the electric wave still can propagate with the elastic wave modes via the constitutive relationship (4). After P is solved, the phase velocity can be defined as
(18)cp=ωP,
and A is the corresponding wave attenuation.
Illustration of equiphase and equiamplitude planes and exponential variation of the amplitude along the phase propagation direction.
4. Results and Discussions
In order to discuss the problem in greater detail and to find out the effects of the rotation speed Ω of the body, propagation angle θ, attenuation angle γ on the phase speed cp, and attenuation coefficient A of the inhomogeneous wave, we have computed them by taking the following piezoelectric material parameters in Table 1. All the physical constants are rewritten with the help of Voigt notation, whose rule is that the subscripts of a tensor are transformed by the rule {11→1,22→2,33→3,23→4,31→5,12→6}.
Material properties of Ba2NaNb5O15 crystal.
Cij (GPa)
C11
C22
C33
C44
C55
C66
C12
C13
C23
239
247
135
65
66
76
104
50
52
eij (C·m−2)
e31
e32
e33
e24
e31
e31
−4.0
−0.3
4.3
3.4
−0.4
4.3
ϵij
ϵ11/ϵ0
ϵ22/ϵ0
ϵ33/ϵ0
ϵ0=8.854×10-12 F/m
222
227
32
ρ (kg/m3)
5,700
For convenience, a parameter Ki can be defined as
(19)Ki=ω|Ω|,Ω=Ω1e1+Ω2e2+Ω3e3,
which is used to discuss the effects of rotation speed vector on the phase velocity and attenuation. Also the direction of rotation speed vector along x3, x1 will be considered and compared in the following. The wave frequency ω here is set to be 2π×106 1/second.
(I) The Phase Velocity. Figures 3 and 4 illustrate the phase velocity of QT1 wave when the piezoelectric body rotation about the x3 and x1 axes with varied rotation speeds and γ=0, respectively. The data show that the rising rotation speed leads to declining phase velocities. Because of the anisotropic property of piezoelectric body, the phase velocity performs differently at different propagation angles. It can be seen that there is a sharp drop in phase velocity at Ki=1; that is, the rotation speed is equal to the wave frequency; at the same time, the rotation direction influences the velocities. When Ki is below 1 or the rotation speed is more than wave frequency, the velocity slope is larger than when Ki is above 1. It is found that the attenuation angle γ almost does not influence the phase velocity.
Phase velocity of quasitransverse wave (QT1) versus propagation angle θ ranging from 0° to 360° with γ=0 and varied Ki, when Ω=Ω3e3.
Phase velocity of quasitransverse wave (QT1) versus propagation angle θ ranging from 0° to 360° with γ=0 and varied Ki, when Ω=Ω1e1.
The cases of quasitransverse wave QT2 are found in Figures 5 and 6, which is similar with QT1, except that the phase velocity of Ki=1 is larger than that of other values of Ki at most propagation angles.
Phase velocity of quasitransverse wave (QT2) versus propagation angle θ ranging from 0° to 360° with γ=0 and varied Ki, when Ω=Ω3e3.
Phase velocity of quasitransverse wave (QT2) versus propagation angle θ ranging from 0° to 360° with γ=0 and varied Ki, when Ω=Ω1e1.
Figures 7 and 8 illustrate the quasilongitudinal wave (QL) velocities along x3, x1 axes. The data suggest that there is no quasilongitudinal QL wave in the case of Ki=1 (when rotation speed is equal to wave frequency). Rather than QT1, 2, the phase velocities of Ki=0.1, 0.01 and Ki = infinite, 1000, 100, 10 are very close.
Phase velocity of quasilongitudinal wave (QL) versus propagation angle θ ranging from 0° to 360° with γ=0 and varied Ki, when Ω=Ω3e3.
Phase velocity of quasitransverse wave (QL) versus propagation angle θ ranging from 0° to 360° with γ=0 and varied Ki, when Ω=Ω1e1.
(II) The Attenuation. Instead of phase velocity, the attenuation angle γ influences the wave attenuation notably, which can be demonstrated in Figures 9–14. These figures imply that only large rotation speed or small Ki can affect wave attenuation, and attenuation angle can amplify such effect significantly. Therefore, the cases of Ki=0.01 are taken to demonstrate the attenuation angle influences.
Attenuation of quasitransverse wave (QT1) versus propagation angle angle θ ranging from 0° to 180° with γ=0°, 60°, when Ω3=2π×108 or Ki=0.01.
Attenuation of quasitransverse wave (QT1) versus propagation angle θ ranging from 0° to 180° with γ=0°, 60°, when Ω1=2π×108 or Ki=0.01.
Attenuation of quasitransverse wave (QT2) versus propagation angle θ ranging from 0° to 180° with γ=0°, 60°, when Ω3=2π×108 or Ki=0.01.
Attenuation of quasitransverse wave (QT2) versus propagation angle θ ranging from 0° to 180° with γ=0°, 60°, when Ω1=2π×108 or Ki=0.01.
Attenuation of quasilongitidinal wave (QL) versus propagation angle θ ranging from 0° to 180° with γ=0°, 60°, when Ω3=2π×108 or Ki=0.01.
Attenuation of quasilongitidinal wave (QL) versus propagation angle θ ranging from 0° to 180° with γ=0°, 60°, when Ω1=2π×108 or Ki=0.01.
Figures 9 and 10 depict the wave attenuations of QT1 wave mode. It is seen that when attenuation angle γ=0 of homogeneous wave, no attenuation is found to exist for any rotation speed; only when the attenuation angle is above zero, that is γ=60, the wave attenuation gains sharply around propagation angle θ=140 for piezoelectric body rotating around x3 and 0, 180 for rotating around x1. With the same value of rotation speed, the rotation around x1 axis shows more impact on the wave attenuation than that around x3 axis.
Turning back to QT2 wave, it is revealed that there is slight oscillating when attenuation angle γ=0. Likewise, nonzero attenuation angle plays large roles, which are shown in Figures 11 and 12.
Further, the attenuations of quasilongitudinal wave (QL) are shown in Figures 13 and 14. The largest increase in attenuation is found at propagation angle θ=0, 180 for rotating about x3 axis. The attenuation of γ=60 jumps when propagation angle is above 120, when the body rotates about x1.
5. Conclusions
In the framework of inhomogeneous wave theory, the propagation of waves in rotating piezoelectric solid has been analyzed firstly. Throughout this paper, four independent parameters are used to study the waves of rotating body, which are rotation speed value and its direction, propagation angle, and attenuation angle. The obtained results demonstrate that the rotation speed influences the wave characteristics significantly. It is found that Ki=1 is the critical point, that is, when rotation speed is equal to the wave frequency, around which the phase velocity varies substantially; at this point, no quasilongitudinal wave exists. Even though no damping parameters are taken into account, still tiny wave attenuation can be induced by the rotation speed. The large rising attenuations are attributable to attenuation angle. If there is no rotation, the attenuations of QT1, QT2, and QL waves are found to be zero. It is also noted that the wave velocity, as well as attenuation, behaves differently for different directions of rotation speed vector.
Acknowledgment
This research is partially sponsored by HPU Grant no. B2010-80 (648473), given to the first author.
HustonR. L.Wave propagation in rotating elastic media196423575576HustonR. L.In-plane vibration of spinning disks19653815191520SchoenbergM.CensorD.Elastic waves in rotating media19733111151252-s2.0-0015616722CensorD.SchoenbergM.Two dimensional wave problems in rotating elastic media19732714014142-s2.0-3425042712010.1007/BF00382503LaoB. Y.Gyroscopic effect in surface acoustic wavesProceeding of the Ultrasonics Symposium198068769110.1109/ULTSYM.1980.19748BeraR. K.Propagation of waves in random rotating infinite magneto-thermo-visco-elastic medium1998369851022-s2.0-0032204512WauerJ.Waves in rotating conducting piezoelectric media199910626266362-s2.0-003280735710.1121/1.427082DestradeM.SaccomandiG.Some results on finite amplitude elastic waves propagating in rotating media20041731–419312-s2.0-1114422683310.1007/s00707-004-0185-xAuriaultJ. L.Body wave propagation in rotating elastic media200431121272-s2.0-034576186210.1016/j.mechrescom.2003.07.002DestradeM.Surface acoustic waves in rotating orthorhombic crystals200446020426536652-s2.0-154261697610.1098/rspa.2003.1192TingT. C. T.Surface waves in a rotating anisotropic elastic half-space20044043293462-s2.0-434469563910.1016/j.wavemoti.2003.10.005AuriaultJ.-L.Acoustics of rotating deformable saturated porous media20056122352372-s2.0-2044449907210.1007/s11242-004-8214-xYangJ.A review of analyses related to vibrations of rotating piezoelectric bodies and gyroscopes20055256987062-s2.0-2384453619710.1109/TUFFC.2005.1503958SinghJ.TomarS. K.Plane waves in a rotating micropolar porous elastic solid200710270749060749072-s2.0-3534881309710.1063/1.2784973SharmaJ. N.GroverD.Body wave propagation in rotating thermoelastic media20093667157212-s2.0-6724908719010.1016/j.mechrescom.2009.03.005KumarR.RupenderR.Effect of rotation in magneto-micropolar thermoelastic medium due to mechanical and thermal sources2009414161916332-s2.0-6734912802910.1016/j.chaos.2008.07.004BiryukovS. V.SchmidtH.WeihnachtM.Gyroscopic effect for SAW in common piezoelectric crystalsProceedings of the IEEE International Ultrasonics Symposium (IUS '09)2009Rome, Italy2133213610.1109/ULTSYM.2009.5441457SharmaJ. N.GroverD.KaurD.Mathematical modelling and analysis of bulk waves in rotating generalized thermoelastic media with voids2011357339634072-s2.0-7995318692710.1016/j.apm.2011.01.014WegertH.ReindlL. M.RuileW.MayerA. P.On the Coriolis effect in acoustic waveguides2012131537943801PrasadR.MukhopadhyayS.Effects of rotation on harmonic plane waves under two-temperature thermoelasticity2012351110371055KothariS.MukhopadhyayS.Study of harmonic plane waves in rotating thermoelastic media of type III2012178824839Abd-AllaA. M.YahyaG. A.Thermal stresses in infinite circular cylinder subjected to rotation201233810591078YuanX.KuangZ.The inhomogeneous waves in pyroelectrics20103331721862-s2.0-7704912107610.1080/01495730903454728YuanX.The energy process of pyroelectric medium20103354134262-s2.0-7795111653510.1080/01495731003733094KuangZ. B.YuanX.Reflection and transmission of waves in pyroelectric and piezoelectric materials20113306111111202-s2.0-7865070849410.1016/j.jsv.2010.09.026YuanX.ZhuZ. H.Reflection and refraction of plane waves at interface between two piezoelectric media20122231225092521