Next Article in Journal
Deposition of Very-Low-Hydrogen-Containing Silicon at a Low Temperature Using Very-High-Frequency (162 MHz) SiH4 Plasma
Next Article in Special Issue
Self-Calibration Technique with Lightweight Algorithm for Thermal Drift Compensation in MEMS Accelerometers
Previous Article in Journal
Low Cost, Easily-Assembled Centrifugal Buoyancy-Based Emulsification and Digital PCR
Previous Article in Special Issue
Biaxial Angular Acceleration Sensor with Rotational-Symmetric Spiral Channels and MEMS Piezoresistive Cantilevers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Optimization of a Novel MEMS Tuning Fork Gyroscope Microstructure

1
School of Mechanical Science and Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
2
Quality Engineering Center, China Aero Polytechnology Establishment, Beijing 100028, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(2), 172; https://doi.org/10.3390/mi13020172
Submission received: 14 December 2021 / Revised: 13 January 2022 / Accepted: 23 January 2022 / Published: 24 January 2022
(This article belongs to the Special Issue MEMS Mechanical Sensing Device)

Abstract

:
This paper presents the design and optimization of a novel MEMS tuning fork gyroscope microstructure. In order to improve the mechanical sensitivity of the gyroscope, much research has been carried out in areas such as mode matching, improving the quality factor, etc. This paper focuses on the analysis of mode shape, and effectively optimizes the decoupling structure and size of the gyroscope. In terms of structural design, the vibration performance of the proposed structure was compared with other typical structures. It was found that slotting in the middle of the base improved the transmission efficiency of Coriolis vibration, and opening arc slots between the tines reduced the working modal order and frequency. In terms of size optimization, the Taguchi method was used to optimize the relevant feature sizes of the gyroscope. Compared with the initial structure, the transmission efficiency of Coriolis vibration of the optimized gyroscope was improved by about 18%, and the working modal frequency was reduced by about 2.7 kHz. Improvement of these two indicators will further improve the mechanical sensitivity of the gyroscope.

1. Introduction

The application of traditional gyroscopes such as mechanical rotor gyroscopes [1,2], laser gyroscopes [3], fiber optic gyroscopes [4], etc., has been severely restricted due to their large size, high cost and unsuitability for mass production. Benefiting from the small size, low cost, low power consumption and batch fabrication, MEMS gyroscopes have been widely used in the fields of automotive safety, mobile robots, consumer electronics, aerospace navigation, military weapons, etc., over the past few decades [5,6,7,8,9]. According to different vibration structures, MEMS gyroscopes can be divided into gimbal gyroscopes, tuning fork gyroscopes, vibrating ring gyroscopes, solid-state gyroscopes, etc. [10]. Among them, the tuning fork gyroscope is mostly preferred due to its superior performance with regard to its flat structure for mass production, common-mode rejection, and low power consumption [11,12,13]. The sensitivity of the tuning fork gyroscope is one of the main challenges in terms of superior performance [14].
There are generally two ways to improve the sensitivity of tuning fork gyroscopes: peripheral circuit gain and improving mechanical sensitivity. Much research has been conducted on the former [15,16,17,18], however, it is not enough to enhance sensitivity only by circuit gain. Without sufficient mechanical sensitivity, improvements in circuit gain are not very useful. Therefore, many scholars have carried out related research, such as performing modal matching through various techniques [19,20,21,22], improving the quality factor [23,24], and suppressing the zero-rate output drift, etc. For example, Behbahan et al. proposed a systematic post fabrication technique to reduce the modal frequency differences between the n = 2 and n = 3 modes in an axisymmetric resonator [19]. Li et al. presented a leverage mechanism to improve the mechanical sensitivity of a MEMS gyroscope [14]. Xu et al. investigated a tuning fork gyroscope with a polygon-shaped vibration beam, and improved the sensitivity by optimizing the spindle azimuth [25]. Nguyen et al. proposed a z-axis tuning fork gyroscope with a freestanding and anti-phase controlled architecture, which could improve the performance by limiting the influence of air damping [11]. Wu et al. reported a z-axis quartz tuning fork gyroscope based on shear stress detection, and designed a symmetric tapered beam to increase the sensitivity [26]. Wang et al. designed the multiple-beam tuning-fork structure to achieve high quality factors, promoting the sensitivity [13]. To suppress the zero-rate output drift of MEMS tuning fork gyroscopes, Xu et al. proposed a temperature-induced stress release method [27].
This study focuses on the mode shape analysis of the gyroscope to achieve the purpose of optimizing the mechanical sensitivity. For a tuning fork gyroscope to work properly, it includes a drive structure to provide drive motion and a sensing structure to detect Coriolis vibration signals. However, the kinetic coupling between the two structures may interfere with the detection signal and reduce the mechanical sensitivity [28]. In this paper, a novel tuning fork gyroscope decoupling microstructure is proposed, which produces a relative separation between driving motion and sensing motion. Furthermore, by analyzing the mode shape of the model, the structure and size of the gyroscope are optimized, thereby improving the mechanical sensitivity. The material of the gyroscope is quartz, benefiting from good piezoelectric effect, insulation, mechanical properties and high quality factor. Details of the gyroscope structure and analysis of mode shape are provided in Section 2. Section 3 compares the advantages of the proposed structure over other typical structures and optimizes the relevant feature sizes of the gyroscope with the Taguchi method. Section 4 contains the discussion. Section 5 concludes this paper with a summary.

