Calculation of Inductance and Capacitance in Power System Transmission Lines Using Finite Element Analysis Method

Download as pdf or txt
Download as pdf or txt
You are on page 1of 6
At a glance
Powered by AI
The paper proposes using magnetic flux-linkage FEA method to calculate inductance and electric potential FEA method to calculate capacitance for transmission lines with various conductor configurations.

Analytical method provides clear physical interpretation but is not applicable to complex conductor shapes, while FEA method is applicable in general without complex calculations but depends on conductor shape. FEA results agree with analytical method when assumptions are valid.

The proposed FEA method calculates inductance using magnetic flux-linkage and capacitance using electric potential. General equations are deduced from the FEA model as the summation of inductances/capacitances between conductor pairs divided by the number of pairs.

Hazlee Azil ILLIAS1,2, Ab Halim ABU BAKAR2, Hazlie MOKHLIS1,2, Syahirah Abd HALIM1

Faculty of Engineering, University of Malaya (1), UMPEDAC (2)

Calculation of inductance and capacitance in power system


transmission lines using finite element analysis method
Abstract. In power system, every transmission line exhibits many electrical properties. Analytical method has been widely used in determination of
the inductance and capacitance for various transmission line configurations. However, it not applicable in general especially for complicated
conductor arrangements. Therefore, in this work, a magnetic flux-linkage finite element analysis (FEA) method has been proposed to calculate the
inductance in various conductor arrangements in different transmission line configurations while the capacitance is calculated using electric potential
FEA method.
Streszczenie. W artykule opisano metod okrelania indukcyjnoci I pojemnoci linii przesyowej wykorzystujcy analiz strumienia rozproszonego
metoda elementu skoczonego. Rozpatrzono przykady rnych konfiguracji linii. (Obliczenia indukcyjnoci i pojemnoci linii przesyowej z
wykorzystaniem metody elementu skoczonego)

Keywords: Capacitance; inductance; transmission line modelling; finite element methods.


Sowa kluczowe: linia przesytowa, okrelanie pojemnoci, okrelanie indukcyjnoci.

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

flux-linkage method uses less computation spaces and is


more robust than the energy perturbation method.
Therefore, the flux-linkage method has been suggested as
the preferred method for calculation of 2D inductance of
electrical machines.
In this work, the calculations of the inductance and
capacitance of various conductor arrangements in
transmission lines have been performed by the proposed
magnetic flux-linkage and electric potential FEA method.
The calculated inductance and capacitance using the
proposed method have been compared with the analytical
method to validate the proposed method. In the end of the
work, general equations for inductance and capacitance in
transmission lines have been deduced based on the results
obtained through the proposed method. Compared to the
analytical method, the proposed method is not only fast and
simple but can also be used for general conductor
arrangements in power transmission line.
Calculation of inductance and capacitance using
analytical method
General assumptions made for calculating inductance of
transmission lines are:
a) The cross section of transmission line is uniform
b) The current and charge densities are uniform over the
entire length of conductors
c) The conductivity, permittivity and permeability over
the entire length of conductors are constant
For a single cylindrical conductor carrying a current I or
charge q across the cross sectional area, as shown in Fig.
1(a), the inductance, L and capacitance, C between two
points, p1 and p2 in the magnetic and electric field regions
around the conductor are obtained by [1]
(1)
(2)

L 12 / I 2 10 7 ln p2 p1

C q / V12 2 0 / ln p2 / p1

where 12 and V12 are the magnetic flux-linkage and electric


potential between the two points respectively.
A single-phase two-wire conductor of radius r carrying a
current, I1 or charge q1 and the return current I1 or charge
q1, separated by a distance, D is shown in Fig. 1(b). The
inductance, L and capacitance, C per phase per unit length
are calculated using [1, 6]
(3)

L 2 10 7 0.5 ln D / r

PRZEGLD ELEKTROTECHNICZNY (Electrical Review), ISSN 0033-2097, R. 88 NR 10a/2012

C 0 / ln D / r

(4)

