Couplage magnéto-mécanique pour le dimensionnement de machines électriques Thème

Couplage magnéto-mécanique pour le dimensionnement de

machines électriques
1 2 1 1
ABLA DAHIA , FREDERIC DUBAS , M’HAMED BELHADI , GUILLAUME KREBS ,
1 1
LAURENT DANIEL ET CLAUDE MARCHAND
1
: GeePs, Laboratoire Génie électrique et électronique de Paris
UMR CNRS 8507, CentraleSupélec, Univ. Paris-Sud, Sorbonne Universités, UPMC Univ. Paris 06
3
, 11 rue Joliot Curie, Plateau de Moulon, 91192 Gif-sur-Yvette CEDEX, France.
[email protected]
2
: Département ENERGIE, Institut FEMTO-ST, UMR CNRS 6174, UBFC, Belfort, France [email protected]

Résumé : Ce papier présente un couplage magnéto-mécanique pouvant être utilisé dans une
procédure de dimensionnement d'un convertisseur électromécanique. Le modèle basé sur un circuit
magnétique équivalent (CME) non-linéaire et un modèle mécanique par éléments finis (EF) en deux
dimensions (2D) sont détaillés ainsi que la procédure globale. Le dispositif étudié s’apparente à un
actionneur à réluctance variable simplifié.

Mots-clés — Couplage magnéto-mécanique ; circuit magnétique équivalent ; modèle éléments finis ; lois de
comportement.

Les principales caractéristiques géométriques de

1. INTRODUCTION
l’actionneur sont décrites dans le Tableau 1. Les tôles
Au sein des machines électriques, plusieurs sources ferromagnétiques sont constituées d’un alliage fer-
de contraintes mécaniques peuvent affecter le silicium. Le bobinage statorique est constitué de deux
comportement magnétiques des matériaux (e.g., la bobines en série de Nts = 50 spires chacune parcourues
variation de perméabilité magnétique, les pertes fer par un courant Is = 10 A.
[1],…). On peut citer les efforts dus à l’inertie qui
peuvent fragiliser par exemple les ponts magnétiques
dans les rotors des machines synchrones [3]. Des
efforts résiduels peuvent apparaître suite à certains
procédés de fabrication (e.g., l’emboutissage,
l’usinage,…). La température et la découpe mises alors
en jeu modifient les propriétés magnétiques jusqu’à
plusieurs millimètres du bord de coupe [3]. Enfin, les
efforts magnétiques participent eux aussi aux
modifications de l’état de contrainte et donc influent
sur les caractéristiques magnétiques [4].
Du point de vue magnétique, les CMEs non-
linéaires offrent un compromis (temps de
calcul/précision des résultats) appréciable et sont
régulièrement utilisés en optimisation [5]. Nous
proposons donc dans ce papier le couplage d’un tel Fig. 1 Système étudié et son CME non-linéaire en position
de conjonction.
modèle avec un modèle mécanique. L’idée est de tenir
compte du comportement magnéto-mécanique lors
Tableau 1. Principales dimensions de l’actionneur à
d’un processus de dimensionnement. Dans cet objectif réluctance variable.
une structure magnétique simplifiée a été retenue. Sa
représentation magnétique à l’aide d’un CME non- Rayon Rayon Epaisseur Entrefer Ouverture
linéaire est d’abord présentée puis le modèle externe d’alésage de culasse polaire
mécanique associé est décrit. Le schéma de couplage statorique
est donné ainsi que des résultats qui en sont issus. 25,5 mm 10 mm 0,5 mm 46 deg.
50 mm
2. DISPOSITIF ETUDIE
La structure à modéliser magnétiquement et
mécaniquement est présentée sur la Fig. 1. Elle peut
être considérée comme une structure élémentaire
d'actionneur à réluctance variable à double saillance.

La Revue 3EI n°88

Avril 2017 17
 Thème Couplage magnéto-mécanique pour le dimensionnement de machines électriques

= ∙
MMF = ⋮ =−
3. MODELE MAGNETIQUE
avec (5d)
À partir des équations électromagnétiques, le CME
(dont les flux de mailles ψ sont les inconnues), peut avec
être exprimé par [6] :
= ⋱ et = ⋱ (5e)
F-χ· · χ T· =0 (1)

