Spacecraft Formation Reconfiguration using Impulsive Control Input
 Author: Bae Jonghee, Kim Youdan
 Organization: Bae Jonghee; Kim Youdan
 Publish: International Journal Aeronautical and Space Sciences Volume 14, Issue2, p183~192, 30 June 2013

ABSTRACT
This paper presents formation reconfiguration using impulsive control input for spacecraft formation flying. Spacecraft in a formation should change the formation size and/or geometry according to the mission requirements and space environment. To modify the formation radius and geometry with respect to the leader spacecraft, the follower spacecraft generates additional control inputs; the two impulsive control inputs are general control type of the spacecraft system. For the impulsive control input, Lambert’s problem is modified to construct the transfer orbit in relative motion, given two position vectors at the initial and final time. Moreover, the numerical simulation results show the transfer trajectories to resize the formation radius in the radial/alongtrack plane formation and in the alongtrack/crosstrack plane formation. In addition, the maneuver characteristics are described by comparing the differential orbital elements between the reference orbit and transfer orbit in the radial/alongtrack plane formation and alongtrack/crosstrack plane formation.

KEYWORD
Spacecraft Formation Flying , Formation Reconfiguration , Impulsive Control Input , Lambert’s Problem

1. Introduction
Spacecraft formation flying has been widely investigated due to increasing interests in the design of clusters and their application. Multiple spacecraft in formation have many advantages in terms of flexibility, reliability, financial benefits, and high image resolution compared to a single large spacecraft. Autonomous formation flying has been developed for future space missions in many countries. National Aeronautics and Space Administration (NASA) provided the EarthObserving 1 (EO1), StarLight (ST) series, Terrestrial Planet Finder (TPF), MicroArcsecond Xray Imaging Mission (MAXIM), and Magnetospheric Multiscale (MMS) mission using multiple spacecraft [14]. Moreover, European Space Agency (ESA) developed the Cluster II mission, the Dawin, the Xray Evolving Universe Spectroscopy (XEUS) missions, and PROBA3 consisting of multiple spacecraft for various missions [5,6]. The Laser Interferometer Space Antenna (LISA) was designed to develop a spacebased gravitational wave detector by NASA and ESA. During the spacecraft operation, formation reconfiguration may be needed to change the formation radius and/or geometry according to the mission requirements. For formation reconfiguration, the follower spacecraft should modify the orbit with respect to the leader spacecraft by additional control inputs.
Lots of works on the maneuver problem have been solved in an inertial frame. The goals of the traditional maneuver problem are to change the orbit size or shape of a single spacecraft with respect to a planet considering mission requirements, and to minimize the control effort during the maneuver. For spacecraft formation, on the other hand, the formation can be reconfigured by increasing or decreasing the formation radius and/or changing the formation geometry through maneuver in relative motion with respect to the leader spacecraft. This transfer between two orbits can be generally performed by impulsive control input. The traditional maneuver problem for spacecraft orbit transfer has been solved by two impulsive control inputs. The Hohmann transfer is well known as twoburn minimum fuel transfer in circular coplanar orbits [7,8]. For the general twoimpulse maneuver problem between two fixed position vectors with the specified flight time, the transfer orbit can be determined by solving Lambert's problem [810], and therefore the required velocity vectors of the resulting transfer orbit can be calculated at the initial and final time. Many researchers utilized Lambert's problem to obtain the transfer trajectory between two orbits with respect to a planet [1113]. Won considered minimumfuel and minimumtime orbit transfer problem in two coplanar elliptic orbits [11]. Prussing considered the minimumfuel impulsive spacecraft trajectory through the multiplerevolution [12], and Shen and Tsiotras studied the optimal fixedtime, twoimpulse rendezvous problem between two satellites in coplanar circular orbit using multiplerevolution Lambert's problem [13].
Although the traditional maneuver problem has been widely studied, the formation reconfiguration problem is different from that because the formation of multiple spacecraft is expressed by the relative motion dynamics with respect to the leader spacecraft, not a planet. Vaddi et al. studied the formation reconfiguration problem using impulsive control based on Gauss's variational equation [14], and Ketema considered the optimal transfer problem using two impulsive control inputs considering the relative motion of the spacecraft formation [15]. However, the transfer method for the follower spacecraft from one orbit to another with respect to the leader spacecraft in the local frame, where the initial and final positions, velocities, transfer time, and orbital elements of the leader spacecraft are specified, has not been sufficiently investigated. In this study, the traditional spacecraft maneuver method, Lambert’s problem, is modified to design the transfer orbit trajectory of the follower spacecraft, where two relative position vectors and the flight time in the local frame are given. Furthermore, the characteristics of the formation reconfiguration in two formation geometries, radial/alongtrack plane formation and alongtrack/crosstrack plane formation, are analyzed in terms of the orbital elements. The numerical simulation results show that the transfer trajectories are designed to resize the formation radius in two formation geometries.
This paper is organized as follows. In Section 2, the relative motion dynamics are described. Spacecraft formation reconfiguration using impulsive control input is derived in Section 3, and the numerical simulation results and analysis to verify the formation reconfiguration are shown in Section 4. Finally, conclusion is presented in Section 5.
2. Spacecraft Relative Motion Dynamics
The inertial coordinate system and a rotating local coordinate system are defined for spacecraft formation flying. As shown in Fig. 1, the inertial coordinate system {
X, Y, Z } is attached to the center of the Earth, and a rotating localverticallocalhorizontal (LVLH) Cartesian frame {x, y, z } is attached to the center of the leader spacecraft. Thex axis points in the radial direction, they axis points in the alongtrack direction, and thez axis is perpendicular to the orbital plane and points in the direction of the angular momentum vector.denotes the position vector of the leader spacecraft and
is the relative position vector of the follower spacecraft with respect to the leader spacecraft in the LVLH frame. Then, the position of each follower spacecraft within a formation is given by
The nonlinear equation of the relative motion for elliptic reference orbits can be represented as
where
is the angular velocity, and
are external forces acting on the leader and follower spacecraft, respectively. When the external forces only act on a central gravitational field, the nonlinear relative dynamics for elliptic reference orbits can be represented as follows:
where
μ is the gravitational coefficient,θ is the argument of latitude of the leader spacecraft,is the orbital angular velocity, and
u_{x} ,u_{y} , andu_{z} are the control inputs in theradial, alongtrack, and crosstrack direction, respectively. In orbital mechanics, the radius and angular velocity of the leader spacecraft are given as [16]
where
a is the semimajor axis, e is the eccentricity,v is the true anomaly, andis the natural frequency of the reference orbit.
Equation (3) can be expressed using the true anomaly
v as a free variable, becausev of the reference orbit monotonically increases and describes the angular velocity. The following transformation is considered:where ( )' and ( )'' denote the first and second derivatives with respect to
v . With the transformation of Eq. (5) and the assumption thatthe relative motion equation can be expressed as follows:
Using Eq. (7), the maneuver problem for the formation reconfiguration is solved and the initial and final conditions are selected in the true anomaly domain to meet the periodicity condition as presented in Ref. [17].
3. Formation Reconfiguration using Impul sive Control Input
3.1 Lambert’s Problem
Lambert’s theorem states that the flight time from one point to another depends on the semimajor axis of transfer orbit, the distance of the initial and final points of the arc from the center of force, and the length of the chord joining the points. The chord length,
c , between two position vectors, as shown in Fig. 2, is defined by the cosine law as follows:
In Fig. 2,
F is the primary focus which is the location of the Earth,are the position vectors at the initial time
t _{1} and the final timet _{2}, andθ_{T} is the transfer angle between two position vectors defined byThe direction of flight can be defined by the transfer angle; the spacecraft moves along the short way when
θ_{T} ≤π ,, and it moves along the long way whenθ_{T} ？π .To determine the transfer orbit between two position vectors, two foci and semimajor axis should be found. From the definition of the ellipse, the location of the vacant focus and semimajor axis of the transfer orbit can be determined. As shown in Fig. 2, two circles can be drawn having 2
ar _{1} and 2ar _{2} as radii and pointsP _{1} andP _{2} as centers, respectively, since the sum of the distance from the foci to any point on the ellipse equals twice the semimajor axis as 2a =r _{1}+(2ar _{1})=r _{2}+(2ar _{2}). As a result, two intersected points,F ' andF '', are determined as shown in Fig. 2, which are the candidate locations of other focus of the transfer orbit. Using two vacant foci, two elliptic orbits can be obtained as shown in Fig. 2. Note from Fig. 2 that two ellipses have different eccentricities; the ellipse with vacant focusF ' has the smaller eccentricity, and the ellipse with F'' has the larger eccentricity. The location of vacant focus depends on the semimajor axis of the transfer orbit, and therefore the semimajor axis should be determined carefully through iterative calculation for the minimum energy.When the transfer orbit is determined, the velocity vectors of the transfer orbot at the initial and final time,
, as shown in Fig. 3, can be calculation by Kepler's equation. Consequently, the velocity changes Δ
v _{1} and Δv _{2} att _{1} andt _{2} can be obtained aswhere
are velocity vectors of the spacecraft around a planet at the initial and final time. Figure 3 shows the velocity changes using two impulse inputs. The total velocity change can defined as follows:
3.2 Modified Lambert’s Problem for Follower Spacecraft
The transfer orbit for the follower spacecraft connecting two relative position vectors can be solved modifying Lambert’s problem, when the orbital elements of the leader spacecraft and the flight time are known. The geometry of the transfer orbit for the follower spacecraft is illustrated in Fig. 4.
The relative position vectors of the follower spacecraft with respect to the leader spacecraft,
ath the initial time
t _{1} and the final timet _{2} in the LVLH frame, respectively, can be expressed as follows:The position vectors of the follower spacecraft with respect to the Earth at
t _{1} andt _{2} in the ECI (Earth Centered Inertial) frame can be written aswhere
are the position vectors of the leader spacecraft in the ECI frame at
t _{1} andt _{2}, respectively.Let us define the velocity vectors of the follower spacecraft at
t _{1} andt _{2} in the ECI frame as, respectively. Given two position vectors,
the transfer orbit for the follower spacecraft can be determined through Lambert’s problem as explained in the previous section. The transfer angle can be defined by
The chord length is determined as
Using
θ_{T} andc , Lambert’s problem can provide the velocities of the transfer orbit,as shown in Fig. 4. The transfer time for the elliptic orbit is defined from Kepler’s equation as follows [7,8]:
where Δ
t =t _{2}t _{1}, andE _{1} andE _{2} are eccentric anomalies att _{1} andt _{2}. In this study, multiple revolutions are not considered, and then 2πk term in Eq. (18) can be ignored. Equation (18) can be rewritten aswhere Δ
E =E _{2}E _{1}. To solve Kepler's equation, the universal variable,χ , is adopted in this study, which relates energy and angular momentum [8]. The universal variable,χ is defined as follows:In addition, parameters
s _{2} ands _{3} are defined as [8]Using Eqs. (20) and (21), Eq. (19) can be expressed as
With trigonometric identity, sin(
a _{1}a _{2})=sina _{1}cosa _{2} cosa _{1}sina _{2}, Eq. (22) can be rewritten asThe relation between eccentric anomaly
E and true anomalyv is defined by [7,8]After some manipulation with Eq. (24), Eq. (23) can be expressed as follows:
With the trigonometric identity, Eq. (25) is presented as
The radius is defined as
r =a (1e cosE ) in orbital mechanics, and then Eq. (26) can be rewritten aswhere
The resulting
A andY give the velocity vectors of the transfer orbit att _{1} andt _{2} as follows:where
Finally, the velocity changes Δ
V _{1} and ΔV _{2} for the follower spacecraft att _{1} andt _{2} can be obtained as follows:Consequently, the total velocity change can be calculated as
Note that
are velocities with respect to the Earth and expressed in the ECI frame. Thus, the relative velocity vectors of the follower spacecraft with respect to the leader spacecraft can be expressed as follows:
where
are velocity vectors of the leader spacecraft at
t _{1} andt _{2}, respectively, and ω is the angular velocity vector of the leader spacecraft as defined byFinally, the relative velocity vectors in the LVLH frame and the orbital elements of the transfer orbit can be found.
4. Numerical Simulation and Analysis
In this section, numerical simulations are performed for the maneuver problem of formation reconfiguration. Using impulsive control input, two maneuvering cases are considered to verify the formation reconfiguration in relative motion: change of the formation radius,
r _{d}, (i) in the radial/alongtrack plane formation (RAPF) and (ii) in the alongtrack/ crosstrack plane formation (ACPF).4.1 Radial/AlongTrack Plane Formation
In this section, the maneuver for formation reconfiguration in the RAPF is performed using impulsive control input. The orbital elements of the leader spacecraft are summarized in Table 1. To change the formation radius in the RAPF, we consider a situation that the follower spacecraft moves from the initial formation radius of 5,000 m to the final formation radius of 10,000 m. The initial and final conditions are selected to satisfy the periodicity and zero offset condition according to the formation geometry as presented in Ref. [17]. The initial and final conditions of the follower spacecraft in the LVLH frame are summarized in Table 2. The flight time is obtained as 5166.4 seconds given the initial and final conditions as described in Table 2.
Figures 5？7 show the simulation results for maneuvering in the RAPF. Figure 5 shows the trajectories of the leader and follower spacecraft in the ECI frame; the normal line denotes the trajectory of the leader spacecraft and the thick line transfer trajectory of the follower spacecraft. As shown in Fig. 5, it is difficult to understand the relative motion of the follower spacecraft with respect to the leader spacecraft in the ECI frame. Instead, Fig. 6 transfer trajectory of the follower spacecraft in the LVLH frame. As shown in Fig. 6, the follower spacecraft moves from the initial position to the final position on the
xy plane in the LVLH frame. Figure 7 describes the relative velocity history ofwer spacecraft during the maneuver. In Fig. 7, the triangle denotes specified relative velocities at the initial and final time in the LVLH frame given in Table 2. As a result, the difference between the resulting velocity from the modified Lambert’s problem and the specified velocity is equal to the required velocity change for the orbit transfer. Table 3 summarizes the relative velocities of the follower spacecraft and the velocity change, △V , for the reconfiguration maneuver in the RAPF.4.2 AlongTrack/CrossTrack Plane Formation
In this section, the maneuver for formation reconfiguration in the ACPF is performed using impulsive control input. The orbital elements of the leader spacecraft are same as previous section. To change the formation radius in the ACPF, we consider a situation in which the follower spacecraft moves
from the initial formation radius of 5,000 m to the final formation radius of 10,000 m. The initial and final conditions of the follower spacecraft in the LVLH frame are given by
formation geometry in Table 4 as presented in Ref. [17].
Figures 8？9 show the simulation results for formation reconfiguration using impulsive control input in the ACPF. Figure 8 shows the transfer trajectory of the follower spacecraft with respect to the leader spacecraft in the LVLH frame. Figure 9 shows the velocity history of the follower spacecraft during the maneuver. In Fig. 9, the triangle denotes specified velocities at the initial and final time in the LVLH frame given in Table 4. Therefore, the required velocity change can be obtained by computing the difference between the resulting velocity from the modified Lambert’s problem and the specified velocity. Table 5 summarizes the velocities of the follower spacecraft and △
V for the maneuver in the ACPF.4.3 Analysis
In this section, the simulation results for the formation reconfiguration using impulsive control input in the RAPF and ACPF are analyzed. The simulation results are investigated in terms of the difference of the orbital elements between the transfer orbit and reference orbit with respect to the formation geometry.
As described in Sections 4.1 and 4.2, the transfer orbits are generated where the initial and final positions are arbitrarily chosen in the RAPF and ACPF. By comparing with the orbital elements of the leader spacecraft in Table 1, it is seen that the transfer orbits of two formation reconfigurations have almost same orbital elements, as described in Fig. 5; however, some orbital elements are different. Table 6 summarizes the differences of the orbital elements: (i) between the reference orbit and the transfer orbit in the RAPF as described in Sections 4.1, and (ii) between the reference orbit and the transfer orbit in the ACPF in Section 4.2. These results show the maneuver characteristics of relative motion in the RAPF and ACPF. As shown in Table 6, the difference of the semimajor axis for the ACPF case is much larger than that for the RAPF case. This difference leads to the total velocity change △
V for the formation reconfiguration; it can be seen in Tables 3 and 6 that △V for resizing the formation radius in the ACPF is increased compared with △V in the RAPF as △a in the ACPF becomes large. Furthermore, as represented in Table 6, the transferbetween two orbits in RAPF produces a difference in orbital elements including the semimajor axis
a , eccentricitye , and augment of perigeeω , which means that orbit transfer in the RAPF belongs to the coplanar maneuver. On the other hand, in the ACPF, not onlya ,e , andω , but also inclination i and longitude of the ascending node Ω are changed during the maneuver, as shown in Table 6. It means that the orbit transfer in the ACPF is related to the combined maneuver of coplanar and noncoplanar maneuvers.5. Conclusion
A formation reconfiguration using impulsive control input is developed for spacecraft formation flying. To transfer between two position vectors at initial and final time, two impulsive control inputs are generated through modification of the classical Lambert’s problem. The modified method provides the transfer orbit in relative motion for spacecraft formation flying, and therefore the follower spacecraft can change the formation radius and/or geometry with respect to the leader spacecraft. Numerical simulations are performed to solve the maneuver problem for formation reconfiguration in the RAPF and ACPF. These results show the maneuver characteristics of the transfer orbit with respect to the formation type. Through these result, it can be concluded that the transfer orbit for resizing formation radius in the RAPF is related to the coplanar maneuver, and transfer orbit in the ACPF belongs to the combined maneuver of coplanar and noncoplanar maneuvers.

[Fig. 1.] Local Vertical Local Horizontal frame

[Fig. 2.] Geometry for Lambert’s proble

[Fig. 3.] Orbital transfer using two impulsive control inputs

[Fig. 4.] Orbit transfer for follower spacecraft

[Table 1.] Orbital elements of leader spacecraft

[Table 2.] Initial and final condition (RAPF)

[Fig. 5.] Trajectory of two spacecraft in ECI frame (RAPF)

[Fig. 6.] Trajectory of follower spacecraft in LVLH frame (RAPF)

[Fig. 7.] Velocity history of follower spacecraft (RAPF)

[Table 3.] Solution of modified Lambert’s problem (RAPF)

[Table 4.] Initial and final condition (ACPF)

[Fig. 8.] Trajectory of follower spacecraft in LVLH frame (ACPF)

[Fig. 9.] Velocity history of follower spacecraft (RAPF)

[Table 5.] Solution of modified Lambert’s problem (ACPF)

[Table 6.] Difference of orbital elements