Published: 15 August 2015

Periodic solution to a nonlinear oscillator arising in micro electro mechanical system

A. Nikkar1
H. Jodeiri2
Sh. Ahari Abesh Ahmad3
H. Tariverdian4
1Young Researchers and Elite Club, Ahar Branch, Islamic Azad University, Ahar, Iran
2Young Researchers and Elite Club, Tabriz Branch, Islamic Azad University, Tabriz, Iran
3Faculty of Civil Engineering, University of Tabriz, Tabriz, Iran
4Department of Civil Engineering, Sardroud Branch, Islamic Azad University, Sardroud, Iran
Corresponding Author:
H. Jodeiri
Views 33
Reads 16
Downloads 1027

Abstract

In this paper, a periodic solution of nonlinear oscillator arising in the micro-beam based microelectromechanical system (MEMS) has been analytically achieved. The Amplitude Frequency Formulation (AFF) and Max Min Approach (MMA) have been used. In the second method (which is called MMA), an approximate solution of the nonlinear equation can be easily deduced by finding Maximal and minimal solution thresholds of this nonlinear problem. What we understood is that both methods, works properly and scales down the deal of the work. Compare conclusions with the results from fourth order Runge-Kutta method and energy balance method (EBM) shows that obtained results are of high accuracy and convenient.

1. Introduction

Study micro-electro-mechanical systems (MEMS) are batch-fabricated devices and structures at a microscale level [1] which can be defined as a system of small dimensions (less than one cubic centimeter) fulfilling a smart function. Since its inception, MEMS technologies are of tremendous importance in various engineering fields. Because of its small size, low power consumption and high reliability, we have seen many potential applications of MEMS devices which require few mechanical components and small voltage levels for actuation actuators and sensors in aerospace, optical and biomedical engineering [2-4]. MEMS devices are generally classified according to their actuation mechanisms. Actuation mechanisms for MEMS vary depending on the suitability to the application at hand. The most common actuation mechanisms are electrostatic, pneumatic, thermal, and piezoelectric [5]. With reported successes in applications ranging from signal filtering to chemical and mass sensing, electrostatically actuated resonant microbeam devices which form a broad class of MEMS devices due to their simplicity, as they require few mechanical components and small voltage levels for actuation [5] continue to be one of the most widely studied topics of the MEMS community.

Almost every natural event can be modeled with nonlinear equations; this is due to nature of the phenomena in the world. Although it is easy to find the solutions of some problems by using of computers, it is still difficult to solve nonlinear equations either numerically or analytically. Also it is very difficult to obtain the exact answer for these equations. Nonlinearity in MEMS which may cause some difficulties in computations, is made by Electrostatic actuation, large deflections and damping. According to these difficulties, some analytic techniques are available for nonlinear problems, such as MEMS [5-15]. Amplitude Frequency Formulation and Max Min approach (MMA), suggested by J. H. He. [16, 17]. Both methods are simple and effective methods to solve nonlinear oscillatory equations. In this paper, AFF and MMA are used to study a nonlinear oscillator arising in the micro-beam based MEMS where the mid plane stretching effect and distributed electrostatic force are both considered.

2. Problem description

Consider a clamped-clamped microbeam-based electromechanical resonator with uniform thickness h, length l, width b (b5h), effective modulus E-=E/(1-ν)2, Young’s modulus E, Poisson’s ratio v and density q, as shown in Fig. 1 together with coordinate oxyz, where g0 is the initial gap and V the electrostatic load. The equation of motion that governs the transverse deflection is written as:

1
E-I4wx4+ρS2wt2-N-+ES2l0lwx22wx2=qx,t,

where S=bh and N- denotes the tensile or compressive axial load created by the mismatch of both thermal expansion coefficient and crystal lattice period between substrate and the thin film (microbeam). The second term on the right-hand side, q(x,t), denotes the driving force per unit length, resulting from electrostatic excitation [18]. Using the classical beam theory and taking into account of the mid-plane stretching effect as well as the distributed electrostatic force, the following dimensionless equation of motion for the microbeam can be formulated via the Galerkin method [18]:

2
u¨(a1u4+a2u2+a3)+a4u+a5u3+a6u5+a7u7=0,

where:

