Next Article in Journal
Real-Time Vehicle Classification and Tracking Using a Transfer Learning-Improved Deep Learning Network
Next Article in Special Issue
A Wi-Fi Indoor Positioning Method Based on an Integration of EMDT and WKNN
Previous Article in Journal
sEMG and Vibration System Monitoring for Differential Diagnosis in Temporomandibular Joint Disorders
Previous Article in Special Issue
LiDAR- and Radar-Based Robust Vehicle Localization with Confidence Estimation of Matching Results
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Self-Calibration and Compensation Method for Installation Errors of Uniaxial Rotation Module Inertial Navigation System

1
College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou 310027, China
2
School of Aeronautics and Astronautics, Zhejiang University, Hangzhou 310027, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(10), 3812; https://doi.org/10.3390/s22103812
Submission received: 22 April 2022 / Revised: 14 May 2022 / Accepted: 16 May 2022 / Published: 17 May 2022
(This article belongs to the Collection Navigation Systems and Sensors)

Abstract

:
Calibration and compensation techniques are essential to improve the accuracy of the strap-down inertial navigation system. Especially for the new uniaxial rotation module inertial navigation system (URMINS), replacing faulty uniaxial rotation modules introduces installation errors between modules and reduces navigation accuracy. Therefore, it is necessary to calibrate these systems effectively and compensate for the installation error between modules. This paper proposes a new self-calibration and compensation method for installation errors without additional information and equipment. Using the attitude, velocity, and position differences between the two sets of navigation information output from URMINS as measurements, a Kalman filter is constructed and the installation error is estimated. After URMINS is compensated for the installation error, the average of the demodulated redundant information is taken to calculate the carrier’s navigation information. The simulation results show that the proposed method can effectively assess the installation error between modules with an estimation accuracy better than 5”. Experimental results for static navigation show that the accuracy of heading angle and positioning can be improved by 73.12% and 81.19% after the URMINS has compensated for the estimated installation errors. Simulation and experimental results further validate the effectiveness of the proposed self-calibration and compensation method.

1. Introduction

The superiority of a strap-down inertial navigation system (SINS) is that no external reference is needed to determine the vehicle’s attitude, velocity, or position, and no information is transmitted to the environment [1]. Therefore, it is extensively used in vehicles, submarines, missiles, aircraft and ships [2,3,4,5]. A global positioning system (GPS) is usually integrated with SINS to improve the carrier’s navigation and positioning accuracy. However, GPS is unavailable in some environments, such as underwater navigation and indoor localization [6]. Compared with the GPS, the SINS has the advantages of good concealment and strong anti-interference capability and can be applied in GPS-denied environments such as underwater navigation [7,8]. However, the drawback of SINS is that navigation errors accumulate with time owing to the drift-bias error of internal sensors such as gyroscope and accelerometer, which deteriorates the navigation accuracy of SINS [9]. Since the carrier may work in a covert environment, no other information is available to fuse with the SINS to suppress errors. Furthermore, if sensors in the SINS fail, these failed sensors must be replaced and maintained for the system to function correctly. However, reinstalling new sensors introduces errors into the system. The traditional method of calibrating the system for various errors is to take the inertial navigation system back to the laboratory for recalibration [10], which is inconvenient and wastes time. Therefore, the self-calibration and self-compensation techniques of errors in the SINS are of great significance for improving the performance of SINS [11,12].
Nowadays, rotation modulation technology is the primary method to realize the self-compensation of errors in the SINS, which can suppress the bias drift and other slowly changing errors of inertial sensors by periodic rotation [13,14]. The inertial sensors are rotated around one axis or multiple axes [15,16,17]. The uniaxial rotation modulated inertial navigation system (RMINS) has been widely studied and applied [18]. However, the inertial sensor’s constant bias along the rotational axis cannot be suppressed and still influences the navigation accuracy of SINS [19]. The inertial sensors’ constant errors are compensated in bi-axis and tri-axis RMINS, but the reliability of the RMINS is reduced due to the complex mechanics [20,21,22]. To avoid the complex multi-axis mechanism and enhance the accuracy and reliability of SINS, a novel triple rotary unit strap-down inertial navigation system (TRUSINS) was proposed. Each rotating unit of TRUSINS is uniaxial rotation, and the positioning accuracy of TRUSINS is improved by ten times compared to SINS [23]. It can be seen that the application of uniaxial rotation and data redundancy has a great development prospect in improving the accuracy of SINS. To improve the accuracy and configuration flexibility of SINS and simplify the design complexity of rotating shafting, a uniaxial rotation module inertial navigation system (URMINS) composed of several single-axis rotating modules is proposed in this paper. The structure of URMINS is described in Section 2 of this article. This navigation system can be applied in GPS-denied environments, such as outdoor navigation of vehicles, ships, and autonomous underwater vehicles.
However, with structures such as TRUSINS and URMINS, installation errors between modules have a detrimental effect on the system’s accuracy and must be compensated for [24]. The proposed methods used to calibrate errors of SINS are generally divided into direct and indirect methods. The direct calibration method uses the turntable to provide accurate reference information to SINS and estimates various error parameters in the system model by comparing the output of SINS with the reference information. For example, reference [25] used a single-axis turntable rotation to generate a reference angular rate for calibrating the gyro error. However, the calibration accuracy is related to the precision of the turntable [26]. The indirect calibration methods are used to determine error parameters by observing vehicle’s attitude, velocity or position errors. For example, reference [16] proposed an estimation method for calibrating nonorthogonal angle in dual-axis rotational INS, which takes attitude, velocity and altitude errors output by SINS as measurement values and assumes that the carrier is in a stationary state during the calibration process. Another study ([27]) based on the static environment used the attitude angle difference as the measurement value of a Kalman filter to estimate the error parameter to improve the attitude output accuracy. Once the carrier is in motion, the output of SINS contains the carrier’s velocity, attitude, and position information. It is difficult to estimate the error without the assistance of other sensors.
Current calibration methods are usually based on laboratory turntables or external information-generated reference signals to calibrate installation errors in the SINS [28]. In recent years, to calibrate the error of inertial sensors accurately, some novel calibration methods based on machine learning have also been widely studied [29,30]. However, these methods generally use high-precision data to train a network model and then apply the model to low-precision data to improve performance. That is to say, the novel calibration method based on machine learning still needs some external information assistance to obtain the error model of the system or sensor. Therefore, the limitation of the existing approach is that it is difficult to estimate the installation error caused by module replacement under dynamic conditions without external information. For calibrating the SINS’s error in a dynamic environment, a method was proposed to estimate the installation error by using the difference in navigation information output by two sets of SINS installed on the same carrier [31,32,33]. Nevertheless, the uniaxial rotation module (URM) consists of only two gyroscopes and two accelerometers, and a single module cannot output navigation information. Therefore, calibration methods for dynamic conditions cannot be directly applied to the URMINS to estimate installation errors caused by the replacement of a faulty module.
This paper proposes a new method for self-calibration and compensation of URMINS installation errors. Compared with existing methods, this proposed method can calibrate the installation error of URMINS in a dynamic environment without additional information sources or external devices. This approach uses the two sets of attitude, velocity, and position differences from URMINS as measurements. The Kalman filter is applied to estimate the installation error, given the linear URMINS error model and measurement error. The advantage of the Kalman filter is that it uses all measurements, regardless of their accuracy, to estimate system state values by properly weighting these measurements. However, we do not think that the Kalman filter is better than other filters when estimating installation errors. Many researchers have studied improvement algorithms of the traditional Kalman filter so that the performance of the Kalman filter has been improved [34].
The motivation behind and contributions of this paper are summarized as follows. First, to improve the navigation accuracy and system configuration flexibility of pure SINS and simplify the design of the rotating shaft, this paper proposes a URMINS composed of multiple URMs. Each URM can demodulate measurements from two gyroscopes and two accelerometers, with inertial sensor measurements oriented perpendicular to the axis of rotation. Second, to realize the replacement and maintenance of faulty URMs in a dynamic environment without relying on external equipment, this paper proposes a self-calibration and compensation method for the installation error caused by replacing faulty URMs. Using the difference between the two sets of navigation information output by URMINS as the measurement value of the Kalman filter, the installation error is stimulated by the movement of the carrier in the dynamic environment and is estimated. An indexing scheme is designed to suppress the inertial sensor’s constant and slowly changing errors. Installation errors are compensated for during the demodulation of the URMINS information. Both simulation and experimental results show that the method can calibrate the URMINS installation errors and that the navigation accuracy is improved after the installation error are compensated significantly.
The rest of this paper is organized as follows. In Section 2, the configuration of URMINS is presented, and the installation error model is established and analyzed; In Section 3, the self-calibration and compensation scheme are given. Simulation results are provided in Section 4 to prove the method’s validity. Experimental results and analysis are presented in Section 5. Finally, our conclusions are summarized in Section 6.

