Next Article in Journal
Recent Advances and Future Perspectives on Microfluidic Mix-and-Jet Sample Delivery Devices
Next Article in Special Issue
Integrated Resonant Micro/Nano Gravimetric Sensors for Bio/Chemical Detection in Air and Liquid
Previous Article in Journal
Effect of Nanopores on Mechanical Properties of the Shape Memory Alloy
Previous Article in Special Issue
Control of Spring Softening and Hardening in the Squared Daisy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Optimization of the Resonator in a Resonant Accelerometer Based on Mode and Frequency Analysis

School of Mechanical Electronic & Information Engineering, China University of Mining and Technology-Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Micromachines 2021, 12(5), 530; https://doi.org/10.3390/mi12050530
Submission received: 1 April 2021 / Revised: 2 May 2021 / Accepted: 3 May 2021 / Published: 7 May 2021
(This article belongs to the Special Issue Integrated MEMS Resonators)

Abstract

:
This study aims to develop methods to design and optimize the resonator in a resonant accelerometer based on mode and frequency analysis. First, according to the working principle of a resonant accelerometer, the resonator is divided into three parts: beam I, beam II, and beam III. Using Hamilton’s principle, the undamped dynamic control equation and the ordinary differential dynamic equation of the resonant beam are obtained. Moreover, the structural parameters of the accelerometer are designed and optimized by using resonator mode and frequency analysis, then using finite element simulation to verify it. Finally, 1 g acceleration tumbling experiments are built to verify the feasibility of the proposed design and optimization method. The experimental results demonstrate that the proposed accelerometer has a sensitivity of 98 Hz/g, a resolution of 0.917 mg, and a bias stability of 1.323 mg/h. The research findings suggest that according to the resonator mode and frequency analysis, the values of the resonator structural parameters are determined so that the working mode of the resonator is far away from the interference mode and avoids resonance points effectively. The research results are expected to be beneficial for a practical resonant sensor design.

1. Introduction

A resonant accelerometer converts acceleration directly into a frequency, it benefits from the advantages of high measurement accuracy, strong anti-disturbance capability, good stability, wide dynamic range, and quasi-digital output [1,2,3,4,5]. A resonator is a device designed to vibrate at a specified frequency. The specific frequency is one where the system retains input energy with minimum loss, which is one of its natural frequencies [6,7]. A resonator is the basis for combining many devices, such as radio frequency filters and resonant sensors [8,9,10]. In a resonant accelerometer, the resonator is the core sensitive element [11,12]. Thus, the mode and frequency of the resonator affect the performance of the resonant accelerometer greatly [13,14].
In their study on resonator mode, Yang and Huang proposed an analytical model to obtain a higher order mode in a rectangular resonator [15]. Phani et al. studied the modal coupling of the z-axis frame vibration gyro [16]. Guo et al. calculated the modal frequency and quality factor for the square resonant cavity [17]. The effect of dispersive modal coupling to tune the frequency difference between two coupled modes in a wide range was studied by Lu et al. [18]. In terms of resonator frequency, Hassanpour et al. studied the influence of axial force, position, mass ratio, and radius of rotation on the natural frequency and resonator mode [19]. A new method was proposed to improve the resonator’s frequency stability by Du et al. [20]. Huang et al. found that the displacement due to warpage brings the natural frequency closer to the real natural frequency of the resonator [21]. Guzman et al. developed a SiC bridge resonator with an unprecedentedly large tunable frequency range [22]. The aforementioned studies mainly investigated the frequency error of the coupled mode and its effect on the resonator. However, none of the studies examined how the mode influences the structural parameters and sensitivity of the resonator explicitly.
To design the structure parameters of the resonator, Shi et al. studied the gap sizes and the excitation parameters and their effect on the resonator performance by using the uncertainty method [23]. Zotov et al. proposed a dynamic equilibrium method to maximize the quality factor in the defect structure of the resonator [24]. Candler et al. demonstrated the effect of the curved resonant beam gap on the quality factor [25]. Zhou et al. used a multi-objective optimization method to design the resonator structure parameters [26]. As for the resonator sensitivity, the influence of the structure parameters on resonator sensitivity was studied by Shi et al. [27]. Xu et al. designed a new structure with a suspended piezoresistive beam (SPB) to improve sensitivity without sacrificing the resonant frequency [28]. Du et al. used the twinborn diaphragms to improve the resonator sensitivity [29]. Zhao et al. used three weakly coupled resonators to improve sensitivity to stiffness changes [30]. The aforementioned studies mainly focused on the influence of the structure parameters on the resonator sensitivity but did not propose a new structure parameter optimization method to improve the sensitivity.
Equating the resonator to a single resonant beam is an effective way to study the dynamic characteristics of the resonator in resonant sensors. Our previous work has established a dynamic resonant beam model and analyzed its dynamic characteristics [31]. In addition, the stability analysis method and vibration excitation parameter characteristics have been proposed based on a single resonant beam model [32,33].
This simplified model that equates the resonator to a single resonant beam is convenient for analysis, and its force–frequency characteristics can be found easily. Unfortunately, there is an error with the real working conditions. In the present study, the displacement model of the resonant beam is established to reduce this error, and the Hamilton principle is used to establish a partial differential dynamic control equation. The linear mode and frequency are solved, which is the quantitative analysis being used to design the resonator, and a series of experiments being carried out to prove the feasibility of the theoretical analysis method.

2. Theory of the Resonator

2.1. Working Mechanism

