Next Article in Journal
Estimating Compressive and Shear Forces at L5-S1: Exploring the Effects of Load Weight, Asymmetry, and Height Using Optical and Inertial Motion Capture Systems
Previous Article in Journal
Research on a Cross-Domain Few-Shot Adaptive Classification Algorithm Based on Knowledge Distillation Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of a Multi-Position Alignment Scheme

1
School of Automation, Northwestern Polytechnical University, Xi’an 710072, China
2
Shaanxi Aircraft Corporation, Xi’an 723214, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(6), 1938; https://doi.org/10.3390/s24061938
Submission received: 4 February 2024 / Revised: 10 March 2024 / Accepted: 14 March 2024 / Published: 18 March 2024
(This article belongs to the Special Issue Challenges and Future Trends of Inertial Sensors)

Abstract

:
The current new type of inertial navigation system, including rotating inertial navigation systems and three-autonomy inertial navigation systems, has been increasingly widely applied. Benefited by the rotating mechanisms of these inertial navigation systems, alignment accuracy can be significantly enhanced by implementing IMU (Inertial Measurement Unit) rotation during the alignment process. The principle of suppressing initial alignment errors using rotational modulation technology was investigated, and the impact of various component error terms on alignment accuracy of IMU during rotation was analyzed. A corresponding error suppression scheme was designed to overcome the shortcoming of the significant scale factor error of fiber optic gyroscopes, and the research content of this paper is validated through corresponding simulations and experiments. The results indicate that the designed alignment scheme can effectively suppress the gyro scale factor error introduced by angular motion and improve alignment accuracy.

1. Introduction

Initial alignment is one of the key technologies in inertial navigation [1,2], and the inertial navigation system (INS) needs to perform initial alignment before entering the navigation state during every startup. Due to the fact that inertial navigation calculation is an integral process, and initial alignment is actually a problem of determining the initial value of integration [3], the alignment result will directly affect the navigation accuracy of the system [4]. Therefore, INS requires accuracy for initial alignment. It requires reducing the initial alignment error to introduce irreversible effects on navigation solutions.
At present, with the application of new INSs such as rotating INS [5,6,7,8] and three-autonomy (autonomous calibration, autonomous alignment and autonomous test) INS [9], some new alignment schemes have also emerged. Due to the fact that the above-mentioned INS has a rotating mechanism, schemes such as multi-position alignment and rotation alignment have been proposed [10,11].
The multi-position alignment scheme utilizes the error estimation and compensation principle of IMU to improve alignment accuracy [12]. The two-position alignment scheme, which rotates 180° around the azimuth axis, is a widely applied multi-position alignment scheme. Fang et al. [13] addressed the issue of poor accuracy in high-precision airborne POS ground two-position alignment caused by output noise from inertial devices by using an adaptive Kalman filter to improve alignment accuracy. Lai et al. [14] proposed an arithmetic model of a fast two-position initial alignment for INS using unscented Kalman filter. Zhang et al. [15] used a piecewise combined Kalman filter utilizing relative azimuth constraint to improve the alignment precision and to reduce the time consumption for error convergence.
Compared with multi-position alignment, rotation alignment is another approach, which uses the principle of error modulation to suppress the impact of IMU-related errors on alignment accuracy through the periodic rotation of the IMU. Wang and Li [16] used a rotation operation and calculated the error based on the projection relationship of the alignment steady-state error at different positions. By compensating for the error, the alignment accuracy of the system is improved. Wei et al. [17] designed a new initial alignment algorithm for dual-axis rotation alignment to solve the problem of azimuth error divergence caused by traditional alignment algorithms when IMU rotates around the horizontal axis.
Because the additional angular motion is introduced by the above alignment methods, analyses of the impact of inertial component errors on the system under rotational conditions term by term are necessary. In this study, detailed analyses of various gyroscope errors under angular motion were conducted, and the impact of these errors on alignment accuracy was assessed. Owing to the coupling of rotation and gyro scale factor error during the alignment process, new additional errors are introduced. Under the condition of large-scale factor errors in gyroscopes, existing two-position alignment schemes will produce significant errors that cannot be ignored. A new multi-position alignment scheme was designed based on the two-position alignment scheme, which effectively suppresses the impact of gyro scale factor error.
The remainder of this paper is structured as follows: The Section 2 elucidates the coordinate definitions employed in the analysis and delineates the error modes inherent to INS. Subsequently, the Section 3 concentrates on scrutinizing the influence of gyroscope errors on the alignment process under rotational circumstances, and devises a multi-position alignment methodology predicated on these analyses. This is followed by the Section 4, wherein the efficacy of the designed multi-position alignment scheme is empirically validated through mathematical simulations and experimental investigations. Finally, the conclusions drawn from this study are encapsulated in the Section 5.

2. Coordinate System and Error Model of INS

2.1. Coordinate System Definition

