Instrumentation et Méthodes pour le contrôle non destructif

Illustration au contrôle des conduites

de précontrainte des ouvrages d’art

Intervenant : Thierry Bore

Dominique Placko – Frédéric Taillade

séminaire Centre d'Alembert - 21 avril 2010

 Plan de la présentation

-  Généralités sur le Contrôle Non Destructif (CND)

-  Précontrainte extérieure : présentation, défauts

rencontrés et procédure d’auscultation

- Sonde capacitive: fonctionnement et performances

-  Modélisation par DPSM

- Application de la technique à notre problème et modélisation

de plans d’expérience
 Généralités sur le CND
Le CND c’est quoi ?

Les Contrôles Non Destructifs sont un ensemble de méthodes

qui permettent l'évaluation de l'état d'intégrité d'un matériau, de
structures, ou d'ensembles mécaniques, sans les dégrader

Ils peuvent être appliqués pratiquement à tous les stades de vie

d'un élément:
- production
- utilisation
- étude de vieillissement
- rétrofit
- maintenance
- étude de vieillissement

Il existe de nombreuses techniques:

endoscopie, ressuage, ultrasons, radiographie, tomographie rayons X,
courant de foucault, magnétoscopie, thermographie, shearographie ...
 Généralités sur le CND

Quelle que soit la méthode adoptée, on peut représenter la mise en

œuvre d'un système CND suivant le synoptique suivant :

La cible se caractérise par un ensemble de paramètres que l'on va

chercher à estimer afin de former un diagnostic d'intégrité

La mise en œuvre d'un système CND adéquat va permettre de produire un

certain nombre de signaux qui sont fonction des paramètres recherchés

Une étape « d'inversion », plus ou moins compliquée, est bien souvent

nécessaire afin de retrouver les paramètres initiaux de la pièce
 Précontrainte extérieure ?
De nombreux ouvrages utilisent la précontrainte externe, que ce soit
dès leur construction ou comme procédé de renforcement de la
structure
 Précontrainte extérieure ?

Intérieur du pont

conduits

 Conduit accessible donc contrôlable !

 Conduit de précontrainte extérieure

coupe d’un conduit

Conduit: PEHD (diamètre ext: 90mm)

Câbles d’acier

Coulis de ciment
 Défauts constatés

•  Saint Cloud 1998 Viaduc de Saint Cloud 1998

–  un câble de renforcement mis en
place en 1981 est trouvé rompu
par un agent de la Subdivision
Entretien de l’A13

•  Pont-A-Mousson 2005
–  Localisation : zone inclinée proche
d’un déviateur en point haut du
câblage
–  Corrosion avec réduction de Viaduc de Pont-A-Mousson 2005
section pouvant atteindre 50% de
la section initiale
–  Absence de coulis dans la gaine,
pâte blanche
–  Analyse du coulis d’injection
•  Forte teneur en fer
•  Apparition du manganèse et du chrome
(éléments d’alliage de l’acier)

 Corrosion des aciers

 Corrosion des aciers

Présence de vides
Processus de corrosion « normaux » par condensation
d’eau sur l’acier par exemple

Présence de ségrégation/ressuage
Création d’une pile électrochimique acier dans
coulis sain / acier dans phase ségrégée/ressuée

Sensibilité de l’acier vis-à-vis de certains paramètres

Formulations
Eau – ciment – adjuvant Chryso GT

E/C = 0,65 (sans l'eau de l'adjuvant)

Extrait sec du chryso GT: 30%
 pH de la pâte blanche mesurée sur le terrain 13,9 !
E. Blactot &T. Chaussadent JOA 2007
 Procédure d’auscultation
Procédures d’auscultation inappropriées:
-  destructives : ouverture de la gaine et caméra endoscopique
-  coûteuses et peu commodes : rayon gamma
-  extrêmement qualitative : sondage au marteau

Méthode de référence !

Nécessité de développer de nouveaux outil d’investigation

automatisables et délivrant des informations quantitatives
 Sonde Capacitive

électrodes
sonde de 3eme génération
 Fonctionnement de la sonde

Support
Electrodes

Intérieur du conduit:
Matériaux inconnus

conduit (PEHD)

Oscillator
Frequency
C(ε) Circuit
Divider
≈ 65 MHz

≈ 13 kHz
 Performance de la sonde (1/5)

Mesure sur le terrain:

SETRA-RST :
Campagne d’auscultation des
ouvrages en béton précontraint

Pont de Labéraudie à Cahors

Mesures réalisées avec le LR Autun

 Performance de la sonde (2/5)
