The paper experimentally deals with the radial in-plane vibration characteristics of disk-shaped piezoelectric transducers. The radial in-plane motion, which is induced due to Poisson’s ratio in the piezoelectric disk polarized in the thickness direction, was measured by using an in-plane laser vibrometer, and the natural frequencies were measured by using an impedance analyzer. The experimental results have been compared with theoretical predictions obtained by simplified theoretical and finite-element analyses. It appears that the fundamental mode of a piezoelectric disk transducer is a radial mode and its radial displacement distribution from the center to the perimeter is not monotonic but shows maximum slightly apart from the perimeter. The theoretically-calculated fundamental frequencies agree well with the finite-element results for small thickness-to-diameter ratio, and they are accurate within 7 % error for the ratio up to 0.4.
1. Introduction
Piezoelectricity
(1)
is one of the main transduction mechanisms adopted in electromechanical sensors and actuators
(2)
. Piezoelectric transducers are widely used, especially in ultrasonic sensors
(3)
which convert mechanical vibrations into electric signals, and actuators
(4)
which convert electric signal into mechanical vibrations. In many ultrasonic devices, such as flowmeters and liquid-level meters etc., piezoelectric transducers generate and detect mechanical waves
(3)
. Piezoelectric transducers for ultrasonic devices are designed on the basis of the vibration characteristics of resonance
(5)
.
The vibration characteristics of piezoelectric transducers depend on geometric shapes and dimensions as well as material properties and boundary conditions. Most piezoelectric transducers are disk-shaped
(6~9)
or ring-shaped
(10~12)
, and their transduction mechanism is based on the motion in the thickness direction. In some cases the shape is cylindrical or spherical. Piezoelectric cylindrical or spherical transducers with radial polarization have resonance characteristics determined by their thickness and curvature
(13,14)
. Torsional transducers made of piezoelectric disks polarized in the circumferential direction have resonance characteristics determined by thickness
(15)
.
In the piezoelectric disk transducers polarized in the thickness direction, radial motion is induced due to Poisson’s ratio
(16)
. Since the radius is larger than the thickness in a thin disk, the lowest several modes are the radial modes governed mainly by radial boundary conditions. Thickness modes appear in the higher modes. This paper experimentally clarifies the radial in-plane vibration of piezoelectric disk transducers and to verify the theoretical prediction by comparing experimental results with analytical ones.
Theoretical approaches appear in much of the literatures. Reference
(6)
reported the dependence of the vibrational mode on the disk diameter- to-thickness ratio calculated by the finite-element method. Reference
(7)
presented vibration characteristics depending on the thickness-graded material properties. Reference
(8)
showed a generalized form of Hamilton’s principle for coupled electromechnical system and compared impedance curves obtained theoretically and experimentally. Reference
(9)
displayed experimental results obtained by three techniques and compared them with numerical ones. In this literature, mode shapes were shown qualitatively to show vibration distribution. Piezoelectric hollow-disks, so-called ring, were also dealt theoretically in much of literature
(10~12)
.
This paper focuses on the fundamental mode of piezoelectric disk transducers, and presents its radial in-plane vibration motion obtained experimentally by in-plane laser vibrometry. The experimentally- obtained shape of the fundamental mode, which was not available in other literatures, is clearly displayed in this paper. A simplified theoretical analysis
(17)
is briefly described to compare the experimental results with theoretical predictions. Some of the theoretical and experimental approaches were cited from other literature about ring transducers
(18)
.
2. Theoretical Prediction
In order to compare experimental results with analytical ones, theoretical prediction is carried out in advance. Two kinds of analysis are described in this section; one is theoretical analysis and another is finite-element analysis.
- 2.1 Theoretical Analysis
- (1) Problem formulation
A piezoelectric disk transducer is schematically shown in
Fig. 1
, where
l
is the thickness and
a
is the radius. The transducer has uniform electrodes on the top and bottom surfaces (
z
= 0,
l
). The electromechanical relations were well formulated, and the constitutive equations of the piezoelectric disk transducers are expressed as follows
(2)
:
Schematic diagram of a disk transducer
where T, S, D, and E are the matrices of stresses, strains, electric displacements, and electric fields, respectively. In addition, c
E
,
e
, and
εS
are the coefficient matrices of the stiffness with constant electric field, the piezoelectric stress constant, and the permittivity with constant strain, respectively.
The axisymmetric vibration in the piezoelectric disk can be formulated in terms of radial displacement
u
(
r
,
z
,
t
), axial displacement
w
(
r
,
z
,
t
), and electric potential
ϕ
(
r
,
z
,
t
) with cylindrical coordinates (
r
,
θ
,
z
) and time (
t
). Normal strains (
εr
,
εθ
,
εz
), shear strains (
γθz
,
γzr
,
γrθ
), and electric field (
Ez
) are related as follows:
Inserting Eqs. (1)~(3) into Eqs. (4) and (5) yields normal stresses (
σr
,
σθ
,
σz
) and electric displacement
Dz
as follows:
- (2) Solution
The disk transducer, the radius of which is much larger than the thickness, satisfies plane stress conditions in the thickness direction, and thus
σz
= 0. This assumption simplifies the theoretical analysis. Eq. (6c) can be rewritten to express the normal strain
εz
as follows:
Eqs. (6a,b) and (7) are rewritten as follows by inserting Eq. (8):
Here the constants with superscript
p
are as follows:
The equation of motion derived from the force equilibrium in the radial direction is
where
ρ
is the mass density. Inserting Eq. (9) into Eq. (12) yields the following governing equation:
where
in the wave Eq. (13) is the propagation speed of the wave.
When the voltage applied to the electrodes is a harmonic function of time
t
with frequency
ω
, the displacement
w
and the electric potential
ϕ
are regarded as harmonic functions of time with the same frequency. Therefore,
u
(
r
,
z
,
t
) can be expressed through the separation of variables as seen in the following:
Substituting Eq. (14) into Eq. (13) provides the following governing equations in terms of
U
(
r
) :
where
k
( =
𝜔
/
c
)is the wave number.
The solution of the Bessel Eq. (15) has the following form:
where
J
1
and
Y
1
are the Bessel functions of the first and second kinds, respectively, of order 1. The unknown constants
A
and
B
are determined according to the boundary conditions.
- (3) Radial mode characteristics
For the transducer free on the circumferential surface (
r
=
a
), the boundary conditions are the following:
In order to satisfy the boundary condition (17a),
B
= 0 in Eq. (16), and thus the solution becomes
Here
qn
( =
ka
) is the
k
times the radius
a
. Satisfying the boundary condition (17b) with Eq. (6a) in the uniform electric field yields
The natural frequencies of the radial mode are calculated from Eq. (19). In the characteristic Eq. (19), the unknown variable
qn
can be determined by using numerical tools, such as Mathematica
(19)
, and the result is
qn
= 2.07, 5.40, 8.58, … . Frequency
f
has a relation with wave number
k
and wave speed
c
as follows:
In addition to the natural frequencies, the radial mode shapes can also be calculated from Eq. (18).
The piezoelectric disk specimens considered in this paper are shown in
Fig. 2
and listed in
Table 1
. The properties of the constituent PZT-4 and PZT-5A materials are well-known, and they are listed in
Table 2
. The material properties to be inserted into the equations are obtained by converting the properties in
Table 2
. The converted properties are listed in
Table 3
.
Specimens of piezoelectric disks; (a) specimen A, (b) specimen B, (c) specimen C, (d) specimen D
Sizes and material of piezoelectric disk specimens
Sizes and material of piezoelectric disk specimens
Material properties of PZT specimens
Material properties of PZT specimens
Material properties of PZT specimens, coverted from the properties inTable 2
Material properties of PZT specimens, coverted from the properties in Table 2
The natural frequencies calculated from Eq. (20) are listed in
Table 4
. The shapes of radial modes calculated from Eq. (18) are displayed in
Fig. 3
. They will be compared with the experimental results in Section 4.
Natural frequencies of the radial modes calculated for free specimens
Natural frequencies of the radial modes calculated for free specimens
Radial mode shapes obtained by theoretical analysis ; (a) 1st mode, (b) 2nd mode, (c) 3rd mode
- 2.2 Finite-element Analysis
The theoretical analysis described in the previous section has advantages in the expressions to calculate the natural frequencies and mode shapes. However, it has disadvantages in that the analysis is based on some assumptions to simplify the real physical phenomenon. Before comparing the theoretically- calculated results with experimental ones, the theoretical analysis is enhanced by a finite-element analysis.
The natural frequencies and mode shapes are obtained by using the commercial software ANSYS. Modal analysis and harmonic analysis were carried out with the boundary conditions that all surfaces are mechanically free and the electric field is established in the thickness direction by the electrodes on the top and bottom surfaces.
Harmonic analysis results are displayed in
Fig. 4
in the form of impedance curves for the models of four specimens. In the impedance curves, the minimum points represent the resonances of the transducers. The natural frequencies obtained by the finite-element analysis are listed in
Table 5
. The fundamental frequencies are 91.0, 81.4, 49.5, and 39.6 kHz for specimens A, B, C, and D, respectively. They are compared with the theoretically- calculated values in Section 4. Incidentally, the mode shapes obtained by the finite-element modal analysis are displayed in
Fig. 5
. The results of specimen A only are shown in the figure as an example. The shapes of the first three modes in
Fig. 5
show radial modes even though other types of modes were not excluded in the analysis.
Impedance curves obtained by the finite element analysis; (a) specimen A, (b) specimen B, (c) specimen C, (d) specimen D
Natural frequencies obtained by the finite element analysis
Natural frequencies obtained by the finite element analysis
Mode shapes obtained by the finite-element analysis for specimen A; (a) 1st mode, (b) 2nd mode, (c) 3rd mode
3. Experiments
The vibration characteristics of the piezoelectric disk transducers are experimentally determined in this section, and they are compared with theoretical predictions in the next section. The specimens used in the experiments are shown in
Fig. 2
, and their sizes and materials are listed in
Table 1
.
- 3.1 Impedance Analysis
By using an impedance gain/phase analyzer (Agilent Technology 4192A), impedance curves were measured for four kinds of specimens. The measurements were repeated with three pieces of four specimens. The results of impedance curves obtained with one piece of each specimen are displayed in
Fig. 6
. In the impedance curves, the minimum points represent the resonances of the transducers. The natural frequencies experimentally -obtained by the impedance curve measurements are listed in
Table 6
, where the values were averaged from the results of three pieces. The fundamental frequencies are 88.3, 79.6, 48.4, and 38.3 kHz for specimens A, B, C, and D, respectively. They are compared with the theoretically-obtained results in Section 4.
Impedance curves obtained by experiments; (a) specimen A, (b) specimen B, (c) specimen C, (d) specimen D
Natural frequencies experimentally-obtained by the impedance-curve measurement
Natural frequencies experimentally-obtained by the impedance-curve measurement
- 3.2 In-plane Laser Interferometry
The mode shapes of the radial in-plane vibration were measured by using an in-plane laser vibrometer. The apparatus used in the experiment is based on the laser interference, and it consists of an optical sensor head (LSV-065-306F) and a controller (OFV-3320) made by Polytec. This apparatus is to measure the moving velocity of the plane perpendicular to the central line bisecting two converging laser beams as shown in
Fig. 7
. The principle of this apparatus is the Doppler phenomenon, where an optical frequency shift of the back-scattered laser light caused by surface motion is evaluated electronically.
Schematic diagram of the in-plane laser vibrometer
As shown in
Fig. 7
, the laser output from a diode is divided into two beams. These two beams progress with incident angle
ϕ
into the plane to be measured, and the superposition of the beams results in the interference pattern of parallel bright and dark fringes. The fringe spacing Δ
S
is a system constant which depends only on the incident angle
ϕ
and laser wavelength λ as follows:
The controller measures the Doppler frequency
fD
and multiplies Δ
S
. The output signal of the controller is monitored at an oscilloscope or a signal analyzer.
A signal generator (Agilent 33220A) was used to exert electric signals with a constant voltage and variable frequency onto piezoelectric disks. The output signal of the controller was monitored on an oscilloscope (Tektronix TDS3032) and the frequency of the signal resulting at the maximum amplitude was determined as the natural frequency of the disk. The fundamental frequencies of the disks measured in this way are 88.2, 79.4, 47.9, and 38.0 kHz for specimen A, B, C, and D, respectively. These results are listed in
Table 7
and will be compared with other results in Section 4.
Comparison of the fundamental frequencies obtained by calculation, finite-element method, and measurement for three specimens with free boundary conditions
Comparison of the fundamental frequencies obtained by calculation, finite-element method, and measurement for three specimens with free boundary conditions
Vibration amplitude at the natural frequency was measured by using a signal analyzer (B&K 2035). The measured amplitude is RMS voltage, which is proportional to the vibration velocity, and is converted to the vibration displacement. The measurement is carried out along a radius from the center to the perimeter at every 1 mm. The measured vibration amplitude was normalized to the maximum amplitude and displayed in
Fig. 8
for four specimens. In
Fig. 8
, the solid line is the theoretically-calculated result obtained in Section 2.1. The mode shapes displayed in
Fig. 8
are discussed in Section 4.
Radial mode shapes of the fundamental mode measured by the laser vibrometer; (a) specimen A, (b) specimen B, (c) speimen C, (d) specimen D
4. Comparisons and Discussions
Natural frequencies and mode shapes obtained theoretically in Section 2 and experimentally in Section 3 are compared with each other and discussed in this section.
The fundamental frequencies obtained by calculation, the finite-element method, and measurement for four specimens are listed in
Table 7
and displayed in
Fig. 9
. The experimental results agree well with the theoretical predictions. The fundamental frequency of the piezoelectric disk transducer depends on its radius, and the frequency is higher as the radius is smaller. The theoretically- calculated fundamental frequencies are compared with the finite-element results in
Fig. 9
. The theoretically-calculated results agree well with the finite-element results for small thickness- to-diameter ratio, and they are accurate within 7 % error for the ratio up to 0.4.
Comparison of the fundamental frequencies obtained by calculation, finite-element method, and measurements; (a) PZT-4, (b) PZT-5A
The mode shapes obtained by the laser vibrometer are compared with the theoretical predictions in
Fig. 8
. The results agree well with each other. The displacement distribution from the center to the perimeter is not monotonic. The vibration amplitude is maximum slightly apart from the perimeter.
5. Conclusion
The paper presented the radial in-plane vibration characteristics of disk-shaped piezoelectric transducers polarized in the thickness direction. The radial in-plane vibration characteristics of piezoelectric disk transducers were experimentally investigated by measuring natural frequencies and mode shapes. The natural frequencies of piezoelectric transducers were measured by using an impedance analyzer and compared with theoretical predictions. The radial in-plane motion, which was induced due to Poisson’s ratio in the piezoelectric disk polarized in the thickness direction, was measured by using an in-plane laser vibrometer and compared with theoretically obtained results.
The fundamental mode shape of a piezoelectric disk is a radial mode and the radial displacement distribution from the center to the perimeter is not monotonic but shows maximum slightly apart from the perimeter. The fundamental frequencies obtained by the simplified theoretical analysis briefly described in this paper well agreed in small thickness-to-diameter ratio range with other results obtained by finite-element analysis and experiments. In the range of the ratio over 0.4, the theoretically-obtained natural frequencies showed errors more than 7 %. The theoretical analysis may be used to predict the vibration characteristics of thin piezoelectric disk transducers.
‡ Recommended by Editor SungSoo Na
Acknowledgements
This research was supported by the MSIP (Ministry of Science, ICT & Future Planning), Korea, under the Convergence-ITRC (Convergence Information Technology Research Center) support program (NIPA-2014-H0401-14-1005) supervised by the NIPA (National IT Industry Promotion Agency).
BIO
Dae Jong Kim received the B.S. degree in mechanical engineering from Soongil University in 2013. During his stay at Soongsil as a graduate student, he has been working on the vibration characteristics of piezoelectric transducers.
Se Hwan Oh received the B.S. and M.S. degrees in mechanical engineering from Soongil University in 2011 and 2013, respectively. During his stay at Soongsil as a graduate student, he has been working on the vibration characteristics of piezoelectric transducers. Since 2013 he has been working in Hyundai Wia Co. Ltd.
Jin Oh Kim received the B.S. and M.S. degrees in mechanical engineering from Seoul National University in 1981 and 1983, respectively, and the Ph.D. degree from University of Pennsylvania in 1989. Since 1997, he has been with the Faculty of Soongsil University, where he is currently a Professor of mechanical engineering. His research interests are in the areas of ultrasonic sensors and actuators using mechanical vibrations and human-body dynamics for rehabilitation robots.
Ikeda T.
1996
Fundamentals of Piezoelectricity
Oxford University Press
Oxford
Busch-Vishniac I. J.
1999
Electromechanical Sensors and actuators, Chapter 5
Springer
New York
Lynnworth L. C.
1989
Ultrasonic Measurements for Process Control
Academic Press
Boston
Choi S.-B.
,
Han Y.-M.
2010
Piezoelectric Actuators : Control Applications of Smart Materials
Taylor & Francis
Boca Raton
Miu D. K.
1993
Mechatronics: Electromechanics and Contromechanics, Chapter 6
Springer
New York
Kunkel H. A.
,
Locke S.
,
Pikeroen B.
1999
Finite-element Analysis of Vibrational Modes in Piezoelectric Ceramic Disks
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
37
316 -
328
Lee P. C. Y.
,
Yu J.-D.
,
Li X.
,
Shih W.-H.
1999
Piezoelectric Ceramic Disks with Thickness-graded Material Properties
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
46
205 -
215
Ho S.-T.
2007
Modeling of a Disk-type Piezoelectric Transformer
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
54
2110 -
2119
Lin Y.-C.
,
Ma C.-C.
2004
Experimental Measurement and Numerical Analysis on Resonant Characteristics of Piezoelectric Disks with Partial Electrode Designs
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
51
937 -
947
Iula A.
,
Lamberti N.
,
Pappalardo M.
1996
A Model for the Theoretical Characterization of Thin Piezoceramic Rings
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
43
370 -
375
Li H. L.
,
Hu J. H.
,
Chan H. L. W.
2004
Finite Element Analysis on Piezoelectric Ring Transformer
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
51
1247 -
1254
Ho S.-T.
2007
Modeling and Analysis on Ring-type Piezoelectric Transformers
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
54
2376 -
2384
Kim J. O.
,
Lee J. G.
2007
Dynamic Characteristics of Piezoelectric Cylindrical Transducers with Radial Polarization
Journal of Sound and Vibration
300
241 -
249
Kim J. O.
,
Lee J. G.
,
Chun H. Y.
2005
Radial Vibration Characteristics of Spherical Piezoelectric Transducers
Ultrasonics
43
531 -
537
Kim J. O.
,
Kwon O. S.
2003
Vibration Characteristics of Piezoelectric Torsional Transducers
Journal of Sound and Vibration
264
453 -
473
Meitzler A. H.
,
O’Bryan H. M.
,
Tiersten H. F.
1973
Definition and Measurement of Radial Mode Coupling Factors in Piezoelectric Ceramic Materials with Large Variations in Poisson’s Ratio
IEEE Transactions on Sonics and Ultrasonics
SU-20
233 -
239
Kim J. O.
,
Oh S. H.
,
Piao C.
2013
Radial Mode Vibration Characteristics of Disc-type Piezoelectric Transducers
Proceedings of the 15th Asia Pacific Vibration Conferences
1058 -
1061
Piao C.
,
Kim J. O.
2014
In-plane Vibration Characteristics of Piezoelectric Ring Transducers
Transactions of the Korean Society for Noise and Vibration Engineering
24
(10)
780 -
787
Wolfram S.
1999
The Mathematica Book
4th ed.
Wolfram Media Inc.
Champaign