ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

24
ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT Décembre 2020 Projet financé par

Transcript of ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Page 1: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Décembre 2020

Projet financé par

Page 2: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Présentation ESEO 2

Page 3: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

3

Système d’imagerie ultrasonore pour le

contrôle non destructif (CND)

Présentation ESEO

Time of Flight Diffraction experiment (TOFD) Echo experiment ( A-scan)

Page 4: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

4

Système d’imagerie ultrasonore

Présentation ESEO

Page 5: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

5

Système d’imagerie ultrasonore non-linéaire

Présentation ESEO

Page 6: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

6

Système d’imagerie ultrasonore médical

Présentation ESEO

Page 7: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

7

Système d’imagerie ultrasonore médical

Présentation ESEO

Page 8: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

8

Système d’imagerie ultrasonore non-linéaire en CND

Présentation ESEO

Page 9: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

9

Système d’imagerie ultrasonore en CND

Présentation ESEO

Quelle est la commande qui maximise la probabilité de détection des défauts ?

Page 10: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

10

Système d’imagerie ultrasonore non-linéaire

Présentation ESEO

Page 11: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Présentation ESEO 11

Page 12: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Principe de la commande optimale

Présentation ESEO 12

Page 13: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

1. Choix de la fonction de coût J(θ)Objectif : maximiser la probabilité de détection

Signature nonlinéaire des défauts

Contraintes

Extraction des nonlinéairités

en sortie du système :

• Filtrage du signal de sortie

• Codage

• Modélisation (Hammerstein, Volterra…)

13

Mise en œuvre de la commande optimale

Présentation ESEO

Syst.input 𝑋𝐿𝑖𝑛 + 𝑋𝑁𝐿Parameters

(amplitude, frequency, shape ...)

Nonlinear part (high amplitude responce)

Linear part (lowamplitude responce)

ൠmax(𝐸𝑁𝐿) → +∞min(𝐸𝐿𝑖𝑛) → 0

max 𝐽 = 𝑚𝑎𝑥𝐸𝑁𝐿𝐸𝐿𝑖𝑛

Page 14: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

1. Choix de la fonction de coût J(θ)Objectif : maximiser la probabilité de détection

Signature nonlinéaire des défauts

2. Choix des paramètres θ d’entréeSimplicité

Hypothèse physique : résonance…

θ : Fréquence

3. Choix de l’algorithme d’optimisationSi fonction J(θ) convexe algorithme d’optimisation locale

Si fonction J(θ) avec plusieurs optimum locaux métaheuristique

algorithme d’optimisation locale : Simplex downhill

14

Mise en œuvre de la commande optimale

Présentation ESEO

Syst.input 𝑋𝐿𝑖𝑛 + 𝑋𝑁𝐿Parameters

(amplitude, frequency, shape ...)

Nonlinear part (high amplitude responce)

Linear part (lowamplitude responce)

ൠmin(𝐸𝑁𝐿) → +∞max(𝐸𝐿𝑖𝑛) → 0

max 𝐽 = 𝑚𝑎𝑥𝐸𝑁𝐿𝐸𝐿𝑖𝑛

Page 15: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Présentation ESEO 15

Page 16: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Simulation : propagation et crack

Présentation ESEO 16

2D aluminium plate

P-wavespeed [m/s]

6398

S-wavespeed [m/s]

3030

Density[kg/m3]

3122

Sample

dimensions

[100 ; 50]

mm

Crack

dimensions[8 ; 8] mm

• Propagation : Pseudospectral / staggered A-B method, with linear propagation of longitudinal and shear waves.

• Crack hysteresis behavior described by hysteretic elementaryunits (HEU)

Young

modulus

[Pa/m2]

70·109

Poisson

ratio 0.33

Density

[kg/m3] 2966

Str

ain

Stress

𝑃𝑜 ≠ 𝑃𝑐

𝑃𝑜 𝑃𝑐

𝑃𝑜 = 𝑃𝑐

𝑃𝑜 = 𝑃𝑐

Str

ain

StressOpening

Pressure

𝑃𝑜

Closing

Pressure𝑃𝑐

𝜌 𝑃𝑜, 𝑃𝑐

McCall & Guyer, J. Geophys. Res. 99 (1994) 23887–23897

Page 17: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Simulation : capteurs et signaux

Présentation ESEO 17

2D aluminium plate Tx Band-pass filter [fc = 1MHz]

Rx Band-pass filter [fc = 3MHz]

𝑥𝑒𝑥𝑡 𝑡 = 𝐴 ∙ 𝑒𝑥𝑝(−𝑡 − 𝑡0

2

𝑡02

) ∙ sin(𝑡 ∙ 2𝜋 ∙ 𝑓𝑒𝑥𝑡)

Output signal

Page 18: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Système d’imagerie ultrasonore non-linéaire en CND

Présentation ESEO 18

Cas linéaire

Cas non-linéaire

Crack

Page 19: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Simulation : optimisation de la fréquence

Présentation ESEO 19

Frequency empirical sweep Mechanical stress

700:25:1100 kHz 800 kPa (1 kPa = 1000 N/m 2)

Exemple of optimisation for 800 kPa

Page 20: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Simulation : optimisation de la fréquence

Présentation ESEO 20

Mechanical stress

100 : 100: 1500 kPa

The optimal nonlinear ratio and optimal frequency for different stresses

The gain of nonlinear linear ratio between optimal frequencies and central frequency of the transducer

Page 21: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

• Mise en place du démonstrateur

21

Experimentation

Présentation ESEO

Page 22: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Présentation ESEO 22

Page 23: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

• Adaptation de la commande optimale au CND

• Cas paramétrique fréquence

• Application en simulation

• Présentation au Forum Acusticum, Lyon, décembre 2020

• Mise en œuvre du démonstrateur

• En cours

• Problématique : plusieurs sources de non-linéarités

Développement de méthode d’extraction des non-linéarités à

partir de modèle de Hammerstein

Conclusion

Présentation ESEO 23

Page 24: ÉNIGOT, Volodymyr GATSA & Vincent TOURNAT

Projet financé par

[email protected]