Failure Patterns of Transmission Tower-Line System Caused by Landslide Events

Abstract

A transmission tower-line system consists of many towers and cables, which is a long-span structure. Due to topographical restrictions, the structure inevitably passes through wide areas prone to landslide hazards. The landslide causes the failure of the tower foundation and threatens the safety of the power transmission system. Therefore, the aim of this study was to investigate the collapse patterns of a transmission tower-line structure under landslides. The explicit analysis method was used to deal with the nonlinear dynamic response equations of the structure. To prove the effect of the landslides on the collapse patterns of the transmission tower-line structure, two foundation failure cases (Case 1 and Case 2) under landslides were considered. For each case, the displacement responses at tower tops and the progressive collapse behavior were analyzed, which illustrate that the failure patterns of the tower-line structure under landslides are dominated by two factors, i.e., the failure mode of the tower foundation and the pulling effect of cables. The overall tilting and hinge formation in the middle of the tower body are the main collapse patterns of the structure under Case 1 and Case 2, respectively. In addition, because of the pulling of cables, the collapse of one tower induced by the landslides always triggers the local failure of neighboring tower heads.

1. Introduction

Due to the increasing demand for electricity, the malfunctions of the power grid can produce adverse effects on people’s livelihoods and the national economy. For example, the massive outage in the state of Texas in the United States [1] has caused a more than $200 billion economic loss [2] and 70 deaths [3], and the severe electricity shortage in China has caused financial and operational crises for many companies [4]. The tower-line structure is a carrier for the transmission of electric energy, which belongs to the main part of the power grids. For such a long-span structure, the towers are inevitably constructed on high-altitude mountainous areas, where the geohazard, especially the landslide, usually occurs at the steep slopes and ridges of the mountains [5]. Landslides can probably lead to the tilting and even collapse of transmission towers, as shown in Figure 1. A recent report [6] by China Electric Power Research indicates that there are more than 1500 transmission towers that are in the potential threat of landslides. Figure 2 illustrates the field pictures of the landslide located at Tower N18 in the Dan-Gong transmission line in China. Therefore, the investigation into the failure mechanisms of transmission tower-line structures under landslides is significant.

Landslide events can be induced by earthquakes and rainfall and they are a major geomorphological process that influences landscape evolution in mountainous areas [7]. The mechanism of the development of the landslide has been studied by many researchers. Some studies [8,9,10,11] have investigated the effect of various influential factors, such as the slope angle, slope height, geological structure, and rock type, on the landslide spatial distribution. In addition, the landslide travel distance has also been evaluated by the theoretical model [12], empirical model [13,14], and numerical model [15]. Because the essential input parameters in the theoretical and numerical models are difficult to be measured during the landslide movement [16], the two models need to be greatly simplified through many assumptions. The empirical model follows the rules that are calibrated based on the experimental and field-measured data. It can describe the landslide using simple mathematical formulas and therefore it is more popular than the other two above-mentioned models.

The transmission tower-line system consists of many towers, which are connected by a series of cables. For this type of structure, the progressive collapse behavior usually occurs [17,18]. Therefore, some studies have been conducted to investigate the nonlinear dynamic response of transmission tower-line systems under extreme environmental conditions [19,20,21,22,23]. Zheng and Fan [19] investigated the progressive collapse mechanisms of a tower-line system subjected to hurricanes through fully coupled dynamic simulations, showing that the collapsed transmission tower would lead to the cascade collapse of the neighboring towers. Kassem et al. [20] also conducted a comparative study on the failure patterns of a transmission tower under far-field and near-field earthquakes. Li et al. [21] investigated the collapse mechanism and failure modes of the transmission tower under strong wind load, showing that the wind-induced coefficients in the current design code are not accurate enough. Fu et al. [22] comprehensively evaluated the failure modes of a transmission tower under combined wind and rain loads and they provided a critical collapse surface to quickly and conveniently judge the structural failure states at different wind intensities. However, the above-mentioned studies are mainly on the failure patterns of tower-line structures under hurricanes, earthquakes, and so on, whereas the relevant research for the structures subjected to landslides is scarce. Though one study [24] has tried to identify the failure modes of transmission tower foundations under landslides, the consequent collapse patterns of the tower-line structures have not been investigated yet.

In this research, we selected a segment of the transmission tower-line system over the Yangtze River as the study object. The nonlinear dynamic response and the collapse process of the structure under landslides were simulated. To deal with the discontinuous problems involved in the collapse simulations, the MATLAB code based on the explicit computation method was adopted. Two failure cases under landslides were considered. For each case, the structural displacement response and the collapse scenarios at the specified moments were analyzed and the corresponding collapse mechanisms of the transmission tower-line system under landslides were then identified.

