1. Introduction
Piezoelectric structures such as bimorphs and unimorphs are of great research interest for a wide variety of applications including ultrasonic motors [
1], microgrippers [
2], and microcantilever biosensors [
3]. Impedance or admittance curves play a remarkable role in the study of piezoelectric structures, and electromechanical properties like resonance frequency, antiresonance frequency, quality factor, and losses are revealed in those curves. Various methods are used to obtain the impedance and admittance curves theoretically. The typical one is Finite Element Analysis (FEA) which mainly relies on FEA software like ANSYS and ABAQUS, but the software is expensive and losses in piezoelectric materials are ignored. The other method is by way of a Mason equivalent circuits model. The equivalent circuit separates the piezoelectric structure into an electrical port and two acoustic ports through the use of an ideal electromechanical transformer [
4]. It contains information derived from the mathematical statement of the structure, such as the piezoelectric constitutive equation and the motion equation [
5], and gives much faster calculation compared with FEA [
6]. Moreover, Mason’s model only uses one-dimensional assumptions and a more accurate result can be obtained [
7]. Therefore, it is more suitable to use equivalent circuits to present the electromechanical behavior of piezoelectric structures.
Mason’s equivalent circuit of piezoelectric structures has been investigated extensively since its introduction. Mason developed the Butterworth–Van Dyke (BVD) equivalent circuit from a one-port model into a three-port model by considering not only the electrical terminal but also the mechanical terminal of the piezoelectric patch [
8]. Germano [
9] used a simplified Mason’s equivalent circuit to describe the electromechanical and electroacoustic behaviors of a bimorph. Bao et al. [
10] studied the thickness mode of a bilaminar actuator. Yang et al. [
11] modeled a piezoelectric energy harvester composed of a rectangular unimorph. Wang et al. [
12] simulated a bimorph piezoelectric energy harvester with segmented electrodes. The above studies focus on the bending vibration modes of piezoelectric patches and related metal substrates; thus, the motion equations which generate the equivalent circuits are of flexure vibrations of the structures, and do not take the longitudinal vibration modes into consideration because longitudinal vibrations are not the working condition of those bimorphs and unimorphs. Nevertheless, for some particular applications, such as multimode ultrasonic motors, the longitudinal vibration modes of piezoelectric patches play a remarkable role [
13,
14,
15,
16,
17]. For those ultrasonic motors, first-order longitudinal vibration is used together with second- or fourth-order bending vibration to generate an elliptical motion locus at the drive feet. Except for the study of magnetoelectric laminated composites by Dong et al. [
18,
19], little attention has been drawn to the longitudinal vibration modes of piezoelectric structures.
Another aspect paid little attention is the loss of metal substrates in piezoelectric structures. Loss mechanisms of piezoelectric materials have been investigated for a long time and plenty of studies have been conducted. Holland and Uchino extensively described the losses in piezoelectric materials, and found that there are three losses: dielectric, elastic, and piezoelectric losses [
20,
21]. Sherrit et al. [
22] compared the KLM (Krimholtz, Leedom, and Matthae) and Mason’s equivalent circuits including the three losses. Chen et al. [
23] established an equivalent circuit composed of complex material numbers which represent the three losses. Dong et al. [
24] developed Mason’s equivalent circuits with three losses and external loads for different configurations of electrodes. Among those studies, “pure” piezoelectric ceramics are the targets when the equivalent circuits are related with losses and no metal substrates are concerned. This approach simplifies the derivation of motion equations, while for practical applications, the metal substrates are included. Therefore, it is more accurate to take the whole structure as a “composite” when investigating the properties of piezoelectric structures. In addition, since the metal substrates account for large portions of the entire structures, the elastic losses of the metal substrate should not be neglected while the losses in the piezoelectric ceramics are involved.
The aim of this paper is to establish an equivalent circuit for a beamlike piezoelectric structure in longitudinal vibration mode. In the motion equation, elastic loss in the metal substrate and three losses in the piezoelectric patch are built into the related parameters, and the longitudinal vibration mode is dealt with. Based on the equation, an equivalent circuit with four kinds of losses is derived. Using MATLAB, the impedance curve of the structure is obtained by calculating the equivalent circuit. Finally, the effectiveness of the equivalent circuit considering both elastic loss of the metal substrate and the three piezoelectric material losses is verified through experiments.
2. Loss and Motion Equation of the Piezoelectric Structure
Figure 1 shows the configuration of the piezoelectric structure, which consists a slender beam (metal substrate, the length/thickness ratio is 15.9) made of aluminum and a piezoelectric patch (PZT). The whole structure divides into four parts; Sections 0–2 represent the metal substrate. L
0 and L
1 denote the lengths of related sections in Cartesian coordinates. The nomenclature used in this paper is summarized in
Appendix B.
Dam**, which is responsible for dissipation of energy, is one of the key properties that determines the dynamic responses of structures [
25]. Dam** from the oscillatory system results in the decay of amplitude of free vibration [
26]. Losses expressed in complex numbers are among the most suitable indices for describing dam** [
27]. In this paper, we apply the superscript “*” to indicate the complex number parameters of the metal substrate and PZT. It should be noted that the superscript “*” differs from the complex conjugate. For the metal substrate, we use aluminum, and the complex Young’s modulus is defined to describe the properties of aluminum.
where
is the complex Young’s modulus of aluminum,
is the Young’s modulus,
j is the imaginary notation, and
is the loss factor of aluminum, where the subscript “
M” indicates the metal substrate (aluminum).
For piezoelectric materials, heat generation due to losses leads to the degradation of material properties; losses are a major concern for miniaturized devices with high power density. Loss in piezoelectric materials is considered to have three components: dielectric, elastic, and piezoelectric. The tangent functions with superscript “
”—
,
, and
—are used to represent “intensive” dielectric, elastic, and piezoelectric loss factors;
represents the “extensive” elastic loss factor and the subscript “
P” indicates the PZT [
21,
28]. In this paper, we use PZT5 for its relatively large losses and low quality factor. Compared with PZT4, which has small losses and a higher quality factor, the usage of PZT5 can illustrate the great impact of PZT losses when they are considered in the equivalent circuit model.
where
is the stiffness under a constant electric field,
is the dielectric constant under constant stress,
is the compliance under a constant electric field, and
is the piezoelectric constant.
A differential equation of motion is the basis for equivalent circuit establishment. Since this structure is a slender beam and we assume that the boundary condition is free–free, with no external force and only axial stress considered, the standard 3D piezoelectric constitutive equation can be reduced to a 1D form [
29]:
Losses in the piezoelectric patch and metal substrate are included in the derivation, so complex numbers are used and are denoted with a superscript “*”, as mentioned before. Here, D3 is the electric displacement in the z-direction, E3 is the electric field in the z-direction, X1i is the axial stress in the x-direction, and S1 is the axial strain in the x-direction.
Section 0 and PZT in this structure play the main role and act as the excitation source; the motion equation is derived upon it and the corresponding parameters are depicted in
Figure 2. Here,
ρ,
V,
m,
A, and
h stand for the density, volume, weight, cross-sectional area, and thickness, respectively; the related subscripts “
P” and “
M” indicate that those parameters describe the PZT and the metal substrate (aluminum), respectively; and similarly hereinafter. A driving voltage
U is applied on the electrode surface of the PZT,
U0 is the amplitude of voltage, and
E3 =
U/
hP.
Lw stands for the width of the structure.
Using the Lagrange function and the variational principle [
30], the motion equation of the piezoelectric structure under longitudinal vibration mode is obtained:
where
is the displacement along the
x-direction and is a function of both position
x and time
t.
According to the principle of composite materials, we can rewrite Equation (7) as the following:
where
is the density of the PZT and metal substrate (aluminum) as a composite material and
is the composite Young’s modulus as a complex number [
31].
3. Equivalent Circuit with Four Kinds of Losses
Equivalent circuits derive from the mathematical statement of the structure, namely, the motion equation. Although the calculation procedures of the equivalent circuit are developed in Mason’s edited book [
32], the main parameters presented here are in complex numbers. For Section 0 and PZT of this piezoelectric structure, the motion equation is expressed in Equation (8). Since we focus on the longitudinal vibration, a displacement formula for an arbitrary point in Section 0 and PZT is assumed:
where
α and
β are coefficients that will be introduced later. Here,
is the wave number rewritten as a complex number, and
where
ω is the angular frequency.
The velocity formula for an arbitrary point is
The distribution of forces and velocities in Section 0 and PZT is shown in
Figure 3, where
ν stands for velocity and F for force. The subscript “
P” represents the PZT, and “
M” the metal substrate; the other subscript “1” means that the corresponding parameter is located at the position of
x = 0, and “2” means the same but at
x =
L0.
3.1. Equivalent Circuit of PZT
First of all, we conduct the equivalent circuit of the PZT part. According to
Figure 3, we have
By inserting Equation (13) into (14) and calculating coefficients
α and
β, we get
The strain of an arbitrary point in the PZT is
The corresponding stress is
Then, the forces can be expressed as
The following is the calculation of the current flow
I*. The charge
Q* of the electrode surface is
The current flow
I* is
and we assume
Meanwhile, we use the following equations to simplify the expression of Equations (18)–(21).
In Equations (24) and (25), we collect all the real parts of those parameters together and mark them as
and
, which represent resistors (the mechanical energy consumption). The imaginary parts of those parameters are marked as
and
, which stand for reactance (the mechanical energy storage). The complete expressions of
,
,
, and
are listed in
Appendix A. In Equation (26), the real part and imaginary part of
are separated into
C01 and
C02.
By inserting Equations (24)–(26) into Equations (18)–(21) correspondingly, we have the final equations to decide the equivalent circuit of PZT, as follows:
Kirchhoff’s current law and Kirchhoff’s voltage law are applied to build the equivalent circuit of PZT, as shown in
Figure 4, where
,
.
3.2. Equivalent Circuit of Section 0 and Complete Equivalent Circuit
According to
Figure 3, and ignoring the effect of the bonding layer, we have
The strain of an arbitrary point in Section 0 is
The corresponding stress is
The forces can then be expressed as
We use the following equations to simplify the expression of Equations (31) and (32). The complete expressions of
,
,
, and
are listed in
Appendix A and their meanings are the same as introduced before.
By inserting Equations (33) and (34) into Equations (31) and (32) correspondingly, we have the final equations to decide the equivalent circuit of Section 0, as follows:
Then, Kirchhoff’s current law and voltage law are applied to build the equivalent circuit of Section 0, as shown in
Figure 5.
The equivalent circuits of remaining sections have the same derivation procedures as that of Section 0, apart from the different lengths and different directions of forces and velocities, as shown in
Figure 6.
Finally, by combining the above separate equivalent circuits together, we can acquire the complete equivalent circuit of the introduced piezoelectric structure, as shown in
Figure 7. The complete expressions of
,
,
, and
are listed in
Appendix A. This equivalent circuit represents the longitudinal vibration mode of a beam with a PZT covering part of the beam length. The principle of our new equivalent circuit differs from the conventional Mason’s equivalent circuit, for losses in the metal substrate and PZT are integrated, as mentioned in the introduction. The effectiveness of the proposed equivalent circuit will be verified through experiments.
4. Experiment
In order to verify the effectiveness of our equivalent circuit, a prototype of this piezoelectric structure was fabricated and tested using a Polytec PSV-400 laser vibrometer (PolyTec Inc., Waldbronn, Germany). The experimental setup is shown in
Figure 8. The piezoelectric structure is supposed to work at free–free boundary conditions, so it was lightly fixed using expanded polystyrene boards to avoid restriction on its longitudinal vibration. The driving signal was generated by a function generator (33210A, Keysight Technologies, Inc., Santa Rosa, CA, USA) and amplified by a power amplifier (HFVA-42, Nan**g Foneng Technology Industry Co., Ltd., Nan**g, China). The vibrometer measurement system contained a Junction Box, Vibrometer Controller, Data Management System, Scanning Head, computer, and monitor, as shown in
Figure 8.
Figure 9 shows the vibrometer measurement results of this structure. The tested edge surface is perpendicular to the
x-axis (
Figure 1) and demonstrates longitudinal vibration at a frequency of 81.4 kHz.
Aor stands for the original cross-sectional area, while
Ac and
Ae are the contraction and extension of the original cross-sectional area, respectively. This structure was also measured using an Agilent 4294A impedance analyzer (Keysight Technologies, Inc., Santa Rosa, CA, USA); the frequency tested by the impedance analyzer was 81.5 kHz, which is close to the vibrometer result of 81.4 kHz. An impedance curve was obtained and compared with the equivalent circuit result. The geometry and material parameters, including loss factors, are listed in
Table 1 and
Table 2. For the loss factor of aluminum
, since the loss factor decreased with increasing frequency [
33] and this structure was tested at high frequency (above 20 kHz), we assumed that it was equal to 1.0
10
-3 [
34,
35]. The piezoelectric ceramic we used was Haiying P-51 (Haiying Enterprise Group Co., Ltd., Wuxi, Jiangsu, China), a PZT5 ceramic. Compared with PZT4, the loss factors in
,
, and
of the P-51 are much higher; thus, the impact of losses can be clearly seen in the impedance curves in
Figure 10.
Figure 10a demonstrates that the impedance–frequency result of the equivalent circuit considering both the three PZT losses and the elastic loss of the metal substrate (aluminum loss or Al loss) has good agreement with the experimental result. It is notable that the deviations in
Figure 10 are mostly the values of impedances
ZfR and
ZfA in resonance and antiresonance frequencies, respectively, while the frequency values simulated by the equivalent circuit have high accuracy and discrepancies in resonance and antiresonance frequencies are close to 0%. As for the values of impedances, for neglect of the aluminum loss (
Figure 10b) or for neglect of the PZT losses (
Figure 10c), the discrepancies in impedance values become higher compared with the experiment results. Moreover, the neglect of both the PZT losses and the aluminum loss leads to extremely high discrepancies (
Figure 10d). The details are listed in
Table 3. In summary, the consideration of losses (both PZT losses and Al loss) in the equivalent circuit increases its accuracy and thus makes it more effective in the study of piezoelectric structures.
5. Conclusions
This paper describes an equivalent circuit of a piezoelectric structure which consists of a slender aluminum beam and a PZT patch. The longitudinal vibration mode of this structure was studied, and the elastic loss of aluminum and three PZT losses were built into the equivalent circuit. The result of the equivalent circuit, namely, the impedance curve of this structure, has good agreement with the experiment result. The typical frequencies such as resonance and antiresonance simulated by the equivalent circuit are almost the same as the experiment results, while the values of impedances in resonance and antiresonance have percent errors of 4.3678% and 8.4631%, respectively. By contrast, neglect of the elastic loss of aluminum in the equivalent circuit leads to discrepancies of 37.890% and 34.149% in values of impedances, neglect of the PZT losses leads to discrepancies of 65.839% and 182.47%, and neglect of both the PZT losses and the aluminum loss leads to discrepancies of 99.831% and 18,105%. The introduction of PZT losses and aluminum loss to the equivalent circuit dramatically contributed to the result accuracy. This result can be extended to the analysis and design of piezoelectric sensors and actuators. In this paper, only longitudinal vibration is investigated; however, in practical applications, such as in multimode ultrasonic motors, bending or torsional modes are involved. Thus, the next step of this work is to study the equivalent circuits combined with multiple vibration modes on more practical assumptions, and to increase the result accuracy of equivalent circuits containing losses.