Brevet FR3137453A1 — publié le 20240105

PROCÉDÉ D’IDENTIFICATION AUTOMATIQUE D’UNE SOURCE ACOUSTIQUE A PARTIR D’UN SIGNAL ACOUSTIQUE PRODUIT Ce procédé comporte : - une étape de reconstruction (58) du signal acoustique produit lors de l’apparition ou l’évolution du défaut, cette étape de reconstruction comportant la construction d’une estimation S e (w) du signal produit au niveau de la position du défaut à partir d’un signal ultrasonore F j (t) mesuré par un capteur et à l’aide d’une estimation ou d’une mesure d’un produit R j (w)G j (w) qui relie, dans le domaine fréquentiel, le signal ultrasonore F j (w) mesuré au signal acoustique S(w) produit, et - une étape de classification (76) automatique de la source acoustique réalisée à partir de l’estimation du signal acoustique reconstruit. Fig. 2 PROCÉDÉ D’IDENTIFICATION AUTOMATIQUE D’UNE SOURCE ACOUSTIQUE A PARTIR D’UN SIGNAL ACOUSTIQUE PRODUIT L’invention concerne un procédé et un appareil d’identification automatique d’une source acoustique à partir d’un signal acoustique produit par l’apparition ou l’évolution d’un défaut dans une structure. L’invention concerne également un support d’enregistrement d’informations pour la mise en œuvre de ce procédé. De tels procédés d’identification sont utilisés pour identifier le type de défaut qui est apparu dans une structure. La structure est, par exemple, une plaque en plastique renforcée par des fibres ou similaires. Dans ce cas, la source acoustique à identifier est le type de défaut qui est apparu ou qui évolue. Dans ce texte, l’expression « type de défaut » désigne la nature physique du défaut. Par exemple, dans le cas des structures constituées d’une matrice en plastique renforcée par des fibres, il existe principalement quatre types de défaut à savoir : - des casses de fibre, - du délaminage, - de la décohésion de matrice-fibre, et - de la fissuration de la matrice. Ces procédés sont donc utilisés dans le domaine du contrôle non destructif ou du SHM (« STRUCTURAL HEALTH MONITORING »). Des procédés connus d’identification d’une source acoustique à partir du signal acoustique produit par l’apparition ou l’évolution d’un défaut dans la structure comportent les étapes suivantes : - l’instrumentation de la structure avec un ou plusieurs capteurs capables de mesurer des signaux ultrasonores, - la mesure des signaux ultrasonores causés par l’apparition ou l’évolution du défaut, puis - l’extraction de certaines caractéristiques physiques des signaux ultrasonores mesurés puis, - la classification automatique de la source acoustique dans une classe parmi plusieurs classes possibles à partir des caractéristiques physiques extraites. Des exemples de tels procédés connus sont décrits dans les articles suivants : 1) Pashmforoush, F. et Al : « Damage characterization of glass/epoxy composite under three-point bending test using acoustic emission technique », Journal of materials engineering and performance, vol. 21(7), pp 1380-1390, 2012, 2) Morizet, N. et AL : « Classification of acoustic emission signals using wavelets and Random Forests : Application to localized corrosion », Mechanical Systems and Signal Processing, 70, 1026-1037, 2016 3) Hisham A Sawan et AL : « Unsupervised learning for classification of acoustic emission events from tensile and bending experiments with open-hole carbon fiber composite samples », Composites Science and Technology, Volume 107, 2015, Pages 89-97. Les articles numérotés de 1) à 3) ci-dessus sont, dans ce texte, désignés par, respectivement, les références suivantes : Pashmforoush2012, Morizet2016 et Sawan2015. Un tel procédé connu est également décrit dans la demande RU2737235C1. Il est souligné que le procédé décrit dans la demande RU2737235C1 comporte en plus une étape de localisation de la position de la source acoustique avant de réaliser la classification. Dans ce procédé, la position de la source acoustique est utilisée pour déterminer la distance qui sépare chaque capteur de la source acoustique. Ensuite, le type de défaut est identifié à partir d’un ratio entre l’amplitude de certaines composantes fréquentielles du signal ultrasonore mesuré par un capteur et la distance qui sépare ce capteur de la source acoustique. Ces procédés connus d’identification présentent notamment l’avantage qu’il n’est pas nécessaire d’équiper la structure avec un émetteur électronique d’ondes telles que des ondes ultrasonores. En effet, c’est l’onde ultrasonore générée par le défaut lui-même qui est utilisée pour réaliser cette identification. Les procédés connus d’identification fonctionnent correctement pour des structures de petites dimensions. Par contre, ces procédés fonctionnent moins bien avec des structures de grandes dimensions dans lesquelles la distance entre la source acoustique et les capteurs peut être beaucoup plus grande. L’invention vise à proposer un procédé d’identification d’une source acoustique qui présente les mêmes avantages que les procédés connus tout en fonctionnant aussi bien avec des structures de petites dimensions qu’avec des structures de grandes dimensions. Elle a donc pour objet un procédé d’identification automatique d’une source acoustique à partir d’un signal acoustique S(t) produit par l’apparition ou l’évolution d’un défaut dans une structure, ce procédé comportant les étapes suivantes : a) l’instrumentation de la structure en fixant sur cette structure au moins un capteur, puis b) la mesure, par le capteur, d'un signal ultrasonore F j (t) causé par le signal acoustique produit par l’apparition ou l’évolution du défaut dans la structure, ce signal ultrasonore mesuré étant relié, dans le domaine fréquentiel, au signal acoustique produit par la relation suivante : F j (w) = R j (w)G j (w)S(w), où : - j est un indice identifiant le capteur, - w est la pulsation en radian, - F j (w) est la transformée de Fourier du signal ultrasonore F j (t) mesuré dans le domaine temporel par le capteur j, - S(w) est la transformée de Fourier du signal acoustique S(t) produit par l’apparition ou l’évolution du défaut dans la structure, - R j (w) est la réponse, dans le domaine fréquentiel, du capteur j, - G j (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique dans la structure entre la position du défaut et l’emplacement du capteur j, puis c) la classification automatique de la source acoustique dans une classe de sources acoustiques choisie parmi plusieurs classes de sources acoustiques possibles, dans lequel : - le procédé comporte également entre les étapes b) et c), une étape d) de reconstruction du signal acoustique produit lors de l’apparition ou l’évolution du défaut, cette étape de reconstruction comportant la construction d’une estimation S e (w) du signal produit au niveau de la position du défaut à partir du signal ultrasonore F j (t) mesuré par le capteur et à l’aide d’une estimation ou d’une mesure du produit R j (w)G j (w) qui relie, dans le domaine fréquentiel, le signal ultrasonore F j (w) mesuré au signal acoustique S(w) produit, et - la classification automatique est réalisée à partir de l’estimation du signal acoustique obtenue à l’issue de l’étape d). Les modes de réalisation de ce procédé peuvent comporter une ou plusieurs des caractéristiques suivantes : 1) - l’étape a) d’instrumentation comporte la fixation sur la structure d’au moins trois capteurs à des emplacements dont les coordonnées sont connues dans un repère fixe par rapport à la structure, - l’étape b) comporte la mesure, par chacun des capteurs, d'un signal ultrasonore F j (t) respectif causé par le signal acoustique produit par l’apparition ou l’évolution du défaut dans la structure, chacun des signaux ultrasonores mesurés étant relié, dans le domaine fréquentiel, au signal acoustique produit par la relation suivante : F j (w) = R j (w)G j (w)S(w), - le procédé comporte entre les étapes b) et d), une étape de localisation, dans le repère fixe, de la position du défaut qui a produit le signal acoustique à partir des mesures réalisées par les capteurs lors de l’étape b) et des coordonnées connues de ces capteurs dans ce repère fixe, et - l’étape d) de reconstruction comporte l’obtention de l’estimation ou de la mesure du produit R j (w)G j (w) à l’aide de la position du défaut obtenue à l’issue de l’étape de localisation. 2) l’étape d) de reconstruction comporte : - pour chaque capteur, une opération de construction d’une estimation G e,j (w) de la fonction G j (w) de propagation du signal acoustique dans la structure entre la position où le défaut apparaît et l’emplacement de ce capteur, à l’aide de la relation suivante : où Γ est une phase définie par la relation suivante : où : - « arg » est la fonction qui retourne l’argument d’un nombre complexe, et - W j sont des poids prédéfinis, puis - la construction de l’estimation S e (w) du signal S(w) à partir des termes (R j (w)G e,j (w))*F j (w) ou des termes F j (w)/(R j (w)G e,j (w)), où le symbole «(...)* » désigne le conjugué de la fonction complexe entre parenthèses. 2) chaque poids W j est défini par la relation suivante : où : - P s sont les coordonnées dans le repère fixe du défaut obtenues à l’issue de l’étape de localisation, - P j sont les coordonnées dans le repère fixe du capteur j, - c est la célérité à laquelle se propage le signal acoustique à l’intérieur de la structure, et - w est la pulsation en radian. 4) chaque fonction R j (w) est égale à un quelle que soit la valeur de la pulsation w. 5) la fonction G e,j (w) de propagation est une fonction de Green. 6) - avant l’étape d) de reconstruction, le procédé comporte une étape d’apprentissage du produit R j (w)G j (w) des fonctions R j (w) et G j (w) pour différents emplacements possibles du défaut, cette étape d’apprentissage comportant les opérations suivantes : i) pour chaque capteur, la mesure ou l’estimation par simulation numérique d’un signal ultrasonore F j,k (t) mesuré par ce capteur lorsqu’un signal acoustique connu S c,k (t) est appliqué sur la structure en un emplacement P k dont les coordonnées sont connues dans le repère fixe, chacun des signaux ultrasonores F j,k (t) étant relié, dans le domaine fréquentiel, au signal acoustique connu S c,k (t) par la relation suivante : F j,k (w) = R j (w)G j,k (w)S c,k (w), où : - k est un indice identifiant l’emplacement P k dont les coordonnées sont connues, - F j,k (w) est la transformée de Fourier du signal ultrasonore F j,k (t) mesuré ou estimé dans le domaine temporelle par le capteur j lorsque le signal acoustique connu est appliqué à l’emplacement P k, - S c,k (w) est la transformée de Fourier du signal acoustique connu S c,k (t) appliqué à l’emplacement P k , - R j (w) est la réponse, dans le domaine fréquentiel, du capteur j, - G j,k (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique connu S c,k (t) dans la structure entre l’emplacement P k et l’emplacement du capteur j, puis ii) pour chaque capteur, la mémorisation, en association avec le capteur j et les coordonnées de l’emplacement P k , du produit R j (w)G j,k (w) des fonctions R j (w) et G j,k (w), le produit R j (w)G j,k (w) étant obtenu à l’aide de la relation suivante : R j (w)G j,k (w) = F j,k (w)/S c,k (w), iii) la répétition des opérations i) et ii) pour plusieurs emplacements P k possibles, - l’étape d) de reconstruction comporte les opérations suivantes : iv) la sélection des produits R j (w)G j,k (w) mémorisés en association avec l’emplacement P k le plus proche de la position du défaut obtenue à l’issue de l’étape de localisation, puis v) la détermination de l’estimation S e (w) du signal S(w) à partir des termes (R j (w)G j,k (w))*F j (w) ou des termes F j (w)/(R j (w)G j,k (w)), où : - les produits R j (w)G j,k (w) sont ceux sélectionnés lors de l’opération iv), et - le symbole «(...)* » désigne le conjugué de la fonction complexe entre parenthèses. Procédé selon l’une quelconque des revendications précédentes, dans lequel : - le procédé comporte la construction de l’estimation S e (t), dans le domaine temporel, du signal acoustique S(t) en appliquant une transformation de Fourier inverse à l’estimation S e (w), puis - la classification automatique est réalisée à partir de l’estimation S e (t) dans le domaine temporel. L’invention a également pour objet un support d’enregistrement d’informations lisible par un microprocesseur, caractérisé en ce qu’il comporte des instructions non-transitoires pour l’exécution des étapes b) à d) du procédé d’identification ci-dessus, lorsque ces instructions sont exécutées par le microprocesseur. L’invention a également pour objet un appareil d’identification automatique d’une source acoustique à partir d’un signal acoustique S(t) produit par l’apparition ou l’évolution d’un défaut dans une structure, cet appareil comportant : - au moins un capteur fixé sur cette structure, ce capteur étant apte à mesurer un signal ultrasonore F j (t) causé par le signal acoustique produit par l’apparition ou l’évolution du défaut dans la structure, ce signal ultrasonore mesuré étant relié, dans le domaine fréquentiel, au signal acoustique produit par la relation suivante : F j (w) = R j (w)G j (w)S(w), où : - j est un indice identifiant le capteur, - w est la pulsation en radian, - F j (w) est la transformée de Fourier du signal ultrasonore F j (t) mesuré dans le domaine temporel par le capteur j, - S(w) est la transformée de Fourier du signal acoustique S(t) produit par l’apparition ou l’évolution du défaut dans la structure, - R j (w) est la réponse, dans le domaine fréquentiel, du capteur j, - G j (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique dans la structure entre la position où le défaut apparaît et l’emplacement du capteur j, - un calculateur électronique configuré pour exécuter les étapes suivantes : - l’acquisition des mesures du capteur, puis - la classification automatique de la source acoustique dans une classe de sources acoustiques choisie parmi plusieurs classes de sources acoustiques possibles, dans lequel le calculateur électronique est également configuré pour : - avant la classification automatique, exécuter une étape de reconstruction du signal acoustique produit lors de l’apparition ou l’évolution du défaut, cette étape de reconstruction comportant la construction d’une estimation du signal produit au niveau de la position du défaut à partir du signal ultrasonore F j (t) mesuré par le capteur et à l’aide d’une estimation ou d’une mesure du produit R j (w)G j (w) qui relie, dans le domaine fréquentiel, le signal ultrasonore F j (w) mesuré au signal acoustique S(w) produit, et - la classification automatique est réalisée à partir de l’estimation du signal acoustique obtenue à l’issue de l’étape de reconstruction. L’invention sera mieux comprise à la lecture de la description qui va suivre donnée uniquement à titre d’exemple non limitatif et faite en se référant aux dessins sur lesquels : - la est une illustration schématique de l’architecture d’un appareil d’identification automatique d’une source acoustique ; - la est un organigramme d’un procédé d’identification automatique d’une source acoustique à l’aide de l’appareil de la ; - la est un graphe illustrant les résultats expérimentaux obtenus en mettant en œuvre le procédé de la , et - la est un organigramme d’un autre procédé d’identification susceptible d’être mis en œuvre par l’appareil de la . Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. Dans la suite de cette description, les caractéristiques et fonctions bien connues de l’homme du métier ne sont pas décrites en détail. CHAPITRE I : Exemples de modes de réalisation La représente un appareil 2 d’identification automatique d’une source acoustique à partir du signal acoustique produit par l’apparition ou l’évolution d’un défaut dans une structure 4. Dans ce cas, la source acoustique est le défaut qui apparaît ou évolue. Dans ce premier mode de réalisation, la structure 4 peut être toute structure dans laquelle, lorsqu’un défaut apparaît ou évolue, un signal acoustique S(t) est généré par ce défaut. L’essentiel de la puissance du signal S(t) est généralement situé dans la bande des fréquences ultrasonores, c’est-à-dire dans une bande de fréquences allant de 16 kHz à 10 MHz. En général, le spectre de puissance du signal S(t) ne recouvre qu’une partie de la bande des fréquences ultrasonores. La structure 4 est réalisée dans un matériau qui permet à ce signal S(t) de se propager à l’intérieur de cette structure sur une distance suffisante pour pouvoir être mesurée par un capteur et distinguée d’un bruit de fond. A titre d’illustration, ici, la structure 4 est un tube en matériau composite formé d’une matrice en plastique renforcée par des fibres. Ce tube s’étend le long d’un axe 6 de révolution sur une distance de 2 mètres. Sa section transversale est constante et circulaire sur toute sa longueur. Un repère R est attaché sans aucun degré de liberté à la structure 4. La position de chaque point de la structure 4 est repérée par des coordonnées dans le repère R. Ici, le repère R comporte un axe X confondu avec l’axe 6 de révolution et un axe Y perpendiculaire à l’axe X. Par exemple, l’axe Y est vertical et l’axe X est horizontal. Par la suite, étant donné que dans le cas particulier décrit ici, la structure 4 est invariante par toute rotation autour de l’axe 6, la position de chaque point P de la structure 4 est définie par des coordonnées cylindriques (x p , θ p ), où : - x p représente l’abscisse de ce point le long de l’axe X du repère R, et - θ p est l’angle entre un vecteur O k P et une direction parallèle à l’axe Y, où le vecteur O k P est le vecteur, perpendiculaire à l’axe 6 qui s’étend du point O k situé sur l’axe 6 jusqu’au point P de la structure 4. L’appareil 2 comporte : - un calculateur électronique 10, - une interface homme/machine 12 raccordée au calculateur 10, et - N capteurs C 1 à C N raccordés au calculateur 10. Le calculateur 10 comporte un microprocesseur programmable 14 et une mémoire 16. La mémoire 16 comporte les instructions et les données nécessaires à l’exécution du procédé de la ou de la , lorsque ces instructions sont exécutées par le microprocesseur 14. L’interface homme/machine 12 est capable de communiquer, de façon directement intelligible par un être humain, les résultats de la mise en œuvre du procédé d’identification de la ou 4. Par exemple, l’interface 12 comporte un écran. Chaque capteur C j mesure le signal ultrasonore produit par l’apparition ou l’évolution d’un défaut dans la structure 4. Dans ce texte, l’indice j est un nombre entier compris entre 1 et N qui identifie le capteur C j . Le calculateur 10 acquiert les mesures de chaque capteur C j pour les traiter. L’appareil 2 comporte au moins trois capteurs C j et, le plus souvent, au moins quatre ou six ou huit capteurs C j . Ici, à titre d’illustration, le nombre total N de capteur C j est égal à huit. Sur la , seuls les capteurs C 1 , C 2 , C j et C N sont représentés. La représentation des autres capteurs a été remplacée par des pointillés pour simplifier la . Ici, les capteurs C j sont tous structuralement identiques les uns aux autres et se distinguent les uns des autres seulement par leurs positions dans le repère R. Par exemple, les capteurs C j sont des capteurs piézoélectriques dont la bande passante recouvre au moins la bande [100 kHz ; 500 kHz] et, de préférence, la bande [50 kHz ; 0,5 MHz] ou [20 kHz ; 0,5 MHz]. Chaque capteur C j est fixé sur la structure 4 à un emplacement respectif P j où il est capable de mesurer un signal ultrasonore qui se propage dans la structure 4. Les coordonnées (x j ; θ j ) de chaque emplacement P j sont connues et mémorisées dans la mémoire 16. La fréquence d’échantillonnage f e des mesures des capteurs C j est élevée c’est-à-dire typiquement supérieure à 1 MHz. Ici, à titre d’illustration, les emplacements P j sont uniformément répartis le long d’un axe parallèle à l’axe 6. Toutefois, d’autres répartitions des emplacements P j sur la surface de la structure 4 sont possibles. Par exemple, plusieurs des emplacements P j peuvent aussi être répartis le long de la circonférence circulaire de la structure 4. De préférence, la distance entre l’un quelconque des emplacements P j et l’emplacement P j+1 le plus proche, est calculée de manière à assurer une bonne couverture de la structure, à savoir qu’au moins deux capteurs peuvent mesurer un signal significatif pour une source en tout point de la structure. Le fonctionnement de l’appareil 2 va maintenant être décrit en référence au procédé de la . Lors d’une étape 50 d’instrumentation, les capteurs C j sont fixés chacun à un emplacement P j respectif sur la structure 4. Les coordonnées (x j ; θ j ) de chaque position P j sont enregistrées dans la mémoire 16 en association avec l’identifiant j du capteur C j . Débute alors une phase 52 d’utilisation de l’appareil 2 pour détecter et identifier une source acoustique. Lors d’une étape 54, chaque capteur C j mesure le signal ultrasonore F j (t) qui se propage dans la structure 4. Typiquement, l’enregistrement et le traitement du signal par chaque capteur est déclenché par la détection d’une onde par ce capteur, par exemple, par le passage du signal ultrasonore F j (t) au-dessus d’un seuil. L’enregistrement s’arrête lorsque l’onde n’est plus détectée, par exemple lorsque le signal ultrasonore est en-dessous du seuil pendant une durée significative déterminée au préalable. Le seuil et la durée significative enregistrée après la dernière détection de l’onde sont donc choisis de manière à ne pas perdre d’information. La durée enregistrée est typiquement de l’ordre de la milliseconde. A l’inverse, la durée significative après la dernière détection est suffisamment courte et le seuil suffisamment élevé pour que, dans la très grande majorité des cas, un seul défaut apparaît ou évolue pendant l’enregistrement et le bruit ne déclenche pas l’enregistrement. Par exemple, pour cela, généralement, la durée significative après dernière détection est de l’ordre de quelques centaines de µs. En même temps, lors de l’étape 54, les mesures des capteurs C j sont acquises par le calculateur 10 à la fréquence d’échantillonnage f e . La fréquence f e est suffisamment grande pour permettre d’acquérir les signaux F j (t) mesurés en évitant ou en éliminant les phénomènes de repliement de spectres. Par exemple, ici, la fréquence f e est supérieure à 1 MHz. Avantageusement, la fréquence f e est deux fois supérieure à la plus grande borne de la bande passante des capteurs C j . Le signal F j (t) mesuré par chaque capteur C j est relié, dans le domaine fréquentiel, au signal S(t) généré par le défaut par la relation suivante : F j (w) = R j (w)G j (w)S(w), où : - j est un indice identifiant le capteur C j , - w est la pulsation en radian, - F j (w) est la transformée de Fourier du signal ultrasonore F j (t) mesuré dans le domaine temporel par le capteur C j , - S(w) est la transformée de Fourier du signal acoustique S(t) produit par l’apparition ou l’évolution du défaut dans la structure, - R j (w) est la réponse, dans le domaine fréquentiel, du capteur C j , - G j (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique dans la structure entre la position P s où le défaut apparaît et l’emplacement P j du capteur C j . Une fois que les signaux F j (t) ont été acquis, lors d’une étape 56, le calculateur 10 localise la position P s du défaut qui a généré les signaux ultrasonores acquis par les capteurs C j. Cette localisation consiste à déterminer les coordonnées (x s ; θ s ) de la position P s dans le repère R à partir des signaux F j (t) mesurés et acquis et des coordonnées connues des emplacements P j . De plus ici, la célérité c à laquelle se propage le signal S(t) dans la structure 4 est également déterminée lors de l’étape 56. Par exemple ici, pour déterminer les coordonnées (x s ; θ s ) et la célérité c, la méthode de triangulation décrite dans l’article suivant est implémentée : Zhang, F. et AL : « Evaluation of acoustic emission source localization accuracy in concrete structures », Structural Health Monitoring, vol. 19(6), pages 2063-2074, 2020. Plus précisément, les instants t j d’arrivée du signal S(t) au niveau de chaque emplacement P j sont relevés à partir des signaux ultrasonores F j (t) acquis. Par exemple, l’instant t j correspond au premier instant où l’amplitude du signal F j (t) dépasse un seuil prédéterminé. Ensuite, la position P s et la célérité c sont pris égales à la position et à la célérité qui minimise la fonction de coût J(P s ; c) suivante : Pour cela, un algorithme de minimisation de la fonction de coût J(P s ; c) est implémenté. Par exemple, il peut s’agir d’un algorithme tel que l’algorithme de NEWTON ou l’algorithme du SIMPLEXE. Des algorithmes génétiques peuvent aussi être mis en œuvre. Une fois la localisation de la source acoustique terminée, lors d’une étape 58 de reconstruction, le calculateur 10 construit une estimation S e (w) du signal S(w) produit à la position P s par le défaut. Cette étape vise à compenser la propagation du signal S(t) entre la position P s et les emplacements P j . Pour cela, dans ce premier mode de réalisation, le calculateur 10 utilise les signaux F j (t), la position P s et la célérité c déterminées lors de l’étape 56. Plus précisément, dans ce premier mode de réalisation, une méthode de reconstruction dite « aveugle » est mise en œuvre. Cette méthode est qualifiée de « aveugle » car elle ne nécessite pas que chaque fonction G j (w) de propagation du signal S(t) de la position P s jusqu’à la position P j soit apprise au préalable lors d’une étape d’apprentissage. Ces méthodes sont plus connues sous le terme anglais de « Blind source deconvolution »). Ici, la méthode mise en œuvre par le calculateur 10 est décrite dans l’article suivant : Sabra, K. G. et AL : « Blind deconvolution in ocean waveguides using artificial time reversal », The Journal of the Acoustical Society of America, vol. 116(1), pages 262-271, 2004. Par la suite, cet article est désigné par la référence Sabra2004. Selon cette méthode, lors d’une opération 60, le calculateur 10 construit une estimation G e,j (w) de chaque fonction G j (w) de propagation du signal S(t) de la position P s jusqu’à l’emplacement P j . Ici, chaque fonction G e,j (w) de propagation est une fonction de Green. Pour cela, l’estimation G e,j (w) de la fonction G j (w) est prise égale au signal F j (t) normalisé par la norme L 2 de tous les signaux F j (t) et multiplié par une phase estimée Γ. Plus précisément, l’estimation G e,j (w) est construite à l’aide des relations suivantes : où : - i est le nombre imaginaire tel que i 2 =-1, - Γ est une phase définie par la relation suivante : où : - « arg » est la fonction qui retourne l’argument d’un nombre complexe, et - W j sont des poids prédéfinis. De plus, ici, pour rendre le procédé de la indépendant des caractéristiques et des propriétés physiques de la structure 4, les poids W j sont calculés à l’aide de la relation suivante : où : - i est le nombre imaginaire tel que i 2 =-1, - ||P s – P j || est la distance qui sépare les positions P s et P j , - c est la célérité déterminée lors de l’étape 56, et - w est la pulsation en radian. Ensuite, lors d’une opération 62, l’estimation S e (w) est construite par retro-propagation des signaux F j (t) au niveau de la position P s . Pour cela, ici, chaque signal F j (t) est rétro-propagé au niveau de la position P s . La rétro-propagation d’un signal F j (t) consiste à multiplier le signal F j (w) par le conjugué de la fonction complexe de propagation. Autrement dit, le signal F j (w) rétro-propagé au niveau de la position P s est donné par la relation suivante : (R j (w)G e,j (w))*F j (w), où le symbole «(...)* » désigne le conjugué de la fonction complexe R j (w)G e,j (w) entre parenthèses. L’estimation S e (w) est donc construite à partir des termes (R j (w)G e,j (w))*F j (w) et chaque terme (R j (w)G e,j (w))*F j (w) est le produit du conjugué de la fonction R j (w)G e,j (w) par la fonction F j (w). Ici, l’estimation S e (w) est prise égale à la moyenne des N termes (R j (w)G e,j (w))*F j (w). Ainsi, l’estimation S e (w) est construite à l’aide de la relation suivante : De plus, dans ce mode de réalisation, pour simplifier, chaque fonction R j (w) est égale à un quelle que soit la valeur de la pulsation w. Ainsi, l’estimation S e (w) est simplement construite à l’aide de la relation suivante : Une fois l’estimation S e (w) construite, dans ce mode de réalisation, lors d’une étape 70, l’estimation S e (t) du signal S(t) dans le domaine temporel est obtenue, par exemple, en appliquant une transformation de Fourier inverse à l’estimation S e (w). Ensuite, dans ce mode de réalisation, le calculateur 10 exécute une étape 72 de réduction de la dimensionnalité. Une telle étape est également connue sous le terme de « réduction de dimensions ». En effet, les estimations S e (w) et S e (t) sont généralement volumineuses, c’est-à-dire sont souvent composées de plusieurs dizaines de milliers ou centaines de milliers d’échantillons. L’étape 72 vise à réduire la quantité d’informations à traiter tout en conservant les informations essentielles qui vont permettre une identification fiable et robuste du défaut. Dans ce mode de réalisation, l’étape 72 consiste à extraire des caractéristiques physiques propres au signal S(t) à partir des estimations S e (w) et S e (t) précédemment construites. Par exemple, ici, il est procédé comme décrit dans la section 2.2 de l’article Morizet2016 sauf qu’aucune énergie dans une bande supérieure à 200 kHz n’est extraite. Enfin, une étape 76 de classification est exécutée par le calculateur 10. Lors de cette étape 76, le calculateur 10 classe le signal S(t) dans une classe choisie parmi plusieurs classes possibles à partir des caractéristiques extraites lors de l’étape 72. Pour cela, le calculateur 10 exécute une méthode de classification automatique. Par exemple, ici, le calculateur 10 exécute une méthode de classification non-supervisée connue. Cette méthode de classification non-supervisée connue est, par exemple, une méthode de classification automatique utilisant un modèle de mélange gaussien (« Gaussian Mixture Model » en anglais) ou la méthode connue sous le nom anglais de « k-means » et décrite, par exemple, dans l’article Pashmforoush2012. Ici, pour la mise en œuvre de cette méthode de classification, le nombre N k de classes a été déterminé au préalable. Par exemple, pour la structure 4, le nombre N k de classe a été choisi égal à trois. Le nombre N k de classes le plus approprié a ici été déterminé en essayant plusieurs valeurs possibles pour le nombre N k puis en calculant pour chacune de ces valeurs du nombre N k , la valeur de l’index DB (« DAVIES-BOULDIN »). La valeur du nombre N k pour laquelle la valeur de l’index DB est minimale a été retenue. Ce calcul de la valeur de l’index DB est par exemple décrit dans la section 3.2 de l’article Pashmforoush2012. Chaque classe est reliée à une nature physique du défaut qui a produit le signal acoustique classé. Ainsi, le procédé d’identification permet non seulement de détecter l’apparition ou l’évolution d’un défaut mais aussi d’identifier la nature physique de ce défaut. Ce qui a été décrit ci-dessus est répété pour chaque nouvel événement mesuré par les capteurs C j . La illustre sous la forme d’un graphe les résultats obtenus en déformant la structure 4 de façon répétée de manière à faire apparaître des défauts. Sur cette figure, chaque point correspond à une source acoustique et donc à l’apparition ou à l’évolution d’un défaut dans la structure 4. L’axe des ordonnées est gradué dans une unité arbitraire proportionnelle à l’amplitude du signal S e (t) et l’axe des abscisses en Hertz. Pour chaque source acoustique identifiée, l’amplitude de l’estimation S e (t) et le centroïde fréquentiel de l’estimation S e (w) ont été calculés puis cette source acoustique a été placée sur le graphe de la en utilisant la valeur de son centroïde fréquentiel comme abscisse et l’amplitude calculée comme ordonnée. Le centroïde fréquentiel est calculé comme décrit dans la section 6.1.1 de l’article suivant : Geoffroy Peeters : « A large set of audio features for sound description (similarity and classification) in the CUIDADO project », 2004. Le centroïde fréquentiel est exprimé en Hertz. Chaque source acoustique a été classée dans l’une des trois classes possibles en mettant en œuvre le procédé de la . Le graphe de la comporte trois zones 31, 32 et 33 qui entourent chacune les points classés dans une classe respective. Les chevauchements entre ces trois zones 31, 32 et 33 sont très faibles, ce qui illustre le fait que la classification réalisée est robuste. La représente un autre procédé d’identification susceptible d’être mis en œuvre à la place du procédé de la par le calculateur 10. Le procédé de la est identique au procédé de la sauf que : - une étape 80 d’apprentissage de chaque fonction de propagation G j (w) est introduite entre l’étape 50 et la phase d’utilisation, et - la phase 52 d’utilisation est remplacée par une phase 82 d’utilisation. Lors de l’étape 80, un signal acoustique connu S c (t) à large spectre est appliqué à différents emplacements P k de la structure 4. Le signal S c (t) présente un spectre qui recouvre la bande passante de la structure 4, c’est-à-dire qui inclut toutes les fréquences susceptibles de se propager et d’être mesurées dans la structure 4 sans être trop atténuées. Par exemple, ici, le spectre du signal S c (t) recouvre la bande des fréquences ultrasonores. Par exemple, dans ce mode de réalisation, le signal S c (t) appliqué à chaque emplacement P k est une source de Hsu-Nielsen c’est-à-dire une source acoustique produite par la rupture d’une mine de crayon. Cette source étant connue, elle n’est pas décrite plus en détail. Les emplacements P k sont ici uniformément répartis sur la surface de la structure 4. Par exemple, les emplacements P k correspondent aux sommets d’un maillage qui recouvre la surface 4. Typiquement, les tuiles de ce maillage sont identiques ou similaires. La distance la plus courte entre deux emplacements P k contigus est par exemple comprise entre 1 mm et 100 mm ou entre 3 mm et 20 mm. Ici, la distance la plus courte entre deux emplacements P k contigus est comprise entre 5 mm et 10 mm. Le nombre N pk d’emplacements P k est typiquement proportionnel à la surface de la structure 4. Par exemple, le nombre N pk est pris égal au ratio S 4 /S ref , où : - S 4 est la surface de la structure 4, et - S ref est une surface de référence, par exemple, comprise entre 10 -2 m² et 1 m². Les coordonnées dans le repère R de chaque emplacement P k sont connues et mémorisées dans la mémoire 16. A chaque fois que le signal S c (t) est appliqué à un emplacement P k , chaque capteur C j mesure le signal ultrasonore F j,k (w) généré par le signal S c (t) qui se propage dans la structure 4. Le signal ultrasonore F j,k (w) est relié au signal acoustique S c (t) par la relation suivante : F j,k (w) = R j (w)G j,k (w)S c (w), où : - k est un indice identifiant l’emplacement P k dont les coordonnées sont connues, - F j,k (w) est la transformée de Fourier du signal ultrasonore F j,k (t) mesuré dans le domaine temporel par le capteur C j lorsque le signal acoustique connu S c (t) est appliqué à l’emplacement P k, - S c (w) est la transformée de Fourier du signal acoustique connu S c (t) appliqué à l’emplacement P k , - R j (w) est la réponse, dans le domaine fréquentiel, du capteur C j , - G j,k (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique connu S c (t) dans la structure 4 entre l’emplacement P k et l’emplacement du capteur C j . Ensuite, pour chaque emplacement P k où le signal S c (t) a été appliqué, le produit R j (w)G j,k (w) des fonctions R j (w) et G j,k (w), est enregistré dans la mémoire 16 en association avec l’identifiant j du capteur C j et l’emplacement P k où a été appliqué le signal S c (t). Le produit R j (w)G j,k (w) est égal au ratio F j,k (w)/S c (w) du signal ultrasonore F j,k (w) mesuré par le capteur C j et du signal acoustique connu S c (w). La phase 82 est identique à la phase 52 d’utilisation sauf que l’étape 58 de reconstruction est remplacée par une étape 88 de reconstruction. L’étape 88 de reconstruction comporte : - une opération 90 de sélection, dans la mémoire 16, des produits R j (w)G j,k (w) mémorisés, et - une opération 92 de détermination de l’estimation S e (w) du signal S(w). Lors de l’opération 90, le calculateur 10 sélectionne les produits R j (w)G j,k (w) associés à l’emplacement P k le plus proche de la position P s obtenue à l’issue de l’étape 56 de localisation. Lors de l’opération 92, l’estimation S e (w) est construite seulement à partir des produits R j (w)G j,k (w) sélectionnés lors de l’opération 90. Plus précisément, ici, l’estimation S e (w) est construite à partir des termes (R j (w)G j,k (w))*F j (w) de façon similaire à ce qui a été décrit dans le cas de l’opération 62. Par exemple, l’estimation S e (w) est construite à l’aide de la relation suivante : où les produits R j (w)G j,k (w) sont ceux sélectionnés lors de l’opération 90. CHAPITRE II : Variantes Variantes de la reconstruction « aveugle » : Les valeurs des poids W j peuvent être différentes. En particulier, dans un mode de réalisation simplifié, les valeurs des poids W j ne dépendent pas de la position P s du défaut. Par exemple, les poids W j sont des constantes toutes égales à +1. Dans ce cas, l’étape 56 de localisation peut être omise. De plus, lorsque la localisation du défaut est omise, le nombre N de capteurs C j peut être inférieur à trois. Par exemple, le nombre N de capteurs C j est égal à deux. Les valeurs des poids W j peuvent aussi être déterminées en fonction des caractéristiques physiques de la structure 4 comme, par exemple, la dispersion théorique des ondes guidées dans la structure. Dans ce dernier cas, les poids W j dépendent des caractéristiques physiques de la structure de sorte que le procédé d’identification est spécifiquement configuré pour identifier un signal acoustique dans cette structure. En variante, au lieu d’utiliser une rétro-propagation des signaux mesurés par les capteurs pour obtenir l’estimation du signal S e (w), il est possible d’utiliser un filtrage inverse (« inverse filtering » en anglais) comme décrit dans l’article Sabra2004. Dans ce dernier cas, ce sont les termes F j (w)/(R j (w)G e,j (w)) qui sont utilisés pour construire l’estimation S e (w) au lieu d’utiliser les termes (R j (w)G e,j (w))*F j (w). Dans ce cas, l’estimation S e (w) est obtenue à l’aide de la relation suivante : Dans ce dernier cas, de préférence, la bande passante de l’estimation S e (w) est limitée pour exclure les fréquences fortement atténuées, par exemple les fréquences élevées atténuées par la structure 4. Autrement dit, au-delà d’une fréquence de coupure f c prédéterminée, l’estimation S e (w) est nulle. Cela permet de limiter la sensibilité au bruit de l’estimation S e (w). En effet, pour les fréquences supérieures à f c , le produit R j (w)G e,j (w) est proche de zéro, de sorte que le bruit est fortement amplifié au-delà de la fréquence f c si la bande passante de l’estimation S e (w) n’est pas limitée. L’étape de reconstruction peut être réalisée en utilisant d’autres méthodes de reconstruction « aveugle » que celle décrite. De nombreux exemples d’autres méthodes de reconstruction « aveugle » peuvent être trouvés dans le domaine de la déconvolution de source aveugle plus connu sous le terme anglais de « blind source deconvolution ». L’estimation S e (w) peut aussi être construite à partir d’une somme pondérée des termes (R j (w)G e,j (w))*F j (w) ou F j (w)/(R j (w)G e,j (w)). Par exemple, l’estimation S e (w) est construite à l’aide de la relation suivante : où q j est un coefficient de pondération. La somme des coefficients q j est égale à un. Par exemple, plus le capteur C j est proche de la position P s , plus le coefficient q j est grand de manière à donner plus de poids aux mesures des capteurs qui sont les plus proches du défaut. Dans un cas simplifié, seuls les coefficients q j des M capteurs C j les plus proches sont non-nuls et les coefficients q j est (N-M) capteurs C j restants sont nuls, où M est un nombre entier inférieur à N. Par exemple, M est égal à un ou deux ou trois. En variante, la réponse R j (w) n’est pas une constante égale à un et l’estimation S e (w) est construite à l’aide de la relation suivante : Pour cela, si la réponse impulsionnelle de chaque capteur C j est connue, alors R j (w) peut être construit à partir de cette réponse impulsionnelle connue. La réponse R j (w) de chaque capteur C j peut aussi être mesurée. Dans ces derniers cas, la réponse R j (w) n’est généralement pas une constante. Variantes de la reconstruction par apprentissage : Le signal acoustique connu appliqué peut être généré par une autre source acoustique qu’une source de Hsu-Nielsen. Par exemple, le signal acoustique connu peut aussi être généré par un émetteur acoustique électronique. Lors de l’étape d’apprentissage, en variante, ce n’est pas le même signal acoustique connu qui est appliqué à chaque emplacement P k . Dans ce cas, le produit mémorisé en association avec chaque emplacement P k est obtenu à l’aide de la relation suivante : R j (w)G j,k (w) = F j,k (w)/S c,k (w), où S c,k (w) est le signal acoustique appliqué à l’emplacement P k . Comme dans le cas de la reconstruction «aveugle », l’estimation S e (w) peut aussi être construite par filtrage inverse. Elle peut aussi être construite à partir d’une somme pondérée ou non des termes (R j (w)G j ,k (w))*F j (w) ou F j (w)/(R j (w)G j,k (w)). Les données utilisées lors de l’étape 80 d’apprentissage peuvent être issues de simulations numériques. Variantes communes à tous les modes de réalisation : Une étape de pré-traitement des signaux ultrasonores F j (t) mesurés par les capteurs C j peut être réalisée avant l’étape de localisation et/ou avant l’étape de reconstruction afin de réaliser ces étapes sur des signaux ultrasonores pré-traités et non pas sur les signaux ultrasonores brutes directement délivrés par les capteurs. Des exemples de pré-traitements qui peuvent être appliqués aux signaux ultrasonores brutes F j (t) sont par exemple décrits dans la section 2.1 de l’article Morizet2016. D’autres modes de réalisation de l’étape 56 de localisation sont possibles. Par exemple, en variante, la localisation de la position P s du défaut est obtenue en implémentant la méthode décrite dans l’article suivant : Pearson M. R. et AL : « Improved acoustic emission source location during fatigue and impact events in metallic and composite structures », Structural Health Monitoring, vol. 16(4), pages 382-399, 2017. De plus, cet article décrit aussi une autre méthode pour relever les instants t j d’arrivée. Si la célérité c à laquelle se propage le signal S(t) à l’intérieur de la structure 4 est connue, il n’est alors pas nécessaire de la déterminer lors de l’étape 56 de localisation. Par exemple, la célérité c connue est pré-enregistrée dans la mémoire 16. Lorsque l’estimation S e (w) est construite à partir des termes F j (w)/(R j (w)G e,j (w)) ou F j (w)/(R j (w)G j,k (w)), si le dénominateur de ces fractions est proche de zéro dans une plage de fréquences, alors cette plage de fréquences est exclue et n’est pas prise en compte lors de l’étape de classification. Une telle situation peut se rencontrer, par exemple, lorsque la structure atténue fortement la propagation des ondes à des fréquences supérieures à f max . Par exemple, si le produit R j (w)G e,j (w) ou R j (w)G j,k (w) est inférieure à 0,1 au-delà de la fréquence f max , alors l’estimation S e (w) est limitée à la plage de fréquences [0 ; f max ]. Cela permet de travailler dans la plage de fréquences où le rapport signal sur bruit est bon et d’exclure la plage de fréquences où ce rapport signal sur bruit est mauvais. En fin de compte, cela améliore la robustesse du procédé d’identification. D’autres modes de réalisation de l’étape 72 de réduction de la dimensionnalité sont possibles. Par exemple, d’autres caractéristiques que celles citées dans l’article Morizet2016 peuvent être utilisées en plus ou à la place des caractéristiques citées dans cet article. Il est aussi possible d’utiliser un plus petit ou un plus grand nombre de caractéristiques physiques. Il est aussi possible de réduire le nombre d’informations à traiter lors de l’étape de classification par d’autres méthodes que l’extraction de caractéristiques physiques. Par exemple, des méthodes mathématiques comme la méthode d’analyse en composantes principales peuvent être utilisées pour réaliser l’étape de réduction de dimensionnalité. Une brève description de cette méthode peut être trouvée dans la section 3.1 de l’article Pashmforoush2012. L’étape de réduction de dimensionnalité peut aussi être omise. Dans ce cas, aucune réduction du nombre d’échantillons du signal S e (w) ou S e (t) traité n’est réalisée avant l’étape de classification. D’autres méthodes de classification que celle précédemment décrites peuvent être mises en œuvre pour identifier la source acoustique. Par exemple, il est aussi possible d’utiliser la méthode KGA (« K-means Genetic Algorithm » en anglais) ou la méthode WPT (« Wavelet Packet Transform ») toutes les deux décrites dans l’article Pashmfouroush2012. La méthode espérance-maximisation connue sous l’acronyme EM (« expectation-maximization algorithm ») ou encore les méthodes de classification non-supervisées décrites dans l’article Sawan2015 peuvent aussi être utilisées. Des méthodes de classification supervisée peuvent aussi être mises en oeuvre. Par exemple, la méthode « Random Forest », décrite dans l’article Morizet2016 peut être implémentée. Les méthodes de classification supervisée suivantes sont aussi utilisables : k-NN (« k-Nearest Neighbors »), SVM (« Support Vector Machine »), un réseau de neurones artificiels et d’autres. La classification peut être réalisée à partir de l’estimation S e (w) seule. Dans ce cas, l’étape 70 de construction de l’estimation S e (t) peut être omise. A l’inverse, la classification peut aussi être faite seulement à partir du signal S e (t). Autres variantes : La structure n’est pas nécessairement un tube en matériau composite. Ce qui a été décrit ici s’applique à d’autres structures comme par exemple : - des plaques, par exemple en métal ou en matériau composite, - des structures en béton comme des piles de pont. Le microprocesseur 14 peut être un processeur générique, un processeur spécifique, un circuit intégré propre à une application (connu aussi sous le nom anglais d’ASIC pour « Application-Specific Integrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gate Array »). Plusieurs des variantes décrites ici peuvent être combinées dans un même mode de réalisation. CHAPITRE III : Avantages des modes de réalisation décrits Le fait d’utiliser directement l’estimation du signal acoustique S(t) produit au niveau de la position P s pour la classification à la place des signaux ultrasonores mesurés, permet de rendre la classification plus robuste et donc d’augmenter la fiabilité du procédé d’identification du signal acoustique. En particulier, cela permet de diminuer l’influence des distorsions du signal acoustique qui se produisent lorsqu’il se propage dans la structure pour atteindre les emplacements P j où il est mesuré par les capteurs C j . Grâce à cela, les procédés d’identification décrits ici sont applicables à toute structure et, en particulier, à des structures de grandes dimensions dans lesquelles le signal acoustique parcourt une distance importante avant d’atteindre les capteurs C j . Le fait d’utiliser les estimations G e,j (w) des fonctions G j (w) de propagation permet d’améliorer la robustesse de l’identification de la source acoustique sans pour cela nécessiter l’exécution d’une étape d’apprentissage lors de laquelle différentes fonctions G j,k (w) sont mesurées pour différentes positions P k possibles du défaut et pour un signal acoustique connu. De plus, puisque les estimations G e,j (w) sont construites à chaque fois que les signaux F j (t) sont mesurés, ces estimations G e,j (w) suivent et s’adaptent automatiquement aux variations des fonctions G j (w). En effet, les fonctions G j (w) peuvent varier dans le temps notamment en fonction de l’état de la structure et des conditions environnementales dans lesquelles est placée la structure. Le fait d’utiliser des poids W j dont les valeurs sont fonction de la distance qui sépare l’emplacement P j et la position P s , permet d’améliorer la précision de l’estimation S e (w) et donc d’améliorer la fiabilité du procédé décrit. De plus, dans ce cas, les valeurs des poids W j sont indépendantes des caractéristiques de la structure instrumentée. Ainsi, dans ce cas, le procédé d’identification fonctionne avec tout type de structure et pas seulement avec des tuyaux ou des structures planes comme des plaques. Le fait que la réponse, dans le domaine fréquentiel, de chaque capteur soit assimilée à une constante égale à un, simplifie la mise en œuvre du procédé d’identification. L’étape d’apprentissage des produits R j (w)G j,k (w) permet de construire une estimation précise du signal S e (w) et donc d’obtenir un procédé d’identification robuste. Procédé d’identification automatique d’une source acoustique à partir d’un signal acoustique S(t) produit par l’apparition ou l’évolution d’un défaut dans une structure, ce procédé comportant les étapes suivantes : a) l’instrumentation (50) de la structure en fixant sur cette structure au moins un capteur, puis b) la mesure (54), par le capteur, d'un signal ultrasonore F j (t) causé par le signal acoustique produit par l’apparition ou l’évolution du défaut dans la structure, ce signal ultrasonore mesuré étant relié, dans le domaine fréquentiel, au signal acoustique produit par la relation suivante : F j (w) = R j (w)G j (w)S(w), où : - j est un indice identifiant le capteur, - w est la pulsation en radian, - F j (w) est la transformée de Fourier du signal ultrasonore F j (t) mesuré dans le domaine temporel par le capteur j, - S(w) est la transformée de Fourier du signal acoustique S(t) produit par l’apparition ou l’évolution du défaut dans la structure, - R j (w) est la réponse, dans le domaine fréquentiel, du capteur j, - G j (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique dans la structure entre la position du défaut et l’emplacement du capteur j, puis c) la classification automatique (76) de la source acoustique dans une classe de sources acoustiques choisie parmi plusieurs classes de sources acoustiques possibles, caractérisé en ce que : - le procédé comporte également entre les étapes b) et c), une étape d) de reconstruction (58 ; 88) du signal acoustique produit lors de l’apparition ou l’évolution du défaut, cette étape de reconstruction comportant la construction d’une estimation S e (w) du signal produit au niveau de la position du défaut à partir du signal ultrasonore F j (t) mesuré par le capteur et à l’aide d’une estimation ou d’une mesure du produit R j (w)G j (w) qui relie, dans le domaine fréquentiel, le signal ultrasonore F j (w) mesuré au signal acoustique S(w) produit, et - la classification automatique est réalisée à partir de l’estimation du signal acoustique obtenue à l’issue de l’étape d). Procédé selon la revendication 1, dans lequel : - l’étape a) d’instrumentation comporte la fixation sur la structure d’au moins trois capteurs à des emplacements dont les coordonnées sont connues dans un repère fixe par rapport à la structure, - l’étape b) comporte la mesure (54), par chacun des capteurs, d'un signal ultrasonore F j (t) respectif causé par le signal acoustique produit par l’apparition ou l’évolution du défaut dans la structure, chacun des signaux ultrasonores mesurés étant relié, dans le domaine fréquentiel, au signal acoustique produit par la relation suivante : F j (w) = R j (w)G j (w)S(w), - le procédé comporte entre les étapes b) et d), une étape de localisation (56), dans le repère fixe, de la position du défaut qui a produit le signal acoustique à partir des mesures réalisées par les capteurs lors de l’étape b) et des coordonnées connues de ces capteurs dans ce repère fixe, et - l’étape d) de reconstruction comporte l’obtention de l’estimation ou de la mesure du produit R j (w)G j (w) à l’aide de la position du défaut obtenue à l’issue de l’étape de localisation. Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape d) de reconstruction comporte : - pour chaque capteur, une opération (60) de construction d’une estimation G e,j (w) de la fonction G j (w) de propagation du signal acoustique dans la structure entre la position où le défaut apparaît et l’emplacement de ce capteur, à l’aide de la relation suivante : où Γ est une phase définie par la relation suivante : où : - « arg » est la fonction qui retourne l’argument d’un nombre complexe, et - W j sont des poids prédéfinis, puis - la construction (62) de l’estimation S e (w) du signal S(w) à partir des termes (R j (w)G e,j (w))*F j (w) ou des termes F j (w)/(R j (w)G e,j (w)), où le symbole «(...)* » désigne le conjugué de la fonction complexe entre parenthèses. Procédé selon les revendications 2 et 3 prises ensemble, dans lequel chaque poids W j est défini par la relation suivante : où : - P s sont les coordonnées dans le repère fixe du défaut obtenues à l’issue de l’étape de localisation, - P j sont les coordonnées dans le repère fixe du capteur j, - c est la célérité à laquelle se propage le signal acoustique à l’intérieur de la structure, et - w est la pulsation en radian. Procédé selon l’une quelconque des revendications 2 à 4, dans lequel chaque fonction R j (w) est égale à un quelle que soit la valeur de la pulsation w. Procédé selon l’une quelconque des revendications 3 à 5, dans lequel la fonction G e,j (w) de propagation est une fonction de Green. Procédé selon la revendication 2, dans lequel : - avant l’étape d) de reconstruction, le procédé comporte une étape (80) d’apprentissage du produit R j (w)G j (w) des fonctions R j (w) et G j (w) pour différents emplacements possibles du défaut, cette étape d’apprentissage comportant les opérations suivantes : i) pour chaque capteur, la mesure ou l’estimation par simulation numérique d’un signal ultrasonore F j,k (t) mesuré par ce capteur lorsqu’un signal acoustique connu S c,k (t) est appliqué sur la structure en un emplacement P k dont les coordonnées sont connues dans le repère fixe, chacun des signaux ultrasonores F j,k (t) étant relié, dans le domaine fréquentiel, au signal acoustique connu S c,k (t) par la relation suivante : F j,k (w) = R j (w)G j,k (w)S c,k (w), où : - k est un indice identifiant l’emplacement P k dont les coordonnées sont connues, - F j,k (w) est la transformée de Fourier du signal ultrasonore F j,k (t) mesuré ou estimé dans le domaine temporelle par le capteur j lorsque le signal acoustique connu est appliqué à l’emplacement P k, - S c,k (w) est la transformée de Fourier du signal acoustique connu S c,k (t) appliqué à l’emplacement P k , - R j (w) est la réponse, dans le domaine fréquentiel, du capteur j, - G j,k (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique connu S c,k (t) dans la structure entre l’emplacement P k et l’emplacement du capteur j, puis ii) pour chaque capteur, la mémorisation, en association avec le capteur j et les coordonnées de l’emplacement P k , du produit R j (w)G j,k (w) des fonctions R j (w) et G j,k (w), le produit R j (w)G j,k (w) étant obtenu à l’aide de la relation suivante : R j (w)G j,k (w) = F j,k (w)/S c,k (w), iii) la répétition des opérations i) et ii) pour plusieurs emplacements P k possibles, - l’étape d) de reconstruction (88) comporte les opérations suivantes : iv) la sélection (90) des produits R j (w)G j,k (w) mémorisés en association avec l’emplacement P k le plus proche de la position du défaut obtenue à l’issue de l’étape de localisation, puis v) la détermination (92) de l’estimation S e (w) du signal S(w) à partir des termes (R j (w)G j,k (w))*F j (w) ou des termes F j (w)/(R j (w)G j,k (w)), où : - les produits R j (w)G j,k (w) sont ceux sélectionnés lors de l’opération iv), et - le symbole «(...)* » désigne le conjugué de la fonction complexe entre parenthèses. Procédé selon l’une quelconque des revendications précédentes, dans lequel : - le procédé comporte la construction de l’estimation S e (t), dans le domaine temporel, du signal acoustique S(t) en appliquant une transformation de Fourier inverse à l’estimation S e (w), puis - la classification automatique est réalisée à partir de l’estimation S e (t) dans le domaine temporel. Support d’enregistrement d’informations lisible par un microprocesseur, caractérisé en ce qu’il comporte des instructions non-transitoires pour l’exécution des étapes b) à d) d’un procédé d’identification conforme à l’une quelconque des revendications précédentes, lorsque ces instructions sont exécutées par le microprocesseur. Appareil d’identification automatique d’une source acoustique à partir d’un signal acoustique S(t) produit par l’apparition ou l’évolution d’un défaut dans une structure, cet appareil comportant : - au moins un capteur (C j ) fixé sur cette structure, ce capteur étant apte à mesurer un signal ultrasonore F j (t) causé par le signal acoustique produit par l’apparition ou l’évolution du défaut dans la structure, ce signal ultrasonore mesuré étant relié, dans le domaine fréquentiel, au signal acoustique produit par la relation suivante : F j (w) = R j (w)G j (w)S(w), où : - j est un indice identifiant le capteur, - w est la pulsation en radian, - F j (w) est la transformée de Fourier du signal ultrasonore F j (t) mesuré dans le domaine temporel par le capteur j, - S(w) est la transformée de Fourier du signal acoustique S(t) produit par l’apparition ou l’évolution du défaut dans la structure, - R j (w) est la réponse, dans le domaine fréquentiel, du capteur j, - G j (w) est la fonction, dans le domaine fréquentiel, de propagation du signal acoustique dans la structure entre la position où le défaut apparaît et l’emplacement du capteur j, - un calculateur électronique (10) configuré pour exécuter les étapes suivantes : - l’acquisition des mesures du capteur, puis - la classification automatique de la source acoustique dans une classe de sources acoustiques choisie parmi plusieurs classes de sources acoustiques possibles, caractérisé en ce que le calculateur électronique (10) est également configuré pour : - avant la classification automatique, exécuter une étape de reconstruction du signal acoustique produit lors de l’apparition ou l’évolution du défaut, cette étape de reconstruction comportant la construction d’une estimation du signal produit au niveau de la position du défaut à partir du signal ultrasonore F j (t) mesuré par le capteur et à l’aide d’une estimation ou d’une mesure du produit R j (w)G j (w) qui relie, dans le domaine fréquentiel, le signal ultrasonore F j (w) mesuré au signal acoustique S(w) produit, et - la classification automatique est réalisée à partir de l’estimation du signal acoustique obtenue à l’issue de l’étape de reconstruction.