2. Materials and Methods

2.1. Architectural Design

A novel MEMS tuning fork gyroscope microstructure is proposed, as shown in Figure 1. The gyroscope includes two drive tines, two sense tines, and a base with a support beam. The drive tines are used to provide driving motion, while the sense tines are used to detect the Coriolis signal. The base transmits the Coriolis force generated on the drive tines to the sense tines to eliminate the coupling between the driving motion and the sensing motion.
When the gyroscope is in driving mode, as shown in Figure 2a (displacement field distribution; the unit of the analysis result is meters), the drive tines will move in opposite directions in same plane. Once there is an angular velocity input in the Y direction, the drive tines will generate Coriolis vibration in the Z direction. The Coriolis force on the two drive tines is not in the same plane, thus forming a coupling. This coupling will cause the torsion of the base, thereby transmitting the Coriolis vibration on the drive tines to the sense tines, as shown in Figure 2b. The drive tines and sense tines of the gyroscope are relatively independent, where the driving motion and sensing motion do not interfere with each other. Thus, the decoupling microstructure relatively solves the problem of interdigital motion coupling and capacitive parasitic coupling, which is conducive to signal detection.
The conversion of vibration and electrical signals in the tuning fork gyroscope microstructure is based on the piezoelectric effect, thereby quartz crystals with this property are selected as a material for manufacturing. After stimulation by the driving voltage and angular velocity signal, the gyroscope outputs an electrical signal to the external circuit, and vibration characteristics are obtained by analyzing the electrical signal. The overall electrode arrangement of the gyroscope is shown in Figure 3a. Figure 3b,c represent the electrode directions of drive tine 1 and drive tine 2, and Figure 3d,e represent the electrode directions of sense tine 1 and sense tine 2, respectively. To simplify the model, since the actual electrode coating is very thin, its influence is ignored in the subsequent mode shape analysis of the gyroscope.

2.2. Mode Shape

The actual vibration of the system can be understood as the superposition of the various mode shapes, according to the relevant knowledge of mechanical vibration. The system response to the initial stimulus follows Formula (1). ® general response of the system to any stimulus follows Formula (2) [29].
{ q ( t ) } = r = 1 n ( { u ( r ) } T [ m ] { q 0 } cos ω r t + { u ( r ) } T [ m ] { q 0 } 1 ω r sin ω r t ) { u ( r ) }
{ q ( t ) } = r = 1 n { u ( r ) } N 0 r / ω r 2 [ 1 ( ω / ω r ) 2 ] 2 + ( 2 ξ r ω / ω r 2 ) 2 sin ( ω t φ r )
where {q(t)} is the system response, {u(r)} is the r-th modal vector, [m] is the mass matrix, {q0} is initial displacement, {q0} is initial velocity, ωr is the r-th natural frequency, N0r is excitation amplitude, ω is excitation frequency, ξr is damping ratio, and φr is phase of the system response.
Formula (1) shows that in the total vibration of the system, the initial motion conditions of the system determine the degree to which the r-th mode is excited. If {q0} approximates the s-th modal vector, the system will vibrate almost according to the s-th mode. Formula (2) shows that the amplitude of the system response is related to the excitation frequency. Obviously, when we apply excitation to make the gyroscope work in the driving mode while the excitation frequency is equal to the driving modal frequency, the gyroscope will resonate in the driving mode. At this time, it is uncertain whether the sensing vibration aroused by the angular velocity vibrates in the manner of the sensing mode.
To verify the conjecture, three sets of analyses have been conducted, and the analysis models of the gyroscope are shown in Figure 4. Dimensions of all models are in millimeters, and are analyzed by Ansys software. Moreover, the geometric structure of the designed gyroscope is simple and symmetrical, so it is more suitable to select hexahedral elements for meshing. Taking type A as an example, its mesh model is shown in Figure 5. The node A coincides with the endpoint of the right drive tine, and the node B coincides with the endpoint of the right sense tine. Their specific locations are indicated in Figure 5. The analysis index is the vibration state of node A and node B in the actual sensing vibration and sensing mode, and analyzing the sensing vibration with harmonic analysis. In addition, a zero-displacement constraint is applied to the bottom of the gyroscope, and a set of sinusoidal forces with equal amplitude and opposite phase are applied to the front and back surface of the drive tines, respectively. Moreover, the amplitude of the load F is 0.0001 N, while the frequency of it is 200 Hz less than the sensing modal frequency. The application area is 1/3 of the upper part of the drive tine. The loading diagram is shown in Figure 6a. In the sensing mode, a zero-displacement constraint of the bottom is the only load, which is shown in Figure 6b.

