# A Contrast on Conductor Galloping Amplitude Calculated by Three Mathematical Models with Different DOFs

Authors: Bin Liu; KuanJun Zhu; XiaoQin Sun; Bing Huo; XiJun Liu
Journal: Shock and Vibration
 (2014)

---

## Abstract

It is pivotal to find an effective mathematical model revealing the galloping mechanism. And it is important to compare the difference between the existing mathematical models on the conductor galloping. In this paper, the continuum cable model for transmission lines was proposed using the Hamilton principle. Discrete models of one DOF, two DOFs, and three DOFs were derived from the continuum model by using the Garlekin method. And the three models were compared by analyzing the galloping vertical amplitude and torsional angle with different influence factors. The influence factors include wind velocity, flow density, span length, damping ratio, and initial tension. The three-DOF model is more accurate at calculating the galloping characteristics than the other two models, but the one-DOF and two-DOF models can also present the trend of galloping amplitude change from the point view of qualitative analysis. And the change of the galloping amplitude relative to the main factors was also obtained, which is very essential to the antigalloping design applied in the actual engineering.


---

## Body

## 1. Introduction

Conductor galloping is a typical self-excited aeroelastic vibration phenomenon [1]. This conductor motion is characterized by large amplitude (possibly > 10 m) and low frequency (approximately 0.1–3 Hz) [2]. Galloping can cause various kinds of structural and electrical damages for overhead transmission lines [3, 4]. Therefore, it has become a difficult problem for transmission technology and excited more and more attention in the world.Conductor galloping represents a classical motion instability mechanism in the steady flow over a noncircular cross section caused by ice accretion on the conductor. It is pivotal to find an effective mathematical model to reveal the galloping mechanism. The galloping models have been studied extensively over a half century since den Hartog [5] firstly established the vertical galloping mechanism based on a simple single degree-of-freedom (DOF) spring-oscillator model. And the quasisteady hypothesis (QSH) was applied to describe the nonlinear aerodynamic forces, whereafter Nigol and Buchan [6, 7] proposed a galloping torsional mechanism based on a two-DOF torsion-oscillator model coupled by vertical and torsional oscillations. They proved that torsional oscillation sometimes plays a crucial role to excite the galloping phenomenon, especially for the bundled conductors. Then, Yu et al. [8, 9] obtained a three-DOFs oscillator model coupled by vertical, horizontal, and torsional oscillations. They found that horizontal oscillation has an important influence on the galloping critical conditions and galloping characteristics. But these models are all based on the simple spring-oscillation simulation for iced conductor cross section in two-dimensional space. Some important factors, such as span length, initial tension, and sag, cannot be considered in these models.In order to remedy the above disadvantages of the simple oscillator models, the cable modeling method was proposed in recent years. And the models for transmission lines have been studied in many works, which can also be classified as a continuum mathematical model. Luongo et al. [10, 11] and Lee and Perkins [12] established a flexible cable model which only accounts for the static rotation. Yu et al. [9], McConnell and Chang [13], and White et al. [14] employed a cable-beam model, accounting for twisting but not for bending and neglecting the cable initial curvature in defining the torsion strain. Because the bending moment is not taken into account in these models, the moment equilibrium around the tangent to the cable is violated [15]. So a consistent cable-beam model has become a matter of great interest, which can be able to take into account all the factors [16, 17]. Recently, Luongo et al. [18–20] have formulated consistent linear and nonlinear models of cable-beam accounting for both torsion and bending. In order to apply in the actual engineering easily, Qin et al. [21] constructed one-DOF, two-DOF, and three-DOF models from the continuum models of transmission lines by using the Hamilton principle. But the comparison of the three models is not studied, which is very important in the application of these models.In this paper, the continuum cable model for transmission lines was proposed using the Hamilton principle, where the initial displacement, the geometric nonlinearity caused by the large deformation, and the aerodynamic nonlinearity caused by flow were considered. And discrete models of one DOF with vertical motion, two DOFs with vertical-torsional coupling motion, and three DOFs with horizontal-vertical-torsional coupling motion were derived from the continuum model by using the Galerkin method. In order to compare the three models, the galloping amplitude was analyzed by the different factors including wind velocity, flow density, span length, damping ratio, and initial tension. The difference in the three models was proposed by comparing the vertical amplitude and the torsional angle caused by the conductor galloping. And the change of the galloping amplitude relative to the influence factors was also presented, which is much important to the antigalloping design applied in the actual engineering.


## 2. Construction of the Models

### 2.1. Cable Model

