Next Article in Journal
Distributed Architecture to Integrate Sensor Information: Object Recognition for Smart Cities
Previous Article in Journal
Structural Health Monitoring System Based on FBG Sensing Technique for Chinese Ancient Timber Buildings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coupled Dynamic Analysis and Decoupling Optimization Method of the Laser Gyro Inertial Measurement Unit

1
Wuhan National Laboratory for Optoelectronics, Huazhong Institute of Electro-Optics, Wuhan 430223, China
2
School of Electromechanical and Architectural Engineering, Jianghan University, Wuhan 430010, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(1), 111; https://doi.org/10.3390/s20010111
Submission received: 18 August 2019 / Revised: 13 December 2019 / Accepted: 17 December 2019 / Published: 23 December 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
The mechanical dithered ring laser gyro (RLG) effectively overcomes the lock-in effect and ensures the sensitive accuracy of the low angular rate for the gyro. However, in the inertial measurement unit (IMU) system, the dither excitation of three RLGs causes the coupled vibration of the IMU structure, which could seriously limit the measuring accuracy of RLGs. In this paper, the vibration frequency response characteristic of laser gyro IMU is taken as the focus point, and the method of multi-rigid body dynamics is used to establish the dynamic model of IMU suitable for vibration frequency response analysis. On the basis of the model, the multi-degree-of-freedom coupling vibration of IMU with the gyro dither excitation is clearly described. A new IMU dynamic decoupling optimization method is proposed to minimize the coupled vibration frequency response, and compared with the previous optimal design method. The prototype experimental test results show that the coupled vibration of IMU is restrained more effectively by the proposed new method than by the previous optimal design method. Finally, on the basis of this new method, the measuring accuracy of the RLGs in the IMU system is improved, which is quite useful for practical engineering application.

1. Introduction

Gyroscopes are the key angular rate sensors. In the current applications of inertial navigation, two types of gyroscopes dominate, that is, the optoelectronic types including ring laser gyros (RLGs) [1,2] and fiber optic gyros (FOGs) [3,4], based on the Sagnac effect, and solid vibratory types such as hemispherical resonator gyros (HRGs) [5,6], based on the Coriolis force. Among these products, the RLG is the main inertial device owing to the high performance and reliabilities to vibration and temperature, especially for the harsh environments. The FOG senses the angular rate by phase shift between two counter-propagating waves in a fiber coil, and its performance is up to a similar grade to that of the RLG, but it may be more easily affected by external dynamic, thermal, or stress perturbations [7]. The HRG detects the precession angle of the stress wave in a quartz hemispherical resonator so as to sense the angular rate [8]. Compared with the optoelectronic gyros, the HRG has superiority for its high reliability, low power, and good environmental adaptability. On the other hand, the chip-scale optoelectronic gyroscopes will play important roles in the gyroscope market owing to the low costs and miniaturization [9]. The resonant micro-optic gyro (RMOG) is a promising sensor in the field of integrated optical gyros, and more innovative studies will focus on RMOG’s key technology, that is, the silica-on-silicon ring resonator with an ultra-high Q-factor and the monolithic integration based on gyro on a chip [10,11].
Even though more and more types of high-performance gyroscopes are being studied and manufactured, the RLG still plays an extremely important role in the field of high precision inertial navigation owing to its mature manufacturing process. In the family of RLGs, the mechanical dithered RLG is most widely applied, which overcomes the gyro’s lock-in effect by the dither bias system and ensures the high measurement accuracy of the gyro for a low angular rate [12].
However, the mechanical dither of RLG becomes the internal vibration source in the inertial navigation system (INS), especially for the inertial measurement unit (IMU), which is the core component of the INS. Even if the RLGs are orthogonally mounted on the support structure in the IMU, the spatial multi-degree-of-freedom (DOF) coupled dynamic behavior with the excitation of each gyro could still be caused by the asymmetry of the IMU structure, which could seriously affect the accuracy of the gyros in IMU and greatly restrict the performance of INS [13].
In order to eliminate the RLG’s lock-in effect so as to ensure the dither dynamics, many studies have been carried out in the areas of RLG’s dither control methods and the lock-in error compensation technology. In terms of dither control, Wang [14] deeply studied the influence of the dither frequency bias on the reduction of dynamic lock-in width and the accuracy of the gyro. Fan [15] accurately extracted the dither frequency bias through the sensitive angular rate of the gyro, and used it to improve the stability of dither control. In terms of lock-in error compensation, Song [16] proposed a new algorithm to estimate the Sagnac phase errors due to lock-ins, and to compensate the random walk error of the RLG. Subsequently, Fan [17] proposed a novel lock-in error correction method to pick up the lost information and remove the random walk error effectively. The above studies have offered important technical support for the design of the dither bias system, and meanwhile laid the foundations for the structure design of the laser gyro IMU.
The studies of the structural dynamics on the laser gyro IMU can be divided into two categories, specifically, the studies on the calibration of the systematic errors of RLGs in the INS, and the studies on the optimal design of the IMU structures. In terms of the former, the studies are almost focused on the calibration and compensation of the attitude errors of the RLGs in IMU, that is, Zheng [18] employed a multi-position calibration method to estimate the g-dependent bias of the RLG in the dual-axis INS; Kim and Wang [19,20] proposed the calibration technique of gyro sensitive-axis deviation angle under static and linear vibration conditions, respectively; and, recently, Tu [21] accurately analyzed the time-domain response of the gyro deviation angle in a biaxial INS through modal decoupling of the IMU.
In terms of the optimal design of the laser gyro IMU structure, a dynamic-vibration-absorber-type solution was implemented on the inertial sensor assembly by Jamil [22] to improve the response of the IMU structure to the RLG’s dither. Subsequently, many researchers focused on the IMU dither decoupling, and proposed the relevant optimal design principles [23,24,25] such as the design requirements of the overlap between the IMU centroid and the center of the multiple-damper layout, and that between the gyro sensitive axis and the principle inertia axis of IMU.
Although the design principles [23,24,25] are proposed, they are sometimes difficult to be satisfied in practical engineering owing to various constraints such as the size of IMU and the installation mode of printed circuit board (PCB) in actual IMU or INS design. Secondly, these design principles are not available for some new structures of IMU, that is, the skewed redundant INS [26], as the orthogonal layouts of the RLGs in the new types of IMU are always impossible. Furthermore, whatever the studies on the calibration or on the structure optimal design of the laser gyro IMU, the methods of modelling and optimization applied are always based on the traditional Newton–Euler method, and mainly focus on the time domain response of IMU vibration, which is not available for frequency decoupling.
Therefore, frequency-domain-analysis based modelling and optimization methods [27] may become the new key techniques for the vibration decoupling and optimization of the laser gyro IMU. Frequency-domain analytical methods focus on the analysis and optimization of the modal and frequency response for the mechatronic systems, which have been widely used in multiple precision manufacturing fields, that is, the dynamic analysis and optimization of the high-performance machine tools [28], the design of the active vibration isolators for the precision instruments [29], and the NVH (noise, vibration, and harshness) analysis of the vehicle [30], among others.
In this paper, the vibration frequency response characteristics of the laser gyro IMU are taken as the focus point, and the multi-rigid body dynamic model of IMU for the frequency response analysis is established to clearly describe the multi-DOF coupled vibration behavior in IMU under gyro dither excitation. Then, the IMU dynamic decoupling optimization method is proposed to minimize the coupled vibration response of IMU, which overcomes the limitation of the previous IMU optimal design principle constrained by the actual design conditions. The results of the prototype experiment show that the dynamic decoupling optimization method of IMU proposed in this paper could effectively suppress the coupled vibration in IMU and improve the measurement accuracy of the RLGs in the IMU system.
The remainder of this paper is organized as follows. The description of the structural components of laser gyro IMU is presented in Section 2. The dynamic model of laser gyro IMU is constructed and its analysis results are given in Section 3. The optimization design of IMU dynamic decoupling is given in Section 4. The experimental verification of IMU dynamic decoupling is presented in Section 5, and Section 6 provides the conclusion.

2. Structure Components of Laser Gyro IMU

IMU is the core component of the INS, which obtains the navigation information such as speed, position, and attitude of the carrier in real-time through processing and solving the inertial signals provided by the three-axis gyro and accelerometer [31]. A typical laser gyro IMU consists of a sensor supporting structure; an outer housing; vibration isolators (also called dampers); a group of orthogonally distributed measuring devices, such as three RLGs; and accelerometers. The RLGs and the accelerometers are mounted on the sensor supporting structure, and multiple vibration isolators connect the outer housing and the sensor supporting structure. Considering the bearing capacity of a single rubber damper, the IMU always uses an eight-point damping method, that is, the outer housing and the sensor supporting structure are connected by eight vibration isolators [32]. The structure diagram of the laser gyro IMU assembly is shown in Figure 1.

3. Dynamic Modelling and Discussion of Laser Gyro IMU

3.1. Multi-Rigid Body Dynamic Modeling of Laser Gyro IMU

Generally, the mechanical dither frequency of RLG is much lower than the first-order modal frequency of the sensor supporting structure in IMU. Therefore, the sensor supporting structure and the resonance cavity of RLG can be simplified into rigid bodies, and the dithering wheel is equivalent to a torsion spring. In the work of [33], dynamic structures are simplified as multiple rigid bodies, springs, and hinges, and a topological diagram is used to represent the constraint relations of the dynamic model. On that basis, in this paper, the multi-rigid body dynamic model of the IMU structure is developed based on the method of topology analysis of constraint relations. The three elements, such as spring-damper units, multiple rigid bodies, and ideal constrained pairs, are combined with the topology form of constraint relations of IMU to describe the RLG’s dither and multi-point vibration isolation in IMU.
The modeling process is divided into four steps.
Step 1: Definitions of the basic elements of the multi-rigid body system of IMU
The vibration displacements of the centroid of single rigid body relative to the global coordinate system are defined as the generalized coordinates of the dynamic equation: q i = [ x i y i z i α i β i γ i ] T ( i is the serial number of bodies in the system). The generalized coordinates ( α β γ ) set as small angular displacements, which can be described by cardan angles.
The schematic diagram of the IMU dynamic model is shown in Figure 2, and the global coordinate system is defined as { g x y z } , the dither coordinate systems of the three RLGs are { r i x y z }   ( i = 1 ~ 3 ), and the local coordinate systems of eight vibration isolators are { v j x y z }   ( i = 1 ~ 8 ). In this paper, the dynamic equation of IMU is described in the global coordinate system, which is fixed on the ground. The dither coordinate systems of RLGs and the coordinate systems of vibration isolators are fixed on the corresponding conjoined bases (i.e., three gyro cavities or eight vibration isolators), and they will change with the movements of the conjoined bases.
As shown in Figure 2, the whole IMU dynamic model includes four motion units, such as the sensor supporting structure and three gyro cavities. Without considering the constraint relationship, the model contains 6 ×   4 = 24 DOFs. Further, the rotation pair between the gyro cavity and the sensor support structure is defined to limit the other five DOFs between the IMU and each RLG, except the direction of dither. Whereby, the DOFs of model are finally reduced to nine, so the dynamic model is composed of a nine-order differential equation.
Step 2: Dynamic modeling of system without considering constraints
The serial numbers of the sensor supporting structure and the three gyro resonators are sequentially from 1 to 4. The centroid coordinate system of the sensor supporting structure coincides with the global coordinate system. The direction cosine matrix from the global coordinate system { Σ g } to the three gyro local coordinate systems { Σ r i } is defined as A i g r (i = 2~4), then the system mass matrix in the global coordinate system { Σ g } can be written as follows:
M = [ M 1 g M 2 g M 3 g M 4 g ] ,
where M i g = [ ( A i g r ) T 0 0 ( A i g r ) T ] M i r [ A i g r 0 0 A i g r ] , i = 2, 3, 4.
The stiffness matrices of eight vibration isolators in each of the local coordinate systems { Σ v j } can be written as follows:
K v j = d i a g ( k p k r k r 0 0 0 ) , j = 1 ~ 8 .
Equation (2) considers the linear stiffness ( k p k r k r ) of the vibration isolator along each of the three coordinate axes. The structure of the vibration isolator in IMU is given in Figure 3, and the stiffness parameters in the three axes can be checked according to the product manual provided.
Correspondingly, the stiffness matrices of three gyro resonators in each local coordinate system { Σ r i } can be written as follows:
K r i = d i a g ( k x i k y i k z i k r x i k r y i k d ) , i = 2 ~ 4 ,
where the dither stiffness of RLG is k d .
On the basis of the Lagrange method, without constraints, the dynamic equation of IMU can be derived as follows:
M q ¨ + K q = 0 ,
where M is represented in Equation (1), and K is as follows:
K = [ K 11 K 12 K 13 K 14 K 21 K 22 0 0 K 31 0 K 33 0 K 41 0 0 K 44 ] ,
where
K 11 = j = 1 8 ( T j g v ) T ( R j g v ) T K v i R j g v T j g v + i = 2 4 ( T i g r ) T ( R i g r ) T K r i R i g r T i g r ,
K 1 i = ( T i g r ) T ( R i g r ) T K r i R i g r I ,   i = 2 ~ 4 ,
K i i = I ( R i g r ) T K r i R i g r I ,   i = 2 ~ 4 ,
K i 1 = I ( R i g r ) T K r i R i g r T i g r ,   i = 2 ~ 4 .
In Equations (6)–(9), T j g v   (j = 1~8) and T i g r   (i = 2~4) are the vibration displacement transfer matrix from the IMU centroid to the center of each vibration isolator and to the dither center of each RLG, respectively. R j g v   (j = 1~8) and R i g r   (i = 2~4) are the pose transformation matrix from the global coordinate system { Σ g } to the local coordinate system { Σ v } of each vibration isolator and to the local coordinate system { Σ r } of each RLG’s resonator, respectively.
Step 3: Topology analysis of the constraint relation
An additional revolute pair is applied between each RLG’s resonant cavity and the sensor supporting structure. Eight vibration isolators connecting the sensor supporting structure and the ground are equivalent to spring-dampers. The constraint topological form of IMU is shown in Figure 4, in which the serial number 0 is the ground, 1 is the sensor supporting structure of IMU, and 2~4 are RLG’s resonant cavities.
The joint matrix H and the variable separation matrix D of each revolute pair are defined, respectively, as follows:
H = d i a g ( 1 1 1 1 1 0 ) ,
D = [ 0 0 0 0 0 1 ] .
Considering the revolute constraints between the RLGs and IMU, the established constraint matrix of IMU is as follows:
B = [ I 1 0 0 0 P 12 Q 12 0 0 P 13 0 Q 13 0 P 14 0 0 Q 14 ] [ I 1 D D D ] = [ I 1 0 0 0 P 12 Q 12 D 0 0 P 13 0 Q 13 D 0 P 14 0 0 Q 14 D ] ,
where
P 1 i = I ( R i g r ) T H R i g r T i g r ,   i = 2 ~ 4 ,
Q 1 i = I ( R i g r ) T ( I H ) ,   i = 2 ~ 4 .
Step 4: Dynamic modeling of the multi-rigid body system of IMU
Substituting the constraint matrix (12)–(14) into the dynamic Equation (4) of the non-constraint system, and considering the hysteresis damping characteristics of the rubber damper [34] (i.e., the elastic and damping effects of the vibration isolator are F = ( k + i ω c ) x = k ( 1 + i β ) x , where β is the hysteresis damping coefficient), meanwhile considering the quality factor Q of RLG’s dither system, so the dynamic model of IMU can be obtained as follows:
M ^ q ¨ + K ^ ( 1 + i β ) q = 0 ,
K ^ = [ K ^ 11 K ^ 12 K ^ 21 K ^ 22 ] ,
K ^ 11 = i = 1 8 [ k r 0 0 0 k r z i k r y i 0 k p 0 k p z i 0 k p x i 0 0 k r k r y i k r x i 0 0 k p z i k r y i k p z i 2 + k r y i 2 + k d 1 8 k r x i y i k p x i z i k r z i 0 k r x i k r x i y i k r z i 2 + k r x i 2 + k d 2 8 k r z i y i k r y i k p x i 0 k p x i z i k r z i y i k p x i 2 + k r y i 2 + k d 3 8 ] ,
K ^ 12 = [ 0 0 0 0 0 0 0 0 0 k d 1 0 0 0 k d 2 0 0 0 k d 3 ] ,
K ^ 21 = [ 0 0 0 k d 1 0 0 0 0 0 0 k d 2 0 0 0 0 0 0 k d 3 ] ,
K ^ 22 = [ k d 1 0 0 0 k d 2 0 0 0 k d 3 ] ,
M ^ = [ M ^ 11 0 0 M ^ 22 ] ,
M ^ 11 = [ 3 m + m b 3 m + m b m ( z r 2 + z r 3 + z r 4 ) 3 m + m b m ( y r 2 + y r 3 + y r 4 ) m ( z r 2 + z r 3 + z r 4 ) m ( y r 2 + y r 3 + y r 4 ) m ( y r 2 2 + y r 3 2 + y r 4 2 + z r 2 2 + z r 3 2 + z r 4 2 ) + 2 I x + I x x m ( z r 2 + z r 3 + z r 4 ) m ( x r 2 + x r 3 + x r 4 ) I x y m ( x r 2 y r 2 + x r 3 y r 3 + x r 4 y r 4 ) m ( y r 2 + y r 3 + y r 4 ) m ( x r 2 + x r 3 + x r 4 ) I x z m ( x r 2 z r 2 + x r 3 z r 3 + x r 4 z r 4 ) m ( z r 2 + z r 3 + z r 4 ) m ( y r 2 + y r 3 + y r 4 ) m ( x r 2 + x r 3 + x r 4 ) m ( x r 2 + x r 3 + x r 4 ) I x y m ( x r 2 y r 2 + x r 3 y r 3 + x r 4 y r 4 ) I x z m ( x r 2 z r 2 + x r 3 z r 3 + x r 4 z r 4 ) m ( x r 2 2 + x r 3 2 + x r 4 2 + z r 2 2 + z r 3 2 + z r 4 2 ) + I x + I y + I y y I y z m ( y r 2 z r 2 + y r 3 z r 3 + y r 4 z r 4 ) I y z m ( y r 2 z r 2 + y r 3 z r 3 + y r 4 z r 4 ) m ( x r 2 2 + x r 3 2 + x r 4 2 + y r 2 2 + y r 3 2 + y r 4 2 ) + 2 I y + I z z ] ,
M ^ 22 = [ I Z I Z I Z ] .
In Equations (16)–(23), ( x i y i z i )   (i = 1~8) are the positions of eight vibration isolators relative to the centroid of IMU, and ( x r i y r i z r i )   (i = 2~4) are the positions of three RLG’s cavities relative to the centroid of IMU. It should be emphasized that this paper focuses on the frequency response analysis of the coupled vibration behavior of three RLGs in IMU. Therefore, the external generalized force is not considered in the above dynamic model.

3.2. Discussion of Coupled Vibration Behavior of IMU

Through analyzing the stiffness matrix K ^ in Equations (15)–(23), we can find that the stiffness distribution of IMU is related to the positions ( x i   y i   z i ) of the vibration isolators relative to the IMU centroid (i = 1~8), besides the stiffness ( k p k r k r ) of the vibration isolators in three axes and the torsional stiffness ( k d 1 k d 2 k d 3 ) of three RLGs in each dither direction.
Similarly, the mass distribution of the system is related to the mass properties of the sensor supporting structure (including the mass parameter m b and inertia parameters [ I x x , I y y , I z z , I x y , I y z , I z x ] ), the inertia I z of each RLG along the respective dither directions, and the positions ( x r i   y r i   z r i ) (i = 2~4) of three RLGs relative to the centroid of the IMU.
The discussion of the coupled vibration behavior of the IMU under gyro dither is equivalent to the diagonalizing study of the mass matrix M ^ and stiffness matrix K ^ in the dynamic Equations (15)–(23). Therefore, the conditions for vibration decoupling of IMU can be derived as follows:
(1)
The IMU centroid coincides with the equivalent elastic center of the vibration isolators, that is, the eight vibration isolators, are symmetrically distributed with respect to the IMU. At this time, the stiffness matrix K ^ in the model is diagonalized. The specific condition is as follows:
i = 1 8 x i = i = 1 8 y i = i = 1 8 z i = i = 1 8 x i y i = i = 1 8 x i z i = i = 1 8 y i z i = 0 .
(2)
The principle inertia axes in the IMU centroid coordinate system are respectively along the dither axes of each RLG, that is, the non-diagonal terms in the mass matrix are zero. At this time, the inertia matrix M ^ in the model is diagonalized. The specific condition is as follows:
I x y m ( x r 2 y r 2 + x r 3 y r 3 + x r 4 y r 4 ) = I x z m ( x r 2 z r 2 + x r 3 z r 3 + x r 4 z r 4 ) = I y z m ( y r 2 z r 2 + y r 3 z r 3 + y r 4 z r 4 ) = 0 .
The similar decoupling studies can be found in the literatures [23,24,25,26], such as Equations (16)–(18) in the work of [23], Section 3.1 in the work of [25], and Section 4 in the work of [26]. Among the above previous studies, the work of [26] is the relatively recent study on the optimal design of the laser gyro IMU, where a new skewed redundant IMU by employing four RLGs is designed, and the dithered coupling of RLGs in the redundant configuration of IMU is considered. In the work of [26], the methods of increasing the weight of the IMU and allocating the RLGs’ frequencies reasonably are applied to reduce the coupled vibration, but no analytical model representing the coupled dynamics of IMU and relevant decoupling optimization method are involved.
Unlike the work of [26], other papers [23,24,25] propose similar decoupling conditions (i.e., the inertia product in the mass matrix is 0, namely I x y = I x z = I y z = 0 ) as Equations (24) and (25) in this paper.
By comparing the IMU dynamic decoupling conditions in Equations (24) and (25) and the previous decoupling conditions in the works of [23,24,25], some conclusions could be derived as follows:
(1)
The proposed decoupling condition such as Equation (25) is more general, compared with those in the works of [23,24,25]. The reason is that, in the actual IMU structure, the distance from the centroid of the IMU to each RLG is always not 0 (i.e.,   x r i 0 ,   y r i 0 ,   z r i 0 ), in that case, the dynamic decoupling is not achieved just by satisfying the equation I x y = I x z = I y z = 0 .
(2)
On the other hand, the condition of I x y = I x z = I y z = 0 means a fully symmetrical structure of IMU, which is also not easy to be satisfied in the actual design of INS owing to other design requirements (such as mass, size, and layouts of PCBs or sensors).
In summary, no matter which decoupling conditions are considered, such as Equations (24)–(25) in this paper or the similar ones in the works of [23,24,25], they are all based on the diagonalizing of the system stiffness and the mass matrices. However, in order to obtain the minimum-coupled vibration response of the IMU, as well as improve the engineering availability, it is still necessary to analyze the vibration response in the frequency domain, so as to propose a new optimal design method for the IMU structure, which will be focused in the next section in this paper.

4. Optimization Design of IMU Dynamic Decoupling

On the basis of the established IMU dynamic model in Section 3, the IMU coupling frequency response amplitude can be minimized by optimizing the IMU dynamic parameters. The frequency response function matrix under vibration excitation is obtained as H ( ω ) = [ ( I + j G ) K ω 2 M ] 1 [35]. Extracting the lth row and pth column of H ( ω ) , the following can be obtained:
H l p ( ω ) = x l ( ω ) f p ( ω ) = r = 1 N ψ l r ψ p r M r [ ( 1 + j g ) K r M r ω 2 ] .
Equation (26) represents that, if the point p is excited, the frequency response function of point l is obtained. If point p and point l are in different directions, Equation (26) represents the coupling frequency response of the system. At the same time, the formula also describes the mapping relationship between the structural parameters (positional relationship of vibration isolator relative to IMU centroid, the parameters of stiffness or damping, and so on) and dynamic response of the system.
According to Equation (26), with the aid of an intelligent algorithm, such as genetic algorithm, the coupling frequency response in IMU is regarded as the optimized objective, and the positions of vibration isolators relative to the IMU centroid ( ( x i , y i , z i ) in Equation (17)) are the optimized parameters to achieve the IMU dynamic decoupling.

4.1. Initial Parameters

The IMU structure before optimization and the serial number of RLGs and vibration isolators are shown in Figure 5.
In Figure 5, the sensor supporting structure of IMU is designed beforehand according to the decoupling principles of Equations (24) and (25), in order to make sure that the off-diagonal terms in the IMU inertia properties are smaller than the 1/50 of diagonal terms. The vibration isolators are chosen as AM003-15 rubber dampers manufactured by Corp. Lord. The other initial physical properties of the IMU structure are as follows:
(1)
Weight of sensor supporting structure: 9.1 kg;
(2)
Dither frequencies of RLGs: no. ①–③ are 490 Hz, 570 Hz, and 670 Hz, respectively;
(3)
Quality factor of the gyro: 150;
(4)
Stiffness of vibration isolator: axial 58 N/mm, radial 65 N/mm;
(5)
Damping coefficient of vibration isolator: 0.2;
(6)
Positions of the vibration isolators relative to the IMU centroid: as shown in Table 1.

4.2. Design Parameters

According to Equations (15)–(23), the structure parameters to determine the IMU dynamic response include the following: (1) weight of sensor supporting structure; (2) inertia distribution of the sensor supporting structure around its centroid; (3) the axial and radial stiffness of the vibration isolators; and (4) the positions and poses of the vibration isolators relative to the IMU centroid coordinate system.
Owing to the limitation of design factors, such as system layout, size, and installation mode, in this paper, the positions of eight vibration isolators relative to the IMU centroid are regarded as the design parameters. For the further sake of engineering achievability, only the x and z coordinates of each of vibration isolators in IMU centroid coordinate system are optimized. Specifically, the positions of no. 1~4 vibration isolators can be adjusted along the x and z axes of the IMU centroid coordinate system, and the positions of no. 5~8 vibration isolators can only be adjusted along the z axis, as shown in Figure 6.
Figure 6 is along the y-axis of the centroid coordinate system, where the positions of no. 2, 4, 6, and 8 vibration isolators overlap the positions of no. 1, 3, 5, and 7 vibration isolators, respectively. The design parameters in Figure 6 are x i ( i = 1 ~ 4 ) , z i ( i = 1 ~ 8 ) . The initial values of these design parameters are given in Table 1. The range of these parameters are as follows:
(1)
x i ϵ [ 80   100 ] ,   i = 1 ~ 4 , ( u n i t :   m m )   ;
(2)
z i ϵ [ 43   63 ] ,   i = 1 , 2 , 5 , 6 , ( u n i t :   m m )   ;
(3)
z i ϵ [ 51 31 ] ,   i = 3 , 4 , 7 , 8 , ( u n i t :   m m ) .

4.3. Determination of Objective Function

The excitation sources of the frequency response function are set to the torsional excitations around the dither directions of no. 1~3 RLGs to simulate the mechanical dithers of RLGs (i.e., the angular accelerations of the dither cavities of RLGs around the x, y, and z axes of the centroid coordinate system in the Figure 5). The response directions are set to the six DOFs for each of the RLGs.
Under each RLG excitation, and in six DOFs of each RLG, only one response is in the same direction, and the other five responses are in different directions (hereinafter referred to as coupled frequency response), as shown in Table 2. H x r x , H y r y   and   H z r z   in the table represent the responses in the same direction, and the rest are coupled frequency responses.
The optimization objective is set to minimize the ratio of the coupled frequency response to the same direction frequency response, that is, the following:
{ M i n ( H x x H x r x ) , M i n ( H x y H x r x ) , M i n ( H x z H x r x ) , M i n ( H x r y H x r x ) , M i n ( H x r z H x r x ) , M i n ( H y x H y r y ) , M i n ( H y y H y r y ) , M i n ( H y z H y r y ) , M i n ( H y r x H y r y ) , M i n ( H y r z H y r y ) , M i n ( H z x H z r z ) , M i n ( H z y H z r z ) , M i n ( H z z H z r z ) , M i n ( H z r x H z r z ) , M i n ( H z r y H z r z ) ,
According to Equation (27), the final optimization objective is the sum of the above 15 single objectives, which is the following:
F = M i n ( β 1 i = x , y , z j = x , y , z H i j H i r i + β 2 i = x , y , z , i j j = x , y , z H i r j H i r i ) ,
where β 1 and β 2 are the weighting factors.

4.4. Optimization Results and Discussion

The optimized positions of each vibration isolator relative to the IMU centroid are given in Table 3, in which the bold figures are the optimal values.
The frequency response characteristics of optimized IMU structure and the initial structure are compared, as shown in Figure 7. It shows the 3 × 3 array of the frequency response curves. The diagonal and non-diagonal terms are the frequency response curves in the same directions and coupled directions, respectively. The blue curves are for the initial IMU, and the red ones are for the optimized IMU with the optimal positions of vibration isolators. The blue circles in the figure represent the coupled frequencies in each coupled frequency response curve.
According to the IMU structure configuration in Figure 5 and the results of the frequency response analyses in Figure 7, it can be seen that the coupled vibration of the system is well attenuated after optimizing the positions of the vibration isolators relative to the IMU centroid. In Figure 7, the figures in the first row are the responses under a torsional excitation along the x direction, and the resonant peaks of 570 Hz and 670 Hz are the coupled response peaks. The figures in the second row are the response under a torsional excitation along the y direction, and the resonant peaks of 490 Hz and 670 Hz are the coupled response peaks. The figures in the third row are the responses under a torsional excitation along the z direction, and the resonant peaks of 490 Hz and 570 Hz are the coupled response peaks.
The amplitudes of each of the coupled response peaks before and after the IMU structure optimization are shown in Table 4.
As shown in Table 4, the coupled vibration of the system based on the previous IMU design principle is further attenuated by the IMU decoupling optimization method proposed in this paper. The transmissibility of coupled vibration is decreased by 6 dB at least, and by 20 dB at most.

5. Experimental Verification of IMU Dynamic Decoupling Optimization Method

In this section, the vibration response tests and RLGs’ accuracy tests under the dither excitations of RLGs were done to demonstrate the IMU dynamic decoupling optimization method.

5.1. Testing Object

Two sets of sensor supporting structures are manufactured according to the above two design methods, respectively. One is based on the dynamic decoupling conditions of Equations (24) and (25), which is called as the IMU structure before optimization. The other is based on the decoupling optimization method proposed in this paper, which is named as the IMU structure after optimization. The distribution positions of the vibration isolators of the two sets of sensor supporting structures are given in Table 1 and Table 3, respectively, and the other physical properties are almost the same, which are shown in Table 5.
By selecting the same group of RLGs as inertial sensors and the same type of rubber dampers as vibration isolators, two sets of IMU prototypes are assembled as the test objects.
The nominal parameters of the selected three RLGs in the test are given in Table 6, where the accuracy of the RLGs is regarded as the standard deviation derived from the data sequence formed by the mean values of the measuring signals per 100 s in the length of 3 h.

5.2. Testing Scheme

In different IMU structures, the accuracy of three RLGs are tested; meanwhile, the coupled vibration responses of each IMU structure are obtained with the vibration test system. The corresponding measured accuracy of the RLGs and the vibration responses of each IMU structure are compared to verify the effectiveness of the IMU dynamic decoupling optimization proposed in this paper.
The same set of hardware modules to meet the normal operation requirements of each RLG in IMU are used in the test, including the IMU power module (supplying 15 V to each RLG) and signal processing module (resampling and outputting one datum per second for each RLG). Meanwhile, the B&K vibration test system made in Denmark is used to obtain the vibration response of the sensor supporting structure when three RLGs dither simultaneously. The main parameters of B&K vibration test system and the environment of the test are shown in Table 7.
The test configuration is shown in Figure 8 and the signal flow diagram in the test is shown in Figure 9.
The three-axis accelerometers are attached to the vertices of each sensor supporting structure in the IMU, and the serial number of each accelerometer is shown in Figure 10, as well as its coordinate system. Simultaneously, the axes of each accelerometer are respectively required to be along with the sensitive axes of the RLGs in IMU. In Figure 10, the blue arrows are the sensitive axes of the three RLGs in IMU, and the red arrows are the measuring axes of eight accelerometers. This layout method makes that each accelerometer could be easily sensitive to the linear acceleration owing to the dither of three RLGs in IMU. For example, the x-axis of no. 1 accelerometer is easily sensitive to the dither components of 490 Hz and 670 Hz, while the dither components of the frequency response of 570 Hz will be considered as coupled vibration. In this way, the frequency response of each measuring axis of the eight accelerometers includes two frequency components along with the excitation direction (hereinafter referred to as the same direction response), and one frequency component orthogonal to the excitation direction (that is, the coupled response), as shown in Table 8.

5.3. Testing Results

The acceleration responses are collected in the range of 0~800 Hz, and the frequency resolution is set to 0.5 Hz. The total 24-axis vibration frequency responses of eight accelerometers under three gyro dither excitations in each IMU are obtained. For the sake of the length of this paper, it only shows the measured vibration frequency response curves of the no. 1 and no. 8 accelerometers in the three axes, as shown in Figure 11. The results of the resonant peak amplitude at the coupled vibration frequency before and after IMU optimization are shown in Table 9.
The test results in Table 9 show that, after the dynamic decoupling optimization, the 24 groups of acceleration responses in IMU are effectively attenuated, among them, 20 groups are more than 50% attenuated, and the coupled vibration amplitude in the optimized IMU is attenuated by 93.4% compared with the original IMU at most. Therefore, compared with the previous IMU design principle, the coupled vibration could be further suppressed by the dynamic decoupling optimization method proposed in this paper. In the whole optimization process, the positions of the vibration isolators relative to the IMU centroid are only adjusted, which relaxes the strict requirements of the mass distribution of IMU in the previous design principles, such as referred in the works of [19,20,21]. Thus, the method proposed in this paper is more suitable for the practical engineering application.
The measured accuracy of the RLGs in IMU before and after optimization are compared, as shown in Figure 12. After three hours of testing, the accuracy of the RLGs after IMU optimization was improved compared with that before optimization. The reason is that the coupled vibration caused by the simultaneous dither of three RLGs is effectively attenuated after the IMU decoupling optimization, which makes each RLG work more stably in IMU.
The repeatability of the accuracy results of the RLGs in two sets of IMU structures is verified. Through five repeated tests, the comparative results of the measured accuracy of the three RLGs in two sets of IMU structures are shown in Table 10.
From Table 10, except for a few results are shown as small perturbations, most of the measured accuracy results are improved more than 1 × 10−4 °/h, which could be a good promotion for a long-endurance navigation. Through the above repeatability of tests, the effectiveness of the proposed method based on the dynamic modeling and decoupling design of IMU structure is verified.

6. Conclusions

In this paper, the multi-rigid body dynamic model of IMU system is modeled based on the topology analysis method subject to the constraint relations, and the multi-DOF coupled vibration behavior in IMU is clearly discussed under the RLGs’ dither excitation. Simultaneously, on the basis of the IMU dynamic model and focused on the vibration frequency response, the new dynamic decoupling method is proposed to optimize the sensor supporting structure of IMU by adjusting the positions of vibration isolators relative to the centroid of IMU. The testing results show that, compared with the previous optimal design principle through adjusting the mass distribution of IMU, the proposed dynamic decoupling optimization method could further suppress the coupled vibration and improve the measuring accuracy of the RLGs in the IMU system. This paper provides a more effective design idea for IMU to attenuate the coupled dither excitation of RLGs and ensure the accuracy of in IMU.

Author Contributions

The idea and the relevant theoretical derivation in the paper are proposed by F.F., the experimental tests in the paper are achieved by W.Z., and relevant experimental data in the paper are processed by F.F., W.Z. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Chow, W.W.; Gea-Banacloche, J.; Pedrotti, L.M.; Sanders, V.E.; Schleich, W.; Scully, M.O. The ring laser gyro. Rev. Mod. Phys. 1985, 57, 61–104. [Google Scholar] [CrossRef]
  2. Faucheux, M.; Fayoux, D.; Roland, J.J. The ring laser gyro. J. Opt. 1988, 19, 101–115. [Google Scholar] [CrossRef]
  3. Yu, A.; Siddiqui, A.S. Theory of a novel high sensitivity optical fibre gyroscope. IEE Proc. J. 1993, 140, 150–156. [Google Scholar] [CrossRef]
  4. Lef, E.; Re Heacute, C. Fundamentals of the Interferometric Fiber-Optic Gyroscope. Opt. Rev. 1997, 4, 20–27. [Google Scholar] [CrossRef]
  5. Perelyaev, S.E.; Bodunov, B.P.; Bodunov, S.B. Solid-State Wave Gyroscope: A New-Generation Inertial Sensor. In Proceedings of the 2017 24th Saint Petersburg International Conference on Integrated Navigation Systems (ICINS), Saint Petersburg, Russia, 29–31 May 2017. [Google Scholar]
  6. Lynch, D.D. Hemispherical resonator gyro. Inertial Technology for the Future. IEEE Trans. Aerosp. Electron. Syst. 1984, 20, 414–444. [Google Scholar]
  7. Culshaw, B. The optical fibre Sagnac interferometer: An overview of its principles and applications. Meas. Sci. Technol. 2005, 17, 1–16. [Google Scholar] [CrossRef]
  8. Ragot, V.; Remillieux, G. A new control mode greatly improving performance of axisymmetrical vibrating gyroscopes. Gyroscopy Navig. 2011, 2, 229. [Google Scholar] [CrossRef]
  9. Liu, K.; Zhang, W.; Chen, W.; Li, K.; Dai, F.; Cui, F.; Wu, X.; Ma, G.; Xiao, Q. TOPICAL REVIEW: The development of micro-gyroscope technology. J. Micromech. Microeng. 2009, 19. [Google Scholar] [CrossRef]
  10. Dell’Olio, F.; Tatoli, T.; Ciminelli, C.; Armenise, M.N. Recent advances in miniaturized optical gyroscopes. J. Eur. Opt. Soc. Rapid Publ. 2014, 9, 164–167. [Google Scholar] [CrossRef] [Green Version]
  11. Dell’Olio, F.; Indiveri, F.; Innone, F.; Russo, P.D.; Ciminelli, C.; Armenise, M.N. System test of an optoelectronic gyroscope based on a high Q-factor InP ring resonator. Opt. Eng. 2014, 53, 127104. [Google Scholar] [CrossRef]
  12. Tang, J.X. Research and Design for Dither Bias System of Mechanically Dithered Ring Laser Gyroscopes. Master’s Thesis, National Defense Science and Technology University, Changsha, China, 2000. [Google Scholar]
  13. Chen, A.; Li, J.; Chu, Z. Dither Signal Removal of Ring Laser Gyro POS Based on Combined Digital Filter. In Proceedings of the 2012 8th IEEE International Symposium on Instrumentation and Control Technology (ISICT), London, UK, 11–13 July 2012; pp. 178–182. [Google Scholar]
  14. Wang, K.; Gu, Q. Optimization of Dither Amplitude and Frequency of RLG. In Proceedings of the 2002 IEEE Position Location and Navigation Symposium (IEEE Cat. No.02CH37284), Palms Springs, CA, USA, 15–18 April 2002; pp. 277–282. [Google Scholar]
  15. Fan, Z.; Luo, H.; Lu, G.; Hu, S. Direct dither control without external feedback for ring laser gyro. Opt. Laser Technol. 2012, 44, 767–770. [Google Scholar] [CrossRef]
  16. Song, S.-W.; Lee, J.-C.; Hong, S.-K.; Chwa, D. New random walk reduction algorithm in ring laser gyroscopes. J. Opt. 2010, 12, 11. [Google Scholar] [CrossRef]
  17. Fan, Z.; Luo, H.; Hu, S.; Xiao, G. Research on lock-in correction for mechanical dithered ring laser gyro. Opt. Eng. 2011, 50, 1–6. [Google Scholar] [CrossRef]
  18. Zheng, Z.; Han, S.; Zheng, K. An eight-position self-calibration method for a dual-axis rotational Inertial Navigation System. Sens. Actuators A Phys. 2015, 232, 39–48. [Google Scholar] [CrossRef]
  19. Kim, K.; Park, C.G. Drift error analysis caused by RLG dither axis bending. Sens. Actuators A Phys. 2007, 133, 425–430. [Google Scholar] [CrossRef]
  20. Wang, L.; Wu, W.; Li, G.; Pan, X.; Yu, R. Ring Laser Gyro G-Sensitive Misalignment Calibration in Linear Vibration Environments. Sensors 2018, 18, 601. [Google Scholar] [CrossRef] [Green Version]
  21. Tu, Y.; Yang, G.; Cai, Q.; Wang, L.; Yin, H. Dynamical analysis and experimental verification of deviation angles caused by rubber dampers deformation in high precision mechanically dithered RLG dual-axis RINS. Mech. Syst. Signal Process. 2019, 126, 553–567. [Google Scholar] [CrossRef]
  22. Lahham, J.I.; Wigent, D.J.; Coleman, A.L. Tuned Support Structure for Structure-Borne Noise Reduction of Inertial Navigator with Dithered Ring Laser Gyros (RLG). In Proceedings of the IEEE 2000 Position Location and Navigation Symposium (Cat. No.00CH37062), San Diego, CA, USA, 13–16 March 2000; pp. 419–428. [Google Scholar]
  23. Li, G.; Wei, G.; Xie, Y.P.; Zhang, P.F.; Yu, X.D.; Long, X.W. Isotropic Design Method of Suspension System of Dithered RLG Strapdown Inertial Measurement Unit. Adv. Mater. Res. 2014, 1044–1045, 788–797. [Google Scholar] [CrossRef]
  24. Li, J.; Fang, J.; Ge, S.S. Kinetics and Design of a Mechanically Dithered Ring Laser Gyroscope Position and Orientation System. IEEE Trans. Instrum. Meas. 2013, 62, 210–220. [Google Scholar] [CrossRef]
  25. Cheng, J.; Fang, J.; Wu, W.; Li, J. Optimized design method of vibration isolation system in mechanically dithered RLG POS based on motion decoupling. Measurement 2014, 48, 314–324. [Google Scholar] [CrossRef]
  26. Gao, C.; Wang, Q.; Wei, G.; Long, X. Research on the optimal structure configuration of dither RLG used in skewed redundant INS. In Proceedings of the SPIE Defense + Security, Baltimore, MD, USA, 13 May 2016. [Google Scholar]
  27. Choi, K.K.; Nam-Ho, K. Structural Sensitivity Analysis and Optimization 1: Linear Systems; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  28. Stepan, G.; Kiss, A.K.; Ghalamchi, B.; Sopanen, J.; Bachrathy, D. Chatter avoidance in cutting highly flexible workpieces. CIRP Ann. 2017, 66, 377–380. [Google Scholar] [CrossRef]
  29. Mohanty, S.; Dwivedy, S.K. Nonlinear dynamics of piezoelectric-based active nonlinear vibration absorber using time delay acceleration feedback. Nonlinear Dyn. 2019, 98, 1465–1490. [Google Scholar] [CrossRef]
  30. Haris, A.; Motato, E.; Mohammadpour, M.; Theodossiades, S.; Rahnejat, H.; O’Mahony, M.; Vakakis, A.F.; Bergman, L.A.; McFarland, D.M. On the effect of multiple parallel nonlinear absorbers in palliation of torsional response of automotive drivetrain. Int. J. Non-Linear Mech. 2017, 96, 22–35. [Google Scholar] [CrossRef] [Green Version]
  31. Tao, Y.; Rui, X.; Yang, F.; Hao, B. Development of a MRE isolation system for strapdown inertial measurement unit. Mech. Syst. Signal Process. 2019, 117, 553–568. [Google Scholar] [CrossRef]
  32. Denice, M.; Mark, J.; Tazartes, D. Isolator and Assembly Configuration. U.S. Patent US20030167863A1, 29 March 2005. [Google Scholar]
  33. Jiang, W.; Chen, X.; Yan, T. Symbolic formulation of rigid multibody systems for vibration analysis based on matrix transformation. Chin. J. Mech. Eng. 2008, 44, 54–60. [Google Scholar] [CrossRef]
  34. Ghiringhelli, G.L.; Terraneo, M. Analytically driven experimental characterisation of damping in viscoelastic materials. Aerosp. Sci. Technol. 2015, 40, 75–85. [Google Scholar] [CrossRef]
  35. Ewins, D.J. Modal Testing: Theory and Practice; Research Studies Press: Letchworth, UK, 1984. [Google Scholar]
Figure 1. Structure diagram of laser gyro inertial measurement unit (IMU). RLG, ring laser gyro.
Figure 1. Structure diagram of laser gyro inertial measurement unit (IMU). RLG, ring laser gyro.
Sensors 20 00111 g001
Figure 2. The schematic diagram of IMU dynamic model.
Figure 2. The schematic diagram of IMU dynamic model.
Sensors 20 00111 g002
Figure 3. The structure of the vibration isolator.
Figure 3. The structure of the vibration isolator.
Sensors 20 00111 g003
Figure 4. The constraint topological form of IMU.
Figure 4. The constraint topological form of IMU.
Sensors 20 00111 g004
Figure 5. The IMU structure and the serial number of vibration isolators and RLGs.
Figure 5. The IMU structure and the serial number of vibration isolators and RLGs.
Sensors 20 00111 g005
Figure 6. The optimization design parameters of IMU.
Figure 6. The optimization design parameters of IMU.
Sensors 20 00111 g006
Figure 7. Comparison with the frequency response characteristics before and after IMU optimization.
Figure 7. Comparison with the frequency response characteristics before and after IMU optimization.
Sensors 20 00111 g007
Figure 8. The configuration of the experiment.
Figure 8. The configuration of the experiment.
Sensors 20 00111 g008
Figure 9. The signal flow diagram in the IMU test. (PC in the figure means Personal Computer)
Figure 9. The signal flow diagram in the IMU test. (PC in the figure means Personal Computer)
Sensors 20 00111 g009
Figure 10. The layout of the measuring points of the three-axis accelerometers in the vibration test.
Figure 10. The layout of the measuring points of the three-axis accelerometers in the vibration test.
Sensors 20 00111 g010
Figure 11. The comparisons of the measured coupled vibration responses of no. 1 and 8 accelerometers before and after IMU optimization. (a) The measured accelerations of x-axis of no. 1 accelerometer before and after optimization; (b) The measured accelerations of y-axis of no. 1 accelerometer before and after optimization; (c) The measured accelerations of z-axis of no. 1 accelerometer before and after optimization; (d) The measured accelerations of x-axis of no. 8 accelerometer before and after optimization; (e) The measured accelerations of y-axis of no. 8 accelerometer before and after optimization; (f) The measured accelerations of z-axis of no. 8 accelerometer before and after optimization.
Figure 11. The comparisons of the measured coupled vibration responses of no. 1 and 8 accelerometers before and after IMU optimization. (a) The measured accelerations of x-axis of no. 1 accelerometer before and after optimization; (b) The measured accelerations of y-axis of no. 1 accelerometer before and after optimization; (c) The measured accelerations of z-axis of no. 1 accelerometer before and after optimization; (d) The measured accelerations of x-axis of no. 8 accelerometer before and after optimization; (e) The measured accelerations of y-axis of no. 8 accelerometer before and after optimization; (f) The measured accelerations of z-axis of no. 8 accelerometer before and after optimization.
Sensors 20 00111 g011aSensors 20 00111 g011bSensors 20 00111 g011c
Figure 12. The comparisons of the measured accuracy of three RLGs in two sets of IMU structures. (a) The measured accuracy of gyro no. 1 with 490 Hz dither frequency; (b) the measured accuracy of gyro no. 2 with 570 Hz dither frequency; (c) the measured accuracy of gyro no. 3 with 670 Hz dither frequency.
Figure 12. The comparisons of the measured accuracy of three RLGs in two sets of IMU structures. (a) The measured accuracy of gyro no. 1 with 490 Hz dither frequency; (b) the measured accuracy of gyro no. 2 with 570 Hz dither frequency; (c) the measured accuracy of gyro no. 3 with 670 Hz dither frequency.
Sensors 20 00111 g012aSensors 20 00111 g012b
Table 1. The initial positions of vibration isolators relative to the inertial measurement unit (IMU) centroid.
Table 1. The initial positions of vibration isolators relative to the inertial measurement unit (IMU) centroid.
Serial Numberx (mm)y (mm)z (mm)
1100.1−87.662.8
2100.182.462.8
3100.1−87.6−51.2
4100.182.4−51.2
5−88.9−87.662.8
6−88.982.462.8
7−88.9−87.6−51.2
8−88.982.4−51.2
Table 2. The components of excitation and response in IMU.
Table 2. The components of excitation and response in IMU.
Excitation DirectionDirection of Frequency Response
Translational ResponseRevolutional Response
X Y Z X Y Z
Revolution Excitation of X   H x x     H x y     H x z   H x r x   H x r y   H x r z
Revolution Excitation of Y H y x H y y H y z H y r x H y r y H y r z  
Revolution Excitation of Z H z x H z y   H z z   H z r x   H z r y   H z r z  
Table 3. The optimized positions of the vibration isolators relative to the IMU centroid.
Table 3. The optimized positions of the vibration isolators relative to the IMU centroid.
Serial NumberIn X Direction (mm)In Y Direction (mm)In Z Direction (mm)
193.3−87.662.1
291.382.462.1
396.3−87.6−49.5
496.382.4−46.5
5−88.9−87.662.6
6−88.982.466.6
7−88.9−87.6−47.5
8−88.982.4−50.5
Table 4. Comparison with the coupled response peaks of IMU structures before and after optimization.
Table 4. Comparison with the coupled response peaks of IMU structures before and after optimization.
Situation of IMU StructureCoupled Response Peaks in Each Coupled Direction (dB)
H x r y     H x r z   H y r x   H y r z   H z r x   H z r y  
Before IMU optimization−116−102−110−102−91−97
After IMU optimization−122−120−117−122−108−117
Table 5. The physical properties of two sets of sensor supporting structures. RLG, ring laser gyro.
Table 5. The physical properties of two sets of sensor supporting structures. RLG, ring laser gyro.
Physical PropertyIMU Structure before OptimizationIMU Structure after Optimization
MaterialAluminum alloyAluminum alloy
Weight (kg)2.89 kg2.88 kg
Inertia ( kgmm 2 )
(Including three RLGs)
[ 6.56 × 10 4 1.09 × 10 2 6.61 × 10 2 1.09 × 10 2 7.75 × 10 4 3.55 × 10 2 6.61 × 10 2 3.55 × 10 2 7.54 × 10 4 ] [ 6.64 × 10 4 1.24 × 10 2 5.69 × 10 2 1.24 × 10 2 7.84 × 10 4 2.09 × 10 2 5.69 × 10 2 2.09 × 10 2 7.62 × 10 4 ]
Table 6. The nominal parameters of the three RLGs.
Table 6. The nominal parameters of the three RLGs.
Nominal ParametersRLG No. 1RLG No. 2RLG No. 3
Bias drift0.075 °/h0.11 °/h0.106 °/h
Angle random walk 0.00065   ° / h   0.00057   ° / h   0.00041   ° / h  
Accuracy0.0026 °/h0.0027 °/h0.0019 °/h
Dither frequency493 Hz568 Hz667 Hz
Table 7. The main parameters of B&K system and the environment of the test.
Table 7. The main parameters of B&K system and the environment of the test.
Parameters of
B&K System
Analyzing BandwidthDC-51.2 kHz
Dynamic Range160 dB
Input Voltage Range±10 V
A/D24 bit
Sensitivity of Acceleration10 mV/g
Environment
of Test
Location30.392 °N, 114.349 °E
Temperature25 ± 2 °C
Table 8. Corresponding relations of the response components of each measuring point of each accelerometer.
Table 8. Corresponding relations of the response components of each measuring point of each accelerometer.
Measuring PointsMeasuring AxisSame Direction Response Frequency Components (Hz)Coupled Response Frequency Components (Hz)
1~8x490670570
y570670490
z490570670
Table 9. The comparison result of the dynamic decoupling effect before and after IMU optimization.
Table 9. The comparison result of the dynamic decoupling effect before and after IMU optimization.
Number of Measuring PointDirectionResonant Peak Amplitude (m/s2)Reduction Ratio of Amplitude
Before IMU OptimizedAfter IMU Optimized
1X-axis1.9760.20889.5%
Y-axis0.6230.04193.4%
Z-axis0.8460.33854.1%
2X-axis0.2010.17512.9%
Y-axis0.3050.10465.9%
Z-axis0.7040.32653.7%
3X-axis0.6310.35443.9%
Y-axis0.7540.22370.4%
Z-axis0.4510.13969.2%
4X-axis0.6550.31851.5%
Y-axis0.1890.04775.1%
Z-axis0.5940.2459.6%
5X-axis0.5710.16870.6%
Y-axis5.6312.70252%
Z-axis0.4550.25444.2%
6X-axis0.2680.05878.4%
Y-axis0.5390.1474%
Z-axis0.3180.20834.6%
7X-axis1.0340.16683.9%
Y-axis0.4630.10178.2%
Z-axis2.0690.64968.6%
8X-axis1.4550.17687.9%
Y-axis0.4330.11174.4%
Z-axis1.8390.50772.4%
Table 10. The repeatability of the measured accuracy of the three RLGs in two sets of IMU structures.
Table 10. The repeatability of the measured accuracy of the three RLGs in two sets of IMU structures.
Times of TestsSerial Number of RLGMeasured Accuracy (°/h)Trend of Accuracy
Before IMU OptimizedAfter IMU Optimized
112.77 × 10−32.66 × 10−3Improvement
22.83 × 10−32.81 × 10−3Small perturbation
32.34 × 10−31.94 × 10−3Improvement
212.84 × 10−32.65 × 10−3Improvement
22.91 × 10−32.93 × 10−3Small perturbation
32.31 × 10−31.99 × 10−3Improvement
312.71 × 10−32.54 × 10−3Improvement
22.81 × 10−32.69 × 10−3Improvement
32.11 × 10−32.01 × 10−3Improvement
412.83 × 10−32.68 × 10−3Improvement
22.72 × 10−32.75 × 10−3Small perturbation
32.15 × 10−31.98 × 10−3Improvement
512.81 × 10−32.69 × 10−3Improvement
22.88 × 10−32.75 × 10−3Improvement
32.21 × 10−31.93 × 10−3Improvement

Share and Cite

MDPI and ACS Style

Fang, F.; Zeng, W.; Li, Z. Coupled Dynamic Analysis and Decoupling Optimization Method of the Laser Gyro Inertial Measurement Unit. Sensors 2020, 20, 111. https://doi.org/10.3390/s20010111

AMA Style

Fang F, Zeng W, Li Z. Coupled Dynamic Analysis and Decoupling Optimization Method of the Laser Gyro Inertial Measurement Unit. Sensors. 2020; 20(1):111. https://doi.org/10.3390/s20010111

Chicago/Turabian Style

Fang, Fang, Wenhui Zeng, and Zilong Li. 2020. "Coupled Dynamic Analysis and Decoupling Optimization Method of the Laser Gyro Inertial Measurement Unit" Sensors 20, no. 1: 111. https://doi.org/10.3390/s20010111

APA Style

Fang, F., Zeng, W., & Li, Z. (2020). Coupled Dynamic Analysis and Decoupling Optimization Method of the Laser Gyro Inertial Measurement Unit. Sensors, 20(1), 111. https://doi.org/10.3390/s20010111

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