 Open Access
 Authors : Pham Viet Hung , Ngo Quy Tuan
 Paper ID : IJERTV9IS061000
 Volume & Issue : Volume 09, Issue 06 (June 2020)
 Published (First Online): 06072020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Prediction of Wind Induced Vibration of High Voltage Electrical Bundled Conductors based on Time Domain Analysis
Pham Viet Hung
Faculty of Engineering and Technology Hue University, Hue city,
Vietnam
Ngo Quy Tuan
Faculty of Engineering and Technology Hue University, Hue city,
Vietnam
Abstract Study is made to predict wind induced vibration of high voltage electrical bundled conductors in time domain. A finite element model of bundled conductors associated model of nodal wind forces is analyzed by using Direct transient response analysis solver. The comparison of fieldmeasured data and simulation of wind induced vibration is presented to investigate the response developed in time. The result may be used in designing and controlling damages due to wind induced vibration of bundled conductors. A significant feature of this approach is that the time dependent characteristics of unsteady aerodynamic forces and the nonlinearities of both aerodynamics and structural behavior of bundled conductors can be taken into account.
KeywordsElectrical conductor; prediction, windinduced vibration, time domain

INTRODUCTION
High voltage electrical bundled conductors are used to transmit electricity from power plants over many types of terrain. The electrical conductors are long, slender, flexible, and wind sensitive structures. Wind induced vibration of conductors leading to the accidents like the breaking of insulators attachments, spacers and porcelain plates due to the large amplitude vibrations of high voltage electrical bundled conductors in gusty winds is observed in the field (Fig. 1). Any damage leading to interrupt electricity system was very costly and affects strongly to modern life.
To control wind induced vibration response of bundled conductors still depends on field trial and error procedure [1] [2]. Two types of the wind induced vibration response of high voltage electrical bundled conductors are the galloping and the gust response. The characteristics of galloping and gust response are entirely different, so that the way to control them would be different. It is, therefore, indispensable to identify the characteristics of high voltage electrical conductors.
For analysis of wind induced vibration response of electrical conductors developed in time, time domain gust analysis is considered. The feature of this approach is that the time dependent characteristics of unsteady aerodynamic forces and the nonlinearities of both aerodynamics and structural origins can be taken into account in estimating the vibration response electrical conductors. In the timedomain analysis, the time history of turbulence acting on structural nodes is input data. The turbulent field is considered as multivariate spatially correlated random Gaussian processes. It is difficult to obtain the time series of turbulence in the field at all points of the bundled conductor. Therefore, the digital simulation techniques are often used to generate the time history of turbulence. The
decomposition techniques as the Cholesky decomposition [3] and modal one have been applied. The buffeting response analysis of bridges in coupling with the aeroelastic forces in the time domain have been carried out [4]. The proper orthogonal decomposition (POD) has been applied in order reduced modeling of the random wind fields [5] [6]. The stochastic response of structures formulated in the time domain using the covariance proper transformation is very promising because of its direct and capable solutions for nonlinear problems and unsteady aerodynamics.
In this study, a finite element (FE) model of bundled conductors associated model of nodal wind forces is analyzed by using Direct transient response analysis solver. The comparison of fieldmeasured data and simulation of wind induced vibration is presented to investigate the response developed in time. The result may be used in designing and controlling damages due to wind induced vibration of bundled conductors.
Figure 1. Damages of high voltage electrical bundled conductors

MATERIALS AND METHODS

Finite element model of bundled conductor
In order to evaluate dynamic characteristics of high voltage electrical bundled conductor systems, the eigenvalue analysis is performed. The FE model of bundled conductor is created and analyzed by using the computing software [7]. In creating the FE models of bundled conductor, the important part is to find
out the proper static equilibrium configurations of conductors subjected to the body loads. Catenary theory is a popular solution for static equilibrium configuration of conductorss under its selfweight.
D t 1 U 2d n C
b c c D Du
b c c D Du
2
1
f 2ut C C
U D L Dw
2ut
f wt
U
wt
The details of FE models of bundled conductor are indicated in Fig. 2. The modeling of electrical conductors is
Mb t
U 2 d n 2 C f
2 c c M Mu U

CM Mw f
U (2)
most essential in the FE modeling of bundled conductor
where CL , CD , CM
are aerodynamic coefficients at a
systems. In this study, three degreeoffreedom, twonode tube elements with zero bending stiffness apply to the conductors. The spacers in the transmission line systems have a function to maintain the distance between conductors in a bundled conductors and each conductor is clamped at the end of spacer arm. The spacers are modeled by a twonode rigid tube element having the rigid connection with the conductors. The material properties of conductors and the spacer of the electrical bundled conductor system are shown in Table I.
Components
Parameter
Units
Values
Conductor
Elastic modulus
(109N/m2)
90.6
Diameter
(mm)
22.4
Spacing
(m)
0.4
Weight
(kg/m)
1.11
Spacer
Weight
(kg)
9.5
Elastic modulus
(109N/m2)
210
Components
Parameter
Units
Values
Conductor
Elastic modulus
(109N/m2)
90.6
Diameter
(mm)
22.4
Spacing
(m)
0.4
Weight
(kg/m)
1.11
Spacer
Weight
(kg)
9.5
Elastic modulus
(109N/m2)
210
Figure 2. FE models of the of high voltage electrical bundled conductors TABLE I. SPECIFICATION OF THE ELECTRICAL BUNDLED CONDUCTOR


Method of vibration response analysis in time domain
Two different numerical methods can be used for a gust response analysis: modal and direct response analysis methods. The direction method performs a numerical integration of the complete coupled equations of motion. The modal method utilizes the mode shapes of the structure to reduce and uncouple the equations of mtion. The solution is then obtained through the summation of the individual modal responses. In direct response analysis, structural response is computed by solving a set of coupled equations using direct numerical integration. In this study, the direction response method is used.
MultiDOF motion equations of electrical bundled conductor subjected to the buffeting forces is expressed:
balanced angle; CL , CD , CM are firstorder derivatives with respect to angle of attack; dc , nc are the diameter of the bundled conductor and number of bundled conductor; F F L, D, M; u, w is aerodynamic transfer
functions between turbulent components and forces.
In current study, Eq. (1) will be solved by computing software by using Direct transient response analysis method based on FE model associated with nodal turbulent wind force acting on the nodes of models
Direct transient response analysis is the most general method for computing the forced dynamic response. The purpose of a transient response analysis is to compute the behavior of a structure subjected to timevarying excitation. The transient excitation is explicitly defined in the time domain. All of the forces applied to the structure are known at each instant in time. The important results obtained from a transient analysis are typically nodal translation in three directions.
For Direct transient response analysis solver, important information is selected as follows. Nonlinear properties on geometry, large displacement is selected in Nastran bulk data options. Direct solution type is selected in Nastran modal analysis setup. Full length of wind velocity time series is recorded in 600 second with time step of 0.02 second. However, vibration response analyzed in time domain is recommended to simulate in 100 second because of long enough of analysis time, capacity of computer and easily selecting wind velocity and vibration response in comparison. Therefore, in transient time step intervals, 1000 steps, 0.1 second time step and one output interval per one step of loading are selected. Due to the cable type structure, the overall structural damping coefficient is selected to be 1%. Loading case in the analysis is combining of body load and nodal wind force.
Nodal wind force is calculated by using Eq. (2) with their magnitude of nodal timevarying wind force depended on number of applied nodes along the span. Variant turbulent processes u (t) and w (t) are selected from fieldmeasurement. In time domain analysis, nodal wind force varied in time is
Mu (x,t) Cu(x,t) Ku(x,t) Fb (t)
(1)
defined by using time depended function. In the current study,
where u(x, t) u(x, t) u (x, t) are deflection and its first, second order derivatives; M , C , K are mass, damping, stiffness matrices; Fb (t) is fullscale buffeting forces.
The buffeting forces per unit length of the bundled conductor can be expressed and corrected by the frequency dependent admittance functions as follows. (Eq.2).
nodal wind forces are applied uniformly to all nodes at different location of nodes.