The resonant accelerometer structure is shown in Figure 1, consisting mainly of an external frame, a proof mass, two resonators, and four elastic beams. Each of the resonators is actuated at a resonance using drives in this particular implementation. The external acceleration applied to the proof mass along the sensitive axis results in an inertial force transmitted axially to the resonator. The applied axial force results in a shift in the resonant frequency of the resonator. This effect is the same as tuning a guitar string to resonate at different frequencies by varying the tension in the string. The output of the accelerometer is the difference between the output frequency of the two resonators.
In the resonator working mode, the local vibration mode of its upper and lower beam is the shape of the reverse half wave. As shown in Figure 2, the smaller the displacement of the intersection of the structure and the dotted line, the higher the energy of the excitation that is transferred to the working mode. Thus, in order to facilitate the analysis in the working mode, the dynamic characteristics of half of the resonator structure can be analyzed based on the assumption that the intersection points A and B are static.
The resonator model is simplified to be made up of beam I, beam II, and beam III, as shown in Figure 2. As a function of the applied acceleration, the response of the resonant beam to the external axial force is to bend and twist. The main response of its lower order vibration is to bend. The effects of the changing external force result in a bending vibration in the plane of the three beams. In the simplified resonator model, one end of each resonant beam is fixed and the other end is free; beams I and III are kept parallel; and the length, width, and thickness of the beams are variable parameters that can be used to describe different resonator states.

2.2. Vibrational Mode and Frequency

