Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
METHOD FOR DETERMINING A PHYSICAL PARAMETER, IMAGING METHOD, AND DEVICE FOR IMPLEMENTING SAID METHOD
Document Type and Number:
WIPO Patent Application WO/2011/018573
Kind Code:
A1
Abstract:
The invention relates to a method for determining a physical parameter representative of a point P of a plate by means of a device comprising a first receiver for measuring a wave propagating in the plate, and a calculation unit. The method includes the following steps: measuring, by means of the first receiver, a first signal s1 (t) representative of a wave propagating in the plate; defining a closed contour C on the plate surrounding said point P, the contour C being the location on the plate on which either a wave is generated by a first emitter, or the first signal s1 (t) is measured; and determining the physical parameter at point P on the plate by identifying, thanks to at least said first signal, a shape function fshape(f), in the following equation: formula (I), where Wcontour and Gplate are functions of Green, and the shape function fshape(f) is dependent on the frequency / of the wave and on the physical parameter.

Inventors:
ING, Ros Kiri (11 rue de Chaligny, Paris, F-75012, FR)
ETAIX, Nicolas (63 avenue Paul Vaillant Couturier, Gentilly, F-94250, FR)
LEBLANC, Alexandre (38 rue Anatole France, Saint-Nicolas, F-62223, FR)
FINK, Mathias (16 rue Edouard Laferrière, Meudon, F-92190, FR)
Application Number:
FR2010/051618
Publication Date:
February 17, 2011
Filing Date:
July 29, 2010
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE - CNRS - (3 rue Michel Ange, Paris, Paris, F-75016, FR)
ING, Ros Kiri (11 rue de Chaligny, Paris, F-75012, FR)
ETAIX, Nicolas (63 avenue Paul Vaillant Couturier, Gentilly, F-94250, FR)
LEBLANC, Alexandre (38 rue Anatole France, Saint-Nicolas, F-62223, FR)
FINK, Mathias (16 rue Edouard Laferrière, Meudon, F-92190, FR)
International Classes:
G01N29/44; G01B17/02
Attorney, Agent or Firm:
BURBAUD, Eric et al. (Cabinet Plasseraud, 52 rue de la Victoire, Paris Cedex 09, F-75440, FR)
Download PDF:
Claims:
REVENDICATIONS

1. Procédé pour déterminer un paramètre physique représentatif d'un point P d'une plaque, dans lequel le paramètre physique est choisi parmi l'épaisseur de la plaque h, la vitesse de propagation d'une onde dans la plaque VP et le produit Vph de l'épaisseur par la vitesse de propagation d'une onde dans la plaque, et

ledit procédé est mis en œuvre par un dispositif comprenant :

- au moins un premier récepteur (Rl) adapté pour réaliser une mesure d'une onde se propageant dans la plaque, et

- une unité de calcul (CALC) reliée audit premier récepteur (Rl),

ledit procédé étant caractérisé en ce que ledit procédé comprend les étapes suivantes :

- on mesure par ledit premier récepteur (Rl) un premier signal srft) représentatif d'une onde se propageant dans la plaque,

- on définit sur la plaque un contour fermé C entourant ledit point P, le contour C étant le lieu de la plaque sur lequel soit on génère par un premier émetteur

(El) une onde se propageant dans la plaque, soit on mesure par ledit premier récepteur (Rl) ledit premier signal srft) représentatif d'une onde se propageant dans la plaque, et

- on détermine le paramètre physique au point P de la plaque en identifiant grâce au moins audit premier signal srft) une fonction de forme fshaPe(f), dans la relation suivante :

W contour (?) = fshape (f Y* plate (/ ~ Ks )

Wcontour est une fonction de Green représentant l'onde le long du contour C,

Gpiate est une fonction de Green représentant l'onde en un point de vecteur position r de la plaque n'appartenant pas au contour C, par rapport à un point S de vecteur position rs de la plaque représentant une source de l'onde, et

ladite fonction de forme fshapJJ) est dépendante au moins de la fréquence / de l'onde et du paramètre physique, et adaptée à la forme du contour C.

2. Procédé selon la revendication 1, dans lequel la fonction de forme est une fonction de Bessel de première espèce Jo(Z) comprenant des zéros Zn, n étant un entier positif ou nul, ladite fonction de Bessel est fonction d'un paramètre d'échelle a multiplié par la racine carrée de la fréquence / de l'onde, de telle sorte que :

3. Procédé selon la revendication 2, dans lequel le premier récepteur (Rl) est adapté pour mesurer une onde sur le contour C, et ledit procédé comprenant l'étape suivante :

- si le premier signal si(t) a une amplitude inférieure à un seuil prédéterminé pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Zn de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

"

4. Procédé selon la revendication 2, dans lequel le dispositif comprend en outre un premier émetteur (El), et l'un des premier émetteur (El) et du premier récepteur (Rl) étant adapté soit pour générer ou mesurer une onde sur le contour C, l'autre étant adapté soit pour générer ou mesurer une onde en un point de la plaque n'appartenant pas au contour C, et ledit procédé comprenant les étapes suivantes : - on génère dans la plaque une première onde par un premier signal d'émission erft) alimentant le premier émetteur (El ) ,

- on mesure par le premier récepteur (Rl) un premier signal srft) représentatif de ladite première onde émise, et

- si le premier signal ei(t) a une amplitude inférieure à un seuil prédéterminé pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Zn de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

a=——

5. Procédé selon la revendication 2, dans lequel le dispositif comprend en outre un deuxième récepteur (R2), le premier récepteur (Rl) étant adapté pour mesurer une onde sur le contour C et le deuxième récepteur (R2) étant adapté pour mesurer une onde en un point de la plaque n'appartenant pas au contour C, et ledit procédé comprenant les étapes suivantes :

- on mesure simultanément par le premier récepteur (Rl) un premier signal sj(t) représentatif d'une première onde, et par le deuxième récepteur (R2) un deuxième signal S2(t) représentatif de ladite même première onde.

6. Procédé selon la revendication 2, dans lequel le dispositif comprend en outre :

- un deuxième récepteur (R2), et

- un premier émetteur (El),

l'un des premier récepteur (Rl), deuxième récepteur (R2) et premier émetteur (El) étant adapté pour mesurer ou générer une onde sur le contour C et ledit procédé comprenant les étapes suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur (El ) ,

- on mesure simultanément par le premier récepteur

(Rl) un premier signal sj(t) représentatif de ladite première onde émise, et par le deuxième récepteur (R2) un deuxième signal S2(t) représentatif de ladite même première onde émise .

7. Procédé selon la revendication 2, dans lequel le dispositif comprend en outre :

- un premier émetteur (El), et

- un deuxième émetteur (E2),

l'un des premier récepteur (Rl), premier émetteur (El) et deuxième émetteur (E2) étant adapté pour mesurer ou générer une onde sur le contour C et ledit procédé comprenant les étapes suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission erft) alimentant le premier émetteur (El ) ,

- on mesure par le premier récepteur (Rl) un premier signal srft) représentatif de ladite première onde émise,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission β2(t) alimentant le deuxième émetteur (E2) ,

- on mesure par le premier récepteur (Rl) un deuxième signal S2(t) représentatif de ladite deuxième onde émise . 8. Procédé selon la revendication 7, dans lequel le deuxième signal d'émission e2(X) est déphasé de π/2 par rapport au premier signal d'émission erft), et ledit procédé comprenant les étapes suivantes :

- on calcule un signal sommé s(t), qui est la somme du premier signal srft) et d'un deuxième signal S2(t), et

- si le premier signal srft) est en phase avec le signal sommé s(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Zn de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par

α=-^=

9. Procédé selon la revendication 2, dans lequel le dispositif comprend en outre :

- un premier émetteur (El), et

- un deuxième émetteur (E2),

l'un des premier récepteur (Rl), premier émetteur (El) et deuxième émetteur (E2) étant adapté pour mesurer ou générer une onde sur le contour C et ledit procédé comprenant les étapes suivantes :

- on génère simultanément dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur (El), et une deuxième onde par un deuxième signal d'émission e2(X) alimentant le deuxième émetteur (E2), et

- on mesure par le premier récepteur (Rl) un premier signal srft) représentatif de la superposition sur le lieu du premier récepteur (Rl) desdites première et deuxième ondes émises. 10. Procédé selon la revendication 9, dans lequel le deuxième signal d'émission e2(X) est déphasé de π/2 par rapport au premier signal d'émission erft), et ledit procédé comprenant les étapes suivantes :

- si le premier signal srft) est en phase avec le premier signal d'émission eι(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Zn de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par : a =_ ^n

11. Procédé selon l'une des revendications 5 à 7, comprenant les étapes suivantes :

- on calcule un signal sommé s(t), qui est la somme du premier signal srft) et d'un deuxième signal déphasé S2*(t), le deuxième signal déphasé S2* (t) étant égal au deuxième signal S2(t) déphasé de π/2, et

- si le premier signal srft) est en phase avec le signal sommé s(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Zn de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

a=——

12. Procédé selon l'une des revendications 5 à 7, comprenant les étapes suivantes :

- on calcule une première transformée de Fourier Si(f) du premier signal srft) et une deuxième transformée de Fourier ^(/"J du deuxième signal S2(t),

- on calcule une fonction de test ftest(f), qui compare le signe de la partie réelle de la première transformée de Fourier Si(f) au signe de la partie réelle de la deuxième transformée de Fourier ^f/J, et qui affecte une première valeur Vi si les signes sont identique et une deuxième valeur V2 si les signes sont différents :

\ si sign s j{f))) = sign s 2{f))) alors ftest(f) = Vj

