REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE L’ENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

UNIVERSITE DE MOHAMED BOUDIAF – M’SILA

FACULTE DE TECHNOLOGIE DEPARTEMENT DE GENIE ELECTRIQUE

MEMOIRE DE FIN D’ETUDES EN VUE DE L’OBTENTION DU DIPLÔME DE

MASTER EN ELECTROMECANIQUE

SPECIALITE : ELECTROMECANIQUE

THEME

Modélisation et commande d’une machine

asynchrone à double alimentation

Proposé et dirigé par : Présenté par :

. BELOUNIS Ouassila ZIANE Mouloud

SIOUANI Ala eddine

Année Universitaire : 2017 / 2018

N° d’ordre
 REMERCIEMENT

Nous remercions Allah de nous avoir donné la

volonté, la santé et la patience qu'il nous a donnée
durant toutes ces longues années d’étude.

Nous remercions sincèrement notre encadreur

Mme BELOUNIS Ouassila pour avoir d'abord
proposer ce thème et diriger ce modeste travail et
pour son continuel suivi et continu tout au long de la
réalisation et la correction de cette mémoire. Elle
nous a présenté un excellent exemple de la
compétence scientifique, discipline et générosité.

Nous adressons nos plus vifs remerciements aux

membres du jury pour l'honneur qu'ils nous font en
acceptant la charge de juger ce travail.

Enfin, nous ne pourrons terminer ces

remerciements sans remercier notre famille et tous
nos amis pour leurs aides, compréhensions,
encouragements et soutiens qu’elles nous ont
apportés tout le long de nos études.
SOMMAIRE

SOMMAIRE
LISTE DES FIGURES
NOTATION
INTRODUCTION GÉNÉRALE……………………………………………………………..……..1

CHAPITRE I : Modélisation de la machine asynchrone à double

alimentation.
I.1 Introduction………………………………………………………………….……………4

I.2 Modélisation de la machine asynchrone à double alimentation.…………..…………..4

I.2.1 Modèle mathématique de la MADA dans le référentiel générale.…..………………..4

I.2.1.1 Les équations électriques de la MADA………………………………………….4

I.2.1.2 Représentation de la MADA dans le système triphasé et biphasé……………….5

I.2.1.3 Les équations magnétiques sous forme matricielle……………………………….5

I.2.1.4 L'équation mécanique…………………………………………………………….7
I.2.2. Choix de référentiel………………………………………………………………………7

I.2.2.1 Référentiel fixé par rapport au stator…………………………………………......7

I.2.2.2 Référentiel lié au rotor……………………………………………………………7

I.2.2.3 Référentiel lié aux champs tournant………………..………………….………….7

I.2.3 Transformation de PARK………………..………………..………………....……..8

I.2.3.1 Application de la transformation de PARK………………..……………..…...…9

I.3 Représentation du modèle de la MADA sous forme d’état………………..………….11

I.4 Expression du couple……..………………..………………..………………..…………14
I.5 Expressions des puissances active et réactive……..………………..………………….15
I.6 Simulation de la MADA……..…..……………..………………..…………………..... .15

I.6.2 Résultats de simulation fonctionnement moteur (MADA) ………………………….15

I.6.2 Interprétations des résultats………………..……………………..…………….……18

I.6.3 Résultats de simulation fonctionnement générateur (GADA) ……………….…..….18

I.6.4. Interprétations des résultats………………..……………………..………….………21

I.7 Conclusion………..…………………………..………………………..…………………21
SOMMAIRE

Chapitre II : Commande vectorielle des puissances

Active et réactive.
II.1 Introduction………..……………………………...…………………..…..…………....23

II.2 Principe de la commande vectorielle………………………………………………….23

II.3 Etablissement du modèle de la MADA……………………………………………….24

II.3.1 Choix du référentiel pour le modèle de la MADA dans le repère d-q ……...……...25

II.3.2 Relation entre le courant statorique et le courant rotorique………………………...26

II.3.3. Expression des puissances actives et réactives dans le repère synchrone………….26

II.3.4 Expressions des tensions rotoriques en fonction de courants rotoriques……………27

II.3.5 Type des régulateurs utilisés………..………………………..…………………...…29

II.4 Commande indirecte des puissances active et réactive………....……….…………...30

II.4.1 Commande indirecte sans boucles des puissances………..…….………………..…30

II.4.2 Commande en boucle fermée………..………………………..………………….…31

II.4.3 Résultats de simulation………..………………………..………………………...…32

II.4.4 Interprétation des résultats………..………………………..………………………..34

II.5 Onduleur de tension………..………………………..………………………..……..35

II.5.1 Modélisation de l’onduleur à deux niveaux………..………………………..……....35

II.5.1.1 Modèle de l’onduleur à deux niveaux………..………………………..…….…..36

II.5.2 Commande triangulaire-sinusoïdale………..………………………..…….………38

II.5.2.1 Algorithme de commande………..………………………..………………….…39

II.5.2.2 Simulation d’onduleur MLI en générale………..………………………….……39

II.5.3 Résultats de simulation avec onduleur………..………………………..………...…40

II.5.4 Interprétation des résultats………..………………………..………………………..42

II.7 Conclusion………..………………………..………………………..………………..42
SOMMAIRE

Chapitre III : Commande par backstepping des puissances

active et réactive

III.1 Commande par backstepping des puissances actives et réactives……..………….44

III.2 Méthodes de Lyapunov …………………………………………………………….………………44

III.2.1 Première méthode de Lyapunov …………………………………………………...44

III.2.2 Deuxième méthode de Lyapunov …………………………………………………45

III.3 Méthode générale de synthèse récursive par backstepping ……………..……...…45

III.4 Application de la commande par Backstepping à la GADA ……………..………...49

III.4.1 Modele de la génératrice asynchrone à double alimentation …………………...…49

III.4.2 Synthese de la loi de commande ……………………………………………….…51

III.5 Résultats de simulation …………………………………………………………...…53

III.6 Interprétation des résultats ………………………………………………………..…56

III.7 Résultats de simulation avec onduleur …………………………………………………..56

III.8 Interprétation des résultats …………………………………………………………..58

III.9 Etude comparative entre les techniques de commande proposées ………………..59

III.10 Conclusion …………………………………………………………………………..61

CONCLUSION GENERALE ……………………………………………………………...62

BIBLIOGRAPHIE ……………………………………………………………………...…..63

ANNEXES …………………………………………………………………………………...66
LISTE DES FIGURES

CHAPITRE I : MODÉLISATION de la machine DE LA MACHINE ASYNCHRONE

A DOUBLE ALIMENTATION .

Figure. I.1 : Representation de la MADA dans le systéme biphasé et triphasée…..………………………….2

Figure. I.2 : Passage de systéme triphasée au systéme biphasé…..………………………..……………………..6

Figure. I.3 : Le couple électromagnétique et la vitesse mécanique de MADA…..……………………….12

Figure. I.4 : Les composantes du courant statoriques de MADA…..………………………..…………………13

Figure. I.5 : Les composantes du courant rotoriques de MADA…..………………………..………………….13

Figure. I.6 : Les composantes du flux statoriques de MADA…..………………………..……………………..13

Figure. I.7 : Les composantes du flux rotoriques de MADA…..………………………..……………………...14

Figure. I.8 : Les puissances actives et réactives statoriques de MADA…..………………………..…………14

Figure. I.9 : Les courants statoriques et rotoriques triphasée (a,b,c) de MADA…..……………………...14

Figure. I.10 : Les courants statoriques triphasée avec zoom…..………………………..………………………15

Figure. I.11 : Les courants rotoriques triphasée avec zoom…..………………………..……………………….16

