Anti-Symmetric Mode Vibration of Electrostatically Actuated Clamped-Clamped Microbeams for Mass Sensing.

This paper details study of the anti-symmetric response to the symmetrical electrostatic excitation of a Micro-electro-mechanical-systems (MEMS) resonant mass sensor. Under higher order mode excitation, two nonlinear coupled flexural modes to describe MEMS mass sensors are obtained by using Hamilton's principle and Galerkin method. Static analysis is introduced to investigate the effect of added mass on the natural frequency of the resonant sensor. Then, the perturbation method is applied to determine the response and stability of the system for small amplitude vibration. Through bifurcation analysis, the physical conditions of the anti-symmetric mode vibration are obtained. The corresponding stability analysis is carried out. Results show that the added mass can change the bifurcation behaviors of the anti-symmetric mode and affect the voltage and frequency of the bifurcation jump point. Typically, we propose a mass parameter identification method based on the dynamic jump motion of the anti-symmetric mode. Numerical studies are introduced to verify the validity of mass detection method. Finally, the influence of physical parameters on the sensitivity of mass sensor is analyzed. It is found that the DC voltage and mass adsorption position are critical to the sensitivity of the sensor. The results of this paper can be potentially useful in nonlinear mass sensors.

1. Introduction

Doubly clamped microbeam is a common resonant element in Micro-electro-mechanical-systems (MEMS) sensors [1,2,3]. Due to their great potential and unique characteristics, microbeam resonant sensors have the advantages of small, fast, high sensitivity [4,5]. Besides, only a small amount of power is required to operate [6,7]. However, the structure nonlinearity and nonlinear electrostatic force seriously affect the performance of conventional micro-mass sensor [8,9,10]. For example, the nonlinear electrostatic force can cause shifts in their resonant frequency and lead to error of the measured mass [11]. Recently, nonlinear MEMS mass sensors have attracted attention due to their unique advantages. Firstly, the nonlinear parameter identification method of mass sensors can solve the error caused by nonlinear stiffness [12]. Besides, the sensitivity and accuracy of the sensor can be improved by using frequency stability and amplitude jump in nonlinear vibration [13]. In this paper, we study the effect of added mass on anti-symmetric mode vibration and utilize the dynamic jump motion of anti-symmetric mode to propose a new mass detection method.

The mode coupled vibration can introduce rich nonlinear phenomena into the MEMS research and reveal the mechanism of the complex dynamic behaviors [14,15,16,17,18]. Anti-symmetric response can be induced by the mode coupled vibration. Li et al. [17] presented coupled vibration behavior between second order and third order modes caused by the axial stress, which can be used to suppress large amplitude vibration and to reduce the possibility of large deflection. Kirkendall et al. [19] reported multi-stable coupled vibration between resonant modes of an electroelastic crystal plate and introduced a mixed analytical-numerical approach to provide new insight into these complex interactions. Wang et al. [20] proposed a simplified oscillation system, which consists of two beam-shaped cantilevers. The results showed the possibility of doubling the frequency response signal from the low frequency cantilever to the high-frequency cantilever based on this super harmonic synchronization. Okamoto et al. [21] used anti-symmetric vibration in two coupled GaAs oscillators to realize the high-sensitivity charge detection. In contrast to the frequency-shift detection using a single oscillator, coupled vibration allowed a large readout up to the strongly driven nonlinear response regime. Hammad et al. [22,23] presented an analytical model and closed form expressions describing the response of mechanically coupled resonant structures. Mode coupled vibration was utilized to implement an adjustable filter. Younis et al. [24] studied possibility of activating a three-to-one internal resonance between the first and second modes. Besides, Mode coupled vibration can be used to improve the frequency stability of nonlinear systems [13]. Du et al. [25] reported the experimental observations of the internal resonance in a coupled ductile cantilever system from the viewpoint of mass sensing and disclosed a frequency enhancement mechanism. Hajjaj et al. [26] used the nonlinear hardening, softening and veering phenomena to realize a bandpass filter of sharp roll off from the passband to the stopband. In general, mode coupled vibration behavior was gradually applied to improve the performance of resonators and expand the scope of MEMS applications [27].

