Modelling and Simulation in Engineering
Volume 2016 (2016), Article ID 8543616, 6 pages
http://dx.doi.org/10.1155/2016/8543616

Research Article

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

P. Di Maida and G. Bianchi

Dipartimento di Scienze e Metodi dell’Ingegneria (DISMI), Universitá degli Studi di Modena e Reggio Emilia, Via G. Amendola 2, 42122 Reggio Emilia, Italy

Received 4 December 2015; Revised 10 May 2016; Accepted 5 October 2016

Academic Editor: Julius Kaplunov

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.

Abstract.

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 [1–5]. 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 [11–14]. Failure mechanisms of MEMS include cracking [15–18], peeling of the cantilever [19–21], stiction to the substrate [22–24], 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 variables

Then the governing equation for the cantilever reads

whereandare the Casimir and the electrostatic line-load, being

Here, 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) as

where 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 read

Let us define ; thus (4) further reduces to

with and the BCs

It 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 give

It 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) give

which shows that the situation is not relevant in this problem. Besides, it follows

and evaluation at lends

Consideration of the sign for and yields the inequality

Figure 1: Micro-cantilever switch.

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 equations

which 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 derivative

Making use of the derivative matrix (see [33, Chap.7]) we have

and (16) and (17) may be rewritten through the derivative matrix as

where 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 have

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

Figure 2: Displacement at and .

Figure 3: Plot of versus according to the FE approximation: this product should return the constant .

Figure 4: Plot of versus according to the pseudo-spectral approximation: this product should return the constant .

Figure 5: Slope (solid), bending moment (dashed), and shearing force (dotted) near pull-in. Graphs are indistinguishable between the pseudo-spectral and the FE solution.

Figure 6: Line load density on the cantilever beam near pull-in according to the FE approximation ().

Figure 7: Line load density on the cantilever beam near pull-in according to the FE approximation ().

Figure 8: Line load density on the cantilever beam near pull-in according to the pseudo-spectral approximation ().

Figure 9: Line load density on the cantilever beam near pull-in according to the pseudo-spectral approximation ().

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 [34–36] and beam-plates [37, 38] or to include viscoelastic effects [39–42].

Competing Interests

The authors declare that they have no competing interests.

References

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.

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.

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.

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.

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.

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.

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.

J. Kaplunov and A. Nobili, “Multi-parametric analysis of strongly inhomogeneous periodic waveguideswith internal cutoff frequencies,” Mathematical Methods in the Applied Sciences, 2016.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

C. H. Mastrangelo, Suppression of Stiction in MEMS, vol. 605 of MRS Proceedings, Cambridge University Press, Cambridge, UK, 1999.

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.

Z. Yapu, “Stiction and anti-stiction in MEMS and NEMS,” Acta Mechanica Sinica, vol. 19, no. 1, pp. 1–10, 2003.

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.

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.

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.

G. Napoli and A. Nobili, “Mechanically induced Helfrich-Hurault effect in lamellar systems,” Physical Review E, vol. 80, no. 3, Article ID 031710, 2009.

J. P. Boyd, Chebyshev and Fourier Spectral Methods, Dover Publications, 2nd edition, 2000.

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.

A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations, CRC Press, Boca Raton, Fla, USA, 1st edition, 1995.

P. L. Sachdev, Nonlinear Ordinary Differential Equations and Their Applications, Marcel Dekker, 1991.

D. Funaro, Polynomial Approximation of Differential Equations, vol. 8 of Lecture Notes in Physics, Springer, 1992.

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.

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.

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.

R. Barretta and R. Luciano, “Exact solutions of isotropic viscoelastic functionally graded Kirchhoff plates,” Composite Structures, vol. 118, no. 1, pp. 448–454, 2014.

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.

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.

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.

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.

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.