The Spectral Approach of Love and Mindlin-Herrmann Theory in the Dynamical Simulations of the Tower-Cable Interactions under the Wind and Rain Loads

: The paper presents the transmission line dynamic response considering tower-conductor cable coupling under wind and rain load is presented. The cable and the tower are modeled using the Spectral Element Method. The ﬁnite Element Method with different types of elements is used for veriﬁcation. Excitations, such as wind and rain, are researched in the analysis and the results show the variation of the OTL displacements in each case. The tower modulates the dynamic cable response close to the coupling point and has a high inﬂuence in the response of the system. The application of the methodology based on spectral modeling, conﬁrming good feasibility for the overhead transmission line, tower, and cable is the purpose of the study. The new approach of Love and Mindlin-Herrmann theory in the dynamical simulations and tower-cable interactions depending on the distance from the cable connection point is studied. Analysis of the wave number depending on the frequency is also shown.


Introduction
Overhead transmission lines (OTL) are the worldwide standard system used to transport energy over long distances. The simultaneous action of wind and rain on infrastructure facilities such as overhead transmission line is a phenomenon often encountered, both in the case of normal wind and heavy wind, with values exceeding the speeds indicated in the wind standards, e.g., Eurocodes, such as tornadoes or hurricanes, which accompanied by heavy rainfall. Another important aspect is the aerodynamic phenomena that accompany the wind flow around the OTL under the influence of wind and rain as single lines or in groups, particularly turbulence effects. Dynamic analysis of the OTL response allows for the comprehensive identification of phenomena that affect the durability of infrastructure. From the point of view of the conducted analyses, it is worth paying attention to in-situ research, laboratory tests, and numerical simulations. The literature review is presented in this wide context, confirming that the discussed topics are still relevant and that the presented issues require further analysis. The methodology proposed by the authors of this paper and the used tools, supported by the results of the analyzes, lead to the conclusion that the paper is the next step in the investigation of the response of OTL under wind and rain excitation, revealing another picture of OTL behavior, not yet presented in the literature. Confirming the above statement, in the paper [1], an OTL with three spans was analyzed with rain and wind loads. The horizontal component of the load acting on the line was the sum of the horizontal drop-impinging force and the force of the wind. The authors pointed out the need for further research on the drag parameters of the conductor with the excitation of rain and wind, as there was insufficient information in the literature on that damping parameters of a conductor. In [15] a numerical model, taking into consideration the ice shedding was analyzed. Fluid-structure interaction methodology was used by Keyhan et al. [16]. The method yielded a more accurate representation of pressure loads acting on moving conductors than provided by the pseudo-static pressure calculation based on Bernoulli's equation, which is the current approach used in the design. The results based on the proposed method were compared to those obtained using the Bernoulli load model using four natural wind records to perform a nonlinear dynamic analysis. Yin et al. [17] presented the used vibration data, measured from sensors to detect structural damage in the transmission towers. The final load patterns for the tower structures using FEM were presented in the [18]. In [19] authors underlined that the forces from the cables to the tower cannot be neglected in the damage identification analysis. Tian et al. [20] analyzed the broken lines, ice, and wind as an additional load. The analysis was based on an explicit algorithm. The advanced static and dynamic analysis was presented in [21]. Authors in [22] gave practical information on modeling techniques to be used for lattice structures. For the purpose of the experiments, an 8 m long section of a transmission line tower was built in the laboratory, pulled at different levels of solicitation, and left to vibrate freely after the load was suddenly released. Numerical modeling was also conducted and compared to the experimental results. In the paper [23], the results of full-scale measurements and a time series analysis for the wind-induced vibration of a transmission line system in terms of the effects of the coupling motion between a steel tower and conductors were described. The authors emphasized that the coupling response characteristics were prominent in a longitudinal direction and were affected by the modeling manner of the supporting condition of the end of the conductors. The results of the paper [24] presented the differences in the response characteristics, and the peak factors computed from a time-series response were greater than those computed from power spectrum density. In the paper [25], the application of the pendulum damper was analyzed to reduce the vibration coming from the wind. In [26] the influence of a single wire on vibration was investigated. The specific behavior of the OTL system could be observed under heavy wind loads. Hamada and El Damatty [27,28] presented influnece of tornado on OTL. Li [29] carried out a wind tunnel experimental test in the OTL prototype and developed a probabilistic analysis monitoring possible failure in the structure submitted to different periods of wind return. Therefore, fast and accurate techniques have been explored for the dynamic and monitoring analysis of OTL. Computational fluid dynamics, especially in the field of drag and lift forces, is an important source of information supplementing wind tunnel research. Meynen et al. [30], based on Navier-Stokes equations investigated the vibrations due to the flow around the cable. Gołebiowska and Dutkiewicz [31] modeled the flow past the Stockbridge-type damper attached to the overhead transmission line determining aerodynamic drag and lift coefficients and the pressure around them and in their wakes. In [32,33], the element shape functions were obtained from the analytical solution of governing differential equations and the dynamic system solution was written in the frequency domain. In [34] the spectral and finite element method was used for the analysis of cable's damage in relation to frequency response functions. Dutkiewicz et al. [35] investigated the SEM model due to a change in section area and axial force acting on the cable. In [36] authors presented the comparison of the cable's natural frequencies resulting from measurements and numerical simulations. The research included, among others, wave number sensitivity and damping. In [37] the Wittrick-Williams algorithm for solving the transcendental eigenvalue problem of the line was presented. The application of SEM for rods and high-order spectral elements analysis were presented in works [32,[38][39][40]. This work presents a transmission line dynamic analysis with the application of SEM and FEM verification. The purposes of the research are addressed to the major issues: (i) to describe the truss spectral elements using different rods order elements and compare the efficiency of each approach to model the tower; (ii) to use SEM to model the OTL combining truss and cable elements; (iii) to demonstrate the dynamic coupling behavior of the tower and connecting cable; (iv) to predict the response of the OTL with rain and wind loads.

