Memoire Magister 4

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
MINISTERE DE LENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

UNIVERSITE MOULOUD MAMMERI DE TIZI-OUZOU
FACULTE DE GENIE ELECTRIQUE ET DE LINFORMATIQUE

DEPARTEMENT DELECTROTECHNIQUE
MEMOIRE DE MAGISTER
EN ELECTROTECHNIQUE
Option : Entrainements Electriques
Prsent par
Mr. Madjid MEZIANI

Ingnieur dtat en lectrotechnique
Thme
Contribution la modlisation
analytico-numrique des transformateurs
de puissance
Soutenu publiquement le 18/12/2011 devant le jury dexamen :
Mr. Hassane MOHELLEBI
Professeur
Mr. Mhemed RACHEK
Maitre de Confrences A UMMTO
Rapporteur
Mr. Nacer Eddine BENAMROUCHE
Professeur
UMMTO
Examinateur
Mr. Salah HADDAD
Professeur
UMMTO
Examinateur
Mr. Tahar OTMANE CHERIF
Maitre de confrences A
UMMTO
Examinateur
UMMTO
Prsident
Remerciements
Jadresse des remerciements tout particuliers mon promoteur Monsieur Mhemed

RACHEK maitre de confrence A luniversit de Tizi-Ouzou, qui a rendu possible ce
travail de mmoire en me faisant profiter de ses vastes connaissances et de son exprience.
Je remercie profondment Monsieur Hassane MOHELLEBI professeur au
dpartement de llectrotechnique de lUMMTO de mavoir fait lhonneur daccepter de
prsider le jury de soutenance.
Cest pour moi un grand honneur de voire juger ce travail par monsieur Nacer
Eddine BENAMROUCHE professeur au dpartement de llectrotechnique de lUMMTO, je
voudrai lui tmoigner aussi, ma profonde gratitude et ma sincre reconnaissance.
Que Monsieur Salah HADDAD professeur au dpartement de llectrotechnique de
lUMMTO, soit vivement remerci pour lintrt quil accord au sujet et pour nous avoir
fait lhonneur daccepter lvaluation de ce travail.
Mes sincres remerciements vont galement Monsieur Tahar OTMANE CHERIF
maitre de confrences A luniversit de Tizi-Ouzou pour avoir accept de participer
dans le jury de soutenance.
Je tiens vivement remerci Monsieur Arezki DICHE enseignant lUMMTO, pour
ses conseils judicieux et sa sympathie.
Sommaire
Introduction gnrale
Chapitre I : Etat de lart sur la modlisation des transformateurs

I.1 Introduction
I.2 Principaux lments des transformateurs...
I.3 Principe de fonctionnement
I.4 Diffrents types de transformateurs..
I.4.1 Transformateur colonnes...
I.4.1.1 Les enroulements.....
I.4.1.2 Le circuit magntique..
I.4.1.3 Limitations des transformateurs colonnes....
I.4.2 Transformateur cuirass..............
I.4.2.1 Les enroulements
I.4.2.2 Les crans lectrostatiques.
I.4.2.3 Avantages et inconvnients du transformateur cuirass.....
I.4.3 Transformateur ordinaire de phase...
I.4.4 Transformateur immerg dans l'huile et transformateur sec...
I.4.4.1 Transformateur lhuile.....
I.4.4.2 Transformateur sec.....
10
I.5 Transformateurs triphass..
10
I.6 Groupe de couplage................ ...
10
I.6.1 Couplages normaliss...
11
I.6.1.1 Couplage toile toile..................
11
I.6.1.2 Couplage triangle toile
11
I.6.1.3 Couplage toile zigzag.
12
I.6-2 Choix du couplage dans les transformateurs usuels.
12
I.7 Contraintes subies par le transformateur
13
I.8 Rgimes de fonctionnement...
13
I.8.1 Rgime de surcharge
13
I.8.2 Rgime dsquilibr.....
14
I.8.3 Rgime transitoire
14
I.9 Protection des transformateurs contre les surtensions
14
I.10 Diffrentes approches de modlisation des transformateurs
14
I.11 Schmas quivalent en tant que modle...
15
I.11.1 Schma classique
15
I.11.2 Phnomnes lectromagntiques.
16
I.11.2.1 Couplage parfait (transformateur idal)
17
I.11.2.2 Couplage imparfait.................................................
17
I.11.2.3 Pertes cuivre.
17
I.11.2.4 Calcul des rsistances ..
18
I.11.2.5 grandeurs magntiques.
18
I.11.2.6 Energie magntique emmagasine...
19
I.11.2.7 Pertes fer...
19
I.11.3 Phnomnes lectrostatiques
20
I.11.3.1 Pertes dilectriques...
20
I.11.3.2 Calcul des capacits..
20
I.11.3.3 Capacit entre enroulements
21
I.11.3.4 Capacit entre conducteurs
21
I.11.3.5 Capacit entre phases...
22
I.11.3.6 Capacit entre enroulement et cuve..
22
I.12 Modlisation par les mthodes semi-analytiques.
23
I.13 Modle dtaill (fin) bas sur le calcul des inductances et des rsistances.
23
I.13.1 Calcul des selfs inductances..
24
I.13.2 Calcul des inductances mutuelles..
24
I.13.3 Calcul des rsistances
25
I.14 Modlisation par la mthode des lments finis (MEF) 25

I.15 Conclusion. 25
Chapitre II : Modlisation lectromagntique des transformateurs
II.1 Introduction...
26
II.2 Equations de Maxwell..
26
II.3 Lois de comportement des matriaux...
27
II.4 La loi dOhm.....
27
II.5 Relations de passage.....
28
II.6 Equations de Maxwell simplifies.
29
II.7 Modles lectromagntiques
29
II.7.1 Modle magntostatique......
29
II.7.2 Modle magntodynamique
30
II.8 Le choix de la variable dtat...
30
II.9 Formulation en potentiel vecteur magntique...
31
II.10 Conditions de Jauge
31
II.11 Formulation bidimensionnelle en potentiel vecteur magntique..
32
II.11.1 Formulation 2D cartsien.....
32
II.11.2 Formulation 2D axisymtrique....
33
II.12 Les conditions aux limites..
33
II.13 Mthodes de rsolution..
34
II.13.1 Mthodes numriques..
34
II.13.1.1 Principe de la mthode des lments finis.....
34
II.13.2 Mthodes semi-analytiques..
35
II.13.2.1 Principe de la mthode des circuits coupls...
35
II.13.2.2 Interaction entre les lments de circuits
36
II.13.2.3 Alimentation sinusodale.
38
II.14 Conclusion
38
Chapitre III : Modlisation des transformateurs par une mthodologie mixte

Circuits Electriques Magntiquement Coupls (CEMC) et Elments Finis(EF)
III.1 Introduction.
39
III.2 Modle gomtrique
39
III.3 Discrtisation des enroulements..
40
III.4 Formulation en potentiel vecteur.
41
III.4.1 Equations caractristiques en axisymtrique
41
III.5 Expression de potentiel vecteur magntique
44
III.6 Prise en compte du circuit magntique.........................
49
III.6.1 Modlisation du noyau par des densits superficielles de courant
50
III.6.2 Expression intgrale caractristique du noyau.
51
III.7 Calcul analytique des paramtres inductifs du modle...
54
III.7.1 Inductance mutuelle..
55
III.7.2 Inductance propre...
56
III.8 Ncessit du couplage avec la mthode des lments finis..
57
III.8.1 Technique mixte circuits coupls-lments finis.... 58

III.8.2 Formulation lments finis des quations lectromagntiques du noyau... 58
III.8.3 Mthodes de rsolution des systmes algbriques
61
III.9 Conclusion...
64
Chapitre IV: Applications est validations

IV.1 Introduction...
65
IV.2 Description du transformateur..
65
IV.3 Caractristiques du dispositif
66
IV.4 Dimensions gomtriques du transformateur...
66
IV.5 Caractristiques physiques
67
IV.6 Discrtisation de la structure.
68
IV.7 Calcul des pertes.... 71

IV.7.1 Le processus de calcul des pertes dans les enroulements.. 71
IV.7.2 Pertes en court-circuit..........................................
71
IV.7.2.1 Calcul des rsistances et selfs....
72
IV.7.2.2 Calcul des mutuelles et selfs de chaque conducteur 73

IV.7.2.3 Calcul des inductions et 73
IV.7.3 Calcul des pertes vide...
74
IV.8 Interprtation des rsultats..
81
IV.9 Description du second transformateur
83
IV.10 Caractristiques du dispositif.
83
IV.11 Dimensions gomtriques du second transformateur.
83
IV.12 Essai en court-circuit.
85
IV.13 Essai vide...
86
IV.14 Conclusion
86
Conclusion gnrale....
87
Bibliographie
Introduction Gnrale
La modlisation des transformateurs de puissance

beaucoup dattention durant des annes, pour
traditionnellement occup
leur importances
dans les systmes de
puissance, le grand intrt port ltude de ces dernier on permit de raliser un

dveloppement remarquable dans le domaine de la conversion, du transport et de la
distribution de lnergie lectrique. La comprhension des phnomnes lectromagntiques
qui rgissent leur fonctionnement et la dtermination des diffrentes pertes et les paramtres
lectriques constitue une tape dterminante pour une meilleure exploitation de ces
dispositifs. Ce qui permet galement de prdire sont comportement aux diffrentes contraintes
pour viter sa dtrioration, et largir sa dure de vie en vue de minimiser son cot
dexploitation. Il existe deux types de modlisations : la modlisation par calcul de champs
dite numrique ou semi-analytique et la modlisation par circuit lectrique quivalent dite
analytique.
Vu la complexit des caractristiques gomtriques et physiques des systmes
lectromagntiques rels, les mthodes disponibles sont alors des mthodes numriques ou
semi-analytiques qui font appel des techniques de discrtisations. Les outils numriques de
modlisation sont ncessaires pour traiter des problmes lis au fonctionnement du
transformateur en exploitation et aussi leur conception, cela offre une modlisation fine des
phnomnes physiques en permettant daccder aux variables locales et aux variations spatiotemporelles des diffrentes grandeurs. Plusieurs mthodes ont ts dveloppes, chaque
mthode ses avantages et ses inconvnients, ses applications et ses limites. On peut citer la
mthode des volumes finis, mthode des diffrences finis et mthode des lments finis. Le
choix se fait selon plusieurs contraintes physiques et selon la complexit des gomtries des
dispositifs lectrotechniques, les proprits physiques des matriaux, et le cout en capacit
mmoire et en temps de rsolution.
Le but de notre travail est de mettre en uvre une mthode de modlisation, pour
retranscrire le comportement lectromagntique du transformateur de puissance dans les
diffrents tats oprationnels aussi bien en rgime permanent quen rgime transitoire.
On a opt pour la mthode des Circuits Electriques Magntiquement Coupls (CEMC)

qui se base sur des techniques de discrtisations. Une mthode qui offre beaucoup de
possibilits de modlisations locale des transformateurs et permet aussi lvaluation des
paramtres lectriques dans une large gamme de frquences. Cette dernire trouve ces limites
en regardant la gomtrie rel du transformateur et aux dtails (non linarit, hystrsis,
anisotropie de la permabilit magntique et conductivit lectrique du circuit magntique).
La prsence de circuit magntique complique lapproche semi-analytique, la mthode des
lments finis est une excellente alternative pour des gomtries complexes et caractres
physiques non linaires. Mais elle peut conduire des difficults pour la prise en compte de
grand nombre de conducteurs, et plus particulirement en prsence de conducteurs filaires.
Pour cette raison on adopt pour une mthodologie hybride pour la modlisation
lectromagntique qui sappuie sur des formulations analytiques et numriques bases
respectivement sur la thorie des circuits coupls et la mthode des lments finis.
Ce travail ce dcompose en quatre chapitres :
Le premier chapitre expose les diffrents types de transformateurs de puissances, leur
description et constitution ainsi que les mthodes de modlisation. Nous aborderons le
problme de modlisation en utilisant deux mthodes : la premire est base sur le schma
lectrique quivalent global et la deuxime mthode est une approche locale base sur un
calcul numrique ou bien semi-analytique effectu respectivement par la mthode des
lments finis (EF) et la mthode des circuits lectriques magntiquement coupls (CEMC).
Le deuxime chapitre a pour objet la description des phnomnes lectromagntiques
dans le cadre des hypothses de l'lectrotechnique. Le modle mathmatique, que nous allons
laborer comprendra les quations de Maxwell, les lois de comportement des matriaux et les
conditions aux limites. Puis, nous donnons les formulations, magntostatique et
magntodynamique en introduisant le potentiel vecteur magntique.
Le troisime chapitre est consacr lapplication de la mthode des circuits lectriques
magntiquement coupls sur un transformateur de puissance. La considration de contraintes
de conservation en courant et tension permet l'assemblage du systme matriciel
lectromagntique complet ncessaire au calcul de la distribution de densits de courant et
donc des pertes Joule dans les conducteurs.
La problmatique dtude ncessite une mthode mixte ou bien
un couplage
analytico-numrique, cela nous permettra de tirer profit de deux mthodes, afin de pouvoir
retranscrire efficacement le comportement de cet quipement. Le circuit magntique rel sera
pris en charge travers le modle magntique formul numriquement par la mthode des
lments finis.
Le quatrime chapitre est consacr la validation ainsi que lexploitation des rsultats
obtenus avec le modle dvelopp. Ces derniers seront confronts aux mesures
exprimentales donnes par lentreprise Electro-Industries dAzzazga.
Chapitre I
tat de lart sur la modlisation
des transformateurs
Chapitre I
tat de lart sur la modlisation des transformateurs
I.1 Introduction
Les
transformateurs
sont
des
machines
statiques
induction
lectromagntique
gnralement utiliss pour modifier les caractristiques de lnergie lectrique alternative afin
de la rendre aussi commodes que possible tous les stades, de la production, du transport
de la distribution et de lutilisation. Cet appareil est dun emploi absolument universel, en
effet il transforme les signaux des sources de tensions et de courants sinusodaux en signaux
de mme frquence mais de valeurs efficaces gnralement diffrentes.
La comprhension des phnomnes lectromagntiques qui rgissent le fonctionnement
des transformateurs et la dtermination de ses pertes, constitue une tape dterminante pour
lexploitation optimale de ces appareils et galement pour leur conception. La modlisation
des transformateurs est une tche trs difficile due l'introduction des caractristiques du
noyau magntique tel que la saturation, lhystrsis et les pertes par courants de Foucault,
aussi bien que de sa topologie et la configuration de ses enroulements. Ainsi, une varit de
mthodes pour sa modlisation et analyse de sont comportement pour diffrents tats
oprationnels ont t dveloppes. On peut distinguer deux approches de modlisation :
Approche locale oriente vers le calcul numrique ou bien semi-analytique.
Approche globale oriente vers le calcul analytique bas sur un schma quivalent
lectrique ou magntique (rseau de rluctances).
I.2 Principaux lments des transformateurs

D'une manire gnrale, un transformateur est constitu d'un circuit magntique feuillet
et d'un ensemble de bobines spares par des crans lectrostatiques qui entourent des noyaux
magntiques. Chaque bobine formant le milieu conducteur est organise en paquets de spires
et
chaque
spire
tant
constitue
de
brins
lmentaires.
Les
diffrents
types
de transformateurs se distinguent suivant la disposition gomtrique de leurs constituants,

de la forme de leur circuit magntique et du type de refroidissement.
I.3 Principe de fonctionnement
En vertu de la loi de Faraday, lorsquun flux dinduction magntique variable circule dans
le circuit magntique, il induit dans chacun des enroulements une force lectromotrice
proportionnelle dans le temps aux taux de changement ( / ) et au nombre de spire que
comporte cet enroulement.
4
Chapitre I
Lorsque le primaire est aliment par une source alternative, il circule dans le circuit
magntique un flux galement alternatif dont lamplitude dpend du nombre de spires du
primaire et de la tension applique. Ce flux induit dans lenroulement secondaire une tension
proportionnelle au nombre de spires du secondaire. La fermeture du secondaire sur une charge
provoque la circulation du courant secondaire.
Le noyau magntique fournit un chemin de canalisation de flux magntique tel que montr
par la figure (I-1)
I.4 Diffrents types de transformateurs
Il existe diffrents types de transformateurs, mais nous nous intresserons ltude

des transformateurs de puissances de rseaux lectriques car ils prsentent le plus grand
intrt dans le dveloppement de linterconnexion des rseaux.
du point de vue construction, deux principales technologies sont distingues, savoir les
transformateurs colonnes et les transformateurs cuirasss.
I.4.1 Transformateur colonnes
Le transformateur colonnes est constitu de deux enroulements concentriques par phase.

Ces enroulements sont monts sur un noyau ferromagntique qui se referme ses extrmits
via des culasses gnralement de sections circulaires afin d'assurer une bonne canalisation du
flux magntique. Dans cette technologie, ce sont les enroulements qui entourent le circuit
magntique (figure. I.2).
Chapitre I

Les enroulements sont constitus de spires dont le nombre est diffrent pour les enroulements
de haute et de basse tension respectivement primaire/secondaire. Les spires sont elles mmes
subdivises en plusieurs brins mis en parallle. Les conducteurs formant les enroulements
sont en cuivre ou en aluminium pour les puissances plus importantes. Ces matriaux sont ainsi
prfrs d'autres pour leurs bonnes conductivits thermiques. Selon la gamme de puissance
des transformateurs, les conducteurs peuvent tres sous forme de fils massifs, de section
circulaire ou carre, de type mplat, ou encore lamins en fines feuilles [Lef-2006].
I.4.1.2 Le Circuit magntique
Le noyau est compos dun empilage de tles ferromagntiques cristaux orients, isoles
lectriquement entres elles afin de rduire les pertes par courants de Foucault. Les tles de
circuit magntique de type (Fe-Si), grain orient lamin froid, ont une paisseur qui varie
de 0.15mm 0.3mm [Gue-1994]. Des joints enchevtrs et orthogonaux effectuent la liaison
entre les colonnes et les culasses sur les transformateurs de petite puissance. Ces joints qui
constituent une succession de stratifications jouent un rle important dans la performance du
noyau afin d'obtenir une stabilit mcanique plus leve, et diminuent galement le bruit des
vibrations pendant le fonctionnement du transformateur. [Nak-1982] [Oli-2003]. L'espace
dlimit par deux colonnes successives et les culasses est dnomm "fentre magntique" et
accueille donc les bobinages [Lef-2006].
6
Chapitre I
I.4.1.3 Limitations des transformateurs colonnes

Le circuit magntique des transformateurs colonnes est plus grand en volume que celui d'un
type cuirass. Par consquence, le nombre de spires et le rapport volumique entre les
matriaux conducteurs et ferromagntiques sont plus importants. Par ailleurs ce type de
construction qui fait ses preuves pour des puissances relativement faibles (environ jusqu
30 KVA), prsente quelque difficults pour des fortes puissances (plus de 100 KVA) et des
tensions plus leves. Pour des tensions suprieures 220 KV, certains constructeurs ont
abandonn la construction du transformateur colonne au profit de la structure cuirasse.
I.4.2 Transformateur cuirass
Dans cette technologie, le circuit magntique entoure les enroulements forms de bobines
rectangulaires axe horizontal. Le circuit magntique, de section rectangulaire est constitu
de tles poses plat. La cuve assure le calage du circuit magntique et des enroulements.
Ces transformateurs sont utiliss au sein des rseaux de transport et de rpartitions ou les
surtensions transitoires sont frquentes. Dans cet environnement, ils doivent se prmunir des
effets nfastes de ces surtension sur les enroulements. Pour cela des crans sont utiliss afin
de rduire les contraintes lies aux champs lectriques dans les bobinages.
Chapitre I

Dans ce type de construction la bobine qui est en fait une grande spirale rectangulaire
trs plate (galette), contient un certain nombre de spires, chaque spire tant elle mme
constitue dune ou de plusieurs couches de conducteurs. Les galettes sont alternes entre
la haute et la basse tension afin de diminuer les fuites magntiques et le gradient de tension.
En outre, cette topologie offre l'avantage d'engendrer des forces en opposition entre chaque
galette lors de court-circuit.
Galette HT
Conducteurs
Galette BT
Fig. I.4 : Transformateur triphas type cuirass [Kos-1979] [Lef-2006].
I.4.2.2 Les crans lectrostatiques

Des feuilles mtalliques haute rsistivit et de faible paisseur sont employes pour rduire
les contraintes entre les spires. Celles-ci sont insres entre les feuilles isolantes des bobines
haute tension. L'effet capacitif, uniformment distribu de cette manire permet alors de
mieux rpartir la tension sur toute la longueur de la bobine lors d'une onde de choc.
I.4.2.3 Avantages et inconvnients du transformateur cuirass
Le transformateur cuirass possde un circuit magntique plus court, ce qui permet davoir un
courant vide relatif plus faible et ses enroulements sont plus simples car le nombre de spires
est moins grand vu que la section du noyau dans un transformateur cuirass peut tre plus
grande que dans le transformateur colonnes. Mais il prsente aussi certains inconvnients :
ses enroulements sont moins accessibles lagent refroidissant, lexamen et la rparation sont
plus difficiles et demande davantage de matriaux isolants pour la haute tension.
Chapitre I
I.4.3 Transformateur ordinaire de phase

Du point de vue de leurs applications, les transformateurs ordinaires de phase (monophas,
biphas, et triphas), groups dans trois ou cinq colonnes se classent en trois catgories :
Les transformateurs pour les grands rseaux et les grandes centrales, leur puissance varie
de 100 400 MVA.
Les transformateurs pour rseaux de rpartition qui alimentent les lignes moyenne
tension, leur puissance varie de 5 30 MVA.
Les transformateurs de distribution destins lalimentation des utilisateurs de
lnergie clectique en basse tension (380 ou 220V), leur puissance varie de 5 1000 KVA.
I.4.4 Transformateur immerg dans l'huile et transformateur sec
Selon le
type de refroidissement on distingue : les transformateurs lhuile et les
transformateurs secs.
I.4.4.1 Transformateur lhuile
Pour prvenir laction nfaste de lair sur lisolation des bobines et amliorer
le refroidissement du transformateur, on place le noyau magntique avec les enroulements
dans une cuve remplie dhuile minrale. Malgr ces proprits avantageuses, lhuile de
transformateurs a deux dfauts principaux : elle est inflammable et sa vapeur forme avec lair
dans certaines conditions un mlange explosif.
En plus de son rle de rservoir au liquide dilectrique, la cuve assure le maintien mcanique
du circuit magntique et des enroulements. La construction de la cuve gnralement de forme
ovale est lie au calcul thermique du transformateur. Le refroidissement du transformateur est
dautant plus difficile raliser que la puissance du transformateur est grande. La cuve est
tapisse de shunts magntiques, dont on distingue deux types :
Les shunts magntiques forms dun empilement de tles magntiques semblable au

noyau et qui canalise le flux de fuite.
les shunts amagntiques plus conomiques, constitus de plaque de cuivre ou aluminium

ayant pour rle de repousser le flux de fuite.
Chapitre I
I.4.4.2 Transformateur sec
Les transformateurs immergs dans lhuile liquide base minrale ou de silicone sont plus
rpondus pour les plus fortes puissances et les niveaux de tension levs, mais prsentent
des risques de fuite, dincendie et la pollution de lenvironnement. Cela permis de laisser
la place aux transformateurs avec technologie sche, avec des enroulements enrobs
(imprgns). Le systme disolation lectrique est remplac par une rsine (poxyde) et lair.
Des rsines ont t dveloppes pour rsister aux tensions lectriques et aux contraintes
mcaniques, thermiques qui apparaissent dans un transformateur en service. Laspect
favorable est son comportement non inflammable et lger. Ce genre de transformateur
est le plus appropris pour la distribution de l'lectricit en degr lev de sret [Esl-2010].
I.5 Transformateurs triphass
Les transformateurs triphass sont prsents diffrents endroits dans les rseaux lectriques
pour adapter les valeurs efficaces des tensions aux niveaux souhaitables. De faon simplifie,
lnergie lectrique est produite dans les centrales sous des tensions moyennes. Le transport
longue distance exige des hautes tensions afin de limiter les pertes par effet Joule et rduire le
dimensionnement des conducteurs, et lutilisation demande des tensions basses ou moyennes.
Des transformateurs lvateurs sont ncessaires au dpart, et des transformateurs abaisseurs
sont indispensables larrive cot consommateurs.
I.6 Groupe de couplage

Le groupe de couplage dsigne lassociation des couplages des enroulements choisis pour la
haute tension et pour la basse tension. Lenroulement haute tension dun transformateur
triphas de nombre de spires
peut tre connect en toile (symbole Y) ou en triangle
(symbole D). Lenroulement basse tension avec un nombre de spires
peut tre connect en
toile (symbole y), en triangle (symbole d) ou en zigzag (symbole z). Les systmes des
tensions du primaire et du secondaire dun transformateur triphas sont en gnral dphass
avec un angle . Comme cette angle est multiple de 6, on dfinit lindice horaire par :
10
Chapitre I
I.6.1 Couplages normaliss

Trois couplages sont particulirement utiliss en pratique : ce sont les couplages
normalissYy0,Dy11,Yz11.
I.6.1.1 Couplage toile toile
Pour le couplage toile toile, les tensions aux bornes des bobines du primaire est du
secondaire de la mme colonne sont des tensions simples (Fig. I.5-a). Ce qui permet de
/
construire le diagramme vectoriel (Fig. I.5-b). le rapport de transformation est

Le dphasage
de basse tension par rapport la haute tension est nul, ce qui donne un indice
horaire I=0.
Fig. I.5-b : Diagramme vectoriel
Fig. I.5-a: Couplage Yy 0
I.6.1.2 Couplage triangle-toile

Pour le couplage triangle-toile la tension aux bornes dune phase du primaire
est une tension compose, alors que la tension aux bornes de la phase correspondante
du secondaire est une tension simple (Fig. I.6-a). Le rapport de transformation
est
est de
3.
. Le dphasage
entre les tensions du primaire et le secondaire
/6 , ce qui veut dire lindice horaire est I=11. (Fig. I.6 -b)
11
Chapitre I
/6
Fig. I.6-a : Couplage Dy 11

I.6.1.3 Couplage toile zigzag
Pour le couplage toile zigzag (Fig. I.7-a), la tension aux borne dune phase du primaire est
une tension simple alors quune tension aux bornes dune phase du secondaire est la somme
des deux tensions aux bornes de demi bobines ce qui permet de construire le diagramme
vectoriel (Fig. I.7-b). Le rapport de transformation est
Le dphasage
entre les tension est de
3.
/6 , en choisissant la premire dtermination
position positif de langle, ce qui veut dire lindice horaire est I=11.
/2
/2
Fig. I.7-a : Couplage Yz 11
I.6-2 Choix du couplage dans les transformateurs usuels.

La prsence du conducteur neutre dans une distribution basse tension permet de disposer
des deux types de tension ; la tension simple pour lusage domestique usuel, et le systme
triphas de tensions composes pour lusage artisanal ou industriel. De plus il est intressant,
cot haute tension de disposer dun couplage ayant un neutre et de mettre ce neutre, ainsi que
le circuit magntique et les parties mtalliques (cuve) du transformateur au potentiel de la
terre. Ceci permet de rduire les distances disolement des bobines hautes tensions[May2006].
12
Chapitre I
Ce quil faut viter cest davoir le mme couplage au primaire et au secondaire dun
transformateur, ceci vite de transmettre intgralement le dsquilibre ventuel des grandeurs
(courants, tensions) dun cot du transformateur lautre. On voit donc apparatre lintrt
dun couplage Yz offrant un neutre des deux cots et des couplages diffrents.
I.7 Contraintes subies par le transformateur
Les transformateurs de puissance sont soumis plusieurs contraintes qui sont les effets
dactions physiques ou chimiques. Les contraintes normales sont celles qui entrainent une
usure et un vieillissement lent du matriel. Les contraintes anormales tant celles qui
provoquent un vieillissement prmatur ou destruction immdiate du matriel.
Les surtensions de foudre, de manuvre ou de court-circuit peuvent provoquer des dfauts
mcaniques comme les vibrations en rgime permanent ou les efforts lectrodynamiques
sur les bobinages.
Les pertes qui se dgagent dans le noyau et les enroulements du transformateur lors de son
fonctionnement se manifestent par de la chaleur qui induit lchauffement des parties
constitutives de lappareil. Son action se traduit
par une destruction lente des isolants
ou la diminution de leur rigidit dilectrique pouvant par consquent amplifier le risque

davarie du matriel.
I.8 Rgimes de fonctionnement
Les rgimes de fonctionnement des transformateurs sont les rgimes de surcharge,
les rgimes dsquilibrs et les rgimes transitoires. [Gue-1994]
I.8.1 Rgime de surcharge
Les transformateurs dinterconnexion de grandes puissances de rseau lectrique doivent
supporter les variations de la charge du rseau lectrique. Ils sont soumis des surcharges
plus ou moins fortes et sur des dures plus ou moins longues.
Llvation des courants de lignes lors des rgimes de surcharge induisent un accroissement
du flux de fuite, encore amplifi par la saturation du circuit magntique. Les pertes dans les
lments principaux : les bobinages, la cuve et les shunts sont accrues. Ces pertes contribuent
un accroissement de la temprature globale ou locale (points chauds) des lments
constitutifs et de lhuile de refroidissement de lisolation. Ainsi, la raction de dcomposition
des isolants cartonns baigns par lhuile se trouve acclre. ce qui peut diminuer sa dure de
vie et provoquer sa destruction terme.
13
Chapitre I
I.8.2 Rgime dsquilibr

Le rgime dsquilibr correspond un incident sur lune des trois phases du rseau.
La composante homopolaire du courant dans les bobines devient importante, entrainant
un dsquilibre des ampres-tours du circuit magntique et par consquent une augmentation
des pertes supplmentaires dans le fer et le cuivre provoquant des surchauffes locales
considrables.
I.8.3 Rgime transitoire
Les phnomnes transitoires dans le transformateur sont causs dune part par
les enclenchements dautre part par les courts-circuits. Les enclenchements engendrent
des surtensions, qui peuvent dtriorer lisolation. Les courts- circuits crent des surintensits
qui engendrent des chauffements et des efforts lectrodynamiques importants au niveau des
enroulements.
I.9 Protection des transformateurs contre les surtensions
Depuis la construction des premires lignes haute tension le problme de protection
des
transformateurs
contre
les
surtensions
prsentaient
un
intrt
exceptionnel.
Ces protections peuvent tres intrieures ou extrieures. Les mesures de protections internes
sont le renfoncement adquat de lisolation des bobines dentre et dextrmit, o peuvent
apparaitre les plus grands gradients de tension. [Kos-1979]
Les mesures de protections extrieures ont pour but dattnuer londe qui arrive dans le
transformateur en diminuant son amplitude et en la rendant moins raide, au moyen dlments
capables dabsorber une grande partie de lnergie transporte par la surtension en lcoulant
vers la terre. On distingue les deux types rpondus : les parafoudres et les clateurs.
I.10 Diffrentes approches de modlisation des transformateurs
La modlisation des transformateurs consiste tablir une structure mathmatique qui dcrit
lensemble des phnomnes lectromagntiques qui sont rgis par les quations de Maxwell.
Celles-ci constituent un systme dquations aux drives partielles qui lient les phnomnes
magntiques aux phnomnes lectriques, la rsolution de ces quations est intimement lie
la gomtrie et aux caractristiques de la structure tudie.
14
Chapitre I
Selon quon sintresse aux grandeurs locales ou globales diffrentes types de mthodes
peuvent se prsenter :
Mthodes bases sur le calcul numrique
(mthode des Elments Finis) qui fait
intervenir des grandeurs locales lectromagntiques.
Mthodes dites analytiques bases sur un schma quivalent lectrique ou magntique

qui font intervenir des grandeurs globales.
Mthodes dites semi-analytiques bases sur un schma quivalent lectrique fin

qui fait
intervenir des grandeurs locales. La mthode la plus approprie est
la mthode des Circuits Electriques Magntiquement Coupls.

I.11 Schma lectrique quivalent en tant que modle
Un schma quivalent de transformateur s'articule autour d'lments modlisant les
phnomnes magntiques pour le circuit magntique, et les phnomnes lectromagntique
et lectrostatique pour les bobinages, en intgrant le couplage, les inductances, les rsistances
et les capacits. L'laboration d'un schma quivalent suppose de passer d'un problme
gomtrique dcrit dans un espace deux ou trois dimensions (ce que nous appelons un
modle spatial) seulement quelques paramtres caractristiques: les lments localiss d'un
schma lectrique [Rob-1999].
I.11.1 Schma classique
Un schma quivalent complet dcrit bien le fonctionnement du transformateur quon peut
prsenter comme suit :
Fig. I.8 : Schma quivalent dun transformateur rel [Lef-2006]
15
Chapitre I
En dsignant par
primaire de
les rsistances et inductances de fuites des enroulements
spires et du secondaire de
spires ;
rsistance du circuit magntique et
est la ractance de magntisation du circuit magntique correspondant ce flux.

,
sont des capacits parasites qui apparaissent en hautes frquences.
et
Le transformateur est aliment sous une tension variable, dlivrant une tension
et travers par les courants
au primaire et
au secondaire.
I.11.2 Phnomnes lectromagntiques

Il existe un champ magntique dans le bobinage d aux fuites de champ magntique.
Ces lments se traduisent, en plus des rsistances, par deux circuits prenant en compte ces
et
fuites inductances
Dans les enroulements, les phnomnes magntiques peuvent tre rsums par un coefficient
d'inductance liant le flux au courant qui lui donne naissance.
I. 1
I. 2
est le flux commun toutes les spires.