3. Results and Analysis

The actual sensing vibration state is in line with the sensing mode. As shown in Table 1, in the harmonic analysis, the amplitude ratio and phase difference of the displacement between node A and node B are both approximate to the sensing mode. This shows that although the vibration frequency is different from the sensing modal frequency by 200 Hz, when the drive tines vibrate in the Z direction oppositely affected by Coriolis force, the sense tines will also vibrate in the Z direction correspondingly, and the vibration state is almost consistent with the sensing mode. According to Formula (1), at this time, the Coriolis vibration state of the gyroscope is almost completely determined by the sensing mode. Therefore, it is credible to evaluate the transmission of Coriolis vibration by use of the sensing mode, which can greatly reduce the calculation required.
The gyroscope generally works in resonance, and it is known from Table 1 that the displacement of the drive tines affected by Coriolis force has little connection with the shape of the base. Hence, the key to improving the mechanical sensitivity of the gyroscope is in improving the transmission efficiency of the Coriolis signal from the drive tines to the sense tines as much as possible. The transmission efficiency is denoted as ŋ and is defined by the formula given below.
η = Z B Z A
where ZB is the displacement of node B in the sensing mode, and ZA is the displacement of node A.

3.1. Structural Analysis

The mechanical sensitivity of the gyroscope is denoted as S and can be calculated according to Formula (4), where Y 0 is sense tine amplitude, Ω is angular velocity, η is transmission efficiency of the Coriolis vibration, F is driving force, and the driving direction quality factor and resonant frequency are represented by Q x and ω x , respectively. In addition, m is mass of the drive tine, and Q y and   ω y represent the quality factor and resonant frequency of driving direction, respectively.
S = Y 0 Ω = 2 η F Q x m ω x ω y 2 1 ( 1 ( ω x ω y ) 2 ) 2 + ( ω x ω y Q y ) 2
According to the calculation formula of mechanical sensitivity, the gyroscope microstructure can be optimized from the following aspects:
(1)
The transmission efficiency of Coriolis vibration ŋ should be as high as possible;
(2)
The working modal (driving mode and sensing mode) order and frequency should be as low as possible, which can achieve greater response amplitude in the structural vibration. When the excitation frequency is close to the r-th natural frequency ωr, the system response {q(t)} is approximately expressed by the following formula.
{ q ( t ) } { u ( r ) } N 0 r 2 ξ r ω r sin ( ω r t π 2 )
where {u(r)} is the r-th modal vector, N0r is excitation amplitude, and ξr is damping ratio.
It is known from Formula (5) that the lower the working modal frequency is, the greater the response amplitude is, and therefore, the higher is the sensitivity of the gyroscope;
(3)
The driving mode and the sensing mode should be adjacent to reduce the effects of interfering modes.
In addition to the three structures shown in Figure 4, six typical structures are analyzed, which are shown in Figure 7. The mode frequency and Coriolis vibration transmission efficiency results of all models are shown in Table 2.
As seen from Figure 8, slot structure types have different effects on the transmission efficiency. Based on type D, slotting on the sides of the base will reduce the transmission efficiency, such as in type C. On the contrary, slotting in the middle of the base will improve the transmission efficiency, such as in types E, F and G. According to Table 1, opening arc slots between the tines can reduce the working modal order and frequency, such as in type I. Therefore, combining the two advantageous features—slotting in the middle of the base and opening arc slots between the tines—can greatly improve the mechanical sensitivity of the gyroscope. The proposed type A not only has the above characteristics, but additionally, its base protrudes slightly outward, so that the sense tine is also symmetrical in a local area, thereby obtaining better vibration performance.
After modal analysis of different models, it was found that type A, with high symmetry, had the best vibration performance. Its transmission efficiency was the highest, and the working mode order and frequency were also the lowest. In addition, the driving and sensing modes were adjacent to each other. In summation, type A perfectly meets the three preset performance requirements.

3.2. Size Optimization with Taguchi Method

