# Analysis of nonlinear dynamic behavior and pull-in prediction of micro circular plate actuator

## Chin-Chia Liu1, Chien-Hung Liu2

1Department of Industrial Education and Technology, National Changhua University of Education, Bao-Shan Campus: Number 2, Shi-Da Road, Changhua, 500, Taiwan, R.O.C.

2Department of Mechanical Engineering, National Chung Hsing University, Taichung, Taiwan, R.O.C.

1Corresponding author

Journal of Vibroengineering, Vol. 16, Issue 4, 2014, p. 2062-2069.
Received 8 October 2013; received in revised form 11 December 2013; accepted 18 December 2013; published 30 June 2014

Copyright © 2014 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 69
Abstract.

The dynamic behavior of micro circular plate electrostatic devices is not easily analyzed using traditional methods such as perturbation theory or Galerkin approach method due to the complexity of the interactions among the electrostatic coupling effect, the residual stress and the nonlinear electrostatic force. Accordingly, the present study proposes a approach for analyzing the dynamic response of such devices using a hybrid numerical scheme comprising the differential transformation method and the finite difference method. The feasibility of the proposed approach is demonstrated by modeling the dynamic response of a micro circular plate actuated by a DC voltage. The numerical results for the pull-in voltage are found to deviate by no more than 0.27 % from those derived in the literature using various computational methods. Thus, the basic validity of the hybrid numerical scheme is confirmed. Moreover, the effectiveness of a combined DC/AC loading scheme in driving the micro circular actuator is examined. It is shown that the use of an AC actuating voltage in addition to the DC driving voltage provides an effective means of tuning the dynamic response of the micro circular plate.

Keywords: micro circular plate, pull-in voltage, MEMS, hybrid method, differential transformation, finite difference method.

#### 1. Introduction

Micro-electro-mechanical systems (MEMS) devices have found widespread use throughout industry for such applications as accelerometers and pressure sensors [1], micro-scale actuators [2], electrostatic rotary comb actuators [3], and so on. Various schemes have been proposed for driving MEMS devices, including electrostatic actuation [4], thermal actuation, piezoelectric actuation [5] and electromagnetic actuation [6]. Among these various methods, electrostatic actuation is the most commonly preferred due to its simplicity and high efficiency [7]. Actually, in the electrostatic actuation of a micro-structure system, the electrostatic force is produced from the voltages of two electrodes. If the electrostatic force is greater than the elastic restoring force of the micro-structure system, this represents an unstable phenomenon [8], and the two electrodes attract and come into contact with each other suddenly. The critical value of the voltage is defined as the pull-in voltage, which has a tremendous influence on the electrostatic device. For example, the electrostatic device is regarded as an actuator when the operation voltage is greater than the pull-in voltage and the upper electrode can be attracted to the fixed bottom electrode very quickly; therefore, the pull-in voltage limits the operation range of the actuator. The pull-in behavior phenomenon, however, can be used in the design of such components as switches [9] and relays; moreover, it could be used to measure Young’s modulus and residual stress values [8]. Hence, in designing electrostatically-actuated MEMS devices, a thorough understanding of the dynamic behavior of the device structure is required such that the pull-in event can be either avoided or induced, as required. Abdel-Rahman et al [10] proposed a nonlinear model for characterizing the mechanical behavior of electrostatically-actuated micro-beams Hung [11] examined the feasibility of using the leveraged bending and strain-stiffening methods to extend the travel range of analog-tuned electrostatic actuators prior to pull-in.

Various researchers have investigated the use of hybrid DC/AC schemes in driving electrostatically-actuated MEMS devices [12-15]. Consider the lumped model of a parallel plate capacitor shown in Fig. 1, in which the initial gap height between the two plates is denoted by $G$. The dynamic behavior of the system is governed by the micro electrostatic force, the mechanical restoring force and the damping force together govern the dynamic behavior of micro system. The micro electrostatic force is determined by the driving voltage $V\left(t\right)$, which comprises both a DC voltage component and an AC voltage component. The DC component offsets the upper plate to a new equilibrium position, while the AC component causes the plate to vibrate around this newly established equilibrium position. Younis et al [13] proposed a novel RF MEMS switch actuated using a combined DC/AC loading scheme. It was shown that an appropriate specification of the magnitude and frequency of the AC voltage enabled the driving voltage required to induce the pull-in event to be significantly reduced. Younis [14] examined the dynamic behavior of micro-beams subject to combined DC/AC loading and derived analytical expressions for the micro-beam motion under primary resonance conditions. In general, MEMS devices with a circular plate configuration have a better structure flexibility than those with a rectangular plate configuration due to a lack of sharp edges and a lower residual stress after multiple depositions [15].

