Modélisation électromagnétique : détermination

de la corrosion à l’aide d’une bobine d’inductance

Simulations par éléments finis d’une bobine d’inductance..

RAPPORT D’ÉTUDE

2021

© copyright photo
 Le Cerema (Centre d’études et d’expertise sur les risques, l’environnement, la mobilité et
l’aménagement) est un établissement public tourné vers l’appui aux politiques publiques,
placé sous la double tutelle du ministère de la transition écologique et du ministère de la
cohésion des territoires et des relations avec les collectivités territoriales.
Les métiers du Cerema s'organisent autour de 6 grands domaines d'activité complémen-
taires visant à accompagner les acteurs territoriaux dans la réalisation de leurs projets.
6 domaines d’activité :
Expertise et ingénierie territoriale / Bâtiment / Mobilités / Infrastructures de transport
/ Environnement et risques / Mer et littoral

Site web : cerema.fr

2/24 Titre - Sous-titre /mois année

Modélisation électromagnétiques :
détermination des endommagements
dans les câbles

Compléments des simulations par éléments finis, par rapport à la

Norme

Commanditaire :

Auteur :

Responsable du rapport

Prénom NOM – Département xxx xxx – Unité xxx xxx

Tél. : +33(0)X XX XX XX XX – Fax : +33(0)X XX XX XX XX

Courrier : pré[email protected]

Direction - adresse

Historique des versions du document

Version Date Commentaire

Références

N° d’affaire : XXXXXXXXX
Partenaire(s) : XXXXXXXX
Devis n° xx/xxxxx du jj/mm/aaaa

Nom Service Rôle Date Visa

Auteur principal

Contributeur

Titre - Sous-titre /mois année 3/24

 Contributeur

Relecteur

Relecteur

4/24 Titre - Sous-titre /mois année

Résumé de l’étude
Cette étude fait suite aux travaux précédents sur les modélisations des dispositifs d’auscultation ma-
gnéto-inductifs à l’aide des boucles magnétiques et aimants permanents.

5 à 10 mots clés à retenir de l’étude

Statut de communication de l’étude

Les études réalisées par le Cerema sur sa subvention pour charge de service public sont par défaut
indexées et accessibles sur le portail documentaire du Cerema. Toutefois, certaines études à caractère
spécifique peuvent être en accès restreint ou confidentiel. Il est demandé de préciser ci-dessous le
statut de communication de l’étude.

Accès libre : document accessible au public sur internet

Accès restreint : document accessible uniquement aux agents du Cerema

Accès confidentiel : document non accessible

Cette étude est capitalisée sur la plateforme documentaire CeremaDoc, via le dépôt de document :
https://doc.cerema.fr/depot-rapport.aspx

Titre - Sous-titre /mois année 5/24

Contexte et objet de l’étude

Cette étude vient à la suite de l’étude sur l’endommagement des câbles à l’aide des
méthode électro-inductive.

Elle fait aussi suite aux travaux de conception et de fabrication de bobine d’inductance
par l’ex-CECP de la Direction Territoriale Normandie Centre du CEREMA. Elle vise
donc la détectabilité des défauts par la variation d’une caractéristique d’une bobine
alimentée à l’aide d’un courant alternatif. Il s’agit donc de calculer les variations de
l’inductance dans les bobines réceptrices, générés par la présence d’un défaut dans
un câble métallique.

6/24 Titre - Sous-titre /mois année

Titre - Sous-titre /mois année 7/24
SOMMAIRE
1. RAPPEL SUR L’EXPERIMENTATION 11
1.1 Rappel sur la Norme EN 12927 .......................................................................... 11
1.2 Dispositif expérimental ........................................................................................ 11
1.2.1 Principe de la détection de la corrosion :............................................................... 13
1.3 Calcul de l’inductance : .......................................................................................13

2. SIMULATION DE LA BOBINE DE CORROSION 14

2.1 Modélisation géométrique de la bobine avec CST : ............................................14
2.1.1 Résultat de la simulation : ..................................................................................... 14
2.2 Modélisation avec COMSOL: ..............................................................................16
2.2.1 Modélisation géométrique : ................................................................................... 16
2.2.2 Simulation physique du phénomène : ................................................................... 17
2.3 Résultat et analyse de la simulation :..................................................................18

3. 20

4. ERGO. 21
4.1 Titre 1.1 - e.................................................................. Erreur ! Signet non défini.

5. BIBLIOGRAPHIE 22
5.1 Bibliographie .......................................................................................................22
5.2 Glossaire.............................................................................................................22

8/24 Titre - Sous-titre /mois année

Titre - Sous-titre /mois année 9/24
Introduction

Les dispositifs expérimentaux grâce aux bobines d’oxydation, permettent la détection

de la détérioration d’un câble par la variation de la mesure de l’inductance.

Nous ne nous intéresserons dans cette étude qu’aux champs magnétiques et courants
induits et à l’énergie calculées générés dans les bobines.

