Next Article in Journal
Abnormal Ship Behavior Detection Based on AIS Data
Previous Article in Journal
Identification of 3D Lip Shape during Japanese Vowel Pronunciation Using Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FEM Analysis of the Effect of Polarization on the Electromechanical Coupling Factor of Resonators with a Wrap-Around Electrode

Faculty of Mechatronics, Informatics and Interdisciplinary Studies, Technical University of Liberec, Studentská 2, 46117 Liberec, Czech Republic
Appl. Sci. 2022, 12(9), 4634; https://doi.org/10.3390/app12094634
Submission received: 8 April 2022 / Revised: 27 April 2022 / Accepted: 28 April 2022 / Published: 5 May 2022

Abstract

:

Featured Application

The published results show a difference in the size of the electromechanical coupling coefficient in piezoceramic resonators with wrap-around electrodes, compared to resonators with fully electroded bases. The results can be applied in the design of new types of piezoceramic resonators.

Abstract

The efficiency of electromechanical conversion in piezoelectric elements is assessed according to the magnitude of the electromechanical coupling coefficient. For power applications of piezoelectric elements, it is desirable that the efficiency of electrical to mechanical energy conversion be as high as possible. In case of resonators with electrodes incompletely covering their bases, an inhomogeneous electric field is generated, which results in inhomogeneous polarization of the resonator. The resonator will be polarized in some places either in a direction other than the desired one or not polarized at all. The degree of polarization of the resonator is one of the factors affecting the electromechanical coupling coefficient. The aim of this work is to analyze the magnitude of the electromechanical coupling coefficient of resonators with electrodes incompletely covering the bases and to compare the results with the electromechanical coupling coefficient of resonators with fully electroded bases. The physical description is given by the linear piezoelectric equations, the Gaussian equation for the electric field, and by Newton’s law of force. On this basis, a FEM model is developed and used to analyze the electromechanical coupling coefficient. The result of the calculation of the electromechanical coupling cofficient is that for resonators with a wrap-around electrode in studied shape and dimensions there will be a decrease in the electromechanical coupling coefficient compared to a fully electroded resonator of identical dimensions. The presented conclusions are compared with analytically and experimentaly achieved results.

1. Introduction

