FR3065079B1 - Procede et dispositif de sondage ultrasonore pour l'obtention de courbes de dispersion d'un milieu sonde - Google Patents
Procede et dispositif de sondage ultrasonore pour l'obtention de courbes de dispersion d'un milieu sonde Download PDFInfo
- Publication number
- FR3065079B1 FR3065079B1 FR1753103A FR1753103A FR3065079B1 FR 3065079 B1 FR3065079 B1 FR 3065079B1 FR 1753103 A FR1753103 A FR 1753103A FR 1753103 A FR1753103 A FR 1753103A FR 3065079 B1 FR3065079 B1 FR 3065079B1
- Authority
- FR
- France
- Prior art keywords
- interest
- ultrasonic
- matrix
- frequency
- signals
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000006185 dispersion Substances 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims abstract description 59
- 230000002123 temporal effect Effects 0.000 claims abstract description 44
- 230000005540 biological transmission Effects 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 4
- 239000000523 sample Substances 0.000 claims description 17
- 239000013598 vector Substances 0.000 claims description 17
- 230000001902 propagating effect Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 239000007787 solid Substances 0.000 claims description 10
- 238000002604 ultrasonography Methods 0.000 claims description 10
- 238000000926 separation method Methods 0.000 claims description 9
- 239000002033 PVDF binder Substances 0.000 claims description 7
- 229920002981 polyvinylidene fluoride Polymers 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 230000003247 decreasing effect Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 4
- 230000005693 optoelectronics Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 9
- 230000005284 excitation Effects 0.000 description 8
- 238000012360 testing method Methods 0.000 description 7
- 239000002131 composite material Substances 0.000 description 6
- 230000004044 response Effects 0.000 description 5
- 230000000739 chaotic effect Effects 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000004913 activation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000015654 memory Effects 0.000 description 2
- 235000019687 Lamb Nutrition 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/041—Analysing solids on the surface of the material, e.g. using Lamb, Rayleigh or shear waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/11—Analysing solids by measuring attenuation of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/32—Arrangements for suppressing undesired influences, e.g. temperature or pressure variations, compensating for signal noise
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4472—Mathematical theories or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/50—Processing the detected response signal, e.g. electronic circuits specially adapted therefor using auto-correlation techniques or cross-correlation techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/015—Attenuation, scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0231—Composite or layered materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02827—Elastic parameters, strength or force
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0427—Flexural waves, plate waves, e.g. Lamb waves, tuning fork, cantilever
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Ce procédé de sondage ultrasonore comporte les étapes suivantes : commande (100) d'une pluralité d'émetteurs pour M émissions successives d'ondes ultrasonores dans un milieu ; commande (100) de N récepteurs pour recevoir, pour chaque émission successive, N signaux de mesure ; obtention (100) d'une matrice [K(t)] de signaux temporels de taille MxN, chaque coefficient Ki,j(t) de cette matrice représentant le signal de mesure reçu par le j-ème récepteur dû à la i-ème émission ; transformation (102) de la matrice [K(t)] en une matrice [FTK(f)] de signaux fréquentiels sur une pluralité de fréquences temporelles d'intérêt. Pour chaque fréquence temporelle d'intérêt, on obtient (106) une matrice diagonalisable par autocorrélation de la matrice [FTK(f)] et on projette (112) plusieurs signaux de référence à plusieurs fréquences spatiales d'intérêt distinctes sur un sous-espace propre de la matrice diagonalisable obtenue. On obtient alors (114) une représentation bidimensionnelle de courbes de dispersion caractéristiques du milieu à partir de ces projections (112) à chaque fréquence temporelle d'intérêt et à chaque fréquence spatiale d'intérêt.
Description
La présente invention concerne un procédé et un dispositif de sondage ultrasonore d’un milieu pour en obtenir une représentation bidimensionnelle de courbes de dispersion relatives à des modes de propagation ultrasonore dans le milieu.
La propagation d’ondes mécaniques, notamment d’ondes ultrasonores, dans un milieu matériel est par nature dépendante de caractéristiques mécaniques du milieu, comme sa densité, sa géométrie, ses modules d’élasticité, etc. L’analyse de la propagation de ces ondes est donc étudiée dans la littérature scientifique depuis des décennies pour retrouver de façon non invasive et non destructive les propriétés d’un milieu de propagation, que ce soit dans le domaine médical, notamment en échographie, ou dans le domaine du contrôle non destructif de matériaux inertes.
Mais ces ondes présentent la particularité d’être souvent dispersives et multimodales. Cela signifie que leur vitesse de propagation dans le milieu sondé dépend de la fréquence et qu’à une fréquence donnée il existe plusieurs modes se propageant chacun à une vitesse donnée. C’est d’autant plus vrai que le milieu sondé est à structure mince devant la longueur d’onde utilisée puisqu’alors les ondes sont guidées par les surfaces de la structure, devenant ainsi des ondes guidées, voire des ondes de Lamb lorsque les interfaces de la structure sont sensiblement planes et parallèles. Cette particularité rend l’analyse de ces ondes très complexe.
Une approche généralement choisie est d’obtenir une représentation bidimensionnelle de courbes de dispersion représentatives de modes de propagation ultrasonore dans le milieu sondé. Ces courbes de dispersion sont obtenues pour un ensemble de fréquences temporelles d’intérêt, formant une première dimension relevant d’une transformation de Fourier de la dimension temporelle, et pour un ensemble de fréquences spatiales d’intérêt ou équivalent (nombre d’onde, vitesse de phase, vitesse de groupe, etc.), formant une deuxième dimension relevant d’une transformation de Fourier d’une dimension spatiale. Elles sont effectivement représentatives des caractéristiques mécaniques du milieu et il existe des méthodes pour retrouver ces caractéristiques à partir d’une représentation bidimensionnelle des courbes de dispersion dans les deux dimensions précitées, notamment par un calcul d’inversion basé sur un modèle théorique prédéterminé ou par comparaison. Ainsi, plus la représentation bidimensionnelle des courbes de dispersion est précise, plus ces méthodes sont efficaces.
Une représentation bidimensionnelle de courbes de dispersion peut par exemple être obtenue par un calcul de double FFT (de l’anglais « Fast Fourier
Transform») spatio-temporelle, comme enseigné dans l’article de Alleyne et al, intitulé « A two-dimensional Fourier transform method for the measurement of propagating multimode signais », publié dans Journal of the Acoustical Society of America, volume 89, n° 3, pages 1159 à 1168, mars 1991. Selon ce principe, un transducteur est commandé pour l’émission d’une onde ultrasonore impulsionnelle large bande dans le milieu. En réception, N transducteurs spatialement distincts, par exemple N éléments électroacoustiques d’un capteur multi-éléments, sont commandés de manière à recevoir simultanément et pendant une durée prédéterminée N signaux temporels de mesure des ondes se propageant dans le milieu selon plusieurs modes de propagation ultrasonore. Une première transformée de Fourier (temporelle) est appliquée par calcul de FFT indépendamment à chacun des N signaux temporels reçus sur une pluralité de fréquences temporelles d’intérêt. Ensuite, pour chaque fréquence temporelle f de cette pluralité de fréquences d’intérêt, une deuxième transformée de Fourier, cette fois-ci spatiale, est appliquée par calcul de FFT au signal constitué des N valeurs obtenues par les N premières transformées de Fourier à cette fréquence f, sur une pluralité de fréquences spatiales k d’intérêt. Il en résulte une représentation bidimensionnelle de cette double FFT dans l’espace des fréquences temporelles et spatiales, cette double FFT formant elle-même une représentation des courbes de dispersion relatives aux différents modes de propagation ultrasonore dans le milieu sondé. Le calcul est simple et rapide mais très bruité, notamment par l’apparition de lobes secondaires autour des maxima locaux. Ainsi les modes de propagation proches se recouvrent et sont difficiles à distinguer, surtout dans les cas complexes de structures composites en contrôle non destructif.
Une amélioration est apportée dans l’article de Ambrozinski et al, intitulé « Identification of material properties - efficient modelling approach based on guided wave propagation and spatial multiple signal classification », publié dans Structural Control and Health Monitoring, volume 22, pages 969 à 983, 2015. Comme précédemment, un transducteur est commandé pour l’émission d’une onde ultrasonore impulsionnelle dans le milieu. En réception, N transducteurs spatialement distincts, par exemple N éléments électroacoustiques d’un capteur multi-éléments ou N capteurs optoélectroniques d’un interféromètre laser, sont commandés de manière à recevoir simultanément et pendant une durée prédéterminée N signaux temporels de mesure des ondes se propageant dans le milieu selon plusieurs modes de propagation ultrasonore. Une première transformée de Fourier (temporelle) est appliquée par calcul de FFT indépendamment à chacun des N signaux temporels reçus sur une pluralité de fréquences temporelles d’intérêt. Ensuite, selon ce nouvel enseignement, pour chaque fréquence temporelle f de cette pluralité de fréquences d’intérêt, ce n’est pas une deuxième transformée de Fourier spatiale qui est appliquée, mais la projection de signaux de référence sur un sous-espace propre d’une autocorrélation du signal spatial à une dimension constitué des N valeurs obtenues par les N premières transformées de Fourier à cette fréquence temporelle f. Plus précisément, selon le principe de classification multiple de signaux MUSIC (de l’anglais « MUItiple Signal Classification »), le projecteur utilisé est un modèle paramétrique de pseudo-spectre spatial présentant des pôles à facteurs de qualité élevés. Il en résulte une représentation des courbes de dispersion nettement moins bruitée. Mais une certaine imprécision subsiste malgré tout.
Ces bruits et imprécisions peuvent en partie s’expliquer par le fait que les deux méthodes précitées de représentation des courbes de dispersion sont basées sur une seule émission d’onde ultrasonore impulsionnelle. Elles peuvent donc fournir des résultats satisfaisants dans le cas de milieux sondés à faible nombre de modes de propagation par ailleurs suffisamment énergétiques. En revanche, elles atteignent leurs limites dès que la complexité des modes de propagation augmente, incluant des modes de faible énergie.
Une solution est alors apportée par le brevet FR 2 946 753 B1, selon l’un de ses modes de réalisation. Elle consiste, comme l’objet de la présente invention, à proposer un procédé de sondage comportant les étapes suivantes : - commande d’une pluralité de transducteurs d’émission pour M émissions successives d’ondes ultrasonores dans le milieu, - commande de N transducteurs de réception de manière à recevoir pendant une durée prédéterminée, pour chaque émission successive, N signaux temporels de mesure des ondes se propageant dans le milieu selon plusieurs modes de propagation ultrasonore, - obtention d’une matrice [K(t)J de signaux temporels acoustiques de taille MxN, chaque coefficient Kjj(t) de cette matrice représentant le signal de mesure reçu par le j-ème transducteur de réception dû à la i-ème émission, - transformation de la matrice [K(t)j de signaux temporels en une matrice [FTK(f)] de signaux fréquentiels sur une pluralité de fréquences temporelles f d’intérêt, - à partir de la matrice [FTK(f)] de signaux fréquentiels, obtention de la représentation bidimensionnelle des courbes de dispersion.
Plus précisément, dans le document FR 2 946 753 B1, le principe de classification multiple de signaux MUSIC n’étant pas applicable à la matrice [FTK(f)] de signaux fréquentiels puisqu’elle n’est pas diagonalisable, pour chaque fréquence temporelle f d’intérêt il est proposé à la place une décomposition en valeurs singulières de cette matrice [FTK(f)], puis un filtrage par suppression des vecteurs singuliers de réception considérés comme liés à du bruit et la projection d’une base d’ondes planes de fréquences spatiales k différentes (ou de vitesses de phase différentes) sur les vecteurs singuliers de réception conservés. Cette méthode améliore nettement les deux précédentes mais contient encore du bruit affectant le calcul des courbes de dispersion et la séparation des modes de propagation visibles, en particulier lorsqu’elle s’applique au sondage ultrasonore de matériaux composites. L’article de Xu et al, intitulé « Sparse SVD method for high-resolution extraction of the dispersion curves of ultrasonic guided waves », publié dans IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, volume 63, n° 10, pages 1514 à 1524, octobre 2016, propose une amélioration profitant de la propriété parcimonieuse des modes de propagation mais au prix d’une complexité largement accrue. Elle impose notamment des paramètres heuristiques à ajuster par essais/erreurs pour chaque configuration étudiée.
Il peut ainsi être souhaité de prévoir un procédé de sondage ultrasonore pour l’obtention de courbes de dispersion qui permette de s’affranchir d’au moins une partie des problèmes et contraintes précités.
Il est donc proposé un procédé de sondage ultrasonore d’un milieu pour en obtenir une représentation bidimensionnelle de courbes de dispersion relatives à des modes de propagation ultrasonore dans le milieu, comportant les étapes suivantes : - commande d’une pluralité d’émetteurs ultrasonores pour M émissions successives d’ondes ultrasonores dans le milieu, - commande de N récepteur de manière à recevoir pendant une durée prédéterminée, pour chaque émission successive, N signaux temporels de mesure des ondes se propageant dans le milieu selon lesdits modes de propagation ultrasonore, - obtention d’une matrice [K(t)J de signaux temporels acoustiques de taille MxN, chaque coefficient Kjj(t) de cette matrice représentant le signal de mesure reçu par le j-ème récepteur dû à la i-ème émission, - transformation de la matrice [K(t)j de signaux temporels en une matrice [FTK(f)] de signaux fréquentiels sur une pluralité de fréquences temporelles f d’intérêt, - à partir de la matrice [FTK(f)] de signaux fréquentiels, obtention de la représentation bidimensionnelle, dans lequel l’obtention de la représentation bidimensionnelle à partir de la matrice [FTK(f)] de signaux fréquentiels comporte les étapes suivantes : - pour chaque fréquence temporelle f d’intérêt : • obtention d’une matrice diagonalisable par autocorrélation de la matrice [FTK(f)] de signaux fréquentiels, • projections de plusieurs signaux de référence à plusieurs fréquences spatiales k d’intérêt distinctes sur un sous-espace propre de la matrice diagonalisable obtenue, - obtention de la représentation bidimensionnelle des courbes de dispersion à partir desdites projections à chaque fréquence temporelle f d’intérêt et à chaque fréquence spatiale k d’intérêt.
En procédant ainsi, le bruit affectant le calcul des courbes de dispersion et la séparation des modes de propagation visibles est de façon surprenante nettement réduit sans ajout de complexité sensible. Il n’est en particulier pas nécessaire alors d’ajouter des paramètres heuristiques ou empiriques à ajuster pour améliorer la qualité des courbes de dispersion obtenues.
De façon optionnelle, pour chaque fréquence temporelle f d’intérêt : - la matrice obtenue par autocorrélation de la matrice [FTK(f)] de signaux fréquentiels est décomposée en valeurs propres et vecteurs propres, - les vecteurs propres liés aux modes de propagation sont séparés des vecteurs propres liés à du bruit selon un critère de séparation prédéterminé, et - le sous-espace propre sur lequel sont projetés les signaux de référence à fréquences spatiales k distinctes est constitué des vecteurs propres liés au bruit.
De façon optionnelle également, le critère de séparation entre les vecteurs propres liés aux modes de propagation et les vecteurs propres liés au bruit est la détection d’un maximum d’une dérivée seconde d’une courbe décroissante des valeurs propres.
De façon optionnelle également, les projections sont réalisées, selon le principe de classification multiple de signaux MUSIC, à l’aide d’un modèle paramétrique de pseudo-spectre spatial présentant autant de pôles que de vecteurs propres liés aux modes de propagation.
De façon optionnelle également, les projections sont réalisées à l’aide d’un projecteur défini de la façon suivante :
où vn(f) est le vecteur propre d’indice n à la fréquence temporelle d’intérêt f, où P est l’indice du dernier vecteur propre lié aux modes de propagation dans l’ordre décroissant des valeurs propres associées à la fréquence temporelle d’intérêt f, où TFs désigne la transformée de Fourier spatiale, exprimée dans la base des signaux spatiaux de référence ek(f) = exp(jkx), j étant l’imaginaire pur tel que j2 = -1, exp( ) étant la fonction exponentielle, TFsk désignant alors la k-ième composante de cette transformée de Fourier spatiale, et où « | | » est le module vectoriel.
De façon optionnelle également, les projections sont réalisées à l’aide d’un projecteur défini de la façon suivante :
où Àn(f) est la valeur propre d’indice n à la fréquence temporelle d’intérêt f, où vn(f) est le vecteur propre d’indice n à la fréquence temporelle d’intérêt f, où P est l’indice du dernier vecteur propre lié aux modes de propagation dans l’ordre décroissant des valeurs propres associées à la fréquence temporelle d’intérêt f, où ek = exp(jkx), j étant l’imaginaire pur tel que j2 = -1, exp( ) étant la fonction exponentielle, « H » étant le symbole hermitien de la transposée conjuguée vectorielle, et où « | | » est le module vectoriel.
De façon optionnelle également, une représentation bidimensionnelle de niveaux d’énergie des courbes de dispersion est en outre obtenue à chaque fréquence temporelle f d’intérêt et à chaque fréquence spatiale k d’intérêt à partir des valeurs et vecteurs propres liés aux modes de propagation.
Il est également proposé un dispositif de sondage à ultrasons d’un milieu pour en obtenir une représentation bidimensionnelle de courbes de dispersion relatives à des modes de propagation ultrasonore dans le milieu, comportant : - une sonde comprenant une pluralité d’émetteurs ultrasonores,
- une pluralité de récepteurs de signaux temporels de mesure des ondes se propageant dans le milieu selon lesdits modes de propagation ultrasonore, et - des moyens de commande des émetteurs et récepteurs, et de traitement conçus pour mettre en œuvre un procédé de sondage ultrasonore selon l’invention.
De façon optionnelle : - chaque émetteur ultrasonore est un transducteur électroacoustique comportant au moins l’un des éléments de l’ensemble constitué d’un transducteur piézoélectrique PZT, d’un transducteur micro-usiné de type cMUT, pMUT ou mMUT, d’un transducteur à film PVDF et d’un transducteur électromagnétique-acoustique, ou - la sonde comporte : • un objet solide réverbérant présentant une face de contact destinée à être disposée au contact du milieu à sonder, • au moins un transducteur électroacoustique disposé contre une paroi de l’objet solide réverbérant, et les émetteurs ultrasonores sont constitués d’un réseau de points de la face de contact à partir desquels sont émis des signaux ultrasonores résultant d’au moins un signal émis par ledit au moins un transducteur électroacoustique et reconstitué par focalisation ultrasonore adaptative dans l’objet solide réverbérant.
De façon optionnelle également, chaque récepteur est : - un transducteur électroacoustique comportant au moins l’un des éléments de l’ensemble constitué d’un transducteur piézoélectrique PZT, d’un transducteur micro-usiné de type cMUT, pMUT ou mMUT, d’un transducteur à film PVDF et d’un transducteur électromagnétique-acoustique, ou - un transducteur optoélectronique, par exemple conçu sous la forme d’un interféromètre laser. L’invention sera mieux comprise à l’aide de la description qui va suivre, donnée uniquement à titre d’exemple et faite en se référant aux dessins annexés dans lesquels : - la figure 1 représente schématiquement la structure générale d’un dispositif de sondage à ultrasons selon un mode de réalisation de l’invention, - la figure 2 illustre les étapes successives d’un procédé de sondage ultrasonore réalisé à l’aide du dispositif de la figure 1, selon un mode de réalisation de l’invention, - la figure 3 illustre, à l’aide d’un diagramme, une étape de sélection d’un sous-espace propre de projection du procédé de la figure 2, - les figures 4A, 4B et 4C illustrent trois exemples de représentations bidimensionnelles de courbes de dispersion obtenues avec (figure 4B) ou sans (figures 4A et 4C) application du procédé de la figure 2, - les figures 5A et 5B illustrent deux autres exemples de représentations bidimensionnelles de courbes de dispersion obtenues avec (figure 5B) ou sans (figure 5A) application du procédé de la figure 2, et - les figures 6A et 6B illustrent deux exemples de représentations bidimensionnelles de niveaux d’énergie de courbes de dispersion obtenues avec (figure 6B) ou sans (figure 6A) application du procédé de la figure 2. L’installation représentée schématiquement sur la figure 1 comporte un dispositif de sondage à ultrasons 10 couplé à un milieu 12 dont on souhaite obtenir des courbes de dispersion relatives à des modes de propagation ultrasonores pour identifier certaines de ses caractéristiques mécaniques. Le milieu 12 est par exemple une structure mince, c’est-à-dire une structure dont l’épaisseur représente une fraction des longueurs d’ondes moyennes émises par le dispositif 10. Dans ce cas, les ondes se propagent de façon particulièrement dispersive selon plusieurs modes de propagation guidés par les surfaces de la structure. Il peut s’agir d’une structure organique, par exemple une structure osseuse, pour une application médicale, ou d’une structure industrielle inerte pour une application de contrôle non destructif.
Le dispositif 10 comporte une sonde d’émission 14 comprenant une pluralité d’émetteurs ultrasonores Eb ..., E,, ..., EM. Ces émetteurs sont des transducteurs électroacoustiques qui transforment un signal d’excitation électrique en un signal d’excitation ultrasonore apte à se propager dans le volume intérieur du milieu 12. Ils peuvent par exemple être conçus sous la forme de transducteurs piézoélectriques PZT (pour Titano-Zirconate de Plomb), de transducteurs micro-usinés de type cMUT (de l’anglais «capacitive Micromachined Ultrasonic Transducer»), pMUT (de l’anglais « piezoelectric Micromachined Ultrasonic Transducer») ou mMUT (de l’anglais «magnetostrictive Micromachined Ultrasonic Transducer»), de transducteurs à film PVDF (de l’anglais « PolyVinyliDene Fluoride ») ou de transducteurs électromagnétiques-acoustiques. Ils peuvent être collés à une paroi du milieu 12. Ils peuvent aussi être simplement couplés à cette paroi à l’aide d’un gel de transmission ultrasonore.
Selon un autre mode de réalisation envisageable, la sonde 14 peut comporter : - un objet solide réverbérant, généralement qualifié de « cavité chaotique », présentant une face de contact destinée à être disposée au contact du milieu 12 à sonder, et - au moins un transducteur électroacoustique disposé contre une paroi de cet objet solide réverbérant.
Une telle sonde est ainsi proposée dans l’article de Montaldo et al, intitulé « Building three-dimensional images using a time-reversal chaotic cavity », publié dans IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, volume 52, n° 9, pages 1489 à 1497, septembre 2005. L’objet solide réverbérant y est implémenté sous la forme d’une cavité chaotique de forme parallélépipédique tronquée par un usinage sphérique.
Conformément au fonctionnement d’une telle sonde : - une phase d’apprentissage est nécessaire et inclut la mesure, en un point de la face de contact, d’une réponse ultrasonore à un signal d’excitation émis par ledit au moins un transducteur électroacoustique dans l’objet solide réverbérant, et - une phase de sondage du milieu 12, lorsque ce milieu est en contact avec le point considéré de la face de contact, fait suite à la phase d’apprentissage et comporte l’émission par ledit au moins un transducteur électroacoustique d’un signal reconstitué sur la base d’une focalisation ultrasonore adaptative de ladite réponse ultrasonore mesurée, de manière à obtenir l’émission de l’excitation ultrasonore souhaitée dans le milieu 12 à partir du point considéré de la face de contact.
Par linéarité spatiale des méthodes de focalisation ultrasonore adaptative, une émission multi-éléments peut être réalisée à partir de plusieurs points distincts de la surface de contact formant ainsi le réseau de plusieurs émetteurs ultrasonores Eb ..., E,, ..., Em. Une loi de retard peut même être appliquée à cette émission multi-éléments. Il suffit de mesurer plusieurs réponses ultrasonores respectivement aux différents points considérés. Ensuite, la focalisation ultrasonore adaptative peut être combinée pour les mesures issues de chacun des points considérés de sorte que l’émission du signal reconstitué par ledit au moins un transducteur électroacoustique en phase de sondage engendre les excitations ultrasonores souhaitées aux points considérés. L’utilisation d’une sonde d’émission à cavité chaotique réverbérante peut être avantageuse dans certaines applications où des émissions de forte amplitude sont souhaitées, notamment pour provoquer des non linéarités susceptibles de faire apparaître des modes de propagation initialement non excités par la structure. En particulier, la caractérisation d’interfaces entre deux milieux collés est plus sensible aux ondes non linéaires qu’aux ondes linéaires, ceci étant valable pour l’ensemble des non linéarités acoustiques de contact CAN (de l’anglais « Contact Acoustic Nonlinearities »).
Le dispositif 10 comporte en outre une pluralité 16 de récepteurs Rb ..., Rj, ..., Rn de signaux temporels de mesure des ondes se propageant dans le milieu 12 selon plusieurs modes de propagation ultrasonore à identifier. Chacun de ces récepteurs peut être un transducteur électroacoustique qui transforme une onde ultrasonore s’étant propagée dans le volume intérieur du milieu 12 en un signal électrique de mesure. Les récepteurs peuvent ainsi par exemple être conçus, comme les émetteurs, sous la forme de transducteurs piézoélectriques PZT (pour Titano-Zirconate de Plomb), de transducteurs micro-usinés de type cMUT (de l’anglais «capacitive Micromachined Ultrasonic Transducer »), pMUT (de l’anglais « piezoelectric Micromachined Ultrasonic Transducer») ou mMUT (de l’anglais « magnetostrictive Micromachined Ultrasonic Transducer»), de transducteurs à film PVDF (de l’anglais « PolyVinyliDene Fluoride ») ou de transducteurs électromagnétiques-acoustiques. Ils peuvent être collés à une paroi du milieu 12. Ils peuvent aussi être simplement couplés à cette paroi à l’aide d’un gel de transmission ultrasonore.
La pluralité 16 de récepteurs peut être conçue sous la forme d’une sonde ultrasonore comme la sonde d’émission 14. Une seule et même sonde peut aussi servir en émission et réception.
Selon un autre mode de réalisation envisageable, chacun des récepteurs Rb ..., Rj, ..., Rn peut être un transducteur optoélectronique, par exemple conçu sous la forme d’un interféromètre laser.
Le dispositif 10 comporte en outre une unité 18 de commande et de traitement interagissant avec la sonde d’émission 14 et la pluralité 16 de récepteurs. Cette unité 18 peut par exemple être mise en œuvre au moins partiellement dans un dispositif informatique tel qu’un ordinateur comportant un processeur associé à une ou plusieurs mémoires pour le stockage de fichiers de données et de programmes d’ordinateurs dont il exécute les instructions. L’unité 18 de commande et de traitement telle qu’illustrée sur la figure 1 comporte ainsi deux modules fonctionnels au moins en partie constitués d’un ou plusieurs programmes d’ordinateurs. Un premier module fonctionnel 20 de commande des émetteurs et récepteurs comporte un générateur 22 de signaux d’excitation électrique des émetteurs ultrasonores Eb Eh EM, un collecteur 24 de signaux électriques de mesure fournis par les récepteurs Rb Rj, ..., RN et un coordinateur 26 du générateur 22 et du collecteur 24 pour synchroniser les émissions et réceptions. Un deuxième module fonctionnel de traitement 28 comporte un transformateur de Fourier 30, un autocorrélateur 32, un projecteur 34 et un module 36 de reconstitution de courbes de dispersion. On notera par ailleurs que les modules fonctionnels précités de l’unité 18 de commande et de traitement, et leurs éléments constitutifs programmables, sont présentés comme distincts, mais cette distinction est purement fonctionnelle. Ils pourraient tout aussi bien être regroupés selon toutes les combinaisons possibles en un ou plusieurs logiciels. Leurs fonctions pourraient aussi être au moins en partie micro programmées ou micro câblées dans des circuits intégrés dédiés. Ainsi, en variante, le dispositif informatique mettant en œuvre l’unité 18 de commande et de traitement pourrait être remplacé par un dispositif électronique composé uniquement de circuits numériques et de mémoires ou registres (sans programme d’ordinateur) pour la réalisation des mêmes actions.
Conformément aux principes généraux de la présente invention : - le premier module fonctionnel de commande 20 permet une activation des émetteurs ultrasonores Eb ..., E,, ..., EM pour M émissions successives d’ondes ultrasonores dans le milieu 12, une activation des récepteurs Rb ..., Rj, ..., Rn de manière à recevoir pendant une durée prédéterminée, pour chaque émission successive, N signaux temporels de mesure des ondes se propageant dans le milieu 12 et l’obtention d’une matrice [K(t)j de signaux temporels acoustiques de taille MxN, chaque coefficient K,j(t) de cette matrice représentant le signal de mesure reçu par le récepteur Rj dû à la i-ème émission, - le transformateur de Fourier 30 a pour fonction de transformer la matrice [K(t)J de signaux temporels en une matrice [FTK(f)] de signaux fréquentiels sur une pluralité de fréquences temporelles f d’intérêt, - l’autocorrélateur 32 a pour fonction d’obtenir, pour chaque fréquence temporelle f d’intérêt, une matrice diagonalisable par autocorrélation de la matrice [FTK(f)] de signaux fréquentiels, - le projecteur 34 a pour fonction de projeter, pour chaque fréquence temporelle f d’intérêt, plusieurs signaux de référence à plusieurs fréquences spatiales k d’intérêt distinctes sur un sous-espace propre de la matrice diagonalisable obtenue par l’autocorrélateur 32, et - le module 36 de reconstitution de courbes de dispersion a pour fonction d’obtenir une représentation bidimensionnelle de courbes de dispersion à partir desdites projections à chaque fréquence temporelle f d’intérêt et à chaque fréquence spatiale k d’intérêt.
De façon optionnelle, le projecteur 34 peut aussi avantageusement remplir une fonction d’estimation d’une énergie des courbes de dispersion pour chaque fréquence temporelle f d’intérêt et pour chaque fréquence spatiale k d’intérêt.
Les étapes successives d’un procédé de sondage ultrasonore pouvant être réalisé à l’aide du dispositif de sondage ultrasonore 10 de la figure 1 vont maintenant être détaillées en référence à la figure 2.
Au cours d’une étape 100, le premier module fonctionnel de commande 20 commande les séquences d’émissions et de réceptions des émetteurs et récepteurs Et , ..., E,, ..., Em, Ri, , Rj, , Rn pour l’acquisition de la matrice [K(t)j.
Ces séquences sont au nombre de M, et comprennent chacune un tir d’émetteur. En variante et de façon connue en soi, la sonde d’émission 14 pourrait comporter plus de M émetteurs et plusieurs émetteurs adjacents pourraient être sollicités à chaque tir. Après chaque tir, les signaux sont reçus sur l’ensemble des N récepteurs, numérisés et transmis à l’unité 18 de commande et de traitement qui construit la matrice [K(t)J des réponses impulsionnelles, chaque coefficient Kj/t) de cette matrice représentant le signal de mesure reçu par le récepteur Rj en réponse à la i-ème émission, ce signal étant numérisé pour faciliter son traitement ultérieur.
On notera que l’activation de plusieurs émetteurs situés à des endroits différents au contact du milieu 12 permet d’enrichir les signaux enregistrés, y compris lorsque le signal émis par chaque émetteur est le même, parce que le trajet ultrasonore est différent dans le milieu 12 à chaque tir. Ceci est d’autant plus utile que le milieu est dispersif et multimodal. En pratique, on peut utiliser des impulsions gaussiennes et un minimum d’une dizaine de séquences de tirs. D’autres signaux peuvent être considérés, tels qu’une impulsion à distribution de Dirac, un « chirp » (terme anglais désignant un signal pseudopériodique modulé en fréquence autour d’une fréquence porteuse et également modulé en amplitude par une enveloppe dont les variations sont lentes par rapport aux oscillations de phase), ou, si l’on dispose d’une électronique de pilotage programmable disposant d’autant de voies qu’il y a d’émetteurs, un signal codé qui est ensuite décodé juste après réception par le récepteur Rj, les autres étapes restant inchangées.
De façon optionnelle, un filtrage peut être réalisé avec un masque temporel pour se débarrasser de certains éléments perturbant le calcul des courbes de dispersion, tels que le front incident ou l’onde de Rayleigh.
Au cours d’une étape 102, le transformateur de Fourier 30 du deuxième module fonctionnel de traitement 28 effectue une transformée de Fourier discrète de la matrice [K(t)J pour obtenir la matrice [FTK(f)] de signaux fréquentiels sur une pluralité de fréquences temporelles f d’intérêt qui dépend de la longueur temporelle des coefficients Kjj(t) et des bandes passantes des émetteurs et récepteurs. Chaque coefficient FTKjj(f) de la matrice [FTK(f)j représente en effet la transformée de Fourier du coefficient Kjj(t) correspondant de la matrice [K(t)J.
Le procédé passe ensuite à une boucle d’étapes 104, 106, 108, 110, 112 impliquant l’autocorrélateur 32 et le projecteur 34. Un indice est créé et géré à l’étape 104 pour le suivi des fréquences temporelles f d’intérêt. Le procédé sort de la boucle d’étapes dès qu’il est vérifié à l’étape 104 que toutes les fréquences temporelles f d’intérêt ont été traitées. Les étapes 106, 108, 110 et 112 sont exécutées pour chacune de ces fréquences temporelles.
Au cours de l’étape 106 exécutée par l’autocorrélateur 32 pour une fréquence temporelle f d’intérêt donnée, une matrice d’autocorrélation [Rxx(f)j est calculée à partir de la matrice [FTK(f)j. Elle est carrée et de taille NxN. Elle présente en outre la propriété, obtenue par le principe même du calcul d’autocorrélation, d’être diagonalisable.
Concrètement, [Rxx(f)j peut être définie comme la matrice des coefficients de corrélation du produit des moments de Pearson PPMCC (de l’anglais « Pearson Product-Moment Corrélation Coefficients »). Pour un coefficient normalisé d’indices (i, j) de cette matrice :
où Cj.j(f) correspond à la covariance en f de l’émetteur E, par rapport au récepteur Rj.
Cette covariance peut prendre la forme suivante :
où I et J sont des vecteurs correspondant respectivement à la ligne de la matrice [FTK(f)] d’indice i et à la colonne de la matrice [FTK(f)] d’indice j, et E[ ] est l’espérance mathématique.
Au cours de l’étape 108 exécutée indifféremment par l’autocorrélateur 32 ou le projecteur 34 pour une fréquence temporelle f d’intérêt donnée, une diagonalisation de la matrice [Rxx(f)j est exécutée. De façon bien connue, cette diagonalisation consiste à décomposer [Rxx(f)j en valeurs propres et vecteurs propres, c’est-à-dire à calculer les matrices [P(f)j et [D(f)J telles que :
où [P(f)j est la matrice des vecteurs propres Vj(f) et [D(f)J est la matrice diagonale des valeurs propres Àj(f) que l’on range par ordre décroissant, à la fréquence temporelle f.
Au cours de l’étape 110 exécutée indifféremment par l’autocorrélateur 32 ou le projecteur 34 pour une fréquence temporelle f d’intérêt donnée, les vecteurs et valeurs propres liés aux modes de propagation résultant des M émissions sont séparés des vecteurs et valeurs propres liés à du bruit selon un critère de séparation prédéterminé. Cette étape a pour objectif de sélectionner un sous-espace propre de projection. La figure 3 illustre une courbe décroissante de valeurs propres Àj(f) qui peuvent résulter d’une diagonalisation 108. On s’aperçoit que la décroissance suit principalement deux pentes successives. La première, forte, peut être considérée comme concernant le sous-espace des vecteurs et valeurs propres liés aux modes de propagation, alors que la deuxième, nettement plus faible, peut être considérée comme concernant le sous-espace des vecteurs et valeurs propres liés au bruit. Au vu de cette caractéristique, le critère de séparation entre les deux sous-espaces propres précités peut être multiple : il peut être lié à un seuil fixe de valeur normalisée de valeur propre ; il peut être lié à un nombre prédéterminé souhaité de vecteurs et valeurs propres liés aux modes de propagation à détecter ; il peut être lié à la détection d’une rupture de pente par régressions linéaires. D’autres critères de séparation peuvent également être envisagés. De nombreux critères sont enseignés dans la littérature scientifique. Avantageusement, un critère astucieux peut tirer profit de la rupture de pente précitée d’une façon optimale : il existe un point d’inflexion dans la courbe de décroissance des valeurs propres, à l’endroit de la rupture de pente. Ce point d’inflexion peut être déterminé par la détection d’un maximum d’une dérivée seconde de cette courbe.
A l’issue de l’étape 110, il existe un nombre entier P, 0<P<N, indiquant la dimension du sous-espace propre lié aux modes de propagation à la fréquence temporelle f, ou de façon équivalente le nombre de modes de propagation à la fréquence temporelle f si l’on considère que chaque vecteur propre représente un mode de propagation et un seul. La combinaison linéaire des N vecteurs et valeurs propres peut alors être séparée en deux sommes :
la première de ces deux sommes représentant les signaux utiles et la deuxième les signaux liés au bruit.
Au cours de l’étape 112 exécutée par le projecteur 34 pour une fréquence temporelle f d’intérêt donnée, plusieurs signaux de référence à plusieurs fréquences spatiales k d’intérêt distinctes sont projetés sur l’un des deux sous-espaces propres déterminés à l’étape précédente. Ces signaux de référence représentent par exemple plusieurs ondes planes de fréquences spatiales différentes.
Dans un mode de réalisation préféré de l’invention, le sous-espace propre de projection est constitué des vecteurs propres liés au bruit, c’est-à-dire les N-P vecteurs propres vP+1(f), ..., vN(f). Dans ce cas, les projections sont avantageusement réalisées, selon le principe de classification multiple de signaux MUSIC, à l’aide d’un modèle paramétrique de pseudo-spectre spatial présentant autant de pôles que de vecteurs propres liés aux modes de propagation, soit P pôles. A la fréquence temporelle f d’intérêt considérée, le projecteur peut ainsi s’écrire :
où TFs désigne la transformée de Fourier spatiale, exprimée dans la base des signaux spatiaux de référence ek = exp(jkx), j étant dans cette expression l’imagniaire pur tel que j2 = -1, x étant la dimension spatiale dans laquelle s’étendent les récepteurs Rb ..., RN dont la distance inter-éléments est contante, exp( ) étant la fonction exponentielle, TFsk désignant alors la k-ème composante de cette transformée de Fourier spatiale, et où « | | » est le module vectoriel. L’homme du métier saura adapter l’expression de cette projection au cas où la distance inter-éléments entre les récepteurs Rb ..., RN n’est pas constante. Dans ce cas, la base des signaux spatiaux de référence ek = exp(jkx) peut rester la même, les distances inter-éléments variables étant prises en compte dans les valeurs discrètes spatiales de la variable x. La projection ne se fait alors pas à proprement parler à
l’aide d’une transformée de Fourier spatiale selon la terminologie généralement admise.
En variante, le projecteur peut aussi s’écrire :
où « H » est le symbole hermitien de la transposée conjuguée vectorielle. Néanmoins, cette variante donne généralement de moins bons résultats que la précédente. Bien d’autres variantes à la portée de l’homme du métier peuvent aussi être envisagées et privilégiées selon les applications.
De façon optionnelle, le projecteur 34 peut, au cours de cette même étape 112 et pour une fréquence temporelle f d’intérêt donnée, estimer une valeur d’énergie E(f,k) aux différentes fréquences spatiales k d’intérêt précitées, ou plus précisément aux P fréquences spatiales ku, avec u e {1, ..., P}, auxquelles la projection PROJ(f,k) atteint les valeurs maximales liées à ses P pôles. On peut en effet considérer que E(f,k) est nulle en toute valeur de k sauf aux P fréquences spatiales kb ..., kP à partir du moment où l’on accepte l’hypothèse selon laquelle ces P fréquences spatiales sont celles des P modes de propagation à la fréquence temporelle f d’intérêt donnée. Un exemple d’estimation des P valeurs discrètes d’énergie Eu = E(f,ku) va maintenant être détaillé.
On sait tout d’abord que, compte tenu de l’expression de la diagonalisation de la matrice [Rxx(f)j :
Vn e {1, ...,?},[«%%(/)]. vn(f) = , chaque valeur propre !„(/), donc également chaque vecteur propre associé vn(f), correspondant à un mode de propagation.
On estime ensuite la variance du bruit dans le signal que l’on note ow et que l’on peut considérer comme égale à la dernière et plus petite des valeurs propres MO-
On en déduit l’ensemble de relations suivantes :
où eu = [1, expQ'kuX^,..., exp(jkuxN)] et x^ ..., xN sont les coordonnées des N récepteurs R^ ..., RN dans la dimension spatiale d’intérêt.
Le système d’équations précédent peut s’écrire sous forme matricielle :
D’où le résultat recherché, par simple inversion :
Lorsqu’il a été constaté à l’étape 104 que l’ensemble des fréquences temporelles f d’intérêt a été exploré, le procédé passe à une dernière étape 114, exécutée par le module 36 de reconstitution de courbes de dispersion, pour l’obtention d’une représentation bidimensionnelle de courbes de dispersion à partir des projections précitées estimées à chaque fréquence temporelle f d’intérêt et à chaque fréquence spatiale k d’intérêt. Cette représentation peut être tout simplement la représentation en niveaux de gris (en dB par exemple) de la fonction PROJ(f,k). Mais il convient de noter que les fréquences spatiales k d’intérêt, autrement appelées généralement « nombres d’ondes », ne sont pas les seuls paramètres physiques spatiaux selon lesquels les courbes de dispersion peuvent être représentées. Les vitesses de phase, ou même les vitesses de groupe, peuvent être préférées dans certains contextes. Le paramètre de vitesse de phase est par exemple mathématiquement directement lié par la pulsation au nombre d’onde, de sorte que les multiples représentations possibles restent équivalentes tout en étant différentes.
De façon optionnelle, le module 36 de reconstitution de courbes de dispersion fournit également à l’étape 114 une représentation bidimensionnelle de niveaux d’énergie des courbes de dispersion à l’aide des valeurs d’énergies E(f,k) éventuellement calculées à l’étape 112.
Les figures 4A, 4B et 4C illustrent des résultats obtenus par des tests expérimentaux. Ces tests ont été menés avec une sonde émettrice et réceptrice à 128 éléments de fréquence centrale égale à 20 MHz, dont 25 éléments sont utilisés en émission. La figure 4A montre le résultat obtenu à l’aide du procédé enseigné dans l’article cité précédemment de Ambrozinski et al, en choisissant la base des signaux spatiaux de référence ek = exp(jkx) précités pour réaliser la projection MUSIC. La figure 4B montre le résultat obtenu à l’aide du procédé selon la présente invention décrit précédemment. La figure 4C montre le résultat obtenu à l’aide du procédé enseigné dans le document de brevet FR 2 946 753 B1. On voit que les nombreux modes de propagation sont bien plus faciles à distinguer et identifier sur la figure 4B que sur les deux autres. En particulier, les courbes de dispersions visibles
ont quasiment toutes une épaisseur d’un seul pixel sur la figure 4B, ce qui témoigne d’une excellente définition.
Les figures 5A et 5B illustrent d’autres résultats obtenus par des tests expérimentaux. Ces tests ont été menés avec une sonde émettrice et réceptrice à 128 éléments de fréquence centrale égale à 10 MHz, dont 15 éléments sont utilisés en émission. Mais le milieu 12 observé est un assemblage collé de type composite sur composite. En effet, les matériaux composites ne sont généralement pas simples à étudier car ils sont très bruités et l’énergie des modes de propagation se dissipe rapidement. La figure 5A montre le résultat obtenu à l’aide du procédé enseigné dans le document de brevet FR 2 946 753 B1. La figure 5B montre le résultat obtenu à l’aide du procédé selon la présente invention décrit précédemment. On voit que les nombreux modes de propagation sont encore une fois bien plus faciles à distinguer et identifier sur la figure 5B que sur la figure 5A.
Les figures 6A et 6B illustrent d’autres résultats obtenus par des tests expérimentaux. Ces tests ont été également menés avec une sonde émettrice et réceptrice à 128 éléments de fréquence centrale égale à 10 MHz, dont 15 éléments sont utilisés en émission. Mais dans ces résultats, ce sont les niveaux d’énergie à chaque fréquence temporelle f et à chaque fréquence spatiale k qui sont calculés et représentés. La figure 6A montre le résultat obtenu à l’aide du procédé enseigné dans l’article de Alleyne et al précité. La figure 6B montre le résultat obtenu à l’aide du procédé selon la présente invention décrit précédemment. On voit que les modes de propagation sont mieux discriminés sur la figure 6B que sur la figure 6A, de sorte que les éventuels transferts d’énergie d’un mode de propagation à l’autre y sont plus clairement détectables. En particulier, le principe de calcul présenté précédemment, pour estimer les niveaux d’énergie uniquement en des valeurs discrètes de la représentation bidimensionnelle correspondant pour chaque fréquence temporelle f d’intérêt aux modes de propagation détectés, permet de conserver l’excellente résolution des courbes de dispersion obtenues. On dispose alors d’une estimation de l’énergie par mode, ce qui permet de vérifier la répartition énergétique des modes réellement excités dans le milieu sondé en fonction de divers paramètres liés à l’excitation ultrasonore, comme par exemple l’amplitude de l’excitation, ou des conditions aux limites éventuelles.
Il apparaît clairement qu’un dispositif de sondage à ultrasons tel que celui décrit précédemment permet d’obtenir une représentation de courbes de dispersion de qualité améliorée par rapport aux méthodes connues. Cette amélioration est d’autant plus visible que les modes de propagation sont nombreux et difficiles à différencier, en particulier lorsque le milieu sondé est un matériau composite.
On notera par ailleurs que l’invention n’est pas limitée aux modes de réalisation décrits précédemment. Il apparaîtra en effet à l'homme de l'art que diverses modifications peuvent être apportées aux modes de réalisation décrits ci-dessus, à la lumière de l'enseignement qui vient de lui être divulgué. Dans les revendications qui suivent, les termes utilisés ne doivent pas être interprétés comme limitant les revendications aux modes de réalisation exposés dans la présente description, mais doivent être interprétés pour y inclure tous les équivalents que les revendications visent à couvrir du fait de leur formulation et dont la prévision est à la portée de l'homme de l'art en appliquant ses connaissances générales à la mise en œuvre de l'enseignement qui vient de lui être divulgué.
Claims (6)
- REVENDICATIONS1. Procédé de sondage ultrasonore d’un milieu (12) pour en obtenir une représentation bidimensionnelle de courbes de dispersion relatives à des modes de propagation ultrasonore dans le milieu, comportant les étapes suivantes : commande (100) d’une pluralité (14) d’émetteurs ultrasonores (Eb ..., Ej, ..., Em) pour M émissions successives d’ondes ultrasonores dans le milieu (12), commande (100) de N récepteurs (Rb ..., Rj, ..., RN) de manière à recevoir pendant une durée prédéterminée, pour chaque émission successive, N signaux temporels de mesure des ondes se propageant dans le milieu (12) selon lesdits modes de propagation ultrasonore, obtention (100) d’une matrice [K(t)J de signaux temporels acoustiques de taille MxN, chaque coefficient Kjj(t) de cette matrice représentant le signal de mesure reçu par le j-ème récepteur (Rj) dû à la i-ème émission, transformation (102) de la matrice [K(t)j de signaux temporels en une matrice [FTK(f)] de signaux fréquentiels sur une pluralité de fréquences temporelles f d’intérêt, à partir de la matrice [FTK(f)] de signaux fréquentiels, obtention de la représentation bidimensionnelle, caractérisé en ce que l’obtention de la représentation bidimensionnelle à partir de la matrice [FTK(f)] de signaux fréquentiels comporte les étapes suivantes : pour chaque fréquence temporelle f d’intérêt : • obtention (106) d’une matrice diagonalisable par autocorrélation de la matrice [FTK(f)] de signaux fréquentiels, • projections (112) de plusieurs signaux de référence à plusieurs fréquences spatiales k d’intérêt distinctes sur un sous-espace propre de la matrice diagonalisable obtenue, obtention (114) de la représentation bidimensionnelle des courbes de dispersion à partir desdites projections (112) à chaque fréquence temporelle f d’intérêt et à chaque fréquence spatiale k d’intérêt.
- 2. Procédé de sondage ultrasonore selon la revendication 1, dans lequel, pour chaque fréquence temporelle f d’intérêt : la matrice obtenue (106) par autocorrélation de la matrice [FTK(f)] de signaux fréquentiels est décomposée (108) en valeurs propres et vecteurs propres, les vecteurs propres liés aux modes de propagation sont séparés (110) des vecteurs propres liés à du bruit selon un critère de séparation prédéterminé, et le sous-espace propre sur lequel sont projetés (112) les signaux de référence à fréquences spatiales k distinctes est constitué des vecteurs propres liés au bruit.
- 3. Procédé de sondage ultrasonore selon la revendication 2, dans lequel le critère de séparation (110) entre les vecteurs propres liés aux modes de propagation et les vecteurs propres liés au bruit est la détection d’un maximum d’une dérivée seconde d’une courbe décroissante des valeurs propres. 4. Procédé de sondage ultrasonore selon la revendication 2 ou 3, dans lequel les projections (112) sont réalisées, selon le principe de classification multiple de signaux MUSIC, à l’aide d’un modèle paramétrique de pseudo-spectre spatial présentant autant de pôles que de vecteurs propres liés aux modes de propagation. 5. Procédé de sondage ultrasonore selon la revendication 4, dans lequel les projections (112) sont réalisées à l’aide d’un projecteur défini de la façon suivante :où vn(f) est le vecteur propre d’indice n à la fréquence temporelle d’intérêt f, où P est l’indice du dernier vecteur propre lié aux modes de propagation dans l’ordre décroissant des valeurs propres associées à la fréquence temporelle d’intérêt f, où TFs désigne la transformée de Fourier spatiale, exprimée dans la base des signaux spatiaux de référence ek = exp(jkx), j étant l’imaginaire pur tel que j2 = -1, x étant la dimension spatiale dans laquelle s’étendent les récepteurs (F^, ..., Rj, ..., RN), exp( ) étant la fonction exponentielle, TFsk désignant alors la k-ième composante de cette transformée de Fourier spatiale, et où « | | » est le module vectoriel. 6. Procédé de sondage ultrasonore selon la revendication 4, dans lequel les projections (112) sont réalisées à l’aide d’un projecteur défini de la façon suivante :où Àn(f) est la valeur propre d’indice n à la fréquence temporelle d’intérêt f, où vn(f) est le vecteur propre d’indice n à la fréquence temporelle d’intérêt f, où P est l’indice dudernier vecteur propre lié aux modes de propagation dans l’ordre décroissant des valeurs propres associées à la fréquence temporelle d’intérêt f, où ek = exp(jkx), j étant l’imaginaire pur tel que j2 = -1, x étant la dimension spatiale dans laquelle s’étendent les récepteurs (R^ Rj, RN), exp( ) étant la fonction exponentielle, « H » étant le symbole hermitien de la transposée conjuguée vectorielle, et où « | | » est le module vectoriel.
- 7. Procédé de sondage ultrasonore selon l’une quelconque des revendications 2 à 6, dans lequel une représentation bidimensionnelle de niveaux d’énergie des courbes de dispersion est en outre obtenue (114) à chaque fréquence temporelle f d’intérêt et à chaque fréquence spatiale k d’intérêt à partir des valeurs et vecteurs propres liés aux modes de propagation. 8. Dispositif (10) de sondage à ultrasons d’un milieu (12) pour en obtenir une représentation bidimensionnelle de courbes de dispersion relatives à des modes de propagation ultrasonore dans le milieu, comportant : une sonde (14) comprenant une pluralité d’émetteurs ultrasonores (Ei, ..., Ej, ..., Em), une pluralité de récepteurs (R^ ..., Rj, ..., RN) de signaux temporels de mesure des ondes se propageant dans le milieu (12) selon lesdits modes de propagation ultrasonore, et des moyens (18) de commande (20) des émetteurs (E^ ..., E,, ..., EM) et récepteurs (R^ ..., Rj, ..., RN), et de traitement (28) conçus pour mettre en œuvre un procédé de sondage ultrasonore selon l’une quelconque des revendications 1 à 7.
- 9. Dispositif (10) de sondage à ultrasons selon la revendication 8, dans lequel : chaque émetteur ultrasonore (Eb ..., E,, ..., EM) est un transducteur électroacoustique comportant au moins l’un des éléments de l’ensemble constitué d’un transducteur piézoélectrique PZT, d’un transducteur micro-usiné de type cMUT, pMUT ou mMUT, d’un transducteur à film PVDF et d’un transducteur électromagnétique-acoustique, ou la sonde (14) comporte : • un objet solide réverbérant présentant une face de contact destinée à être disposée au contact du milieu (12) à sonder, • au moins un transducteur électroacoustique disposé contre une paroi de l’objet solide réverbérant, et les émetteurs ultrasonores (E^ ..., E,, ..., EM) sont constitués d’un réseau de points de la face de contact à partir desquels sont émis des signaux ultrasonores résultant d’au moins un signal émis par ledit au moins un transducteur électroacoustique et reconstitué par focalisation ultrasonore adaptative dans l’objet solide réverbérant.
- 10. Dispositif (10) de sondage à ultrasons selon la revendication 8 ou 9, dans lequel chaque récepteur (Rb ..., Rj, ..., RN) est : un transducteur électroacoustique comportant au moins l’un des éléments de l’ensemble constitué d’un transducteur piézoélectrique PZT, d’un transducteur micro-usiné de type cMUT, pMUT ou mMUT, d’un transducteur à film PVDF et d’un transducteur électromagnétique-acoustique, ou un transducteur optoélectronique, par exemple conçu sous la forme d’un interféromètre laser.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1753103A FR3065079B1 (fr) | 2017-04-10 | 2017-04-10 | Procede et dispositif de sondage ultrasonore pour l'obtention de courbes de dispersion d'un milieu sonde |
PCT/FR2018/050822 WO2018189451A1 (fr) | 2017-04-10 | 2018-04-03 | Procédé et dispositif de sondage ultrasonore pour l'obtention de courbes de dispersion d'un milieu sonde |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1753103 | 2017-04-10 | ||
FR1753103A FR3065079B1 (fr) | 2017-04-10 | 2017-04-10 | Procede et dispositif de sondage ultrasonore pour l'obtention de courbes de dispersion d'un milieu sonde |
Publications (2)
Publication Number | Publication Date |
---|---|
FR3065079A1 FR3065079A1 (fr) | 2018-10-12 |
FR3065079B1 true FR3065079B1 (fr) | 2019-06-21 |
Family
ID=59031182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1753103A Active FR3065079B1 (fr) | 2017-04-10 | 2017-04-10 | Procede et dispositif de sondage ultrasonore pour l'obtention de courbes de dispersion d'un milieu sonde |
Country Status (2)
Country | Link |
---|---|
FR (1) | FR3065079B1 (fr) |
WO (1) | WO2018189451A1 (fr) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117031540B (zh) * | 2023-08-01 | 2024-02-13 | 北京阳光吉澳能源技术有限公司 | 一种微动探测数据处理方法、***、装置和介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2858099B1 (fr) * | 2003-07-25 | 2006-03-24 | Centre Nat Rech Scient | Procede et dispositif de focalisation d'ondes acoustiques |
FR2946753B1 (fr) | 2009-06-11 | 2011-07-22 | Centre Nat Rech Scient | Procede et dispositif ultrasonores pour caracteriser un milieu |
-
2017
- 2017-04-10 FR FR1753103A patent/FR3065079B1/fr active Active
-
2018
- 2018-04-03 WO PCT/FR2018/050822 patent/WO2018189451A1/fr active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2018189451A1 (fr) | 2018-10-18 |
FR3065079A1 (fr) | 2018-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2878545C (fr) | Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants | |
EP2294452B1 (fr) | Procede et dispositif de sondage par propagation d'ondes | |
EP2440140B1 (fr) | Procede et dispositif ultrasonores pour caracteriser un milieu | |
EP1531733B1 (fr) | Dispositif et procede pour la mesure de l'elasticite d'un organe humain ou animal | |
EP2294401B1 (fr) | Dispositif acoustique de mesure localisée et sans contact des non-linearites élastique et dissipative et de la viscoélasticité | |
Larose et al. | Reconstruction of Rayleigh–Lamb dispersion spectrum based on noise obtained from an air-jet forcing | |
WO2019091705A1 (fr) | Contrôle de sante d'une structure industrielle | |
FR2883982A1 (fr) | Procede et dispositif d'imagerie utilisant des ondes de cisaillement | |
EP1169636A1 (fr) | Procede et dispositif d'imagerie utilisant les ondes de cisaillement | |
EP0751743A1 (fr) | Procede et dispositif d'evaluation et de caracterisation des proprietes des os | |
EP3084416B1 (fr) | Procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage à ultrasons correspondants | |
FR3065079B1 (fr) | Procede et dispositif de sondage ultrasonore pour l'obtention de courbes de dispersion d'un milieu sonde | |
Weaver et al. | The mean and variance of diffuse field correlations in finite bodies | |
FR3065078B1 (fr) | Procede et dispositif de sondage ultrasonore par focalisation adaptative au moyen d'un objet solide reverberant | |
EP1649449B1 (fr) | Procédé et dispositif d'imagerie par ondes acoustiques | |
EP4083659B1 (fr) | Dispositif et procédé de traitement du signal issus d'un ensemble de transducteurs ultrasonores | |
CA3046106A1 (fr) | Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants | |
WO2017202628A1 (fr) | Procede et dispositif de caracterisation de la peau | |
WO2011124822A2 (fr) | Procédé et système d'échographie ultrasonore anharmonique |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20181012 |
|
PLFP | Fee payment |
Year of fee payment: 3 |
|
PLFP | Fee payment |
Year of fee payment: 4 |
|
PLFP | Fee payment |
Year of fee payment: 5 |
|
PLFP | Fee payment |
Year of fee payment: 6 |
|
PLFP | Fee payment |
Year of fee payment: 7 |
|
PLFP | Fee payment |
Year of fee payment: 8 |