Research on Vibration of WindPower Tower Based on Response Spectrum Analysis

Mingliang YANG*, Wenjun MENG**, Zhengyan CHANG***, ShouchengWANG****, Peng WANG***** *School of Machinery and Electronics Engineering, Taiyuan University of Science and Technology, Taiyuan, Shanxi, 030024, China, E-mail: yangmingliang1997@163.com **School of Machinery and Electronics Engineering, Taiyuan University of Science and Technology, Taiyuan, Shanxi, 030024, China, E-mail: tyustmwj@126.com ***School of Machinery and Electronics Engineering, Taiyuan University of Science and Technology, Taiyuan, Shanxi, 030024, China, E-mail: 535061752@qq.com ****Taiyuan Heavy Industry Co., Ltd, Taiyuan, Shanxi, 030000, China, E-mail: 13603535316@139.com *****School of Machinery and Electronics Engineering, Taiyuan University of Science and Technology, Taiyuan, Shanxi, 030024, China, E-mail: 595855655 @qq.com


Introduction
The wind-power tower is a flexible structure the characteristics of 'high, large, long and thin'.The width and depth of the tower is much smaller than its height.The wind load is one of the main design loads, which is divided into two parts, including static wind load and fluctuating wind load.The latter can cause structural vibration [1].For the calculation of wind-induced vibration, the peak response of the structure considered as a linear structure is concerned and the spectral analysis method is used to solve the problem [2][3][4].Based on the analysis of wind load response, the response of the structure with dampers or without dampers is analyzed by the spectrum analysis method.

Transverse wind load response
Under the action of horizontal wind, the tower structure with circular cross section will only cause the lateral wind response.When the airflow rounds the tower structure, a vortex will be formed due to the airflow meeting again behind the structure.So, the maximum acceleration of the structure may occur in the crosswind.
After the air flows tower, flow characteristics at different stages are divided into three ranges according to the Reynolds number [1], as shown in Fig. 1.
Subcritical range is and the random vibration will be produced for structure, irregularly.
Trans critical range is Experiments show that when the vortex shedding frequency reaches the natural frequencies of the structure owing to the increase of wind speed, the frequency is constant in the wind speed region, which is called the locked area.Especially when the vortex shedding frequency is consistent with the natural frequency of the structure, a large resonance response will be generated.
Reynolds number formula is expressed as: .
Where: v is the wind speed.l is the maximum size of the cross-section of wind-power tower. is the dynamic viscosity.For the air, it And etc. From Eq. ( 2), it can be seen that the Reynolds number is proportional to the wind speed.Therefore, the key of the wind power structure design is to check whether the resonance wind speed is a critical cross range or not.Especially when the frequency of vortex shedding is consistent with the structure frequency, it will produce a resonant response of several times larger than that of the static force [1].For the structure in the sub critical range, although it can also occur resonance, the wind speed is smaller.For the structure function being not as serious as the cross critical range, the method of increasing structural stiffness and strength is usually adopted to deal with it.
The vibration period of the wind-power tower is between 0.5 and 2.5s, and the resonance speed is from 40 m/s to 24 m/s.For

Re
, the wind speed can be in the extreme weather condition.

Vibration control of wind-power tower structure
The purpose of vibration control of wind-power tower is to reduce the displacement and acceleration response of the structure [5-10].

Mechanical model of damper
A high-rise structure in the j mass of the structure is installed with a tuned mass damper (TMD).For the structure under the action of the damper, the calculation diagram of the damper and tower structure is given in Fig. 2. The motion equation of the tower structure under the TMD action can be designed as: where: [M], [C] and [K] are the mass, damping and stiffness matrixes of the tower structure, respectively.      are the acceleration, velocity and displacement vectors of the tower structure with respect to the ground, respectively.    t , z f is the wind load vector.  TMD f is the damping load vector.

TMD motion equation
The TMD motion equation is defined as: where: md,cd and kd are the mass, damping and spring stiffness of TMD, respectively.xj s the displacement of the j mass of the structure with respect to the ground.xdis the displacement of TMD relative to the ground.

Mean square value of relative displacement of structure at any point and TMD
After installing TMD, if only the first vibration mode is considered, Eq. ( 6) can be obtained by the Eq. ( 3)and where: v is the displacement of TMD with respect to the j mass, It can be obtained by Eq. ( 5).So, .
The formula (8) and type (9) are written in matrix form: Eq. ( 11) can be simplified as: Frequency response function matrix can be introduced as: Where: The spectral density function of stochastic generalized coordinates is defined as: When the TMD is installed, the mean square value of the displacement of the structure at any point m can be designed as: After the installation of TMD, with no TMD, of the ratio of the displacement mean square value of the structure at any point m is given as: Similarly, the power spectral density of the relative displacement of the TMD mass can be obtained as: The mean square value of the relative displacement of TMD can be expressed as:

Structural damper
The damping device can reduce the effect of the wind vibration response, which is caused by the shortening of the system period due to the additional stiffness of dampers and the reduction of the energy generated by the damping characteristic of the damper [11].
In order to analyze the influence of the parameters of TMD on the wind load of tower structure, for tower structure with and without TMD, the displacement formula of the wind load is expressed by the parameters of TMD, mass, damping and stiffness.

Frequency response function
According to the integral formula of the complex frequency response function, Where: 1 ; From Eq. ( 25), it can be seen that when the structure of TMD is set up, the damping of the tower structure will be increased.The greater the damping, the smaller the response is.The equivalent damping is related to the mass

Parameters selection
The frequency modulation quality controller is composed of four parts: mass block, spring system, damping system and mass block supporting system, as shown in Fig. 3.