sinon ftest(f) = v2

- on recherche des fréquences particulières /„ auxquelles la fonction de test ftest(f) change de valeur, en passant soit de la première valeur Vi à la deuxième valeur V2, soit réciproquement de la deuxième valeur V2 à la première valeur Vi, et

- on calcule le paramètre d'échelle a par : a =_ ^n

13. Procédé selon l'une des revendications 5 à 7, comprenant les étapes suivantes :

- on calcule une première transformée de Fourier Si(f) du premier signal srft) et une deuxième transformée de Fourier ^(/"J du deuxième signal S2(t),

- on calcule une différence de phase Δφ entre la première transformée de Fourier et la deuxième transformée de Fourier, par :

Δφ=φ{S2{f)-S1{f))

- on recherche des fréquences particulières /„ de la différence de phase Δφ, pour lesquelles ladite différence de phase présente un saut entre 0 et π ou entre π et 0, et qui sont proportionnelles au carré des zéros Zn de la fonction de Bessel de première espèce Jo, et

- on calcule le paramètre d'échelle a par

7

a _= ^n

Tn 14. Procédé selon l'une des revendications 5 à 7, comprenant les étapes suivantes :

- on calcule une première transformée de Fourier

Si(f) du premier signal srft) et une deuxième transformée de

Fourier ^(/"J du deuxième signal S2(t),

- on détermine le paramètre d'échelle a de telle sorte que le module de la fonction de forme suivante : où b est un autre paramètre d'échelle, et

1.1 est la fonction module,

s'approche au mieux de :

\s2(fys!(f)\

pour un ensemble de fréquences de test /„.

15. Procédé selon l'une des revendications 5 7, comprenant les étapes suivantes :

- on calcule une première transformée de Fourier Si(f) du premier signal srft) et une deuxième transformée de Fourier S2(T) du deuxième signal S2(t) .

- on détermine le paramètre d'échelle a de telle sorte de la fonction de forme suivante :

où b est un autre paramètre d'échelle, et

φ(.) est la fonction de phase,

s'approche au mieux de :

C(S2(T)ZS1(T))

pour un ensemble de fréquences de test /„.

16. Procédé selon l'une des revendications 2 à 15, dans lequel on détermine le paramètre physique de produit VPh, produit de l'épaisseur par la vitesse de propagation d'une onde dans la plaque, par la formule suivante :

a est le paramètre d'échelle de la fonction de

Bessel, précédemment déterminé, et

R est la longueur d'un segment entre le point P et un point du contour C dans la direction de l'onde. 17. Procédé selon l'une des revendications 2 à 15, dans lequel on détermine le paramètre physique d'épaisseur h de la plaque par la formule suivante :

a est le paramètre d'échelle de la fonction de

Bessel, précédemment déterminé,

R est la longueur d'un segment entre le point P et un point du contour C dans la direction de l'onde, et Vp est la vitesse de propagation connue d'une onde dans le matériau de la plaque.

18. Procédé selon l'une des revendications 2 à 15, dans lequel on détermine le paramètre physique de vitesse de propagation Vp d'une onde dans la plaque par la formule suivante :

a est le paramètre d'échelle de la fonction de

Bessel, précédemment déterminé,

R est la longueur d'un segment entre le point P et un point du contour C dans la direction de l'onde, et

h est l'épaisseur connue de la plaque.

19. Procédé selon l'une des revendications 1 à 18, dans lequel le contour C est sensiblement un cercle de rayon R centré sur le point P.

20. Procédé selon l'une des revendications 1 à 18, dans lequel le contour C est sensiblement une ellipse centrée sur le point P. 21. Procédé selon la revendication 20, dans lequel on détermine au préalable la forme du contour C grâce à :

- un dispositif de test comprenant :

- un premier émetteur (El) adapté pour générer une onde au point P,

- au moins un deuxième émetteur (E2) adapté pour générer une onde sur un contour de test Cn ayant la forme d'une ellipse prédéterminée, n étant un indice entier positif,

- des premier et deuxième récepteurs (Rl, R2 ) adaptés pour mesurer une onde en des points n'appartenant pas au contour de test Cn, et grâce à :

un procédé de test comprenant les étapes de test suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur (El ) ,

- on mesure par le premier récepteur (Rl) un premier signal Su(t) représentatif de ladite première onde émise, et on calcule une première transformée de Fourier Suif) de ce premier signal,

- on mesure par le deuxième récepteur (R2) un deuxième signal Sn{t) représentatif de ladite première onde émise, et on calcule une deuxième transformée de Fourier Snif) de ce deuxième signal,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission β2(t) alimentant le deuxième émetteur (E2 ) ,

- on mesure par le premier récepteur (Rl) un troisième signal S2i(t) représentatif de ladite deuxième onde émise, et on calcule une troisième transformée de Fourier S2i(f) de ce troisième signal,

- on mesure par le deuxième récepteur (R2) un quatrième signal S22(t) représentatif de ladite deuxième onde émise, et on calcule une quatrième transformée de Fourier S22(f) de ce quatrième signal,

- on calcule la fonction de différence de phase suivante :

Δφ(f)=φ(S11(f).S12(f)*)-<p(S21(f).S22(f)*)

* désigne la fonction conjuguée, et

φ(.) est la fonction de phase, et

- la forme de l'ellipse du contour de test Cn correspond à un contour optimal, telle que la première onde se propage et se superpose spatialement sur la plaque sensiblement à la deuxième onde, lorsque la fonction de différence de phase Λφ(f) est minimale pour un ensemble de contours de test. Cn auxquels on applique les étapes de test précédentes .

22. Procédé selon l'une des revendications 1 à 18, dans lequel le contour C est sensiblement un rectangle centré sur le point P.

23. Procédé selon la revendication 22, dans lequel le contour C comporte huit points de contour Cl à C8, et dans lequel lesdits points de contour et le point P forment un maillage rectangle régulier.

24. Procédé selon la revendication 22 et l'une des revendications 16 à 18, dans lequel le segment de longueur R est un segment de longueur moyenne calculé par :

R = dl+d2

2

di est la moitié de la longueur de la plus grande médiane dudit rectangle, et

d,2 est la longueur de la diagonale dudit rectangle.

25. Procédé d'imagerie, dans lequel on construit une image d'une plaque, ladite image comportant une pluralité de pixels, et chaque pixel correspondant à un point de la plaque et représentant un paramètre physique de la plaque audit point de la plaque, et ledit paramètre physique dudit point étant déterminé par le procédé selon l'une des revendications précédentes. 26. Dispositif pour mettre en œuvre le procédé pour déterminer un paramètre physique représentatif d'un point P d'une plaque, dans lequel le paramètre physique est choisi parmi l'épaisseur de la plaque h, la vitesse de propagation d'une onde dans la plaque VP et le produit Vph de l'épaisseur par la vitesse de propagation d'une onde dans la plaque, et ledit dispositif comprenant :

- au moins un premier récepteur (Rl) adapté pour réaliser une mesure d'un premier signal srft) représentatif d'une onde se propageant dans la plaque,

- un contour fermé C défini en entourant ledit point P, le contour C étant le lieu de la plaque sur lequel soit on génère par un premier émetteur (El) une onde se propageant dans la plaque, soit on mesure par ledit premier récepteur (Rl) ledit premier signal srft) représentatif d'une onde se propageant dans la plaque, et

- une unité de calcul (CALC) reliée audit premier récepteur (Rl ) ,

ladite unité de calcul (CALC) étant adaptée pour déterminer le paramètre physique au point P de la plaque en identifiant grâce au moins audit premier signal srft) une fonction de forme fshaPe(f), dans la relation suivante :