2. The Yanzi Landslide

As shown in Figure 3, the Yanzi landslide is located in the power transmission zone (Tower 2) from Wuxia to Xilingxia. The overall terrain is low in the east and north and high in the west and south, where the volume of the landslide body is about 60 × 10⁴ m³ with an area of 60 × 10⁴ m² and an average thickness of 10 m [24]. The landslide body looks like a long tongue, which has a slope of about 25–40 degrees. The material composition of the landslide is characterized by the silty clay mixed with gravel soil. The slipping zone is the contact surface between the overlying deposits and the bedrock, as shown in Figure 4. Since 2002, the Yanzi landslide has shown many creep-related slipping events due to rainfall and many cracks have presented on the surface of the slope. Then, the G209 National Highway was constructed through the slope and a 100 mm crack formed between the transmission tower feet. More information about the Yanzi landslide can be found in Ref. [24].

As shown in Figure 3, the prototype of the transmission tower-line system is a segment of a 500 kV transmission line, which involves three lattice transmission towers with the same structural design and four spans of cables. The tower height is 122 m and the foot-to-foot distance is 20.24 m, as shown in Figure 5, which illustrates the top and front views of the structure. The main legs and diagonal members are steel tubes and the dimensions are also given in Figure 5. The material of the tower members is Q235 steel and the material constants are listed in Table 1. The cables can be classified into conductors and ground wires, whose span length is 800 m. The type of conductors is JL/G1A-400/35 with a diameter of 26.82 mm, a mass per unit length of 1.35 kg/m, and an elastic modulus of 65 GPa. The type of ground wires is JLB40-150 and the corresponding values are 15.75 mm, 0.7 kg/m, and 103.6 GPa, respectively.

3. Wind Load Formation

When the wind travels through the transmission tower-line structure, drag forces on the transmission towers are produced. For convenience in calculating the wind loads, the tower was divided into 12 panels, as illustrated in Figure 5. The drag forces, $F_x$ (X-direction) and $F_y$ (Y-direction), can be expressed as

(1)$F_x=0.5\rho U^2{\cos }^2\alpha C_{dx}{A}_x$

(2)$F_y=0.5\rho U^2{\sin }^2\alpha C_{dy}{A}_y$

where $U$ is the wind speed; $\alpha$ is the angle between wind and panels; $C_{dx}$ and $C_{dy}$ are the drag coefficients of tower members [25]; $A_x$ and $A_y$ are the projected areas of tower members; $\rho$ is the air density.

The drag forces on cables, $F_{xl}$, can be calculated using

(3)$F_{xl}=0.5\rho U^2{\cos }^2\beta C_{dl}{A}_{xl}$

where $\beta$ is the angle between wind speed and cables; $C_{dl}$ is the drag coefficient of cables, whose value is 1.0 [19]; $A_{xl}$ is the projected area of cables.

The wind speed $U$ is the summation of the turbulent component $u$ and the mean component ${\overline{u}}_{\mathrm{m}}$.

(4)$U(x,y,z,t)=u(x,y,z,t)+{\overline{u}}_{\mathrm{m}}(z)$

The mean wind component ${\overline{u}}_{\mathrm{m}}$ varies with the height $z$ and it can be expressed using a normal wind profile.

(5)${\overline{u}}_{\mathrm{m}}(z)=u_{\mathrm{ref}}{(\frac{z}{z_{\mathrm{ref}}})}^{\eta }$

where $u_{\mathrm{ref}}$ is the wind speed at the reference height $z_{\mathrm{ref}}$; $\eta$ is the ground roughness index, taken as 0.14 [26].

For a wind field point, the turbulent wind speed $u$ can be simulated using the Kaimal spectrum [27]. However, to consider the turbulent wind speed correlations between Point i and Point j in a wind field, the cross-spectrum, $S_{ij}(f)$, (Equation (6)) should be used, which can be obtained based on the Kaimal spectrum and correlation function $Co{h}_{ij}$ (Equation (7)).

(6)$S_{ij}(f)=\sqrt{S_{ii}(f)S_{jj}(f)}Co{h}_{ij}(f),$

(7)$Co{h}_{ij}(f)=\exp \left\{-12{\left[{(fr/u_{\mathrm{ref}})}^2+{(0.12r/L_{\mathrm{c}})}^2\right]}^{0.5}\right\}$

where $f$ is the wind speed frequency; $S_{ii}(f)$ is the Kaimal spectrum at Point i; $S_{jj}(f)$ is the Kaimal spectrum at Point j; $r$ is the distance between the two points; $L_{\mathrm{c}}$ is the correlation coefficient between the two points [27]. If there are n field points, the cross-spectrum matrix, $S(f)$, can be yielded.