Wind field simulation
In some cases, without fieldrecorded wind velocity, empirical auto Power Spectral Density (PSD) is very usable. Some empirical formulae can be applied to simulate the longitudinal turbulence u (t) and the vertical one w (t). Von
L t 1 U 2d n C f 2ut C C
f wt
Karman (1948) firstly introduced a practical empirical formula
b 2 c c L Lu U
L D Lw
U
of the auto power spectral density based on ongoing turbulent flow conditions such as mean velocity, intensities of turbulence, the length scale of turbulence and so on. This study
used Kaimals spectra (Eq.3) to simulate the longitudinal turbulence u (t) and the vertical one w (t).
affect the results significantly and therefore the value of k is taken as 15. The parameters are collected in Table II&III.
Su f
a nu 2
u
u
u
u
f (1 b n5c / 3
(3)
In preliminary time domain gust response analysis, along and vertical wind turbulence is selected from a fieldmeasured wind velocity. Based on Eq. (2), time series of nodal wind
Sw f
f
a nu 2
w
w
w
w
(1 b n5c / 3
force are applied uniformly to nodes.
Variable
Description
Units
Value
Air density
kg/m3
1.204
CD,L,M
Drag coefficient single conductor
1.0
Aerodynamic admittance
1
k
Empirical constant in frequency dependant correlation of wind velocity
15
Variable
Description
Units
Value
Air density
kg/m3
1.204
CD,L,M
Drag coefficient single conductor
1.0
Aerodynamic admittance
1
k
Empirical constant in frequency dependant correlation of wind velocity
15
TABLE II. ASSUMED PARAMETERS IN ANALYSIS.
where
n fz
U
is the nondimensional coordinates; f is
frequency variable (Hz); U, z are mean velocity (m/s) and

kU
altitude z (m), respectively;
u
ln z / z0
is friction or
shear velocity (m/s); k, length (m).
z0 are scale factor and roughness
TABLE III. AERODYNAMIC DAMPING RATIOS
Mode
Horizontal
Vertical
Torsion
1
5%
2%
2%
2
4%
2%
1%
3
4%
2%
1%
Mode
Horizontal
Vertical
Torsion
1
5%
2%
2%
2
4%
2%
1%
3
4%
2%
1%
The longitudinal u (t), vertical turbulent components w (t) are considered as two multivariate spatiallycorrelated Gaussian random processes with zero mean. The turbulent fields at n discrete nodes (Fig.3) are expressed as follows:
ut u1t , u2 t ,…,uN
t T
(4)
In this study, variate along and vertical turbulent processes
u (t) and w (t), respectively, are selected from the field
1 2 N
1 2 N
wt w t , w t ,…, w t T
Figure 3. Turbulent wind at n discrete nodes of conductor model


Vibration response simulation in time domain
Model of wind force is shown in Eq. (2). In the equation, there are some parameters assumed in the following analysis. The air density depends on the temperature, pressure and presence of moisture, and is 1.204 kg/m3 at 20oC under the standard atmospheric pressure and in the absence of moisture. In this analysis, its value is taken as 1.204 kg/m3 because the possible change in its value may not affect the results of analysis significantly. The force coefficients for electrical conductors is provided by the company. In this study, an average value of 1.0 is assumed. The aerodynamic admittance is introduced to counter the effect of flow mdification due to the motion of the structure on the force coefficients. Transmission line conductors, being very slender and moving with very small velocity as compared to the velocity of the oncoming wind, do not modify the flow significantly. Therefore, the aerodynamic admittance for the transmission lines is assumed as one. An appropriate value of k , which is
the empirical constant in the frequencydependent correlation of wind velocity, can be determined only through experimentation, but Holmes discusses on its general range from 10 to 20 [8]. The possible change in its value may not
measurement in order to simulate the nodal wind force that is employed as input of time domain gust response analysis. Fieldmeasured large amplitude fluctuation responses corresponding to those turbulences in two directions are also considered in comparison with simulation results. Scheme of time domain gust response analysis is shown in Fig. 4.
It should be noted that both vibration response and turbulent time series used in simulation and comparison extract from corresponding 600 second time series measured in the field. In this study, 100 secondtime series is used in simulation and comparison because time series of 100 seconds is long enough for structural vibration response dominated by studying excitations. In simulation, FE model with huge number of node and element as longspan, large bundledconductor transmission lines, long time series often leads to a big problem with analysis time, etc. Although simulated vibration responses in Fig. 4 is obtained from fieldmeasured turbulences (input data), model of wind force and FE model of the transmission line system may be a bit different from the real ones due to the assumptions. Therefore, it is difficult to get exact vibration response time series as fieldmeasured vibration response time series. The main objective of this analysis is to investigate the way of developing of large amplitude gust response in time. Based on the result, large amplitude of bundleconductors due to gusty wind can be predicted for further works as in designing, structural health monitoring (phase to phase flashovers), etc.
In comparison, some points of discussion can be given as follows. Firstly, the agreement between fieldmeasured response and simulation is unexpected at the beginning of 100 second vibration time series because lack of excitation before 100 second time series in simulation. Secondly, peaktopeak and time period of response selected in a simulation should be considered in same sampling frequency (10Hz) with field measured vibration response that may keep proper comparison.
Figure 4. Scheme of framework in time domain gust response analysis


RESULTS AND DISCUSSIONS

Results of the eigenvalue analyses of FE model
The eigenvalue analyses are then conducted for the FE models of transmission lines with proper static equilibrium configurations, and the results of eigenvalue analysis are summarized in Table IV. It is confirmed, in the case of horizontal mode, the lowest natural frequency corresponds to one loopperspan mode and that the higher natural frequencies are nearly equal to the integral multiple of the lowest one. Results of eigenvalue analysis present natural frequencies and mode shapes that are used to update by fieldmeasured vibration in frequency domain for a realizable FE model of the electrical bundled conductor.
TABLE IV. NATURAL FREQUENCIES AND MODE SHAPES
Mode
Horizontal
Vertical
Torsion
1
0.099 Hz
0.145 Hz
0.131 Hz
2
0.156 Hz
0.190 Hz
0.183 Hz
3
0.197 Hz
0.195 Hz
0.212 Hz

Results of wind field simulation
The fitting result to Kaimals spectra is conducted with a=220; b=45; c=0.95 with k=0.4; z1=0.01m; z2=100m;
u* 1.5 ; =0.2 (terrain condition). The along turbulence spectra comparisons are plotted in Fig.5.
Fig.6 shows the nodalwind time series at a sample grid
point. The meanwind velocity at the referred point is 19.4m/s (top of tower). It could be said that it is impossible to obtain the simulated wind field that is equal exactly to the measuredwind field because of the stochastic procedure. However, with the consideration of the fitting the spectra model to the measured spectra and the terrain condition of the site, the result of the wind simulation is assumed to be appropriate to the wind field.
Figure 5. Fitting result to Kaimals spectra
Figure 6. Wind time series at node

Vibration response analysis in time domain
In this analysis, the timedomain gustresponse analysis subjected to the simulatedwind in previous part is conducted by using the basic computing software solver. The time series of drag wind force, lift wind force and moment wind forces are calculated by using Eq. 2. Based on the length of the elements, the nodalwindforce time series is obtained. For example, by assuming parameters in Table V, the nodalwindforce time series at midspan is shown in Fig.7.
Vibration responses are simulated and compared with field measured vibration data at the same position of middle span. The time series of longitudinal, vertical and torsional compared with the fieldmeasured vibration responses turbulences are presented in Fig.8, Fig. 9 and Fig.10, respectively.
TABLE V. NODAL WIND FORCE PARAMETERS
(kg/m3)
U
(m/s)
dC
(mm)
nC
CD, CL, CM
CD, CL, CM
1.204
18.9
0.04
4
1
0.76, 0.05, 0.003
0.1, 0.006, 0.0002
Figure 7. Nodal wind force of subconductor at mid span

E+02

1.E+01
PSD (m2/Hz)
PSD (m2/Hz)
1.E+00
1.E01
1.E02
1.E03
1.E04
1.E05
1.E06

Field measurement Simulation
0.01 0.1 1 10
Frequency (Hz)
Figure 11. PSD of field and simulation fluctuation horizontal responses

E+02
Figure 8. Time series of Horizontal vibration response
1.E+01
PSD (m2/Hz)
PSD (m2/Hz)
1.E+00
1.E01
1.E02
1.E03
1.E04
1.E05
1.E06

Field measurement Simulation
Figure 9. Time series of Vertical vibration response
Figure 10. Time series of Torsional vibration response
In the figure, time series of simulated vibration responses of
3 directions seems to fit properly to time series of field measured vibration responses. The MPPA of 3 directions is similar between the field response and simulated response. The shape of response amplitudes is similar in some period of time series. Especially, the torsional response in Fig.10 shows very good agreement.
To check the agreement of horizontal and vertical response, PSD comparison is calculated and presented in Fig. 11 and Fig.12, respectively. In the figures, both horizontal and vertical PSDs show shapes look quite similar and close peak responses in spectra around 0.1 Hz of horizontal PSD and 0.2 Hz of vertical PSD as obtained in the results of eigenvalue analysis.
0.01 0.1 1 10
Frequency (Hz)
Figure 12. PSD of field and simulation fluctuation vertical responses


CONCLUSIONS AND REMARKS
The study was conducted to predict the wind induced vibration response of high voltage electrical bundled conductors in time domain based on a comparison between the fields measured data and simulated vibration response analysis. The results lead to the following conclusions:
To eal with the spatial windfield simulation for all nodes of the electrical conductors model, of which the nodal coordinates change in the two dimensions, the grid of finite locations of the windtime series that covers all of the nodes of the conductor bundle is used.
By doing so, the wind time series at any coordinate can be easily obtained. Additionally, the grid of nodal wind time series makes the advance NX Nastran solver become more convenient due to interpolating the coordinates of the structural and aerodynamic model.
FE model of the high voltage electrical bundled conductors associated model of the spatial wind is analyzed in the NX Nastran computing software by using the basic solver, the nonlineardirect transient response analysis.
The simulatedresponse time series seem to be understandable. However, the basic solver is inconvenient for the aerodynamic analysis, especially for setting the damping matrix. Additionally, the structural node will change in time due to the continuousmovement of conductor in space leading change of the attacked angle as well as force coefficient.
In fact, it is not easy to simulate exactly the spatialwind that can induce the vibration response of the transmission lines as observed in the field. However, by using the timedomain gustresponse analysis, the timedependent characteristics of unsteady aerodynamic forces and nonlinearities of the both aerodynamic and structural origins can be taken into account in characterizing the type of wind induced vibration response of high voltage electrical bundled conductors. The result may be used in designing and controlling damages due to wind induced vibration of bundled conductors.
REFERENCES

Ohkuma T., Marukawa H., 1999. Galloping of overhead transmission lines in gusty wind. First International Conference on Advances in Structural Engineering and Mechanics, Seoul (Korea), pp. 71 76.

CIGRE. 2007. State of the Art of Conductor Galloping. Technical Brochure 322. Task Force B2.11.06.

D. M. Hawkins and W. J. R. Eplett, The Cholesky factorization of the inverse correlation or covariance matrix in multiple regression, Technometrics, vol. 24, no. 3, pp. 191198, 1982.

Chen, X., Matsumoto, M., Kareem, A. (2000), Time domain flutter and buffeting response analysis of bridges, J. Engrg. Mech., ASCE 126(1), pp. 716.

Solari, G., Tubino, F. (2005), Gust response of long span bridges by double modal transformation, The 4th Euro. and Afri. Conf. on Wind Eng. (EACWE4), Praque, Czech.

Matsumoto, M., Shirato, H., Le, T.H. (2007), Time and frequency domain gust response prediction using proper orthogonal decomposition, The 12th Intl Conference on Wind Engineering (ICWE12), Cairne Australia, July pp. 16.

FEMAP with NX Nastran. User guide. 2010

Holme J.D., 2001. Wind loading of structures. London Spon Press

Matsumoto, M., Chen, X. (1996), Time domain analytical buffeting responses for longspan bridges, Proc. 14th National Symposium on Wind Engineering, JAWE pp. 515520.

Matsumoto, M., Shirato, H., Le, T.H. (2007), Time and frequency domain gust response prediction using proper orthogonal decomposition, The 12th Intl Conference on Wind Engineering (ICWE12), Cairne Australia, July pp. 16.