La modélisation numérique est faite grâce au logiciel de simulation par élément fini,
Comsol Multiphysics, qui permet de résoudre des problèmes en 2D et 3D par la mé-
thode des éléments finis.

10/24 Titre - Sous-titre /mois année

1. RAPPEL SUR L’EXPERIMENTATION

L’expérimentation consiste à utiliser le dispositif mis au point par l’ex CECP du DERDI
Grand-Quevilly.

1.1 Rappel sur la Norme EN 12927

Il s’agit d’une norme de « Prescriptions de sécurité des installations à câbles transpor-

tant des personnes », notamment au niveau des câbles

Par souci de sécurité et de prévention par rapport aux risques de ruptures des ponts
suspendus, les câbles de suspentes doivent être contrôles et inspecter.

1.2 Dispositif expérimental

Le câble expérimental testé est placé à l’intérieur de la bobine de 96 spires et parcouru

par un courant de 100 mA d’intensité.

La bobine suivante est conçue pour des câbles de 70 mm de diamètre.

Titre - Sous-titre /mois année 11/24

Figure 1 : la bobine d'induction fermée pour des câbles de 70 mm.

La bobine est composée de 96 spires et est alimentée par un courant sinusoïdal alter-
natif.

Figure 2 : bobine composée de 96 spires.

12/24 Titre - Sous-titre /mois année

Cette bobine d’inductance, ou aussi appelée, bobine de corrosion a été fabriquée par
l’ex Centre d’Etude de Prototype du Cerema NC.

1.2.1 Principe de la détection de la corrosion :

L’inductance, présente dans les bobines à courant alternatif, est une grandeur qui per-
met de quantifier l’opposition de la bobineau passage du courant, on parle aussi de
réactance inductive. Le défaut est détecté lorsque la valeur de l’inductance varie par
rapport à la valeur d’origine.
Notons que la forme, le diamètre de sa section ainsi que le nombre de spires font varier
la valeur de l’inductance.

1.3 Calcul de l’inductance :

Dans le cas d’une bobine comprenant N spites, à vide, de longueur l et de surface S, l’expres-
sion de l’inductance peut être simplifiée par rapport aux formules mathématiques complexes
la définissant:
𝑁²𝑆
𝐿 = µ0
𝑙
L = inductance est en Henry

µ0 = 4π x 10-7 Tm/A est la perméabilité magnétique du vide.

l est la longueur du « noyau » de la bobine en mètre (m).
N le nombre de spire
S la section de la bobine en mètre carré (m²).

Si la bobine contient un noyau à l’intérieur, comme un câble ou une barre en acier :

𝑁²𝑆
𝐿 = µ0 µ𝑟
. 𝑙
µr est la perméabilité relative du matériau magnétique dans la bobine, le noyau.
Application numérique :

Titre - Sous-titre /mois année 13/24

Nous pouvons ainsi calculer les valeurs de l’inductance pour la bobine modélisée par la suite,
avec : - le nombre de spire est de N = 19
- la longueur de la bobine l = 0,23 m
- le diamètre de la bobine d’inductance D = 0,025 m
- la valeur de la permittivité relative de l’acier du câble µr = 4000.
Les valeurs de l’inductance à vide L vide et avec le câble L câble sont :
L vide = 9,68 . 10-7 H
L câble = 38,72 10-4 H.
Nous devons comparer ces valeurs avec celles données par nos simulations.

2. SIMULATION DE LA BOBINE DE CORROSION

2.1 Modélisation géométrique de la bobine avec CST :

La boucle émettrice est simulée avec un courant électrique sinusoïdal, dont l’intensité,
de valeur efficace de 100 mA, oscille avec une fréquence sinusoïdale de 10 KHz.

Cette simulation s’est effectuée avec deux modèles fournies par CST mais qui ne cor-
respondent pas aux dimensions de notre bobine.

2.1.1 Résultat de la simulation :

14/24 Titre - Sous-titre /mois année

Figure 3: résultat de la simulation, calcul et visualisation du courant J avec un diamètre de 0,02

Figure 4 : Courbe de la variation de l’Inductance en fonction du diamètre du câble dans la bobine.

On constate que les points sont placés exactement sur une courbe d’interpolation po-
lynomiale du troisième degré.
Les valeurs de l’inductance calculées par CST sont plus élevées que les valeurs théo-
riques. Cela est sûrement dû aux différences de dimensions entre les bobines.

Titre - Sous-titre /mois année 15/24

2.2 Modélisation avec COMSOL:

2.2.1 Modélisation géométrique :

Le dispositif expérimental comprend une bobine d’induction alimenté par un courant

alternatif sinusoïdal.

Géométriquement comme un câble plein de 0, 045 m diamètre et de 0,8 m de lon-

gueur.

Figure 2: Modélisation géométrique de la bobine d’inductance à vide câble.

Figure 2: Modélisation géométrique du câble, de la bobine réceptrice et du défaut.

