FR 3049076 A1 20170922 FR 1658049 A 20160830 1 PROCEDE DE CONTROLE D'UN AERONEF SUIVEUR PAR RAPPORT A DES TOURBILLONS GENERES PAR UN AERONEF MENEUR DOMAINE TECHNIQUE La présente invention concerne un procédé de contrôle d'un aéronef suiveur permettant d'estimer les coefficients des moments générés sur un aéronef suiveur par des tourbillons générés par un aéronef meneur et de commander le positionnement de l'aéronef suiveur par rapport à ces tourbillons, ainsi qu'un aéronef suiveur pour lequel un tel procédé de contrôle est mis en oeuvre. ETAT DE LA TECHNIQUE ANTERIEURE Un aéronef en vol génère dans son sillage deux tourbillons, avec un tourbillon partant de chacune de ses ailes du fait de la différence de pression entre l'intrados et l'extrados de l'aile, et de la déflexion de l'écoulement d'air vers le bas qui en résulte. Ces tourbillons sont des tourbillons contrarotatifs et sont caractérisés par un champ de vent globalement ascendant à l'extérieur des tourbillons et globalement descendant entre les tourbillons. A partir des ailes, les tourbillons tendent d'abord à se rapprocher l'un de l'autre, puis à maintenir une distance plus ou moins constante entre eux tout en perdant de l'altitude par rapport à l'altitude à laquelle ils ont été générés. Du fait de cette configuration des tourbillons, il est intéressant, pour un aéronef suiveur qui suit l'aéronef meneur générant les tourbillons, de pouvoir profiter des vents ascendants afin de réduire sa consommation de carburant. EXPOSE DE L'INVENTION Un objet de la présente invention est de proposer un procédé de contrôle qui permet d'estimer les coefficients des moments générés sur un aéronef suiveur par des tourbillons générés par un aéronef meneur afin de contrôler la trajectoire de l'aéronef suiveur dans lequel est implanté un tel procédé de contrôle afin de pouvoir bénéficier des vents ascendants des tourbillons générés par l'aéronef meneur. A cet effet, est proposé un procédé de contrôle de la trajectoire d'un aéronef, dit aéronef suiveur, présentant un processeur et un système de navigation, par rapport à la trajectoire d'un aéronef, dit aéronef meneur, situé devant l'aéronef suiveur et générant des tourbillons, le procédé de contrôle comportant : 3049076 2 - une étape de recueil au cours de laquelle le processeur recueille des données relatives à l'aéronef meneur et à l'aéronef suiveur, - une étape d'estimation au cours de laquelle le processeur estime, à partir des données ainsi recueillies, des coefficients des moments latéraux perturbateurs générés 5 par les tourbillons, - une étape de calcul au cours de laquelle le processeur calcule, à partir des coefficients des moments latéraux perturbateurs ainsi estimés, un écart entre une position actuelle de l'aéronef suiveur et une position à atteindre pour l'aéronef suiveur, la position à atteindre correspondant à une position prédéterminée où 10 l'aéronef suiveur utilise au mieux l'effet de tourbillons, - une étape de détermination au cours de laquelle le système de navigation détermine des commandes à appliquer à des gouvernes de l'aéronef suiveur pour annuler l'écart ainsi calculé, - une étape de commande au cours de laquelle le système de navigation 15 commande les gouvernes de l'aéronef suiveur selon les commandes ainsi déterminées, et - un bouclage sur l'étape de recueil. Un tel procédé permet ainsi à l'aéronef suiveur de maintenir sa position en fonction des tourbillons générés par l'aéronef meneur afin de pouvoir bénéficier des 20 vents ascendants tout en mettant à jour régulièrement les données relatives à l'aéronef meneur. Avantageusement, le coefficient de moments de roulis Cig est donné par la formule Clg = Cl - (Claero + Clsurf) où Cl est un terme d'inertie, CI', est un terme aérodynamique, et Clsurf est un 25 terme de surfaces de contrôle. L'invention propose également un aéronef, dit aéronef suiveur, comportant un système de navigation et un processeur, le processeur étant configuré pour : - recueillir des données relatives à un aéronef, dit aéronef meneur, situé devant l'aéronef suiveur et générant des tourbillons, ainsi que des données relatives à 30 l'aéronef suiveur, - estimer, à partir des données ainsi recueillies, des coefficients des moments latéraux perturbateurs générés par les tourbillons, 3049076 3 - calculer, à partir des coefficients des moments latéraux perturbateurs ainsi estimés, un écart entre une position actuelle de l'aéronef suiveur et une position à atteindre pour l'aéronef suiveur, la position à atteindre correspondant à une position où l'aéronef suiveur utilise au mieux l'effet des tourbillons, 5 - transmettre l'écart au système de navigation, et le système de navigation étant configuré pour : - recevoir l'écart ainsi transmis, - déterminer des commandes à appliquer à des gouvernes de l'aéronef suiveur pour annuler l'écart, 10 - commander les gouvernes de l'aéronef suiveur selon les commandes ainsi déterminées. BREVE DESCRIPTION DES DESSINS Les caractéristiques de l'invention mentionnées ci-dessus, ainsi que d'autres, apparaîtront plus clairement à la lecture de la description suivante d'un exemple de 15 réalisation, ladite description étant faite en relation avec les dessins joints, parmi lesquels : la Fig. 1 montre une vue de dessus d'un exemple de navigation entre un aéronef meneur et un aéronef suiveur selon l'invention, la Fig. 2 représente l'architecture d'un système de contrôle d'un aéronef mettant 20 en oeuvre un procédé de contrôle selon l'invention, la Fig. 3 est un organigramme d'un procédé de contrôle selon l'invention, la Fig. 4 est une vue de derrière de l'aéronef suiveur et des tourbillons générés par l'aéronef meneur, et la Fig. 5 montre la relation entre le coefficient de moment de roulis et les 25 distances latérales. EXPOSE DETAILLE DE MODES DE REALISATION La Fig. 1 montre en vue de dessus un aéronef meneur 10 avec deux ailes 12a-b qui est en vol et qui génère au niveau de chacune des ailes 12a-b, un tourbillon 14a-b (respectivement bâbord - tribord) représenté par son centre. 30 La Fig. 1 montre également un aéronef suiveur 20 selon l'invention qui suit l'aéronef meneur 10 et qui comporte un système de contrôle 50 permettant la mise en oeuvre d'un procédé de contrôle qui est décrit ci-après. 3049076 4 Les positions de l'aéronef meneur 10 et de l'aéronef suiveur 20 sont, de préférence, déterminées dans un repère géographique (latitude, longitude, altitude) et relevées par un GPS du système de navigation de l'aéronef correspondant, mais un autre repère peut également être utilisé et il faut alors mettre en oeuvre des 5 changements de repère. Localement et par convention, on appelle X l'axe longitudinal du fuselage de l'aéronef 10, 20 orienté positivement dans le sens d'avancement de l'aéronef et passant par le centre de gravité de l'aéronef 10, 20, on appelle Y l'axe horizontal qui passe par le centre de gravité de l'aéronef 10, 20 et qui est horizontal, et Z l'axe vertical qui 10 passe par le centre de gravité de l'aéronef 10, 20 et qui est vertical, ces trois directions X, Y et Z étant orthogonales entre elles et formant un repère orthonormé ayant pour origine le centre de gravité O. La Fig. 1 montre également la trajectoire géographique 30 de l'aéronef meneur 10. 15 En référence avec la Fig. 2, le système de contrôle 50 comprend, reliés par un bus de communication 51: un processeur 52 ou CPU (« Central Processing Unit » en anglais), une mémoire vive RAM 53 (« Random Access Memory » en anglais), une mémoire morte ROM 54 (« Read Only Memory » en anglais), une unité de stockage 55 telle qu'un disque dur ou un lecteur de support de stockage, au moins une interface 20 de communication 56, permettant au système de contrôle 50 de communiquer avec des modules extérieurs, tels que des modules de commande des gouvernes de l'aéronef suiveur 20, un système de navigation, une base de données, un récepteur de données, des capteurs. Le processeur 52 est capable d'exécuter des instructions chargées dans la RAM 25 53 à partir de la ROM 54, de l'unité de stockage 55, d'une mémoire externe (non représentée), d'un support de stockage (tel qu'une carte SD), ou d'un réseau de communication. Lorsque l'équipement est mis sous tension, le processeur 52 est capable de lire de la RAM 53 des instructions et de les exécuter. Ces instructions forment un programme d'ordinateur causant la mise en oeuvre, par le processeur 52, 30 de tout ou partie des algorithmes et étapes du procédé de contrôle décrits ci-après. Tout ou partie des algorithmes et étapes du procédé de contrôle décrits ci-après peuvent être implémentés sous forme logicielle par exécution d'un ensemble d'instructions par une machine programmable, par exemple un DSP (« Digital Signal Processor » en anglais) ou un microcontrôleur, ou être implémenté sous forme 3049076 5 matérielle par une machine ou un composant dédié, par exemple un FPGA (« Field-Programmable Gate Array » en anglais) ou un ASIC (« Application-Specific Integrated Circuit » en anglais). La Fig. 3 montre les étapes d'un procédé de contrôle 500 de la trajectoire de 5 l'aéronef suiveur 20 par rapport à la trajectoire de l'aéronef meneur 10 et des tourbillons 14a-b générés par l'avion meneur 10 et qui comporte : - une étape de recueil 502 au cours de laquelle le processeur 52 recueille des données relatives à l'aéronef suiveur 20 et à l'aéronef meneur 10, - une étape d'estimation 504 au cours de laquelle le processeur 52 estime, à 10 partir des données ainsi recueillies, des coefficients des moments latéraux perturbateurs générés par les tourbillons 14a-b, - une étape de calcul 506 au cours de laquelle le processeur 52 calcule, à partir des coefficients des moments latéraux perturbateurs ainsi estimés, un écart entre une position actuelle de l'aéronef suiveur 20 et une position à atteindre pour l'aéronef 15 suiveur 20, la position à atteindre correspondant à une position prédéterminée que l'aéronef suiveur 20 doit occuper pour utiliser au mieux l'effet des tourbillons 14a-b, cette position à atteindre étant également appelée position désirée, - une étape de détermination 508 au cours de laquelle un système de navigation de l'aéronef suiveur 20 détermine les commandes à appliquer aux gouvernes de 20 l'aéronef suiveur 20 pour annuler l'écart ainsi calculé et ainsi approcher l'aéronef suiveur 20 de la position désirée, et - une étape de commande 510 au cours de laquelle le système de navigation commande les gouvernes de l'aéronef suiveur 20 selon les commandes ainsi déterminées. 25 Le but est d'utiliser au mieux le tourbillon 14a-b, ce qui correspond globalement à une conciliation entre une réduction de la consommation de carburant visée et les impacts non-désirés de la turbulence de sillage sur l'aéronef suiveur 20 et pouvant créer un inconfort des passagers. Après l'étape de commande 510, le procédé de contrôle 500 boucle (512) sur 30 l'étape de recueil 502 afin que les calculs postérieurs s'appuient sur des données actualisées et pour garantir que l'aéronef suiveur 20 suive toujours une trajectoire optimisée. L'étape de recueil 502 consiste à recueillir les données nécessaires aux calculs qui doivent être effectués au cours des étapes postérieures. Ces données peuvent être 3049076 6 disponibles dans une base de données de l'aéronef suiveur 20 ou mesurées par des capteurs présents sur l'aéronef suiveur 20. Les données relatives à l'aéronef meneur 10 peuvent être issues d'une base de données de l'aéronef suiveur 20, d'une transmission de données reçues de l'aéronef meneur 10 ou de tout autre moyen 5 approprié. Ces données sont : lxx : le moment d'inertie autour de l'axe X de l'aéronef suiveur 20, Izz : le moment d'inertie autour de l'axe Z de l'aéronef suiveur 20, Ixz : le produit d'inertie autour des axes X, Z de l'aéronef suiveur 20, 10 p, r : les accélérations de roulis et de lacet de l'aéronef suiveur 20, q : la pression dynamique qui s'exerce sur la voilure de l'aéronef suiveur 20, S : la surface aérodynamique de référence de l'aéronef suiveur 20, ici la surface de la voilure de l'aéronef suiveur 20, b : l'envergure de l'aéronef suiveur 20, 15 f3 : l'angle de dérapage de l'aéronef suiveur 20, p : la vitesse de roulis de l'aéronef suiveur 20, VF : la vitesse vraie de l'aéronef suiveur 20, r : la vitesse de lacet de l'aéronef suiveur 20, C113 : le gradient aérodynamique de roulis dû au dérapage de l'aéronef suiveur 20 20, Clp : le gradient aérodynamique de roulis dû à la vitesse de roulis de l'aéronef suiveur 20, : le gradient aérodynamique de roulis dû à la vitesse de lacet de l'aéronef suiveur 20, 25 CnR: le gradient aérodynamique de lacet dû au dérapage de l'aéronef suiveur 20, C,ip : le gradient aérodynamique de lacet dû à la vitesse de roulis de l'aéronef suiveur 20, C, : le gradient aérodynamique de lacet dû à la vitesse de lacet de l'aéronef suiveur 20, 30 Nail : le nombre d'ailerons de l'aéronef suiveur 20, Nsp : le nombre de spoilers de l'aéronef suiveur 20, Nrud : le nombre de gouvernes de direction de l'aéronef suiveur 20, : l'efficacité aérodynamique en roulis due à la déflection de l'aileron i de l'aéronef suiveur 20, 3049076 7 C IÔsp,1 : l'efficacité aérodynamique en roulis due à la déflection du spoiler i de l'aéronef suiveur 20, Cisr,1 : l'efficacité aérodynamique en roulis due à la déflection de la gouverne de direction i de l'aéronef suiveur 20, 5 CnÔa,1 : l'efficacité aérodynamique en lacet due à la déflection de l'aileron i de l'aéronef suiveur 20, C'Ôsp,1 : l'efficacité aérodynamique en lacet due à la déflection du spoiler i de l'aéronef suiveur 20, Cnsr,1 : l'efficacité aérodynamique en lacet due à la déflection de la gouverne de 10 direction i de l'aéronef suiveur 20, bai : la déflection de l'aileron numéro i de l'aéronef suiveur 20, 8spi : la déflection du spoiler numéro i de l'aéronef suiveur 20, Sri : la déflection de la gouverne numéro i de l'aéronef suiveur 20, bL : l'envergure de l'aéronef meneur 10, 15 mL : la masse de l'aéronef meneur 10, VL : la vitesse vraie de l'aéronef meneur 10, p : la densité de l'air, g : la pesanteur, : la corde aérodynamique moyenne de l'aéronef suiveur 20, 20 c(1) : la corde en un point localisé au niveau de la corde aérodynamique moyenne et situé à une distance 1 de l'axe X de l'aéronef suiveur 20, 1 : la distance entre un point sur la voilure de l'aéronef suiveur 20 au niveau de la corde aérodynamique moyenne et le milieu de la voilure au même niveau, dl : l'élément de longueur de la voilure AZ de l'aéronef suiveur 20, 25 Cza : le gradient de coefficient de portance par rapport à l'incidence de l'aéronef suiveur 20, et ai : l'incidence locale sur la voilure de l'aéronef suiveur 20. L'étape d'estimation 504 consiste à utiliser des équations analytiques et des modèles tabulés ou venant de bases de données pour estimer la turbulence. Pour ce 30 faire, un modèle de référence est utilisé. Ce modèle de référence comporte : - des termes d'inertie C1 et Cn qui peuvent être calculés par exemple par les formules : C1= qSb 1'p -1't et 3049076 8 Izzt IxzP cri= qSb - des termes aérodynamiques Ciaero et Cnaero qui peuvent être calculés par exemple par les formules : Claero pb rb = 13 + Cl - + Cl P 2VF 2VF et pb rb Cnaero = Cig 13 + Cn - + Cn ZEF, 5 P 2VF r 2VF, 5 - des termes de surfaces de contrôle Clsurf et Cnsurf qui peuvent être calculés par exemple par les formules : Nad Nsp Nrud C1SUrfClda,iacti +C1,55p,i asPC1drar i=i i=i i=i et Nazi Nsp Nrud Cnsurf = Cnôa,i ai Cndsp,i sP C nsr,i Sri i=1 i=1 i=1 Le modèle de référence permet alors de déterminer les coefficients des moments latéraux perturbateurs de roulis Cig et de lacet Cng qui sont calculés en retranchant la 10 contribution des termes aérodynamiques et de surfaces de contrôle aux termes d'inertie, ce qui correspond aux formules : C1 = C1 (C laero C lsurf) et Cn9 = Cn - (Cnaero Cnsurf) La relation entre le moment de roulis et le coefficient de moment de roulis est par exemple : 15 M 1 = -2. p - s -v2 .1. 1 ou A 11= q-S-l-Ci avec MI: Moment de roulis (N.m), p: la densité de l'air (kg/m3), V : la vitesse de l'air (m/s), S : la surface de référence (m2), 20 1: le bras de levier (m), CI: le coefficient de moment de roulis (SU), q: la pression dynamique (Pa). 3049076 9 L'étape de calcul 506 consiste à déterminer l'écart entre la position actuelle de l'aéronef suiveur 20 et la position à atteindre pour l'aéronef suiveur 20 à partir des coefficients des moments latéraux perturbateurs estimés précédemment. La Fig. 4 montre l'aéronef suiveur 20 et les tourbillons 14a-b de l'aéronef 5 meneur 10. La variable AY est la distance latérale et la variable AZ est la distance verticale et elles définissent l'écart entre la position actuelle de l'aéronef suiveur 20 et la position actuelle du centre du tourbillon 14b le plus proche. Comme cela est expliqué ci-après, on cherche d'abord à déterminer AY et AZ, la 10 position à atteindre, et enfin l'écart entre la position actuelle et la position à atteindre. Dans la mesure où le coefficient de moment de roulis Cig présente une plus grande amplitude que le coefficient de moment de lacet Cng au voisinage d'un tourbillon 14a-b, on utilise une relation faisant intervenir le coefficient de moment de roulis Cig. 15 Le coefficient de moment de roulis Cig s'obtient en intégrant le moment latéral local généré par la portance locale le long de la voilure et la portance locale est liée à l'incidence locale de l'écoulement aérodynamique local sur la voilure. Ceci se traduit par l'équation : b/2 - S - E - CI = [d - c(1) - dl - Czaa, - 1] g -b/2 En aérodynamique des tourbillons de sillage, le champ de vent, globalement 20 ascendant à l'extérieur des tourbillons en fonction de la circulation Fo de l'air dans les tourbillons, est donné par la formule : ro wzo 2n-b, où by est la distance séparant les deux tourbillons 14a-b générés par l'aéronef meneur 10 et qui est donnée par la formule : n- by = -4 LIL La circulation Fo de l'air dans les tourbillons de sillage est calculée à partir de 25 données de l'avion meneur 10 selon la formule : mL. .9 p.17L.bL L'effet d'un tourbillon à une distance 'R' du centre du tourbillon est donné par exemple par le modèle Burnham-Hallock qui constitue un modèle de décroissance spatiale de la circulation et qui est donné par la formule : R2 r(R) ro R2 + R2 ro 3049076 10 où 12, représente le rayon de coeur du tourbillon et fait partie des paramètres du modèle déterminés de manière empirique et prenant typiquement la valeur Rc = 0.04 x LIL La relation liant la vitesse ascendante du tourbillon et la circulation est donnée par la formule : r(R) ro R2 wz(R) = - = 2itR 2itR R2 + l'? 5 En faisant une approximation aux petits angles, le champ spatial des incidences locales est défini par la formule : wz (R) ro R ai(R) = VF = 27FVF R2 + l'? Le coefficient de moment de roulis a l'expression spatiale suivante : CIro ça fb/2 [ R ] = c(l) l dl g 2n-VF S - E Lb/2 R2 + l'? La Fig. 4 montre que la distance 'R' dépend de la distance latérale AY et de la distance verticale AZ de l'avion suiveur 20 par rapport au tourbillon ici 14b, de son 10 attitude latérale cp et de la position '1' du point considéré sur la voilure. L'intégration sur toute la voilure permet donc d'obtenir la relation désirée entre le coefficient de moment de roulis généré par le tourbillon 14b et la position de l'avion suiveur 20 par rapport audit tourbillon 14b. La relation est ainsi donnée par la formule Cig = f (3,17, ,AZ, cp). 15 La détermination de AY et AZ s'effectue à partir de la mesure du coefficient de moment de roulis Cig. Cette détermination peut s'effectuer de différentes manières, soit en adoptant une résolution par valeurs absolues, soit en adoptant une résolution par valeurs relatives. La Fig. 5 montre la relation entre le coefficient de moment de roulis et les 20 distances latérales AY et AZ. La Fig. 5 montre les deux tourbillons 14a-b en vue de derrière. La courbe 302 représente la variation du coefficient de moment de roulis en fonction de la distance latérale AY et la courbe 304 représente la variation du coefficient de moment de roulis en fonction de la distance verticale AZ. Le point A et le point B sont les points de portance virtuelle maximale et le 25 point C est un point où le coefficient de moment de roulis est nul. La résolution par valeurs absolues consiste à affecter une valeur de distance latérale AY à une valeur de coefficient de moment de roulis. A cette fin, l'aéronef suiveur 20 est placé dans une zone dite « bijective » où la relation Cl, = f (AY, 3,z) est bijective et où l'on considère que (p=0 car on fait une convergence latérale quasi- 3049076 11 statique. Ce placement est obtenu par exemple par un modèle de transport de tourbillons. La zone bijective 306 est la zone hachurée qui se situe au-dessus et à l'extérieur du tourbillon 14b considéré et la limite de la zone bijective 306 se situe lorsque le coefficient de moment de roulis passe par un minimum, ici le point A. 5 Dans un exemple, le modèle de transport de tourbillons consiste à calculer la vitesse de descente W, des centres des tourbillons générés par l'aéronef meneur 10, par exemple à partir d'une relation du type : mL.g .n, w = , où : 271-.p.VL.b. mL est la masse de l'aéronef meneur 10 (en kg), 10 g est l'accélération de la pesanteur à l'altitude considérée (en m/s2), p est la densité de l'air au point de vol de l'aéronef meneur 10 (en kg.m3), VL est la vitesse air vraie de l'aéronef meneur 10 (en m.s-1), b, est l'espacement entre les centres des deux tourbillons (en m), et nz est le facteur de charge vertical subi par l'aéronef meneur 10 (en g). 15 L'espacement b, (en mètres) entre les centres des deux tourbillons peut être calculé par exemple par la relation : 7-t- b =-*b v 4 L où : bL est l'envergure de l'aéronef meneur 10 (en mètres). Dans d'autres modes de réalisation, b, peut être envoyé par l'aéronef meneur 10 20 à l'aéronef suiveur 20 ou autres. Les positions des centres des tourbillons au cours du temps sont déterminées des positions de l'aéronef meneur 10 et de la vitesse de descente des tourbillons. Dans la zone bijective 306, la connaissance du coefficient de moment de roulis Cig permet la détermination sans ambiguïté de la distance latérale AY et de la distance 25 verticale AZ de l'aéronef suiveur 20 par rapport au tourbillon 14b le plus proche. La résolution par valeurs relatives consiste à rechercher une distance verticale AZ optimale pour une distance latérale AY fixée, puis à rechercher une distance latérale AY avec une distance verticale AZ fixée à la distance verticale AZ optimale. La recherche de la distance verticale AZ optimale consiste à placer l'aéronef 30 suiveur 20, représenté ici par une croix, loin du tourbillon 14b tout en restant dans sa zone d'influence, et à exploiter la propriété de symétrie du coefficient de moment de roulis par rapport au plan défini par AZ = AZ optimal qui correspond au point de portance virtuelle maximale 'A'. 3049076 12 Ainsi, pour une distance latérale AY fixée et en partant par exemple de la zone supérieure du tourbillon 14b, une première valeur du coefficient de moment de roulis Cigi est relevée, ainsi que l'altitude h1 correspondante. Ensuite l'aéronef suiveur 20 descend en altitude et le coefficient de moment de 5 roulis Cig est relevé au cours de cette descente. Au cours de cette descente, le coefficient de moment de roulis Cig augmente puis diminue jusqu'à passer par la même valeur Cigi pour une autre altitude h2 qui est relevée. Du fait de la symétrie, la distance verticale AZ optimale est calculée comme le milieu du segment [h1, h2] et est donnée par la formule : h1 + h2 àZ optimale = 2 10 La recherche de la distance latérale AY avec une distance verticale AZ fixée à la distance verticale AZ optimale se base sur un coefficient de moment de roulis lui-même basé sur un modèle et qui s'écrit sous la forme de la formule : = k X (AI') où cr (AY) est le modèle dépendant de la distance latérale AY et 'k' une constante et est issu du calcul avec crd(AY) = f(AY,AZ=AZoptimale,(p=0). 15 La formule présente donc deux inconnues 'k' et AY et il est alors possible en utilisant un estimateur approprié d'estimer les valeurs de ces deux inconnues en corrélant les variations du coefficient de moment de roulis Cig avec celles de AY. Pour ce faire, il faut une mesure de la position de l'aéronef suiveur 20, fournie par exemple par un GPS embarqué et une information de déplacement du tourbillon 20 14b, fournie par exemple par le modèle de transport de tourbillons décrit ci-dessus. L'introduction de ces éléments dans un estimateur, tel qu'un estimateur de Kalman par exemple, permet d'estimer AY. Les équations d'état et de mesure utilisées par un tel estimateur sont données par les formules : 25 Le point O est considéré comme un point d'origine O correspondant à un point éloigné du tourbillon 14b où le coefficient de moment de roulis est nul. Selon l'axe latéral, la distance entre le point d'origine O et l'aéronef suiveur 20 est notée Y et la distance latérale AY correspond à la distance entre l'aéronef suiveur 20 et le centre du tourbillon 14b. {[,:1 = [ 01 [,1 io] LYyl Cig = k x (AY) 3049076 13 De la même manière, la distance entre le point d'origine O et le point de portance virtuelle maximale 'A' est Yo et la distance entre le point 'A' et le centre du tourbillon 14b est 'd'. La distance latérale 4Y peut alors s'exprimer par la formule : A Y = Yo - Y + d. 5 La distance 'd' est connue et ne dépend que de l'envergure de l'aéronef suiveur 20. Cette distance est donc facilement modélisable ou peut être issue d'une base de données. La distance 'Y' peut être fournie par exemple par un GPS. La vitesse latérale Y de l'aéronef suiveur 20 par rapport au point d'origine O 10 peut également être fournie par un système de positionnement comme un GPS. La vitesse 3:17 est fournie par un modèle de transport de tourbillons. L'estimateur de Kalman permet alors l'estimation de 'k' et de 4Y. Ainsi, connaissant la distance latérale 4Y et la signature du vortex sur le moment de roulis, la position optimale 'A' de l'aéronef suiveur 20 peut être 15 déterminée. A partir de la connaissance de la position à atteindre 'A' et de la position actuelle de l'aéronef suiveur 20 par rapport au centre du tourbillon 14b, l'écart entre ces deux positions peut être déterminé pour alimenter le système de navigation au cours de l'étape de détermination 508. 20 Ainsi, le processeur 52 calcule l'écart et les lois de pilotage et les mouvements de l'aéronef suiveur 20 ainsi commandés tendent à faire converger cet écart vers 0 avec le temps, afin que la position de l'aéronef suiveur 20 coïncide avec la position désirée. L'aéronef suiveur 20 comporte ainsi un processeur 52 configuré pour mettre en 25 oeuvre l'étape de recueil 502, l'étape d'estimation 504 et l'étape de calcul 506 ainsi que pour transmettre l'écart au système de navigation qui est configuré pour recevoir l'écart ainsi transmis, et mettre en oeuvre l'étape de détermination 508 et l'étape de commande 510. Le processeur 52 est ainsi configuré pour : 30 - recueillir des données relatives à un aéronef meneur 10 ainsi que des données relatives à l'aéronef suiveur 20, - estimer, à partir des données ainsi recueillies, des coefficients des moments latéraux perturbateurs générés par les tourbillons 14a-b, 3049076 14 - calculer, à partir des coefficients des moments latéraux perturbateurs ainsi estimés, un écart entre la position actuelle de l'aéronef suiveur 20 et une position à atteindre pour l'aéronef suiveur 20, - transmettre l'écart au système de navigation, et 5 le système de navigation est configuré pour : - recevoir l'écart ainsi transmis, - déterminer des commandes à appliquer à des gouvernes de l'aéronef suiveur 20 pour annuler l'écart, - commander les gouvernes de l'aéronef suiveur 20 selon les commandes ainsi 10 déterminées. Le système de navigation et le processeur peuvent être deux éléments distincts ou être combinés en un seul composant. Les étapes qui consistent à déterminer 508 les commandes et à commander 510 les gouvernes sont des étapes qui sont classiquement mises en oeuvre dans un aéronef. 15 Par exemple, lorsqu'une destination à atteindre est entrée dans son système de navigation, le système de navigation calcule la différence entre la destination à atteindre et la position actuelle et la trajectoire à suivre pour atteindre cette destination et donc annuler la différence, ainsi que les commandes à appliquer aux gouvernes pour diriger l'aéronef. 20 L'étape de bouclage 512 permet de mettre à jour régulièrement les différentes données relatives à l'aéronef meneur 10 et à l'aéronef suiveur 20 et donc les différents calculs afin que le système de navigation maintienne l'aéronef suiveur 20 en position pendant que l'aéronef meneur 10 et l'aéronef suiveur 20 se déplacent dans des conditions atmosphériques variables. 25 REVENDICATIONS1) Procédé de contrôle (500) de la trajectoire d'un aéronef, dit aéronef suiveur (20) présentant un processeur (52) et un système de navigation, par rapport à la trajectoire d'un aéronef, dit aéronef meneur (10) situé devant l'aéronef suiveur et générant des tourbillons, le procédé de contrôle (500) comportant : - une étape de recueil (502) au cours de laquelle le processeur (52) recueille des données relatives à l'aéronef meneur (10) et à l'aéronef suiveur (20), - une étape d'estimation (504) au cours de laquelle le processeur (52) estime, à partir des données ainsi recueillies, des coefficients des moments latéraux perturbateurs générés par les tourbillons (14a-b), - une étape de calcul (506) au cours de laquelle le processeur (52) calcule, à partir des coefficients des moments latéraux perturbateurs ainsi estimés, un écart entre une position actuelle de l'aéronef suiveur (20) et une position à atteindre pour l'aéronef suiveur (20), la position à atteindre correspondant à une position prédéterminée où l'aéronef suiveur (20) utilise au mieux l'effet de tourbillons (14a-b), - une étape de détermination (508) au cours de laquelle le système de navigation détermine des commandes à appliquer à des gouvernes de l'aéronef suiveur (20) pour annuler l'écart ainsi calculé, - une étape de commande (510) au cours de laquelle le système de navigation 20 commande les gouvernes de l'aéronef suiveur (20) selon les commandes ainsi déterminées, et - un bouclage (512) sur l'étape de recueil (502). 2) Procédé de contrôle (500) selon la revendication 1, caractérisé en ce que le coefficient de moments de roulis Cig est donné par la formule C19 = C1 - (c,,,', + c,s.,'f) 25 où C1 est un terme d'inertie, Cia',, est un terme aérodynamique, et Clsurf est un terme de surfaces de contrôle. 3) Aéronef, dit aéronef suiveur (20), comportant un système de navigation et un processeur (52), le processeur (52) étant configuré pour : 3049076 16 - recueillir des données relatives à un aéronef, dit aéronef meneur (10), situé devant l'aéronef suiveur (20) et générant des tourbillons (14a-b), ainsi que des données relatives à l'aéronef suiveur (20), - estimer, à partir des données ainsi recueillies, des coefficients des moments 5 latéraux perturbateurs générés par les tourbillons (14a-b), - calculer, à partir des coefficients des moments latéraux perturbateurs ainsi estimés, un écart entre une position actuelle de l'aéronef suiveur (20) et une position à atteindre pour l'aéronef suiveur (20), la position à atteindre correspondant à une position où l'aéronef suiveur (20) utilise au mieux l'effet des tourbillons (14a-b), 10 - transmettre l'écart au système de navigation, et le système de navigation étant configuré pour : - recevoir l'écart ainsi transmis, - déterminer des commandes à appliquer à des gouvernes de l'aéronef suiveur (20) pour annuler l'écart, 15 - commander les gouvernes de l'aéronef suiveur (20) selon les commandes ainsi déterminées.