Next Article in Journal
Progressive Line Processing of Kernel RX Anomaly Detection Algorithm for Hyperspectral Imagery
Next Article in Special Issue
Design and Validation of a 150 MHz HFFQCM Sensor for Bio-Sensing Applications
Previous Article in Journal
EPPRD: An Efficient Privacy-Preserving Power Requirement and Distribution Aggregation Scheme for a Smart Grid
Previous Article in Special Issue
Pinhole Zone Plate Lens for Ultrasound Focusing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Semi-Analytical Solution for the Thickness-Vibration of Centrally Partially-Electroded Circular AT-Cut Quartz Resonators

State Key Laboratory of Mechanics and Control of Mechanical Structures/College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(8), 1820; https://doi.org/10.3390/s17081820
Submission received: 28 June 2017 / Revised: 26 July 2017 / Accepted: 2 August 2017 / Published: 7 August 2017
(This article belongs to the Special Issue Acoustic Wave Resonator-Based Sensors)

Abstract

:
Vibration frequencies and modes for the thickness-shear vibrations of infinite partially-electroded circular AT-cut quartz plates are obtained by solving the two-dimensional (2D) scalar differential equation derived by Tiersten and Smythe. The Mathieu and modified Mathieu equations are derived from the governing equation using the coordinate transform and the collocation method is used to deal with the boundary conditions. Solutions of the resonant frequencies and trapped modes are validated by those results obtained from COMSOL software. The current study provides a theoretical way for figuring out the vibration analysis of circular quartz resonators.

1. Introduction

Nowadays, acoustic wave resonators are widely used for frequency generation and control in telecommunications, as well as mass and acceleration sensors. Quartz is the most widely used crystal for resonators, due to its piezoelectric effects. The desire for better resonators has increased with the rapid development of electric devices, which calls for thorough studies on the vibration of the quartz plate. In most applications, quartz crystal plate works with thickness-shear (TS) vibrations [1]. However, because of the complicated properties of quartz crystal, such as anisotropy and electro-mechanical coupling, it is difficult to obtain an exact three-dimensional analytical solution. In order to solve the problem, Mindlin and his co-authors firstly developed a two-dimensional plate theory for solving the vibrations of elastic and piezoelectric plates [2,3,4,5], which is especially suitable for analyzing low-order modes of thickness-shear vibrations. Many researchers have adopted Mindlin’s theory to solve vibration of quartz plates in various situations [6,7,8,9,10,11]. Lee [12] established another effective two-dimensional plate theory using trigonometric expansion to express the electric potential and displacement, by which he tried to analyze the stretch and forced vibrations of crystal plates [13,14]. Furthermore, Tiersten and Smythe suggested a two-dimensional scalar differential equation for high-order modes in AT-cut [15] and Stress Compensated-cut (SC-cut) [16] quartz plates, which proved to be quite accurate yet simple. In many of the two-dimensional studies, the plates were assumed to be infinitely large in the x 1 and x 3 directions, and the effects induced by boundaries were considered by in-plane wavenumbers in the x 1 and x 3 directions, which are small compared with the wavenumber in the x 2 direction. This simplification has been shown to be effective in analyzing actual plate vibrations [17,18,19]. In this paper, we will adopt this two-dimensional scalar differential equation for the analysis of various TS modes.
For most quartz resonators, the shape of the crystal plates or the electrodes on the top and bottom of the plates is either rectangular or circular. It is true that rectangular geometry is widely used for general miniaturized quartz resonators. For some specific applications, such as base stations, which do not have strong limitation on space, a circular design of the quartz plate and/or the electrodes on it can show better performance for energy trapping resonance characteristics than a rectangular plate. It is well known that electrode corners in rectangular design can cause electric field concentration, which is associated with degrading effects or even failure of the materials involved. That problem disappears or is very much reduced when a circular design is used, since in this way the electric field concentration can be avoided with a minimal difference in manufacturing. In the case of rectangular plates, the two-dimensional scalar differential equations can be easily decoupled and solved in the Cartesian coordinates, and most of theoretical and numerical studies on TS vibrations mentioned above focused on rectangular plates. In the case of circular resonators, however, the variables in the two-dimensional scalar differential equations cannot be separated easily due to material anisotropy. Therefore, it is necessary to develop a semi-analytical method for circular quartz resonators.
Due to the difficulty of theoretical study on circular resonators, some researchers have turned to numerical solutions. Yong et al. [20] firstly applied Finite Element Method (FEM) to solve the thickness-shear motions of circular crystal plates, followed by Wang et al. [21]. Furthermore, Liu et al. [22,23] developed a new differential quadrature hierarchical FEM method, which can improve accuracy and calculation efficiency. Generally speaking, the FEM approach requires a large amount of time for computation, especially the post process of extracting resonance frequencies and modes is extremely complicated. For theoretical deduction of TS vibration modes, Tiersten proposed a perturbation method, where an incremental item due to anisotropy is added to a nearby isotropic solution [15]. Wang et al. [24] applied the method to solve TS vibration problems of circular plates in order to obtain a truncated 2D equation, while Yang et al. [25] discussed the vibration modes of a circular crystal plate with transversely varying thickness. More recently, He et al. solved the two-dimensional scalar differential equation for a fully-unelectroded finite circular AT-cut quartz plate by coordinate transform and variable separation in elliptical coordinates, by which the original equations were rewritten into the Mathieu and modified Mathieu equations [26]; their article showed better results than the previous perturbation method.
For the actual application of resonators, the electrode usually covers part of the surface on either the top or the bottom, so a partially-electroded model is more accurate in describing an actual resonator than the fully-electroded model. In this paper, the TS vibrations of infinite circular AT-cut quartz plates partially covered with circular electrodes are examined in detail. In Section 2, Tiersten’s two-dimensional scalar differential equations are introduced and transformed into elliptical coordinates for decoupling. By introducing boundary and constitutive conditions, the infinite plate's TS vibration frequencies and mode shapes are calculated theoretically, with the numerical models illustrated and discussed in Section 3. In Section 4, FEM simulation by the commercial software COMSOL is described. The FEM solutions are compared with the theoretical results and show good agreement. A thorough understanding of TS vibrations for partially-electroded circular quartz plates can be a theoretical guide for the design, calibration and manufacturing optimization of circular quartz resonators.

