Peronnet

%0$503"5%&-6/*7&34*5%&506-064&
Institut National Polytechnique

de Toulouse (INP Toulouse)

Gnie mcanique, mcanique
des matriaux

Elodie PERONNET

jeudi 4 octobre 2012

Identification exprimentale du comportement d'un fuselage composite :
mesures de champs
dtection de dfauts par

Mcanique, Energtique, Gnie civil et Procds (MEGeP)

Laboratoire Gnie de Production (Ecole Nationale d'Ingnieurs de Tarbes)

Moussa KARAMA Professeur des Universits, Ecole Nationale d'Ingnieurs de Tarbes
Sbastien MISTOU Matre de Confrences, Ecole Nationale d'Ingnieurs de Tarbes

Alain GERARD Professeur
des Universits, Universit de Bordeaux
Abdelghani SAOUAB Professeur
des Universits, Universit du Havre
M :
Florent EYMA Matre de Confrences, IUT de Tarbes
Abdelouhed LAKSIMI Professeur des Universits, UT de Compigne
Mostapha TARFAOUI Matre de confrences - HDR, ENSTA Bretagne
Hlne WELEMANE Matre de Confrences, Ecole Nationale d'Ingnieurs de Tarbes
Serge CRISTANTE Ingnieur DAHER SOCATA SAS Invit
Remerciements
Ces travaux de thse ont t raliss au sein de lquipe Mcanique des Matriaux, des
Structures et des Procds (M2SP) du Laboratoire Gnie de Production (LGP) de lcole
Nationale dIngnieurs de Tarbes (ENIT). Ils ont t financs, dans le cadre du programme
FUSCOMP, par DAHER-SOCATA et par le Conseil Gnral des Hautes Pyrnes.
Je tiens tout dabord remercier lensemble des membres du jury qui ont aliment les
discussions grce leurs remarques et questions constructives, et ont permis denrichir ce
manuscrit. Je remercie plus particulirement Mostapha Tarfaoui davoir prsid ce jury,
ainsi que Alain Grard et Abdelghani Saouab davoir accept de rapporter mon mmoire
de thse.
Je remercie galement Moussa Karama, mon directeur de thse, de sa disponibilit

lorsque ncessaire. Jadresse de remerciements Sbastien Mistou, mon co-directeur de
thse, pour mavoir soutenu tout au long de ces travaux, dans les bons comme dans les
moments difficiles et de mavoir cadr et recadr quand ctait ncessaire. Merci aussi
Hlne Welemane et Florent Eyma de stre impliqus dans la ralisation de ces travaux
de thse.
Je veux galement remercier Marie-Laetitia Pastor et Olivier Dalverny qui ont eu un

rle cl dans le dveloppement des deux thmatiques de cette thse.
Je remercie Serge Cristante de stre toujours souci de la valorisation industrielle de

cette thse.
Jaimerais aussi remercier lensemble des membres du CRCC et de Nimitech qui mont
accueilli avec bienveillance et aid la fabrication des cinq prouvettes talons. Je vou-
drais remercier, plus particulirement, Nathalie Coureau pour son implication dans la
ralisation de lensemble des plaques talons.
Jadresse aussi des remerciements Christophe Pescay et Marianne Perrin pour les
nombreuses mesures ultrasonores ralises. Je remercie galement, Richard Huillery de
Thermoconcept et Sylvain Gnot de TomoAdour pour leur collaboration.
Je passe ensuite une ddicace spciale tous les doctorants que jai eu le plaisir
de ctoyer durant ces trois ans Tarbes et toute la troupe du Laboratoire Gnie de
Production. Et je remercie, plus particulirement, mes colloques de bureau (Benji et
Damien) ainsi que mes voisins de palier (Valentin et Marina) pour mavoir support et
soutenu.
Table des matires
Table des matires i
Table des figures v
Liste des tableaux xi
Introduction 1
1 Contrle de dfauts par mesures de champs acoustiques, densimtriques,

thermiques 5
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2 Mesures de champs acoustiques : les ultrasons . . . . . . . . . . . . . . . . 8
1.2.1 Principe physique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2.2 Principe de contrle . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.2.3 Bilan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.3 Mesures de champs densimtriques : radiographie rayons X et tomogra-
phie rayons X . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.3.2 Principe de contrle . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.3.3 Bilan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.4 Mesures de champs thermiques : la thermographie infrarouge . . . . . . . . 29
1.4.2 Principe du contrle . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.4.3 Bilan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
1.5 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2 Matriaux et protocoles exprimentaux 41

2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.2 Prsentation des chantillons . . . . . . . . . . . . . . . . . . . . . . . . . . 42
ii TABLE DES MATIRES
2.2.1 Eprouvette 1 : matriau monolithique avec inserts Tflon . . . . . 43

2.2.2 Eprouvette 2 : matriau monolithique avec trous fond plat . . . . 44
2.2.3 Eprouvette 3 : matriau sandwich avec inserts Tflon (peaux de
0.66 mm) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.2.4 Eprouvette 4 : matriau sandwich avec inserts Tflon (peaux de
1.08 mm) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.2.5 Eprouvette 5 : matriau multi-structures . . . . . . . . . . . . . . . 50
2.3 Protocoles exprimentaux . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.3.1 Mesures de champs acoustiques : les ultrasons . . . . . . . . . . . . 54
2.3.2 Mesures de champs densimtriques . . . . . . . . . . . . . . . . . . 61
2.3.3 Mesures de champs thermiques : thermographie infrarouge . . . . . 66
2.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
3 Etude comparative sur prouvettes talons 75

3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
3.2 Principes danalyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
3.2.1 Dimensionnement des dfauts par traitement dimage . . . . . . . . 77
3.2.2 Critres danalyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
3.3 Rsultats exprimentaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
3.3.1 Etalons monolithiques avec inserts Tflon . . . . . . . . . . . . . 80
3.3.2 Etalon monolithique avec trous fond plat . . . . . . . . . . . . . . 88
3.3.3 Etalons sandwich avec inserts Tflon . . . . . . . . . . . . . . . . 94
3.3.4 Synthse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
3.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
4 Identification des dfauts par mesures de champs cinmatiques 111

4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
4.2 Mthodes de mesures de champs . . . . . . . . . . . . . . . . . . . . . . . . 113
4.2.1 Mesures surfaciques (2D) . . . . . . . . . . . . . . . . . . . . . . . . 113
4.2.2 Mesures surfaciques et hors-plan (2D 1/2) . . . . . . . . . . . . . . 115
4.2.3 Mesures volumiques (3D) . . . . . . . . . . . . . . . . . . . . . . . . 118
4.3 Mthodes didentification du comportement . . . . . . . . . . . . . . . . . 119
4.3.1 Mthodes itratives . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
4.3.2 Mthodes non itratives . . . . . . . . . . . . . . . . . . . . . . . . 123
4.4 Mthodologie didentification multi-zones . . . . . . . . . . . . . . . . . . . 125
4.5 Identification des proprits dun matriau orthotrope plan . . . . . . . . . 128
TABLE DES MATIRES iii
4.5.1 Modles de rfrence et de recalage . . . . . . . . . . . . . . . . . . 129

4.5.2 Description du fonctionnement de lalgorithme didentification . . . 131
4.5.3 Rsultats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
4.6 Identification des proprits dune zone dgrade . . . . . . . . . . . . . . . 135
4.6.1 Cration de la cartographie derreurs . . . . . . . . . . . . . . . . . 136
4.6.2 Estimation de la taille de la zone dgrade . . . . . . . . . . . . . . 139
4.6.3 Identification des proprits de la zone dgrade . . . . . . . . . . . 141
4.7 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144
Conclusion 147
iv TABLE DES MATIRES
Table des figures
1.1 Principe gnral dun contrle non destructif . . . . . . . . . . . . . . . . 6

1.2 Intensit acoustique sur laxe du faisceau (D est le diamtre de la source
utilise et la longueur donde) (Dumont-Fillon [71]) . . . . . . . . . . . . 9
1.3 Types dondes et modes de propagation des ultrasons dans un solide (Dumont-
Fillon [71]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4 Ondes de volume : (a) ondes de compression, (b) ondes de cisaillement
(Lefebvre et al. [133]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.5 Ondes et contraintes induites : (a) onde volumique de compression, (b)
onde volumique de cisaillement, (c) onde de surface (Mistou [143]) . . . . . 11
1.6 Mthode cho ultrasonore par immersion ou par contact avec sabot (Cher-
faoui [55], Lambert et al. [121]) . . . . . . . . . . . . . . . . . . . . . . . . 14
1.7 Mthode cho ultrasonore par contact (Framezelle et al. [79]) . . . . . . . 14
1.8 Diffrents modes de prsentation de linspection ultrasonore (Dumont-Fillon
[71]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.9 Contrle en mode transmission (Chaix et al. [53]) . . . . . . . . . . . . . . 16
1.10 Principe de dtection par transmission (Lacroix [119]) . . . . . . . . . . . . 17
1.11 Spectre lectromagntique (Garnier [83]) . . . . . . . . . . . . . . . . . . . 20
1.12 Interaction entre le rayonnement et la matire lors de la traverse dun
matriau (Lecomte [130]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
1.13 Evolution du coefficient dattnuation dun milieu en fonction de lnergie
E des photons . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
1.14 Principe de contrle par rayons X (Lecomte [130]) . . . . . . . . . . . . . 22
1.15 Contraste dun dfaut sur limage radiante (Lecomte [130]) . . . . . . . . . 23
1.16 Schma dun tube rayons X (Badel [24], Le-Bras [125], Lecomte [130]) . . 24
1.17 Schma de principe du contrle par tomographie X : du faisceau jusquaux
coupes (Carmona [50]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
1.18 Schma dune coupe (Lenoir [134]) . . . . . . . . . . . . . . . . . . . . . . 25
vi TABLE DES FIGURES
1.19 Gnrations de scanners (Lenoir [134]) . . . . . . . . . . . . . . . . . . . . 26

1.20 Principe de la mthode de reconstruction dimages tomographiques (Lenoir
[134]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.21 Bilan radiatif (Pajani [155]) . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.22 Types de sollicitations thermiques (Maillard et al. [135]) . . . . . . . . . . 33
1.23 Correspondance entre les mthodes de thermographie IR et leur signal res-
pectif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
1.24 Thermographie infrarouge en face avant . . . . . . . . . . . . . . . . . . . 34
1.25 Thermographie infrarouge en face arrire . . . . . . . . . . . . . . . . . . . 35
1.26 Schma de cblage de la TIR . . . . . . . . . . . . . . . . . . . . . . . . . 36
1.27 Schma de cblage de la TIRM . . . . . . . . . . . . . . . . . . . . . . . . 36
1.28 TIRM ultrasonore (Gleiter et al. [89]) . . . . . . . . . . . . . . . . . . . . . 38
1.29 Schma de cblage de la TIRP . . . . . . . . . . . . . . . . . . . . . . . . 38
2.1 Photo du matriau monolithique avec inserts Tflon . . . . . . . . . . . 43

2.2 Vue de dessus du matriau monolithique avec inserts Tflon - Position-
nement des inserts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
2.3 Matriau monolithique avec inserts Tflon - Coupe schmatique dans
lpaisseur au niveau dune range dinserts de mme diamtre . . . . . . . 44
2.4 Photo du matriau monolithique avec trous fond plat . . . . . . . . . . . 45
2.5 Vue de dessus du matriau monolithique avec trous fond plat - Position-
nement des trous . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2.6 Matriau monolithique avec trous fond plat - Coupe schmatique dans
lpaisseur au niveau dune range de trous de mme diamtre . . . . . . . 46
2.7 Photo du matriau sandwich avec inserts Tflon (peaux de 0.66 mm) . . 46
2.8 Vue de dessus du matriau sandwich avec inserts Tflon (peaux de 0.66
mm) - Positionnement des inserts Tflon . . . . . . . . . . . . . . . . . . 47
2.9 Matriau sandwich avec inserts Tflon (peaux de 0.66 mm) - Coupe sch-
matique au niveau dune range dinserts de mme diamtre . . . . . . . . 47
2.10 Photo du matriau sandwich avec inserts Tflon (peaux de 1.08 mm) . . 48
2.11 Vue de dessus du matriau sandwich avec inserts Tflon (peaux de 1.08
mm) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.12 Matriau sandwich avec inserts Tflon (peaux de 1.08 mm) - Coupe sch-
matique dans lpaisseur au niveau dune range dinserts de mme diamtre 49
2.13 Photo du matriau multi-structures . . . . . . . . . . . . . . . . . . . . . 50
2.14 Master en rsine usinable . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
TABLE DES FIGURES vii
2.15 Moule rigide en verre-poxy . . . . . . . . . . . . . . . . . . . . . . . . . . 51

2.16 Vue de dessus du matriau multi-structures - Positionnement des inserts
Tflon . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
2.17 Matriau multi-structures - Coupes schmatiques dans lpaisseur au ni-
veau dune range dinserts de mme diamtre : (a) partie monolithique
fine, (b) zone sandwich, (c) partie paisse . . . . . . . . . . . . . . . . . . . 53
2.18 Mise en place du film micro perfor . . . . . . . . . . . . . . . . . . . . . . 53
2.19 Capteur multi-lments quip du sabot . . . . . . . . . . . . . . . . . . . 55
2.20 Dispositif exprimental du contrle par cho ultrasonore par contact . . . 55
2.21 Systme de fixation de la structure inspecte (matriau multi-structures) . 56
2.22 A-scan ralis en cho ultrasonore par contact (phase dtalonnage, sans
dfaut) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
2.23 C-Scan (matriau monolithique avec insert Tflon ) : (a) en amplitude,
(b) en paisseur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
2.24 Pice et capteur immergs dans la cuve . . . . . . . . . . . . . . . . . . . 58
2.25 Dispositif exprimental de lcho ultrasonore par immersion . . . . . . . . 59
2.26 Dispositif de fixation de lprouvette . . . . . . . . . . . . . . . . . . . . . 59
2.27 Montage de transmission ultrasonore par immersion . . . . . . . . . . . . 60
2.28 Positionnement de la porte A . . . . . . . . . . . . . . . . . . . . . . . . . 61
2.29 C-Scan en amplitude (sur une partie du matriau multi-structures) . . . . 61
2.30 Montage exprimental du contrle par radiographie X . . . . . . . . . . . 61
2.31 Dispositif exprimental ddi au contrle par radiographie X utilis par
Daher Socata sur le site de Tarbes . . . . . . . . . . . . . . . . . . . . . . 62
2.32 Radiographie X (matriau sandwich avec peau de 0.66 mm) . . . . . . . . 63
2.33 Scanner SOMATOM Sensation Open . . . . . . . . . . . . . . . . . . . . . 64
2.34 Schma de fonctionnement du scanner mdical SIEMENS . . . . . . . . . 64
2.35 Reconstructions par tomographie X (matriau sandwich avec peaux de 0.66
mm) : (a) Coupe 2D reconstruite paralllement la surface moyenne de la
pice, (b) : reconstruction volumique . . . . . . . . . . . . . . . . . . . . . 65
2.36 Dispositif exprimental en configuration face avant . . . . . . . . . . . . . 67
2.37 Dispositif exprimental en configuration face arrire . . . . . . . . . . . . 67
2.38 Mesures de champ thermique (en relatif) ralises par TIR (matriau mo-
nolithique avec inserts Tflon ) . . . . . . . . . . . . . . . . . . . . . . . . 69
2.39 Cartographie 2D des mesures de champs thermiques obtenus par TIRM
(matriau monolithique avec inserts Tflon ) : (a) amplitude, (b) phase . 70
2.40 Dispositif exprimental de la TIRP en configuration face avant . . . . . . 71
viii TABLE DES FIGURES
2.41 Dispositif exprimental de la TIRP en configuration face arrire . . . . . . 72

2.42 Cartographie 2D des mesures de champs thermiques par TIRP (matriau
monolithique avec inserts Tflon ) . . . . . . . . . . . . . . . . . . . . . . 73
3.1 Procdure de traitement dimage automatique . . . . . . . . . . . . . . . . 78

3.2 Prcision de dtection en fonction de la profondeur des dfauts - Etalons
monolithiques : (a) matriau monolithique avec inserts Tflon , (b) prou-
vette multi-structure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
3.3 Prcision de dtection en fonction du diamtre des dfauts - Etalons mono-
lithiques : (a) matriau monolithique avec inserts Tflon , (b) prouvette
multi-structure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
3.4 Prcision de dtection en fonction de la profondeur des dfauts - Matriau
monolithique avec trous fond plat . . . . . . . . . . . . . . . . . . . . . 91
3.5 Prcision de dtection en fonction du diamtre de dfauts - Matriau mo-
nolithique avec trous fond plat . . . . . . . . . . . . . . . . . . . . . . . 92
sandwich avec peaux de 0.66 mm . . . . . . . . . . . . . . . . . . . . . . . 99
sandwich avec peaux de 1.08 mm . . . . . . . . . . . . . . . . . . . . . . . 99
3.8 Prcision de dtection en fonction de la profondeur des dfauts - Zones
sandwichs de lprouvette multi-structure . . . . . . . . . . . . . . . . . . 100
3.9 Prcision de dtection en fonction du diamtre de dfauts - Matriau sand-
wich avec peaux de 0.66 mm . . . . . . . . . . . . . . . . . . . . . . . . . 100
3.10 Prcision de dtection en fonction du diamtre de dfauts - Matriau sand-
wich avec peaux de 1.08 mm . . . . . . . . . . . . . . . . . . . . . . . . . 101
3.11 Prcision de dtection en fonction du diamtre de dfauts - Zones sandwichs
de lprouvette multi-structure . . . . . . . . . . . . . . . . . . . . . . . . 101
3.12 Prcision moyenne de dtection - Structures monolithiques . . . . . . . . . 104
3.13 Prcision moyenne de dtection - Structures sandwichs . . . . . . . . . . . 106
3.14 Prcision moyenne de dtection - Dfauts de 6 mm de diamtre . . . . . . 107
4.1 Sytme possdant deux camras (Fazzini [76]) . . . . . . . . . . . . . . . . 116

4.2 Principe stroscopique (Fazzini [76]) . . . . . . . . . . . . . . . . . . . . . 116
4.3 Images de gauche et droite et reconstruction 3D du cylindre (Fazzini [76]) 117
4.4 Dtermination du champ de dplacement 3D par stro-corrlation (Fazzini
[76]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
4.5 Principe de la mthode didentification FEMU . . . . . . . . . . . . . . . 121
TABLE DES FIGURES ix
4.6 Problme considr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126

4.7 Stratgie didentification multi-zones . . . . . . . . . . . . . . . . . . . . . 126
4.8 Structure multi-zones de rfrence au niveau de la zone dintrt : (a) prin-
cipe, (b) modle exprimental numrique associ (vue de la face arrire,
les zooms montrent les diffrents plis de chaque zone) . . . . . . . . . . . . 127
4.9 Orientation de la symtrie matrielle du matriau orthotrope . . . . . . . . 129
4.10 Identification des proprits orthotropes - Zone dintrt . . . . . . . . . . 130
4.11 Schma de fonctionnement de lalgorithme didentification . . . . . . . . . 131
4.12 Courbes de convergence du cas FUSCOMP . . . . . . . . . . . . . . . . . . 135
4.13 Calcul derreur - Zone dintrt . . . . . . . . . . . . . . . . . . . . . . . . 136
4.14 Evolution des composantes (Ux , Uy ) du champ de dplacement U et des
composantes (xx , yy , xy ) du champ de dformations dans la base (x, y)
en face avant de la structure de rfrence le long dune ligne AB traversant
les zones saine et dgrade . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
4.15 Cartographie derreur en dformation (rfrence vs matriau monozone sain)139
4.16 Identification de la taille de la zone dgrade : (a) coupes ralises sur la
carte derreur, (b) et (c) : volution des niveaux de gris et dfinition des
seuils . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
4.17 Intgration de la zone dgrade . . . . . . . . . . . . . . . . . . . . . . . . 141
4.18 Identification des proprits dgrades - Zone dintrt . . . . . . . . . . . 141
4.19 Cartographie derreur en dformation (rfrence vs structure multizones
base sur la thorie de stratifis) . . . . . . . . . . . . . . . . . . . . . . . . 142
4.20 Courbes de convergence de lidentification multi-zone . . . . . . . . . . . . 144
x TABLE DES FIGURES
Liste des tableaux
3.1 Rsultats du traitement dimage par famille de mesures de champs . . . . 79

3.2 Sensibilit de dtection - Matriau monolithique avec inserts Tflon . . . 81
3.3 Sensibilit de dtection - Zones monolithiques de lprouvette multi-structure
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
3.4 Prcision moyenne de dtection - Etalons monolithiques . . . . . . . . . . 84
3.5 Sensibilit de dtection - Matriau monolithique avec trous fond plat . . 89
3.6 Prcision moyenne de dtection - Matriau monolithique avec trous fond
plat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
3.7 Sensibilit de dtection - Matriau sandwich avec peaux de 0.66 mm . . . 95
3.8 Sensibilit de dtection - Matriau sandwich avec peaux de 1.08 mm . . . 95
3.9 Sensibilit de dtection - Zones sandwichs de lprouvette multi-structure . 96
3.10 Prcision moyenne de dtection - Etalons sandwichs . . . . . . . . . . . . 98
3.11 Sensibilit de dtection - Structures monolithiques . . . . . . . . . . . . . . 104
3.12 Sensibilit de dtection - Structures sandwichs . . . . . . . . . . . . . . . 105
3.13 Sensibilit de dtection - Dfauts de 6 mm de diamtre . . . . . . . . . . . 107
4.1 Donnes de base et rsultats des problmes directs et indirects (Kleiner-

mann [118]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
4.2 Proprits planes du tissu Hexcel HexForce G0926 dans ses axes principaux130
4.3 Identification de matriaux orthotropes par procdure FEMU (le calcul de
chaque paramtre selon la thorie classique des stratifis est donn entre
parenthses) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
4.4 Thorie des stratifis : empilement [45/0/45] . . . . . . . . . . . . . . . 142
xii LISTE DES TABLEAUX
Introduction
Le dveloppement de nouveaux procds de fabrication pour les composites constitue

lheure actuelle une problmatique cruciale et ouverte, dont les enjeux concernent la fois
loptimisation des mthodes et tapes de mise en oeuvre mais galement celle des mat-
riaux et structures du point de vue de leurs performances mcaniques. Dans ce contexte, le
procd dinfusion de rsine liquide (Liquid Resin Infusion) prsente de nombreux atouts
dont la possibilit de fabriquer des pices de formes complexes et des panneaux entiers de
fuselage en matriaux composites, rduisant considrablement les tapes dassemblages et
donc les temps de production. Paralllement, la rsistance mcanique de ces pices savre
particulirement intressante grce notamment la possibilit de raliser des structures
intgres (raidisseurs par exemple).
Labellis par le ple de comptitivit Aerospace Valley, le projet de Recherche & Dve-
loppement FUSCOMP (FUSelage COMPosite) vise la ralisation dun test de dmons-
trateur de fuselage composite produit par LRI. Ces travaux runissent principalement
quatre partenaires de la Rgion Midi-Pyrnes : la socit DAHER-SOCATA situe
Louey, le Laboratoire Gnie de Production de lEcole Nationale dIngnieurs de Tarbes,
lInstitut Universitaire Technologique de Tarbes (dpendant de lInstitut Clment Ader
de Toulouse) et le Centre de Ressources et de Comptences Composites (CRCC) du Lyce
Jean Dupuy de Tarbes.
Le thme principal de la thse prsente ici est lidentification de dfauts internes au

sein de matriaux composites constituant le fuselage. Il sagit dun sujet primordial pour
le dveloppement dun procd de fabrication puisque la validation de la sant matire
des pices est porte par de nombreuses normes aronautiques, dont les standards AITM
(Airbus Industry Test Method - Mthode dessai dAIRBUS Industrie), et correspond
une tape obligatoire du processus qualit permettant de rpondre aux exigences du client
avant livraison.
2 Introduction
Dans ce cadre, deux objectifs majeurs ont t viss pour ce travail de recherche. Le pre-
mier porte sur lvaluation des mthodes de Contrle Non Destructif (CND) sur ces types
de matriaux, incluant des stratifis monolithiques et des sandwichs. DAHER-SOCATA
ayant la volont dacqurir de nouvelles comptences dans ce domaine, cette premire
tape a consist en une analyse comparative de diffrentes techniques juges intressantes
pour les applications industrielles. Prcisment, les mthodes tudies sont toutes bases
sur des mesures de champs capables de mesures globales sous forme de cartographies, ce
qui est un avantage certain du point de vue de la maintenance des structures. Ces tech-
niques se distinguent par ailleurs par leur rapidit dexcution, leur capacit contrler
des pices gomtrie complexe et leur utilisation le plus souvent sans contact.
Dans la continuit de cette analyse, le second enjeu abord dans ce travail porte
sur lidentification paramtrique des matriaux composites plans partir de mesures de
champs cinmatiques. Lide ce stade est denrichir la caractrisation fournie par le
CND au travers dune quantification des proprits lastiques. Sappuyant sur certains
dveloppements existants au sein du laboratoire, une procdure didentification globale
puis locale par recalage de modles lments finis a t mise en place de manire pouvoir
capter le caractre htrogne dune structure prsentant un dfaut et dterminer les
paramtres de chacune des zones (saine et dgrade).
La rdaction de ce mmoire sorganise autour de quatre parties.
Le premier chapitre est ddi la revue bibliographique des mthodes de CND exis-
tantes et bases sur des mesures de champs. Trois grandes familles sont notamment trai-
tes, fournissant des mesures de champs acoustiques, densimtriques et thermiques sur
structures composites. Pour chacune des mthodes abordes dans lanalyse comparative,
des prsentations du principe physique et du principe de contrle sont ralises. Un bilan
des avantages et inconvnients de chaque famille est galement propos.
La prsentation des matriaux et des protocoles exprimentaux employs fait lobjet du

deuxime chapitre. Outre les dtails sur la constitution et la fabrication des prouvettes
employes ainsi que sur leurs dfauts, cette partie expose prcisment la dmarche des
mesures ralises de manire mieux comprendre les rsultats obtenus par la suite.
Le troisime chapitre est consacr ltude comparative en elle-mme. A partir de

critres danalyse dfinis en lien avec les objectifs industriels, les diffrentes mthodes de
contrle utilises sont values sur les diffrents matriaux et les diffrents types de dfauts.
Introduction 3
Une synthse gnrale par type de structures vient conclure cette partie en indiquant les
tendances observes pour les trois familles de mesures de champs.
Lidentification des dfauts par mesures de champs cinmatiques est aborde dans
le dernier chapitre. Aprs une synthse bibliographique sur les mthodes de mesures de
champs et les techniques didentification les plus rpandues, la procdure didentification
sur structure composite multi-zones est dtaille de manire apprhender les verrous
scientifiques levs dans le cas prsent (matriaux anisotropes et macroscopiquement ht-
rognes). Un certain nombre de rsultats sont enfin prsents pour illustrer chacune des
diffrentes tapes de la mthodologie.
4 Introduction
Chapitre 1
Contrle de dfauts par mesures de

champs acoustiques, densimtriques,
thermiques
Sommaire
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2 Mesures de champs acoustiques : les ultrasons . . . . . . . . . 8
1.2.1 Principe physique . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2.2 Principe de contrle . . . . . . . . . . . . . . . . . . . . . . . . 13
1.2.3 Bilan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.3 Mesures de champs densimtriques : radiographie rayons
X et tomographie rayons X . . . . . . . . . . . . . . . . . . . 19
1.3.2 Principe de contrle . . . . . . . . . . . . . . . . . . . . . . . . 22
1.3.3 Bilan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.4 Mesures de champs thermiques : la thermographie infrarouge 29
1.4.2 Principe du contrle . . . . . . . . . . . . . . . . . . . . . . . . 32
1.4.3 Bilan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
1.5 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
1. Contrle de dfauts par mesures de champs acoustiques,
6 densimtriques, thermiques
1.1 Introduction
Les Contrles Non Destructifs (CND) ont pour vocation de fournir des informations
concernant la sant matire des pices contrles sans altrer ni modifier leurs proprits
physico-chimiques ainsi que mcaniques. Ils valuent les dfauts dune pice durant son
cycle de vie (durant le processus de fabrication, aprs fabrication, en maintenance) afin
de dfinir si celle-ci est toujours utilisable suivant les critres mcaniques associs lors de
sa conception. Cette dmarche dinvestigation est un lment majeur du contrle qualit
des pices dans lindustrie aronautique.
Gnralement, une technique de CND se dcompose en trois phases (Boro-Djordjevic

[44], Cawley [52], Cherfaoui [55], Framezelle et al. [79], Hsu [107], figure 1.1) :
la sollicitation, souvent de nature mcanique, lectromagntique ou thermique,
la perturbation, dfinie par la pice et les anomalies quelle contient, et qui traduit
linteraction avec la sollicitation,
la rvlation, assure soit par lil soit par des capteurs associs un systme
imageur, traduisant la raction entre lexcitation et la perturbation.
Figure 1.1 Principe gnral dun contrle non destructif
En ce qui concerne les structures composites, les principaux dfauts que lon cherche
identifier sont dune part les imperfections dues au procd de fabrication (porosits,
mauvaises imprgnations, fibres non alignes,..) et dautre part, les dcohsions de matire
1.1. Introduction 7
induits en service (rupture de fibres, fissures, dlaminages, dcollements fibres-matrice,

impacts,..) (Adams and Cawley [1], Colin [62], Geier [85]).
Parmi les techniques de CND destines rpondre cette problmatique, on distingue
gnralement deux grandes classes :
les contrles surfaciques tels que le contrle visuel (Cherfaoui [55], Perdijon [162]),
les courants de Foucault (Gros [101], Mook et al. [149]) ou le ressuage (Al-Maghribi
[14], Duplessis-Kergomard [72], Geier [85]) ; ces mthodes visent dtecter des in-
dications ou des dfauts dbouchants mais leur application au cas des composites
reste limite dans lindustrie du fait en particulier de perturbations lies la couleur
des matriaux concerns ou la contamination par le fluide dexamen/de ressuage
(Tsao and Hocheng [192]) ;
les contrles volumiques savrent quant eux plus performants car ils permettent
daccder galement des dfauts internes en prcisant ventuellement leur profon-
deur ; outre la shearographie (Hung [108], Hung et al. [109], Muzet et al. [152]) ou
lmission acoustique (Aggelis et al. [3], Benmedakhene et al. [35], Dzenis and Qian
[74], Mizutani et al. [144], Wevers [195]), on trouve notamment dans cette catgo-
rie diffrentes mthodes bases sur des mesures de champs tels que les ultrasons,
les analyses par rayons X ou la thermographie infrarouge. Labsence de contact, la
richesse des informations obtenues (mesures globales), la rapidit dexcution ainsi
que ladaptation linspection des pices gomtries complexes sont autant dar-
guments pour expliquer leur plein essor (Cawley [52], Grediac [93]).
Dans cette premire partie du travail, nous avons focalis notre tude sur trois types
dacquisition :
la mesure de champs acoustiques (contrle ultrasonore), qui est actuellement la
technique la plus employe dans le secteur aronautique,
la mesure de champs densimtriques, incluant la radiographie rayons X et la
tomographie rayons X,
la mesure de champs thermiques par thermographie infrarouge.
Toutes ces mesures sont bases sur le mme principe gnral :

1. lmission par une source dun flux dondes dont la nature diffre suivant la mesure
concerne,
2. la rception par un dtecteur de ce flux propag par la pice analyse ; deux op-
tions sont ici possibles : soit une mesure en rflexion (on dit aussi contrle face
avant) si le dtecteur se trouve du mme ct que lmetteur 1 , soit une mesure en

transmission (ou contrle face arrire) si le dtecteur est situ du ct oppos
lmetteur,
3. lobtention dune mesure de champ par cartographie 2D ou 3D selon les techniques.
On se propose donc prsent de dcrire les mthodes prcites en prcisant pour cha-
cune leur principe physique ainsi que leurs diffrentes configurations dutilisation dans
lindustrie.
1.2 Mesures de champs acoustiques : les ultrasons

La mesure de champ acoustique tudie ici est celle du contrle ultrasonore des pices
qui correspond la mthodologie danalyse non destructive la plus rpandue actuellement
dans lindustrie aronautique (Adams and Cawley [1], Adams and Drinkwater [2]). Par
le recul exprimental et industriel cumul depuis des annes, la pratique de cette tech-
nique est trs normalise (AITM6-0013 [5], AITM6-0016 [7], AITM6-4002 [8], AITM6-
4010 [10], AITM6-4014 [11], AITM6-0011 [4], AITM6-4005 [9]).
Suivant la dmarche gnrale expose prcdemment, le contrle ultrasonore suit les

tapes suivantes (Dumont-Fillon [71], Perdijon [161]) :
lmission dun flux dondes acoustiques (volumiques ou superficielles),
la rception dun flux acoustique par un dtecteur, en face avant ou arrire,
lobtention dune mesure de champ acoustique en cartographie 2D (qui permet
didentifier la projection en surface des dfauts), dont on peut tirer ventuellement
une reprsentation 3D (par combinaison des vues 2D).
1.2.1 Principe physique

1.2.1.1 Le positionnement optimal du traducteur
Le traducteur ultrasonore est le capteur qui permet dmettre et/ou de recevoir les
ondes ultrasonores. Il est principalement constitu par le transducteur, soit llment
pizo-lectrique traduisant les pressions mcaniques induites par les ondes ultrasonores
en signal lectrique. Un traducteur est caractris par sa frquence, le type donde quil
gnre et la forme de son faisceau acoustique.
1. dans ce cas, le mme capteur peut ventuellement assurer les rles dmetteur et de rcepteur
1.2. Mesures de champs acoustiques : les ultrasons 9
La frquence choisie doit satisfaire deux condions ncessaires au respect de lhypothse

des milieux continus lors du contrle non destructif (Castaings [51]) :
le milieu solide doit tre homogne lchelle des longueurs donde utilises. Pour
cela, il suffit que la plus petite longueur donde des modes gnrs soit suprieure
aux dimensions des htrognits de la structure. Si le milieu nest pas homogne,
il est dispersif et la propagation des ondes gnres est rgie par des phnomnes de
rfraction.
le milieu doit pouvoir tre assimil un milieu semi-infini, cest--dire que les chos
des divers modes gnrs doivent tre parfaitement spars, aprs transmission de
lchantillon, pour permettre la mesure du temps de vol (temps de passage de londe
dans lchantillon).
Le faisceau acoustique du traducteur se dcompose classiquement en deux parties :

le champ proche (dit zone de Fresnel) et le champ lointain (dit zone de Fraunhofer).
De manire saffranchir de linstabilit du champ proche tout en ayant le maximum
2
dintensit acoustique, il convient de positionner le capteur D4 de la surface de la pice
inspecte (figure 1.2).
Figure 1.2 Intensit acoustique sur laxe du faisceau (D est le diamtre de la source
utilise et la longueur donde) (Dumont-Fillon [71])
1.2.1.2 Les diffrents types dondes
Suivant lorientation du capteur (do implicitement celle du transducteur), plusieurs

types dondes peuvent tre gnres (Colin [62], Dumont-Fillon [71], Lefebvre et al. [133],
figure 1.3) :
les traducteurs droits produisent des ondes de volume longitudinales,
les traducteurs dangles permettent de produire des ondes de volume transversales,

des ondes de surface ainsi que des ondes de plaque.
Figure 1.3 Types dondes et modes de propagation des ultrasons dans un solide
(Dumont-Fillon [71])
Les ondes de volume longitudinales dites ondes de compression se caractrisent par un

dplacement des particules lmentaires parallle au dplacement de londe (figure 1.4-a).
Ces ondes font intervenir des contraintes de traction-compression (figure 1.5-a) et se pro-
pagent aussi bien dans les liquides, les gaz et les solides. Elles permettent de dtecter des
dfauts parallles la surface de dtection, ce qui savre particulirement pertinent dans
le cas de composites multicouches puisque les dfauts y sont majoritairement parallles
aux plis.
Dans le cas des ondes de volume transversales dites ondes de cisaillement, le dpla-
cement des particules lmentaires se fait perpendiculairement au dplacement de londe
(figure 1.4-b). Celles-ci induisent des contraintes de cisaillement-torsion (figure 1.5-b) et
ne peuvent tre transmises par des liquides ou des gaz. Leur utilisation est adapte la
dtection de dfauts non parallles la surface de dtection tels que des kissing bonds
(Adams and Drinkwater [2]).
(a)
(b)
Figure 1.4 Ondes de volume : (a) ondes de compression, (b) ondes de cisaillement
(Lefebvre et al. [133])
Les ondes de surface dites ondes de Rayleigh rsultent quant elles de la combinaison
dune onde longitudinale et dune onde transversale en quadrature se propageant la
surface libre du matriau. Le dplacement des particules lmentaires suit une orbite
elliptique dont le grand axe est perpendiculaire la surface inspecte et dont le petit axe
est parallle la direction de propagation de londe qui correspond au contour de la surface
contrle (onde guide). Ces ondes ne prennent naissance qu la surface du matriau et
permettent de dtecter des dfauts trs proches de la surface mais non dbouchants ou
qui se situeraient dans la zone de Fresnel de certains traducteurs (Lefebvre et al. [133]).
(a)
(b)
(c)
Figure 1.5 Ondes et contraintes induites : (a) onde volumique de compression, (b)
onde volumique de cisaillement, (c) onde de surface (Mistou [143])
On distingue enfin les ondes de plaque dites ondes de Lamb qui se propagent parall-
lement aux surfaces tudies (onde guide galement) et qui rsultent de la superposition
de deux ondes longitudinales et de deux ondes transversales. Ces ondes sont utilises
pour contrler des plaques de faible paisseur (quelques longueurs dondes au maximum)
puisque leur caractristique est de se propager dans lpaisseur dune paroi (Cherfaoui
[55]). Elles sont ainsi adaptes la dtection de dfauts se situant au niveau de la colle
entre 2 plis ou entre un pli et une me car les ondes de plaques sont galement sensibles
aux proprits mcaniques des adhsifs (Adams and Drinkwater [2]), lvaluation de d-
laminages produits par un impact ainsi que de dcollements de semelle domga (Yashiro
et al. [198]).
1.2.1.3 Lattnuation ultrasonore
Les composites stratifis sont gnralement dots de proprits dabsorption non ngli-
geables qui sont dues au comportement viscolastique des matrices. De plus, leur structure
multi-couche entraine des effets de diffraction des ondes dans la structure interne (inter-
faces solide-solide imparfaites, porosits, fibres, fissures et autres dfauts). La combinaison
de ces deux phnomnes entrane lattnuation de lamplitude de londe ultrasonore lors
de sa propagation dans le composite, soit une perte dans lamplitude du signal de sortie
par rapport celle du signal dentre.
Lanalyse de lattnuation ultrasonore est utilise dans les techniques de contrle non
destructif par ultrasons (cho ultrasonore, transmission ultrasonore). Dans le cas par
exemple de porosits se situant dans un joint de colle entre 2 plis ou entre un pli et
une mousse, la comparaison de lattnuation mesure sur la pice par rapport celle
obtenue sur des prouvettes talons non poreuses permet dvaluer quantitativement le
taux volumique de vides (Boyer [46], Ledru [131], Phani [165]). Il est important de noter
que la rsolution de la dtection par contrle ultrasonore samliore en augmentant la
frquence de londe ultrasonore, et en consquence lpaisseur de la structure pouvant
tre contrle diminue (Bar-Cohen et al. [28]).
A partir des principes gnraux sur les ondes ultrasonores qui viennent dtre voqus,
on se propose prsent de dtailler les mthodologies de contrle par ultrasons les plus
utilises sur les structures composites.
1.2.2 Principe de contrle

1.2.2.1 Echo ultrasonore : contrle en face avant
Cette mthode est trs approprie au contrle des pices ne possdant quune surface
accessible puisquelle ne ncessite lutilisation que dun seul capteur jouant la fois le
rle dmetteur et de rcepteur. Deux typologies de capteurs ultrasonores peuvent tre
employs (Framezelle et al. [79]) :
les capteurs mono-lments contenant un seul lment pizo-lectrique mettant et
recevant londe et qui sont les plus employs dans lindustrie,
les capteurs multi-lments appels aussi phased array, composs de plusieurs l-
ments individuels pouvant tre pilots indpendamment ; ils permettent une focalisa-
tion lectronique du faisceau ultrasonore diminuant les zones mortes dinvestigation.
Dun point de vue pratique, une couche dagent couplant est place entre le capteur et
la pice de faon guider le passage de londe. Ce couplant peut tre de nature diverse
(Bar-Cohen et al. [28], Mouritz et al. [151], Zani [199]) :
un corps gras, un gel ou de leau (ventuellement additionne dun agent mouillant)
en petite quantit positionne sur la pice ; on parle alors de mthode par contact,
de leau (ventuellement additionne dun agent mouillant) en grande quantit, par
immersion partielle (capteur et une partie de la pice) ou totale (capteur et pice)
dans une cuve.
Lpaisseur de lagent de couplage doit absolument rester constante tout au long du
contrle afin dassurer une inspection complte et homogne de la zone balaye par le
traducteur sur la structure. Lorsque la taille de la cuve le permet, ceci justifie le recours
prfrentiel aux mesures en immersion qui bnficient dautre part de systmes de guidage
automatique. Dans certains cas, il est ncessaire de placer en plus un sabot entre la pice
et le traducteur afin dloigner lentre du signal dans la pice de la zone de Fresnel, et
ainsi permettre un contrle optimal ds les premiers millimtres.
Lanalyse des dfauts est ensuite base sur le principe du radar : les ondes ultraso-
nores mises dans la pice sont rflchies sur les interfaces principales du matriau (entre
et fond) ainsi que sur les htrognits internes produisant des chos par rupture dim-
pdance acoustique. La distance entre les diffrents chos par rapport lcho dentre
dans la pice est relie directement aux temps de propagation aller-retour de londe dans
lpaisseur de celle-ci. Le rsultat mesur se prsente alors sous forme dun chogramme
(galement appel A-scan) traduisant lamplitude (en dB) des chos en fonction de la
distance traverse par londe (Mouritz et al. [151]) (figure 1.6) :
lcho dmission reprsentant le lieu dmission du train dondes,

lcho dentre 1 reprsentant la pntration du train dondes dans la pice inspec-
te ; si le traducteur nest utilis quavec un couplant acoustique de faible paisseur
(mthode par contact), lcho dentre est confondu avec lcho dmission dans la
pice (figure 1.7).
lcho de dfaut 2 reprsentant la rflexion de londe sur limperfection (si elle existe)
(figure 1.6),
lcho de fond 3, toujours situ en dernier dans lchogramme, reprsentant la r-
flexion de londe sur le fond de la pice (figure 1.6).
Figure 1.6 Mthode cho ultrasonore par immersion ou par contact avec sabot (Cher-
faoui [55], Lambert et al. [121])
Figure 1.7 Mthode cho ultrasonore par contact (Framezelle et al. [79])
Lobservation de ces chos permet de localiser les dfauts internes prsents dans une
structure mais galement den raliser des mesures dimensionnelles. Un traitement du
signal peut ventuellement tre ralis pour amliorer lanalyse, comme par exemple pour
la mthode -6dB, qui consiste faire un crtage des amplitudes pour enlever une partie
du bruit de mesure. La taille des dfauts est ensuite estime par la comparaison entre
les amplitudes des chos de fond et de dfaut ; la profondeur des dfauts est quant elle
dtermine par le temps de vol entre lcho dentre (ou de fond) dans la pice et lcho du
dfaut (Bar-Cohen et al. [28]) (figure 1.6). Pour cela, il faut bien sr talonner avant tout
test lappareil de mesure avec une plaque talon reprsentative de la structure inspecte.
Cette plaque est compose de dfauts artificiels ayant pour but de reprsenter au plus prs
les dfauts devant tre identifis aussi bien en terme de nature, quen terme de dimension.
Lutilisation dune plaque talon permet de raliser les diffrents rglages ncessaires la
visualisation des dfauts signals comme critiques pour la structure contrle.
Il existe traditionnellement deux cartographies 2D interprtes dun contrle ultraso-
nore, savoir le B-scan sous la forme dune coupe dans lpaisseur de la pice (figure
1.8) et le C-scan sous la forme dune projection en surface des dfauts prsents dans la
pice contrle (Bar-Cohen et al. [28], Perdijon [161]) (figure 1.8). Le couplage de ces
trois vues permet davoir une visualisation 3D de la rpartition des dfauts dans la pice
contrle. On peut enfin signaler que lemploi de sondes multi-lments permet daccder
des vues sous des angles diffrents (par exemple B-scan angulaire appel aussi le S-scan
(Fidahoussen et al. [78])).
Figure 1.8 Diffrents modes de prsentation de linspection ultrasonore (Dumont-Fillon

[71])
Lcho ultrasonore permet didentifier essentiellement les dfauts suivants dans les
composites dus :
la fabrication : porosits, dcollements, variations du taux dimprgnation (Boro-
Djordjevic [44], Ledru [131], Wong and Ron [196]),
et/ou leur utilisation en service : dlaminages, impacts (Boro-Djordjevic [44], Mouritz
et al. [151], Wong and Ron [196]).
1.2.2.2 Transmission ultrasonore : contrle en face avant et face arrire
La transmission ultrasonore utilise deux traducteurs, lun jouant le rle dmetteur et

lautre de rcepteur. Cette technique se base sur des mcanismes de transmission dondes
ou de diffusion et rflexion suivant la configuration choisie (Chaix et al. [53], Cherfaoui
[55], Framezelle et al. [79], Zani [199]) :
en mode transmission, metteur et rcepteur sont placs de part et dautre de la
pice inspecte (figure 1.9) ; cette approche est la plus classique car la plus simple
mettre en oeuvre et gnralement plus prcise ;
en mode mission-rception, les deux traducteurs sont du mme ct de la pice.
Figure 1.9 Contrle en mode transmission (Chaix et al. [53])
Dans le cas de cette technique, la dtection des dfauts se ralise par mesure de latt-
nuation ou coupure directe du faisceau ultrasonore suivant les propits dabsorption du
dfaut intercept (figure 1.10). Gnralement, ltalonnage de lappareil de test se ralise
vide ou avec un matriau sain. Avant chaque test, il convient de sassurer que les cap-
teurs soient correctement positionns pour que le signal soit mis et rceptionn de faon
optimale. Une fois le test ralis, les rsultats sont bien sr visualisables par les 3 vues
(A-scan, B-scan et C-scan) dj dtailles dans la partie prcdente.
Figure 1.10 Principe de dtection par transmission (Lacroix [119])
Les difficults lies cette technique sont les modifications de londe ultrasonore qui
ne sont pas dues la prsence des dfauts et qui, de fait, perturbent lanalyse : problme
li lagent couplant, dsalignement des deux traducteurs, changement de structure dans
la pice contrle par exemple. Par consquent, cette technique est gnralement associe
aux techniques dcho ultrasonore par contact ou immersion.
Lagent couplant utilis pour la transmission ultrasonore peut tre un corps gras, de
leau (en immersion ou par jet deau) ou de lair. Le recours au jet deau est trs frquent
dans le milieu industriel, associe gnralement un systme automatique de dplacement
des traducteurs assurant une liaison acoustique optimale (Hsu [107], Robert and Green
[176]). Toutefois, face la contamination potentielle des surfaces de la structure par ce
moyen (endommagement, corrosion) et autre perturbation induite par les changements de
milieux dans la propagation de londe, des mthodes sans contact et sans milieu couplant
tendent se dvelopper. On peut citer par exemple les trois suivantes qui peuvent tre
ventuellement couples entre elles (Cawley [52], Clorennec [61], Eyma et al. [75], Kazys
et al. [116], Robert and Green [176], Ruzek et al. [179], Zani [199]) :
les ultrasons couplage air, qui sont particulirement adapts au contrle des
structures hydrophobes ou sensibles leau telles que composites sandwichs (Hsu
[107], Kazys et al. [116]) ; on peut noter que les capteurs coupls lair sont adap-
ts la gnration et la dtection dondes de Lamb en mode mission-rception
(Framezelle et al. [79], Hsu [107]),
les ultrasons lasers, gnres par un laser et dont la dtection est souvent ralise par
une sonde interfromtrique comme le systme de Fabry-Perot (Cawley [52], Robert
and Green [176], Voillaume and Campagne [194]) ; dans cette dmarche, la prsence
dune surface rflchissante savre ncessaire (par exemple, un film mince peut tre
coll cette fin en surface de composites en carbone-poxy) ;
ultrasons avec traducteurs magnto-acoustiques.
Ces nouvelles technologies sont de plus en plus convoites dans les inspections sur struc-
tures composites (Cawley [52], Zani [199]). Outre les avantages prcits, elles minimisent
en effet les bruits de mesures (impliquant une bonne reproductibilit), permettent lins-
pection de pices chaudes ou gomtrie complexe comportant de multiples courbures
ainsi quun contrle sur la ligne de fabrication.
Parmi les dfauts sur matriaux composites pouvant tre identifis par la transmission
ultrasonore, on peut citer galement :
des dfauts de fabrication : porosits, dcollements entre lme et les peaux de
matriaux sandwichs Hsu [107],
et/ou leur utilisation en service : dlaminages, fissures Hsu [107], Mouritz et al.
[151].
1.2.3 Bilan
Le contrle ultrasonore est aujourdhui la technique la plus utilise et la plus recon-
nue dans lindustrie aronautique au travers des standards AITM. Parmi ses avantages
principaux, il faut citer :
mme si les reprsentations planes sont les plus couramment employes, la possibilit
de dtecter des dfauts dans le volume de la pice (par combinaison des vues 2D),
une grande sensibilit de dtection des dfauts plans, grce lutilisation trs r-
pandue des ondes longitudinales incidence normale,
lautomatisation possible du contrle permettant dobtenir une liaison acoustique
optimale,
le caractre portable de certains dispositifs exprimentaux.
1.3. Mesures de champs densimtriques : radiographie rayons X et
tomographie rayons X 19
Toutefois, cette technique prsente quelques inconvnients dont notamment une mise
en oeuvre difficile et dlicate, le besoin le plus souvent de lagent couplant ainsi que la
ncessit de raliser une plaque talon reprsentative de la pice inspecte et des dfauts
recherchs. Par ailleurs, il existe certaines contraintes dapplication dues lemploi de
traducteurs spcifiques selon le matriau test, son paisseur, la forme de la structure
ainsi que de son tat de surface. Enfin, il convient de signaler la sensibilit de la mthode
par rapport au type dondes choisi (les rsultats obtenus pour chaque type diffrent selon
lorientation de la surface du dfaut) ainsi qu leur frquence.
1.3 Mesures de champs densimtriques : radiographie

rayons X et tomographie rayons X
Depuis de nombreuses annes, les mesures de champs densimtriques sont employes
dans le milieu industriel (Babout [22], Badel [24], Lecomte [130]). Dans lindustrie aro-
nautique, sa pratique est normalise, entre autres par les normes AITM (AITM6-7002
[12], AITM6-7006 [13]). Au fil des annes, beaucoup de dveloppements ont t raliss
sur les mesures de champs densimtriques comme moyen de CND pour les structures com-
posites, sous limpulsion notamment du secteur mdical. En particulier, la radiographie X
est devenue maintenant incontournable (Sreelekshmi et al. [187]) alors que la tomographie
X, qui constitue une extension 3D de la radiographie, est en plein essor (Badel [24]).
Les mesures de champs densimtriques sont bases sur :
lenvoi dun flux dondes lectromagntiques (rayons X) au travers de la pice,
la rception du flux transmis par un ou des dtecteurs (contrle en face arrire),
lobtention dune mesure de champ en cartographie 2D (radiographie X) ou 3D aprs
traitement dimages et reconstruction (tomographie X).
Compte tenu de leur base physique commune, on fait en premier lieu une prsentation
conjointe des grandes lignes de ces deux techniques avant daborder leurs spcificits
propres.

1.3.1.1 Le rayonnement X
Les rayons X sont des rayonnements lectromagntiques de trs forte nergie, de

quelques keV plusieurs MeV, cest--dire de longueur donde trs faible (de 1010
1014 m) par rapport au rayonnement visible (voir figure 1.11). Lintraction des rayons
X avec la matire conduit la production dlectrons responsables dionisations dans le

milieu travers. Ils sont classs dans le domaine des rayonnements ionisants, mais ils le
sont en fait indirectement. Leurs effets dans le milieu parcouru proviennent des ionisa-
tions provoques par les particules secondaires. Ce sont ces effets qui sont utiliss pour la
dtection et la visualisation des rayonnements.
Figure 1.11 Spectre lectromagntique (Garnier [83])
1.3.1.2 Loi dattnuation des rayons X
Lorsque des rayons X traversent un milieu, on note une perte par rapport au nombre
de particules incidentes, soit en dautres termes une variation dintensit du faisceau. Ce
phnomne, appel attnuation, rsulte de linteraction des photons avec la matire au
travers de plusieurs processus distincts (Le-Bras [125], Freud [81], Lecomte [130]) :
interaction avec les atomes du milieu (effet photolectrique),

interaction avec les lectrons des atomes (effet Compton),
interaction avec le noyau des atomes (effet Thomson-Rayleigh) ; vis--vis des nergies
utilises dans la radiographie (entre 10 et 150 keV), cet aspect est habituellement
nglig car il ne contribue que trs faiblement lattnuation globale.
Durant une interaction, lnergie du photon incident est compltement ou partiellement

absorbe et transfre un lectron sous forme dnergie cintique. Dans la figure 1.12, les
diffrents parcours et interactions possibles pour les photons X lors de la traverse dun
matriau sont dtaills.
Figure 1.12 Interaction entre le rayonnement et la matire lors de la traverse dun

matriau (Lecomte [130])
La loi dattnuation des rayons X dpend de lpaisseur de lobjet et du coefficient

massique dattnuation du matriau (Babout [22], Badel [24], Bayraktara et al. [32], Le-
Bras [125], Carmona [50], Lecomte [130], Lenoir [134], Bossi et al. [45]). Ce dernier point
reprsente tout lintret de ces mesures dattnuation des radiations X dans le cadre du
CND par la dpendance la densit du matriau. Il sagit ainsi dune mthode de mesure
de champ densimtrique : chaque pixel (radiographie) ou voxel (tomographie) caractrise
en effet la moyenne des densits massiques des diffrents constituants compris dans celui-
ci et permet donc de diffrencier localement les htrognits telles que les vides ou
inclusions. On note que lnergie du faisceau de photons intervient galement dans cette
analyse puisque le coefficient dattnuation du milieu dcrot lorsque lnergie des photons
augmente (figure 1.13).
Figure 1.13 Evolution du coefficient dattnuation dun milieu en fonction de lnergie

E des photons
1.3.2 Principe de contrle

1.3.2.1 Dmarche globale
De manire gnrale, le contrle par rayons X se ralise en transmission (contrle en

face arrire, figure 1.14). La chane de contrle comprend la source de rayons X, le support
de lobjet contrl ainsi que le dtecteur avec ou sans systme imageur.
Figure 1.14 Principe de contrle par rayons X (Lecomte [130])
Les rsultats du contrle par rayons X tant bass sur les diffrences dabsorption de
la pntration de ces radiations dans le matriau inspect, deux paramtres conditionnent
fortement la rsolution des images obtenues et donc la performance de la technique en
terme de dtection :
dune part la capacit du matriau test absorber les radiations X, soit limpor-
tance de son coefficient massique dattnuation ; pour les composites en carbone dont
le coefficient dabsorption est faible la base, il est ainsi recommand dutiliser une
faible nergie de radiation (les tensions les plus couramment employes se trouvent
entre 10 keV et 50 keV mais lintervalle de 10 keV 20 keV est gnralement suf-
fisant) voire demployer un liquide pntrant radiographiquement opaque (comme
par exemple le ttrabromthane) si lon souhaite dtecter des dfauts dbouchants
et connects entre eux (Schilling et al. [181], Wong and Ron [196]),
dautre part le constraste dabsorption entre le matriau dit sain et les dfauts au
sein de celui-ci, soit le rapport entre le nombre N2 de photons X transmis par le
dfaut et le nombre N1 de photons X transmis par son voisinage (voir figure 1.15
pour une pice dpaisseur constante comportant un vide, Lecomte [130]) ; de faon
vidente, le contraste sera dautant meilleur que le dfaut sera pais.
Figure 1.15 Contraste dun dfaut sur limage radiante (Lecomte [130])
Au niveau de la rsolution spatiale des images obtenues (soit la taille des pixels/voxels),
un certain nombre de paramtres entre en jeu et de manire diffrente suivant la configu-
ration matrielle du systme :
la source de rayons X (taille, type, nergie des rayons),
le dtecteur (nature, dimensions, rsolution spatiale, grossissement 2 ),
lchantillon (taille et forme).
Du point de vue du CND, la qualit de la dtection sera lie une rsolution spatiale
de la mesure infrieure la taille habituelle des dfauts recherchs en ultrasons ; cela
conditionnera donc une grande partie des rglages du dispositif exprimental (Lenoir
[134]) ainsi que les temps dacquisition ncessaires.
1.3.2.2 La radiographie rayons X
La radiographie X correspond lapplication bidimensionnelle de la dmarche prcdente.

La cration des rayons X est ralise dans ce cas par un tube rayons X dans lequel
les lectrons sont acclrs sous une diffrence de potentiel importante venant percuter
une cible mtallique gnralement en tungstne (Freud [81], Lecomte [130]) (figure 1.16).
2. Le grossissement est le rapport entre la taille de limage de lchantillon enregistre par le dtecteur
et la taille relle de lchantillon ; plus il est important, meilleure est la rsolution spatiale.
Lors de cette acclration, lnergie des photons va varier de quelque keV lnergie
correspondant la haute tension choisie. Les interactions des lectrons avec la cible sont
de deux types (Le-Bras [125], Lecomte [130]) :
par collision des lectrons avec les lectrons des atomes de la cible,
par collision et freinage des lectrons dans le champ lectrique des noyaux de la
cible.
Figure 1.16 Schma dun tube rayons X (Badel [24], Le-Bras [125], Lecomte [130])
Parmi les diffrents supports de visualisation des rsultats, le plus frquemment utilis
est le film radiographique (film argentique). Ce dernier comprend un support transpa-
rent dont les deux faces sont recouvertes de couches dmulsion comportant des cristaux
dhalognures dargent. Le rayonnement X y produit des actions de diffrents types : tout
dabord larrachement dlectrons et ensuite des ractions chimiques entre ces lectrons
librs et les cristaux de bromure dargent prsents dans le film. Un traitement chimique
de dveloppement transforme compltement ces cristaux modifis en argent mtallique
noir, les autres tant limins. Ce traitement peut tre effectu de faon manuelle ou au-
tomatiquement en machine. Laction du rayonnement X sur le film dpend, entre autres,
de la dure de lexposition et de son intensit (Lecomte [130]).
En ce qui concerne les matriaux composites, les principales applications de la radio-
graphie X concernent la dtection dinclusions comme des inserts mtalliques (Avril et al.
[18]), des manques de colle (Herve et al. [105]) ou dendommagements locaux tels que des
fissures et dlaminages (Bertin et al. [37], Thomas [191]). Nanmoins, les mesures fournies
par cette technique restent une projection 2D sur le dtecteur des dfauts prsents dans le
volume de la structure contrle. Afin de pouvoir complter cette information et localiser
les htrognits en profondeur, il convient de se tourner vers lextension 3D de cette
mthode : la tomographie X.
1.3.2.3 La tomographie rayons X
La tomographie X ou tomodensitomtrie permet une mesure volumique du champ

densimtrique. Cette extension tridimensionnelle sobtient partir de plusieurs images ra-
diographiques 2D appeles coupes tomographiques (figure 1.17) (Babout [22], Bayraktara
et al. [32], Carmona [50], Lenoir [134]). Ces coupes sont acquises par un systme imageur
numrique et sont composes dune matrice de pixels sur une paisseur finie de lobjet.
Lassociation de ces diffrentes coupes (phase de reconstruction) conduit limage volu-
mique compose de voxels dont la taille implique la rsolution spatiale de la mesure (figure
1.18).
Figure 1.17 Schma de principe du contrle par tomographie X : du faisceau jusquaux

coupes (Carmona [50])
Figure 1.18 Schma dune coupe (Lenoir [134])

Classiquement, la mesure par tomographie est effectue laide de scanners pour ap-
plications mdicales ou industrielles. Ces appareils, bass sur le mme principe de cration
des rayons X que prcdemment, peuvent tre distingus selon lagencement des couples
metteur/dtecteur (figure 1.19) (Lenoir [134]) :
La premire gnration de scanners utilise un faisceau en pinceau et un seul d-
tecteur, situs dans le plan de coupe ; la mesure est effectue par translation de
lensemble source/dtecteur sur toute la largeur de lchantillon puis en recommen-
ant lopration aprs rotation de lchantillon ou de lensemble source/dtecteur,
Dans la deuxime gnration, la source subit galement une translation mais se

prsente en ventail avec un faible angle douverture et une ligne de dtecteurs ;
lavantage est dobtenir plus dinformations en moins de temps par rapport un
scanner de premire gnration ; il permet galement dobtenir des donnes plus
prcises que les autres gnrations utilisant une source en ventail mais pour un
temps dacquisition plus long,
Les scanners de troisime et quatrime gnrations utilisent enfin une source et des
dtecteurs plus larges permettant la plupart du temps de couvrir lchantillon en
entier ; la source est en ventail ou en cne, et les dtecteurs sont respectivement
en ligne ou en plan ; ces configurations rduisent considrablement le temps dac-
quisition notamment la configuration cne/plan qui permet de scanner plusieurs
coupes en mme temps ; la diffrence entre les deux gnrations rside dans le fait
que seul lchantillon tourne dans la troisime gnration alors que tout lensemble
source/dtecteur tourne dans la quatrime (si lensemble des dtecteurs est circu-
laire comme dans la figure 1.19) ; les configurations o la rotation est effectue par
lchantillon permettent damliorer la qualit des images car des filtres peuvent tre
facilement placs devant la source (ou les dtecteurs) afin de rduire les artefacts.
Figure 1.19 Gnrations de scanners (Lenoir [134])

La tomographie X peut tre galement ralise laide de synchrotrons. Dans ce cas,

le processus de cration des radiations diffre des scanners mdicaux et industriels et
consiste un changement de trajectoire dans un acclrateur particules (Lenoir [134]).
Les lectrons sont acclrs une vitesse proche de la lumire (lectrons relativistes) puis
dvis par un champ magntique (force de Lorentz) ce qui cre lmission dun faisceau
de photons tangentiellement la courbe suivie par les lectrons. Par rapport aux tubes
rayons X, les nergies mises sont plus importantes, le spectre dnergies dlivrables est
plus large et les faisceaux peuvent tre considrs comme convergents donc plus efficaces,
pour une meilleure rsolution spatiale.
La diversit des angles dincidence des rayons X et la performance des dtecteurs de

cette technique dimagerie conduisent gnralement de meilleurs rsultats que la radio-
graphie X avec films argentiques (Bayraktara et al. [32]). Parmi les facteurs influents sur
la qualit de la mesure par tomographie, on peut ainsi citer le type et la gnration de
scanners utilis ainsi que la rsolution des dtecteurs (soit la taille des cellules photolec-
triques qui enregistrent les rayons transmis). A titre dexemple, les scanners mdicaux
prsentent des rsolutions spatiales de lordre de plusieurs centaines de microns, les scan-
ners industriels de lordre de quelques dizaines de microns et les synchrotrons de lordre du
micron. En terme de rapidit dacquisition, lappareil le plus lent est le scanner industriel
(jusqu plusieurs heures), et les plus rapides sont le scanner mdical et le synchrotron
(quelques minutes).
Il convient enfin de signaler que les rsultats obtenus sont particulirement sensibles
au principe de reconstitution des images. La mesure radiographique de lattnuation des
rayons X travers un objet est une mesure projete pour un angle incident donn. Afin
daccder lattnuation en chaque point du volume, il est ncessaire dassocier un en-
semble de projections sous diffrents angles rgulirement espacs (figure 1.20 (Lenoir
[134])). Lalgorithme utilis pour raliser la reconstitution doit tre capable de corriger les
erreurs engendres par les diffrences entre les conditions relles de mesure et les condi-
tions thoriques de la reconstruction afin de rduire ou supprimer les artfacts qui peuvent
masquer ou modifier des zones dintrt.
Figure 1.20 Principe de la mthode de reconstruction dimages tomographiques (Lenoir

[134])
En ce qui concerne les matriaux composites, la tomographie X est utilise pour quan-
tifier la prsence de dfauts induits par le procd de fabrication ou par des sollicitations
externes, par exemple :
des porosits (Babout [22], Bayraktara et al. [32], Hanke et al. [104], Schilling et al.
[181]),
la distribution de lorientation des fibres (Bernasconi et al. [36]),
des fissurations (Moffat et al. [145], Schilling et al. [181]),
des dlaminages et leurs contours (Bayraktara et al. [32], Geier [85], Hanke et al.
[104], Schilling et al. [181], Tsao and Hocheng [192]).
On notera galement que cette technique permet le suivi de ces dfauts sous contraintes
lorsque la taille et le temps dacquisition des tomographes sont respectivement compatibles
avec la dimension des essais et les vitesses de phnomnes (Babout et al. [23], Carmona
[50], Lenoir [134], Maire et al. [136, 137], Moffat et al. [145], Reifsnider and Williams
[174]).
1.3.3 Bilan
Les mesures de champs densimtriques par radiographie X et tomographie X per-
mettent daccder de manire non intrusive des informations sur la microstructure in-
terne des matriaux, en particulier de rvler toute htrognit ou singularit. Dans le
1.4. Mesures de champs thermiques : la thermographie infrarouge 29
cadre du CND, il sagit donc de techniques particulirement intressantes dont les princi-
paux avantages sont :
leur applicabilit la quasi-totalit des matriaux,
la possibilit de contrler des pices gomtrie complexe,
la dtection volumique de dfauts (inclusions gazeuses, inclusions de matire tran-
gre,...).
Par rapport la radiographie qui noffre quune visualisation 2D du dfaut, la tomo-
graphie a permis des avances notables grce une information dans le volume et la
localisation des imperfections dans lpaisseur des pices. Dans les deux cas, on peut
nanmoins citer comme limitations potentielles les contraintes de protection ncessaires
lors de la ralisation du contrle, la matrise technique indispensable leur utilisation,
le cot des appareillages et leur taille qui ne permet gnralement pas de contrler des
pices de grande envergure.
1.4 Mesures de champs thermiques : la thermographie

infrarouge
Dans le cadre du CND sur structures mcaniques de faible paisseur, une nouvelle
technologie sans contact, rapide en terme dexcution et permettant de raliser des acqui-
sitions de champs en une prise merge depuis quelques annes. Il sagit de la thermogra-
phie infrarouge qui permet daccder des mesures de champs thermiques (temprature
ou proprits thermiques, Bremond [47], Chrysochoos [57], Guibert [102], Lambert et al.
[121]). Cette technique commence tre prsente dans lindustrie puisque sa pratique est
dj normalise par les standards AITM (AITM6-0015 [6]).
Elle est base sur :

lexcitation interne ou externe de la pice,
la rception du flux de chaleur restitu en rponse par celle-ci laide dun dtecteur
infrarouge, en face avant ou arrire,
lobtention dune mesure de champ thermique en cartographie 2D (qui permet
didentifier la projection en surface des dfauts).
1.4.1.1 Rayonnement InfraRouge (IR) et bilan radiatif
Le rayonnement thermique est une mission dondes lectromagntiques caractrises

par des longueurs dondes dites infrarouges comprises entre 3 m et 1000 m (Clerjaud
[60], Garnier [83], Guibert [102], Pastor [159]) (figure 1.11). Le spectre de lumire infra-
rouge se dcompose quant lui en trois domaines :
le domaine infrarouge proche (de 0,759 3 m),

le domaine infrarouge moyen (de 3 5 m),
le domaine infrarouge lointain (de 5 m 1000 m),
les gammes les plus frquemment tudies tant les domaines proche et moyen.
Lorsquun flux thermique est envoy sur la pice (figure 1.21), une partie r du
rayonnement est rflchie sur le corps, une autre partie t est transmise au travers de
celui-ci et une dernire partie a est enfin absorbe. En absorbant ce flux a , le matriau
schauffe et le flux de chaleur e quil met est fonction de sa temprature linstant
considr. Si une sollicitation interne induit cet chauffement, la dissipation produite est
de la mme faon lie sa temprature. Cette transcription en temprature se fait partir
de la connaissance de lmissivit (rapport dmission entre le matriau et un corps de
rfrence dit corps noir 3 ) qui caractrise laptitude du matriau mettre de lnergie par
rayonnement. Lmissivit du matriau, comprise entre 0 et 1, est une proprit thermique
fonction notamment de sa densit, de son tat de surface, de la temprature et de la
longueur donde (Guibert [102], Ibarra-Castanedo [110], Pajani and Audaire [156]).
3. Un corps noir est un matriau idal opaque dmissivit gale 1, qui absorbe intgralement tous les
rayonnements qui lui parviennent, sans aucune rflexion, et indpendemment de la direction dmission
et de la longueur donde.
Figure 1.21 Bilan radiatif (Pajani [155])
1.4.1.2 Capteurs infrarouges
Il existe essentiellement deux types de capteurs du rayonnement infrarouge (Guibert

[102]) :
les capteurs thermiques comme par exemple les microbolomtres, qui sont des r-
sistances ou des jonctions voluant avec la chaleur incidente,
les capteurs photoniques, plus communs du fait de leur grande sensibilit et de leurs
bonnes rsolutions spatiales, qui mesurent directement la densit du flux de photons.
Pour raliser des analyses dimagerie thermique, les capteurs infrarouges sont disposs
en matrice au sein de la camra, chaque capteur formant un pixel de limage. Pour les
mesures infrarouges, des lentilles en silicium ou germanium sont employes car elles sont
transparentes ces longueurs donde (Guibert [102]).
Une mesure de temprature effectue laide dune camra infrarouge demande cer-
taines prcautions. Outre des rflexions optiques parasites et des effets de lenvironnement,
les capteurs infrarouges dans la camra peuvent tre perturbs par la temprature de la
lentille. Ainsi les bords de limage peuvent tre plus chauds (respectivement froids) que
son centre lorsque lobjet observ est plus froid (resp. chaud) que la temprature de la
camra. Pour limiter au maximum ces diffrents bruits, il conviendra par consquent de
travailler dans un laboratoire thermorgul o les vibrations environnantes sont minimi-
ses, de placer la camra et les pices examiner dans celui-ci quelques jours avant le
test et de raliser les essais dans le noir. Lorsque les mesures sont ralises sur de trs
faibles dissipations thermiques, une calibration est par ailleurs recommande avant chaque
srie dexpriences (Guibert [102]). Enfin, quel que soit le sens de contrle adopt (voir
ci-aprs), lobjectif de la camra doit tre positionn paralllement la surface examine
de faon optimiser le champ optique et viter tout effet de distorsion de limage.
1.4.2 Principe du contrle

Le CND par thermographie consiste mesurer le rayonnement infrarouge mis par la
pice suite sa sollicitation et de fournir lvolution temporelle du champ de temprature
ou des proprits thermiques associes. Lexistence dhtrognits dans ce champ traduit
alors des diffrences locales de transmission et de restitution du flux thermique induites
par des variations de densits (soit la prsence de dfauts internes et/ou dbouchants)
ainsi que des diffrences locales de rflexion (soit la variation de ltat de surface et/ou la
prsence de dfauts dbouchants) (Framezelle et al. [79], Larbi [122], Prodjinonto [170]).
Diffrentes mthodologies de contrle par thermographie ont t mises en place pour

analyser la rponse thermique des matriaux. Elles se distinguent notamment par la nature
et le type de signal de la (ou des) sources dexcitation ainsi que par la position respective
de la source et du capteur (sens du contrle). On se propose dvoquer ci-aprs ces deux
aspects puis de passer en revue trois techniques couramment employes.
1.4.2.1 Types de sollicitations
Les dispositifs exprimentaux de contrle comportent une ou plusieurs sources de

sollicitation induisant une dissipation au sein du matriau. Celles-ci peuvent tre de nature
diffrente (figure 1.22) :
en interne : la sollicitation par ultrasons (Framezelle et al. [79], Potet [168]) ou par
courants de Foucault (Maillard et al. [135]),
en externe : par exemple une lampe halogne, une lampe flash, un laser, un inducteur
haute frquence, une convection force dair ou encore une conduction par liquide
(Guibert [102], Pajani and Audaire [156]).
Dans la gamme des sollicitations externes par lampes, le signal mis par la source peut
prendre diffrentes formes (sinusodale, crneau, impulsionnelle, figure 1.23) correspon-
dant des mthodes de thermographie infrarouge (IR) particulires que lon se propose
de dtailler dans la suite.
Figure 1.22 Types de sollicitations thermiques (Maillard et al. [135])
1.4.2.2 Sens du contrle
Contrle en face avant

Pour le contrle en face avant, la source thermique et la camra sont positionnes
du mme ct de la pice (figure 1.24). Le flux mesur par la camra correspond aux
flux rflchis et restitus par la pice. Cette configuration est souvent employe quand les
dfauts sont supposs tre prs de la surface chauffe (Ibarra-Castanedo [110]). Elle se
Figure 1.23 Correspondance entre les mthodes de thermographie IR et leur signal

respectif
rvle galement attractive dans le cas dassemblages ferms puisquelle ne ncessite pas
davoir accs aux deux cts de la pice inspecte.
Figure 1.24 Thermographie infrarouge en face avant
Sajoutant au bruit environnant, il existe dans cette configuration deux autres bruits
de mesure importants :
les rflexions du rayonnement de la source, qui sera dautant plus accentu que la
surface tudie est brillante,
limpact thermique sur la surface au niveau du point focal de la source, qui acclre
la saturation de la pice, et donc celle de limage, rendant trs peu visible les dfauts
localiss cet endroit.
De plus, suivant la gomtrie de la pice contrle ainsi que la taille de la fentre de test
utilise, la notion dencombrement entre la source thermique et la camra IR peut devenir
un problme crucial pour la bonne ralisation du contrle.
Contrle en face arrire

Dans ce cas, la camra IR est positionne du ct oppos la source thermique (figure
1.25) ; la surface mesure est donc oppose la surface chauffe. Dans ce cas, le flux mesur
par la camra correspond aux flux transmis et restitus par la pice. Cette configuration,
assez peu utilise dans lindustrie, est prfrentiellement adopte lorsque les dfauts sont
supposs tre prs du ct de la camra (Ibarra-Castanedo [110]).
Figure 1.25 Thermographie infrarouge en face arrire
Dans cette configuration, un seul autre type de bruit est adjoint leffet de lenvi-
ronnement, savoir les rflexions de la brillance de ltat de surface sur lobjectif de la
camra.
On notera que les trois techniques de thermographie infrarouge qui sont dtailles
maintenant peuvent tre utilises soit en configuration face avant, soit en configuration
face arrire.
1.4.2.3 Thermographie InfraRouge (TIR)
La TIR est caractrise par le signal en crneau de la source thermique. Le schma de

cblage pour effectuer le contrle est assez simple (figure 1.26) : le flux thermique de la
source est gnr par lamplificateur de signal qui permet de faire varier son amplitude
(dans les limites de la puissance de la source thermique et de lamprage de lamplifica-
teur), le contrle (dclenchement du test, acquisition des donnes et arrt) tant gr par
le PC reli la camra IR.
Figure 1.26 Schma de cblage de la TIR
Cette technique permet dobtenir une cartographie 2D du champ thermique sur la

surface observe, dont les htrognits induites par les dfauts internes peuvent tre
dimensionnes. Bien que particulirement rapide et aise mettre en oeuvre, trs peu
dtudes dans la littrature ont t consacres lanalyse de structures composites par le
biais de cette mthode de thermographie IR.
1.4.2.4 Thermographie InfraRouge Module (TIRM)
La TIRM peut galement porter le nom de thermographie infrarouge lock-in, dtec-

tion synchrone ou encore radiomtrie photothermique (Clerjaud [60], Prodjinonto [170]).
Cette technique est caractrise par lutilisation dun signal priodique sinusodal de la
source thermique. Le schma de cblage de cette technique sarticule autour du gnrateur
de signal (figure 1.27) : celui-ci pilote la source thermique via lamplificateur de signal et
permet de faire varier la fois lamplitude mais aussi la frquence de londe envoye ; il
communique par ailleurs via la camra thermique sa frquence de travail au PC qui gre
comme prcdemment le contrle.
Figure 1.27 Schma de cblage de la TIRM

La TIRM conduit lacquisition de deux types de cartographies de dfauts :

celle du module (modification de lamplitude du signal), qui correspond la donne
obtenue par la mthode TIR, et qui permet de dimensionner la taille des dfauts
prsents dans la pice contrle,
celle de la phase (dphasage entre les signaux mis et mesur), qui conduit une
information la fois sur la dimension des dfauts mais galement sur leur profon-
deur dans la pice ; prcisment, on peut dduire de lanalyse de la cartographie de
phase la position relative des dfauts dans la profondeur (caractrisation qualita-
tive) ; lestimation directe (quantitative) de leur localisation respective requiert un
traitement des rsultats laide dalgorithmes spcifiques (Choi et al. [56], Guibert
[102], Ibarra-Castanedo [110]).
Outre des facteurs propres la pice (notamment lpaisseur de lprouvette ou la nature
des dfauts), la capacit de dtection de cette technique est influence par les paramtres
de londe mise (amplitude et frquence) Bai and Wong [27, 26, 25], Horny and Lannoy
[106], Meola et al. [138, 140]. Ces derniers peuvent donc tre modifis de faon capter
linformation voulue, ce qui est un avantage certain pour la mthode. En particulier, une
frquence faible augmente la profondeur de dtection.
Ce fonctionnement gnral de la TIRM a t appliqu laide de diffrentes sources

dexcitation :
partir de sources internes, gnres par exemple par une onde ultrasonore sinu-
sodale traversant la pice (mthode ultrasonique ou vibrothermographie (Gleiter
et al. [89], Potet [168], Zweschper et al. [202]), figure 1.28) ; pour lidentification de
fissures, on peut utiliser galement un champ lectrothermique induit par les cou-
rants de Foucault autour des pointes de celles-ci (Maillard et al. [135], Sakagami
and Kubo [180]) (mthode dite singulire) ;
partir de sources externes telles que les lampes au Xenon (Sakagami and Kubo
[180]) ("mthode disolation" permettant de dtecter des dfauts de type trous
fond plat dans des matriaux tels que la cramique et le graphite) ou plus classi-
quement des lampes halognes (Choi et al. [56]).
Wu (Wu and Busse [197]) prsente une tude comparative entre TIRM avec sources de cha-
leur interne (ultrasonore) et externe (lampes halognes) pour la dtection de dlaminages
dans du stratifi carbone-poxy et dimpacts dans de la cramique. Les rsultats sont
apparemment meilleurs avec la TIRM source interne. Toutefois, ces deux dmarches
semblent complmentaires puisque la TIRM avec lampes halognes est une technique
danalyse globale alors que la TIRM ultrasonore prsente un domaine dinvestigation plus
local.
Figure 1.28 TIRM ultrasonore (Gleiter et al. [89])
1.4.2.5 Thermographie InfraRouge Pulse (TIRP)
La TIRP, galement appele mthode flash associe la thermographie, est carac-

trise par lutilisation dun signal impulsionnel de type Dirac de la source sollicitante
(Prodjinonto [170], Sakagami and Kubo [180]). Selon lpaisseur et la diffusivit de la
pice inpecte, la dure de limpulsion varie de quelques millisecondes quelques secondes
(Ibarra-Castanedo [110], Larbi [122]). Le schma de cblage sapparente celui de la TIR,
la seule diffrence se situant au niveau de la source thermique, en loccurrence des lampes
flash qui ne ncessitent pas damplificateur de signal (figure 1.29).
Figure 1.29 Schma de cblage de la TIRP

Lutilisation de la lampe flash permet de solliciter la pice inspecte sur toutes ses fr-
quences caractristiques en une seule impulsion (Larbi [122]). De nombreux traitements
dimages (par exemple analyse de contraste, traitement statistique, dcomposition en va-
leurs singulires, cartographie de paramtres thermophysiques) sont appliqus aux images
thermiques ralises lors des contrles par TIRP dans le but de mettre en vidence les
dfauts prsents dans la pice inspecte (Larbi [122], Prodjinonto [170]). Les donnes de
sortie gnralement utilises sont des mesures de champs de diffusivit ou conductivit
(Parker et al. [158]) et les dfauts prsents au sein de la structure inspecte sont mis
en vidence partir des variations locales de ces proprits thermiques. Depuis quelques
annes, de multiples amliorations ont t apportes, rendant cette technique particulire-
ment efficace et rapide pour le CND en milieu industriel (Batsale et al. [30, 29], Bendada
et al. [33, 34], David et al. [66], Degiovanni and Laurent [69], Degiovanni [68], Ramond
et al. [172, 171]).
Cette mthode est plutt adapte ltude des matriaux fins (Prodjinonto [170])
(les temps dexposition requis pour des matriaux pais, de lordre de lheure, ne sont
pas compatibles avec un maintien uniforme du flash). Des essais en TIRP ont t rali-
ss sur des structures composites sandwichs avec une me en nid dabeille daluminium
et des structures monolithiques de type CFRP (plastique renforc de fibres de carbone
(Larbi [122])). La TIRP est galement employe dans le cadre du suivi dessai de fatigue
(Steinberger et al. [188]).
Il existe aussi une autre mthode appele thermographie infrarouge de phase pulse
qui est une combinaison des avantages de la thermographie pulse avec sa rapidit dac-
quisition de donnes, et de la thermographie module avec la localisation du dfaut en
profondeur moyennant lexcution de diffrents calculs aprs tests (Clerjaud [60], Guibert
[102], Ibarra-Castanedo and Maldague [113, 112], Ibarra-Castanedo [110]).
1.4.3 Bilan
Les mesures de champs thermiques sont des mesures optiques sans contact prsen-
tant de nombreuses potentialits pour le CND. Du point de vue industriel, on peut citer
comme avantages principaux leur rapidit dexcution, un caractre portable et une mise
en oeuvre facile ainsi que la capacit contrler les pices gomtrie complexe, de
manire globale et/ou locale selon le besoin. Certains inconvnients subsistent lheure
actuelle, parmi lesquels la difficult obtenir une information sur la profondeur du d-

faut, une dtection limite pour des pices composites paisses, une mesure de champs
par visualisation surfacique uniquement.
1.5 Conclusion
Ce premier chapitre a t ddi un tat de lart des techniques industrielles de CND
par mesure de champs, en focalisant notamment sur des applications composites. Trois
familles de mthodes ont t abordes :
les ultrasons, technique la plus utilise actuellement dans laronautique, qui conduit
des mesures de champs acoustiques,
la radiographie X et la tomographie X qui fournissent des mesures de champs den-
simtriques,
la thermographie infrarouge qui permet des mesures de champs thermiques.
Pour chacune delles, leur base physique et les principes gnraux de leur mise en oeuvre
ont t dtaills. Leurs avantages et les limitations ont t galement mis en vidence. Ceci
permet davoir une vision globale de ce qui est pratiqu actuellement dans lindustrie.
Si certaines des techniques voques sont couramment employes et matrises par

DAHER SOCATA, dautres en revanche suscitent son intrt. Lenjeu est alors de pou-
voir mesurer leurs performances respectives sur des configurations de matriaux/dfauts
spcifiques au secteur aronautique et de dvelopper si besoin de nouvelles comptences
en matire de CND. A cette fin, on se propose donc prsent de mener une tude com-
parative sur des mthodes reprsentatives de chacune des trois familles prcites.
Chapitre 2
Matriaux et protocoles exprimentaux
Sommaire
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.2 Prsentation des chantillons . . . . . . . . . . . . . . . . . . . 42
2.2.1 Eprouvette 1 : matriau monolithique avec inserts Tflon . . 43
2.2.2 Eprouvette 2 : matriau monolithique avec trous fond plat . . 44
2.2.3 Eprouvette 3 : matriau sandwich avec inserts Tflon (peaux
de 0.66 mm) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
de 1.08 mm) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.2.5 Eprouvette 5 : matriau multi-structures . . . . . . . . . . . . . 50
2.3 Protocoles exprimentaux . . . . . . . . . . . . . . . . . . . . . 54
2.3.1 Mesures de champs acoustiques : les ultrasons . . . . . . . . . . 54
2.3.2 Mesures de champs densimtriques . . . . . . . . . . . . . . . . 61
2.3.3 Mesures de champs thermiques : thermographie infrarouge . . . 66
2.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
42 2. Matriaux et protocoles exprimentaux
2.1 Introduction
Cette partie est consacre la mise en place dun certain nombre de techniques de CND
parmi les trois familles de mthodes de mesure de champs voques au chapitre prcdent.
La dmarche comparative a t mene sur des structures reprsentatives des constituants
dun fuselage composite, savoir des structures monolithiques et des structures sandwichs.
Chacune des mthodes employes a t teste sur cinq prouvettes talons :

deux prouvettes monolithiques et deux prouvettes sandwichs fabriques au Centre
de Ressources et de Comptences Composites de Tarbes,
une prouvette multi-structure comprenant deux zones monolithiques et une zone
sandwich, fabrique par la socit NIMITECH situe Bagnres-de-Bigorre.
Ces talons sont constitus de deux typologies de dfauts artificiels visant reprsenter
des porosits et/ou dlaminages :
des inserts Tflon insrs lors du process de fabrication,
des trous fond plat excuts aprs fabrication.
Ces dfauts rpondent aux exigences des normes AITM AITM6-0013 [5], AITM6-0016
[7], AITM6-4002 [8], AITM6-4005 [9], AITM6-4010 [10] rgissant notamment le contrle
ultrasonore, et sont frquemment employs dans les tudes de CND sur structures com-
posites (Bai and Wong [25], Batsale et al. [29], Bremond [47], Choi et al. [56], Ibarra-
Castanedo et al. [111], Junyan et al. [115], Larbi [122], Meola et al. [138, 139, 140], Montanini
[148], Rodriguez et al. [177], Sakagami and Kubo [180]).
2.2 Prsentation des chantillons

Chaque talon fabriqu est compos de dfauts de diffrents diamtres diffrentes
profondeurs. Parmi les tailles employes, lensemble des chantillons comporte des inserts
Tflon ou des trous fond plat de 6 mm de diamtre, taille de dfaut rfrence dans
les standards aronautiques. Nanmoins, dautres tailles (de 15 mm 3 mm de diamtre)
sont utilises de manire mieux apprhender les limites de dtection des techniques
employes.
Les diffrents procds de fabrication utiliss permettent lobtention dtats de surface

distincts. Ils ont t rpertoris de la manire suivante :
type A : aspect lisse et brillant,
type B : aspect rugueux, bossel et mat,
type C : aspect rugueux et mat,
2.2. Prsentation des chantillons 43
type D : aspect lisse et mat,

type E : aspect rugueux et brillant.
Ils seront expliqus dans la suite du chapitre.
2.2.1 Eprouvette 1 : matriau monolithique avec inserts Tflon

Ce matriau est une structure monolithique stratifie compose de 8 plis de tissus
carbone Hexcel HexPly G803(tissu)/1457(rsine poxy). Chaque pli est orient 30 et
prsente une paisseur thorique de 0.35 mm (figure 2.1).
Figure 2.1 Photo du matriau monolithique avec inserts Tflon
Cette prouvette a t ralise partir de fibres pr-imprgnes de rsine par moulage

poinon/matrice ou moulage sous vide (dpression/autoclave). Il sagit du procd de
fabrication actuellement le plus rpandu pour les composites dans laronautique [90].
Le protocole exprimental de fabrication est compos de plusieurs tapes :

1. choix du procd/contre moule : la pice peut tre ralise avec un moule et contre-
moule rigides (pressage chaud) ou avec un moule rigide et un contre-moule souple
(dpression/autoclave) ; concernant cette prouvette, un moule et un contre-moule
rigides ont t employs, conduisant un aspect lisse et brillant pour les deux
surfaces de la pice (type A dtat de surface),
2. prparation du moule et du contre-moule (dmoulant),
3. empilement des tissus selon la squence dsire,
4. mise en place du contre-moule,
5. pressage chaud.
Au cours de ltape 4, des disques en Tflon de diffrents diamtres (4, 6, 8, 10, 12
et 15 mm) ont t mis en place diffrentes profondeurs (0.7, 1.4 et 2.1 mm) (figure
2.2) ; prcisment, chaque range horizontale dinserts a t positionne tous les deux plis
(figure 2.3).
Figure 2.2 Vue de dessus du matriau monolithique avec inserts Tflon - Position-
nement des inserts
Figure 2.3 Matriau monolithique avec inserts Tflon - Coupe schmatique dans
lpaisseur au niveau dune range dinserts de mme diamtre
2.2.2 Eprouvette 2 : matriau monolithique avec trous fond plat

Ce matriau stratifi est compos de tissus en carbone unidirectionnels primprgns
Hexcel HexPly T700 (fibres Toray )/M21 (rsine poxy). Cette pice comprend une
variation dpaisseur avec une partie 4 mm sur 50 mm de long et une partie 8 mm

sur 130 mm de long (figures 2.4 et 2.5). La squence dempilement de la partie 4 mm
dpaisseur correspond [0, 45, 90, -45]2,S et celle de la partie 8 mm dpaisseur est
[0, 45, 90, -45]4,S .
Figure 2.4 Photo du matriau monolithique avec trous fond plat
Le procd de fabrication suit globalement la mme dmarche que prcdemment.

Cependant, ltape 1 consiste ici utiliser un moule rigide et un contre-moule souple
(emploi dune bche) sans tissu darrachage, ce qui implique dans ce dernier cas lobtention
dun tat de surface lisse et brillant pour les deux cts de la pice (type A). Ltape 5
est, quant elle, dcompose en deux parties avec dabord une mise sous vide puis une
mise en autoclave.
Les dfauts artificiels sont ici des trous fond plat de diffrents diamtres (2, 3, 6 et
10 mm) jusqu diffrentes profondeurs (0.52, 2.8, 3, 3.48 et 5.2 mm) aprs fabrication de
lprouvette (figures 2.5 et 2.6).
Figure 2.5 Vue de dessus du matriau monolithique avec trous fond plat - Position-
nement des trous
Figure 2.6 Matriau monolithique avec trous fond plat - Coupe schmatique dans
lpaisseur au niveau dune range de trous de mme diamtre

de 0.66 mm)
Le matriau sandwich faible paisseur de peau comprend une me en mousse haute
performance Airex R82 fabriqu par Alcan Composites de 10 mm dpaisseur et de
deux peaux composes de 3 plis de tissus en carbone-verre primprgns Hexcel HexPly
G973(tissu)/913(rsine poxy). Lorientation des plis est 0 pour une paisseur thorique
de 0.22 mm (figure 2.7).
Figure 2.7 Photo du matriau sandwich avec inserts Tflon (peaux de 0.66 mm)
Le procd de fabrication employ sarticule autour des tapes suivantes :

1. prparation du moule rigide ; la face de ce ct du matriau est donc lisse et brillante
(type A, figure 2.9),
2. empilement des tissus selon la squence dsire,
3. mise en place du tissu darrachage (impliquant un aspect rugueux),
4. mise en place du feutre de drainage,
5. mise en place dune bche comme contre-moule souple,

6. compactage des tissus par tirage de vide (dpression) et mise en tuve.
Lors de lempilement des tissus (tape 2), des inserts Tflon de 3 et 6 mm de diamtre
ont t positionns entre tous les plis des peaux ainsi quentre la mousse et le premier
pli de chaque peau (figure 2.9). Du fait de la fine paisseur des plis utiliss, deux inserts
ont t superposs, chaque localisation de dfaut artificiel, afin de simuler au mieux
un dlaminage (figure 2.9). Cette double paisseur, de lordre de grandeur de lpaisseur
des plis, a alors gnr des surpaisseurs sur la surface finale de la pice du ct du
contre-moule souple, qui prsente de fait un aspect rugueux, bossel et mat (type B). Par
consquent, cette solution na pas t retenue par la suite.
Figure 2.8 Vue de dessus du matriau sandwich avec inserts Tflon (peaux de 0.66
mm) - Positionnement des inserts Tflon
Figure 2.9 Matriau sandwich avec inserts Tflon (peaux de 0.66 mm) - Coupe
schmatique au niveau dune range dinserts de mme diamtre

de 1.08 mm)
Cette dernire mono-structure est un matriau sandwich prsentant une me en mousse
haute performance Airex R82 de 15 mm dpaisseur et des peaux constitues de 3 plis
de tissu carbone Hexcel HexForce G0926. La squence dempilement de chaque peau est
[45/0/45] pour une paisseur thorique de 1.08 mm (figure 2.10). Cette prouvette a
t fabrique par infusion de rsine liquide (procd LRI, Liquid Resin Infusion), procd
pour lequel les prformes sches sont imprgnes par la rsine poxy SR1710/SD8822 dans
lassemblage moule/contre-moule avant polymrisation thermique.
Figure 2.10 Photo du matriau sandwich avec inserts Tflon (peaux de 1.08 mm)
Aprs la prparation de la mousse (chanfreinage des bords, perage au laser tous les
5 mm avec des trous de 1 mm de diamtre pour conduire la rsine), un point important
de ce procd consiste prparer la prforme sche selon les diffrentes tapes suivantes :
1. drapage de tissus carbone secs poudrs (peau infrieure), positionnement de lme
en mousse, drapage de tissus carbone secs poudrs (peau suprieure) ; dans le cas
prsent, le drapage a t ralis sur une plaque de verre jouant le rle de moule
rigide plan, amenant de ce ct un aspect de type A,
2. mise en place du tissu darrachage,
3. mise en place du feutre de drainage,
4. implantation du rseau de vide et mise en place de la bche qui sert de contre-moule
souple,
5. prformage des tissus secs poudrs par tirage de vide (dpression) et mise en tuve
(dure 1h 105C avec une rampe de 1C/min).
Pour ltape dinfusion, il convient en premier lieu denlever au montage prcdent

le feutre de drainage et dajouter le rseau des valves dalimentation rsine qui dpend
de la taille de lprouvette. Dans le cas prsent, la rsine choisie permet cette opration
temprature ambiante et donc le suivi de son volution au cours du temps. Une fois
linfusion termine, lensemble est mis en tuve (4h 80C) afin de polymriser la rsine.
Durant cette tape de polymrisation, le rseau de vide est maintenu tandis que celui
dinfusion est clamp.
Durant le drapage, des inserts Tflon de 3 et 6 mm de diamtre ont t insrs entre
chaque pli et entre la mousse et le premier pli (figures 2.11 et 2.12). Contrairement
lprouvette sandwich prcdente, un seul insert a t positionn chaque localisation de
dfauts afin de minimiser la dformation des tissus et doptimiser la planit de ltat de
la surface. Compte tenu de ce protocole, ltat de surface du ct du contre-moule est
rugueux et mat (type C).
Figure 2.11 Vue de dessus du matriau sandwich avec inserts Tflon (peaux de 1.08
mm)
Figure 2.12 Matriau sandwich avec inserts Tflon (peaux de 1.08 mm) - Coupe
schmatique dans lpaisseur au niveau dune range dinserts de mme diamtre
2.2.5 Eprouvette 5 : matriau multi-structures

Ce matriau multi-structures est compose de deux structures monolithiques strati-
fies dpaisseurs diffrentes et dune structure sandwich avec une me en mousse haute
performance Airex R82 de 15 mm dpaisseur (figure 2.13).
Figure 2.13 Photo du matriau multi-structures
Tous les plis de cette structure sont raliss partir de tissus carbone Hexcel HexForce
G0926. La plus mince des parties monolithiques est compose de 6 plis avec la squence
dempilement [45/0/45]2 , la plus paisse comporte 12 plis empils selon la squence
[(45/0/45) (0)6 (45/0/45)]. Les peaux de la partie sandwich sont quant elles com-
poses de 3 plis de ce mme tissu avec la squence dempilement [45/0/45].
Le procd dinfusion LRI a galement t employ pour fabriquer cette prouvette,

selon un protocole exprimental quasi-identique celui de la pice prcdente. Nanmoins,
tant donne sa gomtrie complexe, certaines particularits ont t mises en place, que
lon se propose de dtailler prsent.
Tout dabord, la stratgie choisie ici est lutilisation dun moule et dun contre-moule
rigides : le moule rigide permettra de donner sa forme gomtrique la structure tandis
que le contre-moule rigide (plaque de PMMA) donnera la surface plane de la pice. Pour
raliser le moule rigide, il a fallu au pralable raliser un master en rsine usinable (figure
2.14) sur lequel le moule rigide a t ralis.
Figure 2.14 Master en rsine usinable
Le moule rigide est ralis en deux parties :

la premire partie est du gel coat appliqu directement sur le master permettant
dobtenir, dune part la forme gomtrique de la pice, et dautre part, lors de la
fabrication de la pice, un tat de surface lisse et mat ;
la deuxime partie est un empilement de tissus de verre (moulage au contact, rsine
poxy) directement drapps sur la couche de gel coat ayant atteint ltat B (figure
2.15).
Figure 2.15 Moule rigide en verre-poxy

A partir de ce moule rigide, on ralise ltape de prformage des tissus secs poudrs
de la mme manire que pour lprouvette 4.Durant ltape de prformage, des inserts
Tflon de 3 et 6 mm de diamtre (figure 2.16) ont t positionns tous les deux plis
dans les parties monolithiques et tous les plis dans la structure sandwich, selon la pro-
cdure dtaille prcdemment (figure 2.17). Diffrents prformages ont t raliss afin
doptimiser dune part la forme de la pice avec le respect des angles et des pentes au
niveau de la structure sandwich et, dautre part, de maximiser la rticulation des tissus
poudrs sur lesquels les inserts Tflon ont t positionns. Pour cela, les trois premiers
et trois derniers plis ont t prforms un par un ; concernant le bloc de six plis 0,
ceux-ci ont t prforms tous les deux plis.
Figure 2.16 Vue de dessus du matriau multi-structures - Positionnement des inserts

Tflon
(a)
(b)
(c)
Figure 2.17 Matriau multi-structures - Coupes schmatiques dans lpaisseur au ni-

veau dune range dinserts de mme diamtre : (a) partie monolithique fine, (b) zone
sandwich, (c) partie paisse
En ce qui concerne ltape dinfusion, un film micro perfor a t dispos entre la pice
et le tissu darrachage de faon limiter laspect rugueux obtenu par le contact direct de
la pice avec le tissu darrachage (figure 2.18). De plus, le contre-moule rigide en PMMA
a t usin aux dimensions de lprouvette et perfor selon le positionnement du rseau
dinjection. Une fois linfusion termine, la structure a t mise en tuve (16h 60C)
afin de polymriser le mlange rsine/durcisseur inject (rsine poxy SR1710/SD8822) et
ensuite dmoule. Lensemble du procd de fabrication (prformage et infusion) implique
lobtention du ct moule dun tat de surface lisse et mat du fait du gel coat (type D)
alors que le ct contre-moule est rugueux et brillant grce au film micro perfor (type E,
figure 2.17).
Figure 2.18 Mise en place du film micro perfor

2.3 Protocoles exprimentaux

On se propose dans cette partie de dtailler les diffrentes tapes de mise en oeuvre
des techniques employes en vue de ltude comparative. A la fin de la prsentation de
chacune dentre elles, une illustration de la mesure de champ obtenue sera prsente.
2.3.1 Mesures de champs acoustiques : les ultrasons
Lensemble des contrles ultrasonores ont t raliss au laboratoire IUT-ICA de

Tarbes. Les trois techniques mises en jeu sont lcho ultrasonore (par contact et par
immersion) et la transmission ultrasonore par immersion.
En ce qui concerne le premier aspect de la dmarche dtude, savoir lmission des

ondes acoustiques, il convient de prciser quelles sont dans le cas prsent du type volu-
mique longitudinale. Pour ce qui est de leur rception, le flux sera dtect en configuration
face avant pour lcho ultrasonore et en configuration face arrire pour la transmission
ultrasonore par immersion. Lensemble des contrles ultrasonores ont t effectus avec
des traducteurs de frquence 5 MHz, et ce de faon respecter les usages aronautiques
actuels sur chaque talon fabriqu.
2.3.1.1 Echo ultrasonore par contact : contrle en face avant
Le dispositif exprimental (figure 2.20) est compos des lments suivants :

un traducteur multi-lments 5L64-NW1 (Olympus) quipe dun sabot ; il sagit
dun traducteur ultrasonore de 5 MHz comprenant 64 lments (figure 2.19),
un Omniscan 32 : 128 PR (Olympus) ; cet appareil de contrle ultrasonore ddi
lutilisation des traducteurs haute frquence gre dune part lmission et la d-
tection du flux dondes et, dautre part, ltalonnage du test, lanalyse ainsi que la
visualisation des rsultats,
un Freescan (IXTREM) ; ce systme de codage, sans contact par tracking optique,
calcule les mouvements du traducteur chaque instant grce la localisation de la
pastille argente positionne sur le traducteur (figure 2.19) ; ceci permet de raliser
les cartographies 2D par projection des dfauts en surface.
2.3. Protocoles exprimentaux 55
Figure 2.19 Capteur multi-lments quip du sabot
Figure 2.20 Dispositif exprimental du contrle par cho ultrasonore par contact
La mthode dcho ultrasonore par contact tant manuelle, le systme de fixation de la

structure inspecte doit tre stable, parallle aux capteurs optiques du Freescan et centr
par rapport son champ optique (figures 2.20 et 2.21).
Figure 2.21 Systme de fixation de la structure inspecte (matriau multi-structures)
Une fois ce dispositif exprimental ajust, la structure est positionne sur celui-ci
et recouverte de lagent couplant (constitu ici de 50% de gel couplant et 50% deau).
Ltalonnage du capteur par rapport la structure inspecte est ralis avec lOmniscan
et se dcompose en trois tapes :
1. calcul de la vitesse de propagation de londe dans le matriau inspect : ceci permet
de traduire la mesure en temps de vol en mesure en paisseur, et ainsi den dduire
la profondeur des dfauts en fonction des chos mesurs,
2. positionnement des trois portes de dtection des chos correspondant aux trois types
dchos dtaills dans le chapitre 1 :
la porte I (jaune sur le A-scan de la figure 2.22) place sur lcho dentre et
synchronise sur limpulsion de londe ultrasonore,
la porte A (rouge) positionne sur lcho de fond et synchronise sur le pic dentre
de londe dans la pice,
la porte B (verte) place entre lcho dentre et lcho de fond afin de distinguer
lapparition dun dfaut et synchronise comme la porte A sur le pic dentre de
londe dans la pice.
Ainsi la distance entre lcho dentre (dtect par la porte I) et lcho de fond
(dtect par la porte A) correspond lpaisseur de la zone inspecte, la porte B
localisant la position du dfaut. La hauteur de chaque porte doit tre au minimum
suprieure celle du bruit de mesure,
3. rglage de la sensibilit de dtection de faon obtenir le mme pouvoir de dtection

quelle que soit la loi focale utilise ; ici lcho de fond a t ajust 80% de son
amplitude.
Figure 2.22 A-scan ralis en cho ultrasonore par contact (phase dtalonnage, sans
dfaut)
Les cartographies 2D obtenues aprs contrle sont de deux types :

un C-Scan en amplitude reprsentant la rpartition des amplitudes des chos (figure
2.23-a),
un C-Scan en paisseur permettant de localiser les dfauts dans lpaisseur de la
structure (figure 2.23-b).
(a)
(b)
Figure 2.23 C-Scan (matriau monolithique avec insert Tflon ) : (a) en amplitude,
(b) en paisseur
2.3.1.2 Echo ultrasonore par immersion : contrle en face avant
Le principe de contrle de lcho ultrasonore par immersion est quasiment identique

celui de lcho ultrasonore par contact. Dans cette mthode, la structure inspecte et le
traducteur sont compltement immergs dans une cuve, optimisant ainsi la couche dagent
couplant (figure 2.24).
Figure 2.24 Pice et capteur immergs dans la cuve
Les lments suivants composent le dispositif exprimental (figure 2.25) :

un traducteur multi-lments 5L64-NW1 (Olympus),
un Omniscan 32 : 128 PR (Olympus),
une cuve munie de bras automatiques permettant le positionnement des capteurs
en X et en Y,
un ordinateur avec le logiciel UT-WIN permettant de piloter les bras automatiques
en X et en Y.
Figure 2.25 Dispositif exprimental de lcho ultrasonore par immersion
En comparaison avec la mthode prcdente, la structure inspecte est simplement

pose sur deux cales positionnes sur le rflecteur dondes en PMMA (figure 2.26). La
procdure dtalonnage du capteur par rapport la pice inspecte diffre quant elle
simplement au niveau de la synchronisation de la porte I. Celle-ci doit tre synchronise
sur le pic dentre de londe ultrasonore dans la pice de faon saffranchir dventuels
problmes de paralllisme entre lprouvette et les axes XY.
Figure 2.26 Dispositif de fixation de lprouvette
Une fois lprouvette place et ltalonnage termin, la sonde multi-lments doit tre
positionne au dessus dun coin de la pice contrle qui constituera le dpart de la
trajectoire du contrle. La distance entre le capteur et la pice (distance de focalisation)
est prise de faon sextraire du champ proche du faisceau du traducteur, soit 50 mm
dans le cas prsent. Les cartographies 2D obtenues aprs le contrle sont les mmes que
pour lcho ultrasonore par contact.
2.3.1.3 Transmission ultrasonore par immersion : contrle en face arrire
La transmission ultrasonore par immersion suit les mmes principes dinspection que
lcho ultrasonore par immersion mais selon un contrle en configuration face arrire.
Deux traducteurs mono-lments V309 (frquence de 5 MHz, diamtre de 0,5 pouce) sont
utiliss, lun jouant le rle dmetteur et lautre de rcepteur (figure 2.27).
Figure 2.27 Montage de transmission ultrasonore par immersion
Le support de fixation de lprouvette a pour objectif de la maintenir en position

verticale (figure 2.27). Concernant la procdure dtalonnage, ltape de positionnement
des portes ultrasonores diffre des deux prcdentes techniques. Le principe de contrle
par transmission consiste en effet dtecter labsence ou lattnuation de londe mise
par le premier capteur. De ce fait, le A-scan ne prsente ici que lcho de rception ; seule
donc la porte A doit tre positionne sur cet cho (figure 2.28). La hauteur de cette porte
est dfinie au-dessus du bruit de mesure comme prcdemment. Lorsque lamplitude de
lcho de rception prsente une baisse 50% par rapport ltalonnage ralis sur une
zone sans dfaut, on considre quil y a prsence dun dfaut (mthode -6 dB). La seule
cartographie 2D obtenue avec cette mthode est un C-Scan en amplitude (figure 2.29).
Figure 2.28 Positionnement de la porte A
Figure 2.29 C-Scan en amplitude (sur une partie du matriau multi-structures)
2.3.2 Mesures de champs densimtriques

2.3.2.1 Radiographie X
Tous les contrles de radiographie X ont t raliss par le site de Tarbes du groupe
Daher-Socata. Le montage exprimental utilis est compos du tube rayons X (met-
teur), du film radiographique (rcepteur) ainsi que de lcran de blocage (figures 2.30 et
2.31). Le pupitre de commande permet dajuster lintensit, la tension ainsi que la position
du tube rayons X ainsi que de dclencher le test.
Figure 2.30 Montage exprimental du contrle par radiographie X

Figure 2.31 Dispositif exprimental ddi au contrle par radiographie X utilis par
Daher Socata sur le site de Tarbes
Avant de dmarrer lanalyse, il convient de rgler les deux distances spcifiques la

configuration adopte, savoir la distance source-film (DSF), gale dans notre cas 950
mm, et la distance pice-film (DPF), ici nulle (la pice est pose sur le film). Ensuite,
lintensit et la tension du tube rayons X doivent tre ajustes de faon obtenir un
faisceau de rayons X suffisamment puissant pour optimiser le contrle. En loccurrence,
lintensit du tube retenue est de 19 mA, la tension employe de 20 kV pour les structures
monolithiques et 10 kV pour les structures sandwichs. Une fois le contrle termin, les
ngatifs du film sont ensuite dvelopps pour obtenir la visualisation des rsultats (figure
2.32).
Figure 2.32 Radiographie X (matriau sandwich avec peau de 0.66 mm)
2.3.2.2 Tomographie X
Les contrles par tomographie X ont t pris en charge par la socit Tomo Adour
situe Pau. Un scanner mdical de quatrime gnration adapt la mesure de pices
industrielles a t employ, savoir le SOMATOM Sensation Open de SIEMENS dnergie
maximale 140 keV (figure 2.33). Cet appareil se compose :
dun ordinateur de gestion permettant dexcuter les programmes danimation, de
contrle et de dialogue avec lutilisateur,
dun systme dacquisition assurant la saisie et le prtraitement des mesures datt-
nuation ; celui-ci est constitu du banc mcanique de positionnement et de balayage
de lobjet, du gnrateur de rayonnement X (source en ventail) et de lensemble
de dtection (capteurs cramiques haute dfinition en ligne, figure 2.34) associ
llectronique de numrisation,
du systme danalyse des images et de reconstruction fournissant notamment les
coupes en temps rel et les vues 3D aprs post-traitement.
Du fait dune large couverture permise par la configuration du couple source-dtecteurs
ainsi que par son temps de rotation rapide (environ une seconde), ce scanner permet
lanalyse de volumes importants dans un temps trs court.
Figure 2.33 Scanner SOMATOM Sensation Open
Figure 2.34 Schma de fonctionnement du scanner mdical SIEMENS
La prparation des mesures sarticule autour des points importants suivants :

mise en place et rglage de la position des prouvettes : les pices sont fixes sur un
support en bois dont la faible densit minimise les perturbations avec de ladhsif
double face pour viter tout mouvement lors du dplacement du support mobile ;
les pices sont ensuite positionnes prcisment laide de lasers de faon dfinir
lorigine de lexamen ;
rglage du banc mcanique sur lequel reposent les contraintes spciques du contrle
effectuer : ses paramtres (longueur de la translation, poids,..) dpendent directe-
ment du type et de la nature des pices examiner ;
confinement de la pice o se trouve la machine pour viter toute radiation des

oprateurs,
excution automatique du protocole de pr-chauffe de la machine comprenant no-
tamment ltalonnage,
dfinition de la rsolution spatiale, soit ici une taille de voxel de 0,1 x 0,1 x 0,1
mm3 (en dautres termes, les coupes 2D dacquisition prsentent une paisseur de
0.1 mm).
Les coupes radiographiques sont ralises perpendiculairement la surface moyenne

des pices. A partir de ces rsultats, on peut galement accder des cartographies 2D
dites reconstruites (reslice), cest--dire situes dans des axes diffrents de ceux de la
mesure initiale, notamment dans des plans parallles la surface moyenne des pices
(figure 2.35-a). Les algorithmes de reconstruction permettent, galement, de se ramener
au volume de la pice inspecte (figure 2.35-b).
(a)
(b)
Figure 2.35 Reconstructions par tomographie X (matriau sandwich avec peaux de

0.66 mm) : (a) Coupe 2D reconstruite paralllement la surface moyenne de la pice,
(b) : reconstruction volumique
2.3.3 Mesures de champs thermiques : thermographie infrarouge
Trois techniques de thermographie infrarouge ont t tudies :

les thermographies classique (TIR) et module (TIRM) au laboratoire ENIT-LGP
Tarbes,
la thermographie pulse (TIRP) par lentreprise THERMOCONCEPT situe Ta-
lence et associe au laboratoire TREFLE de lEcole Nationale des Arts et Mtiers
de Bordeaux.
2.3.3.1 Thermographie infrarouge (TIR) : contrle en face avant et face ar-

rire
Pour la thermographie infrarouge classique, le matriel utilis est le suivant :

1 ou 2 lampes halognes 120V 1000W,
un amplificateur de signal (Pulsar),
une camra infrarouge JADE CEDIP rtroffite en FLIR Titanium (rsolution : 20
mK, objectif : MW 50 mm 2.0 Jade) qui joue le rle de dtecteur,
le logiciel ALTAIR (FLIR) permettant la ralisation du film ainsi que sa visualisa-
tion.
Lamplificateur de signal permet de piloter le voltage des lampes halognes faisant varier
la puissance du flux thermique mis par ces dernires. Le choix de cette tension dpend de
la nature du matriau inspect, de son paisseur, de sa conductivit, de son tat de surface
et de la configuration de test adopte. Gnralement, les valeurs retenues sont comprises
entre 90 et 110 V (niveau maximal dutilisation). Si un voltage suprieur savre ncessaire,
on associe alors deux lampes halognes.
Le dispositif exprimental utilis en face avant, pour lequel le flux thermique est mis
et dtect du mme ct de la pice, est illustr la figure 2.36. La configuration face
arrire qui consiste dtecter le flux thermique du ct oppos son mission est quant
elle prsente sur la figure 2.37. Dans ces deux configurations, le positionnement des
metteurs et rcepteurs doit en premier lieu tre ajust afin de minimiser les bruits de
mesure. En particulier, la lampe halogne et lobjectif de la camra doivent tre parallles
la pice et leur point focal respectif (thermique pour lhalogne, optique pour la camra)
doit tre centr sur la zone inspecte (des lasermtres sont utiliss cette fin, figure 2.37).
Figure 2.36 Dispositif exprimental en configuration face avant
Figure 2.37 Dispositif exprimental en configuration face arrire

Cinq tapes doivent ensuite tre ralises avant de procder aux essais :
ltalonnage sur une rfrence (corps noir) ; cette tape seffectue de faon automa-
tique avec la camra employe qui bnficie dun systme de dernire gnration,
la dfinition de lmissivit du matriau tudi ; cette donne est trs sensible
de nombreux paramtres intrinsques mais aussi extrinsques au matriau tudi
(Bodnar [39]) ; en labsence dinformations spciques, nous avons choisi dentrer ici
une valeur constante et uniforme, savoir celle du carbone (0,98),
le rglage de la frquence dacquisition du champ de temprature ; dans le cas pr-
sent, nous avons considr une valeur de 50 Hz, soit 50 images par seconde,
la mesure de la distance prouvette-camra afin dobtenir la rsolution spatiale ;
dans notre cas, cette distance est comprise suivant les essais entre 80 cm et 1 m,
la mise au point de lobjectif.
Une fois ces diffrents rglages effectus, le test se droule selon le protocole suivant :
1. extinction des lumires, la salle de contrle devant tre compltement dans le noir
de faon optimiser le test,
2. simultanment, dmarrage du film sur le logiciel Altair et allumage de la lampe

halogne,
3. chauffage de la pice pendant 10 secondes en moyenne,
4. extinction de la lampe halogne,
5. attente de 20 secondes pour le refroidissement de la pice,
6. arrt du film.
Comme dtaill dans le chapitre 1, le rsultat du test correspond un film de la dure du

contrle (environ 30 secondes) reprsentant lvolution temporelle du champ thermique
au cours du temps. Aprs la phase dacquisition, un traitement des donnes peut savrer
ncessaire. De faon liminer au maximum les effets perturbants pour la mesure (notam-
ment limpact de la temprature ambiante) et amliorer linterprtation des rsultats,
nous avons opt pour une analyse en relatif des donnes en soustrayant chaque instant
limage thermique du point de dpart de lessai. Les rsultats traiteront donc exclusive-
ment des champs de variations de temprature T = T T0 avec T la temprature
linstant t et T0 la temprature initiale, sous la forme de cartographies 2D (figure 2.38).
Figure 2.38 Mesures de champ thermique (en relatif) ralises par TIR (matriau
monolithique avec inserts Tflon )
2.3.3.2 Thermographie infrarouge module (TIRM) : contrle en face avant

et face arrire
En comparaison avec le matriel employ pour la TIR, cette mthode requiert en plus
un gnrateur de signal (HAMEG) qui permet dune part denvoyer un signal sinusodal
via lamplificateur la lampe halogne et, dautre part, de le connecter avec la camra
infrarouge. Le logiciel employ pour la ralisation des films est par ailleurs dans ce cas
ALTAIR-LI (FLIR), capable de reconnatre la frquence de travail du gnrateur ainsi
que le type de signal mis.
Les configurations employes pour cette technique ainsi que les prcautions de posi-
tionnement de la source et de la camra sont les mmes que celles indiques prcdemment.
Le protocole exprimental diffre lgrement :
1. extinction des lumires,

2. simultanment, dmarrage du film sur le logiciel ALTAIR-LI et lancement du signal
sinusodal,
3. attente de larrt simultan du film et du flux thermique.
La dure dun film de TIRM peut varier entre 1 min 20 sec et 2 min 40 sec suivant
lpaisseur des prouvettes testes. Ce film permet daccder des cartographies 2D de
lvolution de lamplitude du champ thermique (figure 2.39-a), et de lvolution du d-

phasage entre lmission et la rception (figure 2.39-b). On notera que lidentification des
dfauts a t faite partir des cartographies de phase qui savrent nettement moins
bruites.
(a)
(b)
Figure 2.39 Cartographie 2D des mesures de champs thermiques obtenus par TIRM
(matriau monolithique avec inserts Tflon ) : (a) amplitude, (b) phase
2.3.3.3 Thermographie infrarouge pulse (TIRP) : contrle en face avant et

face arrire
Cette technique a t mise en oeuvre laide du matriel suivant :

2 lampes flash 230V 1500J,

une camra infrarouge SC7500 (rsolution : 20 mK) (FLIR),
un boitier de dclenchement (THERMOCONCEPT) destin dune part synchro-
niser le lancement du test avec celui du flash de chaque lampe employe et, dautre
part, de trigger la camra infrarouge avec les deux sources thermique utilises,
le logiciel ALTAIR (FLIR), pour la ralisation du film,
le logiciel THERMO-CND (THERMOCONCEPT) qui traite les donnes de manire
fournir une seule image reprsentative de lensemble du film effectu.
La socit THERMOCONCEPT utilise un dispositif exprimental spcifique pour les

configurations en face avant ou arrire (figure 2.40 et 2.41). Dans ces deux cas, le banc
exprimental CND permet le dplacement vertical de la camra, le dplacement horizontal
et vertical des systmes dexcitation ainsi que le positionnement de lchantillon. De plus,
un rflecteur est positionn lextrmit de chaque lampe flash de faon focaliser et
concentrer lmission du flux thermique sur la zone inspecte.
Figure 2.40 Dispositif exprimental de la TIRP en configuration face avant

Figure 2.41 Dispositif exprimental de la TIRP en configuration face arrire
Le protocole exprimental excut est le suivant :
1. extinction des lumires,
2. lancement du film et du flash des lampes par dclenchement synchronis,
3. attente de larrt du film.
La dure dun test varie de 0.18 11 sec suivant lpaisseur des talons contrls. Aprs
traitement du film par le logiciel THERMO-CND, le rsultat obtenu est une unique car-
tographie 2D caractristique des mesures de champs thermiques (figure 2.42).
2.4. Conclusion 73
Figure 2.42 Cartographie 2D des mesures de champs thermiques par TIRP (matriau
monolithique avec inserts Tflon )
2.4 Conclusion
Il tait important dans cette partie de bien dtailler tant la conception et le procd
de fabrication des matriaux et structures faisant lobjet de ltude que les dtails de
mise en oeuvre des huit techniques de CND employes dans chacune des familles de
mesures de champs (trois ultrasonores, deux densimtriques, trois thermiques). Toutes
ces informations seront en effet cruciales pour la comprhension et linterprtation des
rsultats qui vont tre prsents au chapitre suivant.
Chapitre 3
Etude comparative sur prouvettes

talons
Sommaire
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
3.2 Principes danalyse . . . . . . . . . . . . . . . . . . . . . . . . . 77
3.2.1 Dimensionnement des dfauts par traitement dimage . . . . . . 77
3.2.2 Critres danalyse . . . . . . . . . . . . . . . . . . . . . . . . . . 79
3.3 Rsultats exprimentaux . . . . . . . . . . . . . . . . . . . . . . 80
3.3.1 Etalons monolithiques avec inserts Tflon . . . . . . . . . . . 80
3.3.2 Etalon monolithique avec trous fond plat . . . . . . . . . . . 88
3.3.3 Etalons sandwich avec inserts Tflon . . . . . . . . . . . . . . 94
3.3.4 Synthse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
3.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
76 3. Etude comparative sur prouvettes talons
3.1 Introduction
Lenjeu de ce chapitre est de prsenter la dmarche danalyse comparative de tech-
niques de CND. Si certaines tudes dans ce sens existent dans la littrature, elles ne portent
en gnral que sur les variantes de mesures de champs thermiques abordes au chapitre
prcdent (voir par exemple Clerjaud [60], Giorleo and Meola [88], Ibarra-Castanedo et al.
[111], Larbi [122], Pradere et al. [169], Sakagami and Kubo [180]). Loriginalit de ce tra-
vail consiste en une large porte de cette comparaison qui couvre les trois grandes familles
de mthodes de mesure de champs au travers de huit techniques spcifiques mises en
oeuvre sur les cinq prouvettes talons prsentes au chapitre prcdent :
champs acoustiques avec lcho ultrasonore par contact, lcho ultrasonore par im-
mersion et la transmission ultrasonore par immersion,
champs densimtriques avec la radiographie X et la tomographie X,
champs thermiques avec la thermographie Infra-Rouge (TIR), la TIR Module
(TIRM) et la TIR Pulse (TIRP).
A ce stade, lprouvette multi-structure permet une analyse uniformise en terme de
matriaux et de dfauts tandis que les autres prouvettes talons prsentent des dfauts
de diffrents types (inserts ou trous fond plat) et des diamtres plus varis.
Ces mthodes ont gnralement des applications et des notorits diffrentes dans
le domaine du CND sur structures composites pour lindustrie aronautique. DAHER-
SOCATA ayant la volont dlargir ses comptences en la matire, lobjectif principal de
cette tude comparative est de fournir un bilan des capacits de dtection de chacune
des techniques employes en terme de largeur de spectre et de prcision. Outre la mise en
vidence de ladquation de celles-ci selon le type de structures, cette dmarche permettra
plus gnralement de comparer les capacits de techniques rcentes (tomographie X et
mesures de champs thermiques) par rapport aux mthodes plus classiques (mesures de
champs acoustiques et radiographie X).
Ce chapitre sarticule autour de deux parties avec tout dabord lexplication du principe
danalyse des rsultats, puis ltude dtaille de lensemble des rsultats obtenus sur les
cinq structures prsentes prcdemment.
3.2. Principes danalyse 77
3.2 Principes danalyse

Comme cela a t prsent au chapitre prcdent, chacune des mthodes de CND
tudies conduit une visualisation des rsultats par une cartographie de dfauts 2D
ou 3D sous un format image ou assimil, en couleur ou en niveau de gris. De ce fait, le
traitement des rsultats, et notamment leur analyse quantitative, ncessite en gnral un
traitement de ces images.
Le principe danalyse des rsultats se dcompose en deux parties :

lexplication du traitement dimage permettant de dimensionner les dfauts identi-
fis,
le dtail des critres danalyse utiliss sur lensemble des rsultats.
3.2.1 Dimensionnement des dfauts par traitement dimage

Le dimensionnement des dfauts localiss a t ralis par traitement dimage pour
chaque technique de CND :
les rsultats des contrles ultrasonores par des mthodes manuelles de traitement
bases sur les cartographies en paisseur pour les techniques par cho ultrasonore
(figure 2.23-b) et sur les cartographies en amplitude pour la transmission ultrasonore
(figure 2.23-a) : mthodes habituellement utilises dans lindustrie aronautique. On
considrera quil y a prsence dun dfaut quand lamplitude de lcho de fond chute
de 50 % ; on parle galement de mthode -6 dB.
tous les autres rsultats avec une procdure automatique crite en langage Python ;
le choix du langage Python tant li avec lidentification de dfauts par mesures
de champs cinmatiques objet du chapitre suivant et permettant duniformiser le
traitement des donnes.
On se propose dans cette partie de dtailler le fonctionnement de cette dernire procdure
qui a t mise en place afin de sextraire au maximum de la dpendance du rsultat vis-
-vis de loprateur. Il convient de noter que cette dmarche est galement applicable sur
les cartographies ultrasonores. Afin dillustrer les diffrentes tapes de cette procdure,
on se propose justement de lappliquer au C-scan en paisseur du matriau monolithique
avec inserts Tflon obtenu par cho ultrasonore par immersion (figure 2.23-b).
Limage est en premier lieu segmente manuellement de faon isoler chacun des
dfauts identifis (figure 3.1-a). Sur celle-ci, la procdure consiste appliquer les oprations
suivantes :
1. conversion de limage en niveau de gris (figure 3.1-b),
2. recherche du niveau de gris maximum sur labscisse et lordonne de limage (figure

3.1-c),
3. reprage des coordonnes du point identifi ltape prcdente (figure 3.1-d),
4. seuillage de la courbe de niveau de gris reprsentant le dfaut (figure 3.1-e), mi-

hauteur environ, ce traitement est explicit en dtail dans le chapitre suivant,
5. dimensionnement du dfaut correspondant au seuillage dtermin (figure 3.1-f),
6. conversion en mm du diamtre dtermin en pixel par une mise lchelle de limage.
A titre dexemple, le tableau 3.1 rcapitule le type de rsultat obtenu par cette dmarche
pour des techniques de chacune des familles des mthodes de mesures de champs.
(a) dfaut isol (b) Etape 1
(c) Etape 2 (d) Etape 3
(e) Etape 4 (f) Etape 5

Figure 3.1 Procdure de traitement dimage automatique
3.2. Principes danalyse 79
Diamtres
Famille de
Image en Image en rels /
mesure de Seuillage
sortie de test niveau de gris dtec-
champs
ts
Acoustiques
4/4
(cho US par
mm
immersion)
Densimtriques
10/10.84
(radiographie
mm
X)
Thermiques 15/16.53
(TIRM) mm
Table 3.1 Rsultats du traitement dimage par famille de mesures de champs
3.2.2 Critres danalyse

Les huit mthodes de CND ont t testes sur les cinq prouvettes talons dcrites au
chapitre prcdent, ce qui reprsente un traitement de 1160 dfauts. De faon raliser une
comparaison complte et dtaille sous diffrents points de vue, les rsultats exprimen-
taux obtenus ont dabord t analyss par talon puis par type de structure (monolithique
et sandwich). Pour ces tudes, les critres danalyse reposent sur deux aspects importants
du point de vue du CND :
la sensibilit de dtection, soit le nombre de dfauts dtects par la technique ; la
technique juge la plus efficace dans ce domaine est celle qui dtecte le plus grand
nombre de dfauts,
pour les dfauts dtects, la prcision de dtection caractrise par la moyenne des
carts de mesure (en mm) calculs entre la valeur du diamtre du dfaut insr
et celle du dfaut dtecte par le contrle ; la mthode la plus prcise pour ce cri-
tre conduit lcart moyen de mesure le plus petit ; notons que cette prcision sera
apprcie dune part en fonction de la localisation des dfauts en profondeur (la pro-
fondeur des dfauts est dfinie par rapport la source dmission des flux), et dautre
part, en fonction de leur diamtre. Cette dernire analyse apporte notamment des
informations quant la prdisposition des techniques identifier des dfauts consi-
drs comme critiques par le rfrentiel aronautique (de taille suprieure 6 mm
de diamtre (AITM6-4002 [8])) voire galement de plus petite taille.
3.3 Rsultats exprimentaux

Les protocoles exprimentaux concernant les techniques de mesure de champs ther-
miques (TIR, TIRM, TIRP) ont abord les configurations de contrle en face avant et
en face arrire. Toutes les prouvettes fabriques ont t testes selon ces deux modes
mais, afin dallger le propos, seuls les rsultats de la configuration la plus pertinente
pour chaque type de structure sont dtaills dans ce qui suit, prcisment les mesures en
configuration face arrire pour les structures monolithiques et en configuration face avant
pour les structures sandwichs.
Lincapacit de certaines techniques identifier lensemble des dfauts dpend la fois

du matriel de mesure mais galement des chantillons tudis. Cest le couple matriel
utilis/matriau control qui va dfinir la sensibilit et la prcision des rsultats obtenus.
3.3.1 Etalons monolithiques avec inserts Tflon
Cette partie recouvre les rsultats concernant le matriau monolithique avec inserts
Tflon ainsi que les deux zones monolithiques de lprouvette multi-structure.
3.3.1.1 Sensibilit de dtection
Les tableaux 3.2 et 3.3 rcapitulent les analyses de la sensibilit de dtection pour ces
talons monolithiques, respectivement pour lprouvette monolithique et pour les zones
concernes de la structure multi-structure. Trois techniques prsentent une efficacit re-
marquable pour ces deux prouvettes, conduisant la dtection de la totalit des dfauts
(tableau 3.2) :
deux techniques de mesures acoustiques en configuration face avant, savoir les
mthodes dcho ultrasonore par contact et par immersion,
une technique de mesure densimtrique en mode face arrire : la radiographie X.
3.3. Rsultats exprimentaux 81
Table 3.2 Sensibilit de dtection - Matriau monolithique avec inserts Tflon
Table 3.3 Sensibilit de dtection - Zones monolithiques de lprouvette multi-structure

Si lon examine les caractristiques des dfauts non dtects, on remarque que la majo-
rit dentre eux possde un diamtre infrieur au dfaut critique des normes aronautiques
(6 mm de diamtre (AITM6-4002 [8])). De manire gnrale, les techniques en face arrire
sont dautant plus perturbes que les dfauts sont situs en profondeur par rapport la
source dmission. Par exemple pour le matriau monolithique avec inserts Tflon , les
phnomnes de dispersions dondes sont notablement amplifis partir de la profondeur
mdiane (1,4 mm), amenant une rduction de la capacit de dtection des mesures.
Transmission ultrasonore
Bien quefficace pour lprouvette monolithique, cette mthode savre dfaillante pour
les zones monolithiques de lprouvette multi-structure : 52% des dfauts ne sont pas
dtects, dont 84% dinserts de 3 mm de diamtre, et ce quelle que soit la profondeur.
Cette diffrence peut sexpliquer par les tats de surface des structures tudies. En effet,
sur lprouvette n 5 le contrle a t effectu avec le capteur metteur du ct lisse (type
D) et le capteur rcepteur sur la surface rugueuse et brillante (type E). Ainsi, lentre des
ondes ultrasonores est optimise (minimisation des phnomnes de rflexion) mais pas leur
rception. De plus, pour cette technique, la dtection des ondes est ralise uniquement par
le positionnement de la porte A rcuprant lcho de rception. Lors du dimensionnement
-6 dB, la chute damplitude gnre par la prsence de dfauts de 3 mm de diamtre
nest pas suffisante pour permettre lidentification de dfauts.
Tomographie X
En dpit des nombreuses tudes doptimisation dans le traitement des images, lincon-
vnient majeur de la tomographie X reste la prsence ventuelle dartefacts produits lors
de la phase de reconstruction. En effet, ceux-ci soit empchent directement la visibilit
des dfauts, soit perturbent la dtermination prcise de leur contour. Cest ainsi le cas
des 3 dfauts non dtects pour le matriau monolithique (artfact situ 1,4 mm de
profondeur) et de la partie monolithique la plus fine de lprouvette multi-structure. Et
contrairement ce qui est prconis par la norme AITM, la nature des dfauts (Tflon )
nest peut tre pas adquat la tomographie X en rapport la densit et lattnuation
un peu trop proche de la rsine polymre.
TIR
Cette technique dtecte trs difficilement les dfauts de 3 et 4 mm prsents en pro-
fondeur (au-del de 2.1 mm) dans les deux structures tudies. De manire observer le
plus grand nombre de dfauts, le point focal des halognes est positionn au centre de
la zone inspecte. Les dfauts excentrs par rapport ce point sont ds lors difficilement
accessibles. Par ailleurs, la source thermique de TIR met un signal en crneau qui a
tendance se disperser facilement, amenant une sollicitation de faible intensit pour les
dfauts de petite taille localiss en profondeur.
A cette justification, sajoutent dautres perturbations de la propagation du flux ther-
mique associes la composition et au procd de fabrication de lprouvette multi-
structure. La taille du tissage faisant environ 3 mm de large, la vision de larmure est
trs prpondrante sur les images thermiques et du mme ordre de grandeur que la taille
des petits dfauts. De plus, les images du film sont galement trs bruites sur toute
lpaisseur de la structure au niveau des emplacements des grilles destines propager la
rsine lors de linfusion (tat de surface de type E). De ce fait, la dtection des dfauts
localiss en ces endroits est rendue quasiment impossible dans la mesure o le bruit de
mesure est quivalent lcart thermique gnr par chaque dfaut.
TIRM
Les performances de dtection de la TIRM sont bien meilleures que celles de la TIR
grce lemploi du signal priodique basse frquence (sinusodal de frquence 0,05 Hz).
Le signal priodique attnue en effet les phnomnes de dissipation thermique puisque la
pice emmagasine de la chaleur tout au long du contrle et par consquent en transmet
davantage au dtecteur. Dans le mme temps, lutilisation dune onde basse frquence
augmente la capacit de celle-ci se propager dans lpaisseur de la structure inspecte et
diminue le bruit des images des films thermiques. Enfin, les cartographies de phase sont
plus contrastes que les cartographies en amplitude et permettent donc une identification
plus complte des dfauts, notamment des plus profonds.
TIRP
Concernant cette technique, les rsultats obtenus sont quasiment identiques ceux de
la TIR. Comme pour toute technique thermique, le point focal des lampes flash est situ au
centre de la zone contrler de faon raliser un film optimisant la propagation du champ
thermique. Ici, lutilisation dun signal non priodique (impulsion de Dirac) augmente les
phnomnes de dissipation thermique au niveau des extrmits de la structure, rendant
nouveau les petits dfauts en profondeur difficilement identifiables.
3.3.1.2 Analyse de la prcision de dtection
La prcision moyenne de chaque technique est tablie, pour les dfauts dtects, laide
de la moyenne des carts entre la valeur du diamtre rel et celle du diamtre mesur,
quel que soit le diamtre et quelle que que soit sa profondeur. Il est bien entendu que les
techniques capables de dtecter 100% des dfauts sont susceptibles dengendrer des carts
moyens plus important que celles ne dtectant que 40% de dfauts. Lcart moyen calcul
sur un dnominateur commun en terme de dfauts modifie les rsultats (i.e. diminution
des carts moyens pour les techniques dtectant 100% des dfauts) mais le dnominateur
commun est souvent faible, seulement 15% des dfauts environ, et donc non reprsentatif.
Le tableau 3.4 apporte une vue globale de la performance dans ce domaine de chaque
technique. Les mesures de champs densimtriques sont globalement les plus pertinentes
avec un cart moyen infrieur 0.6 mm, suivies de prs par les mthodes par thermogra-
phie : la TIR conduit de meilleurs rsultats pour le matriau monolithique avec inserts
Tflon alors que la radiographie X semble plus adapte lprouvette multi-structure
(carts infrieurs ou gaux 0.3 mm). Dans les deux cas cependant, les techniques les
moins prcises sont les trois mthodes ultrasonores, avec des carts moyens pratiquement
tous suprieurs 1 mm. En particulier, la mthode dcho ultrasonore par immersion
conduit une assez mauvaise estimation sur le cas du matriau monolithique (cart su-
prieur 6 mm), les rsultats ntant pas reprsentatifs et notamment influencs par une
mauvaise compaction de lprouvette, ils ont t supprims des figures suivantes (figures
3.2, 3.3)pour une meilleure lisibilit des autres rsultats.
Table 3.4 Prcision moyenne de dtection - Etalons monolithiques
Afin daller plus loin dans lanalyse, ces rsultats gnraux sont dtaills pour chaque
technique vis--vis des profondeurs des dfauts (figure 3.2) et vis--vis de leur diamtre
(figure 3.3).
(a)
(b)
Figure 3.2 Prcision de dtection en fonction de la profondeur des dfauts - Etalons
monolithiques : (a) matriau monolithique avec inserts Tflon , (b) prouvette multi-
structure
(a)
(b)
Figure 3.3 Prcision de dtection en fonction du diamtre des dfauts - Etalons mono-
lithiques : (a) matriau monolithique avec inserts Tflon , (b) prouvette multi-structure
Mesures de champs acoustiques
Pour chacun des matriaux tests, les carts de dtection de lcho ultrasonore par
contact augmentent sensiblement avec la profondeur des dfauts. Plus le dfaut est proche
de la surface de dtection, plus il y a des risques quil soit confondu avec lcho dentre
ou lcho de fond. Qui plus est, plus le dfaut est profond, plus le chemin parcourir
par londe ultrasonore est important : londe de retour, reue par le traducteur, sera donc
plus attnue et perturbe par les dfauts prsents dans lpaisseur de lprouvette. La
capacit de dtection de la technique est quant elle faiblement impacte par la taille des
dfauts ; notons en particulier quelle se rvle performante pour des dfauts infrieurs
6 mm de diamtre.
Enfin, comme cela a dj t mis en vidence sur la sensibilit de dtection, lanalyse

faite sur la prcision confirme une trs bonne performance de la transmission ultrasonore
par immersion sur le matriau monolithique avec inserts Tflon , dun niveau quasi ana-
logue aux techniques les plus prcises, alors que les rsultats sont plus discutables pour
le cas de lprouvette multi-structure. Cette diffrence de rponse entre les deux pices
sexplique notamment par les tats de surface mis en jeu : le matriau monolithique avec
inserts Tflon possde deux surfaces lisses alors que la plaque multi-structure possde
une surface lisse (o est plac lmetteur) et une surface rugueuse (o se trouve le rcep-
teur). Il semble donc que les rugosits dans cette seconde situation perturbent la sortie
des ondes ultrasonores et engendrent des difficults de dtection. Du point de vue des
diamtres dtects, on soulignera que la prcision de cette technique tend samliorer
globalement avec laugmentation des taille de dfauts.
Mesures de champs densimtriques
Les constats concernant la radiographie X et la tomographie X sont comparables.

Elles conduisent des rsultats prcis, cohrents pour les deux types dprouvettes et
pratiquement constants quelles que soient la profondeur des dfauts et leur taille.
Mesures de champs thermiques
Dans la gamme daccessibilit de la TIR (profondeurs de dfauts infrieures 2.1 mm),

les rsultats obtenus conduisent pour les deux structures une prcision quivalente aux
mthodes bases sur les mesures de champs densimtriques, et ce de manire constante
quelle que soit la profondeur. Pour le matriau monolithique, on note galement une
amlioration avec laugmentation du diamtre des dfauts, cette tendance restant limite
pour lprouvette multi-structure du fait des bruits de mesure sur les films thermiques
voqus lors de lanalyse de la sensibilit de dtection.
Du point de vue de la prcision, la performance de la mthode TIRM est un peu

moins bonne que celle de la TIR. Elle dcrot globalement avec la profondeur des dfauts
ainsi quavec la taille de ceux-ci. Ainsi, si la priodicit du signal permet damliorer la
sensibilit de dtection, elle apporte en revanche des difficults supplmentaires sur le
dimensionnement des dfauts puisque le signal dtect est en perptuel mouvement tout
au long du contrle. De ce fait, limage optimale pour effectuer le dimensionnement est
parfois difficile obtenir et conditionne la performance des rsultats.
Le traitement des images dvelopp par THERMOCONCEPT facilite enfin norm-

ment lanalyse des rsultats de la TIRP puisque le film global est transform en une unique
image reprsentative des proprits thermiques et conduit des prcisions de mesure com-
parables la radiographie X. En terme de tendance, elle suit une volution en fonction
des profondeurs et diamtres pratiquement quivalente la TIRM.
3.3.2 Etalon monolithique avec trous fond plat
Abordons prsent selon la mme mthodologie danalyse le cas du matriau mo-

nolithique avec trous fond plat. Pour cette plaque, la face avant o sont places les
sources dmission correspond la surface plane et la face arrire comporte les trous. La
profondeur des dfauts est indique par rapport la face avant.
Daprs le tableau 3.5, trois mthodes sur huit dtectent la totalit des dfauts :
lcho ultrasonore par immersion (configuration face avant),
les deux mthodes de champs densimtriques, radiographie X et tomographie X
(configurations face arrire).
On peut associer ces mthodes les plus efficaces lcho ultrasonore par contact (confi-
guration face avant) dont le seul dfaut non dtect est le plus petit dfaut (2 mm de
diamtre) le plus proche de la surface de dtection (1 mm de profondeur) pour lequel
lcho de dfaut est difficilement dissociable de lcho dmission.
Table 3.5 Sensibilit de dtection - Matriau monolithique avec trous fond plat
La majorit des dfauts non dtects sont localiss aux profondeurs les plus impor-
tantes (5.2 et 7.48 mm). Examinons dans le dtail le cas des techniques qui prsentent
des limitations pour ce cas de figure.
Transmission ultrasonore par immersion

En transmission ultrasonore, le faisceau dondes doit traverser lpaisseur complte de
la pice inspecte et savre donc considrablement attenu lorsquil arrive au dtecteur
dans le cas de pices paisses. En prsence de dfauts, le signal reu est encore plus faible.
De manire gnrale, la dtection de dfauts de petite ou moyenne taille (infrieure 6
mm) est rendue difficile avec cette technique du fait de la faible attnuation du signal
gnre par ces derniers. La taille critique est donc ici de 3 mm (aucun dfaut de diamtre
infrieur ou gal 3 mm nest dtect) et ce phnomne saccrot dautant plus que les
dfauts sont profonds (les dfauts de 6 mm ne sont plus accessibles au-del de 5.2 mm, le
seul dfaut visible 7.48 mm de profondeur fait 10 mm de diamtre).
TIR
A nouveau, la diffusion du flux thermique des halognes et donc la profondeur din-
vestigation de cette technique est limite dans cette pice paisse. Aucun dfaut nest
ainsi identifi 7.48 mm de profondeur alors quils le sont tous 0.52 mm de profon-
deur. De manire gnrale, la profondeur limite de dtection dun dfaut critique pour
laronautique (6 mm de diamtre (AITM6-4002 [8])) vaut ici 2.8 mm de profondeur. En
comparaison, celle dun dfaut de 10 mm de diamtre est de 5.2 mm de profondeur.
TIRM
Les rsultats de la TIRM sont peu prs quivalents ceux de la TIR. Nanmoins,
on remarque une amlioration notable concernant lidentification des dfauts localiss aux
faibles profondeurs (entre 0.52 et 2.8 mm) pour lesquelles le signal priodique semble plus
adapt la dtection.
TIRP
Bien moins performante que les autres mesures par thermographie, la TIRP ne dtecte
sur cette prouvette que les dfauts proches de la surface (0.52 mm de profondeur). La
puissance du flash thermique ne semble donc pas tre suffisante vis--vis de lpaisseur
inspecte.
Les rsultats concernant la prcision moyenne sur les dfauts dtects sont prsents
la figure 3.6. Lorsquils sont visibles par les diffrentes techniques, les trous fond plat
engendrent davantage dhtrognits dans les champs que les inserts, et du coup sont
facilement dtects. Ainsi, sur les huit mthodes employes,
cinq dentre-elles (cho ultrasonore par contact, cho ultrasonore par immersion,
radiographie X, tomographie X et TIR) prsentent un cart de dtection global
infrieur ou gal 0.5 mm ; en particulier, la tomographie X est la plus efficace ;
les trois autres (transmission ultrasonore par immersion, TIRM, TIRP) ne dpassent
pas un cart de 2 mm.
Table 3.6 Prcision moyenne de dtection - Matriau monolithique avec trous fond
plat
Comme prcdemment, les figures 3.4 et 3.5 illustrent lanalyse des prcisions par
profondeur et par diamtre de dfauts.
Figure 3.4 Prcision de dtection en fonction de la profondeur des dfauts - Matriau

monolithique avec trous fond plat
Figure 3.5 Prcision de dtection en fonction du diamtre de dfauts - Matriau mo-

nolithique avec trous fond plat
Lcho ultrasonore par contact conduit, comme dans le cas des structures monoli-
thiques, de trs bons rsultats et fait preuve dune capacit didentification relativement
constante quelle que soit la profondeur des dfauts et leur taille.
En ce qui concerne lcho ultrasonore par immersion, la courbe de la prcision en

fonction de la profondeur indique de moins bons rsultats dans les zones proches de la
surface de mesure et celles proches de la surface de fond, que lon attribue la mauvaise
identification de lcho dentre ou de lcho de fond par rapport lcho de dfaut.
Paralllement, ses performances augmentent logiquement avec le diamtre des dfauts,
avec une rponse trs satisfaisante pour les dfauts de 6 et 10 mm.
Bien que la prcision samliore avec laugmentation des tailles de dfauts, la trans-
mission ultrasonore est le contrle qui gnre les carts les plus importants, et ce sur la
quasi totalit des profondeurs des dfauts. La particularit de la transmission ultrasonore
par immersion rside dans le placement dune seule porte sur lcho de rception (porte
A). Londe ultrasonore tant fortement attenue par lpaisseur de la pice, lamplitude
de cet cho de rception est faible. De ce fait, il faut ajouter du gain pour amplifier le
signal, et par voie de consquence le bruit de mesure. Cette perturbation, associe la
faible attnuation provoque par les dfauts sur une partie aussi paisse expliquent ces
imprcisions des mesures obtenues.
Comme pour les inserts, les rsultats de radiographie X et tomographie X sont trs
bons avec une bonne stabilit de la prcsion quelle que soit la profondeur et la taille des
trous.
Le nombre de dfauts dtects avec la TIR diminue avec laugmentation de la profon-

deur de ces derniers, les rsultats sont ds lors lgrement fausss et ne permettent pas
lobservation dune vraie tendance (la prcision didentification 5.2 mm ne fait inter-
venir quun dfaut de 10 mm de diamtre). Nanmoins, pour la gamme des profondeurs
comprises entre 0.52 et 2.8 mm, et des tailles de dfauts entre 3 et 10 mm de diamtre,
la dtection est bonne et stable avec des carts infrieurs 1 mm.
Les rsultats de la TIRM sont assez htrognes, tant pour lanalyse des profondeurs
que des tailles. Dans la mesure o un seul dfaut est identifi 1 et 5.2 mm de profondeur et
quaucun dfaut nest dtect 7.48 mm de profondeur, les profondeurs reprsentatives
des performances de cette mthode sont 0.52 et 2.8 mm. Plutt satisfaisante sur des
dfauts proches de la surface dmission (0.5 mm dcart 0.52 mm de profondeur), cette
technique est trs sensible laugmentation de la profondeur (plus de 2 mm dcart 2.8
mm de profondeur). Cet aspect se confirme au niveau des diamtres des dfauts avec le
dfaut de 2 mm parfaitement mesur se situant proche de la surface dmission. Pour les
dfauts de 3 10 mm, la tendance globale est lamlioration avec laugmentation des
tailles, linstar de la TIR. En revanche, si le signal priodique amliore la sensibilit
de dtection par rapport cette dernire, la prcision de la TIRM est globalement plus
faible.
Lanalyse de la TIRP reste limite car seuls sont reprs les dfauts prsents 0.52
mm de profondeur. Qui plus est, le rsultat pour chaque diamtre ne correspond qu un
seul dfaut. La prcision dans ce cas est moyenne (cart lgrement suprieur 1.5 mm),
induite par des phnomnes dovalisation sur les images dautant plus importants que les
dfauts sont grands.
3.3.3 Etalons sandwich avec inserts Tflon

Les rsultats des tests sur les prouvettes sandwichs ainsi que la partie sandwich de la
plaque multi-structure vont tre dtaills maintenant.
Seuls les rsultats des mthodes de mesures ultrasonores et thermiques mens en confi-
guration face avant sont prsents ici. La transmission ultrasonore nest par ailleurs pas
intgre dans lanalyse qui suit car elle na pas permis de dtecter le moindre dfaut. En
effet, lme en mousse tant un isolant acoustique, les ondes ultrasonores mises narrivent
pas traverser les matriaux tudis et tre dtectes par le traducteur rcepteur.
Dautre part, afin de faciliter la comprhension des phnomnes expliqus par la suite,
on rappelle brivement certains aspects essentiels concernant ces matriaux. Chacun des
talons est constitu de fines peaux : lpaisseur des peaux de la premire prouvette est
de 0.66 mm alors que celles des deux autres prouvettes est de 1.08 mm. Par ailleurs, les
tats de surface de chacune des structures conditionne la qualit de la mesure :
le sandwich avec peaux de 0.66 mm possde un tat de surface lisse et brillant du
ct moule (type A) et un tat de surface rugueux, bossel et mat du ct contre-
moule (type B) ; les mesures pour cette prouvette ont t ralises du ct bossel ;
le sandwich avec peaux de 1.08 mm prsente un tat de surface lisse et brillant (type
A) du ct moule et un tat de surfacerugueux (type C) du ct contre-moule ; la
face contrle correspond dans ce cas au ct lisse ;
lprouvette multi-structure a quant elle un tat de surface lisse et mat du ct
moule (type D) et un tat de surface brillant et bossel du ct contre-moule (type
E) ; pour cette dernire, les deux cts ont t tests et lon prsente dans la suite la
moyenne des rsultats obtenus qui savre donc reprsentative dun tat de surface
de nature intermdiaire entre les deux types de cts.
Daprs les tableaux 3.7, 3.8 et 3.9, les mthodes de mesures de champs densimtriques
sont les plus efficaces pour le contrle des structures sandwichs. La radiographie X d-
tecte notamment la totalit des dfauts insrs dans les trois prouvettes tudies. La
tomographie X identifie quant elle lensemble des dfauts prsents sur les prouvettes
sandwichs mais sa mesure est lgrement perturbe par des artefacts pour lprouvette
multi-structure. En ce qui concerne la qualit des mesures ultrasonores et thermiques, on
note des disparits suivant les prouvettes contrles.
Table 3.7 Sensibilit de dtection - Matriau sandwich avec peaux de 0.66 mm
Table 3.8 Sensibilit de dtection - Matriau sandwich avec peaux de 1.08 mm

Table 3.9 Sensibilit de dtection - Zones sandwichs de lprouvette multi-structure
Comme prcdemment, examinons dans le dtail les sensibilits de dtection de cha-

cune des techniques.
Echo ultrasonore par contact

Trois impratifs majeurs doivent tre respects pour la mise en oeuvre du contrle
en cho ultrasonore par contact sur des matriaux sandwichs avec peau de fine paisseur
(infrieure ou proche de 1 mm) :
pouvoir diffrencier lcho dentre dans la pice, lcho de fond et lcho de dfaut,
assurer la visibilit de ces chos dans le A-scan en ayant des chos de fond et de
dfaut damplitude suffisante,
maintenir (manuellement) une liaison acoustique constante et homogne.
Les mauvais rsultats obtenus pour les trois prouvettes, notamment pour la premire,
sont associs des difficults dans ces trois domaines. En effet, le sandwich avec peaux de
0.66 mm possde dune part la peau la plus fine et, dautre part, sa surface de dtection
est rugueuse et bossele. Une paisseur constante dagent couplant ne peut alors tre
assure ; le faisceau dondes tend tre dvi et les portes de dtection (portes I, B
et A) sont difficiles positionner. Les performances sont meilleures sur les deux autres
prouvettes pour lesquelles la technique a t teste sur des cts plus lisses. On notera
en revanche lincapacit de cette technique accder aux dfauts localiss 1.08 mm de
profondeur dans la structure multi-structure. Cette difficult sexplique par la position du
dfaut entre la peau et la mousse.
Echo ultrasonore par immersion

Lcho ultrasonore par immersion dtecte la totalit des dfauts sur la seconde prou-
vette qui a t teste du ct lisse, renforant ainsi la qualit de limmersion par rapport
au contact. En revanche, aucun dfaut nest dtect sur le premier talon pour les mmes
raisons que celles voques prcdemment (peau de fine paisseur, tat de surface bossel).
En ce qui concerne lprouvette multi-structure, son tat de surface moyen est moins
lisse que la seconde prouvette, les seuls dfauts non dtects mesurent 3 mm et sont
localiss prs des transitions vers les parties monolithiques o les effets de bord tendent
diminuer fortement le pouvoir de dtection de la technique.
TIR
Les mesures thermiques en configuration face avant sont de manire gnrale sensibles
deux phnomnes :
la saturation thermique trs rapide des structures suite lmission du flux, qui
rduit le temps de visibilit des dfauts,
limpact thermique du point focal de lhalogne qui minimise la visibilit des dfauts
cet endroit.
Plus la surface est brillante, plus ces aspects sont accentus et affectent la capacit de
dtection de ces mesures. Les rsultats obtenus avec la TIR confirment ce propos : cette
technique dtecte la totalit des dfauts dans la premire structure (tat de surface mat)
et moins de 50% pour les autres talons (tats de surface brillant et intermdiaire). Dans
le premier cas, il faut noter que le doublement des inserts accrot sans doute galement les
carts thermiques et amliorent leur visibilit. Les deux dernires structures, fabriques
partir du mme tissu carbone, sont galement perturbes par le bruit gnr par son ar-
mure vis--vis des dfauts de 3 et 6 mm. Enfin, la diffrence entre les rsultats de celles-ci
(50% pour lprouvette sandwich paisse, 35% pour lprouvette multi-structure daspect
moyen pourtant plus mat) sexplique en grande partie par leffet isolant induit par les nom-
breuses tapes de prformage ncessaires la fabrication de lprouvette multi-structure
(modification physico-chimique) ainsi que par les effets de bords dus la gomtrie de
lprouvette.
TIRM
Les commentaires prcdants pour la TIR peuvent tre tendus la TIRM. Nanmoins,
lutilisation dun signal priodique et le traitement des images sur la cartographie de phase
permettent de diminuer le bruit de mesure li larmure du tissu et donc damliorer la
sensibilit de dtection pour lprouvette multi-structure.
TIRP
A nouveau, les conclusions sur la TIR sappliquent cette dernire mthode par ther-
mographie. Le traitement dimages spcifique ne permet cependant pas de sextraire to-
talement de la perturbation lie larmure du tissu de lprouvette multi-structure.
Le tableau 3.10 rcapitule, pour les dfauts dtects, les carts moyens obtenus pour
chacune des techniques. A nouveau, les mesures de champs densimtriques sont les plus
performantes (cart moyen infrieur 0.5 mm), les mesures thermiques donnent de trs
bons rsultats (avec un cart moyen infrieur 1 mm hormis pour la TIR sur le premier
matriau), et les mesures acoustiques restent un peu en retrait (cart moyen autour de 2
mm).
Table 3.10 Prcision moyenne de dtection - Etalons sandwichs
Lvolution de ces carts de mesure vis--vis des profondeurs (figures 3.6, 3.7 et 3.8)
et des diamtres de dfauts (figures 3.9 , 3.10 et 3.11) compltent lanalyse. Il convient
de noter pour la premire prouvette que les dfauts situs 0.22 mm de profondeur font
tous 6 mm de diamtre, que ceux 0.44 mm de profondeur font tous 3 mm de diamtre et
que les deux diamtres de dfauts sont prsents 0.66 mm de profondeur. Ceci explique
notamment le pic dcart observ 0.44 mm de profondeur sur la figure 3.6 pour de
nombreuses techniques. Pour les deux autres prouvettes, les dfauts de 3 mm et 6 mm
de diamtre sont rpartis quitablement selon la profondeur.

sandwich avec peaux de 0.66 mm

Figure 3.8 Prcision de dtection en fonction de la profondeur des dfauts - Zones

sandwichs de lprouvette multi-structure
Figure 3.9 Prcision de dtection en fonction du diamtre de dfauts - Matriau sand-

wich avec peaux de 0.66 mm
Figure 3.10 Prcision de dtection en fonction du diamtre de dfauts - Matriau

Figure 3.11 Prcision de dtection en fonction du diamtre de dfauts - Zones sandwichs

de lprouvette multi-structure

Compare aux techniques tudies, la prcision de lecho ultrasonore par contact sur
ces structures sandwichs est assez faible. Affecte par un tat de surface bossel, elle
semble galement sensible au procd de fabrication. Le premier talon est fabriqu
partir de tissus pr-imprgns alors que les deux autres sont raliss par infusion de rsine
liquide. Dans ces derniers cas, la prsence de rsine en plus forte quantit entre le dernier
pli et la mousse, et dans une moindre mesure, au sein de la mousse elle-mme doivent
diminuer lattnuation des ondes ultrasonores dans lme. Lidentification de lcho de
fond est beaucoup moins dlicate et donc la dtection des dfauts est plus prcise en
profondeur quen surface pour ces deux talons, contrairement au premier matriau qui
subit pleinement leffet de la mousse. En terme de taille, comme il a pu tre constat
prcdemment, cette mthode est plus prcise sur les dfauts de 6 mm de diamtre que
sur ceux de 3 mm.
Lcho ultrasonore par immersion apporte une rponse en terme de prcision globale-
ment meilleure que la prcdente mthode ultrasonore, grce notamment aux amliora-
tions associes son protocole exprimental. Si lvolution avec la profondeur et la taille
suit la mme tendance que lcho par contact sur le matriau sandwich avec peaux de 1.08
mm, les rsultats obtenus sur ces deux critres pour lprouvette multi-structure sont plus
instables. Une des explications est lutilisation dun procd de fabrication lgrement
modifi entre les deux talons (notamment vis--vis des tapes de prformage).
Les mauvais rsultats obtenus par transmission ultrasonore par immersion sexpliquent
par le manque dadquation entre le matriel utilis (frquence du capteur) et le type de
structure contrle (paisseur de la pice). En effet, actuellement de nombreuses pices
sandwichs aronautiques sont inspectes par transmission ultrasonore par immersion avec
des frquences beaucoup plus faibles (< 1MHz).
Sur ces talons sandwichs, les mesures de radiographie X et tomographie X sont

nouveau trs proches, dune grande prcision et stables par rapport la profondeur din-
vestigation et la taille des dfauts.
Si lon extrait le cas de la premire prouvette o la TIR prsente les mmes spcificits
de prcision que lcho par contact, cette mthode thermique conduit pour les deux autres
talons de bons rsultats. De manire gnrale, la prcision diminue sensiblement avec
la profondeur, de faon plus marque pour la premire prouvette du fait de la mauvaise
dtection des petits dfauts situs en profondeur. L encore, les prcisions les plus grandes
sont obtenues sur des dfauts de plus grand diamtre.
Le signal priodique de la TIRM et le signal impulsionel associ au traitement dimages

spcifique de la TIRP amliorent les rsultats de ces mthodes par rapport la TIR
pour le premier talon. En revanche, les rsultats pour les deux autres prouvettes sont
relativement analogues pour ces trois mthodes. Notons cependant que la TIRP se rvle
la plus stable quelle que soit la profondeur grce une meilleure performance sur la
dtection des dfauts de 3 mm (voir par exemple sur le premier talon, figures 3.6 et 3.9).
3.3.4 Synthse
On se propose pour conclure cette partie dtablir un bilan gnral des performances
des diffrentes techniques, regroupes selon les trois familles de mesures de champs.
Compte tenu des problmatiques industrielles, le point de vue retenu pour ce bilan est
celui dune analyse par types de structures, monolithiques et sandwichs, et vis--vis de la
taille des dfauts. Les moyennes ont t calcules partir des 1160 rsultats individuels et
non partir des tableaux prcdents lis la prcision moyenne de dtection qui portent
sur des chantillons diffrents.
3.3.4.1 Structures monolithiques
Le tableau 3.11 et le graphique 3.12 prsentent un bilan des moyennes obtenues pour
toutes les structures monolithiques tudies (prouvettes monolithiques et prouvette
multi-structure, avec inserts et trous fond plat), respectivement vis--vis de la sensi-
bilit de dtection et de la prcision moyenne. Pour ces matriaux, on rappelle que les
rsultats des mthodes thermiques concernent la configuration face arrire.
En terme de sensibilit de dtection (tableau 3.11), les mthodes acoustiques et den-

simtriques sont globalement les plus performantes sur les structures monolithiques dans
Table 3.11 Sensibilit de dtection - Structures monolithiques
la gamme dpaisseurs considres (0.72 8 mm). En particulier, les techniques dcho

ultrasonore (contact et immersion) et la radiographie X prsentent une efficacit remar-
quable. Les techniques bases sur les mesures de champs thermiques, ralises en face
arrire, identifient quant elles environ la moiti des dfauts en prsence, du fait dune
limitation de dtection en profondeur.
Figure 3.12 Prcision moyenne de dtection - Structures monolithiques
Lanalyse des prcisions moyennes montre en revanche une trs bonne capacit des
mthodes densimtriques et thermiques avec un cart moyen de dtection infrieur 1
mm, voire infrieur 0.5 mm pour la radiographie X et la tomographie X qui sont moins
affectes par la profondeur des dfauts (figure 3.12). Il faut noter en particulier les trs bons
rsultats de la TIR pour ces structures monolithiques, dune prcision globale analogue
la radiographie X. Avec le matriel notre disposition, les mesures de champs acoustiques

savrent gnralement en retrait sur ce critre, mme si lcho ultrasonore par contact
approche la prcision de mesure de la TIRM. Le mauvais rsultat de lcho ultrasonore par
immersion est quant lui fortement accentu par le problme de porosit sur lprouvette
monolithique avec inserts, le niveau de prcision moyenne de cette technique sans inclure
cet talon valant environ 1.68 mm.
3.3.4.2 Structures sandwichs
Une analyse identique est propose pour les structures sandwichs (prouvettes mo-
nolithiques et prouvette multi-structure). Dans ces cas, les techniques thermiques ont
t mises en oeuvre en configuration face avant et la transmission ultrasonore par im-
mersion a chou du fait dune trop forte attnuation du signal par la mousse dans cette
configuration (matriel utilis/matriau contrl).
Table 3.12 Sensibilit de dtection - Structures sandwichs
Y compris sur les structures sandwichs, la sensibilit de dtection la plus performante

est celle des mthodes de mesures de champs densimtriques avec une forte proportion de
dfauts dtects, la radiographie X dtectant en particulier lintgralit des dfauts sur
toutes les prouvettes (tableau 3.12). Les mthodes ultrasonores et thermiques sont quant
elles globalement quivalentes avec en moyenne 60% des dfauts dtects du fait des
perturbations induites par les tats de surface dans cette configuration face avant :
par le relief en mesure ultrasonore, notamment lcho par contact,
par la brillance en mesure thermique.
La trs bonne performance de la radiographie X et de la tomographie X se confirme

au niveau des prcisions moyennes, avec un cart infrieur 0.5 mm (figure 3.13). Outre
la configuration de test en face avant, les contraintes lies ltude des sandwichs (peaux
de fine paisseur, me absorbante) engendrent, en moyenne, pour les autres mthodes
une perte de qualit notable par rapport aux structures monolithiques 1 . Nanmoins,
on observe que les mthodes par thermographie sont un peu moins impactes par ces
phnomnes et conduisent des prcisions qui restent intressantes, de lordre du mm
pour la TIRM et la TIRP.
Figure 3.13 Prcision moyenne de dtection - Structures sandwichs
3.3.4.3 Taille des dfauts
En terme de taille de dfaut, il faut souligner la capacit des mesures de champs den-
simtriques mesurer avec une grande prcision les dfauts critiques pour laronautique
(diamtre 6 mm AITM6-4002 [8] (tableau 3.13 et figure 3.14)) et mme ceux de taille
infrieure (cart infrieur 0.5 mm). Les rsultats des mthodes par thermographie sont
1. Pour lcho ultrasonore par immersion, ce constat est fait par rapport au niveau de prcision a
priori le plus reprsentatif, soit celui obtenu sans inclure ltalon monolithique avec inserts.
galement satisfaisants pour la gamme 6 mm avec un cart gnralement infrieur 1 mm

(tableau 3.13 et figure 3.14). En revanche, il est difficile de conclure de manire gnrale
sur les techniques ultrasonores dans la mesure o les rsultats sont assez disparates selon
les matriaux, le type de dfauts et la technique de CND employe (tableau 3.13 et figure
3.14).
Table 3.13 Sensibilit de dtection - Dfauts de 6 mm de diamtre
Figure 3.14 Prcision moyenne de dtection - Dfauts de 6 mm de diamtre

3.4 Conclusion
Cette tude comparative a permis de comparer des techniques de CND bases sur des
mesures de champs acoustiques, densimtriques et thermiques :
les mesures ultrasonores sont les plus utilises dans le milieu aronautique,
la radiographie X est peu employe dans les structures composites,
la tomographie X et les mesures de champs thermiques connaissent actuellement un
important essor.
Lensemble des dfauts artificiels considrs dans cette tude (inserts Tflon et trous
fond plat) sont en gnral pertinents pour les contrles ultrasonores mais beaucoup
moins reprsentatifs de dfauts de cohsion (fissurations, cavits, dlaminages) du point
de vue de leur densit ou des proprits thermiques. En dpit dun contraste qui nest
donc pas optimal, les mthodes de mesures de champs densimtriques et thermiques ont
rvl sur les matriaux tudis des capacits de dtection au moins quivalentes et sou-
vent meilleures que les mthodes de CND courantes de laronautique. Lanalyse des
sensibilits de dtection et des prcisions dmontre en effet dexcellentes performances
de la radiographie X et de la tomographie X, quels que soient le type de structure (mo-
nolithique ou sandwich), la profondeur ou la taille des dfauts. Les mesures de champs
acoustiques et thermiques montrent enfin une relle complmentarit au regard de len-
semble des rsultats sur les deux types de structures composites. Dautres paramtres
pourraient galement tre intgrs cette analyse, savoir : le temps de mise en oeuvre
et de contrle, non favorable aux techniques acoustiques ; le cot du matriel et des in-
frastructures ncessaires, non favorable aux techniques densimtriques. Enfin, la notion
de contrle face avant ou face arrire est aussi importante prendre en compte, car la
face arrire nest pas toujours accessible ce qui peut tre non favorable aux techniques
densimtriques.
Lintrt de la dtection de dfauts en production comme en maintenance rside dans

la possibilit dune rparation. Se posent alors les questions de la localisation du dfaut
et de la dtermination de sa taille. La localisation est tridimensionnelle, cest--dire quil
faut connatre la profondeur. Pour identifier la profondeur du dfaut, seules les tech-
niques dcho-ultrasonore et de tomographie X sont aujourdhui capables de fournir cette
information (les dveloppements en thermographie IR sont en cours). Dans ce travail,
cette information de profondeur na pu tre identifie que par la tomographie X, avec des
prcisions identiques que celles dans le plan.
On peut enfin conclure sur des choix ou positions qui rsultent de cette tude :
3.4. Conclusion 109
Procds de fabrication : les techniques acoustiques (ultrasons) ne sont pas adaptes

pour des tats de surface rugueux induits par la prsence de tissus darrachage, ou
grilles de drainage de rsine ; les techniques thermiques ne sont pas adaptes pour des
tats de surface brillants induits par la prsence de certains consommables (micro-
perfor, dmoulant) ou dtats de surfaces doutillages trs faible rugosit ; les
techniques densimtriques sont indpendantes de ces tats de surface,
Structures sandwichs : les techniques acoustiques et thermiques ne sont pas adaptes
pour les structures sandwichs, il est souvent ncessaire davoir accs aux deux peaux
et si les peaux sont peu paisses (infrieur au mm) les techniques thermiques sont
les seules adquates ; il faut donc prfrer les techniques densimtriques qui sont trs
bien adaptes ce type de structures,
Position et taille des dfauts dans les structures monolithiques : les techniques acous-
tiques permettent de mesurer des dfauts de diamtre 3 mm ou plus dans un inter-
valle de profondeur compris entre 1 et 5 mm, les techniques thermiques permettent
de mesurer des dfauts de diamtre 3 mm ou plus dans un intervalle de profon-
deur compris entre 0 et 1.5 mm, les techniques densimtriques nont pas de limite
en taille et en profondeur. En gnral les dfauts proches de la surface (quelques
diximes de profondeur) sont les plus difficiles identifier. Il faut enfin noter que la
prcision ntant pas la mme suivant les techniques, un dfaut de diamtre 6 mm
sera identifi comme un dfaut de 7.5 mm par les techniques acoustiques, de 7 mm
par les techniques thermiques et de 6.5 mm par les techniques densimtriques, le
critre dacceptabilit du dfaut sera interprt diffremment.
Chapitre 4
Identification des dfauts par mesures

de champs cinmatiques
Sommaire
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
4.2 Mthodes de mesures de champs . . . . . . . . . . . . . . . . . 113
4.2.1 Mesures surfaciques (2D) . . . . . . . . . . . . . . . . . . . . . 113
4.2.2 Mesures surfaciques et hors-plan (2D 1/2) . . . . . . . . . . . . 115
4.2.3 Mesures volumiques (3D) . . . . . . . . . . . . . . . . . . . . . 118
4.3 Mthodes didentification du comportement . . . . . . . . . . 119
4.3.1 Mthodes itratives . . . . . . . . . . . . . . . . . . . . . . . . . 120
4.3.2 Mthodes non itratives . . . . . . . . . . . . . . . . . . . . . . 123
4.4 Mthodologie didentification multi-zones . . . . . . . . . . . . 125
4.5 Identification des proprits dun matriau orthotrope plan . 128
4.5.1 Modles de rfrence et de recalage . . . . . . . . . . . . . . . . 129
4.5.2 Description du fonctionnement de lalgorithme didentification . 131
4.5.3 Rsultats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
4.6 Identification des proprits dune zone dgrade . . . . . . . 135
4.6.1 Cration de la cartographie derreurs . . . . . . . . . . . . . . . 136
4.6.2 Estimation de la taille de la zone dgrade . . . . . . . . . . . . 139
4.6.3 Identification des proprits de la zone dgrade . . . . . . . . . 141
4.7 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144
112 4. Identification des dfauts par mesures de champs cinmatiques
4.1 Introduction
La dmarche du contrle non destructif consiste identifier des zones de faiblesse au

sein des structures et en mesurer leur tendue. Cette approche, aborde prcdemment
sous langle de diverses techniques exprimentales, permet daccder une localisation
spatiale des zones dgrades et dintroduire des alertes dans les procds de fabrication
et/ou maintenance si la taille des dfauts excde une taille critique dfinie selon les appli-
cations.
On se propose dans cette dernire partie du mmoire de complter ce diagnostic de

la sant matire au travers dune valuation quantitative des proprits de ces zones
dgrades. En plus de laspect dimensionnel des dfauts capt par les mthodes de CND,
il semble en effet intressant et complmentaire de pouvoir dfinir leur taux de dgradation
par rapport aux proprits du matriau sain. Ces deux informations conduiraient ainsi
une caractrisation plus complte, dun intrt crucial pour la prise de dcision par
rapport la pice contrle et de son ventuelle rparation.
A cette fin, notre objectif est de mettre en place une procdure didentification des
proprits lastiques de matriaux macroscopiquement htrognes partir de mesures
de champs cinmatiques. En lien avec les dveloppements de mthodes de mesures telles
que la stro-corrlation dimages numriques, lidentification paramtrique a fait lobjet
dtudes rcentes au sein du laboratoire vis--vis du comportement lasto-plastique de
matriaux mtalliques isotropes (Fazzini [76], Fazzini et al. [77] ). Dans le cadre de ce
travail, la premire difficult est donc dtendre ces travaux au cas des composites propres
au projet FUSCOMP, en loccurrence des matriaux orthotropes plans. A partir de ce
rsultat, une des contributions majeures de ce travail consiste tablir une mthodo-
logie danalyse multi-zones, qui vise dterminer la variation spatiale des proprits
macroscopiques du matriau (zones saines et zones dgrades).
Ce chapitre comprend en premier lieu une prsentation bibliographique des techniques

de mesures de champs cinmatiques ainsi que des principales mthodes didentification.
Aprs une explication des grandes lignes de la dmarche employe dans cette tude, les
procdures spcifiques dans les cas macroscopiquement homognes et macroscopiquement
htrognes seront dtailles. Les performances des algorithmes dvelopps en Python
seront finalement illustres au travers de quelques rsultats.
4.2. Mthodes de mesures de champs 113
4.2 Mthodes de mesures de champs

De nombreuses mthodes optiques dimensionnelles de mesure de champs peuvent tre
rencontres dans la littrature. Avec une sensibilit mcanicienne, nous avons choisi de les
classer en trois familles selon la dimension de mesure ralise. Elles peuvent en effet tre
utilises pour des mesures surfaciques dites 2D , des mesures surfaciques et hors-plan
dites 2D 1/2 ou des mesures volumiques dites 3D.
Au cours de cette revue, une attention particulire sera consacre la stro-corrlation
dimages numriques. Cette technique est en effet employe depuis dix ans environ au sein
du Laboratoire Gnie de Production et les dmarches didentification dj mises au point
par lquipe Mcanique des Matriaux, des Structures et des Procds sappuient sur les
mesures issues de celle-ci.
4.2.1 Mesures surfaciques (2D)

Les mesures surfaciques permettent dobtenir un champ de dplacement 2D sur une
surface plane. Dans un repre (O, x, y, z), les composantes planes du dplacement u et v
sont mesures et elles sont fonction de x, y variables et z0 constante (quation 4.1) :
(
u(x, y, z0 )
(4.1)
v(x, y, z0 )
Le gradient de dplacement F dcrit par lquation 4.2 ne permet pas de calculer le tenseur
de Green-Lagrange. Seules les composantes planes du tenseur des petites dformations
=
sont accessibles (quation 4.3).

Fxx Fxy
F = Fyx Fyy (4.2)

xx xy
1
= FT + F I = xy yy (4.3)

= 2 = = =

Les techniques de mesure surfaciques les plus couramment rencontres sont la photo-
lasticit, linterfromtrie, le moir et les mthodes utilisant lanalyse dimages.
La photolasticimtrie 2D Cest une technique exprimentale danalyse des dfor-

mations et des contraintes. La mthode est base sur la proprit de double rfraction ou
birfringence de certains matriaux transparents comme les rsines poxy ou le polym-

thylmthacrylate (PMMA). Des franges isochromes et isoclines reprsentant la rpartition
de la dformation (diffrence de dformations principales et directions principales) sont
visibles avec laide dun polariscope (Chen [54], Jessop and Harris [114], Morimoto et al.
[150], Patterson and Wang [160]).
Le moir Le moir consiste faire apparatre des franges dinterfrence entre une grille
de rfrence (rseau rgulier de lignes) et une seconde grille, identique la grille de rf-
rence, qui est dpose sur la surface de la pice tudier et qui se dforme avec celle-ci. Les
franges sont alors analyses par le biais dun interfromtre moir qui permet de mesurer
des dplacements et des dformations la surface de la pice (Post et al. [167]).
Les techniques de corrlation - corrlation dimages et mthode des grilles Ces

techniques permettent la dtermination des champs de dplacement et de dformation
la surface dun objet par comparaison de deux images : une image de rfrence de
lobjet ltat initial et une image du mme objet aprs avoir subi une transformation
mcanique (dplacement ou dformation) (Sutton et al. [190]). En marquant la surface
dun motif rgulier (pour la mthode des grilles, Avril [19], Surrel [189]) ou alatoire (pour
la corrlation dimages), chaque point de la surface est identifi dans chaque image. En
corrlant limage de rfrence et limage de lobjet dform, nous obtenons une information
en deux dimensions des dplacements.
Linterfromtrie holographique Cette technique consiste faire apparatre des

franges dinterfrences en superposant des ondes lumineuses dont une au moins est pro-
duite par un hologramme(Franon [80], Rastogi [173], Smigielski [183, 184, 185]), un ho-
logramme tant un enregistrement par interfromtrie de londe lumineuse diffuse ou
diffracte par un objet.
Linterfromtrie de speckle Linterfromtrie de speckle permet lanalyse de franges

dinterfrences en superposant les intensits lumineuses diffuses ou diffractes par un
objet prsentant une surface rugueuse en lclairant par des faisceaux de lumire cohrente,
issues du mme laser, et orientes dans des directions diffrentes par un jeu de miroirs
(Laermann [120], Leendertz and Butters [132], Rastogi [173], Smigielski [186]).
4.2.2 Mesures surfaciques et hors-plan (2D 1/2)

Les mesures surfaciques et hors-plan permettent dobtenir un champ de dplacement
3D sur une surface. Dans un repre (O, x, y, z), les composantes du dplacement u, v et w
sont mesures. Elles sont fonction de x, y et z, cette dernire tant exprime en fonction
de x et y :

u(x, y, z(x, y))

v(x, y, z(x, y)) (4.4)

w(x, y, z(x, y))

Le gradient de dplacement dcrit par lquation 4.5 ne permet daccder quaux compo-
santes planes du tenseur de Green-Lagrange E (quation 4.6).
=

Fxx Fxy
F = Fyx Fyy (4.5)

Fzx Fzy

E xx E xy
1 T
E= F + F I = Exy Eyy (4.6)

= 2 = = =

Dans cette gamme, linterfromtrie, le moir et la stro-corrlation dimages figurent
parmi les mthodes les plus rpandues.
Linterfromtrie holographique et interfromtrie de speckle Ces techniques,

dj dcrites, permettent galement des mesures hors-plan (i.e. en augmentant le nombre
de direction dillumination).
Le moir dombre La technique du moir peut tre utilise pour dterminer le relief
dune pice. Dans ce cas, lombre de la grille de rfrence est projete sur la surface dun
objet par une source ponctuelle. La superposition du rseau et de son ombre montre les
lignes de niveau de lobjet par rapport au plan de rfrence dfini par le plan de la grille
(Breque et al. [48], Post et al. [167]).
Stro-corrlation dimages numriques La stro-corrlation permet deffectuer

des mesures de dplacement plan et hors-plan en combinant la corrlation dimages et
lutilisation de la stroscopie (Cornille [64], Garcia [82], Orteu [153]). Deux types dap-
pariements dimages sont raliss : un appariement temporel dune paire dimages issues
dune mme camra prise deux instants diffrents, et un appariement stroscopique
dune paire dimages prise au mme instant par deux camras diffrentes (figure 4.1).
Figure 4.1 Sytme possdant deux camras (Fazzini [76])
Lutilisation de deux camras en strovision permet en effet daccder une informa-

tion tridimensionnelle. Deux points sur la mme droite projective depuis lune des camras
peuvent ainsi tre distingus par la seconde camra (figure 4.2) et, de faon plus gn-
rale, on peut partir des deux images 2D dun mme point sous deux angles diffrents
dterminer la position tridimensionnelle de celui-ci par triangulation.
(a) une camra (b) deux camras
Figure 4.2 Principe stroscopique (Fazzini [76])

En enregistrant une paire dimages stroscopiques relatives chaque tat de dforma-

tion, la stro-corrlation dimages permet ainsi de mesurer : le champ de dplacements
tridimensionnel des points de sa surface tudie, et par voie de consquence, lvolution
de la gomtrie 3D dune surface (figure 4.3). De ce fait, elle permet aussi, lors de suivis
dessais, de mesurer les champs de dformations de lobjet inspect.
Figure 4.3 Images de gauche et droite et reconstruction 3D du cylindre (Fazzini [76])
Techniquement, le principe de la mthode est rsum la figure 4.4 dans le cas de

deux paires stroscopiques (une avant dformation et une aprs). A ltat initial, une paire
dimages stroscopiques est utilise pour reconstruire lobjet au repos par un appariement
stroscopique (tape 1). Une seconde paire dimages est prise lorsque lobjet est dform
et un appariement temporel est effectu sur les images de la camra de gauche ( ltat
initial et ltat dform) (tape 2). Enfin, un nouvel appariement stroscopique est
ralis sur les images de la camra de droite et la camra de gauche ltat dform
(tape 3).
Figure 4.4 Dtermination du champ de dplacement 3D par stro-corrlation (Fazzini

[76])
4.2.3 Mesures volumiques (3D)

Les mesures volumiques permettent dobtenir un champ de dplacement 3D dans le
volume du systme tudi. Dans un repre (O, x, y, z), les composantes du dplacement
u, v et w sont mesures et elles sont fonction de x, y et z :

u(x, y, z)

v(x, y, z) (4.7)

w(x, y, z)

Le gradient de dplacement F dcrit par lquation 4.8 permet daccder cette fois
lintgralit du tenseur de Green-Lagrange E (quation 4.9).
=

Fxx Fxy Fxz
F = Fyx Fyy Fyz (4.8)

Fzx Fzy Fzz

E xx E xy E xz
1 T
E= F + F I = Exy Eyy Eyz (4.9)

= 2 = = =
Exz Eyz Ezz
On dtaille ci-dessous deux ensembles de techniques de mesure volumiques, la photolas-
ticit et les mthodes utilisant lanalyse dimages.
4.3. Mthodes didentification du comportement 119
La photolasticimtrie 3D par dcoupage optique La photolasticimtrie 3D par

dcoupage optique consiste isoler un feuillet dune pice tridimensionnelle constitue
de matriaux birfringent par deux faisceaux laser plans. Un programme danalyse des
champs lumineux diffuss permet dobtenir une srie de franges isochromes et isoclines.
Cette image est similaire celle obtenue en plaant dans un polariscope rectiligne le mme
feuillet dcoup mcaniquement (Bilek et al. [38], Dupr and Lagarde [73], Germaneau
[86], Zenina [200], Zenina et al. [201]). Contrairement la photolasticimtrie 2D, cette
mthode lavantage dtre non destructive et permet des mesures en temps quasi-rel.
Corrlation volumique La corrlation volumique est une extension en trois dimen-

sions de la corrlation dimage numrique bidimensionnelle (pixels) utilisant des images
numriques de volumes (voxels) gnrs par tomographie rayons X (Bay et al. [31], Bornert
et al. [43], Germaneau [86]) ou des images obtenues par tomographie optique (Doumalin
et al. [70], Germaneau [86]) (mthode consistant isoler un plan du volume de la pice
par un faisceau laser et reconstituer le volume en regroupant les diffrentes coupes).
4.3 Mthodes didentification du comportement
Lidentification de paramtres pilotant un modle de comportement, partir dune

rponse mcanique mesure au cours dun essai statiquement indtermin, constitue un
problme inverse (Bui [49], Decultot [67]). Ce type de problme soppose au problme
direct qui consiste calculer la rponse dun systme mcanique, pour un chargement et
une loi de comportement donne (Latourte [123]) (tableau 4.1).
Le plus souvent, le problme direct est dit bien pos au sens de Hadamard : une solution
existe, elle est unique, elle dpend continment des donnes. Les problmes inverses sont
dit souvent mal poss car ne vrifiant pas au moins une des trois conditions de Hadamard
(Bonnet [40]).
Problme direct Problme indirect

Donnes de base Rsultats Donnes de base Rsultats
gomtrie initiale gomtrie dforme gomtrie initiale gomtrie dforme
maillage contraintes maillage contraintes
chargement dformations chargement dformations
lois de
conditions limites conditions limites
comportement
lois de un rsultat
comportement exprimental
paramtres des lois
Table 4.1 Donnes de base et rsultats des problmes directs et indirects (Kleinermann
[118])
Les mthodes didentification paramtrique peuvent globalement tre classes en deux

catgories :
les techniques itratives, pour la plupart inspires de la rsolution du problme
direct,
les techniques non itratives, en gnral spcifiques aux mesures de champs.
On se propose dans ce qui suit de raliser un inventaire des mthodes les plus courantes
de chacune de ces classes (Avril et al. [20], Grediac [93, 94]).Notons que la technique de
recalage par modles lments finis y sera davantage approfondie dans la mesure o elle
est utilise dans le cadre de ces travaux.
4.3.1 Mthodes itratives

Recalage de modles lments finis (FEMU) La mthode de recalage de modles
lments finis (Finite Element Model Updating method) est sans doute lapproche la plus
intuitive de toutes les mthodes didentification inventories. Elle consiste recaler de fa-
on itrative les paramtres dun modle Elments Finis (EF) afin dapprocher au mieux
des champs cinmatiques dits de rfrence. Dans lidal, la rfrence est bien entendu la
mesure exprimentale mais dans certains cas (mesures non disponibles, mise en place de
procdure) il est possible que les champs de rfrence soient donns par un modle num-
rique jug reprsentatif de la ralit : on parle alors de modle exprimental numrique
(Crouzeix [65], Fazzini [76], Lecompte [128], Pagnacco and Caro-Bretelle [154], figure 4.5).
Le but de la mthode est ainsi de minimiser lcart entre le modle EF identifier et la
rfrence, appel galement fonctionnelle cot ou fonction objectif. Sil existe diffrentes
fonctions cot et diffrentes mthodes de rsolution, la minimisation de la fonction se fait

gnralement par les algorithmes de Gauss-Newton ou de Levenberg-Marquardt (Avril
et al. [20]).
Figure 4.5 Principe de la mthode didentification FEMU
Grce notamment lutilisation des lments finis, cette mthode est flexible et peut
tre utilise avec un nombre limit de donnes (elle nest dailleurs pas restreinte aux
mesures de champs). Elle est ainsi largement applicable des pices de gomtrie et
chargement complexes. De nombreux travaux ont t raliss pour identifier les paramtres
lastiques ou lasto-plastique de matriaux isotropes (Avril and Pierron [21], Claire et al.
[58, 59], Crouzeix [65], Fazzini [76], Pannier [157]), les paramtres lastiques de matriaux
orthotropes (Leclerc et al. [127], Lecompte et al. [129], Molimard et al. [147], Rikards et al.
[175], Vautrin et al. [193]). Elle a galement t employe dans des cas de comportements
non-linaires tels que lendommagement (Geers et al. [84]), llasto-plasticit (Cooreman
[63], Fazzini [76], Fazzini et al. [77]) ou la viscoplasticit (Le-Magorou et al. [126]). On
notera que les calculs itratifs peuvent nanmoins conduire dans certains cas des temps
de calcul importants et quil convient de vrifier la stabilit/unicit de la solution par
rapport aux valeurs initiales.
Erreur en relation de comportement (CEGM) La mthode de lerreur en relation

de comportement (Constitutive Equation Gap Method) repose sur une fonctionnelle cot
construite partir des principes variationnels de la thorie de llasticit (Bonnet and
Pagano [42], Latourte [123], Latourte et al. [124]). Lobjectif est de minimiser une erreur
en relation de comportement entre le dplacement u mesur et la contrainte sur le
volume , dont lexpression est base sur lcriture de la somme des nergies potentielles
et complmentaires :
Z
1
E(u, , C) = ( C : (u)) : C 1 : ( C : (u))dv (4.10)
2
avec C le tenseur dlasticit recherch. Daprs les principes variationnels, cette erreur
E est positive pour tout champ de dplacement u cinmatiquement admissible et tenseur
de contraintes statiquement admissible et nulle pour un problme de mcanique des
solides bien pos. Dans la pratique, avec des champs mesurs, cette minimisation peut
tre partielle (Bonnet and Pagano [42]).
Cette technique, qui peut galement fonctionner avec un nombre limit de mesures, a
t mise en place pour lidentification de comportements lastiques linaires ou non ainsi
que de comportements lastoplastiques (Bonnet et al. [41], Geymonat et al. [87], Latourte
[123], Merzouki et al. [142]).
Ecart la rciprocit (RGM) La mthode dcart la rciprocit (Reciprocity Gap

Method) consiste minimiser une fonctionnelle cot dfinie partir du Principe des
Puissances Virtuelles (PPV) ou partir du thorme de Maxwell-Betti (Andrieux et al.
[16], Andrieux and Constantinescu [17], Bonnet et al. [41]). Lide est alors dcrire ce
principe entre deux tats, lun issu de la mesure u (,,C) et lautre adjoint u ( ,u ,
C ) :
Z Z

: [C C] : dV = (T u T u ) dS (4.11)

avec T le vecteur contrainte sur la frontire extrieure au domaine. Cette mthode ne

requiert ainsi que des donnes de bord (supposes connues dans le problme pos).
Concernant les applications, la mthode de lcart la rciprocit est principalement

utilise pour lidentification de fissures (Andrieux et al. [16], Andrieux and Constantinescu
[17], Bonnet et al. [41]) et de champs de modules lastiques ou de rigidits de flexion
htrognes. Il convient de souligner quil ny a pas lheure actuelle de mise en oeuvre
de cette dmarche partir de donnes exprimentales.
4.3.2 Mthodes non itratives

Mthode des champs virtuels (VFM) Introduite par Grediac [91], la mthode des
champs virtuels (Virtual Fields Method) repose sur lcriture de lquilibre global de
lprouvette teste avec le Principe des Travaux Virtuels (PTV). En labsence de forces
volumiques extrieures et pour un chargement surfacique T sur la surface extrieure
dun volume , ce principe scrit :
Z Z

: d = T u dS (4.12)

pour tout champ de dplacement virtuel cinmatiquement admissible u ( la dfor-
mation associe, Grediac [94], Grediac et al. [95, 98, 99], Meraghni et al. [141], Pannier
[157]).
Le principe de la mthode est de construire un ensemble dquations en introduisant

plusieurs champs virtuels judicieusement choisis de manire faire apparatre les para-
mtres de la loi de comportement = f () (Grediac [91], Grediac et al. [95], Pannier
[157]). Une procdure a t propose par Grediac et al. [99] pour gnrer ces champs
virtuels spciaux dont le choix (nombre et type) dpend de la loi de comportement d-
terminer et va conditionner le degr dindpendance des quations du systme dvelopp.
En particulier, sil y a au moins autant de champs virtuels spciaux que de paramtres
dterminer et si ces paramtres sont crits sous forme polynomiale, le systme dquations
sera linaire et les inconnues seront dtermines directement aprs inversion du systme.
La mthode des champs virtuels permet didentifier des paramtres de lois de com-
portement de matriaux partir de champs cinmatiques htrognes. Elle a t appli-
que lidentification de rigidits de flexion de plaque en statique (Grediac and Vautrin
[100], Grediac [92]) et en dynamique (Grediac and Paris [96], Guo [103]). Elle a gale-
ment permis lidentification dans le plan de proprits lastiques de matriaux (Pierron
et al. [166]). Son emploi reste cependant assez rare pour les comportements non linaires
(Grediac and Pierron [97]) et la construction des champs virtuels spciaux peut sav-
rer trs dlicate dans des cas plus complexes (matriau anisotrope, gomtrie complexe,
chargement multiaxial).
Mthode de lcart lquilibre (EGM) Lobjectif de la mthode de lcart lqui-

libre (Equilibrium Gap Method) est de dterminer un champ de proprits mcaniques
satisfaisant, pour un champ de dplacement donn, lquilibre interne (div() = 0) de
la structure (Amiot et al. [15], Claire et al. [58, 59], Crouzeix [65]). La discrtisation des
quations dquilibre, crites sous une forme faible (entre des lments adjacents pour une
discrtisation lments finis), conduit un systme dquations linaires dont les donnes
sont les champs de dplacement mesurs et les inconnues les champs de modules lastiques.
Les rsultats de cette mthode sont des cartographies de proprits mcaniques.
En pratique, la technique de lcart lquilibre consiste minimiser les rsidus def-

forts nodaux associs aux contrastes de proprits lastiques entre lments voisins (Amiot
et al. [15], Perie et al. [164, 163]) :
{fres } = [K({d})] {umes } {fvol } (4.13)
avec {fres } le vecteur des rsidus defforts nodaux associs au champ de dplacement
mesur {umes }, {d} le champ dendommagement dterminer, [K({d})] la matrice de

rigidit et {fvol } les forces nodales lies au chargement extrieur. Cette mthodologie
a par exemple t mise en place pour lidentification de modules de cisaillement dun
composite orthotrope htrogne (essai sur prouvette biaxiale) (Crouzeix [65]).
Une des grandes limites de cette mthode est lie au maillage sur lequel sapplique la
discrtisation de la formulation qui doit tre bas directement sur les points de mesure.
Cela ne permet pas ltude de gomtries complexes et justifie le recours gnralement
des champs simuls (tests virtuels) comme champs mesurs (Claire et al. [58, 59], Crouzeix
[65], Perie et al. [163, 164], Roux and Hild [178]).
4.4. Mthodologie didentification multi-zones 125
4.4 Mthodologie didentification multi-zones
Comme cela a t voqu prcdemment, la mthode FEMU a t implmente au

sein du laboratoire en lien avec le dveloppement des mthodes optiques de mesures de
champs cinmatiques. Cette dmarche didentification a ainsi t mise en application sur le
cas dun alliage daluminium lastoplastique isotrope partir de mesures exprimentales
ralises par stro-corrlation dimages. On pourra trouver des informations dtailles
ce sujet dans la thse de Fazzini [76].
Vis--vis des objectifs lis ce travail, nous nous sommes donc naturellement appuys
sur le savoir-faire et une partie des dveloppements internes dans ce domaine. La tech-
nique didentification FEMU a en particulier t privilgie, et ce dautant plus que le
contexte de multi-zones limite une grande partie des autres mthodes didentification. Si
la mthode des champs virtuels a t rcemment applique au cas dune plaque composite
endommage (Kim et al. [117]), un certain nombre de problmes de stabilit se posent
encore vis--vis du choix des fonctions polynomiales introduites dans la formulation (Vau-
trin et al. [193]). Le logiciel de calcul par lments finis employ est quant lui Abaqus
(Simulia [182]).
Lenjeu de ce chapitre est donc de prsenter la dmarche didentification permettant

didentifier la variation spatiale des proprits lastiques macroscopiques dune structure
composite. Prcisment, le problme considr ici fait le lien avec le premier axe de la
thse ddi au CND au travers de ltude dune pice comportant deux zones :
lune saine constituant majoritairement la structure,

lautre dgrade reprsentant un dfaut interne localis (figure 4.6).
Les donnes dentre de ce problme sont les champs cinmatiques induits sur cette pice
lors dun chargement mcanique, qui prsentent une htrognit macroscopique spcifi-
quement induite par la prsence du dfaut. A partir de ces informations, la mthodologie
didentification vise alors :
dlimiter spatialement les deux zones,

dterminer pour chacune dentre elles les proprits lastiques effectives du matriau.
Figure 4.6 Problme considr
A linstar de Vautrin et al. [193], la stratgie employe pour rsoudre cette problma-
tique repose sur une analyse globale puis locale qui sarticule autour des grandes tapes
suivantes (figure 4.7) :
1. Identification globale des paramtres lastiques du matriau sain (macroscopique-
ment homogne) sur une prouvette talon,
2. Analyse des carts entre la rfrence et le modle EF de recalage constitu dun
matriau macroscopiquement homogne et sain puis localisation de la zone dgrade
par seuillage au niveau des carts significatifs,
3. Intgration de la zone dgrade dans le modle EF de recalage (matriau macrosco-
piquement htrogne comportant deux zones de proprits distinctes),
4. Identification locale des paramtres lastiques de la zone dgrade.
La dmarche propose ici se distingue cependant des travaux de Vautrin et al. [193] par
ltude dune structure htrogne pleine (ce qui rend lidentification globalement plus
difficile, voir Molimard et al. [146, 147]), dont la zone dgrade est constitu dun matriau
de mme classe (do un constraste de proprits plus faible) et orthotrope.
Figure 4.7 Stratgie didentification multi-zones

4.4. Mthodologie didentification multi-zones 127
Prcisment, le cas industriel faisant lobjet de cette tude concerne un composite

stratifi monolithique en tissu carbone Hexcel HexForce G0926 associs la rsine RTM6
possdant la squence dempilement [45/0/45]2 . La prsence dun dfaut dans cette
structure saine est prise en compte au travers dune singularit gomtrique, savoir un
trou fond plat de 6 mm de diamtre situ en face arrire de cette plaque et faisant 3 plis
dpaisseur. La procdure de recalage concerne une zone dintrt carre situe autour du
trou (voir la figure 4.8-a) et porte sur les donnes de champs cinmatiques obtenus en face
avant qui intgrent de manire continue (la face avant tant plane) lexistence de ce dfaut.
Notons que la taille du dfaut considr a t dicte par les exigences aronautiques ; sa
profondeur a quant elle t choisie de manire obtenir une htrognit suffisamment
significative macroscopiquement sur les dplacements de la face avant tout en conservant
le type de symtrie matrielle du matriau sain (notamment les axes principaux) au niveau
de la zone du dfaut.
(a)
(b)
Figure 4.8 Structure multi-zones de rfrence au niveau de la zone dintrt : (a)

principe, (b) modle exprimental numrique associ (vue de la face arrire, les zooms
montrent les diffrents plis de chaque zone)
Comme cela a t voqu prcdemment, la mthodologie classique par FEMU se base

sur des rsultats exprimentaux pour constituer la rfrence de lidentification. Dans ce
travail, nous avons adopt une dmarche progressive afin de sextraire en premire ap-
proche dun certain nombre de difficults lies aux bruits de mesures et aux discontinuits
locales des champs exprimentaux. Ainsi, la rfrence consiste ici en un jeu de donnes
simules numriquement (modle exprimental numrique), ce qui permet de construire
entirement la procdure didentification multi-zones. Afin de sapprocher au plus prs de
la ralit du comportement du composite tudi, la rfrence sera tablie laide dune
modlisation par lments finis qui tient compte de lempilement des diffrents plis (pro-
prits et orientation, figure 4.8-b).
Lensemble de la dmarche tant pose, on se propose prsent de dtailler les diff-

rentes tapes de sa mise en place.
4.5 Identification des proprits dun matriau ortho-

trope plan
On sintresse dans cette partie lidentification des proprits dun matriau ma-
croscopiquement homogne lastique et anisotrope (tape 1 de la stratgie dcrite la
figure 4.7). Prcisment, il correspond dans le problme tudi au composite stratifi
G0926/RTM6 sain (sans dfaut) qui prsente une symtrie matrielle de type orthotropie
plane.
En rgle gnrale, cette lasticit anisotrope est dcrite par quatre paramtres las-
tiques plans ainsi que par lorientation des axes dorthotropie. Le tissu de carbone utilis
dans le cas dtude est parfaitement quilibr et les orientations des plis considres pour
le stratifi (matriau FUSCOMP et autres validations) sont par la suite exclusivement
0, 45 et 90 ; le repre principal dorthotropie (e1 , e2 , e3 ) correspond alors au repre de
calcul (x, y, z) (ces deux bases tant orthonormes directes, figure 4.9). Lidentification
revient donc dterminer les quatre coefficients lastiques plans suivants (la direction
z = e3 tant perpendiculaire la surface du stratifi) :
les modules dYoung E1 et E2 , associs respectivement aux directions de vecteur
unitaire e1 et e2 ,
le coefficient de Poisson 12 , associ au couple de directions (e1 , e2 ),
le module de cisaillement G12 , associ au couple de directions (e1 , e2 ).
4.5. Identification des proprits dun matriau orthotrope plan 129
Figure 4.9 Orientation de la symtrie matrielle du matriau orthotrope
4.5.1 Modles de rfrence et de recalage

Lidentification des quatres paramtres prcits requiert la gnration dun champ
de contraintes htrognes. Selon (Molimard et al. [146, 147]), lutilisation des champs
cinmatiques induits sur la surface dune prouvette munie dun trou dbouchant lors dun
essai de traction uniaxiale permet lobtention de ces proprits. Dans le cas prsent, les
axes principaux dans le plan de lprouvette de traction correspondent (x = e1 , y = e2 )
et lessai est ralis selon laxe y = e2 (figure 4.9).
La structure de rfrence considre pour cette premire tape est ainsi une prouvette
de traction troue en son centre, dont on dcrit chaque pli laide dune modlisation
3D solide stratifi dans Abaqus. La procdure didentification est illustre par la suite sur
plusieurs cas de matriaux composites bass sur lemploi du mme tissu carbone Hexcel
HexForce G0926 mais qui se distinguent vis--vis de lorientation de leurs plis :
le composite FUSCOMP 6 plis, dempilement [45/0/45]2 ,
un composite not COMP-A 4 plis, dempilement [0/90]s ,
un composite not COMP-B 4 plis, dempilement [45/ 45]s .
Pour chaque pli i de ces structures, les donnes prciser sont donc :
son orientation dfinie par langle i = (x, x0i ) entre le repre de calcul (x, y, z) et le
repre principal du tissu (x0i , yi0 , z) (cette base est galement orthonorme directe),
les proprits principales planes du tissu : les modules dYoung E10 et E20 (E10 = E20
puisque le tissu est quilibr), le coefficient de Poisson 120
et le module de cisaillement
G12 (tableau 4.2).
0
Pour tous ces matriaux, chaque pli est discrtis dans lpaisseur en cinq lments finis de
faon sapprocher au mieux du comportement rel du composite et notamment assurer
la continuit des contraintes de cisaillement transverse aux interfaces entre les plis (figure
4.6-a). Les lments finis utiliss sont des ttradres C3D4 et le modle comporte en tout
470219 lments soit 1880876 noeuds dont 2399 noeuds sur la surface de la zone dintrt.
Proprits lastiques principales du pli

modules dYoung E10 = E20 = 64000 M P a
coefficient de Poisson 0
12 = 0.05
module de cisaillement 0
G12 = 4500 M P a
Table 4.2 Proprits planes du tissu Hexcel HexForce G0926 dans ses axes principaux
(a) rfrence (matriau FUSCOMP) (b) modle EF de recalage

Figure 4.10 Identification des proprits orthotropes - Zone dintrt
Le modle EF de recalage sert quant lui identifier les proprits macroscopiques du

matriau. Il est galement trou en son centre mais construit laide dune modlisation
3D solide classique, cest--dire homogne dans lpaisseur (figure 4.6-b). Deux lments
finis ont galement t placs dans lpaisseur pour assurer la stabilit numrique du
calcul. Les lments finis utiliss sont des ttradres C3D4 et le modle comporte 168891
lments soit 675564 noeuds dont 2400 noeuds sur la surface de la zone dintrt.
Il est enfin important de signaler deux points importants de la dmarche propose :
dune part, les maillages utiliss pour chacun des deux modles (rfrence exp-
rimental numrique et recalage) sont diffrents ; la plus grande gnralit a t
conserve ce niveau afin de prparer dores et dj lextension une rfrence
vritablement exprimentale ;
lidentification a t ralise ici sur une zone dintrt comportant prs de 2400
points (voir figure 4.6), soit en dautres termes en comparant les dplacements de
2400 noeuds en surface du modle EF de recalage avec leurs homologues de la
rfrence (voir bloc PtVoisin ci-aprs) ; par rapport lalgorithme utilis par
Fazzini [76], le nombre de points danalyse a ainsi t fortement augment.
4.5.2 Description du fonctionnement de lalgorithme didentifica-

tion
Le script du programme didentification, lanc partir du logiciel Abaqus et ralis
en langage Python, comporte douze blocs fonctions subdiviss en deux familles de pro-
gramme :
le programme principal FonctionnelleCout qui fait appel lensemble des autres
modules,
les onze autres programmes subordonns au module principal (figure 4.11).
Outre les fonctions classiques dune procdure didentification FEMU, ce programme com-
porte un ensemble dactions spcifiques permettant la comparaison des diffrents maillages
(de la rfrence et du modle EF de reclage).
Figure 4.11 Schma de fonctionnement de lalgorithme didentification
4.5.2.1 Programme principal FonctionnelleCout
Le programme principal est employ chaque itration par la fonction scipy . opti-
mize . leastsq qui est un algorithme doptimisation de Levenberg-Marquardt disponible
dans la librairie Python. Cette partie du programme rgit ensuite lenchanement des onze
autres modules : elle fait appel eux pour leur fonction individuelle et rcupre dans une
variable la modification apporte. Le module FonctionnelleCout renvoie la fonction
cot ncessaire lalgorithme doptimisation. Dans le cas prsent, la fonction cot est
dfinie par : v
u l m n
uX X X
exp 2
J(p) = t unum
ij (p) uij (4.14)
k=1 i=1 j=1
o p est le vecteur des valeurs des paramtres. Le couple (i ; j) correspond aux coordonnes
suivant x et y et le paramtre k indique le step.
4.5.2.2 Programme subordonns
Bloc modifModele Ce premier bloc permet de modifier directement les para-

mtres de la loi de comportement dun matriau dans le fichier modle EF qui est utilis
chaque itration. Dans notre cas, ce bloc actualise donc les modules dYoung E1 et E2 ,
le coefficient de Poisson 12 et le module de cisaillement G12 .
Blocs PointsSimulation et LectureExp Le premier bloc rcupre les co-

ordonnes dans le repre (x, y, z) des points du maillage du modle EF de recalage pour
chaque pas de calcul, le second ralisant une opration similaire sur lobjet rfrence. Ces
points sont extraits dans un tableau n lignes et 6 colonnes, avec n le nombre de points
du maillage : les trois premires colonnes reprsentent les coordonnes des points ltat
non dform, les trois colonnes suivantes correspondant au cas dform.
Blocs Normaliser et Denormaliser La fonction Normaliser permet de

dfinir une variation des paramtres de la loi de comportement comprise entre +1 (borne
haute) et -1 (borne basse) des paramtres. Cette normalisation permet un travail plus
efficace de la fonction doptimisation.
La fonction Denormaliser effectue bien entendu le travail inverse et fournit la fin
de chaque itration de calcul du programme didentification, les valeurs des paramtres
de la loi de comportement.
Bloc FaireLeMenage Cette fonction permet de supprimer les fichiers temporaires

crs par le calcul lments finis.
Les deux blocs suivants permettent dappliquer le processus itratif didentification

sur les mmes points physiques de la structure, bien que les maillages de rfrence et de
recalage soient distincts.
Bloc PtVoisin Ce bloc calcule la distance entre un point du maillage de rfrence

et un point du maillage numrique. Les trois points du maillage numrique les plus proches
dun point exprimental sont reprs. Notons que le maillage EF ttradrique permet
dviter les configurations dalignements de points (amenant des aires dinterpolation
quasi nulles et des divisions par zro dans le calcul dinterpolation).
Bloc InterpolationNoeud Cette partie du programme permet, laide dune

fonction dinterpolation, dvaluer les coordonnes dun point en fonction des trois points
voisins trouvs dans le bloc prcdant. Ces points sont situs sur le maillage non dform
et proches du point de la rfrence concern. Aprs cette tape, une matrice de points
identiques aux points de rfrence est obtenue possdant les coordonnes dformes du
modle numrique de recalage.
Bloc CoordInterpoles Cette fonction permet de prparer les trois variables

matricielles ncessaires au calcul du point numrique cr par linterpolation. Une des
matrices correspond aux coordonnes dun point de la rfrence ltat initial, une autre
aux coordonnes des trois points voisins du modle EF de recalage ltat initial et une
troisime aux coordonnes des trois points voisins de ce dernier ltat dform.
Bloc Ecart Ce bloc calcule en chaque point le vecteur cart entre les points
numriques interpols et les points de rfrence.
Bloc CoupTorchon Cette fonction permet de supprimer tous les points aberrants
lis une erreur de la fonction dinterpolation.
4.5.3 Rsultats
Le tableau 4.3 indique les rsultats issus de la procdure didentification pour les
trois matriaux orthotropes homognes considrs. Les valeurs des paramtres issus de la
thorie classique des stratifis pour chaque empilement sont galement prciss. Pour un
temps de calcul CPU relativement court (environ 30 minutes avec processeur Intel Xeon
E5410 2,33 GHz 4 coeurs), on observe donc une bonne cohrence des rsultats obtenus
par recalage vis--vis de cette reprsentation.
Module Module Module de

Coefficient de
dYoung E1 dYoung E2 cisaillement
Poisson 12
[MPa] [MPa] G12 [MPa]
Matriau 37224 37653 0.4406 19871
FUSCOMP (37849 ) (37849 ) (0.438 ) (21817 )
Matriau 63707 64006 0.0496 4494
COMP-A (64000 ) (64000 ) (0.05 0) (4500 )
Matriau 15885 15877 0.7647 30241
COMP-B (15878 ) (15878 ) (0.764 ) (30476 )
Table 4.3 Identification de matriaux orthotropes par procdure FEMU (le calcul de
chaque paramtre selon la thorie classique des stratifis est donn entre parenthses)
Afin de dmontrer la stabilit de la procdure dveloppe, les volutions des para-

mtres ainsi que de la fonction cot avec le nombre ditrations sont indiques pour le
cas du matriau FUSCOMP (figure 4.12). Les diffrentes courbes se distinguent par les
valeurs initiales employes pour les paramtres. Pour les quatre jeux de dpart employs,
on note une convergence assez rapide vers des valeurs quasi-identiques, une valeur de d-
part loigne ncessitant bien entendu davantage ditrations. Egalement, la fonction cot
obtenue est dans tous les cas proche de zro, ce qui dmontre la bonne convergence de
lalgorithme.
4.6. Identification des proprits dune zone dgrade 135
(a) Module dYoung E1 [MPa] (b) Module dYoung E2 [MPa]
(c) Coefficient de poisson 12 (d) Module de cisaillement G12 [MPa]
(e) Fonction cot

Figure 4.12 Courbes de convergence du cas FUSCOMP
4.6 Identification des proprits dune zone dgrade

Lidentification dun matriau orthotrope macroscopiquement homogne tant ralise,
revenons prsent au problme multi-zones de dpart (figure 4.6). A ce stade, les donnes
considres sont donc :
la rfrence, constitue par le modle exprimental numrique dcrit la figure
4.8, cest--dire un composite (matriau FUSCOMP) comportant un trou fond
plat en face arrire,
le modle EF de recalage dpaisseur constante (plein), dont on cherche identifier

les paramtres dans la zone dgrade pour se rapprocher au plus prs des dplace-
ments observs en face avant de la rfrence.
Lobjectif de cette partie est donc de traiter les tapes 2 4 de la stratgie didentification
multi-zones dcrite la figure 4.7. La plus grande difficult dans ce processus rside en
ltape 2 qui consiste en la localisation de la zone dgrade. Elle requiert en premier lieu
la quantification dune erreur pour distinguer le matriau sain du matriau dgrad puis
la dfinition dun seuil pour prciser sa position spatiale.
4.6.1 Cration de la cartographie derreurs

Les proprits du matriau sain tant connues (tape 1), on cherche mettre en
vidence lcart entre celui-ci et le matriau dgrad. Pour cela, on compare pour des
prouvettes daxes (x = e1 , y = e2 ) soumises un essai de traction simple daxe y = e2 :
la rfrence dcrite ci-dessus,
avec le modle EF de recalage plein, auquel on affecte des proprits macroscopiques
homognes (mono-zone) correspondant au matriau sain.
Compte tenu du dfaut dans la structure de rfrence, il y aura par consquent une
diffrence de rponse entre les deux.

Figure 4.13 Calcul derreur - Zone dintrt
A linstar de la dmarche pour matriaux orthotropes, la structure de rfrence intgre

la stratification du composite FUSCOMP laide dune modlisation 3D solide stratifi
dans Abaqus ; la diffrence porte uniquement sur la prsence du trou fond plat en
face arrire (6 mm de diamtre, profondeur correspondant trois plis, figure 4.13-a). Le
maillage prsente galement les mmes spcificits (5 lments dans lpaisseur de chaque
pli, lments ttradriques C3D4, 545221 lments soit 2180884 noeuds dont 2800 noeuds
en face avant de la zone dintrt).
Le modle EF de recalage suit lui aussi les mmes prrogatives quau paragraphe
4.5.1, avec une modlisation 3D solide mais cette fois dpaisseur constante et continue
(structure pleine). Son maillage est conforme au cas 4.5.1 (2 lments dans lpaisseur,
lments ttradriques C3D4, 174860 lments soit 699440 noeuds dont 2800 noeuds sur
la surface de la zone dintrt). Les proprits saines affectes au matriau qui le constitue
sont celles identifies dans la partie prcdente pour le matriau FUSCOMP (tableau 4.3).
La comparaison des rponses fournies par la rfrence et le modle EF de recalage lors

de lessai de traction est ralise sur une zone dintrt comportant plus de 2800 points.
Prcisment, deux options sont a priori possibles quant lobjet de cette comparaison :
le champ de dplacements ,
le champ de dformations U ,
Afin de guider ce choix, on prsente la figure 4.14, les volutions des composantes
de ces diffrents champs sur la face avant de la rfrence, selon une ligne traversant
les zones saine et dgrade. On observe ainsi que les champs de dformation mettent
beaucoup plus en vidence la singularit gomtrique et offrent ainsi la rponse la plus
reprsentative de la prsence du dfaut. Le calcul portera donc sur cette information.
Dautre part, lensemble des composantes du tenseur de dformation tant affectes, il
semble pertinent pour intgrer de manire globale tous ces phnomnes de recourir la
dformation quivalente au sens de Von Mises dans le calcul derreur :
r
5 2 5 4
equi = xx + 2yy xx yy + 32xy (4.15)
6 6 3
(a) Ligne de mesure AB
(b) Dplacement Ux (c) Dplacement Uy
(d) Dformation xx (e) Dformation yy
(f) Dformation xy
Figure 4.14 Evolution des composantes (Ux , Uy ) du champ de dplacement U et des

composantes (xx , yy , xy ) du champ de dformations dans la base (x, y) en face avant
de la structure de rfrence le long dune ligne AB traversant les zones saine et dgrade
Un algorithme spcifique en langage Python a ainsi t dvelopp pour raliser au-

tomatiquement la cartographie derreur entre le modle de rfrence et le modle EF de
recalage suppos homogne et sain. Ses principales tapes sont les suivantes :
extraction pour les deux structures, des coordonnes de chaque point de leur maillage
ainsi que des trois composantes planes du champ de dformation (11 , 22 , 12 ) en
celui-ci,
calcul de la dformation quivalente au sens de Von Mises en chacun des points des
deux structures,
interpolation des points du maillage numrique sur celui de la structure de rfrence

de faon comparer les mmes points gomtriques de chacune des structures (voir
paragraphe 4.5.2),
calcul dcart entre les deux modles :
ref modele EF
equi equi
ecart = (4.16)
ref
equi
La figure 4.15 illustre une cartographie derreur obtenue selon cette dmarche. Notons
quelle est ralise en 256 niveaux de gris de faon faciliter les traitements dimages par
la suite.
Figure 4.15 Cartographie derreur en dformation (rfrence vs matriau monozone

sain)
4.6.2 Estimation de la taille de la zone dgrade

A partir de la carte derreur, on trace les profils de niveaux de gris obtenus selon deux
coupes correspondant aux directions principales du problme (figure 4.16), le long de la
direction axiale de traction (x = 0) et le long de la direction transversale (y = 0) ; ce
traitement dimage est du mme type que celui utilis dans le chapitre 3 pour lanalyse
CND des dfauts.
On dfinit les hauteurs relles h des profils entre le point haut de la courbe dter-
min automatiquement et la valeur basse identifie manuellement. Reste dfinir le seuil
qui permet destimer la taille de la zone dgrade. Plusieurs choix sont possibles, point
dinflexion, mi-hauteur, ... il a t dcid dans une premire approche de se baser sur le
seuil reprsentatif du diamtre du dfaut. Une zone affecte par le dfaut de proprits
moyennes entre la zone saine et la zone dgrade na pas t prise en compte dans cette
premire approche. Les seuils identifis correspondant au dfaut de 6 mm de diamtre

dans lprouvette sont reprsents en pointills rouges sur les figures 4.16-b et c. Sur la
coupe 4.16-c, le seuil tabli par rapport la taille relle du dfaut reprsente ainsi 67% de
la hauteur relle du profil alors que sur la coupe 4.16-b, il correspond 73%. La mesure
derreur tant effectue sur ltat dform de la pice, cet effet est induit par la traction
effectue selon laxe y qui a tendance accentuer lovalisation du trou fond plat dans
ce sens.
De manire simplifier la dmarche de localisation et de dimensionnement du dfaut,

un niveau mdian entre les deux tats est alors considr qui correspond 70% de la
hauteur de chacune des coupes (en pointills verts sur les figures 4.16-b et c). Ceci permet
dobtenir les coordonnes de quatre points (E, F sur la coupe x = 0 ; G et H sur la coupe
y = 0) par lesquels passent simultanment les deux profils. La forme du dfaut tant
approxime par des fonctions polynomiales, lobtention de ces quatre points permet donc
ce stade dtablir la disposition spatiale de la zone dgrade dans la pice (figure 4.17).
(a)
(b) Coupe sur x = 0 (c) Coupe sur y = 0
Figure 4.16 Identification de la taille de la zone dgrade : (a) coupes ralises sur la
carte derreur, (b) et (c) : volution des niveaux de gris et dfinition des seuils
Figure 4.17 Intgration de la zone dgrade
4.6.3 Identification des proprits de la zone dgrade

Dans cette dernire tape, la structure de rfrence est la mme quau paragraphe
4.6.1 mais il convient de modifier le modle EF de recalage de sorte que (figure 4.18) :
deux zones de proprits distinctes soient dclares selon la position spatiale issue
de ltude prcdente (figure 4.17),
les proprits identifies pour le matriau FUSCOMP (tableau 4.3) soient affectes
la zone saine,
les proprits identifier (quatre paramtres apparents plans lastiques orthotropes)
soient celles de la zone dgrade.
Les paramtres de ce nouveau modle numrique sont analogues aux tudes prcdentes
(modle 3D solide comportant deux matriaux, lments ttradriques C3D4, 178562 l-
ments soit 714248 noeuds dont un peu plus que 2800 noeuds en surface de la zone dint-
rt). Notons simplement que la zone dgrade comporte 400 points alors que la zone saine
en comprend 2400. La comparaison des rponses fournies par la rfrence et le modle EF
de recalage lors de lessai de traction est donc ralise sur une zone dintrt comportant
plus de 2800 points.

Figure 4.18 Identification des proprits dgrades - Zone dintrt
Deux jeux de dpart, extrmes, ont t employs pour lidentification de la zone d-

grade :
jeu 1 : les valeurs du matriau FUSCOMP sain (tableau 4.3),
jeu 2 : les valeurs des paramtres calcules par la thorie classique des stratifis pour
un empilement [45/0/45] (tableau 4.4).
Compte tenu du type de dfaut considr (trou fond plat), les valeurs didentification de
la zone dgrade devraient naturellement sapprocher de celles obtenues par la thorie des
stratifis pour les 3 plis. A la diffrence de lerreur entre la rfrence et le modle EF o
les proprits de la zone dgrade correspondent au premier jeu de dpart (figure 4.15),
lerreur entre la rfrence et le modle EF avec le second jeu comme proprits pour la
zone dgrade est bien quasi nulle sur lessentiel de la zone dintrt (figure 4.19). Les
paramtres issus de la thorie des stratifis constituent donc bien des valeurs cohrentes
pour le rsultat de ce calcul.
Module Module Module de

Coefficient de
dYoung E1 dYoung E2 cisaillement
Poisson 12
[MPa] [MPa] G12 [MPa]
18925 18925 0.438 10909
Table 4.4 Thorie des stratifis : empilement [45/0/45]
Figure 4.19 Cartographie derreur en dformation (rfrence vs structure multizones

base sur la thorie de stratifis)
Nanmoins, lvolution des valeurs des diffrents paramtres identifis en fonction du

nombre ditrations dmontre la non convergence de la rponse didentification selon le
jeu de paramtres de dpart choisi (figure 4.20). Avec le jeu 1, les paramtres restent en
gnral stables, alors quils voluent avec le jeu 2 ; dans le mme temps, les valeurs des
fonctions cot respectives sont faibles.
Ces constats nous incitent penser plusieurs origines potentielles pour ces pro-
blmes. Compte tenu des discontinuits de matire et de contraintes induites par le trou
fond plat dans la rfrence, il est certain tout dabord quune reprsentation bi-zone telle
que suppose la figure 4.8-a constitue une simplification forte qui perturbe sans doute
lidentification paramtrique. Celle-ci prend en compte en effet tous les points de la zone
dintrt de manire identique quils soient dans la zone saine, dans la zone centrale de la
zone dgrade et dans la zone de transition entre les deux. Cela rend donc la procdure
de minimisation et de convergence plus difficile. Par ailleurs, les soucis de localisation du
dfaut voqus prcdemment accentuent dautant plus ce phnomne. Une des solutions
envisages viserait slectionner des points reprsentatifs des proprits recherches,
soit parmi la zone centrale de la zone dgrade (la figure 4.19 confirme bien ladquation
des paramtres issus de la thorie des stratifis dans ce secteur).
Un autre point important concerne la sollicitation applique lprouvette. En lien

avec le reste du travail, nous avons considr dans cette dernire partie un essai de trac-
tion unixiale. Si ce chargement permet en effet de solliciter tous les modules lastiques
plans dune prouvette troue (Molimard et al. [146, 147]), il nest pas forcment le plus
adapt lidentification dune prouvette pleine. Ds lors, pour ce cas macroscopiquement
htrogne, il serait peut-tre prfrable denvisager une (ou plusieurs) sollicitation(s) de
traction bi-axiale qui permettrait de vritablement mettre en oeuvre diffrents modes de
dformation (voir par exemple Lecompte et al. [129]).
(a) Module dYoung E1 [MPa] (b) Module dYoung E2 [MPa]
(c) Coefficient de Poisson 12 (d) Module de cisaillement G12 [MPa]
(e) Fonction cot

Figure 4.20 Courbes de convergence de lidentification multi-zone
4.7 Conclusion
Au travers ce chapitre, une procdure didentification paramtrique FEMU sur struc-
ture composite orthotrope multi-zones a t mise en place. Un certain nombre daspects
importants ont pu tre valids, notamment :
lidentification des paramtres lastiques orthotropes plans dun matriau sain,
4.7. Conclusion 145
la dmarche de localisation et de dimensionnement dune zone dgrade prsentant

des proprits apparentes distinctes.
Ces deux tapes sappuient sur un nombre important de points danalyse et tendent
les travaux existants au contexte anisotrope. Le second point est de plus ralis sur une
prouvette pleine, permettant dtablir long terme une dmarche analogue au CND
partir de champs cinmatiques. La dernire tape concernant lidentification des proprits
de la zone dgrade nest pas compltement aboutie lheure actuelle, mais un certain
nombre de pistes de travail ont t voques afin de rsoudre les difficults de convergence.
Notons enfin que les perspectives majeures de ce travail porteront sur lidentification de
matriaux orthotropes dorientation arbitraire (ce qui supposera lextension la prise
en compte de langle de symtrie matrielle) ainsi qu la prise en compte de champs
cinmatiques exprimentaux comme rfrence de la procdure FEMU (issues notamment
de mesures par stro-corrlation dimages).
Conclusion
Ce travail de recherche a t consacr lidentification de dfauts au sein de structures

composites employes dans le secteur de laronautique. Outre des matriaux fabriqus
classiquement par fibres pr-imprgnes, cette dmarche a port galement sur des com-
posites raliss par infusion de rsine liquide, procd innovant faisant lobjet du projet
FUSCOMP.
Ltude prsente dans ce mmoire a contribu la dfinition de mthodologies per-

mettant la dtection et la caractrisation de dfauts internes, aujourdhui indispensables
pour le dveloppement et la maintenance des arostructures. Lide gnrale a t de
sappuyer sur des donnes de champs qui permettent une analyse globale et gnralement
sans contact de panneaux de fuselage. Prcisment, deux aspects ont t abords ici :
laspect exprimental en ce qui concerne la comparaison de mthodes de contrle non

destructif (CND) par mesures de champs acoustiques, densimtriques et thermiques
sur talons monolithiques et sandwichs ; vis--vis de lobjectif de DAHER-SOCATA
dlargir ses comptences actuelles, il a ainsi t dmontr les perspectives particu-
lirement riches et intressantes de techniques rcentes telles que la tomographie X
ou la thermographie par infrarouge ; les performances de celle-ci savrent en effet
comparables voire meilleures que les mthodes classiques bases notamment sur les
ultrasons, tant sur leur capacit dtecter les dfauts qu en mesurer ltendue ;
laspect numrique vis--vis de lidentification de proprits lastiques de composites

partir de champs cinmatiques ; base sur la mthode de recalage par lments
finis, la procdure mise en place sur des champs simuls permet ce stade dtablir
les modules plans apparents dun matriau orthotrope et de distinguer la prsence
dune zone dgrade de proprits diffrentes.
148 Conclusion gnrale
En perspective ce travail de thse, plusieurs pistes peuvent tre envisages en CND

tout dabord, en terme de technique de contrle, en dehors des mthodes acoustiques et
densimtriques qui sont bien matrises, la thermographie infra-rouge module gagnerait
tre complte par lobtention dune mesure en localisation du dfaut en profondeur
partir du film de phase. Toujours en CND, des tudes de sensibilit doivent tre mises en
place vis--vis des procds de fabrication des composites et des natures des constituants
(armure, fils, densit, poudrage, rsine, ...). Enfin, la nature des dfauts est remettre
en cause, contrairement la spcification AITM ddie principalement aux techniques
ultrasonores, en effet, la thermographie et la tomographie ncessitent des caractristiques
de dfauts diffrentes en terme de proprits thermiques et densimtriques.
En ce qui concerne lidentification partir de mesures de champs cinmatiques, et

plus particulirement lidentification paramtrique multi-zones, des amliorations et enri-
chissements restent tre apports, concernant notamment la caractrisation de la zone
dgrade et surtout lextension de la procdure au cas de champs cinmatiques mesurs
exprimentalement. Plus gnralement enfin, il serait intressant denvisager une passe-
relle entre les deux approches sur le plan industriel, la premire permettant de dtecter
et localiser des dfauts, la seconde donnant des indications sur le niveau de dgradation
associ, dans un but damlioration de laide au diagnostic.
Bibliographie
[1] Adams, R. D. and P. Cawley. 1988. A review of defect types and nondestructive
testing techniques for composites and bonded joints, NDT International , 21(4),
208222.
[2] Adams, R. D. and B. W. Drinkwater. 1997. Nondestructive testing of adhesively

bonded joints, NDT & E International , 30(2), 9398.
[3] Aggelis, D., N.M. Barkoula, T. Matikas, and A. Paipetis. 2012. Acoustic structural
health monitoring of composite materials : damage identification and evaluation in
cross ply laminates using acoustic emission and ultrasonics, Comp. Sci. Tech..
[4] AITM6-0011. 2010. Norme aitm6-0011 : Non destructive inspection of composite

parts.
[5] AITM6-0013. 2008. Norme aitm6-0013 : Evaluation of conventional ultrasonic ins-

pection facilities, equipment and probes.
[6] AITM6-0015. 2008. Norme aitm6-0015 : Thermographic inspection - general.
[7] AITM6-0016. 2008. Norme aitm6-0016 : Evaluation of phased-array ultrasonic ins-

pection facilities, equipment and probes.
[8] AITM6-4002. 2008. Norme aitm6-4002 : Ultrasonic through transmission - inspec-

tion of fibre composite.
[9] AITM6-4005. 2010. Norme aitm6-4005 : Ultrasonic pulse-echo inspection of carbon

fibre plastics.
[10] AITM6-4010. 2008. Norme aitm6-4010 : Inspection of fibre composites by ultrasonic

phased array technique.
[11] AITM6-4014. 2009. Norme aitm6-4014 : Thickness measurement of fiber composites

by ultrasonic pulse-echo technique.
150 BIBLIOGRAPHIE
[12] AITM6-7002. 2007. Norme aitm6-7002 : Radiographic inspection-general.
[13] AITM6-7006. 2008. Norme aitm6-7006 : Computed tomography with x-rays.
[14] Al-Maghribi, A. 2008. Comportement des matriaux composites fibres courtes :

application limpact basse vitesse, Ph.D. thesis, Universit Paul-Sabatier.
[15] Amiot, F., F. Hild, J. N. Pri, and S. Roux. 2009. Mthode de lcart lquilibre,
Tech. rep., Ecole thmatique CNRS, Clermont-Ferrand, France.
[16] Andrieux, S., A. Ben Abda, and H. D. Bui. 1997. Sur lidentification de fissures
planes via le concept de lcart la rciprocit en lasticit, C.R. Acad. Sci. Paris,
324, 14311438.
[17] Andrieux, S. and A. Constantinescu. 2009. Ecart la rciprocit, Tech. rep., Ecole
thmatique CNRS, Clermont-Ferrand, France.
[18] Avril, C., D. Perreux, F. Thiebaud, and B. Reck. 2009. Optimisation du compor-
tement mcanique de structures composites soumises des hyper-acclrations, in
JNC16 , Toulouse, France.
[19] Avril, S. 2002. Application des mthodes de mesure de champs la caractrisation

mcanique de poutre en bton arm rpares par matriaux composites, Ph.D. thesis,
Ecole nationale suprieure des Mines de Saint-Etienne et Universit Jean Monnet.
[20] Avril, S., M. Bonnet, A. S. Bretelle, M. Grdiac, F. Hild, P. Ienny, F. Latourte,

D. Lemosse, S. Pagano, E. Pagnacco, and F. Pierron. 2008. Overview of identification
methods of mechanical parameters based on full-field measurements, Experimental
Mechanics, 48, 381402.
[21] Avril, S. and F. Pierron. 2007. General framework for identification of constitutive
parameters from full field measurements in linear elasticity, International Journal
of Solids and Structures, 44, 49785002.
[22] Babout, L. 2002. Etude par tomographie X et modlisation de lendommagement de

matriaux mtalliques modles, Ph.D. thesis, Institut National des Sciences Appli-
ques de Lyon.
[23] Babout, L., E. Maire, J.Y. Buffire, and R. Fougres. 2001. Characterization by x
ray computed tomography of decohesion, porosity growth and coalescence in model
metal matrix composites, Acta Mater., 49, 20552063.
BIBLIOGRAPHIE 151
[24] Badel, P. 2008. Analyse msoscopique du comportement mcanique des renforts tis-
ss de composites utilisant la tomographie aux rayons X , Ph.D. thesis, Institut Na-
tional des Sciences Appliques de Lyon.
[25] Bai, W. and B. S. Wong. 2001. Evaluation of defects in composite plates under
convective environments using lock-in thermography, Meas. Sci. Technol., 12, 142
150.
[26] Bai, W. and B. S. Wong. 2001. Photothermal models for lock-in thermographic
evaluation of plates with finite thickness under convection conditions, Journal of
applied physics, 89(6), 32753282.
[27] Bai, W. and B.S. Wong. 2000. Non-destructive evaluation of aircraft structure using
lock-in thermography, in SPIE - The International Society for Optical Engineering,
vol. 3994, 3746.
[28] Bar-Cohen, Y., A. K. Mal, C. J. Hellier, W. Plumstead, K. Fowler, R. Grills, G. An-

derews, M. C. Tsao, J. J. Snyder, J. F. Cook, D. A. Aldrich, and R. W. Pepper.
1996. Non destructive evaluation and quality control , ASM International, vol. 17,
chap. Ultrasonic Inspection, 231277.
[29] Batsale, J. C., J. L. Battaglia, and D. Mourand. 2001. Traitement dimages infra-
rouges fortement bruites associes au transfert thermique transitoire dans les mi-
lieux htrognes : application au contrle non destructif de faible cout, in XVme
Congrs Franais de Mcanique, Nancy, France.
[30] Batsale, J.C., D. Maillet, and A. Degiovanni. 1994. Extension de la mthode des
quadriples thermiques laide de transformations intgrales - applications au d-
faut plan bidimensionnel, International Journal of Heat Transfer , 37, 111127.
[31] Bay, B., T. Smith, D. Fyhrie, and M. Saad. 1999. Digital volume correlation :
Three-dimensional strain mapping using x-ray tomography., Experimental Mecha-
nics, 39(3), 217226.
[32] Bayraktara, E., S. D. Antolovich, and C. Bathias. 2008. New developments in non-
destructive controls of the composite materials and applications in manufacturing
engineering, Journal of materials processing technology, 206, 3044.
[33] Bendada, A., J. C. Batsale, A. Dediovanni, and D. Mallet. 1994. Estimation de la

rpartition spatiale dune rsistance thermique dinterface partir de la mesure de
temprature de surface, in Journe dtude SFT-GUT .
152 BIBLIOGRAPHIE
[34] Bendada, A., D. Maillet, J. C. Batsale, and A. Degiovanni. 1998. Reconstruction of

a non-uniform interface thermal resistance by inverse conduction, Inverse Problems
in Engineering, 6, 79123.
[35] Benmedakhene, S., M. Kenane, and M. Benzeggagh. 1999. Initiation and growth of
delamination in glass/epoxy composites subjected to static and dynamic loading by
acoustic emission monitoring, Comp. Sci. Tech, 59, 201208.
[36] Bernasconi, A., F. Cosmi, and D. Dreossi. 2008. Local anisotropy analysis of injec-
tion moulded fibre reinforced polymer composites, Composites Science and Techno-
logy, 68, 25742581.
[37] Bertin, M., F. Touchard, and M. C. Lafarie-Frenot. 2011. Stratifi carbone epoxy
sous sollicitation de fatigue thermomcanique simulant le remplissage rapide de
rservoirs hydrogne : analyse exprimentale de linfluence de la temprature, in
JNC 17 , Poitiers, France.
[38] Bilek, A., J. C. Dupr, F. Bremand, and A. Ouibrahim. 2004. Studies of contact
problems by 3d photoelasticity, comparison with finite element analysis, in Interna-
tional Conference on Experimental Mechanics - ICEM12 .
[39] Bodnar, J. L. 2010. Contrle non destructif par techniques infrarouges : principes et
exemples dapplications, in Colloque Contrles et Mesures Optiques pour lIndustrie
CMOI , Toulouse-Labge, France.
[40] Bonnet, M. 2009. Prsentation gnrale de la problmatique et des mthodes

didentification disponibles, Tech. rep., Ecole thmatique CNRS, Clermont-Ferrand,
France.
[41] Bonnet, M., H. D. Bui, and A. Constantinescu. 2003. Principes variationnels et

exploitation de mesures de champs en lasticit, Mecanique & Industries, 4(6), 687
697.
[42] Bonnet, M. and S. Pagano. 2009. Erreur en relation de comportement, Tech. rep.,
Ecole thmatique CNRS, Clermont-Ferrand, France.
[43] Bornert, M., P. Doumalin, and E. Maire. 2004. Mesure tridimensionnelle des champs
de dformation par corrlation dimages microtomographiques, in Colloque Photo-
mcanique 2004 .
BIBLIOGRAPHIE 153
[44] Boro-Djordjevic, B. 2009. Nondestructive test technology for the composites : Appli-
cation of contemporary non-destructive testing in engineering, in The 10th Interna-
tional Conference of the Slovenian Society for Non-Destructive Testing, Ljubljana,
Slovenia.
[45] Bossi, R., K. Friddell, and A. Lowrey. 1990. Non destructive testing of fiber reinforced
plastic composites - Ch. Computed tomography, London : Elsevier Applied Science
Publishers.
[46] Boyer, Y. 1993. Caractrisation de la porosit de composite carbone-rsine minces

par ltude de la rflexion des ondes ultrasonores, Ph.D. thesis, Universit de tech-
nologie de Compigne.
[47] Bremond, P. 2001. Lock-in thermography : a tool to mesure stress and to detect
defects in aircraft industry, in AIAC , Canberra Australia.
[48] Breque, C., J. C. Dupre, and F. Bremand. 2004. Calibration of a system of projec-
tion moir for relief measuring : biomechanical applications, Optics and Lasers in
Engineering, 41(2), 241 260.
[49] Bui, H. D. 1993. Introduction aux problmes inverses en mcanique des matriaux .
[50] Carmona, V. 2009. Etude de lendommagement de matriaux composites par to-

mographie X et mission acoustique, Ph.D. thesis, Institut National des Sciences
Appliques de Lyon.
[51] Castaings, M. 1993. Propagation ultrasonore dans les milieux stratifis plans consti-
tus de matriaux absorbants et orthotropes, Ph.D. thesis, Universit de Bordeaux
I.
[52] Cawley, P. 2006. Inspection of composites : Current status and challenges, in

ECNDT .
[53] Chaix, J. F., V. Garnier, and G. Corneloup. 2006. Ultrasonic wave propagation
in heterogeneous solid media : Theoretical analysis and experimental validation,
Ultrasonics, 44, 200210.
[54] Chen, T. Y. 2000. Digital photoelasticity, in Photomechanics, 197232.
[55] Cherfaoui, M. 2006. Essais non destructifs chapitre contrle par ultrason, Tech. rep.,
Techniques de lingnieur.
154 BIBLIOGRAPHIE
[56] Choi, M., K. Kang, J. Park, W. Kim, and K. Kim. 2008. Quantitative determination
of a subsurface defect of reference specimen by lock-in infrared thermography, NDT
& E International , 41, 119124.
[57] Chrysochoos, A. 2002. La thermographie infrarouge un outil en puissance pour

tudier le comportement des matriaux, Mcanique & Industries, 3, 314.
[58] Claire, D., F. Hild, and S. Roux. 2004. A finite element formulation to identify
damage fields : The equilibrium gap, Int. J. Num. Meth. Engng., 61(2), 189208.
[59] Claire, D., F. Hild, and S. Roux. 2007. Identification of a damage law by using full
field displacement measurements, Int. J. Damage Mech., 16, 179197.
[60] Clerjaud, L. 2010. Mthode dhtrodynage pour la caractrisation de proprits ther-

mophysiques par thermographie infrarouge dans une large gamme spatiale et tempo-
relle, Ph.D. thesis, Universit de Bordeaux I.
[61] Clorennec, D. 2002. Gnration et dtection optiques dondes guides sur une pice
cylindrique : Application au contrle non destructif sans contact, Ph.D. thesis, Uni-
versit Paris 7.
[62] Colin, N. 1990. Etude comparative de mthodes danalyse de signaux dchographie

ultrasonore dun matriau composite, Ph.D. thesis, Universit de technologie de
Compigne.
[63] Cooreman, S. 2008. Identication of the plastic material behaviour through full-field
displacement measurements and inverse methods, Ph.D. thesis, Universit de Brus-
sel.
[64] Cornille, N. 2005. Accurate 3D Shape and Displacement Measurement using a Scan-
ning Electron Microscope, Ph.D. thesis, Ecole des Mines dAlbi-Carmaux et Uni-
versit de Caroline du Sud.
[65] Crouzeix, L. 2008. Identification de champs de proprits mcaniques de structures

composites partir de mesures de champs de dplacements, Ph.D. thesis, Universit
de Paul Sabatier Toulouse III.
[66] David, D., J. C. Batsale, and J. Gounot. 1998. New sequential method to process
noisy temperature reponse from flash experiment measured by infrared camera,
Review of Scientific Instrument, 69(3), 14371441.
BIBLIOGRAPHIE 155
[67] Decultot, N. 2009. Formage incrmental de tle daluminium : tude du procd

laide de la mesure de champs et identification de modles de comportement, Ph.D.
thesis, Universit de Paul Sabatier Toulouse III.
[68] Degiovanni, A. 1977. Diffusivit et mthode flash, Rev. Gn. Thermique, 185, 420
441.
[69] Degiovanni, A. and M. Laurent. 1986. Une nouvelle technique didentification de la

diffusivit thermique par la mthode flash, Rev. Gn. Phys. Appl., 21, 229237.
[70] Doumalin, P., J. Drouot, and J. C. Dupr. 2004. Mesure du champ de dformation
par corrlation dimages tridimensionnelles : application aux matriaux transpa-
rents, in Colloque Photomcanique.
[71] Dumont-Fillon, J. 1996. Contrle non destructif chapitre ultrasons, Tech. rep., Tech-
niques de lingnieur.
[72] Duplessis-Kergomard, Y. 2007. Etude exprimentale et numrique de la fissuration

intralaminaire et de linitiation du dlaminage de bord libre dans les structures stra-
tifies pseudo-tisses, Ph.D. thesis, Ecole des Mines de Paris.
[73] Dupr, J. and A. Lagarde. 1997. Photoelastic analysis of a three-dimensional speci-

men by optical slicing and digital image processing, Experimental Mechanics, 37(4),
393397.
[74] Dzenis, Y. and J. Qian. 2001. Analysis of microdamage evolution histories in com-
posites, Int. J. Solids Struct., 38, 18311854.
[75] Eyma, F., D. Marlot, L. E. Denaud, and J. N. Felices. 2009. Utilisation dune m-
thode ultrasonore sans contacts pour la recherche de dfauts dans le contreplaqu
de peuplier : manque de matires/dfauts de collage, in COFREND2009 , Mont de
Marsant, France.
[76] Fazzini, M. 2009. Dveloppement de mthodes dintgration des mesures de champs,

Ph.D. thesis, Universit de Paul Sabatier Toulouse III.
[77] Fazzini, M., O. Dalverny, and S. Mistou. 2011. Identification of materials properties
using displacement field measurement, Key Engineering Materials, 482, 5765.
[78] Fidahoussen, A., P. Calmon, and M. Lambert. 2008. Dveloppement de techniques

de reconstitution ultrasonore multi-lments pour la localisation et la caractrisation
de dfauts dans des pices de gomtrie complexe, in COFREND.
156 BIBLIOGRAPHIE
[79] Framezelle, G., F. Lepiller, J. Hatsch, and S. Remond. 2011. Contrle non destructif
de structures composites pour laronautique, in COFREND.
[80] Franon, M. 1986. Loptique moderne et ses dveloppements depuis lapparition du

laser , Paris : CNRS.
[81] Freud, N. 2003. Modlisation et simulation de systmes dimagerie par rayons X ou

gamma, Ph.D. thesis, Institut Nationnal des Sciences Appliques de Lyon.
[82] Garcia, D. 2001. Mesure de formes et de champs de dplacements tridimensionels

par stro-corrlation dimages, Ph.D. thesis, Institut National Polytechnique de
Toulouse.
[83] Garnier, C. 2011. Etude du comportement dynamique des structures composites ra-
lises par LRI : application limpact et la fatigue, Ph.D. thesis, Institut National
Polytechnique de Toulouse.
[84] Geers, M. G. D., I. R. De Borsts, and W. A. M. Brekelmans. 1996. Computing strain

fields from discrete displacement fields in 2d-solids, International Journal of Solids
and Structures, 33(29), 42934307.
[85] Geier, M. 1989. Manuel qualit des composites.
[86] Germaneau, A. 2007. Dveloppement de techniques de mesure dans le volume : pho-

tolasticimtrie 3D par dcoupage optique et corrlation volumique par tomographie
optique et rayon X. Application ltude des effets mcaniques 3D dans les struc-
tures et les biomatriaux , Ph.D. thesis, Universit de Poitiers.
[87] Geymonat, G., F. Hild, and S. Pagano. 2002. Identication of elastic parameters by
displacement field measurement, Comptes Rendus Mecanique, 330(6), 403408.
[88] Giorleo, G. and C. Meola. 2002. Comparison between pulsed and modulated ther-
mography in glas epoxy laminates, NDT & E International , 35, 287292.
[89] Gleiter, A., G. Riegert, T. Zweschper, and G. Busse. 2006. Ultrasound lockin ther-
mography for advanced depth resolved defect selective imaging, in ECNDT2006 .
[90] Goidescu, C. 2011. Caractrisation et modlisation de lendommagement par micro-

fissuration des composites stratifis - Apports des mesures de champs et de lhomo-
gnisation, Ph.D. thesis, Institut National Polytechnique de Toulouse.
BIBLIOGRAPHIE 157
[91] Grediac, M. 1989. Principe des travaux virtuels et identification, Comptes rendus
de lacadmie des sciences, 309, 15.
[92] Grediac, M. 1996. On the direct determination of invariant parameters govering

anisotropic plate bending problems, International Journal of Solids and Structures,
33(27), 39693982.
[93] Grediac, M. 2004. The full-field measurement methods in composite material cha-
racterization : interest and limitations, Composites, 35, 751761.
[94] Grediac, M. 2010. Identication et mesures de champs, in Colloque National MECA-

MAT , Aussois, France.
[95] Grediac, M., S. Avril, and E. Toussaint. 2009. La mthode des champs virtuels,
Tech. rep., Ecole thmatique CNRS, Clermont-Ferrand, France.
[96] Grediac, M. and P. A. Paris. 1996. Direct identification of elastic constants of ani-
sotropic plates by modal analysis : theoretical and numerical aspects, Journal of
Sound and Vibration, 195(3), 401415.
[97] Grediac, M. and F. Pierron. 2006. Applying the virtual fields method to the identifi-
cation of elasto-plastic constitutive parameters, International Journal of Plasticity,
22(4), 602627.
[98] Grediac, M., E. Toussaint, and F. Pierron. 2002. Identification des proprits m-
caniques de matriaux avec la mthode des champs virtuels, une alternative au
recalage par lments finis, C. R. Mecanique, 330, 107112.
[99] Grediac, M., E. Toussaint, and F. Pierron. 2002. Special virtual fields for the direct
determination of material parameters with the virtual fields method, International
Journal of Solids and Structures, 39, 26912705.
[100] Grediac, M. and A. Vautrin. 1990. A new method for determination of bending
rigidities of thin anisotropic plates, Journal of Applied Mechanics, 57, 964968.
[101] Gros, X. E. 1995. An eddy current approach to the detection of damage caused by
low-energy impacts on carbon fibre reinforced materials, Materials & Design, 16(3),
167173.
[102] Guibert, S. 2007. La thermographie infrarouge dtection synchrone applique aux

matriaux composites, Masters thesis, Universit de LAVAL, QUEBEC.
158 BIBLIOGRAPHIE
[103] Guo, B. 2007. Identification simultane des paramtres de rigidit et damortis-

sement de plaques isotropes minces en vibration par la mthode des champs vir-
tuels, Ph.D. thesis, Ecole Nationale Suprieure dArts et Mtiers en Chlon-en-
Champagne.
[104] Hanke, R., T. Fuchs, and N. Uhlmann. 2008. X ray based methods for non des-
tructive testing and material characterization, Nuclear Instruments and Methods in
Physics Research, 591, 1418.
[105] Herve, C., F. Berthelot, and H. Walaszek. 2008. Application des mthodologies
dexamen non destructives pour la fiabilit des structures en matriaux composites,
in COFREND2008 .
[106] Horny, N. and B. Lannoy. 2003. Lock-in thermography with a focal plane array,
Meas. Sci. Technol., 14, 439443.
[107] Hsu, D. K. 2008. Nondestructive inspection of composite structures : methods and

practice, in 17th World Conference on Nondestructive Testing, Shanghai, China.
[108] Hung, Y. Y. 1996. Shearography for non-destructive evaluation of composite struc-

tures, Optics and Lasers in Engineering, 24, 161182.
[109] Hung, Y. Y., Y. S. Chen, S. P. Ng, L. Liu, Y. H. Huang, B. L. Luk, R. W. L.

Ip, C. M. L. Wu, and P. S. Chung. 2009. Review and comparison of shearogra-
phy and active thermography for nondestructive evaluation, Materials Science and
Engineering, 64, 73112.
[110] Ibarra-Castanedo, C. 2005. Quantitative subsurface defect evaluation by pulsed phase

thermography : depth retrieval with the phase, Ph.D. thesis, Universit de LAVAL,
QUEBEC.
[111] Ibarra-Castanedo, C., M. Genest, S. Guibert, J. M. Piau, X. P. V. Maldague, and

A. Bendada. 2007. Inspection of aerospace materials by pulsed thermography, lock-
in thermography and vibrothermography : A comparative study, in SPIE - The
International Society for Optical Engineering, vol. 6541.
[112] Ibarra-Castanedo, C. and X. Maldague. 2004. Defect depth retrieval from pulsed
phase thermographic data on plexiglas and aluminium samples, in SPIE -The in-
ternational Society for optical Engineering, vol. 5405, 348356.
BIBLIOGRAPHIE 159
[113] Ibarra-Castanedo, C. and X. Maldague. 2005. Pulsed phase thermography inversion

procedure using normalized parameters to account for defect size variations, in SPIE
-The international Society for optical Engineering, vol. 5782, 334341.
[114] Jessop, H. T. and F. C. Harris. 1952. Photolasticit, Paris.
[115] Junyan, L., W. Yang, and D. Jingmin. 2010. Research on thermal wave processing of
lock-in thermography based on analyzing image sequences for ndt, Infrared Physics
& Technology, 53, 348357.
[116] Kazys, R., A. Demcenko, E. Z. ukauskas, and L. Mazeika. 2006. Air coupled ultra-
sonic investigation of multi-layered composite materials, Ultrasonics, 44, 819822.
[117] Kim, J. H., F. Pierron, M. R. Wisnom, and K. Syed-Muhamad. 2007. Identification

of the local stiffness reduction of a damaged composite plate using the virtual fields
method, Composites Part A : Applied Science and Manufacturing, 38, 20652075.
[118] Kleinermann, J. P. 2000. Identification paramtrique et optimisation procds mise

forme par problmes inverses, Ph.D. thesis, Universit de Lige.
[119] Lacroix, M. 1986. Essais non destructifs chapitre ultrasons, Tech. rep., Techniques
de lingnieur.
[120] Laermann, K.H. 2000. Optical methods in experimental solid mechanics, New York.
[121] Lambert, A., J. Rivenez, G. Wach, and M. Cherfaoui. 1994. Les contrles non
destructifs gnralits, Tech. rep., Cahiers Formation CETIM.
[122] Larbi, W. B. 2010. Comparaison exprimentale de la thermographie module et de

la thermographie pulse pour lvaluation non destructive des matriaux employs
en arospatial , Masters thesis, Universit de LAVAL.
[123] Latourte, F. 2007. Identification des paramtres dune loi lastoplastique de prager
et calcul de champs de contrainte dans des matriaux htrognes, Ph.D. thesis,
Universit de Montpellier II.
[124] Latourte, F., A. Chrysochoos, S. Pagano, and B. Watrisse. 2008. Elastoplastic be-
havior identification for heterogeneous loadings and materials, Experimental Mecha-
nics, 48(4), 435449.
160 BIBLIOGRAPHIE
[125] Le-Bras, A. 2004. Exploration des potentialits du systme eos pour la caracterisation
mcanique de structures osseuses : application lextrmit suprieure du femur ,
Ph.D. thesis, Ecole Nationale Suprieure Arts et Mtiers.
[126] Le-Magorou, L., F. Bos, and F. Rouger. 2002. Identication of constitutive laws
for wood-based panels by means of an inverse methode, Composites Science and
Technology, 62(4), 591596.
[127] Leclerc, H., J. N. Perie, F. Hild, and S. Roux. 2009. Corrlation dimages numriques
intgre pour identification de proprits mcaniques, in 19 me Congrs Franais
de Mcanique, Marseille, France.
[128] Lecompte, D. 2007. Elastic and Elasto-Plastic Material Parameter Identification

by Inverse Modeling of Static Tests using Digital Image Correlation, Ph.D. thesis,
Universit de Brussel.
[129] Lecompte, D., A. Smits, H. Sol, J. Vantomme, and D. Van Hemelrijck. 2007.
Mixed numerical-experimental technique for orthotropic parameter identication
using biaxial tensile tests on cruciform specimens, International Journal of Solids
and Structures, 44, 16431656.
[130] Lecomte, G. 2005. Analyse dimages radioscopiques et fusion dinformations multi-

modales pour lamlioration du contrle de pices de fonderie, Ph.D. thesis, Institut
National des Sciences Appliques de Lyon.
[131] Ledru, Y. 2009. Etude de la porosit dans les matriaux composites stratifis aro-
nautiques, Ph.D. thesis, Universit de Toulouse.
[132] Leendertz, J. A. and J. N. Butters. 1973. An image-shearing speckle-pattern inter-

ferometer for measuring bending moments, Journal of Physics E : Scientific Instru-
ments, 6(11), 11071110.
[133] Lefebvre, J. P., P. Lasaygue, C. Potel, J. F. De Belleval, and P. Gatignol. 2004.

Acoustique ultrasonore et ses applications, Tech. rep.
[134] Lenoir, N. 2006. Comportement mcanique et rupture dans les roches argileuse tu-
dis par micro-tomographie rayon X , Ph.D. thesis, Universit Joseph Fourier Gre-
noble I.
BIBLIOGRAPHIE 161
[135] Maillard, S., J. Cadith, H. Walaszek, A. Dillenz, and J. L. Bodnar. 2008. La ther-
mographie infrarouge stimule, une nouvelle technique de contrle sur les lignes de
production, in Cofrend .
[136] Maire, E., J. Y. Buffire, L. Salvo, J. J. Blandin, W. Ludwig, and J. M. Letang.

2001. On the application of x-ray tomography in the field of materials science, Adv.
Engng. Sci., 3, 539546.
[137] Maire, E., V. Carmona, J. Courbon, and W. Ludwig. 2007. Fast x-ray tomography
and acoustic emission study of damage in metals during continuous tensile tests,
Acta Mater., 55, 68066815.
[138] Meola, C., G. M. Carlomagno, A. Squillace, and G. Giorleo. 2002. Non-destructive

control of industrial materials by means of lock-in thermography, Meas. Sci. Tech-
nol., 13, 15831590.
[139] Meola, C., G. M. Carlomagno, A. Squillace, and G. Giorleo. 2004. The use of infrared
thermography for nondestructive evaluation of joints, Infrared Physics & Techno-
logy, 46, 9399.
[140] Meola, C., G. M. Carlomagno, A. Squillace, and A. Vitiello. 2006. Non-destructive

evaluation of aerospace materials with lock-in thermography, Engineering Failure
Analysis, 13, 380388.
[141] Meraghni, F., H. Chalal, F. Pierron, and M. Grdiac. 2002. Identification directe du
comportement lastique endommageable de matriaux composites par la mthode
des champs virtuels, Journal de Physique IV , 12, 291300.
[142] Merzouki, T., F. Meraghni, H. Chalal, and T. B. Zineb. 2009. Comparaison de

configurations dessais htrognes en vue damliorer lidentification du compor-
tement lastique orthotrope par erreur en relation de comportement, Mcanique &
Industries, 10, 351364.
[143] Mistou, S. 2000. Dveloppement dun lment fini de coque composite double cour-
bure pour limplmentation dune application oriente objet, Ph.D. thesis, Universit
de Bordeaux I.
[144] Mizutani, Y., K. Nagashima, M. Takemoto, and K. Ono. 2000. Fracture mecha-
nism characterization of cross-ply carbon-fiber composites using acoustic emission
analysis, NDT & E Int., 33, 101110.
162 BIBLIOGRAPHIE
[145] Moffat, A. J., P. Wright, J. Y. Buffire, I. Sinclair, and S. M. Spearing. 2008. Micro-
mechanisms of damage in 0 splits in a [90/0]s composite material using synchrotron
radiation computed tomography, Scripta Materialia, 59, 10431046.
[146] Molimard, J., R. Le Riche, A. Vautrin, and J. R. Lee. 2004. Identification of orthotro-
pic plate stiffness using open hole tensile test, in Xth SEM congress on experimental
mechanics, Costa Mesa (Californie, USA).
[147] Molimard, J., R. Le Riche, A. Vautrin, and J. R. Lee. 2005. Identication of the
four orthotropic plate stinesses using a single open-hole tensile test, Experimental
Mechanics, 45(5), 404411.
[148] Montanini, R. 2010. Quantitative determination of subsurface defects in a reference

specimen made of plexiglas by means of lock-in and pulse phase infrared thermo-
graphy, Infrared Physics & Technology, 53, 363371.
[149] Mook, G., R. Lange, and O. Koeser. 2001. Non-destructive characterisation of

carbon-fibre-reinforced plastics by means of eddy-currents, Composites Science and
Technology, 61, 865873.
[150] Morimoto, Y., Y. Morimoto Jr., and T. Hayashi. 1994. Separation of isochromatics
and isoclinics using fourier transform, Experimental Techniques, 18, 1317.
[151] Mouritz, A. P., C. Townsend, and M. Z. Shah Khan. 2000. Non destructive detection
of fatigue damage in thick composites by pulse-echo ultrasonics, Composites Science
and Technology, 60, 2332.
[152] Muzet, V., P. Blain, and Y. Guillard. 2008. Application de la shearographie la

dtection de fissures sur ouvrages dart, Tech. rep., LRPC de Strasbourg.
[153] Orteu, J. J. 2009. 3D computer vision in experimental mechanics, Optics and Lasers
in Engineering, 47(3-4), 282291.
[154] Pagnacco, E. and A. S. Caro-Bretelle. 2009. Identification par recalage de modles

elments finis, Tech. rep., Ecole thmatique CNRS, Clermont-Ferrand, France.
[155] Pajani, D. 2001. Thermographie - principes et mesure, Tech. rep., Techniques de

lingnieur.
[156] Pajani, D. and L. Audaire. 2001. Thermographie technologies et application, Tech.

rep., Techniques de lingnieur.
BIBLIOGRAPHIE 163
[157] Pannier, Y. 2006. Identification de paramtres lastoplastiques par des essais stati-
quement indtermins : mise en oeuvre exprimentale et validation de la mthode
des champs virtuels, Ph.D. thesis, Ecole Nationale Suprieure dArts et Mtiers de
Centre de Chlons-en-Champagne.
[158] Parker, W. J., R. J. Jenkins, C. P. Butler, and G. L. Abbott. 1961. Flash method
of determining thermal diffusivity, heat capacity and thermal conductivity, Journal
of applied physics, 32(9), 16791684.
[159] Pastor, M. L. 2007. Contribution ltude en fatigue de structures en aluminium

renforces par patchs composites, Ph.D. thesis, Universit Blaise Pascal - Clermont
II.
[160] Patterson, E. A. and Z. F. Wang. 1991. Towards full field automated photoelastic
analysis of complex components, Strain, 27, 4953.
[161] Perdijon, J. 1993. Le contrle non destructif par ultrasons.
[162] Perdijon, J. 2002. Contrle des matriaux .
[163] Perie, J. N., H. Leclerc, S. Roux, and F. Hild. 2009. Digital image correlation and
biaxial test on composite material for anisotropic damage law identification, Inter-
national journal of solids and structures, 46, 23882396.
[164] Perie, J. N., H. Leclerc, S. Roux, and F. Hild. 2009. Identification de lois dendom-
magement anisotrope partir de la corrlation dimages numriques, in Comptes
rendus des JNC 16 , Toulouse, France.
[165] Phani, S. K. M. 2007. Characterization of advanced composites : a nondestructive

approach, Ph.D. thesis, Louisiana State University and Agricultural and Mechanical
College.
[166] Pierron, F., S. Zhavoronok, and M. Grediac. 2000. Identification of the through-
thickness properties of thick laminated tubes using the virtual fields method, Inter-
national Journal of Solids and Structures, 37(32), 44374453.
[167] Post, D., B. Han, and Pr. G. Ifju. 2000. Moir methods for engineering and science
- moir interferometry and shadow moir, in Photomechanics, 151196.
[168] Potet, P. 2004. Connaissez-vous la vibrothermographie ?, Mesures, 768, 6466.

164 BIBLIOGRAPHIE
[169] Pradere, C., L. Clerjaud, J. C. Batsale, and S. Dilhaire. 2011. High speed hetero-
dyne infrared thermography applied to thermal diffusivity identification, Scientific
Instruments, 82(054901), 16.
[170] Prodjinonto, V. 2011. Contribution lconomie dnergie dans le btiment : mesure

de capacit de stockage dynamique dune paroi , Ph.D. thesis, Universit de Bordeaux
I.
[171] Ramond, L., J. C. Batsale, and C. Gobbe. 1999. Mthodologie de mesure de coef-
ficient dchange en milieu poreux - application dun modle deux tempratures,
Int. J. of Thermal Sciences, 38, 250257.
[172] Ramond, L., J. C. Batsale, C. Gobbe, and M. Quintard. 1998. Measurement of

a thermal exchange coefficient in porous media by infrared thermography : two
temperature model application, in 11 th Int. Heat Transfer Conference, Kyong Ju
Korea.
[173] Rastogi, P. K. 2000. Principles of holographic interferometry and speckle metrology,

in Photomechanics, 103150.
[174] Reifsnider, K. L. and R. S. Williams. 1974. Determination of fatigue related heat

emission in composite materials, Exp. Mech., 14, 479485.
[175] Rikards, R., A. Chate, and G. Gailis. 2001. Identification of elastic properties of la-
minates based on experiment design, International Journal of solids and structures,
38, 50975115.
[176] Robert, E. and Jr. Green. 2004. Non contact ultrasonic techniques, Ultrasonics, 42,
916.
[177] Rodriguez, F. L., V. De Paulo Nicolau, H. Oshiro, and D. Tancredi. 2011. Non-
destructive evaluation of composites materials by pulsed-phase thermography :
depth inversion, in COBEM 2011 - 21st Brazilian Congress of Mechanical Engi-
neering.
[178] Roux, S. and F. Hild. 2008. Digital image mechanical identification, Experimental
Mechanics, 48, 495508.
[179] Ruzek, R., R. Lohonka, and J. Jironc. 2006. Ultrasonic C scan and shearography
NDI techniques evaluation of impact defects identification, NDT & E International ,
39, 132142.
BIBLIOGRAPHIE 165
[180] Sakagami, T. and S. Kubo. 2002. Application of pulse heating thermography and
lock-in thermography to quantitative nondestructive evaluations, Infrared Physics
& Technology, 43, 211218.
[181] Schilling, P. J., B. R. Karedla, A. K. Tatiparthi, M. A. Verges, and P. D. Herring-

ton. 2005. X-ray computed microtomography of internal damage in fiber reinforced
polymer matrix composites, Composites Science and Technology, 65, 20712078.
[182] Simulia. 2002. Simulia : http ://www.simulia.com.
[183] Smigielski, P. 1984. Holographie industrielle, Toulouse.
[184] Smigielski, P. 1998. Holographie optique, principes, Tech. Rep. AF 3340, Techniques
de lIngnieur, trait Sciences fondamentales.
[185] Smigielski, P. 1998. Interfromtrie holographique, principes, Tech. Rep. AF 3345,

Techniques de lIngnieur, trait Sciences fondamentales.
[186] Smigielski, P. 2001. Interfromtrie de speckle, Tech. Rep. R 6331, Techniques de

lIngnieur, trait Mesures et contrle.
[187] Sreelekshmi, C., V.R. Ravindran, and V.P. Mahadevan Pillai. 2006. Theoretical
modeling of digital radiograph of composite rocket nozzle for the analysis of defect,
in National Seminar on Non-Destructive Evaluation.
[188] Steinberger, R., T. I. Valadas Leitao, E. Ladstatter, G. Pinter, W. Billinger, and

R.W. Lang. 2006. Infrared thermographic techniques for non-destructive damage
characterization of carbon fibre reinforced polymers during tensile fatigue testing,
International Journal of Fatigue, 28, 13401347.
[189] Surrel, Y. 2000. Fringe analysis, in Photomechanics, 55102.
[190] Sutton, M. A., S. R. McNeill, J. D. Helm, and Y. J. Chao. 2000. Advances in two-
dimensional and three-dimensional computer vision, in Photomechanics, 323372.
[191] Thomas, C. 2011. Etude des mcanismes dendommagement des composites fibres de
carbone / matrice polyamide : application la ralisation de rservoirs de stockage
de gaz sous haute pression de type IV , Ph.D. thesis, Ecole nationale suprieure des
mines de Paris.
166 BIBLIOGRAPHIE
[192] Tsao, C. C. and H. Hocheng. 2005. Computerized tomography and c-scan for mea-
suring delamination in the drilling of composite materials using various drills, In-
ternational Journal of Machine Tools & Manufacture, 45, 12821287.
[193] Vautrin, A., G. Silva, J. Molimard, and R. Le Riche. 2009. Localisation et iden-
tification de lhtrognit de comportement mcanique de matriaux, in 19 me
Congrs Franais de Mcanique, Marseille, France.
[194] Voillaume, H. and B. Campagne. 2008. Les ultrasons laser appliqus aux matriaux
composites aronautiques, in Cofrend .
[195] Wevers, M. 1997. Listening to the sound of materials : acoustic emission for the
analysis of material behaviour, NDT & E Int., 30, 99106.
[196] Wong, B. S. and C. F. M. Ron. 1998. Non destructive testing of fibre reinforced
composites and honeycomb structures, Tech. rep., Nanyang Technological Univer-
sity, Singapore.
[197] Wu, D. and G. Busse. 1998. Lock-in thermography for nondestructive evaluation of
materials, Revue Gnrale de Thermique, 37(8), 693703.
[198] Yashiro, S., J. Takatsubo, and N. Toyama. 2007. An ndt technique for composite
structures using visualized lamb-wave propagation, Composites Science and Tech-
nology, 67, 32023208.
[199] Zani, M. L. 2001. CND par ultrason : le sans contact est prt mais..., Mesures, 734,
6366.
[200] Zenina, A. 1998. Analyse par photolasticimtrie dune pice tridimensionnelle :

modlisations et exprimentations, Ph.D. thesis, Universit de Poitiers.
[201] Zenina, A., J.C. Dupr, and A. Lagarde. 1999. Separation of isochromatic and iso-
clinic patterns of a slice optically isolated in a 3d photoelastic medium, European
Journal of Mechanics - A/Solids, 18(4), 633 640.
[202] Zweschper, Th., A. Dillenz, , and G. Busse. 2000. Ultrasound lockin thermography
- a ndt method for the inspection of aerospace structures, in QIRT2000 .
Doctorat de lUniversit de Toulouse
Dlivr par : Institut National Polytechnique de Toulouse
Ecole doctorale : MEGeP
Spcialit : Gnie mcanique, mcanique des matriaux
Date de soutenance : Jeudi 04 Octobre 2012
Auteur : Elodie PERONNET
Directeurs de thse : Moussa KARAMA / Sbastien MISTOU
Titre : Identification exprimentale du comportement dun fuselage composite : dtection de dfauts

par mesures de champs
Rsum :
Ces travaux concernent lidentification du comportement dde structures composites. Le premier volet
porte sur une tude exprimentale base sur la rvlation de dfauts internes par des techniques de contrle
non destructif (mesures de champs acoustiques, thermiques et densimtriques), visant dterminer leur
position et leur taille. Le deuxime volet prsente une dmarche didentification numrique permettant
dvaluer les caractristiques mcaniques rsiduelles dune structure composite constitue dune zone saine
et dune zone endommage. Cette thse sinscrit dans une dmarche qualit de dveloppement dun outil
daide la dcision visant dune part, la dtection et la localisation des dfauts internes, et dautre part,
lvaluation de leur niveau de dgradation.
Mots-cls : composite, contrle non destructif, mesures de champs, ultrasons, thermographie infrarouge,
tomographie X, identification, recalage lments finis
Title : Experimental identification of the behavior of a composite fuselage : defect detection by full field
measurements
Abstract :
These works are related to the identification of the behavior of composite structure. The first part concerns
an experimental study based on internal defect revelation by non destructive testing methods (acoustic,
thermal, densimetric full field measurements) to determine their location and their size. The second
part presents a numerical identification method to evaluate the four in-plane elastic parameters of an
orthotropic composite structure (sample with a healthy zone and a damaged zone simulated by a flat
bottomed hole). This thesis is part of a quality approach in order to develop of a decision support tool
with firstly, the detection and the localization of internal defects, and secondly, the evaluation of their
damage level.
Keywords : conposite, non destructive testing, full field measurements, ultrasound, infrared thermogra-
phy, X ray tomography, experimental identification, finite element model updating
Laboratoire Gnie de Production de lcole Nationale dIngnieurs de Tarbes
47, Avenue dAzereix, BP 1629, 65016 TARBES Cedex

Peronnet

Transféré par

Droits d'auteur :

Formats disponibles

Peronnet

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Peronnet

Transféré par

Droits d'auteur :

Formats disponibles

%0$503"5%&-6/*7&34*5%&506-064&

Institut National Polytechnique

Je remercie galement Moussa Karama, mon directeur de thse, de sa disponibilit

Je veux galement remercier Marie-Laetitia Pastor et Olivier Dalverny qui ont eu un

Je remercie Serge Cristante de stre toujours souci de la valorisation industrielle de

Table des matires i

Table des figures v

Liste des tableaux xi

1 Contrle de dfauts par mesures de champs acoustiques, densimtriques,

2 Matriaux et protocoles exprimentaux 41

2.2.1 Eprouvette 1 : matriau monolithique avec inserts Tflon . . . . . 43

3 Etude comparative sur prouvettes talons 75

4 Identification des dfauts par mesures de champs cinmatiques 111

4.5.1 Modles de rfrence et de recalage . . . . . . . . . . . . . . . . . . 129

1.1 Principe gnral dun contrle non destructif . . . . . . . . . . . . . . . . 6

1.19 Gnrations de scanners (Lenoir [134]) . . . . . . . . . . . . . . . . . . . . 26

2.1 Photo du matriau monolithique avec inserts Tflon . . . . . . . . . . . 43

2.15 Moule rigide en verre-poxy . . . . . . . . . . . . . . . . . . . . . . . . . . 51

2.41 Dispositif exprimental de la TIRP en configuration face arrire . . . . . . 72

3.1 Procdure de traitement dimage automatique . . . . . . . . . . . . . . . . 78

4.1 Sytme possdant deux camras (Fazzini [76]) . . . . . . . . . . . . . . . . 116

4.6 Problme considr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126

3.1 Rsultats du traitement dimage par famille de mesures de champs . . . . 79

4.1 Donnes de base et rsultats des problmes directs et indirects (Kleiner-

Le dveloppement de nouveaux procds de fabrication pour les composites constitue

Le thme principal de la thse prsente ici est lidentification de dfauts internes au

La rdaction de ce mmoire sorganise autour de quatre parties.

La prsentation des matriaux et des protocoles exprimentaux employs fait lobjet du

Le troisime chapitre est consacr ltude comparative en elle-mme. A partir de

Contrle de dfauts par mesures de

Gnralement, une technique de CND se dcompose en trois phases (Boro-Djordjevic

Figure 1.1 Principe gnral dun contrle non destructif

induits en service (rupture de fibres, fissures, dlaminages, dcollements fibres-matrice,

Toutes ces mesures sont bases sur le mme principe gnral :

avant) si le dtecteur se trouve du mme ct que lmetteur 1 , soit une mesure en

1.2 Mesures de champs acoustiques : les ultrasons

Suivant la dmarche gnrale expose prcdemment, le contrle ultrasonore suit les

1.2.1 Principe physique

La frquence choisie doit satisfaire deux condions ncessaires au respect de lhypothse

Le faisceau acoustique du traducteur se dcompose classiquement en deux parties :

1.2.1.2 Les diffrents types dondes

Suivant lorientation du capteur (do implicitement celle du transducteur), plusieurs

les traducteurs dangles permettent de produire des ondes de volume transversales,

Les ondes de volume longitudinales dites ondes de compression se caractrisent par un

1.2.1.3 Lattnuation ultrasonore

1.2.2 Principe de contrle

lcho dmission reprsentant le lieu dmission du train dondes,

Figure 1.8 Diffrents modes de prsentation de linspection ultrasonore (Dumont-Fillon

1.2.2.2 Transmission ultrasonore : contrle en face avant et face arrire

La transmission ultrasonore utilise deux traducteurs, lun jouant le rle dmetteur et

Figure 1.9 Contrle en mode transmission (Chaix et al. [53])

Figure 1.10 Principe de dtection par transmission (Lacroix [119])

1.3 Mesures de champs densimtriques : radiographie

1.3.1 Principe physique

Les rayons X sont des rayonnements lectromagntiques de trs forte nergie, de

X avec la matire conduit la production dlectrons responsables dionisations dans le

Figure 1.11 Spectre lectromagntique (Garnier [83])

1.3.1.2 Loi dattnuation des rayons X

interaction avec les atomes du milieu (effet photolectrique),

%0$503"5%&-6/7&345%&506-064&