Fig. 1. Lumped model of parallel-plate capacitor

Differential transformation theory was originally proposed by Zhao as a means of solving linear and nonlinear initial value problems in the circuit analysis domain. However, in more recent years, researchers have extended its use to the analysis of a variety of initial value problems in the mechanical engineering field. For example, Chen et al. [16] showed that the combined differential transformation and finite difference method provides a precise and computationally-efficient means of analyzing the nonlinear dynamic behavior of fixed-fixed micro-beams. The same group also used the hybrid method to analyze the nonlinear dynamic response of an electrostatically-actuated micro circular plate subject to both effects of residual stress and a uniform hydrostatic pressure on the upper surface [17-18, 20]. A numerical investigation was performed into the entropy generated within a mixed convection flow with viscous dissipation effects in a parallel-plate vertical channel using differential transformation method by Chen et al. [19].

#### 2. Differential transformation theory

The basic principles of the differential transformation method are introduced below. If $x\left(t\right)$ is an analyzable function in time domain $T$, a definition of the differential transformation of $x$ at $t={t}_{0}$ in the $K$ domain is:

(1)

where $k$ belongs to the set of non-negative integers denoted as the $K$ domain, $X\left(k;{t}_{0}\right)$ is the transformed function in the transformation domain, otherwise called the spectrum of $x\left(t\right)$ at $t={t}_{0}$ in the $K$ domain, $M\left(k\right)$ is the weighting factor, and $q\left(t\right)$ is regarded as a kernel corresponding to $x\left(t\right)$. Both $M\left(k\right)$ and $q\left(t\right)$ are non-zero and $q\left(t\right)$ is an analyzable function in time domain $T$. Therefore, the differential inverse transformation of $X\left(k;{t}_{0}\right)$ can be described as:

(2)

If $M\left(k\right)={H}^{k}/k!$ and $q\left(t\right)=\text{1,}$ where $H$ is the time interval. Let ${t}_{0}=\text{0;}$ Eq. (1) then becomes:

(3)

The differential inverse transformation of $X\left(k\right)$ can then be expressed as below by Eq. (2):

(4)

Substituting Eq. (3) into Eq. (4) gives:

(5)

Eq. (5) can be derived by Taylor series expansion. Therefore, the main basic operation properties of the differential transform are as listed below:

(a) Linearity operation:

(6)
$T\left[\alpha x\left(t\right)+\beta y\left(t\right)\right]=\alpha X\left(k\right)+\beta Y\left(k\right),$

where $T$ denotes the differential transform and $\alpha$ and $\beta$ are any real number.

(b) Differential operation:

(7)
$T\left[\frac{{d}^{n}x\left(t\right)}{d{t}^{n}}\right]=\frac{\left(k+n\right)!}{k!{H}^{n}}X\left(k+n\right),$

where $T$ denotes the differential transform and $n$ is the order of differentiation [16-19].

#### 3. Dynamic modeling of micro circular plate

The analysis performed in this study considers the micro circular plate shown in Fig. 2. As shown, the upper plate is actuated by a driving voltage $V\left(t\right)={V}_{DC}+{V}_{AC}\mathrm{s}\mathrm{i}\mathrm{n}\left(\omega t\right)$, where ${V}_{DC}$ is the DC polarization voltage, ${V}_{AC}$ is the magnitude of AC voltage, and $\omega$ is the excitation frequency. At a certain critical value of the driving voltage, the magnitude of the electrostatic force between the upper and lower plates exceeds the restoring force within the upper, movable plate, and thus a pull-in occurs.

Fig. 2. Schematic illustration showing electrostatic actuation of micro circular plate

The governing equation is therefore given as [12]:

(8)
$\rho h\frac{{\partial }^{2}w}{\partial {t}^{2}}+D\left(\frac{{\partial }^{2}}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial }{\partial r}\right)\left(\frac{{\partial }^{2}w}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial w}{\partial r}\right)-{T}_{r}h\left(\frac{{\partial }^{2}w}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial w}{\partial r}\right)=\frac{{\epsilon }_{0}V{\left(t\right)}^{2}}{2{\left(G-w\right)}^{2}}+{S}_{0},$