Structural Spectral Element
This section introduces the spectral elements formulation of rods and cables used in the overhead transmission model. The rod and cable equations were derived from the exact solutions of each element described from the solution of the wave equation [33].

Truss Spectral Element
Truss spectral element is described in the local coordinate system [41]. It is necessary to transform the local into global coordinates. Figure 1 shows the nodal displacements for a single elementary and Love rods on the left-hand side and Mindlin-Herrmann on the right side. The Love theory that per node, the one-degree of freedom assumption was made for the bar, and the transformation matrix has the following form: in the case of the Mindlin-Herrmann rod and cable structure it is necessary to apply the transformation matrix corresponding to a frame, expressed as The assembling of the matrix into the global dynamic stiffness matrix is given by

Elementary Rod Spectral Element
The theory of the elementary rod considers axial load and longitudinal displacement as represented in Figure 2. The rod's undamped equilibrium equation in the frequency domain is written as where parameter ω is the circular frequency, u is the longitudinal displacement, ρA is the mass per unit length, and EA is the axial rigidity. Hysteretic damping of structure is assumed as E = E(1 + iη), where η is the damping factor and i = √ −1 [32]. The general solution of Equation (4) is assumed to be u(x) = ae −ik(ω)x , where a is the indeterminate amplitude at each frequency. By substituting the general solution it into Equation (4) it gives the longitudinal wavenumber k. For the simple case where all properties are uniform, it is k = ω ρ/E. For the element of length L, the general solution is where a i , i = 1, 2 are constants. The propagation is non-dispersive The boundary condition are at u 1 = u(x = 0), and u 2 = u(x = L) and by applying into Equation (5) Substituting Equation (5) in Equation (7) the constant of boundary conditions vector a is eliminated given the spectral nodal displacement u(x, ω) = g(x, ω)d, where the interpolation function is g(x, ω) = e(x, ω)B −1 (ω) The internal axial force of the bar in spectral form is given by, and the spectral nodal axial forces for the finite elementar rod are defined as, where and the spectral element matrix for the elementar rod element is Or, in an explicit form, it has