2. Configuration of URMINS and Installation Error Modeling

2.1. Configuration of URMINS

This system includes three same URMs and one navigation and power management module, and each module is mounted on the base. The mechanical structure of URMINS is shown in Figure 1. The rotary axes of the URM are perpendicular to each other, forming a Cartesian coordinate system. Each URM comprises two fiber optic gyroscopes (FOGs), two quartz-flexible accelerometers, a frameless torque motor, a RESM angle encoder, two read-heads, a wireless communication unit, and a wireless power transmission unit. The inertial sensors in each rotation framework rotate around the corresponding axis. The wireless data communication (wireless optical communication unit) and wireless power transmission are used in each module to avoid wire entanglement and friction of slip rings. The CAN bus is used for data transmission between modules and between each module and the navigation and power management module. Since each module has the same structure, the modules are interchangeable and convenient for replacement and maintenance.

2.2. The Installation Error Modeling and Analysis

Next, some proper coordinate frames relevant to self-calibration, navigation, and mutual transformations are discussed.
The i-frame is the earth-centered inertial frame. The Earth-centered Earth-fixed frame is the e-frame which rotates along with the earth. The navigation frame, named n-frame, its axes aligned with the east, north, and up directions. The body frame is referred to as the b-frame. The bi-frame ( i = 1 , 2 , 3 ) represents the orthogonal coordinate system obtained after calibration in the laboratory turntable for each module. Bi-frame ( i = 1 , 2 ) represents the equivalent body coordinate system of the combined strap-down inertial navigation system (CSINS), including CSINS1 and CSINS2. The s0-frame is the rotation coordinate frame under the ideal installation conditions at the initial moment. The s-frame is the rotation coordinate frame under the ideal installation conditions at any time. The g-frame and a-frame are the actual installation frames for the gyroscope and accelerometer, respectively. Here, the a-frame and g-frame are nonorthogonal coordinates, while the other frames are orthogonal. The coordinate frames mentioned above belong to the right-hand coordinate system.
To describe the transformation between two coordinate systems clearly, the transformation matrix C t k ( t = a , g , s , s 0 , b ;   k = s , s 0 , b , n , B i ) is introduced to represent the transformation from t-frame to k-frame. Since the internal structure of each module is the same, the direction of the rotation axis between modules is different. Therefore, take URM 1# as an example, and continue to analyze the installation error. In the b1-frame, the measured value of the accelerometer is f b 1 .
f b 1 = C s 0 b 1 C s s 0 C a s f a
And then, the transformation matrices in Equation (1) are written as follows.
C a s = [ 1 β x z 1 β x y 1 β y z 1 1 β y x 1 β z y 1 β z x 1 1 ]  
C s s 0 = [ c o s θ 1 s i n θ 1 0 s i n θ 1 c o s θ 1 0 0 0 1 ]  
C s 0 b 1 = [ 1 μ z 1 μ y 1 μ z 1 1 μ x 1 μ y 1 μ x 1 1 ]  
In the transformation matrix C a s , the installation error angles β z y 1 and β z x 1 are zero since no sensor is mounted along the Z a 1 -axis. The rotation angle of the motor is indicated by the symbol θ i   ( i = 1 , 2 , 3 ) , where i is the module number. The vector f a represents the measurement of the accelerometer installed in URM 1# in the a-frame, where there exists f a = [ f x 1 a   f y 1 a   0 ] T .
There are two types of installation errors in the URMINS, one is within the URM, and the other is between modules. The installation error inside the URM is introduced by the transformation from a nonorthogonal coordinate system, such as a-frame and g-frame, to an orthogonal coordinate system, such as b1-frame, b2-frame, and b3-frame. Take the transformation of the accelerometer’s frame as an example. The modeling of installation errors from the a-frame to the s-frame is shown in Figure 2a. Here, the installation error angle is represented by β i j k (i, j = x , y , z ; k = 1 ,   2 ,   3 ), representing the installation error angle obtained by rotating the i-axis of the a-frame around the j-axis of the s-frame. The installation error from the s0-frame to the b-frame is denoted by the symbol μ k i   ( i = 1 ,   2 ,   3 ;   k = x , y , z ), where i is the module number, and k represents the axis. Installation errors modeling of URM 1# from the s0-frame to b1-frame is shown in Figure 2b. When each module is calibrated on the same base, the b1-frame, the b2-frame, and the b3-frame coincide with the b-frame. Green cylinders represent sensors, including a gyroscope and an accelerometer. In the s0-frame, sensors are installed along the X-axis and the Y-axis. The Z-axis is the rotation axis without installing any sensors. Installation errors within the URM 2# and URM 3# are defined in the same way as URM 1#, and the rotation axes are the Y-axis and the X-axis, respectively.
The accelerometers’ output errors of the X-axis and Y-axis in the b-frame due to the absence of the Z-axis are shown in Equation (5). Here, we only consider the primary term of the installation error.
δ f x 1 b 1 = ( β y z 1 + μ z 1 ) s i n θ 1 f x 1 a ( β x z 1 + μ z 1 ) c o s θ 1 f y 1 a + ( μ y 1 + β y x 1 s i n θ 1 + β x y 1 c o s θ 1 ) f z a δ f y 1 b 1 = ( β y z 1 + μ z 1 ) c o s θ 1 f x 1 a ( β x z 1 + μ z 1 ) s i n θ 1 f y 1 a + ( μ x 1 + β x y 1 s i n θ 1 β y x 1 c o s θ 1 ) f z a
Similarly, the accelerometers’ output errors of the URM 2# and URM 3# can be deduced as follows. The derivation of Equations (6) and (7) is shown in Appendix A.
δ f x 2 b 2 = ( β z y 2 + μ y 2 ) s i n θ 2 f x 2 a + ( μ z 2 + β z x 2 s i n θ 2 β x z 2 c o s θ 2 ) f y a + ( β x y 2 + μ y 2 ) c o s θ 2 f z 2 a δ f z 2 b 2 = ( β z y 2 + μ y 2 ) c o s θ 2 f x 2 a + ( μ x 2 + β x z 2 s i n θ 2 + β z x 2 c o s θ 2 ) f y a ( β x y 2 + μ y 2 ) s i n θ 2 f z 2 a
δ f y 3 b 3 = ( μ z 3 + β z y 3 s i n θ 3 β y z 3 c o s θ 3 ) f x a ( β z x 3 + μ x 3 ) s i n θ 3 f y 3 a ( β y x 3 + μ x 3 ) c o s θ 3 f z 3 a δ f z 3 b 3 = ( μ y 3 + β y z 3 s i n θ 3 β z y 3 c o s θ 3 ) f x a + ( β z x 3 + μ x 3 ) c o s θ 3 f y 3 a ( β y x 3 + μ x 3 ) s i n θ 3 f z 3 a
For the URM 1#, installation errors β x y 1 , β y x 1 ,   μ x 1 ,   φ y 1 , β y z 1 + μ z 1 , β x z 1 + μ z 1 can be calibrated in the laboratory, and the non-orthogonal installation error β i j k will not be changed by module reinstallation, so this paper only considers the variation of orthogonal installation error μ k i between modules. Here, the variation of the installation error angle introduced by the reinstallation of the module is defined as Δ μ k i   ( i = 1 ,   2 ,   3 ;   k = x , y , z ) , and the installation error angle between b2-frame and b3-frame and b1-frame is shown in Figure 3a–b.
Since the installation error coupling terms β y z 1 + μ z 1 and β x z 1 + μ z 1 cannot be separated, they are treated as one term. The error term related to f z a in Equation (5) is not compensated in a single module since there is no sensor mounted in the direction of the rotation axis. Hence, the absence of this measurement f z a in a single module increases the navigation error. Similarly, the installation errors in URM 2# and URM 3# are the same as the analysis method in URM 1#. Therefore, it is necessary to study the self-calibration and compensation method of the installation error Δ μ k i , and a scheme will be given in Section 3.

3. Self-Calibration and Compensation Scheme

3.1. Self-Calibration Scheme

This system comprises three URMs, including six gyroscopes and six accelerometers. Two sets of inertial navigation systems, CSINS1 and CSINS2, can be formed using redundant information. There are four ways to configure two groups of inertial navigation systems when all sensors’ output information is applied. Using the difference in the navigation information of the two inertial navigation systems as a measure of the Kalman filter, the installation error caused by the module replacement is estimated.

3.1.1. Design of Kalman Filter

The angular velocity and specific force measured by the sensors perpendicular to the rotary axis are demodulated to the b-frame and output through the URM. Therefore, six angular velocities and specific forces from three URMs can be obtained in the b-frame, and there are two sets of measurements along with the directions X b , Y b , and Z b . Here, the CSINS1 consists of X b 1 , Y b 1 of URM 1#, and Z b 3 of URM 3#, while X b 2 , Z b 2 of URM 2#, and Y b 3 of URM 3# constitute the CSINS2. The two sets of CSINS calculate and output the navigation information. The difference in navigation information of CSINS1 and CSINS2 is used as the measurement value to estimate the installation errors between modules. Assuming that the b-frame overlapped with the b1-frame of URM 1#, the installation errors Δ μ x 1 , Δ μ y 1 and Δ μ z 1 are zero.
The state variable X is as follows.
X = [ ϕ E 12   ϕ N 12   ϕ U 12   δ V E 12   δ V N 12   δ L 12   δ λ 12   Δ μ x 2   Δ μ y 2   Δ μ z 2   Δ μ x 3   Δ μ y 3   Δ μ z 3 ] T  
The difference between the misalignment angles of CSINS1 and CSINS2 is written as ϕ 12 = [ ϕ E 12   ϕ N 12   ϕ U 12 ] T = ϕ 1 ϕ 2 . The difference between the east and north velocity errors of CSINS1 and CSINS2 is noted as δ V 12 = [ δ V E 12   δ V N 12 ] T = δ V 1 δ V 2 . The difference between the position errors of CSINS1 and CSINS2 is written as δ p 12 = [ δ L 12   δ λ 12 ] T = δ p 1 δ p 2 . The symbol μ k 2   ( k = x , y , z ) represents the installation error angle between URM 2# and URM 1#. The symbol μ k 3   ( k = x , y , z ) represents the installation error angle between URM 3# and URM 1#.
The attitude error equations are as follows.
ϕ ˙ E 12 = ( w U + V E t a n L R N + h ) ϕ N 12 ( w N + V E R N + h ) ϕ U 12 1 R M + h δ V N 12 δ w E 12 n
ϕ ˙ N 12 = ( w U + V E t a n L R N + h ) ϕ E 12 V N R M + h ϕ U 12 + 1 R N + h δ V E 12 w U δ L 12 δ w N 12 n
ϕ ˙ U 12 = ( w N + V E R N + h ) ϕ E 12 + V N R M + h ϕ N 12 + t a n L R N + h δ V E 12 + ( w N + V E s e c 2 L R N + h ) δ L 12 δ w U 12 n
The east and north velocity error equations are shown below.
δ V ˙ E 12 = f U ϕ N 12 + f N ϕ U 12 + V N t a n L R N + h δ V E 12 + ( 2 w U + V E t a n L R N + h ) δ V N 12 + ( 2 V N w N + V E V N s e c 2 L R N + h ) δ L 12 + δ f E 12 n
δ V ˙ N 12 = f U ϕ E 12 f E ϕ U 12 2 ( w U + V E t a n L R N + h ) δ V E 12 ( 2 V E w N + V E 2 s e c 2 L R N + h ) δ L 12 + δ f N 12 n
The positional error equations for longitude and latitude are as follows.
δ L ˙ 12 = 1 R M + h δ V N 12
δ λ ˙ 12 = s e c L R N + h δ V E 12 + V E t a n L s e c L R N + h δ L 12
The symbol w i e n is the earth’s rotation velocity projection under the n-frame, and its expression is w i e n = [ w E   w N   w U ] T . R M and R N are the curvature radius of the meridian and prime vertical, respectively. The symbol f i a n is the specific force measured by the accelerometer in the n-frame, and its expression is f i a n = C b n f i a b = [ f E   f N   f U ] T . The symbol δ w 12 n represents the difference in measurement error under the n-frame caused by gyroscopes of the two CSINSs, and its expression is δ w 12 n = [ δ w E 12 n   δ w N 12 n   δ w U 12 n ] T . The symbol δ f 12 n represents the difference in measurement error under the n-frame caused by the accelerometer measuring the assembly of the two CSINSs, and its expression is δ f 12 n = [ δ f E 12 n   δ f N 12 n   δ f U 12 n ] T . The specific expressions of the difference in sensor measurement error are shown in Equations (16) and (17).
δ w 12 n = C B 1 n δ w 1 B 1 C B 2 n δ w 2 B 2  
δ f 12 n = C B 1 n δ f 1 B 1 C B 2 n δ f 2 B 2  
Take the measurement error of the accelerometer as an example. It can be obtained from Equations (5)–(7). Furthermore, the measurement error of the accelerometer caused by CSINS1 and CSINS2 are presented as δ f 1 B 1 = [ δ f x 1 b   δ f y 1 b   δ f z 3 b ] T and δ f 2 B 2 = [ δ f x 2 b   δ f y 3 b   δ f z 2 b ] T . They are obtained from Equations (18) and (19). Here, there is δ f x 1 b = δ f y 1 b = 0 .
δ f x 2 b = μ y 2 s i n θ 2 f x 2 a μ z 2 f y a + μ y 2 c o s θ 2 f z 2 a δ f z 2 b = μ y 2 c o s θ 2 f x 2 a + μ x 2 f y a μ y 2 s i n θ 2 f z 2 a  
δ f y 3 b = μ z 3 f x a μ x 3 s i n θ 3 f y 3 a μ x 3 c o s θ 3 f z 3 a δ f z 3 b = μ y 3 f x a + μ x 3 c o s θ 3 f y 3 a μ x 3 s i n θ 3 f z 3 a  
In theory, there is C B 1 n = C B 2 n . The missing information in the third axis of each URM is determined by the average of the other two modules, whose expressions are shown in Equations (20)–(22). This method is also applicable to compensation installation errors of URMINS.
f x a = ( f x 1 b 1 + f x 2 b 2 ) / 2
f y a = ( f y 1 b 1 + f y 3 b 3 ) / 2
f z a = ( f z 2 b 2 + f z 3 b 3 ) / 2

3.1.2. Measurement Equation

By analyzing the structure of URMINS, it is obvious that the navigation information output by CSINS1 and CSINS2 comes from the same carrier. Therefore, there are:
A ˜ 1 = A + δ A 1 ,   A ˜ 2 = A + δ A 2 V ˜ 1 n = V n + δ V 1 n ,   V ˜ 2 n = V n + δ V 2 n p ˜ 1 = p + δ p 1 ,   p ˜ 2 = p + δ p 2
The symbols A ,     V n and p, respectively, represent the actual values of the carrier’s attitude, velocity, and position, and the specific expressions are shown in Equation (24).
A = [ θ   γ   ψ ] T , V n = [ V E n   V N n ] T , p = [ L   λ ] T  
The symbols A ˜ 1 , δ A 1 , A ˜ 2 and δ A 2 represent the attitude and attitude error of CSINS1 and CSINS2, respectively. The symbols V ˜ 1 n , δ V 1 n , V ˜ 2 n , δ V 2 n , p ˜ 1 , δ p 1 , p ˜ 2 , δ p 2 have the similar definition as the symbols of attitude.
The difference between the output information of CSINS1 and CSINS2 can be obtained by Equation (25). The acceleration and angular velocity caused by the carrier’s motion are not included in the measured values, so this method can estimate the installation error in a dynamic environment.
δ A ˜ 12 = A ˜ 1 A ˜ 2 ,   δ V ˜ 12 n = V ˜ 1 n V ˜ 2 n , δ p ˜ 12 = p ˜ 1 p ˜ 2  
The measurement equation of the system is as follows.
Z(t)=H(t)X(t)+V(t)
where Z(t) represents the measurement vector of the system. H(t) is the measurement matrix of the URMINS. V(t) is the measurement noise of the system. The detailed expressions of the symbol in Equation (26) are shown in Equations (27)–(30).
Z ( t ) = [ δ A ˜ 12   ;   δ V ˜ 12 n   ;   δ p ˜ 12 ] .
H ( t ) = [ C ϕ A 0 3 × 2 0 2 × 3 I 2 0 2 × 3 0 2 × 2     0 3 × 2 0 3 × 6 0 2 × 2 0 2 × 6 I 2 0 2 × 6 ]  
C ϕ A = [ c o s ψ s i n ψ 0 s i n ψ / c o s θ c o s ψ / c o s θ 0 s i n ψ t a n θ c o s ψ t n a θ 1 ]
I 2 = [ 1 0 0 1 ]
In this paper, to further demonstrate the effectiveness of the Kalman filter algorithm, observability analysis of the state variables in Kalman filter is essential. Here the attitude, velocity, and position errors ϕ 12 , δ V 12 , and δ p 12 are obtained from δ A ˜ 12 , δ V ˜ 12 n , and δ p ˜ 12 , respectively, so these navigation error terms are entirely observable. Analysis of Equations (18) and (19) shows that the installation error term Δ μ k i is related to the measurement value of the accelerometer, the encoder and the third-axis sensor of each module. Therefore, the linear motion, angular motion and design of the rotation scheme of the carrier can stimulate the installation error and improve the observability of the installation error of the system.

3.2. Compensation Scheme

For the installation error caused by replacing modules of URMINS, the block diagram of the installation error compensation method proposed in this paper is shown in Figure 4.
The rotation scheme of a URM needs to be designed, and the motor rotary speed of each module is set to 6°/s to achieve continuous rotation in one direction. Each module’s inertial sensor measurement data are demodulated, and the sensor data, called actual demodulation information, is obtained in the navigation solution module. Firstly, installation errors Δ μ y 2 and Δ μ x 3 are compensated to get the equivalent information f x a , f y a and f z a along the rotation axis of each module. Among them, the calculation method of the equivalent information is as Equations (20)–(22). Equations (31) and (32) show the compensated error terms. Secondly, the installation error between modules associated with URM 2# and URM 3# is compensated by Equations (18) and (19), and the sensor’s data under the b-frame after compensation is obtained. Finally, since the accuracy of the sensors of each module is the same, the average value is taken as the data output of each axis on the b-frame, and the traditional algorithm of SINS is used to calculate the navigation information.
δ f x 2 c = Δ μ y 2 s i n θ 2 f x 2 a + Δ μ y 2 c o s θ 2 f z 2 a δ f z 2 c = Δ μ y 2 c o s θ 2 f x 2 a Δ μ y 2 s i n θ 2 f z 2 a
δ f y 3 c = Δ μ x 3 s i n θ 3 f y 3 a Δ μ x 3 c o s θ 3 f z 3 a δ f z 3 c = Δ μ x 3 c o s θ 3 f y 3 a Δ μ x 3 s i n θ 3 f z 3 a  
We explain the error compensation method in detail. The error compensation process is further introduced by taking the error δ f x 2 b of the accelerometer of URM 2# as an example. The error terms presented by substituting URM 2# are shown in Equation (18). The algorithm in Section 2 can estimate the installation error μ k 2 , and the encoder measures the rotation angle θ 2 . f x 2 a and f z 2 a can be measured by two accelerometers mounted in the URM 2#. The unknown quantity f y a needs to be supplied from URM 1# or URM 3# as shown in Equation (34). Here, the installation error caused by Δ μ x 3 needs to be compensated first, as shown in Equation (33). Furthermore, the error caused by the replacement of URM 2# can be paid by Equation (34). The error compensation method of the remaining axis is similar to the above procedure and will not be described in detail.
f y 3 b = f y 3 b 3 δ f y 3 c
f y a = ( f y 1 b 1 + f y 3 b ) / 2

4. Simulation Results and Analysis

4.1. Simulation of Carrier Motion

Simulation experiments are implemented to verify the effectiveness of the installation error calibration method in Section 3.
The URMINS is mounted on moving carriers to effectively excite the calibration error to calibrate installation errors between modules. We consider the carrier’s linear and angular motion, the simulated carrier’s attitude, velocity, and trajectory change with time, as shown in Figure 5 and Figure 6. In Figure 6, the starting point of the carrier is at the origin of the coordinates.

4.2. Calibration Results of Installation Errors

The carrier motion in Section 4.1 is used as an incentive for self-calibration of the mounting error. The parameters used in the simulation, including inertial sensors’ errors and mounting errors, are shown in Table 1. The method to simulate inertial sensors’ data without noise was described in reference [35]. Here, inertial sensors’ errors are simulated with noise models such as Gaussian white noise, random constant deviation and first-order Markov process. The sensor data used in the simulation are obtained by summing the noise-free raw data with the noisy data. This paper assumes that the b1-frame of URM 1# overlaps with the b-frame, then the installation error terms μ x 1 ,   μ y 1 , and μ z 1 are set to zero. The initial alignment errors of CSINS1 and CSINS2 are the same, as shown in Table 2.
Assuming that the motor of each module is locked in a fixed position, the installation errors of URM 2# and URM 3# obtained by simulation are shown in Figure 7a,b, respectively. With the increase in simulation time, installation errors obtained by simulation present a stable trend. The covariance curves of installation error are shown in Figure 8. As shown in Figure 8, after the simulation time is 600 s, the covariance values of all installation errors are within 5″. Therefore, the installation error can be estimated quickly and accurately using the self-calibration method proposed in this paper.
Furthermore, the effectiveness of the proposed calibration method was verified, and the installation error was simulated five times. The mean value of installation errors was calculated after 600 s, and the mean value and residual of installation errors were presented in Table 3. The maximum value of the residuals of all installation errors is 1.54″, which further verifies the algorithm’s reliability.

4.3. Compensation Results and Analysis

The evaluation method, called true root mean square (TRMS), was used to evaluate the accuracy of the URMINS. The attitude, velocity, and position error curves before and after compensating for the installation error are shown in Figure 9a–c. The simulation data were produced based on MATLAB for one hour under the static condition of the carrier. The motor speed of each module was set to 6°/s. It can be seen from the simulation results of attitude angle errors in Figure 9a that, compared with before compensating for the installation error, the error of the misalignment angles is suppressed after compensating for the installation error. It can be seen from Figure 9b,c that after compensating for the installation error, both the speed errors and the position errors are significantly reduced. The positioning results before and after compensating for URMINS installation errors are shown in Figure 9d to demonstrate the advantages of positioning results after compensating for installation errors. After the installation error of URMINS is compensated, the positioning accuracy is higher.
To verify the validity of the installation error compensation method of the URMINS, the TRMS values of the five simulation results, including the attitude angle error, velocity error and position error, are listed in Table 4. It can be seen from Table 4 that the heading angle accuracy of the URMINS is improved by 2.61%, and the position accuracy of the URMINS is improved by 34.18%. In a word, the installation errors compensation method proposed in this paper can compensate for the installation errors of the URMINS without relying on external information in a dynamic environment. At the same time, the navigation accuracy of the URMINS can be improved.

5. Experimental Results and Analysis

5.1. Experimental Results of Self-Calibration

An experiment with the proposed method was performed using the URMINS fabricated by our laboratory to demonstrate the effectiveness of the self-calibration method in practical applications. In this paper, two modules are used to form URMINS to achieve self-calibration, and the experimental prototype is shown in Figure 10. A computer collected the data output by the URMINS at 100 Hz. The self-calibration method proposed in this paper does not depend on external references, such as turntable and GPS. It only needs the carrier to be dynamic to stimulate the installation error. To verify the effectiveness of the proposed self-calibration algorithm, we installed the URMINS on a dual-axis turntable. The turntable was controlled to provide a sway randomly motion for the URMINS, simulating the angular motion of the carrier. The installation error is estimated by applying a self-calibration algorithm. The two-axis turntable used in the experiment is 2TD-550 fabricated by the Jiujiang Precision Test Technology Research Institute, and the specific parameters of this turntable are shown in Table 5.
Here, the CSINS1 consists of X b 2 , Z b 2 of URM 2#, and Y b 3 of URM 3#, while X b 2 , Z b 3 of URM 3#, and X b 2 of URM 2# constitute the CSINS2. The structure of each module is consistent with that described in Section 2.1. The difference in navigation information of CSINS1 and CSINS2 was used as the measurement value to estimate the installation errors between modules. Assuming that the b-frame overlapped with the b2-frame of URM 2#, the installation error Δ μ k 2 is zero. Reinstall URM 3# and use the dynamic self-calibration method proposed in this paper to estimate the installation error Δ μ k 3 . The motor of each module was controlled to stop at a position where the angle information output by the encoder is zero. The turntable is controlled to sway randomly to simulate the random motion of the carrier. The attitude of the carrier is shown in Figure 11.
The data along the Z b -axis of the carrier were obtained through URM 2# and URM 3# by this experimental prototype, and only the information along the Z b -axis was redundant. These data along the X b -axis and Y b -axis were reused in CSINS1 and CSINS2. Therefore, by further analyzing the error term δ f z 3 b of Equation (19), it can be concluded that installation errors Δ μ x 3 and Δ μ y 3 are fully observable, while Δ μ z 3 is unobservable.
The experimental results of the installation error for URM 3# are shown in Figure 12. Δ μ x 3 and Δ μ y 3 tend to be stable over time, while Δ μ z 3 cannot be estimated. The covariance curve of the installation error is shown in Figure 13, which also proves that Δ μ z 3 is unobservable, and Δ μ x 3 and Δ μ y 3 are gradually converging over time. The mean value of the installation error was calculated after 800 s, and the mean values of Δ μ x 3 and Δ μ y 3 were −544.72” and −871.52”, respectively. Therefore, the experimental results agree with the theoretical analysis, which proves the effectiveness of the proposed self-calibration method.

5.2. Compensation Results and Analysis

Compensate for the installation error of URMINS with the estimated mean values of Δ μ x 3 and Δ μ y 3 . Since Δ μ z 3 cannot be calculated, this value is set to zero. The experimental data were collected for one hour under the static condition of the carrier. The motor’s speed of each module was set to 5.32°/s. From the experimental results of the navigation errors in Figure 14a–c, it is found that the attitude angle error, velocity error and position error are suppressed after compensating for the installation error compared to before compensating for the installation error. The positioning results before and after compensating for the URMINS installation error are shown in Figure 14d to further demonstrate the advantages of the positioning results after compensating for the installation error. The positioning accuracy is higher after compensating for the installation error of URMINS.
To digitally present the improvement effect of navigation accuracy, the experimental results, including the TRMS values of attitude angle error, velocity error, and position error, are statistics in Table 6. Among the navigation error parameters, the accuracy improvement range of the eastward misalignment angle is the lowest at 53.51%, and the suppression ranges of the other navigation error parameters are higher than this value. It can be seen from Table 6 the heading angle accuracy of the URMINS was enhanced by 73.12%, and the position accuracy of the URMINS was improved by 81.19%. In conclusion, the experimental results demonstrate that the installation error compensation method proposed in this paper can compensate for the installation error of URMINS without relying on external information and improve the navigation accuracy.

6. Conclusions

This paper presents a new self-calibration and compensation method for URMINS installation error, which can estimate the installation error caused by the replacement of a faulty module without additional equipment and significantly improve the navigation accuracy of the system. In this paper, the installation error model of URMINS was established and analyzed. Then, three URMs were combined with URMINS, and the system was designed to output two sets of navigation information. The difference in attitude, speed, and position between two groups of CSINS was taken as the measurement value of the Kalman filter, and the installation error between modules was estimated. After demodulation, the average value of redundant information was taken to calculate the carrier’s attitude, velocity, and position. The simulation and experimental results show that the installation errors caused by module replacement can tend to the exact value with this method. Compared with the static navigation results without compensating for the installation error, both simulation and experiment prove that the carrier’s attitude, speed, and position accuracy significantly improve after paying for the installation error. The experimental results show that after compensating for the installation error, the heading angle accuracy and horizontal positioning accuracy of URMINS are improved by 73.12% and 81.19%, respectively. The self-calibration method proposed in this paper can be applied in a real situation where the carrier is in a dynamic environment and the inertial sensor has information redundancy. The estimation of the installation error can be realized without relying on the external reference. In addition, the URMINS is easy to maintain and can suppress errors in pure INSs, which provides convenience for carriers to achieve long-term, high-precision navigation in concealed environments.

Author Contributions

Conceptualization, M.N. and X.S.; methodology, M.N. and K.S.; software, H.M.; validation, M.N. and H.M.; formal analysis, T.H. and K.S.; investigation, M.N.; data curation, H.M.; writing—original draft preparation, M.N.; writing—review and editing, X.S. and T.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Aerospace Science and Technology Fund of China (grant number KM20200075). The authors would appreciate the support and funds.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank Shangru Duan of the College of Biomedical Engineering and Instrument Science of Zhejiang University for assistance with simulation verifications related to this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The accelerometers’ output errors of the X-axis and Z-axis in the b-frame due to the absence of the Y-axis are shown in Equation (6). Here, we only consider the primary term of the installation error. In the b2-frame, the measured value of the accelerometer is f b 2 .
f b 2 = C s 20 b 2 C s 2 s 20 C a 2 s 2 f a 2
And then, the transformation matrices in Equation (A1) are written as follows.
C a 2 s 2 = [ 1 β x z 2 β x y 2 β y z 2 1 β y x 2 β z y 2 β z x 2 1 ]
C s 2 s 20 = [ c o s θ 2 0 s i n θ 2 0 1 0 s i n θ 2 0 c o s θ 2 ]
C s 20 b 2 = [ 1 μ z 2 μ y 2 μ z 2 1 μ x 2 μ y 2 μ x 2 1 ]
In transformation matrix C a 2 s 2 , the installation error angles β y z 2 and β y x 2 are zero since no sensor is mounted along the Y a 2 2 -axis. The rotation angle of the motor is indicated by the symbol θ 2 . The vector f a 2 represents the measurement of the accelerometer installed in URM 2# in the a2-frame, where there exists f a 2 = [ f x 2 a   0   f z 2 a ] T .
The accelerometers’ output errors of the Y-axis and Z-axis in the b-frame due to the absence of the X-axis are shown in Equation (7). Here, we only consider the primary term of the installation error. In the b3-frame, the measured value of the accelerometer is f b 3 .
f b 3 = C s 30 b 3 C s 3 s 30 C a 3 s 3 f a 3
And then, the transformation matrices in Equation (A5) are written as follows.
C a 3 s 3 = [ 1 β x z 3 β x y 3 β y z 3 1 β y x 3 β z y 3 β z x 3 1 ]
C s 3 s 30 = [ 1 0 0 0 c o s θ 3 s i n θ 3 0 s i n θ 3 c o s θ 3 ]
C s 30 b 3 = [ 1 μ z 3 μ y 3 μ z 3 1 μ x 3 μ y 3 μ x 3 1 ]
In transformation matrix C a 3 s 3 , the installation error angles β x z 3 and β x y 3 are zero since no sensor is mounted along the X a 3 3 -axis. The rotation angle of the motor is indicated by the symbol θ 3 . The vector f a 3 represents the measurement of the accelerometer installed in URM 3# in the a3-frame, where there exists f a 3 = [ 0   f y 3 a   f z 3 a ] T .

References

  1. Noureldin, A.; Karamat, T.B.; Georgy, J. Fundamentals of Inertial Navigation, Satellite-Based Positioning and Their Integration; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; ISBN 978-3-642-30465-1. [Google Scholar]
  2. Adams, G.; Gokhale, M. Fiber Optic Gyro Based Precision Navigation for Submarines. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Boston, MA, USA, 19–22 August 2013; p. 4384. [Google Scholar]
  3. Kang, L.; Ye, L.; Song, K. A Fast In-Motion Alignment Algorithm for DVL Aided SINS. Math. Probl. Eng. 2014, 2014, e593692. [Google Scholar] [CrossRef] [Green Version]
  4. Huang, T.; Jiang, H.; Zou, Z.; Ye, L.; Song, K. An Integrated Adaptive Kalman Filter for High-Speed UAVs. Appl. Sci. 2019, 9, 1916. [Google Scholar] [CrossRef] [Green Version]
  5. Zou, Z.; Huang, T.; Ye, L.; Song, K. CNN Based Adaptive Kalman Filter in High-Dynamic Condition for Low-Cost Navigation System on Highspeed UAV. In Proceedings of the 2020 5th Asia-Pacific Conference on Intelligent Robot Systems (ACIRS), Singapore, 17–19 July 2020; pp. 103–108. [Google Scholar]
  6. Shan, G.; Park, B.; Nam, S.; Kim, B.; Roh, B.; Ko, Y.-B. A 3-Dimensional Triangulation Scheme to Improve the Accuracy of Indoor Localization for IoT Services. In Proceedings of the 2015 IEEE Pacific Rim Conference on Communications, Computers and Signal Processing (PACRIM), Victoria, BC, Canada, 24–26 August 2015; pp. 359–363. [Google Scholar]
  7. Liu, Y.; Fan, X.; Lv, C.; Wu, J.; Li, L.; Ding, D. An Innovative Information Fusion Method with Adaptive Kalman Filter for Integrated INS/GPS Navigation of Autonomous Vehicles. Mech. Syst. Signal Process. 2018, 100, 605–616. [Google Scholar] [CrossRef] [Green Version]
  8. Xing, H.; Liu, Y.; Guo, S.; Shi, L.; Hou, X.; Liu, W.; Zhao, Y. A Multi-Sensor Fusion Self-Localization System of a Miniature Underwater Robot in Structured and GPS-Denied Environments. IEEE Sens. J. 2021, 21, 27136–27146. [Google Scholar] [CrossRef]
  9. Tucker, T.; Levinson, E. The AN/WSN-7B Marine Gyrocompass/Navigator. In Proceedings of the 2000 National Technical Meeting of the Institute of Navigation, Anaheim, CA, USA, 26–28 January 2000; pp. 348–357. [Google Scholar]
  10. Bekkeng, J.K. Calibration of a Novel MEMS Inertial Reference Unit. IEEE Trans. Instrum. Meas. 2009, 58, 1967–1974. [Google Scholar] [CrossRef]
  11. Levinson, E.; Majure, R. In Proceedings of the Accuracy Enhancement Techniques Applied to the Marine Ring Laser Inertial Navigator (MARLIN), Anaheim, CA, USA, 20–23 January 1987; pp. 71–80.
  12. Ban, J.; Chen, G.; Meng, Y.; Shu, J. Calibration Method for Misalignment Angles of a Fiber Optic Gyroscope in Single-Axis Rotational Inertial Navigation Systems. Opt. Express 2022, 30, 6487–6499. [Google Scholar] [CrossRef]
  13. Ishibashi, S.; Tsukioka, S.; Sawa, T.; Yoshida, H.; Hyakudome, T.; Tahara, J.; Aoki, T.; Ishikawa, A. The Rotation Control System to Improve the Accuracy of an Inertial Navigation System Installed in an Autonomous Underwater Vehicle. In Proceedings of the 2007 Symposium on Underwater Technology and Workshop on Scientific Use of Submarine Cables and Related Technologies, Tokyo, Japan, 17–20 April 2007; pp. 495–498. [Google Scholar]
  14. Levinson, E.; Ter Horst, J.; Willcocks, M. The Next Generation Marine Inertial Navigator Is Here Now. In Proceedings of the 1994 IEEE Position, Location and Navigation Symposium-PLANS’94, Las Vegas, NV, USA, 11–15 April 1994; pp. 121–127. [Google Scholar]
  15. Li, K.; Chen, Y.; Wang, L. Online Self-Calibration Research of Single-Axis Rotational Inertial Navigation System. Measurement 2018, 129, 633–641. [Google Scholar] [CrossRef]
  16. Deng, Z.; Sun, M.; Wang, B.; Fu, M. Analysis and Calibration of the Nonorthogonal Angle in Dual-Axis Rotational INS. IEEE Trans. Ind. Electron. 2017, 64, 4762–4771. [Google Scholar] [CrossRef]
  17. Gao, P.; Li, K.; Song, T.; Liu, Z. An Accelerometers-Size-Effect Self-Calibration Method for Triaxis Rotational Inertial Navigation System. IEEE Trans. Ind. Electron. 2018, 65, 1655–1664. [Google Scholar] [CrossRef]
  18. Sun, W.; Xu, A.-G.; Che, L.-N.; Gao, Y. Accuracy Improvement of SINS Based on IMU Rotational Motion. IEEE Aerosp. Electron. Syst. Mag. 2012, 27, 4–10. [Google Scholar] [CrossRef]
  19. Li, K.; Gao, P.; Wang, L.; Zhang, Q. Analysis and Improvement of Attitude Output Accuracy in Rotation Inertial Navigation System. Math. Probl. Eng. 2015, 2015, e768174. [Google Scholar] [CrossRef] [Green Version]
  20. Kang, L.; Ye, L.; Song, K.; Zhou, Y. Attitude Heading Reference System Using MEMS Inertial Sensors with Dual-Axis Rotation. Sensors 2014, 14, 18075–18095. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Yuan, B.; Liao, D.; Han, S. Error Compensation of an Optical Gyro INS by Multi-Axis Rotation. Meas. Sci. Technol. 2012, 23, 025102. [Google Scholar] [CrossRef]
  22. Hu, X.; Wang, Z.; Weng, H.; Zhao, X. Self-Calibration of Tri-Axis Rotational Inertial Navigation System Based on Virtual Platform. IEEE Trans. Instrum. Meas. 2021, 70, 1–10. [Google Scholar] [CrossRef]
  23. Chen, G.; Li, K.; Wang, W.; Li, P. A Novel Redundant INS Based on Triple Rotary Inertial Measurement Units. Meas. Sci. Technol. 2016, 27, 105102. [Google Scholar] [CrossRef]
  24. Liu, J.; Chen, G. A Four-Position Calibration Method of Inconsistent Angles between Units of TRUSINS. In Proceedings of the 2019 International Conference on Communications, Information System and Computer Engineering (CISCE), Haikou, China, 5–7 July 2019; pp. 529–532. [Google Scholar]
  25. Syed, Z.F.; Aggarwal, P.; Goodall, C.; Niu, X.; El-Sheimy, N. A New Multi-Position Calibration Method for MEMS Inertial Navigation Systems. Meas. Sci. Technol. 2007, 18, 1897–1907. [Google Scholar] [CrossRef]
  26. Song, N.; Cai, Q.; Yang, G.; Yin, H. Analysis and Calibration of the Mounting Errors between Inertial Measurement Unit and Turntable in Dual-Axis Rotational Inertial Navigation System. Meas. Sci. Technol. 2013, 24, 115002. [Google Scholar] [CrossRef]
  27. Jiang, R.; Yang, G.; Zou, R.; Wang, J.; Li, J. Accurate Compensation of Attitude Angle Error in a Dual-Axis Rotation Inertial Navigation System. Sensors 2017, 17, 615. [Google Scholar] [CrossRef] [Green Version]
  28. Hwang, Y.; Jeong, Y.; Kweon, I.S.; Choi, S. Online Misalignment Estimation of Strapdown Navigation for Land Vehicle under Dynamic Condition. Int. J. Automot. Technol. 2021, 22, 1723–1733. [Google Scholar] [CrossRef]
  29. Peng, H.; Bai, X. Machine Learning Gyro Calibration Method Based on Attitude Estimation. J. Spacecr. Rockets 2021, 58, 1636–1651. [Google Scholar] [CrossRef]
  30. Mahdi, A.E.; Azouz, A.; Abdalla, A.E.; Abosekeen, A. A Machine Learning Approach for an Improved Inertial Navigation System Solution. Sensors 2022, 22, 1687. [Google Scholar] [CrossRef] [PubMed]
  31. Hua, M.; Li, K.; Lv, Y.; Wu, Q. A Dynamic Calibration Method of Installation Misalignment Angles between Two Inertial Navigation Systems. Sensors 2018, 18, 2947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hu, P.; Xu, P.; Chen, B.; Wu, Q. A Self-Calibration Method for the Installation Errors of Rotation Axes Based on the Asynchronous Rotation of Rotational Inertial Navigation Systems. IEEE Trans. Ind. Electron. 2018, 65, 3550–3558. [Google Scholar] [CrossRef]
  33. Sun, C.; Li, K. A Positioning Accuracy Improvement Method by Couple RINSs Information Fusion. IEEE Sens. J. 2021, 21, 19351–19361. [Google Scholar] [CrossRef]
  34. Shen, C.; Zhang, Y.; Tang, J.; Cao, H.; Liu, J. Dual-Optimization for a MEMS-INS/GPS System during GPS Outages Based on the Cubature Kalman Filter and Neural Networks. Mech. Syst. Signal Process. 2019, 133, 106222. [Google Scholar] [CrossRef]
  35. Yan, G.; Wang, J.; Zhou, X. High-Precision Simulator for Strapdown Inertial Navigation Systems Based on Real Dynamics from GNSS and IMU Integration. In Proceedings of the China Satellite Navigation Conference (CSNC), Xi’an, China, 13–15 May 2015; Sun, J., Liu, J., Fan, S., Lu, X., Eds.; Springer: Berlin/Heidelberg, Germany, 2015; Volume 3, pp. 789–799. [Google Scholar]
Figure 1. The mechanical structure of URMINS.
Figure 1. The mechanical structure of URMINS.
Sensors 22 03812 g001
Figure 2. (a) Installation errors from the a-frame to the s-frame; (b) the installation errors of URM 1# from the s0-frame to the b1-frame.
Figure 2. (a) Installation errors from the a-frame to the s-frame; (b) the installation errors of URM 1# from the s0-frame to the b1-frame.
Sensors 22 03812 g002
Figure 3. (a) The installation error of URM 2# from the b2-frame to the b1-frame; (b) the installation error of URM 3# from the b3-frame to the b1-frame.
Figure 3. (a) The installation error of URM 2# from the b2-frame to the b1-frame; (b) the installation error of URM 3# from the b3-frame to the b1-frame.
Sensors 22 03812 g003
Figure 4. A block diagram of compensating for installation errors.
Figure 4. A block diagram of compensating for installation errors.
Sensors 22 03812 g004
Figure 5. (a) The attitude of the carrier; (b) the velocity of the carrier.
Figure 5. (a) The attitude of the carrier; (b) the velocity of the carrier.
Sensors 22 03812 g005
Figure 6. The trajectory of the carrier.
Figure 6. The trajectory of the carrier.
Sensors 22 03812 g006
Figure 7. (a) Simulation results of installation errors for URM 2#; (b) simulation results of installation errors for URM 3#.
Figure 7. (a) Simulation results of installation errors for URM 2#; (b) simulation results of installation errors for URM 3#.
Sensors 22 03812 g007
Figure 8. The covariance of installation errors is obtained by simulation.
Figure 8. The covariance of installation errors is obtained by simulation.
Sensors 22 03812 g008
Figure 9. (a) Simulation results of attitude errors; (b) simulation results of velocity errors; (c) simulation results of position errors; (d) compare the trajectory of the carrier.
Figure 9. (a) Simulation results of attitude errors; (b) simulation results of velocity errors; (c) simulation results of position errors; (d) compare the trajectory of the carrier.
Sensors 22 03812 g009
Figure 10. A prototype is used in the experiment.
Figure 10. A prototype is used in the experiment.
Sensors 22 03812 g010
Figure 11. The attitude of the carrier was simulated by the sway of the turntable.
Figure 11. The attitude of the carrier was simulated by the sway of the turntable.
Sensors 22 03812 g011
Figure 12. Experimental results of installation errors for URM 3#.
Figure 12. Experimental results of installation errors for URM 3#.
Sensors 22 03812 g012
Figure 13. Experimental results of installation errors’ covariance.
Figure 13. Experimental results of installation errors’ covariance.
Sensors 22 03812 g013
Figure 14. (a) Experimental results of attitude errors; (b) experimental results of velocity errors; (c) experimental results of position errors; (d) compare the trajectory of the carrier.
Figure 14. (a) Experimental results of attitude errors; (b) experimental results of velocity errors; (c) experimental results of position errors; (d) compare the trajectory of the carrier.
Sensors 22 03812 g014
Table 1. Set the errors used in the simulation.
Table 1. Set the errors used in the simulation.
Inertial Sensors ErrorsValueInertial Sensors ErrorsValue
Gyroscope constant drift 0.02   ° / h Accelerometer constant bias 50   μ g
Gyroscope scale factor 50   ppm Accelerometer scale factor 50   ppm
Gyroscope angle random walk 0.02 ° / h Accelerometer random walk 5   μ g / Hz
URM 2# orthogonal installation error μ k 2 360 URM 3# orthogonal installation error μ k 3 360
Table 2. Initial alignment errors.
Table 2. Initial alignment errors.
ErrorsValueErrorsValue
Pitch error (″)20Velocity error (m/s)0.01
Roll error (″)20Position error (m)0.1
Yaw error (°)1
Table 3. The installation errors and residuals were obtained by simulation.
Table 3. The installation errors and residuals were obtained by simulation.
μ x 2 ( ) μ y 2 ( ) μ z 2 ( ) μ x 3 ( ) μ y 3 ( ) μ z 3 ( )
1361.86359.19359.08359.73357.42357.62
2360.03359.90358.70360.67359.71359.40
3358.21360.05361.03359.93362.16358.63
4367.11360.27359.82361.17353.90359.95
5359.88359.51359.23360.20359.13359.45
Average361.42359.79359.57360.34358.46359.01
Residual−1.420.210.43−0.341.540.99
Table 4. Simulation results of attitude angle error, velocity error, and position error.
Table 4. Simulation results of attitude angle error, velocity error, and position error.
Navigational
Errors
Before
Compensation
After
Compensation
Improvement Range (%)
ϕ E   ( ° ) 0.070.0448.60
ϕ N   ( ° ) 0.060.0185.29
ϕ U   ( ° ) 1.000.972.61
δ V E   ( m / s ) 4.081.6260.35
δ V N   ( m / s ) 11.159.8711.51
δ L   ( km ) 21.0413.9733.60
δ λ   ( km ) 5.542.8947.78
δ P   ( km ) 21.5914.2134.18
Table 5. Parameters of the two-axis turntable.
Table 5. Parameters of the two-axis turntable.
ParametersValueParametersValue
Angular position error≤2″Angular rate error≤0.002°/s
Rotational error of shafting≤2″Angular rate range of the inner frame ± 0.001 ~ ± 500 ° / s
Shafting perpendicularity error ≤2″Angular rate range of the outer frame ± 0.001 ~ ± 300 ° / s
Table 6. Experimental results of attitude angle error, velocity error, and position error.
Table 6. Experimental results of attitude angle error, velocity error, and position error.
Navigational
Errors
Before
Compensation
After
Compensation
Improvement Range (%)
ϕ E   ( ° ) 0.340.1653.51
ϕ N   ( ° ) 0.860.1484.28
ϕ U   ( ° ) 1.000.2773.12
δ V E   ( m / s ) 65.788.8186.60
δ V N   ( m / s ) 87.6013.4884.61
δ L   ( km ) 154.5426.1583.08
δ λ   ( km ) 57.7318.1268.62
δ P   ( km ) 162.2530.5281.19
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Niu, M.; Ma, H.; Sun, X.; Huang, T.; Song, K. A New Self-Calibration and Compensation Method for Installation Errors of Uniaxial Rotation Module Inertial Navigation System. Sensors 2022, 22, 3812. https://doi.org/10.3390/s22103812

AMA Style

Niu M, Ma H, Sun X, Huang T, Song K. A New Self-Calibration and Compensation Method for Installation Errors of Uniaxial Rotation Module Inertial Navigation System. Sensors. 2022; 22(10):3812. https://doi.org/10.3390/s22103812

Chicago/Turabian Style

Niu, Meng, Hongyu Ma, Xinglin Sun, Tiantian Huang, and Kaichen Song. 2022. "A New Self-Calibration and Compensation Method for Installation Errors of Uniaxial Rotation Module Inertial Navigation System" Sensors 22, no. 10: 3812. https://doi.org/10.3390/s22103812

APA Style

Niu, M., Ma, H., Sun, X., Huang, T., & Song, K. (2022). A New Self-Calibration and Compensation Method for Installation Errors of Uniaxial Rotation Module Inertial Navigation System. Sensors, 22(10), 3812. https://doi.org/10.3390/s22103812

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