Formula sheet

Aircraft Flight Dynamics and Control Systems Formula Sheet

Aerospace engineering formula sheet for aircraft trim, static margin, load factor, turn performance, state-space models, modal damping, control authority, actuator limits, sampling, latency, envelope margins and validation checks.

This formula sheet collects first-pass relationships used to review aircraft flight dynamics and flight-control systems. It focuses on trim, stability, modes, control authority, sensors, actuators, sampling, latency, envelope margins and validation checks. Use it together with aerodynamic data, mass properties, control-law documentation, actuator limits, sensor calibration and flight-test evidence.

The equations below are screening and review tools. Detailed flight dynamics requires validated aerodynamic derivatives, configuration control, structural flexibility, atmosphere data, actuator models, sensor models, software timing and approved test procedures.

Notation

SymbolMeaningTypical unit
Vtrue airspeedm/s
\rhoair density\text{kg/m}^3
\bar{q}dynamic pressurePa
Waircraft weightN
maircraft masskg
Sreference wing area\text{m}^2
\bar{c}mean aerodynamic chordm
bwing spanm
I_x,I_y,I_zmass moments of inertia\text{kg m}^2
\alphaangle of attackrad
\betasideslip anglerad
\phibank anglerad
p,q,rroll, pitch and yaw ratesrad/s
C_L,C_D,C_m,C_l,C_nlift, drag, pitching, rolling and yawing coefficientsdimensionless
\delta_e,\delta_a,\delta_relevator, aileron and rudder deflectionrad
\omega_nnatural frequencyrad/s
\zetadamping ratiodimensionless
f_ssampling frequencyHz
T_dtotal delay or latencys

Dynamic Pressure and Aerodynamic Forces

Dynamic pressure:

\displaystyle \bar{q}=\frac{1}{2}\rho V^2

Lift:

L=\bar{q}SC_L

Drag:

D=\bar{q}SC_D

Pitching moment:

M=\bar{q}S\bar{c}C_m

Rolling moment:

L_r=\bar{q}SbC_l

Yawing moment:

N=\bar{q}SbC_n

Use

These equations are shared with aerodynamics, but flight dynamics uses them as inputs to motion, trim, control authority and envelope decisions. State configuration, Mach number, Reynolds number, center of gravity, mass and control deflections.

Trim Conditions

Level unaccelerated flight:

L\approx W
T\approx D

Pitching moment trim:

C_m=0

Linearized pitching moment model:

\displaystyle C_m=C_{m0}+C_{m_\alpha}\alpha+C_{m_{\delta_e}}\delta_e+C_{m_q}\frac{q\bar{c}}{2V}

For steady trim with q=0:

\displaystyle \delta_e=-\frac{C_{m0}+C_{m_\alpha}\alpha}{C_{m_{\delta_e}}}

Use

Trim depends on center of gravity, configuration, flap setting, thrust line, Mach number, Reynolds number, tail effectiveness and control limits. A trim equation is not a handling-quality assessment by itself.

Static Longitudinal Stability

Static margin:

\displaystyle SM=\frac{x_{np}-x_{cg}}{\bar{c}}

where x_{np} is neutral point and x_{cg} is center-of-gravity position measured in the same direction.

Simplified pitching-moment slope relation:

C_{m_\alpha}\approx -C_{L_\alpha}SM

Static longitudinal stability screen:

C_{m_\alpha}<0

Use

This simplified relation assumes consistent sign convention and a conventional small-disturbance interpretation. Static margin is not enough: control power, stall behavior, aeroelasticity and control-law mode can dominate the real aircraft.

Load Factor and Coordinated Turns

Load factor:

\displaystyle n=\frac{L}{W}

Level coordinated turn load factor:

\displaystyle n=\frac{1}{\cos\phi}

Turn rate:

\displaystyle \Omega=\frac{g\tan\phi}{V}

Turn radius:

\displaystyle R=\frac{V^2}{g\tan\phi}

Stall speed at load factor n:

V_{S,n}=V_{S,1}\sqrt{n}

Use

These relationships assume coordinated level turning and enough thrust, control authority and structural margin. They do not prove that the aircraft can sustain the turn.

Linear State-Space Model

Small-disturbance linear model:

\dot{x}=Ax+Bu

Output model:

y=Cx+Du

Longitudinal state example:

x=\begin{bmatrix}u&w&q&\theta\end{bmatrix}^T

Lateral-directional state example:

x=\begin{bmatrix}\beta&p&r&\phi\end{bmatrix}^T

Control input example:

u_c=\begin{bmatrix}\delta_e&\delta_a&\delta_r&T\end{bmatrix}^T

Use

State-space matrices are condition-specific. A matrix identified at one speed, altitude, mass, center of gravity, configuration or control-law mode should not be reused outside its evidence boundary without justification.

Second-Order Mode Approximation

Standard second-order denominator:

s^2+2\zeta\omega_ns+\omega_n^2=0

Damped natural frequency:

\omega_d=\omega_n\sqrt{1-\zeta^2}

Oscillation period:

\displaystyle T_d^{mode}=\frac{2\pi}{\omega_d}

Approximate 2 percent settling time:

\displaystyle t_s\approx\frac{4}{\zeta\omega_n}

Percent overshoot for a simple underdamped second-order response:

M_p=100e^{-\zeta\pi/\sqrt{1-\zeta^2}}

Use

Short-period, phugoid, Dutch-roll and structural modes are not always clean second-order systems, but this approximation helps screen damping, frequency and pilot or controller workload.

Damping from Flight-Test Peaks

Logarithmic decrement over N cycles:

\displaystyle \delta=\frac{1}{N}\ln\left(\frac{A_1}{A_{N+1}}\right)

Damping ratio:

\displaystyle \zeta=\frac{\delta}{\sqrt{4\pi^2+\delta^2}}

If measured period is T:

\displaystyle \omega_d=\frac{2\pi}{T}
\displaystyle \omega_n=\frac{\omega_d}{\sqrt{1-\zeta^2}}

Use

Use consistent peak measurement, filter settings and flight condition. Do not infer envelope clearance from one noisy decay trace without uncertainty, repeatability and configuration control.

Control Authority and Angular Acceleration

Pitch control moment from elevator:

M_{\delta_e}=\bar{q}S\bar{c}C_{m_{\delta_e}}\delta_e

Pitch angular acceleration:

\displaystyle \dot{q}\approx\frac{M_{\delta_e}}{I_y}

Roll control moment from aileron:

L_{\delta_a}=\bar{q}SbC_{l_{\delta_a}}\delta_a

Roll angular acceleration:

\displaystyle \dot{p}\approx\frac{L_{\delta_a}}{I_x}

Yaw control moment from rudder:

N_{\delta_r}=\bar{q}SbC_{n_{\delta_r}}\delta_r

Yaw angular acceleration:

\displaystyle \dot{r}\approx\frac{N_{\delta_r}}{I_z}

Use

Control authority depends on dynamic pressure, flow attachment, hinge moment, actuator force, structural stiffness and control-law limits. Low-speed authority and high-speed load limits may govern different parts of the envelope.

Actuator Position and Rate Limits

Position saturation:

|\delta|\le\delta_{max}

Rate saturation:

|\dot{\delta}|\le\dot{\delta}_{max}

Minimum slew time:

\displaystyle t_{slew}=\frac{|\Delta\delta|}{\dot{\delta}_{max}}

Command tracking error during rate limiting:

e_\delta(t)=\delta_{cmd}(t)-\delta(t)

Use

An actuator that is fast enough for trim may still be too slow for gust rejection, upset recovery, flutter excitation limits or envelope protection.

Sampling, Filtering and Delay

Nyquist condition:

f_s>2B

Delay phase lag:

\phi_d=-\omega T_d

In degrees at frequency f:

\phi_d=-360fT_d

Approximate remaining phase margin:

PM_{remaining}=PM_{continuous}+\phi_d

where \phi_d is negative.

Use

Flight-control implementation delay includes sensor sampling, filtering, bus transfer, computation, scheduling, command output and actuator update. A control law that has good continuous-time margin can lose robustness after implementation.

Envelope and Margin Checks

Dynamic-pressure margin:

\displaystyle M_q=\frac{\bar{q}_{limit}-\bar{q}_{test}}{\bar{q}_{limit}}

Load-factor margin:

\displaystyle M_n=\frac{n_{limit}-n_{meas}}{n_{limit}}

Control-position margin:

\displaystyle M_\delta=1-\frac{|\delta_{trim}|+|\delta_{maneuver}|}{\delta_{max}}

Actuator-rate margin:

\displaystyle M_{\dot{\delta}}=1-\frac{|\dot{\delta}_{required}|}{\dot{\delta}_{max}}

Use

Margins should be tied to the relevant failure mode. A load-factor margin does not replace control-position margin, phase margin, sensor validity or flutter margin.

Worked Example 1: Elevator Trim Estimate

At a flight condition, use the simplified moment model:

C_m=C_{m0}+C_{m_\alpha}\alpha+C_{m_{\delta_e}}\delta_e

with:

C_{m0}=0.040
C_{m_\alpha}=-0.85\ \text{rad}^{-1}
C_{m_{\delta_e}}=-1.15\ \text{rad}^{-1}
\alpha=5.0^\circ=0.0873\ \text{rad}

Set C_m=0:

\displaystyle \delta_e=-\frac{0.040+(-0.85)(0.0873)}{-1.15}
\delta_e=-0.0298\ \text{rad}=-1.71^\circ

Engineering comment: the sign depends on the deflection convention. The result is a trim estimate, not a control-surface release. Check hinge loads, actuator range, stick force, tail stall, Mach effects and center-of-gravity range.

Worked Example 2: Static Margin from Neutral Point

The neutral point is:

\displaystyle \frac{x_{np}}{\bar{c}}=0.42

The center of gravity is:

\displaystyle \frac{x_{cg}}{\bar{c}}=0.30

Static margin:

SM=0.42-0.30=0.12=12\%

With:

C_{L_\alpha}=5.5\ \text{rad}^{-1}

the simplified pitching-moment slope is:

C_{m_\alpha}\approx-(5.5)(0.12)=-0.66\ \text{rad}^{-1}

Engineering comment: the aircraft is statically stable by this screen. That does not prove acceptable handling quality; damping, elevator authority, configuration changes, maneuvering loads and control-law mode still need evidence.

Worked Example 3: Banked Turn Load and Stall Margin

An aircraft flies a coordinated level turn at:

\phi=45^\circ

Load factor:

\displaystyle n=\frac{1}{\cos45^\circ}=1.414

If wings-level stall speed is:

V_{S,1}=32\ \text{m/s}

then turn stall speed is:

V_{S,n}=32\sqrt{1.414}=38.1\ \text{m/s}

At:

V=70\ \text{m/s}

turn rate is:

\displaystyle \Omega=\frac{9.81\tan45^\circ}{70}=0.140\ \text{rad/s}=8.03^\circ/\text{s}

Turn radius:

\displaystyle R=\frac{70^2}{9.81\tan45^\circ}=499\ \text{m}

Engineering comment: stall margin looks adequate in this simplified calculation, but the maneuver still requires thrust, roll authority, structural load margin, buffet margin and air-data validity.

Worked Example 4: Damping Ratio from Decay Peaks

A flight-test pulse produces a decaying pitch-rate oscillation. The first peak is:

A_1=4.0\ \text{deg/s}

After three cycles:

A_4=2.2\ \text{deg/s}

Logarithmic decrement:

\displaystyle \delta=\frac{1}{3}\ln\left(\frac{4.0}{2.2}\right)=0.199

Damping ratio:

\displaystyle \zeta=\frac{0.199}{\sqrt{4\pi^2+0.199^2}}=0.0317

If measured period is:

T=1.8\ \text{s}

then:

\displaystyle \omega_d=\frac{2\pi}{1.8}=3.49\ \text{rad/s}
\displaystyle \omega_n=\frac{3.49}{\sqrt{1-0.0317^2}}=3.49\ \text{rad/s}

Engineering comment: damping is low. A single calculation should trigger review of test repeatability, sensor filtering, pilot input, mass properties and whether the mode is close to a structural or control-loop interaction.

Worked Example 5: Actuator Rate and Delay Margin

A control law commands an elevator step:

|\Delta\delta|=12^\circ

The actuator rate limit is:

\dot{\delta}_{max}=40^\circ/\text{s}

Minimum slew time:

\displaystyle t_{slew}=\frac{12}{40}=0.30\ \text{s}

If the intended response assumes the movement occurs within:

0.15\ \text{s}

the actuator is rate limited and the simulation assumption is invalid.

Now check implementation delay. Continuous-time phase margin is:

PM_{continuous}=50^\circ

Crossover frequency is:

f_c=3.0\ \text{Hz}

Total implementation latency is:

T_d=25\ \text{ms}=0.025\ \text{s}

Delay phase lag:

\phi_d=-360(3.0)(0.025)=-27^\circ

Remaining phase margin:

PM_{remaining}=50^\circ-27^\circ=23^\circ

Engineering comment: both actuator rate and implementation delay reduce the real aircraft margin. The controller should not be released from the ideal continuous-time model; hardware timing and actuator dynamics must be included.

Validation Checklist

Before accepting a flight-dynamics or flight-control calculation, confirm:

  1. Axis, sign and unit conventions are explicit.
  2. Aerodynamic derivatives match configuration, Mach number, Reynolds number and control deflection range.
  3. Mass, center of gravity and inertia match the tested or released state.
  4. Trim, control-position and actuator-rate margins are checked across the envelope.
  5. Static and dynamic stability are evaluated at representative flight points.
  6. Sampling, filtering, latency and quantization are included in closed-loop margins.
  7. Sensor validity, estimator drift and degraded modes are defined.
  8. Flexible modes and flutter margins are separated from rigid-body modes.
  9. Failure cases state remaining control authority and transition logic.
  10. Simulation predictions are compared with ground, hardware-in-the-loop or flight-test evidence before envelope expansion.

Common Mistakes

Common mistakes include:

  • mixing degrees and radians in derivatives or control deflections;
  • using a stable trim point as proof of dynamic stability;
  • applying aerodynamic derivatives outside the tested envelope;
  • ignoring actuator rate limits in control-law simulation;
  • adding filters to reduce noise without checking phase margin;
  • using air-data values without consistency checks during sensor faults;
  • validating nominal flight while leaving degraded modes untested;
  • treating rigid-body stability as proof that aeroelastic interaction is safe.

Flight dynamics is a closed-loop aircraft problem. The usable result is not one equation; it is a consistent set of forces, moments, states, sensors, actuators, software timing, structural modes and validation evidence.

REF

See also