Shock and Vibration
Volume 2020 (2020), Article ID 8828319, 15 pages
https://doi.org/10.1155/2020/8828319

Research Article

Galloping Stability and Aerodynamic Characteristic of Iced Transmission Line Based on 3-DOF

Xiaohui Liu,¹ Ming Zou,² Chuan Wu,³ Bo Yan,⁴ and Mengqi Cai⁵

¹State Key Laboratory of Mountain Bridge and Tunnel Engineering, Chongqing Jiaotong University, Chongqing 400074, China
²College of Civil Engineering, Chongqing Jiaotong University, Chongqing 400074, China
³State Grid Henan Electric Power Research Institute, Zhengzhou 450052, China
⁴College of Aerospace Engineering, Chongqing University, Chongqing 400044, China
⁵College of Architecture and Civil Engineering, Chengdu University, Chengdu 610106, China

Correspondence should be addressed to Xiaohui Liu; cqdxlxh@126.com

Received 15 May 2020; Revised 3 August 2020; Accepted 28 August 2020; Published 8 September 2020

Academic Editor: Francesco Franco

Copyright © 2020 Xiaohui Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A new calculation method of critical wind speed based on three degrees of freedom (3-DOF) is proposed for galloping problem of iced transmission line. Based on the quasistatic theory, the aerodynamic load of iced transmission line is obtained, which considers the influence of transverse and torsional motion on the relative wind angle of attack. Finally, the equivalent galloping model of 3-DOF iced transmission line is established. At the initial angle of attack, the aerodynamic load is expanded by Taylor, and the unsymmetrical linear aerodynamic coefficient matrix is obtained. The Routh–Hurwitz criterion is used to judge the stability of iced transmission line system, and then the critical wind speed is calculated. The in-plane and out-plane frequencies corresponding to the first-order mode of the transmission line are solved by the analytical method and numerical simulation method. The results obtained by the two methods are compared and verified. The influence of dimensionless transmission line parameter λ on the in-plane and out-of-plane frequencies is discussed. The aerodynamic coefficients of the iced transmission line are measured by wind tunnel test and the aerodynamic characteristics are analyzed. According to the theoretical formula, the critical wind speed is calculated by MATLAB. The critical wind speed determined in this paper is compared with the critical wind speed determined by Den Hartog and Nigol theory. The influences of torsional vibration frequency, ice thickness, and ice shape on critical wind speed are analyzed. The research results of this paper have important theoretical significance for the stability judgment of iced transmission lines.

1. Introduction

Galloping is the self-excited vibration of transmission line under the action of wind load [1]. In winter, the transmission line is easy to be covered with ice, which will not only increase the gravity of the transmission line, but also cause the cross section of the transmission line to become asymmetrical. Icing is not only an important cause of transmission line galloping, but also a difficult problem [2]. Large vibration for a long time will cause serious accidents such as bolt loosening, broken hardware ring, transmission line breaking, tower arm damage, and tower collapse, which will cause huge loss to people’s life and property. Therefore, it is of great engineering significance to study the conditions of galloping [3].

In the field of transmission line galloping research, three major galloping theories are widely accepted internationally. The first is the single-degree-of-freedom vertical galloping theory proposed by the American scholar Hartog [4]. This theory considers the main vibration of transmission line but ignores the influence of horizontal vibration and torsional vibration on galloping. Later, based on this theory, Parkinson [5], Davenport [6], and Iwan [7] developed this theory. The second is the torsional galloping theory proposed by Canadian scholar Nigol et al. [8, 9]. According to this theory, the galloping of transmission line is caused by torsional vibration, and lateral vibration is not considered. The third is the eccentric inertial coupled galloping theory proposed by Japanese scholar Yu et al. [10]. This theory considers the eccentric coupling and ignores the influence of lateral vibration on galloping of transmission lines.

Based on the three galloping theories, experts and scholars all over the world have done a lot of research work. Liu et al. studied the characteristics of free vibration of multispan transmission [11]. McClure et al. proposed a macromodeling method for linear dynamic analysis, focusing on capturing the significant characteristics of the propagation of such impact loads in the online section [12]. Liu et al. analyzed the relationship between the galloping amplitude of the conductor and its influencing factors by using the single-degree-of-freedom model; only the single-degree-of-freedom vibration is considered [13]. John established a quasi-steady aerodynamic model for the vibration of an arbitrary cylinder with two translational degrees of freedom; the coupling vibration between the two degrees of freedom is considered [14]. In ABAQUS software, Yan et al. obtained a torsion degree of freedom cable element used to simulate the icing conductor. The two-dimensional flow around a tetrad conductor covered with different thickness of crescent ice is numerically simulated [15, 16]. The above research is mainly analyzed by numerical simulation method, which can get the results under specific circumstances, and the results are lack of universal applicability. At the same time, some necessary experiments are carried out to verify the results of numerical simulation [15].

In terms of the conductor mode and frequency, Angelo et al. established a nonlinear cable model that can be twisted. The motion is controlled by the classical equation of fully flexible cable and another equation that controls the development of torsion [17, 18]. Gurung et al. proposed a multichannel modal analysis method composed of random decrement method and feature system implementation algorithm to identify galloping [19]. Barbieri et al. tried to establish a method to identify the cable damping of transmission lines, to estimate the system damping matrix with a simple method [20]. Triantafyllou et al. deduced the in-plane modal and natural frequency formulas of the stay cable [21, 22]. Many scholars have studied the influence of gravity components from other perspectives [23, 24].