Fig.1. (a) Current-carrying conductor perpendicular to the plane, (b)


single-phase and (c) three phase-transmission line

A three-phase transmission lines is shown in Fig. 1(c).


The arrangement consists of three conductors of phases a,
b and c, each having a radius r and separated by a distance
Dab, Dac and Dbc between each of them. The inductance, L
and capacitance, C per phase per unit length are defined as
[1, 6-8]
(5)
(6)

ln

L 2 10 7 ln 3 Dab Dac Dbc

C 0

re
1 / 4

Dab Dac Dbc r

(7)

L 2 10 7 ln 3 Dab ' Dac ' Dbc ' re 1/ 4 6 Da1a 2 Db1b 2 Dc1c 2

C 0 ln 3 Dab ' Dac ' Dbc '

r 6 Da1a 2 Db1b 2 Dc1c 2

The partial differential equation (PDE) which is used to


solve the magnetic field problem in the model is

A / 0 r J e

(9)

A three-phase double circuit lines is shown in Fig. 2,


where each conductor carrying a current of relative phase
position a1b1c1-c2b2a2. The inductance per phase unit length
for this transmission line configuration are defined as [1, 6]

(8)

Fig.2. Three-phase double-circuit lines; (a) a transmission line


tower and (b) model of the conductors

where

where 0 is the vacuum permeability, r is the relative


permeability (in Hm-1), A is the magnetic potential or the
-1
-2
flux-linkage (Wbm ) and Je is the current density (Am ).
Free triangular mesh elements are used since it allows
more freedom of meshing [9]. For a higher accuracy in the
calculation, the mesh is refined near the conductor region.
The subdomain settings of the model are detailed in Table
1.
Table I. Subdomain settings in the model geometry for magnetic
field problem
Symbol Value
Subdomain
Parameter
Conductor
Air

External current density


Relative permeability
External current density
Relative permeability

Jc
rc
Ja
ra

-2

1 Am
1
-2
0 Am
1

Dab ' 4 Da1b1 Da1b2 Da 2 b1 Da 2 b2


Dbc ' 4 Db1c1 Db1c 2 Db2 c1 Db2 c 2
Dac ' 4 Da1c1 Da1c 2 Da 2 c1 Da 2 c 2
Proposed Magnetic Flux-Linkage FEA Method to
Calculate the Inductance
The proposed calculation of inductance in transmission
line using FEA method in this work is solved with COMSOL
Multiphysics software. Fig. 3 shows a two-dimensional (2D)
model geometry that has been developed using finite
element analysis (FEA) method to solve the magnetic field
problem in transmission lines. The model geometry consists
of a conductor (radius r) carrying a current perpendicular to
the plane enclosed by surrounding air.

Air
Conductor

Fig.3. 2D model geometry of a single conductor

From the FEA model, the inductance, L between two


points, p1 and p2 external to a conductor is obtained using

PRZEGLD ELEKTROTECHNICZNY (Electrical Review), ISSN 0033-2097, R. 88 NR 10a/2012

279

flux-linkage, divided by the current I through the


conductor,
(10)

p2
L / I B dl
p
1

Sc

where Jc is the current density across the cross-sectional


area, Sc of the conductor. The flux linkage, is obtained by
integration of the magnetic flux density, B along a line
connecting the two points.

conductors divided by the number of different combination


of two conductors,

(14)

Lab Lac Lbc

To further verify eq. (14), the similar method is tested on


a three-phase double-circuit transmission line, as shown in
Fig. 5. The line between two phase conductors is used to
calculate the flux-linkage between the two conductors. The
inductance from each conductor phase is equal to

La1b1 La1b 2 La1c1 La1c 2 La1a 2


La 2 5 La 2b1 La 2b 2 La 2c1 La 2 c 2 La 2 a1
Lb1 1 5 Lb1a1 Lb1b 2 Lb1c1 Lb1c 2 Lb1a 2
Lb 2 1 5 Lb 2 a1 Lb 2 a 2 Lb 2c1 Lb 2 c 2 Lb 2b1
Lc1 1 5 Lc1b1 Lc1b 2 Lc1a1 Lc1c 2 Lc1a 2
Lc 2 1 5 Lc 2b1 Lc 2b 2 Lc 2c1 Lc 2 a 2 Lc 2 a1
La1

