Next Article in Journal
A New Methodology for 3D Target Detection in Automotive Radar Applications
Next Article in Special Issue
An Improved Alignment Method for the Strapdown Inertial Navigation System (SINS)
Previous Article in Journal
Photon Counting Imaging with an Electron-Bombarded Pixel Image Sensor
Previous Article in Special Issue
Sample-Wise Aiding in GPS/INS Ultra-Tight Integration for High-Dynamic, High-Precision Tracking
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Design of a Center Support Quadruple Mass Gyroscope (CSQMG) †

Engineering Research Center for Navigation Technology, Department of Precision Instruments, Tsinghua University, Beijing 100084, China
*
Authors to whom correspondence should be addressed.
This paper is an extended version of our paper entitled “Innovation of Flat Gyro: Center Support Quadruple Mass Gyroscope”. In Proceedings of the 2016 IEEE International Symposium on Inertial Sensors and Systems, Laguna Beach, CA, USA, 22–25 February 2016.
Sensors 2016, 16(5), 613; https://doi.org/10.3390/s16050613
Submission received: 24 March 2016 / Revised: 18 April 2016 / Accepted: 20 April 2016 / Published: 28 April 2016
(This article belongs to the Special Issue Inertial Sensors and Systems 2016)

Abstract

:
This paper reports a more complete description of the design process of the Center Support Quadruple Mass Gyroscope (CSQMG), a gyro expected to provide breakthrough performance for flat structures. The operation of the CSQMG is based on four lumped masses in a circumferential symmetric distribution, oscillating in anti-phase motion, and providing differential signal extraction. With its 4-fold symmetrical axes pattern, the CSQMG achieves a similar operation mode to Hemispherical Resonant Gyroscopes (HRGs). Compared to the conventional flat design, four Y-shaped coupling beams are used in this new pattern in order to adjust mode distribution and enhance the synchronization mechanism of operation modes. For the purpose of obtaining the optimal design of the CSQMG, a kind of applicative optimization flow is developed with a comprehensive derivation of the operation mode coordination, the pseudo mode inhibition, and the lumped mass twisting motion elimination. The experimental characterization of the CSQMG was performed at room temperature, and the center operation frequency is 6.8 kHz after tuning. Experiments show an Allan variance stability 0.12°/h (@100 s) and a white noise level about 0.72°/h/√Hz, which means that the CSQMG possesses great potential to achieve navigation grade performance.

1. Introduction

In view of the fact microelectromechanical systems (MEMS) technology is rapidly developing, it has been possible to create gyroscopes with small size, low weight, low cost, and low power under the premise of high precision. Therefore, MEMS gyros are used not only in high precision applications such as aerospace and military, but also a wide range of consumer market uses including cellphones, automotive, and intelligent wear [1,2].
MEMS vibratory gyroscopes generally utilize Coriolis’ force to sense angular velocity. To operate the CSQMG under drive mode, when an angular velocity is introduced, the moving part of the gyroscope experiences a Coriolis force, and the gyro induces a sense mode oscillation. The input angular velocity can then be measured by detecting the sense mode vibration [3]. Most MEMS gyros are of the vibratory type, and high sensitivity dual-mass Tuning Fork Gyros (TFGs) play an important role. In the tuning fork design, linear acceleration/shock signals are rejected by its symmetric tuning fork architecture [4,5], however, the performance and environmental robustness are limited by the asymmetric sense mode [6,7].
Since the fully symmetric structures with higher performance were identified in recent years, research on micro 3-D symmetric structures has been emphasized. Micro-hemispheric structures have been demonstrated to achieve high Q and high precision potential, while many different manufacturing processes such as glass blowing, silicon etching, and integrated precision machining have been studied to create high precision 3-D hemispherical resonators [8,9,10,11]. As the exploration of new fabrication technology has made the micro HRG become a new research hotspot, the development potential of symmetrical structures has come to be highly regarded at the same time [12].
A new coupling design of MEMS tuning-fork structures, Quad Mass Gyro (QMG) sealed with getter, was proposed to approach gyro compassing and navigation grade performance [13]. The design of the coupling beam, which creates a high frequency separation between anti-phase and in-phase vibratory modes, and balanced operation mode eliminates energy dissipation through the substrate [14,15,16,17]. However, this design causes more folded beams, which introduce more error factors and consequently increase the design and fabrication difficulty [18].
Considering the increasing maturity of the highly batch consistent Silicon-on-Glass/Silicon-on-Insulator (SOG/SOI) planar process, rational MEMS gyro structure pattern with symmetric form is a feasible development direction to achieve breakthrough performance [19]. Recently, a new structure adopting a simpler flexure pattern and a single center support anchor to support four proof masses was introduced [20], which integrates the advantages of TFG (large mass and high resolution) and HRG (fully symmetric and high Q) into a single flat structure (shown in Figure 1). Compared to conventional flat structures, the new design forms a fully symmetric configuration, and adjusts the mode frequency distribution through serpentine beams. Therefore, the fabrication process robustness is enhanced.
Besides, the Y-shaped coupling beam component is introduced to solve the problem of in-phase and anti-phase mode frequency confounding, and this design enhances the synchronization mechanism of operation modes. Because of its simple structure, without introducing complicated structures, the CSQMG achieves similar operation modes to HRGs (0° and 45° vibration modes) by its double tuning fork configuration. Thus the CSQMG can achieve the same environmental robustness as HRGs.
Along with the benefits of the new design, the parameters of a single beam affect the stiffness acting on more than one direction, which means more factors should be considered in the traditional structure optimization processes. At the same time, the lumped masses with two support points can be affected by torsion coupling stiffness. Motivated by the detailed analysis of the structure, a comprehensive derivation of the new design and the function of the serpentine beams is required to show the design rationality. An optimization flow considering multi-factor needs to be put forward in order to simplify the design.
In this paper, a more complete description of the CSQMG is presented. In Section 2, a description of the quadruple mass oscillator architecture with its operation mode is discussed. Section 3 presents a comprehensive derivation of the structure and the function of the serpentine beams, and this model is evaluated by Finite-Element Analysis (FEA). An optimization flow based on the model is concluded, and the comprehensive design parameters of the gyro are given in Section 4. Afterwards, the fabrication process and experimental characterization are presented in Section 5. Section 6 concludes this paper with a discussion and summary of the results.

2. Operation Mode