The Taguchi method is a reliable design method that guides the optimal design of products in complex environments. In this section, the Taguchi method is used to optimize the relevant feature sizes of the gyroscope. The feature size codes of the gyroscope are shown in Figure 9. The height of tines h2, drive-sense interval d1, drive–drive interval d2, width of the support beam w1, and height of the groove h3, are the signal factors. The noise factor is the mesh size of the finite element model. Parameters and their levels are given in Table 3 and Table 4.
The best experimental result is evaluated by the calculation of signal to noise ratio (SNR) in the Taguchi method. The larger the SNR is, the smaller the quality loss and the better the product quality. To obtain the optimum result, three basic categories which are smaller-the-better, larger-the-better and nominal-the-best characteristics are calculated by the formulas given below [30].
(a)
Smaller-the-better characteristic:
SNR = 10 lg ( 1 n i = 1 n y i 2 )
(b)
Larger-the-better characteristic:
SNR = 10 lg ( 1 n i = 1 n 1 y i 2 )
(c)
Nominal-the-best characteristic:
SNR = 10 lg ( 1 n i = 1 n ( y i m ) 2 )
where yi is measured values, and n and m are number of experiments and nominal value respectively.
The orthogonal experiment table was created according to Table 3 and Table 4, and the corresponding finite element model was established to analyze according to different combinations of each factor level. The observed experiment results include the following items:
(1)
The transmission efficiency of Coriolis vibration which meets the larger-the-better characteristic;
(2)
Frequency difference between the driving mode and sensing mode. When the frequency difference between the two is large, the working bandwidth of the gyroscope meets Formula (9) [31]:
B W = 0.54 Δ F
where ΔF is the frequency difference between the driving mode and sensing mode.
Given that the applications of tactical and inertial gyroscope require bandwidths of around 100 Hz [9], the frequency difference ΔF should be maintained at about 200 Hz. For the convenience of calculation, this paper defines the difference between ΔF and 200 Hz as a new variable Δf, which meets the smaller-the-better characteristic and follows the Formula (10) given below:
Δ f = | Δ F 200 |
(3)
The driving modal frequency which meets the smaller-the-better characteristic;
(4)
The drive coupling coefficient which meets the smaller-the-better characteristic. Similar to Section 3.2, the drive coupling coefficient is denoted as λ and defined by the Formula (11) given below.
λ = X B X A
where XB is the displacement of node B in the driving mode, and XA is the displacement of node A.
Since it is preferable for the undesired coupling motion between the drive tines and the sense tines to be separated, it is better to keep the sense tines stationary when the drive tines are vibrating in the opposite direction, and, therefore, the drive coupling coefficient needs to be as low as possible.
According to the different quality characteristics of each result index, the corresponding SNR calculation formula is used for data processing. The relationship between the SNR of different result indexes and the signal factors is shown in Figure 10. Moreover, variance analysis is performed to examine the influence degree of each factor on different outcome indicators, as shown in Figure 11.
The influence degree of every signal factor on each outcome indicator is different. As can be seen from Figure 11, for the indicator (a), the contribution rate of each signal factor is ranked as C > B E > D A , where the effect of factor B and factor C is huge. For indicator (b), the contribution rate of each signal factor is ranked as A > B > C > E > D , where the effect of factor A is powerful. For the indicator (c), the contribution rate of each signal factor is ranked as A ( B C D E ) , where the effect of factor A is significant. For the indicator (d), the contribution rate of each signal factor is ranked as D > A > B > C > E , where the influence degree of each factor is little different.
In the Taguchi method, the factor that has a significant impact is prioritized, and the factor with insignificant effect is comprehensively considered according to actual needs. According to Figure 10a, the levels of factor B and factor C are selected as B1 and C5, respectively. According to Figure 10b,c, the level of factor A is selected as A5. The levels of factor D and factor E are selected as D1 and E4, taking Figure 10a–d into account. Therefore, when considering only the above four performance indicators, the optimal parameter according to the Taguchi method is A5B1C5D1E4. However, the performance indicators that need to be considered in practical application also include miniaturization, impact resistance, processing difficulty, etc. Therefore, it is necessary to further optimize the results obtained by the Taguchi method in accordance with the gyroscope application requirements.
It was found that the modal frequency will be too concentrated when the height of tines h2 is too high. In order to isolate the influence of the interference mode, the level of factor A was optimized from A5 to A4. Combined with the specific structure of the gyroscope, the difference between factor B and factor C cannot be too large, otherwise the two sides of the drive tine will not be sufficiently symmetrical, which will affect the gyroscope mode shape. As known from Figure 10a, within the appropriate range, a smaller value of factor B and bigger value of factor C will result in better vibration performance. Considering miniaturization and the index, the difference between factor B and factor C cannot be too large, therefore, the level of factor B was selected as B1, and the value of factor C was selected as 0.72 mm after further calculation, which was near C2 and was denoted as C2′. Considering the processing difficulty, the level of the factor D was optimized from D1 to D3. Considering the miniaturization, the level of the factor E was optimized from E4 to E2. Therefore, combined with practical application requirements, the further optimization result based on the Taguchi method is A4B1C2′D3E2.
The feature sizes of the gyroscope optimized by the Taguchi method are shown in Figure 12, where the gyroscopic thickness is 0.5 mm.
The feature sizes of the gyroscope before and after Taguchi optimization are shown in Table 5, where the relevant size codes have been shown in Figure 9. The related vibration performance is shown in Table 6.
As can be seen from Table 6, compared with the initial structure, the transmission efficiency ŋ of the optimized gyroscope has improved by about 18%, and the working modal frequency has been reduced by about 2.7 kHz. According to the previous analysis, the improvement of these two indicators will further improve the mechanical sensitivity of the gyroscope.