Beam III is parallel to beam I prior to the deformation of the resonator structure. As shown in Figure 3, o 1 x 1 y 1 is the local coordinate of beam I, which is consolidated at the left end of beam I. o 2 x 2 y 2 is the local coordinate of beam II, which is consolidated at the left end of beam II. o 3 x 3 y 3 is the local coordinate of beam III, which is consolidated at the left end of beam III. In the local coordinate o i x i y i , the lateral displacement of the beam is v i and the axial displacement is u i . In the global coordinate o x y , the displacement of the resonant beam is r i .
The kinetic energy expression of each section of the resonator is as follows:
T i = 0 L i A i 1 2 ρ i | r ˙ i | 2 d A i d x i = 0 L i h T i d x i , i { 1 , 2 , 3 } ,
where the dots indicate derivation with time t, L i is the length of each beam, ρ i is the density of each beam, A i is the cross-sectional area of each beam.
The potential energy expression of each section of the resonator is as follows:
V i = 0 L i A i ε i σ i 2 d A i d x i = 0 L i h V i d x i , i { 1 , 2 , 3 } ,
where ε i is the axial line strain of each beam, σ i is the axial positive stress of each beam. The expressions are as follows:
ε i = u i y i v i , σ i = E ε i ,
Using the Hamilton principle,
t 1 t 2 i = 1 3 ( δ T i δ V i ) d t = 0 ,
The dynamic control equation can be expressed as follows by the energy density function:
2 h T i u ˙ i t 2 h V i u i x i = 0 , 2 h T i v ˙ i t + 3 h V i v i x i 2 = 0 , i { 1 , 2 , 3 } ,
In the equation, the six unknown functions represent the lateral and axial displacements of the three resonant beams in the resonator. Expressions of the 18 boundary conditions corresponding to the motion control in Equation (5) are described as follows:
[ i = 1 3 0 L i 2 h T i u ˙ 1 ( L 1 ) t d x i h V 1 u 1 | x 1 = L 1 ] δ u 1 ( L 1 ) = 0 , h V i u i δ u i | x i = 0 = 0 ,
[ h V 2 u 2 | x 2 = L 2 + 0 L 3 2 h T 3 u ˙ 2 ( L 2 ) t d x 3 ] δ u 2 ( L 2 ) = 0 , h V 3 u 3 δ u 3 | x 3 = L 3 = 0 ,
[ 2 h V 1 v 1 x 1 | x 1 = L 1 i = 1 3 0 L i 2 h T i v ˙ 1 ( L 1 ) t d x i ] δ v 1 ( L 1 ) = 0 , 2 h V i v i x i δ v i | x i = 0 = 0 ,
[ i = 1 3 0 L i 2 h T i v ˙ 1 ( L 1 ) t d x i + h V 1 v 1 δ v 1 | x 1 = L 1 ] δ v 1 ( L 1 ) = 0 , h V i v i δ v i | x i = 0 = 0 ,
[ 2 h V 2 v 2 x 2 | x 2 = L 2 0 L 3 2 h T 3 v ˙ 2 ( L 2 ) t d x 3 ] δ v 2 ( L 2 ) = 0 , 2 h V 3 v 3 x 3 δ v 3 | x 3 = L 3 = 0 ,
[ h V 2 v 2 δ v 2 | x 2 = L 2 + 0 L 3 2 h T 3 v ˙ 2 ( L 2 ) t d x 3 ] δ v 2 ( L 2 ) = 0 , h V 3 v 3 δ v 3 | x 3 = L 3 = 0 , i { 1 , 2 , 3 } ,
The axial displacement in the dynamic control in Equation (5) and its boundary condition expressed in Equations (6)–(11) are ignored considering the removal of undamped free vibration and the external excitation terms. The linear undamped homogeneous dynamic control equation is derived from Equation (5) as follows:
v ¨ 1 + E 1 J 1 ρ 1 A 1 v 1 ( 4 ) = 0 , v ¨ 2 + E 2 J 2 ρ 2 A 2 v 2 ( 4 ) + cos ( α ) v ¨ 1 ( L 1 ) + x 2 v ¨ 1 ( L 1 ) = 0 ,
v ¨ 3 + E 3 J 3 ρ 3 A 3 v 3 ( 4 ) + cos ( α ) [ L 2 v ¨ 1 ( L 1 ) + v ¨ 2 ( L 2 ) ] + x 3 [ v ¨ 1 ( L 1 ) + v ¨ 2 ( L 2 ) ] + v ¨ 1 ( L 1 ) = 0 ,
where the numbers with small brackets represent a partial derivative to spatial coordinates, the apostrophe represents a partial derivative to the spatial coordinates of the function, and the double dot represents the derivative to the time.
Twelve linear boundary conditions of the three unknown functions in Equations (12) and (13) are expressed as follows:
h V 1 v 1 δ v 1 | x 1 = 0 = 0 , [ h V 1 v 1 | x 1 = L 1 + 0 L 2 2 h T 2 v ˙ 1 ( L 1 ) t d x 2 + 0 L 3 2 h T 3 v ˙ 1 ( L 1 ) t d x 3 ] δ v 1 ( L 1 ) = 0 ,
2 h V 1 v 1 x 1 δ v 1 | x 1 = 0 = 0 , [ 2 h V 1 v 1 x 1 | x 1 = L 1 0 L 2 2 h T 2 v ˙ 1 ( L 1 ) t d x 2 0 L 3 2 h T 3 v ˙ 1 ( L 1 ) t d x 3 ] δ v 1 ( L 1 ) = 0 ,
h V 2 v 2 δ v 2 | x 2 = 0 = 0 , [ h V 2 v 2 | x 2 = L 2 + 0 L 3 2 h T 3 v ˙ 2 ( L 2 ) t d x 3 ] δ v 2 ( L 2 ) = 0 ,
2 h V 2 v 2 x 2 δ v 2 | x 2 = 0 = 0 , [ 2 h V 2 v 2 x 2 | x 2 = L 2 0 L 3 2 h T 3 v ˙ 2 ( L 2 ) t d x 3 ] δ v 2 ( L 2 ) = 0 ,
h V 3 v 3 δ v 3 | x 3 = 0 = 0 , 2 h V 3 v 3 x 3 δ v 3 | x 3 = 0 = 0 , h V 3 v 3 δ v 3 | x 3 = L 3 = 0 , 2 h V 3 v 3 x 3 δ v 3 | x 3 = L 3 = 0 ,
The following transformation is introduced:
v 2 = v 2 + cos ( α ) v 1 ( L 1 ) + x 2 v 1 ( L 1 ) ,
v 3 = v 3 + v 1 ( L 1 ) + cos ( α ) [ L 2 v 1 ( L 1 ) + v 2 ( L 2 ) ] + x 3 [ v 1 ( L 1 ) + v 2 ( L 2 ) ] ,
Equations (12) and (13) are substituted into Equation (5) to obtain the following:
v ¨ 1 + E 1 J 1 ρ 1 A 1 v 1 ( 4 ) = 0 , v ¨ 2 + E 2 J 2 ρ 2 A 2 v 2 ( 4 ) = 0 , v ¨ 3 + E 3 J 3 ρ 3 A 3 v 3 ( 4 ) = 0 ,
The solution to Equation (21) can be expressed as follows:
φ i ( x i ) = C i 1 cos ( λ i x i ) + C i 2 sin ( λ i x i ) + C i 3 c h ( λ i x i ) + C i 4 s h ( λ i x i ) , i { 1 , 2 , 3 } ,
Using Equations (14)–(18) and (22), the system frequency characteristic equation can be obtained; then, the frequency eigenvalue can be obtained and the linear frequency can be obtained for each order. By substituting the obtained eigenvalues into Equation (22) and subsequently combining Equations (12) and (13), a modal function can be obtained.
According to the solution of Equation (5) and its corresponding boundary conditions Equations (6)–(11), the two-dimensional (2D) six-degrees-of-freedom beam element is used to divide the resonator. The degrees of freedom of the unit are as follows:
δ i = [ u ς i v ς i ϕ ς i u η i v η i ϕ η i ] T ,
where u ς i , v ς i and ϕ ς i represent the axial degree of freedom, lateral degree of freedom, and section angle of the node ς for unit i, respectively. u η i , v η i , and ϕ η i represent the axial degree of freedom, lateral degree of freedom, and section angle of node η for unit i, respectively. Superscript T stands for transpose.
The displacement of unit i can be written as follows:
[ u i v i ] = N i δ i = N i [ u ς i v ς i ϕ ς i u η i v η i ϕ η i ] T , i { 1 , 2 , 3 , , n } ,
where N i is a 2 × 6 matrix representing the unit functions of the displacement shape, and its expression is as follows:
N i = [ 1 x l i 0 0 ( x l i ) 0 0 0 1 3 ( x l i ) 2 + 2 ( x l i ) 3 ( x l i ) 2 ( x l i ) 2 + ( x l i ) 3 0 3 ( x l i ) 2 2 ( x l i ) 3 ( x l i ) 2 + ( x l i ) 3 ]
The kinetic energy and potential energy of unit i are as follows:
δ T i = [ u ˙ i v ˙ i ] [ δ u ˙ i δ v ˙ i ] d m i = δ ˙ i T ( Π i T N i T ρ i N i d m i Π i ) δ δ ˙ i ,
δ V i = ε i σ i δ ε i d v i = δ i T ( Π i T B i T σ i B i d v i Π i ) δ δ i ,
where the transformation matrix between the local coordinates of the unit and the global coordinate of the resonator is Π i , and where the shape function of the axial strain is B, the expressions are as follows:
Π 1 = [ 0 1 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 1 ] , Π 2 = [ 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 ] , Π 3 = [ 0 1 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 1 ] B i = [ 1 l i y ¯ ( 6 l i 2 12 l i 3 x ¯ ) y ¯ ( 4 l i 6 l i 2 x ¯ ) 1 l i y ¯ ( 12 l i 3 x ¯ 6 l i 2 ) y ¯ ( 2 l i 6 l i 2 x ¯ ) ]
where ( x ¯ , y ¯ ) represents the local coordinates of unit i. To obtain the dynamic control equation of all degrees of freedom of the resonator, unit i can be represented by the positioning matrix and the column vector of the total degrees of freedom; the expression is as follows:
δ i = L i P = L i [ u 0 v 0 φ 0 u n v n φ n ] T ,
where P is the global generalized coordinate array and L i is the positioning matrix, the expressions are as follows:
p = [ u 0 v 0 φ 0 u n v n φ n ] T ,
1 ς η n + 1 L i = [ 0 0 I 3 × 3 0 0 0 0 0 0 0 0 0 0 I 3 × 3 0 0 ] 6 × 3 ( n + 1 ) ,
Substituting Equations (23)–(31) into the Hamiltonian equation,
t 1 t 2 i = 1 n ( δ T i δ V i ) d t = 0 ,
The dynamic control equation for the global generalized coordinate array is obtained as follows:
Μ p ¨ + Κ p = 0 ,
where M and K are the global mass matrix and global stiffness matrix, the expressions are as follows:
M = i = 1 n L i T Π i T N i T N i d m i Π i L i , K = i = 1 n L i T Π i T B i T σ i B i d v i Π i L i ,
When the physical parameters of the resonator, such as geometry and material, are specified, the n-order natural frequency and vibration mode of the resonator can be calculated by Equation (34).

