Table of Contents Author Guidelines Submit a Manuscript
Shock and Vibration
Volume 2015 (2015), Article ID 136026, 11 pages
http://dx.doi.org/10.1155/2015/136026
Research Article

Numerical Investigation of Flapwise-Torsional Vibration Model of a Smart Section Blade with Microtab

1School of Hydraulic, Energy and Power Engineering, Yangzhou University, Yangzhou, Jiangsu 225127, China
2Aerospace Engineering, Embry-Riddle Aeronautical University, Daytona Beach, FL 32114-3900, USA
3Wright State Research Institute, AFRL RQVC, Beavercreek, OH 45431, USA

Received 16 April 2015; Revised 19 July 2015; Accepted 27 July 2015

Academic Editor: Peng Chen

Copyright © 2015 Nailu Li 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.

Abstract

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.

1. Introduction

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) [1]. 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 [4]. The trailing-edge flap has been widely used as the control surface to reject the load and control the aeroelastic instability of the blade [57]. For the recent decades, microtabs, proposed by Yen et al. [8], 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. [9] first proved the effectiveness of the microtab for load control numerically. Baek et al. [10] 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. [13] presented microtab design and wind tunnel tests to verify the use of microtabs for active load control. Maughmer and Bramesfeld [14] 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 [15] 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 [16]. 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 [19], 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 [2022]. 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 [2325]. 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.

Figure 1: Model of the blade section with trailing-edge microtabs.

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 [26]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 [17] 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 [25]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

The state-space equations of the dynamic lift and moment coefficients in (25), (26) can be 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 [20]: 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 [20] 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 [20]

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.

Table 1: Eigenvalues for structural and aeroelastic models with  m/s and .
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 .

Table 2: Eigenvectors of aeroelastic models with  m/s.
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.

Table 3: Eigenvectors of aeroelastic models with  m/s.

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 [27] 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.

Figure 2: Eigenvalues λ of the aeroelastic model with different output controllers for flutter instability case ( m/s).
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.

Figure 3: Eigenvalues λ of the aeroelastic model with different output controllers for divergence instability case ( m/s).
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.

Table 4: Eigenvectors of the feedback system with output 3 for case  1 and case  2.
5.2. Control Input of Microtab Actuation

Since the feedback aeroelastic system is chosen with the appropriate output of , the limitation of the microtab deployment needs to be considered for practice (normally 2% of of , is the chord length of the airfoil section). In this section, three different deployments of the microtab are tested to verify the control capability of the microtab on flutter instability and divergence instability cases. The system is taken as , , with the output . , , are corresponding to the different microtab deployment (2% of , 4% of , and 5% of ).

The flutter instability feedback system is with the control gain and ; the divergence instability feedback system is with the control gain and . The eigenvalues of those two systems with different deployment of microtab are shown in Figure 4. For the divergence instability system, the instability has been improved with the increased microtab deployment. Until the full deployment of 5% of , the microtab has the full controllability on vibration modes. For the flutter instability system, 2% of is good enough to suppress the vibration modes. It indicates that the small microtab deployment can have good applicability on the flutter instability control while divergence instability system needs full deployment actuation of the microtab. However, the microtab performs better on the divergence instability system than the flutter instability system with the full deployment of 5% of , corresponding to the increased mode damping.

Figure 4: Eigenvalues of the vibration modes with different microtab deployments and output 4 for (a) flutter instability case and (b) divergence instability case.
5.3. Control Effect of the Microtab

As indicated in Figures 5 and 6, the microtab control has a good performance on the suppression of both the flutter instability and divergence instability of the airfoil section in a short time. Specifically, for flutter instability case, the oscillation of plunge displacement and torsional deflection can be mostly suppressed in 3 seconds and in 4 seconds separately. For divergence instability case, it takes around 4.5 seconds to control the undamped instability of the aeroelastic system, which also reveals that it is a little harder to stabilize the divergence instability than the flutter instability.

Figure 5: Responses of plunge displacement and torsional deflection without control for (a) flutter instability case or (b) divergence instability case.
Figure 6: Responses of plunge displacement and torsional deflection with microtab control for flutter instability and divergence instability cases.

6. Conclusion

A method for developing the aeroelastic model of the smart blade section with the microtab actuator has been presented considering flapwise and torsional modes, aerodynamic forces, and microtab deployment effects. This method allows for the incorporation of the aerodynamic model approximated to the Theodorsen theory and the aerodynamic effects of the microtab into a smart blade vibration model. This model is consistent with the known aeroelastic effects of a high-order blade section mode, and it is appropriate to study controller development via the microtab actuation for passive vibration control of the smart blade section in classic flutter.

Different output feedback controllers have been investigated, and the controller using the system output of derivative of torsional deflection is confirmed as the most effective controller for passive vibration control via microtab. The control ability of the microtab is also analyzed with different actuation deployments on both flutter and divergence instability cases, indicating that flutter instability control needs less microtab deployment than divergence instability control. The simulation results of the microtab control reveal that both the instability cases can be suppressed successfully in a short time of less than 5 seconds. And it needs a little more time to control the divergence instability than the flutter instability.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors acknowledge support of this work through Grant 14KJB480006 from the Natural Science Foundation of the Colleges and Universities in Jiangsu Province, Grant 2014CXJ028 from School Innovation Fund of Yangzhou University, and Grant DESC0001261 from the Department of Energy, USA. The project is also sponsored by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry.

References

  1. M. O. L. Hansen, J. N. Sørensen, S. Voutsinas, N. Sørensen, and H. A. Madsen, “State of the art in wind turbine aerodynamics and aeroelasticity,” Progress in Aerospace Sciences, vol. 42, no. 4, pp. 285–330, 2006. View at Publisher · View at Google Scholar · View at Scopus
  2. D. W. Lobitz, “Parameter sensitivities affecting the flutter speed of a MW-sized blade,” Journal of Solar Energy Engineering, vol. 127, no. 4, pp. 538–543, 2005. View at Publisher · View at Google Scholar · View at Scopus
  3. A. Staino and B. Basu, “Dynamics and control of vibrations in wind turbines with variable rotor speed,” Engineering Structures, vol. 56, pp. 58–67, 2013. View at Publisher · View at Google Scholar · View at Scopus
  4. L. Librescu and P. Marzocca, “Advances in the linear/nonlinear control of aeroelastic structural systems,” Acta Mechanica, vol. 178, no. 3-4, pp. 147–186, 2005. View at Publisher · View at Google Scholar · View at Scopus
  5. D. W. Lobitz, “Flutter speed predictions for MW-Sized wind turbine blades,” Wind Energy, vol. 7, no. 3, pp. 211–224, 2004. View at Google Scholar
  6. D. G. Wilson, D. E. Berg, M. F. Barone, J. C. Berg, B. R. Resor, and D. W. Lobitz, “Active aerodynamic blade control design for load reduction on large wind turbines,” in Proceedings of the European Wind Energy Conference, Marseille, France, 2009.
  7. M. D. Conner, D. M. Tang, E. H. Dowell, and L. N. Virgin, “Nonlinear behavior of a typical airfoil section with control surface freeplay: a numerical and experimental study,” Journal of Fluids and Structures, vol. 11, no. 1, pp. 89–109, 1997. View at Publisher · View at Google Scholar · View at Scopus
  8. D. T. Yen, C. P. Van Dam, F. Bräeuchle, R. L. Smith, and S. D. Collins, “Active load control and lift enhancement using MEM translational tabs,” in Proceedings of the Fluids Conference and Exhibit, AIAA Paper 2000-2422, Denver, Colo, USA, 2000.
  9. C. P. Van Dam, D. Yen Nakafuji, C. Bauer, D. Chao, and K. Standish, “Computational design and analysis of a microtab based aerodynamic loads control system for lifting surfaces,” in MEMS Components and Applications for Industry, Automobiles, Aerospace, and Communication II, vol. 4981 of Proceedings of SPIE, pp. 25–31, San Jose, Calif, USA, 2002.
  10. P. Baek, M. Gaunaa, N. N. Sørensen, and P. Fuglsang, “Comparative study of distributed active load control concepts for wind turbine blades,” in Science of Making Torque from Wind, pp. 611–617, European Academy of Wind Energy, 2010. View at Google Scholar
  11. D. G. Wilson, D. E. Berg, D. W. Lobitz, and J. R. Zayas, “Optimized active aerodynamic blade control for load alleviation on large wind turbines,” in Proceedings of the AWEA Wind Power Conference & Exhibition, pp. 1–4, Houston, Tex, USA, 2008.
  12. J. P. Baker, K. J. Standish, and C. P. van Dam, “Two-dimensional wind tunnel and computational investigation of a microtab modified airfoil,” Journal of Aircraft, vol. 44, no. 2, pp. 563–572, 2007. View at Publisher · View at Google Scholar · View at Scopus
  13. D. T. Y. Nakafuji, C. P. van Dam, R. L. Smith, and S. D. Collins, “Active load control for airfoils using microtabs,” Journal of Solar Energy Engineering, Transactions of the ASME, vol. 123, no. 4, pp. 282–289, 2001. View at Publisher · View at Google Scholar · View at Scopus
  14. M. D. Maughmer and G. Bramesfeld, “Experimental investigation of gurney flaps,” Journal of Aircraft, vol. 45, no. 6, pp. 2062–2067, 2008. View at Publisher · View at Google Scholar · View at Scopus
  15. K. J. Standish, Aerodynamic analysis of blunt trailing edge airfoils & a microtab-based load control system, University of California, Davis, Calif, USA, 2003.
  16. K. J. Standish and C. P. van Dam, “Computational analysis of a microtab-based aerodynamic load control system for rotor blades,” Journal of the American Helicopter Society, vol. 50, no. 3, pp. 249–258, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. T. Theodorsen, “General theory of aerodynamic instability and the mechanism of flutter,” NACA Report 496, 1935. View at Google Scholar
  18. P. P. Friedmann, “Recent trends in rotary-wing aeroelasticity,” Vertica, vol. 11, no. 1-2, pp. 139–170, 1987. View at Google Scholar
  19. B. G. Van der Wall and J. G. Leishman, “On the influence of time-varying flow velocity on unsteady aerodynamics,” Journal of the American Helicopter Society, vol. 39, no. 4, pp. 25–36, 1994. View at Publisher · View at Google Scholar · View at Scopus
  20. B. S. Kallesøe, “A low-order model for analysing effects of blade fatigue load control,” Wind Energy, vol. 9, no. 5, pp. 421–436, 2006. View at Publisher · View at Google Scholar · View at Scopus
  21. T. Liu, Y. Ren, and X. Yang, “Nonlinear aeroelastic stability analysis of wind turbine blade with bending-bending-twist coupling,” Journal of Fluids and Structures, vol. 42, pp. 488–502, 2013. View at Publisher · View at Google Scholar · View at Scopus
  22. J. W. Larsen and S. R. K. Nielsen, “Non-linear dynamics of wind turbine wings,” International Journal of Non-Linear Mechanics, vol. 41, no. 5, pp. 629–643, 2006. View at Publisher · View at Google Scholar · View at Scopus
  23. J. W. Larsen and S. R. K. Nielsen, “Nonlinear parametric instability of wind turbine wings,” Journal of Sound and Vibration, vol. 299, no. 1-2, pp. 64–82, 2007. View at Publisher · View at Google Scholar · View at Scopus
  24. J. W. Larsen, R. Iwankiewicz, and S. R. K. Nielsen, “Nonlinear stochastic stability analysis of wind turbine wings by Monte Carlo simulations,” Probabilistic Engineering Mechanics, vol. 22, no. 2, pp. 181–193, 2007. View at Publisher · View at Google Scholar · View at Scopus
  25. B. S. Kallesøe, “Equations of motion for a rotor blade, including gravity, pitch action and rotor speed variations,” Wind Energy, vol. 10, no. 3, pp. 209–230, 2007. View at Publisher · View at Google Scholar · View at Scopus
  26. G. Platanitis and T. W. Strganac, “Control of a nonlinear wing section using leading- and trailing-edge surfaces,” Journal of Guidance, Control, and Dynamics, vol. 27, no. 1, pp. 52–58, 2004. View at Publisher · View at Google Scholar · View at Scopus
  27. K. Ogata, Modern Control Engineering, Prentice Hall, New York, NY, USA, 5th edition, 2009.