The inspiration of CSQMG comes from HRGs, especially the n = 2 wine glass mode, as depicted in Figure 2. In order to transform the HRGs into a planar structure, the 0°/90° resonance mode and the −45°/45° resonance mode of the n = 2 wineglass mode are regarded as drive mode and sense mode, respectively. Then the hemisphere structure is simplified to planar structure by the following two steps. Firstly, as for drive mode, the distributed mass is concentrated at the harmonic displacement peak of the hemisphere resonator, forming four lumped masses. The four placement points of the resonator are supported through the bracing structure forming four fixed points (illustrated in Figure 2a,e). Secondly, as for sense mode, considering the tangential motion component of hemisphere resonator mass elements, the resonant elements are turned into four lumped masses alternating towards and away from the central point (illustrated in Figure 2b,f).
Though the latter transformation can’t match the sense mode resonance perfectly, the hemispherical resonator is thus creatively transformed into an in-plane dual tuning fork resonator. It is obvious that, when the quadruple mass resonator works as a gyro, the four lumped masses vibrate along the radius of the symmetry center in the drive mode and vibrate along the tangential orientation under the sense mode. Moreover, in the polar coordinates reference frame, the two adjacent lumped masses oscillate in anti-phase motion, which is the same as tuning fork anti-phase mode. Compared to conventional dual-mass TFGs, the major advantage of this symmetric CSQMG is the 4-fold symmetrical axes pattern. With in-plane dual tuning fork resonator and differential signal extraction, the common mode linear and centrifugal interferences such as acceleration and shock are rejected.
In drive mode, the lumped mass is restrained by the bilateral symmetrical structure and only radial motion is permitted. However, a coupling component along the lumped mass twisting motion potentially exists during the sense mode resonance, and it should be eliminated during the design optimization process.
Accompanied by the transformation, in-phase modes are also inevitable in CSQMG. Instead of vertical mode and tilting mode acting as major pseudo modes of HRGs, “constriction mode” and “twisting mode” become the major pseudo modes of the new structure. In constriction mode, the resonance alternates between contraction and expansion (depicted in Figure 2c,g), while in twisting mode, the resonance alternates between clockwise and counterclockwise (depicted in Figure 2d,h). Although the gyro can also operate in in-phase modes, the common mode signal cannot be eliminated.
The application of the N-shaped beam and the Y-shaped coupling beam enables a vibration of two in-plane quadrature directions with a single resonant beam [15]. Collectively, the four lumped masses, four resonant beams and four fixed points form a circular resonator, and the optimized distribution of each mode’s natural frequency is obtained by changing the design parameters.
The bracing structure is transformed into a single anchor support form in order to achieve low anchor loss [21], like HRGs [9]. As the four placement points of the circle resonator simplify the support form, a symmetric generalized force is exerted on the four fixed points during the sense mode resonance. In the ideal situation, the four points are regarded as zero-force points in the sense mode (Figure 3a). While in drive mode, because of the anti-symmetric generalized force, such an equivalence is no longer applicable, so a mat-shaped bracing structure is designed to eliminate the placement point torque at the four center points of each curved beam. These points become zero-force points in the drive mode (Figure 3b). Therefore, the single center supporting anchor is regarded as a zero-force point in both operation modes.

3. Theoretical Model

As mentioned in Section 2, on account of the N-shaped beam and the Y-shaped coupling beam design, the four lumped masses oscillate synchronously as a single circle resonator in comparison to conventional TFGs. Unlike the usual folded beam that determines stiffness of only one orientation, this design enables a vibration of two in-plane quadrature directions with a single resonant beam. This makes it difficult to establish a second-order system equation for the single resonant mass. As there are no referential design guidelines, the relationship between the geometry and the frequency distribution of the modes needs to be confirmed.

3.1. Derivation of the Flexibility Matrix

Taking account of the 4-fold symmetrical axes pattern, we simplify the analysis as a symmetric model with the following assumptions [18]:
  • The lumped mass is regarded as a rigid body of length 2q, and all of the generalized forces are concentrated in the barycenter of the mass.
  • The mass of the beams is neglected and the placement points are treated as fixed during oscillation.
  • Following the hypotheses of linear strain and small deformation, and the effect of tension and compression stresses is neglected in the analysis process.
  • As the Y-shaped coupling beam and barycenter of the lumped mass are symmetric axes of the structure, we presume an identical symmetric or anti-symmetric orthogonal generalized force is applied on both sides.
The simplified model is still a statically indeterminate system. We solve the system with the statics flexibility method [22,23], remove the constraint at the lumped masses symmetric axes, exert external generalized unit forces FM1, FM2, MM and assume internal generalized unit forces FB1, FB2, MB, as shown in Figure 4. Then the canonical equation of deformation (1) based on the virtual displacement principle is established, and the torque expressions xak~xfk (k ∈ [1, 9] ∩ kN) are listed in Table 1 based on Figure 5.
k = 1 9 1 E k I k 0 l k ( x d k 2 x d k x e k x d k x f k x d k x e k x e k 2 x e k x f k x d k x f k x e k x f k x e k 2 ) d x k ( X p 1 X p 2 X p 3 ) + k = 1 9 1 E k I k 0 l k ( x p k x d k x p k x e k x p k x f k ) d x k = 0
In the canonical Equation (1), Xp1, Xp2, Xp3 (p = a, b, c) represent the internal forces under the action of external force FM1 (p = a), FM2 (p = b), MM (p = c), respectively.
Internal generalized unit forces are obtained by solving the canonical equation. Based on Castigliano’s second theorem, the flexibility matrix derivation of operation modes (anti-phase modes) is formulated, as illustrated in Equation (2):
( x y θ ) = ( S 11 0 0 0 S 22 S 23 0 S 32 S 33 ) ( F x F y M ) S 11 = k = 1 9 1 E k I k 0 l k ( x a k + X a 1 x d k + X a 3 x f k ) 2 d x k S 22 = k = 1 9 1 E k I k 0 l k ( x b k + X b 2 x e k ) 2 d x k S 33 = k = 1 9 1 E k I k 0 l k ( x c k + X c 2 x e k ) 2 d x k S 23 = S 32 = k = 1 9 1 E k I k 0 l k ( x c k + X c 2 x e k ) ( x b k + X b 2 x e k ) d x k
S11 represents the flexibility of anti-phase drive orientation, S22 represents the flexibility of anti-phase sense orientation, and S33 represents the flexibility of lumped mass torsion orientation. It’s noteworthy that S23 is non-vanishing which means a coupling component along the lumped mass twisting motion exists in the sense mode. This is an undesirable characteristic in this structure as well, which should be considered in the subsequent optimization flow.
For in-phase modes, the symmetry of internal forces and external forces interchange. By changing the external force orientations with the internal force orientations, listing torque expressions, and using the same calculation steps, the flexibility matrix derivation of pseudo modes is also built.

3.2. FEM Validation of Theoretical Model

For the purpose of verifying the accuracy of the theoretical model, the geometric parameters obtained by optimization method are ported into FEM software ANSYS, as depicted in Figure 6. To reduce the amount of computation, the model is meshed with two kinds of grid division. Swept meshing with element size less than 0.02 mm is used on the beams, and default values with element size less than 0.2 mm is adopted for the rest. To reduce unnecessary division, the comb fingers are replaced by an evenly distributed mass as well. Static load is applied in accordance with the resonance form of the CSQMG, by means of fixing the single anchor and exerting a unit force on the lumped mass.
The average stiffness error of the four modes is less than 4%, proving the accuracy of the modeling process (Table 2). The main errors of the model come from the connection points of each beams which do not fit the assumptions during the simulation process.

4. Design Optimization Process

Optimal natural frequency allocation of the CSQMG is the primal goal of the design optimization. In consideration of the length and width of the serpentine beams containing 18 variables and the varying degree of the parameter impacts, the analysis of the natural frequency configuration is divided into N-shaped beam and Y-shaped beam, respectively. To alleviate the optimization complexity, the control variables method is used to analyze the two components. Then a targeted independent optimization step to eliminate the coupling component between the sense mode motion and the lumped mass twisting motion is discussed. After the overall structure parameters are confirmed, an optimized design of the comb fingers based on reducing the mechanical noise is performed.

4.1. Optimal Design of the N-Shaped Beam

The model derived in Section 3.1 is used to the analysis of the N-shaped beams. Even taking no account of the effect of the Y-shaped coupling beam, the relationship between the static characteristics and the parameters of the beams are still complicated. The parameters contain a total of 16 variables, namely the width and the length of Beam No. 1 to No. 9 except Beam No. 5 (the Y-shaped coupling beam). Therefore, the model needs to be simplified for convenience of analysis. Further presuppositions are as follows:
  • The parameters of serpentine beams on both sides of the symmetric axes are considered equal.
  • The length of the beams is fixed for subsequent analysis later, as the stiffness is proportional to the length and the third power of the width, and the physical dimension is limited.
  • The width of Beam No. 2 and Beam No. 4 are regarded identical, since the effect of Beam No. 2 and Beam No. 4 are almost the same.
As shown in Figure 7, the simplifying assumptions above reduce the number of variables of the N-shaped beams from sixteen to three, i.e., the width of Beam No. 1 (B1), the width of Beam No. 2 and No. 4 (B2 = B4), and the width of Beam No. 3 (B3). In order to show the variation trend between the N-shaped beams and the frequency distribution of the modes, a Cartesian coordinate system is established using the three variables, and the differences of the stiffness of every two different modes are shown using a color code. The Cartesian coordinate system forms three four-dimensional images, as shown in Figure 8.
Figure 8a shows the stiffness difference between the two operation modes. The red area means the stiffness of drive mode is larger than the sense mode, while the blue area indicates the opposite tendency. The closer to the red region, the larger the difference. As the CSQMG operates under the rate mode, the natural frequency of the drive mode should be set equal to or slightly higher than that of the sense mode. As shown in the figure, with the rise of B2, B3, B4 and the drop of B1, the stiffness difference between the two modes grows. The data points with stiffness differences greater than 0 but less than 50 N/m are marked with black spots as the Optimal Values Area 1 in the coordinate system. As the using of numerical methods limits the number of sample points, the range of black spots forms an approximate three-dimensional surface.
Figure 8b demonstrates the stiffness difference between the in-phase drive mode and the anti-phase drive mode. The red area means the stiffness of the in-phase drive mode is larger than the anti-phase drive mode. Shown from the figure, the stiffness difference is significantly affected by B2 and B4, but only slightly affected by B1 and B3. The difference becomes larger with the increase of B2 and B4. In order to get a gyro with better performance and reduce the impact of pseudo mode, the natural frequency of pseudo modes should be kept away from the operation mode. The data points with stiffness difference values greater than 100 N/m are marked with black spots forming Optimal Values Area 2.
Similar to Figure 8b, Figure 8c shows the stiffness difference between the in-phase sense mode and the anti-phase drive mode. The red area means the stiffness of the anti-phase drive mode is larger than the in-phase sense mode. The variation of B2 and B4 have less influence on stiffness changes, and with the increase of B1 and the decrease of B3 the difference become significant. The data points with stiffness difference values greater than 100 N/m are marked with black spots as the Optimal Values Area 3.
In summary, the variation of the parameter B1, B2, B3, and B4 changes the natural frequency distribution of each mode, and the variable ranges of the parameters are limited by the intersection of the three Optimal Values Area, especially the Optimal Values Area 1. In order to expand the scope, and increase the frequency difference between the operation modes and the pseudo modes, the Y-shaped coupling beam should be further analyzed.

4.2. Optimal Design of the Y-Shaped Coupling Beam

For the purpose of distinguishing each vibratory mode resonance frequency, a new type of Y-shaped coupling beam is put forward with the following advantages: (1) low relative machining error based on a simple structure; (2) in-phase modes separated from anti-phase modes with a single structure; (3) less energy is transforred to the substrate by locating it at the placement point.
Connecting the support frame and the circle resonator at the placement points which act as symmetric axes, the Y-shape beam bears generalized forces from both sides. Considering the 4-fold symmetric axes pattern, an anti-symmetric generalized force is applied on both sides under anti-phase drive mode and in-phase sense mode, while a symmetric generalized force is applied in in-phase drive mode and anti-phase sense mode, as shown in Figure 9a. Intuitively, the design of the Y-shaped coupling beam changes the resonant stiffness of four planar modes.
In order to illustrate the influence of the Y-shaped beam more clearly, Figure 9b is presented. It is observed that by changing the Y-shaped beam parameters which include length increase and width decrease, the frequencies of the anti-phase drive mode and in-phase sense mode are shifted below the anti-phase sense mode. Similarly, with the parametric variation of beams, the stiffness of operation modes matches together, forming a range of optimal values. In addition, a small length value and width value are needed for low elastic damping loss purpose, while shifting above out-of-plane pseudo mode stiffness requires a large value of the parameters. Taking all these conditions into consideration, the optimal value range can be tightened to a narrow range.

4.3. Elimination of the Twisting Motion Coupling

As mentioned in Section 2, the S23 of the flexibility matrix is non-vanishing, which means a coupling component exists between the lumped mass twisting motion and the sense mode motion. This phenomenon indirectly affects the coupling oscillation of drive mode and sense mode during the lumped mass planar motion, which deteriorates the performance of the CSQMG.
An optimization method should be applied to eliminate the coupling. Observing the simplified model, Beam No. 2 and No. 4 are perpendicular to the direction of sense mode vibration, and a forty-five-degree angle is formed between Beam No. 5 and the sense mode vibration. As the width of the beams and the parameters of the Beam No. 5 have been optimized in the above section, an optimization result should be provided via an analysis of the relationship between L2, L4 and the twisting motion coupling stiffness (Figure 10b). With the increasing of L2 and decreasing of L4, the coupling stiffness decreases and intersects the zero plane forming a straight line. This line is the optimal value range of L2 and L4 when the rest of the parameters are determined. Moreover, by means of changing the parameter of L2 and L4 and keeping the summation of these two parameters as a constant, a reasonable allocation is completed without disturbing the frequency distribution sharply.

4.4. FEM Validation of the Optimization Flow

Because chamfers are used at the connection points of each beams, the stress concentrations are reduced and fabrication process consistency is ensured. There is a slight discrepancy between the FEM analysis and the practical design. This discrepancy makes the optimization model unsuitable for the quantitative optimization task. However, the model is still effective in parameterized modeling as a qualitative optimization guide. Optimized analysis of the structure design proceeds through two steps: firstly, iterative calculation and analysis by changing the design of the N-shaped beam and Y-shaped beam are adopted, and these two conditions are taken into account at the same time: (a) the frequency of the drive mode should be slightly higher than the sense mode but no more than 20 Hz (considering about the range of electrostatic tuning) to tolerate the variation in the fabrication process; (b) the frequency of pseudo modes should be separated from the operation mode to more than 150 Hz to reduce the interference of pseudo modes. Secondly, L2 and L4 are changed and the summation of these two parameters is kept constant to find the minimum absolute value of the coupling component between the sense mode motion and the lumped mass twisting motion.
All the parameters of the actual CSQMG structure optimization are shown in Table 3, using the above two steps. By optimizing the width of the beams, the drive mode frequency is 2 Hz higher than the sense mode, the pseudo mode frequency is distinguished from the operation mode to 150 Hz and 434 Hz, respectively. Through optimizing the length of L2 and L4, the rotational coupling of sense mode is reduced from 0.022 N−1 to 4.2 × 10−4 N−1 while the frequency distribution is only slightly changed.

4.5. Optimal Design of Comb Fingers

For the conventional line vibration MEMS gyroscope, the main noise sources are as follows: the mechanical-thermal noise, the preamplifier noise, the digital quantization noise, and the ADC conversion noise. Compared to the mechanical-thermal noise, the digital quantization noise and the ADC conversion noise are a small amount in the CSQMG [24]. The mechanical-thermal noise of the conventional line vibration MEMS gyroscope is:
Ω t h = x n o i s e x sec ( 1 ° / s ) = 90 π A 0 ω n 1 4 ω n 2 k B T m Q 2
where A0 is the amplitude of the drive mode, ωn1 is the natural frequency of the drive mode, ωn2 is the natural frequency of the sense mode, the kB is the Boltzmann constant, T is the Kelvin temperature, m is the quality of lumped masses, and Q2 is the quality factor of the sense mode.
In order to achieve large displacements, which are necessary for low noise operation, 4 × 100 pairs of slide film comb fingers with a 40 μm long stroke are set in the driving direction [25]. Multi-comb fingers result in a maximum driving force of about 16 μN on each lumped mass, when the gyro operates at 10 V. 4 × 52 pairs of squeeze film comb fingers with differential capacitor type leading the sense direction possess over 3.3 pF (each lumped mass) total capacitance and 56 N/m stiffness tuning. The force balance electrodes with 4 × 14 pairs comb fingers guarantee 4 μN balance force at least. Chip size with footprint is 7700 μm × 7700 μm and the design parameters of the CSQMG are listed in Table 4. Substituting the parameters of the CSQMG, such as m, ωn1, ωn2, Q2, A0, into Equation (3), and the mechanical-thermal noise Ωth is calculated to be 0.186°/h/√Hz.

5. Fabrication and Characterization

5.1. Fabrication Process

A MEMS gyro sensitive structure fabrication process based on SOG and inductively coupled plasma deep reactive ion etching (ICP DRIE) with an aspect ratio of about 20 developed by Peking University is shown in Figure 11a [26,27,28]. The starting wafers are 4 inch highly doped silicon wafer produced by OKMETIC (Vantaa, Finland) and a Pyrex 7740 glass wafer. On the Pyrex wafer, a 200 nm thick Ti/Pt/Au layer is patterned by a liftoff process (Figure 11(a1)). On the silicon wafer steps of 20 μm height are etched by DRIE to define the anchor areas and the wafer is doped by phosphorus ion implantation to obtain good Ohmic contacts at the anchor areas (Figure 11(a2)). Afterwards the two wafer are bonded together (Figure 11(a3)). The silicon substrate is then thinned to about 100 μm using KOH (Figure 11(a4)) and released by a second DRIE step (Figure 11(a5)). Finally, the wafers are diced and wire bonded. Figure 11b,c show a SEM image and photograph of the fabricated gyroscope, respectively.

5.2. Frequency Response Curve Analysis

Measurement of the gyro response curve is performed by applying sine excitations on drive combs and reading the signal on the sense combs under a pressure of 30 Pa [3]. Observing from the curve shown in Figure 12a, a total of six resonance peaks appear at 5.7, 6.2, 6.7, 6.8, 7.0 and 7.1 kHz approximately. Considering the influence of processing errors, each resonance peak matches the frequency of in-phase sense mode, pseudo mode 1, anti-phase drive mode (drive mode), anti-phase sense mode (sense mode), pseudo mode 2, and in-phase drive mode, respectively. Comparing to the FEM simulation results, the maximum deviation of the modal frequency is approximately 13.1%, and the deviation of the operation modes are 8.7% and 7.5%. The reductions of the modal frequency are attributed to the overetching during the fabrication process, and a much larger range of compensation of the beams width should be adopted in any subsequent improvement based on the theoretical model analysis. At the same time, due to the frequency response curve analysis, the theoretical model cannot reflect all of the CSQMG’s planar vibration information. Therefore, a mass-spring-damper model with multiple degrees of freedom should be put forward to simplify the CSQMG’s symmetric stiffness analyses, and a system parameter identification basing on experimental data needs to be studied in the future work.

5.3. Experimental Characterization

