Elektrisk systembeskrivelse av kontaktledningsanlegg ver01: Forskjell mellom sideversjoner

Fra Lærebøker i jernbaneteknikk
Hopp til navigering Hopp til søk
(→‎Autotransformator: Rettet opp kildehenvisning slik at den viser til autotransformator ikke åk konsoll i ProArc)
 
(218 mellomliggende versjoner av 2 brukere er ikke vist)
Linje 1: Linje 1:
__NUMBEREDHEADINGS__
__NUMBEREDHEADINGS__
= Generelt =
= 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 vil bli omtalt videre i dette kapittelet:
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 .


<figtable id="tab:Ledere_oversikt">
<figtable id="tab:Ledere_oversikt">
Linje 7: Linje 7:
|+ <caption>Ledere oversikt</caption>
|+ <caption>Ledere oversikt</caption>
|-
|-
! Leder
! Leder !! Forklaring !! Nominell spenning
! Forklaring
! Nominell spenning
|-
|-
| KL
| KL || Kontaktledning, omfatter kontakttråd og bæreline || 15 kV
| Kontaktledningsanlegg, omfatter kontakttråd og bæreline
| 15 kV
|-
|-
| RR
| RR || Kjøreskinner || 0 kV
| Kjøreskinner
| 0 kV
|-
|-
| RL
| RL || Returleder || 0 kV
| Returleder
| 0 kV
|-
|-
| FSL
| LJ || Langsgående jordleder || 0 kV
| Forsterkningsleder
| 15 kV
|-
|-
| PL
| FSL || Forsterkningsleder || 15 kV
| Positivleder (for AT-system)
| normalt + 15 kV
|-
|-
| NL
| PL || Positivleder (for AT-system) || + 15 kV
| Negativleder (for AT-system)
|-
| normalt - 15 kV
| NL || Negativleder (for AT-system) || - 15 kV
|}
|}
</figtable>
</figtable>


I dette kapittelet vil det bli beskrevet en matematisk modell som beskriver hvordan strøm og spenning fordeler seg mellom ulike ledere i et slikt flerledersystem som et kontaktledningsanlegg utgjør. En slik modell kan brukes til å beregne:
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,
* impedansen mellom matestasjon og belastning i kontaktledningssystemet,
* potensial i returkretsen,
* potensial i returkretsen,
Linje 43: Linje 31:
* belastning på enkeltledere og komponenter i kontaktledningsnettet.
* belastning på enkeltledere og komponenter i kontaktledningsnettet.


= Introduksjon: Enkel transmisjonslinje =
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.
== Telegraflikningene ==
 
Den enkleste formen for problemet er ei enkelt linjesløyfe. En kontaktledning med retur i kjøreskinner kan forenklet betraktes på denne måten, hvis man ser bort fra lekkasje til jordsmonn. Telegraflikningen tar utgangspunkt i et svært kort linjesegment d''x'' av linjesløyfa. Linjesegmentet har en seriell resistans ''R'' · d''x'' og reaktans ''X'' · d''x'', og en parallell konduktans ''G'' · d''x'' og susceptans ''B'' · d''x''.
= 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. Ved å anta lineær oppførsel neglisjeres ikke-lineære fenomener, i første rekke magnetisk metning og hysterese i magnetiske materialer som jern.
 
<xr id="fig:LinearModell"/> 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.
 
<figure id="fig:LinearModell">
[[Bilde:Admittansmatrisenettverk 4leder.png|thumb|none|680px| <caption>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.</caption>]]
 
</figure>
 
En lineær elektrisk modell av en komponent eller et system av komponenter kan uttrykkes ved hjelp av en admittansmatrise. Admittansmatrisen angir strømmen som påtrykkes alle systemets terminaler, som funksjon av spenningen ved hver enkelt terminal, der både strømmen og spenningen er viserverdier. Skrevet på matriseform for et system med ''n'' terminaler:
 
<math> \begin{pmatrix} I_1 \\ \vdots \\ I_{n} \end{pmatrix} = \begin{pmatrix} Y_{1,1} & \cdots & Y_{1,n} \\ \vdots & \ddots & \vdots \\ Y_{n,1} & \cdots & Y_{n,n} \end{pmatrix} \cdot \begin{pmatrix} U_1 \\ \vdots \\ U_{n} \end{pmatrix} </math>
 
