Research Article  Open Access
Bacem Zghal, Imen Bel Mabrouk, Lassâad Walha, Kamel Abboudi, Mohamed Haddar, "Analyses of Dynamic Behavior of Vertical Axis Wind Turbine in Transient Regime", Advances in Acoustics and Vibration, vol. 2019, Article ID 7015262, 9 pages, 2019. https://doi.org/10.1155/2019/7015262
Analyses of Dynamic Behavior of Vertical Axis Wind Turbine in Transient Regime
Abstract
In this paper, the dynamic behavior of a onestage bevel gear used in vertical axis wind turbine in transient regime is investigated. Linear dynamic model is simulated by fourteen degrees of freedom. Gear excitation is induced by external and internal sources which are, respectively, the aerodynamic torque caused by the fluctuation of input wind speed in transient regime and the variation of gear mesh stiffness. In this study, the differential equations governing the system motion are solved using an implicit Newmark algorithm. In fact, there are some design parameters, which influence the performance of vertical axis wind turbine. In order to get the appropriate aerodynamic torque, the effect of each parameter is studied in this work. It was found that the rotational speed of the rotor shaft has a significant effect on the aerodynamic torque performance.
1. Introduction
Generally, vertical axis wind turbines (VAWT) have a particular architecture compared with horizontal ones. They are composed of two main parts: the blade rotor in vertical position and a mechanical gear transmission (bevel gear). Vibrations of the aerodynamic part caused by the wind speed variation are transmitted to the other part (gear transmission system) via shaft, gears, and bearing.
In literature, plenty of authors studied the aerodynamic performance of Darrieustype of VAWT. There are two mainly approaches: momentum models and Computational Fluid Dynamics (CFD). The main benefit of momentum models is that their time of resolution is quicker than the other approach. Although the Computational Fluid Dynamics have been a useful design tool for studying the efficiency of wind turbine, the mesh generation in threedimensional analyses needs a lot of time for the simulation. Among analytical models are researches [1–3] based on Actuator Disk and BladeElement Method (BEM) to predict the aerodynamic torque of VAWT. In addition, computational aerodynamics methods such as multiple stream tube method or vortex model [4–7] have been developed to optimize VAWTs performance.
The (CFD) method has been widely used in developing the characteristics of wind turbine (torque fluctuation, power output, and pressure distribution). Jiang et al. [8] employed a commercial CFD simulation for studying the effect of geometrical parameters and airfoil type on the performance of the HDarrieus turbine with fixed pitch angle. Also, a numerical analysis of Hrotor Darrieus turbine is introduced by M.H Mohamed et al. [9]. The developing of torque fluctuation is a very important step for the reason that the instantaneous torque produced by the rotor blade is directly related to both the power generation and the gearbox vibration.
Accordingly, there are many related literatures studying bevel gears transmission. CaiWan Changjian [10] studied the nonlinear dynamic behavior of bevel gear system. Besides, Fujii et al. [11] analyzed the dynamic vibration of straight bevel gear supported by angular bearings and tapered roller. M. Li and H. Y. Hu [12] studied the dynamic analysis of a spiral bevelgeared rotorbearing system. Y. Wang et al. [13] and J.F. Besseling [14] developed a new approach based on finite element theory to model bevel gear systems. Moreover, Driss Yassine et al. [15] present the model of twostage straight bevel gear system excited with only internal excitation which is the periodic fluctuations of the gear meshes’ stiffness.
In this paper, we discuss the impact of some design parameters including number of blades, turbine radius, chord length, blade length, and rotational speed on the aerodynamic torque of the HDarrieus VAWT through analytical approach.
The main objective of the present work is to predict the dynamic behavior of the onestage straight bevel gear system used in vertical axis wind turbine and powered by two main sources of excitation which are the optimum aerodynamic torque selected through parametrical study and the periodic variation of the gear meshes’ stiffness.
2. Theoretical Modeling of VAWT Rotors
In this section, analytical investigation of aerodynamic torque of Darrieus wind turbine is established. The actuator disk theory is chosen for the aerodynamic study of the Darrieustype wind turbine with straight blade. This theory characterizes the turbine as a disc with a discontinuity of pressure in the stream tube of air, which causes a deceleration of the wind speed. Referring to Figure 1, the relative flow velocity can be obtained as follows:The angle of attack is defined as the angle between the resultant air velocity vector and the blade chord. It is expressed as follows:The value of the axial induction factor (a) can be introduced by actuator disk theory. The resultant air velocity is dependent on the induced velocity and the tip speed ratio (TSR) defined asIn this work, the induced velocity is modeled in deterministic form, as a sum of several harmonics [17]:The resulting aerodynamic forces in the blade can be founded by the interpolation of the lift and drag coefficients relative to the symmetrical airfoil used (NACA0018), the angle of attack, and the given Reynolds number.
The tangential and normal forces as function of the azimuth angle θ can be calculated using the bladeelement theory [2].where and are the normal and tangential coefficients, respectively, calculated from the lift and drag coefficients ( and ) using the same theory (bladeelement momentum theory), given by following expressions:The average torque produced by the rotor (n blades) is generated from the average tangential force acting on one blade:The average torque coefficient is calculated byFinally, the power coefficient Cp is estimated from the average torque coefficient:
3. Theoretical Modeling of the OneStage Bevel Gear System
This part investigates the studying of the dynamic behavior of bevel gear system used in vertical axis wind turbine. The main excitations of the onestage bevel gear system are the selected aerodynamic torque estimated through parametrical analysis in addition to the internal mesh stiffness excitation.
The dynamic model of singlestage bevel gear is presented by fourteen degrees of freedom (see Figure 2). This model includes two blocks where the first block is constituted by the Darrieus rotor modeled by the mass linked to the wheel via a first shaft with torsional rigidity ; this block is supported by a bearing.
The second block includes the pinion , the second shaft with torsional rigidity , and the generator modeled by a mass . It is also supported by a bearing. The wheel is connected to the pinion via teeth mesh stiffness. The bearings are modeled by linear springs acting on the lines of action. The mesh stiffness characterizes the elastic deformations managing the relative positions of the two gear wheels; it can modeled by the mesh stiffness .
The proposed dynamic model is modeled by the generalized coordinate vector : are the linear displacements of bearing in each block.
are the angular displacement of bearing and wheel (i=1:2).
The equations of motion describing the dynamic behavior of the model are established using the formalism of Lagrange for each degree of freedom of the system.The total mass matrix is defined by is the average stiffness matrix of the structureThe mesh stiffness matrix can be defined by[k(t)] is the total meshing stiffness of the gear pair. Sha Wei et al. [18] have expressed the total meshing stiffness of the gear pair by a periodic excitation decomposed on Fourier series.
The tooth deflection following the line of action is defined byThe components of the tooth deflection are presented in Table 1.