The electromechanical coupling coefficient indicates the degree to which the piezoelectric element is able to convert the input energy W i n into the output energy W o u t . 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 k e f f , [2], can be determined. This coefficient may be expressed independently of the mode of oscillation as
k e f f 2 = f P 2 f S 2 f P 2 f a 2 f r 2 f a 2 ,
where f S , f P are the series and parallel resonance frequencies expressed from the electrical equivalent circuit, and f r and f a 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 k e f f . 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 c , e , ε . We denote the volume of the continuum as Ω and its boundary as Γ . Behavior of the piezoelectric continuum is in time range 0 , τ governed by two differential equations [16]. Gauss’s law (2) and Newton’s law of motion (3),
· D = x k ( e k i j · S ij + ε k j S · E j ) = 0 ,
· T = ϱ 2 u i t 2 = x j ( c i j k l E · S kl e k i j · E k ) i = 1 , 2 , 3 , x Ω , t 0 , τ ,
where T, D, S, E, u are stress, the vector of electric displacement, strain, electric field, mechanical displacement. Initial and boundary conditions are added:
φ ( 0 ) = φ 0
u ( 0 ) = u 0
φ = φ D on Γ 1 ,
D k n k = D N on Γ 2 ,
u i = u i D i = 1 , 2 , 3 , on Γ 1 ,
T ij n j = t i N i = 1 , 2 , 3 , on Γ 2 .
The right-hand sides φ D , D N , u i D , and t i N 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],
M 0 0 0 U ¨ 0 + K P P T E U Φ = F Q ,
where K is the stiffness matrix, M is the mass matrix, P is the piezoelectric matrix, E is the electric matrix, U is the displacement vector, Φ is the vector of electric potentials, F is the vector of external forces, and Q 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
Φ b = 0 ,
the vector of electric potentials corresponding to the nodes of the mesh at the places of the ungrounded electrode is denoted by Φ n . The system (10) can be written as
M 0 0 0 0 0 0 0 0 U ¨ 0 0 + K P e P n P e T E e P en P n T P en T E n U Φ e Φ n = F 0 Q n ,
where Q n is the vector of electric charges in the nodes belonging to the ungrounded electrode, and Φ e is the vector of electric potentials in the nodes outside the electrode. By eliminating Φ e from the second equation of the system (12) we obtain
M 0 0 0 U ¨ 0 + K ¯ P ¯ P ¯ T E ¯ U Φ n = F Q n ,
where
K ¯ = K P e E e 1 P e T ,
E ¯ = E n P en T E e 1 P en ,
P ¯ = P n P e E e 1 P en .
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., Φ n = 0 . In resonance we can regard F = 0 , the system (13) can be modified to the form
M 0 0 0 U ¨ 0 + K ¯ P ¯ P ¯ T E ¯ U 0 = 0 Q n .
The resonance frequencies are then obtained as solutions of a generalized eigenvalue problem in the form
K ¯ λ M U = 0 .
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
Q n = 0 .
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, F = 0 will apply; the system (13) can be written as
M 0 0 0 U ¨ 0 + K ¯ P ¯ P ¯ T E ¯ U Φ n = 0 0 .
From the second equation we express
Φ n = E ¯ 1 P ¯ T U .
We then obtain the antiresonance frequency as a solution to the generalized problem of eigenvalues in the form
K ^ λ M U = 0 ,
where
K ^ = K ¯ P ¯ E ¯ 1 P ¯ T .
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 k e f f 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 k e f f 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 k e f f of the resonators with wrap-around electrodes with the k e f f of the fully electroded resonators, we conclude the effect of polarization on the electromechanical coupling coefficient.
This approach is applied when computing k e f f 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 α = d / t , where d is the diameter of the disc and t its thickness. We will consider a disc with a dimensional ratio of α 1 , 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 α = 23.7 , 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 1 / 25 .
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:
  • An axially asymmetrical shape of the oscillation that does not reflect reality;
  • Inaccuracies in the calculation of the resonance frequency.
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 M 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 K , 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
D k n k = 0 ,
T ij n j = 0 ,
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.

3. Results

3.1. Calculation of Resonance and Antiresonance Frequencies for Fully Electroded Resonators, Model Verification

