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.


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. 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.

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 x = L 1 . 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. 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.

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.
. By using Hamilton's principle, the equation of motion that governs the transverse deflection ( , ) w x t is written as [17] with the boundary conditions 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.

Physical parameter (units).
Value Length of the beam electrode, L (µm) 150 Thickness of the beam electrode, h (µm) 1 By using Hamilton's principle, the equation of motion that governs the transverse deflection w(x,t) is written as [17] [ρA + δ(x − L 1 )m] ..
with the boundary conditionsŵ where .ŵ = ∂ŵ ∂t andŵ = ∂ŵ ∂x . 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.
Then, the non-dimensional variables are introduced 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 w dc (x), and a dynamic component due to the AC voltage, denoted by w ac (x); that is w = w dc + w ac 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 V dc 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. 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 ac w , the result is The solution of Equation (9) can be expressed as 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 i φ , and integrating the outcome from x = 0 to 1, yield 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 w ac , the result is ..
Due to V dc >> V ac [17], (V dc + V ac cos Ωt) 2 ≈ V 2 dc + 2V dc V ac cos Ωt is obtained. The solution of Equation (9) can be expressed as w ac ( 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 φ i , and integrating the outcome from x = 0 to 1, yield ..
(1−w dc ) 2 , η n = φ 2 n (L 1 /L)m/ρAL. Through Equation (11), we can obtain the resonant frequency where χ n 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. L 1 = 50 µm and L 1 = 75 µm are considered. When L 1 = 75 µm, the added mass has a significant effect on the third modes. However, it has no effect on the second mode. Similarly, when L 1 = 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 L 1 = 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 L 1 = 50 µm.

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

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 To study the dynamic behavior of the anti-symmetric mode, we take w ac ( ..
where the dots indicate the time derivative and the parameters are given in "Appendix A". η 2 represents equivalent added mass when L 1 = 50 µm; η 3 represents equivalent added mass when L 1 = 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.

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 f 3 = O(ε 3 ), 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).

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.

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 L 1 = 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, a 2 = 0 becomes unstable. Then, substituting a 2 = 0 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 L 1 = 75 µm; (3) the added mass can decrease the minimum critical frequency of anti-symmetric mode when L 1 = 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.

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)

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 M < 0 results in the subcritical bifurcation. With the increase of AC voltage, the second order amplitude suddenly appeared. Likewise, the case M > 0 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 L 1 = 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 L 1 = 75 µm.

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

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.
Firstly, L 1 = 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 m = 1 × 10 −6 µg, the supercritical bifurcation voltage V ac1 is obtained. With the increase of the driving voltage, the second order amplitude appears gently. When m = 5 × 10 −6 µg, the bifurcation voltage increases and the supercritical bifurcation is transformed into subcritical bifurcation. The subcritical bifurcation voltage V ac2 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. 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.  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 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.
Similarly, L 1 = 75 µm is considered. From Equation (23), the added mass has no effect on bifurcation type when L 1 = 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.

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 L 1 = 50 µm and L 1 = 75 µm, respectively. Figures 7-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 Substituting Equation (20) into Equation (18) yields where a 3 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, L 1 = 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 δ = −0.55. 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, Figures 14 and 15 show force-amplitude curve and amplitude-frequency response curve in the case of L 1 = 75 µm. With Equations (26) and (27), we obtain mass identification results, as shown in Tables 4 and 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 L 1 = 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.

Number
The true mass m (10 −6 µg)  Table 3. Three groups of mass detection results obtained by the amplitude-frequency response curve

Number
The true mass   Table 4. Three groups of mass detection results obtained by the force-amplitude curve when L 1 = 75 µm, δ = −0.55.

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. Figures 16 and 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 L 1 = 75 µm is improved and the sensitivity of the mass sensor in the case of L 1 = 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.
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.

Conclusion
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

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 L 1 = 75 µm. On the contrary, with the increase of DC voltage, the sensitivity of the mass sensor is suppressed when L 1 = 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.