FR 2932887 A1 20091225 FR 0803542 A 20080624 DISPOSITIF ACOUSTIQUE DE MESURE LOCALISEE ET SANS CONTACT DES NON-LINEARITES ELASTIQUE ET DISSIPATIVE ET DE LA VISCOELASTICITE DOMAINE TECHNIQUE DE L'INVENTION [0001] L'invention se rapporte à un procédé et dispositif acoustique de mesure localisée, sans contact, des non-linéarités élastique et dissipative et de la viscoélasticité dans un milieu soumis à une contrainte acoustique basse fréquence et sondé par des impulsions ou trains d'ondes ultrasonores. [0002] La présente invention concerne le domaine technique de la détection et la quantification du micro-endommagement (micro-fissures) dans des matériaux industriels ou biologiques, notamment le tissu osseux. [0003] L'invention concerne également le domaine technique de la caractérisation d'autres milieux complexes tels que les milieux granulaires consolidés ou non-15 consolidés (roches, sable, béton, ...). ETAT DE LA TECHNIQUE ANTERIEURE [0004] Il est connu de l'état de la technique plusieurs documents décrivant des systèmes de mesure de non-linéarités acoustiques. Ces publications décrivent des méthodes permettant de mesurer les effets non linéaires acoustiques en 20 soumettant le matériau testé à des contraintes mécaniques d'amplitude variable. En effet, l'application d'une contrainte mécanique variable au matériau testé modifie ses propriétés élastiques, et parfois dissipatives. [0005] Une partie de ces méthodes mesure les variations de la vitesse de propagation et de l'amplitude des ondes acoustiques dans un échantillon, induites 25 par une variation quasi-statique de l'état de contrainte mécanique (pression hydrostatique ou contrainte uniaxiale) de cet échantillon. L'autre partie de ces méthodes utilisent l'interaction de deux ondes acoustiques : N. Ichida, T. Sato, M. Linzer; Imaging the nonlinear parameter of the medium, Ultrason. Imaging, Vol. 5, 1983 ; pp. 295-299. Z. Zhu, M.S. Roos, W.N. Cobb, K. Jensen; Determination of the acoustic nonlinearity parameter B/A from phase measurements; J. Acoust. Soc. Am., Vol. 74(5), 1983; pp.1518-1521 ; C. Barrière, D. Royer; Diffraction effects in the parametric interaction of acoustic waves: application to measurements of the nonlinearity parameter B/A in liquids; IEEE Trans. UFFC., Vol. 48 (6), 2001, pp 1706-1715 ; C.M. Sehgal, R.C. Bahn, J.F. Greenleaf; Measurement of the acoustic nonlinearity parameter B/A in human tissues by a thermodynamic method; J. Acoust. Am. Soc., Vol. 76(4), 1984; pp.1023-1029 ; V. Zaitsev, V. Nazarov, V. Gusev et B. Castagnede; Novel nonlinearmodulation acoustic technique for crack detection; NDT&E International, 39 (2006) 184-194; G. Gremaud, M. Bujard, et W. Benoît ; The coupling technique: A two-wave acoustic method for the study of dislocation dynamics ; J. Appl. Phys. 1987 61(5), 1795-1805. [0006] Cependant, de tels systèmes de mesure présentent plusieurs inconvénients. [0007] Les techniques utilisant une machine de test quasi-statique (nécessitant une mise en contact de l'échantillon avec le dispositif) ne sont pas applicables aux milieux fluides, aux gels et aux milieux vivants (applications biomédicales in vivo). [0008] Les méthodes en milieu sec nécessitant un contact (collage par exemple) des transducteurs ultrasonores avec les échantillons, et donc altérant leurs surfaces, ne sont pas compatibles avec certaines applications, notamment les applications biomédicales. [0009] Les méthodes employant une onde sonde ultrasonore émise en continue peuvent générer des ondes stationnaires complexifiant l'analyse des mesures. [0010] Dans le cas d'une interaction colinéaire entre l'onde basse fréquence et l'onde ultrasonore, il peut exister des zones d'ombre acoustique dans le champ basse-fréquence créé par les transducteurs ultrasonores. [0011] Les méthodes employant une interaction colinéaire entre l'onde basse fréquence et l'onde ultrasonore, et où l'onde ultrasonore doit faire un ou plusieurs allers-retours dans l'échantillon, ne sont applicables aux milieux fortement atténuants, tel que le tissu osseux. [0012] Les méthodes employant un dispositif mécanique en contact pour faire varier les contraintes dans l'échantillon, sans utiliser une résonance mécanique de cet échantillon, ont une gamme d'utilisation limitée en termes de vitesse de déformation (au maximum quelques kHz). [0013] Les méthodes employant un dispositif mécanique en contact pour faire varier les contraintes dans l'échantillon, utilisant une résonance mécanique de cet échantillon, peuvent atteindre des fréquences de déformation de quelques dizaines de kHz, cependant la fréquence de résonance mécanique est conditionnée par les dimensions géométriques et les propriétés physiques de l'échantillon. [0014] Les méthodes faisant varier la pression hydrostatique par paliers et nécessitant d'enfermer l'échantillon testé dans une enceinte hermétique remplie d'un liquide sont de conception et d'usage relativement complexes. [0015] Les méthodes n'utilisant pas d'onde acoustique basse fréquence pour générer les variations de contrainte permettent difficilement d'effectuer un test à la fois en compression et en traction (ou compression et détente). Dans le cas des machines de test, le dispositif de maintien altère l'échantillon. Dans le cas des enceintes hermétiques, une diminution de la pression hydrostatique inférieure à la pression atmosphérique est de conception relativement complexe. EXPOSE DE L'INVENTION [0016] La présente invention vise à pallier les inconvénients de l'état de la technique en proposant une méthode et un dispositif sans contact de mesure localisée et instantanée des effets viscoélastiques et des variations non-linéaires élastique et dissipative d'un échantillon soumis à des contraintes extérieures. [0017]Selon un autre aspect, la présente invention permet également d'envisager la mesure des non-linéarités élastique et dissipative et de la viscoélasticité dans des domaines d'applications où le contact avec le système de mesure est non souhaitable tels que l'agro-alimentaire et les applications biomédicales. [0018] L'invention vise par ailleurs à mesurer ces variations instantanées d'élasticité et d'atténuation en fonction de la variation instantanée de la contrainte extérieure, pendant les phases de compression et de traction de l'échantillon, ainsi qu'à quantifier le coefficient de viscoélasticité de l'échantillon analysé. [0019] Pour ce faire, les contraintes extérieures appliquées à l'échantillon sont générées par une onde acoustique basse fréquence dont la longueur d'onde est nettement supérieure à la taille de l'échantillon dans le milieu où se trouve l'échantillon. Les variations de pression hydrostatique ainsi générées compriment et détendent l'échantillon dans les trois directions de l'espace. Simultanément, une série d'impulsions ultrasonores (hautes fréquences) est transmise à travers l'échantillon afin de mesurer les variations de vitesse de propagation et d'atténuation du milieu associées aux variations de la pression hydrostatique. [0020] Plus précisément, l'invention a pour objet un dispositif acoustique de mesure, localisée et sans contact, des micro-endommagements d'un échantillon disposé dans un réservoir, dans lequel des moyens d'émission d'une onde acoustique basse fréquence créent une variation périodique basse fréquence de la pression hydrostatique dans l'échantillon. Ce dispositif comporte également des moyens de génération d'impulsions ultrasonores haute fréquence et une unité d'analyse comportant un module de calcul des modulations du temps de vol et de l'amplitude des impulsions ultrasonores lors du passage de l'onde acoustique basse fréquence dans l'échantillon, de sorte à caractériser les effets non-linéaire élastiques et dissipatifs ainsi que la viscoélasticité de l'échantillon. [0021] Selon des modes de réalisation particuliers, le dispositif comporte : [0022] - un module d'imagerie multi-paramétrique des effets viscoélastiques (tan()) et des effets non-linéaires acoustiques élastique (a, (3 et b) et dissipatif ; [0023] - des moyens de génération de l'onde acoustique basse fréquence comportent un pot vibrant pour la mise en résonance d'un piston dédié ; [0024] - des moyens d'émission des impulsions ultrasonores haute fréquence et des moyens de génération de l'onde acoustique basse fréquence orientés de sorte que les directions de propagation des impulsions ultrasonores et de l'onde acoustique basse fréquence sont perpendiculaires. [0025] L'invention a également pour objet un méthode acoustique de mesure, localisée et sans contact, des micro-endommagements d'un échantillon disposé dans un réservoir, comportant une étape d'émission d'une onde acoustique basse fréquence pour créer une variation périodique basse fréquence de la pression hydrostatique dans l'échantillon, une étape de génération d'impulsions ultrasonores haute fréquence traversant l'échantillon, et une étape de caractérisation, via une unité d'analyse comportant un module de calcul des modulations du temps de vol et de l'amplitude des impulsions ultrasonores lors du passage de l'onde acoustique basse fréquence dans l'échantillon, des effets non-linéaire élastiques et dissipatifs ainsi que de la viscoélasticité de l'échantillon. [0026]Selon des mises en oeuvre particulières: [0027]- l'unité d'analyse réalise un calcul instantané des effets non-linéaires élastiques et dissipatifs ainsi que du paramètre de viscoélasticité de l'échantillon en fonction des phases de traction et de compression de la pression hydrostatique basse fréquence ; [0028]- l'unité d'analyse réalise un calcul, dans le domaine fréquentiel, des non-linéarités élastique et dissipative en fonction de l'amplitude moyenne de la pression hydrostatique basse fréquence, pour une fréquence nulle (ordre 0), une fréquence fondamentale de résonance du piston (ordre 1) et une fréquence double du piston (ordre 2) ; [0029]- la fréquence de l'onde acoustique basse fréquence est comprise en quelques Hz et 100kHz, les impulsions ultrasonores présentent une gamme de fréquence de 20kHz à 100MHz et une cadence de tir environ 10 fois supérieure à la fréquence de l'onde acoustique basse fréquence ; [0030]l'unité d'analyse traite une séquence résultant du moyennage de 10 à 100 tirs de trains d'onde acoustique basse fréquence successifs. [0031] Le fait de générer une variation de la pression hydrostatique à l'intérieur du réservoir à l'aide d'une onde acoustique basse fréquence interagissant de manière synchronisée dans l'échantillon avec une onde � sonde � ultrasonore permet d'obtenir une mesure instantanée des effets viscoélastiques et des effets non linéaires élastique et dissipatif dans l'échantillon. Par ailleurs, l'utilisation d'une onde acoustique pour générer les variations de la pression hydrostatique dans le milieu permet d'étudier en dynamique les effets non-linéaires acoustiques pendant les phases successives de traction et de compression de l'échantillon. [0032] Une hystérésis peut être éventuellement mesurée entre les phases d'augmentation et de diminution de la pression hydrostatique. [0033] L'association d'un piston à un pot vibrant permet de concevoir et réaliser à façon des moyens de génération de l'onde basse fréquence avec une grande latitude pour le choix de la vitesse de déformation. En effet, les paramètres géométriques et structurels du piston déterminent sa fréquence de résonance. [0034] L'utilisation d'une onde acoustique basse fréquence faisant varier les conditions de pression hydrostatique de manière homogène dans les trois dimensions de l'espace permet de � sonder � l'échantillon en faisant varier la direction de propagation de l'onde ultrasonore et étudier ainsi l'anisotropie des non-linéarités viscoélastique et dissipative. [0035] La génération sans contact dans le milieu d'une contrainte basse fréquence en traction et en compression au moyen d'une onde acoustique permet une étude des non-linéarités viscoélastique et dissipative : - dans les milieux fluides et gels non préhensibles dans une machine de test mécanique classique ; - dans les milieux à protéger d'une éventuelle contamination extérieure (milieux biologiques et agro-alimentaires) ; - dans les milieux solides à géométries complexes ; et - pour l'application au vivant (à titre d'exemple, étude du micro-endommagement dans l'os du talon). BREVE DESCRIPTION DES FIGURES [0036] D'autres caractéristiques et avantages de l'invention ressortiront à la lecture qui suit d'un exemple de réalisation détaillé, en référence aux figures annexées qui représentent respectivement : la figure 1, une vue schématique d'un exemple de réalisation d'un dispositif de mesure des paramètres non linéaires élastique et dissipatif et du paramètre viscoélastique d'un échantillon selon l'invention ; la figure 2, une illustration de la pression hydrostatique (onde acoustique basse fréquence) mesurée à l'aide d'un hydrophone situé dans le réservoir dans le même plan que les transducteurs ultrasonores sous le piston basse fréquence, selon l'invention ; la figure 3, une illustration de la modulation du temps de vol des trains d'ondes ultrasonores hautes fréquences, en fonction du temps, se propageant dans un os de talon humain ; 25 la figure 4, une illustration du diagramme instantané de la modulation du temps de vol des impulsions ultrasonores représentant les non-linéarités élastiques en fonction des variations de la contrainte basse fréquence (pression hydrostatique) dans un os de talon ; la figure 5, une illustration du taux de modulation d'amplitude des trains d'ondes ultrasonores hautes fréquences, en fonction du temps, se propageant dans un os de talon humain ; la figure 6, une illustration du diagramme instantané du taux de modulation d'amplitude des impulsions ultrasonores représentant les non-linéarités dissipatives en fonction des variations de la contrainte basse fréquence (pression hydrostatique) dans un os de talon humain ; 15 les figures 7a et 7b, des illustrations de la modulation du temps de vol et du taux de modulation d'amplitude des trains d'ondes ultrasonores hautes fréquences, mesurées dans le domaine fréquentiel, aux ordres 0, 1 et 2, correspondant respectivement à la fréquence nulle, à la 20 fréquence fondamentale de résonance du piston, et à la fréquence double du piston, après propagation dans un os de talon humain. DESCRIPTION DETAILLEE D'UN MODE DE REALISATION [0037] La figure 1 illustre un exemple de réalisation d'un dispositif de mesure 25 instantanée, localisée et sans contact des variations de célérité et d'atténuation d'impulsions ultrasonores hautes fréquences se propageant dans un échantillon (2) soumis à des variations de pression hydrostatique générées par une onde acoustique basse fréquence. [0038] Dans cet exemple, l'échantillon (2) analysé est placé dans un réservoir (4) 30 empli d'un liquide (6) tel que de l'eau. Remplir le réservoir d'eau permet d'assurer une propagation correcte des impulsions ultrasonores. Néanmoins ce liquide peut10 avantageusement être remplacé par un autre fluide liquide pour les matériaux oxydables. [0039] Le dispositif comprend un appareil de synchronisation (8) permettant de générer un signal de synchronisation (8a), de préférence sous forme carrée. Ce signal de synchronisation (8a) permet la synchronisation d'un générateur de signaux (10) et d'un oscilloscope (30). [0040]Le générateur (10) de signal comporte deux voies (10a) et (10b) et permet de délivrer deux signaux synchronisés, de préférence sous forme sinusoïdale. Alternativement, deux générateurs peuvent être utilisés à la place du générateur double voies (10). Dans ce cas, les deux générateurs de signaux sont synchronisés par le signal de synchronisation (8a). [0041]La première voie (10a) de ce générateur (10) est avantageusement reliée à un amplificateur de puissance haute fréquence (12) et permet d'obtenir un signal d'émission électrique haute fréquence (12a). L'amplificateur de puissance haute fréquence (12) est relié à un transducteur ultrasonore d'émission (14) permettant de transformer le signal d'émission électrique haute fréquence (12a), constitué d'impulsions électriques haute fréquence, en un train d'impulsions ultrasonores (14a). [0042]La fréquence du transducteur d'émission (14) sera avantageusement comprise entre 20 kHz et 100 MHz. À titre d'exemple, la fréquence nominale déterminée pour une application dans l'os du talon est égale à 1 MHz. Le transducteur d'émission (14) est également placé dans le réservoir (4) et disposé en vis-à-vis d'un transducteur ultrasonore de réception (16) afin que les impulsions ultrasonores (14a) traversant l'échantillon (2) soient reçues au niveau du transducteur de réception (16). Les impulsions ultrasonores (14a) peuvent être constituées d'une à quelques périodes ultrasonores. [0043] Un signal de réception électrique haute fréquence (18) issu du transducteur ultrasonore de réception (16) est avantageusement transmis à un amplificateur de réception haute fréquence (28), puis numérisé par l'oscilloscope numérique (30). [0044]Les transducteurs d'émission (14) et de réception (16) peuvent être des mono-transducteurs plans ou focalisés, des réseaux multi-éléments annulaires plans ou focalisés, ou des réseaux multi-éléments linéaires plans ou focalisés, mono ou bi-dimentionnel. Le transducteur de réception (16) peut également être un hydrophone. À titre d'exemple, les résultats dans l'os du talon ont été obtenus à l'aide de mono-transducteurs plans. [0045] La deuxième voie (10b) du générateur (10) est, de préférence, reliée à un amplificateur de puissance basse fréquence (20). Cet amplificateur de puissance basse fréquence (20) transmet un signal électrique basse fréquence (20a) amplifié à un pot vibrant (22) apte à mettre en mouvement un piston (24) solidaire de l'axe du pot vibrant (22) et disposé à l'intérieur du réservoir (4). Le pot vibrant (22) permet de générer des vitesses de déformation avantageusement comprises entre 10Hz et 15kHz. À titre d'exemple, des fréquences comprises entre 2 et 5kHz ont été utilisées dans l'application relative à l'os du talon. Pour un fonctionnement idéal, il est nécessaire d'envoyer quelques dizaines de périodes sinusoïdales afin de laisser le piston (24) atteindre sa résonance d'amplitude maximale. [0046] La longueur d'onde basse-fréquence doit être nettement supérieure à la distance entre les transducteurs ultrasonores émetteur (14) et récepteur (16) afin de pouvoir considérer le champ de pression basse fréquence quasi-statique pendant la durée de propagation d'une impulsion ultrasonore (14a) entre les transducteurs ultrasonores (14) et (16). [0047] Le mouvement du piston (24) génère une onde acoustique basse fréquence (24a), de préférence sinusoïdale et présentant avantageusement une direction de propagation orthogonale à la direction de propagation des impulsions ultrasonores (14a). Cette onde acoustique basse fréquence (24a) génère une variation sinusoïdale de la pression hydrostatique dans une région localisée du réservoir (4). En effet, la longueur d'onde basse fréquence (24a) est très grande vis-à-vis d'au moins une des dimensions de l'échantillon (2) sondé. [0048] Alternativement à une interaction orthogonale entre les impulsions ultrasonores (14a) et l'onde acoustique basse fréquence (24a), il est possible de disposer la paire de transducteurs ultrasonores (14) et (16) avec un angle quelconque vis-à-vis de la direction de propagation de l'onde basse fréquence (24a). [0049] La configuration d'une interaction orthogonale entre les directions de propagation des impulsions ultrasonores (14a) et basse fréquence (24a) autorise une étude de l'anisotropie des paramètres non linéaires de l'échantillon (4) dans un plan parallèle au piston (24). Par ailleurs, si on considère un angle quelconque entre le plan des transducteurs haute fréquence (14,16) et la direction de propagation de l'onde acoustique basse fréquence (24a), on accède à une étude de l'anisotropie des paramètres non linéaires dans les trois dimensions de l'espace. [0050] Afin de maximiser les variations d'amplitude de l'onde acoustique basse fréquence (24a), le piston (24) fonctionne de préférence sur son premier mode de résonance. Les choix de la densité et de la rigidité du matériau constitutif du piston (24), ainsi que son diamètre et son épaisseur, déterminent sa fréquence de résonance et permettent d'obtenir des variations de pression hydrostatique de l'ordre de 70 kPa. [0051]À titre d'exemple, un piston en verre de diamètre 140mm et d'épaisseur 6mm permet de générer une vibration à la fréquence 2800Hz et permet d'atteindre des variations d'amplitude de pression hydrostatique de 70kPa à 15mm de la surface du piston. [0052] Les moyens d'émission de l'onde acoustique basse fréquence (24a) pourraient alternativement être composés d'un haut-parleur sous-marin (projecteur acoustique) ou sonar afin d'atteindre des vitesses de déformation supérieures à 15kHz tout en maintenant des niveaux de pression élevés. [0053] Les variations de la pression hydrostatique basse fréquence sont mesurées à l'aide d'un hydrophone (26) placé parallèlement au piston dans le plan des transducteurs haute fréquence (14) et (16). Ainsi, l'onde basse fréquence générée dans le réservoir (4) est mesurée par cet hydrophone (26). Le signal électrique basse fréquence (26a) issu de l'hydrophone (26) est avantageusement transmis à un amplificateur (34) puis communiqué à l'oscilloscope (30) pour numérisation. La forme typique de ce signal électrique basse fréquence (26a) amplifié est illustrée, à titre d'exemple, dans la figure 2. [0054] Les courts trains d'onde haute fréquence (14a) traversent l'échantillon (2) avec une cadence de tirs permettant un échantillonnage correct des variations de la pression hydrostatique dans le réservoir (4). Typiquement pour une onde acoustique basse fréquence (24a) de fréquence 2800 Hz, une cadence des impulsions ultrasonores (14a) d'environ 25 kHz est utilisée, soit environ 10 tirs ultrasonores par période de l'onde acoustique basse fréquence (24a). [0055] La longueur temporelle de la fenêtre d'acquisition des signaux électriques de réception haute fréquence (18) et du signal électrique basse fréquence (26a) est judicieusement déterminée afin d'enregistrer une partie du signal non soumis aux variations de pression basse fréquence et la totalité de la montée en résonance du piston basse fréquence. À titre d'exemple, si on utilise un piston résonant à 2800Hz, une fenêtre de 70ms est nécessaire pour couvrir l'ensemble de l'expérience (figure 2). [0056] Afin d'optimiser le rapport signal sur bruit du signal électrique de réception haute fréquence (18), un moyennage sur N expériences est effectué en envoyant N fois un train d'ondes acoustiques basse fréquence (24a). [0057] Une fréquence d'échantillonnage du signal électrique de réception haute fréquence (18) nettement supérieure (de l'ordre de 10 à 50 fois supérieure) à la fréquence des impulsions ultrasonores (14a) est choisie. À titre d'exemple' on obtient dans l'eau un niveau de bruit de mesure de la modulation du temps de vol de l'ordre de 10-10 secondes pour une fréquence d'échantillonnage de 25MHz, un moyennage N=30 et des impulsions ultrasonores de fréquence 1 MHz. [0058]Le signal électrique de réception haute fréquence (18) et le signal électrique basse fréquence (26a) numérisés par l'oscilloscope (30) sont ensuite transmis à l'unité d'analyse (32) pour stockage et traitement des données. [0059] L'unité d'analyse (32) permet de traiter les informations contenues dans le signal électrique de réception haute fréquence (18) et le signal électrique basse fréquence (26a) préalablement numérisés et synchronisés afin de calculer les non-linéarités élastique et dissipative ainsi que le coefficient de viscoélasticité de l'échantillon (2). [0060] Les non-linéarités élastique et dissipative et la viscoélasticité sont calculées à partir des variations du temps de vol et de l'amplitude du signal électrique de réception haute fréquence (18) en fonction des variations de la pression basse fréquence. [0061] Les variations du temps de vol des courts trains d'ondes ultrasonores (14a), générés pendant les phases de compression et de détente de la pression hydrostatique basse fréquence, sont liées aux variations de la vitesse de propagation ultrasonore dans le milieu (équation 1) : aTOF L ac aPBF C' 9PBF , où Eq. 1 - L est la distance entre les transducteurs ultrasonores émetteur et récepteur (14) et (16) ; - TOF (Time Of Flight) est le temps de propagation des impulsions hautes fréquences (14a) entre les transducteurs ultrasonores (14) et (16) ; - c est la vitesse de propagation de l'onde ultrasonore dans le milieu considéré soumis aux variations de pression basse fréquence ; - Co est la vitesse de propagation de l'onde ultrasonore dans le milieu considéré en l'absence de variations de pression basse fréquence ; [0062] La vitesse de propagation ultrasonore étant reliée au module élastique K et à la densité du milieu considéré po, on peut écrire de la même façon les variations de la célérité c ou du temps de vol (TOF) en fonction des variations de K (équation 2): cTOF L aK ÔPBF 2Co po ÔPBF où Eq.2 - K est le module élastique ; po est la densité du milieu en l'absence de variations de la pression basse 10 fréquence. [0063] L'unité d'analyse (32) comprend un module de calcul des non-linéarités élastiques et les non-linéarités dissipatives de l'échantillon (2) en fonction du temps, telles qu'illustrées sur les figures 3 et 5, respectivement. L'unité d'analyse (32) comporte également un module de calcul du paramètre de viscosité. 15 [0064] Le module de calcul des non-linéarités élastiques applique un algorithme d'intercorrélation entre la première impulsion haute fréquence issue du signal électrique de réception haute fréquence (18) (correspondant à une impulsion ultrasonore (14a) n'ayant pas été soumise aux variations de la pression basse fréquence) et chacune des impulsions haute fréquence suivantes de ce signal 20 électrique de réception (18). L'intercorrélation conduit à l'estimation du délai entre la première impulsion haute fréquence de référence et la Hème impulsion haute fréquence du signal électrique de réception (18). Afin d'optimiser l'estimation de ce délai, une interpolation parabolique du pic d'intercorrélation est implémentée dans le but de s'affranchir de l'échantillonnage numérique. Le résultat in fine conduit à 25 la modulation du temps de vol (Time Of Flight Modulation : TOFM) en fonction du temps (figure 3). [0065] Alternativement à un algorithme d'intercorrélation, un algorithme de démodulation de la phase des signaux issus de l'onde sonde ultrasonore (14a) peut être utilisé dans le cas où cette dernière est envoyée en continu durant le temps d'un tir de train d'ondes basse fréquence (24a). [0066] Pour le module de calcul des non-linéarités dissipatives, trois approches sont possibles : 1. La mesure du taux de variation de l'amplitude pic-pic des impulsions du signal électrique de réception haute fréquence (18) dans le domaine temporel ; 2. La mesure du taux de variation de l'amplitude maximale des impulsions du signal électrique de réception haute fréquence (18) dans le domaine fréquentiel (figure 5); 3. La mesure du taux de variation de la pente d'atténuation fréquentielle. [0067]À l'instar des non-linéarités élastiques, la variation est calculée dans chacun de ces trois cas à partir d'une impulsion du signal électrique de réception haute fréquence (18), signal de référence n'ayant pas subi de variations de pression basse fréquence. Dans le cas de l'approche 3), l'atténuation fréquentielle, appelée couramment BUA (Broadband Ultrasonic Attenuation) dans les applications osseuses, correspond à la pente du rapport des modules des spectres du signal de référence (première impulsion du signal électrique de réception haute fréquence (18)) et de la nème impulsions du signal électrique de réception haute fréquence (18) sur la bande passante considérée. Le résultat in fine conduit au taux de variation en amplitude ou en atténuation en fonction du temps (figure 5). [0068] À partir, d'une part, des fonctions temporelles des non-linéarités élastique et dissipative exposées précédemment et, d'autre part, de la mesure synchronisée des variations de la pression basse fréquence en fonction du temps, il est possible d'analyser et présenter les résultats selon deux modalités : 1) Modalité instantanée : diagrammes instantanés des non-linéarités élastique et dissipative en fonction de l'amplitude de la pression basse fréquence ; 2) Modalité fréquentielle : analyses fréquentielles des non-linéarités élastique et dissipative aux ordres 0, 1 et 2 (respectivement fréquence nulle, fréquence fondamentale de résonance du piston et fréquence double du piston) en fonction de l'amplitude moyenne de la pression basse fréquence. [0069] Dans le cas de la modalité instantanée, chaque point de la mesure de l'amplitude de la pression basse fréquence est synchronisé avec un point de calcul de la modulation du temps de vol (non-linéarité élastique) ou du taux de la modulation d'amplitude (non-linéarité dissipative). Ces mesures synchronisées conduisent naturellement à une représentation de la modulation du temps de vol ou de l'amplitude des impulsions ultrasonores (14a) en fonction de l'amplitude de la pression basse fréquence, et au suivi instantané des non-linéarités durant les phases de compression et de détente dans l'échantillon (2). La fenêtre d'analyse peut être choisie durant la phase de montée en résonance du piston, sur le plateau une fois la résonance établie et/ou durant la phase de relaxation du piston (après l'arrêt du signal électrique (10b). [0070] Si l'échantillon considéré ne présente pas de non-linéarité de type hystérétique (milieux endommagés ou granulaires), et que cependant ce diagramme présente une hystérésis, celui-ci peut conduire à la mesure d'un paramètre de viscoélasticité. Ce dernier est obtenu, par le biais d'une unité de calcul du paramètre de viscoélasticité, à partir de la mesure du déphasage D entre la pression basse fréquence et la modulation de temps de vol des signaux électriques de réception haute fréquence (18). Le coefficient de viscoélasticité (ou facteur de perte) couramment utilisé vaut alors tan((D). [0071] À titre indicatif, les figures 4 et 6 illustrent respectivement un exemple de modulation du temps de vol et un exemple du taux de modulation de l'amplitude des signaux électriques de réception haute fréquence (18) en fonction de l'amplitude de la pression basse fréquence dans le cadre de la mesure d'un os du talon. [0072] À l'instar d'un essai mécanique sur une machine de test quasi-statique, une telle analyse permet de visualiser distinctement et au cours d'un même essai les comportements en traction, pour une variation de pression hydrostatique négative, et en compression, pour une variation de pression hydrostatique positive. [0073] Dans le cas de la modalité fréquentielle, les modulation de temps de vol et d'amplitude des signaux électriques de réception haute fréquence (18) ainsi que l'amplitude de la pression basse fréquence en fonction du temps sont découpées en fenêtres glissantes de tailles égales à 4 ou 5 périodes basse fréquence durant les phases de montée en résonance du piston et de descente. Les fenêtres des deux modulations non-linéaires (élastique et dissipative) sont analysées dans le domaine fréquentiel, où les amplitudes du module de ces spectres sont relevées pour la fréquence nulle (ordre 0), la fréquence fondamentale de résonance du piston (ordre 1) et la fréquence double du piston (ordre 2). [0074] En parallèle, les modules des spectres des fenêtres correspondant à l'amplitude de la pression basse fréquence conduisent au calcul de l'amplitude moyenne de cette pression basse fréquence. Le choix de 4 à 5 périodes basse fréquence résulte du compromis entre une longueur suffisante pour calculer une transformée de Fourier pertinente et une longueur minimale pour pouvoir considérer constante l'amplitude moyenne de la pression basse fréquence. [0075] À titre indicatif, les figures 7a et 7b illustrent respectivement un exemple de modulation du temps de vol et un exemple du taux de modulation de l'amplitude des signaux électriques codés et amplifiés (18), obtenus à partir de la mesure dans le domaine fréquentiel aux ordres 0, 1 et 2, en fonction de l'amplitude moyenne de la pression basse fréquence dans le cadre de la mesure d'un os du talon. [0076] Contrairement à la modalité instantanée, l'information liée aux phases de traction et de compression n'est plus accessible, l'amplitude de la pression basse fréquence étant moyennée sur quelques périodes dans le domaine de Fourier. En revanche, l'analyse fréquentielle conduit à une décomposition des non-linéarités selon différents ordres fréquentiels : constante, fondamental et harmoniques. Par ailleurs, elle permet d'étudier les comportements éventuellement différents lors des phases d'augmentation et de diminution de l'amplitude basse fréquence, lors de la résonance du piston. [0077] Une relation phénoménologique (K.R. McCall et R.A. Guyer; Equation of state and wave propagation in hysteretic nonlinear elastic materials; Journal of Geophysical Research, Vol. 99(B12), 1984; pp.23 887-23 897.) basée sur un développement en série de Taylor jusqu'à l'ordre 2 du module élastique K en fonction de la pression basse fréquence, modifiée pour prendre en compte les effets non-linéaires non classiques ou hystérétiques est utilisée pour extraire ou identifier les paramètres non-linéaires classiques 13, b (ordres 1 et 2) et non classique a selon l'équation : K = Ko + .PBF - b. PBF i + a A(PBF) - PBF.sign PBF , Eq.3 où : Ko - Ko est le module élastique linéaire ; - PBF est la dérivée temporelle de la pression basse fréquence. [0078] Cette identification peut être effectuée à partir d'un ajustement sur le diagramme instantané reliant la modulation de temps de vol (ou élastique) aux variations de la pression basse fréquence. [0079] L'invention n'est pas limitée aux exemples de réalisation décrits et représentés. Il est également possible de prévoir des modes de réalisation pour lesquels le réservoir (4) contiendrait un élément en phase gazeuse et non pas liquide. Par ailleurs, les différents modes de réalisation de l'unité d'analyse (32) ainsi que les résultats obtenus dans ces modes peuvent être combinés de différentes façons sans sortir du cadre de l'invention. REVENDICATIONS1. Dispositif acoustique de mesure, localisée et sans contact, des micro- endommagements d'un échantillon (2) disposé dans un réservoir (4), caractérisé en ce qu'il comporte des moyens d'émission (10, 20, 22, 24) d'une onde acoustique basse fréquence (24a) pour créer une variation périodique basse fréquence de la pression hydrostatique dans l'échantillon (2), des moyens de génération (10, 12, 14) d'impulsions ultrasonores haute fréquence (14a) et une unité d'analyse (32) comportant un module de calcul des modulations du temps de vol et de l'amplitude des impulsions ultrasonores (14a) lors du passage de l'onde acoustique basse fréquence (24a) dans l'échantillon (2), de sorte à caractériser les effets non-linéaire élastiques et dissipatifs ainsi que la viscoélasticité de l'échantillon (2). 2. Dispositif acoustique de mesure selon la revendication 1, comportant un module d'imagerie multi-paramétrique des effets viscoélastiques (tan()) et des effets non-linéaires acoustiques élastique (a, 13 et 8) et dissipatif. 3. Dispositif acoustique de mesure selon la revendication 1 ou 2, dans lequel les moyens de génération (24) de l'onde acoustique basse fréquence (24a) comportent un pot vibrant (20) pour la mise en résonance d'un piston (24) dédié. 4. Dispositif acoustique de mesure selon l'une quelconque des revendications 1 à 3, dans lequel les moyens d'émission (10, 12, 14) des impulsions ultrasonores haute fréquence (14a) et les moyens de génération (10, 20, 22, 24) de l'onde acoustique basse fréquence (24a) sont orientés de sorte que les directions de propagation des impulsions ultrasonores (14a) et de l'onde :30 acoustique basse fréquence (24a) sont perpendiculaires. 5. Méthode acoustique de mesure, localisée et sans contact, des micro-endommagements d'un échantillon (2) disposé dans un réservoir (4), caractérisé en ce qu'il comporte une étape d'émission d'une onde acoustique basse fréquence (24a) pour créer une variation périodique basse fréquence de la pression hydrostatique dans l'échantillon (2), une étape de génération d'impulsions ultrasonores haute fréquence (14a) traversant l'échantillon (2), et une étape de caractérisation, via une unité d'analyse (32) comportant un module de calcul des modulations du temps de vol et de l'amplitude des impulsions ultrasonores (14a) lors du passage de l'onde acoustique basse fréquence (24a) dans l'échantillon (2), des effets non-linéaire élastiques et dissipatifs ainsi que de la viscoélasticité de l'échantillon (2). 6. Méthode acoustique de mesure selon la revendication 5, dans laquelle l'unité d'analyse (32) réalise un calcul instantané des effets non-linéaires élastiques et dissipatifs ainsi que du paramètre de viscoélasticité de l'échantillon (2) en fonction des phases de traction et de compression de la pression hydrostatique basse fréquence. 7. Méthode acoustique de mesure selon la revendication 5, dans laquelle l'unité d'analyse (32) réalise un calcul, dans le domaine fréquentiel, des non-linéarités élastique et dissipative en fonction de l'amplitude moyenne de la pression hydrostatique basse fréquence, pour une fréquence nulle (ordre 0), une fréquence fondamentale de résonance du piston (ordre 1) et une fréquence double du piston (ordre 2) . 8. Méthode acoustique de mesure selon l'une quelconque des revendication 5 à 7, dans laquelle la fréquence de l'onde acoustique basse fréquence (24a) est comprise en quelques Hz et 100kHz, les impulsions ultrasonores (14a) présentent une gamme de fréquence de 20kHz à 100MHz et une cadence de tir environ 10 fois supérieure à la fréquence de l'onde acoustique basse fréquence (24a). 9. Méthode acoustique de mesure selon l'une quelconque des revendications 5 à 8, dans laquelle l'unité d'analyse (32) traite une séquence résultant du moyennage de 10 à 100 tirs de trains d'onde acoustique basse fréquence (24a) successifs.