To sum up, scholars have studied galloping for a long time and obtained some results. However, because galloping is a kind of complex nonlinear vibration involving coupling between fluid and structure, there are still many problems not completely solved. As far as the authors know, most of the calculation models of the conductor galloping mechanism are based on single degree of freedom and two degrees of freedom, so it is not accurate to judge galloping. In an actual situation, galloping of transmission lines is basically a coupling vibration of 3-DOF, and energy will be transferred in the direction of 3-DOF. The galloping models of single degree of freedom and two degrees of freedom cannot meet the actual engineering requirements. Because of the coupling of galloping in horizontal, vertical, and torsional directions, galloping becomes more complex. It is necessary to establish a more accurate 3-DOF galloping model to deduce the galloping determination formula. The existing research focuses on the use of numerical simulation, and the conclusions are not universally applicable. As far as the author knows, considering the changes in wind field and conductor position at different times, the galloping mechanical model of the iced conductor with three degrees of freedom has not been studied. Based on the quasistatic theory, a 3-DOF vibration mechanical model for galloping of the iced conductor is established in this paper. The influence of the conductor structure parameters on the stability of the iced conductor is studied. The critical wind speed of the iced conductor is calculated by MATLAB software based on the stability discriminant.

2. Model Building

2.1. Galloping Mechanical Model

The high-voltage transmission lines that have galloping are mainly in open areas outside the city. In the open area, we assume that the wind speed does not change along the direction of the conductor, so this paper mainly studies the influence of stable wind on galloping. For the open area, we think the shape of ice covering along the line direction is the same. Based on the above assumption, the partial differential equation of continuum vibration can be transformed into the ordinary differential equation of cross section vibration by modal truncation method, which simplifies the complex problems. In this way, a simple and convenient formula can be obtained. The common iced shapes are crescent and sector, and the cross sections of both are symmetrical, so the displacement equation and dynamic equation of the two are the same. In this paper, the crescent ice is selected as an example to establish the mechanical model. For the 3-DOF mechanical model for analyzing galloping stability, as shown in Figure 1, it is assumed that the steady wind flows through the iced conductor from left to right, and the wind speed is U. Assuming that after the wind load is applied, the axis of symmetry of the iced conductor at rest is oZ, and the angle between the axis of symmetry and the horizontal direction is the initial angle of attack α₀, which is an important parameter for studying galloping of iced transmission line. The dynamic torsional angle of the conductor is θ, and the symmetry axis of the iced conductor after vibration is oZ′. The mass of the cross section of the equivalent iced conductor is m. The restraint stiffness in horizontal direction, vertical direction, and torsional direction are Kz, Ky, and Kθ, respectively.

Figure 1: Cross section mechanical model of the iced conductor.

In this paper, the quasistatic theory [25] is used to calculate the wind load on the iced conductor. The relative relationship between the conductor and the wind field is shown in Figure 2. The lift, drag, and torque of the iced conductor are FL, FD, and M, respectively. The lift FL is perpendicular to the relative wind speed Ur, and the drag FD is parallel to the relative wind speed Ur. The radius of the conductor is R. Under wind load, the speed will be generated when the conductor moves along the horizontal direction, where the superscript point indicates the time derivative, so the relative wind speed of the conductor in the horizontal direction is U − . When the conductor vibrates in torsion, the conductor will twist around the center of the conductor section due to the restraint of both ends of the conductor, so the torsion radius of the conductor is R. In actual situation, the torsion angle of the conductor is relatively small. The vertical motion of the conductor will produce velocity , and the torsional motion of the conductor will also produce vertical velocity component , so the relative velocity of the conductor in the vertical direction is .

Figure 2: Relative flow.

Therefore, the actual relative wind speed Ur is

The relative wind attack angle can be calculated by the following formula:

where is the dynamic torsional angle of the conductor and is the amount of change in angle due to the change in relative wind speed. The lift FL, drag FD, and torque M of the iced conductor can be obtained by the following formula:

where ρ is the air density and d are the conductor diameter. CD, CL, and CM are the drag coefficient, lift coefficient, and torque coefficient of the iced conductor, respectively. When solving the vibration equation of the iced conductor, it is necessary to project the lift and drag onto the z-axis and y-axis:

where Cz, Cy, and Cθ are aerodynamic coefficients in the horizontal, vertical, and torsional directions, respectively, which can be expressed as

At the initial wind attack angle , the horizontal, vertical, and torsional aerodynamic coefficients can be Taylor expanded, the linear part of which can be expressed as

where , , . Combining equations (4) and (6), the linear part of the aerodynamic load can be expressed as follows:

where Ca is the load damping matrix, which can be expressed as

Ka is the linear relationship between aerodynamic force and motion displacement, which can be expressed as

At the initial wind angle of attack , , combining equations (1), (2), and (5), the elements of load damping matrix and load stiffness matrix are derived and simplified. The simplified load damping matrix is

The simplified Ka is

If only the vibration of the iced conductor under static load is considered and the inertia coupling term and load stiffness term are ignored, the vibration equation of the iced conductor with 3-DOF can be obtained:

where M is the mass matrix, K is the inherent stiffness matrix, Cd is the inherent damping matrix, and X is the displacement matrix. Each term means

where m is the mass per unit length of the iced conductor and J is the moment of inertia. Ky, Kz, and Kθ are the equivalent restraint stiffness of the iced conductor in vertical, horizontal, and torsional direction, respectively. ζy, ζZ, and ζθ are the damping ratios of the iced conductor in vertical, horizontal, and torsional directions, respectively. ωy, ωz, and ωθ are the circular frequencies of the iced conductor in the vertical, horizontal, and torsional directions, respectively. According to formulas (10) and (12), the damping matrix of the system is

2.2. Stability Condition

When only vertical vibration is considered, according to the Den Hartog galloping criterion, the first term of the diagonal element in formula (14) is the criterion for judging galloping. When the following conditions are met, the system is in a stable state:

When only torsional vibration is considered, according to the Nigol galloping criterion, the third term of the diagonal element in formula (14) is the criterion for judging galloping. When the following conditions are met, the system is in a stable state:

In the actual project, the conductor will vibrate in three directions, and the above two galloping criteria cannot meet the current engineering requirements. Due to the coupling of galloping in horizontal, vertical, and torsional directions, galloping becomes more complex, so it is necessary to establish more accurate criteria for determining the galloping of three degrees of freedom. In this paper, the vibration of the iced conductor in three directions is considered simultaneously. In order to analyze the stability of the iced conductor system, let ; then the characteristic polynomial can be obtained according to formula (12):

The coefficient of β can be obtained by calculating the determinant of the matrix. Let ; then each coefficient of β can be expressed as

If the matrix is a positive definite matrix, then the 3-DOF dynamic system is stable. Use the Routh–Hurwitz criterion to judge the positive definiteness of the matrix. When the following conditions are met, the system is in a stable state:

where is the main determinant of each order of the Hurwitz coefficient matrix. The Hurwitz coefficient matrix Q can be expressed as

2.3. Galloping Modal Analysis of the Iced Conductor

The modal frequency of the conductor is the basis of studying conductor galloping, and the frequency of the conductor vibration is also needed to calculate the critical wind speed of the conductor. Considering that the frequency of the conductor is related to the physical parameters of the conductor, the frequency characteristics of the single span conductor are discussed below. According to [26], the control equations of in-plane and out-of-plane vibration of single span conductor are as follows:

where formula (21a) is the vibration control equation in plane of the conductor and formula (22b) is the vibration control equation out of plane of the conductor. H is the horizontal tension of the conductor under the self-weight. is the dynamic tension. y, z are the vertical and horizontal coordinates of the conductor under the self-weight. is the vertical and horizontal displacement of the conductor under the external load. m is the mass of the conductor per unit length. s is the arc coordinate along the length of the conductor, and t is the time. The solution of formula (21a) is complex, which can be divided into symmetrical mode and antisymmetric mode. According to the solution of the equation, the expression of in-plane mode of the conductor can be obtained as follows:

where A is the normalized modal coefficient, is the dimensionless frequency, f is the frequency, l is the conductor span, and n is the modal order.

The solution of equation (21b) is relatively simple. According to the solution of the equation, the out-of-plane mode of the conductor can be obtained as follows:

According to formulas (21a) and (22), the dimensionless frequency of in-plane symmetrical mode and antisymmetric mode can be determined as

where the comprehensive parameter λ can be expressed as

The comprehensive parameter λ is affected by the tensile rigidity EA, the initial tension H, the mass per unit length m, and the span l of the conductor. The comprehensive parameter λ reflects the physical characteristics of the conductor comprehensively and is an important parameter to study the vibration of the conductor. According to formulas (21b) and (23), it can be determined that the out-of-plane dimensionless frequency is

According to the boundary conditions, the relationship between dimensionless frequency and frequency f can be expressed as

Physical parameters of the iced conductor are given in Table 1.

Table 1: Physical parameters of the iced conductor.

Parameters
Notation
Units
Numerical value

Mass per unit length
m
kg/m³
1.53

Moment of inertia
Jθ
kg m
0.0057

Diameter of bare cable
d
m
27.6 × 10⁻³

Span length
l
m
200

Damping ratio in the y, z direction
ζy, ζz

0.005

Damping ratio in the θ direction
ζθ

0.038

Axial rigidity
EA
N
13.3 × 10⁶

Horizontal tension
H
N
41.65 × 10³

Substituting the physical parameters of the iced conductor into formulas (24)–(26), we can get the change characteristics of dimensionless frequency with the comprehensive parameter λ, as shown in Figure 3.

Figure 3: Dimensionless natural frequencies.