The load vector can be written as and are the aerodynamic torque produced by the rotor and the resisting torque, respectively. is the proportional damping matrixwhere , , , and .
and r_{21} are the halfangle of the base circle of bevel gear (i) of the block (j) and the radius of the sphere which contains two bevel gears, and , respectively (r_{12} = r_{21}).
The angles and u describe the different parameters of the bevel gear [18].
4. Influence of Rotor Configuration on Turbine Performance
The main purpose of this study is to investigate the effect of different design parameters (blade chord, number of blades, and radius turbine) on the aerodynamic torque evaluation of HDarrieus turbine, using analytical approach.
The solidity is defined as an important nondimensional parameter, which influences the selfstarting abilities and determines the applicability of the momentum models. The turbine is able to selfstart for high solidity ( ≥ 0.4) [2]. For straight bladed VAWTs, the solidity is calculated byIn order to study the dynamic behavior of bevel gear system used on vertical axis wind turbine and powered by the aerodynamic torque in transient regime in addition to the periodic variations of mesh stiffness, we used numerical simulation based on implicit method of Newmark.
The parameters of the bevel gear transmission are presented in Table 2. Specifications of the vertical axis wind turbine are shown in Table 3.


The number of blades (n) is an important factor that influences the torque produced. It is well known that a bigger number of blades give rise to the solidity and produce a global torque with small fluctuation. However, increasing the number of blades lead to increase in the turbine drag by increasing the number of connecting shafts. Figure 3 shows the effect of the blade number on the torque in full revolution for a case where the radius and blade chord are maintained constant while the number of blades changes. As can be deduced, 3bladed VAWT perform more efficient than 2 and 4 bladed turbines.
In Figure 4 the effect of blade chord on the torque fluctuation is presented. We remark that torque magnitude increases when the chord length increases. Also, blades with smaller chords require a bigger tip speed ratio to achieve a maximum torque so the selection of chord length affects the selfstarting behavior of Darrieus turbine.
As can be assumed from the torque equation (8), the rotor radius and the blade length have a major contribution in the torque behavior of VAWT turbine. An increase in radius turbine and blade height advances the instantaneous torque as shown in Figures 5 and 6.
Figure 7 shows the total torque evolution versus azimuth angle generated for a complete revolution with different rotation speed. It can be observed that the torque increases considerably with the increase of the rotation speed. Also, the torque fluctuation becomes positive when the rotation speed reaches 178 rev/min. However, for a case of 50 rev/min and 177 rev/min the instantaneous torque presents some fluctuation with negative magnitude which can explain the no selfstarting of NACA0018 airfoil.
As shown in Figures 8 and 9 rising wind speed from 8 to 14 m/s significantly affects the torque which increases considerably. The cyclic distribution of total torque at higher rotational speed (N=180 rev/min) and at small rotational speed (N=50 rev/min) is observed at the same wind speed of 14 m/s.
We clearly see only positive fluctuation of the aerodynamic torque at N=180 tr/min; however, this torque presents negative oscillation at small rotational speed (N=50 tr/min) at fixed wind speed of 14 m/s. Consequently, we conclude the most favorable speed excitation of the considerable Darrieus rotor that corresponds to wind velocity of 14m/s and rotational speed of 180 tr/min which respects the condition of positive torque evolution.
5. Conclusion
This paper presents a threedimensional model of onestage straight bevel gear system used in vertical axis wind turbine. Aerodynamic torque fluctuation and periodic oscillation of the gear meshes’ stiffness are the main sources of excitation for the bevel gear system.
Influence of geometrical parameters of Darrieus rotor has been done in order to select the appropriate parameters. The optimization process has been carried out on the effect of each parameter on the aerodynamic torque produced.
It was found that the aerodynamic torque increases when the chord, the radius, and the height of VAWT rise. However, the best performance is detected for 3bladed VAWT. Finally, the most significant parameter that affects the aerodynamic torque is the rotation speed of the rotor shaft.
Nomenclature
a:  Axial induction factor 
A:  Turbine swept area 
c:  Chord (m) 
:  Proportional damping matrix 
:  Tangential force coefficient 
:  Normal force coefficient 
:  Average torque coefficient 
Cp:  Power coefficient 
, :  Lift and drag coefficients 
:  Tangential force 
:  Normal force 
:  External excitation force 
h:  Blade height (m) 
:  Moment of inertia of the wheel of block i 
, :  Diametrical moment of wheel of block respectively following the X and Yaxes 
:  Gear mesh of stiffness 
:  Time stiffness matrix of the gear mesh fluctuation 
:  Stiffness matrix of the average structure 
:  Shafts torsional stiffness 
, , , and :  Bending and tractioncompression bearing stiffness 
:  Mass matrix 
:  Mass of block i 
:  Vector of geometric parameters of the dynamic model 
n:  Blade number 
N:  Rotational speed of the rotor (tr/min) 
:  Generalized coordinate vector 
R:  Radius of the wind turbine (m) 
:  Resisting torque 
:  Aerodynamic torque 
:  Wind free stream velocity (m/sec) 
Va:  Induced velocity 
Vc, Vn:  Chordal velocity component and the normal velocity component, respectively 
W:  Relative flow velocity 
:  Azimuth angle 
:  Angular velocity (rad/sec) 
:  Air density[kg/m^{3}] 
:  Tip speed ratio 
:  Relative displacement of the contact point along the line of action 
, :  Geometric angle of bevel gear 
:  Bearing displacements in each blocs (i=1:2) 
:  Angular displacement of the bearing following X and Y, respectively 
:  Angular displacement of wheel and gear following Z direction (i=1:2) 
:  Solidity 
:  Angle of attack. 
Data Availability
All results are obtained by simulation with Matlab.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 R. E. Wilson, “Windturbine aerodynamic,” Journal of Industrial Aerodynamics, pp. 357–372, 1980. View at: Google Scholar
 C. Javier, Smallscale vertical axis wind turbine design. Bachelors Thesis, Aeronautical Engineering, Tampere University of Applied Sciences, 2011.
 D. Prathamesh and C. L. Xian, “Numerical study of giromilltype wind turbines with symmetrical and nonsymmetrical airfoils,” Journal of Science and Technology, 2013. View at: Google Scholar
 J. H. Strickland, “The Darrieus turbine: a performance prediction model using multiple streamtube,” Sandia Laboratories Report. SAND, pp. 75–431, 1975, United States of America. View at: Google Scholar
 I. paraschivoiu, “Doublemultiple stream tube model for studying verticalaxis wind turbines,” Journal of Propulsion and Power, vol. 4, no. 4, pp. 370–377, 1988. View at: Publisher Site  Google Scholar
 H. Beri and Y. Yao, “Double multiple stream tube model and numerical analysis of vertical axis wind turbine,” International Journal of Energy and Power Engineering, vol. 3, no. 3, pp. 262–270, 2011. View at: Publisher Site  Google Scholar
 L. Wang, L. Zhang, and N. Zeng, “A potential flow 2D vortex panel model: applications to vertical axis straight blade tidal turbine,” Energy Conversion and Management, vol. 48, no. 2, pp. 454–461, 2007. View at: Publisher Site  Google Scholar
 Z.C. Jiang, Y. Doi, and S.Y. Zhang, “Numerical investigation on the flow and power of smallsized multibladed straight Darrieus wind turbine,” Journal of Zhejiang University SCIENCE A, vol. 8, no. 9, pp. 1414–1421, 2007. View at: Publisher Site  Google Scholar
 M. H. Mohamed, A. M. Ali, and A. A. Hafiz, “CFD analysis for Hrotor Darrieus turbine as a low speed wind energy converter,” Engineering Science and Technology, an International Journal, vol. 18, no. 1, pp. 1–13, 2015. View at: Publisher Site  Google Scholar
 C.W. ChangJian, “Nonlinear dynamic analysis for bevelgear system under nonlinear suspensionbifurcation and chaos,” Applied Mathematical Modelling, vol. 35, no. 7, pp. 3225–3237, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 M. Fujii, Y. Nagasaki, and M. Nohara Trauchi, “Effect of bearing on dynamic behavior of straight bevel gear train,” Transactions of the Japan Society of Mechanical Engineers A, vol. 61, pp. 234–238, 1995. View at: Google Scholar
 M. Li and H. Y. Hu, “Dynamic analysis of a spiral bevelgeared rotorbearing system,” Journal of Sound and Vibration, vol. 259, no. 3, pp. 605–624, 2003. View at: Publisher Site  Google Scholar
 Y. Wang, W. Zhang, and H. Cheung, “A finite element approach to dynamic modeling of flexible spatial compound bar–gear systems,” Mechanism and Machine Theory, vol. 36, no. 4, pp. 469–487, 2001. View at: Publisher Site  Google Scholar
 J. F. Besseling, “Derivatives of deformation parameters for bar elements and their use in bucking and postbucking analysis,” Computer Methods in Applied Mechanics and Engineering, vol. 12, pp. 97–124, 1977. View at: Google Scholar
 D. Yassine, H. Ahmed, W. Lassaad, and H. Mohamed, “Effects of gear mesh fluctuation and defaults on the dynamic behavior of twostage straight bevel system,” Mechanism and Machine Theory, vol. 82, pp. 71–86, 2014. View at: Publisher Site  Google Scholar
 M. Islam, D. Ting, and A. Fartaj, “Aerodynamic models for Darrieustype straightbladed vertical axis wind turbines,” Renewable & Sustainable Energy Reviews, vol. 12, no. 4, pp. 1087–1109, 2008. View at: Publisher Site  Google Scholar
 A. Mirecki, Comparative study of energy conversion system dedicated to a small wind turbine [Ph.D. thesis], Polytechnic National Institute, Toulouse, France, 2005.
 S. Wei, J. Zhao, Q. Han, and F. Chu, “Dynamic response analysis on torsional vibrations of wind turbine geared transmission system with uncertainty,” Journal of Renewable Energy, vol. 78, pp. 60–67, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Bacem Zghal et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.