A MEMS resonant mass sensor mainly realizes the detection by changing the resonant frequency and vibration amplitude of the structure caused by the adsorption of the elastic element of the sensor to the target analyzers [28]. Frequency shift tracking is the most common method [29]. Bouchaala et al. [30] obtained analytical formulations to calculate the induced resonance frequency shifts caused by the added mass. The results indicated that the detection sensitivity increases with the decrease of size. However, with the reduction of size, there are obvious nonlinear effects and complex bifurcation behaviors [31], which seriously affect the dynamic mechanical characteristics and mass detection performance of the sensor. Therefore, nonlinear mass sensors were proposed. Younis et al. [32] utilized the dynamic bifurcations to realize novel methods and functionalities for mass detection. It was noted that bifurcation-based mass detection methods provided for dramatically enhanced sensitivity and less performance deterioration due to measurement noise as compared to frequency shift-based methods [33]. Similarly, Kumar et al. [34] also detailed proof-of-concept experiments on bifurcation-based sensing. Preliminary results revealed the bifurcation-based sensing technique to be a viable alternative to existing resonant sensing methods. Hasan et al. [35] studied the intelligent adjustable threshold pressure switch. When the pressure exceeds the threshold, the system can be induced to produce amplitude jump, realizing the rapid sensing of pressure. However, the dynamic behavior near the bifurcation point is easily disturbed by the ambient noise, which affects the stability of the sensor.

It can be concluded from the above analysis that mode coupled vibration can induce anti-symmetric modes and improve the performance of resonators, which may be beneficial to improve the accuracy and sensitivity of mass sensors [36,37,38]. Recently, most nonlinear mass sensors were realized by using the periodic saddle bifurcation of the resonant system, which greatly increases the sensitivity of the sensor. To the best of our knowledge, there are fewer quantitative results about a general analysis of nonlinear mass sensors by using anti-symmetric mode vibration. Besides, the mechanism of the effect of added mass on anti-symmetric mode vibration is not well understood. Through bifurcation analysis, the influence of added mass on transition mechanism of nonlinear jumping phenomena and complex nonlinear dynamic behaviors can be predicted, which motivates our present work. In this study we exploit the nonlinear jump phenomenon of the anti-symmetric mode to realize mass detection.

The structure of this paper is as follows. In Section 2, the Hamilton’s principle and Galerkin discretization is applied to obtain the equation of motion of a mass sensor. Then static analysis is carried out under different direct current (DC) voltage and added mass. In Section 3, the method of multiple scales is applied to produce an approximate solution. In Section 4, we analyze the physical conditions of the anti-symmetric mode vibration. Meanwhile, the effect of added mass on nonlinear dynamic behavior is considered. In Section 5, we propose the mass detection method by exploring the exploitation of amplitude jump behavior. In Section 6, the influence of DC voltage and mass adsorption position on the sensitivity of the sensor is considered. Finally, summary and conclusions are presented in the last section.

2. Problem Formulation

The bifurcation-based mass detection methods provide for dramatically enhanced sensitivity and less performance deterioration [32]. In our previous work, it was found that the anti-symmetric mode vibration of microbeams can be realized under the symmetrical electrostatic excitation [17]. Here, the dynamic jump behavior of anti-symmetric mode is utilized to realize mass detection. Figure 1 shows the schematic of a resonant mass sensor. The adsorption material is added to the microbeam. Then, a lumped mass m is added at . The added mass can affect the equivalent mass of the system and lead to a shift in the natural frequency of the resonator. The actuation of the microbeam is realized by means of a bias voltage and an alternating current (AC) voltage component.

Figure 1

Schematic of a resonant mass sensor.

By using Hamilton’s principle, the equation of motion that governs the transverse deflection is written as [17] with the boundary conditions where and .