It can be seen from the figure that the dimensionless frequency corresponding to the first out-of-plane symmetric mode φ₁ is 1, and its value does not change with λ. The dimensionless frequency corresponding to the first in-plane symmetric mode ψ₁ increases with the increase of λ, and its value range is 1∼2.86. The dimensionless frequencies corresponding to the second in-plane antisymmetric mode ψ₂ and out-of-plane antisymmetric mode φ₂ are the same, their values do not change with λ, and their values are 2. The dimensionless frequency corresponding to the third out-of-plane symmetric mode φ₃ is 3, and its value does not change with λ. The dimensionless frequency corresponding to the third in-plane symmetric mode ψ₃ increases with the increase of λ, and its value range is 3∼4.92. From formula (24), it can be seen that the comprehensive parameter λ is positively related to the mass per unit length of the conductor m, the span l of the conductor, and the tensile rigidity EA, and the comprehensive parameter λ is negatively related to the initial tension H. From the overall view, the dimensionless frequencies of the first in-plane symmetric mode ψ₁ and the third in-plane symmetric mode ψ₃ decrease with the increase of the initial tension H. The initial tension H has a great influence on the dimensionless frequency. According to the observation and analysis of the actual conductor galloping event, it is found that the mode of galloping is mainly low order. Therefore, this paper selects the in-plane mode as ψ₁ and the out-of-plane mode as φ₁ when the gallop occurs. In order to verify the correctness of the analytical solution of modal frequency, the finite element software ABAQUS is used for numerical simulation. The finite element model of the iced conductor is established in ABAQUS, and the fixed constraints are adopted at both ends of the conductor. Considering that the span of the conductor is far larger than the diameter of the conductor, the influence of bending stiffness is ignored and three-node parabolic truss element is used to simulate the conductor. The conductor model is divided into 200 grids. The equivalent node load of truss element is composed of aerodynamic force and gravity on the conductor. Since galloping of the conductor is a typical nonlinear problem, the stiffness matrix, element unbalance force vector, and nonlinear aerodynamic load vector should be updated; according to the current displacement, after the incremental solution of each load step is completed. In order to reduce the calculation as much as possible and improve the efficiency of galloping analysis, the degree of freedom condensation method is used in this study. After the load is applied, the conductor frequency obtained by numerical simulation is shown in Figure 4.

Figure 4: Frequency corresponding to first-order in-plane and out-of-plane modes.

The dotted curve in Figure 4 represents the frequency corresponding to the first in-plane symmetric mode ψ₁. Its frequency decreases with the increase of λ, and its value is between 0.2 and 0.4 Hz. The square point curve represents the frequency corresponding to the first out-of-plane symmetrical mode φ₁. Its frequency increases first and then decreases with the increase of λ, and its value is between 0.4 and 0.6 Hz. In the entire interval, the in-plane frequency is greater than the out-of-plane frequency. The triangle point curve represents the ratio of the first in-plane frequency to the out-of-plane frequency. In the interval, the frequency ratio rises, and its value is between 1.12 and 2.99. The first in-plane and out-of-plane frequency ratio triangle curve in Figure 4 coincides with the first in-plane and out-of-plane ratio frequency triangle curve in Figure 3, indicating that the analytical solution of frequency coincides with the numerical solution, which proves the correctness and authenticity of the frequency solution.

3. Wind Tunnel Test and Aerodynamic Characteristics Analysis of the Iced Conductor

The previous chapter completes the derivation of the calculation formula of the critical wind speed of the iced conductor and the analysis of the conductor mode. The calculation of the critical wind speed requires the aerodynamic coefficient of the iced conductor. The aerodynamic coefficient provides raw data support for the calculation of the critical wind speed. The analysis of aerodynamic characteristics can explain the cause of the conductor galloping.

3.1. Wind Tunnel Equipment

This test was conducted in a 1.4 m × 1.4 m low-speed wind tunnel at the China Aerodynamics Research and Development Center Low Wind Speed Research Institute. The wind tunnel is shown in Figure 5. The wind tunnel is a direct-flow low wind speed wind tunnel. The cross-sectional area of the test section is 1.8632 m², the cross-sectional shape is a chamfered rectangle, the length of the test section is 2.8 m, and the test wind speed range is 0～65 m/s.

Figure 5: Wind tunnel with low wind speed (1.4 m × 1.4 m).

The PXI system is selected as the data collection system for the test of the iced conductor. Angle control, speed, and pressure control are realized by corresponding engineering control system. The instructions are transmitted by network communication between devices.

3.2. Iced Conductor Model

This test simulates a transmission line located in the Da Feng Ding area of Xi Chang City, Sichuan Province, China. The type of transmission line in this area is XLGJ-400/50, the safety of transmission line factor is 2.5, the maximum operating stress of transmission line is 105.9 MPa, and the average operating stress of transmission line is 66.2 MPa. The diameter of transmission line is 27.6 mm. The material of the conductor test model is duralumin, with a diameter of 27.6 mm and a length of 700 mm. The outer surface is coiled with a rubber tube to simulate the outer surface of the actual transmission line. An aerodynamic balance is installed inside the conductor model. The conductor model is shown in Figure 6.

Figure 6: Conductor model.

The material of the ice model is light wood with a length of 700 mm. During the test, the ice model is pasted on the surface of the conductor model. The iced cross-sectional shape is symmetrical crescent shape and sector shape. The sector shape outer arc angle is 120° and the inner arc angle is 140°. The wind angle of attack is defined as the angle between the counterflow direction and the symmetry axis of the iced section, facing the incoming flow direction, and turning counterclockwise is positive. The initial attack angle is 0°. Each time the aerodynamic coefficient is measured, the conductor model rotates 5°, and the rotation angle range is 0°∼180°. The ice size and attack angle are shown in Figure 7. The influence of torsional vibration on aerodynamic force during galloping is not considered, the use of quasistatic assumptions has little effect on the results. However, the effect of this assumption on the aerodynamic stability of the conductor can be neglected.

