# Investigation on Influences of  Two Discrete Methods on Galloping Characteristics of Iced Quad Bundle Conductors

Authors: Xiaohui Liu; Guangyun Min; Chuan Wu; Mengqi Cai
Journal: Advances in Civil Engineering
 (2020)

---

## Abstract

The partial differential galloping equation of iced quad conductors can be transformed into an ordinary differential galloping equation by two discrete methods: one is a direct discrete method and the other is an indirect discrete method. The two discrete methods are reasonable and effective and have their own advantages and disadvantages, but whether the two different discrete methods would cause the differences in galloping characteristics of the iced quad conductor has not been studied. Based on this concept, this paper studies this problem systematically. Firstly, based on the variational principle for Hamiltonian, the partial differential galloping equation with 3DOFs of the iced quad bundle conductor is derived and then two discrete methods are used to transform the partial differential galloping equation into an ordinary differential galloping equation. One is to use a direct method to transform partial differential galloping equation into an ordinary differential galloping equation, while the other is to use an indirect method to transform partial differential galloping equation into an ordinary differential galloping equation. Secondly, based on the wind tunnel test, the three-component aerodynamic coefficients of each subconductor of the iced quad conductor are obtained, and the equivalent aerodynamic coefficients at the central axis of the quad bundle conductor are obtained by using a reasonable method. Then, the aerodynamic coefficients are fitted by Taylor rules and the aerodynamic coefficients of wind angle of attack which is 55° are used in the analysis of galloping characteristics of the iced quad conductor. Finally, based on the numerical method, the displacement response of the two discrete methods is obtained. By comparing the differences of the displacement response obtained by the two discrete methods, it is found that the two discrete methods have certain influences on the phase, frequency, and amplitude of the iced quad bundle conductor. By comparing the calculation process of these two discrete methods, it can be obtained that the calculation process of the direct discrete method is more complex and the calculation process of the indirect discrete method is simpler. By comparing the calculation results of these two discrete methods, the amplitude obtained by the indirect discrete method is bigger than that obtained by the direct discrete method, especially the amplitude in the torsional direction. The research conclusion of this paper can offer some guidance to civil and electric engineering.


---

## Body

## 1. Introduction

With the rapid development of power grid, the galloping characteristics of iced transmission lines have attracted the attention of most scholars [1–5]. In fact, the span length of transmission lines is far larger than its diameter, which belongs to the flexible cable structure. For cable structure, Irvine firstly studied the statics and linear and nonlinear dynamics of the suspension cable structure and proposed important Irvine parameters [6]. Perkins divided the external excitation of elastic suspension cable into two types, parametric excitation and periodic excitation, and studied the interaction between mode shapes and mode shapes of elastic suspension cable [7]. Rega firstly systematically expounded the mechanical modeling, response, and nonlinear phenomena of the suspension cable from the theoretical analysis, then verified the correctness of the theoretical analysis from the experiment, and found lots of valuable guidance and suggestions for civil and electric engineering [8]. Some scholars also have made many contributions to the study of suspension structure. Lilien and Havard established the finite element model of conductor with 6-DOF joints based on the space curved beam theory, then used the Newton–Raphson method to solve the static nonlinear problem, and found some very meaningful conclusions [9]. Based on the existing research, Yan et al. considered the influences of bending stiffness on the suspension cable and modified the error of bending stiffness not considered in the existing literature [10]. According to the variational principle for Hamiltonian, Xiao et al. deduced the nonlinear dynamic equation with 3DOFs considering bending stiffness of suspension cable, then analyzed the possible internal resonance modes of the cable by the multiscale method, and studied the influences of bending stiffness on internal resonance by the numerical method [11].Although the existing papers are based on the cable structure to analyze the vibration characteristics of iced transmission lines, the galloping of iced transmission lines is different from the vibration of cable. Because the galloping characteristics of transmission lines are related to ice type, ice thickness, and other factors, the galloping of transmission lines is essentially a fluid-structure coupling. vibration. Because the galloping of iced transmission lines is affected by aerodynamic loads, the top priority of studying the galloping characteristics is to obtain the three-dimensional aerodynamic coefficients of transmission lines. Zhao et al. analyzed the actual galloping data of a single conductor, two-bundle conductor, three-bundle conductor, and quad bundle conductor and then obtained some valuable results, which can provide some references for the study of the galloping characteristics of conductor [12]. Kimura et al. made the static wind tunnel test and dynamic wind tunnel test on the crescent iced conductor and then analyzed the differences between static aerodynamic characteristics and dynamic aerodynamic characteristics [13]. Lou et al. investigated the influences of iced thickness, initial iced accretion angle, and subconductor on the aerodynamic properties [14], and Cai et al. measured the aerodynamic coefficients of the iced 8-bundle conductor under different turbulence intensities flow by the wind tunnel experiment [15].In order to study the galloping characteristics of the iced quad conductor, the partial differential galloping equation of iced quad conductor is firstly required; then the partial differential galloping equation should be transformed into an ordinary differential galloping equation by the reasonable method. While according to existing literature, it can be found that some scholars are used to choose the direct discrete method to transform the partial differential galloping equation into ordinary differential galloping equation, while others are used to select the indirect discrete method to transform partial differential galloping equation into ordinary differential galloping equation. These two discrete methods both are correct and effective, but whether these two discrete methods would affect the galloping characteristics of the iced quad bundle conductor has not been studied by anyone. Based on this concept, this paper systematically studies the influences of the two different discrete methods on galloping characteristics of the iced quad bundle conductor. In this paper, the partial differential galloping equation of iced quad conductor with 3DOFs is established firstly, then two different methods are used to transform the partial differential galloping equation into an ordinary differential galloping equation. Then, through the wind tunnel test, the three-component aerodynamic coefficients of the quad bundle conductor are measured, and the equivalent aerodynamic coefficients at the central axis of the quad bundle conductor are obtained by a reasonable method. Finally, the galloping characteristics of the iced quad conductor are studied, and the influences of two different discrete methods on the displacement response of the iced quad conductor are also investigated.