Size parameters and physical properties of resonators are listed in Table 1. The last two terms of Equation (1) represent mid-plane stretching effect and the parallel-plate electric actuation which is composed of DC and AC components. A and I represent the area and moment of inertia of the cross section.

Table 1

Mass sensor parameters and physical properties.

Physical Parameter (units).	Value
Length of the beam electrode, L (µm)	150
Thickness of the beam electrode, h (µm)	1
Width of the beam electrode, b (µm)	10
Gap between the electrodes, d (µm)	1.5
Density of the electrode material, ρ (kg/m3)	2300
Young’s Modulus, E (GPa)	169
Dielectric constant of the medium, ε_0	8.85 × 10−12
Quality factor, Q	12.5

Then, the non-dimensional variables are introduced

Substituting Equation (3) into Equation (1,2), yields the non-dimensional equation of motion of the mass sensor with the boundary conditions

The parameters in Equation (4) are where represents added equivalent mass.

The microbeam deflection under an electric force is composed of a static component due to the DC voltage, denoted by , and a dynamic component due to the AC voltage, denoted by ; that is

Ignoring the time derivatives and the AC forcing term in Equation (4), we can obtain the static deflection of the microbeam

The static displacement of the resonator is very important to the natural frequency of the system. It was found that the nine-order mode discretization can accurately predict the static behavior of the resonator [9]. Here, the Galerkin method is introduced to calculate Equation (8). Meanwhile, the finite element method results are obtained from the software COMSOL (COMSOL Inc., Stockholm, Sweden) by using the Multi-field solver [39]. Figure 2 shows the relationship between midpoint deflections and the DC voltages obtained with the Galerkin method and the finite element method. Results are presented for values of ranging from 0 V to pull-in voltage, where the solid line represents the position of the potential well and the dotted line represents the position of the barrier. Finite element results verify the validity of the theoretical model. It should be noted that the added mass has no effect on the static displacement.

Figure 2

The maximum displacement of the microbeam under different voltages obtained by using the Galerkin method and COMSOL.

Substituting Equation (7) into Equation (4) and using Equation (8) to eliminate the terms representing the equilibrium position, the problem governing the dynamic behavior of the microbeam around the deflected shape is generated. To third-order in , the result is

Due to [17], is obtained.

The solution of Equation (9) can be expressed as , where is the i-th linear undamped mode shape of the straight microbeam. Then, the linear undamped eigenvalue problem is obtained

Substituting Equation (10) into the resulting Equation (9), multiplying by , and integrating the outcome from x = 0 to 1, yield where , .

Through Equation (11), we can obtain the resonant frequency where is the resonant frequency of the n-th order mode.

In our previous work, it was found that the third-order frequency is approximately equal to two times the second order frequency [17]. Hence we study the possibility of activating a 1:2 internal resonance between the second and third modes when the third mode is excited with a higher order excitation. Firstly, the effect of added mass on the second and third natural frequencies is considered, as shown in Figure 3. It should be noted that the effect of different adsorption positions on the natural frequency is very important. L1 = 50 µm and L1 = 75 µm are considered. When L1 = 75 µm, the added mass has a significant effect on the third modes. However, it has no effect on the second mode. Similarly, when L1 = 50 µm, the added mass has a significant effect on the second modes, and it has no effect on the third mode. To explain this phenomenon, the shape diagram of the second and third modes is obtained, as shown in Figure 4. When L1 = 75 µm, the added mass is at the node of the second-order mode. Thus, it has no effect on the second mode. Similarly, the added mass is at the node of the third-order mode under L1 = 50 µm.

Figure 3

Variation of the (a) second and (b) third natural frequencies of the mass sensors with various values of adsorption position and adsorption mass (the solid lines: the theoretical results; the rectangles: COMSOL results).

Figure 4

The second and third modes of the micro-mass sensor in COMSOL.

To study the dynamic behavior of the anti-symmetric mode, we take and obtain that where the dots indicate the time derivative and the parameters are given in “Appendix A”. represents equivalent added mass when L1 = 50 µm; represents equivalent added mass when L1 = 75 µm. When the driving frequency is close to two times the natural frequency of the second order mode, the second order amplitude can produce the bifurcation jump phenomenon [17]. We take advantage of the dynamic jump behavior of the anti-symmetric mode to realize mass detection.