(a)

(b)

Figure 7: Iced size and attack angle. (a) Sector. (b) Crescent.

The thickness of crescent ice is 12 mm, 20 mm, and 28 mm, respectively, and that of sector ice is 18 mm and 25 mm, respectively. The ice model is shown in Figure 8.

Figure 8: Iced model.

3.3. Test Device of the Conductor

In the middle of the iced conductor test device is a conductor model. To ensure the two-dimensional characteristics of the fluid, circular end plates are installed at both ends of the conductor model as shown in Figure 9. A gap of about 5 mm is left between the lower end of the conductor model and the lower end plate to ensure no contact between the model and the lower end plate during the test. The supporting rod of the aerodynamic balance directly passes through the upper end plate and is fixedly connected with the fixed rack of the supporting device. During the test, the attack angle is changed by turning upper and lower turntable synchronously.

Figure 9: Test device of the conductor.

3.4. Test Content and Data Processing

The content of this test is to test the aerodynamic coefficient of the iced conductor model under different attack angles, different wind speeds, different ice shapes, and different ice thickness. The aerodynamic coefficients include drag coefficient, lift coefficient, and torque coefficient. The dimensionless aerodynamic coefficient is defined as follows:

where d is the diameter of the conductor model, d = 26.7 mm, U is the test wind speed, L is the length of the conductor, L = 700 mm, ρ is the air density under the test environment, and ρ = 1.29 kg/m³. The coordinate system is the wind axis coordinate system. Due to the small windward area perpendicular to the airflow, the wind tunnel blockage is very small, and the data will not need to be corrected.

3.5. Test Results and Aerodynamic Characteristics Analysis

3.5.1. Influence of Wind Speed on the Aerodynamic Coefficient

Test the aerodynamic coefficients of the crescent iced conductor with an ice thickness of 12 mm at four wind speeds of 10, 12, 14, and 18 m/s. The test results are shown in Figure 10.

(a)

(b)

(c)

Figure 10: Aerodynamic coefficient of the iced conductor at different wind speeds. (a) Lift coefficient. (b) Drag coefficient. (c) Torque coefficient.

It can be seen from Figure 10 that, under different wind speeds, the curve of aerodynamic coefficient changing with the angle of attack is approximately the same, and the wind speed has little influence on the aerodynamic coefficient of the iced conductor. It can be seen from Figure 10(a) that the curve of the lift coefficient with the angle of attack is approximately sinusoidal. When the angle of attack α = 0°, α = 90°, and α = 180°, the lift coefficient is almost zero. At the same angle of attack, with the increase of wind speed, the change of lift coefficient is very small, and the change is not regular. It can be seen from Figure 10(b) that the curve of the drag coefficient of the conductor with the angle of attack is a half-wave shape. When the angle of attack α = 0°, the windward area of the conductor is the smallest and the drag coefficient is the smallest. When the angle of attack α = 90°, the windward side of the conductor is the largest, and the drag coefficient has a peak value. At the same angle of attack, with the increase of wind speed, the drag coefficient decreases, but the change value is very small. It can be seen from Figure 10(c) that the curve of the conductor torque coefficient changing with the angle of attack is approximately sinusoidal. When the angle of attack α = 0°, α = 150°, and α = 180°, the torque coefficient is almost 0. At the same angle of attack, with the increase of wind speed, the torque coefficient is almost unchanged.

3.5.2. Influence of Ice Thickness on the Aerodynamic Coefficient

At a wind speed of 14 m/s, the aerodynamic coefficients of the crescent iced conductor at three thicknesses of 12, 20, and 28 mm were tested. The test results are shown in Figure 11.

(a)

(b)

(c)

Figure 11: Aerodynamic coefficient of the iced conductor at different ice thicknesses. (a) Lift coefficient. (b) Drag coefficient. (c) Torque coefficient.

It can be seen from Figure 11 that under different icing thickness, the curve of aerodynamic coefficient changing with the angle of attack is approximately the same, but the difference of aerodynamic coefficient is very large, and the icing thickness has a great influence on the aerodynamic coefficient of the conductor. At the same angle of attack, with the increase of ice thickness, the absolute value of aerodynamic coefficient increases.

3.5.3. Influence of Ice Shape on Aerodynamic Coefficient

When the wind speed is 14 m/s, test the aerodynamic coefficient of the iced conductor under the two ice shapes of crescent and sector. In order to make the cross-sectional area of the two ice conductors approximately the same, the ice thickness is 20 mm for crescent and 18 mm for sector. The test results are shown in Figure 12.

(a)

(b)

(c)

Figure 12: Aerodynamic coefficient of the iced conductor at different ice shapes. (a) Lift coefficient. (b) Drag coefficient. (c) Torque coefficient.

It can be seen from Figure 12 that, for different ice shapes, the curves of aerodynamic coefficients changing with angle of attack are different, and the ice shape has a great influence on the aerodynamic coefficients of the conductor. Different ice shapes will lead to different windward areas of the conductor, which will lead to different aerodynamic coefficients. At the same attack angle, the ice shape is different, the aerodynamic coefficients are very different, and the sign may also be different. Sometimes, the aerodynamic coefficients of the crescent ice conductor are larger than those of the sector ice conductor; sometimes, the aerodynamic coefficients of the sector ice conductor are larger than those of the crescent ice conductor; there is no uniform law. For the lift coefficient, the maximum value of the lift coefficient of the crescent ice conductor exactly corresponds to the minimum value of the lift coefficient of the sector ice conductor. For the torque coefficient, the curve of the torque coefficient of the crescent iced conductor and the curve of the torque coefficient of the sector iced conductor have the same trend changing with the attack angle.

4. Calculation and Analysis of Critical Wind Speed

Use MATLAB to write judgment and loop statement programs and substitute the aerodynamic coefficients obtained in chapter 3 and the conductor frequency obtained in chapter 2 into the stability discriminant formulas (15), (16), and (19) to calculate different galloping judgment criteria under various engineering situations critical wind speed. This chapter compares and analyzes the critical wind speeds calculated by different galloping criteria. Set different variables and analyze the influence of the comprehensive conductor parameter, torsion frequency, ice thickness, and ice shape on the critical wind speed. Regarding the torsion frequency of the iced conductor, there is no generally accepted research result. In order to make the result reasonably representative, this paper selects three torsion frequency values for calculation, namely, 0.2 Hz, 0.4 Hz, and 0.6 Hz. The results of the calculation are now analyzed:

4.1. Comparative Analysis of Critical Wind Speeds with Den Hartog and Nigol

According to the calculation results, the critical wind speeds of three different galloping theories at a crescent iced coating thickness of 20 mm, a torsional frequency of 0.4 Hz, and an attack angle of 165° are selected to compare and verify the 3-DOF galloping discrimination proposed in this paper. The result is shown in Figure 13.

Figure 13: Critical wind speeds at different theories.

The black curve represents the critical wind speed determined by Den Hartog theory, the red curve represents the critical wind speed determined by Nigol theory, and the green curve represents the critical wind speed determined by the 3-DOF theory in this paper. It can be seen from Figure 13 that critical wind speed determined by Den Hartog theory increases first and then decreases as the frequency ratio increases, and its value is between 6 and 7.5 m/s. Critical wind speed determined by Nigol theory does not change with the frequency ratio, and its value is 6.2 m/s. It can be seen from formula (16) that Nigol theory has nothing to do with in-plane and out-of-plane frequencies, so its critical wind speed will not change with the ratio of in-plane and out-of-plane frequencies and is a fixed value. Critical wind speed determined by the 3-DOF theory in this paper increases first and then decreases as the frequency ratio increases, and its value is between 4.6 and 6.2 m/s. Overall, the critical wind speed determined by Den Hartog theory is the largest, and the critical wind speed determined by the 3-DOF theory in this paper is the smallest. When the frequency ratio is 1.8, the critical wind speed is the largest. Since the 3-DOF theory in this paper is related to frequencies in three directions, considering the transfer of energy between three degrees of freedom, the critical wind speed is lower, and galloping is more likely to occur.

4.2. Effect of Torsion Frequency on Critical Wind Speed

In order to analyze the influence of torsion frequency on the critical wind speed of the iced conductor, according to the calculation results, the crescent-shaped ice coating thickness is 20 mm, the torsion frequency is 0.2 Hz, 0.4 Hz, and 0.6 Hz, and the angle of attack is 105°. The 3-DOF stability discriminant formula in this paper gives the critical wind speed as shown in Figure 14.

Figure 14: Critical wind speeds at different torsion frequencies.

It can be seen from Figure 14 that the torsional frequency has a greater influence on the critical wind speed. When the torsional frequency is 0.2 Hz and the torsional frequency is 0.6 Hz, the maximum difference in critical wind speed is 2.4 m/s. From the analysis of the three curves, as the torsional frequency decreases, the critical wind speed decreases. When the torsional frequency decreases, the torsional stiffness will also decrease, and conductor galloping is more likely to occur. As a whole, the critical wind speed increases first and then decrease as the frequency ratio increases. When the frequency ratio is 1.8, the critical wind speed reaches its peak value. When the frequency ratio is close to 1 and 3, the critical wind speed is low, and the conductor is easy to gallop.

4.3. Effect of Ice Thickness on Critical Wind Speed

In order to analyze the influence of ice thickness on the critical wind speed, according to the calculation results, when the torsional frequency is 0.4 Hz and the attack angle is 90°, the icing shape of the conductor is crescent, under different ice thicknesses. The critical wind speed determined by the theory of 3-DOF in this paper is used. The results are shown in Figure 15.

Figure 15: Critical wind speeds at different ice thicknesses.

It can be seen from Figure 15 that the thickness of ice has a greater influence on the critical wind speed. When the thickness of ice is 12 mm and 20 mm, the maximum critical wind speed difference reaches 4 m/s. As the thickness of the ice increases, the critical wind speed decreases. As the thickness of the ice increases, the windward surface of the conductor increases, and the aerodynamic coefficients in all three directions increase, making the conductor more prone to galloping.

4.4. Effect of Ice Shape on Critical Wind Speed

When analyzing the influence of ice shape on critical wind speed, in order to make the cross-sectional area of two types of the iced conductors approximately the same, the crescent iced conductor with 20 mm ice thickness and the sector iced conductor with 18 mm ice thickness are selected for comparative analysis, excluding the influence of cross-sectional area on critical wind speed. According to the calculation results, the critical wind speed is determined by the 3-DOF theory in this paper when the torsional frequency is 0.4 Hz and the attack angle is 85°. The results are shown in Figure 16.

Figure 16: Critical wind speeds at different ice shapes.

As can be seen from Figure 16, the critical wind speeds of the two types of the iced conductor have different trends with the ratio of in-plane and out-plane frequency. As the frequency ratio increases, the critical wind speed of the sector iced conductor increases first, then tends to remain unchanged, and finally increases. As the frequency ratio increases, the critical wind speed of the crescent iced conductor increases first and then decreases. In the entire frequency range, sometimes, the critical wind speed of the crescent ice conductor is greater than the critical wind speed of the sector iced conductor, and, sometimes, the critical wind speed of the crescent iced conductor is less than the critical wind speed of the sector iced conductor. There is no uniform rule for the effect of ice shape on the critical wind speed of the conductor galloping, and different conclusions will be drawn for different situations.

5. Conclusion

In this paper, the new mechanical model of the 3-DOF iced conductor is established, and the aerodynamic coefficient is expanded by Taylor at the initial angle of attack. The galloping vibration of the iced conductor under static wind is studied, and the vibration equation of the 3-DOF iced conductor is obtained by neglecting the inertia coupling term and the load stiffness term. According to the Routh–Hurwitz stability criterion, the equilibrium stability conditions of the transmission line system are derived. The equation of in-plane and out-of-plane vibration is established, and the influence of physical parameters of the conductor on modal frequency is studied. In the wind tunnel test, different wind speeds, ice thicknesses, and ice shapes are set. The aerodynamic coefficients in three directions of the iced conductor are measured by wind tunnel test. Using Den Hartog theory, Nigol theory, and 3-DOF theory, the critical wind speed is calculated and compared. The critical wind speeds of different torsional frequencies, ice thicknesses, and ice shapes are calculated by using the 3-DOF stability discriminant. The characteristics of the critical wind speed with the frequency ratio in and out of plane were investigated. The following conclusions were obtained:

(1)The in-plane and out-of-plane frequencies corresponding to the first-order mode of the conductor obtained by the analytical method and the numerical simulation method are similar, which verifies the accuracy of the results. For the iced single conductor, in the case of other physical parameters of transmission line being the same, the change in conductor tension will significantly affect the in-plane and out-of-plane vibration frequencies, and the frequency change will affect the critical wind speed.(2)The 3-DOF galloping judgment criterion obtained by the method in this paper considers the coupling of vertical motion, lateral motion, and torsional motion in more detail, which is more in accordance with the actual engineering situation. According to the calculation results, based on Den Hartog theory, the attack angles that may occur galloping are 15°, 165°, and 170°. Based on Nigol theory, the attack angles that may cause galloping are 10°, 165°, 170°. Based on the 3-DOF stability theory in this paper, the attack angles that may occur galloping are 10°, 55°, 60°, 70°, 80°, 95°, 165°, and 170°; the range of initial angle of attack where galloping occurs is wider than that determined by Den Hartog theory and Nigol theory. The critical wind speed determined by the theory of this paper is lower than the critical wind speed determined by the other two theories. Therefore, it is safer and more reliable to use 3-DOF stability theory in this paper to guide the antigalloping of the transmission line.(3)The frequency ratio of in-plane and out-of-plane has a significant effect on the critical wind speed of galloping. With the increase of frequency ratio, the critical wind speed first increases and then decreases. When the frequency ratio is close to 1 and 3, the critical wind speed is low, and the transmission line is prone to galloping. When the frequency ratio is close to 1.8, the critical wind speed is higher, and the transmission line is not easy to gallop. In the design of transmission line antigalloping, the frequency ratio should be close to 1.8.(4)The torsional frequency has a significant effect on the critical wind speed of galloping of iced transmission line. With the decrease of torsional frequency, the critical wind speed decreases; with the decrease of torsional frequency, the torsional stiffness decreases, and the galloping of the conductor is more likely to happen. Therefore, the transmission line parameters with higher torsional frequency should be selected for antigalloping design.(5)With the increase of ice thickness, the critical wind speed decreases. Under the snowy weather in winter, where the transmission line is easy to be covered with ice, special attention should be paid to the effect of ice on galloping in the design of transmission line antigalloping.(6)For transmission line with different ice coating shapes, the angle of attack affects the critical wind speed of the transmission line.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interests.

Acknowledgments

This study was financially supported by the National Inherent Science Foundation of China and Science (Grants nos. 51308570 and 51808085); Technology Research Project of Chongqing Education Commission (Grant no. KJ201600712182); and Basics and Cutting Edge Project of Chongqing Science and Technology Commission (Grant no. cstc2017jcyjAX0246).

References

F. Wang, K. Du, J. Sun, F. Huang, and Z. Xiong, “Shaking table array tests of an ultra-high-voltage cup-type transmission tower-line system,” Shock and Vibration, vol. 2019, Article ID 2350675, 20 pages, 2019.