Figure. I.12 : Les composantes du courant statoriques de la GADA (d-q…..……………………………...16

Figure. I.13 : Les composantes du courant rotoriques de la GADA (d-q…..………………………..………..16

Figure. I.14 : Les composantes du flux statoriques de la GADA (d-q)………………………………..17

Figure. I.15 : Les composantes du flux rotoriques de la GADA (d-q)……………………...................17

Figure. I.16 : Les puissances actives et réactivse statoriques de la GADA (d-q)……………………..17

Figure. I.17 : Couple électromagnétique et la puissance apparent …………………………………....18

CHAPITRE II: COMMANDE VECTORIELLE DES PUISSANCES ACTIVES ET

REACTIVE

Figure II.1: Analogie entre MCC à excitation séparée et GADA.……………………………………..20

Figure II.2: Référentiel diphasé d q lié au champ tournant statorique…………………………………25

Figure II.3 : Modèle de la MADA pour le contrôle des puissances…………………….......................29

Figure II.4: Schéma de régulation des puissances statoriques Ps et Os ……………………................30

Figure II.5 : Schéma bloc de la commande indirecte en boucle ouverte………………………………31

Figure II.6 : Schéma bloc de la commande indirecte en boucle fer mée…………………………...….31

LISTE DES FIGURES

Figure II.7 : La puissance active et réactive statorique et sontréférence…………….………………...32

Figure II.8 :Le courant rotorique selon l’axe d et q……………………………………………………32

Figure II.9 :Le courant statorique selon l’axe d et q…………………………………………………...33

Figure II.10 : Le flux statorique selon l’axe d et q……………………………………………...…….33

Figure II.11 : Le flux rotorique selon l’axe d et q……………………………………………………...33

Figure II.12 : Courant rotorique triphasés et le couple électromagnétique……………………...…….34

Figure II.13 : Onduleur à deux niveaux utilisant des transistors IGBT………………………………..35

Figure II.14: Modèle équivalent de l’onduleur à deux niveaux………………………………...……...36

Figure II.15: Principe et réponses de la commande MLI sinus-triangle……………………………….39

Figure II.16: la puissance active et réactive et sont référence…………………………………………40

Figure II.17: les courants rotoriques quadrature et directes ………………………………….………..40

Figure II.18 : les courants statoriques quadrature et directes …………………………………………40

Figure II.19 : le flux statorique quadrature et directes ………………………………………….……..41

Figure II.20 : le courant statorique des phases A,B,C…………………………………………………41

Figure II.21 : le courant rotorique des phases A,B,C………………………………………………….41

Figure II.22 : le couple électromagnétique…………………………………………………………….42

CHAPITRE III: COMMANDE PAR BACKSTEPPING DES PUISSANCES ACTIVES

ET REACTIVES :

Figure III.1 : Schéma bloc de la structure de commande par Backstepping…………………………..53

Figure III.2 : Les puissances active et réactive statorique et sont référence……………………53

Figure III.3 : Le courant rotorique selon l’axe d et q…………………………………………54

Figure III.4 : Le flux statorique selon l’axe d et q……………………………………………54

Figure III.5 : Le courant statorique selon l’axe d et q. ………………………………………54

Figure III.6: Courant statorique triphasés avec zoom. ………………………………………55

Figure III.7: Courant rotorique triphasés avec zoom…………………………………………55

Figure III.8: le couple électromagnétique. …………………………………………………...55

Figure III.9 : Les puissances active et réactive statorique et sont référence…………………56

LISTE DES FIGURES

Figure III.10 : Le courant rotorique selon l’axe d et q………………………………………..56

Figure III.11 : Le flux statorique selon l’axe d et q. …………………………………………57

Figure III.12 : Le courant statorique selon l’axe d et q……………………………………….57

Figure III.13 : Courant statorique triphasés avec zoom. ……………………………………..57

Figure III.14 : Courant rotorique triphasés avec zoom. ……………………………………...58

Figure III.15 :Le couple électromagnétique.………………………………………………….58

Figure III.16: La puissance active et réactive statorique pour la (C.V) avec le zoom………..59

Figure III.17: La puissance active et réactive statorique pour la (C.B) avec le zoom………..60
LISTE DES SYMBOLES

MADA : Machine asynchrone à double alimentation.

GADA : Génératrice asynchrone à double alimentation.
M.L.I : Modulation de largeur d’impulsions
fs :Fréquence des grandeurs statorique
fr :Fréquence des grandeurs rotorique
fp :Fréquence de la porteuse
S1, S2, S3 : sont les impulsions de commande de l’onduleur
Vdc La tension de bus continue
a, b, c : Indice correspondents aux trois phases a, b, c
C.V : Commande vectorielle
B,C :Commande par backstepping
s,r : Indice du stator, du rotor
d,q : Indice des composantes orthogonales directs et en quadrature
Rs : Résistance des enroulements statoriques.
Rr : Résistance des enroulements rotoriques
ls : Inductance propre d’une phase statorique
lr : Inductance propre d’une phase rotorique
Tr : Constante de temps rotorique
Ts : Constante de temps statorique
Msr : Inductance mutuelle maximale entre le stator et le rotor
Ms : Inductance mutuelle entre deux phases du stator.
Mr : Inductance mutuelle entre deux phases du rotor
Фs,Фr : Respectivement flux stotorique rotorique
Фsd, Фsq :Les flux statorique directes et en quadrature
Фrd, Фrq :Les flux rotorique directes et en quadrature
isa, isb, isc : Courants triphasés statoriques.
ira, irb, irc : Courants triphasés rotoriques
Vsa,Vsb,Vsc : Tensions simples triphasées statoriques.
Vra, Vrb, Vrc : Tensions rotoriques triphasées
Vsd,Vsq : Tension statorique directes et en quadrature

Vrd,Vrq : Tension rotorique directes et en quadrature

: Angel électrique statorique, rotorique
Wr : Pulsation rotorique.
LISTE DES SYMBOLES

Ws : Pulsation statorique.
Wm : Pulsation mécanique
P : Nombre de paires de pôles.
g : Glissement
Lr : Inductance cyclique d’un phase rotorique
Ls : Inductance cyclique d’un phase statorique
M :L’inductance mutuelle.
[Ap] : Matrice de PARK.
[Ap]-1 : Matrice inverse de PARK.
J : Moment d’inertie.
Cem : Couple électromagnétique.
Cr : Couple résistant.
Ps : La puissance statorique active.
Qs : La puissance statorique réactive.
PI : Proportionnel – Intégrale.
kp : Gain de l’action proportionnel du régulateur PI.
ki : Gain de l’action intégrale du régulateur PI.
INTRODUCTION
GÉNÉRALE
INTRODUCTION GÉNÉRALE

INTRODUCTION GÉNÉRALE

L'électricité est une forme d'énergie. C'est un phénomène énergétique associé à la

mobilité ou au repos de particules chargées positivement ou négativement.
Des phénomènes naturels, tels que la foudre, étaient déjà observés dès l'Antiquité, mais
pendant très longtemps l'électricité a terrifié les hommes qui voyaient en elle une
manifestation de la colère divine ou d'un pouvoir surnaturel.
La production d'électricité permet de mettre à disposition de l'ensemble des
consommateurs un approvisionnement adapté à leurs besoins en énergie électrique, à tout
moment.
La production d'électricité se fait depuis la fin du 19ème siècle à partir de différentes
sources d'énergies primaires [1].

Aujourd'hui, la production peut se faire à partir d’une transformation d’énergie

mécanique (liée au mouvement) en énergie électrique. Les centrales électriques peuvent être
comparées à une dynamo de vélo. Dans ces centrales, l'énergie mécanique est convertie en
énergie électrique mais à plus grande échelle. On peut convertir également d'énergie fossile
(charbon, gaz naturel ou pétrole), d'énergie nucléaire, thermique, hydraulique, solaire, et
d’énergie biomasse, ou encore éolienne en énergie électrique [2].

Tout ça à l’aide des machines électriques, en 1832 Nicolas Constant Pixxi réalise la
première machine électrique à induction, comprenant un aimant qui tournait en face des pôles
d’un électroaimant fixe. Il s’agit d’un générateur à courant alternatif qui peut être redressé
grâce à un commutateur rustique permettant l’inversion de polarité. Cette invention n’aura
cependant pas d’application puisqu’à cette époque seul le courant continu était utilisé de
manière usuelle [3].
Depuis l’ère de l’industrialisation, les chercheurs ont été affronté au "comment
commander les machines électriques à des vitesses variables". Car les entraînements
électriques exigent de plus en plus de hautes performances, une fiabilité accrue, et un coût
réduit.
Au 19éme siècle, ce problème a été résolu avec les moteurs à courant continu, dont on
peut modifier la vitesse de rotation de zéro à la vitesse maximale, et ceci en agissant sur la
tension d’induit en gardant le flux d’induction constant. Cependant, la présence du système

Page 1
INTRODUCTION GÉNÉRALE

balais- collecteur les pénalise : ces machines ne peuvent servir ni dans le domaine des grandes
puissances, ni en milieux corrosifs, s’y ajoute aussi l’entretien que nécessite le collecteur [4].

La machine à double alimentation (MADA) en anglo-saxon "doubly-fed induction

machine " c’est une machine asynchrone triphasée avec un rotor bobiné qui peut être alimenté
par des tensions externes. Le schéma de raccordement le plus typique pour cette machine
consiste à connecter le stator directement au réseau, alors que le rotor est alimenté à travers un
convertisseur de puissance contrôlé. Cette solution est plus attractive pour toutes les
applications où les variations de vitesse sont limitées autour de celle de synchronisme.

L'avantage principal de cette machine, est que le convertisseur de puissance doit traiter
seulement (selon le glissement) une fraction de 20-30 % de toute la puissance de système,
ceci signifie que les pertes dans le convertisseur de puissance peuvent être réduites par rapport
au traitement de toute la puissance par celui-ci; et par conséquent le coût de production. En
plus la MADA en mode génératrice permet de produire de l’énergie électrique à fréquence
constante avec une vitesse mécanique variable. De même qu’elle offre une large gamme
opérationnelle par rapport à la machine asynchrone à cage d’écureuil [5].

L’objectif de notre travail est d’introduire des stratégies de commandes pour améliorer
les performances de la machines asynchrone à double alimentation (MADA) utilisée pour la
production de l’énergie électrique.

Ce travail est structuré comme suit:

 Dans le premier chapitre on présentera le modèle de la machine asynchrone à double

alimentation (MADA), sa mise en équation permettra de simuler sont modèle dans le mode
moteur et générateur, ceci a pour objectif la connaissance du comportement de ce type de
machine dans ses différents régimes de fonctionnement.

 Le deuxième chapitre sera consacré à la commande vectorielle en puissance active et

réactive statorique de la machine asynchrone à double alimentation qui est entraînée par une
vitesse fixe. Le stator de la MADA et alimentée par une source triphasée et le rotor est
connectée à un onduleur triphasé à MLI, commandé par la stratégie triangulo- sinusoïdale.
L’alimentation continue de l’onduleur est supposée constante. Les résultats de simulation par
MATLAB/Simulink seront présentés.

Page 2
INTRODUCTION GÉNÉRALE

 Dans le troisième chapitre nous présenterons la théorie de la commande par Backstteping,

en occurrence, toutes les relations concernant ce mode de réglage, Nous aborderons ensuite la
conception de l’algorithme de commande de la machine asynchrone à double alimentation
(MADA) avec ses différentes étapes. Les performances de la méthode seront justifiées par la
simulation.

Nous terminons ce travail par une conclusion générale et quelques perspectives de recherche
envisagées.

Page 3
CHAPITRE I :
MODÉLISATION DE LA MACHINE
ASYNCHRONE A DOUBLE ALIMENTATION.
Chapitre I Modélisation de la machine asynchrone à double alimentation

I.1. Introduction

La modélisation de la machine électrique est une étape primordiale de son développement.

Les progrès de l’informatique et du génie des logiciels, permettant de réaliser des
modélisations performantes et d’envisager l’optimisation des machines électriques.

Le modèle classique de la machine asynchrone à double alimentation correspondant au

schéma monophasé équivalent, n’est valable qu’en régime sinusoïdale à fréquence fixe au
stator, pour étudier son comportement à une fréquence variable et avec des signaux fournis
par les convertisseurs statiques qui ne sont pas sinusoïdaux, il faut faire appel à un modèle
plus compliqué; on considère que chaque armature triphasée statorique et rotorique peut être
représentée par une armature diphasée équivalente à l’aide de la transformation de PARK ,
chaque armature diphasée est donc représentée par deux enroulements identiques placés sur
deux axes en quadrature et directe, tout en tenant compte des hypothèses simplificatrices
suivantes:

- La machine est de construction parfaitement symétrique;

- Même nombre des phases entre le stator et le rotor;
- Répartition sinusoïdale des champs magnétiques de chaque bobinage, le long de l’entrefer;
- L’absence de saturation dans le circuit magnétique [6].

I.2. Modélisation de la machine asynchrone à double alimentation

I.2.1 Modèle mathématique de la MADA dans le référentiel générale

La machine asynchrone à double alimentation comporte trois bobines statoriques

(As,Bs,Cs) décalées entres elles par un angle de (2π/3) et aussi trois bobines rotorique
identique de répartition et similaire à celles du stator,(Ar, Br , Cr). La vitesse de glissement de
ce champ par rapport au stator est Ωs = ws/p, ou ws désigne la pulsation de réseau
d'alimentation et p le nombre de paires de pôles [7].

Les équations des tensions des phases statoriques et rotoriques qui décrivent le
fonctionnement de la machine s’écrivent comme suit:

I.2.1.1 Les équations électriques de la MADA

A partir de la loi de Faraday, on a :

d
 V  R.i (I.1)
dt

Page 4
Chapitre I Modélisation de la machine asynchrone à double alimentation

On applique cette relation sur l'enroulement triphasé du stator et de rotor on trouve:

Pour le stator :
 d sa
V sa  Rs .i sa  dt
 d sb

V sb  R s .i sb  (I.2)
 dt
V  R .i  d sc
 sc s sc
dt
Pour le rotor :
 d ra
Vra  Rr .ira  dt
 d rb
Vrb  Rr .irb  (I.3)
 dt
V  R .i  d rc
 rc r rc
dt

I.2.1.2 Représentation de la MADA dans le système triphasé et biphasé

La machine asynchrone à double alimentation doit être représentée dans le système

triphasé et biphasé par la figure (I.1)
d
Bs
Br
ds

dr
q
Ar
qs
qr

As

Cs

Cr

Fig. I.1 : Représentation de la MADA dans le système biphasé et triphasé.

I.2.1.3 Les équations magnétiques sous forme matricielle
Les équations magnétiques sont données par les expressions suivantes :

Page 5
Chapitre I Modélisation de la machine asynchrone à double alimentation

 s   l s  M sr  is 
   M  l r   ir 
. (I.4)
 r   sr
Tel que :
  SA    rA 
 S    SB  ;  r     rA  (I.4.1)
 SC   rC 

 S  : Vecteur Flux totale statoriques

 r  : Vecteur Flux totale rotoriques

 ias  iar 
is   ibd  ; ir   ibr  (I.4.2)
 ics  icr 

is  : Vecteur courant traversant les enroulements statoriques.

ir  : Vecteur courant traversant les enroulements rotoriques.
 ls Ms Ms
l s   M s ls M s  (I.5)
 M s Ms l s 

 lr Mr Mr 
l r   M r lr M r 
 M r Mr l r 

(I.6)

Où l s , lr sont respectivement les inductances propres des phases statoriques et rotoriques.

M s , M r sont les inductances mutuelles entre phases statoriques et rotoriques.

M rs  est la matrice des inductances mutuelles ou matrice de couplage stator-rotor donnée par:

 cos  r cos( r  2 / 3) cos( r  2 / 3)

M rs  t 
 M sr   M sr .cos( r  2 / 3) cos  r cos( r  2 / 3) (I.7)
cos( r  2 / 3) cos( r  2 / 3) cos  r 

Cette dernière matrice est nommée inductance mutuelle entre le rotor et le stator.

Page 6
Chapitre I Modélisation de la machine asynchrone à double alimentation

I.2.1.4 L'équation mécanique

d 1
Ce  C r  j  f    .(Ce  Cr )  f (I.8)
dt j
Ce : Le couple électromagnétique,

Cr : Le couple resistant,
f : Le coefficient de frottement visqueux de la MADA,
 : La vitesse de rotation,
j : L’inertie des parties tournantes.
Les équations différentielles décrivent le fonctionnement de cette machine sont en fonction
de θr (système à résolution différentielles).
L'application de la transformation de Park, s'avère nécessaire, cette dernière appliquée aux
courants, tensions, et aux flux pour obtenir des équations différentielles à des coefficients
constants.
I.2.2. Choix de référentiel
Choix du repère (U-V) ce qui rend la transformation de Park attrayante, est que
l’orientation du repère (U-V) peut être quelconque. Il existe trois choix importants, le repère
(U-V) peut être fixé au stator, au rotor ou au champ tournant.

I.2.2.1 Référentiel fixé par rapport au stator

Dans ce référentiel, les axes (U-V) sont immobiles par rapport au stator ce
référentiel est le mieux adapté pour travailler avec les grandeurs instantanées et dont
l’avantage ne nécessite pas une transformation vers le système réel. L’utilisation de ce
systéme permet d’étudier les régimes de démarrage et de freinage des machines à courant
alternatif [8].

I.2.2.2 Référentiel lié au rotor

Dans ce référentiel, les axes (U-V) sont immobiles par rapport au rotor tournant à une
vitesse donc ), L’utilisation de ce référentiel permet d’étudier les
régimes transitoires dans les machines alternatives synchrones et asynchrones avec une
connexion non symétrique des circuits du rotor.

I.2.2.3 Référentiel lié aux champs tournant

Dans ce référentiel, les axes (U-V) sont immobiles par rapport au champ
électromagnétique créé par les enroulements statoriques, d’où : ce référentiel est généralement

Page 7
Chapitre I Modélisation de la machine asynchrone à double alimentation

utilisé dans le but de pouvoir appliquer une commande de vitesse, de couple, ...etc. Puisque
les grandeurs dans ce référentiel sont de forme continue

En générale ce dernier est utilisé pour réaliser le contrôle vectoriel du fait que les
grandeurs de réglage deviennent continues. Le système d'axe lié au champ tournant à
l'avantage qu'avec une alimentation sinusoïdale les grandeurs variables sont traitées comme
des grandeurs continues.
Ainsi que ce référentiel est le seul qui n'introduit pas de simplification dans la
formulation des équations, il est très intéressant dans les problèmes ou la fréquence
d'alimentation est constante ce qui simplifie considérablement les calculs [9].
En plus, le choix d'un tel référentiel approprié pour la modélisation repose sur la stratégie
de commande appliquée et les grandeurs à commander [10].

I.2.3. Transformation de PARK

La transformation de Park consiste à appliquer aux courants, tensions, et aux flux un
changement de variable faisant intervenir l'angle entre l'axe d'une phase (Va) statorique et le
système d'axe (d-q), elle est présenté par la figure (I.2), et définie par :

 
 cos  cos(  2 / 3) cos(  2 / 3) 
2
AP   .  sin   sin(  2 / 3)  sin(  2 / 3) (I.9)
3 1
 1 1 
 2 2 2 

L'angle θ dans la matrice de Park [AP] prend la valeur (θs) pour les grandeurs statoriques et
(θs- θr) pour les grandeurs rotoriques.

2
Le facteur pour conserver la puissance électrique instantanée.
3
D'autre part:

 1 
 cos   sin  
 2
AP 1  2
. cos(  2 / 3)  sin(  2 / 3)
1 
(I.10)
3 2
 1 
cos(  2 / 3)  sin(  2 / 3) 
 2

Page 8
Chapitre I Modélisation de la machine asynchrone à double alimentation

Fig. I.2 : Passage de système triphasé au système biphasé.

I.2.3.1 Application de la transformation de PARK

En appliquant la transformation de Park aux équations de la machine asynchrone à

double alimentation dans le repère quelconque (équation (I.2), (I.3)).
Alors le modèle de la machine est obtenu en tenant compte des composantes homopolaire et
quadrature, directs sous la forme suivante :

AP 1.Vdqo   R.AP 1.idqo  d AP 1. dqo  (I.11)

dt

AP .Vdqo   R.idqo  d  dqo   d .AP 1. dqo  (I.12)

dt  dt 

0  1 0 
AP . .AP    1 0 0. d 
d 1 
(I.13)
 dt   dt 
0 0 0

Tel que :
θ = θs : pour les grandeurs statoriques
θ = θs- θr : pour les grandeurs rotoriques

En remplaçant la relation (I.12) dans (I.13) on obtient le modèle biphasé équivalent suivant :

Page 9
Chapitre I Modélisation de la machine asynchrone à double alimentation

 d d d
Vd  [Link]  dt  dt  q
 d q d

Vq  [Link]   d (I.14)
 dt dt
V  R.i  d o
 o o
dt

Remarque

La composante homopolaire du système (I.14) est de valeur nulle pour un système

équilibré ; à partir de ce qui précède on tire les équations suivantes :

 d sd
Vsd  Rs .i sd    s . sq
dt
 d sq
Vsq  Rs .i sq    s . sd
 dt
 d rd
(I.15)
Vrd  Rr .ird   ( s   m ). rq
 dt
 d rq
Vrq  Rr .irq   ( s   m ). rd
 dt

L'équation (I.15) représente le modèle de la machine asynchrone à double alimentation dans

le repère (d,q) lié au champ tournant.
En appliquant la transformation de Park pour établir les relations entre les flux d'axes (d.q.o)
et ceux des axes (A.B.C) on obtient :

 
  dqos  AP 
. s 

 
  dqor  AP 
. r 
(I.16)

Au stator :

 
  dqos  AP .l s 
 . is   M sr .ir 

       
  dqos  AP .l s . A p . idqos  AP M sr AP  . idqor

1 1 (I.17)

Après un calcul prolongé des équations (I.16) et (I.17), on obtient le système matriciel
suivant :

Page 10
Chapitre I Modélisation de la machine asynchrone à double alimentation

 sd   l s  M s 0 0 3 / 2( M sr ) 0 0  i sd 
   0 ls  M s 0 0 3 / 2( M sr ) 0  i sq 
 sq    
 so   0 0 l  2M s 0 0 0  i so 
 
(I.18)
. 
 rd  3 / 2( M sr ) 0 0 lr  M r 0 0  i rd 
 rq   0 3 / 2( M sr ) 0 0 lr  M r 0  i rq 
    
 rà   0 0 0 0 0 lr  2M r  i ro 

Ou :
Ls = ls – Ms : Inductance cyclique statorique.
Lr = lr – Mr : Inductance cyclique rotorique.
M = 3/2Msr : Inductance mutuelle cyclique entre l'enroulement du stator et
celui du rotor.
Los = ls+2Ms : Inductance homopolaire statorique.
Lor = lr+2Mr : Inductance homopolaire rotorique.

Si on élimine la composante homopolaire du système (I.18) on obtient :

   L i  Mi
 sd s sd rd
   L i  Mi
 qs s sq rq
 (I.19)
 rd  Lr ird  Misd

 rq  Lr irq  Misq


L , L sont les inductances cycliques statoriques et rotoriques de la machine.

s r
M L’inductance mutuelle.
 ,  ,  et  sont les flux statoriques et rotoriques directs et quadrature du systéme
sd qs rd rq

I.3. Représentation du modèle de la MADA sous forme d’état

Dans ce cas, nous considérons les flux statorique (  sd ,  sq ) et les courants rotorique
(ird,irq) comme variable d'état, et les tensions (Vsd,Vsq,Vrd,Vrq) comme variables de commande,
on a donc :

X   A.X  B.U  Avec : X = ( Фrd , Фrq ; ird, irq)

U = (Vsd, Vrq ; Vrd , Vrq )

D'après l'équation (1.19) nous avons :

Page 11
Chapitre I Modélisation de la machine asynchrone à double alimentation

  sd M
isd   .ird
Ls Ls

  sq M
isq   .irq
 Ls Ls
 (I.20)
i  M
 rd  .i sd
 rd Lr Lr

  rq M
ird   .i sq
 Lr Lr

En remplaçant les équations (I.15) et (I.19) dans l'équation (I.20) et après le calcul et la
simplification nous aurons :

d 1 M
 dt  sd   T . sd  ωs sq  T .i rd  Vsd
 s s

d 1 M
 dt  sq   T . sq  ωs sd  T .i rq  Vsq
 s s
 (I.21)
 d i  δi  ω i  α.  β.(ω  ω )  M .V  1 V
 dt rd rd r rq sd s r sq
σL r L s
sd
σL r
rd


 d i  δi  ω i  α.  β.(ω  ω )  M .V  1 V

 dt
rq rq r rd sq s r sd
σL r L s
sd
σL r
rd

On pose :
 Ls
Ts  R
 s

 Lr
Tr 
 Rr
 M2
  1 
 Lr L s

  M
 Ts Lr Ls

  1 ( 1  M
2

 )
 Tr Ts Lr Ls

 M
   L L
 r s

En remplaçant σ, Ts, Tr,  et  dans le système d'équation (I.21) on obtient système sous
forme matricielle est donné par :

Page 12
Chapitre I Modélisation de la machine asynchrone à double alimentation

 -1 M   1 0 0 0 

sd   T    
0
0 
s
Ts
   s   sd   0 1 0

d sq   -1 M  sq   M 1 
      s 0 .    0 0  (I.22)
dt i   Ts Ts
  rd  
i L r
L s
Lr 
rd
i      .(s  r )   r    M 1 
 rq    irq   0  0 
  .(s  r )   r    Lr Ls Lr 

On remarque que :

 -1 M 
 T s 0 
Ts
 s

 -1 M 
A     s Ts
0
Ts  (I.23)
 
    .(s   r )  r 
 
  .(s   r )   r  

 1 0 0 0 
 0 1 0 0 
 
 
B   M 0
1
0  (I.24)
 L r Ls L r 
 M 1 
 
L r 
0 0
 L r Ls

Pour faciliter la simulation, nous avons écris la matrice [A] comme suit :

[A] = [A1 ] + ws .[A2 ] + wr .[A3 ] (I.25)

Telle que :
 -1 M 
0 0
 Ts Ts 
 
 -1 M 
A 1    0
Ts
0
Ts  (I.26)
 
  0  0 
 
 0  0  

Page 13
Chapitre I Modélisation de la machine asynchrone à double alimentation

 0 1 0 0  0 0 0 0
  0 0 
1 0 0 0   A3   
0 0
 A 2    0 β 0 1
0 β 0 0  
  (I.27)  β 0 0 0
(I.28)
 β 0 0 0 

I.4. Expression du couple

On a: Pa = [Link] + [Link]
En remplacent les tensions par leur valeurs, on obtient :
d sd d sq
Pa  R.i sd2  .i sd   s . sq .i sd  Rs .i sq2  .i sq   s . sd .i sq
dt dt
d d sq
Pa  Rs (i sd2  i sq2 )  ( sd .i sd  .i sq )   s ( sd .i sq   sq .i sd ) (I.29)
dt dt
Pa  Pj  Pr  Pel

 Pj  Rs .(i sd2  i sq2 )

 d sq
 d sd
  Pr  .i sd  .i sq (I.30)
 dt dt
 Pel   s .( sd .i sq   sq .i sd )

A partir de la relation qui donne l'énergie électromagnétique en fonction du couple:

Pel = ω[Link] (I.31)

On distingue l'expression du couple :

Ce = P. (Ф[Link]-Ф[Link]) (I.32)

En plus, on peut écrire l'expression du couple en fonction des courants statorique et les flux
rotorique comme suit :
M
Ce  P. ( sq .ird   sd .irq ) (I.33)
Lr
Donc l'équation mécanique devient :


1 M 
  . P. ( sq .ird   sd .irq )  Cr  (I.34)
j Lr 
Avec P le nombre de paires de poles de la MADA.

Page 14
Chapitre I Modélisation de la machine asynchrone à double alimentation

I.5. Expressions des puissances active et réactive

Les expressions de puissances active et réactive sont données par:

Ps  (vsd .isd  vsq .isq )

Qs  (vsq .isd  vsd .isq ) (I.35)

I.6. Simulation de la MADA

La mise sous forme d’état du modèle de la machine asynchrone à double alimentation

permet la simulation de la machine avec ces paramètres qui sont donnés en annexe.
L’objectif de l’étude réalisée dans ce chapitre est d’établir un schéma fonctionnel à partir
des tensions d’alimentation nous permettent de déterminer les grandeurs électrique,
électromagnétique et mécanique en fonction du temps en régime dynamique.
A l’aide de logiciel MATLAB, la simulation de la machine asynchrone à double
alimentation pour les deux types de fonctionnement moteur et générateur à été faite. Le
modèle est basé sur les équations obtenues avec la transformation de Park dans le plan (d-q)
lié aux champs tournant. Les résultats montrés sur les figures ci-dessous sont ceux obtenus
par le modèle d’une machine, 4KW pour le fonctionnement moteur et générateur.

I.6.1. Résultats de simulation fonctionnement moteur (MADA)

Les résultats de simulation présentée dans les figures ci-dessous sont ceux obtenus pour le
modèle d’une machine asynchrone à double alimentation fonctionnement moteur "MADA",
alimentée par une source de tension triphasée au niveau du stator avec une fréquence du
réseau qui est 50 Hz et d’amplitude de 220V.
La simulation a été réalisée pour un démarrage à vide suivi par une variation de charge à
l’instant à t = 3s en appliquant un couple résistant Cr = 5 N.m.
200 400
Couple électromagnétique Cem [Nm]

Vitesse mécanique w [Rad/s]

150
300

100
200
50

100
0

-50 0
0 1 2 3 4 5 0 1 2 3 4 5
Temps(s) Temps(s)

Figure I.3 : Le couple électromagnétique et la vitesse mécanique de MADA.

Page 15
 Chapitre I Modélisation de la machine asynchrone à double alimentation

20 20

Courant statorique quadrature isq[A]

Courant statorique direct isd[A]

0 0

-20 -20

-40 -40

-60 -60

-80 -80
0 1 2 3 4 5 0 1 2 3 4 5
Temps(s) Temps(s)

Figure I.4 : Les composantes du courant statoriques de MADA.

60 80
Courant rotorique quadrature irq[A]
Courant rotorique direct ird[A]

40 60

20 40

0 20

-20 0
0 1 2 3 4 5 0 1 2 3 4 5
Temps(s) Temps(s)
Figure I.5 : Les composantes du courant rotoriques de MADA

0.5 0.5
Flux statorique quadrature fsq [wb]
Flux statorique direct fsd [wb]

0
0

-0.5
-0.5
-1

-1
-1.5

-2 -1.5
0 1 2 3 4 5 0 1 2 3 4 5
Temps(s) Temps(s)

Figure I.6 : Les composantes du flux statoriques de MADA.

Page 16
 Chapitre I Modélisation de la machine asynchrone à double alimentation

0.5 1

Flux rotorique quadrature frq [wb]

Flux rotorique direct frd [wb]

0
0.5

-0.5

0
-1

-1.5 -0.5
0 1 2 3 4 5 0 1 2 3 4 5
Temps(s) Temps(s)
4 Figure I.7 : Les composantes du flux rotoriques
4
de MADA.
x 10 x 10
6 6
Puissence réactive statorique [VAR][w]
Puissence active statorique [w]

5 5

4 4

3 3

2 2

1 1

0 0
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)
Figure I.8 : Les puissances actives et réactives statoriques de MADA.
100 100
Isa Ira
Courant statorique Isabc[A]