3. Perturbation Analysis

The method of multiple scales is introduced to investigate the response of the mass sensor with small amplitude vibration. Then, ε is introduced as a small nondimensional bookkeeping parameter to indicate the significance of each term in the equation of motion. Considering the electrostatic force term , scaling the dissipative terms, we obtain

To describe the 1:2 internal resonance, detuning parameters and are defined by

Finally, the frequency response equation can be derived as

The detailed derivation process is given in “Appendix B”. In this paper, pseudo-trajectory processing method is introduced to solve Equations (17) and (18).

4. The Anti-Symmetric Mode Vibration

In this section, anti-symmetric mode vibration behaviors of MEMS mass sensors are considered. Firstly, we study the physical conditions of the anti-symmetric mode vibration.

4.1. Physical Conditions for Mode Transition

Anti-symmetric modes cannot be directly excited by symmetric driving forces. Mode coupled vibrations can be utilized to induce anti-symmetric modes. When the driving frequency is far away from two times the second natural frequency or the electrostatic excitation is very small, there is no second order vibration [17]. The bifurcation analysis is introduced to obtain the physical conditions of the second order vibration. When the second order amplitude occurs, we can obtain the following equation by Equation (17)

Then, the threshold of third order amplitude can be obtained by solving Equation (19)

Because the displacement of the second order vibration mode at the midpoint of the microbeam is always 0, the added mass has no effect on the critical amplitude when L1 = 75 µm.

When third order amplitude is more than the above threshold, the mode coupled vibration can occur. From Equation (20), the precondition of anti-symmetric mode vibration is obtained.

When the anti-symmetric mode vibration occurs, becomes unstable. Then, substituting into Equation (17) and Equation (18), the physical conditions of anti-symmetric mode vibration under different DC voltage and added mass are obtained, as shown in Figure 5. The unstable region represents the occurrence of anti-symmetric mode vibration. It is noted that: (1) with the increase of DC voltage, the minimum critical frequency of anti-symmetric mode vibration decreases; (2) the added mass cannot affect the minimum critical frequency of anti-symmetric mode when L1 = 75 µm; (3) the added mass can decrease the minimum critical frequency of anti-symmetric mode when L1 = 50 µm; (4) the anti-symmetric mode vibration behavior depends heavily on DC voltage and added mass. Therefore, we can use the bifurcation behavior of anti-symmetric mode vibration to detect the mass.

Figure 5

Physical conditions of anti-symmetric mode vibration under different DC voltage and added mass. (The red dotted line represents the precondition of anti-symmetric mode vibration).

4.2. Stability of the Nontrivial Solution

To further study the bifurcation behavior when the anti-symmetric mode vibration occurs, the stability analysis of the nontrivial solution is introduced. The subcritical bifurcation can lead to unstable branches near the critical points. On the contrary, the supercritical bifurcation can lead to stable branches. To determine the stability of periodic vibration, we ignore the higher-order nonlinear terms and obtain the following equation by Equation (17)

Substituting Equation (20) into Equation (22) yields the discriminant

The case results in the subcritical bifurcation. With the increase of AC voltage, the second order amplitude suddenly appeared. Likewise, the case results in the supercritical bifurcation. With the increase of AC voltage, the small vibration in the second order mode appears.

Figure 6 shows variation of the bifurcation behavior versus , DC voltage and added mass. In the yellow area, the dynamic jump behavior of anti-symmetric mode vibration occurs with the increase of AC voltage. It is found that the increase of the DC voltage makes the dynamic jump behavior occur easily. It is interesting to note that low frequency perturbation parameter is more advantageous to realize the dynamic jump behavior than the high frequency perturbation parameter. Meanwhile, the added mass has obvious influence on the bifurcation behavior of the system when L1 = 50 µm. As the added mass increases, the critical frequency of subcritical bifurcation decreases. However, the added mass has no influence on the critical frequency of subcritical bifurcation when L1 = 75 µm.