The CSQMG is packaged in the metal case, with a vacuum level of 30 Pa [20]. The quality factors of the drive mode and sense mode are measured to be about 8500 and 500, respectively. The key parameters of the CSQMG have been tested by a rate table from ACUITAS GmbH (Altendorf, Switzerland) as shown in Figure 12b. The MEMS and ASIC are assembled in an aluminum case, which is shown in Figure 12b. The gyro outputs are digitally transferred by RS422 protocol and collected by PC.
By using the rate table, the scale factor characterization can be evaluated [29]. This CSQMG’s measurement range is ±300°/s, and the scale factor nonlinearity is 96 ppm (Figure 12c). A 2 h static test was performed at room temperature and the results show that the white noise level is about 0.7185°/h/√Hz and the Allan variance stability is 0.12°/h (@100 s), (Figure 12d). Compared to the conventional TFGs, the CSQMG achieves a great performance improvement on account of its innovative structural design and its optimal design flow. According to the design parameters of the structure, shown in Table 5, the designed value of the mechanical noise level is about four times smaller than the measurement result, which means the main noise of the CSQMG is caused by the circuits, and the measured white noise level could be reduced by improving the performance of the electronics.

6. Conclusions

This paper has proposed and demonstrated a new double tuning fork gyro with a single support and 4-fold symmetric axes. The test results, especially the white noise level, shows the CSQMG has great potential to achieve navigation grade performance. As proved by both the mechanical model analysis and FEM analysis, its simpler coupling beam component, a Y-shaped coupling beam, separates the natural frequency of operation mode from pseudo mode by at least 200 Hz. By optimizing the coupling component between the sense mode motion and the lumped mass twisting motion can be eliminated. Meanwhile, the CSQMG natural frequency distribution is effectively controlled by the optimization flow which is based on a comprehensive derivation of the N-shaped beam, Y-shaped beam, and the effectiveness has been proved by the frequency response curve experiments. The derivation will lay the foundations for the error analysis in further study.
Future works will focus on the circuit optimization, quality factor promotion and MEMS processing improvements. CSQMGs are expected to approach the goal performance of HRGs by increasing the structure thickness. A consistent further study of CSQMG based on rate integrated mode will also be realized.

Acknowledgments

The authors acknowledge the technical and scientific staff of the Institute of Microelectronics, Peking University for the support in fabrication.

Author Contributions