^contour(?) = fshape(fY*plate(/ ~ Ks )

Wcontour est une fonction de Green représentant l'onde le long du contour C,

Gpiate est une fonction de Green représentant l'onde en un point de vecteur position r de la plaque n'appartenant pas au contour C, par rapport à un point S de vecteur position fs de la plaque représentant une source de l'onde, et

ladite fonction de forme fshaPe(f) est dépendante au moins de la fréquence / de l'onde et du paramètre physique, et adaptée à la forme du contour C. 27. Dispositif selon la revendication 26, dans lequel le premier récepteur (Rl) est un vibromètre laser à balayage .

28. Dispositif selon la revendication 27, comprenant en outre un deuxième récepteur (R2), et dans lequel le deuxième récepteur (R2) est réalisé par ledit vibromètre laser à balayage.

Description:
Procédé pour déterminer un paramètre physique, procédé d' imagerie et dispositif pour mettre en œuvre ledit

procédé . La présente invention est relative à un procédé pour déterminer un paramètre physique représentatif d'un point P d'une plaque, un dispositif pour mettre en œuvre ledit procédé et un transducteur pour émettre et/ou recevoir une onde sur ladite plaque.

Plus particulièrement, l'invention concerne un procédé pour déterminer un paramètre physique représentatif d'un point P d'une plaque, dans lequel le paramètre physique est choisi parmi l'épaisseur de la plaque h, la vitesse de propagation d'une onde dans la plaque V P et le produit Vph de l'épaisseur par la vitesse de propagation d'une onde dans la plaque, et ledit procédé est mis en œuvre par un dispositif comprenant :

- au moins un premier récepteur adapté pour réaliser une mesure d'une onde se propageant dans la plaque, et

- une unité de calcul reliée audit premier récepteur .

On connaît des procédés de ce type. Notamment, selon une méthode connue, on génère ou émet une onde de vibration en un premier point de la plaque, on mesure l'onde reçue en un deuxième point de la plaque, et on estime le temps de vol de l'onde entre les premier et deuxième points. La distance séparant les premier et deuxième points permet alors d'estimer la vitesse de propagation de l'onde dans la plaque.

Pour avoir une bonne précision sur cette estimation, la précision temporelle est importante et cette méthode nécessite une onde haute fréquence sensible à l'état de surface de la plaque.

De plus, pour déterminer l'épaisseur de la plaque cette méthode nécessite un modèle de la plaque, qui intègre les couplages et conditions limites, ce qui complexifie la méthode et surtout apporte une incertitude liée aux valeurs des paramètres du modèle utilisé.

La présente invention a pour but de proposer une alternative aux méthodes connues, permettant notamment, d'estimer précisément des paramètres physiques d'une plaque .

A cet effet, le procédé de l'invention est caractérisée en ce qu'il comprend les étapes suivantes :

- on mesure par ledit premier récepteur un premier signal srft) représentatif d'une onde se propageant dans la plaque,

- on définit sur la plaque un contour fermé C entourant ledit point P, le contour C étant le lieu de la plaque sur lequel soit on génère par un premier émetteur une onde se propageant dans la plaque, soit on mesure par ledit premier récepteur ledit premier signal srft) représentatif d'une onde se propageant dans la plaque, et

- on détermine le paramètre physique au point P de la plaque en identifiant grâce au moins audit premier signal srft) une fonction de forme f s ha P e(f), dans la relation suivante :

W contour (^) = fshape (f X^ 'plate (/ ~ KS )

Wcontour est une fonction de Green représentant l'onde le long du contour C,

Gpi a te est une fonction de Green représentant l'onde en un point de vecteur position r de la plaque n'appartenant pas au contour C, par rapport à un point S de vecteur position f s de la plaque représentant une source de l'onde, et

ladite fonction de forme f s ha P e(f) est dépendante au moins de la fréquence / de l'onde et du paramètre physique, et adaptée à la forme du contour C.

Grâce à ces dispositions, il est possible de déterminer simplement et précisément un paramètre physique de la plaque au point P.

Dans divers modes de réalisation du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes :

- la fonction de forme est une fonction de Bessel de première espèce Jo(Z) comprenant des zéros Z n , n étant un entier positif ou nul, ladite fonction de Bessel est fonction d'un paramètre d'échelle a multiplié par la racine carrée de la fréquence / de l'onde, de telle sorte que :

- le premier récepteur est adapté pour mesurer une onde sur le contour C, et ledit procédé comprend l'étape suivante :

- si le premier signal si(t) a une amplitude inférieure à un seuil prédéterminé pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

- le dispositif comprend en outre un premier émetteur, et l'un des premier émetteur et du premier récepteur étant adapté soit pour générer ou mesurer une onde sur le contour C, l'autre étant adapté soit pour générer ou mesurer une onde en un point de la plaque n'appartenant pas au contour C, et ledit procédé comprend les étapes suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur,

- on mesure par le premier récepteur un premier signal si(t) représentatif de ladite première onde émise, et

- si le premier signal si(t) a une amplitude inférieure à un seuil prédéterminé pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

a =-

If n

- le dispositif comprend en outre un deuxième récepteur, le premier récepteur étant adapté pour mesurer une onde sur le contour C et le deuxième récepteur étant adapté pour mesurer une onde en un point de la plaque n'appartenant pas au contour C, et ledit procédé comprend les étapes suivantes :

- on mesure simultanément par le premier récepteur un premier signal srft) représentatif d'une première onde, et par le deuxième récepteur un deuxième signal S 2 (t) représentatif de ladite même première onde ;

- le dispositif comprend en outre :

- un deuxième récepteur, et

- un premier émetteur,

l'un des premier récepteur, deuxième récepteur et premier émetteur étant adapté pour mesurer ou générer une onde sur le contour C et ledit procédé comprend les étapes suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission ej(t) alimentant le premier émetteur,

- on mesure simultanément par le premier récepteur un premier signal srft) représentatif de ladite première onde émise, et par le deuxième récepteur un deuxième signal S 2 (t) représentatif de ladite même première onde émise.

- le dispositif comprend en outre :

- un premier émetteur, et

- un deuxième émetteur,

l'un des premier récepteur, premier émetteur et deuxième émetteur étant adapté pour mesurer ou générer une onde sur le contour C et ledit procédé comprend les étapes suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur,

- on mesure par le premier récepteur (Rl) un premier signal sj(t) représentatif de ladite première onde émise,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission e 2 (t) alimentant le deuxième émetteur,

- on mesure par le premier récepteur un deuxième signal S 2 (t) représentatif de ladite deuxième onde émise ;

- le deuxième signal d'émission e 2 (t) est déphasé de π/2 par rapport au premier signal d'émission erft), et ledit procédé comprend les étapes suivantes :

- on calcule un signal sommé s(t), qui est la somme du premier signal srft) et d'un deuxième signal S 2 (X), et

- si le premier signal srft) est en phase avec le signal sommé s(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

- le dispositif comprend en outre :

- un premier émetteur, et

- un deuxième émetteur,

l'un des premier récepteur, premier émetteur et deuxième émetteur étant adapté pour mesurer ou générer une onde sur le contour C et ledit procédé comprend les étapes suivantes :

- on génère simultanément dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur, et une deuxième onde par un deuxième signal d'émission e 2 ftj alimentant le deuxième émetteur, et - on mesure par le premier récepteur un premier signal srft) représentatif de la superposition sur le lieu du premier récepteur desdites première et deuxième ondes émises ;

- dans le procédé :

- le deuxième signal d'émission β 2 (t) est déphasé de π/2 par rapport au premier signal d'émission eι(t), et dans lequel ledit procédé comprend les étapes suivantes :

- si le premier signal srft) est en phase avec le premier signal d'émission ei(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

- le procédé comprend les étapes suivantes :

- on calcule un signal sommé s(t), qui est la somme du premier signal srft) et d'un deuxième signal déphasé S 2 * (t), le deuxième signal déphasé S 2 * (t) étant égal au deuxième signal S 2 (t) déphasé de π/2, et

- si le premier signal srft) est en phase avec le signal sommé s(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, alors on calcule le paramètre d'échelle a par :

- le procédé comprend les étapes suivantes :

- on calcule une première transformée de Fourier Si(f) du premier signal srft) et une deuxième transformée de Fourier S 2 {f) du deuxième signal S2(t),

- on calcule une fonction de test ft es t(f), qui compare le signe de la partie réelle de la première transformée de Fourier Si(f) au signe de la partie réelle de la deuxième transformée de Fourier ^f/J, et qui affecte une première valeur Vi si les signes sont identique et une deuxième valeur V 2 si les signes sont différents :

- on recherche des fréquences particulières /„ auxquelles la fonction de test f tes t(f) change de valeur, en passant soit de la première valeur Vi à la deuxième valeur

V 2 , soit réciproquement de la deuxième valeur V 2 à la première valeur Vi, et

- on calcule le paramètre d'échelle a par

Z n

a =-

If n

- le procédé comprend les étapes suivantes :

- on calcule une première transformée de Fourier Si(J) du premier signal Si(t) et une deuxième transformée de Fourier S2(J) du deuxième signal S2(t),

- on calcule une différence de phase Δφ entre la première transformée de Fourier et la deuxième transformée de Fourier, par :

Δφ=C(S 2 (T)-S 1 (Z))

- on recherche des fréquences particulières /„ de la différence de phase Δφ, pour lesquelles ladite différence de phase présente un saut entre 0 et π ou entre π et 0, et qui sont proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, et

- on calcule le paramètre d'échelle a par :

a=^= ;

If n

- le procédé comprend les étapes suivantes :

- on calcule une première transformée de Fourier Si(J) du premier signal srft) et une deuxième transformée de Fourier S2(J) du deuxième signal S2(t),

- on détermine le paramètre d'échelle a de telle sorte que le module de la fonction de forme suivante :

où b est un autre paramètre d'échelle, et 1.1 est la fonction module,

s'approche au mieux de :

\s 2 (fys!(f)\

pour un ensemble de fréquences de test /„ ;

- le procédé comprend les étapes suivantes :

- on calcule une première transformée de Fourier Si(f) du premier signal srft) et une deuxième transformée de Fourier S 2 (f) du deuxième signal S2(t),

- on détermine le paramètre d'échelle a de telle sorte de la fonction de forme suivante :

où b est un autre paramètre d'échelle, et

φ(.) est la fonction de phase,

s'approche au mieux de :

C(S 2 (T)ZS 1 (J))

pour un ensemble de fréquences de test /„ ;

- on détermine le paramètre physique de produit V P h, produit de l'épaisseur par la vitesse de propagation d'une onde dans la plaque, par la formule suivante :

a est le paramètre d'échelle de la fonction de Bessel, précédemment déterminé, et

