Time-Optimal Switching Surfaces for Triple Integrator under Full Box Constraints

Résumé technique : Surfaces de commutation optimales en temps pour un triple intégrateur sous contraintes boîte complètes

Énoncé du problème
Cet article traite du problème de commande optimale en temps pour un système triple intégrateur soumis à des contraintes boîte complètes. La dynamique du système fait intervenir des commandes limitées en jerk ($u$) et des contraintes d'état sur la position ($x_1$), la vitesse ($x_2$) et l'accélération ($x_3$). L'objectif est de transférer le système d'un état initial arbitraire $\mathbf{x}_0$ vers un état terminal fixe $\mathbf{x}_f$ en un temps minimal $t_f$.

Bien que la commande optimale en temps pour les double intégrateurs soit bien comprise, le cas du triple intégrateur présente des défis significatifs, particulièrement lorsqu'il s'agit de :

1. Contraintes asymétriques : Où les bornes supérieures et inférieures du jerk et des états ne sont pas égales en magnitude.
2. Conditions aux limites non stationnaires : États initiaux et terminaux arbitraires.
3. Contraintes de position actives : Scénarios où la trajectoire doit toucher ou glisser le long des frontières des contraintes de position ($x_3 = \bar{x}_3$ ou $x_3 = \underline{x}_3$).

Les solutions existantes, telles que l'algorithme commercial Ruckig (version communautaire), échouent souvent sous des contraintes asymétriques ou des contraintes de position actives. Les approches analytiques précédentes ont soit ignoré les contraintes d'état, soit supposé des contraintes symétriques, soit échoué à fournir une caractérisation complète des surfaces de commutation et des conditions actives pour les contraintes de position.

Méthodologie
Les auteurs utilisent le principe du maximum de Pontryagin (PMP) pour caractériser la loi de commande optimale, qui suit une structure Bang-Bang-Singulière (BBS). Le cœur de la méthodologie implique une analyse géométrique de l'espace d'état pour dériver des surfaces de commutation optimales en temps.

1. Loi de commutation augmentée (ASL) : La trajectoire optimale est représentée comme une séquence d'arcs (arcs bang-bang non contraints et arcs contraints où les états sont fixes) et de « marqueurs tangents ». Un marqueur tangent se produit lorsque la contrainte de position devient active de manière tangentielle (c'est-à-dire $x_3 = \bar{x}_3$, $x_2 = 0$, et $x_1 \geq 0$).
2. Construction des surfaces de commutation :
   • Sans contraintes de position : Les auteurs construisent des variétés de commutation par intégration rétrograde depuis l'état terminal $\mathbf{x}_f$. Ils définissent des surfaces correspondant à des ASL spécifiques (par exemple, $\bar{0}\bar{0}\bar{0}$, $\bar{0}\bar{1}\bar{0}\bar{0}$) qui partitionnent l'espace d'état 3D.
   • Avec contraintes de position : L'article dérive les conditions spécifiques selon lesquelles les contraintes de position deviennent actives. Il introduit de nouvelles surfaces de commutation (par exemple, $\bar{0}(\bar{3}, 2)\bar{0}\bar{0}\bar{0}$) qui tiennent compte des trajectoires tangentes à la frontière de position. Le théorème 2 fournit une classification complète de toutes les formes ASL possibles pour le triple intégrateur sous contraintes boîte complètes, établissant que les trajectoires optimales consistent en des séquences d'arcs non contraints et de marqueurs tangents.
3. Implémentation algorithmique :
   • Algorithme 1 : Détermine la direction relative entre $\mathbf{x}_0$ et $\mathbf{x}_f$ en analysant les projections d'ordre 2, gérant les cas où les surfaces de commutation sont multivaluées.
   • Algorithme 2 : Résout le problème sans contraintes de position en connectant les états via des commandes optimales d'ordre 2 et en résolvant des équations polynomiales pour les durées des arcs.
   • Algorithme 3 : Étend la solution aux contraintes boîte complètes. Si la solution non contrainte viole les contraintes de position, l'algorithme connecte récursivement l'état à un marqueur tangent (dérivé du théorème 2) et résout les sous-problèmes.
   • Le système d'équations pour les durées des arcs est résolu en utilisant des bases de Gröbner, réduisant le problème à une seule équation polynomiale de degré $\leq 6$.

Contributions clés

• Caractérisation complète des surfaces de commutation : Pour la première fois, l'article fournit une caractérisation analytique complète des surfaces de commutation optimales en temps pour le triple intégrateur sous contraintes boîte complètes, incluant la structure géométrique des solutions avec contraintes de position actives.
• Déduction des conditions actives : L'article dérive explicitement les conditions selon lesquelles les contraintes de position deviennent actives (la condition de « marqueur tangent »), une lacune précédemment absente de la littérature.
• Théorème 2 (Classification ASL) : Il établit l'ensemble complet des lois de commutation augmentées autorisées pour ce problème, prouvant que les trajectoires optimales suivent des séquences spécifiques impliquant des marqueurs tangents.
• Algorithme efficace : Un algorithme est proposé qui planifie des trajectoires optimales en temps pour des états aux limites arbitraires et des contraintes asymétriques.

Résultats
La méthode proposée a été évaluée par rapport à trois références : Ruckig (version communautaire), CasADi (basé sur IPOPT) et SCP (Programmation Convexe Séquentielle).

• Taux de réussite : La méthode proposée a atteint un taux de réussite de 100 % pour trouver des solutions réalisables à la fois pour les cas symétriques (Cas 1) et asymétriques avec contraintes de position actives (Cas 2). En revanche, Ruckig a échoué dans le Cas 2 en raison de l'asymétrie, et SCP/CasADi ont montré des taux de réussite significativement plus faibles (34,1 % et 97,7 % respectivement dans le Cas 2).
• Efficacité computationnelle : Le temps de calcul moyen par trajectoire était d'environ 10 $\mu$s (spécifiquement $2,38 \times 10^{-6}$ s pour le Cas 1 et $1,25 \times 10^{-5}$ s pour le Cas 2). Cela représente une réduction de 5 ordres de grandeur par rapport aux références basées sur l'optimisation (CasADi : ~0,28 s ; SCP : ~7,28 s).
• Optimalité : La méthode a systématiquement trouvé des solutions globalement optimales. Dans les cas de test où les méthodes basées sur l'optimisation convergaient vers des optima locaux (ASL sous-optimaux), la méthode proposée a réduit le temps terminal jusqu'à 90,9 % (Cas 1) et 73,4 % (Cas 2).

Signification
Les auteurs affirment que ce travail fournit la première caractérisation complète des surfaces de commutation et un algorithme efficace pour le problème de commande optimale en temps d'ordre 3 sous des conditions aux limites non stationnaires et des contraintes de position actives. En s'éloignant de l'optimisation itérative vers une solution géométrique et analytique, l'article démontre qu'un contrôle optimal en temps de haute précision peut être atteint avec des temps de calcul de l'ordre de la microseconde, le rendant hautement adapté aux applications industrielles en temps réel telles que le mouvement robotique et l'usinage à commande numérique où les contraintes sont souvent asymétriques et actives.