Commande Vectorielle de la MASDE par IA
Commande Vectorielle de la MASDE par IA
Mémoire de master
Option : réseaux électriques
Commande Vectorielle de la machine asynchrone double étoile par les techniques de
l’intelligence artificielle
Présenté par : Encadré par:
HAMDI .N
CHEBOUKI BILAL CO-Encadreur:
TAMRABET MOHAMMED EL AMINE DR. DJEBABLA .A
2020/2021
ÉTAT DE L’ART
PLAN DE TRAVAIL : DE LA MASDE Ch:01
01
MODÉLISATION ET
Ch:02 SIMULATION DE LA
MASDE 02
ETUDE ET
SIMULATION DE LA Ch:03
03 MASDE
COMMANDE
VECTORIELLE DE
Ch:04
LA MASDE
04 APPLICATION DE
RESEAUX
NEURONES SUR LA Ch:05
05 MASDE
Introduction :
Ls1,s1 L L
s1 s1,s2 s1,r I
s1
L L I
s2Ls2,s1 s2,s2 s2,r s2 Représentation des enroulements statoriques
r L L L Ir et rotoriques de la MASDE MASD
r,s1 r,s2 r,r
:
Equations mécaniques
j C em C r kf
t
P d d
C em i abc , s 1 L s 1, r iabc , r t iabc , s 2 L s 2 , r t iabc , r t
2 d d
Modèle biphasé :
• Transformation de Park est basée sur la transformation d’un système triphasé d’axes (a, b, c) en
système équivalent biphasé d’axe (x, y) et inversement :
cos
cos 4 /3 cos 2 /3
p ( ) 2
sin sin 4 /3 sin 2 /3
3
1 2 1 2 1 2
• chaque interrupteur est représenté par une paire transistor-diode qui est modélisé par deux
états définis par la fonction de connexion logique :
1 l'interrupteur i est fermé ( conduit, bloqué
f
i =
0 l'interrupteur i est ouvert ( bloqué, conduit)
f
v
as 1 E
2 1 1 1
v b s 1 3 1 2 1 f 2
v 1 1 2
f
cs 1 3
Si v as 1ref v p (t ) f 1 1 si non f 1 0
Si v bs 1ref v p (t ) f 2 1 si non f 2 0
Si v cs 1ref v p (t ) f 3 1 si non f 3 0
Vréf Vp
1
tension (v)
0
-1
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
t (s)
commande f1
1
-1
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
t (s)
500
vas1
-500
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
t (s)
Les figures représente la réponse de la MASDE alimentée par deux onduleurs de tension à commande MLI sinus-
triangle, , suivi de l’application des charge nominale Cr =10 N.m a temps t 3 s
300 50
40
250
30
200
CEM [N*m]
wr [rad/s]
20
150
10
100
0
50
-10
0 -20
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
temps [s] temps [s]
FLUX PHy FLUX PHx
0.7 0.4
0.6
0.2
0.5
0.4 0
P H y [W b r]
0.3
P H x [W b r]
-0.2
0.2
0.1 -0.4
0
-0.6
-0.1
-0.2 -0.8
-0.3
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
-1
temps [s] 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
temps [s]
COURANT IYS
COURANT IYS1
0 10
-5 5
-10 0
I Y S 1 [A ]
IY S [A ]
-15 -5
-20
-10
-25
-15
-30
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 -20
temps [s] 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
temps [s]
COURANT IXS1 COURANT IXS
8 0
-5
4
2
-10
IX S 1 [A ]
IX S [A ]
0
-2
-15
-4
-6
-20
-8
-10 -25
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
temps [s] temps [s]
Commande Vectorielle de la MASDE :
Principe de la Commande Vectorielle:
Le principe de la commande vectorielle (ou commande par flux orienté) consiste à orienter une
des composantes de flux statoriques, rotorique ou de l’entrefer sur un axe du référentiel tournant
à la vitesse ωs .Ce principe de commande d´écoulée, conditionnant le fonctionnement stable du
moteur asynchrone est la principale caractéristique de la commande vectorielle conduisant aux
hautes performances industrielles des entraînements asynchrones. supportant les perturbations
de la charge
• Procède d’orientation du flux :
Orientation du flux rotorique avec les conditions dr r qr 0
Orientation du flux statorique avec les conditions ds s qs 0
Orientation du flux d’entrefer avec les conditions dm m qm 0
• Méthodes de la Commande Vectorielle :
1. Méthode directe
2. Méthode indirecte
• Méthode indirecte :
• Le principe de cette méthode consiste à ne pas utiliser l’amplitude du flux rotorique mais
seulement sa position. Le descriptif “ méthode indirecte” signifie qu’on peut ´éliminer un
estimateur du flux mais elle exige la présence d’un capteur de position du rotor. Cette méthode
est sensible aux variations des paramètres de la machine
résultats de simulation :
Fonctionnement en charge :
• La figure (IV.6) représente l’évolution des caractéristiques de la MASDE avec la régulation de vitesse par
la méthode indirecte, suivi de l’application des charges Cr = 14N.m à t=5s, en imposant la vitesse de
référence = 300tr/mn.
300 60
250 50
200 40
Cem [N*m ]
wr [rad/s]
150 30
20
100
10
50
0
0
-10
-50 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]
COURANT IQS1 COURANT IQS2
16 16
14 14
12 12
10 10
IQ S 2 [A ]
IQ S 1 [A ]
6 6
4 4
2 2
0 0
-2 -2
-4 -4
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.6 0.8
1.4
0.6
1.2
0.4
P H Q [ W b r]
P H D [W b r]
1
0.2
0.8
0
0.6
-0.2
0.4
-0.4
0.2
0 -0.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]
les réseaux neurones :
• Fonctionnement d’un RNA :
• Le Fonctionnement d’un RNA est basé sur Le neurone formel qui calcule la somme de ses
entrées puis cette valeur passe à travers la fonction d'activation pour produire sa sortie.
• En règle générale, le calcul de la valeur de cette fonction peut se comporte en deux étapes :
• Une combinaison linéaire des entrées :
350
300
250
200
Vitesse (rad/s)
150
5 NEURONE
100
7 NEURONE
3 NEURONE
50
-50
0 0.5 1 1.5 2 2.5 3
Temps (s)
Testes de robustesse :
• Afin de tester la robustesse du contrôle proposé, nous avons étudié les variations et l’inversion de
la vitesse, Et les variations des charges résistantes
• résultats de simulation :
1. test 01 : démarrage a vide :
350 350
300 300
250 250
200 200
Couple (N.m )
Vitesse (rad/s)
150 150
100 100
50 50
0 0
-50 -50
0 0.5 1 1.5 2 2.5 3 0 0.5 1 1.5 2 2.5 3
Temps (s) Temps (s)
20 60
isd1
50
15
40
10
30
is d1
isq1 (A)
5
20
0 10
-5
0
-10
-10 0 0.5 1 1.5 2 2.5 3 3.5
0 0.5 1 1.5 2 2.5 3 temps(s) 4
x 10
Temps (s)
30 2
phird
phirq
20
1.5
10
flux rotorique
isa1 (A)
(wb)
-10 0.5
-20
0
-30
-40 -0.5
0 2 4 6 8 10 12 14 16 18
0 0.5 1 1.5 2 2.5 3 temps
x 10
4
Temps (s)
2. test 02 : demarrage a vide suivi des charges nominales :
Temps (Sec) =[0 1.5 1.5 2.5 2.5 ] ; Cr (N.m) = [ 0 0 14 14 0 ]
350
350
300
300
250
250
V ite s s e (ra d /s )
200
Couple (N.m )
200
150
150
100
100
50
50
0
0
-50
0 0.5 1 1.5 2 2.5 3
Temps (s)
-50
0 0.5 1 1.5 2 2.5 3
Temps (s)
2 15
10
1.5
flux rotorique (w b)
is a1 (A )
0
0.5
-5
0 -10
-15
-0.5 0 0.5 1 1.5 2 2.5 3
0 0.5 1 1.5 2 2.5 3 3.5 Temps (s)
Temps (s)
30 20
18
25
16
20
14
15 12
is q 1
10 10
is d 1
8
5
6
0
4
-5 2
0
-10
0 0.5 1 1.5 2 2.5 3 3.5
temps(s) 4
-2
x 10 0 0.5 1 1.5 2 2.5 3 3.5
temps(s) 4
x 10
Pour réaliser ce test, on a inversé le signe de la vitesse de +300 rad/sec à -50 rad /sec à -300 rad/sec à partir de t1
=0 sec et t2=1.5 sec et t3=4 sec . Successivement
300 60
200
40
Vitesse (rad/s)
100
Couple (N.m)
0 20
-100
0
-200
-300 -20
0 1 2 3 4 5 6 0 1 2 3 4 5 6
Temps (s) Temps (s)
15 20
10 15
5 10
isa1 (A)
isq1 (A)
0 5
-5 0
-10 -5
-15 -10
0 1 2 3 4 5 6 0 1 2 3 4 5 6
Temps (s) Temps (s)
2 1.5
1.5 1
phirq (Wb)
phird (Wb)
1 0.5
0.5 0
0 -0.5
-0.5 -1
0 1 2 3 4 5 6 0 1 2 3 4 5 6
Temps (s) Temps (s)
• 4. test 04: variation de la vitesse :
t =[0 1.5 1.5 2.5 2.5 3 3.5] Wref= [300 300 250 250 200 200 300]
60
300
40
Vitesse (rad/s)
Couple (N.m)
200
20
100
0
0 -20
0 0.5 1 1.5 2 2.5 3 3.5 0 0.5 1 1.5 2 2.5 3 3.5
Temps (s) Temps (s)
10 20
5 10
isd1(Wb)
isq1 (A)
0 0
-5 -10
0 0.5 1 1.5 2 2.5 3 3.5 0 0.5 1 1.5 2 2.5 3 3.5
Temps (s) Temps (s)
2 1
1.5 0.5
phird (Wb)
phirq (Wb)
1 0
0.5 -0.5
0 -1
0 0.5 1 1.5 2 2.5 3 3.5 0 0.5 1 1.5 2 2.5 3 3.5
Temps (s) Temps (s)
10
5
isa1 (A)
-5
-10
350 350
WRN
Cem RN
WPI
Cem PI
300 300
250 250
200 200
V ites s e (rad/s )
Couple (N.m )
150 150
100 100
50 50
0
0
-50
-50 0 0.5 1 1.5 2 2.5 3
0 0.5 1 1.5 2 2.5 3
Temps (s)
Temps (s)
2
err-PI
1.5 err-RN
0.5
erre ur d e v ite s s e
0
-0.5
-1
-1.5
-2
Les machines multiphasées sont utilisées dans les applications à haute puissance car
elles se caractérisent par leur efficacité et leur capacité à Répartition de puissance
la commande vectorielle indirect de la machine asynchrone à double étoile basée sur la
régulation neuronale a donné des résultats très performants.
Merci pour votre attention