Love Rod Spectral Element
By assuming the general solution u(x) = ae −ik(ω)x , one can rewrite the equation of motion in the frequency domain as, The wavenumber can be obtained by grouping similar terms in the equation, it has By following a similar procedure described for the elementar rod theory, the wavenumber can be expressed The spectral nodal displacement d, matrix of boundary conditions B, and the interpolation function g(x, ω) for the Love rod are similar to the elementar model presented from Equation (5) to Equation (8). The spectral nodal forces for the Love theory are found to be in matrix form expressed as and by relating spectral nodal displacements with nodal forces using the boundary condition term a, one has F = ΦB −1 d, where the dynamic stiffness matrix for Love rod will be S(ω) = ΦB −1 .

Mindlin-Herrmann Rod Spectral Element
In the analysis each node has two degrees of freedom, with longitudinal u and transversal v displacements, where v = ψ x y, ψ represents the transverse contraction [32]. The two nodes Mindlin-Hermann rod element is presented in Figure 3. The equations of motion for the Mindlin-Herrmann rod is The associated boundary conditions are specified in terms of where I is the inertia of the cross-section, µ, λ, K r1 , and K r2 are a set of coupled equations defined by For the two dependent variables u and ψ it is assumed constant coefficients, the general solutions will be From the general solution and replacing Equation (20) into Equation (19), in a matrix form, one has with the determinant of the matrix equal to zero, it has the characteristic equation expressed by the polynomial where, By solving the polynomial equation of Equation (24), it has the wavenumber As a result of the quadratic terms of the k 2 derived from Mindlin-Herrmann theory, it predicts the two modes related to the real (k 1 ) and imaginary (k 2 ) wavenumber. The general solution for longitudinal and transversal displacements is expressed as where R is relating term of the amplitude given by [38]: the subscript i varies between 1 and 2. The nodal displacements d for the Mindlin-Herrmann rod with length L are given in the form as The coefficients a i , i = 1, . . . , 4 can be found by substituting the element boundary conditions of the element nodal displacement, in a matrix form B will be, where p 1 = e −ik 1 L and p 2 = e −ik 2 L . By replacing the derivative of the displacement solutions presented in Equation (26) in the boundary conditions of Equation (20) where Substituting the vector a of Equation (28) in Equation (29), the dynamic stiffness matrix representing the Mindlin-Herrmann rod is determined likewise in elementar and Love theories.

Cable Spectral Element
The cable model derived from an Euler-Bernoulli beam subject to an axial load is shown in Figure 4. It node has flexural and rotational degrees of freedom and axial tension, shear, and momentum forces. By considering free vibration and hysteretic damping, the spectral form of the equation of motion is and replacing the spectral form of solution as v j = v j e −i(kx−ωt) in Equation (30), the characteristic equation for the wavenumbers k j , will be [32,42] There are four roots, representing two sets of wave mode pairs, being two real and two imaginary roots of the form For the Euler-Bernoulli beam spectral element subjected to an axial load of length L, we assumed that ±k 1 = k and ±k 2 = ik. The traveling waves can be expressed in the form for The spectral nodal displacements and slopes are as follows where, for the boundary condition terms a being expressed as a = B −1 (ω)d. Spectral shear V and moment M forces presented in the cable are defined as by solving the derivative and applying the boundary conditions, the spectral forces are obtained by: and the dynamic stiffness matrix obtained as in the rod elements, so that The displacement is interpolated from the nodal displacement vector d by eliminating the constant vector a, v where the shape function is When the beam is uniform without any sources of discontinuity, it can be represented by a single spectral element with accurate solutions [32].

Distributed Load
The distributed load can be incorporated into the spectral formulation and determined its contribution by multiplying the differential equation for v(x) by the function of the form g(x) and integrating between the limits of the length of the element [32], described by (41) where represents the spatial derivative representation. In this way, the stiffness ratio for the cable with a distributed load is written as The same as in the FEM [43], the elements in SEM can be assembled to form a global structure matrix system [33].

Wind and Rain Loads
The speeds and directions of wind and rain flow change in space and time. In practice, the loads resulting from these phenomena are determined from numerical simulations, field simulations and wind tunnel tests [20,29,30].