4. Discussion

According to the simulation results, it was found that slotting in the middle of the base will improve the transmission efficiency of Coriolis vibration, while slotting on the sides of the base will reduce it. In addition, the transmission efficiency was improved by decreasing drive–sense interval d1 and increasing drive–drive interval d2. We believe that this may be due to the following reasons:
(1)
Slotting on both sides of the base will increase the local torsional stiffness of the edge area of the base, which is not conducive to the transmission of Coriolis vibration. Slotting in the middle of the base will increase the local torsional stiffness of the middle area of the base, thereby hindering the loss of Coriolis vibration towards the middle direction, and playing a guiding and promoting role in the transmission of Coriolis vibration;
(2)
The decrease in drive–sense interval d1 reduces the distance between the drive tine and the sense tine, thereby reducing energy loss in the transmission process and improving the transmission efficiency;
(3)
The increase in drive–drive interval d2 increases the width of the base, thereby reducing the overall torsional stiffness of the base and facilitating the transmission of Coriolis vibration.
In terms of manufacturing, some adjustments were made to the structure and size optimization process of the gyroscope according to manufacturing process issues. For example, the hole pattern of the gyroscope base was as simple and symmetrical as possible, and the process difficulty was comprehensively considered when optimizing the size. In future studies, these considerations will benefit the fabrication and testing of the proposed quartz tuning fork gyroscope.

5. Conclusions

This paper designed a novel MEMS tuning fork gyroscope microstructure and optimized it from the perspective of mode shape, providing a new idea for the design of MEMS gyroscopes. Firstly, the paper demonstrated the correctness for using the mode shape to reflect the vibration state of the gyroscope. Then the finite element method was utilized to compare the vibration performance of the proposed structure with other typical structures, and the advantages of the designed structure were noted. Finally, the Taguchi method was used to optimize the relevant feature sizes of the gyroscope. In short, the research conclusions can be summarized as follows:
(1)
This gyroscope microstructure worked by transmitting vibrations between the drive tine and sense tine, relatively eliminating the kinematic coupling. It was credible to evaluate the transmission efficiency of Coriolis vibration by using the sensing mode, which can greatly reduce the amount of calculation required;
(2)
During gyroscope structural analysis, slotting in the middle of the base improved the transmission efficiency, and opening arc slots between the tines reduced the working modal order and frequency. Moreover, the base protruding slightly outward improved local symmetry to obtain better vibration performance;
(3)
The height of tines had a large influence on the frequency difference of working mode and the frequency of driving mode, while the tine interval had a large effect on the transmission efficiency of Coriolis vibration. Optimized by the Taguchi method, the transmission efficiency was improved by about 18%, and the working modal frequency was reduced by about 2.7 kHz. The improvement of these two indicators will further improve the mechanical sensitivity of the gyroscope microstructure.

Author Contributions

Writing—original draft preparation, C.X.; conceptualization, P.Z. and W.L.; methodology, C.X. and P.Z.; investigation, F.L. and M.Z.; writing—review and editing Y.H.; project administration, funding acquisition, F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China (Number: 52075208 and U20A6004) and the Aeronautical Science Foundation of China (Number 20160279003).

Acknowledgments

