# Full-Scale Reconstruction for Transmission Line Galloping Curves Based on Attitudes Sensors

Authors: Decheng Wu; Shizhong Yang; Lu Tao; Hailin Cao; Xiaoheng Tan; Li Jing; Wenjian Qin
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Caused by strong winds or nonuniform icing, conductor galloping is one of the major hazards that should be monitored in a timely fashion. In this paper, we proposed a new full-scale reconstruction method for transmission line curves based on the attitude sensors that uses only the tangential information and arc-length constraint. Meanwhile, a comparatively low-complexity algorithm is presented:(1) quaternion rotation to generate the tangent vectors; (2) a modified B-spline global interpolation to evaluate the curvature information; and (3) a linear mapping method to correct reconstruction curves. Moreover, sensors placement method is put forward according to the galloping features. Finally, both simulation and experimental results demonstrate that the proposed method can accurately reconstruct the galloping curves in a relatively short time, which fulfills the requirement of the real-time galloping monitoring system.


---

## Body

## 1. Introduction

The overhead transmission line is the primary approach for long-distance electrical energy transmission, which is always challenged by the harsh environmental conditions. Caused by strong winds or nonuniform icing, conductor galloping is a kind of self-excited vibration with low frequency (0Hz-3Hz) and high amplitude (a few meters to tens of meters), and that can cause various types of damage, such as conductor fracture, tower collapse, and large-scale grid paralysis. Hence, galloping prevention is imperative [1–3].Since the famous galloping vertical and torsional mechanism was proposed by Den Hartog [4] and Nigol [5], extensive studies have been performed on the galloping phenomenon. Recently, several researchers have theoretically investigated the galloping mechanism, and mathematical models have been presented to simulate the dynamical behavior of transmission [6–9]. Other authors have attempted to establish equations of motion to describe the observed three-dimensional galloping behavior on iced transmission lines [10, 11]. Moreover, Finite Element Method-based (FEM-based) approaches [12, 13] with various analytical tools, such as ABAQUS [14] and FLUENT [15], have been applied in numerical investigations on aerodynamic galloping problems. Furthermore, full-scale test transmission lines [1, 3] and wind tunnels [16] have been established to test the galloping behavior and antigalloping techniques, but some limitations remain. Although such studies have been performed in both theory and practice, there is no universal theory to explain the galloping mechanism.Galloping data play an important role in the mechanism and effect analysis of dancing prevention. They can be crucial means to verify the antigalloping design theory and determine the galloping prevention plan and effect [17]. Unlike theoretical studies on the galloping mechanism, many authors have focused on the operation-state monitoring for the transmission line using various techniques and sensor devices. The maximum galloping amplitudes [3, 18] and conductor sags [19] are the most common parameters to investigate dancing. Using the differential global positioning system (DGPS) to achieve the real measurement of the conductor position is probably the most direct method for obtaining conductor sag information, but it lacks full-scale galloping information [19]. Video- and image-processing [20] and photography methods [17, 21] have been explored to detect the spatial parameters of transmission lines. However, they are easily affected by rough weather conditions such as dense fog or heavy rain. Moreover, using video cameras for transmission line galloping requires many storages and transmission resources. A new fiber Bragg grating (FBG) sensor was applied to achieve the real-time horizontal tension monitoring of a transmission line, and the galloping amplitude could be calculated from the relationship between the galloping amplitude and the horizontal tension [22]. A new noncontact monitoring technology of operation-sate monitoring for high-voltage transmission lines is presented and uses the correlation between the magnetic-field variation and the electrical and spatial parameters of the transmission line [23]. A sensor module for galloping monitoring was designed using the inertial measurement unit (IMU) and wireless communication, and it can measure the conductor displacements along three directions [24]. Nevertheless, most existing studies are limited to a single point of the single conductor, which is difficult to cover the full-scale transmission line intuitively.Unlike the aforementioned methods, this paper focuses on the full-scale reconstruction of the transmission lines. We aim to reconstruct the galloping curves via tangential information with a length constraint using few attitude sensors, which do not need to know the spatial position of the conductor obtained using DGPS technology. There have been a few studies of shape sensing using attitude sensors, such as Morphosense (a flexible ribbon-like device) [25] and Human wearing vest (with integrated shape sensing sensor array) [26]. The author in [25] proposed a spatial interpolation with a Pythagorean-hodograph, which has better reconstruction performance than does the natural cubic spline method. However, this method requires extensive iterative computations, which is not suitable for a real-time galloping monitoring system. Reference [26] presented a fast computation algorithm for the global shape reconstruction. Since the reconstruction precision is notably limited, many sensors must be used to ensure the reconstruction performance. Similarly, this method is not suitable for transmission line monitoring system because of the requirement of high sensor density. In this paper, we reconstruct the galloping curves by using a few attitude sensors, and a low-complexity algorithm is designed for the curve reconstruction. Using an attitude sensor that contains 9 axes (3 axes of gyroscope, accelerometer, and magnetometer), we can obtain conductor tangential information. Here, we only know the dip angle and length information of the conductor from a few attitude sensors. It should be noted that this problem is not a simple spline interpolation problem because the absolute spatial locations of the sensor nodes are not known. Because all directions of the 3-dimensional unit tangent vectors form a sphere, which is known as the Gaussian sphere, a simple and reasonable assumption of the tangent vectors between consecutive sensors is the geodesic of the given sensors. We use the quaternion rotation operation to fulfill the tangent information between the sensors. The “angular velocity” of the rotation, i.e., the curvature of the transmission line curves, remains unknown. With the curvilinear distances and rotation angles between successive sensors, a modified B-spline global interpolation is introduced to evaluate the curvature of the transmission line. Then, when the global tangent information and length parameter are combined, the desired galloping curves are reconstructed by integral operation. Moreover, due to the integral accumulated error effect, the reconstruction curve may cause a bias to the actual one, and this effect increases with the increase in integral length. According to the prior position information of the transmission tower, a linear mapping method is applied for curve correction, which can fix the reconstruction better to fit the real one. The reconstructed galloping curves can offer a visual method to monitor the galloping and extract galloping parameters such as amplitude, frequency, and waveform.We believe, to the best of our knowledge, that this paper is the first study to attempt to examine the full reconstruction of galloping curves from a set of tangential information with the arc-length constraints.


