Forskjell mellom versjoner av «Elektrisk systembeskrivelse av kontaktledningsanlegg ver01»

Fra Lærebøker i jernbaneteknikk
Hopp til: navigasjon, søk
m (Transmisjonslinje med flere parallelle ledere)
m (Transmisjonslinje med flere parallelle ledere)
Linje 190: Linje 190:
 
<math> \frac{\mathrm{d}}{\mathrm{d}x} \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} = \begin{pmatrix} \mathbf{0} & \mathbf{-Z} \\ \mathbf{-Y} & \mathbf{0} \end{pmatrix} \cdot \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} = \mathbf{A} \cdot \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} </math>
 
<math> \frac{\mathrm{d}}{\mathrm{d}x} \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} = \begin{pmatrix} \mathbf{0} & \mathbf{-Z} \\ \mathbf{-Y} & \mathbf{0} \end{pmatrix} \cdot \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} = \mathbf{A} \cdot \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} </math>
  
Her er systemmatrisen <math> \mathbf{A} </math> en matrise med dimensjoner (2''n'' x 2''n''). Parametermatrisene <math>\mathbf{R}</math>, <math>\mathbf{X}</math>, <math>\mathbf{G}</math> og <math>\mathbf{B}</math> er bestemt ut ifra benyttede ledermaterialer, isolasjonsmaterialer, jordsmonn, systemets geometri og driftsfrekvens. Bestemmelse av parametermatrisene er beskrevet i [[Parametre_for_flerledersystem|''Lenke:'' Parametre for flerledersystem]]
+
Her er systemmatrisen <math> \mathbf{A} </math> en kompleks matrise med dimensjoner (2''n'' x 2''n''). Parametermatrisene <math>\mathbf{R}</math>, <math>\mathbf{X}</math>, <math>\mathbf{G}</math> og <math>\mathbf{B}</math> er bestemt ut ifra benyttede ledermaterialer, isolasjonsmaterialer, jordsmonn, systemets geometri og driftsfrekvens. Bestemmelse av parametermatrisene er beskrevet i [[Parametre_for_flerledersystem|''Lenke:'' Parametre for flerledersystem]]
  
 
Systemmatrisen <math> \mathbf{A} </math> har 2''n'' egenverdier <math> \lambda_i </math> med 2''n'' tilhørende egenvektorer <math> \mathbf{v_i} </math>. Disse må finnes numerisk av datamaskiner for aktuelle parametermatriser. Systemets løsning blir på formen:
 
Systemmatrisen <math> \mathbf{A} </math> har 2''n'' egenverdier <math> \lambda_i </math> med 2''n'' tilhørende egenvektorer <math> \mathbf{v_i} </math>. Disse må finnes numerisk av datamaskiner for aktuelle parametermatriser. Systemets løsning blir på formen:

Revisjonen fra 23. jan. 2018 kl. 12:25

1 Generelt

Kontaktledningsanlegget overfører effekt mellom matestasjonen og traksjonsmateriell og andre belastninger tilknyttet kontaktledningen. I eldre anlegg skjer overføringen i kontaktledning ved 15 kV nominell spenning, med retur i kjøreskinner ved 0 kV. I nyere anlegg er det innført returledere eller AT-system med positivleder og negativleder. Følgende ledere er vanlige i et kontaktledningsanlegg. Det kan finnes ytterligere ledere parallelt med jernbanetraseen som påvirker eller påvirkes av kontaktledningsanlegget .

Tabell 1: Ledere oversikt
Leder Forklaring Nominell spenning
KL Kontaktledningsanlegg, omfatter kontakttråd og bæreline 15 kV
RR Kjøreskinner 0 kV
RL Returleder 0 kV
LJ Langsgående jordleder 0 kV
FSL Forsterkningsleder 15 kV
PL Positivleder (for AT-system) + 15 kV
NL Negativleder (for AT-system) - 15 kV

I dette kapittelet beskrives en matematisk modell som beregner hvordan strøm og spenning fordeler seg mellom ulike ledere i et slikt flerledersystem som et kontaktledningsanlegg utgjør. Denne modellen er det teoretiske grunnlaget for beregning av:

  • impedansen mellom matestasjon og belastning i kontaktledningssystemet,
  • potensial i returkretsen,
  • indusert spenning i ledere som går parallelt med jernbanetraseen,
  • belastning på enkeltledere og komponenter i kontaktledningsnettet.