## 2. Galloping Equation of the Iced Conductor

### 2.1. Mathematical Model and Galloping Equation

The mathematical model of transmission lines is established as shown in Figure1. The two ends of the conductor are constrained by fixed supports. Г1 is the equilibrium configuration under the action of gravity and Г2 is the dynamic configuration under the action of gravity and other external loads.Figure 1
The mathematical model of the conductor.In order to obtain the galloping equation of the conductor, the Cartesian coordinate system has been established and the left support of the conductor has been taken as the origin of the Cartesian coordinate system, which can be seen in Figure1.Because of the influences of iced, the cross-sectional area of the conductor would become noncircular, so its center of mass is not at the center of the circle. The variational principle for Hamiltonian belongs to the energy method, which is essentially the conservation of energy in the system. Although the conductor would rotate eccentrically due to the influences of iced, the variational principle for Hamiltonian would be also suitable for the derivation of the galloping equation of the conductor. Since it can be obtained from the variational principle for Hamiltonian,(1)∫δkv−∫δ∏+∫δw′=0,where δkv is the kinetic energy of the iced conductor, δ∏ is the potential energy of the iced conductor, and δw′ is the virtual work associated with gravity, external, and damping forces.Kinetic energy, potential energy, and virtual work of the conductor are, respectively, as follows:(2a)∫0lkvds=∫0l12mu˙i2ds+12∫0lJθ˙2ds+∫0lszu˙3θ˙ds−∫0lsyu˙2θ˙ds,(2b)∫0l∏ds=∫0lTε+12EAε2ds+12∫0lEI∂2ui∂s22ds+∫0lGIPεθ2ds,(2c)∫0lδw′ds=∑i=13pi−ciu˙i+pθ−cθθ˙,where SZ and Sy are static moments, EI is bending stiffness, EA is tensile stiffness, GIP is torsional stiffness, ε is tensile strain, J is moment of inertia, θ is torsional angle, ui is dynamic displacement, εθ is torsional strain, pi and pθ are external aerodynamic forces, ci and cθ are damping coefficients, T is static tension, and l is span length.Substituting equations (2a)–(2c) into equation (1), can obtain that(3)∂∂sHdxids+∂ui∂s+EAΘ=Jθ¨+Szu¨3−Syu¨2δiθ+Szδi3−Syδi2θ¨+EI∂4ui∂s4−GIp∂2θ∂x1δiθ+mu¨i,where(4)Θ=dxids∂ui∂s+12∂ui∂s2dxids+∂ui∂s.Only considering the vertical (u2), transverse (u3), and torsional (θ) movement of the conductor and ignoring the horizontal (u1) movement of the conductor and then substituting equation (4) into equation (3), the partial differential galloping equation of the conductor in vertical (u2), transverse (u3), and torsional direction (θ) can be obtained as follows:(5a)∂∂xEAdydx∂u2∂x+12∂u2∂x2+∂u3∂x2dydx+∂u2∂x+H∂u2∂x=mu¨2−Syθ¨+EI∂4u2∂x4+f2u˙2−p2,(5b)∂∂xEAdydx∂u2∂x+12∂u2∂x2+∂u3∂x2∂u3∂x+H∂u3∂x=mu¨3+Szθ¨+EI∂4u3∂x4+f3u˙3−p3,(5c)Jθ¨+Szu¨3−Syu¨2−GIp∂2θ∂x+fθθ˙−pθ=0.


### 2.2. Direct Discrete Method

