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 x
ak~x
fk (
k ∈ [1, 9] ∩
k ∈
N) are listed in
Table 1 based on
Figure 5.
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):
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:
where A
0 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 k
B 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, A
0, into Equation (3), and the mechanical-thermal noise Ω
th is calculated to be 0.186°/h/√Hz.
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.