Modelling and Simulation in Engineering

Modelling and Simulation in Engineering / 2016 / Article
Special Issue

Structural Modelling at the Micro-, Meso-, and Nanoscales

View this Special Issue

Research Article | Open Access

Volume 2016 |Article ID 8543616 |

P. Di Maida, G. Bianchi, "Numerical Investigation of Pull-In Instability in a Micro-Switch MEMS Device through the Pseudo-Spectral Method", Modelling and Simulation in Engineering, vol. 2016, Article ID 8543616, 6 pages, 2016.

Numerical Investigation of Pull-In Instability in a Micro-Switch MEMS Device through the Pseudo-Spectral Method

Academic Editor: Julius Kaplunov
Received04 Dec 2015
Revised10 May 2016
Accepted05 Oct 2016
Published08 Nov 2016


A pseudo-spectral approximation is presented to solve the problem of pull-in instability in a cantilever micro-switch. As well known, pull-in instability arises when the acting force reaches a critical threshold beyond which equilibrium is no longer possible. In particular, Coulomb electrostatic force is considered, although the method can be easily generalized to account for fringe as well as Casimir effects. A numerical comparison is presented between a pseudo-spectral and a Finite Element (FE) approximation of the problem, both methods employing the same number of degrees of freedom. It is shown that the pseudo-spectral method appears more effective in accurately approximating the behavior of the cantilever near its tip. This fact is crucial to capturing the threshold voltage on the verge of pull-in. Conversely, the FE approximation presents rapid successions of attracting/repulsing regions along the cantilever, which are not restricted to the near pull-in regime.

1. Introduction

Micro-Electro-Mechanical Systems (MEMS) form a rather diverse and inhomogeneous group of micro-devices aimed at sensing and actuating in a wide array of fields, ranging from mechanical or electronic engineering to chemistry or biology, from micro-mechanics to micro-machining [15]. The manufacturing technology is the common standground for such devices, which heavily relies on the different lithographic techniques borrowed from the technology of micro-electronics. Indeed, MEMS devices are mostly obtained from a silicon substrate. It is observed that MEMS are really “systems” in the sense that they are often made up of several functional parts joint together in the device (like piezo- and magneto-sensors [6]). Among MEMS, micro-switch forms a distinct set with great application potential, with special regard to phase shifters and Radio Frequency MEMS (RFMEMS). They are usually gathered in two groups, namely, capacitor and metal-air-metal switches. Besides, they are further divided according to the actuation method: electrostatic, electrothermal, magnetostatic, and piezoelectric among the most common. A study of magnetoelastic actuated micro-switch is given in [7, 8] for the low-frequency asymptotic analysis of energy scavengers. In this paper, we focus attention on the pull-in instability of a capacitor micro-switch actuated by electrostatic Coulomb force. This particular application has received extensive attention in the literature, owing to the importance of pull-in induced failures in applications. A recent review on the subject can be found in [9]. A theoretical analysis of this problem within the static regime is provided in [10] and references therein. Pull-in voltage in cantilever MEMS have been considered in [1114]. Failure mechanisms of MEMS include cracking [1518], peeling of the cantilever [1921], stiction to the substrate [2224], and temperature [25, 26]. Besides, micropolar theories are often preferred when dealing with micro- and nanodevices to incorporate the scale effect [27, 28]. Spotlight is set on a pseudo-spectral approximation of the problem, which is compared with a Finite Element (FE) solution. Spectral methods belong to the family of Galerkin’s (or Ritz’s) methods [29]. Spectral methods are often divided into two groups, namely, pseudo-spectral or interpolating. The former group enforces the fulfillment of the differential operator at a set of points termed nodes (this is sometimes also named orthogonal collocation). For the latter group, wherein the Galerkin’s method is properly placed, the expansion coefficients are obtained projecting the solution onto the basis set [30].

The paper is structured as follows: Section 2 sets forth the governing equations and the boundary condition for a cantilever. Section 3 introduces the pseudo-spectral method. A numerical comparison with the FE method is illustrated in Section 4. Finally, conclusions are drawn in Section 5.

2. Governing Equations