(15)

Deriving the inductance per phase per unit length as the


average of inductances from each conductor phase, again it
is found that the inductance per phase per unit length
equals to the summation of inductances between different
combinations of two conductors, divided by the number of
different combination of two conductors, where

L
(16)

Fig.4. 2D model geometry of (a) two and (b) three conductors

Fig. 4a shows 2D FEA geometry of two conductors


carrying current density Jc and -Jc while Fig. 4b shows three
0
conductors carrying current densities of Jc, Jc120 and
0
Jc240 . The flux linkage between two conductors, is
obtained by integration of the magnetic flux density, B along
a line connecting the two conductors. For Fig. 4a, the
inductance per unit length has been found to be equal to

L12

(11)

12

I ave

c1

c2

c1

Bc1c 2 dl

S c1 J c 2 S c 2

where Iave is the average current of the two conductors, c1


and c2 and 12 is the integration of the magnetic field, Bc1c2
along the line between the two conductors. It is multiplied
by since 12 is the net flux-linkage due to two conductors.
For three conductors or three-phase transmission line,
as shown in Fig. 4b, the inductance of phases a, b and c in
the FEA model can be calculated using
(12) La

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

La1 La 2 Lb1 Lb2 Lc1 Lc 2

1
15

La1b1 La1b 2 La1c1 La1c 2 La1a 2

La 2b1 La 2b 2 La 2c1 La 2c 2 Lb1b 2


L L L L

b1c 2
b 2 c1
b 2 c 2 Lc1c 2
b1c1

Therefore, using flux-linkage method in finite element


analysis, the inductance per phase per unit length can be
deduced in general as the average of the inductances
between different combinations of two conductors, Lij, or
(17)

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 )

where is the permittivity, E is the electric field and is the


charge density.
The charge in a conductor per unit length and time can
be related to the current density, Jc on the cross-sectional
area of the conductor, S by
(19)

S J c St

For per unit length and time, equals to Jc in eq. (19).


Table II details the subdomain setting to compute the
electric field in the model.
The capacitance, C between two points, p1 and p2
external to a conductor is obtained using the charge, q

PRZEGLD ELEKTROTECHNICZNY (Electrical Review), ISSN 0033-2097, R. 88 NR 10a/2012

divided by the electric potential, V12 between the two points


from the FEA model,

a1

p2
C q / V12 c S E dl
p1

(20)

where E is the electric field. V12 is obtained by integration of


E along the line between the two points.
For two conductors in parallel to each other as in Fig.
4a, the capacitance per unit length is calculated as the
average of charge in conductors 1 and 2 divided by the
electric potential between the two conductors,

C12

(21)

V12

c1

p1

E dl

where Jc1 and Jc2 are the current density in conductors 1


and 2.
In Fig. 4b, for a three-phase transmission line, the
capacitance of phases a, b and c in the FEA model can be
calculated using
(22)

Ca

Cab Cac , Cb 1 2 Cab Cbc , Cc 1 2 Cac Cbc

where Cab, Cac and Cbc are the capacitance between


phases a and b, a and c and b and c which are calculated
using (21). The capacitance per phase per unit length has
been found to be equal to the average of capacitances from
each conductor phase, where

(23)

Ca Cb Cc

By inserting eq. (22) into eq. (23), the capacitance per


phase per unit length equals to a summation of
capacitances between different combination of two
conductors divided by the number of different combination
of two conductors, where

(24)

Cab Cac Cbc

Similar to the inductance calculation, eq. (24) is further


verified on a three-phase double-circuit transmission line,
as shown in Fig. 5. The capacitance from each conductor
phase is equal to

Ca1b1 Ca1b 2 Ca1c1 Ca1c 2 Ca1a 2