## 2. Transmission Line Galloping Monitoring System

The transmission line galloping monitoring system [27] is shown in Figure 1. The attitude sensors are deployed on the transmission line to monitor the angle information at a predefined rate periodically. The sensor nodes disseminate the acquired data in real time or timed to the tower monitoring extension using the wireless network. The transmission line galloping includes the horizontal swing, vertical vibration, and axial torsion, which correspond to displacement in three orthogonal directions (x-, y-, and z-axes). Each attitude sensor can provide three-angle information, i.e., the pitch, roll, and yaw of the current position, and each has a separate controller and a power supply. Also, the number of sensors depends on the span of the transmission line and measurement accuracy, which should meet the requirement that there is not more than one galloping wave between the adjacent sensors. The tower monitoring extension gathers the uploaded information from each node, analyzes the data, and makes the appropriate data processing. The valid output data are sent to the monitoring center via the remote communication technology (Internet). Then, the full-scale transmission line curve is reconstructed in the monitoring center, and the galloping parameters (amplitude, frequency, and waveform) can be analyzed from the reconstructed data.Figure 1
Diagram of the transmission line galloping monitoring system with six sensors.Obviously, the design of a sensor module should consider various factors: (1) electrical properties, e.g., visual corona and radio interference, resistance of short-circuit current rush, indirect lightning stroke, temperature rising property, and impulse stability; (2) electromagnetic compatibility, such as electrostatic discharge, radiated susceptibility, pulse magnetic field immunity, and power frequency magnetic field immunity; (3) climatic protection, such as high-temperature performance, low-temperature performance, and hydrothermal properties; (4) mechanical characteristics, including vibration-resistance performance antifatigue property [27]. Recently, many monitoring sensors have been designed and put into actual test for transmission line galloping monitoring. The authors in [28] have developed a cable inspection robot to acquire light detection and ranging (LiDAR) data. Also, the authors in [24] have designed a wireless sensor module with smooth spherical structure shell, and it contains an inertial measurement unit and wireless communication components. Since this is an original study; this paper focuses on the reconstruction method via tangential information with a length constraint using few attitude sensors.We formulate the problem to design a space curve C(s):(1)Cs=P0+∫0stsds,s∈0,Lswhere P0 is the initial point (location of transmission tower), Ls is the length of the transmission line, and ts is the tangent vector with the arc-length parameter s. Obviously, we must reasonably calculate ts with the dip angle information from the attitude sensors.


## 3. Galloping Curve Reconstruction Method

The curve reconstruction is not merely a classical interpolation problem because we do not know the position of the sensor nodes. Moreover, we cannot independently address the dip information in each dimension because the spatial tangent vector is an entity whose norm is 1. Fortunately, because all three-dimensional tangent vectors lie on the unit sphere [29, 30], we can always find a curve on the sphere to denote the tangent vectors between the arc joining sensor nodes. In this paper, we interpolate the tangent vectors using the quaternion rotation.


### 3.1. Tangential Information with Arc-Length Constraint