Let us consider a micro-cantilever switch device subjected to electrostatic attractive force (Figure 1). The micro-cantilever acts as one armor of a capacitor under the electric potential difference . Let denote the distance between the capacitor armors. We consider a plane problem and introduce the transverse displacement for the cantilever. Let us introduce the dimensionless variablesThen the governing equation for the cantilever readswhereandare the Casimir and the electrostatic line-load, beingHere, stands for the electric potential difference acting between the capacitor armors (in the SI this is expressed in volt, i.e., where N stands for newton, m for meter, and C for electric charge, expressed in Coulomb),is the armors width, is the electric permittivity (in vacuum), andis generally a function ofwhich takes into account the fringe effect. For the sake of illustrating the method, we neglect the Casimir force contribution and assume independent of. Then, we can rewrite the governing equation (2) aswhere prime denotes differentiation with respect to and the following driving parameter is obtained:Under the attractive electrostatic force, it is and the boundary conditions (BCs) for the cantilever readLet us define ; thus (4) further reduces towith and the BCsIt is observed that, integrating and making use of the last BC, it may be deduced that which shows that the shearing force is generally positive and it is zero only at . The same argument may be applied to infer that is generally negative, apart from the point where it is zero, and that is generally negative, although it vanishes at . Consequently, is a monotonic decreasing function of and . The nonlinear fourth order ODE (7) may be integrated once [31, §4.2.1] to giveIt is observed that, in the case , (10) falls into the Emdem-Fowler class of nonlinear ODEs, which, in special cases, may admit closed form solutions [31, 32]. Evaluation atand making use of the BCs (8) givewhich shows that the situation is not relevant in this problem. Besides, it followsand evaluation at lends Consideration of the sign for and yields the inequality

3. Pseudo-Spectral Method

The governing equations (7) may be numerically solved through a pseudo-spectral approach [29, 33]. Accordingly, a -degree polynomial function is adopted to approximate the function on the interval . Since a collocation method is adopted, the polynomial will be uniquely determined enforcing (7) at some predetermined points (nodes). This procedure results in a system of nonlinear algebraic equations, which may be solved through standard methods, such as the iterative Newton method. The Jacobian of the system may be supplemented in closed form to the numerical equation solver.

Let be the -degree polynomial approximation of . The collocation set is defined through the first Gauss-Lobatto points where stands for the th degree Legendre polynomial. Equation (7), evaluated at the interior nodes, yields the system of algebraic equationswhich is then supplemented with the boundary conditions (8)The problem is now rewritten in matrix form. To this aim, let the unknown column vector, the vector of the square, and the vector of the th derivativeMaking use of the derivative matrix (see [33, Chap.7]) we haveand (16) and (17) may be rewritten through the derivative matrix aswhere is the identity matrix. Here, is supplemented with the BCs (8); that is, the first row is set to zero apart from the first entry that is set to; the second, the last-but-one and the last rows are replaced by the first and the last rows of , , and , respectively. Then we haveAt the RHS of system (20) we have the zero vector supplemented with the BCs, namely,

4. Numerical Solution

In this section, the pseudo-spectral approximation is compared with a Finite Element solution, both methods employing the same number of degrees of freedom, which corresponds to for the order of the interpolating polynomial. For the FE solution, we introducenodes, each endowed withdegrees of freedom, that is, one translational and the other rotational. In Figure 2, the system of algebraic equations (20) is solved for different values of the driving parameter (voltage difference) with the FE method. It clearly appears that, for , the cantilever dimensionless tip deflection is very close to ; that is, the cantilever is at the verge of pull-in instability. Figure 3 shows the product , which should be constant along the cantilever and equal to , for the FE approximation. It is evident that the quality of the numerical solution rapidly deteriorates near the cantilever tip, which is exactly where best accuracy is demanded to effectively capture the pull-in threshold. Conversely, Figure 4 plots the product for the pseudo-spectral approximation. The comparison between the two plots is a striking example of the effectiveness of this method in this kind of nonlinear problems. Figure 5 plots the slope, bending moment and shearing force along the cantilever beam, and it illustrates that BCs are well captured by the numerical solution, either FE or pseudo-spectral. The plot is obtained by successive differentiation of the displacement field and it is readily available for the spectral method, where polynomial functions are employed. Conversely, obtaining the corresponding curves for the FE approximation needs some extra care for curve fitting of the nodal displacement is first applied, which is then successively differentiated. Figure 6 plots the applied line load density for the FE solution near pull-in, which corresponds to (the negative of) the electrostatic Coulomb force . It is remarkable that the electrostatic force appears highly oscillatory (note that curve-fitting is employed to get a continuous plot), and it attains unphysical negative values. As a comparison, Figure 7 describes the same behavior for , that is well below the threshold value for pull-in. It appears that the electrostatic force is rather poorly approximated by the method, even far from instability. On the contrary, Figures 8 and 9 illustrate the same graphs as obtained from the pseudo-spectral method. The smoothness of the solution is remarkable. The reason for the superior performance of the pseudo-spectral approximation seems to lie in the fact that it is more robust in dealing with little deviations of the cantilever tip displacement on the verge of contact. Indeed, it is well known that pseudo-spectral approximation guarantees high precision and exponential convergence (under suitable assumptions, see [29]), and this feature proves important in smoothly approximating the highly nonlinear behavior of the electrostatic force.

