FR 3080453 A1 20191025 FR 1853544 A 20180423 Arrière-plan de l'invention L'invention concerne le domaine général de la mécanique. Elle se rapporte plus particulièrement au contrôle non destructif 5 par ultrasons de pièces mécaniques comme par exemple des pièces aéronautiques telles que des billettes en alliage de titane, des disques de soufflantes ou des disques aubagés monoblocs équipant des moteurs d'aéronefs et notamment des turboréacteurs. Aucune limitation n'est attachée toutefois à la nature ni à la forme des pièces mécaniques 10 considérées. De façon connue, on désigne par contrôle non destructif un ensemble de méthodes qui permettent de caractériser l'état d'intégrité et/ou la qualité de structures ou de matériaux sans les dégrader. Le contrôle non destructif a une application privilégiée mais non limitative 15 dans le domaine de l'aéronautique, et plus généralement dans n'importe quel domaine industriel dans lequel les structures dont on souhaite caractériser l'état ou la qualité sont coûteuses et/ou leur fiabilité de fonctionnement est critique. Le contrôle non destructif peut être avantageusement réalisé sur la structure ou le matériau considéré aussi 20 bien en cours de production qu'en cours d'utilisation ou de maintenance. Un certain nombre de pièces mécaniques critiques équipant les moteurs d'aéronef et plus particulièrement les turboréacteurs, sont réalisées en alliage de titane à partir de billettes cylindriques. Afin de s'assurer de l'intégrité de telles pièces mécaniques avant leur montage 25 dans un turboréacteur, il est connu de réaliser sur ces pièces un contrôle non destructif utilisant la propagation d'ondes ultrasonores, et permettant de détecter la présence éventuelle de défauts ou d'indications de défauts dans les pièces contrôlées. Ce contrôle non destructif est généralement réalisé sur les billettes et après forgeage. Il est classique d'utiliser à cet 30 effet un transducteur (i.e. un capteur) composé de multiples éléments piézoélectriques aptes à émettre et à recevoir des ondes ultrasonores. Lors du contrôle non destructif, des ondes ultrasonores sont émises par le transducteur multiélément à une fréquence déterminée, et interagissent avec la microstructure du matériau de la pièce dont on 35 souhaite vérifier l'intégrité. Un « bruit » ultrasonore est mesuré par le transducteur multiélément : ce bruit est en concurrence avec les échos 2 associés aux indications de défauts potentiellement présentes dans la pièce. Plus spécifiquement, chaque élément du transducteur multiélément mesure un signal comprenant plusieurs composantes, à savoir : une composante dite de diffusion simple, représentant les trajets d'onde où l'onde propagée dans la pièce mécanique n'interagit qu'une seule fois avec le matériau de la pièce (i.e. n'est réfléchie qu'une seule fois sur les diffuseurs présents dans le matériau de la pièce, ces diffuseurs pouvant être liés à des défauts présents dans la pièce ou aux différentes faces de la pièce) avant d'être mesurée par les éléments du transducteur multiélément (on parle aussi d'onde rétrodiffusée reçue par les éléments du transducteur multiélément) ; éventuellement une composante dite de diffusion multiple, représentant les trajets d'onde où l'onde propagée dans la pièce mécanique interagit plusieurs fois avec le matériau de la pièce mécanique (i.e. fait l'objet de plusieurs réflexions successives sur les diffuseurs présents dans le matériau) avant d'être mesurée par les éléments du transducteur multiélément ; et une composante de bruit électronique. L'utilisation d'un transducteur multiélément offre avantageusement la possibilité d'appliquer aux signaux mesurés par le transducteur, un filtrage numérique (matriciel) permettant de réduire la composante de bruit électronique et le cas échéant, la composante de diffusion multiple, présente(s) dans ces signaux au profit des composantes reflétant le cas échéant la présence d'indications de défauts dans la pièce mécanique. Ceci permet d'améliorer l'analyse des signaux mesurés observés pour contrôler l'intégrité de la pièce mécanique, et le cas échéant la détection de défauts dans la pièce mécanique. Le document WO 2010/001027 décrit une méthode de filtrage numérique basée sur la séparation de la composante de diffusion simple de la composante de diffusion multiple. Cette méthode consiste à projeter les signaux mesurés par le transducteur multiélément sur un sous-espace caractéristique du régime de diffusion simple et estimé de façon théorique à partir des propriétés de propagation du matériau de la pièce contrôlée (vitesse de propagation de l'onde dans le matériau et nombre d'onde notamment). Cette estimation théorique s'appuie sur l'hypothèse d'une interface plane entre le transducteur multiélément et la pièce mécanique, 3 le transducteur multiélément étant placé en contact direct avec la pièce mécanique. En pratique, les contrôles non destructifs réalisés sur les pièces aéronautiques, et notamment sur les billettes mentionnées 5 précédemment, sont souvent réalisés en immersion afin de permettre une inspection rapide du volume entier des pièces. Dans ce contexte, le transducteur multiélément est déplacé autour de la pièce mécanique immergée et il n'est pas toujours possible de le positionner en contact direct avec la pièce mécanique. En outre, les pièces mécaniques dont on 10 souhaite vérifier l'intégrité ne sont pas toujours planes : elles peuvent au contraire avoir des formes plus complexes, comme par exemple des formes cylindriques, sphériques ou toutes autres, qui ne correspondent pas aux hypothèses du document WO 2010/001027. Il existe donc un besoin d'une technique de filtrage numérique 15 pouvant être utilisée pour le contrôle non destructif par ultrasons d'une pièce mécanique, quelle que soit la forme de cette pièce mécanique et de l'interface entre la pièce mécanique et le transducteur multiélément apte à émettre les ondes ultrasonores. 20 Objet et résumé de l'invention La présente invention répond notamment à ce besoin en proposant un procédé de contrôle non destructif d'une pièce mécanique par propagation d'ondes ultrasonores réalisée au moyen d'un transducteur multiélément placé en vis-à-vis de la pièce mécanique et comprenant une 25 pluralité d'éléments piézoélectriques e(1),e(2),...,e(N), N désignant un entier supérieur à 1, ce procédé de contrôle comprenant : - une étape de mesure comprenant, pour chaque élément e(i) du transducteur, l'émission d'une onde ultrasonore à une fréquence donnée et la mesure, par chaque élément e(j) du transducteur distinct 30 de l'élément e(i), d'un signal temporel noté kij(t) représentatif de l'onde ultrasonore rétrodiffusée reçue par l'élément e(j) ; - une étape de détermination d'au moins une première matrice de composantes temporelles à partir des signaux kij(t) mesurés, i,j=1,...,N, pris sur une durée prédéterminée autour d'une référence 35 temporelle notée T; 4 - une étape de détermination d'au moins une deuxième matrice de composantes fréquentielles correspondant à une fréquence déterminée à partir de la fréquence de l'onde sonore en appliquant une transformée de Fourier aux composantes temporelles de ladite au 5 moins une première matrice ; - une étape de filtrage de ladite au moins une deuxième matrice comprenant une projection de ladite au moins une deuxième matrice de composantes fréquentielles sur un sous-espace de diffusion simple ; et 10 - une étape de vérification de l'intégrité de la pièce mécanique en utilisant ladite au moins une deuxième matrice filtrée. Le procédé de contrôle est remarquable en ce que le sous-espace de diffusion simple utilisé lors de l'étape de filtrage est déterminé au moyen d'un calcul numérique utilisant un algorithme de tracé de rayons. 15 Par exemple, dans un mode particulier de réalisation, l'étape de filtrage de la deuxième matrice comprend : une étape de détermination du sous-espace de diffusion simple, ledit sous-espace de diffusion simple étant défini à partir d'une base de 2N1 matrices notées Em, m=1,...,2N-1, ladite étape de détermination 20 comprenant : o pour chaque couple d'éléments distincts e(i) et e(j) du transducteur, i, j, =1,...,N, une étape de calcul numérique en utilisant un algorithme de tracé de rayons, pour une pluralité D de diffuseurs notés Pd, d=1...,D répartis en divers points de la 25 pièce mécanique et situés à une profondeur déterminée en fonction de la référence temporelle T, d'un temps de propagation tijd d'une onde ultrasonore émise par l'élément e(i), diffractée par le diffuseur Pd et reçue par l'élément e(j), le temps de propagation tijd étant calculé en tenant compte de la 30 forme de la pièce mécanique, du placement du transducteur et des milieux dans lesquels se propage l'onde ultrasonore entre les éléments e(i) et e(j) du transducteur et le diffuseur Pd ; o une étape d'estimation d'une contribution de la diffusion simple dans les composantes fréquentielles de la deuxième matrice à 35 partir des temps de propagation calculés pour chaque couple 5 d'éléments distincts e(i) et e(j) du transducteur, ladite étape d'estimation fournissant une troisième matrice ; o une étape d'extraction des composantes fréquentielles d'une anti-diagonale centrale de la troisième matrice pour construire la 5 base de matrices Em, m=1,...,2N-1 définissant le sous-espace de diffusion simple ; - une étape de projection de ladite au moins une deuxième matrice sur chacune des matrices Em, m=1,..,2N-1 de la base de 2N-1 matrices ; et 10 - une étape d'évaluation de ladite au moins une deuxième matrice filtrée à partir des matrices Em de la base pondérées par les projections obtenues lors de l'étape de projection. Corrélativement, l'invention vise aussi un système de contrôle non destructif d'une pièce mécanique comprenant : 15 - un transducteur multiélément placé en vis-à-vis de la pièce mécanique et apte à émettre des ondes ultrasonores, ledit transducteur comprenant une pluralité d'éléments piézoélectriques e(1),e(2),...,e(N), N désignant un entier supérieur à 1, ledit transducteur étant configuré pour émettre, via chaque élément e(i) du transducteur, une onde 20 ultrasonore à une fréquence donnée, et à mesurer, via chaque élément e(j) du transducteur distinct de l'élément e(i), un signal temporel noté kij(t) représentatif de l'onde ultrasonore rétrodiffusée reçue par l'élément e(j) ; - un dispositif de validation comprenant : 25 o un premier module de détermination, configuré pour déterminer au moins une première matrice de composantes temporelles à partir des signaux kij(t) mesurés, i,j=1,...,N, pris sur une durée prédéterminée autour d'une référence temporelle notée T; o un second module de détermination, configuré pour déterminer 30 au moins une deuxième matrice de composantes fréquentielles correspondant à une fréquence déterminée à partir de la fréquence de l'onde sonore en appliquant une transformée de Fourier aux composantes temporelles de ladite au moins une première matrice ; 35 o un module de filtrage, configuré pour filtrer ladite au moins une deuxième matrice en projetant ladite au moins une deuxième 6 matrice de composantes fréquentielles sur un sous-espace de diffusion simple ; et o un module d'analyse, configuré pour vérifier l'intégrité de la pièce mécanique en utilisant ladite au moins une deuxième 5 matrice filtrée ; le système de contrôle étant remarquable en ce que le module de filtrage est configuré pour déterminer le sous-espace de diffusion simple au moyen d'un calcul numérique utilisant un algorithme de tracé de rayons. Dans un mode particulier de réalisation, le module de filtrage 10 comprend notamment : - un troisième module de détermination, configuré pour déterminer le sous-espace de diffusion simple à partir d'une base de 2N-1 matrices notées Em, m=1,...,2N-1, ledit troisième module de détermination étant configuré pour: 15 o pour chaque couple d'éléments distincts e(i) et e(j) du transducteur, j, j, =1,...,N, calculer numériquement en utilisant un algorithme de tracé de rayons, pour une pluralité D de diffuseurs notés Pd, d=1...,D répartis en divers points de la pièce mécanique et situés à une profondeur déterminée en 20 fonction de la référence temporelle T, un temps de propagation tijd d'une onde ultrasonore émise par l'élément e(i), diffractée par le diffuseur Pd et rétrodiffusée par l'élément e(j), ledit temps de propagation étant calculé en tenant compte de la forme de la pièce mécanique, du placement du transducteur et 25 des milieux dans lesquels se propage l'onde ultrasonore entre les éléments e(i) et e(j) du transducteur et le diffuseur Pd ; o estimer une contribution de la diffusion simple dans les composantes fréquentielles de la deuxième matrice à partir des temps de propagation calculés pour chaque couple d'éléments 30 distincts e(i) et e(j) du transducteur, et obtenir une troisième matrice à partir des contributions estimées ; o extraire des composantes fréquentielles d'une anti-diagonale centrale de la troisième matrice pour construire la base de matrices Em, m=1,...,2N-1 définissant le sous-espace de 35 diffusion simple ; 7 - un module de projection, configuré pour projeter ladite au moins une deuxième matrice sur chacune des matrices Em, rn---1,..,2N-1 de la base de 2N-1 matrices ; et - un module d'évaluation configuré pour évaluer ladite au moins une 5 deuxième matrice filtrée à partir des matrices Em de la base pondérées par les projections obtenues lors de l'étape de projection. L'invention propose un filtrage numérique dans le domaine fréquentiel des signaux mesurés par le transducteur multiélément qui ne s'appuie plus, contrairement à la méthode proposée dans le document WO 10 2010/001027, sur une estimation théorique de ce sous-espace de diffusion simple basée sur un certain nombre d'hypothèses concernant l'interface entre la pièce mécanique et le transducteur multiélément, mais sur un sous-espace de diffusion simple calculé numériquement. Ce calcul numérique utilise un algorithme de tracé (ou de lancer) de rayons (« ray 15 tracing » en anglais), qui permet avantageusement de simuler de façon réaliste la propagation des ondes ultrasonores émises par le transducteur multiélément dans la pièce mécanique. Plus précisément, l'algorithme de tracé de rayons utilisé propose de relier les éléments du transducteur multiélément avec les diffuseurs 20 appartenant à la pièce mécanique par des chemins compatibles avec les lois de la physique et notamment avec les lois de Snell Descartes. Ceci permet de bien caractériser les phénomènes physiques qui se produisent lors de l'émission d'ondes ultrasonores par le transducteur multiélément à destination de la pièce mécanique que l'on cherche à contrôler (ex. 25 pénétration dans la pièce mécanique, réflexion(s) par les diffuseurs appartenant au matériau de la pièce mécanique, puis retour aux éléments du transducteur). En outre, les chemins simulés au moyen de l'algorithme de tracé de rayons tiennent compte de la forme de la pièce physique (au niveau des points d'entrée et de sortie des ondes ultrasonores dans la 30 pièce mécanique). L'utilisation d'un tel algorithme permet donc de tenir compte des conditions dans lesquelles le contrôle non destructif est opéré et par là-même, de pouvoir relâcher les contraintes sur ces conditions, par exemple au niveau de la forme de la pièce mécanique dont on souhaite vérifier 35 l'intégralité, du placement du transducteur multiélément par rapport à la pièce mécanique, des milieux dans lesquels se propage l'onde ultrasonore 8 entre les éléments du transducteur multiélément et les diffuseurs appartenant à la pièce mécanique, etc. Elle conduit à un calcul numérique précis, dont il résulte une amélioration de la qualité du contrôle non destructif opéré et le cas échéant, de la détection de défauts dans la pièce 5 mécanique. L'invention propose donc un filtrage numérique des signaux mesurés par un transducteur multiélément lors d'un contrôle non destructif d'une pièce mécanique qui peut être utilisé quelle que soit la forme de la pièce mécanique et quelle que soit l'interface entre le 10 transducteur multiélément et la pièce mécanique. Il suffit en effet pour utiliser l'invention que le transducteur multiélément soit positionné en vis-à-vis de la pièce mécanique, c'est-à-dire en face de la pièce mécanique (en contact direct ou non) afin que les ondes ultrasonores émises par le transducteur puissent se propager dans la pièce. Le transducteur peut par 15 ailleurs être déplacé autour de la pièce afin de permettre l'inspection de tout le volume de la pièce et ainsi faciliter son contrôle. Le procédé de contrôle non destructif proposé par l'invention n'est donc pas tributaire d'un contact direct entre le transducteur multiélément et la pièce mécanique, ni d'une forme plane de cette dernière. Il peut être aisément 20 utilisé pour le contrôle non destructif de pièces complexes telles que des pièces aéronautiques et en particulier des billettes, des disques aubagés ou des disques de soufflantes tels qu'évoqués précédemment, y compris lorsque ces pièces sont placées en immersion. Dans un mode particulier de réalisation de l'invention, chaque 25 temps de propagation tijd est calculé selon : tijd = minpe tid(Pe) + minps tdj (Ps) où Pe désigne un point d'entrée de coordonnées (xe,ze) de l'onde ultrasonore dans la pièce mécanique et Ps un point de sortie de coordonnées (xs,zs) de l'onde ultrasonore diffractée par le point Xd de coordonnées (xd,zd), avec : tid(x = + ,i(xe - xi)2 + (ze - zi)2 . N I (xd - xe)2 + (zd - ze)2 ) cl c2 30 et \ tdj(x) = I + (xd - xs)2 + (zd - zs)2 -,/ (xs - xj)2 + (zs - zj)2 c2 cl 9 où cl désigne la vitesse de propagation de l'onde ultrasonore dans le milieu séparant le transducteur de la pièce mécanique et c2 désigne la vitesse de propagation de l'onde ultrasonore dans la pièce mécanique. Ce mode de calcul du temps de propagation tient compte à la 5 fois de la forme du transducteur (linéaire, matricielle ou autre), de la nature de l'interface entre le transducteur et la pièce mécanique (configuration de placement par exemple du transducteur par rapport à la pièce mécanique), et de la forme de la pièce mécanique. Dans un mode particulier de réalisation, ladite au moins une 10 deuxième matrice est déterminée à la fréquence de l'onde ultrasonore. En variante, plusieurs deuxièmes matrices sont déterminées à une pluralité de fréquences contenues dans un intervalle prédéterminé défini autour de la fréquence de l'onde ultrasonore. Par ailleurs, dans un autre mode particulier de réalisation, 15 plusieurs premières matrices correspondant à des références temporelles différentes et/ou à des durées prédéterminées différentes sont déterminées. Ces différents modes de réalisation permettent d'améliorer la précision du contrôle non destructif. Les références temporelles 20 considérées correspondent typiquement à des profondeurs dans la pièce mécanique où l'on soupçonne la présence de défauts. Si une telle profondeur est connue, il suffit de se focaliser sur une référence temporelle correspondant au temps de propagation des ondes ultrasonores jusqu'à cette profondeur. Si on ne dispose pas d'indications 25 sur la localisation de défauts éventuels, plusieurs références temporelles et durées prédéterminées autour de ces références temporelles peuvent être envisagées pour améliorer le contrôle effectué. On note qu'aucune limitation n'est attachée à la technique utilisée pour vérifier l'intégrité à proprement parler de la pièce mécanique 30 à partir de la ou des deuxièmes matrices filtrées conformément à l'invention. Par exemple, l'étape de vérification de l'intégrité de la pièce mécanique peut utiliser une technique DORT de diagonalisation de l'opérateur de retournement temporel ou une méthode de focalisation totale, connues en soi. La contribution de l'invention réside davantage sur 35 le filtrage des signaux mesurés par le transducteur multiélément : ce filtrage est amélioré, ce qui permet de diminuer le bruit présent dans les 10 signaux utilisés pour le contrôle non destructif à proprement parler, quelle que soit la méthode utilisée à cet effet. Dans un mode particulier de réalisation, le système de contrôle selon l'invention comprend en outre un bras robotisé sur lequel est placé 5 le transducteur multiélément, ledit bras robotisé étant configuré pour permettre un déplacement du transducteur multiélément autour de la pièce mécanique. Comme mentionné précédemment, grâce à la mobilité du transducteur multiélément autour de la pièce mécanique, il est possible 10 d'observer aisément l'ensemble du volume de la pièce et de proposer un contrôle non destructif plus efficace et plus rapide. Il convient de noter que l'invention s'applique à différents types de transducteur multiélément. Celui-ci peut être linéaire, mais également matriciel ou avoir une forme plus complexe. 15 Brève description des dessins D'autres caractéristiques et avantages de la présente invention ressortiront de la description faite ci-dessous, en référence aux dessins annexés qui en illustrent un exemple de réalisation dépourvu de tout 20 caractère limitatif. Sur les figures : - la figure 1 représente, dans son environnement, un système de contrôle non destructif conforme à l'invention, dans un mode particulier de réalisation ; - la figure 2 représente un exemple d'architecture matérielle d'un 25 dispositif de validation appartenant au système de contrôle non destructif de la figure 1; - la figure 3 représente, sous forme d'ordinogramme, les principales étapes d'un procédé de contrôle non destructif selon l'invention tel qu'il est mis en oeuvre par le système de la figure 1; 30 - la figure 4 représente les principales étapes mises en oeuvre lors du filtrage opéré par le dispositif de validation de la figure 2 conformément à l'invention ; - la figure 5 illustre les propriétés de la matrice inter-élément mesurée par le système de contrôle non destructif de la figure 1 lors de la mise 35 en oeuvre du procédé de contrôle selon l'invention ; et 11 - la figure 6 représente un chemin considéré par l'algorithme de tracé de rayons utilisé lors du procédé de contrôle conformément à l'invention. Description détaillée de l'invention 5 La figure 1 représente, dans son environnement, un système de contrôle non destructif 1 conforme à l'invention, dans un mode particulier de réalisation. Dans l'exemple envisagé à la figure 1, le système de contrôle non destructif 1 est utilisé à titre illustratif pour réaliser le contrôle de 10 l'intégrité d'une pièce mécanique 2 de forme cylindrique. Cette pièce mécanique est par exemple ici une billette en alliage de titane telle qu'évoquée précédemment, et utilisée pour la conception de pièces aéronautiques destinées à équiper un moteur d'aéronef tel qu'un turboréacteur. 15 Bien entendu cet exemple n'est donné qu'à titre illustratif et l'invention s'applique à d'autres pièces mécaniques, quelle que soit leur destination et de forme quelconque (par exemple à un disque aubagé monobloc ou à un disque de soufflante d'un moteur d'aéronef). Conformément à l'invention, le système de contrôle non 20 destructif 1 comprend : - un transducteur multiélément 3 (linéaire ici), comprenant une pluralité d'éléments piézoélectriques e(1),e(2),...,e(N), N désignant un entier supérieur à 1. Chaque élément piézoélectrique e(i) du transducteur est apte à émettre une onde ultrasonore à une fréquence donnée ; chaque 25 élément piézoélectrique e(j), jei, du transducteur est apte à mesurer un signal temporel noté kij(t) représentatif de l'onde ultrasonore émise par l'élément piézoélectrique e(i) et rétrodiffusée, reçue par l'élément eu). Dans le mode de réalisation décrit ici, le transducteur multiélément 3 est piloté par un amplificateur multivoie 4, connu en 30 soi, et non décrit en détail ici. Celui-ci est apte à générer et à amplifier les N signaux (ondes ultrasonores) destinés à être appliqués aux éléments piézoélectriques, et à les appliquer de façon indépendante sur chacun des éléments piézoélectriques ; il est également apte à amplifier les signaux mesurés par chaque élément piézoélectrique ; 35 - un bras robotisé 5 sur lequel est positionné le transducteur multiélément et qui permet de placer le transducteur multiélément en 12 vis-à-vis de la pièce mécanique 2 à tester. Le bras robotisé 5 peut en outre être configuré pour déplacer le transducteur multiélément 3 autour de la pièce mécanique 2 afin de faciliter l'examen de l'ensemble du volume de la pièce ; et 5 - un dispositif de validation 6 de la pièce mécanique 2, conforme à l'invention, apte à utiliser les signaux mesurés par le transducteur multiéléments en réponse aux ondes ultrasonores propagées dans la pièce mécanique 2 pour déterminer si la pièce mécanique 2 contient ou non des indications de défauts. 10 On note qu'aucune limitation n'est attachée à la forme du transducteur multiélément 3. Dans le mode de réalisation décrit ici, celui-ci est linéaire, mais d'autres formes peuvent être envisagées (par exemple carrée, rectangulaire, etc.). Le contrôle non destructif de la pièce mécanique 2 peut en 15 outre être réalisé en immergeant la pièce mécanique 2 dans une cuve de contrôle (non représentée sur la figure 1) contenant un corps liquide, tel que par exemple de l'eau, pour favoriser le transfert de l'énergie ultrasonore du transducteur multiélément vers la pièce mécanique 2. Dans le mode de réalisation décrit ici, le dispositif de validation 20 6 est un ordinateur, dont l'architecture matérielle est illustrée sur la figure 2. Il comprend, notamment un processeur 7, une mémoire vive 8, une mémoire morte 9, une mémoire flash non volatile 10 ainsi que des moyens de communication 11 permettant notamment au dispositif de 25 validation 6 de communiquer avec l'amplificateur multivoie 4 pour récupérer notamment les signaux mesurés par le transducteur multiélément 3. Ces moyens de communication comprennent par exemple un bus de données numériques ou si le dispositif de validation 6 et l'amplificateur multivoie 4 sont reliés via un réseau de télécommunications 30 (local ou autre, filaire ou sans fil, etc.), une carte réseau ou une interface permettant de communiquer sur ce réseau. La mémoire morte 9 du dispositif de validation 6 constitue un support d'enregistrement, lisible par le processeur 7 et sur lequel est enregistré un programme d'ordinateur (ou un ensemble de programmes 35 d'ordinateurs) PROG. Ce programme d'ordinateur PROG comprend une pluralité d'instructions exécutables par le processeur 7, et qui définissent 13 différents modules fonctionnels et logiciels ici du dispositif de validation 6, ces modules fonctionnels s'appuyant sur et/ou commandant les éléments matériels 7-11 du dispositif de validation mentionnés précédemment. Ils comprennent notamment ici, comme illustré sur la figure 1 : 5 - un premier module de détermination 6A, configuré pour déterminer au moins une première matrice de composantes temporelles à partir des signaux kij(t) mesurés par le transducteur multiélément 3, i,j=1,...,N, pris sur une durée prédéterminée notée AT autour d'une référence temporelle notée T; 10 - un second module de détermination 6B, configuré pour déterminer au moins une deuxième matrice de composantes fréquentielles correspondant à une fréquence déterminée à partir de la fréquence des ondes ultrasonores émises par le transducteur multiélément 3, ce module 6B étant configuré pour appliquer à cet effet une transformée 15 de Fourier aux composantes temporelles de ladite au moins une première matrice ; - un module de filtrage 6C, configuré pour filtrer ladite au moins une deuxième matrice en la projetant sur un sous-espace de diffusion simple ; et 20 - un module d'analyse 6D, configuré pour vérifier l'intégrité de la pièce mécanique en utilisant ladite au moins une deuxième matrice filtrée obtenue par le module de filtrage 6C. Le module de filtrage 6C du dispositif de validation 6 comprend lui-même une pluralité de modules, définis par des instructions du 25 programme d'ordinateur PROG, comprenant notamment : - un troisième module de détermination 6C1, configuré pour déterminer le sous-espace de diffusion simple utilisé par le module de filtrage 6C. Ce sous-ensemble de diffusion simple est défini à partir d'une base de 2N-1 matrices notées Em, m=1,...,2N-1. Pour déterminer cette base, le 30 troisième module de détermination 6C1 est configuré pour : o pour chaque couple d'éléments distincts e(i) et e(j) du transducteur multiélément 3, i, j=1,...,N, calculer numériquement en utilisant un algorithme de tracé de rayons, pour une pluralité D de diffuseurs notés Pd, d=1...,D répartis en 35 divers points de la pièce mécanique 2 et situés à une profondeur déterminée en fonction de la référence temporelle T, 14 un temps de propagation tijd d'une onde ultrasonore émise par l'élément e(i), diffractée par le diffuseur Pd et rétrodiffusée par l'élément e(j), ledit temps de propagation étant calculé en tenant compte de la forme de la pièce mécanique, du placement 5 du transducteur et des milieux dans lesquels se propage l'onde ultrasonore entre les éléments e(i) et e(j) du transducteur et le diffuseur Pd; o estimer une contribution de la diffusion simple dans les composantes fréquentielles de la deuxième matrice à partir des 10 temps de propagation calculés pour chaque couple d'éléments distincts e(i) et e(j) du transducteur, et obtenir une troisième matrice à partir des contributions estimées ; o extraire des composantes fréquentielles d'une anti-diagonale centrale de la troisième matrice pour construire la base de 15 matrices Em, m=1,...,2N-1 définissant le sous-espace de diffusion simple ; - un module de projection 6C2, configuré pour projeter ladite au moins une deuxième matrice sur chacune des matrices Em, m=1,..,2N-1 de la base de 2N-1 matrices ; et 20 - un module d'évaluation 6C3 configuré pour évaluer ladite au moins une deuxième matrice filtrée en sommant les 2N-1 matrices Em pondérées respectivement par les projections obtenues par le module de projection 6C2. Les fonctions de ces modules 6A à 6D sont décrites plus en 25 détail maintenant en référence aux étapes du procédé de contrôle non destructif selon l'invention. La figure 3 représente les principales étapes du procédé de contrôle non destructif mis en oeuvre par le système de contrôle non destructif 1 pour contrôler l'intégrité de la pièce mécanique 2, dans un 30 mode particulier de réalisation de l'invention. On suppose que le transducteur multiélément 3 est positionné et orienté par le bras robotisé 5 selon une position POS dite courante dans laquelle il est placé en vis-à-vis de la pièce mécanique 2 (étape E10). Aucune limitation n'est attachée au choix de la position courante POS. 35 Toutefois, préférentiellement, le transducteur multiélément 3 est positionné à une distance inférieure à 1m de la pièce mécanique 2, qui 15 peut être ajustée de sorte à permettre la mesure d'un niveau de signal suffisant avec le transducteur multiélément 3 pour pouvoir réaliser le filtrage numérique conformément à l'invention. Un tel niveau de signal peut être déterminé expérimentalement. Il est intrinsèquement lié ici aux 5 caractéristiques et aux propriétés de l'amplificateur multivoie 4 qui est utilisé pour amplifier les signaux reçus par le transducteur élément 3 de sorte à les rendre exploitables. On évitera également de placer le transducteur multiélément 3 en face d'arêtes saillantes de la pièce mécanique 2 pour éviter que le niveau de signal renvoyé par la pièce 10 mécanique 2 ne soit trop bas. Une fois le transducteur multiélément 3 positionné, l'amplificateur multivoie 4 est activé pour piloter le transducteur multiélément 3 pour réaliser des mesures (étape E20). Plus spécifiquement, chaque élément piézoélectrique e(i), i=1,...,N du 15 transducteur 3 est activé pour émettre chacun à leur tour une onde ultrasonore à destination de la pièce mécanique 2 à une fréquence donnée, notée ici f0. Suite à l'émission d'une onde ultrasonore par un élément piézoélectrique e(i), les autres éléments piézoélectriques e(j) distinct de l'élément e(i) sont chacun activés pour mesurer et enregistrer 20 sur une certaine durée paramétrable un signal temporel, noté kij(t) représentatif de l'onde ultrasonore rétrodiffusée et reçue par chacun de ces éléments e(j). Chaque signal temporel kij(t) ainsi mesuré caractérise donc la réponse du matériau de la pièce mécanique 2 lorsque l'élément piézoélectrique e(i) émet une onde ultrasonore à la fréquence f0 et 25 l'élément e(j) enregistre l'onde ultrasonore rétrodiffusée par le matériau. Les signaux kij(t) ainsi mesurés sont transmis par le transducteur multiélément 3 à l'amplificateur multivoie 4, qui les transmet à son tour au dispositif de validation 6. Ils forment ici les composantes d'une matrice dite inter-élément notée k(t). 30 Le dispositif de validation 6 détermine, via son premier module de détermination 6A, une première matrice notée kn(t) dont les composantes sont des signaux temporels extraits des signaux kij(t) mesurés, i,j=1,...,N, et transmis au dispositif de validation 6 (étape E30). Plus spécifiquement, le premier module de détermination 6A extrait de 35 chaque signal kij(t), une tranche temporelle kijn(t) de durée AT prise 16 autour d'une référence temporelle notée TO, et forme avec les signaux ainsi extraits la première matrice kr°. Autrement dit : kTO(t)=[kiiT0(0]i,i_i,...m Dans le mode de réalisation décrit ici, par souci de 5 simplification, le premier module de détermination 6A détermine une unique matrice kiim(t) dont les signaux sont des signaux temporels de durée AT autour de TO. Toutefois, dans un autre mode de réalisation, le premier module de détermination 6A peut déterminer plusieurs « premières » matrices distinctes, correspondant à des tranches 10 temporelles prises autour de références temporelles distinctes de TO, les durées de ces tranches temporelles pouvant être prises égales à ou être différentes de AT. Puis le dispositif de validation 6, via son deuxième module de détermination 6B, applique une transformée de Fourier aux composantes 15 temporelles de la première matrice kr° de sorte à déterminer une deuxième matrice notée 1 N La base des matrices Em, m=1,. ..,2N-1 ainsi construite par le module de détermination 6C1 du module de filtrage 6C définit le sous-espace de diffusion simple utilisé par le module de filtrage 6C pour filtrer la deuxième matrice 1 Où = Tr(AB*) où Tr(M) désigne la trace de la matrice M et M* 10 désigne la transposée conjuguée de la matrice M. Puis le module de filtrage 6C, via son module d'évaluation 6C3, évalue la deuxième matrice filtrée (notée ci-dessous FK1-°(wc)) en sommant l'ensemble des matrices de la base pondérées respectivement par les projections obtenues sur ces bases, soit (étape E65) : 2N-1 Fe° (wc) = Z pmEm m.1 15 Ceci clôture l'étape de filtrage E60 mise en oeuvre par le module de filtrage 6C du dispositif de validation 6. La matrice filtrée FKT°(a)c) obtenue est ensuite fournie au module d'analyse 6D du dispositif de validation 6 pour analyse et détection le cas échéant, d'indications de défaut dans la pièce mécanique 2 (étape E70 de la figure 3). 20 Si plusieurs matrices ont été déterminées à l'issue de l'étape E50, l'étape E60 est mise en oeuvre sur chacune de ces matrices, et l'ensemble des matrices filtrées est fourni au module d'analyse 6D. Aucune limitation n'est attachée à l'algorithme utilisé par le module d'analyse 6D pour vérifier l'intégrité de la pièce mécanique 2, et 25 détecter le cas échéant la présence d'un défaut dans la pièce mécanique 2 à partir de la ou des matrices filtrées qui lui sont fournies. Ainsi, par exemple, le module d'analyse 6D peut utiliser à cet effet une technique DORT de diagonalisation de l'opérateur de retournement temporel telle que décrite notamment dans le document de 30 C. Prada et al. intitulé « Eigenmodes of the time reversai operator : a 24 solution to selective focusing in multiple-target medias », Wave motion, vol. 20 n°2 pages 151-163, 1994. Il peut procéder à cet effet comme décrit dans le document WO 2010/001027 en considérant un seuil sur la première valeur singulière de la matrice filtrée. 5 En variante, le module d'analyse 6D peut utiliser une méthode de focalisation totale telle que décrite dans le document de C. Holrnes et al. intutiié « post-processing of the full rnatrix of ultrasonic transmitreceive array data for non-destructive evaluation », NDF E. International, vol. 38, pages 701-711, 2005. 10 Ces exemples ne sont toutefois donnés qu'à titre illustratif et d'autres algorithmes peuvent être envisagés pour vérifier l'intégrité de la pièce mécanique 2 à partir de la matrice filtrée FO(coc). Si un défaut est détecté par le module d'analyse 6D lors de l'étape de vérification E70 (réponse oui à l'étape test E80), la pièce 15 mécanique 2 est déclarée non conforme (ou non valide) et rejetée (étape E90). Si aucun défaut n'est détecté (réponse non à l'étape test E80), la pièce mécanique 2 est validée (étape E100), et peut être utilisée conformément à sa destination finale (i.e. montée sur le moteur d'aéronef 20 à laquelle elle est destinée dans l'exemple envisagé ici). L'invention permet donc avantageusement un contrôle non destructif d'une pièce mécanique au moyen d'ondes ultrasonores émises par un transducteur multiéléments, quelle que soit la forme de cette pièce mécanique et quelle que soit la complexité de l'interface entre le 25 transducteur multiélément et la pièce mécanique. Le filtrage réalisé conformément à l'invention permet une réduction substantielle du bruit présent dans les mesures réalisées au moyen du transducteur multiélément ce qui permet d'améliorer l'efficacité du contrôle non destructif opéré ensuite. L'invention peut ainsi s'appliquer REVENDICATIONS1. Procédé de contrôle non destructif d'une pièce mécanique (2) par propagation d'ondes ultrasonores réalisée au moyen d'un transducteur multiélément (3) placé en vis-à-vis de la pièce mécanique et comprenant une pluralité d'éléments piézoélectriques e(1),e(2),...,e(N), N désignant un entier supérieur à 1, ledit procédé de contrôle comprenant : - une étape (E20) de mesure comprenant, pour chaque élément e(i) du transducteur, l'émission d'une onde ultrasonore à une fréquence donnée et la mesure, par chaque élément e(j) du transducteur distinct de l'élément e(I), d'un signal temporel noté kij(t) représentatif de l'onde ultrasonore rétrodiffusée reçue par l'élément e(j) ; - une étape (E30) de détermination d'au moins une première matrice de composantes temporelles à partir des signaux kij(t) mesurés, i,j=1,...,N, pris sur une durée prédéterminée autour d'une référence temporelle notée T; - une étape (E40,E50) de détermination d'au moins une deuxième matrice de composantes fréquentielles correspondant à une fréquence déterminée à partir de la fréquence de l'onde ultrasonore en appliquant une transformée de Fourier aux composantes temporelles de ladite au moins une première matrice ; - une étape (E60) de filtrage de ladite au moins une deuxième matrice comprenant une projection de ladite au moins une deuxième matrice de composantes fréquentielles sur un sous-espace de diffusion simple ; et - une étape de vérification (E70-E100) de l'intégrité de la pièce mécanique en utilisant ladite au moins une deuxième matrice filtrée ; le procédé de contrôle étant caractérisé en ce que le sous-espace de diffusion simple utilisé lors de l'étape de filtrage est déterminé au moyen d'un calcul numérique utilisant un algorithme de tracé de rayons. 2. Procédé de contrôle non destructif selon la revendication 1 dans lequel l'étape de filtrage (E60) de la deuxième matrice comprend : - une étape de détermination (E61-E63) du sous-espace de diffusion simple, ledit sous-espace de diffusion simple étant défini à partir d'une 26 base de 2N-1 matrices notées Em, m=1,...,2N-1, ladite étape de détermination comprenant : o pour chaque couple d'éléments distincts e(i) et e(j) du transducteur, i, j, =1,...,N, une étape de calcul numérique en 5 utilisant un algorithme de tracé de rayons, pour une pluralité D de diffuseurs notés Pd, d=1...,D répartis en divers points de la pièce mécanique et situés à une profondeur déterminée en fonction de la référence temporelle T, d'un temps de propagation tijd d'une onde ultrasonore émise par l'élément e(i), 10 diffractée par le diffuseur Pd et reçue par l'élément e(j), le temps de propagation tijd étant calculé en tenant compte de la forme de la pièce mécanique, du placement du transducteur et des milieux dans lesquels se propage l'onde ultrasonore entre les éléments e(i) et e(j) du transducteur et le diffuseur Pd ; 15 o une étape d'estimation d'une contribution de la diffusion simple dans les composantes fréquentielles de la deuxième matrice à partir des temps de propagation calculés pour chaque couple d'éléments distincts e(i) et e(j) du transducteur, ladite étape d'estimation fournissant une troisième matrice ; 20 o une étape d'extraction des composantes fréquentielles d'une anti-diagonale centrale de la troisième matrice pour construire la base de matrices Em, m=1,...,2N-1 définissant le sous-espace de diffusion simple ; une étape de projection (E64) de ladite au moins une deuxième 25 matrice sur chacune des matrices Em, m=1,..,2N-1 de la base de 2N-1 matrices ; et une étape d'évaluation (E65) de ladite au moins une deuxième matrice filtrée en sommant les matrices Em de la base pondérées par les projections obtenues lors de l'étape de projection. 30 3. Procédé de contrôle selon la revendication 2 dans lequel chaque temps de propagation tijd est calculé selon : tijd = minpe tid(Pe) + minps td es) où Pe désigne un point d'entrée de coordonnées (xe,ze) de l'onde ultrasonore dans la pièce mécanique et Ps un point de sortie de 27 coordonnées (xs,zs) de l'onde ultrasonore diffractée par le point Xd de coordonnées (xd,zd), avec : - \ tid( x) = I + (xe - xi)2 + (ze - zi)2 -,I (xd - xe)2 + (zd - ze)2 cl c2 et . tdj (x) = \ I (xd - xs)2 + (zd - zs)2 + - xj)2 + (zs - zj)2 c2 cl où cl désigne la vitesse de propagation de l'onde ultrasonore dans le 5 milieu séparant le transducteur de la pièce mécanique et c2 désigne la vitesse de propagation de l'onde ultrasonore dans la pièce mécanique. 4. Procédé de contrôle selon l'une quelconque des revendications 1 à 3 dans lequel la pièce mécanique (2) est une billette ou 10 un disque de soufflante ou un disque aubagé monobloc d'un moteur d'aéronef. 5. Procédé de contrôle selon l'une quelconque des revendications 1 à 4 dans lequel plusieurs premières matrices 15 correspondant à des références temporelles différentes et/ou à des durées prédéterminées différentes sont déterminées. 6. Procédé de contrôle selon l'une quelconque des revendications 1 à 5 dans lequel ladite au moins une deuxième matrice est 20 déterminée à la fréquence de l'onde ultrasonore. 7. Procédé de contrôle selon la revendication 6 dans lequel plusieurs deuxièmes matrices sont déterminées à une pluralité de fréquences contenues dans un intervalle prédéterminé défini autour de la 25 fréquence de l'onde ultrasonore. 8. Procédé de contrôle selon l'une quelconque des revendications 1 à 7 dans lequel l'étape de vérification (E70) de l'intégrité de la pièce mécanique utilise une technique DORT de diagonalisation de 30 l'opérateur de retournement temporel ou une méthode de focalisation totale. 28 9. Système de contrôle (1) non destructif d'une pièce mécanique (2) comprenant : - un transducteur rnultiélément (3) placé en vis-à-vis de la pièce 5 mécanique et apte à émettre des ondes ultrasonores, ledit transducteur comprenant une pluralité d'éléments piézoélectriques e(1),e(2),...,e(N), N désignant un entier supérieur à 1, ledit transducteur étant configuré pour émettre, via chaque élément e(i) du transducteur, une onde ultrasonore à une fréquence donnée, et à 10 mesurer, via chaque élément e(j) du transducteur distinct de l'élément e(i), un signal temporel noté kij(t) représentatif de l'onde ultrasonore rétrodiffusée reçue par l'élément e(j) ; - un dispositif de validation (6) comprenant : o un premier module de détermination (6A), configuré pour 15 déterminer au moins une première matrice de composantes temporelles à partir des signaux kij(t) mesurés, i,j=1,...,N, pris sur une durée prédéterminée autour d'une référence temporelle notée T; o un deuxième module de détermination (6B), configuré pour 20 déterminer au moins une deuxième matrice de composantes fréquentielles correspondant à une fréquence déterminée à partir de la fréquence de l'onde ultrasonore en appliquant une transformée de Fourier aux composantes temporelles de ladite au moins une première matrice ; 25 o un module de filtrage (6C), configuré pour filtrer ladite au moins une deuxième matrice en projetant ladite au moins une deuxième matrice de composantes fréquentielles sur un sous-espace de diffusion simple ; et o un module d'analyse (6D), configuré pour vérifier l'intégrité de 30 la pièce mécanique en utilisant ladite au moins une deuxième matrice filtrée ; le système de contrôle étant caractérisé en ce que ledit module de filtrage (6C) est configuré pour déterminer le sous-espace de diffusion simple au moyen d'un calcul numérique utilisant un algorithme de tracé de rayons. 35 29 10. Système de contrôle selon la revendication 9 dans lequel ledit module de filtrage (6C) comprend : un troisième module de détermination (6C1), configuré pour déterminer le sous-espace de diffusion simple à partir d'une base de 5 2N-1 matrices notées Em, m=1,...,2N-1, ledit troisième module de détermination étant configuré pour: o pour chaque couple d'éléments distincts e(i) et e(j) du transducteur, i, j, =1,...,N, calculer numériquement en utilisant un algorithme de tracé de rayons, pour une pluralité D de 10 diffuseurs notés Pd, d=1...,D répartis en divers points de la pièce mécanique et situés à une profondeur déterminée en fonction de la référence temporelle T, un temps de propagation tijd d'une onde ultrasonore émise par l'élément e(i), diffractée par le diffuseur Pd et rétrodiffusée par l'élément e(j), ledit 15 temps de propagation étant calculé en tenant compte de la forme de la pièce mécanique, du placement du transducteur et des milieux dans lesquels se propage l'onde ultrasonore entre les éléments e(i) et e(j) du transducteur et le diffuseur Pd ; o estimer une contribution de la diffusion simple dans les 20 composantes fréquentielles de la deuxième matrice à partir des temps de propagation calculés pour chaque couple d'éléments distincts e(i) et e(j) du transducteur, et obtenir une troisième matrice à partir des contributions estimées ; o extraire des composantes fréquentielles d'une anti-diagonale 25 centrale de la troisième matrice pour construire la base de matrices Em, m=1,...,2N-1 définissant le sous-espace de diffusion simple ; - un module de projection (6C2), configuré pour projeter ladite au moins une deuxième matrice sur chacune des matrices Em, m=1,..,2N-1 de la 30 base de 2N-1 matrices ; et - un module d'évaluation (6C3) configuré pour évaluer ladite au moins une deuxième matrice filtrée en sommant les matrices Em de la base pondérées par les projections obtenues lors de l'étape de projection. 35 11. Système de contrôle selon la revendication 9 ou 10 comprenant en outre un bras robotisé (5) sur lequel est placé le 30 transducteur multiélément, ledit bras robotisé étant configuré pour permettre un déplacement du transducteur multiélément autour de la pièce mécanique. 5 12. Système de contrôle selon l'une quelconque des revendications 9 à 11 dans lequel le transducteur (3) est linéaire.