Mesure sur la génératrice supérieure
Travée centrale câble A trajet longitudinal départ pile RD vers pile RG

Fréquence

Abscisse
 Performance de la sonde (3/5)
Mesure circulaire:
Travée centrale câble A trajet circulaire (2 tours) à 1 m de la pile RD

Vide
Fréquence (Hz)

Abscisse
 Performance de la sonde (4/5)

Ouverture du conduit à posteriori:

 Performance de la sonde (5/5)
Performance de la sonde actuelle:
Permet la détection de défauts:
-  vide
-  pâte blanche
-  distingue un retrait d’un remplissage partiel

Limitations de la sensibilité:

-  en fonction de la permittivité: pas sensible à la présence d’eau

-  spatiale: limitée au voisinage immédiat des électrodes donc
fortement sensible à la nature et la géométrie du conduit
(variation d’épaisseur de la gaine)
 Choix d’une technique de modélisation

Besoin d’un modèle complet pour:

-  prendre en compte la géométrie du conduit et l’existence de

plusieurs milieux à l’intérieur du conduit, développer une
modélisation 3D
-  interpréter les signaux donnés par notre capteur
-  élaborer des stratégies d’inversion qui permettront d'extraire,
à partir des données, des renseignements précis et attendus
par les gestionnaires des ouvrages sur l'état des conduits de
précontraintes

Nous avons choisi d’utiliser un modèle semi analytique,

le modèle DPSM

DPSM = Distributed Point Sources Method

 Principe général de la méthode DPSM

Les surfaces actives (surfaces de transducteurs par

exemple) sont représentées par un ensemble de
sources ponctuelles situées hors de la surface et
dont l’amplitude et le diagramme de rayonnement
sont calculés pour satisfaire les conditions aux
limites du problème
Exemple de source DPSM ultrasonore isochrone

Une source ponctuelle sinusoïdale de fréquence

f=1MHz, vitesse radiale 1 m/s, dans un milieu de
densité 1 et de celérité 1490m/s et rayonnant selon
un diagramme de rayonnement de 240°
Exemple de source DPSM ultrasonore impulsionnelle

Une source ponctuelle impulsionnelle de fréquence

f=1MHz, vitesse radiale 1 m/s, dans un milieu de
densité 1 et de celérité 1490m/s et rayonnant selon
un diagramme de rayonnement de 240°
 La méthode DPSM

Qu'est-ce que la méthode DPSM ?

Il s'agit d'une méthode générique de modélisation

3D des interactions capteur/cible, applicable au
cas où les équations de comportement sont
connues dans les différents milieux
Cette méthode se fonde sur une distribution
spatiale de sources ponctuelles, qui ‘synthétisent’
les grandeurs physiques (champs électriques,
magnétiques, etc) dans les différents milieux
Développée à l'origine pour des problèmes de CND
de nature électromagnétique et ultrasonore, son
utilisation devient peu à peu générique
 Influence des interfaces: illustration

Synthèse par superposition d’un ensemble de

sources DPSM d’un champ de pression acoustique
transmis à la traversée d’une interface
 Modélisation d’une interface
La synthèse des grandeurs réfléchies et transmises au niveau
de l’interface est assurée par des couches de sources situées
de part et d’autre, et calculées pour satisfaire les conditions
aux limites sur un ensemble de points tests, du problème pris
dans son ensemble.

Chaque ‘bulle’ contient une ou plusieurs sources ponctuelles:

charges électrostatiques, éléments de courants ix, iy, iz, etc..
 Illustration du rôle des interfaces

Illustration avec deux milieux séparés par une interface:

Une configuration particulière avec deux milieux identiques permet
de tester les capacités des sources virtuelles à synthétiser les
grandeurs de part et d’autre de l’interface.
 Raccordement aux interfaces

Quelles grandeurs raccorder aux interfaces ?

Pour développer le caractère générique de la méthode, on cherche
à raccorder les grandeurs adaptables à tous les problèmes ( US,
EM, électrostatique, magnétostatique…): le potentiel et sa dérivée,
qu’il soit scalaire ou vectoriel conviennent parfaitement…

électrostatique

ultrasons

électromagnétisme
 Construction de la matrice de résolution (1/2)
Une fois cette étape franchie, on peut faire apparaître les matrices
de couplage entre les différents éléments qui constituent le problème

Les matrices d’auto couplage: les sources d’un objet rayonnant

sur les points tests de cet objet

Les matrices d’inter couplage: les sources d’un objet rayonnant

sur les points tests d’un autre objet

2
1
 Construction de la matrice de résolution (2/2)
