01.03.2014 Aufrufe

Formelsammlung (Stand: 26.11.2013) - EAL Lehrstuhl für ...

Formelsammlung (Stand: 26.11.2013) - EAL Lehrstuhl für ...

Formelsammlung (Stand: 26.11.2013) - EAL Lehrstuhl für ...

MEHR ANZEIGEN
WENIGER ANZEIGEN

Erfolgreiche ePaper selbst erstellen

Machen Sie aus Ihren PDF Publikationen ein blätterbares Flipbook mit unserer einzigartigen Google optimierten e-Paper Software.

Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

1 Trigonometrische Formeln (siehe [1])<br />

Im Folgenden gelte x,y ∈ R (bei entsprechender Einschränkung des Bereiches falls notwendig).<br />

sin(x±y) = sin(x)cos(y)±cos(x)sin(y) (1)<br />

cos(x±y) = cos(x)cos(y)∓sin(x)sin(y) (2)<br />

tan(x±y) = tan(x)±tan(y)<br />

1∓tan(x)tan(y) = sin(x±y)<br />

cos(x±y)<br />

(3)<br />

sin(x)sin(y) = 1 ( )<br />

cos(x−y)−cos(x+y) (4)<br />

2<br />

cos(x)cos(y) = 1 ( )<br />

cos(x−y)+cos(x+y) (5)<br />

2<br />

sin(x)cos(y) = 1 ( )<br />

sin(x−y)+sin(x+y) (6)<br />

2 ( ) ( ) x+y x−y<br />

sin(x)+sin(y) = 2sin sin<br />

(7)<br />

2 2<br />

( ) ( ) x+y x−y<br />

sin(x)−sin(y) = 2cos sin<br />

(8)<br />

2 2<br />

( ) ( ) x+y x−y<br />

cos(x)+cos(y) = 2cos cos<br />

(9)<br />

2 2<br />

( ) ( ) y +x y −x<br />

cos(x)−cos(y) = 2sin sin<br />

(10)<br />

2 2<br />

⎧<br />

( 1<br />

⎨<br />

π<br />

−arctan(x) ,x > 0<br />

arctan = 2<br />

x)<br />

⎩− 2 −arctan(x) ,x < 0. (11)<br />

arctan(−x) = −arctan(x) (12)<br />

Folgende Funktionswerte arctan(x) ergeben sich <strong>für</strong> ausgewählte Argumente x:<br />

x ±∞ ± √ 3 ±1 ± 1 √<br />

3<br />

0<br />

arctan(x) ± π 2<br />

± π 3<br />

± π 4<br />

± π 6<br />

0<br />

Tabelle 1: Ausgewählte Argumente <strong>für</strong> und Funktionswerte von arctan : R → R.<br />

x [ ◦ ] 0 [0 ◦ ]<br />

sin(x) 0<br />

cos(x) 1<br />

tan(x) 0<br />

π π<br />

6 [30◦ ]<br />

4 [45◦ ]<br />

√<br />

1 2<br />

√2<br />

3<br />

√2<br />

2<br />

2 2<br />

1<br />

√<br />

3<br />

1<br />

π π 2π 3π<br />

3 [60◦ ]<br />

2 [90◦ ]<br />

3 [120◦ ]<br />

4 [135◦ ] π [180 ◦ 3π<br />

]<br />

2 [270◦ ]<br />

√ √ √<br />

3 3 2<br />

1<br />

0 −1<br />

2 2 2<br />

1<br />

0 − 1 √<br />

2<br />

− −1 0<br />

2 2 2<br />

√ √<br />

3 ±∞ − 3 −1 0 ±∞<br />

Tabelle 2: Ausgewählte Argumente <strong>für</strong> und Funktionswerte von cos, sin, tan : R → R.<br />

— Seite 1/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

2 Energieeinheiten und Umrechnungsfaktoren<br />

Energieträger Energiegehalt Anmerkung<br />

1 [kg] Steinkohle 8.14 [kWh] –<br />

1 [kg] Rohöl 11.63 [kWh] Benzin: 8.7 [kWh/Liter]; Diesel: 9.8 [kWh/Liter]<br />

1 [m 3 ] Erdgas 8.82 [kWh] –<br />

1 [kg] Holz 4.3 [kWh] (bei 15% Feuchte)<br />

Tabelle 3: Umrechnungsfaktoren verschiedener Energieträger (siehe [2, Tab. 1.2])<br />

[kJ] [kcal] [kWh] [kg] SKE [kg] RÖE [m 3 Erdgas]<br />

1[kJ] = 1000[Ws] 1 0.2388 2.78·10 −4 3.4·10 −5 2.4·10 −5 3.2·10 −5<br />

1[kcal] 4.1868 1 1.163·10 −3 1.43·10 −4 1·10 −4 1.3·10 −4<br />

1[kWh] 3600 860 1 0.123 0.086 0.113<br />