Ca 2 5 Ca 2b1 Ca 2b 2 Ca 2c1 Ca 2c 2 Ca 2 a1
Cb1 1 5 Cb1a1 Cb1b 2 Cb1c1 Cb1c 2 Cb1a 2
Cb 2 1 5 Cb 2 a1 Cb 2 a 2 Cb 2c1 Cb 2c 2 Cb 2b1
Cc1 1 5 Cc1b1 Cc1b 2 Cc1a1 Cc1c 2 Cc1a 2
Cc 2 1 5 Cc 2b1 Cc 2b 2 Cc 2c1 Cc 2 a 2 Cc 2 a1

Ca1

b1

b2
c1

S J c2 S2

1
p2

c2

a2

Fig.5. Three-phase double-circuit transmission lines


Table 2. Subdomain Settings in the Model Geometry for Electric
Field Problem
Symbol Value
Subdomain
Parameter
Conductor
Air

Space charge density


Relative permittivity
Space charge density
Relative permittivity

c
rc
a
ra

-2

Jc Cm
1
-2
0 Cm
1

Therefore, the general equation of capacitance per


phase per unit length in transmission lines for n conductors
using FEA method can be written as

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

The capacitance per phase per unit length has been


found to be equals to a summation of capacitances
between different combinations of two conductors, divided
by the number of different combination of two conductors,
where

C
(26)

Ca1 Ca 2 Cb1 Cb 2 Cc1 Cc 2

1
15

Ca1b1 Ca1b 2 Ca1c1 Ca1c 2 Ca1a 2

Ca 2b1 Ca 2b 2 Ca 2c1 Ca 2c 2 Cb1b 2


C C C C

b1c 2
b 2 c1
b 2 c 2 Cc1c 2
b1c1

Results and Discussion


a) Single conductor
Table 3 compares the inductance between two points
external to a current-carrying conductor (radius of 0.015 m)
obtained using the proposed flux-linkage FEA method and
the analytical method. The inductance is higher when the
distance between two points external to the conductor is
larger. This is due to the inductance is dependent on fluxlinkage between the two points, where a further distance
yields higher inductance. The inductance is approaching
constant values for a longer distance between the two
points because the magnetic fields are nearly zero. For the
inductance calculation, there is a good agreement between
results obtained using the proposed flux-linkage FEA
method and the analytical method. Therefore, this shows
that the proposed FEA method can be used to calculate
inductance between two points external to a currentcarrying conductor.
The capacitance between two points external to a
current-carrying conductor (radius of 0.015 m) obtained

PRZEGLD ELEKTROTECHNICZNY (Electrical Review), ISSN 0033-2097, R. 88 NR 10a/2012

281

using the proposed electric potential FEA method and the


analytical method is shown in Table 4. The capacitance
between the two points decreases with their spacing
because the potential different between them increases.
From this table, the results obtained using the proposed
FEA method and analytical method is in good agreement.
Hence, the proposed FEA method can also be used to
calculate the capacitance between two points external to a
current-carrying conductor.
Table 4. Comparison of capacitance between two points external to
a single conductor between the proposed and analytical methods
Distance between Analytical method Proposed FEA method
two points (m)
(pF/m) [1]
(pF/m)
5
10
15
20
25
30
35
40
45

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

Table 6. Comparison of capacitance for different distance between


two conductors between the proposed and analytical methods
Distance between Analytical method
Proposed FEA
two conductors (m)
(pF/m) [1]
method (pF/m)
10
20
30
40
50
60
70
80

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

c) Three conductors or three-phase transmission line


Table 6 shows comparison between capacitance
obtained from the proposed FEA and analytical methods for
different distance between centres of two current-carrying
conductors (radius of 0.015 m each). It can be seen that
both methods are in good agreement. Therefore, the
proposed electric potential FEA method can be used to

282

calculate the capacitance between two current-carrying


