Next Article in Journal
Recent Developments in (K, Na)NbO3-Based Lead-Free Piezoceramics
Previous Article in Journal
Analysis and Verification of Heat Dissipation Structures Embedded in Substrates in Power Chips Based on Square Frustums Thermal through Silicon Vias
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Temperature Drift Error Estimation Model for Capacitive MEMS Gyros Using Thermal Stress Deformation Analysis

College of Intelligent Systems Science and Engineering, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Micromachines 2024, 15(3), 324; https://doi.org/10.3390/mi15030324
Submission received: 3 November 2023 / Revised: 26 January 2024 / Accepted: 23 February 2024 / Published: 26 February 2024

Abstract

:
Because the conventional Temperature Drift Error (TDE) estimation model for Capacitive MEMS Gyros (CMGs) has inadequate Temperature Correlated Quantities (TCQs) and inaccurate parameter identification to improve their bias stability, its novel model based on thermal stress deformation analysis is presented. Firstly, the TDE of the CMG is traced precisely by analyzing its structural deformation under thermal stress, and more key decisive TCQs are explored, including ambient temperature variation ∆T and its square ∆T2, as well its square root ∆T1/2; then, a novel TDE estimation model is established. Secondly, a Radial Basis Function Neural Network (RBFNN) is applied to identify its parameter accurately, which eliminates local optimums of the conventional model based on a Back-Propagation Neural Network (BPNN) to improve bias stability. By analyzing heat conduction between CMGs and the thermal chamber with heat flux analysis, proper temperature control intervals and reasonable temperature control periods are obtained to form a TDE precise test method to avoid time-consuming and expensive experiments. The novel model is implemented with an adequate TCQ and RBFNN, and the Mean Square Deviation (MSD) is introduced to evaluate its performance. Finally, the conventional model and novel model are compared with bias stability. Compared with the conventional model, the novel one improves CMG’s bias stability by 15% evenly. It estimates TDE more precisely to decouple Si-based materials’ temperature dependence effectively, and CMG’s environmental adaptability is enhanced to widen its application under complex conditions.

1. Introduction

Capacitive MEMS Gyros (CMGs) have so many merits, such as significant miniaturization, high integration, low power consumption, and cost, and they are widely utilized and applied in diverse fields, for example, inertial navigation, attitude detection, and heading guidance [1,2,3,4]. However, because CMGs are manufactured with Si-based materials of temperature dependence, temperature must change the physical properties of Si-based materials, which introduces Temperature Drift Error (TDE) to reduce CMG’s bias stability and limits their applications in some fields [5]. Intrinsically, optimizing the production process of Si-based materials to decouple their temperature dependence is the fundamental way to eliminate TDE [6]. Due to current material process limitations, it is so hard to decouple temperature dependence in the present technological context. In the background of the great repeatability of CMGs, temperature control and mathematical estimation are mainstream methodologies to eliminate TDE [7,8,9]. The first one uses a temperature control system to stabilize ambient temperature, and it has the merits of high temperature control accuracy. In contrast, its demerits, high power consumption, large volume, and complex control, limit it to be used on a large scale [10,11,12]. The second one is based on given models to estimate TDE precisely and calibrate CMG’s outputs in real time. Although it has the merits of low cost and easy implementation, its accuracy depends on precise TDE traceability and parameter identification of TDE estimation models [13,14].
On one hand, TDE traceability is important to enhance CMG’s bias stability [15]. Considering the consistent working principle of capacitive MEMS devices, Temperature Correlated Quantities (TCQs), exciting TDE is universal and invariant, and it ought to be figured out completely to offer references to describe TDE more accurately. Based on many experiments, Maj et al. show ambient temperature variations of 1 °C introduce a sensitivity variation of 1% of the scale factor of a capacitive MEMS accelerometer, which indicates that the direct cause of TDE is ambient temperature [16]. Kim et al. research the influence of ambient temperature on three-dimensional sizes of Si-based materials, and the experimental results point out that Si-based materials deform linearly with ambient temperature [17]. So, a composite model for estimating TDE is built with ambient temperature and TDE. After estimation and compensation, the bias stability of capacitive MEMS devices is improved by 10%. Wen et al. figured out that ambient temperature variation is key to describing structural deformation inside CMGs [18]. Ambient temperature and its variation are utilized as the TCQ of TDE, and the Radial Basis Function Neural Network (RBFNN) is used to estimate TDE precisely; then, CMG’s bias stability is improved by 10%. Bekkeng explores TCQs, ambient temperature, and their variation’s square, and a Kalman filter is applied to estimate TDE [19]. Then, the absolute rate error of MEMS inertial devices is reduced by more than a factor of 10. Chen et al. use linear microstructure thermal analysis to show microstructure deformation of the combs in MEMS capacitive accelerometers and obtain TCQs, ambient temperature variation, and their variation’s square [20]. TDE is estimated with TCQs and RBFNNs, and its bias stability is improved by 60%. But, microstructure deformation is analyzed from an uncomprehensive perspective, and the factors that may cause nonlinear deformation are not studied to confirm whether nonlinear deformation occurs. Shi et al. propose a multi-physical analysis to research TDE [21]. It reveals that the deviation of the thermal expansion coefficient between simulations and experiments is about 1~2 mg/°C, and thermal expansion properties of Si-based materials are nonlinear, which is non-ignorable to TDE traceability. Ambient temperature and its square and ambient temperature variation and its square are explored as TCQs, and a novel TDE estimation model is established with them and a Back Propagation Neural Network (BPNN) using particle swarm optimization plus genetic algorithm. The experimental results show that its bias stability is improved by 16.01% compared with the conventional model. Furthermore, a conclusion is formed that Si-based materials have complex structural deformation under the action of bending moment and elastic modulus, and it further shows that partial structural deformation inside CMGs is nonlinear [22,23,24,25]. So, the key to tracing TDE accurately lies in exploring complete TCQs based on the nonlinear structural deformation of Si-based materials, especially under thermal stress.
On the other hand, precise parameter identification which describes the relationship between TCQs and TDE is another key factor, and it requires the prerequisites that the TDE estimation model is accurate. An accurate TDE estimation model focuses on the relationship between its inputs and outputs, and the bias stability of CMGs is compensated to the target ones with its inputs [26]. Cheng et al. apply a particle swarm optimization algorithm in optimizing support vector machine models to improve bias stability. Small batch data processing methods are used to ensure its modeling in real time [27]. Although it improves capacitive MEMS accelerometers’ bias stability by up to 18.96%, support vector machine models with complicated structures are unsuitable for disposing a lot of experimental results and even reduce in real time. Pan et al. show a TDE estimation model using a wavelet neural network, which improves bias stability after compensation to 10% of that before [28]. Due to the diversity and uncertainty of signal transmission in wavelet neural networks, extensive training is needed to identify its precise parameters. Xu et al. propose a TDE estimation model based on a BPNN, and the nonlinear maximum of TDE is reduced from 3329 ppm to 603 ppm [29]. As a kind of classical artificial neural network, the BPNN has a simple structure, including an input layer, one or more hidden layers, and an output layer, which achieve more outstanding estimation performance to ensure its higher precision and better real time performance, as well as easier implementation in applications compared with other methods [30]. As is known, when BPNNs are trained, their local optimums are prone to occur and may introduce non-optimal estimates of TDE. Wang et al. use a genetic algorithm to assist in training BPNNs to avoid them [31]. The bias stability of the compensated capacitive MEMS accelerometer is less than 0.017% during −10~80 °C, and it is 173 times more precise than the TDE estimation model with polynomial fitting. But, the genetic algorithm has the demerits of probability disorder, reducing its calculation accuracy and real time performance and increasing the training amount of BPNNs and the targeted model’s identification difficulty. So, a proper parameter identification is able to describe the complex nonlinearity and estimate TDE precisely. Additionally, it should have more real time to compensate for CMG’s bias stability with its TCQs to ensure that it outputs angular velocity timely, even in the case that the carrier rotates at a large angular velocity. Compared with the BPNN, the RBFNN is a neural network based on the function approximation method with three neural layers. The function approximation method ensures that it describes the targeted complex nonlinearity as accurately as possible, and three neural layers reduce calculation time to guarantee TDE estimated in real time [32]. So, the RBFNN is a more proper choice compared to other neural networks. However, inaccurate experimental results result in its imprecise implementation to influence its estimation accuracy, so accurate experimental results are an important premise to ensure RBFNN’s performance.
In the paper, TDE traceability for CMGs is explored using thermal stress deformation analysis by simulating structural deformation under diverse conditions, and all-new TCQs are obtained as the key factors to TDE traceability, ambient temperature variation ∆T, and its square ∆T2, as well as its square root ∆T1/2. To improve TDE estimation performance, an RBFNN is used to modify the conventional model based on a BPNN and establish a novel model. Moreover, by analyzing the heat conduction between CMGs and thermal chambers principally, a precise TDE test is established, and two important characteristics are determined: proper temperature control interval and a reasonable temperature control period. The novel model estimates TDE more accurately to decouple the temperature dependence of Si-based materials, which further enhances the environmental adaptability of CMGs.

