Résumé technique : Le problème de la masse initiale maximale dans la conception de trajectoires à faible poussée

Définition du problème

L'article traite d'une lacune dans la théorie du contrôle optimal à faible poussée : l'absence de formalisation explicite et de traitement analytique du problème de la Masse Initiale Maximale (MIM - Maximum Initial Mass). Bien que le domaine soit dominé par les formulations de minimisation de la propulsion (maximisation de la masse finale) et de minimisation du temps de transfert, le problème MIM pose une question distincte : pour une géométrie de transfert donnée, un état final fixé et un temps de transfert donné, quelle est la plus grande masse initiale à partir de laquelle un engin spatial à faible poussée peut accomplir avec succès le transfert ?

Conceptuellement, cela inverse le paradigme standard « étant donné une masse initiale, minimiser le temps/la propulsion » par « étant donné un temps de vol, maximiser la masse initiale ». Bien que cette perspective ait été utilisée de manière informelle dans des compétitions (par exemple, les GTOC) et pour des études de faisabilité de l'exploitation minière d'astéroïdes, elle n'a pas été rigoureusement établie comme une formulation de contrôle optimal autonome possédant ses propres conditions nécessaires et propriétés structurelles.

Méthodologie

Les auteurs utilisent le Principe du Minimum de Pontryagin (PMP) pour dériver les conditions nécessaires de premier ordre pour le problème MIM et les comparer structurellement au problème classique de temps minimum.

1. Formulation :

   • Dynamique : Le système est modélisé à l'aide de la dynamique standard de faible poussée où le vecteur d'état $\mathbf{x}$ (par exemple, les éléments équinociaux modifiés) et la masse $m$ évoluent sous une commande de poussée $u \in [0,1]$ et une direction de poussée $\hat{\mathbf{i}}$.
   • Objectif : Maximiser la masse initiale $m(t_0)$ sous réserve d'un temps de transfert fixe $t_f - t_0$ et de conditions aux limites $\mathbf{x}(t_0) = \mathbf{x}_0, \mathbf{x}(t_f) = \mathbf{x}_f$. Cela est converti en une minimisation de $J = -m(t_0)$.
   • Analyse du Hamiltonien : Les auteurs dérivent le Hamiltonien et sa fonction de commutation associée. Une conclusion clé est que la poussée optimale est toujours $u^*=1$ (poussée pleine) dans le cas non dégénéré, et que la direction de la poussée suit une loi du vecteur primaire ($\hat{\mathbf{i}}^* \propto -\mathbf{B}(\mathbf{x})^T\boldsymbol{\lambda}$). Cela produit une structure de contrôle continue, évitant les profils « bang-bang » discontinus souvent observés dans d'autres formulations, ce qui facilite la convergence rapide pour les méthodes numériques directes et indirectes.

2. Transversalité et conditions de jauge :

   • L'article réexamine de manière critique la condition de transversalité pour les problèmes à temps libre. Il soutient que l'hypothèse courante selon laquelle le Hamiltonien doit être nul au temps final ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) = 0$) pour les problèmes de temps minimum est un choix de jauge plutôt qu'une condition nécessaire dérivée du PMP pour les optima de bord. La véritable condition est une inégalité unidirectionnelle ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) \le 0$).
   • Cette analyse permet la normalisation du vecteur des co-états $\lambda_m(t_0) = 1$ sans violer les conditions nécessaires.

3. Correspondance extrémale :

   • En utilisant les conditions dérivées, les auteurs prouvent le Théorème 1 : Il existe une correspondance biunivoque entre les extrémaux du problème de temps minimum (avec masse initiale fixe) et le problème de masse initiale maximale (avec temps de transfert fixe). Plus précisément, une trajectoire qui est optimale en temps pour une masse donnée est également optimale en masse pour le temps de vol résultant, et vice versa.

4. Application de l'homotopie :

   • Les auteurs appliquent la formulation MIM au problème de temps minimum à multiples révolutions. Ils construisent un chemin d'homotopie en faisant varier le nombre de révolutions ($N_{rev}$) tout en suivant la surface de la masse initiale maximale $M(t_f, N_{rev})$. Cette approche est contrastée avec les homotopies traditionnelles de continuation de poussée, qui souffrent souvent de mauvais conditionnement dans les scénarios à multiples révolutions.

Résultats clés

• Propriétés numériques : Le problème MIM produit un problème de valeur limite bien posé avec sept équations et sept inconnues (pour des points finaux fixes), évitant les problèmes d'indétermination d'échelle courants dans les méthodes de tir indirectes de temps minimum. Le contrôle optimal est strictement à pleine poussée, garantissant une convergence numérique fluide.
• Clarification théorique : L'article clarifie que la condition d'égalité du Hamiltonien terminal n'est pas une condition nécessaire pour l'optimalité temporelle mais un choix de normalisation. Cela résout une « erreur commune » dans la littérature et valide l'utilisation de conditions de jauge alternatives pour améliorer la stabilité numérique.
• Visualisation du paysage de solutions : Dans un transfert de référence de l'orbite de transfert géostationnaire (GTO) vers l'orbite géostationnaire (GEO), l'homotopie MIM a réussi à cartographier le paysage de solutions. La méthode a révélé une surface lisse et différentiable de la masse initiale maximale.
  • Elle a identifié plusieurs trajectoires extrémales, incluant des minima et maxima locaux du temps de vol le long des courbes d'iso-masse.
  • Elle a localisé avec succès la solution globale optimale en temps (1500 kg de masse initiale) et a démontré que certaines solutions précédemment rapportées (par exemple, la solution E du benchmark) correspondent à des maxima locaux de temps de vol, et non des minima.
• Efficacité de l'homotopie : Le chemin de continuation dans le cadre MIM (variation de $N_{rev}$ pour suivre les courbes d'iso-masse) a été observé comme étant lisse et bien conditionné, offrant une alternative plus robuste aux méthodes de continuation de poussée pour naviguer dans les paysages complexes à multiples révolutions.

Signification et revendications

L'article soutient que le problème de la Masse Initiale Maximale mérite d'être élevé au rang de formulation de contrôle optimal de « première classe » dans la conception de trajectoires à faible poussée. Sa signification réside dans trois domaines principaux :

1. Primitive de conception complémentaire : Bien que structurellement lié au problème de temps minimum, la formulation MIM offre une perspective distincte qui fixe le temps et libère la masse. Cela fournit un outil complémentaire pour la conception de missions, particulièrement pour l'analyse préliminaire où de grands ensembles de transferts doivent être évalués.
2. Rigueur analytique : En dérivant formellement les conditions nécessaires et en corrigeant la compréhension des conditions de transversalité dans les problèmes d'optimalité temporelle, l'article fournit une base théorique plus rigoureuse pour l'optimisation à faible poussée.
3. Utilité algorithmique : La formulation MIM permet la construction de chemins d'homotopie lisses qui sont moins sujets au mauvais conditionnement trouvé dans les méthodes de continuation de poussée traditionnelles. Cela permet l'exploration systématique d'espaces de solutions complexes (tels que les transferts à multiples révolutions) et la génération de vastes bases de données de trajectoires optimales pour la modélisation de substitution (surrogate modeling) et l'analyse de l'accessibilité (reachability analysis).

Les auteurs concluent que l'adoption de la formulation MIM comme standard de conception peut améliorer l'efficacité et la robustesse de la conception de missions à faible poussée, particulièrement pour les approches basées sur les données et l'analyse de mission préliminaire.