Abstract

The determination of minimum-fuel or minimum-time relative orbit trajectories represents a classical topic in astrodynamics. This work illustrates some symmetry properties that hold for optimal relative paths and can considerably simplify their determination. The existence of symmetry properties is demonstrated in the presence of certain boundary conditions for the problems of interest, described by the linear Euler-Hill-Clohessy-Wiltshire equations of relative motion. With regard to minimum-fuel paths, the primer vector theory predicts the existence of several powered phases, divided by coast arcs. In general, the optimal thrust sequence and duration depend on the time evolution of the switching function. In contrast, a minimum-time trajectory is composed of a single continuous-thrust phase. The first symmetry property concerns minimum-fuel and minimum-time orbit paths, both in two and in three dimensions. The second symmetry property regards minimum-fuel relative trajectories. Several examples illustrate the usefulness of these properties in determining minimum-time and minimum-fuel relative paths.

1. Introduction

Minimum-time and minimum-propellant-consumption space trajectories have represented a subject of great relevance in astrodynamics since the 1950s. Pioneering works in this field are due to Lawden [1], Leitmann [2], Miele [3], and Edelbaum [4]. More recently, Marec [5], Vinh [6], and Prussing [7, 8] provided relevant contributions, with reference to different types of aerospace missions.

Specifically, relative orbital motion has attracted an increasing interest in the last decades, both for formation flying design and for the purpose of planning spacecraft proximity maneuvers. A very useful mathematical framework for investigating relative orbital motion is represented by the Euler-Hill-Clohessy-Wiltshire (EHCW) equations [9]. They derive from the linear expansion of the dynamics equations that govern position and velocity of a spacecraft relative to a real or virtual vehicle, placed in a neighboring reference circular orbit. In the 1960s Tschauner and Hempel [10] addressed the problem of relative orbital motion with respect to a reference elliptic orbit, whereas Abrahamson and Stern [11] provided an extension to hyperbolic paths. Further works [12, 13] have been devoted to investigating relative orbital motion involving eccentric orbits, through the use of either the true or the eccentric anomaly as the independent variable. Most recently, Melton [14] developed an approximate solution that depends explicitly on time and uses a series expansion in terms of the reference orbit eccentricity. Schaub [15, 16] employed orbit element differences relative to a chief orbit with the intent of describing the relative orbit geometry and in order to establish relative orbits invariant with respect to the Earth oblateness. The coordinated motion of a set of spacecraft is a crucial issue for formation flying, and several works [1723] have been focused on this central topic, with the final aim of minimizing the control effort needed for relative configuration maintenance. However, orbital maneuvers may be needed both for formation flying disposal and, more generally, in the context of proximity maneuvers aimed at rendezvous or orbit rephasing. Minimizing the fuel consumption leads to extending the satellite operational life or to enhancing its operational capabilities. As a pioneering contribution in this field, Prussing [7, 8] concentrated on applying the primer vector theory to two-, three-, and four-impulse relative orbit maneuvers. Carter [24] considered finite-thrust paths, by employing an interesting geometric approach similar to that adopted by Prussing. More recent contributions either deal with constrained impulsive and finite-thrust maneuvers [2529] or are focused on formation flying reconfiguration [30].

In general, the problem of spacecraft trajectory optimization is to be formulated as an optimal control problem, which consists in finding the optimal thrust sequence, magnitude, and direction that minimize propellant consumption or time of flight. Only a few simplified problems are amenable to a closed-form solution, and therefore numerical techniques are usually employed in trajectory optimization problems. Classical (indirect and direct) methods need a first-attempt approximate guess to yield a refined solution, and this can be challenging to provide. A recently introduced methodology, termed indirect heuristic technique [31], employs the necessary conditions for optimality, that is, the Pontryagin minimum principle, the Euler-Lagrange equations, and the Weierstrass-Erdmann conditions, in conjunction with a heuristic optimizer. The latter requires the conversion of the optimal control problem into a parameter optimization problem in order to yield the desired solution. This goal is achieved by employing the necessary conditions for optimality, which allow expressing the control variables as functions of the adjoint variables conjugate to the dynamics equations. This circumstance has the consequence that the initial values of the adjoint variables are treated as the unknown parameters sought by the heuristic solver. In this context, a particular selection of values of the unknown parameters is associated with an individual. A population of individuals, whose number increases as the number of unknown parameters increases, is generated and evolves toward the optimal parameter set. Some recently published papers [32, 33] prove the effectiveness and accuracy of this approach, used also for the illustrative examples reported in this work.

This paper has the purpose of illustrating a pair of symmetry properties that hold for optimal relative space trajectories and can simplify their determination. In the mathematical framework established by the EHCW equations, this work is intended to(a)derive the analytical conditions related to minimum-fuel and minimum-time relative trajectories in a unified perspective,(b)prove a first symmetry property that holds for both minimum-fuel and minimum-time relative paths,(c)report the statement and proof of a second symmetry property related to the optimal thrust sequence and direction of minimum-fuel relative trajectories,(d)illustrate the use of both properties, also pointing out their specific utility.The first symmetry property was employed in [32], where only feasibility of symmetric paths was proven, but not their optimality. The second symmetry property was demonstrated and employed in [33] and is reported in this paper for the sake of completeness. What will be proven is that both properties lead to doubtless computational advantages.

2. Linear Orbit Theory

The relative motion of a spacecraft with respect to a virtual or real target vehicle, placed in a neighboring reference orbit, can be described using linear orbit theory [9], which provides the equations of relative motion, referred to as the Euler-Hill-Clohessy-Wiltshire (EHCW) equations.

The EHCW equations represent the linear expansion of the dynamics equations that govern position and velocity of the pursuing spacecraft () relative to the target vehicle (), which is placed in a reference circular orbit of radius . If is the gravitational parameter of the attracting body, the mean motion of (in this case equal to its constant angular rate) is . Let and represent the inertial position and velocity vectors (with or ) and the relative position of . The EHCW equations are written in terms of the projections of along the right-hand reference frame , where and are, respectively, the unit vectors associated with the target position vector and the specific angular momentum : . The EHCW equations are [9]where can be regarded as a perturbing or thrust acceleration. Letting , (1)–(3) can be rewritten in the form of state equations as whereVector equation (4) is equivalent to six linear differential equations with constant coefficients and admits a general solution composed of two termswhere is the state transition matrix, which obeys the following equation:( is the 6 by 6 identity matrix). Letting , an analytical expression exists for , where and .

The first term of (6) is the general solution of the homogeneous system associated with (6); that is, it is the solution of the EHCW equations in the absence of . If the initial time is set to 0, it is relatively straightforward to show that—omitting the velocity components for the sake of brevity—the homogeneous term (denoted by the subscript “”) can be conveniently rewritten aswhere the constants and the angles and depend on the initial conditions and have the following expressions:Equations (9)–(12) allow identifying some special solutions to the homogeneous EHCW equations of motion. The following solution is of special interest. (a) Relative elliptic motion, if : this special solution corresponds to the absence (in (10)) of the second term, responsible for the drift of with respect to . If also then the ellipse lies in the -plane (the spacecraft motion is two-dimensional).From inspecting (9)–(11) it is apparent that the dynamics of the (out-of-plane) components of position and velocity (i.e., and ż) are completely uncoupled from the dynamics of the remaining (in-plane) components. As an immediate consequence, the analysis and optimization of planar cases (involving only the components and ) can be performed by simply setting and ż = 0 and neglecting the respective equations, which become unnecessary.

3. Optimal Relative Orbit Trajectories

In general, in relative orbit maneuvering two types of optimal paths are sought: (a) trajectories that minimize fuel expenditure in a given time and (b) paths that minimize the time of flight. They are related to specific analytical conditions for optimality, which are the subject of this section.

3.1. Optimal Thrust Programming for Minimum-Fuel Relative Trajectories

This subsection is focused on the problem of minimizing the fuel consumption needed for performing orbital rendezvous or transfers in a specified time . The thrust acceleration (with magnitude denoted by ) represents the control authority over the dynamics of . The pursuing spacecraft employs a time-varying thrust (with magnitude ), and therefore is given by , with constrained to . Moreover, letting (where is the initial spacecraft mass), the following additional (scalar) state equation is introduced:The component is included in the state vector . In (13) is the effective exhaust velocity of the propulsive system, which is assumed constant. The thrust direction is defined by two angles: (out-of-plane angle), constrained to , and (in-plane angle), constrained to . As a result, the thrust acceleration components are