The resonance frequency of the thickness oscillation will be first determined for a resonator whose thickness and diameter are given in Table 1, with fully electroded circular surfaces. The result will be compared with the analytical calculation and the experimentally determined value. The calculation will be performed on the finite element meshes shown in Figure 3. We will compare the results obtained on both types of finite element meshes. The material properties are given in Table 2. We consider the polarization direction to be homogeneous throughout the resonator volume, in the direction of the resonator thickness, i.e., along the z axis, as shown in Figure 3. The calculated frequency spectrum contains the oscillation mode shown in Figure 4.
Half of the resonator cross section is shown, and the center of the resonator is on the left. The shape of the oscillation corresponds to the shape of the thickness oscillation in Section 5.3 in [1]. The calculated resonance frequencies are
f r a = 988 , 961 Hz   for   mesh   variant   ( a ) , f r b = 988 , 915 Hz   for   mesh   variant   ( b ) .
The calculated antiresonance frequencies are
f a a = 1 , 115 , 829 Hz   for   mesh   variant   ( a ) , f a b = 1 , 115 , 712 Hz   for   mesh   variant   ( b ) .
According to [1], the resonance frequency of the fundamental thickness oscillation can be expressed as
f r t = η 2 π c 33 D ϱ ,
where η is the wavenumber, c 33 D is the element of the elastic modulus tensor at constant electric induction, and ϱ is the density. The wavenumber η must satisfy the relevant boundary conditions on the circular surfaces of the resonator and the transcendental equation applies to it:
tg ( η b ) = η t 2 k t 2 .
The electromechanical coupling coefficient for the thickness mode of the oscillation can be expressed as
k t = e 33 c 33 D ε 33 S .
Let us reiterate that the polarization is assumed to be in the direction of the resonator thickness, i.e., in the direction perpendicular to the circular surfaces. Therefore, the components of material tensors with coefficients 33 appear in the relation (28). The resonance frequency of the fundamental mode of the thickness oscillation calculated for said resonator is
f r t = 994 , 270 Hz .
The difference between the value of the resonance frequency f r t calculated analytically and the numerically obtained values of f r a and f r b is 0.6%, which can be considered a very good agreement. We can conclude that the numerical model is thus verified.
Using a similar procedure as with the resonance frequency, the antiresonance frequency can be analytically calculated, see [1]. The value is
f a t = 1 , 121 , 000 Hz .
It can be noted that the values calculated by the finite element method (FEM) model are 0.6% lower in the case of both resonance and antiresonance frequencies. Since the calculation of both values is identical in principle and the same values of material parameters are used, a similar difference can be expected between the values obtained analytically and numerically.
The experimentally determined value of the resonance frequency reported in [1] is
f r e = 1 , 001 , 600 Hz
Compared to f r a , f r b , it differs by approximately 1.3%. Let us now compare the difference between the frequencies f r a , f r b that were calculated on different meshes, as documented in Figure 3. The difference is in the order of 0.001%. The calculation performed on both variants of the meshes yields comparable results. It can be thus concluded that we have a verified model through which we can capture the trends in the resonance frequency changes of the thickness oscillation caused by the electrodes being wrapped around, respectively by the inhomogeneous polarization of the resonator. At the same time, it can be stated that the different discretization of variants (a) and (b) has only a negligible effect on the calculated resonance frequency of the thickness oscillation.
The experimentally determined value of the antiresonance frequency reported in [1] is
f a e = 1 , 112 , 500 Hz .
Compared to f a a , f a b differs by approximately 0.3%.

3.2. Calculation of Resonance and Antiresonance Frequencies of Resonators with Wrap-Around Electrodes

The calculation of the resonance frequency of resonators with a wrap-around electrode will be based on the assumptions made in Section 2.2. In particular, identical finite element meshes will be used for the calculation so that the same discretization error can be achieved and the results from Section 3.1 can be compared with the results obtained for resonators with wrap-around electrodes. The meshes used will be the finite element meshes shown in Figure 3, according to the electrode structure (a) or (b) as shown in Figure 2.
The material properties listed in Table 2 will be used. The polarization and its direction, i.e., its influence on the orientation of the material tensors, are taken from [3]. In this way, the influence of different polarization in the parts with the wrap-around electrode compared to the homogeneous polarization in the whole transducer in the case of a fully electroded resonator is taken into account.
Let us first pay attention to the calculated shape of the oscillation. In the case of a fully electroded resonator, the thickness oscillation is rotationally symmetric along the axis of symmetry of the resonator. For resonators with wrap-around electrodes, the thickness oscillation is symmetric only along the plane of symmetry of the electrodes. The damping of the oscillation can be seen in the area under the wrap-around electrode; see Figure 5, where the oscillation for wrapping-around according to variant (a) is shown. The damping of the oscillation for variant (b) is similar.
The calculated resonance frequencies are
f r a = 1 , 006 , 413 Hz   for   the   electrode   structure   ( a ) , f r b = 1 , 004 , 792 Hz   for   the   electrode   structure   ( b ) .
Compared to the fully electroded resonator, there is an increase in the value of the resonance frequency in both variants of the mesh, namely
  • By 1.8% in the case of the electrode structure (a);
  • By 1.6% in the case of the electrode structure (b).
The calculated antiresonance frequencies are
f a a = 1 , 126 , 303 Hz   for   the   electrode   structure   ( a ) , f a b = 1 , 125 , 329 Hz   for   the   electrode   structure   ( b ) .
and conclusions similar to those for resonance frequencies can be drawn for them.

3.3. Calculation of the Effective Coefficient of Electromechanical Coupling

As mentioned in Section 2.2, the reference element is a circular fully electroded resonator with dimensions according to Table 1.