⋯
χ = et χ = −χ de dimension ×
− ⋯ −
où [R] représente la matrice diagonale des (5f)
réluctances, ["χ" ] la matrice topologique (ou
d'incidence) et ["F" ] le vecteur des MMFs de mailles
dont N = 5 représente le nombre de discrétisation à la
qui est défini par :
moitié du rotor.
F=χ· (2)
4. MODELE MECANIQUE
avec [MMF] le vecteur des MMFs de branches qui
Le domaine Ω étudié (supposé élastique et
dépendent des sources électromagnétiques (e.g., les
aimants permanents ou le bobinage parcouru par des homogène) de frontière dΩ, soumis à l’effet de forces
courant de forme d'onde quelconque). Connaissant [ψ], volumiques Fi et surfaciques fi, est soumis aux
le vecteur des inductions magnétiques [B] peut être conditions d’équilibre définies par :
déterminé par : σij,j +Fi =0 (6)

= ∙ σij +fi =0 (7)

= χ T·
(3)
Les conditions aux limites sont définies sur les
où [ϕ] est le vecteur des flux de branches, [S] la déplacements. La relation liant les déformations ε aux
matrice diagonale des surfaces de passage des flux de déplacements u est définie par :
branches. Le vecteur des champs magnétiques est 1 ∂ui ∂uj
défini selon : εij = + (8)
2 ∂xj ∂xi

= ∙ ∙
Les contraintes sont liées aux déformations par la
(9)
loi d’élasticité linéaire (loi de Hooke):
σij =Cijkl εkl (9)
avec [µr] la matrice diagonale des perméabilités
magnétiques relatives. Le problème d’élasticité est exprimé selon :
Analytiquement, (1) peut être résolue linéairement
de façon itérative par la méthode du point fixe (µr est σij εij dx = F.v dx + f.v ds (10)
constant à chaque itération) [6]. La caractéristique
B(H) du circuit ferromagnétique (courbe
anhystérétique) est décrite par une fonction Une discrétisation par EF est ensuite employée.
mathématique, dite de Marocco [7]. S’appuyant sur un maillage triangulaire, et avec
Dans l'exemple d'étude, le CME est représenté sur comme termes sources les forces discrétisées aux
la Fig. 1 pour la position de conjonction. Pour la nœuds, les inconnues liées aux déplacements sont
position choisie, le modèle est simple et rapide à déterminées. Les contraintes sont ensuite calculées.
résoudre. Les matrices [R], ["χ" ] et [µr] ainsi que le À titre d’illustration (le nombre d’éléments est
vecteur des sources électromagnétiques [MMF], définis choisi délibérément faible) et pour une vitesse de
dans (1) ~ (4), sont représentés par : rotation de 10 000 tr/min, la répartition des efforts
nodaux dus à la rotation (efforts centrifuges) et les
déplacements obtenus sont donnés sur la Fig. 2. La
⋱
répartition des contraintes radiales et tangentielles sont
données sur la Fig. 3. Elles montrent que le rotor est
= (5a) sollicité en bi-traction (valeurs positives des

⋱
contraintes). L’amplitude maximale est atteinte au
centre car les extrémités du rotor vont "avoir tendance"
à s’éloigner, voir Fig. 4.
− − −
χ = (5b)
− − −
χ χ

= =
= ⋱ avec (5c)

La Revue 3EI n°88

18 Avril 2017
 Couplage magnéto-mécanique pour le dimensionnement de machines électriques Thème

5. COMPORTEMENT MAGNETO-
MECANIQUE
Le comportement magnéto-mécanique
anhystérétique est ici pris en compte à l’aide d’une
version simplifiée d’un modèle de comportement
magnéto-élastique multi-échelle[8]. La multi-axialité
des contraintes agissant sur les caractéristiques
magnétiques est bien considérée. Dans le cas d’un
matériau initialement isotrope, le modèle repose sur
trois paramètres :
[N] [m] • l’aimantation à saturation ;
• la constante de déformation de
Fig. 2 Efforts nodaux et déplacements dus à la rotation sur magnétostriction à saturation ;
une partie du rotor.
• et d'un paramètre d’ajustement.
Pour trois niveaux de contraintes uniaxiales (0 MPa
et jusqu'à ±55 MPa) dans une pièce rectangulaire en
Fe-Si, les caractéristiques B(H) ont été calculées au
centre de cette dernière [voir Fig. 5]. On peut noter
qu’en traction les propriétés magnétiques de l’acier
sont améliorées. Pour une même valeur de champ
magnétique et en dehors de la saturation, il y a une
augmentation significative de la perméabilité
magnétique. Néanmoins, dans le cas d’un régime
magnétique variable, les pertes fer risquent d’être plus
importantes.
[Pa] [Pa]
Fig. 3 Contraintes radiales et tangentielles au sein du 6. COUPLAGE MAGNETO-MECANIQUE
rotor. Le schéma décrivant le couplage magnéto-
mécanique est donné sur la Fig. 6. Le premier calcul
magnétique est effectué sans tenir compte des
contraintes mécaniques puis le problème mécanique est
résolu. À partir de la direction du champ et des
contraintes, le modèle de comportement magnéto-
mécanique donne une nouvelle caractéristique
Bi =f(Hi) pour chaque réluctance R4-i et R5-i du rotor
(vecteurs de deux dimensions). Les contraintes sont
calculées au barycentre des zones délimitées par les
réluctances. Si on ne tient compte que des efforts dus à
l’inertie, le calcul mécanique est lancé une seule fois.
Fig. 4 Maillages du rotor, initial et après application des Le modèle magnétique associé au modèle de
déplacements (facteur × 5.104).
comportement est alors exécuté jusqu’à convergence
des résultats.

Contraintes σ
Modèle (MPa)
mécanique
Comportement
Modèle
magnéto-
magnétique
mécanique
Txx Txx Bi=f(Hi, σ=0) {Bi, Hi} Bi=f(Hi, σ)

Fig. 6 Schéma du couplage magnéto-mécanique

Le critère d’arrêt de la boucle magnéto-mécanique

est de considérer un écart minimum entre deux
Fig. 5 Caractéristiques B(H) pour différents niveaux de itérations des niveaux d’induction magnétique locale.
contraintes uniaxiales. Dans le même esprit, il est aussi envisageable de
comparer les caractéristiques Bi=f(Hi) locales.

La Revue 3EI n°88

Avril 2017 19
 Thème Couplage magnéto-mécanique pour le dimensionnement de machines électriques

7. RESULTATS ET COMMENTAIRES modèle mécanique analytique [9]. Pour le cas traité, les
résultats obtenus ont montré qu’il n’y avait peu ou pas
Pour une vitesse de 1000 tr/min, le couplage
de modification globale en termes de flux capté par
magnéto-mécanique est résolu. L’induction magnétique
exemple. Localement et en fonction du niveau de
dans les reluctances R4-1 à R4-5 (celles présentes dans la
contraintes, une modification des caractéristiques
moitié du rotor) vaut 0,78 T et les caractéristiques
magnétiques a été constatée.
Bi=f(Hi) sont présentées sur la Fig. 7. On peut
constater, pour cette vitesse, que le niveau de
contraintes n’est pas assez important pour modifier les
caractéristiques magnétiques du rotor.
Pour une vitesse de 30 000 tr/min [voir Fig. 8], les
contraintes (ici en traction), modifient les
caractéristiques Bi=f(Hi) du rotor de manière favorable,
du point de vue magnétique. La perméabilité est
améliorée à mesure que l’on se rapproche du centre du
rotor à cause du niveau de contrainte (traction) qui
augmente lui aussi. Cependant, le niveau d’induction
magnétique reste pratiquement le même qu’à 1 000
tr/min, d’une part à cause des ampères-tours
consommés dans l’entrefer nettement supérieurs à ceux
consommés dans le fer et d’autre part à la discrétisation
verticale du rotor.
Fig. 7 Caractéristiques Bi=f(Hi) de réluctances R4-i au
À 60 000 tr/min [voir Fig. 9], les caractéristiques rotor pour une vitesse de 1 000 tr/min.
Bi=f(Hi) sont nettement modifiées mais de nouveau il
n’y pas de modification notable de l’induction
magnétique au rotor.
La Fig. 10 montre les caractéristiques Bi=f(Hi) pour
la réluctance la plus proche du centre du rotor
(i.e., R4-5) et pour deux vitesses de rotation. Pour le
point de fonctionnement (0,78 T), la perméabilité est
pratiquement multipliée par deux entre les deux
vitesses.
Concernant les temps de calcul, comme la
géométrie du rotor est simple la résolution du problème
mécanique est courte (entre 15 et 30 sec selon la
densité du maillage mais le programme sous Matlab
peut être optimisé). Si on ne tient compte que des
efforts dus à l’inertie (hors pression magnétique et
contraintes résiduelles dues à l’emboutissage), le
Fig. 8 Caractéristiques Bi=f(Hi) de réluctances R4-i au
problème mécanique peut être lancé une seule fois. Il rotor pour une vitesse de 30 000 tr/min.
en est presque le même pour le problème magnétique
dans ce cas présent, à cause de la discrétisation
magnétique verticale du rotor. Le sens du champ est de
ce fait imposé, et il est possible pour réduire les temps
de calculs de n’exécuter que deux fois le problème
magnétique.

8. CONCLUSIONS
Dans ce papier, le couplage entre un CME non-
linéaire et un modèle mécanique est présenté. L’intérêt
est d’utiliser un tel couplage dans une procédure de
dimensionnement de structure de machine faisant
apparaître des zones fortement contraintes
mécaniquement (pont magnétique, application haute
vitesse…). Une approche plus fine et qui permettrait de
tenir compte de plusieurs positions mécaniques est à
envisager. Il serait alors nécessaire de discrétiser le Fig. 9 Caractéristiques Bi=f(Hi) de réluctances R4-i au
rotor avec des réluctances selon les différentes rotor pour une vitesse de 60 000 tr/min.
directions afin de ternir compte des répartitions des
contraintes et du champ magnétique. Les temps de
calculs pourraient être réduits grâce à l’utilisation d’un

La Revue 3EI n°88

20 Avril 2017
 Couplage magnéto-mécanique pour le dimensionnement de machines électriques Thème

[3] Y. Kai, Y. Tsuchida, T. Todaka, and M.

Enokizono, « Effect of local residual stress in
rotating machine core on vector magnetic
property », in Proc. ICEM, Rome, Italy, 06-08
Sept., 2010.
[4] Z. Han, H. Yang, and Y. Chen, « Investigation
of the rotor mechanical stresses of various
interior permanent magnet motors », in Proc.
ICEMS, Tokyo, Japan, 15-18 Nov., 2009.
[5] M. Hage Hassan, G. Remy, G. Krebs, and C.
Marchand, « Radial output space mapping for
electromechanical systems design », COMPEL,
Fig. 9 Caractéristiques Bi=f(Hi) de la réluctance proche vol. 33, no. 3, pp. 965-975, 2014.
du centre du rotor et pour deux vitesses.
[6] R. Benlamine, F. Dubas, S-A. Randi, D.
Lhotellier, and C. Espanet, « Modeling of an
9. REMERCIEMENTS Axial-Flux Interior PMs Machine for an
Automotive Application Using Magnetic
Ce travail s’inscrit dans le cadre du projet COCTEL Equivalent Circuits », in Proc. ICEMS, Pattaya,
(Conception Optimale des Chaines de Traction Thailand, 25-28 Oct., 2015.
Électrique) piloté par RENAULT-SAS et financé par
l’ADEME (Agence de l’Environnement et de la [7] A. Marrocco, « Analyse numérique de
Maîtrise de l’Énergie). problèmes d'électrotechnique », Ann. Sc. Math.,
vol. 01, pp. 271-296, 1977.
10. REFERENCES [8] L. Bernard, and L. Daniel, « Effect of stress on
magnetic hysteresis losses in a switched
[1] R. Akaki, Y. Takahashi, K. Fujiwara, M. reluctance motor: application to stator and rotor
Matsushita, N. Takahashi, and M. Morita, « shrink-fitting », IEEE Trans. on Magnetics, vol.
Effect of magnetic property in bridge area of 51, no. 9, 7002513, Sept. 2015.
IPM motors on torque characteristics », IEEE
Trans. on Magnetics, vol. 49, no. 5, pp. 2335- [9] F. Chai, Y. Li, P. Liang, and Y. Pei, «
2338, May 2013. Calculation of the maximum mechanical stress
on the rotor of interior permanent magnet
[2] J. Kolehmainen, and J. Ikäheimo, « Motors with synchronous motors », IEEE Trans. on
buried magnets for medium-speed applications Industrial Electronics, in Press, 2016.
», IEEE Trans. on Energy Conversion, vol. 23,
no. 1, pp. 86-91, March 2008.

La Revue 3EI n°88

Avril 2017 21