Modeling of MEMS Resonator Piezoelectric Disc Partially Covered with Electrodes
Ismail Naciri^{1,}^{ *}, Lahoucine Elmaimouni^{1,}^{ *}, Jean-Etienne Lefebvre^{2}, Faniry Emilson Ratolojanahary^{3}, Mohamed Rguiti^{4}, Tadeusz Gryba^{2}
^{1}Laboratoire Sciences Ingénierie et Energie, Energie Renouvelable, Microsystèmes Acoustique et Micanique, Polydisciplinary Faculty of Ouarzazate, Ibn Zohr University, Morocco
^{2}The Institute of Electronics, Microelectronics and Nanotechnology, Opto-Acousto-Electronic Department, University of Valenciennes, France
^{3}Laboratory of applied Physics, Fianarantsoa University, Madagascar
^{4}Laboratoire des Matériaux Céramiques et procédés Associés, Université de Valenciennes, Maubeuge, France
Email address:
To cite this article:
Ismail Naciri, Lahoucine Elmaimouni, Jean-Etienne Lefebvre, Mohamed Rguiti, Faniry Emilson Ratolojanahary, Tadeusz Gryba. Modeling of MEMS Resonator Piezoelectric Disc Partially Covered with Electrodes. American Journal of Mechanics and Applications. Vol. 4, No. 1, 2016, pp. 1-9. doi: 10.11648/j.ajma.20160401.11
Received: August 31, 2016; Accepted: September 26, 2016; Published: October 19, 2016
Abstract: The Legendre polynomial method has been extended to the modeling of MEMS resonator disc partially covered with electrodes. The disc has been divided into two areas: one with electrodes and the other without electrodes. For each area, The Maxwell equations and the piezoelectric constitutive equations of motion are studied and solved to yield a frequency response and electrical behavior of the MEMS resonator applying a semi analytical method based on a Legendre polynomials series and trigonometric functions. However, the method allows incorporating the boundary conditions directly into the governing equations by assuming position-dependent of elastic constants, mass density and delta functions. The alternating electrical source is described by specific terms which are also introduced into the equation of motion. The formalism has been developed which allows for both harmonic and modal analyses. In order to validate our polynomial approach, numerical results are presented such as resonant and anti-resonant frequencies, electric input admittance, electromechanical coupling coefficient and field profiles of fully and partially metallized PZT5A resonator discs. The results obtained were compared with those obtained by an approximated analytical method. The developed software proves to be very efficient to retrieve the contour modes of all orders.
Keywords: MEMS Resonators, Legendre Polynomial Approach, Centralized Metallization, Piezoelectric Resonator Disc, Electrical Admittance, Resonant, Anti-resonant Frequencies
1. Introduction
When a piezoelectric material is subjected to a mechanical strain, electrical charges were generated and conversely. This phenomenon was widely exploited in various engineering applications such as Micro-Electro-Mechanical systems (MEMS) technology. MEMS technology has been obtained significant growth in its field of application such as in electro-optic modulators, ultrasonic detectors, accelerometers, transducers, oscillators, electromechanical sensors and actuators [1-6].
Although, the vibration characteristics of piezoelectric materials are extracted from the piezoelectric constitutive equations, linear piezoelectricity and the Maxwell equations [7-8].
Several methods allow calculating vibration characteristics of piezoelectric devices. In 1967, Eer Nisse [9] presented a vibrational method to analyze the vibrational behavior of piezoelectric disks. C.H.Huang and C.C Ma used Electronic speckle pattern interferometry method to study vibration characteristics for piezoelectric cylinders [10]. Kharouf et al [11] studied the vibrational axisymmetric characteristics of hollow cylinders and disks made from piezoelectric materials by using the Rayleigh-Ritz method. These methods give accurate results especially for simple geometries. In the case of complicated geometries, alternative methods are needed. The finite elements method was used to analyze vibrational modes of ceramic disk composed of PZT5A with the diameter to thicknesses ratio range from 0.2 to 10 [12]. By the same method, Guo et al [13] have been calculated the vibration characteristics of disks composed of PZT5A with diameter to thicknesses ratio ranging from 10 to 20.
Moreover, Legendre polynomial method is one of the most methods used recently to check the electrical response of piezoelectric resonators and gives excellent precision for various structures [14-21]. This method is based on expressing the components of electrical and mechanical displacement by double series of Legendre polynomials.
In this article, polynomial method was extended for studying the frequency spectrum of a piezoelectric PZT5A disc resonator divided into two areas. Centralized area is electroded whereas the annular area is free of electrodes.
2. Formulation of the Problem
As shown in figure 1, the studied geometrical configuration is a homogeneous cylinder of finite length, divided into two areas named N and M. M area contains electrodes from the top and the bottom. H and R are the thickness and the radius of the cylinder respectively.
Figure 1. Structure of a piezoelectric disc partially metallized with electric excitation source .
Studying the structure in cylindrical coordinates . The z-axis coincides with the axis of anisotropy and the origin is placed at the center of the disc.
The vibration characteristics of a piezoelectric materials are studied by following equations [22-23]:
(1)
where are the components of electrical displacements and denotes the stress tensor. are the components of mechanical displacement in the radial, circumferential and axial direction.
Under the assumption of small deformations, the strain mechanical displacements relations in terms of cylindrical coordinates system are expressed by [24].
(2)
Along this paper, we adopt the following change of variables:
Normalized mechanical displacement components: , , .
and .
Equations expressed in (1) can be written as:
(3a)
(3b)
(3c)
(3d)
Taking in account the axisymmetric character of electrical source and studied structure, the variation of wave field does not depend on direction. Moreover, the mechanical displacement in circumferential direction is uncoupled from the electric field. We can only tike in consideration the equations (3a), (3c) and (3d).
Two rectangular window functions and express the boundary conditions related to our structure.
(4)
The equations expressed in (4) are directly incorporated in the constitutive equations of our problem.
The normalized stress components are expressed by
(4a)
(4b)
(4c)
(4d)
The electric displacement components:
(5a)
(5b)
The index j=1, 2 refers to area M and area N respectively.
Both for the metalized and non-metalized area, we express the field of governing equations (4) and their corresponding boundary, symmetry and continuity conditions.
For the metallized area, we have: and with
In this area, the boundary, symmetry and continuity conditions are defined by following relations:
The field of equations (3) becomes:
(6a)
(6b)
(6c)
For the non-metalized area, we have: and with
We express also the boundary, symmetry and continuity conditions defined by following relations:
The field of equations (3) becomes:
(7a)
(7b)
(7c)
Mechanical and electrical displacement components for area M and N are defined respectively as follow:
(8)
(9)
where:
and
are the Legendre polynomials of degree m and n respectively. The expansions coefficients and are in Angstrom and Volts, respectively. The term in the equation of electric potential ensure that. The delta functions, and refers to the boundary and continuity conditions in the studied piezoelectric disc. Those functions are multiplied by different components of the electrical displacement and mechanical stress components.
Substituting the mechanical and electrical displacements (8) and (9) in the field of equations (6) and (7) respectively for area M and area N and multiplying all equations by . j and k runs from 0 to M and N respectively. Governing equations (6) and (7) are integrated over from 0 to for metalized area, from to 1 for non-metalized area while runs from -1 to 1. This formulation allows obtaining a system of equations with unknowns. Obtained system can be rewrite as following form:
(10)
(11)
where:
; ;
Denotes the normalized frequency expressed as follow: .
Based on foregoing formulation, we present simultaneously modal and harmonic analyses.
2.1. Harmonic Analyses
The harmonic analysis allows determining the normalized electric impedance which is expressed by using the expression of displacement current density in our piezoelectric structure.
From the average electrical current that flows through the electrode of surface S, we determine the normalized electric input admittance by following relation:
(12)
Where denote the static capacitance of the disc resonator.
From equation 12, we obtain:
(13)
The matrix element factors and can be expressed as:
2.2. Modal Analyses
Through modal analyses, resonant () and anti-resonant frequencies () are calculated by following relations, respectively.
(14)
(15)
where
and
3. Numerical Results
We start our calculations with properties in table (1), the computer program is established using Matlab software.
The radius and the thickness of the studies resonator disc are H=1mm, R=15mm. The truncation number for our calculation is limited on M=N=12.
3.1. Resonant and Antiresonant Frequencies
The obtained results using polynomial approach are compared with results calculated analytically [25]. Analytical method is presented in appendix A. Tables 2 and 3 show resonant and anti-resonant frequencies for the first five contour modes of a piezoelectric disc partially metallized with a metallization rate 80%. Associated accuracy is calculated by following relations:
Table 4 reveals the electromechanical coupling coefficient of the first five contour modes. This coefficient is calculated for a specific pairs of, by: .
As shown in those tables, the higher error is less than 0.8% which means good agreement was found.
3.2. Normalized Electric Admittance
Our model allows also determining normalized electric input admittance of piezoelectric disc resonator for different rates of metallization. Figure.2 and Figure 3 displays normalized electric input admittance of a MEMS PZT5A resonator as a function of normalized frequency of contour modes for metallization rate 50% and 60% respectively.
3.3. Mechanical and Elctrical Displacements Profiles
For two values of resonant frequencies and, electrical and mechanical displacement profiles are calculated and presented in a 3D configuration for 80% rate of metallization. Figure 4 and Figure 5 shows the displacements profiles in radial and axial direction respectively. In Figure 6, electrical displacement profiles were presented.
Figure 4. Radial components of displacement profiles for and , M=N=15.
Figure 5. Axial components of displacement profiles for kHz and kHz, M=N=15.
Figure 6. Electrical displacement profiles into PZT5A resonator for: kHz and kHz, M=N=15.
We can denote that the rate of metallization is an important factor in controlling the frequency variations of the resonator. Therefore, resonant and anti-resonant frequencies increase with the decrease of the length on the electrode.
4. Conclusion
Along this paper, Legendre polynomial method was successfully extended for studying the vibrational characteristics of a piezoelectric MEMS resonator partially metallized with centralized metallization. As an optimization tool, Legendre polynomial approach provides satisfied results such as resonant and anti-resonant frequencies, normalized electric impedance, electromechanical coupling coefficient, mechanical and electrical displacement profiles in three dimensional configurations. The formulation is easy to implement in which boundary conditions and electric excitation source are directly incorporated in the equation of motion.
Appendix
The analytical method used for calculating normalized admittance is developed for a thin disc (H<<R) made from piezoelectrique material.
where:
where the matrix elements are expressed as follow:
; ;
Where:
; ; ; ;
where and are the Bessel functions of the first kind, and are the elastic compliance constants at constant electrical field. and are the piezoelectric constants, is the Poisson ratio and is the dielectric constant at constant stress.
References