FR 2998702 A1 20140530 FR 1203197 A 20121127 Procédé de déteimination dynamique de la position du point d'arrêt d'un aéronef sur un terrain, et système associé La présente invention concerne un procédé de détermination dynamique de la position du point d'arrêt d'un aéronef sur un terrain, comprenant les étapes suivantes : - mesure dynamique de la position géographique courante de l'aéronef, et de sa vitesse sol courante ; - calcul d'une première longueur représentative d'une phase de vol de l'aéronef, et calcul d'une deuxième longueur représentative d'une phase au sol de l'aéronef ; - calcul de la position du point d'arrêt sur la base de la première longueur, de la deuxième longueur et de la position géographique courante. Un tel procédé est destiné à être mis en oeuvre, lors de l'approche de l'aéronef vers un terrain (par exemple une piste d'atterrissage) et lors de l'atterrissage sur ce terrain. La phase d'atterrissage d'un aéronef est une phase critique du vol. En effet, lors de cette phase, le pilote doit parvenir à arrêter l'aéronef sur la piste, ou doit engager une procédure spécifique d'urgence, s'il estime qu'il ne sera pas en mesure d'arrêter l'aéronef avant l'extrémité de la piste. La procédure spécifique est par exemple une procédure de remise des gaz si la vitesse de l'avion et la longueur de piste restante l'y autorise. À cet égard, les études statistiques des accidents aériens survenus durant la dernière décennie sur des vols commerciaux montrent qu'une majorité des accidents ont lieu lors de l'atterrissage, et qu'environ 18 % des accidents ont pour cause principale une sortie de piste. Une part importante de ces accidents est due, lors de la phase d'atterrissage, à des sorties de pistes en longueur. Plusieurs causes sont isolées pour ce type d'accidents, mais dans la majeure partie des cas, l'équipage avait mal estimé la capacité à faire atterrir l'aéronef sur la piste. Pour pallier ce problème, des systèmes et des procédés ont été développés pour permettre au pilote, lorsqu'il est en vol et/ou au sol, d'estimer s'il sera en mesure de poser l'aéronef en toute sécurité, et d'estimer les capacités de freinage de son avion sur la piste considérée. L'estimation de la position d'arrêt de l'avion, présentée à l'équipage durant l'approche et jusqu'à l'arrêt complet de l'aéronef, est donc une information qui présente un fort intérêt en terme de sécurité. Cet intérêt est renforcé si l'estimation tient compte de l'évolution de la situation en temps réel. Dans un procédé connu, une distance minimale d'atterrissage le long d'une piste, jusqu'à l'arrêt de l'aéronef est calculée à partir de la position courante de l'aéronef, d'un plan de descente standard vers la piste d'atterrissage déterminée, et d'un modèle de performance de l'aéronef. Le calcul de la distance minimale d'atterrissage comprend l'estimation d'une distance parcourue au-dessus de la piste en vol, fondée sur un temps moyen de vol au- dessus de la piste, et l'estimation d'une distance de roulage sur la piste. Ce procédé considère uniquement le plan de descente constant de l'aéronef. Un but de l'invention est donc d'obtenir un procédé de détermination du point d'arrêt d'un aéronef au sol, pouvant être mis en oeuvre lors d'une phase d'approche et lors du roulage, et dont la précision est améliorée. À cet effet, l'invention a pour objet un procédé du type précité, caractérisé en ce que le calcul de la première longueur comprend une phase d'estimation dynamique d'une distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef. Le procédé selon l'invention peut comprendre l'une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toute combinaison techniquement possible : - il comprend une étape de mesure de la pente courante de l'aéronef, la distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef étant calculée en fonction de la pente courante de l'aéronef ; - lors de la phase d'estimation, la distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef est estimée en fonction de la vitesse sol courante de l'aéronef, mesurée lors de l'étape de mesure ; - lors de la phase d'estimation, la distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef est estimée en fonction du facteur de charge nzarr de l'aéronef, estimé durant son arrondi ; - lors de la phase d'estimation, la distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef est estimée par l'équation: gamm.(0-/i 1 D2(t) = fo- [ ,xV s'oi(t)x cos(y + dy)1, dy g x za' - co s(y) ) 2 où g est l'accélération de la pesanteur, 13 est la pente du terrain, nza, est le facteur de charge de l'aéronef estimé durant l'arrondi, Vsol(t) est la vitesse sol courante de l'aéronef, et gamma(t) est la pente courante de l'aéronef ; - la phase d'estimation comporte le calcul d'une hauteur estimée de l'aéronef à laquelle démarre l'arrondi, la distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef étant calculée à partir de la hauteur estimée ; - le calcul de la première longueur comporte une phase d'estimation d'une distance parcourue lors de l'approche de l'aéronef vers le terrain avant l'arrondi et/ou lors d'une première partie de la phase d'arrondi de l'aéronef, la distance parcourue lors de l'approche de l'aéronef vers le terrain avant l'arrondi et/ou lors d'une première partie de la phase d'arrondi de l'aéronef étant calculée en fonction de la pente courante mesurée de l'aéronef, et avantageusement, en fonction de la hauteur courante de l'aéronef par rapport au terrain ; - l'étape de calcul de la deuxième longueur comporte une phase de calcul d'une distance représentative du passage de la décélération courante à la décélération correspondant à un freinage établi de l'aéronef, en fonction de la décélération courante ; - la distance représentative du passage de la décélération courante à la décélération correspondant au freinage établi est calculée en fonction d'une constante Km représentative de la dynamique de l'aéronef entre le toucher du train principal et l'établissement du freinage établi ; - l'étape de calcul de la deuxième longueur comporte une phase de calcul d'une distance de freinage à partir d'un profil de décélération prédéterminé de l'aéronef sur le terrain, en fonction de la vitesse sol courante ; - le profil de décélération prédéterminé est une décélération constante ; - l'étape de calcul du point d'arrêt comporte la détermination de la position géographique du point d'arrêt sur la base de la position géographique courante de l'aéronef, de la première longueur, de la deuxième longueur, et d'un relèvement correspondant à la route de l'aéronef. L'invention a également pour objet un système de détermination dynamique de la position du point d'arrêt d'un aéronef sur un terrain, comprenant : - un ensemble de mesure dynamique de la position géographique courante de l'aéronef, et de sa vitesse sol courante ; - un ensemble de calcul d'une première longueur représentative d'une phase de vol de l'aéronef, et un ensemble de calcul d'une deuxième longueur représentative d'une phase au sol de l'aéronef ; - un ensemble de calcul de la position du point d'arrêt sur la base de la première longueur, de la deuxième longueur et de la position géographique courante ; caractérisé en ce que l'ensemble de calcul de la première longueur comprend des moyens d'estimation d'une distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef. Le système selon l'invention peut comprendre la caractéristique suivante : - un ensemble de mesure de la pente courante de l'aéronef, les moyens d'estimation étant propres à calculer la distance parcourue par l'aéronef lors d'au moins une partie d'une phase d'arrondi de l'aéronef en fonction de la pente courante de l'aéronef. 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 annexés, sur lesquels : - la figure 1 est une vue schématique d'un aéronef comportant un système destiné à la mise en oeuvre du procédé selon l'invention ; - la figure 2 est une vue dans le plan de la trajectoire de l'aéronef lors de son approche vers le sol, avec l'aéronef représenté dans sa phase d'approche avant le début de l'arrondi, pour illustrer un premier mode de réalisation du procédé selon l'invention ; - la figure 3 est une vue analogue à la figure 2, lorsque l'aéronef amorce une phase d'arrondi ; - la figure 4 est un logigramme représentant les différentes étapes du procédé selon l'invention ; - la figure 5 est une vue analogue à la figure 2, pour un deuxième mode de réalisation du procédé selon l'invention ; et - la figure 6 est une vue analogue à la figure 3 pour le deuxième mode de réalisation du procédé selon l'invention. Sur la Figure 2, l'évolution d'un aéronef 12 entre une position courante en vol et une position d'arrêt 16 sur un terrain 14 comporte une phase en vol (pv tant que l'aéronef 12 est en sustentation et jusqu'à son toucher des roues sur le sol et une phase au sol cps, du toucher des roues sur le sol jusqu'à l'arrêt complet de l'aéronef 12. La phase en vol (pv comporte une première phase (pl dans laquelle l'aéronef 12 s'approche du terrain 14 selon une première trajectoire, par exemple une trajectoire rectiligne à pente constante puis une deuxième phase (p2 dite d'arrondi, au cours de laquelle le pilote réalise la transition entre la phase de sustentation et la phase au sol en réduisant la pente de l'aéronef 12 pour diminuer sa vitesse verticale et assurer le confort passager durant cette transition. La phase au sol (ps comporte une première phase (p3 au cours de laquelle la décélération de l'aéronef 12 varie jusqu'à atteindre une valeur correspondant au freinage maximum (par exemple une valeur constante) puis une deuxième phase (p4 au cours de laquelle la décélération de l'aéronef 12 est la décélération associée au freinage maximum, jusqu'à l'arrêt complet de celui-ci à sa position d'arrêt 16. Un premier système 10 de détermination dynamique de la position du point d'arrêt d'un aéronef 12 selon l'invention est illustré schématiquement par la figure 1. En référence à la figure 2, le système 10 est destiné à déterminer, à tout instant t pendant l'atterrissage de l'aéronef 12, vers le terrain 14, la position géographique du point d'arrêt 16 de l'aéronef 12 sur le terrain 14. Comme on le verra plus bas, la position géographique du point d'arrêt 16 est déterminée par ses coordonnées géographiques, par exemple par sa latitude Lpos Arrêt (t) et sa longitude GPos Arrêt (t). Cette détermination est faite avantageusement indépendamment de la position d'une piste d'atterrissage 18 représentée schématiquement sur la figure 2, la piste d'atterrissage 18 présentant un seuil amont 20 et un seuil aval 22. Le système 10 est destiné à fournir, en temps réel, à l'équipage de l'aéronef 12 une indication sur la position estimée du point d'arrêt 16, en fonction des paramètres courants de vol de l'aéronef 12, et notamment de sa position courante (LNc(t) GNc(t)), de sa décélération longitudinale courante par rapport au sol Jx(t), de sa vitesse sol courante Vsol(t), de sa pente sol courante gamma(t), et de sa hauteur courante H(t) par rapport au terrain 14. Outre le système 10, l'aéronef 12 comporte, en référence à la figure 1, un système 30 de mesure et de positionnement, et un système 32 de surveillance de pannes des systèmes avion. Il comprend en outre un dispositif d'affichage 34 d'informations relatives au point d'arrêt 16 de l'aéronef 12 sur le terrain 14. Le système 30 de mesure et de positionnement comporte par exemple des capteurs de mesure de paramètres liés à l'aéronef et à son environnement, tels que la température, la pression ou la vitesse, sa position, son altitude et ses attitudes. Le système 32 est propre à surveiller et à déterminer notamment l'état actuel des systèmes de l'aéronef 12 (notamment les systèmes de freinage), tenant compte des pannes sur ces systèmes. Le système de détermination 10 est raccordé au système de mesure et de positionnement 30 et au système de surveillance 32. Il comporte au moins un processeur 40 et une mémoire 42. Selon l'invention, le système 10 comprend un premier ensemble 52 de calcul d'une première longueur LV(t) représentative de la phase de vol (pv de l'aéronef 12 et un deuxième ensemble 54 de calcul d'une deuxième longueur LS(t) représentative de la phase au sol cf), de l'aéronef 12. Le système 10 comporte en outre un troisième ensemble 56 de calcul du point d'arrêt 16, sur la base de la première longueur LV(t), de la deuxième longueur LS(t), et d'au moins la position géographique courante de l'aéronef 12. Le système 10 est raccordé à une base de données topographique 57 et à une base de données 58 de performances de l'aéronef 12. La base de données 57 comporte des données topographiques relatives au terrain 14 en regard duquel évolue l'aéronef 12. Elle comporte avantageusement des informations relatives aux pistes d'atterrissage 18 présentes sur le terrain 14, et notamment relatives à la position géographique du seuil 20, 22 de chaque piste 18. La base de données de performances 58 comporte des tables relatives aux paramètres de performances de l'aéronef 12 lors de son arrondi (facteurs de charge nzarr) et lors de son freinage au sol. Dans un mode de réalisation particulier, le facteur de charge nzarr est considéré constant lors de l'arrondi. Elle comporte en outre une table de données de performance de freinage établissant l'évolution de la décélération J -x tab de l'aéronef 12 en fonction de sa vitesse sol courante V',(t). Cette évolution est également dépendante de la configuration de l'aéronef (notamment de la position des volets), des performances de freinage de l'aéronef 12, et éventuellement, de l'état de la piste (piste mouillée, piste sèche, piste contaminée). . Dans un mode de réalisation particulier, cette décélération Jx tab est supposée constante en fonction du temps et égale à la décélération constante Jx max correspondant au freinage maximal de l'aéronef 12. La décélération correspondant au freinage maximal de l'aéronef Jx max J est fonction de la configuration de l'aéronef 12, des performances de freinage de l'aéronef 12, et de l'état de la piste (piste mouillée, piste sèche, piste contaminée). La base de données de performances 58 comprend en outre une table de données de constantes Km représentant la dynamique de l'aéronef 12 entre le toucher du train principal d'atterrissage et l'établissement du freinage maximal. Le système de positionnement et de mesure 30 comporte un premier ensemble 60 de mesure dynamique de la position géographique courante (Levc(t) ; GNdt) ; H(t)) de l'aéronef 12, un deuxième ensemble 62 de mesure de la vitesse sol Vsol(t) courante de l'aéronef 12, et un troisième ensemble 64 de mesure de la pente courante gamma(t) de l'aéronef 12. Il comporte en outre un quatrième ensemble 66 de mesure de la décélération courante Jx(t) de l'aéronef 12. Le premier ensemble 60 est raccordé aux capteurs de position. Il est apte à déterminer dynamiquement, sur la base des données reçues des capteurs de position, et sur la base d'informations reçues de la base de données géographique 57, la position courante de l'aéronef, à savoir sa latitude LNc(t), sa longitude GNc(t), et sa hauteur H(t) par rapport au terrain 14. Le deuxième ensemble 62 est également raccordé aux capteurs. Il est apte à déterminer dynamiquement la vitesse sol courante V',(t) de l'aéronef. Le troisième ensemble 64 est raccordé à des capteurs de position. Il est apte à déterminer dynamiquement, sur la base des données reçues des capteurs, la pente sol avion gamma(t). Le quatrième ensemble 66 est aussi raccordé aux capteurs. Il est apte à déterminer dynamiquement la décélération courante Jx(t) de l'aéronef 12. Le premier ensemble de calcul 52 est raccordé aux ensembles de mesure 60 à 66. Dans un premier mode de réalisation illustré sur les figures 2 et 3, il comporte des moyens d'estimation 70 d'une première distance D1(t) correspondant à la distance séparant la projection 701 sur le terrain 14 de l'aéronef 12 en vol, et le plan vertical 703 passant par l'intersection 705 du plan de descente courant 707 avec le terrain 14. Il comporte en outre des moyens 72 d'estimation d'une deuxième distance D2(t) correspondant à la distance séparant le plan vertical 703 du point de toucher des roues 709 de l'aéronef sur le terrain 14. Dans ce mode de réalisation, et par rapport à une trajectoire supposée de l'aéronef 12, tant que l'aéronef 12 n'a pas démarré la phase d'arrondi (p2, la distance D1(t) est représentative de la phase d'approche (pi et d'une première partie de la phase d'arrondi (p2 tandis que la distance D2(t) est représentative de la seconde partie de la phase d'arrondi (p2. Une fois que l'aéronef 12 a amorcé la phase d'arrondi (p2, la distance D1(t) est représentative de la première partie de la phase d'arrondi (p2 tandis que la distance D2(t) est représentative de la seconde partie de la phase d'arrondi cp2. Dans un mode de réalisation particulier, le procédé selon l'invention consiste donc à assimiler D2(t) à la projection d'un demi arc de cercle 0(0 / 2 de l'arrondi, l'arrondi étant supposé s'effectuer suivant une trajectoire circulaire. 0(t) et gamma(t) sont liés par la relation géométrique suivante : 0(t) = gamma(t) - f3 ; p corresponsant à la pente de la piste 18 ((3 étant égal à 0 pour une piste horizontale). Les premiers moyens d'estimation 70 sont propres à calculer dynamiquement la distance D1(t) en fonction de la hauteur courante H(t) de l'aéronef 12 et de sa pente courante gamma(t), par exemple par l'équation (1) ci-dessous : H(t) D1(t) = (1) tan[gamma(t)] En particulier, lorsque la pente courante est faible, l'équation (2) ci-dessous peut être réécrite sous la forme : H(t) D1(t) = r , (2) [gamma(t)] Lorsque la pente courante est très proche de zéro, gamma(t) est alors bornée à une valeur seuil gammamin. Les deuxième moyens d'estimation 72 sont propres à calculer dynamiquement la distance D2(t) en fonction de la vitesse sol courante Vaol(t), de la pente avion courante gamma(t), et du facteur de charge nzarr de l'aéronef 12, obtenu à partir de la base de données de performances 58 et de l'accélération de la pesanteur g, par exemple par l'équation (3) ci-dessous : gamma(t)-13 1D2(t)=J 2 X V2/ (t) X COS( y +dyI dy (3) g x(nza' - cos(y)) 2 En particulier, lorsque l'arrondi est supposé être circulaire, le facteur de charge nzarr est constant. Avec l'hypothèse des petits angles, l'équation (3) peut être simplifiée sous la forme de l'équation suivante (4) : D2(t) = 1 xV2 ,(t)x (gamma(t) - fi) (4) 2 x g x(nzarr La prise en compte effective de l'arrondi offre ainsi une plus grande précision dans l'estimation de la position d'arrêt de l'aéronef que celles des procédés d'estimation de l'art antérieur. Le deuxième ensemble de calcul 54 comporte des troisièmes moyens d'estimation 74 de la distance D3(t) nécessaire au passage, dans la phase de freinage au sol, de la décélération courante Jx(t) à la décélération correspondant à un freinage établi, notamment au freinage maximal J -x frmax- Le deuxième ensemble de calcul 54 comporte en outre des quatrièmes moyens d'estimation 76 de la distance D4(t) parcourue au sol par l'aéronef 12 en freinage établi. Les troisièmes moyens d'estimation 74 sont aptes à calculer dynamiquement la distance D3(t) en fonction de la décélération correspondant à un freinage établi, notamment au freinage maximal J -x frmax, de la décélération courante Jx(t), et d'une constante Kms représentant la dynamique de l'aéronef 12 entre le toucher du train d'atterrissage principal et l'établissement du freinage établi, tels qu'obtenue à partir de la base de données de performances 58. En particulier, la distance D3(t) est estimée par l'équation (5) suivante : D3(t) = Kms x (.1x trmax - Js (t)) (5) Les quatrièmes moyens d'estimation 74 sont aptes à calculer dynamiquement la distance D4(t) en fonction du profil estimé Jx tab de décélération de l'aéronef 12 sur le sol, tel qu'il est tabulé dans la base de données de performances 58. Ce profil estimé .1 -x tab est tabulé en fonction de la vitesse sol courante V',(t). La distance D4(t) est obtenue dynamiquement par une intégration de la fonction correspondant à ce profil J -x tab entre l'instant courant t et le temps d'arrêt t -arrêt de l'aéronef 12. Les moyens 74 comprennent, à cet effet, des moyens de calcul du temps d'arrêt tarrêt de l'aéronef. Ce temps d'arrêt t .arrêt est estimé à l'aide du profil de décélération tabulé Jx(Vsol(t),t) correspondant à la vitesse sol courante Vsol(t). Il est par exemple obtenu par le système d'équations (6) suivant : t _t a(rrêt> 0 1 arrêt x tab Vsol(t),t). d Dans le cas où le profil de décélération est constant, et est par exemple égal à la décélération correspondant au freinage maximal .1 -x frmax, le temps d'arrêt t -arrêt est obtenu par l'équation (7) suivante : Vs° (t) t arrêt = i (7) x fr max La distance D4 est alors obtenue par l'équation (8) suivante : f Jtarrêt [tarrêt D4(t) = x _tab(Vsol(t)5t) - dtl - dt (8) Dans le cas où le profil de décélération est constant, et est par exemple égal à la décélération correspondant au freinage maximal .1 -x frmaxe la distance D4(t) est calculée par l'équation (9) suivante : 2 D4(t) = V"1(t) (9) 2 xx fr max L'estimation de la distance horizontale L(t) séparant la projection au sol V de la position actuelle de l'aéronef de sa position d'arrêt estimée s'exprime ainsi : L(t) = LV (t) + LS(t) = Dl(t) + D2(t) + D3(t) + D4(t) (10) Le troisième ensemble de calcul 56 est apte à calculer les coordonnées géographiques (Lpos are) ; Gpos arrêt(t)) du point d'arrêt 16, en fonction des coordonnées géographiques courantes de l'aéronef 12 (LNc(t) ; GNc(t)), du relèvement x (t), de =0 (6) l'estimation de la distance horizontale L(t) entre l'aéronef 12 et son point d'arrêt 16 et du rayon moyen terrestre RO, avec l'hypothèse d'une terre sphérique. Le relèvements (t) est l'angle formé dans le plan horizontal par la ligne d'un observateur vers un objet par rapport au nord géographique. Il est avantageusement pris égal à la route sol avion, c'est-à-dire l'angle formé par le vecteur vitesse sol de l'avion avec le nord géographique.. La formule suivante est avantageusement utilisée pour calculer les coordonnées géographiques (Lpos arrêt(t) ; Gpos arrêt(t)) P-posArret(t)1 LGposArret (t)i Arc sin[sin(ilx coe Ac (t))x cos(x(t))+ sin(L Ac (t))x cos(11 RO RO sin(x(t))x sin(LR(Ot)) cos(41)- RO x coeAc (t))- cos(x(t))x sieAc (t))x sin(R(0) GAc (0+ Arc tan (11) Le dispositif d'affichage 34 est par exemple disposé dans le cockpit de l'aéronef 12. Il comporte généralement un afficheur 80, tel qu'un écran et un ensemble 82 de gestion de l'affichage sur l'afficheur 80. Il est raccordé au système 10, et notamment au troisième ensemble de calcul 56. L'ensemble de gestion 82 est propre à extraire des informations reçues du système 10, et notamment la position du point d'arrêt 16 estimée dynamiquement par le système 10 pour permettre l'affichage d'une information représentative de ce point d'arrêt 16 sur l'afficheur 80. Cette information peut être un symbole, une ligne, une jauge. L'information est par exemple affichée sur une représentation géographique du terrain 14 présent à l'avant de l'aéronef 12, et/ou sur une représentation d'une piste d'atterrissage 18 et de son environnement sur l'afficheur 80. Un premier procédé de détermination dynamique de la position du point d'arrêt 16 d'un aéronef 12 sur un terrain 14 va maintenant être décrit, en regard de la figure 3. Ce procédé est mis en oeuvre dans un système 10 selon l'invention, tel que décrit plus haut. Le procédé selon l'invention comporte une étape 100 de mesure dynamique de paramètres courants de l'aéronef 12, une étape 102 de calcul d'une première longueur LV(t) représentative d'une phase de vol de l'aéronef 12, une étape 104 de calcul d'une deuxième longueur LS(t) représentative d'une phase au sol de l'aéronef 12, et une étape 106 de calcul de la position du point d'arrêt 16 sur la base de la première longueur LV, de la deuxième longueur LS(t) et de la position géographique courante de l'aéronef 12. Le procédé comporte en outre avantageusement une étape 108 d'affichage d'une information représentative de la position du point d'arrêt 16 sur un afficheur 80 du dispositif d'affichage 34. À l'étape 100, l'ensemble 60 de mesure dynamique de la position géographique courante mesure, sur la base des données reçues des capteurs de position, et sur la base des informations reçues de la base de données géographique 57, la position courante de l'aéronef, à savoir sa latitude LA,c(t), sa longitude GA/c(t), et sa hauteur H(t) par rapport au terrain 14. L'ensemble 62 mesure dynamiquement la vitesse sol courante V',(t) de l'aéronef 12. L'ensemble 64 détermine dynamiquement, la pente avion courante gamma(t). Ainsi, à chaque instant t, la position géographique courante de l'aéronef 12, et sa pente courante gamma(t) sont mesurées. L'étape 102 comporte une première phase 110 d'estimation d'une distance D1(t) parcourue lors de l'approche de l'aéronef 12 vers le terrain 14, et/ou lors d'une première partie de l'arrondi, en fonction de la pente courante mesurée gamma(t) de l'aéronef 12, et, selon l'invention, une deuxième phase 112 d'estimation d'une distance D2(t) parcourue lors d'une deuxième partie de l'arrondi de l'aéronef 12, au-dessus du terrain 14, en fonction de la pente courante gamma(t) de l'aéronef 12. Lors de la phase 110, la distance D1(t) parcourue lors de l'approche de l'aéronef 12 vers le terrain 14, au-dessus du terrain 14, avant d'effectuer un arrondi, et/ou lors d'une première partie de l'arrondi, est calculée en fonction de la hauteur courante H(t) de l'aéronef 12, et de sa pente courante gamma(t) par les premiers moyens d'estimation 70. Ce calcul est effectué par exemple par l'une des équations (1) ou (2) ci-dessus. Puis, lors de la phase 112, la distance D2(t) parcourue en vol par l'aéronef, lors d'une deuxième partie de l'arrondi de l'aéronef 12 est calculée, par les deuxièmes moyens d'estimation 72, en fonction de la vitesse sol courante Vsol(t), de la pente avion courante gamma(t), et du facteur de charge nzarr de l'aéronef 12 lors de l'arrondi, obtenu à partir de la base de données de performances 58, par exemple par l'une des équations (3) ou (4) ci-dessus. La première longueur LV(t) correspond à la somme des distances D1(t) et D2(t) obtenues lors des phases 110, 112. L'étape 104 comporte une troisième phase 114 de calcul de la distance D3(t) nécessaire au passage, dans la phase de freinage au sol, de la décélération courante Jx(t) à la décélération correspondant à un freinage établi -I -x maxi en fonction d'une constante Km, représentant la dynamique de l'aéronef 12 entre le toucher du train d'atterrissage principal et l'établissement du freinage établi. Le calcul de la distance D3(t) est en outre effectué en fonction de la décélération courante Jx(t), et de la décélération correspondant à un freinage établi J -x max. Cette distance D3(t) est calculée par les troisièmes moyens d'estimation 74, en particulier par l'équation (5). L'étape 104 comporte en outre une quatrième phase 116 de calculardr:une distance D4(t) de roulage en freinage établi, sur la base du profil estimé Jx(Vsol(t),t) de décélération de l'aéronef 12 sur le terrain 14, et d'une détermination du temps d'arrêt t de l'aéronef 12. À cet effet, un calcul du temps d'arrêt -arrêt de l'aéronef 12 est effectué par les moyens d'estimation 74, à l'aide de l'une des équations (6) ou (7). La distance D4(t) est obtenue par une intégration du profil estimé Jx(V,01(t),t) de décélération de l'aéronef 12 sur le terrain 14, avantageusement par l'une des équations (8) ou (9). La deuxième longueur LS(t) correspond à la somme des distances D3(t) et D4(t) obtenues lors des phases 114,116. Lors de l'étape 106, le troisième ensemble de calcul 56 calcule les coordonnées géographiques Lpos are) ; Gpos arrêt(t) du point d'arrêt 16, en fonction des coordonnées géographiques courantes LNdt) ; GA/c(t) de l'aéronef 12, du relèvement x (t) pris égal à la route avion, et de la distance horizontale L(t) entre l'aéronef 12 et le point d'arrêt 16, obtenue en effectuant la somme de la première longueur LV(t) et de la deuxième longueur LS(t) obtenues respectivement aux étapes 102 et 104. Comme précisé plus haut, cette position est calculée par une formule classique, avec l'hypothèse d'une terre sphérique. À l'étape 108, l'ensemble de gestion 82 du dispositif d'affichage 34 extrait les informations reçues du système 10, et notamment la position du point d'arrêt 16, et affiche en dynamique sur l'afficheur 80, une information représentative de la position de ce point d'arrêt 16, tel que décrite plus haut. La mise en oeuvre du procédé selon l'invention améliore grandement la précision de l'estimation de la position du point d'arrêt 16 de l'aéronef 12 sur le terrain 14, en tenant compte, en temps réel, de l'évolution de l'aéronef 12. En particulier, l'utilisation de la pente avion à chaque instant améliore l'estimation de la distance D2(t) parcourue lors de l'arrondi de l'aéronef 12 et contribue à augmenter la précision de l'estimation. Dans une variante, le système 10 comporte, en complément ou en remplacement du dispositif d'affichage 34, un dispositif d'alarme propre à émettre une alarme sonore au pilote, par exemple lorsque la position d'arrêt estimée de l'aéronef est située au-delà de la position du seuil de piste 22. Dans une variante, le système 10 recueille des informations dans le système avionique 32 et/ou dans les bases de données 58 pour déterminer des paramètres supplémentaires influençant les performances de freinage de l'avion, tel que le fonctionnement des systèmes avion, en particulier du système de freinage, et/ou le fonctionnement et la configuration des éléments de sustentation. Dans une autre variante, lorsque l'approche vers la piste d'atterrissage 18 n'est pas rectiligne, le système 10 est apte à modéliser la position du point d'arrêt 16, sur la base d'une estimation de la projection courbe de la trajectoire d'approche sur le plan horizontal. Dans encore une autre variante, le système 10 tient compte de l'environnement de l'aéronef 12, et notamment du relief présent sur le terrain 14, en particulier pour corriger la hauteur courante H(t) mesurée et la ramener à la hauteur par rapport au seuil de piste 20. Dans une autre variante, le système 10 comprend des moyens d'estimation en temps réel de l'opportunité pour le pilote de l'aéronef 12 de faire une manoeuvre de remise des gaz, sur la base de la position du point d'arrêt 16 déterminé. Un autre mode de réalisation de l'invention est illustré sur les figures 5 et 6. Dans ce mode de réalisation, les moyens d'estimation 70 estiment une première distance D1(t) correspondant à la distance séparant le projection 701 sur le terrain 14 de l'aéronef en vol et le plan vertical 713 passant par le point 711 de la trajectoire de l'aéronef où celui-ci amorce la phase d'arrondi (P2. Les moyens d'estimation 72 estiment une deuxième distance D2(t) correspondant à la distance séparant le plan vertical 713 du point de toucher des roues 709 de l'aéronef sur le terrain 14. Dans ce mode de réalisation, la distance D1(t) est représentative de la phase d'approche (pi tandis que D2(t) est représentative de toute la phase d'arrondi (p2. Pour ce faire, les moyens d'estimation 70 sont propres à calculer une hauteur AH(t) par rapport à la hauteur du seuil de piste 20 à partir de laquelle l'arrondi sera effectué. La hauteur AH(t) est calculée en fonction de la vitesse sol courante Vs,,i(t), de la pente avion courante gamma(t), et du facteur de charge nzarr de l'aéronef 12, obtenu à partir de la base de données de performances 58, et de l'accélération de la pesanteur g, par exemple par l'équation (12) ci-dessous : o V:sol (t)2 1 AI-1(t) =sin' (y)- dy (12) fgamma(t)--fl g 112 a' - COS cos(y) La distance Dl (t) est ensuite calculée par l'équation suivante : D1(t)= h(t)- AH (13) gamma(t) La distance Dl (t) est écrêtée inférieurement par 0, pour assurer Dl (t)=0 pendant l'arrondi. Ensuite, la distance D2(t) prenant en compte l'arrondi est calculée à chaque instant en fonction de la hauteur AH(t) par l'équation suivante : 2 2 D2(t) = [2. V", (t) 1 sin(gaminal -[3,1e (14) soit g n z a' - cos(gamma(t)) 2 D2(t) = [2 Vsd(t) 2 1sin(gamma(t)I 2o v soi (t) 2 1 sin(y).dyl 2 1 [r g n z a' - cos(gamma(t)) 2 ) jgamma" g n z a'. -cos(y) (15) Les formules de calcul de D3(t) et D4(t) ainsi que de la position d'arrêt restent elles inchangées par rapport au mode de réalisation décrit précédemment. REVENDICATIONS1.- Procédé de détermination dynamique de la position du point d'arrêt (16) d'un aéronef (12) sur un terrain (14), comprenant les étapes suivantes : - mesure (100) dynamique de la position géographique courante de l'aéronef (12), et de sa vitesse sol courante ; - calcul (102) d'une première longueur (LV(t)) représentative d'une phase de vol de l'aéronef (12), et calcul (104) d'une deuxième longueur (LS(t)) représentative d'une phase au sol de l'aéronef (12) ; - calcul (106) de la position du point d'arrêt (16) sur la base de la première longueur (LV(t)), de la deuxième longueur (LS(t)) et de la position géographique courante caractérisé en ce que le calcul de la première longueur (LV(t)) comprend une phase (112) d'estimation dynamique d'une distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12). 2.-. Procédé selon la revendication 1, caractérisé en ce qu'il comporte une étape (100) de mesure de la pente courante de l'aéronef (12), la distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12) étant calculée en fonction de la pente courante de l'aéronef (12). 3.- Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que, lors de la phase d'estimation (112), la distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12) est estimée en fonction de la vitesse sol courante de l'aéronef, mesurée lors d'une étape de mesure (100). 4.- Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que, lors de la phase d'estimation (112), la distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12) est estimée en fonction du facteur de charge nzarr de l'aéronef (12), estimé durant son arrondi. 5.- Procédé selon les revendications 2 à 4, caractérisé en ce que, lors de la phase d'estimation (112), la distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12) est estimée par l'équation: gamma(t)- 1 dy D2(t)= fl g x xV! i(t)xcos(y + -dy za' - cos(y)) ° 2 où g est l'accélération de la pesanteur, 13 est la pente du terrain (14), nzarr est le facteur de charge de l'aéronef (12) estimé durant l'arrondi, V',(t) est la vitesse sol courante de l'aéronef (12), et gamma(t) est la pente courante de l'aéronef (12). 6.- Procédé selon l'une quelconque des revendications 2 à 4, caractérisé en ce que la phase d'estimation (112) comporte le calcul d'une hauteur estimée (EH) del'aéronef (12) à laquelle démarre l'arrondi, la distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12) étant calculée à partir de la hauteur estimée (EH). 7.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le calcul de la première longueur (LV(t)) comporte une phase (110) d'estimation d'une distance (D1(t)) parcourue lors de l'approche de l'aéronef (12) vers le terrain (14) avant l'arrondi et/ou lors d'une première partie de la phase d'arrondi de l'aéronef (12), la distance (D1(t)) parcourue lors de l'approche de l'aéronef (12) vers le terrain (14) avant l'arrondi et/ou lors d'une première partie de la phase d'arrondi de l'aéronef (12) étant calculée en fonction de la pente courante mesurée de l'aéronef (12), et avantageusement, en fonction de la hauteur courante de l'aéronef (12) par rapport au terrain (14). 8.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (104) de calcul de la deuxième longueur (LS(t)) comporte une phase (114) de calcul d'une distance (D3(t)) représentative du passage de la décélération courante à la décélération correspondant à un freinage établi de l'aéronef (12), en fonction de la décélération courante. 9.- Procédé selon la revendication 8, caractérisé en ce que la distance représentative du passage de la décélération courante à la décélération correspondant au freinage établi est calculée en fonction d'une constante Km, représentative de la dynamique de l'aéronef (12) entre le toucher du train principal et l'établissement du freinage établi. 10.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape de calcul (104) de la deuxième longueur (LS(t)) comporte une phase (116) de calcul d'une distance (D4(t)) de freinage à partir d'un profil (.1x tab, 1 de décélération .-- prédéterminé de l'aéronef (12) sur le terrain (14), en fonction de la vitesse sol courante. 11.- Procédé selon la revendication 10, caractérisé en ce que le profil de , décélération prédéterminé (Jxtab) est une décélération constante. 12.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (106) de calcul du point d'arrêt (16) comporte la détermination de la position géographique du point d'arrêt (16) sur la base de la position géographique courante de l'aéronef (12), de la première longueur (LV(t)), de la deuxième longueur (LS(t)), et d'un relèvement (X(t)) correspondant à la route de l'aéronef (12). 13.- Système (10) de détermination dynamique de la position du point d'arrêt (16) d'un aéronef (12) sur un terrain (14), comprenant : - un ensemble (30) de mesure dynamique de la position géographique courante de l'aéronef (12), et de sa vitesse sol courante ;- un ensemble (52) de calcul d'une première longueur (LV(t)) représentative d'une phase de vol de l'aéronef (12), et un ensemble (54) de calcul d'une deuxième longueur (LS(t)) représentative d'une phase au sol de l'aéronef (12) ; - un ensemble (56) de calcul de la position du point d'arrêt (16) sur la base de la première longueur (LV(t)), de la deuxième longueur (LS(t)) et de la position géographique courante ; caractérisé en ce que l'ensemble de calcul (52) de la première longueur (LV(t)) comprend des moyens (72) d'estimation d'une distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12). 14.- Système (10) selon la revendication 13, caractérisé en ce qu'il comporte un ensemble (64) de mesure de la pente courante de l'aéronef (12), les moyens d'estimation (72) étant propres à calculer la distance (D2(t)) parcourue par l'aéronef (12) lors d'au moins une partie d'une phase d'arrondi de l'aéronef (12) en fonction de la pente courante de l'aéronef (12).