Informatics
Mission Analysis
May 7, 2013

SAGAS

A visualization of the best solution found
A visualization of the best solution found

In this trajectory problem we design what is commonly called a deltaV-EGA manouvre to then fly-by Jupiter and reach 50AU. The objective function considered is the overall mission length and has to be minimsed. This creates an MGADSM problem where two more variables need to be added to the decision vector in order to be able to evaluate the keplerian orbit reached after the last fly-by.


The code

  1. MATLAB: use the function sagas.m and pass to it the MGAproblem variable contained in sagas.mat
  2. C++: call the function "double sagas (const std::vector & x, double& DVtot, double& DVonboard)" provided in the GTOPtoolbox.
  3. C++ (PaGMO): use the class pagmo::problem::sagas
  4. Python 2.7 (PyGMO): use PyGMO.problem.sagas().obj_fun(x)

Problem Description

The box bounds on each of the decision vector variable are given below.

State

Variable

Lower Bounds

Upper Bounds

Units

x(1)t070009100MJD2000
x(2)Vinf07km/sec
x(3)u01n/a
x(4)v01n/a
x(5)T1502000days
x(6)T23002000days
x(7)eta10.010.9n/a
x(8)eta20.010.9n/a
x(9)r_p11.057n/a
x(10)r_p28500n/a
x(11)b_incl1-pipirad
x(12)b_incl2-pipirad

We also consider the following two contraints:

  1. deltaV_1 + deltaV_2 < 1.782
  2. deltaV:=deltaV_1 + delta V_2 + delta V_infinity < 6.782 km/s

which capture the limited on board fuel and the launcher performance.

Note that the bound on the departure deltaV is quite large and include a very strong minima at around 1-4 km/sec (1:1 Earth orbit resonance) that often tricks the optimisers. Clearly, by reducing this bound (knowledge-based pruning) one can drastically help any optimiser to locate the correct global optima. As here we are interested in the algorithmic performances, we selected the bounds to create an interesting optimisation problem rather than to simplify the problem as much as possible and thus we included zones of the search space one could prune out by experience.

The objective function is considered to the precision of 0.01 years.


Solutions

The solutions submitted to this problem are listed below in chronological order

Author

Affiliation

Date

Obj Fun

Solution

Algorithm

T., Vinko, D., IzzoEuropean Space AgencyOctober 200518.19 years
x=[7020.49, 5.34817, 1, 0.498915, 788.763, 484.349, 0.4873, 0.01, 1.05, 10.8516, -1.57191, -0.685429]
Differential Evolution
Hamburger icon
Menu
Advanced Concepts Team