REPUBLIQUE DU CAMEROUN REPUBLIC OF CAMEROON

UNIVERSITE DE DOUALA UNIVERSITY OF DOUALA

******** *******
FACULTE DE GENIE INDUSTRIEL FACULTY OF INDUSTRIAL ENGINEERING
[Link] [Link]

Département de Robotique Industrielle (ROI)

NIVEAU 4

OPTION : Ingénierie électrique et électromécanique

Conception Assisté par Ordinateur

(CAO)

GUIDE DE TRAVAUX PRATIQUE

Année académique
2018 / 2019
Introduction
Ce document est un guide simplifié de travaux pratique sur le logiciel de simulation Comsol
multiphysic. Le terme logiciel de simulation désigne l’outil informatique qui permet de simuler
numériquement un modèle. Des exemples de bases traitant du domaine du génie électrique y
sont développés ainsi que les instructions permettant leurs simulations. Ces travaux pratiques
permettront aux lecteurs étudiants d’asseoir certaines notions de base sur la simulation
numérique de phénomènes physiques et multi-physiques. D’acquérir également des aptitudes
certaines dans la prise en main du logiciel Comsol multiphysic.

La version 5.1 de Comsol multiphysic a été utilisée pour tous les exemples illustrés dans le
document car elle est largement suffisante pour les problèmes les plus complexes rencontrés
durant les études de génie.

La démarche de modélisation et de simulation est une approche scientifique qui a pour but de
produire de la connaissance sur des phénomènes naturels ou artificiels. Cette démarche est
utilisée pour comprendre, prédire voire contrôler des objets ou des phénomènes qu’ils soient
déjà existants ou encore au stade de la conception. Le principe est de créer une simplification
(un modèle) du phénomène ou de l’objet étudié, par exemple une maquette réduite, un
prototype ou un modèle numérique. En étudiant le comportement du modèle dans différentes
situations il est possible de répondre aux questions que l’on se pose sur le phénomène initial.
Une telle démarche trouve son intérêt lorsqu’il s’agit d’étudier des phénomènes ou des objets
qui sont difficilement manipulables, observables (de par leur taille, leur dangerosité ou tout
simplement à cause du coût financier des expériences) ou qui ne sont pas encore conçus.

Proposer un modèle et en tirer des résultats exploitables n’est pas une activité anodine. Pour
que les résultats obtenus soient utilisables, il est nécessaire de respecter une certaine
démarche. Depuis plusieurs dizaines d’années, les recherches dans ce domaine ont fait de
l’activité de modélisation et de simulation une science à part entière en proposant théories,
outils et vocabulaire spécifiques.

Ce guide est susceptible d’être modifié et complété régulièrement, au fil des besoins
pouvant apparaitre lors des séances de TP-projets.

1
 Objectif

Tous les problèmes décrits dans ces TP seront traités à l’aide du logiciel COMSOL.
L’objectif de cette séance est une prise en main du logiciel, d’en montrer les possibilités et les
limites dans l’étude numérique des équations aux dérivées partielles. Pour en faire une
utilisation non triviale, il faut comprendre les propriétés des équations aux dérivées partielles
et les propriétés des solutions.
Les problèmes que nous traitons dans ces TP sont simples, leur modélisation ne présente pas
de grandes difficultés. Dans des problèmes plus complexes il faut d’abord construire un
modèle mathématique bien posé puis définir un modèle numérique cohérent et de difficulté
raisonnable avant d’espérer faire une simulation sur un logiciel.
Ce logiciel vous ouvre d’extraordinaires possibilités de faire des expériences, simulées
numériquement, dans tous les domaines de la science : donc explorez, testez, imaginez,
jouez... sans oublier qu’un modèle mathématique n’est pas la réalité

Quelques conseils

 Le logiciel Comsol permet de calculer une approximation par la méthode des éléments
finis de la solution d’un système d’équations aux dérivées partielles très général, en
statique ou en dynamique. Les différents “modèles” ne sont que des particularisations
d’un algorithme général : vocabulaire adapté, choix de valeurs par défaut des
paramètres.
 En principe Comsol permet donc de modéliser toutes les EDP, c’est sa force, mais sa
faiblesse est dans sa généralité : il y a souvent des méthodes plus efficaces,
notamment pour les problèmes hyperboliques pour lesquels la méthode des éléments
finis n’est pas bien adaptée.
 Il est nécessaire d’utiliser l’aide pour comprendre ce que fait réellement une fonction
particulière du logiciel.

2
 Présentation générale du logiciel Comsol

Pour accéder au logiciel, connectez-vous sur une machine, lancez le logiciel en double-
cliquant sur l’icône COMSOL MULTIPHYSIC.

Au démarrage du logiciel la dimension de l’espace dans lequel s’effectuera l’étude est à choisir
(1D, 2D, 2D axisymétrique, 3D,…), puis le navigateur de physique apparait. Pour chacun des
domaines de physique il est précisé quelles sont les variables et quel est le suffixe propre à
cette physique. Ces variables sont toujours suffixées par le nom du modèle auquel elles sont
rattachées. Par exemple Es est la composante en x du champ électrique dans le modèle
« électrostatique ».

Lorsque le ou les domaines physiques ont été choisis, il est question ensuite de spécifier le
domaine d’étude (stationnaire, fréquentiel,...etc.) et de cliquer sur « étude », l’écran général
de Comsol apparait.

D’une manière générale, lors de la création d’un projet il faut parcourir le menu en allant de la
gauche vers la droite. Nous allons donc d’abord construire la géométrie du problème.

On définit ensuite le comportement des domaines en leur associant des propriétés. Ensuite,
on impose des conditions sur certaines frontières, et éventuellement sur certains points. Cette
opération est à refaire pour chaque modèle utilisé dans le projet.

Une fois la physique du problème posée, il faut effectuer le maillage de la géométrie. Il est
possible de faire un raffinement local du maillage (normal, fin, extra fin, extrêmement fin…).Il
ne reste ensuite plus qu’à simuler le projet.

Toutes les inconnues sont alors déterminées. Il est possible de les afficher dans la fenêtre de
post –traitement. L’utilisateur peut choisir quelles variables il désire représenter et sous quelle
forme.

3
 TP 1: Dipôle électrostatique
1-) But
Le but de ce TP est de calculer la distribution du potentiel et du champ électrique d’un dipole
électrostatique.

2-) Rappel théorique

Dipôle électrostatique : L’ensemble de deux charges ponctuelles opposées (+q et –q) à une
distance d.
-q +q
 
N P
-Charge de dipôle : q

-On définit le moment dipolaire électrique par : p  qNP

 Potentiel électrique du dipôle :

On note O le milieu de [NP], a = NP.

On définit l’axe Oz par la droite (NP) orienté par NP .Un point M de l’espace est repéré par
ses coordonnées sphériques r, θ, φ.

Potentiel crée en M :

qa cos  [Link]
V (M )  
4 0 r 2
4 0 r 2 (1) On l’appelle le potentiel du dipôle P en O.

 Champ électrique du dipôle

V 1 V 1 V
E ( M )   grad M V   er  e  e
r r  r sin  

P cos 
V (M ) 
4 0 r 2

3( P.r )r  r 2 P
E(M ) 
4 0 r 5 (2)

4
Cette formule peut ainsi être utilisée dans n’importe quel système de coordonnées.
r  OM ; r  OM

3-) Etapes de résolution dans Comsol Multiphysics

1. Choix du modèle

Pour commencer l’étude, il faut ouvrir le logiciel Comsol multiphysic 5.1 en double-cliquant sur
son icône. Cette action entraine l’ouverture d’une fenêtre appelant à la création d’un nouveau
modèle d’étude. Deux options nous sont proposées, de passer par l’assistant de modélisation
qui nous permettra de construire notre modèle selon nos objectifs ou de prendre un modèle
vierge.

 Faites le choix de « Assistant de modélisation » (1)

 Sélectionnez la dimension de l’espace de travail (2D) (2)

 Observez la fenêtre qui s’affiche, sélectionné dans le navigateur ’’Radio fréquence’’,

le module onde électromagnétique (emw).Puis cliquez sur « ajouter », puis sur
« étude » (3)

5
 Choisir le domaine fréquentiel et démarrer (4)

L’interface de travail du logiciel est maintenant disponible.

6
 2. Géométrie du modèle

La première chose à faire une fois le programme ouvert est de dessiner la géométrie de notre
système.

 On choisit l’unité micromètre au niveau de l’unité de longueur dans la géométrie pour créer
par la suite nos formes géométriques. (5)

 Faire un clic droit sur géométrie et choisir le carré, ensuite, entrer 10um comme longueur
du coté, positionné à (x= -5 ; y= -5) puis cliquez sur construire tous les objets (6)

 Faites de nouveau un clic droit sur géométrie et sélectionner un point et construisez le (7)

7
 3. Introduction des propriétés des domaines

Le but ici est de définir les paramètres physiques qui définissent le domaine : la constante
diélectrique, la conductivité, la permittivité relative etc. On a à définir également les différents
matériaux liés à notre modèle physique étudié.

 Aller sur l’onglet Matériaux, faire un clic droit et choisir « ajouter un nouveau matériau,
cliquez sur ‘’ intégré’’, puis double-cliqué sur « Air » (8)

Le carré construit précédemment est automatiquement défini comme de l’air.

 Cliquez droit sur l’onglet ondes électromagnétiques Domaine fréquentiel (emw),et

choisissez « conditions aux limites de diffusion » (9)

 Allez sur l’interface graphique de dessin, placer le pointeur de la souris et sélectionnez

chacun des côtés du carré puis revenez à nouveau sur l’onglet « ondes
électromagnétique (emw) », cliquez droit, allez sur ‘’points’’ puis choisir ‘’dipôle
électrique ponctuel’’ (11)

 Cliquez ensuite sur le point apparu sur l’interface graphique pour le sélectionner. Au
niveau de la spécification du dipôle prenez ‘’Moment dipolaire’’ et entrez les paramètres
de P (-1 1) (12)

Le point ainsi crée est considéré comme dipôle électrique.

4. Maillage des domaines

Le maillage des domaines est automatique. Il est possible de le raffiner globalement ou de

manière localisée.

8
  Allez sur l’onglet « Maillage », à la taille des éléments choisir ‘’extra fin’’ et cliqué sur
construire toute la séquence. (13)

5. Résolution du problème

Lorsque le problème est complètement défini, il ne reste plus qu’à le résoudre. Il est possible
d’accéder à des menus pour modifier la résolution. De même, lors de la résolution d’un
problème multi- physique, il est possible de ne résoudre que pour une certaine partie des
variables.

 En dessous de l’onglet « Etude », cliquez sur ‘’étape 1 domaine fréquentiel’ ’au niveau
de la valeur de la fréquence entrer : c_const/500[nm] (14)

 Cliquez ensuite sur « calculer » et attendre que le résultat de la simulation apparaisse.

(15)

 Afficher les résultats sous forme de surface de la norme du champ électrique E. (16)

9
  Afficher le résultat pour de [Link] ; Ex^2 ; puis sqrt(Ex^2+Ey^2) ; enfin :
sprt([Link]^2+[Link]^2+[Link]^2). (17)

 Interprétez l’ensemble des résultats ainsi obtenus (18)

---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------

10
TP 2: Potentiel électrique autour de 3 fils de phase (A réaliser par les étudiants)
1-) But
L’objectif de ce TP est de déterminer le champ et le potentiel de trois fils de phase sous
tension et autour de ceux-ci.

2-) Rappel théorique

Considérons les conducteurs en équilibre électrostatique. Un conducteur est en
équilibre électrostatique lorsque la répartition des charges qu’il contient est indépendante du
temps c’est-à-dire qu’aucune charge électrique ne se déplace plus à l’intérieur du conducteur.
Cela signifie que la vitesse moyenne des charges, et en particulier celui des charges libres,
est nulle. Il n’y a donc pas de courant électrique à l’intérieur.

 Champ électrique à l’intérieur d’un conducteur.

A l’équilibre, l’intérieur d’un conducteur sur le plan macroscopique n’est le siège d’aucun
courant électrique. Le champ électrique intérieur est donc nul. On peut donc écrire que :

Eint  0 (1)

 Potentiel électrique d’un conducteur à l’équilibre.

E int   grad .V  0 (2)

On en déduit que le potentiel V est constant dans le volume du conducteur a l’équilibre. On

peut donc écrire :

Vint  cste
(3)

Le conducteur à l’équilibre est donc un volume équipotentiel.

 Champ au voisinage d’un conducteur.

Le champ électrique externe situé au voisinage immédiat d’un conducteur est perpendiculaire
à la surface. Ceci est dû au fait que la surface est une équipotentielle et que les lignes de
champ sont perpendiculaires aux équipotentielles. Ce champ électrique est lié directement à
la densité de charges surfaciques locales.

Puisque le champ électrique interne est nul, si nous appelons E ext le champ électrique externe
immédiat, le théorème de coulomb donne en fonction de la densité superficielle de charges

11
noté σ et de la permittivité du milieu (le vide pour les conducteurs de transport aérien).La valeur

E ext  n
de ce champ est :
0 (4)

Où n est un vecteur unitaire normal au conducteur et dirigé vers l’extérieur.

3-) Etapes de résolution

Air
 Géométrie :

A
Ame
conductrice

B C
Isolant

Figure 2.1 : Conducteurs

Le diamètre des conducteurs est de 20mm

L’épaisseur de l’isolant est de 2mm

 Données :

Permittivité relative de l’aluminium est de Ԑ=1

Permittivité relative de la couche isolante (PVC) Ԑ= 5

Ligne sous tension U= 3KV

 TAF :

1-) En suivant l’organigramme ci-dessous, déterminer à l’aide du logiciel Comsol le potentiel

et le champ électrique autour de la ligne constitué des 3 fils conducteurs.

2-) Afficher les résultats du potentiel électrique des fils, sous forme de surface, puis sous forme
d’isovaleurs.

12
3-) Afficher les résultats de la distribution du champ électrique des 3 fils et vérifier leurs
conformité avec l’équation (1)

4-) Afficher le résultat du calcul du champ suivant la composante x (Ex)

5-) Interprétez l’ensemble des résultats et concluez.

---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------

13
 Début

Choix de la dimension d’étude

2D

Choix de la physique
. AC / DC
…….Electrostatique(es)
………..Stationnaire

Description de la géométrie de la ligne à

. étudier

Choix des matériaux de la ligne à

. étudier

Insertion des paramètres physique :

Conditions aux limites
Potentiel sue les conducteurs

Génération du maillage

Non
Maillage effectué
Oui

Effectuer le calcul

Non
Calcul effectué
Oui
Extraction des résultats

Fin

Figure 2.2: Organigramme de calcul avec Comsol Multiphysics

14
 T P 3 (A réaliser avec l’enseignant)

Au cours de ce TP, nous allons tenter de simuler le comportement d’un

actionneur électrostatique simple. Il s’agit d’une poutre encastrée qui se déforme sous l’action
des forces électrostatique.

Pour cela, nous allons d’abord calculer la distribution des charges et du champ électrique dans
le système puis nous utiliserons un module RDM qui nous donnera la déformation de la poutre
en fonction des forces appliquées sur celle-ci.

1- Choix du modèle physique

Nous allons commencer par l’étude du champ électrostatique.

Lancez donc le module ‘’électrostatique dans le navigateur de physique. Nous allons ici nous
limiter à une étude à 2 dimensions.

 Observez le nom des variables dépendantes et le nom de l’application qui seront

générés

2- Géométrie

La première chose à faire une fois le programme ouvert est de dessiner la

géométrie de notre système.

2mm
1.5mm

1mm

10mm

15mm

 Remarquez que chacun des rectangles que vous dessinez peut être édité. Chaque
figure reçoit un nom par défaut.

15
 3- Introduction des propriétés des domaines
Le but ici est de définir les paramètres physiques qui définissent le domaine :
constante diélectrique, permittivité relative, charges volumique, etc. Les paramètres sont
évidemment liés au modèle physique étudié. Il est également possible d’écarter un des
domaines du modèle étudié.

Une bibliothèque de matériaux est à disposition dans Comsol. Les deux lames sont en
aluminium tandis que le bloc entre les deux sera en verre. Le dernier domaine est constitué
d’air.

