10.1515 - Scjme 2017 0005
10.1515 - Scjme 2017 0005
10.1515 - Scjme 2017 0005
VOL 67 (2017), NO 1, 45 - 54
Abstract: In this contribution, the analysis of ice-shedding from Aluminium Conductor Steel Reinforced (ACSR)
power lines is presented. The impact of the icing position on the overhead power lines, the resulting jump height,
and impact on attachment tension points after ice-shedding is examined. In the numerical simulations the effective
material properties of the ACSR conductor is calculated using the homogenisation method. Numerical analysis of
one power line and double-bundle power lines with icing over the whole range or only on certain sections of single
and double-bundle power lines are performed.
KEYWORDS: ACSR power line, ice-shedding, finite element method, transient analysis
1 Introduction
In cold regions, atmospheric icing is one of the major external loads threatening the
reliability and mechanical integrity of overhead power lines. Ice-shedding form an iced
overhead power line (eq. temperature rise) cause a vertical jump of the power lines and may
lead to flashover if the phase-to-phase or phase-to-tower distance is smaller than the tolerable
insulation distance. Therefore, it is necessary to determine the maximum jump height of an
overhead power line after ice-shedding and to provide a reference for the design of the overhead
power lines [1]. Ice-shedding can also cause dangerous vibrations in power lines, which can
result in mechanical damage of the power line and power line pylons. Therefore, the transient
analysis of ice-shedding is necessary.
Ice-shedding has been investigated with experimental, numerical and theoretical methods by
many authors. Very approximate practical models have been suggested as early as the 1940s
[2]. With the improvement of computational mechanics, numerical simulation methods (over
all the finite element method - FEM) were used to study ice-shedding from the power lines eq.
in [3] the simulation of the dynamic responses of transmission lines with different parameters
after ice-shedding by means of FEM. In [4, 5] the dynamic behavior of bundle conductors and
five-span line section after ice-shedding is numerically simulated. In [1], a new theoretical
method to calculate the jump height of the overhead power line after ice-shedding is presented.
In this paper, the results of the transient analyses of ice-shedding from the ACSR power lines
are presented. The results are calculated using a commercial finite element software ANSYS.
DOI: 10.1515/scjme-2017-0005, Print ISSN 0039-2472, On-line ISSN 2450-5471 2017 SjF STU Bratislava
This gives the power lines an overall high tensile strength. So, the material of the ACSR
conductor is inhomogeneous (Fig. 1) and it is complicated to model as a real conductor.
Therefore, simplified models obtained by the homogenization of material properties are used.
Methods based on homogenization theory (e.g. the mixture rules [6]; self-consistent methods
[7]) have been designed and successfully applied to determine the effective material properties
of heterogeneous materials from the corresponding material behaviour of the constituents (and
of the interfaces between them) and from the geometrical arrangement of the phases. In this
context, the microstructure of the material under consideration is basically taken into account
by the Representative Volume Element (RVE).
The homogenization techniques derived at our department for modelling Functionally
Graded Material (FGM) [8, 9] can also be used for homogenization of the ACSR power lines.
In case of the power line, the material properties vary layer-wise in the radial direction. The
effective elasticity moduli for axial, shear, flexural and torsional loading are calculated from
the condition, that the relevant stiffness of the cross-section compared with real construction
(Fig. 2) is equal to the stiffness of the homogenized cross-section.
z i Ri cos i (2)
Then the quadratic moment of the ith wire cross-sectional area Ai d i2 / 4 according to the
y axis. With respect to the z axis the position can be calculated by equations [8, 10]:
d i 4 d i 2 d i 4 d i 2
I yi zi I zi yi
2 2
(3)
64 4 64 4
46 2017 SjF STU Bratislava Volume 67, No. 1 (2017)
The polar moment of the wire cross-sectional area to the origin of the coordinate system x, y is:
I pi I yi I zi (4)
We assume that the maximal and minimal elasticity modulus for lateral and transversal
bending can be calculated by equations [8, 11]:
n
E I i yi
E LMmax i 1
M Hy zH
E L max n
(5)
I i 1
yi
M Hy
E L min E LMmin
zH
64
2
nFe d Fe E Fe n Al d Al2 E Al (6)
where nFe is the number of steel wires and nAl is the number of aluminium wires. The maximal
elasticity modulus represents the case, when all the wires are fixed together (e.g. after several
years of lifetime), and the minimal elasticity modulus represents the case, when wires can slide
over each other. In practice the effective elasticity modulus for lateral and transversal bending
is assumed as the average value for maximal and minimal elasticity moduli [11]:
E y E L min
M Hy M H
L max
M yH
E L (7)
2
E LMmax
zH
E LMmin
zH
E M zH
L (8)
2
We also assume that the effective Poisson ratio is equal to [8]:
n
i Ai
NH
L i 1
n
(9)
A
i 1
i
where i [-] is the Poisson ratio of the ith wire. The effective mass density for axial beam
vibration is [4]:
n
A k i i i
NH
L i 1
(10)
A
where i [kgm-3] is the material density of each wire and ki [%] is the weight of each material
increment due to stranding.
The nominal design icing load ID is calculated using the formula [11]:
I D K h I R I (11)
where Kh [-] is the height factor of the icing load, γI [-] is the partial icing load factor. In the
calculations of the mechanical stresses, it is needed to calculate the weather load factor z [-]
(include icing), which is calculated using equation [11]:
q1 I D
z (12)
q1
where q1 [Nm-1] is the gravity load of the power line. In practice, calculations mainly utilize the
icing area in N2, so in simulating ice-shedding the icing load is calculated for this case.
When the ice formed on the power line drops off due to changing weather conditions
(temperature, wind) or man-made shocks, it can cause ice-shedding (Fig. 3). The resulting drop
causes the power line to swing upward. According to Fig. 3 it holds:
Δy ymz ymax (13)
where ymax [m] and ymz [m] are maximal deflection with and without icing, respectively. The
maximal jump height ysh [m] of the power lines can be calculated as [11]:
ysh 2Δy (14)
4 Numerical experiment
Transient analyses of the single and double-bundle power lines according to Fig. 4 have been
considered in this article. The symmetric conductor indicated as AlFe 445/74 (Fig. 4b) with a
span length of L = 300 m is used. The maximal deflection (without icing and for ambient
temperature Tamb = 0 °C) of the power line [11] is yA = 7.73 m.
Tab. 2: The change in deflection y and maximal axial forces at the attachment point
dependant on the icing location
single power double-bundle Single power double-bundle
Icing location line power line line power line
[m]
y [m] y [m] Nmax [N] Nmax [N]
– – – 25975.7 25946.4
0 – 300 1.833 1.834 49229.9 49306.2
25 – 275 1.836 1.837 48387.4 48468.6
50 – 250 1.806 1.808 46162.9 46254.5
0 – 200 1.662 1.663 45764.6 45799,8
0 – 150 0.901 0.902 39003.5 39037.1
As we can see in Tab. 2, there are small differences between the results of axial forces and
changes in deflection in single and double-bundle power line. The differences are caused by
spacer dampers placed on the conductors.
After ice-shedding the power line swings up. The maximal jump height ysh [m] of the power
lines are different for every range of the icing on the power lines in span. The results of swing
heights in single and double-bundle power lines are shown in Tab. 3. In the case of an
asymmetrical position of icing (0 – 200 and 0 – 150 m) the point of maximal deflection is
observed (point X in Fig. 4).
After ice-shedding the power line will undergo oscillating motion. This dynamic
phenomenon is damped by naturally by the power line itself and the damping of the
environment. The constants 0.10249 [-] and = 0.0244 [-] for the Rayleigh’s damping were
used and are calculated from the damping ratio used in [14]. Time dependency of a single power
line after ice-shedding for different ranges of ice-shedding is shown in Fig. 5. Time dependency
of axial forces at the attachment point (point B) for a single power line is shown in Fig. 6.
5 Conclusion
In this paper the FEM analysis of ice-shedding from overhead power lines - AlFe 445/74
(single and double-bundle) is presented using the commercial software ANSYS. The swing up
height and axial forces in the attachment points for different ranges of ice-shedding have been
evaluated and discussed.
From the numerical simulation is obvious that maximal axial forces at the attachment (on
insulator) and the maximum jump heights of the power line after ice-shedding of partial icing
are smaller than that in ice shedding over the entire span. It is also difficult to determine partial
icing on the power line. Therefore, it is suitable to use extreme cases of ice-shedding over the
entire span to assess jump height and axial forces at the attachment points of the power line
after ice-shedding.
ACKNOWLEDGEMENT
This work was supported by the Slovak Research and Development Agency under the contract
No. APVV-0246-12 and APVV-14-0613, by Grant Agency VEGA, grant No. 1/0228/14 and
1/0453/15.
REFERENCES
[1] Chuan Wu, Bo Yana, Liang Zhang, Bo Zhang, Qing Li. A method to calculate jump height
of iced transmission lines after ice-shedding. Cold Regions Science and Technology, 2016
(125), 40 – 47
[2] VT. Morgan, DA. Swift Jump height of overhead-line conductors after the sudden release
of ice loads. Proc Inst Electr Eng 1964 (111), 1736 – 1746.
[3] M. Roshan Fekr, G. McClure. Numerical modelling of the dynamic response of ice-
shedding on electric transmission lines, Atmos. Res., 1998 (46), 1 – 11.
[4] L.E. Kollár, M. Farzaneh Vibration of bundled conductors following ice shedding. IEEE
Trans. Power Delivery, 2008, (23), No. 2, 1097 – 1104.
[5] L.E. Kollár, M. Farzaneh, P. Van Dyke. Modeling ice shedding propagation on
transmission lines with or without interphase spacers. IEEE Trans Power Apparatus Syst
2013 (28), 261 – 267.