Fig. 3 Pendant frequency modulation mass damper of tower tube of wind generator
In the vibration control.The inertia force generated by the mass block and the relative motion of the structure is mainly used to achieve the purpose of reducing vibration.The ratio of the mass and the tower structure quality is between 0.02 and 0.005 [4].
The spring system is composed of 4 helical springs.The function is to adjust the frequency of the vibration control, so that it is close to the controlled natural frequency of the structure.
The damping system provides damping for the whole vibration control to achieve the good vibration absorbing effect.The appropriate damping coefficient is chosen from .The mass support system is a hanging hook and the friction between the hook and the mass is very small, which can ensure free sloshing of the mass block.

Response spectrum analysis
The response spectrum analysis is mainly used to determine the dynamic response of wind-power tower structure to random wind load.The result of tower modal analysis is connected with the wind speed spectrum to calculate the displacement and stress of the model.

Response spectrum
Random vibration of wind power tower under dynamic wind load will be produced and the stochastic process simulation should be used in dynamic response analysis.In this paper, the linear filtering method combined with Matlab software is used to simulate the fluctuating wind speed time history of the wind-power tower under the wind load.
Wind speed spectrum for structural analysis davenport spectrum [12].1200 Where: n is the frequency of wind speed in the wind, Hz. x is the turbulence integral scale coefficient.For the wind speed being 11.4m/s and wind-power generator being in the standard condition, the Reynolds number Re is 2.934×106.The structure is not conducted the check calculation in the supercritical range, because the tower wind vortex shedding wind vibration is not formed, obviously.
For the wind speed being 25m/s and wind-power generator being in the cutting condition, the Reynolds number Re is 6.43×106.At this time, the structure may occur the resonance exceeding the critical range.So, the resonance response should be calculated.
For the wind speed being 50m/s and wind turbines being in extreme condition, the Reynolds number Re is 12.87×106.At this time, the structure exceeds the critical range.So, the resonance response should be determined.Using the linear filtering method, under the three conditions including standard operating condition, cutting out condition and limit condition, the fluctuating wind speed time history of wind-power tower structure bearing the wind load is simulated by Matlab.For the limit of speed being 50m/s, the wind velocity spectrum is as shown in Fig. 4 and Fig. 5.
From the above chart, the power spectral density function curve of the numerical simulation is in good agreement with the target spectral density curve.So, the simulation curve of fluctuating wind speed-time history of wheel center has a very high accuracy.The wind speed spectrum curve can be adopted to conduct the finite element wind vibration response analysis to wind-power tower structure.For the wind speed being 50m/s, the displacement response of the structure is presented in Fig. 6.
It can be seen from Table 2 that the maximum displacement and acceleration of structures under extreme climatic conditions are greatly reduced for providing the suspension tuned mass damper.

Conclusion
In this paper, the response of wind power tower structure is studied by the spectral analysis method under three different conditions and the vibration control equipment is used to make the comparison of structures.There are the following main conclusions: 1) The AR method of the linear filtering method is used to simulate the wind speed spectrum of tower structure under three different conditions and the fluctuating wind speed simulation spectrum is always floating around the target spectrum curve and almost symmetrical, which indicate that wind speed meets expectations and has a fairly high accuracy.So, it can be used to analyze the wind-induced vibration response of the wind-power tower.
2) Three dimensional finite element analysis model of wind power tower is built, and the random vibration analysis of the structure is carried out.From the analysis results, for the wind speed being 11.4m/s, the displacement and acceleration of the structure in the supercritical range is smaller and the tower wind vortex shedding wind vibration is not formed, obviously.For the wind speed being 25m/s and wind-power generator being in the cutting condition, the structure may occur the resonance exceeding the critical range.So, the resonance response should be calculated.For the wind speed being 50m/s and wind turbines being in extreme condition, the structure exceeds the critical range.At this time, the structure displacement and acceleration are larger, but the value is less than the cut off working condition.
3) The maximum displacement and acceleration of the structure are greatly reduced under three different wind loads, owing to the suspension tuned mass damper be set.As a flexible structure with the height-to-width ratio being quite large, the wind-power tower structure is easily produce the wind-induced vibration under the action of dynamic wind load.First of all, the probabilistic characteristics of the response of structure and its control system in the frequency domain are analyzed in theory.Secondly, the structural responses of wind load under three conditions, including standard operating condition, cutting condition and limiting condition, are analyzed by the spectral analysis method in practice and the control effect of the damper on the structure is discussed.The results show that the tower structure has different response characteristics under different wind speeds, and the maximum displacement and acceleration of the structure are greatly reduced after the suspension tuned mass damper is set up.

Fig. 1
Fig. 1 Schematic diagram of transverse wind partition vibration will be occurred periodically.

Fig. 2
Fig. 2 Mechanical model of wind turbine tower 3.1.1.Motion equation of tower structure under TMD action

. 1 M 1  1 
 is the ration of TMD mass and structure generalis the generalized mass. is the damping ratio of first modal for the structure.is the first natural frequency of structure.j  is the value of First mode at the point j.
and it is also related to the position of TMD in the structure.

10 v
is the average wind speed at 10 m height from the ground.k is the surface drag coefficient，Corresponding to the values of A, B, C, and D. Four type landforms for k are shown in

Fig. 4 Fig. 5
Fig. 4 Wind conditions the simulated wind speed spectrum of wheel center under 50m/s wind speed

Fig. 6
Fig. 6 Displacement of the tower structure under the 50m/s condition

Table 1
Four type landforms for k 1.

Table 2 .
Table 2 Tower structure response under limit conditions