Bin Zhou conceived the structure form; Tian Zhang, Bin Zhou and Rong Zhang designed the gyro structure, completed the finite element simulation; Tian Zhang, Rong Zhang and Zhiyong Chen derived the mechanical model; Tian Zhang, Bin Zhou and Peng Ying performed the experiments, analyzed the data; Tian Zhang, Bin Zhou, Rong Zhang, Peng Yin and Zhiyong Chen wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Froyum, K.; Goepfert, S.; Henrickson, J.; Thorland, J. Honeywell micro electro mechanical systems (MEMS) inertial measurement unit (IMU). In Proceedings of the 2012 IEEE/ION Position Location and Navigation Symposium (PLANS), Myrtle Beach, SC, USA, 23–26 April 2012; pp. 831–836.
  2. Xia, D.; Yu, C.; Kong, L. The development of micromachined gyroscope structure and circuitry technology. Sensors 2014, 14, 1394–1473. [Google Scholar] [CrossRef] [PubMed]
  3. Acar, C.; Shkel, A. MEMS Vibratory Gyroscopes: Structural Approaches to Improve Robustness; Springer Science & Business Media: Dordrecht, The Netherlands, 2008. [Google Scholar]
  4. Sharma, A.; Zaman, F.M.; Amini, B.V.; Ayazi, F. A high-Q in-Plane SOI Tuning Fork Gyroscope. In Proceedings of the 2004 IEEE Sensors, Vienna, Austria, 24–27 October 2004; pp. 467–470.
  5. Weinberg, M.S. How to Invent (or Not Invent) the First Silicon MEMS Gyroscope. In Proceedings of the 2015 IEEE International Symposium on Inertial Sensors and Systems, Hapuna Beach, HI, USA, 23–26 March 2015; pp. 1–5.
  6. Yoon, S.W.; Li, S.W.; Perkins, N.C.; Najafi, K. Vibration Sensitivity of MEMS Tuning Fork Gyroscopes. In Proceedings of the 2007 IEEE Sensors, Atlanta, GA, USA, 28–31 October 2007; pp. 115–119.
  7. Yoon, S.W.; Lee, S.; Najafi, K. Vibration-induced errors in MEMS tuning fork gyroscopes. Sens. Actuators A Phys. 2012, 180, 32–44. [Google Scholar] [CrossRef]
  8. Mirbeik, A.; Tavassoli, V.; Ayazi, F.; Tavassolian, N. Three-Dimensional Super-Wideband Micro-Antenna for High-Resolution Millimeter-Wave Medical Imaging. In Proceedings of the 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Chicago, IL, USA, 26–30 August 2014; pp. 4999–5002.
  9. Cho, J.Y.; Woo, J.-K.; Yan, J.; Peterson, R.L.; Najafi, K. Fused-Silica Micro Birdbath Resonator Gyroscope (-BRG). J. Microelectromech. Syst. 2014, 23, 66–77. [Google Scholar] [CrossRef]
  10. Shao, P.; Mayberry, C.L.; Gao, X.; Tavassoli, V.; Ayazi, F. A Polysilicon Microhemispherical Resonating Gyroscope. J. Microelectromech. Syst. 2014, 23, 762–764. [Google Scholar] [CrossRef]
  11. Wang, Y.; Shkel, A.M. Near-Zero Power Wake up System Enabled by Ultra High Q-Factor Array of MEMS Resonators. In Proceedings of the 4th International Symposium on Energy Challenges and Mechanics—Working on Small Scales, Aberdeen, UK, 11–13 August 2015.
  12. Tsanh-Hung, S.; Nitzan, S.H.; Taheri-Tehrani, P.; Kline, M.H.; Boser, B.E.; Horsley, D.A. Silicon MEMS Disk Resonator Gyroscope With an Integrated CMOS Analog Front-End. IEEE Sens. J. 2014, 14, 3426–3432. [Google Scholar]
  13. Prikhodko, I.P.; Zotov, S.A.; Trusov, A.A.; Shkel, A.M. Foucault Pendulum on a Chip: Angle Measuring Silicon MEMS Gyroscope. In Proceedings of the 2011 IEEE 24th International Conference on Micro Electro Mechanical Systems (MEMS), Cancun, Mexico, 23–27 January 2011; pp. 161–164.
  14. Trusov, A.A.; Schofield, A.R.; Shkel, A.M. Gyroscope Architecture with Structurally Forced Anti-Phase Drive-Mode and Linearly Coupled Anti-Phase Sense-Mode. In Proceedings of the TRANSDUCERS 2009—2009 International Solid-State Sensors, Actuators and Microsystems Conference, Denver, CO, USA, 21–25 June 2009; pp. 660–663.
  15. Trusov, A.A.; Schofield, A.R.; Shkel, A.M. Micromachined rate gyroscope architecture with ultra-high quality factor and improved mode ordering. Sens. Actuators A Phys. 2011, 165, 26–34. [Google Scholar] [CrossRef]
  16. Trusov, A.; Prikhodko, I.P.; Zotov, S.; 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]
  17. Simon, B.R.; Trusov, A.; Shkel, A.M. Anti-Phase Mode Isolation in Tuning-Fork MEMS Using a Lever Coupling Design. In Proceedings of the 2012 IEEE Sensors, Taibei, Taiwan, 28–31 October 2012; pp. 1–4.
  18. Iyer, S.V. Modeling and Simulation of Non-Idealities in a Z-Axis CMOS-MEMS Gyroscope. Ph.D. Thesis, Carnegie Mellon University, Pittsburgh, PA, USA, April 2003. [Google Scholar]
  19. Trusov, A.; Atikyan, G.; Rozelle, D.; Meyer, A.; Zotov, S.; Simon, B.R.; Shkel, A.M. Flat is Not Dead: Current and Future Performance of Si-MEMS Quad Mass Gyro (QMG) System. In Proceedings of the IEEE/ION Position, Location and Navigation Symposium—PLANS 2014, Monterey, CA, USA, 5–8 May 2014; pp. 252–258.
  20. Zhou, B.; Zhang, T.; Yin, P.; Chen, Z.; Song, M.; Zhang, R. Innovation of Flat Gyro: Center Support Quadruple Mass Gyroscope. In Proceedings of the IEEE International Symposium on Inertial Sensors and Systems, Laguna Beach, CA, USA, 22–25 February 2016; pp. 1–4.
  21. Senkal, D.; Ahamed, M.J.; Ardakani, M.H.A.; Askari, S.; Shkel, A.M. Demonstration of 1 Million-Factor on Microglassblown Wineglass Resonators with Out-of-Plane Electrostatic Transduction. J. Microelectromech. Syst. 2015, 24, 29–37. [Google Scholar] [CrossRef]
  22. Shi, Q.; Wang, S.; Qiu, A.; Xu, Y.; Ji, X. Design Principle of Suspension of MEMS Gyroscope. In Proceedings of the 1st IEEE International Conference on Nano/Micro Engineered and Molecular Systems, Zhuhai, China, 18–21 January 2006; pp. 242–245.
  23. Simon, B.R.; Khan, S.; Trusov, A.A.; Shkel, A.M. Mode Ordering in Tuning Fork Structures with Negative Structural Coupling for Mitigation of Common-Mode g-Sensitivity. In Proceedings of the 2015 IEEE Sensors, Busan, Korea, 1–4 November 2015; pp. 1–4.
  24. Leland, R.P. Mechanical-thermal noise in MEMS gyroscopes. IEEE Sens. J. 2005, 5, 493–500. [Google Scholar] [CrossRef]
  25. Senkal, D.; Efimovskaya, A.; Shkel, A.M. Dual Foucault Pendulum gyroscope. In Proceedings of the 2015 Transducers—2015 18th International Conference on Solid-State Sensors, Actuators and Microsystems (TRANSDUCERS), Anchorage, AK, USA, 21–25 June 2015; pp. 1219–1222.
  26. Ding, H.; Liu, X.; Lin, L.; Chi, X.; Cui, J.; Kraft, M.; Yang, Z.; Yan, G. A high-resolution silicon-on-glass axis gyroscope operating at atmospheric pressure. IEEE Sens. J. 2010, 10, 1066–1074. [Google Scholar] [CrossRef]
  27. Guo, Z.; Yang, Z.; Zhao, Q.; Lin, L.; Ding, H.; Liu, X.; Cui, J.; Xie, H.; Yan, G. A lateral-axis micromachined tuning fork gyroscope with torsional Z-sensing and electrostatic force-balanced driving. J. Micromech. Microeng. 2009, 20, 025007. [Google Scholar] [CrossRef]
  28. Chu, Y.; Dong, J.; Chi, B.; Liu, Y. A Novel Digital Closed Loop MEMS Accelerometer Utilizing a Charge Pump. Sensors 2016, 16, 389. [Google Scholar] [CrossRef] [PubMed]
  29. Han, F.; Sun, B.; Li, L.; Wu, Q. Performance of a sensitive micromachined accelerometer with an electrostatically suspended proof mass. IEEE Sens. J. 2015, 15, 209–217. [Google Scholar]
