5th Global Trajectory Optimisation Competition

22 downloads 0 Views 436KB Size Report
Nov 1, 2010 - GP2). 3359.45. 3122.40. 59827.12. 60120.57. 293.45. 2.90. 293. 6. (2001. GP2). (2001. GP2). 3082.40. 3005.40. 60120.57. 60222.60. 102.03.
5th Global Trajectory Optimisation Competition: Results of Team 5 Dario Izzo, Francesco Biscani, Chit Hong Yam, Luis Felismino Simoes Advanced Concepts Team, European Space Agency, The Netherlands

David di Lorenzo, Bernadetta Addis, Andrea Cassioli, Fabio Schoen, Marco Locatelli Global Optimization Laboratories, University of Florence

November 1, 2010

1 Introduction and Method The solution proposed to this year problem for the fth edition of the Global Trajectory Optimization Competition was obtained by us using a two-step process. 1. We build an ad hoc analytical model for the optimal minimum time transfer from an orbiting object to itself constraining the relative re-encounter velocity to a minimum of 400 m/s. The model returns an analytical estimate for both time of ight and mass spent. The linear model is then used in conjunction to a multiple revolution lambert solver to search for good 'chemical representations' of trajectories. A branch and prune approach is adopted to search the immense tree of resulting possibilities. 2. The chemical trajectories candidates are then turned into low-thrust by keeping epochs and asteroid sequence as initial guess to our global optimization methods for constrained optimization problems. We work directly on a high delity trajectory model constructed using a variant to the sims-anagan transcription and applying monotonic basin hopping as main solving paradigm.

2 Solution Overview Our nal solution accumulates a total primary score of J = 16.0 with a secondary objective τ = 14.19 years - The launch V∞ = 2.35 km/s the nal mass mf = 500 kg (nal penetrator mass accounted). A summary of the complete trajectory is given in the table below. The visualization of the ecliptic projection of all inter asteroid phases is also reported in the gures below.

1

2

phase start asteroid end asteroid start mass end mass start epoch end epoch 1 Earth (2007 UN12) 4000.00 3874.24 59062.81 59241.23 2 (2007 UN12) (2007 UN12) 3834.24 3736.18 59241.23 59366.70 (2007 UN12) (2006 JY26) 3735.18 3487.07 59366.70 59718.72 3 4 (2006 JY26) (2006 JY26) 3447.07 3360.45 59718.72 59827.12 5 (2006 JY26) (2001 GP2) 3359.45 3122.40 59827.12 60120.57 6 (2001 GP2) (2001 GP2) 3082.40 3005.40 60120.57 60222.60 7 (2001 GP2) 225312 (1996 XB27) 3004.40 2706.22 60222.60 60645.66 8 225312 (1996 XB27) 225312 (1996 XB27) 2666.22 2610.63 60645.66 60768.01 9 225312 (1996 XB27) (2003 SM84) 2609.63 2491.05 60768.01 61044.11 10 (2003 SM84) (2003 SM84) 2451.05 2391.03 61044.11 61135.71 11 (2003 SM84) (2006 HE2) 2390.03 2227.29 61135.71 61410.80 12 (2006 HE2) (2006 HE2) 2187.29 2132.20 61410.80 61490.44 13 (2006 HE2) (2007 HC) 2131.20 2006.92 61490.44 61679.56 14 (2007 HC) (2007 HC) 1966.92 1916.45 61679.56 61753.71 15 (2007 HC) (1998 KY26) 1915.45 1734.75 61753.71 62127.99 16 (1998 KY26) (1998 KY26) 1694.75 1649.67 62127.99 62192.70 17 (1998 KY26) (2009 TD17) 1648.67 1552.93 62192.70 62383.79 18 (2009 TD17) (2009 TD17) 1512.93 1473.12 62383.79 62440.28 19 (2009 TD17) (2007 YF) 1472.12 1345.03 62440.28 62647.84 20 (2007 YF) (2007 YF) 1305.03 1271.97 62647.84 62697.49 21 (2007 YF) (1999 AO10) 1270.97 1190.78 62697.49 62949.99 22 (1999 AO10) (1999 AO10) 1150.78 1127.17 62949.99 63002.23 23 (1999 AO10) (2006 RH120) 1126.17 1027.11 63002.23 63303.21 24 (2006 RH120) (2006 RH120) 987.11 964.00 63303.21 63344.09 25 (2006 RH120) (2008 UA202) 963.00 926.63 63344.09 63481.71 26 (2008 UA202) (2008 UA202) 886.63 867.50 63481.71 63522.77 27 (2008 UA202) (2008 HU4) 866.50 770.33 63522.77 63874.80 28 (2008 HU4) (2008 HU4) 730.33 713.20 63874.80 63905.35 29 (2008 HU4) (2004 VJ1) 712.20 657.25 63905.35 64061.71 30 (2004 VJ1) (2004 VJ1) 617.25 603.89 64061.71 64090.16 31 (2004 VJ1) (2003 WT153) 602.89 554.57 64090.16 64225.14 (2003 WT153) (2003 WT153) 514.57 501.00 64225.14 64244.62 32

leg tof tot tof thrust duration 178.43 0.49 140 125.47 0.83 125 352.02 1.80 280 108.41 2.09 108 293.45 2.90 293 102.03 3.18 102 423.06 4.33 337 122.35 4.67 70 276.10 5.42 162 91.60 5.68 71 275.09 6.43 218 79.64 6.65 62 189.12 7.16 150 74.15 7.37 58 374.28 8.39 374 64.71 8.57 50 191.10 9.09 111 56.49 9.25 43 207.56 9.82 164 49.65 9.95 37 252.50 10.64 144 52.23 10.79 24 300.99 11.61 140 40.88 11.72 22 137.62 12.10 32 41.06 12.21 17 352.03 13.17 167 30.55 13.26 15 156.36 13.69 54 28.45 13.76 12 134.99 14.13 59 19.47 14.19 13

