Travaux Pratiques N 3: Calcul Du Champ Dans Un Circuit Magnétique

Université de Souk-Ahras
Faculté de sciences et de la technologie

Département de Génie électrique : Master 1- Machines Electriques
Module: TP Champ magnétique dans les machines électriques
Année universitaire : 2021/2022
Enseignant : H ZAIMEN
Travaux Pratiques N03 : Calcul du champ dans un circuit magnétique
1. Introduction :
Dans ce TP de simulation on va s’intéresser par le calcul de la répartition du champ
magnétique dans un circuit magnétique en utilisant les commandes et fonctions de
PDETOOL-Matlab.
2. Application en 2D : circuit magnétique sans entrefer :

La géométrie en 2D de circuit magnétique est illustrée dans la figure 1 :
5 mm
2 mm
5 mm
16 mm 30 mm
2 mm
2 mm
40 mm
Fig.1 dimensions du circuit magnétique.
3. Paramètres physiques :
Densité de courant Source Js= 5 106 A/m2
Perméabilité magnétique du Fer μFer = 2000. μ0 [H/m]
Perméabilité magnétique des bobines μbobines= 4 π 10-7 H/m
4. Dessin de la géométrie :
PDE-Toolbox permet de modéliser une très grande variété de systèmes électromagnétiques.
Une surface 2D peut être réalisée par l’utilisation des droites et des cercles.
TP Champ magnétique dans les machines électriques 1

Code Matlab Script de la géométrie :

% (1) dessin de la géometrie
% (1-1) Boite d'aire R1
x0=0.1;y0=0.1;
R1=[2 0 x0 0 0 1 0
2 x0 x0 0 y0 1 0
2 x0 0 y0 y0 1 0
2 0 0 y0 0 1 0];
% (1-2) Circuit magnétique R2

L=40e-3; H=30e-3; ep=5e-3; L2=30e-3; H2=20e-3;
x=(x0-L)/2; y=(y0-H)/2;
R2=[2 x x+L y y 2 1
2 x+L x+L y y+H 2 1
2 x+L x y+H y+H 2 1
2 x x y+H y 2 1
% -------------------------------------------
2 x+ep x+ep+L2 y+ep y+ep 1 2
2 x+ep+L2 x+ep+L2 y+ep y+ep+H2 1 2
2 x+ep+L2 x+ep y+ep+H2 y+ep+H2 1 2
2 x+ep x+ep y+ep+H2 y+ep 1 2];
% (1-3) Bobine ALLER R3

Lb=2e-3; Hb=16e-3; ep2=2e-3;
R3=[2 x+ep+ep2 x+ep+ep2+Lb y+ep+ep2 y+ep+ep2 3 1

2 x+ep+ep2+Lb x+ep+ep2+Lb y+ep+ep2 y+ep+ep2+Hb 3 1
2 x+ep+ep2+Lb x+ep+ep2 y+ep+ep2+Hb y+ep+ep2+Hb 3 1
2 x+ep+ep2 x+ep+ep2 y+ep+ep2+Hb y+ep+ep2 3 1];
% (1-4) Bobine RETOUR R4
R4=[2 x-ep2-Lb x-ep2 y+ep+ep2 y+ep+ep2 4 1

2 x-ep2 x-ep2 y+ep+ep2 y+ep+ep2+Hb 4 1
2 x-ep2 x-ep2-Lb y+ep+ep2+Hb y+ep+ep2+Hb 4 1
2 x-ep2-Lb x-ep2-Lb y+ep+ep2+Hb y+ep+ep2 4 1];
Geometry=[R1;R2;R3;R4]';
figure(1),pdegplot(Geometry);
5.Le Maillage :
Le maillage éléments finis EF d’un dispositif est décrit par trois matrices de nœuds, de
connectivités et de triangles (p,e,t) :
% (2) MAILLAGE
[p,e,t]=initmesh(Geometry);
[p,e,t]=refinemesh(Geometry,p,e,t,1);
figure(2),pdeplot(p,e,t);
% Nbre des nœuds et triangles
Nno=size(p,2);Ntr=size(t,2);
% Cordonnées des noeuds
xp=p(1,:)'; yp=p(2,:)';

6. Les coefficients :
On va affecter à chaque élément du maillage les coefficients correspondants (µ, σ et J) qui
diffère d’une région à une autre.
Problème magnétostatique en 2D : A=(0,0,Az) , Js=(0,0,Jsz) et B=(Bx, By, 0) ;
L’équation du champ magnétostatique en termes de potentiel vecteur A :
 1 
rot  rot A z  = J sz 1 
    −.  A z  = J sz  −. ( cu ) +au = f (elliptic PDE )
  
div A = 0
% (3)Les coefficients a,c,f:

