Brevet FR3089287A1 — publié le 20200605

Description Titre de l'invention : Système AHRS hybridé comportant un dispositif de mesure de l'intégrité de l'attitude calculée [0001] Le domaine de l'invention est principalement celui de la navigation aéronautique et, en particulier, (les systèmes de mesures hybridés mettant en oeuvre à la fois des unités de mesure inertielles et des moyens de géolocalisation par satellite. [0002] Les systèmes dits « AHRS », acronyme signifiant « Attitude and Heading Reference System », sont des équipements intégrant des mesures inertielles d'accélérations et de rotations angulaires afin de déterminer les attitudes d'un porteur. On entend par attitude l'orientation du porteur dans l'espace. Cc porteur est généralement un aéronef. Dans cc cas, l'attitude est typiquement définie par les valeurs de roulis, de tangage et de lacet. [0003] En raison des défauts des mesures des différents senseurs, ces attitudes sont erronées et la quantité d'erreur est incrémentale. Dans certaines AHRS telles que celle représentée sur la figure I, on corrige les attitudes intégrées Ai obtenuesaprès intégration des vitesses angulaires VA en les comparant à des attitudes de référence AR qui sont généralement des attitudes dites accélérométriques obtenues par la projection de la gravité terrestre sur des axes horizontaux liés au porteur. L'écart entre les attitudes AHRS et les attitudes de référence est rebouclé avec un certain gain G. Sur cette figure et sur la figure 3, les intégrateurs sont représentés par le symbole de l'intégration figuré par un S allongé. [0004] Afin d'améliorer la précision du système, les mesures inertielles peuvent être fu- sionnées avec des mesures provenant d'un système de navigation par satellites, système dit « GNSS », acronyme signifiant « Global Navigation Satellite System ». [0005] Dans le contexte aéronautique, les données concernant l'attitude de l'appareil sont exigées avec une quantification en temps réel de leur intégrité, c'est-à-dire de la probabilité que l'erreur de mesure soit supérieure à un seuil déterminé prenant en compte les différentes causes de défaillances possibles. Par exemple, une intégrité inférieure à 107par heure est fréquemment exigée pour des mesures d'attitudes d'un aéronef. [0006] Les attitudes issues des AHRS sont fournies sans quantification de leur intégrité. Il est cependant possible d'estimer cette intégrité sous réserve d'utiliser un algorithme voisin de celui des centrales inertielles hautes performances. Cette méthode impose : - Une forte complexité et une charge de calcul élevée, l'intégrité des centrales inertielles étant estimée grâce à l'emploi en parallèle de plusieurs filtres de Kalman. La demande européenne de référence EP 3 018 447 intitulée « Procédé de détermination d'un rayon de protection associé à un paramètre de navigation d'un système hybride de 2 navigation inertielle et système associé » propose une première solution de ce type. Lc brevet US 9 341 718 intitulé « Method and System for providing integrity for hybrid attitude and truc hcading » propose également une solution de ce type. - Des contraintes de performances sur les senseurs inertiels des centrales AHRS, afin de disposer d'une précision autorisant l'utilisation de filtres de Kal man. A titre d'exemple, la dérive des gyromètres doit être inférieure à dix degrés par heure. [0007] Le système hybride selon l'invention ne présente pas les inconvénients précédents. En effet, l'estimation d'attitude est réalisée de manière conventionnelle. Des positions AHRS sont estimées par triple intégration des mesures inertielles AHRS. Ces positions hybrides AHRS sont comparées avec les positions GNSS pour en déduire des corrections. Ces corrections sont appliquées sur les données hybrides. [0008] L'invention porte sur la production d'une mesure de l'intégrité de ces attitudes. Dans le système selon l'invention, l'intégrité est estimée grâce à la connaissance des fonctions de transfert du filtre d'hybridation, de la valeur du rayon de protection associé à la position donnée par le récepteur du système de positionnement par satellites et par le calcul d'une constante majorante issue de la connaissance du système AHRS hybridé. Le rayon de protection est calculé pour une intégrité de 107/h. [0009] Plus précisément, l'invention a pour objet un système de mesure d'attitude d'un porteur, ledit système étant hybridé et comportant au moins : - un système de positionnement par satellites ; - une unité de mesure comportant trois gyromètres et trois accéléromètres ; - une électronique de calcul des informations d'attitude du porteur à partir des informations issues du système de positionnement par satellites et de l'unité de mesure, caractérisé en ce que l'électronique de calcul comporte un moyen de calcul d'un rayon angulaire de protection, c'est-à-dire un rayon angulaire tel que le risque que l'erreur entre l'attitude calculée et l'attitude vraie du porteur ne soit pas comprise à l'intérieur de ce rayon est inférieur à une probabilité donnée, ledit rayon étant égal à la somme de deux contributions, la première contribution étant égale aux erreurs d'attitudes liées à l'unité de mesure, la seconde contribution étant égale aux erreurs induites sur les estimations d'attitude par les erreurs de position dues au système de positionnement par satellite, en prenant en compte l'hypothèse d'une panne du système de positionnement par satellites. [0010] Avantageusement, la première contribution est égale à l'erreur de positionnement maximale obtenue à partir de simulations de Monte-Carlo prenant en compte l'ensemble des erreurs possibles de l'unité de mesure. [0011] Avantageusement, la seconde contribution est égale au produit : [0012] Du gain d'une fonction de transfert de l'électronique de calcul de l'attitude du porteur, ladite fonction de transfert étant déterminée pour une fréquence d'oscillation 3 de la position donnée par le système de positionnement par satellites ; [0013] Par la valeur du rayon de protection du système de positionnement par satellites, ladite valeur étant connue sous la terminologie « HIL », acronyme de « Horizontal Integrity Li mit ». [0014] Avantageusement, l'électronique de calcul de l'attitude comporte trois intégrateurs en série et un filtre correcteur à gains constants, les données d'entrée du premier intégrateur étant les informations de vitesses angulaires du porteur et les données d'entrée du deuxième intégrateur étant les accélérations du porteur et le résultat du premier intégrateur, les données d'entrée du troisième intégrateur étant le résultat du deuxième intégrateur initialisées par les informations de position données par le système de positionnement par satellites, les données de sortie du premier intégrateur étant les informations d'attitude du porteur dite attitude hybridée, les données de sortie du second intégrateur étant les informations de vitesse du porteur dite vitesse hybridée et les données de sortie du troisième intégrateur étant les informations de position du porteur dite position hybridée. [0015] Avantageusement, le rayon angulaire de protection est inférieur à 1 degré. [0016] Avantageusement, l'intégrité associée au rayon de protection angulaire est comprise entre 10 'heure d'utilisation et 10 9/heure d'utilisation. [0017] Avantageusement, le porteur est un aéronef [0018] L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description qui va suivre donnée à titre non limitatif et grâce aux figures annexées parmi lesquelles : [0019] [fig.1] précédemment commentée, représente une centrale AHRS selon l'art antérieur ; [0020] [fig.2] représente un système de mesure d'attitude d'un porteur selon l'invention ; [0021] [fig.3] représente une électronique de calcul des informations d'attitude selon l'invention. [0022] Le système de mesure d'attitude selon l'invention est, comme il a été dit, un système hybride. Il est représenté sur la figure 2. Il est intégré à un porteur, non représenté sur cette figure. Il comporte au moins : - un système 10 de positionnement par satellites GNSS ; - une unité de mesure inertielle « UMI » 20 comportant des gyromètres et des accéléromètres ; - une électronique 30 de calcul des informations hybridées d'attitude A, de vitesse V et de positionnement P du porteur à partir des informations issues du système de positionnement par satellites 10 et de l'unité de mesure 20. [0023] Le système de mesure fonctionne en temps réel. [0024] Les systèmes de positionnement par satellites 10 sont connus de l'homme du métier. Ils fournissent en temps réel les informations de position du porteur. De la même façon, les unités de mesure de la vitesse angulaire et de l'accélération 20 sont connues de l'homme du métier. [0025] Le système de mesure utilise les incréments mesurés par l'unité de mesure inertielle UMI 20, un cap de référence et les sorties conventionnelles du récepteur GNSS 10. La plateforme de calcul temps réel 30, intégrée dans le dispositif réalise les calculs des données hybrides d'attitude A, de vitesse V et de position P. L'estimation d'attitude est réalisée de manière conventionnelle. Les positions AHRS sont estimées par triple intégration des mesures inertielles AHRS. Ces positions hybrides AHRS sont comparées avec les positions GNSS pour en déduire des corrections. Ces corrections sont appliquées sur les données hybrides au moyen de filtres correcteurs 34. [0026] Plus précisément, le calcul des attitudes AHRS hybrides est réalisé dans la plateforme de calcul interne à l'équipement 30 qui est représentée sur la figure 3. Elle comporte une chaîne de trois intégrateurs 31, 32 et 33. L'intégrateur 31 est un intégrateur des vitesses angulaires VA issues de l'unité de mesure inertielle 20. Il fournit des attitudes hybrides AH. L'intégrateur 32 est un intégrateur d'attitudes. Il fournit des vitesses hybrides VH. L'intégrateur 33 est un intégrateur des positions. Il fournit des positions hybrides PH. Le fonctionnement de ces intégrateurs s'apparente à celui des centrales inertielles. La plateforme de calcul inclut également une fonction de calcul de l'accélération de la pesanteur qui est ajoutée aux incréments accélérométriques intégrés, ainsi qu'une fonction de correction de l'accélération de Coriolis, pour prendre en compte le fait que l'intégration des accélérations se fait dans un référentiel non galiléen. Ces fonctions ne sont pas représentées sur la figure 3. [0027] La mécanisation de la plateforme peut être réalisée en axes géographiques ou en axes liés à la plateforme. Les axes géographiques sont, par exemple, référencés Nord, Est, Down ou « NED ». [0028] Chaque position hybride obtenue par intégration des données de l'UMI 20 est comparée à la même position obtenue par le système de positionnement 10 GNSS pour en déduire un écart. Les écarts sont obtenus en trois dimensions selon trois axes, par exemple sous forme d'écarts de latitude, de longitude et d'altitude dans le cas d'une mécanisation NED. [0029] A partir de ces écarts, des corrections sont élaborées par des filtres correcteurs 34 à gains constants, intégrés à la plateforme de calcul. Les gains sont ajustables en fonction du temps de réponse et de la robustesse souhaitée. Des gains élevés conduisent à un temps de réponse court mais une plus grande sensibilité aux erreurs de position GNSS. Des gains petits conduisent à une grande robustesse aux erreurs de position GNSS. [0030] Chaque filtre correcteur est dédié à une voie correspondant à un axe particulier. Par exemple, la première voie correspond à un écart de latitude, la seconde voie à un écart de longitude et la troisième voie à un écart d'altitude. Pour les écarts de latitude ou de longitude, le filtre correcteur produit des corrections dc position, de vitesse, d'attitude et de précession. Pour les écarts d'altitude, le filtre correcteur produit uniquement des écarts d'altitude et de vitesse verticale. [0031] Les gains Ki des filtres correcteurs sont calculés en établissant la fonction de transfert analytique du système. Ils sont définis en fonction de la constante dc temps r du système. Par exemple, dans le cas du filtre correcteur de latitude, les gains sont les suivants : [0032] [Math .1] t3- [0033] [Math.2] K [0034] Rt étant le rayon terrestre, g étant l'accélération de la pesanteur et étant le co- efficient d'amortissement.. Ce coefficient est voisin del. [0035] La valeur de t dimensionne la réponse du système. Une constante de temps r petite conduit à des corrections importantes et un meilleur suivi de l'évolution des écarts de position. Une constante de temps t élevée conduit à des corrections plus faibles et un lissage plus important de l'évolution de l'écart de position. [0036] Un des points délicats de ce type (le système est l'estimation de l'erreur (les données de sortie. Il est cependant important, dans un certain nombre d'applications, de fournir les données d'attitude avec une intégrité connue correspondant à la probabilité que l'erreur de mesure soit supérieure à un seuil déterminé, en prenant en compte les causes de défaillances possibles. [0037] L'objet de l'invention est de borner ces erreurs de mesure de façon à garantir une intégrité donnée. Dans le système selon l'invention, l'intégrité des attitudes est assurée grâce à la production d'un rayon de protection, c'est-à-dire un seuil tel que l'erreur entre l'attitude vraie et l'attitude estimée n'est supérieure à ce seuil qu'en dessous d'une probabilité donnée. Dans un contexte aéronautique, selon les applications, cette probabilité varie entre 10-5/heure d'utilisation et 109/heure d'utilisation. [0038] Ce seuil doit prendre en compte l'ensemble des contributions d'erreurs sur les attitudes AHRS hybrides. On identifie dans ce contexte : - Les erreurs ducs aux capteurs inertiels ; 6 - Les erreurs introduites par les asservissements retenus ; - Les erreurs ducs aux entrées GNSS utilisées pour corriger les attitudes hybrides. Il faut prendre en compte les erreurs « courantes » liées à la précision du GNSS et les erreurs de type « panne du système GNSS ». [0039] On modélise le rayon de protection HIL des attitudes AHRS hybrides sous la forme [0040] [Math.3] HZ H HO' 1-Hfl [0041] Avec : - HO : contribution des erreurs liée aux défauts des capteurs et à l'algorithme. - H 1 : contribution liées aux erreurs GNSS, prenant en compte l'hypothèse d'une panne GNSS. [0042] En raison de la complexité du système, de la variété des comportements des erreurs capteurs, ou des cas d'utilisation du dispositif, un calcul analytique de HO,,, n'est pas réalisable. On l'estime donc à partir de simulations de Monte Carlo sur un modèle logiciel du système. Pour être couvrant vis-à-vis de l'ensemble des combinaisons d'erreurs dues aux capteurs et des différents cas d'utilisation, HO,,,N est choisi égal à la valeur d'erreur d'attitude maximale correspondant au pire cas d'erreur identifié lors des simulations. A titre d'exemple, le nombre de simulations nécessaires à l'estimation de cette erreur peut être de 106. [0043] Le calcul de Hl peut être fait analytiquement. Les attitudes, vitesses et positions hybrides sont calculées par intégration successives des incréments inertiels, puis corrigées grâce au filtre correcteur à gains constants. En conséquence, l'ensemble du système est parfaitement connu et modélisable. Grâce à cette connaissance, il est possible de calculer la fonction de transfert du système H(p). [0044] [Math.4] S(P) 11(P) E( [0045] La fonction de transfert H(p) permet de décrire la sortie S(p) en fonction d'une entrée déterminée E(p). Le module de H(p) donne en particulier le coefficient multiplicateur entre l'amplitude d'une perturbation en entrée et la résultante en sortie. Les perturbations en entrée peuvent être sur les incréments inertiels accélérométriques ou gyrométriques et/ou des erreurs sur la position GNSS utilisée en référence. [0046] Dans le cas de l'AHRS hybridée GNSS, on peut considérer la sortie « erreur des attitudes hybrides » en fonction de l'entrée « erreur sur la position GNSS ». En utilisant le module de cette fonction de transfert, pour une erreur de position ErrGNSS oscillant à une fréquence f, on peut calculer l'impact sur l'erreur d'attitude ErrATT à la même fréquence. En notant H E . la fonction de transfert, on a alors la 7 relation : [0047] [Math.5] Err AMI) = H r'.0.2,3_,,,e(f)' F.r7G,VSSCO [0048] L'erreur sur la position GNSS peut revêtir différentes natures. Elle peut être une erreur « normale et rare » correspondant aux erreurs classiques du GNSS, une erreur oscillante ou une erreur liée à une panne du système GNSS. Dans tous les cas, la position GNSS est fournie avec une mesure de son intégrité ou « Horizontal Integrity Limit» HIL. Le système GNSS garantit que sous une probabilité de panne donnée, l'erreur sur la position GNSS est en dessous de ce rayon. Au-dessus de ce rayon, le système GNSS se déclare en panne. [0049] On utilise alors la fonction de transfert pour calculer la réponse du L 1/4,fi système AHRS hybridé avec en entrée la valeur de HILGNss. La sortie est alors le rayon d'intégrité des attitudes hybrides Rpmco en faisant l'hypothèse que les erreurs GNSS oscillent à unc fréquence donnée. [0050] [Math.6] 14.ksstf) [0051] Pour borner cette erreur, on choisir la fréquence induisant la plus forte réponse, c'est-à-dire le maximum de la fonction de transfert. On peut alors écrire [0052] [Math.7] [0053] Le rayon de protection des attitudes hybrides AHRS est donc la somme de ces deux contributions. [0054] [Math.8] = -fmax(lHenecs, [0055] L'ensemble des calculs nécessaires à la détermination des deux contributions du rayon de protection HIL ne pose pas de problèmes particuliers pour l'homme du métier et nécessitent uniquement des adaptations des logiciels d'une centrale AHRS. [0056] Lc système de mesure d'attitude d'un porteur selon l'invention présente de nombreux avantages détaillés ci-dessous. [0057] Lc système peut s'adapter à tout type de senseurs inertiels de classe AHRS. En par- ticulier, le calcul du rayon de protection peut se faire dans un large spectre de dérives des centrales inertiels, compris entre 100°/h et 0.1°/h. La perte de performance des senseurs inertiels n'influence que peu la performance de la méthode, le HIL estimé étant dépendant essentiellement du HIL GNSS. 8 [0058] Le système selon l'invention permet de calculer des attitudes AHRS intègres avec un rayon de protection de l'ordre du degré pour une intégrité de 107/h à partir de constituants standards existant qui sont une centrale AHRS, un récepteur GNSS, un calculateur électronique. [0059] La méthode de calcul des attitudes AHRS est simple. Elle se hase sur un unique filtre à gains constants. [0060] La méthode de calcul de l'intégrité est également simple et peu conteuse en termes de charge de calcul temps réel. Le calcul d'intégrité dans les systèmes inertiels est généralement traité à l'aide d'algorithmes complexes utilisant plusieurs filtres de Kalman en parallèle. Dans le système selon l'invention, on utilise uniquement la connaissance de la fonction de transfert du filtre à gain constant. Il suffit de multiplier le HIL GNSS par une constante, calculée analytiquement à la conception du filtre correcteur et d'ajouter une constante. Cette solution peut donc être embarquée dans des systèmes peu complexes comportant des capacités de calcul limitées. [0061] Le rayon de protection est majorant. En effet, les hypothèses d'intégrité réalisées sont majorantes par rapport aux conditions réelles dans la mesure où l'on retient le cas le plus défavorable. Ceci assure une grande pertinence à l'intégrité estimée. [0062] Au final, le système selon l'invention est particulièrement bien adaptée à une uti- lisation aéronautique. Dans ce cas, l'architecture d'implantation comprend les éléments suivants : - Une centrale inertielle aéronautique de type AHRS ; - Un récepteur GNSS aéronautique ; - Un calculateur, qui peut être intégré à l'équipement AHRS. 9 [Revendication 1] [Revendication 2] [Revendication 3] [Revendication 4]