Université Abdelhamid Ibn Badis MOSTAGANEM
Faculte Des Sciences et de la Technologie
TP Méthodes Numérique et Optimisation Departement du Genie Electrique
TP N°2 : Etude et simulation en 2D le champ magnétique dans un moteur électrique
bipolaire avec la Méthode des éléments finis
Objectif du TP :
Etude et simulation en 2D le champ magnétique dans un moteur électrique bipolaire c’est-à-dire
trouvez
rouvez le champ magnétique statique induit par les enroulements du stator dans un moteur
électrique bipolaire. Les aimants, les moteurs électriques et les transformateurs sont des domaines
dans lesquels des problèmes liés à la magnétostatique peuvent être rencontrés. La «statique»
implique que le taux de changement
ngement dans le temps est lent.
Pour la simulation on utilise la toolboxe pdetool du MATLAB.
Enoncé du TP
a géométrie :
1- Création de la
La géométriesese compose de trois régions: Deux pièces ferromagnétiques: le stator et le rotor (acier
du transformateur) l'entrefer
'entrefer entre le stator et le rotor ainsi que laa bobine de cuivre d'induit
transportant le courant continu (voir figure 1).
Figure 1
La géométrie de ce moteur électrique est [Link] modèle est une union de cinq cercles et de
deux [Link] réduction au premier quadrant est réalisée par intersection avec un carré.
Pour dessiner la géométrie, entrez les commandes suivantes dans la fenêtre de commande
MATLAB :
pdecirc(0,0,1,'C1')
pdecirc(0,0,0.8,'C2')
pdecirc(0,0,0.6,'C3')
pdecirc(0,0,0.5,'C4')
pdecirc(0,0,0.4,'C5')
pderect([-0.2
0.2 0.2 0.2 0.9],'R1')
0.9],
pderect([-0.1
0.1 0.1 0.2 0.9],'R2')
0.9],
pderect([0 1 0 1],'SQ1')
)
Réduisez le modèle au premier quadrant.
Pour ce faire, entrez= (C1 + C2 + C3 + C4 + C5 + R1 + R2) * SQ1 dans le champ définir la
formule.
TP Methode Numériques et OPT Chagé du TP :BENAOUALI
BENAOUALI Mohamed __ 2021
Master 1ELT_Ind et ER 1
Université Abdelhamid Ibn Badis MOSTAGANEM
Faculte Des Sciences et de la Technologie
TP Méthodes Numérique et Optimisation Departement du Genie Electrique
Supprimez les bordures
rdures de sous-domaine
sous inutiles jusqu'à ce que la géométrie se compose de
quatre sous-domaines:
domaines: le rotor (sous-domaine
(sous 1), le stator (sous-domaine
domaine 2), l'entrefer (sous
(sous-
domaine 3) et la bobine (sous--domaine 4). Voir figure 2
Figure 2
quation : définir l’équation suivante
2- Déclaration de l’équation
B=µH
B = ∇A -∇ . ( (1/ µ)∇A) = J
où B est la densité de flux magnétique, H est l'intensité du champ magnétique, J est la densité de
courant et µ est la perméabilité magnétique du matériau.
A : potentiel de vecteur magnétique.
coefficients la PDE sur chaque sous-domaine
3- Spécifiez les coefficientsde domaine et spécifiez les coefficients
suivants:
Coil (bobine) : µ = 4*pi*10^(-7) H/m, J = 10 A/m2.
- Dans le stator et le rotor, µ est défini par
µmax= 5000, µmin= 200, et c = 0.05 sont des valeurs qui pourraient représenter l'acier dutransformateur.
Stator et rotor: µ = 4*pi*10^(-7)*(5000./(1+0.05*(ux.^2+uy.^2))+200)
7)*(5000./(1+0.05*(ux.^2+uy.^2))+200) H/m,
Où ux.^2+uy.^2 egal à |∇A |2, J = 0 (courant nul).
l’entrefer (Air gap)) : µ = 4*pi*10^(-7) H/m, J = 0.
4- Spécifiez les conditions aux limites.
limites
spécifiez une condition aux limites de Neumann avec g = 0 and q = 0 suivant l’axe des X
X.
Toutes les autres limites ont une condition aux limites Dirichlet avec h = 1 et r = 00.
TP Methode Numériques et OPT Chagé du TP :BENAOUALI
BENAOUALI Mohamed __ 2021
Master 1ELT_Ind et ER 2
Université Abdelhamid Ibn Badis MOSTAGANEM
Faculte Des Sciences et de la Technologie
TP Méthodes Numérique et Optimisation Departement du Genie Electrique
5- Initialisez le maillage en sélectionnantMesh-InitializeMesh.
sélectionnant
6- Résoudre : Choisissez le solveur non linéaire. choisir d'utiliser le solveur adaptatif avec le
solveur nonlinéaire.
Utiliser : Options Axes Limits,
Limits ajustez les limites des axes si né[Link] exemple, utilisez la
case à cocher Auto.
Question :
TP N°2 : Etude et simulation en 2D le champ magnétique
1- Tracer la densité de flux magnétique B à l'aide des flèches et donner un commentaire ?
…………………………………………………………
…………………………………………………………….
……………………………………………………………
2- Tracer la densité de flux magnétique B à l'aide des flèches et des contours et donner un
commentaire ?
………………………………………………………………………………….
………………………………………………………………………..
3- Tracer le champ magnétique H à l'aide des flèches et des contours et donner un commentair
commentaire
…………………………………………………………………………….
………………………………………………………………
4- Tracer le vecteur potentiel magnétique A à l'aide des flèches et des contours et donner un
commentaire ?
………………………………………………………………….
………………………………………………………………..
5- Donner les 5 premiers points de solution ?
………………………………………………………..
6- Donner le nombre d’élément et la Qualité du maillage ?
TP Methode Numériques et OPT Chagé du TP :BENAOUALI
BENAOUALI Mohamed __ 2021
Master 1ELT_Ind et ER 3