Calculation of Inductance and Capacitance in Power System Transmission Lines Using Finite Element Analysis Method
Calculation of Inductance and Capacitance in Power System Transmission Lines Using Finite Element Analysis Method
Calculation of Inductance and Capacitance in Power System Transmission Lines Using Finite Element Analysis Method
Introduction
In power system, transmission lines transfer electrical
energy between power generation unit and distribution
system. There are various configurations of transmission
lines with different phase conductors, shield wires and kV
ratings. Every transmission line exhibits many electrical
properties, where the most common properties are
inductance and capacitance. The capacitance and
inductance in a transmission line depend on the line
configuration itself. These two electrical properties are
important in the development of transmission line models
used in power system analysis [1].
Nowadays, analytical method has been commonly used
to calculate the inductance and capacitance of many
transmission line configurations. The advantage of using
analytical method is the physical interpretation of
inductance and capacitance equations in transmission line
is clear and well defined. However, analytical method is not
applicable in general, especially for complicated conductor.
Therefore, finite element analysis (FEA) method can be
used to calculate the inductance and capacitance of
transmission lines in general without the need of many
complex calculations [2, 3].
There is a work that has reported on modelling of
capacitance and inductance for multi-conductors system
using finite element analysis (FEA) method [3]. The models
consists of two conductors within two different layers of
dielectric material, two conductors within the same dielectric
layers and three conductors within three different dielectric
layers. The inductance and capacitance per unit conductor
length of multi-conductor transmission lines are related to
each other. A matrix of a transmission line capacitance is
obtained from the FEA model and its inverse matrix is then
used to calculate the inductance of the line. The results
obtained for the capacitance and inductance per unit
conductor length using this method agree well with values
calculated in the literature.
Previous works have considered the energyperturbation and flux-linkage methods to determine
inductance values of electrical machines from a twodimensional (2D) FEA model [4, 5]. The numerical
characteristics and equivalence of the energy-perturbation
and flux-linkage methods have been demonstrated through
calculation of the time-varying apparent inductances of the
stator and rotor field windings of a turbine generator under
transient conditions. It was found that both methods yield
inductance values which agree to each other. However, the
278
L 12 / I 2 10 7 ln p2 p1
C q / V12 2 0 / ln p2 / p1
L 2 10 7 0.5 ln D / r
C 0 / ln D / r
(4)
ln
C 0
re
1 / 4
(7)
A / 0 r J e
(9)
(8)
where
Jc
rc
Ja
ra
-2
1 Am
1
-2
0 Am
1
Air
Conductor
279
p2
L / I B dl
p
1
Sc
(14)
(15)
L
(16)
L12
(11)
12
I ave
c1
c2
c1
Bc1c 2 dl
S c1 J c 2 S c 2
1
1
2 Lab Lac , Lb 2 Lab Lbc , Lc 2 Lac Lbc
where Lab, Lac and Lbc are the inductance due to flux-linkage
between phases a and b, a and c and b and c which are
calculated using eq. (11). The inductance per phase per
unit length has been found to be equal to the average of
inductances from each conductor phase, where
(13)
La Lb Lc
Inserting eq. (12) into eq. (13), the inductance per phase
per unit length can be written as a summation of
inductances between different combination of two
280
1
15
b1c 2
b 2 c1
b 2 c 2 Lc1c 2
b1c1
n
2
Lij , j i
n(n 1) i 1, j 1
where i and j are the i-th and j-th conductor and n is the
total number of conductors. The number of different
combinations of two conductors equals to [n(n-1)]/2. Lij is
calculated using eq. (11). Eq. (17) can be used in general
for n number of conductors in transmission lines.
Proposed Electric Potential FEA Method to Calculate
the Capacitance
The partial differential equation (PDE) which is used to
solve the electric field problem in the model is
(18)
(E )
S J c St
a1
p2
C q / V12 c S E dl
p1
(20)
C12
(21)
V12
c1
p1
E dl
Ca
(23)
Ca Cb Cc
(24)
Ca1
b1
b2
c1
S J c2 S2
1
p2
c2
a2
c
rc
a
ra
-2
Jc Cm
1
-2
0 Cm
1
(27)
n
2
Cij , j i
n(n 1) i 1, j 1
where i and j are the i-th and j-th conductor and Cij is
calculated using eq. (20). Eq. (27) is in similar pattern to
that of the general equation for inductance as eq. (17).
Table 3. Comparison of inductance between two points external to
a single conductor between the proposed and analytical methods
Distance between
Analytical method
Proposed FEA
two points (m)
(mH/km) [1]
method (mH/km)
5
10
15
20
25
30
35
40
45
1.1624
1.3008
1.3818
1.4392
1.4838
1.5203
1.5511
1.5778
1.6013
1.1624
1.3008
1.3818
1.4392
1.4838
1.5203
1.5511
1.5778
1.6013
(25)
C
(26)
1
15
b1c 2
b 2 c1
b 2 c 2 Cc1c 2
b1c1
281
9.7345
8.8496
8.4071
7.9646
7.7434
7.5221
7.5218
7.5215
7.5211
11.0619
9.7345
9.2921
8.8496
8.8493
8.8490
8.8488
8.8485
8.8482
b) Two conductors
Comparison between values obtained using the
proposed FEA and analytical methods of inductance as a
function of distance between centres of two current-carrying
conductors (radius of 0.015 m each) is shown in Table 5.
Similar to a single current-carrying conductor, the
inductance increases when the spacing between the two
conductors is larger. The inductance calculated using the
proposed method agrees with the analytical method. This
again shows that the proposed FEA method can be used to
calculate inductance between two current-carrying
conductors.
Table 5. Comparison of inductance for different distance between
two conductors between the proposed and analytical methods
Distance between Analytical method Proposed FEA method
two conductors (m)
(mH/km) [1]
(mH/km)
10
20
30
40
50
60
70
80
1.3505
1.4891
1.5702
1.6277
1.6723
1.7088
1.7396
1.7663
1.3462
1.4730
1.5340
1.5631
1.6472
1.6642
1.6787
1.6863
4.2668
3.8858
3.6572
3.1238
3.4286
3.4284
3.4282
3.4281
4.2667
3.8857
3.6571
3.0857
3.3905
3.3903
3.3901
3.3900
282
1.3505
1.4891
1.5702
1.6277
1.6723
1.7088
1.7396
1.7663
1.3152
1.4415
1.4987
1.5140
1.4743
1.6373
1.6320
1.6235
4.2593
3.8889
3.6111
3.5185
3.4259
3.4257
3.4256
3.4255
4.3519
3.9815
3.7037
3.6111
3.4260
3.4258
3.4257
3.4256
1.3505
1.4891
1.5702
1.6277
1.6723
1.7088
1.7396
1.7663
Two
Arbitrary-shaped
Proposed FEA
method (mH/km)
1.3031
1.4368
1.5151
1.5706
1.6136
1.6488
1.6785
1.7043
4.2668
3.8858
3.6572
3.1238
3.4286
3.4284
3.4282
3.4281
4.4112
4.0173
3.7810
3.2295
3.5446
3.5444
3.5442
3.5441
[4]
r2
r1
I2
I1
D
[5]
[6]
[7]
[8]
283