Lorsque les frquences des courants lis la charge augmentent, le trajet des lignes de fuite
est modifi et diminue par la mme occasion les valeurs de inductances
Sachant quune tension variable, par l'intermdiaire du flux, induit une force lectromotrice
non seulement dans le conducteur qui porte ce courant mais galement dans tout conducteur
plac proximit. Ces forces sont donnes comme suit:
I. 3
I. 4
Un transformateur rel est fondamentalement caractris par le systme dquations

ci-dessous:
I. 5
I. 6
I. 7
16
Chapitre I
La rluctance est l'analogue magntique de la notion de rsistance. Elle vaut classiquement

pour un tube de flux de longueur et de section
constante:
I. 8
I.11.2.1 Couplage parfait (transformateur idal)

La fonction premire du transformateur est d'assurer un couplage entre ses deux
enroulements, couplage qui trouve son origine dans l'existence d'un flux commun ceux-ci.
Le couplage est parfait lorsque la totalit du flux est commun aux deux enroulements. Le
transformateur idal nglige les pertes cuivre, aucune hystrsis et les pertes de courant de
Foucault dans les enroulements, et assume le couplage idal [Dix-2004].
Les courants et tensions sont alors lis par le rapport de transformation, not
et dfini par:
I. 9
I.11.2.2 Couplage imparfait

En pratique, il est impossible de coupler parfaitement deux enroulements: du fait que dans les
dispositifs rels il existe toujours un flux de fuite qui rduit le flux commun. L'importance
de cette rduction est caractrise par le coefficient de couplage dont la valeur dans ce cas est
infrieure l'unit.
I. 10
I.11.2.3 Pertes cuivre
Un schma quivalent doit galement modliser les effets dissipatifs dans les enroulements,
savoir les pertes cuivre, qui sont les pertes par effet Joule dans les enroulements,
elles sexpriment par la formule suivante :
,
I. 11
17
Chapitre I
I.11.2.4 Calcul des rsistances (grandeurs lectriques)

La rsistance de chaque enroulement est en fonction de la longueur de fil
, de sa surface
et de sa rsistivit. Sa relation gnrale est donne par :

.
I. 12
, comme on ne connat
Pour obtenir la rsistivit lectrique du matriau
pas la temprature du cuivre ou laluminium a priori, lutilisateur doit lestimer en fonction

de la norme de lenroulement.
1
I. 13
: Conductivit lectrique du milieu
: est la rsistivit lectrique du matriau 0 C
1.72 10 pourlecuivre
.
2.69 10 pourlaluminium
: est le coefficient de temprature du matriau soit
3.93
4.03
10
10
pourlecuivre
pourlaluminium
La dpendance frquentielle des pertes dans le conducteur de lenroulement est lie

aux courants de Foucault. Ces derniers sont dus la variation dans le temps du champ
magntique. En effet les courants de Foucault engendrent une augmentation des pertes
et une rduction du flux magntique total. Ces pertes peuvent se dcomposer suivant deux
origines : les pertes propres dues leffet de peau par le biais du champ magntique cr par
le conducteur sur lui-mme et celles dues aux effets de proximit par le biais dun champ cre
par les autres conducteurs. Leffet de peau et de proximit sont la source de pertes Joule
additionnelles. Les pertes augmentent alors la temprature des enroulements et par consquent
leurs rsistances.
I.11.2.5 grandeurs magntiques
Ceci est schmatis par une ractance
en parallle aux bornes du primaire du
transformateur. La ractance est lie la rluctance quivalente de circuit magntique par la

relation suivante :
.
2
I. 14
I. 15
18
Chapitre I
I.11.2.6 Energie magntique emmagasine

Lnergie stocke dans linductance magntisante
1
2
1
.
2
a pour expression :
. I. 16
Le thorme dampre appliqu au circuit magntique permet dcrire :

I. 17
Ce qui nous conduit lexpression :
I. 18
est la longueur du circuit magntique,
sa section et
la permabilit magntique
absolue du matriau le constituant, est linduction magntique.

Lnergie stocke dans le circuit magntique est inversement proportionnelle la permabilit
magntique du matriau constituant le circuit magntique. Si on veut limiter cette nergie
on a intrt choisir un matriau de permabilit magntique leve.
I.11.2.7 Pertes fer
Les pertes fer sont souvent modlises par une simple rsistance
en parallle sur
l'inductance de magntisation. Elles rsultent de la variation du flux dans le circuit

magntique et se composent des pertes par courant de Foucault (
hystrsis (
et les pertes par
Pour des faibles frquences les pertes fer sont dfinies par la formule semi-empirique
suivante : [Sho-2004]
I. 19
Dans celle-ci, les coefficients
dpendent des particularits de fabrication de la tle
mais galement de la frquence et de la valeur de l'induction. Les coefficients
reprsentent respectivement le volume du noyau et l'paisseur des tles ferromagntiques,

tandis que
peau
est coefficient de Steinmetz (1.7 2). A hautes frquences, quand l'paisseur de
devient trs infrieure l'paisseur des tles, les pertes par courant de Foucault sont
alors proportionnelles la frquence [Lef-2006].

19
Chapitre I
Dans les matriaux ferromagntiques feuillets les pertes fer sont donnes par la formule :
la rsistivit du matriau,
4 . 10
I. 20
permabilit relative et permabilit du vide
).
Pourleshautesfrquences
1.4 ; Pourlesbassesfrquences
1.4 [Fre-
2009].
I.11.3 Phnomnes lectrostatiques
Dans le transformateur il ya diffrentes capacits distribues le long des enroulements ; entre
conducteurs ; entre enroulement haute et basse tension ; entre bobinages et cran ou masse
mcanique. Ces lments sont la consquence du champ lectrostatique dont leffet est
beaucoup accentu en haute frquence.
I.11.3.1 Pertes dilectriques
Les pertes dilectriques dues des courants de dplacement dans les isolants sont pour
la
plupart
du
de dissipation
temps
ngligeables
par
rapport
aux
autres
pertes.
Le
facteur
estdelordre10 10 . Ces pertes sont importantes lorsque
les lments actifs interviennent, car leur isolation contient une teneur dhumidit leve
[Oli-2003], et peuvent tres modlises par un rseau de rsistances et de condensateurs en
cascade et en parallle [Ahm-92], leur expression est de la forme suivante :
. .
. I. 21
: est la tension en volt.
C : est la capacit (F).
: est la pulsation (rad/s).
I.11.3.2 Calcul des capacits (grandeurs lectrostatiques)

Le calcul des ces capacits dpend de la gomtrie et des caractristiques lectriques des
spires des enroulements (Fig. I.12). [Ven-2009].
Certaines de ces capacits sont calcules avec des formules gomtriques simplifies et
dautres sont calcules avec des formulations semi-empiriques.
20
Chapitre I
Noyau
Enroulement BT
Enroulement HT
Fig. I 8 : Reprsentation de leffet lectrostatique par diffrentes capacits [Ven-2009]
I.11.3.3 Capacit entre enroulements
Les enroulements n'tant rien d'autre que des conducteurs proches ports des potentiels
diffrents, ils introduisent des effets capacitifs dont les principales manifestations sont
des rsonances avec les inductances du transformateur.
En considrant une configuration cylindrique la capacit entre une spire de lenroulement
haute tension et lenroulement basse tension est calcule au moyen : [Ven-2009]
2
I. 22
: Capacit entre enroulement de BT et HT

: Rayon de lenroulement HT
: Permittivit du vide (
10 36
: Rayon de lenroulement BT
: Permittivit relative du milieu considr
I.11.3.4 Capacit entre conducteur

La capacit entre conducteurs dune mme spire ou bien entre conducteurs de spires
diffrentes est approche par lquation suivante :[Ven-2009]
I. 23
21
Chapitre I
I.11.3.5 Capacit entre phases

La formule suivante est applique pour capacits entre de jambes adjacentes.
I. 24
1
Avec
b : la distance entre deux enroulements adjacent (centre-centre),
est le diamtre
de lenroulement.
Fig. I 10 : Dimension entre deux enroulements adjacents et entre lenroulement

et la cuve [Bje-2005]
I.11.3.6 Capacit entre enroulement et cuve
La capacit entre lenroulement et la cuve peut tre calcul partir de lquation I. 24 ,
avec un facteur rduit de 0.25 pour une configuration cylindrique [Bje-2005] :
0.25.
I. 25
1
D est la largeur interne de la cuve.
La capacit entre lextrmit de la phase et la cuve peut tre approche par [Bje-2005] :
Cela pour un facteur rduit
de 0.75 (si la cuve nest pas compltement circulaire).
0.75.
I. 26
22
Chapitre I
I.12 Modlisation par la mthode semi-analytique

La mthode semi-analytique est base sur des techniques de discrtisation numriques en
petits lments, sur lesquelles on applique des formules analytiques. Cette mthode nous
permet davoir un modle de schma lectrique fin et dtaill et le calcul des ces paramtres
est plus prcis. Parmi ces mthodes on peut citer la mthode des circuits lectriques
magntiquement coupls. Elle est dun emploi restreint aux milieux linaires, et dispositif
avec gomtrie simple et homogne [Bje-2005].
I.13 Modle dtaill (fin) bas sur le calcul des inductances et les rsistances
La mthode est base sur un rseau dlments compos des selfs et des inductances
mutuelles. Quel que soit le type de reprsentation, les lments sont calculs en fonction des
caractristiques gomtriques du bobinage ainsi que des caractristiques gomtriques du
noyau magntique. Lobjectif de cette mthode cest dtablir un modle dun circuit
lectrique quivalent dtaill et fin des enroulements illustr sur (Fig. I.11). [Ahm-1992]
[Bje-2005] [Ras-2008]
Fig. I.11 : Modle dun circuit lectrique quivalent

de lenroulement primaire [Bje-2005]
23
Chapitre I
I.13.1 Calcul des selfs inductances

Pour le calcul de la self-inductance d'une spire de section rectangulaire on peut utiliser
une formule approximative base sur la gomtrie des enroulements donne comme suit:
8
ln
2 I. 27
: rayon moyen de la spire

: rayon gomtrique moyen de la spire donn par Grover : [Rah-2003] [Bje-2005] :
2
3
ln
2
3
12
12
25
I. 28
12
Fig. I.12 : Dimension dune spire de section rectangulaire
Les paramtres
et
sont respectivement
le rayon moyen, la hauteur et la largeur
de la spire.
I.13.2 Calcul des inductances mutuelles
Une formule exacte pour la mutuelle inductance entre llment et dsign par
est
donne par Maxwell en utilisant les intgrales elliptiques: [Shi-2004] [Bje-2005]
.
4
I. 29
I. 30
et rayons de llment et
;
integrales lyptiques de 1er et 2eme ordre.
Fig. I.13 : Mutuelle inductance entre deux lments filamentaires

circulaires [Shi-2004]
24
Chapitre I
I.13.3 Calcul des rsistances

La rsistance dune spire lmentaire est calcule au moyen :
I. 31
:rsistivit de llment i ;
et
sont respectivement la longueur et la section de llment.
I.14 Modlisation par la mthode des lments finis (MEF)

La mthode des lments finis est une mthode incontournable lorsque l'on traite
la modlisation des milieux continus. Son domaine d'application trs tendu fait probablement
d'elle le premier outil de conception. Ainsi, applique au domaine de l'lectromagntisme, elle
permet d'tudier des problmes aux formes complexes, avec des matriaux linaires ou non,
des couplages circuits, du mouvementetc. Du fait de sa prise en compte des phnomnes au
niveau local, sans injection de connaissance a priori, elle est la mthode utiliser pour
comprendre finement les phnomnes physiques en jeu dans un dispositif. Lobstacle souvent
rencontr par cette mthode dans beaucoup dapplications est la prise en compte des
conditions aux limites existant linfini : lutilisation dun maillage pour approcher le
problme rel demande un nombre de nuds important et prjudiciable au temps de calcul.
De plus cette mthode trouve ses limites lorsque le domaine dtude prsente de singularits.
I.15 Conclusion
Ce premier chapitre nous a permis de prsenter les principaux types de transformateurs,
leurs rgimes de fonctionnement. Du point de vue modlisation, diffrentes mthodes on ts
dveloppes entre autres le schma lectrique quivalent. Ce modle permet de reproduire
certains comportements lectromagntiques du transformateur. Mais masque en ralit des
phnomnes plus complexes et comme tout modle souffre d'hypothses simplificatrices. Ces
hypothses conduisent cependant des paramtres de reprsentation trop idaliss pour les
applications industrielles actuelle. Les difficults de modlisation ncessitent le recours aux
mthodes numriques ou semi-analytique de discrtisation en petits lments de circuit telles
que la mthode des lments finis, et la mthode des circuits coupls. Ces modles sont les
plus rpandus et utiliss, l'heure actuelle, pour modliser un transformateur. Quel que soit le
type de reprsentation, les lments sont calculs en fonction des caractristiques
gomtriques et dilectriques du bobinage ainsi que des caractristiques gomtriques du
noyau magntique.
Le prochain chapitre sera consacr la modlisation lectromagntique des transformateurs.
25
Chapitre II
Modlisation lectromagntique des
transformateurs
Chapitre II
Modlisation lectromagntique des transformateurs
II.1 Introduction
Nous prsentons dans ce deuxime chapitre, les bases mathmatiques des quations
lectromagntiques pour la modlisation dun problme lectrotechnique, travers le rappel
des quations de Maxwell, ainsi que les diffrentes formulations magntostatiques
et magntodynamiques qui en rsultent. Ces formulations conduisent des quations aux
drivs partielles auxquelles il faut associer des conditions aux limites et des conditions
dinterfaces.
La rsolution analytique de telles quations ncessite gnralement des hypothses
simplificatrices fortes, et ne rend donc pas compte de la ralit des phnomnes physiques.
Dans le souci de prendre en considration la complexit des structures des systmes
lectromagntiques (non linarit, phnomnes coupls, gomtrie 3D), on a recours soit aux
mthodes couples analytico-numrique soit purement numriques.
II.2 Equations de Maxwell
La rpartition spatiale et temporelle des champs magntique et lectrique est donne
par les quations de Maxwell: [Dha-2005] [Ped-2003] [Hen-2004]
.

,-&- .
0
0 )
.
)

#$%

123$
Avec :
: Champ magntique ;A/2>

4
' : Champ lectrique ;V/2>
'
$&
56 7
4
$&

.
.

.

()

)

.
* ' .
(&
&

.
.

(
(
85 7

*4
9

(&
(&
6

II. 1"
II. 2"
II. 3"
II. 4"
: Induction magntique ;T>

)
: Induction lectrique ;A/2B >

56 : Densit de courant de conduction ;1/2B >
: Densit de charge volumique ;A/2C >

: charge contenue dans le volume V ;A>
: Elment de volume ;2C >
& : Le temps ;>

: Elment de surface ;2B >
: Elment de longueur ;2>
Ces quations sont associes aux relations constitutives des matriaux qui sont en gnral,
anisotropes ou non linaires
26
Chapitre II
II.3 Lois de comportement des matriaux
Une
relation
constitutive
dcrit
localement
le
comportement
des
grandeurs
lectromagntiques dans un matriau donn. Elles sont donnes dans le cas le plus gnral

2 -$ &$3
G H'
E
KH , %, M"K '
Dans un milieu dilectrique
2 -$ -&$3
F
E
HO' P'
2 -- -$ &$3
D

-$ &$3
G ) S4 7 )T
E
7 )
T
KS , %, M"K 4
Dans un milieu magntique )
-$ -&$3
F
E
SO4
P4
7 )
T
-- -$ &$3
D)
H : Permitivit lectrique absolue ;F/2>
II. 5"
II. 6"R
II. 7"
II. 8"
II. 9" R
II. 10"
S : Permabilit magntique absolue ;4/2>

