# Stall Flutter Control of a Smart Blade Section Undergoing Asymmetric Limit Oscillations

Authors: Nailu Li; Mark J. Balas; Pourya Nikoueeyan; Hua Yang; Jonathan W. Naughton
Journal: Shock and Vibration
 (2016)

---

## Abstract

Stall flutter is an aeroelastic phenomenon resulting in unwanted oscillatory loads on the blade, such as wind turbine blade, helicopter rotor blade, and other flexible wing blades. Although the stall flutter and related aeroelastic control have been studied theoretically and experimentally, microtab control of asymmetric limit cycle oscillations (LCOs) in stall flutter cases has not been generally investigated. This paper presents an aeroservoelastic model to study the microtab control of the blade section undergoing moderate stall flutter and deep stall flutter separately. The effects of different dynamic stall conditions and the consequent asymmetric LCOs for both stall cases are simulated and analyzed. Then, for the design of the stall flutter controller, the potential sensor signal for the stall flutter, the microtab control capability of the stall flutter, and the control algorithm for the stall flutter are studied. The improvement and the superiority of the proposed adaptive stall flutter controller are shown by comparison with a simple stall flutter controller.


---

## Body

## 1. Introduction

Stall flutter is a nonlinear aeroelastic phenomenon with self-excited oscillation of limited amplitude. Those oscillations are called limit cycle oscillations (LCOs). And stall flutter is a LCO occurring as a result of the dynamic stall, which is an aerodynamic nonlinearity [1]. Two dynamic stall regimes are defined [2]: moderate stall and deep stall. Moderate stall is a phenomenon similar to static stall. Deep stall is characterized by the maximum lift and moment values far in excess of their static counterparts. The oscillations of stall flutter associated with moderate stall or deep stall can result in large oscillatory fatigue loads on the wind turbine blade or helicopter rotor blade and thus reduce blade life span as well as harming the safety and the efficiency of the blade operation. Therefore, there had been an interest in the stall flutter suppression for years.For the past decades, flutter control has been studied by a number of researchers. Yan and Xin [3] studied the oscillatory blowing control of the airfoil flutter by numerical simulation. Nonlinear control methods were reviewed by Kurdila and Akella [4] for high-energy limit cycle oscillations. They concluded that full and partial feedback control or nonlinear model reference adaptive control (MRAC) methodologies were effective to control LCOs. Strganac et al. [5] suggested the adaptive control as the appropriate control strategy for limit cycle oscillations in aeroelastic systems. Those studies mainly emphasized the design of the control algorithms.Recently, there has been an increasing interest in flutter control via the smart blade, which is the blade equipped with the appropriate and deployed actuator [6]. Various kinds of such actuators have been investigated for flutter suppression. Tang and Dowell [7] completed the LCO suppression of a high-aspect ratio wing by a controllable slender body tip mass distribution. Li and Fleeter [8] adopted piezoelectric actuators for active suppression of nonlinear stall flutter, including limit cycle and chaotic and quasi-periodic separated flow induced airfoil vibration. Sun et al. [9] studied stall flutter and flutter suppression control systems with a novel state-space model description. Liu [10] investigated blade pitch actuation in the aeroservoelastic control of stall-induced flutter of wind turbine blade section. In addition, the trailing edge flap has been widely used as the control surface to reject the unwanted oscillatory loads and control the aeroelastic instability of the blade [11–13].However, the flutter control using the microtab is rarely studied. The microtab, earlier proposed by Yen et al. [14], is simply short flat plates (approximately 1–5% of the airfoil chord in height) attached to the trailing edge perpendicular to the chord line on the pressure or suction side of the airfoil. In contrast to many other actuators for addressing blade loading issues that seem to be complex and energy-consuming, the potential for the microtab to regulate blade loads with high efficiency has been numerically and experimentally studied by researchers [15–20]. The author has also made several studies on microtab control. First, the effect of microtab was studied to control small sinusoidal aerodynamic loads on the blade section [21]. Then, the microtab was used to control aerodynamic load of the blade section at low wind speeds associated with operating regions of a typical wind turbine [22, 23]. The microtab controller was also designed to have full controllability on pitch motion [24]. Moreover, the microtab was recently investigated for the control of classic flutter at relatively small angle of attack [25].However, the microtab control of the blade section undergoing dynamic stall has not been studied extensively. In most of the previous works with microtab control application, the static aerodynamic force coefficients have been considered as a baseline for the airfoil loading behavior or a small and smooth change in the unsteady aerodynamic loading to control the aeroelastic response of the blade section only in the attached flow. Nevertheless, it has been proved that, with the occurrence of the dynamic stall, the airfoil section suffers a dramatic change in aerodynamic loads in short time intervals, which will worsen the control effect and even cause the instability of the control system. Thus, it is essential to develop new microtab controller to address stall flutter in different dynamic stall regimes.The purpose of the present work is to investigate a new microtab controller to suppress the asymmetric LCOs in moderate stall flutter and deep stall flutter conditions. The dynamic moderate/deep stall of the blade section will be simulated via the Beddoes-Leishman dynamic stall model and validated by the static test curves. The consequent limit cycle oscillations of the moderate/deep stall flutter will be studied by analyzing the time-domain responses and phase plane results. The design of the stall flutter controller with microtab actuator includes the sensor signal study, the microtab control capability analysis, and the stall flutter control algorithm design. The control results of the proposed controller are compared with the simple stall flutter controller to show its improvement and effectiveness in both dynamic stall conditions.


## 2. Aeroelastic Model of Stall Flutter System

In this study, an aeroelastic model of the blade section shown in Figure1 is used for aeroelastic analysis and control design. This is a typical two-dimensional wing section with two degrees of freedom, that is, the plunge motion and the pitch motion. The plunge motion is defined by the plunge displacement h with the direction of positive downwards, and the pitch motion is defined by the pitch angle θ with the direction of positive nose down. The aerodynamic control device, microtabs, is installed close to the trailing edge of the blade. The microtab deployment height g, is regarded as a control input to the system and acts in a similar way to the gurney flap. b is the mid-chord length. K
h and K
θ are plunge stiffness and pitch stiffness, restricted to the elastic axis and C
h and C
θ are the plunge and pitch damping coefficients. The section is exposed to the aerodynamic forces L and aerodynamic moment M at the aerodynamic centre.Figure 1
Typical two-dimensional blade section.


### 2.1. Dynamic Stall Model

To include the loading from the unsteady dynamic stall phenomena, the Beddoes-Leishman (B-L) model is adopted to account for the effect of the separated flow at the stall region. Since the basic assumption in B-L model of static lift for trailing edge separation is the lift on a flat plate in a potential Kirchhoff flow, the separation pointf
s
t can be defined as follows [26]:(1)
f
s
t
=
2
C
L
s
t
α
C
L
,
α
α
-
α
0
-
1
2
,where C
L
s
t is the static lift coefficient for the airfoil and f
s
t indicates a fully separated flow when it equals zero.Based on the separation point, the fully separated flow is described as(2)
C
L
f
s
=
C
L
s
t
-
C
L
,
α
α
-
α
0
f
s
t
1
-
f
s
t
.Note that there will be a problem whenf
s
t
=
1. The value of C
L
f
s can be approximated as C
L
s
t
(
α
)
/
2 when f
s
t gets close to one.The dynamic behavior of the trailing edge separation can be described by the separation related to the pressure distribution over the airfoil section and the separation due to the dynamics of the boundary layer. For the pressure distribution separation, it is assumed that there is a time lag between the unsteady lift coefficientC
L
P
′and the actual lift coefficient C
L
P. The unsteady lift coefficient is modeled as(3)
C
L
P
′
˙
+
T
p
-
1
C
L
P
′
=
T
p
-
1
C
L
p
t
,where T
p is the time constant for the pressure lag and C
L
P
=
2
π
α
E
-
α
0
+
π
b
θ
/
U
˙
,
α
0 is the angle of attack at zero lift.The overall dynamics of the flow separation can be revealed by the unsteady separation point, which is caused by the boundary layer dynamics, and involves the effect of the pressure distribution. The equation for the unsteady separation pointf
d
y, which makes the separation point lag behind the quasi-steady value f
q
s, is given as(4)
f
˙
d
y
+
T
f
-
1
f
d
y
=
T
f
-
1
f
q
s
t
,where T
f is the time delay in the boundary layer; the equivalent quasi-steady separation point f
q
s
=
f
s
t
(
α
f
), α
f
=
C
L
P
′
/
C
L
,
α
+
α
0, is an equivalent angle of attack that gives the same quasi-steady value of unsteady lift coefficient C
L
P
′.The unsteady aerodynamic forces of the blade section in stall should involve the unsteady forces induced by both attached flow and the separation flow. The unsteady lift curve can be represented by a liner interpolation between the lift for the fully attached flow and the lift for the fully separated flow as(5)
L
s
=
ρ
U
2
b
C
L
0
+
L
1
f
d
y
+
ρ
U
2
b
C
L
f
s
α
E
1
-
f
d
y
,where L
s is the unsteady lift including the effect of trailing edge separation, C
L
0is the static lift coefficients at zero angle of attack, L
1 is the lift for the fully attached flow, and α
E is the effective angle of attack and can be given as (6)
α
E
t
=
b
h
˙
t
U
+
θ
t
+
b
θ
˙
t
U
1
-
A
1
-
A
2
+
x
1
+
x
2
,
(7)
x
1
=
b
1
A
1
1
b
U
∫
0
t
ω
3
/
4
t
′
U
t
′
e
-
b
1
/
b
∫
t
′
t
U
τ
d
τ
d
t
′
,
x
2
=
b
2
A
2
1
b
U
∫
0
t
ω
3
/
4
t
′
U
t
′
e
-
b
2
/
b
∫
t
′
t
U
τ
d
τ
d
t
′
,where x
1 and x
2, yielding the effective downwash, are two state variables of the attached flow; ω
3
/
4 is the downwash at the three-quarter point. A
1, A
2, b
1, and b
2 are constant, indicating two time lags in Wagner function of ϕ
s
=
1
-
A
1
e
-
b
1
s
-
A
2
e
-
b
2
s.The unsteady moment is affected by the trailing edge separation through movement of the pressure center due to the separation. Define the equivalent pressure center as(8)
a
s
t
=
C
M
s
t
-
C
M
0
C
L
s
t
,where C
M
s
t and C
L
s
t are the static lift and moment coefficient and C
M
0 is the moment coefficient at zero lift. Thus, the unsteady moment is given as (9)
M
s
=
2
ρ
U
2
b
2
C
M
0
+
M
1
+
2
ρ
U
2
b
2
L
s
a
s
t
f
d
y
-
a
s
t
f
s
t
α
E
,where the first term is the static moment at zero angle of attack, the second term M
1 is the moment for the fully attached flow, and the last term is the unsteady moment due to the dynamic trailing edge separation.Assume that the dynamic stall model is linearized with small amplitude vibrations about the equilibrium, given by the steady angle of attackα
0. Then the linearized unsteady lift and moment using Taylor expansion about the steady angle of attack are approximated as (10)
L
s
=
ρ
U
2
b
C
L
0
+
L
1
f
0
+
ρ
U
2
b
α
E
d
C
L
f
s
d
α
α
=
α
0
1
-
f
0
+
ρ
U
2
b
C
L
,
α
α
0
-
α
0
-
C
L
f
s
α
0
f
d
y
,
M
s
=
2
ρ
U
2
b
2
C
M
0
+
M
1
-
2
ρ
U
2
b
2
α
E
C
L
0
d
f
s
t
d
α
α
=
α
0
d
a
s
t
d
f
α
=
α
0
+
2
ρ
U
2
b
2
C
L
0
d
a
s
t
d
f
α
=
α
0
f
d
y
,where C
L
0 and C
M
0 are the static lift and moment coefficients at the steady angle of attack. f
0 is the static separation point at the steady angle of attack. The details of the formulation of the attached lift L
1 and moment M
1 can be found in the literature [25].


