Research Article  Open Access
Research on Nonlinear Dynamics of Drive Mode in ZAxis Silicon Microgyroscope
Abstract
A novel drive method is proposed to increase the vibration amplitude of drive mode of silicon vibratory gyroscopes based on the nonlinear dynamics research made in this paper. A nonlinear dynamic equation, considering the nonlinear mechanicalelastic and electrostatic force, was firstly deduced according to the simplified structure schematic of drive mode. Then, experiments were made to characterize the nonlinear amplitudefrequency response characteristic as well as its dependence on drive voltage. Besides, by solving the nondimensional nonlinear dynamics using Matlab and Simulink, the influence of initial status of drive mode, that is, the displacement and the velocity, on the stable vibration amplitude was presented. The amplitudefrequency characteristic was also analyzed by simulation. The result of simulation shows good agreement with that of experiments. When nonlinearity is present in the amplitudefrequency characteristic, there is a wide flat region with larger amplitude than the peak amplitude obtained in the linear characteristic. Finally, a novel drive method was presented to excite the nonlinear vibration of large amplitude in the downwardssweeping characteristic curve. It is numerically proved that this merit can be utilized to enlarge the drive mode’s vibration amplitude so as to enhance silicon gyroscopes’ performance.
1. Introduction
As part of a new generation of inertial sensors, silicon microgyroscopes, used for measurement of angular rate or angle of a rotating carrier, have a lot of merits, such as low power consumption, small volume and weight, high shock resistivity, and suitability for batch fabrication [1–5]. However, the performance of silicon microgyroscopes is lower than conventional mechanical, rotating wheel, fiberoptical, and ring laser gyroscopes [6]. Thus, it is necessary to improve the performance so as to broaden their application market.
Large vibration amplitude of drive mode is one of the conditions for silicon gyroscopes to achieve high performance [7, 8]. Therefore, capacitive actuators of more combs and smaller gaps are designed in microgyroscopes so as to increase the vibration amplitude when the same drive voltage is applied. However, when gaps become small and the vibration amplitude gets large, the nonlinear terms of the force applied on proof mass can become sufficiently high to bend the amplitudefrequency characteristic curve. For the traditional drive methodology, drive mode works at the natural frequency in the linear amplitudefrequency characteristic curve. However, when the amplitudefrequency characteristic curve is nonlinear, the drive methodology for the drive mode of silicon gyroscopes is not yet noticed.
Nonlinear dynamics of MEMS resonators has been studied analytically, numerically, and experimentally by some institutions [9–14]. In these literatures, dynamic models are constructed to match the theoretical characteristic curves with the experimental phenomena of real resonator samples [11, 12]. The models can be solved by different methods, such as the theoretical method of multiple scales [10], the semianalytical harmonic method [9], the numerical integration method [9–11], and the electrical equivalentcircuit modelling method [14]. It is known that the nonlinear phenomena refer to the spring softening and hardening phenomena in the amplitudefrequency curves. In capacitive MEMS resonators, the phenomena mainly result from the nonlinear mechanicalelastic reaction force when the mechanical support system is out of its linear work range, as well as the nonlinear electrostatic actuation force produced when the squeeze gaps of capacitors are small. As the drive mode of gyroscope is considered as a resonator, some more research is made in this paper on the nonlinear characteristics of the gyroscope prototype designed by our research team. Based on the nonlinear amplitudefrequency characteristic, a novel drive method is provided to make drive mode vibrate at large amplitude.
The outline of the paper is as follows. First, in Section 2, the nonlinear dynamics is derived based on the simplified schematic of the drive mode. Then, in Section 3, the experimental amplitudefrequency curves are obtained by using a detect circuit and sweeping the actuation frequency increasingly and decreasingly, separately. Besides, the nonlinear phenomena are also illustrated by solving the dynamics numerically using the Matlab and Simulink model in Section 4. Finally, in Section 5, a novel drive method is proposed for actuating the vibration of large amplitude in the nonlinear amplitudefrequency characteristic curve.
2. Dynamics of Drive Mode
2.1. Structure Schematic for Drive Mode
In the axis silicon microgyroscope designed by our research team, the proof mass is supported by straight beams, as well as Ushaped beams, and driven and detected by capacitors of comb fingers [15–17]. The dimensions of the gyroscope sample, fabricated by SOG process, are around 8 mm × 8 mm × 60 μm and the SEM photos of beams and capacitors are displayed in Figure 1. Figure 1 shows the simplified structure schematic for the drive mode of the microgyroscope, consisting of the proof mass, comb actuators, sense combs, and supporting beams. As seen in the graph, Beam 1 is a fixedguided straight beam, and Beam 2 and Beam 3 are the parallel and the serial combination of Beam 1, respectively. The proof mass is displaced by the electrostatic force produced when differential voltage is applied on drive electrodes. Since system’s stiffness in the sense direction (axis) is much larger than that in the drive direction (axis), the motion of drive mode is simplified as 1DOF translation in the axis, as displayed in (1). is the displacement of the proof mass along the drive direction, and are the electrostatic force and the mechanicalelastic reaction force, and is the damping force with being the damping coefficient:
2.2. MechanicalElastic Reaction Force
In Figure 1, as the displacement of the proof mass in the drive direction is , the deformation of the straight beam in Beams 3 is . Since the equivalent stiffness of Beam 2 is twice the stiffness of Beam 1, the assumption is made that the deformation in Beam 1 is twice that in Beam 2. Thus, the deformation of the straight beam in Beams 1 and 2 is and , respectively. The elastic reaction force applied on proof mass by all Beams 2 is approximated through a cubic polynomial relation, as displayed in (2) [10], where is the elastic modulus of silicon and , , and are the width, depth, and length of a straight beam. Similarly, the reaction force produced by all Beams 3 is shown in (3). Thus, the total mechanicalelastic reaction force applied on the proof mass in the drive direction is deduced in (4), where is defined as the linear mechanical stiffness and as the quadratic mechanical stiffness coefficient. Consider
2.3. Electrostatic Force
Two drive electrodes in Figure 1 are applied with differential voltages and relative to the proof mass, in which is the DC bias, is the amplitude of AC voltage, and is the actuation frequency. N groups of comb fingers are attached to one drive stator. Figure 2 illustrates one group of comb fingers, made up of one moving comb and two fixed combs. As explained in Section 2.1, the moving combs are assumed to vibrate only in the axis. Then, the electrostatic force applied on the proof mass along the axis is derived in (5) [5]. is generated by the overlapped areas with a slide gap of and is generated by overlapped areas with a squeeze gap of . It is easily noticed that is independent of while has a nonlinear relationship with it. Thus, the nonlinear electrostatic force applied on proof mass is relevant to the squeeze gap .
2.4. Nonlinear Dynamics
Substituting (4) and (5) into (1), the nonlinear dynamics for drive mode is displayed in (6), which can be simplified as (7) when is much larger than . Based on (7), the natural frequency of drive mode is approximated in (8), where is the linear electroelastic stiffness and is the quadratic electroelastic stiffness coefficient. It is easily seen that the natural frequency increases with the mechanicalelastic stiffness and decreases with the electroelastic stiffness. Besides, the mechanicalelastic and electroelastic stiffness both cause the nonlinearity between the natural frequency and the displacement. The nondimensional differential dynamics, as seen in (9), is obtained from (6) by the replacement of , , , and . The nondimensional dynamics is used in the following numeric analysis:
3. Experimental Characterization
3.1. Experimental Setup
As presented in Figure 3, the carrier voltage is applied on the anchor, that is, the moving combs, and two differential sense electrodes are connected with their own preamplifiers. The current, and i_{s2}, is integrated and the output voltage, and , passes through the differential amplifier. The process is shown as below: Since the squeeze gap in the sense combs is much larger than that in the drive combs, the sense capacitance between the squeeze faces is ignored and the approximate sense capacitance change is proportional to the displacement ; that is, . Thus, if , the output voltage is Then, is filtered by a highpass filter and displayed by a spectrum analyzer on which the vibration information of drive mode is at twice the actuation frequency.
The practical experimental setup is illustrated in Figure 4. Through the probe station, a packaged gyroscope sample is connected with the detect circuit. During testing, the actuation frequency of the carrier voltage is consecutively tuned by a step of 0.1 Hz increasingly or decreasingly between 3310 Hz and 3330 Hz, and the corresponding amplitude information is recorded by acquiring the value at twice the actuation frequency on the analyzer screen.
3.2. Experimental Results
When V, V, and sweeps up slowly from 3310 Hz to 3330 Hz, the frequency characteristic curve of the drive mode is fitted in red, as shown in Figure 5, using the amplitude data acquired at each frequency step. Similarly, the characteristic curve, while frequency sweeps down, is fitted in black. It is noticed that when the actuation frequency sweeps down, the natural frequency is around 3313 Hz, much lower than that obtained during upwards sweeping, that is, around 3319 Hz. Both frequency response curves bend towards lower frequency.
Experiments are also made on the influence of applied voltage on the nonlinearity of frequency response. In Figure 6, when increases from 5 V to 6 V, the natural frequency of upwardsweeping reduces from 3319 Hz to 3316 Hz and that of downwardsweeping decreases from 3313 Hz to 3308 Hz. Thus, increasing makes the spring softening phenomena more evident. Similarly, as seen in Figure 7, increasing also aggravates the nonlinear softening characteristic.
(a)
(b)
(a)
(b)
4. Numerical Characterization
4.1. Simulink Model
Matlab and Simulink are used to determine the motion of the proof mass as well as its frequency response characteristic. Parameters, which are given in Table 1, are firstly defined in the Matlab code and then transferred to the Simulink model for simulation. Figure 8 gives the Simulink model, in which an SFunction module is used to construct (9), differential voltage is applied, and the response is monitored by oscilloscopes. By simulating the model, the stable solution of displacement and velocity is obtained. Besides, the amplitudefrequency characteristic curve is determined by fitting the stable vibration amplitude acquired at each actuation frequency. This model is also used to determine and compare the response under different initial status and voltage loads. Here, the actuation frequency is swept in two ways. Within the frequency range of interest, one sweeping method is to increase frequency consecutively by a certain step and the other is to decrease it by the same step. It must be pointed out that, at each frequency sweeping step, the initial status of the displacement and velocity as well as the phase of applied voltage is maintained equal to the value at the finish time of the previous frequency step.

