Compute flutter and divergence speeds for a 2-DOF wing section (plunge h + pitch α) in real time. Explore V-g and V-omega diagrams and watch the animated wing cross-section show coupled oscillation behavior.
Wing Section Parameters
Chord length c
Bending frequency fh5.0 Hz
Torsion frequency fθ12.0 Hz
Air density ρ1.225 kg/m³
Mass ratio μ = m/(πρb²)20
Radius of gyration rα0.50
Eccentricity e/b0.20
Results
— m/s
Flutter speed Vf
— m/s
Divergence speed Vd
—
Reduced freq. k
—
Frequency ratio fh/fθ
Wing section animation — coupled plunge h and pitch α oscillation
V-g Diagram (damping vs speed)
V-g Diagram (Flutter Speed)
V-ω Diagram (frequency vs speed)
Vom
Theory & Key Formulas
\(V_f \approx \omega_\theta \cdot b \cdot r_\alpha \sqrt{\dfrac{\mu}{1+(2\pi/k_c)}}\)
\(V_{div}= \omega_\theta \cdot b \cdot \sqrt{\dfrac{\mu r_\alpha^2}{2\pi e/b}}\)
b = c/2 (semi-chord), k = ωb/V (reduced frequency)
What is Aeroelastic Flutter?
🙋
What exactly is "flutter"? I've heard it's a dangerous vibration in airplane wings.
🎓
Basically, it's a self-excited oscillation where the wing extracts energy from the airflow. In practice, the bending and twisting motions of the wing couple together, feeding each other. Above a critical speed—the flutter speed—the oscillations grow exponentially, leading to catastrophic failure in seconds. Try moving the "Bending Frequency" slider in the simulator to see how it immediately changes the predicted flutter speed.
🙋
Wait, really? So it's not just turbulence shaking the wing? What are these two speeds, flutter and divergence, calculated here?
🎓
Right, it's a fundamental instability of the structure itself. The simulator calculates two key speeds. The Flutter Speed (\(V_f\)) is where coupled oscillations become unstable. The Divergence Speed (\(V_{div}\)) is a static instability where the wing twists itself apart—like a strip of paper in strong wind. For instance, in the simulator, if you increase the "Eccentricity (e/b)", you'll see the divergence speed drop, making the wing more susceptible to twisting failure.
🙋
That makes sense. So what do all these parameters like "Mass Ratio (μ)" and "Radius of Gyration (rα)" actually represent for a real wing?
🎓
Great question! They describe the wing's physical properties. The Mass Ratio (μ) compares the wing's mass to the mass of air it displaces—a lightweight wing has a low μ. The Radius of Gyration (rα) tells us how the wing's mass is distributed for twisting. A common case is a wing with heavy engines on the tips, which increases rα. Play with the rα slider; you'll see a higher value generally increases the flutter speed, making the wing more stable against flutter but often more prone to divergence.
Physical Model & Key Equations
The simulator uses a classic 2-degree-of-freedom (2-DOF) model for a wing section, considering pure bending (plunging, h) and pure twisting (pitching, θ) motions. The interaction with aerodynamic forces leads to the characteristic instability speeds. The simplified estimate for the flutter speed is derived from equating the aerodynamic energy input to the structural damping.
$$V_f \approx \omega_\theta \cdot b \cdot r_\alpha \sqrt{\dfrac{\mu}{1+(2\pi/k_c)}}$$
Where:
• \(V_f\) = Flutter Speed
• \(\omega_\theta = 2\pi f_\theta\) = Torsional natural frequency (rad/s)
• \(b = c/2\) = Semi-chord length
• \(r_\alpha\) = Radius of gyration about the elastic axis
• \(\mu = m/(\pi\rho b^2)\) = Mass ratio
• \(k_c\) = A reduced frequency parameter related to the frequency ratio.
The divergence speed represents a static aeroelastic failure where the aerodynamic pitching moment overwhelms the wing's torsional stiffness. It occurs when the twist increases without bound for a small increase in airspeed.
$$V_{div}= \omega_\theta \cdot b \cdot \sqrt{\dfrac{\mu r_\alpha^2}{2\pi e/b}}$$
Where:
• \(V_{div}\) = Divergence Speed
• \(e\) = Distance between the aerodynamic center (lift point) and the elastic axis (twist axis).
• \(e/b\) = Non-dimensional eccentricity. A positive \(e\) (center ahead of axis) is destabilizing. This equation shows why moving the wing's spar (changing e/b) is a critical design choice.
Frequently Asked Questions
Flutter speed is the speed at which coupled bending-torsion vibrations diverge, determined by the point where the damping ratio crosses zero on the V-g diagram. Divergence speed is a static instability phenomenon determined solely by torsional stiffness, occurring when the aerodynamic moment exceeds the elastic restoring force. This tool automatically calculates both and displays markers on the graph.
The V-g diagram plots speed on the horizontal axis and damping ratio g on the vertical axis; instability (flutter) occurs when g exceeds 0. The V-ω diagram shows the change in frequency with speed, and flutter is more likely to occur near the region where bending and torsion modes approach or cross each other. By examining both graphs together, you can accurately determine the stability boundary.
Refer to actual wing design values or experimental data to input mass, stiffness, center of gravity position (static unbalance), and moment of inertia. In particular, a larger static unbalance strengthens the coupling and lowers the flutter speed. By setting appropriate aerodynamic coefficients (such as lift slope) according to the airfoil type, more realistic predictions are possible.
The animation visualizes the coupled behavior of wing cross-section bending (vertical motion) and torsion (rotational motion). As the speed changes, the phase difference and amplitude ratio of the vibrations vary, and near flutter, you can intuitively observe both modes diverging in phase. This helps in visually identifying design issues.
Real-World Applications
Aircraft Wing Design: This calculation is foundational. Engineers must ensure the operating speed (including a safety margin) is below both \(V_f\) and \(V_{div\). For instance, the mass ratio (μ) is minimized using composite materials to push the flutter speed higher, while the eccentricity (e/b) is carefully controlled by placing the main wing spar.
Wind Turbine Blades: Long, flexible blades are susceptible to flutter. CAE tools using this core model help optimize the internal spar layout and material distribution (affecting rα) to prevent instability during extreme gusts or operational shutdowns.
Formula 1 Front Wings: These wings are extremely lightweight and subject to high loads. Flutter analysis ensures the complex multi-element design does not develop oscillatory instabilities that could cause loss of downforce or structural failure at high speed.
Bridge Design (Tacoma Narrows Lesson): While the classic failure was torsional flutter, the principles are analogous. Modern bridge design uses similar aeroelastic stability checks to prevent wind-induced oscillations, often testing scaled sections in wind tunnels based on these equations.
Common Misconceptions and Points to Note
When you start using this tool, there are several points where beginners, especially those new to CAE, often stumble. A major misconception is thinking "the calculation result is the safety margin itself." This simulator is based on the most basic theoretical model for a 2D cross-section. An actual wing is a 3D structure with complex coupling of multiple modes. For example, the flutter speed \(V_f\) calculated here often requires a design margin of over 30% or more for a real aircraft. Consider the tool's results as a means to "understand trends and compare the effects of parameters."
Next, a pitfall in parameter setting: do not casually leave the density ratio \(\mu\) at its default value. This represents the "ratio of wing mass to air mass" and significantly reflects the design philosophy. For instance, flutter characteristics are entirely different between a lightweight glider (small \(\mu\)) and a heavy fighter jet (large \(\mu\)). In actual projects, the golden rule is to use an accurate \(\mu\) calculated from the air density at the anticipated flight altitude.
Finally, the idea that "you only need to look at static divergence or flutter" is dangerous. Try increasing the "eccentricity distance e/b" from 0.2 to 0.4 in the tool. You'll immediately see the divergence speed \(V_{div}\) drop sharply, indicating it can become critical before flutter. In practice, you must check both critical speeds for all flight conditions and focus your design on the lower one. This is the first step in thinking comprehensively about the "aeroelastic trio (flutter, divergence, control reversal)."
Set chord length (slCNum) between 0.5–3 m and elastic axis position (slC, as fraction of chord)
Input structural frequencies: fh (plunge in Hz) and fθ (pitch in Hz), typically 5–25 Hz for transport aircraft wings
Define air density (slRho, kg/m³—sea level ~1.225), mass per span (slMu, kg/m), radius of gyration (slRa, m), and Young's modulus (slE, GPa) for the wing material
Simulator calculates flutter speed Vf (m/s) and divergence speed Vd (m/s), plotting V-g and V-ω diagrams in real time
Worked Example
Transport wing section: chord 2.0 m, elastic axis at 0.3c, fh=8 Hz, fθ=12 Hz, air density 1.225 kg/m³, mass 180 kg/m, radius of gyration 0.6 m, aluminum (E=72 GPa). Reduced frequency k=(ωb)/V where b=chord/2=1.0 m. At sea level cruise (ρ=1.225), simulator yields flutter speed Vf≈145 m/s and divergence speed Vd≈215 m/s. Frequency ratio fh/fθ=0.67 produces typical coalescence behavior on V-g plot near Vf.
Practical Notes
Aeroelastic instability occurs when damping (g-curve) crosses zero; flutter speed Vf marks first coalescence of plunge and pitch modes—critical for certification above Vf
Divergence speed Vd (structural instability) typically exceeds Vf for typical wing mass distributions; if Vd approaches Vf, increase pitch stiffness (fθ)
Reduced frequency k varies inversely with airspeed; subsonic incompressible theory applies below Mach 0.3; for transonic wings, apply Prandtl-Mach correction: replace b with b/√(1−M²)