In the strapdown inertial navigation system, the commonly used coordinate systems include inertial frame, earth frame, geographic frame, navigation frame, body frame, etc. In the rotation inertial navigation system, another IMU frame is added [18,19]. The coordinate system used in the subsequent analysis of this paper is defined as follows:
Inertial frame (i frame)— O e x i y i z i : The origin is located in the center of the earth, and the coordinate axis has no rotation relative to the star. Axis O e z i is traditionally selected along the north pole of the earth. Axis O e x i and axis O e y i are selected in the equatorial plane of the earth and point to two stars in space.
Earth frame (e frame)— O e x e y e z e : The origin is at the center of the earth, and the coordinate axis is firmly connected with the earth. Axis O e z e points to the direction of the earth’s polar axis. Axis O e x e is along the intersection of the Greenwich meridian and the equatorial plane of the earth. The earth frame rotates around axis O e z i at angular velocity ω i e relative to the inertial frame.
Geographic frame (g frame)— O x g y g z g : The coordinate origin is located at the center of gravity of the carrier. Axes o x g and o y g are at the local level, where axis o x g points east and axis o y g points north. Axis o z g points to the sky along the vertical direction.
Navigation frame (n frame)— o x n y n z n : The navigation frame is selected as the navigation reference according to the needs of the system during navigation. This paper uses the geographic frame as the navigation frame.
Body frame (b frame)— o x b y b z b : The body frame is fixed on the carrier. The coordinate origin o of the body frame is located at the center of gravity of the carrier. Axis o x b points to the right along the horizontal axis of the carrier, and axis o y b points forward along the longitudinal axis of the carrier. Axes o z b , o x b and o y b form a right-handed coordinate system.
IMU frame (p frame)— o x p y p z p : The coordinate origin of the IMU frame coincides with that of the body frame, which is a unique frame of the rotating inertial navigation system. Axes are the directions of the sensitive axes of the inertial components. When each rotation angle of the shaft is zero, the p frame coincides with the b frame.

2.2. Error Model of IMU

The error models for the gyroscope and accelerometer are shown in Equations (1) and (2) [20]:
δ ω i p p = δ K g + G ω i p p + ε
δ f i p p = δ K a + A f i p p +
where δ ω i p p = δ ω i p x p δ ω i p y p δ ω i p z p T and δ f i p p = δ f i p x p δ f i p y p δ f i p z p represent the output error of the gyro and accelerometer, respectively; ω i p p = ω i p x p ω i p y p ω i p z p T and f i p p = f i p x p f i p y p f i p z p respectively represent angular velocity and acceleration input. The superscript p indicates that the measurement coordinate system is the p frame, the subscript ip indicates the angular velocity or acceleration of the p frame relative to the i frame, and x, y, z represent the corresponding coordinate axes. ε = ε x ε y ε z T and = x y z T represent the constant drift of gyro and the constant bias of accelerometer respectively. δ K g = d i a g δ K g x δ K g y δ K g z T and δ K a = d i a g δ K a x δ K a y δ K a z T respectively represent the scale factor error of the gyro and accelerometer. Gyro misalignment angle G and accelerometer misalignment angle A represent the included angle between the IMU-sensitive axis and the p frame, which could be expressed as:
G = 0 G x y G x z G y x 0 G y z G z x G z y 0
A = 0 A x y A x z A y x 0 A y z A z x A z y 0

2.3. Error Model of INS

The attitude error equation and velocity error equation of the INS are shown in Equations (5) and (6) [21]:
ϕ ˙ = ϕ × ω i n n + δ ω i n n C p n δ ω i p p
δ V ˙ n = ϕ × f n + δ V n × 2 ω i e n + ω e n n + V n × 2 δ ω i e n + δ ω e n n + C p n δ f i p p
where ϕ = ϕ E ϕ N ϕ U T is attitude error; δ V n = δ V E δ V N δ V U T is velocity error; ω i e n is Earth’s rotational motion, ω e n n is the rotation angular velocity of the n frame, and ω i n n = ω i e n + ω e n n ; f n = f E f N f U T is the specific force output under the navigation frame; V n = V E V N V U T indicates the speed of the carrier; and C p n represents the transformation matrix from the p frame to the n frame.

3. Error Analysis and Scheme Design for Alignment

3.1. Error Analysis for Alignment

In general, the carrier remains stationary during the alignment process. Theoretically, the speed output of the system is zero, but the actual speed output of INS is not zero due to the inertial sensor error. Therefore, the speed output is used as the measurement to inform the alignment process [22]. During the alignment process, the alignment accuracy of the horizontal attitude of the system depends on the equivalent bias of the horizontal accelerometers E and N , and the initial alignment accuracy of the azimuth depends on the equivalent bias of the east accelerometer E and the equivalent drift of the east gyro ε E . Considering that gyro accuracy is the main factor affecting the alignment accuracy, this paper mainly analyzes the effect of gyro error term by term.

3.1.1. Constant Bias

For the two-position scheme, during the alignment process, the IMU basically remains at the horizontal position, and remains stationary at the first position for a period of time. Then the IMU rotates 180° around the azimuth axis to reach the second position, and remains stationary for the same time. Assuming the initial attitude is [0°, 0°, 0°], when the coarse alignment error is ignored, the horizontal gyro bias before and after rotation is shown in Figure 1.
At the first stationary position, ε x t 1 is the easterly component of the x-axis gyro, i.e., the equivalent easterly gyro drift ω presents the rotation angular rate of IMU. At the second stationary position, ε x t 2 is the equivalent easterly gyro drift. During the rotation, the equivalent eastward gyro drift caused by the x-axis gyro is:
ε x t s = 0 t s ε x cos ω t d t t s
Considering the rotation process and the two stationary positions, the equivalent eastward gyro drift caused by the x-axis gyro in the two-position alignment process is:
ε x E = ε x t 1 + ε x t s + ε x t 2 = 0 t 1 ε x d t + 0 t s ε x cos ω t d t 0 t 2 ε x d t t
In Equation (8), the integration of the constant drift of the gyroscope during the stationary position and rotation process can be understood as the accumulation of the effect of the gyro drift over a period of time. By adding the integrals of the three processes and averaging them to the whole alignment time, the average equivalent result of gyro drift in the east direction during the whole alignment process can be obtained. ε x E can be expressed as equivalent eastward gyro drift caused by x-axis gyro bias during two-position alignment. It can be seen from Equation (8) that when the IMU is stationary at two positions for the same time, the equivalent eastward gyro bias caused by ε x is zero.
We can then analyze ε y in the same way. Since the y-axis initially points to the north, the equivalent eastward gyro drift caused by the y-axis gyro in the two-position alignment process is the equivalent component in the rotation process, namely:
ε y E = ε y t s = 0 t s ε y sin ω t d t t 2 ε y ω t
The estimation accuracy of azimuth error caused by horizontal gyro drift can be obtained by combining Equations (8) and (9):
δ ϕ ^ U ε = ε E ω i e cos L = 2 ε y ω t ω i e cos L
where δ ϕ ^ U ε is the estimation accuracy of azimuth error, L is the local latitude, and ω i e is the rotational velocity of earth.
It can be seen from Equation (10) that although the y-axis gyro bias will affect the alignment accuracy, its impact is very small. According to the above analysis process, it can be concluded that the effect of gyro constant bias on azimuth alignment accuracy can be well suppressed by IMU staying in the horizontal symmetric position at the same time.
A simple simulation analysis of single item error for the above theoretical analysis is then carried out. IMU stands still for 295 s at the first position, then rotates 180° around the azimuth axis to the second position within 10 s and stands still for 295 s. Assume that the constant drift of horizontal axis gyro is 0.01°/h, and ignore other component errors and coarse alignment errors. The initial attitude is [0°, 0°, 0°]. The simulation results are shown in Figure 2.
No matter the alignment process of the first position or the second position, the effect of the y-axis gyro drift on the azimuth alignment is not significant, and from the alignment results, there is basically no significant difference. The reason is that the y-axis points to the north at both positions, and the component in the east direction is always small. The effect of the rotation process on the east component is relatively small from the average to the whole alignment period. However, as shown in Figure 2, the effect of x-axis gyro drift is quite remarkable. The azimuth error of the first position is relatively large, but the azimuth error of the second position is gradually reduced. When the alignment time of the two positions is consistent, the azimuth error is basically consistent with the situation without component error. It can be seen that the second position gradually modulates the azimuth error generated in the alignment process of the first position.
Z-axis gyro drift will also cause azimuth alignment error. When the rotation axis of two-position alignment basically coincides with the horizontal plane, in the first position, the azimuth drift introduced by ε z is ε z t 1 , whose minimum variance estimation error is ε z t 1 / 2 ; similarly, in the second position, the minimum variance estimation error introduced by ε z is ε z t 2 / 2 . In the whole two-position alignment process, the azimuth estimation error caused by ε z is:
δ ϕ ^ U ε z = ε z t 2
The azimuth estimation error caused by the z-axis gyroscope with constant drift of 0.01°/h used in the later analysis in the 10-min alignment time is about 6″.
For the rotation scheme, IMU rotates in full circle around the azimuth axis. Similar to the analysis of the two-position scheme in the previous text, the equivalent eastward gyro drift caused by the horizontal gyro in the two-position alignment process is:
ε ω E = 0 t ε x cos ω t ε y sin ω t d t t
The estimation accuracy of azimuth error caused by horizontal gyro drift is:
δ ϕ ^ U ε = 0 t ε x cos ω t ε y sin ω t d t t ω i e cos L
Observing the integral term in the Equation (13), it can be observed that when the IMU rotates around the z-axis for the entire cycle, its sine and cosine integral results are zero, indicating that the bias of the horizontal gyro is basically completely suppressed when rotated for one cycle. To verify the theoretical analysis, a simple simulation is conducted. Assume that the constant drift of horizontal axis gyro is 0.01°/h, and ignore other component errors and coarse alignment errors. The initial attitude is [0°, 0°, 0°]. IMU undergoes reciprocating rotational motion around the z-axis with an angular velocity of 6 °/s. The simulation results are shown in Figure 3.
From Figure 3, it can be seen that the bias of the horizontal gyroscope has been well modulated.

3.1.2. Scale Factor Error

Because the alignment method analyzed in this paper includes the angular motion of IMU, the scale factor error cannot be ignored when some components with relatively poor scaling factor performance such as FOG [23] are used. Next, the influence of scale factor error during alignment is analyzed.
According to the attitude error equation, the attitude error introduced by the gyro scale factor error is:
ϕ ˙ K = C b n C p b [ δ K g ] ω i p p
where C b n is the attitude transfer matrix, and ω i p p is gyroscopic output angular velocity. In the two-position alignment scheme, assuming that the body frame coincides with the navigation frame, take IMU’s positive rotation around the azimuth axis as an example, in Equation (15):
C p b = cos ω t s sin ω t s 0 sin ω t s cos ω t s 0 0 0 1
where t s is the rotation time of IMU, and ω is the rotation angular velocity of IMU. The attitude error caused by the gyro scale factor error can be obtained as:
0 t ϕ ˙ k d t = ω i e cos L sin 2 ω t 2 ω δ K g x δ K g y ω i e cos L sin 2 ω t 4 ω δ K g y δ K g x + t ω i e cos L 2 δ K g x + δ K g y δ K g z ω i e sin L + ω t
Ignoring small quantities related to the rotation of the earth, during IMU rotation along the azimuth axis, the azimuth error caused by the gyro scale factor error can be simplified as:
ϕ K U = δ K g z ω t
According to Equation (17), IMU rotation will cause azimuth error by coupling with the gyro scale factor error. When δ K g z = 30   ppm , the attitude error caused by IMU rotation 180° is about 10″, as shown in Figure 4, which cannot be ignored. Therefore, for gyroscopes with unsatisfactory scale coefficient error, the effect of scale coefficient error of rotating axis gyroscopes should be taken into account when selecting the method of rotating IMU to improve the alignment accuracy. According to Equation (17), to suppress the influence of this error, the IMU needs to rotate the same angle in the opposite direction around the same rotation axis.
For the rotation scheme, when the IMU adopts the alignment scheme of a full cycle reciprocating rotation, according to Equation (17), the azimuth errors caused by scale factor error in the forward and reverse directions are as follows:
ϕ K U + = δ K g z ω t ϕ K U = δ K g z ω t
In this case, the estimation curve of the azimuth error will fluctuate up and down. Taking δ K g z = 30   ppm as an example, when rotating for one cycle, the fluctuation of the azimuth error is 30 ppm × 360° = 38.88″. When the IMU rotation direction changes, the direction of error accumulation also changes. The simulation results are shown in Figure 5.

3.1.3. Misalignment

When the IMU rotates, the attitude error introduced by the misalignment of gyros is:
ϕ ˙ G = C b n C p b G ω i p p = ω i e cos L G x y cos 2 ω t G y x sin 2 ω t + ω i e sin L + ω G x z cos ω t G y z sin ω t ω i e cos L sin 2 ω t G x y + G y x / 2 + ω i e sin L + ω G x z sin ω t + G y z cos ω t G z x ω i e cos L sin ω t + G z y ω i e cos L cos ω t
When the IMU rotates in a positive direction around the azimuth axis, it can be obtained when the small amount related to the Earth’s rotation is ignored ( ϕ G U = 0 ). When the IMU rotates around the azimuth axis, the misalignment of the gyro will not affect the azimuth alignment accuracy. For horizontal attitude alignment accuracy, IMU rotation will also introduce corresponding errors, as shown in Equation (20):
ϕ G E = G x z sin ω t G y z cos ω t ϕ G N = G x z cos ω t G y z sin ω t
The simple simulation analysis is carried out with 5″ misalignment of the gyro, and the results are shown in Figure 6. According to the simulation results, the size of the error is basically the same level as the misalignment of the gyro. The reason why the heading error in the simulation is not zero is that the influence of the Earth’s rotational velocity was ignored in the derivation. Under the premise that the calibration accuracy of the misalignment is high, the importance of this error could be ranked after the scale factor error.
For the rotation scheme, the simulation results are shown in the Figure 7.
Owing to the full rotation of the scheme, the misalignment of the azimuth gyro has also been modulated.
To sum up, the main gyro error items that affect the alignment accuracy include the constant bias of the horizontal gyro, the gyro scale factor error of the rotation axis, and the gyro misalignment related to the rotation axis. Among them, the constant bias is the main influence factor. The gyro scale factor error only introduces the relevant alignment error when the IMU is rotating, and the effect of gyro misalignment is the smallest.

3.2. Alignment Scheme Design

3.2.1. Alignment and Rotation Scheme Design

According to the theoretical analysis, the two-position alignment scheme only considers the effect of gyro constant bias on the alignment accuracy, and ignores the alignment error introduced by other related error items when IMU rotates. Therefore, when designing the multi-position alignment scheme, it is necessary to comprehensively consider the various errors of the components and the modulation effect of the rotation scheme on the various errors. According to the theoretical analysis in Section 2.1 and Section 2.2, the main factors that affect the initial alignment accuracy are the constant bias and the scale factor error of the gyro. To suppress the effect of constant bias, the IMU needs to remain at the same time in the symmetrical position, while to suppress the effect of the gyro scale factor error, the IMU needs to turn the same angle around the positive and negative directions of the rotation axis. Therefore, in this study, we designed a three-position alignment scheme that rotates twice around the azimuth axis. In this scheme, the IMU rotates 180° forward and backward around the azimuth axis to ensure the suppression effect of the gyro scale factor error. It stays at the first position and the third position for the same time, and stays at the second position for twice as long as at the first position to ensure the suppression effect on the constant bias of the gyro. During the alignment process, the z-axis gyro output of the two-position scheme, rotation scheme and three-position schemes is shown in Figure 8.

3.2.2. Alignment Algorithm Design

In the traditional initial alignment scheme, the initial alignment algorithm uses a 12-dimensional state vector Kalman filter for alignment solution. In this paper, due to the introduction of the research on the influence of the gyroscopic scale factor error of the rotating axis on the alignment accuracy, the gyro scale factor error of the rotating axis is introduced into the state vector of the initial alignment filter, namely:
X A = φ E φ N φ U δ V E δ V N δ V U ε x p ε y p ε z p x p y p z p δ K g z T
The state vectors in Equation (21) are attitude error, velocity error, gyroscope bias, accelerometer bias and scale factor error of the z-axis gyroscope.
The state equation and measurement equation are [24,25]:
X ˙ = F X + W
Z = H X + V
where W and V represent the state noise matrix and the measurement noise matrix, respectively. The velocity output calculated by the system is selected as the measurement information. According to the IMU error model and the INS error model, the state transition matrix in the filter can be obtained as:
F = F 1 F 2 F 3 0 7 × 6 0 7 × 6 0 7 × 1 13 × 13
In Equation (24), each sub-matrix is:
F 1 = 0 ω i e sin L ω i e cos L 0 1 R M + h 0 ω i e sin L 0 0 1 R N + h 0 0 ω i e cos L 0 0 tan L R N + h 0 0 0 f U f U 0 2 ω i e sin L 2 ω i e cos L f U 0 f E 2 ω i e sin L 0 0 f N f E 0 2 ω i e cos L 0 0
F 2 = C p n 0 3 × 3 0 3 × 3 C p n
F 3 = C p n 0 0 ω i p z p T 0 1 × 3 T
The measurement matrix can be represented as H = [ 0 3 × 3 I 3 × 3 0 3 × 7 ] ; I 3 × 3 represents a third-order identity matrix, and 0 m × n represents an m × n zero matrix.

4. Simulation and Experimental Analysis

4.1. Simulation Analysis

In order to evaluate the alignment accuracy of the three-position alignment scheme, a comparative simulation experiment was carried out for three alignment schemes: two-position alignment, rotation alignment and three-position alignment. Assume that the initial attitude error is [0.5°, 0.5° 1°] and the position is [108.9° E, 34.2° N, 400 m], and the alignment time is 600 s. IMU error parameter settings in the simulation are shown in Table 1. In total, 20 simulation tests on the three alignment schemes are carried out, and the root mean square curves of the system attitude error during the alignment process are drawn, as shown in Figure 9 and Figure 10.
It can be seen from Figure 9 and Figure 10 that the horizontal attitude error converges rapidly as the IMU rotates in the two-position scheme and three-position scheme. The reason is that the horizontal attitude alignment accuracy is mainly affected by the horizontal equivalent accelerometer bias. As the accelerometer bias is rapidly and accurately estimated, the horizontal attitude error converges rapidly. In the rotation alignment scheme, due to the continuous rotation of IMU, it can be seen from the analysis in Section 3.1.3 that although the horizontal accelerometer bias can be accurately estimated and compensated, the gyro misalignment will cause the horizontal attitude to be sinusoidal. The misalignment of the gyro is also the reason why the horizontal error curve of the two-position scheme and the three-position scheme cannot be further converged. However, in terms of the magnitude of the alignment error caused by misalignment of the gyro, its effect is not as obvious as that of gyro constant bias and gyro scale factor error.
As can be seen from Figure 10, from the final result, the azimuth alignment accuracy of the three-position alignment scheme is similar to that of the rotation alignment scheme, which is better than the two-position alignment scheme. Among them, the azimuth error curve of the two-position alignment scheme fluctuates obviously when the alignment time is about 300 s. According to Section 3.1.2, this is caused by the coupling of the gyro scale factor error and IMU rotation motion, which is consistent with the fluctuation of error curve of the three-position scheme around 180 s and 440 s. But the two rotation directions of the three-position alignment scheme are opposite, which effectively suppresses the influence of the gyro scale coefficient error on the alignment accuracy. The continuous fluctuation of the angular error curve of the rotation alignment scheme is caused by the continuous rotation of the IMU. Since the gyro scale factor error of the rotating axis cannot be accurately estimated and compensated, the continuous rotation of IMU will cause continuous fluctuation of azimuth.
The azimuth accuracy of the rotation alignment scheme at the end of alignment is similar to that of the three-position alignment scheme. However, in some special applications that require a rapid response, such as self-propelled guns, missile launchers, etc., the rotation alignment scheme may encounter some problems. In such cases, the vehicle may be temporarily transferred before the preset alignment time. In the circumstances that the IMU is not at the end of a complete cycle (forward and reverse for one cycle), the adoption of the rotating alignment scheme means that the scale factor error of the gyroscope is not fully modulated, so the alignment results obtained at this time are not ideal. It can also be seen from the azimuth error angle curve of rotation alignment in Figure 10 that the scheme fluctuates greatly during the alignment process. Once the alignment needs to be temporarily terminated in advance, the credibility of the obtained alignment results will be greatly reduced.
A statistical analysis of the standard deviation of the estimation error was conducted. From Figure 11, Figure 12 and Figure 13, it can be seen that the standard deviation curve of the estimation error of accelerometer bias can converge in the vicinity of zero, indicating that the filter has a good estimation effect on accelerometer bias, while the filter has poor estimation effect on gyro bias and z-axis gyro scale factor error. For gyro bias, these schemes all adopt the idea of error modulation to suppress their influence. Due to the inability to estimate the scale factor error of the z-axis gyroscope, the two-position scheme generates additional errors. These analyses are also consistent with the results of the analysis in Section 3.1.

4.2. Test Verification

In order to verify the research content of this paper, the two-position alignment, rotation alignment and three-position alignment schemes were compared and analyzed through experiments. The experimental equipment mainly includes a three-axis rotary fiber optic gyro inertial navigation system, as shown in Figure 14. Inertial component parameters adopted by the system are basically consistent with the simulation parameters. Considering the high accuracy of the system and the difficulty in obtaining the high-precision attitude reference information matched with it, ten experiments were carried out on three alignment schemes respectively, then the standard deviations of the alignment results were calculated, and the statistical results were compared and analyzed. The results are shown in Table 2.
It can be seen from Table 2 that the horizontal alignment accuracy of the three alignment schemes is basically the same. The reason is that the three alignment schemes can accurately estimate the bias of the horizontal accelerometer, which is the main error that affects the accuracy of horizontal alignment. The estimation result of the bias of the horizontal accelerometer is shown in Figure 15. From the perspective of azimuth alignment accuracy, the three-position alignment scheme designed in this paper can reach the same accuracy level as the rotation alignment scheme, which is significantly improved compared with the two-position alignment scheme. According to the analysis in Section 2 of this paper, the reason is that in the three-position alignment scheme, the equivalent biases of the horizontal gyros in the second position and the first and third positions are modulated. The two-position scheme and the rotation scheme can also achieve the corresponding error modulation. However, the three-position scheme suppresses the influence of the gyro scale factor error on azimuth alignment accuracy through two angular movements with opposite directions and the same rotation angle, while the two-position alignment scheme only rotates once, and the azimuth error introduced by the gyro scale factor error cannot be modulated. In the rotation alignment scheme, the coupling between the continuous rotation of IMU and the error of the gyro scale factor causes the continuous fluctuation of the azimuth error. The above three schemes cannot estimate the gyro scale factor error of the rotating axis, as shown in Figure 16. Therefore, the influence of this error term can only be suppressed by error modulation. The above experimental results are basically consistent with the theoretical analysis and simulation results of this paper, which verifies the effectiveness of the alignment error suppression theory and the alignment scheme designed in this paper.

5. Conclusions

This paper investigates the influencing factors and principles for error suppression in multi-position initial alignment accuracy. Using the classic two-position alignment scheme, an analysis is conducted on the primary component error terms and their respective magnitudes impacting alignment accuracy during IMU rotation. Notably, the bias of gyroscopes and the gyro scale factor error emerge as principal factors influencing alignment accuracy within the two-position alignment paradigm. Building upon an elucidation of the suppression principle of IMU rotation on initial alignment errors, a novel three-position initial alignment scheme is proposed, accompanied by a detailed analysis of its error suppression mechanism. Simulation and experimental outcomes corroborate the efficacy of the proposed scheme in effectively mitigating the adverse effects of gyro scale factor error and gyro bias on initial alignment accuracy. Furthermore, the error suppression strategies outlined herein offer valuable insights for the design of initial alignment schemes tailored to diverse devices operating under corresponding conditions.

Author Contributions

Conceptualization, B.G. and D.W.; methodology, B.G. and Q.F.; software, Z.L.; validation, Z.L.; formal analysis, B.G.; data curation, Z.L.; writing—original draft preparation, B.G.; writing—review and editing, B.G. and D.W.; project administration, Q.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Aviation Science Foundation of China (2023Z022053001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

Author Zhongping Liu was employed by the company Shaanxi Aircraft Corporation. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Chang, Z.; Zhang, Z.; Zhou, Z.; Chen, H.; Zhao, J. Alignment method of rotating SINS with swinging base. Syst. Eng. Electron. 2019, 41, 2087–2091. [Google Scholar]
  2. Emel’yantsev, G.; Stepanov, O.; Stepanov, A.; Blazhnov, B.; Dranitsyna, E.; Evstifeev, M.; Eliseev, D.; Volynskiy, D. Integrated GNNS/IMU Gyrocompass with Rotating IMU. Development and Test Results. Remote Sens. 2020, 12, 3736. [Google Scholar] [CrossRef]
  3. Lin, Y.; Miao, L.; Zhou, Z. A high-accuracy initial alignment method based on backtracking process for strapdown inertial navigation system. Measurement 2022, 201, 111712. [Google Scholar] [CrossRef]
  4. Sun, W.; Li, S.; Li, R. Analysis of the Initial Alignment of Rotary Fiber Inertial Navigation System. Chin. J. Sens. Actuators 2014, 27, 1082–1087. [Google Scholar]
  5. Chen, J.; Cheng, X.; Li, J. Rapid alignment algorithm for rotary INS in inertial frame. Transducer Microsyst. Technol. 2017, 36, 142–145. [Google Scholar]
  6. Zhang, Q.; Wang, L.; Liu, Z.; Feng, P. An Accurate Calibration Method Based on Velocity in a Rotational Inertial Navigation System. Sensors 2015, 15, 18443–18458. [Google Scholar] [CrossRef] [PubMed]
  7. Li, J.; Su, L.; Wang, F.; Li, K.; Zhang, L. An Improved Online Fast Self-Calibration Method for Dual-Axis RINS Based on Backtracking Scheme. Sensors 2022, 22, 5036. [Google Scholar] [CrossRef] [PubMed]
  8. Lu, Y.; Wang, W.; Liu, Y.; Guo, Z. An Improved Rotational Modulation Scheme for Tri-Axis Rotational Inertial Navigation System (RINS) with Fiber Optic Gyro (FOG). Appl. Sci. 2023, 13, 8394. [Google Scholar] [CrossRef]
  9. Liu, B.; Wei, S.; Lu, J.; Wang, J.; Su, G. Fast Self-Alignment Technology for Hybrid Inertial Navigation Systems Based on a New Two-Position Analytic Method. IEEE Trans. Ind. Electron. 2020, 67, 3226–3235. [Google Scholar] [CrossRef]
  10. Xu, C.; Miao, L.; Zhou, Z.; Lin, Y. A rapid and high-accuracy rotation alignment method based on bidirectional process for dual-axis rotational inertial navigation system. Meas. Sci. Technol. 2020, 31, 125106. [Google Scholar] [CrossRef]
  11. Wang, L.; Wu, W.; Wei, G.; Pan, X.; Lian, J. Navigation Information Fusion in a Redundant Marine Rotational Inertial Navigation System Configuration. J. Navig. 2018, 71, 1531–1552. [Google Scholar] [CrossRef]
  12. Mao, S.; Hu, P.; Ding, X.; Tan, J. Parameter correction method for dual position-sensitive-detector-based unit. Appl. Opt. 2016, 55, 4073–4078. [Google Scholar] [CrossRef]
  13. Fang, J.; Han, X.; Yang, F. Two-position alignment method of airborne position and orientation system. J. Chin. Inert. Technol. 2013, 21, 318–323. [Google Scholar]
  14. Lai, J.; Xiong, J.; Liu, J.; Jiang, B. Improved arithmetic of two-position fast initial alignment for sins using unscented kalman filter. Int. J. Innov. Comput. Inf. Control. 2012, 8, 2929–2940. [Google Scholar]
  15. Zhang, L.; Wu, W.; Wang, M. Rapid SINS Two-Position Ground Alignment Scheme Based on Piecewise Combined Kalman Filter and Azimuth Constraint Information. Sensors 2019, 19, 1125. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, D.; Li, G. The improved method of the SINS initial alignment based on single-axis rotation. J. Harbin Eng. Univ. 2014, 35, 925–930. [Google Scholar]
  17. Wei, D.; Li, S.; Fu, Q. A New Initial Alignment Scheme for Dual-Axis Rotational Inertial Navigation System. IEEE Trans. Instrum. Meas. 2022, 71, 8503010. [Google Scholar] [CrossRef]
  18. Wang, M.; Wang, L.; Han, H. Research on Innovative Self-Calibration Strategy for Error Parameters of Dual-Axis RINS. IEEE Sens. J. 2023, 21, 27918–27927. [Google Scholar] [CrossRef]
  19. Seo, Y.-B.; Yu, H.; Ryu, K.; Lee, I.; Oh, J.; Kim, C.; Lee, S.J.; Park, C. Analysis of Gyro Bias Depending on the Position of Inertial Measurement Unit in Rotational Inertial Navigation Systems. Sensors 2022, 22, 8355. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, L.; Liu, W.; Lian, J.; Wu, M. Depression of ring laser gyro drift using rotation modulation technique. J. Chin. Inert. Technol. 2011, 19, 515–520. [Google Scholar]
  21. Wang, B.; Weng, H.; Wang, Z.; Zhang, Y.; Hu, X. Calibration method for two-axis indexing RLG-SINS with speed damping. J. Chin. Inert. Technol. 2019, 41, 421–425. [Google Scholar]
  22. Yan, G.; Weng, J.; Bai, L.; Qin, Y. Initial in-movement alignment and position determination based on inertial reference frame. Syst. Eng. Electron. 2011, 33, 618–621. [Google Scholar]
  23. Yu, H.; Wang, W.; Huang, L. Improved performance of scale factory linearity in closed-loop IFOG. J. Chin. Inert. Technol. 2007, 15, 449–451. [Google Scholar]
  24. Fu, Q. Key Technologies for Vehicular Positioning and Orientation System. Ph.D. Thesis, Northwestern Polytechnical University, Xi’an, China, 2014. [Google Scholar]
  25. Dranitsyna, E.V. IMU Calibration Using Sins Navigation Solution: Selection of the Rate Table Motion Scenario. In Proceedings of the 24th Saint Petersburg International Conference on Integrated Navigation Systems, ICINS 2017, St. Petersburg, Russia, 29–31 May 2017. [Google Scholar]
Figure 1. Rotation diagram of ε x and ε y during two-position alignment.
Figure 1. Rotation diagram of ε x and ε y during two-position alignment.
Sensors 24 01938 g001
Figure 2. Influence of x-axis gyro bias on azimuth error (two-position scheme).
Figure 2. Influence of x-axis gyro bias on azimuth error (two-position scheme).
Sensors 24 01938 g002
Figure 3. Influence of horizontal gyro bias on azimuth error (rotation scheme).
Figure 3. Influence of horizontal gyro bias on azimuth error (rotation scheme).
Sensors 24 01938 g003
Figure 4. Influence of the scale factor error of azimuth gyro bias on the azimuth error (two-position scheme).
Figure 4. Influence of the scale factor error of azimuth gyro bias on the azimuth error (two-position scheme).
Sensors 24 01938 g004
Figure 5. Influence of the scale factor error of azimuth gyro bias on the azimuth error (rotation scheme).
Figure 5. Influence of the scale factor error of azimuth gyro bias on the azimuth error (rotation scheme).
Sensors 24 01938 g005
Figure 6. Influence of the misalignment of the azimuth gyro on the azimuth error (two-position scheme).
Figure 6. Influence of the misalignment of the azimuth gyro on the azimuth error (two-position scheme).
Sensors 24 01938 g006
Figure 7. Influence of the misalignment of the azimuth gyro on the azimuth error (rotation scheme).
Figure 7. Influence of the misalignment of the azimuth gyro on the azimuth error (rotation scheme).
Sensors 24 01938 g007
Figure 8. Schematic diagram of three alignment schemes.
Figure 8. Schematic diagram of three alignment schemes.
Sensors 24 01938 g008
Figure 9. Curves of the horizontal attitude error of three alignment schemes.
Figure 9. Curves of the horizontal attitude error of three alignment schemes.
Sensors 24 01938 g009
Figure 10. Curves of the azimuth error of three alignment schemes.
Figure 10. Curves of the azimuth error of three alignment schemes.
Sensors 24 01938 g010
Figure 11. Standard deviation of the estimation error of gyro bias.
Figure 11. Standard deviation of the estimation error of gyro bias.
Sensors 24 01938 g011
Figure 12. Standard deviation of the estimation error of accelerometer bias.
Figure 12. Standard deviation of the estimation error of accelerometer bias.
Sensors 24 01938 g012
Figure 13. Standard deviation of the estimation error of the scale factor error of Z-axis gyro.
Figure 13. Standard deviation of the estimation error of the scale factor error of Z-axis gyro.
Sensors 24 01938 g013
Figure 14. Rotational inertial navigation system (RINS) for the experiment.
Figure 14. Rotational inertial navigation system (RINS) for the experiment.
Sensors 24 01938 g014
Figure 15. Estimation curve of horizontal accelerometer bias.
Figure 15. Estimation curve of horizontal accelerometer bias.
Sensors 24 01938 g015
Figure 16. Estimation curve of gyroscope scale factor error of the rotating axis.
Figure 16. Estimation curve of gyroscope scale factor error of the rotating axis.
Sensors 24 01938 g016
Table 1. IMU error parameter settings in the simulation.
Table 1. IMU error parameter settings in the simulation.
Error TermsParameters
GyroConstant bias0.002°/h
Random walk0.0003°/√h
Scale factor error30 ppm
Misalignment2″
Driving noise (first-order Markov process)0.0003°/h
Correlated time3000 s
AccelerometerConstant bias20 μg
Random walk5 μg/√Hz
Scale factor error30 ppm
Misalignment2″
Table 2. Standard deviation statistics of the alignment results of three alignment schemes.
Table 2. Standard deviation statistics of the alignment results of three alignment schemes.
Attitude AngleTwo-Position SchemeRotation SchemeThree-Position Scheme
pitch4.59″2.19″3.51″
roll1.32″1.58″2.03″
yaw20.32″14.63″12.83″
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guan, B.; Liu, Z.; Wei, D.; Fu, Q. Design of a Multi-Position Alignment Scheme. Sensors 2024, 24, 1938. https://doi.org/10.3390/s24061938

AMA Style

Guan B, Liu Z, Wei D, Fu Q. Design of a Multi-Position Alignment Scheme. Sensors. 2024; 24(6):1938. https://doi.org/10.3390/s24061938

Chicago/Turabian Style

Guan, Bofan, Zhongping Liu, Dong Wei, and Qiangwen Fu. 2024. "Design of a Multi-Position Alignment Scheme" Sensors 24, no. 6: 1938. https://doi.org/10.3390/s24061938

APA Style

Guan, B., Liu, Z., Wei, D., & Fu, Q. (2024). Design of a Multi-Position Alignment Scheme. Sensors, 24(6), 1938. https://doi.org/10.3390/s24061938

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