The formulation of the optimal control problem is thus the following: determine the optimal thrust magnitude and direction that minimize fuel consumption; that is,where represents the control vector. Constraints for the problem at hand are equations of motion (4) and (13) and boundary conditions of the formwhere the subscripts “” and “” denote the initial and final values of the respective variable.

A Hamiltonian and a function of the terminal conditions are introduced, with the intent of deriving the necessary conditions for optimality [34]where and are the adjoint variables conjugate to the dynamics equations and to the boundary conditions, respectively. First, the Pontryagin minimum principle [34] assumes the form(the symbol “” denotes the optimal time history of the respective control variable). As the term is nonnegative, for (18) can be rewritten asand yields the following relations: Equations (20) are meaningless during coast intervals. After inserting (20), (18) becomesThe term in parentheses is usually referred to as the switching function , because it determines the optimal sequence of thrust and coast intervals. In fact, under the assumption that does not vanish over finite time intervals (i.e., excluding singular-control arcs), (21) yields the optimal thrust magnitudeThis means that minimum-fuel trajectories are composed of maximum-thrust intervals and coast arcs. Their optimal sequence is established by , according to (22). Moreover, the adjoint (or costate) equations for are [34]The six components of can be expressed in closed form, after deriving the state transition matrix for the linear system (23) (first equation),The matrix fulfills a differential equation formally identical to (7), with in place of . Letting , assumes the following form:It is worth remarking that the adjoint vector has the same analytical form along thrust intervals and coast arcs. The remaining necessary conditions regard the initial and final times, and . Since these conditions are problem-dependent, only their general form is reported,where the subscripts “” and “” refer, respectively, to and . As is unspecified, for all of the cases that are being considered.

3.2. Optimal Thrust Programming for Minimum-Time Relative Trajectories

In this subsection the problem of determining minimum-time rendezvous and transfer trajectories in the Euler-Hill frame is considered. The thrust acceleration represents again the control authority over the dynamics of and has components given by (5). The state vector includes the same variables introduced in the preceding subsection; the initial time is specified (and set to 0), whereas the final time is unspecified and is to be minimized.

The formulation of the optimal control problem is thus the following: determine the optimal thrust magnitude and direction that minimize the time of flight; that is,where represents the control vector. Constraints for the problem at hand are equations of motion (4) and (13) and boundary conditions of form (16).

A Hamiltonian and a function of the terminal conditions are introduced, with the purpose of deriving the necessary conditions for optimality,The Pontryagin minimum principle assumes the formAs the term is nonnegative, (29) yields again relations (20). After inserting (20), (29) becomesDue to (28), the adjoint equations are the same as those reported in (23), with related boundary conditions (26), which yield again , because is unspecified. The latter condition, in conjunction with the adjoint equation for , implies that . As a result, the term in parentheses in (30) is always negative, and this leads to determining the optimal thrust magnitude that minimizes the Hamiltonian; that is,Moreover, as is unspecified, the transversality condition holds,

4. Symmetry Properties

This section addresses a pair of symmetry properties that hold in the presence of certain boundary conditions:(a)symmetry of two-dimensional and three-dimensional minimum-fuel and minimum-time paths,(b)symmetry of the optimal thrust direction and sequence for two-dimensional minimum-fuel trajectories.

4.1. First Symmetry Property

The first interesting symmetry property regards minimum-fuel and minimum-time trajectories in the Euler-Hill frame. In both cases, the adjoint variables have the same (analytical) expressions and determine the optimal thrust direction (cf. (20)), as well as the optimal thrust sequence of minimum-fuel paths (cf. (22)).

Proposition 1. Let be the optimal solution of problem A, which is in determining either the minimum-fuel (I) or the minimum-time (II) relative trajectory, while holding equations of motion (4) and (13) and the boundary conditions The set is the optimal solution to problem B, formulated either as (I) or as (II), with boundary conditionsprovided that

Proof. By definition, the solution meets the necessary conditions for optimality, that is, (20) and (23)–(26), in conjunction with (22) for objective (I) and (31) and (32) for objective (II). Omitting the superscript “” for the elements of solution henceforward, (26) yield Now, if for problem the adjoint vectors , , and are assumed to be given bythen adjoint equations (23) and related boundary conditions (26) are fulfilled by , , and . In fact, the adjoint equations, written for , reduce to those holding for . In addition, due to (36), the boundary conditions for imply satisfaction of the corresponding boundary conditions for ,With regard to the control variables, (22) (for objective (I)) and (31) (for objective (II)) in conjunction with (20) and (38) yieldDue to the previous relations on and , if , then the adjoint equation for and the related boundary condition () are both satisfied. Insertion of the control variables in the state equations (1) leads to their fulfillment, because these equations reduce to the corresponding equations that hold for solution . In conclusion, since all the necessary conditions for optimality are met, is the optimal solution for problem , formulated either as minimum-fuel (I) or as minimum-time (II) problem.

4.2. Examples of Symmetric Minimum-Time Relative Trajectories

Two examples that fulfill the conditions stated in Proposition 1 are reported in this subsection, for the purpose of illustrating the previously described symmetry property. The reference circular orbit around the Earth has an altitude of 400 km. The thrust magnitude is set to 0.01 m/sec2, whereas the effective exhaust velocity is 25 km/sec.(a)Transfer between two coplanar relative elliptic orbits (Figure 1): the problem consists in determining (i) the departure point from the periodic orbit (with semimajor axis of 20 km) and (ii) the thrust direction that minimize the time needed to inject the pursuing spacecraft in a different coplanar relative elliptic orbit (with semimajor axis of 2 km).(b)Three-dimensional rendezvous from a relative elliptic orbit (Figure 2): the problem consists in determining (i) the departure point from the periodic orbit and (ii) the thrust direction that minimizes the time needed to rendezvous with the target, denoted by in Figure 2(b). The initial relative elliptic orbit is associated with the position coordinates (4)–(6), with ,  km,  km, and . The condition corresponds to relative elliptic orbits (as previously remarked), implies symmetry of the relative orbit with respect to the target location, whereas the values and and the angular relation define, respectively, the ellipse size and spatial orientation.Further details pertaining to computational aspects are reported in [32]. It is apparent that the utility of Proposition 1 is in determining immediately the second optimal path and the related optimal control time history, once the first optimal trajectory has been found. This may be relevant in the presence of specific practical constraints (e.g., if the transfer must be performed at the earliest occurrence, one can choose between the two options represented by trajectory 1 and path 2, with no difference in terms of time of flight and propellant expenditure).

4.3. Second Symmetry Property

A further symmetry property concerns minimum-fuel two-dimensional trajectories (i.e., , , , and   ). As a preliminary step, from inspection of (20) and (22) it is apparent that the thrust direction and magnitude depend only on the adjoint variables . It is relatively straightforward to rewrite their analytical expressions aswhere , , and depend on , , , and (the subscript “0” corresponds to ). Under the assumption that , the geometrical locus associated with is either an ellipse (if ) or a cycloid (if ). Figure 3 illustrates a cycloid in an orbital period and shows the optimal control direction, determined using (20).

If the effective exhaust velocity is very high, then the mass depletion rate is modest, and the following approximation can be adopted: . This implies that   . As a result, (13) becomes unnecessary, as well as the respective adjoint equation for (cf. (23)). Under this assumption the switching function assumes the form Hence, in the -plane the condition for thrust phases to occur () corresponds to the cycloid (or ellipse) arcs that lie outside the unit circle. Intersections of the cycloid with the unit circle correspond to the switching times between thrust intervals and coast arcs. In an orbital period, their maximum number equals 6 (for a cycloid) or 4 (for an ellipse). This means that at most four powered phases, separated by three coast arcs, can occur. This property has been well known for impulsive relative paths since the 1960s [7, 8] and holds also for finite-thrust trajectories, due to the analytic form of the adjoint variables and , as also remarked by Carter [24].

Proposition 2. With regard to two-dimensional minimum-fuel paths, if the maximal thrust acceleration is constant, that is,   , and ifthen (i) the optimal thrust sequence is symmetric with respect to the intermediate time and (ii) the optimal thrust angle fulfills the following relation:

Proof. As a first step, the time variable is defined as . Then, the state components and are easily obtained from (6):It is worth noticing that in this subsection the subscript “0” corresponds to the intermediate time (i.e., ). Equations (45) are evaluated at the initial time and at the final time ; then, they are inserted in (43), to yieldAfter combining (46), the following relation is obtained:Due to the geometrical properties of the -locus, (47) is satisfied only if in the interval (i.e., ) the portion of cycloid (or ellipse, in the degenerate case ) is symmetrical with respect to the -axis. This implies also that the following two properties hold: (a) the function is odd and (b) the function is even. Properties (a) and (b) are equivalent to (44) (statement (ii)). In addition, symmetry of the -locus implies also the symmetry of the optimal thrust sequence (statement (i)).

Due to the analytical form of (appearing in (25)), it is also immediate to recognize that (a) is even and (b) is odd. As a result,

As an immediate consequence of the previously mentioned symmetry property, the following integral equalities can be easily derived, by considering the properties of odd and even functions of time:These relations can be profitably employed to obtain the values of and at , that is, and . With regard to , the use of the related analytic expression deriving from (6) (evaluated at and ), together with (50) and (51), yields Then, the corresponding expression for , together with (49)–(51), leads to the following relationship:Equations (52) and (53) prove that the values of and at the intermediate time depend only on the respective initial and final values, as well as on the time of flight .

Lastly, due to the analytic form of the adjoint variables, (48) yield useful relations among the initial values of the costate variables. In fact, letting , and can be rewritten asLetting and , after inserting these expressions in (48), one obtainsEquations (55) can be employed to derive two initial costate variables (e.g., and ) as functions of the remaining ones.

In conclusion, if the state components are specified at both the final and the initial time and boundary conditions (43) hold, then the numerical solution of the relative trajectory optimization problem may be sought(a)in the time interval , instead of ,(b)by using two initial values of the adjoint variables as unknown parameters, instead of four unknown parameters ,(c)by assuming (52) and (53) as the only boundary conditions, without the need for enforcing the conditions on and .Points (b) and (c) imply that the parametric problem is reduced to half of its original size, whereas point (a) means that each trajectory is to be computed up to the intermediate time during the numerical solution process. These circumstances have very positive effects on computational efficiency. In fact, if the indirect heuristic method is used, halving the number of individuals and the propagation time (for each trajectory) implies reducing the computational runtime by a factor 4.

4.4. Example of Symmetric Minimum-Fuel Relative Trajectory

The optimal relative trajectory presented in this subsection is very similar to the four-thrust-arc transfer treated by Carter [24]. The reference orbit is a low Earth circular orbit at an altitude of 435 km, while m/sec2 and  min. The boundary conditions areFor this case special symmetry conditions (43) and (44) are satisfied and the optimal trajectory is composed of four powered phases. Figures 4 and 5 illustrate the optimal trajectory and the related velocity components, thrust pointing angle, and switching function. Further details pertaining to computational aspects are reported in [33].

5. Conclusion

This work is focused on a pair of symmetry properties that hold for optimal relative orbit trajectories. The existence of symmetry properties is proven in the presence of special boundary conditions for the problems of interest. First, both minimum-time and minimum-fuel paths are proven to enjoy a property that leads to defining the optimal trajectory (and the respective optimal control law) that is symmetrical with respect to a given, already determined optimal path. Minimum-fuel relative trajectories are composed of multiple powered phases (at most 4 in an orbital period), divided by coast arcs. The optimal thrust sequence, duration, and direction depend on the time evolution of the switching function and satisfy a second interesting symmetry property, provided that a pair of elementary boundary conditions holds. This allows halving the number of boundary conditions and unknown parameters with respect to the original formulation of the problem. Some examples illustrate the properties at hand, pointing out their doubtless interest and usefulness for the purpose of simplifying the determination of optimal relative orbit maneuvers.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.