FR 2728339 A1 19960621 FR 9415050 A 19941214 La présente invention concerne un procédé pour estimer les biais gyrométriques de mesures effectuées par au moins un gyromètre monté sur un aéronef, ainsi qu'un dispositif pour la mise en oeuvre de ce procédé. Dans le cadre de la présente invention, on entend par biais des erreurs de mesure d'un capteur apparaissant lors de chaque mesure et dues par exemple à un mauvais étalonnage ou à un défaut de précision dudit capteur, et donc par biais gyrométriques on entend des erreurs de mesure ainsi définies d'un gyromètre. Bien que non exclusivement, la présente invention sera plus particulièrement appropriée à être mise en oeuvre sur des missiles lancés à partir du sol et munis d'une centrale inertielle. De façon connue, une telle centrale inertielle est pourvue de capteurs, à savoir essentiellement des accéléromètres et des gyromètres, destinés à effectuer des mesures utilisées par ladite centrale inertielle, en parti- culier pour déterminer la position du missile en vol. En raison de l'imperfection de ces capteurs entraînant l'appa- rition de biais dans les mesures effectuées, la position du missile estimée par la centrale inertielle diffère, rapide- ment après le lancement du missile, de la position effec- tive. Il est alors nécessaire, si on veut éviter un échec de la mission du missile, d'effectuer un recalage de la posi-25 tion dudit missile. Or, un tel recalage est généralement difficile à mettre en oeuvre et de plus très coûteux. Une solution pour éviter d'avoir à effectuer un tel recalage est d'utiliser des capteurs à très haute performance, dont les éventuels biais sont négligeables. Ceci concerne surtout30 les gyromètres, puisque les biais gyrométriques entraînent les erreurs les plus importantes dans l'estimation de la position du missile. Toutefois, de tels gyromètres à très haute performance sont extrêmement coûteux, de sorte que cette solution n'est guère envisageable, surtout pour des aéronefs ou des missiles munis d'une centrale inertielle de faible coût. La présente invention a pour objet de remédier à ces in- convénients. Elle concerne un procédé pour estimer les biais gyrométriques de mesures effectuées par un ensemble de gyromètres comportant au moins un gyromètre et monté sur un aéronef. A cette fin, selon l'invention, ledit procédé est remarqua- ble en ce que: - on immobilise ledit aéronef en un endroit déterminé sur la terre; - on détermine les angles d'assiette, de roulis et de cap dudit aéronef ainsi immobilisé; -on détermine la latitude dudit endroit déterminé; -on calcule, à partir de ladite latitude et desdits angles d'assiette, de roulis et de cap, la vitesse de rotation de la terre dans un référentiel de mesure lié audit aéronef; - on mesure la vitesse de rotation de la terre dans ledit référentiel de mesure au moyen dudit ensemble de gyromè- tres; et - on effectue la soustraction entre la vitesse de rotation calculée et la vitesse de rotation mesurée, de manière à obtenir les valeurs estimées desdits biais gyrométriques. Ainsi, grâce à l'invention, on est en mesure de déterminer les biais gyrométriques de mesures effectuées par des gyromètres, ce qui permet de corriger lesdites mesures de30 manière à obtenir des valeurs gyrométriques extrêmement précises. Par conséquent, on peut utiliser des gyromètres de performance moyenne, et donc de faible coût, tout en obte- nant des résultats très précis (après la correction des biais gyrométriques déterminés par le procédé conforme à l'invention). La présente invention est donc particulièrement appropriée à être mise en oeuvre sur des aéronefs comportant des centra- les inertielles de faible coût ou sur des aéronefs, en particulier des missiles, à longue portée présentant un temps de vol et une distance de vol élevés puisque, dans ce dernier cas, même si l'aéronef est muni d'une centrale inertielle très performante, une légère erreur dans les mesures gyrométriques qui doivent être répétées sur une longue durée peut provoquer une erreur importante, par exemple dans l'estimation de la position de l'aéronef, ce qui peut entraîner des conséquences très dommageables. Or, grâce à l'invention, de telles erreurs préjudiciables, même si elles ne sont que minimes, peuvent être facilement corrigées. Avantageusement, les angles d'assiette et de roulis sont déterminés à partir de mesures accélérométriques effectuées sur l'aéronef, de préférence au moyen d'accéléromètres de la centrale inertielle de l'aéronef. Dans le cas o le référentiel de mesure correspond à un référentiel OXYZ dans lequel: - l'axe OX correspond à l'axe de l'aéronef, - l'axe OY est perpendiculaire audit axe OX suivant les ailes dudit aéronef, et - l'axe OZ est perpendiculaire auxdits axes OX et OY, on détermine de façon avantageuse les angles d'assiette e et de roulis * à partir d'au moins deux des trois relations suivantes: rx = -g.sine.(l+fx) ry = g. cose.sino.(l+fy) rz = g.cose.coso.(l+fz) dans lesquelles: - rx, ry et rz représentent des mesures accélérométriques effectuées par des accéléromètres respectivement selon les axes OX, OY et OZ, - _ représente l'accélération de la pesanteur, et - fx, fy et fz représentent des facteurs d'échelle desdits accéléromètres, correspondant à une erreur multiplicative à la différence des biais qui correspondent à une erreur additionnelle. On notera que l'angle de cap ne peut pas être déterminé à partir des mesures accélérométriques précisées précédemment puisqu'il est défini dans un plan qui est perpendiculaire à l'accélération de la pesanteur. Aussi, de préférence, on détermine ledit angle de cap à partir de mesures effectuées par un moyen de mesure appro- prié externe à la centrale inertielle de l'aéronef, ce qui permet de plus d'utiliser un moyen de mesure non embarqué à bord de l'aéronef et pouvant être réutilisé ultérieurement. En outre, de façon avantageuse, on calcule la vitesse de rotation wcal de la terre dans le référentiel de mesure à partir de la relation: / cose.cosp.w.cosÀ + sine.w.sinA wcal = ( cosî.sine.sin0 - sinp. coso).w.cosÀ - cose.sins.w.sinX (coso.sine.cos + sin0.sin>).w.cosÀ cosb.cose4.sinX dans laquelle: - À correspond à la latitude dudit point déterminé, - e, * et T représentent respectivement les angles d'as- siette, de roulis et de cap, et - X représente la vitesse de rotation de la terre autour de son axe de rotation. Par ailleurs, avantageusement, ledit aéronef, par exemple un missile, est immobilisé pour la mise en oeuvre de l'inven- tion à son poste de lancement, ce qui permet de déterminer les biais gyrométriques juste avant le lancement après la mise en service des capteurs nécessaires à cette mise en oeuvre. La présente invention concerne également un dispositif associé audit ensemble de gyromètres comportant au moins un gyromètre et monté sur l'aéronef, pour la mise en oeuvre du procédé précité. Avantageusement, le dispositif conforme à l'invention est associé à un système auxiliaire susceptible de déterminer les angles d'assiette et de roulis de l'aéronef, et compor- te: - un premier moyen de mesure déterminant la latitude dudit endroit déterminé de la terre; - un second moyen de mesure déterminant l'angle de cap dudit aéronef; - un premier calculateur calculant la vitesse de rotation de la terre à partir d'informations reçues dudit système auxiliaire et desdits premier et second moyens de mesure; et - un second calculateur déterminant les biais gyrométriques à partir d'informations reçues dudit premier calculateur et dudit ensemble de gyromètres. De façon avantageuse, ledit système auxiliaire et ledit ensemble de gyromètres font partie de la centrale inertielle30 de l'aéronef. De plus, lesdits premier et second calculateurs sont regrou- pés de préférence dans un calculateur unique. Par ailleurs, les valeurs des biais gyrométriques estimées par ledit dispositif peuvent être stockées avantageusement dans une mémoire intégrée dans ladite centrale inertielle, ce qui permet à ladite centrale inertielle de corriger directement les mesures gyrométriques ultérieures effectuées par ledit ensemble de gyromètres. Les figures du dessin annexé feront bien comprendre comment l'invention peut être réalisée. Sur ces figures, des réfé- rences identiques désignent des éléments semblables. La figure 1 est le schéma synoptique d'un dispositif confor- me à l'invention. La figure 2 illustre de façon schématique dans un plan vertical la position initiale immobile de l'aéronef. La figure 3 illustre schématiquement dans un plan horizontal la position de la figure 2. La figure 4 illustre un référentiel utilisé pour la mise en oeuvre de l'invention et défini par rapport à la terre. Le dispositif 1 conforme à l'invention et représenté schéma- tiquement sur la figure 1 est destiné à estimer les biais gyrométriques de mesures effectuées par des gyromètres d'un ensemble de gyromètres non représenté et faisant partie de la centrale inertielle 2 d'un aéronef A, en particulier un25 missile. On sait que les biais gyrométriques correspondent & des erreurs concernant les mesures gyrométriques, dues par exemple à un défaut d'étalonnage des gyromètres et apparaissant lors de chaque mesure. De tels biais gyrométri- ques entraînent des erreurs dans les calculs effectués à partir de ces mesures par la centrale inertielle 2, ce qui peut évidemment entraîner des conséquences très dommagea- bles, en particulier un échec de la mission à effectuer. La présente invention a pour objet d'estimer ces biais gyrométriques afin de pouvoir les corriger dans les mesures effectuées par les gyromètres, lesdites mesures gyrométri- ques correspondant de façon connue à des mesures de vitesses angulaires. A cet effet, selon l'invention, on immobilise ledit aéronef A par rapport à la terre T, par exemple sur une rampe ou un tube de lancement non représentée, de manière que ledit aéronef A ne soit soumis qu'à la rotation de la terre T. A partir d'une telle position immobile, on détermine les biais gyrométriques en effectuant, conformément à l'inven- tion, et de préférence juste avant le lancement de l'aéronef A, les opérations suivantes: - on calcule la vitesse de rotation w cal de la terre T dans un référentiel de mesure R1 lié à l'aéronef A et décrit ci-après; - on mesure, à l'aide des gyromètres de ladite centrale inertielle 2, la vitesse de rotation mes de la terre T mes dans ledit référentiel de mesure R1; et - on effectue la soustraction entre la vitesse de rotation calculée w cal et la vitesse de rotation mesurée w mes' de manière à obtenir les valeurs estimées desdits biais gyrométriques. Ledit référentiel de mesure R1 est formé d'axes OX, OY et OZ tels que: - l'axe OX correspond à l'axe de l'aéronef A; - l'axe OY est perpendiculaire audit axe OX suivant les ailes A1 et A2 dudit aéronef A; et - l'axe OZ est perpendiculaire à OX et à OY, tel que repré- senté sur les figures 2 et 3. Selon l'invention, pour calculer la vitesse de rotation de la terre T, on définit un référentiel de base R2 représenté sur la figure 4 et défini par trois axes V, Nm et Em tels que: - V est un axe parallèle au gradient du champ de gravité terrestre non représenté, orienté vers le centre C de la terre T, et donnant la verticale locale; -Nm est un axe horizontal donnant le nord magnétique; et - Em est un axe horizontal donnant l'est magnétique. Sur la figure 4, on a représenté le référentiel de base R2 en un point P correspondant à l'endroit du globe terrestre T o se trouve l'aéronef A au moment des mesures. La latitude de ce point P est désignée par À. Par ailleurs, on sait que la terre T tourne autour de son axe de rotation CC' à une vitesse w, avecw = 150 par heure. Aussi, dans ledit référentiel de base R2, la vitesse de rotation wl de la terre T s'écrit: w.CosX Wl = 0 -w.sink On peut déduire de cette expression la vitesse de rotation mcal de la terre T dans le référentiel de mesure R1, à partir de la relation: cal= M.wl mcal M étant la matrice de passage dudit référentiel de base R2 audit référentiel de mesure R1. De façon connue, la matrice M s'écrit: (cose.cosi I cose.sin I - sine M = cosP.aine.sin0-sin4.cos> sin>.sine.sin+cos. cos I cose.sino cos0.sine.cosP+sin0.sin I cos0.sine.sinp-sin0.cos coso. coue avec: - e l'angle d'assiette de l'aéronef A, correspondant à l'angle entre l'axe OX dudit aéronef A et le plan horizon- tal défini par les axes Em et Nm, tel que représenté sur la figure 2; - * l'angle de roulis de l'aéronef A; et - l'angle de cap, c'est-à-dire l'angle par rapport à la direction Nm indiquant le Nord de la projection de l'axe OX sur le plan horizontal formé des axes Nm et Em, comme représenté sur la figure 3. Dans le cas d'un missile lancé à partir du sol au moyen d'une rampe de lancement non représentée, l'angle d'assiette e correspond à l'angle entre cette rampe de lancement et le20 plan horizontal (Nm, Em) lorsque l'axe OX de l'aéronef est parallèle au plan de lancement de ladite rampe. Selon l'invention, les angles e et 0 sont déterminés à partir de mesures accélérométriques rx, ry et rz effectuées par des accéléromètres non représentés de la centrale25 inertielle 2, respectivement suivant les axes OX, OY et OZ. De façon connue, ces mesures accélérométriques vérifient les relations suivantes: rx = -g.sine.(1+fx) ry = g.cose.sino.(l+fy) rz = g. cose.coso.(l+fz) dans lesquelles: - fx, fy et fz représentent des facteurs d'échelle desdits accéléromètres, et - q correspond à l'accélération de la pesanteur. On remarquera que ces mesures rx, ry et rz peuvent éventuel- lement être corrigées par des biais accélérométriques si ces derniers sont connus. On notera de plus que, pour l'ensemble des mesures à la fois accélérométriques et gyrométriques effectuées lors de la mise en oeuvre de la présente invention, on réalise de préférence un filtrage de ces mesures, par exemple pendant quelques secondes, pour s'affranchir d'éventuels bruits susceptibles d'avoir perturbé les capteurs utilisés. Ceci est possible en pratique puisque les calculs effectués à partir de ces mesures ne doivent pas être réalisés immé- diatement en temps réel, mais peuvent être effectuées sans problème quelques secondes après les mesures. Les angles e et 0 peuvent donc être déduits à partir d'au moins deux des mesures accélérométriques rx, ry, rz préci- sées cidessus. En revanche, l'angle de cap t, tel que défini précédemment, ne peut pas être obtenu à partir de ces mesures accéléromé- triques, puisqu'il est formé dans un plan (Nm, Em) qui est perpendiculaire à l'accélération de la pesanteur colinéaire à l'axe V. Aussi, selon l'invention, on détermine l'angle de cap i & l'aide d'un moyen de mesure approprié 3 de type connu, non décrit plus en détail dans la présente description. En plus dudit moyen de mesure 3, le dispositif 1 conforme à l'invention comporte un moyen de mesure 4 pour déterminer la latitude X du point P, ainsi que des calculateurs 5 et 6. Ledit calculateur 5 relié à la centrale inertielle 2 et auxdits moyens de mesure 3 et 4, respectivement par l'inter- médiaire de liaisons 7, 8 et 9, reçoit de ces derniers les valeurs des angles 8, ç et Tp, ainsi que la valeur de la latitude À. A partir de ces informations, ledit calculateur calcule la vitesse de rotation w cal de la terre T dans le référentiel de mesure R1 à partir de la relation: cose.cosi.w.co9X + sine.w.aini wcal = (cos.sine.sinO - sin4.cos9).w. cosX - cose.sino.w.sink (cosO.sine.cos9 + sinO.sinT).W.cosu - cos%.cose. w.sinX) Le calculateur 5 transmet ensuite le résultat obtenu au calculateur 6 par l'intermédiaire d'une liaison 10. Ce dernier calculateur 6, qui reçoit de plus, par l'intermé- diaire d'une liaison 11, la vitesse de rotation de la terre wmes mesurée par les gyromètres de la centrale inertielle 2, moes calcule les biais gyrométriques bo, be et bip des angles 0, e et T, en effectuant la soustraction entre mmes et w cal' Lesdits biais gyrométriques vérifient en effet la relation suivante, à des facteurs d'échelle gyrométriques près: mes - cal b y Le calculateur 6 peut transmettre le résultat obtenu par l'intermédiaire d'une liaison 12, par exemple à un système utilisateur non représenté ou à une mémoire également non représentée reliée à la centrale inertielle 2 et stockant ledit résultat, de manière à permettre la correction ulté- rieure de mesures gyrométriques de la centrale inertielle 2 au moyen des biais gyrométriques bo, be et b4 ainsi stockés. Dans un mode de réalisation préféré, les calculateurs 5 et 6 sont regroupés dans un calculateur unique 13. Ainsi, grâce à l'invention, on peut utiliser des gyromètres à performance moyenne et de faible coût, puisqu'on est capable de corriger les erreurs de mesure qui sont, de façon connue, constituées essentiellement des biais gyrométriques déterminés par le procédé conforme à l'invention. La présente invention est par conséquent particulièrement avantageuse: pour des aéronefs, par exemple des missiles, à centrale inertielle de faible coût, puisque les gyromètres sont les capteurs les plus coûteux de la centrale inertielle, ou -pour des aéronefs à distance et à temps de vol élevés puisque, dans ce cas, les erreurs de mesure même minimes sont susceptibles d'être très dommageables car les mesures doivent être répétées sur une longue durée. Des simulations ont par ailleurs montré que la précision des biais gyrométriques déterminés par le procédé conforme à l'invention dépend de la position initiale (immobile) de l'aéronef A, c'est-à-dire de l'angle d'assiette e (par exemple l'angle de la rampe de lancement d'un missile par rapport à l'horizontale), de l'angle de roulis 0 et de l'angle de cap i dudit aéronef A. Ces simulations ont montré que, pour obtenir la précision la plus élevée possible, il convient d'utiliser un angle d'assiette e important, par exemple 45 , et des angles de roulis * et de cap! les plus faibles possibles. On a pu vérifier en outre que les éventuelles erreurs de mesures des moyens de mesure 3 et 4 du type connu, dues à un défaut de précision et concernant donc l'angle de cap p et la latitude À, n'influent pas sur la précision des biais gyrométriques déterminés conformément5 à l'invention. REVENDICATIONS 1. Procédé pour estimer les biais gyrométriques de mesures effectuées par un ensemble de gyromètres comportant au moins un gyromètre et monté sur un aéronef (A), caractérisé en ce que: - on immobilise ledit aéronef (A) en un endroit déterminé (P) sur la terre (T); - on détermine les angles d'assiette (8), de roulis ( obtenir les valeurs estimées desdits biais gyrométriques. , 2. Procédé selon la revendication 1, caractérisé en ce que les angles d'assiette (e) et de roulis (4) sont déterminés à partir de mesures accélérométriques effectuées sur l'aéronef (A). * 3. Procédé selon la revendication 2, ledit référentiel de mesure (R1) correspondant à un référentiel OXYZ dans le- quel: - l'axe OX correspond à l'axe de l'aéronef (A), - l'axe OY est perpendiculaire audit axe OX suivant les ailes (Al, A2) dudit aéronef (A), et - l'axe OZ est perpendiculaire auxdits axes OX et OY, caractérisé en ce qu'on détermine les angles d'assiette e et de roulis ç à partir d'au moins deux des trois relations suivantes: rx = -g.sine. (l+fx) ry = g.cose.sinO.(l+fy) rz = g.cose.coso.(l+fz) dans lesquelles: - rx, ry et rz représentent des mesures accélérométriques effectuées par des accéléromètres respectivement selon les axes OX, OY et OZ, - g représente l'accélération de la pesanteur, et - fx, fy et fz représentent des facteurs d'échelle desdits accéléromètres. 4. Procédé selon l'une des revendications 2 ou 3, caractérisé en ce que lesdites mesures accélérométriques sont effectuées par des accéléromètres d'une centrale inertielle (2) montée sur l'aéronef (A). 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'on détermine l'angle de cap () à partir de mesures effectuées par un moyen de mesure appro- prié (3) externe à la centrale inertielle (2) de l'aéronef (A). 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'on calcule la vitesse de rotation wcal de la terre (T) dans le référentiel de mesure (R1) à partir de la relation: cose. cosT.W.cosA + sine.w.sink wcal = (cose.sine.sin0 - sin4.cos0).w.cosÀ cose.sin.w.sinX \(cosçb.sine.cos + sini.sin).w.cos À - co,.cosee.u.-inX dans laquelle: - À correspond à la latitude dudit point déterminé (P), - e, X et t représentent respectivement les angles d'as- siette, de roulis et de cap, et - w représente la vitesse de rotation de la terre (T) autour de son axe de rotation (CC'). À 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que ledit aéronef (A) est immobilisé à son poste de lancement. 8. Dispositif pour la mise en oeuvre du procédé spécifié sous l'une quelconque des revendications 1 à 7, ledit dispositif (1) étant associé audit ensemble de gyromètres comportant au moins un gyromètre et monté sur l'aéronef (A), caractérisé en ce qu'il est associé à un système auxiliaire susceptible de déterminer les angles d'assiette (e) et de roulis (0) de l'aéronef (A), et en ce qu'il comporte: - un premier moyen de mesure (3) déterminant la latitude (X) dudit endroit (P) déterminé de la terre (T); - un second moyen de mesure (4) déterminant l'angle de cap (*) dudit aéronef (A); - un premier calculateur (5) calculant la vitesse de rota- tion de la terre à partir d'informations reçues dudit système auxiliaire et desdits premier et second moyens de mesure (3, 4); et - un second calculateur (6) déterminant les biais gyrométri- ques à partir d'informations reçues dudit premier calcula- teur (5) et dudit ensemble de gyromètres. 9. Dispositif selon la revendication 8, caractérisé en ce que ledit système auxiliaire et ledit ensemble de gyromètres font partie de la centrale inertielle30 (2) de l'aéronef (A). 10. Dispositif selon la revendication 9, caractérisé en ce que les valeurs des biais gyrométriques estimées par ledit dispositif (1) sont stockées dans une mémoire intégrée dans ladite centrale inertielle (2). 11. Dispositif selon l'une quelconque des revendications 8 à , caractérisé en ce que lesdits premier et second calculateurs (5, 6) sont regroupés dans un calculateur unique (13).