Figure 6

Variation of the bifurcation behavior versus , DC voltage and add mass. (The solid line represents the critical condition of anti-symmetric mode vibration).

4.3. Dynamic Analysis

To further study the effect of added mass on anti-symmetric mode vibration behaviors, dynamic behaviors of the anti-symmetric mode under different adsorption position and mass are introduced.

L1 = 50 µm

Firstly, L1 = 50 µm is considered. As shown in Figure 6, point A and point B represents supercritical bifurcation and subcritical bifurcation, respectively. Then, Figure 7 shows the force-amplitude curves corresponding to point A and point B. When μg, the supercritical bifurcation voltage Vac1 is obtained. With the increase of the driving voltage, the second order amplitude appears gently. When μg, the bifurcation voltage increases and the supercritical bifurcation is transformed into subcritical bifurcation. The subcritical bifurcation voltage Vac2 is obtained and the second order amplitude is suddenly generated with the increase of the driving voltage. To validate the above analysis, long-time integration (LTI) of Equation (14) is used to obtain some numerical solutions (discrete points), compared with the analytical solution derived from the method of multiple scales.

Figure 7

Comparison of the force-amplitude curves obtained by theoretical method (line) and numerical method (rectangle) corresponding to A-B in Figure 6 when L1 = 50 µm (solid line: stable; dashed line: unstable).

Then, Figure 8 shows the effect of added mass on amplitude-frequency response curve. There are four critical frequencies for anti-symmetric mode vibration behavior without considering the added mass. P1, P2, P3 and P4 indicate the bifurcation frequency of the anti-symmetric mode vibration as shown in Figure 5c. It is found that: (1) when the frequency is less than P1, no anti-symmetric mode vibration occurs; (2) when the frequency is between P1 and P2, the anti-symmetric mode vibration occurs and there is only one stable periodic solution in the system; (3) when the frequency is between P2 and P3, the trivial solution of the second order amplitude becomes stable. There are two stable and one unstable periodic solutions in the system and the second order amplitude depends on the initial conditions; (4) when the frequency is between P3 and P4, the anti-symmetric mode vibration occurs and there are two stable periodic solutions in the system; (5) when the frequency is more than P4, the anti-symmetric mode vibration behavior disappears. Then, there are two critical frequencies for anti-symmetric mode vibration behavior with considering added mass. When the frequency is between F1 and F2 (Figure 5d), the anti-symmetric mode vibration occurs and there are two stable periodic solutions in the system. F1 represents the supercritical bifurcation type and F2 represents the subcritical bifurcation type. As the drive frequency decreases, the second order amplitude can be suddenly generated when the frequency is equal to F2. From Figure 8, it is found that the added mass can change the number of critical frequencies for anti-symmetric mode vibration behavior. Meanwhile, the added mass reduces the subcritical bifurcation frequency.

Figure 8

Amplitude-frequency response curves of anti-symmetric mode when m = 0 µg and m = 5 × 10−6 µg (line: theoretical method; rectangle: numerical method).

L1 = 75 µm

Similarly, L1 = 75 µm is considered. From Equation (23), the added mass has no effect on bifurcation type when L1 = 75 µm. However, the added mass can affect the bifurcation voltage and frequency. Figure 9 shows the force-amplitude curves obtained by pseudo-trajectory processing method (line) and long-time integration method. With the increases of added mass, the subcritical bifurcation voltage increases and the supercritical bifurcation voltage decreases. Then, amplitude-frequency response curves of the system considering the different added masses are introduced, as shown in Figure 10. There are two critical frequencies for anti-symmetric mode vibration behavior. The increase of the added mass shifts down the forced frequency response, ΔΩ being the subcritical bifurcation frequency shift.

Figure 9

Comparison of the force-amplitude curves with considering the different added mass when L1 = 75 µm.

Figure 10