T. Ishihara and S. Oka, “A numerical study of the aerodynamic characteristics of ice-accreted transmission lines,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 177, pp. 60–68, 2018.

C. B. Gurung, H. Yamaguchi, T. Yukino et al., “Identification of large amplitude wind-induced vibration of ice-accreted transmission lines based on field observed data,” Engineering Structures, vol. 24, no. 2, pp. 179–188, 2002.

J. P. D. Hartog, “Transmission line vibration due to sleet,” Transactions of the American Institute of Electrical Engineers, vol. 51, no. 4, pp. 1074–1076, 1932.

G. V. Parkinson and N. P. H. Brooks, “On the aeroelastic instability of bluff cylinders,” Journal of Applied Mechanics, vol. 28, no. 2, pp. 252–258, 1961.

M. Novak and A. G. Davenport, “Aeroelastic instability of prisms in turbulent flow,” Journal of the Engineering Mechanics Division, vol. 96, no. 1, pp. 17–39, 1970.

W. D. Iwan, “Galloping oscillations of hysteretic structures,” Journal of the Engineering Mechanics Division, vol. 99, no. 6, pp. 1129–1146, 1973.

O. Nigol and P. Buchan, “Conductor galloping part I-Den Hartog mechanism,” IEEE Transactions on Power Apparatus and Systems, vol. PAS-100, no. 2, pp. 699–707, 1981.

O. Nigol and P. Buchan, “Conductor galloping-part II torsional mechanism,” IEEE Transactions on Power Apparatus and Systems, vol. PAS-100, no. 2, pp. 708–720, 1981.

P. Yu, A. H. Shah, and N. Popplewell, “Inertially coupled galloping of iced conductors,” Journal of Applied Mechanics, vol. 59, no. 1, pp. 140–145, 1992.

X. Liu, L. Liu, M. Cai, and B. Yan, “Free vibration of transmission lines with multiple insulator strings using refined models,” Applied Mathematical Modelling, vol. 67, pp. 252–282, 2019.

G. McClure and M. Lapointe, “Modeling the structural dynamic response of overhead transmission lines,” Computers & Structures, vol. 81, no. 8–11, pp. 825–834, 2003.

Y. Liu, Z. Fu, X. H. Yang et al., “Galloping amplitude analysis and observation on full-scale test overhead transmission lines,” in Proceedings of the 2015 IEEE Power and Energy Society General Meeting, pp. 1–5, Denver, CO, USA, July 2015.

J. H. G. Macdonald and G. L. Larose, “Two-degree-of-freedom inclined cable galloping-part 1: general formulation and solution for perfectly tuned system,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 96, no. 3, pp. 291–307, 2008.

J. Hu, B. Yan, S. Zhou, and H. Zhang, “Numerical investigation on galloping of iced quad bundle conductors,” IEEE Transactions on Power Delivery, vol. 27, no. 2, pp. 784–792, 2012.

M. Cai, B. Yan, X. Lu, and L. Zhou, “Numerical simulation of aerodynamic coefficients of iced-quad bundle conductors,” IEEE Transactions on Power Delivery, vol. 30, no. 4, pp. 1669–1676, 2015.

A. Luongo, D. Zulli, G. Piccardo et al., “On the effect of twist angle on nonlinear galloping of suspended cables,” Computers & Structures, vol. 87, no. 15-16, pp. 1003–1014, 2009.

A. Luongo, D. Zulli, G. Piccardo et al., “Analytical and numerical approaches to nonlinear galloping of internally resonant suspended cables,” Journal of Sound and Vibration, vol. 315, no. 3, pp. 375–393, 2008.

C. B. Gurung, H. Yamaguchi, T. Yukino et al., “Identification and characterization of galloping of Tsuruga test line based on multi-channel modal analysis of field data,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 91, no. 7, pp. 903–924, 2003.

N. Barbieri, O. H. de Souza Júnior, R. Barbieri et al., “Dynamical analysis of transmission line cables: part 2-damping estimation,” Mechanical Systems and Signal Processing, vol. 18, no. 3, pp. 671–681, 2004.

M. S. Triantafyllou, “The dynamics of taut inclined cables,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 37, no. 3, pp. 421–440, 1984.

M. S. Triantafyllou and L. Grinfogel, “Natural frequencies and modes of inclined cables,” Journal of Structural Engineering, vol. 112, no. 1, pp. 139–148, 1986.

Q. Wu, K. Takahashi, and S. Nakamura, “Formulae for frequencies and modes of in-plane vibrations of small-sag inclined cables,” Journal of Sound and Vibration, vol. 279, no. 3–5, pp. 1155–1169, 2005.

M. Zhang, G. Zhao, J. Li et al., “Nonlinear dynamic analysis of high-voltage overhead transmission lines,” Shock and Vibration, vol. 2018, Article ID 1247523, 35 pages, 2018.

Y. M. Desai, P. Yu, N. Popplewell, and A. H. Shah, “Finite element modelling of transmission line galloping,” Computers & Structures, vol. 57, no. 3, pp. 407–420, 1995.

G. Rega and R. Alaggio, “Experimental unfolding of the nonlinear dynamics of a cable-mass suspended system around a divergence-Hopf bifurcation,” Journal of Sound and Vibration, vol. 322, no. 3, pp. 581–611, 2009.