For systemet i <xr id="fig:LinearModell"/> har vektorene <math> \mathbf{I} </math> og <math> \mathbf{U} </math> n elementer hver, en for hver terminal, mens admittansmatrisen <math> \mathbf{Y} </math> har (n x n) elementer.
 
Hver komponent i systemet kan også uttrykkes ved sin admittansmatrise. Systemets admittansmatrise er summen av admittansmatrisene for alle komponentene. Elementer i admittansmatrisen som mellom terminaler som ikke er sammenkoplet, har admittansverdien 0.
 
Under forutsetning av at en komponent eller et system har minst en jordforbindelse er den tilhørende admittansmatrisen <math> \mathbf{Y} </math> for komponenten eller systemet inverterbar. Det er da mulig å uttrykke alle systemets spenninger som funksjon av påtrykket strøm i alle terminaler, via den inverse admittansmatrisen, som kalles impedansmatrisen. De aller fleste av terminalene blir ikke påtrykket strøm, og i de resterende terminalene 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_{n} \end{pmatrix} = \begin{pmatrix} Z_{1,1} & \cdots & Z_{1,n} \\ \vdots & \ddots & \vdots \\ Z_{n,1} & \cdots & Z_{n,n} \end{pmatrix} \cdot \begin{pmatrix} I_1 \\ \vdots \\ I_{n} \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,
* temperaturavhengighet for parametre som ohmske motstander,
* metning og hysterese i transformatorer, og
* aktivt styrte belastninger, for eksempel tog.
 
