Keyhole Performans
Keyhole Performans
ABSTRACT
The common two-axis azimuth-elevation gimbaled pedestal has a full-hemispheric, horizon-to-zenith field of regard.
This pedestal has no kinematic difficulties at low elevation angles. In this position, the line-of-sight of the mounted
sensor is perpendicular to both the azimuth and elevation gimbal axes, which thus provide two orthogonal degrees of
freedom. However, as the line-of-sight approaches zenith, the sensor axis nears alignment with the azimuth axis. The
azimuth axis thus loses its ability to move the line-of-sight orthogonally to the sweep of the elevation axis. This
condition is known as gimbal lock and the position range in which dynamic difficulties occur is called keyhole. The
keyhole region is a solid cone centered around the zenith axis. The onset of dynamics difficulties is a continuum from
horizon to zenith, and as such defining the keyhole region is arbitrary. However, dynamic difficulties become rapidly
pronounced somewhere between 70 and 80 degrees, so it is generally agreed that the keyhole region starts in this range.
This paper provides a comprehensive analysis of the keyhole region. While performance problems at keyhole are well
known (high torque, acceleration, and speed requirements), certain dynamic effects actually reduce in keyhole, such that
for some systems the range of worst-case performance is actually outside the keyhole region. Gimbal geometry is
introduced and pointing equations derived using vector methodology. Kinematic equations are then developed, with a
focus on the requirements of line-of-sight stabilization for vehicle-mounted systems. Examples and simulation results are
provided to illuminate the issue.
Keywords: gimbals, pedestals, keyhole, gimbal lock, gyroscope, pointing, stabilization, tracking
1. INTRODUCTION
The two-axis, elevation-on-azimuth, “az-el” gimbaled pedestal is widely used to pointing and stabilize a wide variety of
optical and RF systems, including telescopes, cameras, night-vision systems, lasers, and communication and radar
antennas. Az-el pedestals are found in all environments, including ground-fixed, ground-mobile, shipboard, airborne,
and spacecraft applications. The kinematic and dynamic difficulties of these pedestals at high-elevation angles, known as
the keyhole or gimbal lock is well known. While az-el pedestals are feasible in some shipboard applications, such as
littoral water craft or low-performance consumer products, the wide-use of three-axis pedestals in shipboard applications
is in fact a direct consequence of the keyhole limitations of the two-axis az-el gimbal. The keyhole problem is so well
known that it is tempting to assume that the worst-case performance of a two-axis gimbal will always be in the keyhole
region. Surprisingly, there are circumstances in which the worst-case performance of a system using an az-el gimbaled
pedestal will occur away from the keyhole region.
The kinematic, dynamic, and control equations for a two-axis gimbal are presented in this paper. Performance in
stabilized mobile applications is investigated, as it is the stabilization application that provides the most variety of
performance scenarios. The effects of pointing geometry, vehicle motion, friction, motor current and voltage limits,
direct and indirect drive configurations, and noise is investigated and discussed. Results are shown analytically or with
simulation results as best illuminates the issue.
A
d& 1 = ω × d 1
A D
(1)
The first term in the vector cross product is the angular velocity of elevation gimbal D within local-level reference frame
A. This quantity can be generated using the addition theorem of angular velocity:
A
ω D = Aω B + B ω C + C ω D (2)
The first term on the right is the angular velocity of the vehicle B in the local-level reference frame A, which can be
expressed as follows:
A
ω B = ω 1 b1 + ω 2 b 2 + ω 3 b 3 (3)
The scalar components ω1, ω2, and ω3 are often referred to as the body rates of the vehicle. These angular rates,
sometimes designated in the literature as P, Q, and R, are the quantities measured by a strap-down triad of gyroscopes as
part of an attitude and heading reference system (AHRS), inertial reference unit (IRU), inertial navigation system (INS),
or other inertial measurement unit (IMU).
The last two terms on the right side of equation 2 are the angular velocities of the azimuth and elevation gimbals C and
D in the base and azimuth reference frames B and C respectively. As both gimbals move about axes fixed in the
referenced reference frame, these quantities are simple angular velocities, readily expressed as:
B
ω C = θ&c 3 (4)
and:
C
ω D = ϕ& d 2 (5)
The cross-product of equation 1 is readily handled by expressing equation 2 in the vectors of the line-of-sight reference
frame D. The substitutions are found by inspection of the vector triad geometry:
Stabilization is achieved when the time-rate-of-change of the line-of-sight vector, given by equation 11, is equal to zero.
Equation 11 is uniformly zero only when the d2 and d3 scalar components are zero.
3. STABILIZATION CONTROL
Comparison of equations 10 and 11 yields an alternate definition of stabilization: The line-of-sight d1 is stabilized when
the d2 and d3 scalar components of the angular velocity of line-of-sight reference frame D (which here is the elevation
gimbal) are zero. This definition is suggestive of the most common servo configuration for stabilization. Specifically,
gyroscopes are mounted to the elevation gimbal so that they are aligned in the d2 and d3 directions. Control loops are
then implemented using the gyroscopes as feedback, with motors on the elevation and azimuth axis. The control loops
act to hold the gyroscope outputs to zero, thus satisfying the kinematic conditions for stabilization.
Vector d2 is aligned with the elevation axis, and the gyroscope aligned in this direction is called the elevation gyroscope
with no ambiguity. The vector d3, which is perpendicular to both d1 and d2, is aligned with the azimuth gimbal axis only
when the elevation angle is zero. The gyroscope aligned in the d3 direction is commonly called the azimuth gyroscope. A
more technically accurate but less commonly used name is cross-elevation gyroscope. While gyroscope names are not
critical, it is important however to distinguish that angular rate sensed by this gyroscope (specifically, cross-elevation
rate, which is the d3 component of equation 10) is not the same as the inertial rate of the azimuth gimbal.
3.1 Elevation
The block diagram for the kinematics of elevation stabilization control is shown in Figure 1. The primary, “T=Jα”
dynamics of the elevation gimbal are included. Secondary dynamics, such as imbalance, centripetal or product-of-inertia
terms, are not shown on the diagram, but would be included if significant in the torque term. The dynamic equation is
formulated here in inertial terms, such that the output of the acceleration-to-rate integrator is inertial rate. For the
elevation axis, this term, designated in the figure as “Elevation LOS Inertial Rate”, is sensed directly by the elevation
gyroscope. The other terms in the d2 component of equation 10 are shown at the bottom of the diagram and are used
within the model to generate the elevation angle as shown.
Figure 2. Block diagram model of azimuth and cross-elevation kinematics and dynamics
Figure 3. Simplified model of cross-elevation (azimuth) rate loop with roll-coupling input ω1.
G 2 (ω ) = sin(ϕ ) G RC (ω ) (14)
Equation 13 indicates that, for any given frequency, the roll coupling gain is zero at horizon (φ = 0) and maximum in the
keyhole region (φ approaching 90°). However, note that sine increases quickly from zero as elevation angle increases
from horizon. The rise is rapid enough that roll coupling becomes a significant error source at angles well below the
keyhole region.
The following example is provided for clarification of this effect. A rate loop with PI compensation and a 15 Hz open-
loop crossover has an loop gain of approximately 60 at 1 Hz, which corresponds to a gain GRC(ω) of -36 dB. The roll
coupling gain as a function of elevation angle for this system is calculated using Eq. 14 and shown in Figure 4. The
curve rises quickly, with roll coupling gain within 3 dB of its maximum value at 45 degrees and within 1 dB at 60
degrees.
Figure 4. Roll coupling gain as a function of elevation angle for example system.
Figure 5. Model of cross-elevation (azimuth) rate loop expanded to show secant correction
term, normalized azimuth gimbal inertia, and azimuth gimbal acceleration αAZ.
The increase in acceleration requirements at keyhole in roll coupling conditions can be shown by generating the transfer
function from ω1 to azimuth acceleration αAZ. Observation of the loop in Figure 5 shows this transfer function,
designated here as G4(s), to be:
α AZ ⎛ s G 3 ( s) ⎞
G 4 (s) = = tan(ϕ )⎜ ⎟ (15)
ω1 ⎜ s + G 3 ( s) ⎟
⎝ ⎠
At any given frequency the second term is constant, indicating that torque requirements to stabilize in roll coupling
conditions increase with the tangent of the elevation angle. Similar analysis shows that the torque requirements in
response to a commanded rate ωCMD increase with the secant of the elevation angle. These two curves are shown in
Figure 6.
Figures 7a and 7b. X-Y error plots for example system in ship’s motion with no current limits. Figure 7a on the left is at
low elevation angles (30 degrees) and the Figure 7b on the right is for elevation angles near keyhole (80 degrees).
The sign function makes the response of this system definitely non-linear, but under any conditions friction is an
unwelcome torque disturbance that degrades the level of stabilization performance. Though the system itself is non-
linear, the degree to which friction disturbs the cross-elevation rate can be investigated by ignoring (for now) the non-
linear feedback path from ωAZ to friction and using the same linear analysis technique as before. In this way, the transfer
function G5(s) from the friction input TF to the cross-elevation rate ωXEL can be derived from observation of the block
diagram in Figure 9:
ω XEL cos(ϕ )
G 5 (s) = = (16)
TF s + G 3 (s)
Transfer function G5(s) characterizes the torque disturbance rejection capabilities of the system. Like roll coupling,
torque disturbances produce an undesirable effect, so we would like the gain of the torque disturbance transfer function
G5(s) to be as low as possible. Disregarding the cosine term, the gain of G5(s) is made small by making the gain of G3(s)
large. That is, increasing the control gain of the stabilization loop improves torque disturbance rejection by driving the
gain of the torque disturbance rejection curve down.
The effect of the cosine term in Eq. 16 is similar but opposite to the effect of the sine term in Eq. 15. Using the same
formulations as Eq. 13 and 14, and the baseline torque disturbance gain GTD as the gain value of 1/(1+ G3(s)) at a given
frequency ω yields:
1
GTD ≡ (17)
1 + G3 (ω )
The gain of the torque disturbance rejection curve is expressed as:
Equation 18 indicates that the magnitude of the torque disturbance rejection curve (for any given frequency) is maximum
at horizon (φ = 0) and diminishes rapidly in the keyhole region (φ approaching 90°). This means that the sensitivity of
the system to torque disturbances decreases with increasing elevation angle. Note however that cosine does not decrease
rapidly from unity as elevation angle increases from horizon. As such, the sensitivity to torque disturbances is relatively
flat until the keyhole region is approached.
Figure 10. Torque disturbance gain as a function of elevation angle for example system.
Note that, while Eq. (16) is formulated as the transfer function from the friction torque TF to the cross-elevation rate
ωXEL, it is in fact independent of the value of TF. As it is, the torque disturbance transfer function is valuable tool to
investigate the effects of any torque disturbance TD that acts on the azimuth gimbal. This fact is used in the next section
to investigate the effect of other torque disturbances.
Note also that, though the torque disturbance transfer function is independent of fiction level, it is in fact the combined
effects of torque disturbance gain and friction level and frequency content that determine actual stabilization
performance. A common measure of the friction level in a system is the friction-to-inertia ratio. The friction-to-inertia
ratio has units of angular acceleration, usually radians-per-second-squared though degrees-per-second-square are
sometimes used. The friction-to-inertia ratio is a measure of the amount of acceleration the friction level can produce,
and in that regard friction-to-inertia ratio is a method of normalized the friction level to allow direct comparisons of the
effect of friction in systems of different sizes and different inertias. Note that in the normalized inertia formulation used
here, the friction level KF is the friction-to-inertia ratio. Large gimbals, even with moment-bearing, preloaded duplex
bearings, typically have friction-to-inertia levels of under 1.0 rad/s2. Small gimbals, such as found in miniature missile
seekers or the polarization mechanisms of stabilized antennas, can have friction-to-inertia ratios in the 20 to 40 range or
higher. This wide range plays into later discussions of total system stabilization performance.
Figure 11. Indirect-drive torque disturbance model for a belt-driven, co-rotating, motor-on-base azimuth drive
4.5.2 Noise
By far the largest source of noise in stabilization loops is sourced within the feedback gyroscope. To a great extent, the
cost of gyroscopes is in inverse to the noise level. Generally speaking, any gyroscope used for gimbal stabilization will
have a noise floor far in excess of the noise floors of any other electrical component or process in the loop. A model to
investigate the effects of noise as a function of elevation angle in azimuth loops is shown in Figure 12.
The torque disturbance TD from the gyroscope noise εXEL is embedded within the motor torque TM and is calculated from
the figure to be:
ω XEL G 3 (s)
= (20)
ε XEL s + G 3 ( s)
Equation 20 indicates that the disturbance to the line-of-sight from gyroscope noise is independent of the elevation
angle. This is a somewhat non-intuitive, especially if you have witnessed the audible growling and clearly visible
shaking of the azimuth gimbal of a high-gain system as the elevation angle increases into the keyhole region. However,
the increase in azimuth motion from noise amplification of the secant correction term is indeed offset exactly by the
reduction in mechanical gain of the azimuth gimbal to cross-elevation rate.
Secant amplification of gyroscope noise is, however, anything but benign. The noise amplification increases motor
current requirements, heats the motor, and reduces the voltage and current dynamic range headroom. It can indeed lead
to a reduction in stabilization performance, either through gradual means, such as early saturation of the azimuth motor,
or through hard-limiting such as servo-amplifier foldback.
Take, for example, a small satellite communications antenna mounted to an azimuth-elevation pedestal designed to
operate on on-road and off-road ground vehicles within the continental United States. Maximum satellite elevations in
the United States are in the approximate range of 35 to 55 degrees elevation. The roll coupling gains are not excessive
for these angles, so it is quite feasible to provide an azimuth drive with the overhead needed to avoid current limiting and
motor saturation. If this system has a high friction-to-inertia ratio and/or a high-ratio indirect drive, it is very likely that
maximum stabilization errors will occur at the lower elevation angles. This is especially true given that ground vehicles,
especially in off-road conditions, tend to undergo frequent yaw reversals, thus adding high-frequency content to the
azimuth friction and further increasing the effects of friction and indirect-drive coupling on line-of-sight stabilization.
ACKNOWLEDGEMENTS
The author wishes to thank Rohit Murthy of General Dynamics Satcom Technologies, Richardson, Texas, for his
encouragement and contributions during the preparation and review of this manuscript.
[1] Rue, A., “Stabilization of precision electro-optical pointing and tracking systems,” IEEE Trans. Aerosp. Electron.
Syst., AES-5(5), 805–819 (1969). (Errata: Nov., 1970, pp. 855–857).
[2] Hilkert, J. M., “Inertially stabilized platform technology,” IEEE Contr. Syst. Mag., 28(1), 26–46 (2008).
[3] Oslon, A.R., “Shipboard Stabilization of Large Antennas”, ASME Paper 65-WA/AUT-8 (1965).
[4] Kane, T.R. and Levinson, D.A, [Dynamics, Theory and Applications], McGraw-Hill, New York (1985).
[5] DeBruin, J.C. and Johnson, D., “Line-of-sight reference frames: A unified approach to plane-mirror kinematics,”
Proc. SPIE 1697, 111–129 (1992).
[6] Ellison, B. and Richi, J., “Inertial stabilization of periscopic sights band driven three axle gimbal”, Proc. SPIE 389,
107-119 (1983).
[7] Masten, M., “Inertially stabilized platforms for imaging optical systems,” IEEE Contr. Syst. Mag., 28(1), 47–64
(2008).