Abstract

In this paper, the dynamic behavior of a one-stage 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 Darrieus-type 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 three-dimensional analyses needs a lot of time for the simulation. Among analytical models are researches [13] based on Actuator Disk and Blade-Element Method (BEM) to predict the aerodynamic torque of VAWT. In addition, computational aerodynamics methods such as multiple stream tube method or vortex model [47] 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 H-Darrieus turbine with fixed pitch angle. Also, a numerical analysis of H-rotor 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. Cai-Wan Chang-jian [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 bevel-geared rotor-bearing 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 two-stage 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 H-Darrieus VAWT through analytical approach.

The main objective of the present work is to predict the dynamic behavior of the one-stage 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 Darrieus-type 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 blade-element theory [2].where and are the normal and tangential coefficients, respectively, calculated from the lift and drag coefficients ( and ) using the same theory (blade-element 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 One-Stage 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 one-stage bevel gear system are the selected aerodynamic torque estimated through parametrical analysis in addition to the internal mesh stiffness excitation.

The dynamic model of single-stage 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 r21 are the half-angle 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 (r12 = r21).

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 H-Darrieus turbine, using analytical approach.

The solidity is defined as an important nondimensional parameter, which influences the self-starting abilities and determines the applicability of the momentum models. The turbine is able to self-start 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, 3-bladed 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 self-starting 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 self-starting 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 three-dimensional model of one-stage 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 3-bladed 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 Y-axes
:Gear mesh of stiffness
:Time stiffness matrix of the gear mesh fluctuation
:Stiffness matrix of the average structure
:Shafts torsional stiffness
, , ,   and  :Bending and traction-compression 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/m3]
: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.