Remarquez que si aucune propriété concernant le matériau utilisé n’est pertinente pour le
modèle étudié les valeurs précédentes seront conservées.

On considèrera dans notre exemple qu’il n’y’a pas de charges d’espace.

 Explorez les différents onglets pour voir les différentes options proposées.

4- Introduction des conditions aux limites

Pour résoudre les EDPs, il est nécessaire de définir un certains nombres de conditions
aux limites cohérentes.

A chaque frontière de domaine est associé un segment orienté. Il est possible de donner des
valeurs de conditions aux limites pour chacun de ces segments. Les options seront différentes
suivant qu’il s’agisse ou non d’une frontière intérieure au système.

Attention : Remarquez déjà que comme les frontières sont orientées, les normales aux
frontières ne seront pas nécessairement des normales extérieures.

Les conditions doivent être choisies de manière à ce que la lame supérieure soit portée à un
potentiel de 10V et la lame inférieure soit portée à la masse.

 Choisissez de manière pertinente les différentes conditions aux limites. Vous pourrez
les changer plus tard, de manière à voir leur influence.

5- Maillage des domaines

Le maillage des domaines est automatique. Il est possible de le raffiner globalement ou

de manière localisée.

 Réalisez le maillage automatique du système et affinez globalement ‘’fin’’ une fois

6- Résolution du problème

Lorsque le problème est complètement défini, il ne reste plus qu’à le résoudre.

16
Il est possible d’accéder à des menus pour paramétriser la résolution. De même, lors de la
résolution d’un problème multiphysique, il est possible de ne résoudre que pour une certaine
partie des variables.

 Essayez de résoudre le problème. S’il y’a des problèmes de convergence, il est

possible que le problème soit mal posé (conditions aux limites).

7- Post-traitement : analyse des résultats

Lorsque le programme a convergé vers une solution, il bascule automatiquement vers

la fenêtre de post-traitement, qui permet d’analyser les résultats.

Il est possible de les visualiser sous forme de remplissage de couleur, de lignes de courant,
de champ de flèches, d’isovaleurs, etc…

 Essayez de visualiser les lignes de champ électrique, les équipotentielles, les charges
de surface qui apparaissent sur les frontières. Vous pouvez éventuellement changer
les types de conditions aux limites pour en coir l’influence.

 Cette représentation vous parait-elle réaliste ? Que se passe-t-il si vous décidez de

modifier la taille de la grande boite ?

 Remarquez qu’on a ici accès au nom de certaines variables, ce qui peut constituer un
aide-mémoire utile et facilement accessible.

8- Enrichissement du modèle

Nous allons à présent enrichir le modèle en incluant la déformation de la lame supérieure

sous l’action des forces électrostatiques. Pour cela nous allons utiliser une formulation de
résistance des matériaux. L’intérêt étant ici d’ajouter une nouvelle physique et non de
remplacer le modèle électrostatique.

 Lancez le navigateur de physique depuis le menu multiphysique pour ajouter une

physique. Sélectionnez la physique « Structural mechanics/plane stress/Statis
analysis » et puis cliquez sur ‘’Add’’

La géométrie utilisée sera bien entendu la même que pour l’étape précédente. On peut donc
directement passer à la paramétrisation du domaine.

17
 9- Introduction des propriétés des domaines

Dans la barre de titre de l’application vous voyez l’indication du modèle actif

actuellement (ps pour plane stress,es pour électrostatique)

Lorsque vous êtes dans le mode « plane stress », il faut entrer les paramètres relatifs aux
différents domaines, comme pour le modèle électrostatique. Remarquez que les paramètres
ne sont évidemment plus les mêmes.

Dans notre exemple, nous nous intéresserons aux contraintes et déformations sur les parties
solides. Le domaine contenant l’air ne doit donc plus être inclus dans ce modèle.

 Introduisez les paramètres des matériaux pour les différents domaines

 Comment feriez-vous pour tenir compte du poids propre des objets ?

10- Introduction aux conditions aux limites

 Regardez les différents types de conditions aux limites que vous pouvez introduire

 Le côté gauche et le dessous de la lame du bas sont supposés être fixes. La lame du