(8)$S(f)={\left[S_{ij}(f)\right]}_{n\times n},$

Then, the turbulent wind speed time history at Point i, which considers the influences from other field points, can be written as

(9)$u_i(t)={\displaystyle \sum _{j=1}^i{\displaystyle \sum _{k=1}^{N_f}\left|H_{ij}(2\mathsf{\pi }{f}_k)\right|}}\sqrt{4\mathsf{\pi }\Delta f}\times \cos \left[2\mathsf{\pi }{f}_{k}t-{\theta }_{ij}(2\mathsf{\pi }{f}_k)+{\phi }_{jk}\right]$

where $N_f$ is the number of sampling points; ${\phi }_{jk}$ is the phase angle that satisfies the uniform distribution; $H_{ij}$ is the lower triangular matrix that can be computed via the Cholesky decomposition of $S(f)$ [28]; ${\theta }_{ij}(2\mathsf{\pi }{f}_k){=\tan }^{-1}\left[\frac{\mathrm{Im}{H}_{ij}(2\mathsf{\pi }{f}_k)}{\mathrm{Re}{H}_{ij}(2\mathsf{\pi }{f}_k)}\right]$.

4. Dynamic Response Analysis Method of the Tower-Line Structure

Figure 6 provides the discretization model of the transmission tower-line structure, which is represented by a number of mass particles and member elements. The total mass of the structure is modeled by mass particles. The member element describes the topological relationship between mass particles. The element model developed by the researchers in [29] was adopted to simulate tower members, which can account for the yielding in tension, the buckling in compression, and the critical force degradation due to the Bauschinger effect [30]. For the cable members, the $\epsilon -\sigma$ relationship is shown in Figure 7. It can be found that when the cables are tensioned, the material behaves in a bilinear way. The slope of stage OA is $E_1$, and the slope of the yielding stage AB is $E_{\mathrm{Y}}$. The value of $E_Y/E_1$ is 0.05 in this study [31]. The cable member element in Figure 7 assumes that the cables cannot experience any compression actions, which means that the internal forces become zero when the cables are compressed.

The dynamic response of structures is actually a function of time and space. The motion response of mass particles can be regarded as a series of sequential motions in small time steps [32]. The motion of mass particles in a small time step is called a path unit, which can be calculated based on Newton’s second law [33]. For the mass particle i, the expression can be expressed as

(10)$\left[\begin{array}{ccc}{m}_i& 0& 0\\ 0& m_i& 0\\ 0& 0& m_i\end{array}\right]\left[\begin{array}{c}{\ddot{d}}_x\\ {\ddot{d}}_y\\ {\ddot{d}}_z\end{array}\right]=\left[\begin{array}{c}{F}_{x,\mathrm{int}}\\ F_{y,\mathrm{int}}\\ F_{z,\mathrm{int}}\end{array}\right]+\left[\begin{array}{c}{F}_{x,\mathrm{ext}}\\ F_{y,\mathrm{ext}}\\ F_{z,\mathrm{ext}}\end{array}\right]+\left[\begin{array}{c}{F}_{x,\mathrm{damp}}\\ F_{x,\mathrm{damp}}\\ F_{x,\mathrm{damp}}\end{array}\right]$

where $m_i$ is the mass of the particle i, which is the half sum of member elements connected to the particle; $d={[\begin{array}{ccc}{d}_x& d_y& d_z\end{array}]}^{\mathrm{T}}$ and $\ddot{d}={[\begin{array}{ccc}{\ddot{d}}_x& {\ddot{d}}_y& {\ddot{d}}_z\end{array}]}^{\mathrm{T}}$ are the vectors of displacement and acceleration of the particle i; ${\left[\begin{array}{ccc}{F}_{x,\mathrm{int}}& F_{y,\mathrm{int}}& F_{z,\mathrm{int}}\end{array}\right]}^{\mathrm{T}}$ is the vector of internal force in member elements, which can be obtained through the model of member elements; ${\left[\begin{array}{ccc}{F}_{x,\mathrm{ext}}& F_{y,\mathrm{ext}}& F_{z,\mathrm{ext}}\end{array}\right]}^{\mathrm{T}}$ is the vector of external force exerted on the particle i, which depends on the wind speed (Equations (1)–(3)) and landslide conditions used in this study. Finally, ${\left[\begin{array}{ccc}{F}_{x,\mathrm{damp}}& F_{y,\mathrm{damp}}& F_{z,\mathrm{damp}}\end{array}\right]}^{\mathrm{T}}$ is the vector of damping force, which can be further written as