R est la longueur d'un segment entre le point P et un point du contour C dans la direction de l'onde ;

- on détermine le paramètre physique d'épaisseur h de la plaque par la formule suivante :

a est le paramètre d'échelle de la fonction de Bessel, précédemment déterminé,

R est la longueur d'un segment entre le point P et un point du contour C dans la direction de l'onde, et

Vp est la vitesse de propagation connue d'une onde dans le matériau de la plaque ;

- on détermine le paramètre physique de vitesse de propagation Vp d'une onde dans la plaque par la formule suivante :

a est le paramètre d'échelle de la fonction de Bessel, précédemment déterminé,

R est la longueur d'un segment entre le point P et un point du contour C dans la direction de l'onde, et

h est l'épaisseur connue de la plaque ;

- le contour C est sensiblement un cercle de rayon

R centré sur le point P ;

- le contour C est sensiblement une ellipse centrée sur le point P ;

- on détermine au préalable la forme du contour C grâce à :

- un dispositif de test comprenant :

- un premier émetteur adapté pour générer une onde au point P,

- au moins un deuxième émetteur adapté pour générer une onde sur un contour de test C n ayant la forme d'une ellipse prédéterminée, n étant un indice entier positif,

- des premier et deuxième récepteurs adaptés pour mesurer une onde en des points n'appartenant pas au contour de test C n , et grâce à :

- un procédé de test comprenant les étapes de test suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur,

- on mesure par le premier récepteur un premier signal Su(t) représentatif de ladite première onde émise, et on calcule une première transformée de Fourier Su(J) de ce premier signal,

- on mesure par le deuxième récepteur un deuxième signal Sn(t) représentatif de ladite première onde émise, et on calcule une deuxième transformée de Fourier Sn(J) de ce deuxième signal,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission ez(t) alimentant le deuxième émetteur,

- on mesure par le premier récepteur (Rl) un troisième signal S 2 i(t) représentatif de ladite deuxième onde émise, et on calcule une troisième transformée de Fourier S21O) de ce troisième signal,

- on mesure par le deuxième récepteur un quatrième signal S 22 (t) représentatif de ladite deuxième onde émise, et on calcule une quatrième transformée de Fourier S 22 O) de ce quatrième signal,

- on calcule la fonction de différence de phase suivante :

Aφ{f)= φiSjjif\S 12 (f)*)-φ(S 21 (f).S 22 (f)*)

* désigne la fonction conjuguée, et

φ(.) est la fonction de phase, et

- la forme de l'ellipse du contour de test C n correspond à un contour optimal, telle que la première onde se propage et se superpose spatialement sur la plaque sensiblement à la deuxième onde, lorsque la fonction de différence de phase Δφ(f) est minimale pour un ensemble de contours de test. C n auxquels on applique les étapes de test précédentes ;

- le contour C est sensiblement un rectangle centré sur le point P ;

- le contour C comporte huit points de contour Cl à C8, et dans lequel lesdits points de contour et le point P forment un maillage rectangle régulier ; - le segment de longueur R est un segment de longueur moyenne calculé par :

R _ d / +d 2

2

di est la moitié de la longueur de la plus grande médiane dudit rectangle, et

d 2 est la longueur de la diagonale dudit rectangle.

L'invention se rapporte également à un procédé d'imagerie, dans lequel on construit une image d'une plaque, ladite image comportant une pluralité de pixels, et chaque pixel correspondant à un point de la plaque et représentant un paramètre physique de la plaque audit point de la plaque, et ledit paramètre physique dudit point étant déterminé par le procédé défini précédemment.

L'invention se rapporte également à un dispositif pour mettre en œuvre le procédé pour déterminer un paramètre physique représentatif d'un point P d'une plaque selon l'une des revendications précédentes, dans lequel le paramètre physique est choisi parmi l'épaisseur de la plaque h, la vitesse de propagation d'une onde dans la plaque V P et le produit Vph de l'épaisseur par la vitesse de propagation d'une onde dans la plaque, et ledit dispositif comprenant :

- au moins un premier récepteur adapté pour réaliser une mesure d'un premier signal srft) représentatif d'une onde se propageant dans la plaque,

- un contour fermé C défini en entourant ledit point P, le contour C étant le lieu de la plaque sur lequel soit on génère par un premier émetteur une onde se propageant dans la plaque, soit on mesure par ledit premier récepteur ledit premier signal srft) représentatif d'une onde se propageant dans la plaque, et

- une unité de calcul reliée audit premier récepteur ladite unité de calcul étant adaptée pour déterminer le paramètre physique au point P de la plaque en identifiant grâce au moins audit premier signal srft) une fonction de forme f s ha P e(f), dans la relation suivante :

^contour(?) = fshape(/)°plate(? ~ ?S )

Wcontour est une fonction de Green représentant l'onde le long du contour C,

Gpi a te est une fonction de Green représentant l'onde en un point de vecteur position r de la plaque n'appartenant pas au contour C, par rapport à un point S de vecteur position f s de la plaque représentant une source de l'onde, et

ladite fonction de forme f s ha P e(f) est dépendante au moins de la fréquence / de l'onde et du paramètre physique, et adaptée à la forme du contour C.

Dans divers modes de réalisation du dispositif selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes :

- le premier récepteur est un vibromètre laser à balayage ;

- le dispositif comprend en outre un deuxième récepteur, et dans lequel le deuxième récepteur est réalisé par ledit vibromètre laser à balayage.

D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'un de ses modes de réalisation, donné à titre d'exemple non limitatif, en regard des dessins joints.

Sur les dessins :

- la figure 1 est une vue d'une plaque sur laquelle on peut mettre en œuvre l'invention ;

- la figure 2 représente une fonction de Bessel de première espèce Jo(Z) ;

- la figure 3 est un premier mode de réalisation d'un dispositif selon l'invention ; - les figures 4a et 4b sont un deuxième mode de réalisation d'un dispositif selon l'invention ;

- la figure 5 est un troisième mode de réalisation d'un dispositif selon l'invention ;

- la figure 6a et 6b sont un quatrième mode de réalisation d'un dispositif selon l'invention ;

- la figure 7a, 7b et 7c sont un cinquième mode de réalisation d'un dispositif selon l'invention ;

- la figure 8 représente un premier transducteur adapté pour mettre en œuvre l'invention ;

- la figure 9 représente un deuxième transducteur adapté pour mettre en œuvre l'invention ;

- les figures 10 à 12 illustrent des variantes du procédé selon l'invention ;

- la figure 13 représente une image réalisée à l'aide du procédé selon l'invention.

Dans la suite du document, le terme « vibration » sera compris comme pouvant désigner une onde vibratoire, une onde acoustique, ou une onde ultrasonore. L'onde en question a une fréquence par exemple comprise entre 100 Hz et 50 kHz, et de préférence entre 1000 Hz et 20 kHz, de sorte que des matériels peu coûteux permettent de mesurer une telle onde.

L'invention concerne un procédé pour déterminer un paramètre physique représentatif d'un point P d'une plaque 1, dans lequel le paramètre physique est choisi notamment parmi l'épaisseur de la plaque h, la vitesse de propagation d'une onde dans la plaque V P et le produit Vph de l'épaisseur par la vitesse de propagation d'une onde dans la plaque.

Le procédé est mis en œuvre par un dispositif comprenant :

- au moins un récepteur Rl adapté pour réaliser au moins une mesure d'une onde sur ladite plaque, et

- une unité de calcul CALC reliée audit récepteur, adaptée pour obtenir ladite mesure dudit récepteur, et adaptée pour déterminer ledit paramètre physique au point P à partir de ladite mesure.

On définit sur la plaque un point P correspondant au lieu où ledit paramètre physique est à déterminer et un contour fermé C entourant ledit point P.

Le procédé est basé sur l'émission par un émetteur d'une onde sur le contour C ou réciproquement la réception par un récepteur d'une onde sur le contour C. Le principe de réciprocité de la propagation d'une onde acoustique ou vibratoire dans une structure, conduit à plusieurs modes de réalisation de ce procédé, dans lesquels des éléments sont soit des émetteurs, soit des récepteurs.

L'onde peut être une onde vibratoire, acoustique ou ultrasonore. Elle se propage dans la plaque ou sur une surface de ladite plaque. Cette onde peut être mesurée par un récepteur, ou générée par un émetteur.

Un émetteur ou récepteur d'une onde peut être un transducteur, tel que par exemple un dispositif à matériau piézoélectrique fixé sur la plaque. Dans un mode récepteur, le transducteur transforme un déplacement, une déformation, une contrainte ou une pression en une tension électrique, représentant une mesure dudit déplacement, de ladite déformation, de ladite contrainte ou de ladite pression. La tension électrique peut être transformée par un convertisseur analogique-numérique, pour fournir une valeur numérique à une unité de calcul. Réciproquement, dans un mode émetteur, le transducteur transforme une tension électrique en un déplacement, une déformation, une contrainte ou une pression. La tension électrique peut être produite par un convertisseur numérique-analogique d'une unité de calcul, et éventuellement suivie d'un amplificateur de tension.

Alternativement, l'émetteur ou le récepteur peut être sans contact avec la plaque. Par exemple, un transducteur électromagnétique ou un laser de forte puissance fonctionnant en régime impulsionnel sont utilisables. Réciproquement, on peut utiliser un récepteur optique, tel qu'un vibromètre laser, adapté pour mesurer sans contact et à distance une vibration d'un point ou d'une multitude de points sur une plaque.

