Brevet FR2789172A1 — publié le 20000804

APPAREIL À GYROMTRES ET ACCELEROMETRES POUR LA DETERMINATION DES A'TTITUDES D'UN AERODYNE L'invention est relative à un appareil de détermina- tion des attitudes d'un engin volant. Pour piloter un aérodyne, il est nécessaire de connaître ses attitudes. Les attitudes sont par définition les trois angles d'Euler y, 0, A, ou, respectivement, cap, assiette, et incli- naison, qui définissent l'orientation du repère lié à l'aérodyne par rapport au repère géographique local. Le repère (ou trièdre) géographique local est un trièdre orthonormé dont un axe est suivant la verticale locale (vers le bas) et les deux autres axes suivant les directions nord et est locales. Le repère lié à l'aérodyne est un trièdre orthonormé dont l'axe X, ou axe de roulis, est colinéaire au fuselage de l'aérodyne, l'axe Y, ou axe de tangage, est perpendiculaire à son plan de symétrie, et l'axe Z, ou axe de lacet, est inclus dans le plan de symétrie. Les attitudes sont déterminées à l'aide d'un disposi- tif gyrométrique dit à composants liés, constitué d'au moins trois gyromètres solidaires entre eux et solidaires de la struc- ture de l'aéronef. Chacun de ceux-ci fournit une composante, suivant son axe, du vecteur Q de rotation instantanée de l'aé- rodyne par rapport à un référentiel inertiel. Ce dispositif per- met de mesurer les trois composantes de Q dans le repère X, Y, Z lié à l'aérodyne. Ces composantes sont, respectivement, le roulis p suivant l'axe X, le tangage q suivant l'axe Y et le lacet r suivant l'axe Z. Ces trois composantes sont intégrées, par exemple par la méthode dite des quaternions. L'intégration fournit des qua- ternions d'attitudes, d'o l'on déduit, avec l'orientation ini- tiale de l'aérodyne, une matrice d'attitudes puis les angles d'Euler. Les signaux de mesure délivrés par les gyromètres pré- sentent des défauts appelés dérives. Comme ces signaux sont intégrés, les erreurs sur les attitudes augmentent avec le temps. C'est pourquoi, habituellement, pour déterminer les atti- tudes, on fait appel, en plus des gyromètres, à des accéléro- mètres qui sont utilisés pour corriger les résultats fournis par l'intégration des signaux des gyromètres. Cette correction s'ef- fectue dans une boucle d'asservissement. Les accéléromètres solidaires du dispositif gyromé- trique à composants liés peuvent fournir, grâce à la pesanteur, les attitudes en assiette et inclinaison de l'aérodyne; toute- fois les signaux qu'ils délivrent sont entachés d'erreurs qui rendent les résultats difficilement exploitables tels quels. Par contre, la combinaison, dans une boucle d'asservissement, des signaux fournis par des gyromètres et des accéléromètres permet de compenser les dérives des gyromètres tout en conservant l'a- vantage de ces derniers qui est de fournir des résultats avec un bruit relativement faible à court terme. Dans des appareils connus, la comparaison entre les signaux gyrométriques et accélérométriques constitue une opé- ration complexe car les signaux des accéléromètres sont intégrés de façon à représenter des vitesses et ces dernières sont com- parées à des références; cette comparaison sert à élaborer la correction du dispositif gyrométrique. L'intégration fait appel aux attitudes issues du dispositif gyrométrique. L'invention fournit un appareil dans lequel les cor- rections sont effectuées de façon simple et fiable. L'appareil selon l'invention est caractérisé en ce que la comparaison est effectuée entre, d'une part, un vecteur de direction verticale fourni par le dispositif gyrométrique et, d'autre part, un vecteur de direction verticale fourni par les accéléromètres. De préférence, la comparaison consiste à effec- tuer le produit vectoriel de ces deux vecteurs. Le vecteur de direction verticale fourni par les accé- léromètres est colinéaire à l'accélération de la pesanteur, durant les phases de vol rectiligne non accéléré. Ainsi, il n'est pas nécessaire d'effectuer une intégration des signaux accélérométriques. Par ailleurs, dans le dispositif gyromé- trique, le vecteur de direction verticale constitue la dernière colonne de la matrice d'attitudes obtenue avec ce dispositif. En outre, les deux vecteurs de direction verticale peuvent être exprimés dans un même repère lié à l'aérodyne. On rappelle ici que la matrice d'attitudes est une matrice 3x3 dans laquelle la première colonne représente la direction du nord dans le repère XYZ de l'aérodyne, et les seconde et troisième colonnes représentent respectivement les directions est et verticale, également en repère XYZ. Il est d'ailleurs à noter que le produit vectoriel est isotrope, c'est-à-dire indépendant du repère. Quand on fait appel à un produit vectoriel, les cor- rections à apporter sont aussi d'une grande simplicité car le module du produit vectoriel est proportionnel au sinus de l'an- gle de la rotation qui fait coïncider le vecteur de direction verticale fourni par le dispositif gyrométrique avec le vecteur de direction verticale fourni par le dispositif accéléromé- trique. De plus, la direction du produit vectoriel donne l'axe de cette rotation. Si la comparaison consiste à effectuer un produit vec- toriel, des dispositions supplémentaires peuvent être néces- saires. En effet, un produit vectoriel étant, par définition, perpendiculaire à chacun des vecteurs du produit, sa composante est nulle suivant la direction verticale, et si l'on veut esti- mer la composante verticale de la dérive du dispositif gyromé- trique, il est nécessaire de disposer de données supplémen- taires. Dans un mode de réalisation de l'invention, pour obte- nir une estimation de la composante verticale de la dérive du dispositif gyrométrique, on effectue une projection orthogonale du vecteur! de composantes p, q r, qui est fourni par les gyromètres, après compensation des dérives (par l'asservissement mentionné ci-dessus), sur la direction verticale fournie par le dispositif gyrométrique. Pour obtenir cette projection, il suffit d'effectuer le produit scalaire du vecteur O par le vecteur unitaire Ug de la direction verticale fournie par le dispositif gyrométrique. Ce vecteur projeté peut être utilisé, dans un asservissement à retour intégral, pour compenser la composante verticale de la dérive du vecteur Q fourni par le dispositif gyrométrique. Le produit scalaire O.U g représente la dérive du gyromètre en direction verticale seulement si l'avion est en phase de vol rectiligne, c'est-à-dire qu'il ne tourne pas autour de la direction verticale. Ainsi, dans un mode de réalisation, on détecte les phases de virage et on ouvre la boucle d'as- servissement de lacet (permettant de compenser la composante verticale des dérives du dispositif) pendant les phases de virage. De même, le dispositif accélérométrique ne fournit un vecteur représentant l'accélération de la pesanteur que pendant les phases de vol rectiligne non accéléré. Il est donc pré- férable d'ouvrir au moins certaines boucles d'asservissement du dispositif gyrométrique pendant les phases de virage et/ou d'ac- célération. On ouvre la boucle d'asservissement de roulis dans les phases de virage et on ouvre la boucle d'asservissement de tangage pendant les phases d'accélération. Lors des virages accélérés, on ouvre de préférence la boucle de roulis, la boucle de tangage restant fermée. Dans une variante, pour obtenir une estimation de la composante verticale de la dérive du dispositif gyrométrique, au lieu d'utiliser un asservissement en lacet, on fait appel à un asservissement en cap, c'est-à-dire qu'on asservit la direction du nord exprimée en repère aérodyne, fourni par le dispositif gyrométrique, sur le nord magnétique fourni par un capteur externe tel qu'un magnétomètre. De préférence, l'écart de cap est représenté par un vecteur dont la direction correspond à la verticale du dispo- sitif gyrométrique de façon qu'il n'entre pas en conflit avec le produit vectoriel mentionné ci-dessus qui, par définition, est perpendiculaire à la direction verticale fournie par le dispo- sitif gyrométrique. A cet effet, on peut utiliser le produit vectoriel d'un vecteur Um représentant la direction du nord magnétique (fourni par un magnétomètre) par un vecteur bl représentant la direction du nord fourni par le dispositif gyro- métrique dans la première colonne de la matrice d'attitudes. Ce produit vectoriel est alors projeté sur la direction verticale b3 fournie par le dispositif gyrométrique. L'invention prévoit ainsi un appareil de détermination d'attitudes, notamment l'assiette et l'inclinaison d'un aéro- dyne, comportant: un dispositif gyrométrique à composants liés fournis- sant les composantes du vecteur de rotation instantanée de l'aé- rodyne dans un repère lié à l'aérodyne, des moyens pour calculer, à partir des mesures gyromé- triques, une matrice d'attitudes définissant l'orientation de l'aérodyne par rapport au repère géographique local, des accéléromètres solidaires du dispositif gyromé- trique, et des moyens pour comparer des données calculées à par- tir des accéléromètres à des données calculées à partir des gyromètres afin de fournir des corrections pour compenser les erreurs ou dérives du dispositif gyrométrique; cet appareil est caractérisé en ce que la comparaison est effectuée entre un vecteur Ug de direction verticale fourni par le dispositif gyrométrique et un vecteur Ua de direction verticale fourni par les accéléromètres. L'appareil comporte, dans une réalisation, des moyens pour calculer, à partir des mesures accélérométriques, la direc- tion de la verticale dans un repère lié à l'aérodyne. Les vec- teurs de direction verticale sont, par exemple, définis par leurs coordonnées dans un repère lié à l'aérodyne. Le vecteur de direction verticale fourni par le dispo- sitif gyrométrique peut être obtenu à partir de la dernière colonne de la matrice d'attitudes. De préférence, les vecteurs de direction verticale ont une longueur unitaire. Dans une réalisation, les moyens de comparaison effec- tuent le produit vectoriel des deux vecteurs de direction ver- ticale. L'appareil comporte une boucle d'asservissement de type proportionnel qui, de préférence, présente un gain limité à une valeur maximale prédéterminée. Dans un mode de réalisation, l'appareil comprend une boucle d'asservissement à intégrateur pour compenser les dérives. Quand on effectue le produit vectoriel des deux vec- teurs de direction verticale, il est avantageux que l'appareil comporte un moyen pour effectuer une correction suivant la direction verticale. Dans ce cas, ce dernier moyen comprend, par exemple, un moyen de calcul de la projection, sur la direction verticale, du vecteur de rotation instantanée fourni par le dispositif gyrométrique, la direction verticale sur laquelle est effectuée la projection étant fournie par le dispositif gyro- métrique, et le vecteur de rotation instantanée qui est projeté étant corrigé par la boucle d'asservissement. Selon un mode de réalisation, le vecteur de direction verticale est fourni, d'une part, par les accéléromètres et, d'autre part, à partir de la vitesse de l'aérodyne délivrée, par exemple, par des moyens anémométriques. Dans ce dernier cas, la direction verticale est, par exemple, déterminée à partir de l'équation suivante: Fd dV + AV- FC m dt m Fd étant la force de la pesanteur, de direction ver- ticale, V le vecteur de vitesse de l'aérodyne déterminé par les moyens anémométriques, Q le vecteur de rotation instantanée de l'aérodyne déterminé par le dispositif gyrométrique, et Fc la résultante des forces de contact déterminée par les accélé- romètres, et m la masse de l'aérodyne. Dans une réalisation, l'appareil comporte un moyen d'asservissement du cap gyrométrique sur le cap magnétique. Dans ce dernier cas, et lorsqu'on effectue le produit vectoriel des vecteurs de direction verticale, l'asservissement du cap peut être déterminé à partir du vecteur de correction satisfaisant à l'équation suivante: Fc = [(Um A bllb3lb3 = (m.b2b3 équation dans laquelle Um est un vecteur représentant la direc- tion du champ magnétique terrestre, de préférence d'amplitude unité, et les vecteurs b1, b2 et b3 sont trois vecteurs dont les coordonnées sont fournies par les colonnes de la matrice d'attitudes délivrée par le dispositif gyrométrique. Selon une réalisation, l'appareil comprend des boucles fermées d'asservissement des attitudes utilisant les signaux de correction fournis par les moyens de comparaison, et des moyens pour ouvrir au moins une boucle en dehors des phases de vol rec- tiligne et à vitesse constante de l'aérodyne. Dans ce dernier cas, l'ouverture des boucles d'asser- vissement peut intervenir quand l'écart entre le vecteur de direction verticale fourni par le dispositif gyrométrique et le vecteur de direction verticale fourni par les accéléromètres dépasse une valeur prédéterminée. En cas de virage, par exemple, une boucle d'asservissement en roulis reste ouverte, alors qu'une boucle d'asservissement correspondante en tangage reste fermée. Le taux de virage peut être utilisé pour distinguer entre des phases de vol rectiligne et des phases de vol dif- férentes et une boucle d'asservissement en roulis peut être fer- mée quand le taux de virage est inférieur en valeur absolue à un seuil déterminé. Quand on prévoit une boucle pour compenser les dérives, on peut aussi prévoir une boucle d'asservissement pour assurer une correction des mesures gyrométriques suivant la direction verticale et cette boucle d'asservissement comprend un intégrateur qui, de préférence, est commun à la boucle compen- sant les dérives. Dans une réalisation de cette dernière disposition, l'intégrateur comporte des moyens permettant de limiter la pente moyenne de son signal de sortie et des moyens pour que, à chaque instant, le signal de sortie présente un écart inférieur à un seuil par rapport à un signal à variation linéaire ayant une pente déterminée. D'autres caractéristiques et avantages de l'invention apparaîtront avec la description de certains de ses modes de réalisation, celle-ci étant effectuée en se référant aux dessins ciannexés sur lesquels: la figure 1 est un schéma d'un boîtier, avec ses axes, d'un appareil selon l'invention, la figure 2 est un schéma de boucle d'asservissement d'un appareil conforme à l'invention, la figure 3 est un schéma analogue à celui de la figure 2, mais pour une variante, les figures 4, 5 et 6 sont des diagrammes illustrant des modes de fonctionnement, conformes à l'invention, des cir- cuits d'asservissement représentés sur la figure 2 ou la figure 3, la figure 7 est un schéma d'un mode de réalisation d'intégrateur utilisé dans le mode de réalisation représenté sur la figure 2 ou la figure 3, la figure 8 est un diagramme donnant, à titre d'exem- ple, un mode de fonctionnement de l'intégrateur représenté sur la figure 7, et la figure 9 correspond à une variante de la figure 3. L'exemple de réalisation de l'invention que l'on va décrire en relation avec les figures se rapporte à un instrument ou appareil pour déterminer les attitudes d'un avion qui com- prend au moins trois gyromètres (non représentés) dont les com- posants ont une position fixe par rapport à l'avion. Le dispo- sitif gyrométrique est, dans ces conditions, dit à "composants liés". Autrement dit, les gyromètres sont des capteurs dont les axes ont des positions bien définies par rapport aux axes de l'avion. L'appareil comprend aussi des accéléromètres (non mon- trés) pour corriger les erreurs ou dérives des gyromètres. On sait, en effet, que les gyromètres fournissent des données (les attitudes d'un avion) à faibles bruits mais présentant des dérives, ou biais, entraînant des erreurs croissant avec le temps. Pour compenser les dérives, on fait appel à des accé- léromètres qui donnent les composantes de l'accélération de l'aérodyne correspondant à la résultante des forces de contact s'appliquant à cet aérodyne. Par "forces de contact", on entend la poussée et les forces aérodynamiques (traînée et portance). A partir de ces composantes fournies par les accéléromètres, on peut en déduire l'accélération de la pesanteur, c'est-à-dire la direction verticale. En effet: L'équation générale de la dynamique s'écrit: Fd+ = c d Yc (1) m m Dans cette formule Fd est la somme des forces dites "à distance" qui s'appliquent à l'aérodyne, c'est-à-dire la pesan- teur et l'inertie, Ec est la somme des forces de contact men- tionnées ci-dessus, yc est l'accélération cinématique de l'aé- rodyne et m la masse de ce dernier. Les forces d'inertie sont les forces de Coriolis qui sont ici négligeables. Dans le cas d'un vol rectiligne à vitesse uniforme, Yc = 0. Les accéléromètres fournissent alors les composantes de l'accélération de la pesanteur. Les composantes de g sont déterminées par rapport aux axes des capteurs que constituent les accéléromètres. Comme ces capteurs sont fixés au bâti de l'avion, il est aisé de déterminer les composantes de g dans le repère de l'avion. Les accéléromètres fournissent donc une indi- cation sur l'orientation de l'avion. Par ailleurs, on peut observer que les accéléromètres ne fournissent une composante du vecteur g que lors de phases de vol rectiligne uniforme de l'avion. Dans les autres phases, on peut connaître la direction verticale, à condition de dispo- ser d'autres données, comme on le verra plus loin. Dans l'exemple, on prévoit trois gyromètres disposés suivant un trièdre trirectangle colinéaire aux axes du boîtier 10 de l'appareil (figure 1). On prévoit donc un gyromètre selon l'axe x, un selon l'axe y et un selon l'axe z. Par contre, on prévoit seulement deux accéléromètres, l'un selon l'axe x et l'autre selon l'axe y. En effet, pour déterminer la direction du vecteur g d'accélération de la pesanteur, deux composantes du vecteur y = g suffisent, car la norme g du vecteur g étant connue, on peut en déduire la com- posante yz selon l'axe Z à partir de cette norme et des com- posantes Yx et yy grâce à la formule suivante: d2 _2 2 Yz = g- Yx - y (2) Dans les appareils connus, pour corriger les données gyrométriques par les données fournies par les accéléromètres, on prévoit une boucle d'asservissement dans laquelle le vecteur ) fourni par les gyromètres est intégré pour fournir les atti- tudes. On intègre les accélérations en utilisant les attitudes gyrométriques pour effectuer des projections de changement de repère et on compare les vitesses ainsi obtenues à une référence pour en déduire les corrections qui sont utilisées pour cor- riger, en boucle fermée, les données gyrométriques. Contrairement aux traitements effectués par les appa- reils classiques dans lesquels on intègre les données fournies par les accéléromètres, l'invention consiste à comparer directe- ment le vecteur de direction verticale fourni par les accéléro- mètres au vecteur de direction verticale fourni par la matrice d'attitudes que délivre le dispositif gyrométrique. Dans le mode de réalisation préféré de l'invention, cette comparaison est obtenue en élaborant le produit vectoriel des deux vecteurs unitaires Ua (direction verticale fournie par les accéléromètres) et Ug (direction verticale fournie par le dispositif gyrométrique) de direction verticale, c'est-à- dire: E: = Ua A Ug (3) L'appareil selon l'invention est donc particulièrement simple car, d'une part, il n'est pas nécessaire d'intégrer les données fournies par les accéléromètres et, d'autre part, les vecteurs Ua et Ug sont dans le même repère qui est celui des capteurs, c'est-à-dire du boîtier 10. En outre, le vecteur E est exploitable pour effectuer la correction puisque sa norme est sine, a étant l'angle entre les vecteurs Ug et Ua. La correction consiste à faire tourner le vecteur Ug d'un angle a de façon à le ramener dans la direc- tion du vecteur Ua. On utilise ce vecteur c dans une boucle d'asservissement qui fait converger le vecteur Ug vers le vec- teur Ua. Cette convergence s'effectue avec une vitesse de rota- tion égale à la norme du vecteur g et donc proportionnelle (si l'angle a est petit) à l'écart angulaire entre les vecteurs Ug et 0a. De plus, la correction est isotrope, c'est-à-dire indépendante du repère dans lequel on se trouve. Toutefois, on effectuera les calculs d'asservissement dans le repère capteur. Le vecteur Ug est fourni par la dernière colonne de la matrice d'attitudes du dispositif gyrométrique. On va maintenant décrire en relation avec la figure 2 un circuit d'asservissement d'un appareil conforme à l'inven- tion. Ce circuit d'asservissement comporte: un moyen de calcul d'intégration, ou intégrateur 12, fournissant un quater- nion d'attitudes, un moyen 14 pour extraire la matrice B d'at- titudes des quaternions, et un bloc 16 qui extrait de la matrice B la dernière colonne qui représente le vecteur Ug, c'est-à- dire le vecteur unitaire de direction verticale dans le repère XYZ du boîtier 10. Le vecteur Ug fourni par le bloc 16 est appliqué sur la première entrée d'un bloc 18 de produit vectoriel dont la seconde entrée reçoit le vecteur vertical unitaire Ua fourni par le dispositif accélérométrique. Ainsi, le bloc 18 délivre sur sa sortie le vecteur e. Pour faire converger le vecteur Ug vers le vecteur Uar on prévoit une première boucle d'asservissement proportion- nel comportant un amplificateur (ou atténuateur) 20 de coeffi- cient k1. Cet asservissement proportionnel est complété par un asservissement à intégrateur 22 permettant d'estimer et de com- penser les dérives du dispositif gyrométrique. On prévoit donc une seconde boucle avec un second amplificateur (ou atténuateur) 24 de coefficient k2. Le signal de la première boucle est fourni sur une première entrée (moins) d'un soustracteur 26 dont la sortie est reliée à l'entrée de l'intégrateur 12 et qui comporte une entrée + (plus) 28 recevant le vecteur Qg fourni par les trois gyromètres et une seconde entrée connectée à la sortie de l'intégrateur 22 en série avec l'amplificateur 24. Dans chacune des boucles d'asservissement, en série avec les amplificateurs 20 et 24, on prévoit des interrupteurs, respectivement 32 et 34, qui sont ouverts pendant les phases dynamiques de vol. En effet, pendant ces phases, les accéléro- mètres ne fournissent plus les composantes de l'accélération de la pesanteur. Cependant, pendant ces ouvertures des boucles d'asservissement, on continue à compenser les dérives des gyro- mètres grâce à la sortie de l'intégrateur 22. Il est à noter que le schéma représenté sur la figure 2 est une représentation vectorielle. Chaque ligne fléchée représente une information vectorielle de dimension 3. Il y a donc, en pratique, trois boucles d'asservissement, une pour chaque composante. Le repère le plus approprié est le repère (X, Y, Z) lié à l'aérodyne. Pour détecter les phases dynamiques pendant lesquelles les boucles d'asservissement doivent être ouvertes, on compare les composantes du vecteur ú à des seuils. Si les seuils sont dépassés, on ouvre la ou les boucles. Comme on le verra plus loin, une boucle d'asservissement dans une direction peut être ouverte, alors que d'autres boucles d'asservissement restent fermées. Quoiqu'il en soit, il est nécessaire de sélectionner soigneusement les instants des ouvertures afin de ne pas fausser l'asservissement. Il faut également éviter que le vecteur Ug dérive de façon excessive pendant les ouvertures de boucles. On décrira plus loin des moyens permettant d'effectuer les ouvertures des boucles d'asservissement à des instants correctement choisis, ainsi que des sécurités permettant d'éviter les divergences du vecteur Ug pendant les ouvertures des boucles. On va maintenant décrire, en relation avec la figure 3, un circuit d'asservissement qui utilise des boucles d'asser- vissement analogues à celles représentées sur la figure 2 et, en plus, une autre boucle d'asservissement permettant de compenser la composante verticale de dérive dQ. * En effet, le produit vectoriel E étant, par défini- tion, perpendiculaire au vecteur Ug, sa composante suivant la direction verticale est nulle. Ce produit vectoriel ne permet donc pas d'estimer la composante verticale de la dérive dO g si l'avion conserve une attitude constante. Cette compensation est nécessaire pour pouvoir obtenir correctement les trois composantes de la rotation instantanée Q. En outre, selon un aspect de l'invention, la composante ver- ticale de la rotation instantanée constituant le taux de virage, cette composante sert à distinguer les phases rectilignes des phases de virages pour ouvrir, à des moments correctement choi- sis, certaines des boucles d'asservissement. Dans le mode de réalisation représenté sur la figure 3, on effectue une projection orthogonale du vecteur Qc, c'est- à-dire le vecteur de rotation instantanée compensé des dérives estimées, sur la direction Ug. A cet effet, on effectue le pro- duit scalaire Qc.qug, et on utilise le vecteur (Qc.Ug).Ug pour effectuer l'asservissement de lacet. Pour éviter toute confusion, on a noté Og le vecteur d'accélération angulaire sans compensation et Qc le même vec- teur avec compensation. Sur la figure 3, les éléments correspondant à ceux de la figure 2 portent les mêmes chiffres de références. Dans cette réalisation, le vecteur Qg est appliqué à l'entrée + d'un soustracteur 40 dont la sortie est reliée à l'entrée + d'un autre soustracteur 42 dont l'entrée - reçoit la sortie de l'amplificateur 20 (kl). Par ailleurs, l'entrée - du soustracteur 40 est reliée à la sortie d'un additionneur 44 par l'intermédiaire de l'intégrateur 22. Une première entrée de l'additionneur 44 est reliée à la sortie de l'amplificateur 24 (k2), tandis que la seconde en- trée de cet additionneur est reliée à la sortie d'un amplifi- cateur 46 de gain k3 se trouvant dans la boucle d'asservissement de lacet. Cette boucle d'asservissement de lacet comporte un bloc 48 qui calcule les composantes du vecteur (QC.Ug).Ug. A cet effet, ce bloc 48 présente une entrée 50 recevant le vecteur Qc (Qg compensé des dérives, du fait de la boucle d'asservis- sement à amplificateur 24 et intégrateur 22), et une seconde entrée 5recevant le vecteur!g fourni par la sortie du bloc 16. La sortie du bloc 48 est reliée à l'entrée de l'amplifi- cateur 46 par l'intermédiaire d'un interrupteur 54. Étant donné que la boucle d'asservissement à amplifi- cateur 46 doit compenser une dérive, elle comporte un intégra- teur, et ce dernier est l'intégrateur 22, qui fait également partie de la boucle à amplificateur 24. Le produit scalaire Qc-.Ug représente la composante verticale de la dérive du dispositif gyrométrique à la condition que l'avion ne tourne pas autour de la verticale. C'est pourquoi l'interrupteur 54 est ouvert pendant les virages. On va maintenant décrire, notamment en relation avec les figures 4 à 6, un certain nombre de dispositions permettant de commander correctement l'ouverture des interrupteurs 32, 34 et 54 des circuits d'asservissement pendant certaines phases transitoires pour lesquelles ces circuits (ou une partie d'entre eux) ne peuvent pas fonctionner correctement. On considérera tout d'abord l'interrupteur 32 (figures 2 et 3) du circuit d'asservissement de type proportionnel. Cet interrupteur 32 doit être ouvert pendant les accé- lérations ou les virages, car, dans ce cas, le vecteur Ua fourni par les accéléromètres ne représente pas la direction verticale réelle. Pour détecter ces phases, on détecte les écarts angulaires entre les vecteurs Ug et Ua qui est repré- senté par les composantes du vecteur = Ua ^ Ug. Pour limiter les ouvertures de l'interrupteur 32, on confère une caractéristique de saturation à l'amplificateur 20. Cette caractéristique est représentée sur la figure 4, qui est un diagramme montrant, pour une composante (x par exemple), la variation de la correction c portée en ordonnées, de l'ampli- ficateur 20 en fonction de la composante Ex du vecteur E, cette composante cx étant portée en abscisses. Cette variation est représentée par une ligne 60 présentant: une zone linéaire de pente K1 représentée par un segment de droite 62 passant par l'origine O, des zones de saturation représentées par les seg- ments de droites 64 et 66 parallèles à l'axe des abscisses, et des zones de coupure au-delà des abscisses exc et -xc. Lorsque l'écart sx dépasse, en valeur absolue, la valeur Exc, l'interrup- teur 32 est ouvert. La saturation limite le signal appliqué sur l'entrée - du soustracteur 42. Cette limite du signal correspond au maximum admissible pour la compensation des dérives. Une telle caracté- ristique est un compromis entre une correction proportionnelle et la boucle ouverte. Cette caractéristique de saturation atté- nue l'effet des erreurs de E sur la précision de Ug, sans avoir à couper l'asservissement. Pour commander l'ouverture de l'interrupteur 34, on compare les normes de g et de Qc à des seuils pour détecter les évolutions de l'aérodyne, lorsqu'on passe en phase de virage ou d'accélération. L'interrupteur 54 est ouvert lorsque le taux de virage dépasse un seuil. Ce taux de virage est représenté par le pro- duit scalaire Qc.Ug. On va maintenant décrire, en relation avec les dia- grammes des figures 5 et 6, une disposition qui permet de limi- ter le risque de divergence pendant les ouvertures de l'inter- rupteur 32X lors des phases de virage. Cette disposition consiste à maintenir fermé l'inter- rupteur 32y à partir du moment o l'interrupteur 32X est ouvert, quel que soit sy, en remplaçant, dans ce cas, le profil de l'am- plificateur 20y tel que décrit avec la figure 4 par celui décrit avec la figure 5, c'est-à-dire avec une pente plus faible et continue, sans seuil. Autrement dit, le gain kly de l'amplificateur 20y pré- sente une valeur plus faible pendant les phases de virages. Dans le diagramme de la figure 5, on a porté en abscisses la compo- sante gy du vecteur -, et, en ordonnées, le signal ry en sortie de l'amplificateur 20y. La ligne 60y représentée en traits interrompus correspond à la ligne 60 de la figure 4. Cette variation de ry en fonction de sy intervient pendant les phases de vol rectiligne sans accélération, l'interrupteur 32X de la boucle d'asservissement roulis étant fermée dans ce cas. Par contre, pendant les phases de virage, l'interrupteur 32X est ouvert et l'interrupteur 32y reste fermé. Mais, comme représenté par la ligne 72y, la caractéristique de l'amplificateur 20y est modifiée dans cette situation: la pente de la droite 72y est inférieure à la pente de la droite 62y. Avec cette disposition, lorsque l'avion tourne autour de la verticale, les erreurs d'attitude en tangage et en roulis convergent vers zéro suivant une spirale comme représentée par la figure 6, si la vitesse de l'aérodyne reste constante. Ainsi, en régime permanent, on obtient des erreurs nulles en tangage et roulis. Sur le diagramme de la figure 6, on a porté en ordon- nées l'erreur ô4 (roulis) et, en abscisses, l'erreur 80 (tan- gage). La courbe en spirale 74 converge vers des erreurs nulles en (roulis) et en 0 (tangage). Dans le cas d'un virage, s'il se produit des accélé- rations ou décélérations longitudinales (axe x), la référence de verticale Ua fournie par les accéléromètres devient fausse aussi en tangage. Cette erreur AO en tangage sur le vecteur de direction verticale Ua se traduit par une erreur sur le roulis Q. En régime permanent, cette erreur est proportionnelle à AO, K1 et 1/Q, Q étant le taux de virage et par une erreur nulle pour 0. La réduction du gain de l'amplificateur rend ainsi le dispositif moins sensible aux erreurs de correction dues aux accélérations longitudinales. On va maintenant décrire une autre disposition qui permet de fiabiliser la distinction entre les phases rectilignes et les phases de virage de façon à maintenir fermé l'interrup- teur 32X lorsqu'on détecte une phase de vol rectiligne. On a vu que le maintien à l'état ouvert ou fermé des interrupteurs 32 dépend de la comparaison de chaque composante de È à un seuil. On a vu aussi qu'on pouvait, pour distinguer entre les phases de vol rectiligne et les phases de virages, utiliser aussi le produit scalaire Qc. Ug qui représente le taux de virage. Selon l'invention, on interdit l'ouverture de 32x lorsque úx dépasse son seuil de coupure si le taux de virage calculé indique que l'on n'est pas en virage. Dans ce cas, on remplace le profil de l'amplificateur 20y tel que représenté sur la figure 4 par un profil linéaire comme la droite 72 sur la figure 5. Toutefois, la pente de cette droite est supérieure à la pente de la droite 62. De cette manière, on ramène la direc- tion verticale du dispositif gyrométrique dans la direction ver- ticale vraie 0a fournie par les accéléromètres, dans le cas o, en sortie de virage, úx serait sorti de sa plage de capture (zone -Ec +Ec sur le diagramme de la figure 4) à la suite d'une divergence pendant la phase en boucle ouverte. On va maintenant décrire, en relation avec les figures 7 et 8, un mode de réalisation de l'intégrateur 22 utilisé pour estimer les dérives. On rappelle ici que l'intégrateur 22 doit pouvoir com- penser les dérives en boucle ouverte, c'est-à-dire quand les asservissements sont inhibés. L'intégrateur 22 tient compte des caractéristiques des dérives du dispositif gyrométrique. En particulier, l'évolution des dérives en fonction du temps étant limitée par une pente maximale caractéristique des gyromètres, on limite également la pente de sortie de l'intégrateur 22. La figure 7 représente un système échantillonné, de -1 période dt. Chaque bloc Z est un retard, ce qui signifie que la sortie est en retard d'une période par rapport à l'entrée. Dans l'exemple représenté sur la figure 7, l'intégra- teur 22 comporte deux étages, respectivement 90 et 92. Le pre- mier étage comporte un retard 94 dont l'entrée est reliée à la sortie d'un additionneur 96 dont une entrée 961 constitue l'en- trée de l'intégrateur 22 et dont l'autre entrée 962 est reliée à la sortie 222 de l'intégrateur 22. La sortie du retard 94 est reliée à l'entrée + d'un soustracteur 98 dont la sortie est reliée à l'entrée d'un satu- rateur d'écart 100 limitant les signaux entre deux valeurs +Emax et - Emax de signes opposés. La sortie du saturateur 100 est reliée à une entrée d'un additionneur 102 dont la sortie est reliée à la sortie 222 de l'intégrateur 22. Le second étage 92 comporte aussi un retard 104 dont la sortie est reliée à la seconde entrée de l'additionneur 102 et dont l'entrée est reliée à la sortie d'un additionneur 106 dont une entrée 106l1 est reliée à la sortie d'un saturateur de pente 108 dont l'entrée est reliée à la sortie du soustracteur 98 en amont du saturateur d'écart 100. La sortie du retard 104 est également connectée à l'entrée - du soustracteur 98, et à la deuxième entrée de l'ad- ditionneur 106. Le saturateur de pente 108 permet de limiter la pente du second étage 92 à deux pentes maximales de signes opposés Pmax, -Pmax- Le premier étage 90 intègre directement le signal d'entrée. Sa sortie est limitée par le saturateur 100 dans une zone de largeur fixée dont le centre est le signal de sortie du second étage 92. Le signal de sortie du second étage évolue selon deux pentes maximales opposées Pmax et -Pmax déterminées par le satu- rateur 108. Le second étage 92 sert à limiter la pente moyenne du signal intégré, tandis que la sortie de l'intégrateur 22 est fournie par le premier étage. L'évolution du signal de sortie est représentée sur la figure 8. On voit que le signal de sortie 105 est confiné entre deux lignes brisées 107 et 109 dont tous les segments sont parallèles entre eux. Chacun des segments des lignes 107 et 109 ayant alternativement une pente Pmax et une pente -Pmax. Tant que l'écart entre les signaux des deux étages 90 et 92 reste inférieur à l'écart maximum 2Emax (parties 110, 112, etc. de la courbe 105), l'estimateur a un fonctionnement linéaire, c'est-à- dire qu'il fonctionne comme un simple intégrateur. On conserve ainsi les propriétés dynamiques de l'asservissement autour de la valeur zéro. Par contre, lorsque l'écart dépasse l'écart maximum déterminé par le saturateur 100, la sortie de l'estimateur est "bridée" afin que son signal n'évolue pas plus que la pente maximale autorisée Pmax. Ce bridage est représenté par les seg- ments 111, 113, etc. sur la figure 8. En résumé, l'intégrateur 22 représenté sur la figure 7 permet de limiter la vitesse moyenne d'évolution du signal d'es- timation et non de limiter la vitesse instantanée de ce signal d'estimation. Ces caractéristiques permettent de s'affranchir des difficultés qui pourraient résulter de vols avec de nombreux virages et donc avec des ouvertures fréquentes de boucles d'as- servissement. Dans tous les exemples décrits ci-dessus, on a indiqué qu'il était nécessaire de prendre des précautions lors des phases de vol non rectilignes et/ou accélérées, car la correc- tion fournie par les accéléromètres ne peut être utilisée pen- dant ces phases. En effet, dans le cas de vol en virage ou accé- léré, les mesures accélérométriques ne fournissent plus la ver- ticale vraie liée à la pesanteur mais la verticale apparente liée à la portance de l'aérodyne. Cependant, dans l'équation (1) mentionnée cidessus, le vecteur Yc d'accélération cinématique peut être déterminé à partir de la vitesse V et de la rotation instantanée Q selon la formule classique ci-dessous: dv Yc = + Q A V (4) dt Dans cette formule, les vecteurs V et Q ainsi que d9 l'accélération - sont considérés en repère aérodyne. dt On voit ainsi que si l'on connaît le vecteur V en repère aérodyne, on sait déterminer Yc, car Q est fourni par les gyromètres. Autrement dit, si le vecteur V peut être déter- miné, les données gyrométriques peuvent être corrigées en boucle fermée, même lors de virages ou d'accélérations. La vitesse de l'aérodyne, V, peut être déterminée par un anémomètre. Ainsi, quand l'aérodyne vole sans dérapage à incidence nulle ou faible, la vitesse en repère aérodyne est un vecteur V dont la composante Vx est la vitesse fournie par l'anémomètre et les composantes Vy et Vz sont nulles (étant donné qu'il n'y a pas de dérapage et que l'incidence est nulle). Un anémomètre peut ainsi être utilisé soit pour main- tenir constamment fermées les boucles d'asservissement, soit pour limiter les périodes d'ouverture de ces boucles. On va maintenant décrire, notamment en relation avec la figure 9, une disposition qui peut être utilisée en combinai- son avec l'utilisation de la vitesse fournie par un anémomètre pour corriger les erreurs et dérives des gyromètres. Cette disposition part de la constatation qu'en phase de vol rectiligne, la référence de verticale est très sensible aux erreurs sur le vecteur Q, c'est-à-dire aux erreurs d'esti- mation de dérives. En d'autres termes, la précision sur le vec- teur de direction verticale dépend fortement des erreurs sur le vecteur Q; en outre, cette précision se dégrade quand la vitesse augmente. L'estimation de la dérive de lacet, qui a été décrite avec la figure 3, n'est pas suffisamment précise pour permettre de compenser cette altération de précision. Ainsi, selon une autre disposition de l'invention, on remplace l'asservissement de lacet par un asservissement du cap gyrométrique sur le cap magnétique fourni par un magnétomètre. Cet asservissement est analogue à celui des attitudes, mais ne nécessite pas d'ouverture de boucle d'asservissement pendant les virages. Cet asservissement en cap permet, bien entendu, de limiter les erreurs de cap et améliore l'estimation de la dérive de lacet. On améliore la précision du calcul de la verticale vraie dans l'équation (4) lorsqu'on utilise la vitesse anémo- métrique. On peut aussi améliorer la détection des phases recti- lignes de vol. Enfin, les contraintes de qualité sur le gyro- mètre de lacet sont moins sévères, ce qui permet l'utilisation d'un gyromètre de lacet plus économique. Un mode de réalisation de cette disposition de l'in- vention consiste à fournir un vecteur de correction de cap sc de direction verticale afin que cette correction n'entre pas en conflit avec celle fournie par le vecteur E mentionné ci-dessus qui est par construction, dans le plan horizontal (produit vec- toriel des directions verticales fournies par le dispositif gyrométrique et par les accéléromètres). L'écart de cap peut être obtenu grâce à un magnéto- mètre ou appareil analogue qui fournit les composantes du champ magnétique terrestre B. Dans ces conditions, l'écart de cap est calculé selon la formule ci-dessous: EC = [(Um A b1îb31b3 = ( m.b2})3 (5) Dans cette formule, le vecteur Um est le vecteur uni- taire dans la direction du nord magnétique et b1, b2 et b3 sont les trois vecteurs fournis par chacune des trois colonnes de la matrice d'attitudes B avec: bll b12 b13 B = b21 b22 b23 = [1b23] (6) b31 b32 b33 Dans ce cas, le vecteur b1 est le vecteur d'attitude dans la direction nord, le vecteur b2 est le vecteur d'attitude dans la direction est et le vecteur b3 est le vecteur d'at- titude en direction verticale. L'équation (5) signifie que l'on effectue le produit vectoriel de Um par b1 et que l'on projette le vecteur résul- tat sur la verticale donnée par b3. La figure 9 représente, de façon simplifiée, les bou- cles d'asservissement utilisées dans ce mode de réalisation o on effectue une correction de cap. Cette figure 9 est analogue à celle de la figure 3. Toutefois, comme déjà indiqué ci-dessus, l'asservissement de lacet est remplacé par l'asservissement en cap. On voit que l'intégrateur 22 est commun aux boucles d'as- servissement d'attitude et de cap. On a représenté par un bloc l'organe qui détermine les vecteurs gc d'écart d'attitudes (obtenu à partir du produit vectoriel des vecteurs de direction verticale fournis, d'une part, par le dispositif gyrométrique et, d'autre part, par les accéléromètres) et cc d'écart de cap. Ce vecteur Ec est obtenu à partir de la matrice d'attitudes et du vecteur B de champ magnétique terrestre. On prévoit, dans l'asservissement en cap, une correc- tion de type proportionnel avec un gain a (amplificateur 122) et une correction de type intégral de gain P (amplificateur 124), cette dernière correction faisant appel au même intégrateur 22 que l'asservissement en attitude de gain k2. Quel que soit le mode de réalisation, l'appareil com- porte, de façon classique, des moyens de réinitialisation habi- tuellement dénommés "caging". Cette réinitialisation est déclen- chée de façon manuelle par le pilote pour faire converger rapi- dement les attitudes afin que la direction verticale déterminée par le dispositif gyrométrique corresponde à la verticale réelle. Cette fonction peut être utilisée en vol à vitesse constante rectiligne, ou dans n'importe quelle phase de vol, si on combine les données fournies par un anémomètre, ou analogue, avec les données fournies par les accéléromètres. Pour assurer cette fonction de réinitialisation, dans un mode de réalisation de l'invention, les gains des amplifica- teurs des boucles d'asservissement proportionnel sont majorés pendant un temps déterminé, par exemple trois minutes, sur les trois axes et, pendant cette durée, les boucles d'asservissement sont maintenues à l'état fermé et les saturations sont inhibées. Les intégrateurs étant bridés, les boucles d'asservissement per- mettent de faire converger rapidement les erreurs d'attitudes afin qu'elles puissent rentrer dans la plage de capture de ces asservissements. REVENDICATIONS 1. Appareil de détermination d'attitudes, notamment l'assiette (0) et l'inclinaison (W) d'un aérodyne, comportant: un dispositif gyrométrique à composants liés fournis- sant les composantes du vecteur (Q) de rotation instantanée de l'aérodyne dans un repère lié à l'aérodyne, des moyens pour calculer, à partir des mesures gyromé- triques, une matrice d'attitudes (b1, b2, b3) définissant l'orientation de l'aérodyne par rapport au repère géographique local, des accéléromètres solidaires du dispositif gyromé- trique, et des moyens pour comparer des données calculées à par- tir des accéléromètres à des données calculées à partir des gyromètres afin de fournir des corrections pour compenser les erreurs ou dérives du dispositif gyrométrique, caractérisé en ce que la comparaison est effectuée entre un vecteur Ug de direction verticale fourni par le dispo- sitif gyrométrique et un vecteur Ua de direction verticale fourni par les accéléromètres. 2. Appareil selon la revendication 1, caractérisé en ce qu'il comporte des moyens pour calculer, à partir des mesures accélérométriques, la direction de la verticale dans un repère lié à l'aérodyne. 3. Appareil selon la revendication 2, caractérisé en ce que les vecteurs de direction verticale sont définis par leurs coordonnées dans un repère lié à l'aérodyne. 4. Appareil selon la revendication 1, 2 ou 3, carac- térisé en ce que, le vecteur de direction verticale fourni par le dispositif gyrométrique est obtenu à partir de la dernière colonne (b3) de la matrice d'attitudes. 5. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce que les vecteurs de direction verticale ont une longueur unitaire. 6. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce que les moyens de comparaison effectuent le produit vectoriel (g) des deux vecteurs de direc- tion verticale. 7. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une boucle d'as- servissement de type proportionnel. 8. Appareil selon la revendication 7, caractérisé en ce que la boucle d'asservissement de type proportionnel présente un gain limité à une valeur maximale prédéterminée (64, 66). 9. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend une boucle d'as- servissement à intégrateur (22) pour compenser les dérives. 10. Appareil selon la revendication 6, caractérisé en ce qu'il comporte un moyen pour effectuer une correction suivant la direction verticale. 11. Appareil selon la revendication 10, caractérisé en ce que le moyen pour effectuer une correction suivant la direc- tion verticale comprend un moyen de calcul de la projection, sur la direction verticale, du vecteur (Q) de rotation instantanée fourni par le dispositif gyrométrique, la direction verticale sur laquelle est effectuée la projection étant fournie par le dispositif gyrométrique, et le vecteur de rotation instantanée qui est projeté étant corrigé par la boucle d'asservissement. 12. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce que le vecteur de direction ver- ticale est fourni, d'une part, par les accéléromètres et, d'autre part, à partir de la vitesse de l'aérodyne délivrée, par exemple, par des moyens anémométriques. 13. Appareil selon la revendication 12, caractérisé en ce que la direction verticale est déterminée à partir de l'équa- tion suivante: Fd dV FC m d + A Vm m dt m Fd étant la force de la pesanteur, de direction verticale, V le vecteur de vitesse de l'aérodyne déterminé par les moyens anémo- métriques, Q le vecteur de rotation instantanée de l'aérodyne déterminé par le dispositif gyrométrique, et Fc la résultante des forces de contact déterminée par les accéléromètres, et m la masse de l'aérodyne. 14. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte un moyen d'asser- vissement du cap gyrométrique sur le cap magnétique. 15 Appareil selon les revendications 6 et 14, carac- térisé en ce que l'asservissement du cap est déterminé à partir du vecteur de correction (Ec) satisfaisant à l'équation suivante: C = [(Um ^A bîlb31b3 = (jm'b2b3 (5) équation dans laquelle Um est un vecteur représentant la direc- tion du nord magnétique, de préférence d'amplitude unité, et les vecteurs b1, b2 et b3 sont trois vecteurs dont les coordonnées sont fournies par les colonnes de la matrice d'attitudes déli- vrée par le dispositif gyrométrique. 16. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend des boucles fer- mées d'asservissement des attitudes utilisant les signaux de correction fournis par les moyens de comparaison, et des moyens pour ouvrir au moins une boucle en dehors des phases de vol rec- tiligne et à vitesse constante de l'aérodyne. 17. Appareil selon la revendication 16, caractérisé en ce que l'ouverture des boucles d'asservissement intervient quand l'écart entre le vecteur de direction verticale fourni par le dispositif gyrométrique et le vecteur de direction verticale fourni par les accéléromètres dépasse une valeur prédéterminée (Exc) - 18. Appareil selon la revendication 16 ou 17, carac- térisé en ce qu'une boucle d'asservissement en roulis étant ouverte pendant un virage, une boucle d'asservissement corres- pondante en tangage reste fermée pendant ce virage. 19. Appareil selon l'une quelconque des revendications 16 à 18, caractérisé en ce que, le taux de virage étant utilisé pour distinguer entre des phases de vol rectiligne et des phases de vol différentes, une boucle d'asservissement en roulis est fermée quand le taux de virage est inférieur en valeur absolue à un seuil déterminé. 20. Appareil selon les revendications 9 et 10, carac- térisé en ce qu'il comporte une boucle d'asservissement pour assurer une correction des mesures gyrométriques suivant la direction verticale et en ce que cette boucle d'asservissement comprend un intégrateur (22) commun à la boucle compensant les dérives. 21. Appareil selon la revendication 20, caractérisé en ce que l'intégrateur (22) comporte des moyens (108) permettant de limiter la pente moyenne de son signal de sortie et des moyens (100) pour que, à chaque instant, le signal de sortie présente un écart inférieur à un seuil par rapport à un signal à variation linéaire ayant une pente déterminée.