Historically, quaternions were conceived by Hamilton as extensions to complex numbers. A quaternionq is defined by(2)q=w+xi+yj+zk=x,y,z,w=v,wwhere v=x,y,z∈R3, w∈R, and i, j, and k are different imaginary units with i2=j2=k2=ijk=-1. More detailed properties of quaternions can be found in [31].Given a unit quaternionq∈S3, i.e., Nq=qq∗=1, particularly if q=cos⁡θ+v0sin⁡θ∈S3 for an angle θ and an arbitrary vector t=x,y,z,0, the 3D rotation operation is defined as(3)Qt=qtq-1where Qt denotes the rotation by 2θ around the axis v0.For example, letv0 be a unit vector with inclination (or elevation) angle α and azimuth angle β. Then, q=sin⁡αcos⁡βsin⁡θ,sin⁡αsin⁡βsin⁡θ,cos⁡αsin⁡θ,cos⁡θ, and t2 is the vector after the rotation operation; i.e., t2=Qt0, as shown in Figure 2. Essentially, the rotation Qt=qtq-1 is divided into two steps. First, t1 is obtained from a θ rotation around the axis v0 from t0; i.e., t1=qt0. Second, t2 is obtained from a θ rotation around the axis v0 from t1; i.e., t2=t1q-1.Figure 2
Geometry of rotation for a given axisv0 and vector t0.Now, we return to the problem ofts interpolation. By exploring the attitude sensor, we can obtain the 3-dimensional angle information, which can be converted into tangential information. The tangent vectors are unit vectors that lie on the Gaussian sphere.A simple and possible answer is that these unknown tangential coefficients lie on the great arc between pointst0 and t1, as shown in Figure 3. This arc is named the geodesic curve and is in S2 instead of S3, which is given by(4)γθt=qθtt0,θt∈0,θ0,1(5)qθt=v0sin⁡θt,cos⁡θt=exp⁡v0θtwhere θt0,t1=cos-1t0∙t1∈0∘,180∘ is the angle between vectors t0 and t1. v0=t0×t1/t0×t1 is a unit vector in the direction of the cross-product of t0 and t1.Figure 3
Geodesic curve for given tangent vectorst0 and t1.According to the arc-length constraint,θt is a variable that is associated with the arc-length parameters s, i.e., θt(s), and it satisfies θt0=0 and θtlt0,t1=θt0,t1, where lt0,t1 is the arc-length between contiguous sensor nodes. Now, we obtain t(s) between the first two sensor nodes as follows:(6)ts=γθts=qθtst0,s∈0,lt0,t1.


### 3.2. Curvature Information Using B-Spline Interpolation

From (6), we obtain the curvature κs:(7)κs=dtsds=q′θtst0θt′s=q′θtst0·θt′s=θt′s,s∈0,l0,1.Sinceq′(θt(s)) and t0 are unit quaternions, the multiplication q′(θt(s))t0=1. Therefore, the curvature κs depends on θt′(s) or θt(s).A simple example isθt(s)=θ0,1/l0,1s, s∈0,l0,1. Apparently, κs=θt′(s)=θ0,1/l0,1 is a constant, and the resulting curve (1) is a circular arc, which is not suitable for the situation of the transmission galloping curve.Now, the calculation ofts is transformed into the problem of evaluating θt(s). Since the transmission line is an entire conductor, there is an internal connection between successive sensor nodes. Thus, a modified form of B-spline interpolation is introduced to estimate the curvature information.Whenn+1 attitude sensors are placed, attitude sensors are placed, the length of the transmission line in each interval li,i+1 is determinate. Meanwhile, the angles θi,i+1 between successive sensors can be calculated using θi,i+1=cos-1ti∙ti+1, 0≤i≤n. We cumulatively reform the rotation angles and arc-lengths as follows: (8)s0=0,s1=l0,1,⋯,sn=∑i=0n-1li,i+1(9)θ0=0,θ1=θ0,1,⋯,θn=∑i=0n-1θi,i+1.We obtain the data pointsθ=θ0,θ1,⋯,θn with the corresponding parameter s=s0,s1,⋯,sn, and the global curve of θ(s) is obtained by the classical B-spline interpolation, as shown in Figure 4.Figure 4
Diagram of the B-spline global interpolation curve (a)θ(s); (b) θt(s).
(a) (b)Then, we calculate each segment of the rotation angle information as(10)θi,i+1s=θs-θi,for0≤i≤n-1.Then, we revise (6) in the global frame, and the tangent curve ts is reformulated by(11)ts=∏i=0nexp⁡viθtst0,s∈0,sn(12)vi=ti×ti+1ti×ti+1.By combining (1) and (11) with (12), the desired transmission line galloping curve can be reconstructed.


### 3.3. Reconstruction Curve Correction

In the actual situation, the results will be affected by the measurements. The errors from the sensors will have a direct impact on the shape of the reconstruction curves. Moreover, the error increases at the end of the reconstruction curve because of the integral accumulated error. Because the location of the two transmission towers is determined, we correct the curve using the linear mapping(13)xc=plx-p⌢lxl×x+xyc=ply-p⌢lyl×x+yzc=plz-p⌢lzl×x+zwhere pl and p⌢l denote the transmission tower positions in reality and reconstruction, respectively; l is the span of the transmission tower; x is the span distance to the first transmission tower, 0≤x≤l; and x,y,z and xc,yc,zc are the three-dimensional coordinates of the reconstruction and correction, respectively. The endpoint of the reconstruction curve clearly matches the transmission tower after the correction. In addition, if the reconstruction of the suspension point is far from the real position, this situation is regarded as transmission line fault, which will be explained in the next section. The reconstruction curve after correction is shown in Figure 5.Figure 5
Reconstructions of the transmission lines: galloping curves withNs=8 sensors (red square), direct reconstruction curve (yellow dotted line), and correction curve (pink dash line).


## 3.1. Tangential Information with Arc-Length Constraint

Historically, quaternions were conceived by Hamilton as extensions to complex numbers. A quaternionq is defined by(2)q=w+xi+yj+zk=x,y,z,w=v,wwhere v=x,y,z∈R3, w∈R, and i, j, and k are different imaginary units with i2=j2=k2=ijk=-1. More detailed properties of quaternions can be found in [31].Given a unit quaternionq∈S3, i.e., Nq=qq∗=1, particularly if q=cos⁡θ+v0sin⁡θ∈S3 for an angle θ and an arbitrary vector t=x,y,z,0, the 3D rotation operation is defined as(3)Qt=qtq-1where Qt denotes the rotation by 2θ around the axis v0.For example, letv0 be a unit vector with inclination (or elevation) angle α and azimuth angle β. Then, q=sin⁡αcos⁡βsin⁡θ,sin⁡αsin⁡βsin⁡θ,cos⁡αsin⁡θ,cos⁡θ, and t2 is the vector after the rotation operation; i.e., t2=Qt0, as shown in Figure 2. Essentially, the rotation Qt=qtq-1 is divided into two steps. First, t1 is obtained from a θ rotation around the axis v0 from t0; i.e., t1=qt0. Second, t2 is obtained from a θ rotation around the axis v0 from t1; i.e., t2=t1q-1.Figure 2
Geometry of rotation for a given axisv0 and vector t0.Now, we return to the problem ofts interpolation. By exploring the attitude sensor, we can obtain the 3-dimensional angle information, which can be converted into tangential information. The tangent vectors are unit vectors that lie on the Gaussian sphere.A simple and possible answer is that these unknown tangential coefficients lie on the great arc between pointst0 and t1, as shown in Figure 3. This arc is named the geodesic curve and is in S2 instead of S3, which is given by(4)γθt=qθtt0,θt∈0,θ0,1(5)qθt=v0sin⁡θt,cos⁡θt=exp⁡v0θtwhere θt0,t1=cos-1t0∙t1∈0∘,180∘ is the angle between vectors t0 and t1. v0=t0×t1/t0×t1 is a unit vector in the direction of the cross-product of t0 and t1.Figure 3
Geodesic curve for given tangent vectorst0 and t1.According to the arc-length constraint,θt is a variable that is associated with the arc-length parameters s, i.e., θt(s), and it satisfies θt0=0 and θtlt0,t1=θt0,t1, where lt0,t1 is the arc-length between contiguous sensor nodes. Now, we obtain t(s) between the first two sensor nodes as follows:(6)ts=γθts=qθtst0,s∈0,lt0,t1.


## 3.2. Curvature Information Using B-Spline Interpolation

From (6), we obtain the curvature κs:(7)κs=dtsds=q′θtst0θt′s=q′θtst0·θt′s=θt′s,s∈0,l0,1.Sinceq′(θt(s)) and t0 are unit quaternions, the multiplication q′(θt(s))t0=1. Therefore, the curvature κs depends on θt′(s) or θt(s).A simple example isθt(s)=θ0,1/l0,1s, s∈0,l0,1. Apparently, κs=θt′(s)=θ0,1/l0,1 is a constant, and the resulting curve (1) is a circular arc, which is not suitable for the situation of the transmission galloping curve.Now, the calculation ofts is transformed into the problem of evaluating θt(s). Since the transmission line is an entire conductor, there is an internal connection between successive sensor nodes. Thus, a modified form of B-spline interpolation is introduced to estimate the curvature information.Whenn+1 attitude sensors are placed, attitude sensors are placed, the length of the transmission line in each interval li,i+1 is determinate. Meanwhile, the angles θi,i+1 between successive sensors can be calculated using θi,i+1=cos-1ti∙ti+1, 0≤i≤n. We cumulatively reform the rotation angles and arc-lengths as follows: (8)s0=0,s1=l0,1,⋯,sn=∑i=0n-1li,i+1(9)θ0=0,θ1=θ0,1,⋯,θn=∑i=0n-1θi,i+1.We obtain the data pointsθ=θ0,θ1,⋯,θn with the corresponding parameter s=s0,s1,⋯,sn, and the global curve of θ(s) is obtained by the classical B-spline interpolation, as shown in Figure 4.Figure 4
Diagram of the B-spline global interpolation curve (a)θ(s); (b) θt(s).
(a) (b)Then, we calculate each segment of the rotation angle information as(10)θi,i+1s=θs-θi,for0≤i≤n-1.Then, we revise (6) in the global frame, and the tangent curve ts is reformulated by(11)ts=∏i=0nexp⁡viθtst0,s∈0,sn(12)vi=ti×ti+1ti×ti+1.By combining (1) and (11) with (12), the desired transmission line galloping curve can be reconstructed.


## 3.3. Reconstruction Curve Correction

In the actual situation, the results will be affected by the measurements. The errors from the sensors will have a direct impact on the shape of the reconstruction curves. Moreover, the error increases at the end of the reconstruction curve because of the integral accumulated error. Because the location of the two transmission towers is determined, we correct the curve using the linear mapping(13)xc=plx-p⌢lxl×x+xyc=ply-p⌢lyl×x+yzc=plz-p⌢lzl×x+zwhere pl and p⌢l denote the transmission tower positions in reality and reconstruction, respectively; l is the span of the transmission tower; x is the span distance to the first transmission tower, 0≤x≤l; and x,y,z and xc,yc,zc are the three-dimensional coordinates of the reconstruction and correction, respectively. The endpoint of the reconstruction curve clearly matches the transmission tower after the correction. In addition, if the reconstruction of the suspension point is far from the real position, this situation is regarded as transmission line fault, which will be explained in the next section. The reconstruction curve after correction is shown in Figure 5.Figure 5
Reconstructions of the transmission lines: galloping curves withNs=8 sensors (red square), direct reconstruction curve (yellow dotted line), and correction curve (pink dash line).


## 4. Sensors Placement and Galloping Alert

Not only the number of attitude sensors but also the placement position will affect the reconstruction performance of galloping curves. To capture the galloping waves with the least number of sensors, the placement of the sensors should meet two conditions. (1) Ensure the angle variation between the two adjacent sensors is less than 180 degrees and (2) ensure this variation is monotonic.Easy to illustrate, we take the two-dimensional curve as an example (three-dimensional case is likewise). The tangent vectors are unit vectors that formed a unit circle, as shown in Figure6(a). The variation of the tangent vector can be interrupted as the rotation on the circle. In order to capture every waveform, the variation angle θi,i+1 is not more than a half circle; i.e., θi,i+1≤180∘, which is calculated as (14)θi,i+1=cos-1ti∙ti+1∈0∘,180∘.Figure 6
Illustration of the sensors installment. (a) The two-dimensional tangent vectors and (b) corresponding reconstruction curves.
(a) (b)Meanwhile, the variation should be monotonic; otherwise, partial curvature information will be lost. As shown in Figure6(b), if the sensor node is installed at the position t1′ instead of t1, then curvature information from t1 to t1′ is missing. Meanwhile, we lost a galloping waveform between adjacent sensor nodes.In summary, sensors should be ideally mounted on the inflection point of the curve. However, this condition is quite critical in the real application. We could not determine the position of the inflection point since transmission line galloping or change the position of sensors which had been placed on the transmission line. The other way to improve the reconstruction precision is to increase the number of sensors. For example, if we add another sensor node betweent1 and t1′, then, we can get a better curve fitting to the real one.The placement of sensors should match the galloping feature, i.e., the conduct sag [19] is easy ice-coating and galloping, and the place closer to the suspension ends is more fixed. Therefore, more sensors should be installed close to the conduct sag to capture more galloping information. In the study, a simple evenly installment method was proposed. We place the sensors evenly along the span direction, which can be mainly in line with the suspension characteristics of the transmission line. Moreover, as mentioned before, more sensors could get a better curve fitting to the real galloping curves. Figure 7 shows a case for sensor placements with different number of sensors. As we can see, Ns=6 sensors mounted on the inflection point can achieve a fine performance. However, the performance gets worse when Ns increased to 8 with the even placement. But then Ns rises to 10 or more; the reconstruction can achieve acceptable accuracy.Figure 7
The reconstruction curves with different numbers of sensors. (a)Ns=6 mounted on the inflection point. (b), (c), and (d) Evenly installed along the span with Ns=8, Ns=10, and Ns=12, respectively.
(a) (b) (c) (d)With a proper number of sensors, the full-scale transmission line curve could be well reconstructed in the monitoring center. We can watch the reconstructed curves directly and can extract the galloping parameters (amplitude, frequency, and waveform) from the reconstructed data. It will be beneficial for the state monitoring and safety improvement of power delivery systems. For the case that the power transmission tower is slanted or collapsed, or the transmission line is fractured, as shown in Figure8, the angle information will change accordingly, furthermore, the reconstruction curve will be greatly different (especially for the conductor sag and the suspension ends). These anomalous phenomena will be reflected directly in the reconstructed curves, so the workers can discover failure and repair the fault in the first time.Figure 8
Illustration of the galloping alert for (a) tower collapse and (b) line fracture.
(a) (b)


## 5. Simulation and Analysis

### 5.1. Simulation Setup

In this experiment, we used the single-degree-of-freedom (SDOF) model which only considers the vertical vibration to test the effectiveness of the method. The overhead transmission line curve with vertical vibration can be approximately formulated [3] as(15)yx,t=hlx-γxl-x2σ0cos⁡β+A0sin⁡nπlxsin⁡ωtwhere h (m) and l(m) are the height difference and span of the crossing tower, respectively. γ (MPa/m) is the specific load, which is defined as the load in the unit length and cross-sectional area. σ0 (MPa) is the horizontal axial stress; β=arctan⁡(h/l) is the angle of elevation difference. A0 (m) denotes the galloping amplitude, n is the galloping loops, and ω(rad) is the galloping frequency. Using the transmission line parameters of the Great Span Project in Hubei province, China, the corresponding parameters are as follows: l=1055 m, h=19.5 m, γ=42.6×10-3MPa/m, σ0=101.45MPa, A0=10 m, n=5, and ω=1.678 rad. In the simulation, we take the first five time slots of transmission line data as instances, as shown in Figure 9.Figure 9
Illustration of SDOF galloping model in the first five time slots.We set two evaluating indicators to assess the method performance.eL indicates the length error between real and reconstructed lengths of the transmission line and is expressed as(16)eL=Ls-LrLs×100%where Ls and Lr are the length of the transmission line and reconstruction length, respectively. eP evaluates the reconstruction performance and is given by (17)eP=1Ls1Ns∑i=1Nspi-p⌢i22×100%where Ns is the number of sensors and pi and p⌢i are the positions of the i-th attitude sensor in reality and reconstruction, respectively. To some extent, eP represents the shape difference between the actual curves and the reconstructions.


### 5.2. Computation Unitpl and Computational Complexity

The reconstruction performance is substantially dependent on the computation unit. In the simulation,Ns was fixed to 10 and the computation unit pl was chosen as pl=[1,5,10,15,20,25](m), and the reconstruction performance as shown in Figure 10. With pl increasing, eL rises sharply and eP also shows a trend of slow growth. Therefore, pl has a high impact on the length reconstruction. The bigger integrator unit leads to increased cumulative error. Meanwhile, the computation unit will directly decide the computational complexity. In this paper, all calculations were performed on a PC with a 3 GHz CPU (Intel Pentium G3220), and the corresponding average time consumption was t=[0.903,0.166,0.086,0.059,0.046,0.039](s). It is clear that, with the minor unit, the reconstruction will achieve higher precision but also higher computational complexity. It should be mentioned that, even with the minor unit pl=1 m, the average reconstruction time is less than one second, which meets the demand of the real-time application.Figure 10
Reconstruction performance versus computation unitpl.


### 5.3. Number of SensorsNs

Figure11 shows the reconstruction performance (error bar) versus different numbers of sensors. The computation unit pl was fixed to 5 m and Ns was varied from 7 to 12 at the step of 1. It shows that Ns almost does not affect the length reconstruction eL (as shown in blue line) and has a slow decline trend in the average position error eP. This result is reasonable because more sensors lead to higher accuracy. When the sensors are twice than the galloping loops; i.e., Ns=10, the average eP has a sharp fall from 0.24% to 0.19%. Also, it can be observed that low bound of eP with Ns=7 is even lower than the one with Ns=11. This situation is because that placement with Ns=7 may be closed to the inflection point for some galloping time slots.Figure 11
Reconstruction performance versus different numbers of sensorsNs.


### 5.4. Curve Correction