Le problème initial, complexe, est réduit à un ensemble de
matrices élémentaires décrivant les relations entre les différents
objets et permettant d’évaluer les grandeurs que l’on souhaitera
raccorder
Ensuite, on agence l’ensemble de ces matrices pour former une
matrice globale de résolution M

auto couplage

inter couplage

La matrice M relie le vecteur contenant l’ensemble des valeurs des

sources DPSM J (inconnues à ce stade), à un vecteur de conditions
aux limites C, intrinsèques ou imposées par l’utilisateur:
C=M*J
 Formalisme DPSM en électrostatique
Pour un point cible et Ns sources:

Pour Np points cibles et Ns sources:

matrice de couplage
pour le potentiel

matrice de couplage
pour la composante
selon z du champ
Electrostatique: exemple de géométrie avec Interfaces
Condensateur plan dans lequel est inséré deux diélectriques::

sources α1

points test a1
sourcesβ2

axe z points test b1

sources δ1 =
points test d2 points test b2

sources γ2

axe y
 Electrostatique: calcul de la Matrice globale
Surfaces actives: Expressions des conditions limites ‘utilisateurs’

fixé par
l’utilisateur
Interfaces: Expressions des conditions limites ‘intrinsèques’:
Continuité du potentiel V
Continuité de la composante normale de D

Formulation matricielle du problème:

INCONNUES DU PROBLEME
 Electrostatique

potentiel et composante Ez
calculés dans le plan y0z
 Electrostatique

potentiel et lignes de champ

calculés dans le plan y0z
 Les avantages “concurrentiels”
Quels sont les atouts de la méthode DPSM ?

Rester proche de la physique du problème et

obtenir un modèle semi analytique sous forme
matricielle
Absence de maillage volumique (3D) , donc...
Temps de calcul court: possibilité de réaliser des
calculs paramétrés (animations, temps quasi réel,
modélisation de plans d’expériences)
Possibilité de séparer les effets des sources liées aux
différents objets (extinction des sources du capteur
inducteur par exemple)

Simulation multiphysique simultanée et

impulsionnelle possible
 Exemple modélisation DPSM (1/3)
Exemple électromagnétisme (1/2)
Air – sphère diélectrique (εr = 15), f= 500 MHz
Mode isochrone

spire de courant
dans un plan xOy

I(t)
V(t)

milieu 1
milieu ambiant
champ magnétique
induit

milieu 2
obstacle sphérique
milieu diélectrique
 Exemple modélisation DPSM (2/3)
Exemple électromagnétisme (2/2)
CND par Courant de Foucault

Inclusion sphérique dans un matériau conducteur

f=4.5 kHz
 Exemple modélisation DPSM (3/3)
Exemple en acoustique

Interface air – eau, f=200 kHz

Mode impulsionnel
 Illustration de la ‘séparation des effets’
Eau –bulle d’air, f=1 MHz

Visualisation pression ‘totale’’ Visualisation onde rétro diffusée

et onde à l’intérieur de la bulle
 Retour au problème
Quelle géométrie modéliser ?
Problème rencontré lors de l’injection::
 Application de la technique DPSM à notre problème

Le laboratoire dispose de corps d’épreuve à défaut maîtrisés:

Conduit sans toron métallique:
l’objectif est pour l’instant de
détecter la présence de pâte
Blanche et de vide

Un étalonnage permet d’obtenir

la capacité en fonction
de la fréquence mesurée

Conduit auquel on applique la méthode:

Maillage DPSM du conduit
 Représentation des grandeurs

Représentation des grandeurs dans

un plan transverse

Les propriétés des milieux étant

identiques: permet d’illustrer la
continuité aux interfaces
 Calcul de la capacité

Soient deux conducteurs en influence:

L’équilibre de deux conducteurs en influence est défini par:

les cij sont des coefficients ne dépendant que de la configuration

géométrique du problème: en utilisant des excitations différentes
on peut calculer les 4 coefficients et donc la capacité
 Capacité en fonction de la position
Calculs de capacité réalisés en rotation
Allure des courbes obtenues : signature des milieux
 Comparaison mesures/simulations

Mesures réalisées sur le conduit d’essai à z = 500 mm

Modélisation aux EF (réalisées avec un logiciel commercial) : modélisation 2D
 Utilisation du modèle dans un plan d’expérience
Réalisation d’une étude paramétrique
OBJECTIF : comprendre l’influence de la nature des matériaux à
l’intérieur du conduit sur la capacité

Matériau inconnu, εr =X

Matériau inconnu, εr =Y

ciment

gaine (PEHD)

On réalise une série de calculs