(11)$\left[\begin{array}{c}{F}_{x,\mathrm{damp}}\\ F_{x,\mathrm{damp}}\\ F_{x,\mathrm{damp}}\end{array}\right]=-\mu \left[\begin{array}{ccc}{m}_i& 0& 0\\ 0& m_i& 0\\ 0& 0& m_i\end{array}\right]\left[\begin{array}{c}{\ddot{d}}_x\\ {\ddot{d}}_y\\ {\ddot{d}}_z\end{array}\right]$

where $\mu$ is the damping coefficient, which is the same as the definition in Ref. [34]; $\dot{d}={[\begin{array}{ccc}{\dot{d}}_x& {\dot{d}}_y& {\dot{d}}_z\end{array}]}^{\mathrm{T}}$ is the velocity vector of the particle i.

To avoid the complex iteration calculations and the formation of a global stiffness matrix, Equation (10) can be solved by an explicit central difference method [35]. For this, the acceleration and velocity can be written as

(12)$\ddot{d}=\frac{1}{\Delta t^2}\left({\ddot{d}}_{n+1}-2{\ddot{d}}_n+{\ddot{d}}_{n-1}\right)$

(13)$\dot{d}=\frac{1}{2\Delta t}\left(d_{n+1}-d_{n-1}\right)$

where $d_{n+1}$, $d_n$, and $d_{n-1}$ are displacement vectors of a particle at steps n + 1, n, and n-1, respectively; $\Delta t$ is the time step, taken as 10⁻⁵ s [36]. By substituting Equations (11)–(13) into Equation (10), the displacement response of the structure can be obtained.

Because the time step in a path unit is pretty small, the deformation and internal force increment of member elements in the path unit are linear and small [37]. The material properties of elements, the number and mass of particles, and the topological relationship can be regarded to be not changed in the path unit, but they can be changed between path units. Therefore, the total deformation of member elements could be nonlinear and very large throughout the whole analysis process, and even the discontinuous issues of member fractures can also be dealt with between path units. The above-mentioned calculation method has been programmed using MATLAB [38], which was already validated by comparison with other numerical models and was found to be accurate with an error within 4.04% [33].

5. Collapse Analysis of the Tower-Line System under Wind and Landslide Action

5.1. Numerical Model of the Tower-Line Structure

As illustrated in Figure 3, the investigated segment includes three transmission towers connected by four spans of cables. MATLAB was used to create the numerical model of the transmission tower-line structure, as shown in Figure 8. The members in the numerical model are connected by hinges. The feet of Towers 1 and 3 (P1–P4 and P9–P12) are assumed to be fixed. However, the boundary condition at the feet of Tower 2 (P5–P8) varies with the landslide cases defined in Section 5.2. The boundary conditions at P13–P16 and P17–P20 are modeled using end springs and their stiffness is defined by the elastic characteristic of the transmission towers [39]. According to previous studies [40,41], the damping of the steel tower is 2% and the damping of cables is 20%. For each simulation, the wind loads are always applied to the structure, whereas the landslide actions are exerted on the structure over 5–15 s.

5.2. Load Cases

The landslide would lead to the foundation failure of transmission towers. In Figure 8, the P5-P8 are feet of Tower 2. We investigated two failure cases of the tower foundation on the landslide, as illustrated in Figure 9. For Case 1, P7 and P8 of Tower 2 are located in the sliding zone (Figure 9a); for Case 2, P6, P7, and P8 of Tower 2 are located in the sliding zone (Figure 9b). We assumed that part of the tower foundation is on a potential landslide. Once the slide occurs, the part of the tower foundation that sits on the potential landslide could lose its bearing capacity and the corresponding tower feet could lose their suddenly supports [24].

5.3. Progressive Collapse under Case 1

For Case 1, the displacement time histories at the tops of Towers 2 and 3 are shown in Figure 10. The collapse scenarios for the transmission tower-line system at the special moments are presented in Figure 11, Figure 12, Figure 13 and Figure 14. The displacement response of Tower 1 is not presented here since the response curves show similar features to those of Tower 3. The coordinate system of OXYZ is defined in Figure 8. In Figure 10, it can be found that the displacements in the X-, Y-, and Z-direction are very small at 5 s because there are only wind loads and dead weight applied to the structure for 1–5 s. After that, the feet of Tower 2 at P7 and P8 start losing their supports due to the landslide (Figure 9a) and the tower body gradually tilts in the X-direction, as shown in Figure 11. Accordingly, the displacement at the tower top also sees an increment in the X- and Z-direction, as illustrated in Figure 10a. It also can be observed in Figure 11 that the cables between transmission towers are also lifted by wind loads, which produce insignificant influences on the response of towers. As the tilt level of the tower body increases, the slack cables in the initial state (Figure 8) have straightened at t = 12 s, as shown in Figure 12. Due to the pulling effect by cables, the neighboring Towers 1 and 3 show the bending shape in the Y-direction at t = 14 s (Figure 13). The bending shape at the head of Towers 1 and 3 becomes more obvious at t = 14.3 s (Figure 14). This conclusion can also be drawn from the increment of the Y-displacement in Figure 10b. Then, the displacement at the top of Tower 3 shows a sharp decline in three directions and the duration from the start of the bending to the end of the total failure (Figure 10b) at the tower head is less than 1 s. This indicates that the collapse process of the neighboring transmission towers would be greatly accelerated by the pulling of the cables.