D'autre part, un émetteur ou un récepteur sur un contour C peut être réalisé de multiples manières.

Selon une première possibilité, on utilise des transducteurs piézoélectriques. On fixe ou place sur la plaque un nombre T prédéterminé de transducteurs T 1 , i étant un indice entier positif compris entre 1 et T , placés régulièrement ou non le long du contour C, comme cela est représenté en figure 8. Par exemple, T peut être égal à 3, 4, 8 ou 16. On comprendra que plus le contour C a une longueur importante, plus le nombre T de transducteurs devra être élevé pour obtenir un signal ou générer une onde équivalente à un transducteur continu le long du contour C. Chaque transducteur piézoélectrique comporte une anode et une cathode. Tous les transducteurs sont alors connectés entre eux en parallèle, c'est-à-dire que toutes les anodes sont reliées entre elles par un premier conducteur 12, et toutes les cathodes sont reliées entre elles par un deuxième conducteur 13. L'ensemble des transducteurs T 1 forme donc un transducteur unique, alimenté par uniquement deux conducteurs 12, 13. Il peut alors être connecté à un dispositif de la même manière qu'un unique transducteur, tel que le transducteur P alimenté par deux autres conducteurs 10, 11. En variante (non représentée), au lieu d'être connectés entre eux en parallèle, les transducteurs pourraient être montés en série, ce qui donnerait un résultat équivalent.

Selon une deuxième possibilité représentée en figure 9, on peut également réaliser un transducteur adapté pour mesurer une onde le long d'un contour C, à l'aide d'un matériau polymère piézoélectrique, tel que du polyfluorure de vinylidène (PVDF) . Ce matériau est souple et peut être mis en forme sur un film adhésif 2, adapté pour être collé sur la plaque 1 par sa face inférieure. Le film 2 comprend donc une pastille P en PVDF entouré d'un contour C circulaire en PVDF. La pastille P est reliée à une première borne 20 par un conducteur 10, et à une deuxième borne 21 par un conducteur 11. Le contour C est relié à ladite première borne 20 et à une troisième borne 22 par un troisième conducteur 12. Les bornes de connexion sont par exemple situées sur la périphérie du film 2, et présentent chacune une face conductrice supérieure du côté de la face supérieure du film 2 et adaptée à être connectée à un dispositif .

Un tel dispositif intègre donc un premier et deuxième récepteur Rl, R2 tel que cela peut être souhaité pour mettre en œuvre le procédé de l'invention. De plus, ce dispositif peut être peu coûteux et facile à mettre en œuvre .

Selon une troisième possibilité, on utilise un vibromètre (non représenté) . Le vibromètre fait un balayage du contour C avec un pas prédéterminé pour produire la mesure de chaque point sur ledit contour C. La mesure sur le contour C sera alors calculée numériquement en sommant les signaux de chaque point.

Le même vibromètre peut être utilisé pour effectuer des mesures d'un ou plusieurs autres points sur la plaque n'appartenant pas au contour C, de telle sorte qu'un seul dispositif de mesure sera utilisé pour réaliser l'ensemble des mesures nécessaires.

Lorsqu'il est énoncé dans cette demande de brevet que l'on émet ou que l'on reçoit une onde sur un contour C, on comprendra que l'on peut positionner sur la plaque un nombre prédéterminé de transducteurs connectés entre eux pour mesurer ou générer cette onde sur l'ensemble du contour C, ou que l'on peut utiliser un capteur sans contact, tel qu'un vibromètre à balayage, pour réaliser cette fonction, ou tout autre moyen connu.

Selon une quatrième possibilité (non représentée) , on utilise des transducteurs piézoélectriques du type de ceux de la première possibilité de la figure 8, par exemple en céramique, montés sur un film souple, par exemple en plastique, de sorte que les transducteurs T 1 sont placés à des positions prédéterminées les uns par rapport aux autres en formant un contour C fermé, autour d'un transducteur P à l'intérieur dudit contour C. Le nombre T de transducteurs peut être de 3, 4, 8 ou 16, ou plus si le contour a une longueur importante.

Le film peut avoir une face revêtue d'un adhésif pour fixer directement tous les transducteurs sur la plaque. L'adhésif ne sera cependant pas présent sur tout le film, mais uniquement localisé sous les transducteurs et le film doit avoir une faible élasticité pour que l'ensemble ne modifie pas localement la réponse vibratoire de la plaque .

En variante, les transducteurs peuvent avoir une face revêtue d'un adhésif pour être directement fixé sur la plaque .

Les transducteurs T 1 sont reliés entre eux par des conducteurs électriques souples formés sur le film, et selon un même principe que pour le cas de la première possibilité. L'ensemble des transducteurs T 1 , du transducteur P, du film et des conducteurs souples forme donc un produit comprenant un premier et deuxième récepteur Rl, R2 pour mettre en œuvre le procédé. De plus, cet ensemble est prêt à être placé sur une plaque rapidement et facilement . Les fondements théoriques permettant de comprendre les divers modes de réalisations du dispositif et du procédé de l'invention sont décrits ci-dessous. On se reportera à la figure 1 représentant une plaque 1 et permettant d'illustrer cette partie.

Par une approche mécanique de plaque mince dans le domaine fréquentiel, le rayonnement d'un point de la plaque est régi par l'équation suivante :

r vecteur position du point de la plaque, de préférence en coordonnées polaires,

(O pulsation de l'onde où avec / fréquence de l 'onde,

p : la masse volumique du matériau de la plaque, et

h : l'épaisseur de la plaque,

S: fonction de Dirac représentant une source ponctuelle centrée à l'origine des coordonnées. et D = Ehr

12(1-^) '

E est le module de Young du matériau de la plaque, et

(J est le coefficient de Poisson du matériau. Pour une plaque infinie, la solution de cette équation est une fonction de Green Gf ree pour un point de coordonnée r et avec une source vibratoire ou acoustique placé au point S de la plaque et de coordonnées r s :

£ 4 = , (4)

H 0 est la fonction de Hankel de première espèce.

Pour une plaque de dimension finie, l'onde vibratoire résulte également de l'interférence avec de multiples ondes réfléchies sur les bords de la plaque, de telle sorte que la fonction de Green G p ι ate pour le point de coordonnées r sur une plaque de dimension finie, peut s ' écrire :

G plate { \r - r s \}) ( 5 ) où

C refl est une fonction qui représente uniquement les réflexions sur les bords de la plaque. Cette fonction dépend de la fonction de Green Gf ree sur une plaque de dimension infinie. Cette fonction est linéaire. Le rayonnement du contour C peut être calculé par une somme de rayonnements ponctuels le long de ce contour, chacun étant calculé par la formulation précédente. Dans le cas d'un contour C sous la forme d'un cercle, l'utilisation du théorème d'addition des harmoniques cylindriques permet d'obtenir la fonction de Green W cιrc ι e pour un contour C ayant la forme cercle sur une plaque de dimension finie :

W circle ( F) = 2πAJ 0 , (6) où

A est une amplitude, et

Jo est une fonction de Bessel de première espèce dépendant produit d'un nombre d'onde k et du rayon du cercle R.

En outre :

où :

Vp est le paramètre de vitesse de propagation de l'onde dans la plaque,

/ est la fréquence de l'onde, et

h est l'épaisseur de la plaque.

On obtient donc le produit kR suivant :

Ainsi, le produit Vph, produit de l'épaisseur de plaque h par la vitesse de propagation de l'onde dans la plaque Vp, s'écrit : (ut?

La fonction de Bessel Jo(Z) est une fonction oscillante, représentée en figure 2. Cette fonction s'annule ou présente des zéros ou racines pour des valeurs d'abscisse Z n particulières, n étant un indice entier positif ou nul. Les cinq premiers zéros peuvent être notés

Zo, Zi, X2, Z3, Z4, Z5 et ils ont pour valeurs approchées :

Z 0 ≈ 2,4048

Z 1 ≈ 5,5201

Z 2 ≈ 8,6537

(10) Z 3 ≈ 117915

Z 4 ≈ 14,9309

Z 5 ≈ 18,0711

Les zéros de la fonction de Bessel Jo sont espacés l'un de l'autre d'une manière périodique, de telle sorte que, connaissant la fréquence d'une onde, on peut déterminer des produits JcR correspondant à chaque zéro Z n , et en déduire le produit Vph, produit de l'épaisseur h par la vitesse de propagation de l'onde dans la plaque Vp. On en déduit donc un paramètre physique de la plaque.

Divers modes de réalisation du dispositif sont possibles, chacun ayant plusieurs variantes.

Dans un premier mode de réalisation du dispositif représenté en figure 3, le dispositif comporte un seul récepteur Rl. Ce récepteur ou capteur de vibration Rl est adapté pour mesurer une onde sur le contour C de la plaque.

Le contour C est éventuellement un cercle de rayon R centré sur un point P. Ce premier mode de réalisation ne comporte pas d'émetteur. C'est donc un dispositif passif, qui utilise le bruit et/ou les vibrations de l'environnement du dispositif .

Dans un deuxième mode de réalisation du dispositif, le dispositif comporte un seul récepteur Rl et un seul émetteur El. Si le récepteur ou capteur Rl est adapté pour mesurer une onde sur le contour C de la plaque, l'émetteur El est adapté pour générer une onde en un point quelconque de la plaque n'appartenant pas au contour C (figure 4a). Réciproquement, si l'émetteur El est adapté pour générer une onde sur le contour C de la plaque, le récepteur Rl est adapté pour mesurer une onde en un point quelconque de la plaque n'appartenant pas au contour C (Figure 4b) . Le contour C est éventuellement un cercle de rayon R centré sur un point P. Ce deuxième mode de réalisation comporte un émetteur, c'est donc un dispositif actif.

Dans un troisième mode de réalisation du dispositif, le dispositif comporte deux récepteurs Rl, R2, mais pas d'émetteur. La figure 5 montre un exemple de dispositif de ce troisième mode de réalisation. Un premier récepteur Rl est adapté pour mesurer une onde sur le contour C de la plaque. Le contour C est éventuellement un cercle de rayon R centré sur un point P. Un deuxième récepteur R2 est adapté pour mesurer une onde en un point quelconque de la plaque n'appartenant pas au contour C. Notamment, ce deuxième récepteur R2 peut être adapté pour mesurer une onde au point P ou à un point intérieur au contour C, c'est-à-dire entouré par le contour C, ou à un point extérieur au contour C. Ce troisième mode de réalisation ne comporte pas d'émetteur, c'est donc un dispositif passif.

Dans un quatrième mode de réalisation du dispositif, le dispositif comporte deux récepteurs Rl, R2 et un émetteur El. Les figures 6a et 6b montrent un exemple de dispositif de ce quatrième mode de réalisation. L'émetteur El, les récepteurs Rl, R2 peuvent être adaptés pour mesurer ou générer une onde, indifféremment sur le contour C ou sur un premier point quelconque de la plaque ou sur un deuxième point quelconque de la plaque. Le contour C est éventuellement un cercle de rayon R centré sur un point P. Ce quatrième mode de réalisation comporte un émetteur, c'est donc un dispositif actif. Dans un cinquième mode de réalisation du dispositif, le dispositif comporte deux émetteurs El, E2 et un récepteur Rl . Les figures 7a, 7b et 7c montrent un exemple de dispositif de ce cinquième mode de réalisation. Le récepteur Rl, les émetteurs El, E2 peuvent être adaptés pour mesurer ou générer une onde, indifféremment sur le contour C ou sur un premier point quelconque de la plaque ou sur un deuxième point quelconque de la plaque. Le contour C est éventuellement un cercle de rayon R centré sur un point P. Ce cinquième mode de réalisation comporte deux émetteurs, c'est donc un dispositif actif.

Divers modes de réalisation du procédé sont possibles, chacun étant adapté à un ou plusieurs des modes de réalisation du dispositif. Ces modes de réalisations du procédé sont décrits ci-dessous.

Dans un premier mode de réalisation du procédé adapté notamment aux premier et deuxième modes de réalisation du dispositif comprenant uniquement un seul récepteur Rl, ledit procédé pour déterminer le paramètre physique comprend alors les étapes suivantes :

- on mesure par le récepteur Rl sur le contour C un signal srft) représentatif de la propagation d'une onde dans la plaque.

