PLECS 3.4 Online Help

Transmission Line (3ph)

Purpose

3-phase transmission line

Library

Electrical / Passive Components

Description

pict

This component implements a three-phase transmission line. A transmission line is characterized by a uniform distribution of inductances, resistances and neutral capacitances along the line. In multi-wire lines, there are also uniformly distributed mutual inductances and coupling capacitances.

The user has the choice between two different implementations: one with series-connected pi sections of lumped elements and another one with distributed parameters based on traveling wave theory. A stiff solver is recommended for simulating models containing this component.

Pi-Section Line

In many cases, the uniformly distributed parameters of a transmission line can be approximated by a series of pi sections consisting of lumped inductors, capacitors and resistors. The figure below illustrates a single pi section exemplified for a 2-phase line. Depending on the desired fidelity at higher frequencies, the number of series-connected pi sections can be configured.

pict

Let l   be the length of the line and n   the number of pi sections representing the line. The inductance L  , the resistance R  , the neutral capacitance CN   as well as the coupling capacitances Cij   and mutual inductances Lij   of the discrete elements can then be calculated from their per-unit-length counterparts L ′ , R′ , CN′  , C ′ij   and L′ij   using the following equations:

L = lL′,    R = lR ′,    C  = -lC′
    n           n        N   n  N

      l             l
Cij = n-C′ij,   Lij = n-L′ij

It is possible to specify the parameters for each phase individually in order to model asymmetric lines. In this case, the parameters must be provided in vector format. Otherwise, the parameter can be a scalar assigning the same value to all phases.

Distributed Parameter Line

The implementation of a distributed parameter line is based on the traveling wave theory, which describes the time delay phenomenon. This approach is numerically more efficient due to the absence of numerous state variables and should be used in large models.

Modeling asymmetric lines is not supported, therefore all parameters need to be scalar.

Single-Phase Lossless Line

pict

Consider a lossless transmission line with inductance L′  and capacitance C ′  per unit length. At a certain point x   along the total length d  , the relation between the line voltage and current can be described with partial differential equations:

- ∂e-= L ′⋅ ∂i
 ∂x       ∂t

  ∂i       ∂e
- ∂x-= C ′⋅ ∂t-

Since a wave entering the sending end "s" of the line must remain unchanged when it arrives at the receiving end "r" (and vice versa), the following expression is derived:

is = 1-⋅es(t)- Ish
    Z

    1-
ir = Z ⋅er(t)- Irh

where

      1
Ish = Z-⋅er(t- τ)+ ir(t- τ)

I  =  1-⋅e(t- τ)+ i(t- τ)
 rh   Z  s         s

with surge impedance     ∘ L′-
Z =   C′  and travel time       √----
τ = d ⋅ L′C ′  . This model can be represented by a two-port equivalent circuit, where the electrical conditions at port "s" are transferred after a time delay τ   to port "r' via the controlled current source Irh   .

pict

Approximation of Series Resistance

pict

Since the shunt conductance is usually negligible, the series resistance is responsible for the major part of the power losses. Such series resistance can be approximated by three lumped resistors, two of which with the value R-
 4   are placed at both ends of the line while one with the value R-
2   is placed in the middle. R = R ′ ⋅d   is the total series resistance of the line. After aggregation and substitution, the original expression of the two equivalent current source becomes:

I  = 1+-h-⋅( 1-⋅e(t- τ)+ h⋅i(t- τ))+ 1--h-( 1-⋅e(t- τ)+ h⋅i(t- τ))
 sh     2    ZR   r          r          2   ZR  s           s

Irh = 1+-h-⋅( 1-⋅es(t- τ)+ h⋅is(t- τ))+ 1--h-( 1-⋅er(t- τ)+ h⋅ir(t- τ))
       2    ZR                         2   ZR

with          R-
ZR = Z +  4   and       R
h = ZZ-+R4
       4    .

Three-Phase Line

pict

The differential equations of a 3-phase system with vector variables ⃗e = [ea,eb,ec]T   , ⃗i = [ia,ib,ic]T   can be expressed as:

  ∂⃗e-   ′ ∂⃗i   ′ ⃗
- ∂x = L ⋅∂t + R ⋅i

- ∂⃗i-= C ′ ⋅ ∂⃗e
  ∂x       ∂t

Under the assumption of symmetrical phase parameters, the per unit length inductance, capacitance and resistance can be written in matrix form:

    ⌊  L′S  L ′M  L ′M ⌋
