See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/280881618

Méthodes de calcul des sources pour des bobines non-maillées avec la

formulation éléments finis A − V

Conference Paper · June 2015

CITATIONS READS

0 144

7 authors, including:

Pauline Ferrouillat Christophe Guérin

Altair Engineering France Altair Engineering
19 PUBLICATIONS   10 CITATIONS    93 PUBLICATIONS   795 CITATIONS   

SEE PROFILE SEE PROFILE

Gerard Meunier Brahim Ramdane

French National Centre for Scientific Research G2Elab
491 PUBLICATIONS   5,216 CITATIONS    79 PUBLICATIONS   298 CITATIONS   

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

On the theory of bing bang View project

FLUX 10.4 View project

All content following this page was uploaded by Pauline Ferrouillat on 11 August 2015.

The user has requested enhancement of the downloaded file.

 Méthodes de calcul des sources pour des bobines non-maillées avec la
formulation éléments nis A − V
P. Ferrouillat1,2 , C. Guérin2 , G. Meunier1 , B. Ramdane1 , P. Dular3 , P. Labie1 et D. Dupuy2
1
Université Grenoble Alpes, CNRS, G2Elab, Grenoble, France
2
CEDRAT, Meylan, France
3
Université de Liège, Liège, Belgique

E-mail : pauline.ferrouillat@g2elab.grenoble-inp.fr

Résumé  Dans cet article, diérents calculs du du maillage, Nl et Wk les fonctions de forme nodale et
champ électromagnétique source sont décrits pour d'arête. Les termes de droite (3b)(3c) sont déduits de
la formulation éléments nis A − V avec des bobines
Z
Wk ·Js dΩ, avec Js la densité de courant source telle
non-maillées. Ces diérents calculs sont appliqués Ω
sur une machine asynchrone pour dénir la meilleure que Js = rot Hs . Le terme (3c) avec la condition n×Hs
approche. Cette source doit pouvoir ensuite être uti- est nul quand Γ est à l'inni. Mais quand il est intégré
lisée sur un domaine réduit xe ou en mouvement. sur un domaine réduit, il ne s'annule pas.

I. INTRODUCTION
III. CALCUL DES TERMES SOURCES

Les formulations en T−φ et en T−φr sont largement

utilisées. Cependant, ces formulations requièrent l'intro- Pour des bobines non-maillées, le champ magné-
duction de coupures articielles [1], contrairement à la tique source Hs et le potentiel vecteur magnétique
formulation en A − V [2]. En général, les bobines sont source As sont calculés analytiquement avec la loi de
maillées avec la formulation en A − V , mais dans cet Biot-Savart [4] tels qu'avec ν0 la réluctivité du vide :
article, la prise en compte des bobines non-maillées est
proposée. Diérents calculs de source pour des bobines Hs = ν0 rot As . (5)
non-maillées sont décrits pour la formulation en A − V
[3]. Ces calculs sont étudiés sur une machine asynchrone Or, As et Hs doivent être évalués à chaque fois qu'ils
pour dénir la meilleure source. Enn, cette source est sont utilisés. Pour éviter de les recalculer, trois calculs
calculée par domaine, xe ou en mouvement. de termes sources sont présentés dans cette partie.
Le champ magnétique source discrétisé Hj est inter-
polé à partir d'Hs avec la norme de L2 (Ω) telle que [3] :
II. FORMULATION EN A−V
Z Z
Les potentiels vecteur magnétique A et scalaire élec- Wk · Hj dΩ = Wk · Hs dΩ. (6)
trique V sont dénis en magnétodynamique tels que [2] : Ω Ω

B = rot A, (1) L'espace des éléments d'arêtes est H(curl, Ω), et