phases 2 & 3 1.0

0.5

0.5

y (AU)

y (AU)

phases 1 & 2 1.0

0.0

0.5

0.0

0.5

1.0

1.0

1.5

1.0

0.5

0.0 x (AU)

0.5

1.5

1.0

1.5

phases 3 & 4

1.0

0.5

0.0 x (AU)

0.5

1.5

1.0

phases 4 & 5

1.5

1.0 1.0 0.5

y (AU)

y (AU)

0.5 0.0

0.0 0.5 0.5 1.0 1.5

1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

1.0

phases 5 & 6

1.5

1.5

1.0

0.5

x (AU)

0.0

0.5

1.0

phases 6 & 7

1.5 1.0

1.0

0.5 y (AU)

y (AU)

0.5 0.0

0.0 0.5 0.5

1.0

1.0

1.5

1.0

x (AU)

0.0

0.5

1.5

1.0

phases 7 & 8

1.5 1.0

1.0

0.5

0.5

0.0

0.5

1.0

1.0

1.5

1.0

0.5

0.0 x (AU)

1.0

0.5

1.0

1.5

1.5

3

0.5

0.0 x (AU)

0.5

1.0

1.5

0.5

1.0

1.5

phases 8 & 9

0.0

0.5

1.5

1.5

1.5

y (AU)

y (AU)

0.5

1.5

1.0

0.5

0.0 x (AU)

phases 9 & 10

1.0

1.0

0.5

0.5

0.0

0.0

0.5

0.5

1.0

1.0

1.5

1.5

1.0

0.5

0.0 x (AU)

phases 10 & 11

1.5

y (AU)

y (AU)

1.5

0.5

1.5

1.5

1.0

1.5

1.0

0.5

0.0 x (AU)

0.5

1.5

1.0

(i) phases 11 & 12

1.5

phases 12 & 13

1.5

1.0

1.0

0.5 y (AU)

y (AU)

0.5 0.0

0.0 0.5 0.5

1.0 1.5

1.5

1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

1.0

phases 13 & 14

1.5

1.5

1.0

0.5

0.5

1.0

1.5

0.5

1.0

1.5

y (AU)

0.5 y (AU)

1.0

0.5

phases 14 & 15

1.5

1.0

0.0 x (AU)

0.0

0.0

0.5

0.5

1.0

1.5

1.0

0.5

0.0 x (AU)

0.5

1.0

1.0

1.5

phases 15 & 16

1.5

1.5

1.0

0.5 y (AU)

0.5 y (AU)

1.0

0.0 x (AU)

phases 16 & 17

1.5

1.0

0.5

0.0

0.0

0.5

0.5

1.0

1.5

1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

4

1.0

1.5

1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

phases 17 & 18

1.5

phases 18 & 19

1.5

0.5

0.5 y (AU)

1.0

y (AU)

1.0

0.0

0.0

0.5

0.5

1.0

1.5

1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

1.0

phases 19 & 20

1.5

1.5

1.0

0.5

0.5

0.5

1.0

0.5

1.0

y (AU)

0.5 y (AU)

1.0

0.0

phases 20 & 21

1.5

1.0

x (AU)

0.0

0.0

0.5

0.5

1.0

1.5

1.0

0.5

x (AU)

0.0

0.5

1.0

1.0

1.5

1.0

0.5

phases 21 & 22

1.5

x (AU)

0.0

phases 22 & 23 1.0

1.0 0.5

y (AU)

y (AU)

0.5 0.0

0.0 0.5 0.5 1.0 1.0

1.5

1.0

0.5

x (AU)

0.0

0.5

1.5

1.0

1.0

0.5

1.0

0.5

0.5

0.0

0.5

1.0

1.5

0.5

1.0

1.5

0.0

0.5

1.0 1.5

0.5

phases 24 & 25

1.0

y (AU)

y (AU)

phases 23 & 24

0.0 x (AU)

1.0 1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

5

1.5

1.0

0.5

0.0 x (AU)

phases 26 & 27 1.0

0.5

0.5

y (AU)

y (AU)

phases 25 & 26 1.0

0.0

0.5

0.5

1.0 1.5

0.0

1.0 1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

1.5

1.0

0.5

1.0

0.5

0.5

0.0

0.5

1.0

1.5

0.5

1.0

1.5

0.0

0.5

1.0 1.5

0.5

phases 28 & 29

1.0

y (AU)

y (AU)

phases 27 & 28

0.0 x (AU)

1.0 1.0

0.5

0.0 x (AU)

0.5

1.0

1.5

1.5

phases 29 & 30

1.0

0.5

0.0 x (AU)

phases 30 & 31

1.5

1.0 1.0 0.5

y (AU) 0.0

0.5 0.5 1.0 1.5

1.0

0.5

0.0 x (AU)

0.5

1.0

1.0

1.5

1.0

0.5

phases 31 & 32

1.5

1.0

0.5 y (AU)

y (AU)

0.5 0.0

0.0

0.5

1.0

1.0

0.5

0.0

x (AU)

6

0.5

1.0

1.5

0.0

x (AU)

0.5

1.0

1.5

complete trajectory

1.5 1.0

y (AU)

0.5 0.0 0.5 1.0 1.5

1.5

1.0

0.5

0.0 x (AU)

7

0.5

1.0

1.5

1.5 1.4

distance (AU)

1.3 1.2 1.1 1.0 0.9 0.8 0.7 59000

60000

61000

62000 Epoch (MJD)

8

63000

64000

65000