Optimal Control for Proximity Operations and Docking

International Journal of Aeronautical and Space Sciences.
2010.
Sep,
11(3):
206-220

- Published : September 15, 2010

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

This paper proposes optimal control techniques for determining translational and rotational maneuvers that facilitate proximity operations and docking. Two candidate controllers that provide translational motion are compared. A state-dependent Riccati equation controller is formulated from nonlinear relative motion dynamics, and a linear quadratic tracking controller is formulated from linearized relative motion. A linear quadratic Gaussian controller using star trackers to provide quaternion measurements is designed for precision attitude maneuvering. The attitude maneuvers are evaluated for different final axis alignment geometries that depend on the approach distance. A six degrees-of-freedom simulation demonstrates that the controllers successfully perform proximity operations that meet the conditions for docking.
Standard International Space Station (ISS) approach (adaptedfrom Walker et al. [2005]).
This trajectory satisfies the many constraints on visual and sensor visibility, plume impingement pressures and contamination, propellant consumption, and other factors(Cloutier, 1997). The reaction control system (RCS) of the Shuttle orbiter is used to provide thrusters for this approach trajectory. The 38 Primary Reaction Control System (PRCS)thrusters are arrayed around the Shuttle orbiter. The final orbit of the Shuttle orbiter’s rendezvous profile targets a point 183 m (600 ft) below the ISS along the R-bar direction. The Shuttle orbiter crew begins manual trajectory control at a distance of 610 m (2000 ft). The negative R-bar direction control is activated at 305 m (100 ft) to provide plume protection by inhibiting thrusters that fire toward the ISS. At the 183 m (600 ft) point, the Shuttle orbiter begins an 11.5 minute positive pitch automatic maneuver to the final ISS approach attitude.A simultaneous, manual V-bar translation of 0.3 m/s (1 ft/s)accomplishes a slow transition from the 183 m (600 ft) R-bar departure point to a final approach corridor along the ISS V-bar. In the LVLH frame, this transition appears as a gradual spiral from 183 m (600 ft) along the R-bar to approximately 350 ft (107 m) along the V-bar. From the arrival point along the V-bar, the Shuttle orbiter slowly approaches the docking port at a rate of 0.06 to 0.03 m/sec (0.2 to 0.1 ft/sec).
To effect translational control, the state-dependent Riccati equation (SDRE) (Cimen, 2008; Cloutier, 1997; Stansbery and Cloutier, 2000; Sznaier and Suarez, 2001) control and linear quadratic tracking (LQT) (Alba-Flores and Barbieri,2006; Budiyono and Wibowo, 2007; Lewis and Syrmos,1995; Naidu, 2003) control with the free-final state are both considered, and their control results are compared and discussed. To provide attitude control, a linear quadratic Gaussian (LQG) controller is used, which is a combination of a linear quadratic regulator (LQR) (Lewis and Syrmos, 1995;Naidu, 2003) and a Kalman filter. The LQR technique is a well-known and accepted theory upon which many modern controllers are based. However, most dynamic systems requiring control are nonlinear. The SDRE method entails factorization of the nonlinear dynamics into the state vector and the product of a matrix-valued function that depends on the state itself. With this factorization, the SDRE algorithm fully captures the nonlinearities of the system, expressing the nonlinear system as a non-unique linear structure known as a state-dependent coefficient (SDC) matrix that minimizes a nonlinear performance index having a quadratic-like structure. A state-dependent Riccati equation using SDC matrices is then solved in real time to yield the suboptimal control law. The SDRE control is based on the nonlinear system whereas LQT is based on the linear system. As an alternative approach, an LQT controller with the free-final state that maintains the output as closely as possible to the desired or reference output with minimum control energy is also considered.
CW and body-fixed reference frame.
In the above,
ax , ay ,
and
az
are unit vectors in the spacecraft local orbital frame, A, and the axes,
e
x
,
e
y
and
ez
, are unit vectors in the CW frame, E. The +
a
x
axis is referred to as the V-bar axis, the +ay axis as the H-bar axis, and the
+az
axis is termed the R-bar axis. The axes,
bx , by,
and
bz ,
are unit vectors in the spacecraft body-fixed frame used to describe the attitude with respect to the inertial frame.
where
R_{e}
is the radius of the Earth, i is the inclination, and θ is the argument of the latitude. The relative effect of the earth’s oblateness then becomes:
where the subscripts,
c
and
t
, denote the chaser and the target, respectively. The radius, inclination, and argument of the latitude of the target are assumed to be known from the on-board navigation system of the target vehicle. The radius,inclination, and argument of the latitude of the chaser are then used to relate the position and velocity vectors through Lee and Pernicka (2009)
where C is the 3-1-3 rotation sequence
C=C_{3}(θ_{t})C_{1}(i_{t})C_{3}(Ω_{t})
,
η
is the relative position, and
is the relative velocity. The ω×term is the cross product matrix. The perturbing acceleration in the CW frame due to aerodynamic drag is computed by expressing the acceleration in terms of the ECI frame(Madonna, 1997; Vallado and McClain, 2001) as
where ρ is the atmospheric density, which often is a difficult parameter to determine. The ballistic coefficient is c
_{D}
A/
m
,another measure of a spacecraft’s susceptibility to drag effects, and
V_{rel}
is the velocity vector in the ECI coordinates relative to the rotating atmosphere. The relative effect of atmospheric drag in the CW frame is then given by:
in which the drag is computed for both the chaser and the target orbit. Thus, the sum of the relative effects of Earth’s oblateness and drag becomes:
Consequently, the equations of motion for the SDRE control formulation become:
In the above,
x, y,
and
z
represent the relative position of the chaser spacecraft with respect to the target spacecraft, r
_{t}
refers to the scalar radius of the target from the center of Earth,r
_{c}
refers to the scalar radius of the chaser from the center of Earth, μ is the gravitational parameter, and
f
represents the true anomaly of the target. The terms, F
_{x},
F
_{y}
, and F
_{z,}
are the control forces, and
m(t)
is the time-varying mass of the chaser spacecraft due to the propellant mass consumption.The general linear relative equations of motion used for the LQT controller with the free-final state are derived from Eq.(10) by neglecting the perturbation effects and the higherorder terms (Schaub and Junkins, 2003) to yield
where
I_{c}
is the inertia matrix of the chaser. The orientation of the body-fixed frame, B, of the chaser with respect to the spacecraft’s local orbital frame,
A
, of the chaser is described by the direction cosine matrix,
C^{B/A}
, as indicated below.
The orientation of the local orbital frame,
A
, of the target spacecraft with respect to the CW frame,
C
, is described by the direction cosine matrix, C
A/C
, such that:
The direction cosine matrix, C
^{B/A}
, can be expressed using successive rotations with the inertial frame,
N
, through:
The angular velocity of the chaser,
ω=ω^{B/N}
, and
a
z
can be expressed in terms of the basis vector of the body-fixed frame, B, of the chaser as:
The gravity-gradient torque matrix becomes:
A full description of the rotational motion of a rigid spacecraft requires both kinematic and dynamic equations of motion.For most modern spacecraft applications, quaternion kinematics (Lefferts et al., 1993) are preferred. The quaternion kinematic equations are given by:
where
The adopted quaternion is defined by:
where is defined as [q
_{1}
q
_{2}
q
_{3}
]
^{T}
=
e
sin(ϑ/2), and q4=cos(ϑ/2),where
e
is the axis of rotation and ϑ is the angle of rotation.Euler's rotational equation of motion, including the gradient torque, is given by:
where
г_{c}
is the applied control torque and where
T_{d}
is the external disturbance torque, which is modeled by Gaussian-noise.
In the above,
x
_{r}(t
)∈R
^{6}
,
f
∈R
^{6}
, B∈R
^{6×3}
, R(x)∈R
^{3×3}
,
u
_{Fs}
∈R
^{3}
;Also,
x
_{r}^{(t)}
is the reference or desired state vector provided by the guidance scheme based on the CW state transition matrix (Sznaier and Suarez, 2001) and the straight-line V-bar approach (Fehse, 2003). Assume
f(x)
is continuously differentiable in
x
and
f(0)=0
; B is a constant matrix. The weighting matrices, Q(
x
)≥0 and R(
x
)>0, for all
x. u
_{Fs}
(x) is the control force exerted by the SDRE controller in the CW frame. The translational maneuvering of the spacecraft is accomplished through the use of the control force. The nonlinear equations of the spacecraft dynamics in Eq. (8) are written in the linear-like state-space form of Eq. (22) through the SDC parameterization. Then, the matrices are given by:
The SDRE method for obtaining a suboptimal locally asymptotically stabilizing solution of problem Eqs. (21, 22)is as follows.
Construct the following nonlinear feedback controller equation.
The resulting SDRE-controlled trajectory becomes the solution of the quasi-linear closed-loop dynamics.
Then, the state-feedback gain for minimizing Eq. (21) is:
The state weight matrix for the performance index in Eq. (21)is given by:
The control weight matrix is given by:
The weight matrices used from the initial time are readjusted at the steady-state conditions in order to reduce the steadystate tracking error The additional control forces are then generated. This tuning of the weight matrices is very important in order to generate suitable control forces with no thruster saturation. The readjusted weight matrix,
Q(x)
,is then 1,000 times the original
Q(x
(
t_{0}
)), and the readjusted weight matrix, R(
x
), is 1/100 times the original R(
x
(
t_{0}
)).
where
u
_{FL}
∈R
^{3}
is the control force in the CW frame from the LQT controller. The objective of the LQT controller is to control the system in Eq. (30) such that
x
(
t
) tracks the reference state,
x^{r}(t )
, as closely as possible during the interval [
t_{0}, t_{f}
] with minimum control effort with the quadratic performance index chosen as follows.
Further, the boundary conditions are defined as
x(t _{0})=x_{0}
and free
x
(
t_{f}
). The weight matrices, F(
t_{f}
) and
Q(t )
, are symmetric positive semi-definite matrices, and
R(t
) is a symmetric positive definite matrix. When
F(t^{f}
) becomes large, the free final state becomes a fixed final state. The optimal control law consists of the sum of two components given by:
where the symmetric, positive definite matrix, P(
t
)∈R
^{6×6}
,is the solution of the nonlinear, matrix differential Riccati equation (DRE), and the first term is a full-state feedback with Kalman gain. Thus,
with the final condition being
The vector,
g
(
t
)∈R
^{6}
, is the solution of the linear, nonhomogeneous vector differential equation,
with the final condition being
Whereas the SDRE is solved online, the DRE and the nonhomogeneous vector differential equations are solved offline before control is performed. The optimal state is the solution of the linear state equation,
When the linearized equations of the spacecraft dynamics in Eq. (9) are adopted, the time-varying system matrix, A(
t
),becomes:
The applied weight matrices, Q(
t
) and R(
t
), in Eq. (31) are equal to the matrices in Eqs. (28) and (29). The weight matrix for the final state,
F(t_{f})
, is set to a relatively larger value than the value of
Q(t)
to achieve the effect of the fixed-final state control; it is given by:
where
x
_{a}
includes the quaternion and angular velocity of the chaser. The general form is given by:
where
u_{г}= г_{c}R^{3}
is the control torque. The goal is to drive the state to zero with minimum control energy. The regulator problem is then formulated with the performance index,
where the subscript, “a,” is used here to differentiate from translational maneuvers. The reference or desired quaternion,
q
_{d}
, is defined that also obeys the kinematic equation,
where
ω
_{d}
is the desired angular velocity vector and
q
_{d}
is assumed to be provided by on-board navigation in the target spacecraft. The error quaternion is defined as:
where ⊗ denotes quaternion multiplication. Also, the quaternion inverse is defined by
q
_{d}
^{-1}
=[-ρ
_{d}
q
_{4}
]. This work adopts the convention of Lefferts et al. (1993), who multiply the quaternions in the same order as the attitude matrix multiplication. Then, δ
ρ
and δq
_{4}
can be shown to be given by the following.
As δ
ρ
approaches zero, the actual quaternion approaches the desired quaternion. Some properties of the quaternion error measurement are described below, making it possible to develop an exact linearization of the error dynamics for a formulation. Assume that the closed-loop dynamics have the linear form (Bach and Paielli, 1993; Paielli and Bach, 1993),
where L
_{1}
and L
_{2}
are 3×3 gain matrices. These matrices can be determined using an LQR approach starting with the following.
In the above,
L=[L_{1} L_{2}]
. The state-space formulation of Eq.(48) is given by:
where
The state weight matrix for the performance index in Eq. (42) is given by:
The control weight matrix is given by:
For the generation of effective control torque commands that avoid actuator saturation with good attitude tracking,the weight matrix,
Q_{a}(x
_{a}
(t
_{0}
)), is also readjusted at steady state as is the case with SDRE control. The ARE is solved to produce a constant gain matrix. The goal of this control is to determine a control torque,
г
_{c}
, in Eq. (19) that achieves the desired closed-loop dynamics given by Eq. (46). Toward this end, two time derivatives of Eq. (45a) are first taken and then substituted into Eq. (46) (Crassidis and Junkins, 2004)yielding:
Taking the time derivative of Eq. (16) leads to:
where the identity, Ω
^{2}
(
ω)=-(ω^{T} ω
)I
_{4×4}
, is used. An identical expression for the desired quaternion is also given as:
where
can be derived from Euler’s equations of motion using Eq. (19) (neglecting the unmodeled torque,
T
_{d}
).Substituting Eq. (16) into Eq. (52) gives:
Substituting Eq. (16) and Eq. (54) into Eq. (51) and solving for
г
_{c}
yields:
For precision attitude control, precise attitude sensor sensors are adopted and combined with LQR control. The goal of the EKF application is the estimation of the gyro biases along all three axes and the attitude quaternions of the chaser spacecraft. The sensors used here include three axis gyros and star trackers whose output is the attitude quaternion referenced to J2000 inertial coordinates. For this sensor, a widely used model is given by the following (Crassidis and Junkins, 2004).
In the above,
η
_{v}
and
η
_{u}
are independent zero-mean Gaussian white noise processes with covariances usually given by σ
_{v}
^{2}
I
_{3×3}
and σ
_{u}
^{2}
I
_{3×3}
, respectively, with σ
_{v}
= √10×10
^{-7}
rad/sec
^{1/2}
and σ
_{v}
= √10×10
^{-10}
rad/sec
^{3/2}
.
^{12}
β is a bias vector, and
is the measured observation. A combined quaternion from two star trackers is used as the measurement. To generate synthetic measurements, the following model is used:
where
q
_{m}
is the quaternion measurement,
q
is the truth,and
v
is the measurement noise, which is assumed to be a zero-mean Gaussian noise process with a covariance of 0.001I
_{3×3 }
deg
^{2}
. The measured quaternion is normalized to ensure a normalized measurement. An error quaternion between the measured quaternion and the estimated quaternion of the chaser,
is used for measurement in the filter. This is computed using the error quaternion,
For small angles, the vector portion of the quaternion is approximately equal to half angles so that δ
ρ
?
α
/2 and δq
_{4}
=1.
The estimated quaternions of the chaser are fed to the control torque,
г
_{c}
, in Eq. (54) as
q
. The implementation of the LQGtype control system for the chaser’s rotational maneuver is shown in
Fig. 3
.
Implementation of a linear quadratic gaussian-type control system
Geometry of axial alignment.
The quaternions of the target,
q
_{t}
, were found through the useof the quaternion kinematic equations. The target angularvelocity was assumed to be constant (consistent with anEarth-pointing attitude) during the scenario. The inertia matrix of the target was assumed as follows (Fehse, 2003).
The inertia matrix of the chaser was also assumed constantand taken as follows (Nagata et al., 2001).
The chaser was initially located at the relative position, [609.6213.36 0.1] m, with a relative velocity of [0.01 -0.43 0.01] m/swith respect to the target. The initial quaternion of the chaserwas [0.2473 0.4123 0.8651 0.1426], which corresponded to aninitial attitude of -30 degrees of pitch rotation with respectto the CW frame of the chaser. The initial biases along eachgyro axis were set as 1 deg/hr (Naidu, 2003). The dockingtarget location (Olszewski, 1990) was given as [27.30 12.71-2.74] m in the CW frame. The flight segments were dividedinto a total of six stages that comprised of sub-segments. Thefirst sub-segment started from the initial point, S1, wherethe final mid-course correction maneuver was executed andtargeted S2, a point located 183 m along R-Bar. The secondsub-segment started from S2 and continued to S3. At S2,the Shuttle started a positive pitch maneuver to the LOSof the target V-bar; the third was the first station-keepingto capture the LOS of the target; the fourth was the V-barapproach toward the target V-bar; the fifth was the secondstation-keeping location to decrease the approach velocityand ensure capture of the LOS; and the final sub-segmentwas the straight-line V-bar approach intended to accomplishsoft docking. The simulation lasted for 68 minutes fromthe holding point, S1, to the target docking port with a stepsize of 0.1 seconds. Whereas the translational control wasinitialized at the beginning of each sub-segment, the attitudecontrol was initialized twice to execute axial alignments. Thefirst attitude was initialized at S1, and the second attitudewas initialized at S2. Since the H-bar relative motion in thisscenario was much smaller than the R-bar and V-bar relativemotions, a planar motion perspective was used to betterillustrate the results. The results of translational maneuveringwere independently determined by both SDRE control andLQT control. The data were then plotted together from theholding point to the target docking port.
Proximity operation trajectories.
Out-of-plane relative motion.
Fly-around maneuver: S1 to S2.
After the second sub-segment was completed, a five minuteperiod of station-keeping was executed to ensure proper LOSalignment by nullifying the arrival velocity generated by theCW terminal guidance (Fehse, 2003; Prussing and Conway,1993), as shown in
Fig. 9
. After the fly-around maneuver fromS2 to S3 was complete, the chaser needed to align its attitude
Fly-around maneuver: S2 to S3.
First station-keeping subsegment.
with the target LOS to the docking port during the stationkeeping period. Five minutes was considered an adequate amount of time to align with the LOS during the stationkeeping phase. Once the LOS was acquired during the station-keeping phase, the chaser prepared for the straight V-bar final approach phase and additional control forces were generated so that the chaser would not drift. After the first station-keeping period was complete, a V-bar hopping approach along the LOS was executed for 25 minutes, as shown in
Fig. 10
. The nominal V-bar hopping trajectory and the transfer time were also determined by the CW terminal guidance scheme, which typically provides a fuel-efficient approach.
The chaser then executed a second period of stationkeeping prior to the straight LOS final approach, shown in
Fig. 11
, for five minutes. Unlike the first station-keeping period, the chaser’s arrival velocity was opposite to that of the arrival velocity from the second R-bar maneuver. During the second station-keeping period, the chaser slowed the arrival velocity, again aligned the LOS to the docking port,and prepared for the straight LOS final approach. After the second station-keeping phase was completed, the chaser spacecraft completed the straight-line V-bar approach to dock with the target, as shown in
Fig. 12
. A straight-line V-bar guidance approach with very slow constant velocity was
V-bar approach.
Second station-keeping sub-segment.
specified to ensure safety and avoid an unacceptably high physical impact on the target.
According to the docking constraints (Pearson, 1989) in use for the Shuttle and the ISS, the lateral docking tolerance is a maximum of 0.330 m (13 inches), the lateral velocity tolerance is 0.0457 m/s (0.15 ft/s), and the closing velocity tolerance is 0.0914 m/s (0.30 ft/s). It is especially important to note the position tracking error and approach range rate to monitor these conditions. In order to facilitate an effective assessment of whether or not these conditions are met, the final straight-line V-bar approach over the last 10 minutes was magnified. After control initialization in each sub-segment, the tracking error gradually decreased.This decrease in the tracking error is thus acceptable for achieving soft docking. As the translational maneuver was stabilized, the controlled state tracked the reference state more closely. The SDRE control achieved this accuracy as a result of adjusting the weight matrices in the final approach phase. However, LQT control with the free-final state could achieve very precise position tracking at the final stage by the use of a large weight-matrix, F(
tf
), for the final state in Eq. (40) with no changes in the other weight matrices. By applying this large weight matrix, F(
tf
), LQT control with the free-final state could produce the result using fixed-final state control (Lewis and Syrmos, 1995).
Figure 13
a shows the
Straight-line V-bar approach.
Position tracking error history.
position tracking error during the entire simulation. While the approach trajectory using the SDRE controller shows a constant straight-line trajectory whose error is less than 0.02 m, the one resulting from LQT control with the free-final
Approach range rate history.
Performance in translational maneuvering in the final straight-line V-bar approach
state shows a variation from about 0.14 m to 4×10
^{-3}
m at the final time, as shown in
Fig. 13
b. The approach range rate using SDRE control, shown in
Fig. 14
a, converged to 0.0077 m/s, whereas the approach range rate using LQT control with the free-final state, shown in
Fig. 14
b, dropped to 0.0049 m/s. The performance regarding translational maneuvering in the final approach is listed in
Table 1
.
The control force histories produced by the two controllers are shown in
Fig. 15
. The figure also shows that the additional control forces that resulted from the readjustment
Control force histories
Propellant mass consumption
of the weight matrices in the SDRE controller increased impulsively before the final straight-line approach. LQT control with the free-final state shows that the chaser spacecraft more precisely approached each destination(shown in
Fig. 1
as S2, S3, and the docking port) than the SDRE controller did. The translational maneuvers using the SDRE controller by adjusting the weight matrices in the straight-line LOS final approach could track the reference state provided by the straight-line guidance scheme and avoid thruster saturation. The tracking error was effectively reduced. However, additional control forces were required in straight-line maneuvering using the SDRE controller. Unlike LQT control with the free-final state, the effects of relative perturbations to the nonlinear system were considered in the SDRE approach in an attempt to control the translational maneuvering more precisely. For more precise translational maneuvers, the SDRE controller using nonlinear relative motion dynamics including relative perturbations can be closer to true system control than LQT control with the freefinal state. As the control forces and torques are generated,propellant mass is consumed, resulting in variations in the total mass and moment of inertia of the chaser. Figure 16 shows the propellant mass consumption as a result of the applied control forces and control torques. This study assumed that the locations of all RCSs were known so that the propellant mass consumed by the control torques could be computed. The total is the sum of the mass consumed by both the control forces and the control torques. This study
Targe t and chaser quaternion histories
Euler angle histories.
assumed that the locations of all RCSs were known so that the propellant mass consumed by the control torques could be computed. The total is the sum of the mass consumed by both the control forces and the control torques.
Euler angle error history
Angular rate histories
Euler angle histories of the target and the chaser expressed with the 3-2-1 rotational sequence of the chaser, as converted from quaternions, are shown in
Fig. 18
.
Figure 19 shows the Euler angle error history between the target and the chaser spacecraft for axial alignments that are composed of 30 and 90 degrees of pitch rotation. The target is orbiting in a near-circular 350 × 450 km altitude Earthpointing orbit. The target angular rate can then be expressed as
ω_{t}
=[0
n
0]
^{T}
in the adopted body-frame, where n is the mean motion of the target and its value is 0.065 deg/s. This attitude control can be achieved by using the chaser’s precise estimation of the attitude with the sensors previously listed and by adjusting the weight matrices. The chaser spacecraft then maintains the desired attitude by using the control torque continuously until the terminal time, leading to roll, pitch, and yaw errors of less than 0.1 degrees. Figure 20 shows the nominal target angular rate that the chaser can track for the two different types of axial alignment, which are illustrated in
Fig. 4
. For effective and rapid axial alignment,the weight matrix was adjusted during the steady-state interval, as was done for SDRE control. Even in the presence of unmodeled disturbance torques in the controller and uncertainty regarding the moment of inertia, the controller
Gravity-gradient and external torque histories
Control torque history.
succeeded in maintaining an acceptable attitude tracking error. The rotational maneuvers thus could meet the alignment condition for the docking phase. The performance regarding rotational maneuver in the final straight-line V-bar approach is listed in
Table 2
.
Performance regarding rotational maneuver in the finalstraight-line V-bar approach LQG: linear quadratic tracking.
The upper plot in
Fig. 21
shows the gravity-gradient torque history acting on the chaser, and the lower plot shows the external disturbance torque history acting on the chaser,which was simulated by white Gaussian-noise with mean[10 10 10]T, both of which were added to Eq. (19). Figure 22 shows the applied control torque history. Along with the two required attitude changes, the large, initial, disturbing control torques were applied and the response was then reduced to nearly zero.

Optimal control
;
Proximity operations and docking
;
State-dependent Riccati equation controller
;
Linear quadratic tracking controller
;
Linear quadratic Gaussian controller

1. Introduction

Autonomous rendezvous and docking of vehicles are essential for future space exploration of the Moon, Mars, and beyond, and for the supply and repair of vehicles such as the International Space Station (ISS). Spacecraft rendezvous and docking dates back to the manned US Gemini and Apollo programs and the unmanned Russian Cosmos missions of the late 1960s. In all the US missions, the human pilots were in the vehicle control loop during rendezvous and docking,whereas the Russians used a primarily automated approach,giving the pilots a supervisory role. The Gemini and Apollo programs developed the initial concept for rendezvous and capture. The Shuttle has demonstrated that these can be performed for various low Earth orbit (LEO) missions.Recent programs, such as Demonstration of Autonomous Rendezvous Technology (DART) (Rumford, 2002; Zimpfer et al., 2005), XSS-11 (Zimpfer et al., 2005), Orbital Express(Zimpfer et al., 2005), Automatic Transfer Vehicle (ATV)(Gonnaud and Pascal, 1999; Zimpfer et al., 2005), H-2 Transfer Vehicle (HTV) (Zimpfer et al., 2005), and the Hubble Robotic Servicing and Deorbit Mission (HRSDM)(Zimpfer et al., 2005), have been proposed and developed to demonstrate many of the technologies required for exploration that entails rendezvous and docking. Such programs have spurred significant developments in the area of autonomous rendezvous and capture (AR&C). However,the DART mission in 2005 experienced some failures.A 2006 National Aeronautics and Space Administration(NASA) report (NASA, 2006) provides an overview of the DART mishap investigation. The lessons learned from the mishap will hopefully facilitate the future development of autonomous capabilities.
Proximity operations and docking are critical phases of a rendezvous mission because highly precise translational and rotational maneuvers are required. AR&C capabilities will continue to be important for the successful execution of the various phases of a typical rendezvous mission, including the homing phase, the closing phase, proximity operations,and the final translational approach. The last two phases are critical from the perspective of mission completion and safety (Wertz and Bell, 2003). These phases are characterized by a chaser undergoing controlled motion that tracks a predetermined reference trajectory toward the docking port of the target. Like proximity operations, docking maneuvers have been accomplished in manned missions since the early days of the space program. Autonomous proximity operations are required for many future missions; however,further research and development is necessary if they are to be accomplished routinely. These requirements have led to the development and evaluation of several relative navigation sensors, such as the video guidance sensor (VGS),the global positioning system (GPS), light detection and ranging (LiDars) sensors, the laser dynamics range imager(LDRI), optical sensors, star trackers, inertial navigation systems (INS), and others. Any of these sensors can provide the basis for accurate relative positions and attitudes.Various navigation sensors are used to estimate the relative state information for feedback to an automated rendezvous and docking operations controller. Of course, autonomous spacecraft rendezvous and docking requires highly precise translational and rotational maneuvers that are integrated with precise sensors.
The motion of the chaser with respect to the target spacecraft is modeled by nonlinear relative equations of motion including the oblateness and aerodynamic drag.The attitude dynamics are modeled by the well-known Euler rotational equations of motion, including the gravitygradient torque and attitude kinematic equations using quaternions. Optimal control techniques are applied to the Shuttle orbiter’s manual phase flight segment covering proximity operations and docking, as shown in Fig. 1. The Shuttle program is scheduled to be retired by 2010; however,the Shuttle orbiter’s standard ISS approach technique may be extended or applied to other programs. The intent here is to develop suitable control algorithms that can facilitate autonomous proximity operations and docking. The trajectory profile design in the manual flight phase of the Shuttle is highly dependent on the payload configuration (Olszewski,1990). Payload attitude control and susceptibility to plume impingement are primary drivers in the final approach. The preferred technique for preventing plume impingement is the use of the V-bar guidance method during the final approach (Pearson, 1989), which guides the approach of the Space Shuttle along a velocity vector toward a target such as the ISS. This well-known approach is a form of pursuit guidance, and has been thoroughly investigated dating back to the Gemini program (Pearson, 1989). To initiate the V-bar approach, the active vehicle nulls the orthogonal relative velocities along the V-bar direction and accelerates to the desired closing rate along the V-bar direction, which is now along the line of sight to the target. The V-bar final approach is desirable because it is relatively fuel efficient. In addition, the constant orientation of the earth’s horizon provides a good reference for piloting, and closing rates can be easily and immediately nulled with subsequent station-keeping should some Shuttle or payload system anomaly occur (Pearson,1989). The trajectory of the Shuttle to the ISS in use since the STS-102 mission in 2001 is by default the starting point for the design of the lower surface-inspection maneuver. Figure 1 illustrates the following approach trajectory in the rotating local vertical local horizontal (LVLH) frame, centered at the ISS center of mass.
Lager Image

2. Dynamic Modeling

- 2.1 Coordinate systems

The coordinate systems used are LVLH frames centered on the target and the chaser spacecraft, and an orthogonal body-fixed frame at the center of mass and the earthcentered inertial (ECI) frame, N, as shown in
Fig. .2
The LVLH frame is sometimes referred to as the CW frame (Fehse, 2003;Prussing and Conway, 1993; Schaub and Junkins, 2003), E,with the x-axis being directed radially outward along the local vertical, the y-axis along the direction of motion or velocity direction, and the z-axis normal to the reference orbit plane.In a rendezvous mission, the motion of the chaser spacecraft is commonly described relative to the target spacecraft.Instead of the CW frame E, the spacecraft local orbital frame A is adopted here to describe its motion. This frame is related to the CW frame E.
Lager Image

Lager Image

- 2.2 Relative motion dynamics

Two sets of equations of motion describing the chaser motion relative to the target in the LVLH frame are defined here. The first set includes the exact nonlinear relative equations of motion including the relative oblateness effect,J2, and the aerodynamic drag. The second set includes the linear relative equations of motions with no perturbations.The first set was used for the SDRE control formulation. The second set was used for the LQT control formulation with the free-final state. Among the many sources of perturbations,the earth’s oblateness and aerodynamic drag in the LEO are dominant. The relative effects of Earth’s oblateness and aerodynamic drag are included in the exact nonlinear relative equations for more precise dynamic modeling. In the CW frame, E, the perturbing acceleration due to oblateness is described by Prussing and Conway (1993):
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 2.3 Attitude dynamics and kinematics

The rotational motion of the chaser expressed in the bodyfixed frame is described by the well-known Euler’s equations of motion. Like the perturbing accelerations in relative motion dynamics, the attitude dynamics also experience disturbing torques such as the torque due to aerodynamic drag, magnetic field torque, and gravity-gradient torque due to the asymmetry of the spacecraft. This study models only the gravity-gradient torque. The gravity-gradient torque due to asymmetry of the body, expressed using the local orbital frame, A, is given in vector/dyadic form as (Wie, 1998):
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

3. Control Law Formulation

In this section, three optimal control techniques are derived for translational and rotational maneuvering. The SDRE tracking controller used with the nonlinear relative motion dynamics and the LQT controller with the free-final state used with the linearized relative dynamics are both derived for translational maneuvers. By using thrusters for translational and rotational control, both can be uncoupled to a high degree of approximation. However, a disturbance torque generated by thruster firing is considered as an unmodeled disturbance torque. The two controllers have different purposes but are executed simultaneously. The state vector for the translational maneuvers is given by:
Lager Image

- 3.1 The SDRE tracking formulation for translational maneuvers

Motivated by the LQR method, which is characterized by the solution to the algebraic Riccati equation (ARE),SDRE feedback control is an extended linearization (SDC)control method that provides an approach similar to that of a nonlinear regulation problem. For reference trajectory tracking, the regulator problem must be recast as a tracking problem. The goal is to drive the error between the reference and the output to zero with minimum control energy. The tracking problem is then formulated with the performance index as follows.
Lager Image

Lager Image

Lager Image

Lager Image

- i) Use the direct parameter method to bring the nonlinear equation into the SDC form as in Eq. (21).
- ii) Solve the SDRE to obtain P(x)≥0,

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 3.2 LQT formulation for relative translational motion

Tracking systems need to track a desired trajectory in some optimal sense. LQT control with the free-final state is developed to maintain the output as closely as possible to the desired output with minimum control acceleration. A linear, observable system from Eq. (9) is given as:
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 3.3 LQG formulation for rotational maneuvers

This section describes some properties of quaternions that make it possible to realize an exact linearization of the error dynamics formulation. This study uses this linearized equation to take advantage of the simplified equation of motion to determine the precise attitude control. The control law formulation using LQR is combined with the extended Kalman filter (EKF), which leads to an LQG-type control system. The state vector for attitude control is given by:
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

4. Numerical Results and Analysis

This study tested the controllers developed here in a Shuttle-like scenario in which the proximity operations and docking have historically been performed by the crew in the manual phase flight segment, as shown in
Fig. 1.
A six degrees-of-freedom simulation and a passive target (ISS) were created to demonstrate the performance of the controllers. The ultimate objective was to have the chaser (Shuttle) docking port approach the target docking port leading to a soft docking with the desired attitude. There were two types of axial alignment during this phase. The first one was the simplest case in which both vehicles had the same direction of axial alignment, as illustrated by
Figs. 4
a and b. The second one was the case in which the chaser had to execute a pitch rotation of +90 degrees, as illustrated by
Figs. 4
a and c.
Lager Image

Lager Image

Lager Image

- 4.1 Results of translational maneuvering

Figures 5
and
6
show proximity operation trajectories andout-of plane motions for the entire duration of simulation. Using the CW guidance scheme and final straight-lineapproach guidance with subcentimeter-per-secondvelocities, the reference trajectories were autonomouslydirected to the controllers. The transfer time and the nominaltrajectory for the R-bar relative motion were determined bythe CW terminal guidance scheme (Fehse, 2003; Prussingand Conway, 1993).
Figures 7
and
8
show the fly-aroundmaneuver from S1 to S2 for 13 minutes, and the fly-aroundmaneuver from S2 to S3 for 10 minutes, respectively. As soonas each sub-segment was complete, the newly initializedcontrol was applied sequentially to the next sub-segment.
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Performance in translational maneuvering in the final straight-line V-bar approach

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 4.2 Rotational maneuver

Because of varying mass, the moment of inertia is not constant. The robustness of the rotational controller was evaluated for uncertainties in the moment of inertia. The uncertainty was quantified by adding 30 percent of the initial moment of inertia to this initial value (Xin et al., 2004).Moreover, external disturbances, which were not modeled in the design of the controller, were added to the Euler rotational equation of the chaser in Eq. (19). All of the results to follow were compiled after incorporating the moment of inertia uncertainty and the external disturbances. The first rotational maneuver, illustrated by a comparison of
Figs. 4
a and b, was executed to align the chaser body axis with the target body axis in the first segment. The second rotational maneuver,illustrated by comparing
Figs. 4
a and c, was performed up to the final sub-segment after the first rotational maneuver was complete.
Figure 17
shows that the chaser exhibited good attitude tracking during the entire simulation. To more clearly illustrate the geometry of the angular motion, the
Lager Image