en faisant varier la valeur de X et Y
 Sensibilité aux permittivités (1/2)

Représentation de la capacité pour une position en fonction de X et Y:

Permet de retrouver des résultats sur la sensibilité de la

sonde:
pas ou peu sensible à l’eau (εr = 80 pour l’eau pure)
 Sensibilité aux permittivités (2/2)

Réalisation d’un ‘fit’ des courbes

Les paramètres de notre ‘fit’ constituent les nouveaux paramètres

de notre problème direct: on pourra tenter de l’inverser
 Estimation des permittivités

Bonne décorrélation des estimations selon X et Y

Robustesse vis-à-vis du bruit additif gaussien

Néanmoins, erreur pour les faibles valeurs de Y

 Utilisation du modèle dans un plan d’expérience
Réalisation d’une étude paramétrique
OBJECTIF : comprendre l’influence de l’épaisseur du menisque
d’air et de la couche de pâte

air : épaisseur e1

pâte blanche: épaisseur e2

ciment

gaine (PEHD)

On réalise une série de calculs

en faisant varier la valeur de e1 et e2
 Sensibilité aux épaisseurs

Représentation de la capacité pour une position en fonction de

e1 et e2 :
 Estimation des épaisseurs

Bonne décorrélation des estimations selon e1 et e2

Robustesse vis-à-vis du bruit additif gaussien

Estimation de e1 (épaisseur de pâte) de moins bonne

qualité
 Conclusion

Pour la partie instrumentation

Identification des phénomènes de corrosion dans

les conduits de précontrainte extérieure
Développement d’une sonde capacitive pour
inspecter les conduits
La sonde capacitive nous permet d’identifier la
présence d’un vide ou la présence de pâte blanche

L’efficacité de la sonde capacitive a été démontrée

sur le terrain: utilisée par les équipes d’inspection
 Conclusion

Pour la partie modélisation

Bonne adéquation entre la modélisation et les

mesures
Possibilité de simuler des plans d’expérience

Pour la partie information & traitement du signal

Estimation des données internes de la gaine via une

mesure de capacité possible & robuste
 Perspectives

Prochaines étapes
Estimation des données internes de la gaine à partir
de données obtenues en laboratoire ou sur le
terrain

Enrichir l’information délivrée par la sonde pour

estimer plus de paramètres

-Sonde multi électrode ?

-Mesures multi fréquence ?
-Multi physique ?
 Perspectives pour la modélisation DPSM

Perspective d’évolution de la méthode DPSM:

Application de la méthode en mécanique des fluides

Devellopper la méthode en thermique pour des

applications dans le domaine du bâtiment

Coupler la méthode avec des “mailleurs” de type

Industriel (CATIA par exemple)
Merci de votre attention
 Eléments de bibiographie
- D. PLACKO, Ecole Normale Supérieure, Cachan, France, [Link], University of Arizona, USA:
"DPSM - Distributed Point Source Method for Modeling Engineering Problems"
372 pages, juin 2007, éditeur John Wiley & Sons

- D. PLACKO, Ecole Normale Supérieure, Cachan, France, T. KUNDU, University of Arizona, USA:
"Modeling of Ultrasonic Field by Distributed Point Source Method"
Chapitre 2 de l'ouvrage: "Ultrasonic Nondestructive Evaluation: Engineering and Biological Material Characterization",
pp 143-202
Coordonateur: [Link], éditeur: CRC press, USA, janvier 2004

- [Link], T. KUNDU, University of Arizona, USA, [Link], Ecole Normale Supérieure, Cachan, France:
"Modeling of the Ultrasonic Field of Two Transducers Immersed in a Homogeneous Fluid using Distributed Point Source
Method"
Chapitre 5 de l'ouvrage: " Advanced Ultrasonic Methods for Material and Structure Inspection", pp 159-187.
Coordonateur: [Link], éditeur Hermes Science Publishing, UK, janvier 2007

- [Link], [Link], [Link]:

"Procédé pour évaluer une grandeur physique representative d'une interaction entre une onde et un obstacle".
Dépôt de Brevet N° 02 14108 ENS Cachan/CNRS/Université d'Arizona (novembre 2002). Extension PCT (novembre
2003).

[Link], [Link], [Link], [Link]:

"Procédé universel de modélisation des interactions entre au moins une onde et au moins un objet, la surface de chaque
objet définissant une interface entre au moins deux milieux".
Dépôt de Brevet français N° 05/13219 (décembre 2005) et extension PCT N° 63831 ENS Cachan/CNRS/Université
d'Arizona, (décembre 2006). Extension Canada et Etats-Unis en cours (janvier 2009).