c=zeros(1,Ntr); a=zeros(1,Ntr); f=zeros(1,Ntr);
mu0=4*pi*1e-7;
muFer=2000*mu0; muB=mu0; Js=5*1e6;
for ii=1:Ntr % Boucle sur les élements finis

if t(4,ii)==1 % Region 1: boite d'air
c(1,ii)=1/mu0;
a(1,ii)=0;
f(1,ii)=0;
elseif t(4,ii)==2 % Region 2: Fer
c(1,ii)=1/muFer;
a(1,ii)=0;
f(1,ii)=0;
elseif t(4,ii)==3 % Region 3: Bobine aller
c(1,ii)=1/muB;
a(1,ii)=0;
f(1,ii)=+Js;
elseif t(4,ii)==4 % Region 4: Bobine retour
c(1,ii)=1/muB;
a(1,ii)=0;
f(1,ii)=-Js;
end
end
7.Construction des matrices globales-Assemblage :

Les matrices élémentaires sont assemblées afin d’obtenir les matrices globales du système
algébrique : [M]. [A] = [F].
%(4) Assemblage M & F ---

[M1,M2,F]=assema(p,t,c,a,f);
M=M1+M2;

8. Les conditions aux limites CAL:

On va imposer une condition de type Dirichlet sur les 04 limites du domaine d’étude : Az=0.
%(5) CAL..magnetiques ------------
%..... Limites magnetiques.. 4 lignes..........
L1=find(p(2,:)==0 & p(1,:)<=x0);
L2=find(p(2,:)<=y0 & p(1,:)==x0);
L3=find(p(2,:)==y0 & p(1,:)<=x0);
L4=find(p(2,:)<=y0 & p(1,:)==0);
nox=[L1 L2 L3 L4]; % tous les noeuds des limites
%..... Conditions de type Dirichlet..Az=0........

%..... M(nox,nox)=1 & M(nox,:)=0 & Az(nox)=0 ....
for n=1:length(nox)
M(nox(n),:)=0; M(:,nox(n))=0; M(nox(n),nox(n))=1;
Az(nox(n))=0;
F=F-M(:,nox(n))*Az(nox(n));
end
9. La résolution du système algébrique [M].[Az] = [F]:

La résolution est représentée par un vecteur qui donne les valeurs du potentiel vecteur
magnétique Az dans tous les nœuds de maillage :
%(6)Résolution M*A=F ------
Az=inv(M)*F;
10.Tracer les lignes équipotentielles (lignes du champ) :
%(7) Equipotential lines ----

figure(3);
pdeplot(p,e,t,'xydata',Az,'xystyle','off','contour','on','levels',50,'color
bar','on')
11. Calcul de l’induction magnétique B :
%(8) Calcul de B ----

[dA_dx,dA_dy]=pdegrad(p,t,Az);
Bx=dA_dy; By=-dA_dx;
B=sqrt(Bx.^2+By.^2);

Travaux Pratiques N 3: Calcul Du Champ Dans Un Circuit Magnétique

Transféré par

Droits d'auteur :

Formats disponibles

Travaux Pratiques N 3: Calcul Du Champ Dans Un Circuit Magnétique

Transféré par

Informations du document

Titre original

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

Travaux Pratiques N 3: Calcul Du Champ Dans Un Circuit Magnétique

Transféré par

Droits d'auteur :

Formats disponibles

Université de Souk-Ahras

Faculté de sciences et de la technologie

Travaux Pratiques N03 : Calcul du champ dans un circuit magnétique

2. Application en 2D : circuit magnétique sans entrefer :

Densité de courant Source Js= 5 106 A/m2

Perméabilité magnétique du Fer μFer = 2000. μ0 [H/m]

Perméabilité magnétique des bobines μbobines= 4 π 10-7 H/m

TP Champ magnétique dans les machines électriques 1

Code Matlab Script de la géométrie :

% (1-2) Circuit magnétique R2

% (1-3) Bobine ALLER R3

R3=[2 x+ep+ep2 x+ep+ep2+Lb y+ep+ep2 y+ep+ep2 3 1

% (1-4) Bobine RETOUR R4

R4=[2 x-ep2-Lb x-ep2 y+ep+ep2 y+ep+ep2 4 1

TP Champ magnétique dans les machines électriques 2

% (3)Les coefficients a,c,f:

for ii=1:Ntr % Boucle sur les élements finis

7.Construction des matrices globales-Assemblage :

%(4) Assemblage M & F ---

TP Champ magnétique dans les machines électriques 3

8. Les conditions aux limites CAL:

%..... Conditions de type Dirichlet..Az=0........

9. La résolution du système algébrique [M].[Az] = [F]:

10.Tracer les lignes équipotentielles (lignes du champ) :

%(7) Equipotential lines ----

11. Calcul de l’induction magnétique B :

%(8) Calcul de B ----

TP Champ magnétique dans les machines électriques 4

Vous aimerez peut-être aussi