Printer Friendly

A wire electrode inside parallelepipedicaly-shaped ground inhomogeneity: comparison of two solutions.

I. Introduction

The pillars realized in practice (parts of the overhead power network or the lightning protection system) necessarily include a corresponding grounding system. Such system is usually formed of a basic star-, rectangularly- or circularly-shaped electrode [1], [2] connected to an iron armature of the concrete pillar foundation (usually of the shape of parallelepiped [3]), which can be treated as a second part of the grounding system. Official publications as [1] usually neglect influence of the concrete foundation. One simple procedure for approximate modelling of the described influence was proposed in [4] and provides reduction of the analysed grounding system, which includes concrete foundation (ground inhomogeneity), into a problem of a grounding system in the homogeneous ground, which can be solved using well known procedures (in this case Method of Moments (MoM)). The approach is based on approximation of parallelepipedical concrete domain with circle cross-section cylinder and analysis of the current density field. Level of the described influence depends on concrete's and ground's electrical parameters, which can have different values depending on the ground's structure, humidity, etc. The described procedure has been applied in [5]-[7].

In this paper, the resistance of the single wire electrode inside parallelepipedically-shaped domain surrounded by homogenous ground is determined and compared with available results obtained with hybrid boundary elements method (HBEM). This problem was often analysed [8], [9]. The HBEM was recently proposed in [10], based on the idea to substitute boundary segments with total charges placed in free space at the centres of boundary segments. This approach includes using of corresponding Green's function for electric scalar potential and point matching method (PMM) [11] for matching values of potential and boundary condition for normal components of the electric field. It has been applied for electrostatic as well as magnetostatic problems [12]. The results which are used for comparison in this paper are directly adopted from [13] where the method has been applied for analysis of quasi-stationary grounding systems with complex conductivity of multilayered media.

II. The Procedure for Approximating Single Wire Conductor Inside Parallelepiped