The experts represented by Lou et al. [16, 17] generally transformed the partial differential galloping equation into an ordinary differential galloping equation by direct discrete method; that is, equation (4) is directly substituted into equation (3).The detailed steps of the direct discrete method are given below. The dynamic displacements can be expressed as(6a)u2x,t=ψ2xq2t,(6b)u3x,t=ψ3xq3t,(6c)θx,t=ψθxqθt,where ψ2, ψ3, and ψθ are mode shapes and q2, q3, and qθ are vibration functions.Substituting equations (6a)–(6c) into equations (5a)–(5c), then it can be obtained that(7a)∂∂xHψ2′q2+EAℓy′+ψ2′q2=mψ2q¨2−Syψθq¨θ+EIψ′‴q2+f2ψ2q˙2−p2,(7b)∂∂xHψ3′q3+EAℓψ3′q3=mψ3q¨3+Szψθq¨θ+EIψ3′‴q3+f3ψ3q˙3−p3,(7c)Jψθq¨θ−Szψ3q¨3+Syψ2q¨2=GIpψθ′qθ−fθψθq˙θ+pθ,where(8)ℓ=y′ψ2′q2+12ψ′22q22+ψ′33q32.Using the Galerkin method for equation (7a), then it can be obtained that(9a)b1q¨2+b2q˙2+b4q2+b5q22+b6q32+b8q23+b9q2q32+b10q¨θ=pp2.Using the Galerkin method for equation (7b), then it can be obtained that(9b)c1q¨3+c2q˙3+c3q3+c5q2q3+c6q22q3+c7q33+c8q¨θ=pp3.Using the Galerkin method for equation (7c), then it can be obtained that(9c)d1q¨θ+d2q˙θ+d3q¨3+d4q¨2+d5qθ=pp4.Equations (9a)–(9c) are the three-degree-of-freedom coupling ordinary differential galloping equation of iced conductor obtained by the direct discrete method, and the expression of the coefficients involved in equations (9a)–(9c) is shown in the appendix.


### 2.3. Indirect Discrete Method

Some experts [18, 19] thought that the axial galloping of the iced conductor is very weak, so the influences of its axial inertia force can be ignored. Based on this concept, equation (4) can be averaged along the span length; that is,(10)Θm=1l∫0ldxids∂ui∂s+12∂ui∂s2dxids+∂ui∂sdx.Substituting equation (10) into equations (5a)–(5c), then it can be obtained that(11a)∂∂xEAl∫0ldydx∂u2∂x+12∂u2∂x2+∂u3∂x2dxdydx+∂u2∂x+H∂u2∂x=mu¨2−Syθ¨+EI∂4u2∂x4+f2u˙2−p2,(11b)∂∂xEAl∫0ldydx∂u2∂x+12∂u2∂x2+∂u3∂x2dx∂u3∂x+H∂u3∂x=mu¨3+Szθ¨+EI∂4u3∂x4+f3u˙3−p3,(11c)Jθ¨+Szu¨3−Syu¨2−GIp∂2θ∂x+fθθ˙−pθ=0.Substituting equations (6a)–(6c) into equation (11a) and using the Galerkin method, then it can be obtained that(12a)b1′q¨2+b2′q˙2+b4′q2+b5′q22+b6′q32+b8′q23+b9′q2q32+b10′q¨θ=pp2′.Substituting equations (6a)–(6c) into equation (11b) and using the Galerkin method, then it can be obtained that(12b)c1′q¨3+c2′q˙3+c3′q3+c5′q2q3+c6′q22q3+c7′q33+c8′q¨θ=pp3′.Substituting equations (6a)–(6c) into equation (11c) and using the Galerkin method, then it can be obtained that(12c)d1′q¨θ+d2′q˙θ+d3′q¨3+d4′q¨2+d5′qθ=pp4′.Equations (12a)–(12c) are the three-degree-of-freedom coupling ordinary differential galloping equation of the iced conductor obtained by the indirect discrete method.


### 2.4. Aerodynamic Loads on the Iced Conductor

The mathematical model of aerodynamic load on the iced conductor is established as shown in Figure2. The iced model is crescent, FL is aerodynamic lift, FD is aerodynamic drag, U is true wind speed, Ur is relative wind speed, and D is the diameter of the conductor.Figure 2
Cross-sectional model of the iced conductor and relative flow. (a) Cross-sectional model of the iced conductor. (b) Relative flow in the quasi-steady assumption.
(a)(b)p2 is aerodynamic load in the vertical direction, p3 is aerodynamic load in the transverse direction, pθ is the aerodynamic load in the torsional direction; then,(13a)p2=FLcosαt−FDsinαt=0.5ρU2DCyα,(13b)p3=FLsinαt+FDcosαt=0.5ρU2DCzα,(13c)pθ=0.5ρU2D2CMα,where Cy(α), Cz(α), and CM(α) are aerodynamic coefficients, α is instantaneous wind angle of attack, and ρ is air density. The aerodynamic coefficients can be fitted by a cubic curve; that is,(14a)Cyα=α1α+α2α2+α3α3,(14b)Czα=β1α+β2α2+β3α3,(14c)CMα=γ1α+γ2α2+γ3α3,where αi,βi, and γi are the three-component coefficients determined by wind tunnel test, which are related to the initial wind angle of attack, torsional angle, and relative wind speed.According to equations (13a)–(13c) and equations (14a)–(14c), then using the Galerkin method, it can be obtained that(15a)pp2=pp2′t=∫0l12ρU2D−α1Uψ2xq˙2t+α2U2ψ2xq˙2t2−α3U3ψ2xq˙2t3ψ2xdx,(15b)pp3=pp3′t=∫0l12ρU2D−β1Uψ2xq˙2t+β2U2ψ2xq˙2t2−β3U3ψ2xq˙2t3ψ3xdx,(15c)pp4=pp4′t=∫0l12ρU2D2−γ1Uψ2xq˙2t+γ2U2ψ2xq˙2t2−γ3U3ψ2xq˙2t3ψ3xdx.By substituting equations (15a)–(15c) into equations (9a)–(9c), the ordinary differential equation with considering the influences of aerodynamic loads on the iced conductor under direct discrete method can be obtained as follows:(16)b1q¨2+b2+b∗q˙2+b4q2+b5q22+b6q32+b8q23+b9q2q32+b10q¨θ+b11∗q˙22+b12∗q˙23=0,c1q¨3+c2q˙3+c3q3+c5q2q3+c6q22q3+c7q33+c8q¨θ+c9∗q˙2+c10∗q˙22+c11∗q˙23=0,d1q¨θ+d2q˙θ+d3q¨3+d4q¨2+d5qθ+d6∗q˙2+d7∗q˙22+d8∗q˙23=0,where the expression of the coefficients involved in equation (16) is shown in the appendix.By substituting equations (15a)–(15c) into equations (12a)–(12c), the ordinary differential equation with considering the influences of aerodynamic loads on the iced conductor under indirect discrete method can be obtained as follows:(17)b1′q¨2+b2′+b∗q˙2+b4′q2+b5′q22+b6′q32+b8′q23+b9′q2q32+b10′q¨θ+b11∗q˙22+b12∗q˙23=0,c1′q¨3+c2′q˙3+c3′q3+c5′q2q3+c6′q22q3+c7′q33+c8′q¨θ+c9∗q˙2+c10∗q˙22+c11∗q˙23=0,d1′q¨θ+d2′q˙θ+d3′q¨3+d4′q¨2+d5′qθ+d6∗q˙2+d7∗q˙22+d8∗q˙23=0,where the expression of the coefficients involved in equation (17) is shown in the appendix.