4.2. Influence of Initial Status
When drive voltage is applied with V, V, and = 3315 Hz, the phase plane portrait of (9) is illustrated in Figure 9, in which the arrow tells the progressing trend of the nondimensional displacement and velocity along with time [18]. Besides, three curves in the graph refer to the equivalent transient response trajectories under different initial status. It is observed that the nondimensional stable amplitude of displacement is around 0.1 when initial status is , while it is around 0.85 when . In addition, the trajectory is divergent when . It is easily concluded that the stable motion of drive mode is dependent on the initial motion status of the proof mass when voltage is applied as here.
4.3. AmplitudeFrequency Response
The Simulink model in Figure 8 and the parameters in Table 1 are reused to obtain the characteristic curve of frequency response while voltage is applied as V, V, and the actuation frequency is increased and decreased consecutively by a step size of 0.5 Hz between 3305 Hz and 3330 Hz. Besides, the initial status at the start frequency is . The corresponding stable vibration amplitude at each frequency step is collected to fit characteristic curves.
Figure 10 gives the simulated characteristic curves, in which the pink solid line and the blue dashed line refer to the responses while frequency sweeps up and down, respectively. It is noticed that two natural frequencies and appear in frequency response curves. The curves both bend towards lower frequency, manifesting high nonlinearity. During sweeping up, when the actuation frequency gets near , the amplitude rushes to the peak and then descends gradually. Similarly, during sweeping down, the amplitude ascends slowly to the peak and then drops dramatically at the frequency . It is also seen that two natural frequencies are different with , and their corresponding peak amplitude is not the same with . Besides, when Hz, the obtained downwardssweeping amplitude is over 7 times larger than the upwardssweeping amplitude and even larger than . In the downwardssweeping curve, there is a flat region with relatively large amplitude before is reached. Thus, the nonlinear downwardssweeping amplitudefrequency characteristic curve is very attractive for performance enhancement of silicon gyroscopes.
By using the same Simulink model, some investigation is made on the dependence of the amplitudefrequency characteristic on applied voltage. Simulations are made by changing parameters and separately. The simulated results are shown in Figures 11 and 12. When and increase, the two natural frequencies obtained by different sweeping methods both descend but their corresponding peak amplitude ascends. Besides, the decrement of is much larger than that of. It can be concluded that the frequency response nonlinearity can be adjusted by and and the flat region of the downwardssweeping curve grows larger when nonlinearity is intensified.
(a)
(b)
(a)
(b)
5. Excitation of Nonlinear Vibration
According to Figure 9, when the initial status of vibration is and differential voltages and are applied on drive combs, the vibration of the drive mode will be stabilized at relatively large amplitude, which is especially attractive for improving the scale factor of silicon gyroscopes. Therefore, Figure 13 illustrates how to excite the nonlinear vibration with larger amplitude. A triangular impulse is firstly loaded to start the motion of drive mode and the vibration displacement and velocity are recorded. As shown in Figure 14, the impulse intensity is adjusted to make the motion status appear at some time. Thus, the differential drive voltage is loaded at this time. The nonlinear response, as illustrated in Figure 15, is obtained by simulating the model in this way. It is noticed that when changes consecutively from 3315 Hz and 3314 Hz to 3316 Hz, the stable displacement amplitude varies from 0.848 and 0.867 to 0.823 and the phase relative to the normalized sine voltage is maintained around −108.6 degrees.
According to the numerical and experimental characteristics, the nonlinear vibration with larger amplitude can be actuated by decreasing the drive frequency consecutively with a certain frequency step and time interval. Thus, the nonlinear vibration in the downward sweeping curve can be excited by monitoring or gradually accumulating the status of vibration. Finally, the novel drive methodology makes drive mode vibrate at large amplitude and it is robust with no drops of amplitude and phase when the drive frequency fluctuates.
6. Discussion
When differential voltages and are applied on drive stators with V and V, the amplitudefrequency characteristic curves obtained by experiments and simulations are in good agreement. It is displayed that the frequency characteristic curves bend towards lower frequency with the natural frequency much smaller than the initial value. Besides, the peak vibration amplitude and the natural frequency, obtained by different sweeping methods, are not the same.
According to (8), it is known that and lead to reduction of the natural frequency while raises it. In addition, makes the frequency response curves bend towards lower frequency while makes them bend towards higher frequency. Thus, the spring softening phenomena in Figures 5 and 10 imply the relationship between the nonlinear stiffness coefficients in (8); that is, .
Based on the definition of and , increasing or decreasing makes electroelastic stiffness larger, leading to more evident spring softening phenomena. Besides, according to (8), the nonlinear stiffness is also increased by displacement enlargement. Thus, increasing and , which makes the vibration amplitude of the proof mass larger, aggravates spring softening phenomena. The same influence of applied voltage on the amplitudefrequency characteristic is indicated by the experimental results in Figures 6 and 7 as well as the numerical analysis in Figures 11 and 12.
In the nonlinear amplitudefrequency characteristic curves, though the vibration amplitude and the natural frequency are multivalued, the flat region of downwardssweeping curves with extremely large amplitude is useful for improving gyroscope’s performance, such as sensitivity and robustness.
It is numerically proved that the vibration of the drive mode ends with different stable amplitude under different initial motion status, when applied voltage is large enough. Thus, the nonlinear vibration with larger amplitude can be excited by loading drive voltage at the specific status or decreasing the drive frequency gradually by a certain frequency and time step. The obtained amplitude in the downward sweeping curve is much larger than the amplitude in the linear frequency response curve. Besides, the amplitude and the phase of vibration are less influenced by the fluctuation of the actuation frequency. It can be concluded that this novel drive methodology is beneficial for performance improvement of silicon vibratory gyroscopes.
7. Conclusion
In this paper, the nonlinear frequency softening characteristic of the drive mode in the axis silicon microgyroscope designed by our research team is determined by analytical, experimental, and numerical methods, results of which show good agreement qualitatively. It is illustrated that the spring softening phenomena of the amplitudefrequency characteristic appear when applied voltage and the vibration amplitude are large enough. Besides, it is also proved that applied voltage plays a significant role in the nonlinear vibration of the gyroscope, so the nonlinearity can be tuned by drive voltage. Moreover, under nonlinear vibration, the wide flat region with large amplitude in the downwardssweeping amplitudefrequency characteristic curve is attractive for improving gyroscope’s performance. The nonlinear vibration of large amplitude can be excited by applying drive voltage at a certain motion status or gradually decreasing the actuation frequency. Besides, the large amplitude and the phase of vibration are maintained when the actuation frequency fluctuates. Thus, for silicon microgyroscopes, it can be a promising drive method to make the drive mode vibrate at the large amplitude in the flat region of the nonlinear amplitudefrequency characteristic curve.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work is supported by preresearch fund Project of 9140A09011313JW06119.
References
 S. Wang, L. Huang, and B. Yang, Weiguanxing yibiao yu weixitong, Weapon Industry Press, Beijing, China, 2011.
 N. Yazdi, F. Ayazi, and K. Najafi, “Micromachined inertial sensors,” Proceedings of the IEEE, vol. 86, no. 8, pp. 1640–1658, 1998. View at: Publisher Site  Google Scholar
 A. M. Shkel, C. Acar, and C. Painter, “Two types of micromachined vibratory gyroscopes,” in Proceedings of the 4th IEEE Conference on Sensors, pp. 531–536, Irvine, Calif, USA, November 2005. View at: Publisher Site  Google Scholar
 A. M. Shkel, “Micromachined gyroscopes: challenges, design solutions, and opportunities,” in Smart Structures and Materials: Smart Electronics and MEMS, vol. 4334 of Proceedings of SPIE, pp. 74–85, Newport Beach, Calif, USA, August 2001. View at: Publisher Site  Google Scholar
 W. A. Clark, Micromachined vibratory rate gyroscopes [Ph.D. thesis in EngineeringElectrical Engineering and Computer Science], University of California, Berkeley, Calif, USA, 1997.
 W. Geiger, J. Bartholomeyczik, U. Breng et al., “MEMS IMU for AHRS applications,” in Proceedings of the IEEE/ION Position, Location and Navigation Symposium (PLANS '08), pp. 225–231, May 2008. View at: Publisher Site  Google Scholar
 A. Sharma, F. M. Zaman, B. V. Amini, and F. Ayazi, “A highQ inplane SOI tuning fork gyroscope,” Proceedings of IEEE Sensors, vol. 1, pp. 467–470, 2004. View at: Publisher Site  Google Scholar
 A. Sharma, M. F. Zaman, and F. Ayazi, “A sub0.2°/hr bias drift micromechanical silicon gyroscope with automatic CMOS modematching,” IEEE Journal of SolidState Circuits, vol. 44, no. 5, pp. 1593–1608, 2009. View at: Publisher Site  Google Scholar
 F. Braghin, F. Resta, E. Leo, and G. Spinola, “Nonlinear dynamics of vibrating MEMS,” Sensors and Actuators, A: Physical, vol. 134, no. 1, pp. 98–108, 2007. View at: Publisher Site  Google Scholar
 A. M. Elshurafa, K. Khirallah, H. H. Tawfik, A. Emira, A. K. S. Abdel Aziz, and S. M. Sedky, “Nonlinear dynamics of spring softening and hardening in foldedmems comb drive resonators,” Journal of Microelectromechanical Systems, vol. 20, no. 4, Article ID 5772897, pp. 943–958, 2011. View at: Publisher Site  Google Scholar
 R. M. C. Mestrom, R. H. B. Fey, K. L. Phan, and H. Nijmeijer, “Simulations and experiments of hardening and softening resonances in a clampedclamped beam MEMS resonator,” Sensors and Actuators, A: Physical, vol. 162, no. 2, pp. 225–234, 2010. View at: Publisher Site  Google Scholar
 R. M. C. Mestrom, R. H. B. Fey, J. T. M. van Beek, K. L. Phan, and H. Nijmeijer, “Modelling the dynamics of a MEMS resonator: simulations and experiments,” Sensors and Actuators A: Physical, vol. 142, no. 1, pp. 306–315, 2008. View at: Publisher Site  Google Scholar
 M. Palaniapan and L. Khine, “Nonlinear behavior of SOI freefree micromechanical beam resonator,” Sensors and Actuators, A: Physical, vol. 142, no. 1, pp. 203–210, 2008. View at: Publisher Site  Google Scholar
 T. Veijola and T. Mattila, “Modeling of nonlinear micromechanical resonators and their simulation with the harmonicbalance method,” International Journal of RF and Microwave ComputerAided Engineering, vol. 11, no. 5, pp. 310–321, 2001. View at: Publisher Site  Google Scholar
 B. Yang, B. Zhou, S. Wang, L. Huang, and Y. Yin, “A quadrature error and offset error suppression circuitry for Silicon MicroGyroscope,” in Proceedings of the 3rd IEEE International Conference on Nano/Micro Engineered and Molecular Systems (NEMS '08), pp. 422–426, January 2008. View at: Publisher Site  Google Scholar
 Y. Yin, S. Wang, C. Wang, and P. Sheng, “Structural scheme design and simulation of structuredecoupled dualmass MEMS gyroscope,” Journal of Southeast University, vol. 38, no. 5, pp. 918–922, 2008. View at: Google Scholar
 Y. Ni, H. Li, L. Huang, and L. Zhao, “Design and test of quadrature correction structure for silicon microgyroscope,” Journal of Southeast University: Natural Science Edition, vol. 43, pp. 1227–1231, 2013. View at: Google Scholar
 W. Zhang, G. Meng, and D. Chen, “Stability, nonlinearity and reliability of electrostatically actuated MEMS devices,” Sensors, vol. 7, no. 5, pp. 760–796, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Lu Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.