| name | flutter-analysis |
| description | Flutter analysis — classical flutter (bending-torsion), panel flutter, Theodorsen function, p-k method, V-g method, flutter speed, divergence, control surface flutter, FAR 25.629, aeroelastic margin, GVT correlation. |
| metadata | {"priority":7,"promptSignals":{"phrases":["flutter analysis","aeroelastic flutter","panel flutter","Theodorsen flutter","FAR 25.629","flutter speed prediction"],"minScore":3}} |
Flutter Analysis — Complete Skill
Flutter Fundamentals
Flutter: dynamic aeroelastic instability; coupling of aerodynamic, elastic, and inertial forces → growing oscillations
Types:
- Classical (bending-torsion) flutter: two modes couple through aerodynamic forces
- Panel flutter: thin panel in supersonic flow → limit-cycle oscillation
- Control surface flutter: control surface hinge moment coupling
- Stall flutter: nonlinear — occurs in stall conditions (blades, rotors)
- Buffet: forced response from separated flow (not true instability)
Flutter margin requirement (FAR 25.629):
U_flutter ≥ 1.15 × V_D (dive speed) [where V_D = design dive speed; V_flutter must have 15% margin above dive speed]
Must demonstrate freedom from flutter, divergence, control reversal at 1.15 × V_D
Theodorsen Theory (2D Airfoil)
2D Section flutter equations (pitch-plunge):
[M]{ẍ} + [C]{ẋ} + [K - q × A(k)]{x} = {0}
M = structural mass matrix; K = structural stiffness; q = dynamic pressure; A(k) = aerodynamic matrix
k = ωb/V = reduced frequency; b = semi-chord; V = airspeed; ω = flutter frequency
Theodorsen lift deficiency function C(k):
C(k) = F(k) + iG(k) = H₁⁽²⁾(k) / [H₁⁽²⁾(k) + iH₀⁽²⁾(k)] [H = Hankel functions of second kind]
F(k): real part; G(k): imaginary part
k → 0 (quasi-steady): C(k) → 1.0; k → ∞ (high freq): C(k) → 0.5
Aerodynamic lift and moment per unit span:
L = 2πρV²b × C(k) × (α + ḣ/V + b(0.5-a)α̇/V) + ρπb² × (V α̇ + ḧ - baα̈)
M = ... [full expression in Fung, "Introduction to Aeroelasticity"]
Flutter determinant:
det[K - ω² M - q A(k)] = 0 → complex equation; find q and ω satisfying it
Solution Methods
V-g (k) Method
Procedure:
- Assume reduced frequency k (sweep over k = 0.01–1.0)
- Solve eigenvalue problem: [K/q + i g_artificial] I - ω² M/q + A(k)] {x} = 0
- For each k: extract eigenvalue λ = (ω/V)² and artificial damping g = 2ζ
- Plot g vs. V for each mode; flutter occurs where g = 0
V-g plot interpretation:
Mode 1 (bending): g < 0 at all speeds → stable
Mode 2 (torsion): g = 0 at V = V_flutter → onset of flutter
Modes coalesce in frequency as flutter speed approached (frequency coalescence)
Limitation: V-g method gives exact result only at k corresponding to actual flutter point; must iterate
p-k Method (More Accurate)
Eigenvalue problem:
[p² M + p C + K - q A(p)] {x} = {0}
p = eigenvalue; p = σ + iω (complex frequency); σ < 0: stable; σ > 0: unstable (flutter)
Procedure:
- Guess k; compute A(k)
- Solve eigenvalue for p; extract new k = Im(p)×b/V
- Iterate until k_assumed = k_extracted
- Sweep V; find V where Re(p) > 0
More physical than V-g: p represents actual damping; correct at all speeds (not just flutter point)
Structural Model Requirements
GVT (Ground Vibration Test)
Required before flutter analysis:
- Measure mode shapes, frequencies, and modal damping for all critical modes
- Correlation criterion: MAC (Modal Assurance Criterion) ≥ 0.9 for critical modes
MAC = |{Φ_FEA}ᵀ{Φ_test}|² / (|Φ_FEA|² × |Φ_test|²)
- Frequency correlation: |f_FEA - f_test| / f_test ≤ 3% for critical modes
Critical modes for flutter:
Wing bending (1B), wing torsion (1T), antisymmetric bending, control surface rotation
Typically include 10–20 modes below 5 × flutter frequency
Structural Model (FEM → Flutter)
Normal modes extraction:
[K - ω² M]{Φ} = {0}; use first N modes as basis
Modal mass matrix: [m] = [Φ]ᵀ [M] [Φ] (diagonal if mass-normalized)
Modal stiffness: [k] = [ω²] (diagonal)
Generalized aerodynamic force (GAF) matrix:
[Q(k)] = [Φ]ᵀ [A_full(k)] [Φ] [computed from DLM or CFD]
Aerodynamic Methods
Doublet Lattice Method (DLM) — Linear Subsonic
Standard method for subsonic flutter (M < 0.7 practical limit):
Lifting surface divided into panels; doublet lattice singularities satisfy no-penetration BC
Computes [Q(k)] directly; implemented in NASTRAN (CAERO1/CAERO2 + PAERO cards)
NASTRAN flutter solution:
SOL 145 (FLUTTER = PKNL or PKS)
MKAERO1: define Mach numbers and reduced frequencies
FLFACT: airspeed sweep
Supersonic Panel Methods (M > 1.3)
Mach box method: panel aerodynamics using supersonic flow theory
ZONA51 / modified strip theory: captures spanwise effects at supersonic speeds
CFD-based GAF: unsteady RANS + structural modes → high-fidelity GAF matrix
Control Surface Flutter
Hinge moment coupling:
Control surface mass unbalance → aerodynamic hinge moment → flutter at lower speed than wing alone
Mass balance: add mass ahead of hinge line → reduce unbalanced moment
Static margin check:
x_cg_control_surface < x_hinge [CG of control surface must be at or forward of hinge]
Mass balance weight: m_balance = -I_unbalance / (x_balance - x_hinge)
Control surface frequency requirement:
ω_control_surface > ω_flutter_mode × 2.0 (decoupling criterion; prevents subcritical flutter coupling)
Panel Flutter (Supersonic)
Dynamic pressure parameter:
λ = q × a³ / D [D = panel bending stiffness = E × t³ / (12(1-ν²)); a = panel length in flow direction]
Critical dynamic pressure:
λ_cr = [μ/(μ + M × Im(k_flutter))] × (π/a)⁴ [complex; from panel flutter theory]
Simplified: λ_cr ≈ 500 (simply supported panel, M > 1.5, μ = mass ratio = ρ a/m_panel)
Limit cycle oscillation (LCO):
Post-flutter panel may oscillate at finite amplitude (nonlinear hardening) vs. catastrophic divergence
LCO amplitude depends on panel tension and geometric nonlinearity
Flutter Speed Estimate
Quick estimate — bending-torsion flutter:
V_f ≈ k_flutter × b × √(GJ / (ρ × chord³)) [k_flutter ≈ 0.5–0.8; depends on structural coupling]
Or: V_f ≈ V_divergence / √2 (Typical relation for classical flutter vs. divergence)
Divergence speed:
q_D = (dC_L/dα) × (e/S_tor) [e = aerodynamic-elastic axis offset; S_tor = torsional stiffness]
V_D = √(2 q_D / ρ)
Standards
| Standard | Scope |
|---|
| FAR 25.629 | Aeroelastic stability requirements (commercial transport) |
| MIL-A-8870 | Flutter, vibration, and divergence (military aircraft) |
| EASA CS-25.629 | European aeroelastic stability requirements |
| NASTRAN SOL 145 | NASTRAN flutter solution documentation |
| AIA/EUROCONTROL | Flutter test procedures |
Output
Provide: analysis method (Theodorsen/DLM/p-k/V-g), aircraft configuration (wing planform, t/c [%], sweep [°]), structural modes included (N modes; lowest 3 frequencies [Hz]), GVT correlation (MAC and frequency error [%]), aerodynamic method (DLM/Mach box/CFD-GAF), flutter speed V_flutter [KEAS] and flutter frequency ω_flutter [Hz], flutter mode description (bending-torsion coupling), divergence speed V_div [KEAS], design dive speed V_D [KEAS], flutter margin V_flutter / V_D (must be ≥ 1.15 per FAR 25.629), control surface mass balance status, and applicable standard (FAR 25.629, MIL-A-8870).
