Brevet FR2818383A1 — publié le 20020621

ANEMOMETRE POUR HELICOPTER L'invention concerne la mesure de la vitesse d'un hélicoptère par rapport à l'air. La mesure de la vitesse des aéronefs se fait en général à partir de centrales anémobarométriques qui permettent de calculer la vitesse de l'aéronef par rapport à l'air ambiant à partir de mesures de pression statique, de pression totale et de température : à pression statique et température statique constante, la pression totale mesurée croît et décroît comme la vitesse selon une loi connue. Ces centrales sont bien adaptées aux avions sur lesquels on sait placer des sondes de mesure de pression et température, mais elles sont mal adaptées aux hélicoptères, à cause des basses vitesses, des fortes composantes de déplacement latéral, et des perturbations de mesure de pression dues aux flux d'air du rotor. C'est pourquoi on a déjà proposé de placer des sondes de pression sur les pales du rotor et d'observer les variations de pression. Si l'hélicoptère ne se déplace pas par rapport à l'air ambiant, la vitesse de la sonde par rapport à l'air est constante (le vecteur vitesse tourne seulement avec le rotor) ; la pression mesurée par la sonde est constante et sa valeur est liée uniquement à la vitesse de rotation du rotor. Si l'hélicoptère se déplace à vitesse Vair par rapport à l'air, l'amplitude de la vitesse instantanée de la sonde par rapport à l'air varie sinusoïdalement autour d'une valeur moyenne qui représente la vitesse de rotation des pales ; l'amplitude de la sinusoïde représente la vitesse d'avancement de l'hélicoptère par rapport à l'air ; et la phase représente la direction d'avancement par rapport à une direction de référence de l'hélicoptère : le maximum de la sinusoïde de vitesse se produit au moment où la sonde est tournée vers la direction d'avance de l'hélicoptère, tandis que le minimum se produit au moment où la sonde est tournée dans la direction opposée. La pression totale mesurée par la sonde ne varie pas linéairement avec la vitesse. Mais elle suit une loi de variation connue, classique en anémobarométrie. Lorsque la vitesse varie sinusoïdalement, la pression suit alors une courbe de variation qu'on peut qualifier de sinusoïde déformée. Si la variation de pression suivait une loi exactement sinusoïdale, identique à la variation sinusoïdale de la vitesse de la sonde, il serait facile de déterminer l'amplitude et la phase de cette sinusoïde et d'en déduire l'amplitude et la direction du vecteur vitesse de l'hélicoptère. Mais ce n'est pas le cas du fait que la sinusoïde est déformée. On a déjà proposé dans le passé d'effectuer une analyse par transformée de Fourier de la courbe de pression mesurée en bout de pale, pour en tirer des valeurs numériques qui permettent de remonter jusqu'à la valeur et à la direction du vecteur vitesse par rapport à l'air. Mais ce procédé réclame des calculs lourds et donc des moyens électroniques sophistiqués pour obtenir en temps réel une indication précise du vecteur vitesse. L'invention a pour but de proposer une solution particulièrement simple pour déterminer l'amplitude et la direction du vecteur vitesse d'un hélicoptère. Selon l'invention, on propose un procédé de détermination de vitesse d'hélicoptère utilisant une sonde de pression montée sur une pale du rotor, ce procédé consistant à : - obtenir N échantillons de mesure de pression Pi répartis sur l'ensemble d'une rotation complète de la pale, - produire un premier signal binaire changeant d'état au passage de la pale par des orientations remarquables situées à 0'eut 1800 d'une orientation de référence de la pale, et un deuxième signal binaire en quadrature de phase avec le premier et changeant d'état à 90'eut 270', et effectuer les sommes suivantes S1 et S2 pendant une rotation complète : - S1 : somme des échantillons obtenus pendant un premier état du premier signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du premier signal ; - S2 : somme des échantillons obtenus pendant un premier état du deuxième signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du deuxième signal binaire, - calculer une valeur numérique représentant arctg (S2/S1), le résultat de ce calcul déterminant une phase remarquable correspondant à une direction de la vitesse cherchée. Pour déterminer ensuite l'amplitude de la vitesse, on prévoit de stocker les N échantillons de pression utilisés au cours de la détermination de la phase remarquable ainsi que les instants de mesure correspondants ; on additionne une valeur n/2 et une valeur 3n/2 à la phase (modulo 2), pour obtenir deux autres phases remarquables Pmax et cumin représentant les instants de mesure théoriques du maximum et du minimum de la courbe de valeurs de pression mesurées au cours de la rotation, et on recherche dans les échantillons stockés les valeurs d'échantillons de pression Pmax et Pmin correspondant à ces instants. Enfin, on calcule la vitesse instantanée Vmax et Vmin à partir de ces valeurs de pression, et on détermine l'amplitude de la vitesse à partir de Vmax et Vmin : la vitesse Vair est en principe la moitié de la différence Vmax-Vmin. L'invention propose en outre un système de mesure de vitesse comportant les moyens nécessaires à la mise en oeuvre du procédé qui vient d'être décrit dans ses grandes lignes. Le système de mesure'de vitesse d'hélicoptère par rapport à l'air, comprend : - une sonde de pression placée sur une pale du rotor de l'hélicoptère, cette sonde fournissant des signaux représentant la pression mesurée par la sonde, - des moyens pour obtenir des échantillons de mesure de pression répartis au cours d'une rotation complète de la pale, - des moyens pour produire un premier signal binaire changeant d'état au passage de la pale par des orientations remarquables situées à 0'eut 1800 d'une orientation de référence de la pale, et un deuxième signal binaire en quadrature de phase avec le premier et changeant d'état à 900 et 270 , - un premier accumulateur numérique pour calculer la somme S1 des échantillons obtenus pendant un premier état du premier signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du premier signal ; - un deuxième accumulateur numérique pour calculer la somme S2 des échantillons obtenus pendant un premier état du deuxième signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du deuxième signal binaire, - des moyens pour calculer une valeur numérique représentant arctg (S2/S1), le résultat de ce calcul déterminant une phase remarquable correspondant à une direction de la vitesse cherchée. Le système comprend de préférence également, pour la détermination de l'amplitude de la vitesse, les moyens suivants : - des moyens pour stocker les échantillons de pression en relation avec les instants de mesure correspondants, - des moyens pour additionner une valeur jar/2 et une valeur 3/2 à la phase La sonde de pression peut être une sonde fournissant des impulsions dont la périodicité est fonction de la pression, et il est alors prévu de préférence un circuit de détermination de la valeur de cette périodicité, ce circuit fournissant des valeurs d'échantillons de mesure de pression représentés par la valeur de la période déterminée. Ce circuit peut être un simple circuit de comptage d'impulsions à fréquence beaucoup plus élevée que celle des impulsions fournies par la sonde, ce circuit comptant le nombre d'impulsions de fréquence élevée émises pendant une ou plusieurs périodes d'impulsions issues de la sonde. Les accumulateurs qui calculent les sommes S1 et S2 peuvent alors recevoir simplement pendant ce temps la fréquence de comptage de ce compteur, et compter ou décompter les impulsions reçues à fréquence élevée selon l'état des premier et deuxième signaux binaires. Les accumulateurs sont remis à zéro au bout d'une rotation complète de la pale, et c'est alors qu'est effectué le calcul d'arctg (S2/S1). D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit et qui est faite en référence aux dessins annexés dans lesquels : - la figure 1 représente une hélicoptère dont une des pales de sustentation porte une sonde de pression, la pale étant dans une position de référence par rapport à la structure de l'hélicoptère ; - la figure 2 représente l'hélicoptère avec la pale dans une position quelconque au cours de sa rotation ; - la figure 3 représente la sinusoïde de vitesse de la sonde ; - la figure 4 représente la sinusoïde déformée de pression mesurée par la sonde ; - la figure 5 représente le principe de corrélation permettant de déterminer la phase de la sinusoïde déformée ; - la figure 6 représente un schéma bloc d'un circuit de mise en oeuvre préférentiel du procédé de mesure de vitesse selon l'invention. Sur la figure 1, on voit l'hélicoptère 10 avec quatre pales de sustentation entraînées par un rotor, parmi lesquelles une pale 12 porte une sonde de pression 14 ; cette pale sera considérée comme pale de référence. Les autres pales pourraient aussi porter des sondes de pression, et dans ce cas on pourrait combiner les informations issues des différentes sondes pour accroître la précision ou la fiabilité du système. De même, il pourrait y avoir plusieurs sondes sur une pale. Dans ce qui suit on considérera qu'il n'y a qu'une sonde de pression sur une seule pale puisque cela suffit dans le principe à la détermination de vitesse de l'hélicoptère. La sonde peut être située en bout de pale ou dans une partie intermédiaire. En bout de pale, elle doit alors être plus résistante aux accélérations fortes, et notamment à la force centrifuge. La sonde est tournée de préférence vers l'avant de la pale, dans le sens de rotation de la pale (perpendiculairement à la longueur de la pale par conséquent) pour mesurer la pression d'air maximale due à la vitesse de rotation de la pale plutôt que la dépression à l'arrière de la pale. Le sens de rotation est le sens des aiguilles d'une montre sur les figures. La sonde doit être également placée dans un plan correspondant à l'assiette de référence de l'hélicoptère, de manière que la vitesse d'hélicoptère calculée soit référencée par rapport à cette assiette. A la figure 1, la pale 12 est représentée dans une position de référence utilisée dans le calcul de vitesse. Par souci de simplicité, on prend ici comme référence la position de la pale 12 dans l'axe longitudinal de la structure de l'hélicoptère et vers l'arrière, la vitesse de l'hélicoptère par rapport à l'air étant définie par rapport à cet axe longitudinal. A la figure 2, on a représenté la pale 12 en rotation à vitesse n supposée constante, la pale étant dans une position quelconque définie à un instant t par l'angle nt de la pale par rapport à sa position de référence (on suppose qu'à t=O, la pale est dans la position de référence). Une flèche désignée par Va, r représente la vitesse de l'hélicoptère par rapport à l'air. Elle fait un angle avec l'axe longitudinal de référence de l'hélicoptère. La détermination de la vitesse comprend la détermination de l'angle et de l'amplitude Vair. La vitesse linéaire instantanée de la sonde par rapport à l'hélicoptère est un vecteur Vi d'amplitude LQ, si L est la distance entre la sonde et l'axe du rotor, et de direction #t+ 7r/2. Mais la vitesse de la sonde par rapport à l'air, sachant que l'hélicoptère avance avec la vitesse Vair dans la direction définie par l'angle ()), est un vecteur V'i qui est la combinaison vectorielle de Vi et de Vair. Ce vecteur combiné V'i peut être décomposé en deux composantes, l'une dans la direction longitudinale de la pale, sans effet sur la sonde de pression, et l'autre, désignée par V"i dans la direction perpendiculaire à la pale, direction où la sonde est le plus sensible à la pression. La pression vue à l'instant t par la sonde est celle qui résulte de la composante de vitesse instantanée V"i = Ç2L-Vairsin (nt-). Cette composante de vitesse varie donc selon une courbe sinusoïdale oscillant avec une amplitude de +Vair à-Vair autour d'une valeur moyenne nL et avec une phase (déterminée par rapport à la position de référence de la pale) qui est l'angle # représentant la direction d'avancement de l'hélicoptère. Cette courbe est représentée à la figure 3. Si la mesure de pression par la sonde permettait de calculer directement cette courbe, la mesure de vitesse consisterait simplement à repérer l'amplitude de l'oscillation et sa phase pour connaître le vecteur vitesse Va, r en amplitude et direction. Toutefois les sondes de pression ne fournissent pas une vitesse mais une pression, et la formule mathématique (connue et classique en anémobarométrie) qui permet de passer de la vitesse à la pression n'est pas une formule linéaire. C'est une fonction du type Pt = Ps (l +0, 2M2) 7/2 où Pt est la pression totale mesurée par la sonde, Ps est la pression statique atmosphérique (en dehors de tout mouvement), M est la vitesse de la sonde exprimée en nombre de Mach, c'est-à-dire en fraction de la vitesse du son dans l'air, la vitesse du son dépendant elle-même de la température de l'air. De plus, la sonde elle-même a des caractéristiques qui ne sont pas forcément linéaires. Par exemple, une sonde de type capacitif, à membrane résonante, produisant un signal de sortie sous forme d'une valeur de fréquence de résonance, fournit des échantillons de mesure de pression qui ne varient pas linéairement en fonction de la pression. Typiquement, on peut utiliser une sonde fournissant une fréquence de sortie F = Fo [ (P+Po)/2Po] , P étant la pression, Po une pression au repos et Fo la fréquence correspondante. Po dépend en pratique de la température du capteur. La pression mesurée par la sonde, varie dans ces conditions selon une sinusoïde déformée représentée à la figure 4, dont les maxima et minima se situent aux mêmes endroits que ceux de la sinusoïde de vitesse, mais dont la phase Dans ce qui suit on fera le raisonnement mathématique en parlant d'échantillons de pression, que ces échantillons de pression soient véritablement une pression ou une fréquence représentant cette pression, ou encore la période de cette fréquence, ou tout autre signal électrique représentant la pression, pourvu que ce signal varie de manière monotone avec la pression. Le raisonnement est en effet le même dans tous les cas, en effectuant cependant les transpositions nécessaires si le signal varie en sens inverse de la pression : par exemple si les échantillons de pression sont établis sous forme d'un comptage de la période de signal et non pas de la fréquence de signal dans le cas de la sonde mentionnée ci-dessus, il est évident que les maxima de la période coïncident avec les minima de pression (et de vitesse) de la sonde, et vice-versa. On pourrait bien sûr calculer point par point toutes les valeurs de vitesse V"i correspondant à toutes les pressions Pi mesurées, à partir du calcul inverse de celui qui a été indiqué ci-dessus, puis, ayant retrouvé une forme sinusoïdale, déterminer assez facilement l'amplitude et la direction de la vitesse Vair. Mais cela demanderait des moyens de calcul très puissants. Un traitement par transformée de Fourier sur les échantillons de pression obtenus au cours d'une rotation complète permet également de remonter à la phase et à la vitesse Vair mais demande également des moyens de calcul puissants. La présente invention propose un moyen beaucoup plus simple qui repose sur le fait que bien que la sinusoïde soit déformée, elle est fortement corrélée avec la forme de la sinusoïde de vitesse qui a même période et même phase, les maxima et minima de la sinusoïde déformée coïncidant avec ceux de la sinusoïde de vitesse (ou éventuellement intervertis) ; on devrait donc pouvoir, par une corrélation simple, retrouver la phase de la sinusoïde de vitesse correspondant à la sinusoïde de pression ou à la sinusoïde des échantillons de signal quel qu'il soit. La sonde fournit donc à chaque tour de rotor N échantillons périodiques de mesure de pression dont les valeurs sont stockées en mémoire et associées chacune à la valeur de l'instant de mesure (ce stockage est utile à la mesure d'amplitude mais non à la mesure de la phase) ; cet instant est exprimé par exemple sous forme d'un angle de phase 2k7t/N où k représente le numéro de l'échantillon de 1 à N. Par exemple, c'est simplement le numéro d'échantillon k qui sert d'adresse dans une mémoire où on stocke la valeur de l'échantillon à cette adresse. On fait la corrélation de la manière la plus simple possible, consistant alors simplement à créer deux formes d'onde carrées, alternant chacune entre deux valeurs binaires (qu'on désignera par simplicité par +1 et-1). Ces deux formes d'onde carrée ont pour période la période de rotation du rotor, sont synchronisées sur cette rotation, et sont déphasées de 900 l'une par rapport à l'autre. La première forme d'onde a par exemple la valeur +1 pendant que l'angle de rotation Dt est supérieur ou égal à 0"eut inférieur à 180 , puis la valeur-1 tant que l'angle est supérieur ou égal à 1800 et inférieur à 3600. La deuxième forme d'onde prend alors la valeur-1 de 00 à 900 et de 2700 à 3600 et la valeur +1 de 900 à 2700. Les valeurs numériques des échantillons de mesure de pression sont multipliés par +1 ou par-1 selon la valeur binaire de la forme d'onde, et additionnés de manière à former les sommes algébriques suivantes : Somme S1 : somme des échantillons obtenus pour un premier état de la première forme d'onde (+1 en principe), diminuée de la somme des échantillons obtenus pour le deuxième état de la première forme d'onde Somme S2 : somme des échantillons obtenus pour un premier état de la deuxième forme d'onde (+1 en principe), diminuée de la somme des échantillons obtenus pour le deuxième état de la deuxième forme d'onde. La figure 5 représente les formes d'onde carrées F1 et F2, calées en phase sur la rotation du rotor, en regard de la courbe de valeurs de pression transmises par la sonde. La théorie et l'expérimentation montrent qu'on peut tirer l'angle L'angle Etant donné que la fonction arctg admet deux solutions séparées de 1800 pour l'angle, on lève l'ambiguïté sur la valeur de l'angle par la connaissance du signe de la somme S1 ou S2 par exemple (positif si l'angle est inférieur à 180 ). Par conséquent, on considère que la valeur arctg (S2/S1) représente le déphasage qu'aurait la sinusoïde de vitesse correspondant à la courbe de pression mesurée par la sonde, ceci sans avoir besoin de recalculer les points de la sinusoïde de vitesse. Les valeurs de phase indiquées ci-dessus pour les signaux binaires F1 et F2 concernent le cas où les échantillons de pression varient comme la vitesse de la sonde. S'ils varient en sens inverse, on les modifiera, par exemple en inversant la phase de F1 qui aurait alors la valeur-1 de 0 à 1800. Connaissant maintenant cet angle de déphasage , représentant la direction du vecteur vitesse de l'hélicoptère par rapport à l'air, on calcule les angles Pmax = 4 > +n/2 (modulo 2n) et +3n/2 (modulo 2n) auxquels se situent le maximum et le minimum de la courbe de pression puisque ces maximum et minimum coïncident avec ceux de la sinusoïde de vitesse. Et on va chercher dans la table des valeurs d'échantillons les valeurs de pression qui correspondaient aux instants de mesure pour ces angles. En pratique, on exprime la phase calculée 4 > sous forme d'une numéro d'échantillon ko et on exprime les positions de (I) et et (min sous la forme ko+N/4 et k0 + 3N/4 modulo N. Ces nombres représentent l'adresse à laquelle on va chercher les échantillons représentant la pression maximale Pmax et la pression minimale Pmin. Si ces échantillons sont stockés sous forme d'une valeur de fréquence d'impulsions ou sous forme d'une valeur'de période, on effectue évidemment le calcul qui convient pour trouver la pression. Dans le cas de la sonde capacitive, la formule simplifiée de la fréquence en fonction de la pression a été indiquée. Toutefois, la valeur de pression de référence dépend de la température et un capteur de température de la sonde est en général associé à la sonde pour permettre le calcul de la pression. Ayant les valeurs de pression maximale et minimale, on calcule alors les vitesses air correspondantes par la formule inverse de la formule Pt=Ps (1 +0, 2M pour aboutir à deux vitesses Vi"max et Vi"min. La différence de ces vitesses est égale à deux fois la vitesse de l'hélicoptère par rapport à l'air. Vair = (Vi"max-Vi"min)/2 Le calcul de vitesse fait intervenir classiquement la pression statique et la température de l'air, qui doivent être connues par ailleurs. La mise en oeuvre de l'invention peut se faire à l'aide de la circuiterie représentée sur la figure 6. On a illustré l'invention dans le cas d'une sonde de pression fournissant des impulsions dont la fréquence (donc la période) varie en fonction de la pression. On pourrait utiliser pour constituer des échantillons de mesure de pression, un simple comptage des impulsions pendant une durée déterminée, et se servir du comptage pour effectuer les sommes algébriques S1 et S2 définies précédemment. Cependant, pour plus de précision lorsque la fréquence des impulsions de la sonde n'est pas très élevée (quelques dizaines de kilohertz), on préfère utiliser comme valeur d'échantillons un comptage de la durée de la période des impulsions fournies par la sonde. Ce comptage est fait avec une fréquence de comptage fixe beaucoup plus grande que la fréquence des impulsions fournies par la sonde. On peut compter la période séparant deux impulsions successives, ou mieux la durée de plusieurs impulsions successives en provenance de la sonde. Une horloge générale HG, synchronisée sur la rotation du rotor qui en général est régulière, définit les N instants de mesure, par exemple sous forme d'impulsions à fréquence NQ. Elle définit aussi les signaux binaires F1 et F2 servant à la corrélation. A un instant de mesure donné, t=k/NQ on détermine le début de la première impulsion de sonde qui arrive, et on déclenche le comptage d'un compteur CPT1 qui compte à partir de zéro à une fréquence f, fixe et très supérieure à celle des impulsions de la sonde (par exemple f = 20 MHz). On arrête le comptage à l'arrivée de l'impulsion suivante, ce qui donne un contenu de compteur représentant directement cette période, donc l'inverse de la fréquence, donc une valeur représentant la pression. Ou bien encore, on arrête le comptage au bout de plusieurs impulsions (toujours le même nombre bien sûr, défini par exemple par un compteur CPT2) ce qui donne un compte représentant un somme de plusieurs périodes, représentant également la pression. La valeur du compteur en fin de comptage est stockée dans une mémoire MEM, à une adresse k correspondant à l'instant de mesure. Cette adresse est incrémentée à chaque nouvel échantillon par l'horloge générale HG qui définit les instants de mesure. Une valeur stockée peut par ailleurs être lue ultérieurement en appliquant une adresse de lecture correspondant à un instant de mesure déterminé. La mémoire est complètement remplie à chaque tour de rotor en remplaçant les anciennes valeurs. Des mémoires tampons pour le remplissage de la mémoire, ou un fonctionnement avec deux mémoires travaillant en alternance en lecture et écriture à chaque tour est bien entendu prévu si la durée nécessaire au calcul de vitesse exige de conserver le contenu de la mémoire d'un tour avant de la remplir pour le tour suivant. Les créneaux binaires F1 et F2 contrôlent chacun un accumulateur respectif ACC1, ACC2. Les accumulateurs sont initialisés à zéro lors du passage rotor par sa position de référence. Ils additionnent les échantillons au fur et à mesure de leur arrivée, en positif ou en négatif selon l'état des signaux binaires F1, F2 qui les commandent. Ils calculent donc les sommes S1 et S2 qui sont disponibles dans les accumulateurs à la fin d'un tour pour être utilisées dans un circuit de calcul de arctg (S2/S1). Les accumulateurs pourraient être des additionneurs/soustracteurs qui reçevraient les échantillons du compteur CPT1 au fur et à mesure qu'ils arrivent et les additionneraient ou les soustrairaient en fonction de l'état de F1 et F2. Mais, comme cela est représenté sur la figure 5, il est plus simple de prévoir que les accumulateurs ACC1 et ACC2 sont de simples compteurs/décompteurs qui reçoivent directement les impulsions de comptage du compteur CPT1 (en principe à fréquence f ou sous-multiple de f) tant que le compteur compte et qui les accumulent algébriquement en comptant ou décomptant ces impulsions selon l'état de F1 et F2. A la fin d'un tour de rotor, les accumulateurs contiennent les sommes S1 et S2 et les transfèrent à un circuit de calcul de arctg (S2/S1) avant d'être remis à zéro pour un nouveau tour. Le circuit de calcul de = arctg (S2/S1) fournit directement en sortie la direction du vecteur vitesse de l'hélicoptère. De plus, il fournit une information à un circuit de calcul de +90 et +270 . Ce dernier fournit à la mémoire MEM, en lecture, les adresses ko +N/4 et ko+3N/4 correspondant à ces phases. Les valeurs d'échantillons stockés sont extraits de la mémoire et transmis à un circuit de calcul de pression Pmin et Pmax à partir des valeurs d'échantillons stockés. Le circuit de calcul de Pmin et Pmax utilise les données caractéristiques de la sonde de pression pour déduire à partir des nombres représentant les périodes d'impulsions issues de la sonde des valeurs numériques de pression. Ce circuit peut faire intervenir une mesure de température de la sonde, et dans ce cas une information de température de la sonde est fournie par celle-ci (cas des sondes de pression associées à une sonde de température locale). Enfin, à partir des valeurs de pression Pmin et Pmax, un circuit de calcul de Vmin et Vmax détermine Vmin et Vmax en fonction des formule anémobarométriques indiquées ci-dessus. Ce calcul fait en principe intervenir une mesure de pression statique et une valeur de température de l'air (température statique ou température totale selon les définitions en usage dans l'anémobarométrie). Les valeurs de pression statique et de température, obtenues par ailleurs, sont donc transmises à ce circuit. Le circuit de calcul de Vmin et Vmax fournit la valeur Vair qui est la moitié de la différence entre Vmin et Vmax. REVENDICATIONS 1. Procédé de détermination de vitesse d'hélicoptère utilisant une sonde de pression (14) montée sur une pale (12) du rotor, ce procédé consistant à : - obtenir N échantillons de mesure de pression Pi répartis sur l'ensemble d'une rotation complète de la pale, - produire un premier signal binaire (F1) changeant d'état au passage de la pale par des orientations remarquables situées à 0 et 1800 d'une orientation de référence de la pale, et un deuxième signal binaire (F2) en quadrature de phase avec le premier et changeant d'état à 90"eut 270", et effectuer les sommes suivantes S1 et S2 pendant une rotation complète : - S1 : somme des échantillons obtenus pendant un premier état du premier signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du premier signal ; - S2 : somme des échantillons obtenus pendant un premier état du deuxième signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du deuxième signal binaire, - calculer une valeur numérique représentant arctg (S2/S1), le résultat de ce calcul déterminant une phase remarquable correspondant à une direction de la vitesse cherchée. 2. Procédé selon la revendication 1, caractérisé en ce que, pour déterminer ensuite l'amplitude de la vitesse, on prévoit de stocker les N échantillons de pression utilisés au cours de la détermination de la phase remarquable ainsi que les instants de mesure correspondants ; on additionne une valeur 7c/2 et une valeur 3/2 à la phase (modulo 2n), pour obtenir deux autres phases remarquables Pmax et Omin représentant les instants de mesure théoriques du maximum et du minimum de la courbe de valeurs de pression mesurées au cours de la rotation, on recherche dans les échantillons stockés les valeurs d'échantillons de pression Pmax et Pmin correspondant à ces instants, et enfin on calcule la vitesse instantanée Vmax et Vmin à partir de ces valeurs de pression, en fonction d'une relation vitesse/pression, et on détermine l'amplitude de la vitesse de l'hélicoptère à partir de Vmax et Vmin. 3. Procédé selon la revendication 2, caractérisé en ce que la relation vitesse/pression est la suivante : Pt = Ps (1+0, 2M2) 7/2 où Pt est la pression totale mesurée par la sonde, Ps est la pression statique atmosphérique, M est la vitesse exprimée en nombre de Mach. 4. Système de mesure de vitesse d'hélicoptère comprenant une sonde de pression (14) placée sur une pale (12) du rotor de l'hélicoptère, cette sonde fournissant des signaux représentant la pression (Pi) mesurée par la sonde, et des moyens pour obtenir des échantillons de mesure de pression répartis au cours d'une rotation complète de la pale, caractérisé en ce qu'il comprend en outre : - des moyens (HG) pour produire un premier signal binaire (S1) changeant d'état au passage de la pale par des orientations remarquables situées à 0 et 1800 d'une orientation de référence de la pale, et un deuxième signal binaire (F2) en quadrature de phase avec le premier et changeant d'état à 900 et 2700, - un premier accumulateur numérique (ACC1) contrôlé par le premier signal binaire pour calculer la somme S1 des échantillons obtenus pendant un premier état du premier signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du premier signal ; - un deuxième accumulateur numérique (ACC2) contrôlé par le deuxième signal binaire pour calculer la somme S2 des échantillons obtenus pendant un premier état du deuxième signal binaire diminuée de la somme des échantillons obtenus pendant le deuxième état du deuxième signal binaire, - des moyens pour calculer une valeur numérique représentant arctg (S2/S1), le résultat de ce calcul déterminant une phase remarquable correspondant à une direction de la vitesse cherchée. 5. Système selon la revendication 4, caractérisé en ce qu'il comprend également, pour la détermination de l'amplitude de la vitesse : -des moyens (MEM) pour stocker les échantillons de pression en relation avec les instants de mesure correspondants, - des moyens pour additionner une valeur 7 : 12 et une valeur 3Z/2 à la phase 6. Système selon l'une des revendications 4 et 5, caractérisé en ce que la sonde de pression est une sonde fournissant des impulsions dont la périodicité est fonction de la pression, et en ce qu'il est prévu un circuit (CPT1) de comptage d'impulsions à fréquence beaucoup plus élevée que celle des impulsions fournies par la sonde, ce circuit comptant le nombre d'impulsions de fréquence élevée émises pendant une ou plusieurs périodes d'impulsions issues de la sonde, et le contenu du compteur représentant une valeur d'échantillon de pression. 7. Système selon la revendication 6, caractérisé en ce que les accumulateurs (ACC1, ACC2) sont des compteurs/décompteurs qui reçoivent les impulsions d'incrémentation du circuit de comptage d'impulsions (CPT1) et qui s'incrémentent en comptage ou en décomptage selon l'état du signal binaire qui les contrôle, les sommes S1 et S2 étant établies au bout d'un tour complet du rotor et les accumulateurs étant remis à zéro ensuite.