Upgrade to Pro — share decks privately, control downloads, hide ads and more …

Manœuvre planaire d'un vaisseau spatial - Le ph...

Avatar for GdR MOA 2015 GdR MOA 2015
December 03, 2015

Manœuvre planaire d'un vaisseau spatial - Le phénomène de chattering et la solution sous-optimale

by J. Zhu

Avatar for GdR MOA 2015

GdR MOA 2015

December 03, 2015
Tweet

More Decks by GdR MOA 2015

Other Decks in Science

Transcript

  1. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Manœuvre planaire d’un vaisseau spatial - Ph´ enom` ene de chattering et solution sous-optimale Jiamin ZHU Laboratoire Jacques-Louis Lions Collaboration: E. TR´ ELAT, M. CERF [email protected] 3/12/2015
  2. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Table of content 1 Formulation du probl` eme MTTP 2 Analysis des ext´ emales 3 Solution sous-optimale 4 Conclusions
  3. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Mod` ele du probl` eme • ´ Etat x ∈ R4: x = (x1, x2, x3, x4) := (vx , vy , θ, ω) • Contrˆ ole u ∈ R: u = −µ/µmax ∈ [−1, 1] • Coordonn´ ees de r´ ef´ erence: x0 ∧ y0 • Coordonn´ ees de corps: ea ∧ ec
  4. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Mod` ele du probl` eme • Syst` eme ` a contrˆ ole affine ˙ x = f0(x) + uf1(x) • Champs de vecteurs lisses f0 =(a cos x3 − cx1x2) ∂ ∂x1 + (a sin x3 + cx2 1 − g0) ∂ ∂x2 + (x4 − cx1) ∂ ∂x3 , f1 =b ∂ ∂x4 . o` u a = T/m, b = Tlµmax /I, c = 1/r sont des constantes positives.
  5. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Formulation du probl` eme MTTP Amener le syst` eme ˙ x = f0(x) + uf1(x) de x(0) = x0 (condition initiale) ` a x(tf ) ∈ M1 (condition finale) en temps minimal tf sous la contrainte de contrˆ ole u(t) ∈ [−1, 1], o` u M1 ={(x1, x2, x3, x4) ∈ R4| x2(tf ) − x1(tf ) tan x3f = 0, x3f = x3(tf ), x4f = x4(tf )}
  6. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Table of content 1 Formulation du probl` eme MTTP 2 Analysis des ext´ emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering 3 Solution sous-optimale 4 Conclusions
  7. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Application de Principe du Maximum de Pontryagin • Vecteur adjoint p = (p1, p2, p3, p4), r´ eel p0 ≤ 0 • ´ Equation des extr´ emales ˙ x = ∂H ∂p , ˙ p = − ∂H ∂x ⇒ (x(·), p(·), p0, u(·)) une extr´ emale • Hamiltonien: H := H0 + uH1 = p, f0(x) + u p, f1(x) • Condition de maximisation: H(x(t), p(t), p0, u(t)) = max −1≤v(t)≤1 H(x(t), p(t), p0, v(t)) • Condition de transversalit´ e (p1(tf ), p2(tf )) ⊥ ea(tf ) ⇒ p1(tf ) = −p2(tf ) tan x3f
  8. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Extr´ emales r´ eguli` ere et singuli` ere • Hamiltonien: H := H0 + uH1 = p, f0(x) + u p, f1(x) • Fonction de commutation (switching function): H1(t) = p(t), f1(x(t)) • Condition de maximisation: H(x(t), p(t), p0, u(t)) = max −1≤v(t)≤1 H(x(t), p(t), p0, v(t)) • Contrˆ ole r´ egulier (Bang-Bang): u(t) = signH1(t), H1(t) = 0
  9. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Extr´ emales r´ eguli` ere et singuli` ere • Fonction de commutation(switching function): H1(t) = p(t), f1(x(t)) • Contrˆ ole singulier: H1(t) ≡ 0, t ∈ I ⊂ [0, tf ] • Ordre du contrˆ ole singulier: un contrˆ ole singulier d’ordre k apparaˆ ıt lin´ eairement dans la d´ eriv´ ee d’ordre 2k de la fonction H1 • Condition g´ en´ eralis´ ee de Legendre-Clebsch (condition n´ ecessaire d’optimalit´ e d’ordre k): (−1)k ∂ ∂u d2kH1 dt2k ≤ 0, • Contrˆ ole singulier d’ordre un nombre pair ⇒ Chattering (H.J. Kelley, R.E. Kopp, H.G. Moyer, ect.)
  10. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Ph´ enom` ene de chattering (Ph´ enom` ene de Fuller) • Jonction singuli` ere: jonction d’un arc singulier et d’un arc r´ egulier • Contrˆ ole chattering: contrˆ ole commute infiniment sur un intervalle compact pr` es d’une jonction singuli` ere Figure : Une illustration du ph´ enom` ene de chattering.
  11. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Jonction singuli` ere et chattering • Corollaire (jonction): For the problem (MTTP), any optimal singular arc can- not be connected with a nontrivial bang arc. We must then have chattering, in the following sense. Let u be an optimal control, solution of (MTTP), and assume that u is singular on the sub-interval (t1, t2) ⊂ [0, tf ] and is non-singular else- where. If t1 > 0 (resp., if t2 < tf ) then, for every > 0, the control u switchings an infinite number of times over the time interval [t1 − , t1] (resp., on [t2, t2 + ]). • Preuve par l’absurde
  12. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Portrait des solutions proches de la surface singuli` ere • Surface singuli` ere S := {(x, p)|H1 ≡ 0} • Corollaire (portrait des solutions) For the problem (MTTP), there exist two subbundles Σ+ and Σ− having the singular surface S as a base, and two fibers N+ and N− of dimension two filled by chattering solutions. • Preuve: “Theorem on Bundles with Chattering Fibers” (M.I.Zelikin, V.F. Borisov)
  13. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Optimalit´ e de chattering - M´ ethode de M.I. Zelikin et V.F. Borisov • Changement de variable: (x, p) ∈ R2n → (z, w) ∈ R2n ⇒ syst` eme semi-canonique • Construire un sous vari´ et´ e Lagrangienne (i.e. form´ e par les extr´ emales optimales) • Condition suffisante d’optimalit´ e: La projection canonique π∗ : (z, w) → x restreinte ` a ces extr´ emales est r´ eguli` ere - Identique ` a la th´ eorie des temps conjugu´ es
  14. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Optimalit´ e de Chattering - Syst` eme semi-canonique • Supposer que le contˆ ole singulier est d’ordre 2 • z ∈ R4: zi = H(i−1) 1 + o(H(i−1) 1 ), i = 1, · · · , 4 t.q. ˙ zi = zi+1+fi (z, w), i = 1, 2, 3, ˙ z4 = α(w)+uβ(w)+f4(z, w) o` u les fi sont infinit´ esimaux de zi+1, i = 1, 2, 3, 4 • Compl´ eter avec w ∈ R2n−4: ˙ w = F(z, w, u) • (z, w) = (z(x, p), w(x, p)) ∈ R2n t.q. det( D(z, w) D(x, p) ) = 0 • Surface singuli` ere S = {(x, p)|H1 ≡ 0} ⇒ S = {(z, w)|z = 0}
  15. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Optimalit´ e de chattering - Syst` eme semi-canonique - Nouvelle variable (z, w) du probl` eme MTTP:                z1 = p4, z2 = p(1) 4 , z3 = p(2) 4 , z4 = p(3) 4 + ac sin x3p3 , w1 = a(p1 sin x3 + p2 cos x3) , w2 = a sin x3(−x4 + 3cx1)p2 − a cos x3(−x4 + cx1 + cx2)p1 , w3 = p2/p1 , w4 = x1 + x2 . - Consid´ erons la domaine o` u la condition g´ en´ eralis´ ee de Legendre-Clebsch est v´ erifi´ ee (−1)k ∂ ∂u d2kH1 dt2k ≤ 0, ⇒ w1/w3 ≥ 0
  16. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Optimalit´ e de chattering - Sous vari´ et´ e Lagrangien • Surface singuli` ere: S = {(z, w)|z = 0} • S0 := S ∩ {H = 0} • Projection canonique: π∗ : (z, w) → x • M∗ 1 := {x|x2 = Ψ1(x1)} ∩ π∗(S0) • Sous vari´ et´ ed N: associer un point x ∈ M∗ 1 ` a un point (x, p(x)) ∈ S0 et satisfaire la condition de transversali´ e (dimension 1) ⇒ N est en co-dimension 6 + 1 et est un lift de M∗ 1 . • Ψ1(·) peut ˆ etre choisie t.q. N est transversal aux extr´ emales singuli` eres dans S.
  17. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Optimalit´ e de chattering - Sous vari´ et´ e Lagrangien • S0 := S ∩ {H = 0} • Projection canonique: π∗ : (z, w) → x • M∗ 1 := {x|x2 = Ψ1(x1)} ∩ π∗(S0) • Sous vari´ et´ e N: un lift de M∗ 1 • Sous vari´ et´ e N∗: constitu´ e des extr´ emales singuli` eres situ´ ees transverses ` a N ⇒ N∗ est une sous vari´ et´ e Lagrangienne
  18. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Optimalit´ e de chattering - Condition suffisante • Surface de commutation Γ = {(z, w)|z1 = 0} • Notons V = (v1, v2, v3, vk) • v1 vecteur tangent ` a N • v2 vitesse des extr´ emales singuli` eres en un point dans S • v3 vitesse des extr´ emales r´ eguli` eres en un point dans S • vk d´ eriv´ ee de la projection π∗ du vecteur ∂/∂z3 : dπ∗(∂/∂z3 ) • Th´ eor` eme (optimalit´ e local du chattering) Les solutions chattering du MTTP sont optimales locale- ment pour la C0-topologie si det V = 0 • Preuve: les conditions suffisantes pour la r´ egularit´ e de π∗ propos´ ees par M.I. Zelikin et V.F. Borisov
  19. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Table of content 1 Formulation du probl` eme MTTP 2 Analysis des ext´ emales 3 Solution sous-optimale Pr´ ediction Solutions sous-optimales 4 Conclusions
  20. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Pr´ ediction • Structure de la solution optimale (dynamique rapide-lente) • Supposons: • x30 ∈ [−π/2, π/2], t ∈ [0, tf ] • x3f ∈ [−π/2, π/2], (x3 (t) − x3f ) ∈ [−π/2, π/2], t ∈ [0, tf ] • manœuvre “counter-clockwise”: x3 (t) ∈ [x30, x3f + π/2] • manœuvre “clockwise”: x3 (t) ∈ [x3f − π/2, x30 ] • Manœuvre “counter-clockwise”/“clockwise”: x2f /x1f ≥ x20/x10 / x2f /x1f < x20/x10 • Rappel: (x1, x2) = (vx , vy )
  21. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Pr´ ediction • Soit t2 − t1 = 0 1er arc r´ egulier associ´ e ` a u = ±1 ∓1 ⇒ (x3(0), x4(0)) = (x30, x40), x3(t1) = x∗ 3 , x4(t1) = x∗ 4 ⇒ ∆1(x1, x2) 2nd arc r´ egulier associ´ e ` a u = ∓1 ±1 ⇒ (x3(t2), x4(t2)) = (x∗ 3 , x∗ 4 ), x3(tf ) = x3f , x4(tf ) = x4f ⇒ ∆2(x1, x2) • ⇒ (x1(tf ), x2(tf )) = (x10, x20) + ∆1(x1, x2) + ∆2(x1, x2) • Rappel: condition initiale x(0) = x0 condition finale x(tf ) ∈ M1 o` u M1 ={(x1, x2, x3, x4) ∈ R4| x2(tf ) − x1(tf ) tan x3f = 0, x3f = x3(tf ), x4f = x4(tf )}
  22. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Pr´ ediction • Il n’existe pas d’arc singulier si x2(tf ) − x1(tf ) tan x3f ≥ 0, pour les manœuvres “counter-clockwise”, et si x2(tf ) − x1(tf ) tan x3f ≤ 0, pour les manœuvres “clockwise” • Rappel: condition finale x(tf ) ∈ M1 o` u M1 ={(x1, x2, x3, x4) ∈ R4| x2(tf ) − x1(tf ) tan x3f = 0, x3f = x3(tf ), x4f = x4(tf )}
  23. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Pr´ editction • Remarque: si la condition x2(tf ) − x1(tf ) tan x3f ≥ 0, ou x2(tf ) − x1(tf ) tan x3f ≤ 0, n’est pas v´ erifi´ ee, un nombre plus important d’arcs r´ eguliers peut apparaˆ ıtre
  24. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Resolution num´ erique - Pr´ ediction • Rappel: (x1(tf ), x2(tf )) = (x10, x20) + ∆1(x1, x2) + ∆2(x1, x2) SC := x2(tf ) − x1(tf ) tan x3f SC ≥ 0 ⇒ pas d’arc singulier (cas “counter-clockwise”) SC ≤ 0 ⇒ pas d’arc singulier (cas “clockwise”) • Influence de v0: Soit (x10, x20) := v0 (cos x30, sin x30) ⇒ R´ esoudre SC (v0) = 0 ⇒ v0 = ¯ vup (“counter-clockwise”) ou v0 = ¯ vdown (“clockwise”) ⇒ Si v0 ≤ ¯ vup (“counter-clockwise”) ou v0 ≤ ¯ vdown (“clockwise”) ⇒ pas d’arc singulier
  25. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Solutions sous-optimales (SODM) • Deux m´ ethodes num´ eriques principales • M´ ethodes indirectes: PMP ⇒ probl` eme de tir le contrˆ ole chattering n’est plus continu par morceaux ⇒ ´ echec de l’int´ egration num´ erique ⇒ ne sont pas appropri´ ees • M´ ethodes directes: Discr´ etisation ⇒ programmation non-lin´ eaire • Solution sous-optimale (SODM): m´ ethode directe + ap- proximation par contrˆ ole continu par morceaux • Outils: IPOPT + AMPL
  26. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Convergence de la SODM • Nombre de pas de discr´ etisation: N • Contrˆ ole u(t) converge vers u(t)∗ quand N augmente Figure : Contrˆ ole u(t) avec N different.
  27. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Convergence de la SODM • Nombre de pas de discr´ etisation: N • Temps de manœuvre tf converge vers t∗ f quand N aug- mente. Figure : Temps de manœuvre tf en fonction de N different. • Conjecture: notons tN f , uN, xN correspondant ` a N lim N→+∞ tN f = t∗ f , lim N→+∞ uN = u∗, lim N→+∞ xN = x∗. o` u t∗ f , u∗, x∗ correspondent ` a la solution r´ eelle
  28. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Solutions sous-optimales (SOFSN) • Solution sous-optimale avec un nombre de commutation fix´ e (SOFSN) • Fixer le nombre de commutation: (K − 1) < +∞ ⇒ trajectoire est divis´ ee par K intervalles • Variables: contrˆ ole (constant) sur chaque intervalle + dur´ ee de temps entre deux temps de commutation adja- cent • Trois parties: arc r´ egulier + arc singulier + arc r´ egulier • Outils: IPOPT + AMPL
  29. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Exemple num´ erique - SOFSN • Soit K = 5 u = (1, −1, 0.0724, −1, 1), ζ = (9.6156, 10.4054, 23.3685, 9.9737, 9.0723)s, tf = 62.4355s • Soit K = 7 u =(1, −1, 1, −0.0055, 1, −1, 1), ζ =(9.6350, 11.7880, 2.2523, 16.1259, 2.2205, 11.3196, 9.0889)s, tf = 62.4302s • Comparaison de la SODM et la SOFSN: • SODM avec N = 200: t∗ f = 62.42s ⇒ un peu moins • SOFSN: t∗ f = 62.43s
  30. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Table of content 1 Formulation du probl` eme MTTP 2 Analysis des ext´ emales 3 Solution sous-optimale 4 Conclusions
  31. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions Conclusions • L’existence et l’optimalit´ e locale des extr´ emales de chatter- ing • Les conditions pour lesquelles il n‘y a pas d’arc singulier • Les solutions sous-optimales
  32. Jiamin ZHU Formulation du probl` eme MTTP Analysis des ext´

    emales Application de Principe du Maximum de Pontryagin Existence du chattering Optimalit´ e du chattering Solution sous-optimale Pr´ ediction Solutions sous-optimales Conclusions MERCI DE VOTRE ATTENTION