Figure1(a) shows the schematic diagram of the transmission line under initial tension T0. l is the span length. u, v, and w are the longitudinal, vertical, and horizontal displacements of an arbitrary point A in the line, respectively. The static balance state of the transmission line is selected as the initial position. The displacement of an arbitrary infinitesimal element dx in the line is described schematically in Figure 1(b). AB is the initial position with the static arc length ds0 and A′′B′′ is the movement position with the changed arc length ds taken from the line at time Δt. Suppose that the displacements of the start and the end in longitudinal (x), horizontal (z), and vertical (y) directions are u(x,t), w(x,t), andv(x,t) and u+uxdx, w+wxdx, and v+vxdx, respectively, where the right subscript “x” indicates differentiation with respect to coordinate x.Model of the transmission line: (a) cable model and (b) infinitesimal element.
(a)
(b)The initial arc length vector of the line infinitesimal element can be described as(1)ds0⃑=dxi⃑+dyj⃑+dzk⃑,
where i⃑, j⃑, and k⃑ are the unit vectors of x, y, and z directions.Then, arc length vector after movement can be derived as(2)ds⃑=ds0⃑+BB′⃑-AA′⃑+B′B′′⃑-A′A′′⃑=(dx+uxdx)i⃑+(dy+vxdx)j⃑+wxdxk⃑,
where k⃑ is the unit vector of z direction.Using arc length formula, the arc length expressions of the start and end in the movement can be described as(3)ds0=1+y0x2dx,ds=(1+ux)2+(y0x+vx)2+wx2dx,
where y0x is the initial configuration of the line under the initial tension. And y0x can be expressed simply using catenary shape function as [22]
(4)y0=-2σ0γshγx2σ0shγ(l-x)2σ0,
where σ0 is the axial stress at the lowest position of the line and γ is combined load in the line.Suppose thatdy0x2 can be omitted compared to 1 and ux can be omitted compared to vx and wx. Using power series expansion and omitting the terms higher than fourth order, the line element strain can be derived from the expressions (1)–(4):
(5)ε=ds-ds0ds0=y0xvx-y0x3vx+14y0x5vx-y0x2vx2+12wx2-12y0x2wx2+38y0x4vx2-18vx4+116y0x2vx4-12y0xvx3+14y0x3vx3+116y0x2wx4-12y0xvxwx2+14y0x3vxwx2+18y0x2vx2wx2-18wx4+12vx2+18y0x4wx2-14vx2wx2.


### 2.2. Torsional Motion

It is well known that the torsional displacements have a considerable effect on the instability of the conductor galloping. Especially for bundled conductors, many observations and studies verified that the torsional motion is an integral part of the galloping phenomenon and the initiation of this vibration is exactly based on the coupling between the translational motion and the torsional motion. So the torsional motion is also an essential factor to analyze the conductor galloping. Then, the torsional motion is considered as in Figure2 [21].Coordinate system: (- -) original location, (—) location after torsion.
(a)
(b)Assuming thate1 is the eccentricity in z direction and e2 is the eccentricity in y direction, then the displacements of the center of mass are
(6)v=v-+e1sinθ+e2cos⁡θ,w=w-+e1cos⁡θ-e2sinθ,
where v- and w- are the displacements of y direction and z direction, respectively.Forθ is small, (6) can be simplified as
(7)v=v-+e1θ+e2,w=w-+e1-e2θ.And the torsional strain can be defined as(8)εθ=∂θ∂x.


### 2.3. Aerodynamic Force Model

Figure3 shows a typical shape of a crescent-shaped ice accretion on the conductor surface, which is observed frequently in the case of freezing rain. In this figure, U is the absolute wind velocity perpendicular to the conductor axis in the horizontal direction. FD, FL, and M are the aerodynamic lift, drag, and moment due to wind action, respectively. θ is the torsion angle due to the moment on the conductor. θ0 is the initial static off-set angle, namely, the initial ice accretion position on the conductor. α1 is the instantaneous angle of attack with respect to the wind, which has the expression
(9)α1=tanv˙+Rθ˙cos⁡θU-ω˙+Rθ˙sinθ≈v˙+Rθ˙U,
where a dot superscript indicates differentiation with respect to time t, v˙ is the conductor velocity in the vertical direction, and θ˙ is the angular velocity of the conductor.Figure 3
The aerodynamic model of iced conductor.Then the dynamic angle of attackα can be expressed based on the quasisteady theory [23]:
(10)α=θ0+θ-v˙+Rθ˙U.This quasisteady theory assumes that the aerodynamic forces acting on the iced conductor at any instant in motion are identical to that on the stationary conductor under the same flow condition. And the use of the quasisteady theory is justified because the frequencies of conductor galloping are much lower than the vortex shedding frequencies. So the curves of the aerodynamic coefficients are obtained by the wind tunnel tests based on the stationary iced conductors at different angles of incidences. In order to establish the aerodynamic model conveniently, the polynomial curve is applied to fit onto the force coefficients versus the dynamic angle of attack by using the experimental data. Then drag coefficientCD, lift coefficient CL, and moment coefficient CM can be expressed as
(11)CL=∑i=1kc1iαi,CD=∑i=1kc2iαi,CM=∑i=1kc3iαi,
where c1i, c2i, and c3i are the fitting coefficients of the lift, drag, and moment coefficient curves, respectively, and k is the order of the polynomial. In this study, without loss of generality, k=3 is selected for the subsequent analysis. Actually, in the major existing literatures, the third-order polynomial fitting formula is most widely used.Then the aerodynamic forces have the following expressions:(12)FD=12ρU2DCD,FL=12ρU2DCL,M=12ρU2D2CM,
where ρ is the density of the flow and D is the reference length of the iced cross section that is usually substituted by the conductor diameter. And the aerodynamic forces in y and z directions are
(13)Fy=FLcos⁡α1-FDsinα1,Fz=FLsinα1+FDcos⁡α1.


