Recent Advances on Methods and Applications of Nonlinear Differential Equations
View this Special IssueResearch Article  Open Access
Solution for Nonlinear ThreeDimensional Intercept Problem with Minimum Energy
Abstract
Classical orbit intercept applications are commonly formulated and solved as Lamberttype problems, where the timeofflight (TOF) is prescribed. For general threedimensional intercept problems, selecting a meaningful TOF is often a difficult and an iterative process. This work overcomes this limitation of classical Lambert’s problem by reformulating the intercept problem in terms of a minimumenergy application, which then generates both the desired initial interceptor velocity and the TOF for the minimumenergy transfer. The optimization problem is formulated by using the classical Lagrangian and coefficients, which map initial position and velocity vectors to future times, and a universal time variable . A NewtonRaphson iteration algorithm is introduced for iteratively solving the problem. A generalized problem formulation is introduced for minimizing the TOF as part of the optimization problem. Several examples are presented, and the results are compared with the Hohmann transfer solution approaches. The resulting minimumenergy intercept solution algorithm is expected to be broadly useful as a starting iterative for applications spanning: targeting, rendezvous, interplanetary trajectory design, and so on.
1. Introduction
A fundamental problem of astrodynamics is concerned with computing intercept trajectories or interplanetary mission orbit for objects in space [1, 2]. These calculations are often performed assuming a predetermined timeofflight (TOF). This is the wellknown Lambert’s problem [3–6]. Applications of Lambert’s problem are common in interplanetary trajectory design, spacecraft intercept, rendezvous, ballistic missile targeting, and so on. These problems are formulated by specifying the initial position vectors of an interceptor and target satellite. When the TOF is specified, the initial velocity vector for the interceptor is an unknown implicit function of the local gravity field and can only be recovered by a successive approximation strategy. Other authors have considered alternative problem formulations for a specified TOF that have included minimumfuel and multipleimpulse strategies [7, 8]. This work overcomes the limitations of these approaches by introducing a minimumenergy approach that simultaneously generates both the TOF and the initial velocity vector for the interceptor.
For the special case of coplanar orbits, the Hohmann transfer algorithm generates a twoimpulse minimumenergy orbit transfer by using tangential burns [3–5]. This technique provides a reference orbit transfer for various space applications. For direct applications of the Hohmann transfer to interplanetary orbit transfer, the position vectors of the target planet and initial departure planet are specified assuming a prescribed TOF. When the spacecraft is far from the initial position, one must be alert to the possibility that a multiorbit maneuver may be required.
Clearly, the TOF is a critical parameter for various applications. Once a TOF is determined, the rest of the procedure is solved readily by the solution of Lambert’s problem. This work addresses the problem that there are no adequate methods available for determining a TOF, especially, in general threedimensional (3D) cases. The problem of finding an optimal TOF only becomes well defined when one specifies a minimization criteria. To this end, a minimumenergy version of classical interceptor problem is formulated for recovering the TOF for a 3D orbit transfer. The results of this calculation are useful as a reference value for interplanetary trajectory design, spacecraft intercept, rendezvous, ballistic missile targeting, and so on. Of course, one can also bound the range of achievable transfer trajectory times by solving for the minimum TOF consistent with the maximum energy that can be generated. Yielding a mission design space that spans the range of TOF consists in the range.
The design goal for the optimization problem is to simultaneously recover the required initial interceptor velocity and the TOF for the intercept. The mathematical advantage of this approach is that the problem has a unique optimal solution, rather than the family of solutions that characterize the classical Lambert’s problem. Mathematically, the problem is defined by a constrained optimization algorithm. Particular care is exercised in formulating the problem for handling the nearparabolic orbits that arise in intercept applications. Analytically, this is handled in a comprehensive way by introducing a universal variable that permits a single TOF equation to be developed that is valid for all conic orbits.
This work is organized in three sections. First, Kepler’s equation is used to define the TOF equation. This is followed by a description of the universal variable used for the problem formulation. For completeness, Lambert’s problem is briefly described. Second, the minimumenergy problem for the intercept problem is introduced and solved. Third, simulation results are presented which compare the TOF obtained for an interplanetary trajectory design with a trajectory developed using the Hohmann transfer methodology and interceptor design solution approaches.
2. Mathematical Review
A fundamental approach for determining the TOF for spacecraft starts with Kepler’s equation that is given by whereis the mean anomaly,anddenote the eccentric anomaly and the eccentricity, respectively,is the time of periapsis passage,is the TOF,is the mean motion defined as,denotes the gravitational constant, andis the semimajor axis of orbit.
Asthe solution for Kepler’s equation becomes more difficult to obtain. This problem is overcome by introducing the universal variable given by [3] whereis the position of spacecraft. As shown in [3], by introducing the universal variable defined by (2), one can express Kepler’s equation and the radial spacecraft coordinate in the following form: whereis assumed to be zero without loss of generality andandare the initial position and velocity vectors of spacecraft, respectively. These necessary conditions describe the position and velocity of an orbiting object as a function of time. If the value of the universal variable from (3) is known, the position of the spacecraft at that time is evaluated. Even though (3) is transcendental in, a Newton’s iteration technique is used to successfully solve forwhen the TOF,, is given.
Assuming that there are no external forces, then the four vectors,,, andare assumed to be governed by Keplerian motion. To computeandin terms of,, and, the position and velocity vectors of spacecraft at timeare described as [9] where,,, andare scalar timedependent constants, which are subject to the following constraint: where
The energy minimum form of Lambert’s problem is solved by introducing the classical Lagrangian coefficients and universal variable in the problem necessary conditions.
3. TimeofFlight for MinimumEnergy Orbit Transfer
The major objective in this paper is to compute (i) the TOF and (ii) the initial velocity for an interceptor object for two arbitrary given position vectors so that the transfer orbit energy is a minimum. This problem differs from the classical Lambert’s problem, which fixes a TOF and only recovers a solution for the initial velocity for the interceptor object, given initial and final position vectors of spacecraft.
The optimization problem is formulated by assuming that a target and an interceptor exist in arbitrary orbits, respectively. The problem geometry is illustrated in Figure 1, whereanddenote the initial position and velocity vectors of the target, respectively, andandrepresent the initial position and velocity vectors of interceptor, respectively.
The unknowns for the problem are the TOF and initial velocity correction for the interceptor. The goal of the trajectory optimization is to reduce the displacement position vector locating the interceptor relative to the target to zero values at the TOF, while minimizing the orbit energy of the interceptor. The problem is formulated as a nonlinear optimization problem.
4. Constrained Optimization Problem
For given,, and, findandby minimizing the performance index defined as the interceptor’s orbit energy,, defined as subject to whereand(9) denote the universal variables for the target and the interceptor, respectively, andis the TOF to be determined.
The displacement vectors for the target and interceptor are expressed usingandas follows:
As a constraint vector,, (3) forandis rewritten as Note that the augmented variables to be obtained are,,, and.
5. Optimal Necessary Conditions
Since the second term of the energy is constant, it does not affect the performance index so that the index is redefined, without loss of generality, as [10]
The Hamiltonian is formed by appending the constraints of (9) with Lagrange multipliers as follows: whereand. To minimize the performance index with respect to the augmented variables, the necessary conditions provide the following [11]: whererepresents a row vector,is a matrix (refer to the Appendix for detail derivation), and
The necessary conditions of (14)–(17) are simplified by the following manipulations. First, from (15), the Lagrange multiplieris obtained as
Second, substituting (19) into (17) yields which can be solved for, leading to
Third, by using (14), the Lagrange multiplierfollows as
Collecting the Lagrange multiplier solutions from (19) and (22), introducing the results into (16), one obtains
This equation is further simplified by recalling the terminal constraint, leading to
Substituting (21) into (3) yields the final necessary condition required for finding the TOF for the intercept problem: where the new matrix is defined for simplicity as
Satisfaction of the new equation implies that the interceptor can meet the target with minimum energy in a time provided by the computed TOF, not by a predetermined time.
6. Summary
The approach for obtaining the nonlinear 3D intercept problem by using the classical Lagrangianandcoefficients is summarized. Computing the TOF and the interceptor’s initial velocity is the goal to meet the requirement that the final distance between the two spacecrafts becomes zero. Firstly, one can simply obtainandusing the following: where the orbital energy is given by
Then, one can find the universal variables, the initial velocity for the interceptor, and the TOF using the following equations: where the semimajor axisof the interceptor can be iteratively computed with estimated. The NewtonRaphson iteration algorithm is applied to solve the previous equations. Next, one can compute all of theandexpressions using (7). Then, (5) is applied to obtain the final position and velocity vectors.
There are many feasible performance indices to specify a TOF. For example, consider the candidate performance index whereis nonnegative weight. By adding the time as one part of the performance index, the TOF is expected to be shortened with respect to the variation of. In a similar manner with the minimumenergy procedure in the previous section, the optimization solution to this problem is readily determined. The partial derivative of the Hamiltonianwith respect tois given by
Finally, a costeffective equation weighted to the time is obtained as
Numerical convergences based on different methods and their overall computational cost depend on the chosen parameterization, the initial guess, and the numerical technique used for solving the resulting equation. It is known that singularities exist when solving Lambert’s problem that prevent some algorithms from converging for particular cases or make convergence extremely slow. For example, Lambert’s method fails when the transfer angle is 180 degrees [9]. Therefore, the features of the suggested method must be analyzed. However, this is out of the scope of this paper, which is focused on approaches to determine the TOF and initial velocity of the interceptor.
7. Application Examples
The specification of a TOF for an intercept problem is generally not unique, and a family of solutions are possible when the initial trust level is variable. As a result, the process of determining a useful TOF requires experimentation and iteration. The minimumenergy optimization approach of this work finds a unique value for the TOF. The solution for the intercept problem simultaneously determines the initial interceptor velocity vector and TOF. Numerical examples are presented that compare and contrast the classical Hohmann transfer with the proposed method.
Let us briefly review the Hohmann transfer and compare the minimumenergy problem with it. The geometry of the Hohmann transfer is illustrated in Figure 2. The distances of the departure and arrival orbits are denoted asand, respectively.
The semimajor axis for the elliptic orbit and the energy are given by [3] and, then, the departure velocity of the transfer orbit is readily obtained as
Since the velocity of the departure orbit is given by the velocity change for the Hohmann transfer is calculated as and the TOF of the Hohmann transfer is written as whereis obtained from (33). Two circular orbits are assumed with the radii ofkm andkm, respectively. Then, the velocity of the departure circular orbit iskm/s, and the remaining parameters for the Hohmann transfer are obtained assec,km/s, andkm/s. To navigate to the final position of the arrival orbit by the Hohmann transfer, the initial position and velocity vectors are assumed to be given by
Four cases including the Hohmann case are analyzed. Initial velocities, velocity changes, and TOF obtained by the solution of the proposed minimumenergy problem are arranged in Table 1. Also, the initial positions of target spacecraft and their resultant transfer trajectories are displayed in Figure 3.