### 2.2. Governing Equations of the System

With the dynamic stall modeled in the previous section, the equations of the pitch and plunge motions of the smart blade section undergoing the stall dynamics are as follows:(11)
m
m
b
x
θ
m
b
x
θ
I
θ
h
¨
θ
¨
+
C
h
0
0
C
θ
h
˙
θ
˙
+
K
h
0
0
K
θ
h
θ
=
-
L
s
M
s
+
-
L
g
M
g
,where m is the mass of the airfoil section, x
θ is the dimensionless distance between the elastic axis and the mass center, and I
θ is the inertia moment about the elastic axis.Definez
=
h
θ
T, F
=
-
L
s
M
s
T, and F
g
=
-
L
g
M
g
T; (11) of the structural model can be rewritten in the form as (12)
M
0
z
¨
+
C
0
z
˙
+
K
0
z
=
F
+
F
g
,where F is the unsteady aerodynamic forces, F
g is the aerodynamic forces due to control surface, L
g
=
ρ
U
2
b
C
L
,
g
g and M
g
=
2
ρ
U
2
b
2
C
M
,
g
g represent the change in the airfoil lift and pitching moment due to the microtab deployment, and C
L
,
g and C
M
,
g are the corresponding aerodynamic coefficients. These coefficients were calculated by numerical simulation using the commercial CFD code, FLUENT.DefiningX
a as the vector of x
1 and x
2 and X
s as the vector of C
L
P
′ and f
d
y, (10) can be rewritten in the matrix form as (13)
F
=
-
L
s
M
s
=
D
1
z
+
D
2
z
˙
+
C
1
X
a
+
C
2
X
s
+
C
0
,where (14)
D
1
=
0
-
ρ
U
2
b
a
-
C
L
,
α
f
0
+
a
-
C
L
,
α
0
f
s
1
-
f
0
0
2
ρ
U
2
b
2
a
-
C
M
,
α
0
s
t
-
a
-
C
L
0
f
α
0
s
t
a
α
0
s
t
,
D
2
=
-
ρ
U
a
-
b
2
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
-
ρ
U
a
-
b
2
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
+
π
a
-
2
ρ
U
a
-
b
3
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
2
ρ
U
a
-
b
3
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
-
π
2
a
-
,
C
1
=
-
ρ
U
2
b
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
-
ρ
U
2
b
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
2
ρ
U
2
b
2
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
2
ρ
U
2
b
2
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
,
C
2
=
0
-
ρ
U
2
b
C
L
,
α
α
0
-
α
0
-
C
L
f
s
α
0
0
2
ρ
U
2
b
2
C
L
0
a
α
0
s
t
,
C
0
=
-
ρ
U
2
b
C
L
0
2
ρ
U
2
b
2
C
M
0
,
a
-
=
1
-
A
1
-
A
2
,
C
L
,
α
0
f
s
=
d
C
L
f
s
d
α
α
=
α
0
,
C
M
,
α
0
s
t
=
d
C
M
s
t
d
α
α
=
α
0
,
f
α
0
s
t
=
d
f
s
t
d
α
α
=
α
0
,
a
α
0
s
t
=
d
a
s
t
d
f
α
=
α
0
.With the definition ofX
a and X
s, (3), (4), (6), and (7) of the aerodynamic stall model can be rewritten in the state-space forms as (15)
X
˙
s
=
A
3
X
a
+
A
4
X
s
+
B
3
z
+
B
4
z
˙
,
X
˙
a
=
A
1
X
a
+
A
2
X
s
+
B
1
z
+
B
2
z
˙
,where(16)
A
1
=
-
T
1
-
1
0
0
-
T
2
-
1
,
A
2
=
0
0
0
0
A
3
=
T
p
-
1
C
L
,
α
T
p
-
1
C
L
,
α
0
0
,
A
4
=
-
T
p
-
1
0
T
f
-
1
C
L
,
α
f
α
0
s
t
-
T
f
-
1
,
B
1
=
0
T
1
-
1
A
1
0
T
2
-
1
A
2
,
B
2
=
T
1
-
1
A
1
b
U
T
1
-
1
A
1
b
U
T
2
-
1
A
2
b
U
T
2
-
1
A
2
b
U
,
B
3
=
0
a
-
T
p
-
1
C
L
,
α
0
0
,
B
4
=
a
-
b
T
p
-
1
C
L
,
α
U
a
-
b
T
p
-
1
C
L
,
α
U
+
π
T
P
-
1
b
U
0
0
.The aeroservoelastic model is a combination of the structural model in (12) and the aerodynamic model in (15). Defining the state variables of the system as q
=
z
z
˙
X
a
X
s
T, the governing equations of the aeroservoelastic model can be written in state-space form as (17)
z
˙
z
¨
X
˙
a
X
˙
s
=
0
2
×
2
I
2
×
2
0
2
×
2
0
2
×
2
-
M
0
-
1
K
0
-
D
1
-
M
0
-
1
C
0
-
D
2
M
0
-
1
C
1
M
0
-
1
C
2
B
1
B
2
A
1
A
2
B
3
B
4
A
3
A
4
︸
A
z
z
˙
X
a
X
s
+
0
2
×
1
G
0
2
×
1
0
2
×
1
︸
B
g
,where(18)
G
=
-
ρ
U
2
b
I
θ
C
L
,
g
-
2
ρ
U
2
b
3
m
x
θ
C
M
,
g
m
I
θ
-
m
2
b
2
x
θ
2
ρ
U
2
b
2
m
x
θ
C
L
,
g
+
2
ρ
U
2
b
2
m
C
M
,
g
m
I
θ
-
m
2
b
2
x
θ
2
.Note that (17) is a linearized equation, with the state matrix A, containing the corresponding matrices from the aerodynamic stall model and the structural model.


## 2.1. Dynamic Stall Model

To include the loading from the unsteady dynamic stall phenomena, the Beddoes-Leishman (B-L) model is adopted to account for the effect of the separated flow at the stall region. Since the basic assumption in B-L model of static lift for trailing edge separation is the lift on a flat plate in a potential Kirchhoff flow, the separation pointf
s
t can be defined as follows [26]:(1)
f
s
t
=
2
C
L
s
t
α
C
L
,
α
α
-
α
0
-
1
2
,where C
L
s
t is the static lift coefficient for the airfoil and f
s
t indicates a fully separated flow when it equals zero.Based on the separation point, the fully separated flow is described as(2)
C
L
f
s
=
C
L
s
t
-
C
L
,
α
α
-
α
0
f
s
t
1
-
f
s
t
.Note that there will be a problem whenf
s
t
=
1. The value of C
L
f
s can be approximated as C
L
s
t
(
α
)
/
2 when f
s
t gets close to one.The dynamic behavior of the trailing edge separation can be described by the separation related to the pressure distribution over the airfoil section and the separation due to the dynamics of the boundary layer. For the pressure distribution separation, it is assumed that there is a time lag between the unsteady lift coefficientC
L
P
′and the actual lift coefficient C
L
P. The unsteady lift coefficient is modeled as(3)
C
L
P
′
˙
+
T
p
-
1
C
L
P
′
=
T
p
-
1
C
L
p
t
,where T
p is the time constant for the pressure lag and C
L
P
=
2
π
α
E
-
α
0
+
π
b
θ
/
U
˙
,
α
0 is the angle of attack at zero lift.The overall dynamics of the flow separation can be revealed by the unsteady separation point, which is caused by the boundary layer dynamics, and involves the effect of the pressure distribution. The equation for the unsteady separation pointf
d
y, which makes the separation point lag behind the quasi-steady value f
q
s, is given as(4)
f
˙
d
y
+
T
f
-
1
f
d
y
=
T
f
-
1
f
q
s
t
,where T
f is the time delay in the boundary layer; the equivalent quasi-steady separation point f
q
s
=
f
s
t
(
α
f
), α
f
=
C
L
P
′
/
C
L
,
α
+
α
0, is an equivalent angle of attack that gives the same quasi-steady value of unsteady lift coefficient C
L
P
′.The unsteady aerodynamic forces of the blade section in stall should involve the unsteady forces induced by both attached flow and the separation flow. The unsteady lift curve can be represented by a liner interpolation between the lift for the fully attached flow and the lift for the fully separated flow as(5)
L
s
=
ρ
U
2
b
C
L
0
+
L
1
f
d
y
+
ρ
U
2
b
C
L
f
s
α
E
1
-
f
d
y
,where L
s is the unsteady lift including the effect of trailing edge separation, C
L
0is the static lift coefficients at zero angle of attack, L
1 is the lift for the fully attached flow, and α
E is the effective angle of attack and can be given as (6)
α
E
t
=
b
h
˙
t
U
+
θ
t
+
b
θ
˙
t
U
1
-
A
1
-
A
2
+
x
1
+
x
2
,
(7)
x
1
=
b
1
A
1
1
b
U
∫
0
t
ω
3
/
4
t
′
U
t
′
e
-
b
1
/
b
∫
t
′
t
U
τ
d
τ
d
t
′
,
x
2
=
b
2
A
2
1
b
U
∫
0
t
ω
3
/
4
t
′
U
t
′
e
-
b
2
/
b
∫
t
′
t
U
τ
d
τ
d
t
′
,where x
1 and x
2, yielding the effective downwash, are two state variables of the attached flow; ω
3
/
4 is the downwash at the three-quarter point. A
1, A
2, b
1, and b
2 are constant, indicating two time lags in Wagner function of ϕ
s
=
1
-
A
1
e
-
b
1
s
-
A
2
e
-
b
2
s.The unsteady moment is affected by the trailing edge separation through movement of the pressure center due to the separation. Define the equivalent pressure center as(8)
a
s
t
=
C
M
s
t
-
C
M
0
C
L
s
t
,where C
M
s
t and C
L
s
t are the static lift and moment coefficient and C
M
0 is the moment coefficient at zero lift. Thus, the unsteady moment is given as (9)
M
s
=
2
ρ
U
2
b
2
C
M
0
+
M
1
+
2
ρ
U
2
b
2
L
s
a
s
t
f
d
y
-
a
s
t
f
s
t
α
E
,where the first term is the static moment at zero angle of attack, the second term M
1 is the moment for the fully attached flow, and the last term is the unsteady moment due to the dynamic trailing edge separation.Assume that the dynamic stall model is linearized with small amplitude vibrations about the equilibrium, given by the steady angle of attackα
0. Then the linearized unsteady lift and moment using Taylor expansion about the steady angle of attack are approximated as (10)
L
s
=
ρ
U
2
b
C
L
0
+
L
1
f
0
+
ρ
U
2
b
α
E
d
C
L
f
s
d
α
α
=
α
0
1
-
f
0
+
ρ
U
2
b
C
L
,
α
α
0
-
α
0
-
C
L
f
s
α
0
f
d
y
,
M
s
=
2
ρ
U
2
b
2
C
M
0
+
M
1
-
2
ρ
U
2
b
2
α
E
C
L
0
d
f
s
t
d
α
α
=
α
0
d
a
s
t
d
f
α
=
α
0
+
2
ρ
U
2
b
2
C
L
0
d
a
s
t
d
f
α
=
α
0
f
d
y
,where C
L
0 and C
M
0 are the static lift and moment coefficients at the steady angle of attack. f
0 is the static separation point at the steady angle of attack. The details of the formulation of the attached lift L
1 and moment M
1 can be found in the literature [25].


## 2.2. Governing Equations of the System

With the dynamic stall modeled in the previous section, the equations of the pitch and plunge motions of the smart blade section undergoing the stall dynamics are as follows:(11)
m
m
b
x
θ
m
b
x
θ
I
θ
h
¨
θ
¨
+
C
h
0
0
C
θ
h
˙
θ
˙
+
K
h
0
0
K
θ
h
θ
=
-
L
s
M
s
+
-
L
g
M
g
,where m is the mass of the airfoil section, x
θ is the dimensionless distance between the elastic axis and the mass center, and I
θ is the inertia moment about the elastic axis.Definez
=
h
θ
T, F
=
-
L
s
M
s
T, and F
g
=
-
L
g
M
g
T; (11) of the structural model can be rewritten in the form as (12)
M
0
z
¨
+
C
0
z
˙
+
K
0
z
=
F
+
F
g
,where F is the unsteady aerodynamic forces, F
g is the aerodynamic forces due to control surface, L
g
=
ρ
U
2
b
C
L
,
g
g and M
g
=
2
ρ
U
2
b
2
C
M
,
g
g represent the change in the airfoil lift and pitching moment due to the microtab deployment, and C
L
,
g and C
M
,
g are the corresponding aerodynamic coefficients. These coefficients were calculated by numerical simulation using the commercial CFD code, FLUENT.DefiningX
a as the vector of x
1 and x
2 and X
s as the vector of C
L
P
′ and f
d
y, (10) can be rewritten in the matrix form as (13)
F
=
-
L
s
M
s
=
D
1
z
+
D
2
z
˙
+
C
1
X
a
+
C
2
X
s
+
C
0
,where (14)
D
1
=
0
-
ρ
U
2
b
a
-
C
L
,
α
f
0
+
a
-
C
L
,
α
0
f
s
1
-
f
0
0
2
ρ
U
2
b
2
a
-
C
M
,
α
0
s
t
-
a
-
C
L
0
f
α
0
s
t
a
α
0
s
t
,
D
2
=
-
ρ
U
a
-
b
2
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
-
ρ
U
a
-
b
2
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
+
π
a
-
2
ρ
U
a
-
b
3
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
2
ρ
U
a
-
b
3
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
-
π
2
a
-
,
C
1
=
-
ρ
U
2
b
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
-
ρ
U
2
b
C
L
,
α
f
0
+
C
L
,
α
0
f
s
1
-
f
0
2
ρ
U
2
b
2
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
2
ρ
U
2
b
2
C
M
,
α
0
s
t
-
C
L
0
f
α
0
s
t
a
α
0
s
t
,
C
2
=
0
-
ρ
U
2
b
C
L
,
α
α
0
-
α
0
-
C
L
f
s
α
0
0
2
ρ
U
2
b
2
C
L
0
a
α
0
s
t
,
C
0
=
-
ρ
U
2
b
C
L
0
2
ρ
U
2
b
2
C
M
0
,
a
-
=
1
-
A
1
-
A
2
,
C
L
,
α
0
f
s
=
d
C
L
f
s
d
α
α
=
α
0
,
C
M
,
α
0
s
t
=
d
C
M
s
t
d
α
α
=
α
0
,
f
α
0
s
t
=
d
f
s
t
d
α
α
=
α
0
,
a
α
0
s
t
=
d
a
s
t
d
f
α
=
α
0
.With the definition ofX
a and X
s, (3), (4), (6), and (7) of the aerodynamic stall model can be rewritten in the state-space forms as (15)
X
˙
s
=
A
3
X
a
+
A
4
X
s
+
B
3
z
+
B
4
z
˙
,
X
˙
a
=
A
1
X
a
+
A
2
X
s
+
B
1
z
+
B
2
z
˙
,where(16)
A
1
=
-
T
1
-
1
0
0
-
T
2
-
1
,
A
2
=
0
0
0
0
A
3
=
T
p
-
1
C
L
,
α
T
p
-
1
C
L
,
α
0
0
,
A
4
=
-
T
p
-
1
0
T
f
-
1
C
L
,
α
f
α
0
s
t
-
T
f
-
1
,
B
1
=
0
T
1
-
1
A
1
0
T
2
-
1
A
2
,
B
2
=
T
1
-
1
A
1
b
U
T
1
-
1
A
1
b
U
T
2
-
1
A
2
b
U
T
2
-
1
A
2
b
U
,
B
3
=
0
a
-
T
p
-
1
C
L
,
α
0
0
,
B
4
=
a
-
b
T
p
-
1
C
L
,
α
U
a
-
b
T
p
-
1
C
L
,
α
U
+
π
T
P
-
1
b
U
0
0
.The aeroservoelastic model is a combination of the structural model in (12) and the aerodynamic model in (15). Defining the state variables of the system as q
=
z
z
˙
X
a
X
s
T, the governing equations of the aeroservoelastic model can be written in state-space form as (17)
z
˙
z
¨
X
˙
a
X
˙
s
=
0
2
×
2
I
2
×
2
0
2
×
2
0
2
×
2
-
M
0
-
1
K
0
-
D
1
-
M
0
-
1
C
0
-
D
2
M
0
-
1
C
1
M
0
-
1
C
2
B
1
B
2
A
1
A
2
B
3
B
4
A
3
A
4
︸
A
z
z
˙
X
a
X
s
+
0
2
×
1
G
0
2
×
1
0
2
×
1
︸
B
g
,where(18)
G
=
-
ρ
U
2
b
I
θ
C
L
,
g
-
2
ρ
U
2
b
3
m
x
θ
C
M
,
g
m
I
θ
-
m
2
b
2
x
θ
2
ρ
U
2
b
2
m
x
θ
C
L
,
g
+
2
ρ
U
2
b
2
m
C
M
,
g
m
I
θ
-
m
2
b
2
x
θ
2
.Note that (17) is a linearized equation, with the state matrix A, containing the corresponding matrices from the aerodynamic stall model and the structural model.


## 3. Stall Conditions and Stall Flutter Analysis

The purpose of this section is to investigate the effect of different dynamic stall conditions and the consequent aeroelastic behavior of the section undergoing stall flutter. Naughton et al. [27] found that the mean angle of attack is one of the important parameters that affect the behavior of the airfoil in stall significantly. Clearly, the extent of the dynamic stall phenomenon can also affect the aeroelastic behavior of the airfoil. Thus, the study focuses on two different cases of dynamic stall regime. The first case is the moderate stall condition, which holds the mean angle of attack right before the static stall angle (α
s
=
12
°). The second case is the deep stall condition, which has the mean angle of attack beyond the static stall angle. In deep dynamic stall condition, large levels of hysteresis exist in the lift and moment coefficient curves, as opposed to the moderate dynamic stall case where small hysteresis in aerodynamic coefficients is seen. The aeroelastic behavior of the airfoil induced by both cases will be discussed. In order to investigate aeroelastic stability, the following model parameters are used [28]: m
=
50 kg, K
h
=
2884 N/m, C
h
=
7.54 kg/s, K
θ
=
700, C
θ
=
0.23 kgm2/s, I
θ
=
0.4 m2kg, ρ
=
1.2 kg/m3, and b
=
0.5 m.


### 3.1.  Moderate Dynamic Stall

Before investigating the aeroelastic response of the airfoil under moderate stall condition, the moderate stall behavior of the blade section needs to be well predicted by the simulation model. The moderate stall aerodynamic curves are simulated based on static experimental data.Figure2 shows hysteresis in the dynamic lift and moment curves. The dynamic lift loop has small overshoot based on the static curves, the dynamic moment loop is well surrounding the static curves, and the dynamic curves follow the trend of the static curves as indicated in [26]. For the dynamic lift curve, it is clear that the moderate stall is indicated by a gentle drop between max. point and reattached point and a small time lag between static stall point and max. point. For the dynamic moment curve, there is also a relatively smooth drop between max. point and reattached point. The max. point before the static point is caused by the dynamic moment shape that follows the trend of static curve before the static stall point.Figure 2
Moderate dynamic stall simulated curve and static experimental data. (a) Lift force coefficient and (b) moment coefficient.
(a)
(b)The pitch and plunge responses at moderate stall flutter are shown in Figure3, where the responses tend to damp at first 30 seconds, and they are stabilized after 40 seconds with asymmetric LCOs (e.g., plunge displacement LCO with magnitude from 0 m to −0.5 m and pitch angle LOC with magnitude from −5 deg. to −10 deg.). This phenomenon is convinced by the experimental results in [1], which indicated the occurrence of only asymmetric LCOs with either positive or negative pitching centering.Figure 3
Pitch and plunge responses at moderate stall flutter.


### 3.2. Deep Dynamic Stall

The deep stall aerodynamic curves are also simulated based on static experimental data. Figure4 shows the deep stall motions of the simulated unsteady lift and moment coefficients at angle of attack α
=
13.5
°
+
3.5
°
s
i
n
⁡
(
10
t
). The deep stall curves are also well predicted by the proposed model to capture the trend of static curves in the neighborhood as expected. In contrast with the moderate stall in Figure 2, the characteristics of deep stall dynamic are revealed by a relatively sharp drop between max. point and reattached point due to the strong separation of the flow and more drastic hysteresis with significant difference between up-stroke and down-stroke phase of dynamic curves. In addition, the maximum lift and moment point are far beyond the static point with a large time lag.Figure 4
Deep dynamic stall simulated curve and static experimental data. (a) Lift force coefficient and (b) moment coefficient.
(a)
(b)The stall phenomena in Figures2 and 4 can also be confirmed by the model validation results in [26], where the lift loop near max. point has a very different shape owing to the separation effects, and the lift curve has expected overshoot, while moment curves revolve around the static curves.Figure5 shows the pitch and plunge responses at deep stall flutter, where the asymmetric LCOs are also observed and their magnitudes are invariant after 25 seconds, resulting in magnitude of the plunge from 0 m to −1 m and magnitude of the pitch from −5 deg. to −15 deg. Compared with the moderate stall flutter, the deep stall can further decrease the damping ratio of the pitch and plunge modes, leading to relatively larger LCOs in a shorter time.Figure 5
Pitch and plunge responses at deep stall flutter.Theθ
-
θ
˙ phase plane is illustrated in Figure 6 for both moderate stall flutter case and deep stall flutter case. Two asymmetric limit cycle oscillations are observed. It is revealed that the asymmetric LCOs grow in size, and pitching center moves towards the left half of the plane as the section goes through more significant stall hysteresis, which can bring severe oscillations. This phase analysis result is consistent with the results of time-domain responses in Figures 3 and 5.Figure 6
Phase plane ofθ and θ
˙ of the system at moderate/deep stall flutter.


## 3.1.  Moderate Dynamic Stall

Before investigating the aeroelastic response of the airfoil under moderate stall condition, the moderate stall behavior of the blade section needs to be well predicted by the simulation model. The moderate stall aerodynamic curves are simulated based on static experimental data.Figure2 shows hysteresis in the dynamic lift and moment curves. The dynamic lift loop has small overshoot based on the static curves, the dynamic moment loop is well surrounding the static curves, and the dynamic curves follow the trend of the static curves as indicated in [26]. For the dynamic lift curve, it is clear that the moderate stall is indicated by a gentle drop between max. point and reattached point and a small time lag between static stall point and max. point. For the dynamic moment curve, there is also a relatively smooth drop between max. point and reattached point. The max. point before the static point is caused by the dynamic moment shape that follows the trend of static curve before the static stall point.Figure 2
Moderate dynamic stall simulated curve and static experimental data. (a) Lift force coefficient and (b) moment coefficient.
(a)
(b)The pitch and plunge responses at moderate stall flutter are shown in Figure3, where the responses tend to damp at first 30 seconds, and they are stabilized after 40 seconds with asymmetric LCOs (e.g., plunge displacement LCO with magnitude from 0 m to −0.5 m and pitch angle LOC with magnitude from −5 deg. to −10 deg.). This phenomenon is convinced by the experimental results in [1], which indicated the occurrence of only asymmetric LCOs with either positive or negative pitching centering.Figure 3
Pitch and plunge responses at moderate stall flutter.


## 3.2. Deep Dynamic Stall

The deep stall aerodynamic curves are also simulated based on static experimental data. Figure4 shows the deep stall motions of the simulated unsteady lift and moment coefficients at angle of attack α
=
13.5
°
+
3.5
°
s
i
n
⁡
(
10
t
). The deep stall curves are also well predicted by the proposed model to capture the trend of static curves in the neighborhood as expected. In contrast with the moderate stall in Figure 2, the characteristics of deep stall dynamic are revealed by a relatively sharp drop between max. point and reattached point due to the strong separation of the flow and more drastic hysteresis with significant difference between up-stroke and down-stroke phase of dynamic curves. In addition, the maximum lift and moment point are far beyond the static point with a large time lag.Figure 4
Deep dynamic stall simulated curve and static experimental data. (a) Lift force coefficient and (b) moment coefficient.
(a)
(b)The stall phenomena in Figures2 and 4 can also be confirmed by the model validation results in [26], where the lift loop near max. point has a very different shape owing to the separation effects, and the lift curve has expected overshoot, while moment curves revolve around the static curves.Figure5 shows the pitch and plunge responses at deep stall flutter, where the asymmetric LCOs are also observed and their magnitudes are invariant after 25 seconds, resulting in magnitude of the plunge from 0 m to −1 m and magnitude of the pitch from −5 deg. to −15 deg. Compared with the moderate stall flutter, the deep stall can further decrease the damping ratio of the pitch and plunge modes, leading to relatively larger LCOs in a shorter time.Figure 5
Pitch and plunge responses at deep stall flutter.Theθ
-
θ
˙ phase plane is illustrated in Figure 6 for both moderate stall flutter case and deep stall flutter case. Two asymmetric limit cycle oscillations are observed. It is revealed that the asymmetric LCOs grow in size, and pitching center moves towards the left half of the plane as the section goes through more significant stall hysteresis, which can bring severe oscillations. This phase analysis result is consistent with the results of time-domain responses in Figures 3 and 5.Figure 6
Phase plane ofθ and θ
˙ of the system at moderate/deep stall flutter.


## 4. Stall Flutter Control

The purpose of this section is to design the effective controller to suppress the limit cycle oscillations of both moderate stall flutter and deep stall flutter, discussed in the previous section. Since the stall flutter occurs with the unsteady time-varying aerodynamic forces, which is induced by the time-varying change of angle of attackα
0
(
t
), the aeroservoelastic system in (17) can be rewritten as (19)
X
˙
=
A
t
X
+
B
u
y
=
C
X
,where A
(
t
) is the state matrix, which holds the coefficients in (17), u is the control input, that is, the microtab deployment of the actuator, and the output sensor signal y can be determined by the output matrix C.


### 4.1. The Test of Former Aeroelastic Controllers

Before the study of the control algorithm for stall flutter control, two former controllers that have been developed by the author in classic flutter control study [25] and pitch control study [24] are first tested on the stall flutter case. Those controllers are also developed with the actuator of the microtab. The aim is to investigate the effectiveness of two former aeroelastic controllers on the new control case, that is, the stall flutter control.Controller One. It is the classic flutter controller u
=
G
e
y, y
=
θ
˙. This controller is an output feedback controller, which was shown to be effective in classic flutter suppression of the airfoil section.Controller Two. It is the adaptive pitch controller u
=
G
e
y (y
=
θ
+
0.1
θ
˙, G
˙
e
=
-
γ
e
y
2). This adaptive control was proved to have full controllability of pitch mode and partial controllability of plunge mode.Those two aeroelastic controllers are adopted to control the stall flutter of the blade section separately. The control results are shown in Figure7 for moderate stall case and in Figure 8 for deep stall case. It can be found that the stall flutter oscillations are not effectively suppressed via the adaptive pitch controller in 5 seconds even though there is tendency of decaying for both stall cases. At the same time, even though the classic flutter controller works for the pitch mode and performs better than the adaptive pitch controller on the plunge mode, the oscillations of the plunge displacement cannot be totally suppressed in 5 seconds as they do for classic flutter case. The better control performance of pitching than plunging is due to control design with system output of pitching signal, as indicated in controller one and controller two. Therefore, those former flutter controllers are not good enough to eliminate the stall flutter effectively in a short time. Stall flutter controller needs to be designed with the investigation of new sensor signal and new control algorithm.Figure 7
System responses with classic flutter controller and adaptive pitch controller at moderate stall flutter case.Figure 8
System responses with classic flutter controller and adaptive pitch controller at deep flutter case.


### 4.2. The Design of Stall Flutter Controller

In order to develop the effective stall flutter controller for both stall cases, there are several aspects that need to be investigated: first, the effective sensor signals (or the control input signals for the controller) that can control the aeroelastic behavior of the blade section undergoing dynamic stall significantly, second, the controllability of the microtab (the control capability of the actuator), and third, the powerful control algorithm that can suppress the asymmetric LCOs of the stall flutter in a short time.


#### 4.2.1. Sensor Signal Effect on Stall Flutter Control

To study the sensor signal effect on the control of stall flutter, the eigenvalues of the closed-loop system with a static output feedback of the sensor signal are analyzed. The state matrix of the closed-loop system is taken asA
c
=
A
+
B
G
e
C, A
=
A
t
|
α
0
t
=
11, for moderate stall case and A
=
A
t
|
α
0
t
=
13.5 for deep stall case; y
=
C
q is the sensor signal. Note that the controllability matrix can also be computed to indicate the sensor signal effect, but here the eigenvalue analysis makes a more vivid explanation. The LCOs of the stall flutter indicate the imaginary stability of open-loop system.For moderate stall flutter case, Figure9 illustrates the eigenvalues of plunge and pitch modes in closed-loop systems with different sensor signals. It is revealed that the sensor of h, the sensor of θ, and the sensor of θ
˙ all lose the controllability of the plunge mode in stall flutter, even though they can have the full controllability of the pitch mode. Nevertheless, both the sensor of h
˙ and the sensor of h
˙
+
θ
˙ can suppress the LCOs with no imaginary parts in plunge mode and all negative values in pitch mode. Among those two sensor signals, the sensor of h
˙
+
θ
˙ is more influential in control of moderate stall flutter. Note that, for consideration of the sensor number limitation or sensor cost in practice, the single sensor signal of h
˙ can be chosen for implementation. For deep stall flutter case, the eigenvalues results shown in Figure 10 are similar with those for moderate stall flutter case with little less negative values in pitch mode.Figure 9
Closed-loop eigenvalues of the system with different sensor signals for moderate stall case.Figure 10
Closed-loop eigenvalues of the system with different sensor signals for deep stall case.


#### 4.2.2. Microtab Control Analysis on Stall Flutter

Since the potential sensor signals are studied for both stall flutter cases, the limitation of the microtab deployment also needs to be considered for practice. The airfoil type is DU97W300-10, and the location of the microtab is at 100% chord. The deployment of the microtab is saturated to 1%–5% of chord. The aerodynamic coefficients are adopted with respect to the deployment of the microtab. The static test curves of aerodynamic coefficients without/with microtab and with different microtab deployments are detailed in [29]. In this section, four different limitations of the microtab deployment are tested to verify the control capability of the microtab on LCOs of moderate/deep stall flutter.The damping ratios of plunge and pitch modes in the feedback system with microtab limitations (1% of c, 2% of c, 4% of c, and 5% of c) are illustrated in Figure11. It is clear that both modes have increasing damping ratios (a sharp increase in plunge damping ratio and a moderate increase in pitch damping ratio) with the growing microtab deployments, indicating that the LCOs in both stall flutter cases can be successfully suppressed within the limitations of the microtab. And the microtab has more control ability on plunge mode than pitch mode.Figure 11
Damping ratios versus microtab limitation for moderate/deep stall flutter cases.For the control of plunge motion, small microtab limitation (e.g., 1-2% of c) leads to the decaying LCOs with small damping ratios less than one, and full microtab deployment (5% of c) brings completely damped responses with damping ratio of 1 for both stall flutter cases. However, the microtab deployment of 4% of c has a weaker impact on deep stall flutter with less damping ratio than that of moderate stall flutter case, revealing that it is a little harder to suppress the deep stall flutter compared to moderate stall flutter. For the control of pitch motion, although there is a small increase in damping ratio, the microtab seems to have similar effect on the increase of the damping ratio for both stall flutter cases. Thus, microtab is proved as the potential and effective actuator for the stall flutter control, with powerful deployment of 4% of c for moderate stall case and 5% of c deployment for deep stall case.


#### 4.2.3. Stall Flutter Controller

The control objective is to minimize the asymmetric LCOs in the stall flutter. Since the former aeroelastic controllers do not have capability in stall flutter control, the controller needs to be improved with new control strategies or new effective sensor signal for the stall flutter application. In this section, two different stall flutter controllers are developed based on different control algorithms and the new potential sensor signal for stall flutter control.Simple Stall Flutter Controller. First, the simple output feedback control strategy used for classic flutter control is adopted to design the stall flutter controller along with the new potential sensor signal in Section 4.2.1. Then, the simple stall flutter controller is given as (20)
u
=
G
e
y
,
y
=
h
˙
+
θ
˙
.Note that, on this condition, the stall flutter control system is taken as a single input and single output system.Adaptive Stall Flutter Controller. To improve the control capability of each vibration mode, the stall flutter control system is treated as a single-input two-output system, and the adaptive control is used as the control algorithm. The block diagram of the adaptive stall flutter control strategy is illustrated in Figure 12. Here, the output of the system is defined as the vector y
=
h
˙
θ
˙
T. Since the system is time-varying, the control gain should also be designed to accommodate with the time-varying change of system parameters.Figure 12
Block diagram of the adaptive stall flutter control strategy.The adaptive stall flutter controller is presented as(21)
u
=
G
e
t
y
,
G
e
t
=
g
1
t
g
2
t
,where G
e
(
t
) is the time-varying adaptive control gain, which can be updated by the adaptive law as (22)
g
1
t
=
r
1
h
˙
2
+
r
2
h
˙
θ
˙
,
g
2
t
=
r
1
h
˙
θ
˙
+
r
2
θ
˙
2
,where r
1 and r
2 are scalar. The adaptive law can regulate the change of adaptive control gain to make plunge displacement h and pitch angle θ asymptotically stable.Since the value ofr
1 and r
2 can impact the performance of the adaptive law, four strategies of r
1 and r
2 are studied for the stall flutter control: case 1, r
1
=
-
1 and r
2
=
-
1, case 2, r
1
=
-
1 and r
2
=
-
0.1, case 3, r
1
=
-
0.1 and r
2
=
-
1, and case 4, r
1
=
-
0.5 and r
2
=
-
0.5. First, the control responses of adaptive stall flutter controllers with case 1, case 2, and case 3 are compared in Figure 13. It can be revealed that the controller with all control strategy cases can suppress the LCOs of plunge and pitch modes in 5 seconds as expected. However, the control strategy of case 2 leads to high-frequency pitch oscillations in first 2 seconds for moderate stall case and even as long as 4 seconds for deep stall case because case 2 lacks control power in pitch mode with r
1
≪
r
1. At the same time, case 3 performs well in pitch oscillation suppression, while it performs worse in plunge mode among all three cases because of r
1
≪
r
1. Only case 1 which has equivalent weigh on both modes (r
1
=
r
2) can suppress all LCOs effectively. The unsatisfying control results of case 2 and case 3 are due to the fact that there is not very strong coupling between pitching and plunging under dynamic stall condition. Hence, the individual control of either plunging or pitching cannot suppress them all effectively at the same time. In order to further investigate the effect of r
1 and r
2, the controller with case 4 is compared with that with case 1 in Figure 14. It is shown that the decrease of value of r
1
,
r
1 brings the similar overshoot of the control response but with a longer settling time, especially for the plunge control. All in all, the adaptive stall flutter control with control strategy of case 1 is the most potential one for LCOs suppression for both stall flutter cases.Figure 13
Adaptive stall flutter control responses with different control strategies for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter caseFigure 14
Adaptive stall flutter control responses with different equal-weigh control strategies for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter case


### 4.3. Control Results

To show the superiority of the proposed adaptive stall flutter controller, it is compared with the simple stall flutter controller in (20). The results of control performance with those controllers are shown in Figure 15. In addition, the effects of adaptive stall flutter control on moderate stall flutter and deep stall flutter are illustrated in Figure 16.Figure 15
Response of the closed-loop system with simple stall flutter controller and adaptive stall flutter controller for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter caseFigure 16
Control results of the closed-loop system for moderate/deep stall flutter cases.


#### 4.3.1. Simple Stall Flutter Controller versus Adaptive Stall Flutter Controller

Figure15 shows the control results of the simple flutter controller and the adaptive stall flutter controller. With respect to the classic flutter controller in Section 4.1, the simple flutter controller has improved control performance on plunge mode with the new sensor signal. Specifically, the plunge displacement has a quick decay of the oscillations at the first several simulation intervals with a relatively short settling time and small magnitude of the oscillations, which confirms the effectiveness of the new sensor signal in the stall flutter control. In other words, the sensor signal has a significant impact on the improvement of the control effect on different flutter cases even with the same control algorithm.Although the simple flutter controller can effectively suppress the LCOs and stabilize the plunged mode in about 4 seconds and the pitch mode at the same time, the dynamic performance is not very satisfying with lots of small oscillations around the steady state after 1 second, which indicate that the utilization of only new sensor signal is not powerful enough to eliminate all vibrations in a short time for both modes. In contrast, the adaptive stall flutter controller shows the superiority with the new control algorithm and the new sensor signal. As indicated in Figure15, for moderate stall case, the proposed adaptive stall flutter controller greatly improves the control performance of pitch and plunge modes. Both vibration responses can converge to the steady state in 2 seconds, and there are no significant oscillations after 1.5 seconds to show the great dynamic performance. However, the static stall flutter controller is not capable of suppressing all oscillations in a short time even though it can stabilize the system with the decaying responses for all stall cases. For deep stall flutter with stronger dynamic stall, simple stall flutter controller performs even worse with larger oscillations in vibration modes, while the adaptive stall flutter control can still suppress the stall flutter quickly with good dynamic performance. Thus, the adaptive stall flutter controller outperforms the static stall flutter controller.For the physical meaning, first, the poor control performance of the classic flutter controller on pitching control is due to the weak coupling between pitching and plunging on stall flutter cases. Then, this weak coupling is strengthened at specific operation point by the simple stall flutter controller with new sensor signal, which is the most potential sensor signal in stall flutter control. However, as the parameters of the system are varying from one operation point to another, the coupling dynamic is also changing from time to time. On this condition, the simple stall flutter controller cannot adapt to the coupling variance with the constant control gain. Nevertheless, the proposed adaptive stall flutter controller can adapt to the varying coupling conditions at different operation points and adjust the control input by the adaptive law. Overall, there is improvement of the proposed adaptive stall flutter controller on the following aspects: (1) involving both of plunge acceleration and pitch rate in system output signal instead of only pitch rate for the output feedback control, (2) adding the cross term ofh
˙
θ
˙in the adaptive control law in (22), and (3) distributing the equal weigh on the control of plunged mode and pitch mode with r
1
=
-
1 and r
2
=
-
1 in (22). All those new control strategies can strengthen the coupling between pitching and plunging along the whole system trajectory and further improve the overall control performance. Hence, the proposed adaptive stall flutter controller is the effective and robust controller for the stall flutter suppression.


#### 4.3.2. Moderate Stall Flutter Control versus Deep Stall Flutter Control

Figure16 illustrates the control responses of the adaptive stall flutter control of both moderate and deep stall flutter cases. For both stall flutter cases, the proposed adaptive stall flutter controller has the full controllability of both vibration modes, and it can suppress all LCOs with good dynamic performance in a short time within the saturation of the microtab.However, the control effect in moderate stall flutter seems a little better than that of deep stall flutter. It can be found clearly that, compared with the moderate flutter case, the system responses of the deep stall flutter system tend to have a relatively larger overshoot, larger steady state error, and longer settling time, and the system needs more actuation of the microtab at around 1 second, when the system responses begin to converge. Therefore, it is more difficult to suppress the deep stall flutter than the moderate stall flutter. Note that the responses of the microtab deployment can quickly converge to zero after 1.5 seconds, which also confirms the fast actuation and the control capability of the microtab.The control results can be verified by the results of microtab control analysis on stall flutter. As indicated in Figure11, the damping ratio of both plunge and pitch modes in stall flutter can be increased by the actuator of microtab with more significant improvement of plunge damping ratio than pitch damping ratio, which confirms the better control performance of plunge displacement with more smooth response in Figure 16. At the same time, Figure 11 shows that damping ratios for moderate stall flutter can be improved in a better way than that for deep stall flutter, which verifies the conclusion of control results demonstrating that deep stall flutter control is tougher than the moderate stall flutter control.The control results in Figure16 are reasonable and can also be confirmed by physical explanation. Since the sharp change in aerodynamic loads for deep stall case leads to severe change of the coupling of vibration modes, the controller is more reluctant to the rapid change of the system dynamic and it will definitely need more power from the actuator to bond the strong coupling between vibration modes and gain the steady conditions.


## 4.1. The Test of Former Aeroelastic Controllers

Before the study of the control algorithm for stall flutter control, two former controllers that have been developed by the author in classic flutter control study [25] and pitch control study [24] are first tested on the stall flutter case. Those controllers are also developed with the actuator of the microtab. The aim is to investigate the effectiveness of two former aeroelastic controllers on the new control case, that is, the stall flutter control.Controller One. It is the classic flutter controller u
=
G
e
y, y
=
θ
˙. This controller is an output feedback controller, which was shown to be effective in classic flutter suppression of the airfoil section.Controller Two. It is the adaptive pitch controller u
=
G
e
y (y
=
θ
+
0.1
θ
˙, G
˙
e
=
-
γ
e
y
2). This adaptive control was proved to have full controllability of pitch mode and partial controllability of plunge mode.Those two aeroelastic controllers are adopted to control the stall flutter of the blade section separately. The control results are shown in Figure7 for moderate stall case and in Figure 8 for deep stall case. It can be found that the stall flutter oscillations are not effectively suppressed via the adaptive pitch controller in 5 seconds even though there is tendency of decaying for both stall cases. At the same time, even though the classic flutter controller works for the pitch mode and performs better than the adaptive pitch controller on the plunge mode, the oscillations of the plunge displacement cannot be totally suppressed in 5 seconds as they do for classic flutter case. The better control performance of pitching than plunging is due to control design with system output of pitching signal, as indicated in controller one and controller two. Therefore, those former flutter controllers are not good enough to eliminate the stall flutter effectively in a short time. Stall flutter controller needs to be designed with the investigation of new sensor signal and new control algorithm.Figure 7
System responses with classic flutter controller and adaptive pitch controller at moderate stall flutter case.Figure 8
System responses with classic flutter controller and adaptive pitch controller at deep flutter case.


## 4.2. The Design of Stall Flutter Controller

In order to develop the effective stall flutter controller for both stall cases, there are several aspects that need to be investigated: first, the effective sensor signals (or the control input signals for the controller) that can control the aeroelastic behavior of the blade section undergoing dynamic stall significantly, second, the controllability of the microtab (the control capability of the actuator), and third, the powerful control algorithm that can suppress the asymmetric LCOs of the stall flutter in a short time.


### 4.2.1. Sensor Signal Effect on Stall Flutter Control

To study the sensor signal effect on the control of stall flutter, the eigenvalues of the closed-loop system with a static output feedback of the sensor signal are analyzed. The state matrix of the closed-loop system is taken asA
c
=
A
+
B
G
e
C, A
=
A
t
|
α
0
t
=
11, for moderate stall case and A
=
A
t
|
α
0
t
=
13.5 for deep stall case; y
=
C
q is the sensor signal. Note that the controllability matrix can also be computed to indicate the sensor signal effect, but here the eigenvalue analysis makes a more vivid explanation. The LCOs of the stall flutter indicate the imaginary stability of open-loop system.For moderate stall flutter case, Figure9 illustrates the eigenvalues of plunge and pitch modes in closed-loop systems with different sensor signals. It is revealed that the sensor of h, the sensor of θ, and the sensor of θ
˙ all lose the controllability of the plunge mode in stall flutter, even though they can have the full controllability of the pitch mode. Nevertheless, both the sensor of h
˙ and the sensor of h
˙
+
θ
˙ can suppress the LCOs with no imaginary parts in plunge mode and all negative values in pitch mode. Among those two sensor signals, the sensor of h
˙
+
θ
˙ is more influential in control of moderate stall flutter. Note that, for consideration of the sensor number limitation or sensor cost in practice, the single sensor signal of h
˙ can be chosen for implementation. For deep stall flutter case, the eigenvalues results shown in Figure 10 are similar with those for moderate stall flutter case with little less negative values in pitch mode.Figure 9
Closed-loop eigenvalues of the system with different sensor signals for moderate stall case.Figure 10
Closed-loop eigenvalues of the system with different sensor signals for deep stall case.


### 4.2.2. Microtab Control Analysis on Stall Flutter

Since the potential sensor signals are studied for both stall flutter cases, the limitation of the microtab deployment also needs to be considered for practice. The airfoil type is DU97W300-10, and the location of the microtab is at 100% chord. The deployment of the microtab is saturated to 1%–5% of chord. The aerodynamic coefficients are adopted with respect to the deployment of the microtab. The static test curves of aerodynamic coefficients without/with microtab and with different microtab deployments are detailed in [29]. In this section, four different limitations of the microtab deployment are tested to verify the control capability of the microtab on LCOs of moderate/deep stall flutter.The damping ratios of plunge and pitch modes in the feedback system with microtab limitations (1% of c, 2% of c, 4% of c, and 5% of c) are illustrated in Figure11. It is clear that both modes have increasing damping ratios (a sharp increase in plunge damping ratio and a moderate increase in pitch damping ratio) with the growing microtab deployments, indicating that the LCOs in both stall flutter cases can be successfully suppressed within the limitations of the microtab. And the microtab has more control ability on plunge mode than pitch mode.Figure 11
Damping ratios versus microtab limitation for moderate/deep stall flutter cases.For the control of plunge motion, small microtab limitation (e.g., 1-2% of c) leads to the decaying LCOs with small damping ratios less than one, and full microtab deployment (5% of c) brings completely damped responses with damping ratio of 1 for both stall flutter cases. However, the microtab deployment of 4% of c has a weaker impact on deep stall flutter with less damping ratio than that of moderate stall flutter case, revealing that it is a little harder to suppress the deep stall flutter compared to moderate stall flutter. For the control of pitch motion, although there is a small increase in damping ratio, the microtab seems to have similar effect on the increase of the damping ratio for both stall flutter cases. Thus, microtab is proved as the potential and effective actuator for the stall flutter control, with powerful deployment of 4% of c for moderate stall case and 5% of c deployment for deep stall case.


### 4.2.3. Stall Flutter Controller

The control objective is to minimize the asymmetric LCOs in the stall flutter. Since the former aeroelastic controllers do not have capability in stall flutter control, the controller needs to be improved with new control strategies or new effective sensor signal for the stall flutter application. In this section, two different stall flutter controllers are developed based on different control algorithms and the new potential sensor signal for stall flutter control.Simple Stall Flutter Controller. First, the simple output feedback control strategy used for classic flutter control is adopted to design the stall flutter controller along with the new potential sensor signal in Section 4.2.1. Then, the simple stall flutter controller is given as (20)
u
=
G
e
y
,
y
=
h
˙
+
θ
˙
.Note that, on this condition, the stall flutter control system is taken as a single input and single output system.Adaptive Stall Flutter Controller. To improve the control capability of each vibration mode, the stall flutter control system is treated as a single-input two-output system, and the adaptive control is used as the control algorithm. The block diagram of the adaptive stall flutter control strategy is illustrated in Figure 12. Here, the output of the system is defined as the vector y
=
h
˙
θ
˙
T. Since the system is time-varying, the control gain should also be designed to accommodate with the time-varying change of system parameters.Figure 12
Block diagram of the adaptive stall flutter control strategy.The adaptive stall flutter controller is presented as(21)
u
=
G
e
t
y
,
G
e
t
=
g
1
t
g
2
t
,where G
e
(
t
) is the time-varying adaptive control gain, which can be updated by the adaptive law as (22)
g
1
t
=
r
1
h
˙
2
+
r
2
h
˙
θ
˙
,
g
2
t
=
r
1
h
˙
θ
˙
+
r
2
θ
˙
2
,where r
1 and r
2 are scalar. The adaptive law can regulate the change of adaptive control gain to make plunge displacement h and pitch angle θ asymptotically stable.Since the value ofr
1 and r
2 can impact the performance of the adaptive law, four strategies of r
1 and r
2 are studied for the stall flutter control: case 1, r
1
=
-
1 and r
2
=
-
1, case 2, r
1
=
-
1 and r
2
=
-
0.1, case 3, r
1
=
-
0.1 and r
2
=
-
1, and case 4, r
1
=
-
0.5 and r
2
=
-
0.5. First, the control responses of adaptive stall flutter controllers with case 1, case 2, and case 3 are compared in Figure 13. It can be revealed that the controller with all control strategy cases can suppress the LCOs of plunge and pitch modes in 5 seconds as expected. However, the control strategy of case 2 leads to high-frequency pitch oscillations in first 2 seconds for moderate stall case and even as long as 4 seconds for deep stall case because case 2 lacks control power in pitch mode with r
1
≪
r
1. At the same time, case 3 performs well in pitch oscillation suppression, while it performs worse in plunge mode among all three cases because of r
1
≪
r
1. Only case 1 which has equivalent weigh on both modes (r
1
=
r
2) can suppress all LCOs effectively. The unsatisfying control results of case 2 and case 3 are due to the fact that there is not very strong coupling between pitching and plunging under dynamic stall condition. Hence, the individual control of either plunging or pitching cannot suppress them all effectively at the same time. In order to further investigate the effect of r
1 and r
2, the controller with case 4 is compared with that with case 1 in Figure 14. It is shown that the decrease of value of r
1
,
r
1 brings the similar overshoot of the control response but with a longer settling time, especially for the plunge control. All in all, the adaptive stall flutter control with control strategy of case 1 is the most potential one for LCOs suppression for both stall flutter cases.Figure 13
Adaptive stall flutter control responses with different control strategies for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter caseFigure 14
Adaptive stall flutter control responses with different equal-weigh control strategies for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter case


## 4.2.1. Sensor Signal Effect on Stall Flutter Control

To study the sensor signal effect on the control of stall flutter, the eigenvalues of the closed-loop system with a static output feedback of the sensor signal are analyzed. The state matrix of the closed-loop system is taken asA
c
=
A
+
B
G
e
C, A
=
A
t
|
α
0
t
=
11, for moderate stall case and A
=
A
t
|
α
0
t
=
13.5 for deep stall case; y
=
C
q is the sensor signal. Note that the controllability matrix can also be computed to indicate the sensor signal effect, but here the eigenvalue analysis makes a more vivid explanation. The LCOs of the stall flutter indicate the imaginary stability of open-loop system.For moderate stall flutter case, Figure9 illustrates the eigenvalues of plunge and pitch modes in closed-loop systems with different sensor signals. It is revealed that the sensor of h, the sensor of θ, and the sensor of θ
˙ all lose the controllability of the plunge mode in stall flutter, even though they can have the full controllability of the pitch mode. Nevertheless, both the sensor of h
˙ and the sensor of h
˙
+
θ
˙ can suppress the LCOs with no imaginary parts in plunge mode and all negative values in pitch mode. Among those two sensor signals, the sensor of h
˙
+
θ
˙ is more influential in control of moderate stall flutter. Note that, for consideration of the sensor number limitation or sensor cost in practice, the single sensor signal of h
˙ can be chosen for implementation. For deep stall flutter case, the eigenvalues results shown in Figure 10 are similar with those for moderate stall flutter case with little less negative values in pitch mode.Figure 9
Closed-loop eigenvalues of the system with different sensor signals for moderate stall case.Figure 10
Closed-loop eigenvalues of the system with different sensor signals for deep stall case.


## 4.2.2. Microtab Control Analysis on Stall Flutter

Since the potential sensor signals are studied for both stall flutter cases, the limitation of the microtab deployment also needs to be considered for practice. The airfoil type is DU97W300-10, and the location of the microtab is at 100% chord. The deployment of the microtab is saturated to 1%–5% of chord. The aerodynamic coefficients are adopted with respect to the deployment of the microtab. The static test curves of aerodynamic coefficients without/with microtab and with different microtab deployments are detailed in [29]. In this section, four different limitations of the microtab deployment are tested to verify the control capability of the microtab on LCOs of moderate/deep stall flutter.The damping ratios of plunge and pitch modes in the feedback system with microtab limitations (1% of c, 2% of c, 4% of c, and 5% of c) are illustrated in Figure11. It is clear that both modes have increasing damping ratios (a sharp increase in plunge damping ratio and a moderate increase in pitch damping ratio) with the growing microtab deployments, indicating that the LCOs in both stall flutter cases can be successfully suppressed within the limitations of the microtab. And the microtab has more control ability on plunge mode than pitch mode.Figure 11
Damping ratios versus microtab limitation for moderate/deep stall flutter cases.For the control of plunge motion, small microtab limitation (e.g., 1-2% of c) leads to the decaying LCOs with small damping ratios less than one, and full microtab deployment (5% of c) brings completely damped responses with damping ratio of 1 for both stall flutter cases. However, the microtab deployment of 4% of c has a weaker impact on deep stall flutter with less damping ratio than that of moderate stall flutter case, revealing that it is a little harder to suppress the deep stall flutter compared to moderate stall flutter. For the control of pitch motion, although there is a small increase in damping ratio, the microtab seems to have similar effect on the increase of the damping ratio for both stall flutter cases. Thus, microtab is proved as the potential and effective actuator for the stall flutter control, with powerful deployment of 4% of c for moderate stall case and 5% of c deployment for deep stall case.


## 4.2.3. Stall Flutter Controller

The control objective is to minimize the asymmetric LCOs in the stall flutter. Since the former aeroelastic controllers do not have capability in stall flutter control, the controller needs to be improved with new control strategies or new effective sensor signal for the stall flutter application. In this section, two different stall flutter controllers are developed based on different control algorithms and the new potential sensor signal for stall flutter control.Simple Stall Flutter Controller. First, the simple output feedback control strategy used for classic flutter control is adopted to design the stall flutter controller along with the new potential sensor signal in Section 4.2.1. Then, the simple stall flutter controller is given as (20)
u
=
G
e
y
,
y
=
h
˙
+
θ
˙
.Note that, on this condition, the stall flutter control system is taken as a single input and single output system.Adaptive Stall Flutter Controller. To improve the control capability of each vibration mode, the stall flutter control system is treated as a single-input two-output system, and the adaptive control is used as the control algorithm. The block diagram of the adaptive stall flutter control strategy is illustrated in Figure 12. Here, the output of the system is defined as the vector y
=
h
˙
θ
˙
T. Since the system is time-varying, the control gain should also be designed to accommodate with the time-varying change of system parameters.Figure 12
Block diagram of the adaptive stall flutter control strategy.The adaptive stall flutter controller is presented as(21)
u
=
G
e
t
y
,
G
e
t
=
g
1
t
g
2
t
,where G
e
(
t
) is the time-varying adaptive control gain, which can be updated by the adaptive law as (22)
g
1
t
=
r
1
h
˙
2
+
r
2
h
˙
θ
˙
,
g
2
t
=
r
1
h
˙
θ
˙
+
r
2
θ
˙
2
,where r
1 and r
2 are scalar. The adaptive law can regulate the change of adaptive control gain to make plunge displacement h and pitch angle θ asymptotically stable.Since the value ofr
1 and r
2 can impact the performance of the adaptive law, four strategies of r
1 and r
2 are studied for the stall flutter control: case 1, r
1
=
-
1 and r
2
=
-
1, case 2, r
1
=
-
1 and r
2
=
-
0.1, case 3, r
1
=
-
0.1 and r
2
=
-
1, and case 4, r
1
=
-
0.5 and r
2
=
-
0.5. First, the control responses of adaptive stall flutter controllers with case 1, case 2, and case 3 are compared in Figure 13. It can be revealed that the controller with all control strategy cases can suppress the LCOs of plunge and pitch modes in 5 seconds as expected. However, the control strategy of case 2 leads to high-frequency pitch oscillations in first 2 seconds for moderate stall case and even as long as 4 seconds for deep stall case because case 2 lacks control power in pitch mode with r
1
≪
r
1. At the same time, case 3 performs well in pitch oscillation suppression, while it performs worse in plunge mode among all three cases because of r
1
≪
r
1. Only case 1 which has equivalent weigh on both modes (r
1
=
r
2) can suppress all LCOs effectively. The unsatisfying control results of case 2 and case 3 are due to the fact that there is not very strong coupling between pitching and plunging under dynamic stall condition. Hence, the individual control of either plunging or pitching cannot suppress them all effectively at the same time. In order to further investigate the effect of r
1 and r
2, the controller with case 4 is compared with that with case 1 in Figure 14. It is shown that the decrease of value of r
1
,
r
1 brings the similar overshoot of the control response but with a longer settling time, especially for the plunge control. All in all, the adaptive stall flutter control with control strategy of case 1 is the most potential one for LCOs suppression for both stall flutter cases.Figure 13
Adaptive stall flutter control responses with different control strategies for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter caseFigure 14
Adaptive stall flutter control responses with different equal-weigh control strategies for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter case


## 4.3. Control Results

To show the superiority of the proposed adaptive stall flutter controller, it is compared with the simple stall flutter controller in (20). The results of control performance with those controllers are shown in Figure 15. In addition, the effects of adaptive stall flutter control on moderate stall flutter and deep stall flutter are illustrated in Figure 16.Figure 15
Response of the closed-loop system with simple stall flutter controller and adaptive stall flutter controller for (a) moderate stall flutter case and (b) deep stall flutter case.
(a)
Moderate stall flutter case
(b)
Deep stall flutter caseFigure 16
Control results of the closed-loop system for moderate/deep stall flutter cases.


### 4.3.1. Simple Stall Flutter Controller versus Adaptive Stall Flutter Controller

Figure15 shows the control results of the simple flutter controller and the adaptive stall flutter controller. With respect to the classic flutter controller in Section 4.1, the simple flutter controller has improved control performance on plunge mode with the new sensor signal. Specifically, the plunge displacement has a quick decay of the oscillations at the first several simulation intervals with a relatively short settling time and small magnitude of the oscillations, which confirms the effectiveness of the new sensor signal in the stall flutter control. In other words, the sensor signal has a significant impact on the improvement of the control effect on different flutter cases even with the same control algorithm.Although the simple flutter controller can effectively suppress the LCOs and stabilize the plunged mode in about 4 seconds and the pitch mode at the same time, the dynamic performance is not very satisfying with lots of small oscillations around the steady state after 1 second, which indicate that the utilization of only new sensor signal is not powerful enough to eliminate all vibrations in a short time for both modes. In contrast, the adaptive stall flutter controller shows the superiority with the new control algorithm and the new sensor signal. As indicated in Figure15, for moderate stall case, the proposed adaptive stall flutter controller greatly improves the control performance of pitch and plunge modes. Both vibration responses can converge to the steady state in 2 seconds, and there are no significant oscillations after 1.5 seconds to show the great dynamic performance. However, the static stall flutter controller is not capable of suppressing all oscillations in a short time even though it can stabilize the system with the decaying responses for all stall cases. For deep stall flutter with stronger dynamic stall, simple stall flutter controller performs even worse with larger oscillations in vibration modes, while the adaptive stall flutter control can still suppress the stall flutter quickly with good dynamic performance. Thus, the adaptive stall flutter controller outperforms the static stall flutter controller.For the physical meaning, first, the poor control performance of the classic flutter controller on pitching control is due to the weak coupling between pitching and plunging on stall flutter cases. Then, this weak coupling is strengthened at specific operation point by the simple stall flutter controller with new sensor signal, which is the most potential sensor signal in stall flutter control. However, as the parameters of the system are varying from one operation point to another, the coupling dynamic is also changing from time to time. On this condition, the simple stall flutter controller cannot adapt to the coupling variance with the constant control gain. Nevertheless, the proposed adaptive stall flutter controller can adapt to the varying coupling conditions at different operation points and adjust the control input by the adaptive law. Overall, there is improvement of the proposed adaptive stall flutter controller on the following aspects: (1) involving both of plunge acceleration and pitch rate in system output signal instead of only pitch rate for the output feedback control, (2) adding the cross term ofh
˙
θ
˙in the adaptive control law in (22), and (3) distributing the equal weigh on the control of plunged mode and pitch mode with r
1
=
-
1 and r
2
=
-
1 in (22). All those new control strategies can strengthen the coupling between pitching and plunging along the whole system trajectory and further improve the overall control performance. Hence, the proposed adaptive stall flutter controller is the effective and robust controller for the stall flutter suppression.


### 4.3.2. Moderate Stall Flutter Control versus Deep Stall Flutter Control

Figure16 illustrates the control responses of the adaptive stall flutter control of both moderate and deep stall flutter cases. For both stall flutter cases, the proposed adaptive stall flutter controller has the full controllability of both vibration modes, and it can suppress all LCOs with good dynamic performance in a short time within the saturation of the microtab.However, the control effect in moderate stall flutter seems a little better than that of deep stall flutter. It can be found clearly that, compared with the moderate flutter case, the system responses of the deep stall flutter system tend to have a relatively larger overshoot, larger steady state error, and longer settling time, and the system needs more actuation of the microtab at around 1 second, when the system responses begin to converge. Therefore, it is more difficult to suppress the deep stall flutter than the moderate stall flutter. Note that the responses of the microtab deployment can quickly converge to zero after 1.5 seconds, which also confirms the fast actuation and the control capability of the microtab.The control results can be verified by the results of microtab control analysis on stall flutter. As indicated in Figure11, the damping ratio of both plunge and pitch modes in stall flutter can be increased by the actuator of microtab with more significant improvement of plunge damping ratio than pitch damping ratio, which confirms the better control performance of plunge displacement with more smooth response in Figure 16. At the same time, Figure 11 shows that damping ratios for moderate stall flutter can be improved in a better way than that for deep stall flutter, which verifies the conclusion of control results demonstrating that deep stall flutter control is tougher than the moderate stall flutter control.The control results in Figure16 are reasonable and can also be confirmed by physical explanation. Since the sharp change in aerodynamic loads for deep stall case leads to severe change of the coupling of vibration modes, the controller is more reluctant to the rapid change of the system dynamic and it will definitely need more power from the actuator to bond the strong coupling between vibration modes and gain the steady conditions.


## 4.3.1. Simple Stall Flutter Controller versus Adaptive Stall Flutter Controller

Figure15 shows the control results of the simple flutter controller and the adaptive stall flutter controller. With respect to the classic flutter controller in Section 4.1, the simple flutter controller has improved control performance on plunge mode with the new sensor signal. Specifically, the plunge displacement has a quick decay of the oscillations at the first several simulation intervals with a relatively short settling time and small magnitude of the oscillations, which confirms the effectiveness of the new sensor signal in the stall flutter control. In other words, the sensor signal has a significant impact on the improvement of the control effect on different flutter cases even with the same control algorithm.Although the simple flutter controller can effectively suppress the LCOs and stabilize the plunged mode in about 4 seconds and the pitch mode at the same time, the dynamic performance is not very satisfying with lots of small oscillations around the steady state after 1 second, which indicate that the utilization of only new sensor signal is not powerful enough to eliminate all vibrations in a short time for both modes. In contrast, the adaptive stall flutter controller shows the superiority with the new control algorithm and the new sensor signal. As indicated in Figure15, for moderate stall case, the proposed adaptive stall flutter controller greatly improves the control performance of pitch and plunge modes. Both vibration responses can converge to the steady state in 2 seconds, and there are no significant oscillations after 1.5 seconds to show the great dynamic performance. However, the static stall flutter controller is not capable of suppressing all oscillations in a short time even though it can stabilize the system with the decaying responses for all stall cases. For deep stall flutter with stronger dynamic stall, simple stall flutter controller performs even worse with larger oscillations in vibration modes, while the adaptive stall flutter control can still suppress the stall flutter quickly with good dynamic performance. Thus, the adaptive stall flutter controller outperforms the static stall flutter controller.For the physical meaning, first, the poor control performance of the classic flutter controller on pitching control is due to the weak coupling between pitching and plunging on stall flutter cases. Then, this weak coupling is strengthened at specific operation point by the simple stall flutter controller with new sensor signal, which is the most potential sensor signal in stall flutter control. However, as the parameters of the system are varying from one operation point to another, the coupling dynamic is also changing from time to time. On this condition, the simple stall flutter controller cannot adapt to the coupling variance with the constant control gain. Nevertheless, the proposed adaptive stall flutter controller can adapt to the varying coupling conditions at different operation points and adjust the control input by the adaptive law. Overall, there is improvement of the proposed adaptive stall flutter controller on the following aspects: (1) involving both of plunge acceleration and pitch rate in system output signal instead of only pitch rate for the output feedback control, (2) adding the cross term ofh
˙
θ
˙in the adaptive control law in (22), and (3) distributing the equal weigh on the control of plunged mode and pitch mode with r
1
=
-
1 and r
2
=
-
1 in (22). All those new control strategies can strengthen the coupling between pitching and plunging along the whole system trajectory and further improve the overall control performance. Hence, the proposed adaptive stall flutter controller is the effective and robust controller for the stall flutter suppression.


## 4.3.2. Moderate Stall Flutter Control versus Deep Stall Flutter Control

Figure16 illustrates the control responses of the adaptive stall flutter control of both moderate and deep stall flutter cases. For both stall flutter cases, the proposed adaptive stall flutter controller has the full controllability of both vibration modes, and it can suppress all LCOs with good dynamic performance in a short time within the saturation of the microtab.However, the control effect in moderate stall flutter seems a little better than that of deep stall flutter. It can be found clearly that, compared with the moderate flutter case, the system responses of the deep stall flutter system tend to have a relatively larger overshoot, larger steady state error, and longer settling time, and the system needs more actuation of the microtab at around 1 second, when the system responses begin to converge. Therefore, it is more difficult to suppress the deep stall flutter than the moderate stall flutter. Note that the responses of the microtab deployment can quickly converge to zero after 1.5 seconds, which also confirms the fast actuation and the control capability of the microtab.The control results can be verified by the results of microtab control analysis on stall flutter. As indicated in Figure11, the damping ratio of both plunge and pitch modes in stall flutter can be increased by the actuator of microtab with more significant improvement of plunge damping ratio than pitch damping ratio, which confirms the better control performance of plunge displacement with more smooth response in Figure 16. At the same time, Figure 11 shows that damping ratios for moderate stall flutter can be improved in a better way than that for deep stall flutter, which verifies the conclusion of control results demonstrating that deep stall flutter control is tougher than the moderate stall flutter control.The control results in Figure16 are reasonable and can also be confirmed by physical explanation. Since the sharp change in aerodynamic loads for deep stall case leads to severe change of the coupling of vibration modes, the controller is more reluctant to the rapid change of the system dynamic and it will definitely need more power from the actuator to bond the strong coupling between vibration modes and gain the steady conditions.


## 5. Conclusion

In this study, the control of asymmetric LCOs of the moderate/deep stall flutter is investigated. The moderate/deep dynamic stall conditions of the blade section are first simulated based on experimental data, and the corresponding moderate/deep stall flutters are analyzed separately. The adaptive stall flutter controller is designed and studied for effective moderate/deep stall flutter control. The superiority of the proposed controller is proved by the comparison with other aeroelastic controllers. The conclusions can be made based on the results as follows.(1) The proposed aeroelastic model is capable of predicting the different dynamic stall conditions of the section with the simulation results of dynamic lift owing to expected overshoot and dynamic moment revolving around the static curves. The stall flutter analysis shows that, compared with the moderate stall flutter, the deep stall flutter can further decrease the damping ratio of the pitch and plunge modes, leading to relatively large asymmetric LCOs in a shorter time.(2) The passive microtab control is rarely studied in the field of moderate/deep stall flutter. For the design of effective stall flutter controller via microtab, several aspects are investigated. The sensor signal study shows that the sensor ofh
˙
+
θ
˙ is more influential in control of stall flutter, and the signal sensor of h
˙ also works for LCOs suppression for both plunge and pitch mode. Moreover, microtab is proved as the potential and effective actuator for the stall flutter control, with a sharp increase of damping ratio in plunge mode and a moderate increase of damping ratio in pitch mode. Furthermore, the control strategy study reveals that the equivalent weigh on the gains in adaptive law of the adaptive stall flutter controller leads to more satisfactory control performances at both stall flutter cases.(3) The superiority of the proposed adaptive stall flutter controller is proved by the comparison with the simple stall flutter controller. The simulation results illustrate that the proposed stall flutter controller performs better in LCOs suppression with great dynamic performance. In addition, the deep stall flutter is more difficult to control than the moderate stall flutter.The future work includes the implementation of the proposed control method on the 3D wind turbine blade. For 3D geometries, the vibration control is more difficult and complicated since the vibration modes increase to three degrees of freedom motions. The aeroservoelastic model for the 3D wind turbine blade should be established for the control design. In order to solve the control problem for 3D blade, the adaptive controller should be redesigned for all three vibration modes, and new control strategies in the control law should be adopted to achieve the best control performance.



---

*Source: 5096128-2016-03-20.xml*