Wind Load
The action of the wind in contact with a body generates aerodynamic and friction forces. The difference in pressure downstream and upstream of the structures causes drag pressure, which is considerably greater than the viscous drag of the air. Wind speed changes in different altitudes and varies in time [44]. It is assumed that the average wind speed is stable in the period of 10 min to 2 h as demonstrated in [44], in this interval, the instantaneous speed can be obtained by, where V m is the average speed and v the speed variation, those parameters will change from geographic location, e.g in Brazil the forces due to wind admits the average speed remains constant for a period of up to 10 min [45]. However, fluctuations in the wind can introduce important fluctuations in the direction of the wind acting in slender structures, such as in cables. The basic wind speed V b is defined as the wind speed referred to a return period of 50 years [46], at the height of 10 m in relation to the ground, referring to an integration period of 10 min measured in terrain with a degree of roughness B type.
In transmission line analysis the design speed V p is computed by relating the correction factor V b adopted according to the degree of roughness of the terrain, the wind action time interval, and the height of the obstacle, so that where h is the wind action height level, K r roughness coefficient,K d relationship between integration period and terrain roughness, the n the correction factor. The dynamic reference pressure is given by for ρ v = 1.293 (1+0.00367)λ 16000+64λ−κ 16000+64λ+κ being the air density, λ the temperature [ • C ], and κ the altitude of the region. The wind action on the cable generates a perpendicular force, and at the point of attachment to each support is given by where C a represents the drag coefficient equal to 1 determined by the number of Reynolds, α is the effectiveness factor, d the diameter of the cable in meters, L length of span considered in meters and θ the angle of incidence of the wind (≤90 • ). The technical standard NBR6123 [45] defines the action of winds on buildings, with the global force in the wind direction as F v = C a P v A, where A is the area under the effect of the wind load, and the dynamic pressure P v is obtained as

Rain Load
The mathematical expression of the rain load considers the raindrops' motion and the law of conservation of the moment [3]. It is assumed that the raindrop splashing is neglected in the process of the raindrop impact. The trajectory of the raindrop equation for vertical (y), cross-wind(z) and horizontal along-wind (x), directions are given by [2] where m = 4πr 3 3 ρ w is the raindrop mass, U, V, W are the components of wind speed, ρ w is the water density, ρ v the air density, µ the air viscosity, r is the raindrop radius and R ey the Reynolds number given as By replacing Equation (49) in (48), it has In the wind field of the atmospheric boundary layer, raindrops are influenced by the wind flow. After a simplification of the drag coefficients in the all directions, the components C ax , C ay , and C az are equal to the sphere drag coefficient, approximated by 0.47 [47]. The motion of the drops is due to the inertia [3]. The variables characterizing the rain load in the vertical and horizontal speed direction are the specific capture rate N(d) related to the raindrop diameter d r , and the speed ration γ(h, d r , α). In short, γ(h, d r , α) is proportional to d r of the raindrop and to the power law exponent α, while inversely proportional to its height h, expressed as The relationship between the intensity of rain in the horizontal and vertical direction is given by the ratio N(d) which is also equal to the rate of rain flow per unit area at the same orientation, given by where V h is the horizontal velocity and is the vertical terminal velocity. For V term greater than V h , the vertical rain velocity will be greater than the velocity in the horizontal direction. For rain, with a drop diameter, the intensity is defined by the variable V drop is the velocity that is perpendicular to the plane of the raindrop spectrum and n(d) is the number of raindrops with diameter per unity of volume. Substituting V term (d) and R(d) into Equation (52) the catch ratio can be expressed as where the spectrum in the horizontal and vertical direction are n h = n v = n 0 e (−Λd) , receptively [48]. For n 0 = 8 × 10 3 in 1/(m 3 mm) and Λ = 4.1R −0.21 h in 1/mm, R h is the intensity of the horizontal rain. In the process of colliding the raindrop with the structure, its speed V h becomes zero. Since the collision time is based on the moment conservation law, it has the average impact force F(τ) of the drop over a time interval τ is The rain load with a specific diameter acting on the structure in a unitary volume must adopt the spectrum in the vertical plane so that F d = F(τ)n v (d).
Therefore, the load related to rain acting on the structure is obtained by [3]: where ∀ is the volume of the flow acting on the structure, and the rain pressure is given by: In situ, it is typically that rain follows by wind, hence there is the necessity to study simultaneous these events to better understand the behavior of the structure [1]. The wind and rain loads acting simultaneously are developed as the sum of forces yields (60)