where ${\epsilon }_{0}$, $h$ and $G$ are represent the permittivity of free space, the thickness of the micro circular plate, and the initial gap height between the upper and lower plates, respectively. In addition, $V\left(t\right)$ is the voltage between the two plates, (i.e., ), $\rho$ is the density of the micro circular plate, and $w$ is the transverse deflection of the micro circular plate at a distance $r$ from the center of the plate. Note that $w$ is a function only of the position $r$ and the time $t$, i.e., $w=w\left(r,t\right)$. In other words, the symmetry transverse deflection of micro circular plate is irrelevant to polar coordinate $\theta$. Finally, ${S}_{0}$ is the uniform hydrostatic pressure acting on the upper surface of the plate, ${T}_{r}$ is the residual stress within the plate, and $D$ is the flexural rigidity of the plate, i.e.:

(9)
$D=\frac{E{h}^{3}}{12\left(1-{\upsilon }^{2}\right)},$

where $\upsilon$ and are the Poisson ratio and Young’s modulus of the upper circular plate, respectively. The boundary conditions for the governing equation of motion of the micro circular plate for are defined as follows:

(10)

where $R$ is the radius of the micro circular plate. (Note that Eq. (10) assumes that the plate is clamped at its edges). Finally, the initial condition is defined as:

(11)
$w\left(r,0\right)=\frac{\partial w\left(r,0\right)}{\partial t}=0.$

#### 4. Application of hybrid method to solution of governing equation of motion

For analytical convenience, let the the transverse displacement of the micro circular plate $w$ be normalized with respect to the initial gap between the electrodes. Furthermore, let the radial distance position $r$ be normalized with respect to the radius of the micro circular plate, and let time $t$ be normalized with respect to a time constant ${T}_{n}$, where ${T}_{n}$ is defined as ${T}_{n}=\sqrt{\rho h{R}^{4}/D}$. Finally, let the excitation frequency $\omega$ be normalized by taking the product of $\omega$ and the time constant ${T}_{n}$, i.e.:

(12)

Let the following parameters be defined:

(13)

Substituting Eqs. (12) and (13) into Eqs. (8), (10) and (11), the dimensionless governing equation of motion for the micro circular plate is obtained as:

(14)
$\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {\stackrel{-}{t}}^{2}}+\frac{{\partial }^{4}\stackrel{-}{w}}{\partial {\stackrel{-}{r}}^{4}}+\frac{2}{\stackrel{-}{r}}\frac{{\partial }^{3}\stackrel{-}{w}}{\partial {\stackrel{-}{r}}^{3}}-\frac{1}{{\stackrel{-}{r}}^{2}}\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {\stackrel{-}{r}}^{2}}+\frac{1}{{\stackrel{-}{r}}^{3}}\frac{\partial \stackrel{-}{w}}{\partial \stackrel{-}{r}}-{\stackrel{-}{T}}_{r}\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {\stackrel{-}{r}}^{2}}-{\stackrel{-}{T}}_{r}\frac{1}{\stackrel{-}{r}}\frac{\partial \stackrel{-}{w}}{\partial \stackrel{-}{r}}=\frac{{\alpha }_{1}{\left[{V}_{DC}+{V}_{AC}\mathrm{sin}\left(\stackrel{-}{\omega }\stackrel{-}{t}\right)\right]}^{2}}{{\left(1-\stackrel{-}{w}\right)}^{2}}+\stackrel{-}{{S}_{0}}.$

The corresponding dimensionless boundary conditions are given as follows:

(15)

Finally, the initial condition is equal to:

(16)
$\stackrel{-}{w}\left(\stackrel{-}{r},0\right)=\frac{\partial \stackrel{-}{w}\left(\stackrel{-}{r},0\right)}{\partial \stackrel{-}{t}}=0.$

The nonlinear governing equation of motion for the micro circular plate (Eq. (14)) is solved using the hybrid differential transformation and finite difference method described in [16-19]. The solution procedure commences by discretizing the equation of motion with respect to the time domain $t$ using the differential transformation method. In accordance with the hybrid numerical method, the transverse displacement of the micro plate is then discretized spatially in the radial direction using the finite difference approximation method based on fourth-order and second-order accurate central difference formulae.