The authors would like to thank the China Aero Polytechnology Establishment for the help to the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Belyanin, L.N.; Vu, D.K. Nutation damping of dynamically tuned gyroscopes. In Proceedings of the MATEC Web of Conferences, Tomsk, Russia, 18–20 April 2017; p. 1008. [Google Scholar]
  2. Bencze, W.J.; Brumley, R.W.; Eglington, M.L.; Buchman, S. Precision electrostatic suspension system for the gravity probe B relativity mission’s science gyroscopes. In Proceedings of the SICE 2003 Annual Conference, Fukui, Japan, 4–6 August 2003; pp. 2726–2731. [Google Scholar]
  3. Di Virgilio, A.; Schreiber, K.U.; Gebauer, A.; Wells, J.P.; Tartaglia, A.; Belfi, J.; Beverini, N.; Ortolan, A. A laser gyroscope system to detect the Gravito-Magnetic effect on Earth. Int. J. Mod. Phys. D 2010, 19, 2331–2343. [Google Scholar] [CrossRef] [Green Version]
  4. Usui, R.; Ohno, A. Recent progress of fiber optic gyroscopes and applications at JAE. In Proceedings of the 2002 15th Optical Fiber Sensors Conference Technical Digest, Portland, OR, USA, 10 May 2002; pp. 11–14. [Google Scholar]
  5. Xu, Q.; Xiao, D.; Hou, Z.; Zhuo, M.; Li, W.; Xu, X.; Wu, X. A novel high-sensitivity butterfly gyroscope driven by horizontal driving force. IEEE Sens. J. 2019, 19, 2064–2071. [Google Scholar]
  6. Xia, D.; Yu, C.; Kong, L. The development of micromachined gyroscope structure and circuitry technology. Sensors 2014, 14, 1394–1473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Guo, X.; Sun, C.; Wang, P.; Huang, L. Vision sensor and dual MEMS gyroscope integrated system for attitude determination on moving base. Rev. Sci. Instrum. 2018, 89, 015002. [Google Scholar] [CrossRef]
  8. Shkel, A.M. Micromachined gyroscopes: Challenges, design solutions, and opportunities. Smart Structures and Materials 2001: Smart Electronics and MEMS. Int. Soc. Opt. Photonics 2001, 4334, 74–85. [Google Scholar]
  9. Yazdi, N.; Ayazi, F.; Najafi, K. Micromachined inertial sensors. Proc. IEEE 1998, 86, 1640–1659. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, K.; Zhang, W.; Chen, W.; Li, K.; Dai, F.; Cui, F.; Wu, X.; Ma, G.; Xiao, Q. The development of micro-gyroscope technology. J. Micromech. Microeng. 2009, 19, 113001. [Google Scholar] [CrossRef]
  11. Nguyen, M.N.; Ha, N.S.; Nguyen, L.Q.; Chu, H.M.; Vu, H.N. Z-Axis micromachined tuning fork gyroscope with low air damping. Micromachines 2017, 8, 42. [Google Scholar] [CrossRef]
  12. Weinberg, M.S.; Kourepenis, A. Error sources in in-plane silicon tuning-fork MEMS gyroscopes. J. Microelectromechanical Syst. 2006, 15, 479–491. [Google Scholar] [CrossRef]
  13. Wang, R.; Cheng, P.; Xie, F.; Young, D.; Hao, Z. A multiple-beam tuning-fork gyroscope with high quality factors. Sens. Actuators A Phys. 2011, 166, 22–33. [Google Scholar] [CrossRef]
  14. Li, Z.; Gao, S.; Jin, L.; Liu, H.; Guan, Y.; Peng, S. Design and mechanical sensitivity analysis of a MEMS tuning fork gyroscope with an anchored leverage mechanism. Sensors 2019, 19, 3455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Zhao, Y.; Zhao, J.; Wang, X.; Xia, G.M.; Shi, Q.; Qiu, A.P.; Xu, Y.P. A sub-0.1°/h bias-instability split-mode MEMS gyroscope with CMOS readout circuit. IEEE J. Solid-State Circuits 2018, 53, 2636–2650. [Google Scholar] [CrossRef]
  16. Yang, C.; Li, H. Digital control system for the MEMS tuning fork gyroscope based on synchronous integral demodulato. IEEE Sens. J. 2015, 15, 5755–5764. [Google Scholar] [CrossRef]
  17. Zhang, Z.; Feng, L.; Sun, Y. Study on digital correlation demodulation technology of micro quartz tuning fork gyroscope. In Proceedings of the Conference on Advanced Sensor Systems and Applications IV, Beijing, China, 9 November 2010. [Google Scholar]
  18. Prikhodko, I.P.; Trusov, A.A.; Shkel, A.M. Compensation of drifts in high-Q MEMS gyroscopes using temperature self-sensing. Sens. Actuators A Phys. 2013, 201, 517–524. [Google Scholar] [CrossRef] [Green Version]
  19. Behbahani, A.H.; M’Closkey, R.T. Multimodal tuning of an axisymmetric resonator. J. Dyn. Syst. Meas. Control. 2019, 141, 091010. [Google Scholar] [CrossRef] [Green Version]
  20. Howard, H.G.; Behbahani, A.H.; Gibson, J.S.; M’Closkey, R.T. Adaptive cancellation of parasitic coupling. J. Microelectromechanical Syst. 2018, 27, 844–853. [Google Scholar]
  21. Li, C.; Yang, B.; Guo, X.; Wu, L. A digital calibration technique of MEMS gyroscope for closed-loop mode-matching control. Micromachines 2019, 10, 496. [Google Scholar] [CrossRef] [Green Version]
  22. Ge, H.H.; Behbahani, A.H.; M’Closkey, R.T. MEMS gyro drift compensation using multiple rate measurements derived from a single resonator. In Proceedings of the 2018 IEEE/ION Position, Location and Navigation Symposium (PLANS), Monterey, CA, USA, 23–26 April 2018; pp. 288–293. [Google Scholar]
  23. Trusov, A.A.; Prikhodko, I.P.; Zotov, S.A.; Shkel, A.M. Low-Dissipation silicon tuning fork gyroscopes for rate and whole angle measurements. IEEE Sens. J. 2011, 11, 2763–2770. [Google Scholar] [CrossRef]
  24. Bernstein, J.; Bancu, M.; Cook, E.; Henry, T.; Kwok, P.; Nyinjee, T.; Perlin, G.; Teynor, B.; Weinberg, M. Diamond hemispherical resonator fabrication by isotropic glass etch. In Proceedings of the Solid-State Sensors, Actuators, Microsystems Workshop, Hilton Head Island, SC, USA, 8–12 June 2014; pp. 273–276. [Google Scholar]
  25. Xu, Q.; Hou, Z.; Kuang, Y.; Miao, T.; Ou, F.; Zhuo, M.; Xiao, D.; Wu, X. A tuning fork gyroscope with a polygon-shaped vibration beam. Micromachines 2019, 10, 813. [Google Scholar] [CrossRef] [Green Version]
  26. Wu, X.; Xie, L.; Xing, J.; Dong, P.; Wang, H.; Su, J. A z-axis quartz tuning fork micromachined gyroscope based on shear stress detection. IEEE Sens. J. 2011, 12, 1246–1252. [Google Scholar] [CrossRef]
  27. Xu, P.; Wei, Z.; Jia, L.; Zhao, Y.; Han, G.; Si, C.; Ning, J.; Yang, F. ZRO drift reduction of MEMS gyroscopes via internal and packaging stress release. Micromachines 2021, 12, 1329. [Google Scholar] [CrossRef] [PubMed]
  28. Bai, B.; Li, C.; Zhao, Y. Development of an effective method to reduce mechanical coupling error in a micro quartz tuning fork gyroscope. Sens. Actuators A Phys. 2021, 332, 113058. [Google Scholar] [CrossRef]
  29. Bishop, R.E.D.; Johnson, D.C. The Mechanics of Vibration; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  30. Canel, T.; Zeren, M.; Sınmazçelik, T. Laser parameters optimization of surface treating of Al 6082-T6 with Taguchi method. Opt. Laser Technol. 2019, 120, 105714. [Google Scholar] [CrossRef]
  31. Bao, M.H. Micro mechanical transducers-pressure sensors, accelerometers and gyroscopes. In Handbook of Sensors and Actuators; Elsevier: Amsterdam, The Netherlands, 2000; pp. 139–198. [Google Scholar]