Numerical Analyses and Discussion
The numerical analysis examined the results regarding the tower, cable, and cabletower coupled model under the wind and rain loads. The span of the transmission line of 100 m is shown in Figure 5. The tower is based on the structure proposed by [19] which is a self-supporting tower composed of rods of the same cross-sectional areas and parameters: modulus of elasticity E = 210 GPa, a density of ρ = 7860 kg/m 3 , the cross-section area A = 0.01 m 2 , and Poisson coefficient equals to 0.3. The tower is 26 m high and contains 48 rods in total with 22 connection nodes. The cable was made of aluminum with the properties E = 74 GPa, ρ = 2700 kg/m 3 , η = 0.01, area of A = 50 × 10 −5 m 2 , length L = 100 m and a tension load of 27 kN. The suspension string was omitted, therefore, this article does not analyze the issue of the swing of the cable. In the numerical model, the truss and the cable elements are used. In Figure 5, the analyzed points on structures are shown: point (a) is located on the tower, (b) is the coupling point, point (c) is on the cable, 1 m from the coupling point, point (d) is located in the middle of the cable. The assumed parameters of wind and rain load are: base wind speed of 20 m/s, air density of 1225 kg/m 3 , raindrops density of 1000 kg/m 3 , correction rate n = 12, terrain category B, roughness category K r = 1, K d = 1.6, temperature 25 • C, drag coefficient C a = 1, power-law coefficient α = 0.65, rain intensity of 150 mm/h.

Tower Dynamic Analysis
The overhead transmission tower is modeled with the rod and beam structural elements as a small lateral rotation in vibration conditions. Different types of rods are used for modeling the tower and verifying the vibration response, as well as to investigate the computational cost of each approach. Figure 6 shows the dispersion diagrams for a single elementary and Love rods on the left-hand side and Mindlin-Herrmann on the right side with properties described in Section 5. Elementary rod has a non-dispersive behavior in the whole frequency range. The Love bar has a non-dispersive behavior up to 100 kHz and from this frequency onward a dispersive one, caused by the lateral contraction that is considered in this theory. The Mindlin-Herrmann dispersion diagram shows the first (k 1 ) waveform that presents a purely real behavior. The second mode (k 2 ) shows a complex waveform with a cutoff frequency ω c around 140 kHz. It has an imaginary value and therefore, evanescent character. Above ω c it presents a real waveform with propagating waves.   The time processing increases as the complexity of the theory, which is shown in Table 1. The analysis with the application of the Mindlin-Herrmann model increases the computation time in 107% and the Love model in 53%, in comparison with an elementary rod. Therefore, to select the rod model used in the tower, it is important to consider the frequency range used, the deformation component conditions, and computational time. The transmission tower is modeled with the three-rod theory and compared to the FEM model aiming to verify the application and lead to the rod selection. For the FEM tower analysis, it is used the ANSYS software and Link 180 elements. Figure 8 shows the tower's receptance obtained for the excitation with unitary force and measured at point (a) that is shown on Figure 5. The tower modeled with elementary (El), Love, and 0.01 m mesh size FEM elements has a similar dynamic response in the frequency of 0 to 200 Hz. The FEM analysis required a larger number of elements in the mesh compared to SEM models. The resonance frequency for Mindlin-Herrmann (M-H) elements differs from the third mode compared to the other models. The mode shapes are the same. Mindlin-Herrmann 1.4402 Figure 8. Tower receptance at point (a) as located in Figure 5, with elementary, FEM, Love, and Mindlin-Herrmann models Table 2 shows the first five frequencies of resonance. The results have good approximation among the theories, although FEM needed a refined mesh to converge with SEM and elementary rod. Table 3 shows the mesh discretization ratio of the tower modeled by FEM related to the resonant frequencies. By comparing SEM single rod and FEM, the best outcoming of FEM is with an element sizing 0.01 m. Figure 9 brings the FEM convergence analysis related to the percentage error, where it is shown that the error decreases with the increase of the mesh discretizations. The convergence percentage error was calculated by subtracting the SEM (ω S ) resonance frequency from the FEM (ω F ) and divide by (ω F ), as [(ω F − ω S )/ω F ] × 100.  The computational time for processing is shown in Table 4 and the simulations were performed using a Dell microcomputer with an Intel (R) Core (TM) i7 processor and 8 GB RAM. Note that the elementary rod had the shortest processing time. Comparing the SEM and FEM solutions, it was demonstrated the mesh convergence was proportional to the computational time. The SEM proved to be an accurate method with a low computational cost for all approaches. Because the results with the different SEM rod theories show small variations in the analyzed frequency band and economical cost involved in the analysis, the elementary theory shows to be the best choice for the SEM rod, and it will be used in the next studies.