1[kg] SKE 29308 7000 8.14 1 0.7 0.923<br />

1[kg] RÖE 41868 10000 11.63 1.428 1 1.319<br />

1[m 3 Erdgas] 31736 7580 8.816 1.083 0.758 1<br />

Tabelle 4: Umrechnungsfaktoren zwischen verschiedenen Energieeinheiten (siehe [3, Tab. 1.1]) mit den<br />

Abkürzungen [kJ]: Kilojoule, [Ws]: Wattsekunde, [kcal]: Kilokalorie, [kWh]: Kilowattstunde,<br />

[SKE]: Steinkohleeinheit, RÖE: Rohöleinheit und [ m 3] : Kubikmeter (Volumen).<br />

Vorsatz Symbol Wert Vorsatz Symbol Wert<br />

Milli m 10 −3 (Tausendstel) Kilo k 10 3 (Tausend)<br />

Mikro µ 10 −6 (Millionstel) Mega M 10 6 (Million)<br />

Nano n 10 −9 (Milliardstel) Giga G 10 9 (Milliarde)<br />

Piko p 10 −12 (Billionstel) Tera T 10 12 (Billion)<br />

Femto f 10 −15 (Billiardstel) Peta P 10 15 (Billiarde)<br />

Atto a 10 −18 (Trillionstel) Exa E 10 18 (Trillion)<br />

Tabelle 5: Vorsätze, Symbole und Faktoren<br />

3 Wechsel- und Drehstromsysteme (siehe [4])<br />

Im Folgenden sei t 0 ≥ 0 ein beliebiger Zeitpunkt und x: R ≥0 → R ein periodisches Signal mit<br />

Periodendauer T > 0 [s] (d.h. x(t) = x(t+T) <strong>für</strong> alle t ≥ 0) und Amplitude ˆx > 0.<br />

• Gleichwert oder arithmetischer Mittelwert:<br />

X := 1 T<br />

t∫<br />

0 +T<br />

t 0<br />

x(τ)dτ = 1 T<br />

∫ T<br />

0<br />

x(τ)dτ (13)<br />

— Seite 2/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

• Gleichrichtwert:<br />

X DC := 1 T<br />

t∫<br />

0 +T<br />

|x(τ)|dτ = 1 T<br />

∫ T<br />

|x(τ)|dτ (14)<br />

(<strong>für</strong> sinus- oder cosinusförmige Signale gilt X DC =<br />

t 0<br />

0<br />

2<br />

πˆx ≈ 0.637ˆx)<br />

• Effektivwert (engl. root-mean-square/RMS value) oder quadratischer Mittelwert:<br />

X eff := √ 1 t∫<br />

0 +T<br />

√ 1 ∫ T<br />

x(τ)<br />

T T<br />

2 dτ (15)<br />

t 0<br />

x(τ) 2 dτ =<br />

(<strong>für</strong> sinus- oder cosinusförmige Signale gilt X eff = ˆx √<br />

2<br />

≈ 0.707ˆx)<br />

• Momentanleistung im Wechselstromsystem mit Stromi(t) [A] und Spannungu(t) [V]<br />

0<br />

∀t ≥ 0: p ∼ (t) := u(t)i(t) [W] (16)<br />

• Momentanleistung im Drehstromsystem mit Strangströmeni abc (t) = (i a (t), i b (t), i c (t)) ⊤ [A] 3<br />

und Strangspannung u abc (t) = (u a (t), u b (t), u c (t)) ⊤ [V] 3 (jeweils vektoriell zusammengefasst;<br />

drei Stränge mit den Bezeichnungen a,b,c)<br />

∀t ≥ 0: p 3∼ (t) := u abc (t) ⊤ i abc (t) = u a (t)i a (t)+u b (t)i b (t)+u c (t)i c (t)[W] (17)<br />

• Mittlere Leistung <strong>für</strong> periodische Ströme und Spannungen mit Periodendauer T > 0<br />

P ∼ := 1 T<br />

∫ T<br />

0<br />

p ∼ (τ)dτ bzw. P 3∼ := 1 T<br />

∫ T<br />

0<br />

p 3∼ (τ)dτ (18)<br />

• Scheinleistung <strong>für</strong> sinusförmige Ströme und Spannungen mit Periodendauer T > 0<br />

S ∼ := U eff I eff [VA] bzw. S 3∼ := 3U eff I eff [VA] (19)<br />

die Effektivwerte der Stranggrößen sind (Sternschal-<br />

wobei U eff = U eff und I eff = I eff<br />

tung).<br />

• Wirkleistung <strong>für</strong> sinusförmige Ströme und Spannungen mit Periodendauer T > 0<br />

P ∼ := S ∼ cos(ϕ) = U eff I eff cos(ϕ) [W] bzw.<br />

P 3∼ := S 3∼ cos(ϕ) = 3U eff I eff cos(ϕ) [W]<br />

(20)<br />

• Blindleistung <strong>für</strong> sinusförmige Ströme und Spannungen mit Periodendauer T > 0<br />

Q ∼ := S ∼ sin(ϕ) = U eff I eff sin(ϕ) [var] bzw.<br />

Q 3∼ := S 3∼ sin(ϕ) = 3U eff I eff sin(ϕ) [var]<br />

(21)<br />

— Seite 3/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

4 Grundlagen linearer Regelungstechnik (siehe [1, 5, 6, 7])<br />

• Lösungsformel <strong>für</strong> ax 2 +bx+c = 0<br />

x 1,2 = −b (<br />

1±<br />

2a<br />

√<br />

1−4 ac<br />

b 2 )<br />

<strong>für</strong> a,b,c ∈ R (22)<br />

• Rechenregeln <strong>für</strong> Logarithmus zur Basis x<br />

log x (a·b c ) = log x (a)+clog x (b) <strong>für</strong> a,b,c,x ∈ R (23)<br />

• Komplexe Rechnung<br />

s = σ +j ·ω = |s|exp(j∠s) ∈ C <strong>für</strong> σ,ω ∈ R (24)<br />

wobei s ∗ := σ −j ·ω, |s| = √ s ∗ s = √ R{s} 2 +I{s} 2 und tan∠s = I{s}<br />

R{s} . Dann<br />

∠s = arctan<br />

( ) I{s}<br />

+<br />

R{s}<br />

⎧<br />

⎪⎨ 0 ,I{s} ≥ 0∧R{s} ≥ 0<br />

π ,(I{s} < 0∨I{s} > 0)∧R{s} ≤ 0<br />

⎪⎩<br />

2π ,I{s} ≤ 0∧R{s} ≥ 0.<br />

(25)<br />

• Laplace-Transformation x(s) = ∫ ∞<br />

0<br />

x(t)exp(−st)dt (oder kurz x(t) ❞ x(s))<br />

Für a,b ∈ R und a ≠ b gilt:<br />

ẋ(t) ❞ sx(s)−x(0+) ∗ (26)<br />

ẍ(t) ❞ s 2 x(s)−sx(0+)−ẋ(0+) ∗ (27)<br />

{<br />

1 ,t ≥ T<br />

σ(t−T) = ❞ e −sT<br />

(28)<br />

0 ,t < T s<br />

t n ,n ∈ N ❞ n!<br />

0 (29)<br />

s n+1<br />

1<br />

a−b (exp(−bt)−exp(−at)) ❞ 1<br />

(30)<br />

(s+a)(s+b)<br />

exp(−at)cos(bt) ❞ <br />

1<br />

b exp(−at)sin(bt) ❞ 1<br />

s+a<br />

(s+a) 2 +b 2 (31)<br />

(32)<br />

(s+a) 2 +b 2<br />

1 ❞ 1<br />

(<br />

)<br />

exp(−at)−exp(−bt)+(a−b)texp(−bt) (33)<br />

(s+a)(s+b) 2 (a−b) 2<br />

• Faltungsregel <strong>für</strong> Impulsantworten f(t) ❞ F(s) und g(t) ❞ G(s)<br />

f(t)∗g(t) :=<br />

∫ t<br />

0<br />

f(τ)g(t−τ)dτ ❞ F(s)G(s) (34)<br />

— Seite 4/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

• Additivität von Betrag (in [dB]) und Phase im Bode-Diagramm<br />

F(jω) = F 1 (jω)·...·F n (jω)<br />

= |F 1 (jω)|exp(j∠F 1 (jω))·...·|F n (jω)|exp(j∠F n (jω)) (35)<br />

|F(jω)| dB = 20log(|F 1 (jω)|·...·|F n (jω)|)<br />

= 20log(|F 1 (jω)|)+...+20log(|F n (jω)|)<br />

= |F 1 (jω)| dB +...+|F n (jω)| dB (36)<br />

∠F(jω) = ∠F 1 (jω)+...+∠F n (jω) (37)<br />

• <strong>Stand</strong>ardreglerstrukturen<br />

P-Regler:<br />

u(t) = V R e(t) ❞ F P (s) = u(s)<br />

e(s) = V R (38)<br />

PI-Regler:<br />

u(t) = V R<br />

(<br />

e(t)+ 1 T n<br />

∫<br />

)<br />

e(t)dt<br />

❞ F PI (s) = u(s) (<br />

e(s) = V R 1+ 1 ) ( ) 1+sTn<br />

= V R<br />

sT n sT n<br />

(39)<br />

PD-Regler:<br />

u(t) = V R (e(t)+T v ė(t)) ❞ F PD (s) = u(s)<br />

e(s) = V R(1+sT v ) (40)<br />

PID-Regler:<br />

u(t) = V R<br />

(e(t)+T v ė(t)+<br />

∫ ) e(t)<br />

dt<br />

T n<br />

❞ F PID (s) = u(s) (<br />

e(s) = V R 1+sT v + 1 )<br />

sT n<br />

( ) 1+sTn +s 2 T v T n<br />

= V R<br />

sT n<br />

(41)<br />

• Anfangs- und Endwertsätze <strong>für</strong> F(s) = y(s)<br />

u(s) = Z(s)<br />

N(s)<br />

Wenn die Endwerte lim t→0+ y(t), lim t→∞ y(t) und lim t→0 ẏ(t) existieren und endlich sind,<br />

dann gilt<br />

lim<br />

t→0+<br />

lim<br />

t→∞<br />

y(t) = lim(sF(s)u(s)) <strong>für</strong> deg(Z) < deg(N) (42)<br />

s→∞<br />

y(t) = lim(sF(s)u(s)) (43)<br />

s→0<br />

limẏ(t) = lim(s 2 F(s)u(s)) (44)<br />

t→0 s→∞<br />

wobei deg(Z) und deg(N) die Ordnung des Zähler- bzw. Nennerpolynoms beschreiben.<br />

• Stabilität von linearen Regelkreisen<br />

Ein Regelkreis der Ordnung m,n ∈ N 0 , V S ∈ R, c 0 ,...,c m−1 ∈ R und a 0 ,...,a n−1 ∈ R<br />

— Seite 5/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

mit der Übertragungsfunktion<br />

F S (s) = y(s)<br />

u(s) := V Z(s)<br />

S<br />

N(s) = V S<br />

c 0 +c 1 s+···+c m−1 s m−1 +s m<br />

a 0 +a 1 s+···+a n−1 s n−1 +s n mit m ≤ n (45)<br />

zwischen Eingang u(s) und Ausgang y(s) ist (exponentiell) stabil (d.h. Systemantwort<br />

klingt ab), wenn alle Pole λ i von F S (s), d.h. alle Nullstellen des Nennerpolynoms N(λ i ) =<br />

0 negativen Realteil R{λ i } < 0 <strong>für</strong> alle i = 1,...,n besitzen.<br />

• Routh-Hurwitz Stabilitätskriterium <strong>für</strong> lineare Regelkreise<br />

Dazu untersucht man das charakteristische Polynom n-ter Ordnung des LTI Systems<br />

N(s) = a 0 +a 1 s+a 2 s 2 +···+a n s n , a 0 ,...,a n ∈ R. (46)<br />

Das charakteristische Polynom entspricht dem Nennerpolynom der Übertragungsfunktion<br />

F S (s) in (44). Es gilt: N(s) ist ein Hurwitz-Polynom (d.h. System ist exponentiell stabil),<br />

(i) dann sind alle Koeffizienten a i > 0 in N(s) (notwendige Bedingung, d.h. nicht unbedingt<br />

ausreichend!);<br />

(ii) genau dann, wenn der Koeffizient a n > 0 und alle nordwestlichen Hurwitz-Unterdeterminanten<br />

D i > 0 <strong>für</strong> i = 1,...,n−1 (notwendige & hinreichende Bedingung) 1<br />

Die Unterdeterminanten D i entstehen aus den Determinanten der entsprechenden (i,i)-<br />

Untermatrizen in der linken oberen (“nordwestlichen”) Ecke der Koeffizienten-Matrix<br />

⎡<br />

M n =<br />

⎢<br />

⎣<br />

⎤<br />

a n−1 a n−3 a n−5 ... a n−2n+3 a n−2n+1<br />

a n a n−2 a n−4 ... a n−2n+4 a n−2n+2<br />

0 a n−1 a n−3 ... a n−2n+5 a n−2n+3<br />

0 a n a n−2 ... a n−2n+6 a n−2n+4<br />

.<br />

.<br />

.. .<br />

⎥<br />

0 0 ... ... a 1 0 ⎦<br />

0 0 ... ... a 2 a 0<br />

∈ R n×n . (47)<br />

Zu untersuchen sind D 1 = a n−1 und D 2 =<br />

∣ a ∣<br />

n−1 a n−3 ∣∣∣<br />

, etc. Es gilt a<br />

a n a −k = 0 <strong>für</strong> k > 0.<br />

n−2<br />

• Kausalität (Realisierbarkeit)<br />

Ein lineares dynamisches System F S (s) = y(s) := V u(s) S Z(s) = V N(s) S c 0+c 1 s+···+c m−1 s m−1 +s m<br />

a 0 +a 1 s+···+a n−1<br />

wird<br />

s n−1 +s n<br />

kausal genannt, wenn m ≤ n. Die Systemantwort y(t) eines kausalen Systems hängt<br />

lediglich vom (vorangegangenen) Verlauf der Eingangsgröße u(τ) mit 0 ≤ τ ≤ t ab.<br />

• Zustandsdarstellung eines LTI 1 Systems (Regelungsnormalform)<br />

1 Falls alle Koeffizienten a i > 0 kann über das Liénard-Chipart-Kriterium die Anzahl der zu untersuchenden<br />

Determinaten reduziert werden. Es müssen lediglich die Determinanten D i mit ungeradem Index i = 1,3,5,...<br />

oder geradem Index i = 2,4,6,... auf Positivität geprüft werden. Dieser Sachverhalt basiert auf der linearen<br />

Abhängigkeit der Determinaten <strong>für</strong> a i > 0 <strong>für</strong> alle i = 1,2,3,....<br />

1 linearen, zeit-invarianten (engl. linear, time-invariant)<br />

— Seite 6/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

Zustandsgleichung (Vektordifferentialgleichung) und Ausgangsgleichung eines dynamischen<br />

Systems n-ter Ordnung<br />

}<br />

