FR 2717149 A1 19950915 FR 9402852 A 19940311 l L'invention concerne le contrôle d'attitude, notamment en roulis- lacet, à l'aide de générateurs de moment magnétique, d'un satellite destiné à être stabilisé selon trois axes sur une orbite terrestre. Elle vise de préférence,5 mais non exclusivement, le cas o l'orbite a une inclinaison intermédiaire (typiquement entre environ 15 et environ 75o). L'utilisation de bobines magnétiques pour participer à la stabilisation en attitude d'un satellite en orbite autour de la Terre a déjà été proposée dans de10 nombreux documents parmi lesquels on peut citer: - Reusable Reentry Satellite (RRS) System Design Study Appendix E: Attitude Control System Study Ithaco Incorporated, Ithaca, New York NASA STAR Technical report issue 04 - Magnetic attitude control system for spinning small spacecraft par PAL, PARIMAL, SELBY, VAUGHN H. Ithaco Incorporated, Ithaca New York NASA IAA Conference paper issue 10 - Development of magnetic attitude control system of ETS-III par IMAI F., SHIMIZU J. National space Development Agency of Japan, Tokyo NASA IAA Conference paper, issue 18 An elementary magnetic attitude control system par STICKLER A.C., ALFRIEND K.T. Ithaco Incorporated, Ithaca, New York - Naval Research Laboratory, Washington D.C. AIAA, Mechanics and control flight conference, Anaheim, California/August 5-9, 1974 - Magnetic roll/yaw attitude control of a momentum biased near polar satellite par BRUBERLE E., REINDLET A. Messerschmitt Bâlkow Blohm GmbH, Ottobrunn, Germany IFAC/ESA symposium on automatic control in space, Noordwijkerhout, The Netherlands 5-9 July 1982 - Autonomous spacecraft attitude control using magnetic torquing only par MUSSER K.L., EBERT W.L. The Johns Hopkins University Applied physics laboratory, Laurel Maryland In NASA, Goddard Space Flight Center, Flight Mechanics/Estimation Theory Symposium, 1989 pp 23-38. Dans ces documents, soit les bobines sont utilisées à de simples fins d'amortissement, soit elles n'interviennent que pour désaturer les roues cinétiques qui assurent le contrôle d'attitude, soit elles interviennent en tout ou rien (soit, en variante, en tout "positif", en tout "négatif" ou en rien; les bobines servent alors à désaturer les roues), soit combinent des valeurs du champ magnétique en roulis et des dérivées de valeurs du champ magnétique en tangage, soit se placent dans des hypothèses très restric- tives (par exemple condition BO = 0). Le document de MUSSER et EBERT propose quant à lui un procédé de contrôle d'atti- tude par bobines magnétiques, mettant en oeuvre une modula- tion du courant circulant dans les bobines; cette solution est tout à fait performante mais présente l'inconvénient de nécessiter de disposer d'un modèle complet du champ terrestre et connaître avec précision l'orbite et la position du satellite; d'autre part, elle ne peut pas tenir compte des écarts du champ magnétique terrestre réel par rapport à ce modèle, ni de faire face au moindre orage magnétique. On peut également citer les documents suivants: - US-4.916.622 qui nécessite une roue d'inertie capable de pivoter, - US-4.489.383 qui se limite à des orbites fortement inclinées et n'active les bobines magnétiques que de manière intermittente (sur un quart d'orbite sur deux), - US-4.807.835 qui nécessite l'emploi d'un mât de gravité et n'utilise des bobines magnétiques qu'à des fins d'amortissement, - FR-2.498.155 qui se limite à amortir la nutation, - FR-2.319.150 qui se limite au cas d'orbites synchrones, - FR-2.525.359, qui concerne la désaturation d'un moment cinétique embarqué dans le plan roulis- lacet, - FR-2.525.360, qui nécessite la détermination d'un couple extérieur perturbateur que l'on veut compenser, - FR-2.525.361, qui fait intervenir un signal de polarité magnétique et un organe de lissage, - EP-0.544.242 qui vise à compenser des perturbations magnétiques, - US-5.047.945 qui nécessite la connaissance d'un modèle du champ magnétique terrestre et vise la compensation du champ magnétique résiduel du satellite. L'invention a pour objet de pallier les incon- vénients précités grâce à un procédé et un dispositif de contrôle d'attitude, au moins en roulis-lacet, utilisant l'interaction de générateur(s) de moment magnétique avec le champ magnétique terrestre, permettant d'obtenir une bonne précision de pointage avec toutefois une faible masse, un faible coût, une grande durée de vie grâce à sa simplicité et à sa fiabilité, tout en étant capables de s'adapter à des variations de ce champ magnétique. L'invention propose à cet effet un procédé de contrôle de l'attitude d'un satellite devant être stabilisé selon des axes de roulis, de tangage et de lacet, selon lequel, pour le contrôle en roulis-lacet: - on détermine l'orientation instantanée dans le plan roulis-lacet de la direction instantanée a définie par le produit vectoriel unitaire Y de l'axe de tangage avec le champ magnétique terrestre local B, - on détermine des grandeurs d'attitude en roulis-lacet du satellite, - on détermine l'amplitude et le sens d'un couple de correction Cmxz à appliquer au satellite selon ladite direction, cette amplitude étant déterminée à partir desdites grandeurs d'attitude en roulis-lacet selon une loi de commande prédéterminée dont les gains sont élaborés de façon10 adaptative en sorte de respecter des conditions prédétermi- nées sur les pôles du système modélisant la dynamique du satellite dans le plan roulis-lacet, - on génère dans un plan perpendiculaire à cette direction instantanée un moment magnétique d'amplitude propre à générer selon ladite direction ledit couple de correction Cxz par interaction avec le champ magnétique terrestre, grâce à quoi ce couple de correction du roulis-lacet n'a aucune composante en tangage. Selon des dispositions préférées de l'invention: - les conditions prédéterminées sur les pôles du système sont des conditions de sensible constance au cours du temps, - lesdits gains sont élaborés par interpolation à partir de données stockées à bord et d'au moins l'orientation de la composante Bxz du champ magnétique terrestre instantané dans la plan roulis-lacet, - lesdits gains sont élaborés à partir de coefficients de FOURIER stockés à bord et de la valeur instantanée de l'orientation de la composante Bxz du champ magnétique terrestre local dans le plan roulis-lacet, - la loi de commande est du type proportionnel/dérivée, les grandeurs d'attitude sont les angles et les vitesses angulai- res en roulis-lacet, et les gains sont au nombre de quatre, - la loi de commande est du type proportionnel/intégral, les grandeurs d'attitude sont les angles en roulis-lacet et leurs intégrales et les gains sont au nombre de quatre, - la loi de commande est du type proportionnel/intégral/dérivée, les grandeurs d'attitude sont les angles en roulis-lacet et leurs vitesses et leurs intégrales, et les gains sont au nombre de six, - les pôles auxquels sont appliquées les conditions prédé- terminées sont caractérisés par les grandeurs f et w, telles que le polynôme caractéristique du système puisse s'écrire10 sous la forme d'un produit de n termes: x22 X + 2 i (aiX + wi o i est compris entre 1 (inclus) et n (inclus), - les pôles sont choisis en sorte que les grandeurs wi, sont supérieures d'au moins un ordre de grandeur à la pulsation orbitale, - on détermine l'amplitude de la composante Bxz et on détermine l'amplitude du moment magnétique, noté My, par division de l'amplitude du couple de correction par cette amplitude de la composante Bxz, et on génère ce moment magnétique My selon l'axe de tangage, - on détermine l'amplitude IBI du champ magnétique terrestre B, et on génère le moment magnétique, noté Mb, perpendiculai- rement au champ magnétique terrestre, tel que: b BAnCmxz Mb-= IBI2 - en outre, on contrôle l'attitude en tangage par modifica- tion de la vitesse de rotation d'une roue d'inertie en tangage, - on génère des couples en tangage de désaturation de la roue d'inertie par activation d'au moins un générateur de moment magnétique générant un moment magnétique dans le plan roulis- lacet, - on contrôle l'attitude en tangage par génération dans le plan roulis-lacet d'un moment magnétique d'amplitude détermi- née en fonction des mesures d'attitude en tangage, - on génère ce moment magnétique dans le plan roulis-lacet selon ladite direction instantanée, on détermine les grandeurs d'attitude par détection à l'aide de détecteurs terrestre, solaire et/ou stellaire, de gyroscopes ou de gyromètres, et/ou estimation à l'aide de filtres estimateurs, - on détermine certaines ou moins des grandeurs d'attitude par comparaison de mesures, selon un ou plusieurs axes, du champ magnétique terrestre instantané avec un modèle de ce champ magnétique terrestre, - on détermine le champ magnétique terrestre au moins en roulis-lacet par mesure à l'aide d'au moins un magnétomètre à un ou plusieurs axes de mesure, - on détermine au moins l'orientation du champ magnétique terrestre en roulis-lacet à partir de signaux reçus du sol, - on détermine au moins l'orientation du champ magnétique terrestre local en roulis-lacet à partir d'un modèle de ce champ, - le satellite suit une orbite d'inclinaison comprise entre et 75 environ. L'invention propose également un satellite stabilisé en roulis, en tangage et en lacet, comportant un dispositif de contrôle d'attitude en roulis-lacet comportant: - un premier dispositif de détermination d'au moins l'orientation du champ magnétique terrestre local instantané en roulis-lacet Bxz, - un deuxième dispositif de détermination d'attitude pour fournir des signaux d'attitude au moins en roulis et en lacet, - un calculateur connecté aux premier et second dispositifs, contenant une loi de commande adaptative à gains variables donnant une valeur de couple de correction en roulis/lacet à appliquer perpendiculairement à l'axe de tangage et à l'orientation instantanée de Bxz à partir des signaux d'attitude, et une loi donnant lesdits gains varia- bles en fonction d'au moins l'orientation de Bxz, - un générateur de moment magnétique adapté à générer, dans le plan de l'axe de tangage et de Bxz, un moment magnétique propre à générer par interaction avec le champ magnétique terrestre ladite valeur de couple à appli- quer. Les dispositions préférées précitées s'appli- quent, mutatis mutandis à ce satellite et à son dispositif de contrôle d'attitude. On appréciera que l'invention: - permet de réaliser le contrôle du roulis-lacet sans exercer de couple en tangage, - laisse la possibilité de contrôler aussi le tangage par actionnement magnétique, - ne nécessite pas de moment cinétique pivotant, - ne se limite pas ni aux orbites fortement inclinées, ni aux orbites faiblement inclinées, - n'a pas besoin de mât de gravité, - n'a pas besoin ni de roues de réaction, ni de roue cinéti- que dont l'axe de rotation a une composante non nulle dans le plan roulis-lacet, ni de tuyères, - n'a pas besoin de prévoir un retard avant d'appliquer la commande, - n'a pas besoin de déterminer un éventuel moment perturba- teur extérieur, et ne se limite pas à la compensation de couples perturbateurs, - n'exige pas la prévision du champ magnétique, ni donc la connaissance d'un modèle de ce champ magnétique sur toute l'orbite (la connaissance ponctuelle de ce champ suffit), - n'exige pas la connaissance à bord de la position du satellite sur son orbite. Des objets, caractéristiques et avantages de l'invention ressortent de la description qui suit, donnée à titre d'exemple non limitatif, en regard des dessins annexes sur lesquels: - la figure 1 est une présentation dans l'espace d'un repère lié à un satellite et de divers vecteurs interve- nant dans le contrôle d'attitude au moins en roulis-lacet,5 conformément A l'invention, - la figure 2 est un schéma de principe du procédé de contrôle d'attitude en roulis-lacet de l'inven- tion, - la figure 3 est un groupe de six graphiques représentant les évolutions en fonction du temps des angles de roulis, tangage et lacet et des coordonnées du vecteur rotation instantanée du satellite, en simulation selon l'invention et en l'absence de tout couple perturbateur, - la figure 4 est un groupe de quatre graphiques représentant les évolutions en fonction du temps des quatre gains en roulis, lacet correspondant à la simulation de la figure 3, - la figure 5 est un groupe de six graphiques analogue à la figure 3, mais selon une simulation prenant en compte des couples perturbateurs, - la figure 6 est un groupe de quatre graphiques analogue à celui de la figure 4, montrant les évolutions des gains adaptatifs, - la figure 7 est un groupe de trois graphiques montrant les composantes selon les trois axes des couples perturbateurs pris en considération, d'origine aérodynamique et solaire, et - la figure 8 est un groupe de trois graphiques montrant les évolutions au cours du temps (le long de l'orbite) des composantes selon les trois axes du champ magnétique terrestre instantané que subit le satellite. La figure 1 représente un repère XYZ lié au satellite S que l'on veut stabiliser en attitude, au moins en roulis-lacet. Les axes X, Y et Z ont leurs significations habituelles, en ce sens que: - X est l'axe de roulis, - Y est l'axe de tangage, normalement perpendiculaire au plan orbital, c'est-à-dire au plan roulis-lacet, - Z est l'axe de lacet, normalement dirigé vers le centre de la Terre, avec les précisions suivantes: * le repère X Y Z est direct * X, qui est normalement contenu dans le plan roulis- lacet, est orienté dans le même sens que la vitesse instantanée du satellite sur son orbite T; X est norma- lement tangent à l'orbite lorsque celle-ci est circu- laire. On se préoccupe dans ce qui suit, dans un premier temps, du contrôle d'attitude en roulis-lacet. On appelle e et t les angles de roulis et de lacet, c'est-à-dire les erreurs angulaires que peut présenter l'attitude du satellite autour des axes X et Z par rapport à son attitude de consigne définie ci-dessus par rapport à l'orbite et à la Terre. Le contrôle d'attitude en roulis-lacet est, selon l'invention, réalisé au moyen d'une bobine magnétique orientée en tangage, c'est-à-dire d'axe actif parallèle à l'axe Y, ou plus généralement au moyen d'un générateur de moment magnétique dont l'axe magnétique actif, éventuellement variable, est dans le plan de cet axe Y et du champ magnéti- que terrestre B. Soit My le moment magnétique généré par cette bobine (il est donc orienté selon Y) et Mb le moment magnéti- que généré dans le plan (Y, B) par ce générateur de moment magnétique. On s'intéresse d'abord au cas de la bobine orientée en tangage. L'interaction de My avec la composante en tangage du champ magnétique terrestre local se traduit par un couple nul puisque, comme on le sait, l'interaction entre deux champs ou moments magnétiques est un couple égal au produit vectoriel de ces champs/moments. Par contre My interagit avec la composante Bxz du champ magnétique terrestre dans le plan roulis-lacet, en sorte de générer magnétiquement un couple Cmxz, parallèle au plan roulis-lacet, perpendiculaire à My et Bxz et de module5 égal au produit des modules de ces derniers puisque ceux-ci sont perpendiculaires entre eux. Si l'on s'intéresse maintenant au cas général d'un moment magnétique Mb généré dans le plan (Y, B) on constate que l'interaction de ce moment magnétique avec Bxz10 aboutit à un couple ayant la même direction que précédemment mais on comprend qu'en outre l'interaction de Mb avec la composante en tangage By du champ magnétique terrestre aboutit aussi à un couple ayant cette direction. Plus généralement l'interaction de Mb avec B est un couple Cmxz15 donné par: Cmxz = MbAB Un cas particulier de cette solution générale est le cas précité o le moment magnétique est parallèle à Y (il suffit alors d'une seule bobine) Un autre cas particulier intéressant est celui o l'on choisit de donner à Mb une orientation, dans le plan (Y, B), perpendiculaire à B (l'efficacité est alors maximale); le module de Cmxz est alors donné par le produit des modules de Mb et de B. Un autre cas consiste par exemple à fixer Mb dans le plan roulis-lacet, c'est-à-dire parallèle à Bxz (seules des bobines (ou équivalents) en roulis-lacet sont alors requises). Soit a l'angle que fait avec l'axe de roulis X la direction instantanée d définie par le produit vectoriel du vecteur unitaire Y de l'axe de tangage avec le champ magnéti- que terrestre local; c'est aussi l'angle que fait Bxz avec l'axe de lacet Z. Il est important de noter que l'orientation de Bxz détermine l'orientation de Cmxz; si Bx et Bz sont les composantes en roulis et en lacet (axes X et Z du satellite) du champ magnétique terrestre local et si on compte positive- ment a autour de -Y, on a5 a = Arc tg (-Bx/Bz) Ainsi que cela apparaîtra plus loin, l'invention prend en considération, dans les équations de contrôle d'attitude au moins en roulis-lacet, le fait que tous les couples réalisables dans le plan roulis-lacet à un instant10 donné par interaction de Mb avec le champ magnétique terres- tre sont nécessairement dans la direction a définie ci- dessus par l'angle a (les couples peuvent bien sûr être négatifs et donc être en sens inverse de c). La deuxième équation de la dynamique, appliquée au système constitué par le satellite pris dans son environ- nement permet d'écrire (en écriture vectorielle) en tenant compte du fait que les efforts extérieurs sont la somme d'un couple perturbateur Cp (en principe inconnu) et d'un couple de commande Cm (à déterminer): I + H + AIQ + QAH = Cp + Cm o Q est le vecteur rotation instantanée du satellite dans son repère, I est la matrice d'inertie du satellite, H est le moment cinétique (éventuel) embarqué au sein du satellite; il s'agit par exemple d'une roue cinétique d'orientation fixe, orientée de préférence selon l'axe de tangage. On sait que l'équation dynamique vectorielle précitée à trois composantes peut se décomposer en une équation dynamique matricielle à deux composantes dans le plan roulis-lacet (correspondant à un système de deux équations scalaires) et une équation en tangage. On s'inté- resse dans ce qui suit A l'équation dans le plan roulis- lacet; celle-ci peut s'écrire, en linéarisant au premier ordre les petits angles, sous la forme: AX+ BX+ DX = Cpxz + Cmxz ou: Cpxz est la composante du couple perturbateur dans le plan roulis-lacet, Cmxz est la composante du couple de contrôle dans le plan roulis-lacet (l'indice "m" rappelant ici que, comme indiqué ci-dessous, ce couple est d'origine magnétique), X est le vecteur colonne formé des angles t et $ de roulis et de lacet, A, B et D sont des matrices de coefficients caractérisant le satellite, notamment en inertie. Compte tenu de ce que l'on connaît l'orientation dans laquelle on exerce le couple Cmxz de contrôle on peut écrire ce couple Cmxz en projection dans le plan roulis-lacet sous la forme: Cmxz = C. [sin " o C est un scalaire (correspondant au module du vecteur Cmxz de la figure 1). On peut choisir a priori la forme de la loi de commande donnant Cmxz en fonction de X, ainsi que de sa dérivée X et/ou de son intégrale fx. Une commande de type proportionnel/dérivée revient A écrire: C = Kc.X + Lc.X o Kc est une matrice ligne de gains k, et k2 et o Lc est une matrice ligne de gains Ql et e2, l'indice 1 correspondant aux gains en roulis et l'indice 2 correspondant aux gains en lacet. Si l'on choisissait une loi de commande de type proportionnel/intégrale, on écrirait: C = Kc.X + Nc.fX oû les matrices Kc et Nc seraient de même deux matrices lignes de gains k1, k2, n1, n2. Si on choisit une loi de commande de type PID (proportionnel/intégrale/dérivée) on ferait intervenir six coefficients k1, k2, Ql, 2, n1 et n2. Bien entendu la loi de commande est d'autant plus simple que le nombre de coefficients ou gains est faible, en pratique égal à 4. Dans le cas de la loi de commande proportionnel/- dérivée on peut écrire l'équation dynamique précitée sous la forme: A X + BX + D'X = Cpxz en intégrant dans B' et D' les coefficients de X et de X dans l'expression du couple Cmxz. Si l'on choisit une loi de commande faisant intervenir l'intégrale de X, on aurait dans le membre gauche de l'équation dynamique un quatrième terme par exemple noté F. Cela explique que, pour des raisons de simplicité de l'équation, on préfère utiliser ici une loi de type propor- tionnel/dérivée. Selon un aspect important de l'invention, on fait varier les divers gains au cours du temps selon un critère prédéterminé approprié, tenant compte des valeurs instanta- nées des composantes en roulis-lacet du champ magnétique terrestre, ou de a et du module de Bxz défini précédemment. C'est pourquoi on peut dire que la loi de commande utilisée est une loi adaptative. De manière préférée, l'élaboration des gains se fait à partir de choix sur les pôles de l'équation dynamique du système, c'est-à-dire sur les racines de son équation caractéristique. Ainsi qu'on le sait, l'équation caractéristique du système est un polynôme dont l'ordre n dépend de l'ordre de l'équation différentielle considérée (égal à 2 si la loi de commande ne comporte que des termes proportionnels dérivés, égal à 3 si la loi de commande comporte un terme10 intégral) et du nombre de composantes du vecteur X considéré (2 si l'on ne s'intéresse qu'au système dans le plan roulis- lacet). En fait n est le produit des valeurs précitées. Il est à noter que le nombre n, c'est-à-dire le nombre de pôles ou de racines, peut être différent du nombre de gains de la loi de commande (cas d'une loi de commande de type propor- tionnel/intégral ou dévrivée/intégral). Ces pôles sont généralement des valeurs complexes, mais il existe toutefois des grandeurs intermédiaires qui ont l'avantage d'être toujours réelles et d'avoir une signification physique. Plus précisément, au moins dans le cas de la loi de commande de type proportionnel/dérivée appliqué au contrôle roulis-lacet, on peut écrire le polynôme caractéris- tique sous la forme: P= (x2 + 2(C olx +)2) (x2 + 2C2 2X + 2) Les grandeurs intermédiaires sont des pulsations W1 et w2 et des amortissements fj et 2' Dans le cas o il y a six pôles, on utilise des grandeurs supplémentaires w3 et 3. Les pulsations w1 et w2 seront choisies d'autant plus grandes que l'on veut réagir vite (commande "nerveuse"), tandis que les amortissements f 1 et 2 seront choisis d'autant plus grands que l'on veut amortir rapidement les mouvements. En particulier, on peut donner à ces grandeurs des valeurs constantes au cours du temps. En variante, on peut les faire évoluer en fonction de la position du satellite sur son orbite. Pour assurer une bonne stabilité du système lors- que a varie, les pulsations wl et w2 sont de préférence choisies très supérieures à la pulsation orbitale o (par exemple 10 à 100 fois plus grands). Cela conduit à un bon découplage en fréquence de l'évolution des gains de commande d'une part et de l'évolution de l'attitude du satellite d'autre part. A partir de la donnée de ces pulsations et amortissements, donc des pôles, on peut déterminer les coefficients du polynôme caractéristique et en déduire les valeurs instantanées à donner aux gains kl, k2, e1 et Q2 (voire les autres gains s'il y en a plus que 4), ce qui donne la valeur instantanée C à réaliser dans la direction a par la relation C = kl.q + 11.4 + k2.- + 12.4 puis, en fonction de la direction du plan (Y, B) dans laquelle on a choisi de générer le moment magnétique, en déduire l'amplitude de ce dernier. Ainsi que cela est schématisé à la figure 2, le contrôle d'attitude en roulis-lacet comporte donc les étapes suivantes, avec une loi de commande proportionnel/dérivé: - mesure en i des valeurs instantanées de Bx et de Bz (ou directement de Bxz) à l'aide de magnétomètres de tout type connu approprié, et identification de a, détermination en 2, à partir de la connaissance des caractéristiques dynamiques du satellite (A, B, D) et de la donnée des caractéristiques imposées au contrôle (contraintes imposées à f, f2, 1, w2 ou aux racines proprement dites du polynôme caractéristique P), des valeurs instantanées des gains de kl, k2, e1, Q2, c'est-à- dire des matrices lignes de la loi de commande, - détermination en 2 également, à partir de la mesure en 3 par tout système connu approprié de détection d'attitude en roulis-lacet (détecteurs terrestre, solaire voire stellaire et/ou gyroscopes et/ou filtres estimateurs voire réception depuis une source extérieure) des angles et des vitesses à corriger en roulis et en lacet (c'est- à-dire à partir de la connaissance de X et X) du couple Cmxz à obtenir en roulis-lacet (couple de consigne), détermination en 4 du moment magnétique My à générer selon Y, ou plus généralement de Mb, - et réalisation en 4 de ce moment qui, par interaction en avec le champ magnétique terrestre B, génère un couple réel, d'ou une réaction en attitude du satellite détectée en 3 et venant ainsi fermer la boucle du système, via la dynamique du satellite 6. Ce cycle d'opérations avec observation des valeurs instantanées de Bxz, et des nouvelles valeurs d'attitude compte tenu de la réaction du satellite au moment My généré et interagissant avec le champ magnétique terrestre assure le contr6le de l'attitude du satellite. La détermination de la projection du champ magnétique terrestre instantané (ou courant) dans le plan roulis-lacet peut être faite par mesure (voir ci-dessus) mais aussi à l'aide d'un modèle mathématique ou par combinaison de mesures et de ce modèle. Il est rappelé que la direction de cette projection est variable et fait un tour par orbite. Il en est de même pour la direction d dans laquelle il est possible de générer un couple magnétique purement dans le plan roulis-lacet. Il est à noter que cette direction est unique dans le cas des orbites inclinées par rapport à l'équateur magnétique, et que dans le cas d'une orbite non inclinée par rapport à cet équateur on peut réaliser un couple purement dans le plan roulis-lacet dans n'importe quelle direction de ce plan, l'invention n'est donc pas limitée à ce cas particu- lier. Pour déterminer les gains de contrôle en roulis, en lacet, en vitesse de roulis et en vitesse de lacet (dans le cas d'une loi de contrôle de type proportionnel/dérivée, mais ce qui suit se généralise aisément à toute autre loi) on identifie les coefficients du polynôme caractéristique du système contrôlé, calculés avec les quatre gains de retour comme paramètres, avec les coefficients du polynôme caracté-10 ristique désiré, correspondant aux pôles choisis, puis on résoud le système de quatre équations ainsi obtenues, avec pour inconnues les quatre gains de contrôle. En pratique les gains peuvent être calculés à bord par interpolation de plusieurs valeurs de gains stockées à bord, précalculées au sol pour plusieurs valeurs de a. Les gains peuvent aussi être calculés au sol et transmis au satellite en vol. On peut également calculer au sol lesprincipaux coefficients de FOURIER de chacun des quatre gains en fonction de a et les stocker à bord, dans ce cas, on peut alors reconstituer chaque gain à bord en fonction de l'angle a courant et des coefficients stockés (en pratique les cinq premiers coefficients suffisent). Pour ce qui est du contrôle complémentaire en tangage, celui-ci est par exemple réalisé à l'aide d'une roue cinétique en tangage (cas envisagé ci-dessus) dont la vitesse est modifiée pour compenser les variations d'attitude en tangage. En temps utile, voire en permanence, la désatura- tion de cette roue peut être réalisée à l'aide de bobines magnétiques agissant dans le plan roulis-lacet (par exemple une bobine en roulis et une bobine en lacet) en sorte de générer, par interaction avec la composante en roulis-lacet de champ magnétique terrestre, le couple Cy de désaturation souhaité. En variante, on peut laisser la roue d'inertie (si elle existe, ce qui est de préférence le cas pour assurer une stabilisation gyroscopique de l'axe de tangage) tourner à une vitesse constante et générer les couples de correction en tangage à l'aide des bobines magnétiques précitées disposées dans le plan roulis-lacet. Pour assurer une efficacité maximale de ce contrôle ou de cette désaturation en tangage les bobines magnétiques du plan roulis-lacet sont de préférence comman- dées en sorte de générer un moment magnétique qui est10 perpendiculaire à Bxz; il existe toutefois une possibilité d'interaction avec le champ magnétique en tangage qui provoque un couple perturbateur dans le plan roulis-lacet. Le contrôle ou la désaturation en tangage perturbe donc le roulis-lacet, mais le contrôle du roulis-lacet ne perturbe pas le tangage, ainsi il n'y a pas de possibilité de rebou- clage des perturbations susceptible de déstabiliser le système complet (roulis-lacet et tangage). La loi de commande pour le contrôle en tangage peut être de tout type (proportionnel/dérivée, intégral/déri- vée ou proportionnel/intégral/dérivée) avec des gains éventuellement constants; toutefois, de manière préférée, les gains sont variables en sorte de respecter une contrainte sur les pôles (ou les pulsation et amortissement), par exemple une contrainte de constance. Des résultats de simulation sont présentés aux figures 3 à 6 qui correspondent aux données numériques suivantes: -altitude 900 km inclinaison de l'orbite 650 - w 21/6000 - période orbitale 6000 secondes w et w2 0,004 et 0,01 rd/sec 1 et 2 1 et 1,1 - moment cinétique en Y 30 N.m.s - désaturation du moment cinétique en tangage par des bobines en X et Z - matrice d'inertie diagonale du satellite - éventuels couples perturbateurs selon figure 7 - modèle de champ magnétique d'ordre 8 (8 harmo niques sphériques) - conditions initiales: = 0,5 = 10 0o = 1,5 La comparaison des figures 3 et 4 (sans couples perturbateurs) aux figures 5 et 6 (avec couples perturba- teurs) montre la capacité du procédé de contrôle à supporter des couples perturbateurs importants. Vis à vis des solutions antérieures, l'invention élargit considérablement le champ d'application des systèmes de contrôle d'attitude magnétique: aux orbites ni quasi-polaires, ni quasi- équatoriales, - aux satellites naturellement très instables, en particu- lier à cause du gradient de gravité (par exemple à cause d'un terme croisé roulis-lacet important dans la matrice d'inertie). Ainsi, l'invention permet: - d'éviter l'adjonction de roues cinétiques ou à réaction dans le plan roulis-lacet, - de se passer d'un système propulsif pénalisant en terme de masse et limitant la durée de vie, - de se dispenser d'une stabilisation passive par gradient de gravité, également pénalisante en masse, moins précise, et difficilement applicable aux moyens et gros satellites, - d'éviter les interruptions de commande deux quarts d'orbite sur 4 qui laissent dériver l'attitude du satellite 50 % du temps, dégradant ainsi le pointage, - de ne demander que des calculs simples (permet de se dispenser d'un modèle embarqué de champ magnétique), - de s'adapter au champ magnétique local réel, conférant ainsi une grande robustesse, - de maîtriser les caractéristiques dynamiques du système et ainsi garantir une excellente stabilité et une bonne précision de pointage malgré des perturbations importan- tes. On peut résumer que le système représentant la dynamique des axes roulis et lacet d'un satellite en orbite inclinée par rapport à l'équateur magnétique et un générateur de moment magnétique capable de produire un couple de contrôle par interaction avec le champ magnétique, est instationnaire du fait de l'évolution du champ magnétique en axes satellites lorsque ce dernier parcourt son orbite. L'invention permet de contrôler ledit système instationnaire en fixant avantageusement ses pôles en tout point de l'orbite, donc en lui donnant des caractéristiques dynamiques - telles que la rapidité ou l'amortissement - à volonté en tout point de l'orbite (dans la limite correspon- dant à la saturation du générateur de moment magnétique, bien sûr). En particulier, on pourra choisir des pôles constants tout au long de l'orbite. Les pulsations correspondant aux pôles du système pourront être choisies (non restrictif) au moins quatre fois plus grandes que la pulsation orbitale, et les amortissements correspondants choisis de l'ordre de l'unité. REVENDICATIONS 1. Procédé de contrôle de l'attitude d'un satellite devant être stabilisé selon des axes de roulis, de tangage et de lacet, selon lequel, pour le contrôle en5 roulis-lacet: - on détermine l'orientation instantanée dans le plan roulis-lacet de la direction instantanée c définie par le produit vectoriel unitaire Y de l'axe de tangage avec le champ magnétique terrestre local B, - on détermine des grandeurs d'attitude en roulis-lacet du satellite, - on détermine l'amplitude et le sens d'un couple de correction Cmxz à appliquer au satellite selon ladite direction, cette amplitude étant déterminée à partir desdites grandeurs d'attitude en roulis-lacet selon une loi de commande prédéterminée dont les gains sont élaborés de façon adaptative en sorte de respecter des conditions prédétermi- nées sur les pôles du système modélisant la dynamique du satellite dans le plan roulis-lacet, - on génère dans un plan perpendiculaire à cette direction instantanée un moment magnétique d'amplitude propre à générer selon ladite direction ledit couple de correction Cmxz par interaction avec le champ magnétique terrestre, grâce à quoi ce couple de correction du roulis-lacet n'a aucune composante en tangage. 2. Procédé selon la revendication 1, caractérisé en ce que les conditions prédéterminées sur les pôles du système sont des conditions de sensible constance au cours du temps. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que lesdits gains sont élaborés par interpolation à partir de données stockées à bord et d'au moins l'orientation de la composante Bxz du champ magnétique terrestre instantané dans le plan roulis-lacet. 4. Procédé selon la revendication 1 ou 2, caractérisé en ce que lesdits gains sont élaborés à partir de coefficients de FOURIER stockés à bord et de la valeur instantanée de l'orientation de la composante Bxz du champ magnétique terrestre local dans le plan roulis-lacet. 5. Procédé selon l'une quelconque des revendica- tions 1 à 4, caractérisé en ce que la loi de commande est du type proportionnel/dérivée, les grandeurs d'attitude sont les angles et les vitesses angulaires en roulis-lacet, et les gains sont au nombre de quatre. 6. Procédé selon l'une quelconque des revendica- tions 1 à 4, caractérisé en ce que la loi de commande est du type proportionnel/intégral, les grandeurs d'attitude sont les angles en roulis-lacet et leurs intégrales et les gains sont au nombre de quatre. 7. Procédé selon l'une quelconque des revendica- tions 1 à 4, caractérisé en ce que la loi de commande est du type proportionnel/intégral/dérivée, les grandeurs d'attitude sont les angles en roulis-lacet et leurs vitesses et leurs intégrales, et les gains sont au nombre de six. 8. Procédé selon l'une quelconque des revendica- tions 1 à 7, caractérisé en ce que les pôles auxquels sont appliquées les conditions prédéterminées sont caractérisés par les grandeurs f et wi telles que le polynôme caractéris- tique du système puisse s'écrire sous la forme d'un produit de n termes: X2 x + 2i iX + 2i ou i est compris entre 1 (inclus) et n (inclus). 9. Procédé selon la revendication 8, caractérisé en ce que les pôles sont choisis en sorte que les grandeurs wi sont supérieures d'au moins un ordre de grandeur à la pulsation orbitale. 10. Procédé selon l'une quelconque des revendica- tions 1 à 9, caractérisé en ce qu'on détermine l'amplitude de la composante Bxz et on détermine l'amplitude du moment magnétique, noté My, par division de l'amplitude du couple de5 correction par cette amplitude de la composante Bxz, et on génère ce moment magnétique My selon l'axe de tangage. 11. Procédé selon l'une quelconque des revendica- tions 1 à 9, caractérisé en ce qu'on détermine l'amplitude IBI du champ magnétique terrestre B, et on génère le moment magnétique, noté Mb, perpendiculairement au champ magnétique terrestre, tel que: BAYEnxz Mb ACmxz IB12 12. Procédé selon l'une quelconque des revendica- tions 1 à 11, caractérisé en ce qu'en outre on contrôle l'attitude en tangage par modification de la vitesse de rotation d'une roue d'inertie en tangage. 13. Procédé selon la revendication 12, caracté- risé en ce qu'on génère des couples en tangage de désatura- tion de la roue d'inertie par activation d'au moins un générateur de moment magnétique générant un moment magnétique dans le plan roulis-lacet. 14. Procédé selon l'une quelconque des revendica- tions 1 à 11, caractérisé en ce qu'on contrôle l'attitude en tangage par génération dans le plan roulis-lacet d'un moment magnétique d'amplitude déterminée en fonction des mesures d'attitude en tangage. 15. Procédé selon la revendication 13 ou la revendication 14, caractérisé en ce qu'on génère ce moment magnétique dans le plan roulis- lacet selon ladite direction instantanée. 16. Procédé selon l'une quelconque des revendica- tions 1 à 15, caractérisé en ce que l'on détermine les grandeurs d'attitude par détection à l'aide de détecteurs terrestre, solaire et/ou stellaire, de gyroscopes ou de gyromètres, et/ou estimation à l'aide de filtres estimateurs. 17. Procédé selon l'une quelconque des revendications 1 à 15, caractérisé en ce qu'on détermine certaines au moins des grandeurs d'attitude par comparaison de mesures, selon un ou plusieurs axes, du champ magnétique terrestre instantané avec un modèle de ce champ magnétique terrestre. 18. Procédé selon l'une quelconque des revendica- tions 1 à 17, caractérisé en ce qu'on détermine le champ magnétique terrestre au moins en roulis-lacet par mesure à l'aide d'au moins un magnétomètre à un ou plusieurs axes de mesure. 19. Procédé selon l'une quelconque des revendica- tions 1 à 17, caractérisé en ce qu'on détermine au moins l'orientation du champ magnétique terrestre en roulis-lacet à partir de signaux reçus du sol. 20. Procédé selon l'une quelconque des revendica- tions 1 à 17, caractérisé en ce qu'on détermine au moins l'orientation du champ magnétique terrestre local en roulis- lacet à partir d'un modèle de ce champ. 21. Procédé selon l'une quelconque des revendica- tions 1 à 20, caractérisé en ce que le satellite suit une orbite d'inclinaison comprise entre 15 et 75 environ. 22. Satellite stabilisé en roulis, en tangage et en lacet, comportant un dispositif de contrôle d'attitude en roulis-lacet comportant: - un premier dispositif de détermination d'au moins l'orientation du champ magnétique terrestre local instantané en roulis-lacet Bxz, - un deuxième dispositif de détermination d'attitude pour fournir des signaux d'attitude au moins en roulis et en lacet, - un calculateur connecté aux premier et second dispositifs, contenant une loi de commande adaptative à gains variables donnant une valeur de couple de correction en roulis/lacet à appliquer perpendiculairement à l'axe de tangage et à l'orientation instantanée de Bxz à partir des signaux d'attitude, et une loi donnant lesdits gains varia- bles en fonction d'au moins l'orientation de Bxz, - un générateur de moment magnétique adapté à générer dans le plan de l'axe de tangage et de Bxz un moment magnétique propre à générer par interaction avec le champ magnétique terrestre ladite valeur de couple à appliquer. 23. Satellite stabilisé selon la revendication 22, caractérisé en ce que le satellite comporte en outre un système de contrôle d'attitude en tangage comportant une roue d'inertie en tangage à vitesse variable. 24. Satellite stabilisé selon la revendication 23, caractérisé en ce que le système de contrôle d'attitude en tangage comporte un générateur de moment magnétique en roulis-lacet adapté à générer avec la composante en roulis- lacet du champ magnétique terrestre un couple de désaturation en tangage. 25. Satellite stabilisé selon la revendication 22, caractérisé en ce que le dispositif de contrôle d'atti- tude en roulis-lacet fait partie d'un dispositif de contrôle d'attitude en roulis, tangage et lacet, le calculateur contenant en outre une loi de contrôle d'attitude en tangage déterminant un couple de contrôle d'attitude en tangage, et le générateur de moment magnétique étant adapté à générer un moment magnétique engendrant, par interaction avec le champ magnétique terrestre, le couple de correction en roulis-lacet et le couple de correction en tangage.