Le signal srft) présente une amplitude nulle pour certaines fréquences, notamment pour des antirésonances de la structure de la plaque, mais également pour des fréquences particulières de la fonction de forme.

Dans le cas d'un contour sous la forme d'un cercle et d'un matériau sensiblement isotrope et d'après les équations (6) et (8), la fonction de forme est une fonction de Bessel de première espèce Jo. Cette fonction de Bessel est fonction d'un paramètre d'échelle a multiplié par la racine carrée de la fréquence / de l'onde :

el le s ' annule pour les zéros Z n , Z n = CIyJf Dans les autres cas dans lesquels le matériau n'est pas isotrope ou le contour C n'est pas un cercle, on peut déterminer une fonction de forme numériquement, ayant également des zéros pour certaines fréquences particulières.

On considère un ensemble de fréquences de test /„, n étant un indice entier positif ou nul compris entre zéro et N, N également entier positif, N pouvant par exemple être égal à cinq. Les fréquences de tests /„ sont définies comme proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo. Si pour les fréquences de test /„, le signal srft) a une amplitude faible ou nulle, et par exemple, inférieure à un seuil prédéterminé S, alors ces fréquences de test /„ correspondent bien aux zéros de la fonction de Bessel Jo et on peut calculer le paramètre d'échelle a de la fonction de Bessel par :

a = r — pour tout n entre zéro et N.

On peut alors calculer un paramètre physique.

Notamment d'après l'équation (9), le produit Vph, produit de l'épaisseur h par la vitesse de propagation Vp d'une onde dans la plaque au point P peut être calculé par :

pour n compris entre 0 et N,

où /„ est la fréquence de test d'indice n de 1 'ensemble,

Z n est le zéro d'indice n de la fonction de Bessel de première espèce Jo, ledit zéro Z n correspondant à ladite fréquence de test /„ de même indice, et

R est le rayon du contour C.

Si on connait la valeur de la vitesse de propagation Vp d'une onde dans le matériau de la plaque, on peut calculer l'épaisseur de la plaque au point P par :

pour n compris entre 0 et N.

Si on connait la valeur de l'épaisseur de la plaque h, on peut calculer la vitesse de propagation d'une onde dans la plaque par :

pour n compris entre 0 et N.

Selon un deuxième mode de réalisation du procédé adapté notamment aux troisième et quatrième modes de réalisation du dispositif de l'invention comprenant au moins deux récepteurs Rl, R2 et éventuellement un émetteur

El, l'un d'entre eux étant sur le contour C (figures 5a, 6a et 6b) , ledit procédé pour déterminer le paramètre physique comprend alors les étapes suivantes :

- on mesure par le premier récepteur Rl un premier signal srft),

- on mesure par le deuxième récepteur R2 un deuxième signal S2(t),

- on déphase de π/2 le deuxième signal.

Par conséquent, si le premier signal est du type

alors d'après l'équation 6, le deuxième signal déphasé s*2(t) peut s'écrire :

s* 2 {ή= AJ 0 {kR)sin{2φ)

Soit s(t) la somme de si(t) et s*2(t) .

En posant, tan(φ) = AJo(JcR), on obtient :

(A _cos{2≠-φ)

Λ V) ~ ( \

cosyφ)

Par conséquent, pour les zéros de la fonction de Bessel de première espèce Jo, tan(φ) = 0. Donc, φ = 0. Dans ces conditions, le signal sommé s(t) est en phase avec le premier signal srft).

Pour déterminer si un signal est en phase avec un autre, toute technique pourra être employée, dans le domaine temporel ou fréquentiel.

Le deuxième mode de réalisation du procédé comprend alors également une étape dans laquelle des fréquences de tests /„ sont définies comme proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo. Si pour les fréquences de test /„, le signal srft) sensiblement en phase avec un signal sommé s(t) correspondant à la somme du premier signal et du deuxième signal déphasé de π/2, alors ces fréquences de test /„ correspondent bien aux zéros de la fonction de Bessel Jo.

On peut alors calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus.

Ce deuxième mode de réalisation du procédé est également utilisable avec le cinquième mode de réalisation du dispositif (figures 7a, 7b et 7c) comprenant deux émetteurs El, E2 et un seul récepteur Rl.

Dans ce cas :

- on génère dans la plaque une première onde par un premier signal d'émission erft) alimentant le premier émetteur El,

- on mesure par le premier récepteur Rl, un premier signal srft) représentatif de ladite première onde émise,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission e 2 (X) alimentant le deuxième émetteur E2,

- on mesure par le premier récepteur Rl, un deuxième signal S 2 (t) représentatif de ladite deuxième onde émise .

Ensuite, on déphase de π/2 le deuxième signal S 2 (t) pour former un deuxième signal déphasé s* 2 (t) . La suite du procédé est alors identique à ce qui est décrit ci-dessus.

Les deux variantes décrites du deuxième mode de réalisation du procédé utilisent donc le déphasage du deuxième signal à la réception.

Selon un troisième mode de réalisation du procédé, on réalise un déphasage à l'émission. Ce troisième mode de réalisation du procédé est notamment adapté au cinquième mode de réalisation du dispositif comprenant deux émetteurs El, E2 et un récepteur Rl (figures 7a, 7b et 7c) . Le procédé pour déterminer le paramètre physique comprend alors les étapes suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission erft) alimentant le premier émetteur El,

- on mesure par le premier récepteur Rl, un premier signal srft) représentatif de ladite première onde émise,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission e 2 (X) alimentant le deuxième émetteur E2,

- on mesure par le premier récepteur Rl, un deuxième signal S 2 (t) représentatif de ladite deuxième onde émise .

On calcule alors un signal sommé s(t), qui est la somme du premier signal srft) et d'un deuxième signal S 2 (t) :

s{t)= Sl {t)+s 2 {t).

