Modélisation d'un aérogénérateur double étoile
Modélisation d'un aérogénérateur double étoile
THEME
Présenté par :
YOUNES SAHRAOUI
BELMOSTEFA ISSAME
I
a). Equation de tension ................................................................................................................... 18
b). Equations des flux ...................................................................................................................... 18
c). Expression du couple électromagnétique ...................................................................................... 19
d). Equation mécanique .................................................................................................................. 20
II.3 TRANSFORMATION DE PARK ............................................................................................................ 20
II.4 Choix du référentiel........................................................................................................................... 22
II.4.1. Référentiel lié au stator ............................................................................................................. 22
II.4.2. Référentiel lié au rotor............................................................................................................... 22
II.4.3. Référentiel lie au champ tournant ............................................................................................. 22
II.4.4. Modèle de la MASDE dans le repère de Park (d,q) lie au champ tournant ................................. 22
II.5 REPRESENTATION D’ETAT DE LA MASDE ............................................................................................ 23
II.5.1. Puissance absorbée et couple électromagnétique...................................................................... 25
II.6 MODELISATION DE LA TURBINE EOLIENNE ........................................................................................ 26
II.6.1. Modélisation de la vitesse du vent ............................................................................................. 27
II.6.2. Modélisation de la turbine ......................................................................................................... 27
II.6.3. Génération du couple éolien ...................................................................................................... 28
II.6.4. Modèle du multiplicateur .......................................................................................................... 29
II.6.5. Equations dynamique de l’arbre de la génératrice ..................................................................... 29
II.7 CONCLUSION..................................................................................................................................... 30
CHAPITRE III :LA COMMANDE VECTORIELLE DE LA GASDE ............................................................................. 32
III.1. INTRODUCTION ................................................................................................................................. 32
III.2. COMMANDE VECTORIELLE ................................................................................................................ 32
III.2.1. Choix d’orientation du flux ......................................................................................................... 33
III.3. DIFFERENTES METHODES DE LA COMMANDE VECTORIELLE .............................................................. 34
III.3.1. Méthode de commande directe ................................................................................................. 34
III.3.2. Méthode de commande indirecte .............................................................................................. 34
III.4. COMMANDE VECTORIELLE INDIRECTE SANS REGLAGE DE VITESSE .................................................... 34
III.5. APPLICATION DE LA COMMANDE INDIRECTE SUR LA MASDE ............................................................. 37
III.5.1. Commande vectorielle indirecte avec régulation de vitesse ....................................................... 37
III.6. CONTROLE DE LA TURBINE ................................................................................................................ 38
III.6.1. Maximisation de la puissance éolienne ...................................................................................... 40
III.6.2. Modélisation de l’onduleur MLI ................................................................................................. 41
Commande de la tension par MLI sinus-triangle .................................................................................... 44
III.7. RESULTATS DE SIMULATION ET INTERPRETATION.............................................................................. 45
Interprétationdesrésultatsobtenus ........................................................................................................ 47
III.8. CONCLUSION..................................................................................................................................... 47
II
CONCLUSION GÉNÉRALE .................................................................................... Error! Bookmark not defined.
Bibliographie ................................................................................................................................................ 50
III
الملخص
الهدف من تطبيق تقنية، ( في نظام رياح متذبذب السرعةGASDE) يتناول هذا العمل التحكم الشعاعي غير مباشر لمولد مزدوج النجم غير متزامن
ومن هنا قمنا في الجزء األول بمنقاشه عامة، ( المطبقIFOC) ( وتحسين أداء التحكم غير المباشرGASDE) التحكم هو تنظيم السرعة التي تولدها
(ومن اجلPWM) . (التي يتم تغذيتهاعن طريق محوالت الجهدMASDE ) وأما فالجزء الثاني فتطرفناإلىنمذجة الماكنة رياضيا، حول طاقة الرياح
و في األخير توضيح فعالية التحكم غير المباشر، تقليدي لتنظيم السرعةPI و عليه فقد تم تطبيق التحكم الشعاعي غير مباشر وتم استخدام منظم، تحكم فيها
.من خالل نتائج المحاكاة
كلمات مفتاحية
.PI منظم، )FOC( تحكمالشعاعي غير مباشر، ، نظام رياح متغير السرعة،مولد مزدوج النجم غير متزامن
Résumé
Ce travail traite la commande vectorielle d'un générateur asynchrone à double étoile (MASBE) dans un système de
vent oscillant. L'objectif de l'application de la technologie de contrôle est de réguler la vitesse nominale (MASDE) et
d'améliorer les performances du couple appliqué (FOC), et à partir de cela dans la première partie, nous avons une
discussion générale sur l'énergie éolienne, et dans la deuxième partie, nous référons à une description de (MASDE), qui
est obtenue grâce à des transformateurs de tension. (MLI) Afin de modéliser, MASDE, nous avons appliqué un contrôle
vectoriel dans la direction du flux, et un régulateur PI classique a été utilisé pour réguler la vitesse. L'efficacité du
contrôle indirect de la transmission a été démontrée par Pendant les résultats de la simulation.
Mots clés
Génératrice asynchrone à double étoile, système éolienne à vitesse variable, Algorithme MPPT , Onduleur de
tension, commande vectorielle (FOC), régulateur PI .
Abstract
This work deals the vector control of a double-stator asynchronous generator (DSAG) in an oscillating wind
system. The objective of the application of the control technology is to regulate the rated speed (DSAG) and to improve
the performance of thevector control (FOC), and from this In the first part, we have a general discussion about wind
energy, and in the second part, we refer to a description of (DSAG), which is obtained through voltage transformers.
(PWM) In order to model, MASDE we applied vector control through the direction of flow, and a recursive PI
classiquewas used to regulate the speed. The effectiveness of indirect transmission control was demonstrated by During
the simulation results.
Key words
Double star asynchronous generator, Variable speed wind system, Technical optimization, MPPT, Voltage inverter
vector control, Algorithm vector control (FOC), PI regulator.
I
LISTE DES FIGURES
Fig.I.6. Principauxélémentsd’unsystèmeéolien.……………………………………..........………..8
Fig. I.14. Flux d’air sur un profil de pale " stall ".……………………………………………….12
[Link]èmeéolienbasésurlamachineasynchronedoublealimentation………………….……13
Fig. [Link] aérodynamique en fonction de ratio de vitesse 𝛌pour une éolienne de 1.5
MW. ………………………………………………………………………………………………….29
II
Fig. II.6. Schéma bloc du modèle de la turbine. …………………………………………………30
Fig. III.7. Schéma bloc du modèle de la turbine avec le régulateur PI. ………………………...39
Fig. III.9. Alimentation de la MASDE par deux onduleurs de tension triphasé[Link] commande
des courants. ………………………………………………………………………………………...42
III
INTRODUCTION GENERALE
La production d'énergie sera un défi important dans les années à venir. En effet, les besoins
énergétiques des sociétés industrialisées ne cessent d'augmenter. De plus, les pays en développement
auront besoin d'une quantité croissante d'énergie pour atteindre leurs objectifs.
Développement. Aujourd'hui, une grande partie de la production mondiale d'énergie est
garantie. À base de sources fossiles La consommation de ces ressources entraîne des émissions de
gaz à effet de serre, et donc une augmentation des polluants un risque supplémentaire provient de la
une consommation excessive de ressources naturelles épuise les réserves de ressources naturelles de
ce type.
L'énergie de manière dangereuse pour les générations futures Les énergies renouvelables telles
que l'éthanol, l'énergie solaire, la géothermie et l'énergie résiduelle (énergie issue de la combustion
des déchets) sont des ressources utilisables indéfiniment sans polluer, permettant aux générations
futures de vivre sur la même planète qu'aujourd'hui.
Parmi ces énergies renouvelables l'énergie éolienne, qui est une forme d'électricité propre et
abondante. Aujourd'hui, l'avancement de la technologie éolienne à fréquence constante et à vitesse
variable est au centre d'une transition réussie et globale vers cette énergie.
Dans ce but, nous présentons nos recherches sur l'énergie éolienne, et c'est ce qui a poussé les
chercheurs à mener des études afin d'améliorer l'efficacité de la conversion électromécanique et de
l'énergie produite.
Plusieurs machines multi-étages sont apparues au siècle dernier que les chercheurs se sont
intéressés à étudier et à contrôler,La machine asynchrone double étoile (MSADE), où notre objectif
est de savoir comment contrôler les régulateurs PI conventionnels .
Le premier chapitre
Est consacré à des statistiques sur le développement de l’énergie Éolienne puis nous allons
intéresser à la construction et le fonctionnement de ces aérogénérateurs, suivi par l’étude des
différentes types des éoliennes et des génératrices utilisées et les éléments principaux dans la chaine
de conversionélectromagnétique ,enfin la présentation les avantages, les inconvénients l’énergie
Éolienne.
Le deuxième chapitre
Est consacré à la modalisation de la MASDE alimentée par deux Onduleurs de tension, en
utilisant la transformation du Park (passage triphasé-biphasé).
1
Le troisième chapitre
Est consacré la présentation du principe et l’application de la commande vectorielle à la GASDE. Les
résultats de simulation avec le régulateur classique.
2
CHAPITRE I
I.1. INTRODUCTION
L'énergie éolienne est énergies renouvelables largement utilisées comme alternative aux
énergies fossiles, c'est une énergie abondante et renouvelable que l'on retrouve dans toutes les
régions, mais son abondance varie d'un endroit à l'autre. C'est une énergie propre et renouvelable qui
ne produit pas d'émissions telles que des gaz à effet de serre pendant l'exploitation, et elle nécessite
des surfaces variables selon la taille de l'usine et le type de tours utilisées.
Réfléchir à l'utilisation de l'énergie éolienne aujourd'hui est devenu très important afin de la
développer de plus en plus, en raison de son impact positif sur l'environnement en réduisant le
pourcentage de pollution, car elle est considérée comme respectueuse de l'environnement.
L'objectif de ce chapitre est de présenter un aperçu général sur l’énergie éolienne, son
historique et la croissance de l’exploitation de cette énergie, en présentant la description de l’éolienne
et son fonctionnement puis les différents types des éoliennes et leur variété selon
la position de l'arbre et la situation par rapport au vent ,le coefficient de puissance et les techniques de
régulation de la puissance en Pitch ou Stall Control. En termina ce chapitre par les avantages et les
inconvénients de l’énergie éolienne.
I.2. LA SITUATION ACTUELLE DE L’ENERGIE EOLIENNE
Le GWEC (Global Wind Énergie Council) a publié le 4 avril son Global Wind Report 2022,
Voici ci-après quelques grandes données chiffrées de l’organisme défendant les intérêts de la filière
éolienne au niveau international.[w1].
En 2021, la filière éolienne a connu, en matière de nouvelles installations, « la 2e meilleure année de
son histoire » au niveau mondial (après 2020) : 93,6 GW de nouvelles capacités éoliennes ont été
déployées, portant la puissance du parc mondial en fin d’année à 837 GW (soit 12% de plus qu’à fin
2020).[w1].
Sur ce total, le GWEC répertorie 72,5 GW de nouvelles éoliennes terrestres, soit une chute de 18%
par rapport au niveau de 2020 en raison d'un ralentissement dans les deux principaux pays de
développement de la filière : la Chine et les États-Unis (la fin de tarifs d’achats garantis en Chine en
2021 a causé une chute de 39% des installations dans ce pays par rapport à l'année 2020). [w1].
3
L’éolien offshore a pour sa part connu de loin la meilleure année de son histoire, avec 21,1 GW de
nouvelles installations, soit près de « trois fois plus que l’année précédente » (portant à 57 GW la
puissance cumulée des éoliennes offshore dans le monde à fin 2021) [w1].
4
I.4. ENERGIE EOLIENNE EN ALGERIE
Les ressources éoliennes en Algérie varient considérablement d'un endroit à l'autre. Ceci est
principalement dû à une topographie et un climat très diversifiés. En effet, notre vaste pays est divisé
en deux régions géographiques distinctes : la Méditerranée Nord se caractérise par un littoral de 1200
km et une région montagneuse, par l'Atlas Tell et l'Atlas Entre eux se trouvent des plaines et des
plateaux au climat continental. De son côté, le Sud est caractérisé par le Sahara. [w2].
Les cartes des vents développées par plusieurs chercheurs montrent que le Sud-Est est
caractérisé par des vitesses plus élevées que le Nord, notamment au Sud-
[Link] chercheurs du premier atlas éolien algérien établi avec la vitesse la plus
élevée de l'ordre de grandeur du 6 m/s et situé dans la région d'Adrar. Ces résultats, réalisés sur 10
ans de mesure, sont basés sur des cartes de vent développées par des chercheurs DER. «Centre de
Développement des Energies Renouvelables »
Les travaux effectués récemment, ont mis en évidence l’existence de sites ventés dans d’autres
régions du Sud. Les régions de Tamanrasset, Djanet et In Salah disposeraient d’un potentiel éolien
exploitable. Le gisement éolien en Algérie est donc toujours en cours d’évaluation.
Après avoir fait des tests sur 63 sites étudiés à l’intérieur du pays et 24 dans les pays voisins, les
résultats des chercheurs du centre du Développement sont comme suite:
On constate que le printemps est la saison la plus active suivie par l’été au sud et l’hiver au nord,
avec un vent stable.
La vitesse annuelle du vent varie du 1.2 à 6.3 m/s à 10 m de hauteur.
Les sites situés dans le sud de l'Algérie sont plus élevés avec un maximum enregistré à Adrar
de 6,3 m/s, de Hassi Ramel avec 6,1 m/s.
La densité annuelle varie selon les régions, la plus grande valeur est enregistrée à Hassi R'mel avec
295w/𝒎𝟐 .
En conclusion, on peut classer la carte des vents en Algérie selon les zones suivantes
(ClassementPNL« Pacific Northwest Laboratory »):[w2].
5
Fig.I.2. Carte des vents de l’Algérie à 10 m d’altitude (m/s).
Ce sont des machines appelées aérogénérateurs ou plus communément éoliennes, ce sont des
systèmes permettant de récupérer l’énergie cinétique du vent et de la convertir en énergie mécanique
pour produire le plus souvent de l’électricité [Hoc13], par l’intermédiaire d’une génératrice.
Cette conversion se fait en deux étapes :
Au niveau de la turbine, qui extrait une partie de l’énergie cinétique du vent disponible pour la
convertir en énergie mécanique.
Au niveau de la génératrice, qui reçoit de l’énergie mécanique et les convertis en énergie électrique
qui la transmit ensuite sur les réseaux électrique. [LOU16].
Il existe deux principaux types d'éoliennes qui se diffèrent essentiellement dans leur organe
capteur d’énergie à savoir l’aéroturbine. En effet, selon la disposition de la turbine par rapport au sol
on obtient :[Lat12]
Eoliennes à axe vertical.
Eoliennes à axe horizontal.
L’axe de rotation de ce type d’aérogénérateur est vertical par rapport au sol et perpendiculaire à
la direction du vent. Ce type de turbine peut recevoir le vent de n’importe quelle direction, ce qui
rend inutile tout dispositif d’orientation. Le générateur et le multiplicateur sont disposés au niveau du
sol, ce qui est plus simple et donc économique. La maintenance du système est également simplifiée
du fait qu’elle se fait au sol. La majorité des éoliennes à axe vertical tourne à faible vitesse, ce qui est
très pénalisant dans les applications de génération d’électricité avec connexion au réseau car la boite
de vitesse doit permettre une importante démultiplication. Le faible rendement et la quantité du vent
réduite qu’elles reçoivent au niveau du sol, constituent les principaux inconvénients de ce type
d’éolienne. La figure (I-4) montre un exemple d’éoliennes à axe vertical.[CHE15].
7
b) Eoliennes à axe horizontal
Une turbine à axe de rotation horizontal demeure face au vent, comme les hélices des avions et
des moulins à vent. Elle est fixée au sommet d’une tour, ce qui lui permet de capter une quantité plus
importante d’énergie éolienne. La plupart des éoliennes installées sont à axe horizontal. Ce choix
présente plusieurs avantages, comme la faible vitesse d’amorçage (démarrage). Toutefois, la boite de
vitesses et la machine électrique doivent être installées en haut de la tour, ce qui pose des problèmes
mécaniques et économiques. Selon son nombre de pales, l’éolienne à axe horizontal est dite mono-
pale, bipales, tripales ou multiples.[CHE15].
Lafigure(I.6)montrelesprincipauxélémentsd’unsystèmeéolienquisont
8
Le rotoretses pales;
Lemoyeu;
L’arbre;
La coquilledelanacelle;
Laboitedevitesseou multiplicateur;
Lefreinagedesécurité;
Lagénératrice;
Lemât outour;
Lanacelle.
I.7. FONCTIONNEMENT D’UN SYSTEME EOLIEN
Tout d'abord examinons ses principes de base si le vent qui souffle peut faire tourner les ailes
nous obtiendrons de l'électricité par le générateur qui est raccordé mais comment le vent fait-il
tourner les ailes ! regardant de plus près la lame, la lame dispose d'un grand profil aérodynamique
formé de plusieurs tailles et forme de la base à la pointe.(Fig.I.7)
La technologie simple de la structure aérodynamique fait tourner les pales de l'éolienne cela
signifie qu'une force de poussée est créé par le déplacement d'un fluide sur la surface d'une aile de
cette façon l'éolienne atteint la rotation normale à laquelle nous sommes habitués.(Fig.I.8)
Tout comme dans un train en mouvement les choses se passent de façon relative la pâle de
l'éolienne en mouvement subit également le vent de façon relative pour la palle en mouvement la
vitesse relative du vent et celle indiquée sur l'illustration par conséquent la palle de l'éolienne est
positionné de manière incliné afin de s'aligner avec la vitesse relative du vent, cetterelation est donner
par :
𝑉𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑒 = 𝑉𝑣𝑒𝑛𝑡 − 𝑉𝑝𝑎𝑙𝑒
(I.1)
On note que :
𝑽𝒘𝒊𝒏𝒅 =𝑽𝒗𝒆𝒏𝒕
𝑽𝒃𝒍𝒂𝒅𝒆=𝑽𝑷𝒂𝒍𝒆
9
[Link] vitesse relative du vent.
Au final, on va rassembler toutes ces étapes évoquées précédemment dans cette figure, d'où on
note que tous les composants sont connectés au final avec une génératrice c'est le cas utilisé dans
notre travail. Turbine
R
Vent Multiplicateur
Génératrice
Pmec
Arbre lent
v
t Arbre rapide
mec
L’objectif de la régulation mécanique est double, d’une part de protéger l’éolienne contre le vent fort
et d’autre part délimiter la puissance. En effet, la turbine éolienne est dimensionnée pour fournir une
puissance nominale pn à une vitesse de vent nominale Vn , au-delà de cette vitesse, les paramètres de
la turbine doivent évoluer afin de fournir la puissance nominale et de ne pas produire au-delà d’une
vitesse maximale Vm qui pourrait endommager la turbine. [ZER15]
On peut définir quatre zones de fonctionnement :
Zone I : où P = 0, la turbine ne fonctionne pas.
Zone II : Dans laquelle la puissance fournie surl'arbre dépend de la vitesse du vent.
Zone III : Où la vitesse de rotation est maintenueconstante et où la puissance P fournie reste égale à
Pn.
Zone IV : Dans laquelle le système de sûreté defonctionnement arrête le transfert de l’énergie.
10
FigI.11. Zones de fonctionnement d’une éolienne à vitesse variable.
Pour contrôler la puissance aérodynamique capturée et la limiter lorsque la vitesse du vent est très
élevée, on peut utiliser l’une des deux techniques de contrôle connues surtout dans les grandes
turbines sous les noms de "Pitch Control" et "Stall Control"[TOU 10].
C’est un contrôle par variation de l’angle de calage des pales (pas des pales) autour de leur axe
longitudinal en fonction de la vitesse du vent et de la puissance active à extraire.
En agissant par le Pitch Control, le coefficient de puissance𝐂𝐩 ne varie plus en fonction de λmais en
fonction de la vitesse du vent à vitesse de rotation constante, ce qui va permettre de maintenir la
puissance constante dans la zone III à sa valeur maximale:[TOU10].
On peut noter que pour modifier l’angle de calage dans le Pitch Control dans un système de
régulation automatique de vitesse de rotation on a besoin d’un régulateur électronique avec des
masselottes centrifuges ou d’un régulateur électronique avec des servomoteurs, ce qui va permettre :
Le démarrage à une vitesse du vent 𝐕𝐝 réduite.
La régulation de la vitesse angulaire de la turbine Ω𝑡 pour une vitesse du vent supérieure à la vitesse
nominale𝐕𝐧 .
L’optimisation du régime de conversion de l’énergie quand la vitesse du vent évolue entre 𝐕𝐝 et 𝐕𝐧 .
La protection de l’éolienne contre les grandes vitesses du vent (>𝐕𝐦𝐚𝐱 ) par la mise en drapeau des
pales sans nécessité d’un système puissant freinage pour les arrêts d’urgence.
11
Fig.I.12. Variation de l'angle de calage Fig.I.13. Influence de l’angle de calage
La plupart des éoliennes connectées au réseau électrique nécessitent une vitesse de rotation fixe
pour des raisons de cohérence de fréquence avec le réseau. Le système de limitation de vitesse le plus
simple et le moins coûteux est un système de limitation naturelle (intrinsèque à la forme de la pale)
dit "stall". Il utilise le phénomène de décrochage aérodynamique. Lorsque l’angle d’incidence devient
important, c’est à dire lorsque la vitesse du vent dépasse sa valeur nominale 𝑉𝑛 , l’aspiration créée par
le profil de la pale n’est plus optimale ce qui entraîne des turbulences à la surface de la pale et par
conséquent une baisse du coefficient de puissance. Ceci empêche alors une augmentation de la
vitesse de rotation.[Fré03]
Fig. I.14. Flux d’air sur un profil de pale " stall ".
Ce système est simple et relativement fiable mais il manque de précision car il dépend de la
masse volumique de l'air et de la rugosité des pales donc de leur état de propreté. Il peut, dans
certains cas, être amélioré en autorisant une légère rotation de la pale sur elle-même (système "stall
actif") permettant ainsi de maximiser l’énergie captée pour les faibles vitesses de vent. Pour les fortes
vitesses de vent, la pale est inclinée de façon à diminuer l’angle de calage et renforcer ainsi l’effet
"stall" de la pale. La répercussion des variations de vitesse de vent sur le couple mécanique fournie
par l’éolienne est ainsi moins importante.[Fré03].
12
I.9.1 MachineSynchrone
Cetypedemachine(MS)esttrèsutilisédanslaplupartdesprocédéstraditionnelsdeproductiondel’élect
ricité,notammentdansceuxdetrèsgrandepuissance(centralesthermique,hydrauliqueounucléaires).Lesgé
nérateurssynchronesutilisésdansledomaineéolien,ceuxde500KWà2MWsontbienpluschersquelesgénér
ateursasynchronesdelamême taille. De plus, lorsque ce type de machine est directement connecté au
réseau, sa vitesse
derotationestfixeetproportionnelleàlafréquenceduré[Link]équencedecettegranderigiditédelacon
nexiongénérateur-
réseau,lesfluctuationsducouplecaptéparl’aérogénérateursepropagentjusqu’àlapuissanceélectriqueprod
uite.C’estpourquoilesmachinessynchronesnesontpasutiliséesdanslesaérogénérateursdirectementconne
ctésauré[Link]éeslorsqu’ellessontconnectéesauréseauparl’intermédiairedeconv
ertisseursde puissance.
Ledéveloppementdesmatériauxmagnétiquesapermislaconstructiondemachinessynchrones à
aimants permanents à des coûts qui deviennent compétitifs. Les machines de ce typesont à grand
nombre depôles et permettent de développer descouples mécaniquesconsidé[Link]
figuresuivanteillustreun systèmebasésur la machine synchrone.
[Link]èmeéolienbasésurlamachine synchrone.
I.9.2 Machinesasynchronesàdoublealimentation(MADA)
13
fréquence(électronique de puissance) peut être réduit. Ces machines sont un peu plus complexes que
desmachines asynchrones à cage avec lesquelles ont en commun la nécessité d’un multiplicateur
devitesse. Leur robustesse est légèrement diminuée par la présence de système à bagues et balais,
maislebénéficedufonctionnementàvitessevariableestunavantagesuffisantpourquedetrèsnombreuxfabri
cants utilisent cetypedemachines.
Fig.I.16. Systèmeéolienbasésurlamachineasynchronedoublealimentation.
14
I.9.3 Machineasynchronedoublestator(doubleétoile)
Avec:
𝑃𝑚g:Puissancetransmise parlagénératrice;
Lavitesseétantliée aunombredepairesde pôles parcette relation
fi𝑠=𝜔⁄𝑃(I.13)
A une faible vitesse correspond un nombre de paires de pôles élevé. D’autre part, on a unstator de
forte puissance correspondant à une vitesse élevée et donc à faible nombre de paires depôlespourles
vitesses de vent élevées.
15
I-10 AVANTAGES ET INCONVENIENTS DE L’ENERGIE EOLIENNE
I-11 CONCLUSION
Dans ce chapitre, on a exposé, en première partie quelques généralités sur l’énergie éolienne, son
historique et son évolution. Puis, une description d’éolienne, son fonctionnement et un bref aperçu
sur les différents types des éoliennes. On a présenté la variation de la puissance extraite du vent en
fonction de la vitesse de rotation par rapport à la vitesse du vent. On a donné aussi le principe des
deux techniques de réglage de puissance "Pitch Control" et "Stall Control". Les machines électriques
utilisées dans la conversion aérodynamique sont aussi présentées avec quelques spécifications
propres liées aux structures d’alimentation différentes .A la fin de ce chapitre, on a mentionné
quelques avantages liés à la croissance de l’énergie éolienne et quelques inconvénients liés à
l’empêchement de son avancement.
16
Compte tenu des recherches récentes sur la machine asynchrone double étoile (MASDE) et vu ses
avantages, il est très intéressant de pouvoir étudier cette dernière qui sera l’un des objectifs du chapitre suivant
où on décrit le modèle de simulation du système développé.
17
CHAPITRE II
II.1 INTRODUCTION
16
Fig. II.1-Représentation des enroulements de la MASDE.
Les grandeurs relatives à la première étoile (stator 1) sont notées par l'indice s 1
, tandis que celles correspondant à la deuxième étoile (stator 2) sont notées par l'indice
s 2 (stator 2). Les phases de la première étoile sont s (a1,), s (b1,) s (c1,) et les phases
de la deuxième étoile sont s (a2,), s (b2,) s (c2,). Les phases du rotor sont r (a,), r (b,) r
(c) (c,) L'angle de séparation entre les deux étoiles est appelé. La position du rotor
(phase r (a,)) par rapport à l'étoile 1 (phase s (a1,)) est représentée par 1. θ 2 La
position du rotor par rapport à l'étoile 2 est déterminée par les équations suivantes.
𝜃1 = 𝛺𝑚𝑒𝑐 𝑡 + 𝜃0 (II. 1)
𝜃2 = 𝜃1 − 𝛼 (II. 2)
Tels que :
𝛺𝑚𝑒𝑐 : Vitesse mécanique du rotor.
𝜃0 : Position initiale du rotor parrapport à l’étoile1.
Les deux enroulements statorique sont alimentés chacune par un système
triphasé équilibré des courants entraînant la création d’un champ tournant dans le long
de l’entrefer.
17
a). Equation de tension
La combinaison de la loi d’Ohm et de la loi de Lenz permet d’´écrire les
relations suivantes
𝑑
[𝑉𝑠1 ] = [𝑅𝑠1 ][𝐼𝑠1 ] + [𝛷 ] ( II. 3)
𝑑𝑡 𝑠1
𝑑
[𝑉𝑠2 ] = [𝑅𝑠2 ][𝐼𝑠2 ] + [𝛷𝑠2 ] (II. 4)
𝑑𝑡
𝑑
[𝑉𝑟 ] = [𝑅𝑟 ][𝐼𝑟 ] + [𝛷𝑟 ] (II. 5)
𝑑𝑡
Tell que :
[R S1 ],[R S2 ],[R r ]: Matrices des résistances statoriques (étoile 1 et 2) et rotoriques
respectivement .
Les vecteurs des tensions statorique, courant et flux sont définit par :
[𝑉𝑠1 ] = [ 𝜈𝑎𝑠1 , 𝜈𝑏𝑠1 , 𝜈𝑐𝑠1 ]𝑇
[𝐼𝑠1 ] = [𝑖𝑎𝑠1 , 𝑖𝑏𝑠1 , 𝑖𝑐𝑠1 ]𝑇 (II.6)
[𝛷𝑠1 ] = [𝜙𝑎𝑠1 , 𝜙𝑏𝑠1 , 𝜙𝑐𝑠1 ]𝑇
[𝑉𝑠2 ] = [ 𝜈𝑎𝑠2 , 𝜈𝑏𝑠2 , 𝜈𝑐𝑠2 ]𝑇
[𝐼𝑠2 ] = [𝑖𝑎𝑠2 , 𝑖𝑏𝑠2 , 𝑖𝑐𝑠2 ]𝑇 (II.7)
[𝛷𝑠2 ] = [𝜙𝑎𝑠2 , 𝜙𝑏𝑠2 , 𝜙𝑐𝑠2 ]𝑇
Les vecteurs de tensions, courant et flux rétorque sont donnés par
[𝑉𝑟 ] = [ 𝜈𝑎𝑟 𝜈𝑏𝑟 𝜈𝑐𝑟 ]𝑇
[𝐼𝑟 ] = [𝑖𝑎𝑟 𝑖𝑏𝑟 𝑖𝑐𝑟 ]𝑇 (II.8)
[𝛷𝑟 ] = [𝜙𝑎𝑟 𝜙𝑏𝑟 𝜙𝑐𝑟 ]𝑇
18
𝐿𝑠1 + 𝐿𝑚𝑠 − 𝐿𝑚𝑠 ⁄2 − 𝐿𝑚𝑠 ⁄2
[𝐿𝑠1,𝑠1 ] = [ − 𝐿𝑚𝑠 ⁄2 𝐿𝑠1 + 𝐿𝑚𝑠 − 𝐿𝑚𝑠 ⁄2 ](II.12)
− 𝐿𝑚𝑠 ⁄2 − 𝐿𝑚𝑠 ⁄2 𝐿𝑠1 + 𝐿𝑚𝑠
2𝜋 4𝜋
𝑐𝑜𝑠 𝛼 𝑐𝑜𝑠 (𝛼 + ) 𝑐𝑜𝑠 (𝛼 + )
3 3
4𝜋 2𝜋
[𝑀𝑠1,𝑠2 ] = 𝐿𝑚𝑠 𝑐𝑜𝑠 (𝛼 + ) 𝑐𝑜𝑠 𝛼 𝑐𝑜𝑠 (𝛼 + ) (II.15)
3 3
2𝜋 4𝜋
[𝑐𝑜𝑠 (𝛼 + 3
) 𝑐𝑜𝑠 (𝛼 + 3
) 𝑐𝑜𝑠 𝛼 ]
2𝜋 4𝜋
𝑐𝑜𝑠 𝜃1 𝑐𝑜𝑠 (𝜃1 + ) 𝑐𝑜𝑠 (𝜃1 + )
3 3
4𝜋 2𝜋
[𝑀𝑠1,𝑟 ] = 𝑀𝑠𝑟 𝑐𝑜𝑠 (𝜃1 + ) 𝑐𝑜𝑠 𝜃1 𝑐𝑜𝑠 (𝜃1 + ) (II.16)
3 3
2𝜋 4𝜋
[𝑐𝑜𝑠 (𝜃1 + 3
) 𝑐𝑜𝑠 (𝜃1 + 3
) 𝑐𝑜𝑠 𝜃1 ]
2π 4π
cos θ2 cos (θ2 + ) cos (θ2 + )
3 3
4π 2π
[Ms2,r ] = 𝑀𝑠𝑟 cos (θ2 + ) cos θ2 cos (θ2 + ) (II.17)
3 3
2π 4π
[cos (θ2 + 3
) cos (θ2 + 3
) cos θ2 ]
Où :
T T T
[Ms2,s1 ]=[Ms1,s2 ] ;[Mr,s1 ]=[Ms1,r ] ;[Mr,s2 ]=[Ms2,r ] .
2
Avec : 𝐿𝑚𝑠 = 𝐿𝑚𝑟 = 𝐿𝑠𝑟 = 3 𝐿𝑚
19
Tel que :
1
𝑤 = 2 ([𝑖𝑠1 ]𝑡 [𝜑𝑠1 ] + [𝑖𝑠2 ]𝑡 [𝜑𝑠2 ] + [𝑖𝑠3 ]𝑡 [𝜑𝑠3 ]) (II.19)
Les deux transformations sont représentées par les deux équations suivantes
[𝐺𝑑𝑞𝑜 ] = [𝐴][𝐺𝑎𝑏𝑐 ](II.24)
[𝐺𝑎𝑏𝑐 ] = [𝐴−1 ][𝐺𝑑𝑞𝑜 ](II.25)
20
Avec :
[Gabc ] : Est le vecteur assemblé des grandeurs du système triphasé équilibre.
[Gdqo ] : Est le vecteur assemble des grandeurs du système biphasé.
De plus les composantes des flux statorique et rétorque sont exprimées par les
équations suivantes.
Stator :
21
𝜙𝑑𝑠1 = 𝐿𝑠1 𝑖𝑑𝑠1 + 𝐿𝑚 (𝑖𝑑𝑠1 + 𝑖𝑑𝑠2 + 𝑖𝑑𝑟 )
𝜙𝑞𝑠1 = 𝐿𝑠1 𝑖𝑞𝑠1 + 𝐿𝑚 (𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 + 𝑖𝑞𝑟 )
(II.28)
𝜙𝑑𝑠2 = 𝐿𝑠2 𝑖𝑑𝑠2 + 𝐿𝑚 (𝑖𝑑𝑠1 + 𝑖𝑑𝑠2 + 𝑖𝑑𝑟 )
𝜙𝑞𝑠2 = 𝐿𝑠2 𝑖𝑞𝑠2 + 𝐿𝑚 (𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 + 𝑖𝑞𝑟 )
Rotor :
𝜙𝑑𝑟 = 𝐿𝑟 𝑖𝑑𝑟 + 𝐿𝑚 (𝑖𝑑𝑠1 + 𝑖𝑑𝑠2 + 𝑖𝑑𝑟 )
(II.29)
𝜙𝑞𝑟 = 𝐿𝑟 𝑖𝑞𝑟 + 𝐿𝑚 (𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 + 𝑖𝑞𝑟 )
22
Stator :
𝑑
𝑣𝑑𝑠1 = 𝑅𝑠1 𝑖𝑑𝑠1 + 𝜙𝑑𝑠1 − 𝜔𝑠 𝜙𝑞𝑠1
𝑑𝑡
𝑑
𝑣𝑞𝑠1 = 𝑅𝑠1 𝑖𝑞𝑠1 + 𝑑𝑡 𝜙𝑞𝑠1 + 𝜔𝑠 𝜙𝑑𝑠1
𝑑
(II.30)
𝑣𝑑𝑠2 = 𝑅𝑠2 𝑖𝑑𝑠2 + 𝑑𝑡 𝜙𝑑𝑠2 − 𝜔𝑠 𝜙𝑞𝑠2
𝑑
𝑣𝑞𝑠2 = 𝑅𝑠2 𝑖𝑞𝑠2 + 𝑑𝑡 𝜙𝑞𝑠2 + 𝜔𝑠 𝜙𝑑𝑠2
Rotor :
𝑑
0 = 𝑅𝑟 𝑖𝑑𝑟 + 𝑑𝑡 𝜙𝑑𝑟 − (𝜔𝑠− 𝜔𝑟 )𝜙𝑞𝑟
𝑑
(II.31)
0 = 𝑅𝑟 𝑖𝑞𝑟 + 𝑑𝑡 𝜙𝑞𝑟 + (𝜔𝑠− 𝜔𝑟 )𝜙𝑑𝑟
Le choix des variables d’état dépend des objectifs lies soit à la commande soit à
l’estimation.
Pour le modèle complet, la vitesse mécanique est une variable d’état aussi. Dans ce
travail le vecteur d’état est choisi comme suit :
T
X =[ϕds1 , ϕqs1 , ϕds2 , ϕqs2 , ϕdr , ϕqr ] .
Et cela à fin d’avoir les caractéristiques dynamiques et statiques de la MASDE, ainsi
sa commande.
Par ailleurs, le flux magnétisant 𝜙𝑚 total est la somme des deux flux
magnétisants direct 𝜙𝑚𝑑 et quadratique 𝜙𝑚𝑞 , d’où :
𝜙𝑚 = √𝜙 2 𝑚𝑑 + 𝜙 2 𝑚𝑞 (II.32)
Les deux expressions des flux magnétisants en fonction des courants statorique et
rétorques sont exprimées par.
𝜙𝑚𝑑 = 𝐿𝑚 (𝑖𝑑𝑠1 + 𝑖𝑑𝑠2 + 𝑖𝑑𝑟 )
(II.33)
𝜙𝑚𝑞 = 𝐿𝑚 (𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 + 𝑖𝑞𝑟 )
En introduisant les expressions des flux magnétisants (II.29) dans les relations du
flux, on obtient le système d’équation suivant.
Stator :
𝜙𝑑𝑠1 = 𝐿𝑠1 𝑖𝑑𝑠1 + 𝜙𝑚𝑑
𝜙𝑞𝑠1 = 𝐿𝑠1 𝑖𝑞𝑠1 + 𝜙𝑚𝑞
(II.34)
𝜙𝑑𝑠2 = 𝐿𝑠2 𝑖𝑑𝑠2 + 𝜙𝑚𝑑
𝜙𝑞𝑠2 = 𝐿𝑠2 𝑖𝑞𝑠2 + 𝜙𝑚𝑞
Rotor:
23
𝜙𝑑𝑟 = 𝐿𝑟 𝑖𝑞𝑟 + 𝜙𝑚𝑑
(II.35)
𝜙𝑞𝑟 = 𝐿𝑟 𝑖𝑞𝑟 + 𝜙𝑚𝑞
De plus, les expressions des courants peuvent être exprimées de cette façon :
Stator :
𝑖𝑑𝑠1 = (𝜙𝑑𝑠1 − 𝜙𝑚𝑑 )⁄𝐿𝑠1
𝑖𝑞𝑠1 = (𝜙𝑞𝑠1 − 𝜙𝑚𝑞 )⁄𝐿𝑠1
(II.36)
𝑖𝑑𝑠2 = (𝜙𝑑𝑠2 − 𝜙𝑚𝑑 )⁄𝐿𝑠2
𝑖𝑞𝑠2 = (𝜙𝑞𝑠2 − 𝜙𝑚𝑞 )⁄𝐿𝑠2
𝐑𝐨𝐭𝐨𝐫:
𝑖𝑑𝑟 = (𝜙𝑑𝑟 − 𝜙𝑚𝑑 )⁄𝐿𝑟
(II.37)
𝑖𝑞𝑟 = (𝜙𝑞𝑟 − 𝜙𝑚𝑞 )⁄𝐿𝑟
24
1 0 0 0
0 1 0 0
0 0 1 0
B= (II.41)
0 0 0 1
0 0 0 0
[0 0 0 0]
Tel que :
𝐿
𝑇𝑠1 = 𝑅𝑠1 ∶ Constante de temps statorique de la première étoile ;
𝑠1
𝐿
𝑇𝑠2 = 𝑅𝑠2 ∶ Constante de temps statorique de la deuxième étoile ;
𝑠2
𝐿𝑟
𝑇𝑟 =𝑅 : Constante de temps rétorque.
𝑟
En remplaçant les flux donnés par (II.29) et (II.30) dans (II.45), on aura .
𝐶𝑒𝑚 = 𝑝𝐿𝑚 [(𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 )𝑖𝑑𝑟 − (𝑖𝑑𝑠1 + 𝑖𝑑𝑠2 )𝑖𝑞𝑟 ](II.46)
De plus à partir des équations des flux rétorques (𝜑𝑑𝑟 𝑒𝑡𝜑𝑞𝑟 ) exprimées par on tire
les expressions des courants rétorques comme suite :
25
1
𝑖𝑑𝑟 = 𝐿 [𝜙𝑑𝑟 − 𝐿𝑚 (𝑖𝑑𝑠1 + 𝑖𝑑𝑠2 )](II.47)
𝑚 +𝐿𝑟
1
𝑖𝑞𝑟 = 𝐿 [𝜙𝑞𝑟 − 𝐿𝑚 (𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 )](II.48)
𝑚 +𝐿𝑟
26
II.6.1. Modélisation de la vitesse du vent
Les aérogénérateurs existent grâce au vent. Le point d'entrée du système est à l'étude.
Une partie de l'énergie cinétique du vent est récupérée, puis convertie en énergie
mécanique via la turbine, et enfin en énergie électrique via le générateur. Les
propriétés dynamiques du vent sont donc cruciales pour l'étude de l'ensemble du
système de conversion d'énergie, car la puissance éolienne varie avec le cube de la
vitesse du vent dans des conditions idéales. La vitesse du vent est un vecteur
tridimensionnel.
Cependant, la direction du vecteur vitesse du vent considérée dans ce modèle est
limitée aux dimensions horizontales. En conséquence, le modèle comportemental des
ventes peut être grandement simplifié. La vitesse du vent est généralement représentée
à l'aide d'une fonction scalaire qui évolue dans le temps.
𝑉 = 𝑓(𝑡)(II.50)
Cette fonction du vent, peut être décomposée en une composante moyenne, variant
lentement, et des fluctuations
27
II.6.3. Génération du couple éolien
L'aspect général d'une turbine à axe horizontal, ainsi que le profil des pales, sont
visibles sur la Figure (II.4). convient également de noter que la turbine utilisée dans
cette étude dispose d'un moyen de contrôler l'angle de calage (Pitch Control).
28
Fig. II.5 coefficient aérodynamique en fonction de ratio de vitesse 𝛌pour une
éolienne de 1.5 MW.
A partir des relevés réalisés sur cette éolienne de 1.5 MW, l’expression du coefficient
de puissance est interpolée sous la forme suivante.
𝜋(𝜆+0.1)
𝐶𝑝 (𝜆, 𝛽 ) = (0.5 − 0.0167(𝛽 − 2)) 𝑠𝑖𝑛 [18.5−0.3(𝛽−2)] − 0.00184(𝜆 −
3)(𝛽 − 2)(II.54)
Connaissant la vitesse de la turbine, le couple aérodynamique (couple de la
turbine éolienne est donc directement déterminé par.
𝑃𝑇 1
𝐶𝑎𝑒𝑟 = = 2𝛺 𝐶𝑝 𝜌𝑆𝑉 3 (II.55)
𝛺𝑡 𝑡
On écrit donc :
1
𝐶𝑎𝑒𝑟 = 2 𝐶𝑡 𝜌𝑆𝑉 2 (II.57)
29
obtient le modèle définit par la relation suivante :
𝑑𝛺𝑟
𝐽 + 𝐷𝛺𝑟 = 𝐶𝑟 − 𝐶𝑒𝑚 (II.60)
𝑑𝑡
Avec :
𝐽
𝐽 = 𝐺𝑡2 + 𝐽𝑒 (II.61)
𝐷
𝐷 = 𝐺 2𝑡 + 𝐷𝑒 (II.62)
Basée sur les équations présentées précédemment, la figure (II.6) peut définir un
modèle physique de la turbine ayant comme entrées l’angle de calage des pales, la
vitesse du vent et le couple électromagnétique fourni par la génératrice.
II.7 CONCLUSION
30
CHAPITRE III:LA COMMANDE VECTORIELLE DE LA GASDE
III.1. INTRODUCTION
Après avoir étudié et modélisé le système pour l'énergie éolienne, nous allons dans ce dernier
consacrée à basée sur la commande vectorielle de la GASDE l'orientation du flux rétorqué, puis on
va modéliser l'alimentation. Electronique du stator constitue de deux convertisseurs commandés par
MLI l'explication de la stratégie de commande de la turbine associée d'une boucle de régulation de
vitesse basée sur régulateur PI classique par mettant d'optimiser l'extraction maximale de l'énergie du
vent
Le but de commodevectoriel, également appelé contrôle d'orientation du flux, est d'aligner une
ou plusieurs composantes du flux (statoique, rotorque ou entrefer) sur un axe de référence qui tourne
avec la vitesse (Fig. III.1 (a)). Ce principe de commande nouvellement développé, qui assure le
fonctionnement stable des machines asynchrones, est la caractéristique clé de la commande
vectorielle, qui conduit à des performances industrielles élevées des machines asynchrones (machines
à papier, laminoirs, traction électrique et générateurs, etc.) qui peuvent supporter fluctuations de
charge
32
En effet dans une machine à courant continu le courant d’induit 𝐼𝑎 contrôle le couple et le courant
inducteur 𝐼𝑓 contrôle le flux. L’expression du couple électromagnétique de la MCC est donne par :
𝐶𝑒𝑚 = 𝐾𝜙𝑓 𝐼𝑎 = 𝐾 ′ 𝐼𝑓 𝐼𝑎 (III.1)
Avec :
𝜑𝑓 : Flux imposé par le courant d’excitation 𝐼𝑓 ;
𝐼𝑎 : Courant d’induit ;
𝐾, 𝐾 ′ : Constantes.
Pour la machine asynchrone double étoile, la commande vectorielle consiste à réaliser un
découplage des grandeurs du couple électromagnétique et du flux rétorqué. Pour cela, il faut choisir
une loi de commande et un système d’axes assurant le découplage du flux et du couple.
Pour la MASDE, on opte pour le choix de l’orientation du flux rétorqué (III.2), car cela permet
d’obtenir une variation de vitesse où le flux et le couple électromagnétique sont indépendamment
commandés à travers les courants statorique direct et quadratique, respectivement.
Les lois de la commande sont obtenues à partir des équations de la MASDE liées au champ
tournant par orientation du flux rétorque. En considérant comme grandeurs de références le flux
rétorque et le couple et en exprimant que
𝜙 𝑑𝑟 = 𝜙𝑟∗
(III.5)
𝜙𝑞𝑟 = 0
La relation finale du couple électromagnétique devient alors :
33
∗ 𝐿𝑚
𝐶𝑒𝑚 = 𝑝𝐿 𝜙𝑟 (𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 )(III.6)
𝑚 +𝐿𝑟
D’après l’équation (II.6), on constate que le couple électromagnétique résulte de l’interaction d’un
terme de flux et d’un terme de courant. Cette expression rappelle le couple de la machine à courant
continu à excitation séparée.
La commande vectorielle de la MASDE peut être soit directe ou indirecte. Dans ce qui suit on
va décrire d’une façon les deux méthodes de la commande vectorielle.
Les lois de commande sont issues des équations qui régissent le fonctionnement de la machine
double étoile asynchrone en fonction du champ tournant et de l'orientation du rotor. Figure(III.2)
représente le schéma de principe de la commande orientée flux du MASDE électromagnétique, en
utilisant le flux rotorique et la paire comme points de référence.
34
Fig. III.2 Schéma de principe de la commande à flux orienté (FOC).
En exprimant que :
𝜙𝑑𝑟 = 𝜙𝑟∗
𝜙𝑞𝑟 = 0 (III.7)
En remplaçant (II.7) dans les équations des tensions rotoriques (I.40et I.41), on obtient.
𝑅𝑟 𝑖𝑑𝑟 = 0 ⟹ 𝑖𝑑𝑟 = 0(III.8)
∗
∗ 𝜔𝑔𝑙 𝜙𝑟∗
𝑅𝑟 𝑖𝑞𝑟 + 𝜔𝑔𝑙 𝜙𝑟∗ = 0 ⟹ 𝑖𝑞𝑟 = − (III.9)
𝑅𝑟
𝜙𝑟∗ 𝐿𝑚
𝑖𝑑𝑟 = 𝐿 −𝐿 (𝑖𝑑𝑠1 + 𝑖𝑑𝑠2 )(III.10)
𝑚 +𝐿𝑟 𝑚 +𝐿𝑟
𝐿𝑚
𝑖𝑞𝑟 = − 𝐿 (𝑖𝑞𝑠1 + 𝑖𝑞𝑠2 )(III.11)
𝑚 +𝐿𝑟
En introduisant (II.) et (II.11) dans le système d’équations des flux statoriques (II.28), on aura.
𝜙𝑑𝑠1 = 𝜆1 ∙ 𝑖𝑑𝑠1 + 𝐿𝑟 𝛽𝑖𝑑𝑠2 + 𝛽 ∙ 𝜙𝑟∗
𝜙𝑞𝑠1 = 𝜆1 ∙ 𝑖𝑞𝑠1 + 𝐿𝑟 𝛽𝑖𝑞𝑠2
(III.12)
𝜙𝑑𝑠2 = 𝜆2 ∙ 𝑖𝑑𝑠2 + 𝐿𝑟 𝛽𝑖𝑑𝑠1 + 𝛽 ∙ 𝜙𝑟∗
𝜙𝑞𝑠2 = 𝜆2 ∙ 𝑖𝑞𝑠2 + 𝐿𝑟 𝛽𝑖𝑞𝑠1
Avec
𝐿𝑚
𝛽= 𝑒𝑡 𝜆1,2 = 𝐿𝑠1,2 + 𝛽𝐿𝑟
𝐿𝑚 +𝐿𝑟
35
En remplaçant (III.14) et (III.12) dans le système d’équations des tensions statoriques(II.30), on
obtient :
∗ ∗
𝑣𝑑𝑠1 = 𝑅𝑠1 ∙ 𝑖𝑑𝑠1 + 𝐿𝑠1 𝑠𝑖𝑑𝑠2 − 𝜔𝑠∗ (𝐿𝑠1 ∙ 𝑖𝑞𝑠1 + 𝜏𝑟 ∙ 𝜙𝑟∗ ∙ 𝜔𝑔𝑙 )
∗ ∗ ∗
𝑣𝑞𝑠1 = 𝑅𝑠1 ∙ 𝑖𝑞𝑠1 + 𝐿𝑠1 𝑠𝑖𝑞𝑠2 + 𝜔𝑠 (𝐿𝑠1 ∙ 𝑖𝑑𝑠1 + 𝜙𝑟 )
∗ ∗
(III.15)
𝑣𝑑𝑠2 = 𝑅𝑠2 ∙ 𝑖𝑑𝑠2 + 𝐿𝑠2 𝑠𝑖𝑑𝑠1 − 𝜔𝑠∗ (𝐿𝑠2 ∙ 𝑖𝑞𝑠2 + 𝜏𝑟 ∙ 𝜙𝑟∗ ∙ 𝜔𝑔𝑙 )
∗ ∗ ∗
𝑣𝑞𝑠2 = 𝑅𝑠2 ∙ 𝑖𝑞𝑠2 + 𝐿𝑠2 𝑠𝑖𝑞𝑠1 + 𝜔𝑠 (𝐿𝑠2 ∙ 𝑖𝑑𝑠2 + 𝜙𝑟
Avec :
𝐿𝑟 ∗
𝜏𝑟 = 𝑒𝑡 𝜔𝑔𝑙 = 𝜔𝑠∗ − 𝜔𝑟
𝑅𝑟
∗ ∗ ∗ ∗
Le système d’équations électriques (II.14) montre que les tensions ( 𝑣𝑑𝑠1 , 𝑣𝑞𝑠1 , 𝑣𝑑𝑠2 , 𝑣𝑞𝑠2 )
influent au même temps sur les composantes des courants statoriques directes et en quadratures
(𝑖𝑑𝑠1 , 𝑖𝑞𝑠1 , 𝑖𝑑𝑠2 , 𝑖𝑞𝑠2 ), donc sur le flux et sur le couple. Il est alors nécessaire de réaliser un
découplage. Cela en définissant des nouvelles variables (𝑣𝑑𝑠1𝑟 , 𝑣𝑞𝑠1𝑟 , 𝑣𝑑𝑠2𝑟 , 𝑣𝑞𝑠2𝑟 ) n’agissent
respectivement que sur (𝑖𝑑𝑠1 , 𝑖𝑞𝑠1 , 𝑖𝑑𝑠2 , 𝑖𝑞𝑠2 ) tel que.
Afin de compenser l’erreur introduite lors de découplage, les tensions statoriques de références
à flux constant sont exprimées par :
∗
𝑣𝑑𝑠1 = 𝑣𝑑𝑠1𝑟 − 𝑣𝑑𝑠1𝑐
∗
𝑣𝑞𝑠1 = 𝑣𝑞𝑠1𝑟 + 𝑣𝑞𝑠1𝑐
∗ (III.19)
𝑣𝑑𝑠2 = 𝑣𝑑𝑠2𝑟 − 𝑣𝑑𝑠2𝑐
∗
𝑣𝑞𝑠2 = 𝑣𝑞𝑠2𝑟 + 𝑣𝑞𝑠2𝑐
36
Où:
∗
𝑣𝑑𝑠1𝑐 = 𝜔𝑠∗ (𝐿𝑠1 ∙ 𝑖𝑞𝑠1 + 𝜏𝑟 ∙ 𝜙𝑟∗ ∙ 𝜔𝑔𝑙 )
∗ ∗
𝑣𝑞𝑠1𝑐 = 𝜔𝑠 (𝐿𝑠1 ∙ 𝑖𝑑𝑠1 + 𝜙𝑟 )
∗
(III.20
𝑣𝑑𝑠2𝑐 = 𝜔𝑠∗ (𝐿𝑠2 ∙ 𝑖𝑞𝑠2 + 𝜏𝑟 ∙ 𝜙𝑟∗ ∙ 𝜔𝑔𝑙 )
𝑣𝑞𝑠2𝑐 = 𝜔𝑠∗ (𝐿𝑠2 ∙ 𝑖𝑑𝑠2 + 𝜙𝑟∗ )
37
Ce défluxage est schématisé par la figure suivante.
Pour des vitesses de rotor inférieures ou égales à la vitesse nominale de la machine, le flux est
généralement maintenu constant à sa valeur nominale ; pour des vitesses plus élevées, le flux diminue
au fur et à mesure que la vitesse augmente pour limiter la tension sur les roulements de la machine.
Ce diagramme décrit le schéma de réglage de la vitesse MASDE, qui comprend un bloc de défluxage
:
38
Fig. III.6 Allure de performance d’une turbine éolienne.
La partie ascendante de la courbe (zone I) représente la zone MPPT précédant le point PO idéal
; dans cette zone, on s'efforce d'extraire le maximum de puissance de la turbine et de fonctionner au
PO. La moitié inférieure de la courbe (zone II) décrit le fonctionnement avec des contraintes où est
supérieur à et Cp est dégradé, permettant de limiter la puissance générée par la turbine.
La vitesse de rotation de la turbine est contrôlée par action sur le couplage électromagnétique
de la génératrice, selon un schéma de fonctionnement du modèle de turbine. La vitesse du vent est
considérée comme un facteur perturbateur de ce système (voir figure II.7).
39
III.6.1. Maximisation de la puissance éolienne
L'algorithme MPPT classique est contrôlé par paires en ajustant la vitesse de rotation mesurée.
Il s'agit de déterminer le couple de référence pour chaque vitesse de rotation en gardant à l'esprit que
les conditions sont idéales ; les vitesses et couple sont donnés comme suit :
𝑉𝑣𝑒𝑛𝑡 𝜆𝑜𝑝𝑡
𝛺𝑜𝑝𝑡 = (III.22)
𝑅
∗ 2
𝐶𝑒𝑚 = 𝐶𝑡_𝑜𝑝𝑡 = 𝐾𝑜𝑝𝑡 𝛺𝑜𝑝𝑡 (III.23)
Avec :
𝑜𝑝𝑡
1 𝐶𝑝
𝐾𝑜𝑝𝑡 = − 3 𝜌𝜋𝑅5
2 𝜆𝑜𝑝𝑡
Le contrôle du système dépend de la zone de fonctionnement.
Zone I : l’usage de l’algorithme MPPT nous permet de déterminer la puissance mécanique de
référence comme suit .
∗
𝑃𝑚𝑒𝑐 = 𝑃𝑚𝑒𝑐_𝑜𝑝𝑡 = 𝐶𝑒𝑚 Ω𝑜𝑝𝑡 (III.24)
Zone II : afin de limiter les contraintes électriques sur la génératrice la puissance mécanique de
référence est maintenue à sa valeur nominale alors on a :
𝑃𝑚𝑒𝑐 = 𝑃𝑛 (III.25)
Dans le présent travail on s’intéresse à la zone I.
Le schéma du bloc de commande est représenté par la figure (III.8), sachant que :
∗ ∗ ∗ ∗
𝑖𝑑𝑠1 = 𝑖𝑞𝑠1 𝑒𝑡 𝑖𝑑𝑠2 = 𝑖𝑞𝑠2
40
Fig. III.8. Représentation schématique du bloc de commande des courants.
41
Deux onduleurs triphasés symétriques sont utilisés pour alimenter le MASDE, comme le
montre ce schéma.
Pour modéliser les deux onduleurs de la figure précédente on commence par la modélisation de
l’onduleur triphasé :
On note. D’après la figure (III.10) que 𝑣𝑎 , 𝑣𝑏 𝑒𝑡 𝑣𝑐 sont les tensions de phases mesurées par
rapport au neutre N de la charge triphasée symétrique. (Les enroulements statoriques de l’étoile 1 de
la MASDE par exemple).
Chaque interrupteur (transistor + diode) supposé idéal. (𝑑𝑜𝑛𝑐 𝐾𝑖𝑗 𝑖 = 1,2 𝑜𝑢3 𝑗 = 1 𝑜𝑢 2) ; On
peut établir les relations suivantes :
42
2𝑣10 − 2𝑣𝑎 + 𝑣𝑏 + 𝑣𝑐 − 𝑣20 − 𝑣30 = 0 (III.29)
La somme des courants (𝑖𝑎 , 𝑖𝑏 , 𝑖𝑐 ) doit être nulle car la charge est à neutre isolé et symétrique. Même
chose pour les tensions des phases (𝑣𝑎 + 𝑣𝑏 + 𝑣𝑐 = 0). Donc, dans l’équation (I.29), on peut
remplacer 𝑣𝑏 + 𝑣𝑐 par −𝑣𝑎 et on tire :
Les tensions de branches 𝑣𝑗0 peuvent êtres égales à E ou à 0 selon l’état des interrupteurs 𝐾𝑖𝑗 .
𝑣𝑎𝑠1 2 −1 −1 𝐹1
𝐸
𝑣
[ 𝑏𝑠1 ] = [−1 2 −1] [𝐹2 ] (III.31)
3
𝑣𝑐𝑠1 −1 −1 2 𝐹3
𝑣𝑎𝑠1 2 −1 −1 𝐹1
𝐸
[𝑣𝑏𝑠1 ] = [−1 2 −1] [𝐹2 ] (III.32)
3
𝑣𝑐𝑠1 −1 −1 2 𝐹3
𝑣𝑎𝑠2 2 −1 −1 𝐹4
𝐸
[𝑣𝑏𝑠2 ] = [−1 2 −1] [𝐹5 ](III.33)
3
𝑣𝑐𝑠2 −1 −1 2 𝐹6
Pour déterminer les moments où les interrupteurs se ferment et s'ouvrent, utilisez la technique
de modulation (MLI) du grand triangle sinusoïdal d'impulsion (inverseur de tension) ou la commande
de courant par la technique MLI à hystérésis (inverseur de courant). L'utilisation de fonctions de
référence ou de modulation distingue la majorité des stratégies de modulation d'impulsions à grande
échelle (MLI) pour les onduleurs. L'utilisation de fonctions de modulation non sinusoïdales permet
d'améliorer significativement les performances du convertisseur en termes de puissance :
augmentation du contenu fondamental et harmonique dans la tension de sortie. Cependant, ces
méthodes rendent difficile, voire impossible, la mise en œuvre d'une commande en temps réel .
43
Commande de la tension par MLI sinus-triangle
Cette méthode consiste à calculer les intersections d’une tension de référence sinusoïdale et
d’une tension de modulation triangulaire. Les six signaux des références pour les deux onduleurs sont
donnés par.
Avec : 𝑘 = 𝑎, 𝑏 𝑜𝑢 𝑐 , 𝑗 = 1 , 2 𝑜𝑢 3 .
𝑡
𝑉𝑝𝑚 [4 𝑇 − 1] 𝑠𝑖 0 ≤ 𝑡 ≤ 𝑇𝑝 ⁄2
𝑝
𝑉𝑝 (𝑡) = { (III.36)
𝑡
𝑉𝑝𝑚 [−4 +3] 𝑠𝑖 𝑇𝑝 ⁄2 ≤ 𝑡 ≤ 𝑇𝑝
𝑇𝑝
𝑓
𝑚 = 𝑓𝑝: L’indice de modulation ;
𝑉
𝑟 = 𝑉 𝑚 : Coefficient de réglage.
𝑝𝑚
Où :
𝑇𝑝 : Le période de modulation
𝑓𝑝 : La fréquence de modulation ;
44
III.7. RESULTATS DE SIMULATION ET INTERPRETATION
A l'aide du logiciel MATLAB, un programme numérique a été créé pour simuler l'évolution des
caractéristiques du système éolien étudié, avec les paramètres listés en annexe. Ce système est relié
électriquement à un puissant réseau triphasé stable et équilibré.
Au début de la simulation, le système éolien est supposé fonctionner à sa position optimale,
telle que : à une vitesse de vent spécifique, le coefficient de puissance, et l'éolienne extrait la
puissance potentielle maximale du vent. Pendant une durée de 15 minutes, le vent change de vitesse
selon une rampe vers une valeur différente, comme indiqué sur le schéma (III.13).
12.3
Vitesse du vent (m/s)
12.2
12.1
12
11.9
0 2 4 6 8 10 12 14 16 18
Temps (s)
145 0
Couple électromagnétique (KN.m)
Vitesse de la GASDE (rad/s)
135 -1.5
-2
130
-2.5
125 -3
0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18
Temps(s) Temps (s)
45
9.005 0.5241
Coeficient de puissance Cp
Vitesse spécifique lambda
9.0048 0.5241
9.0046 0.5241
9.0044 0.5241
9.0042 0.5241
9.004 0.5241
0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18
Temps(s) Temps (s)
1000 1000
0 0
-500 -500
-1000 -1000
0 2 4 6 8 10 12 14 16 18 6 6.01 6.02 6.03 6.04 6.05 6.06
Temps(s) Temps(s)
1000 1000
ias1(A) ias1(A)
ias2 (A) ias2 (A)
Courant statorique
500 500
Zoom de Courant
0 0
-500 -500
-1000 -1000
0 2 4 6 8 10 12 14 16 18 6 6.01 6.02 6.03 6.04 6.05 6.06
Temps(s) Temps(s)
46
5
x 10
2
1
0
0.5 -2
-4
0
-6
-0.5 -8
dr -10
-1
qr -12
-1.5 -14
0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18
Temps (s) Time(s)
[Link]. 11 .Résultats de simulation de l’aérogénérateur asynchrone double étoile et à flux rotorique orienté, muni
d’un régulateur de vitesse PI classique (Profilde vent constant).
Interprétationdesrésultatsobtenus
La puissance active éolienne suit sa référence optimale, cette allure est conforme aussi à celle
ducoupleélectromagnétiquede la GASDE ;
La vitesse de la GASDE est l’image du vent entraînant l’éolienne, elle suit
convenablement saréférence;
Le déphasage de 180 entre la tension et le courant de phase statorique , tradui une
production de la puissance active uniquement au stator comme illustre les figures
des puissances ;
Les allures du couple électromagnétique de la GASDE et de sa référence, sont
pratiquementidentiques, mais différentes de l’allure du profil de la vitesse du vent à
cause du couple dynamiquedûàl’inertiel ;
L allure des composantes du flux rotorique montre une bonne orientation du flux
garantissante une commande vectorielle decouplee de la GASDE.
III.8. CONCLUSION
Dans ce chapitre à été consacré à la commande par orientation du flux rotorique de la MASDE,
l’objectif était la régulation de la vitesse, Pour cela on applique la méthode vectorielle indirecte. Cette
méthode elle est nécessite seulement la connaissance de la position de flux.
L’avantage de la régulation par les algorithmes de réglage classique tels que les régulateurs PI
est assure de bonnes performances dynamiques et statiques et une meilleure poursuite de la
trajectoire dans le cas du système invariable dans le temps mais face à cette avantage, on a un
inconvénient c’est la sensibilité aux incertitudes paramétrique de la machine et leur variation donc
dans ce cas on est obligatoire de cherche une autre technique de régulation plus efficace.
47
CONCLUSION GENERALE
La machine asynchrone double étoile tient une place très importante dans le monde industriel
d’aujourd’hui grâce à sa robustesse et à sa simplicité de construction. Elle intéresse également le
monde de la recherche de l'industrie car, bien qu'elle soit plus chère que la structure triphasée, elle
présente une fiabilité plus élevée et offre plus d’options pour la réalisation des différentes
commandes en moteur ou bien en géné[Link] sur la base de ces avantages , on choisi les MASDE.
Le travail présenté dans ce mémoire la commande vectorielle indirecte appliqué sur une
machine asynchrone double étoile en fonctionnement génératrice à base d’un système de conversion
de l’énergie éolienne alimentée par deux onduleurs de tension par la conception de régulateur
classique de type PI utilisé pour la régulation de la vitesse.
Nous avons tout d’abord entamé ce travail par un rappel sur aperçu dynamique sur les énergies
renouvelables (et en particulier l’énergie éolienne) en Algérie, nous avons décrit et modélisé
mathématiquement la turbine à vitesse variable et les composants associés, la puissance du vent capté
et la puissance électrique convertie. Nous avons ensuite présenté la modélisation mathématique de la
génératrice asynchrone double étoile (GASDE) dans le repère naturel (a, b, c) et dans le repère
biphasé (d,q) de Park. Comme il est fait d’habitude, la complexité de ce modèle est relativement
réduite en utilisant un certain nombre d’hypothèses simplificatrices et en effectuant un changement
de repère par l’application de la transformation de Park. Après la phase modélisation de la chaîne
éolienne, nous nous sommes intéressés à l’aspect commande. Pour maitriser la difficulté de son
réglage on utilisée la commande vectorielle indirecte par orientation du flux. On a munie à cette
technique de commande vectorielle d’un régulateur de vitesse de type PI qui a permis au système
éolien d’obtenir une bonne performance en termes de poursuite du point de fonctionnement optimal à
maximum de puissance.
Les simulations sont réalisées en sachant que l'approche adoptée pour le captage d’air au
niveau de la turbine est la technique MPPT. L'outil logiciel utilisé durant la réalisation des
simulations est Matlab
L’utilisation d’éoliennes pour produire de l’énergie électrique signifie qu’il n’y a pas de
pollution de l’aire car rien n’est rejeté dans l’aire après la transformation électrique.
48
L’utilisation de commande vectorielle pour le réglage de tension est rapide.
L’utilisation de la méthode indirecte est plus attractive et plus utilisée.
L’utilisation de régulateur PI il donne une bonne poursuite, mais un dépassement et un temps
de réponse relativement grands.
49
BIBLIOGRAPHIE
50