T : induction rmanente ;X>
)
II.4 La loi dOhm
Cest la loi exprimant, dans un milieu conducteur la densit de courant 5 en fonction
du champ lectrique ' , la loi sapplique tout milieu prsentant une conductivit lectrique.
5, Z'
Y
5, Z' 7 5
2 ,-,&$
2 ,-,&$ , $,
5 : Densit de courant source ;1/2B >
27
II. 11"
II. 12"
Chapitre II
II.5 Relations de passage
par des relations de continuit. Soient deux milieux [\ et [B , et soit -\B le vecteur normal la
Les champs lectriques et magntiques des diffrents milieux sont lis linterface
surface dirig du milieu [\ vers le milieu [B . (Figure II.1)
La frontire entre [\ et [B
[\
S\
[B
-\B
H\
Z\
SB
HB ZB
Fig. II-1 : Interface entre deux milieux [Doi-2007]
due aux
La discontinuit du milieu de la composante tangentielle du champ magntique 4
.
densits de courants surfaciques ]
\ 4
B P ^ - ]

O4
II. 13"
O'\ 'B P ^ - 0
II. 14"
\
B P. -
O
II. 15"
\ )
B P. - 0
O)
II. 16"
La conservation de la composante tangentielle du champ lectrique '
due aux charges surfaciques

La discontinuit de la composante du dplacement lectrique
.

La conservation de la composante normale de linduction magntique )
28
Chapitre II
II.6 Equations de Maxwell simplifies

Les hypothses simplificatrices gnralement poses dans le traitement des quations de
modlisations des problmes Electrotechniques sont :
Dans le domaine de llectrotechnique, les frquences de travail sont infrieures aux

radios frquences qui sont de lordre 1012 Hz, ce qui nous permet de ngliger les courants
(& devant les courants de conductions 5 dans lquation II. 4",

de dplacement 5_ (
,
le systme est quasi-stationnaire.
Les charges despace dans les conducteurs sont ngliges dans lquation II. 1", ainsi
pour les conducteurs le rapport (H/Z) est trs faible, de lordre de 10a\b a\ . ce qui est
ngligeable devant celui quon considre pour les frquences infrieures 10\B Hz.
En considrant les hypothses simplificatrices prcdentes, les quations de Maxwell

0

II. 17"
Z' 7 5
4
$&
II. 20"
scrivent sous la forme suivante :

(&
' ()
$&
II. 18"
0
)
II. 19"
II.7 Modles lectromagntiques

Sur la base des quations de Maxwell simplifies II. 17 20" relatives au domaine de
llectrotechnique nous dfinirons les modles magntostatique et magntodynamique.
II.7.1 Modle magntostatique
/(&) est nul et les champs
par des sources de courant indpendantes du temps. Le terme (()
Dans le cas de la magntostatique, on supose que le champ magntique est produit
sont dcoupls.
lectrique ' et magntique )
Le systme d'quations rgissant les phnomnnes magntostatiques s'crit :

0
)
5

$&4
II. 21"
II. 22"
A ce systme, on peut ajouter la loi de comportement des matriaux magntiques II. 8", ainsi
que les conditions aux limites.
29
Chapitre II
II.7.2 Modle magntodynamique
Dans le cas magntodynamique, les sources de tension ou de courant varient en fonction

du temps, cela veut dire que les courants induits rsultants feront le couplage du champ
lectrique
et
magntique.
Le
systme
dquations
rgissant
les
phnomnes
magntodynamiques scrit :
'
$&

()
(&
II. 23"
Z' 7 5
4
$&
II. 24"
Il faut ajouter ce systme les lois de comportement des matriaux lectriquesII. 5)

et magntiquesII. 8" et les conditions aux limites.
II.8 Choix de la variable dtat
Pour traiter les quations magntostatique et magntodynamique diffrentes formulations

mathmatiques du problme peuvent tres obtenues, soit en utilisant les champs magntiques
ou lectriques comme variables principales, ou en introduisant de nouvelles variables que sont
les potentiels lectriques ou magntiques. A chacune de ces grandeurs magntiques et
lectriques correspond une quation aux drives partielles.
Les formulations en champs prsentent un inconvnient majeur, qui est la discontinuit aux
interfaces et lutilisation dlments darrtes dans la discrtisation. De plus les vecteurs ont
gnralement trois composantes, augmentant par-l le nombre dinconnues, ce qui entraine un
temps de rsolution plus important. Pour remdier ces inconvnients on fait appel aux
potentiels vecteur/scalaire lectriques ou magntiques.
Il existe plusieurs types de potentiels, parmi les quels, on peut citer le potentiel vecteur
magntique et le potentiel scalaire lectrique. La rsolution des quations en potentiel vecteur
magntique,
savre
tre
particulirement
intressante
dans
des
configurations
axisymtriques, car seule la composante azimutale est non nulle. La condition de Jauge de
continuit qui assure lunicit du potentiel vecteur magntique A est naturellement vrifie
[Lef-2005] [Doi-2007].
30
Chapitre II
II.9 Formulation en potentiel vecteur magntique
issu de
Cette formulation est base sur lexistence dun potentiel vecteur magntique A
lquationII. 3", tel que :

$&

)
1
La combinaison des quations II. 2" et II. 25" conduit :

^ 8' 7
c
(1
90
(&
II. 25"
II. 26"
La relation II. 26" nous permet de dduire quil existe un potentiel lectrique scalaire V tel
que :
'
(1
"
d$
(&
partir des quations II. 4", II. 25" et II. 27" nous pouvons crire :

GZ 8(19 Z.
d$ i&- 2d-&%-2i
(&
E
E
(1
1

9
i&- 2d-&%-2i - $,
Z
8

$& g $&1h
(&
S
F

i&- 2d-&&&i , $,
E Z. d$
E
0
i&- 2d-&&&i - $,
D
II. 27"
II. 28"
II. 29"R
II. 30"
II. 31"
Lquation II. 28" reprsente la formulation magntodynamique en terme de 1 tenant compte
du terme source.
II.10 Conditions de Jauge

Les champs lectromagntiques dfinis par lensemble des quations de Maxwell
ne sont pas uniques. En effet les champs divergence sont dfinis un rotationnel prs et les
champs rotationnel sont dfinis un gradient prs. Il convient donc dimposer une condition
dquations II. 28"-II. 31" . Cette condition, appele condition de Jauge est gnralement
supplmentaire afin dassurer lunicit de la solution issue de la rsolution du systme
exprime sous deux formes [Doi-2007].
31
Chapitre II
1" 0
1" 7 l
jd A2k
(
0
(&
jd $-&M
II. 32"
II. 33"
II.11 Formulation bidimensionnelle en potentiel vecteur magntique

Le terme gauche de lquation II. 28" scrit :
Om d$
Om $&
1P n
$&
1P 7 d$
Om n
1P
Y
1
m
S
II. 34"
En adoptant la condition de Jauge de Coulomb donne parII. 32", lquation II. 34" devient :
Om d$
Om $&
1P n
$&
1P Z 8
(1
9 7 5
(&
Si on considre que la source est harmonique sinusodale lquation II. 35" scrit :
Om d$
n
1P 7 Zjo1 5
II. 35"
II. 36"
II.11.1 Formulation 2D cartsien

Dans le cas dun problme 2D cartsien , %" , le potentiel vecteur magntique 1 et la densit
de courant 5 nayant quune seule composante suivant oz. Ces derniers sont de la forme :
1 0,0, 1p "
et 5 0,0, 5p ".
En coordonnes cartsiens nous avons :

(
1 qq
$&
(
0
r
(
(%
0
s
(q
q
(M
1p
II. 37"
, r , s : Vecteurs unitaires.

$&1 t
(1p "
(1p "
u 7 t
u r
(%
(
Om $&
1P t
$&
(
(1p
(
(1p "
gm.
h gm.
hu p
(
(
(M
(%
32
II. 38"
II. 39"
Chapitre II
Lquation II. 36" devient :

(
(1p
(
(1p
gm.
h
gm.
h Z. j. o. 1p 7 5p
(
(%
(
(%
II. 40"
II.11.2 Formulation 2D axisymtrique
Dans le cas dun problme 2D axisymtrique$, M", le potentiel vecteur magntique 1 et la

densit de courant 5 nayant quune seule composante azimutale.
1 0, 1v , 0" et 5 0, 5v , 0".
En coordonnes cylindriques axisymtriques nous avons :

T
1 (
1 . q
$&
$ ($
0
$v
(
(w
$1v
p
(
q
(M
0
II. 41"
( m ($1v "
(
(1v "
g .
h gm.
hu v
($ $
(M
($
(M
II. 42"
(
(1v "
( m ($1v "
gm.
h g .
h Z. j. o. 1v 7 5v
(M
(M
($ $
($
II. 43"
Om $&
1P t
$&
Lquation II. 36" devient :

II.12 Les conditions aux limites

Les valeurs dun champ 1 la frontire du domaine satisfait principalement deux relations
simples. Le champ peut tre soit de valeur connue, soit normal la surface, ce qui de traduit
par :
R1|y 1z
(1
R | 1z
(- y
II. 44"
,-&- $,{&
II. 45"
,-&- }2--
Les conditions de Dirichlet et de Neumann sont dites homognes si les valeurs imposes 1z
sont nulles. Elles sont dites non homognes dans le cas contraire.
33
Chapitre II
II.13 Mthodes de rsolution

II.13.1 Mthodes numriques
Lutilisation des mthodes numriques de discrtisation pour la rsolution des quations
mathmatiques tablies, consiste ramener la rsolution des quations aux drives partielles
dans le domaine dtude, compte tenu des conditions aux limites, celle dun systme
dquations algbriques dont la solution donne les valeurs et la distribution des grandeurs
recherches. Parmi ces mthodes on trouve la mthode aux diffrences finies, la mthode des
volumes finis et la mthode des lments finis. La mthode des lments finis est sans doute
la plus utilise car elle est mieux adapte pour traiter les gomtries complexes et les milieux
non linaires.
II.13.1.1 Principe de la mthode des lments finis
Dans les formulations prsentes, il sagit de rsoudre des quations aux drives partielles
sur le domaine , auxquelles sont associes des conditions aux limites sur la frontire . La
mthode des lments finis ne sapplique pas directement pour la rsolution des quations aux
drives partielles, mais une formulation intgrale du problme, en utilisant lune des deux
approches suivantes :
La mthode variationnelle exige la connaissance au pralable de la fonctionnelle dnergie

du systme tudi, elle sexprime par :
.
II. 46"
#1" .

6
: Lagrangien construit partir de la diffrence entre lnergie cintique et
lnergie potentielle du systme tudier.
La mthode des rsidus pondrs ou mthode projective qui consiste minimiser le rsidu
induit par lapproximation de la fonction inconnue. Cest une mthode plus gnrale, son
application ne ncessite pas la connaissance de la fonctionnelle dnergie du systme,
elle traite directement lquation aux drives partielles quon veut rsoudre.
Lquation rsoudre tant :
(#1"
(.1" (.1"
(.1"
0

0
(1\
(1B
(1
(1
-: Nombre de nuds du domaine dtude

1 : Inconnu au nud du domaine.
34
II. 47"
Chapitre II
Linconnue A du problme est approxime sur les lments de maillage par la fonction
dinterpolation r telle que :
1r
r . 1r
r
r : dpend des coordonnes des nuds correspondant chaque lment de maillage.
II. 48"
II.13.2 Mthodes semi-analytiques

Parmi les mthodes intgrales dites semi-analytiques, on distingue la mthode des circuits
lectriques magntiquement coupls (CEMC), et la mthode des intgrales de frontires.
Contrairement aux modles bass sur les mthodes purement numriques qui incluent lair
environnant dans le domaine de rsolution, les mthodes semi-analytiques saffranchissent de
la discrtisation de lair. La mthode des circuits coupls est une mthode de rsolution semi
analytique qui sapplique quaux milieux conducteurs. Applique notamment pour ltude des
systmes induction structure asymtrique, dont le fonctionnement est bas sur les courants
de Foucault [Meu-2002].
II.13.2.1 Principe de la mthode des circuits coupls
La mthode des circuits lectriques magntiquement coupls consiste subdiviser
le conducteur en lments de circuits circulaires et coaxiaux fictifs, dont la forme
et les dimensions sont telles que la densit de courant peut tre suppose constante dans
chacun de ces lments. La subdivision peut se faire en mailles non rgulires selon la pice
modlise, zone de forte variation des grandeurs physiques, effet de peau important.
Axe de symtrie
\
B
Elments de
discrtisation - de la
spire }
Fig. II.2: Dcoupage des conducteurs en lments [Lef-2006]

35
Spire lmentaire
Chapitre II
III.13.2.2 Interaction entre les lments de circuits lectriques

Chaque lment du circuit lectrique est reprsent par sa rsistance lectrique et son
inductance, ces lments forment entre eux un systme de circuits mutuellement coupls
[Mao-1996] [Che-1994] [Doi-2007]. La discrtisation par la mthode des circuits coupls est
reprsente dans la (Fig.II.3), o tout conducteur d'une bobine est forme de N spires
principales en srie et chacune de ces spires est constitues de n spires lmentaires en
parallles.
\\
\\
\
B\ B\
\\
B\
B\,\
3$"1"
\
\B
\B
B
BB
BB BB
\\
B
3$"2"
\ \
\
B

3$"}"
Fig. II.3: Schma lectrique quivalent pour le primaire [Lef-2006]
Le systme rsultant est alors rgi par - } quations, - } inconnues qui sont
- courants et les } tensions . Chaque spire lmentaire est soumise une tension
identique la tension la quelle est soumise la spire principale.
Lapplication de la loi de Kirchhoff au schma quivalent de la (Fig. II.3) nous permet

dcrire :
Dans la spire 1
\\
s
G \ \\ \\ 7 \\
7
\\,s
&
&
s\
\
E
.
E
.
F
.
E

s
\
E 7
7
\,s
\
\ \
\
D
&
&
s\
\
36
II. 49" R
Chapitre II
Dans la spire 2
\B
s
G B \B \B 7 \B
7
\B,s
&
&
s\
\
E
.
E
.
F
.
E

s
B
E 7
7
B,s
B
B B
B
D
&
&
s\
\
II. 50" R
Dans la spire N
\
s
G \ \ 7 \
7
\,s
&
&
s\
\
E
.
E
.
F
.
E
s
E 7 7
,s
D
&
&
s\
\
II. 51" R
l dsigne la spire lmentaire et est la spire principale et lensemble des quations de circuit
II. 49" II. 51" conduit au systme dquations diffrentielles linaire suivant :
;>;> 7 ;>
(
;> ;>
(&
;> : Vecteur des courants dans les spires filiformes
II. 52"
;> : Vecteur de tension dans les spires filiformes
;> : Matrice diagonale, ses lments reprsentent les rsistances des spires lmentaires
;> : Matrice inductance pleine, llment de la diagonale reprsente linductance propre
de la spire lmentaire et llment hors diagonale r reprsente linductance mutuelle
entre la spire et la spire j.
37
Chapitre II
III.13.2.3 Alimentation sinusodale

Dans le cas des phnomnes lectromagntiques sinusodaux, le systme dquations obtenu
est linaire et coefficients complexes, il est de forme :
;>;> ;>
II. 53"
Les diffrents oprateurs de diffrentiation temporelle sont remplacs par le coefficient

jo .
;>

7 jo
s

II. 54"
;> : Matrice impdance de dimension -} -}"

: Rsistance de la spire ;
: Inductance propre de la spire
s : Mutuelle inductance entre la spire et la spire l

II.14 Conclusion
Les quations gnrales rgissant les phnomnes lectromagntiques sont dcrites par des
quations aux drives partielles. La rsolution de ces quations permet de connaitre
lvolution du champ en tout point de lespace. Nous nous somme intresss dune faon
particulire aux modles magntostatique et magntodynamique qui sont obtenus laide du
potentiel vecteur magntique. Ces quations peuvent tes rsolues par des mthodes
analytiques ou mthodes numriques ou bien semi-analytiques. Parmi ces mthodes on
retiendra la mthode des circuits lectriques magntiquement coupls, et un modle
numrique bas sur une approche lments finis qui prend en compte la non linarit et
lanisotropie des matriaux, les grandes lignes de ces deux mthodes ont ts prsentes.
Dans le chapitre qui suit nous dtaillerons la mthode des circuits lectriques magntiquement
coupls et lappliquerons notre dispositif.
38
Chapitre III
Modlisation des transformateurs
par une mthodologie mixte
Circuits Electriques Magntiquement
Coupls (CEMC) et Elments
Finis(EF)
Chapitre III
Modlisation des transformateurs par une mthodologie mixte

Circuits Electriques Magntiquement Coupls (CEMC)
et lments finis (EF)
III.1 Introduction
La mthode des circuits lectriques magntiquement coupls est base sur une discrtisation
numrique du domaine conducteur en lments de circuit. Chaque lment ou spire
lmentaire est reprsent par sa rsistance et son inductance. Pour cela les bobinages HT et
BT sont discrtiss en lments lectriques correspondant une portion (spires lmentaires)
auxquels sapplique des expressions analytiques des grandeurs lectromagntiques, dduites
de la loi de Biot et Savart. Le calcul de diffrents effets rsistifs et inductifs, propres et
mutuels entre ces multiples spires et ventuellement le circuit magntique est exprim partir
de lexpression analytique du potentiel vecteur magntiqueA.
III.2 Modle gomtrique
Notre dispositif est un transformateur triphas, colonnes, pour lequel on modlise
le bobinage dune colonne latrale. Le dispositif est ainsi approch par un systme
axisymtrique est ayant une symtrie de rvolution, permettant de prendre en considration
une partie du domaine dtude tel que reprsent par la Fig. III.1.
Enroulement HT
Enroulement BT
Noyau magntique
Circuit magntique
Phase A
Phase B
Secondaire
Phase C
Zoom
Air
Fig. III.1.a : Architecture du transformateur

[Hoc-2009]
39
Primaire
Fig. III.1.b : Domaine dtude
Chapitre III

III.3 Discrtisation des enroulements

Classiquement il existe deux modles de conducteurs : les conducteurs filamentaires et les
conducteurs massifs. Cependant un nouveau type de conducteur trs efficace est appliqu
aux transformateurs, il sagit de feuilles denroulements, rsistant aux efforts axiaux et aux
tempratures leves dues aux courts circuits [Ger-2001].
On peut distinguer deux sortes de discrtisations : dans le cas ou lpaisseur de peau est trs
faible (hautes frquences), on ne discrtise que le pourtour de la section du conducteur
(Fig.III.2), et dans le cas contraire en (basses frquences), on discrtise toute la section
(Fig.III.3).
Section circulaire
Section rectangulaire
Fig.III.2 : Discrtisation linique de la frontire du conducteur
Section rectangulaire
Section circulaire
Fig.III.3 : Discrtisation totale de la section du conducteur [Mao-1996]

Si lpaisseur de peau est grande par rapport aux dimensions des conducteurs et pour des
frquences gales 50 Hertz, les densits de courant sont pratiquement uniformes comme
cest le cas pour les rgions filaires pour lesquels toute la section du conducteur est
discrtise. Par contre lorsque lpaisseur de peau est petite par rapport aux dimensions des
conducteurs, dans ce cas les pertes produites par les flux de fuite dominent, on ne discrtise
pratiquement que le pourtour de ces derniers.
40
Chapitre III

III.4 Formulation en potentiel vecteur magntique

Cette mthode consiste formuler la loi dOhm pour chaque lment, de faon aboutir
un systme dquations diffrentielles, faisant intervenir les chutes de tension rsistives et
inductives correspondant dune part linductance propre et dautre part aux mutuelles
inductances avec tous les autres lments.
Pour la mise en uvre de cette mthode, des hypothses ont t faites :
La gomtrie tudie est symtrie de rvolution.
Lhlicit des enroulements est nglige.
Les matriaux utiliss sont proprits physiques isotropes et linaires.
Les modles diffrent selon la manire dintroduire la source, le premier suppose

une source de courant et le second prend comme source une tension. Chacun de ces modles
est formul suivant la frquence dexcitation, des caractristiques lectriques, magntiques
et gomtriques du problme. Le modle dexcitation en tension pour objectif
la dtermination de la rpartition de la densit de courant circulant dans chaque spire.
III.4.1 quations caractristiques en axisymtrique
Pour dcrire le phnomne lectromagntique, nous choisissons de ramener les quations

de Maxwell un systme de deux quations exprimant le couplage entre la densit du courant
et le potentiel vecteur magntique, les quations de Maxwell mnent au systme suivant :
quations II. 36 et III. 1 [Mao-2006] [Mao-2009] [Esl-2010]
III. 1
densitdecourantinduit
41
Chapitre III

Spire
Fig.III. 4: Coupe transversale et radiale dune spire [Mao-2006]

Le fait que le courant circule le long du primtre moyen de la spire on ramne le potentiel
scalaire lectrique la variation angulaire tel que :
1
III. 2
La variation de la tension suivant le rayon est nulle lquation III. 2 devient :

1
III. 3
La tensionU au borne de la spire est donne par :

2 .
III. 4
En utilisant le thorme de stockes (quation III. 1 ) on peut crire :

1
2 .
2 . III. 5
A partir des quations III. 2 , III. 3 , III. 4 , III. 5
on obtient lquation gnrale
pour un lment de circuit.
2 .
2 .
III. 6
En rgime harmonique III. 6 scrit :

2 .
2 .
III. 7
42
Chapitre III

On reconnat lexpression de la tension aux bornes dune spire lmentaire utilise

dans la formulation de la mthode des circuits coupls, exprime en termes de densit
de courant et du potentiel vecteur . Le terme 2 .
rsistive, alors que le terme
correspond la chute de tension
reprsente la chute de tension inductive.
2 .
Le potentiel vecteur magntique comprendrait la contribution du courant de la spire

elle-mme (inductance propre), ainsi que la contribution des courants circulant dans les autres
spires (inductances mutuelles), cela est reprsent sur la (Fig.III.5) :
Fig.III.5: Schma lectrique quivalent des deux enroulements HT et BT

Lexpression III. 8 , peut tre crite pour toute spire lmentaire
de la spire
avec la contribution des courants circulant dans les spires lmentaires de la spire comme
suit : [Lef-2005]
2
Nombre de spire des enroulements primaire /secondaire
Nombre de discrtisation par spire primaire / secondaire*

43
III. 8
Chapitre III

Spire
Spire
Fig.III.6 : Reprsentation des spires lmentaires
et
[Leo-1992]
III.5 Expression du potentiel vecteur magntique

Nous considrons le cas de la (fig. III.7), dans lair, qui exprime le potentiel cre par la spire
de contour
contour
parcourue par un courant
au point .
est la distance entre lments du
Fig. III.7: Potentiel vecteur magntique cre par une spire filiforme
[Mao-1996] [Cha-2003]
44
Chapitre III

Le potentiel vecteur
du champ lectromagntique entre deux nuds metteur et rcepteur
est dfini par l'expression loi de Biot et Savart donne par lquation: [Mao-2006]
.
III. 9
Lapproximation du courant est constant dans la section de la spire :

,
Dans
. .
III. 10
une configuration axisymtrique, le potentiel vecteur magntique et la densit
de courant lectrique sont rduits leurs composants angulaires. [Leo-1992] [Mao-1996]

,
:
III. 11
au moyen de :
2
|
On dfini
III. 12
2
. cos III. 13
Le rarrangement de lquation III. 11 nous permet dcrire
III. 14
Avec :
4
III. 15
2 ,
2
45
Chapitre III

Lintgrale III. 14 peut scrire :
2. 2
III. 16
En considrant lidentit [Mao-1996] :

1
Le potentiel vecteur magntique sexprime par lintermdiaire des intgrales elliptiques de

premier et de deuxime espce
2.
1
2
1
8
2 .
III. 18
comme suit [Rez-1992] :
1
2
2.
III. 19
Correspond alors au potentiel vecteur magntique au sein de la boucle lmentaire
rceptrice
de la spire , cr par la densit de courant
mettrice de la spire (Fig.III.6).

46
1;
peuvent tres calculs par une mthode des gomtries arithmtiques, qui consiste
calculer une sries de deux coefficients
III. 17
1
et
et
dans la boucle lmentaire
Chapitre III

Dans le cas ou lpaisseur de peau
est trs petite devant la largeur des spires lmentaires
on intgre par rapport llment de tangente , en tenant compte de la dcroissance

exponentielle du courant la distance
de la normale
do [Che-1995] [Che-1996] :
III. 20
Lquation III. 19 peut se ramener :
2.
III. 21
. .
. .
III. 22
On pose
.
. .
4
. .
4
III. 23
Lquation III. 22 devient : [Mao-1996]
III. 24
La tension totale aux bornes de la bobine est la somme des tensions
aux bornes de chaque
spire. Le nombre de tensions inconnues est la somme des tensions de deux bobines (
[Lef-2005] [Lef-2006].
47
Chapitre III

III. 25
La conservation de courants dune spire lautre et exprime par lquation :

.
III. 26
Dans le but de diminuer la taille du systme dquations rsoudre et comme on ne

sintresse quaux densits de courant, on limine les tensions lmentaires, par consquent
on rduit le systme matriciel.
Daprs
III. 24 , III. 25 , III. 26
on obtient un systme de
,quations
inconnues de la forme :
.
.
.
.
.
.
..
.
.
.
.
.
.
0
0
.
.
.
.
.
..
III. 27
0
0
.
.
.
.
.
: Matrice reprsentant linfluence entre spires du primaire.

: Matrice reprsentant linfluence entre spires du secondaire.
48
Chapitre III

A partir de la loi de conservation du courant III. 26 nous construisons les matrices

de dimension
qui reprsentent respectivement la matrice sections spires
primaire et secondaire tel que : [Doi-2007]
0 0

III. 28
0 . . .
0 0

III. 29
0 . . . 0
A partir de lquation III. 24 , et les quations de conservation de tension (III.25)

et du courant (III.26)
nous obtenons les matrices
et
caractrisant
leffet des conducteurs sur eux-mmes qui dfinissent le couplage entre les spires
lmentaires
des
conducteurs,
. La matrice
1 ,
respectivement
de
dimension
1 ,
est forme en assemblant les deux matrices
et
tel que : [Doi-2007]
La matrice
III. 30
est forme en assemblant les deux matrices
et
III. 31
: Matrice reprsente linteraction entre spires du primaire et celles du secondaire de

dimension
1 ,
: Matrice transpose de mme dimension que
III.6 Prise en compte du circuit magntique

Lutilisation de circuits magntiques (tle ferromagntique, ferrites pour des frquences
plus leves) peut prsenter un certain nombre davantages, notamment lorsquil sagit
de concentrer localement le champ magntique ou de dvier celui-ci vers une rgion
bien prcise. La prsence du circuit magntique modifie le potentiel vecteur magntique.
49
Chapitre III

Certain matriaux tels que lacier, le Nickel ou le Cobalt se caractrisent par une permabilit
magntique dpendante du champ magntique. Ces matriaux sont fortement utiliss
dans lindustrie, il est donc important de prendre en compte leur spcificit.
Nous considrons que la saturation du circuit magntique nest pas atteinte ; ainsi la
permabilit magntique est suppose constante. Les courants de Foucault sont donc ngligs
dans le circuit. Ce dernier se comporte donc comme un matriau homogne et isotrope.
Lensemble de ces considrations nous permet de modliser toute partie de la culasse
par des courants fictifs surfaciques. [Dur-1968].
III.6.1 Modlisation du noyau par des densits superficielles de courant
La mthode des circuits lectriques magntiquement coupls base sur une discrtisation
de linducteur en spires lmentaires et une discrtisation en segments du circuit magntique,
donnant ainsi un maillage des diffrents domaines conducteurs. A chaque spire est associe
dune part lquation magntodynamique (II.28) et les lois de Kirchhoff pour les
enroulements, et dautre part les quations de Fredholm de 2me espce pour la culasse
magntique. Les inconnues tant les densits de courant dans chaque lment de
discrtisations des enroulements et du noyau magntique. [Ern-1987]
Llment ferromagntique, peut tre remplac par des densits de courant fictives
surfaciques et ventuellement volumiques. Considrons le circuit magntique (noyau) et une
distribution de courant vhicul par le conducteur caractrise par une densit de courant
Nous
supposons
que
le
milieu
est une conductivit lectrique
magntique
une
permabilit
0 . Le milieu extrieur au circuit magntique tant

tant la permabilit de vide.
de lair, il a une permabilit
La permabilit magntique du noyau est importante par rapport celle de lair, et le champ
magntique est perpendiculaire la surface du noyau [Leo-1992].
L'induction dans le circuit magntique se dcompose en un terme source
par la densit de courant source
gnr
cre par lenroulement primaire, et un terme de
magntisation
III. 32
Cest linduction magntique
rsultante qui agit sur le milieu pour produire lintensit
daimantation .
50
Chapitre III

Fig.III.8 : Schma de principe du systme bobine culasse [Doi-2007]

Le milieu magntique produisant linduction magntique
se comporte comme la
superposition de deux distributions de courant. Une distribution volumique fictive

distribution de courant superficiel fictive
et une
, circulant respectivement lintrieur de la
culasse et sur sa surface, avec :

III. 33
III. 34
tant la surface oriente positivement lextrieur en tout point de la surface. Lhypothse

dhomognit de la permabilit magntique de la culasse implique que
0, donc
est nulle, ainsi la culasse est reprsente uniquement par une distribution de courants
superficiels
sur sa surface [Ern-1987] [Ern-1990] [Joa-2004].
III.6.2 Expression intgrale caractristique du noyau

Le noyau de gomtrie axisymtrique est donc modlis par distribution linique de courants
fictifs. A partir des conditions aux limites linterface du noyau, lquation intgrale de
Fredholm de 2me espce en
, peut tre utilise pour calculer linduction magntique en
chaque segment de la surface du noyau tel que :[Ern-1987]
51
Chapitre III
1
2

.
1
1
segments de longueur , sur les quels les densits de
Nous discrtisons le noyau en

courants fictifs
III. 35
sont constantes. La superposition de lensemble des lments liniques du
noyau, et tenant compte de leffet des bobines permet dexprimer la relation magntique au
point du noyau comme suit : III. 36
1
2
1
1
Segment du noyau
Magntisation propre
Segment du noyau
Magntisation mutuelle
III. 36
Segment de linducteur
Source
Pour le calcul de III. 36 nous devons dterminer linduction magntique cre par une spire
circulaire en tout point de la discrtisation
du noyau. Pour ce faire, nous utilisons
lexpression du potentiel vecteur magntique cr par une spire.

Ainsi, les composantes de linduction magntique lmentaire
et
drives de
lquation III. 21 scrivent :

. .
4 . .
III. 37
0 III. 38
.
4 .
III. 39
Avec
52
Chapitre III

Noyau ferromagntique
Axe de symtrie
Enroulement BT
Secondaire
Enroulement HT
Primaire
Rayon
_
_
_
_
Fig. III.9: Discrtisation du systme bobine et noyau [Leo-1992] [Doi-2007]

Finalement lquation de la tension aux bornes dune spire lmentaire est complte par
lquation de Fredholm en incluant la contribution du noyau ferromagntique. Cela augmente
le nombre dinconnues, ainsi alourdissement de la matrice globale sachant que le cot en
place mmoire dpend du nombre de discrtisations.
III. 40
Ainsi, le champ lectromagntique dans le circuit magntique peut tre calcul partir des
quations intgrales de Fredholm de deuxime espce.
53
Chapitre III

Le systme matriciel final est donn comme suit :

0
0
.
.
.
.
.
.
.
.
.
.
.
.
.
..
.
.
.
.
.
.
..
.
.
..
0
0
III. 41
.
.
.
.
.
.
.
..
0
.
.
0
;
;
sont des matrices pleines, elles proviennent de lquation de Fredholm de
deuxime espce III. 36 .
, elle reprsente linfluence entre les spires du noyau.
: Matrice de dimension
* , elle reprsente linteraction entre spires du

est la matrice transpose de mme dimension que la
primaire et celles du noyau ;
matrice
.
* , elle reprsente linteraction entre spires du

secondaire et celles du noyau ;
matrice
est la matrice transpose de mme dimension que la
III.7 Calcul analytique des paramtres inductifs du modle

Plusieurs techniques sont proposes pour le calcul des paramtres du modle
savoir les rsistances et les inductances. On peut classer celles-ci en trois grandes classes
la mthode base sur la gomtrie et les proprits des matriaux ; mthode dnergie
emmagasine dans le systme et la mthode base sur une formulation en potentiel vecteur
magntique.
54
Chapitre III

III.7.1 Inductance mutuelle
Daprs le thorme de Stockes le flux dinduction magntique est donn par lquation
suivante : [Bab-2002] [Cha-2003] [Ven-2009]
.
III. 42
En utilisant lexpression analytique du potentiel vecteur, on dduit les expressions des

mutuelles
:
.
III. 43
L'inductance mutuelle est donne par l'expression I. 29 : [Ven-2009]

On pose :
2
III. 44
III. 45
Lexpression de la mutuelle se simplifie comme suit :
III. 46
Dans un systme (solnode massif) nous subdivisons lenroulement en mailles filamentaires,

la section de l'enroulement primaire est divise en
secondaire en
cellules, et celui de l'enroulement
cellules, chaque cellule contient un filament, et on suppose que la densit de
courant dans l'enroulement est uniforme, de sorte que tous les courants du filament soient
gaux. L'inductance mutuelle de cette configuration est obtenue par la superposition des
inductances mutuelles de chaque couple compos par les enroulements (spires filamentaires)
[Bab-2004] [Bon-1997] [Bab-2003].
55
Chapitre III

.
.
III. 47
: Nombre de spires de lenroulement primaire.

est le nombre de spires de lenroulement secondaire.
est La mutuelle entre deux lments lmentaires et .
III.7.2 Inductance propre

Le calcul de linductance propre dune spire circulaire ne peut tre entrepris sans dfinir
certains paramtres gomtriques qui sont : le rayon moyen de la spire ; le rayon interne de
la spire
; le diamtre du fil constituant la spire et le rapport k dfinis comme k= .
Fig. III.10 : Dfinition gomtrique dune spire lmentaire

Linductance propre, est la somme de deux termes. Un premier interne , que nous
est associ au flux traversant la section du conducteur, un second
dnommerons
externe , traduisant le flux partiel embrass par le fil. [Fou-1979]

III. 48
Linductance interne dune bobine est associe la distribution de la densit dnergie interne
aux fils constituant la bobine.
2
III. 49
est lnergie stocke dans le volume
des conducteurs, et
56
le courant efficace total.
Chapitre III

En accord avec la loi dAmpre, applique des conducteurs de section circulaire, le champ
dinduction magntique dans les conducteurs est donn par :
III. 50
Avec III. 49 et III. 50 ,
devient :
III. 51
: est la longueur totale de la spire.

On peut assimiler linductance propre extrieure la mutuelle entre la fibre moyenne
et une fibre priphrique quelconque [Fou-1979], do :
ln
2 III. 52
Lexpression de linductance propre est donne par

ln
7
III. 53
4
III.8 Ncessit du couplage avec la mthode des lments finis
La mthode des circuits lectriques magntiquement coupls qui ne discrtise que les parties
actives (inducteur, charge amagntique, etc.), peut tre avantageusement couple une
mthode numrique de reprsentation des milieux non linaires. Le noyau ferromagntique
peut tre modlis par une mthode analytique base sur l'introduction de densits fictives et
superficielles de courant pour un cas linaire isotrope est homogne, vue les caractristiques
physiques non linaires du noyau magntique la ncessit dun couplage avec une mthode
numrique est invitable. Le couplage tire profit des avantages des deux mthodes, cela nous
offre lavantage de limiter le domaine de rsolution de ces mthodes aux seuls domaines
conducteurs. Cela engendre la rduction de maillage (non maillage de lair) et un gain en
espace mmoire et en temps de rsolution.
57
Chapitre III

III.8.1 Technique mixte circuits coupls-lments finis

On applique la mthode des circuits coupls uniquement dans les enroulements primaire et
secondaire et sur la frontire du domaine, et on lui associe la mthode des lments finis
lintrieur du domaine
Fig. III.1
Noyau
Enroulement BT
Enroulement HT
Fig. III.1 : Discrtisation des enroulements et du noyau

III.8.2 Formulation lments finis des quations lectromagntiques du noyau
Le domaine dtude est dcompos en un nombre fini dlments triangulaires qui forment
le maillage. La valeur du potentiel vecteur est dtermine sur tous les sommets (les nuds)
des triangles. En utilisant des fonctions dinterpolation appropries, la solution sera
dtermine sur tout point sommet de llment.
Pour transformer un systme dquations aux drives partielles en une formulation intgrale,
on utilise soit la formulation des rsidus pondrs ou encore la mthode variationnelle. On se
limitera dans ce qui suit prsenter la mthode projective de Galerkine [Sab-1986] [Dha2005].
58
Chapitre III

En remplaant le terme
par dans lquation magntodynamique donne par II. 43 , et
en annulant le terme source on aura :
0 III. 54
. . .
La formulation rsidus pondrs relative aux problmes magntodynamiques harmoniques

exprims par lquation prcdente sera :
.
. . .
0 III. 55
Lapplication du thorme de Green nous permet dcrire

.
. . .
0 III. 56
Utilisant la procdure de Galerkine qui consiste choisir des fonctions dinterpolation

identiques aux fonctions de pondration,
est approxim par : [Per-1995]
Avec :
: Nombre de nuds du domaine de rsolution.
: Inconnues en chaque nud du domaine.
Lorsquon introduit lexpression de discrtisation de linconnue dans la formulation lments
finis cette dernire devient :
.
. .
.
.
0 III. 57
59
Chapitre III

En crivant III. 57 pour tout nud du domaine on aboutit au systme matriciel suivant
.
0 III. 58
Les coefficients des diffrentes matrices sont fournis par :

.
III. 59
: Fonction de forme associe au nud du domaine.

.
. .
III. 60
Le couplage hybride CEMC et EF se fera travers les valeurs du potentiel vecteur

magntique
pos comme tant conditions aux limites de type Dirichlet.
On obtient un systme de
quations qui scrit :
III. 61
Avec :
III. 62
: correspond aux valeurs du potentiel sur les frontires du noyau.

A partir
III. 24 , III. 25 , III. 26 , III. 58 , III. 61 On aura la matrice globale, les
inconnues sont les termes du potentiel vecteur magntique sur la frontire et lintrieur du
noyau (milieu).
60
Chapitre III
0 0

0
.
0
.
.
. III. 63
0
0
.
.
0
: Matrice pleine, qui caractrise leffet du primaire sur la frontire du noyau.

: Matrice pleine, qui caractrise leffet du primaire sur le noyau.
: Matrice pleine, qui caractrise leffet du secondaire sur la frontire du noyau

: Matrice pleine, qui caractrise leffet du secondaire sur le noyau
: Matrice creuse, qui caractrise leffet de potentiel vecteur sur la frontire du noyau
: Matrice creuse, qui caractrise leffet de potentiel vecteur magntique dans le noyau
III.8.3 Mthodes de rsolution des systmes algbriques

Gnralement on peut classer les mthodes de rsolution en deux catgories : les mthodes
itratives et les mthodes directes. Les matrices qui proviennent de la mthode des lments
finis sont creuses par contre les matrices qui proviennent de la mthode des circuits coupls
sont pleines.
61
Chapitre III

Il existe des cas ou les structures du systme linaire ne peuvent tres rsolues par les
mthodes directes. C'est par exemple le cas des systmes ou la matrice
est trs creuse. C'est
la raison pour laquelle, on prfre utiliser des mthodes itratives.

Les mthodes itratives donnent, en thorie, la solution dun systme linaire aprs
un nombre fini ditrations. A chaque pas, elles ncessitent le calcul du rsidu du systme.
Le modle algbrique obtenu sera rsolu par un processus itratif, on suppose a priori
la valeur initiale de linconnu (fournie par les conditions initiales).
La solution une itration donne est obtenue partir de la solution prcdente, on arrte les
itrations lorsque la diffrence relative entre deux solutions successives
et
est
infrieure une valeur donne tel que :
Le modle dquations algbriques scrit :
III. 64
III. 65
III. 66
III. 67
Le systme dquations devient :
III. 68
III. 69
III. 70
III. 71
62
Chapitre III

Nous utilisons lalgorithme (Fig. III.2) pour la rsolution de notre problme
Initialisation
et
Calcul de
de lquation III. 68
Connaissant
et
on calcul de
de lquation III. 69
=
=
Connaissant
et
on calcul
de lquation III. 70
Connaissant
et
on calcul
de lquation III. 71
Non
Teste de convergence
Oui
Affiche des rsultats et calcul
des pertes Joule et pertes fer
Fig. III.2 : Algorithme de rsolution

63
Chapitre III

III.9 Conclusion
Nous avons dtaill dans ce chapitre la mthode des circuits lectriques magntiquement
coupls, applique un dispositif axisymtrique reprsent par un transformateur. Les
expressions introduites permettent de calculer les inductances et les rsistances et dterminer
la rpartition des courants inducteurs, cela en appliquant le modle excitation en tension .
La modlisation du circuit magntique est faite analytiquement par lutilisation de densits
superficielles et fictives de courant, ce dernier se comporte comme un matriau linaire,
isotrope et homogne utiles pour le fonctionnement vide, en court circuit et au voisinage
dun point en charge en supposant que sa permabilit magntique est constante.
Lapplication de la mthode des circuits coupls nous permettra de dterminer les valeurs de
la densit de courant dans chaque spire. Cette mthode semi-analytique est limite aux
dispositifs axisymtriques de milieux linaires. Son couplage ventuel avec une mthode
numrique permet de considrer la non linarit, pour cette raison nous avons effectu un
couplage hybride de deux mthodes analytique et numrique. Un modle 2D lments finis en
termes de potentiel vecteur magntique pour le calcul du champ magntique dans la partie
ferromagntique (non linaire), combin une mthode semi-analytique (Circuit Electriques
Magntiquement Coupls) dans les enroulements primaire et secondaire.
Finalement ce modle de couplage est appliqu deux modles de transformateurs de
puissance pour valider lefficacit de la mthode propose.
64
Chapitre IV
Applications et validations
Chapitre IV
IV.1 Introduction
Le but de ce chapitre est de valider les codes de calcul implment travers la simulation du
fonctionnement en rgime statique proprits physiques linaires du transformateur. Les
modles circuits lectriques magntiquement coupls et lments finis mis en uvre ainsi que
la
procdure
de
leur
couplage,
sont
implments
sous
environnement
Matlab.
Limplmentation des codes de calcul porte sur la discrtisation-maillage, calcul des diverses
matrices inductances et rsistance, des intgrants lments finis, et de lalgorithme de
couplage hybride (CEMC-EF).
En premier lieu nous dfinissons les caractristiques gomtriques et physiques du
transformateur dtude. En second lieu, nous prsenterons les simulations ralises en rgime
vide et en court-circuit en vue de connaitre dune part les niveaux de valeurs des diverses
grandeurs courants, potentiel vecteur magntique, induction magntique, paramtres rsistifs
et inductifs, et dautres part dterminer les pertes joules et les pertes fers respectivement dans
les enroulements primaire et secondaire et dans le circuit magntique.
Les rsultats de simulation du comportement lectromagntique relatifs au deux modles
semi-analytique (circuits lectriques magntiquement coupls) et hybrides CEMC-EF, aussi
bien en rgime de fonctionnement vide quen court-circuit seront confronts aux mesures
exprimentales fournies par lentreprise Electro-industries (ex- ENEL Azzazga, Tizi-Ouzou).
IV.2 Description du transformateur
Le transformateur de distribution utilis est fabriqu et commercialis par lentreprise Electroindustries [E.E.I-2008]. Il sagit dun transformateur triphas de distribution trois colonnes
refroidissement par huile. Les enroulements haute tension (HT) sont connects en triangle et
comportent chacun 3000 spires rparties sur 14 couches, et les enroulements basse tension
(BT) connects en toile comportent chacun 22 spires rparties sur deux couches.
65
Chapitre IV
Caractristiques
Valeur
Unit
Puissance apparente (Sn)
630
kVA
Frquence
50
Hz
Tension nominale de lenroulement HT (Un)
30
kV
Tension nominale de lenroulement BT (Un)
0.4
kV
Groupe de couplage
Dyn 11
IV.4 Dimensions gomtriques du transformateur

Composante
Unit
198
mm
Diamtre intrieur de lenroulement BT
Diamtre extrieur de lenroulement BT
239
mm
Diamtre intrieur de lenroulement HT
265
mm
356
mm
hauteur des enroulements
436
mm
Epaisseur de lenroulement BT
20.5
mm
43
mm
180
mm
hauteur du noyau
810
mm
Diamtre du fil nu
15
mm
Diamtre du fil avec isolation
16
mm
Epaisseur de lenroulement HT
Diamtre du noyau
Section de fil
7.5
66
valeur
mm
Chapitre IV

_
_
Fig. IV.1 : Dimensionnement du transformateur
IV.5 caractristiques physiques
Milieu
Caractristiques
Physiques
Air
Enroulement
Noyau
Conductivit lectrique
0,
10
10
4
5.9 10
67
Permabilit magntique
1,
10 ,
1,
.
=1
Chapitre IV
IV.6 Discrtisation de la structure

Nous considrons le transformateur dans la condition quilibre triphase. Par consquent
nous pouvons seulement calculer une phase. Cependant, le calcul peut tre prolong pour
couvrir les trois phases. [Liu-1996]
Zoom
Enroulement BT
Enroulement HT
Noyau magntique
Enroulement hlicodale
Fig. IV.2 : Modle dune phase utilis en simulation [Jan-2009]
Lorsque la permabilit magntique relative est gale lunit
1 , toute la surface du
circuit magntique peut tre discrtise tel que donn par la figure (IV.3).
68
Chapitre IV
400
300
200
100
Zoom
0
-100
-200
-300
-400
100
200
Fig. IV.3 : Discrtisation dune colonne du transformateur
Lorsque la permabilit magntique est diffrente de lunit, la mthode des circuits

lectriques magntiquement coupls est applique a travers la discrtisation des enroulements
et de la frontire du noyau magntique. Les enroulements sont discrtiss par des spires
lmentaires circulaires, et la surface du noyau par 400 segments. La densit de courant est
calcule au centre de chaque lment de discrtisation linique du noyau ferromagntique du
transformateur.
69
Chapitre IV
400
300
10
200
2 10
100
-100
5.9 10
-200
4 10
-300
1
0
-400
0
100
200
Fig. IV.4 : Discrtisation dune colonne du transformateur par la CEMC

La seconde approche plus gnrale est issue du couplage circuits lectriques magntiquement
coupls dans les enroulements et lutilisation de la mthode des lments finis dans le noyau.
400
120
300
110
100
90
200
80
70
60
100
100
110
120
130
140
150
160
170
0
360
340
-100
320
300
280
-200
260
-20
20
40
60
80
100
120
-300
-400
-400
-300
-200
-100
100
200
300
400
500
600
Fig. V.5 : Maillage mixte circuits lectriques magntiquement coupls et

lments finis (CEMC- EF)
70
Chapitre IV
IV.7 Calcul des pertes

Cette partie rpertorie les pertes qui ont pour origines des phnomnes lectromagntiques
dans le transformateur, et appliqu les modles mis en uvre en vue de calculer ces dernires.
Notre intrt portera sur la dtermination des pertes joules gnres dans les enroulements est
des pertes fer dans le circuit magntique.
Pour un rgime magntodynamique une frquence donne du rseau, et des valeurs
constantes des paramtres lectriques des enroulements, les pertes joules dpendent
principalement des courants eux-mmes lis aux rgimes de fonctionnement ( vide, en courtcircuit, et en charge). Les pertes fer frquence et permabilit magntique donnes sont
fixes par les valeurs des inductions magntiques dans le noyau. Les nivaux dinductions
magntiques sont implicitement lis aux courants dans les enroulements.
IV.7 .1 Processus de calcul des pertes totales
Le calcul des pertes totales se fait selon les tapes suivantes :
1) Calcul des courants primaire et secondaire
2) Calcul de la rsistance de chaque enroulement
3) Estimation des pertes Joule.
4) Calcul des selfs inductances
des enroulements et des mutuelles
entre les
conducteurs.
5) Calcul de linduction magntique axiale
et radiale
de chaque conducteur des deux
enroulements en utilisant les courants primaire et secondaire calculs.

6) Calcul des pertes par courant de Foucault.
7) Estimation des pertes totales.
IV.7.2 Pertes en court-circuit
cest la tension rduite quil faut appliquer au primaire lorsque le secondaire est en court
circuit pour que le courant nominal circule au primaire. Sa valeur est comprise entre 4 et 5%
de la tension nominale.
La puissance absorbe lors du test en court circuit correspond aux pertes Joule dues au
courant de charge transitant dans les enroulements, donne par lexpression suivante:
IV. 1
71
Chapitre IV
sont respectivement la rsistance quivalente de lenroulement primaire
et lenroulement secondaire, calcules comme suit :

.
: est la partie diagonale des matrices
et
et
Les pertes lies la frquence dans les conducteurs se produisent comme courants
de Foucault dus aux champs magntiques variables dans le temps. Les courants de
Foucault causent une augmentation de pertes et une rduction de la quantit nette de flux
magntique.
Les pertes totales en court-circuit en (Watt) sont donnes par la formule suivante
IV. 2
Tension de court-circuit
0V;
nulle
= 5 % de la tension nominale, et la tension au secondaire est
2000.
Essai en court circuit
Courant primaire
Courant secondaire
Modle
Modle
Mesures
Modle
Modle
Mesures
CEMC
CEMC-EF
Electro-
CEMC
CEMC-EF
Electro-
Industries
= 9.646 A
= 9.299 A
=10.15 A
Industries
= 947.53 A
=930.79A
= 909 A
IV.7.2.1 Calcul des rsistances et selfs

Les inductances
sont calcules comme suit :
: Partie imaginaire de la matrice.

rsistance du primaire
Modle CEMC
Mesures ENEL
17.22
17.33
self du primaire
0.5
Inductance mutuelle (primaire-secondaire)
4.49mH
rsistance du secondaire
Modle CEMC
Mesures ENEL
1.35
1.42
self du secondaire
2.1
72
Chapitre IV
IV.7.2.2 Calcul des mutuelles et selfs de chaque conducteur

Le calcul des mutuelles entres cinq points de lenroulement primaire et secondaire est donn
dans le tableau suivant :
Inductance
pri/sec
= 0.211
=0.210
=0.208
=0.206
0.204
Inductance
pri/pri
=0.126
=0.710
=0.594
=0.526
0.478
Inductance
sec/sec
= 0.097
=0.221
=0.139
=0.961
0.693
IV.7.2.3 Calcul des inductions

La formule pour les pertes
et
par courant de Foucault peut tre donne par: [Shi-2002]
[Lin-2009] [Jan-2009]:
Pour un seul conducteur de rayon

.
IV. 3
Les pertes pour
conducteurs sont donnes par : [Shi-2002]
IV. 4
32
Fig. IV.6 : Effet de peau dans un conducteur rond [Shi-2002]

73
Chapitre IV
Enroulement
Enroulement primaire
Enroulement secondaire
Point (m)
A (1,2)
0 .952
0.124
0.090
0.099
B (1,3)
0.475
0.031
0.041
0.026
C (1,4)
0.316
0.014
0.025
0.012
2
32
2
32
IV. 5
IV. 6
IV. 7
Les pertes totales en court circuit sont donnes par le tableau suivant :
Pertes Joule
Pertes par courant de Foucault
8.441
Pertes en court-circuit
1.104
9.545
Pertes Joule totales au primaire et au secondaire

Modle CEMC
Modle CEMC-EF
Mesures
Electro-Industries
8.441
8.084
7.75
Effet de peau nglig dans

les conducteurs
IV.7.3 Calcul des pertes vide

Dans ce cas on suppose que le matriau magntique a un comportement linaire caractris
par une permabilit magntique relative constante calcule dans la partie linaire de la
courbe de premire aimantation. Dans ces conditions la permabilit relative vaut
74
2000.
Chapitre IV
Les pertes vide
se composent des pertes par courant de Foucault et les pertes par

est petit par rapport sa valeur en pleine charge,
hystrsis. Comme le courant vide
les pertes Joule sont ngliges dans les enroulements, ainsi :

IV. 8
Ces pertes vide sont appeles aussi pertes fer qui sont conditionnes par la valeur de
linduction la surface de la tle magntique et la gamme de frquence utilise.
Lexpression IV. 8 Peut scrire :
.
Les pertes dilectriques (
IV. 9
) sont pour la plupart du temps ngligeables par
rapport aux autres pertes.

.
Lexpression IV. 9 presente ce quon appelle une dcomposition des pertes. La premire
intgrale reprsente les pertes Joule (pertes par courant de Foucault) et la seconde les pertes
magntiques.
Les pertes par courant de Foucault sont calcules au moyen de : [Liu-1996]
1,6. .
1.77
IV. 10
(Dtermine daprs Fig.IV.7)
: La masse de la tle magntique
: Valeur Max de linduction magntique.

: Densit de masse
: paisseur dune feuille de la tle magntique
75
Chapitre IV
2
Inductionmagntique aux points de discritisation (T)
1.8
Induction magntique (T)
1.6
1.4
1.2
0.8
0.6
1.8
1.6
1.4
1.2
0.8
0.6
0.4
0.4
-0.5
0
Diretion axiale (m)
0.5
0.02
0.04
0.06
0.08
Direction radiale (m) Nc=400
Fig. IV.7 : Induction dans le circuit magntique vide
La rsistance quivalente correspondant aux pertes par courant de Foucault dans le

noyau est donne par la formule suivante :
IV. 11
4.44
IV. 12
est la fem induite et
le flux magntique maximal dans le noyau.
Daprs les quations IV. 10 et IV. 11 , IV. 12 on peut crire :

4.44
1.6
861.8 10 IV. 13
: Section du noyau
76
Chapitre IV
Les pertes dans le noyau magntique sont proportionnelle au carre de la

frquence IV. 10 et au carre de la tension induite au primaire. La rsistance
qui
prsente ces pertes est constante. [Ell-2003]

Les pertes totales dans le circuit magntique peuvent tres approches par : [Ell-2003]
IV. 14
dcrivent respectivement le noyau et la culasse

,
sont la masse du noyau et de la culasse.
: Facteur de construction.
: 7.48 10
5.2 10
1.79.
La tension au primaire
400.
Essai vide
Courants primaire vide
Modle CEMC
Modle CEMC-EF
Mesures
Electro-Industries
11.91
12.15
12.85
Pertes vide
Modle CEMC
Modle CEMC-EF
Mesures
Electro-Industries
1.433
1.316
77
1.365
Chapitre IV

6
x 10
5.3581
6
Densit de courant primaire (A/m2)
Densit de courant du secondaire en court circuit (A/m)
5.3581
x 10
5.3581
5.3581
5.358
5.358
5.358
5.358
5.358
5.3579
0
0.13
5.3579
-1
-0.5
0
0.5
1
1.5
Longueur de la spire cot secondaire (m)
Fig. IV.8 : Evolution de la densit de courant au

secondaire en court circuit
0.14
0.15
0.16
0.17
0.18
Longueur de la spire cot primaire (m)
Fig. IV.9 : Evolution de la densit de courant

au primaire en court circuit
Essai vide
6
90
Dens it de c ourant au prim aire v ide (A /m 2)
D ens it de c ourant c ot s ec ondaire v ide (A /m 2)
100
80
70
60
50
40
30
20
10
0
-1
-0.5
0.5
0
0.13
1.5
Longueur de la spire cot secondaire (m)
0.14
0.15
0.16
0.17
Longueur de la spire cot primaire (m)
Fig. IV.10 : Evolution de la densit de courant

au secondaire vide
78
x 10
Fig. IV.11 : Evolution de la densit de

courant au secondaire vide
Chapitre IV
Induction magntique radiale aux points de discrtisations ( T )
1.8
Induction magntique (T)
1.6
1.4
1.2
0.8
0.6
0.4
-0.5
0
Diretion axiale (m)
0.5
-3
x 10
4
3
2
1
0
-1
-2
-3
-4
-5
0.05
0.1
Distance de point de discritisation

au centre du noyau (m)
Fig. IV.12 : Allure de linduction magntique axiale et

radiale sur la frontire vide
La rpartition du potentiel vecteur magntique dans le noyau pour diffrentes hauteurs selon
(oz) et suivant la frontire sont illustres par les figures (IV.13), (IV.14), (IV.15), (IV.16).
Evolution du potentiel vecteur magntique vide
-5
Module du potential vecteur magntique [T.m]
x 10
Hauteur z=370 mm
Hauteur z=0 mm
Hauteur z=-370 mm
8
7
6
5
4
3
2
1
0
0
10
20
30
40
50
60
Direction radiale en [mm]
70
80
90
Fig. IV. 13 : Evolutions radiales du module du potentiel vecteur magntique

dans le noyau pour diffrentes hauteurs axiales
79
Chapitre IV
-4
1.6
x 10
Abscisse r=90 mm
1.4
1.2
1
0.8
0.6
0.4
0.2
0
-500
-400
-300
-200
-100
0
100
200
Direction axiale en [mm]
300
400
500
Fig. IV. 14 : Evolution axiale du module du potentiel vecteur magntique

sur la frontire du noyau
Evolution du potentiel vecteur magntique en court circuit
-5
x 10
Hauteur z=370 mm
Hauteur z=0 mm
Hauteur z=-370 mm
1.8
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0
10
20
30
40
50
60
Direction radiale en [mm]
70
80
90
Fig. IV.15 : Evolutions radiales du module du potentiel vecteur magntique

dans le noyau pour diffrentes hauteurs axiales
80
Chapitre IV
-4
x 10
Abscisse r=90 mm
0
-500
-400
-300
-200
-100
0
100
200
Direction axiale en [mm]
300
400
500
Fig. IV. 16 : Evolution axiale du module du potentiel vecteur magntique

sur la frontire du noyau
IV.8 Interprtation des rsultats
Les rsultats pour la premire application qui concerne la discrtisation du noyau pour une
permabilit relative
correspond une forte saturation du noyau, le circuit
magntique ne canalise pas le champ magntique et le flux de fuite devient important ce qui
engendre une augmentation des pertes dans le circuit magntique.
Cela nous a conduis choisir une autre valeur de permabilit diffrente de lunit pour le
circuit magntique.
Ltude des grandeurs globales que sont la rsistance, et courant montre la similarit entre la
mthode CEMC et le calcul exprimental.
Pour lessai en court circuit linductance de fuite totale est trs majoritairement associe
l'enroulement primaire (quelques
contre quelques
pour l'enroulement secondaire).
Les flux de fuite se referment dans l'air et la prise en compte du noyau a alors tendance faire
diminuer l'inductance de l'enroulement extrieur tout en augmentant celle de l'enroulement
intrieur. Linductance totale diminue pour ce transformateur en prsence du circuit
ferromagntique.
81
Chapitre IV
Pour juger la fiabilit des modles axisymtriques une confrontation avec des mesures
fournies par Electro-Industries, savre ncessaire afin darriv une valuation prcise des
performances des modles ainsi proposs.
Lcart entre les valeurs des courants obtenue par le calcul axisymtrique et celle donnes par
Electro-Industries induit une erreur aux primaire de 5% avec la mthode des circuits
lectriques magntiquement coupls et 8.3 avec le couplage hybride, et une erreur au
secondaire de 4% avec la CEMC et 2,3 avec le couplage hybride. Des erreurs relatives
gnralement acceptables.
Les pertes par courant de Foucault sajoutent aux pertes Joule donnant ainsi les pertes totales
en court- circuit dune valeur de 9.54 KW.
En se basant sur les valeurs de rsistance et des courants calculs, nous aurons 8.441KW de
pertes Joule avec la CEMC et de 8.084KW avec le couplage hybride. Comparant la valeur
donne par Electro-industries des erreurs relatives associes sont respectivement de 8.2% et
4.1%. Le modle hybride donne pratiquement de meilleurs rsultats.
En essai vide, Lcart entre les valeurs n'est que denviron 5% avec le couplage hybride et
7% avec la CEMC cela entrane une faible erreur dans le calcul. Par consquent les pertes
vide sont lgrement dfrentes de celles mesures par Electro-industries. Lerreur relative
commise est inferieure la valeur admissible qui est de 15%.
Les figures (IV.8, IV.9) donnent la rpartition des densits de courant dans lenroulement
primaire et secondaire en court circuit.
La densit de courant
est
uniforme dans lenroulement secondaire cela revient la
profondeur de pntration des courants induit qui est grande par rapport au rayon du
conducteur. Effet de peau nglig.
Par contre dans lenroulement secondaire la densit du courant est non uniforme cela
sexplique par leffet de proximit entre spires et entre couches. L'effet de proximit
caractrisant l'impact des effets induits est troitement corrl au nombre de spires
conductrices et le nombre de couches d'enroulement composant le bobinage complet.
Les figures (IV.10, IV.11) donnent la rpartition des densits de courant dans
lenroulement primaire vide. Pratiquement on trouve un courant nul dans le secondaire.
La prsence du noyau ferromagntique modifie cependant le trajet de la rpartition des
densits de courant et par consquent vient dajouter un effet mutuel sur les conducteurs.
La figure (IV.12) donne la rpartition des inductions magntiques radiale et axiale sur le
contour du noyau. La valeur maximale du module de cette dernire est de 1,77 Tesla.
82
Chapitre IV
La prcision de la solution du modle circuit coupl dpend grandement du nombre de

dcoupage liniques en segments du circuit magntique.
Les figures (IV.14-16) reprsentent la rpartition du module de potentiel vecteur magntique
sur la frontire dcoulant du modle coupl (EF-CEMC).
Le potentiel vecteur magntique prend la valeur max la hauteur axiale z=0, de 4,9
10 Tmen court circuit et dune valeur de1.5 10 Tm vide et diminue puis sannule
aux extrmits du noyau.
A lintrieur du noyau le potentiel vecteur magntique diminue en augmentant la hauteur
axiale (z) et en sloignant de lintervalle entre noyau et enroulements. Cela revient la haute
attraction de flux de fuite d la permabilit magntique du circuit magntique. et prend la
valeur max 1.9 10 Tm en court circuit et de 8 10 Tm vide, figures (IV.13-15).
IV.9 Description du second transformateur

Transformateur
triphas
trois
colonnes
de
puissance
de315kVA,11.5/0.4KV.
Les enroulements (HT) comportent 1882 spires, les enroulements (BT) comportant 36 spires
une seule couche [Bah-1999].
Caractristiques
Valeur
Unit
Puissance apparente (Sn)
315
kVA
Frquence
50
Hz
Tension nominale de lenroulement HT (Un)
11.5
kV
Tension nominale de lenroulement BT (Un)
0.4
kV
Nombre de spires de lenroulement BT
36 (une couche)
IV.11 Dimensions gomtriques du second transformateur
Composante
Unit
184
mm
Diamtre intrieur de lenroulement BT
Diamtre extrieur de lenroulement BT
220
mm
242
mm
83
Valeur
Chapitre IV
300
mm
580
mm
590
mm
18
mm
29
mm
Diamtre du noyau
160
mm
hauteur du noyau
630
mm
Section du conducteur
227.5
mm
Epaisseur de lenroulement BT
Epaisseur de lenroulement HT
Section du conducteur
1.2
Poids du noyau
294
Kg
Poids de la culasse
280
Kg
300
200
100
-100
-200
-300
0
100
200
Fig. IV.18 : Discrtisation du dispositif
84
mm
Chapitre IV
IV.12 Essai en court-circuit

Tension de court-circuit
= 4 % de la tension nominal
Courant primaire
Modle
Modle
CEMC
CEMC-EF
= 9.59 A
=10.37 A
Courant secondaire
Calcul
Modle
Modle
CEMC
CEMC-EF
= 492.38 A
=443.37 A
Calcul
[Bah-1999]
= 455 A
Pertes Joule totales au primaire et au secondaire

Modle CEMC
Modle CEMC-EF
3.30
Les selfs
Calcul [Bah-1999]
3.176
et linductance mutuelle
3.45
sont calculs comme suit :
rsistance du primaire
rsistance du secondaire
5.167
1.7
self du primaire
self du secondaire
0.27
3.2
Inductance mutuelle primaire-secondaire
3.319
Le calcul des mutuelles est reprsent par le tableau suivant :
Inductance
pri/sec
= 0. 197
= 0.196
= 0.193
= 0.189
0.184
Inductance
pri/pri
= 0.115
= 0.573
= 0.468
= 0.407
0.363
Inductance
sec/sec
= 0.097
= 0.217
= 0.140
= 0.100
0.089
85
Chapitre IV
IV.13 Essai vide
La tension au primaire
la tension au secondaire
11500;
la permabilit magntique
250 ;
400
Essai vide
Courants primaire vide
Modle CEMC
Modle CEMC-EF
8.99
Calcul [Bah-1999]
8.354
Pertes vide
Modle CEMC
Modle CEMC-EF
1.06
0.735
Calcul [Bah-1999]
0.66
IV.14 Conclusion
Dans ce chapitre nous avons procd aux validations de lapproche semi-analytique pour un
circuit magntique isotrope et linaire, et une approche de couplage hybride pour un circuit
non linaire.
La rsolution de lalgorithme de couplage implment sous environnement Matlab nous
conduit lestimation des diffrentes valeurs des densits de courant sources et induits dans
les enroulements ainsi que le potentiel vecteur magntique dans le circuit magntique et sur la
frontire de ce dernier.
Les diffrentes pertes vide et les pertes en court circuit dans les parties actives du
transformateur sont en concordance avec celles donnes par Electro-Industrie.
La mthode hybride de modlisation fournit pratiquement de meilleurs rsultats pour
lanalyse du comportement du transformateur et peut tre effectivement utilise pour prvoir
les performances du transformateur .
La validation de deux mthodes dveloppes est effectue sur deux modles de
transformateurs de puissances 630 et 315 KVA.
86
Conclusion gnrale
Ce travail nous a permit de faire une tude sur le comportement lectromagntique
dun transformateur de puissance en rgime statique, et dvaluer les paramtres lectriques
de ce dernier aussi bien que lestimation des pertes fer et les pertes dans le cuivre travers un
modle semi-analytique (CEMC) et un modle de couplage hybride (CEMC-EF).
Nous nous sommes intresss en premier lieu une approche semi-analytique
(CEMC), La mthode des circuits lectriques magntiquement coupls a montr dexcellentes
performances modliser les structures lectromagntiques. La raison en est sa bonne
capacit prendre en compte les gomtries axisymtriques comprenant un grand nombre de
conducteurs et systmes de grande dimension, de plus sa facilit dutilisation et lestimation
des pertes par courant de Foucault dans les enroulements. Le nombre de discrtisation nest
limit que par la capacit mmoire des micro-ordinateurs.
Une confrontation des rsultats en terme de grandeurs globales (rsistances, courant)
obtenus par simulation et ceux issus des exprimentations fournis par lentreprise ElectroIndustrie vrifie la concordance des rsultats obtenus. Ce qui nous permet de conclure la
validit du code de calcul labor. Nanmoins, ces modlisations sont imparfaites et
prsentent des limitations.
La CEMC ne prend pas en compte la non linarit du circuit magntique ce qui veux
dire elle traite mal les matriaux magntiques et une partie des pertes totales est nglige. Les
valeurs du champ magntique sont obtenues la surface du noyau exactement en chaque
point de discrtisation. A cet effet nous avons labor un code de calcul hybride (CEMC-EF)
permettant partir des ces valeurs, d'tendre au circuit magntique complet et rel (maillage),
et un calcul prcis des pertes dans le circuit ferromagntique. Ce dernier est moins coteux en
temps de calcul et en place mmoire, tout en garantissant une meilleure prcision. Ce travail
ouvre des perspectives pour une modlisation fine du problme lectromagntique.
La validit de lapproche nous permettra de lappliquer pour ltude du rgime
transitoire et la variation des grandeurs lectromagntiques en fonction du temps suite une
surintensit interne ou externe, et galement tre prolong aux frquences suprieures
87
Bibliographie
Ouvrages
[Dur-1968] E. Durand Magntostatique Masson -1968.
[Fou-1979] G. Fournet Electromagntisme partir des quations locales Masson-1979
ISBN 2-7462-0548-3.
[Sab-1986] G.Sabonnadire J.Louis La mthode des lments finis-du modle la CAO
Edition Hermes, Paris, 1986.
[Sho-2004] T.A. Short Electric power distribution handbook Schenectady, New York,
2004 by CRC Press LLC, ISBN 0-8493-1791-6.
[Dha-2005] G.Dhatt D.Touzot E.Lefranois Mthode des lments finis Edition Hermes,
Lavoisier, Paris, 2005.
[May-2006] Pierre. May Aide mmoire- lectrotechnique Dunod,Paris, 2006, ISBN 210
049578 X.
Thses de doctorat et mmoires de magister
[Ahm-1992] Ahmad Ahmad Contribution a la modlisation des transformateurs de
puissance et de leur comportement en haute frquence Ecole centrale de Lyon, Thse de
Doctorat en Gnie Electrique 1992.
[Gue-1994] Christophe Guerin Dtermination des pertes par courants de Foucault
dans les cuves des transformateurs-modlisation des rgions minces et prise en compte
de la saturation des matriaux magntiques en rgime harmonique Thse de doctorat Institut
National Polytechnique de Grenoble, 1994.
[Che-1994] Abdellah Chentouf Contribution a la modlisation lectrique, magntique
et thermique dun applicateur de plasma inductif haute frquence Ecole Centrale de Nantes
IUT de Saint-Nazaire, Thse de doctorat en dynamique des fluides et des transferts, 1994.
[Mao-1996] Bachir Maouche Etude de dveloppement semi-analytique de lquation
de diffusion lectromagntique avec thermes de dplacement dans le cas de dispositifs
axisymtrique excit en courant ou en tension Mmoire de Magister en gnie lectrique
Universit de Bejaia, 1996.
[Liu-1996] Ye Liu Finite element modeling of the frequency characteristics of
transformers thesis, for the Degree of Master of Applied Science, Department of Electrical
and Computer Engineering, University of Toronto,1996.
[Le-1999]
Yvonnick Le Menach Contribution la modlisation numrique

tridimensionnelle des systmes lectrotechniques Thse de Doctorat de lUniversit de
Lille1, Anne 1999.
[Rob-1999] Frederic Robert Modlisation et simulation de transformateurs pour
alimentations dcoupage Thse de Doctorat en Sciences Appliques Universit Libre
de Bruxelles, Aot 1999.
[Dix-2004] Brandon Dixon, B.S.E.E Pulsed power conditioning with strongly coupled
transformers Thesis, Graduate Faculty of Texas Tech University, Master of Science
in Electrical Engineering May, 2004.
[Joa-2004] Michal Joan Modlisation des paramtres R et L de matriels lectriques
bobins par la mthode des lments finis 3D Institut National Polytechnique de Grenoble
Thse de Doctorat en Gnie Electrique 2004.
[Hen-2004] T.Henneron Contribution la prise en compte des grandeurs globales dans les
problmes dlectromagntisme rsolus avec la mthode des lments finis Thse de
Doctorat de lUniversit de Lille1, Anne 2004.
[Bje-2005] Eilert Bjerkan High frequency modeling of power transformers stresses and
diagnostics Doctoral Thesis for the degree of Doctor Ingenior, NTNU Norwegian
University of Science and technology, Trondheim, May 2005.
[Lef-2006] Anthony Lefevre Contribution a la modlisation lectrique, lectromagntique
et thermique des transformateurs application l'tude de l'chauffement sur charges
non linaires Thse de Doctorat de lUniversit de Nantes Anne 2006.
[Doi-2007] Vincent Doirat Contribution la modlisation de systmes de contrles
non destructifs par courants de Foucault-application la caractrisation physique
et dimensionnelle de matriaux de laronautique Thse de Doctorat de lUniversit
de Nantes Anne 2007.
Articles de revues
[Nak-1982] T. Nakata ; N. Takahashi ; Y. Kawase magnetic performance of step-lap joints
in distribution transformer cores IEEE Transactions On Magnetics, Vol. Mag-18, No. 6,
November 1982.
[Ern-1987] Roland.Ernst Annie.Gagnoud Isabelle.Leclercq Etude du comportement
dun circuit magntique dans un systme de chauffage par induction revue gnrale de
llectricit, N 9, (1987) 123-130.
[Pic-1989] L. Pichon M. Ayoub A.Razek Comparaison de techniques mixtes pour le calcul

de courants de Foucault en gomtrie axisymtrique Revue Phys Appl, 24 (1989)
pp 1049 -1056.
[Ern-1990] Roland.Ernst Annie.Gagnoud Isabelle.Leclercq Comportement dun induit
magntique dans un systme de chauffage par induction revue gnrale de llectricit, N 8,
(1990).
[Leo-1992] Francisco de Leon Adam Semlyen Efficient calculation of elementary
parameters of transformers Transaction on Power Delivery, Vol 7 No 1, pp 376-383
January 1992.
[Rez-1992] A.Rezzoug J.P.Caron F.M.Sargos Analytical calculations of flux induction
and forces of thick coils with finite length IEEE Transactions On Magnetics, Vol. 28, No. 5,
pp 2250-2252, September 1992.
[Per-1995] J.F.Pern S.N.Yeh Calculating the current distribution in power trasformer
windings using finite element analysis with circuit constraints IEE- Proc-Sci.Meas.Technol.
Vol 14, No 3, pp 231-236, May 1995.
[Che-1995] A.Chentouf J.Fouladgar G.Develey A simplified method for calculation
of the impedance of an induction plasma IEEE Transactions on Magnetic, Vol 31, No 3, pp
2100-2103, may 1995.
[Che-1996] A.Chentouf J.Fouladgar G.Develey Calculation of the impedance of an
induction plasma with a thick - walled cold crucible IEEE Transactions On Magnetics, Vol.
12, No. 3, pp 2030-2033, May 1996.
[Bon-1997] Ki.Bong.Kim Enrico.Levi Zivan.Zabar Leo.Birenbaum Mutual Inductance
of noncoaxial circular coils with constant current density IEEE Transactions
On Magnetics, Vol. 33, No. 5, pp 4303-4309, September 1997.
[Ger-2001] Herbert.De.Gersem Kay.Hameyer A Finite Element Model for Foil Winding
Simulation IEEE Transactions on Magnetics, Vol. 37, No. 5, pp 3427-3432, September
2001.
[Shi-2002] Wang Shishan, JI ShengChang, II Yanming The Study of Eddy Current Losses
in Coaxially Insulated Windings of Power Transformer International Conference on power
system technology , Vol.3 , No. pp 1392- 1395, 2002.
[Bab-2002] Slobodan Babic, Slaven Kincic New and fast procedures for calculating
the mutual inductance of coaxial circular coils (circular coildisk coil) IEEE Transactions
On Magnetics, Vol. 38, No. 5, pp 2367-3-2369, September 2002.
[Bab-2003] Slobodan Babic Cevdet Akyel Sheppard J.Salon New procedures

for calculating the mutual inductance of the system: filamentary circular coil-massive Circular
Solenoid IEEE Transactions On Magnetics, Vol. 39, No. 3, pp 1131-1134, May 2003.
[Cha-2003] H.L.Chan K.W.E.Cheng D.Sutanto Calculation of inductances of high
frequency air-core transformers with superconductor windings for DC-DC converters IEE
Proc.-Electr: Power Appl, Vol. 150, No. 4, pp 447-454, July 2003.
[Ell-2003] Mohammed Elleuch, Michel Poloujadoff, Analytical Model of Iron Losses in
Power Transformers IEEE Transactions On Magnetics, Vol. 39, No. 2, pp 973-980, Mar
2003.
[Rah-2003] Ebrahim.Rahimpour Jochen.Christian Kurt.Feser Hossein.Mohseni Transfer
Function method to diagnose axial displacement and radial deformation of transformer
windings IEEE Transactions On Power Delivery, Vol. 18, No. 2, pp 493-505, April 2003.
[Oli-2003] Juan Carlos Olivares; Yilu Liu; Jose M. Canedo; Rafael Escarela-Perez ; Johan
Driesen; Pablo Moreno Reducing Losses in Distribution Transformers IEEE Transactions
on Power Delivery, vol. 18, No. 3, pp 821-826, July 2003.
[Bab-2004] Slobodan Babic Sheppard Salon Cevdet Akyel The mutual inductance of two
thin coaxial disk coils in air IEEE Transactions On Magnetics, Vol. 40, No 2, pp 822-825
March 2004.
[Shi-2004] Yoshikazu.Shibuya Shigeto.Fujita High frequency model of transformer
winding Electrical Engineering in Japan, Vol.146, No, 3, 2004. Translate from Denki
Gakkai Ronbunshi, Vol 123-B, No 2, pp 201-207, February 2003.
[Lef-2005]
Anthony.Lefevre
Laurence.Migeville
Java.Fouladgar
Guy.Olivier
3D
computation of transformers overheating under nonlinear loads IEEE Transaction

On Magnetics, vol.41, No.5, pp 1564-1567, May 2005.
[Li-2005] Xiaosong.Li Qiaofu.Chen Jianbo.Sun Yu.Zhang Guzong.Long Analysis
of magnetic field and circulating current for HTS transformer windings IEEE Transactions
On Applied Superconductivity, Vol-15, No 3, pp 3808-3813, September 2005.
[Mao-2006] B.Maouche, M.Feliachi, N.Khenfer A half-analytical formulation for
the impedance variation in axisymmetrical modeling of eddy current non destructive testing
Eur. Phys. J. Appl. Phys. 33, 5967 (2006).
[Ras-2008] V. Rashtchi H. Shayeghi M. Mahdavi A. Kimiyaghalam E. Rahimpour Using
an improved PSO algorithm for parameter identification of transformer detailed model
International Journal of Electrical Power and Energy Systems Engineering - Summer 2008.
[Fre-2009] W. Frelin L. Berthet M. Petit J. C. Vannier Transformer Winding Losses

Evaluation when Supplying Non Linear Load Author manuscript, published in "The 44th
International Universities" Power Engineering Conference - UPEC 2009, Glasgow :United
Kingdom (2009).
[Jan-2009] arko Janic Zvonimir Valkovic eljko tih Helical windings magnetic field
in power transformers Electr Eng (2009) 91:161166 published online: 9 October 2009
Springer-Verlag.
[Hoc-2009] Labar Hocine, Rekik Badri, Bounaya Kamel, Kelaiaia Mounia Samira
Transformer Diagnosis Based on Coupled Circuits Method Modeling World Academy of
Science, Engineering and Technology 54 2009.
[Lin-2009] A. Linsuo Zeng B. Liangxian Zhang C. Dexin Xie D. Shan Bai E. Yan Wu
F. Yuanhai Xia Separation Calculation of High-Voltage Dry type Smoothing Reactor under
Harmonic Magnetic Fields IEEE Transaction On Magnetics, International Conference On
Electrical Machines and Systems (ICEMS), pp 1-4, 15-18 Nov, 2009.
[Mao-2009] Bachir.Maouche Rezak.Alkama Mouloud.Feliachi Semi-analytical calculation
of the impedance of a differential sensor for eddy current non-destructive testing NDT&E
International 42 (2009) 573 580-Science Direct.
[Ven-2009] Vicente.Venegas Jose.L.Guardado Serguei.G.Maximov Enrique.Melgoza
A computer model for surge distribution studies in transformer windings IEEE, 978-14244-3861-7/09, pp 451-457, 2009.
[Esl-2010] M. Eslamian B. Vahidi S.H. Hosseinian Combined analytical and FEM
methods for parameters calculation of detailed model for dry-type transformer Simulation
Modelling Practice and Theory 18 (2010) 390403-ScienceDirect.
Rfrences de transformateurs
[E.E.I-2008] Procs verbale de contrle de transformateur, Entreprise Electro-industries (exENEL Azzazga, Tizi-Ouzou).
[Bah-1999] Bahmane.Omar Chaouche.Reda Chami.Abdelhakim Processus de ralisation et
calcul dun transformateur mmoire de fin dtude DEUA en lectrotechnique, Universit
Mouloud Mammeri de Tizi-Ouzou.
Contribution la modlisation analytico-numrique

des transformateurs de puissance
Rsum: Le travail prsent dans ce mmoire concerne ltude et la mise en uvre dun
modle prdictif pour retranscrire le comportement lectromagntique du transformateur
de puissance. Un modle bas sur la mthode des circuits lectriques magntiquement
coupls permet la modlisation des enroulements primaire et secondaire, avec prise en
compte de leffet de peau et de proximit. Par contre la modlisation du noyau magntique
sappuis sur des hypothses simplificatrices, en considrant son matriau homogne, isotrope
et linaire qui est reprsent par des densits fictives de courant surfacique. Pour remdier
aux insuffisances de cette mthode semi-analytique et la difficult de modlisation du circuit
ferromagntique, nous avons adopt une mthode mixte dit un couplage hybride (CEMC-EF),
la mthode des circuits lectriques magntiquement coupls est applique a travers la
discrtisation des enroulements et le noyau magntique sera pris en charge par la mthode
des lments finis. Ce couplage permet de prendre en compte les non linarits et la variation
des grandeurs physiques des matriaux ferromagntiques cela nous offre lavantage dtudier
notre transformateur aussi en rgime statique quen rgime transitoire.
Mots cls : Transformateur de puissance, Circuits Electriques Magntiquement Coupls,
Couplage CEMC-MEF, Elments Finis.

Memoire Magister 4

Transféré par

Droits d'auteur :

Formats disponibles

Memoire Magister 4

Transféré par

Informations du document

Description originale:

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

Memoire Magister 4

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE LENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

FACULTE DE GENIE ELECTRIQUE ET DE LINFORMATIQUE

Mr. Madjid MEZIANI

Mr. Mhemed RACHEK

Maitre de Confrences A UMMTO

Mr. Nacer Eddine BENAMROUCHE

Mr. Salah HADDAD

Mr. Tahar OTMANE CHERIF

Jadresse des remerciements tout particuliers mon promoteur Monsieur Mhemed

Chapitre I : Etat de lart sur la modlisation des transformateurs

I.2 Principaux lments des transformateurs...

I.3 Principe de fonctionnement

I.4 Diffrents types de transformateurs..

I.4.1 Transformateur colonnes...

I.4.1.1 Les enroulements.....

I.4.1.2 Le circuit magntique..

I.4.1.3 Limitations des transformateurs colonnes....

I.4.2 Transformateur cuirass..............

I.4.2.1 Les enroulements

I.4.2.2 Les crans lectrostatiques.

I.4.2.3 Avantages et inconvnients du transformateur cuirass.....

I.4.3 Transformateur ordinaire de phase...

I.4.4 Transformateur immerg dans l'huile et transformateur sec...

I.4.4.1 Transformateur lhuile.....

I.4.4.2 Transformateur sec.....

I.5 Transformateurs triphass..

I.6 Groupe de couplage................ ...

I.6.1 Couplages normaliss...

I.6.1.1 Couplage toile toile..................

I.6.1.2 Couplage triangle toile

I.6.1.3 Couplage toile zigzag.

I.6-2 Choix du couplage dans les transformateurs usuels.

I.7 Contraintes subies par le transformateur

I.8 Rgimes de fonctionnement...

I.8.1 Rgime de surcharge

I.8.2 Rgime dsquilibr.....

I.8.3 Rgime transitoire

I.9 Protection des transformateurs contre les surtensions

I.10 Diffrentes approches de modlisation des transformateurs

I.11 Schmas quivalent en tant que modle...

I.11.1 Schma classique

I.11.2 Phnomnes lectromagntiques.

I.11.2.1 Couplage parfait (transformateur idal)

I.11.2.2 Couplage imparfait.................................................

I.11.2.3 Pertes cuivre.

I.11.2.4 Calcul des rsistances ..

I.11.2.5 grandeurs magntiques.

I.11.2.6 Energie magntique emmagasine...

I.11.2.7 Pertes fer...

I.11.3 Phnomnes lectrostatiques

I.11.3.1 Pertes dilectriques...

I.11.3.2 Calcul des capacits..

I.11.3.3 Capacit entre enroulements

I.11.3.4 Capacit entre conducteurs

I.11.3.5 Capacit entre phases...

I.11.3.6 Capacit entre enroulement et cuve..

I.12 Modlisation par les mthodes semi-analytiques.

I.13.1 Calcul des selfs inductances..