Si le premier signal srft) est en phase avec le signal sommé s(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, alors ces fréquences correspondent aux fréquences particulières recherchées.

On peut alors calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus.

En variante, les émissions des premier et deuxième émetteurs sont effectuées en simultané. Dans ce cas, le procédé pour déterminer le paramètre physique comprend alors les étapes suivantes :

- on génère simultanément dans la plaque une première onde par un premier signal d'émission erft) alimentant le premier émetteur El, et une deuxième onde par un deuxième signal d'émission β 2 (t) alimentant le deuxième émetteur E2, et

- on mesure par le premier récepteur Rl un premier signal srft) représentatif de la superposition sur le lieu du premier récepteur Rl desdites première et deuxième ondes émises.

Le deuxième signal d'émission e 2 W est déphasé de π/2 par rapport au premier signal d'émission ei(t).

Si le premier signal srft) est en phase avec le premier signal d'émission eι(t), pour un ensemble de fréquences de test /„, les fréquences de test /„ étant proportionnelles au carré des zéros Z n de la fonction de Bessel de première espèce Jo, alors ces fréquences correspondent aux fréquences particulières recherchées.

On peut alors calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus.

Selon un quatrième mode de réalisation du procédé, adapté notamment aux troisième et quatrième modes de réalisation du dispositif de l'invention comprenant au moins deux récepteurs Rl, R2 et éventuellement un émetteur El, l'un d'entre eux étant sur le contour C (figures 5, 6a et 6b) , ledit procédé pour déterminer le paramètre physique comprend les étapes suivantes :

- on mesure par le premier récepteur Rl un premier signal Sι(t), et on calcule une première transformée de

Fourier Si(f) de ce premier signal,

- on mesure par le deuxième récepteur R2 un deuxième signal S 2 (t), et on calcule une deuxième transformée de Fourier S 2 Φ de ce deuxième signal.

Le premier signal srft) et le deuxième signal S 2 (t) sont en phase ou en opposition de phase pour des fréquences particulières correspondant aux abscisses pour lesquelles la fonction de Bessel Jo présente un zéro.

Ces fréquences peuvent alors être identifiées : - soit directement en comparant le signe des parties réelles des transformées de Fourier, - soit indirectement en calculant une différence de phase .

Dans le premier cas, le signe de la partie réelle de la première transformée de Fourier Si(f) est comparé au signe de la partie réelle de la deuxième transformée de

Fourier ^ 2 (X). Il faut observer les bandes de fréquences dans lesquelles les signes sont identiques et les bandes de fréquences dans lesquelles les signes sont opposés. Les fréquences de transition entre ces bandes de fréquences permettent d'identifier les fréquences particulières des zéros de la fonction de Bessel.

On calcule, par exemple, une fonction de test f tes t(f) suivante qui prenant une première valeur V / si les signes sont identique et une deuxième valeur V 2 si les signes sont différents :

\ si sign(^(S j(f))) = sign(^(S 2 (f))) alors f test (f) = V j

sinon f test {f) = V 2

où Vl et V2 peuvent prendre des valeurs quelconques différentes. Par exemple, Vi = I et V 2 = 0.

On recherche des fréquences particulières /„ auxquelles la fonction de test f tes t(f) change de valeur, en passant soit de la première valeur V / à la deuxième valeur

V 2 , soit réciproquement de la deuxième valeur V 2 à la première valeur V / .

On peut alors calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus.

Dans le deuxième cas, on calcule une différence de phase Δφ entre la première transformée de Fourier Si(f) et la deuxième transformée de Fourier S 2 (f), par :

Aφ= φ(S 2 (f)-S 1 (f)) ;

La différence de phase Δφ présente alors des sauts de phase entre 0 et π, ou entre π et 0, pour les fréquences particulières recherchées.

On recherche des fréquences particulières /„ de la différence de phase Δφ, pour lesquelles ladite différence de phase Δφ présente un tel saut, lesdites fréquences particulières /„ étant proportionnelles au carré des zéros Z n de la fonction de première espèce Jo.

Pour détecter un saut sur une fonction, telle que la différence de phase, toute technique pourra être employée. Notamment, on pourra détecter le franchissement d'un seuil intermédiaire, proche de π/2, dans le sens croissant ou décroissant, avec ou sans hystérésis.

Une fois les fréquences identifiées, il est alors possible de calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus.

Ce quatrième mode de réalisation du procédé est également utilisable avec le cinquième mode de réalisation du dispositif (figures 7a, 7b et 7c) comprenant deux émetteurs El, E2 et un seul récepteur Rl.

Dans ce cas :

- on génère dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur El,

- on mesure par le premier récepteur Rl, un premier signal srft) représentatif de ladite première onde émise,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission β 2 (t) alimentant le deuxième émetteur E2,

- on mesure par le premier récepteur Rl, un deuxième signal S 2 (t) représentatif de ladite deuxième onde émise .

On calcule une première transformée de Fourier Si(f) du premier signal sj(t) et une deuxième transformée de

Fourier S2Φ du deuxième signal S2(t) .

De même que pour le troisième mode de réalisation du procédé, le premier signal srft) et le deuxième signal S 2 (t) sont en phase ou en opposition de phase pour des fréquences particulières correspondant aux abscisses pour lesquelles la fonction de Bessel Jo présente un zéro.

Dans la suite du procédé, on identifie de même les fréquences particulière soit directement en comparant le signe de la partie réelle de la première transformée de Fourier au signe de la partie réelle de la deuxième transformée de Fourier, soit indirectement en calculant une différence de phase Δφ, la suite du procédé étant identique.

Ayant déterminé les fréquences particulières, on peut calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus. Selon un cinquième mode de réalisation du procédé adapté notamment aux troisième et quatrième modes de réalisation du dispositif de l'invention comprenant au moins deux récepteurs Rl, R2 et éventuellement un émetteur El, l'un d'entre eux étant sur le contour C (figures 5, 6a et 6b) , ledit procédé comprend alors les étapes suivantes :

- on mesure par le récepteur Rl un premier signal Si(t), et on calcule une première transformée de Fourier Si(f) de ce premier signal,

- on mesure par le récepteur R2 un deuxième signal S 2 (t), et on calcule une deuxième transformée de Fourier S 2 (J) de ce deuxième signal.

Lorsque le contour C est un cercle et le matériau de la plaque un matériau sensiblement isotrope, on peut appliquer l'équation 6 pour identifier sur un ensemble de fréquences de test /„ la forme de la fonction de Bessel de première espèce Jo.

Cette identification peut être effectuée :

- soit par les modules,

- soit par les phases.

Dans le premier cas, on recherche les paramètres a et b d'une fonction paramétrique bJoψ-Jf] s 'approchant au mieux de pour un ensemble de fréquences de test /„, 1.1 étant la fonction module.

Dans le deuxième cas, on recherche les paramètres a et b d'une fonction paramétrique (fψJoψ-^jf )) s 'approchant au mieux de φ (S2(f)/ Si(J)) pour un ensemble de fréquences de test fnr φ(-) étant la fonction phase.

Une fois cette identification du paramètre d'échelle a étant effectuée, on a par l'équation 8, le lien entre l'abscisse a-yjf de la fonction de Bessel J 0 et le paramètre physique de la plaque.

On peut calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus.

Ce cinquième mode de réalisation du procédé est également utilisable avec le cinquième mode de réalisation du dispositif (figures 7a, 7b et 7c) comprenant deux émetteurs El, E2 et un seul récepteur Rl.

Dans ce cas :

- on génère dans la plaque une première onde par un premier signal d'émission ei(t) alimentant le premier émetteur El,

- on mesure par le premier récepteur Rl, un premier signal si(t) représentatif de ladite première onde émise,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission β2(t) alimentant le deuxième émetteur E2,

- on mesure par le premier récepteur Rl, un deuxième signal S2(t) représentatif de ladite deuxième onde émise .

On calcule une première transformée de Fourier Si(f) du premier signal si(t) et une deuxième transformée de Fourier S2Φ du deuxième signal S2(t) .

Dans la suite, lorsque le contour C est un cercle et lorsque le matériau de la plaque un matériau sensiblement isotrope, on peut également appliquer l'équation 6 pour identifier sur un ensemble de fréquences de test /„ la forme de la fonction de Bessel de première espèce Jo.

Cette identification est effectuée soit sur les modules, soit sur les phases, comme précédemment, pour obtenir un paramètre d'échelle a.

On peut calculer un paramètre physique de la plaque, avec les formules 11 à 13 explicitées ci-dessus.

Ainsi, on peut appliquer un ou plusieurs modes de réalisation du procédé à chaque mode de réalisation du dispositif, et dans tous les modes de réalisation du procédé on peut déterminer un paramètre physique de la plaque .

Les procédés précédents peuvent être appliqués également si le matériau de la plaque est anisotrope. Dans ce cas, le contour C aura une forme différente d'un cercle.

Dans un premier cas, la vitesse de propagation d'une onde dépend de la direction selon une loi elliptique, du type :

