Open access peer-reviewed chapter

# High Voltage Transmission Line Vibration: Using MATLAB to Implement the Finite Element Model of a Wind-Induced Power-Line Conductor Vibration

By Chiemela Onunka and Evans Eshiemogie Ojo

Submitted: September 6th 2017Reviewed: February 11th 2018Published: September 19th 2018

DOI: 10.5772/intechopen.75186

## Abstract

Wind-induced vibration affects the performance and structural integrity of high voltage transmission lines. The finite element method (FEM) is employed to investigate wind-induced vibration in MATLAB. First, the FEM model was used to develop the equation of motion of the power line conductor. In addition, dampers, conditions for damping, free and forced vibrations of the overhead conductor were considered in the FEM model. Wind-induced experiments were conducted in the laboratory using an actual overhead power conductor. The developed FEM models were simulated in the MATLAB computing environment. The results from the MATLAB simulation, finite element and experimental recordings were compared in order to evaluate the efficacy of models simulated in MATLAB and developed using the FEM.

### Keywords

• aeolian vibration
• power conductor damping
• resonant frequency
• MATLAB

## 1. Introduction

The availability and use of electrical power in the society is crucial in the development and growth of the society. Power generated from power stations is transmitted using high voltage transmission lines. The transmission of power from the point of generation to the point of use requires complex network of high voltage lines, systems and components . High voltage conductors are usually subjected to vibration and the vulnerabilities of the power lines to vibration can lead to fatigue failure. Thus, power loading determination and control on the power grid can influence the integrity of the transmission network. High voltage conductor vibration is very difficult to model due to the fact that the responses exhibit a non-linear behaviour. There has been concerted effort to try and predict the conductor response as a result of aeolian vibration. Evaluations of conductor vibration caused by aeolian forces have been investigated be several researchers [2, 3, 4, 5, 6].

Recent researches developed models to investigate wind-induced vibration using nonlinear time history, expert systems, the concept of principal modes, aero-elastic and bending stiffness [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. The investigation of vortex formation and the phenomenon of wind-induced vibration was done using the concept fluid–solid dynamic excitation [18, 19, 20]. The models were used to determine how wind loading influenced the oscillation of transmission lines. This form of investigation was done by experimental studies carried out in a wind tunnel [21, 22]. The outcomes of these experiments were used to determine conductor input loading. Several models developed by various researchers can be used to determine conductor damping and also the placement of vibration absorbers on the line conductors to curtail the effect of cable mechanical oscillation [23, 24, 25, 26, 27]. Based on the various models that have been developed by researchers as indicated in the first and second paragraphs, there is a need to further analyze wind-induced vibration using finite element method (FEM) in MATLAB.

The design, construction and maintenance of power transmission network requires adequate understanding of the system dynamics that occurs when subjected to vortex induced vibration . Various analysis can be conducted using techniques that suit certain objectives. System integrity in high voltage transmission lines is of paramount importance. MATLAB is a multi-model simulation environment used for numerical computing. It can be integrated with physical hardware or systems in order to determine real-time performance, characteristics and behavior. MATLAB also provides a platform for special hardware in loop simulations . These functionalities amongst others are vital in determining various characteristics and behaviors in high voltage transmission lines.

High voltage transmission lines and grid can experience vulnerabilities such as vibration, electromagnetic transients, fatigue, transmission loss, switching surges, conductor sag fluctuation [30, 31]. When the conductor experiences vibration, the transmission lines experience high amplitudes of vibrations from wind forces and can lead to fatigue of the transmission lines . The use of systems simulation and analysis provides the platform to understand the response of the transmission conductor. The responses considered in the chapter include transmission line excitation through wind loading, conductor properties such as damping and damper placement used in mitigating the vibration.

The chapter discussed the development and implementation of a wind-induced high voltage transmission line vibration using finite element method (FEM) in MATLAB. The sections in the chapter discussed the development of transmission line equation of motion, the solution to the equation of motion, free and forced vibration of the transmission line, dampers and conductor self-damping, FEM MATLAB setup and implementation, simulation of FEM models. The chapter also discussed results from FEM models, simulation and experimental investigation. The chapter is focused towards the development of a finite element method and its implementation on the MATLAB software. The developed finite element method (FEM) was based on the concept of the simply supported beam model and it was used in modeling the transverse vibration of power line conductors. The results from the FEM were then compared with results from the analytical model and results obtained from experimental studies documented in . The results from MATLAB simulations from the finite element models and experimental results were compared in order to determine the accuracy of the models. The developed FEM was then used as the means to verify the effect of varying the conductor axial tension on the natural frequencies of the conductors.

## 2. Transmission line equation of motion (EOM)

The transverse displacement of high voltage transmission line conductor is generally caused by wind loading. This form of vibration with small displacement is known as aeolian vibration and it is a source of concern to the power lines reliability. One the vulnerabilities is that it can cause fatigue failure of the transmission lines. Conductors are example of continuous or distributed systems and modeling its mechanical vibration can either be as a beam or taut string. In [18, 19], it was ascertained that modeling a conductor as a beam is more accurate than modeling it as a taut string due to the effect of the bending stiffness. Hence, in line with the above, the conductor transverse vibration was modeled as a beam, simply supported or pinned at both ends. The distributed loading on the conductor is replaced by effective point load that can effectively have the same resultant effect as that of the actual distributed load.

The high voltage transmission line equation of motion was formulated by assuming that power conductors can modeled as beams with fixed ends. The following assumptions were considered :

• The power conductor is uniform along its length and it is slender

• The power conductor is a solid with cylindrical body having both linear and homogeneous physical properties throughout its cross-sectional area

• The power conductor has a symmetrical plane which acts as the plane of vibration such that there is the decoupling of translational and rotational motion.

The assumptions were based on beam theory. In considering the power conductor as a beam, sagged by a tensile force S, being acted upon by a concentrated wind load fxt, with cross-sectional area A, density ρ, flexural rigidity EI, displaced at a distance of xafter time. In Eq. (1), the high voltage transmission line equation of motion is expressed as:

fxt=EI4yxtx4S2yxtx2+ρA2yxtt2E1

For x01,t>0. The boundary conditions are expressed and indicated in Eqs. (2) and (3):

y0t=20tx2=0E2
ylt=2ltx2=0E3

The initial conditions at t=0are indicated in Eq. (4), Eq. (5) and expressed as:

yx0=yoxE4
ẏx0=ẏxE5

Introducing the mass per unit length of the power conductor, the new equation of motion indicated in Eq. (6) is expressed as:

fxt=EI4yxtx4S2yxtx2+m2yxtt2E6

In order to derive a possible solution, the model was simplified using dimensionless functions and Dirac delta functions. In Eqs. (7)(12), the variables are expressed in dimensionless form and expressed as:

Y=yxtDE7
X=xLE8
τ=tfE9
Ip=Df2gE10
Sp=SDγL2E11
MpEIDγL4E12

Eq. (13) indicates the revised equation of motion and it is expressed as:

Mp.4YX4Sp2YX2+Ip2Yτ2=1γFXτ+nδXXnFnτE13

Where γrepresents the power conductor weight per unit length and grepresents gravitational constant. XnδXXnrepresents the Dirac delta function, FXτdenotes the net transverse force per unit length acting on the conductor and Fnτdenotes the nth concentrated force acting transversely on the conductor.

## 3. Solution to the EOM

The general solution to the high voltage transmission line equation of motion was derived using Euler-Bernoulli equation. The particular solution to the equation of motion was derived using a product of two functions. The two functions were first separated using the principle of variable separation as expressed in Eq. (14) :

Yxt=XxTtE14

Where Xxis the normalized function representing the mode shape of the equation of motion. The normalized function ensures that orthogonality condition was satisfied in the derivation of the EOM model solution. Applying the normalized function in the EOM yields Eqs. (15) and (16):

EIX////xSX//xω2ρAXx=0E15
T¨t+ω2Tt=0E16

Where X////x=d4ydx4, X//x=d2ydx2, T¨t=d2ydt2and ω2is a constant that equates xand t. Assuming that Xx=ZeΨx, the model is expressed in Eq. (17) as:

ZeΨxEIΨ4SΨ2ρAω2=0E17

Considering that ZeΨx0, hence EIΨ4SΨ2+ρAω2=0. The general solution of the Euler-Bernoulli equation which represents the solution to the equation of the motion of the transmission line is expressed in Eqs. (18) and (19) as :

Ω2,Ψ2=S±S24EIρAω22EIE18
Ω,Ψ=±S±S2+4EI(ρAω2)2EIE19

The values of Ωand Ψrepresents the general solution of the equation of motion. The practical implication of the derived solution is that it represents the transverse vibration of the high voltage transmission line. The derived solution has infinite number of solutions and the solution is indexed to accommodate all the possible solutions from the model. The indexed solution is expressed in Eqs. (20) and (21) as:

Ωn=S2EI+S22EI2+mL2πfn2EIE20
Ψn=S2EI+S22EI2+mL2πfn2EIE21

Where ωn=2πfnand for n=1,2,3,.

In Eqs. (22)(24), the infinite natural frequencies of the power conductor were derived while considering that the mode shape is the same as a pinned-pinned beam eigenfunction model with no external force. Hence,

Ynxt=sinnπxlcosωntE22
EIl4sinnπxlcosωntSl2sinnπxlcosωnt+ρAωnsinnπxlcosωnt=0E23
sinnπxlcosωntEIρAl4+SρAl2ωn2=0E24

The natural frequency of the power conductor in rad/s is expressed in Eqs. (25) and (26) as:

ωn2=l2S+l4EIE25
ωn=L2SmL1+L2EISE26

The natural frequency in Hz is expressed in Eq. (27) as:

Fn=12πL2SmL1+L2EISE27

## 4. Free vibration of power conductor

The self-damping model of the power conductor provided the basis to analyze free vibration experienced by the conductor. Free vibration occurs when the forcing function causing the power conductor to vibrate become zero. Hence the equation of motion is expressed in Eq. (28) as :

EI4yxtx4S2yxtx2+βI5yxtx4t+Cyxtt+ρA2yxtt2=0E28

Applying the principle of separation of variable to the equation of motion yields Eq. (29):

EIX////xTtSX//xTt+βIX///xṪt+CXxṪt+ρATtXx=0E29

Integrating the eigenfunction Xnx=sinnπxlin the model yields Eqs. (30)(33):

EIl4sinnπxlTtSlsinnπxlTt+βIl4sinnπxlṪt+CsinnπxlṪt+ρAsinnπxlT¨t=0E30
sinnπxlEIl4Tt+Sl2Ṫt+βIl4Tt+CṪt+ρAT¨t=0E31
ρAT¨t+βIl4+CṪt+Sl2+EIlTt=0E32
T¨t+βIρAl4+CρAṪt+SρAl2+EIρAl4Tt=0E33

Considering that the vibration model represents a multi-degree vibration system. The natural frequency of the power conductor is determined and expressed in Eqs. (34) and (35):

ωn2=SρAl2+EIρAl4E34
2ξωn2=βIρAl4+CρAE35

The temporal solution to the free vibration model is expressed in Eq. (36) as:

Tn=A1eξnωntsinωdt+ϕE36

The solution can also be represented in Eq. (37) and expressed as:

Tn=eξnωntB1sinωdt+B2cosωdtE37

Where the damped frequency of the power conductor is expressed in Eq. (38) as:

ωd=ωn1ξ2E38

The system response is expressed in Eq. (39) as:

yxt=n=1A1eξnωntsinωdt+φsinnπxlE39

The response can also be represented in Eq. (40) and expressed as:

yxt=n=1eξnωntB1sinωdt+B2cosωdsinnπxlE40

## 5. Forced vibration of power conductor

High voltage transmission lines are exposed to loading from the wind. The actual system representation through system simulation strategy considers a case of distributed load through the span of the conductor. In order to simplify simulations, the external force acting on the conductor is represented as a point load. In Eqs. (41)(43), the equation of motion is solved with an excitation force in order to evaluate the actual response of high voltage transmission lines under aeolian vibration . Hence,

EI4yxtx4S2yxtx2+βI5yxtx4t+Cyxtt+ρA2yxtt2=fxtE41
sinnπxlEIl4Tt+Sl2Ṫt+βIl4Tt+CṪt+ρAT¨t=FsinωdrtE42
T¨t+βIρAl4+CρAṪt+SρAl2+EIρAl4Tt=FsinωdrtE43

Expressing the model as a multi-degree system yields Eq. (44):

Tt=Aeζωntsinωdt+ϕ+XcosωtθE44

The natural frequency of the power conductor under forced vibration is expressed in Eqs. (45) and (46) as:

ωn2=SρAl2+EIρAl4E45
2ξωn2=βIρAl4+CρAE46

The solution to the equation of motion under forced vibration is expressed in Eq. (47) as:

yxt=sinnπxlAeζωtsinsinωt+ϕ+XcosωtθE47

## 6. Conductor self-damping and dampers

The influence of external and internal damping mechanisms was considered in the conductor vibration model. The factors considered included the following [1, 32]:

• The power conductor inter-strand motion and fluid damping. This is proportional to the conductor velocity and represented as viscous damping in the model.

• The rate of strain in the power conductor. This proportional to the internal damping of the power conductor.

The high voltage transmission line damped model is expressed in Eq. (48) as:

EI4yxtx4S2yxtx2+βI5yxtx4t+Cyxtt+ρA2yxtt2=fxtE48

Where Cand βrepresent damping constants. In the presence of axial load, viscous air damping, strain rate damping or Kelvin-Voigt damping, high voltage transmission line integrity can be managed.

There are various types of dampers that can be used to reduce vibration. The dampers are excited by the vibration of the power conductor and the vibration of their masses connected by the massager cable help to damp out energy. Stockbridge dampers are commonly installed on high voltage transmission lines to reduce aeolian vibrations. Stockbridge dampers can be symmetrical or asymmetrical in their design. An example of dampers installed on high voltage transmission lines is shown in Figure 1 . The design of Stockbridge dampers follows the principle of cantilever beams with mass at the free ends. The contribution of dampers to power conductor vibration mitigation is to lower the severity of the vibration to a level that might prevent failure to the line.

## 7. FEM MATLAB model setup, formulation and implementation

In order to implement the conductor model in MATLAB environment, finite element analysis formulation was done as function of the physical state of power transmission line conductor. The models developed using finite element analysis can then be implemented in MATLAB. The FEM model enables the analysis of the dynamic behavior and response of power line conductor to the dynamic forces of wind . Consider a power transmission line subjected to dynamic aeolian vibration as an assembly of thin strands having distributed mass and elasticity. The physical model can be represented used partial differential equations. Each strand in the transmission line experiences axial, bending and torsional loads from the wind [34, 35]. The accurate representation of each factor is critical in the determination of the dynamic behavior of power transmission lines [36, 37]. Euler-Bernoulli curved beam theory was used to formulate the finite element model of power transmission lines.

Consider a power transmission line experiencing a vertical force, curvature and an axial force has an axial displacement modeled in Eq. (49) as :

uxy=u0x+vRyθx+uxRE49

Where θxrepresents the rotation of the power line due to flexural effect, Rrepresents the radius of rotation, yrepresents the distance from the axis of rotation to the centroidal axis of the conductor or transverse displacement, vrepresents tangential displacement and uxrepresents the axial displacement of the power lines. The shape function for power transmission lines having rotation, bending and axial motion components is modeled using discretization techniques and represented in Eqs. (50)(52) as:

us=b0+b1sE50
vs=c0+c1S+c2S2+c3S3E51
θ=dv(S)ds=c1+2c2S+3c3S2E52

The solution to the discretization of the models yields Eq. (53) to Eq. (59):

uvθ=N10N20N3000N4000N4N5N60N5N6u1v1θ1u2v2θ2E53

Where

N1=121ζE54
N2=121+ζE55
N3=1423ζ+ζ3E56
N4=141ζζ2ζ3E57
N5=142+3ζζ3E58
N6=141ζ+ζ2+ζ3E59

The power line matrix model contains the strand stiffness K, mass matrix Mand the load vector F. They are expressed in Eq. (60) as:

K=12NaTEANaδζ+12NBTEINBδζ+12NBTTNB,TδζE60

Where Arepresents the cross-sectional area of the power line, Erepresents the young modulus of the power line material, Irepresents polar moment of area, Trepresents the kinetic energy of the system. The matrix is modeled in Eq. (61).

M=12u̇TρAu̇+12ġTρAv̇E61

Where ρrepresents the density of the power line material and Eq. (62) indicates external excitation.

δW=12FδuE62

The power line conductor model is constructed using Euler-Bernoulli theories and summarized in Eq. (63) as:

M11M12M21M22u¨v¨+K11K12K21K22uv=F1F2E63

The finite element analysis follows a step by step numerical computation in the MATLAB environment as documented in [38, 39]. The dynamic response analysis assumes continuous displacement, velocity and acceleration [40, 41]. The numerical integration technique utilized was based on Newmark integration method. The compact form of the high voltage transmission line model is expressed in Eqs. (64)(69) as :

My¨+Cẏ+Ky=FE64
Ks+1=Ks+1+a3ms+1E65
Fs,s+1=Fs+1+ms+1ays+aẏs+ay¨sE66
a3=2γΔt2E67
a4=2γΔtE68
a5=1γ1E69

The initial conditions are expressed in Eq. (70) as:

y¨0=M1F0K1y0E70

The acceleration vector is expressed in Eqs. (71)(74) as:

y¨s+1=a3ys+1ysa4ẏa5y¨sE71
y¨s+1=y¨s+a2y¨s=a1y¨s+1E72
a1=αΔtE73
a2=1αΔtE74

## 8. FEM MATLAB model implementation strategy

In order to test the validity of the models discussed earlier using MATLAB, an aluminum power conductor with a steel core having a total diameter of 35.56 mm and having an ultimate tensile strength of 51.51kN was used in setting up the MATLAB simulation. Further physical properties of the power cable are shown in Table 1 . The power conductor had a minimum bending stiffness EIminof 8.66 Nm2 and maximum bending stiffness EImaxof 433 Nm2. The wholistic finite element models where implemented in MATLAB using strategy expressed in Figure 2 .

Strand layerStrand materialDiameter (mm)No. of strandsPitch per length (cm)Lay direction
Layer 0Steel2.251
Layer 1Aluminum3.38616.1Left hand lay
Layer 2Aluminum3.381222.2Right hand lay

### Table 1.

Power transmission conductor physical properties.

The inputs in the MATLAB algorithm were bending and axial loads, the cross-sectional area of the power conductor, strand radius and strand material type. The type of analysis which can be either static or dynamic was also specified as part initial and boundary conditions. Also included in the algorithm was to specify if the computation focuses on local vibration of the power conductor or the global vibration model.

## 9. Experimental investigation of conductor vibration

MATLAB code was written for the FEM and this was used to model the dynamic analysis of the problem of conductor vibration. To validate the FEM model an experimental study was conducted at the Vibration and Research Testing Centre (VRTC) situated at the University of KwaZulu-Natal which comprises of apparatus similar to that shown in Figure 3 . The sweep tests (resonance search) were carried out and the test results were used to obtain natural frequencies and the modes of vibration for a Pelican conductor. The frequency range for the Pelican conductor was between 5 and 50 Hz and testing was done for three axial tensions of 20, 25, 30 and 35% of its ultimate tensile strength (UTS). The experimental results obtained were used to validate the developed FEM model. The comparison between results from the experimental data, FEM and the theoretical model for the three different axial tensions for high voltage conductors are reported in the next section.

## 10. Simulation and experimental results

The results from the MATLAB simulations were compared with results from the finite element models (FEM) and experimental recordings. These are shown in Figures 4 7 . The results were compared in terms of the natural frequency of vibration or vertical displacement of the power conductor.

## 11. Conclusion

The results showed that the implementation of the derived models in MATLAB provided a reliable strategy in the determination of the wind-induced dynamic properties of high voltage transmission lines. The results from MATLAB simulation, finite element method and experimental recordings were similar in values and showed similar trend. MATLAB as an environment can be used as a reliable simulation tool to implement and analyze high voltage conductor dynamics. The parameters obtained from the results, to some degree of accuracy can be used to predict the response of conductors due to aeolian vibration caused by wind loading.

## More

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Chiemela Onunka and Evans Eshiemogie Ojo (September 19th 2018). High Voltage Transmission Line Vibration: Using MATLAB to Implement the Finite Element Model of a Wind-Induced Power-Line Conductor Vibration, MATLAB - Professional Applications in Power System, Ali Saghafinia, IntechOpen, DOI: 10.5772/intechopen.75186. Available from:

### Related Content

Next chapter

#### Introductory Chapter: MATLAB Applications in Power System

By Ali Saghafinia

#### Recent Developments on Power Inverters

Edited by Ali Saghafinia

First chapter

#### Two‐Phase Inverters with Minimum Switching Devices

By Branislav Dobrucky, Tomas Laskody and Roman Konarik

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books