Qi Wang , Changsong Yang Shaoen Wu and Yuxiang Wang
Abstract: Transfer alignment is used to initialize SINS (Strapdown Inertial Navigation System) in motion. Lever-arm effect compensation is studied existing in an AUV(Autonomous Underwater Vehicle) before launched from the mother ship. The AUV is equipped with SINS, Doppler Velocity Log, depth sensor and other navigation sensors.The lever arm will cause large error on the transfer alignment between master inertial navigation system and slave inertial navigation system, especially in big ship situations.This paper presents a novel method that can effectively estimate and compensate the flexural lever arm between the main inertial navigation system mounted on the mother ship and the slave inertial navigation system equipped on the AUV. The nonlinear measurement equation of angular rate is derived based on three successive rotations of the body frame of the master inertial navigation system. Nonlinear filter is utilized as the nonlinear estimator for its capability of non-linear approximation. Observability analysis was conducted on the SINS state vector based on singular value decomposition method.State equation of SINS was adopted as the system state equation. Simulation experiments were conducted and results showed that the proposed method can estimate the flexural lever arm more accurately, the precision of transfer alignment was improved and alignment time was shortened accordingly.
Keywords: Transfer alignment, SINS, lever-arm effect, AUV, observability analysis.
Transfer alignment is the process of initializing the position, velocity and attitude of a slave INS, using the data supplied by another INS known as the master inertial navigation system [John and Leondes (1972)]. As the initial attitude errors cause the navigation errors to increase much rapidly than the initial velocity and position errors [Cheng, Wang and Liu (2014)], thus the relative attitude error of the slave INS with respect to the master INS is a major error source to result in the position error growth after launching the inertial guided weapons [Zhu and Cheng (2013)]. A digital filter method was applied in compensating the lever arm effect [Xu and Wan (1994)]. Special attention has been paid to the estimation of the slave INS attitude estimation. In the case of shipborne slave INS,both the time taken to achieve an alignment and the accuracy which can be achieved are highly dependent on the angular motion of the ship, and there are still some difficulties in aligning slave INS with the mater reference [Kassas and Humphreys (2014)]. For example,the application of velocity matching [Sun, Xu and Liu (2016)] alone is greatly limited in the shipborne alignment because ships are unable to make such maneuver as an aircraft do. In addition, the dynamic state of different points [Zhu, Cheng and Wang (2016)]. The location of the slave INS is not usually equipped at the swaying center, resulting in a lever effect on the system, leading to an impact on the initial alignment. The optimal estimation theory [Li,Wang and Liu (2014)] was introduced in the master-slave inertial navigation system supposing the master inertial navigation system was already calibrated and initialized accurately. The velocity error between the master inertial navigation system and the slave inertial navigation system can be regarded as the measurements in the system. An angle rate matching method [Wang, Gu, Ma et al. (2017); Zhou, Wang, Wu et al. (2017)] was applied in estimating the angle error. The velocity vector projecting on the inertial frame differences can be used in the angle error between the two frames requiring the angle movement Butterworth filter was designed according to the real acceleration and the leverarm acceleration distributed in the different frequency zones restraining the bad acceleration and compensating the lever arm effect [Chen, Chen and Wu (2017)]. A SINS error model under swaying base was constructed [Xiong, Peng, Wang et al. (2015)],analyzing the ship-board weapon self-initialization procedure [Wang, Shen and Ma (2008)]and proposed a new lever-arm compensation method under the base of traditional lever-arm compensation method. The novel method [Chen, Liu, Cheng et al. (2017); Fu, Huang, Ren et al. (2017)] was proved more accurate and with high precision with simulation experiments. An improved fast transfer alignment algorithm was proposed [Xia, Zhang and Lei (2017); Liu, Xu, Zhao et al. (2014)] and derived the construction method of velocity and z axis misalignment measurement [Fu, Wu, Wang et al. (2017); Sun and Gu (2017)].The reduced-order transfer alignment filtering model [Guo, Xu and Li (2017); Hong, Lee,Chun et al.(2005)] was presented and the performance of the reduced order filtering model is validated by simulation.
The remainder of this paper is organized as follows: the description of novel initial alignment error models is presented in Section II. The Observability analysis of the transfer alignment is given in Section III. Simulation experiments are conducted and the simulation results are given in Section IV. The conclusion can be found in Section V.
Firstly, define relevant basic frames:
i frame-earth centered inertial frame;
e frame-earth centered earth fixed frame;
g frame-east-north-up geographic frame;
m frame-body frame of master INS;
s frame-body frame of slave INS;
Assuming that L represents the latitude,λrepresents longitude,ωieis the earth rotation angular velocity, the local acceleration due to the actual gravity is denoted by g,φis defined as the misalignment angle.(ObXbYbZb)stands for the body frame,(OeXeYeZe)stands for the earth frame,(OnXnYnZn)stands for the navigation frame,(OiXiYiZi)represents the inertial frame. The subscripts b, e, i, n denote the corresponding components in each frame.
Suppose (OiXiYiZi)represents the inertial frame,(ObXbYbZb)stands for the body frame.Obis the swaying center of the ship, also is known as the gravitational center of the body.
Figure 1: Lever arm effect sketch
The position of the center of gravity is calculated according to the general design of the load distribution, assuming that the gravity center is fixed and the master INS installation position is in coincidence with Ob, the slave INS accelerometer is fixed point P mounted on the carrier coordinates, position vector coordinate origin, P point position relative to the inertial coordinate vector at P, relative to the position vector vector coordinate origin.Obviously, they have the following relationship:
The position of the gravity center is usually calculated according to the load distribution of design, which is always fixed and coincides with the master inertial navigation system equipping position. The accelerometers of slave inertial navigation system are equipped at the fixed point p in the body frame. The position vector of the origin of the carrier coordinate system is, the position vector of p to inertial frame origin is, the position vector of p to the carrier frame origin is, abiding the following rules:
The differential of the upper formulae to time can be obtained:
The differential of the upper formulae to time can be obtained:
According to the relative micro quotient principle of vector differential, it can be obtained:
The same reason can be obtained:
The linear acceleration of p to inertial frame can be expressed as:
The carrier structure is rigid in the study of lever arm effect. Other methods are needed to compensate the error caused by the flexible deformation. P is fixed to the carrier frame here.
The installation point should be in the swing center of the carrier in the ideal situation,=0, the lever arm effect did not exist in the situation. The master inertial navigation is equipped at the carrier swing center, while the slave inertial navigation system cannot satisfy the needs in transfer alignment. Lever arm effect cannot be neglected in real application. The latter two components are caused by the lever arm effect, sensed by the slave inertial navigation system, not sensed by the master inertial navigation system.stands for the lever arm acceleration, the standard equation of the lever arm effect error is,
For
The swing angle rate and angle rate differential in the carrier frame to the navigation frame are as follows,
For
Suppose the vehicle did not move linearly, that is, the above formulae can be simplified as
The lever arm effect acceleration can be expressed in the navigation frame as follows
Suppose
Suppose
And for
Can be derived that
Before AUV entering the water, transfer alignment is conducted before integrated navigation period. Transfer alignment uses the navigation information of position,velocity and attitude from master inertial navigation system equipped on the big ship.Observability analysis was conducted on the SINS state vector based on SVD method.State equation of SINS was adopted as the system state equation. Observability analysis was conducted under different vehicle motion with SINS/DVL integration and different matching methods under three-axis swaying mode.
Suppose the vehicle was navigating at a certain depth, the navigation velocity was 10 m/s,navigation depth was 20 meters. The initial latitude L=32othe initial longitude λ=118o,the gyro constant drift was1o/h ,the gyro random drift was 1o/h, the accelerometer constant drift was50 μg ,the accelerometer random drift was 50 μg, sampling time was 1 s.Observability analyses of velocity matching were conducted under different vehicle motions. Simulation experiments were conducted under stationary base, dual-axis swaying, three-axis swaying, acceleration and turning motion with velocity matching.The measurement information is the vehicle velocity (north velocity and east velocity)and velocity from the master inertial navigation system.
The motion process of initial alignment on moving base can be designed as follows, the gyro drift is estimated under three-axis swaying and linear acceleration. The accelerometer bias is estimated under linear acceleration and turning, the angle error can be estimated under linear constant navigation.
The observability degree of the state variables under three-axis swaying with velocity,heading angle error and position matching is presented in Tab. 1.
Table 1: Observability degree under three-axis swaying with velocity, heading angle error and position matching
From Tab. 1, eastern velocity error, northern velocity error, up heading angle error,longitude error, latitude error, eastern gyro drift, northern gyro drift, up gyro drift are high in observability degree and the filter effects are good. Eastern accelerometer bias and northern accelerometer bias are low in observability degree and the filter effects are not good.
From the analysis above, the velocity error can be well estimated with velocity matching while heading angle and position error cannot be well estimated. The heading angle error can be better estimated with velocity and heading matching, while the position error cannot be well estimated. The heading angle error, velocity error, position error can be well estimated with velocity, heading and position matching.
When the ship is under mooring condition, the linear acceleration and thelinear velocity is zero, the error equation of the inertial system can be reduced as follows,
n'is the navigation frame, b is the body frame,δVis the body velocity error,φis the vehicle attitude error,is the earth rotation rate on the navigation frame,is the computational error,gn'is the gravitational acceleration on the navigation frame.is the lever arm acceleration,?bis the acceleration constant drift,εbis the gyro constant drift,rbis the lever arm length on the body frame,is the strapdown matrix.
In the strapdown inertial navigation system, Kalman filter is applied in the estimation of lever arm length. The accelerometer error and the gyro drift gyro are expanded as the states in the Kalman filter. The system state equation is as follows,
State variables and the noise vector are as follows,
Coefficient matrix
φis the latitude,Ris the earth radius,is the gyro angle e velocity,is the angle acceleration.
The measurement variable isand, the measurement equation is
In the simulation experiments, the performance of each transfer alignment algorithm is evaluated. The ship speed is about 30knots. The angular motions of the ship are generated as follows
Suppose the initial longitude is 118°, the initial latitude is 32°, the initial height is 0. The initial longitude error is 5’’, the initial latitude error is 5’’, and the initial height error is 5’’. The initial eastern velocity is 10m/s, the initial northern velocity is 10m/s, the initial eastern velocity error is 0.1m/s, the initial northern velocity is 0.1m/s, the initial heading H0is 45°, the initial pitching angle is 0°, the initial rolling angle is 0°. Initial heading,pitching and rolling angle error is 10′′,10′′,30′′respectively. Accelerometer constant drift and random drift is 100ug. The gyro constant drift and random drift is 0.1°/h. The swing amplitude of heading, pitching and rolling is 14°, 9°and 12° respectively. The swing period of heading, pitching and rolling is 6 s, 8 s, and 10 s respectively. The ship swing model which is suffering the wind and tide under the mooring condition is as follows,
The lever arm length is [2m2m1m]T, DVL velocity error is 0.1 m/s, the DVL error coefficient is 0.0001.
This paper presents an efficient transfer alignment approach for swaying base big ship navigation system. A novel algorithm associated with the transfer alignment is employed to obtain the mathematical platform for the navigation system. Considering the environmental disturbances and the sensor drift as the main error source, a nonlinear filter approach is applied to the system to reduce the lever arm effect on the acceleration measurement. Observability analysis is conducted to different vehicle motion. Simulation experiments were conducted and the results showed that the novel method is able to improve the rapidity and precision of transfer alignment, overcoming the lever arm effect and disturbances existing in the host inertial navigation system and the slave inertial navigation system in the application of big ship navigation initial alignment.
Acknowledgement:This work is funded by Natural Science Foundation of Jiangsu Province under Grant BK20160955, a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions and Science Research Foundation of Nanjing University of Information Science and Technology under Grant 20110430.Open Foundation of Jiangsu Key Laboratory of Meteorological Observation and Information Processing (KDXS1304), Open Foundation of Jiangsu Key Laboratory of Ocean Dynamic Remote Sensing and Acoustics (KHYS1405).
Computers Materials&Continua2019年4期