Pogoni u robotici Pogoni u robotici - FER-a - Sveučilište u Zagrebu
Pogoni u robotici Pogoni u robotici - FER-a - Sveučilište u Zagrebu
Pogoni u robotici Pogoni u robotici - FER-a - Sveučilište u Zagrebu
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
<strong>Pogoni</strong> u <strong>robotici</strong><br />
• Planiranje trajektorije � generiranje reference<br />
(pozicije, brzine, ubrzanja) uz uvažavanje ograničenja<br />
svakog zgloba<br />
• Upravljanje zglobom ≡ upravljanje motorom u zglobu i<br />
radnim mehanizmom pričvršćenim na osovinu motora �<br />
upravljanje slijednim sustavom<br />
• Osnovne veličine s kojima se upravlja: položaj, brzina,<br />
moment (tj. struja), sila dodira<br />
• Mehanički prijenos (reduktor) � smanjuje brzinu zgloba,<br />
a povećava moment (konstantna snaga P=MΩ)<br />
• Iz dinamičkog modela znamo dinamička obilježja zgloba<br />
robota � promjenljivi moment inercije J i koeficijent<br />
viskoznog trenja B, prisutno je trenje i utjecaj sile teže<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
<strong>Pogoni</strong> u <strong>robotici</strong><br />
• Osnovni elementi pogona zgloba robota<br />
– Slijedni regulator (servo regulator)<br />
– Slijedno pojačalo (servo pojačalo snage)<br />
– Slijedni motor (servomotor) s radnim mehanizmom<br />
(člankom robota)<br />
– Mjerni član (enkoder, rezolver, rijeđe potenciometar i<br />
tahogenerator)<br />
1<br />
2
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Vrste pogona u <strong>robotici</strong><br />
• Pneumatski pogoni<br />
– Jeftini, najčešće se ugrađuju u translacijski zglob<br />
– Koriste se u uvjetima gdje može doći do eksplozije<br />
– Loša strana što zahtijevaju stalno održavanje i što pod<br />
djelovanjem tereta ne drže položaj (razlog je u stlačivosti<br />
zraka), razmjerno male snage<br />
• Hidraulički pogoni<br />
– Skuplji, najčešće se ugrađuju u translacijski zglob<br />
– Velike snage (manipulacija s težim predmetima) i dosta su<br />
brzi (brzina vrha alata oko 1m/s)<br />
– Loša strana je potreba za stalnim održavanjem i problemi<br />
s brtvljenjem (curenjem ulja), javlja se i buka<br />
• Električki pogoni<br />
– Dominantno u uporabi, odlične upravljačke karakteristike<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Električki pogoni u <strong>robotici</strong><br />
• Istosmjerni motori s permanentnim magnetima<br />
– Pogodni su za upravljanje<br />
– Snage su reda veličine 100 W do 10 kW<br />
– Nominalna brzina vrtnje motora iznosi 1000 – 3000 min−1 ,<br />
– Nedostatak je u kolektoru, koji određuje vijek trajanja i<br />
maksimalnu struju armature u prijelaznom procesu (2<br />
– 2.5 puta veća od nominalne vrijednosti)<br />
– Nominalno (maksimalno) ubrzanje motora i radnog<br />
mehanizma određeno je omjerom nominalnog<br />
(maksimalnog) momenta motora i momenta inercije<br />
motora i radnog mehanizma<br />
– Elektromehanička vremenska konstanta istosmjernih<br />
motora standardne izvedbe reda je veličine 15 – 150 ms,<br />
a nominalno ubrzanje iznosi (1 – 1.5)103 rad/s2 3<br />
4
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Električki pogoni u <strong>robotici</strong><br />
• Istosmjerni motori specijalne izvedbe<br />
– imaju mnogo veće ubrzanje od strojeva standardne izvedbe,<br />
– smanjenje djelovanja kolektora (izvedbe s glatkim, tj. nenazubljenim<br />
rotorom) - poboljšava se hlađenje vodiča armature i smanjuje se<br />
induktivitet armature � poboljšana komutacijska svojstva - prijelazno<br />
strujno preopterećenje motora 8–10 puta veće od nominalnog<br />
• Dodatno povećanje brzine odziva motora s glatkim rotorom postiže<br />
se izvedbom izduženog rotora manjeg promjera � na taj se<br />
način postiže da motori s glatkim rotorom imaju za red veličine<br />
veće ubrzanje od motora standardne izvedbe.<br />
– smanjenje momenta inercije rotora - rotor kao šuplji cilindar ili disk od<br />
nemagnetskog materijala (namot armature obično štampan ili štancan),<br />
– zamjena kolektora (mehaničkog komutatora) elektroničkim<br />
komutatorom (elektronički komutirani, odnosno beskolektorski motori)<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Električki pogoni u <strong>robotici</strong><br />
• Sinkroni motori s permanentnim magnetima<br />
– Beskontaktni, odnosno elektronički komutirani istosmjerni motori �<br />
motor je bez kolektora, ima elektronički komutator i davač položaja<br />
rotora, koji generira signale za upravljanje elektroničkim komutatorom<br />
– ovisno o rasporedu namota armature statora, motori mogu imati<br />
trapezni (BLDC) ili sinusni (PMSM) oblik protuelektromotorne sile<br />
• Asinkroni motori<br />
– Imaju manje dimenzije i težinu, jeftiniji su od istosmjernih strojeva iste<br />
snage i jednostavniji za održavanje,<br />
– Upravljanje asinkronim motorima je složenije, pa su asinkroni pogoni<br />
nešto skuplji od istosmjernih pogona,<br />
– Metode upravljanja – vektorsko, neposredno upravljanje momentom.<br />
• Koračni motori<br />
– pogodni za neposredno digitalno upravljanje – velika brzina odziva i<br />
točnost pozicioniranja � moguć rad u otvorenoj i zatvorenoj petlji<br />
5<br />
6
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Istosmjerni motori za primjenu u<br />
slijednim sustavima robota<br />
Matematički model za istosmjerni motor:<br />
U općem slučaju vrijedi za bilo koji zglob:<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Istosmjerni motori za primjenu u<br />
slijednim sustavima robota<br />
M t<br />
7<br />
8
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Istosmjerni motori za primjenu u<br />
slijednim sustavima robota<br />
Može se zanemariti<br />
Aperiodski odziv brzine vrtnje<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Istosmjerni motori za primjenu u<br />
slijednim sustavima robota<br />
Kombinacijom i preuređenjem prijenosnih funkcija na način da na jednu stranu dođu nezavisne<br />
varijable U a i M s , dobije se:<br />
Drugi član na lijevoj strani može se preurediti množeći brojnik i nazivnik s I a (K a =1/R a ):<br />
• Omjer viskoznog trenja BΩ prema razvijenom momentu KI a nije nikad veći od 0.1<br />
• Pad napona na otporu armature I a R a znatno je manji po iznosu od susjednog člana KΩ (tj.<br />
protuelektromotorne sile E) � zato je taj član za dodatni red veličine manji u odnosu na E,<br />
pa se može zanemariti u analizi dinamičkog ponašanja.<br />
• Ako se zanemareni član i uzme u obzir, rezultati analize dinamičkog ponašanja neće se<br />
bitno promijeniti, jer se u matematičkom smislu jedino vrijednost koeficijenta K malo<br />
povećava.<br />
9<br />
10
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Istosmjerni motori za primjenu u<br />
slijednim sustavima robota<br />
Ovo je najčešći oblik prijenosne funkcije istosmjernog servomotora u robotičkoj literaturi<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Elektronički komutirani istosmjerni<br />
motori u slijednim sustavima robota<br />
• U osnovnoj izvedbi ovaj motor proizvodi sinusoidalnu<br />
protuelektromotornu silu i zbog toga zahtijeva<br />
sinusoidalne struje statora da bi proizveo konstantan<br />
moment motora.<br />
• Rotirajući magnetski tok u zračnom rasporu<br />
sinkronog motora s permanentnim magnetima<br />
predstavlja osnovnu razliku u odnosu na stacionarni<br />
magnetski tok istosmjernog motora.<br />
• Permanentni magnet na rotoru osigurava konstantnu<br />
uzbudu, dok se u pogledu ostalog ovaj motor ne<br />
razlikuje od klasičnog sinkronog motora.<br />
• Transformacijom jednadžbi sinkronog stroja iz RST<br />
koordinatnog sustava u d, q koordinate sustava čvrsto<br />
vezanog za položaj rotora, sve impedancije sinusno<br />
ovisne o zakretu poprimaju oblik konstanti � najčešći<br />
oblik upravljanja – vektorsko upravljanje.<br />
q<br />
d<br />
11<br />
Vektor magnetskog toka Ψ m smješten<br />
je na d-osi, tj. u smjeru osi rotora.<br />
Komponenta struje statora i q koja<br />
razvija moment motora okomita je na<br />
tok Ψ m i položena je u smjeru q-osi.<br />
12
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Elektronički komutirani istosmjerni<br />
motori u slijednim sustavima robota<br />
• Cilj vektorskog upravljanja je održavanje komponente statorske struje i d što bliže nuli,<br />
kako bi bio ispunjen uvjet za struju statora |i s | = i q .<br />
• Za vektorski upravljan elektronički komutirani motor poželjno je odrediti matematički<br />
model u d, q prostoru, a polaznu osnovu za to predstavlja klasični model sinkronog motora.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Elektronički komutirani istosmjerni<br />
motori u slijednim sustavima robota<br />
• Prevođenje iz RST u d,q prostor pomoću Parkove transformacije:<br />
• Ako se statorske naponske jednadžbe izraze u d, q sustavu, koji je čvrsto vezan za položaj<br />
rotora, tada one poprimaju oblik:<br />
13<br />
približno ispunjeno zahvaljujući vektorskom upravljanju<br />
Uz uvjet da je i d =0,<br />
14
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Elektronički komutirani istosmjerni<br />
motori u slijednim sustavima robota<br />
Uz uvjet da je i d =0,<br />
Primjenom Laplaceove transformacije:<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Prijenos pomoću reduktora<br />
Omjer mehaničkog prijenosa i mp definiran je kao omjer mehaničke brzine vrtnje<br />
motora ω m i brzine zgloba ω i :<br />
znatno smanjena ovisnost elektromehaničke vremenske konstante o položaju robota,<br />
15<br />
16
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Načini upravljanja slijednim sustavima<br />
robotskih mehanizama<br />
Mogući oblik digitalnog sustava upravljanja robotskim mehanizmom:<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Načini upravljanja slijednim sustavima<br />
robotskih mehanizama<br />
• Dominiraju kaskadni sustavi upravljanja<br />
• Dva osnovna pristupa:<br />
• Prvi se pristup sastoji u upravljanju varijablama zgloba pomoću podređene<br />
regulacijske petlje struje, odnosno momenta robotske ruke, pa se moment<br />
motora M m u tom slučaju uzima kao ulazna veličina u zglob.<br />
• U ovakvoj se strukturi čitava regulacijska petlja struje može predstaviti kao<br />
ekvivalentnim strujnim pojačalom, koje kao izlaz daje armaturnu struju.<br />
17<br />
18
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Načini upravljanja slijednim sustavima<br />
robotskih mehanizama<br />
• Dominiraju kaskadni sustavi upravljanja<br />
• Dva osnovna pristupa:<br />
• Drugi se pristup sastoji u upravljanju varijablama zgloba pomoću podređene<br />
regulacijske petlje brzine vrtnje pa se robotski zglob uzima kao uzrok<br />
promjene momenta tereta koji djeluje na osovinu motora<br />
• U ovakvoj se strukturi cijela petlja brzine vrtnje može predstaviti<br />
ekvivalentnim naponskim pojačalom, koje kao izlaz daje promjenljivi napon<br />
armature motora.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Načini upravljanja slijednim sustavima<br />
robotskih mehanizama<br />
• Prvi se pristup koristi kada izvršni dio ruke može doći u kontakt s drugim<br />
objektima koji sprečavaju gibanje duž koordinate, jer tada je podesno<br />
regulirati moment<br />
• Drugi se pristup koristi kada se robotska ruka može slobodno gibati duž<br />
neke koordinate (npr. roboti za bojenje), jer tada je prikladno regulirati brzinu<br />
gibanja.<br />
• Treba primijetiti da je moguće regulirati ili brzinu ili moment, ali ne i oboje<br />
zajedno.<br />
• Pri podređenom upravljanju brzinom vrtnje, izlaz iz mikroprocesorskog<br />
regulatora položaja proporcionalan je zadanoj brzini vrtnje, a u slučaju<br />
podređenog upravljanja momentom, zadanom momentu.<br />
• Tijekom analize dinamičkog ponašanja digitalnih slijednih sustava robotskih<br />
mehanizama koja slijedi efekti diskretizacije su zanemareni, jer se smatra da<br />
je frekvencija diskretizacije znatno manja od presječne frekvencije amplitudno<br />
frekvencijske karakteristike otvorenog kruga razmatranih sustava.<br />
19<br />
20
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Primjeri iz prakse: Projektiranje surlastog robota<br />
Zahtjevi:<br />
• male dimenzije elemenata<br />
• veliki moment izlazne osovine<br />
• visoka preciznost pozicioniranja<br />
• mogućnosti lakog programiranja<br />
• standardni tip brze komunikacije<br />
(npr. CAN sabirnica)<br />
Karakteristike robota:<br />
• sastavljen od osnovnih robotskih modula s 2 SSG<br />
• mogućnost naknadnog dodavanja modula (npr. do 12 SSG)<br />
• opremanje robota s dodatnim senzorima (sonar, kamera, ispitne<br />
sonde) ili alatom (laserski nož, hvataljka)<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odabrano rješenje nakon pretrage<br />
na Internetu<br />
21<br />
22
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljački uređaj Faulhaber MCDC 3003/06 C<br />
Karakteristike: Kaskada P regulatora pozicije, PI regulatora brzine<br />
vrtnje i PWM regulatora struje armature<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
(Robot Kiwi)<br />
Karakteristike: Kaskada<br />
vanjskog regulatora pozicije, PI<br />
regulatora brzine vrtnje i PWM<br />
regulatora struje armature uz<br />
moguće konfiguriranje kaskade<br />
(pomoću DIP1 preklopke)<br />
Maxon 4-Q-DC Servoamplifier ADS<br />
23<br />
24
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Aerotech A3200 Servo controller (robot Hrid)<br />
Karakteristike: Kaskada<br />
PI regulatora pozicije, PI<br />
regulatora brzine vrtnje i<br />
PWM regulatora struje<br />
armature uz moguće<br />
konfiguriranje kaskade te<br />
mogućnosti kompenzacije<br />
i predupravljanja<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Aerotech A3200 Servo controller (robot Hrid)<br />
Potpuno digitalizirani<br />
sustav upravljanja<br />
25<br />
26
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Sinteza nominalnog slijednog regulatora<br />
• Zna se da parametri J u i B ovise o položaju, odnosno brzinama zglobova<br />
robota i da su stoga promjenljivi parametri čije vrijednosti su iz intervala<br />
[J umin , J umax ], odnosno [B min , B max ],<br />
• Za provođenje sinteze slijednog regulatora potrebno je odrediti granice<br />
intervala [J umin , J umax ] i [B min , B max ], te definirati nominalne vrijednosti<br />
parametara J un i B n za svaki robotski zglob � Podloga za određivanje<br />
nominalnih parametara slijednih regulatora pojedinih zglobova robota<br />
• Određivanje parametara � primjena postupaka identifikacije (npr. rekurzivna<br />
metoda najmanjih kvadrata - RLS)<br />
zglob1<br />
zglob2<br />
J umax J un J umin<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Rezultati identifikacije prvog zgloba 4-osnog<br />
SCARA robota Kiwi<br />
manji J, brži odziv - identifikacija u skladu s očekivanjima<br />
27<br />
28
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Metode sinteze slijednih sustava<br />
• Sinteza sustava obuhvaća izbor strukture i parametara<br />
korekcijskih elemenata (regulatora) u svrhu postizanja<br />
zadanih pokazatelja kvalitete upravljanja<br />
• Osnovni elementi sustava (objekt upravljanja, izvršni uređaj,<br />
pojačalo i mjerni elementi) – odabrani prema tehničkim<br />
zahtjevima<br />
• Pretpostavka da su parametri svih osnovnih elemenata<br />
sustava stacionarni<br />
• Odabir linearnog regulatora<br />
– U regulacijskim petljama položaja dominantno se koristi<br />
P i PI regulator<br />
– U regulacijskim petljama brzine vrtnje dominantno se<br />
koristi PI regulator<br />
– U regulacijskim petljama armaturne struje dominantno<br />
se koristi PI regulator, PWM regulator i regulator s<br />
histerezom<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Metode sinteze slijednih sustava<br />
• Podešavanje parametara regulatora – djeluje se na pogrešku<br />
regulacije δ us i nadvišenje σ m. Ako se pri tom ne može<br />
ostvariti potrebna brzina odziva, uvode se unutarnje povratne<br />
veze, tj. primjenjuje se kaskadni sustav regulacije<br />
• U pogledu zahtijevane točnosti, općenito vrijedi δ us ↓ � K o ↑,<br />
ali povećanje K o može uzrokovati nestabilnost<br />
• Zato se postavlja dodatni zahtjev na korekcijski element �<br />
osiguranje amplitudnog i faznog osiguranja<br />
• Za sintezu parametara mogu se primijeniti direktne i<br />
indirektne metode<br />
29<br />
30
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Metode sinteze slijednih sustava<br />
• Direktne metode:<br />
– analitičke metode – pogodne samo za analizu jednostavnijih<br />
sustava automatskog upravljanja, a sasvim nepogodne za<br />
sintezu sustava<br />
– metode simuliranja sustava na računalu – daju neposredan<br />
uvid u dinamičko ponašanje konkretne strukture i parametara<br />
sustava. Može se organizirati iterativni postupak optimiranja<br />
sustava, kojim se uz odabrani kriterij optimiranja dobiju<br />
parametri sustava<br />
• Indirektne metode:<br />
– frekvencijske metode – pokazale su se vrlo pogodnim za<br />
inženjersku primjenu<br />
– metoda mjesta korijena – dosta nepogodna za sintezu<br />
sustava višeg reda<br />
– integralne metode – za jednostavnije sustave gdje se zadani<br />
pokazatelji kvalitete nastoje postići minimiziranjem nekog<br />
integralnog kriterija (IAE, ITAE, ISE, ITSE)<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Frekvencijske metode sinteze PI regulatora<br />
• Tehnički optimum:<br />
– Određivanje integralne vremenske konstante T I – cilj je odrediti<br />
T I koja kompenzira najveću (dominantnu) vremensku<br />
konstantu u regulacijskom krugu<br />
– Određivanje koeficijenta pojačanja K R – cilj je za tako određeni<br />
T I odrediti K R s kojim se ostvaruje zadano fazno osiguranje<br />
• Simetrični optimum:<br />
– Primijenjuje se na astatičke sustave prvog reda – parametri PI<br />
regulatora određuju se tako da Bodeova frekvencijska<br />
karakteristika na presječnoj frekvenciji ω c ima nagib -1 te da<br />
su frekvencije loma lijevo i desno od ω c jednako udaljene<br />
(simetrične) od ω c � širina pojasa određuje fazno osiguranje<br />
– Izbor ω I = 1/T I – frekvencija loma lijevo od ω c<br />
– Izbor K R – određuje presječnu frekvenciju ω c<br />
• Pitanje: Mogu li se oba načina sinteze parametara PI regulatora<br />
primijeniti na pozicijske slijedne sustave?<br />
31<br />
32
Lo<br />
Lz<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Frekvencijske metode sinteze PI regulatora<br />
ωc<br />
Presječna<br />
frekvencija<br />
otvorenog<br />
regulacijskog<br />
sustava<br />
ωg<br />
Granična frekvencija<br />
zatvorenog<br />
regulacijskog<br />
sustava ω g =ω p<br />
LAFK<br />
otvoreni krug<br />
LAFK<br />
zavoreni krug<br />
-3 dB<br />
ω<br />
ω<br />
Pitanje: Kakvo je ponašanje pozicijskog slijednog<br />
sustava s PI regulatorom podešenim<br />
prema “tehničkom optimumu”?<br />
Odgovor: Nestabilno!<br />
Osnovno o frekvencijskim metodama sinteze:<br />
• Pogodno je za sintezu koristiti Bodeov prikaz<br />
(aproksimaciju s pravcima) frekvencijskih<br />
karakteristika otvorenog slijednog sustava<br />
• Osim željenog faznog osiguranja, neophodno je<br />
uspostaviti barem približnu vezu između<br />
frekvencijskih karakteristika otvorenog sustava i<br />
pokazatelja kvalitete prijelazne pojave<br />
zatvorenog sustava<br />
• Iz teorije je poznato da nagib i nivo logaritamske<br />
amplitudno-frekvencijske karakteristike na niskim<br />
frekvencijama određuje pogrešku u stacionarnom<br />
stanju, dok oblik te karakteristike na srednjim<br />
frekvencijama određuje dinamičko ponašanje<br />
sustava<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Frekvencijske metode sinteze - inženjerski pristup<br />
tipični slijedni sustav – primijetiti da je astatički (nagib -1 na<br />
niskim frekvencijama)<br />
• Nomogrami za tipične slijedne sustave - mogu se odrediti pokazatelji<br />
odziva na skokovitu promjenu reference (h m , ω c t m , ω c t p , ω p /ω c � A m , ω m )<br />
• Vrlo često je ispunjeno 1≤ ω 3 /ω c ≤ 4 � normirano vrijeme maksimuma<br />
ω c t m ≈ 3 i malo se mijenja s L 1 i ω 3 � t m ≈ 3/ω c (može poslužiti za<br />
procjenu presječne frekvencije sustava iz vremena maksimuma).<br />
• Iz nomograma slijedi da je ω p (ω g ) ≈ ω c (vidi prethodnu sliku)<br />
33<br />
34
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Frekvencijske metode sinteze – opće preporuke<br />
• Karakteristike nekorigiranog sustava L o (ω) i karakteristika korigiranog sustava L ok (ω)<br />
trebaju se podudarati u što je moguće većem rasponu frekvencija. U suprotnom,<br />
realizacija korekcijskog elementa postaje složena.<br />
• U području niskih frekvencija nagib karakteristika treba iznositi -r·20dB/dek<br />
ω=1 � L o (ω) = 20 log K o , iznos za r ovisi o zahtjevima na δ us , δ uk i δ ud<br />
• Nagib karakteristike na frekvenciji presjeka ω c poželjno je da iznosi -20dB/dek<br />
(postiže se potrebno fazno osiguranje). Širina pojasa s tim nagibom treba iznositi 2-4<br />
oktave (pokriva nadvišenja u području σ m = 0 – 30 %)<br />
Određivanje parametara serijskih korekcijskih elemenata:<br />
Lok ( ω) = Lo( ω) + Lk(<br />
ω)<br />
ϕ ( ω) = ϕ ( ω) + ϕ ( ω)<br />
ok o k<br />
• Karakteristika korekcijskog elementa:<br />
Lok ( ω) = Lok ( ω) −Lo(<br />
ω)<br />
ϕ ( ω) = ϕ ( ω) −ϕ<br />
( ω)<br />
k ok o<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Zadovoljavajuća korekcija često nije moguća samo promjenom pojačanja:<br />
• K o ↑ dva puta � pomak LAFK za 6 dB � promjena ω c � promjena F.O. i A.O. (tj. desi se<br />
određena promjena pokazatelja kvalitete)<br />
• U nekim slučajevima zadana kvaliteta postiže se serijskim korekcijskim elementom oblika<br />
(lead-lag kompenzator, skraćeno L-L)<br />
1+<br />
τ Ts 1+<br />
T1ks Gs () = Kk = Kk<br />
1+ Ts 1+<br />
T s<br />
2k<br />
odnosno,<br />
G( s)<br />
= KC<br />
( s − z0<br />
)<br />
( s − p )<br />
gdje je: p0 = 1/T, z0 = 1/τT, and KC = Kkτ. Određivanje parametara L-L korekcijskog elementa:<br />
• Kk se određuje iz zahtjeva na pogrešku u stacionarnom stanju<br />
• Ako treba smanjiti σm i tm , onda treba povećati F.O. (γ) i ωc � τ > 1, odnosno T1k >T2k , tj.<br />
kompenzator treba imati pozitivnu fazno-frekvencijsku karakteristiku (postaje element s<br />
faznim prethođenjem)<br />
• Pretpostavka: Kk = 1, T1k = 1 (ϕm – maksimalna vrijednost koju dostiže ϕk (ω))<br />
τ = T 1k /T 2k<br />
ϕ m [ o ]<br />
10<br />
55<br />
• Približna procjena djelovanja DT 1 i PT 1 elemenata na fazu: ±45 o /dek, ±15 o /okt<br />
5<br />
42<br />
2<br />
19<br />
0<br />
35<br />
36
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• U slučaju kad je τ > 1, iznos nule zo je manji od iznosa pola po . Zato L-L kompenzator<br />
pomiče krivulju mjesta korijena prema lijevoj poluravnini. To rezultira povećanim F.O. i<br />
bržim odzivom sustava.<br />
• U frekvencijskom području kompenzator dodaje pozitivnu fazu na području frekvencija<br />
1/τT do 1/T, odnosno 1/T1k do 1/T2k . Iz tablice je vidljivo da ovisno o parametru τ,<br />
maksimum dodane faze je u pravilu manji od 90o . Za veće vrijednosti potrebno je koristiti<br />
dva ili više kompenzatora.<br />
• Maksimalni iznos faze dodan je na središnjoj frekvenciji:<br />
1<br />
ω n =<br />
T τ<br />
• Maksimalni Iznos dodane faze može se izračunati iz vrijednosti parametra τ<br />
τ −1<br />
sin<br />
φ =<br />
τ + 1<br />
• Dodani pozitivni iznos faze povećava F.O. i tako povećava stabilnost sustava. L-L<br />
kompenzator s faznim prethođenjem projektira se određivanjem τ iz potrebnog iznosa<br />
faze da se zadovolje zahtjevi na F.O, te određivanjem T, odnosno T 2k , tako da se<br />
potrebna promjena faze doda na mjestu nove presječne frekvencije – može se koristiti<br />
prethodno spomenuta procjena ±45 o /dek, ±15 o /okt<br />
• Nedostatak korekcije uz τ > 1 je povećanje pojačanja na visokim frekvencijama (dakle,<br />
šuma i smetnji – pri skokovitoj promjeni ulazne veličine pojačanje za ω>> iznosi τ·K k (dok<br />
je K k za ω=0)<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Ako je kvaliteta prijelazne pojave zadovoljavajuća (σ m , t m ), a potrebno je povećati<br />
točnost u stacionarnom stanju, onda se koristi L-L kompenzator s faznim zaostajanjem<br />
� τ < 1, odnosno T 2k >T 1k . S obzirom da je iznos nule z o veći od iznosa pola p o ,L-L<br />
kompenzator pomiče krivulju mjesta korijena udesno bliže ishodištu fazne ranine, što<br />
nije poželjno. Zato u pravilu z o i p o moraju biti smješteni tako da bi što manje mijenjali<br />
dinamiku sustava, a da bi se poboljšao odziv u stacionarnom stanju.<br />
• Na visokim frekvencijama L-L kompenzator s faznim zaostajanjem ima jedinično<br />
pojačanje, a na niskim frekvencijama pojačanje je iznosa z 0 /p 0 > 1, pa je time za toliko<br />
manja i pogreška u stacionarnom stanju.<br />
• Pretpostavka: K k = 1, T 1k = 1 (ϕ m – maksimalna vrijednost koju dostiže ϕ k (ω))<br />
τ = T 1k /T 2k<br />
ϕ m [ o ]<br />
ok<br />
o<br />
1<br />
0.1<br />
-55<br />
2k 1k<br />
c<br />
0.2<br />
-42<br />
3<br />
0.5<br />
-19<br />
37<br />
38
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
•<br />
•<br />
•<br />
Primjena L-L na astatičke sustave – da se smanji kinetička pogreška δuk �Kok > Ko σm nepromijenjeno � γ nepromijenjeno<br />
tm nepromijenjeno � ωc nepromijenjeno<br />
1<br />
ω 1k<br />
=<br />
T1k<br />
ωc<br />
≤<br />
10<br />
osigurava se da L-L ne djeluje na F.O.<br />
• Ako je r=0, a postoji uvjet δus = 0, tada korekcijski element mora sadržavati integralno<br />
ponašanje:<br />
⎛ 1 ⎞ ⎛ 1 ⎞ 1+<br />
Ts I<br />
GR() s = Gk() s = Kk ⎜1+ ⎟= KR⎜1+ ⎟=<br />
KR<br />
⎝ Ts k ⎠ ⎝ Ts I ⎠ Ts I<br />
• Proporcionalna komponenta dodaje se zato da se isključi integralno djelovanje u<br />
području srednjih frekvencija (oko ω c ).<br />
• PI djelovanjem se uz isto γ, postiže veći ω c nego samo I djelovanjem � brži odziv uz<br />
isto nadvišenje (manji t m )<br />
• Određivanje integralne vremenske konstante:<br />
– Statički (brzinski) slijedni sustavi r=0 � T I = T max � tehnički optimum<br />
– Astatički (pozicijski) r=1 � AFK simetrična oko ω c � simetrični optimum<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Povezivanje frekvencijske i vremenske domene<br />
• Ocjena pokazatelja kvalitete slijednih sustava, koji su u pravilu sustavi višeg reda,<br />
može se približno izvršiti na temelju ocjene pokazatelja kvalitete referentnog<br />
sustava drugog reda – ideja je uspostaviti vezu između nadvišenja σ m u odzivu<br />
zatvorenog kruga i faznog osiguranja γ kod frekv. karakteristika otvorenog kruga<br />
• Za referentni sustav drugog reda vrijede sljedeće relacije:<br />
– za Bodeov prikaz aproksimiran pravcima:<br />
[ ]<br />
o<br />
σm% + γ ⎡<br />
⎣<br />
⎤<br />
⎦ ≈63<br />
– za točan Bodeov prikaz:<br />
o<br />
σm[ % ] + γ ⎡<br />
⎣<br />
⎤<br />
⎦<br />
≈70<br />
• Ove relacije mogu se koristiti za procjenu nadvišenja sustava drugog i višeg reda,<br />
odnosno za procjenu potrebnog faznog osiguranja, da bi odziv sustava imao<br />
zadano nadvišenje<br />
• Relacije za t m i σ m mogu se koristiti u postupku sinteze za određivanje<br />
frekvencijske karakteristike otvorenog sustava u području srednjih frekvencija<br />
prema pokazateljima kvalitete prijelazne pojave<br />
39<br />
40
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Komentari na tehnički i simetrični optimum<br />
Zadane pokazatelje kvalitete želi se postići sintezom parametara PI regulatora K R i T I<br />
– slučaj upravljanja brzinom vrtnje, strujom armature – upravljanje<br />
statičkim sustavom koji ne sadrži niti jedan element s integralnim<br />
djelovanjem � otvoreni krug ima konačno pojačanje L(0) na niskim<br />
frekvencijama (nagib 0). Dodavanjem PI regulatora se želi to<br />
promijeniti (dodaje se integralno djelovanje i dodatna nula u otvoreni<br />
krug.<br />
– slučaj upravljanja položajem, silom dodira – upravljanje astatičkim<br />
sustavom koji sadrži jedan element s integralnim djelovanjem �<br />
otvoreni krug ima beskonačno pojačanje L(0) na niskim frekvencijama<br />
(nagib -1). Dodavanjem PI regulatora dodaje se još jedno integralno<br />
djelovanje i dodatna nula u otvoreni krug.<br />
U slučaju tehničkog optimuma kada je T I = T k =T max =T 1 , prijenosne funkcije<br />
nekorigiranog i korigiranog sustava 2. reda imaju oblik:<br />
K<br />
o<br />
ok<br />
G o(<br />
s)<br />
=<br />
Gok<br />
( s)<br />
= Gk<br />
( s)<br />
⋅Go<br />
( s)<br />
=<br />
Kok<br />
=<br />
( 1+<br />
T1s)(<br />
1+<br />
T2s)<br />
s(<br />
1+<br />
T2s)<br />
Iako samo po sebi vidljivo, vrijedi istaknuti da je regulirani slijedni sustav statički<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Bodeov prikaz:<br />
20 log KR/TI<br />
20 log Ko<br />
20 log KR/TI<br />
L(ω)<br />
Lo<br />
ω=1<br />
Lk(ω)<br />
ω=1<br />
Lok<br />
-1<br />
ω1=1/T1<br />
ωI=1/TI<br />
-1<br />
ωc<br />
ωck<br />
20 log KR<br />
K<br />
ωck > ωc<br />
ω2=1/T2<br />
-2<br />
ω<br />
K RK<br />
T<br />
ω<br />
I<br />
o<br />
41<br />
42
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Ako je ω ck = ω 2 /2 � γ ≈ 60 o , pa nadvišenje iznosi σ m ≈ 64 - γ = 4%<br />
• Za sustav 2.reda � ξ = 0.707 (σ m = 4,3%) � Takav način projektiranja regulatora se<br />
u klasičnoj literaturi naziva projektiranje prema tehničkom optimumu<br />
• Tehnički optimum u proširenom smislu – pojačanje regulatora u skladu sa zadanim<br />
σ m , odnosno γ (relacija γ ≈ 63 – σ m )<br />
• Postiže se presječna frekvencija ω ct (t mt ≈ 3.3/ω ct )<br />
• Realni slijedni sustav n ≥3:<br />
• Projektiranje prema tehničkom optimumu:<br />
Ko<br />
Go() s =<br />
(1 + Ts 1 )(1 + Ts 2 )(1 + Ts 3 )<br />
G () s = G () s G () s =<br />
K<br />
K<br />
K K<br />
=<br />
ok R o<br />
ok k o<br />
s(1 + Ts 2 )(1 + Ts 3 )<br />
ok<br />
TI<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Bodeov prikaz:<br />
20 log KR/TI<br />
20 log Ko<br />
20 log KR/TI<br />
L(ω)<br />
Lo<br />
ω=1<br />
Lk(ω)<br />
ω=1<br />
Lok<br />
-1<br />
ω1=1/T1<br />
ωI=1/TI<br />
-1<br />
ωc<br />
ωck<br />
ωck > ωc<br />
20 log KR<br />
ω2<br />
-2<br />
ω3<br />
-3<br />
ω<br />
ω<br />
43<br />
44
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Za sustave višeg reda (n ≥3) � određivanje koeficijenta pojačanja regulatora<br />
(korekcijskog elementa) prema približnoj relaciji za određivanje γ iz zadanog<br />
nadvišenja σ m :<br />
γ[ o ] ≈ 58 - σ m<br />
• Također približna relacija za određivanje vremena maksimuma:<br />
t mt ≈ 2.9/ω ct<br />
• Sažetak: Projektiranje PI regulatora prema tehničkom optimumu može se<br />
primijeniti samo na statičke objekte upravljanja, jer astatičke sustave čini a priori<br />
nestabilnima (zbog dva integratora u otvorenom krugu faza na niskim<br />
frekvencijama kreće od -180 o )<br />
Kako stoje stvari s određivanjem K R i T I prema simetričnom optimumu, kada se želi<br />
dobiti simetričnu LAFK oko ω c ?<br />
• Klasična definicija kaže da se ovaj optimum primjenjuje na astatički sustav prvog<br />
reda, tj. za r=1. Zašto za r=1? Zato, jer tehnički optimum nije primjenjiv, pa treba<br />
primijeniti neki drugi kriterij.<br />
• Zašto ne za r=2,3,..? Jasno, zbog problema nestabilnosti.<br />
• Za r=1 u otvorenom krugu su dva integratora pa treba izborom parametara K R i<br />
T I osigurati potrebnu stabilnost sustava, tj. potrebno F.O. na ω c .<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Projektiranje prema simetričnom optimumu za astatički sustav (r=1):<br />
Ko<br />
Go() s =<br />
s(1 + Ts 1 )<br />
Najjednostavniji slučaj – inače se na T1 gleda<br />
kao na dominantnu vremensku konstantu Tmax Kok (1 + Ts I )<br />
Gok () s = Gk() s Go() s = 2<br />
s (1 + Ts) KRKo Kok<br />
=<br />
T<br />
• AFK korigiranog sustava je simetrična, ako su ω I =1/T I i ω 1 =1/T 1 jednako udaljene od<br />
ω c :<br />
ω ω<br />
=<br />
ω ω<br />
I cs =<br />
cs<br />
I<br />
• Širina područja nagiba (-1) iznosi:<br />
• Faza na frekvenciji ω cs iznosi:<br />
a<br />
ω T ω ω<br />
ω ω ω<br />
1 I 1 cs<br />
= = ⋅ =<br />
I T1<br />
cs I<br />
a<br />
2<br />
1<br />
o o<br />
1<br />
ϕ( ωcs ) =− 180 + arctgTIωcs− arctgT1ωcs =− 180 + arctga−arctg a<br />
I<br />
45<br />
46
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Bodeov prikaz:<br />
20 log KR/TI<br />
20 log Ko<br />
20 log KR/TI<br />
L(ω)<br />
Lo<br />
ω=1<br />
Lk(ω)<br />
ω=1<br />
Lok<br />
-1<br />
-2<br />
-1<br />
ωc<br />
ωI=1/TI<br />
ωI=1/TI<br />
ωcs > ωc<br />
ωcs<br />
ω1=1/T1<br />
20 log KR<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
1<br />
a -<br />
o<br />
180 ( ) arctg arctg arctg a −<br />
γ = + ϕ ωcs<br />
= a − = = arctg<br />
a 1<br />
1+<br />
a ⋅<br />
2a<br />
a<br />
• Fazno osiguranje iznosi:<br />
1+ sinγ<br />
a =<br />
cosγ<br />
• Približne vrijednosti γ i σ m u odnosu na razne vrijednosti koeficijenta a:<br />
a<br />
γ [ o ]<br />
σ m [%]<br />
• Približne relacije za γ i σ m (vrijede za a > 1):<br />
γ<br />
γ ≈15 a, odnosno a≈<br />
15<br />
64 −σ<br />
m<br />
σm≈64 −γ ≈64 −15 a, odnosno a≈<br />
15<br />
2<br />
30<br />
34<br />
-2<br />
2<br />
1 a 1<br />
3<br />
45<br />
19<br />
4<br />
60<br />
4<br />
ω<br />
ω<br />
47<br />
48
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
σ m [%]<br />
γ [ o ]<br />
a<br />
a 2<br />
10<br />
54<br />
3.6<br />
12.96<br />
U pravilu se kod sinteze PI regulatora prema simetričnom optimumu projektira<br />
sustav za nadvišenje veće od zadanog (σm > σmz ), npr. σm = k σmz , gdje je k=2, 3, 4,<br />
5, pri čemu vrijedi ograničenje σm ≤ 60%. Zadano nadvišenje ostvaruje se<br />
dodavanjem filtra na ulaz referentnog signala:<br />
1<br />
Gf() s =<br />
1+<br />
Ts<br />
f<br />
20<br />
44<br />
2.93<br />
8.6<br />
30<br />
34<br />
2.26<br />
5.14<br />
Primjer: Neka je sustav projektiran za nadvišenje σ mω =30%. Za taj slučaj vrijede<br />
odnosi prikazani u tablici:<br />
σ mz [%]<br />
σ m /σ mz<br />
T f<br />
30<br />
1<br />
0<br />
20<br />
2/3<br />
t m /4<br />
10<br />
1/3<br />
t m /2<br />
5<br />
1/6<br />
2t m /3<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Pseudoderivativni regulator s predupravljanjem (PDFF)<br />
• PDFF regulator može zamijeniti standardne PI regulatore u slučajevima kad su<br />
parametri PI regulatora dobro podešeni u odnosu na smetnju, ali zahtijevaju dodatno<br />
nezavisno podešenje u odnosu na referencu.<br />
• Struktura PDFF regulatora s nazivima ulaza i izlaza sukladnim nazivima blokova za<br />
PLC:<br />
PV ( s)<br />
⎛ KP ⎞<br />
OUT ( s)<br />
= = P ⋅ KP ⋅ SP(<br />
s)<br />
− KP ⋅ PV ( s)<br />
+ ( SP PV )<br />
G ( s)<br />
⎜<br />
TI s ⎟ −<br />
p<br />
⎝ ⎠<br />
SP P⋅ KP<br />
OvisnooparametruP, struktura<br />
regulatora se mijenja i ima različit<br />
utjecaj na dinamičko ponašanje<br />
SP-PV<br />
KP<br />
TI s<br />
+<br />
+<br />
-<br />
OUT<br />
slijednog sustava<br />
Za P=1 � PI regulator<br />
PV KP<br />
⎛ K ⎞ R<br />
Go() s ⎜P⋅ KR<br />
+ ⎟<br />
PV () s<br />
TIs Gz() s = =<br />
⎝ ⎠<br />
SP() s ⎛ K ⎞ R<br />
1 + Go( s) ⎜KR + ⎟<br />
⎝ TIs⎠ 49<br />
50
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Treba primijetiti da podesivi parametar P ne utječe na polove zatvorenog sustava, ali<br />
zato utječe na nulu:<br />
⎛PTs ⋅ I + 1⎞<br />
KR⋅Go() s ⎜ ⎟<br />
PV () s<br />
TIs Gz() s = =<br />
⎝ ⎠<br />
SP() s ⎛Ts I + 1⎞<br />
1 + KR⋅Go( s)<br />
⎜ ⎟<br />
⎝ TIs ⎠<br />
Budući da parametar P utječe na slijedni sustav putem predupravljanja, on ne utječe<br />
na prijenosnu funkciju zatvorenog kruga s obzirom na smetnju. S druge strane,<br />
njegov utjecaj na nulu zatvorenog kruga neposredno utječe na dinamiku zatvorenog<br />
kruga u odnosu na referencu. To znači da PDFF regulator omogućava nezavisno<br />
podešavanje dinamike odziva slijednog sustava na promjenu reference.<br />
Kod praktične realizacije, isto kao kod PI regulatora, treba izvršiti provjeru izlaza u<br />
odnosu na postavljena ograničenja, te u slučaju ulaska i zadržavanja izlaza<br />
regulatora u limitu, potrebno je zaustaviti integriranje (anti-windup).<br />
Preporuka: Odrediti K R i T I kao za standardni PI regulator koristeći simetrični (za<br />
astatičke procese) ili tehnički optimum (za statičke procese), zatim mijenjati P dok se<br />
ne dobiju zadovoljavajući odzivi.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Primjer: PDFF regulacija brzine vrtnje istosmjernog motora (Matlab+Simulink)<br />
regulator<br />
+ pojačalo<br />
snage<br />
Parametri određeni prema tehničkom optimumu: K R =5.8, T I =0.37 s<br />
Simulacijski eksperimenti: P=0; 0,25; 0.5; 0,75; 1.<br />
Promjena reference: ΔU r =0.65 V (ΔΩ=10 s -1 )<br />
Promjena poremećajne veličine (momenta tereta): ΔM t =1 Nm<br />
51<br />
52
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Promjena reference<br />
Promjena tereta<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Sustavi upravljanja zglobom robota uz<br />
upravljanje momentom<br />
• Kaskadna struktura s vanjskom petljom po položaju i unutarnjom petljom po<br />
struji (momentu)<br />
• Osnovni cilj upravljanja kretanjem robotske ruke predstavlja ostvarivanje<br />
potrebnog momenta u svakom zglobu da se pokrene manipulirani teret i da<br />
se pri tom svladaju svi momenti trenja, sile teže i dinamički momenti uslijed<br />
inercije, koji se suprotstavljaju gibanju.<br />
• Osnovna shema upravljanja sadrži P ili PD regulator položaja (vidi sliku)<br />
• PD regulator položaja se u praksi realizira kao paralelni rad P i D regulatora<br />
položaja, odnosno P regulatora položaja i P regulatora brzine<br />
strujna petlja<br />
53<br />
zanemareno je viskozno trenje<br />
54
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Analiza statičkog i dinamičkog ponašanja<br />
• prijenosna funkcija zatvorene regulacijske petlje ima oblik:<br />
• gdje je:<br />
• Iz karakteristične jednadžbe<br />
• slijedi da je<br />
Koeficijent prigušenja ovisi o<br />
ukupnom momentu inercije J u , koji<br />
se mijenja s promjenom položaja<br />
robotske ruke kao i radnog tereta.<br />
Za dobro projektirani manipulator,<br />
J u se mijenja u odnosu 10 : 1.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Analiza statičkog i dinamičkog ponašanja<br />
• Zbog promjenljivog J u postoji problem određivanja parametara<br />
nominalnog regulatora (vezanih uz nominalnu vrijednost J un )<br />
• Obično se odabire da je J un = J sr<br />
• U literaturi se mogu pronaći preporuke da se za J u = J sr regulacijska<br />
petlja podesi na kritičnu vrijednost, tj. ξ =1.<br />
• Tada vrijedi:<br />
• Budući da je ξ sr = 1, koeficijent prigušenja za najmanji moment inercije<br />
(J sr /3) je ξ = 1.73, a za najveći (3J sr ) je ξ = 0.58.<br />
• Prirodna frekvencija neprigušenih oscilacija za J sr označena je s ω no<br />
• Shodno tome, frekvencija neprigušenih oscilacija iznosi ω n = 1.73 ω no<br />
za najmanji, a ω n = 0.58 ω no za najveći moment inercije.<br />
• Manji moment inercije uzrokuje brži odziv iveći koeficijent prigušenja,<br />
a time i manje nadvišenje u odzivu.<br />
• Treba uočiti da uvijek postoji nadvišenje zbog D djelovanja regulatora<br />
55<br />
56
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odzivi slijednog sustava zgloba s nominalnim regulatorom položaja<br />
uz regulaciju momenta<br />
σ m=6%<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Ovisnost promjene momenta motora o promjenama referentne<br />
vrijednosti položaja i momenta tereta može se izraziti sljedećom<br />
prijenosnom funkcijom:<br />
• Poremećajni moment M s sadrži dinamičke i statičke komponente<br />
� u stacionarnom stanju djeluju samo statičke komponente<br />
(prvenstveno moment sile teže M g), pa vrijedi M m=M g<br />
• Postojanje pogreške položaja u stacionarnom stanju (s=0)<br />
• Opisani nedostaci osnovne regulacijske petlje predstavljaju<br />
zapreku pri primjenama robota koje zahtijevaju preciznost<br />
pozicioniranja i čim bolja svojstva dinamičkog momenta motora.<br />
57<br />
58
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljačka petlja po momentu s<br />
kompenzacijskim proširenjem<br />
• Estimirana vrijednost momenta inercije uvedena je kao programirano pojačanje J ue.<br />
• Ovo pojačanje može se uvesti na način, kao što je prikazano na slici dolje, tj.<br />
neposredno iza PD regulatora, ili pak samo u bloku predupravljanja ubrzanjem (tj.<br />
kao J ues 2 ).<br />
• Estimirana vrijednost statičkog momenta ovisnog o djelovanju sile teže M ge kao<br />
dodatni korekcijski signal na izlazu regulatora.<br />
• Predupravljanje ubrzanjem, dodano sa svrhom da se poboljša točnost u dobivanju<br />
zadanog dinamičkog momenta M d.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljačka petlja po momentu s<br />
kompenzacijskim proširenjem<br />
• Nakon uvođenja navedenih elemenata, jednadžba slijednog regulatora<br />
poprima oblik:<br />
• Jednadžba zatvorenog regulacijskog kruga položaja poprima oblik:<br />
• Iz karakteristične jednadžbe sustava može se odrediti koeficijent<br />
prigušenja i prirodna frekvencija neprigušenih oscilacija:<br />
M d<br />
59<br />
60
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljačka petlja po momentu s<br />
kompenzacijskim proširenjem<br />
• Uz uvjet (prikladan izbor vrijednosti pojačanja<br />
strujne petlje)<br />
• Jednadžba zatvorenog regulacijskog kruga položaja poprima<br />
oblik:<br />
• pogreška položaja u stacionarnom stanju pri jediničnoj<br />
skokovitoj promjeni referentne veličine ima oblik:<br />
• Uz uvjet da je M ge =M g i J ue =J u , slijedi da je razvijeni moment<br />
idealizirani slučaj koji<br />
realno nije ispunjen<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljačka petlja po momentu s<br />
kompenzacijskim proširenjem<br />
• Ako se definira omjer estimirane i stvarne vrijednosti c= J ue/J u, može se promatrati<br />
utjecaj odstupanja njegove realne vrijednosti (c≠1) na odziv položaja zgloba robota<br />
• Karakteristično je da pri skokovitoj promjeni reference Δθ odzivi položaja poprimaju<br />
u početku vrijednost θ(0) = cΔθ (to je tako jer u lineariziranom modelu derivativno<br />
djelovanje nema limita) � kad je stvarna vrijednost manja od estimirane javlja se<br />
nadvišenje koje može biti jako veliko (čak i u realnim uvjetima uz uvažavanje svih<br />
ograničenja u slijednom sustavu)<br />
• Postoji potreba za postupkom podešenja nominalnog slijednog regulatora koji će<br />
riješiti ili ublažiti problem mogućeg nadvišenja u odzivu položaja � predlaže se tzv.<br />
min-max postupak određivanja parametara slijednog regulatora K p i K d<br />
– Uz fiksno J ue � najbolje ponašanje uz J ue =J min (c≤1)<br />
– Minimalni koeficijent prigušenja ξ min dobije se uz J ue =J max<br />
– Ideja je da se postavi ξ min > 1 tako da su sve moguće vrijednosti prigušenja uvijek<br />
veće od 1<br />
– Da bi se ovaj postupak mogao primijeniti treba dobro poznavati (J min , J max )<br />
61<br />
62
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu<br />
• Nedostaci koncepta linearnih metoda upravljanja položajem uz<br />
upravljanje momentom<br />
– zbog ograničenih mogućnosti digitalnih računala nemoguće je<br />
eliminirati utjecaj poremećajne veličine M s (ne mogu se točno i<br />
pravovremeno estimirati sve komponente od M s)<br />
– traže se učinkovitiji algoritmi koji se zasnivaju na drugim načelima<br />
korekcije poremećajne veličine<br />
• Robusno upravljanje<br />
– Eliminacija pogreške nastale u sustavu uslijed djelovanja<br />
poremećajne veličine predstavlja osnovni cilj razrade robusnog<br />
algoritma upravljanja.<br />
• Adaptivno upravljanje<br />
– Uklanjanje razlika u kvaliteti odziva sustava nastalih zbog<br />
promjena parametara sustava predstavlja osnovni cilj kod razrade<br />
adaptivnog algoritma upravljanja.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
• Dinamička analiza sustava upravljanja položajem uz upravljanje<br />
momentom uz korištenje istosmjernih motora<br />
• Promjene parametara zgloba robota<br />
• Zadatak je robusnog i adaptivnog algoritma upravljanja poništiti<br />
djelovanje veličine h, tj. bilo kakve promjene parametara ili bilo<br />
kakvog poremećaja, koji se javlja u sustavu ili djeluje izvana.<br />
63<br />
64
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
• Moguć način kompenzacije h je dodavanje korekcijskog signala na izlaz iz<br />
regulatora:<br />
v je dodatni korekcijski signal<br />
• Treba primijetiti da uvođenje dodatnog upravljačkog signala v podrazumijeva<br />
potpunu kompenzaciju djelovanja sile teže, što povlači za sobom ispuštanje<br />
potrebe za dodatnim signalom estimiranog momenta sile teže (tj. M ge = 0).<br />
• Stoga je nominalni slijedni regulator opisan sljedećom jednadžbom:<br />
• Slijedi da je f upravo onaj zakretni moment kojeg motor treba ostvariti da bi sustav<br />
zadržao zadano, tj. nominalno dinamičko ponašanje<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
• Uz uvjet J ue=J u=J un slijedi da jednadžba signala razlike sustava ima<br />
oblik:<br />
Uz uvjet K pK=1<br />
• Slijedi da korekcijski signal v ima oblik:<br />
• Proračun dodatnog korekcijskog signala v može započeti sljedećim<br />
estimacijskim postupkom:<br />
On se zasniva na ideji da je v upravo ona korekcija<br />
koja osigurava nominalno dinamičko ponašanje<br />
65<br />
66
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
• Blokovska shema upravljanja<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
• Robusnost upravljanja zasniva se na postojanju lokalne petlje<br />
beskonačnog pojačanja čiji ulaz i izlaz su ograničeni realnim<br />
ograničenjima slijednog sustava – ostvaren je uvjet za ostvarenje<br />
točnosti u stacionarnom stanju<br />
67<br />
68
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
Prikaz sheme robusnog upravljanja u obliku sheme adaptivnog upravljanja - signalna<br />
adaptacija<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
Prikaz sheme robusnog upravljanja u obliku sheme adaptivnog upravljanja -<br />
parametarska adaptacija<br />
+ m s<br />
69<br />
oblik v isti je kao i kod signalne adaptacije, jer<br />
se zbog robusnosti m s može zanemariti<br />
70
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Robusno i adaptivno upravljanje položajem uz<br />
upravljanje po momentu – Hsia metoda<br />
Praktična realizacija Hsia regulatora<br />
– Limit izlaza iz regulatora (limit struje)<br />
– Realni derivator u adaptacijskom mehanizmu<br />
– Korištenje signala povratne veze brzine vrtnje<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Regulacijske petlje robotskih sustava CNC tipa<br />
• Blokovska shema slijednog sustava CNC tipa<br />
– Kaskadna regulacija s dvije regulacijske petlje brzine vrtnje<br />
– Vanjska petlja s integralnim regulatorom ili PI regulatorom<br />
– Unutarnja petlja s proporcionalnim regulatorom<br />
• Uz referentni signal brzine vrtnje šalje se i referentni signal<br />
položaja θ r upotrebom interpolatora iz CNC nadzornog računala.<br />
• U svakom intervalu diskretizacije uspoređuje se referentni signal<br />
položaja sa stvarnim položajem, koji je izražen s brojem impulsa<br />
iz inkrementalnog davača impulsa.<br />
• Kada su obje vrijednosti jednake, tada se zglob nalazi u zadanom<br />
položaju i zaustavlja se referentni signal brzine vrtnje.<br />
71<br />
72
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Regulacijske petlje robotskih sustava CNC tipa<br />
• Generiranje reference brzine vrtnje<br />
• U t 1 razlika između zadanog i stvarnog<br />
položaja iznosi e (vidi sliku desno), a regulator<br />
ju otklanja u intervalu između trenutaka t 1 i t 2<br />
• Zadovoljavajuće ponašanje robotskih sustava<br />
CNC tipa dobije se za:<br />
• Vrijednost pojačanja K o podešava se za<br />
maksimalnu očekivanu vrijednost<br />
elektromehaničke vremenske konstante T m ,<br />
koja ovisi o maksimalno očekivanom iznosu<br />
ukupnog momenta inercije J u .<br />
• Za manje vrijednosti ukupnog momenta<br />
inercije koeficijent prigušenja će poprimiti<br />
vrijednosti veće od 0.7.<br />
• Za male momente inercije koeficijenti<br />
prigušenja mogu značajno porasti �<br />
posljedica je veoma trom prijelazni proces i<br />
usporeni radni ciklus robota.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira manipulatora<br />
• Kretanje manipulatora može biti bez ograničenja ili s ograničenjima<br />
• Primjeri ograničenog kretanja:<br />
– gibanje po površini stola (pisanje, brušenje, rezanje)<br />
– guranje i povlačenje pretinaca<br />
– različiti zadaci montaže (umetanje, zavrtanje)<br />
• Razlozi za uvođenje regulacije sile dodira:<br />
– Robotski poslovi (npr. proces brušenja) kod kojih upravljanje po trajektoriji gibanja<br />
nije dovoljno, jer se alat kojim se radi troši pa se nakon nekog vremena, ukoliko se<br />
upravlja samo po trajektoriji gibanja, izgubi kontakt između predmeta i alata,<br />
– Nepoznavanje karakteristika obrađivanog predmeta (npr. hrapavosti i oblika ruba<br />
brušenog predmeta).<br />
– Robotski poslovi za koje je trajektorija zadana uz pretpostavljeno poznavanje radne<br />
okoline robotskog sustava u bilo kojem trenutku gibanja, a kod kojih nepredviđene<br />
promjene u radnoj okolici mogu prouzročiti neplanirani kontakt robota s okolinom.<br />
– Ako postoji samo upravljanje po položaju, tada zbog realno postojećih pogrešaka u<br />
planiranju željene operacije te zbog konačne točnosti regulacijske petlje položaja,<br />
manipulator može udariti i oštetiti sebe ili elemente radne okoline � može se<br />
zaključiti da je u tim slučajevima sustav suviše krut (tj. nedovoljno podatan) –<br />
rješenje može biti dodavanje regulacijskog kruga sile dodira.<br />
73<br />
74
Podatn<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Pisanje po ploči Okretanje ručice<br />
x<br />
z<br />
y<br />
xc<br />
ωc<br />
zc<br />
yc<br />
Prirodna<br />
ograničenja<br />
vx =<br />
f y<br />
vz =<br />
ω x<br />
ω y<br />
τ z<br />
0<br />
= 0<br />
0<br />
= 0<br />
= 0<br />
= 0<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Umetanje klina u rupu<br />
xc<br />
zc<br />
yc<br />
Prirodna<br />
ograničenja<br />
vx =<br />
v y<br />
f z<br />
ω x<br />
ω y<br />
τ z<br />
0<br />
= 0<br />
= 0<br />
= 0<br />
= 0<br />
= 0<br />
Umjetna<br />
ograničenja<br />
f x<br />
= 0<br />
0<br />
z= v<br />
= 0<br />
= 0<br />
= 0<br />
vy =<br />
v d<br />
τ x<br />
τ y<br />
ω z<br />
Umjetna<br />
ograničenja<br />
f x<br />
v y<br />
f z<br />
τ x<br />
τ y<br />
ω d<br />
= 0<br />
= 0<br />
= 0<br />
= 0<br />
= 0<br />
z= ω<br />
75<br />
76
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira manipulatora<br />
• Senzori sile (ugrađuju se u prihvatnice robotske ruke)<br />
– registriraju kontakt i iznos sile kontakta, nakon čega se ovisno o robotskom<br />
sustavu, donosi odluka o akciji � akcija može biti trenutno zaustavljanje<br />
gibanja robota ili isključenje upravljanja po trajektoriji gibanja i uključenje<br />
upravljanja po sili.<br />
• Podatnost kao mjera mekoće dodira<br />
– Pasivna podatnost - podrazumijeva se svojstvo samog manipulatora (npr.<br />
prihvatnica presvučena gumom),<br />
– Aktivna podatnost - podatnost koja se podržava programskom podrškom, tj.<br />
računalom. Prednosti aktivne podatnosti su mogućnost preprogramiranja i<br />
prilagođavanja raznim primjenama.<br />
• Dualnost upravljanja silom i upravljanja položajem<br />
– Ako su osigurani uvjeti slobodnog gibanja po svim osima � dovoljna je samo<br />
regulacija položaja.<br />
– Ako su uvjeti ograničenja prisutni po svim osima � jedino je moguća regulacija sile.<br />
– Između dviju krajnosti, kombinacijom istovremene regulacije položaja i sile, postiže<br />
se tražena podatnost manipulatora.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira manipulatora<br />
Dva osnovna pristupa upravljanju<br />
sile dodira:<br />
• hibridno upravljanje položajem<br />
i silom dodira<br />
– upravljanje položajem, koje<br />
se odvija u vremenu kada<br />
gibanje nije ograničeno<br />
okolinom i upravljanje silom<br />
kada je gibanje ograničeno<br />
okolinom. Upravljanje<br />
položajem odvija se u<br />
tangencijalnoj ravnini, a<br />
upravljanje silom u ravnini<br />
okomitoj na sučeljenu<br />
površinu<br />
• upravljanje impedancijom<br />
manipulatora.<br />
77<br />
78
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira robotskog mehanizma<br />
s jednim stupnjem slobode gibanja<br />
• Senzor sile dodira - montiran<br />
na mjestu korijena šake<br />
• U režimu slobodnog gibanja<br />
po svim osima, mehanizam je<br />
upravljan po položaju<br />
•Uslučaju dodira s okolinom,<br />
upravljan je po sili dodira<br />
• Zglob pogoni vektorski<br />
upravljan čoperski napajan<br />
elektronički komutirani motor<br />
• Kaskadna regulacija – P<br />
regulator položaja/sile, PI<br />
regulator brzine vrtnje i PWM<br />
regulator struje<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Model dodira manipulatora i radne okoline (1 stupanj<br />
slobode gibanja uz translacijsko gibanje vrha alata)<br />
79<br />
80
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Model dodira manipulatora i radne okoline (1 stupanj<br />
slobode gibanja uz translacijsko gibanje vrha alata)<br />
Profil predmeta, tj. poremećajna veličina<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Matematički model regulacijske petlje sile dodira<br />
• Regulacijska petlja struje<br />
– Fazne struje statora reguliraju se pomoću regulatora s impulsno širinskom<br />
modulacijom faznih napona statora<br />
– Tipične vrijednosti frekvencije prekidanja čopera su u rasponu 5 − 20 kHz,<br />
čime se postiže gotovo trenutačno upravljanje strujom.<br />
– Zbog toga se zatvorena strujna petlja vektorski upravljanog sinkronog<br />
motora može opisati sljedećom prijenosnom funkcijom:<br />
• Regulacijska petlja brzine vrtnje<br />
– PI regulator brzine vrtnje:<br />
81<br />
82
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Blokovska shema regulacijska petlje brzine vrtnje<br />
• Prijenosna funkcija otvorene regulacijske petlje brzine vrtnje (bez regulatora)<br />
parametri svedeni na<br />
osovinu motora<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
s<br />
83<br />
• Zadani su sljedeći pokazatelji kvalitete<br />
odziva regulacijske petlje brzine vrtnje:<br />
– σ m =20%, t m =15 ms<br />
• Sinteza parametara regulatora brzine<br />
vrtnje primjenom metode krivulje mjesta<br />
korijena: K R=15, T I=0.03 s<br />
• Pol u ishodištu i pol 1/T M vrlo su blizu<br />
jedan drugoga pa se uslijed odabranog<br />
mjerila teško razlikuju. Iz slike se može<br />
vidjeti da je utjecaj nule gotovo eliminiran<br />
vrlo bliskim polom, tako da je dinamika<br />
sustava određena dominantnim parom<br />
konjugirano kompleksnih polova. Preostali<br />
pol sustava nalazi se vrlo daleko od<br />
ishodišta i njegovo djelovanje je<br />
zanemarivo. Zanemaren je i utjecaj 1/T zi.<br />
84
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odziv signala povratne veze brzine vrtnje pri<br />
skokovitoj promjeni reference.<br />
• Zadani su sljedeći pokazatelji kvalitete<br />
odziva regulacijske petlje brzine vrtnje:<br />
– σ m =20%, t m =15 ms<br />
• Sinteza parametara regulatora brzine<br />
vrtnje primjenom metode krivulje mjesta<br />
korijena: K R=15, T I=0.03 s<br />
• Pol u ishodištu i pol 1/T M vrlo su blizu<br />
jedan drugoga pa se uslijed odabranog<br />
mjerila teško razlikuju. Iz slike se može<br />
vidjeti da je utjecaj nule gotovo eliminiran<br />
vrlo bliskim polom, tako da je dinamika<br />
sustava određena dominantnim parom<br />
konjugirano kompleksnih polova. Preostali<br />
pol sustava nalazi se vrlo daleko od<br />
ishodišta i njegovo djelovanje je<br />
zanemarivo. Zanemaren je i utjecaj 1/T zi.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Postupak sinteze P regulatora sile dodira<br />
85<br />
86
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Postupak sinteze P regulatora sile dodira<br />
• Blokovska shema otvorene regulacijske petlje sile dodira.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Prijenosna funkcija otvorene regulacijske petlje sile dodira uz K F = 1 (bez regulatora)<br />
• Prijenosna funkcija otvorene regulacijske petlje sile dodira nakon<br />
uvrštavanja prijenosnih funkcija G 1 do G 4:<br />
Astatički sustav tipa IT 6<br />
87<br />
88
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Krivulja mjesta korijena regulacijske petlje sile<br />
dodira (gumena površina).<br />
• Zadani su sljedeći pokazatelji kvalitete<br />
odziva regulacijske petlje sile dodira:<br />
– σ m =25%, t m =50 ms<br />
• Sinteza parametara regulatora sile dodira<br />
primjenom metode krivulje mjesta<br />
korijena: K F=0.03<br />
• Polovi, koji su vrlo udaljeni od ishodišta,<br />
zanemareni su zbog neznatnog djelovanja<br />
na dinamičko ponašanje sustava.<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odziv signala sile dodira pri skokovitoj promjeni<br />
reference (gumena površina).<br />
89<br />
Odziv signala povratne veze sile dodira pri skokovitoj<br />
promjeni reference (gumena površina).<br />
90
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
• Primjena metoda robusnog i adaptivnog upravljanja � adaptacija<br />
koeficijenta pojačanja regulatora sile K F<br />
Nestabilan odziv sile za isti K F =0.03<br />
pri promjeni površine (aluminij).<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira preko upravljanja<br />
impedancijom manipulatora<br />
• Impedancija je fizički sustav koji je u stanju prihvatiti ulaznu pobudu<br />
(nalog za gibanje) i odgovoriti s ostvarenim izlaznim učinkom (npr.<br />
silom).<br />
• Robotski manipulator ponaša se kao impedancija, a shodno tome,<br />
radna okolica djeluje kao admitancija, tj. posjeduje suprotna<br />
svojstva.<br />
• Impedancija robotskog manipulatora sastoji se od otpornosti<br />
(trenja), kapacitivnosti (elastičnosti) i inercije (mase), koji su<br />
određeni mehaničkom konstrukcijom manipulatora.<br />
• Parametri impedancije definiraju pasivnu podatnost manipulatora<br />
� npr. umetanjem gumenih dijelova u konstrukciju manipulatora<br />
mijenjat će se pasivna podatnost.<br />
• Upravljani robotski manipulator ponaša se kao aktivna impedancija<br />
� Kod upravljanja impedancijom robota parametri regulatora sile<br />
dodira podešavaju se tako da se uz zadanu konstrukciju robotskog<br />
manipulatora dobije željena podatnost.<br />
91<br />
92
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira preko upravljanja<br />
impedancijom manipulatora<br />
Svrha upravljanja impedancijom je dobivanje<br />
dinamičkog odnosa između pozicije i sile, a ne<br />
upravljanja svakom od tih varijabli posebno:<br />
d<br />
Fa = K s ( xr<br />
− x)<br />
+ bs<br />
( xr<br />
− x)<br />
dt<br />
F a - sila kojom djeluje vrh alata na površinu predmeta,<br />
x - stvarna pozicija vrha alata,<br />
x r - željena pozicija vrha alata,<br />
K s - koeficijent krutosti alata,<br />
b s - koeficijent viskoznosti alata,<br />
x p - površina predmeta,<br />
K e - koeficijent krutosti predmeta.<br />
K s i b s se postavljaju programski<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira preko upravljanja<br />
impedancijom manipulatora<br />
Fp = K e ( x p − x)<br />
F p - sila koju predmet osjeća,<br />
K e - koeficijent krutosti predmeta (određuje se<br />
eksperimentalno),<br />
x p - površina predmeta,<br />
x - pozicija do koje je vrh alata prodro kroz predmet<br />
Položaj vrha alata koji prodire u predmet:<br />
x = x p + Δxa<br />
+ Δx<br />
Δxa<br />
Δx<br />
p<br />
p<br />
- promjena položaja alata prouzrokovana<br />
djelovanjem sile robota na predmet<br />
- promjena profila predmeta prouzrokovana<br />
djelovanjem sile robota na predmet<br />
93<br />
94
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Upravljanje silom dodira preko upravljanja<br />
impedancijom manipulatora<br />
Δx<br />
a<br />
F<br />
=<br />
K s<br />
F<br />
Δx<br />
p =<br />
K e<br />
U stacionarnom stanju (kada je zadana sila<br />
postignuta) pozicija vrha alata se računa na<br />
sljedeći način:<br />
K s xr<br />
+ K e x p<br />
xs<br />
=<br />
K + K<br />
Fr<br />
x =<br />
K e<br />
+<br />
s<br />
e<br />
Mi robotu zadajemo položaj vrha alata x r, a ona<br />
se dobiva na sljedeći način:<br />
K e K e<br />
K e + K s K e<br />
= x − x + x xr = x − x p<br />
K K<br />
K K<br />
xr p<br />
s s<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
U stacionarnom stanju sila je stalna i željeni položaj vrha alata treba<br />
odrediti prema sljedećem izrazu:<br />
K + K<br />
x = ⋅ F + x = K F + x<br />
e s<br />
r<br />
Ke⋅Ks r p Σ r p<br />
Osnovni princip upravljanja impedancijom se svodi na pretpostavku<br />
da se površina okoline u trenutku dodira manipulatora elastično<br />
deformira � pod prepostavkom da je na raspolaganju senzor<br />
položaja dovoljne rezolucije može se mjeriti «prodor» u okolinu �<br />
na temelju tog pomaka uz poznavanje krutosti okoline K e zna se<br />
iznos sile kojom manipulator djeluje na okolinu.<br />
Upravljački zakon upravljanja impedancijom robota ima oblik:<br />
F Ks x Bs dt<br />
s<br />
x<br />
s<br />
p<br />
95<br />
0 x2≤xp dΔx ⎧⎪<br />
= Δ + pri čemu vrijedi: F = ⎨<br />
⎪⎩<br />
K ( x − x ) x > x<br />
e 2 p 2 p<br />
96
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Primjer: Upravljanje silom dodira troosnog planarnog<br />
manipulatora preko upravljanja impedancijom<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Primjer: Upravljanje silom dodira troosnog planarnog<br />
manipulatora preko upravljanja impedancijom<br />
97<br />
98
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Primjer: Upravljanje silom dodira troosnog planarnog<br />
manipulatora preko upravljanja impedancijom<br />
Zakon upravljanja impedancijom robota ima oblik:<br />
dΔd F = KsΔ d + Bs dt<br />
F<br />
pri čemu vrijedi:<br />
⎧⎪ 0 d ≤ dp<br />
= ⎨<br />
⎪⎩<br />
Ke( d − dp) d > dp<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
99<br />
100
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Shema u Matlabu sustava upravljanja pozicijom, brzinom vrtnje i strujom<br />
armature za sva tri motora u zglobovima planarnog manipulatora<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odziv sile dodira pri upravljanju impedancijom<br />
bs<br />
K s<br />
K e<br />
= 100<br />
= 10000<br />
101<br />
=<br />
6000N<br />
/ m<br />
102
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odziv sile dodira pri upravljanju impedancijom – veći K s<br />
bs<br />
= 100<br />
K = 100000<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odziv sile dodira pri upravljanju impedancijom – veći b s<br />
s<br />
K e<br />
b s =<br />
K s<br />
= 6000N<br />
/ m<br />
500<br />
= 10000<br />
103<br />
K =<br />
6000 N / m<br />
e<br />
104
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
Odziv sile dodira pri upravljanju impedancijom – tvrđi predmet<br />
bs<br />
K s<br />
= 100<br />
<strong>Sveučilište</strong> u <strong>Zagrebu</strong>, Fakultet elektrotehnike i računarstva Osnove robotike<br />
Zavod za automatiku i računalno inženjerstvo Predavanja<br />
= 10000<br />
K = 12000 N / m<br />
Odziv sile dodira pri upravljanju impedancijom – tvrđi predmet pomaknut za 1 mm<br />
e<br />
bs<br />
K s<br />
= 100<br />
= 10000<br />
105<br />
K =<br />
12000 N / m<br />
e<br />
106