Reguleringsteknikk: gratis øving, teori og oppgaver

1. Laplace og overføringsfunksjoner

Hva handler det om?

I reguleringsteknikk beskriver vi hvordan et system (en motor, en temperaturprosess, en robotarm) reagerer på et pådrag ved å transformere differensialligningen til Laplace-planet. Der blir differensiallikninger til algebra, og systemet kan beskrives med én kompakt brøk: overføringsfunksjonen $G(s)$. Polene og nullpunktene til $G(s)$ forteller nesten alt du trenger å vite om systemets oppførsel – om det er stabilt, hvor fort det reagerer, og om det svinger.

Målet med enheten er å kunne lese en overføringsfunksjon og si noe kvalitativt og kvantitativt om sprangresponsen uten å løse differensialligningen i tidsplanet hver gang.

Begreper og formler

• Overføringsfunksjon: $G(s) = Y(s)/U(s)$, forholdet mellom Laplace-transformert utgang og inngang med null initialbetingelser.
• Poler: verdiene av $s$ som gjør nevneren (den karakteristiske ligningen) null. Nullpunkter: verdiene som gjør telleren null.
• Stabilitet: systemet er stabilt hvis og bare hvis alle polene har negativ realdel (ligger i venstre halvplan).
• Førsteordens system: $G(s) = \dfrac{K}{\tau s+1}$, med pol i $s=-1/\tau$. Sluttverdi ved enhetssprang: $G(0)=K$. Etter én tidskonstant $\tau$ er responsen på 63 % av sluttverdien.
• Andreordens standardform: $G(s) = \dfrac{\omega_n^2}{s^2+2\zeta\omega_ns+\omega_n^2}$, med naturlig frekvens $\omega_n$ og dempingsforhold $\zeta$.
• Oversving (for $0<\zeta<1$): $M_p = e^{-\pi\zeta/\sqrt{1-\zeta^2}}$. Topptid $t_p = \pi/\omega_d$ og innsvingningstid (2 %) $t_s \approx 4/(\zeta\omega_n)$, der $\omega_d=\omega_n\sqrt{1-\zeta^2}$ er den dempede (svingende) frekvensen.
• Dødtid $T$ gir et ekstra ledd $e^{-sT}$ i overføringsfunksjonen: en ren forsinkelse før utgangen reagerer.

Slik løser du oppgavene

1. Identifiser om systemet er første- eller andreordens, og les av $K$, $\tau$ (eller $\omega_n$, $\zeta$) direkte fra formen på $G(s)$.
2. Finn polene ved å sette nevneren lik null, og sjekk fortegnet på realdelen for stabilitet.
3. Finn sluttverdien ved et enhetssprang med $G(0)$ (sett $s=0$).
4. Bruk standardformlene for oversving, topptid og innsvingningstid for andreordens systemer i stedet for å løse differensialligningen.
5. Sjekk svaret opp mot fysikken: et system uten integrator og med positiv $K$ skal gå mot en endelig, positiv sluttverdi.

Eksempel

Et system har overføringsfunksjonen $G(s) = \dfrac{36}{s^2+3,6s+36}$. Finn dempingsforhold, sluttverdi ved enhetssprang og oversving.

1. Sammenlign med standardformen: $\omega_n^2=36$, så $\omega_n=6$ rad/s.
2. $2\zeta\omega_n=3,6$ gir $\zeta = 3,6/(2\cdot6) = 0,3$.
3. Sluttverdi: $G(0)=36/36=1$.
4. Oversving: $M_p=e^{-\pi\cdot0,3/\sqrt{1-0,3^2}}\approx e^{-0,988}\approx0,372$, altså ca. 37 %.

Svar: $\zeta=0,3$, sluttverdi $=1$, oversving $\approx37$ %.

Vanlige feil

• Å tro at nullpunkter bestemmer stabiliteten – det er polene som avgjør stabilitet, ikke nullpunktene.
• Å glemme å sette $s=0$ riktig når man finner sluttverdien (gjelder bare for stabile systemer uten integrator).
• Å blande sammen $\omega_n$ og den dempede frekvensen $\omega_d=\omega_n\sqrt{1-\zeta^2}$ – de er ikke det samme når $\zeta>0$.
• Å tro at dødtid $e^{-sT}$ er det samme som en pol – den endrer fasen, ikke antall poler i en rasjonal tilnærming.

Begreper i denne delen

Øv på laplace og overføringsfunksjoner i appen →

2. PID-regulering

Hva handler det om?

En PID-regulator er den mest brukte reguleringsalgoritmen i industrien: den sammenligner en målt verdi med en referanse, og beregner et pådrag ut fra avviket. De tre leddene – proporsjonal (P), integral (I) og derivat (D) – har hver sin jobb, og god regulering handler om å balansere dem mot hverandre: rask respons, lite oversving, null stasjonært avvik og robusthet mot støy.

Enheten handler om å forstå hva hvert ledd gjør, hvordan de kombineres i en tilbakekoblet sløyfe, og hvordan man kan finne fornuftige startverdier for forsterkningene.

Begreper og formler

• Ideell PID-regulator: $C(s) = K_p + \dfrac{K_i}{s} + K_ds$, der $K_i=K_p/T_i$ og $K_d=K_pT_d$ ($T_i$ = integrasjonstid, $T_d$ = derivasjonstid).
• P-leddet gir rask respons, men lar et stasjonært avvik stå igjen for de fleste prosesser.
• I-leddet integrerer avviket over tid og fjerner stasjonært avvik, men kan gjøre systemet tregere og mer svingete.
• D-leddet reagerer på endringstakten til avviket og demper oversving, men forsterker målestøy og filtreres derfor ofte.
• Lukket sløyfe med forover-forsterkning $G$ og tilbakekobling $H$: $\dfrac{Y}{R} = \dfrac{G}{1+GH}$ (negativ tilbakekobling).
• Stasjonært avvik ved sprang: $e_\infty = \dfrac{1}{1+K_pG(0)}$ for ren P-regulering; med et I-ledd blir $e_\infty=0$ for et sprang.
• Stasjonært avvik ved rampe (type 1-system, PI-regulering av $G(s)=K/(\tau s+1)$): hastighetskonstanten $K_v = K_pK/T_i$, og $e_\infty = 1/K_v$.
• Ziegler–Nichols (svingemetoden): finn kritisk forsterkning $K_u$ og periode $T_u$ (ren P-regulering helt til vedvarende svingning). PID-tabellen: $K_p=0,6K_u$, $T_i=T_u/2$, $T_d=T_u/8$.
• Integratoroppvikling (windup): I-leddet fortsetter å vokse mens pådraget er i metning, og gir store oversving når avviket snur.

Slik løser du oppgavene

1. Identifiser hvilket ledd (P, I eller D) spørsmålet handler om, og hva slags feil eller virkning det gir.
2. For stasjonært avvik: sjekk om sløyfen har et I-ledd (type 1) – da er avviket null for et sprang, men kan være endelig for en rampe.
3. For lukket-sløyfe-uttrykk: sett opp $G/(1+GH)$ og forenkle før du finner poler eller sluttverdi.
4. For Ziegler–Nichols: sett $K_u$ og $T_u$ rett inn i tabellformlene – ikke bland sammen $K_p$, $T_i$ og $T_d$.
5. Sjekk om svaret gir mening: økt $K_p$ skal gi raskere, men mer urolig respons; et I-ledd skal fjerne stasjonært avvik ved sprang.

Eksempel

En prosess $G(s) = \dfrac{2}{s+1}$ reguleres med en PI-regulator med $K_p=3$ og $T_i=0,5$ s. Referansen er en rampe med helning 1 (enhet/s). Finn det stasjonære avviket.