I praktisk bruk er regnemodellen programmert inn i et regneprogram der numeriske parametre for beregningen angis og resultatet presenteres i form av grafer og tabeller.

2 Lineær modell

Et kontaktledningsanlegg har en tilnærmet lineær oppførsel og kan derfor for mange formål modelleres som et lineært system.

Figur 1 viser skjematisk en linje med serielle og parallelle komponenter for et system med fire ledere. Serielle komponenter har koplingspunkter til en eller flere av lederne på hver side av komponenten, mens parallelle komponenter kun har en kopling til hver leder.

Admittansmatrisenettverk 4leder.png
Figur 1: Linje med 5 noder, 4 ledere og 6 lineære kretselementer. Strekningen har 18 terminaler: 2 for node T og 4 for hver av de andre nodene.

En lineær elektrisk modell av en komponent eller et system kan uttrykkes ved hjelp av en admittansmatrise. Admittansmatrisen angir strømmen som påtrykkes alle systemets terminaler, som funksjon av spenningen ved hver enkelt node, der både strømmen og spenningen er viserverdier. Skrevet på matriseform:

[math] \begin{pmatrix} I_1 \\ \vdots \\ I_{18} \end{pmatrix} = \begin{pmatrix} Y_{1,1} & \cdots & Y_{1,18} \\ \vdots & \ddots & \vdots \\ Y_{18,1} & \cdots & Y_{18,18} \end{pmatrix} \cdot \begin{pmatrix} U_1 \\ \vdots \\ U_{18} \end{pmatrix} [/math]

For systemet i Figur 1 har vektorene [math] \mathbf{I} [/math] og [math] \mathbf{U} [/math] 18 elementer hver, mens admittansmatrisen [math] \mathbf{Y} [/math] har (18 x 18) elementer.

Hver komponent i systemet kan også uttrykkes ved sin admittansmatrise. Systemets admittansmatrise er summen av admittansmatrisene for alle komponentene.

Under forutsetning av at systemet har minst en jordforbindelse er admittansmatrisen [math] \mathbf{Y} [/math] inverterbar. Det er da mulig å uttrykke alle systemets spenninger som funksjon av påtrykket strøm i alle systemets terminaler, via den inverse admittansmatrisen, som kalles impedansmatrisen. De aller fleste av systemets terminaler blir ikke påtrykket strøm, og i de resterende koplingspunktene er strømmen ofte kjent. Derfor er det ofte av interesse å beregne spenning på bakgrunn av kjente verdier for strøm.

[math] \begin{pmatrix} U_1 \\ \vdots \\ U_{18} \end{pmatrix} = \begin{pmatrix} Z_{1,1} & \cdots & Z_{1,18} \\ \vdots & \ddots & \vdots \\ Z_{18,1} & \cdots & Z_{18,18} \end{pmatrix} \cdot \begin{pmatrix} I_1 \\ \vdots \\ I_{18} \end{pmatrix} [/math]

Ved å sette inn belastninger fra tog eller kortslutninger ved aktuelle terminaler får man som resultat spenningen ved samtlige andre terminaler i systemet. Dette kan så benyttes til videre å beregne strøm- og spenningsprofil for hver enkelt komponent.

Forhold som ikke er lineære, og som derfor må tilnærmes i en lineær modell, er:

  • magnetiske egenskaper for stålskinner, som gjør stålets indre reaktans ulineært avhengig av strøm og frekvens,
  • strømfortrengning i ledere, som gjør den elektriske motstanden i en leder ulineært avhengig av frekvens, og
  • aktivt styrte belastninger, for eksempel tog.

I de kommende avsnittene beskrives hvordan admittansmatrisen etableres for de ulike komponenttypene i et kontaktledningssystem:

3 Transmisjonslinje

3.1 Enkel linjesløyfe

I Script: Enkel linjesløyfe er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.

Den enkleste formen for problemet er ei enkel linjesløyfe. Det finnes mye teori som beskriver problemet, og problemet har kjent analytisk løsning. Det tas i avsnittet en rask gjennomgang av den grunnleggende teorien, fordi den samme framgangsmåten benyttes for analyse av et flerledersystem.

En kontaktledning med retur i kjøreskinner kan forenklet betraktes som ei linjesløyfe, hvis man ser bort fra lekkasje til jordsmonn. Telegraflikningen tar utgangspunkt i et svært kort segment [math]\mathrm{d}x[/math] av linjesløyfa. Linjesegmentet har en seriell resistans [math]R \cdot \mathrm{d}x[/math] og reaktans [math] X \cdot \mathrm{d}x = \omega \cdot L \cdot \mathrm{d}x[/math], og en parallell konduktans [math]G \cdot \mathrm{d}x [/math] og susceptans [math]B \cdot \mathrm{d}x = \omega \cdot C \cdot \mathrm{d}x[/math]. Her er L, C og [math]\omega[/math] henholdsvis linjas induktans, kapasitans og vinkelfrekvens.