Courant rotorique Irabc[A]

Isb Irb
50 50
Isc Irc

0 0

-50 -50

-100 -100
0 1 2 3 4 5 0 1 2 3 4 5
Temps(s) Temps(s)
Figure I.9 : Les courants statoriques et rotoriques triphasés (a,b,c) de MADA.

Page 17
 Chapitre I Modélisation de la machine asynchrone à double alimentation

I.6.2. Interprétations des résultats

La modélisation et la simulation de la MADA ont permis de déterminer et d’analyser les

différentes caractéristiques du couple, de la vitesse, le flux et des courants. On remarque que
toutes les grandeurs passent par un régime transitoire avant l’établissement du régime
permanent. Dans le régime transitoire, la machine consomme des courants importants. Après
le régime transitoire les courants diminuent et à l’instant de l’application de la charge ils
augmentent de nouveau.
Au démarrage, le couple électromagnétique présente un pic important et des oscillations.
Lors de l’application d’un couple résistant se traduit par une augmentation du couple
électromagnétique, augmentation de l’amplitude des courants statoriques et une diminution de
la vitesse à cause d’un glissement plus important qui se répercute sur les courants rotoriques
(vérification de la relation (fr =gfs)).

I.6.3 Résultats de simulation fonctionnement générateur (GADA)

