flutter-analysis

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:

1. Assume reduced frequency k (sweep over k = 0.01–1.0)
2. Solve eigenvalue problem: [K/q + i g_artificial] I - ω² M/q + A(k)] {x} = 0
3. For each k: extract eigenvalue λ = (ω/V)² and artificial damping g = 2ζ
4. 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:

1. Guess k; compute A(k)
2. Solve eigenvalue for p; extract new k = Im(p)×b/V
3. Iterate until k_assumed = k_extracted
4. 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).