V p 2 [ [_cos 2 θ Ilv p 2 x + sin 2 θ I/v p 2 y 1 J = 1

X est un axe dans la direction du grand axe de

1 'ellipse,

Y est un axe dans la direction du petit axe de

1 'ellipse,

X et Y étant des axes orthogonaux,

θ est l'angle de la direction de la propagation de l'onde par rapport à un axe X,

V px est la vitesse de propagation de l'onde selon l'axe X,

V py est la vitesse de propagation de l'onde selon l'axe Y.

L'équation (6) peut alors s'écrire pour un contour en forme d'ellipse :

W e iu pse ( r) = 2πAJ 0 {kR)G plate ( \r - f s

Si l'on utilise deux récepteurs (Rl, R2 ) à des points de la plaque de coordonnées r} et F 2 , n'appartenant pas au contour en forme d'ellipse, on peut écrire deux fois la précédente relation, et en déduire :

^W ^ (?) W W ^ ( F 2 )* =|2^J 0 (^)| 2 G ptoe (|F 1 -FjG ptoe (|F 2 -Fj* où * désigne la fonction conjuguée.

Par conséquent, la phase de W' e Mp Se {rj)W e ιiip Se {r2)* doit être égale à la phase de G plate -r s \)G plate .

On en déduit un procédé de test dans lequel on utilise un ou plusieurs émetteurs (E2) adaptés pour générer une onde sur un contour de test C n ayant une forme en ellipse prédéterminée, et n indice entier positif, et on applique sur chacun d'eux un procédé de test destiné à déterminé le meilleur contour de test C n , c'est-à-dire la forme de l'ellipse, et comportant les étapes suivantes :

- on génère dans la plaque une première onde par un premier signal d'émission ej(t) alimentant le premier émetteur (El ) ,

- on mesure par le premier récepteur (Rl) un premier signal Su(t) représentatif de ladite première onde émise, et on calcule une première transformée de Fourier Sn(f) de ce premier signal,

- on mesure par le deuxième récepteur (R2) un deuxième signal Sn{t) représentatif de ladite première onde émise, et on calcule une deuxième transformée de Fourier Sn(I) de ce deuxième signal,

- on génère dans la plaque une deuxième onde par un deuxième signal d'émission e 2 (X) alimentant le deuxième émetteur (E2) ,

- on mesure par le premier récepteur (Rl) un troisième signal S 2 i(t) représentatif de ladite deuxième onde émise, et on calcule une troisième transformée de Fourier S2i(f) de ce troisième signal,

- on mesure par le deuxième récepteur (R2) un quatrième signal S 22 (t) représentatif de ladite deuxième onde émise, et on calcule une quatrième transformée de Fourier S22H) de ce quatrième signal,

- on calcule la fonction de différence de phase suivante :

Aφ{f)=φ{S 11 {f).S 12 {f)*)-φ{S 21 {f).S 22 {f)*)

où * désigne la fonction conjuguée, et φ(.) est la fonction de phase.

Alors, la forme de l'ellipse du contour de test C n correspond à un contour optimal, telle que la première onde se propage et se superpose spatialement sur la plaque sensiblement à la deuxième onde, lorsque la fonction de différence de phase Δφ(f) est minimale pour un ensemble de contours de test. C n auxquels on applique les étapes de test précédentes .

Grâce à ce procédé de test, on détermine la forme optimale du contour C à utiliser dans le procédé de l'invention pour lequel toutes les équations établies pour un cercle sont maintenant utilisable pour l'ellipse.

Dans un cas plus général, dans lequel la vitesse de propagation d'une onde dépend de la direction selon une loi ayant une forme prédéterminée, le contour C à utiliser dans tous les modes de réalisation du procédé de l'invention aura cette même forme prédéterminée, pour se ramener au cas du contour C circulaire, comme cela a été fait dans le cas de l'ellipse. Notamment, l'équation 6 sera vérifiée et la fonction de forme utilisée sera une fonction de Bessel de première espèce Jo.

Autrement dit, la forme du contour C idéal peut être déterminé par la variation angulaire de la vitesse de phase du premier mode de flexion de la plaque. Pour une plaque isotrope, le contour C est circulaire. Pour une plaque orthotrope, le contour C est elliptique. Pour une plaque quelconque, une analyse du premier mode de flexion peut permettre de prédéterminer la forme idéale du contour C à utiliser.

Le contour C peut également prendre une forme ne correspondant pas au profil des vitesses dans le matériau de la plaque.

Par exemple, le contour C peut être un rectangle tel que présenté en figure 10, centré sur le point P, P étant le lieu où le paramètre physique est estimé. Le contour C rectangle comprend un nombre prédéterminé de points de contour Cj, j étant un indice entier positif compris entre 1 et U. U est par exemple égal à huit, de telle sorte que les huit points de contour Cj sont positionnés dans les quatre coins et au milieu de chaque côté du rectangle.

Le procédé comprend les étapes suivantes :

- on mesure en chaque point du contour C-,, un signal de contour s/t) représentatif de l'onde en chacun de ces points de contour ;

- on calcule par toute technique d'interpolation des signaux S k (t) représentatif de l'onde sur des points virtuels CI k positionnés le long d'un contour virtuel CI à l'intérieur du contour C rectangle. Notamment, le contour virtuel CI peut avoir la forme d'un cercle de rayon R prédéterminé, par exemple inférieur à la moitié du plus petit côté du rectangle ;

- on calcule la somme des signaux S k (t) représentatif de l'onde sur des points virtuels CI k pour estimer un premier signal sj(t) le long du contour virtuel CI, et on calcule une première transformée de Fourrier Si(f) dudit premier signal.

On peut alors appliquer un des procédés décrits précédemment, appliqué au contour virtuel CI de rayon R. Notamment :

- on mesure au point P, un deuxième signal S 2 (t) représentatif de l'onde au point P, et on calcule une deuxième transformée de Fourrier ^f/J dudit deuxième signal ;

- on détermine un paramètre d'échelle a de telle sorte que le module de la fonction de Bessel de première espèce Jo suivante :

où b est un autre paramètre d'échelle, et

1.1 est la fonction module,

s'approche au mieux de : \s 2 (fys!(f)\

pour un ensemble de fréquences de test /„.

On détermine alors un paramètre physique comme cela a déjà été décrit.

Comme montré en figure 11, le contour C peut avoir une forme carrée.

Selon une première variante appliquée à un contour C de forme carré, on considère que les points de contour Cj sont proches d'un cercle CI de rayon R calculé par :

di est la moitié de la longueur de la plus grande médiane du contour carré, et

d, 2 est la longueur de la diagonale du contour carré. On peut alors appliquer un des procédés décrits précédemment, appliqué au cercle CI avec le rayon calculé ci-dessus .

Selon une deuxième variante représentée en figure 12 et appliquée à un contour C de forme carré, on considère que les points de contour Cj sont inclus dans deux cercles :

- un premier cercle CIl de rayon di, passant par les points de contour Cj situés au milieu des côtés du carré, et

- un deuxième cercle CI2 de rayon d 2 , passant par les points de contour Cj situés dans les coins du carré.

L'équation (6) s'écrit maintenant sous la forme de deux équations :

W ardel (r) = 2πAJ 0 (Ud 1 )G plate (r), et

W arcle2 (r) = 2πAJ 0 (kd 2 )G plate (r) .

On compare alors par un des procédé décrit ci- dessus :

W ard J 0 ikd 2 ) + W^ 2 J 0 (Id 1 ) , et

J 0 (M 2 )J 0 (M 1 )G^r)

pour déterminer un paramètre physique de la plaque au point P. En outre, le procédé utilisant des points de contour positionnés sur un contour C en forme de rectangle, peut être avantageusement mis en oeuvre avec un vibromètre à balayage. Le vibromètre à balayage fournira les mesures vibratoires de l'onde se propageant sur la plaque 1 pour un ensemble de points, répartis sur un quadrillage matriciel de la plaque 1.

Ainsi, on peut calculer une image de la plaque représentant le paramètre physique, en utilisant successivement chaque point du quadrillage comme un point P de référence où le paramètre physique est à déterminer, et les autres points entourant immédiatement ce dernier point comme points appartenant à un contour C fermé.

L'image comporte une pluralité de pixels. Chaque pixel correspond :

- à un point de la plaque, par exemple un point mesuré par un vibromètre à balayage, et

- représente un paramètre physique de la plaque audit point de la plaque, déterminé par un des procédé décrits ci-dessus.

Notamment, il est possible par exemple de fournir une image de l'épaisseur d'une plaque, à distance et sans contact, ladite image ayant une précision spatiale égale à la distance entre les points mesurés. Une telle image permet donc de détecter, déterminer et localiser une différence d'épaisseur de la plaque.

La figure 13 représente un exemple d'une telle image réalisée sur une plaque en aluminium de 4 mm d'épaisseur ayant une zone de 100x100 mm usinée pour avoir une épaisseur de 3,5 mm. La distance entre chaque pixel de l'image est de 4 mm.

De tels produits et procédés peuvent être mis en œuvre pour mesures des épaisseurs de plaques, tôles sur de grandes structures (coque de bateau, fuselage d'avion, réservoir, bâtiments) ou de petites structures.

Ils peuvent être mis en œuvre sur des plaques planes, courbées ou des tubes.

Dans le cas de plaques courbées ou de tubes, les ondes se propagent avec des vitesses variant avec la direction de propagation. Le dispositif aura alors avantageusement un contour C fermé de forme elliptique adaptée à la courbure de la structure, comme dans le cas d'une plaque plane constituée d'un matériau anisotrope.

De tels produits et procédés ont de nombreuses applications industrielles :

- contrôle de l'épaisseur de structures, de type plaques, tôles, tubes ;

- contrôle d'une épaisseur d'un dépôt sur ces structures, tel que le dépôt de tartre dans des canalisations de transport d'eau ;

- contrôle de l'apparition de défaut dans ces structures, de 1 'endommagement ou du vieillissement de ces structures .