ẋ(t) = Ax(t)+bu(t) ,x(0) = x 0 ∈ R n<br />

y(t) = c ⊤ (48)<br />

x(t)<br />

mit<br />

– dem Zustandsvektor x(t) = ( x 1 (t), ..., x n (t) ) ⊤<br />

∈ R (z.B. Motorstrom & -drehzahl)<br />

– der Stellgröße u(t) ∈ R (z.B. Motorspannung oder -moment)<br />

– der Systemmatrix in Regelungsnormalform (RNF)<br />

⎡<br />

⎤<br />

0 1 0 ... ... 0<br />

0 0 1 0 ... 0<br />

.<br />

A =<br />

. .. . .. . .. .<br />

∈ R n×n (49)<br />

⎢ 0 ... ... 0 1 0<br />

⎥<br />

⎣ 0 ... ... ... 0 1 ⎦<br />

−a 0 −a 1 ... ... ... −a n−1<br />

– dem Steuer-/Einkoppelvektor b = ( 0, ..., 0, V S<br />

) ⊤<br />

∈ R<br />

n<br />

– dem Auskoppel-/Ausgangsvektor c ⊤ = ( c 0 , c 1 , ..., c n−1<br />

)<br />

∈ R<br />

n<br />

• Übertragungsfunktion (aus Regelungsnormalform)<br />

Darstellung im Laplace-Bereich mit s = σ +iω ∈ C ergibt<br />

F S (s) = y(s)<br />

u(s) = c⊤ (sI n −A) −1 c 0 +c 1 s+···+c n−1 s n−1<br />

b = V S<br />

(50)<br />

a 0 +a 1 s+···+a n−1 s n−1 +s n<br />

• Stabilität eines LTI Systems in Zustandsdarstellung (47)<br />

Ein System der Form (47) ist<br />

– <strong>für</strong> jeden Anfangswert x 0 ∈ R n (exponentiell) stabil und<br />

– <strong>für</strong> jeden Anfangswertx 0 ∈ R n und jeden beschränkten Eingangu(·) bounded-input,<br />

bounded-output (BIBO) stabil,<br />

wenn alle Eigenwerte λ 1 ,...,λ n ∈ C der Matrix A negativen Realteil besitzen, d.h.<br />

R{λ i } < 0 <strong>für</strong> alle i ∈ {1,...,n}. Die Eigenwerte können durch Nullsetzen der charakteristischen<br />

[ Gleichung χ A (s) := det(sI n − A) = 0 bestimmt werden. Hierbei entspricht<br />

1<br />

.<br />

I n = .. ∈ R<br />

1]<br />

n×n der Einheitsmatrix der Ordnung n. Die Eigenwerte λ i sind identisch<br />

mit den Polen der Übertragungsfunktion (49).<br />

— Seite 7/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

5 Grundlagen zu Drehfeldmaschinen (siehe [5] und [6])<br />

Wichtige Grundlage zum Verständnis von Drehfeldmaschinen ist die Zeigertheorie (engl. space<br />

vector theory). Hierzu siehe Abbildung 1.<br />

u<br />

v<br />

w<br />

q<br />

q ′<br />

i q′<br />

s<br />

i β s<br />

β<br />

ω s = ˙φ s<br />

i s = i s s<br />

d<br />

ω k = ˙φ k<br />

i d s<br />

b<br />

i q s<br />

b<br />

φ s<br />

φ k<br />

i d′<br />

s<br />

d ′<br />

ω r = ˙φ r<br />

c<br />

b r<br />

c r<br />

a r<br />

a<br />

b r<br />

c<br />

c r<br />

i α s<br />

a r<br />

a<br />

φ r<br />

α<br />

Abbildung 1: Zeigertheorie: Maschine mit Anschlussklemmen u, v, w, Stator-Wicklungen a, b, c und<br />

Rotor-Wicklungen a r , b r , c r (links) und unterschiedliche Koordinatensysteme (rechts):<br />

• 3-phasiges Koordinatensystem (a,b,c),<br />