Amplitude-frequency response curve of the anti-symmetric mode considering the different added mass when L1 = 75 µm (line: theoretical method; rectangle: numerical method).

When L1 = 50 µm and L1 = 75 µm, the added mass can increase the inertial forces of the second order vibration mode and the third order vibration mode, respectively. Through Figure 7, Figure 8, Figure 9 and Figure 10, the effect of added mass on anti-symmetric mode vibration behavior can be summarized as follows: (1) the added mass can adjust the number of bifurcation points near the origin; (2) the supercritical bifurcation near the origin can be transformed into subcritical bifurcation by the added mass; (3) the added mass can reduce the bifurcation frequency of the dynamic jump motion near the origin and increase the bifurcation voltage of the dynamic jump motion near the origin.

5. Mass Detection Method

In Section 4, we study the effect of added mass on anti-symmetric mode vibration behavior in detail. Subcritical bifurcation can lead to the dynamic jump motion, which greatly improves the sensitivity of the sensor. It is found that the added mass has an important influence on the bifurcation voltage and frequency of the dynamic jump motion. Then, we utilize of amplitude jump behavior to realize mass detection in MEMS.

First of all, the principle of mass detection is given as follow:

(1) It is found that when the driving voltage increases and the driving frequency decreases, the dynamic jump motion of the anti-symmetric mode can occur. Through the force-amplitude curve and amplitude-frequency response curve, we can obtain the bifurcation frequency and voltage.

(2) Then, we consider the mass identification formula in the case of L1 = 50 µm and L1 = 75 µm, respectively.

L1 = 50 µm

From Figure 7, Figure 8, Figure 9 and Figure 10, the third order amplitude is very small when the subcritical bifurcation occurs. Thus, we can ignore the nonlinearity of the third order mode and obtain by Equation (18)

Then, substituting Equation (24) into Equation (19) yields

L1 = 75 µm

Substituting Equation (20) into Equation (18) yields where a3 can be obtained by Equation (20).

(3) Finally, we can obtain the added mass m with dimensional transformation.

The numerical studies are introduced to prove the mass identification method. Firstly, L1 = 50 µm is considered. Figure 11 show force-amplitude curve obtained by sweeping up the AC voltage. As the added mass increases, the bifurcation voltage increases when . With Equations (25) and (27), we obtain three kinds of parameter identification results, as shown in Table 2. The results show that the mass detection method presented in this paper can identify the added mass. Figure 12 show amplitude-frequency response curve of the system considering the different added mass obtained by sweeping down the frequency. As the drive frequency decreases, there are jump phenomena in the second order amplitude. Besides, the added mass reduces the bifurcation frequency of the jump point. Then, swept harmonic responses for midpoint displacement are introduced to further verify the bifurcation jump phenomenon, as shown in Figure 13. Here, we also obtain three kinds of parameter identification results by using Figure 12, as shown in Table 3. Similarly, Figure 14 and Figure 15 show force-amplitude curve and amplitude-frequency response curve in the case of L1 = 75 µm. With Equations (26) and (27), we obtain mass identification results, as shown in Table 4 and Table 5. There is a small error between the identification result and the real result. There are two main error sources: (1) The nonlinear stiffness terms of the third order modes are ignored when L1 = 50 µm; (2) the accuracy of numerical studies depends on the density of discrete points. An insufficient number of discrete points may lead to identification error.

Figure 11

Comparison of the force-amplitude curves obtained by long-time integration method considering the different added mass when L1 = 50 µm.

Table 2

Three groups of mass detection results obtained by the force-amplitude curve when , .

Number	The True Mass m (10−6 µg)	Bifurcation Voltage (V)	Identification Results m (10−6 µg)	Error (10−6 µg)
1	1	0.537	0.995	0.005
2	3	0.774	2.833	0.167
3	5	1.032	4.624	0.376

Figure 12

The amplitude-frequency response curve of the system considering the different added mass when L1 = 50 µm.

Figure 13

Swept harmonic responses for midpoint displacement when V and V.

Table 3