1. Hastighetskonstanten: $K_v = K_pK/T_i = 3\cdot2/0,5 = 12$.
2. Stasjonært avvik: $e_\infty = 1/K_v = 1/12 \approx 0,0833$.

Svar: $e_\infty \approx 0,083$ (enheter).

Vanlige feil

• Å tro at et I-ledd fjerner stasjonært avvik for *alle* typer referanser – det fjerner avviket for et sprang, men ikke nødvendigvis for en rampe.
• Å øke $K_p$ mye for å bli kvitt stasjonært avvik i stedet for å legge til et I-ledd – det gir mer oversving og kan gjøre systemet ustabilt.
• Å forveksle $T_i$ (integrasjonstid, stor $T_i$ = svak I-virkning) med $K_i$ (der stor $K_i$ = sterk I-virkning) – de er hverandres inverse skalert med $K_p$.
• Å glemme anti-windup i en regulator med pådragsmetning, slik at I-leddet "låser seg" og gir store oversving.

Begreper i denne delen

Øv på pid-regulering i appen →

3. Frekvensanalyse

Hva handler det om?

Bodediagrammet viser hvordan et system svarer på sinusformede innsignaler ved forskjellige frekvenser: hvor mye amplituden forsterkes eller dempes (i desibel), og hvor mye fasen forskyves. Denne frekvensresponsen brukes til å vurdere hvor stabilt et lukket reguleringssystem blir – uten å måtte finne polene eksplisitt – og til å forstå hvor mye du kan skru opp forsterkningen før systemet begynner å svinge ukontrollert.

Fase- og forsterkningsmargin er de viktigste "sikkerhetsmarginene" en regulatoringeniør leser av et Bodediagram før en regulator settes i drift.

Begreper og formler

• Forsterkning i desibel: $|G|_{dB} = 20\log_{10}|G|$. Omvendt: $|G| = 10^{|G|_{dB}/20}$.
• En enkelt pol ($1/(\tau s+1)$) gir asymptotisk $-20$ dB/dekade over knekkfrekvensen $\omega_c = 1/\tau$; et enkelt nullpunkt gir $+20$ dB/dekade.
• Fasebidrag fra en enkelt pol går gradvis fra $0^\circ$ til $-90^\circ$ rundt knekkfrekvensen; fra et nullpunkt fra $0^\circ$ til $+90^\circ$.
• En ren dødtid $e^{-sT}$ endrer ikke amplituden, men gir en faseforsinkelse $\varphi = -\omega T$ (radianer) som vokser lineært med frekvensen.
• Gjennomskjæringsfrekvens $\omega_c$: der $|L(j\omega)| = 1$ (0 dB), med $L$ = sløyfeforsterkningen $= GH$.
• Fasemargin: $PM = 180^\circ + \angle L(j\omega_c)$ ved gjennomskjæringsfrekvensen. Vanlig mål: 45–60°.
• Forsterkningsmargin: hvor mange dB $|L|$ kan økes med før systemet blir ustabilt, avlest der fasen er $-180^\circ$.
• Båndbredde: frekvensen der amplituden i det lukkede systemet har falt 3 dB fra lavfrekvensverdien.

Slik løser du oppgavene

1. Finn knekkfrekvensene til hver pol/nullpunkt ($\omega=1/\tau$) og resonnér deg fram til asymptotene.
2. Legg sammen amplitude- og fasebidrag fra hvert ledd (poler, nullpunkter, eventuell dødtid) for å få total amplitude og fase.
3. Finn gjennomskjæringsfrekvensen $\omega_c$ der $|L|=1$, og les av fasen der for å regne ut fasemarginen.
4. For dødtid: regn fasebidraget som $-\omega T$ i radianer, og gjør om til grader ved å multiplisere med $180/\pi$.
5. Sjekk at svaret er rimelig: en stabil sløyfe bør ha positiv fase- og forsterkningsmargin.

Eksempel

En sløyfe har $L(s) = \dfrac{4}{s}e^{-sT}$ med dødtid $T = 0,05$ s. Finn fasemarginen.

1. Gjennomskjæringsfrekvens: $|L(j\omega)| = 4/\omega = 1$ gir $\omega_c = 4$ rad/s.
2. Fase fra integratoren: $-90^\circ$.
3. Fase fra dødtiden: $\varphi_T = -\omega_cT\cdot(180/\pi) = -4\cdot0,05\cdot(180/\pi) \approx -11,46^\circ$.
4. Total fase: $-90^\circ - 11,46^\circ = -101,46^\circ$.
5. Fasemargin: $PM = 180^\circ - 101,46^\circ \approx 78,5^\circ$.

Svar: $PM \approx 78,5^\circ$ (god margin).

Vanlige feil

• Å tro at dødtid endrer amplituden – den endrer bare fasen, aldri $|G|$.
• Å lese av forsterkningsmarginen ved gjennomskjæringsfrekvensen i stedet for der fasen er $-180^\circ$ (og omvendt for fasemarginen).
• Å glemme å regne om fasebidraget fra dødtid fra radianer til grader.
• Å tro at høyere båndbredde alltid er bedre – det gir raskere respons, men også mer følsomhet for støy.

Begreper i denne delen

Øv på frekvensanalyse i appen →

4. Stegrespons og førsteordens systemer

Hva handler det om?

Mange prosesser reagerer på samme måte når du endrer inngangen brått: en vannkoker som slås på, en tank som fylles, en motor som får spenning. Utgangen vokser raskt først og flater så ut mot en ny verdi. Dette er et førsteordens system, og det beskrives fullt ut av to tall: forsterkningen $K$ og tidskonstanten $\tau$.

Begreper og formler

• Overføringsfunksjon:

• Stegrespons: når inngangen hopper med $A$ ved $t = 0$, blir utgangen

• Sluttverdi: $y(\infty) = KA$. Forsterkningen $K$ er endringen i utgang delt på endringen i inngang når alt har stabilisert seg.
• Tidskonstanten $\tau$ er tiden det tar å nå 63,2 % av endringen. Etter $3\tau$ er du på 95 %, etter $4\tau$ på 98 % og etter $5\tau$ på 99,3 %.
• Stigetid (fra 10 % til 90 %): $t_r = \tau\ln 9 \approx 2{,}2\,\tau$.
• RC-krets: $\tau = RC$. En termisk prosess: $\tau$ = termisk kapasitans ganger termisk motstand.
• Tid for å nå en andel $p$ av sluttverdien: $t = -\tau\ln(1 - p)$.

Slik løser du oppgavene

1. Les av $K$ og $\tau$ fra $G(s)$. Pass på at nevneren står på formen $\tau s + 1$ (del på konstantleddet om nødvendig).
2. Sluttverdi: $KA$. Verdi ved tiden $t$: $KA(1 - e^{-t/\tau})$.
3. Tid til en viss andel: $t = -\tau\ln(1-p)$.

Eksempel

$G(s) = \dfrac{5}{2s + 1}$, og inngangen hopper fra 0 til 3.

1. $K = 5$ og $\tau = 2$ s.
2. Sluttverdien er $5\cdot 3 = 15$.
3. Etter 2 s er utgangen $15\,(1 - e^{-1}) \approx 9{,}48$, altså 63,2 % av veien.
4. Innsvingningstid til 98 %: $4\tau = 8$ s.

Vanlige feil

• Å lese av feil tidskonstant når nevneren ikke er normalisert: $\dfrac{10}{4s + 2} = \dfrac{5}{2s + 1}$ gir $\tau = 2$, ikke 4.
• Å glemme å gange med steghøyden $A$.
• Å blande 63 % (én tidskonstant) og 98 % (fire tidskonstanter).

Begreper i denne delen

Øv på stegrespons og førsteordens systemer i appen →