Les résultats de simulation présentée sur les figures ci-dessous sont ceux obtenus pour le
modèle d’une machine asynchrone à double alimentation fonctionnement générateur
"GADA", entraînée à une vitesse fixe égale à1440tr/min, alimentée directement par deux
sources de tension triphasée parfait. L’une au niveau du stator avec une fréquence du réseau
qui est 50Hz et d’amplitude de (220 V *√ ;100*pi rad/s;50Hz ) et l’autre au niveau du rotor
avec une amplitude de (12*√ V;4*pi rad/s;2Hz ) et une fréquence égale à la fréquence
rotorique.

100 10
Isa
Courant statorique Isabc [A]

Isb
50 5
Isc

0 0

-50 -5

-100 -10
0 0.2 0.4 0.6 0.8 1 0.16 0.17 0.18 0.19
Temps(s)

Figure I.10 : Les courants statoriques triphasés avec zoom.

Page 18
 Chapitre I Modélisation de la machine asynchrone à double alimentation

100 5
Ira
Courant rotorique Irabc [A]

50 Irb
Irc

0 0

-50

-100 -5
0 0.2 0.4 0.6 0.8 1 0.4 0.6 0.8 1
Temps(s)

Figure I.11 : Les courants rotoriques triphasés avec zoom.

100 100
Courant statorique quadrature Isq [A]
Courant statorique direct Isd [A]

80
50
60

40 0

20
-50
0

-20 -100
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps(s) Temps(s)

Figure I.12 : Les composantes du courant statoriques de la GADA (d-

q).
100 100
Courant rotorique quadrature Irq [A]
Courant rotorique direct Ird [A]

50
50

0
-50

-50 -100
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps(s) Temps(s)

Figure I.13 : Les composantes du courant rotoriques de la GADA (d-q).

Page 19
 Chapitre I Modélisation de la machine asynchrone à double alimentation

0.5 0.5

Flux statorique quadrature fsq [wb]

Flux statorique direct fsd [wb]

0
0

-0.5
-0.5
-1

-1
-1.5

-2 -1.5
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps(s) Temps(s)
Figure I.14 : Les composantes du flux statoriques de la GADA (d-q).
0.5 0.2
Flux rotorique quadrature frq [wb]
Flux rotorique direct frd [wb]

0
0

-0.2
-0.5
-0.4

-1
-0.6

-1.5 -0.8
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps(s) Temps(s)

Figure I.15 : Les composantes du flux rotoriques de la GADA (d-q).

4 4
x 10 x 10
3 4
Puissence réactive statorique [VAR][w]
Puissence active statorique [w]

2 3

1 2

0 1

-1 0

-2 -1
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps (s)
Figure I.16 : Les puissances actives et réactive statoriques de la GADA (d-q).
Page 20
 Chapitre I Modélisation de la machine asynchrone à double alimentation

4
x 10
50 3.5
Couple éléctromgnétique Cem [N.m]

Puissence apparent statorique

0
2.5

-50 2

1.5
-100
1
-150
0.5

-200 0
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps (s)
Temps(s)
Figure I.17 : Couple électromagnétique et la puissance apparente.

I.6.4. Interprétations des résultats :

Les résultats de simulation obtenus montrent l’importance des courants statoriques et

rotoriques pendant le démarrage, après un temps égal à environs 0.02 s, ils se stabilisent et
prennent leurs formes contenu avec une fréquence de 50 Hz pour les courants statoriques. On
remarque que l’augmentation au niveau de la tension rotorique implique d’augmentions au
niveau des grandeurs de la machine, et on peut signaler aussi que la tension et la fréquence
rotoriques influence sur les grandeurs de la génératrice, et influence sur la puissance au niveau
de stator. Le couple électromagnétique est négatif atteint au démarrage une valeur maximal
(−158 Nm).

I.7. Conclusion
L’étude de la modélisation nous a permis d’élaborer un modèle mathématique simple
et adéquat pour commander la machine asynchrone à double alimentation.
Le modèle de la MADA dans le référentiel lié au champ tournant nous a permis
l’analyse de cette machine en régime temporel. Ce choix de référentiel repose sur la
stratégie de commande appliquée et les grandeurs à commander.
Dans ce chapitre, nous avons vu que la machine asynchrone à double alimentation
triphasée a été ramenée à une machine biphasée équivalente à l’aide d’une transformation
de PARK.

Page 21
Chapitre I Modélisation de la machine asynchrone à double alimentation

Le modèle de la "MADA" est un système à équations différentielles dont les coefficients sont
des fonctions périodiques du temps, la transformation de Park nous a permis de le simplifier.

Ce modèle sera exploité pour étudier la stratégie de commande de la machine asynchrone à

double alimentation en puissance active et réactive au chapitre suivant en fonctionnement
générateur, qui est consacré à la commande vectorielle et l’utilisation de ces lois par un
systéme d'orientation de flux statorique.

Page 22
CHAPITRE II :
Commande vectorielle des puissances
actives et réactives
Chapitre II Commande vectorielle des puissances active et réactive

II.1 Introduction

La commande vectorielle par orientation du flux présente une solution attractive pour
réaliser de meilleures performances dans les applications à vitesse variable pour le cas de la
machine asynchrone double alimentée aussi bien en fonctionnement générateur que moteur.
[11]

Il ya deux méthodes pour effectuer la commande en puissance de cette machine, la

Méthode qui consiste à négliger les termes de couplage et à mettre en place un régulateur
indépendant sur chaque axe pour contrôler indépendamment les puissances active et réactive.
Cette méthode sera appelée méthode directe car les régulateurs de puissance contrôlent
directement les tensions rotoriques de la machine. La deuxième méthode consiste à tenir
compte des termes de couplage et à les compenser en effectuant un système comportant deux
boucles permettant de contrôler les puissances et les courants rotoriques. Cette méthode
appelée méthode indirecte. [12]

Le but de ce chapitre est de présenter la technique de la commande vectorielle appliquée

à la machine asynchrone à double alimentation. Quelques résultats de simulation seront
présentés aussi pour montrer l’amélioration des performances dynamiques de la machine.
Ensuite, on calcule le régulateur nécessaire pour la commande, et enfin on termine par les
résultats de la simulation qui signifier la validation de commande.

II.2 Principe de la commande vectorielle

L’expression (I.33) du couple électromagnétique de la machine asynchrone à la forme d’une

différence de couple produit par deux machines à courant continu(en excitation séparée)
tournant en sens inverse l’une par rapport à l’autre. L’orientation du flux revient donc à
annuler l’effet de la machine en inversion de sorte que l’expression du couple
électromagnétique résultant soit de la forme de celle d’une machine à courant continu à
excitation séparée.

Dans une machine à courant continu à excitation séparée, le couple est donné par la relation :
Cem = [Link] .If (II.1)

Où Ia est le courant d’induit, If est le courant d’excitation, K est une constante. En maintenant
le courant If constant par le circuit de réglage de l’inducteur, le couple dépend uniquement du
courant Ia [13].

Page 23
Chapitre II Commande vectorielle des puissances active et réactive

Afin d’obtenir un contrôle analogue à celui de la machine à courant continu à excitation

séparée, est analogue au courant d’excitation, tandis que le courant est analogue au
courant d’induit. Par conséquent, les deux composantes et sont mutuellement
découplées.

La méthode du flux orienté consiste à choisir un système d’axes (d, q) où l’un de ces axes
coïncide avec la direction désirée du flux qui peut être rotorique, statorique ou d’entrefer. Ce
repère nous permet d’avoir deux composantes directe et en quadrature, l’une de flux
statorique et l’autre de courant rotorique. On peut dire que La puissance active dépend
uniquement du courant rotorique Irq , et la puissance réactive dépend uniquement du courant
Ird [14] (voire la figure II.1).

GADA

Figure II.1: Analogie entre MCC à excitation séparée et GADA.

II.3 Etablissement du modèle de la GADA

Nous utilisons la modélisation diphasée de la machine. On oriente le repère (d-q) afin que
l'axe ‘d’ soit aligné sur le flux statorique.
On rappelle d’abord le système équations différentielles de la génératrice décrivant la
machine asynchrone dans un repère lié au champ tournant qui est donné par :

Page 24
Chapitre II Commande vectorielle des puissances active et réactive

 d sd
Vsd  Rs .isd    s . sq
dt
 d sq
Vsq  Rs .isq    s . sd
 dt (II.2)
 d rd
Vrd  Rr .ird   ( s   m ). rq
 dt
 d rq
Vrq  Rr .irq   ( s   m ). rd
 dt

II.3.1. Choix du référentiel pour le modèle de la GADA dans le repère d-q :

En choisissant le référentiel diphasé (d-q) lié au champ tournant et à l’aide de la figure

II.2. en plaçant le vecteur flux statorique  S sur l’axe d, on peut écrire :

  
 sd s
  0 (II.3)
 sq

Vsd = 0
 (II.4)
Vsq = Vs   S  s

Figure II.2: Référentiel diphasé d q lié au champ tournant statorique.

En utilisant les simplifications précédentes, peuvent se simplifier sous la forme suivante :

Vsd  Rs .isd
V  Rs .isq   s . sd
 sq
 d dr
Vrd  Rr .ird   ( s   m ). rd (II.5)
 dt
 d rq
Vrq  Rr .irq   ( s   m ). rd
dt

Page 25
Chapitre II Commande vectorielle des puissances active et réactive

Si on néglige la résistance Rs qui est une hypothèse souvent acceptée pour les machines de
grande puissance utilisée dans la production d’énergie éolienne, les équations des tensions de
la machine se réduisent à la forme suivante [16].

Vsd 0
V   s . s  Vs
 sq
 d rd
Vrd  Rr .ird   ( s   m ). rq (II.6)
 dt
 d rq
Vrq  Rr .irq   ( s   m ). rd
dt

De la même manière que pour les tensions, les équations des flux deviennent :

  L I  Mi
 sd s s sd rd
0  Ls I sq  Mi rq


 rd  Lr I rd  Mi sd
(II.7)

 L I  Mi
 rq r rq sq

II.3.2. Relation entre le courant statorique et le courant rotorique :

A partir des équations des composantes directes et quadrature du flux statorique (II.7),
nous pouvons écrire les équations liant les courants statoriques aux courants rotoriques
comme suit :
ìï F
ïï isd = sd - M .ird
ïï Ls Ls
í
ïï M (II.8)
ïï isq = - .irq
ïî Ls

II.3.3. Expression des puissances actives et réactives dans le repère synchrone

Dans un repère diphasé quelconque les puissances actives et réactives sont exprimées par :


 PS  Vsd isd  Vsq isq
 (II.9)
QS  Vsq isd  Vsd isq


L’application des hypothèses choisies ainsi que le système d’axes nous permet d’établir la
simplification suivante :

Page 26
Chapitre II Commande vectorielle des puissances active et réactive

 PS  Vs isq

QS  Vs isd
Ainsi en remplaçant isd et isq par leurs expressions en fonction des courants rotoriques :

 M
 PS  V s irq
 Ls
 (II.10)
Q  Vs  s  Vs M i
 S Ls Ls
rd

De l’expression (II.6) nous pouvons écrire :

Vs
s  (II.11)
s

D’où d’après (II.10)alors :

 M
P
 S  V s
irq
 Ls
 (II.12)
Q  V s  s  Vs M i
2

 S  s .Ls Ls
rd

Il ressort de l’expression (II.12)que le contrôle des puissances active et réactive au stator est
découplé. En effet, avec une inductance magnétisante constante et un réseau électrique
puissant, la puissance active sera directement proportionnelle au courant rotorique D’axe q, et
2
la puissance réactive proportionnelle au courant rotorique d’axe d à une constante près Vs
 s Ls

imposée par le réseau[14].

II.3.4. Expressions des tensions rotoriques en fonction de courants rotoriques
En remplaçant les courants statoriques direct et en quadrature par leurs expressions dans les
équations (II.7), les expressions des flux rotoriques peuvent être écrites sous la forme
suivante :

  M2 M .Vs
  
 rd  rL    ird 
  Ls   s .Ls (II.13)

   L  M   i
2

 rq  r Ls 
rq
 

Page 27
Chapitre II Commande vectorielle des puissances active et réactive

Avec  r  g s en injectant ces équations dans les expressions des tensions rotoriques :

  M 2  dird  M2 
Vrd  Rr .ird   Lr    g. s  Lr    irq
  Ls  dt  Ls 
 (II.14)
V  R .i   L  M   dirq  g.  L  M   i  g. M .Vs
2 2

 rq r rq  r L  dt s r
Ls 
rd
  s   .Ls

En régime permanent, les termes faisant intervenir les dérivées des courants rotoriques
diphasés disparaissent, nous pouvons donc écrire :

  M2
Vrd  Rr .ird  g. s  Lr 
   irq
  Ls 
 (II.15)
V  R .i  g.  L  M   i  g. M .Vs
2

s r
 rq Ls 
r rq rd
  .Ls

En appliquant la transformation de Laplace aux équations (II.14), on obtient :

   M2   M2 
Vrd   Rr   Lr  .P ird  g. s  Lr    irq
   Ls    Ls 
 (II.16)
   M2   M2 M .Vs
V   R  
 L   .P  i  g .  L    i  g
 Ls    Ls 
rq r r rq s r rq
    .Ls

Avec

Vrd et V rq : Sont les composantes déphasées des tensions rotoriques à imposer à la machine
pour obtenir les courants rotoriques désirés.

( ): Est le terme de couplage entre les deux axes. Une synthèse adéquate des

régulateurs dans la boucle de commande permettra de les compenser.

M .Vs : représente une force électromotrice dépendante de la vitesse de rotation.

g
.Ls

A l’aide des systèmes d’équations (II.12) et (II.16), nous pouvons élaborer le modèle pour
le contrôle des puissances tel qu’il est présenté par le schéma bloc ci-dessous :

