SADR AJ - emp.etf.rsemp.etf.rs/radovi/Diplomski/branko_kovacevic/branko_kovacevic.pdf · elektromagnetna indukcija, zbog koje je potrebna razlika u brzinama izmeu rotora i obrtnog
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Oblast kojom se ovaj diplomski rad bavi su asinhroni motori i njihova primena ubrzinskim servomehanizmima. U prvom delu ovog rada bi e prikazan analiti ki modelasinhronog motora i ostalih komponenti koje se pojavljuju kada se asinhroni motor ( AM )koristi u spomenutoj primeni, a te komponente su trofazni tranzistorski invertor ( radi emo satrofaznim asinhronim motorom ), histerezisni strujni regulator, indirektna vektorska kontrola,brzinski regulator.
Kori enje AM u servomehanizmima povezano je sa odre enim te ko ama kojeproizilaze iz njegovih osobina . Osobina koja dovodi do tih te ko a je slede a : pri konstantnoju estanosti i amplitudi napona napajanja brzina rotora zavisi od momenta optere enja, toiziskuje komplikovane algoritme upravljanja u slu ajevima kada elimo da koristimo asinhronimotor u brzinskim ili pozicionim servomehanizmima, odnosno tamo gde se tra i preciznakontrola brzine i/ili polo aja. Ova pojava je posledica principa rada asinhronog motora, a to jeelektromagnetna indukcija, zbog koje je potrebna razlika u brzinama izme u rotora i obrtnogmagnetnog polja koje generi e stator da bi postojao elektromagnetni moment. Elektronika kojarealizuje pomenute algoritme ranije je bila skupa i to je ote avalo kori enje asinronih motora uovakve svrhe, ali se oni danas sa pojeftinjenjem elektronskih komponenti i kori enjem ra unarau realizaciji algoritama upravljanja sve vi e koriste.
Prvi zadatak koji ovaj diplomski rad treba da obavi je da se pored analiti kog prikazaAM tako e analiti ki predstavi i postupak koji omogu ava kori enje AM u servomehanizmima,a to je vektorska kontrola, u na em slu aju indirektna.
Drugi zadatak koji treba da se uradi je primena dobijenih rezultata, odnosno izlo eneteorije. Mi ne emo raditi na konkretnom asinhronom motoru, ve emo modelirati sve opisanekomponente na ra unaru, u programskom paketu SIMULINK. Da emo model jednogkonkretnog asinhronog motora, a to je Sever ZK-80. Napomenu emo koja se sve zanemarenjavr e pri modeliranju asinhronog motora : sve nelinearnosti kod motora, od kojih je najva nijanelinearna karakteristika magne enja, zatim se zanemaruju histerezisni gubici, povr inski efekat,gubici usled vihornih struja i zanemaruje se energija akumulirana u elektri nom polju.
Modeliranje na ra unaru se svodi na predstavljanje relacija izvedenih za AM u prvomdelu ovog rada uz pomo blok eme. Tako e e biti dat model invertora sa histerezisnimstrujnim regulatorom, model indirektnog vektorskog kontrolera, model brzinskog regulatora i
3
model inkrementalnog enkodera. Na kraju je sve ove komponente ( module ) potrebno povezatiu vektorski kontrolisan elektromotorni pogon sa brzinskom regulacijom.
Nakon pravljenja modela u SIMULINK-u mo e da se pokrene i simulacija, koja namomogu ava da pribli no vidimo kakvo bi bilo pona anje realnih komponenti. Rezultati tihsimulacija bi e dati ovde u vidu grafika karakteristi nih fizi kih veli ina.
Re ‘modularni’ u nazivu diplomskog rada poti e od toga to SIMULINK omogu avada se komplikovane blok eme, kao to je model AM, zamene jednim blokom, koji je potpunoekvivalentan toj emi, a to omogu ava pravljenje veoma komplikovanih modela bez gubljenjapreglednosti.
2. OSNOVNE TEORIJSKE POSTAVKE
U ovom poglavlju bi e re i o analiti kom modelu asinhronog motora do kojeg se dolaziuz navedena zanemarenja, zatim e biti re i o nekim transformacijama koje nam poma u daprika emo taj model u raznim koordinatnim sistemima i o relativizaciji jedna ina. Nakon togaemo predstaviti princip indirektne vektorske kontrole, koja omugu ava primenu asinhronih
motora u servomehanizmima i na kraju emo prikazati postupak odre ivanja parametarabrzinskog PI regulatora na osnovu eljenog pona anja sistema.
2.1. Model asinhronog motora
Po to analiziramo trofazni asinhroni motor, njegov stator ima tri faze ije su osenormalne na osu simetrije motora, a me usoban polo aj im je kao na slici 1. Ozna avamo ihslovima a, b, c. Na rotoru se, zbog same njegove konstrukcije, ne mogu razlikovati faze.Me utim, mo emo formalno usvojiti postojanje nekakvih faza ije ose imaju isti me usobnipolo aj kao i faze statora, ali koje rotiraju zajedno sa rotorom, tako da su nepomi ne u odnosuna namotaje rotora. Ozna i emo te faze sa A, B, C. I one su prikazane na slici 1. Ugao θm je
4
Slika 1. Raspored faza na statoru i rotoru
ugao za koji se obrne rotor od po etka obrtanja. Na slici 1. su ozna ene i ose α, α’, β i β’ ijie smisao biti obja njen u slede em odeljku. Analiti ki model asinhronog motora sastoji se od
nekoliko grupa jedna ina : jedna ine naponskog balansa za stator i rotor, definicione jedna ineflukseva, mehani ka jedna ina koja opisuje zavisnost promene ugaone brzine od rezultuju egmomenta i jedna ina koja odre uje aktivni, elektromagnetni moment.
Jedna ine naponskog balansa za stator pi u se za svaku fazu pojedina no i glase :
u R id
dta s aa= +
Ψ ,
u R id
dtb s bb= +
Ψ ,
u R id
dtc s cc= +
Ψ ,
(2.1)
pri emu je R s ( otpornost namotaja statora ) ista za sve namotaje zbog simetri nosti motora.Jedna ine naponskog balansa za rotor mo emo pisati u istom obliku kao za stator u
onom koordinatnom sistemu u kom se namotaji, odnosno ice rotora stvarno i nalaze, tj. ukoordinatnom sistemu koji se obr e zajedno sa rotorom. Izbor faza A, B, C zadovoljava tajuslov i mo emo pisati :
0 = = +u R id
dtA r AAΨ
,
0 = = +u R id
dtB r BBΨ
,
5
0 = = +u R id
dtC r CCΨ
,
(2.2)
gde je R r otpornost namotaja rotora. Naponi su jednaki 0 jer su namotaji rotora kratko spojeni. Jedna ine koje odre uju flukseve glase :
Ψa s a ab b ac c aA A aB B aC CL i M i M i M i M i M i= + + + + + ,
Ψb ab a s b bc c bA A bB B bC CM i L i M i M i M i M i= + + + + + ,
Ψc ac a bc b s c cA A cB B cC CM i M i L i M i M i M i= + + + + + ,
ΨA aA a bA b cA c r A AB B AC CM i M i M i L i M i M i= + + + + + ,
ΨB aB a bB b cB c AB A r B BC CM i M i M i M i L i M i= + + + + + ,
ΨC aC a bC b cC c AC A BC B r CM i M i M i M i M i L i= + + + + + ,(2.3)
pri emu koeficijenti imaju slede e zna enje : Ls i Lr su respektivno sopstvene induktivnosti
statorskih i rotorskih namotaja; koeficijentiM xy , me usobne induktivnosti namotaja statora, su
jednaki -1/2 Ls zbog me usobnog polo aja namotaja, koeficijenti M XY su analogno tome
jednaki -1/2 Lr a M xY su me usobne induktivnosti statorskih i rotorskih namotaja i promenljivisu u vremenu ako se rotor obr e, a prostoperiodi ni su ako se rotor obr e konstantnombrzinom. Njihova maksimalna vrednost je induktivnost magne enja i ozna ava emo je sa Lm .Jedan od razloga nepodesnosti ovog modela za analizu motora je vremenska promenljivostkoeficijenata M xY .
Mehani ka jedna ina glasi :
Jd
dtM Mm
em opt
ω= − ,
(2.4)
gde je J moment inercije motora, ωm ugaona brzina rotora, M em i M opt su respektivno aktivni
elektromagnetni i opteretni moment. U ovoj jedna ini smo zanemarili frikciju ije uzimanje u
obzir bi na desnoj strani jedna ine dodalo lan -Kf ωm . Izraz za M em emo izvesti kasnije.
Jo jedan razlog to rad sa ovim jedna inama nije pogodan je to je u u ua b c+ + = 0 pajedna ine naponskog balansa za stator, a to va i i za rotor, nisu nezavisne. Jedna mo e da seizrazi preko druge dve. To zna i da u tim jedna inama imamo zapravo dve koordinate stanja.Zato je pogodno pisati jedna ine u Dekartovom pravougaonom koordinatnom sistemu, jer onima dve koordinate koje su me usobno nezavisne.
2.2. Klarkova transformacija
6
Klarkova transformacija predstavlja prelaz sa tri faze, koje su realne, na dve fiktivnefaze, ije se ose nalaze pod uglom od π/2. Za po etak emo usvojiti odvojene parove faza zastator i rotor. Za stator usvajamo faze α i β koje su nepomi ne u prostoru. Osa faze α neka sepoklapa sa osom faze a. Za rotor usvajamo faze α’ i β’ koje se obr u zajedno sa rotorom
ugaonom brzinom ωm. Neka se osa faze α’ poklapa sa osom faze A. Sve ove faze tako e supredstavljene na slici 1. Klarkova transformacija predstavlja i prelaz u suprotnom smeru, tj. sadve faze na tri. Ozna imo sa V bilo koju od veli ina ( napon, struja, fluks ) koje karakteri usvaku fazu ponaosob. Na osnovu slike 1. lako se dolazi do slede ih veza , sa tim da je dodatkoeficijent k iji emo smisao naknadno objasniti :
V k V V V k VV V
a b c ab c
απ π
= + +⎛⎝⎜
⎞⎠⎟
= − −⎛⎝⎜
⎞⎠⎟
cos cos2
3
4
3 2 2 ,
( )V k V V k V Vb c b cβπ π= −⎛
⎝⎜⎞⎠⎟
= −cos cos6 6
3
2 ,
V k V V V k VV V
A B C AB C
απ π
' cos cos= + +⎛⎝⎜
⎞⎠⎟
= − −⎛⎝⎜
⎞⎠⎟
2
3
4
3 2 2 ,
( )V k V V k V VB C B Cβπ π
' cos cos= −⎛⎝⎜
⎞⎠⎟
= −6 6
3
2 .
(2.5)
Uvo enjem ovih fiktivnih faza umesto stvarnih znatno se olak ava analiza asinhronogmotora, a samim tim se upro ava i sinteza raznih upravlja kih algoritama. Da bismo objasniliulogu koeficijenta k, posmatra emo izraz za ulaznu snagu :
P = u i u i u ia a b b c c+ + .(2.6)
Ako izra unamo na osnovu u u ua b c+ + = 0 , i i ia b c+ + = 0 i relacija (2.5) izraz
u i u iα α β β+ , dobi emo da je
u i u iα α β β+ = 2
32k ( )u i u i u i k Pa a b b c c+ + = 2
32 .
(2.7)
Tako e je
V k VV V
kVab c
aα = − −⎛⎝⎜
⎞⎠⎟
=2 2
3
2 .
(2.8)
Iz relacija (2.7) i (2.8) mo emo do i do slede ih zaklju aka : ako elimo da namprimenjena transformacija bude invarijantna po snazi, odnosno da bude P
7
= u i u iα α β β+ , skalira emo jedna ine (2.5) sa koeficijentom k = 2 3/ ; ako elimo
invarijantnost faznih veli ina, tj. da amplitude napona uα i u a budu jednake, uze emo k = 2/3.
Mo e da bude i k = 1, ali tada bi izraz za ulaznu snagu bio P = 2/3 ( u i u iα α β β+ ).
Jedna ine naponskog balansa za stator u α, β sistemu glase :
u R id
dtsα αα= +
Ψ ,
u R id
dtsβ ββ= +
Ψ ,
(2.9)
a za rotor u α’, β’ sistemu :
0 = = +u R id
dtrα αα
' ''Ψ
,
0 = = +u R id
dtrβ ββ
' ''Ψ .
(2.10)
Jedna ine fluksnih obuhvata su :
Ψα α αα α αβ β= + +L i M i M is ' ' ' ' ,
Ψβ β βα α ββ β= + +L i M i M is ' ' ' ' ,
Ψα αα α βα β α' ' ' '= + +M i M i L ir ,
Ψβ αβ α ββ β β' ' ' '= + +M i M i L ir ,
(2.11)
pri emu su koeficijenti M xy' opet promenljivi u vremenu ako se motor obr e.
2.3. Parkova transformacija αβ/dq
Do sada izlo eni modeli asinhronog motora imaju nekoliko mana. Jedna od njih je to sejedna ine za stator i rotor ne pi u u jedinstvenom koordinatnom sistemu, a posledica toga je dase u jedna inama fluksnih obuhvata javljaju vremenski promenljivi koeficijenti. Drugiproblem je to su u stacionarnom stanju sve veli ine prostoperiodi ne, to ote ava sintezu bilokakvog regulatora, a samim tim i primenu motora u servomehanizmima. Prvi problem seotklanja prebacivanjem jedna ina za rotor iz α’β’ u αβ sistem. To se posti e na osnovu slike 1.U αβ sistemu uve emo i dodatni indeks koji mo e biti s ili r i ozna ava da li se ta veli inaodnosi na stator ili rotor. Veze izme u veli ina u ova dva sistema su :
8
0 = = −u u ur m mα α βθ θ' 'cos sin ,
0 = = +u u ur m mβ α βθ θ' 'sin cos ,
(2.12)
0
0
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ =
−⎡
⎣⎢
⎤
⎦⎥
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥
u
u
u
uT
u
ur
r
m m
m m
α
β
α
β
α
β
θ θθ θ
cos sin
sin cos'
'
'
'
.
(2.13)
Analogne jedna ine se mogu pisati za flukseve i struje rotora :
ΨΨ
ΨΨ
α
β
α
β
r
r
T⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥
'
'
,
(2.14)
i
iT
i
ir
r
α
β
α
β
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥
'
'
.
(2.15)
U vezi matrice T treba primetiti da je T TT = −1 . Jedna ine (2.10) mo emo predstaviti umatri nom obliku :
0
0
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥
u
uR
i
id
dtr
α
β
α
β
α
β
'
'
'
'
'
'
ΨΨ
(2.16)
i pomno i emo ovu matri nu jedna inu sa leve strane matricom T :
0
0
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥
u
uT
u
uR T
i
iT
d
dtr
rr
α
β
α
β
α
β
α
β
'
'
'
'
'
'
ΨΨ
,
(2.17)
u
uR
i
iT
d
dtT
r
rr
r
r
T r
r
α
β
α
β
α
β
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥
⎛
⎝⎜⎜
⎞
⎠⎟⎟
ΨΨ
,
(2.18)
( )u
uR
i
iT
d
dtT T T
d
dtr
rr
r
r
T r
r
T r
r
α
β
α
β
α
β
α
β
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥
ΨΨ
ΨΨ
,
(2.19)
9
u
uR
i
id
dtr
rr
r
r
r
rm
r
r
α
β
α
β
α
β
α
βω⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥ +
−⎡
⎣⎢
⎤
⎦⎥
⎡
⎣
⎢⎢⎢
⎤
⎦
⎥⎥⎥
ΨΨ
ΨΨ
0 1
1 0 .
(2.20)
I najzad, u skalarnom obliku prethodna jedna ina glasi :
0 = = + +u R id
dtr r rr
m rα αα
βωΨΨ ,
0 = = + −u R id
dtr r rr
m rβ ββ
αωΨ
Ψ .
(2.21)
Jedna ine naponskog balansa za stator su ve napisane u αβ sistemu, a mi emo ihponoviti uz dodavanje indeksa s :
u R id
dts s ss
α αα= +
Ψ ,
u R id
dts s ss
β ββ= +
Ψ .
(2.22)
Jedna ine fluksnih obuhvata sada, po to se pi u u jedinstvenom koordinatnom sistemuza stator i rotor, ne sadr e vremenski promenljive koeficijente :
Ψα α αs s s m rL i L i= + , Ψβ β βs s s m rL i L i= + ,
Ψα α αr r r m sL i L i= + , Ψβ β βr r r m sL i L i= + .
(2.23)
Pogodnost rada u αβ koordinatnom sistemu je u tome to u jedna inama fluksnihobuhvata imamo konstantne parametre. Mana analize modela asinhronog motora u αβ sistemuje to su u stacionarnom stanju, kada se rotor obr e konstantnom brzinom, sve karakteristi neveli ine na statoru i rotoru prostoperiodi ne. Ugaona u estanost flukseva, struja i napona nastatoru i rotoru u αβ sistemu je ω s , a realnih struja i flukseva na rotoru ω ω ωs m k− = inaziva se ugaona u estanost klizanja. To to je stacionarno stanje prostoperiodi no ote avasintezu regulatora brzine, a to nam je jedan od ciljeva. Sinteza bi bila mnogo lak a ako bi ustacionarnom stanju sve veli ine bile konstantne.
Da bi se ovo ostvarilo, potrebno je usvojiti novi koordinatni sistem. Po to magnetnopolje rotira u estano u napona napajanja ω s , komponente fluksa e biti konstantne u onomkoordinatnom sistemu koji rotira istom tom u estano u. Analogno tome e i komponentenapona i struje i na statoru i na rotoru u tom koordinatnom sistemu , u stacionarnom stanju, bitikonstantne. Usvoji emo jedan takav koordinatni sistem i ozna i emo ga sa dq, kao na slici 2.Na slici 2. ugao θ r je ugao koji osa d zaklapa sa osom α, a ostale oznake su poznate od ranije.
10
Cilj nam je da model asinhronog motora prika emo u ovom koordinatnom sistemu. Postupak jeisti kao u slu aju prelaska sa α’β’ na αβ sistem :
u
u
u
uT
u
uds
qs
r r
r r
s
s
s
s
⎡
⎣⎢
⎤
⎦⎥ =
−⎡
⎣⎢
⎤
⎦⎥
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥
cos sin
sin cos
θ θθ θ
α
β
α
β1 ,
(2.24)
i
iT
i
ids
qs
s
s
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥1
α
β ,
(2.25)
Slika 2. Uvo enje dq koordinatnog sistema
ΨΨ
ΨΨ
ds
qs
s
s
T⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥1
α
β ,
(2.26)
( ) ( )( ) ( )
0
0 2
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ =
− −− − −
⎡
⎣⎢
⎤
⎦⎥
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥
u
u
u
uT
u
udr
qr
r m r m
r m r m
cos sin
sin cos'
'
'
'
θ θ θ θθ θ θ θ
α
β
α
β ,
(2.27)
i
iT
i
idr
qr
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥2
α
β
'
'
,
(2.28)
ΨΨ
ΨΨ
dr
qr
T⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥2
α
β
'
'
.
(2.29)Matrice T1 i T2 imaju istu osobinu kao i matrica T, tj. T TT
1 11= − i T TT
2 21= − .
Postupkom potpuno analognim malo as sprovedenom ( jedna ine (2.17) - (2.21) ) dobija se :
11
u
uR
i
id
dtds
qss
ds
qs
ds
qss
ds
qs
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥ +
−⎡
⎣⎢
⎤
⎦⎥
⎡
⎣⎢
⎤
⎦⎥
ΨΨ
ΨΨ
ω 0 1
1 0 ,
(2.30)
0
0
0 1
1 0
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ =
⎡
⎣⎢
⎤
⎦⎥ +
⎡
⎣⎢
⎤
⎦⎥ +
−⎡
⎣⎢
⎤
⎦⎥
⎡
⎣⎢
⎤
⎦⎥
u
uR
i
id
dtdr
qrr
dr
qr
dr
qrk
dr
qr
ΨΨ
ΨΨ
ω .
(2.31)
Jedna ine fluksnih obuhvata imaju isti oblik kao i u αβ sistemu :
Ψds s ds m drL i L i= + , Ψqs s qs m qrL i L i= + ,
Ψdr r dr m dsL i L i= + , Ψqr r qr m qsL i L i= + .
(2.32)
Sada emo izvesti i izraz za elektromagnetni moment. Pretpostavimo da smo primenili
Klarkovu transformaciju invarijantnu po snazi, tj. sa koeficijentom k = 2 3/ . Tada je :
[ ]P u i u i u i u i u i u ui
ia a b b c c s s s s s s
s
s
= + + = + =⎡
⎣⎢
⎤
⎦⎥ =α α β β α β
α
β
[ ]=⎡
⎣⎢
⎤
⎦⎥u u T T
i
ids qsT ds
qs1 1 ,
(2.33)
P u i u ids ds qs qs= + ,
(2.34)
P R id
dti R i
d
dtis ds
dss qs ds s qs
qs
s ds qs= + −⎛⎝⎜
⎞⎠⎟
+ + +⎛⎝⎜
⎞⎠⎟ =
ΨΨ
ΨΨω ω
( )= + + + −R i id
dti
d
dti is s ds
dsqs
qss ds qs qs ds
2 Ψ ΨΨ Ψω .
(2.35) Proanalizira emo ovaj izraz. lan R is s
2 predstavlja D ulove gubitke u namotajima
statora, lan id
dti
d
dtdsds
qsqsΨ Ψ
+ promenu energije akumulirane u magnetnom polju, pa prema
tome izraz ( )ω s ds qs qs dsi iΨ Ψ− mora predstavljati snagu koja se predaje rotoru. Ta snaga se
naziva snaga obrtnog polja i ozna ava se Pob :
( )P i iob s ds qs qs ds= −ω Ψ Ψ .
(2.36)
12
Elektromagnetni moment se dobija kad se snaga obrtnog polja podeli sa kru nomu estano u napajanja ω s :
MP
i iemob
sds qs qs ds= = −
ωΨ Ψ .
(2.37)
Ova formula va i za motore sa jednim parom polova. Ako postoji vi e pari polova imoment je ve i, i to proporcionalan broju pari polova :
( )M p i iem ds qs qs ds= −Ψ Ψ .
(2.38)
Elektromagnetni moment se mo e izraziti i preko komponenti fluksa rotora :
M pL
L
L
Li
L
Li
pL
L
L L
Li i L i i
L L
Li i L i i
emm
r
r
mds qs
r
mqs ds
m
r
r s
mds qs r dr qs
r s
mds qs r qr ds
= −⎛⎝⎜
⎞⎠⎟ =
= + − −⎛⎝⎜
⎞⎠⎟ =
Ψ Ψ
( ) ( )= + − − = −pL
LL i i L i i L i i L i i p
L
Li im
rr dr qs m ds qs m ds qs r qr ds
m
rdr qs qr dsΨ Ψ .
(2.39)Sada je kompletiran model asinhronog motora u dq koordinatnom sistemu. Mogu e je
izraziti elektromagnetni moment i preko veli ina u αβ sistemu. Na osnovu slike 2. imamo :
Ψ Ψ Ψds s r s r= +α βθ θcos sin , Ψ Ψ Ψqr r r r r= − +α βθ θsin cos ,
i i ids s r s r= +α βθ θcos sin , i i iqs s r s r= − +α βθ θsin cos .
(2.40)
Zamenom (2.40) u (2.39) dobija se :
( )M pL
Li iem
m
rr s r s= −Ψ Ψα β β α .
(2.41)
2.4. Relativizacija jedna ina
Relativizacija jedna ina je posledica relativizacije fizi kih veli ina, odnosno posmatranjafizi kih veli ina ne u svojim realnim vrednostima nego u odnosu na neke odre ene konstantnevrednosti tih fizi kih veli ina. Te odre ene konstantne vrednosti nazivaju se bazne i ozna avajuse indeksom B, kao na primer ΨB B BU I, , . Trenutna vrednost fizi ke veli ine predstavlja sebezdimenzionim brojem koji je jednak koli niku te trenutne vrednosti i odgovaraju e baznevrednosti.
13
Za svaki motor postoje ograni enja u vidu maksimalnog fluksa, napona i trajne vrednostistruje sa kojima mo e raditi. Fluks je ograni en zbog zasi enja feromagnetnog jezgra, napon jeograni en zbog proboja izolacije a struja zbog zagrevanja provodnika. Postoje i nominalnevrednosti za koje je motor predvi en da radi. One su manje ili jednake od najve ih dozvoljenihvrednosti.
Bazne vrednosti se mogu izabrati proizvoljno, mada je uobi ajeno da se za fluks, napon istruju nominalne vrednosti uzimaju za bazne. Bazne vrednosti za ostale veli ine koje sepojavljuju u analizi motora mogu se izraziti preko ve usvojenih baznih vrednosti za drugefizi ke veli ine.
Relativne vrednosti se ozna avaju dodavanjem oznake [p.u.] ( od engleskog per unit )iza oznake fizi ke veli ine. Zbog sa etijeg zapisa mi emo u daljem tekstu za relativne jedinicekoristiti isko ena ( italik ) slova, npr. Ψ, u, i .
Nave emo najva nije nominalne i usvojene bazne vrednosti za konkretan model motorakojim emo se baviti :
I I AB nom= = 2 1. , cos .ϕ n = 0 75 - nominalni faktor snage ,
(2.42)Relativizacija jedna ina se vr i tako to se svaka fizi ka veli ina zameni proizvodom te
fizi ke veli ine izra ene u relativnim jedinicama i njene bazne vrednosti, npr. i Ids ds B= i . Zatimse sve bazne vrednosti grupi u. One e se ili sve skratiti ostavljaju i jedna inu u potpuno istomobliku ili e ostaviti neki konstantan koeficijent. Evo dva primera za relativizaciju jedna ina :
Ψds s ds m drL i L i= + ,(2.43)
Ψ ds B s B ds B m B dr BL I L IΨ = +L i L i ,(2.44)
( )Ψ dsB B
Bs ds m dr s ds m dr
L I= + = +
ΨL i L i L i L i .
(2.45)
Kao to se vidi jedna ina je relativizacijom o uvala potpuno isti oblik. Evo jednogprimera gde to nije slu aj, tj. gde se pojavljuje dodatni koeficijent :
14
d
dtu R is
s s s
Ψαα α= − ,
(2.46)
ΨΒd
dtU R Is
s B s s B B
Ψ αα α= −u R i ,
(2.47)
( ) ( )d
dt
Us Bs s s B s s s
Ψ αα α α αω= − = −
ΨΒ
u R i u R i .
(2.48)
Kao to se vidi, pojavio se koeficijent ω B koji je posledica toga to za vreme, koje jetako e fizi ka veli ina, nismo usvojili baznu vrednost nego ga posmatramo u realnoj vrednosti.
2.5. Indirektna vektorska kontrola
Kada elimo da koristimo asinhroni motor u servomehanizmima, odnosno kada jepotrebna brza i kvalitetna regulacija polo aja i/ili brzine, neophodan uslov za to je veoma brzaregulacija elektromagnetnog momenta. Razlog za to je to u npr. pozicionom servomehanizmupostoje tri petlje regulacije ( videti sliku 3. ) : spolja nja reguli e ugao a upravlja ka veli ina jeugaona brzina, srednja kontura reguli e brzinu a upravlja ka veli ina je moment, a unutra njakontura reguli e moment. Poznato je iz teorije regulacije da unutra nja kontura mora bitinajbr a, a spolja nja najsporija.
θref Δθ REG. ω * Δω REG. M* REG. M 1 ω 1 θ
- POLO• AJA - BRZINE MOMENTA Js s
Slika 3. Principijelan izgled pozicionog servomehanizma
Ta razlika u brzinama, odnosno u irinama propusnih opsega mora biti bar tolika da oddve susedne konture regulacije dinamika one unutra nje prakti no ne uti e na dinamiku onespolja nje, odnosno da je unutra nja kontura za spolja nju transparentna. Time se posti erasprezanje pojedinih kontura regulacije i omogu ava se nezavisna sinteza regulatora.
Potreban uslov za brzu regulaciju momenta je brza regulacija struje, ali to nije dovoljanuslov, jer moment nije slika neke komponente struje, ve zavisi od nekoliko komponenti strujekao i od fluksa ( relacije (2.38), (2.39), (2.41) ). Brza regulacija struje se ostvaruje napajanjemmotora iz trofaznog tranzistorskog invertora kojim upravlja histerezisni strujni regulator. Naslici 4. prikazani su invertor i strujni regulator za jednu fazu statora.
15
Slika 4. Trofazni invertor sa histerezisnim strujnim regulatorom za jednu fazu statora
Postupak koji na osnovu eljenog momenta M∗ sa slike 3. generi e reference za strujei i ia b c
∗ ∗ ∗, , i koji obezbe uje da se brza regulacija struje prenese na moment naziva se vektorskakontrola, a mo e biti direktna ( bez senzora na osovini motora ) i indirektna ( sa senzorom naosovini ). Uloga senzora je da meri ugao osovine motora θ m . Mi emo se baviti indirektnomvektorskom kontrolom.
Indirektna vektorska kontrola ( IVK ) se bazira na modelu asinhronog motora u dqkoordinatnom sistemu, koji je dat jedna inama (2.4), (2.30) - (2.32) i (2.39). Po to statorskestruje reguli emo velikom brzinom, mo e se smatrati da je motor strujno napajan. Komponentestruje statora su tad upravlja ke promenljive, a ne promenljive stanja. Zato u razmatranje neuzimamo jedna ine naponskog balansa za stator, ve samo za rotor. Za promenljive stanjauze emo Ψdr i Ψqr :
0 = = + −u R id
dtdr r drdr
k qr
ΨΨω ,
(2.49)
0 = = + +u R id
dtqr r qrqr
k dr
ΨΨω ,
(2.50)
0 = = −⎛⎝⎜
⎞⎠⎟ + −∗u R
L
L
Li
d
dtdr rdr
r
m
rds
drk qr
Ψ ΨΨω ,
(2.51)
0 = = −⎛⎝⎜
⎞⎠⎟ + +∗u R
L
L
Li
d
dtqr rqr
r
m
rqs
qrk dr
Ψ ΨΨω ,
(2.52)
16
d
dt T
L
Tidr dr
rk qr
m
rds
Ψ ΨΨ= − + + ∗ω ,
(2.53)
d
dt T
L
Tiqr
k drqr
r
m
rqs
ΨΨ
Ψ= − − + ∗ω ,
(2.54)
pri emu je T L Rr r r= / vremenska konstanta rotora. Karakteristi ne u estanosti ovog sistema
su s T jr k1 2 1, /= − ± ω a to ne obezbe uje dovoljno brzu regulaciju momenta, koji zavisi od
fluksa.Jedna ina za elektromagnetni moment glasi :
( )M pL
Li iem
m
rdr qs qr ds= −Ψ Ψ .
(2.55)
Ako bismo postigli da bude Ψqr = 0 , tada bi moment bio proporcionalan proizvodu
Ψdr qsi , a po to je u toku rada motora fluks naj e e konstantan, moment bi prakti no bio slika
struje iqs i bila bi obezbe ena dovoljno brza regulacija momenta. Ako stavimo Ψqr = 0 u
jedna ine (2.53) i (2.54) dobijamo :
d
dt T
L
Tidr dr
r
m
rds
Ψ Ψ+ = ∗ ,
(2.56)
ω ωk drm
rqs k
m qs
dr r
L
Ti
L i
TΨ
Ψ= ⇒ =∗
∗
.
(2.57)
Ako je ugaona u estanost klizanja odre ena relacijom (2.57) tada iz (2.54) sledi :
d
dt Tqr qr
r
Ψ Ψ+ = 0 ,
(2.58)
to zna i da i ako postoji po etna vrednost Ψqr ≠ 0 , ona vremenom nestaje i bi e Ψqr = 0 . U
tom slu aju pona anje motora mo e se predstaviti emom sa slike 5. Vidimo da se motor pona akao motor jednosmerne struje sa nezavisnom pobudom : struja ids odre uje fluks, a struja iqs
odre uje moment. Fluks je uglavnom konstantan i jednak nominalnom fluksu, a do potrebe zanjegovim menjanjem dolazi pri radu na brzinama ve im od nominalne. Motor generi eelektromotornu silu proporcionalnu fluksu i ugaonoj brzini, a napon napajanja, koji je ograni en,mora biti ve i od te elektromotorne sile da bi se mogla injektovati struja statora. Zbog toga je
17
potrebno smanjivati fluks na brzinama ve im od nominalne i to po zakonu 1/ω. To je rad uoblasti slabljenja polja.
ids∗
Lm Ψdr
1+ sTr
iqs∗
pL m M em
Lr
Slika 5. Model motora kada je zadovoljena relacija (2.57)
Funkcija bloka koji realizuje IVK je slede a : on na osnovu eljenog fluksa odre ujestruju ids
∗ , a na osnovu eljenog momenta i znaju i kako se menja fluks odre uje se struja iqs∗ .
Ovaj blok tako e obezbe uje da bude zadovoljena relacija (2.57). Senzor na osovini motorameri ugao θ m , a na osnovu njega i ω k se iz relacije
θ θ ωr m k
t
dt= + ∫0
(2.59)
dobija ugao θ r koji je potreban da bi se na osnovu ids∗ i iqs
∗ dobile struje i sα∗ i i sβ
∗ . Mo emo
smatrati da strujni regulator mo e da reguli e te struje zbog njihovih jednostavnih aritmeti kihveza sa realnim strujama i a
∗ , i b∗ , i c
∗ . Na osnovu dosada njih razmatranja mo e se nacrtatiema bloka IVK :
Parkova tr.Ψr
∗ 1 ids
∗ dq i sα
∗
Lm αβ i sβ∗
1 θ r
1+sTr
Mem∗
Lr iqs
∗
Lm ωk ∫ θm
pLm Lr
Slika 6. ema bloka za indirektnu vektorsku kontrolu
2.6. Korisni ki model sprege IVK+AM
U ovom odeljku bi e prikazano kako korisnik “vidi” spregu bloka IVK i asinhronogmotora. Ta sprega naravno uklju uje i trofazni invertor sa histerezisnim strujnim regulatorom.Motor se modelira emom sa slike 5. jer radi u sprezi sa IVK, a IVK blok je predstavljen emomsa slike 6. Blokovi bi se spojili tako to bi se struje ids
∗ i iqs∗ iz bloka IVK vodile na ulaz eme sa
18
slike 5. Rezultuju i model, koji podrazumeva da je regulacija struje idealna, je veomajednostavan i prikazan je na slici 7.
Po to je na krajnji cilj simulacija rada motora u brzinskom servomehanizmu, potrebnoje da utvrdimo kakav regulator emo koristiti za regulaciju brzine. U tu svrhu e nam poslu itimodel sa slike 7. Ako bismo se odlu ili za regulator sa isto proporcionalnim dejstvom, tada bizbog toga to postojanje nekog konstantnog poreme aja u vidu momenta optere enja M opt
zahteva postojanje elektromagnetnog momenta iste vrednosti u stacionarnom stanju i na osnovuM kem p
∗ = Δω postojala i gre ka u brzini u stacionarnom stanju, to se ne mo e dozvoliti.
Mopt
Mem∗
+ 1 ω m 1 θ m
Js s
Ψr∗ 1 Ψdr
1+sTr
Slika 7. Korisni ki model sprege IVK+AM
Kori enje regulatora sa samo integralnim dejstvom bi prouzrokovalo nestabilnostsistema, pa emo koristiti PI regulator. Fluks i moment koji se mogu dovesti na ulaz modela saslike 7. su ograni eni. Fluks zavisi od ugaone brzine na ranije obja njen na in, a moment koji semo e razviti tako e zavisi od ugaone brzine i parametara motora. Korisno je utvrditi koliki jemaksimalni moment koji motor mo e da razvije na nekoj ugaonoj brzini i inkorporirati toograni enje u regulator, koji bi u suprotnom mogao da generi e eljeni moment M em
∗ koji jeve i od ostvarljivog momenta to bi dovelo do navijanja regulatora i zasi enja aktuatora, uovom slu aju invertora. Da bismo do tog rezultata do li u to jednostavnijem oblikupretpostavi emo da je blok IVK postigao Ψqr = 0 , da je otpor R s u jedna inama naponskog
balansa za stator (2.30) zanemarljivo mali i da je u dq koordinatnom sistemu postignutostacionarno stanje. Jedna ine naponskog balansa za stator se svode na :
u ds s qs= −ω Ψ , u qs s ds= ω Ψ ,
(2.60)
( )u L i L ids s s ds m dr= − +ω , ( )u L i L iqs s s qs m qr= +ω .
(2.61)
Na osnovu Ψqr = 0 , jedna ine (2.49) i pretpostavljenog stacionarnog stanja dobija se
i dr = 0 , a na osnovu izraza za Ψqr u (2.32) se dobija da je
iL
Liqr
m
rqs= − .
(2.62)
19
Sada nam jedna ine (2.61) daju :
u L ids s s ds= −ω , u LL
Li L iqs s s
m
rqs s qse
= −⎛⎝⎜
⎞⎠⎟ =ω ω γ
2
,
(2.63)
pri emu se veli ina Leγ naziva ekvivalentnom rasipnom induktivno u motora. Odavde se
dobija izraz za napon statora u dq sistemu :
( ) ( )u L i L idq s s ds qse= +ω γ
2 2 .
(2.64)
Ovu jedna inu emo relativizovati :
( ) ( )u L i L idq B s B B B s ds qsU L Ie
= +ω ω γ2 2
,
(2.65)
( ) ( )u L i L idq s s ds qse= +ω 2 2
γ .
(2.66)
Po to je U UB nom= iUnom je maksimalna efektivna vrednost napona jedne faze statora,
sledi da je u u umax max maxa b c= = = 2 . Odatle sledi da je u /
maxs = 3 2 , a zatim se na osnovu toga
to je primenjena Klarkova transformacija invarijantna po snazi dobija da je umaxdq = 3 . Sli no
se i za fluks dobija da je Ψ drmax= 3 . Na osnovu u
maxdq = 3 i (2.66) dobija se :
( )iL
L iL
L
Lmaxqss
s dss
s
mdr
e e
= − = −⎛⎝⎜
⎞⎠⎟
1 3 1 32
2
2
2
γ γω ωΨ ,
(2.67)
a elektromagnetni moment dat je izrazom
M pL
Liem
m
rdr qs= Ψ .
(2.68)
i vidimo da je moment slika struje iqs jer je fluks Ψdr konstantan. Fluks Ψ dr = 3 na ugaonim
u estanostima manjim od neke na kojoj emo zapo eti slabljenje polja. Iz relacije (2.67) se vidida se slabljenje polja mora zapo eti na u estanosti manjoj od ω nom da bi se uop te mogla
injektovati struja iqs , odnosno da bi potkorena veli ina bila uvek ve a od 0 . Za potrebe na eg
modela u SIMULINK-u emo uzeti da je ω sgr= 0 9. i na ve im u estanostima emo smanjivati
fluks obrnuto proporcionalno sa ω s .
20
Iz relacije (2.67) se vidi da se za male vrednosti ω s dobijaju veoma velike vrednosti zastruju, to nije realno jer smo zanemarili omsku otpornost namotaja statora. Tako e, du aprimena tolike struje koja je mnogostruko ve a od nominalne bi mogla da prouzrokujepregorevanje namotaja kod realnog motora. Zbog ovih razloga mi emo za struju iqs koristiti
slede e ograni enje : za ω ωs sgr≤ emo uzeti i .
maxqs const= i jednako onoj vrednosti koja se iz
formule (2.67) dobija za ω ωs sgr= . Za ω ωs sgr
> koristi emo formulu (2.67) . Kona an oblik
promene maksimalne dozvoljene struje i momenta se dobija kada se i zakon promene fluksauzme u obzir :
iL
L
L L.
L
Lmaxqs
s
s
m
s
me gr e
= −⎛⎝⎜
⎞⎠⎟ = −
⎛⎝⎜
⎞⎠⎟
3 1 312346
2
2 2
γ γω za ω ωs sgr
≤ = 0 9. , (2.69)
iL
L
L
.
L.
L
Lmaxqss
s
m s s
s
me e
= −⎛⎝⎜
⎞⎠⎟ = −
⎛⎝⎜
⎞⎠⎟
1 33
0 9 31 0 9
2
2 2
γ γω ω ω za ω ωs sgr
> . (2.70)
Da bismo dobili ograni enje momenta prvo emo relativizovati jedna inu (2.68) :
ML
Liem B
m
rdr qs B BM p I= Ψ Ψ ,
(2.71)
ML
Li
L
Liem
m
rdr qs
B B
B
m
rdr qs
B B B
B B n n
pI
Mp
I
U I cos= =Ψ Ψ
Ψ Ψ ωϕ η3
,
(2.72)
ML
Liem
n n
m
rdr qs
p
cos=
3 ϕ ηΨ .
(2.73)
Sada za maksimalni moment, uz pomo relacija (2.69) i (2.70), dobijamo :
ML
Li
L
L L.
L
L
max max maxemn n
m
rdr qs
n n
m
r
s
m
p
cos
p
cose
= =
= −⎛⎝⎜
⎞⎠⎟ =
3
33
312346
2
ϕ η
ϕ η γ
Ψ
= −⎛⎝⎜
⎞⎠⎟ =p
cosconst
n n
m
r
s
meϕ η γ
L
L L.
L
L.
112346
2
za ω ωs sgr≤ ,
(2.74)
21
ML
Li
L
L
.
L.
L
L
max max maxemn n
m
rdr qs
n n
m
r s s
s
m
p
cos
p
cose
= =
= −⎛⎝⎜
⎞⎠⎟ =
3
33
0 9 31 0 9
2
ϕ η
ϕ η γ
Ψ
ω ω
= −⎛⎝⎜
⎞⎠⎟
p
cos n n
m
r s
s
m seϕ η γ
L
L
.
L.
L
L~
0 9 11 0 9
12
2
2ω ω za ω ωs sgr
> .
(2.75)
Za konkretan motor, osim podataka datih u (2.42), poznati su i slede i parametri :
p = 2 , L Lr s= = 2 , L .m = 1916 , R .s = 0 038726 , R .r = 0 04364 ,(2.76)
pa se ubacivanjem ovih vrednosti u formule (2.69), (2.70), (2.74), (2.75) dobija :
i.
..maxqs
s
ss
=≤
>
⎧⎨⎪
⎩⎪
4 0 936
0 9
za
za
ω
ωω , M
.
.maxem
s
s
s=
≤
>
⎧⎨⎪
⎩⎪
8 0 9
0 92
.1 za 6.56
za
ω
ωω .
(2.77)
2.7. Odre ivanje parametara digitalnog PI regulatora brzine
Za regulator smo usvojili PI regulator i njegova funkcija prenosa u z-domenu je
( ) ( )( )
( )( ) ( )W z
M z
z
M z
z zK
K
zR
em em
ref m
PI= =
−= +
−
∗ ∗
−Δω ω ω$ 1 1 ,
(2.78)
gde su ( )ω ref z i ( )$ω m z respektivno z-transformacije referentne ugaone brzine i izmerene
ugaone brzine a KP i KI su respektivno koeficijenti proporcionalnog i integralnog dejstva.
Pretpostavi emo da je regulacija momenta idealna, tj. da je ( ) ( )M z M zem em∗ = , to zna i da je
elektromagnetni moment koji motor razvija konstantan u toku jedne periode odabiranja. Naosnovu toga se i ugaona brzina motora u toku jedne periode odabiranja menja linearno, ako jeopteretni moment konstantan. Ako uvedemo za periodu odabiranja oznaku T, mo emo pisatislede e relacije :
( ) ( ) ( )ω ωm m emk kT
JM k= − + −1 1 ,
(2.79)
22
( ) ( ) ( ) ( )θ θ
ω ωm m
m mk k Tk k
= − ++ −
11
2 .
(2.80)
Inkrementalni enkoder meri ugao θ m a prate a elektronika odre uje ugaonu brzinu :
( ) ( ) ( )$ωθ θ
mm mk
k k
T=
− − 1 .
(2.81)
Prelaskom u z-domen dobija se :
( ) ( ) ( ) ( ) ( )ω ω ωm m em m emz z zT
Jz M z z
T
J zM z= + ⇒ =
−− −1 1 1
1 ,
(2.82)
( ) ( ) ( ) ( ) ( ) ( )θ θ θω ω
ωm mm m
m mz z z Tz z z
zT z
zz= +
+⇒ = +
−−
−1
1
2 2
1
1 ,
(2.83)
( ) ( ) ( ) ( )$ωθ θ
θmm m
mzz z z
T
z
Tz=
−= −− −1 11
.
(2.84)
Na osnovu ovih relacija mo emo dati model celog sistema regulacije u z-domenu :
Mopt
ωref + Δω KP (z-1)+KI z Mem T 1 ωm T z+1 θm z-1 z-1 + J z-1 2 z-1 Tz
$ω m
Slika 8. Model sistema u z-domenu
Funkcija povratnog prenosa je
( ) ( ) ( )( )( )
W zK z K z
z
T
J z
T z
z
z
Tz
T
J
z K z K K z
z zP
P I P P I=− +
− −+−
− =+ − +
−
1
1
1
1 2
1
1
1
2
1
12
. (2.85)
Da bismo odredili koeficijente KP i KI potrebno je usvojiti kriterijum kojim emooceniti kvalitet prelaznog re ima ovog sistema. Za kriterijum emo usvojiti da sistem sazatvorenom povratnom spregom nema kompleksnih polova. Po to je sistem tre eg reda, tajkriterijum se analiti ki mo e izraziti kao
23
( ) ( )f z z= − < <σ σ30 1 , ,
(2.86)
gde je f(z) karakteristi ni polinom sistema. Na osnovu (2.85) karakteristi ni polinom sistema,ije su nule jednake nulama funkcije ( )1+ W zP , dat je izrazom
( ) ( ) ( )( )f z z zT
Jz K z K K zP P I= − + + − + =1
21
2
= + + −⎛⎝⎜
⎞⎠⎟
+ +⎛⎝⎜
⎞⎠⎟
−z zK T
J
K T
Jz
K T
J
K T
JP I I P3 2
2 22 1
2 2 .
(2.87)
Ekvivalentiranjem desnih strana jedna ina (2.86) i (2.87) dobija se sistem tri jedna ine sa trinepoznate ije je re enje :
σ = − =4 1 0 58743 . , K T
JP
20 2= . ,
K T
JI
20 035= . ,
(2.88)
i prema tome, sa poznatim J se u zavisnosti od T odre uju KP i KI . Me utim, na i modeli uSIMULINK-u e biti u relativnim jedinicama, a to va i i za regulator. Zna i da treba odreditikoeficijente k P i k I koji figuri u u relativizovanoj jedna ini :
( ) ( )M em PIz k
k z
zz= +
−⎛⎝⎜
⎞⎠⎟1
Δω .
(2.89)
Ove koeficijente emo odrediti relativizovanjem jedna ine (2.78) :
( ) ( )M em B PI
Bz M KK z
zz= +
−⎛⎝⎜
⎞⎠⎟1
Δω ω ,
(2.90)
( ) ( )M emP B
B
I B
B
zK
M
K
M
z
zz= +
−⎛⎝⎜
⎞⎠⎟
ω ω1
Δω .
(2.91)
Upore uju i (2.89) i (2.91) dobijamo :
24
kK
MPP B
B
=ω
, kK
MII B
B
=ω
.
(2.92)Sada na osnovu (2.88) imamo :
kT
J
MPB
B
= 0 4. ω , k
T
J
MIB
B
= 0 07. ω .
(2.93)
Koeficijent J
MB
B
ω ima dimenziju vremena a javlja se i kod relativizovanja mehani ke
jedna ine sistema :
( )Jd
dtM
J
M
d
dtBm
B em optB
B
mem optω ωω ω
= − ⇒ = −M M M M .
(2.94)
Radi kra eg zapisa ozna i emo taj koeficijent sa Tmeh i zva emo ga mehani ka vremenskakonstanta, pa imamo
Td
dtmehm
em opt
ω= −M M .
(2.95)
Tako e iz jedna ina (2.93) imamo :
kT
TPmeh= 0 4. , k
T
TImeh= 0 07. .
(2.96)
Kod konkretnog motora kojim emo se baviti je T smeh = 1 .
3. MODELIRANJE U PROGRAMSKOM PAKETU SIMULINK
U ovom poglavlju bi e prikazani modeli u SIMULINK-u svih komponenti koje su bilerazmatrane u prethodnom poglavlju. Bi e dat model asinhronog motora, model trofaznogtranzistorskog invertora sa strujnim regulatorom, model bloka za IVK sa modelominkrementalnog enkodera, zatim model brzinskog regulatora i kola koje meri brzinu motora.
25
Zatim e biti prikazan model kompletnog elektromotornog pogona sa brzinskom regulacijomkoji sadr i sve navedene komponente. Na kraju e biti dat primer re avanja konkretnogin enjerskog problema uz pomo ovih modela.
Pravljenje modela u SIMULINK-u je jednostavno i odvija se grafi kim putem. Postojimno tvo blokova koji su podeljeni u biblioteke : Discrete, Linear, Nonlinear, Sources, Sinks itd.Svaki blok ima svoj dialog box, gde se nalazi kra e obja njenje funkcije bloka, a ve ina blokovaima parametre koji se ovde unose. Npr. za blokove SUM i PRODUCT ( sabira i mno a )jedini parametar je broj ulaza, za integrator to je po etna vrednost itd.
Dve zna ajne osobine SIMULINK-a su grupisanje ( Group ) i maskiranje ( Mask )blokova. Grupisanje je vizuelna zamena grupe blokova jednim blokom koji je tipa Subsystem.Time se posti e vizuelno pojednostavljenje komplikovanih modela. U jednom modelu mo e bitivi e Subsystem blokova, a grupisanje mo e biti i na vi e nivoa, tj. da jedan Subsystem blok usebi sadr i druge Subsystem blokove. Maskiranje ide korak dalje i omogu ava kreiranje dialogbox-a za taj jedan blok koji zamenjuje grupu blokova. U tom dialog box-u se unose sviparametri koji su potrebni za grupisane blokove. Mogu e je kreirati i izgled novog bloka.
Svi modeli i sve veli ine e biti u relativnim jedinicama, osim ugla θ m koji e biti urealnoj vrednosti zbog ve e ta nosti u ra unu.
3.1. Model asinhronog motora
Model asinhronog motora je u αβ sistemu, u relativnim jedinicama i dat je na slici 9.Dobija se iz analiti kog modela u αβ sistemu, koji je dat relacijama (2.21) - (2.23), (2.41) i(2.95). Jedna ina (2.95) je ve u relativnim jedinicama, a i ostale jedna ine treba relativizovati.Definicione jedna ine flukseva se relativizacijom ne menjaju ( videti (2.45) ), a ovde emonavesti relativizovane jedna ine (2.21) i (2.22) u eksplicitnom obliku ( izvod fluksa izdvojen nalevoj strani jedna ine - pogledati (2.48) ), kao i relativizovanu jedna inu (2.41) :
( )d
dts
B s s s
Ψ αα αω= −u R i , ( )d
dts
B s s s
Ψ ββ βω= −u R i ,
( )d
dtr
B r r m r
ΨΨωα
α βω= − −R i , ( )d
dtr
B r r m r
ΨΨωβ
β αω= − +R i ,
(3.1)
( )ML
Li iem
n n
m
rr s r s
p
cos= −
3 ϕ η α β β αΨ Ψ .
(3.2)
Prve dve od relacija (3.1) direktno su unesene u model ( blokovi G4, G5, Sum, Sum7,G11, G12 ), a druge dve emo modifikovati tako to emo izraziti struje rotora prekokomponenti fluksa rotora i struja statora, na osnovu definicija komponenti fluksa rotora, pa sedobijaju slede i izrazi :
d
dtr
Br
rr
r m
rs m r
ΨΨ Ψωα
α α βω= − + −⎛⎝⎜
⎞⎠⎟
R
L
R L
Li ,
26
d
dtr
Br
rr
r m
rs m r
ΨΨ Ψωβ
β β αω= − + +⎛⎝⎜
⎞⎠⎟
R
L
R L
Li .
(3.3)
Ove dve relacije su direktno inkorporirane u model ( blokovi Sub, Prod8, Prod9, Sum8, Sum9,G15 i G16 ). Treba jo objasniti kako se ra unaju komponente struje statora. One se izra avajupreko komponenti fluksa rotora i statora eliminisanjem komponenti struje rotora iz definicijaflukseva. Dobijaju se relacije :
iL
L
Lα
γα αs s
m
rr
e
= −⎛⎝⎜
⎞⎠⎟
1 Ψ Ψ ,
iL
L
Lβγ
β βs sm
rr
e
= −⎛⎝⎜
⎞⎠⎟
1 Ψ Ψ .
(3.4)
Dve poslednje relacije tako e su direktno unesene u model, a u njihovom predstavljanjuu estvuju blokovi G, G3, Sum1, Sum2, G1 i G2.
Za odre ivanje elektromagnetnog momenta koristi se relacija (3.2), koja je u modelupredstavljena blokovima G, G3 i Sub1.
Ulazi modela su uαs , uβs i M opt , a izlazi Ψ αs , Ψ βs , iαs , iβs , ω m , θ m i M em . Bitno
je skrenuti pa nju da su sve ulazne i izlazne veli ine u relativnim jedinicama, osim ugla θ m , kojije u stvarnoj vrednosti.
teta_mM_opt
+-
Sum7
+-
Sum
-K-
G12
-K-
G11
-K-
G13
-K-
G3
1/s
Int1
*
Prod9
*
Prod8
1/s
Int4
Sub-K-
G16
-K-
G15
+-+
Sum9
--+
Sum8
rs
G4
rs
G5
Sub1
+-
Sum2
-K-
G1
-K-
G2
-K-
G
+-
Sum1
1/s
Int3
1/s
Int2
1/s
Int
f_alfa_r
f_beta_r
i_alfa_s
i_beta_su_beta_s
u_alfa_s
omega_m
M_em
f_alfa_s
f_beta_s
Slika 9. Model asinhronog motora u SIMULINK-u
27
Va no je napomenuti da izlaz ω m iz modela motora treba podeliti sa p da bi se dobilastvarna brzina obrtanja, po to ve i broj pari polova smanjuje brzinu obrtanja pri istoj u estanostinapajanja.
Po etna stanja svih integratora su 0, a poja anja Gain blokova za one za koje se ne vididirektno sa slike su : za G i G3 poja anje je L / Lm r , za G1 i G2 je 1 / Lγ e
, a za blokove G11-
G13 i G15-G16 poja anje iznosi ω nom .U modelu motora postoje i dva Subsystem bloka iji izgled je dat na slikama 10. i 11.
1/Tmeh
G9
1/s
Int5
+-
Sum5
km
G8
+-
Sum6
1
out_1
*
Prod6
*
Prod5
4
in_4
M_em
5
in_5
2
out_2M_opt
3
in_3
2
in_2
1
in_1
Slika 10. Podsistem Sub1 koji ra una M em i ω m
U podsistemu sa slike 10. vrednost konstante k M je data izrazom ( )p cos n n/ 3 ϕ η a
uvedena je zbog kra eg zapisa i zbog toga da bi model motora imao manji broj parametara.Blokovi ozna eni sa IN su ulazi, a sa OUT izlazi podsistema.
rr/lr
G14
2
out_2
rr*lm/lr
G6
1
out_1
2
in_2
1
in_1
3
in_3
4out_4
rr*lm/lr
G7
rr/lr
G13
3
out_3
4
in_4
f_beta_r f_alfa_r
i_alfa_si_beta_s
Slika 11. Podsistem Sub
Primenom postupaka grupisanja i potom maskiranja ceo model asinhronog motora sezamenjuje jednim blokom, kojeg emo nazvati ASINHRONI MOTOR, i u ijem se dialog box-uunose svi potrebni parametri. Na slici 12. dat je izgled tog bloka.
ASINHRONI
MOTOR
ASMOT
f_alfa_s
f_beta_su_beta_s
u_alfa_s
M_opt
i_alfa_s
i_beta_s
teta_m
omega_m
M_em
Slika 12. Maskirani model asinhronog motora
28
Uz pomo ovog modela emo prikazati kako se pona a asinhroni motor priklju en namre ni napon, tj. na napon konstantne u estanosti i amplitude, i sa opteretnim momentom kojije skokovita funkcija u vremenu. Simulira se priklju enje motora u mre u u trenutku t=0. Naslici 13. je prikazan ovaj eksperiment u SIMULINK-u.
100*pi
Constant1ASINHRONI
MOTOR
ASMOT
*
Product
*
Product1sin(u)
Fcn1
sqrt(3)
Constant
cos(u)
Fcn
1/s
Integrator
omega_s teta_s
M_opt
u_alfa_s
u_beta_s
M_em
omega_m
teta_m
i_beta_s
f_beta_s
i_alfa_s
f_alfa_s
Step Input
Step Input1
+
-
Sum
Slika 13. Simulacija priklju enja asinhronog motora na mre u
Da bismo dobili grafike promene izlaznih veli ina u modelu motora potrebno je na izlazemotora povezati neke od blokova iz Sinks biblioteke, a to su razni Scope blokovi. Oni su naslici 13. izostavljeni jer nemaju uticaj na sam tok eksperimenta. Na slici 14. prikazani su graficipromene nekoliko zna ajnih veli ina. Trajanje eksperimenta je 3 sec.
0 1 2 3-5
0
5
10
Elektromagnetni moment
t
M_em
0 1 2 30
0.2
0.4
0.6
0.8
1
Ugaona brzina
t
omega_m
0 1 2 30
2
4
6
8
10
Struja statora
t
i_rms
0 1 2 30
1
2
3Fluks statora
t
f_rms
Slika 14. Grafici nekih veli ina u eksperimentu sa slike 13.
Promena opteretnog momenta, koji je poreme aj u ovom eksperimentu, je data na slici15. Opteretni moment naglo menja vrednost sa 0 na 2M nom , a zatim opet pada na 0. Ovde je
29
bitno napomenuti da je, na osnovu relacije za M nom u (2.42) i toga to je p=2 nominalni
moment sa kojim motor radi zapravo jednak 2M nom .
0 0.5 1 1.5 2 2.5 3-1
-0.5
0
0.5
1
1.5
2
2.5
3Opteretni moment
t
M_opt
Slika 15. Opteretni moment u eksperimentu sa slike 13.
Iz ovog eksperimenta se vidi, a to smo i ranije napomenuli, da ugaona brzina motora prikonstantnoj u estanosti i amplitudi napona napajanja zavisi od momenta optere enja.
3.2. Model trofaznog invertora sa histerezisnim strujnim regulatorom
Princip rada sprege trofaznog invertora sa histerezisnim strujnim regulatorom je sasvimjednostavan. Upore uju se zadate vrednosti sa izmerenim vrednostima komponenti strujestatora, pa se u zavisnosti od toga koja je vrednost ve a i da li je ta razlika ve a od nekog pragana odre enu fazu statora dovodi pozitivan ili negativan konstantan napon. U praksi su merljivestruje i i ia b c, , pa je potrebno regulatoru dovesti zadate vrednosti struje i i ia b c
∗ ∗ ∗, , . Me utim, na
model motora je u αβ sistemu, a tako e i blok IVK generi e eljene vrednosti struja i is sα β∗ ∗, .
Zato i model na slici 16. ima ulazne i izlazne veli ine u αβ sistemu, a to predstavlja vrlo maloodstupanje od prakse zbog jednostavnih aritmeti kih veza izme u ova dva sistema struja. Ipak,po to su struje i i ia b c, , te koje elimo da odr avamo u odre enim granicama, unutar modela seprelazi na a, b, c sistem, odre uje se na koju fazu statora se dovodi kakav napon, a zatim se vr ipovratak na αβ sistem. Zna i, primenjuje se Klarkova transformacija u oba smera, kori enjem
koeficijenta k = 2 3/ . Transformacija u jednom smeru data je relacijama (2.5) ,
30
i_alfa_s*
i_beta_s*
i_beta_s
i_alfa_s
delta_i_b
delta_i_c
Relay
Relay2
1/2
Gain4
u_c1/2
Gain5
u_b
sqrt(2/3)
Gain1
delta_i_a u_a
Relay1
-+
Sum2
1/sqrt(6)
Gain2
u_alfa_s
+
-
Sum5
+--Sum6
sqrt(2/3)
Gain6
sqrt(2)
Gain9
u_beta_s+-
Sum1
+-
Sum
1/sqrt(2)
Gain3--
Sum3
Slika 16. Model trofaznog invertora sa histerezisnim strujnim regulatoromu drugom smeru relacije se dobijaju na osnovu (2.5) i V V Va b c+ + = 0 , pri emu V mo e bitioznaka i struje i napona :
i ia = 2
3α , i i ib = − +
1
6
1
2α β , i i ic = − −
1
6
1
2α β .
(3.5)
Jedan va an parametar koji se ne vidi direktno na emi sa slike 16. je korak histerezisa h,pod kojim se podrazumeva polovina irine histerezisa. Sva tri Relay bloka imaju isti korak. Akoje na primer h=0.05 to zna i da e svaka od struja i i ia b c, , mo i da odstupa od zadatihvrednosti najvi e za 5% od nominalne vrednosti struje, koja je ujedno i bazna vrednost struje.Grupisanjem blokova sa slike 16. i maskiranjem tako dobijenog Subsystem bloka dobija se jedanblok koji predstavlja invertor sa regulatorom. On je prikazan na slici 17. a u njegovom dialogbox-u mo e se definisati vrednost parametra h.
INVERTOR
SA STRUJNIM
REGULATOROM
INVERTOR
u_alfa_s
u_beta_s
i_alfa_s*
i_alfa_s
i_beta_s*
i_beta_s
Slika 17. Blok kojim se zamenjuje ema sa slike 16.
3.3. Model indirektnog vektorskog kontrolera
31
lm*rr/lr
Gain2
om_nomGain5
1/sIntegrator
+ +Sum
teta_r^
sin(u) Fcn1cos(u)Fcn2
i_ds*
+
+
Sum1
+
-
Sum2
*
*
*
*
M_em*
f_r*1/lm
Gain
+ -Sum3
lr/(lm*km)
Gain1
f_r
omega_k
teta_m^
*
i_alfa_s*
i_beta_s*
lm Gain3
*i_qs*
s1
Integrator1
rr*om_nom/lrGain4
1/uFcn
Slika 18. Model indirektnog vektorskog kontrolera
Model bloka za IVK u SIMULINK-u se jednostavno dobija iz eme bloka IVK sa slike6. Razlika je u tome to je model sa slike 18. u relativnim jedinicama, pa se pojavljuju i nekekonstante, odnosno poja anja. Tako e je blok sa funkcijom prenosa I reda u s-domenu iz emesa slike 6. predstavljen u modelu sa vi e blokova : Gain3, Sum3, Integrator1, Gain4. Razlogtome je to se na ovaj na in postavljanjem po etne vrednosti integratora mo e postaviti po etnavrednost za fluks f_r0. Ako se postavi da je po etna vrednost izlaza integratora npr.
( ) ( )J r r nom0 0 2= . L / R ω , bi e f_r0=0.2 . U stvarnosti fluks rotora u po etku ima vrednost 0.
Me utim, po to je relacija (2.77) koja defini e maksimalni ostvarljivi moment i koja se koristi ubrzinskom regulatoru izvedena na osnovu pretpostavljenog stacionarnog stanja, ona ne uzima uobzir injenicu da je fluks u po etku 0. Ako bismo postavili za po etnu vrednost fluksa umodelu neku malu vrednost, i regulator postavi maksimalnu vrednost momenta M em
∗ , to bi
prouzrokovalo veliku vrednost za struju iqs∗ , koja se realno ne mo e dobiti. Zato je bolje uzeti
f_r0= 3 , to e ograni iti struju iqs∗ . U kasnijem radu modela ovo nema zna aja. Na slici 19.
prikazan je maskiran model bloka IVK, a jedan od parametara u njegovom dialog box-u je if_r0.
INDIREKTNA
VEKTORSKA
KONTROLA
IVK
i_alfa_s*
i_beta_s*
f_r*
M_em*
teta_m^
Slika 19. Maskirani model bloka IVK
32
3.3.1. Model inkrementalnog enkodera ( IE )
IE meri ugao koji zauzima rotor motora u odre enim kvantima, tj. pokazivanje IE semenja za diskretne vrednosti ugla koji zauzima rotor, a i sam rezultat je diskretnog tipa. Takavna in merenja se jednostavno predstavlja jednim blokom u SIMULINK-u, a to je Quantizerblok. Njegov parametar je kvant, odnosno rezolucija merenja. Ulaz je stvarna pozicija rotora, aizlaz je izmerena, zaokru ena vrednost. Na slici 20. su prikazani model IE koji se sastoji odjednog bloka i blok koji se dobija maskiranjem tog modela. Parametar kod maskiranog bloka jebroj proreza N na IE, a korak kvantizacije, tj. rezolucija merenja se ra una kao 2π / N . Ukonkretnom slu aju je N=1024.
Quantizer
teta_m^teta_mINKREMENTALNI
ENKODER
IE
Slika 20. Model IE i blok koji se dobija njegovim maskiranjem
Na jednom primeru emo pokazati kako rade u sprezi IVK, trofazni invertor sa strujnimregulatorom i asinhroni motor. Na ulaz M em
∗ bloka IVK emo dovesti povorku impulsa , a naulaz f_r* konstantnu vrednost. Za moment inercije motora emo postaviti neku veoma velikuvrednost, teorijski beskona nu, da se motor ne bi obrtao. To je zbog toga to najve i momentkoji se mo e ostvariti opada sa brzinom obrtanja, a na ovaj na in uvek je mogu e ostvaritizadati moment. Posmatra emo kako elektromagnetni moment koji motor razvija prati zadatimoment. ema u SIMULINK-u data je na slici 21. Svi parametri motora su zadr ali vrednostiod ranije, samo moment inercije, odnosno mehani ka vremenska konstanta Tmeh ima jako velikuvrednost. Moment optere enja je 0, odnosno motor radi u praznom hodu. Zadati moment menjavrednost izme u nule i nominalnog momenta. Eksperiment traje 4 sec.
INDIREKTNA
VEKTORSKA
KONTROLA
IVK
PulseGenerator
sqrt(3)
Constant1
M_em*
f_r*
INKREMENTALNI
ENKODER
IE
ASINHRONI
MOTOR
ASMOT
teta_m
i_beta_s
i_alfa_su_alfa_s
u_beta_s
M_opt
i_alfa_s*
i_beta_s*
0
Constant M_em
INVERTOR
SA STRUJNIM
REGULATOROM
INVERTOR
33
Slika 21. Eksperiment sa zadatim momentom u vidu povorke impulsa
Na slede im slikama prikazani su grafici promene svih zna ajnih veli ina u ovom eksperimentu.
Slika 22. Grafici zadatog i stvarnog momenta u eksperimentu sa slike 21.Vidimo da na samom po etku ne mo e da se razvije nominalni moment, a to je zbog
toga to smo postavljanjem po etne vrednosti f_r0= 3 ograni ili struju iqs∗ . Oscilacije u
momentu poti u od visokofrekventnih oscilacija u struji, a one su posledica rada strujnogregulatora.U estanost tih oscilacija je u estanost komutacije i zavisi od koraka histerezisa h kojimi odre ujemo. to je taj korak ve i, amplituda oscilacija je ve a, a njihova u estanost je manja.U konkretnom slu aju je postavljeno h=0.05 .
34
Slika 23. Grafici zadatih i stvarnih vrednosti α i β komponente struje statoraVidimo da su zadate vrednosti komponenti struja statora konstantne kada je M em
∗ = 0 .
To je i logi no jer je tada iqs∗ = 0 , pa je i ω k = 0 , a po to se rotor ne obr e sledi ω s = 0 .
Slika 24. Talasni oblici komponenti fluksa statora
3.4. Model digitalnog PI regulatora
35
Model digitalnog PI regulatora dat je na slici 25. Switch blokovi na izlaz propu tajugornji ulaz ako je srednji ulaz ve i od nekog praga ( koji je za oba Switch bloka 0.9 ), ina e sepropu ta donji ulaz.
omega_m*
omega_m^+
-
Sum
z1
Unit Delay1
Tmeh*(kp+ki)/T
Gain
Tmeh*kp/T
Gain1
Switch1
6.56/u^2
Fcn1
8.1
Constant1
f_r*
Switch
sqrt(3)
Constant
Zero-OrderHold
0.9*sqrt(3)/u
FcnAbs
Abs
M_em*+
-Sum1
+
+Sum3 OGRANIC
z1
Unit Delay2
Slika 25. Model brzinskog PI regulatora u SIMULINK-uDigitalni PI regulator smo razmatrali u odeljku 2.7. U odeljku 2.6. je izra unata
maksimalna vrednost momenta koji se mo e razviti na datoj brzini obrtanja i to ograni enjeemo inkorporirati u model regulatora. Da bismo to uradili, funkciju prenosa regulatora emo
predstaviti u obliku :
( ) ( )( ) ( )M em P Iz k z kz
z= − +−
−−1
1
11
1Δω .
(3.6)
Po to deo prenosne funkcije van zagrade predstavlja diskretni integrator, deo u zagradi zajednosa ( )Δω z predstavlja inkrement upravlja ke veli ine. Ograni enje momenta se ostvaruje tako
to se diskretni integrator predstavi preko povratne sprege, to se vidi u modelu regulatora naslici 25. U direktnoj grani se postavi ograni ava koji spre ava integraljenje ako je dostignutamaksimalna vrednost momenta. eljeni fluks se odre uje na osnovu eljene brzine obrtanjamotora. Pod koeficijentima kp i ki koji se pojavljuju u modelu regulatora u SIMULINK-u sepodrazumevaju konstante 0.4 i 0.07 koje su odre ene u odeljku 2.7.
Ovaj model sadr i jedan Subsystem blok koji je nazvan OGRANIC i iji je sadr ajprikazan na slici 26. Prvi ulaz je maksimalni moment koji se mo e razviti za trenutnu brzinuobrtanja, a drugi ulaz je eljeni moment iju vrednost generi e regulator. Blok OGRANIC naizlaz propu ta onu vrednost koja je po apsolutnoj vrednosti manja. Prag Switch bloka je 1.
36
1/u
Fcn1
abs(u)
Fcn2
Switch1
*
1
in_1
2
in_2
sgn(u)
Fcn3
1
out_1
*
M_max
M_em*
Slika 26. Podsistem OGRANIC u modelu regulatora
Da bi regulator mogao da radi, potrebno je da ima informaciju o ugaonoj brzini motora.Tu informaciju daje mu blok- ema sa slike 27. na iji se ulaz dovodi izmereni ugao sa IE. Utrenucima odabiranja uzimaju se odbirci ugla, koji su u realnim vrednostima, a brzina koja je urelativnim jedinicama odre uje se tako to se razlika dve uzastopno izmerene vrednosti uglapodeli periodom odabiranja i nominalnom kru nom u estano u. Ovu blok- emu emo tako emaskirati, a na slici 27. je prikazan i tako dobijen blok. Njegovi parametri su nominalna kru nau estanost i perioda odabiranja, koja je ista kao i za regulator.
teta_m^
Zero-OrderHold
z1
Unit Delay
+
-Sum
1/(T*om_nom)
Gain
omega_m^MERENJE
BRZINE
MERBRZ
Slika 27. Blok- ema mera a brzine i njegov maskirani blok
Kao to e se videti u slede em odeljku, prelazni proces u eksperimentu u kome pratimokako brzina prati referentnu brzinu traje nekoliko desetih delova sekunde, pa emo, da biobezbedili dovoljan kvalitet prelaznog procesa, za periodu odabiranja uzeti T=0.01sec.3.5. Model elektromotornog pogona
Kada se sve dosad opisane komponente pove u u celinu dobije se vektorski kontrolisanelektromotorni pogon sa brzinskom regulacijom. Ulaz takvog jednog pogona je eljena brzinaobrtanja motora, moment optere enja je poreme aj, a va nije izlazne veli ine su stvarna brzinaobrtanja, ugao koji zauzima rotor, elektromagnetni moment koji se razvija, struje statora irotora itd. Model takvog elektromotornog pogona prikazan je na slici 28. i sastoji se od svih dosad napravljenih modela.
37
M_em
omega_m
teta_m
f_beta_s
f_alfa_s
i_beta_s
i_alfa_s
M_opt
u_beta_s
u_alfa_s
ASINHRONI
MOTOR
ASMOT
INKREMENTALNI
ENKODER
IE
INVERTOR
SA STRUJNIM
REGULATOROM
INVERTOR
i_alfa_s*
i_beta_s*
INDIREKTNA
VEKTORSKA
KONTROLA
IVK
M_em*
f_r*
DIGITALNI
PI REGULATOR
REG
omega_m*
MERENJE
BRZINE
MERBRZ
omega_m^
teta_m^
Slika 28. Model vektorski kontrolisanog elektromotornog pogona sa brzinskom regulacijom
Da bismo videli kako ovaj pogon funkcioni e, mo emo izvesti dve vrste eksperimenta.Jedan eksperiment je posmatranje kako pogon reaguje na promenu reference brzine, odnosnokako brzina motora prati zadatu brzinu. U tom eksperimentu opteretni moment je konstantan.Drugi eksperiment se odnosi na pona anje pogona u slu ajevima kada treba odr ati odre enubrzinu motora pri promeni opteretnog momenta. Tu je referentna brzina konstantna. Parametridigitalnog PI regulatora su pode eni tako da sistem u zatvorenoj povratnoj sprezi ima realnepolove, tako da bi trebalo o ekivati neoscilatoran talasni oblik odziva elektromagnetnogmomenta i brzine motora na referencu i poreme aj koji su odsko nog tipa. Za praksu je veomabitno da moment u prelaznom procesu ne menja znak, jer bi u suprotnom zbog neidealne spregemotora i optere enja do lo do pojave mrtvog hoda, odnosno do razlike polo aja motora ioptere enja.
Prvi eksperiment koji emo sprovesti je posmatranje pona anja ovog pogona u slu ajuskokovite promene referentne brzine. Na slici 29. prikazani su grafici referentne, tj. zadatebrzine obrtanja i stvarne brzine, sa podse anjem da je zbog dva para polova potrebno stvarnubrzinu podeliti sa dva. U ovom eksperimentu ulazna veli ina naglo menja vrednost u nekimtrenucima. Ona se analiti ki mo e prikazati kao algebarski zbir nekoliko step-funkcija kojeimaju razli ite amplitude skokova i razli ite trenutke u kojima dolazi do promene vrednosti, alise to u eksperimentu mo e uraditi jednostavnije. Na ulaz se stavi Constant blok, koji na po etkueksperimenta ima vrednost 1. U trenutku u kom elimo da ulazna veli ina promeni vrednostpreko menija Simulation-Pause pauziramo simulaciju, zatim promenimo vrednost izlazaConstant bloka, i onda u meniju sa Simulation-Continue nastavimo simulaciju. Efekat je isti, aovako imamo potpunu slobodu u pogledu trenutka promene izlaza i amplitude te promene.
38
0 1 2 3 4 5
-2
-1
0
1
2
Referenca brzine
t
0 1 2 3 4 5
-2
-1
0
1
2
Stvarna brzina
t
Slika 29. Grafici promene zadate i stvarne brzine motora
Na slici 29. mo emo uo iti da je prelazni proces dosta br i za manje brzine obrtanja, a toje zbog toga to se na manjim brzinama mo e razviti ve i moment.
Slika 30. Grafici promene zadatog i ostvarenog elektromagnetnog momentaNa slici 30. dati su grafici promene zadatog ( generi e ga regulator ) i ostvarenog
elektromagnetnog momenta. Mo e se uo iti kako se zadati moment smanjuje sa pove anjembrzine i prelaskom u oblast slabljenja polja, a to nije zbog zakona upravljanja regulatora vezbog ograni enja momenta koje je inkorporirano u regulator.
39
Na slici 31. je prikazano kako se u ovom eksperimentu menjaju struja i fluks statora.Ove dve veli ine su definisane kao :
i i is s s= +α β2 2 , Ψ Ψ Ψs s s= +α β
2 2 .
(3.7)
Slika 31. Talasni oblici struje i fluksa statora u eksperimentupra enja reference brzine
Grafik promene struje sli an je onom za moment sa tom razlikom to struja ne menjaznak. Osim toga, u stacionarnom stanju struja statora nije jednaka nuli, a to je zbog toga topostoji komponenta struje ids koja generi e fluks. Komponenta iqs tada jeste jednaka nuli. Na
grafiku fluksa jasno mo e da se vidi da se pri radu u oblasti slabljenja polja fluks smanjuje sapove anjem brzine.
Drugi eksperiment se odnosi na pona anje pogona u slu aju naglih promena momentaoptere enja, kada je referentna brzina konstantna i jednaka nominalnoj brzini. Napomenu emoda je rad na nominalnoj brzini, po to slabljenje polja zapo injemo na 0 9. ω nom , tako e rad uoblasti slabljenja polja. Eksperiment obuhvata zaletanje neoptere enog motora do nominalnebrzine, zatim naglo optere ivanje motora momentom 2M nom za koji smo rekli da je zapravonominalni moment sa kojim motor radi, i na kraju naglo ukidanje optere enja. Eksperiment trajeukupno 3sec.
Na slici 32. prikazani su grafici promene brzine motora i opteretnog momenta. U samomeksperimentu za opteretni moment, isto kao i za referentnu brzinu u prethodnom eksperimentu,mo emo da koristimo Constant blok iju vrednost menjamo kad pauziramo simulaciju. Nagrafiku brzine motora se vidi da dolazi do poreme aja u trenucima kada se naglo uklju uje iliisklju uje optere enje, me utim regulator ove poreme aje veoma brzo poni tava, i to bezoscilacija u brzini motora.
40
0 0.5 1 1.5 2 2.5 3-0.5
0
0.5
1
1.5Brzina motora
t
0 0.5 1 1.5 2 2.5 3
0
1
2
Opteretni moment
t
Slika 32. Grafici brzine obrtanja motora i opteretnog momenta ueksperimentu sa konstantnom referentnom brzinom
Na slici 33. mo e se uo iti da, nakon to se motor naglo optereti, elektromagnetnimoment dobija vrednosti ve e od opteretnog momenta, kako bi se nadoknadio gubitak u brzini,a zatim se izjedna ava sa opteretnim momentom.
Slika 33. Talasni oblici zadatog i ostvarenog momentaNa slici 34. su dati grafici struje i fluksa statora u ovom eksperimentu. Struja i fluks su
definisani na isti na in kao i u prethodnom eksperimentu, relacijama (3.7). Struja je opet sli na
41
momentu, sa tim to ne pada na nulu u stacionarnom stanju, a fluks se pri optere ivanju motoraprakti no ne menja.
Slika 34. Talasni oblici struje i fluksa statora
3.6. Primena napravljenih modela u re avanju konkretnih in enjerskih problema
Videli smo kako se modeliranje u SIMULINK-u mo e upotrebiti za re avanje unapredzadatih problema. Ako korisnik eli da napravi svoj model u kojem bi bili kori eni modeli kojisu opisani ovde i koji su prilo eni u fajlovima navedenim u Dodatku, postupak e biti obja njenu slede ih nekoliko re enica. Nakon startovanja MATLAB-a i zatim SIMULINK-a, u menijuFile se odabere New, to zna i da korisnik eli da pravi novi model. Otvara se novi prozor ukojem e se taj model nalaziti. Ako korisnik eli da se u njegovom modelu nalazi modelasinhronog motora, on e zatim otvoriti fajl asmotm.m u kojem se nalazi taj model umaskiranom obliku. Primenom drag’n’drop tehnike se blok ASINHRONI MOTOR dovla i ukorisnikov prozor, i zatim se zatvara fajl asmotm.m, jer on vi e nije potreban. Na isti na in se ukorisnikov model ubacuju i drugi ovde ponu eni modeli. U svoj model korisnik mo e da ubaci ielementarne SIMULINK blokove, tako to e u SIMULINK prozoru otvoriti biblioteku u kojojse nalazi eljeni blok i onda e opet primenom drag’n’drop tehnike dovu i taj blok u svoj model.Da bi se pratio tok simulacije, u model trebaju da se ubace razni Scope blokovi, a tako e iblokovi iz Sources biblioteke koji e predstavljati ulazne veli ine.
Recimo da u eksperimentu sa elektromotornim pogonom u kom se menja referentnabrzina korisnik eli da vidi ta se de ava sa brzinom i strujom statora ako se u jednom trenutkunaglo promeni moment inercije motora, sa vrednosti J na vrednost J J'= 3 . Korisnik da binapravio model treba da uradi slede e stvari : iz fajla elpog.m treba da dovu e modelelektromotornog pogona, a iz biblioteka SIMULINK-a dovu e slede e blokove : Constant,Scope, ToWorkspace, Product, Sum, Fcn. Zatim se blokovi povezuju na eljeni na in. BlokSubsystem se pravi tako to se zaokru ivanjem selektuje grupa blokova, a zatim se u menijuOptions odabere Group. Tako se dobija model prikazan na slici 35.
42
1
Constant1 DIGITALNI
PI REGULATOR
REG
RMS1
ASINHRONI
MOTOR
ASMOT
2
Constant
INKREMENTALNI
ENKODER
IE
i_s
To Workspace6
strujaINVERTOR
SA STRUJNIM
REGULATOROM
INVERTOR
INDIREKTNA
VEKTORSKA
KONTROLA
IVK
MERENJE
BRZINE
MERBRZ
brzina
omega_m
To Workspace1
Slika 35. Eksperiment sa elektromotornim pogonom u kome menjamomoment inercije motora
Blokovi nazvani ‘struja’ i ‘brzina’ na slici 35. su Scope blokovi. Po to se njihov sadr ajne mo e sa uvati da bi kasnije bio preba en u neki dokument i od tampan na tampa u, koristese ToWorkspace blokovi ija je uloga da odbirke veli ine koja se dovodi na njihov ulaz u viduniza prebace u MATLAB-ov radni prostor, gde se dalje mogu obra ivati, odnosno gde se moguprikazati u vidu grafika koji se mo e prebaciti u neki dokument.
Na slici 35. blok tipa Subsystem nazvan RMS1 ima funkciju da na osnovu dvekomponente struje statora ra una iz relacije (3.7) struju statora. Njegov izgled je prikazan naslici 36.
2
in_2
1
in_1 +
+
Sum
*
Product1
*
Product
sqrt(u)
Fcn
1
out_1
i_alfa_s
i_beta_si_s
Slika 36. Izgled podsistema RMS1
Pre nego to se pokrene simulacija, potrebno je dati adekvatne vrednosti parametrimasimulacije, od kojih su najva niji minimalni i maksimalni korak simulacije i metod integracije.Minimalni korak simulacije treba da bude dovoljno mali da se ne propuste zna ajni doga aji usistemu koji se modelira. U ovom slu aju vremenski su kriti ne komutacije invertora, jer imajuveoma visoku u estanost. Minimalni korak simulacije mora biti manji od periode komutacije. Ueksperimentu je uzeto T smin = 5μ . U slu aju izra enih diskontinuiteta, odnosno nelinearnosti usistemu, a takav je i ovaj zbog invertora, ni maksimalni korak simulacije ne sme da bude veliki.Ovde je uzeto da bude jednak minimalnom koraku. Metod integracije se odre uje na osnovutoga kakav je model. Linearnim sistemima, tj. modelima prilago en je Linsim metod, sistemima
43
sa sporom i brzom dinamikom odgovaraju Adams i Gear metodi, a nelinearnim i sistemima sadiskontinuitetima Runge-Kutta metodi.
Simulacija se startuje izborom Start u Simulation meniju. Eksperiment e trajati 4sec. Utrenutku t=0 uklju ujemo pogon, u trenutku t=1 naglo menjamo referentnu brzinu sa ω nom na0, u trenutku t=2 menjamo moment inercije motora, odnosno mehani ku vremensku konstantuTmeh ( proporcionalna je momentu inercije ), a u trenutku t=3 referentnu brzinu menjamo sa 0
na ω nom . Motor e sve vreme biti optere en nominalnim momentom.Promene parametara za vreme simulacije vr e se na na in opisan u prethodnom odeljku.
U glavnom meniju se pod Simulation odabere Pause, a zatim se promeni neki parametar, takoto se otvori dialog box odgovaraju eg bloka i unese nova vrednost. Kada se unesu sve eljene
promene, simulacija se nastavlja biranjem opcije Continue u Simulation meniju. Za referentnubrzinu emo koristiti Constant blok i na opisani na in emo u pojedinim trenucima menjanjemvrednosti konstante koja je parametar ovog bloka simulirati promenu referentne brzine.Promena parametara se mo e izvesti i dok simulacija traje, bez njenog pauziranja. Pauziranjesimulacije je neophodno ako elimo istovremeno da promenimo vi e parametara, po to sampostupak njihove promene zahteva neko vreme.
Promenu momenta inercije emo ubaciti izme u dve promene referentne brzine, kako bivideli efekat promene ovog parametra motora na stacionarno stanje kao i na prelazni proces.Kao veli ine od interesa posmatra emo brzinu obrtanja motora i struju statora koja je dataizrazom (3.7). Na slici 37. prikazani su rezultati dobijeni ovim eksperimentom. Vidi se dapromena momenta inercije motora u stacionarnom stanju ne uti e na brzinu niti na strujustatora. To se moglo i o ekivati po to se iz mehani ke jedna ine, relacije (2.4) i (2.94), vidi damoment inercije ima uticaja samo kada postoji promena brzine, odnosno u prelaznim re imima.
Slika 37. Grafici struje statora i brzine motora
Sa druge strane, vidi se da je prelazni proces znatno lo iji, sa mnogo du im trajanjem i save im preskokom. Tako e je lo iji sa stanovi ta potro nje energije, jer je du e vreme potrebnoinjektovati struju statora. To je posledica razde enosti regulatora, koji nije prilago enpromenjenom momentu inercije motora.
44
4. ZAKLJU AK
Videli smo da je modeliranje i simulacija raznih dinami kih sistema, ne samo asinhronihmotora, u programskom paketu SIMULINK veoma jednostavno, jer se sve odvija grafi kimputem, u Windows okru enju. Tako e injenica da je SIMULINK zapravo dodatak MATLAB-u obezbe uje veliku mogu nost analize i prikazivanja rezultata dobijenih iz modela uSIMULINK-u uz pomo MATLAB-ovih funkcija. Zna i da je obezbe ena i vizuelna lako a ijednostavnost, a i numeri ka snaga.
to se ti e simulacija, odnosno pu tanja u rad modela u SIMULINK-u, obezbe en jeirok dijapazon algoritama za numeri ku integraciju : dva Runge-Kutta algoritma, Adams, Gear,
Euler, Linsim, Adams/Gear. Linsim metod je prilago en linearnim sistemima, Adams i Gearmetodi su prilago eni sistemima koji imaju i brzu i sporu dimaniku, ali nemaju diskontinuitete,dok su Runge-Kutta metodi pogodni za izrazito nelinearne sisteme koji imaju diskontinuitete.Na korisniku je da proceni koji bi metod najvi e odgovarao za simulaciju sa njegovim modelima.
SIMULINK tako e pru a veliku mogu nost eksperimentisanja za vreme simulacija, jerse mogu menjati parametri blokova. Jedino to se ne mo e raditi je da se menja strukturamodela za vreme simulacije.
Tako e dve zna ajne osobine SIMULINK-a, o kojima smo ve govorili, su grupisanje imaskiranje blokova. Maskiranje blokova omogu ava korisniku da kreira nove blokove, sasvojim parametrima. Korisnik mo e da kreira i njihov izgled, koji mo e da bude ak i u vidujednostavnih crte a, kako izgledaju i elementarni blokovi SIMULINK-a. Novi blokovi mogu dase sa uvaju i da se koriste u raznim modelima, to zna i da je obezbe ena modularnost, a toomogu ava pravljenje veoma slo enih modela.
45
DODATAK
Uz ovaj diplomski rad prilo ena je i disketa na kojoj se nalazi tekst diplomskog rada uvidu arhiviranog Word document fajla, kao i m-fajlovi ( ekstenzija ‘m’ je od MATLAB ) kojisadr e sve modele koji su za potrebe ovog diplomskog rada napravljeni u SIMULINK-u. Zanekoliko modela prilo ena su po dva fajla : jedan fajl sadr i model pre maskiranja, koji seobi no sastoji od mno tva elementarnih SIMULINK blokova; drugi fajl sadr i maskirani model,koji se sastoji od samo jednog bloka. Bitno je napomenuti slede e : fajl koji sadr i model premaskiranja nije pogodan za kori enje ( pokretanje simulacije ), ve za upoznavanje sastrukturom modela, a to je zbog toga to su vrednosti parametara blokova zamenjene njihovimimenima kako bi se onda u dialog box-u maskiranog bloka unosile vrednosti tih parametara.Zna i za eksperimente se koriste oni fajlovi koji sadr e maskirane modele.
Samo po jedan fajl prilo en je za eksperimente, a to su eme sa slika 13. , 21. i 35. , kaoi za model elektromotornog pogona sa slike 28.
Na disketi su prilo eni slede i m-fajlovi :
asmot.m - model asinhronog motora sa slike 9.asmotm.m - maskirani model asinhronog motora sa slike 12.invert.m - model trofaznog invertora sa histerezisnim strujnim regulatorom - slika 16.invertm.m - maskirani model invertora sa regulatorom - slika 17.ivk.m - model indirektnog vektorskog kontrolera - slika 18.ivkm.m - maskirani model indirektnog vektorskog kontrolera - slika 19.inkenk.m - model inkrementalnog enkodera - leva polovina slike 20.inkenkm.m - maskirani model IE - desna polovina slike 20.regul.m - model digitalnog PI regulatora - slika 25.regulm.m - maskirani model regulatoramerbrz.m - model mera a brzine - leva polovina slike 27.merbrzm.m - maskirani model mera a brzine - desna polovina slike 27.elpog.m - model elektromotornog pogona sa slike 28.exper1.m - eksperiment sa slike 13.exper2.m - eksperiment sa slike 21.exper3.m - eksperiment sa slike 35.
46
LITERATURA
1. “The Student Edition of SIMULINK - User’s Guide” , The Math Works Inc.2. urkovi Nenad, “Estimacija brzine asinhronog motora bez senzora na
osovini” - diplomski rad3. Mirjana Ne i , “Primena delta modulacije u upravljanju elektromotornog pogona