3. Design of the Resonator

Resonator materials are assumed to be silicon or steel. The relevant parameters for silicon and steel materials are shown in Table 1. The structural parameters for the resonant beam are shown in Figure 4. L1: L2: L3 = 1: PL: 1, B1: B2: B3 = 1: PB: 1, where PL is the length ratio of the three beams, PB is the width ratio of the three beams, and the physical parameters of the three beams are the same. The relationship between frequency and the length, as well as the width, can be obtained by the numerical model Equation (34). The frequency at any point can be compared by using the theoretical model Equations (14)–(18) and (22).
As shown in Figure 5 and Figure 6, the solid lines were obtained by the numerical model Equation (34), and the asterisks were derived from theoretical model Equations (14)–(18) and (22). Figure 5 and Figure 6 show the relationship between the resonant frequencies of the first four resonator modes and the fixed-length ratio PL or the width ratio PB under the condition of silicon or steel materials. As shown in Figure 5, when the width ratio of the three beams is fixed in 1, the resonant frequency of each mode and the frequency difference between the modes decrease with the increase of length ratio. As shown in Figure 6, when the length ratio of the three beams is fixed in 6, the natural frequency of each mode and the frequency difference between the modes increase with the increase of width ratio.
Thus, the working mode of the resonator can be quantitatively selected to keep it away from the interference mode. If the second-order mode is selected as the working mode of the resonator, the working mode of the corresponding single beam is the first-order mode. As demonstrated by the results of the simulation and analysis shown in Figure 5 and Figure 6, PL between 5 and 9 was selected, and PB was selected larger than 0.5, which can effectively keep the working mode away from the interference mode.
Depending on the relationship between the higher order mode frequency and the resonator parameters, it was found that there is a multiple relationship between the higher order mode frequency and a certain resonator parameter. When the structural parameter was adjusted to make the second-order mode frequency close to two or three times the first-order mode frequency, the internal resonant conditions between the first two modes were 1:2 or 1:3. Figure 7 and Figure 8 show the internal resonant frequency of the resonator with the silicon and steel material, which were assumed to have a different fixed-length and fixed-width ratios. As shown in Figure 7, under the fixed-width ratio of 1, the second-order mode resonant frequency and the double frequency of the first-order mode resonant frequency intersect at point P1, where the PL is approximately 6.5, indicating that there is a 2:1 internal resonance between the second-order mode and the first-order mode. Under the fixed-width ratio of 1, the third-order mode resonant frequency and the double second-order mode resonant frequency intersect at point P2, where the PL is approximately 9.3, indicating that the resonator parameters have a 2:1 internal resonance between the third-order mode and the second-order mode. Likewise, as can be seen in Figure 8, with the fixed-length ratio of 7, the second-order mode resonant frequency and the double frequency of the first-order mode resonant frequency intersect at the point where the PB is approximately 0.7. It can be seen from the results of the analysis in Figure 7 and Figure 8 that the PL between 6.5 and 9.3 is selected, and the PB smaller than 0.7 is selected, which can avoid the internal resonance points effectively.
Based on the above analysis, physical parameters have little influence on the performance of the accelerometer; hence, the structural parameters of the resonator should be considered essential when optimizing the resonator. According to the mode and frequency analysis, when optimizing resonator structural parameters, the values of PL and PB are determined by theoretical calculation and numerical simulation so that the working mode of the resonator is far away from the interference mode. In addition, the selected PL and PB values of the resonator were verified to avoid the resonance points effectively.

4. Optimization of the Resonator

4.1. Boundary Conditions

There are two kinds of boundary conditions of a resonator structure, one is single ended and the other is double ended. Figure 9 shows the first-order mode to fourth-order mode of the resonator structure in the case of single-ended fixed support (fixed at the left end, free at the right end). Among them, the second-order mode is the required working mode, but the free end is displaced, so that the excitation energy dissipates to the sensitive structure. When the axial acceleration is perceived, it will form a coupling with the inertia force, which affects the measurement accuracy.
Figure 10 shows the first-order mode to fourth-order mode of the resonator structure in the case of double-ended fixed support, and the first-order mode is the working mode. Under the same geometrical parameters, compared with the boundary of single-ended fixed support, the resonant beam with double-ended fixed support is conducive to signal detection. There is no displacement at both ends of the structure, and the energy is concentrated in the middle section, which is in line with the modal shape and design goal required by the ideal resonator. The structural diagram of the resonator is shown in Figure 11. L1, L2, and L3 represent the length of beam I, beam II, and beam III; B1, B2, and B3 represent the width of beam I, beam II, and beam III; h represents the thickness of the resonator; L4 and b represent the length and width of the beam root, respectively.

4.2. Structural Parameters

Based on the above mode and frequency analysis, the PL between 5 and 9 was selected, and the PB larger than 0.5 was selected, which can keep the working mode away from the interference mode effectively. The PL between 6.5 and 9.3 was selected, and PB smaller than 0.7 was selected, which can avoid an internal resonance point effectively. Therefore, the length ratio PL was determined to be 8 and the width ratio was determined to be 0.6.
Table 2 shows the resonant beam and beam root parameters and the corresponding stiffness of the three sets of resonators.
According to the above three groups of structural parameters, the first four modes of each structure were obtained by simulation as shown in Figure 12, Figure 13 and Figure 14. The relationship between the degree of working frequency away from the interference frequency and the resonator root stiffness under different resonator root parameters are shown in Figure 15. The working mode frequency of the resonator was 3901.4 Hz according to f = 3.75 L 2 E I ρ b h . As shown in Figure 15, the beam root parameters of the resonator were selected from beam root 3, in which the working mode frequencies of the resonator are far away from the interference mode frequencies distinctly.
The whole structure of the resonant accelerometer was simulated according to the designed parameters, and the results are as shown in Figure 16, which indicates that the working mode of the accelerometer was far away from the interference mode. The simulated results reflect the effectiveness of the designed resonator parameters.

5. Experiments and Discussion

The designed parameters of the resonant accelerometer are shown in Table 3. The resonant accelerometer sample was fabricated from steel materials. In the actual fabrication process, in order to reduce the impact of the fabrication technology precision on the performance of the accelerometer, the micromachining method was adapted to achieve the precise control of the thickness, side wall perpendicularity, smoothness and other process parameters, to ensure the structural dimension accuracy of the accelerometer and improve the sensitivity of the accelerometer. The vibration characteristics of the resonant accelerometer was tested in a series of 1 g acceleration tumbling experiments.
During the accelerometer tumbling experiments, the input acceleration changed according to the sine rule. The corresponding output should also change according to the sine rule. Due to the multiple disturbances, the output was a periodic function and did not change entirely according to the sine rule. The actual output periodic function was broken down by the Fourier series, and the coefficients of the Fourier series could be converted into the coefficients of the output model.
The output model equation can be expressed in the following form:
U = K 0 + K 1 A i + K 2 A i 2 + K 3 A i 3 ,
where U is the output of the accelerometer, K 0 is the partial value, K 1 is the sensitivity, K 2 is the nonlinear coefficient of the second order, and K 3 is the nonlinear coefficient of the third order.
The experimental device shown in Figure 17 is composed of an accelerometer, rate table, and signal sampling device. Initially, vibration isolation and anti-tilt equipment were used to improve the experimental accuracy. The accelerometer was fixed on the rate table, and the observed resonant frequency of the resonator was measured using the signal processing device under the designed rotation angle of the rate table. When there were no acceleration and rotation angles applied to the device, two resonators were driven in to vibrate in their natural resonant frequency. Subsequently, an external rotation angle was applied to the device resulting in a frequency modulated oscillator output that was collected and processed by a signal processing device. Theoretically, the linear output model of the accelerometer can be achieved by measuring the modulated frequency output when the accelerometer is at 1 g. In fact, due to the need for accurate measurement, some special rotation angles were selected for multi-point measurement. In the experiments, the accelerometer device was tested on a rate table providing 36 rotation angles between 0° and 360°. The rate table started from the mechanical zero of the accelerometers and then it was turned clockwise and one way precisely to 0°, 10°, 20°, 30°, 40°, 45°, 50° …, 120°, 130°, 135°, 140° …, 220°, 225°, 230°, …, 310°, 315°, 320°, …, 360°. For each angle within 20 s, 5 output frequencies were selected, then the average of the 5 output frequencies were taken as the output. Finally, the rate table was turned counterclockwise and one way precisely to 360°, 350°, 340°, …, 20°, 10°, 0°, the average of 5 frequencies was also selected out of each angle within 20 s. To improve the accuracy of the measurement, the time of residence at each angle position was equal. The above experiments were repeated four times, and the average of the four measurements was obtained as the output.
The output model of the accelerometer can be obtained by processing the test data as follows:
U = 1.553 + 98.229 A i + 1.889 A i 2 + 1.978 A i 3 ,
Figure 18 shows that there is a linear relationship between the output frequencies of the tested accelerometer and the acceleration. Figure 19 shows the sensitivity testing curve of the resonant accelerometer sample. Equation (36) shows that the sensitivity of the resonant accelerometer was 98 Hz/g. At the same time, the accelerometer was installed on a two-axis rotating table surface. When the rotating table was leveled, the surface was locked after turning a small angle around the tilt axis. The rotation axis rotated and modulated at a certain angular velocity to subdivide the gravitational acceleration component. After sampling the output of the accelerometer, FFT low-pass filtering was performed. The resolution of the accelerometer was estimated to be 0.917 mg. In addition, the accelerometer was installed under a normal temperature condition; the power was turned on and off repeatedly for multiple sampling measurements after the installation stress was fully released; the bias stability of the accelerometer was estimated to be 1.323 mg/h. The above test results demonstrated that the accelerometer designed according to the proposed mode and frequency analysis method is feasible. In the future work, we will analyze the influence of error sources on the sensitivity of the resonant accelerometer and the solutions to further enhance the sensitivity. We will also consider a new acceleration sensitive mechanism and improve the sensitivity of the accelerometer to achieve an ultra-sensitive accelerometer.