5.4. Progressive Collapse Simulation under Case 2

For Case 2, the displacement time histories at the tops of Towers 2 and 3 are shown in Figure 15. The collapse scenarios for the transmission tower-line system at the special moments are presented in Figure 16, Figure 17 and Figure 18. Similarly, the displacement response in the first 5 s is not obvious and therefore the curves are not plotted in Figure 15. Case 2 (Figure 9b) has one less support, namely P6, than Case 1 (Figure 9a), which means the feet of Tower 2 at P6, P7, and P8 would lose their support due to the landslide under Case 2. As illustrated in Figure 15a, the displacement response at the top of Tower 2 not only shows variations in the X- and Z-direction similar to these in Case 1 (Figure 10a), but also shows obvious changes in the Y-direction. At t = 9 s, it can be observed that the collapse pattern under Case 2 (Figure 16) is greatly different from that in Case 1, which presents great deformation along the tower body along with the tilting behavior in Case 1 (Figure 11). This can be explained by the fact that the cable tensions on both sides of Tower 2 are no longer equal when the displacement in the Y-direction occurs. As unbalanced tensions become larger and larger with the increment of the Y-displacement, a hinge forms in the middle of the tower at t = 13 s, as shown in Figure 17. After that, the head of Tower 2 drops greatly due to its own gravity and the displacement response curve in the Z-direction also sees an obvious decline (Figure 15a). In addition, as illustrated in Figure 18, a local failure also occurs in the upper part of the neighboring Tower 3, which is also due to the pulling of cables.

6. Discussion

The explicit central difference method was used to solve the dynamic response equation of the structure. It does not require the formation of a global stiffness and therefore it is convenient to process the strong geometry and material nonlinearities involved in the collapse of the tower-line structure caused by landslides. It should be noted that the explicit central difference method is conditionally convergent and the time step in Equations (12) and (13) should be carefully checked. For the finite element of the tower-line system, the time step should be smaller than the critical time step, $\Delta t_{\mathrm{c}}=T_{\mathrm{n}}/\pi$ ($T_{\mathrm{n}}$ is the minimum natural vibration period of tower members). Based on the method, the collapse simulations of the tower-line system subjected to landslides have been performed to illustrate that the collapse patterns are greatly affected by the failure modes of the tower foundation and the pulling of cables. On one hand, the tower foundation failure modes determine the number and the location of tower feet that lose the supports under landslides, which affects the azimuthal angle as the tower body tilts under landslides. On the other hand, the contribution of the cable pulling to the structural collapse varies with the azimuthal angles. For Case 1, the tilting azimuthal angle of Tower 2 with respect to the longitudinal direction is 90 degrees and the pulling of cables at both sides of Tower 2 counteracts; for case 2, the corresponding azimuthal angle is around 135 degrees, for which the pulling of cables between Towers 2 and Tower 3 is greatly different from that between Towers 1 and Tower 2. Therefore, the collapse scenarios under Case 1 are distinguished from those under Case 2. However, the above simulations do not consider the interactions between the tower foundation and soil during landslides, which possibly affect the failure patterns of the tower-line system. Future work will include the mentioned interactions to improve the accuracy of numerical results.

7. Conclusions

In this study, we analyzed the displacement response at tower tops and the structural collapse scenarios at special moments to investigate the failure mechanism of the tower-line system under landslides. The following conclusions were obtained:

The collapse failure patterns of the transmission tower-line structure vary with the tower foundation failure modes induced by landslides. For Case 1, two tower feet lose the supports under the landslide, and the overall tilting is the main collapse pattern for the tower on the landslide. For Case 2, three tower feet lose the supports under the landslide. The displacement response in the longitudinal direction is no longer zero, which leads to the different cable tensions on both sides of the tower. A hinge forming in the middle of the tower is the main collapse pattern in this case, which is very different from the collapse pattern in Case 1.

Considering that the tower-line system includes many towers that are connected by cables, the failure of one tower would cause the cascade collapse of the neighboring towers in the transmission line due to the pulling effect of cables. For Cases 1 and 2, the neighboring towers share similar collapse patterns, which are represented by the local failure close to the tower head.

Author Contributions