∂A non L2 (Ω). Cependant, l'interpolation d'Hj avec
E = − − grad V, (2) la norme ou semi-norme d'H(curl, Ω) nécessite des
∂t
sources impossibles à calculer analytiquement pour des
avec B l'induction magnétique, E le champ électrique.
bobines non-maillées. Dans la suite, le potentiel vecteur
Le potentiel A est discrétisé avec les éléments d'arêtes,
source discrétisé Aj est interpolé sur le bon espace.
alors que V est discrétisé avec les éléments nodaux.
Soit le domaine Ω de frontière Γ. Soit n la normale La pré-résolution d'Aj interpolée à partir d'As et
à Γ. Soit Ωb le domaine des bobines. En appliquant la d'Hs avec la norme d'H(curl, Ω) est telle que [5] :
méthode de Galerkin sur les équations de Maxwell avec Z Z
les potentiels A (1) et V (2), la formulation en A−V en rot Wk · rot Aj dΩ + Wk · Aj dΩ =
magnétodynamique sur le domaine Ω est telle que [3] : ZΩ ZΩ
1
Z  
∂A
 rot Wk · Hs dΩ + Wk · As dΩ. (7)
νrotWk · rotA + σWk · + gradV dΩ (3a) Ω ν0 Ω
Ω ∂t
Alors que la pré-résolution d'Aj interpolée à partir d'Hs
Z
= rot Wk · Hs dΩ (3b) avec la semi-norme d'H(curl, Ω) est telle que [3] :
Ω
Z
+ Wk · (n × Hs ) dΓ, ∀ k = 1 · · · na , (3c)
Z Z
1
Γ rot Wk · rot Aj dΩ = rot Wk · Hs dΩ. (8)
Z   Ω Ω ν0
∂A
σ grad Nl · + grad V dΩ = 0, ∀ l = 1 · · · nn ,
Ω ∂t Les équations (6)(7)(8) sont les pré-résolutions des
(4) termes sources. La formulation en A−V avec ces termes
avec Hs le champ magnétique source, ν la réluctivité, σ sources est résolue avec les équations (3)-(4) en rempla-
la conductivité, nn et na le nombre de n÷uds et d'arêtes çant Hs dans le terme (3b)(3c) par Hj ou ν0 rotAj (5).
 IV. CALCUL DE LA SOURCE SUR UN

DOMAINE RÉDUIT

Soit Ωred un domaine réduit tel que Ωb ⊂ Ωred ⊂ Ω,

avec sa frontière Γred telle que Ωb n'est pas en contact
avec Γred . En général, les machines ont une partie du
domaine en mouvement, Ωr , et une autre xe, Ωs . Ω Ωred
Quand Ωr bouge, la continuité des sources doit être as-
surée à l'interface entre Ωr et Ωs . Pour ne pas l'assurer, (a) (b) Maillage (c) Temps de calcul
Géométrie avec symétrie
la source est calculée sur un domaine réduit Ωr ou Ωs . complète et périodicité
Si la source est calculée sur Ωred , les termes (3b)(3c)
de la formulation en A − V ne sont plus intégrés sur Ω Fig. 1. Géométrie, maillage et temps de calcul pour la machine

et Γ, mais sur Ωred et Γred . La condition n × Hs du asynchrone

terme (3c) n'est pas nulle, car Γred n'est pas l'inni. Le
terme (3c) est simplié en une intégrale volumique.
Les fonctions de forme d'arêtes Wk sont divisées
telles que Wk = Wωk + Wγk avec Wωk pour les arêtes
dans Ωred mais pas sur Γred , et Wγk pour les arêtes (a) Hj avec la (b) Aj avec la (c) Aj avec la (d) Aj avec la
sur Γred . En remplaçant Wk par Wωk , le terme (3c) norme L2 (Ω) norme semi-norme semi-norme
est nul, parce que Wωk est nulle sur Γred , le domaine H(curl, Ω) H(curl, Ω) H(curl, Ωred )
d'intégration de (3c). En remplaçant Wk par Wγk , les
Fig. 2. Isovaleurs de B pour la machine asynchrone
termes (3b)(3c) sont nuls ensemble, parce que le do-
maine des bobines Ωb est Z inclus dans Ωred sans contact
entre Γred et Ωb , donc Wγk · Js dΩ = 0. Secondement, la meilleure source, Aj interpolée avec
Ωred la semi-norme d'H(curl, Ω), est calculée sur un domaine
Finalement, l'équation (3) de la formulation en A−V réduit Ωred , égal au stator et à l'air de la machine. Les
avec le calcul de la source sur Ωred devient telle que : isovaleurs de B sur la Fig. 2 et les temps de calculs sur
Z  
∂A
 la Fig. 1c sont similaires que la source Aj soit calcu-
