1. Introduction
The electromechanical coupling coefficient indicates the degree to which the piezoelectric element is able to convert the input energy
into the output energy
. When the element is in actuator mode, the input energy is electrical energy, while the output energy is mechanical energy, and vice versa in sensor mode. The magnitude of the electromechanical coupling coefficient can be determined either theoretically based on knowledge of the material constants (static value), or calculated from the measured series and parallel resonance frequencies (dynamic value) [
1].
In a simple way, the properties of a resonator can be described by an equivalent electrical circuit with lumped parameters. It describes the resistor impedance in the vicinity of resonance and antiresonance frequencies [
2]. From the equivalent electrical circuit of a piezoelectric resonator, the effective electromechanical coupling coefficient
, [
2], can be determined. This coefficient may be expressed independently of the mode of oscillation as
where
,
are the series and parallel resonance frequencies expressed from the electrical equivalent circuit, and
and
are the resonance and antiresonance frequencies meeting the requirement of zero imaginary impedance component.
It is possible to show the relationship between the effective electromechanical coupling coefficient and the electromechanical coupling coefficients of certain types of oscillations of resonators of certain shapes [
2], an example being the thickness mode of oscillation of circular fully electroded resonators. In a general case, however, the dependence cannot be simply defined and one must make do with the general relation (
1). This is also the case for circular resonators with wrap-around electrodes. Because of the axially asymmetrical electrodes, the magnitude of the electromechanical coupling coefficient cannot be expressed analytically. As shown in [
3], the wrap-around electrodes cause inhomogeneous polarization in the resonator; moreover, in the area where the electrodes are wrapped-around, polarization is not achieved at all. It can be therefore expected that the electromechanical coupling coefficient will be lower for resonators with wrap-around electrodes than for resonators of the same dimensions with fully electroded bases. In the presented work, we will verify this hypothesis computationally and validate it by comparison with experimental results.
In order for the obtained results for resonators with fully electroded bases to be comparable with the results calculated for resonators with wrap-around electrodes, we will use the relation (
1) to calculate the coefficient of electromechanical coupling. It will therefore be necessary to determine the resonance and antiresonance frequencies of the oscillation for which we determine
. The vibration to be studied will be the fundamental thickness oscillation.
The study of the resonance characteristics of piezoelectric resonators with incomplete or asymmetric electrodes is already the subject of [
4,
5], or more, recently [
6,
7]. In the publications [
8,
9], the analytical solution is compared with theexperimental results. However, in all these publications, the authors focus on resonators where a certain form of symmetry can be seen in the electrode surfaces, e.g., circular electrodes with different diameters but a common axis of circles, or rectangular electrodes on a resonator of rectangular plan, where two planes of symmetry can be found.
In the available literature, the author has not found a relevant study concerning circular, partially electroded resonators and effects of this feature on the resonance frequency of the thickness oscillation. Neither for resonators electroded symmetrically according to the axis of symmetry of the disc, nor for asymmetrical electroding similar to the wrap-around electrodes. Neither experimental nor theoretical study has been found. However, mention should be made of the work of [
10], where incomplete or axially asymmetric electroding is studied, but not for thickness oscillation and not in the electrode configurations usable for comparison with wrap-around electrodes. Similar results are reported in [
11], where the application of symmetrical electroding on a circular resonator in order to determine the Poisson number of piezoelectric ceramics is presented. However, the electrode variants and types of excited oscillations presented in his work are not relevant to the case studied here.
In publications [
12,
13,
14,
15], the finite element method is used as a tool to determine the resonance frequency, antiresonance frequency, and electromechanical coupling coefficient. The change in the electromechanical coupling coefficient with the size of the area covered by the electrode in a circular resonator is also examined, but only configurations with electrodes symmetric to the resonator symmetry and almost exclusively for radial oscillation are studied. None of the works mentioned above deals with the issue of changing the electromechanical coupling coefficient of resonators with electrodes not symmetrical with respect to the axis of symmetry of the resonator; resonators with a wrap-around electrode can be included among them. At the same time, the thickness mode of the resonator oscillation is considered in none of these works. The presented work is therefore unique in its analysis of the influence of asymmetric electrodes on the electromechanical coupling coefficient of the resonator thickness oscillation.
2. Materials and Methods
Consider the piezoelectric continuum with density
, characterized by material tensors
. We denote the volume of the continuum as
and its boundary as
. Behavior of the piezoelectric continuum is in time range
governed by two differential equations [
16]. Gauss’s law (
2) and Newton’s law of motion (
3),
where
T,
D,
S,
E,
u are stress, the vector of electric displacement, strain, electric field, mechanical displacement. Initial and boundary conditions are added:
The right-hand sides
,
,
, and
represent excitation of the continuum. If we discretize the Equations (
2) and (
3) by the finite element method, this set of equations will have the form [
16,
17],
where
is the stiffness matrix,
is the mass matrix,
is the piezoelectric matrix,
is the electric matrix,
is the displacement vector,
is the vector of electric potentials,
is the vector of external forces, and
is the electric charge vector. When a piezoelectric resonator is excited by an electric voltage applied to its electrodes, the determining quantity is the difference in electric potentials between the electrodes. According to [
16,
17], we shall consider one of the electrodes grounded; the vector of electric potentials corresponding to the nodes of the finite element mesh at the location of the grounded electrode is
the vector of electric potentials corresponding to the nodes of the mesh at the places of the ungrounded electrode is denoted by
. The system (
10) can be written as
where
is the vector of electric charges in the nodes belonging to the ungrounded electrode, and
is the vector of electric potentials in the nodes outside the electrode. By eliminating
from the second equation of the system (
12) we obtain
where
At resonance, we assume [
16,
17], that the electrical impedance takes the minimum value and the electrodes are short circuited (that is, both electrodes have the same electrical potential), i.e.,
. In resonance we can regard
, the system (
13) can be modified to the form
The resonance frequencies are then obtained as solutions of a generalized eigenvalue problem in the form
When calculating the antiresonance frequency, we use the system of Equation (
13) again. In accordance with [
18,
19], we consider disconnected electrodes, thus the electric charge on the ungrounded electrode meets the condition
However, the part of the system of equations related to the grounded electrode was eliminated due to the condition (
11), see the transition from (
10) to (
12). At the same time,
will apply; the system (
13) can be written as
From the second equation we express
We then obtain the antiresonance frequency as a solution to the generalized problem of eigenvalues in the form
where
Knowing the values of the resonance frequencies, calculated according to the relation (
22), and the values of the antiresonance frequencies according to the relation (
22), the values of the effective electromechanical coupling coefficients can be obtained by substituting in the relation (
1).
2.1. Task Assignment
We calculate of the fundamental thickness oscillation for resonators with wrap-around electrodes, where we expect inhomogeneous polarization due to electrodes being wrapped around. We compare the result with the calculation of thickness oscillation of a resonator of identical dimensions but fully electroded on both circular surfaces, where homogeneous polarization can be expected throughout the resonator volume. By comparing the calculated of the resonators with wrap-around electrodes with the of the fully electroded resonators, we conclude the effect of polarization on the electromechanical coupling coefficient.
This approach is applied when computing
of two types of piezoelectric transducers, labeled (a) and (b), shown in
Figure 1. The transducers differ in the shape of the electrode wrap-around. In case of (a), the overlay is in the shape of a circular segment, in case of (b), it is in the shape of intersecting circles.
To make the results comparable, we will consider both piezoelectric transducers of the same size. The criterion of equivalence for the wrap-around electrodes needs to be established. We will consider equivalent the electrodes that will have the same area of wrap-around electrode in both cases.
Let
, where
d is the diameter of the disc and
t its thickness. We will consider a disc with a dimensional ratio of
, where the basic thickness mode of the oscillation is significantly distinct. In the monograph [
1] (chapter 5.3), the frequency spectrum of a resonator with entire electrodes on both circular surfaces with ratio
, made of NCE51 piezoelectric ceramics, is analyzed. We will thus work with a resonator that maintains the dimensions of resonator from chapter 5.3 of the monograph [
1], except that it will not be fully electroded on both circular surfaces. A diagram of the dimensions of the transducers is shown in
Figure 2. The dimensions of the transducers are given in
Table 1.
The dimensions of the wrap-around are chosen so that the area of the wrapped around electrodes is identical in both cases. The ratio of the wrap-around to the total circular area is approximately .
We consider the transducers made of NCE51 piezoelectric ceramics, again in accordance with Chapter 5.3 in [
1]. The values of the elements of the above-mentioned tensors are given in
Table 2 and are taken from [
1,
20]. In the table, permittivity is given as relative permittivity.
The polarization of the above-mentioned resonators with a wrap-around electrode is analyzed in the work by [
3]. We will consider the polarization directions as calculated in [
3]. Thus, we know the polarization directions of piezoelectric ceramics and consequently the macroscopic piezoelectric properties that are determined by polarization.
2.2. Geometry and Finite Element Mesh
The electromechanical coupling coefficient is calculated using the relation (
1), i.e., based on the knowledge of the resonance and antiresonance frequency of the respective oscillation. The resonance and antiresonance frequencies are determined by the procedure in the introduction of
Section 2. The thickness oscillation that is symmetric to the resonator symmetry axis will be examined. The resonator has the cylindrical shape, so the axis of symmetry is the axis of the cylinder representing the resonator. It is therefore appropriate that the chosen discretization of the resonator area reflects the symmetry of the oscillation and is therefore homogeneous. The use of inhomogeneous grids may lead to:
We can say that this is a form of discretization error. By using a finite element mesh, the edges of which will have a comparable size in the whole discretized area, we eliminate the mentioned errors in the obtained solution. At the same time, identical finite element meshes will be used to calculate the resonance and antiresonance frequencies of resonators with wrap-around electrodes and to calculate the resonance and antiresonance frequencies of a fully electroded resonator. The use of identical meshes has the following benefits:
The same discretization error can be expected for the calculation of the resonance frequency of a resonator with fully electroded circular surfaces, as well as for resonators with wrap-around electrodes;
A very similar numerical error in the calculation of resonance frequencies can be expected for all analyzed cases.
Under these circumstances, the results of the calculation of the resonance frequencies obtained for each case can be compared. The differences between the results will be negligible because:
Resonators fully electroded on circular surfaces and resonators with wrap-around electrodes have the same diameter and thickness;
Both resonator types are made of the same material and its density determining the mass matrix is the same;
With respect to the condition (
11) the boundary conditions are also the same.
The only difference between the analysed variants is due to the anisotropy of the elastic coefficients tensor. In the case of a resonator fully electroded on circular surfaces, the homogeneous polarizing electric field results in an almost perfect orientation of the individual ceramic grains in the direction of the applied electric field, i.e., in the direction perpendicular to the thickness of the resonator.
In the case of resonators with wrap-around electrodes, we will consider the orientation of the ceramic grains in the directions calculated in the work [
3]. Thus, the difference of the individual analyzed variants will be apparent only in those elements of the stiffness matrix
, which will correspond to the finite elements in which the material properties differ.
The thickness oscillation we are analyzing is symmetric along the resonator axis. Therefore, the geometry corresponding to the symmetric half of the resonator can be used for the calculation, taking into account the boundary conditions
in the plane of symmetry. The finite element meshes used are shown in
Figure 3. They have a sufficient number of nodes for the calculation of the fundamental thickness oscillation over the base diameter, both variants (a) and (b) contain a total of about 14,000 nodes.
4. Discussion
It can be hypothesized that for the resonators with the wrap-around electrode in the shape and dimensions according to the scheme in
Figure 2 there will be an increase in the resonance frequency compared to the fully electroded resonators of the same dimensions. The differences between a fully electroded resonator and a resonator with the wrap-around electrode, which are also very plausible causes of the increase in resonance frequency, are:
Inhomogeneous polarization, especially in the area of electrode wrap-around, resulting in inhomogeneous material properties;
Inhomogeneous distribution of the excitation electric field in the resonator.
One may also ask why in the case of the mesh, or—more precisely—of the wrap-around shape according to variant (a), the growth of the resonance frequency is about 0.2% higher than in the case of the wrap-around according to variant (b). A discretization error cannot be ruled out even though the difference in the calculated resonance frequency between variants (a) and (b) is less than an order of magnitude lower than a tenth of a percent in the case of the calculated resonance frequency.
However, we can hypothesize that the difference in the calculated frequency is caused by the different amount of polarized volume of the transducer in the direction of its thickness. This, in accordance with the results given in [
3], is
This reasoning is consistent with the conclusions made, for example, in [
12]. The effect of electrode size on the resonance frequency of symmetrically electroded discs is studied there. The increase in resonance frequency with decreasing electrode area is shown. Compared with the problem of the resonance frequency of a thickness oscillation resonator with a wrap-around electrode studied here, there are two major differences:
In [
12], axially symmetrically electroded discs are studied, in contrast to resonators with a wrap-around electrode, which are symmetric only in one plane;
In [
12] the radial oscillations of the disc are studied, in contrast to the thickness oscillations studied in resonators with a wrap-around electrode.
For longitudinal oscillations of piezoelectric rods and thickness oscillations of plates, a relation for the electromechanical coupling coefficient valid for this type of oscillations [
1,
2], can be derived. The condition is full electroding of the opposite circular bases. For the thickness oscillations of the plate with full electrodes [
2], the following applies
As mentioned in the introduction to this chapter, the relationship between
and the electromechanical coupling coefficient of a given type of oscillation can be expressed. In the case of the thickness oscillation and its electromechanical coupling coefficient
, the following is true [
2],
For comparison, let us calculate the coefficient
for the fully electroded resonator as well as for resonators with a wrap-around electrode. The calculation will be performed, although we know that for resonators with a wrap-around electrode it will only be approximately valid. The results of the calculation are summarized in
Table 3.
The calculated values of are higher than the corresponding values of in all of the considered cases, but the trend is the same for both quantities. In the case of resonators with wrap-around electrodes, is lower than in the case of a fully electroded resonator.
5. Conclusions
Based on the above, it can be hypothesized that for resonators with a wrap-around electrode in the shape and dimensions according to the scheme in
Figure 2 there will be a decrease in the electromechanical coupling coefficient compared to a fully electroded resonator of identical dimensions. The hypothesis is based on the comparison of the results obtained by the same method, which is the FEM model calculation for both the resonators with a wrap-around electrode and the fully electroded resonator. The reasons that lead to the decrease in the electromechanical coupling coefficient are likely to lie in the differences between the fully electroded resonator and the resonators with a wrap-around electrode; in particular, the uneven polarization of the resonator due to the shape and placement of the electrodes in resonators with wrap-around electrodes.
In the available literature, the author has not found a study of the effect of partial electroding, asymmetrical along the resonator axis, on the electromechanical coupling coefficient. Mention may be made of the works of [
8,
12,
14,
15,
21]. Here, however, the cases of partial electroding symmetrical along the axis of symmetry of the disc and for the case of radial oscillations are studied.
It can be concluded that the claim of a decrease in the electromechanical coupling coefficient during partial electrode position is consistent with the results published in [
12]. In this work, the graph in Figure 10 shows the dependence of the electromechanical coupling coefficient on the size of the electrodes. However, these findings apply to a symmetrically electroded disc and radial type of oscillation. Thus, the hypothesis made in the present paper cannot be based on the results published in [
12].
The results presented here can be verified experimentally; for the reasons given in
Section 1 it is not possible to determine an analytical solution. The experimental determination of the frequency spectra of thickness oscillating resonators with a wrap-around electrode is the subject of the technical report of [
22]. The report contains the values of the electromechanical coupling coefficient calculated from the measured resonance and antiresonance frequencies of resonators with double-sided electroding and resonators with a wrap-around electrode. The electromechanical coupling coefficient of the thickness oscillation is determined here to be lower in most cases for resonators with a wrap-around electrode than for the fully electroded resonators. Although the measurements were carried out on a small number of samples, this result supports the proposed hypothesis of a decrease in the electromechanical coupling coefficient of resonators with a wrap-around electrode.