16/24 Titre - Sous-titre /mois année

Figure 2: Modélisation géométrique de la bobine avec le câble comme noyau.

2.2.2 Simulation physique du phénomène :

Nous utilisons le modèle Magnétostatique quasi-statique, qui permet de simuler des

phénomènes magnétiques avec des courants alternatifs à basse fréquence.

Le dispositif d’excitation, une bobine émettrice à basse fréquence de 1 kHz, permettant

de magnétiser le câble, est modélisé par un dipôle magnétique, situé au centre du
câble.

Le câble est en acier et sa magnétisation suit une loi de type B = B(H). nous utilisons
les courbes dont disposent les propriétés des matériaux magnétisables dans Comsol.
La courbe B – H est généralement utilisé pour décrire le degré de magnétisation qu’un
matériau ferromagnétique obtient en réponse d’un champ magnétique appliqué.
En d’autres mots, la courbe caractérise la perméabilité µ d’un matériau magnétique
non linéaire, définit par :
:
2 ∗ 𝑊𝑖𝑛𝑡𝑀𝑎𝑔𝑛
𝐿=
𝐼𝑐𝑜𝑖𝑙 ²

L est l’inductance s’exprime en Henry (H),

WintMagn est la densité de l’énergie magnétique en J,

Titre - Sous-titre /mois année 17/24

et Icoil est la norme du courant exprimé en A/m², parcourant la boucle.

L’algorithme permet de décrire la courbe d’un matériau magnétique non-linéaire

jusqu’à la saturation suivant la valeur de saturation désirée, ici 2,4 T, suivant une équa-
tion du type :
𝐵(𝐻) = µ0 µ𝑟 .H
, où B est la densité de flux magnétique, appelé couramment champ magnétique, en
T et H est l’excitation magnétisation en A/m ,
μ0 , la perméabilité magnétique du vide vaut μ0 = 4π.10-7 en H.m-1 et
µr , la perméabilité relative est ici prise égale à µr = 4000.
C’est cette valeur de µr qui est utilisée pour la saturation.

2.3 Résultat et analyse de la simulation :

Pour un diamètre variable de la barre en acier inséré au centre dans le dispositif avec
la bobine, nous avons évalué les valeurs de l’inductance :

Pour un diamètre de D = 0,02

Répartition du champ magnétique à l’intérieur du tube :

18/24 Titre - Sous-titre /mois année

Figure 2: Valeur de l’énergie magnétique dans la bobine.

Les valeurs trouvées par la simulation sont sur une courbe d’interpolation polynomiale
du troisième degré.

Titre - Sous-titre /mois année 19/24

20/24 Titre - Sous-titre /mois année
3. CONCLUSION.

Les valeurs trouvées par les simulations et le calcul théorique sont assez différentes. Pour
CST, cela est dû au fait que les dimensions de la bobine et du câble n’ont pas pu être respec-
tées, notre période d’évaluation du logiciel CST étant trop courte. Il convient donc que nous
puissions redimensionner la bobine et le câble comme avec l’expérimentation.

Titre - Sous-titre /mois année 21/24

4. BIBLIOGRAPHIE

4.1 Bibliographie
Pour citer un ouvrage
(auteur) Nom, Prénom. Titre de l’ouvrage. Lieu d’édition. Éditeur, année de publication,
nombre de pages.
Pour citer un article papier
« Titre de l’article », Titre de la revue, année de publication/numéro de tome (si numérotation
annuelle, ex : 2012/4) ou numéro continu (ex : n° 100), p. xx-xx
Pour citer un article disponible en ligne
« Titre de l’article », Titre de la revue, année de publication/numéro de tome (si numérotation
annuelle, ex : 2012/4) ou numéro continu (ex : n° 100), p. xx-xx. DOI s’il s’agit d’un article
scientifique qui en dispose. URL : lien vers la page
Pour citer un site web

COMSOL Multiphysics 5.0 « Modeling of a 3D Inductor »

(auteur) Nom, Prénom. Titre de la page web ou du site web consulté. DOI s’il s’agit d’un article
scientifique qui en dispose. URL : lien vers la page, [consulté le jj/mm/aaaa]
Pour citer une base de données
Titre de la base de données. Editeur, URL si la base est consultable en ligne : lien vers la
page, [Consulté le jj/mm/aaaa]

4.2 Glossaire
DtecTV : Cerema, Direction technique Territoires & ville
DtecREM : Cerema, Direction technique Risque, eau et mer
DtecITM : Cerema, Direction technique Infrastructures, transport & matériaux

22/24 Titre - Sous-titre /mois année

 Centre d’études et d’expertise sur les risques, l’environnement, la mobilité et l’aménagement
Siège social : Cité des mobilités - 25 avenue François Mitterrand - CS 92 803 - 69674 Bron Cedex -
Tél. : +33 (0)4 72 14 30 30 – www.cerema.fr