The idea to reduce problem of a vertical conductor (which models armature's conductors system) inside the concrete foundation on the problem of wire electrode having equivalent length and cross-section radius placed in homogeneous ground, resulted with the procedure proposed and in details explained in [4]. The basics of this procedure (with additionally included description of the complex function theory application) are briefly presented in the text that follows.

Firstly, using the complex function theory, the system of vertical conductors having a circular cross-section of radius [r.sub.0] and placed on the circle of radius a, can be replaced by a single wire vertical conductor of circular cross-section of radius [a.sub.C], Fig. 1 [14]. The complex function

[] = [C.sub.1]1n[] + [C.sub.2], [] = r[e.sup.j[theta]], (1)

is applied on analysis of the cage system formed from N conductors placed in linear homogeneous media of electrical conductivity a and having leakage current of density per unit length [I.sub.leak]. Now, the complex potential is

[] = -([I.sub.leak] / 2[pi][sigma])[N.summation over (n = 1)]1n([] - [[].sub.n]) + C. (2)

In previous expression are [[].sub.n] = a[e.sup.j(n - 1)[alpha]] and [alpha] = 2[pi]/N, while C is constant depending on the referent level of zero potential. Since is

[[].sup.N] - [a.sup.N] = [N.[product] over (n = 1)][] - [[].sub.n], (3)

it obtains

[] = -([I.sub.leak]/2[pi][sigma])1n([[].sup.N] - [a.sup.N]) + C. (4)

The electrical scalar potental is

[phi] = Re([]) = -[[I.sub.leak]/(4[pi][sigma]N)]1n([r.sup.2N] + [a.sup.2N] -2[(ra).sup.N] cos N[theta]) + C. (5)

Applying previous expression for determining the potential of the single conductor at the point A (Fig. 1) defined with x = a + [r.sub.0], y = 0 = 0, i.e. r = a + [r.sub.0], [theta] = 0, obtains


In (6), with [I.sub.tot] = [NI.sub.leak] the density per unit length of total leakage current from the cage system is labeled. Using condition [r.sub.0] << 2a[pi]/N one can write

[(a + [r.sub.0]).sup.N] [approximately equal to] [a.sup.N] + [Na.sup.N - 1][r.sub.0] (7)

and consequently, the potential of the single cage conductors' system is approximately


If instead of the cage the single cylindrical electrode of cross-section having radii [a.sub.C] , and leakage current of density per unit length [I.sub.tot] is observed, the following expression can be formed

U = -[[I.sub.tot]/ (2[pi][sigma])] 1n([a.sub.C]) + C. (9)

Comparing expressions (8) and (9) for the equivalent radius of the conductors' system from Fig. 1 obtains [4]

[a.sub.C] = [a.sup.N][square root of ([Nr.sub.0]/a.]

In the second step, Fig. 2, adopting that the vertical conductor has a semi-spherically shaped end, applying the equation of continuity and determining the electrical scalar potential in the surrounding ground, the resistance of the vertical electrode is obtained as ((5) in [4])


where [[sigma].sub.1] and [[sigma].sub.C] are electrical conductivities of the ground and the concrete, respectively, and the rest of parameters in (11) can be understood from Fig. 2. Applying the same procedure on a single electrode of equivalent length [l.sub.e] = [K.sub.e][l.sub.c] and of cross-section radius [a.sub.e] = [K.sub.e][a.sub.C] , placed in a homogeneous ground of electrical conductivity [[sigma].sub.1], the electrode's resistance is obtained ((6) in [4])

[] = (1/[K.sub.e])(1/2[pi][[sigma].sub.1][l.sub.C])1n([l.sub.C]/[a.sub.C]), (12)

and consequently, the constant K e can be calculated using the expression ((7) in [4])

[K.sup.-1.sub.e] = [[sigma].sub.1]/[[sigma].sub.C] + (1 - [[sigma].sub.1]/[[sigma].sub.C])[1n([l.sub.C]/b)] / 1n(1 + [1.sub.c]/[a.sub.C]). (13)

The general expression for calculating the electric scalar potential in the vicinity (point defined with field vector [??]) of the single wire electrode of the length [l.sub.e] from Fig. 1 is


In previous expression, with I(s ') is labeled unknown longitudinal current distribution along conductor axis (s' coincide with z axis), K([??], s') = (1/[r.sub.1]) + (1/[r.sub.2])is the kernel, while [r.sub.1] and [r.sub.2] are distances from the conductor element, i.e. its image in the flat mirror from the point where the potential is determined, respectively. The antenna model is used and the longitudinal current is assumed in polynomial form having unknown coefficients [I.sub.m], m = 0,1, ..., M [15]

I(s') = [M.summation over (m = 0)][I.sub.k][(s'/[l.sub.e]).sup.m] (15)

Taking in consideration a fact that quasi-stationary regime is assumed, and after applying partial integration, the expression for calculating the electric scalar potential gets the form


Obviously, there are in total M + 1 unknown coefficients. One equation can be obtained from the condition (Eq. (9c) in [4]):


where [I.sub.leak](s') is leakage current density per unit length.

The rest M equations are formed matching potential value (16) in the points at the conductors' surface [11]


Since is I(s' = 0) = [I.sub.0] = [I.sub.g], the electrode resistance is

[R.sub.g] = U/[I.sub.g] = U/[I.sub.0]. (19)

III. The Application of Hybrid Boundary Elements Method

The application of the hybrid boundary elements method (HBEM) is presented in Fig. 3. Unknown total charges distribution at the boundary surface total charges are modelled with point charges [] n = 1, ... N, placed at the points [] and [], n = 1, ... N. (In this case with r = [square root of [chi square] + [y.sup.2] is labeled radial coordinate of the cylindrical coordinate system).

The longitudinal current distribution is assumed in polynomial form (15). Potential of the system is expressed as [13]


where are:

[k.sub.1n.sup.2] = [4rr.sub.t n] / [(r + [r.sub.t n).sup.2] + [(z + [z.sub.t n]).sup.2], (21)

[k.sub.2n.sup.2] = [4rr.sub.t n] / [(r + [r.sub.t n).sup.2] + [(z + [z.sub.t n]).sup.2] (22)



is a complete elliptic integral of the first kind.

The system of charges and currents in Fig. 3 and potential expression (20) is formed using image theory in the flat mirror and symmetrical features of the structure from Fig. 3. Total M equations (the number of unknown coefficients in leakage current distribution (17)) are obtained matching potential value at the conductors' surface, and N equations (total number of point charges) is formed based on the relation between normal electric field component and surface charge density [13]. Since feeding current can be obtained as


The electrode resistance can be calculated as

[R.sub.g] = 1/[G.sub.g] = U/[I.sub.g] (25)

IV. Numerical Results

Based on the described models and procedures, the single resistance of wire electrode having length [l.sub.C] = 1m and cross-section radius [a.sub.C] = 1m placed in parallelepiped concrete cylinder of electrical conductivity [[sigma].sub.c] and surrounded by ground of electrical conductivity [[sigma].sub.1] = 0.001S/m is determined. The other parameters values are [b.sub.s] = [b.sub.1] = 0.3m. The results obtained with procedure of approximating wire conductor inside parallelepiped (AWCP) are compared with those ones available from [13], which are calculated applying HBEM. Parameter M in (15) is M = 6. Normalized values of the conductance obtained using both described procedures are presented in Table I. It can be noticed that matching between results obtained using described procedures increases with decreasing of electrical conductivity of parallelepiped. The graphs obtained based on the values from Table I are shown in Fig. 4.

V. Conclusions

The resistance of the single wire electrode inside parallelepipedically-shaped domain surrounded with homogeneous ground is determined in the paper. The results obtained applying one recently proposed procedure for approximating wire electrode inside parallelepiped with one equivalent single wire electrode placed in homogeneous ground are presented. In this way obtained values for the resistance of the observed system are compared with recently published results, realized applying hybrid boundary elements method. Agreement of the results improves with decreasing electrical conductivity (i.e. increasing electrical resistivity) of the parallelepipedically-shaped domain. The procedure of approximating wire conductor inside parallelepiped is more "engineering" oriented and simpler for application. Since this kind of a problem is very realistic model of the armature conductors inside concrete pillar foundation, one can conclude that for large values of the concretes' electrical resistivity, the procedure based on application of HBEM gives sufficiently correct results for the resistance of the observed structure. 10.5755/j01.eee.19.5.1845

Manuscript received June 08, 2012; accepted February 16, 2013.

This paper is supported in part by project grant HI44004 (2011-2014) financed by the Ministry of Education and Science of Republic of Serbia.


[1] Grounding of the Power Network Pillars, Technical Recommendation of Electric Power Industry of Serbia, TR No. 9. Belgrade, Serbia, 2010. (in Serbian). [Online]. Available:

[2] J. Nahman, "Digital Calculation of Earthing Systems in Nonuniform Soil", Electrical Engineering, Springer-Verlag, vol. 62, no. 1, pp. 19-24, 1980. [Online]. Available: 10.1007/BF01578113

[3] T. Radojicic, Concrete Construction. Nis, Serbia, Prosveta, 1981. (in Serbian).

[4] N. N. Cvetkovic, P. D. Rancic, "Model for Analyzing a Pillar Grounding System with Cylindrically-shaped Concrete Foundation", Electromagnetics, Taylor & Francis, vol. 29, no. 2, pp. 151-164, 2009. [Online]. Available:

[5] P. D. Rancic, N. N. Cvetkovic, "Pillar Star-shaped Grounding System," Przeglqd Elektroctechniczny (Electrical Review), SIGMANOT Company, vol. 86, no. 5, pp. 100-104, 2010. [Online]. Avaliable: issue.php?lang=0&num=05/2010

[6] P. D. Rancic, N. N. Cvetkovic, "Square shaped electrode as pillar grounding system", in Proc. of the 9th Int. Conf. TELSIKS, October 07-09, 2009, Nis, Serbia, vol. 2 pp. 351-354 [Online]. Available: TELSKS.2009.5339524

[7] N. N. Cvetkovic, D. D. Krstic, "One Model of Vertical Armature Construction inside Concrete Foundation and Ring Electrodes as Pillar Grounding System", International Review on Modelling and Simulations (I.R.E.M.O.S.), Praise Worthy Prize S.r.l., vol. 4, no. 5, pp. 2724-2733, 2011.

[8] F. Dawalibi, D. Mukhedkar, D. Bensted, "Measured and Computed Current Densities in Buried Ground Conductors", IEEE Transactions on Power Apparatus and Systems, vol. 100, no. 8, pp. 4083-4092, 1981. [Online]. Available:

[9] T. Takashima, T. Nakae, R. Ishibahsi, "High frequency characteristics of impedances to ground and field distributions of ground electrodes", IEEE Transaction on Power Apparatus and Systems, vol. 100, no. 4, pp. 1893-1900, 1980. [Online]. Avaliable: 10.1109/TPAS.1981.316532

[10] N. B. Raicevic, S. R. Aleksic, "One Method for Electric Field Determination in the Vicinity of Infinitely Thin Electrode Shells", Engineering Analysis with Boundary Elements, Elsevier, vol. 34, no. 2, pp. 97-104, 2010. [Online]. Available:

[11] R. Harrington, Field Computation by Moment Methods. New York, The Macmillan Company, 1969.

[12] N. B. Raicevic, S. S. Ilic, S. R. Aleksic, "A Hybrid Boundary Element Method for Multilayer Electrostatic and Magnetostatic Problems", Electromagnetics, Taylor & Francis, vol. 30, no. 6, pp. 507-524, 2010. [Online]. Available:

[13] S. S. Ilic, N. B. Raicevic, S. R. Aleksic, "Application of new hybrid boundary element method on grounding systems", in Proc. of the 14th International IGTE Symposium, Graz, Austria, 2010, 20-22. September 2010, pp. 56-61.

[14] J. V. Surutka, Electromagnetics. Belgrade, Serbia, Gradjevinska knjiga, 1989, p. 158. (in Serbian).

[15] B. Popovic, "Polynomial Approximation of Current Along Thin Symmetrical Cylindrical Dipole", in Proc. of the Institution of Electrical Engineers, The Institute of Engineering and Technology, vol. 117, no. 5, pp. 873-878, 1970. [Online]. Available:

N. N. Cvetkovic (1)

(1) University of Nis, Faculty of Electronic Engineering Aleksandra Medvedeva 14, 18000 Nis, Serbia


[sigma]c(S/m)     [G.sub.g]/(4[pi][[sigma].sub.1]lc)

                        AWCP             HBEM

5 * [10.sup.-2]        0.2420           0.3798
[10.sup.-2]            0.2229           0.2864
5 * [10.sup.-3]        0.2028           0.2412
[10.sup.-3]            0.1180           0.1172
5 * [10.sup.-4]        0.0774           0.0726
[10-.sup.4]            0.0207           0.0181
COPYRIGHT 2013 Kaunas University of Technology, Faculty of Telecommunications and Electronics
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2013 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Cvetkovic, N.N.
Publication:Elektronika ir Elektrotechnika
Article Type:Report
Date:May 1, 2013
Previous Article:A particle filter track-before-detect algorithm for multi-radar system.
Next Article:Direct torque control of induction drive using fuzzy controller.

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters