FR 3067889 A1 20181221 FR 1755343 A 20170614 Procédé de mesure de vitesse de rotation d'un arbre de turbomachine tenant compte du bruit de mesure DOMAINE TECHNIQUE DE L'INVENTION La présente invention concerne la mesure précise de régime d'un moteur d'aéronef. Il s'agit plus spécifiquement de tenir compte de l'existence de bruits parasites venant s'ajouter de manière imprévisible au signal utile, lors de la 10 mesure de la vitesse de rotation d'un arbre de moteur par capteur à roue phonique. PRESENTATION GENERALE DE L'INVENTION 15 Une détection précise de la vitesse de rotation d'un arbre de moteur d'aéronef est nécessaire, dans la mesure où cette information peut commander la coupure ou la régulation de l'alimentation du moteur, par exemple si la vitesse de rotation ou l'accélération dépassent une valeur de sécurité prédéterminée. 20 Il est bien connu de mesurer la vitesse de rotation d'un arbre à l'aide d'un capteur à reluctance variable associé à une roue phonique. On peut parler invariablement de capteur à roue phonique, ou de capteur à réluctance variable. La roue phonique est agencée de manière à être solidaire de 25 l'arbre. Le capteur crée, à l'aide d'un bobinage, un champ magnétique qui se referme soit dans l'air à l'espace entre deux dents de la roue phonique, soit sur une dent de la roue phonique. Ainsi, le flux magnétique mesuré varie en fonction du passage des dents de la roue phonique, de manière proportionnelle à la vitesse de rotation de la roue phonique. La fréquence 30 de la tension alternative ainsi générée est égale à la fréquence de passage des dents de la roue phonique, elle-même illustrative de la vitesse de 2 rotation de l'arbre, et l'amplitude de cette tension alternative est fonction de l'entrefer et de la fréquence du signal. La Déposante a décrit, dans la demande de brevet français publiée sous le numéro WO 2014/207369, un écrou denté pouvant être vissé autour d'un arbre de moteur d'aéronef, par 5 exemple pour solidariser mécaniquement cet arbre avec un roulement à bille, dont les dents peuvent coopérer avec un capteur de vitesse de rotation fonctionnant par ce principe. Au moins une partie des dents de l'écrou comportent en outre des évidements, préservant la fonction de serrage de l'écrou tout en fournissant un plus grand nombre d'intervalles de déficit dans la pièce pour permettre une mesure précise de vitesse de rotation de l'arbre. Des systèmes de mesure de vitesse de rotation à roue phonique peuvent notamment être utilisés pour mesurer le régime dit « Ni » de rotation du corps basse pression d'un turboréacteur à double corps. Sur des moteurs « direct drive » (entraînement direct de la soufflante par l'arbre basse pression), ce régime Ni est aussi le régime de la soufflante. Un capteur à roue phonique peut également permettre de mesurer le régime N2 de rotation du corps haute pression du turboréacteur double corps, ou la vitesse de rotation d'un autre élément tournant du moteur. Le signal sinusoïdal alternatif acquis aux bornes d'un capteur à reluctance variable peut être transformé, par un circuit électronique placé en amont du calculateur moteur, en signal créneau sur lequel le calculateur moteur peut réaliser une mesure de fréquence. La transformation du signal sinusoïdal en signal créneau est notamment réalisée par un trigger de Schmidt, avec un seuil de redémarrage et un seuil de réarmement prédéterminés. Le signal créneau en sortie de trigger de Schmidt commute vers une valeur haute lorsque la tension d'entrée dépasse le seuil de réarmement, et commute vers une valeur basse lorsque la tension d'entrée passe en-deçà du seuil de redémarrage. 3 Or, le signal créneau obtenu pour un capteur de régime à roue phonique peut être très perturbé à cause de bruits parasites. Des parties du capteur peuvent subir des perturbations électromagnétiques, dues à des vibrations mécaniques de certaines parties de moteur. On peut observer des 5 événements ponctuels de bruit, mais également des signaux de bruit du type alternatif se superposant au signal représentatif de la rotation de la roue phonique (signal utile). A ce titre, la figure 1 illustre, sur un premier axe, la tension Sv aux bornes du capteur à roue phonique au cours du temps, dans le cas d'un signal altéré par un bruit parasite de type alternatif 10 de fréquence supérieure à celle du signal utile, et illustre, sur un deuxième axe, le signal créneau Sc obtenu après crénelage du signal de tension Sv, en sortie de trigger. L'amplitude pic à pic du bruit pouvant être localement supérieure à la différence de tension entre le seuil de redémarrage UR et le seuil de réarmement UA, le bruit alternatif peut provoquer des commutations 15 parasites du signal créneau. Un troisième axe illustre les variations d'un signal Cp de comptage des périodes successives du signal créneau Sc au fil du temps. La temporisation illustrée par la courbe de ce troisième axe est relancée à chaque commutation descendante du signal créneau, du niveau haut vers le niveau bas. Par exemple, au cours de la période Pi illustrée sur 20 la figure, il ne se produit pas de commutation parasite due au bruit. En revanche, sur la période suivante P2, le signal Sv détecté franchit le seuil de réarmement UA, puis redescend en-deçà du seuil de redémarrage UR, du fait du bruit alternatif parasite superposé au signal utile. La durée écoulée entre ces deux passages de seuils de trigger est bien inférieure à la période 25 Pi, ce qui provoque des commutations rapides successives du signal créneau Sc, et la détection, en lieu et place d'une unique période « complète » P2 qui aurait été obtenue sans commutation parasite, de deux périodes successives P2a et P2b « morcelées » du signal créneau Sc. La mesure de vitesse de rotation de l'arbre est alors perturbée. 30 4 Une solution connue consiste à filtrer le bruit parasite en amont de l'électronique qui produit le signal crénelé. Il s'agit d'une solution bas niveau, mise en oeuvre en amont du calculateur moteur. Il est par exemple fréquent de filtrer, par un filtrage passe-bas, les bruits dont la fréquence est 5 supérieure à celle du signal utile. Cependant, le filtrage passe-bas est susceptible de déphaser la mesure obtenue. Une mesure déphasée peut s'avérer gênante pour d'autres fonctions que peut remplir le signal acquis aux bornes du capteur à reluctance variable, par exemple une fonction d'équilibrage moteur, si la position de l'arbre est prise comme référence 10 pour le contrôle d'autres éléments du moteur. De plus, le filtrage passe-bas peut réduire l'amplitude du signal. A seuils de redémarrage et de réarmement constants et à amplitude réduite du signal de mesure de vitesse de rotation, la précision de la mesure de fréquence obtenue est réduite. Ce problème est d'autant plus pertinent en bas régime moteur (par 15 exemple au démarrage du moteur), où le rapport de l'amplitude pic à pic d'un bruit parasite de type alternatif sur l'amplitude pic à pic du signal utile peut être plus important que pour un régime moteur plus élevé. Par ailleurs, il existe des solutions software, mises en oeuvre au niveau de 20 le système applicatif de l'interface numérique du moteur (par exemple, un système de type FADEC), réalisant un post-traitement de la tension sinusoïdale acquise aux bornes d'un capteur à reluctance variable. Cependant, aucune des solutions existantes n'est précisément adaptée à la problématique de traitement d'un signal bruité par un bruit alternatif de 25 haute fréquence. Notamment, des solutions de l'art antérieur qui reposent sur un simple moyennage de mesures de période acquises à des temps antérieurs, pour déterminer la période du passage des dents de la roue phonique au temps 30 courant, ne peuvent être efficaces qu'à condition d'acquérir un nombre suffisant de périodes « complètes » non bruitées. 5 PRESENTATION GENERALE DE L'INVENTION La présente invention résout les problèmes mentionnés ci-avant par un 5 post-traitement des valeurs de période détectées sur le signal crénelé, ce traitement étant implémenté au niveau de l'applicatif software de l'interface numérique du moteur, pour restituer la « période utile » (c'est-à-dire la période du signal utile qui serait obtenu sans bruit parasite) à partir d'un ensemble d'échantillons de périodes du signal crénelé enregistrés dans une 10 mémoire de type buffer circulaire. Ce post-traitement repose sur l'observation qu'en présence de commutations parasites telles que celles mises en exergue sur le signal Sc de la figure 1 présentée ci-avant, on peut obtenir, en réalisant la somme de plusieurs périodes consécutives du signal Sc, une valeur très proche de celle de la période qui serait observée pour 15 le signal utile. Le procédé de l'invention permet de retrouver les bonnes sommes à effectuer, en réalisant une pluralité de sommes d'échantillons de période consécutifs et en comparant ces sommes à des seuils de vraisemblance, calculés en fonction des mesures de période et de régime réalisées précédemment. 20 Un tel traitement ne nécessite plus de disposer d'échantillons de période correspondant à des périodes complètes du signal utile. Même si l'on dispose seulement d'échantillons de périodes « morcelées » du fait de commutations parasites obtenues sur le signal créneau, le procédé de 25 l'invention est en mesure de retrouver une période proche de celle qui serait obtenue pour le signal utile. A ce titre, l'invention vise, selon un premier aspect, un procédé de mesure de vitesse de rotation d'un arbre d'un moteur d'aéronef, comprenant des 30 étapes de : 6 acquisition d'un signal alternatif de détection de vitesse de rotation, obtenu aux bornes d'un capteur à roue phonique agencé à proximité de l'arbre, transformation dudit signal alternatif en signal créneau, le signal créneau commutant vers un niveau bas de tension, respectivement vers un niveau 5 haut de tension, en correspondance avec le passage de la tension du signal alternatif en-dessous d'un seuil de redémarrage, respectivement au-dessus d'un seuil de réarmement, enregistrement, dans une mémoire de type buffer circulaire, d'une pluralité d'échantillons de période de signal créneau, et comparaison de chacun 10 desdits échantillons à une borne inférieure de période et une borne supérieure de période, pour déterminer des échantillons dits valides dont la valeur est comprise entre ces deux bornes; si le nombre d'échantillons valides déterminés est supérieur à un premier seuil, détermination à partir desdits échantillons valides d'une période utile 15 du signal créneau, et au moins à la condition que le nombre d'échantillons valides soit inférieur au premier seuil, détermination de la période utile du signal par les sous-étapes suivantes : - calcul d'une pluralité de sommes d'au moins deux échantillons, 20 détermination parmi cette pluralité de sommes, de sommes valides d'échantillons qui sont comprises entre les bornes inférieure et supérieure de période, - calcul de la moyenne d'un ensemble constitué des échantillons valides et des sommes déterminées valides d'échantillons, la période utile 25 du signal créneau étant prise égale à ladite moyenne. De manière avantageuse, ce procédé peut présenter les caractéristiques supplémentaires suivantes : 30 - si la somme du nombre d'échantillons valides et du nombre de sommes valides d'échantillons est inférieure à un deuxième seuil, la période 7 utile du signal créneau est prise égale à une période utile du signal créneau obtenue lors d'une itération précédente du procédé ; - les bornes inférieure et supérieure de période sont calculées à partir 5 d'une valeur de période consolidée, et d'une valeur de variation maximale de période ; - la période consolidée est alors prise égale à une période utile de signal créneau obtenue à un instant antérieur ; 10 - la valeur de variation maximale de période est alors obtenue à l'aide d'une loi de gradient maximal de vitesse de rotation, où le gradient maximal à un instant donné est fonction de la vitesse de rotation ; 15 - durant une phase d'initialisation de la mesure, phase pendant laquelle le nombre d'échantillons valides de période ne dépasse pas un troisième seuil, ou alors pendant laquelle les échantillons de période détectés sur le signal créneau sont dispersés entre eux au-delà d'un pourcentage seuil prédéterminé, une période utile du signal créneau est 20 estimée par un modèle numérique théorique du régime ; - la période utile du signal créneau est alors déterminée pendant la phase d'initialisation par un modèle numérique du régime en fonction de paramètres thermodynamiques du moteur, de préférence la valeur du 25 régime de corps haute pression, la pression atmosphérique, le calage stator du compresseur. Selon un deuxième aspect, l'invention concerne un calculateur de moteur de turbomachine mettant en oeuvre un procédé de mesure de vitesse de 30 rotation d'arbre comme ci-avant. 8 L'invention concerne enfin selon un troisième aspect un produit programme d'ordinateur, destiné à être utilisé dans un système applicatif d'un calculateur de moteur de turbomachine, dont des instructions de code permettent la mise en oeuvre d'un procédé comme ci-avant. 5 PRESENTATION GENERALE DES FIGURES La figure 1 a déjà été décrite ci-avant. 10 D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, accompagnée des dessins annexés supplémentaires suivants : La figure 2 illustre, sur un exemple de signal bruité en sortie de capteur à 15 roue phonique, des détections erronées de période de signal créneau. La figure 3 représente un signal de tension théorique résultant de la superposition d'un signal utile et d'un bruit de type alternatif sinusoïdal, pour illustrer le principe de l'addition d'échantillons de période. 20 La figure 4 représente de manière schématique un mode de réalisation d'un algorithme de traitement d'échantillons de période acquis dans une mémoire circulaire pour recalculer la période utile. 25 La figure 5 représente un mode de réalisation d'un algorithme de détermination d'intervalle de vraisemblance au sein de l'algorithme de la figure 4. La figure 6 illustre les résultats obtenus, avec l'algorithme illustré en figure 30 4, sur un signal de vitesse de rotation d'arbre, à partir d'une phase 9 d'initialisation de la mesure et sur une durée comprenant une plage temporelle restreinte Z. La figure 7 représente à échelle de temps plus réduite, sur la plage Z, les 5 signaux de la figure 6. DESCRIPTION DETAILLEE D'UN MODE DE REALISATION Dans toute la description ci-après, on acquiert un signal Sv alternatif aux 10 bornes d'un capteur à roue phonique, de détection de vitesse de rotation de l'arbre du corps basse pression (BP) d'un turboréacteur double corps. En aval du capteur, et en amont d'un Electronic Engine Controller (EEC) tel qu'un FADEC3, le signal Sv est transformé en un signal créneau Sc. Pour ce faire, le signal Sv est ici filtré par un filtre RC, puis écrêté, puis à nouveau 15 filtré, et enfin crénelé par un trigger de Schmidt de seuil de redémarrage UR à O V et de seuil de réarmement UA à 0.232 V. Des périodes de passage de dents de roue phonique sont mesurées sur le signal créneau Sc, en considérant les durées écoulées entre deux fronts 20 descendants du signal créneau. Le capteur à roue phonique mesure ici le régime Ni d'un moteur. Toutefois, les mêmes modalités de traitement peuvent tout aussi avantageusement être utilisées pour la mesure de tout autre régime d'élément tournant du moteur, comme par exemple pour la mesure du régime N2. 25 Or, comme vu ci-avant en relation à la figure 1, le signal créneau Sc peut présenter des commutations qui ne correspondent pas à des fronts montants ou descendants pour la partie utile du signal de détection de régime, mais qui sont au contraire provoquées par une variation de tension 30 d'un bruit alternatif parasite (partie non utile du signal). Ce phénomène est illustré en figure 2. Sur un premier axe, une tension théorique Sv aux 10 bornes d'un capteur à roue phonique est représentée en fonction du temps t. On a aussi représenté sur le même axe la partie utile Su du signal, obtenue de manière théorique en soustrayant au signal Sv la contribution du bruit alternatif parasite. Sur un deuxième axe temporel, on a superposé un signal 5 créneau associé au signal Sv, et un signal créneau associé au signal utile théorique Su, commutant entre un niveau bas Ui et un niveau haut U2. Il y a commutation de ces signaux créneaux vers le niveau haut U2 lorsque la tension du signal Sv passe au-delà du seuil UA, et commutation de ces signaux créneaux vers le niveau bas Ui lorsque la tension du signal Sv 10 passe en-deçà du seuil UR. On observe, sur le deuxième axe, que le signal bruité Sv donne après crénelage un front montant et un front descendant correspondant à la partie 20', puis un front montant et un front descendant correspondant à la partie 15 20', tandis que le signal utile Su donne après crénelage un seul front montant et un seul front descendant correspondant à la partie 20. En outre, les fronts montant et descendant de la partie 20' ne sont pas en correspondance avec les fronts montant et descendant de la partie 20. Ainsi, la détection de période (par exemple la période entre deux fronts 20 descendants successifs de signal créneau) est erronée pour le signal bruité Sv, par rapport à la période obtenue pour le signal utile Su. Or, c'est la période du signal utile Su qui est représentative du régime Ni utile. Dans tout ce qui va suivre, on décrit un procédé de traitement 10 appliqué 25 à une pluralité de valeurs T1,... ,T de période acquises sur le signal créneau Sc. Le procédé 10 est mis en oeuvre par un algorithme implémenté dans le système applicatif du calculateur électronique de régulation du moteur (EEC). Il permet de reconstituer une période consolidée Pt à partir des valeurs ,Tn. Pour rappel, une période de signal créneau Sc est ici 30 obtenue entre deux fronts descendants consécutifs (du niveau de tension haut U2 vers le niveau de tension bas Ui) de ce signal créneau. Il n'est pas 11 nécessaire de supposer à l'avance que ces périodes sont « complètes », c'est-à-dire associées au signal utile théorique de détection de régime hors bruit parasite, ou « morcelées », c'est-à-dire calculées entre des commutations parasites de signal Sc liées au bruit. Le procédé 10 permet 5 de reconstituer une t-ième période Pt réalisant une bonne approximation de la période du signal utile théorique Su. Le procédé 10 est ici mis en oeuvre de manière numérique, au niveau du logiciel applicatif, ou système applicatif, de l'interface numérique EEC d'un moteur d'aéronef. 10 La figure 3 permet d'éclairer, sur un exemple théorique simplifié, les objectifs et le principe du procédé 10. Sur un premier axe, on a représenté au cours du temps un exemple de signal de tension Sv, obtenu en sommant un signal utile théorique Su sinusoïdal de basse fréquence, également représenté sur ce premier axe, et un signal de bruit de plus haute fréquence. 15 On a également représenté le seuil de redémarrage UR et le seuil de réarmement UA d'un trigger de Schmidt utilisé pour créneler le signal Sv et le signal théorique Su, ces derniers signaux correspondant respectivement au signal créneau Sc, représenté en parallèle sur un deuxième axe, et au signal créneau théorique Sc-r, représenté en parallèle sur un troisième axe. 20 On constate que la période Pt qui serait obtenue pour le signal ScT se superpose partiellement à une première période Ti effectivement obtenue pour le signal Sc, puis se superpose complètement à deux périodes suivantes T2 et T3, puis se superpose partiellement à une quatrième période non référencée sur la figure 3. 25 Or, en additionnant les périodes consécutives Ti, T2 et T3, on obtient une valeur de période très proche de la valeur Pt, cette dernière constituant la période utile pour la mesure de fréquence. Une interprétation de ce résultat est que les deux fronts descendants de Sc marquant respectivement le 30 début de la période T1 et la fin de la période T3 correspondent approximativement à des fronts descendants consécutifs qui seraient 12 observés pour le signal Sc-r, donnant la période Pt. Si l'on effectue une analogie entre cet exemple théorique et un cas réel de mesure de signal de détection acquis aux bornes d'un capteur à roue phonique, le signal effectivement observable aux bornes du capteur serait le signal Sv, donnant 5 le signal créneau Sc. Un objectif du procédé 10 est alors d'obtenir une bonne approximation de la période Pt obtenue avec le signal créneau ScT, la période Pt n'étant pas directement observable. On décrit ci-après les étapes du procédé 10 de traitement de valeurs de 10 période d'un signal créneau Sc acquis à partir d'un signal Sv de détection de vitesse de rotation pris aux bornes d'un capteur à roue phonique, le procédé 10 étant ici, pour rappel, mis en oeuvre par le logiciel applicatif, ou système applicatif, de l'interface numérique du calculateur électronique de régulation du moteur. Les étapes du procédé 10 sont ci-après décrites en 15 relation à la figure 4. A une étape 100, l'acquisition des valeurs d'échantillons de période de signal créneau s'effectue de la manière suivante. Un premier compteur détecte les fronts descendants d'un signal créneau d'échantillonnage de 20 fréquence élevée par rapport à celle attendue pour le signal créneau. Par exemple, la carte de mesure des périodes de passage de dents de roue phonique peut avoir une période d'échantillonnage de l'ordre de 0.25.10-6 s. Un deuxième compteur détecte les fronts descendants du signal créneau Sc. Les valeurs de durée écoulée entre deux fronts descendants 25 consécutifs du signal créneau Sc peuvent ainsi être enregistrées dans une mémoire tournante B de type buffer circulaire, pouvant conserver en mémoire un certain nombre de mesures de période. Il est ici préférable d'implémenter le buffer B dans le circuit électronique de mesure de la période, au lieu d'une implémentation numérique dans le système applicatif 30 de l'interface numérique du calculateur électronique du moteur, car la période de calcul de ce calculateur moteur est typiquement de l'ordre des 13 0.015 s et ne permet pas la mesure précise de la période de passage des dents de la roue phonique (dont la vitesse peut atteindre les 7000 tours par minute, à raison de 60 dents). Le buffer circulaire B peut ensuite fournir les valeurs Ti,... ,Tri, ainsi que les dates de mesure de ces échantillons, au 5 système applicatif. Le nombre, ici noté n, d'échantillons est pris supérieur à un nombre minimal d'échantillons nécessaire pour assurer une bonne précision de la mesure. Le nombre minimal peut correspondre, par exemple, au nombre de dents 10 de la roue phonique. En effectuant, à des étapes ultérieures reprises ci- après du procédé 10, une moyenne d'un nombre d'échantillons consécutifs supérieur au nombre de dents de la roue phonique, on gomme d'éventuelles disparités entre les mesures qui seraient liées à des inhomogénéités de géométrie des dents de roue phonique. La fréquence 15 attendue des bruits parasites dont on cherche à réduire l'influence intervient également dans la détermination du nombre minimal d'échantillons du buffer B: plus le bruit est de haute fréquence par rapport au signal utile, plus de périodes « morcelées » pourront être nécessaires pour reconstituer, selon le principe décrit en relation à la figure 3, une période utile de passage 20 de dent de la roue phonique. A une étape 110, les échantillons sont ensuite triés dans l'ordre du plus ancien échantillon au plus récent. Cette étape 110 permet, aux étapes ultérieures, de réaliser des sommes de périodes consécutives dans le 25 temps. On suppose ici, pour plus de commodité, que les échantillons triés conservent leur ordre de Ti à Tn après le tri. En sortie d'étape 110, le système applicatif dispose donc d'un vecteur constitué de n échantillons de périodes consécutives. 14 A la suite de l'étape 110, dans le présent mode de réalisation, on distingue deux cas de figure, en fonction de la valeur d'un signal binaire d'initialisation qui indique si la mesure de période est à un état initialisé ou non. En effet, comme il sera décrit ci-après, aux étapes suivantes de 5 reconstitution de période consolidée à partir de l'étape 120, une valeur de période consolidée obtenue à un pas de calcul précédent est prise en entrée. Ces étapes nécessitent donc qu'au moins une valeur de période consolidée jugée satisfaisante ait été acquise. 10 Dans le présent mode de réalisation, on effectue un test sur les échantillons de période triés dans le buffer B. En alternative, on pourrait également effectuer un test pour l'initialisation sur des échantillons de période non triés. On teste si un certain nombre d'échantillons parmi les échantillons (par exemple, dix échantillons), ou en alternative tous 15 les échantillons, sont suffisamment cohérents entre eux, c'est-à-dire si leurs valeurs sont faiblement dispersées par rapport à un pourcentage seuil de dispersion D prédéterminé. Par exemple, si le test d'initialisation nécessite que trois échantillons de 20 période soient peu dispersés entre eux, et pour une valeur du seuil de dispersion D égale à 5%, on vérifie si l'équation logique 100* P1-P2 REVENDICATIONS1. Procédé de mesure de vitesse de rotation d'un arbre d'un moteur d'aéronef, comprenant les étapes suivantes : - acquisition d'un signal alternatif (Sv) de détection de vitesse de rotation, obtenu aux bornes d'un capteur à roue phonique agencé à proximité de l'arbre ; - transformation dudit signal alternatif (Sv) en signal créneau (Sc), le signal créneau (Sc) commutant vers un niveau bas de tension (1.11), respectivement vers un niveau haut de tension (U2), en correspondance avec le passage de la tension du signal alternatif (Sv) en-dessous d'un seuil de redémarrage (UR), respectivement au-dessus d'un seuil de réarmement (UA); - enregistrement, dans une mémoire de type buffer circulaire (B), d'une pluralité d'échantillons ,T,i) de période de signal créneau (Sc), - comparaison de chacun desdits échantillons à une borne inférieure de période (Pt min) et une borne supérieure de période (Pt max), pour déterminer des échantillons dits valides dont la valeur est comprise entre ces deux bornes; - si le nombre d'échantillons valides déterminés est supérieur à un premier seuil (Ai), détermination à partir desdits échantillons valides d'une période utile (Pt) du signal créneau (Sc), caractérisé en ce qu'au moins à la condition que le nombre d'échantillons valides soit inférieur au premier seuil (Ai) la détermination de la période utile (Pt) du signal créneau (Sc) met en oeuvre les étapes suivantes : - calcul d'une pluralité de sommes d'au moins deux échantillons ,Sin,S23,... ,S(n-i)n) - détermination parmi cette pluralité de sommes, de sommes valides d'échantillons qui sont comprises entre les bornes inférieure et supérieure de période (P t min, Pt max), 25 calcul de la moyenne d'un ensemble constitué des échantillons valides et des sommes déterminées valides d'échantillons, la période utile (Pt) du signal créneau étant prise égale à ladite moyenne. 5 2. Procédé de mesure de vitesse selon la revendication 1, dans lequel si la somme du nombre d'échantillons valides et du nombre de sommes valides d'échantillons est inférieure à un deuxième seuil (A2), la période utile (Pt) du signal créneau est prise égale à une période utile du signal créneau 10 (Pt_i) obtenue lors d'une itération précédente du procédé. 3. Procédé de mesure de vitesse selon l'une des revendications précédentes, dans lequel les bornes inférieure et supérieure (Pt min, Pt max) de période sont calculées à partir d'une valeur de période consolidée, et 15 d'une valeur de variation maximale de période. 4. Procédé de mesure selon la revendication 3, dans lequel la période consolidée est prise égale à une période utile de signal créneau obtenue à un instant antérieur (Pt-i). 20 5. Procédé de mesure selon l'une des revendications 3 ou 4, dans lequel la valeur de variation maximale de période est obtenue à l'aide d'une loi de gradient maximal de vitesse de rotation (GNimax), où le gradient maximal à un instant donné est fonction de la vitesse de rotation. 25 6. Procédé de mesure selon l'une quelconque des revendications précédentes, dans lequel durant une phase d'initialisation de la mesure, phase pendant laquelle le nombre d'échantillons valides de période ne dépasse pas un troisième seuil (A3), une période utile du signal créneau est 30 estimée par un modèle numérique théorique de période (Pth). 26 7. Procédé de mesure selon l'une quelconque des revendications précédentes, dans lequel durant une phase d'initialisation de la mesure, pendant laquelle les échantillons de période détectés sur le signal créneau (Sc) sont dispersés entre eux au-delà d'un seuil prédéterminé (D), une 5 période utile du signal créneau est estimée par un modèle numérique théorique de période (Pth). 8. Procédé de mesure selon la revendication 6 ou 7, où la période utile du signal créneau est déterminée pendant la phase d'initialisation par un 10 modèle numérique en fonction de paramètres thermodynamiques du moteur, de préférence pris parmi la valeur du régime de corps haute pression ou la valeur du régime d'un autre élément tournant du moteur, la pression atmosphérique, le calage stator du compresseur. 15 9. Calculateur d'un moteur d'aéronef mettant en oeuvre un procédé de mesure de vitesse de rotation d'arbre selon l'une quelconque des revendications précédentes. 10.Produit programme d'ordinateur, destiné à être utilisé dans un 20 système applicatif d'un calculateur de moteur de turbomachine, dont des instructions de code permettent la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 8.