Continuation from a flat to a round Earth model in the coplanar orbit transfer problem
|
|
- Louise George
- 5 years ago
- Views:
Transcription
1 Continuation from a flat to a round Earth model in the coplanar orbit transfer problem M. Cerf 1, T. Haberkorn, Emmanuel Trélat 1 1 EADS Astrium, les Mureaux 2 MAPMO, Université d Orléans First Industrial Workshop, SADCO 2011, March 2nd
2 The coplanar orbit transfer problem Spherical Earth Central gravitational field g(r) = µ r 2 System in cylindrical coordinates ṙ(t) = v(t) sin γ(t) ϕ(t) = v(t) cos γ(t) r(t) v(t) = g(r(t)) sin γ(t) + Tmax m(t) u 1(t) v(t) γ(t) = r(t) g(r(t)) «cos γ(t) + Tmax v(t) m(t)v(t) u 2(t) ṁ(t) = βt max u(t) Thrust: T (t) = u(t)t max (T max large: strong thrust) Control: u(t) = (u 1 (t), u 2 (t)) satisfying u(t) = p u 1 (t) 2 + u 2 (t) 2 1
3 The coplanar orbit transfer problem Initial conditions r(0) = r 0, ϕ(0) = ϕ 0, v(0) = v 0, γ(0) = γ 0, m(0) = m 0, Final conditions a point of a specified orbit: r(t f ) = r f, v(t f ) = v f, γ(t f ) = γ f, or an elliptic orbit of energy K f < 0 and eccentricity e f : ξ Kf = v(t f ) 2 µ 2 r(t f ) K f = 0, ξ ef = sin 2 γ + 1 r(t f )v(t f ) 2 «2 cos 2 γ ef 2 = 0. µ (orientation of the final orbit not prescribed: ϕ(t f ) free; in other words: argument of the final perigee free) Optimization criterion max m(t f ) (note that t f has to be fixed)
4 Application of the Pontryagin Maximum Principle Hamiltonian H(q, p, p 0 v, u) = p r v sin γ + p ϕ cos γ + pv r + p γ v r g(r) v g(r) sin γ + Tmax m u 1 ««cos γ + Tmax mv u 2 p mβt max u, «Extremal equations q(t) = H p (q(t), p(t), p0, u(t)), ṗ(t) = H q (q(t), p(t), p0, u(t)), Maximization condition H(q(t), p(t), p 0, u(t)) = max w 1 H(q(t), p(t), p0, w)
5 Application of the Pontryagin Maximum Principle Hamiltonian H(q, p, p 0 v, u) = p r v sin γ + p ϕ cos γ + pv r + p γ v r g(r) v g(r) sin γ + Tmax m u 1 ««cos γ + Tmax mv u 2 p mβt max u, «Maximization condition leads to u(t) = (u 1 (t), u 2 (t)) = (0, 0) whenever Φ(t) < 0 u 1 (t) = r p v (t) p v (t) 2 + pγ (t)2 v(t) 2, u 2 (t) = v(t) r p γ(t) p v (t) 2 + pγ (t)2 v(t) 2 whenever Φ(t) > 0 where s Φ(t) = 1 p v (t) 2 + pγ(t)2 βpm(t) (switching function) m(t) v(t) 2
6 Application of the Pontryagin Maximum Principle Hamiltonian H(q, p, p 0 v, u) = p r v sin γ + p ϕ cos γ + pv r + p γ v r g(r) v g(r) sin γ + Tmax m u 1 ««cos γ + Tmax mv u 2 p mβt max u, «Transversality conditions case of a fixed point of a specified orbit: p ϕ(t f ) = 0, p m(t f ) = p 0 case of an orbit of given energy and eccentricity: r ξ Kf (p γ v ξ ef p v γξ ef ) + v ξ Kf (p r γξ ef p γ r ξ ef ) = 0 Remark p 0 0 (no abnormal) p 0 = 1 no singular arc (Bonnard - Caillau - Faubourg - Gergaud - Haberkorn - Noailles - Trélat)
7 Shooting method Find a zero of r(t f, p 0 ) r f ξ Kf (p 0 ) 1 v(t f, p 0 ) v f ξ ef (p S(t f, p 0 ) = B γ(t f, p 0 ) γ f p ϕ(t f, p 0 ) A or 0 ) B p ϕ(t f, p 0 ) A, p m(t f, p 0 ) 1 p m(t f, p 0 ) 1 Main problem: how to make the shooting method converge? initialization of the shooting method discontinuities of the optimal control
8 Shooting method Main problem: how to make the shooting method converge? initialization of the shooting method discontinuities of the optimal control Several methods: use first a direct method to provide a good initial guess, e.g. AMPL combined with IPOPT: R. Fourer, D.M. Gay, B.W. Kernighan, AMPL: A modeling language for mathematical programming, Duxbury Press, Brooks-Cole Publishing Company (1993). A. Wächter, L.T. Biegler On the implementation of an interior-point lter line- search algorithm for large-scale nonlinear programming, Mathematical Programming 106 (2006), but usual flaws of direct methods (computationally demanding, lack of numerical precision).
9 Shooting method Main problem: how to make the shooting method converge? initialization of the shooting method discontinuities of the optimal control Several methods: use the impulse transfer solution to provide a good initial guess: P. Augros, R. Delage, L. Perrot, Computation of optimal coplanar orbit transfers, AIAA but valid only for nearly circular initial and final orbits. See also: J. Gergaud, T. Haberkorn, Orbital transfer: some links between the low-thrust and the impulse cases, Acta Astronautica 60, no. 6-9 (2007), L.W. Neustadt, A general theory of minimum-fuel space trajectories, SIAM Journal on Control 3, no. 2 (1965),
10 Shooting method Main problem: how to make the shooting method converge? initialization of the shooting method discontinuities of the optimal control Several methods: multiple shooting method parameterized by the number of thrust arcs: H. J. Oberle, K. Taubert, Existence and multiple solutions of the minimum-fuel orbit transfer problem, J. Optim. Theory Appl. 95 (1997),
11 Shooting method Main problem: how to make the shooting method converge? initialization of the shooting method discontinuities of the optimal control Several methods: differential or simplicial continuation method linking the minimization of the L 2 -norm of the control to the minimization of the fuel consumption: J. Gergaud, T. Haberkorn, P. Martinon, Low thrust minimum fuel orbital transfer: an homotopic approach, J. Guidance Cont. Dyn. 27, 6 (2004), P. Martinon, J. Gergaud, Using switching detection and variational equations for the shooting method, Optimal Cont. Appl. Methods 28, no. 2 (2007), but not adapted for high-thrust transfer.
12 Flattening the Earth Observation: Solving the optimal control problem for a flat Earth model with constant gravity is simple and algorithmically very efficient. In view of that: Continuation from this simple model to the initial round Earth model.
13 Simplified flat Earth model System ẋ(t) = v x (t) ḣ(t) = v h (t) v x (t) = Tmax ux (t) m(t) v h (t) = Tmax m(t) u h(t) g 0 max m(t f ) t f free ṁ(t) = βt max qu x (t) 2 + u h (t) 2 Control Control (u x ( ), u h ( )) such that u x ( ) 2 + u h ( ) 2 1 initial conditions: x(0) = x 0, h(0) = h 0, v x (0) = v x0, v h (0) = v h0, m(0) = m 0 final conditions: h(t f ) = h f, v x (t f ) = v xf, v h (t f ) = 0
14 Modified flat Earth model Idea: mapping circular orbits to horizontal trajectories 8 >< >: x = r ϕ h = r r T v x = v cos γ v h = v sin γ 8 >< >: r = r T + h x ϕ = qr T +h v = vx 2 + vh 2 γ = arctan v h v x «««ux cos γ sin γ u1 = u h sin γ cos γ u 2 x f (v xf, v hf) (r T + h f)ϕ f x f v f γ f v h0 h f v 0 γ 0 h f v x0 h 0 h 0 r T ϕ f
15 Modified flat Earth model Plugging this change of coordinates into the initial round Earth model: leads to... ṙ(t) = v(t) sin γ(t) ϕ(t) = v(t) cos γ(t) r(t) v(t) = g(r(t)) sin γ(t) + Tmax m(t) u 1(t) v(t) γ(t) = r(t) g(r(t)) «cos γ(t) + Tmax v(t) m(t)v(t) u 2(t) ṁ(t) = βt max u(t)
16 Modified flat Earth model Modified flat Earth model x(t) ẋ(t) = v x (t) + v h (t) r T + h(t) ḣ(t) = v h (t) v x (t) = Tmax m(t) ux (t) vx (t)v h(t) r T + h(t) v h (t) = Tmax m(t) u vx (t)2 h(t) g(r T + h(t)) + r T + h(t) ṁ(t) = βt max u(t) Differences with the simplified flat Earth model (with constant gravity): the term in green: variable (usual) gravity. the terms in red: correcting terms allowing the existence of horizontal (periodic up to translation in x) trajectories with no thrust.
17 Continuation procedure Simplified flat Earth model (with constant gravity) continuation modified flat Earth model: procedure x(t) ẋ(t) = v x (t) + λ 2 v h (t) r T + h(t) ḣ(t) = v h (t) v x (t) = Tmax m(t) ux (t) λ v x (t)v h (t) 2 r T + h(t) v h (t) = Tmax m(t) u h(t) µ (r T +λ 1 h(t)) 2 + λ v x (t) 2 2 r T + h(t) 2 parameters: 0 λ λ 2 1 ṁ(t) = βt max qu x (t) 2 + u h (t) 2 λ 1 = λ 2 = 0: simplified flat Earth model with constant gravity λ 1 = 1, λ 2 = 0: simplified flat Earth model with usual gravity λ 1 = λ 2 = 1: modified flat Earth model (equivalent to usual round Earth)
18 Continuation procedure Simplified flat Earth model (with constant gravity) continuation modified flat Earth model: procedure x(t) ẋ(t) = v x (t) + λ 2 v h (t) r T + h(t) ḣ(t) = v h (t) v x (t) = Tmax m(t) ux (t) λ v x (t)v h (t) 2 r T + h(t) v h (t) = Tmax m(t) u h(t) µ (r T +λ 1 h(t)) 2 + λ v x (t) 2 2 r T + h(t) 2 parameters: 0 λ λ 2 1 ṁ(t) = βt max qu x (t) 2 + u h (t) 2 Two-parameters family of optimal control problems: (OCP) λ1,λ 2
19 Continuation procedure (OCP) 0,0 flat Earth model, constant gravity (linear) continuation on λ 1 [0, 1] final time t f free (OCP) 1,0 flat Earth model, usual gravity (linear) continuation on λ 2 [0, 1] final time t f fixed (OCP) 1,1 modified flat Earth model (equivalent to round Earth model) Application of the PMP to (OCP) λ1,λ 2 series of shooting problems.
20 Change of coordinates Remark: Once the continuation process has converged, we obtain the initial adjoint vector for (OCP) 1,1 in the modified coordinates. To recover the adjoint vector in the usual cylindrical coordinates, we use the general fact: Lemma Change of coordinates x 1 = φ(x) and u 1 = ψ(u) dynamics f 1 (x 1, u 1 ) = dφ(x).f (φ 1 (x 1 ), ψ 1 (u 1 )) and for the adjoint vectors: p 1 ( ) = t dφ(x( )) 1 p( ). Here, this yields: p r = x r T + h px + p h p ϕ = (r T + h)p x p v = cos γ p vx + sin γ p vh p γ = v( sin γ p vx + cos γ p vh ).
21 Analysis of the flat Earth model System ẋ = v x Initial conditions Final conditions ḣ = v h x(0) = x 0 x(t f ) free t f free v x = Tmax m ux v h = Tmax m u h g 0 ṁ = βt max qu 2 x + u 2 h h(0) = h 0 v x (0) = v x0 v h (0) = v h0 m(0) = m 0 h(t f ) = h f v x (t f ) = v xf v h (t f ) = 0 m(t f ) free max m(t f ) Theorem If h f > h 0 + v2 h0 2g 0, then the optimal trajectory is a succession of at most two arcs, and the thrust u( ) T max is either constant on [0, t f ] and equal to T max, or of the type T max 0, or of the type 0 T max.
22 Analysis of the flat Earth model Main ideas of the proof: Application of the PMP The switching function Φ = 1 m q p 2 v x + p 2 v h βp m satisfies: Φ = p h p vh q m pv 2 x + pv 2 h Φ = β u m Φ has at most one minimum Φ m qp 2 v x + p 2 vh Φ 2 + strategies T max, T max 0, 0 T max, or T max 0 T max The strategy T max 0 T max cannot occur ph 2 q m pv 2 x + pv 2 h
23 Shooting method in the flat Earth model A priori, we have: 5 unknowns p h, p vx, p vh (0), p m(0), and t f 5 equations h(t f ) = h f, v x (t f ) = v xf, v h (t f ) = 0, p m(t f ) = 1, H(t f ) = 0 but using several tricks and some system analysis, the shooting method can be simplified to: 3 unknowns p vx, p vh (0), and the first switching time t 1 3 equations h(t f ) = h f, v x (t f ) = v xf, v h (t 1 ) + g 0 t 1 = g 0 p vh (0) very easy and efficient (instantaneous) algorithm and the initialization of the shooting method is automatic (CV for any initial adjoint vector) automatic tool for initializing the continuation procedure
24 Numerical simulations T max = 180 kn Isp = 450 s λ λ Initial conditions ϕ 0 = 0 (SSO) h 0 = 200 km v 0 = 5.5 km/s γ 0 = 2 deg m 0 = kg λ p h p vx λ p vh p m Final conditions h f = 800 km v f = 7.5 km/s γ f = 0 deg (nearly circular final orbit) Evolution of the shooting function unknowns (p h, p vx, p vh, p m) (abscissa) with respect to homotopic parameter λ 2 (ordinate) continuous but not C 1 path: λ , λ 2 0.8, and λ : 0 λ : T max λ 2 0.8: T max 0 T max 0.8 λ : T max λ 2 1: T max 0 T max
25 Numerical simulations 2 x x 104 h (km) x (km) v x (km/s) v h (km/s) t (s) u x u h u m (kg) t (s) Trajectory and control strategy of (OCP) 1,0 (dashed) and (OCP) 1,1 (plain). t f 1483s Remark In the case of a final orbit (no injecting point): additional continuation on transversality conditions.
26 Numerical simulations Comparison with a direct method: Heun (RK2) discretization with N points combination of AMPL with IPOPT needs however a careful initial guess Continuation method 3 seconds: (OCP) 0,0 : instantaneous from (OCP) 0,0 to (OCP) 1,0 : 0.5 second from (OCP) 1,0 to (OCP) 1,1 : 2.5 seconds Accuracy: Direct method N = 100: 5 seconds N = 1000: 165 seconds Accuracy: 10 6
27 Conclusion Algorithmic procedure to solve the problem of minimization of fuel consumption for the coplanar orbit transfer problem by shooting method approach Does not require any careful initial guess Open questions Is this procedure systematically efficient, for any possible coplanar orbit transfer? Extension to 3D
Continuation from a flat to a round Earth model in the coplanar orbit transfer problem
Continuation from a flat to a round Earth model in the coplanar orbit transfer problem Max Cerf, Thomas Haberkorn, Emmanuel Trélat To cite this version: Max Cerf, Thomas Haberkorn, Emmanuel Trélat. Continuation
More informationLow Thrust Minimum-Fuel Orbital Transfer: A Homotopic Approach
Low Thrust Minimum-Fuel Orbital Transfer: A Homotopic Approach 3rd International Workshop on Astrodynamics Tools and Techniques ESA, DLR, CNES ESTEC, Noordwijk Joseph Gergaud and Thomas Haberkorn 2 5 October
More informationOptimal control and applications to aerospace problems
Optimal control and applications to aerospace problems E. Trélat Université Pierre et Marie Curie (Paris 6), Laboratoire J.-L. Lions and Institut Universitaire de France 6th European Congress of Mathematics,
More informationContinuation from a flat to a round Earth model in the coplanar orbit transfer problem
OPTIMAL CONTROL APPLICATIONS AND METHODS Optim. Control Appl. Meth. 212; 33:654 675 Published online 28 September 211 in Wiley Online Library (wileyonlinelibrary.com)..116 Continuation from a flat to a
More informationRemarks on Quadratic Hamiltonians in Spaceflight Mechanics
Remarks on Quadratic Hamiltonians in Spaceflight Mechanics Bernard Bonnard 1, Jean-Baptiste Caillau 2, and Romain Dujol 2 1 Institut de mathématiques de Bourgogne, CNRS, Dijon, France, bernard.bonnard@u-bourgogne.fr
More informationCOUPLED OPTIMIZATION OF LAUNCHER AND ALL-ELECTRIC SATELLITE TRAJECTORIES
COUPLED OPTIMIZATION OF LAUNCHER AND ALL-ELECTRIC SATELLITE TRAJECTORIES M. Verlet (1), B. Slama (1), S. Reynaud (1), and M. Cerf (1) (1) Airbus Defence and Space, 66 Route de Verneuil, 78133 Les Mureaux,
More informationSECOND ORDER OPTIMALITY CONDITIONS WITH APPLICATIONS. B. Bonnard. J.-B. Caillau
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 27 pp. 45 54 SECOND ORDER OPTIMALITY CONDITIONS WITH APPLICATIONS B. Bonnard Institut de Mathématiques (UMR CNRS 5584)
More informationLow thrust minimum-fuel orbital transfer: an homotopic approach
Low thrust minimum-fuel orbital transfer: an homotopic approach T. Haberkorn, P. Martinon, and J. Gergaud Abstract: We describe in this paper the study of an earth orbital transfer with a low thrust (typically
More informationAUTONOMOUS AND ROBUST RENDEZVOUS GUIDANCE ON ELLIPTICAL ORBIT SUBJECT TO J 2 PERTURBATION.
AUTONOMOUS AND ROBUST RENDEZVOUS GUIDANCE ON ELLIPTICAL ORBIT SUBJECT TO J 2 PERTURBATION Emmanuel GOGIBUS (1), Hervé CHARBONNEL (2), Patrick DELPY (3) (1) Astrium Space Transportation, 66 route de Verneuil,
More informationOptimal control theory and some applications to aerospace problems
Optimal control theory and some applications to aerospace problems Emmanuel Trélat Abstract. In this proceedings article we first shortly report on some classical techniques of nonlinear optimal control
More informationAnalysis of optimal strategies for soft landing on the Moon from lunar parking orbits
Analysis of optimal strategies for soft landing on the Moon from lunar parking orbits R V Ramanan and Madan Lal Aerospace Flight Dynamics Group, Vikram Sarabhai Space Centre, Thiruvananthapuram 695 022,
More informationLYAPUNOV-BASED ELLIPTICAL TO CIRCULAR PLANAR ORBIT TRANSFERS IN LEVI-CIVITA COORDINATES
(Preprint) AAS 2-66 LYAPUNOV-BASED ELLIPTICAL TO CIRCULAR PLANAR ORBIT TRANSFERS IN LEVI-CIVITA COORDINATES Sonia Hernandez, Maruthi R. Akella, and Cesar A. Ocampo INTRODUCTION We consider planar orbit
More informationHOMOTOPY METHOD FOR MINIMUM CONSUMPTION ORBIT TRANSFER PROBLEM
ESAIM: COCV April 26, Vol. 2, 294 3 DOI:.5/cocv:263 ESAIM: Control, Optimisation and Calculus of Variations www.edpsciences.org/cocv HOMOTOPY METHOD FOR MINIMUM CONSUMPTION ORBIT TRANSFER PROBLEM Joseph
More informationOptimization techniques for autonomous underwater vehicles: a practical point of view
Optimization techniques for autonomous underwater vehicles: a practical point of view M. Chyba, T. Haberkorn Department of Mathematics University of Hawaii, Honolulu, HI 96822 Email: mchyba@math.hawaii.edu,
More informationTheory and Applications of Constrained Optimal Control Proble
Theory and Applications of Constrained Optimal Control Problems with Delays PART 1 : Mixed Control State Constraints Helmut Maurer 1, Laurenz Göllmann 2 1 Institut für Numerische und Angewandte Mathematik,
More informationEXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING. Dilmurat Azimov
EXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING Dilmurat Azimov University of Hawaii at Manoa 254 Dole Street, Holmes 22A Phone: (88)-956-2863, E-mail: azimov@hawaii.edu
More informationSpacecraft Trajectory Optimization. Edited. Bruce A. Conway. Published by Cambridge University Press 2010
Spacecraft Trajectory Optimization Edited by Bruce A. Conway Published by Cambridge University Press 2010 Cambridge University Press 2010 ISBN 978-0-521-51850-5 2 Primer Vector Theory and Applications
More informationGoddard s Problem Steven Finch
Goddard s Problem Steven Finch June 12, 2015 Arocketliftsoff vertically at time =0. Let ( ) be the mass of the rocket (payload and fuel) and ( ) be the altitude. We wish to choose the thrust ( ) and a
More informationAsteroid Rendezvous Missions
Asteroid Rendezvous Missions An Application of Optimal Control G. Patterson, G. Picot, and S. Brelsford University of Hawai`i at Manoa 25 June, 215 1/5 G. Patterson, G. Picot, and S. Brelsford Asteroid
More informationOPTIMIZING PERIAPSIS-RAISE MANEUVERS USING LOW-THRUST PROPULSION
AAS 8-298 OPTIMIZING PERIAPSIS-RAISE MANEUVERS USING LOW-THRUST PROPULSION Brenton J. Duffy and David F. Chichka This study considers the optimal control problem of maximizing the raise in the periapsis
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationMinimum Fuel Round Trip from a L 2 Earth-Moon Halo Orbit to Asteroid 2006 RH 120
Minimum Fuel Round Trip from a L Earth-Moon Halo Orbit to Asteroid 6 Monique Chyba, Thomas Haberkorn and Robert Jedicke Abstract The goal of this paper is to design a spacecraft round trip transfer from
More informationLimit solutions for control systems
Limit solutions for control systems M. Soledad Aronna Escola de Matemática Aplicada, FGV-Rio Oktobermat XV, 19 e 20 de Outubro de 2017, PUC-Rio 1 Introduction & motivation 2 Limit solutions 3 Commutative
More informationTutorial on Control and State Constrained Optimal Control Problems
Tutorial on Control and State Constrained Optimal Control Problems To cite this version:. blems. SADCO Summer School 211 - Optimal Control, Sep 211, London, United Kingdom. HAL Id: inria-629518
More information: low-thrust transfer software, optimal control problem, averaging techniques.
J. Fourcade S. Geffroy R.Epenoy Centre National d Etudes Spatiales 8 avenue Edouard Belin 4 Toulouse cedex 4 France e-mail : Jean.Fourcade@cnes.fr Sophie.Geffroy@cnes.fr Richard.Epenoy@cnes.fr Low thrust
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationNumerical Methods for Optimal Control Problems. Part I: Hamilton-Jacobi-Bellman Equations and Pontryagin Minimum Principle
Numerical Methods for Optimal Control Problems. Part I: Hamilton-Jacobi-Bellman Equations and Pontryagin Minimum Principle Ph.D. course in OPTIMAL CONTROL Emiliano Cristiani (IAC CNR) e.cristiani@iac.cnr.it
More informationOptimal control and applications to aerospace: some results and challenges
Optimal control and applications to aerospace: some results and challenges E. Trélat Abstract This article surveys the classical techniques of nonlinear optimal control such as the Pontryagin Maximum Principle
More informationSOLVING NONLINEAR OPTIMAL CONTROL PROBLEMS WITH STATE AND CONTROL DELAYS BY SHOOTING METHODS COMBINED WITH NUMERICAL CONTINUATION ON THE DELAYS
1 2 3 4 5 6 7 8 9 1 11 12 13 14 15 16 17 18 19 2 21 22 23 24 25 26 27 28 29 3 31 32 33 34 35 36 SOLVING NONLINEAR OPTIMAL CONTROL PROBLEMS WITH STATE AND CONTROL DELAYS BY SHOOTING METHODS COMBINED WITH
More informationS12 PHY321: Practice Final
S12 PHY321: Practice Final Contextual information Damped harmonic oscillator equation: ẍ + 2βẋ + ω0x 2 = 0 ( ) ( General solution: x(t) = e [A βt 1 exp β2 ω0t 2 + A 2 exp )] β 2 ω0t 2 Driven harmonic oscillator
More informationMinimum fuel round trip from a L_2 Earth-Moon Halo orbit to Asteroid 2006 RH_120
Minimum fuel round trip from a L_ Earth-Moon Halo orbit to Asteroid 6 RH_ Monique Chyba, Thomas Haberkorn, Robert Jedicke To cite this version: Monique Chyba, Thomas Haberkorn, Robert Jedicke. Minimum
More informationSecond Order Analysis of Control-affine Problems with Control and State Constraints
& state constraint Second Order Analysis of Control-affine Problems with Control and State Constraints María Soledad Aronna Escola de Matemática Aplicada, Fundação Getulio Vargas, Rio de Janeiro, Brazil
More informationAscent Phase Trajectory Optimization for a Hypersonic Vehicle Using Nonlinear Programming
Ascent Phase Trajectory Optimization for a Hypersonic Vehicle Using Nonlinear Programming H.M. Prasanna,D.Ghose, M.S. Bhat, C. Bhattacharyya, and J. Umakant 3 Department of Aerospace Engineering Department
More informationResults found by the CNES team (team #4)
3 rd Global Trajectory Optimisation Competition (GTOC3) organized by the Aerospace Propulsion Group of the Dipartimento di Energetica at Politecnico di Torino Results found by the CNES team (team #4) Presented
More informationSingular arcs in the generalized Goddard s Problem
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE arxiv:math/739v [math.oc] 3 Mar 27 Singular arcs in the generalized Goddard s Problem apport de recherche Frédéric Bonnans Pierre Martinon
More informationwhere s is the horizontal range, in this case 50 yards. Note that as long as the initial speed of the bullet is great enough to let it hit a target at
1 PHYS 31 Homework Assignment Due Friday, 13 September 00 1. A monkey is hanging from a tree limbat height h above the ground. A hunter 50 yards away from the base of the tree sees him, raises his gun
More informationL 1 -minimization for mechanical systems
L 1 -minimization for mechanical systems Z. Chen J.-B. Caillau Y. Chitour November 215 Abstract Second order systems whose drift is defined by the gradient of a given potential are considered, and minimization
More informationOptimal Control - Homework Exercise 3
Optimal Control - Homework Exercise 3 December 17, 2010 In this exercise two different problems will be considered, first the so called Zermelo problem the problem is to steer a boat in streaming water,
More informationOptimization problems for controlled mechanical systems: Bridging the gap between theory and application
0 Optimization problems for controlled mechanical systems: Bridging the gap between theory and application M. Chyba a, T. Haberkorn b, R.N. Smith c a University of Hawaii, Honolulu, HI 96822 USA Mathematics,
More informationSOLUTION OF OPTIMAL CONTINUOUS LOW-THRUST TRANSFER USING LIE TRANSFORMS
AAS 3-93 SOLUTION OF OPTIMAL CONTINUOUS LOW-THRUST TRANSFER USING LIE TRANSFORMS M. Sanjurjo-Rivo, * D. J. Scheeres, M. Lara and J. Peláez This paper addresses the problem of optimal constant continuous
More informationOn the Maximal Orbit Transfer Problem Marian Mureşan
The Mathematica Journal On the Maximal Orbit Transfer Problem Marian Mureşan Assume that a spacecraft is in a circular orbit and consider the problem of finding the largest possible circular orbit to which
More informationRendezvous Missions to Temporarily-Captured Near Earth Asteroids
Rendezvous Missions to Temporarily-Captured Near Earth Asteroids S Brelsford, Monique Chyba, T Haberkorn, G Patterson To cite this version: S Brelsford, Monique Chyba, T Haberkorn, G Patterson. Rendezvous
More informationAveraging and optimal control of elliptic Keplerian orbits with low propulsion
Systems & Control Letters 55 (6) 755 76 www.elsevier.com/locate/sysconle Averaging and optimal control of elliptic Keplerian orbits with low propulsion Bernard Bonnard a,, Jean-Baptiste Caillau b, Romain
More informationMinimum Time Ascent Phase Trajectory Optimization using Steepest Descent Method
IJCTA, 9(39), 2016, pp. 71-76 International Science Press Closed Loop Control of Soft Switched Forward Converter Using Intelligent Controller 71 Minimum Time Ascent Phase Trajectory Optimization using
More informationBernard Bonnard 1, Jean-Baptiste Caillau 2 and Emmanuel Trélat 3
ESAIM: COCV Vol. 13, N o 2, 27, pp. 27 236 DOI: 1.151/cocv:2712 ESAIM: Control, Optimisation and Calculus of Variations www.edpsciences.org/cocv SECOND ORDER OPTIMALITY CONDITIONS IN THE SMOOTH CASE AND
More informationISTS-2017-d-114 ISSFD
ISTS-27-d-4 ISSFD-27-4 A Minimum-Fuel Fixed-Time Low-Thrust Rendezvous Solved with the Switching Systems Theory By Clément GAZZINO, ) Denis ARZELIER, ) Luca CERRI, 2) Damiana LOSA, 3) Christophe LOUEMBET,
More informationOrbit Transfer Optimization for Multiple Asteroid Flybys
Orbit Transfer Optimization for Multiple Asteroid Flybys Bruno Victorino Sarli 1 and Yasuhiro Kawakatsu 2 1 Department of Space and Astronautical Science, The Graduate University for Advance Studies, Sagamihara,
More informationSub-Riemannian geometry in groups of diffeomorphisms and shape spaces
Sub-Riemannian geometry in groups of diffeomorphisms and shape spaces Sylvain Arguillère, Emmanuel Trélat (Paris 6), Alain Trouvé (ENS Cachan), Laurent May 2013 Plan 1 Sub-Riemannian geometry 2 Right-invariant
More informationOrbital Mechanics. John E. Prussing & Bruce A.Conway. Published by Oxford University Press 2013
Orbital Mechanics by John E. Prussing & Bruce A.Conway Published by Oxford University Press 2013 Copyright 2013, 1993 by Oxford University Press ISBN: 978-0-19-983770-0 Chapter 8 Continuous-Thrust Orbit
More informationarxiv: v1 [math.oc] 24 Sep 2015
Coupled attitude and trajectory optimization for a launcher tilting maneuver arxiv:159.7519v1 [math.oc] 24 Sep 215 Jiamin Zhu Emmanuel Trélat Max Cerf Abstract We study the minimum time control problem
More informationOptimal Configuration of Tetrahedral Spacecraft Formations 1
The Journal of the Astronautical Sciences, Vol. 55, No 2, April June 2007, pp. 141 169 Optimal Configuration of Tetrahedral Spacecraft Formations 1 Geoffrey T. Huntington, 2 David Benson, 3 and Anil V.
More informationarxiv: v2 [math.oc] 25 Jul 2013
A SHOOTING ALGORITHM FOR OPTIMAL CONTROL PROBLEMS WITH SINGULAR ARCS 1,2 M. SOLEDAD ARONNA, J. FRÉDÉRIC BONNANS, AND PIERRE MARTINON arxiv:1206.0839v2 [math.oc] 25 Jul 2013 Abstract. In this article, we
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationarxiv: v1 [math.oc] 15 Dec 2015
Pontryagin maximum principle for optimal sampled-data control problems. Loïc Bourdin Emmanuel Trélat arxiv:52.4797v [math.oc] 5 Dec 25 Université de Limoges, Institut de recherche XLIM, CNRS UMR 7252,
More informationLow-Thrust Trajectories to the Moon
3rd WSEAS International Conference on APPLIED and THEORETICAL MECHANICS, Spain, December 14-16, 7 257 Low-Thrust Trajectories to the Moon ANTONIO F. B. A. PRADO Space Mechanics and Control Division INPE
More informationTime-minimal orbital transfers to temporarily-captured natural Earth satellites
Time-minimal orbital transfers to temporarily-captured natural Earth satellites Monique Chyba, Mikael Granvik,Robert Jedicke,Geoff Patterson,Gautier Picot, Jeremie Vaubaillon Abstract In this paper we
More informationLow-Thrust Trajectory Optimization with No Initial Guess
Low-Thrust Trajectory Optimization with No Initial Guess By Nathan L. Parrish 1) and Daniel J. Scheeres 1) 1) Colorado Center for Astrodynamics Research, University of Colorado, Boulder, USA (Received
More informationOptimal Trajectory Control of an Occulter-Based Planet-Finding Telescope
Optimal Trajectory Control of an Occulter-Based Planet-Finding Telescope Egemen Kolemen and N. Jeremy Kasdin Princeton University, Princeton, NJ, 08544, USA The optimal configuration of a satellite formation
More informationHAPPY BIRTHDAY TONY! Chicago Blochfest June 30, 2105
HAPPY BIRTHDAY TONY! 1 CLEBSCH OPTIMAL CONTROL Tudor S. Ratiu Section de Mathématiques Ecole Polytechnique Fédérale de Lausanne, Switzerland tudor.ratiu@epfl.ch 2 PLAN OF THE PRESENTATION Symmetric representation
More informationAn application of PL continuation methods to singular arcs problems
An application of PL continuation methods to singular arcs problems Pierre Martinon and Joseph Gergaud ENSEEIHT-IRIT, UMR CNRS 555, 2 rue Camichel; F-37 Toulouse martinon,gergaud@enseeiht.fr Summary. Among
More informationSingular boundary value problems
Department of Mathematics, Faculty of Science, Palacký University, 17. listopadu 12, 771 46 Olomouc, Czech Republic, e-mail: irena.rachunkova@upol.cz Malá Morávka, May 2012 Overview of our common research
More informationFeedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field
Vol. 4, No. 4, 23 Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Ashraf H. Owis Department of Astronomy, Space and Meteorology, Faculty of Science, Cairo University Department
More informationMotion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space
Motion in Space MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background Suppose the position vector of a moving object is given by r(t) = f (t), g(t), h(t), Background
More informationExtremal Trajectories for Bounded Velocity Mobile Robots
Extremal Trajectories for Bounded Velocity Mobile Robots Devin J. Balkcom and Matthew T. Mason Abstract Previous work [3, 6, 9, 8, 7, 1] has presented the time optimal trajectories for three classes of
More informationThe optimization of electric propulsion (EP) trajectories for interplanetary missions is a quite difficult
Indirect Optimization Method for Low-Thrust Interplanetary Trajectories IEPC-27-356 Presented at the 3 th International Electric Propulsion Conference, Florence, Italy Lorenzo Casalino, Guido Colasurdo
More information7.2 Maximization of the Range of a Rocket
138 CHAPTER 7. SOME APPLICATIONS The counterintuitive answer that a supersonic aircraft must dive first in order to climb to a iven altitude in minimum time was first discovered by Walter Denham and Art
More informationPhysics 170 Week 9 Lecture 2
Physics 170 Week 9 Lecture 2 http://www.phas.ubc.ca/ gordonws/170 Physics 170 Week 9 Lecture 2 1 Textbook Chapter 1: Section 1.6 Physics 170 Week 9 Lecture 2 2 Learning Goals: We will solve an example
More informationThe Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 2 Dianne P. O Leary c 2008 The Plan
More informationPeriodic solutions of weakly coupled superlinear systems
Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize
More information3 Space curvilinear motion, motion in non-inertial frames
3 Space curvilinear motion, motion in non-inertial frames 3.1 In-class problem A rocket of initial mass m i is fired vertically up from earth and accelerates until its fuel is exhausted. The residual mass
More informationClassical Numerical Methods to Solve Optimal Control Problems
Lecture 26 Classical Numerical Methods to Solve Optimal Control Problems Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Necessary Conditions
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More informationGravity Assisted Maneuvers for Asteroids using Solar Electric Propulsion
Gravity Assisted Maneuvers for Asteroids using Solar Electric Propulsion Denilson P. S. dos Santos, Antônio F. Bertachini de A. Prado, Division of Space Mechanics and Control INPE C.P. 515, 17-310 São
More informationLesson 7. Luis Anchordoqui. Physics 168. Tuesday, October 10, 17
Lesson 7 Physics 168 1 Eruption of a large volcano on Jupiter s moon When volcano erupts speed of effluence exceeds escape speed of Io and so a stream of particles is projected into space Material in stream
More informationSolving Optimal Continuous Thrust Rendezvous Problems with Generating Functions
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 29, No. 2, March April 2006 Solving Optimal Continuous Thrust Rendezvous Problems with Generating Functions Chandeok Park, Vincent Guibout, and Daniel J.
More informationOrbital Transfer Trajectory Optimization. James K Whiting
Orbital Transfer Trajectory Optimization by James K Whiting Submitted to the Department of Aeronautical and Astronautical Engineering in partial fulfillment of the requirements for the degree of Master
More informationOptimal Control based Time Optimal Low Thrust Orbit Raising
Optimal Control based Time Optimal Low Thrust Orbit Raising Deepak Gaur 1, M. S. Prasad 2 1 M. Tech. (Avionics), Amity Institute of Space Science and Technology, Amity University, Noida, U.P., India 2
More informationSolving Optimal Continuous Thrust Rendezvous Problems with Generating Functions
AIAA Guidance, Navigation, and Control Conference and Ehibit 15-18 August 25, San Francisco, California AIAA 25-677 Solving Optimal Continuous Thrust Rendezvous Problems with Generating Functions Chandeok
More informationHalo Orbit Mission Correction Maneuvers Using Optimal Control
Halo Orbit Mission Correction Maneuvers Using Optimal Control Radu Serban and Linda Petzold (UCSB) Wang Koon, Jerrold Marsden and Shane Ross (Caltech) Martin Lo and Roby Wilson (JPL) Wang Sang Koon Control
More informationShoot-1.1 Package - User Guide
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE Shoot-. Package - User Guide Pierre Martinon Joseph Gergaud inria-287879, version 2-2 Sep 28 N 354 Juin 28 Thème NUM apport technique ISSN
More informationCalculus of Variations Summer Term 2015
Calculus of Variations Summer Term 2015 Lecture 14 Universität des Saarlandes 24. Juni 2015 c Daria Apushkinskaya (UdS) Calculus of variations lecture 14 24. Juni 2015 1 / 20 Purpose of Lesson Purpose
More informationAbout the minimal time crisis problem and applications
About the minimal time crisis problem and applications Térence Bayen, Alain Rapaport To cite this version: Térence Bayen, Alain Rapaport. About the minimal time crisis problem and applications. IFAC Word
More informationSUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS. Sandro da Silva Fernandes. Cleverson Maranhão Porto Marinho
SUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS Sandro da Silva Fernandes Instituto Tecnológico de Aeronáutica, São José dos Campos - 12228-900 - SP-Brazil, (+55) (12) 3947-5953 sandro@ita.br Cleverson
More informationproblem score possible Total 60
Math 32 A: Midterm 1, Oct. 24, 2008 Name: ID Number: Section: problem score possible 1. 10 2. 10 3. 10 4. 10 5. 10 6. 10 Total 60 1 1. Determine whether the following points are coplanar: a. P = (8, 14,
More informationDirect and indirect methods for optimal control problems and applications in engineering
Direct and indirect methods for optimal control problems and applications in engineering Matthias Gerdts Computational Optimisation Group School of Mathematics The University of Birmingham gerdtsm@maths.bham.ac.uk
More informationBasic Ascent Performance Analyses
Basic Ascent Performance Analyses Ascent Mission Requirements Ideal Burnout Solution Constant & Average Gravity Models Gravity Loss Concept Effect of Drag on Ascent Performance Drag Profile Approximation
More informationOptimal Trajectories of Anti-satellite Missiles
Licensed under the Creative Commons Attribution Noncommercial No Derivatives Communications to SIMAI Congress, ISSN 1827-9015, Vol. 3 (2009) 265 (12pp) DOI: 10.1685/CSC09265 Optimal Trajectories of Anti-satellite
More informationIntroduction to the Optimal Control Software GPOPS II
Introduction to the Optimal Control Software GPOPS II Anil V. Rao Department of Mechanical and Aerospace Engineering University of Florida Gainesville, FL 32611-625 Tutorial on GPOPS II NSF CBMS Workshop
More informationOptimal Control Theory SF 2852
Optimal Control Theory SF 2852 Johan Karlsson Optimization and Systems Theory, Department of Mathematics Royal Institute of Technology (KTH) Spring 2017 Optimal Control Theory SF 2852 1 Course information
More informationTime-optimal control of a 3-level quantum system and its generalization to an n-level system
Proceedings of the 7 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 7 Time-optimal control of a 3-level quantum system and its generalization to an n-level
More informationMinimum-Time Trajectory Optimization of Low-Thrust Earth-Orbit Transfers with Eclipsing
Minimum-Time Trajectory Optimization of Low-Thrust Earth-Orbit Transfers with Eclipsing Kathryn F. Graham Anil V. Rao University of Florida Gainesville, FL 32611-625 Abstract The problem of determining
More informationProject Optimal Control: Optimal Control with State Variable Inequality Constraints (SVIC)
Project Optimal Control: Optimal Control with State Variable Inequality Constraints (SVIC) Mahdi Ghazaei, Meike Stemmann Automatic Control LTH Lund University Problem Formulation Find the maximum of the
More informationOrbital Motion in Schwarzschild Geometry
Physics 4 Lecture 29 Orbital Motion in Schwarzschild Geometry Lecture 29 Physics 4 Classical Mechanics II November 9th, 2007 We have seen, through the study of the weak field solutions of Einstein s equation
More informationDirect Methods. Moritz Diehl. Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U.
Direct Methods Moritz Diehl Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U. Leuven Belgium Overview Direct Single Shooting Direct Collocation Direct Multiple
More informationV Requirements for a Gun Assisted Launch to Circular Orbit
V Requirements for a Gun Assisted Launch to Circular Orbit Gerry Flanagan The Alna Space Program May 12, 2011 Introduction and Assumptions An earth-based gun can be used to send a projectile into space,
More informationConstrained shape spaces of infinite dimension
Constrained shape spaces of infinite dimension Sylvain Arguillère, Emmanuel Trélat (Paris 6), Alain Trouvé (ENS Cachan), Laurent May 2013 Goal: Adding constraints to shapes in order to better fit the modeled
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.it.edu 16.323 Principles of Optial Control Spring 2008 For inforation about citing these aterials or our Ters of Use, visit: http://ocw.it.edu/ters. 16.323 Lecture 10 Singular
More informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18
More informationOn the optimality of Dubins paths across heterogeneous terrain
On the optimality of Dubins paths across heterogeneous terrain Ricardo G. Sanfelice and Emilio Frazzoli {sricardo,frazzoli}@mit.edu Laboratory for Information and Decision Systems Massachusetts Institute
More informationDifferential continuation for regular optimal control problems
Optimization Methods & Software Vol. 27, No. 2, April 212, 177 196 Differential continuation for regular optimal control problems Downloaded by [Bibl de L'Univ de Bourgogne], [J.-B. Caillau] at 5:57 4
More information