2. Methodology

2.1. Precise TDE Traceability

2.1.1. Conventional TDE Estimation Model

CMGs are miniaturized devices manufactured with Si-based materials and include mass, driving circuit, sensing circuit, and substrate. With microelectromechanical technology, they are assembled as microelectromechanical units. Figure 1 shows their basic principles.
Where kD and kS are the springs of the driving circuit and sensing circuit, respectively, and CD and CS are their corresponding capacitors. According to Figure 1, the sensing circuit and driving circuit have interspersed combs, and driving combs and sensing combs are movable. Under the action of driving combs, sensing combs move along the driving axis. Because of air gaps between driving combs, they can be abstracted as plate capacitors. One capacitance of sensing combs C0 is shown with the plate capacitor’s formula as follows:
C 0 = ε 4 π k S 0 d 0
where ε is the relative dielectric constant, k is the electrostatic force constant, S0 is the overlap area between the sensing combs, and d0 is its plate distance. In Figure 1, when the carriers rotate, sensing combs displace away from their original position under the action of the Coriolis force, and their displacements are related to the carriers’ angular velocity. The higher their angular velocity is, the greater air gaps between the sensing combs change, and the more significantly their capacitance varies. Because sensing combs have a common end, they are connected in parallel in terms of electrical form. Assuming that sensing combs have 2n combs, there are n air gaps, and capacitance variation ∆C measured by CMGs is shown as follows:
Δ C = i = 1 n ( C i C 0 ) = n [ ε 4 π k S 0 ( d 0 - Δ d ) ε 4 π k S 0 d 0 ]
where ∆d is plate distance variation and Ci is the capacitance between each air gap. So, the carriers’ rotation is precisely sensed with capacitance variation. Also, sensing and driving circuits’ stiffness determines CMG’s structural consistency. Considering that their stiffness is so related to Si-based materials’ stiffness, ambient temperature is considered the root cause [18]. Sensing and driving circuits deform as the ambient temperature varies inevitably and the capacitance errors of sensing combs are excited, as are the TDE ∆E in CMG’s output. The conventional TDE estimation model describes CMGs’ structural deformation using linear analytical methods and shows that ∆E is determined by ambient temperature variation ∆T and its square ∆T2. So, the conventional model is shown as follows:
Δ E c a p ( Δ T , Δ T 2 )

2.1.2. A Novel TDE Estimation Model with Thermal Stress Deformation Analysis

According to Figure 1, the mass, driving circuit, and sensing circuit are fabricated on the substrate, and they displace together as a whole with the substrate. Even if it deforms as ambient temperature, their relative distances between each other remain relatively stable, especially the air gaps. So, the substrate’s structural deformation has no influence on CMG’s performance, and it can be negligible. In contrast, the structural deformation of sensing combs changes air gaps between them and worsens the accuracy of capacitance variation. So, their structural deformation is the key factor to be considered comprehensively.
Thermal stress is an internal stress inside objects. When ambient temperature varies, the stress generated by the objects due to external and mutual constraints makes its structure unable to fully expand and contract freely. It makes CMG’s structure deform under the excitation of ambient temperature variation and changes the consistency of air gaps between the sensing combs, and the accuracy of capacitance measurements are reduced directly. Then, the influence of thermal stress on the capacitance measurement of sensing combs is analyzed using thermal stress deformation analysis in diverse conditions.
(a)
Ambient temperature T = T0 and angular velocity ω = ω0
When ambient temperature stabilizes at T0, CMG’s structure stays stable. If the carriers rotate at ω = ω0, sensing combs displace in the y-axis under the Coriolis force. Figure 2 shows the structure of sensing combs when T = T0 and ω = ω0.
Where a0 is the comb thickness of sensing combs, b0 is the length of overlap area between the sensing combs, c0 is their width, and S0 = b0c0. Because ambient temperature stays still at T0, CMG’s structure remains unchanging, and the capacitance of the sensing combs is constant. According to (2), ∆C measured by CMGs is shown as follows:
Δ C = i = 1 n ( C i C 0 ) = n ε 4 π k ( b 0 c 0 d 0 - Δ d b 0 c 0 d 0 )
(b)
Ambient temperature T = T1 and angular velocity ω = ω0
When T = T0, it is unmistakable for sensing combs to deform in three dimensions because of the temperature dependence of Si-based materials. According to Thermal Expansion Theory, structural deformation includes expanding or contracting, which changes the air gaps of sensing combs and their capacitance. So, the multi-physical field simulation software COMSOL Multiphysics 6.0 is used to show its structural deformation. Because CMG’s structure is at a micrometer scale, about 1~100 μm, specific physical quantities of CMG’s structure are set at a common value [33,34]. The thermal expansion coefficient of Si-based materials is 2.4 × 10−6/°C, its elastic modulus is 175 GPa, and its Poisson’s ratio is 0.28. Figure 3 shows the structure of sensing combs under thermal stress when T = T1 and ω = ω0 in simulations with COMSOL Multiphysics 6.0.
According to Figure 3, the combs have complex structural deformations under thermal stress, including linear and nonlinear. For non-connecting ends of the sensing combs, which have structural characteristics of long beams of Si-based materials, they deform freely in three dimensions following the linear thermal expansion formula. After deforming, its structure is still a long beam whose three-dimensional sizes scale up or down proportionally as the ambient temperature varies [35]. In contrast, thermal stress from the connecting ends acts on sensing combs, limiting their free structural deformation, and the bending moment from thermal stress brings local nonlinear deformation. According to (4), the capacitance variation of sensing combs is related to the overlap area of the sensing combs and their plate distance. According to Figure 3, some characteristics of structural deformation are obtained as follows:
  • Thermal stress from the connecting ends limits sensing combs to deform laterally, like its width and thickness, and it is useless for longitudinal structural deformation at all. So, its length deforms freely and can be calculated with a linear thermal expansion formula.
  • Thermal stress from the connecting ends and free structural deformation of the non-connecting ends make sensing combs’ length vary in a curved line. Given that its width is bigger than its thickness in normal conditions, thermal stress to the width is also bigger than that to its thickness, which means the curve radian of the thickness varies more quickly than the width. Considering that sensing combs are arranged in a positive and negative way, the overlap area is represented as an approixmate rectangle whose length and width should be calculated with the linear thermal expansion formula.
  • According to Figure 3, because the the curve radian of the thickness of the overlap area varies much more quickly than its width, the distance between the sensing combs is shown as a “wide-narrow-wide” pattern. It causes plate distance to vary nonlinearly. In a word, it should be calculated with a nonlinear expression.