Conceptualization, H.Y. and H.-D.Z.; methodology, H.L. and H.-D.Z.; software, Z.-Q.Z., G.-F.Z. and H.-D.Z.; validation, H.Y., H.L., and D.-H.W.; formal analysis, H.Y.; investigation, H.L.; resources, D.-H.W.; data curation, Z.-Q.Z.; writing—original draft preparation, H.-D.Z.; writing—review and editing, H.Y.; visualization, G.-F.Z.; supervision, D.-H.W. and H.-D.Z.; funding acquisition, D.-H.W. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. The grid accidents induced by landslides: (a) tilting; (b) collapse.

Figure 2. The landslide located at tower N18 in the Dan-Gong transmission line: (a) general view of the landslide, with the yellow rectangle locating the position of the tower; (b) landslide at the foundation of the tower; (c) landslide at one leg of the tower (photographed by authors).

Figure 3. River-crossing transmission tower-line system (Google Maps, 17 August 2022).

Figure 4. Geological section of Yanzi landslide.

Figure 5. Top and front views of Tower 2.

Figure 6. The discretization model of the tower-line system.

Figure 7. The [Forumla omitted. See PDF.] relationship of cables in this study (OA: Elastic stage; AB: Yielding stage).

Figure 8. Numerical model of the structure.

Figure 9. Two failure cases of towers on the landslide. (a) Case 1; (b) Case 2.

Figure 10. Displacement response at tower tops under Case 1. (a) Tower 2; (b) Tower 3.

Figure 11. Behavior of the transmission tower-line structure at t = 9 s under Case 1.

Figure 12. Behavior of the transmission tower-line structure at t = 12 s under Case 1.

Figure 13. Behavior of the transmission tower-line structure at t = 14 s under Case 1.

Figure 14. Behavior of the transmission tower-line structure at t = 14.3 s under Case 1.

Figure 15. Displacement response at tower tops under Case 2. (a) Tower 2; (b) Tower 3.

Figure 16. Behavior of the transmission tower-line structure at t = 9 s under Case 2.

Figure 17. Behavior of the transmission tower-line structure at t = 13 s under Case 2.

Figure 18. Behavior of the transmission tower-line structure at t = 13.8 s under Case 2.

Table 1

Material constants of tower members.

Parameter	Value
Yield stress (MPa)	235
Elastic modulus (GPa)	206
Mass density (kg/m3)	7800
Poisson’s ratio	0.3

References