#### 5. Numerical results and discussion

In performing the comparison, two models of the micro circular actuator are considered (see Table 1). As shown in Table 2, the pull-in voltage of the Model 1 micro circular plate is found to be 363 V when using the hybrid differential transformation and finite difference method. The equivalent values of the pull-in voltage computed using Osternberg’s model [18] and the reduced-order model [18] are 362 V and 364 V, respectively. In other words, the predicted value of the pull-voltage obtained using the hybrid method deviates by just 0.27 % from the results presented in [18].

Table 1. Material and geometry parameters of micro circular plate models

 Parameters Model 1 Model 2 Young’s modulus ($E$ / GPa) 169 130 Poisson’s Ratio ($\nu$) 0.3 0.23 Density ($\rho$ / Kg/m3) 2.33×103 2.33×103 Permittivity of free space (${\epsilon }_{0}$ / F/m) 8.8542×10-12 8.8542×10-12 Thickness of circular plate ($h$ / μm) 20 1 Initial gap ($G$ / μm) 1 1 Radius of circular plate ($R$ / μm) 250 200

Table 2. Comparison of present analytical results and literature results for pull-in voltages

 Analytical results Deviation Hybrid method (H.M.) (Model 1) Osterberg’s model [15] Reduced-order model [15] $∆{e}_{1}$ (%) $\text{∆}{\text{e}}_{\text{2}}$ (%) Pull-in Voltage (V) 363 362 364 0.27 0.27

Figure 3 shows the variation of the center-point deflection of the Model 2 circular plate over time given AC voltages of 1.61 V, 1.62 V and 1.63 V, respectively. Note that the DC voltage is equal to 4.5 V, the residual stress is equal to ${T}_{r}=\text{0}$, and the AC voltage frequency is equal to $\stackrel{-}{\omega }=1$. The results clearly show that for an AC voltage lower voltage of 1.62 V, the micro circular plate oscillates in a stable manner about the equilibrium deflection point. However, when the AC voltage is increased to 1.63 V, the pull-in phenomenon occurs, and the micro circular plate makes transient contact with the lower electrode.

Fig. 3. Variation of dimensionless center-point displacement over time for AC voltages of 1.61 V, 1.62 V and 1.63 V, respectively, (note that DC voltage $=$ 4.5 V and $\stackrel{-}{\omega }=1$)

Figure 4 shows the deformation over time of the Model 2 micro circular plate given an AC voltage of 0.5 V, a dimensionless excitation frequency of $\stackrel{-}{\omega }=\text{1}/\text{2}$, and a DC voltage of 2 V. As stated previously, the transverse displacement of the micro circular plate depends only on radial position and the time. In other words, the displacement is independent of the polar coordinate $\theta$. As a result, the micro circular plate oscillates in a stable fashion.

Fig. 4. Deformation of micro circular plate given DC voltage of 2 V and AC voltage of 0.5 V, (note that $\stackrel{-}{\omega }=1/2$)

#### 6. Conclusions

This study has analyzed the nonlinear dynamic behavior of an electrostatically-actuated micro circular plate using a hybrid numerical scheme comprising the differential transformation method and the finite difference method. The validity of the proposed scheme has been confirmed by comparing the predicted value of the pull-in voltage for the plate with the results presented in the literature. Finally, it has been shown that the use of an AC actuating voltage with an appropriate magnitude in addition to the DC driving voltage provides an effective means of tuning the dynamic response of the micro circular plate. Overall, the numerical results presented in this study show that the hybrid differential transformation method and finite difference method provides an accurate and computationally-efficient means of analyzing the nonlinear dynamic behavior of the micro circular plate structures used in many of today’s MEMS-based actuator systems.

#### Acknowledgements

The authors gratefully acknowledge the financial support provided to this study by the National Science Council of Taiwan under Grant Number NSC 101-2221-E-018-007.

#### References