conductors in a transmission line.
The proposed method has also been used to calculate
the inductance of a three-phase transmission line,
consisting of three conductors (radius of 0.015 m each)
arranged symmetrically. Referring to Table 7, the
inductance values calculated using the proposed FEA
method seems to be in agreement with the values
calculated using the analytical method. The behaviour of
inductance as a function of spacing between three
conductors is similar to a two-conductor transmission line.
The results of capacitance per phase per unit length
using the proposed FEA and analytical methods for a threephase transmission line are shown in Table 8. This again
shows the calculated capacitance values using the
proposed FEA method are in good agreement with the
analytical method.
Table 7. Comparison of inductance for different distance between
three phase conductors between the proposed and analytical
methods
Distance between Analytical method
Proposed FEA
two conductors (m)
(mH/km) [1]
method (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.3152
1.4415
1.4987
1.5140
1.4743
1.6373
1.6320
1.6235

Table 8. Comparison of capacitance for different distance between


three phase conductors between the proposed and analytical
methods
Distance between Analytical method
Proposed FEA
two conductors (m)
(pF/m) [1]
method (pF/m)
10
20
30
40
50
60
70
80

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

Table 9. Calculated Inductance for Three-Phase Double-Circuit


Lines
Proposed
Parameters
Analytical
FEA
method
method
y1,
x1, m x2, m x3, m
y2, m (mH/km) [1] (mH/km)
m
5.5
8.25
6.25
3.5
3.25
0.5692
0.5687
8.25 12.375 9.375 5.25 4.875
0.6097
0.6076
11.0 16.5
12.5
7.0
6.5
0.6385
0.6321
13.75 20.625 15.625 8.75 8.125
0.6608
0.6452
Table 10. Calculated capacitance for three-phase double-circuit
lines
Propose
Parameters
Analytical
d FEA
method
method
x1, m x2, m x3, m
y1,
y2, m
(pF/m) [1]
(pF/m)
m
5.5
8.25
6.25
3.5
3.25
10.219
9.8289
8.25 12.375 9.375 5.25 4.875
9.5102
9.1442
11.0 16.5
12.5
7.0
6.5
9.0642
8.6711
13.75 20.625 15.625 8.75 8.125
8.7461
8.2617

d) Three-phase double-circuit transmission line


For three-phase double-circuit lines, the parameters that
have been used in Fig. 2 and the results calculated using
the proposed FEA and the analytical method are shown in
Table 9. The conductor radius, r that has been used in the

PRZEGLD ELEKTROTECHNICZNY (Electrical Review), ISSN 0033-2097, R. 88 NR 10a/2012

calculation is 0.015 m. It has been found that the


inductance values calculated using all methods are within
reasonable agreement to each other. Therefore, the
proposed FEA method can be said reasonable in
calculating the inductance of three-phase double-circuit
transmission lines.
Table 10 shows the calculated capacitance using the
proposed FEA and analytical methods for a three-phase
double-circuit line with model parameters similar to those
that have been used in Table 9. The results clearly show
that the values calculated using both methods are within
reasonable agreement to each other. Thus, the proposed
FEA method is acceptable for capacitance calculation of a
three-phase double-circuit transmission line.
Table 11. Calculated Inductance for
Conductors
Distance between Analytical method
two conductors (m)
(mH/km) [1]
10
20
30
40
50
60
70
80

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

Table 12. Calculated Capacitance for Two Arbitrary-shaped


Conductors
Distance between Analytical method
Proposed FEA
two conductors (m)
(pF/m) [1]
method (pF/m)
10
20
30
40
50
60
70
80

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

the analytical method assumes that the conductor is a


cylindrical, which simplifies the calculation. Whereas, the
inductance and capacitance calculated using the FEA
method depend on the shape of the conductor, which
influences the magnetic and electric field distributions
around the conductor.
Conclusions
Calculations of the inductance using magnetic fluxlinkage and the capacitance using electric potential in finite
element analysis (FEA) method for various transmission
line configurations have been proposed in this work. The
inductance and capacitance calculated using the proposed
FEA method has been validated through comparison with
the analytical method. It has been found that both results
from the proposed FEA and analytical methods are within
reasonable agreement to each other. From the FEA model,
general equations of inductance and capacitance per phase
per unit length for multi-conductors in power system
transmission lines have been deduced as the summation of
the inductances and capacitances between two conductors
divided by the number of different combinations of two
conductors respectively. The calculations of inductance and
capacitance using the proposed FEA method are simple,
straight forward and is applicable in general, especially for
complicated conductor arrangements and arbitrary
conductor shapes.
This work was supported by University of Malaya, Malaysia
through UMRG Research Grant (Grant no. RG135/11AET).
REFERENCES
[1]
[2]
[3]

