FLUTTER ANALYSIS OF A FIXED-WING AIRCRAFT

Scope of the Project:

 Study of Flutter.
 Influence of mass-axis and flexural-axis on the flutter speed.
 Solving the Flutter equations.
 Methods of solution “Eigenvalue Solution and ‘k’-Method”.
 MATLAB code for Eigenvalue and K-Method.
 Analysis in Ansys FSI.
Literature Study:
 M.R.Moosavi,A.R.Naddaf Oskouei,A.Khelil “Flutter of Subsonic
Wing” Elsevier Thin-Walled Structures 43 (2005) 617-627.
 S.A.Fazelzadeh,A.Mazidi,H.Kalantari “Bending-torsional flutter
of wings with an attached mass subjected to a follower force”
Elsevier Journal of Sound and Vibrations 323(2009) 148-162.

Work Done:
 Taken General form of Aeroelastic Equations.
 Solving the Equations by Eigenvalue Solution and K-Method in
MATLAB.
 Mass and Flexural axis effect the flutter speed.

Work Done after 2nd Review:

 Effect of M𝜃̇ for Rigid and Flexible Wing.
 Analysis in Ansys by One-way FSI.
 Flutter speed without damping.
 Flutter Speed Prediction.
 Flutter Speed Prediction:
𝑎11 𝑎12 𝑞1̈ 𝑏11 𝑏12 𝑞1̇ 𝑐11 𝑐12 𝑒11 0 𝑞1 0
[ ]{ ̈ }+ 𝑉[ ] { } + (𝑉 ∗ 𝑉 [ ]+[ ]) { } = { }
𝑎21 𝑎22 𝑞2 𝑏21 𝑏22 𝑞2̇ 𝑐21 𝑐22 0 𝑒22 𝑞2 0

Assume a solution of the form

𝑞1 𝑞1
{ } = { } 𝑒 𝜆𝑡
𝑞2 𝑞2
and also make sub:

𝑥𝑉 2 = 𝑒11 , 𝑒22 = 𝜇𝑒11 = 𝜇𝑥𝑉 2

𝑎 𝜆2 + 𝑏11 𝑉𝜆 + (𝑐11 + 𝑥)𝑉 2 𝑎12 𝜆2 + 𝑏12 𝑉𝜆 + 𝑐12 𝑉 2

| 11 |=0
𝑎21 𝜆2 + 𝑏21 𝑉𝜆 + 𝑐21 𝑉 2 𝑎22 𝜆2 + 𝑏22 𝑉𝜆 + (𝑐22 + 𝜇𝑥)𝑉 2

Solving the determinant gives the quartic equation

𝑏4 𝜆4 + 𝑏3 𝜆3 + 𝑏2 𝜆2 + 𝑏1 𝜆 + 𝑏0 = 0

Where, 𝑏4 = 𝑎11 𝑎22 − 𝑎21 𝑎12 ,

𝑏3 = (𝑎11 𝑏22 + 𝑏11 𝑎22 − 𝑎21 𝑏12 − 𝑎12 𝑏21 )𝑉,

𝑏2 = [(𝜇𝑎11 + 𝑎22 )𝑥 + (𝑎11 𝑐22 + 𝑏11 𝑏22 + 𝑐11 𝑎22 − 𝑎21 𝑐12 − 𝑏12 𝑏21 − 𝑐21 𝑎12 )]𝑉 2 = (𝑝1 𝑥 + 𝑝0 )𝑉 2 ,

𝑏1 = [(𝜇𝑏11 + 𝑏22 )𝑥 + (𝑏11 𝑐22 + 𝑐11 𝑏22 − 𝑏21 𝑐12 − 𝑐21 𝑏12 )]𝑉 3 = (𝑞1 𝑥 + 𝑞0 )𝑉 3 ,

𝑏0 = [𝜇𝑥 2 + (𝑐22 + 𝜇𝑐11 )𝑥 + 𝑐11 𝑐22 − 𝑐12 𝑐21 ]𝑉 4 = (𝑟2 𝑥 2 + 𝑟1 𝑥 + 𝑟0 )𝑉 4

By sub the above equations in determinant equation we get:

(𝑏4 𝑞12 − 𝑏3 𝑞1 𝑝1 + 𝑏32 𝑟2 )𝑥 2 + (2𝑏4 𝑞1 𝑞0 − 𝑏3 𝑞0 𝑝1 − 𝑏3 𝑞1 𝑝0 + 𝑏32 𝑟1 )𝑥 + (𝑏4 𝑞02 − 𝑏3 𝑞0 𝑝0 + 𝑏32 𝑟0 ) = 0

By solving the about quartic equation we will get two critical flutter speeds

𝑉1 = 122.1236 𝑚⁄𝑠 𝑎𝑛𝑑 𝑉2 = 468.6894 𝑚⁄𝑠

 Mtheta Variation for Rigid wing
1

0.5

0
0 20 40 60 80 100 120 140 160 180 200
-0.5
Mtheta

-1

-1.5

-2

-2.5
Flutter Speed

Mtheta Variation for Flexible Wing

0.5

0
0 10 20 30 40 50 60 70 80 90
-0.5
Mtheta

-1

-1.5

-2

-2.5
Flutter Speed
 6

4
Freq (Hz)

0
0 50 100 150
Air Speed (m/s)

100
Damping Ratio (%)

-100

0 50 100 150
Air Speed (m/s)

Zero Damping Flexible Wing

12

10

8
Freq (Hz)

0
0 50 100 150 200 250 300
Air Speed (m/s)

40

30

20
Damping Ratio (%)

10

-10

-20

-30

-40
0 50 100 150 200 250 300
Air Speed (m/s)

Zero Damping Rigid Wing

 11

10

9
Freq (Hz)
8

4
0 20 40 60 80 100 120 140 160 180 200
Air Speed (m/s)

40

30
Damping Ratio (%)

20

10

-10
0 20 40 60 80 100 120 140 160 180 200
Air Speed (m/s)

Damping Rigid Wing:Flutter speed at 154m/s ,frequency of 8.2Hz

3
Freq (Hz)

0
0 50 100 150
Air Speed (m/s)

100

50
Damping Ratio (%)

-50

-100

-150
0 50 100 150
Air Speed (m/s)

Damping Flexible Wing:Flutter speed at 45m/s ,frequency of 4.5Hz

Flutter speed at 154m/s in ansys ,Frequency of 7.75Hz.

Name: Paramasivam Suresh

Reg No: UR13AE044
Project Guide: Mr. Daniel Antony

Signature of Student Signature of Guide