As Figure10 shows, when pl is small (with pl=1m), the correction operation has limited effect for the reason that the reconstruction accuracy is relatively high. When pl goes larger, correction operation can significantly reduce the length error eL and also has slight improvement on the position error eP. Moreover, as shown in Figure 11, it shows that correction operation can effectively adjust the length error induced by the computational unit pl=5m. Also it can successfully decrease eP when sensors are not enough installed. When the sensors are more than twice of the galloping loops, the effect of correction operation is reduced since the algorithm achieves an acceptable reconstruction accuracy.


## 5.1. Simulation Setup

In this experiment, we used the single-degree-of-freedom (SDOF) model which only considers the vertical vibration to test the effectiveness of the method. The overhead transmission line curve with vertical vibration can be approximately formulated [3] as(15)yx,t=hlx-γxl-x2σ0cos⁡β+A0sin⁡nπlxsin⁡ωtwhere h (m) and l(m) are the height difference and span of the crossing tower, respectively. γ (MPa/m) is the specific load, which is defined as the load in the unit length and cross-sectional area. σ0 (MPa) is the horizontal axial stress; β=arctan⁡(h/l) is the angle of elevation difference. A0 (m) denotes the galloping amplitude, n is the galloping loops, and ω(rad) is the galloping frequency. Using the transmission line parameters of the Great Span Project in Hubei province, China, the corresponding parameters are as follows: l=1055 m, h=19.5 m, γ=42.6×10-3MPa/m, σ0=101.45MPa, A0=10 m, n=5, and ω=1.678 rad. In the simulation, we take the first five time slots of transmission line data as instances, as shown in Figure 9.Figure 9
Illustration of SDOF galloping model in the first five time slots.We set two evaluating indicators to assess the method performance.eL indicates the length error between real and reconstructed lengths of the transmission line and is expressed as(16)eL=Ls-LrLs×100%where Ls and Lr are the length of the transmission line and reconstruction length, respectively. eP evaluates the reconstruction performance and is given by (17)eP=1Ls1Ns∑i=1Nspi-p⌢i22×100%where Ns is the number of sensors and pi and p⌢i are the positions of the i-th attitude sensor in reality and reconstruction, respectively. To some extent, eP represents the shape difference between the actual curves and the reconstructions.


## 5.2. Computation Unitpl and Computational Complexity

The reconstruction performance is substantially dependent on the computation unit. In the simulation,Ns was fixed to 10 and the computation unit pl was chosen as pl=[1,5,10,15,20,25](m), and the reconstruction performance as shown in Figure 10. With pl increasing, eL rises sharply and eP also shows a trend of slow growth. Therefore, pl has a high impact on the length reconstruction. The bigger integrator unit leads to increased cumulative error. Meanwhile, the computation unit will directly decide the computational complexity. In this paper, all calculations were performed on a PC with a 3 GHz CPU (Intel Pentium G3220), and the corresponding average time consumption was t=[0.903,0.166,0.086,0.059,0.046,0.039](s). It is clear that, with the minor unit, the reconstruction will achieve higher precision but also higher computational complexity. It should be mentioned that, even with the minor unit pl=1 m, the average reconstruction time is less than one second, which meets the demand of the real-time application.Figure 10
Reconstruction performance versus computation unitpl.


## 5.3. Number of SensorsNs

Figure11 shows the reconstruction performance (error bar) versus different numbers of sensors. The computation unit pl was fixed to 5 m and Ns was varied from 7 to 12 at the step of 1. It shows that Ns almost does not affect the length reconstruction eL (as shown in blue line) and has a slow decline trend in the average position error eP. This result is reasonable because more sensors lead to higher accuracy. When the sensors are twice than the galloping loops; i.e., Ns=10, the average eP has a sharp fall from 0.24% to 0.19%. Also, it can be observed that low bound of eP with Ns=7 is even lower than the one with Ns=11. This situation is because that placement with Ns=7 may be closed to the inflection point for some galloping time slots.Figure 11
Reconstruction performance versus different numbers of sensorsNs.


## 5.4. Curve Correction

As Figure10 shows, when pl is small (with pl=1m), the correction operation has limited effect for the reason that the reconstruction accuracy is relatively high. When pl goes larger, correction operation can significantly reduce the length error eL and also has slight improvement on the position error eP. Moreover, as shown in Figure 11, it shows that correction operation can effectively adjust the length error induced by the computational unit pl=5m. Also it can successfully decrease eP when sensors are not enough installed. When the sensors are more than twice of the galloping loops, the effect of correction operation is reduced since the algorithm achieves an acceptable reconstruction accuracy.


## 6. Experimental Test

