Hand Calculations logo Hand Calculations All help pages ▾

Dutch Roll Wn⚠ unverified

Aerospace / Stability · Compute the approximate dutch-roll natural frequency

Parameters

InputSymbolUnitDefaultDescription
Cn_betaCn_β1.0Yawing-moment coefficient derivative with respect to sideslip angle (per radian, non-dimensional)
qqPa1.0Dynamic pressure
SSm^21.0Reference wing area
IzzIzzkg*m^21.0Yaw-axis moment of inertia
bbm1.0Wing span
OutputSymbolUnitDescription
resultwnrad/sDutch-roll natural frequency, in radians per second (rad/s)

The science & history

Understanding the Parameters

Derivation (Approaching a Proof)

Model the dutch roll, to leading order, as a single-degree-of-freedom yaw oscillation — a torsional mass-spring system about the yaw axis. Newton's second law for rotation (Torque) about the yaw axis:

$$I_{zz}\,\ddot\psi = N,$$

where $\psi$ is the yaw (heading) perturbation and $N$ is the yawing moment. For small sideslip, the restoring yawing moment is the directional stiffness times the sideslip; taking the sideslip to track the yaw perturbation ($\beta \approx -\psi$ for this simplified mode), the aerodynamic moment is

$$N = q\,S\,b\,C_n = q\,S\,b\,C_{n\beta}\,\beta \approx -\,q\,S\,b\,C_{n\beta}\,\psi.$$

Substituting gives the equation of a simple harmonic oscillator:

$$I_{zz}\,\ddot\psi + \big(q\,S\,b\,C_{n\beta}\big)\psi = 0 \quad\Longrightarrow\quad \ddot\psi + \frac{q\,S\,b\,C_{n\beta}}{I_{zz}}\,\psi = 0.$$

The natural frequency of $\ddot\psi + \omega_n^2\psi = 0$ is $\omega_n = \sqrt{k/I}$ with the effective stiffness $k = q\,S\,b\,C_{n\beta}$:

$$\omega_n = \sqrt{\frac{q\,S\,b\,C_{n\beta}}{I_{zz}}}.$$

The registry carries an extra factor of $\tfrac12$ inside the root (a convention/normalisation choice, giving $\sqrt{C_{n\beta}qSb/(2I_{zz})}$). The damping would come from the yaw-rate derivative $C_{nr}$ and the roll coupling — a second equation this single-DOF model omits, which is why only the frequency emerges here.

Dimensional check. $C_{n\beta}$ is per radian (treated as dimensionless), $q\,S\,b = \text{Pa}\cdot\text{m}^2 \cdot\text{m} = \text{N}\cdot\text{m}$ (a moment), and $I_{zz} = \text{kg}\cdot\text{m}^2$: $$\frac{q\,S\,b}{I_{zz}} = \frac{\text{N}\cdot\text{m}}{\text{kg}\cdot\text{m}^2} = \frac{\text{kg}\cdot\text{m}^2/\text{s}^2}{\text{kg}\cdot\text{m}^2} = \frac{1}{\text{s}^2},$$ and the square root gives $\text{s}^{-1} = \text{rad/s}$ ✓.

History and Development

Related Concepts: Weathercock Stability, Dihedral Effect, Spiral Mode Time Constant, Roll Mode Time Constant, Phugoid Damping, Natural Frequency mass-spring, Damping Ratio, Dynamic Pressure

Notes: Dutch roll = coupled yaw–roll wallowing oscillation (lateral-directional analogue of the short period). Modelled here as a yaw-only mass-spring: stiffness $qSbC_{n\beta}$ (directional stability), inertia $I_{zz}$ → $\omega_n=\sqrt{qSbC_{n\beta}/I_{zz}}$ (registry adds a $\tfrac12$). Faster at high $q$ (low altitude/high speed). Gives frequency, not damping — the damping (the real problem, from $C_{nr}$ + roll coupling) is why jets carry yaw dampers. $C_{n\beta}$ must be $>0$.

← Back to the workspace  ·  All help pages  ·  Getting started