νrotWk · rotA + σWk · + gradV dΩ (9a) lée sur Ω (Fig. 2d) ou Ωred (Fig. 2c). Finalement, la
Ω ∂t
Z source Aj avec la formulation en A − V peut être cal-
= rot Wωk · Hs dΩ, ∀ k = 1 · · · na . (9b) culée sur Ωred , avec des résultats similaires à Ω.
Ωred

La source Hs de l'équation (9) peut être remplacée par VI. CONCLUSION

les sources de la partie précédente, Hj ou ν0 rot Aj .
La meilleure source pour la formulation en A−V avec
des bobines non-maillées est le potentiel vecteur discré-
V. APPLICATION
tisé Aj interpolé avec la semi-norme d'H(curl, Ω). De
Une machine asynchrone est étudiée sur une version plus, il est possible de ne pas calculer la source Aj sur
non-commerciale du logiciel Flux R en régime harmo- le domaine entier Ω pour le limiter à un domaine ré-
nique avec la formulation en A − V (9)(4), en rem- duit Ωred , soit xe ou en mouvement, sans nécessité
∂ d'assurer la continuité d'Aj à l'interface Γred . Cette
plaçant par jω , avec ω la pulsation. La géométrie source ne sera pas recalculée à chaque nouvelle position.
∂t
et le maillage sont représentés sur les Fig. 1a et 1b,
et détaillés dans [3]. Sur la Fig. 1c, les temps de pré- RÉFÉRENCES
résolutions des sources avec le gradient conjugué et de
la résolution de la formulation en A − V avec un sol- [1] J.-C. Vérité, Calculation of multivalued potentials in exterior

veur direct sont représentés. Sur la Fig. 2, les isovaleurs regions, IEEE Trans. Magn., vol. MAG-23, no. 3, pp. 1881
1887, May 1987.
de B sont représentées sur le plan xy à z = 0 pour deux
[2] O. Biro, K. Preis, and K. R. Richter, On the use of the ma-
encoches du stator. gnetic vector potential in the nodal and edge nite element
Premièrement, les diérents calculs de source sont analysis of 3d magnetostatic problems, IEEE Trans. Magn.,
étudiés sur la machine asynchrone. Sur la Fig. 2a, il y a vol. 32, no. 3, pp. 651654, May 1996.

des taches anormales pour les isovaleurs de B au-dessus [3] P. Ferrouillat, C. Guérin, G. Meunier, B. Ramdane, P. Labie,

des encoches du stator. Donc, la source Hj (Fig. 2a) and D. Dupuy, Computation of source for non-meshed coils
with a-v formulation using edge elements, IEEE Transac-
donne des mauvais résultats pour les isovaleurs de B sur tions on Magnetics, En cours de publication.
le domaine Ω, alors qu'avec Aj (Fig. 2b et 2c) elle donne [4] L. K. Urankar, Vector potential and magnetic eld of current-
des bons résultats. Sur la Fig. 1c, le plus petit temps de carrying nite arc segment in analytical form, part iii : Exact
calcul pour la source Aj sur Ω est celui de l'interpolation computation for rectangular cross section, IEEE Trans.
Magn., vol. MAG-18, no. 6, pp. 18601867, November 1982.
avec la semi-norme d'H(curl, Ω), environ 1 heure. Fina-
[5] Z. Wang, T. Henneron, N. Nemitz, J.-C. Mipo, and F. Piriou,
lement, la meilleure source pour la formulation en A−V
Electromagnetic eld projection on nite element overlap-
avec des bobines non-maillées est Aj interpolée à partir ping domains, IEEE Trans. Magn., vol. 49, no. 4, pp. 1290
d'Hs avec la semi-norme d'H(curl, Ω). 1298, April 2013.

View publication stats