Université de Souk-Ahras
Faculté de sciences et de la technologie
Département de Génie électrique : Master 2 : Machines électriques
Module : TP_CAO des machines électriques
Année universitaire : 2023/2024
Enseignant : H ZAIMEN
TP5 : CAO d’une machine synchrone à aimants permanents MSAP
1. Modélisation électromagnétique de la MSAP sous FEMM :
Avant d’aller faire une étude de conception et optimisation da la MSAP, nous allons nous intéresser d’abord
par le comportement électromagnétique la machine : calcul de champ, calcul de couple, ainsi qu’une étude
paramétrique.
→La détermination du couple Cem et de l’induction dans l’entrefer Be peut s’effectuer à partir de différentes
méthodes. Les méthodes numériques comme la MEF, MDF ou MVF permettent de déterminer les grandeurs
magnétiques avec une bonne précision pour des géométries complexes. Pour cet objet, en utilisant un logiciel
de calcul du champ par MEF :
OctaveFEMM.
1.1 Modèles géométrique et
physique de la machine :
Il s’agit d’une machine synchrone à
aimants permanents déposés en surface du
rotor et à aimantation radiale.
La modélisation électromagnétique
s’effectue en 2 dimensions dans la vue en
coupe représentée sur la Fig.1 :
Le Tab.1 suivant récapitule tous les paramètres
géométriques et physiques de notre MSAP :
Fig.1 : Présentation de la machine synchrone à
aimants (un quart de la machine=1 pôle).
Tab.1 : Propriétés physiques et géométriques.
TP5 : CAO d’une machine synchrone à aimants permanents MSAP 1/7
Université de Souk-Ahras
Faculté de sciences et de la technologie
Département de Génie électrique : Master 2 : Machines électriques
Module : TP_CAO des machines électriques
Année universitaire : 2023/2024
Enseignant : H ZAIMEN
1.2 Modèles de la machine sous FEMM :
Les figures 2 et 3 montrent respectivement le modèle en 2D de notre machine sous le logiciel FEMM, et le
maillage éléments finis triangulaires EF.
1.3 Bobinage du stator :
L’enroulement dans toute machine est la pièce responsable de la création de champ
électromagnétique EM. Tout enroulement est caractérisé par les paramètres suivants :
Notre bobinage triphasé du stator est réalisé à pas diamétral (ouverture des bobines y=π électrique). Le
schéma d’enroulement du stator est illustré par la figure 4 :
TP5 : CAO d’une machine synchrone à aimants permanents MSAP 2/7
Université de Souk-Ahras
Faculté de sciences et de la technologie
Département de Génie électrique : Master 2 : Machines électriques
Module : TP_CAO des machines électriques
Année universitaire : 2023/2024
Enseignant : H ZAIMEN
1.4 Bloc matériaux :
▪ Culasse :
Pour la culasse de stator et rotor on utilise le
"M-19 Steel" qui se caractérise par la
courbe de magnétisation montrée en Fig 5 :
▪ Aimants permanents :
L’aimant de type "Néodyme-Fer-Bore" qui
se caractérise par les grandeurs :
►Champ coercitif : Hc= 106 A/m ;
►Perméabilité : µr= 1,07.
Fig. 5 courbe d’aimantation de ferromagnétique
▪ Cuivre : utilisé pour la culasse.
Pour les bobines du stator, on choisit le conducteur de type "10 AWG" dont la conductivité σ =58
MS/m
1.4.1 Commandes pour bibliothèque des matériaux :
Ferromagnétique :
mi_getmaterial('M-19 Steel'); % culasse stator et rotor
Aimants:
mi_addmaterial('NdFeB', mu_r , mu_r , Hc , 0 , 0.667 , 0,0,0,0,0,0); % PM
Air:
mi_getmaterial('Air'); % air
Bobines:
mi_getmaterial('10 AWG'); % cuivre
Fig. 6: Bibliothèque des matériaux.
TP5 : CAO d’une machine synchrone à aimants permanents MSAP 3/7
Université de Souk-Ahras
Faculté de sciences et de la technologie
Département de Génie électrique : Master 2 : Machines électriques
Module : TP_CAO des machines électriques
Année universitaire : 2023/2024
Enseignant : H ZAIMEN
1.4.2 Insérer les matériaux aux régions :
mi_addblocklabel(x,y);
mi_selectlabel(x,y);
mi_setblockprop(' blockname ', automesh, meshsize, ' incircuit ',
magdir, group,turns);
Exemples :
1/ insertion d’un matériau ferromagnétique :
mi_addblocklabel(x,y);
mi_selectlabel(x,y);
mi_setblockprop('M-19 Steel', 0, 0, 'none', 0, 1,0);
2/ insertion d’un aimant:
mi_addblocklabel(x,y);
mi_selectlabel(x,y);
mi_setblockprop('NdFeB', 0, 0, 'none', 30, 1,0);
1.5 Bloc circuit électriques :
Pour créer un circuit électrique on utilise de la commande suivante :
mi_addcircprop(' circuitname', i, circuittype)
circuittype : 0 for a parallel-connected circuit and 1 for a series-connected circuit,
Exemples: création des 3 phases du stator A, B et C (aller + retour) :
mi_addcircprop('A+',0,1);
mi_addcircprop('B+',0,1);
mi_addcircprop('C+',0,1);
mi_addcircprop('A-',0,1);
mi_addcircprop('B-',0,1);
mi_addcircprop('C-',0,1);
1.5.1 Insertion des circuits électriques des phases A, B et C :
mi_selectlabel(x,y);
mi_setblockprop(' blockname ', automesh, meshsize, ' incircuit ',
magdir, group,turns);
Exemples: insertion de la phase A aller :
mi_selectlabel(x,y);
mi_setblockprop('10 AWG', 0, 0, 'A+', 0, 2,Ns);
TP5 : CAO d’une machine synchrone à aimants permanents MSAP 4/7
Université de Souk-Ahras
Faculté de sciences et de la technologie
Département de Génie électrique : Master 2 : Machines électriques
Module : TP_CAO des machines électriques
Année universitaire : 2023/2024
Enseignant : H ZAIMEN
1.6 Conditions aux limites :
Pour insérer une CAL on utilise la commande suivante :
m_addboundprop (‘propname’, A0, A1, A2, Phi, Mu, Sig, c0, c1, BdryFormat): ajouter une nouvelle
condition aux limites avec le nom ‘propname’.
Exemple : insertion CAL de type Dirichlet: A=0 (sur la surface externe du stator) :
mi_addboundprop('A', 0, 0, 0, 0, 0, 0, 0, 0, 0)
mi_selectarcsegment(Rs_ext*cos(10*pi/180),Rs_ext*sin(10*pi/180))
mi_selectarcsegment(Rs_ext*cos(-10*pi/180),Rs_ext*sin(-10*pi/180))
mi_setarcsegmentprop(9,'A',0,2) ;
1.7 Enregistrer:
mi_saveas('[Link]') % Enregister le programme FEMM
1.8 Exécution et résultats :
mi_analyze : exécute le solveur magnétique ;
mi_loadsolution : charger et afficher la solution correspondant à la géométrie actuelle ;
1.9 Calcul des grandeurs I, V et Φ:
La commande mo_getcircuitproperties(’circuit’) permet de donner pour un circuit les
grandeurs suivantes :
- Current : courant transporter par le circuit ;
- Volts: chute de tension à travers le circuit ;
- Flux: flux balayé par le circuit : flux de phase par exemple.
1.10 Calcul d’autres grandeurs EM :
La commande mo_blockintegral(type)permet de calculer pour un bloc sélectionné les grandeurs
suivantes :
TP5 : CAO d’une machine synchrone à aimants permanents MSAP 5/7
Université de Souk-Ahras
Faculté de sciences et de la technologie
Département de Génie électrique : Master 2 : Machines électriques
Module : TP_CAO des machines électriques
Année universitaire : 2023/2024
Enseignant : H ZAIMEN
2. Travail demandé :
2.1 Calcul des lignes du champ (lignes équipotentielles) :
→ Tracer les lignes du champ pour les types d’aimants suivants :
NdFeB: μr=1,07 et Hc=106 A/m,
Alnico 6 : μr=3,3 et Hc= 59907 A/m,
Sm-Co: μr= 1,103 et Hc= 772000 A/m.
Commande utilisée : mi_addmaterial.
2.2 Induction magnétique :
→ Tracer la distribution de l’induction magnétique pour les types d’aimants suivants :
NdFeB: μr=1,07 et Hc=106 A/m,
Alnico 6 : μr=3,3 et Hc= 59907 A/m,
Sm-Co: μr= 1,103 et Hc= 772000 A/m.
→ Commenter les résultats trouvés en indiquant la valeur max de l’induction magnétique.
2.3 Etude paramétrique :
2.3.1 Influence de l’ouverture des aimants (θ_aim):
→Discuter l’influence de l’ouverture des aimants (θ_aim) sur le comportement magnétique de notre
machine (p=2 et aimant NdFeB), on choisit 04 états d’ouvertures d’aimants :
θ_aim= 30o, θ_aim= 45o, θ_aim= 60o, θ_aim= 75o.
→Tracer la répartition de l’induction magnétique pour les 04 cas d’ouverture d’aimants ;
→Tracer la courbe Bmax en fonction de θ_aim, commenter les résultats.
2.4 Flux de phase :
Le rotor fait un demi-tour : θ_méc=1800, et on choisit un pas de rotation pas=10o :
→Pour p=2 et aimant NdFeB, Tracer les formes d’ondes de flux par phase (a,b et c) : Φa, Φb et Φc.
Commandes utilisées : mi_addmaterial.
mo_getcircuitproperties
% Rotating
mi_selectgroup(1); % groupe (1) c'est le rotor
mi_moverotate(0, 0, theta(ii));
% Calcul le flux de phase --exemple de la phase A aller + retour
fA1= mo_getcircuitproperties('A+') ;
fA2= mo_getcircuitproperties('A-') ;
FluxA1(ii) = fA1(1,3);
FluxA2(ii) = fA2(1,3);
commenter les résultats.
TP5 : CAO d’une machine synchrone à aimants permanents MSAP 6/7
Université de Souk-Ahras
Faculté de sciences et de la technologie
Département de Génie électrique : Master 2 : Machines électriques
Module : TP_CAO des machines électriques
Année universitaire : 2023/2024
Enseignant : H ZAIMEN
2.5 Induction magnétique dans l’entrefer :
Tracer la courbe de l’induction magnétique B dans l’entrefer (B en fonction de : θ_méc) pour :
NdFeB: μr=1,07 et Hc=106 A/m,
Alnico 6 : μr=3,3 et Hc= 59907 A/m,
Sm-Co: μr= 1,103 et Hc= 772000 A/m.
Commandes utilisées : mo_getb.
commenter les résultats.
TP5 : CAO d’une machine synchrone à aimants permanents MSAP 7/7