FR 3034535 A1 20161007 FR 1552736 A 20150331 1 DOMAINE TECHNIQUE La présente invention appartient au domaine du contrôle d'attitude d'engins spatiaux, tels que des satellites, et concerne plus particulièrement un procédé et un système de contrôle d'attitude d'un engin spatial en rotation sur lui-même avec un moment cinétique total initial non nul. Par « contrôle d'attitude », on entend ici plus particulièrement modifier l'orientation de l'engin spatial par rapport à l'axe dudit moment cinétique total initial, dit « axe cinétique », c'est-à-dire aligner, dans un repère associé à la géométrie de l'engin spatial, dit « repère engin », ledit axe cinétique avec un axe prédéterminé en repère engin. ÉTAT DE LA TECHNIQUE Pour modifier l'orientation d'un engin spatial par rapport à l'axe cinétique, il est connu : - d'arrêter la rotation de l'engin spatial en repère inertiel en absorbant ledit moment cinétique total au moyen de volants d'inertie (roues à réaction, actionneurs gyroscopiques), - de placer, au moyen desdits volants d'inertie, l'engin spatial dans l'orientation choisie par rapport à l'axe dudit moment cinétique total, - de transférer le moment cinétique total stocké dans les volants d'inertie à l'engin spatial dans l'orientation choisie. Une telle approche suppose toutefois que le moment cinétique total initial se trouve dans la capacité d'absorption des volants d'inertie, ce qui ne peut pas toujours être garanti. Notamment, le moment cinétique total initial, transféré à un satellite au moment de la séparation avec un lanceur dudit satellite, est généralement trop important pour être absorbé par les volants d'inertie dudit satellite. A titre d'exemple, le moment cinétique total initial d'un satellite après lancement peut être de l'ordre de 500 N.m.s à 1000 N.m.s, alors que la capacité des volants d'inertie embarqués dans un satellite est généralement de l'ordre de 50 N.m.s à 100 N.m.s. En outre, quand bien même les volants d'inertie seraient dimensionnés de sorte à présenter une capacité de 1000 N.m.s, des conditions particulièrement défavorables de séparation pourraient toujours transférer au satellite un moment cinétique initial supérieur à 1000 3034535 2 N.m.s, en dehors de la capacité d'absorption des volants d'inertie. C'est pourquoi, de nos jours, la modification de l'orientation d'un satellite, par rapport à l'axe cinétique après lancement, est généralement précédée d'une diminution dudit moment cinétique total au moyen de 5 propulseurs chimiques, permettant de ramener ledit moment cinétique total dans la capacité d'absorption des volants d'inertie. Il est actuellement envisagé que les futurs satellites ne soient plus équipés de propulseurs chimiques, mais uniquement de propulseurs électriques (plasmiques). Or les propulseurs électriques ne peuvent pas se 10 substituer aux propulseurs chimiques, pour réduire le moment cinétique total d'un satellite après séparation, tant que l'autonomie électrique dudit satellite n'est pas assurée. Or l'autonomie électrique du satellite ne peut être assurée qu'en plaçant ledit satellite dans une orientation adaptée. EXPOSÉ DE L'INVENTION 15 La présente invention a pour objectif de remédier à tout ou partie des limitations des solutions de l'art antérieur, notamment celles exposées ci-avant, en proposant une solution qui permette de modifier l'orientation d'un engin spatial par rapport à l'axe cinétique au moyen uniquement de volants d'inertie, y compris lorsque ledit moment cinétique total initial est supérieur à la capacité 20 d'absorption des volants d'inertie. A cet effet, et selon un premier aspect, l'invention concerne un procédé de contrôle d'attitude d'un engin spatial en rotation sur lui-même avec un moment cinétique total 11ToT non nul, ledit engin spatial comportant un ensemble de volants d'inertie adaptés à former un moment cinétique interne 25 HACT d'axe quelconque dans un repère engin. Le procédé de contrôle d'attitude comporte une étape d'alignement de l'axe dudit moment cinétique total HTOT avec un axe principal d'inertie de l'engin spatial, au cours de laquelle les volants d'inertie sont commandés de sorte à former un moment cinétique interne HAcT tel que l'expression suivante : HACT X J-1(HTOT 0 r1HTOT) 30 est négative pendant toute la durée de l'étape d'alignement avec l'axe principal d'inertie si ledit axe principal d'inertie visé est l'axe d'inertie maximale de l'engin spatial, 3034535 3 - est positive pendant toute la durée de l'étape d'alignement avec l'axe principal d'inertie si ledit axe principal d'inertie visé est l'axe d'inertie minimale de l'engin spatial, expression dans laquelle J est la matrice d'inertie de l'engin spatial en repère 5 engin, l'opérateur x est le produit scalaire entre deux vecteurs et l'opérateur 0 est le produit vectoriel entre deux vecteurs. En effet, les inventeurs ont constaté que de telles dispositions permettent toujours d'aligner l'axe du moment cinétique total avec l'axe principal d'inertie visé (soit l'axe d'inertie maximale soit l'axe d'inertie 10 minimale). Dans le cas où l'axe principal d'inertie visé est l'axe d'inertie maximale, de telles dispositions permettent en outre d'accélérer significativement la convergence par rapport au cas où aucun contrôle n'est effectué (un engin spatial soumis à un moment cinétique total non nul ayant, du fait de dissipations internes d'énergie et en l'absence de contrôle, tendance à 15 converger lentement vers une orientation dans laquelle l'axe dudit moment cinétique total est aligné avec l'axe d'inertie maximale). En outre, le signe de l'expression HAcr X J-1(1TOT °J-111TOT) est entièrement déterminé par la direction et le sens du moment cinétique interne HAcT, et est indépendant du module dudit moment cinétique interne HAcT. Par 20 conséquent, le procédé de contrôle d'attitude objet de l'invention peut être mis en oeuvre y compris lorsque le moment cinétique total Hun' initial est supérieur à la capacité d'absorption des volants d'inertie de l'engin spatial. En d'autres termes, on peut toujours assurer que le moment cinétique interne HAcT respecte la capacité maximale de stockage de moment cinétique des volants 25 d'inertie et la capacité maximale de formation de couple des volants d'inertie. Dans des modes particuliers de mise en oeuvre, le procédé de contrôle d'attitude peut comporter en outre l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles. 30 Dans des modes particuliers de mise en oeuvre, les volants d'inertie sont commandés de sorte à former un moment cinétique interne HAcT pour lequel un angle 0 entre ledit moment cinétique interne HAcT et le vecteur 3034535 4 1-1(1-1Tcyr 0 1-11froT) vérifie, pendant toute la durée de l'étape d'alignement avec l'axe principal d'inertie, l'expression suivante : Icosel > 0.9 De telles dispositions permettent d'accélérer significativement la 5 convergence vers l'axe principal d'inertie visé, étant entendu que la direction du vecteur 1-1(FIT0T 0 1-111T0T) est la direction optimale pour avoir la convergence la plus rapide vers l'axe principal d'inertie visé. Dans des modes particuliers de mise en oeuvre, les volants d'inertie sont commandés, pendant toute la durée de l'étape d'alignement avec l'axe 10 principal d'inertie, de sorte à former un moment cinétique interne HAcT : HACT = Kv - U expression dans laquelle Kv est un paramètre scalaire négatif si l'axe principal d'inertie visé est l'axe d'inertie maximale de l'engin spatial ou positif si ledit axe principal d'inertie visé est l'axe d'inertie minimale dudit engin spatial, et U correspond au vecteur unitaire : III l(Hurr 01 1HT0T)11 15 Dans des modes particuliers de mise en oeuvre, le procédé de contrôle d'attitude comporte, après l'étape d'alignement avec l'axe principal d'inertie, une étape d'alignement avec un axe X prédéterminé en repère engin, au cours de laquelle les volants d'inertie sont commandés de sorte à asservir les composantes du moment cinétique interne HAcT suivant des axes Y, Z, 20 transverses à l'axe X, sur des consignes respectivement hy et hz déterminées en fonction des composantes de la vitesse de rotation de l'engin spatial suivant lesdits axes Y, Z : - la consigne hy du moment cinétique interne HAcT suivant l'axe Y est déterminée suivant une loi de commande de type proportionnelle- 25 intégrale à partir de la composante r de la vitesse de rotation suivant l'axe Z, - la consigne hz du moment cinétique interne HAcT suivant l'axe Z est déterminée suivant une loi de commande de type proportionnelle-intégrale à partir de la composante q de la vitesse de rotation = I 1(11ToT 0 i 1Hr:1a) U 3034535 5 suivant l'axe Y. De telles dispositions permettent de réduire la nutation éventuelle de l'engin spatial autour de l'axe principal d'inertie visé. En effet, la loi de commande des volants d'inertie au cours de l'étape d'alignement avec l'axe 5 principal d'inertie est non linéaire permet d'aligner l'axe cinétique avec l'axe d'inertie maximale ou avec l'axe d'inertie minimale quelle que soit la capacité de formation de moment cinétique et de couple des volants d'inertie de l'engin spatial. Toutefois, le vecteur 1-1(Flpyr 0 IlifroT) - tend vers zéro à l'approche de l'axe principal d'inertie visé, de sorte que toute erreur sur la connaissance 10 de la matrice d'inertie J peut se traduire par une nutation résiduelle autour dudit axe principal d'inertie visé. La loi de commande de type proportionnelle-intégrale, linéaire, permet d'assurer la convergence vers un axe X, proche de (dans la limite de la capacité des volants d'inertie) ou confondu avec l'axe principal d'inertie visé, et d'amortir la nutation jusqu'à une valeur nulle autour 15 dudit axe X. Dans des modes particuliers de mise en oeuvre, les consignes hy et hz du moment cinétique interne HAcT suivant les axes respectivement Y et Z sont reliées aux composantes respectivement r et q de la vitesse de rotation de l'engin spatial par les fonctions de transfert suivantes, exprimées dans le 20 domaine de Laplace : { hy = Kz (1 + (j') ) r wS y hz = -Ky (1 + -s ) q expressions dans lesquelles : - s est la variable de Laplace, - Ky et Kz sont des paramètres scalaires de même signe constants au cours du temps, 25 - (1)Y et wz sont des paramètres scalaires positifs constants au cours du temps. Dans des modes particuliers de mise en oeuvre, l'étape d'alignement avec l'axe X est exécutée à moment cinétique total 11ToT constant en repère inertiel. 30 Dans des modes particuliers de mise en oeuvre, l'étape d'alignement 3034535 6 avec l'axe principal d'inertie est exécutée à moment cinétique total HTOT constant en repère inertiel. Dans des modes particuliers de mise en oeuvre, le procédé de contrôle d'attitude comporte, après l'étape d'alignement avec l'axe principal d'inertie, 5 une étape de modification du moment cinétique total HT0T en repère inertiel au moyen de propulseurs de l'engin spatial. Dans des modes particuliers de mise en oeuvre, l'étape d'alignement avec l'axe principal d'inertie est exécutée avant toute modification du moment cinétique total Hun, en repère inertiel. 10 Selon un second aspect, l'invention concerne un produit programme d'ordinateur comportant un ensemble d'instructions de code de programme qui, lorsqu'elles sont exécutées par un processeur, mettent en oeuvre un procédé de contrôle d'attitude selon l'un quelconque des modes de mise en oeuvre de l'invention. 15 Selon un troisième aspect, l'invention concerne un système de contrôle d'attitude d'un engin spatial en rotation sur lui-même avec un moment cinétique total Hun, non nul, ledit engin spatial comportant un ensemble de volants d'inertie adaptés à former un moment cinétique interne HAcT d'axe quelconque dans un repère engin, ledit système de contrôle d'attitude 20 comportant des moyens configurés pour mettre en oeuvre un procédé de contrôle d'attitude selon l'un quelconque des modes de mise en oeuvre de l'invention. PRÉSENTATION DES FIGURES L'invention sera mieux comprise à la lecture de la description suivante, 25 donnée à titre d'exemple nullement limitatif, et faite en se référant aux figures qui représentent : - Figure 1 : une représentation schématique d'un satellite comportant un système de contrôle d'attitude, - Figure 2 : un diagramme illustrant les principales étapes d'un mode 30 particulier de mise en oeuvre d'un procédé de contrôle d'attitude, - Figure 3 : un diagramme illustrant les principales étapes d'un mode préféré de mise en oeuvre d'un procédé de contrôle d'attitude, - Figure 4 : des courbes illustrant les performances d'un procédé de 3034535 7 contrôle d'attitude selon l'invention. Dans ces figures, des références identiques d'une figure à une autre désignent des éléments identiques ou analogues. Pour des raisons de clarté des figures, les éléments représentés ne sont pas nécessairement à l'échelle, 5 sauf mention contraire. DESCRIPTION DÉTAILLÉE DE MODES DE RÉALISATION La présente invention concerne le contrôle d'attitude d'un engin spatial en rotation sur lui-même avec un moment cinétique total Hun' initial non nul. L'axe du moment cinétique total HTOT est désigné ci-après par « axe 10 cinétique » . Par « contrôle d'attitude », on entend ici modifier au moins l'orientation de l'engin spatial par rapport à l'axe cinétique, c'est-à-dire aligner, dans un repère associé à la géométrie de l'engin spatial, dit « repère engin », ledit axe cinétique avec un axe prédéterminé en repère engin. 15 Dans la suite de la description, on considère de manière non limitative le cas où l'engin spatial est un satellite 10 et où le moment cinétique total HTOT initial correspond au moment cinétique transféré audit satellite 10 lors de la séparation avec un lanceur dudit satellite. Par exemple, le satellite 10 a été placé en orbite GTO (« Geostationary Transfer Orbit ») par le lanceur et est 20 destiné à effectuer sa mission en orbite GEO (« Geostationary Orbit »). Rien n'exclut cependant, suivant d'autres exemples, de considérer d'autres types d'engins spatiaux (navette spatiale, sonde, etc.). En outre, l'invention est applicable de manière plus générale à tout engin spatial dont on souhaite modifier l'orientation de l'axe cinétique en repère engin et quelle que 25 soit l'orbite en cours ou finale dudit engin spatial (LEO ou « Low Earth Orbit », MEO ou « Medium Earth Orbit », GTO, GEO, etc.). La figure 1 représente un mode particulier de réalisation d'un satellite 10 comportant un système de contrôle d'attitude. En pratique, le système de contrôle d'attitude peut comporter également d'autres éléments que ceux 30 représentés sur la figure 1, qui sortent du cadre de l'invention. Dans l'exemple non limitatif illustré par la figure 1, le système de contrôle d'attitude comporte un ensemble de volants d'inertie 20 adapté à former un moment cinétique interne HAcT d'axe quelconque en repère engin. 3034535 8 Tel qu'illustré par la figure 1, le système de contrôle d'attitude comporte par exemple au moins trois volants d'inertie 20, tels que des roues de réaction et/ou des actionneurs gyroscopiques, adaptés à former un moment cinétique interne HAcT d'axe quelconque en repère engin. Par exemple, le système de 5 contrôle d'attitude comporte trois roues de réaction d'axes de vecteurs unitaires respectifs linéairement indépendants. Le système de contrôle d'attitude du satellite 10 comporte également un dispositif 30 de commande. Le dispositif 30 de commande contrôle l'attitude du satellite 10 et commande à cet effet, notamment, les volants d'inertie 20. 10 Le dispositif 30 de commande comporte par exemple au moins un processeur et au moins une mémoire électronique dans laquelle est mémorisé un produit programme d'ordinateur, sous la forme d'un ensemble d'instructions de code de programme à exécuter pour mettre en oeuvre les différentes étapes d'un procédé 50 de contrôle d'attitude du satellite 10. Dans une variante, le 15 dispositif 30 de commande comporte également un ou des circuits logiques programmables, de type FPGA, PLD, etc., et/ou circuits intégrés spécialisés (ASIC) adaptés à mettre en oeuvre tout ou partie desdites étapes du procédé 50 de contrôle d'attitude du satellite 10. En d'autres termes, le dispositif 30 de commande comporte un ensemble de moyens configurés de façon logicielle 20 (produit programme d'ordinateur spécifique) et/ou matérielle (FPGA, PLD, ASIC, etc.) pour mettre en oeuvre les différentes étapes d'un procédé 50 de contrôle d'attitude décrit ci-après. Dans l'exemple illustré par la figure 1, le dispositif 30 de commande est embarqué dans le satellite 10. De manière plus générale, le dispositif 30 de 25 commande peut être embarqué dans le satellite 10 ou dans un ou plusieurs équipements distants dudit satellite 10, notamment des équipements terrestres (les volants d'inertie 20 étant le cas échéant commandés à distance). Rien n'exclut en outre, suivant d'autres exemples, d'avoir un dispositif 30 de commande distribué entre le satellite 10 et un ou plusieurs autres équipements 30 distants dudit satellite 10. Dans le cas où le dispositif 30 de commande est au moins partiellement embarqué dans un équipement distant du satellite 10, ledit équipement distant et le satellite 10 comportent des moyens conventionnels respectifs de communication à distance. 3034535 9 Dans son principe, un procédé 50 de contrôle d'attitude selon l'invention comporte principalement une étape 52 d'alignement de l'axe dudit moment cinétique total 11ToT avec un axe principal d'inertie du satellite 10, qui met une loi de commande non linéaire. L'axe principal d'inertie visé est soit 5 l'axe d'inertie maximale (« flat spin ») du satellite 10, soit l'axe d'inertie minimale dudit satellite 10. Au cours de l'étape 52 d'alignement avec l'axe principal d'inertie, le dispositif 30 de commande pilote les volants d'inertie 20 de sorte à former un moment cinétique interne HAcT tel que l'expression suivante : HACT X J-1(u 0 J-1HTOT) 10 est de signe constant pendant toute la durée de ladite étape 52 d'alignement avec l'axe principal d'inertie. Dans l'expression précédente, J est la matrice d'inertie du satellite 10 en repère engin, l'opérateur x est le produit scalaire entre deux vecteurs et l'opérateur CD est le produit vectoriel entre deux vecteurs. Plus particulièrement, 15 les volants d'inertie 20 sont commandés de sorte que l'expression précédente : - est négative pendant toute la durée de l'étape 52 d'alignement avec l'axe principal d'inertie si ledit axe principal d'inertie visé est l'axe d'inertie maximale du satellite 10, - est positive pendant toute la durée de l'étape 52 d'alignement avec 20 l'axe principal d'inertie si ledit axe principal d'inertie visé est l'axe d'inertie minimale du satellite 10. Le signe de l'expression HAcr X J-1(HTOT litroT) est entièrement déterminé par la direction et le sens du moment cinétique interne HAcT, et est indépendant du module dudit moment cinétique interne HAcT. Par conséquent, 25 le procédé 50 de contrôle d'attitude peut être mis en oeuvre y compris lorsque le moment cinétique total Hun' initial est supérieur à la capacité d'absorption des volants d'inertie de l'engin spatial. En d'autres termes, on peut toujours assurer que le moment cinétique interne HAcT, à former au cours de l'étape 52 d'alignement avec l'axe principal d'inertie, respecte la capacité maximale de 30 formation de moment cinétique et de couple des volants d'inertie 20. Le système de contrôle d'attitude comporte par exemple un dispositif 303 4 53 5 10 de mesure (non représenté sur les figures) adapté à mesurer la vitesse de rotation inertielle du satellite 10, à partir de laquelle le dispositif 30 de commande peut déterminer de manière conventionnelle le moment cinétique total Hun, initial du satellite 10. 5 Ensuite, le dispositif 30 de commande détermine par exemple les composantes du moment cinétique total Hun, dans le repère engin à l'instant considéré, et en déduit les conditions que le moment cinétique interne HAcT à former doit vérifier pour assurer que le signe de l'expression HAcr X ri (11TOT 0 1-11froT) est négatif si l'axe principal d'inertie visé est l'axe d'inertie 10 maximale, positif si l'axe principal d'inertie visé est l'axe d'inertie minimale. Le dispositif 30 de commande détermine alors une consigne adaptée de moment cinétique interne HAcT, ainsi qu'une commande correspondante de couple à former par lesdits volants d'inertie 20 pour obtenir ladite consigne du moment cinétique interne HAcT. La valeur réelle du moment cinétique interne HAcT est 15 par exemple déterminée au moyen de mesures des vitesses de rotation respectives des volants d'inertie 20, et l'écart entre la valeur réelle et la consigne du moment cinétique interne HAcT est par exemple utilisé, de manière conventionnelle, pour actualiser les commandes de couple à former par lesdits volants d'inertie 20. Ces différentes étapes sont itérées au cours du temps afin 20 de tenir compte de la rotation du satellite 10 (et des variations des composantes du moment cinétique total Hun, dans le repère engin) pour assurer que le signe de l'expression HAcT X J-1(1TOT 0 riliTOT) reste constant pendant toute la durée de l'étape 52 d'alignement de l'axe cinétique avec l'axe principal d'inertie visé. 25 Tel qu'indiqué précédemment, la direction du vecteur 1-1 (1-1Tcyr 01-11froT) est la direction optimale pour avoir la convergence la plus rapide vers l'axe principal d'inertie visé. Dans des modes préférés de mise en oeuvre, les volants d'inertie 20 sont donc commandés de sorte à former un moment cinétique interne HAcT 30 dont la direction, pendant toute la durée de l'étape 52 d'alignement avec l'axe principal d'inertie, n'est pas trop éloignée de celle du vecteur 1-1 (1-1Tcyr 01-11froT)- Plus particulièrement, les volants d'inertie 20 sont 303 4 53 5 11 commandés de sorte que la valeur absolue du cosinus d'un angle 0 entre ledit moment cinétique interne HAcT et le vecteur 1-1(-1Tcyr 1-11froT) est, pendant toute la durée de l'étape 52 d'alignement avec l'axe principal d'inertie, supérieure à 0.9 (lcosel > 0.9), voire supérieure à 0.98 (lcosel > 0.98). 5 De préférence, les volants d'inertie 20 sont commandés, pendant toute la durée de l'étape 52 d'alignement avec l'axe principal d'inertie, de sorte à former un moment cinétique interne HAcT sensiblement de même direction que celle du vecteur J 10-1ToT 0 J-11-1T0T) HACT = KvU expression dans laquelle Kv est un paramètre scalaire et U correspond à la 10 direction du vecteur J(H- 1.--ToT 01-11-1T0T) = I 1(11ToT 111ToT) U III 1-(HT0T 1HT0T)11 Le paramètre scalaire Kv est soit négatif si l'axe principal d'inertie visé est l'axe d'inertie maximale du satellite 10, soit positif si ledit axe principal d'inertie visé est l'axe d'inertie minimale dudit satellite 10. Le paramètre scalaire Kv est constant ou variable au cours du temps, 15 et sa valeur est déterminée de sorte à assurer que le moment cinétique interne HAcT, à former au cours de l'étape 52 d'alignement avec l'axe principal d'inertie, respecte toujours la capacité maximale de formation de moment cinétique et de couple des volants d'inertie 20. Par exemple, la valeur du paramètre scalaire Kv vérifie les expressions suivantes : Tmax 'Kyi il(v1 Hmax 20 expressions dans lesquelles : - Tmax correspond à la capacité maximale de formation de couple des volants d'inertie 20, - Û correspond à la dérivée temporelle de la direction U du vecteur J-10- froT J-1 IITOT) 25 - Hmax correspond à la capacité maximale de formation de moment cinétique des volants d'inertie 20. La valeur du paramètre du paramètre scalaire Kv peut également être 3034535 12 limitée lorsque la direction du moment cinétique total Hun, se rapproche de l'axe principal d'inertie. Par exemple, lorsque l'angle cp entre la direction du moment cinétique total Hun, et l'axe principal d'inertie devient inférieur à 15°, il est possible d'appliquer un facteur simp/sin(15°) qui décroît progressivement à 5 mesure que l'angle cp diminue. La figure 2 représente un mode préféré de mise en oeuvre d'un procédé 50 de contrôle d'attitude du satellite 10. Tel qu'illustré par la figure 2, le procédé 50 de contrôle d'attitude illustré comporte, après l'étape 52 d'alignement avec l'axe principal d'inertie 10 visé (axe d'inertie maximale ou axe d'inertie minimale), une étape 54 d'alignement de l'axe cinétique avec un axe X prédéterminé en repère engin, au cours de laquelle on utilise une loi de commande différente de celle utilisée au cours de l'étape 52 d'alignement avec l'axe principal d'inertie visé, en l'occurrence une loi de commande linéaire de type proportionnelle-intégrale. 15 Plus particulièrement, au cours de l'étape 54 d'alignement avec l'axe X, les volants d'inertie 20 sont commandés de sorte à asservir les composantes du moment cinétique interne HAcT suivant des axes Y, Z, transverses à l'axe X et formant avec celui-ci un repère engin particulier, sur des consignes respectivement hy et hz déterminées en fonction des 20 composantes de la vitesse de rotation inertielle du satellite 10 suivant lesdits axes Y, Z: - la consigne hy du moment cinétique interne HAcT suivant l'axe Y est déterminée suivant une loi de commande de type proportionnelle-intégrale à partir de la composante r de la vitesse de rotation 25 suivant l'axe Z, - la consigne hz du moment cinétique interne HAcT suivant l'axe Z est déterminée suivant une loi de commande de type proportionnelle-intégrale à partir de la composante q de la vitesse de rotation suivant l'axe Y. 30 Tel qu'illustré par la figure 2, l'axe X visé est de préférence confondu avec l'axe principal d'inertie visé, auquel cas la loi de commande de type proportionnelle-intégrale permet d'amortir la nutation jusqu'à une valeur nulle 3034535 13 autour dudit axe principal d'inertie visé. Toutefois, la loi de commande de type proportionnelle-intégrale permet également d'assurer la convergence, à partir d'un axe cinétique sensiblement aligné avec l'axe principal d'inertie visé, vers un axe X distinct dudit axe principal d'inertie, et d'amortir la nutation jusqu'à 5 une valeur nulle autour dudit axe X. Le cas échéant, l'axe X doit cependant être suffisamment proche de l'axe principal d'inertie pour que l'alignement de l'axe cinétique avec ledit axe X puisse être effectué au moyen des volants d'inertie 20. L'axe X peut être considéré comme suffisamment proche de l'axe principal d'inertie lorsque les produits d'inertie Ixy et Ixz avec l'axe X vérifient les 10 expressions suivantes : PX& II HTOT II G IIHmaxil Ix expressions dans lesquelles Ix correspond à l'inertie suivant l'axe X. Toutefois, il est préférable de prévoir une marge supplémentaire pour assurer que l'alignement de l'axe cinétique avec ledit axe X puisse effectivement être effectué au moyen des volants d'inertie 20 : 11xy I - IIFIT0T II 1), par exemple égal à deux (p = 2). Dans des modes particuliers de mise en oeuvre, les consignes hy et hz du moment cinétique interne HAcT suivant les axes respectivement Y et Z sont reliées aux composantes respectivement r et q de la vitesse de rotation du 20 satellite 10 par les fonctions de transfert suivantes, exprimées dans le domaine de Laplace : { hy = Kz (1 + (j') )r ws y hz = -Ky (1 + -s ) q expressions dans lesquelles : G IIHmaxil Ix Ilxzl ' Il HTOT II 3034535 14 - s est la variable de Laplace, - Ky et Kz sont des paramètres scalaires sans dimension de même signe constants au cours du temps, - wy et wz sont des paramètres scalaires positifs, homogènes à des 5 pulsations (s-1), constants au cours du temps. L'utilisation du domaine de Laplace dans le domaine des systèmes d'asservissement est tout à fait conventionnelle et considérée comme connue de l'homme de l'art. En outre, l'ajustement des paramètres scalaires KY, Kz, wy et wz, étant données les expressions ci-dessus des consignes hy et hz du 10 moment cinétique interne HAcT, est considéré à la portée de l'homme de l'art. Tel qu'indiqué précédemment, le procédé 50 de contrôle d'attitude permet d'aligner l'axe cinétique avec un axe principal d'inertie (axe d'inertie maximale ou axe d'inertie minimale) quelle que soit la capacité de formation de moment cinétique et de couple des volants d'inertie 20 du satellite 10. 15 Par conséquent, l'étape 52 d'alignement de l'axe cinétique avec l'axe principal d'inertie visé peut être exécutée avant toute modification du moment cinétique total HTOT en repère inertiel. De même, l'étape 54 d'alignement de l'axe cinétique avec un axe X confondu avec l'axe principal d'inertie visé ou proche de celui-ci, peut également être exécutée avant toute modification du 20 moment cinétique total HTOT en repère inertiel. Rien n'exclut cependant, suivant d'autres exemples, de modifier le moment cinétique total HTOT avant et/ou pendant d'aligner l'axe cinétique avec un axe prédéterminé en repère engin. La figure 3 représente un mode préféré de mise en oeuvre d'un procédé 50 de contrôle d'attitude du satellite 10. Tel qu'illustré par la figure 3, le 25 procédé 50 de contrôle d'attitude comporte, après l'étape 52 d'alignement avec l'axe principal d'inertie visé, une étape 56 de modification du moment cinétique total HTOT en repère inertiel au moyen de propulseurs (non représentés sur les figures) du satellite 10, commandés par le dispositif 30 de commande. L'étape 56 de modification du moment cinétique total HTOT au moyen de propulseurs du 30 satellite 10 est conventionnelle et est considérée comme connue de l'homme de l'art. Par exemple, l'étape 56 de modification du moment cinétique total HTOT peut comporter la réduction dudit moment cinétique total HTOT jusqu'à ce qu'il 3034535 15 soit dans la capacité d'absorption des volants d'inertie 20. Dans l'exemple non limitatif illustré par la figure 3, le procédé 50 de contrôle d'attitude comporte également l'étape 54 d'alignement avec un axe X, qui met en oeuvre une loi de commande linéaire pour amortir la nutation, qui 5 est par exemple également exécutée avant l'étape 56 de modification du moment cinétique total Hun, en repère inertiel. Dans l'exemple illustré par la figure 3, l'axe X est considéré confondu avec l'axe principal d'inertie visé. De préférence, les propulseurs utilisés au cours de l'étape 56 de modification du moment cinétique total HTOT sont des propulseurs électriques. 10 Avantageusement, l'axe principal d'inertie visé (ou, le cas échéant, l'axe X visé) est un axe suivant lequel le satellite 10 peut assurer son autonomie électrique au moyen de générateurs solaires (non représentés sur les figures), de sorte qu'il est possible d'utiliser des propulseurs électriques pour modifier le moment cinétique total Hun, en repère inertiel. 15 Dans la suite de la description, on se place de manière non limitative, à titre d'exemple, dans le cas où les générateurs solaires du satellite 10 sont mobiles en rotation autour d'un axe de rotation sensiblement orthogonal à l'axe d'inertie maximale dudit satellite 10 (ce sera généralement le cas si les générateurs solaires sont agencés de part et d'autre d'un corps dudit satellite 20 10). L'orientation desdits générateurs solaires autour de l'axe de rotation est commandée au moyen de mécanismes d'entrainement. Dans ce cas, l'axe principal d'inertie visé au cours de l'étape 52 d'alignement (et, le cas échéant, l'axe X visé au cours de l'étape 54 d'alignement) correspond à l'axe d'inertie maximale dudit satellite 10. A l'issue de cette étape, l'axe cinétique est donc 25 sensiblement orthogonal à l'axe de rotation des générateurs solaires. On peut montrer qu'il est alors toujours possible de trouver une orientation des générateurs solaires permettant d'assurer que l'ensoleillement moyen de surfaces photosensibles desdits générateurs solaires, sur la durée d'une rotation complète du satellite 10 sur lui-même, est supérieur à une valeur seuil 30 prédéfinie, par exemple supérieur à 30%. En plaçant les générateurs solaires dans une telle orientation, il est alors possible d'assurer l'autonomie électrique du satellite 10, de sorte que des propulseurs électriques peuvent être mis en oeuvre pour modifier le moment cinétique total Hun' du satellite 10. 3034535 16 La description ci-avant illustre clairement que par ses différentes caractéristiques et leurs avantages, la présente invention atteint les objectifs qu'elle s'était fixés. En particulier, l'invention permet de modifier l'orientation d'un satellite 10 par rapport à un axe cinétique au moyen uniquement de 5 volants d'inertie 20, y compris lorsque le moment cinétique total 11ToT initial est supérieur à la capacité d'absorption desdits volants d'inertie 20, grâce à une loi de commande non linéaire. La figure 4 représente des résultats de simulation illustrant les performances d'un procédé 50 de contrôle d'attitude selon l'invention au cours 10 de l'étape 52 d'alignement de l'axe du moment cinétique total HTOT avec l'axe principal d'inertie d'un satellite 10. Les inerties du satellite 10 considéré sont les suivantes : - Ix = 70000 kg.m2, - ly = 10000 kg.m2, 15 - lz = 50000 kg.m2. de sorte que l'axe principal d'inertie dudit satellite 10 est sensiblement l'axe X. En outre, la capacité maximale Hmax en moment cinétique des volants d'inertie 20 du satellite 10 considéré est égale à 25 N.m.s, et la capacité maximale Tmax en couple est égale à 0.2 N.m. 20 La partie a) de la figure 4 représente les vitesses de rotation ,Qx, l/y, fz (exprimées en rad/s) du satellite 10. On constate que le moment cinétique total 11ToT initial est proche de l'axe Y. En outre, le moment cinétique total 11ToT initial considéré est de 250 N-ms, de sorte qu'il est très supérieur à la capacité maximale Hmax en moment cinétique des volants d'inertie 20 du satellite 10. 25 La partie b) de la figure 4 représente les consignes hx, hy, hz du moment cinétique interne HAcT au cours de l'étape 52 d'alignement. Tel qu'illustré par la partie a) de la figure 4, l'axe du moment cinétique total 11ToT converge bien vers l'axe X du satellite 10, et sans que le moment cinétique interne HACT ne dépasse la capacité maximale Hmax de formation de 30 moment cinétique des volants d'inertie 20 du satellite 10, comme le montre la partie b) de la figure 4. REVENDICATIONS1 - Procédé (50) de contrôle d'attitude d'un engin spatial (10) en rotation sur lui-même avec un moment cinétique total 11ToT non nul, ledit engin spatial (10) comportant un ensemble de volants d'inertie (20) adaptés à former un moment cinétique interne HAcT d'axe quelconque dans un repère engin, caractérisé en ce qu'il comporte une étape (52) d'alignement de l'axe dudit moment cinétique total 11ToT avec un axe principal d'inertie de l'engin spatial (10), au cours de laquelle les volants d'inertie (20) sont commandés de sorte à former un moment cinétique interne HAcT tel que l'expression suivante : HACT X J-1(u 0 J-1uTOT) - est négative pendant toute la durée de l'étape (52) d'alignement avec l'axe principal d'inertie si ledit axe principal d'inertie visé est l'axe d'inertie maximale de l'engin spatial (10), - est positive pendant toute la durée de l'étape (52) d'alignement avec l'axe principal d'inertie si ledit axe principal d'inertie visé est l'axe d'inertie minimale de l'engin spatial (10), expression dans laquelle J est la matrice d'inertie de l'engin spatial (10), l'opérateur x est le produit scalaire entre deux vecteurs et l'opérateur CD est le produit vectoriel entre deux vecteurs. 2 - Procédé (50) selon la revendication 1, dans lequel les volants d'inertie (20) sont commandés de sorte à former un moment cinétique interne HAcT pour lequel un angle 0 entre ledit moment cinétique interne HAcT et le vecteur J 1(1froT rittroT) vérifie, pendant toute la durée de l'étape (52) d'alignement avec l'axe principal d'inertie, l'expression suivante : Icosel > 0.9 3 - Procédé (50) selon la revendication 2, dans lequel les volants d'inertie (20) sont commandés, pendant toute la durée de l'étape (52) d'alignement avec l'axe principal d'inertie, de sorte à former un moment cinétique interne HACT HACT = Kv - U expression dans laquelle Kv est un paramètre scalaire négatif si l'axe 3034535 18 principal d'inertie visé est l'axe d'inertie maximale de l'engin spatial (10) ou positif si ledit axe principal d'inertie visé est l'axe d'inertie minimale dudit engin spatial, et U correspond au vecteur unitaire : = 1-1(HT0T 01-111r0T) U 111-1(-froT 01 -111T0T) II 4 - Procédé (50) selon l'une des revendications précédentes, comportant, 5 après l'étape (52) d'alignement avec l'axe principal d'inertie, une étape (54) d'alignement avec un axe X prédéterminé en repère engin, au cours de laquelle les volants d'inertie (20) sont commandés de sorte à asservir les composantes du moment cinétique interne HAcT suivant des axes Y, Z, transverses à l'axe X, sur des consignes respectivement hy et hz 10 déterminées en fonction des composantes de la vitesse de rotation de l'engin spatial (10) suivant lesdits axes Y, Z : - la consigne hy du moment cinétique interne HAcT suivant l'axe Y est déterminée suivant une loi de commande de type proportionnelle-intégrale à partir de la composante r de la vitesse de rotation 15 suivant l'axe Z, - la consigne hz du moment cinétique interne HAcT suivant l'axe Z est déterminée suivant une loi de commande de type proportionnelle-intégrale à partir de la composante q de la vitesse de rotation suivant l'axe Y. 20 5 - Procédé (50) selon la revendication 4, dans lequel les consignes hy et hz du moment cinétique interne HAcT suivant les axes respectivement Y et Z sont reliées aux composantes respectivement r et q de la vitesse de rotation par les fonctions de transfert suivantes, exprimées dans le domaine de Laplace : { hy = Kz (1 + -wsz) r hz = -Ky (1 + (1`7) q 25 expressions dans lesquelles : - s est la variable de Laplace, - Ky et Kz sont des paramètres scalaires de même signe constants au cours du temps, 3034535 19 - wy et wz sont des paramètres scalaires positifs constants au cours du temps. 6 - Procédé (50) selon l'une des revendications 4 à 5, dans lequel l'étape (54) d'alignement avec l'axe X est exécutée à moment cinétique total HTOT 5 constant en repère inertiel. 7 - Procédé (50) selon l'une des revendications précédentes, dans lequel l'étape (52) d'alignement avec l'axe principal d'inertie est exécutée à moment cinétique total 11ToT constant en repère inertiel. 8 - Procédé (50) selon l'une des revendications précédentes, comportant, 10 après l'étape (52) d'alignement avec l'axe principal d'inertie, une étape (56) de modification du moment cinétique total 11ToT en repère inertiel au moyen de propulseurs de l'engin spatial (10). 9 - Produit programme d'ordinateur caractérisé en ce qu'il comporte un ensemble d'instructions de code de programme qui, lorsqu'elles sont 15 exécutées par un processeur, mettent en oeuvre un procédé (50) de contrôle d'attitude selon l'une des revendications précédentes. 10- Système de contrôle d'attitude d'un engin spatial (10) en rotation sur lui-même avec un moment cinétique total 11ToT non nul, ledit engin spatial comportant un ensemble de volants d'inertie (20) adaptés à former un 20 moment cinétique interne HAcT d'axe quelconque dans un repère engin, caractérisé en ce qu'il comporte des moyens configurés pour commander les volants d'inertie (20) conformément à l'une des revendications 1 à 8.