Elektrisk systembeskrivelse av kontaktledningsanlegg ver01 Vedlegg script: Forskjell mellom sideversjoner
Hopp til navigering
Hopp til søk
Linje 143: | Linje 143: | ||
yk = rk + %i * xk | yk = rk + %i * xk | ||
ym = gm + %i * bm | ym = gm + %i * bm | ||
Yat = zeros(3,3) | |||
Yat(1,1) = yk+ym | Yat(1,1) = yk+ym | ||
Yat(2,2) = yk+ym | Yat(2,2) = yk+ym |
Sideversjonen fra 21. jul. 2017 kl. 10:50
__NUMBEREDHEADINGS__
Generelt
I dette vedlegget presenteres funksjoner srevet i Scilab som gjør de beregningene som er beskrevet i Lenke: Elektrisk systembeskrivelse. Koden er testet med Scilab 6.0.0.
Lineær modell
Ingen script tilgjengelig
Transmisjonslinje
Enkel linjesløyfe
Teoretisk beskrivelse: Linjesløyfe
Funksjon: admLoop Beskrivelse: Beregner admittansmatrisa for ei linjesløyfe Kalles ved: Y = admLoop(rLoop,xLoop,gLoop,bLoop,l)
function Y = admLoop(r,x,g,b,l) z = r + %i * x y = g + %i * b Z0 = sqrt(z/y) gam = sqrt(z*y) Y = (1/(Z0*sinh(gam*l))) * [cosh(gam*l) -1 ; -1 sinh(gam*l)] endfunction
Variabel | Type | Enhet | Type | Beskrivelse |
---|---|---|---|---|
Y | (2x2) kompleks matrise | S | Resultat | Admittansmatrise for linjesløyfe |
r | reell skalar | Ω/km | Input | Spesifikk serieresistans |
x | reell skalar | Ω/km | Input | Spesifikk seriereaktans |
g | reell skalar | S/km | Input | Spesifikk parallell konduktans |
b | reell skalar | S/km | Input | Spesifikk parallell susceptans |
l | reell skalar | km | Input | Linjesløyfas lengde |
Transmisjonslinje med flere parallelle ledere
Teoretisk beskrivelse: Transmisjonslinje
Funksjon: admLine Beskrivelse: Beregner admittansmatrisa for en transmisjonslinje med n parallelle ledere Kalles ved: Y = admLine(r,x,g,b,l)
function Y=admLine(r,x,g,b,l) n = size(r,1) nul = zeros(n,n) z = r + %i*x y = g + %i*b A = [nul , -z ; -y , nul] [M,fi] = spec(A) // Finds eigenvector matrix (M) and eigenvalue matrix (fi) for i=1:2*n fi(i,i) = exp(fi(i,i)*l) end fi = M*fi*inv(M) fi11 = fi(1:n,1:n) fi12 = fi(1:n,n+1:2*n) fi21 = fi(n+1:2*n,1:n) fi22 = fi(n+1:2*n,n+1:2*n) fi12inv = inv(fi12) Y = [-fi12inv*fi11 , fi12inv ; (fi22*fi12inv*fi11-fi21) , -fi22*fi12inv] endfunction
Variabel | Type | Enhet | Type | Beskrivelse |
---|---|---|---|---|
Y | (2nx2n) kompleks matrise | S | Resultat | Admittansmatrise for flerlederlinje |
r | (nxn) reell matrise | Ω/km | Input | Spesifikk serieresistans |
x | (nxn) reell matrise | Ω/km | Input | Spesifikk seriereaktans |
g | (nxn) reell matrise | S/km | Input | Spesifikk parallell konduktans |
b | (nxn) reell matrise | S/km | Input | Spesifikk parallell susceptans |
l | (nxn) reell matrise | km | Input | Linjas lengde |
Sugetransformator
Teoretisk beskrivelse: Sugetransformator
Funksjon: admBoosterTransformer Beskrivelse: Beregner admittansmatrisa for en sugetransformator Kalles ved: Y = admBoosterTransformer(yk,ym)
function Y=admBoosterTransformer(rk,xk,gm,bm) yk = (rk + %i * xk)^(-1) ym = gm + %i * bm ytmp = [yk+ym,-yk;-yk,yk] Y([1:2],[1:2]) = ytmp Y([1:2],[3:4]) = -ytmp Y([3:4],[1:2]) = -ytmp Y([3:4],[3:4]) = ytmp endfunction
Variabel | Type | Enhet | Type | Beskrivelse |
---|---|---|---|---|
Y | (4x4) kompleks matrise | S | Resultat | Admittansmatrise for sugetransformator |
rk | reell skalar | Ω | Input | Kortslutningsresistans |
xk | reell skalar | Ω | Input | Kortslutningsreaktans |
gm | reell skalar | S | Input | Magnetiseringskonduktans |
bm | reell skalar | S | Input | Magnetiseringssusceptans |
Seksjonering
Teoretisk beskrivelse: Seksjonering
Funksjon: admSeriesImpedance Beskrivelse: Beregner admittansmatrisa for en seksjonering Kalles ved: Y = admSeriesImpedance(g,b)
function Y=admSeriesImpedance(g,b) ytmp = g + %i*b Y(1,1) = ytmp Y(1,2) = -ytmp Y(2,1) = -ytmp Y(2,2) = ytmp endfunction
Variabel | Type | Enhet | Type | Beskrivelse |
---|---|---|---|---|
Y | (2x2) kompleks matrise | S | Resultat | Admittansmatrise for serieimpedans eller seksjonering |
g | reell skalar | S | Input | Konduktans |
b | reell skalar | S | Input | Susceptans |
Kommentar: For en ren seksjonering er g og b lik 0 S, og resultatet blir en nullmatrise. |
Autotransformator
Teoretisk beskrivelse: Autotransformator
Funksjon: admAutoTransformer Beskrivelse: Beregner admittansmatrisa for en autotransformator Kalles ved: Y = admAutoTransformer(rk,xk,gm,bm)
function Y=admAutoTransformer(rk,xk,gm,bm) yk = rk + %i * xk ym = gm + %i * bm Yat = zeros(3,3) Yat(1,1) = yk+ym Yat(2,2) = yk+ym Yat(1,2) = yk-ym Yat(2,1) = yk-ym Yat(3,[1:2]) = -2.0*yk Yat([1:2],3) = -2.0*yk Yat(3,3) = 4.0*yk endfunction