Based on that, the length of the overlap area b1 with linear thermal expansion formula can be shown as follows:
b 1 = b 0 + Δ b + Δ b = b 0 + α b 0 Δ T + α b 0 Δ T = b 0 ( 1 + 2 α Δ T )
where α is the thermal expansion coefficient of Si-based materials, ∆b′ is the length deformation of a sensing comb, ∆b″ is that of the next one, and ∆T is ambient temperature variation. Then, the width of the overlap area with linear thermal expansion formula c1 is shown as follows:
c 1 = c 0 + 2 Δ c = c 0 + 2 α c 0 Δ T = c 0 ( 1 + 2 α Δ T )
where ∆c is width deformation. According to Figure 3, the structural deformation of sensing combs is near the connecting end, which influences the in-plate distance of the overlap area but not its area. In a word, the area of the overlap area can be seen as constant, but its plate distance varies at any point of sensing combs. Figure 4 shows the structural deformation of the long beam of Si-based materials under thermal stress.
Given that the sensing combs are axisymmetric, a three-dimensional positive half-axis coordinate system is established on a comb for simplified calculation in Figure 4. Additionally, the bending moment from thermal stress is used to describe nonlinear structural deformation. So, its bending moment M at any position of the sensing combs can be calculated as follows [36]:
M = 0 a 0 2 α E Δ T b 0 ( a 0 2 x ) d x = α E b 0 a 0 2 8 Δ T
where E is the elastic modulus of Si-based materials. According to Figure 4, the thickness of sensing combs a1 under its bending moment is shown as follows:
a 1 = a 0 + 2 Δ a = a 0 + 2 M E I x 2 2 = a 0 + 12 α Δ T a 0 x 2 x ( 0 , b 0 )
where ∆a is thickness deformation, I is the inertia moment of the combs, and I = b0(a0/2)3/12. So, the plate distance under its bending moment d1 is shown as follows:
d 1 = d 0 Δ a P Δ a N = d 0 12 α Δ T a 0 x 2 [ 12 α Δ T a 0 ( b 0 x ) 2 ] = d 0 12 α Δ T a 0 ( 2 x 2 2 b 0 x + b 0 2 )
where ∆aP is the thickness deformation of sensing combs arranged in a positive way and ∆aN is the thickness deformation arranged in a negative way. If the carriers rotate at ω = ω0, the capacitance in the sensing circuit ∆C′ is shown as follows:
Δ C = i = 1 n ( C i C 0 ) = n 0 b 0 ( ε 4 π k b 1 c 1 d 1 Δ d ε 4 π k b 1 c 1 d 1 ) d x   = n ε 4 π k b 0 c 0 ( 1 + 2 α Δ T ) ( 1 + 2 α Δ T ) 0 b 0 [ 1 d 0 12 α Δ T a 0 ( 2 x 2 2 b 0 x + b 0 2 ) Δ d 1 d 0 12 α Δ T a 0 ( 2 x 2 2 b 0 x + b 0 2 ) ] d x   = n ε b 0 c 0 4 π k ( 1 + 2 α Δ T ) 2 1 24 α Δ T a 0 0 b 0 [ 1 ( x b 0 2 ) 2 + d 0 Δ d 24 α Δ T a 0 b 0 2 4 1 ( x b 0 2 ) 2 + d 0 24 α Δ T a 0 b 0 2 4 ] d x   = n ε b 0 c 0 4 π k ( 1 + 2 α Δ T ) 2 1 4 6 α Δ T a 0 ( 1 d 0 Δ d 24 α Δ T a 0 b 0 2 4 ln | b 0 2 + d 0 Δ d 24 α Δ T a 0 b 0 2 4 b 0 2 + d 0 Δ d 24 α Δ T a 0 b 0 2 4 | 1 d 0 24 α Δ T a 0 b 0 2 4 ln | b 0 2 + d 0 24 α Δ T a 0 b 0 2 4 b 0 2 + d 0 24 α Δ T a 0 b 0 2 4 | )
In (2), the capacitance errors measured in sensing circuit ∆Ecap′ can be expressed as follows:
Δ E c a p = Δ C Δ C   = n ε b 0 c 0 4 π k ( 1 + 2 α Δ T ) 2 1 4 6 α Δ T a 0 ( 1 d 0 Δ d 24 α Δ T a 0 b 0 2 4 ln | b 0 2 + d 0 Δ d 24 α Δ T a 0 b 0 2 4 b 0 2 + d 0 Δ d 24 α Δ T a 0 b 0 2 4 | 1 d 0 24 α Δ T a 0 b 0 2 4 ln | b 0 2 + d 0 24 α Δ T a 0 b 0 2 4 b 0 2 + d 0 24 α Δ T a 0 b 0 2 4 | ) n ε 4 π k ( b 0 c 0 d 0 - Δ d - b 0 c 0 d 0 )   = f 1 ( Δ T , Δ T 2 ) f 2 ( Δ T ) f 3 ( Δ T ) f ( Δ T , Δ T 2 , Δ T )
According to (11), there is a deviation between CMG’s output and its reference when ambient temperature varies. (11) can be expressed as the product of three uncorrelated functions that cannot be simplified further, and it is related to ambient temperature variation ∆T and its square ∆T2, as well as its square root ∆T1/2. So, the conventional model is modified to a novel model, which is shown as follows:
Δ E c a p f ( Δ T , Δ T 2 , Δ T 1 / 2 )

2.2. Precise Parameter Identification for the Novel TDE Precise Estimation Model

Based on precise TDE traceability and a novel model for CMGs, identifying its parameters is another important factor. The more precisely the novel model is identified, the less ambient temperature influences TDE, and the better the bias stability of CMGs. So, testing TDE accurately is the key prerequisite, and a proper TDE test method is essential.

2.2.1. TDE Accurate Acquisition Methodology

(a)
Heat conduction solutions
To ensure the experimental results’ reliability, heat from ambient temperature should be conducted to CMGs efficiently. Heat conduction solutions are applied on CMGs, like thermal grease and thermal conductive rubber, which keep the ambient temperature the same as the CMG and reduce the deviation between experimental results and their theoretical values.
(b)
Precise temperature measurement system
Precise ambient temperature represents CMG’s environmental adaptability comprehensively, so it is necessary to apply a precise temperature measurement system. It needs be installed closely on the CMG to reduce the temperature gradient. Its measuring accuracy needs to be over 2 times more accurate than ambient temperature variation, and its measuring frequency needs to be higher than CMG’s output frequency for complete experimental results.
(c)
Proper temperature control interval
TDE is a statistical characteristic deviation represented based on mathematical models. Once CMGs are manufactured, their mechanical properties are fixed, and their TDE can be expressed with bias, scale factor, and random error. So, CMGs have similar environmental adaptability but unique TDEs, and from their datasheets, TDEs can be roughly estimated as follows:
Δ E = γ Δ T + β Δ T
where ∆E′ are roughly estimated TDEs, γ is the parameter “Zero-rate level change vs. temperature”, and β is the parameter “Sensitivity change vs. temperature”. Theoretically, ∆E′ is smaller than TDE because some non-statistical errors are ignored and ∆E′ ≤ ∆E. CMG’s sensitivity ∆ES determines the measured minimum of angular velocity, and higher sensitivity brings more accurate angular velocity. When ambient temperature deteriorates suddenly during CMG work, it is probable that TDE is much bigger than ∆ES, even covering up CMG’s sensitivity to measure inaccurate angular velocity. So, it is essential to vary ambient temperature at a reasonable temperature variation rate. To test TDE precisely, set ∆ES ≈ ∆E′, and the temperature control interval can be further expressed as follows:
Δ T Δ E S γ + β
(d)
Reasonable temperature control period
According to heat conduction theory, a long enough heat conduction period conducts the heat from position A to position B completely, which ensures that ambient temperature variation at position A is the same as position B. An insufficient heat conduction period causes an incomplete heat conduction process, and the TDE of the CMG is tested imprecisely. Usually, a thermal chamber is used to test TDE. It is designed as a cube of front-door-open and -closed insulation. Also, it integrates a temperature control system with six temperature control units on each plane, and they can be further simplified according to the targeted control efficiency. There is a rate table inside to offer the referenced angular velocity. To reflect the ambient temperature of the CMG better, temperature sensors are set on it. Figure 5 shows a schematic diagram of the test platform.
In Figure 5, the heat is transferred through the thermal chamber’s inner wall to control ambient temperature inside the targets, and the heat from temperature control units uniformly transfers to the CMG. Then, the center of the rate table is taken as a reference position, and it is also the central position of the thermal chamber. Given that the sizes of the thermal chamber are L1 × L2 × L3, the farther the position is away from the inner wall, the longer the period of heat transfers. Under the action of temperature control units, CMGs obtain the heat, and their environmental adaptability changes. The position where the CMG is located is the final area where the ambient temperature stays stable. According to energy conservation law and Newton’s law of cooling, the heat flux density equation in three dimensions is shown as follows:
q = k h ( u x i + u y j + u z k )
where u is the ambient temperature of any position in the thermal chamber at moment t, q is its heat flux density, kh is the heat conductivity coefficient, u / x is the spatial varying rate of ambient temperature on the x-axis, u / y is the spatial varying rate of the ambient temperature on the y-axis, u / z is the spatial varying rate of the ambient temperature on the z-axis, i is the unit vector on the x-axis, j is the unit vector on the y-axis, and k is the unit vector on the z-axis. The heat on the x-axis heats CMGs, and according to (15), a heat transfer equation can be established as follows:
Q x = [ ( q x ) A ( q x ) A + Δ A ] Δ y Δ z Δ t = [ ( k h u x ) A ( k h u x ) A + Δ A ] Δ y Δ z Δ t = k h 2 u x 2 Δ x Δ y Δ z Δ t
where Qx is the heat at position A on the x-axis, (qx)A is its heat flux density, (qx)A+ΔA is its heat flux density of its close position A + ΔA near position A on the x-axis, Δy is the width of the plane where the heat conducts, Δz is the height of the plane where the heat conducts, and Δt is the period when the heat conducts. Considering that the CMG is heated on the x-axis, y-axis, and z-axis, a heat transfer equation can be established on the y-axis and z-axis as follows:
Q y = k h 2 u y 2 Δ x Δ y Δ z Δ t Q z = k h 2 u z 2 Δ x Δ y Δ z Δ t
where Qy is the heat at position A on the y-axis and Qz is the heat of position A on the z-axis. Based on (16) and (17), the heat at position A in three dimensions can be deduced as follows:
Q = Q x + Q y + Q z = k h Δ x Δ y Δ z Δ t ( 2 u x 2 i + 2 u y 2 j + 2 u z 2 k ) + F ( x , y , z , t ) Δ x Δ y Δ z Δ t v
where F(x,y,z,t) is the heat flux density of the potential heat sources related to its position and time, v is the 3D unit vector matrix, and v = [ i j k ] T . When the thermal chamber works in a period ts, according to the specific heat capacity formula, we can obtain the following equation:
Q = [ k h Δ x Δ y Δ z ( 2 u x 2 i + 2 u y 2 j + 2 u z 2 k ) + F ( x , y , z , t ) Δ x Δ y Δ z v ] t s = C m Δ T v
where C is the specific heat capacity of air in a thermal chamber in a closed condition, m is its mass, ∆T is temperature variation, and ∆T = |Tb-T0|. T0 is the initial temperature and Tb is the final temperature. Then, (19) can be transformed as follows:
[ k h Δ x Δ y Δ z ( 2 u x 2 i + 2 u y 2 j + 2 u z 2 k ) + F ( x , y , z , t ) Δ x Δ y Δ z v ] t s = C ρ Δ x Δ y Δ z Δ T v
where ρ is the air density inside the thermal chamber. Considering that different positions in the heat conduction pathway have different heat flux densities in the heat conduction process, the further the heat is conducted, the smaller the heat flux density becomes, and the lower heat transfer efficiency stays. The CMG is a sensor, and its heat is relatively insignificant to the thermal chamber, so F(x,y,z,t) can be ignored in magnitude. Then, (20) can be deduced as follows:
( 2 u x 2 i + 2 u y 2 j + 2 u z 2 k ) = C ρ Δ T k h t s v
To calculate total heat along the heat conduction pathway, (21) is integrated on the x-axis, y-axis, and z-axis separately, and then we obtain the following equation:
0 L 1 0 L 1 2 u x 2 d x d x + 0 L 2 0 L 2 2 u y 2 d y d y + 0 L 3 0 L 3 2 u z 2 d z d z = Δ T C ρ k h t s ( 0 L 2 0 L 3 d y d z + 0 L 1 0 L 3 d x d z + 0 L 1 0 L 2 d x d y )
After integration, (22) is transformed and shown as follows:
T x + T y + T z = Δ T C ρ k h t s ( L 2 L 3 + L 1 L 3 + L 1 L 2 )
where Tx is the control target of the temperature control unit on the x-axis, Ty is on the y-axis, and Tz is on the z-axis. Then, ts can be deduced as follows:
t s = C ρ k h ( L 2 L 3 + L 1 L 3 + L 1 L 2 ) ( T x + T y + T z ) Δ T
Usually, a thermal chamber is designed as a cube, and its size is L × L × L. Moreover, the control target on the x-axis, y-axis, and z-axis is set the same as Tb, and (24) can be simplified further as follows:
t s = C ρ L 2 k h T b ( T b T 0 )
Then, the period for heat conduction from the inner wall to the center of the thermal chamber can be calculated with (25). To guarantee that the thermal chamber is heated uniformly, the temperature control period is set as tstp. Taking L3GD20H as an example, it is chosen to test TDE randomly. From its datasheet, ∆ES = 8.75 mdps/digit, γ = ±0.04 dps/°C, and its operating temperature range is −40~85 °C. After dimensional transformation, β is obtained as follows:
β = 2 % / ° C × 8.75 mdps × [ ( 85   ° C ) ( 40   ° C ) ] 2 = 10.9375 mdps / ° C
According to (14), temperature control interval ∆T is shown as follows:
Δ T 8.75 mdps / digit 0.04 dps / ° C   × ( 245 dps / 2000 dps ) + 10.9375 mdps / ° C   0.55   ° C
To simplify the test steps, ∆T = 0.5 °C. A thermal chamber SET-Z-021 is utilized to test L3GD20H, and its parameters are C = 1.005 kJ/(kg×K), kh = 0.0267 W/m°C, L = 0.6 m, and ρ = 1.293 kg/m3. Given that CMG’s operating temperature range is −40~85 °C, (25) is obtained as follows:
t s 25.766   s
So, temperature control units take 25.766 s as the temperature control period to vary the temperature control interval at 0.5 °C. To simplify test steps and reserve an allowance for stable heat conduction, tp = 30 s. L3GD20H is tested in temperature experiments, and its temperature is obtained with a precise temperature measurement system of measuring accuracy of ±0.03 °C and measuring frequency of 10 Hz. Thus, a temperature experiment is designed as follows:
  • CMG L3GD20H is installed on the rate table, its measuring direction is parallel to the rate table, and the referenced true value is the angular velocity of the rate table.
  • Temperature sensors of the precise temperature measurement system are attached close to L3GD20H, the wireless data transmission module transmits the experimental results, and the PC is prepared to receive its temperature TG and its output DG.
  • Cool the thermal chamber to a minimum operating temperature of −40 °C and keep TG and DG recording for 0.5 h when the ambient temperature stays stable.
  • Heat the thermal chamber to a maximum operating temperature of 85 °C at a rate of 60 °C/h, 0.5 °C per 30 s. When TG goes up to 85 °C, stop the test when it stays stable for 0.5 h.
  • Redo steps (2) to (4) five times and record them as the experimental results.
Figure 6 shows a flowchart of the temperature experiment and experimental results.

2.2.2. Implementation of Novel TDE Precise Estimation Model Based on an RBFNN

Under the premise of building a novel TDE precise compensation model and obtaining TDE and TCQs, identifying its precise parameters is another key to estimating the TDE of the CMG accurately. As shown in Figure 6, the ambient temperature goes up from −40 °C to 85 °C, and at the very beginning, the ambient temperature stays stable for 0.5 h, and it is set as the reference ambient temperature. Moreover, the reference output of the CMG is 0 dps/s because ω s = 0 . ∆T, ∆T2, ∆T1/2, and TDE deduced from Figure 6 are shown in Figure 7.
As shown in Figure 7, when ambient temperature T varies, ∆T and ∆T2 vary in a similar trend. Because of the small numerical amplitude of ambient temperature variation at the beginning, ∆T1/2 has the same trend with a small range of numerical fluctuations. Under the excitation of ∆T and ∆T2, as well as ∆T1/2, TDE has an approximate varying trend. Therefore, it shows that there is some relevance among ∆T, ∆T2, ∆T1/2, and TDE. Figure 8 shows the corresponding relationship among them.
As shown in Figure 8, there is a complex nonlinear relationship among ∆T, ∆T2, ∆T1/2, and TDE, and it is so difficult to estimate precisely in real time. So, it is essential to use a nonlinear model of multiple inputs and a single output as well as high accuracy to fit the nonlinearity. The RBFNN is a feedforward neural network based on the function approximation method with a single hidden layer. It consists of neurons and neural layers. Neurons are applied as basic computing units, and neural layers are utilized as a computing framework. Its neural layers include an input layer, a hidden layer, and an output layer, and neurons are distributed in diverse layers. The neurons in the hidden layer and output layer have kernel functions, and inputs in the RBFNN are transmitted between the input layer and the output layer by neurons and kernel functions. By means of neurons and neural layers, as well as kernel functions, the output in the RBFNN perfectly represents or approximately approaches the targets. Usually, it uses the Gaussian function as the kernel function in the hidden layer and the purelin function as the activation function in the output layer, and they are shown as follows:
y = e x c j 2 / 2 σ j 2         y = p u r e l i n ( x ) = a x + b
where cj is the center of the Gaussian function and δj is its width, a and b are constant, x is the input of the functions, and y are their outputs. Moreover, the RBFNN has two merits as follows:
  • Owing to RBFNN’s mathematical principles, its calculation results are optimal in global scope to avoid local minimums, even in flat areas where the error gradient approximates to zero.
  • From the Kolmogorov theorem, a three-layer forward network is able to approach any continuous function with any target accuracy [18]. The RBFNN has the structure of an input layer, a hidden layer, and an output layer, and it can represent the targeted nonlinearity with any accuracy.
Figure 9 shows the structure of the RBFNN for the complex nonlinear relationship among ∆T, ∆T2, ∆T1/2, and TDE.
Where Ii(i = 1⋯3) is the ith neuron of the input layer, Hi(i = 1⋯K) is the ith neuron of the hidden layer, and O1 is the neuron of the output layer. Based on that, the RBFNN is a good choice to represent the nonlinear relationship among ∆T, ∆T2, ∆T1/2, and TDE accurately. So, (12) can also be deduced as follows:
Δ E c a p = R B F N N ( Δ T , Δ T 2 , Δ T 1 / 2 )
Then, the parameters of the novel model should be identified as follows:
  • Two temperature experiments are conducted, and the experimental results are recorded. One of them is a training sample set, and the other one is a verification sample set.
  • Based on sample data of the training sample set, TDE is calculated by subtracting the reference value of CMGs from the sample data of their actual outputs one by one. ∆T is calculated by subtracting the reference temperature of CMGs from the sample data of their actual temperature one by one. Then, ∆T2 is obtained by multiplying itself, and ∆T1/2 is obtained by calculating the square root of ∆T.
  • The RBFNN uses ∆T, ∆T2, and ∆T1/2 as its inputs and TDE as its output. It is trained with mathematical tools until the differences between its outputs and targeted TDE meet the requirements.
  • The compensated results of CMGs are calculated from subtraction between the actual outputs of CMGs and their corresponding estimated outputs of RBFNNs.
The sample data of the training sample set are determined by MEMS Gyros’ output frequency under the target operating conditions. The higher MEMS Gyros’ output frequency is, the more quickly MEMS Gyros outputs, and the higher the sample data of the training sample set. Then, (30) is trained with the experimental results, and its parameters are identified precisely. Figure 10 shows the primary data and those compensated by the conventional model and novel model in five experiments.
To check its estimation performance, the Mean Square Deviation (MSD) is introduced as follows:
M S D = M S E ( x x )
where x is the evaluated sample, x′ is its reference, and the MSE is the mean square error algorithm. The MSD is an index that indicates the dispersion degree of the evaluated sample and its reference. The less the MSD is, the smaller the dispersion degree stays. The bias stability of CMGs is the evaluated sample. So, the MSD can be implemented with programming software based on mathematical principles. The MSD of the bias stability of its primary data is MSD1, the bias stability of those compensated by the conventional model is MSD2, and the bias stability of those compensated by the novel model is MSD3, which are shown in Table 1. An MSD improvement of the conventional model Q1 and the novel model Q2 are shown as follows:
Q 1 = M S D 2 M S D 1 Q 2 = M S D 3 M S D 1
In Figure 10, CMGs compensated by the conventional and novel models perform stably with small fluctuations around the reference, but CMGs compensated by the novel model run more stably and have smaller fluctuations, which shows that the novel model can estimate TDE more accurately. In Table 1, although the conventional and novel models reduce the MSD of the CMG effectively, the novel model reduces the MSD of the CMG less, and the MSD after compensation is enhanced by four orders of magnitude. Moreover, CMG’s bias stability is enhanced better by 10.5% evenly than the conventional model. So, the novel model can estimate TDE more accurately to improve the bias stability of the CMG significantly.

3. Experiments and Analysis

To verify the estimation universality and the robustness of the novel model further, CMG I3G4250D is taken as a test object. It is chosen randomly and tested in five verification experiments. The conventional model is based on the least square method and ∆T (Conventional 1), the conventional one is based on the BPNN and ∆T and ∆T2 (Conventional 2), and the novel one is based on the RBFNN and ∆T, ∆T2, and ∆T1/2 (Novel). They are established on the x-axis, y-axis, and z-axis, and their bias stabilities are compared. To ensure experimental results’ reliability, the referenced angular velocity in three axes is set as ωrefx = 25 dps, ωrefy = 5 dps, and ωrefz = 53 dps randomly. Figure 11 compares the experimental results in the experiments.
According to (31), the MSD of the bias stability of its primary data is MSD1, the bias stability of those compensated by Conventional 1 is MSD4, the bias stability of those compensated by Conventional 2 is MSD5, and the bias stability of those compensated by Novel is MSD6. They are shown as follows:
M S D 1 = M S E ( k k i ) M S D 4 = M S E [ k ( Δ T ) k i ] M S D 5 = M S E [ k ( Δ T , Δ T 2 ) k i ] M S D 6 = M S E [ k ( Δ T , Δ T 2 , Δ T 1 / 2 ) k i ]
where k is the primary data of I3G4250D in three axes; k(∆T) are those compensated by Conventional 1 in three axes; k(∆T, ∆T2) are those compensated by Conventional 2 in three axes; k(∆T, ∆T2, ∆T1/2) are those compensated by Novel in three axes; and ki′(i = x, y, z) is the referenced angular velocity in three axes. In (31), the MSD improvement of Conventional 1 is Q3, the MSD improvement of Conventional 2 is Q4, and the MSD improvement of Novel is Q5. MSDs in five experiments are shown in Table 2, Table 3, Table 4, Table 5 and Table 6.
Q 3 = M S D 4 M S D 1 Q 4 = M S D 5 M S D 1 Q 5 = M S D 6 M S D 1
According to Figure 11, CMGs compensated by Conventional 1, Conventional 2, and Novel perform stably with small fluctuations around the references, and the bias stability of CMGs compensated by the novel model is reduced the most significantly. It shows that the novel model estimates TDE more accurately and effectively, which ensures CMGs perform more stably and reliably. According to Table 2, Table 3, Table 4, Table 5 and Table 6, the MSDs of CMGs compensated by Novel are smaller than Conventional 1 and Conventional 2, and its MSD after compensation is improved about two orders of magnitude at minimum. Although their difference is small, there is a performance improvement of the test data in a statistical sense. So, the primary data’s performance must be improved at a certain moment; it may be much better or a little better, and its actual performance must be improved more than the MSD. Furthermore, Novel and Conventional 2 can enhance CMG’s bias stability better by 15% than Conventional 1, and Novel can enhance CMG’s bias stability more perfectly by 5% than Conventional 2. Hence, the novel model is able to estimate TDE more accurately and enhance the bias stability of CMGs markedly. By comparison, the conventional model has higher TDE estimation accuracy and more perfect universality and robustness, which guarantees CMGs to perform stably and accurately for a long time in diverse complicated conditions.

4. Conclusions

In this paper, TDE traceability for CMGs is explored adequately with thermal stress deformation analysis by simulating structural deformation in diverse conditions, and ∆T,T2, and ∆T1/2 are deduced as more key decisive TCQs to establish the novel TDE precise estimation model. To increase TDE estimation accuracy, an RBFNN is applied to implement the novel model and avoid the local optimums of BPNNs. By analyzing heat conduction between CMGs and the thermal chamber, a precise TDE test is established with heat flux analysis. Two key parameters, the proper temperature control interval and the reasonable temperature control period, are obtained to avoid time-consuming and expensive experiments. The MSD is used to evaluate TDE estimation accuracy, and the bias stability of CMGs compensated by the novel model is improved by about 15% compared with the conventional ones. The novel model can estimate TDE more precisely to decouple the temperature dependence of Si-based materials, which enhances the environmental adaptability of CMGs to widen their application in different complicated conditions. This study provides a universal and robust means of estimating the TDE of CMGs in diverse applications, optimizing the TDE precise estimation model to perform precisely and in real time in fields, such as industrial production, smart agriculture, infrastructure construction, and so on, offering a proper TDE precise test method without the need for time-consuming and expensive experiments.

Author Contributions

Conceptualization, B.Q. and J.C.; methodology, B.Q. and Z.W.; validation, B.Q. and Z.W.; formal analysis, B.Q. and J.C.; investigation, B.Q., J.C. and Z.W.; resources, J.C.; data curation, C.J. (Chao Jiang) and C.J. (Chun Jia); writing—original draft preparation, B.Q.; writing—review and editing, J.C.; visualization, Z.W. and C.J. (Chao Jiang); supervision, J.C.; project administration, B.Q. and J.C.; funding acquisition, B.Q. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was jointly funded by the National Natural Science Foundation of China (Nos. 62003108, 62073093, 61773132, 61803115, 62003109), the National Key Research and Development Program (No. 2021YFB3901300), the 145 High-tech Ship Innovation Project sponsored by the Chinese Ministry of Industry and Information Technology, the Heilongjiang Province Research Science Fund for Excellent Young Scholars (No. YQ2020F009), and the Fundamental Research Funds for Central Universities (Nos. 3072019CF0401, 3072020CFJ0402, 3072020CFT0403).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Verma, P.; Agrawal, P.; Gopal, R.; Arya, S.K. Parametric Sensitivity Analysis of a 2-DOF Drive and 1-DOF Sense Modes MEMS Gyro-Accelerometer Structure. Adv. Sci. Lett. 2015, 20, 1495–1498. [Google Scholar] [CrossRef]
  2. Shen, X.; Yuan, D.; Chang, R.; Jin, W. A Nonlinear Observer for Attitude Estimation of Vehicle-Mounted Satcom-on-the-Move. IEEE Sens. J. 2019, 19, 8057–8066. [Google Scholar] [CrossRef]
  3. Cai, Q.; Zhao, F.J.; Kang, Q.; Luo, Z.Q.; Hu, D.; Liu, J.W.; Cao, H.L. A Novel Parallel Processing Model for Noise Reduction and Temperature Compensation of MEMS Gyroscope. Micromachines 2021, 12, 1285. [Google Scholar] [CrossRef]
  4. Yu, Y.Y.; Luo, H.; Chen, B.Y.; Tao, J.; Feng, Z.H.; Zhang, H.; Guo, W.L.; Zhang, D.H. MEMS Gyroscopes Based on Acoustic Sagnac Effect. Micromachines 2017, 8, 2. [Google Scholar] [CrossRef]
  5. Yang, Y.; Liu, Y.; Liu, Y.H.; Zhao, X.D. Temperature Compensation of MEMS Gyroscope Based on Support Vector Machine Optimized by GA. In Proceedings of the 2019 IEEE Symposium Series on Computational Intelligence, Xiamen, China, 6–9 December 2019. [Google Scholar]
  6. Omar, A.; Elshennawy, A.; AbdelAzim, M.; Ismail, A.H. Analyzing the Impact of Phase Errors in Quadrature Cancellation Techniques for MEMS Capacitive Gyroscopes. In Proceedings of the 2018 17th IEEE SENSORS Conference, New Delhi, India, 28–31 October 2018. [Google Scholar]
  7. Guan, R.; He, C.H.; Liu, D.; Zhao, Q.; Yang, Z.; Yan, G. A temperature control system used for improving resonant frequency drift of MEMS gyroscopes. In Proceedings of the 10th IEEE International Conference on Nano/Micro Engineered and Molecular Systems, Xi’an, China, 7–11 April 2015. [Google Scholar]
  8. Maspero, F.; Delachanal, S.; Berthelot, A.; Joet, L.; Langfelder, G.; Hentz, S. Quarter-mm2 High Dynamic Range Silicon Capacitive Accelerometer With a 3D Process. IEEE Sens. J. 2020, 20, 689–699. [Google Scholar] [CrossRef]
  9. Wang, X.; Bleiker, S.J.; Edinger, P.; Errando-Herranz, C.; Roxhed, N.; Stemme, G.; Gylfason, K.B.; Niklaus, F. Wafer-Level Vacuum Sealing by Transfer Bonding of Silicon Caps for Small Footprint and Ultra-Thin MEMS Packages. J. Microelectromechanical Syst. 2019, 28, 460–471. [Google Scholar] [CrossRef]
  10. Castellanos-Sahagun, E.; Alvarez, J. Temperature-Temperature cascade control of Binary Batch Distillation Columns. In Proceedings of the 2013 European Control Conference, ETH Zurich, Zurich, Switzerland, 17–19 July 2013. [Google Scholar]
  11. Wang, G.L.; Yan, W.W.; Chen, S.H.; Zhang, X.; Shao, H.H. Multivariable constrained predictive control of main steam temperature in ultra-supercritical coal-fired power unit. J. Energy Inst. 2015, 88, 181–187. [Google Scholar] [CrossRef]
  12. Xu, S.; Hashimoto, S.; Jiang, W.; Jiang, Y.Q.; Izaki, K.; Kihara, T.; Ikeda, R. Slow Mode-Based Control Method for Multi-Point Temperature Control System. Processes 2019, 7, 533. [Google Scholar] [CrossRef]
  13. Wu, X.Y.; Wang, X.N.; He, G.G. A Fuzzy Self-tuning Temperature PID Control Algorithms for 3D Bio-printing Temperature Control System. In Proceedings of the 2020 32nd Chinese Control And Decision Conference, Hefei, China, 22–24 August 2020. [Google Scholar]
  14. Dong, H.; Zhang, W.; Li, M.X.; Zeng, Z.B.; Cao, D.; Li, X.P. High-precision air temperature control considering both hardware elements and controller design. Case Stud. Therm. Eng. 2022, 37, 102290. [Google Scholar] [CrossRef]
  15. Krysko, V.A.; Awrejcewicz, J.; Yakovleva, T.V.; Kirichenko, A.V.; Szymanowska, O.; Krysko, V.A. Mathematical modeling of MEMS elements subjected to external forces, temperature and noise, taking account of coupling of temperature and deformation fields as well as a nonhomogenous material structure. Commun. Nonlinear Sci. Numer. Simul. 2019, 72, 39–58. [Google Scholar] [CrossRef]
  16. Maj, C.; Szermer, M.; Zajac, P.; Amrozik, P. Analytical modelling of MEMS Z-axis comb-drive accelerometer. In Proceedings of the 2019 20th International Conference on Thermal, Mechanical and Multi-Physics Simulation and Experiments in Microelectronics and Microsystems (EuroSimE), Hannover, Germany, 24–27 March 2019. [Google Scholar]
  17. Kim, B.; Hopcroft, M.A.; Candler, R.N.; Jha, C.M.; Agarwal, M.; Melamud, R.; Chandorkar, S.A.; Yama, G.; Kenny, T.W. Temperature Dependence of Quality Factor in MEMS Resonators. J. Microelectromechanical Syst. 2008, 17, 755–766. [Google Scholar] [CrossRef]
  18. Qi, B.; Wen, F.Z.; Liu, F.M.; Cheng, J.H. Modification of MTEA-Based Temperature Drift Error Compensation Model for MEMS-Gyros. Sensors 2020, 20, 2906. [Google Scholar] [CrossRef] [PubMed]
  19. Bekkeng, J.K. Calibration of a Novel MEMS Inertial Reference Unit. IEEE Trans. Instrum. Meas. 2009, 58, 1967–1974. [Google Scholar] [CrossRef]
  20. Qi, B.; Chen, J.Y.; Tian, S.S. A Novel Temperature Drift Error Precise Estimation Model of MEMS Capacitive Accelerometers Using MTAM. In Proceedings of the 2022 International Conference on Guidance, Navigation and Control, ICGNC 2022: Advances in Guidance, Navigation and Control, Harbin, China, 5–7 August 2022. [Google Scholar] [CrossRef]
  21. Qi, B.; Shi, S.S.; Zhao, L.; Cheng, J.H. A Novel Temperature Drift Error Precise Estimation Model for MEMS Accelerometers Using Microstructure Thermal Analysis. Micromachines 2022, 13, 835. [Google Scholar] [CrossRef] [PubMed]
  22. Dai, G.; Li, M.; He, X.; Du, L.; Shao, B.; Su, W. Thermal drift analysis using a multiphysics model of bulk silicon MEMS capacitive accelerometer. Sens. Actuators A-Phys. 2011, 172, 369–378. [Google Scholar] [CrossRef]
  23. Kinoshita, T.; Kawakami, T.; Hori, T.; Matsumoto, K.; Kohara, S.; Orii, Y.; Yamada, F.; Kada, M. Thermal Stresses of Through Silicon Vias and Si Chips in Three Dimensional System in Package. J. Electron. Packag. 2013, 134, 020903. [Google Scholar] [CrossRef]
  24. Hong, C.C. Thermal bending analysis of shear-deformable laminated anisotropic plates by the GDQ method. Mech. Res. Commun. 2009, 36, 804–810. [Google Scholar] [CrossRef]
  25. Ab Rahim, R.; Bais, B.; Majlis, B.Y.; Fareed, S. Design Optimization of MEMS Dual-Leg Shaped Piezoresistive Microcantilever. In Proceedings of the 2013 IEEE Regional Symposium on Micro and Nanoelectronics, Daerah Langkawi, Malaysia, 25–27 September 2013. [Google Scholar]
  26. Basarab, M.; Giani, A.; Combette, P. Thermal Accelerometer Simulation by the R-Functions Method. Appl. Sci. 2020, 10, 8373. [Google Scholar] [CrossRef]
  27. Cheng, J.; Fang, J.; Wu, W.; Li, J. Temperature drift modeling and compensation of RLG based on PSO tuning SVM. Measurement 2014, 55, 246–254. [Google Scholar] [CrossRef]
  28. Pan, Y.J.; Li, L.L.; Ren, C.H.; Luo, H.L. Study on the compensation for a quartz accelerometer based on a wavelet neural network. Meas. Sci. Technol. 2010, 21, 105202. [Google Scholar] [CrossRef]
  29. Xu, D.; Yang, Z.; Zhao, H.; Zhou, X. A temperature compensation method for MEMS accelerometer based on LM-BP neural network. In Proceedings of the 2016 IEEE SENSORS, Orlando, FL, USA, 30 October–3 November 2016. [Google Scholar]
  30. Wang, C.C.; Lin, C.J. Dynamic Analysis and Machine Learning Prediction of a Nonuniform Slot Air Bearing System. ASME J. Comput. Nonlinear Dyn. 2023, 18, 011007. [Google Scholar] [CrossRef]
  31. 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. Actuators A-Phys. 2020, 316, 112393. [Google Scholar] [CrossRef]
  32. Wang, C.C.; Kuo, P.H.; Chen, G.Y. Machine Learning Prediction of Turning Precision Using Optimized XGBoost Model. Appl. Sci. 2022, 12, 7739. [Google Scholar] [CrossRef]
  33. Kou, Z.W.; Kong, Z.; Jing, G.L.; Cui, X.M.; Ge, L.; Lu, P.Z. Design of Capacitive MEMS Ring Solid-state Vibrating Gyroscope Interface Circuit. In Proceedings of the 2022 IEEE 6th Advanced Information Technology, Electronic and Automation Control Conference (IAEAC), Beijing, China, 3–5 October 2022. [Google Scholar]
  34. Oouchi, A. Plastic molded package technology for MEMS sensor evolution of MEMS sensor package. In Proceedings of the 2014 International Conference on Electronics Packaging (ICEP), Toyama, Japan, 23–25 April 2014. [Google Scholar]
  35. Datskos, P.G.; Rajic, S.; Datskou, I. Photoinduced and Thermal Stress in Silicon Microcantilevers. Appl. Phys. Lett. 1998, 73, 2319–2321. [Google Scholar] [CrossRef]
  36. Lu, K.H.; Ryu, S.K.; Zhao, Q.; Hummler, K.; Im, J.; Huang, R.; Ho, P.S. Temperature-dependent Thermal Stress Determination for Through-Silicon-Vias (TSVs) by Combining Bending Beam Technique with Finite Element Analysis. In Proceedings of the 2011 IEEE 61st Electronic Components and Technology Conference, Lake Buena Vista, FL, USA, 31 May–3 June 2011. [Google Scholar]
Figure 1. Basic principle of the CMG. (a) Physical design; (b) working principle.
Figure 1. Basic principle of the CMG. (a) Physical design; (b) working principle.
Micromachines 15 00324 g001
Figure 2. Partial view of the structure of sensing combs when T = T0 and ω = ω0. (a) Top view of sensing combs; (b) side view of sensing combs.
Figure 2. Partial view of the structure of sensing combs when T = T0 and ω = ω0. (a) Top view of sensing combs; (b) side view of sensing combs.
Micromachines 15 00324 g002
Figure 3. Partial view of the structure of sensing combs under thermal stress with COMSOL Multiphysics 6.0 when T = T1 and ω = ω0. (a) Top view in simulations; (b) side view in simulations.
Figure 3. Partial view of the structure of sensing combs under thermal stress with COMSOL Multiphysics 6.0 when T = T1 and ω = ω0. (a) Top view in simulations; (b) side view in simulations.
Micromachines 15 00324 g003
Figure 4. Structural deformation of long beams with and without thermal stress. (a) Top view of the long beam; (b) side view of the long beam.
Figure 4. Structural deformation of long beams with and without thermal stress. (a) Top view of the long beam; (b) side view of the long beam.
Micromachines 15 00324 g004
Figure 5. A schematic diagram of the test platform. (a) Front view of the CMG installation; (b) vertical view of the CMG installation; (c) field test of the CMG.
Figure 5. A schematic diagram of the test platform. (a) Front view of the CMG installation; (b) vertical view of the CMG installation; (c) field test of the CMG.
Micromachines 15 00324 g005aMicromachines 15 00324 g005b
Figure 6. Flowchart of the temperature experiment and experimental results. (a) Experimental flow; (b) experimental results.
Figure 6. Flowchart of the temperature experiment and experimental results. (a) Experimental flow; (b) experimental results.
Micromachines 15 00324 g006
Figure 7. Ambient temperature variation and its square, as well as its square root and TDE.
Figure 7. Ambient temperature variation and its square, as well as its square root and TDE.
Micromachines 15 00324 g007
Figure 8. The relationship among ∆T, ∆T2, ∆T1/2, and TDE.
Figure 8. The relationship among ∆T, ∆T2, ∆T1/2, and TDE.
Micromachines 15 00324 g008
Figure 9. The structure of the RBFNN for the complex nonlinear relationship.
Figure 9. The structure of the RBFNN for the complex nonlinear relationship.
Micromachines 15 00324 g009
Figure 10. Comparison of the primary data of CMGs and their compensated outputs in five experiments. (a) First experiment; (b) second experiment; (c) third experiment; (d) fourth experiment; (e) fifth experiment.
Figure 10. Comparison of the primary data of CMGs and their compensated outputs in five experiments. (a) First experiment; (b) second experiment; (c) third experiment; (d) fourth experiment; (e) fifth experiment.
Micromachines 15 00324 g010
Figure 11. Comparison of the experimental results in five verification experiments. (a) First experiment; (b) second experiment; (c) third experiment; (d) fourth experiment; (e) fifth experiment.
Figure 11. Comparison of the experimental results in five verification experiments. (a) First experiment; (b) second experiment; (c) third experiment; (d) fourth experiment; (e) fifth experiment.
Micromachines 15 00324 g011aMicromachines 15 00324 g011b
Table 1. Performance comparison between the conventional and novel models in experiments.
Table 1. Performance comparison between the conventional and novel models in experiments.
MSD1MSD2MSD3Q1Q2|Q2Q1|/Q1 (%)
First experiment9.55523.5555 × 10−23.2172 × 10−23.7210 × 10−33.3670 × 10−39.51%
Second experiment9.55793.5957 × 10−23.1756 × 10−23.7620 × 10−33.3224 × 10−311.68%
Third experiment9.55673.4492 × 10−23.1616 × 10−23.6092 × 10−33.3083 × 10−38.34%
Fourth experiment9.55903.5734 × 10−23.1339 × 10−23.7383 × 10−33.2784 × 10−312.30%
Fifth experiment9.55073.5442 × 10−23.1660 × 10−23.7109 × 10−33.3150 × 10−310.67%
Table 2. MSDs of the experimental results in the first experiment.
Table 2. MSDs of the experimental results in the first experiment.
MSD1MSD4MSD5MSD6Q3Q4Q5|Q4Q3|/Q3 (%)|Q5Q3|/Q3 (%)|Q5Q4|/Q4 (%)
x-axis4.77940.02160.02140.01934.519 × 10−34.478 × 10−34.038 × 10−30.93%10.65%9.81%
y-axis8.57680.02800.02330.02143.265 × 10−32.717 × 10−32.495 × 10−316.79%23.57%8.15%
z-axis0.72580.02700.02200.02103.720 × 10−23.031 × 10−22.893 × 10−218.52%22.22%4.55%
Table 3. MSDs of the experimental results in the second experiment.
Table 3. MSDs of the experimental results in the second experiment.
MSD1MSD4MSD5MSD6Q3Q4Q5|Q4Q3|/Q3 (%)|Q5Q3|/Q3 (%)|Q5Q4|/Q4 (%)
x-axis3.85450.01940.01920.01845.033 × 10−34.981 × 10−34.774 × 10−31.03%5.15%4.17%
y-axis6.82650.02560.02190.02083.750 × 10−33.208 × 10−33.047 × 10−314.53%18.75%5.02%
z-axis0.65820.02710.02110.02024.117 × 10−23.206 × 10−23.069 × 10−222.14%25.46%4.27%
Table 4. MSDs of the experimental results in the third experiment.
Table 4. MSDs of the experimental results in the third experiment.
MSD1MSD4MSD5MSD6Q3Q4Q5|Q4Q3|/Q3 (%)|Q5Q3|/Q3 (%)|Q5Q4|/Q4 (%)
x-axis4.77790.02100.02020.01904.395 × 10−34.227 × 10−33.976 × 10−33.81%9.52%5.94%
y-axis8.59620.02890.02400.02213.362 × 10−32.792 × 10−32.571 × 10−316.96%23.53%7.92%
z-axis0.72860.02800.02270.02143.843 × 10−23.116 × 10−22.937 × 10−224.91%26.35%1.92%
Table 5. MSDs of the experimental results in the fourth experiment.
Table 5. MSDs of the experimental results in the fourth experiment.
MSD1MSD4MSD5MSD6Q3Q4Q5|Q4Q3|/Q3 (%)|Q5Q3|/Q3 (%)|Q5Q4|/Q4 (%)
x-axis3.85750.02000.01990.01895.185 × 10−35.159 × 10−34.900 × 10−30.51%5.50%5.03%
y-axis6.82080.02560.02210.02113.753 × 10−33.240 × 10−33.094 × 10−313.67%17.58%4.52%
z-axis0.65380.02730.02070.02004.176 × 10−23.166 × 10−23.059 × 10−224.18%26.74%3.38%
Table 6. MSDs of the experimental results in the fifth experiment.
Table 6. MSDs of the experimental results in the fifth experiment.
MSD1MSD4MSD5MSD6Q3Q4Q5|Q4Q3|/Q3 (%)|Q5Q3|/Q3 (%)|Q5Q4|/Q4 (%)
x-axis3.85380.01980.01930.01895.138 × 10−35.008 × 10−34.904 × 10−32.53%4.55%2.07%
y-axis6.83390.02560.02170.02123.746 × 10−33.175 × 10−33.102 × 10−315.23%17.19%2.30%
z-axis0.65880.02740.02090.02034.159 × 10−23.172 × 10−23.081 × 10−223.72%25.91%2.87%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qi, B.; Cheng, J.; Wang, Z.; Jiang, C.; Jia, C. A Novel Temperature Drift Error Estimation Model for Capacitive MEMS Gyros Using Thermal Stress Deformation Analysis. Micromachines 2024, 15, 324. https://doi.org/10.3390/mi15030324

AMA Style

Qi B, Cheng J, Wang Z, Jiang C, Jia C. A Novel Temperature Drift Error Estimation Model for Capacitive MEMS Gyros Using Thermal Stress Deformation Analysis. Micromachines. 2024; 15(3):324. https://doi.org/10.3390/mi15030324

Chicago/Turabian Style

Qi, Bing, Jianhua Cheng, Zili Wang, Chao Jiang, and Chun Jia. 2024. "A Novel Temperature Drift Error Estimation Model for Capacitive MEMS Gyros Using Thermal Stress Deformation Analysis" Micromachines 15, no. 3: 324. https://doi.org/10.3390/mi15030324

APA Style

Qi, B., Cheng, J., Wang, Z., Jiang, C., & Jia, C. (2024). A Novel Temperature Drift Error Estimation Model for Capacitive MEMS Gyros Using Thermal Stress Deformation Analysis. Micromachines, 15(3), 324. https://doi.org/10.3390/mi15030324

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