Figure 1. Optical photograph of a Center Support Quadruple Mass Gyroscope, with four lumped masses connected through the N-shape-beams and the Y-shape-beams, supported by a single center support anchor.
Figure 1. Optical photograph of a Center Support Quadruple Mass Gyroscope, with four lumped masses connected through the N-shape-beams and the Y-shape-beams, supported by a single center support anchor.
Sensors 16 00613 g001
Figure 2. FEA shows the operation modes and pseudo modes of CSQMG, comparing to the resonant mode of the HRGs. In the operating modes of CSQMG, also known as anti-phase modes, the four lumped masses vibrate along the radial of the symmetry center in the drive mode, and vibrate along the tangential orientation under the sense mode. Moreover, the two adjacent lumped masses oscillate in anti-phase motion. Corresponding to the operation modes, the two in-phase modes are regarded as the pseudo modes.
Figure 2. FEA shows the operation modes and pseudo modes of CSQMG, comparing to the resonant mode of the HRGs. In the operating modes of CSQMG, also known as anti-phase modes, the four lumped masses vibrate along the radial of the symmetry center in the drive mode, and vibrate along the tangential orientation under the sense mode. Moreover, the two adjacent lumped masses oscillate in anti-phase motion. Corresponding to the operation modes, the two in-phase modes are regarded as the pseudo modes.
Sensors 16 00613 g002
Figure 3. FEA shows the zero-force points of the operation modes. (a) Symmetric generalized force is exerted on the four fixed point during sense mode resonance, and the four placement points of the circular resonator are location points; (b) A mat-shaped bracing structure is designed to eliminate the placement point torque at the four center points of each curved beam, and these points become zero-force points in the drive mode.
Figure 3. FEA shows the zero-force points of the operation modes. (a) Symmetric generalized force is exerted on the four fixed point during sense mode resonance, and the four placement points of the circular resonator are location points; (b) A mat-shaped bracing structure is designed to eliminate the placement point torque at the four center points of each curved beam, and these points become zero-force points in the drive mode.
Sensors 16 00613 g003
Figure 4. Simplified mechanics model of a quarter of the circle resonator. (a) Simulating the operation mode under external forces, and exerting external generalized unit forces FM1, FM2, MM. Numbering the microbeams from No. 1 to No. 9. Identifying the axes of symmetric position with a dotted line; (b) Naming the accompanied internal generalized unit forces FB1, FB2, and MB, when the circle resonator is under external forces.
Figure 4. Simplified mechanics model of a quarter of the circle resonator. (a) Simulating the operation mode under external forces, and exerting external generalized unit forces FM1, FM2, MM. Numbering the microbeams from No. 1 to No. 9. Identifying the axes of symmetric position with a dotted line; (b) Naming the accompanied internal generalized unit forces FB1, FB2, and MB, when the circle resonator is under external forces.
Sensors 16 00613 g004
Figure 5. The bending moment diagram of the circle resonator simplified mechanics model after removing the constraints at the lumped masses symmetric axes. (a) Simulating external forces under the anti-phase drive mode; (b) Simulating external forces under the anti-phase sense mode; (c) Simulating external forces under lumped mass twisting motion coupling with anti-phase sense mode; (d) Assuming internal forces along the tangential orientation; (e) Assuming internal forces along the radial of the symmetric center; (f) Assuming internal forces along the lumped mass twisting motion.
Figure 5. The bending moment diagram of the circle resonator simplified mechanics model after removing the constraints at the lumped masses symmetric axes. (a) Simulating external forces under the anti-phase drive mode; (b) Simulating external forces under the anti-phase sense mode; (c) Simulating external forces under lumped mass twisting motion coupling with anti-phase sense mode; (d) Assuming internal forces along the tangential orientation; (e) Assuming internal forces along the radial of the symmetric center; (f) Assuming internal forces along the lumped mass twisting motion.
Sensors 16 00613 g005
Figure 6. The finite element model of CSQMG with 35 μm chamfer combining two different mesh size. The comb fingers are replaced by evenly distributed mass. Swept meshing with element size less than 0.02 mm is used on the beams.
Figure 6. The finite element model of CSQMG with 35 μm chamfer combining two different mesh size. The comb fingers are replaced by evenly distributed mass. Swept meshing with element size less than 0.02 mm is used on the beams.
Sensors 16 00613 g006
Figure 7. The further simplified model for optimal N-shaped beam design. We presuppose the parameters of the N-shaped beams on both sides of the symmetric axes are exactly the same.
Figure 7. The further simplified model for optimal N-shaped beam design. We presuppose the parameters of the N-shaped beams on both sides of the symmetric axes are exactly the same.
Sensors 16 00613 g007
Figure 8. The relationship of geometric parameters B1,B2,B3,B4 and stiffness differences. (a) The stiffness difference between the two operation modes; (b) The stiffness difference between the in-phase drive mode and the anti-phase drive mode; (c) The stiffness difference between the in-phase sense mode and the anti-phase drive mode.
Figure 8. The relationship of geometric parameters B1,B2,B3,B4 and stiffness differences. (a) The stiffness difference between the two operation modes; (b) The stiffness difference between the in-phase drive mode and the anti-phase drive mode; (c) The stiffness difference between the in-phase sense mode and the anti-phase drive mode.
Sensors 16 00613 g008
Figure 9. (a) The colored elastic strain patterns of Y-shaped coulping beam under anti-phase drive mode (a1); anti-phase sense mode (a2); in-phase drive mode (a3); and in-phase sense mode (a4); (b) The relationship of the Y-shaped beam geometric parameters and four main mode stiffness. Resonant stiffnesses between in-phase modes and anti-phase modes have been distinguished effectively by increasing the Y-shaped beam length and decreasing its width.
Figure 9. (a) The colored elastic strain patterns of Y-shaped coulping beam under anti-phase drive mode (a1); anti-phase sense mode (a2); in-phase drive mode (a3); and in-phase sense mode (a4); (b) The relationship of the Y-shaped beam geometric parameters and four main mode stiffness. Resonant stiffnesses between in-phase modes and anti-phase modes have been distinguished effectively by increasing the Y-shaped beam length and decreasing its width.
Sensors 16 00613 g009
Figure 10. (a) Coupling component along the lumped mass twisting motion exists in the sense mode, which causes the coupling of the drive mode oscillation and sense mode oscillation; (b) The relationship of CSQMG’s geometric parameter L2, L4 and stiffness of twisting motion coupling. The twisting motion coupling can be eliminated by a reasonable allocation.
Figure 10. (a) Coupling component along the lumped mass twisting motion exists in the sense mode, which causes the coupling of the drive mode oscillation and sense mode oscillation; (b) The relationship of CSQMG’s geometric parameter L2, L4 and stiffness of twisting motion coupling. The twisting motion coupling can be eliminated by a reasonable allocation.
Sensors 16 00613 g010
Figure 11. (a) Fabrication process of the CSQMG: (a1) Interconnection forming; (a2) Step etching and ion implantation; (a3) Anodic bonding; (a4) Silicon thining; (a5) Structure releasing; (b) SEM of CSQMG force balance electrodes. The gaps between the adjacent comb fingers are 4 μm and 16 μm; (c) Photograph of the CSQMG die in a metal case ready for experimental characterization.
Figure 11. (a) Fabrication process of the CSQMG: (a1) Interconnection forming; (a2) Step etching and ion implantation; (a3) Anodic bonding; (a4) Silicon thining; (a5) Structure releasing; (b) SEM of CSQMG force balance electrodes. The gaps between the adjacent comb fingers are 4 μm and 16 μm; (c) Photograph of the CSQMG die in a metal case ready for experimental characterization.
Sensors 16 00613 g011
Figure 12. (a) Frequency response curve of the CSQMG. Sine excitations are applied on the drive combs and the signal read on sense combs, then all planar modes are reflected on the curve by vibration coupling; (b) Photograph of the Motion Dynamics 3V-AB 1-Axis Rate-Table (left) and fully packaged CSQMG (right); (c) Measured scale factor of the CSQMG. Above is the fitting result, below is the residual error. The scale factor is 20209.2LSB/(°/s) under ±300°/s measurement range, and nonlinearity is 96 ppm; (d) Measured Allan variance stability of the CSQMG. The white noise level is about 0.7185°/h/√Hz. The Allan variance stability is 0.12°/h (100 s) at room temperature.
Figure 12. (a) Frequency response curve of the CSQMG. Sine excitations are applied on the drive combs and the signal read on sense combs, then all planar modes are reflected on the curve by vibration coupling; (b) Photograph of the Motion Dynamics 3V-AB 1-Axis Rate-Table (left) and fully packaged CSQMG (right); (c) Measured scale factor of the CSQMG. Above is the fitting result, below is the residual error. The scale factor is 20209.2LSB/(°/s) under ±300°/s measurement range, and nonlinearity is 96 ppm; (d) Measured Allan variance stability of the CSQMG. The white noise level is about 0.7185°/h/√Hz. The Allan variance stability is 0.12°/h (100 s) at room temperature.
Sensors 16 00613 g012
Table 1. Torque expressions of xak~xfk.
Table 1. Torque expressions of xak~xfk.
FM1 (a)FM2 (b)MM (c)FB1 (d)FB2 (e)MB (f)
1 x a 1 = F M 1 2 ( q + x 1 ) x b 1 = 0 x c 1 = M M 2 x d 1 = 0 x e 1 = q + x 1 x f 1 = 1
2 x a 2 = F M 1 2 ( q + l 1 ) x b 2 = F M 2 2 x 2 x c 2 = M M 2 x d 2 = x 2 x e 2 = q + l 1 x f 2 = 1
3 x a 3 = F M 1 2 ( q + l 1 + x 3 ) x b 3 = F M 2 2 l 2 x c 3 = M M 2 x d 3 = l 2 x e 3 = q + l 1 + x 3 x f 3 = 1
4 x a 4 = F M 1 2 ( q + l 1 + l 3 ) x b 4 = F M 2 2 ( l 2 x 4 ) x c 4 = M M 2 x d 4 = l 2 x 4 x e 4 = q + l 1 + l 3 x f 4 = 1
5 x a 5 = F M 1 2 ( q + l 1 + l 3 x 5 2 ) + F M 1 2 ( q + l 9 + l 7 x 5 2 ) x b 5 = 0 x c 5 = 0 x d 5 = ( l 2 l 4 x 5 2 ) + ( l 8 l 6 x 5 2 ) x e 5 = 0 x f 5 = 2
6 x a 6 = F M 1 2 ( q + l 9 + l 7 ) x b 6 = F M 2 2 ( l 8 x 6 ) x c 6 = M M 2 x d 6 = ( l 8 x 6 ) x e 4 = q + l 9 + l 7 x f 6 = 1
7 x a 7 = F M 1 2 ( q + l 9 + x 7 ) x b 7 = F M 2 2 l 8 x c 7 = M M 2 x d 7 = l 8 x e 3 = q + l 9 + x 7 x f 7 = 1
8 x a 8 = F M 1 2 ( q + l 9 ) x b 8 = F M 2 2 x 8 x c 8 = M M 2 x d 8 = x 8 x e 2 = q + l 9 x f 8 = 1
9 x a 9 = F M 1 2 ( q + x 9 ) x b 9 = 0 x c 9 = M M 2 x d 9 = 0 x e 1 = q + x 9 x f 9 = 1
Table 2. Stiffness error comparison.
Table 2. Stiffness error comparison.
Design Geometric ParametersBeam NumberLength/μmWidth/μm
1754.570
277634.6
3430.555
484443
511150
ModeStiffness/N/m (Calculation)Stiffness/N/m (ANSYS)Percentage Error
In-phase Drive897.6904.980.82%
Anti-phase Drive834.3834.03−0.03%
Anti-phase Sense781.3866.939.88%
In-Phase Sense667.6703.635.12%
Table 3. Frequency distribution of each mode.
Table 3. Frequency distribution of each mode.
Parameter NameInitial ParametersAfter Mode OptimizationAfter Rotational Coupling Optimization
Beam Length (mm)L10.680.680.68
L20.7550.7960.755
L30.43050.43050.4305
L40.8440.8440.885
L50.10.10.1
Beam Width (mm)B10.080.0720.072
B2 = B40.040.450.45
B30.050.050.05
B50.0240.0270.027
Mode Frequency (Hz)In-Phase Sen581665956558
Pseudo 1620569326923
Anti-phase Dri693073537357
Anti-phase Sen660573647359
Pseudo 2701575557509
In-phase Dri707176457570
Rotational Coupling (1/N)0.0450.0224.2 × 10−4
Table 4. Design parameters of the CSQMG.
Table 4. Design parameters of the CSQMG.
Design ParametersStructureSize With Pads7700 μm × 7700 μm
Thickness80 μm
OscillatorSingle Lumped Mass0.58 mg
Beam Width31 μm~80 μm
Beam Length100 μm~825 μm
Comb FingerComb Drive Gap3.5 μm
Comb Sense Gap4 μm~16 μm
Detection Capacitance3.3 pF in total (dC/dx = 1.2 μF/m)
Mechanical Noise Level0.186°/h/√Hz
Table 5. Measured results of the CSQMG.
Table 5. Measured results of the CSQMG.
Measured ResultsOperating Frequency (After Tuning)6.8 kHz
Q factor(20 °C & ≤ 30 Pa)8500/500 on drive/sense mode
Scale Factor20,209.2 LSB/(°/s)
Scale Factor Nonlinearity96 ppm
Measurement Range±300°/s
White Noise Level0.7185°/h/√Hz
Allan Variance Stability0.12°/h (@100 s)

Share and Cite

MDPI and ACS Style

Zhang, T.; Zhou, B.; Yin, P.; Chen, Z.; Zhang, R. Optimal Design of a Center Support Quadruple Mass Gyroscope (CSQMG). Sensors 2016, 16, 613. https://doi.org/10.3390/s16050613

AMA Style

Zhang T, Zhou B, Yin P, Chen Z, Zhang R. Optimal Design of a Center Support Quadruple Mass Gyroscope (CSQMG). Sensors. 2016; 16(5):613. https://doi.org/10.3390/s16050613

Chicago/Turabian Style

Zhang, Tian, Bin Zhou, Peng Yin, Zhiyong Chen, and Rong Zhang. 2016. "Optimal Design of a Center Support Quadruple Mass Gyroscope (CSQMG)" Sensors 16, no. 5: 613. https://doi.org/10.3390/s16050613

APA Style

Zhang, T., Zhou, B., Yin, P., Chen, Z., & Zhang, R. (2016). Optimal Design of a Center Support Quadruple Mass Gyroscope (CSQMG). Sensors, 16(5), 613. https://doi.org/10.3390/s16050613

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