Cable under Punctual and Distributed Load
The cable is made of aluminum with the mechanical properties described in the previous subsection. The cable dispersion diagram demonstrates the direct relation of the tensile load with the propagation wave type. Increasing the tension, the wave propagation behavior changes from dispersive for a low tensile load to non-dispersive as shown in Figure 10.  Table 5 shows the cable resonance frequencies obtained during the analysis: analytical [49], SEM, and FEM with different mesh discretization. Presented resonance frequencies are for the 100 m cable span with fixed ends. The frequency analysis was done in the static equilibrium configuration under self-weight. The SEM presents the same frequency values compared to the analytical solution. FEM solutions were obtained with 10, 20, 50, and 100 elements in the mesh. For FEM, with the increase in the number of elements, the accuracy increases together with the processing costs. Table 5. Cable resonance frequencies obtained with the spectral element (SEM), analytical solution(An), and finite element method (FEM) calculated with different mesh discretization.  Figure 11 shows the receptance response obtained at the left end of the cable due to a unitary excitation imposed at the same measured point for SEM and FEM with 20 elements in the mesh. The resonance frequencies and mode shape up to 10 Hz are similar, after 10 Hz the resonances have a small difference. In the next case studies, the results of the analysis with the SEM are presented. The influence of the direction of the wind and rain acting on the cable was analyzed.  Figure 12a presents the receptance function of the cable exposed to the wind-rain load acting at 45 • and 90 • angle. Figure 12b shows the results from analysis of the cable that is excited by the wind, rain, and wind-rain as the punctual forces. The basic wind speed was assumed as 20 m/s, corresponding to the category of land as urban areas and lands with many tall trees.

Distributed Load
The cable dynamic response changed with the distributed load. Figure 13a-d shows the receptance FRFs of the cable subjected to a tensile load of 27 KN, excited by (a) distributed unitary load; (b) distributed wind excitation with a base speed of 20 ms; (c) distributed rain load with an intensity of 150 mm/h, and (d) distributed wind and rain. The analysis estimated the receptance over a meter of cable from the support and considered distributed load according to the formulation presented in Section 3.
The cable receptance presented a significant difference when it considered the distributed load compared to the punctual load. Distributed load affects the amplitude of the response and the vibration of the structure itself. However, by comparing the wind and rain-distributed loads with unitary, the difference occurs only on the response amplitudes. The following analyzes assumed distributed loading with wind and rain loads simultaneously.

