Ring Laser Gyro G-Sensitive Misalignment Calibration in Linear Vibration Environments

The ring laser gyro (RLG) dither axis will bend and exhibit errors due to the specific forces acting on the instrument, which are known as g-sensitive misalignments of the gyros. The g-sensitive misalignments of the RLG triad will cause severe attitude error in vibration or maneuver environments where large-amplitude specific forces and angular rates coexist. However, g-sensitive misalignments are usually ignored when calibrating the strapdown inertial navigation system (SINS). This paper proposes a novel method to calibrate the g-sensitive misalignments of an RLG triad in linear vibration environments. With the SINS is attached to a linear vibration bench through outer rubber dampers, rocking of the SINS can occur when the linear vibration is performed on the SINS. Therefore, linear vibration environments can be created to simulate the harsh environment during aircraft flight. By analyzing the mathematical model of g-sensitive misalignments, the relationship between attitude errors and specific forces as well as angular rates is established, whereby a calibration scheme with approximately optimal observations is designed. Vibration experiments are conducted to calibrate g-sensitive misalignments of the RLG triad. Vibration tests also show that SINS velocity error decreases significantly after g-sensitive misalignments compensation.


Introduction
Strapdown inertial navigation systems (SINS) has been widely used in many applications. In order to improve the SINS accuracy, many techniques have been adopted, such as the zero velocity updated (ZUPT) technique [1], the data processing technique [2], and the Kalman filter technique [3,4]. A SINS based on Micro-Electro-Mechanical Syetem (MEMS) has also been used in some highly dynamic environments [5]. In high accuracy applications, the ring laser gyro (RLG) has become a common instrument in SINS for spacecraft and other aerospace systems [6]. In these applications, SINS may be employed in harsh working environments, e.g., vibration environments, maneuver motions with large-amplitude specific forces and angular rates, or other highly dynamic environments. The laser block of an RLG is mounted on the inertial measurement unit (IMU) base through a dither suspension structure. However, because the transverse stiffness of the RLG dither axis is limited [7], the RLG dither axis will bend and exhibit errors due to the specific forces acting on the instrument, which are known as g-sensitive misalignments of the gyros. Namely, the laser block of the RLG will tilt under the influence of specific forces. The g-sensitive error is usually taken into consideration when calibrating MEMS gyros [8,9]. However, the calibration of g-sensitive misalignments regarding high-accuracy RLGs is usually ignored for general applications [10,11]. In aerospace applications, due to the presence of great acceleration and angular rate, g-sensitive misalignments of the RLG triad will introduce equivalent gyro drift, causing further attitude error and velocity error. Hence, g-sensitive misalignments compensation must be taken into consideration in the SINS navigation software.
In order to reduce the influence of g-sensitive misalignments in a vibration environment, parameter optimization has been applied to the RLG dither suspension structure [12], and the vibration isolation structure of the SINS has been elaborately designed [13]. However, the vibration tests performed on the laser SINS for a Cyclone-4 launch vehicle still reveal that the performance degrades regardless of careful vibration isolation design, and equivalent gyro drift occurs in the SINS [14]. In References [15,16], a Kalman filter and a multi-position calibration method are adopted to estimate the RLG g-dependent biases; however, the method to calibrate g-sensitive misalignments of the RLG triad is not discussed. The influence of g-sensitive misalignment errors upon the estimation accuracy of SINS error parameters is verified via numerical simulations in Reference [17]. The drift errors caused by RLG dither axis bending are analyzed in Reference [18]; however, only theoretical and simulation results are discussed. The calibration and compensation methods are not given in Reference [18]. A calibration method is proposed in Reference [19], in which a sequence of coning tests on a high-precision three-axis turntable are conducted to estimate the g-sensitive misalignments of the RLGs caused by gravity; however, the procedures are very complicated and time-consuming. The calibration method proposed in Reference [18] is based on static tests and is only suitable in low-g (less than 1 g) environments. In addition, actual calibration experiments and vibration tests are not conducted.
This paper presents a novel method to calibrate the g-sensitive misalignments of an RLG triad in linear vibration environments. The g-sensitive misalignments of the RLG triad is deduced and the mechanism is analyzed. An equivalent gyro drift model is established to express the error caused by g-sensitive misalignments, which builds a relationship between the equivalent gyro drift and the specific forces as well as the angular rate. Further, the attitude error caused by g-sensitive misalignments is analyzed, whereby a calibration method based on vibration experiments is designed. To fully excite the errors caused by g-sensitive misalignments, vibration experiments are elaborately designed. With the SINS attached to a linear vibration bench through outer rubber dampers, rocking of the SINS can occur when the linear vibration bench is employed. Therefore, linear vibration environments can be created to simulate the harsh environment during aircraft flight. In this condition, g-sensitive misalignments of the RLG triad will cause severe attitude errors. By analyzing the mathematical model of g-sensitive misalignments of an RLG triad, the attitude error observation equations are deduced. With error parameter sensitivity maximized, a calibration scheme with approximately optimal observations is designed. Experiments are also conducted to validate the method proposed in the paper.

Mathematical and Model
Because of the dither suspension structure, the RLG dither axis will bend and exhibit errors due to the specific forces acting on the instrument, which are known as g-sensitive misalignments of the gyros. The bending of X RLG dither axis under specific forces, depicted in Figure 1, is taken as an example to explain the mechanism of g-sensitive misalignments. As shown in Figure 1, the X RLG body coordinate frame-b g x is defined as follows: g xx -axis coincides with the X gyro sensitive axis, namely, the b x -axis of the SINS body frame-b; g xy and g xz -axes are parallel to the b y and b z -axes of the SINS body frame, respectively. The X RLG body coordinate frame-b g x and SINS body frame-b are right-hand orthogonal coordinate frames. They are nominally parallel. In general, the framex g b is parallel to the SINS body frame-b. However, if the specific forces act on the X RLG, the aforementioned two coordinate frames will no longer be parallel. There will be a small angle between them. A rotation vector denoted by x g  can be used to express the rotation from the frame-b to the framex g b , and is given by: where xz  and xy  respectively denote the X gyro-sensitive axis misalignment toward z b and y b -axes of the SINS due to specific forces acting on the X RLG. The quantity of x g  is a small angle.
The angular rate of X RLG expressed in the framex g b is given by: is the angular rate of X RLG expressed in the frame- ib ω is the true angular rate of the SINS expressed in the frameb , and C is the direction cosine matrix from the frameb to the frame- Based on small angle assumption, the direction cosine matrix from the frameb to the frame- C , can be expressed as Hence, Equation (2) can be rewritten as: where 3 I denotes a unit matrix, and In fact, only the first component of Equation (4)  . Hence, we obtain: In general, the frame-b g x is parallel to the SINS body frame-b. However, if the specific forces act on the X RLG, the aforementioned two coordinate frames will no longer be parallel. There will be a small angle between them. A rotation vector denoted by θ g x can be used to express the rotation from the frame-b to the frame-b g x , and is given by: where θ xz and θ xy respectively denote the X gyro-sensitive axis misalignment toward b z and b y -axes of the SINS due to specific forces acting on the X RLG. The quantity of θ g x is a small angle. The angular rate of X RLG expressed in the frame-b g x is given by: where ω b gx ib gx is the angular rate of X RLG expressed in the frame-b g x , ω b ib is the true angular rate of the SINS expressed in the frame-b, and C b gx b is the direction cosine matrix from the frame-b to the frame-b g x .
Based on small angle assumption, the direction cosine matrix from the frame-b to the frame-b g x , i.e., C b gx b , can be expressed as Hence, Equation (2) can be rewritten as: where I 3 denotes a unit matrix, and θ g x × denotes the skew symmetric matrix of vector θ g x . Rewriting Equation (3) in components form leads to another identical expression: In fact, only the first component of Equation (4) can be sensed by X RLG. The term ω b iby θ xy + ω b ibz θ xz denotes the equivalent gyro drift caused by the g-sensitive misalignments of X RLG. To clarify, we use ω b ibx to represent the angular rate sensed by X RLG, i.e., the first component of ω b gx ib gx . Hence, we obtain: Analogously, the angular rates sensed by Y RLG and Z RLG are respectively given by: where ω b iby and ω b ibz respectively denote the angular rate sensed by Y RLG and Z RLG, and θ ij (i = x, y, z; j = x, y, z; i = j) denotes the i gyro misalignment toward the j-axis due to specific forces acting on the i gyro.
According to Equations (5)-(7), the mathematical model for the g-sensitive misalignments of an RLG triad can be represented as: with where δω b ib denotes the equivalent gyro drift of the RLG triad expressed in the body frame-b, ω b ib denotes the angular rate sensed by the RLG triad, and M denotes the g-sensitive misalignments matrix of the RLG triad.
As shown in References [18,19], the g-sensitive misalignments in Equation (9) are proportional to the amount of specific force acting on the RLG, which can be expressed as: where f b i (i = x, y, z) denotes the i-axis specific force sensed by accelerometer triad, and τ ijk (i = x, y, z; j = x, y, z; k = x, y, z; i = j) denotes the flexure coefficient which accounts for i gyro misalignment toward the j-axis due to specific force along k-axis, whose unit is rad/(m/s 2 ) or arc-second/g (1 g ≈ 9.8 m/s 2 ).
Substituting Equations (9) and (10) into Equation (8) leads to: where ω b ib can be expressed in components form as (11) is the equivalent gyro drift induced by the g-sensitive misalignments of an RLG triad. Twelve flexure coefficients notated as a vector τ are contained in Equation (11), and τ is given by: The task of this paper is to estimate these parameters and compensate the corresponding errors. Substituting Equation (12) into Equation (11) and rearranging Equation (11) leads to: with

Calibration Method
In this section, a calibration algorithm based on vibration tests to estimate unknown parameters in Equation (13) is introduced, followed by a calibration scheme designed with approximately optimal observations.

Calibration Algorithm
Based on the attitude update algorithm, the direction cosine matrix expressing the rotation from the body frame-b to the inertial frame-i can be updated by the gyro output ω b ib . The attitude update algorithm is given by: where t 0 is the start time of vibration, and the inertial frame-i is formed by fixing the b-frame at the time t 0 in the inertial space, namely, i b(t 0 ). If there is no vibration acting on the SINS, the calculated attitude according to Equation (15) is without error. However, due to the linear vibration acting on the SINS, g-sensitive misalignments will occur. As a result, the attitude calculated by Equation (15) will exhibit great error introduced by equivalent gyro drift. The time derivative of attitude error is given by: where φ i denotes the attitude error expressed in the inertial frame-i. Substituting Equation (13) into Equation (16), taking the integral of Equation (16) from time t 0 to time t k , and rearranging it leads to: where t k is the end time of the vibration. Here, the initial attitude error in Equation (17) is ignorable, namely, φ i (t 0 ) ≈ 0 3×1 . According to Equation (17), τ can be estimated with attitude error observations. Calibration experiments can be conducted through following steps:

•
Step 1: The SINS is first at rest for a certain time to carry out static self-alignment. The static self-alignment phase provides the initial attitude, which is essential information to make the SINS work. The initial attitude is not influenced by g-sensitive misalignments.

•
Step 2: The linear vibration is applied to the SINS following the static period, and lasts a few minutes, e.g., 8-10 min. The SINS works on inertial navigation mode. Based on the initial attitude calculated by Step 1, the SINS attitude is updated to the end of vibration. During this period, by an elaborately designed fixed device, the attitude error of SINS is stimulated under the influence of g-sensitive misalignments.

•
Step 3: The SINS returns to the stationary state and performs static self-alignment again. This phase is only employed to provide the correct attitude reference of SINS. Based on the attitude reference calculated in Step 3, the error of attitude calculated by Step 2 can be observed and is given by: is the direction cosine matrix calculated by inertial navigation, and C b i (t k ) denotes the true direction cosine matrix at time t k calculated by Equation (19).
is calculated by the alignment process in Step 1, and C n(t k ) n(t 0 ) is calculated by Equation (20) with being the Earth's rotation rate, ω ie the Earth rate, and L the geographic latitude. Note that C n(t k ) n(t 0 ) is determined by Earth's rotation rate and vibration time, and there exists an explicit relationship with them.
At least 12 different attitude error observations are required to solve the unknown parameters in Equation (17), so several vibration experiments are needed. Figure 2 shows the calibration method process. Considering the limitation of the north-finding accuracy, only the horizontal components of the attitude error observations in Equation (18) are used. Equation (17) can be solved by least squares as follows: where the subscript 1, 2, ..., n denotes the vibration experiment number.
Sensors 2018, 18, 601 6 of 13 where ()  (17), so several vibration experiments are needed. Figure 2 shows the calibration method process. Considering the limitation of the north-finding accuracy, only the horizontal components of the attitude error observations in Equation (18) are used. Equation (17) can be solved by least squares as follows:  where the subscript 1, 2, …, n denotes the vibration experiment number.
The SINS firstly carries out static selfalignment to provide initial attitude reference The linear vibration is applied to the SINS and lasts 10 minutes. The attitude error is stimulated The SINS returns to the static state and performs realignment to provide attitude error observation If the observations are sufficient, the parameters can be solved. If not, repeat the previous steps

Optimal Observations
As mentioned in Reference [20], a good observation scheme not only allows all of the unknown parameters to be estimated, but also improves the accuracy of estimation. It is important to maximize the sensitivity of the attitude error observations with respect to these parameters.
Because the influence of the Earth's rotation rate is ignorable in a short time, i.e., C n(t 0 ) n(t) ≈ I 3 , and the attitude change induced by linear vibration satisfies small-angle approximation, i.e., (17) can be simplified as: Taking the partial derivatives of Equation (23) with respect to the components of τ leads to the Jacobian matrix, given as: We take X RLG as an example to illustrate how to obtain the optimal observations for its flexure coefficients. The first four columns of Equation (24) are the partial derivatives of φ i (t k ) with respect to τ xyx , τ xyy , τ xzx , and τ xzz .
Intuitively, if the linear vibration directions are respectively along X-axes as shown in Figure 3a,c, both the values of Equation (25) will reach their maximums accordingly. Specifically, they are approximately optimal observations for τ xyx and τ xzx . If the linear vibration direction is along the diagonal of the YZ plane as shown in Figure 3b, both the values of Equation (26) will reach their maximums. Namely, this is an approximately optimal observation for τ xyy and τ xzz .

Optimal Observations
As mentioned in Reference [20], a good observation scheme not only allows all of the unknown parameters to be estimated, but also improves the accuracy of estimation. It is important to maximize the sensitivity of the attitude error observations with respect to these parameters.
Because the influence of the Earth's rotation rate is ignorable in a short time, i.e.,  (17) can be simplified as: Taking the partial derivatives of Equation (23) with respect to the components of  leads to the Jacobian matrix, given as: We take X RLG as an example to illustrate how to obtain the optimal observations for its flexure coefficients. The first four columns of Equation (24) are the partial derivatives of () Intuitively, if the linear vibration directions are respectively along X-axes as shown in Figure 3a

Experiment Results
Linear vibration experiments are conducted to calibrate the g-sensitive misalignments of the RLG triad of a high-precision navigation-grade SINS. The SINS is specially designed for the aircraft. The SINS includes three mechanically dithered gyroscopes with a bias stability of 0.005 deg/h (1 ) and three quartz accelerometers with a bias stability of 20 μg (1 ). Before the vibration experiments, we first calibrated the scale factor constants, misalignments, and biases of gyroscopes and accelerometers in the laboratory. The scale factor of the RLG is very stable, with an error of less than 1 ppm, and even as small as 0.1 ppm [21]. Compared with the errors caused by g-sensitive misalignments, the scale factor error of the RLG is so small that it can be ignored. In addition, the size effect parameters and nonlinear coefficients including square coefficients and cross-coupling coefficients of accelerometers should also be calibrated in advance.
As shown in Figure 4, the SINS is directly mounted to a fixture, and the fixture is mounted to the vibration bench through four outer rubber dampers. The rubber dampers are located under the fixture. The fixture is L-shaped, and has two planes including a horizontal plane and a vertical plane. By fixing the SINS to the horizontal and vertical planes of the fixture, we can perform a vibration test on each axis of the SINS. In addition, because the inertial navigation system is a strapdown inertial navigation system, there is no special requirement for SINS installation accuracy on the vibration bench. The vibration bench is a THV710A linear vibration table, which is electro-dynamic. The sine trust and random trust of the vibration bench are 50.0 kN and 50.0 kN (Root Mean Square, RMS), respectively. The frequency range of the vibration bench is from 1 to 2500 Hz. The natural frequency of the selected outer rubber dampers is approximately in the bandwidth of 10~20 Hz. Compared with

Experiment Results
Linear vibration experiments are conducted to calibrate the g-sensitive misalignments of the RLG triad of a high-precision navigation-grade SINS. The SINS is specially designed for the aircraft. The SINS includes three mechanically dithered gyroscopes with a bias stability of 0.005 deg/h (1σ) and three quartz accelerometers with a bias stability of 20 µg (1σ). Before the vibration experiments, we first calibrated the scale factor constants, misalignments, and biases of gyroscopes and accelerometers in the laboratory. The scale factor of the RLG is very stable, with an error of less than 1 ppm, and even as small as 0.1 ppm [21]. Compared with the errors caused by g-sensitive misalignments, the scale factor error of the RLG is so small that it can be ignored. In addition, the size effect parameters and nonlinear coefficients including square coefficients and cross-coupling coefficients of accelerometers should also be calibrated in advance.
As shown in Figure 4, the SINS is directly mounted to a fixture, and the fixture is mounted to the vibration bench through four outer rubber dampers. The rubber dampers are located under the fixture. The fixture is L-shaped, and has two planes including a horizontal plane and a vertical plane. By fixing the SINS to the horizontal and vertical planes of the fixture, we can perform a vibration test on each axis of the SINS. In addition, because the inertial navigation system is a strapdown inertial navigation system, there is no special requirement for SINS installation accuracy on the vibration bench. The vibration bench is a THV710A linear vibration table, which is electro-dynamic. The sine trust and random trust of the vibration bench are 50.0 kN and 50.0 kN (Root Mean Square, RMS), respectively. The frequency range of the vibration bench is from 1 to 2500 Hz. The natural frequency of the selected outer rubber dampers is approximately in the bandwidth of 10~20 Hz. Compared with the SINS's own inner vibration isolators, the selected outer rubber dampers are far softer, so there will be a large-amplitude angular rate sensed by the RLG triad due to the vibration acceleration acting on the SINS. Vibration experiments are performed on the SINS with approximately optimal observations, as shown in Figure 3. In order to increase the product of specific force and angular rate regarding parameter sensitivity, the amplitude and frequency of vibration are set as 1.5 g and 20 Hz, which is close to the natural frequency of the outer rubber dampers. the SINS's own inner vibration isolators, the selected outer rubber dampers are far softer, so there will be a large-amplitude angular rate sensed by the RLG triad due to the vibration acceleration acting on the SINS. Vibration experiments are performed on the SINS with approximately optimal observations, as shown in Figure 3. In order to increase the product of specific force and angular rate regarding parameter sensitivity, the amplitude and frequency of vibration are set as 1.5 g and 20 Hz, which is close to the natural frequency of the outer rubber dampers.  Table 1. Taking yxy  as an example, whose value is 0.5108 arc-second/g, this means that 1 g specific force along the Y-axis causes about 0.5 arc-second of misalignment angle of the Y RLG dither axis toward the X-axis. Note that coefficients ijj  (i = x, y, z; j = x, y, z; i  j) account for the lateral anisoelasticity of the i RLG dither axis, whose values are close to each other. In addition, the sign of yxx  and yzz  is positive rather than negative due to the installation relationship of Y RLG.  To verify the effect of g-sensitive misalignments compensation, a linear vibration test is conducted. The SINS first carries out the alignment process for 15 min. Then the linear vibration along the Y-axis in the XY plane is applied to the SINS, lasting about 10 min. The SINS works on pure inertial navigation mode during the vibration. Figure 5 Table 1. Taking τ yxy as an example, whose value is 0.5108 arc-second/g, this means that 1 g specific force along the Y-axis causes about 0.5 arc-second of misalignment angle of the Y RLG dither axis toward the X-axis. Note that coefficients τ ijj (i = x, y, z; j = x, y, z; i = j) account for the lateral anisoelasticity of the i RLG dither axis, whose values are close to each other. In addition, the sign of τ yxx and τ yzz is positive rather than negative due to the installation relationship of Y RLG. To verify the effect of g-sensitive misalignments compensation, a linear vibration test is conducted. The SINS first carries out the alignment process for 15 min. Then the linear vibration along the Y-axis in the XY plane is applied to the SINS, lasting about 10 min. The SINS works on pure inertial navigation mode during the vibration. Figure 5 shows the details of f b y and ω b ibx within 0.5 s. The curves in this figure are plotted according to the IMU gyro output and accelerometer output. By inspection of the figure, we can obtain the frequency, amplitude, and phase information of f b y and ω b ibx . The frequency of f b y and ω b ibx is 20 Hz. The amplitudes of f b y and ω b ibx are about 1 g and 36 deg/s. Due to the damping effect of the outer rubber dampers, there is a phase difference between f b y and ω b ibx , whose value is about 3π/5. Hence, the representations of f b y and ω b ibx are given by: In this condition, the equivalent gyro drift regarding f b y and ω b ibx is given by: The direct component of Equation (28) is the constant equivalent gyro drift. The value is about 0.05 deg/h, which is a great value for high-accuracy SINS. The equivalent gyro drift will further introduce great velocity error, as shown in Figure 6. After g-sensitive misalignments compensation, the velocity error significantly decreases. The percentage of velocity error decrease is about 61%.
In this condition, the equivalent gyro drift regarding b y f and b ibx  is given by: The direct component of Equation (28) is the constant equivalent gyro drift. The value is about 0.05 deg/h, which is a great value for high-accuracy SINS. The equivalent gyro drift will further introduce great velocity error, as shown in Figure 6. After g-sensitive misalignments compensation, the velocity error significantly decreases. The percentage of velocity error decrease is about 61%. In order to further verify the effect of g-sensitive misalignments compensation, other 27 sets of vibration experiments including fixed frequency vibration and random vibration are performed on the SINS. The amplitude and frequency of fixed frequency vibration are listed in Table 2. For random vibration, 0.7 g root mean square (RMS) vibration acceleration is set with constant power spectral density (PSD) in a frequency range from 10 to 60 Hz. The vibration direction is along the diagonal of the SINS, which makes the acceleration act on the SINS' multiple axes at the same time. The SINS first carries out the alignment process for 15 min before the vibration test. Each vibration experiment lasts about 10 min. For the high-accuracy SINS, compared with the error caused by g-sensitive misalignments in vibration, the influence of other error sources, e.g., gyro bias and accelerometer bias, can be ignored within 10 min, which is much less than the Schuler period, i.e., 84.4 min. So, the major errors during the vibration for 10 min are the g-sensitive misalignments. The choice of vibration The direct component of Equation (28) is the constant equivalent gyro drift. The value is about 0.05 deg/h, which is a great value for high-accuracy SINS. The equivalent gyro drift will further introduce great velocity error, as shown in Figure 6. After g-sensitive misalignments compensation, the velocity error significantly decreases. The percentage of velocity error decrease is about 61%. In order to further verify the effect of g-sensitive misalignments compensation, other 27 sets of vibration experiments including fixed frequency vibration and random vibration are performed on the SINS. The amplitude and frequency of fixed frequency vibration are listed in Table 2. For random vibration, 0.7 g root mean square (RMS) vibration acceleration is set with constant power spectral density (PSD) in a frequency range from 10 to 60 Hz. The vibration direction is along the diagonal of the SINS, which makes the acceleration act on the SINS' multiple axes at the same time. The SINS first carries out the alignment process for 15 min before the vibration test. Each vibration experiment lasts about 10 min. For the high-accuracy SINS, compared with the error caused by g-sensitive misalignments in vibration, the influence of other error sources, e.g., gyro bias and accelerometer bias, can be ignored within 10 min, which is much less than the Schuler period, i.e., 84.4 min. So, the major errors during the vibration for 10 min are the g-sensitive misalignments. The choice of vibration  In order to further verify the effect of g-sensitive misalignments compensation, other 27 sets of vibration experiments including fixed frequency vibration and random vibration are performed on the SINS. The amplitude and frequency of fixed frequency vibration are listed in Table 2. For random vibration, 0.7 g root mean square (RMS) vibration acceleration is set with constant power spectral density (PSD) in a frequency range from 10 to 60 Hz. The vibration direction is along the diagonal of the SINS, which makes the acceleration act on the SINS' multiple axes at the same time. The SINS first carries out the alignment process for 15 min before the vibration test. Each vibration experiment lasts about 10 min. For the high-accuracy SINS, compared with the error caused by g-sensitive misalignments in vibration, the influence of other error sources, e.g., gyro bias and accelerometer bias, can be ignored within 10 min, which is much less than the Schuler period, i.e., 84.4 min. So, the major errors during the vibration for 10 min are the g-sensitive misalignments. The choice of vibration time is for error analysis. The SINS works on pure inertial navigation mode during the vibration. The evaluation index is the percentage of velocity error decrease after g-sensitive misalignments compensation. Figure 7 depicts the velocity error curves in one of the random vibration tests. After g-sensitive misalignments compensation, the velocity error decreases by 66.1%. time is for error analysis. The SINS works on pure inertial navigation mode during the vibration. The evaluation index is the percentage of velocity error decrease after g-sensitive misalignments compensation. Figure 7 depicts the velocity error curves in one of the random vibration tests. After g-sensitive misalignments compensation, the velocity error decreases by 66.1%. All of the percentages of velocity error decrease corresponding with the 27 sets of vibration experiments are listed in Table 2. A corresponding statistic pie chart is shown in Figure 8. The statistical results show that the proportion of velocity error decreasing by more than 30% in the total results reach about 82% after g-sensitive misalignments compensation. As shown in Table 2 and Figure 8, though the proportion of velocity error decreasing by less than 30% is about 18%, which is affected by random measurement error to a certain extent, the systematic errors caused by g-sensitive misalignments can be compensated for to a large extent.  All of the percentages of velocity error decrease corresponding with the 27 sets of vibration experiments are listed in Table 2. A corresponding statistic pie chart is shown in Figure 8. The statistical results show that the proportion of velocity error decreasing by more than 30% in the total results reach about 82% after g-sensitive misalignments compensation. As shown in Table 2 and Figure 8, though the proportion of velocity error decreasing by less than 30% is about 18%, which is affected by random measurement error to a certain extent, the systematic errors caused by g-sensitive misalignments can be compensated for to a large extent. time is for error analysis. The SINS works on pure inertial navigation mode during the vibration. The evaluation index is the percentage of velocity error decrease after g-sensitive misalignments compensation. Figure 7 depicts the velocity error curves in one of the random vibration tests. After g-sensitive misalignments compensation, the velocity error decreases by 66.1%. All of the percentages of velocity error decrease corresponding with the 27 sets of vibration experiments are listed in Table 2. A corresponding statistic pie chart is shown in Figure 8. The statistical results show that the proportion of velocity error decreasing by more than 30% in the total results reach about 82% after g-sensitive misalignments compensation. As shown in Table 2 and Figure 8, though the proportion of velocity error decreasing by less than 30% is about 18%, which is affected by random measurement error to a certain extent, the systematic errors caused by g-sensitive misalignments can be compensated for to a large extent.

Conclusions
In this paper, a novel method to calibrate g-sensitive misalignments of an RLG triad in linear vibration environments is proposed to improve the SINS accuracy for applications in spacecraft and other aerospace systems. Compared with the traditional calibration method in a gravitational field, the method presented in this paper is based on linear vibration experiments to simulate the harsh environment during aircraft flight. The derived mathematical model of g-sensitive misalignments can be used to relate attitude error to specific force and angular rate. With approximately optimal observations, g-sensitive misalignments are obtained by least squares estimations. Vibration experiments show the validity of the proposed method.