3.3.1. Calculation of k e f f for Fully Electroded Resonators

We calculate k e f f for fully electroded resonators from the resonance and antiresonance frequencies obtained
  • Computationally analytically;
  • Experimentally;
  • Through the FEM model for mesh variants (a) and (b) as shown in Figure 3.
Analytical calculation of the resonance and antiresonance frequencies is performed according to the relations in [1]. The calculated values are
f r t = 994 , 270 Hz , f a t = 1 , 121 , 000 Hz .
By substituting these values into the relation (1) we obtain the value of the effective coefficient of the electromechanical coupling
k e f f t = 0.4619 .
The experimentally determined values are taken from [1]
  • Resonance frequency f r e = 1 , 001 , 600 Hz ;
  • Antiresonance frequency f a e = 1 , 112 , 500 Hz .
By substituting these values into the relation (1) we obtain the effective coefficient of the electromechanical coupling
k e f f e = 0.4352 .
The values of the resonance and antiresonance frequencies for the finite element mesh variants (a) and (b) according to Figure 3 were determined in Section 3.1. By substituting into the (1) relation, we obtain the value of the effective coefficient of the electromechanical coupling
k e f f a = 0.4631
for the mesh variant (a) and
k e f f b = 0.4630
for the mesh variant (b).
The obtained result is consistent with the fact that the values of resonance and antiresonance frequencies are both shifted by 0.6% compared to the analytical solution. The calculation of k e f f given by the relation (1) can then be expected to be almost identical in value for both the analytical and numerical solutions.

3.3.2. Calculation of k e f f for Resonators with Wrap-Around Electrodes

The values of the resonance and antiresonance frequencies for both variants of the mesh according to Figure 3 were calculated in Section 3.2. By substituting into the relation (1) we obtain the value of the effective coefficient of the electromechanical coupling
k e f f a = 0.4490
for the wrap-around electrode according to the variant (a) and
k e f f b = 0.4503
for the wrap-around electrode according to the variant (b). The electromechanical coupling coefficient calculated for the fully electroded resonator in Section 3.3.1 is equal to 0.46 for both electrode structures. It can be thus concluded that the calculated value of the electromechanical coupling coefficient is lower in the case of resonators with wrap-around electrodes.

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
  • A total of 88% of resonator volume in configuration (a);
  • A total of 90% of resonator volume in configuration (b).
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
k t 2 = π 2 f r f a cotg π 2 f r f a .
As mentioned in the introduction to this chapter, the relationship between k e f f 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 k t , the following is true [2],
k t 2 = π 2 1 k e f f cotg π 2 1 k e f f .
For comparison, let us calculate the coefficient k t 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 k t are higher than the corresponding values of k e f f in all of the considered cases, but the trend is the same for both quantities. In the case of resonators with wrap-around electrodes, k t 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.

Funding

This research was funded by the Technical University of Liberec.

Data Availability Statement

The data presented in this study are available on request from the author. The data are not publicly available.

Conflicts of Interest

The author declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
FEMFinite element method