6. Conclusions

  • According to the working principle of a resonant accelerometer, the resonator was divided into beam I, beam II, and beam III. The shape and displacement of the resonator was described in detail. Using Hamilton’s principle, the undamped dynamic control equation and the corresponding boundary conditions of the resonant beam were obtained. Furthermore, the ordinary differential dynamic equation of the resonant beam with n degrees of freedom was obtained to verify the correctness of the theoretical analysis.
  • The structural parameters of the accelerometer were designed and optimized by using resonator mode and frequency analysis; then, they were verified by using a finite element simulation. As the simulation results demonstrated, the working mode of the resonant accelerometer is far away from the interference mode and avoids resonance points effectively.
  • From the 1 g acceleration tumbling experiments results of the resonant accelerometer sample, the sensitivity is 98 Hz/g, the resolution is 0.917 mg, and the bias stability is 1.323 mg/h. This shows that the proposed design and optimization method of the resonant accelerometer based on resonator mode and frequency analysis is effective and feasible.
  • This research provides an important and novel method for the design and optimization of a resonant accelerometer, which may inspire the use of resonant accelerometers for a variety of motion sensing applications ranging from inertial navigation to vibration monitoring. We will consider a new acceleration sensitive mechanism to achieve an ultra-sensitive accelerometer in the future.

Author Contributions