Overhead Transmission Line Tower-Cable Coupling
The overhead transmission line composes of a tower coupled to the cable as shown in Figure 5. The FRFs of the system were obtained at the points 'a' to d due to a unitary excitation force imposed at point 'a' (red dot). Figure 14a-c are the FRFs measured at the point 'a' (red dot) horizontal and vertical direction in the tower and b (yellow dot) the coupling point, respectively. The receptance response of those three measurements captured only the tower dynamic resonance frequencies and modes, even in the coupling point where there is the cable attached. The FRF obtained on the cable from a meter of the coupling point is presented in Figure 14d, the interaction of the cable and the tower is observed, and noticed the dynamic behavior of the tower modulates the cable's response. FRF calculated at the cable middle length, and the end-span are shown in Figure 14e,f, where the cable dynamics increased, but are still modulated by the tower vibration. For the analyzed problem, the range of frequencies from 0 to 200 Hz is significant. The fatigue failures in the transmission line occur close to the connection coupling point. Figure 15 shows the cable's receptances obtained from the coupling point (c) to 1 m (d) for every 0.1 m distance. At the coupling point, FRF is dominant. In the other points, the cable's response contribution increased gradually but was still modulated by the tower modes. The tower's response has a great contribution to the OTL vibrational response. As shown in the literature, cable fractures usually occur between 0.6 to 1 m from the coupling point. In this region, there is an increasing influence of the cable inducing local vibration fatigue, as demonstrated in Figure 15.  Table 6 shows the first six frequencies calculated in the tower at the connection point, and in the cable at 1 m and 50 m from the connection. One can notice the interaction between the tower and the cable by observing the resonance frequencies. The first frequency within one meter of the connection coincides with the tower, and the others are cable frequencies far from the tower one. In the middle of the cable (50 m from the connection), the first six frequencies are low. The influence of the excitation position is presented in Figure 16. It assumed a unitary force imposed at the coupling point 'b' (black curves), a meter from the connection in 'c' (red curves), and at the cable middle length in 'd' (blue curves). The receptances shown in Figure 16a-e were calculated in the tower, coupling point, in the cable at 1 m and 50 m, and at the end of the cable. The receptances obtained due to the excitation at point 'b' held a significant influence from the tower in all cases. In contrast, the receptance derived due to the excitation at point 'b' presented a notable influence of the cable responses and high modal density. However, the tower dynamics modulated the curves. By excitation the system in the middle of the cable, the receptances demonstrated a greater influence of the cable with greater amplitudes in its responses. Therefore, to minimize the influence of tower vibration on the cable response, it is recommended to select the excitation and measurement points over a meter of the connection.

Distributed Load
FRFs obtained over the tower and at the connection point shown in Figure 18a-c assumed a unitary distributed load over the cable. The responses maintain the amplitude of displacement combined with the dynamic tower response, and it is predominantly influenced by cable vibration. The FRFs obtained in the cable, Figure 18d-f exhibited a low influence of the tower in its dynamic behavior. In Figure 18d close to the coupling the FRFs increased in amplitude, and in Figure 18e,f related to the furthest points from the coupling decreased the amplitude.  Figure 19 presents the FRFs estimated using a distributed wind-rain load excitation over the overhead transmission conductors. By comparing Figure 19 to Figure 18, a similar pattern can be observed with changes only in amplitude due to the increased load. For the excitation with the wind and rain phenomena acting separately, the responses follow the same curve pattern, with changes only in amplitude.

Conclusions
The action of the rain and wind excitation acting at the same time on the tower and the line is performed in the paper. The behavior of the tower-cable system is presented as excited by the rain and wind, which is more sophisticated than the wind load only. The average impact force of the drop is considered. The purpose of this analysis is to compare the results for Mindlin-Herrmann, Love, and the elementary elements. As models have the wave character, for the purposes of comparison, the propagation of the wave number depending on the frequency is shown. The impact of the load variation causes high displacement amplitudes of the system. It is also notable, that there are no significant variations in the FRFs behavior in the case of distributed loads. Effects of the OTL dynamic coupling responses are evident with changes in the frequencies of a modulation of the cable's FRF by the tower response. It is worth mentioning that it is essential to analyze the interaction of the structures since both presents the influence of the coupling in the responses. The results prove that, closer to the connection point, the couple vibrations combine the tower and cable response, which means that these elements influence each other near this connection, but far from this point this interaction decreases. The results show that the Poisson coefficient does not affect the tower's behavior. Therefore, the elementary bar theory in the analysis is further assumed. This is important to emphasize that results of SEM simulations are close to the real vibrations of the cable. Simulation validation was performed and discussed in [37].
The proposed methodology of the spectral and finite element methods is suitable for cable and tower analysis. The method used in the present work shows good feasibility for this type of structure. The transmission tower modeled by the SEM with the three types of rod obtains similar responses to FEM, which has a large mesh in its discretization to achieve a good result. For the cable analysis, the comparison is made among the SEM, FEM, and the analytical solution, demonstrating the good accuracy of SEM.
In the next research, the influence of the shedding of the ice, shock load, and eccentricities of the load at joints of the tower will be analyzed.