dessus est soumise aux forces électrostatiques. Introduisez les conditions aux limites
pertinentes pour décrire cette situation.

11- Analyse des résultats

 Résolvez le modèle multiphysique et regardez les résultats

 Affichez les déformations de la surface et des frontières

 Affichez les contraintes dans les matériaux

 Que se passe-t-il si vous augmentez ou diminuez la raideur du bloc entre les deux
lames ?

 Affichez le champ électrostatique et les charges. Cette répartirions est-elle correcte ?

Pourquoi ?

18
 TP 4 : Condensateur Plan (A réaliser par l’étudiant)
1-) But
Calcul de la répartition du champ électrostatique du condensateur plan

2-) Rappel théorique

On considère un condensateur plan constitué de deux surfaces conductrices (S) planes et
parallèles dont les dimensions sont grandes par rapport à la distance e qui les sépare.

On considère donc que les deux surfaces conductrices (S) sont deux plans infinis parallèles
portant les charges opposées : +Q = σ.S et -Q = - σ.S

Le champ créé par un plan infini uniformément chargé en surface (σ) est :


E n
2 0

Pour un point M quelconque placé entre les armatures du condensateur, on obtient d'après le
théorème de superposition :

 
E (M )  E1  E2  ex  (e x )
2 0 2 0


E (M )  ex
0
19
 3-) Etapes de résolution

1. Choix du modèle

Nous allons faire le choix d’une étude 2D, et lancez dans le navigateur de modèle (AC/DC) le
module électrostatique (es).L’étude se fera en domaine stationnaire.

2. Géométrie

Dessinez la géométrie du système (condensateur) à partir des données spécifiées suivants

r3

r1

r2
c1

c2

Figure 3.1 : Modèle géométrique du condensateur

Rectangle 1 Rectangle 2 Rectangle 3

H = 1m H = 1m H = 0.9m

L = 1m L = 1m L = 0.14m

Centre (0.1 0) Centre (-0.1 0) Centre (0 0)

Cercle 1 Cercle 2

R1= 0.03m R2= 0.03m

O1 (-0.1 -0.5) O2 (0.1 -0.5)

20
 3. Introduction des propriétés des domaines

Matériaux : Entrez les matériaux des différents composants

Air : pour la zone d’étude (c3)

Glass (quartz) : pour le diélectrique entre les deux armatures (r3)

Aluminium 3003-H18: Pour les deux armatures du condensateur.

Pour résoudre les EDP, il est nécessaire de définir un certain nombre de conditions aux limites
cohérentes. Attribuez les différents potentiels électriques aux deux armatures du condensateur

V1 = 0.5V ; V2 = -0.5V

21
4. Maillage
Choisissez le maillage extrêmement fin et le construire
5. Résolution
Calculez et attendez que le programme ait convergé vers une solution.

1- Afficher ensuite le résultat sous forme d’isovaleurs

2- Afficher les lignes de courants

3- Calculez la norme du déplacement électrique

4- Interprétez l’ensemble des résultats, montrez l’influence des effets de bord et concluez

----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------

22
 TP 5
1-) But
L’objectif est de calculer la répartition du champ et du potentiel le long de l’isolateur composite.

2-) Rappel théorique

La tenue diélectrique des matériels d'isolation est tributaire des répartitions du champ
et du potentiel à leurs surfaces. Nous allons considérer les équations et les conditions qui
gouvernent ces répartitions. Nous prenons comme hypothèse que les isolants qui entrent en
jeu ici (air, ailettes, ...) possèdent une permittivité de valeur constante, une résistivité
volumique nulle et une densité volumique de charge libre nulle.

Les équations qui gouvernent la répartition du potentiel et du champ électrique, dans un

milieu donné, dérivent des équations de Maxwell. Celles-ci sont formées par quatre équations
aux dérivées partielles qui lient les phénomènes magnétiques caractérisés par le champ
magnétique H et l'induction magnétique B aux phénomènes électriques caractérisés par le
champ électrique E et l'induction électrique D. Ces quatre équations sont :

B
rot E  
t (1)

divD  v (2)

D
rot H  j 
t (3)

divB  0 (4)
Avec ρv la densité volumique de charges.
Outre ces quatre équations, il y a les relations constitutives qui lient D à E, J à E et B à H:

D  E (5)

j  E (6)

B  H (7)
Où ε, µ, et σ représentent respectivement la permittivité électrique, la susceptibilité magné-
tique, et la conductivité électrique du milieu.
Lorsque l'on applique une différence de potentiel alternative, la dérivation des grandeurs
électriques par rapport au temps revient à les multiplier par la quantité jw où j est l'unité ima-
ginaire et w la pulsation de la source. En tenant compte des relations constitutives, les équa-
tions (1) et (3) deviennent :

23
 rot E   jwB   jw. H (8)

 
div  E  v
(9)

rot H  E  jwE (10)

Comme le champ magnétique est faible, dans les équipements H. T considérés, le second
membre de l'équation (8) peut être négligé, ce qui donne :

rot E  0 (11)
On peut donc découpler les équations (9) et (11) qui gouvernent les grandeurs électriques
de celles qui gèrent les grandeurs magnétiques. L'équation (9) n'est alors rien d'autre que la
forme locale du théorème de Gauss qui traduit la conservation de la charge. L'équation (11)

permet de dire que le champ E dérive d'un potentiel  avec :

E   grad . (12)
L'introduction de la relation (12) dans l'équation (9) donne :

div  .grad .   v
(13)
On obtient ainsi l'équation qui gouverne la répartition du potentiel dans un milieu de permittivité
absolue ε avec une densité de charge volumique de charge ρv.
Pour des matériaux isolants, généralement utilisés dans les appareillages H. T, la den-
sité volumique de charge ρv est nulle. Avec ces hypothèses, on obtient l'équation de Laplace
(14) qui gouverne la répartition du potentiel dans les matériaux isolants.

 2  0

(14)

3-) Etapes de résolution

1. Choix du modèle

Prenez un modèle vierge, une fois le programme lancé créer un composant 2D à partir duquel
vous construirez plus tard un modèle 3D de l’isolateur. A ce moment-là vous ajouterez un
composant 3D avec pour physique le module électrostatique.

2. Géométrie

La procédure pour la construction du modèle géométrique 3D de l’isolateur est dans un

premier temps de dessinez la géométrie 2D (voir fig. 1) et ensuite d’en éditer un modèle 3D
(voir fig. 2)
24
 0.02m

0.2m

0.10m

3. Introduction des propriétés des domaines

Haute tension (HT) : 10KV

Air permittivité relative : Ԑ=1

Silicone permittivité : Ԑ= 11.7

4. Maillage
Choisissez un maillage ‘plus fin’ et construisez la séquence.

25
 5. Résolution du problème
Calculez

 Afficher ensuite le résultat sous forme de surface

 Commentez les résultats obtenus.
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
-----------------------------------------------------

26
Conclusion

L’approche de modélisation et simulation par le biais du logiciel Comsol a pour

but d’accroitre la connaissance que possède un observateur sur un phénomène donné.
L’intérêt de cette démarche réside dans le fait de simplifier l’objet d´étude en créant un modèle.
Ce modèle est beaucoup plus simple à manipuler, à observer que le phénomène réel
modélisé. Les conditions expérimentales sont toujours maitrisées et une même expérience
peut être rejouée autant de fois que nécessaire. L’enjeu de cette méthode consiste à créer un
modèle suffisamment représentatif du phénomène étudié pour que les réponses apportées
par le modèle puissent s’appliquer au phénomène réel étudié.

La méthodologie associée à la modélisation et à la simulation est relativement bien structurée.

Cependant, il est important de noter qu’elle est liée au cadre de pensée du modélisateur et
aux questions posées sur le phénomène réel étudié. Ainsi, un même phénomène peut être
représenté de manières différentes et plusieurs modèles et simulateurs hétérogènes peuvent
coexister. Un même phénomène peut être représenté par plusieurs modèles décrits dans des
formalismes différents et simulés grâce à des politiques d’exécutions différentes. De même le
paradigme de modélisation va changer en passant de l’étude de systèmes non-complexes (qui
présentent des mécanismes de causalité linéaires) aux systèmes complexes.

27