L′ = ⌈ L ′M  L′S  L ′M ⌉
      L ′M  L ′M   L′S

     ⌊                                    ⌋
       C′E + 2⋅C ′K     - C ′K       - C ′K
C ′ = ⌈   - C ′K    C′E + 2⋅C ′K     - C ′K   ⌉
          - C ′K       - C ′K    C′E + 2⋅C ′K

     ⌊   ′        ⌋
  ′  ⌈ R    0′  0  ⌉
R  =    0  R   0′
        0   0  R

The presence of off-diagonal elements (mutual inductance and coupling capacitance) in the matrix make it difficult to solve the equation system. However, this can be overcome with the help of modal transformations. If the differential equations are multiplied by a transformation matrix T   on the left side

      ∂⃗e-       ′   -1      ∂⃗i        ′   -1     ⃗
- (T ⋅∂x) = (T ⋅L ⋅T  )⋅(T ⋅∂t)+ (T ⋅R ⋅T   )⋅(T ⋅i)

- (T ⋅ ∂⃗i) = (T ⋅L′ ⋅T -1)⋅(T ⋅ ∂⃗e)
      ∂x                    ∂t

with

    ⌊ 1  1    1 ⌋
T = ⌈ 1  - 2  1 ⌉
      1  1   - 2

the off-diagonal elements of the inductance, capacitance and resistance matrix can be eliminated:

                    ⌊             ⌋
                      L′u   0   0
L ′mod = T ⋅L ′ ⋅T-1 = ⌈ 0  L ′v  0  ⌉
                       0   0  L ′w

                    ⌊   ′         ⌋
 ′          ′  - 1  ⌈ C u  0′  0  ⌉
C mod = T ⋅C ⋅T   =    0  C v  0′
                       0   0   Cw

                    ⌊            ⌋
                       R′  0   0
R ′mod = T ⋅R′ ⋅T -1 = ⌈ 0 R ′  0 ⌉ = R ′
                       0   0  R ′

Thus the original system, in which the three phases are coupled, has been converted to three decoupled systems in the modal domain (denoted as u   , v   , w   ). They can be treated separately in the same way as the single-phase system.

pict

The simulation output in the modal domain should be eventually transformed back into the phase domain via the inverse of the matrix T ′ .

Parameters

Self inductance per unit length
The series self inductance L ′S   per unit length. If the length l   is specified in meters (m) the unit of L ′S   is henries per meter (H/m).

For a pi-section line, the self inductance can be specified individually per phase by providing a 3-element vector of the form L′ = [L ′1 L′2  L′3]  .

Mutual inductance per unit length
The series mutual inductance L ′M   per unit length. If the length l   is specified in meters (m) the unit of L ′M   is henries per meter (H/m).

In a pi-section line, the mutual inductances M ′ij    between the i-th and j-th phase can be specified individually by providing a 3-element vector M ′ = [M1′2 M ′13 M ′23]   containing the upper triangular coupling matrix.

Resistance per unit length
The series resistance   ′
R  per unit length. If the length l   is specified in meters (m) the unit of  ′
R  is ohms per meter (_O_  /m).

For a pi-section line, the parameter can be a vector of the form   ′    ′   ′   ′
R  = [R 1  R2  R3]  .

Neutral capacitance per unit length
The line-to-neutral capacitance C ′N    per unit length. If the length l   is specified in meters (m) the unit of C′N   is farads per meter (F/m).

For a pi-section line, this parameter can be a vector of the form C ′N = [C′N1 C ′N2 C ′N1]  .

Coupling capacitance per unit length
The line-to-line capacitance C′C   per unit length. If the length l   is specified in meters (m) the unit of C ′C   is farads per meter (F/m).

In a pi-section line, the coupling capacitance C′ij   between the i-th and j-th phase can be specified individually by providing a 3-element vector C ′ = [C ′12 C1′3  C′23]   containing the upper triangular coupling matrix.

Length
The length l   of the line. The unit of l   must match the units  ′
LS  ,   ′
L M  ,   ′
R ,  ′
CN   and   ′
CC   are based on.
Number of pi sections
Number of sections used to model the transmission line. The default is 3  . This parameter only affects the pi-section implementation.

Reference

H. Dommel: "Digital Computer Solution of Electromagnetic Transients in Single and Multiple Networks", IEEE Transactions on Power Apparatus and Systems, Vol. PAS88, No. 4, April, 1969