5. Conclusions

In this paper, the pseudo-spectral method is adopted to numerically solve the problem of pull-in instability in a cantilever beam. The beam constitutes one armor of a capacitor, the other armor being represented by a grounded flat surface. Although only Coulomb electrostatic force is considered, the method is easily extended to deal with the fringe effect and the Casimir force. Results may be especially relevant for cantilever micro-switches in MEMS. It is shown that the pseudo-spectral method compares very favorably with an equivalent Finite Element approximation, equivalency being constituted by an equal number of degrees of freedom in the methods. In particular, good approximation for the cantilever deflection near its tip is crucial to capturing the threshold voltage on the verge of pull-in. Indeed, poor approximation leads to very unphysical oscillatory attraction/repulsion forces along the cantilever. It is further shown that the oscillatory behavior is not restricted to the near pull-in regime. Finally, it is emphasized that both methods exactly satisfy the boundary conditions (BCs). It is remarked that the present analysis can be extended to incorporate functionally graded cantilevers [3436] and beam-plates [37, 38] or to include viscoelastic effects [3942].

Competing Interests

The authors declare that they have no competing interests.


  1. V. A. Salomoni, C. E. Majorana, G. M. Giannuzzi, and A. Miliozzi, “Thermal-fluid flow within innovative heat storage concrete systems for solar power plants,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 18, no. 7-8, pp. 969–999, 2008. View at: Publisher Site | Google Scholar
  2. A. Nobili, L. Lanzoni, and A. M. Tarantino, “Experimental investigation and monitoring of a polypropylene-based fiber reinforced concrete road pavement,” Construction and Building Materials, vol. 47, pp. 888–895, 2013. View at: Publisher Site | Google Scholar
  3. V. A. Salomoni, C. E. Majorana, B. Pomaro, G. Xotta, and F. Gramegna, “Macroscale and mesoscale analysis of concrete as a multiphase material for biological shields against nuclear radiation,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 38, no. 5, pp. 518–535, 2014. View at: Publisher Site | Google Scholar
  4. G. Dinelli, G. Belz, C. E. Majorana, and B. A. Schrefler, “Experimental investigation on the use of fly ash for lightweight precast structural elements,” Materials and Structures, vol. 29, no. 194, pp. 632–638, 1996. View at: Google Scholar
  5. P. Bisegna and R. Luciano, “Bounds on the overall properties of composites with debonded frictionless interfaces,” Mechanics of Materials, vol. 28, no. 1–4, pp. 23–32, 1998. View at: Publisher Site | Google Scholar
  6. A. Nobili and A. M. Tarantino, “Magnetostriction of a hard ferromagnetic and elastic thin-film structure,” Mathematics and Mechanics of Solids, vol. 13, no. 2, pp. 95–123, 2008. View at: Publisher Site | Google Scholar | MathSciNet
  7. A. Kudaibergenov, A. Nobili, and L. Prikazchikova, “On low-frequency vibrations of a composite string with contrast properties for energy scavenging fabric devices,” Journal of Mechanics of Materials and Structures, vol. 11, no. 3, pp. 231–243, 2016. View at: Publisher Site | Google Scholar
  8. J. Kaplunov and A. Nobili, “Multi-parametric analysis of strongly inhomogeneous periodic waveguideswith internal cutoff frequencies,” Mathematical Methods in the Applied Sciences, 2016. View at: Publisher Site | Google Scholar
  9. W.-M. Zhang, H. Yan, Z.-K. Peng, and G. Meng, “Electrostatic pull-in instability in MEMS/NEMS: a review,” Sensors and Actuators A: Physical, vol. 214, pp. 187–218, 2014. View at: Publisher Site | Google Scholar
  10. Y. Zhang and Y.-P. Zhao, “Static study of cantilever beam stiction under electrostatic force influence,” Acta Mechanica Solida Sinica, vol. 17, no. 2, pp. 104–112, 2004. View at: Publisher Site | Google Scholar
  11. L. C. Wei, A. B. Mohammad, and N. M. Kassim, “Analytical modeling for determination of pull-in voltage for an electrostatic actuated MEMS cantilever beam,” in Proceedings of the 5th IEEE International Conference on Semiconductor Electronics (ICSE '02), pp. 233–238, IEEE, Penang, Malaysia, December 2002. View at: Google Scholar
  12. S. Chowdhury, M. Ahmadi, and W. C. Miller, “A closed-form model for the pull-in voltage of electrostatically actuated cantilever beams,” Journal of Micromechanics and Microengineering, vol. 15, no. 4, pp. 756–763, 2005. View at: Publisher Site | Google Scholar
  13. H. Sadeghian, G. Rezazadeh, and P. M. Osterberg, “Application of the generalized differential quadrature method to the study of pull-in phenomena of MEMS switches,” Journal of Microelectromechanical Systems, vol. 16, no. 6, pp. 1334–1340, 2007. View at: Publisher Site | Google Scholar
  14. A. Ramezani, A. Alasty, and J. Akbari, “Closed-form solutions of the pull-in instability in nano-cantilevers under electrostatic and intermolecular surface forces,” International Journal of Solids and Structures, vol. 44, no. 14-15, pp. 4925–4941, 2007. View at: Publisher Site | Google Scholar
  15. A. M. Tarantino, “On the finite motions generated by a mode I propagating crack,” Journal of Elasticity, vol. 57, no. 2, pp. 85–103, 1999. View at: Publisher Site | Google Scholar | MathSciNet
  16. A. M. Tarantino, “Nonlinear fracture mechanics for an elastic Bell material,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 50, no. 3, pp. 435–456, 1997. View at: Publisher Site | Google Scholar | MathSciNet
  17. R. Luciano and J. R. Willis, “Hashin-Shtrikman based FE analysis of the elastic behaviour of finite random composite bodies,” International Journal of Fracture, vol. 137, no. 1–4, pp. 261–273, 2006. View at: Publisher Site | Google Scholar
  18. A. Nobili, E. Radi, and L. Lanzoni, “A cracked infinite Kirchhoff plate supported by a two-parameter elastic foundation,” Journal of the European Ceramic Society, vol. 34, no. 11, pp. 2737–2744, 2014. View at: Publisher Site | Google Scholar
  19. L. Lanzoni and E. Radi, “Thermally induced deformations in a partially coated elastic layer,” International Journal of Solids and Structures, vol. 46, no. 6, pp. 1402–1412, 2009. View at: Publisher Site | Google Scholar
  20. V. Guidi, L. Lanzoni, and A. Mazzolari, “Patterning and modeling of mechanically bent silicon plates deformed through coactive stresses,” Thin Solid Films, vol. 520, no. 3, pp. 1074–1079, 2011. View at: Publisher Site | Google Scholar
  21. N. Tullini, A. Tralli, and L. Lanzoni, “Intefacial shear stress analysis of bar and thin film bonded to 2D elastic substrate using a coupled FE-BIE method,” Finite Elements in Analysis and Design, vol. 55, pp. 42–51, 2012. View at: Publisher Site | Google Scholar | MathSciNet
  22. C. H. Mastrangelo, Suppression of Stiction in MEMS, vol. 605 of MRS Proceedings, Cambridge University Press, Cambridge, UK, 1999.
  23. W. Merlijn Van Spengen, R. Puers, and I. De Wolf, “A physical model to predict stiction in MEMS,” Journal of Micromechanics and Microengineering, vol. 12, no. 5, pp. 702–713, 2002. View at: Publisher Site | Google Scholar
  24. Z. Yapu, “Stiction and anti-stiction in MEMS and NEMS,” Acta Mechanica Sinica, vol. 19, no. 1, pp. 1–10, 2003. View at: Publisher Site | Google Scholar
  25. G. Xotta, G. Mazzucco, V. A. Salomoni, C. E. Majorana, and K. J. Willam, “Composite behavior of concrete materials under high temperatures,” International Journal of Solids and Structures, vol. 64, pp. 86–99, 2015. View at: Publisher Site | Google Scholar
  26. F. Marotti de Sciarra and M. Salerno, “On thermodynamic functions in thermoelasticity without energy dissipation,” European Journal of Mechanics—A: Solids, vol. 46, pp. 84–95, 2014. View at: Google Scholar
  27. A. Nobili, “On the generalization of the Timoshenko beam model based on the micropolar linear theory: static case,” Mathematical Problems in Engineering, vol. 2015, Article ID 914357, 8 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet
  28. G. Napoli and A. Nobili, “Mechanically induced Helfrich-Hurault effect in lamellar systems,” Physical Review E, vol. 80, no. 3, Article ID 031710, 2009. View at: Publisher Site | Google Scholar
  29. J. P. Boyd, Chebyshev and Fourier Spectral Methods, Dover Publications, 2nd edition, 2000.
  30. F. Greco and R. Luciano, “A theoretical and numerical stability analysis for composite micro-structures by using homogenization theory,” Composites—Part B: Engineering, vol. 42, no. 3, pp. 382–401, 2011. View at: Publisher Site | Google Scholar
  31. A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations, CRC Press, Boca Raton, Fla, USA, 1st edition, 1995.
  32. P. L. Sachdev, Nonlinear Ordinary Differential Equations and Their Applications, Marcel Dekker, 1991. View at: MathSciNet
  33. D. Funaro, Polynomial Approximation of Differential Equations, vol. 8 of Lecture Notes in Physics, Springer, 1992. View at: MathSciNet
  34. R. Barretta, L. Feo, R. Luciano, and F. Marotti de Sciarra, “Variational formulations for functionally graded nonlocal Bernoulli-Euler nanobeams,” Composite Structures, vol. 129, pp. 80–89, 2015. View at: Publisher Site | Google Scholar
  35. R. Barretta, L. Feo, R. Luciano, and F. Marotti de Sciarra, “A gradient Eringen model for functionally graded nanorods,” Composite Structures, vol. 131, pp. 1124–1131, 2015. View at: Publisher Site | Google Scholar
  36. R. Barretta, L. Feo, and R. Luciano, “Some closed-form solutions of functionally graded beams undergoing nonuniform torsion,” Composite Structures, vol. 123, pp. 132–136, 2015. View at: Publisher Site | Google Scholar
  37. R. Barretta and R. Luciano, “Exact solutions of isotropic viscoelastic functionally graded Kirchhoff plates,” Composite Structures, vol. 118, no. 1, pp. 448–454, 2014. View at: Publisher Site | Google Scholar
  38. A. Apuzzo, R. Barretta, and R. Luciano, “Some analytical solutions of functionally graded Kirchhoff plates,” Composites Part B: Engineering, vol. 68, pp. 266–269, 2015. View at: Publisher Site | Google Scholar
  39. L. Dezi, G. Menditto, and A. M. Tarantino, “Homogeneous structures subjected to repeated structural system changes,” Journal of Engineering Mechanics, vol. 116, no. 8, pp. 1723–1732, 1990. View at: Publisher Site | Google Scholar
  40. L. Dezi and A. M. Tarantino, “Time-dependent analysis of concrete structures with a variable structural system,” ACI Materials Journal, vol. 88, no. 3, pp. 320–324, 1991. View at: Google Scholar
  41. L. Dezi, G. Menditto, and A. M. Tarantino, “Viscoelastic heterogeneous structures with variable structural system,” Journal of Engineering Mechanics, vol. 119, no. 2, pp. 238–250, 1993. View at: Publisher Site | Google Scholar
  42. R. Barretta, L. Feo, and R. Luciano, “Torsion of functionally graded nonlocal viscoelastic circular nanobeams,” Composites Part B: Engineering, vol. 72, pp. 217–222, 2015. View at: Publisher Site | Google Scholar

Copyright © 2016 P. Di Maida and G. Bianchi. 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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles