Research Article  Open Access
Asteroid Rendezvous Mission Design Using Multiobjective Particle Swarm Optimization
Abstract
A new preliminary trajectory design method for asteroid rendezvous mission using multiobjective optimization techniques is proposed. This method can overcome the disadvantages of the widely employed PorkChop method. The multiobjective integrated launch window and multiimpulse transfer trajectory design model is formulated, which employes minimumfuel cost and minimumtime transfer as two objective functions. The multiobjective particle swarm optimization (MOPSO) is employed to locate the Pareto solution. The optimization results of two different asteroid mission designs show that the proposed approach can effectively and efficiently demonstrate the relations among the mission characteristic parameters such as launch time, transfer time, propellant cost, and number of maneuvers, which will provide very useful reference for practical asteroid mission design. Compared with the PCP method, the proposed approach is demonstrated to be able to provide much more easily used results, obtain better propellantoptimal solutions, and have much better efficiency. The MOPSO shows a very competitive performance with respect to the NSGAII and the SPEAII; besides a proposed boundary constraint optimization strategy is testified to be able to improve its performance.
1. Introduction
The optimization of interplanetary trajectories to an asteroid continues to arouse a great deal of interest [1–4]. Although the lowthrust propulsion is employed in asteroid rendezvous missions, an impulsive trajectory is always assumed for preliminary mission design and optimization, in which the access feasibility is evaluated, the launch window is determined, the gravityassist maneuvers scheme is planned, and so forth. In a preliminary asteroid rendezvous trajectory design, the PorkChop method is widely used. This method employs the twoimpulse algorithms including the optimal twoimpulse noncoplanar transfer [1, 2], the classical Hohmann transfer [3], and the twoimpulse Lambert algorithm [4] to design the transfer trajectory. The total corresponding to different departure time and arrival time is then obtained and the contours of minimum total are plotted to assist the designer to find the best launch window and transfer trajectory.
This type of method is very intuitionistic and easily executed. However, only the twoimpulse trajectory is investigated in this method, and as demonstrated by Lawden’s theory [5, 6], the twoimpulse trajectory is not the propellantoptimal solution under most conditions. Besides, this method is in essence of an exhaustive searching method, its computation cost increases exponentially as the search space increases, and much human intervention is required as no orderliness exists in most of the contours.
In this paper, a new preliminary trajectory design method for asteroid rendezvous mission using multiobjective optimization techniques is proposed. This method can overcome the disadvantages of the PorkChop method. The multiobjective integrated launch window and multiimpulse transfer trajectory design model is formulated, which employs minimumfuel cost and minimumtime transfer as two objective functions. In this model, the Earth departure date, hyperbolic velocity, and the interplanetary transfer impulses are all chosen optimization design variables. The multiobjective particle swarm algorithm is employed to locate the Pareto solution, by which the relationships characteristics among the overall mission parameters can be effectively revealed.
The particle swarm optimization (PSO) algorithm firstly introduced by Kennedy and Eberhart [7] in 1995 has been successfully applied in many fields of research. However, it is only recently applied to aerospace trajectories optimization [8–11]. The singleobjective PSO algorithm has been testified as one successful spacecraft trajectory optimizer in designing impulsive interplanetary trajectories [8], lowthrust trajectories for asteroid exploration [9], Lyapunov and Halo (periodic) orbits [10], and multipleburn rendezvous trajectories [11].
The implementation of the PSO algorithm adopts a population of particles, whose behavior is affected by either the best local (i.e., within a certain neighborhood) or the best global individual. The relative simplicity of PSO and the fact that it is a populationbased technique have made it a natural candidate to be extended for multiobjective optimization. In a survey paper in 2006 on multiobjective particle swarm optimization (MOPSO) [12], it was reported that there were currently over twentyfive different proposals of MOPSO reported in the specialized literature. The studies on MOPSO remain a very active area of research, and the MOPSO has been successfully applied to many practical multiobjective optimization problems, recently applied to robotics [13], industrial management [14], and chemical engineering [15]. It also has been applied in the domain of aerospace including airfoil shape optimization, complex physics/shape optimization, and multidisplinary design optimization [16]. However, it was seldom applied in solving multiobjective spacecraft trajectory optimization problems. The recent studies show that the singleobjective PSO algorithm is one effective, reliable, and accurate spacecraft trajectory optimizer [8–11]. In this study, we will show that the MOPSO could be one successful optimizer for multiobjective spacecraft trajectory optimization problems.
In summary, the main contribution of this paper is twofold. (1) A novel asteroid rendezvous mission design method using the multiobjective techniques is proposed. Compared with the current widely employed PorkChop method, the proposed approach is demonstrated to be able to provide much more easily used results, obtain better propellantoptimal solutions, and have much better efficiency. (2) As far as we know, it is the first time to apply the MOPSO to spacecraft trajectory optimization. The MOPSO proves to be quite effective in finding the Paretooptimal solutions to asteroid rendezvous multiobjective optimization problems. The MOPSO shows a very competitive performance with respect to two highly competitive multiobjective evolutionary algorithms: the nondominated sorting genetic algorithmII (NSGAII) [17] and the strength Pareto evolutionary algorithmII (SPEAII) [18].
2. Optimization Problem
2.1. Asteroid Rendezvous Design Problem
The interplanetary transfer trajectory is always divided into three different segments, that is, planet departure segment, heliocentric transfer segment, and capture segment by using the concept of influencing sphere. As the influencing sphere of the planet is much smaller than that of the sun, the flight path and flight time of the departure and capture segments are much small compared with that of the heliocentric transfer. Therefore, for a preliminary trajectory design of an asteroid exploration mission, the design emphasis is firstly focused on the heliocentric transfer, while the planet departure is assumed as instantaneous process with an impulsive maneuver and the capture segment is omitted.
This paper studies the asteroid rendezvous problem departure from the Earth. represents the spacecraft launch date which is to be designed, and , , and define the departure hyperbolic velocity according to the formulas
Let , be the heliocentric position and velocity of the Earth at ; then the initial position and velocity vector of the spacecraft, and , entering the heliocentric transfer trajectory are defined as
The heliocentric transfer trajectory is modeled by a twobody dynamic model with the following governing equations: where , , and are position components of spacecraft, , , and are velocity components of spacecraft, , , and are acceleration component of spacecraft, and is the gravity parameter of the sun.
The thrust acceleration can be approximated as impulses where is the time where an impulse is applied, , and they are all the variables to be designed. The final state conditions are defined where is the rendezvous time of spacecraft with the target asteroid and and are the position and velocity vector of the target asteroid at the rendezvous time.
2.2. FeasibleSolution Iteration Model
In order to avoid dealing with the equality constraints described as (5), the Lambert algorithm is employed to establish the infeasible iteration optimization model. The chosen independent variables, that is, the optimization variables, are impulse times and the first impulses
The last two impulses are determined by solving the Lambert problem constrained by (5). In this feasible iteration model, each evaluation of the objective function produces a feasible solution that satisfies implicitly the rendezvous conditions. Detail on this multiimpulse rendezvous optimization model using the Lambert algorithm can be found in [19, 20].
2.3. Multiobjective Optimization Model
The total velocity characteristic is chosen as the first objective function
The heliocentric transfer time is chosen as the second objective function
The constraint on the time of impulse is considered. The general constraint on is
As the transfer time is one of the objective functions, is chosen as an optimization variable. The total optimization variables include four parts: , , , , , , and . In order to improve optimization performance, the variablescaling method is imposed on . Let . Therefore, the final optimization variables are
3. Optimization Algorithms
3.1. Introduction of MOEA
A general multiobjective optimization problem is to find the design variables that optimize a vector objective function over the feasible design space. The objective functions are the quantities that the designer wishes to minimize, maximize, or attain a certain value. The problem formulation in standard form for a minimization is given here, which is similar for the other cases: where
For the multiobjective asteroid rendezvous design problem, the two objective functions are described by (7) and (8), the constraints are described by (9), and the optimization variables are described as (10).
The classical optimization method for a multiobjective optimization problem is the weighting method. In recent years, the multiobjective evolutionary algorithms (MOEA) have been greatly investigated in the domain of multiobjective optimization. There are many variants of MOEA reported in the literature; a recent survey on MOEA and their application in aeronautical and aerospace engineering has been made in [16].
In the study, except for the MOPSO, we also test two other mostly popular MOEA. The first is the NSGAII algorithm which is proposed by Deb et al. [17]. This algorithm uses the idea of transforming the objectives to a single fitness measure by the creation of a number of fronts, sorted according to nondomination. During the fitness assignment, the first front is created as the set of solutions that is not dominated by any solutions in the population. These solutions are given the highest fitness and temporarily removed from the population, then a second nondominated front consisting of the solutions that are now nondominated is built and assigned the secondhighest fitness, and so forth. This is repeated until each of the solutions has been assigned a fitness. After each front has been created, its members are assigned crowding distances (normalized distance to closest neighbors in the front in the objective space) later to be used for niching. The NSGAII has been successfully applied in spacecraft trajectory optimization, for example, in designing a threeobjective impulse rendezvous problem [20, 21] and a twoobjective robust rendezvous problem with considering uncertainty [22].
The second is the SPEAII proposed by Zitzler et al. [18]. It uses an archive containing nondominated solutions previously found (the socalled external nondominated set). At each generation, nondominated individuals are copied to the external nondominated set, removing the dominated solutions. For each individual in this external set, a strength value is computed. Pareto dominance is adopted to ensure that the solutions are properly distributed along the Pareto front. It also uses a nearest neighbor density estimation technique and a finegrained fitness assignment strategy which guide the search more efficiently.
3.2. Brief Description of the MOPSO
The MOPSO applied in this study is the algorithm proposed by Pulido and Coello [23], which was competitive against the most popular MOEA such as the NSGAII, the PAES, and other MOPSO on typical benchmark problems, under common performance metrics [23]. The source code of the MOPSO is available from the EMOO repository located at http://delta.cs.cinvestav.mx/~ccoello/EMOO.
The MOPSO is based on the use of Pareto ranking and a subdivision of decision variable space into several subswarms which is done using clustering techniques. The complete execution process of this algorithm can be divided into three stages: initialization, flight, and generation of results [23].
At the first stage, every swarm is initialized. Each swarm creates and initializes its own particles and generates the leaders set among the particle swarm set by using Pareto ranking. In the second stage, it firstly performs the execution of the flight of every swarm; next, it applies a clustering algorithm to group the guide particles. This is performed until reaching a total of GMax iterations. The execution of the flight of each swarm can be seen as an entire PSO process (with the difference that it will only optimize a specific region of the search space). First, each particle will select a leader to which it will follow. At the same time, each particle will try to outperform its leader and to update its position. If the updated particle is not dominated by any member of the leaders set, then it will become a new leader. The execution of the swarm will start again until a total of sgmax iterations are reached. Once all the swarms have finished theirs flights, a clustering algorithm takes the control by grouping the closest particle guides into swarms. Each resulting group will be assigned to a different swarm. The third and final stage will present all the nondominated solutions found.
Details of this algorithm can be found in [23] and its pseudocode code is shown in Algorithm 1.

The MOPSO algorithm requires the following parameters: (1) GMax: the total number of generations that the algorithm will be executed; (2) sgmax: the number of internal generations that the particles of each swarm will run before sharing their leaders; (3) : the total number of particles; and (4) : the number of particle groups.
3.3. Constraints Optimization Method
The multiobjective asteroid mission design problem is a highly constrained problem whose constraints are described in (9). The simulation experiments show that these constraints strongly influence the convergence. In the MOPSO, the constraints are always handled in checking Pareto dominance [23, 24]. When we compare two individuals, we first check their feasibility. If they are both feasible, then the comparison is done using Pareto dominance. If one is feasible and the other is infeasible, the feasible individual wins. If both are infeasible, then the individual with the lower amount of total constraint violation wins.
The total constraint is calculated by making use of a nondifferentiable penalty function. For the general constrained problem in (10)–(13), the penalized, total constraint function is
In the present work is used as a penalty coefficient.
For most optimization problems, each design variable has its own upper and low values. Thus, a strategy to maintain the particles within the search space in case they go beyond their boundaries is necessary for the MOPSO algorithm.
In [23, 24], the following strategy was employed. When a decision variable goes beyond its boundaries, the decision variable takes the value of its corresponding boundary (either the lower or the upper boundary). The strategy will be likely effective when the Pareto solutions are located in the boundary of variables.
However, our simulation experiments show that this strategy is not very effective in solving our problems. Therefore, another simple strategy is proposed. When a decision variable goes beyond its boundaries, the decision variable takes a random value from its feasible design space. The probabilitybased disposal to boundary constraint could enrich the diversity of swarm flight.
The MOPSO with these two different boundary constraint optimization strategies is, respectively, called as MOPSOI and MOPSOII. Their performance is compared in Section 4.3.
3.4. Algorithm Assessment Metrics
In order to allow a quantitative assessment of the performance of the MSOPO, we adopted the following two metrics.
The first one is the epsilon indicator [25]. Given a reference Pareto front (ideally the true Pareto front, if available), the epsilon indicator measures the minimum amount necessary to translate all the points of the found Pareto front to weakly dominate the reference set. The epsilon indicator has two types, that is, the additive type and the multiplicative type, and the multiplicative type is used here.
The second one is the hypervolume indicator [26]. The hypervolume of a set of solutions measures the size of the portion of objective space that is dominated by those solutions collectively. Generally, hypervolume is favored because it captures in a single scalar both the closeness of the solutions to the optimal set and, to some extent, the spread of the solutions across objective space.
The statistical test chosen for result evaluation is the MannWhitney test [27]. This is a nonparametric rankbased test that can be used to compare two independent sets of sampled data. It outputs values that estimate the probability of a failure to reject the null hypothesis of the study question. Here the values can be interpreted as the probability that the performance of one algorithm is different (superior or junior) with statistical significance to that of the other.
4. Simulation Results
4.1. Problem Configuration
In order to testify the effectiveness of the proposed method, two different asteroid mission designs are illustrated. Table 1 lists the orbit elements of the two asteroids. The upper and lower space of the optimization variables are provided in Table 2. Three different cases with the number of impulses of 2, 3, and 4 are, respectively, tested for each asteroid mission.


