Modélisation de la machine double étoile
Modélisation de la machine double étoile
I.1-INTRODUCTION :
Lamodél
i
sat
iond’
unsy
stèmequel
conqueesti
ndi
spens
abl
elor
squenousv
oul
ons
une commande particulière de ce systèmeetl
’
étudedur
égi
met
ransi
t
oir
e.
Lamachi
neasy
nchr
oneàdoubl
eét
oil
e(MASDE)n’
estpasuns
yst
èmesi
mpl
e,car
de nombreux phénomènes compliqués interviennent dans son fonctionnement,
commel
asat
urat
ion,l
’
eff
etdepeau….
etc.
[2]
Cependant
,nousn’
all
onspast
eni
rcompt
edecesphénomènes
,card’
unepar
t,l
eur
f
ormul
ati
on mat
hémat
ique est di
f
fici
l
e, d’
aut
re par
t,l
eur i
nci
dence sur l
e
comportement de la machine est considérée comme négligeable dans certaines
condi
ti
[Link] per
metd’
obt
eni
r des équations simples, qui traduisent
fidèlement le fonctionnement de la machine.
Le modèle de la machine asynchrone à double étoile est un système de neuf
équations différentielles dont les coefficients sont des fonctions périodiques du
temps. La résolution d’
unt
elsy
stèmeestdi
ff
ici
l
emêmeav
ecl
’
uti
l
isat
iondel
’
out
il
informatique.
L’
uti
l
isat
iond’
unet
ransf
ormat
iondi
t
edePAR K,quiestunchangementconv
enabl
e
desv
ari
abl
es,per
metdedét
our
nercet
te di
f
fi
cul
téetd’
obt
eni
runmodèl
efaci
l
ement
exploitable.
Dans ce chapitre, nous présentons la modélisation de la machine asynchrone à
double étoile basée sur la résolution des équations régissant leur fonctionnement en
régime dynamique.
L’
étudedel
amodél
i
sat
iondesmac
hinesél
ect
ri
quespr
ésent
eungr
and intérêt en
raison des impératifs techniques, conduisent à concevoir et exploiter le matériel au
voisinage de leurs limites technologiques.
Electromécanique -1- M’
sila 2006
Chapitre I Modélisation de la machine asynchrone double étoile
I.2-Description de la machine asynchrone double étoile (MASDE):
La machine asynchrone double étoile dont nous allons étudier la mise en équation
correspond à la structure de principe représentée par la figure (I-1) [3].
Bs1
Bs2
Ar
Br
r
As2
Cs1 As1
Cr
Cs2
Figure (I-1): Représentation des enroulements de la MASDE.
La MASDE est composé de deux enroulements au stator fixe: As1, Bs1, Cs1 et As2,
Bs2, Cs2 qui sont déphasés de
=30° entre eux, et de trois phases rotoriques
mobiles: Ar, Br, Cr.
Dans des encoches régulièrement réparties sur la face interne de chaque
enroulement sont logés trois enroulements identiques à p paires de pôles. Leurs
ax
es sontdécal
és ent
re eux d’
un angl
e él
ect
ri
que égalà 2(/3. Ces enroulements
constituent le primaire de la machine.
Le rotor est une armature cylindrique séparée des deux enroulements du stator par un entrefer
de faible é
pai
sse
ur,e
tdontl
’enr
oul
eme
ntc
ompor
tet
roi
spha
sesi
dent
ique
sdé
cal
éede2(
/3.
Ces enroulements représentent le secondaire de la MASDE.
Electromécanique -2- M’
sil
a2006
Chapitre I Modélisation de la machine asynchrone double étoile
La structure de ce rotor est supposée dans notre cas être un rotor à cage (barre
conductrice en aluminium intégrée aux tôles ferromagnétiques).
I.2.1-Principe de fonctionnement de la MASDE:
Les courants triphasés de fréquence «fs» al
i
ment
antl
el’
enr
oul
ement1dust
atorde
la machine, donnent naissance à un champ tournant à la vitesse de synchronisme
«Ns» [3]
Tel que:
fs
N s [tr/s]
p
Electromécanique -4- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Sachant que :
Vsa1
Vsa 2
Vs , abc 1
Vsb1
Vs , abc 2
Vsb 2
,
Vsc1
Vsc 2
De même :
isa1
isa 2
ira
isb1
[is,abc1] ,
isb 2
[is,abc2] ,
irb
[ir,abc]
i sc 1
i sc 2
i rc
Electromécanique -5- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
1 1
(Lr Lmr) Lmr Lmr
2 2
Lrr = Lmr
2
1 1
(Lr Lmr) Lmr
2
1 1
Lmr Lmr (Lr Lmr)
2 2
2π 4π
Lms cos( γ) Lms cos(γ ) Lms cos(γ )
3 3
Ls1s2= Lms cos (γ )
4π
Lms cos( γ) Lms cos(γ )
2π
3 3
2π 4π
Lms cos( γ ) Lms cos (γ ) Lms cos(γ )
3 3
2π 4π
Lsr cos( θr) r
Lsr cos(θ ) Lsr cos(θr )
3 3
Ls1r = Lsr cos (θ
r
4π
) Lsr cos( θr) Lsr cos(θr
2π
)
3 3
2π 4π
Lsr cos( θ r ) r
Lsr cos (θ ) Lsr cos(θr)
3 3
2π 4π
Lsr cos( θr-γ ) Lsr cos(θ r γ ) Lsr cos(θ r γ )
3 3
Ls2r = Lsr cos (θ
r γ
4π
) Lsr cos( θ r γ) Lrs cos(θ r γ
2π
)
3 3
2π 4π
Lsr cos( θ r γ ) Lsr cos (θ r γ ) Lsr cos(θ r γ )
3 3
Les matrices: [Ls2,s1]=[Ls1,s2]t, [Lr,s1]=[Ls1,r]t, [Lr,s2]=[Ls2,r]t
Ls1 , Ls2 , Lr : Inductances de fuites des deux enroulements statoriques et du rotor.
Lms , Lmr ,Lsr : Inductances mutuelles .
Lms= Lmr= Lsr =2/3 Lm.
I.2.3.2- Equation mécanique:
Lecoupl
eél
ect
romagnét
iqueestdonnéparl
’
expr
essi
on génér
ales
uiv
ant
e:
I t . .
Cem L .
I
p d
(I-4)
2
d
Electromécanique -6- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Avec:
L s 1 , s 1 L s1, s 2 L s1, r
is, abc1
L
Ls2, r
is , abc2
L s 2 , s 1 Ls2, s2 ; [I ]
Lr,r
L r , s 1
Lr, s2 ir , abc
Les sous matrices: Ls1,s1 , Ls2,s2, Ls1,s1, Ls2,s1, Lr,r, contiennent des termes constants.
On à donc
après simplification:
p d d
Cem = ( [is,abc1] [Ls1,r] [ir,abc] +[is,abc2] [Ls1,r] [ir,abc] ) (I-5)
2 dθ dθ
L’
équat
ionmécani
quedel
amac
hines’
écr
it
:
d r
J C em C r K f r (I-6)
dt
La résol
uti
onanal
yt
iquedesneuféquat
ionsdi
ff
érent
iel
l
esav
ecl
’
équat
ionducoupl
e
et de la vitesse est très compliquée car certains coefficients font intervenir des
fonctions sinusoïdales dues au mouvement de rotation du rotor. De ce fait, une
simplification adéquat
e s’
i
mpos
e,quiestl
atr
ans
for
mat
ion de c
e sy
stème parl
a
méthode de PARK.
I.3 Transformation de PARK :
Cette transformation consiste à transformer un système triphasé (a,b,c) en un
système diphasé équivalent (d,q),af
ind’
obt
eni
runmodèl
emathématique plus
simple que
q le modèle physique du système .[2],[3],[5]
Vsb1
d q
sb1
Vsb2 isq1
isb1 ra Vsa2
isa2 Vsq1
r
sb2 isb2 ira
r s2 sa2 isq2
Vsa1 Vsq2
s1 isa1
irq
sa1
irb irc
rb d
isc1 rc ird isd2 isd1
Vsc1 Vsd2 Vsd1
Electromécanique -7- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Sa matrice inverse:
1
cos( θs 1 ) - sin( θs 1 )
2
1 2 2π 2π 1
P( θs 1 ) cos( θs 1 - ) - sin( θs 1 - )
3 3 3 2
2π 2π 1
cos( θs 1
) s1
- sin( θ )
3 3 2
Pou
rl’
enr
oul
eme
nt2s
tat
ori
que
s:
2π 2 π
cos(θs 2) cos(θs 2 - ) cos(θs 2 )
3 3
2
2π 2 π
P(θs 2 ) sin(θs 2 ) - sin(θs 2 - ) i2 )
- sin(θ
3 3 3
1 1 1
2 2 2
Sa matrice inverse:
1
cos( θs 2 ) - sin( θs 2 )
2
2 2π 2π 1
P( θs 2 ) 1 cos( θs 2 - ) - sin( θs 2 - )
3 3 3 2
2π 2π 1
cos( θs 2 ) i2
- sin( θ )
3 3 2
Electromécanique -8- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Pour le rotor :
2π 2π
cos( θr ) cos( θr - ) cos( θr )
3 3
2 2π 2π
P( θr )
sin( θr ) - sin( θr ) - sin( θr )
3 3 3
1 1 1
2 2 2
Sa matrice inverse:
1
cos( θr ) - sin( θr )
2
2 2π 2π 1
P( θr ) 1 cos( θr - ) - sin( θr )
3 3 3 2
2π 2π 1
cos( θr
) - sin( θr )
3 3 2
dφsd 1
Vsd 1 [Link] 1 ωs
.φsq1
dt
dφsq 1
Vsq 1 [Link] 1 ωs
.φsd 1
dt
dφsd 2
Vsd 2 rs [Link] 2 ωs
.φsq 2 (I-7)
dt
dφsq 2
Vsq 2 rs [Link] 2 ωs
.φsd 2
dt
dφrd
Vrd [Link] (ωs- ωr
).
φrq
dt
dφrq
Vrq [Link] (ωs- ωr
).
φrd
dt
Electromécanique -9- M’
si
la2006
Pe=Vsa1.isa1+Vsb1.isb1+Vsc1.isc1+ Vsa2.isa2+Vsb2.isb2+Vsc2.isc2
= Vsd1.isd1+ Vsd2.isd2+ Vsq1.isq1+ Vsq2.isq2 (I-10)
En remplaçant (I-7) dans (I-10) on aura :
dφsd1 dφsd 2
Pe=( rs1.isd1 ωs
.φ
sq1 ).i
sd1+( rs [Link] 2 ωs
.φsq 2 ).isd2
dt dt
dφsq1 dφsq2
+( rs1.isq1 ωs
.φsd 1 ).isq1+ ( rs 2.isq2 ωs
.φ
sd 2 ).isq2
dt dt
2 2 2 2
Pe={rs1.(isd1) +rs2.(isd2) +rs1.(isq1) +rs2.(isq1) }
dφsd 1 dφsd 2 dφsq1 dφsq 2
+{ isd1+ isd2+ isq1+ isq2} + {s( φsd 1 .isq1+φsd2 .isq2 - φsq1 .isd1
dt dt dt dt
Electromécanique -10- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Donc:
Cem =p.(φsd1 .isq1+ φsd2 .isq2 - φsq1 .isd1 -φsq2 .isq2) (I-13)
Electromécanique -11- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Tel que:
Avec:
md=Lm. (isd1+isd2+ird)
mq =Lm. (isq1+isq2+irq)
Lescour
ant
sst
ator
iquesetr
otor
iquess’
écr
ivental
orscommesui
t
:
sd1 md
isd1=
Ls1
sd 2 md
isd2=
Ls 2
sq1 mq
isq1=
Ls1
sq 2 mq
isq2= (I-18)
Ls 2
rd md
ird=
Lr
rq mq
irq=
Lr
Sachant que
mq et md s’
écr
iventenf
onct
i
ondesf
luxdel
amachi
necommesui
t:
Electromécanique -12- M’
si
la2006
d sd 1 rs 1
V sd 1 ( sd 1 md ) ωs.sq 1
dt Ls 1
d sq 1 rs 1
Vsq 1 ( sq 1 mq) ωs.sd 1
dt Ls 1
dsd 2 rs 2
Vsd 2 (sd 2 md) ωs
.sq 2 (I-19)
dt Ls 2
d sq 2 rs 2
Vsq 2 ( sq 2 mq) ωs.sd 2
dt Ls 2
d rd rr
V rd ( rd md ) (ωs ωr
)rq
.
dt Lr
d rq rr
V rq ( rq mq ) (ωs ωr
)rd
.
dt Lr
Nousmet
tonsl
esy
stèmed’
équations (I-18) s
ousf
ormed’
unsy
stème d’
équat
ions
d’
état
.
X AX BU
X : c
Avec : Ve teurd’
éta
t
Ve
U : c
teurd’
entr
ée
X= [ φsd 1 , φsd 2 , φsq1 , φsq 2 , φrd , φrq ] t
Electromécanique -13- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
1 0 0 0 Vsd 1
0
0
1 0
Vsd 2
0 0 1 0 ; U
B Vsq 1
0
0 0 1
0 0 0 0 Vsq 2
0 0 0 0
I.5-Simulation numérique:
La figure (I.3) représente les performances de la conduite de la machine asynchrone
t
ri
phaséeàdoubl
est
atorl
orsd’
undémar
rageàv
ide.
La figure (I.4) représente les performances de la conduite de la machine asynchrone
t
ri
phaséeàdoubl
est
atorl
orsd’
undémar
rageav
ecappl
i
cat
i
ond’
uncoupl
ede
charge.
Figure (I-3) :Compor
tementdel
aMASDEl
orsd’
undémar
rageàv
ide.
Electromécanique -14 M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Figure (I-4) : Compor
tementdel aMASDEl or sdel ’appl
i
cat
iond’
uncoupl
ede
charge entre deux instants.
Electromécanique -15- M’
si
la2006
Chapitre I Modélisation de la machine asynchrone double étoile
Electromécanique -16- M’
si
la2006
Electromécanique -17- M’
si
la2006