ControlLab základy PID regulace

Tohle je statická kopie kapitoly pro vyhledávače. Interaktivní verze má animované obrázky, kontrolní otázky a tlačítka, která příklad načtou do kalkulačky.

Ladění: pravidla a jejich meze

Ladicí pravidlo je tabulka, kterou někdo nafitoval na rodinu soustav. Není to návrh. Je to odhad, který tě dostane do správné dekády — a to je na reálné soustavě většina práce.

Ziegler–Nichols metodou mezní citlivosti (1942)

Nejstarší a pořád nejužitečnější, protože nepotřebuje model. Postup je experiment:

  1. Vypni I i D, nech jen P.
  2. Zvyšuj Kp, dokud se smyčka nerozkmitá ustáleně — kmity se ani nezvětšují, ani nezmenšují.
  3. To zesílení je Ku, perioda kmitů Tu.
  4. Odečti nastavení z tabulky.
0.010200.01.02.0čas [s]y0.6·KᵤKᵤ1.35·Kᵤ
Tři zesílení na téže soustavě. Pod Ku kmity odeznívají, nad Ku rostou, a přesně při Ku se drží. Ta hranice není náhoda: je to přesně frekvence, na které má soustava fázi −180°.
TypKpTiTd
P0,50·Ku
PI0,45·KuTu/1,2
PID (klasický)0,60·KuTu/2Tu/8
malý překmit0,33·KuTu/2Tu/3
bez překmitu0,20·KuTu/2Tu/3

Ta zkouška je nebezpečná. Uvádíš stroj záměrně na hranici stability. Na ohřevu vody nevadí, na robotickém rameni nebo v tlakové nádobě je to nepřijatelné. Bezpečnější variantou je relé test: místo zesilování se dočasně nasadí dvoupolohový regulátor, který rozkmitá soustavu s omezenou amplitudou, a Ku se dopočítá.

Ziegler–Nichols z přechodové charakteristiky

Nekmitá se, jen se skočí vstupem a odečte K, T a L. Pravidlo se pak píše normovaně přes a = K·L/T:

a·Kp = 1,2  ·  Ti = 2L  ·  Td = 0,5LZN, přechodová charakteristika, PID

AMIGO — Ziegler-Nichols po šedesáti letech

Åström a Hägglund udělali totéž co ZN, ale s jednou zásadní změnou: přidali požadavek na robustnost (Ms ≤ 1,4). Výsledek je znatelně klidnější a bezpečnější:

Kp = 1K(0,2 + 0,45 TL)  ·  Ti = (0,4L + 0,8T)(L + 0,1T) L  ·  Td = 0,5·L·T0,3L + TAMIGO PID (Åström & Hägglund, rov. 7.7)

SIMC — jedno kolečko, kterým se to dolaďuje

Skogestadovo pravidlo má jediný parametr: τc, tedy jak rychlá má uzavřená smyčka být. Výchozí volba τc = L je „rychlé, ale ještě robustní“; větší hodnota je klidnější.

Kp = TK(τc + L)  ·  Ti = min(T, 4(τc + L))SIMC PI (Skogestad 2003)

To min() je ta část, kterou si všichni odpustí — a bez ní vyjde u pomalé soustavy integrační čas tak dlouhý, že smyčka nikdy nic nesrovná.

0.05.0100.01.02.0wčas [s]yZN classicZN stepAMIGO PIDSIMC PI
ZN classic: Kp = 2,60, Ti = 0,75 s, Td = 0,19 s  ·  překmit 46 %  ·  Ms 2,69
Přepínej pravidla a měň obtížnost soustavy jezdcem τ. Zvýrazněná čára je vybrané pravidlo, ostatní zůstávají čárkovaně pro srovnání. Sleduj, jak Ziegler–Nichols překmitává skoro pořád, a jak se to zhoršuje s rostoucím τ — a jak AMIGO drží Ms kolem 1,4 napříč celým rozsahem.

Co v tomhle nástroji schválně není

Cohen–Coon. Je to známé pravidlo, ale zdroje, které se podařilo dohledat, si navzájem odporují v koeficientech — a ladicí pravidlo, které nejde ověřit, je horší než žádné. Radši ho tady nenajdeš, než aby ti nástroj vydával za ověřené číslo něco, co jsem opsal z obrázku. (Stejný důvod, proč GearLab nepočítá ohyb v patě zubu.)

Ruční doladění, když pravidlo nestačí

  1. Začni pravidlem — dostane tě do správné dekády.
  2. Sleduj tvar kmitání, ne jeho velikost: rychlé = Kp, pomalé = Ti, chlupatá akční veličina = D.
  3. Měň jeden parametr, po dvojnásobcích a polovinách. Změna o 10 % je nepozorovatelná a jen tě zmate.
  4. Zkoušej odezvu na poruchu, ne jen na skok žádané hodnoty. Většina smyček tráví život odoláváním poruchám a skok dostane jednou za směnu.

Vzorce v této kapitole

Ku — mezní zesílení
Ku = 1/|P(jω₁₈₀)| [—] Ziegler & Nichols 1942
Tu — perioda mezního kmitání
Tu = 2π/ω₁₈₀ [s] Ziegler & Nichols 1942
ZN — Ziegler–Nichols PID z mezního kmitání
Kp = 0,6·Ku, Ti = Tu/2, Td = Tu/8 [—] Trans. ASME 64 (1942) 759
ZN-a — Ziegler–Nichols z přechodové charakteristiky
a·Kp = 1,2, Ti = 2L, Td = 0,5L, a = K·L/T [—] ZN 1942, ve tvaru dle Åströma
AMIGO — AMIGO PID
Kp = (1/K)(0,2 + 0,45·T/L); Ti = L(0,4L+0,8T)/(L+0,1T); Td = 0,5LT/(0,3L+T) [—] Åström & Hägglund, Advanced PID Control, eq. 7.7
SIMC — SIMC PI
Kp = T/(K(τc+L)), Ti = min(T, 4(τc+L)), τc = L [—] Skogestad, J. Process Control 13 (2003) 291