Telegraflikningen.png
Figur 2: Telegraflikningene - Kretsskjema for et linjesegment dx

Parametrene kan skrives om slik at impedansen blir [math]Z = R + j X [/math] og admittansen blir [math]Y = G + j B[/math]. Parametrene kan finnes ved beregning eller ved måling.

Det serielle spenningsfallet dU over dette linjesegmentet er gitt av Ohms lov:

(i) [math] \frac{\mathrm{d}U}{\mathrm{d}x} = - Z \cdot I[/math]

Strømmen som lekker gjennom admittansen [math]Y[/math] utgjør forskjellen i strøm over linjesegmentet. Denne strømmen er proporsjonal med spenningen:

(ii) [math] \frac{\mathrm{d}I}{\mathrm{d}x} = - Y \cdot U [/math]

Likningene (i) og (ii) danner et koplet likningssett som i litteraturen kalles for telegraflikningene. Likningssettet kan ordnes med matriserepresentasjon på følgende måte:

[math] \frac{\mathrm{d}}{\mathrm{d}x} \begin{pmatrix} U \\ I \end{pmatrix} = \begin{pmatrix} 0 & -Z \\ -Y & 0 \end{pmatrix} \cdot \begin{pmatrix} U \\ I \end{pmatrix} = \mathbf{A} \cdot \begin{pmatrix} U \\ I \end{pmatrix} [/math]

Her har systemmatrisen A dimensjonene (2 x 2). Løsningen på et slikt likningssett er beskrivet i flere lærebøker i lineæralgebra, for eksempel i Referanse [1]. En rask innføring er gitt i Wikipedia.

Egenverdiene [math] \lambda [/math] til systemet er:

[math] \lambda_1 = +\sqrt{Z \cdot Y} = \gamma [/math]

[math] \lambda_2 = -\sqrt{Z \cdot Y} = -\gamma [/math]

der

[math] \gamma = \sqrt{Z \cdot Y} [/math] er linjas forplantningskonstant.

De tilhørende egenvektorene kan finnes til å være:

[math] \mathbf{v_1} = \begin{pmatrix} 1 \\ -\sqrt{\tfrac{Y}{Z}} \end{pmatrix} = \begin{pmatrix} 1 \\ -\tfrac{1}{Z_0} \end{pmatrix} [/math]

[math] \mathbf{v_2} = \begin{pmatrix} 1 \\ \sqrt{\tfrac{Y}{Z}} \end{pmatrix} = \begin{pmatrix} 1 \\ \tfrac{1}{Z_0} \end{pmatrix} [/math]

der

[math] Z_0 = \sqrt{\tfrac{Z}{Y}} [/math] er linjas karakteristiske impedans.

Forplantningskonstanten [math] \gamma [/math] er en kompleks størrelse og et mål på hvordan en strøm/spenning dempes og forandrer fase langs ei transmisjonslinje. Den karakteristiske impedansen [math] Z_0 [/math] er et mål på sammenhengen mellom strøm og spenning i ei transmisjonslinje.

Egenvektorene [math] \mathbf{v_1} [/math] og [math] \mathbf{v_1} [/math] ordnes så i en egenvektormatrise [math] \mathbf{V} [/math]:

[math] \mathbf{V} = \begin{pmatrix} \mathbf{v_1} & \mathbf{v_2} \end{pmatrix} = \begin{pmatrix} 1 & 1 \\ - \tfrac{1}{Z_0} & \tfrac{1}{Z_0} \end{pmatrix} [/math]

Egenvektormatrisen [math] \mathbf{V} [/math] har følgende inverse:

[math] \mathbf{V^{-1}} = \frac{1}{2} \begin{pmatrix} 1 & - Z_0 \\ 1 & Z_0 \end{pmatrix} [/math]

Løsningen er gitt av følgende uttrykk:

[math] \begin{pmatrix} U \\ I \end{pmatrix}(x) = \mathbf{V} \cdot e^{\mathbf{\Lambda} x} \cdot \mathbf{V^{-1}} \cdot \begin{pmatrix} U_0 \\ I_0 \end{pmatrix} = \mathbf{\Phi}(x) \cdot \begin{pmatrix} U_0 \\ I_0 \end{pmatrix} [/math]

