Polytech Projets Ge

p>. !https://forge.clermont-universite.fr/attachments/download/16442/logo_20140122173631_20140122173645.jpg!

h1=. P13A08 Modélisation et simulation d’un moteur roue

p=. !https://forge.clermont-universite.fr/attachments/download/16443/moteur_20140122195528_20140122200104.png!

p=. Polytech' Clermont-Ferrand

Génie Electrique

Projet 2013 GE5A

Sujet: P13A08

*Entreprise / Client* :

Département Génie Physique de Polytech' Clermont-Ferrand représenté par M. Lionel BATIER

*Responsable Projet* : M. Rafik SMAALI

*Tuteur industriel* : M. Pascal FICKINGER

*Tuteur technique* : M. Christophe PASQUIER

*Equipe de projet* :

GE2: Jaouad FARIS / Mohamed OUAMEUR / Thomas ROUSSEAU.

GE3: Jaouad FARIS / Mohamed OUAMEUR.

---

h2=. Plan

---

[[1. Résumé]] 

[[2. Abstract]] 

[[3. Introduction]] 

[[4. Présentation du Sujet]] 

[[5. Cahier des Charges]] 

[[6. Développement]] 

p(((. [[1. Problématiques]] 

[[2. Faisabilité]] 

[[3. Etude Théorique]] 

[[4. Démarches]] 

[[5. Solutions]] 

[[7. Travail complémentaire]] 

p(((. [[1. Etude thermique]]

[[2. Etude d'une proposition]] 

[[8. Gestion de Projet]]

p(((. [[1. W.B.S.]] 

[[2. Gantt]] 

[[9. Notes d'application]] 

p(((. [[1. sujet 1]]

[[2. sujet 2]] 

[[10. Bilan]] 

p(((. [[1. Etat d'avancement]] 

[[2. Analyse Critique]] 

[[3. Perspectives]] 

[[11. Bibliographie]]

---

h2=. Résumé

---

h2. Ce projet intitulé "Modélisation, simulation et optimisation d’un moteur roue" a pour objectif de créer et valider un modèle du moteur roue. Ce moteur a pour mission de faire avancer la voiture solaire Bélénos développée par le département Génie Physique de Polytech Clermont-Ferrand.

h2. Le projet proposé par le département Génie physique vise à contribuer au succès de la construction de la voiture Bélénos en montrant les performances du moteur roue en simulation.

h2. Proposé en 2012/2013, le projet n’a pas abouti au résultat attendu par faute de temps, puisque le moteur a été modélisé uniquement en 2D.

h2. L'objectif de cette année est de créer un nouveau modèle en 3D et optimiser le dimensionnement de la machine afin d’obtenir un rendement optimale.

h2. MOTS CLES: Energie renouvelable, Champ magnétique, Machine sans balais, Modélisation, Simulation, Force de Laplace, Puissance, Couple.

*%{color:red}Haut de page%*

--- 

h2=. Résumé

---

h2. This project entitled “Modelling, simulation and optimization of a motor-wheel” aims to create and to validate a model of the motor-wheel. This engine has for mission to make the solar car Bélénos move; it is developed by the department of Engineering Physics of Polytech Clermont-Ferrand.

h2. The project proposed by the department of Engineering Physics aims to contribute to the success of the construction of the car Bélénos by showing the performances of the motor-wheel in simulation.

h2. Proposed in 2012/2013, the project had not the result waited by, for lack of time, because the engine was only modelled in 2D.

h2. The objective of this year is to create a new model in 3D and to optimize the sizing of the engine to obtain an optimal efficiency.

h2. KEYWORDS: Renewable Energy, Magnetic field, Brushless machine, Modeling, Simulation, Laplace Force, Power, Torque.

*%{color:red}Haut de page%*

---

h2=. Introduction

---

h2. Dans le cadre d’un enjeu environnemental grandissant, notre client, le département Génie Physique a développé une voiture solaire « Bélénos ». En utilisant les énergies renouvelables, ce système d'innovation est un produit répondant strictement aux besoins de se développer de manière propre et durable. Afin de résoudre une partie des problèmes de pollution et de crise énergétique, ce projet propose une solution dans les domaines de l'électrotechnique, l'électromagnétique, l'électronique de puissance etc.

p=. !https://forge.clermont-universite.fr/attachments/download/16444/voiture_belenos_20140122194537_20140122195423.jpeg!

h2=. Figure 01: Voiture Bélénos développée par le Génie Physique de Polytech Clermont-Ferrand

*%{color:red}Haut de page%*

---

h2=. Présentation du Sujet

---

Le projet intitulé modélisation, simulation et optimisation d’un moteur roue a pour objectif de créer et valider un modèle du moteur roue. Ce moteur a pour mission de faire avancer la voiture solaire Bélénos développée par le département Génie Physique de Polytech’ Clermont-Ferrand en transformant l’énergie électrique en énergie mécanique. Afin d’améliorer les performances du moteur, les modifications des paramètres de construction seront faites après la simulation à l’aide du logiciel *COMSOL MULTIPHYSICS* basé sur la méthode des éléments finis.

*%{color:red}Haut de page%*

---

h2=. Cahier des Charges

---

h2<. +*%{color:red}Etude Théorique%*+ 

Pour leur voiture « Bélénos », le département du génie Physique a besoin d’un moteur roue « brushless » (absence de balais) à flux axial *%{color:red}[Voir figure 02 et 03]%* avec rotor disque et la génération du champ magnétique axial pour maximiser le rendement.

(L'absence de balais et de collecteur se traduit par une durée de vie plus longue qui rend une meilleure fiabilité de dispositif).

p=. !https://forge.clermont-universite.fr/attachments/download/16445/flux_20130509193030_20130509193042.jpeg!

h3=. Figure 02: Flux axial et flux radial

p=. !https://forge.clermont-universite.fr/attachments/download/16446/coupe_20131223163520_20131223163752.gif!

h3=. Figure 03: Coupe de modèle 3D mécanique du moteur roue

h2<. +*%{color:red}2-Caractéristiques Électriques du moteur%*+ 

h2. Pour pouvoir simuler et tester les différents paramètres du moteur il est demandé de respecter les caractéristiques électriques suivantes:

h2(((. - Couple 0< C <90 N.m

- Vitesse de rotation 0 < n < 1061 tr/min

- Tension d'alimentation 99 < U < 155 V

- Courant nominal 12 < I < 30 A

Toutes ces caractéristiques correspondent aux valeurs que le moteur ne doit pas dépasser lorsque la voiture franchit une pente de 15%. C'est à partir de ces valeurs maximales d’utilisation du moteur que nous réaliserons nos simulations

h2<. +*%{color:red}3-Caractéristiques Mécaniques du moteur%*+ 

Pour modéliser le moteur il faudra respecter ses caractéristiques géométriques et ses caractéristiques mécaniques.

p=. !https://forge.clermont-universite.fr/attachments/download/16447/trapeze_20130508205217_20130508205335.jpeg!

h3=. Figure 04: Forme trapézoïdale des aimants

p=. !https://forge.clermont-universite.fr/attachments/download/16448/rayon_20130508205217_20130508205248.jpeg!

h3=. Figure 05 : Positionnement des rayons

h2<. +*%{color:red}4-Objectif%*+

h2. Le projet a pour but de trouver le meilleur rendement possible avec ce type de moteur. En effet un laboratoire Australien à réussi à obtenir un rendement de 98% sur un moteur de même géométrie. Il est demandé de se rapprocher de ce rendement.

Pour ce faire, il est possible de modifier 2 paramètres sur la géométrie du moteur:

h2(((. - L'entre Rotor (en bleue)

- L'épaisseur de l'acier (en rouge) 

*%{color:red}Haut de page%*

---

h2=. Développement

---

h2. *Problématique*

h2. En raison des enjeux environnementaux actuels concernant les problèmes de pollution, la nécessité d’utiliser des énergies renouvelables, et la nécessité de faire des économies d’énergie avec des appareils à haut rendement , le département du Génie Physique a développé une voiture solaire baptisée « Bélénos » et ils veulent augmenter le rendement du moteur roue utilisé.

En utilisant les énergies renouvelables, ce système d'innovation est le produit répondant strictement aux besoins de se développer de manière propre et durable. Ce projet concerne une solution qui appliquent les technologies dans les domaines de l'électrotechnique, l'électromagnétique, l'électronique de puissance, etc.

h2. Le travail à effectuer, est une optimisation du moteur roue. Par des moyens logiciels, cette étude agit directement sur les dimensions et les caractéristiques physiques de cette machine électrique.

h2. Grâce au panneau solaire et au moteur roue, l'énergie solaire est transformée en énergie électrique puis en énergie mécanique qui se traduit par un mouvement linéaire.

---

h2. *Faisabilité*

h2. Dans des situations complexes comme celle-ci, les performances en simulation sont calculées en résolvant les équations de la magnétostatique à l’aide de méthodes numériques telles les différences finies ou les éléments finis.

h2. La méthode des éléments finis est utilisée pour résoudre des problèmes de physiques en résolvant numériquement des équations aux dérivées partielles (EDP). Ces équations peuvent donc représenter le comportement dynamique de certains systèmes physiques même très complexes qui sont continus et décrit par une équation aux dérivées partielles linéaire. Cette méthode permet donc de résoudre de manière discrète une EDP avec une solution approchée suffisamment fiable en mettant en place un algorithme mathématique. C’est cet algorithme qui recherche une solution approchée d’une EDP sur un domaine compact avec conditions aux bords (conditions limites) et/ou à l’intérieur du compact. Ces conditions aux bords permettant d'assurer existence et unicité d'une solution.

p=. !https://forge.clermont-universite.fr/attachments/download/16449/uu_20140122200819_20140122201015.png!

h3=. Figure 06. La méthode des éléments finis 

h2. L’utilisation de cette méthode implique de choisir un maillage de la structure étudiée. Il s’agit d’un découpage de l’espace selon un maillage. D’après le contexte, plusieurs formes d’un maillage sont possibles, telles que maillage carré ou triangulaire. Et il est à noter que quelle que soit la forme du maillage, plus ce maillage est resserré, plus la solution que l’on obtient sera précise et proche du résultat réel. Enfin, le calcul numérique devient plus complexe quand on resserre le maillage.

p=. !https://forge.clermont-universite.fr/attachments/download/16450/oo_20140122200819_20140122201032.png!

h3=. Figure 07. L’illustration d’un maillage

---

h2. Etude Théorique

---

h2. Pour le calcul de rendement les principales équations utilisées sont :

h2(((. *La puissance absorbée :*

h2. Dans cette phase, puisque on fixe notre étude dans le régime statique, on calcule le rendement avec le rapport de la puissance électromagnétique, sur la puissance absorbée, cette dernière est fixée, pour un courant fixé à 12A pour le régime nominal, et 30A pour avoir une idée sur les valeurs maximales, et une tension à 155V. Et si on prend en considération le coefficient de convertisseur statique (l’onduleur) qui adapte l’alimentation du moteur, on peut écrire la puissance absorbée par le moteur :

h2=. Pa = Ubatt*Ibatt* rend_onduleur

h2(((. *La puissance électromagnétique :*

h2. La puissance utile qui fait fonctionner le moteur s’exprime en fonction du couple électromagnétique et la vitesse de rotation. Elle est fixée pour une fréquence de fonctionnement normale. Pour toute l’étude statique faite, la fréquence choisie est 500Hz.

Pour calculer le couple électromagnétique, la méthodologie consiste à calculer le couple à partir de la force électromagnétique (ou la force de Lorenz ou la force mécanique), ou bien à partir de l’induction de l’entrefer.

L’expression de cette puissance est donc :

h2=. Pem=Cem*vitesse_angulaire

h2(((. *Le couple électromagnétique :*

h2. Comme toute machine électromagnétique, on note que le moteur roue est réversible, si on alimente le stator d'une dynamo excitée par une source séparée, Chaque brin actif placé dans un champ et parcouru par un courant va donc être soumis à une force électromagnétique dont le sens est donné par la règle des 3 doigts de la main droite. Les composantes tangentielles de toutes ces forces produisent un couple qui détermine la rotation du moteur. Cette machine qui reçoit de l'énergie électrique et produit de l'énergie mécanique est un moteur de puissance électromagnétique Pem.

Pour calculer ce couple, il suffit de connaitre la force électromagnétique exercée dans chaque point de brin multipliée par la distance entre ce point et l’axe du moteur roue.

p=. !https://forge.clermont-universite.fr/attachments/download/16451/trapeze_20130508205217_20130508205335.jpeg!

h3=. Figure 08. La représentation des forces électromagnétiques

h2. Avec :

h2=. !https://forge.clermont-universite.fr/attachments/download/16452/kk_20140122200819_20140122201238.png!

h2(((. *La force électromagnétique :*

h2. La force de Lorentz, ou force électromagnétique, est la force que va subir un point des fils de bobinage dans un champ électromagnétique.

C'est la principale manifestation de l'interaction électromagnétique. La force de Lorentz, appliquée dans diverses situations, induit l'ensemble des interactions électriques et magnétiques observées sur le moteur. Le sens de cette force est déterminé par la règle des 3 doigts de la main droite.

p=. !https://forge.clermont-universite.fr/attachments/download/16453/dd_20140122201600_20140122201824.gif!

h3=. Figure 09. La règle des 3 doigts de la main droite liant force, champ et courant

h2. Le champ électromagnétique exerce la force suivante sur des points des fils de bobinage

p=. !https://forge.clermont-universite.fr/attachments/download/16454/dd45_20140122204316_20140122204443.png!

h2. E et B sont respectivement le champ électrique et le champ magnétique pris au point. On peut distinguer deux contributions à cette force :

F_elec=qE : qui est la force électrique ;

F_mag=qv ⋀ B : qui est la force magnétique.

h2. En regardant cette force dans un repère orthonormé, on peut analyser le sens de F comme ci-dessous :

p=. !https://forge.clermont-universite.fr/attachments/download/16455/dv_20140122201600_20140122201912.png!

h3=. Figure 10. Représentation des forces électromagnétiques

h2. En outre, dans un niveau macroscopique, cette force porte un autre nom c’est celui de la force de Laplace qui est la résultante de l’action de la force de Lorentz sur tous les points chargés dans notre cas.

Tant la force de Lorentz que la force électromagnétique, elles peuvent être appelées « La force de Laplace ».

Ou plus simplement dans notre cas, avec une densité de courant J, qu’on modifié pour avoir soit 12 ou soit 30A, on peut écrire la force : *F=J*B*

---

h2. *Démarches*

---

h2. La démarche de ce projet se divise en trois phases. La première phase consiste à construire une modélisation en 3D sous le logiciel COMSOL avec des géométries basiques dont on ne respecte pas forcément la forme réelle des aimants, la deuxième phase est celle d’une modélisation partielle pour avoir des idées générales sur les performances, et la dernière phase est celle d’une modélisation complète, avec des simulations des géométries à chaque phase.

En parallèle on travaille sur un programme d’optimisation des paramètres géométriques sous le logiciel MATLAB, avant de mettre en évidence les performances trouvées.

h2. L’algorithme d'optimisation utilisé est un algorithme génétique, c’est un algorithme évolutionniste. Le but est d'obtenir une solution approchée au problème d'optimisation, notamment parce qu'il n'existe pas de méthode exacte sous COMSOL pour le résoudre en un temps raisonnable. Ce type d’algorithme utilise la notion de sélection naturelle et l'applique à une population (nombre d’itérations) de solutions potentielles au problème donné. La solution est approchée par « bonds » successifs, comme dans une procédure de séparation et évaluation, à ceci près que ce sont des formules qui sont recherchées et non plus directement des valeurs.

h2(((. On peut schématiquement résumer cette partie par le système bouclé suivant :

p=. !https://forge.clermont-universite.fr/attachments/download/16456/dzdfff_20140122201600_20140122203401.png!

+*%{color:red}1-Modélisation et Simulation%*+ 

h2. La modélisation et la simulation du moteur seront faites dans un environnement logiciel qui utilise la méthode des éléments finis. La solution proposée est le logiciel COMSOL MULTIPHYSICS.

p=. !https://forge.clermont-universite.fr/attachments/download/16457/kl_20140122213739_20140122215831.jpeg!

h3=. Figure 11. Le logiciel COMSOL MULTIPHYSICS

h2. COMSOL MULTIPHYSICS (anciennement appelé : Toolbox 1.0, Femlab) est un logiciel de simulation numérique développé et commercialisé par la société COMSOL. Ce logiciel est un outil de résolution d’EDP basé sur la méthode des éléments finis. Sa particularité est de disposer d’une base de données d’équations qui permet de modéliser et simuler différents phénomène physiques (multi-physiques) comme le comportement d’un fluide arrivant à grande vitesse sur un obstacle, la déformation d’une structure métallique ou encore les performances d’une génératrice. COMSOL MULTIPHYSICS est multiplateforme, il peut fonctionner sous Windows, Mac, GNU-Linux etc. Il utilise une interface graphique qui nous permet de travailler en plusieurs dimensions (1D, 2D ou 3D) et contient la plupart des équations. Cependant, on peut également entrer manuellement des EDP spécifiques. Il permet également de faire des couplages entre différents domaines comme par exemple entre les domaines électrique, magnétique et/ou thermique, etc. Grâce à l’interface graphique, l’observation de différentes grandeurs (couple, tension, force etc.) devient visuelle.

h2. Afin de modéliser au mieux notre moteur nous avons procédé par étapes. Premièrement, on modélise une partie du moteur en 2D. Ainsi on prend en main le logiciel et on voit les résultats auxquels on doit s’attendre en 3D. Puis on modélise une partie du moteur en 3D.

La modélisation du moteur, comme pour tout autre système, consiste en 4 étapes :

h2(((. - Définition de l’étude voulue

- La création de la géométrie du moteur

- La définition des propriétés des matériaux constituant le moteur

- La définition des limites de la géométrie créée

h2. Une fois que la partie du moteur modélisée, il faut la simuler, c’est-à-dire qu’il faut résoudre les équations qui définissent le problème pour avoir un résultat significatif. Cette étape de simulation se fait en 3 étapes :

h2(((. - Le maillage de la structure

- La définition du type d’analyse et du solveur

- L’analyse des résultats

+*%{color:red}2-Optimisation%*+ 

h2. Dans cette partie l’objectif consiste à optimiser la géométrie du moteur en régime statique, afin de pouvoir délivrer une géométrie optimale permettant d’avoir le meilleur rendement de puissance possible et cela en modifiant des paramètres géométriques.

h2. Etant donné que nous devons faire varier différents paramètres ; l’entre-rotor et l’épaisseur d’acier, il est nécessaire d’utiliser une méthode qui nous permet d’avoir un couple optimal de ces deux paramètres, correspondant au rendement maximal du moteur. La solution proposée est l’utilisation du logiciel MATLAB.

p=. !https://forge.clermont-universite.fr/attachments/download/16458/zd_20140122213739_20140122215857.gif!

h3=. Figure 12. Le logo MATLAB 

h2. MATLAB (« matrix laboratory ») est un langage de programmation de quatrième génération et un environnement de développement, il est utilisé à des fins de calcul numérique. Développé par la société The MathWorks?, ce logiciel permet la manipulation des matrices, l’affichage des courbes et des données, de mettre en œuvre des algorithmes, la création des interfaces utilisateurs, et peut s’interfacer avec d’autres langages comme le C, C++ et Java. Les utilisateurs de MATLAB sont de milieux très différents comme l’ingénierie, comme le cas de notre projet les sciences et l’économie dans un contexte aussi bien industriel que pour la recherche. MATLAB peut s’utiliser seul ou bien avec des toolbox (« boîte à outils »).

L’existence d’une interface entre MATLAB et COMSOL, ainsi que notre connaissance de MATLAB grâce à son utilisation dans notre formation, nous ont poussé à faire le choix d’utiliser ce dernier.

---

h2. *Solutions*

h2. +*%{color:red}1-Modélisation%*+ 

h2((.  *+1.1-Géométrie Basique+*

h2. La géométrie basique servira de base pour la compréhension des simulations en 3D

p=. !https://forge.clermont-universite.fr/attachments/download/16459/geobasique_20130508212731_20130508213049.jpeg!

h3=. Figure 13. modélisation de la géométrie basique

h2((.  *+1.2-Géométrie partielle réelle+*

h2. La géométrie partielle réelle servira de base pour toutes les simulations lorsque l'on devra jouer sur les différents paramètres (entre-rotor et épaisseur de l'acier) pour trouver les caractéristiques optimales

p=. !https://forge.clermont-universite.fr/attachments/download/16460/geopartielle_20140123112214_20140123112236.gif!

h3=. Figure 14. modélisation de la géométrie partielle réelle

h2((.  *+1.3-Géométrie complète+*

h2. La géométrie complète servira à comparer les résultats entre les paramètres initiaux (entre-rotor = 9mm ; épaisseur de l'acier = 6mm) et les paramètres optimaux

p=. !https://forge.clermont-universite.fr/attachments/download/16461/geocomplte_20130508212731_20130508213113.jpeg!

h3=. Figure 15. modélisation de la géométrie partielle

h2. +*%{color:red}3-Simulation%*+ 

h2((.  *1.1-Géométrie Basique+*

h2. L’étude faite comme 1ère étape de la prise en main du logiciel. On a représenté ici la densité de flux magnétique, les lignes du champ magnétique.

p=. !https://forge.clermont-universite.fr/attachments/download/16462/dez_20140122201600_20140122203932.png!

h3=. Figure 16. une partie du moteur en 2D

h2((.  *+1.2-Géométrie partielle réelle+*

h2. On retrouve la même saturation qu’on 2D. On retrouve également les lignes du champ mais cette fois, moins nette qu’en 2D. Ces aberrations sont dues à la taille du maillage qu’on a définie pour la simulation. Cette taille de maillage est la meilleure qu’on puisse avoir étant donnée la quantité de mémoire que cela demande à l’ordinateur utilisé.

p=. !!

h3=. Figure 17. une partie du moteur en 3D