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 V. Let t denote the distance between the capacitor armors. We consider a plane problem and introduce the transverse displacement y(x) for the cantilever. Let us introduce the dimensionless variables(1)ξ=xL,u=yt.Then the governing equation for the cantilever reads(2)EItL4d4udξ4=qc+qe,whereqcandqeare the Casimir and the electrostatic line-load, being(3)qe=αϵ0bV2t21-u2.Here, V stands for the electric potential difference acting between the capacitor armors (in the SI this is expressed in volt, i.e., V=Nm/C where N stands for newton, m for meter, and C for electric charge, expressed in Coulomb),bis the armors width, ϵ0=8.854·10-12C2N-1m-2 is the electric permittivity (in vacuum), andαis generally a function ofuwhich takes into account the fringe effect. For the sake of illustrating the method, we neglect the Casimir force contribution and assume α independent ofu. Then, we can rewrite the governing equation (2) as(4)uiv=A1-u2,where prime denotes differentiation with respect to ξ and the following driving parameter is obtained:(5)A=αϵ0bL4t3EIV2≥0.Under the attractive electrostatic force, it is 0≤u(ξ)≤umax≤1 and the boundary conditions (BCs) for the cantilever read(6)u0=u′0=u′′1=u′′′1=0.Let us define v(ξ)=1-u(ξ); thus (4) further reduces to(7)viv=-Av2,with 0≤vmin≤v(ξ)≤1 and the BCs(8)v0=1,v′0=v′′1=v′′′1=0.It is observed that, integrating and making use of the last BC, it may be deduced that (9)v′′′ξ=∫ξ1Av2τdτ,which shows that the shearing force is generally positive and it is zero only at ξ=1. The same argument may be applied to infer that v′′ is generally negative, apart from the point ξ=1 where it is zero, and that v′ is generally negative, although it vanishes at ξ=0. Consequently, v is a monotonic decreasing function of ξ and v(1)=vmin. The nonlinear fourth order ODE (7) may be integrated once [31, §4.2.1] to give(10)2v′v′′′-v′′2=2Av+43B.It is observed that, in the case B=0, (10) falls into the Emdem-Fowler class of nonlinear ODEs, which, in special cases, may admit closed form solutions [31, 32]. Evaluation atξ=1and making use of the BCs (8) give(11)B=-32Avmin,which shows that the situation B=0 is not relevant in this problem. Besides, it follows(12)2v′v′′′-v′′2=2A1v-1vmin≤0,and evaluation at ξ=0 lends (13)v′′02=2A1vmin-1.Consideration of the sign for v′ and v′′′ yields the inequality(14)v′′≤2A1vmin-1v.
Micro-cantilever switch.
3. Pseudo-Spectral Method
The governing equations (7) may be numerically solved through a pseudo-spectral approach [29, 33]. Accordingly, a n-degree polynomial function pn(ξ) is adopted to approximate the function u(ξ) on the interval [0,1]. Since a collocation method is adopted, the polynomial will be uniquely determined enforcing (7) at some n+1 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 pn(ξ)∈Pn be the n-degree polynomial approximation of v(ξ). The collocation set is defined through the first n+1 Gauss-Lobatto points (15)GLn=ηi such that η0=-1,ηn=+1,pn′ηi=0,i=1,…,n-1,where pn stands for the nth degree Legendre polynomial. Equation (7), evaluated at the interior nodesηi, yields the system of (n-1) algebraic equations(16)pn2ηipnivηi+A=0,i=1,…,n-1,which is then supplemented with the boundary conditions (8)(17)pnη0=1,pn′η0=pn′′ηn=pn′′′ηn=0.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 mth derivative(18)x=pnηi,y=pn2ηi,xm=pnmηi,i=0,…,n,m∈N.Making use of the (n+1)×(n+1) derivative matrix D (see [33, Chap.7]) we have(19)x′=Dx,x′′=DDx=D2x,…and (16) and (17) may be rewritten through the derivative matrix as(20)yD4~x+AId=o~,where Id is the identity matrix. Here, D4~ is D4supplemented with the BCs (8); that is, the first row is set to zero apart from the first entry that is set to1; the second, the last-but-one and the last rows are replaced by the first and the last rows of D, D2, and D3, respectively. Then we have(21)D4~=10⋯0D0,0⋯D0,nD2,04⋯D2,n4⋮⋱Dn-2,04⋯Dn-2,n4Dn-1,02⋯Dn-1,n2Dn,03Dn,13⋯Dn,n3.At the RHS of system (20) we have the zero vector o supplemented with the BCs, namely,(22)o~=10⋮0.
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 n=20 for the order of the interpolating polynomial. For the FE solution, we introduce10nodes, each endowed with2degrees 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) A with the FE method. It clearly appears that, for A=20, the cantilever dimensionless tip deflection is very close to 1; that is, the cantilever is at the verge of pull-in instability. Figure 3 shows the product -v(iv)v2, which should be constant along the cantilever and equal to A, 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 -v(iv)v2 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 v(iv) for the FE solution near pull-in, which corresponds to (the negative of) the electrostatic Coulomb force -qe. 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 A=12, 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.
Displacement u at A=12,14,16,18 and A=20.
Plot of -v2v(iv) versus ξ according to the FE approximation: this product should return the constant A.
Plot of -v2v(iv) versus ξ according to the pseudo-spectral approximation: this product should return the constant A.
Slope (solid), bending moment (dashed), and shearing force (dotted) near pull-in. Graphs are indistinguishable between the pseudo-spectral and the FE solution.
Line load density qe on the cantilever beam near pull-in according to the FE approximation (A=20).
Line load density qe on the cantilever beam near pull-in according to the FE approximation (A=12).
Line load density qe on the cantilever beam near pull-in according to the pseudo-spectral approximation (A=20).
Line load density qe on the cantilever beam near pull-in according to the pseudo-spectral approximation (A=12).
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.
SalomoniV. A.MajoranaC. E.GiannuzziG. M.MiliozziA.Thermal-fluid flow within innovative heat storage concrete systems for solar power plants2008187-896999910.1108/096155308108990512-s2.0-57149146711NobiliA.LanzoniL.TarantinoA. M.Experimental investigation and monitoring of a polypropylene-based fiber reinforced concrete road pavement20134788889510.1016/j.conbuildmat.2013.05.0772-s2.0-84879215645SalomoniV. A.MajoranaC. E.PomaroB.XottaG.GramegnaF.Macroscale and mesoscale analysis of concrete as a multiphase material for biological shields against nuclear radiation201438551853510.1002/nag.22222-s2.0-84896718922DinelliG.BelzG.MajoranaC. E.SchreflerB. A.Experimental investigation on the use of fly ash for lightweight precast structural elements1996291946326382-s2.0-2142802142BisegnaP.LucianoR.Bounds on the overall properties of composites with debonded frictionless interfaces1998281–4233210.1016/s0167-6636(97)00046-x2-s2.0-0032067994NobiliA.TarantinoA. M.Magnetostriction of a hard ferromagnetic and elastic thin-film structure20081329512310.1177/1081286506073716MR23879812-s2.0-39549087052KudaibergenovA.NobiliA.PrikazchikovaL.On low-frequency vibrations of a composite string with contrast properties for energy scavenging fabric devices201611323124310.2140/jomms.2016.11.n3p3KaplunovJ.NobiliA.Multi-parametric analysis of strongly inhomogeneous periodic waveguideswith internal cutoff frequencies201610.1002/mma.3900ZhangW.-M.YanH.PengZ.-K.MengG.Electrostatic pull-in instability in MEMS/NEMS: a review201421418721810.1016/j.sna.2014.04.0252-s2.0-84901353876ZhangY.ZhaoY.-P.Static study of cantilever beam stiction under electrostatic force influence200417210411210.1007/s10338-004-0413-z2-s2.0-4544304317WeiL. C.MohammadA. B.KassimN. M.Analytical modeling for determination of pull-in voltage for an electrostatic actuated MEMS cantilever beamProceedings of the 5th IEEE International Conference on Semiconductor Electronics (ICSE '02)December 2002Penang, MalaysiaIEEE2332382-s2.0-84864652090ChowdhuryS.AhmadiM.MillerW. C.A closed-form model for the pull-in voltage of electrostatically actuated cantilever beams200515475676310.1088/0960-1317/15/4/0122-s2.0-17444404224SadeghianH.RezazadehG.OsterbergP. M.Application of the generalized differential quadrature method to the study of pull-in phenomena of MEMS switches20071661334134010.1109/JMEMS.2007.9092372-s2.0-36949015590RamezaniA.AlastyA.AkbariJ.Closed-form solutions of the pull-in instability in nano-cantilevers under electrostatic and intermolecular surface forces20074414-154925494110.1016/j.ijsolstr.2006.12.0152-s2.0-34047148779TarantinoA. M.On the finite motions generated by a mode I propagating crack19995728510310.1023/a:1007673212904MR17932922-s2.0-0033289811TarantinoA. M.Nonlinear fracture mechanics for an elastic Bell material199750343545610.1093/qjmam/50.3.435MR1468726LucianoR.WillisJ. R.Hashin-Shtrikman based FE analysis of the elastic behaviour of finite random composite bodies20061371–426127310.1007/s10704-005-3067-z2-s2.0-33644848396NobiliA.RadiE.LanzoniL.A cracked infinite Kirchhoff plate supported by a two-parameter elastic foundation201434112737274410.1016/j.jeurceramsoc.2013.12.0292-s2.0-84901284573LanzoniL.RadiE.Thermally induced deformations in a partially coated elastic layer20094661402141210.1016/j.ijsolstr.2008.11.0042-s2.0-59349105859GuidiV.LanzoniL.MazzolariA.Patterning and modeling of mechanically bent silicon plates deformed through coactive stresses201152031074107910.1016/j.tsf.2011.09.0082-s2.0-81855185303TulliniN.TralliA.LanzoniL.Intefacial shear stress analysis of bar and thin film bonded to 2D elastic substrate using a coupled FE-BIE method201255425110.1016/j.finel.2012.02.006MR29094722-s2.0-84858125732MastrangeloC. H.1999605Cambridge, UKCambridge University PressMRS ProceedingsMerlijn Van SpengenW.PuersR.De WolfI.A physical model to predict stiction in MEMS200212570271310.1088/0960-1317/12/5/3292-s2.0-0036732328YapuZ.Stiction and anti-stiction in MEMS and NEMS200319111010.1007/bf024874482-s2.0-0037872055XottaG.MazzuccoG.SalomoniV. A.MajoranaC. E.WillamK. J.Composite behavior of concrete materials under high temperatures201564869910.1016/j.ijsolstr.2015.03.0162-s2.0-84929504380Marotti de SciarraF.SalernoM.On thermodynamic functions in thermoelasticity without energy dissipation2014468495NobiliA.On the generalization of the Timoshenko beam model based on the micropolar linear theory: static case20152015891435710.1155/2015/914357MR33354832-s2.0-84928242570NapoliG.NobiliA.Mechanically induced Helfrich-Hurault effect in lamellar systems200980303171010.1103/physreve.80.0317102-s2.0-70349901339BoydJ. P.20002ndDover PublicationsGrecoF.LucianoR.A theoretical and numerical stability analysis for composite micro-structures by using homogenization theory201142338240110.1016/j.compositesb.2010.12.0062-s2.0-79951555558PolyaninA. D.ZaitsevV. F.19951stBoca Raton, Fla, USACRC PressSachdevP. L.1991Marcel DekkerMR1109918FunaroD.19928SpringerLecture Notes in PhysicsMR1176949BarrettaR.FeoL.LucianoR.Marotti de SciarraF.Variational formulations for functionally graded nonlocal Bernoulli-Euler nanobeams2015129808910.1016/j.compstruct.2015.03.0332-s2.0-84928102803BarrettaR.FeoL.LucianoR.Marotti de SciarraF.A gradient Eringen model for functionally graded nanorods20151311124113110.1016/j.compstruct.2015.06.0772-s2.0-84938151438BarrettaR.FeoL.LucianoR.Some closed-form solutions of functionally graded beams undergoing nonuniform torsion201512313213610.1016/j.compstruct.2014.12.0272-s2.0-84920997588BarrettaR.LucianoR.Exact solutions of isotropic viscoelastic functionally graded Kirchhoff plates2014118144845410.1016/j.compstruct.2014.07.0442-s2.0-84910667854ApuzzoA.BarrettaR.LucianoR.Some analytical solutions of functionally graded Kirchhoff plates20156826626910.1016/j.compositesb.2014.08.0482-s2.0-84907601496DeziL.MendittoG.TarantinoA. M.Homogeneous structures subjected to repeated structural system changes199011681723173210.1061/(ASCE)0733-9399(1990)116:8(1723)2-s2.0-0025466342DeziL.TarantinoA. M.Time-dependent analysis of concrete structures with a variable structural system19918833203242-s2.0-0026151680DeziL.MendittoG.TarantinoA. M.Viscoelastic heterogeneous structures with variable structural system1993119223825010.1061/(ASCE)0733-9399(1993)119:2(238)2-s2.0-0027543374BarrettaR.FeoL.LucianoR.Torsion of functionally graded nonlocal viscoelastic circular nanobeams20157221722210.1016/j.compositesb.2014.12.0182-s2.0-84921726407