References

  1. Pustka, M. Spektrální Vlastnosti Kruhových Piezokeramických Rezonátorů, 1st ed.; VÚTS: Liberec, Czech Republic, 2019. [Google Scholar]
  2. Erhart, J.; Půlpán, P.; Pustka, M. Piezoelectric Ceramic Resonators, 1st ed.; Springer: Cham, Germany, 2017. [Google Scholar]
  3. Novák, J. FEM Analysis of Piezoelectric Resonator Polarization Process. Appl. Sci. 2021, 11, 12045. [Google Scholar] [CrossRef]
  4. Schnabel, P. Dispersion of thickness vibrations of piezoceramic disc resonators. IEEE Trans. Sonics Ultrason. 1978, 25, 16–23. [Google Scholar] [CrossRef]
  5. Sekimoto, H. Analysis of trapped energy resonators with circular electrodes. IEEE Trans. Sonics Ultrason. 1984, 31, 664–669. [Google Scholar] [CrossRef]
  6. Karlash, V.L. Effect of split or partial electrodes on the forced vibrations of bar-type piezoceramic transducers. Int. Appl. Mech. 2016, 52, 535–546. [Google Scholar] [CrossRef]
  7. Zhao, Z.; Wang, B.; Qian, Z.; Yang, J. Trapped-Energy Thickness-Extensional Mode of a Partially Electroded ZnO Thin-Film Resonator. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2018, 65, 1669–1679. [Google Scholar] [CrossRef]
  8. Erhart, J.; Tutu, S. Effective Electromechanical Coupling for the Partially Electroded Ceramic Resonators of Different Geometries. Ann. Dunarea Jos Univ. Galati Fascicle IX Metall. Mater. Sci. 2015, 2, 7–16. [Google Scholar]
  9. Bezverkhy, O.; Zinchuk, L.; Karlash, V. Modeling and experimental verification of the thin multi-electrode piezoceramic bars’ forced vibrations. Math. Model. Comput. 2016, 3, 1–11. [Google Scholar] [CrossRef] [Green Version]
  10. Stefan, O. Obrysové kmity kruhových keramických resonátorů. Československý Časopis pro Fyziku 1970, 20, 113–122. [Google Scholar]
  11. Onoe, M. Contour Vibrations of Isotropic Circular Plates. J. Acoust. Soc. Am. 1956, 28, 1158–1162. [Google Scholar] [CrossRef]
  12. Wu, L.; Chure, M.C.; Chen, Y.C.; Wu, K.K.; Che, B.H. Electrode Size and Dimensional Ratio Effect on the Resonant Characteristics of Piezoelectric Ceramic Disk. In Ceramic Materials, Progress in Modern Ceramics; Shi, F., Ed.; IntechOpen: London, UK, 2012. [Google Scholar]
  13. Kunkel, H.A.; Locke, S.; Pikeroen, B. Finite-Element Analysis of Vibrational Modes in Piezoelectric Ceramic Disks. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 1990, 37, 316–328. [Google Scholar] [CrossRef] [PubMed]
  14. Ivina, N.F. Analysis of the natural vibrations of circular piezoceramic plates with partial electrodes. Acoust. Phys. 2001, 47, 714–720. [Google Scholar] [CrossRef]
  15. Huang, C.H.; Lin, Y.C.; Ma, C.C. Theoretical Analysis and Experimental Measurement for Resonant Vibration of Piezoceramic Circular Plates. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2004, 51, 12–24. [Google Scholar] [CrossRef] [PubMed]
  16. Guo, N.; Cawley, P.; Hitchins, D. The Finite Element Analysis of the Vibration Characteristics of Piezoelectric Discs. J. Sound Vib. 1992, 159, 115–138. [Google Scholar] [CrossRef]
  17. Toftekaer, J.F.; Benjeddou, A.; Hogsberg, J. General numerical implementation of a new piezoelectric shunt tuning method based on the effective electromechanical coupling coefficient. Mech. Adv. Mater. Struct. 2019, 27, 1908–1922. [Google Scholar] [CrossRef]
  18. Rálek, P. Spectral Analysis of Discretized Model of Piezoelectric Resonator. Ph.D. Thesis, Technical University of Liberec, Liberec, Czech Republic, 2006. [Google Scholar]
  19. Trindade, M.A.; Benjeddou, A. Effective Electromechanical Coupling Coefficients of Piezoelectric Adaptive Structures: Critical Evaluation and Optimization. Mech. Adv. Mater. Struct. 2009, 16, 210–223. [Google Scholar] [CrossRef]
  20. Noliac Piezoceramic Material NCE51. Available online: http://www.noliac.com/products/materials/nce51/ (accessed on 1 November 2021).
  21. Rogacheva, N.N. The Dependence of the Electromechanical Coupling Coefficient of Piezoelectric Elements on the Position and Size of the Electrodes. J. Appl. Maths. Mech. 2001, 65, 317–326. [Google Scholar] [CrossRef]
  22. Nováková, K.; Pustka, M.; Novák, J. Frekvenční Spektrum Tloušťkově Kmitajících Rezonátorů s přEtaženou Elektrodou; Technical Report No. KMS-3062-01, Piezoceram, s. r. o., Libřice; Technical University of Liberec: Liberec, Czech Republic, 2006. [Google Scholar]
Figure 1. Piezoelectric transducers with wrap-around electrodes. Wrapping of the electrode from the bottom to the top surface is apparent, (a) wrapping in the shape of a circular segment; (b) wrapping in the shape of an intersection of circles.
Figure 1. Piezoelectric transducers with wrap-around electrodes. Wrapping of the electrode from the bottom to the top surface is apparent, (a) wrapping in the shape of a circular segment; (b) wrapping in the shape of an intersection of circles.
Applsci 12 04634 g001
Figure 2. Diagram of the dimensions of the analyzed transducers, (a) in the shape of circular segment, (b) in the shape of an intersection of circles.
Figure 2. Diagram of the dimensions of the analyzed transducers, (a) in the shape of circular segment, (b) in the shape of an intersection of circles.
Applsci 12 04634 g002
Figure 3. Finite element mesh corresponding to the symmetric halves of the resonators, (a) in the shape of circular segment, (b) in the shape of an intersection of circles.
Figure 3. Finite element mesh corresponding to the symmetric halves of the resonators, (a) in the shape of circular segment, (b) in the shape of an intersection of circles.
Applsci 12 04634 g003
Figure 4. Oscillation shape and frequency for mesh variant (a), for mesh variant (b).
Figure 4. Oscillation shape and frequency for mesh variant (a), for mesh variant (b).
Applsci 12 04634 g004
Figure 5. Shape of the oscillation in the electrode symmetry plane, for the electrode structure (a).
Figure 5. Shape of the oscillation in the electrode symmetry plane, for the electrode structure (a).
Applsci 12 04634 g005
Table 1. Dimensions of the transducers.
Table 1. Dimensions of the transducers.
d (mm)47.9t (mm)2.02
s 1 (mm)4.66 r 1 (mm)9.50
s 2 (mm)2.00 r 2 (mm)7.50
Table 2. Material properties of NCE51 ceramics.
Table 2. Material properties of NCE51 ceramics.
ε 11 S ( 1 ) 906 c 12 E ( Pa ) 2.04 × 10 10
ε 33 S ( 1 ) 823 c 13 E ( Pa ) 8.85 × 10 10
e 31 ( C / m 2 ) −4.80 c 33 E ( Pa ) 17.2 × 10 10
e 33 ( C / m 2 ) 17.2 c 44 E ( Pa ) 2.04 × 10 10
e 15 ( C / m 2 ) 13.7 c 66 E ( Pa ) 2.30 × 10 10
c 11 E ( Pa ) 12.9 × 1010 ϱ ( kg / m 3 ) 7850
Table 3. Calculated values k t .
Table 3. Calculated values k t .
Fully ElectrodedWrap-Around (a)Wrap-Around (b)
k t 0.50130.48670.4858
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Novák, J. FEM Analysis of the Effect of Polarization on the Electromechanical Coupling Factor of Resonators with a Wrap-Around Electrode. Appl. Sci. 2022, 12, 4634. https://doi.org/10.3390/app12094634

AMA Style

Novák J. FEM Analysis of the Effect of Polarization on the Electromechanical Coupling Factor of Resonators with a Wrap-Around Electrode. Applied Sciences. 2022; 12(9):4634. https://doi.org/10.3390/app12094634

Chicago/Turabian Style

Novák, Josef. 2022. "FEM Analysis of the Effect of Polarization on the Electromechanical Coupling Factor of Resonators with a Wrap-Around Electrode" Applied Sciences 12, no. 9: 4634. https://doi.org/10.3390/app12094634

APA Style

Novák, J. (2022). FEM Analysis of the Effect of Polarization on the Electromechanical Coupling Factor of Resonators with a Wrap-Around Electrode. Applied Sciences, 12(9), 4634. https://doi.org/10.3390/app12094634

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop