Research Article | Open Access
Numerical Investigation of Flapwise-Torsional Vibration Model of a Smart Section Blade with Microtab
This study presents a method to develop an aeroelastic model of a smart section blade equipped with microtab. The model is suitable for potential passive vibration control study of the blade section in classic flutter. Equations of the model are described by the nondimensional flapwise and torsional vibration modes coupled with the aerodynamic model based on the Theodorsen theory and aerodynamic effects of the microtab based on the wind tunnel experimental data. The aeroelastic model is validated using numerical data available in the literature and then utilized to analyze the microtab control capability on flutter instability case and divergence instability case. The effectiveness of the microtab is investigated with the scenarios of different output controllers and actuation deployments for both instability cases. The numerical results show that the microtab can effectively suppress both vibration modes with the appropriate choice of the output feedback controller.
Blade vibration becomes one of the primary concerns as the increased-size blade leads to more flexible and long blades maintaining stiffness and minimizing mass. Due to the aerodynamic effects combined with the blade structural movement, the flutter happens when the wind speed reaches the decreased flutter speed due to increasing structural flexibility of the large blade (e.g., 50 m blades and above) . The typical two degrees of freedom (DOE) of the blade caused by flutter involve torsion (torsional mode) and translation (flapwise mode). There are several factors that can affect the blade flutter instability, such as the centre of mass relative to the centre of the elastic axis in the blade sections, and the ratio between the torsional frequency and the flapwise frequency [2, 3]. The classic flutter of the blade section occurs when the wind speed reaches the flutter speed, leading to the unstable oscillations (flutter instability). With increasing wind velocity, large amplitude of vibration without damping (divergence instability) can happen suddenly. Those flutters will greatly decrease the life of the blade without control. Therefore, there is an increasing interest in the flutter suppression via the smart blade, which is the blade equipped with the appropriate actuator, such as trailing-edge flap and microtabs . The trailing-edge flap has been widely used as the control surface to reject the load and control the aeroelastic instability of the blade [5–7]. For the recent decades, microtabs, proposed by Yen et al. , attract attention as an aerodynamic-effective, cost-efficient, quick-actuation, and light-weight control surface.
The potential of the microtab to regulate blade loads has been numerically and experimentally studied by numerous researchers. Van Dam et al.  first proved the effectiveness of the microtab for load control numerically. Baek et al.  compared the load rejection effect of microtabs with that of flaps via proportional controller. The microtabs were also demonstrated as effective load controlling devices responding to a steady state aerodynamic model [11, 12]. Nakafuji et al.  presented microtab design and wind tunnel tests to verify the use of microtabs for active load control. Maughmer and Bramesfeld  carried out experimental investigation of fixed deployment of microtab as gurney flap. They found that a deployment of the microtab from lower surface increased the amount of lift while the deployment on the upper surface acted in the opposite way. Standish  computed complicated computational studies of tab height and tab location on the pressure side and suction side of the S809 airfoil with the behavior as gurney flap, which is simply a short flat plate attached to the trailing-edge perpendicular to the chord line on the pressure or suction side of the airfoil. He then extended the work to rotorcraft applications . However, most of those works mainly focused on microtab in flow field or microtabs as active aerodynamic load controller. The microtab effect on the passive flutter suppression of the blade has not generally been investigated.
For the purpose of microtab control analysis, a suitable, accurate, and computational-efficient aeroelastic model is needed. For classic flutter study, Theodorsen’s theory has been widely used in modeling the two-dimensional incompressible flow for both fixed wing and rotating wing [17, 18]. Later, for the purpose of analyzing the aeroelasticity using time-marching solution techniques, the aerodynamic behavior of the blade section has been modified in state-space form in time-domain. Different approximations to the Theodorsen theory also have been developed , and Van der Wall and Leishman concluded that, for moderate reduced frequencies, those theories give the same results as the Theodorsen theory. In this paper, one of the approximations is used to provide the aerodynamic model in time-domain for analysis in time-domain instead of frequency domain.
Moreover, the mathematical equations of the blade motions and the consequent computation of the solution are usually complicated and computationally expensive [20–22]. However, for classic flutter of the blade, several studies showed that the dominating coupling of the structural modes with the aerodynamic forces is torsional mode and flapwise mode, and there is rarely coupling of the edgewise mode [23–25]. Therefore, the edgewise mode can be neglected to obtain a more efficient aeroelastic model. Although the blade motions and aerodynamic forces on the blade section have been studied, the aeroelastic model of a smart blade section with the coupling of blade motion, aerodynamic load, and microtab effect has not yet been investigated generally.
This paper demonstrates how the microtab actuator installed on the trailing-edge can affect the aeroelastic stability of the smart blade section in classic flutter. The contribution for this type of model is twofold: it is a computational-efficient aeroelastic model that involves the interaction between aerodynamic forces, structural motion, and microtab actuation for the smart section blade study; the model with the microtabs can be used to evaluate the potential sensors and control strategies for passive vibration control of the blade. The objective of this paper is to gain insight into the control capability of the microtabs on the aeroelasticity of the smart section blade, without complex high degree of freedom modeling or detailed computational-fluid-dynamics modeling. More importantly, this paper offers the preliminary results to design the advanced control algorithms for passive smart blade vibration control. The following sections present the structural model of the smart blade section, describe aerodynamic model applied to the smart section blade, summarize the state-space equations of the aeroelastic model, and analyze the passive control effect of the microtabs on two cases of aeroelastic instability with different choice of control inputs.
2. Structural Model
In this study, the smart blade section with the microtab actuator is modeled based on aerodynamic characteristics of the airfoil UA97W300 and the microtabs installed on the trailing-edge of the blade (90% chord). The smart blade section is shown in Figure 1.
Considering the two-degree-of-freedom (DOE) blade section model with structural linearity and aerodynamic nonlinearity, the differential equations of blade motions of flapwise-torsional oscillations can be written as where is the mass, is the flapwise deflection, is the torsional deflection, is the dimensionless static unbalance about the elastic axis, is the semichord, , are the flapwise and torsional damping coefficients, is the mass moment of inertia about the elastic axis, and , are the flapwise and torsional stiffness coefficient. and are the generalized force and moment acting on the blade section.
The physical generalized forces and moments involve aerodynamic force/moment on the blade, the aerodynamic effects of the microtab actuation, and the supplemental terms regarding the asymmetry of the airfoil. The aerodynamic force and moment are caused by the aerodynamic loads that vary periodically as the blade rotates with constant rotor speed. The microtabs can affect the aerodynamic lift/moment acting on the blade by the deployment on the pressure side or the suction side of the blade with the height (Gurney flap). As illustrated in Figure 1, presents the microtab deployment on the pressure side and indicates the deployment on the suction side. The additional terms describe the static aerodynamic force and moment of the blade when the angle of attack is zero. Since the testing airfoil UA97W300 is an asymmetrical airfoil, the extra terms have nonzero coefficients when symmetric airfoil leads to no extra terms. The generalized force and moment can be expressed aswhere is air density, , are lift and moment coefficients with respect to aerodynamic loads, , are lift and moment coefficients for microtab deployment, and , are lift and moment coefficients at zero angle of attack for the airfoil. The apparent wind speed is .
Define , , , , , and , where , are uncoupled flapwise and torsional natural frequencies and , are the critical damping factors in bending and torsion. Define the dimensionless quantities as , , and , where is the dimensionless wind velocity, is the dimensionless flapwise deflection, and is the ratio between flapwise and torsional natural frequencies. The dimensionless motions of the equations of (1) can be rewritten aswhere where is the dimensionless radius of gyration with respect to elastic axis, is the Mach number, and is the dimensionless mass. The , can be computed using experimental static data of the airfoil at zero angle of attack. , are derived from curve fittings of lift force/moment data with respect to gurney flap (the deployment of the microtab) from wind tunnel test. , are aerodynamic lift/moment that can be obtained from aerodynamic model.
Let the state variables for the dimensionless structural model be , , , , and ; (3) can be rewritten aswhereDefine the control input , ; the state-space equations of the blade model with microtab actuator in (5) can be expressed aswhere
3. Aerodynamic Model
3.1. Aerodynamic Model Based on the Theodorsen Theory
In this study, the aeroelastic stability is affected by the aerodynamic force in linear region for classic flutter research. An aerodynamic model approximated to the Theodorsen theory  is used to compute the unsteady aerodynamic force for classic vibration of two-dimensional airfoil in the attached flow. Taking into account Theodorsen theory, the unsteady lift on the blade section is given as where is the apparent wind speed, is the chord length, is the geometric angle of attack between chord and free-stream flow, and is the downwash at the three-quarter point, related to the motions of the blade section. is an indicial function and is the dimensionless time-scale. The unsteady lift is a combination of noncirculatory part due to the acceleration of the air mass and the circulatory part with Duhamel’s integral for the memory effect of previously shed vorticity.
can be approximated as where , , , and are constant, indicating two time lags.
is described as
Then, the effective downwash at the three-quarter point can be expressed as
Substituting (10) and (11) into (12) giveswhere the state variables can be written as where . And is the solution to the differential equations aswhere the angle of attack at three-quarter . Although the equation is linear in aerodynamic state variables , the structural states are coupled nonlinearly with through and . This equation describes the effective downwash determining the dynamics of unsteady lift coefficient for attached flow.
With the aerodynamic state variables , , the effect angle of attack can be written as
The unsteady lift for attached flow can be rewritten as where is the angle of attack at zero lift.
Let ; the dynamic lift coefficient in (7) can be written as where the first term represents the dynamics of the attached flow, the second term describes the dynamics of the separated flow, and the last is the added mass term related to the rate of torsional deflection.
The equivalent pressure center can be obtained by the static lift and moment coefficients as
The dynamic moment coefficient in (7) can be written as where is the static moment coefficient, and the other term is the added mass term related to the rate of torsional deflection.
3.2. Linearization of the Nonlinear Aerodynamic Model
For the purpose of stability analysis and control design, the nonlinear aerodynamic model in (15), (18), and (20) needs to be linearized. Assume the small perturbations about the equilibrium as where is the small perturbation of the dynamic state variables, is far more less than one, is the steady state of the variables, and is the steady angle of attack.
Substituting (27), (28) into (15), (20), and (22) gives the steady states of the equations aswith the -order equations aswhere the time constant , , with . Consequently, the nonlinear unsteady lift/moment coefficients in (23), (26) are linearized aswhere , are static lift and moment coefficients at the steady state angle of attack. , are Taylor expansion coefficients with respect to effective angle of attack and separate point separately in (31) and are given as
Let , , and define the state variables of the aerodynamic model, . The state-space equations of the aerodynamic model in (30) are written as
3.3. State-Space Representation of the Equations of the Fluid-Structural System
The state variables of the system are defined as Define , , , , , . Combining the structural model in (7) and the aerodynamic model in (28), (29), the state-space equations of the system are as follows:
4. Aeroelastic Analysis
In order to investigate the aeroelastic stability of the system and the effect of the microtab on vibration suppression of the system, the following model parameters are used as : kg, N/m, , , , , , and . For these parameters, the flapwise vibration mode and the torsional vibration mode have the undamped natural frequencies of 6.28rad/s (1Hz) and 50.24rad/s (8Hz), which are compared to the values of 1.1Hz and 8Hz modal natural frequencies in  for flapwise and torsional deflection separately. The , are computed by numerical simulation with commercial CFD code in FLUENT 6.3 software based on Reynolds-averaged Navier-Stokes equations based on wind tunnel experimental data. Those lift/moment coefficients are taken for nondimensional gurney flap for microtab deployment.
The equations of motions for the aeroelastic model in (31) can be rewritten as the following form:where , , and are derived from structural model, , , , and are derived from aerodynamic model, is the control input matrix, related to the actuation of the microtab, and is constant matrix. It can be seen that the system state matrix has structural components that are coupled with aerodynamic components. The system input matrix and output matrix only have structural components. The aerodynamic force influence on the stability of the blade section is analyzed by looking at the eigenvalues of the aeroelastic system in (32). The complex conjugate eigenvalues with positive real part indicate flutter instability, and the eigenvalues with only positive real part reveal divergence instability. The specific frequency and the damp ratio of the damping mode can be computed by the eigenvalues. Before the aeroelastic stability and control analysis is studied, the established model needs to be validated.
4.1. Model Validation by Comparison of the Results with the Study of Kallesøe 
The vibration modes of classic flutter for the blade section have been investigated by Kallesøe. In his study, the numerical results were given for the blade with flapwise, edgewise, and torsional motions at the low angle of attack for classic flutter study. The edgewise mode was indicated to have very small coupling with the aerodynamic forces, and the combined torsional and flapwise motion is primarily responsible for the instability, known as classic flutter. The flutter mode changes from well damped (flutter instability) to highly unstable (divergence instability) within a relatively small increase in air speed.
The flapwise mode and torsional mode are also carried out in this study in Section 2. When the apparent wind velocity does not reach the flutter speed, the system is stable with the preflutter speed of m/s. On this condition, the structural mode has the system matrix asWith the angle of attack and m/s, the aeroelastic model has the system matrix as
For the structural model, the homogenous part of equation of motion can be written as , with the state vector ; for the aeroelastic model, the homogenous part of equation of motion is , with the state vector . The eigenvalues of two models are shown in Table 1. The structural model has two eigenvalues predominately related to torsional mode (8.17 Hz) and other two mainly associated with flapwise mode (1 Hz), all with small negative real part, corresponding to small structural damping (). The computed nondimensional resonance frequencies () can be converted back to resonance frequencies by , . The aeroelastic model has four eigenvalues more than the structural model, representing four heavily damped aerodynamic states. The torsional mode is a little affected by the effect of aerodynamic force while flapwise mode has more negative real part compared with that in structural model, corresponding to larger flapwise damping. The results are consistent with the model for blade fatigue load control developed by Kallesøe, which showed the same system dynamic for the structural and aeroelastic models in Table 1.
4.2. Aeroelastic Analysis
There are mainly two kinds of unstable system response of the aeroelastic system. One is the flutter instability, which occurs when the apparent wind velocity reaches the flutter speed and the system begins to flutter with oscillations. The other one is the divergence instability when the system is unstable and heavily damped with no oscillations. In this section, both of those two instability cases are investigated for the aeroelastic model.
4.2.1. Flutter Instability
The eigenvectors of the aeroelastic model at flutter speed m/s are shown at Table 2. The eigenvectors associated with the unstable eigenvalues of and reveal that the torsional mode is flutter instability with the imaginary parts and flapwise mode is largely drawn towards the imaginary axis by the inclusion of the aerodynamic forces. These results are reasonable based on the fact that the aerodynamic flow acts mainly in the attached flow region with the low angle of attack. The other two eigenvalues of −104.31 and −15.73 are mainly for the attached flow, which are also coupled with .
4.2.2. Divergence Instability
For m/s, the eigenvectors of the aeroelastic model are given in Table 3. The divergence instability of torsional modes are shown by the eigenvalues of 7.96 and 0.0018 with no imaginary parts, and the corresponding eigenvectors also reveal the coupling of the mode and the attached flow. The eigenvectors of four aerodynamic eigenvalues indicate that, on the condition of divergence instability, the coupling of the aerodynamic force and structural model is even stronger than that on flutter instability condition.
5. Control Analysis
5.1. Control Strategies for Microtab Actuation
In order to implement the microtab actuation effectively for blade vibration control, the sensors need to be chosen carefully to design the controller. Here, four strategies of adopting sensors are investigated for microtab control on flutter instability and divergence instability cases separately. For appropriate output feedback gain , the feedback aeroelastic system can be written as . The different output controllers are with the output 1 of , output 2 of , output 3 of , and output 4 of . In practice, those output signals can be sensed by meter sensor of , angular sensor of , and rate gyro of . For the controllability and observability, since the output feedback controller is used for microtab, the output controllability criteria  are used to evaluate the controllability and observability of the outputs at the same time. The output controllability matrix is computed with the () matrices in (32) and (33). For different output choice, the is taken as , , , and . The aeroelastic model with microtab and four outputs are (), (), (), and () separately, and the corresponding output controllability matrices are all solved with the rank, which is equal to the row number of , , which satisfies the requirement of the criteria. Therefore, the aeroelastic model for microtab is completely output controllable.
5.1.1. Flutter Instability Case (Case 1)
When the flutter instability occurs at m/s, the eigenvalues of the feedback system with the same feedback gain and four different output controllers are shown in Figure 2. The choice of output 2 has no effect on suppressing the vibration modes (, ); output 1 and output 3 can control the torsional mode (, ) at the cost of losing the controllability of flapwise mode (, ); output 4 (derivative of torsional deflection) has the full controllability of both modes (, , ). It reveals that is the most powerful sensor signal for microtab control on flutter instability case. The failure of output 3 shows the stronger impact of than that of in the signal combination, and an adjustment of the weight of in output 3 can retain the controllability.
5.1.2. Divergence Instability Case (Case 2)
On the divergence instability condition of m/s, the eigenvalues of the feedback system with four different sensor strategies and feedback gain are shown in Figure 3. The aeroelastic stability has not been improved by microtab control associated with output 1 (, , ) and output 2 (, , , ). However, both output 3 (, , ) and output 4 (, , ) can stabilize the system, specifically that output 4 has better controllability on the modes than output 3. It indicates that is also the key parameter for microtab control in divergence instability case, and, more importantly, overweighs on the control performance, corresponding to the stable feedback system with output 3.
5.1.3. Compare Flutter Instability and Divergence Instability Cases
Compare the results of flutter instability case with divergence instability case; it reveals that , for both cases, the output of and the output of are incapable of controlling the blade vibration via the microtab, and the output of has the ability of suppressing the vibration modes; the difference of the results of case (flutter instability) and case (divergence instability) lies in the effect of output 3: for case , the feedback system with output 3 can eliminate the system instability while this is not the case for output 3 on case . In order to investigate this phenomenon, the eigenvectors of the system with output 3 for both cases are shown in Table 4. It indicates that, for flutter instability (case ), the coupling of flapwise mode () and torsional mode () is much weaker than that of divergence instability (case ), so that the flapwise deflection is reluctant to be suppressed by the control of , which is the most effective signal in the vibration control; although the feedback systems with output 4 can be stabilized for both cases, the divergence instability system in case () needs more control gain than the flutter instability system in case (), which shows that divergence instability is more difficult to control than the flutter instability. Note that is already good enough to stabilize the system in flutter instability case while it lacks the control ability for the divergence instability case, which needs more power from the controller.