I de kommende avsnittene beskrives hvordan admittansmatrisen etableres for de ulike komponenttypene i et kontaktledningssystem:
* Seriell komponent: [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01#Transmisjonslinje|Transmisjonslinje]]
* Seriell komponent: [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01#Sugetransformator|Sugetransformator]]
* Seriell komponent: [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01#Seksjonering|Seksjonering eller brudd]]
* Parallell komponent: [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01#Autotransformator|Autotransformator]]
* Parallell komponent: [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01#Jordforbindelse|Jordforbindelse]]
* Parallell komponent: [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01#Kilde og belastning|Kilde og belastning]]
 
= Transmisjonslinje =
== Generelt ==
Modellering av transmisjonslinje er gjort for:
* [[Elektrisk systembeskrivelse av kontaktledningsanlegg ver01#Enkel linjesløyfe|Enkel linjesløyfe]]: Modellering av en enkelt sløyfe med fram- og tilbakeleder.
* [[Elektrisk systembeskrivelse av kontaktledningsanlegg ver01#Transmisjonslinje med flere parallelle ledere|Transmisjonslinje med flere parallelle ledere]]: Modellering av system med flere parallelle ledere og jord.
 
== Enkel linjesløyfe ==
I [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01_Vedlegg_script#Enkel_linjesl.C3.B8yfe|''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.


<figure id="fig:Telegraflikningene">
<figure id="fig:Telegraflikningene">
[[Bilde:Telegraflikningen.png|thumb|none|480px|]]
[[Bilde:Telegraflikningen.png|thumb|none|460px| <caption>Telegraflikningene - Kretsskjema for et linjesegment d''x''</caption>]]
<caption>Telegraflikningene: Kretsskjema for et linjesegment d''x''</caption>
 
</figure>
</figure>


Parametrene kan skrives om om slik at impedansen blir <math>Z = R + j X </math> og admittansen blir <math>Y = G + j B</math>
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 d''U'' over dette linjesegmentet er gitt av Ohms lov:
Det serielle spenningsfallet d''U'' over dette linjesegmentet er gitt av Ohms lov:
Linje 64: Linje 102:
Likningene ''(i)'' og ''(ii)'' danner et koplet likningssett som i litteraturen kalles for telegraflikningene. Likningssettet kan ordnes med matriserepresentasjon på følgende måte:
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} </math>
<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>


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 [https://en.wikipedia.org/wiki/Matrix_differential_equation Wikipedia].
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 [https://en.wikipedia.org/wiki/Matrix_differential_equation Wikipedia].


Egenverdiene <math> \gamma </math> til systemet er:
Egenverdiene <math> \lambda </math> til systemmatrisen '''A''' er:


<math> \lambda_1 = +\sqrt{Z \cdot Y} = \gamma </math>
<math> \lambda_1 = +\sqrt{Z \cdot Y} = \gamma </math>
Linje 76: Linje 114:
der  
der  


<math> \gamma = \sqrt{Z \cdot Y} </math> er linjens transmisjonskonstant.
<math> \gamma = \sqrt{Z \cdot Y} </math> er linjas forplantningskonstant.


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


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


<math> v_2 = \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
der


<math> Z_0 = \sqrt{\tfrac{Z}{Y}} </math> er linjens karakteristiske impedans.
<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_2} </math> ordnes så i en egenvektormatrise <math> \mathbf{V} </math>:


Transmisjonskonstanten <math> \gamma </math> er en kompleks størrelse og et mål på hvordan en strøm/spenning dempes og forandrer fase langs en transmisjonslinje. Den karakteristiske impedansen <math> Z_0 </math> er et mål på sammenhengen mellom strøm og spenning i en transmisjonslinje.
<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:
Løsningen er gitt av følgende uttrykk:


<math> \begin{pmatrix} U \\ I \end{pmatrix} = A \cdot \mathbf{v_1} \cdot e^{\lambda_1 x} + B \cdot \mathbf{v_2} \cdot e^{\lambda_2 x} </math>
<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: <br>
<math> U_0 </math> og <math> I_0 </math> er strøm og spenning ved posisjonen <math> x = 0 </math>, <br>
<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 <br>
<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.
 
<figure id="fig:Admittansbeskrivelse_2porter">
[[Bilde:Admittansbeskrivelse_2porter.png|thumb|none|480px| <caption>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></caption>]]
 
</figure>
 
For ei linje av lengde <math> l </math> som vist i <xr id="fig:Admittansbeskrivelse_2porter" /> 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>
 
== Transmisjonslinje med flere parallelle ledere ==
 
=== Beskrivelse av transmisjonslikningene for flerledersytem ===
I [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01_Vedlegg_script#Transmisjonslinje_med_flere_parallelle_ledere|''Script:'' Transmisjonslinje]] er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet. <br>
For parametrisering av modellen henvises det til [[Parametre_for_flerledersystem|''Lenke:'' Parametre for flerledersystem]].
 
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>
 
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]].
 
På kompakt matriseform kan disse likningene skrives:
 
(i)  <math> \frac{\mathrm{d}}{\mathrm{d}x} \mathbf{U} = \mathbf{-Z} \cdot \mathbf{I} </math>
 
(ii)  <math> \frac{\mathrm{d}}{\mathrm{d}x} \mathbf{I} = \mathbf{-Y} \cdot \mathbf{U} </math>
 
Det er to vanlige tilnærminger for å løse dette likningssettet. De beskrives i avsnittene under.
* Metode 1 er enklest å forstå, men den er utsatt for numerisk ustabilitet ettersom den innebærer at man må invertere en stor matrise som nesten er singulær (determinanten har lav verdi).
* Metode 2 har bedre numerisk stabilitet, men har større abstraksjonsnivå og krever derfor dypere forståelse.
 
=== Metode 1: Sett med første ordens koplete differensiallikninger ===
Det to likningene (i) og (ii) kan kombineres til et felles sett med første ordens koplete differensiallikninger med en systemmatrise '''A''':
 
<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 (2''n'' x 2''n'').
 
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 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: <br>
<math> \mathbf{U_0} </math> og <math> \mathbf{I_0} </math> er strøm og spenning ved posisjonen <math> x = 0 </math>, <br>
<math> e^{\mathbf{\Lambda} x} </math> er en diagonal (2''n'' x 2''n'')-matrise med diagonale elementer <math>e^{\lambda_i x}</math>,<br>
<math> \mathbf{V} </math> er (2''n'' x 2''n'') egenvektormatrisen der kolonne nr. ''i'' er egenvektoren <math> \mathbf{v_i} </math> til egenverdien <math> \lambda_i </math> , og <br>
<math> \mathbf{\Phi}(x) </math> er den resulterende løsningsmatrisen med dimensjoner (2''n'' x 2''n'').
 
<figure id="fig:Admittansbeskrivelse_nlinjer">
[[Bilde:Admittansbeskrivelse_n_linjer.png|thumb|none|400px| <caption>Et linjesegment med lengde <math> l </math> har lineære karakteristikker og kan beskrives komplett ved en 2''n''x2''n'' løsningsmatrise <math> \mathbf{\Phi_l} </math> eller admittansmatrise <math> \mathbf{Y_l} </math></caption>]]
</figure>
 
På samme måte som for den enkle linjesløyfa kan man for ei linje med lengde <math> l </math> som vist i <xr id="fig:Admittansbeskrivelse_nlinjer" /> 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,12}^{-1}} \\
\mathbf{Y_{l,22}} & =  - \mathbf{\Phi_{l,22}} \cdot \mathbf{\Phi_{l,12}^{-1}} \end{align}</math>
 
=== Metode 2: Andre ordens differensiallikning, modalanalyse ===
De to likningene (i) og (ii) kan kombineres ved at (i) først deriveres:
 
<math> \frac{\mathrm{d}^2 \mathbf{U}}{\mathrm{d}x^2} = - \mathbf{Z} \cdot \frac{\mathrm{d}\mathbf{I}}{\mathrm{d}x} </math>
 
Deretter settes uttrykket for <math>\frac{\mathrm{d}\mathbf{I}}{\mathrm{d}x}</math> fra likning (ii) inn i likning (i):
 
<math> \frac{\mathrm{d}^2 \mathbf{U}}{\mathrm{d}x^2} = - \mathbf{Z} \cdot \left( - \mathbf{Y} \cdot \mathbf{U} \right) = \mathbf{Z} \cdot \mathbf{Y} \cdot \mathbf{U}</math>
 
Tilsvarende operasjon kan gjøres for strøm, der resultatet blir
 
<math> \frac{\mathrm{d}^2 \mathbf{I}}{\mathrm{d}x^2} = \mathbf{Y} \cdot \mathbf{Z} \cdot \mathbf{I}</math>
 
Løsningen på disse ekvivalente settene av andre ordens differensiallikninger finnes ved såkalt modalanalyse der likningene transformeres til n uavhengige andre ordens differensiallikninger, som kan løses direkte.
 
De to matrisene '''Z·Y''' og '''Y·Z''' har de samme egenverdiene λ<sub>i</sub>, men ulike egenvektormatriser '''T<sub>v</sub>''' og '''T<sub>i</sub>'''. Det kan vises at man kan finne en modal problemformulering ved hjelp av egenverdiene γ<sub>i</sub> egenvektormatrisene '''T<sub>v</sub>''' og '''T<sub>i</sub>''', der differensiallikningen kan løses direkte for hver mode. Løsningene kombineres deretter tilbake til fysiske admittanser i en admittansmatrise ved den inverse transformasjonen.
 
<math> \mathbf{Z_m} = \mathbf{T_v^{-1}} \cdot \mathbf{Z} \cdot \mathbf{T_i} </math>
 
<math> \mathbf{Y_m} = \mathbf{T_i^{-1}} \cdot \mathbf{Y} \cdot \mathbf{T_v} </math>
 
For hver mode finner man en karakteristisk impedans:
 
<math> Z_{cmi} = \sqrt{\frac{Z_m(i,i)}{Y_m(i,i)}} </math>
 
og en forplantningskonstant γ<sub>mi</sub>:
 
<math> \gamma_{mi} = \sqrt{\lambda_i} </math>
 
For hver mode finner man admittansmatrisen for ei linje med lengde l:
 
<math> \mathbf{Y_{m,l}} = \frac{1}{Z_{cmi} \sinh{(\gamma_{mi} \cdot l)}} \begin{pmatrix} \cosh{(\gamma_{mi} \cdot l)} & -1 \\ -1 & \cosh{(\gamma_{mi} \cdot l)} \end{pmatrix} = \begin{pmatrix} y_{mi,11} & y_{mi,12} \\ y_{mi,21} & y_{mi,22} \end{pmatrix} </math>
 
Til slutt transformeres hvert element i impedansmatrisen for alle modene tilbake til fysiske admittanser. Her kan det utnyttes at de to diagonale elementene er like, og at de to ikke-diagonale elementene er like, slik at transformasjonen kun må gjøres for to elementer. De modale elementene ordnes i diagonale matriser <math>Y_{m,ij} </math> der i og j er 1 eller 2.
 
<math>\mathbf{Y_{l,11}} = \mathbf{Y_{l,22}} = \mathbf{T_i} \cdot \mathbf{Y_{m,11}} \cdot \mathbf{T_v^{-1}} </math>
 
<math>\mathbf{Y_{l,12}} = \mathbf{Y_{l,21}} = \mathbf{T_i} \cdot \mathbf{Y_{m,12}} \cdot \mathbf{T_v^{-1}} </math>
 
Linjas endelige admittansmatrise kan til slutt settes sammen til:
 
<math>\mathbf{Y_l} = \begin{pmatrix} \mathbf{Y_{11}} & \mathbf{Y_{12}} \\ \mathbf{Y_{21}} & \mathbf{Y_{22}} \end{pmatrix} </math>
 
= Sugetransformator =
I [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01_Vedlegg_script#Sugetransformator|''Script:'' Sugetransformator]] er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.
 
Ekvivalentskjema for en sugetransformator kan tegnes som vist i <xr id="fig:Sugetransformator_skjema"/>, der y<sub>m</sub> er magnetiseringsadmittansen og y<sub>k</sub> er kortslutningsadmittansen.
 
<figure id="fig:Sugetransformator_skjema">
[[Bilde:Sugetransformator_ekvivalentskjema.png|thumb|none|400px| <caption>Ekvivalentskjema for sugetransformator. TODO: Endre fra "Node 1" og "Node 2" til "Node R" og "Node S"</caption>]]
</figure>
 
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>


<math> U = A \cdot e^{\gamma x} + B \cdot e^{-\gamma x} </math>
For å finne verdiene for parametrene y<sub>m</sub> og y<sub>k</sub> 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 y<sub>m</sub> settes til 0 S, men ved analyse av brudd i returkretsen og ved modellering av metning ved høy belastning og ved kortslutning får denne betydning.


<math> I = - A \cdot \tfrac{1}{Z_0} \cdot e^{\gamma x} + B \cdot \tfrac{1}{Z_0} \cdot e^{-\gamma x} </math>
<figtable id="tab:Data_for_sugetransformatorer">
{| class="wikitable" style="text-align:left"
|+ <caption>Data for sugetransformatorer</caption>
|-
! Parameter !! Symbol !! Enhet !! Sugetransformator <br> 600 A !! Sugetransformator <br> 600 A !! Sugetransformator <br> 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 || 5 759 || 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) || y<sub>m</sub> || mS || 2,18 || 6,55 || 4,21
|-
| Magnetiseringskonduktans (beregnet) || g<sub>m</sub> || mS || 1,07 || 4,64 || 2,84
|-
| Magnetiseringsreaktans (beregnet) || b<sub>m</sub> || mS || 1,90 || 4,61 || 3,11
|-
| Produsent || || || ABB || Møre Trafo || Møre trafo
|-
| Kilde || || || EH-800131 || EH-800207 || EH-800208
|-
! Kortslutningsadmittans !! y<sub>k</sub> !! S !! 13,08 - ''j'' 25,42 !! 12,84 - ''j'' 29,65 !! 14,62 - ''j'' 32,87
|-
! Magnetiseringsadmittans !! y<sub>m</sub> !! mS !! 1,07 + ''j'' 1,90 !! 4,64 + ''j'' 4,61 !! 2,84 + ''j'' 3,11
|}
</figtable>


Tallverdiene for konstantene A og B er gitt av grenseverdibetingelser.
= Autotransformator =
 
I [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01_Vedlegg_script#Autotransformator|''Script:'' Autotransformator]] er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.
 
En autotransformator kan tegnes som vist i <xr id="fig:Autotransformator_skjema"/>.
 
<figure id="fig:Autotransformator_skjema">
[[Bilde:Autotransformator ekvivalentskjema.png|thumb|none|400px| <caption>Ekvivalentskjema for autotransformator.</caption>]]
</figure>
 
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 ''y<sub>k</sub>'' og ''y<sub>m</sub>'' 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).
 
<figtable id="tab:Data_for_autotransformatorer">
{| class="wikitable" style="text-align:left"
|+ <caption>Data for autotransformatorer</caption>
|-
! Parameter !! Symbol !! Enhet !! Autotransformator <br> 5 MVA <br> 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) || g<sub>m</sub> || µS || 2,3875
|-
| Magnetiseringssusceptans <br> antatt lik 2,5 · g<sub>m</sub> || b<sub>m</sub> || µS || 5,9688
|-
| Produsent || || ||  ABB
|-
| Kilde || || || EH-800403
|-
! Kortslutningsadmittans !! y<sub>k</sub> !! S !! 0,6082 - ''j'' 0.9735
|-
! Magnetiseringsadmittans !! y<sub>m</sub> !! µS !! 2,3875 + ''j'' 5,9688
|}
</figtable>
 
= Seksjonering =
I [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01_Vedlegg_script#Seksjonering|''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 <xr id="fig:Seksjonering_skjema"/>.
 
<figure id="fig:Seksjonering_skjema">
[[Bilde:Seksjonering ekvivalentskjema.png|thumb|none|380px| <caption>Ekvivalentskjema for seksjonering.</caption>]]
</figure>
 
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.
 
= 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
 
''i<sub>R</sub>'' er strømmen i jordforbindelsen tilknyttet tilkoplingspunkt ''R''.
 
''y<sub>RE</sub>'' er admittansen mellom tilkoplingspunkt ''R'' og jordpotensialet ''e''.
 
''u<sub>R</sub>'' er spenningen mot jord ved tilkoplingspunkt ''R''.
 
Verdien for ''y<sub>RE</sub>'' 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.
 
= 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.
 
<figure id="fig:Kilde">
[[Bilde:Systembeskrivelse kilde.png|thumb|none|300px| <caption>Kretsskjema for en kilde modellert som en Norton-ekvivalent</caption>]]
</figure>
 
== Kilde ==
I [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01_Vedlegg_script#Kilde|''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.
 
''i<sub>K</sub>'' er den Norton-ekvivalente strømmen til kilden.
''y<sub>K</sub>'' er den Norton-ekvivalente admittansen til kilden.
 
''i<sub>K</sub>'' og ''y<sub>K</sub>'' bestemmes på følgende måte:
 
- ''y<sub>K</sub>'' er kildeadmittansen for den aktuelle kilden. Kildeadmittansen kan velges som en statisk verdi (regulatorkarakteristikk) eller en dynamisk verdi (admittans i maskiner og transformatorer) avhengig av hva som skal beregnes.
 
- Tomgangsspenningen u<sub>0</sub> velges, vanligvis til 16 500 V i norsk banestrømforsyning.
 
- Den Norton-ekvivalente strømmen i<sub>K</sub> settes lik <math> u_0 \cdot y_k </math>. Ved ''y<sub>K</sub>'' = -''j'' 100 S og ''u<sub>0</sub>'' = 16 500 V blir ''i<sub>K<sub>'' = -''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>
 
== Belastning ==
 
I [[Elektrisk_systembeskrivelse_av_kontaktledningsanlegg_ver01_Vedlegg_script#Belastning|''Script:'' Belastning]] er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.
 
Belastninger kan modelleres på flere måter:
# 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.
# Som en admittans ''y<sub>B</sub>''. 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.
# 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.
# 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 ''i<sub>B</sub>'' 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>


= Referanser =
= Referanser =
[1] Edwards, Penney: ''Elementary Linear Algebra'', Pearson, 1987. ISBN 9780132582605.
[1] Edwards, Penney: ''Elementary Linear Algebra'', Pearson, 1987. ISBN 9780132582605.
== Admittanser ==
Admittansen

Siste sideversjon per 28. jun. 2022 kl. 11:47

__NUMBEREDHEADINGS__

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 .

<figtable id="tab:Ledere_oversikt">

Ledere oversikt
Leder Forklaring Nominell spenning
KL Kontaktledning, 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

</figtable>

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.

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. Ved å anta lineær oppførsel neglisjeres ikke-lineære fenomener, i første rekke magnetisk metning og hysterese i magnetiske materialer som jern.

<xr id="fig:LinearModell"/> 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.

<figure id="fig:LinearModell">

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.

</figure>

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

For systemet i <xr id="fig:LinearModell"/> har vektorene og n elementer hver, en for hver terminal, mens admittansmatrisen har (n x n) elementer.

Hver komponent i systemet kan også uttrykkes ved sin admittansmatrise. Systemets admittansmatrise er summen av admittansmatrisene for alle komponentene. Elementer i admittansmatrisen som mellom terminaler som ikke er sammenkoplet, har admittansverdien 0.

Under forutsetning av at en komponent eller et system har minst en jordforbindelse er den tilhørende admittansmatrisen for komponenten eller systemet inverterbar. Det er da mulig å uttrykke alle systemets spenninger som funksjon av påtrykket strøm i alle terminaler, via den inverse admittansmatrisen, som kalles impedansmatrisen. De aller fleste av terminalene blir ikke påtrykket strøm, og i de resterende terminalene er strømmen ofte kjent. Derfor er det ofte av interesse å beregne spenning på bakgrunn av kjente verdier for strøm.

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,
  • temperaturavhengighet for parametre som ohmske motstander,
  • metning og hysterese i transformatorer, og
  • aktivt styrte belastninger, for eksempel tog.

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

Transmisjonslinje

Generelt

Modellering av transmisjonslinje er gjort for:

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 av linjesløyfa. Linjesegmentet har en seriell resistans og reaktans , og en parallell konduktans og susceptans . Her er L, C og henholdsvis linjas induktans, kapasitans og vinkelfrekvens.

<figure id="fig:Telegraflikningene">

Telegraflikningene - Kretsskjema for et linjesegment dx

</figure>

Parametrene kan skrives om slik at impedansen blir og admittansen blir . Parametrene kan finnes ved beregning eller ved måling.

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

(i)

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

(ii)

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

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 til systemmatrisen A er:

der

er linjas forplantningskonstant.

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

der

er linjas karakteristiske impedans.

Forplantningskonstanten 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 er et mål på sammenhengen mellom strøm og spenning i ei transmisjonslinje.

Egenvektorene og ordnes så i en egenvektormatrise :

Egenvektormatrisen har følgende inverse:

Løsningen er gitt av følgende uttrykk:

der:
og er strøm og spenning ved posisjonen ,
er en diagonal 2x2-matrise med diagonale elementer og , og
er den resulterende løsningsmatrisen.

Skrevet fullt ut er løsningsmatrisen:

Løsningen gir en komplett beskrivelse av strøm og spenning i alle punkter langs ei transmisjonslinje, dersom strømmen og spenningen ved posisjon er kjent.

<figure id="fig:Admittansbeskrivelse_2porter">

Et linjesegment med lengde har lineære karakteristikker og kan beskrives komplett ved en 2x2 løsningsmatrise eller admittansmatrise

</figure>

For ei linje av lengde som vist i <xr id="fig:Admittansbeskrivelse_2porter" /> mellom punktene og kan man finne sammenhengen mellom strøm og spenning ved de to linjeendene ved:

der

. Grunnen til det negative fortegnet for 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:

Man kan finne at matrisen er:

Transmisjonslinje med flere parallelle ledere

Beskrivelse av transmisjonslikningene for flerledersytem

I Script: Transmisjonslinje er det angitt en Scilab-funksjon som utfører beregningene i dette avsnittet.
For parametrisering av modellen henvises det til Lenke: Parametre for flerledersystem.

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:

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:

Parametermatrisene , , og er bestemt ut ifra benyttede ledermaterialer, isolasjonsmaterialer, jordsmonn, systemets geometri og driftsfrekvens. Bestemmelse av parametermatrisene er beskrevet i Lenke: Parametre for flerledersystem.

På kompakt matriseform kan disse likningene skrives:

(i)

(ii)

Det er to vanlige tilnærminger for å løse dette likningssettet. De beskrives i avsnittene under.

  • Metode 1 er enklest å forstå, men den er utsatt for numerisk ustabilitet ettersom den innebærer at man må invertere en stor matrise som nesten er singulær (determinanten har lav verdi).
  • Metode 2 har bedre numerisk stabilitet, men har større abstraksjonsnivå og krever derfor dypere forståelse.

Metode 1: Sett med første ordens koplete differensiallikninger

Det to likningene (i) og (ii) kan kombineres til et felles sett med første ordens koplete differensiallikninger med en systemmatrise A:

Her er systemmatrisen en kompleks matrise med dimensjoner (2n x 2n).

Systemmatrisen har 2n egenverdier med 2n tilhørende egenvektorer . Disse må finnes numerisk for aktuelle parametermatriser. Systemets løsning blir på formen:

der:
og er strøm og spenning ved posisjonen ,
er en diagonal (2n x 2n)-matrise med diagonale elementer ,
er (2n x 2n) egenvektormatrisen der kolonne nr. i er egenvektoren til egenverdien , og
er den resulterende løsningsmatrisen med dimensjoner (2n x 2n).

<figure id="fig:Admittansbeskrivelse_nlinjer">

Et linjesegment med lengde har lineære karakteristikker og kan beskrives komplett ved en 2nx2n løsningsmatrise eller admittansmatrise

</figure>

På samme måte som for den enkle linjesløyfa kan man for ei linje med lengde som vist i <xr id="fig:Admittansbeskrivelse_nlinjer" /> mellom nodene og finne sammenhengen mellom strøm og spenning ved de to linjeendende ved:

For å finne admittansmatrisebeskrivelsen deles løsningsmatrisen inn i fire (n x n) submatriser på følgende måte:

En kan da finne at beskrivelsen på admittansmatriseform blir:

der

Metode 2: Andre ordens differensiallikning, modalanalyse

De to likningene (i) og (ii) kan kombineres ved at (i) først deriveres:

Deretter settes uttrykket for fra likning (ii) inn i likning (i):

Tilsvarende operasjon kan gjøres for strøm, der resultatet blir

Løsningen på disse ekvivalente settene av andre ordens differensiallikninger finnes ved såkalt modalanalyse der likningene transformeres til n uavhengige andre ordens differensiallikninger, som kan løses direkte.

De to matrisene Z·Y og Y·Z har de samme egenverdiene λi, men ulike egenvektormatriser Tv og Ti. Det kan vises at man kan finne en modal problemformulering ved hjelp av egenverdiene γi egenvektormatrisene Tv og Ti, der differensiallikningen kan løses direkte for hver mode. Løsningene kombineres deretter tilbake til fysiske admittanser i en admittansmatrise ved den inverse transformasjonen.

For hver mode finner man en karakteristisk impedans:

og en forplantningskonstant γmi:

For hver mode finner man admittansmatrisen for ei linje med lengde l:

Til slutt transformeres hvert element i impedansmatrisen for alle modene tilbake til fysiske admittanser. Her kan det utnyttes at de to diagonale elementene er like, og at de to ikke-diagonale elementene er like, slik at transformasjonen kun må gjøres for to elementer. De modale elementene ordnes i diagonale matriser der i og j er 1 eller 2.

Linjas endelige admittansmatrise kan til slutt settes sammen til:

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 <xr id="fig:Sugetransformator_skjema"/>, der ym er magnetiseringsadmittansen og yk er kortslutningsadmittansen.

<figure id="fig:Sugetransformator_skjema">

Ekvivalentskjema for sugetransformator. TODO: Endre fra "Node 1" og "Node 2" til "Node R" og "Node S"

</figure>

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

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, men ved analyse av brudd i returkretsen og ved modellering av metning ved høy belastning og ved kortslutning får denne betydning.

<figtable id="tab:Data_for_sugetransformatorer">

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 5 759 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 EH-800131 EH-800207 EH-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

</figtable>

Autotransformator

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

En autotransformator kan tegnes som vist i <xr id="fig:Autotransformator_skjema"/>.

<figure id="fig:Autotransformator_skjema">

Ekvivalentskjema for autotransformator.

</figure>

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

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).

<figtable id="tab:Data_for_autotransformatorer">

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 EH-800403
Kortslutningsadmittans yk S 0,6082 - j 0.9735
Magnetiseringsadmittans ym µS 2,3875 + j 5,9688

</figtable>

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 <xr id="fig:Seksjonering_skjema"/>.

<figure id="fig:Seksjonering_skjema">

Ekvivalentskjema for seksjonering.

</figure>

Admittansmatrisen for en slik seriell admittans er:

For en ren seksjonering er motstanden svært stor, og admittansen y kan settes til S.

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:

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:

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.

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.

<figure id="fig:Kilde">

Kretsskjema for en kilde modellert som en Norton-ekvivalent

</figure>

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 er kildeadmittansen for den aktuelle kilden. Kildeadmittansen kan velges som en statisk verdi (regulatorkarakteristikk) eller en dynamisk verdi (admittans i maskiner og transformatorer) avhengig av hva som skal beregnes.

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

- Den Norton-ekvivalente strømmen iK settes lik . Ved yK = -j 100 S og u0 = 16 500 V blir iK = -j 16 500 kA.

Admittansmatrisen for en kilde er:

Strømvektoren for kilden blir:

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:

Referanser

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