der:
[math] U_0 [/math] og [math] I_0 [/math] er strøm og spenning ved posisjonen [math] x = 0 [/math],
[math] e^{\mathbf{\Lambda} x} [/math] er en diagonal 2x2-matrise med diagonale elementer [math]e^{\gamma x}[/math] og [math]e^{-\gamma x}[/math] , og
[math] \mathbf{\Phi}(x) [/math] er den resulterende løsningsmatrisen.

Skrevet fullt ut er løsningsmatrisen:

[math] \mathbf{\Phi}(x) = \begin{pmatrix} \cosh{(\gamma x)} & -Z_0 \sinh({\gamma x)} \\ - \tfrac{1}{Z_0} \sinh{(\gamma x)} & \cosh{(\gamma x)} \end{pmatrix} [/math]

Løsningen gir en komplett beskrivelse av strøm og spenning i alle punkter langs ei transmisjonslinje, dersom strømmen og spenningen ved posisjon [math] x = 0 [/math] er kjent.

Admittansbeskrivelse 2porter.png
Figur 3: Et linjesegment med lengde [math] l [/math] har lineære karakteristikker og kan beskrives komplett ved en 2x2 løsningsmatrise [math] \mathbf{\Phi_l} [/math] eller admittansmatrise [math] \mathbf{Y_l} [/math]

For ei linje av lengde [math] l [/math] som vist i Figur 3 mellom punktene [math] R [/math] og [math] S [/math] kan man finne sammenhengen mellom strøm og spenning ved de to linjeendene ved:

[math] \begin{pmatrix} U_S \\ -I_S \end{pmatrix} = \mathbf{\Phi}(l) \cdot \begin{pmatrix} U_R \\ I_R \end{pmatrix} = \mathbf{\Phi_l} \cdot \begin{pmatrix} U_R \\ I_R \end{pmatrix} [/math]

der

[math] U_R = U(0), \; U_S = U(l), \; I_R = I(0), \; I_S = -I(l) [/math]. Grunnen til det negative fortegnet for [math] I_S [/math] er at vi nå definerer strømmen som positiv inn i linja ved begge linjeender.

Man vil i mange tilfeller kjenne spenningen ved de to linjeendene, og ønske å beregne strømmen. Da kan likningssettet skrives om på admittansmatriseform slik at:

[math] \begin{pmatrix} I_R \\ I_S \end{pmatrix} = \mathbf{Y_l} \cdot \begin{pmatrix} U_R \\ U_S \end{pmatrix} [/math]

Man kan finne at matrisen [math] \mathbf{Y_l} [/math] er:

[math] \mathbf{Y_l} = \frac{1}{Z_0 \sinh{(\gamma l)}} \begin{pmatrix} \cosh{(\gamma l)} & -1 \\ -1 & \cosh{(\gamma l)} \end{pmatrix} [/math]

3.2 Transmisjonslinje med flere parallelle ledere

I Script: Transmisjonslinje er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.

Ved flere parallelle ledere blir framgangsmåten svært lik den som er skissert for ei linjesløyfe, men matematikken blir så komplisert at den ikke kan gjøres analytisk. Den må derfor i stedet gjennomføres numerisk, vanligvis ved hjelp av en datamaskin.

I tillegg til selvimpedans for hver enkelt leder vil det også være en gjensidig impedans mellom alle lederne. Den gjensidige impedansen gir et tillegg til spenningen over en leder som følge av strømmen i hver av de andre lederne. For et system med n (1, 2, ... n) ledere kan dette beskrives med matriser på følgende måte:

[math] \frac{\mathrm{d}}{\mathrm{d}x} \begin{pmatrix} U_1 \\ \vdots \\ U_n \end{pmatrix} = -\left[ \begin{pmatrix} R_{11} & \cdots & R_{1n} \\ \vdots & \ddots & \vdots \\ R_{n1} & \cdots & R_{nn} \end{pmatrix} + j \cdot \begin{pmatrix} X_{11} & \cdots & X_{1n} \\ \vdots & \ddots & \vdots \\ X_{n1} & \cdots & X_{nn} \end{pmatrix} \right] \cdot \begin{pmatrix} I_1 \\ \vdots \\ I_n \end{pmatrix} = -\begin{pmatrix} Z_{11} & \cdots & Z_{1n} \\ \vdots & \ddots & \vdots \\ Z_{n1} & \cdots & Z_{nn} \end{pmatrix} \cdot \begin{pmatrix} I_1 \\ \vdots \\ I_n \end{pmatrix} [/math]