Lager Image

Lager Image

Lager Image

Performance regarding rotational maneuver in the finalstraight-line V-bar approachLQG: linear quadratic tracking.

Lager Image

Conclusions

This study designed a nonlinear SDRE control technique and an LQT control technique for translational maneuvers between spacecraft. The adopted SDRE control technique was designed for nonlinear relative motion dynamics, including Earth’s oblateness and aerodynamic drag perturbations.Through the SDC parameterization, a linear-like closedform structure was achieved for the nonlinear system control problem. For the tracking command, the controller was designed without increasing the state dimensions, unlike the SDRE integral servo controller formulation. The tracking results using SDRE control were successfully achieved by readjusting the weight matrices at the steady-state interval to avoid saturation. The LQT control designed for the linearized system shows that the free-final state becomes the fixed final state when the weight matrix for the final state is large. This result is very effective in meeting the required final state condition. The LQT control results can be used as a reference in designing the SDRE controller. However, for more precise control, SDRE control is preferable because the controller was designed for the nonlinear system including the relative perturbations. An LQG-type control technique was designed for the rotational maneuvers. A combination of an LQR controller and EKF estimation using star trackers and threeaxis gyros constitutes the LQG-type control technique. The weight matrices were readjusted similarly for the LQG-type control as for SDRE control to decrease the attitude tracking error to within the desired accuracy. The tracking error was maintained even in the presence of disturbance torques and uncertainties in the moments of inertia, demonstrating the robustness of this LQG-type controller. The Shuttle crew’s manual flight segment was chosen to simulate and evaluate the control techniques applied in this study. A six degrees-of freedom simulation demonstrated that the adopted control techniques can successfully conduct proximity operations and meet the conditions for the docking phase. The control techniques can also be applied to other programs that require autonomous proximity operations composed of many subphases,including docking. For executing autonomous and precise proximity operations, including successful docking,the integration of these controllers with accurate state estimation using high-precision sensors is needed.
Alba-Flores R
,
Barbieri E
2006
Real-time infinite horizon linear-quadratic tracking controller for vibration quenching in flexible beams
IEEE International Conference on Systems Man and Cybernetics
Taipei Taiwan
38 -
43

Bach R
,
Paielli R
1993
Linearization of attitudecontrol error dynamics
IEEE Transactions on Automatic Control
38
1521 -
1525
** DOI : 10.1109/9.241567**

Budiyono A
,
Wibowo S. S
2007
Optimal tracking controller design for a small scale helicopter
Journal of Bionic Engineering
4
271 -
280
** DOI : 10.1016/S1672-6529(07)60041-9**

Cimen T
2008
State-Dependent Riccati Equation(SDRE) control: a survey
Proceedings of the 17th International Federation of Automatic Control (IFAC) World Congress
Seoul Korea
3761 -
3775

Cloutier J. R
1997
State-dependent Riccati equation techniques: an overview
Proceedings of the American Control Conference
Albuquerque NM
932 -
936

Crassidis J. L
,
Junkins J. L
2004
Optimal Estimation of Dynamic Systems
Chapman & Hall/CRC
Boca Raton

Fehse W
2003
Automated Rendezvous and Docking of Spacecraft
Cambridge Univeristy Press
Cambridge

Gonnaud J. L
,
Pascal V
1999
ATV guidancenavigation and control for rendezvous with ISS.
Proceedings of the 4th ESA International Conference on Spacecraft Guidance Navigation and Control Systems
Noordwijk the Netherlands
501 -

Lee D
,
Pernicka H
2009
Vision-based relative stateestimation using the unscented Kalman filter. Spaceflight Mechanics Meeting
Savannah GA
Paper No. 09-166

Lefferts E. J
,
Markley F. L
,
Shuster M. D
1993
Kalman filtering for spacecraft attitude estimation
Journal of Guidance Control and Dynamics
5
417 -
429

Lewis F. L
,
Syrmos V. L
1995
Optimal Control
2nd ed
Wiley
New York

Madonna R. G
1997
Orbital Mechanics
Original ed
Krieger Pub. Co
Malabar FL
105 -
107

Nagata T. J
,
Modi V
,
matsuo H
2001
Dynamics and control of flexible multibody systems: Part II: simulation code and parametric studies with nonlinear control
Acta Astronautica
49
595 -
610

Naidu D. S
2003
Optimal Control Systems
CRC Press
Boca Raton FL

National Aeronautics and Space Administration (NASA).
Overview of the DART mishap investigation results (for public release)
(May 15 2006)

Olszewski O. J
1990
Automated terminal guidance for a shuttle rendezvous to space station freedom
AIAA GuidanceNavigation and Control Conference
Portland OR
377 -
388

Paielli R. A
,
Bach R. E
1993
Attitude control with realization of linear error dynamics
Journal of GuidanceControl and Dynamics
16
182 -
189
** DOI : 10.2514/3.11444**

Pearson D. J
1989
Shuttle rendezvous and proximity operations
Proceedings of the CNES International Symposium on Space Dynamics
Paris
833 -
851

Prussing J. E
,
Conway B. A
1993
Orbital Mechanics
Oxford University Press
New York

Rumford T
2002
Demonstration of Autonomous Rendezvous Technology (DART) project summary
CoreTechnologies for Space Systems Conference Colorado SpringsCO.

Schaub H
,
Junkins J. L
2003
Analytical Mechanics of Space Systems
American Institute of Aeronautics and Astronautics
Reston VA

Stansbery D. T
,
Cloutier J. R
2000
Position and attitude control of a spacecraft using the state-dependent Riccati equation technique
Proceedings of the American Control Conference
Chicago IL
1867 -
1871

Sznaier M
,
Suarez R
2001
Suboptimal control of constrained nonlinear systems via receding horizon state dependent Riccati equations
Proceedings of the 40th IEEE Conference on Decision and Control
Orlando FL
3832 -
3837

Vallado D. A
,
McClain W. D
2001
Analytical Mechanics of Space Systems
2nd ed
Kluwer Academic Publishers
Boston

Walker S. R
,
LoPresti J. A
,
Schrock M. B
2005
Space shuttle Rbar pitch maneuver for thermal protection system inspection
AIAA Guidance Navigation and Control Conference and Exhibit
San Francisco CA
AIAA 2005-5983

Wertz J. R
,
Bell R
2003
Autonomous rendezvous and docking technologies-status and prospects
Space Systems Technology and Operations Conference SPIE AeroSense Symposium
Orlando FL
Paper No. 5088-5083

Wie B
1998
Fundamentals of Astrodynamics and Applications
American Institute of Aeronautics and Astronautics
RestonVA
365 -
369

Xin M
,
Balakrishnan S. N
,
Stansbery D. T
2004
Spacecraft position and attitude control with θ-D technique
42nd AIAA Aerospace Sciences Meeting and Exhibit
RenoNV
AIAA 2004-2540

Zimpfer D
,
Kachmar P
,
Tuohy S
2005
Autonomous rendezvous capture and in-space assembly: past present and future
1st Space Exploration Conference: Continuing the Voyage of Discovery
Orlando FL
AIAA 2005-2523

Citing 'Optimal Control for Proximity Operations and Docking
'

@article{ HGJHC0_2010_v11n3_206}
,title={Optimal Control for Proximity Operations and Docking}
,volume={3}
, url={http://dx.doi.org/10.5139/IJASS.2010.11.3.206}, DOI={10.5139/IJASS.2010.11.3.206}
, number= {3}
, journal={International Journal of Aeronautical and Space Sciences}
, publisher={The Korean Society for Aeronautical & Space Sciences}
, author={Lee, Daero
and
Pernicka, Henry}
, year={2010}
, month={Sep}