Pitch, roll, yaw, stability margins, dynamic modes, performance equations, and trajectory analysis. The complete guide to how aircraft and rockets are controlled and how their motion is predicted.
Every aircraft motion (from a gentle bank to a high-g break turn) is described as a combination of rotations about three axes. Understanding these axes and the control surfaces that command them is the first step in flight mechanics.
Rotation about the lateral Y-axis, running through both wingtips. Pitching up increases AoA and generates more lift, but also more induced drag. Exceed the critical AoA and the wing stalls. The horizontal stabiliser provides pitch stability; the elevator commands pitch changes. The distance between the aircraft's centre of gravity and the neutral point (static margin) determines how much stabilising moment the tail must generate.
Rotation about the longitudinal X-axis, running nose to tail. Ailerons deflect differentially (one down (increased lift), one up (decreased lift)) creating asymmetric lift that rolls the aircraft. Dihedral (wings angled upward from root to tip) provides positive roll stability: a disturbance into a bank creates a restoring moment. High-speed roll also generates adverse yaw from differential induced drag, requiring coordinated rudder input.
Rotation about the normal Z-axis, running top to bottom through the aircraft centre. Controlled by the rudder on the vertical fin. The fin acts as a weathervane providing directional stability (tendency to align with the relative wind). Ailerons cause adverse yaw. The rising wing generates more induced drag, pulling the nose toward the outside of the turn. Coordinated rudder input eliminates this, keeping the slip indicator centred.
Roll generates adverse yaw. Yaw can excite Dutch roll or spiral divergence. All three axes interact, flight control computers decouple them for the pilot.
Stability determines whether an aircraft naturally returns to equilibrium after a disturbance. It is the fundamental property that separates a controllable aircraft from an uncontrollable one.
An aircraft is statically stable if, following a disturbance, the initial aerodynamic forces and moments tend to restore it to its original attitude. For pitch: a nose-up gust increases AoA, which must generate a nose-down pitching moment to be stable.
The neutral point (NP) is where an increase in AoA produces zero net pitching moment change. The static margin is the distance CG to NP as a percentage of mean aerodynamic chord (MAC). CG must be forward of NP for stability. Typical airliners: 5 to 15% MAC static margin. Modern fighters may fly at zero or negative margin for manoeuvrability, requiring active FBW stabilisation.
A statically stable aircraft can still be dynamically unstable. The oscillation may grow over time rather than damp. The four key dynamic modes are:
Phugoid, long-period exchange of kinetic and potential energy at roughly constant angle of attack. Lanchester's approximation gives T ≈ π√2·V/g, so the period is set by airspeed: about 60 to 120 s for a transport, shorter on a light aircraft. Damping is ζ ≈ 1/(√2·L/D), so a cleaner aircraft has a less damped phugoid, not a longer one.
Short-period, fast (1 to 5 s) pitch oscillation at near-constant speed. Must be well damped; MIL-STD-1797 Level 1 wants ζ between about 0.35 and 1.3. Dangerous if underdamped.
Dutch roll, coupled yaw-roll oscillation, aggravated by sweep and by high dihedral effect relative to directional stability. Suppressed by the yaw damper.
Spiral mode, slow non-oscillatory bank divergence. Time constant typically 200 to 400 s on a transport, so the pilot corrects it without noticing. A mildly divergent spiral is normal and acceptable; it is traded against Dutch roll damping.
Performance analysis answers the operational questions: how far can it fly, how high, how fast, and how quickly does it climb? These quantities are derived directly from the aerodynamic and propulsion parameters.
The Breguet range equation describes how far an aircraft can fly on a given fuel load as a function of aerodynamic efficiency and propulsion efficiency. Range improves with aerodynamic efficiency and with low specific fuel consumption, but the exact condition depends on the powerplant. A propeller aircraft maximises range at maximum L/D. A jet maximises V·(L/D), which for a parabolic polar occurs at CD0 = CDi/3, about 1.32× the minimum-drag speed. Saying "jets cruise at max L/D" is the single most common error here: they deliberately cruise faster than that, trading a little aerodynamic efficiency for a lot of distance per hour. Either way the speed scales as √(W/S), so a heavier aircraft flies faster for the same coefficient.
Endurance (maximum time aloft) peaks at a different, lower speed than range: minimum power required for a prop, minimum drag for a jet.
Holds L/D and SFC constant through a cruise-climb at fixed CL. Real cruise goes in steps, so treat the answer as a ceiling and subtract climb, descent, reserves and routing. Watch the units too: with SFC in kg/(N·s) the fuel weight flow is g₀·SFC·T, so the equation is R = V·(L/D)·ln(W₀/W₁)/(g₀·SFC). Leave g₀ out and your range comes out 9.81 times too big.
Rate of climb RC = (T − D) · V / W, excess thrust times velocity divided by weight. Maximum RC occurs at the speed where excess power is greatest. As altitude increases, air density drops and engine thrust decreases, reducing excess power. The absolute ceiling is the altitude where rate of climb reaches zero. The service ceiling is where RC = 100 ft/min (the operational limit).
Turn performance is governed by the load factor n = L/W, which is dimensionless: a 60° banked level turn requires n = 2, not "2g". In level flight n = 1/cos φ. Stall speed in a turn rises as √n, so a 60° bank raises it by 41% over the straight-and-level value. Turn radius is R = V²/(g√(n²−1)), which is why a tight turn needs speed and load factor, and why both are bounded by the V-n diagram.
Trajectory analysis integrates the equations of motion to predict where an aircraft or rocket will be at any given time. For rockets, this is the critical calculation that determines whether the payload reaches orbit.
A rigid aircraft has 6 degrees of freedom, three translational (surge, heave, sway) and three rotational (pitch, roll, yaw). The full 6-DOF equations of motion are a coupled set of 12 first-order ODEs. They can be linearised around a trim point to produce separate longitudinal and lateral-directional state-space models:
Longitudinal: states are u (forward speed), w (heave velocity), q (pitch rate), θ (pitch angle). Input: elevator δe.
Lateral-directional: states are v (side velocity), p (roll rate), r (yaw rate), φ (bank angle), β (sideslip). Inputs: aileron δa, rudder δr.
A rocket trajectory integrates thrust, gravity, and aerodynamic drag simultaneously. The equations of motion are:
V̇ = (T·cosα − D)/m − g·sinγ
γ̇ = (T·sinα + L)/(m·V) − (g/V)·cosγ
ẋ = V·cosγ · ḣ = V·sinγ · ṁ = −T/(Isp·g₀)
where T is thrust, D is drag, L is any aerodynamic normal force, m is instantaneous mass, α is the thrust angle relative to the velocity vector, and γ is the flight path angle. Writing the equations in the velocity and flight-path frame rather than in x and y keeps drag correctly aligned with the velocity vector, which is where it physically acts; resolving drag into separate x and y components with the flight-path angle is a common way to double-count it. MATLAB's ode45 is the standard tool for integrating these equations numerically. The Tsiolkovsky equation gives the final ΔV; the trajectory integration gives altitude, downrange distance, and dynamic pressure as functions of time.
One stage, constant exhaust velocity, thrust along the velocity vector, and nothing else acting: no gravity, no drag, no steering. That missing 1.5 to 1.7 km/s is why a launch needs 9.3 to 9.5 km/s to reach an orbit moving at 7.8, and why you integrate the trajectory above instead. The g₀ is a unit conversion, 9.80665 m/s², not local gravity, so altitude does not touch it.
The aerodynamic forces (lift, drag, pitching moments) that flight mechanics must balance. Start here first.
Thrust curves, specific impulse, and how engine choice affects trajectory performance and climb rate.
The Tsiolkovsky equation, 6-DOF rocket trajectory, and how orbital insertion is achieved from the same equations.
Simulate the aerodynamic forces that appear in your flight mechanics equations: Cp distribution, CL, CD.
State-space models, eigenvalue analysis, LQR control design, and step response simulation: all in working MATLAB code.
Apply trajectory equations live: calculate delta-v, check stability margins, and simulate apogee in the rocket builder.
SheCodes Lab teaches Python and C++ from scratch, side by side, free, no experience needed. Includes an engineering module covering NumPy, pandas, ISA models, cost index, and flight data analysis. The same tools used to build the calculators on this site.
shecodeslab.com →