Figure 1. Structure of decoupling MEMS tuning fork gyroscope.
Figure 1. Structure of decoupling MEMS tuning fork gyroscope.
Micromachines 13 00172 g001
Figure 2. Working mode of the MEMS gyroscope: (a) driving mode (b) sensing mode.
Figure 2. Working mode of the MEMS gyroscope: (a) driving mode (b) sensing mode.
Micromachines 13 00172 g002
Figure 3. Schematic of the gyroscope electrodes arrangement: (a) overall electrode arrangement, (b) electrode directions of drive tine 1, (c) electrode directions of drive tine 2, (d) electrode directions of sense tine 1, (e) electrode directions of sense tine 2.
Figure 3. Schematic of the gyroscope electrodes arrangement: (a) overall electrode arrangement, (b) electrode directions of drive tine 1, (c) electrode directions of drive tine 2, (d) electrode directions of sense tine 1, (e) electrode directions of sense tine 2.
Micromachines 13 00172 g003
Figure 4. Model of the mode test: (a) type A, (b) type B, and (c) type C.
Figure 4. Model of the mode test: (a) type A, (b) type B, and (c) type C.
Micromachines 13 00172 g004
Figure 5. Mesh model of type A.
Figure 5. Mesh model of type A.
Micromachines 13 00172 g005
Figure 6. Load of the finite element model: (a) harmonic analysis, and (b) modal analysis.
Figure 6. Load of the finite element model: (a) harmonic analysis, and (b) modal analysis.
Micromachines 13 00172 g006
Figure 7. Six other typical models of the gyroscope: (a) type D: initial structure, (b) type E: single hole in the base, (c) type F: double hole in the base, (d) type G: spindle beam between two holes, (e) type H: single hole in horizontal direction, and (f) type I: opening arc slots between the tines.
Figure 7. Six other typical models of the gyroscope: (a) type D: initial structure, (b) type E: single hole in the base, (c) type F: double hole in the base, (d) type G: spindle beam between two holes, (e) type H: single hole in horizontal direction, and (f) type I: opening arc slots between the tines.
Micromachines 13 00172 g007
Figure 8. Bar chart of the transmission efficiency ŋ of various typical structures.
Figure 8. Bar chart of the transmission efficiency ŋ of various typical structures.
Micromachines 13 00172 g008
Figure 9. Feature size codes of the gyroscope.
Figure 9. Feature size codes of the gyroscope.
Micromachines 13 00172 g009
Figure 10. SNR effect diagram of different result index: (a) SNR effect diagram of the transmission efficiency, (b) SNR effect diagram of the frequency difference, (c) SNR effect diagram of the driving modal frequency, and (d) SNR effect diagram of the drive coupling coefficient.
Figure 10. SNR effect diagram of different result index: (a) SNR effect diagram of the transmission efficiency, (b) SNR effect diagram of the frequency difference, (c) SNR effect diagram of the driving modal frequency, and (d) SNR effect diagram of the drive coupling coefficient.
Micromachines 13 00172 g010aMicromachines 13 00172 g010b
Figure 11. The contribution rate of each signal factor to the SNR of different outcome indicators: (a) the contribution ratio of each signal factor to the SNR of the ŋ, (b) the contribution ratio of each signal factor to the SNR of the Δf, (c) the contribution ratio of each signal factor to the SNR of the driving modal frequency, and (d) the contribution ratio of each signal factor to the SNR of the λ.
Figure 11. The contribution rate of each signal factor to the SNR of different outcome indicators: (a) the contribution ratio of each signal factor to the SNR of the ŋ, (b) the contribution ratio of each signal factor to the SNR of the Δf, (c) the contribution ratio of each signal factor to the SNR of the driving modal frequency, and (d) the contribution ratio of each signal factor to the SNR of the λ.
Micromachines 13 00172 g011
Figure 12. The feature sizes of the optimized gyroscope.
Figure 12. The feature sizes of the optimized gyroscope.
Micromachines 13 00172 g012
Table 1. The finite element results of the mode shape.
Table 1. The finite element results of the mode shape.
Type Harmonic AnalysisModal Analysis
Amplitude A (μm)Amplitude B (μm)Phase A (°) Phase B (°) AB/AA Δ φ   ( ° ) AB/AA Δ φ   ( ° )
A24.141111.8684−0.622179.4070.4916180.0290.4806180
B23.04549.2162−0.750179.2730.3999180.0230.3983180
C23.61778.4889−0.765179.2600.3594180.0250.3561180
Table 2. Vibration performance of different typical structures (# represents driving mode frequency, * represents sensing mode frequency).
Table 2. Vibration performance of different typical structures (# represents driving mode frequency, * represents sensing mode frequency).
TypeModal Frequency (Hz)ŋ
Fifth OrderSixth OrderSeventh OrderEighth Order
A17,302.9 (#)17,369.6 (*)17,713.9--0.4806
B--18,460.318,658.5 (#)19,475.1 (*)0.3983
C--18,676.919,198.2 (#)19,639.6 (*)0.3561
D--19,343.119,384.1 (#)19,777.7 (*)0.3857
E--19,106.6 (#)19,264.819,669.4 (*)0.4285
F--18,565.9 (#)19,171.619,436.1 (*)0.4391
G--18,598.8 (#)19,172.719,436.7 (*)0.4385
H--19,186.2 (#)19,305.619,762.5 (*)0.3941
I--17,394.5 (#)17,460.5 (*)17,961.30.4548
Table 3. Levels of the signal factors.
Table 3. Levels of the signal factors.
LevelsABCDE
h2 (mm)d1 (mm)d2 (mm)w1 (mm)h3 (mm)
13.50.60.660.121
240.70.760.161.5
34.50.80.860.22
450.90.960.242.5
55.511.060.283
Table 4. Levels of the noise factor.
Table 4. Levels of the noise factor.
LevelsF
Mesh Size (mm)
10.06
20.08
Table 5. The feature sizes of the gyroscope before and after Taguchi optimization.
Table 5. The feature sizes of the gyroscope before and after Taguchi optimization.
Feature Sizesd0d1d2h1h2th3h4w1w2
Initial Value (mm)0.460.80.863.54.50.520.50.21.19
Optimized Value (mm)0.460.60.722.850.51.50.50.21.02
Table 6. The vibration performance of the gyroscope before and after Taguchi optimization.
Table 6. The vibration performance of the gyroscope before and after Taguchi optimization.
Drive Modal OrderSense Modal OrderDrive Modal Frequency (Hz)Sense Modal Frequency (Hz)Frequency Difference ΔF (Hz)The Transmission Efficiency ŋ
Initial6516,832.916,623.2209.70.4632
Optimized5613,907.814,101.5193.70.5466
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xiong, C.; Zeng, P.; Lv, W.; Lu, F.; Zhang, M.; Huang, Y.; Zhu, F. Design and Optimization of a Novel MEMS Tuning Fork Gyroscope Microstructure. Micromachines 2022, 13, 172. https://doi.org/10.3390/mi13020172

AMA Style

Xiong C, Zeng P, Lv W, Lu F, Zhang M, Huang Y, Zhu F. Design and Optimization of a Novel MEMS Tuning Fork Gyroscope Microstructure. Micromachines. 2022; 13(2):172. https://doi.org/10.3390/mi13020172

Chicago/Turabian Style

Xiong, Chuanguo, Pengjun Zeng, Weishan Lv, Fengming Lu, Ming Zhang, Yuhua Huang, and Fulong Zhu. 2022. "Design and Optimization of a Novel MEMS Tuning Fork Gyroscope Microstructure" Micromachines 13, no. 2: 172. https://doi.org/10.3390/mi13020172

APA Style

Xiong, C., Zeng, P., Lv, W., Lu, F., Zhang, M., Huang, Y., & Zhu, F. (2022). Design and Optimization of a Novel MEMS Tuning Fork Gyroscope Microstructure. Micromachines, 13(2), 172. https://doi.org/10.3390/mi13020172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop