Next Article in Journal
Towards Single Camera Human 3D-Kinematics
Next Article in Special Issue
A Scalable Distributed Control Algorithm for Bearing-Only Passive UAV Formation Maintenance
Previous Article in Journal
Wireless Sensing Technology Combined with Facial Expression to Realize Multimodal Emotion Recognition
Previous Article in Special Issue
A Context-Aware Smartphone-Based 3D Indoor Positioning Using Pedestrian Dead Reckoning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Rapid Self-Alignment Strategy for a Launch Vehicle on an Offshore Launching Platform

1
School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
2
Institute of Modern Optics, Nankai University, Tianjin 300350, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(1), 339; https://doi.org/10.3390/s23010339
Submission received: 3 November 2022 / Revised: 23 December 2022 / Accepted: 24 December 2022 / Published: 28 December 2022

Abstract

:
To reduce the impact of offshore launching platform motion and swaying on the self-alignment accuracy of a launch vehicle, a rapid self-alignment strategy, which involves an optimal combination of anti-swaying coarse alignment (ASCA), backtracking navigation, and reverse Kalman filtering is proposed. During the entire alignment process, the data provided by the strapdown inertial navigation system (SINS) are stored and then applied to forward and backtrack self-alignment. This work elaborates the basic principles of coarse alignment and then analyzes the influence of ASCA time on alignment accuracy. An error model was built for the reverse fine alignment system. The coarse alignment was carried out based on the above work, then the state of the alignment system was retraced using the reverse inertial navigation solution and reverse Kalman filtering with the proposed strategy. A cycle-index control function was designed to approximate strict backtracking navigation. Finally, the attitude error was compensated for after the completion of the first and the last forward navigation. To demonstrate the effectiveness of the proposed strategy, numerical simulations were carried out in a scenario of launch vehicle motion and swaying. The proposed strategy can maximize the utilization of SINS data and hence improve the alignment accuracy and further reduce the alignment time. The results show that the fully autonomous alignment technology of the SINS can replace the complex optical aiming system and realize the determination of the initial attitude of a launch vehicle before launch.

1. Introduction

The initial alignment (IA) accuracy of a launch vehicle (LV) directly affects the navigation accuracy, and the initial alignment time directly affects the launch vehicle test and launch process. Due to the influence of disturbance factors, such as wind and waves at sea [1], offshore launches of launch vehicles are subject to more severe swaying than land-based launches. It is necessary to solve the self-alignment problem of the launch vehicle under motion and swaying conditions and achieve considerable accuracy in a very short launch period. The alignment is one of the most important stages of a navigation system [2]. The optical sighting or self-sighting algorithm used for initial alignment of land-based launch vehicles cannot be directly applied to the initial alignment at sea [3]. In case of land-based launch, the optical sighting equipment is fixed and the launch vehicle is slightly swayed. If the launch is on an offshore launching platform, both the optical aiming equipment and the launch vehicle are swaying at the same time, it will be impossible to calibrate. How to propose a high-precision initial alignment method suitable for sea launch in response to the characteristics of the marine motion environment is one of the technical difficulties for the launch of launch vehicles [3].
At present, optical sighting is usually used at home and abroad for azimuth reference transfer, and filter compensation or traditional transfer alignment method is used to complete the horizontal reference transfer [3]. The fundamental principle of the common method for initial alignment of sea launches of launch vehicles is shown in Figure 1. The principle of this method is simple and reliable, but the application cost is high, and the accuracy is slightly worse.
Conventional self-alignment methods are performed sequentially in chronological order, and data processing is often one-off; the data of the previous stage are not used in the later stage, and there is no storage requirement, resulting in inadequate information mining [4]. In contrast, the computing power of current navigation computers is strong enough and the capacity is large enough to save the SINS data, which can then be repeatedly processed and analyzed, providing another method of improving the accuracy and speed of an LV’s self-alignment at sea.
The coarse alignment (CA) process mainly accomplishes the coarse estimation of the unknown attitude angle [5]. Qin et al. [6] proposed a coarse alignment algorithm based on gravitational acceleration, which effectively isolated the angular swaying interference by introducing the assumption of a solidification coordinate system (CS) and tracking the attitude change in real time according to the gyro output. Bai et al. [7] demonstrated the effectiveness of the gravitational acceleration-based alignment algorithm through verification experiments. The measured angular velocity of the earth’s rotation and gravitational acceleration can be used in the coarse alignment stage. Swaying is a perturbation situation that mainly affects orientation changes [8]. With the help of the idea of decomposing the inertial coordinate system, the gravity vectors in the inertial space at different moments are not co-linear, which can be solved by the double-vector attitude determination [9,10] or multi-vector attitude determination algorithm, which isolates the swaying interference [11,12,13]. To improve the accuracy and efficiency of latitude and attitude estimation by fully utilizing inertial measurements, an anti-swaying self-alignment method based on the inertia tensor of gravity observation vector was presented in Reference [2].
Coarse alignment results can be improved by applying fine alignment (FA) at the completion of the coarse alignment phase. Fine alignment improves the accuracy of the fusion process by using external sensors or information [14,15]. Tsukerman et al. [16] derived an analytic assessment of the steady-state characteristics of the fine alignment procedure.
To improve the alignment performance by taking full advantage of all sensor data in a limited amount of time [17,18,19], a real-time information-reusing alignment strategy was proposed for rotating an inertial navigation system in Reference [20]. Li et al. [21] designed a backtracking navigation scheme for INS (Inertial Navigation System)/DVL (Doppler Velocity Log) integration. Yan et al. [22] proposed an inter-travel alignment scheme based on reverse data processing, which organically combined the initial alignment of the moving base with the position navigation through the reverse navigation algorithm, which first carried out the initial alignment of the moving base attitude and the sensor sampling data storage of the SINS (Strapdown Inertial Navigation System) gyrocompass and then used the reverse and forward dead reckoning algorithms to obtain the position information, which solved the previous problem that the position needed to be corrected using landmark points after inter-travel alignment. Liang et al. [23] applied the reverse navigation algorithm to the static base alignment of the SINS, which effectively shortened the alignment time and improved the alignment accuracy. Wang et al. [24] introduced the reverse Kalman filtering technique, which processed the collected historical data several times using the forward and reverse Kalman filtering algorithms and completed the reverse fine alignment and forward combined navigation on the basis of the coarse alignment results. In Reference [25], aiming at the problem of long alignment time for traditional methods, a forward-forward backtracking initial alignment method under swaying base conditions was proposed to achieve the purpose of rapid initial alignment, which did not consider linear vibration.
The incremental output SINS is taken as the research object. Enlightened by [19], a self-alignment strategy based on the reverse data processing technology is proposed, which uses pure mathematical analysis to complete the alignment of azimuth and horizontal directions simultaneously without the need for complex optical transmission equipment, thus achieving the purpose of shortening the alignment time and improving the alignment accuracy. The main innovations of this paper are as follows:
  • Reverse Kalman filtering algorithm is designed, and the backtracking navigation is completed by the combination of the reverse navigation solution and reverse Kalman filtering;
  • A cycle-index control function is designed for approximate strict backtracking navigation;
  • In addition to the attitude correction at the initial moment (IM), attitude error compensation is completed twice during the subsequent alignment.
Compared to the existing alignment methods, firstly, using of the gravity projection along the inertial system, with reverse Kalman filter technique, could avoid the effects of disturbed accelerates and angular velocities on the SINS alignment. Secondly, attitude error compensation, to a certain extent, can achieve real-time tracking, which is the most obvious difference from the existing methods. Then, the designed cycle-index control function enables the control of the number of forward and reverse navigation solutions, thus controlling the amount of computation. Simulations and experiments show that the strategy proposed is effective. In addition, the strategy can be applied to engineering practice, as it effectively improves alignment accuracy and reduces alignment time compared with the traditional alignment methods.

2. Problem Description

Coordinate systems and Attitude Angles are described as follows:
Geocentric Inertial CS ( i system): o i x i y i z i . The origin of the CS is located at the center ( o i ) of the earth. The o i x i -axis points to the vernal equinox, and the o i z i -axis is along the earth’s axis of rotation. The o i y i -axis and the o i x i -axis and o i z i -axis form a right-handed CS.
Earth CS ( e system): ο i x e y e z e . The CS is fixed to the earth, and the origin is at the center ( o i ) of the earth. The ο i x e -axis passes through the intersection of the prime meridian and the equator, the ο i z e -axis passes through the north pole and the ο i y e -axis passes through the intersection of the meridian 90°E and the equator.
Geographic CS ( g system): ο b x g y g z g . The origin of the CS is located at the center of mass of the carrier, and the ο b x g -axis, ο b y g -axis and ο b z g -axis point to the east, north and up, respectively.
Navigation CS ( n system): ο n x n y n z n . The CS is the reference CS for the navigation solution, and the geographic CS is used as the navigation CS in this paper.
Calculation of Navigation CS ( n system): ο n x n y n z n . The CS is the navigation CS established according to the navigation parameters obtained from the inertial navigation solution.
Inertial Measurement Unit (IMU) CS ( s system): ο s x s y s z s . The origin of the CS is located at the center of mass of the IMU, and the three axes point to the direction of the IMU sensitive axis. It is assumed that the body CS is rigidly connected to the IMU CS.
IMU CS at the IM ( s 0 system) coincides with the IMU CS ( s system) at the instance of the IA, and there is no rotation relative to inertial space.
Navigation Inertial CS at the IM ( n 0 system) coincides with the navigation CS ( n system) at the instance of the IA, and there is no rotation relative to inertial space.
Launch CS ( n t system): ο t x t y t z t . The origin of the CS is located at the center of the launch point, the ο t y t -axis is horizontally forward along the ballistic direction, the ο t z t -axis is vertically upward, the ο t x t -axis is perpendicular to the ballistic plane to the right and the o t y t -axis and o t z t -axis form a right-handed CS. Obviously the ο t y t z t -plane is the ballistic plane. The angle between the ballistic plane (or ο t y t -axis) and the local geographic north direction is usually called the launch angle, which is recorded as A 0 .
The Inertial Navigation System (IMU) is installed on the LV, and its axial definition is shown in Figure 2. When the LV is “lying down” horizontally, the three axes are along the lateral, longitudinal and vertical directions respectively, i.e., “right-front-up”. Euler angles are defined as “132”:
Pitch angle ( θ ): After a coordinate translation so that the origin of the navigation CS coincides with the origin of the IMU CS, the pitch angle is the angle between the projection line of the longitudinal axis ( ο b y b -axis) of the LV in the ballistic plane and the ο n y n -axis. The angle is positive if the longitudinal axis of the offshore platform points above the horizontal plane; otherwise, it is negative, and the value range is θ [ 0 , 180 ] .
Yaw angle ( ψ ): After a coordinate translation so that the origin of the navigation CS coincides with the origin of the IMU CS, the yaw angle is the angle between the longitudinal axis of the LV and the ballistic plane. Looking at the angle plane, if the rotation from the ο y n -axis to the ο y b -axis is counterclockwise, the angle is positive; otherwise, it is negative, and the value range is ψ [ 90 , 90 ] .
Roll angle ( γ ): After a coordinate translation so that the origin of the navigation CS coincides with the origin of the IMU CS, the roll angle is the angle between the vertical axis of the LV and the plumb plane containing the longitudinal axis of the LV. Looking forward along the longitudinal axis from the tail of the LV, if the ο z b -axis is on the right side of the vertical plane, the angle is positive; otherwise, it is negative. The value range is γ [ 90 , 90 ] .
The conversion relationship applied throughout the article is described as follows: C B N is the representation of the attitude of B system with respect to N system. It is equivalent to say that it is the transformation matrix from N system to B system.
The differential equation for the specific force of SINS in a launch vehicle can be described as:
ν ˙ e n n ( t ) = C b n ( t ) f s f b ( t ) ( 2 ω i e n ( t ) + ω e n n ( t ) ) × ν e n n ( t ) + g n ( t )
where f s f b ( t ) is the ideal output real-time specific force of the accelerometer; ω i e n ( t ) represents the real-time angular rate of rotation of the navigation system due to Earth’s rotation; ω e n n ( t ) represents the n system rotation of the SINS moving near the Earth’s surface due to the curvature of the Earth’s surface; g n ( t ) represents the gravity vector in the n system [27].
Assuming the offshore platform is moored. Influenced by the sea state, the motion state of the launching platform can be decomposed into swaying motion and linear motion. Under such conditions, the initialization of the SINS in a launch vehicle will be obtained by the self-alignment strategy designed later.

3. Design of Self-Alignment Method

To realize the rapid self-alignment process of the LV launched on the offshore platform, this research established a backtracking system error model, and a reverse Kalman filtering algorithm was designed. The self-alignment process is divided into three parts: ASCA, backtracking navigation FA and forward FA.

3.1. Self-Alignment Strategy of LV

The designed self-alignment strategy of the LV launched on the offshore platform is shown in Figure 3. In Figure 3, C ^ b , C n represents the initial attitude direction cosine matrix calculated at the end of the CA. C ^ b , R n represents the initial attitude direction cosine matrix calculated at the end of the backtracking FA. C ^ b n ( k ) represents the initial attitude direction cosine matrix obtained from each backtracking FA estimation.
Firstly, at the IM t 0 of the alignment, the CA process is started by the ASCA method while storing IMU data in the mooring condition. The CA process is completed at the moment t 1 . The computed attitude matrix C ^ b , C n , which is obtained by the CA, is used as the initial value of the first backtracking FA, and the backtracking FA uses the collected IMU data information to return the state of the LV to the IM t 0 . At the same time, the attitude information estimated by the BN is combined with the known velocity information at the starting point to correct the inertial navigation. The first backtracking FA process is completed, and the attitude information with a certain accuracy is obtained. Immediately after the fast forward navigation to return the system state to the t 1 moments, the attitude errors are compensated using the data stored from t 1 to t 2 moment. Then, the backtracking-forward cyclic FA starts from moment t 2 . When the last forward FA returns to moment t 2 , the data stored after moment t 2 are used to catch up with the normal real-time computation until the end of the alignment. The navigation computer maintains high speed computing during the whole process. The entire FA process is completed, and higher accuracy attitude information can be obtained.
The BN algorithm is formally identical to the forward navigation algorithm. As long as the forward saved gyro data and the angular rate of Earth rotation, etc., are inverted and the forward final value is set to the backtracking initial value, BN can be realized. During the backtracking solution, all sensor data must be stored. Therefore, in the IA, if the data processed by the backtracking continues to be used as a new set of data, or if the forward and backtracking attitude solutions are combined and the data are used repeatedly, this is equivalent to increasing the quantity of data.

3.2. Anti-Swaying Coarse Alignment Algorithm

Through the idea of decomposition of inertial solidification CS, the initial attitude matrix C s n is decomposed into
C s n ( t ) = C n 0 n ( t ) C s 0 n 0 C s s 0 ( t )
C n 0 n represents the transformation matrix from the navigation inertial CS ( n 0 -system) at the IM to the navigation CS ( n -system); C s 0 n 0 represents the transformation matrix from the IMU CS ( s 0 -system) at the IM to the navigation inertial CS ( n 0 -system) at the IM; C s s 0 represents the transformation matrix from the IMU CS ( s -system) to the IMU CS ( s 0 system) at the IM.
According to the angular velocity of the earth’s rotation and the time interval between the two moments, the relationship between the two moments of the navigation CS can be obtained; according to the projection of the gravitational acceleration vector at the two moments in the navigation CS at the IM of alignment, and the projection of the specific force measurement value in the initial IMU CS, the conversion relationship between the IMU CS and the navigation CS at the IM can be solved by the double-vector attitude determination algorithm. From the projection and conversion relationship represented by the above three matrices, the conversion relationship between the IMU CS and the navigation CS at the end of the alignment can be determined, and then the CA is completed.

3.3. Backtracking Navigation Algorithm

In order to realize the sampling and utilization of the output signals of the gyroscope and accelerometer without an omission in time, the system collects the angular incremental output of the gyroscope and the velocity incremental output of the accelerometer within a defined time interval, rather than the angular velocity output and the specific force output at discrete time points. The attitude, velocity and position update equations of the incremental SINS inertial guide can be generally expressed as follows.

3.3.1. Forward Navigation [27]

(1)
Attitude update
C b f ( m ) n f ( m ) = C n f ( m 1 ) n f ( m ) C b f ( m 1 ) n f ( m 1 ) C b f ( m ) b f ( m 1 )
C n f ( m 1 ) n f ( m ) = M R V T ( T m ω i n ( m 1 / 2 ) n f ) ,     C b f ( m ) b f ( m 1 ) = M R V ( ϕ i b ( m ) b f )
where C b f ( m 1 ) n f ( m 1 ) and C b f ( m ) n f ( m ) represent the SINS attitude matrix at t m 1 and t m ( m = 1 , 2 , 3 , n m ), respectively. For the navigation update period [ t m 1 , t m ] , T m = t m t m 1 . t ( m 1 / 2 ) = ( t m + t m 1 ) / 2 denotes the midpoint moment. ω i n ( m 1 / 2 ) n f represents the rotation of the n system with respect to the i system at the moment of the midpoint in the forward navigation. The subscript ( m 1 / 2 ) indicates that they are the variables at the midpoint moment. Δ θ m 1 denotes the angular increment in the time interval [ t m 1 t ( m 1 / 2 ) ] , and Δ θ m 2 denotes the angular increment in the time interval [ t ( m 1 / 2 ) t m ] . The equivalent rotation vector can be obtained using the two-subsample conical error compensation algorithm:
ϕ i b ( m ) b f = ( Δ θ m 1 + Δ θ m 2 ) + 2 3 Δ θ m 1 × Δ θ m 2
M R V can be obtained from the equivalent rotation vector ϕ ( ϕ = | ϕ | ):
M R V = I + sin ϕ ϕ ( ϕ × ) + 1 cos ϕ ϕ 2 ( ϕ × ) 2
(2)
Velocity update
v m n f ( m ) = v m 1 n f ( m 1 ) + Δ v s f ( m ) n f + Δ v c o r / g ( m ) n f
where Δ v s f ( m ) n f and Δ v c o r / g ( m ) n f are the specific force velocity increment of the navigation system and the velocity increment of the detrimental acceleration in the time period T m , respectively. Δ v r o t ( m ) b f ( m 1 ) is the rotational error compensation amount of the speed, which is caused by the rotation change of the specific force direction in space during the solution time period; Δ v s c u l ( m ) b f ( m 1 ) is the paddling error compensation amount.
Δ v c o r / g ( m ) n f [ 2 ω i e ( m 1 / 2 ) n f + ω e n ( m 1 / 2 ) n f ] × v m 1 / 2 n f T m + g m 1 / 2 n f T m , Δ v s f ( m ) n f = [ I T m 2 ( ω i n ( m 1 / 2 ) n f × ) ] C b f ( m 1 ) n f ( m 1 ) Δ v m + C b f ( m 1 ) n f ( m 1 ) ( Δ v r o t ( m ) b f ( m 1 ) + Δ v s c u l ( m ) b f ( m 1 ) ) , Δ v r o t ( m ) b f ( m 1 ) = 1 2 Δ θ m × Δ v m ,   Δ v s c u l ( m ) b f ( m 1 ) = 2 3 ( Δ θ m 1 × Δ v m 2 + Δ v m 1 × Δ θ m 2 ) .
where Δ θ m = Δ θ m 1 + Δ θ m 2 , Δ v m = Δ v m 1 + Δ v m 2 , Δ v m 1 denotes the velocity increment in the time interval [ t m 1 t ( m 1 / 2 ) ] , and Δ v m 2 denotes the velocity increment in the time interval [ t ( m 1 / 2 ) t m ] .
(3)
Position update
p m = p m 1 + M p v ( m 1 / 2 ) ( v m 1 n f ( m 1 ) + v m n f ( m ) ) T 2
where M p v = [ 0 1 / R M h 0 sec L / R N h 0 0 0 0 1 ] , M p v ( m 1 / 2 ) can be obtained by linear extrapolation [28].

3.3.2. Backtracking Navigation