## 2.1. Mathematical Model and Galloping Equation

The mathematical model of transmission lines is established as shown in Figure1. The two ends of the conductor are constrained by fixed supports. Г1 is the equilibrium configuration under the action of gravity and Г2 is the dynamic configuration under the action of gravity and other external loads.Figure 1
The mathematical model of the conductor.In order to obtain the galloping equation of the conductor, the Cartesian coordinate system has been established and the left support of the conductor has been taken as the origin of the Cartesian coordinate system, which can be seen in Figure1.Because of the influences of iced, the cross-sectional area of the conductor would become noncircular, so its center of mass is not at the center of the circle. The variational principle for Hamiltonian belongs to the energy method, which is essentially the conservation of energy in the system. Although the conductor would rotate eccentrically due to the influences of iced, the variational principle for Hamiltonian would be also suitable for the derivation of the galloping equation of the conductor. Since it can be obtained from the variational principle for Hamiltonian,(1)∫δkv−∫δ∏+∫δw′=0,where δkv is the kinetic energy of the iced conductor, δ∏ is the potential energy of the iced conductor, and δw′ is the virtual work associated with gravity, external, and damping forces.Kinetic energy, potential energy, and virtual work of the conductor are, respectively, as follows:(2a)∫0lkvds=∫0l12mu˙i2ds+12∫0lJθ˙2ds+∫0lszu˙3θ˙ds−∫0lsyu˙2θ˙ds,(2b)∫0l∏ds=∫0lTε+12EAε2ds+12∫0lEI∂2ui∂s22ds+∫0lGIPεθ2ds,(2c)∫0lδw′ds=∑i=13pi−ciu˙i+pθ−cθθ˙,where SZ and Sy are static moments, EI is bending stiffness, EA is tensile stiffness, GIP is torsional stiffness, ε is tensile strain, J is moment of inertia, θ is torsional angle, ui is dynamic displacement, εθ is torsional strain, pi and pθ are external aerodynamic forces, ci and cθ are damping coefficients, T is static tension, and l is span length.Substituting equations (2a)–(2c) into equation (1), can obtain that(3)∂∂sHdxids+∂ui∂s+EAΘ=Jθ¨+Szu¨3−Syu¨2δiθ+Szδi3−Syδi2θ¨+EI∂4ui∂s4−GIp∂2θ∂x1δiθ+mu¨i,where(4)Θ=dxids∂ui∂s+12∂ui∂s2dxids+∂ui∂s.Only considering the vertical (u2), transverse (u3), and torsional (θ) movement of the conductor and ignoring the horizontal (u1) movement of the conductor and then substituting equation (4) into equation (3), the partial differential galloping equation of the conductor in vertical (u2), transverse (u3), and torsional direction (θ) can be obtained as follows:(5a)∂∂xEAdydx∂u2∂x+12∂u2∂x2+∂u3∂x2dydx+∂u2∂x+H∂u2∂x=mu¨2−Syθ¨+EI∂4u2∂x4+f2u˙2−p2,(5b)∂∂xEAdydx∂u2∂x+12∂u2∂x2+∂u3∂x2∂u3∂x+H∂u3∂x=mu¨3+Szθ¨+EI∂4u3∂x4+f3u˙3−p3,(5c)Jθ¨+Szu¨3−Syu¨2−GIp∂2θ∂x+fθθ˙−pθ=0.


## 2.2. Direct Discrete Method

The experts represented by Lou et al. [16, 17] generally transformed the partial differential galloping equation into an ordinary differential galloping equation by direct discrete method; that is, equation (4) is directly substituted into equation (3).The detailed steps of the direct discrete method are given below. The dynamic displacements can be expressed as(6a)u2x,t=ψ2xq2t,(6b)u3x,t=ψ3xq3t,(6c)θx,t=ψθxqθt,where ψ2, ψ3, and ψθ are mode shapes and q2, q3, and qθ are vibration functions.Substituting equations (6a)–(6c) into equations (5a)–(5c), then it can be obtained that(7a)∂∂xHψ2′q2+EAℓy′+ψ2′q2=mψ2q¨2−Syψθq¨θ+EIψ′‴q2+f2ψ2q˙2−p2,(7b)∂∂xHψ3′q3+EAℓψ3′q3=mψ3q¨3+Szψθq¨θ+EIψ3′‴q3+f3ψ3q˙3−p3,(7c)Jψθq¨θ−Szψ3q¨3+Syψ2q¨2=GIpψθ′qθ−fθψθq˙θ+pθ,where(8)ℓ=y′ψ2′q2+12ψ′22q22+ψ′33q32.Using the Galerkin method for equation (7a), then it can be obtained that(9a)b1q¨2+b2q˙2+b4q2+b5q22+b6q32+b8q23+b9q2q32+b10q¨θ=pp2.Using the Galerkin method for equation (7b), then it can be obtained that(9b)c1q¨3+c2q˙3+c3q3+c5q2q3+c6q22q3+c7q33+c8q¨θ=pp3.Using the Galerkin method for equation (7c), then it can be obtained that(9c)d1q¨θ+d2q˙θ+d3q¨3+d4q¨2+d5qθ=pp4.Equations (9a)–(9c) are the three-degree-of-freedom coupling ordinary differential galloping equation of iced conductor obtained by the direct discrete method, and the expression of the coefficients involved in equations (9a)–(9c) is shown in the appendix.


## 2.3. Indirect Discrete Method

Some experts [18, 19] thought that the axial galloping of the iced conductor is very weak, so the influences of its axial inertia force can be ignored. Based on this concept, equation (4) can be averaged along the span length; that is,(10)Θm=1l∫0ldxids∂ui∂s+12∂ui∂s2dxids+∂ui∂sdx.Substituting equation (10) into equations (5a)–(5c), then it can be obtained that(11a)∂∂xEAl∫0ldydx∂u2∂x+12∂u2∂x2+∂u3∂x2dxdydx+∂u2∂x+H∂u2∂x=mu¨2−Syθ¨+EI∂4u2∂x4+f2u˙2−p2,(11b)∂∂xEAl∫0ldydx∂u2∂x+12∂u2∂x2+∂u3∂x2dx∂u3∂x+H∂u3∂x=mu¨3+Szθ¨+EI∂4u3∂x4+f3u˙3−p3,(11c)Jθ¨+Szu¨3−Syu¨2−GIp∂2θ∂x+fθθ˙−pθ=0.Substituting equations (6a)–(6c) into equation (11a) and using the Galerkin method, then it can be obtained that(12a)b1′q¨2+b2′q˙2+b4′q2+b5′q22+b6′q32+b8′q23+b9′q2q32+b10′q¨θ=pp2′.Substituting equations (6a)–(6c) into equation (11b) and using the Galerkin method, then it can be obtained that(12b)c1′q¨3+c2′q˙3+c3′q3+c5′q2q3+c6′q22q3+c7′q33+c8′q¨θ=pp3′.Substituting equations (6a)–(6c) into equation (11c) and using the Galerkin method, then it can be obtained that(12c)d1′q¨θ+d2′q˙θ+d3′q¨3+d4′q¨2+d5′qθ=pp4′.Equations (12a)–(12c) are the three-degree-of-freedom coupling ordinary differential galloping equation of the iced conductor obtained by the indirect discrete method.


## 2.4. Aerodynamic Loads on the Iced Conductor

The mathematical model of aerodynamic load on the iced conductor is established as shown in Figure2. The iced model is crescent, FL is aerodynamic lift, FD is aerodynamic drag, U is true wind speed, Ur is relative wind speed, and D is the diameter of the conductor.Figure 2
Cross-sectional model of the iced conductor and relative flow. (a) Cross-sectional model of the iced conductor. (b) Relative flow in the quasi-steady assumption.
(a)(b)p2 is aerodynamic load in the vertical direction, p3 is aerodynamic load in the transverse direction, pθ is the aerodynamic load in the torsional direction; then,(13a)p2=FLcosαt−FDsinαt=0.5ρU2DCyα,(13b)p3=FLsinαt+FDcosαt=0.5ρU2DCzα,(13c)pθ=0.5ρU2D2CMα,where Cy(α), Cz(α), and CM(α) are aerodynamic coefficients, α is instantaneous wind angle of attack, and ρ is air density. The aerodynamic coefficients can be fitted by a cubic curve; that is,(14a)Cyα=α1α+α2α2+α3α3,(14b)Czα=β1α+β2α2+β3α3,(14c)CMα=γ1α+γ2α2+γ3α3,where αi,βi, and γi are the three-component coefficients determined by wind tunnel test, which are related to the initial wind angle of attack, torsional angle, and relative wind speed.According to equations (13a)–(13c) and equations (14a)–(14c), then using the Galerkin method, it can be obtained that(15a)pp2=pp2′t=∫0l12ρU2D−α1Uψ2xq˙2t+α2U2ψ2xq˙2t2−α3U3ψ2xq˙2t3ψ2xdx,(15b)pp3=pp3′t=∫0l12ρU2D−β1Uψ2xq˙2t+β2U2ψ2xq˙2t2−β3U3ψ2xq˙2t3ψ3xdx,(15c)pp4=pp4′t=∫0l12ρU2D2−γ1Uψ2xq˙2t+γ2U2ψ2xq˙2t2−γ3U3ψ2xq˙2t3ψ3xdx.By substituting equations (15a)–(15c) into equations (9a)–(9c), the ordinary differential equation with considering the influences of aerodynamic loads on the iced conductor under direct discrete method can be obtained as follows:(16)b1q¨2+b2+b∗q˙2+b4q2+b5q22+b6q32+b8q23+b9q2q32+b10q¨θ+b11∗q˙22+b12∗q˙23=0,c1q¨3+c2q˙3+c3q3+c5q2q3+c6q22q3+c7q33+c8q¨θ+c9∗q˙2+c10∗q˙22+c11∗q˙23=0,d1q¨θ+d2q˙θ+d3q¨3+d4q¨2+d5qθ+d6∗q˙2+d7∗q˙22+d8∗q˙23=0,where the expression of the coefficients involved in equation (16) is shown in the appendix.By substituting equations (15a)–(15c) into equations (12a)–(12c), the ordinary differential equation with considering the influences of aerodynamic loads on the iced conductor under indirect discrete method can be obtained as follows:(17)b1′q¨2+b2′+b∗q˙2+b4′q2+b5′q22+b6′q32+b8′q23+b9′q2q32+b10′q¨θ+b11∗q˙22+b12∗q˙23=0,c1′q¨3+c2′q˙3+c3′q3+c5′q2q3+c6′q22q3+c7′q33+c8′q¨θ+c9∗q˙2+c10∗q˙22+c11∗q˙23=0,d1′q¨θ+d2′q˙θ+d3′q¨3+d4′q¨2+d5′qθ+d6∗q˙2+d7∗q˙22+d8∗q˙23=0,where the expression of the coefficients involved in equation (17) is shown in the appendix.


## 3. Wind Tunnel Tests

In order to obtain the three-component aerodynamic coefficients of the crescent iced quad bundle conductor, the wind tunnel test was carried out in China Aerodynamics Research and Development center. The quasi-static method is used to test the aerodynamic coefficients of a cross-sectional model of the iced conductor. The picture of 1.4 m × 1.4 m wind tunnel with low wind speed is shown in Figure3, the picture of the iced model is shown in Figure 4, and the picture of the force-measuring model is shown in Figure 5.Figure 3
Wind tunnel with low wind speed (1.4 m × 1.4 m).Figure 4
Iced model.Figure 5
Force-measuring model.After measuring the aerodynamic loads, the top motor of force-measuring model would drive force-measuring model to rotate 5°; then, the aerodynamic load is measured again and the range of measurement is 0°∼180°. The aerodynamic coefficients of the iced quad bundle conductor measured by the wind tunnel test include aerodynamic drag, lift, and moment coefficients; then, the dimensionless aerodynamic parameters are defined as follows:(18)CD=FD0.5ρU2LD,CL=FL0.5ρU2LD,CM=ML0.5ρU2LD2,where L is the effective length of the conductor, FD is the lift load of the conductor, FL is the drag load, and ML is the torque of the conductor, respectively.According to the wind tunnel test, the aerodynamic coefficients of crescent iced quad bundle conductor with 12 mm iced thickness under the wind speed of 10 m/s, 12 m/s, 14 m/s, and 18 m/s are shown in Figure6.Figure 6
Aerodynamic coefficients under various wind velocities. (a) Aerodynamic coefficients of subconductor 1. (b) Aerodynamic coefficients of conductor 2. (c) Aerodynamic coefficients of subconductor 3. (d) Aerodynamic coefficients of subconductor 4.
(a)(b)(c)(d)Form Figure6, when α is less than 40°, the curves of the aerodynamic lift coefficients show an upward trend. When α is in the range of 40°∼120°, the curves of the aerodynamic lift coefficients show a downward trend. When α is in the range of 120°∼160°, the curves of the aerodynamic lift coefficients have another downing range, but its downward trend is more obvious than that of α being in the range of 40°∼120°. When α is in the range of 160°∼180°, the curves of aerodynamic lift coefficients show an upward trend again.The aerodynamic drag coefficients of subconductor 1 are affected by the wake of subconductor 2 whenα is 135°. The aerodynamic drag coefficients of subconductor 3 drop suddenly influenced by the wake of subconductor 2 when α is near 45°. The aerodynamic drag coefficients of subconductor 4 are affected by the wake effect of subconductor 1, subconductor 2, and subconductor 3 when α is about 45°, 90°, and 135°, respectively.Whenα is less than 45°, the curves of the aerodynamic torque coefficients increase with α. When α is greater than 45°, the curves of the aerodynamic torque coefficients decrease with α.Because the subconductor of the iced quad bundle would be effect by wake effects, it is not reasonable to choose the aerodynamic coefficients of one subconductor to study the galloping characteristics of the iced quad bundle conductor. In practical civil engineering, the galloping of the quad bundle conductor usually is expressed as the whole vibration of the quad bundle conductor, so the aerodynamic coefficients at the central axis of the quad bundle conductor are obtained by choosing a reasonable method in this paper.In order to obtain the three-component aerodynamic coefficients at the central axis of the quad bundle, the equivalent aerodynamic lift coefficients and drag coefficients at the central axis of the quad bundle are defined as follows:(19a)CLN=1N∑i=1NCLi,(19b)CDN=1N∑i=1NCDi,where CLi and CDi, respectively, represent the aerodynamic lift coefficients and aerodynamic drag coefficients of each subconductor of the quad bundle conductor and N represents the number of each subconductor.The equivalent aerodynamic moment coefficients at the central axis of the quad bundle are defined as follows:(19c)CMN=1NCM,SN+CM,DN+CM,LN,where the physical meaning of CM,SN, CM,DN, and CM,LN is the contribution of the torque, drag, and lift of the subconductor to the quad bundle conductor.According to equations (19a)–(19c), it can be obtained that three-component aerodynamic coefficients at the central axis of the quad bundle under the wind speed of 18 m/s are shown in Figure 7.Figure 7
Aerodynamic coefficients at the central axis of the quad bundle.It can be seen from Figure7 that the curve of the equivalent aerodynamic lift coefficient of the quad bundle has a negative slope when α is in the ranges of 40°–120° and 160°–180° and it would be easy to gallop. Selecting the three-component aerodynamic coefficients when α is 55°, then according to Taylor’ rules, it can be obtained that(20a)Cyα=Cy|α=0+∂Cy∂α|α=0α+∂2Cy∂α2|α=0α2+∂3Cy∂α3|α=0α3+Οα4,(20b)Czα=Cz|α=0+∂Cz∂α|α=0α+∂2Cz∂α2|α=0α2+∂3Cz∂α3|α=0α3+Οα4,(20c)CMα=CM|α=0+∂CM∂α|α=0α+∂2CM∂α2|α=0α2+∂3CM∂α3|α=0α3+Οα4,where the right items of equation (20a) can be expressed as(21a)Cy|α=0=CLcosα+CDsinα=CL|α=0,(21b)∂Cy∂α|α=0=∂CL∂αcosα−CLsinα+∂CD∂αsinα+CDcosα|α=0=∂CL∂α+CD,(21c)∂2Cy∂α2|α=0=∂2CL∂α2cosα−2∂CL∂αsinα−CLcosα+∂2CD∂α2cosα+2∂CD∂αcosα−CDsinα|α=0=∂2CL∂α2−CL+2∂CD∂α,(21d)∂3Cy∂α3|α=0=∂3CL∂α3cosα−3∂2CL∂α2sinα−3∂CL∂αcosα+CLsinα+∂3CD∂α3sinα+3∂2CD∂α2cosα−3∂CD∂αsinα−CDcosα|α=0=∂3CL∂α3−3∂CL∂α+3∂2CD∂α2−CD.For the constant term in equations (21a)–(21d) is a fixed value, which would not affect the galloping characteristics of the iced quad bundle conductor, then by ignoring the constant term and substituting equations (21a)–(21d) into equation (20a), it can be obtained that(22a)Cy=∂CL∂α+CDα+∂2CL∂α2−CL+2∂CD∂αα2+∂3CL∂α3−3∂CL∂α+3∂2CD∂α2−CDα3+Οα4.In the same way, the aerodynamic coefficients in the transverse and torsional direction can be expressed as(22b)Cz=∂CD∂α+CLα+∂2CD∂α2−CD−2∂CL∂αα2+∂3CD∂α3−3∂CD∂α−3∂2CL∂α2+CLα3+Οα4,(22c)Czα=∂CM∂αα+∂2CM∂α2α2+∂3CM∂α3α3+Οα4.Let the coefficients of the primary term, the coefficient of the secondary term, and the coefficient of the tertiary term ofα in equation (22a) be α1, α2, and α3, respectively; that is,(23)∂CL∂α+CD=α1,∂2CL∂α2−CL+2∂CD∂α=α2,∂3CL∂α3−3∂CL∂α+3∂2CD∂α2−CD=α3.Then,Cy(α) can be expressed as(24a)Cy=α1α+α2α2+α3α3.In the same way,CZ(α) and CM(α) can be expressed as(24b)Cz=β1α+β2α2+β3α3,(24c)CM=γ1α+γ2α2+γ3α3.Finally, it can be obtained that(25a)Cy=−0.96060α−1.40716α2+97.62315α3,(25b)Cz=2.39795α−7.00826α2−119.95612α3,(25c)CM=2.6744α+22.30868α2−132.0081α3.


## 4. Numerical Analysis

The galloping characteristics of the iced quad bundle conductor are investigated and the galloping characteristics of the iced quad bundle conductor obtained by these two discrete methods will be compared in this paper. Since the physical parameters of the iced quad bundle conductor are involved, the parameters of the quad bundle conductor are listed in Table1 [20].Table 1
Physical parameters of the conductor.
ParametersNotationUnitsFor lineAxial rigidityAE106 N13.30Torsional rigidityGIPNm2rad−1101DiameterD10−3 m18.80TensionH103 N21.73Damping (y)ζy10–20.067Damping (z)ζz10–30.067Damping (θ)ζθ10–20.253Mass per unit lengthMkg m−11.53The moment of inertiaJ10−4 kg m2m−157.02The static momentSykg m m−10.000459The static momentSZkg m m−1−0.000145By substituting the physical parameters into equations (16) and (17) and combining the three-component aerodynamic coefficients in equation (25), the displacement response of the iced quad bundle conductor under direct discrete method and indirect discrete method can be obtained by using the Runge–Kutta method, as shown in Figures 8 and 9.Figure 8
Displacement response under direct discrete method. (a) Vertical displacement. (b) Transverse displacement. (c) Torsional displacement.
(a)(b)(c)Figure 9
Displacement response under indirect discrete method. (a) Vertical displacement. (b) Transverse displacement. (c) Torsional displacement.
(a)(b)(c)From Figures8 and 9, it can be seen that the iced quad bundle conductor gradually tends to be stable at 400 s under these two discrete methods; that is, the differences in the discrete method would not change the time when the iced quad bundle conductor reaches the steady state. It also can be found that the amplitudes obtained by these two discrete methods are different to some extent, and the most obvious differences are amplitude in the torsional direction.In order to compare the differences of the displacement response obtained by these two discrete methods more clearly, selecting the data as the iced quad bundle conductor is in the steady state to obtain the local curves of displacement response, which are shown in Figure10.Figure 10
Local displacement. (a) Vertical displacement, (b) Transverse displacement, (c) Torsional displacement.
(a)(b)(c)It can be seen from Figure10 that, in the transverse and torsional directions, the amplitude of the indirect discrete method is greater than that of the direct discrete method, and the differences of the amplitude in the torsional direction are particularly obvious. The frequency, phase, and amplitude of the iced quad bundle conductor obtained by the two discrete methods are different to some extent. Different discrete methods may lead to some errors between the theoretical results and the actual conditions, which is a noteworthy phenomenon.When the iced quad bundle conductor is galloping, the amplitude is one of the most concerned parameters for the studying of galloping characteristics. Because the amplitude in the torsion direction obtained by the two discrete methods has obvious differences, in order to compare the amplitude in the torsion direction obtained by the two discrete methods more clearly, amplitude diagrams below are obtained which are shown in Figure11.Figure 11
Amplitude diagrams. (a) Direct discrete method. (b) Indirect discrete method.
(a)(b)Form Figure11, the maximum amplitude in the torsion direction obtained by the direct discrete method is 0.02680 m, while the maximum amplitude in the torsion direction obtained by the indirect discrete method is 0.06300 m. This is because the indirect discrete method ignores the axial inertial force of the iced quad bundle conductor, which leads to the increase of the amplitude in the torsional direction compared with the amplitude in the torsional direction obtained by the direct discrete method.


## 5. Conclusions

Firstly, based on the variational principle for Hamiltonian, the partial differential galloping equation with 3DOFs of the iced quad bundle conductor is derived; then two discrete methods are used to transform the partial differential galloping equation into an ordinary differential galloping equation. Secondly, based on the wind tunnel test, the three-component aerodynamic coefficients of each subconductor of the iced quad conductor are obtained; then, the equivalent aerodynamic coefficients at the central axis of the quad bundle conductor are obtained by using a reasonable method. Finally, based on the numerical method, the displacement responses of the two discrete methods are obtained. The following conclusions can be listed:(1)
Based on the wind tunnel test, it can obtain that each subconductor of the iced quad bundle conductor would be effect by wake effects, so it is not reasonable to choose the aerodynamic coefficients of one subconductor to investigation on the galloping characteristics of the iced quad bundle conductor. For studying the galloping characteristics of the iced quad bundle conductor, three-component aerodynamic coefficients at the central axis of the quad bundle should be obtained firstly.(2)
Based on the numerical method, it can be obtained that the amplitude of the indirect discrete method is greater than that of the direct discrete method, and the differences of the amplitude in the torsional direction are particularly obvious. The frequency and phase of the galloping of the iced quad bundle conductor obtained by the two discrete methods are also different to some extent. This is because the indirect discrete method ignores the axial inertial force of the iced quad bundle conductor, which leads to the increase of the amplitude in the torsional direction compared with the amplitude in the torsional direction obtained by the direct discrete method.(3)
Based on the fact that the calculation process of the partial differential galloping equation is transformed into an ordinary differential galloping equation, it can be obtained that the calculation process of the direct discrete method is more complex and the calculation process of the indirect discrete method is simpler. However, based on the calculation results of the galloping characteristics, the amplitude obtained by the indirect discrete method is bigger than that obtained by the direct discrete method, especially the amplitude in the torsional direction. This is because the indirect discrete method ignores the axial inertial force of the iced quad bundle conductor since the amplitude in the torsional direction obtained by it is worthy of further discussion.



---

*Source: 8818728-2020-11-05.xml*