Abstract

This paper focuses on chaos suppression strategy of a microresonator actuated by two symmetrical electrodes. Dynamic behavior of this system under the case where the origin is the only stable equilibrium is investigated first. Numerical simulations reveal that system may exhibit chaotic motion under certain excitation conditions. Then, bifurcation diagrams versus amplitude or frequency of AC excitation are drawn to grasp system dynamics nearby its natural frequency. Results show that the vibration is complex and may exhibit period-doubling bifurcation, chaotic motion, or dynamic pull-in instability. For the suppression of chaos, a novel control algorithm, based on an integer-order nonsingular fast terminal sliding mode and a fractional-order switching law, is proposed. Fractional Lyapunov Stability Theorem is used to guarantee the asymptotic stability of the system. Finally, numerical results with both fractional-order and integer-order control laws show that our proposed control law is effective in controlling chaos with system uncertainties and external disturbances.

1. Introduction

Microelectromechanical systems (MEMS) have been developed rapidly in the field of sensing and actuating ‎[14]. However, as one of the most undesirable dynamic behaviors, chaotic motion exists extensively in dynamic MEMS, such as microresonators ‎[5], comb-drivers ‎[6], or atomic force microscopes (AFM) ‎[7]. In many cases, the existence of chaos is highly undesirable as it can lead to poor performance ‎[8, 9]. Therefore, suppression of chaos of microcomponents has attracted much attention and is studied by many groups. For instance, applying Melnikov method, Ashhab et al. ‎[10] derived the chaotic region of a microcantilever sample system with a proportional and derivative controller and finally eliminated the possibility of chaos. Yamasue and Hikihara ‎[11] introduced a time delayed feedback control method to suppress the chaotic motion of a microcantilever sensor. An optimal linear control method was used to destroy the chaotic motion of a time-varying capacitor ‎[12]. In the work by Wang et al. ‎[9], a robust algorithm composed of sliding model control and a backstepping feedback was proposed and used to realize the suppression of chaos of an AFM system. Moreover, a nonlinear mechanical model for MEMS-based storage devices was developed and two control strategies: static output feedback control and geometric nonlinear control, were applied to eliminate the chaotic motion ‎[13]. Tusset et al. ‎[14] used three control strategies: state dependent Riccati equation control, optimal linear feedback control, and fuzzy sliding mode control, to sufficiently suppress the chaotic behavior of a comb-driver. Seleim et al. ‎[15] investigated the dynamics of a close-loop controlled electrostatic MEMS resonator and finally determined the optimal operating regions beyond chaos.

Recently, chaotic motion of a doubly clamped microbeam actuated by two symmetrical electrodes was investigated by Haghighi and Markazi ‎[16]. As the existences of cubic stiffness and electrostatic force, system may have different number of equilibrium points. They paid their attention to the vibration case where the origin is an unstable saddle point. Results showed that homoclinic bifurcation could lead to system chaos. Then, a robust adaptive fuzzy control algorithm was applied to control the chaotic motion and hence stabilize it into a high-amplitude oscillation state. What follows are some other control strategies, such as fuzzy sliding mode control ‎[17] and second-order fast terminal sliding mode control ‎[8], that were also used for the suppression of chaos. What is more, Aghababa ‎[18] investigated the chaotic motion of a fractional-order model of this MEMS resonator and introduced a novel fractional finite-time controller to suppress the chaotic behavior with system uncertainty and external disturbance. Although a fractional-order MEMS system is chosen as a research model, this work indicates the possibility of applying fractional-order controllers to suppress chaotic motion in MEMS resonator.

Fractional calculus has a long history ‎[19]. However, for many years, this branch of science had not been applied in physics and engineering. Until recent decades, the investigation of fractional calculus has been growing very rapidly, among which design of fractional-order controller has become an attractive research topic ‎[20, 21]. Meanwhile, sliding mode control is an efficient approach in the application of linear or nonlinear control. For the foregoing reasons, there are more and more new control methods composed of fractional calculus and sliding mode control ‎[22, 23], and their application fields have been gradually expanded from integer-order dynamic systems to fractional-order ones ‎[2426]. Besides, compared with traditional integer-order control, fractional-order control can offer more degrees of freedom to designers to meet a predefined set of performance criteria ‎[27]. If the fractional-order selection is appropriate, it can exhibit faster convergence and stronger robustness than that of integer-order one ‎[28].

So far, to the best of our knowledge, there is still no chaotic analysis of this doubly clamped microbeam-based resonator when system includes only one stable origin equilibrium. Although this vibration situation had been thought to have no chaotic motion under small disturbances ‎[29], we found its chaos phenomenon as the increase of AC voltage during working pattern. Besides, there are fewer fractional-order sliding mode control results about chaos suppression of this MEMS resonator. In this paper, a simulation case where the origin is the only stable center is first carried out to investigate the existence of chaos. Then, bifurcation diagram is used to study the chaotic behavior nearby its primary resonance frequency. For the suppression of chaos, a fractional-order nonsingular fast terminal sliding mode control method (FNFTSMC), based on an integer-order nonsingular fast terminal sliding mode ‎[30] and a fractional-order switching law ‎[22], is proposed and, then, the effectiveness and robustness in controlling chaotic motion are investigated with system uncertainty and external and stochastic disturbances (usually existing in practice).

The structure of this paper is as follows. In Section 2, some basics of fractional calculus are briefly reviewed. In Section 3, dynamic model is introduced ‎[16]. Chaotic analysis under the case where the origin is the only stable equilibrium is then carried out. In Section 4, a FNFTSNC method is proposed and its asymptotical stability is then proved via the Fractional Lyapunov Stability Theorem. In Section 5, simulation results are presented to verify the effectiveness of the proposed control scheme. Finally, Section 6 presents the conclusion.

2. Some Basics of Fractional Calculus

In this section, some basic definitions and preliminaries of fractional calculus are presented first. Then, a numerical approximation of fractional operators in the frequency-domain is given.

2.1. Definitions and Preliminaries

There is a uniform formula for the fractional integral of function with th-order which is defined as ‎[19]where , and are the limits of the operation, and is the well-known Gamma function.

Two important and commonly used definitions of the fractional derivatives are the Riemann-Liouville and the Caputo fractional derivatives. In this paper, the Riemann-Liouville definition is adopted and denoted as , which is defined as ‎[19]where , .

Property 1 ‎(see [19]). For the Riemann-Liouville derivative, the following equality holds:where , , .

Property 2 ‎(see [31]). The following equality holds for the Riemann-Liouville derivative:where .

Property 3. According to Properties 1 and 2, if and , one has

Property 4 ‎(see [22]). If , one haswhere is the sign function.

Theorem 1 ‎(see [26, 32, 33]). Let be an equilibrium point for the nonautonomous fractional-order systemwhere satisfies the Lipschitz condition with Lipschitz constant and . Assume that there exists a Lyapunov function satisfyingwhere , , , and are positive constants and denotes an arbitrary norm. Then the equilibrium point of system (7) is asymptotically stable.

2.2. Numerical Approximation

In practice, it is difficult to obtain the exact solution of fractional differential equations and, in most cases, is unknown. These make numerical approximations be widely used. Among these numerical methods, Oustaloup approximation is one of the most commonly used ways to obtain fractional operators ‎[34]. In 2006, Xue et al. ‎[35] proposed a refined Oustaloup filter, which has a very good approximation of the fractional operators in a specified frequency range of order . This filter is given bywhere is the order of fractional operation; and can be computed fromwhere and are weighted parameters. It is expected that a good approximation is obtained with and .

3. System Description

As is shown in Figure 1, an electrostatically actuated microresonator is simplified as a one degree of freedom system ‎[16]. is the vertical displacement of the microbeam; is the initial gap width. The resonator is deflected by DC voltage and actuated by AC voltage , where and are the amplitudes of DC and AC voltages and is the angular frequency of AC voltage. Displacement and velocity of this resonator can be identified through the measurement of output voltage in practice ‎[36].

The electrostatic force loaded on the microresonator can be expressed as ‎[36]where is the capacitance over the gap when , is the dielectric constant of the air, andis the overlapped area between the microbeam and the fixed electrode. The governing equation of motion of the resonator can be expressed as ‎[16]where the prime represents the derivative with respect to time and , , , and are the effective mass and damping, linear, and cubic stiffness, respectively.

For convenience, one can introduce the following nondimensional variables:where is the natural angular frequency of the microbeam. Then, (12) can be rewritten in the following nondimensional form (without any simplification):where the overdot represents the derivative with respect to nondimensional time .

In our previous study ‎[37, 38], we find that the number of equilibrium points of this system is determined by both and . Under the situation with , this system has only one stable center at the origin when and two stable centers at either side of the origin when . Taking as an example, the system has two stable centers when . Under some certain AC excitations, system may exhibit chaotic motion ‎[16]. In this paper, is considered here to investigate the chaotic behavior of this microresonator when system only has one stable origin center.

Numerical simulations are carried out as follows. Nondimensional parameters are assumed to be , , and . Potential energy and phase portrait of Hamiltonian system of (14) are shown in Figure 2. In this case, system contains only one stable origin equilibrium point . Besides, two symmetrical unstable equilibrium points are located at either side of it. Theoretical maximum vibration amplitude is approximate to . Generally, microresonator is driven to vibrate around its natural frequency ‎[39, 40]. Through some calculation, one can derive that the linearized equivalent natural angular frequency of this microresonator is in nondimensional form. Therefore, AC excitation angular frequency should be selected in the neighborhood of . For case study, excitation parameters are selected as and . Numerical solutions of system (14) are calculated by removing 5000 periods () of AC excitation and retaining the last 128 periods. Initial conditions are assumed to be . Finally, phase portrait and time history diagram of (14) are shown in Figure 3. The red lines in Figure 3(a) are two heteroclinic orbits, representing stable and unstable boundary. Inner region is stable beyond pull-in instability. In Figure 3(b), two red lines indicate the available motion region of the microresonator. From Figure 3, one can notice that system vibration is quite complex. Detailed numerical analyses (Figure 4) via Poincare maps, power spectral density (PSD), and maximum Lyapunov exponents (MLE) demonstrate that the microresonator indeed undergoes chaotic motion. Thus, system (14) may exhibit complex dynamics such as chaos during working pattern (namely, the resonator is excited near its natural frequency).

In what follows, bifurcation diagram is drawn to analyze the chaotic behavior of this resonator, which is obtained based on Poincare section taken by sampling the data at each period of AC excitation. We consider the following system parameters ‎[41]:  kg, μN/μm, μN/μm3, μm, and  F. A quality factor is assumed to be , which is related to the damping coefficient by [40]. Taking  V yields . Then, the dimensional equivalent natural angular frequency is approximate to 4 × 105 rad/s (63.66 kHz). Therefore, the following analysis mainly focuses on the vibration nearby its natural frequency.

In Figure 5, dynamic behavior of the system is depicted as the variation of AC voltage amplitude. As the increase of , the resonator exhibits complex dynamics like Period 2 motion, Period 3 motion, period-doubling bifurcation, chaotic motion and so on. When 7.59  V, the system is almost chaotic. When  V, dynamic pull-in occurs. Next, the effect of frequency of AC excitation on chaotic motion is studied. Taking  Vyields a Period 1 motion (Figure 5). In Figure 6, the vibration characteristic is complex as well as when the angular frequency of AC voltage varies from 3 × 105 rad/s (47.75 kHz) to 5 × 105 rad/s (79.58 kHz), especially between 4 × 105 rad/s (63.66 kHz) and 4.4 × 105 rad/s (70.03 kHz). When varies from 4.03 × 105 rad/s to 4.12 × 105 rad/s, system dynamic changes between Period 1 and chaotic motion. In domain [4.2 × 105 rad/s, 4.4 × 105 rad/s], system is always chaotic.

The above analysis indicates that nearby the natural frequency, the vibration of the system is complex and may exhibit chaotic motion. In the following section, a FNFTSMC method is proposed to cope with this chaotic behavior.

4. Control Method

Consider the following integer-order system:where is a known function, and are uncertainty term and external disturbance, and is the presented control law.

If the tracking errors of the controlled system are defined as and , where is the desired trajectory, then the error-dynamic system can be given as

Assumption 2. Uncertainty and disturbance are usually unknown and assumed to be bounded as follows:where is a positive constant.

Sliding mode control design usually consists of two steps: (a) sliding surface design; (b) effective control law design. In the first step, a uniform integer-order sliding surface is introduced as follows ‎[30]:where and are positive constants, , .