[4]

r2

r1

I2

I1
D

[5]

[6]
[7]

Fig.6. Two conductors of arbitrary shape parallel to each other

e) Inductance and capacitance for arbitrary conductor


shape
Fig. 6 shows two conductors of arbitrary shape with a
cross sectional area of A1 and A2, carrying a current, I1 and
I2, separated by a distance D. For the calculation of
inductance and capacitance, the analytical method
assumes that the conductor is a cylindrical shape with a
radius of r, as shown in eqs. (3) and (4). Thus, this is the
limitation of the analytical method as it is not applicable to
general conductor shape in transmission lines [10]. In Fig.
6, if the radius of each conductor is taken as the average
radius (here, the average radius equals to 0.015 m), if r1 =
r2 and I2 = -I1, the calculated inductance and capacitance
using the proposed FEA and analytical method are shown
in Tables 11 and 12. The parameters used are the same as
in Tables 1 and 2. From these tables, it is found that the
results are slightly different between the two methods,
unlike those obtained in Tables 5 and 6. This may be due to

[8]

H. Saadat, Power System Analysis, 2nd ed.: McGraw-Hill


(2004).
S. Cristina, M. Feliziani, A finite element technique for
multiconductor cable parameters calculation, IEEE T. Magn. 25
(1989) 2986-2988.
S. M. Musa, M. N. O. Sadiku, Application of the finite element
method in calculating the capacitance and inductance of
multiconductor transmission lines, IEEE Southeastcon. (2008)
300-304.
R. Escarela-Perez, E. Campero-Littlewood, M. A. ArjonaLopez, A. Laureano-Cruces, Comparison of two techniques for
two-dimensional finite-element inductance computation of
electrical machines, IEE Proc. - Electric Power Applications.
152 (2005) 855-861.
E. Melgoza, R. Escarela-Perez, J. Alvarez-Ramirez, FiniteElement Inductance Computation in 2-D Eddy-Current
Systems Using Sensitivity Analysis, IEEE T. Energy Conver.
25 (2010) 690-697.
J. D. Glover, M. S. Sarma, T. J. Overbye, Power System
Analysis and Design, 4th ed.: Cengage Learning, USA (2004).
D. Labridis, V. Hatziathanassiou, Finite element computation
of field, forces and inductances in underground SF6 insulated
cables using a coupled magneto-thermal formulation, IEEE T.
Magn. 30 (1994) 1407-1415.
R. Lucas, S. Talukdar, Advances in Finite Element
Techniques for Calculating Cable Resistances and
Inductances, IEEE T. Power Ap. Syst. PAS-97 (1978) 875883.

Authors: dr. Hazlee Azil Illias, Department of Electrical


Engineering, Faculty of Engineering, University of Malaya, 50603
Kuala Lumpur, Malaysia, E-mail: [email protected]; dr. Ab Halim
Abu Bakar, University of Malaya Power Energy Dedicated
Advanced Centre (UMPEDAC), University of Malaya, 50603 Kuala
Lumpur, Malaysia, E-mail: [email protected]; dr. Hazlie Mokhlis,
Department of Electrical Engineering, Faculty of Engineering,
University of Malaya, 50603 Kuala Lumpur, Malaysia, E-mail:
[email protected]; ms. Syahirah Abd Halim, Department of
Electrical Engineering, Faculty of Engineering, University of
Malaya,
50603
Kuala
Lumpur,
Malaysia,
E-mail:
[email protected].

PRZEGLD ELEKTROTECHNICZNY (Electrical Review), ISSN 0033-2097, R. 88 NR 10a/2012

283

You might also like