Since this is an original work and limited by the experimental facility, a primary test platform in Figure12(a) was built to verify the effectiveness of the proposed method. The span of the two hanging poles (simulated as transmission towers) is 9.2 m, and the length of the mooring rope (simulated as the transmission line) is 9.5 m. Eight sensor nodes (shown as red blocks) were evenly installed, in span direction, to capture the dip information (pitch angle and yaw angle). Since we ignore the twist of the test line, Roll angle information is not used in the experiment. The direction of the two hanging poles was fixed to they-axis, and the horizontal and vertical directions were set up as thex-axis andz-axis, respectively.Figure 12
Test platform to simulate transmission line galloping: (a) experimental platform; (b) structure of the sensor node.
(a) (b)As shown in Figure12(b), a sensor node was designed that consisted of an attitude sensor (MPU 9250), a wireless communication device (RF2401), microcontroller (STM32F072RBT6), and a power supply module. The attitude sensor was used to obtain the dip information. The microcontroller served as an interface between the attitude sensor and the wireless communication device and performed network-related tasks such as clock signal buffering and data retransmission.Four states were set up to simulate the galloping curves: initial rest state, vertical vibration, horizontal swing, and tortuosity. In order to obtain the spatial position information of the sensor nodes, all the states were tested under static status. Meanwhile, two precision grids were set topl=0.01 m and pl=0.1 m, and the reconstruction points Ns=Ls/pl= 950 and 95 were used for reconstruction. The reconstruction performance was shown in Figure 13 (with pl=0.01 m) and Table 1.Table 1
Reconstruction performance of four experimental states.
Time(s) 3.02 (pl=0.01 m) 0.49 (pl=0.1 m) Correction No Yes No Yes Error (%) e L e P e L e P e L e P e L e P (a) 0.316 0.714 0.349 0.331 1.276 0.627 0.875 0.463 (b) 0.316 1.364 0.511 0.439 1.274 0.472 1.595 0.565 (c) 0.316 0.491 0.198 0.502 1.276 0.770 0.696 0.601 (d) 0.316 0.615 0.489 0.304 1.274 0.480 0.787 0.431 Mean 0.316 0.796 0.386 0.394 1.275 0.587 0.988 0.515Figure 13
Experimental results of four different states.Ns=8 sensors (red square) were explored on the test line. The reconstruction curves before and after the correction are denoted by the solid blue line and the dotted yellow line, respectively: (a) initial rest state; (b) vertical vibration; (c) horizontal swing; (d) tortuosity.
(a) (b) (c) (d)Because of the limitation of the experimental instrument, we observe that the real experimental result is worse than the simulated one. Nevertheless, the experimental results show that the reconstruction curves are reasonably consistent with the actual shape of the test lines. Comparing the direct reconstruction results and the corrected one, we observe that the correction operation increases the length erroreL slightly but decreases the position error eP markedly. To some extent, eP represents the shape difference between the original curves and the reconstructions. Thus, the reconstruction is closer to the original curve with smaller eP, and the correction operation enhances the reconstruction curves by making them much closer to the actual one.A comparison of the results with different precision grids clearly shows thatpl directly affects the reconstruction of the arc-length: eL=0.316% with pl=0.01m and eL=1.275% with pl=0.1m. However, pl only slightly affects the reconstruction position error: eP=0.394% with pl=0.01m and eP=0.515% with pl=0.1m. The time consumption sharply decreases from 3.02 s to 0.49 s with a small increase in eP. In other words, the proposed method costs less than 0.13s (with pl=0.1 m) for one reconstruction, which is suitable for real-time monitoring of galloping systems.


## 7. Conclusions

This paper presents a new method for the full-scale reconstruction of transmission line galloping curves with few attitude sensors; the method uses only a set of tangential information with arc-length constraints instead of the spatial positions. With a small amount of dip information data from the attitude sensors, the galloping curves are reconstructed in the monitoring center, from which the transmission line galloping status during the monitoring time can be directly observed. Meanwhile, a low-complexity algorithm is introduced for the real-time reconstruction. We used the quaternion rotation to fulfill the tangent vectors between adjacent sensors, and the B-spline interpolation was applied to obtain the global curvature information. Moreover, a linear mapping was used to correct the reconstruction curves. In addition, we discussed the sensor installment and proposed a placement that sensors, with twice more than the galloping loops, were fixed evenly along the span direction. SDOF galloping model and an experimental platform were set up to verify the accuracy and applicability of the method. Both simulation and experimental results demonstrate that the proposed method can accurately reconstruct the galloping curves. Further, the average time cost of the proposed method is less than one second, which is suitable for real-time monitoring demands. This method can provide full-scale galloping trace information for power transmission line monitoring, which will be very helpful for the state monitoring and safety improvement of power delivery systems.



---

*Source: 1095842-2018-06-27.xml*