4.2. Pareto Fronts Analysis
From our experiments, this multiobjective asteroid rendezvous design problem is very difficult to be solved, and the true Pareto fronts of this multiobjective problem are not known. In order to obtain the Pareto fronts as possible close to the true ones, the MOPSO is executed with a much larger number of function evaluations. The parameters of the MOPSO are , , , and .
Considering the stochastic characteristic of the MOPSO, 10 independent runs for each test case are completed. All the Pareto solutions of the 10 independent runs are compared, and the repeated and nonPareto solutions are deleted, and the revised Pareto solutions are selected as the final solutions. In the following examples, the figured Pareto fronts are all obtained in the same method. The MOPSO with different boundary constraint optimization methods, that is, the MOPSOI and MOPSOII, is both tested. The performance comparisons between the MOPSOI and MOPSOII will be analyzed in the next section.
Figure 1 compares the Pareto solution front for twoimpulse, threeimpulse, and fourimpulse of the Asteroid 1999YR14 mission. Figure 2 compares this for the Asteroid 2340 mission. The tradeoffs between the total and the transfer time are clearly demonstrated by Figures 1 and 2, which will be useful for a mission designer. Several inherent peculiarities regarding the optimal multiobjective asteroid rendezvous trajectories have been observed from the optimization results. For the second asteroid, the corresponding total reduces evidently as the transfer time increases when the transfer time is in the range of day, but small change in other ranges; therefore the Pareto front concentrates on a narrow range of 40 days. For the first asteroid, the corresponding total reduces evidently as the transfer time increases in a large range of day; therefore the Pareto solutions distribute in a much larger space, but the Pareto front is discontinuous in small range as demonstrated in Figure 1.
(a) Produced by MOPSOI
(b) Produced by MOPSOII
(a) Produced by MOPSOI
(b) Produced by MOPSOII
In order to explain this phenomenon, the propellantoptimal solutions are obtained by using the approach employed in [28]. In this optimization, the transfer time is fixed and is calculated every 5 days in the transfer time of range days. The relations between the transfer time and of the propellantoptimal solutions are illustrated in Figure 3. Some transfer time ranges with muchhigher propellant cost are located; for example, the is about 9.5 km/s for the transfer time of 300 days, while being only 6.2 km/s for the transfer time of 280 days. The former solution is larger than the latter solution by about 50% in ; besides its transfer time is also larger than the latter. This issue can explain why the Pareto fronts of twoimpulse solutions discontinue near the transfer time of 300 days.
The influence of different number of impulses on the total has been also demonstrated by Figures 1 and 2. It is obvious that a threeimpulse trajectory and a fourimpulse trajectory cost less propellant in compared with a twoimpulse trajectory. It is necessary to increase the number of impulse to find a much better trajectory in the preliminary asteroid rendezvous mission design.
4.3. Performance Analysis of the MOPSO
4.3.1. Boundary Constraint Optimization Comparisons
As seen from Figures 1 and 2, the boundary constraint optimization method of the MOPSO has a great influence on the Pareto fronts. Graphically, the MOPSOII has much better performance in diversity, and the spread of its solutions found is much larger than that of the MOPSOI. This can easily explain that probabilitybased disposal to boundary constraint can enrich the diversity of swarm flight.
It is not easily to determine which one is much closer to the true Pareto fronts from Figure 1, and the twoimpulse comparison case is shown in Figure 4(a), while it is easy to see from Figure 2 that the Pareto fronts obtained by MOPSOI are much closer to the true Pareto fronts, and the threeimpulse comparison case is redrawn in Figure 4(b) to demonstrate this clearly.
(a) Asteroid 1, twoimpulse
(b) Asteroid 2, threeimpulse
Therefore, the quantitative metrics are calculated and the statistical results are provided in Table 3. As seen from Table 3, the MOPSOII has better performance with respect to the hybervolume indicator (larger value means better) for the two cases. While, for the epsilon indicator (smaller value means better), the MOPSOII is much better in case of Asteroid 1, twoimpulse problem and the MOPSOI are much better in case of Asteroid 2, threeimpulse.

Through the comparisons provided in Table 2 and Figure 4, our proposed boundary constraint demonstrates much better performance in total. Besides, these comparisons show that the same algorithm will have different performance in solving the same type of problem with different configurations.
4.3.2. Algorithm Parameters Analysis
Our simulation experiments show that and are the two main parameters affecting the performance of the MOPSO in solving this complex multiobjective problem. In order to quantitatively evaluate their influence , is, respectively, set as 400, 200, and 100, is also, respectively, set as 400, 200, and 100, and a total of nine groups of parameters (in Table 4) are furthered tested. Other parameters are the same as used in Section 4.2.

Ten independent runs for the MOPSOII with each group of parameters in solving the threeimpulse, Asteroid 1 problem, are executed.
Figure 5(a) compares the Pareto fronts of Cases 1, 5, and 9, Figure 5(b) compares that of Cases 2, 5, and 8, and Figure 5(c) compares that of Cases 7, 8, and 9. These comparisons are to show the influence of the total number of function evaluations. It is clearly seen from Figure 5 that the larger number of function evaluations evidently increases the performance of the MOPSO. Thus, a larger number of swarm size and iteration are necessary for obtaining the optimal Pareto fronts for this practical multiobjective optimization problem.
(a) Cases 1, 5, and 9
(b) Cases 2, 5, and 8
(c) Cases 7, 8, and 9
Figure 6(a) compares the Pareto fronts of Cases 2 and 4, Figure 6(b) compares that of Cases 6 and 8, and Figure 6(c) compares that of Cases 3, 5, and 7. These comparisons are to show the influence of and with the same number of function evaluations. Evidently, we cannot determine which case is better by the graphical results provided in Figure 6. Therefore, the quantitative metrics are calculated and the statistical results are provided in Tables 5, 6, 7, and 8.



 
Bold indicates epsilon; the other is hybervolume. 
(a) Cases 2, 4
(b) Cases 6, 8
(c) Cases 3, 5, and 7
The comparisons between Case 2 and Case 4 show that larger size of swarm can obtain better average epsilon indicator with a value of 67.76% but worse average hybervolume indicator with a value of 14.05%. The comparisons between Case 6 and Case 8 show that larger size of swarm both increase average performance in epsilon and hybervolume, with a value of 47.27% and 7.57%. By comparing Cases 3, 5, and 7, Case 5 demonstrates the best performance in total. However, the improvements are not evident as the values are all a large value.
Although the result obtained from this experiment seems to be inconclusive on which parameter is the best, we can argue that the MOPSO’s performance is not sensitive to the main PSO parameters, under the condition that the total number of function evaluation retains a large value. We have obtained competitive results in most cases without paying special attentions on PSO parameters.
4.4. Comparison of the MOPSO with Other MOEA
To demonstrate the performance of the employed MOPSO, two other popular algorithms are tested and compared: the NSGAII and the SPEAII.
In the following examples, the total number of function evaluations was set to 20000 for all the algorithms compared. The NSGAII and the SPEAII were run using a population size of 200, a maximum number of generations of 100. The MOPSOI and MOPSOII used 200 particles, a maximum number of generations of 20, a maximum number of generations per swarm of 5, and a total of 5 swarms. The source code of the NSGAII, and SPEAII provided in the EMOO repository is also employed in this study.
The two test problems are the same to those employed in Section 4.3 to compare boundary constraint optimization. Ten independent runs for each test case are executed. The final Pareto fronts are compared in Figure 7 and the statistical results with respect to epsilon and hybervolume indicators are provided in Tables 9, 10, and 11.

 
Bold indicates epsilon; the other is hybervolume. 
 
Bold indicates epsilon; the other is hybervolume. 
(a) Asteroid 1, twoimpulse
(b) Asteroid 2, threeimpulse
Evidently seen from Figure 7(a), the MOPSOI produces better Pareto fronts than the NSGAII and SPEAII for the first test problem. The statistical results of the two indicators also support this. The value with respect to epsilon is calculated as 0.58% and 1.13%, respectively, and 0.73% and 0.17% for hybervolume. This shows that the MOPSOI is absolutely super to the NSGAII and SPEAII.
Also evidently seen from Figure 7(b), the MOPSOI produces much closer Pareto fronts than the NSGAII and SPEAII for the second test problem. However, its Pareto fronts are a little narrower, which may result in that its epsilon and hybervolume indicator are junior to the NSGAII and SPEAII. Furthermore, the fact that the Pareto fronts are highly discontinuous would make these metrics irrelevant. In the view of practical spacecraft mission design, obtaining much closer Pareto fronts from which the designer can choose one single satisfying solution for engineering design would be much desirable performance for a multiobjective optimizer. Thus, we can say that the MOPSOI is a little better than the NSGAII and SPEAII; even the statistical results on quantitative metrics do not support this.
In conclusion, the MOPSO showed a better performance with respect to the NSGAII and the SPEAII for the two test problems. For other test problems, the MOPSO is not always better, but its performance is very competitive.
4.5. Comparison with PCP Method
The widely used PorkChop (PCP) method is also tested here for comparison. The is calculated every 1 day in both Earth departure time and asteroid arrival time, and the twoimpulse Lambert algorithm is employed. The contours of time of flight corresponding to the optimal transfers for Asteroid 1 are illustrated in Figure 8. For convenience, only the solutions with a less than 18 km/s and Earth departure and asteroid arrive times in (MJD2000) are presented.
(a) Asteroid 1, twoimpulse
(b) Asteroid 2, twoimpulse
We analyze the comparisons between the proposed multiobjective optimization approach and the PCP method in the following three aspects.
4.5.1. Uses of the Results
In the PCP method, the contours of the are plotted to assist the designer to find the best launch window and transfer trajectory. However, the contours have many local peaks and it is not very convenient to determine which solution is the best one. Besides, the size of the region containing the global minimum is a very important factor in finding real trajectories in any given year. This is not easily observed from the PorkChop figures.
In contrast to the PCP method, the proposed multiobjective optimization approach can provide friendly the designer of this information. The relations between total characteristic velocity, transfer time, and earth departure time of those obtained Paretooptimal solutions for two test cases are provided in Figure 8. From Figure 8, it can be clearly observed that the earth departure time for all the optimal rendezvous trajectories concentrates on one arrow domain for these two cases. The optimized Earth launch window and its size can be easily determined through the multiobjective optimization design, which will provide very useful information for engineering design.
4.5.2. Global Optimality
The minimumpropellant solution searched by the PCP method for the first asteroid mission is calculated with a of 6015.3 m/s and it is located in Figure 9 (* denotes). As seen in Figure 1, the MOPSO has located a set of Pareto solutions with a about this value. This issue can also demonstrate the global convergence ability of the MOPSO. Besides, the PCP method obtains only the twoimpulse trajectory. As is well known, the twoimpulse trajectory is not the propellantoptimal solution under most conditions, and increasing the number of impulses will reduce the propellant cost. Our proposed approach can design the multiimpulse trajectory; thus it locates better solution than the PCP method.
4.5.3. Efficiency
The PCP is in essence of an exhaustive searching method. In our test, the search space for departure time and arrive time is day. With the search step of 1 day, the PCP method calculates a total number of trajectories, while only trajectories are calculated in the proposed method for the most highly cost case. The proposed approach improves the calculation efficiency by about 40 times. Besides, the PCP method needs much human intention to determine the final solutions, while the proposed approach is an automated search method.
5. Conclusions
The paper formulates the asteroid rendezvous preliminary trajectory design as a multiobjective optimization problem and employs the multiobjective particle swarm optimization (MOPSO) algorithm to locate the Paretooptimal solution set. Compared with the widely employed PorkChop method, the proposed approach is demonstrated to be able to provide much more easily used results, obtain better propellantoptimal solutions, and have much better efficiency. The results show that the proposed approach can effectively and efficiently demonstrate the relations among the mission characteristic parameters such as launch time, transfer time, propellant cost, and number of maneuvers, which will provide useful reference for practical asteroid mission design. The MOPSO proves to be quite effective in finding the Paretooptimal solutions and its performance can be improved by a proposed boundary constraint optimization strategy. The MOPSO is found to be very competitive with respect to two highly competitive multiobjective evolutionary algorithms: the NSGAII and the SPEAII.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (no. 11222215), the 973 Project (no. 2013CB733100), and the Hunan Provincial Natural Science Foundation of China (no. 13JJ1001).
References
 N. D. Hulkower, C. O. Lau, and D. F. Bender, “Optimum twoimpulse transfer for preliminary interplanetary trajectory design,” Journal of Guidance, Control, and Dynamics, vol. 7, no. 4, pp. 458–461, 1984. View at: Google Scholar
 C. O. Lau and N. D. Hulkower, “Accessibility of nearearth asteroids,” Journal of Guidance, Control, and Dynamics, vol. 10, no. 3, pp. 225–232, 1987. View at: Google Scholar
 E. Perozzi, A. Rossi, and G. B. Valsecchi, “Basic targeting strategies for rendezvous and flyby missions to the nearearth asteroids,” Planetary and Space Science, vol. 49, no. 1, pp. 3–22, 2001. View at: Publisher Site  Google Scholar
 D. Qiao, H. T. Cui, and P. Y. Cui, “Evaluating accessibility of nearearth asteroids via earth gravity assists,” Journal of Guidance, Control, and Dynamics, vol. 29, no. 2, pp. 502–505, 2006. View at: Publisher Site  Google Scholar
 D. F. Lawden, Optimal Trajectories for Space Navigation, Butterworths, London, UK, 1963. View at: MathSciNet
 J. E. Pussing, “Primer vector theory and applications,” in Spacecraft Trajectory Optimization, B. A. Conway, Ed., pp. 16–36, Cambridge University Press, New York, NY, USA, 2010. View at: Google Scholar
 J. Kennedy and R. Eberhart, “Particle swarm optimization,” in Proceedings of the IEEE International Conference on Neural Networks, vol. 4, pp. 1942–1948, IEEE, Piscataway, NJ, USA, December 1995. View at: Google Scholar
 C. R. Bessette and D. B. Spencer, “Identifying optimal interplanetary trajectories through a genetic approach,” in Proceedings of the AIAA/AAS Astrodynamics Specialist Conference and Exhibit, AIAA Paper 20066306, pp. 809–826, Keystone, Colo, USA, August 2006. View at: Google Scholar
 K. Zhu, F. Jiang, J. Li, and H. Baoyin, “Trajectory optimization of multiasteroids exploration with low thrust,” Transactions of the Japan Society for Aeronautical and Space Sciences, vol. 52, no. 175, pp. 47–54, 2009. View at: Publisher Site  Google Scholar
 M. Pontani and B. A. Conway, “Particle swarm optimization applied to space trajectories,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 5, pp. 1429–1441, 2010. View at: Publisher Site  Google Scholar
 M. Pontani, P. Ghosh, and B. A. Conway, “Particle swarm optimization of multipleburn rendezvous trajectories,” Journal of Guidance, Control, and Dynamics, vol. 35, no. 4, pp. 1192–1207, 2012. View at: Publisher Site  Google Scholar
 M. ReyesSierra and C. A. C. Coello, “Multiobjective particle swarm optimizers: a survey of the stateoftheart,” International Journal of Computational Intelligence Research, vol. 2, no. 3, pp. 287–308, 2006. View at: Google Scholar  MathSciNet
 K. B. Lee and J. H. Kim, “Multiobjective particle swarm optimization with preferencebased sort and its application to path following footstep optimization for humanoid robots,” IEEE Transactions on Evolutionary Computation, vol. 17, no. 6, pp. 755–766, 2013. View at: Google Scholar
 X. Y. Shao, W. Q. Liu, Q. Liu, and C. Y. Zhang, “Hybrid discrete particle swarm optimization for multiobjective flexible jobshop scheduling problem,” International Journal of Advanced Manufacturing Technology, vol. 67, no. 9–12, pp. 2885–2901, 2013. View at: Google Scholar
 M. Shokrian and K. A. High, “Application of a multi objective multileader particle swarm optimization algorithm on NLP and MINLP problems,” Computers & Chemical Engineering, vol. 60, pp. 57–75, 2014. View at: Publisher Site  Google Scholar
 A. AriasMontaño, C. A. C. Coello, and E. MezuraMontes, “Multiobjective evolutionary algorithms in aeronautical and aerospace engineering,” IEEE Transactions on Evolutionary Computation, vol. 16, no. 5, pp. 662–694, 2012. View at: Publisher Site  Google Scholar
 K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGAII,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002. View at: Publisher Site  Google Scholar
 E. Zitzler, M. Laumanns, and L. Thiele, “SPEA2: improving the strength Pareto evolutionary algorithm,” in Evolutionary Methods for Design, Optimization and Control with Applications to Industrial Problems, K. Giannakoglou, D. T. Tsahalis, J. Périaux, K. D. Papailiou, and T. Fogarty, Eds., pp. 95–100, Springer, Berlin, Germany, 2002. View at: Google Scholar
 S. P. Hughes, L. M. Mailhe, and J. J. Guzman, “A comparison of trajectory optimization methods for the impulsive minimum fuel rendezvous problem,” in Guidance and Control 2003, I. J. Gravseth and R. D. Culp, Eds., vol. 113 of Advances in the Astronautical Sciences, pp. 85–104, 2003. View at: Google Scholar
 Y. Z. Luo, G. J. Tang, and Y. J. Lei, “Optimal multiobjective nonlinear impulsive rendezvous,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 4, pp. 994–1002, 2007. View at: Publisher Site  Google Scholar
 Y. Z. Luo, G. J. Tang, and Y. J. Lei, “Optimal multiobjective linearized impulsive rendezvous,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 2, pp. 383–389, 2007. View at: Publisher Site  Google Scholar
 Y. Z. Luo, Z. Yang, and H. N. Li, “Robust optimization of nonlinear impulsive rendezvous with uncertainty,” Science China: Physics, Mechanics & Astronomy, vol. 57, no. 3, pp. 1–10, 2014. View at: Publisher Site  Google Scholar
 G. T. Pulido and C. A. C. Coello, “Using clustering techniques to improve the performance of a particle swarm optimizer,” in Genetic and Evolutionary Computation—GECCO 2004, vol. 3102 of Lecture Notes in Computer Science, pp. 225–237, Springer, Berlin, Germany, 2004. View at: Publisher Site  Google Scholar
 C. A. C. Coello, G. T. Pulido, and M. S. Lechuga, “Handling multiple objectives with particle swarm optimization,” IEEE Transactions on Evolutionary Computation, vol. 8, no. 3, pp. 256–279, 2004. View at: Publisher Site  Google Scholar
 E. Zitzler, L. Thiele, M. Laumanns, C. M. Fonseca, and V. G. da Fonseca, “Performance assessment of multiobjective optimizers: an analysis and review,” IEEE Transactions on Evolutionary Computation, vol. 7, no. 2, pp. 117–132, 2003. View at: Publisher Site  Google Scholar
 E. Zitzler and L. Thiele, “Multiobjective optimization using evolutionary algorithms—a comparative case study,” in Parallel Problem Solving from Nature—PPSN V, A. E. Eiben, Ed., vol. 1498 of Lecture Notes in Computer Science, pp. 292–301, Springer, Amsterdam, The Netherlands, 1998. View at: Publisher Site  Google Scholar
 M. Hollander and D. A. Wolfe, Nonparametric Statistical Methods, Wiley Series in Probability and Statistics, John Wiley & Sons, Hoboken, NJ, USA, 2nd edition, 1999. View at: MathSciNet
 Y. Z. Luo, J. Zhang, H. Y. Li, and G. J. Tang, “Interactive optimization approach for optimal impulsive rendezvous using primer vector and evolutionary algorithms,” Acta Astronautica, vol. 67, no. 34, pp. 396–405, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Yazhong Luo and Lini Zhou. 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.