To obtain the inertial navigation attitude, velocity and position update equations for the backtracking solution, Equations (3)–(7) are transformed by shifting the terms, by taking the reverse of the angular velocity of ground rotation, the angular velocity of rotation of the navigation CS relative to the earth CS, the acceleration of gravity, the angular increment of the gyro output and the velocity increment of the angular velocity meter output to obtain:
(1)
Attitude update during backtracking
C b r ( q ) n r ( q ) = C n r ( q 1 ) n r ( q ) C b r ( q 1 ) n r ( q 1 ) C b r ( q ) b r ( q 1 )
C n r ( q 1 ) n r ( q ) = M R V ( T q ω i n ( q 1 / 2 ) n r ) , C b r ( q ) b r ( q 1 ) = M R V ( ϕ i b ( q ) b r )
ϕ i b ( q ) b r = ( Δ θ q 2 + Δ θ q 1 ) + 2 3 Δ θ q 1 × Δ θ q 2
where C b r ( q 1 ) n r ( q 1 ) and C b r ( q ) n r ( q ) represent the SINS attitude matrix at t q 1 and t q ( q = 1 , 2 , 3 , n q ), respectively. For the backtracking navigation update period [ t q 1 , t q ] , T q = t q t q 1 . t ( q 1 / 2 ) = ( t q + t q 1 ) / 2 denotes the midpoint moment. ω i n ( q 1 / 2 ) n r represents the rotation of the n system with respect to the i system at the moment of the midpoint in the backtracking navigation. The subscript ( q 1 / 2 ) indicates that they are the variables at the midpoint moment. Δ θ q 1 denotes the angular increment in the time interval [ t q 1 t ( q 1 / 2 ) ] , and Δ θ q 2 denotes the angular increment in the time interval [ t ( q 1 / 2 ) t q ] .
(2)
Velocity update during backtracking
v q n r ( q ) = v q 1 n r ( q 1 ) + Δ v s f ( q ) n r + Δ v c o r / g ( q ) n r
where Δ v s f ( m ) n r and Δ v c o r / g ( m ) n r are the specific force velocity increment of the backtracking navigation system and the velocity increment of the detrimental acceleration in the time period T q , respectively. Δ v r o t ( q ) b r ( q 1 ) is the rotational error compensation amount of the speed, which is caused by the rotation change of the specific force direction in space during the solution time period; Δ v s c u l ( q ) b r ( q 1 ) is the paddling error compensation amount in the backtracking navigation system.
Δ v c o r / g ( q ) n r ( 2 ω i e ( q 1 / 2 ) n r + ω e n ( q 1 / 2 ) n r ) × v q 1 / 2 n r T q + g q 1 / 2 n r T q , Δ v s f ( q ) n r = [ I + T q 2 ( ω i n ( q 1 / 2 ) n r × ) ] C b r ( q ) n r ( q ) Δ v q + C b r ( q ) n r ( q ) ( Δ v r o t ( q ) b r ( q 1 ) + Δ v s c u l ( q ) b r ( q 1 ) ) , Δ v r o t ( q ) b r ( q 1 ) = 1 2 Δ v q × Δ θ q , Δ v s c u l ( q ) b r ( q 1 ) = 2 3 ( Δ v q 1 × Δ θ q 2 + Δ θ q 1 × Δ v q 2 ) .
where Δ θ q = Δ θ q 1 + Δ θ q 2 , Δ v q = Δ v q 1 + Δ v q 2 , Δ v q 1 denotes the velocity increment in the time interval [ t q 1 t ( q 1 / 2 ) ] , and Δ v q 2 denotes the velocity increment in the time interval [ t ( q 1 / 2 ) t q ] .
(3)
Position update during backtracking
p q = p q 1 M p v ( q 1 / 2 ) ( v q 1 n r ( q 1 ) + v q n r ( q ) ) T 2
Through the above derivation, it can be seen that the attitude, velocity and position can at the same time be basically consistent during the forward and backward navigation process, the same piece of navigation data can be solved repeatedly and the known information at the starting point can be used for correction to improve the alignment accuracy.

3.4. Reverse Kalman Filtering

The error equation of the BN system is deduced according to the inertial navigation error equation [25].
The actual calculated attitude matrix C ˜ b r n r has a small angle of attitude error angle with the ideal attitude matrix C b r n r , then C ˜ b r n r can be expressed as
C ˜ b r n r = [ I ( ϕ × ) ] C b r n r
The calculated attitude matrix C ˜ b r n r is generally solved by the following differential equation:
C ˜ ˙ b r n r = C ˜ b r n r ( ω ˜ i b b r × ) ( ω ˜ i n n r × ) C ˜ b r n r
From the above, it is clear that the BN process requires the inversion of the gyro output angular velocity ω i b b f and the rotation angular velocity ω i n n f of the navigation system, i.e.,
ω ˜ i b b r = ω i b b f δ ω i b b f
ω ˜ i n n f = ω i n n f δ ω i n n f
Substituting the differentiation of Equations (15), (16) and (13) into Equation (14), we get
( ϕ ˙ × ) C b r n r + ( I ( ϕ × ) ) [ C b r n r ( ω i b b f × ) ( ω i n n f × ) C b r n r ] = ( I ϕ × ) C b r n r [ ( ω i b b f δ ω i b b f ) × ] [ ( ω i n n f δ ω i n n f ) × ] ( I ϕ × ) C b r n r
Multiplying both sides of Equation (17) simultaneously right by C n r b r and omitting the second-order small quantity,
( ϕ ˙ × ) = [ ( ω i n n f × ) ( ϕ × ) ( ϕ × ) ( ω i n n f × ) ] ( δ ω i n n f × ) + C b r n r ( δ ω i b b f × ) C n r b r
The inertial navigation attitude error equation for BN is obtained by simplifying Equation (18),
ϕ ˙ = ( ω i n n f × ϕ ) δ ω i n n f + δ ω i b b f
According to the velocity update equation for BN, the differential equation of inertial navigation velocity for BN can be expressed as
ν ˙ n r = C b r n r ( f n f ) [ 2 ( ω i e n f ) + ( ω e n n f ) ] × ν n r + ( g n f )
The actual calculated velocity ν ˜ n r is obtained by Equation (21).
ν ˜ ˙ n r = C ˜ b r n r f ˜ n r ( 2 ω ˜ i e n r + ω ˜ e n n r ) × v ˜ n r + g ˜ n r
where f ˜ b r = f b f δ f b f , ω ˜ i e n r = ω i e n f δ ω i e n f , ω ˜ e n n r = ω e n n f δ ω e n n f , g ˜ n r = g n f δ g n f .
The velocity error δ ν n r can be obtained by subtracting Equations (21) and (20)
δ ν ˙ n r = [ C ˜ b r n r f ˜ n r C b r n r ( f n f ) ] { ( 2 ω ˜ i e n r + ω ˜ e n n r ) × v ˜ n r [ 2 ( ω i e n f ) + ( ω e n n f ) ] × v ˜ n r } + [ g ˜ n r ( g n f ) ]
Substituting Equation (13) into Equation (22) and expanding it gives:
δ ν ˙ n r = [ ( I ϕ × ) C b r n r ( f b f δ f b f ) C b r n r ( f b f ) ] { [ 2 ( ω i e n f δ ω i e n f ) + ( ω e n n f δ ω e n n f ) ] × ( v n r + δ v n r ) [ 2 ( ω i e n f ) + ( ω e n n f ) ] × v ˜ n r } δ g n f
The Equation (23) is simplified by omitting the second-order small quantity error and the gravity error δ g n f , and the velocity error equation for the BN is obtained as follows:
δ ν ˙ n r = C b r n r f b f × ϕ v n r × ( 2 δ ω i e n f + δ ω e n n f ) + ( 2 ω i e n f + ω e n n f ) × δ v n r δ f b f
The position update differential equation for BN can be written as:
L ˙ = ν N n / ( R M + h )
λ ˙ = ν E n sec L / ( R N + h )
h ˙ = ν U
The position error equation for the BN is obtained from the position update differential equation Equations (25)–(27) as follows:
δ L ˙ = δ ν N / ( R M + h ) + ν N δ h / ( R M + h ) 2
δ λ ˙ = δ ν E sec L / ( R N + h ) ν E δ L sec L tan L / ( R N + h ) + ν E δ h sec L / ( R N + h ) 2
δ h ˙ = δ ν U
Integrating the error equation of the BN derived above, the corresponding systematic error equation of the BN system is as follows:
{ ϕ ˙ = ω i n n f × ϕ δ ω i n n f + C b f n f ε b f δ v ˙ = f n f × ϕ v n r × ( 2 δ ω i e n f + δ ω e n n f ) + ( 2 ω i e n f + ω e n n f ) × δ v C b f n f b f δ p ˙ = M p v δ v + M p p δ p ε ˙ n f = [ 0 0 0 ] T ˙ n f = [ 0 0 0 ] T
Define:
M 1 = [ 0 0 0 ω i e sin L 0 0 ω i e cos L 0 0 ] , M 2 = [ 0 1 R M h 0 1 R N h 0 0 tan L R N h 0 0 ] , M 3 = [ 0 0 v N n R M h 2 0 0 v E n R N h 2 v E n sec 2 L R N h 0 v E n tan L R N h 2 ] ,   δ ω i e n f = M 1 δ p ,   δ ω e n n f = M 2 δ v + M 3 δ p ,   M a a = ( ω i n n f × ) ,   M a p = M 1 M 3 ,   M v a = ( f n f × ) , M v v = v n r × M 2 + ( ( 2 ω i e n f + ω e n n f ) × ) ,   M v p = v n r × ( 2 M 1 + M 3 ) , M p v = [ 0 1 R M h 0 sec L R N h 0 0 0 0 1 ] ,   M p p = [ 0 0 v N n R M h 2 v E n sec L tan L R N h 2 0 v E n sec L R N h 2 0 0 0 ]
where ϕ , δ ν and δ p represent the attitude error, velocity error and position error of the SINS in the three directions of East, North and Up, respectively; ε and are the gyroscopic drift and the accelerometer zero error, respectively.
According to Equation (31), the reverse Kalman filtering model of the navigation system can be deduced as:
X ˙ r = F r X + G r W Z r = H r X + V
where X r is a 14-dimensional state variable; Z r is the observation vector; F r and H r represent the system matrix and observation matrix, respectively; G r is the noise allocation matrix; W and V represent the system noise and observation noise, respectively.
State variables:
X = [ ϕ T δ v T δ p T ( ε b ) T ( b ) T ] T
System matrix in reverse Kalman filtering:
F r = [ M a a M a v M a p C b n 0 3 × 3 M v a M v v M v p 0 3 × 3 C b n 0 3 × 3 M p v M p p 0 3 × 3 0 3 × 3 0 6 × 3 0 6 × 3 0 6 × 3 0 6 × 3 0 6 × 3 ]
Observation vector:
Z r = v S I N S
where v S I N S denotes the velocity component obtained from the SINS solution.
Observation Matrix:
H r = [ 0 3 × 3 I 3 × 3 0 3 × 9 ]
According to the above equations, more accurate attitude information can be obtained by retrospectively solving the system state through the backtracking Kalman filtering algorithm.
Remark 1.
The continuous-time reverse Kalman filter model derived above is discretized, and then the corresponding terms are replaced by Δ θ q and Δ v q .
Δ θ q = t q 1 t q ω i b b r ( t ) d t ,   Δ v q = t q 1 t q f s f b r ( t ) d t .
For a launch vehicle, the changes of ω i b b r and f s f b r are small in a short period of time [ t q 1 , t q ] , so it is generally considered that the quadratic functions of Δ θ q and Δ v q are slow quantities of time. The calculation error caused is generally small, and the relatively simple gradient integration method can be used to obtain
Δ θ q = t q 1 t q ω i b b r ( t ) d t ( ω i b b r ( t q 1 ) + ω i b b r ( t q ) ) T q 2 ω i b ( q 1 / 2 ) b r T q   Δ v q = t q 1 t q f s f b r ( t ) d t ( f s f b r ( t q 1 ) + f s f b r ( t q ) ) T q 2 f s f ( q 1 / 2 ) b r T q
where ω i b ( q 1 / 2 ) b r denotes the output angular velocity of the gyroscope at the moment of the midpoint t ( q 1 / 2 ) ; f s f ( q 1 / 2 ) b r denotes the output specific force of the accelerometer at the moment of the midpoint t ( q 1 / 2 ) .

3.5. Cycle-Index Control Function

For incremental SINS, repeated forward and backtracking analysis to calculate the stored data is beneficial to improving the navigation accuracy. In this process, in order to simplify the calculation, some second-order small quantities are often omitted, or small angle approximation is adopted, so that the error of each iteration is accumulated. When the number of forward and backtracking is small, the effect of error on the results is not obvious. When the number of forward and backtracking is large and reaches a certain number of times, the error is amplified to the same multiple as the number of times due to repeated calculations, resulting in oscillation convergence or even the drift of the result not converging.
For the case where the swaying amplitude is small, the backtracking process approximates the strict reverse derivation process, and the solution is more accurate. Multiple iterations of forward and backtracking solutions can guarantee the accuracy and avoid the accumulation of calculation errors. At this point, the cycle-index control function is established as follows:
q = [ α t s i m t c o a r ]
where the whole total sampling time is recorded as t s i m , the first backtracking time is recorded as t c o a r , and α is the frequency control coefficient.

3.6. Attitude Error Compensation

During the self-alignment process of the LV launched on the offshore platform, when returning to the starting point of the BN, the velocity and attitude angle change due to the influence of waves, etc., so attitude error compensation is required.
The first attitude error compensation: after obtaining the relatively attitude matrix at moment t 1 , the attitude update is solved with its initial value using the SINS data stored after moment t 1 . It is updated until the latest sampled data moment t 2 , thereby obtaining the attitude matrix estimated at moment t 2 .
The second attitude error compensation: the cyclic solution starts at the moment t 2 , while the sampled data after t 2 are stored. After the last forward navigation to obtain the exact matrix at time t 2 , the SINS data stored after moment t 2 are used as the initial value for attitude update solving. Until the IA ends, the latest attitude matrix is obtained.

4. Experiments and Results Analysis

Remark 2.
The self-alignment of the LV launched on OP is obviously affected by the line vibration. The oscillation in the convergence section of the estimated curve reflects the characteristics of the motion and swaying of the LV. For the influence of line vibration disturbance, considering that the period of the gravitational acceleration vector is very small with inertia, low-pass filtering and other methods can be used, but this is not within the scope of this study.
Remark 3.
The orientation attitude of the LV’s self-alignment is the roll angle rather than the yaw angle.
Remark 4.
The form of attitude matrix in the table is: [pitch, yaw, roll].
The simulation process simulates the self-alignment process of the LV in the mooring condition. In the Xisha and towards the equatorial sea area, the swaying of the LV and the launch platform disturbed by the waves is approximately sinusoidal with equal amplitude oscillation in sea state level 4, the swaying amplitude does not exceed 2.5 ° and the swing period is 4 s 14 s .
The hardware environment for simulation was an Intel(R) Core(TM) i5-8400 CPU with 2.80 GHz, 8.00 GB RAM and Windows 10 operating system.
The inertial sensor error parameters [29] were set as shown in Table 1. The data output rate of the inertial sensor was set to 200 Hz.
The swaying parameters [30] are shown in Table 2.
In the mooring condition, in addition to swaying, there are also lateral, longitudinal and vertical linear velocities. The linear velocity parameters [30] are shown in Table 3.
The installation height of the SINS in the LV does not exceed 60 m. Combining the swaying of the offshore launch platform with the linear velocity, the trajectory of the LV in the relative coordinate system is shown in Figure 4 and it is bounded by the scope. It has a displacement of approximately one meter, which is characteristic of the engineering application environment.
The simulation scheme is as follows:
Scheme A: The conventional method uses the first 300 s of data for ASCA, and then uses the next 360 s of data for conventional Kalman filtering FA.
Scheme B: ASCA uses the first 300 s of data, the first backtracking FA uses the first 300 s of data, the first forward FA uses the first 330 s of data, repeated backtracking and forward FA uses the first 330 s of data; and the final forward FA uses the first 360 s of data.
The time taken for scheme A and B is shown in Figure 5.
The time taken for scheme A is 660 s, while Scheme B takes 360 s. Therefore, the proposed scheme shortens the time by 45.5% compared with Scheme A.

4.1. Estimation of the Misalignment Angle at the First Cycle

Using the strategy proposed in this paper, a stochastic simulation was performed. Figure 6 shows the estimated curves of the misalignment angles during the first backtracking navigation FA and the first forward navigation FA. The arrow in Figure 6 indicates the estimated direction.
After using the reverse Kalman filtering back to the IM and correcting the SINS with the estimated results, most of the attitude misalignment angles are corrected, and the filter curve converges during the forward alignment process, which can guarantee the alignment accuracy.

4.2. Analysis of the Last Forward Navigation Alignment Results

The CA with an alignment time of 300 s and the FA with an alignment time of 360 s were studied to simulate the conventional self-alignment based on scheme A and the self-alignment strategy based on scheme B, and the results of attitude error obtained for scheme A and scheme B are shown in Figure 7.
A stochastic simulation was performed. In the conventional approach, the roll angle error enters the error band [ 1 , 1 0 ] after completing the alignment, while the roll angle error enters the error band [ 0 , 3 ] in the proposed method. Since the proposed method utilizes the SINS data from the CA stage in the FA, more information is obtained at the same time. The alignment time of the proposed strategy depends on the convergence time of the FA filter, and the CA no longer occupies the alignment time; the SINS data of the CA stage are used twice with less storage and calculation, which not only increases the quantity of data but also effectively shortens the alignment time.
The convergence segment of the simulation results from 300 s to 360 s was selected for analysis. The statistics in Table 4 are for the convergence segment.
As can be seen from Table 4, the convergence quality of the proposed strategy is significantly better than that of the conventional method. The proposed strategy makes use of the SINS data in the CA stage from the beginning and obtains more information in a limited time with better results; with the increase in the frequency of backtracking, there is no obvious adjustment in the alignment process of the method. Through analysis, the algorithm can achieve accuracy after one or two backtracking processes.

4.3. Simulation Analysis of Multiple Reverse Navigation

Under the above simulation conditions, four of the reverse processes were selected and the alignment error curves are shown in Figure 8. Clearly, the orientation attitude angle is continuously corrected as the number of reversals increases.

4.4. Monte Carlo Simulation Results and Discussion

The statistics of the results of 50 simulations under the alignment conditions of scheme A and scheme B are shown in Figure 9.
The statistical characterization of the converged three-axis attitude angle errors was carried out, and the corresponding characterization is shown in Table 5.
As can be seen from Table 5, the alignment accuracy of scheme B is obviously better than that of scheme A, its attitude angle error is less than 3.9′ and the average value is within 0.82′. The standard deviations of pitch error, roll error and yaw error are 0.0853′, 1.4095′ and 0.0123′, respectively. It can be seen that the stability of the convergence section of the attitude angle error in the horizontal direction is improved, and the stability of the convergence section of the roll angle error decreases.

4.5. The Effect of CA Time on Alignment Accuracy of the Proposed Method

At present, there are few studies on the BN that focus on the effect of CA time in alignment accuracy. Based on the proposed method in this paper, 50 Monte Carlo simulations were conducted using the CA time of 120, 240, 300 and 330 s. The statistical analysis of the alignment accuracy results is shown in Table 6.
In general, the ASCA algorithm has a relatively weak anti-interference ability within a short period of time, and if the CA time is extended appropriately, good CA results can often be achieved. However, it can be seen from Table 6 that appropriately extending the CA time and keeping the FA time unchanged does not necessarily improve the final alignment accuracy, which is related to the accumulation of errors during the backtracking.

4.6. Validation Test

To compare the proposed strategy with conventional optical aiming method, the laser IMU is mounted on a six-degree-of-freedom turntable, and the relationship between the IMU orientation and the frame angle of the turntable is determined using the prism on the IMU to establish the reference orientation. The update frequency of IMU’s raw data is set to 100 Hz. The bias stability of gyroscopes and accelerometers are better than 0.01°/h and 20 µg, respectively. Following the laboratory confidentiality requirements, we will not show the physical picture.
Different attitude commands were executed on the turntable to simulate the shaking conditions during the launch of an LV, and the alignment results obtained are shown in Table 7.
As can be seen from Table 7, the experimental results under various conditions are basically consistent. The proposed alignment strategy can track the sway of the SINS well, and the alignment accuracy is within 3′.

5. Conclusions and Outlook

In this paper, a rapid self-alignment strategy is proposed for the characteristics of the maritime motion environment, which can complete the initial alignment work without setting up the aiming equipment and can solve the problem that the traditional alignment algorithm is not applicable in the maritime initial alignment system model. The strategy can effectively improve the initial alignment accuracy and meet the requirements of a launch vehicle on an offshore launching platform. The results show that the proposed fully autonomous alignment algorithm can be used to achieve the initial alignment of the SINS before the launch of the LV. To meet the demand of low cost and fast launch of LVs, the fully autonomous alignment technology will be gradually applied to LVs to replace the optical aiming system.

Author Contributions

Conceptualization, R.M. and T.Z.; data curation, R.M. and T.Z.; formal analysis, R.M. and T.Z.; funding acquisition, R.M.; investigation, T.Z.; methodology, T.Z. and S.L.; project administration, R.M.; resources, T.Z. and S.L.; software, T.Z. and S.L.; validation, R.M., T.Z. and S.L.; writing—original draft, T.Z.; writing—review and editing, R.M., T.Z. and S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fourth Batch of Pre-research on Manned Spaceflight, grant number 18123060201.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Most of the data and models generated or used during the research appear in this publication.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
SINSStrapdown Inertial Navigation System
IMUInertial Measurement Unit
RLGRing Laser Gyroscope
FOGFiber-Optic Gyroscope
MCMonte Carlo
LVLaunch Vehicle
ASCAAnti-Swaying Coarse Alignment
BNBacktracking Navigation
RKFReverse Kalman Filtering
CICFCycle-Index Control Function
OLPOffshore launching Platform
FAFine Alignment
CACoarse Alignment
IMInitial Moment
CSCoordinate system
IAInitial Alignment

References

  1. Fu, M.Y.; Zheng, X.; Deng, Z.H. Transfer Alignment Theory and Applications, 1st ed.; Science Press: Beijing, China, 2016; pp. 1–2. [Google Scholar]
  2. Zhang, H.L.; Zhou, Y.L.; Huang, T.C.; Wang, L. Anti-Swaying Self-Alignment Method Based on Observation Vector Inertia for SINS under Latitude Uncertainty. Measurement 2021, 177, 109250. [Google Scholar] [CrossRef]
  3. Wang, S.Y.; Zhang, Q.; Guo, T. Research on the Application of Adaptive Transfer Alignment Technology in the Initial Alignment on an Offshore Launching Platform. China CIO News. 2021. Available online: https://qikan.cqvip.com/Qikan/Article/ReadIndex?id=7105510250&info=eUjW%2fnPsSZYKMRAjEt0qYn1ByP74B25tKphQ%2f9urBA6BiRRTmRl1gA%3d%3d (accessed on 1 October 2022).
  4. Wen, Z.Y.; Yang, G.L.; Cai, Q.Z.; Sun, Y.D. Odometer Aided SINS in-Motion Alignment Method Based on Backtracking Scheme for Large Misalignment Angles. IEEE Access 2020, 8, 7937–7948. [Google Scholar] [CrossRef]
  5. Silva, F.O.; Hemerly, E.M.; Filho, C.L. Error Analysis of Analytical Coarse Alignment Formulations for Stationary SINS. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 1777–1796. [Google Scholar] [CrossRef]
  6. Qin, Y.Y.; Yan, G.M.; Gu, D.Q.; Zheng, J.B.; Qin, Y.; Yan, G.; Gu, D. A Clever Way of SINS Coarse Alignment despite Rocking Ship. J. Northwest. Univ. 2005, 23, 681–683. [Google Scholar]
  7. Bai, L.; Qin, Y.Y.; Zhu, X.Y. Validation of Self-alignment Algorithm of SINS in a Swaying Vehicle. Fire Control Command Control 2010, 35, 164–167. [Google Scholar]
  8. Liu, X.X.; Liu, X.J.; Song, Q.; Yang, Y.; Liu, Y.T.; Wang, L.H. A novel self-alignment method for SINS based on three vectors of gravitational apparent motion in inertial frame. Measurement 2015, 62, 47–62. [Google Scholar] [CrossRef]
  9. Yan, G.M.; Qin, Y.Y.; Wei, Y.X.; Zhang, L.C.; Yan, W.S. New Initial Alignment Algorithm for SINS on Moving Base. Syst. Eng. Electron. 2009, 31, 634–637. [Google Scholar]
  10. Yan, G.M.; Weng, J.; Bai, L. Initial in-Movement Alignment and Position Determination Based on Inertial Reference Frame. Syst. Eng. Electron. 2011, 33, 618–623. [Google Scholar]
  11. Wu, M.P.; Wu, Y.X.; Hu, X.P.; Hu, D.W. Optimization-based Alignment for Inertial Navigation Systems: Theory and Algorithm. Aerosp. Sci. Technol. 2011, 15, 1–17. [Google Scholar] [CrossRef]
  12. Guo, Y.S.; Fu, M.Y.; Deng, Z.H.; Zhuang, G.C.; Deng, L.; Automation, S.O. Application of Quaternion Estimator Algorithm Dedicated on Alignment of Swaying and Moving Carrier. J. Chin. Inertial Technol. 2017, 25, 182–185. [Google Scholar]
  13. Chang, Z.J.; Zhang, Z.L.; Zhou, Z.F.; Chen, H.; Zhao, J.Y. Alignment Method of Rotating SINS with Swinging Base. Syst. Eng. Electron. 2019, 41, 2087–2091. [Google Scholar]
  14. Groves, P.D. Principles of GNSS, Inertial and Multisensor Integrated Navigation Systems, 2nd ed.; Artech House: Norwood, MA, USA, 2013. [Google Scholar]
  15. Wang, H.L.; Li, G.Z.; Shi, Z.Y. Overview of Initial Alignment Method for Strap down Inertial Navigation System. In Proceedings of the Advances in Materials, Machinery, Electrical Engineering (AMMEE 2017), Tianjin, China, 10–11 June 2017; Volume 114, pp. 149–154. [Google Scholar]
  16. Tsukerman, A.; Klein, I. Analytic Evaluation of Fine Alignment for Velocity Aided INS. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 376–384. [Google Scholar] [CrossRef]
  17. Chang, L.B.; Qin, F.; Li, A. A Novel Backtracking Scheme for Attitude Determination-Based Initial Alignment. IEEE Trans. Autom. Sci. Eng. 2015, 12, 384–390. [Google Scholar] [CrossRef]
  18. Chang, L.B.; Hu, B.Q.; Li, Y. Backtracking Integration for Fast Attitude Determination-based Initial Alignment. IEEE Trans. Instrum. Meas. 2015, 64, 795–803. [Google Scholar] [CrossRef]
  19. Gao, W.; Ben, Y.Y.; Zhang, X.; Li, Q.; Yu, F. Rapid Fine Strapdown INS Alignment Method under Marine Mooring Condition. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 2887–2896. [Google Scholar] [CrossRef]
  20. Fu, Q.; Li, S.; Liu, Y.; Wu, F. Information-Reusing Alignment Technology for Rotating Inertial Navigation System. Aerosp. Sci. Technol. 2020, 99, 105747. [Google Scholar] [CrossRef]
  21. Li, W.L.; Wu, W.Q.; Wang, J.L.; Wu, M.P. A Novel Backtracking Navigation Scheme for Autonomous Underwater Vehicles. Measurement 2014, 47, 496–504. [Google Scholar] [CrossRef]
  22. Yan, G.M.; Yan, W.S.; Xu, D.M. Inverse navigation algorithm and its application to initial alignment of strapdown moving base. In Proceedings of the 27th China Control Conference, Kunming, China, 16–18 July 2008. [Google Scholar]
  23. Liang, L.Y.; Liu, Q.J.; Yang, Y. The SINS Self-Alignment Technology Based on Reverse Navigation Algorithm. Aerosp. Control 2015, 33, 49–53. [Google Scholar]
  24. Wang, D.S.; He, G.Y.; Jiang, X.H. In-motion Alignment Scheme Based on Reverse Kalman Filter. J. Chin. Inertial Technol. 2020, 28, 721–728. [Google Scholar]
  25. Xu, X.; Xu, D.C. A Forward-Forward Backtracking Alignment Method for SINS on Swaying Base. J. Astronaut. 2019, 40, 965–971. [Google Scholar]
  26. Yan, G.M.; Weng, J. Strapdown Inertial Navigation Algorithm and Principle of Integrated Navigation, 1st ed.; Northwestern Polytechnic University Press: Xi’an, China, 2019; pp. 78–95. [Google Scholar]
  27. Jiao, Y.B.; Li, J.; Ma, X.H.; Feng, K.Q.; Guo, X.T.; Wei, X.K.; Feng, Y.J.; Zhang, C.M.; Wang, J.Q. An Anti-Turbulence Self-Alignment Method for SINS under Unknown Latitude Conditions. Sensors 2022, 22, 4686. [Google Scholar] [CrossRef] [PubMed]
  28. Gu, D.Q.; Qin, Y.Y. Research on the Optimal Two-sample Scrolling Compensation Algorithm for SINS Position Updating. J. Astronaut. 2005, 26, 562–566. [Google Scholar]
  29. Yan, G.M.; Li, S.H.; Qin, Y.Y. Inertial Instrument Testing and Data Analysis, 1st ed.; National Defense Industry Press: Beijing, China, 2015; pp. 148–195. [Google Scholar]
  30. Qin, Y.Y. Inertial Navigation, 2nd ed.; Science Press: Beijing, China, 2015; pp. 316–340. [Google Scholar]
Figure 1. Fundamentals of initial alignment of launch vehicles.
Figure 1. Fundamentals of initial alignment of launch vehicles.
Sensors 23 00339 g001
Figure 2. Launch CS and IMU CS of LV [26]. (a) Launch CS and IMU CS; (b) Euler angle definition.
Figure 2. Launch CS and IMU CS of LV [26]. (a) Launch CS and IMU CS; (b) Euler angle definition.
Sensors 23 00339 g002
Figure 3. The self-alignment strategy of LV launched on the offshore platform.
Figure 3. The self-alignment strategy of LV launched on the offshore platform.
Sensors 23 00339 g003
Figure 4. The trajectory of the LV swaying with the offshore launch platform.
Figure 4. The trajectory of the LV swaying with the offshore launch platform.
Sensors 23 00339 g004
Figure 5. The time taken for schemes A and B.
Figure 5. The time taken for schemes A and B.
Sensors 23 00339 g005
Figure 6. Estimated curves of the misalignment angles during the first backtracking navigation FA and the first forward navigation FA.
Figure 6. Estimated curves of the misalignment angles during the first backtracking navigation FA and the first forward navigation FA.
Sensors 23 00339 g006
Figure 7. Three-axis attitude error curves for the FA of scheme A and the last FA of scheme B.
Figure 7. Three-axis attitude error curves for the FA of scheme A and the last FA of scheme B.
Sensors 23 00339 g007
Figure 8. Curves of attitude errors.
Figure 8. Curves of attitude errors.
Sensors 23 00339 g008
Figure 9. Alignment result statistics for 50 simulations: (a) Scheme A; (b) Scheme B.
Figure 9. Alignment result statistics for 50 simulations: (a) Scheme A; (b) Scheme B.
Sensors 23 00339 g009
Table 1. Sensor errors settings.
Table 1. Sensor errors settings.
Error TermsGyroscopeAccelerometer
3-Axis (1 σ ) 3-Axis (1 σ )
Scale factor
repeatability deviation
0.000050.000037
Zero bias stability deviation 0.01 ° / h 20   μ g
Random walk coefficient 0.005 ° / h 1.4   μ g / H z
Table 2. Swaying parameter settings.
Table 2. Swaying parameter settings.
PitchRollYaw
Swing amplitude (°)1.752.51.25
Swing period (s)567
Table 3. Linear velocity parameter settings.
Table 3. Linear velocity parameter settings.
LaterallongitudinalVertical
Amplitude (m/s)0.20.030.02
Period (s)786
Table 4. Simulation results.
Table 4. Simulation results.
SchemeConvergence Analysis
MEAN ( ) STD ( )/(1 σ ) MAX ( ) MIN ( )
A[0.4063,
0.1269,
5.8172]
[0.4846,
0.9554,
2.2690]
[1.2277,
1.6750,
10.2059]
[−0.5708,
−1.4288,
1.0795]
B[−0.0415,
0.0487,
1.2252]
[0.4695,
0.3474,
0.6638]
[0.6912,
0.5822,
2.2717]
[−0.7575,
−0.4713,
0.1863]
Table 5. Statistics for simulation result.
Table 5. Statistics for simulation result.
SchemeStatistics
MEAN ( ) STD ( )/(1 σ ) MAX ( ) MIN ( )
A[0.5065,
−0.0448,
2.8971]
[0.0681,
0.0424,
2.2384]
[0.5677,
0.1502,
7.0392]
[0.3206,
−0.0448,
−1.0826]
B[0.2644,
−0.0169,
−0.8103]
[0.0853,
0.0123,
1.4095]
[0.3048,
0.0184,
3.5274]
[−0.0229,
−0.0395,
−3.8959]
C[−0.0364,
0.0928,
8.7820]
[0.0097,
0.0114,
3.9653]
[−0.0171,
0.1188,
17.0969]
[−0.0687
−0.0613,
−1.6983]
Table 6. Statistics for simulation results.
Table 6. Statistics for simulation results.
t C o a r s e Statistics
MEAN ( ) STD ( )/(1 σ ) MAX ( ) MIN ( )
120 s[0.1516,
−0.0169,
−0.1848]
[0.0633,
0.0145,
2.5233]
[0.2716,
0.0064,
8.5242]
[−0.0067,
−0.0545,
−6.3868]
240 s[0.2075,
−0.0193,
2.1955]
[0.1000,
0.0127,
3.8766]
[0.3662,
0.0087,
9.8579]
[0.0145,
−0.0462,
−5.4194]
300 s[0.2644,
−0.0169,
−0.8103]
[0.0853,
0.0123,
1.4095]
[0.3048,
0.0184,
3.5274]
[−0.0229,
−0.0395,
−3.8959]
330 s[0.1752,
−0.0156,
0.3078]
[0.0728,
0.0131,
2.5647]
[0.3086,
0.0148,
6.5900]
[−0.0090,
−0.0457,
−7.2000]
Table 7. Results of highly precise alignment of laser gyroscope inertial measurement unit dynamic experiment.
Table 7. Results of highly precise alignment of laser gyroscope inertial measurement unit dynamic experiment.
No.Turntable Position (°)Environmental ConditionsOptical Aiming Azimuth (°)Azimuth Error of the Proposed Strategy ( )
Inner FrameMiddle FrameOuter FrameAngular MovementLinear Velocity
1900.03156.7Amplitude: [0.35° 0.5° 0.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
8.986−0.520
2900.03146.7Amplitude: [0.35° 0.5° 0.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
−1.9730.880
3901.03156.7Amplitude: [0.35° 0.5° 0.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
8.9861.790
4901.03146.7Amplitude: [0.35° 0.5° 0.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
−1.973−1.360
5900.03156.7Amplitude: [1.75° 2.5° 1.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
9.772−0.970
6900.03146.7Amplitude: [1.75° 2.5° 1.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
−2.1381.990
7901.03156.7Amplitude: [1.75° 2.5° 1.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
9.7722.170
8901.03146.7Amplitude: [1.75° 2.5° 1.25°]
Period: [5 s 6 s 7 s]
Amplitude: [0.2 0.03 0.02]
Period: [7 s 8 s 6 s]
−2.138−1.780
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

Mu, R.; Zhang, T.; Li, S. A Rapid Self-Alignment Strategy for a Launch Vehicle on an Offshore Launching Platform. Sensors 2023, 23, 339. https://doi.org/10.3390/s23010339

AMA Style

Mu R, Zhang T, Li S. A Rapid Self-Alignment Strategy for a Launch Vehicle on an Offshore Launching Platform. Sensors. 2023; 23(1):339. https://doi.org/10.3390/s23010339

Chicago/Turabian Style

Mu, Rongjun, Tengfei Zhang, and Shoupeng Li. 2023. "A Rapid Self-Alignment Strategy for a Launch Vehicle on an Offshore Launching Platform" Sensors 23, no. 1: 339. https://doi.org/10.3390/s23010339

APA Style

Mu, R., Zhang, T., & Li, S. (2023). A Rapid Self-Alignment Strategy for a Launch Vehicle on an Offshore Launching Platform. Sensors, 23(1), 339. https://doi.org/10.3390/s23010339

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