FR 3090852 A1 20200626 FR 1873496 A 20181220 Description Titre de l'invention : OPTIMISATIOn D'UN MODELE PARAMETRIQUE DE PERFORMANCES D'UN AERONEF Domaine de l'invention [0001] L'invention concerne le domaine de l'avionique en général. L'invention concerne en particulier des procédés et des systèmes pour optimiser tout ou partie des opérations d'un aéronef. Technique antérieure [0002] Les approches connues de l'état de la technique (e.g. WO 2017042166, ou US9290262) s'appuient généralement sur un jeu de données modélisant les performances dudit aéronef. Ces données de performances, bien qu'étant connues des fabricants d'aéronefs, ne sont généralement pas publiques, ce qui représente un frein considérable au développement des méthodes d'optimisation des opérations aériennes. [0003] Par ailleurs, lorsqu'elles sont divulguées, ces données décrivent le comportement d'un type d'appareil (de manière générique) et non d'un aéronef en particulier (i.e. pas de données spécifiques). Enfin, les modèles avions peuvent ne pas être libres de droits. [0004] Le cas échéant, les producteurs de systèmes d'optimisation d'opérations aériennes doivent utiliser les données de performances fournies par les fabricants d'aéronefs (le cas échéant) et se retrouvent donc dans une situation de dépendance qui leur est désavantageuse. [0005] Dans l'état de la technique actuel, certains travaux comme « BADA » (acronyme de « Base of Aircraft Data » en anglais) ou «safety-line » par EUCASS (« European Conference for Aeronautics and Space Sciences ») présentent des limitations. Le modèle BADA est limité en termes de poussée et de traînée. Le modèle EUCASS ne s'applique qu'aux avions propulsés par turboréacteurs (pilotés en N1). [0006] Les optimisations actuelles utilisent différentes techniques. En particulier, certaines approches privilégient des méthodes d'estimation paramétriques non linéaires (par exemple par moindres carrés). Ces techniques reposent sur des algorithmes de descente de gradient qui présentent l'inconvénient de converger, parfois, vers des minimas locaux, ce qui conduit à une mauvaise estimation des paramètres du modèle. [0007] Les problèmes techniques en aéronautique font généralement intervenir beaucoup de paramètres différents et les méthodes paramétriques non-linéaires convergent alors peu_ mal ou pas du tout. Cette absence de convergence (ou de convergence vers des minimas locaux) s'observe en effet fréquemment lorsque les processus modélisés sont de grandes dimensions (et suivent des modèles différents en fonction des points d'observation). 2 100081 Les méthodes fondées sur des modèles intégrant des équations physiques de l'aéronef sont généralement dépendantes de la qualité de la modélisation et de la connaissance du comportement réel de l'aéronef. De fait, ces méthodes ne sont généralement pas robustes à la variabilité du comportement réel de chaque aéronef, comparé à celui d'un aéronef « moyen » (modélisé). [0009] TI existe un besoin industriel pour des procédés et des systèmes avancés pour l'optimisation de tout ou partie des opérations d'un aéronef. Résumé de l'invention [0010] Il est divulgué des procédés et des systèmes pour l'optimisation du vol d'un aéronef. La trajectoire est partitionnée en segments, chacun des segments étant régi par des jeux d'équations distincts, qui sont fonctions du mode de poussée du moteur et du guidage vertical (montée, croisière ou descente). En supposant deux modèles, aérodynamique et de régime moteur, des données issues d'enregistrements du vol sont reçues et un nombre de paramètres d'un moteur d'optimisation paramétrique est itérativement déterminé par application d'un calcul des moindres carrés jusqu'à satisfaction d'un critère de minimalité prédéfini. Lc moteur d'optimisation paramétrique sert ensuite à prédire le point de trajectoire suivant un point donné. Des aspects de logiciel et de système (c.g. FMS et/ou EFB) sont décrits. [0011] Avantageusement, le procédé selon l'invention s'appuie sur un modèle de poussée (et de traînée) plus proche de la réalité que celui mis en oeuvre dans le modèle BADA. [0012] Avantageusement, l'invention permet de déterminer un modèle de performances d'un aéronef, et ce indépendamment de son constructeur. [0013] Avantageusement, l'invention permet d'assurer la convergence d'une méthode d'estimation par moindres carrés sur un problème multi-modèles hybride de grande dimension. [0014] L'invention pourra avantageusement être mise en oeuvre dans différents métiers de l'avionique, et en particulier dans les systèmes destinés à l'assistance au pilotage d'un aéronef [0015] En particulier, l'invention peut être déployée dans un calculateur de type FMS ou dans un ensemble de système interconnectant le FMS avec un EFB. Plus spécifiquement, les applications potentielles concernent le calcul de trajectoires, l'assistance à un avionneur pour l'établissement des perfonnances avions, l'optimisation des opérations aériennes d'une compagnie, la simulation de vol ou encore la mise au point de systèmes avioniques au sens large. Description des figures 100161 Différents aspects et avantages de l'invention vont apparaitre en appui de la des- cription d'un mode préféré d'implémentation de l'invention mais non limitatif, avec 3 référence aux figures ci-dessous : [0017] La [fig.1] illustre des exemples d'étapes du procédé selon un mode de réalisation de l' invention ; [0018] La [fig.2] illustre un exemple d'optimisation paramétrique ; Description détaillée de l'invention [0019] La [fig.11 illustre l'influence du mode de poussée et du guidage vertical sur les pa- ramètres SEP, FF et Ni. [0020] L'acronyme SEP 111 pour Specific Excess Power désigne l'énergie disponible pour la montée de l'aéronef, i.e. la capacité de monter de l'avion divisé par la masse (ce paramètre n'est pas constant) [0021] L'acronyme FF pour Fuel Flow désigne la variation de la masse de carburant. [0022] L'acronyme NI fait référence à la vitesse de rotation du premier étage du moteur, celui qui influe le plus sur la consommation de carburant. La puissance disponible est très liée à cette vitesse NI. [0023] Les paramètres SEP, FF et NI sont étroitement liés. [0024] Une façon de quantifier ou de qualifier les relations entre SEP, FF et NI consiste à formuler des modèles (des équations) modélisant ces relations. [0025] Les relations entre les paramètres FF, NI et SEP sont multiples niais dépendent es- sentiellement du mode de pilotage de l'avion, qui se décompose en « mode de poussée » (THR pour « Thrust » III) d'une part (soit adaptée 1111 soit fixe 1111) et en « guidage vertical » 112 d'autre part (c.g montée, croisière, ou descente). [0026] Ces modes de poussée et de guidage vertical définissent N modes de pilotage du moteur utilisés en production (cc sont les modes du pilote automatique c.g. en mode dit Open Accel, IDLE, « energy sharing », moteur au ralenti, avec ou sans sortie des aérofreins, etc). Par exemple, un mode « open acccl » correspond à une poussée à fond pour atteindre et tenir des paramètres cibles prédéfinis. Ces Nt modes correspondent à autant de modèles ou de jeux d'équations. [0027] Chaque mode induit donc une manière de prédire le point de trajectoire suivant (à partir d'un mode [0028] La [fig.2] illustre des exemples d'étapes du procédé selon l'invention. [0029] Dans un mode de réalisation, le procédé comprend l'étape consistant à déterminer un modèle paramétrique comprenant une pluralité de phases de vol décrivant la mécanique du vol de l'avion, et ce à partir d'enregistrements obtenus lors de vols (réels) commerciaux. [0030] Dans un mode de réalisation, le modèle paramétrique est déterminé en utilisant la méthode des moindres carrés. La méthode des moindres carrés permet de comparer des données expérimentales, généralement entachées d'erreurs de mesure, à un modèle ma- thématique censé décrire ces données. Le modèle d'erreur est généralement gaussien (loi normale). Si le modèle d'erreur est non gaussien, le procédé peut recourir à la méthode du maximum de vraisemblance (dont la méthode des moindres carrés est un cas particulier). [0031] Dans un mode de réalisation, le modèle paramétrique est déterminé en utilisant la méthode des moindres carrés récursifs (en anglais « Recursive least squares » acronyme RLS), qui est un filtre adaptatif (minimisant une fonction d'erreur ou les moindres carrés pondérés). [0032] Dans un mode de réalisation, le modèle paramétrique est déterminé en utilisant la méthode des doubles moindres carrés. Cet estimateur permet notamment d'estimer un modèle de régression linéaire à variables instrumentales. [0033] Dans un mode de réalisation, le vol de l'aéronef comprend des phases de vol distinctes, obéissant à des lois différentes selon les segments de vol le long de la trajectoire: les performances de la méthode des moindres carrés sont avantageusement améliorées s'il est tenu compte des informations de poussée et de guidage vertical. [0034] Les phases de vol comprennent notamment les phases de décollage, montée, croisière, changement de niveau, et atterrissage. [0035] Dans un mode de réalisation, le procédé comprend une étape consistant à déterminer une estimation par moindres carrés (ou selon des variantes) d'un modèle de performances avion, cette estimation se faisant spécifiquement en plusieurs « passages » de complexité croissante, chaque passage utilisant 0 un jeu de données enrichi d'une phase de vol par rapport au passage précédent ; et ii) le résultat du passage précédent comme point d'initialisation (la première passe d'initialisation étant faite avec une phase de montée). [0036] L'invention peut utiliser des modèlesparamétriques prédéfinis (i.e. utilisation d'équations). [0037] Plusieurs modèles sont possibles. Dans la description qui suit, deux modèles sont décrits (aérodynamique et moteur), qui sont nécessaires et suffisants pour déterminer un modèle de performances. [0038] Les modèles spécifiés ci-après conjuguent avantageusement réalisme avec rapidité de calcul. [0039] Modèle aérodynamique [0040] Dans un mode de réalisation, l'aérodynamique est modélisée de la façon suivante : [0041] [Math.1] CD CL, M) = E = 1T 15_ oa LiCL [0042] dans lequel, aià représente des coefficients ou des fonctions, C,, (Est) et C,, (drag) éta- blissent le lien entre la portance, l'incidence et la trainée, et M (Mach) désigne le nombre de Mach (qui exprime le rapport de la vitesse locale d'un fluide à la vitesse du son dans ce même fluide). Ce modèle est générique. [0043] Modèle moteur [0044] Dans un mode de réalisation, le régime moteur est modélisé par une fonction affine de type : [0045] [Math.2] TSP-axz +19-xM-Hc dans lequel TSP , a, b et c représentent des coefficients, ip désigne la pression fonction de l'altitude et M désigne le nombre de Mach. 100461 Plus généralement, dans d'autres modes de réalisation, le régime moteur est modélisé par une fonction polynomiale (à coefficients réels ou complexes). 100471 Le motoriste peut disposer d'équations plus complètes ; la (rainée et la poussée peuvent notamment être des fonctions (configuration aérodynamique). 100481 La 1fig.21 illustre un exemple d'optimisation paramétrique selon l'invention, notamment les étapes de raffinement par ajout successif des phases de vol. 100491 A l'étape 210 (qui est préalable, hors-ligne, d'initialisation), le procédé selon l'invention consiste à segmenter chaque vol en un nombre n phase de phases de vol modélisables par un ensemble E d'équations prédéfinies dépendant d'un nombre p_exp de variables explicatives et d'un nombre n par de paramètres et fournissant un vecteur d'observation de taille q mes. [0050] Les variables explicatives prennent soit des valeurs numériques continues, soit des variables discrètes décrivant les différents sous-processus suivis par le processus modélisé. [0051] A l'étape 220, à partir d'un ensemble de fichiers (provenant par exemple d'enregistrements de vols commerciaux, c'est-à-dire de données réelles), un tableau d'observations est déterminé, lequel comprend un nombre In obs de lignes et p_exp + mes colonnes et contenant les points associés à la phase de montée de l'aéronef [0052] A l'étape 230, le procédé selon l'invention comprend une étape consistant à dé- terminer (par exemple de façon arbitraire, mais pas toujours afin d'accélérer les calculs) un nombre npar de paramètres initiaux et à calculer le critère des moindres carrés du tableau à l'étape 220. [0053] Le critère des moindres carrés peut s'obtenir de la façon décrite ci-après. [0054] Pour chaque ligne du tableau, le procédé détermine- à l'aide de l'ensemble E d'équations définies à l'étape 1 les obs composantes (Y1...Yrn obs) du vecteur Y= E(XI,..Xp exp,P1,...Pn_par). Puis le procédé détermine la norme euclidienne N du vecteur [Y- Yobs], les composantes du vecteur Yobs étant les m obs dernières valeurs de la ligne i du tableau. Il est enfin calculé la moyenne de N sur l'ensemble des lignes 6 du tableau. [0055] A l'étape suivante, le résultat ainsi obtenu est fourni à un moteur d'optimisation para- métrique 240, qui en déduira une itérée des n_par paramètres initiaux et répètera les étapes 230 et 240 jusqu'à satisfaction d'un critère de minimal ité 241, en donnant, à la fin de cc processus, un ensemble de n_par paramètres, dits finaux ou optimisés. [0056] A l'étape 250, les étapes 220, 230 et 240 sont répétées selon des variations pré- définies. [0057] Selon une variation, le tableau d'observations de l'étape 220 contient les points de montée et les points de descente de l'aéronef et les paramètres initiaux de l'étape 230 sont les paramètres finaux obtenus à l'étape 240. [0058] Selon une autre variation, les étapes 220, 230 et 240 sont réitérées, dans lesquelles le tableau d'observations à l'étape 220 contient les points de montée, les points de descente et les points de croisière de l'aéronef, les paramètres initiaux de l'étape 230 étant les paramètres finaux obtenus à l'étape 240. [0059] Dans un mode de réalisation, le tableau d'observations comprend des informations relatives au mode de poussée et au guidage vertical. [0060] L'invention peut s'implémenter à partir d'éléments matériel et/ou logiciel. Elle peut être disponible en tant que produit programme d'ordinateur sur un support lisible par ordinateur. [0061] Dans une variante de réalisation, une ou plusieurs étapes du procédé peuvent être mises en oeuvre au sein d'un calculateur au sol. [0062] Avantageusement, si le calculateur se trouve à bord d'un aéronef et est connecté à un enregistreur de paramètres, les données de vol peuvent être utilisées dans une architecture temps réel pour améliorer la connaissance des performances avion en temps réel. [0063] Dans une variante de réalisation, une ou plusieurs étapes du procédé selon l'invention sont mises en oeuvre sous forme de programme d'ordinateur hébergé sur un ordinateur portable de type « EFB » (acronyme d'« Electronic Flight Bag). Dans une variante de réalisation, une ou plusieurs étapes du procédé peuvent être mises en oeuvre au sein d'un calculateur de type « FMS » (acronyme de « Flight Management System ») ou dans une fonction FM d'un calculateur de vol. [0064] Il est décrit un procédé pour l'optimisation du vol d'un aéronef, le vol dudit aéronef étant segmenté en un nombre n_phase de phases de vol modélisables par un ensemble E d'équations prédéfinies dépendant d'un nombre p_exp de variables explicatives et d'un nombre n_par de paramètres et fournissant un vecteur d'observation de taille q_mes, le procédé comprenant les étapes consistant à : recevoir une pluralité d'enregistrements de vols et déterminer un tableau comportant un nombre ,n_obs de lignes et p_exp + q_mes colonnes; déterminer un nombre n_par de paramètres et 7 calculer le critère des moindres carrés du tableau ainsi déterminé ; déterminer n_par paramètres dits finaux, par itération jusqu'à satisfaction d'un critère de minimalité prédéfini dans un moteur d'optimisation paramétrique. [0065] Les variables explicatives peuvent être des valeurs numériques continues et/ou des variables discrètes décrivant les différents sous-processus suivis par le processus modélisé. [0066] Dans un mode de réalisation, l'étape de calcul par moindres carrés comprenant les étapes consistant à : - calculer pour chaque ligne du tableau, à de l'ensemble E d'équations définies, les m obs composantes (Y1 ...Yrn ()lis) du vecteur Y= E(X1,..Xp exp,P1,...Pn par) ; - déterminer la norme euclidienne N du vecteur Y - Yobs, les composantes du vecteur Yohs étant les m obs dernières valeurs de la ligne 1 du tableau et déterminer la moyenne de N sur l'ensemble des lignes du tableau. [0067] Dans un mode de réalisation, le tableau d'observations comprenant des informations relatives au mode de poussée et au guidage vertical. Ces données (e.g. à fond, ralenti, montée, descente, etc) servent aux intégrateurs (prédiction). [0068] Dans un mode de réalisation, l'aérodynamique est modélisée par un modèle reliant la portance, l'incidence et la traînée, et M qui désigne le nombre de Mach. Le nombre de Mach exprime le rapport de la vitesse locale d'un fluide à la vitesse du son dans ce même fluide. [0069] Dans un mode de réalisation, le régime moteur couplé au modèle aérodynamique est modélisé par une fonction polynomiale. [0070] Dans un mode de réalisation, les paramètres initiaux sont les paramètres finaux. [0071] Dans un mode de réalisation, le critère des moindres carrés est un critère des moindres carrés récursifs. [0072] 11 est décrit un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer une ou plusieurs des étapes du procédé, lorsque ledit programme est exécuté sur un ordinateur. [0073] 11 est décrit un système pour la mise en oeuvre une ou plusieurs des étapes du procédé, le système comprenant un ou plusieurs systèmes avioniques. [0074] Dans un mode de réalisation, les calculs sont effectués au sol. [0075] Dans un mode de réalisation, le système comprend en outre un ou plusieurs systèmes non-avioniques de type sacoche de vol électronique EFB. [0076] Dans un mode de réalisation, le système comprend des systèmes embarqués de calcul et de stockage pour le traitement temps réel des performances de l'aéronef. 8 [Revendication 1] [Revendication 2] [Revendication 3] [Revendication 4] [Revendication 5]