### 2.4. Model with Three DOFs

Based on Sections2.1, 2.2, and 2.3, the model of conductor galloping with three DOFs (out of plan) can be constructed by using the Hamilton principle [24], which is
(14)∫t1t2(δT-δV+δW)dt=0,
where δ is the first order variational operator, T and V are the total kinetic and strain energies, respectively, and W is the work done by the nonconservative forces. And, based on the description in Sections 2.1 and 2.2, the T, V, and W are given by
(15)V=∫0l(T0+12EAε)(ds-ds0)+∫0l12GJεθ2ds,T=∫0l12mv˙2ds0+∫0l12mw˙2ds0+∫0l12Iθ˙2ds0,W=∫0l(Fy-cyv˙)vds0+∫0l(Fz-czw˙)wds0+∫0l(M-cθθ˙)θds0,
where m is the mass per unit length, E is the elastic modulus, A is the cross section area, G is the shear modulus, J is the moment of inertia, I is polar moment of inertia, and cy, cz, and cθ are the damping coefficients.Substituting Equations (1)–(13) and (15) into (14) and omitting axial inertial force, the line motion equations can be obtained. Then, the continuum motion equations can be approximated by the Galerkin procedure. In the displacement expansion, the use of the eigenfunctions is often the best choice in the aspects of convergence and accuracy [11, 19, 25]. But in order to concentrate on the study of the comparison of the different models, a simplified Galerkin procedure is presented by using a sine series of the symmetric in-plane modes as assumed shape functions [26]. So the displacements v, w, and θ can be expressed as
(16)v(x,t)=qv(t)sinπlx,w(x,t)=qw(t)sinπlx,θ(x,t)=qθ(t)sinπlx.Then, employing the Galerkin procedure, the motion equations of conductor galloping can be discretized to a three-DOF model. And neglecting the terms higher than third order, the expression is derived as(17)b101q¨v+(b102+cy)q˙v+b103qv=F1(qv,q˙v,qw,q˙w,qθ,q˙θ),b201q¨w+(b202+cz)q˙w+b203qw=F2(qv,q˙v,qw,q˙w,qθ,q˙θ),b301q¨θ+(b302+cθ)q˙θ+b303qθ=F3(qv,q˙v,qw,q˙w,qθ,q˙θ),
where F1–F3 are the functions relative to qv, q˙v, qw, q˙w, qθ, and q˙θ in which the maximum order is not greater than three and b101–b103, b201–b203, and b301–b303 are the integration constants.


### 2.5. Model with Two DOFs (y-θ Direction)

The motion in the horizontal direction is much smaller than in the vertical direction for the conductor galloping movement trajectory from the field observations in transmission lines. So, in some studies, the horizontal motion is neglected in the galloping model and only the two DOFs of vertical and torsion are considered.As a result, when the horizontal DOF is ignored, the expression (5) is changed into
(18)ε=ds-ds0ds0=y0xvx-y0x3vx+14y0x5vx+12vx2-18vx4-y0x2vx2+38y0x4vx2-12y0xvx3+14y0x3vx3+116y0x2vx4.And (17) is changed into
(19)b1q¨v+(b2+cy)q˙v+b3qv=F4(qv,q˙v,qθ,q˙θ),c1q¨θ+(c2+cθ)q˙θ+c3qθ=F5(qv,q˙v,qθ,q˙θ),
where F4 and F5 are the functions relative to qv, q˙v, qθ, and q˙θ in which the maximum order is not greater than three and b1–b3 and c1–c3 are the integration constants.


### 2.6. Model with One DOF (y Direction)

In some studies, only the vertical vibration is considered in order to discuss the galloping characteristics more conveniently for the more complex nonlinear dynamic problem after the conductor galloping has appeared. Actually, the main damages caused by the conductor galloping are attributed to the large amplitude vibration in the vertical direction. So when only the vertical DOF is considered, the expression (10) of the dynamic angle of attack is changed into
(20)α=θ0-v˙Uz.And (17) is changed into
(21)b01q¨v+(b02+cy)q˙v+b03qv=b04qv2+b05qv3+b06q˙v2+b07q˙v3,
where b01–b07 are the integration constants.