Page 28
Chapitre II Commande vectorielle des puissances active et réactive

Figure II.3: Modèle de la MADA pour le contrôle des puissances.

Ce schéma fait apparaître des fonctions de transfert du premier ordre pour les deux axes liant
les tensions rotoriques aux puissances actives et réactives statoriques. Il montre également que
nous pouvons mettre en place une commande vectorielle étant donné qu'à l'influence du
couplage près, chaque axe peut être commandé indépendamment avec chacun son propre
régulateur. Les grandeurs de référence pour ces régulateurs seront : la puissance active pour
l’axe q rotorique et la puissance réactive pour l’axe d rotorique [14].

Donc la composante quadratique de tension du rotor Vqr contrôle la puissance active (le
couple électromagnétique) et la composante directe Vdr contrôle la puissance réactive
échangée entre le stator et le réseau. On retrouve ainsi le découplage si avantageux des
moteurs à courant continu [15].

II.3.5 Type des régulateurs utilisés

Nous allons maintenant procéder à la synthèse de régulateurs nécessaires à la réalisation de la

commande en puissance active et réactive de la MADA. Le régulateur Proportionnel Intégral
(PI) est nécessaire pour assurer que les puissances active et réactive du MADA suivent leurs
consignes, et reste le plus communément utilisé pour la commande de ces machines en
fonctionnement génératrice, ainsi que dans de nombreux systéme de régulation industriel. les
régulateurs de chaque axe ont pour rôle d’annuler l’écart entre les puissances actives et
réactives de référence et la puissance active et réactive mesurée [14].

Page 29
Chapitre II Commande vectorielle des puissances active et réactive

Ce type de régulateur est utilisé pour assurer une meilleure robustesse vis-à-vis des
perturbations internes ou externes.

Les boucles de régulation des puissances active Ps et réactive Qs peuvent se présenter par le
schéma bloc de la figure II.4:

Figure II.4 Schéma de régulation des puissances statoriques Ps et Qs.

II.4 Commande indirecte des puissances active et réactive

Cette méthode consiste à estimer les valeurs des tensions rotoriques à partir des valeurs des
puissances actives et réactives. Ainsi, les courants rotoriques seront régulés indirectement.

II.4.1 Commande indirecte sans boucles des puissances :

Cette commande est basée sur la régulation indirecte des courants rotoriques mesurées qui
seront asservis avec les courants de références déjà déterminés par les puissances statoriques
de références imposés à la machine[16].

Donc on peut dire que les courants rotoriques irq et ird, sont respectivement les images de la
puissance active statorique Ps et la puissance réactive statorique Qs, doivent poursuivre leurs
courants de références [17].

A partir des expressions de la puissance active et réactive statoriques du système (II.12) on

déduit les références des courants rotoriques direct et quadrature suivant les relations.

ìï LS
ïï I rq = - Ps
ïï ref
VS M ref
í (II.17)
ïï LS Vs
ïï I rdref = - F sref + .QS
ïî VS M ws .M

Ces courants seront utilisés comme références à la place des références sur les puissances
active et réactive, on aboutit alors au schéma bloc ci-dessous.

Page 30
Chapitre II Commande vectorielle des puissances active et réactive

Figure II.5 : Schéma bloc de la commande indirecte en boucle ouverte.

Cette configuration consiste à réguler les courants par une boucle dans chaque axe, mais Elle
n’est pas toujours fiable pour les puissances actives et réactives, car aux temps qu’il y a Une
instabilité du réseau va donc provoquer une erreur sur le suivi des consignes de ces
puissances.

II.4.2 Commande en boucle fermée

Ce type de structure (Figure II.6) est consisté pour réguler les puissances de manière optimale,
et il est composé de deux boucles, la première est celle du courant tel qu’effectuée dans la
section précédente, et la deuxième boucle est celle des puissances. Ainsi on aboutit à une
régulation de la puissance à partir de consignes que l’on impose.

Figure II.6 : Schéma bloc de la commande indirecte en boucle fermée.

Page 31
Chapitre II Commande vectorielle des puissances active et réactive

II.4.3 Résultats de simulation :

Les résultats présentés ci-dessous sont obtenus en appliquant les lois de commande en
puissance active et réactive que nous venons de présenter sur le modèle de la machine
asynchrone à double alimentation, sans réglage de vitesse avec un démarrage à vide puis avec
application d’une puissance active P= -3000W entre t=1s et t=3s, et une puissance réactive
Q=0VAR.
Les figures ci-dessous montrent les performances de la commande vectorielle en puissance
actives et réactives statorique appliquée à la" GADA" et entraînée à une vitesse fixe :
1440tr/min.
1500 1500
Psmes Qsmes
1000 Psref Qsref
Puissance réactive statorique [VAR] 1000
500
Puissance active statorique (w)

0
500
-500

-1000 0

-1500
-500
-2000

-2500
-1000
-3000

-3500 -1500
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.7 : La puissance active et réactive statorique et sont référence

20 7

6
Courant rotorique quadrature (A)

Courant rotorique directes (A)

15
5

4
10

5
2

1
0
0
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.8 : Le courant rotorique selon l’axe d et q.

Page 32
Chapitre II Commande vectorielle des puissances active et réactive

20 20

15 15

courant statorique quadrature (A)

courant statorique directes (A)

10 10

5 5

0 0

-5 -5

-10 -10

-15 -15
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.9 : Le courant statorique selon l’axe d et q.

2 1
Flux statorique quadrature (wb)
Flux statorique directes (wb)

1.5 0.5

1 0

0.5 -0.5

0 -1
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.10 : II.12: Le flux statorique selon l’axe d et q.

4 4

3 3
Flux rotorique quadrature (wb)

Flux rotorique directes (wb)

2 2

1 1

0 0

-1 -1

-2 -2
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.11 : II.12: Le flux rotorique selon l’axe d et q.

Page 33
Chapitre II Commande vectorielle des puissances active et réactive
20 20

15 15

couplr électromagnétique (N*m)

10
10
courant rotorique Ir (A)

5
5
0
0
-5
-5
-10
-10
-15
-15
-20

-20 -25
0 1 2 3 4 5
0 1 2 3 4 5
Temps (s)
Temps (s)

Figure II.12 : Courant rotorique triphasés et le couple électromagnétique.

II.4.4. Interprétation des résultats :

Pour évaluer et tester la technique de contrôle indirecte des puissances active et réactive
par des régulateurs PI, une étude de simulation a été effectuée sous l'environnement
MATLAB. Dans cette section, on est amené à représenter toutes les figures de simulation qui
nous permettent d'évaluer les performances apportées par le système de réglage considéré
pour un fonctionnement à vitesse constante. D'après ces résultats on constate que :
La simulation de cette commande a été faite en imposant les puissances active et réactive
statoriques : Pref et Qref. Pendant que la machine est entraînée à vitesse de rotation constante.
Nous remarquons que les grandeurs de référence ont été suivies par la "GADA" pour les
deux puissances active et réactive. On remarque des oscillations aux instants d’application des
échelons de puissance dues au couplage entre les deux axes de la machine et on remarque
aussi que le couple électromagnétique dépend directement de la puissance active et ceci
traduit par sa forme identique de celle de la puissance active Ps. Dans ce cas, on peut conclure
que la puissance active est une conséquence du couple électromagnétique; alors que la
puissance réactive est une conséquence de l'excitation du circuit rotorique.

Page 34
Chapitre II Commande vectorielle des puissances active et réactive

II.5 Onduleur de tension

L’onduleur est un convertisseur statique d’énergie permettant d’onduler la tension du
bus continu et la fournir au bobinage du rotor. La tension de sortie d’un onduleur a une forme
d’onde périodique qui n’est pas sinusoïdale, mais qui peut être très proche de la forme d’onde
souhaitée avec une fréquence désirée.
Cette hypothèse nous permet de faciliter la commande du convertisseur. Ainsi la
bidirectionnalité de la puissance entre le rotor et le réseau est possible à travers un niveau de
tension du bus continu supérieur à la valeur crête de la tension ondulée
Toujours par soucis de simplification, l’onduleur utilisé est à deux niveaux (Figure II.13)
car dans le cadre de cette étude nous avons pour but d’amélioreront la qualité des tensions
ondulées par la diminution des harmoniques [18].
II.5.1 Modélisation de l’onduleur à deux niveaux :

Les progrès réalisés ces dernières années dans la commande des onduleurs MLI à deux
niveaux de tension pilotant en vitesse variable des machines asynchrones de moyenne et
grande puissance ont permis d’améliorer sensiblement ses performances [19].
Le rotor de la GADA est connecté au bus continu par l'intermédiaire d'un onduleur de tension
(Figure II.14). La tension de sortie de ce dernier est contrôlée par une technique de
modulation de largeur d’impulsion (MLI) qui permet le réglage simultané de la fréquence et
de la tension de sortie de l’onduleur .Un onduleur de tension triphasé est constitué de trois
cellules (bras) de commutation de deux interrupteurs. Pour chaque interrupteur est formé par
un transistor (ou un thyristor) monté en tête bêche avec une diode de récupération. Pour
assurer la continuité des courants alternatifs et éviter le court-circuitage de la source [20].

Figure II.13 : Onduleur à deux niveaux utilisant des transistors IGBT.

Page 35
Chapitre II Commande vectorielle des puissances active et réactive

II.5.1.1 Modèle de l’onduleur à deux niveaux:

La symétrie des onduleurs à deux niveaux permet leur modélisation par bras. Après avoir
modéliser chaque couple transistor – diode par un seul interrupteur bidirectionnel « ki » Ki
(i=1, 2, 3, 4, 5,6). Il devient possible de déduire un modèle complet de l’onduleur :

Figure II.14 : Modèle équivalent de l’onduleur à deux niveaux

Pour exprimer les tensions de lignes en fonction de la tension dans l’étage continu et de
l’état des commutateurs, les variables S1 , S2 et S3 se doivent d’être définies en fonction de
l’état des commutations dans les trois branches :

Branche 1 :
S1 = 0 Si k1 est ouvert et k4 est fermé

S1 = 1 Si k1 est fermé et k4 est ouvert.

Branche 2 : (II.18)

S2 = 0 Si k2 est ouvert et k5 est fermé;
S2 = 1 ; Si k2 est fermé et k5 est ouvert.

Branche3 :
S3 = 0 Si k3 est ouvert et k6 est fermé;


S3 = 1 Si k3 est fermé et k6 est ouvert.

Pour simplifier l’étude, on supposera que :

 la commutation des interrupteurs est instantanée.

 la chute de tension aux bornes des interrupteurs est négligeable.
 le système triphasé obtenu à la sortie de l’onduleur est un système triphasé équilibré
en tension ne contenant que les harmoniques impairs différents de trois.

Page 36
Chapitre II Commande vectorielle des puissances active et réactive

Les tensions composées à la sortie du convertisseur s’expriment alors par :

V A  V AO  VON

VB  VBO  VON (II.19)
V  V  V
 C CO ON

Par addition on a : VA  VB  VC  VB  VAO  VBO  VCO  [Link]

1
Donc : V AO  VBO  VCO  [Link]  0  [Link]   (V AO  VBO  VCO ) (II.20)
3

On remplace (II.20) dans (II.19), on aura le système suivant :

Or, si on considère que les tensions sont équilibrées on peut déduire les expressions des
tensions en lignes par rapport aux tensions composées :

 3 1 1
V A  2 V AO  3 VBO  3 VCO

 1 2 1
Vb   V AO  VBO  VCO (II.21)
 3 3 3
 1 1 2
Vc   3 V AO  3 VBO  3 VCO


Ainsi l’onduleur est pris en compte dans les simulations par l’intermédiaire de l’équation
classique alors, on peut écrire le système (II.21) sous la forme matricielle suivante:

V A  2 - 1 - 1  V Ao 
V   .- 1 2 - 1
1 V  (II.22)
 B 3    Bo 
VC  - 1 - 1 2 VCo 

Avec :

 E
 V AO  .S1
2

 E
VBO  .S 2 (II.23)
 2
 E
VCO  .S 3
2

Page 37
Chapitre II Commande vectorielle des puissances active et réactive

S1 , S 2 , S 3 sont les impulsions de commande de l’onduleur

Et à l’aide de l’équation (II.18) on remplace (II.23)dans (II.22) on aura le système suivant :

V A  2 - 1 - 1   S1 
V   .- 1 2 - 1  S 
E
(II.23)
 B 6    2
VC  - 1 - 1 2  S 3 

Le système (II.23) représente le modèle mathématique de l'onduleur triphasé à MLI.

II.5.2 Commande triangulaire-sinusoïdale :

La fonction MLI (Modulation en Largeur d'Impulsion) ou PWM (Pulse Width

Modulation) joue le rôle d'interface entre la partie commande d’un variateur de vitesse et
la machine électrique associée. Cette fonction agit sur l'onduleur de tension (ou de
courant) de la partie puissance du variateur et joue un rôle essentiel avec des
conséquences sur toutes les performances du système [21].

La MLI permet de former pour chaque alternance de la tension de sortie plusieurs créneaux.
Pour obtenir cela, on adopte une certaine technique de commande des interrupteurs de
l’[Link] stratégie triangulaire-sinusoïdale est l’une des techniques les plus simple. Elle
consiste à compare les trois tension de référence, calculées par une transformation inverse de
Parke à un signal (porteuse) triangulaire d’amplitude fixe et de pulsation notamment
supérieure à celle des tensions de références [22]. La stratégie triangulo- sinusoïdale est
caractérisée par deux paramètres qui sont :

 L’indice de modulation m qui représente le rapport de la fréquence de la porteuse à la

fréquence de référence et le coefficient de réglage en tension qui représente le rapport
fp
de l’amplitude de l’onde de référence à la valeur de crête de la porteuse. m 
f

 Les instants de commutation sont obtenus par la comparaison de trois tension de

référence de fréquence f à l’onde triangulaire de haute fréquence fp [20]. r  Vref
UP

L’onde en impulsion est meilleure que l’onde rectangulaire si [23]:

f porteuse  f réference

Page 38
Chapitre II Commande vectorielle des puissances active et réactive

II.5.2.1 Algorithme de commande

Algorithme de commande de la stratégie triangulosinusoidale pour un onduleur à deux

niveaux pour un bras k peut être résumé en 2 étapes :

 Etape 1

Vrefk  U p  Vk  Vdc
 (II.24)
Vrefk  U p  Vk  Vdc

Tel que Vc est la tension du bus continue.

 Etape 2

Vk  Vdc  S k  1
 (II.25)
Vk  Vdc  S k  0

II.5.2.2 Simulation d’onduleur MLI en générale

m=12
500
Tension [V]

Signaux de
0
comparaison
(référence et
-500
0 0.005 0.01 0.015 0.02
Temps(s)
Impulsion de commande

Impulsion de 0
commande du premier
transistor. -1
0 0.005 0.01 0.015 0.02
Tension [V]
500
Tension [V]

Tension composée de
0
sortie de l’onduleur
(E).
-500
0 0.005 0.01 0.015 0.02
Temps(s)

Figure II.15: Principe et réponses de la commande MLI sinus-triangle.

Page 39
Chapitre II Commande vectorielle des puissances active et réactive

II.5.3. Résultats
4
de simulation avec onduleur :
x 10
1 12000
Puissence active statorique [w]

Qsmes
Psmes 10000 Qsref

Psref

Puissence réacive statorique [VAR]

8000
0.5
6000

4000
0
2000

0
-0.5 -2000

-4000

-1 -6000
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.16: La puissance active et réactive et sont référence

30 40
le courant rotorique quadrature [A]
le courant rotorique directes [A]

20
20
10

0 0

-10
-20
-20

-30 -40
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.17 : Les courants rotoriques quadrature et directes

60 30
le courant sttorique quadrature [A]
le courant statorique directes [A]

20
40
10

20 0

-10
0
-20

-20 -30
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)
Figure II.18 : Les courants statoriques quadrature et directes

Page 40
Chapitre II Commande vectorielle des puissances active et réactive

1 2

Flux statorique quadrature (wb)

Flux statorique directes (wb)

0.5 1.5

0 1

-0.5 0.5

-1 0
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure II.19 : Le flux statorique quadrature et directes

40 40
le courant statorique [A]

le courant statorique [A]

20 20

0 0

-20 -20

-40 -40
0 1 2 3 4 5 0 0.2 0.4 0.6 0.8 1
Temps (s) Temps [s]
Figure II.20 : Le courant statorique des phases A,B,C

30

20 6
le courant rotorique [A]

le courant rotorique [A]

4
10
2
0 0
-10 -2
-4
-20
-6
-30
0 1 2 3 4 5 3.4 3.6 3.8 4 4.2 4.4
Temps (s) Temps (s)

Figure II.21 : Le courant rotorique des phases A,B,C

Page 41
Chapitre II Commande vectorielle des puissances active et réactive

50

le couple électromagnétique 0

-50
0 1 2 3 4 5
Temps (s)

Figure II.22 : Le couple électromagnétique

II.5.4. Interprétation des résultats :

Les résultats de simulation montre que toutes les grandeurs et similaires à ceux trouvés
avec une alimentation continus parfaite. Un bon suivi de consigne pour les puissances actives
et réactives statorique est constaté.
Les courants du réseau suivent parfaitement leurs références. La tension à la sortie du
onduleur coté réseau est à deux niveaux ou son amplitude dépend de la valeur de la tension du
bus continu qui est 300 V.

II.6. Conclusion :

Ce chapitre nous a permis d’étudier et d’appliquer la commande vectorielle pour la régulation

de puissance active et réactive statorique de la génératrice asynchrone à double alimentation.
Le choix sur l’orientation du flux a été pris en orientant le flux statorique selon l’axe ‘d’. La
méthode du flux orienté est appliquée depuis quelques années à la MADA reste la méthode la
plus répondue. En effet, celle-ci nous permet non seulement de simplifier le modèle de la
machine mais aussi de découpler la régulation des puissances actives et réactives. A partir de
la simulation numérique, on a constaté qu’effectivement la technique d’orientation du flux
statorique permet de découpler le flux et les puissances de sorte que la composante directe du
courant rotorique contrôle la puissance réactive, et la composante en quadrature contrôle la
puissance active. Ceci nous permet d’obtenir des performances dynamiques élevées similaires
à celle de la MCC.

Page 42
Chapitre II Commande vectorielle des puissances active et réactive

Dans ce chapitre, on a présente les résultats de simulation dans le cas sans onduleur et dans le
cas de l’association de l’onduleur avec la GADA qui est entraînée par une vitesse fixe.
Pour pallier au problème de la sensibilité de la commande aux perturbations et aux
incertitudes paramétriques (telle que la résistance de la charge non linéaire …), nous allons
introduire un deuxième type de commande non linéaire, en l’occurrence, la commande par
backstepping. Cette dernière fera l’objet du chapitre suivant.

Page 43
CHAPITRE III:
COMMANDE PAR BACKSTEPPING DES

PUISSANCES ACTIVES ET REACTIVES

Chapitre III Commande par backstepping des puissances active et réactive

III.1 Introduction

Depuis quelques années, beaucoup de progrès ont été faits dans le domaine de la
commande des systèmes non linéaires. La technique du backstepping a fait son apparition
dans les années 1990 par P. Kokotovic [24] et fait partie de ces nouvelles percées dans ce
domaine. Elle est basée sur la décomposition du système entier de commande, qui est
généralement multi-variable (MIMO) et d’ordre élevé en une cascade de sous-systèmes de
commande du premier ordre [25].

Le terme backstepping signifie que l’on remonte la chaîne d’intégrateurs jusqu’à la

commande, en construisant pas à pas la fonction de Lyapunov et la commande par retour
d’état. Cette méthode de backstepping est donc une méthodologie de synthèse récursive
conjointe d’une fonction de Lyapunov et de la loi de commande associée. Elle met à profit
les relations causales successives pour les construire de manière itérative et systématique [26].

Dans ce chapitre, il sera question d’exposer le principe de la technique du

Backstepping avec application aux systèmes du second et du troisième ordre afin de se
familiariser avec cette technique en vue de l’appliquer à la commande du générateur
asynchrone à double alimentation.

III.2 Méthodes de Lyapunov

La commande des systèmes non linéaire s’appuie sur deux approches possibles. La
première vise à linéariser le système à commander, La deuxième approche consiste à trouver
une Fonction de Commande de Lyapunov garantissant certaines performances pour le
système en boucle fermée.

III.2.1 Première méthode de Lyapunov

Cette méthode permet d’analyser la stabilité, d’un système à partir de l’étude de la stabilité
locale par linéarisation de la dynamique autour d'un point d'équilibre. Cette méthode est d'une
importance limitée, car elle ne permet d'étudier que la stabilité locale et ne donne pas
d’information sur le domaine de stabilité globale. De plus, dû aux approximations du premier
degré (linéarisation), il n'est pas possible de tenir compte de tous les types de phénomènes
non-linéaires. En fait, l’étude locale est surtout intéressante pour justifier ou non la poursuite
de l’étude de la stabilité [27].

Page 44
Chapitre III Commande par backstepping des puissances active et réactive

III.2.2 Deuxième méthode de Lyapunov

Cette méthode est basée sur le concept d'énergie dans un système. Le principe de cette
méthode consiste à analyser la stabilité du système, sans même résoudre les équations
différentielles non linéaires qui le régissent. La stabilité dépend uniquement de l'étude des
variations (signe de la dérivée) de l'énergie, ou d’une fonction qui lui est équivalente, le long
de la trajectoire du système [27].

III.3 Méthode générale de synthèse récursive par backstepping

L’étude de la stabilité d'un système caractérisé par un vecteur d'état x consiste alors à
chercher une fonction V(x) (représentative de l'énergie) de signe défini, dont la dérivée dV/dt
est semi définie et de signe opposé dans le même domaine.

Les systèmes dans ce cas doivent être sous forme triangulaire. La forme générale du
système à analyser est donnée par :

 x  f ( x )  g ( x ) x ,
 1 1 1 0 1 2


 x 2  f 2 ( x1 , x 2 )  g 1( x1 , x2 ) x3 , (III.1)

 x n  f n ( x1 ,........., x n )  g n ( x1 ,............xn )u,


Avec :

x1 x1 , x 2 .........,x n T  R

Le système étant d’ordre n, la mise en œuvre s’effectue en n étapes.

Etape 1: On commence par la première équation du système (III.1) La première référence

désirée est notée :

( x1 ) d a 0  y r (III.2)

Ce qui conduit à l’erreur de régulation suivante :

e1  x1 a 0 (III.3)

Page 45
Chapitre III Commande par backstepping des puissances active et réactive

Ainsi sa dérivée est :

  
e1  x1 a 0 (III.4)


 f1 ( x1 )  g 0 ( x1 ) x2  a 0

Pour un tel système, nous construisons d’abord la fonction de Lyapunov V1 sous une forme
quadratique

1 2
V1  e1 (III.5)
2

Sa dérivée temporelle est :

(III.6)
V 1  e1 e1  e1  f1 ( x1 )  g 0( x1) x2  a 0 
  

 

Assurait la stabilité Pour la dynamique de (III.4) Pour cela, prenons x2  a1 , telle que :


f1 ( x1 )  g 0 ( x1 ) x2  a 0  k1e1 (III.7)

Où k1  0 est une constante de conception.

Ainsi, la loi de commande pour le système (III.4) sera donnée par :

1   
a1   k1e1  a 0  f1 ( x1 ) (III.8)
g 0 ( x1 )  
 

V 1  k1e12  0 (III.9)

Etape 2 : la nouvelle référence désirée sera la variable de commande pour le sous-système

précédent

( x2 ) d  a1 (III.10)

L’erreur de régulation c’est :

Page 46
Chapitre III Commande par backstepping des puissances active et réactive

e2  x2 a1 (III.11)

Sa dérivée est :
  
e 2  x 2 a1 (III.12)


 f 2 ( x1 , x2 )  g 1( x1 , x2 ) x3  a1
La fonction de Lyapunov étendue pour le système (III.12) est :

1
V2  V1  e22 (III.13)
2


1 2 2
2
e1 e2 
Dont la dérivée est :
  
V 2  V 1  e2 e 2 (III.14)


 k1e12  e2  f 2 ( x1 , x2 )  g 1( x1 , x2 ) x3  a 0 


 

Pour la stabilité de système dynamiquement en choisie x3  a 2 , telle que :


f 2 ( x1 , x2 )  g 1( x1 , x2 ) x3  a 0  k1e1 (III.15)

Où k1  0 est une constante de conception.

1   
a2   k 2 e2  a1  f 2 ( x1 , x2 ) (III.16)
g1 ( x1 , x2 )  
 
Avec :

      
g 0 ( x1 )  k 2 e 2  a 0  f ' 2 ( x1 , x2 )   k1e1  a 0  f1 ( x1 ) g '1 ( x1 )
    
a1      (III.17)
2
g 0 ( x1 )

Cette choix implique que :

Page 47
Chapitre III Commande par backstepping des puissances active et réactive


V 2  k1e12  k 2 e22  0 (III.18)

Etape n: De la même façon, pour cette étape la référence à suivre sera :

( xn ) d a n1 (III.19)

On a l’erreur de régulation
en  xn a n1 (III.20)

Sa dérivée est :
  
e n  x n a n1 (III.21)


 f n ( x1 ,......xn )  g n ( x1 ,......xn )u  a n1
La fonction de Lyapunov étendue est :
1
V n  V 1  V2  ........ en2 (III.22)
2


1 2
2

e1  .... en2 
Dont la dérivée est :
  
V n  V 1  ...en e n (III.23)

 k e  ...  en  f n ( x1 ,...xn )  g n ( x1 ,...xn )u  a n1 

 
2

 
1 1

Dans cette dernière étape, on est arrivé à déduire la loi de commande pour le système entier.
Un bon choix doit satisfaire :


f n ( x1 ,......xn )  g n ( x1 ,......xn )u  a n1  k n en (III.24)

Où k n  0 est une constante de conception.

Ainsi, la loi de commande pour le système entier sera donnée par :

  
1 
u  k n en  a n1  f n ( x1 ,...,x2 )
g n ( x1 ....,x2 )   (III.25)
 

Page 48
Chapitre III Commande par backstepping des puissances active et réactive

Ce qui garantit la négativité de la dérivée de la fonction de Lyapunov étendue :


V n  k1e12        k 2 e22  0 (III.26)

III.4 Application de la commande par Backstepping à la GADA

La commande backstepping appliquée à la machine asynchrone à double alimentation

est conçue de telle façon à garder la même structure générale d’une commande vectorielle de
flux, afin de garantir de bonnes performances tout en assurant une régulation et une limitation
des courants.

Cette approche nous permet de déterminer les composantes des tensions de commande
d'alimentation de la machine en assurant la stabilité globale par la théorie de Lyapunov. La
commande hybride ainsi obtenue permet d’assurer le suivi et le contrôle indépendamment de
puissances active et réactive générées, en assurant à chaque fois des dynamiques stables pour
les erreurs entre les grandeurs de références et celles réelles des puissances [16].

III.4.1 Modèle de la génératrice asynchrone à double alimentation

La technique backstepping est basé sur le contrôle vectoriel comme nous avons fait dans le
précédent chapitre.

Dans ce cas, la commande sera appliquée dans le repère (d,q) en considérant l’orientation du
flux statorique suivant :

  
 sd s
  0 (III.27)

 sq

Dans cette partie, on considère le nouveau modèle de la MADA dans le repère ( d,q) donné
par :

Vsd 0
V   s . sq  Vs
 sq
 d rd
Vrd  Rr .ird    r . rq (III.28)
 dt
 d rq
Vrq  Rr .irq    r . rd
dt
 Relations entre puissances statoriques et courants rotoriques

Page 49
Chapitre III Commande par backstepping des puissances active et réactive

 ref M ref
 PS  Vs L i rq
 s
 2
(III.29)
Q ref  VS  Vs M i ref


S
 s .Ls Ls rd

 ref Ls
i rq   V .M .Ps
ref

 s
 (III.30)
i ref   Vs  L s .Qref
 rd  s .M Vs .M
s

 Equations des dérivées de courants rotoriques de références

 ref
Ls
 i rq   .Psref
 Vs .M
 (III.31)
i ref   L s .Qref
 rd Vs .M
s

 Equations des dérivées de courants rotoriques

 d i rd
 dt  V - R .i  g.w .L . .i .
1
L r .
rd t rd s r rq

 (III.32)
 d i rq   V - R .i  g.w .L . .i - g.w . MVs . 1
  rq t rq s r rd s
ws .Ls  L r .
 dt 
III.4.2 Synthèse de la loi de commande

La combinaison du contrôle vectoriel et du backstepping nous permet d’éliminer les

régulateurs PI dans le contrôle vectoriel par un algorithme de contrôle. Le calcul des
grandeurs de références des courants et des tensions de commande se fait sur la base des
fonctions de Lyapunov [16].

Cette commande est basée sur deux étapes qui sont données comme suit :

Cette étape consiste à identifier les erreurs « e1 » « e2 » , qui représentent respectivement

l’erreur entre la puissance active statorique réelle « Ps » et la puissance de référence « Psref »

qui sera directement proportionnelle au courant rotorique d’axe q ,ainsi que la puissance

Page 50
Chapitre III Commande par backstepping des puissances active et réactive

réactive statorique « Qs » et celui de référence Qsref qui sera proportionnelle au courant

rotorique d’axe d.

e1  (i rq  irq )
ref

 (III.33)
e 2  (i rd  ird )
ref

La dérivée de cette erreur est donnée par :

e1  (i ref  i rq )
  

 rq
   (III.34)
e 2  (i ref
rd
 i rd )

La première fonction de Lyapunov est définie par :

1
v  (e12  e22 ) (III.35)
2
Pour avoir la dérivée de la fonction de Lyapunov nulle, il faut choisir les erreurs comme suit :

 

e1   K 1 .E1 et e 2   K 2 .E2

Alors la dérivée de la fonction de Lyapunov par rapport au temps devient :

v  ( K1 .e 2  K 2 .e22 )
1
, K1 0 K2  0 (III.36)

 

A partir des équations (III.31),(III.32) et (III.34) la dérivée des erreurs « e1 » ,« e 2 »

devient:

  L s  ref 1 MVs 
e1   ( .Ps )  .(Vrq - R t .i rq  [Link] .L r . .i rd - [Link] . ) 
  Vs
.M L r
. ws
.L s 
(III.37)

e 2   ( L s .Q  ref )  1 .V - R .i  g.w .L . .i 
  V .M s L r .
rd t rd s r rq 
  s 

   L s ref 1 1 MV 
e1   ( .P s )  .Vrq  .( - R t .i rq  [Link] .L r . .i rd - [Link] . s ) 
  Vs .M L r . L r . ws .L s 
 (III.38)
e   ( L s .Q ref )  1 .V  1 . - R .i  g.w .L . .i 
 2  V .M s L . rd L . t rd s r rq 
  s r r 
Page 51
Chapitre III Commande par backstepping des puissances active et réactive

  L s  ref 1 1 MV 
 K 1 .e1   ( .Ps )  .Vrq  .( - R t .i rq  [Link] .L r . .i rd - [Link] . s )  (III.39)
  Vs .M L r . L r . ws .L s 

 K .e   ( L s .Q ref )  1 .V  1 . - R .i  g.w .L . .i 
 2 2  V .M s L r .
rd
L r .
t rd s r rq 
  s 

  L MV 
Vrq   L r . .( s .P sref  K 1 .e1 )  Vrq - R t .i rq  [Link] .L r . .i rd - [Link] . s  (III.40)
  Vs .M ws .Ls 

 
V   L r . ( s .Q sref  K 2 .e2 )  Vrd - R t .i rd  [Link] .L r . .i rq 
L
 rd
  Vs .M 

RESEAU
[Link]
[Link]

[Link] Backstepping ond GADA

ML
[Link]
control I

[Link] [Link]
ird..mes
+ +
[Link]

Figure III.1 : Schéma bloc de la structure de commande par Backstepping.

III.5 Résultats de simulation :

Les résultats présentés ci-dessous sont obtenus en appliquant les lois de commande en
puissance active et réactive que nous venons de présenter sur le modèle de la machine
asynchrone à double alimentation, avec application d’une puissance active P= -3000W entre
t=1s et t=3s, et une puissance réactive Q=1000VAR entre t=2s et t=4s.

Page 52
 Chapitre III Commande par backstepping des puissances active et réactive

1500 2000
Psmes Osmes
1000 Psref Qsref

500 1500

Puissance réactive statorique [VAR]

Puissanve active statorique (w)

-500 1000

-1000

-1500 500

-2000

-2500 0

-3000

-3500 -500
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure III.2 : Les puissances active et réactive statorique et sont référence

12 10
Idrmes Iqrmes
9
Idrref Iqrref
10
8
le courant rotorique quadrature (A)
le courant rotorique directes (A)

8 7

6
6
5

4 4

3
2
2
0
1

-2 0
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure III.3 : Le courant rotorique selon l’axe d et q.

1 2

0.8 1.8

0.6 1.6
Flux statorique quadrature (wb)
Flux statorique directes (wb)

0.4 1.4

0.2 1.2

0 1

-0.2 0.8

-0.4 0.6

-0.6 0.4

-0.8 0.2

-1 0
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure III.4 : Le flux statorique selon l’axe d et q.

Page 53
 Chapitre III Commande par backstepping des puissances active et réactive

8 8

6
6
4

le courant statorique quadrature (A)

le courant statorique directes (A)

4
2
2
0

0 -2

-4
-2
-6
-4
-8
-6
-10

-8 -12
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure III.5 : Le courant statorique selon l’axe d et q.

10

8 1
6
le courant statorique (A)

4
le coutant statorique (A)

0.5
2

0 0
-2

-4 -0.5
-6

-8 -1
-10
0 1 2 3 4 5 3.8 3.9 4 4.1 4.2 4.3
Temps (s) Temps (s)

Figure III.6: Courant statorique triphasés avec zoom.

15 3

10 2
le courant rotorique (A)

1
le coutant rotorique (A)

0
0

-1
-5

-2
-10

-3
-15
0 1 2 3 4 5 2.4 2.6 2.8 3 3.2 3.4 3.6
Temps (s) Temps (s)

Figure III.7: Courant rotorique triphasés avec zoom.

Page 54
Chapitre III Commande par backstepping des puissances active et réactive

15

10

le couple éléctromagnétique [N*m]

5

-5

-10

-15

-20

-25
0 1 2 3 4 5
Temps (s)

Figure III.8: le couple électromagnétique.

III.6 Interprétation des résultats :

Les figures ci-dessus représentent les différentes courbes obtenues par la commande par
backstepping des puissances active et réactive de la GADA.
La figure III. 2 montre la fiabilité et la robustesse de cette commande qui est donnée un
bon poursuivre de puissance active et réactive statorique sa référence. En commande par
backstepping les réponses sont sans dépassement, rapide en régime transitoire et l’erreur
statique tend vers zéro. Dans notre cas, le réglage des puissances active et réactive par cette
commande apporte une amélioration remarquable par rapport aux régulateurs classiques PI.
Cette technique a permis d’obtenir un découplage parfait entres les deux composantes de la
puissance statorique. La figure III.3 représente les deux composants du courant rotorique,
dont on peut déduire le découplage entre la commande des deux puissances (la puissance
active est l’image du courant Iqr et la puissance réactrice du courant Idr).

III.7 Résultats de simulation avec onduleur :

10000 12000
Psmes Qsmes
8000 Psref 10000 Qsref
Puissance réactive statorique [VAR]

6000
Puissance active statorique (w)

8000

4000 6000

2000
4000

0
2000

-2000
0
-4000
-2000
-6000
-4000
-8000
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure III.9 : Les puissances active et réactive statorique et sont référence

Page 55
Chapitre III Commande par backstepping des puissances active et réactive

20 30
Idrmes Iqrmes
15 Idrref Iqrref
20
10

lecourant rotorique quadrature (A)

lecourant rotorique directes (A)

10
5

0 0

-5 -10

-10
-20
-15
-30
-20

-25 -40
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure III.10 : Le courant rotorique selon l’axe d et q.

1 2

0.8 1.8

0.6 1.6
Flux statorique quadrature (wb)
Flux statorique directes (wb)

0.4 1.4

0.2 1.2

0 1

-0.2 0.8

-0.4 0.6

-0.6 0.4

-0.8 0.2

-1 0
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)

Figure III.11 : Le flux statorique selon l’axe d et q.

35 40

30
30
25
le courant statorique quadrature (A)
le courant statorique directes (A)

20 20

15
10
10
0
5

0 -10

-5
-20
-10

-15 -30
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s)

Figure III.12 : Le courant statorique selon l’axe d et q.

Page 56
Chapitre III Commande par backstepping des puissances active et réactive
10

8
1
6

le courant statorique (A)

4
le coutant statorique (A)

0.5
2

0
0
-2

-4
-0.5
-6

-8
-1
-10
0 1 2 3 4 5 3.8 3.9 4 4.1 4.2 4.3
Temps (s) Temps (s)

Figure III.13 : Courant statorique triphasés avec zoom.

15 3

10 2
le courant rotorique (A)
le coutant rotorique (A)

5 1

0 0

-1
-5

-2
-10

-3
-15
0 1 2 3 4 5 2.4 2.6 2.8 3 3.2 3.4 3.6
Temps (s) Temps (s)
Figure III.14 : Courant rotorique triphasés avec zoom.

15

10
le couple éléctromagnétique [N*m]

-5

-10

-15

-20

-25
0 1 2 3 4 5
Temps (s)

Figure III.15 : Le couple électromagnétique.

Page 57
Chapitre III Commande par backstepping des puissances active et réactive

III.8 Interprétation des résultats

Les figures ci-dessus représentent les différentes courbes obtenues par la commande
backstepping des puissances active et réactive de la GADA qui est entraînée par une vitesse
fixe. Lorsque l'onduleur rotorique est pris en compte dans les simulations. Les perturbations
induites par l'onduleur ne permettent pas de distinguer des différences de comportements
entre les régulateurs en termes de temps de réponse. Toutefois, on peut constater une absence
de dépassements importants et une bonne régularité des signaux de puissance dans le cas de la
commande par backstepping. La présence de l'onduleur et du modèle de la GADA ont permis
d'effectuer une simulation dans des conditions proches de celles d'un système réel.

1500 1500
Psmes Psmes
1000 Psref 1000 Psref

500 500
Puissance active statorique (w)

Puissance active statorique (w)

0 0

-500 -500

-1000 -1000

-1500 -1500

-2000 -2000

-2500 -2500

-3000 C.V -3000

-3500 -3500
0 1 2 3 4 5 0 0.5 1 1.5
Temps (s) Temps (s)
1500 500
Qsmes Qsmes
Qsref 400 Qsref
1000
Puissance réactive statorique [VAR]

Puissance réactive statorique [VAR]

300

200
500
100

0 0

-100
-500
-200

-300
-1000
-400
C.V
-1500 -500
0 1 2 3 4 5 0 0.2 0.4 0.6 0.8 1
Temps (s) Temps (s)

Figure III.16: La puissance active et réactive statorique pour la (C.V)

avec le zoom.

Page 58
Chapitre III Commande par backstepping des puissances active et réactive

III.9 Etude comparative entre les techniques de commande proposées

Le but dans cette partie est de réaliser une comparaison entre les différentes commandes que
nous avons présentées dans ce chapitre et le précédent.

Les résultats que nous avons vus par simulation numérique de la génératrice asynchrone à
double alimentation sous Matlab, cette génératrice est pour l’utilisation de produire
l’électricité. nous permettent d’effectuer une comparaison entre la commande vectorielle
(C.V) et la commande par backstepping (C,B) pour la régulation des puissances active et
réactive statorique de la génératrice asynchrone à double alimentation (GADA).

Les figure (III.16) et (III.17) représentes respectivement les résultats de simulation

d’obtenus pour la variation de la puissance active et réactive avec vitesse constante.

1500 1500
Psmes Psmes
1000 Psref 1000 Psref

500 500
Puissanve active statorique (w)

Puissanve active statorique (w)

0 0

-500 -500

-1000 -1000

-1500 -1500

-2000 -2000

-2500 -2500

-3000 C.B -3000

-3500 -3500
0 1 2 3 4 5 0.5 1 1.5
2000 Temps (s) 500 Temps (s)
Osmes Osmes
Qsref 400 Qsref

1500 300
C.B
Puissance réactive statorique [VAR]