1. Wu, D.; Zheng, X.; Xu, Y.; Olsen, D.; Xia, B.; Singh, C.; Xie, L. An open-source extendable model and corrective measure assessment of the 2021 texas power outage. Adv. Appl. Energy; 2021; 4, 100056. [DOI: https://dx.doi.org/10.1016/j.adapen.2021.100056]

2. Ivanova, I. Texas Winter Storm Costs Could Top $200 Billion—More Than Hurricanes Harvey and Ike. 2021; Available online: https://www.cbsnews.com/news/texas-winter-storm-uri-costs (accessed on 1 August 2022).

3. Tengle, J.; Renault, M. Suspected Hypothermia Deaths in Homes Mount in Texas. 2021; Available online: https://apnews.com/article/houston-hypothermia-weather-conroe-texas-8323ab5f5c1612e632f7f2e6c2c20358 (accessed on 1 August 2022).

4. Chang, J. What Happened behind the Nationwide Power Cuts?. 2021; Available online: https://m.bjnews.com.cn/detail/163266456914638.html (accessed on 1 August 2022).

5. Liu, S.; Yin, K.; Zhou, C.; Gui, L.; Liang, X.; Lin, W.; Zhao, B. Susceptibility Assessment for Landslide Initiated along Power Transmission Lines. Remote Sens.; 2021; 13, 5068. [DOI: https://dx.doi.org/10.3390/rs13245068]

6. State Grid Equipment Department and China Electric Power Research Institute. Investigation Report on Geological Hazards of Power Transmission Lines (Yearly Report), 2017–2021; State Grid Equipment Department and China Electric Power Research Institute: Beijing, China, 2021.

7. Roering, J.J.; Kirchner, J.W.; Dietrich, W.E. Characterizing structural and lithologic controls on deep-seated landsliding: Implications for topographic relief and landscape evolution in the Oregon Coast Range, USA. Geol. Soc. Am. Bull.; 2005; 117, pp. 654-668. [DOI: https://dx.doi.org/10.1130/B25567.1]

8. Xu, C.; Dai, F.C.; Yao, X.; Chen, J.; Tu, X.B.; Sun, Y.; Wang, Z.Y. GIS-based landslide susceptibility assessment using analytical hierarchy process in Wenchuan earthquake region. Chin. J. Rock Mech. Eng.; 2009; 28, pp. 3978-3985.

9. Chigira, M.; Wu, X.; Inokuchi, T.; Wang, G. Landslides induced by the 2008 Wenchuan earthquake, Sichuan, China. Geomorphology; 2010; 118, pp. 225-238. [DOI: https://dx.doi.org/10.1016/j.geomorph.2010.01.003]

10. Qi, S.; Xu, Q.; Lan, H.; Zhang, B.; Liu, J. Spatial distribution analysis of landslides triggered by 2008.5. 12 Wenchuan Earthquake, China. Eng. Geol.; 2010; 116, pp. 95-108. [DOI: https://dx.doi.org/10.1016/j.enggeo.2010.07.011]

11. Guo, D.; Hamada, M. Qualitative and quantitative analysis on landslide influential factors during Wenchuan earthquake: A case study in Wenchuan County. Eng. Geol.; 2013; 152, pp. 202-209. [DOI: https://dx.doi.org/10.1016/j.enggeo.2012.10.012]

12. TAKAHASHI, T. Estimation of potential debris flows and their hazardous zones: Soft countermeasures for a disaster. Nat. Disaster Sci.; 1981; 3, pp. 57-89.

13. Rickenmann, D. Empirical relationships for debris flows. Nat. Hazards; 1999; 19, pp. 47-77. [DOI: https://dx.doi.org/10.1023/A:1008064220727]

14. Rickenmann, D. Runout prediction methods. Debris-Flow Hazards and Related Phenomena; Springer: Berlin/Heidelberg, Germany, 2005; pp. 305-324.

15. McDougall, S.D.; Hungr, O. Objectives for the development of an integrated three-dimensional continuum model for the analysis of landslide runout. Debris-Flow Hazards Mitigation: Mechanics, Prediction, and Assessment, Proceedings of the 3rd International DFHM Conference, Davos, Switzerland, Millpress 10–12 September 2003; Rickenmann, D.; Chen, C.L. Ios Press Inc.: Amsterdam, The Netherlands, 2003; pp. 481-490.

16. Guo, D.; Hamada, M.; He, C.; Wang, Y.; Zou, Y. An empirical model for landslide travel distance prediction in Wenchuan earthquake area. Landslides; 2014; 11, pp. 281-291. [DOI: https://dx.doi.org/10.1007/s10346-013-0444-y]

17. Tian, L.; Bi, W.; Liu, J.; Dong, X.; Xin, A. Seismic and collapse analysis of a UHV transmission tower-line system under cross-fault ground motions. Earthq. Struct.; 2020; 19, pp. 445-457.

18. Tian, L.; Pan, H.; Qiu, C.; Ma, R.; Yu, Q. Wind-induced collapse analysis of long-span transmission tower–line system considering the member buckling effect. Adv. Struct. Eng.; 2019; 22, pp. 30-41. [DOI: https://dx.doi.org/10.1177/1369433218774961]

19. Zheng, H.D.; Fan, J. Progressive collapse analysis of a truss transmission tower-line system subjected to downburst loading. J. Constr. Steel Res.; 2022; 188, 107044. [DOI: https://dx.doi.org/10.1016/j.jcsr.2021.107044]

20. Kassem, M.M.; Beddu, S.; Qi Min, W.; Tan, C.G.; Mohamed Nazri, F. Quantification of the Seismic Behavior of a Steel Transmission Tower Subjected to Single and Repeated Seismic Excitations Using Vulnerability Function and Collapse Margin Ratio. Appl. Sci.; 2022; 12, 1984. [DOI: https://dx.doi.org/10.3390/app12041984]

21. Li, J.; Gao, F.; Wang, L.; Ren, Y.; Liu, C.; Yang, A.; Yan, Z.; Jiang, T.; Li, C. Collapse Mechanism of Transmission Tower Subjected to Strong Wind Load and Dynamic Response of Tower-Line System. Energies; 2022; 15, 3925. [DOI: https://dx.doi.org/10.3390/en15113925]

22. Fu, X.; Li, H.N.; Li, G.; Dong, Z.Q. Fragility analysis of a transmission tower under combined wind and rain loads. J. Wind. Eng. Ind. Aerodyn.; 2020; 199, 104098. [DOI: https://dx.doi.org/10.1016/j.jweia.2020.104098]

23. Guan, R.; Xiang, C.; Jia, Z. Anti-Wind Experiments and Damage Prediction of Transmission Tower under Typhoon Conditions in Coastal Areas. Energies; 2022; 15, 3372. [DOI: https://dx.doi.org/10.3390/en15093372]

24. Zhou, Y.; Sheng, Q.; Chen, J.; Li, N.; Fu, X.; Zhou, Y. The failure mode of transmission tower foundation on the landslide under heavy rainfall: A case study on a 500-kV transmission tower foundation on the Yanzi landslide in Badong, China. Bull. Eng. Geol. Environ.; 2022; 81, 125. [DOI: https://dx.doi.org/10.1007/s10064-022-02628-9]

25. Wang, F.Y.; Xu, Y.L.; Qu, W.L. Multi-scale failure analysis of transmission towers under downburst loading. Int. J. Struct. Stab. Dyn.; 2018; 18, 1850029. [DOI: https://dx.doi.org/10.1142/S0219455418500293]

26. Zheng, H.D.; Zhang, B.; Wang, S.B.; Zhou, G.Z. Effects of Blade Numbers on Wind-Induced Fatigue Lives of Straight-Bladed Vertical-Axis Wind-Turbine Tower Bases. Metals; 2022; 12, 321. [DOI: https://dx.doi.org/10.3390/met12020321]

27. International Electrotechnical Commission. Wind Turbines-Part 1: Design Requirements. IEC 614001 Ed. 3; International Electrotechnical Commission: Geneva, Switzerland, 2006.

28. Higham, N.J. Analysis of the Cholesky Decomposition of a Semi-Definite Matrix; Oxford University Press: Oxford, UK, 1990.

29. Zheng, H.D.; Fan, J.; Long, X.H. Analysis of the seismic collapse of a high-rise power transmission tower structure. J. Constr. Steel Res.; 2017; 134, pp. 180-193. [DOI: https://dx.doi.org/10.1016/j.jcsr.2017.03.005]

30. Zheng, H.D.; Fan, J. Phenomenological hysteretic model for steel braces including inelastic postbuckling and low-cycle fatigue prediction. J. Struct. Eng.; 2019; 145, 04019037. [DOI: https://dx.doi.org/10.1061/(ASCE)ST.1943-541X.0002319]

31. Yang, S.C.; Hong, H.P. Nonlinear inelastic responses of transmission tower-line system under downburst wind. Eng. Struct.; 2016; 123, pp. 490-500. [DOI: https://dx.doi.org/10.1016/j.engstruct.2016.05.047]

32. Long, X.; Wang, W.; Fan, J. Collapse analysis of transmission tower subjected to earthquake ground motion. Model. Simul. Eng.; 2018; 2018, 2687561. [DOI: https://dx.doi.org/10.1155/2018/2687561]

33. Zheng, H.D.; Fan, J. Analysis of the progressive collapse of space truss structures during earthquakes based on a physical theory hysteretic model. Thin-Walled Struct.; 2018; 123, pp. 70-81. [DOI: https://dx.doi.org/10.1016/j.tws.2017.10.051]

34. Lewis, W.J.; Jones, M.S.; Rushton, K.R. Dynamic relaxation analysis of the non-linear static response of pretensioned cable roofs. Comput. Struct.; 1984; 18, pp. 989-997. [DOI: https://dx.doi.org/10.1016/0045-7949(84)90142-1]

35. Zheng, X.Y.; Zheng, H.D.; Lei, Y.; Chen, H. Nonlinear stochastic responses of a newly developed floating wind-solar-aquaculture system. Ocean. Eng.; 2021; 241, 110055. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2021.110055]

36. Zheng, H.D.; Zheng, X.Y.; Lei, Y.; Li, D.A.; Ci, X. Experimental validation on the dynamic response of a novel floater uniting a vertical-axis wind turbine with a steel fishing cage. Ocean. Eng.; 2022; 243, 110257.

37. Duan, Y.F.; Wang, S.M.; Yau, J.D. Vector form intrinsic finite element method for analysis of train–bridge interaction problems considering the coach-coupler effect. Int. J. Struct. Stab. Dyn.; 2019; 19, 1950014. [DOI: https://dx.doi.org/10.1142/S0219455419500147]

38. Higham, D.J.; Higham, N.J. MATLAB Guide; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2016.

39. Veletsos, A.S.; Darbre, G.R. Dynamic stiffness of parabolic cables. Earthq. Eng. Struct. Dyn.; 1983; 11, pp. 367-401. [DOI: https://dx.doi.org/10.1002/eqe.4290110306]

40. Loredo-Souza, A.M.; Davenport, A.G. The effects of high winds on transmission lines. J. Wind. Eng. Ind. Aerodyn.; 1998; 74, pp. 987-994. [DOI: https://dx.doi.org/10.1016/S0167-6105(98)00090-7]

41. Takeuchi, M.; Maeda, J.; Ishida, N. Aerodynamic damping properties of two transmission towers estimated by combining several identification methods. J. Wind. Eng. Ind. Aerodyn.; 2010; 98, pp. 872-880. [DOI: https://dx.doi.org/10.1016/j.jweia.2010.09.001]