American Journal of Mechanics and Applications
Volume 4, Issue 1, November 2016, Pages: 1-9

Modeling of MEMS Resonator Piezoelectric Disc Partially Covered with Electrodes

Ismail Naciri1, *, Lahoucine Elmaimouni1, *, Jean-Etienne Lefebvre2, Faniry Emilson Ratolojanahary3, Mohamed Rguiti4, Tadeusz Gryba2

1Laboratoire Sciences Ingénierie et Energie, Energie Renouvelable, Microsystèmes Acoustique et Micanique, Polydisciplinary Faculty of Ouarzazate, Ibn Zohr University, Morocco

2The Institute of Electronics, Microelectronics and Nanotechnology, Opto-Acousto-Electronic Department, University of Valenciennes, France

3Laboratory of applied Physics, Fianarantsoa University, Madagascar

4Laboratoire des Matériaux Céramiques et procédés Associés, Université de Valenciennes, Maubeuge, France

Email address:

(I. Naciri)
(L. Elmaimouni)

*Corresponding author

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]:


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].


Along this paper, we adopt the following change of variables:

Normalized mechanical displacement components: , , .

 and .

Equations expressed in (1) can be written as:





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.


The equations expressed in (4) are directly incorporated in the constitutive equations of our problem.

The normalized stress components are expressed by





The electric displacement components:



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:




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:




Mechanical and electrical displacement components for area M and N are defined respectively as follow:





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:




; ;

 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:


Where  denote the static capacitance of the disc resonator.

From equation 12, we obtain:


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.





3. Numerical Results

We start our calculations with properties in table (1), the computer program is established using Matlab software.

Table 1. Physical properties of PZT5A.

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: .

Table 2. Normalized resonant frequencies of the first five contour modes, M=N=15.

Table 3. Normalized anti-resonant frequencies of the first five contour modes, M=N=15.

Table 4. Electromechanical coupling coefficient of the first five contour mode, M=N=15.

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.

Figure 2. Normalized electric input admittance of PZT5A resonator as a function of normalized frequency of contour modes with. metallization rate of 50%, M=N=15.

Figure 3. Normalized electric input admittance of PZT5A resonator as a function of normalized frequency of contour modes with metallization rate of 60%, M=N=15.

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.


The analytical method used for calculating normalized admittance is developed for a thin disc (H<<R) made from piezoelectrique material.


 where the matrix elements are expressed as follow:

; ;


; ; ; ;

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.


  1. M. Lutz, A. Partridge, P. Gupta, N. Buchan, E. Klaassen, J. McDonald, K. Petersen. "MEMS oscillators for high volume commercial applications", Int. S. State Sen Actuat. Microsyst. (14th Conf. Transducers), pp. 49-52, 2007.
  2. M. H. Tsai, Y. C. Liu, C.M. Sun, C. Wang, C. W. Cheng, W. Fang, "3-Axis CMOS-MEMS accelerometer with vertically integrated fully-differential sensing electrodes", 16th Int. Solid-State S. Actu. Micr. Conf. (TRANSDUCERS), pp. 811-814, 2011.
  3. M. Rinaldi, C. Zuniga, G. Piazza, "Ss-DNA functionalized array of AlN contourmode NEMS resonant sensors with single CMOS multiplexed oscillator for subppb detection of volatile organic chemicals", IEEE 24th Int. Conf. MEMS, , pp. 976–979, 2011
  4. T.W. Secord, H.H. Asada, A variable stiffness PZT actuator having tunable resonant frequencies, IEEE Trans. Robot. Vol.26 (6), pp. 993–1005, 2010.
  5. D.W. Greve, J.J. Neumann, I.J. Oppenheim, S.P. Pessiki, D. Ozevin, Robust capacitive MEMS ultrasonics transducers for liquid immersion, IEEE Symp. Ultrason. Vol.1, pp. 581–584, 2003.
  6. D.K. Agrawal, P. Thiruvenkatanathan, J. Yan, A.A. Seshia, Electrically coupled MEMS oscillators, in: Joint Conf. IEEE Int. Freq. Control Eur. Freq. Time Forum (FCS), pp. 1-5, 2011.
  7. H. F. Tiersten, Linear piezoelectric plate vibration, Plenum, New York, 1969.
  8. IEEE Standard on Piezoelectricity, ANSI-IEEE Std. 176, IEEE New York, 1987.
  9. E. P Eer Nisse, "Variational Method for Electroelastic Vibration Analysis," IEEE Trans. Sonics. Ultra. Vol. 14(4), pp. 153–160, 1967.
  10. C. H. Huang and C. C. Ma, "Vibration Characteristics for Piezoelectric Cylinders Using Amplitude- Fluctuation Electronic Speckle Pattern Interferometry" AIAA JL Vol. 36(12), 1998.
  11. Kharouf, N, and Heyliger, P. R., "Axisymmetric Free Vibration of Homogeneous and Laminated Piezoelectric Cylinders," J. Sound. Vib., Vol. 174(4), pp. 539–561, 1994.
  12. H. A Kunkel, S. Locke, and B. Pikeroen, "Finite-Element Analysis of Vibrational Modes in Piezoelectric Ceramics Disks," IEEE. Trans. Ultrason. Ferr. Freq. Cont, Vol. 37(4), pp. 316–328, 1990.
  13. N. Guo, P. Cawley and D. Hitchings, "The Finite Element Analysis of the Vibration Characteristics of Piezoelectric Disks," J Sound. Vib., Vol. 159(1), pp. 115-138, 1992.
  14. L. Elmaimouni, J. E Lefebvre, V. Zhang, and T. Gryba, "Guided waves in radially graded cylinders: a polynomial approach", NDT & E Int, Vol. 38, pp.344-353, 2005.
  15. J. E. Lefebvre, V. Zhang, J. Gazalet, T. Gryba and V. Sadaune, "Acoustic Waves Propagation in Continuous Functionally Graded Plates: An Extension of the Legendre Polynomial Approach"IEEE. Trans. Ultrason. Ferr. Freq. Cont, Vol. 48(5), pp. 1332-1340, 2001.
  16. J. Yu, J. E. Lefebvre and L. Elmaimouni, Toroidal wave in multilayered spherical curved plates, Journal of Sound and Vibration, Vol. 332 (11), pp. 2816-2830, 2013.
  17. L. Elmaimouni, J. E. Lefebvre, F. E. Ratolojanahary, A. Raherison, T. Gryba and J. Carlier, "Modal analysis and harmonic response of resonators: an extension of a mapped orthogonal functions technique" Wave Motion, Vol. 48(1), pp. 93-104, 2011.
  18. L. Elmaimouni, J. E. Lefebvre, F. E. Ratolojanahary, A. Raherison, T. Gryba and J. Carlier, "Modal analysis and harmonic response of resonators: an extension of a mapped orthogonal functions technique" Wave Motion, Vol. 48(1), pp. 93-104 (2011).
  19. L. Elmaimouni, J. E. Lefebvre, F. E. Ratolojanahary A. Raherison, B. Bahani and T. Gryba, "Polynomial approach modeling of resonator piezoelectric disc" Key Engineering Materials, Dynamics of the structures and Non Destructive testing, Vol. 1294 (482), pp. 11-20, 2011.
  20. P. M. Rabotovao, F. E. Ratolojanahary, J. E. Lefebvre, A. Raherison, L. Elmaimouni, T. Gryba, and J. G. Yu, "Modeling of high contrast partially electroded resonators by means of a polynomial approach", J. Applied Physics, Vol. 114 (12), pp.124502, 2013.
  21. L. Elmaimouni, F. E. Ratolojanahary, J. E. Lefebvre, J.G. Yu, A. Raherison and T. Gryba, "Modeling of MEMS resonator piezoelectric disc by means of an equicharge current source method", Ultra, Vol. 53(7), pp. 1270-1279, 2013.
  22. B. A Auld, "Acoustic Fields and Waves in Solids" Krieger Publishing Company, Malabar, Florida, 1990.
  23. D. Royer et E Dieulesaint, Ondes élastiques dans les solides, Masson, Paris 1994.
  24. J. Yu, J. E. Lefebvre, Y. Guo, and L. Elmaimouni "Wave Propagartion in the Circumferential Direction of General Multilayered Piezoelectric Cylindrical Plates". IEEE. Trans. Ultrason. Ferr. Freq. Cont, Vol. 59, pp. 0885-3010, 2012.
  25. N. Guo, "The vibration characteristics of piezoelectric discs", Dissertation, Department of Mechanical Engineering, Imperial College of Science, Technology and Medicine, London, 1989.

Article Tools
Follow on us
Science Publishing Group
NEW YORK, NY 10018
Tel: (001)347-688-8931