As shown in Table 1, the result for the Hohmann case is nearly identical to the output from the classical approach in (35)–(38) with a small numerical error. It proves that the proposed approach provides optimal solutions we are looking for. Moreover, it is obvious that it gives the minimum velocity change, which is tangential with the trajectory, compared with the other cases. If the target is positioned at case 1, 2, or 3, relative to the interceptor’s initial position and it is required to start the orbit transfer mission at this time, it would be a great advantage to have a reference minimumenergy trajectory to accomplish the mission. Fortunately, the results in Table 1 can be utilized, since they represent the minimum velocity in each case. This means that there are no more efficient trajectories in these cases than the transfer orbit listed in Table 1. When the target is positioned forward compared to the Hohmann transfer, the phase angle, sometimes called flightpathangle, at departure should be negative to meet the optimal trajectory requirement. When the target, on the contrary, is positioned backward, the flightpathangle should be positive.
Even if a circular orbit is selected for the comprehensive analysis by comparing with the Hohmann transfer, the application of the proposed approach is not limited. Therefore, an illustrative example in Figure 4 is conducted to demonstrate the performance of space maneuver of the interceptor. There are two arbitrary elliptic orbits, which are not coplanar. The initial positions of the target and interceptor orbit are depicted in Figure 4. By solving the nonlinear 3D intercept problem through minimizing energy, the TOF and initial velocity is computed. Then, applying minimumenergy velocity obtained for the interceptor confirms that the final distance between the two orbits is zero at the computed TOF. Note that the problems formulated by the universal variable andandexpressions in this paper are solved in 3D space for supporting the design of arbitrary intercept problems with minimum energy.
The intercept problem is easily generalized by introducing a time weighting factor in the definition of the optimal control problem, which allows a systematic exploration of the intercept design space as a function of the allowed transfer energy. In this example, the outer orbit is selected as a target orbit to be captured as illustrated in Figure 5. The initial position of the interceptor is at the inner orbit. The TOF obtained from the optimization problem is shortened whenincreases, and the results are illustrated in Figure 5. As expected, the longest TOF is obtained whenis zero. Ifapproaches one in this simulation case, the interceptor can hit the target in a very short time. It means that shortening the TOF is the optimal solution to minimize the chosen cost function.
8. Conclusions
A general algorithm is presented for generalizing the classical Lambert’s transfer problem, where the determination of timeofflight (TOF) for a spacecraft intercept, in arbitrary threedimensional orbit, is addressed. A constrained optimization technique is introduced and iteratively solved to find both the TOF and the initial intercept velocity vector. The proposed algorithm provides a benchmark minimumenergy solution that provides an optimal reference trajectory. A significant advantage of this approach is that the TOF is unique when the energy is minimized. This implies that the interceptor with lower energy than the evaluated minimum energy cannot meet the target. Numerical results are presented, and they compare the intercept solutions with those obtained using the classical Hohmann transfer technique. The proposed algorithm is expected to be broadly useful for all classes of intercept problem that have a Lambertlike character.
Appendix
The Partial Derivatives
The partial derivativesandwith respect toand, respectively, are given by
Applying (4) and using (7), the partial derivatives with respect to the universal variable are written as
The orbit energy has the following relationship:
Sinceis a constant in this case, differentiating both sides with respect toyields
The partial derivative ofwith respect tois readily written as
The row vector, partial derivative ofwith respect to, is given by where
Next, the partial derivative ofwith respect tofollows as where
The matrix,, in (26), consists of the combination of two vectors and one scaled identity matrix. Ifis not zero, the matrix would have a full rank of three. The position and velocity vectors in orbit are in general not parallel. Since the condition thatis zero means thatandare parallel, this is impossible in orbit. It could guarantee the existence of the inverse of the matrix.
Acknowledgment
This study was supported by research fund from Chosun University, 2013.
References
 M. Hawkins, Y. Guo, and B. Wie, “Spacecraft guidance algorithms for asteroid intercept and rendezvous missions,” International Journal of Aeronautical and Space Sciences, vol. 13, no. 2, pp. 154–169, 2012. View at: Google Scholar
 T. S. No, J. M. Lee, G. E. Jeon, D. Lee, and G. Kim, “A study on earthmoon transfer orbti design,” International Journal of Aeronautical and Space Sciences, vol. 13, no. 1, pp. 106–116, 2012. View at: Google Scholar
 R. R. Bate, D. D. Mueller, and J. E. White, Fundamentals of Astrodynamics, Dover Books on Aeronautical Engineering, ch 4–7, Dover, Mineola, NY, USA, 1st edition, 1971.
 D. A. Vallado and W. D. McClain, Fundamentals of Astrodynamics and Applications, Microcosm Press/Springer, El Segundo, Calif, USA, 2007. View at: MathSciNet
 V. A. Chobotov, Orbital Mechanics, AIAA Education, ch 4, AIAA, Reston, Va, USA, 3rd edition, 2002.
 A. Prado and R. A. Broucke, “The minimum deltaV Lambert’s problem,” Sociedade Brasileira de Automática, vol. 7, no. 2, pp. 84–90, 1996. View at: Google Scholar
 J. E. Prussing and J.H. Chiu, “Optimal multipleimpulse timefixed rendezvous between circular orbits,” Journal of Guidance, Control, and Dynamics, vol. 9, no. 1, pp. 17–22, 1986. View at: Google Scholar
 H. Shen and P. Tsiotras, “Optimal twoimpulse rendezvous using multiplerevolution lambert solutions,” Journal of Guidance, Control, and Dynamics, vol. 26, no. 1, pp. 50–62, 2003. View at: Google Scholar
 R. H. Battin and R. M. Vaughan, “An elegant Lambert algorithm,” Journal of Guidance, Control, and Dynamics, vol. 7, no. 6, pp. 662–670, 1984. View at: Google Scholar
 H. Leeghim and B. A. Jaroux, “Energyoptimal solution to the lambert problem,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 3, pp. 1008–1009, 2010. View at: Publisher Site  Google Scholar
 A. E. Bryson, Jr. and Y. C. Ho, Applied Optimal Control: Optimization, Estimation and Control, ch 1, Taylor & Francis, 1975. View at: MathSciNet
Copyright
Copyright © 2013 Henzeh Leeghim et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.