På tilsvarende måte kan det være en konduktiv eller kapasitiv lekkasje mellom hver leder i et ledningssystem. Dette kan beskrives på matriseform på følgende måte for de samme lederne:

[math] \frac{\mathrm{d}}{\mathrm{d}x} \begin{pmatrix} I_1 \\ \vdots \\ I_n \end{pmatrix} = -\left[ \begin{pmatrix} G_{11} & \cdots & G_{1n} \\ \vdots & \ddots & \vdots \\ G_{n1} & \cdots & G_{nn} \end{pmatrix} + j \cdot \begin{pmatrix} B_{11} & \cdots & B_{1n} \\ \vdots & \ddots & \vdots \\ B_{n1} & \cdots & B_{nn} \end{pmatrix} \right] \cdot \begin{pmatrix} U_1 \\ \vdots \\ U_n \end{pmatrix} = -\begin{pmatrix} Y_{11} & \cdots & Y_{1n} \\ \vdots & \ddots & \vdots \\ Y_{n1} & \cdots & Y_{nn} \end{pmatrix} \cdot \begin{pmatrix} U_1 \\ \vdots \\ U_n \end{pmatrix} [/math]

På kompakt matriseform kan dette skrives:

[math] \frac{\mathrm{d}}{\mathrm{d}x} \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} = \begin{pmatrix} \mathbf{0} & \mathbf{-Z} \\ \mathbf{-Y} & \mathbf{0} \end{pmatrix} \cdot \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} = \mathbf{A} \cdot \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix} [/math]

Her er systemmatrisen [math] \mathbf{A} [/math] en kompleks matrise med dimensjoner (2n x 2n). Parametermatrisene [math]\mathbf{R}[/math], [math]\mathbf{X}[/math], [math]\mathbf{G}[/math] og [math]\mathbf{B}[/math] er bestemt ut ifra benyttede ledermaterialer, isolasjonsmaterialer, jordsmonn, systemets geometri og driftsfrekvens. Bestemmelse av parametermatrisene er beskrevet i Lenke: Parametre for flerledersystem

Systemmatrisen [math] \mathbf{A} [/math] har 2n egenverdier [math] \lambda_i [/math] med 2n tilhørende egenvektorer [math] \mathbf{v_i} [/math]. Disse må finnes numerisk av datamaskiner for aktuelle parametermatriser. Systemets løsning blir på formen:

[math] \begin{pmatrix} \mathbf{U} \\ \mathbf{I} \end{pmatrix}(x) = \mathbf{V} \cdot e^{\mathbf{\Lambda} x} \cdot \mathbf{V^{-1}} \cdot \begin{pmatrix} \mathbf{U_0} \\ \mathbf{I_0} \end{pmatrix} = \mathbf{\Phi}(x) \cdot \begin{pmatrix} \mathbf{U_0} \\ \mathbf{I_0} \end{pmatrix} [/math]

der:
[math] \mathbf{U_0} [/math] og [math] \mathbf{I_0} [/math] er strøm og spenning ved posisjonen [math] x = 0 [/math],
[math] e^{\mathbf{\Lambda} x} [/math] er en diagonal (2n x 2n)-matrise med diagonale elementer [math]e^{\lambda_i x}[/math],
[math] \mathbf{V} [/math] er (2n x 2n) egenvektormatrisen der kolonne nr. i er egenvektoren [math] \mathbf{v_i} [/math] til egenverdien [math] \lambda_i [/math] , og
[math] \mathbf{\Phi}(x) [/math] er den resulterende løsningsmatrisen med dimensjoner (2n x 2n).

Admittansbeskrivelse n linjer.png
Figur 4: Et linjesegment med lengde [math] l [/math] har lineære karakteristikker og kan beskrives komplett ved en 2nx2n løsningsmatrise [math] \mathbf{\Phi_l} [/math] eller admittansmatrise [math] \mathbf{Y_l} [/math]

På samme måte som for den enkle linjesløyfa kan man for ei linje med lengde [math] l [/math] som vist i Figur 4 mellom nodene [math] R [/math] og [math] S [/math] finne sammenhengen mellom strøm og spenning ved de to linjeendende ved:

[math] \begin{pmatrix} \mathbf{U_S} \\ -\mathbf{I_S} \end{pmatrix} = \mathbf{\Phi}(l) \cdot \begin{pmatrix} \mathbf{U_R} \\ \mathbf{I_R} \end{pmatrix} = \mathbf{\Phi_l} \cdot \begin{pmatrix} \mathbf{U_R} \\ \mathbf{I_R} \end{pmatrix} [/math]

For å finne admittansmatrisebeskrivelsen deles løsningsmatrisen [math] \mathbf{\Phi_l} [/math] inn i fire (n x n) submatriser på følgende måte:

[math] \begin{pmatrix} \mathbf{U_S} \\ -\mathbf{I_S} \end{pmatrix} = \begin{pmatrix} \mathbf{\Phi_{l,11}} & \mathbf{\Phi_{l,12}} \\ \mathbf{\Phi_{l,21}} & \mathbf{\Phi_{l,22}} \end{pmatrix} \cdot \begin{pmatrix} \mathbf{U_R} \\ \mathbf{I_R} \end{pmatrix} [/math]

En kan da finne at beskrivelsen på admittansmatriseform [math] \mathbf{Y_l} [/math] blir:

[math] \begin{pmatrix} \mathbf{I_R} \\ \mathbf{I_S} \end{pmatrix} = \mathbf{Y_l} \cdot \begin{pmatrix} \mathbf{U_R} \\ \mathbf{U_S} \end{pmatrix} = \begin{pmatrix} \mathbf{Y_{l,11}} & \mathbf{Y_{l,12}} \\ \mathbf{Y_{l,21}} & \mathbf{Y_{l,22}} \end{pmatrix} \cdot \begin{pmatrix} \mathbf{U_R} \\ \mathbf{U_S} \end{pmatrix} [/math]

der

[math] \begin{align} \mathbf{Y_{l,11}} & = - \mathbf{\Phi_{l,12}^{-1}} \cdot \mathbf{\Phi_{l,11}} \\ \mathbf{Y_{l,12}} & = \mathbf{\Phi_{l,12}^{-1}} \\ \mathbf{Y_{l,21}} & = \mathbf{\Phi_{l,22}} \cdot \mathbf{\Phi_{l,12}^{-1}} \cdot \mathbf{\Phi_{l,11}} - \mathbf{\Phi_{l,21}} \\ \mathbf{Y_{l,22}} & = - \mathbf{\Phi_{l,22}} \cdot \mathbf{\Phi_{l,12}^{-1}} \end{align}[/math]

4 Sugetransformator

I Script: Sugetransformator er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.

Ekvivalentskjema for en sugetransformator kan tegnes som vist i Figur 5, der ym er magnetiseringsadmittansen og yk er kortslutningsadmittansen.

Sugetransformator ekvivalentskjema.png
Figur 5: Ekvivalentskjema for sugetransformator. TODO: Endre fra "Node 1" og "Node 2" til "Node R" og "Node S"

Basert på dette skjemaet kan en admittansmatrise etableres som følger:

[math] \begin{pmatrix} I_{R-KL} \\ I_{R-RR} \\ I_{S-KL} \\ I_{S-RR} \end{pmatrix} = \begin{pmatrix} (y_m + y_k) & -y_k & -(y_m + y_k) & y_k \\ -y_k & y_k & y_k & -y_k \\ -(y_m + y_k) & y_k & (y_m + y_k) & -y_k \\ y_k & -y_k & -y_k & y_k \end{pmatrix} \cdot \begin{pmatrix} U_{R-KL} \\ U_{R-RR} \\ U_{S-KL} \\ U_{S-RR} \end{pmatrix} [/math]

For å finne verdiene for parametrene ym og yk leser man datablad for en sugetransformator. I Norge benyttes to typer sugetransformator: 600 A og 800 A, og dataene for disse er gjengitt under. For de fleste praktiske formål kan ym settes til 0 S, kun ved analyse av brudd i returkretsen kan denne få betydning.

