FR 2988868 A1 20131004 FR 1252895 A 20120330 L'invention concerne les drones à voilure tournante tels que les quadricoptères et analogues. Ces drones sont pourvus de rotors multiples entraînés par des moteurs respectifs commandables de manière différenciée afin de piloter le drone en attitude et en vitesse. Un exemple typique d'un tel drone est le AR.Drone de Parrot SA, Paris, France, qui est un quadricoptère équipé d'une série de capteurs (accéléromètres, gyromètres trois axes, altimètre), d'une caméra frontale captant une image de la scène vers laquelle est dirigé le drone, et d'une caméra de visée verticale captant une image du terrain survolé. Les WO 2010/061099 A2 et EP 2 364 757 Al (Parrot SA) décrivent un tel drone ainsi que son principe de pilotage par l'intermédiaire d'un téléphone ou baladeur multimédia à écran tactile et accéléromètre intégré, par exemple un téléphone cellulaire de type iPhone ou un baladeur ou une tablette multimédia de type iPod Touch ou iPad (marques déposées de Apple Inc., USA). Plus précisément, dans ce mode de fonctionnement que l'on appellera par la suite "mode contrôlé", le drone est piloté par l'utilisateur au moyen de signaux émis par le détecteur d'inclinaisons de l'appareil, inclinaisons qui seront répliquées par le drone : par exemple, pour faire avancer le drone l'utilisateur incline son appareil selon son axe de tangage, et pour déporter le drone à droite ou à gauche il incline ce même appareil par rapport à son axe de roulis. De cette manière, si le drone est commandé de manière à s'incliner ou "plonger" vers le bas (inclinaison suivant un angle de tangage), il progressera vers l'avant avec une vitesse d'autant plus élevée que l'inclinaison sera importante ; inversement s'il est commandé de manière à se "cabrer" dans le sens opposé, sa vitesse ralentira progressivement puis s'inversera en repartant vers l'arrière. De la même façon, pour une commande d'inclinaison suivant un axe de roulis le drone penchera à droite ou à gauche, provoquant un déplacement linéaire en translation ho- rizontale vers la droite ou vers la gauche. L'utilisateur dispose d'autres commandes affichées sur l'écran tactile, notamment de "montée/descente" (commande des gaz) et de " pivotement à droite/pivotement à gauche" (pivotement du drone autour de son axe de lacet). Le drone est également pourvu d'une commande de passage en point fixe : dans ce mode de fonctionnement que l'on appellera par la suite "mode autopiloté", lorsque l'utilisateur lâche toutes les commandes de son appareil de télécommande, le drone s'immobilise en point fixe et s'y stabi- lise de façon entièrement automatique. Le problème général de l'invention est le pilotage d'un tel drone dans le vent. On appellera "vent" un déplacement d'air parallèlement au sol, constant dans l'espace et dans le temps (les rafales de vent ne seront pas considé- rées). Un tel "vent" est alors complètement déterminé par ses composan- tes de vitesse horizontale dans un repère terrestre donné, ou bien par sa norme et sa direction dans un tel repère. Le vent représente une perturbation qui, par ses caractéristiques, affecte le drone lorsqu'il est en vol. Ces perturbations sont particulièrement sen- sibles dans le cas des "microdrones" très légers, tels que le AR.Drone précité, dont la masse est de quelques centaines de grammes seulement. La prise au vent d'un tel drone est donc particulièrement sensible, même par vent modéré. Le drone s'appuie pour se déplacer sur une masse d'air mobile en fonc- tion du vent, tandis qu'il mesure par rapport au sol ses composantes hori- zontales de vitesse (par exemple au moyen d'une caméra pointant vers le bas et donnant une image du sol), d'angles (par les accéléromètres et gyromètres embarqués) et de forces extérieures qui lui sont appliquées (par les accéléromètres embarqués). En l'absence de vent le drone peut garder une vitesse nulle et donc une position fixe en maintenant une assiette nulle. Il est donc primordial de mesurer avec précision cet état neutre. Or, les accéléromètres possèdent par construction un biais (que l'on supposera constant sur le court terme) qui doit être estimé et compensé pour utiliser correctement la mesure d'accélération dans la reconstruction des angles. Les accéléromètres étant entre autres utilisés en inclinomètres, une erreur d'estimation des biais des accéléromètres entraîne un biais sur la mesure d'angle, néfaste à la tenue du point fixe. De même ce biais de capteur influe sur la mesure des forces du vent réalisé par l'accéléromètre. En présence de vent, le drone doit s'incliner pour compenser le vent et maintenir le point fixe (asservi sur une image stable de la caméra verticale) : son assiette ne sera donc pas nulle et l'évaluation du biais sera modifiée. Ceci n'empêche pas de maintenir le drone en point fixe (asser- vissement par la caméra), mais va avoir une incidence particulièrement sensible dans deux cas - en point fixe lors d'un mouvement de lacet seul (rotation autour d'un axe vertical, sans déplacement avant/arrière ni droite/gauche) : comme le drone était penché dans une direction précise pour contrer le vent, dès qu'il commence à tourner son inclinaison ne correspond plus à la direction de compensation du vent mais va faire un angle avec cette dernière : de ce fait, le drone va "partir dans le vent" et ne tiendra plus son point fixe ; - en pilotage contrôlé, les déplacements dans le sens du vent sont ac- centués, tandis que les déplacements dans le sens opposé du vent sont freinés, voire rendus impossibles. L'un des buts de l'invention est de proposer un procédé d'estimation du vent permettant de distinguer entre biais de construction de l'accéléromètre et biais de mesure induit par le vent. Un autre but de l'invention est de proposer un procédé de compensation du vent permettant de pallier les deux phénomènes que l'on vient de décrire, c'est-à-dire i) en vol autopiloté, permettre au drone de maintenir son point fixe dans le vent lors d'un mouvement de lacet et ii), en mode commandé, de procurer un contrôle du déplacement qui soit identique dans tous les sens, que ce soit dans le sens du vent ou dans le sens opposé au vent, de façon transparente pour le pilote. Le FR 2 926 637 Al décrit un système d'estimation de la vitesse d'un drone par rapport à l'air, par application d'un modèle exprimant les relations mécaniques entre les forces agissant sur le drone et ses vitesses et accélérations en translation et en rotation. La vitesse par rapport à l'air est notamment estimée à partir de la force de traînée aérodynamique s'opposant au déplacement du drone, qui est indirectement mesurée et utilisée pour obtenir une information sur la vitesse par rapport à l'air. Pour mesurer la vitesse du vent, c'est-à-dire la différence entre i) la vi- tesse du drone par rapport à l'air et ii) la vitesse du drone par rapport au sol, ce document propose d'utiliser un récepteur GPS embarqué mesurant la vitesse du drone par rapport au sol. Cette manière de procéder requiert toutefois la présence d'un tel récepteur dans le drone, et nécessite qu'un signal GPS soit disponible et que sa précision soit suffisante (visibilité d'un nombre suffisant de satellites). En tout état de cause, la précision obtenue ne procure pas une estimation de la vitesse du vent suffisamment fine pour évaluer correctement la part due au vent du biais de mesure du capteur accélérométrique, nécessaire pour compenser les inconvénients de pilotage dans le vent exposés plus haut. Pour résoudre ces problèmes, et atteindre les buts exposés plus haut, l'invention propose une autre technique d'estimation de la vitesse du vent et du biais du capteur. L'invention propose également une technique de compensation permettant d'utiliser les résultats de cette estimation (ou d'une estimation opérée d'une autre manière) pour corriger les défauts de pilotage dans le vent exposés plus haut, aussi bien en mode autopiloté (point fixe) qu'en mode contrôlé (commandes de pilotage appliquées par un utilisateur). Plus précisément, l'invention propose un procédé de pilotage d'un drone à voilure tournante à rotors multiples entraînés par des moteurs respectifs commandables pour piloter le drone en attitude et en vitesse comprenant, de manière en elle-même connue, la génération de consignes angulaires et l'application de ces consignes à une boucle de contrôle des moteurs du drone, ces consignes étant aptes à commander l'attitude du drone selon des axes de tangage et de roulis. De façon caractéristique, le procédé de l'invention comprend en outre les étapes suivantes : l'établissement d'au moins un modèle dynamique du drone, décrivant des composantes de vitesse horizontale du drone en fonction des coefficients de traînée et de la masse du drone, des angle d'Euler caractérisant l'attitude du drone par rapport à un repère terrestre absolu, ainsi que de la rotation du drone autour d'un axe vertical ; la me- sure de la force de traînée aérodynamique du drone ; la mesure de la vitesse relative du drone par rapport au sol ; et l'application audit modèle dynamique du drone desdites mesures de force de traînée aérodynamique et de vitesse relative par rapport au sol, de manière à produire une estimation des composantes horizontales de vitesse du vent latéral. La mesure de la force de traînée aérodynamique du drone peut notamment être dérivée d'une mesure d'accélération du drone. L'étape d'application des composantes horizontales de vitesse de vent latéral au modèle dynamique du drone comprend très avantageusement l'utilisation d'un filtre prédictif de Kalman, notamment un filtre à six états, ces états comprenant : deux composantes de la vitesse de déplacement du drone par rapport au sol, exprimées dans un repère lié au drone ; deux composantes de la vitesse du vent par rapport au sol, exprimées dans un repère terrestre absolu ; et deux composantes horizontales de biais d'ac- céléromètre du drone. Avantageusement, l'établissement du modèle dynamique du drone comprend l'établissement de deux modèles distincts, avec : i) un modèle dynamique du drone en vol, et ii) un modèle dynamique du drone au sol, ces modèles étant utilisés sélectivement en fonction de l'état du drone. Avantageusement, le modèle dynamique du drone en vol peut utiliser une mesure d'accélération du drone comme mesure proportionnelle à la vitesse de l'air dans le repère du drone. Le modèle dynamique du drone en vol peut notamment être du type : (00 0 o) 0 0) 0 0) 0 (0 0 0 \O 0 (00 o (oo o. X= X+ C, Cx in in C C , v Y -(Oz - _::.... 0 représentant l'évolution de la valeur initiale Xdu vecteur d'état, l'entrée U du système étant la gravité U = g, G et Cy étant les coefficients de frottement aérodynamique du drone selon ses axes u et y, m étant la masse du drone, ço, 0 et yr étant les angles d'Euler (respectivement de roulis, de tangage et de lacet) caractérisant l'attitude du drone par rapport au repère NED, coz étant le taux de rotation autour de l'axe w (mouvement de rotation en lacet), et R'v étant la matrice de rotation associée à l'angle v, en dimension 2. Le modèle dynamique du drone au sol peut notamment être du type : où T est un paramètre temporel de décroissance progressive de l'estimation de la vitesse du vent. L'invention propose également une étape de compensation de l'effet du vent latéral sur le positionnement et les déplacements du drone, par géné- ration de consignes correctrices, fonction des composantes horizontales estimées de vitesse du vent latéral, et combinaison de ces consignes correctrices aux consignes angulaires de tangage et de roulis appliquées à la boucle de contrôle des moteurs du drone. Les consignes correctrices peuvent notamment être du type : = Cx - COS (liewind 1P drone) l) - d_comp s OPnd - dron Vivind étant le module de la vitesse estimée du vent, Cx étant le coefficient de traînée du drone, tifwind étant l'angle de cap de la direction du vent et Vdt'one étant l'angle de cap de la direction du drone. L'étape de compensation peut notamment comprendre la définition pour le drone d'une assiette de référence en boucle ouverte. Lorsque le drone est en vol contrôlé par un utilisateur, les consignes correctrices sont combinées aux consignes de pilotage appliquées par l'utilisateur. Lorsque le drone est en vol stationnaire autopiloté, ces consignes correctrices sont combinées à des consignes de stabilisation de point fixe produites en ré- ponse à une mesure de la vitesse horizontale du drone par rapport au sol. (0 0 k-0 0 (0 0 0 ',41 0, 0 = (-T 0 0" 0 0 0 On va maintenant décrire un exemple de mise en oeuvre du dispositif de l'invention, en référence aux dessins annexés où les mêmes références numériques désignent d'une figure à l'autre des éléments identiques ou fonctionnellement semblables. La Figure 1 est une vue d'ensemble montrant le drone et l'appareil de télécommande associé permettant son pilotage à distance. La Figure 2 est un schéma par blocs des différents organes de contrôle, d'asservissement et de pilotage assisté du drone. Les Figures 3 et 4 sont des relevés comparatifs, respectivement du module de la vitesse du vent et de la direction du vent, des valeurs simultanément estimées et mesurées dans des conditions réelles, montrant la pertinence du procédé d'estimation selon l'invention. La Figure 5 est une représentation schématique fonctionnelle illustrant les boucles d'estimation et de compensation du vent, pour une configuration de vol du drone en point fixe autopiloté. La Figure 6 est une représentation schématique fonctionnelle illustrant la boucle d'estimation du vent, pour une configuration de vol du drone en vol 20 contrôlé. 0 Sur la Figure 1, la référence 10 désigne de façon générale un drone, par exemple un quadricoptère tel que le AR.Drone de Parrot SA, Paris, 25 France, décrit dans les WO 2010/061099 A2 et EP 2 364 757 Al précités, ainsi que dans le FR 2 915 569 Al (qui décrit notamment le système de contrôle avec gyromètres et accéléromètres utilisé par le drone). Le drone 10 comporte quatre rotors coplanaires 12 dont les moteurs sont pilotés indépendamment par un système intégré de navigation et de 30 contrôle d'attitude. Il est pourvu d'une première caméra 14 à visée frontale permettant d'obtenir une image de la scène vers laquelle est orienté le drone, par exemple une caméra grand angle à capteur CMOS. Le drone comporte également une seconde caméra à visée verticale (non visible sur la Figure 1) pointant vers le bas, apte à capter des images suc- cessives du terrain survolé et utilisée notamment pour évaluer la vitesse du drone par rapport au sol. Des capteurs inertiels (accéléromètres et gyromètres) permettent de mesurer avec une certaine précision les vitesses angulaires et les angles d'attitude du drone, c'est-à-dire les angles d'Euler décrivant l'inclinaison du drone. De façon générale, par "inclinaison", on entendra l'inclinaison du drone par rapport à un plan horizontal d'un repère terrestre fixe, étant entendu que les deux composantes longitudinale et transversale de la vitesse horizontale sont intimement liées à l'inclinaison suivant les deux axes respectifs de tangage et de roulis. Deux repères peuvent être considérés dans la suite, selon le cas : - un repère local {u,v,w} lié au corps du drone : le drone, bien que fortement symétrique par construction, comporte un avant et un arrière, et l'on considérera la position de la caméra comme pointant vers l'avant, définissant ainsi l'axe u. ; l'axe y est perpendiculaire à u dans le plan moyen du drone, et l'axe w est l'axe vertical dirigé vers le sol ; - un repère absolu de type NED (North East Down), qui est un repère terrestre fixe, notamment un repère géographique {XNED,YNED,ZNED} l'axe XNED correspond par exemple à la direction du Nord magnétique, la direction YNED est la direction parallèle au plan du sol perpendiculaire au Nord géographique (c'est-à-dire l'Est géographique), et la direction ZNED est perpendiculaire au plan du sol et orientée vers le bas. Un télémètre à ultrasons disposé sous le drone et un capteur barométrique embarqué fournissent par ailleurs des mesures qui, combinées, don- nent une estimation de l'altitude du drone par rapport au sol. En ce qui concerne la vitesse linéaire de translation du drone dans un plan horizontal par rapport au sol, celle-ci est évaluée par analyse de l'image fournie par la caméra à visée verticale du drone en combinaison avec les données accélérométriques, grâce à un logiciel qui estime le dé- placement de la scène captée par la caméra d'une image à la suivante et applique à ce déplacement estimé un facteur d'échelle fonction de l'altitude mesurée. Cette technique est décrite en détail dans le EP 2 400 460 Al (Parrot SA), auquel on pourra se référer pour de plus amples détails. Le drone 10 est piloté par un appareil de télécommande distant 16 pourvu d'un écran tactile 18 affichant l'image embarquée par la caméra frontale 14, avec en superposition un certain nombre de symboles permettant l'activation de commandes de pilotage par simple contact du doigt 20 d'un utilisateur sur l'écran tactile 18. L'appareil 16 est pourvu de moyens de liaison radio avec le drone pour l'échange bidirectionnel de données du drone 10 vers l'appareil 16, notamment pour la transmission de l'image captée par la caméra 14, et de l'appareil 16 vers le drone 10 pour l'envoi de commandes de pilotage. Cette liaison peut être par exemple de type réseau local Wi-Fi (IEEE 802.11) ou Bluetooth (marques déposées). L'appareil 16 est également pourvu de capteurs d'inclinaison permettant de contrôler l'attitude du drone en imprimant à l'appareil des inclinaisons correspondantes. Comme indiqué en introduction, l'appareil de télécommande 16 est avantageusement constitué par un téléphone ou baladeur multimedia à écran tactile et accéléromètre intégré, par exemple un téléphone cellulaire de type iPhone, un baladeur de type iPod Touch ou une tablette multimédia de type iPad, qui sont des appareils incorporant les divers organes de contrôle nécessaires à l'affichage et à la détection des commandes de pilotage, à la visualisation de l'image captée par la caméra frontale, et à l'échange bidirectionnel de données avec le drone par liaison Wi-Fi ou Bluetooth. Le pilotage du drone 10 consiste à faire évoluer celui-ci en commandant les moteurs de façon différenciée pour engendrer, selon le cas, des mouvements de : a) rotation autour d'un axe de tangage 22, pour faire avancer ou reculer le drone ; et/ou b) rotation autour d'un axe de roulis 24, pour décaler le drone vers la droite ou vers la gauche ; et/ou c) rotation autour d'un axe de lacet 26, pour faire pivoter vers la droite ou vers la gauche l'axe principal du drone ; et/ou d) translation vers le bas 28 ou vers le haut 30 par changement du ré- gime des gaz, de manière à respectivement réduire ou augmenter l'altitude du drone. Lorsque ces commandes de pilotage sont appliquées par l'utilisateur à partir de l'appareil de télécommande 16 (configuration de "pilotage contrô- lé", commandé par un utilisateur), les commandes a) et b) de pivotement autour des axes de tangage 22 et de roulis 24 sont obtenues par des inclinaisons de l'appareil 16 respectivement autour de son axe longitudinal 32 et de son axe transversal 34 : par exemple, pour faire avancer le drone il suffit d'incliner l'appareil vers l'avant en le penchant autour de l'axe 32, pour le déporter à droite il suffit d'incliner l'appareil en le penchant autour de l'axe 34 vers la droite, etc. Les commandes c) et d), quant à elles, résultent d'actions appliquées par contact du doigt 20 de l'utilisateur (généralement le doigt droit) sur des zones spécifiques correspondantes de l'écran tactile 18. On pourra se référer aux WO 2010/061099 A2 et EP 2 364 757 Al préci- tés pour plus de détails sur ces aspects du système. Le drone possède également un système automatique et autonome de stabilisation en vol stationnaire (configuration de "point fixe" autopiloté), activé notamment dès que l'utilisateur retire son doigt de l'écran tactile de l'appareil, ou automatiquement à la fin de la phase de décollage, ou en- core en cas d'interruption de la liaison radio entre l'appareil et le drone. Le drone passe alors à un état de sustentation où il sera immobilisé et stabilisé automatiquement, sans aucune intervention de l'utilisateur. Contrôle et asservissement des moteurs du drone On va maintenant exposer la manière dont sont élaborées les consignes de pilotage des moteurs du drone, en référence à la Figure 2 qui est un schéma par blocs fonctionnels des différents organes de contrôle et d'as- servissement du drone. On notera que, bien que ces schémas soient présentés sous forme de circuits interconnectés, la mise en oeuvre des différentes fonctions est essentiellement logicielle, cette représentation n'étant qu'illustrative. De façon générale, comme illustré Figure 2, le système de pilotage impli- que plusieurs boucles imbriquées pour le contrôle de la vitesse horizon- tale, de la vitesse angulaire et de l'attitude du drone, en plus des variations d'altitude. La boucle la plus centrale, qui est la boucle 100 de contrôle de la vitesse angulaire, utilise d'une part les signaux fournis par des gyromètres 102 et d'autre part une référence constituée par des consignes de vitesse angu- laire 104. Ces différentes informations sont appliquées en entrée d'un étage 106 de correction de la vitesse angulaire, qui pilote lui-même un étage 108 de contrôle des moteurs 110 afin de commander séparément le régime des différents moteurs pour corriger la vitesse angulaire du drone par l'action combinée des rotors entraînés par ces moteurs. La boucle 100 de contrôle de la vitesse angulaire est imbriquée dans une boucle 112 de contrôle d'attitude, qui opère à partir des indications fournies par les gyromètres 102, par des accéléromètres 114 et par un magnétomètre 116 donnant l'orientation absolue du drone dans un repère géomagnétique terrestre. Les données issues de ces différents capteurs sont appliquées à un étage 118 d'estimation d'attitude de type PI (Proportionnel-Intégrateur). L'étage 118 produit une estimation de l'attitude réelle du drone, appliquée à un étage 120 de correction d'attitude qui compare l'attitude réelle à des consignes d'angle générées par un circuit 122 à par- tir de commandes directement appliquées par l'utilisateur 124 (configura- tion de "pilotage contrôlé) et/ou de données générées en interne par le pilote automatique du drone via le circuit 126 de correction de vitesse horizontale pour le maintien du drone en position immobile (configuration de "point fixe" autopiloté). Les consignes éventuellement corrigées appli- quées au circuit 120 et comparées à l'attitude réelle du drone sont trans- mises par le circuit 120 au circuit 104 pour commander les moteurs de manière appropriée. En résumé, à partir de l'erreur entre la consigne (appliquée par l'utilisateur et/ou générée en interne) et la mesure de l'angle (donnée par le circuit d'estimation d'attitude 118), la boucle 112 de contrôle d'attitude calcule une consigne en vitesse angulaire à l'aide du correcteur PI du circuit 120. La boucle 100 de contrôle de la vitesse angulaire calcule ensuite la différence entre la consigne de vitesse angulaire précédente et la vitesse angulaire effectivement mesurée par les gyromètres 102. La boucle calcule à partir de ces informations les différentes consignes de vitesse de rota- tion (et donc de force ascensionnelle), qui sont envoyées aux moteurs 110 pour réaliser la manoeuvre demandée par l'utilisateur et/ou planifiée par le pilote automatique du drone. De façon caractéristique, le circuit 122 de calcul des consignes d'angle reçoit également des corrections délivrées par un circuit 128 d'estimation et de compensation du vent, que l'on décrira plus bas en détail. La boucle de contrôle de vitesse horizontale 130 utilise une caméra vidéo verticale 132 et des capteurs d'altitude 134 (télémètre à ultrasons combi- né à un capteur barométrique). Un circuit 136 assure le traitement des images produites par la caméra verticale 132, en combinaison avec les signaux de l'accéléromètre 114 et du circuit d'estimation d'attitude 118, pour estimer au moyen d'un circuit 138 les composantes de la vitesse ho- rizontale selon les deux axes de tangage et de roulis du drone. Les vites- ses horizontales estimées sont corrigées par l'estimation de vitesse verticale donnée par un circuit 140 et par une estimation de l'altitude donnée par le circuit 142 à partir des informations des capteurs 134. En ce qui concerne les déplacements verticaux du drone, l'utilisateur 124 applique des commandes à un circuit 144 de calcul de consignes d'altitude, consi- gnes qui sont appliquées à un circuit 146 de calcul de consigne de vitesse ascensionnelle Vz via le circuit de correction d'altitude 148 recevant l'estimation d'altitude donnée par le circuit 142. La vitesse ascensionnelle Vz calculée est appliquée à un circuit 150 qui compare cette vitesse de consigne à la vitesse correspondante estimée par le circuit 140 et modifie en conséquence les données de commande des moteurs (circuit 108) de manière à augmenter ou réduire la vitesse de rotation simultanément sur tous les moteurs de façon à minimiser l'écart entre vitesse ascensionnelle de consigne et vitesse ascensionnelle mesurée. Élaboration des consignes de contrôle d'attitude du drone On va maintenant exposer de façon générale (et en l'absence de vent) la manière dont sont élaborées les consignes de contrôle du drone au moyen des circuits que l'on vient de décrire. On désignera par : - u la composante de la vitesse de translation horizontale dans la direction de l'axe principal du drone (le long de l'axe de roulis 24) ; - y la composante de la vitesse de translation horizontale dans la direction transversale le long de l'axe de tangage 22) ; et - w la vitesse de translation verticale, le tout dans le repère {u,v,w} lié au drone, indépendant de l'inclinaison de celui-ci par rapport au repère NED de référence. Chacune des quatre hélices i du drone (i = 1 ... 4) exerce un couple Fi et une force de poussée ascensionnelle F1 proportionnelle au carré de la vi- tesse de rotation a), du moteur : tri = On applique la relation fondamentale de la dynamique, en projection dans le référentiel mobile du drone, ce qui - en l'absence de vent - donne les trois équations suivantes : ù = (rv - qw) gsin0 Cxu = (pw - ru) + gsimpcos0 + Cyv (qu - pv) + gcoscpcos0 - Lin (Équations 1-3) p, q et r étant les vitesses angulaires suivant les trois axes, g étant l'accélération de la pesanteur, cp et 0 étant les deux angles définissant l'inclinaison du drone par rapport à l'horizontale (angles d'Euler : si correspond à une rotation autour de l'axe ZNED du repère absolu NED, 0 correspond à une rotation autour de l'axe YNED tourné de W, et (r) correspond à une rotation autour de l'axe u), C, et Cy étant les coefficients de résistance à l'avancement (reflétant les forces de frottement subies par le drone) dans les deux axes horizontaux, a étant un coefficient liant la poussée et la vitesse ascensionnelle à la vi- tesse de rotation w, et m étant la masse du drone. On applique de la même façon le théorème du moment dynamique au système, toujours en projection dans le référentiel mobile, ce qui conduit aux trois équations suivantes : + gr(1., - ly) = la(coz lyq + prOx 1.j = la(wï w) ILr + pq(ly i) = 1)(w - + w - (14 (Équations 4-6) I,, ly et I, étant des paramètres représentatifs du coefficient d'inertie du drone dans les trois axes, et 1 étant la distance séparant le moteur de son centre de gravité. Dans ces équations, le premier terme du membre de gauche correspond au moment dynamique du système, le second terme représente la contri- bution au moment dynamique des forces de Coriolis, et le membre de droite correspond aux moments exercés par les forces ascensionnelles F, et les couples F, créés par les hélices de chacun des rotors. Enfin, on démontre la relation suivante faisant intervenir les trois angles d'Euler cp, 0 et y : p + (osimp + rcos(p)tan0 é = ocoscp rsimp osincp + rcosq) :11 cos() (Équations 7-9) Le comportement du système est donc décrit au total par neuf équations à neuf inconnues. Au voisinage du point d'équilibre, avec le drone en sustentation à l'hori- zontale (vitesse et inclinaison nulles), on aura : V = 0 (f) 0 Les Équations 1-9 deviennent : 4 a2 CO , p = = 0, =- coi, 1 g Et donc, au voisinage du point d'équilibre : 1 \ling20 Si l'on pose wi = coi - coo, avec i = 1... 4, et que l'on linéarise à l'ordre 1 les équations précédentes autour du point d'équilibre, on aboutit au système d'équations linéarisées suivant : -9(0) - Cx * u = g(0) + CY * 2awo w (w1+ + W3 + / 4) 111, icuoo x lcutio = (w1 ) bw0 r - (w + w3 - W2 z = P = = r (Équations 10-18) Estimation du vent dans le plan horizontal L'invention propose une technique d'estimation du vent basée sur les équations dynamiques du drone en vol. Ce modèle dynamique, dont on donnera plus bas les équations, sera utilisé avec un estimateur d'état de type "filtre de Kalman", qui est un filtre à réponse impulsionnelle infinie qui estime les états d'un système dynamique (le drone dans le cas présent) à partir d'une série de mesures appliquées en entrée. On trouvera les prin- cipes généraux de cette technique par exemple dans R.E. Kalman, A new Approach to Linear Filtering and Prediction Problems, Transactions of the ASME - Journal of Basic Engineering, Vol. 82 (1960). Dans le cas présent, le filtre de Kalman utilisera et estimera six états, à savoir : - les composantes horizontales VD'fsx et VD/sy de la vitesse de dépla- cement du drone par rapport au sol (D/S), exprimées dans le repère local {u,v,w} lié au corps du drone, 1414) W4) - les composantes horizontales VA/ et VA sy de la vitesse de l'air par rapport au sol (A/S), exprimées dans le repère terrestre absolu NED fixe, et - les biais B'et B de l'accéléromètre sur ses axes X et Y. Le vecteur d'état s'écrit donc : I VA/ s; I VAls yi Bx By L'évolution des six états de l'estimateur du filtre de Kalman peut être prédite à l'aide du modèle dynamique décrivant le comportement du drone. Les valeurs mesurées de la vitesse par rapport au sol et de l'accélération du drone serviront à recaler les prédictions données par cet estimateur, permettant ainsi i) d'améliorer la précision d'estimation des grandeurs di- rectement mesurées (fonction de filtrage de l'estimateur), comme les vitesses VO -4;. et VDJS de déplacement du drone par rapport au sol, et ii) d'estimer les grandeurs inaccessibles à une mesure directe (fonction d'estimateur d'état), comme la vitesse et la direction du vent correspon- dant aux composantes VA/sx et VAisy . On démontre que le modèle dynamique en vol du drone peut s'exprimer sous la forme suivante : VD is ( VD /s y CX -(Oz 171 (0' 0 0 .)( représentant l'évolution de la valeur initiale Xdu vecteur d'état, l'entrée Udu système étant la gravité U= g, et Cy étant les coefficients de frottement aérodynamique du drone selon ses axes u et y, m étant la masse du drone, ço, 0 et y/ étant les angles d'Euler (respectivement de roulis, de tangage et de lacet) caractérisant l'attitude du drone par rapport au repère NED, (oz étant le taux de rotation autour de l'axe w (mouvement de rotation en lacet), et R4, étant la matrice de rotation associée à l'angle yr, en dimension 2. Les termes nuls des deux sous-matrices proviennent du fait que, par hy- pothèse, le vent et les biais de l'accéléromètre sont considérés comme constants sur le court terme. On notera que ceci n'empêche pas l'estimation de ces états car i) la vitesse du vent est couplée dynamiquement à la vitesse du drone par rapport au sol, elle-même mesurée et ii) le biais de l'accéléromètre est mesuré au travers de la mesure d'accélération. On notera par ailleurs que les deux matrices de la relation ci-dessus ne sont pas constantes (ce qui correspond à un filtre de Kalman non stationnaire), étant fonction de l'attitude du drone et de sa vitesse de rotation autour de l'axe de lacet ; mais ces paramètres sont supposés connus de l'estimateur, dans la mesure où ils sont en pratique mesurés ou estimés dans le cadre du contrôle d'attitude du drone. Dans la mesure où l'estimateur utilisé travaille en temps discret (par itérations successives), il est possible d'utiliser une version discrétisée du modèle dynamique en pratique, par approximation au premier ordre (relation de type "état suivant Xt+gt en fonction de l'état actuel Xt"), par approxima- tion au premier ordre : Xt-±st = Xt. St = X t + + BU). 61- Xt+st = (id + A St-)Xt + 6t. B. U D'où : Ad = 1d ± A. cît Bd = Ôt.B Subsidiairement, il est possible de prévoir un modèle dynamique différent représentant le comportement du drone lorsque celui-ci est au sol. En effet, le drone état posé, le vent n'a alors pas d'influence sur lui. Pour permettre aux estimations de revenir lentement à zéro, notamment pour ne pas déstabiliser les estimateurs de biais, on utilisera un modèle dynamique au sol du type : 0 0 0 0) 0 0 CI 0) -1 -1 (') 0 T 0 ( 1 0 Comme l'estimateur n'est dans ce cas pas en mesure de donner une estimation fiable de la vitesse du vent, ceci permet de forcer la décroissance de la vitesse du vent précédemment mesurée progressivement vers zéro par un modèle du premier ordre de paramètre r, temps caractéristique de cette décroissance, par exemple T = 1 seconde. On va maintenant décrire les mesures utilisées par l'estimateur. En ce qui concerne l'accélération, celle-ci est donnée par : 711 Y m. (I (C x ( C x (0 z iii ACCines C Z 0 Le dernier terme de cette relation sert à tenir compte du fait qu'un accéléromètre mesure les forces non gravitationnelles qu'il subit. On en tire ici implicitement parti : en effet, si l'on applique les Équations 10 et 11 au système de mesure accélérométrique, ce dernier ne mesure que les for- ces de frottements aérodynamiques ; l'accéléromètre mesure donc une grandeur proportionnelle à la vitesse du vent d'un facteur C,/m ou Cy/m selon l'axe considéré (au biais de construction près). La relation ci-dessus peut être réécrite en fonction du vecteur d'état : - lorsque le drone est en vol : cernesi,- (Accn, 'Accrnesz - et, de la même façon, lorsque le drone est au sol (posé) : o)) sin ( 1 -sin Op), Ac t, 0\ 0 mes UO 0) 0 0 En ce qui concerne la vitesse, celle-ci est mesurée dans le repère {u,v,w} lié au drone grâce à la caméra embarquée à visée verticale donnant une image du sol. Cette image est analysée par des algorithmes tels que ceux décrits dans le EP 2 400 460 Al précité, qui estime le déplacement de la scène captée par la caméra d'une image suivante et applique à ce dépla- cement estimé un facteur d'échelle fonction de l'altitude, elle-même esti- mée par fusion des données produites par un télémètre à ultrasons et par un capteur barométrique. La vitesse du drone exprimée dans son repère étant un élément du vecteur d'état, l'écriture de la mesure en fonction du vecteur d'état est di- recte : 0\ ° 0 ,0 °)),v + (°). u 0, Les Figures 3 et 4 sont des relevés comparatifs, respectivement du module de la vitesse du vent et de la direction du vent, des valeurs simultanément estimées et mesurées dans des conditions réelles. Sur ces figures, on a représenté en trait plein la mesure estimée du vent (vitesse, ou orientation exprimée en angle de cap), et en tiretés les valeurs réellement mesurées par un anémomètre et une girouette, respectivement. Ces figures permettent de vérifier l'excellente qualité de l'estimation obtenue par le procédé de l'invention, atteinte après quelques dizaines de secondes seulement, temps nécessaire à la stabilisation du filtre. Les mesures à t = 470-480 s et t > 560 s correspondent à un état du drone posé au sol, configuration qui ne permet pas de produire une esti- mation de la vitesse du vent. Toutefois, comme on l'a indiqué plus haut, l'utilisation d'un modèle dynamique spécifique pour le drone au sol permet de forcer la décroissance de la vitesse du vent vers zéro, évitant ainsi toute déstabilisation des estimateurs de biais des accéléromètres. Utilisation de l'estimation de la vitesse du vent - compensation du vent L'estimation du vent a plusieurs utilités : - pour réaliser une compensation permettant de stabiliser le drone contre le vent, en mode contrôlé ou en mode autopiloté ; - pour réaliser une fusion de la vitesse et de l'accélération mesurées te- nant compte du vent estimé, pour le contrôle du point fixe autopiloté ; - pour estimer le biais de construction de l'accéléromètre fourni à l'estimateur d'attitude, permettant ainsi d'obtenir une estimation de l'attitude efficace, non perturbée par le vent même en présence de ce dernier. On va maintenant décrire un procédé de compensation du vent. Ce procédé de compensation du vent est avantageusement mis en oeuvre à partir de l'estimation obtenue de la manière que l'on a décrite plus haut, mais de façon non limitative, dans la mesure où il peut être mis en oeuvre avec d'autres techniques d'estimation du vent. La compensation du vent consiste à déterminer pour le drone une assiette de référence autour de laquelle il va s'asservir, soit en mode de vol contrôlé soit en mode de point fixe autopiloté : - en vol contrôlé, la référence en boucle ouverte sera ajoutée aux consi- gnes du pilote, de manière à soutenir le drone face au vent et inverse- ment le freiner dans le sens du vent, et - lors du passage au point fixe, la stabilisation du drone pourra être atteinte plus rapidement, la compensation étant aussi réactive que l'estimation. La Figure 5 est une représentation schématique fonctionnelle illustrant les boucles d'estimation et de compensation du vent suivant l'axe u, pour une configuration de vol du drone en point fixe autopiloté. L'asservissement est identique selon l'axe y, les termes en x ou en 8 étant substitués par des termes en y ou en ço et la fonction cosinus étant remplacée par la fonction sinus, respectivement. La Figure 6 est homologue de la Figure 5, pour une configuration de vol du drone en mode contrôlé. Sur ces deux figures 5 et 6, on retrouve les blocs 124 (commandes utilisa- teur appliquées en mode contrôlé), 126 (correction du point fixe en mode autopiloté) et 128 (estimation du vent à partir des mesures d'attitude et de vitesse, l'estimation du vent incluant également la décorrélation des biais des accéléromètres) illustrés sur le schéma d'ensemble de la Figure 2. La boucle de compensation en point fixe 126 comprend un étage intégra- teur 154 à gain 156 et un étage proportionnel 158. La consigne en sortie de l'étage additionneur 152 est appliquée par l'intermédiaire d'un étage 160 limitant l'excursion de la consigne effectivement appliquée à la boucle de contrôle des moteurs du drone. Le terme de compensation appliqué à l'additionneur 152 pour chacun des axes u et y est : Cx ev,._ _comp Vwfnd. COS 01 - P dro ) (/) wi n d _comp = nd - SinOliwind 11) drone)- Vwind étant le module de la vitesse estimée du vent, Cx étant le coefficient de traînée du drone, vwind étant l'angle de cap de la direction du vent et /1drone étant l'angle de cap de la direction du drone. REVENDICATIONS1. Un procédé de pilotage d'un drone (10) à voilure tournante à rotors multiples (12) entraînés par des moteurs respectifs (110) commandables pour piloter le drone en attitude et en vitesse, comprenant : - la génération de consignes angulaires (0, g)) et l'application de ces consignes à une boucle (120) de contrôle des moteurs du drone, ces consignes étant aptes à commander l'attitude du drone selon des axes de tangage (22) et de roulis (24), caractérisé en ce qu'il comprend en outre : - l'établissement d'au moins un modèle dynamique du drone, décrivant des composantes de vitesse horizontale du drone en fonction des coef- ficients de traînée et de la masse du drone, des angle d'Euler caractérisant l'attitude du drone par rapport à un repère terrestre absolu, ainsi que de la rotation du drone autour d'un axe vertical ; - la mesure de la force de traînée aérodynamique du drone ; - la mesure de la vitesse relative du drone par rapport au sol ; et - l'application audit modèle dynamique du drone desdites mesures de force de traînée aérodynamique et de vitesse relative par rapport au sol, de manière à produire une estimation des composantes horizontales de vitesse du vent latéral. 2. Le procédé de la revendication 1, dans lequel ladite mesure de la force de traînée aérodynamique du drone est dérivée d'une mesure d'accélération du drone. 3. Le procédé de la revendication 1, dans lequel ladite étape d'application des composantes horizontales de vitesse de vent latéral au modèle dynamique du drone comprend l'utilisation d'un filtre prédictif de Kalman. 4. Le procédé de la revendication 3, dans lequel le filtre prédictif de Kal- man est un filtre à six états, ces états comprenant : - deux composantes de la vitesse de déplacement du drone par rapport au sol, exprimées dans un repère lié au drone,- deux composantes de la vitesse du vent par rapport au sol, exprimées dans un repère terrestre absolu, et - deux composantes horizontales de biais d'accéléromètre du drone. 5. Le procédé de la revendication 1, dans lequel l'établissement d'au moins un modèle dynamique du drone comprend l'établissement de deux modèles distincts, avec : i) un modèle dynamique du drone en vol, et ii) un modèle dynamique du drone au sol, ces modèles étant utilisés sélectivement en fonction de l'état du drone. 6. Le procédé de la revendication 1, dans lequel le modèle dynamique du drone en vol utilise une mesure d'accélération du drone comme mesure proportionnelle à la vitesse de l'air dans le repère du drone. 7. Le procédé de la revendication 5, dans lequel le modèle dynamique du drone en vol est du type : C tc,,, m (Oz (0 0\ ' -sin (0 Cy 0) sin (I)) - (oz - 0 - 0 0' 0 0 0\ (00) 0 0) 0 0 '0 0/ (0\ (\, 0 0 0 o ') o 0) Co 0 0) ) '0) 0 0\ 'k:' représentant l'évolution de la valeur initiale Xdu vecteur d'état, l'entrée U du système étant la gravité U= g, C., et G étant les coefficients de frottement aérodynamique du drone selon ses axes u et y, m étant la masse du drone, 0 et u étant les angles d'Euler (respectivement de roulis, de tangage et de lacet) caractérisant l'attitude du drone par rapport au repère NED, co, étant le taux de rotation autour de l'axe w (mouvement de rotation en lacet), et Rw étant la matrice de rotation associée à l'angle v, en dimension 2. 8. Le procédé de la revendication 5, dans lequel le modèle dynamique du drone au sol est du type : 0 (0 0 \,0 -1 (°0 - (00 T 0 0 0) X + °) 0) 0 0 0 0' où -c est un paramètre temporel de décroissance progressive de l'estima- tion de la vitesse du vent. 9. Le procédé de la revendication 1, comprenant en outre une étape de compensation de l'effet du vent latéral sur le positionnement et les déplacements du drone, par : - génération de consignes correctrices, fonction des composantes hori- zontales estimées de vitesse du vent latéral, et - combinaison de ces consignes correctrices aux consignes angulaires de tangage et de roulis appliquées à la boucle de contrôle des moteurs du drone. 10. Le procédé de la revendication 9, dans lequel les consignes correctrices sont du type : cos ?iddrone) Cx qnd_comp sin( d drone_ Vwind étant le module de la vitesse estimée du vent, Cx étant le coeffi- cient de traînée du drone, v,,,'d étant l'angle de cap de la direction du vent et Vciione étant l'angle de cap de la direction du drone. 11. Le procédé de la revendication 9, dans lequel l'étape de compensation comprend en outre la définition pour le drone d'une assiette de référence en boucle ouverte. 12. Le procédé de la revendication 9, dans lequel lorsque le drone est en vol contrôlé par un utilisateur, les consignes correctrices (Owind_co,np, ÇOwitid_comp) sont combinées aux consignes de pilotage (0 pilote, 0'12ilote) appliquées par l'utilisateur. 13. Le procédé de la revendication 9, dans lequel lorsque le drone est en vol stationnaire autopiloté, les consignes correctrices (OH,ind_comp, ÇOwincl_comp) sont combinées à des consignes de stabilisation de point fixe (0p, p) produites en réponse à une mesure de la vitesse horizontale du drone par rapport au sol ( VD rzX, VD/sy ).