Conceptualization, Y.L. and B.J.; methodology, Y.L.; software, B.J.; validation, Y.L., B.J. and M.Z.; formal analysis, Y.L.; investigation, M.Z.; resources, Y.L.; data curation, B.J.; writing—original draft preparation, B.J.; writing—review and editing, B.J.; visualization, M.Z.; supervision, M.Z.; project administration, F.Y.; funding acquisition, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China with grant No. 61973308 and the Fundamental Research Funds for the Central Universities of China with grant No. 800015Z1160.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, J.; Su, Y.; Shi, Q.; Qiu, A.P. Microelectromechanical resonant accelerometer designed with a high sensitivity. Sensors 2015, 15, 30293–30310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Yang, B.; Wang, X.; Dai, B.; Liu, X. A new z-axis resonant micro-accelerometer based on electrostatic stiffness. Sensors 2015, 15, 687–702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, B.; Zhao, Y.; Li, C.; Cheng, R.; Sun, D.; Wang, S. A differential resonant accelerometer with low cross-interference and temperature drift. Sensors 2017, 17, 178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mustafazade, A.; Pandit, M.; Zhao, C.; Sobreviela, G.; Du, Z.; Steinmann, P.; Seshia, A.A. A vibrating beam MEMS accelerometer for gravity and seismic measurements. Sci. Rep. 2020, 10, 10415. [Google Scholar] [CrossRef]
  5. Han, C.; Li, C.; Zhao, Y.; Li, B. High-Stability Quartz Resonant Accelerometer with Micro-Leverages. J. Microelectromech. Syst. 2021, 30, 184–192. [Google Scholar] [CrossRef]
  6. Han, G.; Fu, X.; Xu, L. Free vibration for a micro resonant pressure sensor with cross-type resonator. Adv. Mech. Eng. 2021, 13, 1687814020987336. [Google Scholar] [CrossRef]
  7. Jiang, B.; Huang, S.; Zhang, J.; Su, Y. Analysis of Frequency Drift of Silicon MEMS Resonator with Temperature. Micromachines 2021, 12, 26. [Google Scholar] [CrossRef]
  8. Carr, A.R.; Patel, Y.H.; Neff, C.R.; Charkhabi, S.; Kallmyer, N.E.; Angus, H.F.; Reuel, N.F. Sweat monitoring beneath garments using passive, wireless resonant sensors interfaced with laser-ablated microfluidics. Npj Digit. Med. 2020, 3, 62. [Google Scholar] [CrossRef]
  9. Kochhar, A.; Galanko, M.E.; Soliman, M.; Abdelsalam, H.; Colombo, L.; Lin, Y.C.; Fedder, G.K. Resonant microelectromechanical receiver. J. Microelectromech. Syst. 2019, 28, 327–343. [Google Scholar] [CrossRef]
  10. Liang, W.; Ilchenko, V.S.; Savchenkov, A.A.; Dale, E.; Eliyahu, D.; Matsko, A.B.; Maleki, L. Resonant microphotonic gyroscope. Optica 2017, 4, 114–117. [Google Scholar] [CrossRef]
  11. Xu, L.; Wang, S.; Jiang, Z.; Wei, X. Programmable synchronization enhanced MEMS resonant accelerometer. Microsyst. Nanoeng. 2020, 6, 63. [Google Scholar] [CrossRef]
  12. Wang, S.; Zhu, W.; Shen, Y.; Ren, J.; Gu, H.; Wei, X. Temperature compensation for MEMS resonant accelerometer based on genetic algorithm optimized backpropagation neural network. Sens. Actuator A Phys. 2020, 316, 112393. [Google Scholar] [CrossRef]
  13. OlssonIII, R.H.; Wojciechowski, K.E.; Baker, M.S.; Tuck, M.R.; Fleming, J.G. Post-CMOS-compatible aluminum nitride resonant MEMS accelerometers. J. Microelectromech. Syst. 2009, 18, 671–678. [Google Scholar] [CrossRef]
  14. Azgin, K.; Valdevit, L. The effects of tine coupling and geometrical imperfections on the response of DETF resonators. J. Micromech. Microeng. 2013, 23, 125011. [Google Scholar] [CrossRef]
  15. Yang, Y.D.; Huang, Y.Z. Mode analysis and Q-factor enhancement due to mode coupling in rectangular resonators. IEEE J. Quantum Electron. 2007, 43, 497–502. [Google Scholar] [CrossRef]
  16. Phani, A.S.; Seshia, A.A.; Palaniapan, M.; Howe, R.T.; Yasaitis, J. Modal coupling in micromechanical vibratory rate gyroscopes. IEEE Sens. J. 2006, 6, 1144–1152. [Google Scholar] [CrossRef]
  17. Guo, W.H.; Huang, Y.Z.; Lu, Q.Y.; Yu, L.J. Modes in square resonators. IEEE J. Quantum Electron. 2003, 39, 1563–1566. [Google Scholar]
  18. Lu, K.; Li, Q.; Zhou, X.; Song, G.; Wu, K.; Zhuo, M.; Xiao, D. Modal Coupling Effect in a Novel Nonlinear Micromechanical Resonator. Micromachines 2020, 11, 472. [Google Scholar] [CrossRef] [PubMed]
  19. Hassanpour, P.A.; Cleghorn, W.L.; Esmailzadeh, E.; Mills, J.K. Vibration analysis of micro-machined beam-type resonators. J. Sound Vibr. 2007, 308, 287–301. [Google Scholar] [CrossRef]
  20. Du, J.; Wei, X.; Ren, J.; Wang, J.; Huan, R. Micromechanical vibration absorber for frequency stability improvement of DETF oscillator. J. Micromech. Microeng. 2019, 29, 045005. [Google Scholar] [CrossRef]
  21. Huang, B.S.; Chang-Chien, W.T.; Hsieh, F.H.; Chou, Y.F.; Chang, C.O. Analysis of the natural frequency of a quartz double-end tuning fork with a new deformation model. J. Micromech. Microeng. 2016, 26, 065006. [Google Scholar] [CrossRef]
  22. Guzman, P.; Dinh, T.; Phan, H.P.; Joy, A.P.; Qamar, A.; Bahreyni, B.; Dao, D.V. Highly-doped SiC resonator with ultra-large tuning frequency range by Joule heating effect. Mater. Des. 2020, 194, 108922. [Google Scholar] [CrossRef]
  23. Shi, H.; Fan, S.; Zhang, Y. Design and optimization of slit-resonant beam in a MEMS pressure sensor based on uncertainty analysis. Microsyst. Technol. 2017, 23, 5545–5559. [Google Scholar] [CrossRef]
  24. Zotov, S.A.; Simon, B.R.; Prikhodko, I.P.; Trusov, A.A.; Shkel, A.M. Quality factor maximization through dynamic balancing of tuning fork resonator. IEEE Sens. J. 2014, 14, 2706–2714. [Google Scholar] [CrossRef]
  25. Candler, R.N.; Duwel, A.; Varghese, M.; Chandorkar, S.A.; Hopcroft, M.A.; Park, W.T.; Kenny, T.W. Impact of geometry on thermoelastic dissipation in micromechanical resonant beams. J. Microelectromech. Syst. 2006, 15, 927–934. [Google Scholar] [CrossRef] [Green Version]
  26. Zhou, X.; Xiao, D.; Hou, Z.; Li, Q.; Wu, Y.; Wu, X. Influences of the structure parameters on sensitivity and Brownian noise of the disk resonator gyroscope. J. Microelectromech. Syst. 2017, 26, 519–527. [Google Scholar] [CrossRef]
  27. Shi, H.; Fan, S.; Li, W. Design and optimization of DETF resonator based on uncertainty analysis in a micro-accelerometer. Microsyst. Technol. 2018, 24, 2025–2034. [Google Scholar] [CrossRef]
  28. Xu, Y.; Zhao, L.; Jiang, Z.; Ding, J.; Peng, N.; Zhao, Y. A novel piezoresistive accelerometer with SPBs to improve the tradeoff between the sensitivity and the resonant frequency. Sensors 2016, 16, 210. [Google Scholar] [CrossRef] [Green Version]
  29. Du, X.; Wang, L.; Li, A.; Wang, L.; Sun, D. High accuracy resonant pressure sensor with balanced-mass DETF resonator and twinborn diaphragms. J. Microelectromech. Syst. 2016, 26, 235–245. [Google Scholar] [CrossRef]
  30. Zhao, C.; Wood, G.S.; Xie, J.; Chang, H.; Pu, S.H.; Kraft, M. A three degree-of-freedom weakly coupled resonator sensor with enhanced stiffness sensitivity. J. Microelectromech. Syst. 2015, 25, 38–51. [Google Scholar] [CrossRef]
  31. Yan, L.; Shangchun, F.; Zhanshe, G.; Jing, L.; Le, C. Study of dynamic characteristics of resonators for MEMS resonant vibratory gyroscopes. Microsyst. Technol. 2012, 18, 639–647. [Google Scholar] [CrossRef]
  32. Jing, L.; Shang-Chun, F.; Yan, L.; Zhan-She, G. Dynamic stability of parametrically-excited linear resonant beams under periodic axial force. Chin. Phys. B 2012, 21, 110401. [Google Scholar]
  33. Wang, H.; Li, Y. Dynamic analysis of the resonator for resonant accelerometer. Sens. Rev. 2018, 38, 478–485. [Google Scholar] [CrossRef]