Remark 3. Expression (18) does not have singularity problem and, for any given initial condition , the system state converges to in finite time ‎[30].

Sliding mode control law can be decomposed into two parts, equivalent control law and switching law. The former one, deduced according to predesigned sliding surface by the application of sliding mode control theory, guarantees system state on sliding surface, while the latter one constrains the state to sliding surface subsequently. Firstly, without any consideration of uncertainty and disturbance, the equivalent control law can be given asNext, two different forms of switching law will be discussed as follows.

4.1. Integer-Order Nonsingular Fast Terminal Sliding Mode Control (INFTSMC)

A traditional integer-order switching law can be described as ‎[30]where , , .

Theorem 4 ‎(see [30]). If control law is applied to the error-dynamic system (16), then the tracking error will converge to zero in finite time if .

Proof. Consider the smooth Lyapunov function asTaking the derivative of (21) with respect to time and using (16) and (18), one hasConsequently, the tracking error will converge to zero.

4.2. Fractional-Order Nonsingular Fast Terminal Sliding Mode Control (FNFTSMC)

Before introducing the fractional-order switching law, a fractional-order reaching law must be described in advance. The expression of reaching law can be given by the following ‎[22]:where .

According to Property 3, the following operation holds:The above equation is equivalent toDue to , . The reaching dynamics governed by (25) will create a stronger push from both sides of the switching manifold ‎[22].

According to the above analysis, a fractional-order switching law can be proposed and defined aswhere .

In our study, a FNFTSMC method based on an integer-order sliding mode and a fractional-order switching law is proposed. The feasibility will be discussed as follows.

Theorem 5. Fractional-order control law is effective in controlling system (16), if .

Although the error-dynamic system (16) is integer-order, the proposed sliding mode control is fractional-order. Therefore, it is more suitable to prove the occurrence of the sliding motion via the Fractional Lyapunov Stability Theorem (Theorem 1) ‎[41].

Proof. Select a nonsmooth Lyapunov function asThen, one hasAccording to Property 4 and (25), one has the following ‎[22]:Then (28) can be simplified asAs the sliding surface is convergent ‎[30], it is reasonable to assume that is bounded in, where is a positive constant. Then, (30) can be written asTherefore, satisfies fractional-order Lyapunov stability theorem (Theorem 1). The proof is completed.

Remark 6 ‎(see [42]). Nonsmooth Lyapunov functions can be applied to show the finite-time convergence of a system.

The signum function of the switching law is one main inducement of chattering. When the system state hits the sliding surface, chattering will occur because of the discontinuity of signum function. To reduce chattering, a saturation function is often used instead of signum function. Here, a new form of saturation function is introduced as follows ‎[43]:where , is the boundary layer, and .

5. Simulation Results

In order to demonstrate the effectiveness of the proposed control method, simulation results are provided in this section. Numerical simulations based on (14) are carried out using the fourth-order Runge-Kutta method with 0.001 time-step size.

Control flow chart is shown in Figure 7. Error-dynamic signals ( and ) can be calculated through the comparison of system signals ( and ) and desired trajectories ( and ). Then, fractional-order controller can be obtained through two steps. Integer-order operator can be used to calculate equivalent control law (19) and fractional-order switching law (26) can be calculated through a fractional-order filter (9). Meanwhile, traditional integer-order controller can be calculated just through an integer-order operator. In practice, displacement and velocity signals of the microresonator can be obtained by transformation. Through calculation, control can now be obtained and loaded on this microresonator in the form of external force, like electrostatic force (electrode is acted by feedback voltage), piezoresponse force (microbeam attached by piezoelectric layers), and so forth. As our main focus is on the control strategy design and its application in the field of chaos suppression, detailed implementation is not discussed here anymore.

Assume  V and  rad/s; other parameters are the same as Section 3. Under these parameters, the resonator evolves into chaotic state. The control target is to let the motion of the resonator track the desired trajectory . In practice, system uncertainties and external disturbances are usually the most common disturbances in MEMS. Thus, to verify the effectiveness and robustness of the proposed method, the uncertainty and disturbance are assumed to be as follows: (a) Case : and ; (b) Case : and , where stochastic parameter follows the standard uniform distribution on the open interval . The controller is activated at . Besides, a small stochastic disturbance is added onto the system when . follows the standard uniform distribution on the open interval . Control parameters in these two cases are given as , , , , , , , , , , , , , , and .

Simulation results are shown in Figures 8 and 9. The black line represents the trajectory of the traditional INFTSMC with signum function switching law ‎[30], and the blue line represents the trajectory of our proposed FNFTSMC with novel saturation function switching law. From these two figures, one can notice that(1)compared with the control method by L. Yang and J. Yang ‎[30], our proposed control method is effective in eliminating the chaotic motion;(2)under this set of control parameters, the trajectories of the proposed FNFTSMC converge slightly faster than that of INFTSMC;(3)the novel saturation function is effective in reducing chattering;(4)the proposed method has strong robustness to the system (random) uncertainty and disturbance.

Appropriate fractional-order selection can make FNFTSMC exhibit faster convergence and stronger robustness ‎[22]. Therefore, our next focus is mainly on a qualitative discussion about the effect of fractional-order on control performance.

Standard deviation can reflect the degree of dispersion of a data set, so this statistic can be used to describe the degree of chattering of control signal . The larger the standard deviation of input is, the more dramatic the chattering is. Meanwhile, a nondimensional statistical function can be defined to approximately estimate the control energy cost. It should be noticed that no matter what the sign of the work done by the controller is, the control system is always energy-intensive. Thus, a mean power to describe the energy cost can be defined aswhere represents the total calculated time from control in action to the end, is the discrete point number when control is on, and is equal to the total length of discrete points. When control is activated, control signal can be calculated according to the difference between the actual trajectory and the desired trajectory . Under the effect of excitation and control input , the trajectory can be driven from to in one time-step. Thus, the work done by the control is approximate as . That is to say, the control energy cost is . By that analogy, one can obtain (33).

For convenience, stochastic disturbance is not considered. Control parameters are the same as Case . Assume the control is activated at and the total simulation time is . Control responses under different fractional-order are shown in Figure 10. It is obvious from the figure that the trajectories converge faster with the increase of fractional-order . FNFTSMC are more sensitive than INFTSMC when control is on. And it can also achieve a better “soft landing” to the desired trajectory. However, the increase of can increase the possibility of chattering and power cost (Table 1). Anyway, compared with integer-order control, both the standard deviation of input and mean power of FNFTSMC are much less than those of INFTSMC. During control parameters selection, we should make reasonable trade-offs between convergence speed, chattering, and power cost. However, this is beyond the scope of our research.

6. Conclusions

This paper presents the chaotic dynamics of a doubly clamped microbeam actuated by two symmetrical electrodes. Different from the research by Haghighi and Markazi ‎[16], this study mainly focuses on system dynamics when this MEMS resonator includes only one stable origin center. Through simulation results, we find, under primary resonance condition, that the resonator may exhibit chaotic behaviors as well. For the suppression of chaos, a novel control strategy, composed of an integer-order nonsingular fast terminal sliding surface with a fractional-order switching law, is proposed and shows its effectiveness in controlling chaotic motion. Besides, the existence of uncertainty and external and stochastic disturbance also verifies its strong robustness. Moreover, a qualitative analysis on the selection of fractional-order is also carried out. Two statistical functions, standard deviation of input and mean power, are used to estimate the degree of chattering and energy cost. Simulation results (Figure 10 and Table 1) reveal that the increase of fractional-order is accompanied by relatively higher chattering and bigger energy cost, but much faster convergence.

This study develops a fractional-order sliding mode control method for chaos control of MEMS resonator. The feasibility is numerically verified by comparing it with an integer-order sliding mode control ‎[30]. However, control properties with different control parameters still need to be investigated in order to find a group of optimal parameters. Future work could focus on control parameter optimization by applying some optimization theories.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work was supported by the National Natural Science Foundation of China (Grants nos. 11372210, 11302149), Tianjin Research Program of Application Foundation and Advanced Technology (Grant no. 15JCQNJC05000), Innovation Team Training Plan of Tianjin Universities and Colleges (TD12-5043), Tianjin College Foundation of Development of Science and Technology (20130408), and Development Foundation of Tianjin University of Technology and Education (KJ11-06).