1 INTRODUCTION
In the 21st century, the global water crisis and environmental issues resulted in an accelerated growth in the deployment of increasingly large onshore and offshore wind turbines, with increasingly high towers, in search of technically feasible and economically viable winds. Therefore, the technical development, trade and installation of wind turbines in the world has advanced fastly, so that the generation of energy from thermoelectric, nuclear and hydroelectric power plants has been complemented and/or replaced by the production of such equipment (^{Hau, 2006}).
The global growth of installed wind capacity in recent years has established an annual record of 63,633 MW in 2015; resulting in a cumulative installed wind capacity of 539,581 MW at the end of 2017 (^{GWEC, 2018}). In Brazil, wind power production reached 14.34 GW of installed capacity in 568 wind farms and more than 7,000 wind turbines in 12 states, with the states of the Northeast Region, where the records of attendance exceed 70% of the energy produced in this region, which account for most of Brazilian production (September data according to ^{Abeeólica, 2018}). This shows that the cost of wind energy is the most competitive and viable amongst the renewable energy sources compared to the conventional ones (^{Burton et al., 2011} and ^{Manwell et al., 2009}).
For this reason, Brazilian political, social and technical issues have been studied in order to make feasible and further develop the use of this type of energy to produce electricity (^{Juárez et al., 2014}). Specifically, the state of Pernambuco has carried out, over the last years, a series of actions, with the objective of fomenting the renewable energy sector (^{Atlas..., 2017}). Among these, we highlight the attraction of industries to the Suape Wind Pole, that resulted in the installation of a group of companies to manufacture wind turbines, tubular towers, large blades and steel forged. For example, the wind power plants of Santa Brigida located in Caétes, Pedra and Paranatama cities of Agreste region of Pernambuco and the wind power plants of Santo Estevão in Araripina city of Chapada region of Araripe; totalizing more than 600 towers of 80 m high for support wind turbines in operation in the state, evidencing the growth of this type of energy.
In view of the vast potential of wind power that can be used in the state of Pernambuco, the size of the wind turbines to be used in the future tends to be increasing and more powerful, making it necessary to install these equipments under the action of more intense and continuous winds and thus causing the dimensions of the towers of these wind turbines to be continuously increased. Particularly, the height of the tower is an essential parameter for capturing stable highaltitude winds; however, the cost of this, which can exceed 20% of the wind turbine overall cost (^{Hau, 2006} and ^{Yoshida, 2006}), makes the height increase a disadvantage. In addition, the transport and the assembly of the tower become more expensive. Therefore, the design and installation of larger wind turbines makes it necessary to have higher and higher towers that require more elaborate structural analysis, stability and dynamics, resulting in a more complex design.
Thus, increasing the towers height, the vibration effects on the wind turbine components are greatly increased, so that in the last decades researchers and manufacturers of wind turbine components have sought solutions to mitigate the vibrations from the operation of mechanical components and wind and seismic actions (^{Menezes et al., 2018}). Even in situations where these structures can withstand such dynamic actions, without necessarily suffering structural damage, the fatigue effects of the constituent materials should not be overlooked.
To achieve the balance between safety and economic efficiency in large wind turbine design, the structural control to mitigate excessive structure responses is a viable option. In this way, passive, active, hybrid and semiactive control systems are valuable alternatives for the control of vibrations of large wind turbine towers. Even considering the studies already carried out and the numerous practical applications of structural control in bridges, in telecommunications towers and in tall buildings, the research and application of control systems for wind turbine towers are relatively recent topics.
The use of Tuned Mass Dampers (TMD) has been the object of analysis carried out by several authors (^{Murtagh et al., 2008}; ^{Lackner and Rotea, 2011}; ^{Stewart and Lackner, 2013}; ^{Shzu et al., 2015}; ^{Colherinhas et al., 2015}; ^{Tong et al., 2015}; and ^{Zuo et al., 2017}). However, despite the studies made, there are fields of analysis that have not yet been fully approached, such as design of a control device with satisfactory performance and with clearly demonstrated execution conditions.
Hence, this paper proposes to analyze the dynamic behavior and design of a vibration control device for a 120 m high steel tubular tower of a 3.2 MW Horizontal Axis Wind Turbine (HAWT). Thus, the Hybrid Mass Damper (HMD), proposed in this research, is presented as a differentiated solution to the previously mentioned studies, since, instead of being allocated to the nacelle (^{Hu and He, 2017}), it was designed to be installed internally to the wind turbine tower (near the top). In this way, the internal space of the nacelle is saved for the installation of the electricmechanical components of the wind turbine (^{Zuo et al., 2017}). This has the advantage of facilitating access to the control equipment, contributing to improve the maintenance conditions.
The insensitivity to orientation (^{Stewart and Lackner, 2014}), combined with the restrained displacement of the HMD, suggests that it could just as effectively be located at the top of the tower rather than in the nacelle, which may be more convenient, because the HMD would not yaw with the nacelle. It would essentially experience all orientations, which does not appear to impact the load reductions.
It is also highlighted the development of constructive model of the DVA developed, in which its main components are described, in order to show the viability of design and execution of the vibration control system designed for the tower analyzed.
2 VIBRATION CONTROL THEORY
2.1 Passive control vibration
For the vibration control methodology used in this paper, the theoretical development of the TMD proposed by ^{Den Hartog (1947}) was initially used for systems with two degrees of freedom, one of which is necessary to describe the movement of the main mass (the structure to be controlled) and the second one referring to the secondary mass (mass of the vibration absorber). The TMD is a device consisting of a mass, a spring (or a combination of springs) and a damper (or a damper assembly) that are attached to the main structure in order to reduce/control its dynamic response. The frequency of the absorber is tuned in a certain frequency of the main structure. When this frequency is excited, the absorber vibrates out of phase in relation to the movement of the structure, dissipating the energy contained in the main structure by the movement of the absorber.
The absorber parameters (mass
Writing the motion equations of the system and applying the complex numbers for resolution when the main structure is subjected to a harmonic excitation, we obtain:
which is the steadystate amplitude of the mass
which is the steadystate amplitude of the mass
which,
where,
With the addition of the absorber to the system, it is intended to reduce the amplitude peaks of the main mass to the lowest possible values. The amplitude peaks of the main structure are infinite when the damping of the absorber is zero (resonance) or infinite. However, there is an intermediate value of damping for which these peaks are minimal, since the work performed by the damping force is maximum.
According to ^{Den Hartog (1947}), by evaluating the values
which is the tuning required as a function of the absorber mass.
Finally, the optimal damping values can be calculated by substituting (9) into (1), deriving the resulting expression from
Thus, by using the damping value above, there is an almost horizontal tangent at each of the fixed points. The stiffness value of the tuned absorber can be calculated using Equations (3) and (4):
The design values of the vibration absorber coefficients are calculated by considering that the tower movement is governed predominantly but independently by the first two natural modes of vibration (the first two flexural modes having practically equal vibration frequencies). Thus, the absorber is independently tuned to the two natural frequencies whose modes define movements according to the first two vibration modes of the tower. The optimal values of the absorber coefficients are then obtained as if these two modes of vibration of the tower were represented by two independent systems of one degree of freedom (dof) each, in order to make feasible the application of Den Hartog's theory.
The TMD is preferably allocated where the deflections of the structure are larger; in this case, the TMD motion is associated with an added translational dof near the tower top. Thus, the mass, stiffness, and damping matrices of the coupled system (towerabsorber) come from the addition of the aforementioned dof.
Therefore, by means of the modal mass, stiffness, and damping values, the coefficients corresponding to the modes of vibration of the main system to be controlled are obtained. Then, starting from the choice of mass ratio
2.2 Hybrid vibration control
The translational hybrid vibration absorber of two degrees (Figure 1b) is defined similarly to the dynamic passive vibration absorber with two degrees of freedom of translation. Hydraulic actuators were added to the system in order to convert the TMD into a Hybrid Mass Damper (HMD), which is applied to the tower top to obtain a hybrid control system.
This model is idealized as part of a full state feedback control system, that is, the actuators apply forces of magnitudes and directions preestablished by the controller. These forces are generated from the signals received from the sensors positioned in a way measure all model state variables. Thus, for a tower model with 8 beam finite elements that already shows good accuracy in relation to the structure of the excessively discrete tower, with respect to the most significant global vibration modes for dynamic final tower response, it is necessary to position the sensors every 15 m of tower length, for example.
Working in the field of analysis and design of modern control systems, it is convenient to represent the differential equations that generate the movement of a dynamic system in the states space (^{Moutinho, 2007}). This analysis involves three types of variables in its modeling, namely: state variables, input variables and output variables. The state variables of a dynamic system are quantities whose set of values determines the state of the system (^{Junkins and Kim, 1993}). Input variables are quantities related to external actions or control forces/moments applied to the structure. The output variables are related to the input or state variables that can be measured, since in many practical situations, not all states or inputs are available for measurement.
Considering the static condensation (^{Clough and Penzien, 2003}; ^{Chopra, 2012} and ^{Humar, 2002}) of the dynamic system that represents the behavior of the tower, in order to dynamically study the horizontal translation dof, the dynamic motion equation of the structureDVA takes the following form:
where:
In the hybrid vibration absorber is considered that external forces are not applied, but only a control force imposed by the actuator. Thus, the vector
the control force vector
and the matrix
in which, terms
Finally, for passing to the state space, we isolate the vector
obtaining the STATE EQUATION of the system:
which will be solved by applying a linear transformation diagonalizing the state matrix
and,
Finally, we have the output equation that defines the output vector
where:
2.3 Active vibration control
An Active Mass Damper (AMD) was implemented to evaluate the behavior of the system controlled by a purely active absorber and to have information of the energy expenditure necessary to perform the control of the tower by this type of absorber. The AMD of two degrees of freedom consists of a system composed of a mass guided exclusively by hydraulic actuators, without springs and dampers. Thus, from the HMD, the stiffness and damping matrices are obtained by assigning null values to the parameters of stiffness and damping of the absorber (that is,
2.4 Linear optimal control
In active and hybrid control systems is necessary to determine the values of the input variables that will be applied by the actuators. In this way, the design of stable control systems based on quadratic performance indexes (or cost function) will be considered in this item. We intend to choose the control vector
where
where
where:
The second term on the righthand side of Equation (24) accounts for the energy expenditure of the control signals, which will be a determining factor for the operation of the control system, because the intensity of the control force will be limited by the capacity of the hydraulic actuators.
To solve the optimization problem, we substitute (23) into (17):
in which, we assume that the matrix
According to ^{Ogata (1997}), we use the Lyapunov approach to solve this optimization problem:
where
replacing Equation (26) without the vector of external forces that may be associated with external disturbances acting on the system, and observing that this equation must be true for any
By ^{Lyapunov's second method (Lyapunov, 1992}), if
In order to obtain the solution to the quadratic optimal control problem, we write the matrix
where
a function of
which is the optimal gain matrix required. Finally, the
that is called the matrix Ricatti equation. Its resolution is obtained by the process of diagonalization, as presented in ^{Kwakernaak and Sivan (1972}).
3 TOWER MODELING
A S355J2 steel tubular tower (design by ^{Lima et al., 2018}; Figure 2) supporting an SWT3.2113 (^{Siemens, 2014}) wind turbune was considered, according to the characteristics specified in Table 1.
Parameters  

IEC (International Electrotechnical Commission) Class  IIA 
Rated power (MW)  3.2 
Rotor diameter (m)  113.0 
Blade length (m)  55.0 
Swept area (m^{2})  10,000 
Hub height (m)  79.5  142.0 (used 122.5 m) 
Power regulation  Pitch regulated 
Annual output at 8.5 m/s  14,402 MWh 
Nacelle weight (ton)  78 
Rotor weight (ton)  67 
The modeling of the tower was done using the Finite Element Method  MEF (^{Cook et al., 2002}), considering nonlinear geometric elastic behavior. Initially, a code was implemented in Mathcad v.15 (^{Mathcad, 2012}), in which the tower was modeled as a structure clamped at the base and free at the top (cantilevered), using 8 dof beam finite elements (4 dof per node: axial and transverse translations, flexional and torsional rotations). An analysis was made of the required level of discretization (convergence analysis). It was noticed that a discretization with about 16 beam elements is already sufficiently accurate to the modal analysis necessary for the design of the control device proposed.
Moreover, a finite element model was created in the software ANSYS r.14.5 (2012), in which the tower was also clamped in the base (Figure 3a) and discretized by 7272 shell elements with 4 node and 6 dof per node (Figure 3b). The nacelle was modeled with tetrahedral solid finite elements, designated by SOLID187 (Figure 3c), with 10 nodes and 3 translational dof per node; representing it by a uniform mass. The main reason that led to the use of two finite element models was the need to evaluate the reliability and accuracy of the numerical results obtained.
The results of the modal analysis performed by ^{Lima (2018}) for both the shell and beam finite element model are similar for the global vibration modes, mainly with respect to the 1st mode of vibration of the tower that is more susceptible to excitation. Therefore, the beam FE model can be used as a representative of the dynamic behavior of the tower when it is subjected to the harmonic excitations resonant to its fundamental frequency.
For excitation modeling, the terms amplitudes of the vector of forces applied horizontally to the tower were calculated using the parameters established by ^{ABNT (1988}): fluctuating component of the wind was established according to the difference between the peak velocity pressure and the main velocity pressure; and a horizontal force from the nacellerotor system applied to the top of the tower (according to ^{Asibor et al., 2015}).
The fluctuating component of the wind and the horizontal force transmitted by the rotor to the top were adopted according to a harmonic excitation considered resonant to the fundamental mode of vibration of the tower. In this way, it was considered that the tower vibrates around the equilibrium position determined by the action coming from the main wind speed (^{Franco, 1993}). The amplitude of the force applied to the tower top has an order of magnitude 100 times greater than the other amplitudes that corroborates with the excitation of its fundamental mode. It is also commented that the random effect of the wind action in the determination of the force vector was not considered, but only the intensity of the wind action that corresponds to its fluctuating component (wind gusts).
4 ANALYSIS AND RESULTS
4.1 Passive vibration control of tower
In order to design a dynamic vibration absorber with two degrees of freedom that is able to control the oscillations according to the lateral directions X and Z, it is assumed that the movement according to the direction X of the tower is governed by the first mode of vibration (1st flexional mode in the XYplane) and that the movement according to the Z direction is governed by the second mode (1st flexional mode in the YZplane). This fact is possible because it is a symmetrical structure, resulting in the first noncoupled modes of vibration (flexional XY and YZ). Thus, the absorber is designed to directly control the first two flexural modes of vibration of the tower, but having its utility directed to the control of the tower as a whole, since the vectorial composition of the stiffness and damping of the absorber allows the control in any direction transverse to the tower. Moreover, exclusively for the purposes of designing the absorber, the modal dampings associated with the 1st vibration modes of the tower are neglected in order to allow the application of Den Hartog's theory.
The modal mass and stiffness from the 1st tower vibration modes are given in Table 2. For the design of the absorber, the mass of the DVA is equal to 3% of the total mass of the tower (that is equal to 8.534 E5 kg), differently than some authors suggest that it is 0.5 to 1.0% (^{Housner et al., 1997}). The tower is a large structure, so a greater mass percentage is needed to attenuate its vibrations without there being an amplitude of displacement exacerbated of the DVA, because there is a space limitation at the top of the tower (the top diameter of the analyzed model is 3.5 m). Thus, the DVA mass was adjusted in order to maintain its vibration amplitude limited to the space within the tubular structure of the tower.
With the values of the modal mass of the tower and the mass of the absorber, the parameters
Structure type  Mass (kg)  Stiffness (N/m)  Damping (kg/s) 

Tower  1.37263 E5  4.79351 E5   
DVA  2.56014 E4  6.35067 E4  1.95797 E4 
The absorber mass of 25601.4 kg is composed of a steel box (1.00 x 1.00 x 2.25 m) filled with lead. Lead was used as fill material because of its high specific mass (11340 kg/m^{3}) and, consequently, reduction of the volume of the control device. Thus, given the parameters of the vibration absorber, the stiffness, mass, and damping matrices of the vibrationcoupled towerabsorber system are assembled.
From the spectral matrix
The function of the passive vibration absorber is to detune the fundamental frequency of the tower with respect to the proposed excitation frequency, since the fundamental frequency of the tower without vibration absorber, which is equal to 0.29745 Hz, becomes be equal to 0.22810 Hz, for the case of the towerDVA assembly. In addition, the 2nd vibration frequency (0.31889 Hz) is slightly higher than the fundamental frequency of the tower without absorber. However, such detuning is done by reducing the natural frequency of the structure, i.e., leaving the tower even more flexible. This low flexibility of the tower has been the subject of research that seeks to develop the socalled soft towers, in order to make them more competitive solutions in the market, since, great savings could be reached in case the dogma of great distances between eigenfrequency and exciting frequency is disregarded (^{Lange and Elberg, 2018}).
4.2 Hybrid vibration control of tower
It is possible to improve the dynamic behavior of the tower (to reduce the intensity and number of vibration cycles) when implementing the hybrid vibration control. Thus, with the passive vibration absorber being designed for the two translational degrees of freedom according to the X and Z directions, the same parameters for the dynamic hybrid vibration absorber are used and the design is complemented with the determination of the optimal control variables to be introduced by the hydraulic actuators.
For the determination of the system state feedback gain matrix
The control ponderation
4.3 Active vibration control of tower
The procedure for obtaining the gain matrix
4.4 Parametric analysis of control systems
In order to compare the performances of the passive, hybrid and active control systems, simulations were made of the response of the controlled tower, the output of the vibration absorber and the control forces as a function of the mass ratio
For plotting the graphs presented, the mass ratio was varied and the other parameters of the absorber, such as stiffness and damping, were recalculated for each mass ratio value, in order to maintain the tuning of the frequency of the absorber with the frequency of the tower, according to Table 3.
Mass ratio

Modal mass ratio

Frequency ratio

DVA stiffness

DVA damping

1.0  6.217  0.941  2.64153 E4  4.44880 E3 
2.0  12.434  0.889  4.71495 E4  1.15540 E4 
3.0  18.651  0.843  6.35067 E4  1.95797 E4 
4.0  24.869  0.801  7.64536 E4  2.79219 E4 
5.0  31.086  0.763  8.67169 E4  3.62791 E4 
7.5  46.629  0.682  1.03960 E5  5.63376 E4 
10.0  62.171  0.617  1.13317 E5  7.45715 E4 
12.5  77.714  0.563  1.17953 E5  9.08481 E4 
15.0  93.257  0.517  1.19692 E5  1.05310 E5 
17.5  108.800  0.479  1.19625 E5  1.18167 E5 
20.0  124.343  0.446  1.18427 E5  1.29632 E5 
In the graph shown in Figure 5a, the modal mass ratio
In Figure 6a is shown the graph of the displacement r.m.s. (root mean square) ratio between the values of the tower top controlled and without control, for the three types of control studied (purely passive, hybrid and purely active), as a function of mass ratio
Additionally, it can be observed that the passive control system presents better performance for a mass ratio around 7.5%, which loses efficiency for higher mass ratios. However, the efficiency in the control of displacements of the tower top, for the hybrid and active systems, increases with the growth of the mass ratio until it reaches a horizontal asymptotic behavior from a mass ratio around 10%. Nonetheless, for values of mass ratio with this value or higher, the feasibility of using the vibration absorber is hampered, if not impossible.
Another important aspect in this parametric analysis, which is comparative among the proposed control systems, refers to the energy required to achieve a certain level of control. It is possible to observe that the hybrid and active systems have similar behavior, both with respect to the tower response (Figure 6a) and with respect to the vibration absorber response (Figure 6b). In Figure 6a, the control behavior of the hybrid system is better than that of the active system because, on average, the r.m.s. displacement of the tower top is 24.3% lower when the hybrid control is used. In addition, Figure 6b shows that, on average, the ratio between displacement r.m.s. of the absorber and tower top controlled in the hybrid system is 26.7% higher than that of the active system.
Finally, comparing the data of the graph of Figure 7, it is observed that the hybrid system requires a control force of a smaller magnitude in relation to a purely active system, because the curves have practically the same behavior, being only separated vertically, on average, by a value of 1.846 E4 N. The previous comments establish and confirm some advantages of the hybrid control system in relation to the purely active system.
4.5 Transient responses
The harmonic excitation to which the tower is subjected was applied for a time of only 5 s, from which the excitation ceases and the towerDVA assembly vibrates freely (Figure 8 and Figure 9). This excitation simulates an intensity of a tipic wind gust. Table 4 shows the tower top maximum and r.m.s. displacements of uncontrolled and controlled systems; as well as the percentage reduction of such displacements.
Without control  Passive control  Hybrid control  Active control  
Maximum (m)  r.m.s. (m)  Maximum (m)  r.m.s. (m)  Maximum (m)  r.m.s. (m)  Maximum (m)  r.m.s. (m) 
0.56497  0.23126  0.23190  0.07451  0.15757  0.05463  0.19671  0.06559 
Percentage reduction of r.m.s. (%)  
Without control → Passive control  Passive control → Hybrid control  Without control → Active control  
67.78  26.68  71.64 
Figure 8 shows the graph of uncontrolled and passively controlled top of tower displacements, from which we can observe a high reduction of such displacement values; as shown in Table 4. There is a reduction of 67.78 and 71.64% in the r.m.s. of the tower top displacement in relation to the uncontrolled case with the introduction of the passive and active DVA, respectively, in the tower structural system.
For comparisons of passively and hybridly controlled tower top displacements, the graphs of Figure 9a were plotted. It is observed that the hybrid system performs better than the purely passive system, reducing the r.m.s. of the tower top displacements by 26.68% in relation to the displacement of the top of the passively controlled tower. In addition, it is noticed that the vibration stops earlier (Figure 9a) and that the vibration peak is smaller (Table 4) when using the hybrid control system.
It is shown in Figure 9b the graph of the displacements of the vibration absorber for tower equipped with passive and hybrid control. It should be noted that the absorber response becomes more rapid with the introduction of the actuators from the hybrid and active control systems, making these systems more suitable for situations where rapid control action is required for efficiency and robustness of the control system (as in the case of control of seismically excited structures).
Finally, the graph of the control force acting on the dynamic vibration absorber with reaction at the top of the tower is shown in Figure 10. By means of the values of the control forces, one has an idea of the energy required for the control of the tower by the hybrid system. It is important to note that there is a lag (about 1.3 s) between the peaks of the control force graphs (Figure 10) and the vibration absorber displacement (Figure 9b), so that the peak of the control force excites the vibration absorber, and it subsequently moves toward its displacement peak to control the structure of the tower.
4.6 Steadystate responses
In this section, the structure of the towerDVA submitted to stacionary resonant harmonic excitation is analyzed, for which an analysis time of 100 s is used to calculate the r.m.s. values. Table 5 shows the tower top maximum and r.m.s. displacements of uncontrolled, passive and active controlled system.
In Figure 11 is shown the graph of the uncontrolled and passively controlled top of the tower, from which we can observe the high efficiency of the passive system in the control of the two 1st vibration modes (XY and YZ) of the tower, in view that the vibration absorber has been tuned to such frequencies. According to the data of Table 5, with the use of the passive and active systems, it is observed that there is a reduction of 93.87 and 95.26% in the r.m.s. of the tower top displacement in relation to the case without control, respectively.
Without control  Passive control  Hybrid control  Active control  
Maximum (m)  r.m.s. (m)  Maximum (m)  r.m.s. (m)  Maximum (m)  r.m.s. (m)  Maximum (m)  r.m.s. (m)  
6.02837  2.91768  0.25723  0.17884  0.15757  0.10828  0.19768  0.13833  
Percentage reduction of r.m.s. (%)  
Without control → Passive control  Passive control → Hybrid control  Without control → Active control  
93.87  39.46  95.26 
For the comparison of the passive and hybrid control systems, the graph of passively controlled and hybridly controlled tower displacements (Figure 12a) and Table 5 data are analyzed; in which it is possible to improve the performance of the control system, since the r.m.s. value of tower top displacements is reduced in 39.46% compared to the tower top displacements controlled passively.
In this way, the introduction of active control improves the robustness of the system in situations that require to refine/improve the vibration control. On the other hand, the passive control considerably reduces the displacements, so that the control system works well in situations of loss of electric power and, consequently, impossibility of actuation of the active system. However, when the structure is subjected to a detuned excitation in relation to the prevailing vibration frequency of the tower, the active control actuators begin to correct or improve/refine the behavior of the purely passive control, since the active control also works well for a frequency of excitation not tuned at any of the modal frequencies of the uncontrolled tower structure. Such performance of hybrid and active systems may become, in certain cases of practical application, indispensable for the proper functioning and safety of the tower structure and its users.
Corroborating with the above findings, in Figure 12b is shown the graph of the displacements of the vibration absorber for tower equipped with passive and hybrid control, in which the most effective action of the hybrid control is observed, since the dynamic vibration absorber responds more rapidly and significantly to the action imposed on the structure of the towerDVA assembly. Finally, in Figure 13 is shown the graph of the control force exerted by the actuators that act on the dynamic vibration absorber and have a reaction at the top of the tower.
The values of the control forces provide a perception of the energy required for the control of the tower by the hybrid system. In addition, the maximum control force (between 20.0 and 30.0 kN) is perfectly possible to be applied using an electronically controlled hydraulic pump with a working pressure between 1000 and 1500 bar, coupled to a piston with a crosssectional area of 200 cm^{2} that highlights the technical feasibility of the proposed control system.
4.7 DVA constructive model
In this item, a sketch of the control mechanism developed in this study and based on the device proposed by ^{Accioly (2006}) is presented, according Figure 14. Then, the main components of this mechanism are described and enumerated, aiming to show the viability of design and execution of the vibration control system in tubular towers for wind turbines developed in this paper.
1  steel housing for control mass support: external dimensions 1,00 x 1,00 x 2,25 m;
2  lead mass control that fills the steel housing (nº 1);
3, 4, 5, 6, 7, 8  viscous dampers (nº 3 e 6), hydraulic actuators (nº 4 e 7) e linear springs (nº 5 e 8) positioned in direction X;
9, 10, 11, 12, 13, 14  viscous dampers (nº 9 e 12), hydraulic actuators (nº 10 e 13) e linear springs (nº 11 e 14) positioned in direction Z;
15a,b  hinges: six hinges (nº 15a: free rotation around the Z axis) linking the components nº 3, 4, 5, 6, 7, and 8 to the element nº 24; and other six hinges (nº 15b: free rotation around the X axis) joining the components nº 9, 10, 11, 12, 13, and 14 to the element nº 24;
16a,b  hinges: six hinges (nº 16a: free rotation around the Y axis) linking the components nº 3, 4, and 5 to the element nº 18, and the components nº 6, 7, and 8 to the element nº 19. And other six hinges (nº 16b: free rotation around the Y axis) joining the components nº 9, 10, and 11 to the element nº 21, and the components nº 12, 13 e 14 to the element nº 22;
17, 18, 19  linear slide guide positioned along the Z axis: allows only the translation movement of the runners (nº 18 e nº 19) in the longitudinal direction of the rail nº 17. The rail is fixed to the component nº 23, that, in turn, is fixed to the support housing (nº 1) to the control mass (nº 2);
20, 21, 22  linear slide guide positioned along the X axis: allows only the translation movement of the runners (nº 21 e nº 22) in the longitudinal direction of the rail nº 20. The rail is fixed to the component nº 23, that, in turn, is fixed to the support housing (nº 1) to the control mass (nº 2);
23  structure for fixing rails: it is a reinforcement structure composed of metal bars with cold Ushaped sections, which serves as a support for the rails of linear slide guides (nº 17 e nº 20). Then, this structure is attached to the control mass (n° 2) support housing (nº 1);
24  tubular structure of the tower: where the hinges nº 15a,b are fixed.
25  king pin: metal shaft supporting the control mass (nº 2). The king pin is connected to the control mass by means of the component nº 26;
26  king pin fixing system: set of metal plates fixed to the kingpin (nº 25) and to the control mass (nº 2);
27  combined bearing: it provides rotation about the king pin's Y axis (nº 25) and the control mass (nº 2). Such a bearing is capable of withstanding axial and radial stresses. Although the control device designed in this work intends to control translations in the X and Z directions, the construction of the vibration absorber with this type of bearing allows the use of the same to the control of torsional rotations (rotation about the Y axis);
28  metal cylinder: which is attached to the component nº 29, it has the function of behaving the combined bearing (nº 27);
29 metallic base: it is manufactured using tubes of rectangular crosssection and has the function of joining the metallic cylinder (nº 28) to the component nº 30;
30  movable upper metallic base: it is manufactured using tubes of rectangular section. At its upper part, it is fixed to the metallic base (nº 29) and, in the bottom part, is fixed to the components nº 31a,b,c,d;
31a,b,c,d; 32a,b  linear slide guide positioned along the Z axis: only allows translation movement of the upper metallic base (nº 30) on the direction of the rails nº 32a,b. The latter are fixed to the component nº 33. The components 31a,b,c,d represent the stakes that are fixed to the element nº 30;
33  movable bottom metallic base: it is manufactured using tubes of rectangular section. At its upper part, it is fixed to the metallic base (nº 29) and, in the bottom part, is fixed to the components nº 31a,b,c,d;
34a,b,c,d; 35a,b  linear slide guide positioned along the X axis only allows translation movement of the bottom metallic base (nº 33)on the direction of the rails nº 35a,b. The latter are fixed to the component nº 36. The components 34a,b,c,d represent the runners that are fixed to the element nº 33;
36  rectangular tubular beam: on it are attached the rails nº 35a,b;
37  counterbalance: formed by a lead mass used to balance the control mass (nº 2), keeping the center of mass thereof coinciding with the center of the tubular tower mass;
38  stopblock with rubber: safety device for limiting the movements of the DVA and for absorbing the impact of the runners nº 31a,b,c,d on the rails nº 32a,b and of the runners nº 34a,b,c,d on the rails nº 35a,b. That is, it serves in the case of an unforeseen loading situation that would lead to the displacement of the runners beyond the limits of the rails.
5 CONCLUSIONS
After analysis developed in this paper, it is observed that the vibration absorbers proposed, both passive and hybrid, presented efficiency in the tower control vibration when the towerabsorber assembly is submitted to a harmonic action resonant to the 1st tower mode of vibration. Additionally, through the constructive scheme of the absorber exposed, the executive viability of the vibration control system for widescale horizontal axis wind turbine tubular towers is confirmed.
The results of the modal analysis performed for both the beam and shell finite element model are similar, especially in relation to the 1st vibration mode of the tower. Therefore, the model represented by beam FE can be used as representative of the tower dynamic behavior when it is submitted to the resonant excitations to its fundamental frequency. With this, the action of the floating part of the wind was modeled by resonant harmonic forces, considering the tendency that there is the excitation of the fundamental way of the tower. Since the increasing intensity of the wind forces with height relative to ground level, the increase of the horizontal displacements to the tower with respect to the height in the 1st mode of vibration and the mass concentrated in the top (nacelle), of the same order of magnitude of the tower mass, are factors that contribute to excitation of the 1st mode of vibration of the tower.
Three types of absorbers were then comparatively studied: tuned mass damper (TMD), active mass damper (AMD) and hybrid mass damper (HMD). The HMD, the main contribution of this paper, reached excellent levels of vibration reduction for the tower subjected to harmonic actions, in a transient and permanent (stationary) regime, resonating with the first vibration mode of the structure without absorber.
In view of the obtained results, it is verified that the purely passive system is quite efficient in the control of vibrations of the excited tower according to the 1st vibration mode for which the absorber is tuned. Hybrid and active systems have been able to significantly improve passive system performance. It can be noted that the absorber response becomes faster with the introduction of hydraulic actuators of the hybrid control system, making this system suitable for situations where rapid control action is required for effectiveness and robustness of the system.
It is concluded, therefore, that the results of this work involve contributions of immediate practical interest, since the research develops subsidies for structural analysis and vibration control of the towers for multimegawatt wind turbines to be implanted in Brazil considering the recent growth of the country's wind potential.