2. Theoretical Deduction

2.1. Governing Equations

An infinite AT-cut circular quartz plate with two circular identical electrodes placed concentrically on both the top and bottom was studied, as shown in Figure 1. The plate has a diameter of 2 R 1 and thickness of 2 h ; each electrode has a diameter of 2 R 0 and thickness of h e . We set the plate radius R 1 = + for simplicity. The vibrations of the electroded region and the unelectroded region are governed by different equations.
The thickness-shear vibration consists of different modes. The displacement u 1 n ( x 1 , x 3 , t ) for the nth-order thickness-shear vibration of AT-cut quartz plate is defined by [5]:
u 1 ( x 1 , x 2 , x 3 , t ) = n = 1 , 3 , 5 u 1 n ( x 1 , x 3 , t ) sin ( n π x 2 2 h ) ,
where n = 1 represents the fundamental thickness-shear mode, and larger values of n represent higher-order overtone modes. According to [15], the items u 1 n representing vibrations of electroded and unelectroded regions are governed by different two-dimensional scalar differential equations.
The scalar equation for the electroded region is shown as:
M n 2 u 1 n x 1 2 + c 55 2 u 1 n x 3 2 + ρ ( ω 2 ω ¯ 2 ) u 1 n = 0 ,
While the scalar equation for the unelectroded region has the similar form as:
M n 2 u 1 n x 1 2 + c 55 2 u 1 n x 3 2 + ρ ( ω 2 ω 2 ) u 1 n = 0 ,
where
M n = c 11 + ( c 12 + c 66 ) r + 4 ( r c ¯ 66 c 66 ) ( r c 22 + c 12 ) c 22 n π κ cot κ n π 2 ,
c ¯ 66 = c 66 + e 26 2 ε 22 ,   κ = ( c ¯ 66 c 22 ) 1 2 ,   r =   c 12 + c 22 c ¯ 66 c 22
c ^ 66 = c ¯ 66 ( 1 8 k ¯ 26 2 n 2 π 2 2 R ) ,   k ¯ 26 2 = e 26 2 c ¯ 66 ε 22 ,   R = 2 ρ h ρ h ,
ω ¯ = n 2 π 2 4 c ^ 66 ρ h 2 ,   ω = n 2 π 2 4 c ¯ 66 ρ h 2 ,
in which ρ is the mass density of quartz; 2hand ρ are the thickness and the density of the electrodes, respectively; cpq, eip and εij are the compact forms of elastic, piezoelectric and dielectric constants; and ω and ω ¯ are the fundamental TS frequencies for infinite unelectroded and electroded plates respectively, as seen in [18].
To deal with the governing equations we used a method similar to that adopted by He et al. [26]. For both the electroded and unelectroded regions, we introduced the following coordinate transformation:
{ x 1 = λ y 1 x 3 = μ y 3 ,
where λ and μ are defined as:
λ = 1 ω M n ρ ,    μ = 1 ω c 55 ρ ,
and Equations (2) and (3) separately become:
2 u 1 n y 1 2 + 2 u 1 n y 3 2 + ( ω 2 ω 2 ω ¯ 2 ω 2 ) u 1 n = 0 ,
2 u 1 n y 1 2 + 2 u 1 n y 3 2 + ( ω 2 ω 2 ω 2 ω 2 ) u 1 n = 0 ,
The boundary of the electrode can be defined by:
x 1 2 + x 3 2 = R 0 2 ,
After the coordinate transformation the circular boundary becomes an ellipse shown as:
λ 2 y 1 2 + μ 2 y 3 2 = R 1 2 ,
and we introduce the elliptical coordinate in order to deal with the boundary conditions:
{ y 1 = c sinh ξ sin η y 3 = c cosh ξ cos η ,
where c is half the focal distance of the ellipse which represents the electrode boundary in elliptical coordinates; and ξ is the radial variable; while η is the angular variable. Then Equations (10) and (11) are transformed into Equations (15) and (16) respectively:
2 u 1 n ξ 2 + 2 u 1 n η 2 + 2 q 1 ( cosh 2 ξ cos 2 η ) u 1 n = 0 ,
2 u 1 n ξ 2 + 2 u 1 n η 2 + 2 q 2 ( cosh 2 ξ cos 2 η ) u 1 n = 0 ,
in which:
{ q 1 = c 2 4 ( ω 2 ω 2 ω ¯ 2 ω 2 ) q 2 = c 2 4 ( ω 2 ω 2 ω 2 ω 2 ) ,
We now perform a variable separation to u 1 n :
u 1 n ( ξ , η ) = u ( ξ ) v ( η ) ,
Equations (15) and (16) can be rewritten as:
{ d 2 v η 2 + [ λ 2 q i cos ( 2 η ) ] v ( η ) = 0 d 2 u ξ 2 + [ 2 q i cosh ( 2 ξ ) λ ] u ( ξ ) = 0 ,
where i = 1 , 2 ; and λ is the separation constant which is related to q . The first and second equations of Equation (19) are known as the Mathieu equation and modified Mathieu equation, and the subscript i = 1 , 2 represents the equation used for the electroded and unelectroded parts, respectively.

2.2. Mathieu and Modified Mathieu Equations

The exact solutions of the Mathieu and modified Mathieu equations can be obtained from [27]. For the Mathieu equation, there are four kinds of periodic solutions with the period 2 π , representing angular distributions, as shown by:
{ c e 2 m ( η , q ) = k = 0 A 2 k 2 m ( q ) cos 2 k η c e 2 m + 1 ( η , q ) = k = 0 A 2 k + 1 2 m + 1 ( q ) cos ( 2 k + 1 ) η s e 2 m + 1 ( η , q ) = k = 0 B 2 k + 1 2 m + 1 ( q ) sin ( 2 k + 1 ) η s e 2 m + 2 ( η , q ) = k = 0 B 2 k + 2 2 m + 2 ( q ) sin ( 2 k + 2 ) η ,
Here Equation (20) is called the Mathieu function, in which c e 2 m and c e 2 m + 1 represent symmetric solutions, while s e 2 m + 2 and s e 2 m + 1 represent anti-symmetric solutions. The parameters m = 0 , 1 , 2 represent solution orders. A k m and B k m are coefficients determined by q .
For the modified Mathieu equation there are also four kinds of solutions, representing radial distributions, as shown by:
{ C e 2 m ( ξ , q ) = k = 0 A 2 k 2 m ( q ) cosh 2 k ξ C e 2 m + 1 ( ξ , q ) = k = 0 A 2 k + 1 2 m + 1 ( q ) cosh ( 2 k + 1 ) ξ S e 2 m + 1 ( ξ , q ) = k = 0 B 2 k + 1 2 m + 1 ( q ) sinh ( 2 k + 1 ) ξ S e 2 m + 2 ( ξ , q ) = k = 0 B 2 k + 2 2 m + 2 ( q ) sinh ( 2 k + 2 ) ξ ,
where m = 0 , 1 , 2 also represents the order of the solutions. Equation (21) is called the modified Mathieu function in which C e 2 m and C e 2 m + 1 represent symmetric solutions while S e 2 m + 2 and S e 2 m + 1 represent anti-symmetric solutions. A k m and B k m are coefficients, the same as those in Equation (20). All the modified Mathieu functions can be expanded with Bessel functions of the first and second kinds shown as:
C e 2 m ( ξ , q ) = M c 2 m J ( ξ , q ) = 1 A 0 2 m k = 0 ( 1 ) k + m A 2 k 2 m ( q ) J k ( v 1 ) Z k J ( v 2 ) , C e 2 m + 1 ( ξ , q ) = M c 2 m + 1 J ( ξ , q ) = 1 A 1 2 m + 1 k = 0 ( 1 ) k + m A 2 k + 1 2 m + 1 ( q ) × [ J k ( v 1 ) Z k + 1 J ( v 2 ) + Z k J ( v 2 ) J k + 1 ( v 1 ) ] , S e 2 m + 1 ( ξ , q ) = M s 2 m + 1 J ( ξ , q ) = 1 B 1 2 m + 1 k = 0 ( 1 ) k + m B 2 k + 1 2 m + 1 ( q ) × [ J k ( v 1 ) Z k + 1 J ( v 2 ) Z k J ( v 2 ) J k + 1 ( v 1 ) ] , S e 2 m + 2 ( ξ , q ) = M s 2 m + 2 J ( ξ , q ) = 1 B 2 2 m + 2 k = 0 ( 1 ) k + m B 2 k + 2 2 m + 2 ( q ) × [ J k ( v 1 ) Z k + 1 J ( v 2 ) Z k J ( v 2 ) J k + 2 ( v 1 ) ] ,
where v 1 = q exp ( ξ ) ; v 2 = q exp ( ξ ) ; and J = 1 , 2 . A m k and B m k are mentioned in Equation (21). When J = 1 , Z k J ( v ) = J k ( v ) , which is a Bessel function of the first kind, and Equation (22) becomes the first kind of modified Mathieu function. When J = 2 , Z k J ( v ) = Y k ( v ) , which is a Bessel function of the second kind, and Equation (22) becomes the second kind of modified Mathieu function. Similar to the definition of the Hankel function, there are also the third and fourth kinds of modified Mathieu functions defined as:
{ M c m ( 3 ) , ( 4 ) ( ξ , q ) = M c m ( 1 ) ( ξ , q ) ± i M c m ( 2 ) ( ξ , q ) M s m ( 3 ) , ( 4 ) ( ξ , q ) = M s m ( 1 ) ( ξ , q ) ± i M s m ( 2 ) ( ξ , q ) ,
Different modified Mathieu functions are used in different kinds of solutions representing different vibration modes [27], which we will discuss in the following sections.

2.3. Vibration Modes and Displacement Solutions

Because of the features of the Mathieu and modified Mathieu functions, the vibration modes can be divided into symmetric and anti-symmetric modes.
For symmetric modes, the displacement u 1 n in the electroded region can be written as
u 1 n ( η , ξ ) = m = 0 A m c e m ( η , q 1 ) M c m ( 1 ) ( ξ , q 1 ) ,
In the unelectroded region, the displacement u 1 n can be written as
u 1 n ( η , ξ ) = m = 0 B m c e m ( η , q 2 ) M c m ( 3 ) ( ξ , q 2 ) ,
For the anti-symmetric modes, the displacement for the electroded region can be written as:
u 1 n ( η , ξ ) = m = 0 A m s e m ( η , q 1 ) M s m ( 1 ) ( ξ , q 1 ) ,
In the unelectroded region the displacement is:
u 1 n ( η , ξ ) = m = 0 B m s e m ( η , q 2 ) M s m ( 3 ) ( ξ , q 2 ) ,
where A m , B m , A m and B m are undetermined constants. In Equation (24) and (26) we chose the first kind of modified Mathieu function, which represents a standing wave because the electroded region is bounded and there could be a displacement in the middle of the region. In Equations (25) and (27), the third kind of modified Mathieu function, which represents a divergent wave, was chosen because the unelectroded region extends to infinity and the waves can propagate outwards.

2.4. Boundary Conditions

We take the symmetric modes as examples. The far-field radiation condition of the unelectroded region has been expressed in [4] by:
u 1 n = 0 ,    ξ = ξ ,
On the interface between the electroded and the unelectroded regions in x 1 x 3 plane, the displacement and stress continuity conditions can be expressed as:
[ u 1 n ξ ] = 0 [ u 1 n ] = 0 ,    ξ = ξ 0 ,
Equation (28) is taken into account first. The substitution of Equation (25) into Equation (28) results in the following equation:
m = 0 B m c e m ( η , q 2 ) M c m ( 3 ) ( ξ , q 2 ) = 0 ,
According to the characteristics of the third kind of modified Mathieu function, M c m ( 3 ) approaches zero when ξ = ξ and Equation (30) is satisfied automatically.
With respect to the continuity conditions, we can substitute Equations (24) and (25) into Equation (29) to get:
m = 0 A m c e m ( η , q 1 ) M c m ( 1 ) ( ξ 0 , q 1 ) m = 0 B m c e m ( η , q 2 ) M c m ( 3 ) ( ξ 0 , q 2 ) = 0 , m = 0 A m c e m ( η , q 1 ) M c m ( 1 ) ( ξ 0 , q 1 ) m = 0 B m c e m ( η , q 2 ) M c m ( 3 ) ( ξ 0 , q 2 ) = 0 ,

2.5. The Collocation Method

In Equations (28) and (29), the displacements are written as a sum of infinite series of terms and we took several of them as the approximate displacements, and the result converged to the accurate value as more items are considered.
Here we adopted the collocation method, i.e., the continuous equations were satisfied at a series of collocation points on the interface contour between the electroded and the unelectroded parts of the structure in the x 1 x 3 plane.
We now considered the symmetric modes of the infinite plate. Due to the features of the Mathieu functions and the symmetry of the model, only half the boundary in the first and fourth quadrants needed to be examined for the displacement and stress continuation. After choosing collocation points we found linear homogeneous equations set from Equation (31) for the undetermined constants A m and B m . To ensure the existence of the nontrivial solutions A m and B m , the determinant of the coefficient matrix should be zero, which can determine the value of ω.
In order to use the collocation method, the convergence of the vibration frequency should be first examined to ensure validity.
For example, an infinite plate with an electrode attached on each of the two major surfaces (in the x 1 x 3 plane) has a thickness of 2h = 0.1 mm, the electrode has a radius of R0 = 0.8 mm, and the mass ratio between the electrode and the plate is taken to be R = 2 ρ h / ρ h = 0.05 . The material constants of AT-cut quartz can be obtained from [28]. When solving the Mathieu equation we took the first thirty terms of each series to conduct the calculation. Two to nine points on half the boundary were chosen to calculate the frequency ω of the fundamental mode when n = 1 , and the result is shown in Figure 2:
It can be seen from Figure 2 that when the number of the chosen points increased to five, the frequency calculated became steady and convergence was ensured.
We then checked the displacement and stress continuity of the interface between the electroded and the unelectroded parts on the elliptical coordinate ξ = ξ 0 . We arbitrarily chose boundary points between the collocation points for examination. For Example 1, with eight collocation points chosen on half of the interface, and the first eight terms of the series for the calculation of the displacement and stress remaining, we obtained the frequency of the fundamental TS mode and the expression of the displacement and the stress. The displacement and stress of two arbitrarily-chosen points on the interface between the collocation points (e.g., (ξ0,π/4) and (ξ0,5π/6)) are shown in Table 1.
Comparing the data in Table 1, we found that the displacement was continuous at the arbitrarily-chosen points besides the collocated points on the boundary, and the corresponding stress difference across the boundary was also acceptable. Thus, we concluded that the collocation method can be used to deal with the boundary conditions under question.

3. Numerical Results and Discussion

In this section, we will discuss the displacement distribution of the TS vibrations for Example 1.
The displacement field for the fundamental TS mode, u 1 1 (n = 1), is shown in Figure 3a, while the displacement field for the 3rd-order overtone TS mode, u 1 3 ( n = 3 ), is shown in Figure 3b. Comparing Figure 3a,b, we found that the shapes of the distributions were almost the same, but that the displacement u 1 n . was more concentrated in the electroded region when n = 3 than when n = 1 . Due to the material anisotropy of the quartz crystal plate, the concentration area of the vibration in the middle of the plate was an ellipse instead of a circle, where the long axis lies in the x 1 direction and the short axis in the x 3 direction. Meanwhile, the frequency for the third thickness-shear overtone was about triple that of the fundamental mode.
Four other thickness-shear modes with n = 1 are shown in Figure 4. Among these, Figure 4a,b are modes symmetric to axis- x 3 , while Figure 4c,d are anti-symmetric modes. Figure 4a–d are not the fundamental modes for n = 1 and they have more nodal lines under the electrodes, which will cause charge cancellation and are unwanted spurious modes. In fact, the mode in Figure 3a is the most useful one.
Furthermore, by observing the vibration in the electroded and unelectroded areas, it was obvious that the vibration amplitude was low in the unelectroded area and most of the vibration energy was confined to the electroded area. This phenomenon is called energy trapping and it is used for device packaging.

4. Verification in the COMSOL Package

4.1. Model, Equation and Boundary Conditions

We then used the COMSOL Multiphysics software to solve the problem in order to verify the numerical results solved by the Mathieu functions in Section 3.
Since the plate is unbounded, we needed to use the infinite element domain whose geometry is set cylindrical. We created an infinite two-dimensional model for the crystal plate in COMSOL. The electroded zone was within the circle at radius R 0 = 0.8 mm, and the outer zone was the unelectroded zone. The outer side of the unelectroded zone tended to infinity.
In addition, the two circles outside the electrode were introduced for mesh, the radius of which needed to be big enough to make the model accurate, yet small enough to be economic in computation. After examination, we determined the appropriate size for the outer radiuses of the two rings to be 1.8 mm and 2.5 mm. We used extra fine free triangles in the model to ensure that the results were accurate enough; the meshed model in the x 1 x 3 plane is shown in Figure 5.
The problem in this article was mode analysis by two-dimensional equations, which is an eigenvalue problem that can be solved by a user-defined Partial Differential Equation (PDE) interface. Coefficient form PDE, which provides a general interface for specifying and solving many well-known PDEs in the coefficient form, was chosen from the PDE interfaces in module mathematics shown as:
λ 2 e a 2 u t 2 λ d a u t + ( c u α u + γ ) + β u + a u = f ,
where u is the unknown scalar to be solved; t and f represent the time and source term, respectively; λ is the eigenvalue; and e a ,   d a ,   c ,   α ,   γ ,   β   and   a are coefficients defined by users in which c is a matrix because of the material anisotropy. This is for a 2D domain, = [ / x , / y ] .
To match with the coefficient form PDE, Equations (2) and (3) were rewritten in the following form:
f 2 ( 4 π 2 ρ ) u 1 n + M n 2 u 1 n x 1 2 + c 55 2 u 1 n x 3 2 ω ¯ 2 u 1 n = 0 ,
f 2 ( 4 π 2 ρ ) u 1 n + M n 2 u 1 n x 1 2 + c 55 2 u 1 n x 3 2 ω 2 u 1 n = 0 ,
where f is the characteristic frequency. Comparing Equations (33) and (34) with Equation (32), we found the user-defined parameters listed as Equations (35) and (36), respectively:
e a = 4 π 2 ρ , c = [ M n 0 0 c 55 ] , a = ρ ω ¯ 2 , d a = β = γ = 0 ,
e a = 4 π 2 ρ , c = [ M n 0 0 c 55 ] , a = ρ ω 2 , d a = β = γ = 0 ,
The continuous conditions between the electroded and unelectroded regions were satisfied automatically. At the outermost circle boundaries, the displacements were set as zero to simulate the radiation condition of Equation (28) and suppress other vibrating modes.
We set the number of eigenvalues that we were interested in and the search ranges according to the theoretical values obtained above. After that, eigenvalues and mode displacements were able to be obtained automatically.

4.2. Simulation Results

COMSOL was used to analyze the model to verify the numerical results and we obtained the same TS vibration modes as those listed in Figure 2 and Figure 3. Both the vibration modes and frequencies will be compared in this part.
The results of the fundamental and third overtone TS modes obtained through FEM simulation are shown in Figure 6. Meanwhile, two spurious modes are shown in Figure 7.
The mode frequencies obtained from the theoretical analysis and FEM simulation are also listed in Table 2, where the first and last rows show the fundamental and third order overtone TS modes, respectively, and the four rows in the middle illustrate four typical spurious modes whose frequencies are close to that of fundamental mode ( n = 1 ). We can see that for each mode the error of frequency is less than 100 ppm, which is acceptable in the analysis and design of resonators. That validates the correctness of our theoretical analysis based on the coordinate transform and the Mathieu function series expansion.

5. Conclusions

The exact frequencies of thickness-shear and thickness-twist modes along with mode shapes for infinite, centrally partially-electroded circular AT-cut quartz plates were obtained through theoretical analysis, by the application of coordinate transform and using Mathieu and modified Mathieu functions. The results were compared with 2D FEM simulations by COMSOL software and good agreement was achieved, illustrating the correctness of our theoretical method. From the mode shapes obtained it was easy to find that all the vibrations were concentrated in the electroded region and there was little vibration in the outer, unelectroded area. This phenomenon is called energy trapping and is used for resonator packaging, which does not affect TS vibrations. Our results are of theoretical importance for the design of circular quartz resonators.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 11502108, 11232007, 51405225), the Program for New Century Excellent Talents in Universities (No. NCET-12-0625), the Natural Science Foundation of Jiangsu Province (No. BK20140037, BK20140808), the Fundamental Research Funds for Central Universities (No. NE2013101), and a project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Author Contributions

B.W. derived the formula and analyzed the numerical results; X.D. and X.Z. did all the numerical calculations and plotted all the figures and tables; Z.Q. proposed the study, corrected all the numerical results and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tiersten, H.F. Thickness vibrations of piezoelectric plates. J. Acoust. Soc. Am. 1963, 35, 53–58. [Google Scholar] [CrossRef]
  2. Mindlin, R.D. High frequency vibrations of crystal plates. Q. Appl. Math. 1961, 19, 51–61. [Google Scholar] [CrossRef]
  3. Mindlin, R.D. Forced Thickness-Shear and Flexural Vibrations of Piezoelectric Crystal Plates. J. Appl. Phys. 1952, 23, 83–88. [Google Scholar] [CrossRef]
  4. Tiersten, H.F.; Mindlin, R.D. Forced vibrations of piezoelectric crystal plates. Q. Appl. Math. 1962, 20, 107–119. [Google Scholar] [CrossRef]
  5. Mindlin, R.D. An Introduction to the Mathematical Theory of Vibrations of Elastic Plates; Yang, J., Ed.; World Scientific: Singapore, December 2006. [Google Scholar]
  6. Wang, J.; Zhao, W.H. The determination of the optimal length of crystal blanks in quartz crystal resonators. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2005, 52, 2023–2030. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, J.; Zhao, W.H.; Du, J.K. The determination of electrical parameters of quartz crystal resonators with the consideration of dissipation. Ultrasonics 2006, 44, e869–e873. [Google Scholar] [CrossRef] [PubMed]
  8. Zhang, C.L.; Chen, W.Q.; Yang, J.S. Electrically forced vibration of a rectangular piezoelectric plate of monoclinic crystals. Int. J. Appl. Electromagn. Mech. 2009, 31, 207–218. [Google Scholar]
  9. Wang, J.N.; Hu, Y.T.; Yang, J.S. Frequency spectra of at-cut quartz plates with electrodes of unequal thickness. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2010, 57, 1146–1151. [Google Scholar] [CrossRef] [PubMed]
  10. Chen, G.J.; Wu, R.X.; Wang, J.; Du, J.K.; Yang, J.S. Five-mode frequency spectra of x3-dependent modes in AT-cut quartz resonators. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2012, 59, 811–816. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, J.; Yang, J.S. Higher-order theories of piezoelectric plates and applications. Appl. Mech. Rev. 2000, 53, 87–99. [Google Scholar] [CrossRef]
  12. Lee, P.C.Y.; Nikodem, Z. An approximate theory for high-frequency vibrations of elastic plates. Int. J. Solids Struct. 1972, 8, 581–612. [Google Scholar] [CrossRef]
  13. Lee, P.C.Y.; Nakazawa, M.; Hou, J.P. Extensional vibrations of rectangular crystal plates. J. Appl. Phys. 1982, 53, 4081–4087. [Google Scholar] [CrossRef]
  14. Lee, P.C.Y.; Syngellakis, S.; Hou, J.P. A two-dimensional theory for high-frequency vibrations of piezoelectric crystal plates with or without electrodes. J. Appl. Phys. 1987, 61, 1249–1262. [Google Scholar] [CrossRef]
  15. Tiersten, H.F.; Smythe, R.C. Coupled thickness-shear and thickness-twist vibrations of unelectroded AT-cut quartz plates. J. Acoust. Soc. Am. 1985, 78, 1684–1689. [Google Scholar] [CrossRef]
  16. Stevens, D.S.; Tiersten, H.F. An analysis of doubly rotated quartz resonators utilizing essentially thickness modes with transverse variation. J. Acoust. Soc. Am. 1986, 79, 1811–1826. [Google Scholar] [CrossRef]
  17. Tiersten, H.F. Analysis of trapped energy resonators operating in overtones of coupled thickness-shear and thickness-twist. Symp. Freq. Control 1975, 76, 71–75. [Google Scholar]
  18. Tiersten, H.F.; Smythe, R.C. An analysis of contoured crystal resonators operating in overtones of coupled thickness shear and thickness twist. J. Acoust. Soc. Am. 1979, 65, 1455–1460. [Google Scholar] [CrossRef]
  19. He, H.J.; Liu, J.X.; Yang, J.S. Effects of mismatched electrodes on an AT-cut quartz resonator. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2012, 59, 281–286. [Google Scholar] [PubMed]
  20. Steward, J.T.; Yong, Y.K. Thickness-shear mode shapes and mass-frequency influence surface of a circular and electroded AT-cut quartz resonator. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1992, 39, 609–617. [Google Scholar]
  21. Wang, J.; Chen, L.P.; Du, J.K.; Hu, Y.T.; Li, G.Q. Finite element analysis of nonlinear thickness-shear vibrations of AT-cut quartz crystal plates. In Proceedings of the 2010 IEEE International Frequency Control Symposium, Newport Beach, CA, USA, 1–4 June 2010; pp. 392–396. [Google Scholar]
  22. Liu, B.; Xing, Y.F. Thickness-shear vibration analysis of rectangular quartz plates by a differential quadrature finite element method. AIP Conf. Proc. 2014, 41, 41–44. [Google Scholar]
  23. Liu, B.; Xing, Y.F.; Wang, W.; Yu, W.D. Thickness-shear vibration analysis of circular quartz crystal plates by a differential quadrature hierarchical finite element method. Compos. Struct. 2015, 131, 1073–1080. [Google Scholar] [CrossRef]
  24. Wang, J.; Chen, H.; Ma, T.F.; Du, J.K.; Yi, L.J. An analysis of thickness-shear vibrations of an annular plate with the Mindlin plate equations. In Proceedings of the 2015 Joint Conference of the IEEE International Frequency Control Symposium and European Frequency and Time Forum, Denver, CO, USA, 12–16 April 2015; pp. 402–405. [Google Scholar]
  25. Yang, J.S.; Tiersten, H.F. Effects of free edge on the vibration characteristics of a contoured, beveled quartz resonator. Int. J. Appl. Electromagn. Mech. 2014, 46, 649–661. [Google Scholar]
  26. He, H.J.; Yang, J.S.; Jiang, Q. Thickness-shear and thickness-twist vibrations of circular AT-cut quartz resonators. Acta Mech. Solida Sin. 2013, 26, 245–254. [Google Scholar] [CrossRef]
  27. Zhang, S.J.; Jin, J.M. Computation of Special Functions; John Wiley: New York, NY, USA, 1996. [Google Scholar]
  28. Tiersten, H.F. Linear Piezoelectric Plate Vibrations: Elements of the Linear Theory of Piezoelectricity and the Vibrations Piezoelectric Plates; Springer: New York, NY, USA, 1 January 1969. [Google Scholar]
Figure 1. A centrally partially-electroded circular AT-cut quartz plate.
Figure 1. A centrally partially-electroded circular AT-cut quartz plate.
Sensors 17 01820 g001
Figure 2. The number of the chosen points and the frequencies calculated.
Figure 2. The number of the chosen points and the frequencies calculated.
Sensors 17 01820 g002
Figure 3. (a) Fundamental thickness-shear mode for n = 1, R = 5%, f = 1.586632 × 107 Hz. (b) The third overtone thickness-shear mode n = 3, R = 5%, f = 4.731520 × 107 Hz.
Figure 3. (a) Fundamental thickness-shear mode for n = 1, R = 5%, f = 1.586632 × 107 Hz. (b) The third overtone thickness-shear mode n = 3, R = 5%, f = 4.731520 × 107 Hz.
Sensors 17 01820 g003
Figure 4. Thickness-shear modes for n = 1, R = 5%. (a) f = 1.606972 × 107 Hz; (b) f = 1.634720 × 107 Hz; (c) f = 1.614183 × 107 Hz; (d) f = 1.639936 × 107 Hz.
Figure 4. Thickness-shear modes for n = 1, R = 5%. (a) f = 1.606972 × 107 Hz; (b) f = 1.634720 × 107 Hz; (c) f = 1.614183 × 107 Hz; (d) f = 1.639936 × 107 Hz.
Sensors 17 01820 g004
Figure 5. The meshed infinite electrode model.
Figure 5. The meshed infinite electrode model.
Sensors 17 01820 g005
Figure 6. Fundamental and third overtone TS modes by FEM simulation in COMSOL (a) n = 1, f = 1.586634 × 10 7   Hz (b) n = 3, f = 4.731538 × 10 7   Hz .
Figure 6. Fundamental and third overtone TS modes by FEM simulation in COMSOL (a) n = 1, f = 1.586634 × 10 7   Hz (b) n = 3, f = 4.731538 × 10 7   Hz .
Sensors 17 01820 g006
Figure 7. Two spurious modes by FEM simulation in COMSOL (a) f = 1.606981 × 10 7   Hz (b) f = 1.614231 × 10 7   Hz .
Figure 7. Two spurious modes by FEM simulation in COMSOL (a) f = 1.606981 × 10 7   Hz (b) f = 1.614231 × 10 7   Hz .
Sensors 17 01820 g007
Table 1. The normalized displacement and stress of selected points on the boundary by solutions from electroded and unelectroded regions, respectively.
Table 1. The normalized displacement and stress of selected points on the boundary by solutions from electroded and unelectroded regions, respectively.
Results PointsDisplacementsStresses
ElectrodedUnelectrodedElectrodedUnelectroded
( ξ 0 , π / 4 ) 0.16450.1645−0.7959−0.7958
( ξ 0 , 5 π / 6 ) 0.15210.1521−0.7002−0.7002
Table 2. Comparison of TS mode frequencies obtained by theoretical analysis and FEM simulation.
Table 2. Comparison of TS mode frequencies obtained by theoretical analysis and FEM simulation.
TS ModesTheoretical Results (MHz)FEM Results by COMSOL (MHz)Error (ppm)
Fundamental n = 115.8663215.866351.84
Spurious mode 116.0697216.069815.6
Spurious mode 216.1418316.1423129
Spurious mode 316.3472016.3475521
Spurious mode 416.3993616.3998731
Overtone n = 347.3152047.315383.8

Share and Cite

MDPI and ACS Style

Wang, B.; Dai, X.; Zhao, X.; Qian, Z. A Semi-Analytical Solution for the Thickness-Vibration of Centrally Partially-Electroded Circular AT-Cut Quartz Resonators. Sensors 2017, 17, 1820. https://doi.org/10.3390/s17081820

AMA Style

Wang B, Dai X, Zhao X, Qian Z. A Semi-Analytical Solution for the Thickness-Vibration of Centrally Partially-Electroded Circular AT-Cut Quartz Resonators. Sensors. 2017; 17(8):1820. https://doi.org/10.3390/s17081820

Chicago/Turabian Style

Wang, Bin, Xiaoyun Dai, Xintao Zhao, and Zhenghua Qian. 2017. "A Semi-Analytical Solution for the Thickness-Vibration of Centrally Partially-Electroded Circular AT-Cut Quartz Resonators" Sensors 17, no. 8: 1820. https://doi.org/10.3390/s17081820

APA Style

Wang, B., Dai, X., Zhao, X., & Qian, Z. (2017). A Semi-Analytical Solution for the Thickness-Vibration of Centrally Partially-Electroded Circular AT-Cut Quartz Resonators. Sensors, 17(8), 1820. https://doi.org/10.3390/s17081820

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