## 2.1. Cable Model

Figure1(a) shows the schematic diagram of the transmission line under initial tension T0. l is the span length. u, v, and w are the longitudinal, vertical, and horizontal displacements of an arbitrary point A in the line, respectively. The static balance state of the transmission line is selected as the initial position. The displacement of an arbitrary infinitesimal element dx in the line is described schematically in Figure 1(b). AB is the initial position with the static arc length ds0 and A′′B′′ is the movement position with the changed arc length ds taken from the line at time Δt. Suppose that the displacements of the start and the end in longitudinal (x), horizontal (z), and vertical (y) directions are u(x,t), w(x,t), andv(x,t) and u+uxdx, w+wxdx, and v+vxdx, respectively, where the right subscript “x” indicates differentiation with respect to coordinate x.Model of the transmission line: (a) cable model and (b) infinitesimal element.
(a)
(b)The initial arc length vector of the line infinitesimal element can be described as(1)ds0⃑=dxi⃑+dyj⃑+dzk⃑,
where i⃑, j⃑, and k⃑ are the unit vectors of x, y, and z directions.Then, arc length vector after movement can be derived as(2)ds⃑=ds0⃑+BB′⃑-AA′⃑+B′B′′⃑-A′A′′⃑=(dx+uxdx)i⃑+(dy+vxdx)j⃑+wxdxk⃑,
where k⃑ is the unit vector of z direction.Using arc length formula, the arc length expressions of the start and end in the movement can be described as(3)ds0=1+y0x2dx,ds=(1+ux)2+(y0x+vx)2+wx2dx,
where y0x is the initial configuration of the line under the initial tension. And y0x can be expressed simply using catenary shape function as [22]
(4)y0=-2σ0γshγx2σ0shγ(l-x)2σ0,
where σ0 is the axial stress at the lowest position of the line and γ is combined load in the line.Suppose thatdy0x2 can be omitted compared to 1 and ux can be omitted compared to vx and wx. Using power series expansion and omitting the terms higher than fourth order, the line element strain can be derived from the expressions (1)–(4):
(5)ε=ds-ds0ds0=y0xvx-y0x3vx+14y0x5vx-y0x2vx2+12wx2-12y0x2wx2+38y0x4vx2-18vx4+116y0x2vx4-12y0xvx3+14y0x3vx3+116y0x2wx4-12y0xvxwx2+14y0x3vxwx2+18y0x2vx2wx2-18wx4+12vx2+18y0x4wx2-14vx2wx2.


## 2.2. Torsional Motion

It is well known that the torsional displacements have a considerable effect on the instability of the conductor galloping. Especially for bundled conductors, many observations and studies verified that the torsional motion is an integral part of the galloping phenomenon and the initiation of this vibration is exactly based on the coupling between the translational motion and the torsional motion. So the torsional motion is also an essential factor to analyze the conductor galloping. Then, the torsional motion is considered as in Figure2 [21].Coordinate system: (- -) original location, (—) location after torsion.
(a)
(b)Assuming thate1 is the eccentricity in z direction and e2 is the eccentricity in y direction, then the displacements of the center of mass are
(6)v=v-+e1sinθ+e2cos⁡θ,w=w-+e1cos⁡θ-e2sinθ,
where v- and w- are the displacements of y direction and z direction, respectively.Forθ is small, (6) can be simplified as
(7)v=v-+e1θ+e2,w=w-+e1-e2θ.And the torsional strain can be defined as(8)εθ=∂θ∂x.


## 2.3. Aerodynamic Force Model

Figure3 shows a typical shape of a crescent-shaped ice accretion on the conductor surface, which is observed frequently in the case of freezing rain. In this figure, U is the absolute wind velocity perpendicular to the conductor axis in the horizontal direction. FD, FL, and M are the aerodynamic lift, drag, and moment due to wind action, respectively. θ is the torsion angle due to the moment on the conductor. θ0 is the initial static off-set angle, namely, the initial ice accretion position on the conductor. α1 is the instantaneous angle of attack with respect to the wind, which has the expression
(9)α1=tanv˙+Rθ˙cos⁡θU-ω˙+Rθ˙sinθ≈v˙+Rθ˙U,
where a dot superscript indicates differentiation with respect to time t, v˙ is the conductor velocity in the vertical direction, and θ˙ is the angular velocity of the conductor.Figure 3
The aerodynamic model of iced conductor.Then the dynamic angle of attackα can be expressed based on the quasisteady theory [23]:
(10)α=θ0+θ-v˙+Rθ˙U.This quasisteady theory assumes that the aerodynamic forces acting on the iced conductor at any instant in motion are identical to that on the stationary conductor under the same flow condition. And the use of the quasisteady theory is justified because the frequencies of conductor galloping are much lower than the vortex shedding frequencies. So the curves of the aerodynamic coefficients are obtained by the wind tunnel tests based on the stationary iced conductors at different angles of incidences. In order to establish the aerodynamic model conveniently, the polynomial curve is applied to fit onto the force coefficients versus the dynamic angle of attack by using the experimental data. Then drag coefficientCD, lift coefficient CL, and moment coefficient CM can be expressed as
(11)CL=∑i=1kc1iαi,CD=∑i=1kc2iαi,CM=∑i=1kc3iαi,
where c1i, c2i, and c3i are the fitting coefficients of the lift, drag, and moment coefficient curves, respectively, and k is the order of the polynomial. In this study, without loss of generality, k=3 is selected for the subsequent analysis. Actually, in the major existing literatures, the third-order polynomial fitting formula is most widely used.Then the aerodynamic forces have the following expressions:(12)FD=12ρU2DCD,FL=12ρU2DCL,M=12ρU2D2CM,
where ρ is the density of the flow and D is the reference length of the iced cross section that is usually substituted by the conductor diameter. And the aerodynamic forces in y and z directions are
(13)Fy=FLcos⁡α1-FDsinα1,Fz=FLsinα1+FDcos⁡α1.


## 2.4. Model with Three DOFs

Based on Sections2.1, 2.2, and 2.3, the model of conductor galloping with three DOFs (out of plan) can be constructed by using the Hamilton principle [24], which is
(14)∫t1t2(δT-δV+δW)dt=0,
where δ is the first order variational operator, T and V are the total kinetic and strain energies, respectively, and W is the work done by the nonconservative forces. And, based on the description in Sections 2.1 and 2.2, the T, V, and W are given by
(15)V=∫0l(T0+12EAε)(ds-ds0)+∫0l12GJεθ2ds,T=∫0l12mv˙2ds0+∫0l12mw˙2ds0+∫0l12Iθ˙2ds0,W=∫0l(Fy-cyv˙)vds0+∫0l(Fz-czw˙)wds0+∫0l(M-cθθ˙)θds0,
where m is the mass per unit length, E is the elastic modulus, A is the cross section area, G is the shear modulus, J is the moment of inertia, I is polar moment of inertia, and cy, cz, and cθ are the damping coefficients.Substituting Equations (1)–(13) and (15) into (14) and omitting axial inertial force, the line motion equations can be obtained. Then, the continuum motion equations can be approximated by the Galerkin procedure. In the displacement expansion, the use of the eigenfunctions is often the best choice in the aspects of convergence and accuracy [11, 19, 25]. But in order to concentrate on the study of the comparison of the different models, a simplified Galerkin procedure is presented by using a sine series of the symmetric in-plane modes as assumed shape functions [26]. So the displacements v, w, and θ can be expressed as
(16)v(x,t)=qv(t)sinπlx,w(x,t)=qw(t)sinπlx,θ(x,t)=qθ(t)sinπlx.Then, employing the Galerkin procedure, the motion equations of conductor galloping can be discretized to a three-DOF model. And neglecting the terms higher than third order, the expression is derived as(17)b101q¨v+(b102+cy)q˙v+b103qv=F1(qv,q˙v,qw,q˙w,qθ,q˙θ),b201q¨w+(b202+cz)q˙w+b203qw=F2(qv,q˙v,qw,q˙w,qθ,q˙θ),b301q¨θ+(b302+cθ)q˙θ+b303qθ=F3(qv,q˙v,qw,q˙w,qθ,q˙θ),
where F1–F3 are the functions relative to qv, q˙v, qw, q˙w, qθ, and q˙θ in which the maximum order is not greater than three and b101–b103, b201–b203, and b301–b303 are the integration constants.


## 2.5. Model with Two DOFs (y-θ Direction)

The motion in the horizontal direction is much smaller than in the vertical direction for the conductor galloping movement trajectory from the field observations in transmission lines. So, in some studies, the horizontal motion is neglected in the galloping model and only the two DOFs of vertical and torsion are considered.As a result, when the horizontal DOF is ignored, the expression (5) is changed into
(18)ε=ds-ds0ds0=y0xvx-y0x3vx+14y0x5vx+12vx2-18vx4-y0x2vx2+38y0x4vx2-12y0xvx3+14y0x3vx3+116y0x2vx4.And (17) is changed into
(19)b1q¨v+(b2+cy)q˙v+b3qv=F4(qv,q˙v,qθ,q˙θ),c1q¨θ+(c2+cθ)q˙θ+c3qθ=F5(qv,q˙v,qθ,q˙θ),
where F4 and F5 are the functions relative to qv, q˙v, qθ, and q˙θ in which the maximum order is not greater than three and b1–b3 and c1–c3 are the integration constants.


## 2.6. Model with One DOF (y Direction)

In some studies, only the vertical vibration is considered in order to discuss the galloping characteristics more conveniently for the more complex nonlinear dynamic problem after the conductor galloping has appeared. Actually, the main damages caused by the conductor galloping are attributed to the large amplitude vibration in the vertical direction. So when only the vertical DOF is considered, the expression (10) of the dynamic angle of attack is changed into
(20)α=θ0-v˙Uz.And (17) is changed into
(21)b01q¨v+(b02+cy)q˙v+b03qv=b04qv2+b05qv3+b06q˙v2+b07q˙v3,
where b01–b07 are the integration constants.


## 3. Contrast Analysis of Different Models on Galloping Amplitude

In this section, without loss of generality, a typical conductor is selected to analyze the galloping characteristics by using the different models under the different factors. And parameters of the line are tabulated in Table1. The fitting lift coefficients are c11=-0.1667, c12=-4.0547, and c13=8.3581. The drag coefficients are c21=-0.8605, c22=0.8325, and c23=1.7815. And the moment coefficients are c31=-0.7272, c32=0.2935, and c33=5.9704.Table 1
Parameter employed in this analysis.
Parameters
Notation
Units
Data
Wind velocity
U
m/s
9
Elastic modulus
E
GPa
47.8
Polar moment of inertia
I
Kg·m (10−3)
0.3334
Diameter of conductor
D
mm
28.6
Horizontal initial tension
T
0
kN
30
Span length
l
m
125.88
Cross-sectional area
A
mm2
423.24
Mass of unit length
m
kg/m
2.3789
Torsional stiffness
G
J
N·m2/rad
101
Initial ice angle
θ
°
40
Damping ratio
ξ
y
0.04
ξ
z
0.02
ξ
θ
0.06There are numerous methods to analyze the nonlinear dynamic equations, which can be divided into two classes, qualitative methods and quantitative methods. Qualitative methods include Multiple Scales Perturbation Method, Harmonic Balance Method, and Cell Mapping. Numerical integration method which is used to solve the galloping equations can be considered as a quantitative method, which may obtain more accurate results [27].The focus of this study is on the comparison of the difference of three models. So only a direct numerical integration method used to obtain the galloping amplitude is selected in this study. And the numerical integration is based on the fourth-order Runge-Kutta method. After specifying initial conditions, the integration is carried out until the system reaches a steady state, determined by checking the constancy of the peak-to-peak amplitudes over several consecutive periods.


### 3.1. Influence of the Wind Velocity on the Galloping Amplitude

Figure4 shows the galloping amplitudes qv and qθ versus the horizontal wind velocity U based on three different models. From Figure 4(a), it is obvious that there are the same trends of the galloping amplitude relative to the wind velocity variation based on the three models. After the galloping phenomenon starts, the galloping amplitude increases as the wind velocity increases. At the same wind velocity, the galloping amplitude calculated by the one-DOF model is much larger than the other models, and the amplitude calculated by the two-DOF model is the smallest. From the results, we can see that the coupling of the horizontal and torsional motions with the vertical motion will reduce the galloping amplitude under the condition of the same wind velocity. And the amplitude calculated by the three-DOF model is larger than the two-DOF model because of the effect of the horizontal motion.The influence of the velocity on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)From Figure4(b), torsional angles obtained by two-DOF and three-DOF models will all increase as the wind velocity increases. And the torsional angles will tend to be equal to the much higher wind velocity. So it has an obvious influence of the horizontal motion to the torsional angle in the lower wind velocity.From the amplitude responses obtained by the purely numerical procedure, it can be approximately identified as the amplitude corresponding to the leading frequencies in 1 : 2 ratio under the given wind velocity, which is consistent with some other literatures [11]. In addition, it is not difficult to obtain that the squared ratio between the transversal frequency and the torsional frequency is small, which is in accord with the torsional theory for the conductor galloping phenomenon. Figure 5 presents the amplitude-frequency spectrum of numerical solutions when U = 9 m/s.The amplitude-frequency spectrum of numerical solutions (U=9 m/s).
(a)
(b)


### 3.2. Influence of the Fluid Density on the Galloping Amplitude

Although the air density has not an obvious change under the actual conditions, the influence of the flow density on the galloping amplitude is also studied in this paper. And the results are shown in Figure6. Figure 6(a) presents the vertical amplitude versus the flow density.The influence of the flow density on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)In the given flow density interval[1,5] kg/m3, the galloping amplitudes have all increased with the increasing of the flow density based on the three models. In the same flow density, the amplitude calculated by the one-DOF model is the largest. And the amplitudes obtained by two-DOF model and three-DOF model have a little gap. So the torsional motion has an obvious influence to the galloping amplitude and the horizontal motion has not an obvious influence on the change of the flow density.Figure6(b) gives the torsional angle versus the flow density. From the results, we can see that the torsional angles decrease with the increasing of the flow density. And, at the same flow density, the torsional angle calculated by the two-DOF model is all but equal to the three-DOF model though the angle from the two-DOF model is slightly larger than the three-DOF model.It must be clarified that the actual air density is approximately 1.29 kg/m3. In this study, in order to analyze the effect of fluid density on the galloping amplitude in theory, the value of parameter ρ is given in an interval [1,5] kg/m3 as the other parameters. So the other values of the parameter ρ used in this analysis have no practical significance.


### 3.3. Influence of the Span Length on the Galloping Amplitude

Figure7 shows the galloping vertical amplitude qv and qθ versus the span length L based on three different models. In this study, the span length interval is selected as [50,260] m, in which it is familiar that the galloping modal is presented with one loop. From Figure 7(a), we can see that the galloping amplitudes have all increased with the increasing of the span length based on the three models. At the same span length, the galloping amplitude calculated by the one-DOF model is much larger than the other models, and the amplitude calculated by the two-DOF model is the smallest. Figure 7(b) presents the torsional angle versus the span length. From the results, we can see that the torsional angles obtained by the 2D model and the 3D model have the same results by and large.The influence of the span length on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)


### 3.4. Influence of the Damping Ratio on the Galloping Amplitude

Figure8 shows the galloping amplitudes qv and qθ versus the damping ratios ζy and ζθ based on three different models. Figure 8(a) gives the vertical amplitude versus the vertical damping ratio. In the given damping ration interval [0,0.08], the galloping amplitudes have all decreased rapidly with the increasing of the damping ratio based on the three models. So it has obvious influence of the damping to the galloping amplitude. And the amplitude obtained by the three-DOF model has the most gradual decrease with the increasing of the damping ratio. Figure 8(b) presents the torsional angle versus the torsional damping ratio. From the results, we can see that the torsional angles have an obvious decrease with the increasing of the damping ratio. And at the same damping ratio, the torsional angle calculated by the two-DOF model is all but equal to the three-DOF model though the angle obtained by the three-DOF model is slightly larger than the two-DOF model.The influence of the damping ratio on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)


### 3.5. Influence of the Initial Tension on the Galloping Amplitude

Figure9 shows the galloping amplitudes qv and qθ versus the initial tension T based on three different models. Figure 9(a) presents the vertical amplitude versus the initial tension. In the given tension interval [30,90] kN, the galloping amplitude will decrease slightly with the increasing of the initial tension. And at the same tension, the galloping amplitude calculated by the one-DOF model is much larger than the other models, and the amplitude calculated by the two-DOF model is the smallest.The influence of the initial tension on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)Figure9(b) gives the torsional angle versus the initial tension. From the results, we can see that it has an obvious difference between the two-DOF model and the three-DOF model. In the given tension interval, it has no obvious change of the initial tension to the torsional angle calculated by two-DOF model. But the torsional angle obtained by three-DOF model has a remarkable decrease as the tension increases. So the horizontal motion has an obvious influence on the torsional motion on the conductor galloping under the condition of different tensions.


## 3.1. Influence of the Wind Velocity on the Galloping Amplitude

Figure4 shows the galloping amplitudes qv and qθ versus the horizontal wind velocity U based on three different models. From Figure 4(a), it is obvious that there are the same trends of the galloping amplitude relative to the wind velocity variation based on the three models. After the galloping phenomenon starts, the galloping amplitude increases as the wind velocity increases. At the same wind velocity, the galloping amplitude calculated by the one-DOF model is much larger than the other models, and the amplitude calculated by the two-DOF model is the smallest. From the results, we can see that the coupling of the horizontal and torsional motions with the vertical motion will reduce the galloping amplitude under the condition of the same wind velocity. And the amplitude calculated by the three-DOF model is larger than the two-DOF model because of the effect of the horizontal motion.The influence of the velocity on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)From Figure4(b), torsional angles obtained by two-DOF and three-DOF models will all increase as the wind velocity increases. And the torsional angles will tend to be equal to the much higher wind velocity. So it has an obvious influence of the horizontal motion to the torsional angle in the lower wind velocity.From the amplitude responses obtained by the purely numerical procedure, it can be approximately identified as the amplitude corresponding to the leading frequencies in 1 : 2 ratio under the given wind velocity, which is consistent with some other literatures [11]. In addition, it is not difficult to obtain that the squared ratio between the transversal frequency and the torsional frequency is small, which is in accord with the torsional theory for the conductor galloping phenomenon. Figure 5 presents the amplitude-frequency spectrum of numerical solutions when U = 9 m/s.The amplitude-frequency spectrum of numerical solutions (U=9 m/s).
(a)
(b)


## 3.2. Influence of the Fluid Density on the Galloping Amplitude

Although the air density has not an obvious change under the actual conditions, the influence of the flow density on the galloping amplitude is also studied in this paper. And the results are shown in Figure6. Figure 6(a) presents the vertical amplitude versus the flow density.The influence of the flow density on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)In the given flow density interval[1,5] kg/m3, the galloping amplitudes have all increased with the increasing of the flow density based on the three models. In the same flow density, the amplitude calculated by the one-DOF model is the largest. And the amplitudes obtained by two-DOF model and three-DOF model have a little gap. So the torsional motion has an obvious influence to the galloping amplitude and the horizontal motion has not an obvious influence on the change of the flow density.Figure6(b) gives the torsional angle versus the flow density. From the results, we can see that the torsional angles decrease with the increasing of the flow density. And, at the same flow density, the torsional angle calculated by the two-DOF model is all but equal to the three-DOF model though the angle from the two-DOF model is slightly larger than the three-DOF model.It must be clarified that the actual air density is approximately 1.29 kg/m3. In this study, in order to analyze the effect of fluid density on the galloping amplitude in theory, the value of parameter ρ is given in an interval [1,5] kg/m3 as the other parameters. So the other values of the parameter ρ used in this analysis have no practical significance.


## 3.3. Influence of the Span Length on the Galloping Amplitude

Figure7 shows the galloping vertical amplitude qv and qθ versus the span length L based on three different models. In this study, the span length interval is selected as [50,260] m, in which it is familiar that the galloping modal is presented with one loop. From Figure 7(a), we can see that the galloping amplitudes have all increased with the increasing of the span length based on the three models. At the same span length, the galloping amplitude calculated by the one-DOF model is much larger than the other models, and the amplitude calculated by the two-DOF model is the smallest. Figure 7(b) presents the torsional angle versus the span length. From the results, we can see that the torsional angles obtained by the 2D model and the 3D model have the same results by and large.The influence of the span length on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)


## 3.4. Influence of the Damping Ratio on the Galloping Amplitude

Figure8 shows the galloping amplitudes qv and qθ versus the damping ratios ζy and ζθ based on three different models. Figure 8(a) gives the vertical amplitude versus the vertical damping ratio. In the given damping ration interval [0,0.08], the galloping amplitudes have all decreased rapidly with the increasing of the damping ratio based on the three models. So it has obvious influence of the damping to the galloping amplitude. And the amplitude obtained by the three-DOF model has the most gradual decrease with the increasing of the damping ratio. Figure 8(b) presents the torsional angle versus the torsional damping ratio. From the results, we can see that the torsional angles have an obvious decrease with the increasing of the damping ratio. And at the same damping ratio, the torsional angle calculated by the two-DOF model is all but equal to the three-DOF model though the angle obtained by the three-DOF model is slightly larger than the two-DOF model.The influence of the damping ratio on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)


## 3.5. Influence of the Initial Tension on the Galloping Amplitude

Figure9 shows the galloping amplitudes qv and qθ versus the initial tension T based on three different models. Figure 9(a) presents the vertical amplitude versus the initial tension. In the given tension interval [30,90] kN, the galloping amplitude will decrease slightly with the increasing of the initial tension. And at the same tension, the galloping amplitude calculated by the one-DOF model is much larger than the other models, and the amplitude calculated by the two-DOF model is the smallest.The influence of the initial tension on the galloping amplitude: (a) vertical amplitude and (b) torsional amplitude.
(a)
(b)Figure9(b) gives the torsional angle versus the initial tension. From the results, we can see that it has an obvious difference between the two-DOF model and the three-DOF model. In the given tension interval, it has no obvious change of the initial tension to the torsional angle calculated by two-DOF model. But the torsional angle obtained by three-DOF model has a remarkable decrease as the tension increases. So the horizontal motion has an obvious influence on the torsional motion on the conductor galloping under the condition of different tensions.


## 4. Conclusions

In this paper, one-DOF, two-DOF, and three-DOF models on the iced conductor galloping were proposed based on the continuum cable model. The contrast with the three mathematical models of different DOFs was also studied by analyzing the galloping amplitude obtained under different factors. The obtained results could be summarized as follows.(1)
These three models can all present the trend of galloping amplitude change relative to the factors from the point view of qualitative analysis.(2)
At the same value of the influence factor, the galloping vertical amplitude calculated by the one-DOF model is much larger than the other models. So the accessorial coupling of the horizontal and torsional motions with the vertical motion will reduce the galloping vertical amplitude. And the vertical amplitude obtained by the three-DOF model is larger than the two-DOF model due to the effect of the horizontal motion.(3)
The torsional motion amplitude on the conductor galloping has an obvious influence on the change of the given factors. It is more sensitive for the torsional angle obtained by the three-DOF model than the two-DOF model with the change of the initial tension. And there is no obvious difference between the two models for other influence factors.(4)
It is more conspicuous for the wind velocity and the damping ratio to influence the galloping vertical amplitude. Actually, the vertical amplitude is the most important characteristic on the conductor galloping to cause the damages of transmission lines. So it will be an effective method to control the galloping phenomenon by modifying the aerodynamic property and the damping ratio of the conductors in the antigalloping design.



---

*Source: 781304-2014-07-21.xml*