A Transfer Alignment Algorithm Study Based On Actual Flight Test Data From A Tactical Air-To-Ground Weapon Launch
A Transfer Alignment Algorithm Study Based On Actual Flight Test Data From A Tactical Air-To-Ground Weapon Launch
--
Abstract This paper provides results from a joint study attitude navigation states. The following is a discussion
of the T/A process and the hardware requiredto perform
of various Transfer/Alignment (T/A) techniques
including velocity matching, integrated velocity a transfer/alignment.
matching, and doubly integrated Velocity matching
approaches. The goal is to determine how the various Transfer Alignment
algorithms perform when driven with actual flight test
data. The study addresses key T/A Kalman filter issues A T/A method is the process of transferring the high
such as the effect of wing flexure vibration on accuracy navigation information from an aircraft INS to
instrument error estimation and rapid alignment a weapon INS which is typically located on a pylon
concepts. A description of how a T/A Kalman filter is station of the aircraft‘swing or on a launcher located in
implemented and used in inertially guided weapons is a remote weapons bay. A direct transfer of aircraft INS
included. The data used in this study was acquired from navigation states (position, velocity, and attitude) would
the Air Force Operational Concept Demonstration not account for the fact that the weapon is not pointed in
(OCD) Program. The OCD effort was performed by the the exact same direction as the aircraft. In other words
Air-To-Surface Weapons System Program Mice (SPO) the aircmft sensors and the weapon sensors are not at the
at Eglin AFB in late 1992 and early 1993. The Eglin same tilts (pitch, roll, and heading) with respect to the
SPO initiated the flight test program to demonstrate the local vertical coordinate frame. Not compensating for
effectiveness of GPS/INS integration in a tactical these tilt differences would lead to erroneous velocity
weapon. The test was conducted using a Honeywell and position estimates, resulting in undesired weapon
Integrated Flight Management Unit (IFMU) and an navigation performance. This limitation in weapon
Interstate Electronics 5-channel P-code GPS receiver in accuracy can be greatly reduced if the weapon has its
a GBU-15 airframe on board a block 40 F-16 aircraft. own INS and access to aircraft INS data. A T/A Kalman
Telemetry data from this OCD weapon was transmitted filter is typically nm on the weapon CPU and receives
continuously from the beginning of weapon alignment aircraft navigation information via a 1553 interface bus.
to final impact on the ground. The OCD telemetry The T/A Kalman filter utilizes the measurements from
stream contained INS and GPS data from both the the aircraft INS and a weapon INS to estimate the
aircraft and weapon during captive carriage T/A misalignment between the two systems (both with
maneuvers. After launch, the weapon data alone was respect to the local vertical), and, as a residual effect, to
transmitted during the free flight to the target. These calibrate the weapon INS sensors.
messages provide the data set used in the T/A study.
Objectives
INTRODUCI~ON
There are many issues concerning the T/A process.
Today’s users of aircraft weapons are requiring that After a minimal discussion of the Kalman filter as it
munitions be “smart.“ One genre of smart weapons applies to the alignment process, this paper will focus on
includes those capable of navigating to a predetermined the more signifkant issues concerning T/A algorithm
target coordinate autonomously. The navigation design including how long a weapon should align, the
performance of such a weapon is dependent on the effect of vibration on sensor calibration, the effect of
proper initialization of the vehicle navigation states aircraft to weapon lever arm miscalculations, and the
which include position, velocity, and attitude. The T/A choice of T/A measmment methods.
approach to initializing a weapon is the most accurate
technique presently used to initialize the weapon
431
0-7803-1435-2/94 $3.00 0 1994 IEEE
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.
KALMANF ~ T ECONCEP~S
R APPLIED estimate the majority of the alignment error from the
TO TRANSFERALIGNMENT input velocity differences. To illustrate this point Figure
1 shows a filter's ability to estimate the misalignment
For the purposes of this study, it is assumed that the angle between the aircraft INS and weapon INS using
reader is familiar with the standard Kalman filter only a 6-state velocity matching filter. These estimates
equations [l]. For stability reasons, a Joseph's form are compared to the OCD 15-state integrated velocity
Kalman filter is implemented and a state correction matching filter that was initialized at t=O seconds and
vector is applied to close the loop [2]. The challenge is has stabilized. It is worth noting that the pitch and roll
to apply the Joseph's form Kalman filter equations to the Euler angle misalignments are observable within 5-10
T/A problem. When a Kalman filter is applied to the seconds due to the influence of the gravity vector. The
T/A process, there are two major benefits which result; azimuth Euler angle misalignment became observable at
alignment error estimation and sensor calibration. t=440 seconds when the aircraft executed a 60 degree
heading change. All subsequent misalignment angle
Alignment Error Estimation plots will compare the T/A study results with the OCD
stabilized estimates.
The most important elements of the T/A Kalman filter
are found in the 6 by 6 element sub-matrix of the system
dynamics matrix composed of three tilt <Yx,Y, ,Yz ) O M
Euler Anpk Uhallpnmnnl Estlmatn : OCO n O-!jtolr VU flw
error equations and three velocity (6V,, 6V, ,SV,) error Note 0'3 flrta Itorled at 1-0 and ha8 rmchmd rlmdy rtaIn.
0.03.
equations. The sub matrix of the system dynamics modinp
matrix is given below.
0 -T3 T2 I 0 0 0.01
[T3 0 -T,i 0 0 0
-------------,--------------
-T2T,
A ( r ) X ( r ) = 0 -a,-5 Z2
O i O
ay I 0
0
a,
-ag
o
a,
-a,
o
I Z,
1 -z, Z,
O -"I
o
-0ML
420 440 460 480 Hx) 520 540
I
5M)
where: nc
Figure I.
A (t ) = KF State Dynamics Matrix
X ( t ) = KF Error State Vector Sensor Calibration
Y = Misalignment Error Estimation Vector
6 V = Veiociv Error Estimation Vector Sensor calibration can be achieved by adding sensor
a , , a y, a , = Local Level Acceleration Vector error states to the Kalman filter. It would be ideal to
model all of the known sensor errors in the filter;
T,=-(p*-R,) however, such an endeavor reduces computational
T2 = P y + Q y capability and requires more memory and throughput in
T3 = - a , the weapon CPU. The alternative is to model only the
z, = -p, - 2 Q x most significant sensor error states and compensate for
2, = p y +2Ry the unmodeled states by increasing the appropriate
parameters in the Kalman filter process noise matrix.
z, = 2R, The sensor error states modeled in the Kalman filter
p, ,p, = x and y Local Level Craft Rate in radfsec used for this study include accelerometer bias and gyro
R, , a yR,
, = Local Level Earth Rates in radfsec bias.
432
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.
velocity estimates of two navigation systems. The filters make up the plant. The model presented below
implemented for this study incorporate the difference represents the doubly integrated velocity matching plant
between aimaft INS velocity measurements and with the associated m r state vector. This plant was
weapon INS velocity measurements. modified to support analysis of the velocity and
m et
m velocity matching techniques. Detailed
There are three methods [l] of processing the velocity infonnation on these techniques are presented in
measurements recently evaluated through simulation for reference [l].
the Naval Weapon Center (NWC). They are the
velocity matching 0, integrated velocity matching
0, and doubly integrated velocity matching
(DIVM) techniques.
TESTINOTRANSFTWALIGNMENT
433
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.
Modem researchers have become interested in alternate good thing, then to integrate twice should be even
measurements such as angle rate matching, integrated better. The goal is to further reduce the random noise
velocity matching, and doubly integrated velocity on the measurements, that should allow more accurate
matching. In this study, the angle rate matching estimation of tilt mrs and sensor biases as presented in
approach is not investigated due to the complexity of [l]. The doubly integrated velocity measurement is
derivation and implementation. Research is being given below.
performed under an Air Force contract with TASC (The
Analytical Sciences Corporation) to investigate the
effectiveness of angle rate matching [5]. dDj,=j
T
2-0
I
T
t=o
Vm*dt - IT IT
f=O z=o
Va*dt
Velocity matching is the standard approach of
differencing the velocity measurements of the aircraft mi,,= Doubly Integrated Velocity
and the weapon. The input to the Kalman filter is given
below. Matching Measurement
AV, = V,,,- V, Figures 2.3 and 4 illustrate the performance of the three
V,,,= Missile Velocity alignment methods (VM. IVM, DIVM) employing a 1
second filter cycle time implemented for this study.
V, = Aircraji Velocity Hgure 2 compares the pitch error estimates for the three
AV,,, = Velocity Matching Measurement methods and includes the OCD real-time estimate as a
baseline. The three methods agree roughly within 1-2
The inputs are processed by the Kalman filter and milliradians. The velocity matching method shows an
estimates of the tilt errors and accelerometer and gyro initial overshoot which is quickly corrected. This
bias are estimated. The argument over which phenomena was not seen on most of the data runs and is
measurement type to use centers around the ability to attributed to initialization errors. A steady state
estimate accelerometer and gyro bias. estimate of the pitch misalignment for all three
approaches is achieved in approximately five seconds.
Integrated velocity matching is similar to velocity There is no discernible difference between the three
matching with the exception that the measurement is the methods that would suggest one is better than the others
difference between the integral of the aircraft and the for pitch estimation.
weapon velocity estimates. In effect, the measurement
becomes the difference of position changes estimated
over the measurement cycle time. The advantage of
differencing the integral of the velocity measurement is
that the random noise which corrupts the velocity
measurements will be attenuated due to the integration
process. A reduced measurement noise should improve
filter performance in the estimation of instnunent errors
such as accelerometer and gyro bias. The integrated
velocity measurement equation is given below.
434
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.
velocity matching converges approximately three to five Measurement Update Rates
seconds later. The three methods converge on the OCD
estimate of misalignment angle without any visible The choice of the measurement technique is influenced
differences. by the T/A filter update rate. The velocity malchmg
technique degrades at the longer filter cycle times (3-5
seconds) [l], therefore a common filter update rate of
one second was chosen to evaluate the three
measurement techniques. The OCD program which is
used as a baseline, implemented a five second update
rate with integrated velocity as the measurement.
k v e r Arm Corrections
VibrationEffects
435
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.
introduction of wing flexure vibration tends to impede
the transfer alignment process. The estimation of sensor
errors such as accelerometer and gyro bias is hindered
by modeled aircraft wing dynamics such as flexure
and vibration. The promised benefits of the integrated
measurement techniques are that the input to the T/A
KaLman filter has reduced noise properties. This
reduction in noise should result in more accurate
estimation of sensor biases. Figures 5 and 6 illustrate
the vibration experienced by a weapon system before
and after launch as measured by the weapon
accelerometers. Figure 7 shows the acceleration
characteristics through launch, including the peak -1 L A
588 588 570 572 574 576 578 580
acceleration experienced when the aircraft ejector
n
cartridges are fired.
Figure 7.
0.1 -
0.05 -
4.15t
i 4t I il -I
-02
-025 I
t I
592 594 596 598 600 602 604
I C
Figure 6.
0 5 1 0 1 5 2 0 2 5 3 0 3 5 4 0 z J 5 0
Hz
Figure 8.
436
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.
4-
3.5-
3-
2.5-
2-
1.5
0
1-
0
-
5 1 0 1 5 2
Figure 9.
5
Hx
437
_.. - _^.I. .
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.
the filters were all found to be sensitive to horizontal [6] Timothy F. Elbert, Paul F. Torrey, David M. Gaskill
accelerations at the initialization point. The integrated "INS/GPS OCD Final Report", Contract Number:
approaches exhibited greater sensitivity initial FO8635-91-C-OOO2, Sverdrup Technology, Inc., P.O.
conditions leading to larger overshoots of the steady Box 1935 Eglin AFB, FL 32542-5000, October 1993.
state values which led to a bias in the final steady state
estimate.
438
Authorized licensed use limited to: J.R.D. Tata Memorial Library Indian Institute of Science Bengaluru. Downloaded on February 20,2025 at 22:54:38 UTC from IEEE Xplore. Restrictions apply.