The dynamic problems of a microactuator with a single edge crack are numerically formulated using radial basis functions. The microactuator model incorporates the taper ratio, electrode shapes, and crack length, all of which govern the dynamic behavior of microactuators. To optimize the design of a microactuator, many characteristics of various shaped cantilevers and curved electrodes are also investigated.
1. Introduction
Microelectromechanical systems exploit microscale effects to extend the range applications of actuators, accelerometers, angular rate sensors, and other devices. Elucidation of the dynamic mechanism of electrostatic microactuators contributes markedly to their design. Legtenberg et al. [1] investigated the dynamic behavior of active joints for different electrostatic actuator designs and proposed the idea of using a curved electrode to improve pull-in performance. Electrostatic actuators are widely applied in microelectromechanical systems. Electrostatic microactuator devices have a high operating frequency and low power consumption. Hong et al. [2] studied the influence of the dimensions and stress of such a device on fatigue endurance when an external force was applied to a normal microcantilever beam and a notch cantilever beam. They performed analysis simulation indicating that the stress was maximal at the fixed end. Their results showed that a deep notch in a specimen concentrate stress and thus promote specimen failure. Mehdaoui et al. [3] presented the vertical cointegration of AlSi MEMS tunable capacitors and Cu inductors in tunable LC blocks. Etxeberria and Gracia [4] proposed tunable MEMS volume capacitors for high-voltage applications. Liu et al. [5] presented actuation by electrostatic repulsion established by nonvolatile charge injection. Gallant and Wood [6] investigated how fabrication techniques affect the performance of widely tunable micromachined capacitors. Borwick III et al. [7] analyzed a high Q microelectromechanical capacitor with large tuning range for RF filter systems. Harsh et al. [8] studied the design and realization of a flip-chip integrated microelectromechanical tunable capacitor. Osterberg et al. [9, 10] proposed a one-dimensional model and a three-dimensional model for analyzing electrostatically deformed diaphragms. Their results revealed that the electrostatic deformation calculated using the one-dimensional model is close to that obtained using a three-dimensional model. Gilbert et al. [11] analyzed the three-dimensional coupled electromechanics of microelectromechanical systems using a CoSolve-EM simulation algorithm. Elwenspoek et al. [12] studied the dynamic behavior of active joints for various electrostatic actuator designs. Shi et al. [13] presented the combination of an exterior boundary element method for analyzing electrostatics and a finite-element method for analyzing elasticity to evaluate the effect of coupling between the electrostatic force and the elastic deformation. Gretillat et al. [14] employed the three-dimensional MEMCAD and finite-element method programs to simulate the dynamics of a nonlinear actuator, taking into account squeeze-film damping. Hung and Senturia [15] developed leveraged bending and strain-stiffening methods to increase the maximum travel distance before the pull-in of electrostatic actuators. Chan et al. [16] measured the pull-in voltage and capacitance-voltage characteristic and performed two-dimensional simulations that included the electrical effects of fringing fields and finite-beam thickness to determine the material properties of electrostatic microactuators. Li and Aluru [17] developed a mixed-regime approach for combining linear and nonlinear theories to analyze large microelectromechanical structure deformations at large applied voltages. Their results demonstrated that electrostatic actuators can undergo large deformation at certain driving voltages. Chyuan et al. [18–20] established the validity and accuracy of the dual boundary element method and employed it to elucidate the effect of a variation in gap size on the levitation of a microelectromechanical comb drive. Qiao et al. [21] presented a suspension beam called a two beam to realize a parallel-plate actuator with an extended working range, but without the disadvantages of complex control circuit and high actuation voltage. In this investigation, radial basis functions are adopted to analyze how cantilever shape, damping, cracks, and electrode shape affect dynamic behavior in electrostatic actuator systems. The radial basis function scheme is applied to formulate the electrostatic field problems in matrix form. The integrity and computational accuracy of radial basis functions are demonstrated with reference to various case studies. To the author’s knowledge, very few published investigates have presented a vibration analysis of a cantilever electrostatic microactuator with an edge crack using radial basis functions.
2. Radial Basis Function
A radial basis function is a real-valued function whose value depends on the distance from an origin. Kansa [22, 23] studied a given function or partial derivatives of a function with respect to a coordinate direction which is expressed as a linear weighted sum of all functional values at all mesh points along the direction that was initiated based on the concept of radial basis function. In their algorithm, the node distribution was completely unstructured. Wang and Liu [24] proposed a point interpolation meshless method that was based on radial basis functions and incorporated the Galerkin weak form for solving partial differential equations. Elfelsoufi and Azrar [25] investigated the buckling, flutter, and vibration of beams using radial basis functions. Hon et al. [26] used radial basis functions for function fitting and solving partial differential equations using global nodes and collocation procedures. Liu et al. [27] constructed shape functions with the delta function property based on radial and polynomial basis functions. In this study, shape functions are constructed using radial basis functions. A radial basis function can be expressed as follows [28, 29]:
(2.1)Bi(x)=(x-xi)2+c2,
where c is a constant. The radial basis function is typically used to develop the functional approximations of the following form [28, 29]:
(2.2)v(x,t)=∑i=1Nai(t)Bi(x),
where ai is the coefficient to be determined. The microcantilever deflection v(x,t) denotes a sum of N radial basis functions, each associated with a different center xi. The domain contains N collocation points. Although this nonlinear equation of the electrostatic microactuator does not have an analytical solution, numerical approaches can be adopted to solve it. These nonlinear partial differential equations are obtained numerically using the radial basis function approach, which does not require a mesh.
3. Dynamic Behavior of Clamped-Free Microactuators
Figure 1 displays the geometry of an electrostatic actuator with an edge crack close to its fixed end. Variable t0 is the thickness of a microactuator at x=0, t1 denotes the tip thickness of a microactuator at x=L, and L is the length of the microbeam. An electrostatic force, introduced by the difference between the driving voltage of the curved electrode and that of the cantilever, pulls the cantilever toward the curved electrode. The electrostatic force is approximately proportional to the inverse of the square of the distance between the curved electrode and the shaped cantilever. The equation of motion of an electrostatic microactuator with an edge crack near the fixed end can be derived as [1, 30]
(3.1)ρA∂2v∂t2+D∂2∂x2(EI∂3v∂x2∂t)+∂2∂x2(EI∂2v∂x2)=ε0bV-2(1+0.65((d+S-αxt0/2L-v)/b))2(d+S-αxt0/2L-v)2,
where S represents the shape of the curved electrode, and is given by a polynomial such that S=(xδ/L)n, δ is the gap between the tip of a curved electrode at x=L and the tip of a micro cantilever at x=L, α is defined as the ratio (t1-t0)/t0, and n is the polynomial order of the shape of the electrode. The electrode shape varies with the value of n. V- is the driving voltage; E is Young’s modulus of the actuator material; ε0 is the dielectric constant of air, ε0=8.85×10-12; b is the width of the microactuator; d is the initial gap, as displayed in Figure 1. The dielectric layer prevents short circuits. The cross-sectional area of a microactuator is A(x)=bt0(1+αx/L), and I(x) is the moment of inertia of the cross-sectional area of a microactuator and is given by I(x)=I0(1+αx/L)3 and I0=bt03/12. The Kelvin-Voigt damping force D(∂2/∂x2)(EI(∂3v/∂t∂x2)) is assumed to model resistance to the actuator strain velocity. D is the Kelvin-Voigt damping coefficient. To ensure generality, Kelvin-Voigt damping effects are considered in the formulation of the equations of motion [31, 32]. Equation (3.1) depicts the fringing effects of the electrical field. The dynamic characteristics of edge-cracked beams are of considerable importance in many designs. The flexibility G caused by a crack of depth a- can be determined using Broek’s approximation [33] to be
(3.2)(1-μ2)KI2E=(Pb)22bdGda-,
where KI is the stress intensity factor under mode I loading; μ is Poisson’s ratio; Pb denotes the bending moment at the crack; G is the flexibility of the micro cantilever. The magnitude of the stress intensity factor can be determined using Tada’s formula [34], as(3.3)KI=6Pbbt02πrt0FI(r),
where(3.4)FI(r)=2πrtan(πr2)0.923+0.199(1-sin(πr/2))4cos(πr/2),r=a-t0.
Substituting the stress intensity factor KI into (3.2), yields(3.5)G=6(1-μ2)t0∫0rπrFI2(r)drEI0.
Since(3.6)kT=1G,
the bending stiffness kT of the cracked section of a micro cantilever can be expressed as (3.7)kT=EI06(1-μ2)t0∫0rπrFI2(r)dr.
A crack can be represented as a spring of zero length and zero mass. The boundary conditions associated with (3.1) for a microactuator with an edge crack near its fixed end are given by
(3.8)v(0,t)=0,EI∂2v(0,t)∂x2=kT∂v(0,t)∂x,EI∂v2(L,t)∂x2=0,∂∂x(EI∂v2(L,t)∂x2)=0.
This nonlinear equation does not have an analytical solution; however, numerical approaches can be utilized to solve it. The radial basis function approach is adopted to solve numerically these nonlinear partial differential equations. In the radial basis function approach, (2.2) is substituted into (3.1). The equation of motion of a fixed-free microbeam can be rearranged into a formula based on the radial basis function approach: (3.9)[∂2EI(xi)∂x2∂2B1(xi)∂x2∂2EI(xi)∂x2∂2B2(xi)∂x2⋯∂2EI(xi)∂x2∂2BN(xi)∂x2][a1a2⋯aN]T+[2∂EI(xi)∂x∂3B1(xi)∂x32∂EI(xi)∂x∂3B2(xi)∂x3⋯2∂EI(xi)∂x∂3BN(xi)∂x3][a1a2⋯aN]T+[EI(xi)∂4B1(xi)∂x4EI(xi)∂4B2(xi)∂x4⋯EI(xi)∂4BN(xi)∂x4][a1a2⋯aN]T+[D∂2EI(xi)∂x2∂2B1(xi)∂x2D∂2EI(xi)∂x2∂2B2(xi)∂x2⋯D∂2EI(xi)∂x2∂2BN(xi)∂x2][∂a1∂t∂a2∂t⋯∂aN∂t]T+[2D∂EI(xi)∂x∂3B1(xi)∂x32D∂EI(xi)∂x∂3B2(xi)∂x3⋯2D∂EI(xi)∂x∂3BN(xi)∂x3][∂a1∂t∂a2∂t⋯∂aN∂t]T+[DEI(xi)∂4B1(xi)∂x4DEI(xi)∂4B2(xi)∂x4⋯DEI(xi)∂4BN(xi)∂x4][∂a1∂t∂a2∂t⋯∂aN∂t]T+[ρA(xi)B1(xi)ρA(xi)B2(xi)⋯ρA(xi)BN(xi)][∂2a1∂t2∂2a2∂t2⋯∂2aN∂t2]T=[ε0bV-2(1+0.65((d+S-αxt0/2L-∑j=1NajBj(xi))/b))B1(xi)2(d+S-αxt0/2L-∑j=1NajBj(xi))2ε0bV-2(1+0.65((d+S-αxt0/2L-∑j=1NajBj(xi))/b))B2(xi)2(d+S-αxt0/2L-∑j=1NajBj(xi))2⋯ε0bV-2(1+0.65((d+S-αxt0/2L-∑j=1NajBj(xi))/b))BN(xi)2(d+S-αxt0/2L-∑j=1NajBj(xi))2]×[a1a2⋯aN]Tfori=1,2,…,N.
Based on the radial basis function approach, the boundary conditions of a clamped-free microactuator with an edge crack can be rearranged into matrix forms as
(3.10)[B1(x1)B2(x1)⋯BN(x1)][a1a2⋯aN]T=[0],[EI(x1)∂2B1(x1)∂x2EI(x1)∂2B2(x1)∂x2⋯EI(x1)∂2BN(x1)∂x2][a1a2⋯aN]T=[kT∂B1(x1)∂xkT∂B2(x1)∂x⋯kT∂BN(x1)∂x][a1a2⋯aN]T,[EI(xN)∂2B1(xN)∂x2EI(xN)∂2B2(xN)∂x2⋯EI(xN)∂2BN(xN)∂x2][a1a2⋯aN]T=[0],[EI(xN)∂3B1(xN)∂x3EI(xN)∂3B2(xN)∂x3⋯EI(xN)∂3BN(xN)∂x3][a1a2⋯aN]T+[∂EI(xN)∂x∂2B1(xN)∂x2∂EI(xN)∂x∂2B2(xN)∂x2⋯∂EI(xN)∂x∂2BN(xN)∂x2]×[a1a2⋯aN]T=[0].
Schematic of a clamped-free curved electrode microactuator with an edge crack.
4. Numerical Results
The following figures summarize the results thus obtained. Figure 2 shows the frequencies of a clamped-free curved electrode microactuator. The material and geometric parameters of the actuator considered herein are E=150GPa, δ=30μm, b=5μm, t0=2μm, r=0, α=0, and d=2μm [1, 35]. The figure plots the analytical solutions and the numerical results obtained using radial basis function approach. Numerical results indicate that the estimated frequencies remain stable even when only fifteen collocation points are considered. They also suggest that the frequencies calculated using the radial basis function approach are extremely close to the exact solutions. Figure 3 shows the frequencies of a clamped-free curved electrode microactuator for various crack depths. The frequency falls, as the crack depth increases. The crack depth significantly affects the frequencies of the micro cantilever. Computational results solved using radial basis function approach are compared with numerical results obtained using the differential quadrature method. Figure 4 compares the tip deflections of an actuating electrode for various driving voltages and electrode shapes. Changing the electrode shape in an electrostatic microactuator is an effective technique for varying the electrostatic force distribution therein. Numerical and measured results suggest that the tip deflections calculated using the radial basis function approach are in good agreement with published experimental results [1]. Numerical results demonstrate that the pull-in voltage declines gradually, as the value of n increases. Figure 5 plots the tip responses of the microactuator with different α values. The tip deflections of the micro cantilever drop, as α increases. As α increases, the applied voltage required to cause a particular deflection of the tip of the micro cantilever increases. The cantilever shapes substantially influence the pull-in behavior of microactuators. Figure 6 plots the tip responses of the microactuator for various values of r. Notably, the tip deflection increases with crack depth. The depth of the crack significantly affects the tip response. Figure 7 shows the variation of the tip responses of the actuator with D. Because of recent advances in stably responding and high performance actuator structures, the enhancement of damping has become a very significant issue. The numerical results in this example show that internal damping can significantly affect the dynamic behavior of the actuator system. Strong residual vibration occurs in a system with a zero internal damping coefficient.
Frequencies of the clamped-free curved electrode microactuator obtained using radial basis function approach at various collocation points.
Frequencies of clamped-free curved electrode microactuators for various values of r.
Comparison of tip deflections of variously shaped actuating electrodes.
Tip responses of microactuators for different values of α.
Tip responses of microactuators for various values of r.
Tip responses of microactuators for various values of D.
5. Conclusions
This work examines the radial basis functions for dynamic problems of an electrostatic actuator with a crack. The effects of internal damping, electrode shape, edge cracking, and cantilever shape on the pull-in behavior of electrostatic microstructures are investigated. The frequency of the microcantilever declines, as the crack depth increases. The value of radial basis functions that describe the pull-in behavior of a microactuator with an edge crack is determined.
LegtenbergR.GilbertJ.SenturiaS. D.ElwenspoekM.Electrostatic curved electrode actuators1997632572652-s2.0-003123512010.1109/84.623115HongH.HungJ. N.GuuY. H.Various fatigue testing of polycrystalline silicon microcantilever beam in bending2008476525652612-s2.0-5504909963910.1143/JJAP.47.5256MehdaouiA.PisaniM. B.FritschiR.AnceyP.IonescuA. M.Vertical co-integration of AlSi MEMS tunable capacitors and Cu inductors for tunable LC blocks2007845–8136913732-s2.0-3424789713810.1016/j.mee.2007.01.222EtxeberriaJ. A.GraciaF. J.Tunable MEMS volume capacitors for high voltage applications2007845–8139313972-s2.0-3424764441510.1016/j.mee.2007.01.063LiuZ.KimM.ShenN. Y. M.KanE. C.Actuation by electrostatic repulsion by nonvolatile charge injection200511912362442-s2.0-1714442499710.1016/j.sna.2004.03.036GallantA. J.WoodD.The role of fabrication techniques on the performance of widely tunable micromachined capacitors20041101–34234312-s2.0-164256650410.1016/j.sna.2003.07.006BorwickR. L.IIIStuparP. A.DeNataleJ.AndersonR.TsaiC.GarrettK.ErlandsonR.A high Q, large tuning range MEMS capacitor for RF filter systems20031031-233412-s2.0-003743902110.1016/S0924-4247(02)00316-3HarshK. F.SuB.ZhangW.BrightV. M.LeeY. C.Realization and design considerations of a flip-chip integrated MEMS tunable capacitor20008021081182-s2.0-003390690710.1016/S0924-4247(99)00255-1OsterbergP.YieH.CaiX.WhiteJ.SenturiaS.Self-consistent simulation and modeling of electrostatically deformed diaphragms199428322-s2.0-0027989013OsterbergP. M.SenturiaS. D.M-test: a test chip for MEMS material property measurement using electrostatically actuated test structures1997621071182-s2.0-003116899010.1109/84.585788GilbertJ. R.LegtenbergR.SenturiaS. D.3D coupled electro-mechanics for MEMS: applications of CoSolve-EMProceedings of the IEEE Micro Electro Mechanical Systems ConferenceFebruary 19951221272-s2.0-0029194943ElwenspoekM.WeustinkM.LegtenbergR.Static and dynamic properties of active jointsProceedings of the 8th International Conference on Solid-State Sensors and Actuators and EurosensorsJune 19954124152-s2.0-0029519371ShiF.RameshP.MukherjeeS.Simulation methods for micro-electro-mechanical structures (MEMS) with application to a microtweezer19955657697832-s2.0-002963471410.1016/0045-7949(95)00007-4GretillatM. A.YangY. J.HungE. S.RabinovichV.AnanthasureshG. K.de RooijN. F.SenturiaS. D.Nonlinear electromechanical behavior of an electrostatic microrelayProceedings of the International Conference on Solid-State Sensors and ActuatorsJune 1997114111442-s2.0-0030655325HungE. S.SenturiaS. D.Extending the travel range of analog-tuned electrostatic actuators1999844975052-s2.0-003331619410.1109/84.809065ChanE. K.GarikipatiK.DuttonR. W.Characterization of contact electromechanics through capacitance-voltage measurements and simulations1999822082172-s2.0-003314867410.1109/84.767117LiG.AluruN. R.Linear, nonlinear and mixed-regime analysis of electrostatic MEMS20019032782912-s2.0-0035880218ChyuanS. W.LiaoY. S.ChenJ. T.An efficient method for solving electrostatic problems20035352582-s2.0-003754551910.1109/MCISE.2003.1196307ChyuanS. W.LiaoY. S.ChenJ. T.Computational study of variations in gap size for the electrostatic levitating force of MEMS device using dual BEM20043597397482-s2.0-394306633810.1016/j.mejo.2004.06.005LiaoY. S.ChyuanS. W.ChenJ. T.Efficaciously modeling the exterior electrostatic problems with singularity for electron devices2004205253410.1109/MCD.2004.1343246QiaoD. Y.YuanW. Z.LiX. Y.A two-beam method for extending the working range of electrostatic parallel-plate micro-actuators20076542562622-s2.0-3384612939710.1016/j.elstat.2006.09.001KansaE. J.Multiquadrics—a scattered data approximation scheme with applications to computational fluid-dynamics-I surface approximations and partial derivative estimates1990198-91271452-s2.0-0025229330ZBL0692.76003KansaE. J.Multiquadrics—a scattered data approximation scheme with applications to computational fluid-dynamics-II solutions to parabolic, hyperbolic and elliptic partial differential equations1990198-91471612-s2.0-0025210711ZBL0850.76048WangJ. G.LiuG. R.A point interpolation meshless method based on radial basis functions20025411162316482-s2.0-003714377510.1002/nme.489ZBL1098.74741ElfelsoufiZ.AzrarL.Buckling, flutter and vibration analyses of beams by integral equation formulations20058331-32263226492-s2.0-2734445320310.1016/j.compstruc.2005.04.001HonY. C.LuM. W.XueW. M.ZhuY. M.Multiquadric method for the numerical solution of a biphasic mixture model1997882-31531752-s2.0-000102385710.1016/S0096-3003(96)00309-8ZBL0910.76059LiuG. R.ZhaoX.DaiK. Y.ZhongZ. H.LiG. Y.HanX.Static and free vibration analysis of laminated composite plates using the conforming radial point interpolation method20086823543662-s2.0-3744903487910.1016/j.compscitech.2007.07.014WrightG. B.2003University of Colorado at BoulderBuhmannM. D.2003New York, NY, USACambridge University PressChenK. S.OuK. S.Development and verification of 2D dynamic electromechanical coupling solver for micro-electrostatic-actuator applications200713614034112-s2.0-3424713622710.1016/j.sna.2006.10.042CloughR. W.PenzienJ.1975New York, NY, USAMcGraw-HillRaoS. S.1990New York, NY, USAAddison-WesleyBroekD.1986Leiden, The NetherlandsMartinus NijhoffTadaH.ParisP. C.IrwinG. R.2000Professional EngineeringLiuY.LiewK. M.HonY. C.ZhangX.Numerical simulation and analysis of an electroactuated beam using a radial basis function2005146116311712-s2.0-2834444093710.1088/0964-1726/14/6/009