Tabell 2: Data for sugetransformatorer
Parameter Symbol Enhet Sugetransformator
600 A
Sugetransformator
600 A
Sugetransformator
800 A
Merkeytelse (datablad) Sr kVA 55 55 95
Merkestrøm (datablad) Ir A 600 600 800
Merkespenning (datablad) Ur V 91,7 91,7 118,75
Merkeimpedans (beregnet) Zr Ω 0,1528 0,1528 0,1484
Kortslutningstap (målt) Pk W 1 099 4 429 7 012
Tomgangstap (målt) Po W 9 39 40
Tomgangsstrøm (målt) Io A 0,2 0,6 (antatt) 0,5 (antatt)
Kortslutningsimpedans (målt) ez / ek % 22,89 20,24 18,2
Kortslutningsresistans (beregnet) er % 10,47 8,05 7,38
Kortslutningsreaktans (beregnet) ex % 20,36 18,57 16,64
Magnetiseringsadmittans (beregnet) ym mS 2,18 6,55 4,21
Magnetiseringskonduktans (beregnet) gm mS 1,07 4,64 2,84
Magnetiseringsreaktans (beregnet) bm mS 1,90 4,61 3,11
Produsent ABB Møre Trafo Møre trafo
Kilde EK.800131 EK.800207 EK.800208
Kortslutningsadmittans yk S 13,08 - j 25,42 12,84 - j 29,65 14,62 - j 32,87
Magnetiseringsadmittans ym mS 1,07 + j 1,90 4,64 + j 4,61 2,84 + j 3,11

5 Seksjonering

I Script: Seksjonering er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.

En seksjonering kan tegnes som en seriell admittans som vist i Figur 6.

Seksjonering ekvivalentskjema.png
Figur 6: Ekvivalentskjema for seksjonering.

Admittansmatrisen for en slik seriell admittans er:

[math] \begin{pmatrix} I_{R} \\ I_{S} \end{pmatrix} = \begin{pmatrix} y & -y \\ -y & y \end{pmatrix} \cdot \begin{pmatrix} U_{R} \\ U_{S} \end{pmatrix} [/math]

For en ren seksjonering er motstanden svært stor, og admittansen y kan settes til [math] y = 0 [/math] S.

6 Autotransformator

I Script: Autotransformator er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.

En autotransformator kan tegnes som vist i Figur 7.

Autotransformator ekvivalentskjema.png
Figur 7: Ekvivalentskjema for autotransformator.

Den tilhørende admittansmatrisen kan finnes ut fra skjemaet til å være:

[math] \begin{pmatrix} I_{NL} \\ I_{PL} \\ I_{RR} \end{pmatrix} = \begin{pmatrix} (y_k + y_m) & (y_k - y_m) & (-2 \cdot y_k) \\ (y_k - y_m) & (y_k + y_m) & (-2 \cdot y_k) \\ (-2 \cdot y_k) & (-2 \cdot y_k) & (4 \cdot y_k) \end{pmatrix} \cdot \begin{pmatrix} U_{NL} \\ U_{PL} \\ U_{RR} \end{pmatrix} [/math]

Parametrene yk og ym kan beregnes ut ifra transformatorens datablad. Merk: For krestsskjemaet og admittansmatrisen som er angitt over, må verdiene autotransformatorens høyspentvikling brukes. Datablad for en 5 MVA autotransformator fra ABB angir (kilde: EB.800106).

Tabell 3: Data for autotransformatorer
Parameter Symbol Enhet Autotransformator
5 MVA
høyspentvikling
Merkeytelse (datablad) Sr kVA 5000
Merkespenning (datablad) Ur V 33 000
Merkestrøm (beregnet) Ir A 151,5
Merkeimpedans (beregnet) Zr Ω 217,8
Kortslutningstap (datablad) Pk W 10 600
Tomgangstap (datablad) Po W 2 600
Tomgangsstrøm Io A Ikke angitt
Kortslutningsimpedans (datablad) ez / ek % 0,4
Kortslutningsimpedans (beregnet) ez / ek Ω 0,8712
Kortslutningsresistans (beregnet) er Ω 0,4617
Kortslutningsreaktans (beregnet) ex Ω 0,7388
Magnetiseringskonduktans (beregnet) gm µS 2,3875
Magnetiseringssusceptans
antatt lik 2,5 · gm
bm µS 5,9688
Produsent ABB
Kilde EK.800106
Kortslutningsadmittans yk S 0,6082 - j 0.9735
Magnetiseringsadmittans ym µS 2,3875 + j 5,9688

7 Jordforbindelse

En jordforbindelse er en forbindelse mellom et tilkoplingspunkt i systembeskrivelsen (normalt tilkoplet returkretsen), og fjernt jordsmonn. Jordforbindelsen omfatter overgangsmotstanden i en jordelektrode. Admittansbeskrivelsen av en jordforbindelse er:

[math] i_{R}= y_{RE} \cdot u_R [/math]

der

iR er strømmen i jordforbindelsen tilknyttet tilkoplingspunkt R.

yRE er admittansen mellom tilkoplingspunkt R og jordpotensialet e.

uR er spenningen mot jord ved tilkoplingspunkt R.

Verdien for yRE i et tilkoplingspunkt er den inverse overgangsmotstanden:

[math] y_{RE}= \frac{1}{r_{RE}} [/math]

Verdien for overgangsmotstand må måles eller anslås ved beregning. For modellering er det antatt at matestasjoner har en admittans på 0,1 S, og for AT-kiosker er det antatt en admittans på 0,02 S for hver AT-kiosk.

8 Kilde og belastning

Ved kilder og belastninger kan strøm styres mellom ulike ledere i systemet. Ved en kilde går strøm inn i kontaktledningen og ut av kjøreskinnene. Motsatt går strøm ut av kontaktledningen og inn i kjøreskinnene ved belastninger.

Systembeskrivelse kilde.png
Figur 8: Kretsskjema for en kilde modellert som en Norton-ekvivalent

8.1 Kilde

I Script: Kilde er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.

Kilder modelleres som Norton-ekvivalenter som er koplet mellom kontaktledningen og kjøreskinnene.

iK er den Norton-ekvivalente strømmen til kilden. yK er den Norton-ekvivalente admittansen til kilden.

iK og yK bestemmes på følgende måte:

- yK settes til en vilkårlig stor negativ susceptans. Vanligvis velges en størrelse rundt yK = -j 100 S.

- Tomgangsspenningen u0 velges, vanligvis til 16 500 V i norsk banestrømforsyning.

- Den Norton-ekvivalente strømmen iK settes lik [math] u_0 \cdot y_k [/math]. Ved yK = -j 100 S og u0 = 16 500 V blir iK = -j 16 500 kA.

Admittansmatrisen for en kilde er:

[math] \begin{pmatrix} I_{KL} \\ I_{RR} \end{pmatrix} = \begin{pmatrix} y_K & -y_K \\ -y_K & y_K \end{pmatrix} \cdot \begin{pmatrix} U_{KL} \\ U_{RR} \end{pmatrix} [/math]

Strømvektoren for kilden blir:

[math] \begin{pmatrix} I_{KL} \\ I_{RR} \end{pmatrix} = \begin{pmatrix} i_K \\ -i_K \end{pmatrix} [/math]

8.2 Belastning

I Script: Belastning er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.

Belastninger kan modelleres på flere måter:

  1. Som en ren strømkilde. Denne metoden gir god kontroll på belastningsstrømmens amplitude, men dårlig kontroll på fasevinkelen. Fasevinkelen må eventuelt bestemmes iterativt for å styre hvor stor belastningsstrømmens fasevinkel blir.
  2. Som en admittans yB. Denne metoden gir god kontroll på belastsningsstrømmens fasevinkel, men den gir dårlig kontroll på amplituden. Admittansens størrelse må eventuelt bestemmes iterativt for å styre hvor stor belastningsstrømmens amplitude blir.
  3. Som en Norton-ekvivalent, dvs. en kombinert admittans og strømkilde. Denne metoden kan brukes for å styre spenningen ved belastningen, men gir ganske komplisert beregning.
  4. Som en fast aktiv effekt P og reaktiv effekt Q. Denne metoden krever at det gjennomføres en lastflytanalyse (for eksempel med Newton-Raphson-iterasjon). Det krever omfattende programmering å implementere dette, og regnetiden kan bli stor.

Dersom det er flere belastninger samtidig i samme beregning, må det enten gjennomføres en lastflytanalyse, eller så må man akseptere at belastningens fasevinkel ikke blir riktig. En funksjon for lastflytanalyse er ikke utarbeidet her, og det er valgt å bruke metode nummer 1 (konstant strømbelastning) for angivelse av belastningens størrelse.

Belastningsstrømmen iB settes lik den strømmen som belastningen skal trekke i den aktuelle noden. Ved beregning av spenningen ved belastningen, justeres fasevinkelen iterativt slik at strømmens fasevinkel blir lik fasevinkelen til spenningen ved belastningen. Vanligvis blir fasevinkelen svært nær 0°, slik at et tilnærmet resultat kan oppnås ved å sette fasevinkelen til 0°.

Strømvektoren for belastningen blir:

[math] \begin{pmatrix} I_{KL} \\ I_{RR} \end{pmatrix} = \begin{pmatrix} -i_B \\ i_B \end{pmatrix} [/math]

9 Referanser

[1] Edwards, Penney: Elementary Linear Algebra, Pearson, 1987. ISBN 9780132582605.