Three groups of mass detection results obtained by the amplitude-frequency response curve when , .

Number	The True Mass m (10−6 µg)	Bifurcation Frequency (100 kHz)	Identification Results m (10−6 µg)	Error (10−6 µg)
1	1	21.598	0.949	0.051
2	3	21.591	2.772	0.228
3	5	21.585	4.542	0.458

Figure 14

Comparison of the force-amplitude curves obtained by long-time integration method considering the different added mass when L1 = 75 µm.

Figure 15

Amplitude-frequency response curve of the system considering the different added mass when L1 = 75 µm.

Table 4

Three groups of mass detection results obtained by the force-amplitude curve when , .

Number	The True Mass m (10−6 µg)	Bifurcation Voltage (V)	Identification Results m (10−6 µg)	Error (10−6 µg)
1	1	0.482	1.191	0.191
2	3	0.538	3.219	0.219
3	5	0.602	5.342	0.342

Table 5

Three groups of mass detection results obtained by the amplitude-frequency response curve when , .

Number	The True Mass m (10−6 µg)	Bifurcation Frequency (100 kHz)	Identification Results m (10−6 µg)	Error (10−6 µg)
1	1	21.601	0.955	0.045
2	3	21.597	3.112	0.112
3	5	21.592	4.273	0.727

6. Sensor Sensitivity

In this section, considering different added positions and DC voltages, we study the sensitivity of mass sensors. Variation of bifurcation voltage and bifurcation frequency caused by the added mass is introduced to characterize the sensitivity of the sensor. Figure 16 and Figure 17 show the effect of added mass at different positions on bifurcation voltage and bifurcation frequency. As the DC voltage increases, the sensitivity of the mass sensor in the case of L1 = 75 µm is improved and the sensitivity of the mass sensor in the case of L1 = 50 µm is suppressed. Thus, we need to select a reasonable DC drive voltage and adsorption position to improve the sensitivity of the mass sensor.

Figure 16

The effect of added mass at different positions on AC voltage of bifurcation jump points.

Figure 17

The effect of added mass at different positions on frequency of bifurcation jump points.

7. Conclusions

This paper studies the anti-symmetric mode vibration behavior of MEMS mass sensors, obtains the influence mechanism of added mass on anti-symmetric mode vibration behavior and utilizes the dynamic jump motion of anti-symmetric mode vibration to realize the mass detection.

Firstly, the Hamilton’s principle and Galerkin discretization are applied to obtain two degrees of freedom equations of the resonant structure. The results show that the added mass at different adsorption positions has an important effect on the natural frequency of the resonator. Through bifurcation analysis, physical conditions of anti-symmetric mode vibration behavior are obtained. The anti-symmetric mode vibration behavior depends heavily on DC voltage and added mass. Besides, it is found that the added mass can change the bifurcation type of the resonator and affect the bifurcation voltage and bifurcation frequency of the anti-symmetric mode vibration.

The dynamic jump motion is introduced to realize mass detection, which greatly improves the sensitivity of the mass sensor. The parameter identification formula of mass detection is deduced and numerical studies are used to verify mass parameter identification method. Finally, the sensitivity of the sensor is analyzed. With the increase of DC voltage, the sensitivity of the mass sensor is improved when L1 = 75 µm. On the contrary, with the increase of DC voltage, the sensitivity of the mass sensor is suppressed when L1 = 50 µm. Therefore, we need to adjust the DC drive voltage of the resonator according to different mass adsorption positions. The framework presented here provides theoretical support for nonlinear mass sensors. It should be emphasized that all the theoretical results in this paper are compared with numerical results, which guarantees the accuracy of our whole investigations.

It is worthy to mention that the proposed mass detection method discussed in this paper need to be investigated for their stability to external disturbances. Fork bifurcation occurs near the critical voltage and frequency. Hence, the stability of operation prior to mass detection must be ensured to prevent accidental bifurcation jump phenomenon due to noises or disturbances. This can be studied by conducting global dynamic analysis to track the basin of attraction of the stable solution.