FR 2988836 A1 20131004 FR 1200925 A 20120328 Procédé de détermination d'une masse estimée d'un aéronef et système correspondant La présente invention concerne un procédé de détermination d'une masse estimée d'un aéronef à un instant de détermination choisi. Le pilotage d'un aéronef repose sur la connaissance de paramètres de vol de celui-ci, telles que sa vitesse par rapport à l'air ambiant, son altitude et son incidence. Ces paramètres sont déterminés au moyen de sondes situées sur la carlingue de l'aéronef. De manière connue, ces sondes comprennent des capteurs de pression statique, des sondes de Pitot pour la mesure de la pression totale, des sondes d'incidence montées sur un dispositif à palette ou pneumatique, et des sondes de température totale. Ces sondes sont alors raccordées à des moyens de détermination de la grandeur correspondante. Notamment, un anémomètre détermine la vitesse de l'aéronef par rapport à l'air à partir des mesures de pressions totale et statique, et un altimètre détermine l'altitude de l'aéronef à partir de mesures de la pression statique. Ces mesures sont ensuite regroupées et affichées sur un dispositif d'affichage qui constitue une source d'information centrale à partir de laquelle le pilotage de l'aéronef est effectué. De manière connue, les sondes d'incidence et de pression se présentent sous la forme de palettes et de tubes protubérants de la peau de l'aéronef. Elles sont ainsi exposées à des facteurs météorologiques ou mécaniques qui peuvent en altérer le fonctionnement, en particulier par obstruction d'orifices de ces sondes par du givre ou des poussières ou insectes, ou par blocage des dispositifs à palette. De telles défaillances conduisent à la production de mesures erronées, et notamment à l'affichage d'incidences, de vitesses et/ou d'altitudes fausses, pouvant conduire le pilote à effectuer des manoeuvres inappropriées. Par exemple, de fausses informations de pilotage peuvent conduire au décrochage de l'aéronef ou à la perte de son contrôle du fait d'une vitesse trop élevée. Afin de minimiser les conséquences de telles défaillances, les régulations en aéronautique imposent aux constructeurs de disposer de moyens de mesures de ces caractéristiques redondants. Ainsi, les aéronefs comprennent généralement au moins une sonde de secours identique à chaque sonde susceptible de faillir. Toutefois, cette solution ne donne pas entière satisfaction. En effet, les sondes de secours existantes sont pour la plupart de type protubérant, et par conséquent présentent les mêmes risques de défaillance que les sondes qu'elles sont destinées à suppléer. Ainsi, en cas de défaillance des sondes de mesure, aucune mesure fiable n'est fournie au pilote. En outre, le pilote ne dispose d'aucun moyen pour vérifier la fiabilité des informations fournies par les sondes et ne dispose d'aucune information alternative fiable. Notamment, le pilote ne connaît pas d'estimation fiable de la masse de l'aéronef en vol, estimation qui permettrait le calcul de caractéristiques de vol de l'aéronef indépendamment de mesures de capteurs défaillants. Un but de l'invention est donc de déterminer une masse estimée de l'aéronef, qui soit fiable et précise. A cette fin, l'invention a pour objet un procédé du type précité, caractérisé en ce qu'il comprend : - la détermination d'une première masse de l'aéronef, à partir de caractéristiques de l'aéronef déterminées avant ou pendant le décollage de l'aéronef, - la détermination, lors du vol de l'aéronef, d'une deuxième masse de l'aéronef, à partir d'une équation de sustentation de l'aéronef exprimant la masse de l'aéronef en fonction d'une information représentative du facteur de charge de l'aéronef, d'un coefficient de portance de l'aéronef, d'une information de vitesse de l'aéronef et d'une pression statique d'une masse d'air traversée par l'aéronef, déterminées à partir de mesures de capteurs de l'aéronef lors du vol de l'aéronef, - l'évaluation de la masse estimée audit instant de détermination, en fonction desdites première et deuxième masses. Suivant des modes particuliers de réalisation, le procédé de détermination comporte l'une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : - lesdites première et deuxième masses sont des masses de l'aéronef à un instant initial au parking, et l'étape d'évaluation de ladite masse estimée comprend la détermination d'une masse estimée au parking audit instant initial en fonction desdites première et deuxième masses, la détermination d'une masse de carburant consommé entre ledit instant initial et ledit instant de détermination, ladite masse estimée étant évaluée en soustrayant ladite masse de carburant consommé de ladite masse estimée au parking ; - la détermination de ladite première masse comprend l'estimation d'une première masse initiale de l'aéronef, réalisée avant le décollage, l'estimation d'une deuxième et d'une troisième masses initiales de l'aéronef à partir de caractéristiques de vol de l'aéronef déterminées lors du décollage, ladite première masse étant déterminée en fonction desdites première, deuxième et troisième masses initiales ; - ladite première masse est déterminée à partir d'une médiane desdites première, deuxième et troisième masses initiales ; - la détermination de la première masse de l'aéronef comprend l'estimation d'une première masse initiale de l'aéronef à partir d'une estimation d'une masse de base de l'aéronef, d'une masse de carburant de l'aéronef au parking et d'une masse du chargement de l'aéronef ; - la détermination de la première masse de l'aéronef comprend l'estimation d'une deuxième masse initiale de l'aéronef, ladite deuxième masse initiale étant déduite à partir d'une valeur de poussée des moteurs, d'une valeur de traînée de l'aéronef et d'une valeur d'accélération de l'aéronef déterminées pendant le décollage de l'aéronef, à un instant choisi, notamment 3 secondes après l'envol de l'aéronef ; - la détermination de la première masse de l'aéronef comprend l'estimation d'une troisième masse initiale de l'aéronef, à partir d'une équation de sustentation de l'aéronef exprimant la masse de l'aéronef en fonction d'une information représentative du facteur de charge de l'aéronef, d'un coefficient de portance de l'aéronef, d'une information de vitesse de l'aéronef et d'une pression statique d'une masse d'air traversée par l'aéronef, déterminées pendant le décollage de l'aéronef, à un instant choisi, notamment 3 secondes après l'envol de l'aéronef ; - la détermination de ladite deuxième masse de l'aéronef comprend : - la détermination d'une pluralité de masses instantanées de l'aéronef à des instants de vol distincts, lesdites masses instantanées étant déterminées à partir de ladite équation de sustentation de l'aéronef, - la détermination à partir desdites masses instantanées et de masses de carburant consommé entre un instant initial et lesdits instants de vol, d'une pluralité de quatrièmes masses initiales de l'aéronef audit instant initial, - la détermination, à partir desdites quatrièmes masses initiales, de ladite deuxième masse de l'aéronef ; - la détermination de ladite deuxième masse de l'aéronef comporte le calcul d'une moyenne pondérée desdites quatrièmes masses initiales en attribuant à chaque quatrième masse initiales un coefficient de pondération fonction d'une crédibilité accordée à la valeur de cette masses initiale, sous la forme : t' fx,p(t)m4i(t) msp t2 où trisp(tv) désigne la masse de l'aéronef à l'instant initial déterminée à l'instant de détermination tv, chaque m4'(t) désigne une quatrième masse initiale déterminée à un instant t, et chaque x,,,(t) est un coefficient de pondération attribué à la quatrième masse initiale m4'(t) ; - le coefficient de pondération attribué à une quatrième masse initiale déterminée à un instant t est égal à un lorsqu'à l'instant t, l'aéronef a sensiblement un vol rectiligne horizontal selon des critères prédéterminés, ou à zéro dans le cas contraire ; - ladite masse estimée est évaluée à partir d'une moyenne pondérée desdites première et deuxième masses ; - ladite moyenne pondérée desdites première et deuxième masses est réalisée en attribuant à chacune desdites première et deuxième masses un coefficient de pondération, le coefficient de pondération attribué à ladite deuxième masse étant proportionnel à la durée de vol rectiligne uniforme de l'aéronef depuis un instant initial ; - le coefficient de pondération attribué à ladite deuxième masse est exprimé sous la forme : fx.(t) p(tv). _t1 où p(tv) désigne le coefficient de pondération attribué à ladite deuxième masse à l'instant de détermination tv, t1 est l'instant initial, etx,(t) une fonction égale à un lorsqu'à l'instant t, l'aéronef a sensiblement un vol rectiligne uniforme, ou à zéro dans le cas contraire. L'invention a également pour objet un système de détermination d'une masse estimée d'un aéronef à un instant de détermination choisi, caractérisé en ce qu'il comprend : - des moyens configurés pour déterminer une première masse de l'aéronef, à partir de caractéristiques de l'aéronef déterminées avant ou pendant le décollage de l'aéronef, - des moyens configurés pour déterminer, lors du vol de l'aéronef, une deuxième masse de l'aéronef, à partir d'une équation de sustentation de l'aéronef exprimant la masse de l'aéronef en fonction d'une information représentative du facteur de charge de tv fx.(t) t2 l'aéronef, d'un coefficient de portance de l'aéronef, d'une information de vitesse de l'aéronef et d'une pression statique d'une masse d'air traversée par l'aéronef, déterminées à partir de mesures de capteurs de l'aéronef lors du vol de l'aéronef, - des moyens configurés pour évaluer la masse estimée audit instant de détermination, en fonction desdites première et deuxième masses. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins suivants, sur lesquels : - la figure 1 illustre de manière schématique un aéronef dans des conditions de vol classiques, auquel est appliquée l'invention ; - la figure 2 est un schéma synoptique général représentant le système selon l'invention ; - la figure 3 est un schéma synoptique illustrant le procédé selon l'invention ; - la figure 4 est un schéma synoptique illustrant les étapes mises en oeuvre pour évaluer la fiabilité de mesures d'incidence selon un mode de réalisation ; - la figure 5 est un schéma synoptique représentant des étapes mises en oeuvre pour évaluer la fiabilité de mesures de pression selon un mode de réalisation ; - la figure 6 illustre un mode de représentation par un dispositif d'affichage auxiliaire d'un premier ensemble d'informations relatives au vol ; - la figure 7 illustre un mode de représentation par un dispositif d'affichage auxiliaire d'un deuxième ensemble d'informations relatives au vol ; - la figure 8 illustre un mode de représentation par un dispositif d'affichage auxiliaire d'un troisième ensemble d'informations relatives au vol. La Figure 1 représente, de manière schématique, un aéronef 1 en vol auquel est appliqué le procédé selon l'invention. L'aéronef 1 est représenté sur la Figure 1 par son seul centre de gravité. Son axe longitudinal est orienté selon un axe X, qui forme avec l'horizontale A un angle 0 appelé assiette de l'aéronef. Il se déplace par rapport à l'air selon un vecteur vitesse Vair , qui forme avec l'horizontale A un angle vair appelé pente air de l'aéronef. L'angle a entre l'axe longitudinal X de l'aéronef 1 et son vecteur vitesse est appelé angle d'incidence. Ces angles vérifient ainsi la relation : 0 = a+ ya,,. . La vitesse propre Vp de l'aéronef, qui est sa vitesse par rapport à l'air dans un plan horizontal, est liée à sa vitesse Vs par rapport au sol dans ce plan horizontal par le ...--^-- triangle des vitesses selon la relation : Vp = Vs - W , où W désigne le vecteur vitesse du vent dans le plan horizontal. Le système 2 de détermination de paramètres de l'aéronef 1 selon l'invention est représenté de manière schématique sur la Figure 2. Le système 2 est propre à déterminer des caractéristiques de l'aéronef 1 ou de l'air ambiant lors du vol, et à évaluer la fiabilité des mesures de capteurs de l'aéronef 1. Ce système 2 comprend un calculateur central 3 muni d'une mémoire 4. La mémoire 4 comprend notamment des tables de correspondance et des abaques fournissant la valeur d'une caractéristique de vol de l'aéronef en fonction de la valeur d'une ou plusieurs autres caractéristiques. En particulier, la mémoire 4 comprend une table d'atmosphère type donnant l'altitude pression de l'aéronef, notée Zp, en fonction de la pression statique autour de l'aéronef. Cette altitude pression est définie comme l'altitude qu'aurait un aéronef dans une atmosphère standard s'il était à cette pression statique. La mémoire 4 comprend également une table de correspondance fournissant la valeur d'un coefficient de portance de l'aéronef à partir de valeurs du nombre de Mach M et d'incidence a de l'aéronef. Cette table de correspondance est issue d'essais en vol préalablement réalisés. La mémoire 4 comprend par ailleurs des abaques, telles qu'un abaque permettant de déterminer la masse de l'aéronef 1 en fonction de la portance et de la traînée pendant la phase initiale de décollage, et de l'accélération en résultant. Ce système 2 comprend par ailleurs plusieurs capteurs 5, et notamment un capteur de pression statique 5a, un capteur de pression totale 5b, un capteur de température 5c et un capteur d'incidence 5d. Le capteur de pression statique 5a est propre à mesurer la pression statique Ps, c'est-à-dire la pression atmosphérique au niveau de l'aéronef. Le capteur de pression totale 5b est par exemple une sonde Pitot. Il est propre à mesurer la pression totale PT, somme de la pression dynamique Pdyn et de la pression statique Ps. Le capteur de température 5c est propre à mesurer la température totale, notée TAT pour « Total Air Temperature », correspondant à la température d'impact de l'air dans la sonde destinée à la mesure. Cette température est constante sur tout le filet d'air qui arrive sur la sonde. La température totale est liée à la température statique de l'air, notée SAT pour « Static Air Temperature », et correspondant à la température que mesurerait un thermomètre à l'arrêt dans la masse d'air, par la relation : (1) SAT = TAT TAT H -1 A42) (1+0.2M2) 2 où M désigne le nombre de Mach de l'aéronef et y 1.4 est le coefficient adiabatique de l'air. Le capteur d'incidence 5d est propre à déterminer l'incidence a de l'aéronef 1. Tous les capteurs sont reliés au calculateur 3. Le système 2 de détermination comprend également un altimètre 6, un anémomètre 7, un machmètre 8 et un accéléromètre 9. L'altimètre 6 est relié au capteur de pression statique 5a et au calculateur 3. Il est propre à déterminer l'altitude pression Zp de l'aéronef 1 par rapport à un niveau de référence à partir de la mesure de la pression statique Ps. A cette fin, l'altimètre 6 utilise par exemple une table d'atmosphère type comprenant des valeurs tabulées de pressions statiques en fonction de l'altitude. L'anémomètre 7 est relié aux capteurs de pression statique 5a et totale 5b et au calculateur 3. Il est propre à déterminer, à partir des pressions statique Ps et totale PT, la pression dynamique Pdyn = PT - Ps , et à déduire de cette pression dynamique la vitesse conventionnelle Vc de l'aéronef par rapport à l'air, à partir de la relation : 3,5 Pdyn ' v 2 -1 (2) 101325 1+ 0,2c \ 661,471 où 101325 correspond à la pression atmosphérique au sol, 661,471 est la vitesse du son au sol en noeuds. La vitesse conventionnelle est ainsi donnée par : = 661,471 Cette vitesse conventionnelles Vc est celle qui aurait produit la même pression dynamique Pdyn en volant dans l'atmosphère standard au sol. L'équivalent de vitesse EV peut être déduit de l'altitude pression Zp et du nombre de Mach M par : 2 1 poEV2 = -1 pa- M 2 2 2 12 Pd" +1\ 3'5 -1 101325 0,2 (3) (4) où p0 désigne la masse volumique de l'air au sol, a la vitesse du son et p la masse volumique de l'air . Lorsque le nombre de Mach de l'aéronef est faible (M Le calculateur 3 est propre à déterminer, à partir d'une valeur de Mach M de l'aéronef, sa vitesse par rapport à l'air, appelée vitesse vraie et notée TAS pour « True Air Speed », suivant la relation : TAS = Ae'lyR - SA T (5) où y est le coefficient adiabatique de l'air et R la constante universelle des gaz parfaits. Le calculateur 3 est également propre à déterminer la norme de la vitesse propre Vp de l'aéronef, norme de la projection horizontale de la vitesse vraie TAS, à partir de la relation : TAS2 = +Vz2(GPS) (6) Le système 2 de détermination comprend par ailleurs des moyens 11 de détermination de la quantité de carburant contenue dans les réservoirs de l'aéronef 1. Les moyens 11 sont reliés au calculateur 3. Ces moyens comprennent par exemple des jauges de réservoirs et des débitmètres. Les jauges de réservoirs sont propres à mesurer la quantité, notamment la masse de carburant contenue dans chaque réservoir, permettant au calculateur 3 de déterminer la masse de carburant dans l'aéronef 1. Les débitmètres sont propres à mesurer le débit massique de carburant fourni à chaque moteur, permettant au calculateur 3 de déduire la masse FU de carburant consommée, et à partir d'une mesure initiale de la masse de carburant, la masse de carburant restant. L'aéronef 1 comprend des moyens 12 d'interface homme-machine. Ces moyens 12 sont reliés au système 2 de détermination. Ils sont propres à présenter des informations à destination de l'équipage et notamment du pilote, et à recevoir des instructions du pilote destinées en particulier au système 2 de détermination. L'aéronef 1 comprend ainsi des instruments de bord, propres à présenter au pilote des informations relatives au vol de l'aéronef 1. Notamment, ces instruments comprennent des moyens d'affichage classiques, propres à afficher les caractéristiques du vol issues des mesures de pression et d'incidence par les capteurs de pression et d'incidence. L'aéronef 1 comprend également un dispositif d'affichage auxiliaire 14. Le dispositif d'affichage auxiliaire 14 comprend des moyens 15 d'affichage sélectif de différentes évaluations de caractéristiques du vol, en fonction d'un état de crédibilité affecté aux mesures des capteurs, et notamment des capteurs d'incidence et de pression. Le dispositif d'affichage auxiliaire 14 est complété en outre par des moyens 16 d'affichage de messages d'alerte destinés à informer le pilote de l'état de crédibilité affecté aux mesures des capteurs, et en particulier à l'avertir lorsque les mesures d'un ou plusieurs capteurs ne sont pas fiables. Ces messages peuvent être des signaux visuels et/ou sonores, par exemple des signaux lumineux. Avantageusement, ces signaux sont des messages textuels ou symboliques et sont inclus dans le dispositif d'affichage des pannes de l'avionique résidente. L'aéronef 1 comprend par ailleurs une interface 17 de saisie, par exemple des boutons de commande et un clavier, permettant au pilote de donner des instructions au système de détermination ou d'entrer des valeurs numériques de paramètres de vol. On a représenté sur la Figure 3 un exemple de mise en oeuvre du procédé selon un mode de réalisation de l'invention, pour la surveillance de la crédibilité d'informations de pilotage fournies par les sondes et appareils de mesure de cet aéronef, et en particulier des valeurs de pression et d'incidence mesurées. Le procédé comprend une étape 21 de détermination d'une masse au parking initiale mpio de l'aéronef 1, en vue de déterminer une masse estimée de l'aéronef 1 à tout instant lors de son vol. Cette masse au parking initiale mpio est élaborée en estimant, selon différentes méthodes, plusieurs masses mo au parking de l'aéronef 1 et en déterminant la masse au parking initiale mpio de l'aéronef 1 à partir des estimations mpi ainsi obtenues. L'étape 21 comprend ainsi avantageusement la détermination de trois masses au parking mpl, ma, ma, et le choix comme valeur de masse au parking initiale mpio de la valeur médiane de ces trois masses (10= 1, 2 ou 3). L'étape 21 comprend une phase 23 de détermination de la première masse au parking mpi de l'aéronef 1, à un instant initial t1 avant le décollage de l'aéronef 1. Cette masse au parking est déterminée par estimation et sommation de la masse de base de l'aéronef 1, la masse de carburant et la masse du chargement. Cette masse n'est généralement pas exacte. Notamment, la mesure de la masse de carburant contenue dans les réservoirs n'est pas une mesure exacte. On estime l'erreur faite sur cette masse mpi à quelques pourcents. En outre, la masse de base et la masse du chargement peuvent subir des erreurs dues à la mesure, à la transmission d'informations, au pilote et/ou à la saisie dans le système. L'étape 21 comprend une phase 25 de détermination de la deuxième masse au parking mp2 de l'aéronef 1, réalisée lors du décollage, à un instant t2 choisi par exemple 3 secondes après la fin du signal d'enfoncement des trains "Weight off wheels". Cette deuxième masse au parking mp2 est déterminée à partir de caractéristiques de l'aéronef à l'instant t2, notamment la poussée nominale des moteurs au décollage dans des conditions de tarage, de température et de pression saisies par le pilote, et la traînée théorique de l'aéronef dans la configuration "train sorti + volets décollage" dans des conditions mesurées de vitesse et d'incidence de l'aéronef. La phase 25 comprend la détermination d'une masse m2 de l'aéronef 1 à l'instant t2 lors du décollage, à laquelle est ajoutée la masse de carburant FU(t2) consommée depuis le décollage. La phase 25 comprend à cette fin la détermination des forces de poussée F et de traînée T exercées sur l'aéronef 1 à cet instant t2. La phase 25 comprend en particulier la mesure de la température totale TAT par le capteur de température 5c et la détermination de la température statique SAT par le calculateur 3 à partir de cette mesure. Elle comprend en outre la mesure de l'incidence a par le capteur d'incidence 5d, la mesure de la pression statique Ps par le capteur de pression statique 5a, de l'accélération longitudinale Jx et selon l'axe Z J2 de l'aéronef 1 par l'accéléromètre 9, et d'une information de vitesse de l'aéronef.. Cette information de vitesse est avantageusement un nombre de Mach indépendant des pressions mesurées et de l'incidence a, noté M**. Ce nombre de Mach M** est déterminé à partir du triangle des vitesses, selon la relation : M (t2)= 11 (7) GS +1/(GPS)k> -171/(t2) SAT(t2) dans laquelle : - GS et Vz(GPS) sont respectivement la vitesse horizontale de l'aéronef 1 par rapport au sol et la vitesse de l'aéronef 1 selon l'axe vertical B, dérivée de l'altitude pression GPS par rapport au temps, déterminées à partir de mesure du capteurs GPS 10 à l'instant t2 ; - ire est un vecteur unitaire ascendant parallèle à l'axe vertical B ; - W(t2 ) est un vecteur vitesse du vent au décollage saisi par le pilote ; et - SAT(t2)est la température statique au décollage. A cet instant t2, les roues de l'aéronef 1 ne sont plus en contact avec le sol, de telle sorte que la force de traînée T est déterminée à partir des valeurs du nombre de Mach M** et d'incidence a déterminées. La force de traînée T est par exemple déterminée à partir de l'équation de propulsion de l'aéronef donnée par : T = 0,7SPsM**2 x(ce,M* * ,conf) (8) où Cx désigne une estimation du coefficient de traînée de l'aéronef, fonction de l'incidence a, du nombre de Mach M** et de la configuration de vol de l'aéronef 1, notée « conf », relative notamment aux volets sortis (ici configuration premier segment). Le coefficient de traînée ax est par exemple estimé à partir d'une table de correspondance stockée dans la mémoire 4 du calculateur 3, fournissant une estimation de la valeur de ax à partir du nombre de Mach M** et de l'incidence a de l'aéronef 1. Cependant, ce coefficient de traînée dépend faiblement du nombre de Mach, et il peut être avantageux, pour simplifier le procédé, de déterminer ax en fonction de l'incidence a et de la configuration seulement, voire de fixer la valeur de ax . La force de poussée F, due aux moteurs, est une fonction connue de la température et de la pression. La masse m2 de l'aéronef 1 au décollage est liée à l'accélération Jx, à la traînée T et à la force de poussée F par la relation : F +T . x x X M2 (9) où Fx et Tx sont les projections de la force de poussée F et de la traînée T respectivement sur l'axe X. La masse m2 est par exemple déterminée à partir de l'accélération, de la traînée et de la force de poussée, au moyen d'un abaque préalablement établi et stocké dans la mémoire 4. L'abaque est établi en déterminant par pesée de manière précise la masse exacte de l'aéronef pour une série de vols d'essais, puis en mesurant pour chaque vol d'essai l'accélération Jx correspondante, avec différentes configurations de l'aéronef. Cette accélération est mesurée avec un accéléromètre. La deuxième masse au parking mp2 est déterminée à partir de la masse au décollage m2 et d'une mesure par les débitmètres de la masse de carburant consommée entre les instants t1 et t2. La deuxième masse au parking mp2 est ainsi égale à : mp2 = m2 (t2 )+ FU(t2) (10) où FU(t2)désigne la masse de carburant consommé entre les instants t1 et t2. L'étape 21 comprend en outre une phase de détermination 26 de la troisième masse au parking mp3 de l'aéronef 1, également réalisée lors du décollage. La détermination 26 de la troisième masse au parking mp3 comprend la détermination d'une masse m3 de l'aéronef 1 à l'instant t2 lors du décollage, à laquelle est ajoutée la masse de carburant FU(t2) consommée depuis le décollage. La masse m3 est calculée selon une équation de sustentation de l'aéronef mettant en correspondance sa masse effective à l'instant t2, le facteur de charge nz, l'incidence mesurée a, le nombre de Mach estimé M**, et la pression statique mesurée Ps. L'équation de sustentation de l'aéronef est donnée de manière générale par : nzm g =0,7SPsM25z(a,M,conf) (11) où m est la masse de l'aéronef, S représente une surface de référence de l'aéronef, nz le facteur de charge de l'aéronef 1 selon l'axe Z, M le nombre de Mach de l'aéronef et Ez un coefficient de portance de l'aéronef 1 selon l'axe Z. La phase 26 comprend ainsi la détermination du facteur de charge nz, de la pression statique Ps, du nombre de Mach et du coefficient de portance de l'aéronef 1 à l'instant t2. Le facteur de charge nz à l'instant t2 est alors déterminé par le calculateur 3 à partir de l'expression : nz(t2)=Jz(t2 ) g où Jz est l'accélération selon l'axe Z déterminée par l'accéléromètre 9 à l'instant t2. (12) Le coefficient de portance az (t2) estimé à l'instant t2 est une projection selon l'axe Z perpendiculaire à l'axe longitudinal de l'aéronef 1 du coefficient de portance C, selon un axe orthogonal au vecteur vitesse de l'aéronef 1, et du coefficient de traînée C, parallèlement à ce vecteur vitesse. Le coefficient de portance az vérifie donc : az = sin a + Cz cos a . Il comprend en général majoritairement un terme de sustentation. Le coefficient de portance az est fonction de l'incidence a, du nombre de Mach M et de la configuration de vol de l'aéronef 1, relative notamment aux volets sortis. A l'instant t2 considéré, il s'agit de la configuration premier segment (train sorti et braquage de volets choisi pour le décollage). Le coefficient de portance az est par exemple estimé à partir d'une table de correspondance stockée dans la mémoire 4 du calculateur 3, fournissant une estimation de la valeur de az à partir du nombre de Mach et de l'incidence a de l'aéronef 1. Cette table de correspondance est préalablement déterminée en effectuant une série de vols d'essais pour un modèle d'aéronef donné. Différentes configurations d'incidence, de facteur de masse, de charge et de Mach, et différentes configurations d'aéronefs sont balayées pour déterminer dans chaque cas le coefficient az . Le coefficient de portance est ainsi estimé au moyen de la table de correspondance à partir du nombre de Mach M**(t2) déterminé à l'instant t2 et d'une valeur d'incidence a(t2) mesurée par le capteur d'incidence à cet instant t2, dans la configuration premier segment. Il est ainsi noté az**(t2) : az (t2)= az(a(t2),M (t2 ),ter segment). La masse m3 de l'aéronef 1 à l'instant t2 est ainsi déduite de l'équation de sustentation selon l'expression : m3 (t2 ) (t2) (13) nz(t2)g La masse au parking mp3 est déduite par ajout de la masse de carburant consommé FU(t2) entre l'instant t1 et l'instant t2, déterminée à partir de mesures des débitmètres. mp3 = m3 (t2 )+ FU(t2) (14) Les trois masses au parking mpi, mp2 et mp3 ainsi déterminées sont généralement différentes. L'étape 21 comprend ainsi une phase 27 de détermination d'une valeur de masse au parking initiale de l'aéronef 1, notée mp,o, à partir des trois masses au parking mpi, mp2 et mp3 préalablement déterminées. Cette masse au parking mpio est avantageusement égale à la médiane des masses mp1, mp2 et mp3. Ainsi, si une de ces masses a une valeur aberrante, cette valeur aberrante est exclue et n'intervient pas dans la valeur de la masse au parking mpio. Lors du vol de l'aéronef 1, des paramètres de vol sont déterminés. Ces paramètres comprennent des caractéristiques de l'air ambiant et du vol de l'aéronef 1, et sont déterminés de manière continue ou périodique à des instant ty du vol, lors d'une étape 28. Cette étape 28 comprend une phase de détermination 29 de caractéristiques de l'air ambiant, notamment de la pression statique Ps et des températures statique SAT et totale TAT. L'étape 28 comprend également une phase de détermination 31 de caractéristiques de vol de l'aéronef 1, et en particulier de son incidence cc par le capteur d'incidence 5d, de l'altitude pression Zp par l'altimètre 6, du nombre de Mach Ma par le machmètre 8 et de l'accélération J par l'accéléromètre 9, à l'instant t,, considéré. Le nombre de Mach Ma constitue une deuxième information de vitesse, en plus de l'équivalent de vitesse EV. L'étape 28 comprend par ailleurs une phase de détermination 33 du coefficient de portance az et du facteur de charge nz de l'aéronef 1 selon l'axe Z à partir des valeurs du nombre de Mach Ma, d'incidence a et d'accélération J déterminées lors de la phase 31. Comme précédemment décrit, le facteur de charge nz est déterminé à partir de l'expression : g nz =Jz (15) où Jz est l'accélération selon l'axe Z déterminée lors de la phase 31. Par ailleurs, le coefficient de portance Ez est estimé à partir de la table de correspondance stockée dans la mémoire 4 du calculateur 3, fournissant une estimation de la valeur de az à partir des valeurs du nombre de Mach Ma et d'incidence a déterminées lors de la phase 31. L'étape 28 comprend en outre une phase de détermination 39 d'une masse estimée tri de l'aéronef 1 lors de son vol. La détermination 39 de la masse estimée m comprend une phase de détermination 41 d'une quatrième masse au parking de l'aéronef 1, lors du vol de l'aéronef 1, à partir de l'équation de sustentation de l'aéronef. Cette quatrième masse est ainsi appelée masse de sustentation msp. Cette masse de sustentation est déterminée en estimant la masse instantanée de l'aéronef 1 à plusieurs instants t distincts lors du vol, avantageusement en continu, en déduisant à partir de chacune des masses instantanées ainsi évaluées une masse au parking de l'aéronef 1, et en effectuant une moyenne pondérée de ces masses au parking. La masse instantanée de l'aéronef 1 à chaque instant t, notée m4(t), est déduite de l'équation de sustentation à partir de la pression statique Ps mesurée lors de la phase 29, du nombre de Mach Ma déterminé lors de la phase 31, et du coefficient de portance az et du facteur de charge nz déterminés lors de la phase 33, selon l'expression : ma(t)= 07SPS Me Z (a, 26 (a Ma , conf) nzg Une masse au parking m4' est déduite par ajout de la masse de carburant consommée FU(t) entre l'instant t1 et l'instant t, déterminée à partir de mesures des débitmètres. Cependant, cette valeur de masse au parking m4' dépend de mesures qui pourraient être erronées, notamment en cas de turbulence, ou de défaillances d'instruments de mesure. On détermine ainsi lors de la phase 41 la masse de sustentation msp comme une moyenne pondérée des masses au parking m4' déterminées entre l'instant t2 de décollage et l'instant tv, cette moyenne étant réalisée seulement sur les instants t lors desquels les mesures sont considérées comme fiables. La masse de sustentation rnsp est ainsi déterminée à un instant t, à partir de l'expression : tv fxn, (t)(m4(t)+ FU(t)) me ) = t2 t (17) fxn, (t) t2 zn.,(t) désigne ici une fonction caractéristique de masse prenant la valeur 1 ou 0 selon des critères géométriques et dynamiques indépendants visant à caractériser le vol rectiligne horizontal de l'aéronef 1. La valeur de la fonction xr,,(t) à chaque instant t est déterminée par le calculateur 3. Notamment, lorsque l'aéronef 1 est dans une phase de vol stabilisée et calme, c'est-à- dire lorsque la vitesse verticale Vz(GPS) issue du capteur GPS 10 est inférieure à un seuil prédéterminé, la valeur absolue de l'inclinaison (issue des centrales à inertie) est inférieure à 5°, le facteur de charge n z est voisin de 1 et une énergie du facteur de charge, définie comme sa variance, est proche de zéro, xm(t) prend la valeur 1. Au contraire, (16) lorsque le vol est agité ou non stabilisé en ligne droite, ou lorsque le calculateur 3 détecte une défaillance d'un capteur ou instrument de mesure, comme décrit ci-après, il fixe la valeurxm (t) = 0 . La masse de sustentation msp n'a donc pas une valeur figée, mais est ajustée lors du vol. La masse estimée in de l'aéronef 1 est alors déterminée lors d'une étape 43 à partir des valeurs de masse au parking mp,0 et msp. La masse estimée Fil est ainsi du type fri = f(moo,m,p). La masse estimée t7) est par exemple une somme pondérée de la masse médiane moo déterminé lors de l'étape 21 et de la masse de sustentation msp. Les coefficients de pondération sont par exemple calculés à partir de la fonction zn.,(t), de telle sorte que le coefficient attribué à la masse de sustentation soit proportionnel à la fraction de temps lors de laquelle cette masse a pu être déterminée. La masse estimée tri est dans cet exemple donnée par : fr1(t,)=[(1- p(t,))mp,o + p(ty)rnsp(tv)J-FU(t,) (18) où la fonction p(t,)est définie par : tv x (t) fm p(t,)=. " t, -t, Ainsi, lorsque le vol est stabilisé depuis longtemps, p(t,) tend vers 1 et la masse estimée in tend vers la masse de sustentation msp dont on a soustrait la masse de carburant consommé. On remarque en outre que tout biais dans la valeur FU(t,), résultant par exemple d'une erreur affectant la débitmétrie de l'aéronef, n'a d'influence que sur le terme - p(tv))moo, qui est d'autant plus faible que p(tv) est proche de 1, c'est-à-dire que la majorité du vol s'effectue en croisière. Ainsi, lorsque le vol est stabilisé depuis longtemps, un biais dans la valeur FU(tv) n'aura quasiment aucune incidence sur la valeur de la masse estimée tri. L'étape 28 comprend par ailleurs une phase de détermination 45 de paramètres de vols auxiliaires, destinés à permettre le calcul de tests des paramètres de vols issus des capteurs ou instruments de mesure, et à remplacer ces paramètres en cas d'identification d'une défaillance d'un ou plusieurs de ces capteurs ou instruments. Ces paramètres auxiliaires sont donc déterminés de manière continue ou périodique à partir (19) de l'instant t2 de décollage de l'aéronef. L'indice n de ces paramètres correspond à la nième valeur déterminée de ces paramètres lorsque la détermination est périodique. La phase 45 comprend notamment la détermination 47 de paramètres de vol à partir de données issues du capteur GPS 10, la détermination 49 d'une première information de vitesse auxiliaire de l'aéronef 1, indépendante des pressions mesurées, la détermination 50 d'une température statique auxiliaire et l'estimation 51 de la vitesse du vent. La connaissance de l'altitude pression de l'aéronef (ou de la pression statique environnant l'aéronef) est d'importance primordiale pour tout procédé de surveillance analytique des paramètres de vol, visant en particulier à détecter une corruption éventuelle des valeurs de pressions statique et totale mesurées. Ainsi, selon l'invention, une information d'altitude issue d'un capteur de position géographique est en permanence disponible lors du vol de l'aéronef 1. Cette information d'altitude est avantageusement issue du capteur GPS 10, l'altitude étant exprimée de manière classique au-dessus du géoïde de référence WGS 84. La phase 47 comprend ainsi la détermination d'une altitude mesurée par le GPS et notée ZGPS, d'un estimateur d'altitude pression, également appelé altitude pression reconstituée et noté Zp**, et d'une pression statique reconstituée notée Ps**, déterminée à partir de l'altitude pression reconstituée Zp** au moyen d'une table d'atmosphère standard. Deux modes sont à distinguer. Selon un premier mode, appelé mode primaire, l'altitude pression reconstituée Zp** est déterminée en retranchant de l'altitude GPS ZGPS un terme correctif d'altitude afin de tenir compte du décalage entre l'atmosphère standard et l'atmosphère réelle de vol. Ce décalage peut être dû à une pression au niveau du sol inférieure, respectivement supérieure, à la pression standard de 1013 bar, générant un décalage des courbes isobares vers le bas, respectivement vers le haut, et à un gradient de température non standard entre le sol et l'aéronef, modifiant l'espacement des isobares. Ce décalage est déterminé à des intervalles DT constants, par exemple égaux à 30 secondes, par comparaison entre l'altitude pression reconstituée Zp** et l'altitude pression Zp issue de l'altimètre 6. Ces écarts constituent une suite de valeurs AZ(/), où j est un entier naturel. Les instants successifs de mesure de ces écarts seront notés t(j). Ainsi : LZ(j) = ZGPS(t(i))- Z p(t(i)) (20) Le terme correctif d'altitude est calculé à partir des écarts AZ(i) comme une suite AZ**(j) donnée par la formule AZ**(j) = AZ(i), où i est un indice inférieur ou égal à j, choisi selon l'algorithme suivant : - Si aucune incohérence des paramètres anémo-clino-barométriques n'a été déclarée, c'est-à-dire si les mesures des capteurs de pression statique 5a, de pression totale 5b et d'incidence 5d n'ont pas été jugés non fiables, et si AZWE [AZ(j -1)- ST;AZ(j -1)+ ST] l'indice i est égal à l'indice j, c'est-à-dire A**Z(j)= AZ(j). A-Z est donc ainsi automatiquement recalé. ST est un seuil de tolérance. La vérification de l'appartenance de AZ(j) à l'intervalle [AZ(j -1)- ST;AZ(j -1)+ ST] revient à analyser analytiquement les sauts de AZ(j) de manière à refuser un recalage aberrant en cas de corruption de l'altitude pression Zp(t(j)). Le seuil de tolérance ST peut par exemple prendre la valeur de 100 ft + 5%IVzI. Le premier terme de 100 pieds est destiné à couvrir un gradient exceptionnel d'isobare en altitude (générateur du vent géostrophique) ainsi que la fluctuation statistique de la valeur d'altitude issue des GPS ; le second terme égal à 5 pourcents de la vitesse verticale Vz(GPS) est destiné à couvrir, en cas de changement d'altitude, un écart de température d'au plus par rapport au standard 30t en valeur absolue. ST peut avantageusement être adapté suite à l'analyse de résultats d'essais en vol. A titre d'exemple, un vent géostrophique de 200 noeuds au pôle correspond à un gradient d'isobare de 40 pieds en 30 secondes de vol. - Si aucune incohérence des paramètres anémo-clino-barométriques n'a été déclarée, et si 32(j)e [AZ(j -1)- ST; AZ(j -1)+ ST], la valeur de l'indice i reste égale à sa valeur précédente, et A**Z n'est pas recalé. De plus, le processus de recalage de la suite d'écarts d'altitude AZ**(/) sur la suite AZ est arrêté jusqu'à une relance éventuelle par le pilote comme cela sera décrit plus bas. - Lorsque au contraire les paramètres anémo-clino-barométriques déterminés à partir des capteurs de pression statique 5a, de pression totale 5b et d'incidence 5d sont jugés incohérents donc non fiables, comme décrit ci-après, le processus de recalage de la suite d'écarts d'altitude AZ**(1) sur la suite AZ est arrêté. Ainsi, la valeur I est figée au dernier indice "pertinent", c'est-à-dire correspondant à un instant passé du vol pour lequel l'incohérence ne s'était pas encore manifestée. Ainsi i En outre, lorsqu'une incohérence entre les paramètres anémo-clino-barométriques issus de capteurs est détectée, l'arrêt du recalage automatique permet de ne plus prendre en compte des valeurs de pression statique et totale douteuses ; en effet, le mauvais fonctionnement des sondes de pression, qu'elles soient statiques ou totales, a un impact direct sur la fiabilité de l'altitude pression élaborée par le calculateur 3. L'altitude pression issue du capteur GPS est affichée par le dispositif d'affichage auxiliaire 14 lorsque les valeurs de pression statique et totale sont jugées non fiables. Le pilote peut alors choisir d'afficher soit l'altitude au-dessus du géoïde de référence ZGPS (il s'agit alors d'un calage « GEO» de l'altitude au-dessus du géoïde de référence), soit l'altitude pression reconstituée Zp** (il s'agit alors d'un calage « Standard» de l'altitude). L'arrêt du recalage automatique de AZ** à cause d'une au moins des deux conditions d'arrêt décrites ci-dessus pourrait nuire, à long terme, à la précision de l'altitude pression reconstituée Zp**. La commande de recalage permet d'éviter une telle imprécision et de relancer le recalage automatique, lorsque le pilote dispose d'informations lui permettant d'accréditer la validité des mesures de pression statique. Cependant, si aucune information disponible à bord ne vient accréditer la validité de ces mesures de pression statique, l'altitude GPS de l'aéronef 1 au-dessus du géoïde et la pression statique environnant l'aéronef 1 sont évaluées selon un mode secondaire, en fonction du QNH local. Le QNH est un code international destiné au calage de l'altimètre de telle sorte que celui-ci indique l'altitude topographique du terrain où est élaboré le QNH lorsque l'aéronef se trouve au sol sur ce terrain. Le QNH local, obtenu ponctuellement par prévision météorologique, ou bien par des moyens radiocom ou datalink, est alors inséré dans ce mode secondaire par le pilote, et le calculateur 3 extrait la pression Pst statique théorique dans une table d'atmosphère standard, correspondant à une altitude pression située à une altitude ZGPS au-dessus de l'isobare du QNH. Dans ce mode secondaire, la pression statique reconsituée Ps** est égale à la pression statique théorique Pst. Il est à noter que cette estimation de la pression statique est d'autant meilleure que l'avion vole près de la surface du sol ou de la mer. Avantageusement, l'entrée dans ce mode secondaire de reconstitution de la pression statique peut se faire par appui sur la même commande que la commande de relance de AZ** en mode primaire. Le duplexage de cette commande peut être réalisé implicitement par l'affichage d'un calage altimétrique par le pilote sur le dispositif d'affichage auxiliaire 14 : si le calage demandé est « standard », l'actionnement de la commande de recalage relancera le processus de recalage automatique selon le mode primaire ; si le calage est de type manuel (QNH affiché), typiquement sous la surface de transition, la commande de recalage activera le mode secondaire de reconstitution de la pression statique. La première vitesse auxiliaire de l'aéronef 1 déterminée lors de la phase 49 est indépendante de la pression statique Ps issue du capteur de pression statique 5a. Elle est avantageusement extrapolée en temps réel à partir de la valeur courante de l'incidence mesurée. Cette vitesse est par exemple un nombre de Mach, appelé Mach fort et noté Ce nombre de Mach est déterminé périodiquement à partir de l'équation de sustentation, sous la forme d'une série par récurrence, selon la relation : Mn = 0.7S P:* a; (tv 2 (22) ryng où : - n est un indice temporel correspondant à un instant tv, - Ps** est la pression statique reconstituée déterminée lors de la phase 47 à l'instant tv, nz est le facteur de charge déterminé lors de la phase 31 à l'instant tv, in est la masse estimée lors de la phase 39 à l'instant tv, et Cz (tv )= az(oe(t.,),Mn*_1,conf) est un coefficient de portance estimé à l'instant tv. Le coefficient de portance estimé a; (tv) est déterminé à partir d'une table de correspondance, fournissant la valeur de a; (t, ) à partir de la valeur de l'incidence a à l'instant tv déterminée lors de la phase 31, et du dernier Mach fort déterminé, c'est-à-dire le Mach fort Mn*_1 d'indice n-1. Ce coefficient de portance estimé Cz (tv) diffère ainsi du coefficient de portance az déterminé lors de la phase 33 en ce qu'il ne fait pas intervenir le nombre de Mach issu du Machmètre 8 à partir de la pression statique mesurée par le capteur de pression statique 5a et de la pression totale mesurée par le capteur de pression totale 5b. Le Mach fort Mn* ne dépend donc pas des pressions mesurées à l'instant tv. La série Mn est initialisée à partir d'un premier terme M; . Ce premier terme est une valeur initiale fiable, provenant de mesures indépendantes des capteurs de pressions. La fréquence de détermination du Mach fort est par exemple comprise entre 1 et 10Hz, de préférence 4Hz. Le Mach fort dépend cependant de l'incidence a. En cas de défaillance du capteur d'incidence 5d détectée, cette vitesse Mn* n'est donc plus calculée. Le Mach fort Mn* ainsi calculé permet de déduire une estimation de température statique SAT; indépendante des mesures de pression provenant des capteurs 5a ou 5b, un estimateur d'écart de température standard AISA,* , estimant l'écart de la température SAT; par rapport à la température statique SATed attendue à l'altitude pression reconstituée Zp** de l'aéronef, et une estimation de la vitesse horizontale du vent. La température statique SAT; est obtenue lors de la phase 50 à partir de la température totale TAT issue du capteur 5c, supposée fiable, et du Mach fort Mn* , en remplaçant dans l'équation (1) le nombre de Mach M par le Mach fort Mn* suivant l'expression : TAT(tv) TAT(tv) SAT; = ± y-1 A ) (1+ 0.2M'*2 ) 2 'vin L'estimateur d'écart de température standard AISA,, estimant l'écart de la températureSAT, par rapport à la température statique SATed attendue à l'altitude pression GPS Zp** de l'aéronef, est déterminé selon l'expression AISAn* = SATstd -SATn* (24) la température SATsta étant par exemple déterminée à partir d'une table d'atmosphère standard. Cet estimateur d'écart de température standard n'est pas utilisé tant que l'incidence mesurée par le capteur d'incidence 5d est jugée fiable et que le Mach fort peut être calculé. Cependant, il permet, lorsque cette incidence n'est plus jugée fiable, d'estimer une température statique en corrigeant une température statique SATstd obtenue à partir de l'altitude pression GPS Zp**, comme décrit ci-après. La vitesse horizontale du vent par rapport au sol est déterminée lors de la phase 51 au moyen d'un estimateur de vent horizontal Wn , évaluant périodiquement la vitesse du vent en force et en direction. Cet estimateur de vent est avantageusement déterminé à partir du triangle des vitesses selon lequel la vitesse Vs de l'aéronef par rapport au sol, dans un plan horizontal, est la somme de sa vitesse propre Vp et de la vitesse du vent, soit : Wn = Vs - Vp (25) (23) La vitesse horizontale Vs de l'aéronef par rapport au sol, c'est-à-dire sa vitesse dans un plan horizontal au sol, est estimée à partir de mesures du capteur GPS 10 par exemple. Ainsi Vs = GS . La vitesse propre Vp de l'aéronef est la composante horizontale de sa vitesse par rapport à l'air, c'est-à-dire la composante horizontale de sa vitesse vraie. Une vitesse vraie de l'aéronef, appelée vitesse vraie forte et notée TAS,* , est déterminée à partir du Mach fort Mh , suivant une relation analogue à la relation (5) : TAS, = A/1,* VyR SAT,* (26) En négligeant la composante verticale de la vitesse du vent, la vitesse de l'aéronef 1 selon l'axe vertical B par rapport au sol est égale à sa vitesse verticale par rapport à l'air. Ainsi, une valeur de la norme de la vitesse propre, appelée vitesse propre forte et notée V; , est déterminée à partir de la vitesse vraie forte à partir de la relation : TASh2 =V;2 (GPS)2 (27) où Vz (GPS) est la vitesse de l'aéronef 1 selon l'axe vertical B par rapport au sol, déterminée à partir de mesures du capteur GPS 10 par exemple. La vitesse propre Vp de l'aéronef est donc estimée par : 117, = uh\iTAS,*2 - V, (GPS)2 (28) où le vecteur uh est un vecteur unitaire horizontal portant la projection horizontale du vecteur vitesse par rapport à l'air, estimé par l'intermédiaire de l'incidence a. A partir de l'hypothèse selon laquelle le dérapage (i.e. l'angle entre le vent relatif et le cap avion) est nul et l'incidence a de l'aéronef donnée par le capteur d'incidence 5d est fiable, l'estimateur de vent W; horizontal est ainsi déterminé à partir de l'expression : = GS NIM '*2 .(yR SAT,* )- Vz (GPS)2 (29) Lors du vol de l'aéronef 1, le calculateur 3 met en oeuvre plusieurs tests permettant de déterminer un état de crédibilité des valeurs de paramètres de vol déterminées au moyen de mesures de capteurs. En particulier, ces tests permettent de déterminer la cohérence des informations de vitesse, de pression et d'incidence, telles que le nombre de Mach Ma, la pression statique Ps et l'incidence oc. Le procédé comprend ainsi une étape 70 de détermination d'un état de crédibilité des informations de pression statique, de pression totale et d'incidence, au moyen d'un premier test. Ce test 70 est un test de surveillance analytique, mis en oeuvre en continu lors du vol, tant qu'un dysfonctionnement n'a pas été détecté Le test 70 est basé sur l'équation de sustentation. Il est réalisé en déterminant si les paramètres de vol mesurés vérifient effectivement cette équation de sustentation. Cette équation fait en effet intervenir la pression statique Ps, le nombre de Mach Ma et l'incidence a, via le coefficient de portance en axes avion C. . Le test 70 comprend ainsi une phase 72 de détermination de la quantité : n,(t, + 0k7i(tv)g T1(ti, = (30) 0,7SPs(tv)Ma2 (tv )Ez (0e(tv ), Ma (tv m'If) dans laquelle les paramètres m, Ps, Ma, Cz sont ceux déterminés lors de l'étape 28 à un instant tv immédiatement antérieur à l'instant tv de la détermination 72. L'indice k indique qu'il s'agit du kième test 70 mis en oeuvre par le calculateur 3. nz est déterminé à un instant (tv + 0), où 0 est un déphasage temporel déterminé une fois pour toutes dans tout ou partie du domaine de vol, afin de maximiser statistiquement une corrélation entre l'incidence a(tv) et le facteur de charge nz (tv + 0) sur l'aéronef considéré, comme décrit par la suite. L'équation de sustentation est vérifiée si T1=1 en l'absence de toute erreur dans les mesures et dans la loi C. . Ainsi, si les valeurs de la pression statique Ps, du nombre de Mach Ma et de l'incidence a issues de mesures des capteurs sont correctes et si la loi az est correcte, la quantité T1 doit être égale à 1. Un léger écart par rapport à cette valeur peut néanmoins être toléré, en raison notamment du bruit des mesures, de l'écart de la masse estimée th par rapport à la masse réelle de l'aéronef, et de l'imprécision potentielle de la valeur du coefficient de portance C. . Le test 70 comprend ainsi une étape 73 lors de laquelle le calculateur 3 vérifie l'appartenance de la valeur T1 à un intervalle restreint autour de la valeur 1, noté ]1- e,;1 + [. ei est un nombre prédéterminé définissant l'écart toléré, par exemple égal à 0.1, correspondant à une erreur de 10%. Si T1 e ]1- ;1 + e, [, le premier test 70 est positif, et les informations de vitesse, de pression et d'incidence sont jugées de crédibilité optimale 74. Dans l'état de crédibilité optimale 74, les paramètres de vol utilisés directement ou indirectement pour le pilotage de l'aéronef sont donc : - l'altitude pression Z,, ; - la température statique SAT ; - le nombre de Mach Ma ; - la vitesse vraie TAS; - la vitesse indiquée IAS. La vitesse vraie de l'aéronef est déterminée à partir du Mach Ma, suivant la relation : TAS = Ms-M - SAT (31) La vitesse indiquée de l'aéronef, notée IAS pour « Indicated Airspeed », est la vitesse directement issue des mesures de pressions, traduites par l'anémomètre 7. Au cours d'une longue phase de vol stabilisée, la fonction p(t) converge vers la valeur 1, de sorte que la masse estimée m(t) converge vers la masse me (t)- FU(t), qui converge elle-même vers la masse m4(t) issue de l'équation de sustentation. Ainsi, la quantité Ti (tvk ) converge naturellement vers la valeur 1 lors d'un vol stabilisé en croisière. Ceci est obtenu même si la loi az comporte un biais ou si la mesure du capteur d'incidence est biaisée. Dans l'état de crédibilité optimale 74, le dispositif d'affichage auxiliaire 12 affiche un ensemble de caractéristiques de vol issues des mesures de capteurs de l'aéronef 1, déterminées notamment à partir des pressions statique et totale et de l'incidence mesurées. Par exemple, le dispositif d'affichage auxiliaire 12 affiche l'altitude pression Z,,, le nombre de Mach Ma et la vitesse indiquée IAS. Le dispositif d'affichage auxiliaire 12 affiche également la pente air ya , r et l'incidence a mesurées. Le dispositif d'affichage auxiliaire 12 affiche en outre une vitesse indiquée de décrochage de l'aéronef 1, notée IASa, sur la base du facteur de charge nz courant. Il s'agit donc de la vitesse IAS minimale qui peut être atteinte par l'aéronef 1 sans décrocher à la valeur du facteur de charge courante. Cette vitesse de décrochage IASa est déterminée par le calculateur 3 en fonction de l'incidence as de décrochage, de la vitesse indiquée IAS et de l'incidence a courantes, selon la relation : IASs =1AS a- ao as - (32) où ao est l'incidence de portance nulle dans la configuration actuelle de l'aéronef. Elle est généralement négative, sa valeur absolue pouvant atteindre quelques degrés au braquage de volets le plus élevé. Cette incidence oto est tabulée suivant la configuration de l'aéronef. Si au contraire Ti « - ;1 + e, [ , T1 correspondant à un instant t vki , le test 70 est négatif, ce qui signifie que l'une au moins des informations de pression statique, de pression totale ou d'incidence est erronée. Le calculateur 3 estime alors que les mesures de vitesse, d'altitude barométrique et d'incidence ne sont plus fiables, et active dans une étape de transition 75 un état de crédibilité faible 76. Tant que cet état de crédibilité faible 76 est activé, les paramètres de vol déterminés par le calculateur 3 à un instant tv et destinés à l'aide au pilotage sont indépendants des valeurs de pression statique, de pression totale et d'incidence issues de capteurs à cet instant tv, même s'ils peuvent dépendre de valeurs antérieures mesurées dans l'état de crédibilité optimale. Lors de l'étape de transition 75, le calculateur 3 fige la valeur de l'écart de température standard AISA,* et du décalage d'altitude pression 32- (j) aux dernières valeurs fiables connues. Ces dernières valeurs fiables connues sont celles déterminées un certain laps de temps avant le test 70 précédent, à un instant noté tf : tf est le dernier instant auquel les paramètres de vol ont été considérés comme fiables. On choisira par exemple tf = tvkl - 30s . Notamment, la mesure de température totale TAT par le capteur de température totale 5c est jugée non fiable, de telle sorte que le recalage de l'écart de température standard AISA,* à partir des relations (23) et (24) ne peut être effectuée. L'écart de température standard AISA,* prend donc la valeur constante A/SAtf** = AISAn* (tf). La température statique est alors estimée à partir l'altitude pression reconstituée Zp** selon la relation : SAT- =SATstd(Zp)+ AISAtf- (33) SAT- est appelée température statique faible. La température statique SAT,* n'est quant à elle plus déterminée. Le décalage d'altitude pression AZ* (j), utilisé pour la détermination 47 de l'altitude pression reconstituée Zp** , prend par ailleurs la valeur constante AZ; = AZ (tf) (la valeur de i est figée, comme décrit ci-dessus). Cependant, comme décrit précédemment, cette valeur 32ff peut être modifiée à tout moment par action manuelle du pilote sur la commande de relance, par exemple en basse altitude avant l'atterrissage si le pilote vérifie que son altitude pneumatique indiquée au calage QNH est cohérente avec son altitude géoïde donnée par un instrument de secours ou sa hauteur radiosonde augmentée de l'altitude du terrain. En outre, la valeur de l'incidence n'est pas fiable, de telle sorte que l'estimateur de vent ne peut plus être déterminé selon la relation (29) ci-dessus. Cet estimateur, alors noté W , prend alors une valeur constante Wt*f* = Wn*(tf). Cependant, la valeur de l'estimateur de vent peut être mise à jour ponctuellement, de manière automatique, lorsque la trajectoire de l'aéronef présente une évolution suffisante, selon un critère décrit ci-après. La mise à jour ou recalage de l'estimateur de vent horizontal W en crédibilité faible est réalisée suivant une méthode similaire à la méthode de détermination de l'estimateur de vent horizontal W; en crédibilité forte, en substituant à l'incidence « mesurée par le capteur d'incidence une incidence reconstituée notée aR, et à la vitesse vraie forte TAS' une vitesse vraie estimée notée TÂS . L'évaluation de l'incidence reconstituée aR est réalisée à partir de l'équation de sustentation de l'aéronef, en restant dans un domaine de vol dans lequel le coefficient de portance de l'aéronef peut être exprimé comme une fonction affine de l'incidence aR, dont les coefficients sont prédéterminés, selon l'expression : d; (oeR m) = /1,(A4)+ aR (34) valable dans un domaine restreint d'incidence [ao(M),ai (M)]. Ainsi, l'incidence reconstituée est donnée par : a az -1144) 2(M) (35) = R ,u(M) = 0,7SP;M-2,u(M) ,u(M) De même, la dérivée de l'incidence reconstituée aR est exprimée en fonction de la dérivée du facteur de charge nz, par : mgc?R, =(36) 0,7SPS,M..2,u(M)nz Les valeurs de 2(M) et de,u(M) sont déterminées à partir des tables donnant la valeur du coefficient de portance en fonction de l'incidence et du nombre de Mach. La vitesse vraie estimée TÂS est alors estimée comme un quotient entre une accélération de l'aéronef et une vitesse angulaire de cet aéronef, sur la base de l'incidence reconstituée a R et de sa dérivée. Cette estimation est valable sous réserve que l'évolution de la trajectoire de l'aéronef soit suffisamment nette, selon le critère décrit ci-après, et en supposant que le dérapage est nul. Le critère utilisé pour estimer si l'évolution de la trajectoire est suffisante repose sur la formule suivante qui exprime, par le biais de la relation fondamentale de la dynamique, que la projection de l'accélération de l'avion dans le référentiel terrestre, projetée dans le plan de l'aéronef, est supérieure à une valeur déterminée, avantageusement 0.4g : gn+g > 0,40 g (37) où Fi est le facteur de charge et g n + g l'accélération de l'aéronef. Ainsi, l'évolution de la trajectoire de l'aéronef est jugée suffisante si : (cos 0 sin çoi)2 + (n z - cos çoi cos 0)2 > 0,16 (38) où Odésigne l'assiette de l'aéronef et çoi son angle de gîte. On remarque notamment qu'il suffit d'avoir (cos Osin vi )2 > 0,16 , donc d'incliner suffisamment l'avion (quelles que soient la poussée, la traînée et la portance) pour obtenir une évolution de trajectoire "suffisante". Typiquement, lorsque l'aéronef 1 est en trajectoire proche du vol rectiligne, le calcul de la vitesse vraie estimée sera déclaré invalide et l'estimateur de vent horizontal W ne sera pas mis à jour. A l'inverse, lorsque l'évolution de la trajectoire est estimée suffisante, la vitesse vraie estimée TÂS est estimée, ce qui permet de recaler l'estimateur de vent horizontal W selon l'expression : W =GS-tThVTÂS2 -1/z(GPS)2 (39) où vh est un vecteur unitaire horizontal portant la projection horizontale du vecteur vitesse air estimé par l'intermédiaire de l'incidence reconstituée aR. Lors de l'étape de transition 75, le calculateur 3 fixe par ailleurs %m (t) = 0 , de sorte que les valeurs de pression statique, de Mach et d'incidence ne sont plus utilisées pour le calcul de la masse estimée. Ainsi, la masse estimée in à un instant tv > tf est donnée par : iii(tv). f )- FU[tf;tv] (40) où FU[tf;ti désigne la masse de carburant consommée entre tfet tv. L'estimation de la vitesse horizontale du vent et de l'écart de température standard permet, en complément de l'altitude pression reconstituée Zp**, de calculer des estimations de vitesse de l'aéronef 1 à partir de sa vitesse GS par rapport au sol. Ces estimations sont par exemple une estimation du nombre de Mach de l'aéronef, appelé Mach faible et noté M-, et une estimation d'un équivalent de vitesse, appelé équivalent de vitesse faible et noté EV**. La vitesse GS de l'aéronef 1 par rapport au sol est par exemple déterminée par le calculateur 3 à partir du capteur GPS 10. Elle peut également être déterminée à partir de l'accélération de l'aéronef telle que mesurée par des centrales à inertie IRS (Inertial Reference System) de l'aéronef 1, éventuellement en hybridant les informations reçues du capteur GPS et des centrales à inertie. La vitesse GS de l'aéronef 1 et l'estimation de la vitesse horizontale du vent W , égale à la valeur figée We ou à une valeur recalée, permettent de déterminer un vecteur vitesse propre de l'aéronef, appelé vitesse propre faible et noté Vp donné par : Vp- = GS-W (41) La vitesse verticale Vz (GPS) de l'aéronef 1 selon l'axe vertical B par rapport à l'air est déterminée à partir de mesures du capteur GPS 10, en négligeant le vent vertical. Une vitesse vraie de l'aéronef, appelée vitesse vraie faible et notée TAS**, est alors déterminée comme une norme du vecteur vitesse de l'aéronef par rapport à l'air, somme de la vitesse propre et de la vitesse verticale, selon la relation : TAS- (t). GS+Vz(GPS)Z-W- (42) dans laquelle : GS est la vitesse horizontale de l'aéronef 1 par rapport au sol à l'instant t ; Vz(GPS) désigne la vitesse de l'aéronef 1 selon l'axe vertical B à l'instant t ; k> est un vecteur unitaire ascendant parallèle à l'axe vertical B ; W- est l'estimateur de vent figé ou recalé. En fonction de l'altitude pression estimée Z*p* et de l'écart de température standard A/SAtf** figé, on estime par ailleurs la température statique environnante SAT =SATed(Zp**)+ 3.1SA; , et la vitesse estimée du son -N/yR - SAT- .** Le Mach faible M** est ensuite déterminé à partir de la vitesse vraie faible et de la vitesse du son estimée, par la relation : M > tf)-- iTAS (t) yR - c'est-à-dire GS+Vz PS)ie - tiVi7 yR SAT- Le Mach faible M** est ainsi indépendant de toute valeur de pression mesurée postérieurement à tf. Le calcul périodique du Mach fort Mn* est interrompu, car ce calcul dépend de l'incidence a. En outre, un équivalent de vitesse faible, noté EV**, est déterminé à partir du Mach faible M**et de l'altitude pression reconstituée Zp- selon la relation : 1 -1 p0EV**2 = 2 '-- *2 (45) 2 L'équivalent de vitesse faible EV** est ainsi obtenu à partir du Mach faible M** et de la pression statique reconstituée P's- par la relation : EV- = 24Ds M-, la pression Po statique reconstituée Ps étant elle-même déterminée à partir de l'altitude pression reconstituée Zp** selon une table de l'atmosphère standard. Dans l'état de crédibilité faible 76, les paramètres de vol utilisés directement ou indirectement pour le pilotage de l'aéronef, à un instant t,, sont donc : - l'altitude-pression GPS Zp** = Z GPS tf ; - la température statique SAT'rt =SATstd(Zp)+ 3.ISAtf- ; - l'estimateur de vent horizontal W figé à la valeur Wt*f* ou recalé ; - le Mach faible M , - la vitesse vraie faible TAS ; - l'équivalent de vitesse faible EV**. M**(t>tf)= > tf (43) (44) Ainsi, en crédibilité faible, les seules mesures en temps réel utilisées sont celle issues du capteur GPS. Pour affiner le test 70 et discriminer les erreurs sur la pression statique Ps ou sur la pression totale PT, dues à un dysfonctionnement des sondes, d'une erreur sur la mesure de l'incidence a, due à un dysfonctionnement du capteur d'incidence 5d, le calculateur 3 effectue un deuxième test 78. Ce test 78 est effectué dans la suite du test 70, pendant un délai d limité après t , par exemple de 5 secondes. Ce deuxième test 78 est destiné à déterminer la crédibilité de l'incidence a telle que mesurée lors de la phase 31. Ce test 78 d'incidence est réalisé en déterminant si l'équation de sustentation est vérifiée, en remplaçant dans cette équation la pression statique Ps mesurée par le capteur de pression statique 5a et le nombre de Mach Ma déterminé par le machmètre 8 à partir de mesures des capteurs de pression statique et totale, par des paramètres de vols auxiliaires correspondants. Ces paramètres de vols auxiliaires correspondants sont la pression statique reconstituée Ps**, déterminée conformément à la phase 45, et le Mach faible M** déterminé lors de l'étape 76. En outre, le coefficient de portance, noté az** , est estimé à partir de la valeur de l'incidence a à l'instant tv déterminée lors de la phase 31, et du Mach faible M** . Le test 78 comprend la détermination 80 du rapport T2 vérifiant : r T2 (tvl + d > tv > tvkl)= nz (t v + 0)riv)g (46) 0,7SP; M** (tv )6;* (tv où az- (tv)=az(a(tv),M- (tv),conf). Cette étape 80 est suivie d'une étape 82 lors de laquelle le calculateur 3 vérifie l'appartenance de T2 à l'intervalle 11-- 62;1+ e2[ , où 82 est un paramètres de tolérance prédéterminé, par exemple égal à el. Le rapport T2 ne dépend pas de la pression statique et du nombre de Mach mesurés. Ainsi, un écart de T2 par rapport à la valeur 1 théorique attendue incrimine spécifiquement la mesure de l'incidence a, car il s'agit de la seule mesure postérieure à tyk1 autre que des mesures issues du capteurs GPS, qui sont quant à elles supposées fiables. Si T2 }1- e2;1+ e2[, le calculateur 3 identifie une défaillance du capteur d'incidence 5d, n'excluant pas une défaillance des capteurs de pression. Le calculateur 3 maintient ainsi l'état de crédibilité faible 76. Le dispositif d'affichage auxiliaire 14 affiche alors des caractéristiques de vol auxiliaires ne faisant pas intervenir les mesures issues des capteurs de pression 5a, 5b et d'incidence 5d, ces mesures étant jugées peu fiables. Ainsi, chacune des caractéristiques de vol faisant intervenir une de ces mesures dans l'état de crédibilité optimale est remplacée par une caractéristique homologue, i.e. représentative de la même caractéristique de vol, sur le dispositif d'affichage auxiliaire 14. Par exemple, l'altitude pression Zp, le nombre de Mach Ma et la vitesse indiquée IAS sont remplacés sur le dispositif d'affichage auxiliaire 12 par des caractéristiques homologues, qui sont respectivement l'altitude pression reconstituée Zp** ou l'altitude ZGPS, selon le choix du pilote, le nombre de Mach faible M**, et l'équivalent de vitesse faible EV**. En outre, la pente air yak et l'incidence a ne sont plus affichées, et disparaissent donc du dispositif d'affichage auxiliaire 14. Par ailleurs, la vitesse indiquée de décrochage IASS est remplacée par un équivalent de vitesse de décrochage faible EV;* sur la base du facteur de charge n, courant. L'équivalent de vitesse de décrochage faible EVs** est déterminé par le calculateur 3 en fonction du facteur de charge n, courant et d'un équivalent de vitesse de décrochage EVso sous une accélération de 1g, selon la relation : EV;* = EVsor.i); (47) L'équivalent de vitesse de décrochage EVso est obtenu à partir d'une table de décrochage en fonction de la masse estimée if; de l'aéronef 1. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité faible. Ce message indique ainsi au pilote que les mesures des capteurs de pression et d'incidence ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Si au contraire T2 E 11- e2;1+ e2{, et sous réserve qu'un test complémentaire évaluant une cohérence dynamique de l'incidence mesurée avec le facteur de charge, détaillé ci-après, n'ait pas fourni un résultat négatif pendant une durée prédéterminée, par exemple 5 minutes, le calculateur 3 exclut une défaillance du capteur d'incidence 5d, ce qui implique une défaillance d'au moins un des capteurs de pression statique et totale. Le calculateur 3 lève ainsi l'état de crédibilité faible, et active dans une étape 83 un état de crédibilité forte 84. Lors de la transition 83, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité forte. Ce message indique que les mesures des capteurs de pression ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Dans cet état de crédibilité forte 84, les paramètres de vol déterminés à un instant tv sont indépendants des valeurs de pression statique et totale issues de capteurs à cet instant t', mais peuvent dépendre de l'incidence issue du capteur d'incidence 5d. Lors de la transition 83 d'un état de crédibilité faible 76 à un état de crédibilité forte 84, le calculateur 3 réinitialise le calcul de la suite des Mach forts, interrompu dans l'état de crédibilité faible, à partir d'un terme initial M; correspondant à la dernière valeur du nombre de Mach faible avant la transition 83. Le terme initial M; est ainsi remis à jour à chaque fois que l'état de crédibilité forte 84 est activé. La température statique SAT; est à nouveau évaluée, suivant l'expression : SAT- = , TAT(ty) - + 0.0.2M;,2)(48) Par ailleurs, le calculateur 3 reprend la détermination de l'estimateur d'écart de température standard A/SA,* , interrompue dans l'état de crédibilité faible en raison de l'indisponibilité du nombre de Mach fort Mn . Cet estimateur est déterminé conformément à l'étape 50 ci-dessus. , Le calculateur 3 reprend également le calcul de l'estimateur de vent 1/V' , évaluant périodiquement la vitesse du vent par rapport au sol, conformément à la phase 51 ci-dessus. Ce calcul fait en effet appel au nombre de Mach fort et à la température statique forte SAT; , à nouveau disponibles. La vitesse vraie forte TAS,* est déterminée à partir du nombre de Mach fort Mn* et de la température statique SAT,* selon l'expression (23) ci-dessus. En outre, un équivalent de vitesse fort, noté EV; , est déterminé à partir du Mach fort Mn* et de l'altitude pression reconstituée Zp** selon la relation : *2 1 -2 poEV, = -2 yPs M,*2 (49) Ainsi, les paramètres de vol utilisés directement ou indirectement pour le pilotage de l'aéronef 1 en crédibilité forte, sont : - l'altitude pression GPS Zp** - la température statique SAT; ; - l'incidence mesurée a ; ZGPS 32tf ; - le Mach fort Mn* ; - la vitesse vraie forte TAS,* ; - l'équivalent de vitesse fort EV,* . Dans l'état de crédibilité forte 84, le dispositif d'affichage auxiliaire 14 affiche des caractéristiques de vol auxiliaires, ne faisant pas intervenir les mesures issues des capteurs de pression. Ainsi, chacune des caractéristiques de vol faisant intervenir une mesure d'un capteur de pression dans l'état de crédibilité optimale est remplacée par une caractéristique homologue sur le dispositif d'affichage auxiliaire 14, indépendante de toute mesure de capteurs de pression réalisée dans l'état de crédibilité forte. Par exemple, l'altitude pression Zp est remplacée sur le dispositif d'affichage auxiliaire 14 par une caractéristique homologue, qui est l'altitude pression reconstituée Zp** ou l'altitude ZGPS, selon le choix du pilote. Le nombre de Mach Ma est remplacé par le nombre de Mach fort M,7* , et la vitesse indiquée IAS est remplacée par l'équivalent de vitesse fort EV,* . La pente air yak et l'incidence a mesurées restent affichées. Le dispositif d'affichage auxiliaire 12 affiche en outre un équivalent de vitesse de décrochage fort noté EV; , déterminé par le calculateur 3 de manière analogue à la vitesse indiqué de décrochage IASs, en fonction de l'équivalent de vitesse fort EV,* , de l'incidence as de décrochage et de l'incidence a courantes, selon la relation : \I EV; = EV, a - ao (50) as - ao La détermination de ces caractéristiques de vol et de l'air ambiant permet ainsi de fournir des informations au pilote de l'aéronef tout en s'affranchissant des pressions statique et totale potentiellement erronées mesurées par les capteurs. Comme précédemment décrit, dans l'état de crédibilité forte, la température totale mesurée par le capteur 5c et l'incidence mesurée par le capteur 5d sont jugées fiables, ce qui permet, grâce au calcul de l'altitude pression reconstituée Zp** , de déterminer un nombre de Mach (Mach fort) et une vitesse vraie (vitesse vraie forte) de l'aéronef. Ces informations de vitesse permettent de déterminer une image de l'atmosphère, et en particulier un estimateur de vent horizontal et un estimateur d'écart de température standard. A l'inverse, dans l'état de crédibilité faible, la température totale mesurée par le capteur 5c et l'incidence mesurée par le capteur 5d sont jugées non fiables. L'altitude pression reconstituée Z*p* , l'estimateur de vent horizontal et l'estimateur d'écart de température standard sont alors figés aux dernières valeurs fiables connues, l'altitude pression reconstituée Zp** et l'estimateur de vent horizontal pouvant cependant être recalés. En outre, la valeur figée de l'écart de température standard AISA,* (telle que déterminés dans l'état de crédibilité forte) permet d'estimer une température statique à partir de l'altitude pression reconstituée Z. . La connaissance de cette image simplifiée de l'atmosphère permet alors de déterminer un nombre de Mach (Mach faible) et un équivalent de vitesse (équivalent de vitesse faible) de l'aéronef. Quel que soit l'état de crédibilité forte ou faible en vigueur dans le système, le calculateur 3 met en oeuvre des tests complémentaires pour évaluer, de manière périodique, la crédibilité des informations fournies par les capteurs de pression et/ou d'incidence, détecter une éventuelle évolution de cette crédibilité et confirmer ou modifier un état de crédibilité faible ou fort. Dans l'état de crédibilité faible 76, les mesures issues des capteurs de pression 5a et 5b et d'incidence 5d sont jugées peu fiables. Pour réévaluer la crédibilité de ces mesures, le calculateur 3 évalue de manière répétée la fiabilité de la mesure d'incidence cc. Si cette mesure est jugée non fiable, le calculateur 3 maintient l'état de crédibilité faible. Si au contraire cette mesure est jugée fiable, le calculateur 3 lève l'état de crédibilité faible, et active un état de crédibilité forte 84, dans lequel seules les mesures de pression sont jugées peu fiables. Dans l'état de crédibilité forte 84, les mesures issues des capteurs de pression 5a, 5b sont jugées peu fiables, tandis que les mesures du capteur d'incidence 5d sont estimées fiables. Si les mesures d'incidence sont jugées fiables lors de l'évaluation par les tests complémentaires, le calculateur 3 maintient l'état de crédibilité forte 84. Si au contraire les mesures d'incidence sont jugées peu fiables lors de l'évaluation par les tests complémentaires, le calculateur 3 active l'état de crédibilité faible 76. Les étapes mises en oeuvre par le calculateur 3 pour évaluer la fiabilité des mesures du capteur d'incidence 5d dans l'état de crédibilité forte 84 ou faible 76 sont illustrées sur la Figure 4. Cette réévaluation repose sur deux tests de crédibilité de l'incidence. Un premier test 90 mesure la cohérence dynamique de l'incidence mesurée oc avec le facteur de charge nz. Une modification de l'incidence a entraîne une variation du facteur de charge nz, avec un léger retard noté 0. Si l'on considère un échantillon de temps, lors duquel l'incidence varie de da, suffisamment court pour que les valeurs m et PSMa2 puissent être considérés comme constantes, la dérivation de l'équation de sustentation sur cet échantillon de temps donne : f 0,7SPsMa2 paz mg Da ) ) Si l'on considère de plus que la grandeur aa aC est constante sur l'échantillon de temps, on obtient dnz = K.da , où K est une constante. Par conséquent, si la valeur de l'incidence mesurée par le capteur d'incidence 5d est juste, cette incidence doit être une fonction affine du facteur de charge nz sur un échantillon de temps court. Le test 90 repose ainsi sur l'évaluation de la covariance entre l'incidence a et le facteur de charge nz, mesurant la corrélation entre ces deux paramètres. Le test 90 comprend une étape 92 de détermination de l'applicabilité du test 90, une étape 94 de calcul d'un coefficient de corrélation entre le facteur de charge et l'incidence, et une étape 96 de comparaison de ce coefficient de corrélation à une valeur seuil prédéterminée. Le test 90 ne peut être effectué que si l'incidence et le facteur de charge varient suffisamment sur l'échantillon de temps considéré. La détermination 92 de l'applicabilité du test 90 comprend l'évaluation d'une variance du facteur de charge nz, mesurant la dispersion du facteur de charge n, autour de sa moyenne sur un intervalle de temps choisi, et appelé énergie empirique brute. L'intervalle de temps est par exemple égal à 2.5 secondes. L'énergie empirique brute en, du facteur de charge est alors donnée par : 1 100 f e2nz 100 = - nz to 40 0 - nz (52) où nz , moyenne du facteur de charge n, sur l'intervalle [to + 0; to + 2.5 + 0], est donné par : 1 100 nz Enz(to+-k 40 +0) (53) 100 k=1 0 est le déphasage entre le facteur de charge et l'incidence et to un instant de mesure choisi. 0 est par exemple égal à 0.01 seconde. Ce déphasage 0 est dépendant dnz = da (51) d'un type d'avion. Il est déterminé une fois pour toutes dans tout ou partie du domaine de vol afin de maximiser en moyenne, dans la partie du domaine de vol considérée, la corrélation entre a(t) et n,(t+0). Lors de l'étape 92, le calculateur 3 compare ensuite la valeur de l'énergie empirique brute en, du facteur de charge à une valeur seuil minimale en, min prédéterminée. La valeur seuil en, min est par exemple comprise entre 0,01 et 0,02, avantageusement égale à 0,015. Si enz Puis, dans l'étape 96, le calculateur 3 compare le coefficient de corrélation déterminé lors de l'étape 94 à un seuil de corrélation prédéterminé Coffrait, . Si Corro(nz,a) Ce test 90 fournit ainsi une information sûre de la vivacité de l'incidence et de sa cohérence temporelle avec le facteur de charge. Cependant, lorsque ce test est positif, un décalage permanent entre la valeur d'incidence mesurée et la valeur réelle n'est pas exclu. La réévaluation de la fiabilité du capteur d'incidence 5d comprend ainsi un deuxième test 100 d'incidence, vérifiant la cohérence entre l'incidence a mesurée et l'assiette Ode l'aéronef, en vol rectiligne horizontal. Ce test 100 comprend une étape 102 de détermination de la pertinence du test 100 et une étape 106 de comparaison de l'incidence a mesurée et de l'assiette O de l'aéronef lorsque le test est pertinent. Le test 100 peut être effectué lorsque l'aéronef vole en palier rectiligne, c'est-à-dire lorsque la pente air yair de l'aéronef et son inclinaison go sont voisines de zéro. L'étape 102 comprend ainsi la comparaison de la pente air vair de l'aéronef et de son inclinaison çpà des valeurs seuil prédéterminées. La pente air yair de l'aéronef est évaluée à partir de mesures issues du capteur GPS. Le calculateur 3 détermine à cette fin une vitesse verticale ±p** de l'aéronef, dérivée de l'altitude pression reconstituée, et un nombre de Mach de l'aéronef, qui peut être soit le Mach faible M** dans l'état de crédibilité faible 76, soit le Mach fort Mn dans l'état de 2; zPI crédibilité forte 84. Le calculateur 3 détermine alors le rapport Pa = ou Pa =, M* proportionnel au sinus de la pente air yaw, et compare ce rapport à une valeur seuil Pmax- La tolérance est par exemple égale à 0,4°. Par ailleurs, la centrale à inertie fournit au calculateur 3 la valeur de l'inclinaison go de l'aéronef. Le calculateur 3 compare alors cette inclinaison à une valeur seuil yomax, par exemple égale à 5°. Si ço > ç,,,ax et/ou Pa > P'x, le test 100 n'est pas poursuivi. Le calculateur 3 réitère alors l'étape 102 jusqu'à ce que le test 100 puisse être effectué. Si çomax et P Pmax, le test 100 peut être poursuivi et le calculateur 3 passe à l'étape 106. Lorsque q' Lors de l'étape 106, le calculateur 3 évalue la fiabilité du capteur d'incidence 5d en comparant l'assiette O, déterminée au moyen de la centrale à inertie, à l'incidence a mesurée par le capteur d'incidence 5d. Par exemple, le calculateur 3 détermine la différence (19- a)et vérifie que cette différence est proche de la valeur 0, c'est-à-dire est qu'elle est inférieure à un seuil e3 en valeur absolue, où e3 est un nombre prédéterminé définissant l'écart toléré. On prendra par exemple e3 =1° . Si I0 -- > e3, le test 100 est négatif. Le calculateur 3 réitère alors l'étape 80 tant qu'un résultat positif n'a pas été obtenu. Si IO -al5e3, le test 100 est positif. Ce test 100 permet de valider la valeur absolue de l'incidence en vol rectiligne horizontal. Il permet ainsi d'éliminer des erreurs constantes sur la valeur de l'incidence. Les tests 90 et 100 sont répétés en permanence par le calculateur 3. Les deux tests 90 et 100 ne sont pas toujours en mesure de donner un résultat de manière concomitante. En effet, le test 90 de corrélation entre l'incidence et le facteur de charge est applicable lorsque le facteur de charge varie, tandis que le test 100 de comparaison entre l'incidence et l'assiette est applicable lorsque l'aéronef est en palier rectiligne. Les résultats positifs à ces deux tests ne sont donc pas nécessairement obtenus ensemble. Ainsi, lorsque le calculateur 3 obtient un résultat positif à l'un des tests 90, 100, il active 107 un état intermédiaire et mémorise ce résultat pendant une durée prédéterminé At, en attente d'un résultat positif ou négatif à l'autre de ces tests 100, 9. La durée At est par exemple égale à 5 min (300 s). Si un résultat positif est obtenu pour le deuxième de ces tests 100, 90 pendant cette durée At, le calculateur 3 considère que les mesures du capteur d'incidence 5d sont fiables. A noter que si ces résultats sont obtenus durant l'intervalle kvki,tvki + di, c'est-à- dire pendant le délai d après un test 70 négatif, ils ne sont pas pris en compte. Si l'état de crédibilité en cours est l'état de crédibilité forte 84, le calculateur maintient cet état. Si l'état de crédibilité en cours est l'état de crédibilité faible 76, le calculateur 3 lève cet état de crédibilité faible 76 et active l'état de crédibilité forte 84, précédemment décrit. Notamment, la transition entre ces deux états est similaire à la transition 83 décrite précédemment. Le test permanent 70 est alors réactivé dès l'entrée en crédibilité forte. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche des caractéristiques de vol auxiliaires adaptées à cet état de crédibilité forte 84. Notamment, l'altitude pression GPS Zp** reste affichée, mais le nombre de Mach faible M** est remplacé par le nombre de Mach fort Mn* , et l'équivalent de vitesse faible EV** est remplacé par l'équivalent de vitesse fort EV,.,* . En outre, la pente air yak et l'incidence a sont à nouveau affichées par le dispositif d'affichage auxiliaire. Enfin, l'équivalent de vitesse de décrochage fort EV; est à nouveau déterminé à partir de mesures d'incidence issues du capteur d'incidence. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité forte. Ce message indique au pilote que les mesures des capteurs de pression 5a, 5b ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Si aucun résultat positif n'est obtenu pour le deuxième de ces tests 100, 90 pendant la durée At, le calculateur 3 considère que les mesures du capteur d'incidence 5d ne sont pas fiables. Si l'état de crédibilité en cours est l'état de crédibilité faible 76, le calculateur maintient cet état de crédibilité faible 76 et réitère les tests 90 et 100. Si l'état de crédibilité en cours est l'état de crédibilité forte 84, le calculateur 3 active l'état de crédibilité faible 76. Notamment, lors de la transition 75 vers l'état de crédibilité faible 76, les valeurs de l'estimateur de vent Wn* , de l'écart de température standard BISA' et du décalage d'altitude pression AZ (j) sont figées aux dernières valeurs fiables connues, qui sont par exemple celles déterminées à un instant tf antérieur au test 90 ou 100 négatif. L'altitude pression reconstituée Zp** et l'estimateur de vent horizontal peuvent cependant être recalés comme décrit ci-dessus. Le dispositif d'affichage auxiliaire 14 affiche par ailleurs des caractéristiques de vol auxiliaires adaptées à cet état de crédibilité faible 76. Notamment, l'altitude pression GPS 4.* reste affichée, mais le nombre de Mach fort Mn* est remplacé par le nombre de Mach faible M**, et l'équivalent de vitesse fort EVn* est remplacée par l'équivalent de vitesse faible EV**. En outre, la pente air ya'- et l'incidence a ne sont plus affichées, et disparaissent du dispositif d'affichage auxiliaire 14. Enfin, l'équivalent de vitesse de décrochage faible EV;* est obtenu à partir d'une table de décrochage en fonction du facteur de charge nz courant et d'un équivalent de vitesse de décrochage EVso sous une accélération de 1g. L'équivalent de vitesse de décrochage EVso est obtenu à partir d'une table de décrochage en fonction de la masse estimée in de l'aéronef 1. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité faible. Ce message indique au pilote que les mesures des capteurs de pression et d'incidence ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. La combinaison des tests 90 et 100, mesurant respectivement la corrélation entre l'incidence et le facteur de charge et la cohérence entre l'incidence et l'assiette en vol horizontal, permet ainsi d'évaluer la fiabilité des mesures d'incidence par le capteur d'incidence 5d. Comme précédemment décrit, l'état de crédibilité faible 76 est activé lorsque le premier test 70 et le deuxième test 78 s'avèrent successivement négatifs. L'état de crédibilité faible 76 est également activé de manière systématique au décollage, et pendant chaque changement de configuration de l'aéronef 1. Après le décollage, l'état de crédibilité faible 76 est levé et l'état de crédibilité forte 84 activé dès que la configuration « deuxième segment » est atteinte, sous réserve que le test 90 de corrélation entre l'incidence et le facteur de charge ne soit pas négatif. De même, dès que le changement de configuration de l'aéronef 1 est achevé, l'état de crédibilité faible 76 est levé et l'état de crédibilité forte 84 établi, sous réserve que le test 90 de corrélation entre l'incidence et le facteur de charge ne soit pas négatif. Par ailleurs, dans l'état de crédibilité forte 84, le calculateur 3 réévalue la fiabilité des capteurs de pression 5a, 5b. Ces tests sont mis en oeuvre de manière permanente dans l'état de crédibilité forte 84. Pour réévaluer la fiabilité des capteurs de pression dans l'état de crédibilité forte, le calculateur 3 met en oeuvre deux tests 110, 112, comme représenté sur la Figure 5. Dans un premier test 110 des capteurs de pression, le calculateur 3 compare le nombre de Mach Ma obtenu par le machmètre 8 à un instant donné à partir des valeurs de pression statique et totale au nombre de Mach fort Mn obtenu au même instant, qui ne dépend pas de ces pressions. En effet, si les valeurs de pression mesurées sont correctes, ces deux vitesses devraient être sensiblement égales. Ainsi, lors du test 110, le calculateur 3 détermine le nombre de Mach Ma conformément à la phase 31 ci-dessus, et le nombre de Mach fort Mn , conformément à l'étape 49 ci-dessus. Le calculateur 3 détermine alors le rapportM et vérifie que ce n rapport est proche de la valeur 1, c'est-à-dire est compris dans un intervalle /1- e4 ;1 + e4 [ , où 64 « 1 , est un nombre prédéterminé définissant l'écart toléré. . M Si - e4;1+ 64 [, le test 110 est négatif. Les valeurs des pressions mesurées Mn sont jugées non fiables, et l'état de crédibilité forte est maintenu. Si M- , -A/1 e e4;1 + 64 [, le test est positif. n Le deuxième test 112 des capteurs de pression est complémentaire au premier test 110. Lors de ce test 112, le calculateur 3 compare l'altitude pression Zp obtenue à un instant donné à partir de la mesure de la pression statique à l'altitude pression GPS Zp** obtenue au même instant, qui ne dépend pas de la mesure de la pression statique. Si la valeur de pression statique mesurée est correcte, ces deux altitudes pression devraient être sensiblement égales. Ainsi, lors du test 112, l'altimètre 6 fournit au calculateur 3 l'altitude pression Zp conformément à la phase 31 ci-dessus, et le calculateur 3 détermine l'altitude pression GPS Zp**, conformément à l'étape 47 ci-dessus. Le calculateur 3 détermine alors la différence z; -Zp, et vérifie que ce rapport est proche de zéro, c'est-à-dire est compris dans un intervalle ]- e5;+ e5 [ , où e5 est un nombre prédéterminé définissant l'écart toléré. Si z; - Zp 65;+ 65 [ , le test 112 est négatif. La valeur de la pression statique mesurée est jugée non fiable, et l'état de crédibilité forte est maintenu. Si Zp -Zp E 65 ; + 65 [ , le test 112 est positif. Dans une étape 113, le calculateur détermine si deux tests 110 et 112 sensiblement concomitants, par exemple effectués dans un intervalle de temps de l'ordre de la seconde, sont positifs. Si c'est le cas, les mesures des capteurs de pression sont jugées fiables. Le calculateur 3 lève donc l'état de crédibilité forte et active l'état de crédibilité optimale 74, décrit ci-dessus. Dans le cas contraire, l'état de crédibilité forte 84 est maintenu. Ainsi, lors d'un changement de configuration ou au décollage, une fois l'état de crédibilité faible levé et l'état de crédibilité forte établi, l'état de crédibilité optimale peut être atteint si les tests 110 et 112 s'avèrent positifs. Ces deux tests 110 et 112 permettent ainsi de réévaluer la crédibilité des mesures de pression lors du vol, lorsque ces mesures ont été jugées peu crédibles, et de fournir au pilote des informations issues de ces mesures lorsque les capteurs de pression fonctionnent à nouveau. L'utilisation de deux tests complémentaires garantit par ailleurs une meilleure sécurité. On a représenté sur la Figure 6 un dispositif d'affichage auxiliaire 14 selon un mode de réalisation de l'invention. Comme précédemment décrit, le dispositif d'affichage auxiliaire 14 est relié au système 2 et en particulier au calculateur 3 de l'aéronef 1, et reçoit des informations et des ordres de commande de celui-ci. Le dispositif 14 d'affichage auxiliaire comprend des moyens 15, propres à afficher des évaluations des valeurs de vitesse, d'altitude, d'incidence et de pente air, choisies en fonction d'estimations de la fiabilité des capteurs d'incidence et de pression et transmises par le calculateur 3. Ces moyens 15 forment un dispositif de visualisation. Le dispositif d'affichage auxiliaire 14 affiche, dans l'état de crédibilité optimale, des caractéristiques de vol dites « principales » de l'aéronef, telles qu'une vitesse, une altitude et un nombre de Mach, obtenues à partir de mesures des capteurs de pression statique 5a, de pression totale 5b et d'incidence 5d. Lorsque les mesures d'au moins un de ces capteurs sont jugées non fiables, chaque caractéristique de vol « principale » obtenue à partir de mesures de ce capteur est remplacée sur le dispositif d'affichage par une caractéristique de vol « auxiliaire », homologue à la caractéristique principale remplacée, i.e. représentative de la même caractéristique de vol, et indépendante de toute mesure du capteur jugé non fiable effectuée lorsque les mesures de ce capteur sont jugées non fiables. Ainsi, la vitesse indiquée IAS est remplacée par une vitesse (l'équivalent de vitesse forte EV,* ou faible EV- ), l'altitude pression Zp est remplacée par une altitude (l'altitude pression GPS Zp** ou l'altitude au-dessus du géoïde de référence ZGPS), et le nombre de Mach Ma est remplacé par un nombre de Mach (nombre de Mach fort Mn ou faible M**). Les Figure 6, 7 et 8 représentent ainsi des informations projetées sur le dispositif 15, affichées sous la forme de symboles, lors de la mise en oeuvre du procédé précédemment décrit, lorsque l'aéronef 1 est dans l'état de crédibilité optimale 74, dans l'état de crédibilité forte 84 et dans l'état de crédibilité faible 76 respectivement. Ces symboles comportent un indicateur 122 de commande de vol, affichant un symbole 124 de maquette de l'aéronef, occupant une position constante, qui matérialise une projection à l'infini de l'axe longitudinal X de l'aéronef 1, et une ligne d'horizon artificiel 126, au centre d'une échelle graduée de pente 128. La position de la ligne 126 par rapport au symbole 124 représente l'assiette Ode l'aéronef 1, cette assiette étant indiquée sur l'échelle graduée de pente 128 en regard de la ligne d'horizon artificiel 126. L'indicateur 122 comprend également un symbole 130 de vecteur vitesse, indiquant la direction du vecteur vitesse de l'aéronef par rapport à l'air. Ce symbole 130 est disposé le long de l'échelle graduée de pente 128. L'écart sur l'échelle graduée de pente 128 entre le symbole 124 de maquette et le symbole 130 de vecteur vitesse est égal à l'incidence a de l'aéronef 1. Avantageusement, les symboles de maquette 124 et de vitesse 130 ne sont affichés que dans les états de crédibilité optimale 74 ou forte 84. La pente air et l'incidence de l'aéronef 1 ne sont ainsi pas affichées dans l'état de crédibilité faible 76. Le dispositif 15 affiche par ailleurs un indicateur de vitesse 132. Cet indicateur 132 comprend une échelle de vitesse 134 graduée, représentée sous la forme d'un segment s'étendant entre deux points fixes, et un symbole de vitesse 136, disposé en regard de l'échelle de vitesse 134. Le symbole de vitesse 136 indique une vitesse de l'aéronef 1. Le symbole 136 a par exemple une forme de chevron. Un symbole 138, indiquant sous forme numérique la valeur de la vitesse, est ajouté en regard du symbole de vitesse 136. La vitesse indiquée dépend de l'état de crédibilité établi par le calculateur 3. La vitesse affichée au moyen du symbole de vitesse 136 est par exemple la vitesse indiquée IAS dans l'état de crédibilité optimale 74, l'équivalent de vitesse forte EV,* dans l'état de crédibilité forte 84 et l'équivalent de vitesse faible EV- dans l'état de crédibilité faible 76. Un symbole 141, disposé en bas de l'échelle de vitesse 134 graduée, indique un nombre de Mach de l'aéronef. Le type de nombre de Mach affiché dépend de l'état de crédibilité courant. Dans l'état de crédibilité optimale 74, il s'agit du nombre de Mach Ma issu du machmètre 8, dans l'état de crédibilité forte 84, le nombre de Mach fort Mn* est affiché, et dans l'état de crédibilité faible 76, le nombre de Mach faible M** est affiché. Les symboles 138 et 141 sont munis d'une identification représentative du type de vitesse et de nombre de Mach affiché, avantageusement un code couleur, de manière à informer le pilote du type de vitesse dont il dispose. Le code couleur choisi indique de manière intuitive quels types de vitesse et de nombre de Mach sont affichés. Par exemple, les symboles 138 et 141 sont verts encadrés de vert dans l'état de crédibilité optimale 74, i.e. quand la vitesse indiquée IAS et le nombre de Mach Ma sont affichés, verts encadrés de jaune dans l'état de crédibilité forte 84, i.e. quand l'équivalent de vitesse forte EV,* et le nombre de Mach fort M* sont affichés, et jaunes encadrés de jaune dans l'état de crédibilité faible 76, i.e. quand l'équivalent de vitesse faible EV- et le nombre de Mach faible M** sont affichés. Comme représenté sur la Figure 6, dans l'état de crédibilité optimale, le dispositif 15 affiche en outre, superposé à l'indicateur de vitesse 132, un ruban 139a indiquant la vitesse indiquée IASS de décrochage de l'aéronef 1. Le ruban 139a est par exemple de couleur rouge. Dans l'état de crédibilité forte 84, un ruban 139a identique indique l'équivalent de vitesse de décrochage fort EV; . Dans l'état de crédibilité faible 76, le ruban 139a est remplacé par un ruban rouge discontinu 139b, comme représenté sur la Figure 8. Le ruban 139b indique l'équivalent de vitesse de décrochage faible EVS* décrit ci-dessus. Le dispositif 15 affiche également un indicateur d'altitude 140. Cet indicateur 140 comprend une échelle d'altitude 142 graduée, représentée sous la forme d'un segment s'étendant entre deux points fixes, et un symbole d'altitude 144, disposé en regard de l'échelle d'altitude 142. Le symbole d'altitude 144 indique une altitude de l'aéronef 1. Le symbole 144 a par exemple une forme de chevron. Un symbole 146, indiquant sous forme numérique la valeur de l'altitude, est ajouté en regard du symbole d'altitude 144. L'altitude indiquée dépend de l'état de crédibilité établi par le calculateur 3. Dans l'état de crédibilité optimale 74, l'altitude affichée au moyen du symbole d'altitude 144 est par exemple l'altitude pression Zp ou, au choix du pilote, une altitude corrigée du calage, notée Zc. Le calage est par exemple un calage QNH ou QFE. Le QFE est un code international destiné au calage de l'altimètre par rapport à un terrain donné de telle sorte que celui-ci indique une altitude nulle lorsque l'aéronef se trouve au sol sur ce terrain. Dans l'état de crédibilité forte 84 ou faible 76, l'altitude affichée au moyen du symbole d'altitude 144 est l'altitude pression GPS Zp** ou, au choix du pilote, l'altitude au-dessus du géoïde de référence ZGPS, donnée par le capteur GPS. Ces choix entre Zp ou Zc et Zp** ou ZGps peuvent être effectués au moyen de l'interface de saisie 17, par exemple par actionnement d'un bouton dédié. Le dispositif 15 indique alors, dans un emplacement réservé 150, sur quelle référence l'altitude est calée. Par exemple, dans l'état de crédibilité optimale 74, si l'altitude affichée est l'altitude pression Zp, le calage « STD » est affiché dans l'emplacement 150, comme représenté sur la Figure 6. Si l'altitude affichée est une altitude corrigée du calage Zc, « QNH » ou « QFE » est affiché dans l'emplacement 150, selon le calage choisi. Dans l'état de crédibilité forte 84 ou faible 76, si l'altitude pression GPS Zp** est affichée, le calage « STD » est affiché dans l'emplacement 150, comme représenté sur les Figure 7. Si l'altitude ZGps au-dessus du géoïde de référence est affichée, le calage «GEO » est affiché dans l'emplacement 150, comme représenté sur la Figure 8. Le symbole 146 est muni d'une identification représentative du type d'altitude affichée, avantageusement un code couleur. Le code couleur choisi indique de manière intuitive quel type de vitesse est affiché. Par exemple, le symbole 146 est vert encadré de vert dans l'état de crédibilité optimale 74, i.e. quand l'altitude Zp ou Zc est affichée, vert encadré de jaune lorsque l'hybridation lâche de Zp** sur Zp est en cours, et jaune encadré de jaune lorsque l'hybridation lâche de Zp** par Zp est arrêtée. Le symbole 146 encadré de jaune signifie donc que l'altitude pression GPS Zp** ou l'altitude au-dessus du géoïde de référence ZGps est affichée, et la couleur du symbole lui-même indique si le processus d'hybridation lâche de Zp** sur Zp est en cours ou non. Le dispositif d'affichage auxiliaire 14 est complété en outre par des moyens 16 d'affichage de messages d'alerte, comme décrit précédemment. Ces moyens 16 comprennent par exemple une fenêtre d'affichage 152, propre à afficher des messages textuels destinés à informer le pilote de l'état de crédibilité affecté aux mesures des capteurs, et de la fiabilité des mesures affichées par les instruments de bord classiques. Cette fenêtre 152 est par exemple intégrée à un système d'alerte de l'équipage (ou CAS pour Crew Alerting System). Par exemple, dans l'état de crédibilité forte, la fenêtre 152 affiche un message indiquant au pilote que les mesures des capteurs de pression 5a, 5b ne sont pas fiables, et un message indiquant que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Ces messages sont par exemple du type « STATIC &/or total pressure input fail » et « M/IAS unreliable (caution) » respectivement, comme illustré sur la Figure 7. Dans l'état de crédibilité faible 76, la fenêtre 152 affiche un message indiquant au pilote que les mesures des capteurs de pression 5a, 5b et d'incidence 5d ne sont pas fiables, et un message indiquant que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Ces messages sont par exemple du type « AOA input fait » et « M/IAS unreliable (caution) » respectivement, comme illustré sur la Figure 8.Le dispositif d'affichage auxiliaire selon l'invention permet donc de fournir au pilote, à tout moment, des informations de pilotage fiables. En effet, lorsque les mesures d'au moins un capteur parmi les capteùrs de pression statique (5a), de pression totale (5b) et d'incidence (5d) sont jugées non fiables, chaque caractéristique de vol déterminée à partir de mesures de ces capteurs est remplacée sur le dispositif d'affichage auxiliaire par une caractéristique de vol auxiliaire. Chaque caractéristique de vol auxiliaire est indépendante de toute mesure de capteur effectuée à un instant auquel ce capteur est jugé non fiable. Les caractéristiques de vol auxiliaires peuvent cependant dépendre de mesures du ou de chaque capteur jugé non fiable effectuées préalablement, lorsque ce ou ces capteurs étaient jugé(s) fiables. Le dispositif d'affichage auxiliaire selon l'invention permet également de prévenir l'équipage en cas de suspicion de défaillance d'un ou plusieurs capteurs. Le procédé et le système selon l'invention permettent ainsi de fournir à l'équipage d'un aéronef une évaluation de la fiabilité des caractéristiques de vol issues de mesures de capteurs de l'aéronef, et de l'alerter en cas de possible défaillance. En outre, ce procédé et ce système permettent de réévaluer en permanence la fiabilité des capteurs pour détecter une éventuelle remise en état de fonctionnement de capteurs jugés non fiables ou au contraire une défaillance d'un capteur auparavant jugé fiable. Le procédé et le système selon l'invention permettent également de fournir à l'équipage des caractéristiques de vol auxiliaires indépendantes des mesures des capteurs jugés non fiables. Ainsi, même en cas de défaillance de capteurs de pression et/ou d'incidence, l'équipage dispose de valeurs fiables des caractéristiques de vol. Il devra être compris que les exemples de réalisation présentés ci-dessus ne sont pas limitatifs. Notamment, la vérification de l'équation de sustentation lors des étapes 70 et 78 peut être réalisée selon une autre forme que celle précédemment décrite. Notamment, elle peut faire intervenir un autre type d'information de vitesse que le nombre de Mach Ma. En outre, les caractéristiques de vol auxiliaires peuvent être déterminées selon d'autres expressions, et éventuellement à partir d'autres types de capteurs qu'un capteur GPS, par exemple un autre type de capteur de position par satellites, un radar Doppler, ou une centrale à inertie. Le coefficient de portance peut également être projeté sur un autre axe que l'axe Z représenté, par exemple selon un axe perpendiculaire au vecteur vitesse. Par ailleurs, le dispositif d'affichage auxiliaire peut afficher d'autres informations auxiliaires, telles que la température statique ou l'estimateur de vent horizontal. Selon une variante, l'aéronef 1 comprend plusieurs capteurs du même type redondants, par exemple plusieurs capteurs de pressions statique et totale et/ou plusieurs capteurs d'incidence. Selon cette variante, le calculateur teste par exemple la fiabilité de chacun des capteurs d'un même type séparément. Ainsi, si un seul capteur d'un type donné est jugé fiable, le calculateur 3 utilise les mesures issues de ce capteur, et avertit le pilote d'une défaillance du ou des autre(s) capteur(s) du même type. Si tous les capteurs d'un même type sont non fiables, le calculateur 3 active un état de crédibilité forte ou faible suivant les cas, conformément au procédé décrit ci-dessus. Bien entendu, d'autres modes de réalisation peuvent être envisagés. En outre, les caractéristiques techniques des modes de réalisation et variantes mentionnées ci-dessus peuvent être combinées entre elles. REVENDICATIONS1.- Procédé de détermination d'une masse estimée (m) d'un aéronef (1) à un instant de détermination (tv) choisi, caractérisé en ce qu'il comprend : - la détermination (21) d'une première masse (mpio) de l'aéronef (1), à partir de caractéristiques de l'aéronef (1) déterminées avant ou pendant le décollage de l'aéronef (1), - la détermination (41), lors du vol de l'aéronef (1), d'une deuxième masse (msp) de l'aéronef (1), à partir d'une équation de sustentation de l'aéronef (1) exprimant la masse de l'aéronef (1) en fonction d'une information représentative du facteur de charge (ni) de l'aéronef (1), d'un coefficient de portance ( ) de l'aéronef (1), d'une information de vitesse (Ma) de l'aéronef (1) et d'une pression statique (Ps) d'une masse d'air traversée par l'aéronef (1), déterminées à partir de mesures de capteurs (5a, 5b, 5c, 5d) de l'aéronef (1) lors du vol de l'aéronef (1), - l'évaluation (39) de la masse estimée (m) audit instant de détermination (tu), en fonction desdites première (mplo) et deuxième (rnsp) masses. 2.- Procédé d'estimation selon la revendication 1, caractérisé en ce que lesdites première (mpio) et deuxième (map) masses sont des masses de l'aéronef (1) à un instant initial (t1) au parking, et en ce que l'étape (39) d'évaluation de ladite masse estimée (tri) comprend la détermination d'une masse estimée au parking audit instant initial (t,) en fonction desdites première (mpio) et deuxième (msp) masses, la détermination d'une masse (FU) de carburant consommé entre ledit instant initial (t,) et ledit instant de détermination (tv), ladite masse estimée ( tri ) étant évaluée en soustrayant ladite masse (FU) de carburant consommé de ladite masse estimée au parking. 3.- Procédé d'estimation selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que la détermination (21) de ladite première masse (m00) comprend l'estimation (23) d'une première masse initiale (mp1) de l'aéronef (1), réalisée avant le décollage, l'estimation (25, 26) d'une deuxième (mp2) et d'une troisième (mp3) masses initiales (me, me) de l'aéronef (1) à partir de caractéristiques de vol de l'aéronef (1) déterminées lors du décollage, ladite première masse (mpio) étant déterminée (27) en fonction desdites première, deuxième et troisième masses initiales (mpl, mp2, me). 4.- Procédé d'estimation selon la revendication 3, caractérisé en ce que ladite première masse (m00) est déterminée à partir d'une médiane desdites première, deuxième et troisième masses initiales (mpl, me, me). 5.- Procédé d'estimation selon l'une quelconque des revendications précédentes, caractérisé en ce que la détermination (21) de la première masse (mplo) de l'aéronef (1) comprend l'estimation (23) d'une première masse initiale (mp1) de l'aéronef (1) à partir d'une estimation d'une masse de base de l'aéronef (1), d'une masse de carburant de l'aéronef (1) au parking et d'une masse du chargement de l'aéronef (1). 6.- Procédé d'estimation selon l'une quelconque des revendications précédentes, caractérisé en ce que la détermination (21) de la première masse (mp,o) de l'aéronef (1) comprend l'estimation (25) d'une deuxième masse initiale (mp2) de l'aéronef (1), ladite deuxième masse initiale (mp2) étant déduite à partir d'une valeur de poussée (Fx) des moteurs, d'une valeur de traînée (Tx) de l'aéronef (1) et d'une valeur d'accélération (Jx) de l'aéronef (1) déterminées pendant le décollage de l'aéronef (1), à un instant (t2) choisi, notamment 3 secondes après l'envol de l'aéronef (1). 7.- Procédé d'estimation selon l'une quelconque des revendications précédentes, caractérisé en ce que la détermination (21) de la première masse (mp,o) de l'aéronef (1) comprend l'estimation (26) d'une troisième masse initiale (mp3) de l'aéronef (1), à partir d'une équation de sustentation de l'aéronef (1) exprimant la masse de l'aéronef (1) en fonction d'une information représentative du facteur de charge (ni) de l'aéronef (1), d'un coefficient de portance (az* ) de l'aéronef (1), d'une information de vitesse (M**) de l'aéronef (1) et d'une pression statique (Ps) d'une masse d'air traversée par l'aéronef (1), déterminées pendant le décollage de l'aéronef (1), à un instant (t2) choisi, notamment 3 secondes après l'envol de l'aéronef (1). 8.- Procédé d'estimation selon l'une quelconque des revendications précédentes, caractérisé en ce que la détermination (41) de ladite deuxième masse (msp) de l'aéronef (1) comprend : - la détermination d'une pluralité de masses instantanées (m4) de l'aéronef (1) à des instants de vol (t) distincts, lesdites masses instantanées (m4) étant déterminées à partir de ladite équation de sustentation de l'aéronef (1), - la détermination à partir desdites masses instantanées et de masses de carburant consommé entre un instant initial (t1) et lesdits instants de vol (t), d'une pluralité de quatrièmes masses initiales (m4') de l'aéronef (1) audit instant initial (t1), - la détermination, à partir desdites quatrièmes masses initiales (m4'), de ladite deuxième masse (msp) de l'aéronef (1). 9.- Procédé d'estimation selon la revendication 8, caractérisé en ce que la détermination de ladite deuxième masse (msp) de l'aéronef (1) comporte le calcul d'une moyenne pondérée desdites quatrièmes masses initiales (m4') en attribuant à chaquequatrième masse initiales (m4') un coefficient de pondération fonction d'une crédibilité accordée à la valeur de cette masses initiale (m4'), sous la forme : t, fx,(t)m41(t) msp ) = t2 , n, (t) fx t, où rnsp(t, ) désigne la masse de l'aéronef à l'instant initial (t1) déterminée à l'instant de détermination tv, chaque m4'(t) désigne une quatrième masse initiale déterminée à un instant t, et chaquexm(t) est un coefficient de pondération attribué à la quatrième masse initiale m4'(t). 10.- Procédé d'estimation selon la revendication 9, caractérisé en ce que le coefficient de pondération (xm(t)) attribué à une quatrième masse initiale (m4'(t)) déterminée à un instant t est égal à un lorsqu'à l'instant t, l'aéronef (1) a sensiblement un vol rectiligne horizontal selon des critères prédéterminés, ou à zéro dans le cas contraire. 11.- Procédé de détermination selon l'une quelconque des revendications précédentes, caractérisé en ce que ladite masse estimée (ir-i) est évaluée à partir d'une moyenne pondérée desdites première (mp,o) et deuxième (msp) masses. 12.- Procédé d'estimation selon la revendication 11, caractérisé en ce que ladite moyenne pondérée desdites première et deuxième masses est réalisée en attribuant à chacune desdites première (mp,o)et deuxième (msp) masses un coefficient de pondération, le coefficient de pondération attribué à ladite deuxième masse (msp) étant proportionnel à la durée de vol rectiligne uniforme de l'aéronef depuis un instant initial (t1). 13.- Procédé d'estimation selon la revendication 12, caractérisé en ce que le coefficient de pondération attribué à ladite deuxième masse (msp) est exprimé sous la forme : tv iZm(t) p(tv)= _ti , où p(ti désigne le coefficient de pondération attribué à ladite deuxième masse (msp) à l'instant de détermination tv, t1 est l'instant initial, etzn(t) une fonction égale à un lorsqu'à l'instant t, l'aéronef (1) a sensiblement un vol rectiligne uniforme, ou à zéro dans le cas contraire. 14.- Système de détermination d'une masse estimée (in) ) d'un aéronef (1) à un instant de détermination (tv) choisi, caractérisé en ce qu'il comprend :- des moyens configurés pour déterminer une première masse (mp,o) de l'aéronef (1), à partir de caractéristiques de l'aéronef (1) déterminées avant ou pendant le décollage de l'aéronef (1), - des moyens configurés pour déterminer, lors du vol de l'aéronef (1), une deuxième masse (msp) de l'aéronef (1), à partir d'une équation de sustentation de l'aéronef (1) exprimant la masse de l'aéronef (1) en fonction d'une information représentative du facteur de charge (ni) de l'aéronef (1), d'un coefficient de portance (az ) de l'aéronef (1), d'une information de vitesse (Ma) de l'aéronef (1) et d'une pression statique (Ps) d'une masse d'air traversée par l'aéronef (1), déterminées à partir de mesures de capteurs (5a, 5b, 5c, 5d) de l'aéronef (1) lors du vol de l'aéronef (1), - des moyens configurés pour évaluer la masse estimée (i5) audit instant de détermination (tu), en fonction desdites première (mp,o) et deuxième (msp) masses.15