Figure 1. Schematic of the resonant accelerometer.
Figure 1. Schematic of the resonant accelerometer.
Micromachines 12 00530 g001
Figure 2. Simplified resonator implementation.
Figure 2. Simplified resonator implementation.
Micromachines 12 00530 g002
Figure 3. Schematic of the displacement of the resonant beam.
Figure 3. Schematic of the displacement of the resonant beam.
Micromachines 12 00530 g003
Figure 4. Structural parameters of the resonant beam.
Figure 4. Structural parameters of the resonant beam.
Micromachines 12 00530 g004
Figure 5. Theoretical and numerical frequencies of the resonator at a fixed-width ratio: (a) silicon; (b) steel.
Figure 5. Theoretical and numerical frequencies of the resonator at a fixed-width ratio: (a) silicon; (b) steel.
Micromachines 12 00530 g005
Figure 6. Theoretical and numerical frequencies of the resonator at a fixed-length ratio: (a) silicon; (b) steel.
Figure 6. Theoretical and numerical frequencies of the resonator at a fixed-length ratio: (a) silicon; (b) steel.
Micromachines 12 00530 g006
Figure 7. Internal resonant frequency of the resonator at a fixed-width ratio: (a) silicon; (b) steel.
Figure 7. Internal resonant frequency of the resonator at a fixed-width ratio: (a) silicon; (b) steel.
Micromachines 12 00530 g007
Figure 8. Internal resonant frequency of the resonator at a fixed-length ratio: (a) silicon; (b) steel.
Figure 8. Internal resonant frequency of the resonator at a fixed-length ratio: (a) silicon; (b) steel.
Micromachines 12 00530 g008
Figure 9. The mode diagram of resonator structure with single-ended fixed support: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Figure 9. The mode diagram of resonator structure with single-ended fixed support: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Micromachines 12 00530 g009
Figure 10. The mode diagram of the resonator structure with double-ended fixed support: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Figure 10. The mode diagram of the resonator structure with double-ended fixed support: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Micromachines 12 00530 g010
Figure 11. Resonator structure diagram.
Figure 11. Resonator structure diagram.
Micromachines 12 00530 g011
Figure 12. Beam root 1: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Figure 12. Beam root 1: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Micromachines 12 00530 g012
Figure 13. Beam root 2: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Figure 13. Beam root 2: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Micromachines 12 00530 g013
Figure 14. Beam root 3: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Figure 14. Beam root 3: (a) first-order mode; (b) second-order mode; (c) third-order mode; (d) fourth-order mode.
Micromachines 12 00530 g014
Figure 15. Stiffness frequency of three resonators.
Figure 15. Stiffness frequency of three resonators.
Micromachines 12 00530 g015
Figure 16. Mode analysis of accelerometer: (a) interference mode; (b) working mode.
Figure 16. Mode analysis of accelerometer: (a) interference mode; (b) working mode.
Micromachines 12 00530 g016
Figure 17. Schematic of the accelerometer 1 g tumbling experiment.
Figure 17. Schematic of the accelerometer 1 g tumbling experiment.
Micromachines 12 00530 g017
Figure 18. Acceleration frequency of two resonators.
Figure 18. Acceleration frequency of two resonators.
Micromachines 12 00530 g018
Figure 19. Sensitivity of the resonant accelerometer.
Figure 19. Sensitivity of the resonant accelerometer.
Micromachines 12 00530 g019
Table 1. Physical parameters of the resonator.
Table 1. Physical parameters of the resonator.
Density/(Kg/m3)Young’s Modulus/(GPa)Poisson Ratio
Silicon23001700.278
Steel78002110.3
Table 2. Parameters of the resonant beam and beam root.
Table 2. Parameters of the resonant beam and beam root.
ParametersResonant BeamBeam Root 1Beam Root 2Beam Root 3
Length (mm)20222
Width (mm)0.30.511.5
Thickness (mm)
Stiffness (N/m)
5
3375
5
1.56 × 107
5
1.25 × 108
5
4.22 × 108
Table 3. Parameters of the resonant accelerometer.
Table 3. Parameters of the resonant accelerometer.
ParametersValueUnit
Density7.8 × 103kg/m3
Young’s modulus2 × 1011Pa
External frame (length × width × thickness)100 × 100 × 5mm3
Poisson ratio0.3/
Inner frame (length × width × thickness)80 × 80 × 5mm3
Proof mass (length × width × thickness)60 × 50 × 5mm3
Resonant beam (length × width × thickness)20 × 0.3 × 5mm3
Beam root (length × width × thickness)2 × 1.5 × 5mm3
Axial force1.2N
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Jin, B.; Zhao, M.; Yang, F. Design and Optimization of the Resonator in a Resonant Accelerometer Based on Mode and Frequency Analysis. Micromachines 2021, 12, 530. https://doi.org/10.3390/mi12050530

AMA Style

Li Y, Jin B, Zhao M, Yang F. Design and Optimization of the Resonator in a Resonant Accelerometer Based on Mode and Frequency Analysis. Micromachines. 2021; 12(5):530. https://doi.org/10.3390/mi12050530

Chicago/Turabian Style

Li, Yan, Biao Jin, Mengyu Zhao, and Fuling Yang. 2021. "Design and Optimization of the Resonator in a Resonant Accelerometer Based on Mode and Frequency Analysis" Micromachines 12, no. 5: 530. https://doi.org/10.3390/mi12050530

APA Style

Li, Y., Jin, B., Zhao, M., & Yang, F. (2021). Design and Optimization of the Resonator in a Resonant Accelerometer Based on Mode and Frequency Analysis. Micromachines, 12(5), 530. https://doi.org/10.3390/mi12050530

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