3
a1=01ϕ6dζ,
4
a2=-201ϕ4dζ,
5
a3=01ϕ2dζ,
6
a4=01(ϕ''''ϕ-Nϕ''ϕ-V2ϕ)dζ,
7
a5=-01(2ϕ''''ϕ3-2Nϕ''ϕ3+αϕ''ϕ01(ϕ')2dζ)dζ,
8
a6=01(ϕ''''ϕ5-Nϕ''ϕ5+2αϕ''ϕ301(ϕ')2dζ)dζ,
9
a7=-01αϕ''ϕ501(ϕ')2dζdζ,

and

10
ϕ=16ξ21-ξ2.

Here, an overdot () represents differentiation with respect to the time variable, while a prime (') demonstrates the partial differentiation with respect to the coordinate variable ξ.

3. Basic idea of amplitude frequency formulation (AFF)

To illustrate the basic idea of AFF method, the following nonlinear oscillator is considered:

11
u¨+Nu¨,u˙,u,t=0, u0=A, u˙0=0.

For a generalized nonlinear oscillator in Eq. (11), we use two trial functions:

12
u1=Acosω1t, ω1=1,
13
u2=Acosω2t, ω2=ω,

where ω is assumed to be the frequency of the nonlinear oscillator, Eq. (11). The residuals are:

14
R1t=-cost+fAcost,
15
R2t=-ω2cosωt+fAcosωt.

Fig. 1Schematics of a fixed microbeam-based electromechanical resonator [18]

Schematics of a fixed microbeam-based electromechanical resonator [18]

We introduce R11 and R22 defined as:

16
R11t=4T10T14R1tcostdt, T1=2π,
17
R22t=4T20T24R2tcosωtdt, T2=2πω.

Applying He’s frequency-amplitude formulation, we have:

18
ω2=ω12R22(t)-ω22R11(t)R22(t)-R11(t),

where ω1= 1 and ω2=ω.

4. Basic idea of max-min approach

In this section, MMA is applied to discuss the problem. To illustrate the basic idea of MMA method, the following nonlinear oscillator is considered:

19
u¨+Nu¨,u˙,u,t=0, u0=A, u˙0=0.

In order to the fact that small parameters or linear terms are not the requirements of MMA, Eq. (18) can be approximately solved by using the MMA. Considering a, b, c and d as the real numbers:

20
ab<x<dc.

Then:

21
ab<ma+ndmb+nc<dc,

where m and n are weighting factors and x is a rough approximation of:

22
x=ma+ndmb+nc.

Eq. (19) can be rewritten in the following form:

23
u¨+ufu¨,u˙,u,t=0.

And the frequency can be identified as follows:

24
ab<ϖ2=ma+ndmb+nc<dc.

Then:

25
u¨+ϖ2u=u¨+Nu¨,u˙,u,t+ρu¨,u˙,u,t,
26
ρu¨,u˙,u,t=0,

here ϖ can be obtained by substitution of Acosϖt as initial assumption into Eq. (25).

5. Implementation of AFF

In this section, in order to assess the accuracy of He’s amplitude-frequency formulation for solving nonlinear governing equation of motion and to compare it with the other solutions, in this section we consider this method. According to He’s amplitude-frequency formulation, we choose two trial functions u1(t)=Acost and u2(t)=Acosωt, where ω is assumed to be the frequency of the nonlinear oscillator Eq. (2). Substituting the trial functions into Eq. (10) results in, respectively, the following residuals:

27
R1(t)=-(a1A4cos(t)4+a2A2cos(t)2+a3)×Acos(t)+a4Acos(t)+a5A3cos(t)3
+a6A5cos(t)5+a7A7cos(t)7,
28
R2(t)=-(a1A4cos(ωt)4+a2A2cos(ωt)2+a3)×Acos(ωt)2ω2+a4Acos(ωt)
+a5A3cos(ωt)3+a6A5cos(ωt)5+a7A7cos(ωt)7.

We introduce R11 and R22 defined as:

29
R11(t)=4T10T14R1(t)cos(t)dt=1128A×(40a6A4-40a1A4-48a2A2+48a5A4
+35a7A6+64a4-64a3), T1=2π,
30
R22t=4T20T24R2tcosωtdt=1128A×(48a5A2+40a6A4+35a7A6+64a4
-64ω2a3-40A4ω2a1-48A2ω2a2), T2=2πω.

Applying He’s frequency-amplitude formulation, we have:

31
ω2=ω12R22(t)-ω22R11(t)R22(t)-R11(t),

where ω1= 1 and ω2=ω. Finally, the amplitude-depended frequency can be approximated as:

32
ωAFF=2464a4+48a5A2+40a6A4+35a7A65a1A4+6a2A2+8a3.

Its approximate solution reads:

33
ut=Acos2464a4+48a5A2+40a6A4+35a7A65a1A4+6a2A2+8a3t.

6. Application of MMA

In this section, the MMA is applied to solve nonlinear sector. Eq. (2) can be rewritten in the following form:

34
u¨+a1a3u3u¨+a2a3uu¨+a4a3+a5a3u2+a6a3u4+a7a3u6u=0.

In addition, boundary conditions are considered. Here, if ut=Acos(ϖt) be chosen as a trial function, then:

35
a4a3<a1a3u3u¨+a2a3uu¨+a4a3+a5a3u2+a6a3u4+a7a3u6
<-a1A4ϖ2-a2A2ϖ2+a4+a5A2+a6A4+a7A6a3,

where a4 is the minimum and -a1A4ϖ2-a2A2ϖ2+a4+a5A2+a6A4+a7A6 is the maximum value of a1a3u3u¨+a2a3uu¨+a4a3+a5a3u2+a6a3u4+a7a3u6.

Here the frequency value can be identified as follows:

36
ϖ2=nm+n-a1A4ϖ2-a2A2ϖ2+a4+a5A2+a6A4+a7A6a3+mm+na4a3
=k(-a1A4ϖ2-a2A2ϖ2+a5A2+a6A4+a7A6)a3+a4a3,

where m and n are weighting factors, and k=n/(m+n). So the frequency can be approximated as follows:

37
ϖ=k-a1A4ϖ2-a2A2ϖ2+a5A2+a6A4+a7A6a3+a4a3.

Then Eq. (24) can be rewritten as follows:

38
u¨+ϖ2u=u¨+a1a3u3u¨+a2a3uu¨+a4a3+a5a3u2+a6a3u4+a7a3u6u+ρ,

where:

39
ρ=a3a1a3u3u¨+a2a3uu¨+a4a3+a5a3u2+a6a3u4+a7a3u6u
+-ka1A4ϖ2-ka2A2ϖ2+a4+ka5A2+ka7A6u.

Combining u(t)=Acosϖt with ρ, and using Fourier expansion series gives:

40
ρϖt=n=0δ2n+1cos2n+1ϖt=δ1cosϖt+δ3cos3ϖt+...δ1cosϖt,
41
δ1=4π0π2ρϕcosϕdϕ=8315π(-105ka1A4ϖ2-105ka2A2ϖ2+210a4
+105ka5A2+105ka6A4+105ka7A6-72a1A4ϖ2-84a2A2ϖ2+84a5A2
+72a6A4+64a7A6).

In order to avoid secular term, the following equation can be set:

42
δ1=0.

Then the following equation can be obtained:

43
k=2105-36a1A4ϖ2-42a2A2ϖ+105a4+42a5A2+36a6A4+32a7A6A2-A2a1ϖ2-a2ϖ2+a5+A2a6+A4a7.

Substituting Eq. (42) into Eq. (36), yields:

44
ϖMMA=3384a5A2+64a7A6+72a6A4+105a424a1A4+35a3+28a2A2.

The answers of AFF method and max-min approach are plotted in Fig. 2. For comparison, results of energy balance [18] and numerical fourth-order Runge-Kutta method are plotted in these figures too. Herein the values of parameters are taken as N= 15, α= 10.

Fig. 2The results of MMA, AFF, EBM and Runge-Kutta method for A= 0.2 and 0.4, N= 15, α= 10

The results of MMA, AFF, EBM and Runge-Kutta method for A= 0.2 and 0.4, N= 15, α= 10
The results of MMA, AFF, EBM and Runge-Kutta method for A= 0.2 and 0.4, N= 15, α= 10

7. Conclusion

In this paper, the main purpose was to illustrate the application of amplitude frequency formulation (AFF) and max-min approach (MMA) in solving nonlinear oscillator arising in the microbeam based MEMS. Also, the capabilities and facile applications of these methods have been demonstrated in comparison with the numerical solution and energy balance method. By intensifying the value of voltage and initial condition A, it has been considered that the disagreement between the answers obtained by AFF, MMA and Runge-Kutta is very insignificant which indicates that these methods provide highly precise answers for nonlinear equations.

References

  • Lee S., Ramadoss R., Buck M., Bright V. M., Gupta K. C., Lee Y. C. Reliability testing of flexible printed circuit-based RF MEMS capacitive switches. Microelectronics Reliability, Vol. 44, Issue 2, 2004, p. 245-250.
  • Rhoads J. F., Shaw S. W., Turner K. L. The nonlinear response of resonant microbeam systems with purely-parametric electrostatic actuation. Journal of Micromechanics and Microengineering, Vol. 16, Issue 5, 2006, p. 890-899.
  • Kuehne I., Frey A., Marinkovic D., Eckstein G., Seidel H. Power MEMS – a capacitive vibration-to-electrical energy converter with built-in voltage. Sensors and Actuators A: Physical, Vol. 142, 2008, p. 263-269.
  • Czarnecki P., Rottenberg X., et al. Effect of substrate charging on the reliability of capacitive RF MEMS switches. Sensors and Actuators A: Physical, Vol. 154, 2009, p. 261-268.
  • Tocchio A., Langfelder G., Longoni A., Lasalandra E. In-plane and out-of-plane MEMS motion sensors based on fringe capacitances. Procedia Engineering, Vol. 5, 2010, p. 1392-1395.
  • Taghipor R., Akhlaghi T., Nikkar A. Explicit solution of the large amplitude transverse vibrations of a flexible string under constant tension. Latin American Journal of Solids and Structures, Vol. 11, 2014, p. 545-555.
  • Reza A., Sedighi H. M. Nonlinear vertical vibration of tension leg platforms with homotopy analysis method. Advances in Applied Mathematics and Mechanics, Vol. 7, 2015, p. 357-368.
  • Ghaffarzadeh H., Nikkar A. Explicit solution to the large deformation of a cantilever beam under point load at the free tip using the variational iteration method-II. Journal of Mechanical Science and Technology, Vol. 27, Issue 11, 2013, p. 3433-3438.
  • Nikkar A., Bagheri S., Saravi M. Dynamic model of large amplitude vibration of a uniform cantilever beam carrying an intermediate lumped mass and rotary inertia. Latin American Journal of Solids and Structures, Vol. 11, 2014, p. 320-329.
  • Sedighi H. M., Daneshmand F. Static and dynamic pull-in instability of multi-walled carbon nanotube probes by He’s iteration perturbation method. Journal of Mechanical Science and Technology, Vol. 28, Issue 9, 2013, p. 3459-3469.
  • Bagheri S., Nikkar A., Ghaffarzadeh H. Study of nonlinear vibration of Euler-Bernoulli beams by using analytical approximate techniques. Latin American Journal of Solids and Structures, Vol. 11, 2014, p. 157-168.
  • Bagheri S., Nikkar A. Higher order explicit solutions for nonlinear dynamic model of column buckling using variational approach and variational iteration algorithm-II. Journal of Mechanical Science and Technology, Vol. 28, Issue 11, 2014, p. 4605-4611.
  • Khan Y., Taghipour R., Falahian M., Nikkar A. A new approach to modified regularized long wave equation. Neural Computing and Applications, Vol. 23, 2013, p. 1335-1341.
  • Sedighi H. M., Daneshmand F. Nonlinear transversely vibrating beams by the homotopy perturbation method with an auxiliary term. Journal of Applied and Computational Mechanics, Vol. 1, Issue 1, 2014, p. 1-9.
  • Sedighi H. M., Daneshmand F., Abadyan M. Dynamic instability analysis of electrostatic functionally graded doubly-clamped nano-actuators. Composite Structures, Vol. 124, 2015, p. 55-64.
  • He J. H. Max-min approach to nonlinear oscillators. International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 9, Issue 2, 2008, p. 207-210.
  • He J. H. An improved amplitude-frequency formulation for nonlinear oscillators. International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 9, 2008, p. 211-212.
  • Fu Y., Zhang J., Wan L. Application of the energy balance method to a nonlinear oscillator arising in the microelectromechanical system (MEMS). Current Applied Physics, Vol. 11, 2011, p. 482-485.

About this article

Received
13 March 2015
Accepted
08 June 2015
Published
15 August 2015
SUBJECTS
Chaos, nonlinear dynamics and applications
Keywords
nonlinear oscillation
amplitude frequency formulation (AFF)
max min approach (MMA)
micro electro mechanical system