Linee elettriche
Si parte dal modello della linea a parametri distribuiti:{{:dinamica_e_controllo:figura65.png|}}
//Fig. 65 circuito equivalente di un concio di linea//
Equazioni di equilibrio
Le equazioni di equilibrio che descrivono il circuito di fig. 65 sono:Equazione alla maglia:
V - (V + dV) = - dV = (r I + l {dI}/{dt}) dx
(1)Equazione al nodo "A":
I - (I + dI) = - dI = (g V + c dV/dT) dx(2)
nella (1) e nella (2), in particolare, i differenziali di I e di V sono espressi in ragione del tempo e dello spazio:
dV = {partial V (x,t)}/{partial x} dx
(3)
dI = {partial I (x,t)}/{partial x} dx
(4)Pertanto la (1) e la (2) possono essere espresse come:
- dV/dx = r I + l {dI}/{dt}
(5)
- dI/dx = g V + c dV/dT(6)
Soluzione delle equazioni differenziali di equilibrio
Per risolvere più agevolmente le suscritte equazioni differenziali, in quanto lineari, si può ricorrere alla trasformata di Laplace in luogo dei metodi risolutivi classici.
V (x,t) doubleright V (x,s)
I (x,t) doubleright I (x,s)
La (5) e la (6) possono pertanto essere riscritte nel modo seguente:
- {partial V _{(x,s)}}/{partial x} = r i + S l i_{(x,s)} = (r + l s) i_{(x,s)}
(7)
- {partial i _{(x,s)}}/{partial x} = g V + S c V_{(x,s)} = (g + c s) V_{(x,s)}
(8)Derivando la (7) e la (8) rispetto a x:
delim{lbrace}{
matrix{2}{1}{
{- {partial^2 i _{(x,s)}}/{partial x^2} = (g + c s) {partial V _{(x,s)}}/{partial x}}
{- {partial^2 V _{(x,s)}}/{partial x^2} = (r + l s) {partial i _{(x,s)}}/{partial x}}
}
}{}
(9)Sostituendo la (7) e la (8), rispettivamente nella prima e nella seconda delle (9), si ottiene:
delim{lbrace}{
matrix{2}{1}{
{- {partial^2 i _{(x,s)}}/{partial x^2} = - (g + c s) (r + l s) i_{(x,s)}}
{- {partial^2 V _{(x,s)}}/{partial x^2} = - (r + l s) (g + c s) V_{(x,s)}}
}
}{}
(10)Si definisce la seguente costante:
k = sqrt{(r + l s) (g + c s)}
(11)La quale, sostituita nella (10), porge:
delim{lbrace}{
matrix{2}{1}{
/{partial x^2} = k^2 i_{(x,s)}}
/{partial x^2} = k^2 V_{(x,s)}}
}
}{}
(12)Integrale generale
Risolvendo il sistema di equazioni differenziali (12) si ottiene la seguente soluzione generale:
V(x,s) = M e^kx + N e^-kx(13)
Di forma analoga è l'espressione della i(x,s) Derivando la (13) si ottiene:
{partial V _{(x,s)}}/{partial x} = M k e^kx - N k e^-kx = - (r + l s) i _{(x,s)}
(14)risolvendo la (14) rispetto a i(x,s) si ottiene:
i _{(x,s)} = - k/{r + l s} (M e^kx - N e^-kx) = - sqrt{{g + c s}/{r + ls}} (M e^kx - N e^-kx)
(15)Si osseva che la grandezza
sqrt{{g + c s}/{r + ls}}
ha le dimensioni di un'impedenza. Pertanto si può porre:
Z_0 = sqrt{{g + c s}/{r + ls}}
(16)
FIXME (rifare i conti perchè forse sono invertiti numeratore e denominatore)[Z0] = [Ω], prende il nome di //impedenza caratteristica della linea//. Sebbene singolarmente le grandezze che la compongono dipendano dalla lunghezza della linea (
R = r l_0 ; L = l l_0 ; G = g l_0 ; C = c l_0), in virtù di semplificazioni danno origine ad una grandezza che __non dipende dalla lunghezza della linea__.
i (x,s) = - 1/Z_0 (M e^kx - N e^-kx)(17)
La costante k, manipolata opportunamente, può essere riscritta nella seguente forma:
k = sqrt{(r + l s) (g + c s)} = sqrt{lc} sqrt{(s + r/l) (s + g/c)} = sqrt{lc} sqrt{(s + alpha + beta) (s + alpha - beta)}
(18)dove:
alpha = 1/2 (r/l + g/c)= costante di attenuazione (19)
beta = 1/2 (r/l - g/c)= costante di fase (20)
Condizioni al contorno
Al fine di calcolare il valore dei parametri M ed N della (17) si considerano le condizioni al contorno, ovvero la tensione e la corrente ad un'estremità della linea:
V(0,s) = V_p (s) = M + N(21)
i(0,s) = I_p (s) = - 1/Z_0 (M - N)(22)
Mettendo a sistema:
delim{[}{
matrix{2}{1}{
{V_p}
{I_p}
}
}{]} = delim{[}{
matrix{2}{2}{
{1}{1}
{-1/Z_0}{1/Z_0}
}
}{]}delim{[}{
matrix{2}{1}{
{M}
{N}
}
}{]}
(23)Risolvendo il sistema (23) si ottiene:
M = {V_p - Z_0 I_p}/2
(24)
N = {V_p + Z_0 I_p}/2
(24)Sostituendo i valori ottenuti nella (17) e nella omessa espressione della v(x,s) si ottiene:
v(x,s) = {V_p - Z_0 I_p}/2 e^kx + {V_p + Z_0 I_p}/2 e^-kx
(25)
i(x,s) = -1/Z_0 ({V_p - Z_0 I_p}/2 e^kx - {V_p + Z_0 I_p}/2 e^-kx)
(26)Applicando la formula di Eulero:
cosh(kx) = {e^kx + e^-kx}/2
(27)la (25) e la (26) possono essere riscritte nel seguente modo:
delim{lbrace}{
matrix{2}{1}{
{v(x,s) = V_p cosh(kx) - Z_0 I_p senh(kx)}
{i(x,s) = -1/Z_0 V_p senh(kx) + I_p cosh(kx)}
}
}{}
(29)Notazione fasoriale
Nella trattazione di linee con perdite in regime sinusoidale, ci si può avvalere della notazione fasoriale: s = jω I queste ipotesi la costante k e l'impedenza caratteristica divengono:
k = sqrt{(r + iota omega l) (g + iota omega c)}
(30)
Z_0 = sqrt{{r + iota omega l}/{g + iota omega c}}
(31)
Linea come doppio bipolo
Occorre anzitutto precisare che la convenzione adottata per i segni nella presente trattazione ha una variante rispetto all'approccio secondo i sistemi elettrici (v fig. 66){{:dinamica_e_controllo:figura66.png|}}
//Fig. 66: differenze tra convenzioni//
Sia "a" la lunghezza della linea:
A = D = cosh(k a)(32)
B = Z_0 senh(k a)(33)
C = 1/Z_0 senh(k a)(34)
Il sistema che descrive la linea come doppio bipolo è pertanto:
delim{lbrace}{
matrix{2}{1}{
{V_p = A V_a + B I_a}
{I_p = C V_a + D I_a}
}
}{}
(35)Essendo le costanti A, B, C e D funzioni iperboliche, si preferisce, al fine di linearizzare, esprimerne lo sviluppo in serie di Taylor:
A = D = cosh(k a) approx 1 + {(k a)^2}/2 + {(k a)^2}/{4 !} = 1 + ZY/2 + {Z^2 Y^2}/{24}
(36)
B = Z_0 senh(k a) approx sqrt{Z/Y} (k a + {(k a)^3}/{3 !} + {(k a)^5}/{5 !})
= sqrt{Z/Y} sqrt{Z Y} (1 + {Z^2 Y^2}/{6} + {Z^4 Y^4}/{120})
(37)
C = 1/Z senh(k a) = Y (1 + {Z^2 Y^2}/6 + {Z^4 Y^4}/120)
(38)Per le linee più comuni (che tipicamente non sono più lunghe di poche centinaia di km) vale l'approssimazione al primo ordine:
delim{lbrace}{
matrix{3}{1}{
{A = D = 1}
{B = Z}
{C = Y}
}
}{}
(39)Linea senza perdite in regime sinusoidale
Si considera, in particolare, il caso di linea senza perdite esercita in regime sinusoidale. In queste ipotesi la (29) può essere riscritta in funzione di jω:
V (x) = V_p cosh(iota omega sqrt{l c} x) - Z_0 I_p senh(iota omega sqrt{lc} x)
(40)
I (x) = - 1/Z_0 V_p senh(iota omega sqrt{lc} x) + I_p cosh(iota omega sqrt{lc} x)
(41)Pertanto le grandezze di arrivo in fuzione di quelle di partenza si ottengono sostituendo la lunghezza "a" al posto della coordinata x nella (40) e nella (41)
V_a = V_p cos(omega sqrt{l c} a) - iota Z_0 I_p sen(omega sqrt{lc} a)
(41)
V_p = V_a cos(omega sqrt{l c} a) + iota Z_0 I_a sen(omega sqrt{lc} a)
(41 bis)
I_a = - 1/Z_0 V_p senh(iota omega sqrt{lc} a) + I_p cosh(iota omega sqrt{lc} a)
(42)Si definisce un nuova grandezza:
lambda = 1/{f sqrt{lc}} = {2 pi}/{omega sqrt{lc}}
(43)Che prende il nome di //lunghezza d'onda della linea//.
definendo:
Z_a = {V_a}/{I_a}
(44)Sostituendo la (44) nella (41 bis) si può esprimere la tensione di partenza come funzione della sola tensione di arrivo:
V_p = V_a (cos({2 pi a}/lambda) + iota {Z_0}/{Z_a} sen({2 pi a}/lambda))
(45)invertendo la (45) si ottiene:
V_a = {V_p}/{cos({2 pi a}/lambda) + iota {Z_0}/{Z_a} sen({2 pi a}/lambda)}
(46)Comportamento a vuoto
Nell'ipotesi in cui la linea sia aperta (Za → ∞) la (46) assume la seguente forma:
V_a = {V_p}/{cos({2 pi a}/{lambda})}
(47)l'andamento della (47) è riportato nel seguente grafico:
{{:dinamica_e_controllo:figura68.png|}}
//Fig. 68: Andamento della tensione di arrivo di una linea a vuoto in funzione della lunghezza//
Esiste pertanto una lunghezza per la quale, a vuoto, la tensione di arrivo è infinita:
lim{a right lambda/4}{V_a} = infty
(48)ovviamente nella realtà cadono alcune delle ipotesi fatte (ad es. linea non dissipativa), per la quale la tensione di arrivo, sebbene sia di valore molto elevato, non è comunque infinita.
Impedenza caratteristica
Se la linea viene invece chiusa su un carico di entità pari alla impedenza caratteristica della linea si ha che:Zl = Z0
V_a = {V_p}/{cos({2 pi a}/lambda) + iota sen({2 pi a}/lambda)} = V_p e^{- iota {2 pi a}/{lambda}}
(49)La tensione di arrivo ha pertanto la stessa ampiezza della tensione di partenza e presenta soltanto uno sfasamento.
Generica condizione di carico
Il luogo dei punti sul piano di Gauss descritto dalla (46) è un'ellisse centrata nell'origine (v. fig. 69){{:dinamica_e_controllo:figura69.png|}}
//Fig. 69: luogo dei punti del denominatore di Va//
Si noti che ad ongi angolo sul grafico di fig. 69 corrisponde una certa lunghezza della linea (fissata una certa lunghezza d'onda).