1. Bosc J. M., Guo Y., Sarihan V., Lee T. Accelerated life testing for micro-machined chemical sensors. IEEE Transactions on Reliability, Vol. 47, Issue 2, 1998, p. 135-141. [Search CrossRef]
2. Zhu Y., Espinosa H. D. Effect of temperature on capacitive RF MEMS switch performance – a coupled-field analysis. Journal of Micromechanics and Microengineering, Vol. 14, Issue 8, 2004, p. 1270-1279. [Search CrossRef]
3. Maluf N. I., Reay R. J., Kovacs G. T. A. High-voltage devices and circuits fabricated using foundry CMOS for use with electrostatic MEM actuators. Sensors and Actuators A, Vol. 52, 1996, p. 187-192. [Search CrossRef]
4. Rezazadeh G., Tahmasebi A., Zubstov M. Application of piezoelectric layers in electrostatic MEM actuators: controlling of pull-in voltage. Microsystem Technologies, Vol. 12, Issue 12, 2006, p. 1163-1170. [Search CrossRef]
5. Luharuka R., LeBlanc S., Bintoro J. S., Berthelot Y. H., Hesketh P. J. Simulated and experimental dynamic response characterization of an electromagnetic microvalve. Sensors and Actuators A, Vol. 143, 2008, p. 399-408. [Search CrossRef]
6. Tong P., Huang W. Large deflection of thin plates in pressure sensor applications. Journal of Applied Mechanics, Vol. 69, Issue 6, 2002, p. 785-789. [Search CrossRef]
7. Varadan V. M., Vinoy K. J., Jose K. A. RF MMES and their applications. Wiley, New York, 2003. [Search CrossRef]
8. Osterberg P. M., Senturia S. D. M-test: a test chip for MEMS material property measurement using electrostatically actuated test structures. Journal of Microelectromechanical Systems, Vol. 6, Issue 2, 1997, p. 107-118. [Search CrossRef]
9. Wong J. E., Lang J. H., Schmidt M. A. An electrostatically actuated MEMS switch for power applications. The IEEE MEMS Conference, Miyazaki, Japan, 2000, p. 633-638. [Search CrossRef]
10. Abdel-Rahman E. M., Younis M. I., Nayfeh A. H. Characterization of the mechanical behavior of an electrically actuated microbeam. Journal of Micromechanics and Microengineering, Vol. 12, Issue 6, 2002, p. 759-766. [Search CrossRef]
11. Hung E. S., Senturia S. D. Extending the travel range of analog-tuned electrostatic actuators. Journal of Microelectromeehanieal Systems, Vol. 8, Issue 4, 1999, p. 497-505. [Search CrossRef]
12. Nayfeh A. H., Younis M. I., Abdel-Rahman E. M. Reduced-order models for MEMS applications. Nonlinear Dynamics, Vol. 41, Issue 1-3, 2005, p. 211-236. [Search CrossRef]
13. Younis M. I., Abdel-Rahman E. M., Nayfeh A. H. Dynamic simulations of a novel RF MEMS switch. NSTI-Nanotech, Vol. 2, 2004, p. 287-290. [Search CrossRef]
14. Younis M. I. Investigation of the mechanical behavior of microbeam-based MEMS devices. Master of Science Thesis, Virginia Polytechnic Institute and State University, 2001. [Search CrossRef]
15. Vogl G. W. V., Nayfeh A. A reduced-order model for electrically actuated clamped circular plates. Journal of Micromechanics and Microengineering, Vol. 15, Issue 4, 2005, p. 684-690. [Search CrossRef]
16. Chen C. K., Lai H. Y., Liu C. C. Application of hybrid differential transformation / finite difference method to nonlinear analysis of micro fixed-fixed beam. Microsystem Technologies, Vol. 15, Issue 6, 2009, p. 813-820. [Search CrossRef]
17. Chen C. K., Lai H. Y., Liu C. C. Nonlinear micro circular plate analysis using hybrid differential transformation / finite difference method. CMES: Computer Modeling in Engineering & Sciences, Vol. 40, Issue 2, 2009, p. 155-174. [Search CrossRef]
18. Liu C. C., Yang S. C., Chen C. K. Nonlinear dynamic analysis of micro cantilever beam under electrostatic loading. Journal of Mechanics, Vol. 28, 2012, p. 559-566. [Search CrossRef]
19. Chen C. K., Lai H. Y., Liu C. C. Numerical analysis of entropy generation in mixed convection flow with viscous dissipation effects in vertical channel. International Communications in Heat and Mass Transfer, Vol. 38, 2011, p. 285-290. [Search CrossRef]
20. Liu C. C., Chen C. K. Modeling and simulation of nonlinear micro-electromechanical circular plate. Smart Science, Vol. 1, Issue 1, 2013, p. 59-63. [Search CrossRef]