Commandes non linéaires des machines asynchrones
Commandes non linéaires des machines asynchrones
UNIVERSITE DE 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 GENIE ELECTRIQUE
SPECIALITE: COMMANDE DES SYSTEMES ELECTRIQUES
THEME
N° d’ordre : 86
A l’issue de cette étude, je remercie ALLAH clément et
miséricordieux tout puissant qui m’a donné tant de courage, de
volonté, de patience, de force, de sagesse et d’abnégation pour
mener à terme ce mémoire.
SAHRAOUI Mebarek
~I~
Dédicace
Je dédie ce modeste travail à :
La personne la plus chère dans le monde, la lumière de
ma vie, la source de tendresse, celle qui a sacrifié et a souffert
les plus belles années de sa vie pour me voir un jour réussir :
♥♥ Ma chère mère ♥♥
SAHRAOUI Mebarek
~ II ~
CHAPITRE I
Figure (I-01) Machine asynchrone. 04
Figure (I-02) Modèle d une machine asynchrone triphasée. 05
Figure (I-03) Passage de système triphasé au système biphasé 09
Figure (I-04) Définition des axes réel de la MAS par rapport au référentiel (α β) 12
Figure (I-05) Définition des axes réel de la MAS par rapport au référentiel (x y) 13
Figure (I-06) Définition des axes réels de la MAS par rapport au référentiel(d q) 13
Figure (I-07) Représentation d’un estimateur de flux 14
Figure (I-08) Résultat de simulation du modele de la MAS à vide 18
Figure (I-09) Résultat de simulation du modele de la MAS en charge 17
CHAPITRE II
Figure (II-01) Schéma de principe de la linéarisation entrée sortie 24
Figure (II-02) Schéma bloc du système linéarise 26
Figure (II-03) Dynamique du système linéarise 27
Figure (II-04) Test de robustesse avec variation de la charge 31
Figure (II-05) Test de robustesse avec variation de la vitesse 32
Figure (II-06) Test de robustesse avec variation de la résistance rotorique 33
CHAPITRE III
Figure (III-01) Théorie de la commande par mode de glissement 36
Figure (III-02) Existence du mode de glissement 38
Figure (III-03) Fonction « SAT » 40
Figure (III-04) Fonction « SMOOTH » 41
Figure (III-05) Représentation de la fonction « Sign » 42
Figure (III-06) schéma block de simulation 49
Figure (III-07) Test de robustesse avec variation de la charge 50
Figure (III-08) Test de robustesse avec variation de la vitesse 51
Figure (III-09) Test de robustesse avec variation de la résistance rotorique 52
Figure (III-10) Les surfaces de glissement, du flux et de la vitesse 53
Table des notations et symboles
obs
A ,B ,C Indices correspondants au trois phases A, B, C
Angle d’observation de la matrice de PARK
Angle électrique rotorique
Angle électrique statorique.
Angle électrique de glissement.
V Tension.
I Courant.
Flux.
Résistance statorique
Résistance rotorique
inductance cyclique statorique et rotorique par phase.
Inductance mutuelle maximale entre une phase de stator et une phase
rotor
Inductance mutuelle cyclique
Couple électromagnétique.
Couple résistant
Moment d’inertie de la partie tournante.
Coefficient de frottement visqueux
� Vitesse mécanique.
� vitesse de rotation statorique
� vitesse de rotation rotorique
Maximum de l’inductance mutuelle entre phase du stator et la phase
Correspondante du rotor
Nombre de paires de pôles.
Coefficient de frottement visqueux.
Terme de couple de frottement visqueux
Constante de temps statorique et rotorique.
Coefficient de dispersion de blondel.
( Axes correspondant au référentiel lie au stator.
( Axes correspondant au référentiel lie au champ tournant.
( Axes correspondant au référentiel lie au stator
(u ,v) Axes du système biphasé
P( Matrice de PARK.
P( Matrice de PARK inverse
[A] matrice d’évolution d’état du système
[B] matrice de système de commande.
[U] matrice de commande
S Surface de glissement
V Fonction de Lyapunov
Glissement
SSV System a structure variable
CSV Commande a structure variable
Sommaire
INTRODUCTION GENERALE ....................................................................................................02
CHAPITRE I : Modélisation de la machine asynchrone
I.1 Introduction ...................................................................................................................................03
I.2 Génialité sur les machines Asynchrones triphasés..................................................................................03
I.2.1 Définition de la machine asynchrone ............................................................................................... 03
I.2.2 Description du moteur asynchrone................................................................................................... 04
I.2.3 Modele généralise triphasé de [Link] ............................................................................................ 05
I.3 Modélisation du moteur asynchrone .......................................................................................................05
I.3.1 Hypothèses simplificatrices .............................................................................................................05
I.3.2 Equations électriques .......................................................................................................................06
I.3.3 Equations magnétiques ....................................................................................................................07
I.3.4 Equation mécanique ......................................................................................................................... 08
I.4 Transformation de PARK . .....................................................................................................................08
I.5 Application de la transformation de Park au modèle de la MAS ............................................................10
I.5.1 Equations électriques .......................................................................................................................10
I.5.2 Equations magnétiques ....................................................................................................................11
I.5.3 Equations mécaniques ......................................................................................................................11
I.6 Choix du référentiel ................................................................................................................................11
I.6.1 Référentiel immobile lie au stator .................................................................................................... 12
I.6.2 Référentiel immobile lie au rotor ..................................................................................................... 12
I.6.3 Référentiel immobile lie au champs tournant .................................................................................. 13
I.7 Estimateur de flux ...................................................................................................................................14
I.8 Modélisation de l’onduleur de tension ....................................................................................................15
I.8.1 Choix des composants ......................................................................................................................15
I.9 Représentation d’état du modèle de la MAS...........................................................................................16
I.10 Présentation du système simulé ............................................................................................................17
I.10.1 Résultats de simulation à vide et en charge ................................................................................... 18
I.10.2 Interprétation des résultats .............................................................................................................19
I.11Conclusion .............................................................................................................................................20
CHAPITRE II : Commande par linéarisation entrée-sortie du moteur asynchrone
II.1 Introduction ...........................................................................................................................................21
II.2 Historique de la commande ...................................................................................................................22
II.3 Avantage de la théorie entrée-sortie ......................................................................................................22
II.4 Outils mathématiques ............................................................................................................................22
II.4.1 Gradient .......................................................................................................................................... 22
II.4.2 Dérivée de Lie.................................................................................................................................23
II.5 La linéarisation entrée-sortie .................................................................................................................23
II.5.1 Principe de la technique de linéarisation entrées-sorties ................................................................23
II.5.2 Système à une entrée et une sortie .................................................................................................. 25
II.5.3 Système à plusieurs entrées et plusieurs sorties .............................................................................. 25
II.6 Application au modèle de la machine asynchrone.................................................................................27
II.6.1 Modélisation de la machine asynchrone ......................................................................................... 28
II.6.2 Choix des grandeurs de sortie .........................................................................................................28
II.6.3 Calcul du degré relatif ....................................................................................................................28
II.6.4 Linéarisation du système ................................................................................................................29
II.6.5 Synthèse des régulateurs .................................................................................................................30
II.7 Simulation et Interprétation de résultat ..................................................................................................30
II.7.1 Test de robustesse avec variation de la charge .............................................................................. 30
II.7.2 Test de robustesse avec variation de la vitesse .............................................................................. 32
II.7.3 Test de robustesse avec variation de la résistance rotorique..........................................................33
II.8 Conclusion .............................................................................................................................................34
CHAPITRE III : Commande par mode glissant
III.1 Introduction ..........................................................................................................................................35
III.2 Bref historique ......................................................................................................................................35
III.3 Théorie de la commande par mode glissement ....................................................................................36
III.3.1-Système à structure variable..........................................................................................................37
III.3.-2 Objectif de la commande par mode glissant ................................................................................ 37
III.4 Notion de base ......................................................................................................................................37
III.4.1 Conception de la commande par mode de glissement ...................................................................37
III.4.2 Bases mathématiques de la commande à structure variable ..........................................................37
III.4.3 Condition d’existence du mode glissant ........................................................................................ 38
III.4.4 Choix de la surface de glissement ................................................................................................ 39
III.5 La fonction de LYAPUNOV ................................................................................................................39
III.6 Phénomène de chattering.....................................................................................................................40
III.6.1 Fonction sat ................................................................................................................................... 40
III.6.2 fonction smooth ............................................................................................................................. 40
III.7 Domaines D’Application de la commande par mode glissant .............................................................41
III.8 Avantages et inconvénients de la commande par mode de glissement ................................................41
III.9 Expression analytique de la commande ...............................................................................................42
III.10 Commande non-linéaire par mode glissant ........................................................................................43
III.10.1 Modèle du moteur asynchrone .................................................................................................... 43
III.10.2 Choix des surfaces de glissements .............................................................................................. 44
III.10.3 Conditions de convergence..........................................................................................................45
III.10.4 Loi de commande par mode glissant ........................................................................................... 45
III.11 Bloc de simulation ..............................................................................................................................48
III.12 Simulation et interprétation du Résultat ...........................................................................................49
III.12.1 Test de robustesse avec variation de la charge ............................................................................ 49
III.12.2 Test de robustesse avec variation de la vitesse ................................... Erreur ! Signet non défini.1
III.12.3 Test de robustesse avec variation de la résistance rotorique ............. Erreur ! Signet non défini.2
III.13 Conclusion ...............................................................................................Erreur ! Signet non défini.3
CONCLUSION GENERALE ................................................................ Erreur ! Signet non défini.4
Introduction Générale
Introduction Générale
Les technologies modernes des systèmes d’entraînement exigent de plus en plus un contrôle précis
et continu de la vitesse et du couple, tout en garantissant la stabilité, la rapidité et le rendement le
plus élevé possible.
Le moteur à courant continu, a satisfait une partie de ces exigences mais il est pourvu des balais
frottant sur le collecteur à lames, ce qui limite la puissance et la vitesse maximale et présente des
difficultés de maintenance et des interruptions de fonctionnement. Alors que le prix des machines
électriques varie peu, celui des composantes électroniques et microinformatique baisse
constamment, de telle façon que la part du variateur dans le coût d’un entraînement à vitesse
variable diminue. Pour toutes ces raisons, le moteur à courant continu trouve des limites
d’utilisation, le moteur à courant alternatif à tendance de le remplacer dans nombreuses
applications.
Actuellement la machine asynchrone est de plus en plus utilisée, elle présente l’avantage d’être
robuste, peut coûteuse. Mais malgré tous les avantages cités précédemment, la commande des
machines asynchrones pose de problèmes du fait que son modèle de base est non linéaire et
fortement couplé. Aussi ce qui complique ce modèle, c’est que les paramètres du moteur
asynchrone sont connus approximativement et peuvent varier avec le temps.
Les recherches ont donné un essor appréciable à la commande non-linéaire [16]. Parmi les
techniques de commande à grande performance assurant un découplage global entre les sorties à
commander quelque soit les profils de trajectoires imposées à la machine, on trouve la commande
non linéaire développée par "Isodori" , cette commande présente l'avantage de pouvoir commander
séparément le flux et la vitesse par un retour d'état non linéaire, chacun de ces sous systèmes
représente une boucle indépendante de commande d'une variable donnée (vitesse, couple, flux,
etc.).
Toutefois, cette technique de commande précitée présente relativement une certaine sensibilité liée
aux variations paramétriques. En effet, elle est dépendante directement du modèle de la machine, la
robustesse de son algorithme de commande est remise en question. Plutôt, il est important d'utiliser
des méthodes de contrôle robuste, soit linéaire ou non linéaire, telle que la commande par mode
glissant.
La commande par mode glissant qui a été développée en Union soviétique il y a plus de 30 ans,
fait partie de ces méthodes de commandes robustes. Elle possède des avantages incontestables pour
1
Introduction Générale
Dans le deuxième chapitre, nous rappelons la théorie de la commande par linéarisation entrée-sortie
et nous faisons la synthèse de cette technique sur le modèle de la machine asynchrone. Une
réalisation sur Sumilink qui nous permet de confirmer la validité de cette méthode
Dans le troisième chapitre, on traitera des rappels de notions de base sur la théorie de la commande
par mode de glissement, après on passera à l’application de cette commande sur la machine
asynchrone. ainsi une présentation des résultats de simulation sera donnée pour voir les
performances de cette technique et enfin une conclusion générale sur les résultats obtenus sera
donnée à la fin de ce travail.
2
CHAPITRE I
Modélisation
de la Machine
Asynchrone
CHAPITRE I Modélisation de la machine asynchrone
I.1. Introduction :
La modélisation des machines électriques est d’une grande importance aussi bien pour
le concepteur que pour l’automaticien. Au niveau de la conception,les équations de maxwell
offrent une analyse parfaite du comportement de la machine alors qu’au niveau de la commande,
un modèle basé sur les équations électriques est suffisant pour sa synthèse.
Il faut noter cependant qu’il existe trois approches concernant la modélisation des
machinesélectriques [1]:
- La modélisation de PARK basée sur les équations électriques des machines avec un certain
nombre d’hypothèses simplificatrices conventionnelles; ce modèle est bien adapté à l’élaboration
d’algorithme de commande.
-La modélisation par réseau de permanence ou le circuit magnétique est modélisée par un circuit
électrique équivalent. Elle peut être utilisée pour tester la robustesse des algorithmes de
commande, comme elle peut contribuer à l’estimation des paramètres de lamachine.
-La modélisation par éléments finis, c’est la modélisation la plus précise, elle est surtout utilisée
pour la conception des machines entre autres pour leurs dimensionnements.
Dans ce chapitre, nous présenterons le modèle mathématique triphasé de la machine asynchrone
(MAS) avec ses hypothèses simplificatrices et sa transformation dans le systèmebiphasé.
On désigne par " " le nombre de pairs de pôles de la machine. Lorsque le rotor tourne à une
vitesse différente de (asynchrone), l’application de la loi de Faraday aux enroulements
rotorique montre que ceux-ci deviennent le siège d’un système de forces électromotrices
triphasées engendrant elles-mêmes trois courants rotorique, d’après la loi de Lenz, ces derniers
s’opposent à la cause qui leur a donné naissance, c’est-à-dire la vitesse relative de l’induction
tournante statorique par rapport au [Link] va entraîner le rotor vers la poursuite du champ et
essayer de le faire tourner à la même vitesse )inférieur , il n’est pas au synchronisme
-L’additivité du flux.
-La constance des inductances propres.
-La loi de variation sinusoïdale des inductances mutuelles entre les enroulements statoriques
et rotoriques en fonction de l’angle électrique entre leurs axes magnétiques[4].
I.3.2Equations électriques
Les équationsde tension des trois phases statoriques et des trois phases rotoriques
sont:
(I-01)
(I-02)
En désignant par:
[ ] [ ][ ] [ ] (I-03)
[ =[ ]. [ ]+ (I-04)
Pour le rotor:
[ ] [ ][ ] [ ] (I-05)
[ =[ ][ ]+ (I-06)
I.3.3Equations magnétiques
Les relations entre les flux et les courants du stator et du rotor s’écrivent comme suit :
Pour le stator :
[ ] =[ ][ ]+[ ][ ] (I-07)
Pour le rotor :
[ ] =[ ][ ]+[ ][ ] (I-08)
Ou:
[ ] (I-09)
[ [ ] (I-10)
Ainsi:
[ ]= [ ]ᵀ= (I-11)
[ ]
Où :
+ [ ][ ] (I-12)
(I-13)
Avec :
:Moment d’inertie des masses tournantes.
:Couple résistant imposé à l’arbre de la machine.
Ω : Vitesse rotorique.
: Couple électromagnétique.
: Coefficient de frottement visqueux.
(I-14)
√ √ √
[ ]
[ ] (I-15)
Tel que :
Avec :
=∫ : est une position quelconque d’observation entre les systèmes d’axes biphasé
par rapport au système d’axes triphasé.
Pour le passage du système triphasé vers le système biphasé, on a les équivalents suivants :
L’équivalent de tension :
[ ] =[P( )] [ ] (I-16)
L’équivalent de courant :
[ ] =[P( )] [ ] (I-17)
L’équivalent du flux :
[ ] = [P ( )] [ ] (I-18)
Où :
[ ] = [p ( ) [ ] (I-19)
[ ] = [p ( ) [ ] (I-20)
[ ] = [p ( ) (I-21)
– (I-22)
(I-23)
(I-24)
(I-25)
= +M (I-26)
= +M (I-27)
= +M (I-28)
= +M (I-29)
Avec:
I.5.3Equations mécaniques
Le couple électromécanique représenté dans l’équation (I-12) devient :
= ( − ) (I-30)
(I-31)
Figure (I-04): Définition des axes réel de la MAS par rapport au référentiel (α,β)
(I-32)
C’est le repère le mieux adapté pour travailler avec les grandeurs instantanées, il possède
des tensions et des courants réelles et peut être utilisé pour étudier les régimesde démarrage et de
freinage des machines à courant alternatif
(I-33)
Figure (I-05) : Définition des axes réels de la MAS par rapport au référentiel (x ,y)
(I-34)
Ce système est utilisé pour étudier les processus transitoires dans les machines synchrones
et asynchrones.
I.6.3Référentiel immobile au champ tournant
(I-35)
q
d
Figure (I-06) : Définition des axes réels de la MAS par rapport auréférentiel (d
q)
(I-36)
Avec :
(I-37)
[Link] de flux
A partir des deux dernières équations du système(I-32) on réalise un estimateur de flux suivant:
{ (I-38)
Pour un onduleur triphasé, les commandes des interrupteurs d’un bras sont complémentaires.
Donc pour chaque bras il y a deux états indépendants. Ces deus états peuvent être considérés
comme une grandeur booléenne.
• =1 : interrupteur du demi- bras haut (a,b ou c) fermé.
• =0 : interrupteur du demi-bras bas (a,b ou c) fermé
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.
• La charge triphasée; est équilibrée en étoile avec neutre isolé.
Pour les tensions simples on a :
(I-39)
{ (I-40)
{ (I-41)
sont les tensions simples de la machine et est la tension fictive entre le neutre
de la MAS et le point fictif d'indice ’’o’’
(I-42)
(I-43)
(I-44)
(I-45)
[ ]= [ ][ ] (I-46)
[A] = (I-48)
[ ]
et la matrice de commande :
B= (I-49)
[ ]
: Coefficient de dispersion
La simulation des modèles des machines électriques permet d’analyser les résultats des systèmes
avec beaucoup de précision.
Dans notre cas, toutes les simulations sont effectuées par le logiciel «SIMULINK» sous
«MATLAB», il offre également plusieurs méthodes de résolution numérique des équations
différentielles. Il comporte des modules qui peuvent réaliser des fonctions comme
la sommation, l’intégrale, la dérivée...etc
20 200
10 150
vitesse(rad/s)
isa(A)
0 100
-10 50
-20 0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
t(s) t(s)
1.5 25
1 20
0.5 15
couple(N.m)
flux(web)
0 10
-0.5 5
-1 0
-1.5 -5
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
t(s) t(s)
20 2
10 1
flux(web)
isa(A)
0 0
-10 -1
-20 -2
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
t(s) t(s)
30 200
20 150
vitesse(rad/s)
couple(N.m)
10 100
0 50
-10 0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
t(s) t(s)
1.9. Conclusion
Dans ce travail, nous avons établi le modèle mathématique de la machine asynchrone alimentée
en tension. La complexité de ce modèle a été réduite par un certain nombre d'hypothèses
simplificatrices et par l'application de la transformation de PARK qui ramène le modèle triphasé
en un modèle diphasé équivalent.
Cette modélisation nous a permis essentiellement de voir les caractéristiques de la machine
asynchrone, de valider ainsi le modèle de Parkpar simulations effectuées sous MATLAB.
Nous avons remarqué que la vitesse de rotation du moteur chute lorsque on augmente la
charge, ceci nécessite sa régulation en particulier lorsqu'il s'agit d'un processus industriel
exigeant une vitesse constante indépendamment de la variation de la charge.
II.1 Introduction
L'intérêt d'améliorer les performances des systèmes commandés conduit à des
modélisations de plus en plus précises. Mais, si un modèle rend compte du comportement d’un
système dans une large plage de fonctionnement, il est malheureusement le plus souvent non-
linéaire et les outils fondamentaux de synthèse de lois de commande dans le domaine linéaire
deviennent insuffisants ce qui oblige d'avoir recours à de nouvelles méthodes.
En effet, plusieurs techniques de commande ont été développées (commande à flux orienté,
commande prédictive etc.). Les commandes conçus restent, cependant, très dépendant du point
de fonctionnement considéré et surtout du taux de variation des non linéarités du modèle du
moteur. Si ces non-linéarités deviennent très prononcées, la linéarisation de premier ordre n'est
donc plus valables et les commandes perdent automatiquement leurs performances et peuvent
rendre le système instable. Cette situation a poussé les chercheurs à s'orienter de plus en plus
vers les techniques de la commande non-linéaire basée sur la théorie de la géométrie
différentielle. Cette commande a été introduite principalement pour remédier aux problèmes
rencontrés avec la commande linéaire. Parmi ces méthodes, on trouve la technique de
linéarisation au sens des entrées - sorties. Son principe consiste à trouver une transformation qui
permet de compenser les non-linéarités du modèle et ainsi rendre la relation entre la sortie d'un
système et son entrée complètement linéaire. Ainsi dans la théorie E/S, on cherche une
description moins détaillée du système, qui met délibérément à part de l’information pour
produire des résultats plus généreux et plus faciles d’utilisation. Ceci présente bien évidemment
des avantages pour l’étude de la robustesse des systèmes de contrôle, dans les quels seuls une
structure approximée du modèle est disponible.[17]
Dans ce chapitre, Nous présentons le principe de la technique de linéarisation au sens des entrées
sorties. Nous allons concevoir un contrôleur non-linéaire basé sur cette technique afin de réguler
la vitesse et le flux du moteur asynchrone. A la fin du chapitre, nous allons montrer les
limitations et les faiblesses de la commande non-linéaire lorsqu'on est en présence d'incertitudes
dans les paramètres du moteur.
h(x)défini par .D'une façon similaire, le gradient d'un champ de vecteur f(x) est défini
= .
= h = [ ]
De plus, si g est un autre champ de vecteurs alors la fonction scalaire est donnée par:
Le but de cette technique est de transformer le système multi entrée non linéaire en un
système linéaire en utilisant un retour d'état linéarisante avec découplage entrée-sortie. De là ,on
pourra appliquer la théorie des systèmes linéaire.
Soit un système d'ordre n, multi entrée et multi sortie, décrit par la représentation d'état non
linéaire suivante :
(II-01)
Avec :
u et y : vecteur de commande (d'entrée ) et vecteur de sortie , respectivement.
ƒ(x), g(x) : champs de vecteur de fonctions continues.
h(x) : vecteur de sortie.
Les éléments des champs vectoriels ƒ, g et h sont des fonctions lisses. Si l'on considère le cas des
systèmes avec m entrée et m sortie, on cherche un bouclage statique de la forme :
tel que le comportement entrée- sortie du système (II-01) après bouclage soit linéaire et
découplé. Ainsi on obtient un ensemble de m sous systèmes mono sortie indépendants où les
entrées du sous-système i n'affectent pas la sortie et réciproquement.
Avec :
v : nouvelle variable de commande du système linéaire, v∈ .
β : Matrice non singulière de dimension m×m.
α : Vecteur de dimension m×1.
On utilise souvent des notions de la géométrie différentielle telle que la dérivée de Lie.
La nouvelle commande v permet de ramener le comportement entrée-sortie du système, définit
par l'équation (II-01) à celui d'un système linéaire, par différentiation des sorties du système
jusqu'à l'apparition des anciennes commandes en utilisant la dérivéede Lie.
∑ ( ) j=1,2,..p (II-04)
Où :
sont les dérivées de Lie de respectivement dans la direction de f et g.
est le nombre de dérivées nécessaires pour qu'au moins, une des entrées apparaisse dans
l'expression (II-04) et il est connu sous le nom du degré relatif correspondant à la sortie .
Le degré relatif total (r) est défini comme étant la somme de tous les degrés relatifs obtenus à
l'aide de (II-04) et doit être inférieur ou égal à l'ordre du système:
∑
On dit que le système (II-02) a pour degré relatif (r) s'il vérifie:
0 < k< ,1 ,1
Dans le cas où le degré relatif total est égal à 1'ordre du système, on est en présence d'une
linéarisation au sens des entrées-états. Si par contre le degré relatif total est strictement inférieur
à l'ordre du système, la linéarisation est dite linéarisation au sens des entrées- sorties. Pour
trouver l'expression de la loi linéarisante u qui permet de rendre la relation linéaire entre l'entrée
et la sortie, on récrit I’ expression (II-04) sous forme matricielle:
[ ] (II-05)
(x) = [ ] (II-06)
D(x)= (II-07)
[ ]
[ ] [ ] [ ]
.(- (x)+ ∑
[ ]
En remplaçant (II-08) dans (II-02), le système équivalent devient linéaire et totalement découplé
de la forme:
(II-09)
(II-11)
� (II-12)
[ ( ) ]
Posons :
., , ., , .
(II-14)
[ ]
Afin d'atteindre les objectifs de contrôle qui consistent en la commande des deux grandeurs
flux et vitesse, on choisit les variables de sorties suivantes :
Y=[ ]= (II-15)
�
La sortie (II-16)
La sortie (II-17)
Avec :
(II-18)
(II-19)
Où : . . .
Le degré relatif associé aux grandeurs de sortie et sont respectivement =2 et =2.
L'ordre de la dynamique interne associé à cette linéarisation est égale à :
n: étant l'ordre du système à contrôler (n= 5).
Cette dynamique interne est inobservable à partir des sorties choisies par la loi de
commande. Ainsi il devient impératif de vérifier sa stabilité. Dans le choix de la position du flux
rotorique par rapport au stator a été effectué comme coordonnée. La position est une grandeur
bornée entre- ce qui garantie la stabilité de la dynamique inobservable.
II.6.4-Linéarisation du système :
Pour linéariser la dynamique entrée-sortie du moteur, on considère seulement les dérivées
des sorties. On obtient :
[ ]= (II-20)
Avec :
E(x)= =( ) (II-21)
Donc :
[ ] 2 (II-22)
Le déterminant de la matrice E(x) est toujours différent de zéro, donc E(x) est une matrice
inversible. Ainsi, la loi de commande pour la linéarisation et le découplage entrée-sortie par
retour d'état du moteur set donnée par :
(II-23)
{ (II-24)
Le système découplé peut être représenté par deux intégrateurs d'ordre '2'.
( ) (II-25)
(� � ) )+ (II-26)
La présence des ondulations sur la réponse du couple électromagnétique est liée aux
harmoniques de courant.
1.03 1.4
ZOOM de flux
1.02
1.2
1.01
1
1
flux(web)
flux(web)
0.99 0.8
0.98
0.6
0.97
0.4
0.96
0.95 0.2
0.94
0
2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
15 14
ce
12 cr
10
10
5
couple(N.m)
8
isa(A)
0 6
4
-5
2
-10
0
-15 -2
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
200
v
180
v ref
160
140
vitesse(rad/s)
120
100
80
60
40
20
-20
0 1 2 3 4 5 6 7 8 9 10
temps(s)
15 1.4
1.2
10
1
5
flux(web)
0.8
isa(A)
0
0.6
-5
0.4
-10
0.2
-15 0
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
160 14
140 v 12 ce
v ref cr
120 10
100 8
vitesse(rad/s)
couple(N.m)
80 6
60 4
40 2
20 0
0 -2
-20 -4
-40 -6
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
1.12
ZOOM de flux
1.1
1.08
1.06
flux(web)
1.04
1.02
0.98
0.96
0.94
0.92
1 2 3 4 5 6 7 8 9 10
temps(s)
10 1.2
1
5
flux(web)
isa(A)
0.8
0
0.6
-5
0.4
-10
0.2
-15 0
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
200 14
v ce
180
v ref 12 cr
160
10
140
vitesse(rad/s)
couple((N.m)
120 8
100
6
80
60 4
40
2
20
0
0
-20 -2
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
1.06
ZOOM de flux
1.04
1.02
flux(web)
0.98
0.96
0.94
0.92
2 3 4 5 6 7 8 9
temps(s)
II.8- CONCLUSION
Dans ce chapitre, nous avons présenté une étude théorique sur les notions appropries à la
commande non linéaire par le biais de la linéarisation entrée-sortie et son application sur le
modèle de la machine. Nous avons vu que le modèle réduit de la machine asynchrone, dans le
référentiel lié au stator avec la vitesse et le flux comme grandeurs de sorties, admet une
linéarisation entrée-sortie. Les résultats de simulation obtenus montrent que cette commande
permet de découpler le réglage de la vitesse et du flux. De même, le modèle ainsi linéaire est
remplie au moyen de la technique de régulation approprie aux systèmes linéaires telle que
placement de pôles.
III.1 Introduction :
Dans le domaine de la commande des machines électriques, les travaux de la
recherches’orientent de plus en plus vers l’application des techniques de commandes robustes.
Cestechniques évoluent avec l’évolution des calculateurs numériques etde l’électronique de
puissance. Ceci permet d’aboutir à des processus industriels de hautesperformances. Nous
pouvons citer à titre d’exemple la commande floue, la commandeadaptative, la commande par
mode de glissement...etc. L’intérêt récent accordé à cettedernière est dû essentiellement à la
disponibilité des interrupteurs à fréquence de commutationélevée et des microprocesseurs de
plus en plus performants [2].
Le mode de glissement est un mode de fonctionnements particulier des systèmes à structure
variable. La commande par mode de glissement est une méthode de réglage danslaquelle la
structure est commutée de manière à ce que le vecteur d’état "x" suit une trajectoire S(x)=0 dans
l’espace d’état [8].Cette technique de commande s’adapte aux variations des paramètres du
procédé àcommander. Dans ce chapitre, nous allons appliquer la commande par mode glissant à
la machine asynchrone dans le but d'obtenir de bonnes performances même en présence des
perturbations.
avec une constante [Link] ces approches ne fournissent aucune garantie quant à la
convergence exacte au mode glissant et présente un dilemme chattering- robustesse.
Control de véhicule
Moteur électrique
Equation du mouvement
Découplage
Invariance
Equation d existence
III.4Notion de base
Avant d’entamer les méthodes utilisées dans la synthèse de contrôle par mode de glissement,
nous allons expliquer quelques notions de bases du régime glissant [11].
III.4.1-Conception de la commande par mode de glissement
La conception des régulateurs par les modes glissants prend en charge les problèmes de
stabilité et des performances désirées d’une façon systématique. La mise en œuvre de cette
méthode de commande nécessite principalement trois étapes
1-Le chois de la surface ;
2-L’établissement des conditions d’existence de la convergence ;
3-La détermination de la loi de commande.
III.4.2-Bases mathématiques de la commande à structure variable
La modélisation mathématique de la commande à structure variable conduit à des équations
différentielles de la forme :
(III-01)
Avec :
: vecteur d'état,
f : vecteur en fonctions de x et t
CSE. Université de M’sila 37
CHAPITRE III commande par mode glissant
{ i=1,2….m ; (III-3)
Dans le cas multi variables le régime glissant s’effectue sur une variété S(x)=0 de dimension (n×m) ,
S( x)=0 a �=� �
b e
�=� �
Les commutations ont lieu sur les droites décalées parallèlement de ±ΔS .Une trajectoire avec
h
trajectoire qui respecte cette condition. La fréquence de commutation est donc infiniment élevée.
[6].
{ (III-04)
S(x)=[( � (III-05)
Avec :
� = – : écart de la variable à régler.
: Constante positive qui interprète la bande passante du contrôle désiré.
r : degré relatif.
L’objectif de la commande est de maintenir la surface à zéro. Cette dernière est une
équation différentielle linéaire dont l’unique solution est e(x) [Link] un choix convenable du
paramètre, ceci revient à un problème de poursuitede trajectoire qui est équivalant à une
linéarisation exacte de l’écart tout en respectant la condition de convergence [7][10][13].
(III-06)
Et sa dérivée par :
(III-07)
Pour que la fonction de Lyapunov décroisse, il suffit d’assurer que sa dérivée est négative. Ceci
est vérifie si :
(III-08)
Cette équation montre que le carré de la distance vers la surface, mesuré par S 2 ( x) , diminue tout
le temps, contraignant la trajectoire du système à se diriger vers la surface des deux cotés .Cette
condition suppose un régime glissant idéal [14].
III.6 Phénomène de chattering
Le phénomène de Chattering est provoqué par une commutation non infiniment rapide de
la commande quand les techniques des modes glissants sont utilisé[Link] phénomène est
indésirable car il ajoute au spectre de la commande des composantes de haute fréquence. Le
broutement (phénomène de Chattering) peut être réduit en remplaçant la fonction« sign » par une
fonction de saturation adéquate qui filtre les hautes fré[Link] donne ci-dessous des
exemples de ces fonctions [15][4].
III.6 1 La fonction SAT :
�
� { � (III-09)
� | |
SAT ( S )
1
S
1
Figure (III-03):Fonction«SAT »
| |
SMOOTH (S )
1
1
Parmi les solutions proposées à ces problèmes on peut citer la commande par mode glissant à
bande limite qui consiste à remplacer la fonction de commutation dans la commande par une
fonction de saturation (hystérésis, arc tangente,… etc.)[16]
III.9-Expression analytique de la commande
Nous nous intéressons au calcul de la commande équivalente et par la suite au calcul de
la commande attractive du système défini dans l’espace d’état par l’équation (III-01).
ẋ(t) f (x, t) g(x, t)u(t) (III-10)
Le vecteur u est composé de deux grandeurs : et soit :
Nous avons :
{ + (III-12)
{ } { } (III-13)
= (III-14)
La solution la plus simple est de choisir sous la forme de relais (Figure III-4) .Dans ce
cas, la commande s’écrit comme suit :
(III-16)
+k
S( X )
-k
Où le facteur est toujours négatif pour la classe des systèmes que nous considérons.
� �
� (III-18)
�
{ ( )
ẋ=f(x)+g(x).u(t) (III-19)
Avec :
U Vs Vs T
x is is r r r x1 x2 x3 x4 x5
T T
(III-20)
Les variables x sont composés de deux états électriques (is, is) et deux états magnétiques
(r, r ) et un état représentant la vitesse de rotation du rotor Ωr, gouverné par une équation
mécanique.
f et g sont des fonctions continues en fonction des variables x.
x1 x1 T x3 p x4 x5 Vs
x2 x2 x4 p x5 x3 Vs
r
Tr (III-21)
x3 x1 x3 p x4 x5
M 1
M
Tr Tr
x4 T x2 T x4 p x3 x5
1
x5 ( x2 x3 x1 x4 ) C r
r r
J
Avec :
1 1
L T ( )
Rs
s r
M
1 1
Ls
1 (III-22)
pM
J L
r
� � (III-23)
Si on introduit l’erreur en vitesse :
� (III-24)
On aura
= ( )+( ) (III-25)
Régulation du Flux Φr
� +� (III-26)
Avec l’erreur en flux :
CSE. Université de M’sila 44
CHAPITRE III commande par mode glissant
� (III-27)
On aura :
( ) (III-28)
Où :
K1 et K2 sont des gains positifs.
Le choix du gain K se fait de façon que la valeur désirée soit imposée à la sortie dechaque
régulateur
III.10.3 Conditions de convergence
Pour que les variables choisis convergent vers leurs valeurs de référence il faut que les deux
surfaces de glissement soient nulles.
S1 0 ( r r ref ) K1 ( r r ref ) 0.
d
S2 0
dt (III-29)
( r r ref ) K2 ( r r ref ) 0.
d
dt
Par conséquent pour une surface de glissement nulle (S1=0) la vitesse mécanique et le flux
rotorique convergeant exponentiellement vers leurs références. Alors pour suivre r r ef et r r ef , il
par mode glissant. Les états considérés pour le contrôle du moteur asynchrone sont la vitesse
et le flux rotorique ( r , r ) .
Pour des raisons de commodité de calcul on prendra :
r r2 r2 r2 x32 x42 (III-30)
2 x x 2 x x
r 3 3 4 4 (III-31)
Les surfaces de glissement seront alors :
S1 K1 ( x5 r r ef ) ( x 5
r r ef ).
S2 K 2 (r r r ef ) (r r r ef ).
(III-32)
J
x5 (x2 x3 x3 x2 x1 x4 x4 x1)
C r
2 M (x x x x ) 2 (x2 x2)
J
(III-34)
r 1 3 2 4 3 4
Tr Tr
2.M (x x x x x x x x )
r 1 3 3 1 2 4 4 2
Tr
2 2M
( (x1 x3 x2 x4) (x3 x4 ))
2
2 2
-
Tr Tr Tr
- K1 r r ef r r ef ( x3 Vs x4 Vs )
r
(III-35)
S 2 ( 2 M K ) 2 M ( M h ( K 3 ) h
2 Tr
2 r 3 2 1
2
Tr Tr Tr Tr
p x5 h2 ) K 2 r r ef r r ef M ( x4 Vs x3 Vs )
Tr
Tels que :
h1 x1 x3 x2 x4
h2 x2 x3 x1 x4
(III-36)
h3 x12 x22
Avec :
C p x5 (h1 r )
J r
1 K1
F1 ( K 1 ) h 2
Tr
r r ef
K1 1
- r r ef
(III-38)
2 M
F2 ( K 2 ) r M ( h3 ( K 2 ) h1
M 3
Tr Tr Tr
T K
p x5 h2 ) r 2 r r ef
Tr
2 2 r r ef
Sous la forme condensé (III-37) sera de la forme :
S F D V (III-39)
Où :
F1 x4 x3
F ; D
F M x M x
2 3 4
S v Sgn(S) (III-40)
En égalisant (III-39) et (III-40) on aura :
V D 1 v Sgn( S ) D 1 F (III-41)
Dans ce chapitre la loi de commande est donnée par :
V Veq Vc (III-42)
Avec :
Veq D 1 F
Vc D 1 v Sgn(s)
(III-43)
Vc v1 0 Sgn(S1)
D 1
V
0 v2 Sgn(S2)
(III-46)
c
Pour que la loi de commutation intervient de la loi de commande globale il faut choisir v1 et v2
suffisamment grand.
v1 F1
v2 F 2
(III-47)
D
La loi de commande existe à condition que la matrice D sera inversible :
2
M r 0 (III-48)
Alors :
x4
M r
x3
D
1
r
(III-49)
r M r
x 3 x4
15 1.4
1.2
10
1
5
flux(web)
0.8
isa(A)
0
0.6
-5
0.4
-10
0.2
-15 0
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
200 14
v ce
180
v ref 12 cr
160
10
140
vitesse(rad/s)
couple(N.m)
120 8
100
6
80
60 4
40 2
20
0
0
-20 -2
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
15 1.4
10 1.2
1
5
flux(web)
isa(A)
0.8
0
0.6
-5
0.4
-10
0.2
-15
0 1 2 3 4 5 6 7 8 9 10 0
0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
200 14
v ce
12
v ref cr
150 10
8
vitesse(rad/s)
couple(N.m)
100 6
2
50
0
-2
0
-4
-6
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
1.2
10
1
5
flux(web)
0.8
isa(A)
0
0.6
-5
0.4
-10 0.2
-15 0
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
14 200
ce v
180
12 cr v ref
160
10 140
vitesse(rad/s)
couple(N.m)
120
8
100
6
80
4 60
40
2
20
0
0
-2 -20
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
temps(s) temps(s)
Z OOM
1
isa(A)
-1
-2
-3
-4
-5
1000 500
500
Surface de vitesse
Surface de flux
0 -500
-500 -1000
-1000 -1500
0 2 4 6 8 10 0 2 4 6 8 10
temps (s) temps (s)
III.13 Conclusion
CONCLUSION GENERALE
La commande d'une MAS peut se faire suivant plusieurs techniques, chacune d'elles offre des
performances dynamiques et statiques bien définies avec des limites d'applications.
Le problème se pose dans le choix de telle ou telle méthode. Le recours à une méthode ou à l'autre
se fait normalement en fonction des contraintes du cahier de charge, aux quelles s'ajoutent parfois
les exigences nouvelles de l'économie d'énergie et de l'économie du matériel qui devraient être
prises en compte. L'objectif principal de ce mémoire était la synthèse des lois de commande par
non linéaire pour la commande de la machine asynchrone et allant de paire entre limitation
d’énergie aux bornes du stator et performances automatiques optimales. Deux commandes sont
considérées: la linéarisation entrée-sortie et la commande par mode glissant
La commande non linéaire par linéarisation au sens des entrées sorties a donnée des résultats
satisfaisants pour tout ce qui concerne la poursuite des variations de la vitesse de consigne, mais
avec des changements dans le couple résistant le flux présente une erreur statique.
La commande par mode glissant est plus robuste vis à vis des variations du couple de charge ou
des variation de la résistance rotorique mais souffre du broutement.
Enfin, de l’étude des deux techniques du contrôle du moteur asynchrone triphasé, on peut conclure
que la technique de commande par mode glissant a conduit à des bonnes performances, dans
beaucoup de cas, on obtient une meilleure qualité de réglage par rapport à la commande non
linéaire, donc la commande par mode glissant offre certains avantages :
1* Une robustesse par rapport aux variations des paramètres du système.
2* Un dynamique très performant "temps de réponse acceptable et erreur stationnaire pratiquement
nulle".
Les performances de cette technique dépendent d’un choix convenable des coefficients de la
surface du glissement et la rapidité de la réponse dépend du couple maximal que peut donner à la
machine.
54
Annexe(A)
Paramètres de la machine
Puissance nominale
Vitesse nominale
Résistance statorique
Résistance rotorique
[3] LAIEB Fadila ; GUERROUCHE Hakima " Etude comparative de la commande par mode
glissant et la commande vectorielle d'une machine asynchrone" Mémoire d’Ingéniorat
Université de M’SILA,2004.
[4] ACHOUR Abdelmalak ; AICHOUCHE Azedine " Réglage par mode glissant d’un machine
asynchrone" Mémoire d’Ingéniorat Université de M’SILA,2005.
[5] BELKAIBECH Ali ; MAALMI Ali " Commande a base d’un régulateur RST de la machine
a induction" Mémoire d’Ingéniorat Université de M’SILA,2006.
[6] YAHIAOUI Omar ; LABZA Tahar " Etude et simulation d’association onduleur de tension
-moteur asynchrone et réglage par mode glissant" Mémoire d’Ingéniorat Université de
M’SILA,2007.
[7] IHADDADENE Kerim ; CHAIBI Madani ; ABDELATIF Billal " Etude comparative entre
la commande pi et celle par mode glissant d’un moteur asynchrone triphasé" Mémoire
d’Ingéniorat Université de M’SILA,2009.
[8] KAHIL Zoubir ; LATRECHE Houari " Commandes par retour d'état non linéaire d'un
moteur asynchrone " Mémoire d’Ingéniorat Université de M’SILA , JUIN 2006.
[9] KEDJOUTI Seddik ; MESSAOUDI Chouaib " Réglage des puissances active et réactive de
la génératrice asynchrone a double alimentation par des régulateurs en mode glissant " Mémoire
d’Ingéniorat Université de M’SILA , JUIN 2005.
[13] AMGHAR Billal ; HEDID Abdelkader " Commande par Fuzzy Sliding Mode d'un
variateur de vitesse hydrodynamique " Mémoire d’Ingéniorat Ecole Nationale
Polytechnique2007.
[14] BELKHEiR Med elamine ; HAMDI Mohamed ; ROUMAIR Tariq " Synthèse de lois de
commande par floue mode glissement avec une surface de glissement non linéaire appliques a
système Hydraulique " Mémoire d’Ingéniorat Université de M’SILA,2010
[16] DJERIOUI Ali ; KHELIFI Faris "Commande et observation par mode glissant des
convertisseurs multicellulaires séries " Mémoire d’Ingéniorat Université de M’SILA,2009
[18] DEBBI Mouna ; BAKHTI Sarah "Commande vectorielle d'une machine asynchrone a cage
en présence des défauts rotoriques " Mémoire d’Ingéniorat Université de M’SILA,2006
MEMOIRE DE FIN D'ETUDES EN VUE DE L'OBTENTION DU DIPÔLME
DE MASTER EN GENIE ELECTRIQUE
SPECIALITE : COMMANDE DES SYSTEMES ELECTRIQUES
Thème :
Commandes non linéaires de la machine asynchrone
Résume
Les machines asynchrones constituent la plus grande partie des moteurs construits jusqu'à
nos jours. La lois de commande classique peut être insuffisante car elle est moins robuste et
moins précise. Pour cela nous avons utilisé d'autres techniques de l'automatique moderne dans
l'objectif est d'améliorer les performances et la robustesse de la machine même en présence des
perturbations et des incertitudes.
Parmi ces méthodes on a commandé la machine asynchrone, alimentée par un onduleur de
tension, par deux types de commandes différentes; la première commande est la linéarisation
entrée-sortie et la seconde utilise le mode glissant. Ces méthodes sont appliquées sur le modèle
de Park dans le référentiel lié au stator. Les résultats de simulations sont satisfaisants.
Mots clé :
Machine asynchrone, mode glissant, linéarisation entrée-sortie , commandes non linéaires
N°d’ordre 086