Puissance réactive statorique [VAR]

200

1000 100

500 -100

-200

0 -300

-400

-500 -500
0 1 2 3 4 5 0 0.2 0.4 0.6 0.8 1
Temps (s) Temps (s)

Figure III.16: La puissance active et réactive statorique pour la (C.B) avec le zoom.

Page 59
Chapitre III Commande par backstepping des puissances active et réactive

Ces résultats montrent que les performances de poursuite de la puissance active et réactive
sont satisfaisantes pour les deux types de commande, mais les résultats obtenus par la
commande par backstepping obtenus sont meilleures par apport de la commande vectoriel.

Remarque

On remarque que la commande par backstepping avait un temps de réponse très faible que la
commande vectoriel, en plus la commande vectorielle prend des perturbations très
importants au régime transitoire au démarrage (dépassements), par contre de la commande
par backstepping est très faible presque nul.

Finalement on peut conclure que le contrôle des puissances par backstepping présente des
meilleures performances que celles de la commande vectorielle.

III.10 Conclusion

Dans ce chapitre nous avons vu une nouvelle loi de commande par backstepping appliquée
au système de production d’énergie d’une génératrice asynchrone à double alimentation cette
commande basée sur le principe de la commande vectorielle de la machine à induction, et
permet d’obtenir de façons constructives une loi de commande qui assure la stabilité de la
machine asynchrone à double alimentation en boucle fermée via une fonction de Lyapunov
nous sommes intéressé à cette technique de commande, où nous avons présenté les principales
définitions: « la stabilité » et « l’équilibre » des systèmes non linéaires, la méthode d’analyse
de Lyapunov qui est l’élément essentiel de la technique du backstepping; aussi, nous avons
présenté la synthèse et la technique de cette commande par la méthode directe de Lyapunov.
Enfin, nous avons calculé une loi de commande backstepping pour commander
indépendamment les puissances produites par la GADA à flux orienté. Les résultats de
simulations numériques obtenus montrent l’intérêt croissant d’une telle commande dans les
systèmes électriques. A travers les caractéristiques de réponse, on observe des bonnes
performances c'est-à-dire qu’elle augmente la fiabilité de nos systéme.

Nous avons présenté, une étude comparative des résultats obtenus par les deux commandes
(commande vectorielle par orientation de flux statorique et la commande par backstepping).
Ces résultats obtenus dans le cadre de cette application en simulation mettent en valeur les
performances de chaque technique de commande (dépassement, rapidité…).

Page 60
Chapitre III Commande par backstepping des puissances active et réactive

Finalement, cette commande présente un algorithme de contrôle simple qui a l'avantage

d'être facilement implantable dans un calculateur.

Page 61
CONCLUSION
GÉNÉRALE
 CONCLUSION GÉNÉRALE

CONCLUSION GÉNÉRALE

Le but de ce travail était la modélisation et la commande la machine asynchrone à double

alimentation (MADA).

Nous avons développé la modélisation de cette machine en utilisant la transformation de

PARK. Une étude comparative entre les deux types de fonctionnement Moteur et Générateur
a été faite et ceci a pour l’objectif la connaissance du comportement de ce type de machine
dans ses différents régimes de fonctionnement. Les résultats de simulation de ce modèle sur
Matlab nous confirment la possibilité de double alimentation par la stabilité des différentes
grandeurs de la machine étudiée pour un fonctionnement générateur que moteur.

La génératrice asynchrone été destinée pour la production d’énergie électrique dans

plusieurs domaine, il est donc très utile de raisonner en terme de puissance. Pour cela, la
commande vectorielle élaborée dans le deuxième chapitre est en puissance active et réactive
statoriques. Dans ce chapitre, on a présenté l’étude théorique de la commande vectorielle, ou
on a exprimé les grandeurs statoriques de la GADA en fonction des grandeurs rotoriques
ayant pour but le pilotage de la GADA par le rotor. Le bon suivi des consignes pour les deux
puissances statoriques "active et réactive" par les puissances réelles débitées par le stator de la
machine a montré l’efficacité de la commande appliquée à priori la stratégie de régulation en
MLI par un onduleur à deux niveaux alimenté par une source continue supposée parfaite.

Dans le dernier chapitre, on a appliqué la technique de commande par backstepping, Cette

commande élimine l'utilisation des régulateurs en les remplaçant par des fonctions de
Lyapunov qui doivent être vérifiées en tout moment. L’idée de base est l’application de
Lyapunov pour la machine asynchrone à double alimentation. A l’issue de ce travail, nous
exposons une synthèse des résultats et de simulation obtenus afin d'établir les perspectives
futures envisagées. Cependant, la conception de nouvelles méthodes de commande par
backstepping. Elles sont plus simples et plus facile à mettre en œuvre, et présentent des
propriétés de stabilité.

On applique cette commande pour contrôler indépendamment des puissances actives et

réactives générées par la machine asynchrone découplée par orientation du flux. On peut donc
conclure que la technique de commande backstepping utilisée pour la commande du la GADA
a conduit à des bonnes performances, dans beaucoup de cas, on obtient une qualité de réglage
meilleur par rapport à la commande vectorielle.

Page 62
BIBLIOGRAPHIE
BIBLIOGRAPHIE

BIBLIOGRAPHIE

[1] Site Internet, [Link].

[2] GHENISSA Hocine; « Etude comparative de l’alimentation en électricité d’un site isolé
par voie classique et systéme photovoltaïque » Mémoire de Master, Université de
Biskra.2014.

[3] HADJSAÏD Nourredine , SABONNADIERE Jean-Claude ; « Histoire de l’électricité : de

Thales à la consommation du 21e siècle »; Encyclopedia de l’énergie », Document
électronique; Version [Link] 2015.

[4] AMOUCHAS Nabil , BOUTALEB Amine, « Commande vectorielle d’une machine

asynchrone à double alimentation par régulateur PI », Mémoire de Master, Université de
M’Hamed Bougara, Boumerdes, 2016.

[5] Mohamed Adjoudj, Mohamed Abid , Abdelghani Aissaoui, Youcef Ramdani , Houria
Bounoua, « Commande par mode glissant d’une machine asynchrone à double alimentation
montée dans une éolienne », Laboratoire IRECOM Faculté des sciences de l’ingénieur,
Université Djillali Liabes, Algérie,2010.

[6] BOYETTE Armand, « Contrôle - Commande d'un générateur asynchrone à double

alimentation avec systéme de stockage pour la production éolienne », Thèse doctorat,
l’Université de Henri Poincaré, Nancy 1, 2006.

[7] POITIERS Frédéric, « Etude et commande de génératrices asynchrones pour l’utilisation

de l’énergie éolienne - Machine asynchrone à cage autonome - Machine asynchrone à double
alimentation reliée au réseau », Thèse doctorat, l'Ecole polytechnique de l'Université de
Nantes, 2003.

[8] BOMMÉ Edouard, « Modélisation et optimisation des machines électriques discoïdes à

double entrefer », Thèse doctorat de l’institut polytechnique de Grenoble, 2009.

[9] A. Tamaarat, A. Benakcha, A. Menacer, « Commande des puissances active et reactive

dans d’une chaîne de conversion basée sur une machine asynchrone à double alimentation »
Laboratoire LGEB, Université de Biskra, Algérie,2013.

[10] Rezgui Salah eddine, « Commande de la Machine électrique en environment », Mémoire

de Magister, Université de Université Mentouri de Constantine, 2009.

Page 63
BIBLIOGRAPHIE

[11] DRID Saïd, « Contribution à la modélisation et à la commande robuste d’une machine à

induction double alimentée à flux orienté avec optimisation de la structure d’alimentation :
théorie et experimentation», Thèse doctorat, l’Université de Batna, Algérie , 2005.

[12] M. Allam, B. Dehiba, M. Abid, [Link], [Link], « Etude comparative entre la

commande vectorielle directe et indirecte de la machine asynchrone à double alimentation
(MADA) dédiée à une application éolienne », Laboratoire IRECOM (Interaction Réseau
Convertisseur Machines), Université de Sidi Bel Abbès, Algerie, 2014.

[13] BOUKHELIFA Akkila, « Les élements d’optimisation du pilotage d’une machine

asynchrone en vue d’un control vectoriel», Thèse doctorat, l’École National Polytechnique
ENP, Algérie, 2007.

[14] ROUABHI Riyadh, « Étude et commande d’une machine asynchrone à double

alimentationapplication : energie éolienne » , Mémoire de Magister, Université Ferhat Abbas-
Sétif, 2012.

[15] KHETTACHE Laid. « Etude et commande d’un systéme éolienne à base d’une machine
électrique double Alimentée » , Mémoire de Magister, Université de Batna, 2007.

[16] ROUABHI Riyadh, « Contrôle des puissances générées par un système éolien à vitesse
variable basé sur une machine asynchrone double Alimentée », thèse de Doctorat,
L’université de Sétif. 2016.

[17] BEN NADJI Hocine , BELLOUL Abdallah, « Commande de la puissance active et

réactive par la méthode de mode glissant », Mémoire de Master, Université Echahid Hamma
LakhdarEl-Oued, 2015.

[18] BENCHERIF Yacine Abderrahmane, « Modélisation et commande d’une machine

asynchrone à double alimentation pour la production de l’énergie éolienne » mémoire
d'ingénieur d'état en Electrotechnique, l’École National Polytechnique ENP, , 2008

[19] DEPERNET Daniel, « Optimisation de la commande d’un onduleur MLI à trois niveaux
de tension pour machine asynchrone », thèse de Doctorat, L’université de Reims Champagne-
Ardenne, 1995.

[20] MEZERZI Imad, « Contribution à la modélisation et la commande d’une éolienne à base

d’une génératrice asynchrone à double alimentation», Mémoire de Master, Université de
M’sila, 2013.

Page 64
BIBLIOGRAPHIE

[21] CAPITANEANU Stefan Laurentiu, « Optimisation de la fonction MLI d’un onduleur de

tension deux-niveaux », L’université de L’institut national Polytechnique de Toulouse, 2002.

[22] K. Idjdarene, D. Rékioua , T. Rekioua , A. Tounzi, « Contrôle d’une éolienne en

fonctionnement autonome basée sur une génératrice asynchrone », Laboratoire
d’électrotechnique et d’électronique de puissance, Université de Lille, France, 2007.

[23] NOUADJA Mohammed et CHACHA Lokman, « Commande par mode glissant

appliquée à un système électro-énergétique », Mémoire de Master, Université de Kasdi
Merbah ourgla, 2017.

[24] Tarek Ahmed-Ali , « Quelques notions sur la commande non linéaire par modes glissants
et la commande non linéaire de type Backstepping », Ecole national supérieure des ingénieurs
de Caen, 2011.

[25] RAFANOTSIMIVA Liva Falisoa. « Etude de commande non lineaires pour réseau
électriques –APPLICATION A UN SYSTÈME SMIB », Thèse de Doctorat, L’université
d’antsiranana école supérieure Polytechnique, 2013.

[26] BENSETTI Nour_el_islam et BOUKHACHBA Mimouna, « Commande par

backstepping de la machine asynchrone alimentée par un onduleur à trois niveaux »,Mémoire
de Master, Université de Kasdi Merbah ourgla, 2017.

[27] BADRAOUI Mourad, BEREKSI REGUIG Meriem, « Etude et modélisation du moteur

synchrone simple et double étoile avec application de la commande backstepping à la
machine synchrone », Mémoire de Master, Université de Abou BekrBelkaïd – Tlemcen, 2016.

Page 65
ANNEXES
ANNEXES

Les paramétres de la machine à induction double Alimentée

Valeurs nominales :
Puissance nominale : Pn =4 kW.
Tension nominale : V/U= 220/380 – 50 Hz.
Courant nominale : i/I=15/8.6A
Lavitesse nominale : Ωn= 1440 tr/min.
Nombre de paires de pôles : p = 2.

Paramètres de la MADA:
Résistance statorique: : Rs = 1.2 Ω.
Résistance rotorique: : Rr = 1.8 Ω.
Inductance cyclique statorique: : Ls = 0.1554 H.
Inductance cyclique rotorique: : Lr = 0.1568 H.
Inductance cyclique mutuelle: : Lm = 0.15 H.
Le nombre de paires de pôles est: : p = 2.

Constantes mécaniques de la MADA:

Moment d’inertie : J mach =0.2 kg.m².
Coefficient de frottement visqueux : f=0.001 N.m.s/rd.

Page 66
 MEMOIRE DE FIN D'ETUDES EN VUE DU L'OBTENTION DU DIPLÔME
MASTER EN GENIE ELECTRIQUE

Spécialité: Electromécanique
Propose et dirige par: Mme. BELOUNIS Ouassila

Présente par: ZIANE Mouloud & SIOUANI Ala eddine

Thème :

Modélisation et commande d’une machine

asynchrone à double alimentation

Résumé :
Ce travail a pour objectif la modélisation et la commande d’une machine asynchrone à
double alimentation. La modélisation dans le repère de Park de cette machine « GADA » a été
présentée et simuler en utilisant MATLAB, ainsi que sa commande vectorielle, et en fin on a
présenté une stratégie de commande par backstepping appliquée au système .Cette technique
trouve sa plus forte justification au problème d’utilisation d’une loi de commande non linéaire
robuste aux incertitudes du modèle. L’objectif est d’appliquer cette commande pour contrôler
indépendamment les puissances active et réactive générées par la machine asynchrone
découplée par orientation du flux statorique.

Mots clés :
Génératrice asynchrone à double alimentation « GADA », commande vectorielle, Onduleur
triphasé à deux niveaux, contrôle des puissances, réglage par Backstepping.

N° d’ordre :