• statorfestes s-Koordinatensystem (α,β),<br />

• rotorfestes r-Koordinatensystem (d ′ ,q ′ ) und<br />

• beliebiges k-Koordinatensystem (d,q)<br />

und Statorstrom i s mit Länge ‖i s ‖ = ‖i s s‖ =<br />

√<br />

(i α s) 2 +(i β s) 2 = ‖is‖ r = ∥ ∥i s<br />

k ∥ mit entsprechenden<br />

Komponenten (z.B. i α s und i β s im Stator-Koordinatensystem).<br />

• Clarke-Transformation von Stranggrößen x abc in Statorgrößen x s = T C x abc mit<br />

T C : R 3 → R 3 ,<br />

⎡<br />

⎛ ⎞ ⎛ ⎞ 1 − 1<br />

x a x α<br />

⎝x b ⎠ ↦→ ⎝x β ⎠ := 2 2<br />

− 1 ⎤<br />

2 ⎛ ⎞<br />

⎢0 √ 3<br />

x c x 0 3 ⎣ 2<br />

− √ x a<br />

3<br />

⎥⎝x b ⎠ (51)<br />

2 ⎦<br />

x c<br />

} {{ } } {{ } 1 1 1<br />

=:x abc =:x s }<br />

2 2<br />

{{<br />

2<br />

}<br />

=:T C ∈R 3×3<br />

[Analogie zu D. Schröder: a 0 := e j0◦ = 1, a 1 := e j120◦ , a 2 := e j240◦ ]<br />

— Seite 8/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

wobei x ∈ {ψ r , ψ s , u r , i s , ...}. T C ist regulär mit inverser Matrix<br />

⎡ ⎤<br />

1 0 1<br />

C = ⎢<br />

⎣− 1 √<br />

3<br />

2 2<br />

1⎥<br />

⎦ , d.h. xabc = T −1<br />

− 1 2<br />

− √ 3<br />

2<br />

1<br />

T −1<br />

C xs . (52)<br />

Oft wird die Nullkomponente x 0 vernachlässigt (z.B. gilt i 0 = 0 = i a +i b +i c bei Sternschaltung),<br />

dann vereinfachen sich die Clarke-Transformationsmatrizen zu<br />

⎡ ⎤<br />

[ ]<br />

1 0<br />

T C = 2 1 −<br />

1<br />

− 1 2 2<br />

√<br />

3 0 3<br />

− √ ∈ R 2×3 und T −1<br />

3<br />

C = ⎢<br />

⎣− 1 √<br />

3<br />

⎥<br />

2 2 ⎦ ∈ R3×2 .<br />

2 2<br />

− 1 − √ 3<br />

2 2<br />

• Park-Transformation von Statorgrößen x s in beliebig umlaufendes Koordinatensystem<br />

(z.B. rotorfestes (d ′ ,q ′ )-KoSy oder rotorflussfestes (d,q)-KoSy)<br />

Zur Vereinfachung wird die Nullkomponente (zero-sequence) x 0 in Statorgrößen vernachlässigt<br />

[ ]<br />

cos(φ) −sin(φ)<br />

T P : R → R 2×2 , φ ↦→ T P (φ) :=<br />

[Analogie zu D. Schröder: e jφ ]<br />

sin(φ) cos(φ)<br />

(53)<br />

wobei T P (φ) regulär <strong>für</strong> alle φ ∈ R mit inverser Matrix<br />

[ ]<br />

cos(φ) sin(φ)<br />

T P (φ) −1 = = T P (φ) ⊤ = T P (−φ). [Analogie zu D. Schröder: e −jφ ]<br />

−sin(φ) cos(φ)<br />

(54)<br />

Es gilt entsprechend<br />

[ ]<br />

cos(φ) −sin(φ)<br />

T P (−φ) −1 = = T P (−φ) ⊤ = T P (φ). (55)<br />

sin(φ) cos(φ)<br />

Mithilfe der trigonometrischen Sätze folgt<br />

∀φ 1 ,φ 2 ∈ R: T P (φ 1 )·T P (φ 2 ) = T P (φ 1 +φ 2 ). [Analogie zu D. Schröder: e j(φ 1+φ 2 ) ]<br />

(56)<br />

Für φ = π/2 gilt<br />

[ ]<br />

( π<br />

) 0 −1<br />

J := T P = ,<br />

2 1 0<br />

[Analogie zu D. Schröder: e j π 2 ] (57)<br />

was einer (positiven) Drehung um π/2 eines Vektors x ∈ R 2 entspricht.<br />

— Seite 9/12 —


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

Die Matrizen J und T P (φ) kommutieren, d.h.<br />

[ ]<br />

−sin(φ) −cos(φ)<br />

∀φ ∈ R: JT P (φ) = = T P (φ)J, (58)<br />

cos(φ) −sin(φ)<br />

Für d dt<br />

und<br />

φ =: ω gilt<br />

[Analogie zu D. Schröder: e j(π 2 +φ) = e j(φ+π 2 ) = je jφ ].<br />

[ ]<br />

Ṫ P (φ) := d −sin(φ) −cos(φ)<br />

dt T P(φ) = ω = ωJT P (φ) = ωT P (φ)J (59)<br />

cos(φ) −sin(φ)<br />

[Analogie zu D. Schröder:<br />

Ṫ P (φ) −1 := d dt T P(φ) −1 = ω<br />

[Analogie zu D. Schröder:<br />

[<br />

−sin(φ) cos(φ)<br />

d<br />

dt ejφ = jωe jφ = e j(π 2 +φ) = ωe j(φ+π 2 ) ]<br />

−cos(φ) −sin(φ)<br />

]<br />

d<br />

dt e−jφ = −jωe −jφ ]<br />

= −ωJT P (φ) −1 = −ωT P (φ) −1 J<br />

Mit der Konvention φ k = φ lässt sich mithilfe der Park-Transformation T P (φ k ) ausgehend<br />

von einem statorfesten Koordinatensystem (α,β) in ein (beliebiges) umlaufendes<br />

k-Koordinatensystem (d,q) (Superskript x k = (x d , x q ) ⊤ ) übergegangen werden<br />

(60)<br />

x k = T P (−φ k )x s = T P (φ k ) −1 x s =⇒ x s = T P (φ k )x k , (61)<br />

z.B. zur Rotorfluss-Orientierung (d.h. ψr q = ˙ψ r q = 0 und ψr d = ‖ψr‖ r = ∥ ∥<br />

∥ψ k r<br />

∥, der<br />

Rotorfluss im k-Koordinatensystem liegt ‘exakt’ auf der d-Achse).<br />

— Seite 10/12 —


Optimierungstabelle<br />

y ref<br />

′<br />

y<br />

ref,0<br />

′<br />

✻<br />

1<br />

Führungsgröße<br />

z<br />

z 0<br />

✻<br />

1<br />

Störgröße<br />

— Seite 11/12 —<br />

z F S<br />

y ′ { }} {<br />

y ref ref<br />

✲ ✲ ❡ ✲ ✲ ❄<br />

F R<br />

❡−<br />

y<br />

F G<br />

✲F Sσ<br />

✲F S1<br />

✲F S2<br />

✲<br />

✻−<br />

Führungsglättung<br />

Regler Strecke<br />

Strecke<br />

Regler<br />

y<br />

y ′ ref,0<br />

y ∞<br />

y ′ ref,0<br />

0<br />

✻<br />

1<br />

0<br />

t an<br />

t aus<br />

✲<br />

t<br />

y max<br />

y ′ ref,0 ±2%<br />

✲<br />

t<br />

0<br />

y<br />

y max<br />

V S z 0 ✻<br />

y<br />

V S z 0 ∞<br />

V S z 0<br />

0<br />

❩<br />

❩ Wendetangente<br />

❩<br />

❩<br />

Verhalten bei Sprung der<br />

Einstellung<br />

Führungsgröße y ref (bzw. y<br />

ref ′ ) Störgröße z Nr.<br />

Günstiger<br />

Opt.<br />

Nr. Typ F S Typ F R t an t aus y max y ∞ T ers t an 1 y max 1 y ∞<br />

Bereich Krit. T n V R T v T G (±2%)<br />

T σ T σ y<br />

ref,0<br />

′ y<br />

ref,0<br />

′ T σ T σ V S z 0 V S z 0<br />

V S<br />

1<br />

1<br />

1 PT 1 beliebig I V R BO — — — 4,7 8,4 1,04 1 2 6,3 0,64 0 1<br />

1+sT σ s 2T σ V S<br />

( )<br />

T 1<br />

2<br />

≫ 1 P V R BO — T 1<br />

— — (4,7) (8,4) 1,04 y ∞ V R V S<br />

1 1<br />

T σ 2T σ V S y ′ 2 (4,7) ≈<br />

2<br />

ref,0<br />

1+V R V S 1+V R V S 1+V R V S<br />

√<br />

V S T 1 1+sT n T 1<br />

T1 0,5...1,2<br />

3 PT 2 > 1 PI V R BO T 1 — — 4,7 8,4 1,04 1 2 5,5<br />

0 3<br />

(1+sT 1 )(1+sT σ ) T σ sT n 2T σ V S T σ T 1 /T σ<br />

T 1 1+sT n T 1<br />

— 3,1 ... 4,7 8,4 ... 16,5 1,04 ... 1,43 1 — 1,2...1,6<br />

4<br />

≥ 4 PI V R SO 4T σ — ≈ 10<br />

0 4<br />

T σ sT n 2T σ V S 0..4Tσ 4,7 ... 7,6 8,4 ... 13,3 1,04 ... 1,08 1 2 ... 4 T 1 /T σ<br />

( )<br />

T 1<br />

5<br />

≫ 1 PD V R (1+sT v ) BO — T 1<br />

T 2 — (4,7) (8,4) 1,04 y ∞ V R V S<br />

T σ 2T σ V S y ′ 2 4+ T 2 1 1<br />

≈<br />

5<br />

ref,0<br />

1+V R V S T σ 1+V R V S 1+V R V S<br />

√<br />

V S T 1 (1+sT n )(1+sT v ) T 1<br />

T 1 T 2 0,5...0,75<br />

6 PT 3 > 1 PID V R BO T 1 T 2 — 4,7 8,4 1,04 1 2 4,4 √ √<br />

(1+sT 1 )(1+sT 2 )(1+sT σ ) T σ sT n 2T σ V S T 2 0 6<br />

σ T1 /T σ T2 /T σ<br />

√<br />

T 1 (1+sT n )(1+sT v ) T 1<br />

— 3,1 ... 4,7 8,4 ... 16,5 1,04 ... 1,43 1 — T2 1,4...1,8<br />

7 T 2 > T σ ≥ 4 PID V R SO 4T σ T 2 ≈ 10 4 √ 0 7<br />

T σ sT n 2T σ V S 0..4Tσ 4,7 ... 7,6 8,4 ... 13,3 1,04 ... 1,08 1 2 ... 4 T σ T 1 /T σ T2 /T σ<br />

T 1<br />

8<br />

≫ 1 P V R BO — T 1<br />

— — 4,7 8,4 1,04 1 2 (4,7) ≈ 1 1<br />

8<br />

V<br />

V S T σ 2T σ V S V R V S V R V S S<br />

IT 1<br />

sT 1 (1+sT σ ) 1+sT n T 1<br />

— 3,1 16,5 1,43 1 —<br />

1,6<br />

9 beliebig PI V R SO 4T σ — 10<br />

0 9<br />

sT n 2T σ V S 4Tσ 7,6 13,3 1,08 1 4<br />

T 1 /T σ<br />

T 1<br />

10<br />

≫ 1 PD V R (1+sT v ) BO — T 1<br />

T 2 — 4,7 8,4 1,04 1 2 4+ T 2<br />

≈ 1 1<br />

10<br />

V<br />

V S T σ 2T σ V S T σ V R V S V R V S S<br />

IT 2<br />

√<br />

sT 1 (1+sT 2 )(1+sT σ ) (1+sT n )(1+sT v ) T 1<br />

— 3,1 16,5 1,43 1 — T2 1,8<br />

11<br />

beliebig PID V R SO 4T σ T 2 ≈ 10 4 √ 0 11<br />

T 2 > T σ sT n 2T σ V S 4T σ 7,6 13,3 1,08 1 4 T σ T 1 /T σ T2 /T σ<br />

Abbildung 2: Optimierungstabelle <strong>für</strong> exakt bekannte Strecken mit Ordnung n ≤ 3 und rein reellen Polen<br />

(T 1 , T 2 ...große Zeitkonstanten, T σ ...kleine (Ersatz-)Zeitkonstante).<br />

t an<br />

✲<br />

t<br />

✲<br />

t<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

Regelung von regenerativen Energiesystemen <strong>Formelsammlung</strong> (WS13/14)


Regelung von regenerativen Energiesystemen<br />

<strong>Formelsammlung</strong> (WS13/14)<br />

Technische Universität München, Christoph Hackl <strong>Stand</strong> 26. 11. 2013<br />

Literatur<br />

[1] L. Råde, B. Westergren, and P. Vachenauer, Springers mathematische Formeln: Taschenbuch<br />

<strong>für</strong> Ingenieure, Naturwissenschaftler, Informatiker, Wirtschaftswissenschaftler.<br />

Springer-Verlag, 3. ed., 2000.<br />

[2] K. Mertens, Photovoltaik. Hanser Verlag, 2013.<br />

[3] V. Quaschning, Regenerative Energiesysteme. Hanser Verlag, 2011.<br />

[4] R. Kories and H. Schmidt-Walter, Taschenbuch der Elektrotechnik. Verlag Harri Deutsch,<br />

1998.<br />

[5] D. Schröder, Elektrische Antriebe — Grundlagen (3.,erw. Auflage). Springer-Verlag, 2007.<br />

[6] D. Schröder, Elektrische Antriebe - Regelung von Antriebssystemen (3., bearb. Auflage).<br />

Berlin: Springer-Verlag, 2009.<br />

[7] D. Hinrichsen and A. Pritchard, Mathematical Systems Theory I — Modelling, State Space<br />

Analysis, Stability and Robustness. No. 48 in Texts in Applied Mathematics, Berlin:<br />

Springer-Verlag, 2005.<br />

— Seite 12/12 —

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!