draft

(Note: the US spelling is maneuver, but it took me years of muscle memory training to spell manoeuvre correctly, and I won’t be losing it easily.)

Acceleration, Manoeuvres, and Aircraft Loading

The preceding analysis has been constrained to steady flight - that is, with zero acceleration. For the cases of climb, the climb rate was assumed to be steady.

To understand unconstrained aircraft manoeuvres requires an understanding of accelerated flight.

Manoeuvre will be broken down into horizontal (e.g., flat or banked turns) and vertical manoeuvres (e.g., loops, pull-ups), comprising curvilinear motion. Such manoeuvres are the result of a force perpendicular to the flight path, giving a normal acceleration.

All of the manoevures discussed are the result of a variation in lift, which can be large. Consider that the dynamic pressure rises with the square of the forward speed, so a five-fold speed increase results in twenty-five times the aerodynamic forces.

Before discussing manoeuvres, a means to represent the allowable amount of load on an aircraft will be introduced.

Load Factor

The load that can be safely taken through an aircraft dictates the load limits on an aircraft - for this the load factor is introduced as a non-dimensional measure of the load variation.

n=LW(1)n=\frac{L}{W}\tag{1}

There are two structural limits defined for aircraft:

  • Limit Loads, nln_{l} are the loads at which plastic deformation will occur. At flight with 1<n<n11<n<n_1, elastic structural deformation will occur on parts of the aircraft, and the parts will return to the design or equilibrium position once the loads are removed

    If flight occurs at n>nln>n_l, the aircraft will require inspection and likely replacement of parts.

  • Ultimate Loads, nun_{u} are the loads at which failure will occur. At flight with n>nun>n_u, parts of the aircraft will break

The numerical value of the load factor is an exercise in structural analysis, whereby the loads and their paths are applied to a model of the aircraft, and a determination of nln_l and nun_u can be made.

In general, nln_l and nun_u will be defined with a safety factor of 50% - such that plastic structural deformation may not actually occur until nl1.5n_l*1.5. This should not be considered a margin to ‘play with’, however!

Values of the load factor

FAR 23 (Federal Aviation Requirement) dictates the minimum load factor required for three categories of aircraft:

  • Commuter Aircraft
  • Utility Aircraft
  • Aerobatic Aircraft

The minimum load factors are defined as (since this is from FAR 23, WW is defined in lbflbf)

Aircraft CategoryMinimum positive load factor, n+n+Minimum negative load factor, nn-
Normal/Commutersmaller of 2.1+24000W+10002.1+\frac{24000}{W+1000} or 3.83.80.4n+0.4n+
Utility4.44.40.4n+0.4n+
Aerobatic6.06.00.5n+0.5n+

Often pilots will talk about load factor in terms of “g”s - for straight and level flight, 1gg feels like regular earth gravity, whilst an n=2n=2 or 2g2g manoeuvre makes the occupants feel twice as heavy.

Since the load factor represents the amount of lift being produced, it is an easy and instructive means to define the structural limits of an aircraft - e.g.,e.g., the loads for which the wings will break off.

V-n diagram

The values of the limit loads, and ultimate loads are presented on a graph vs. airspeed. For the following, the loads will be presented against equivalent airspeed because that enables a single graph to be plotted, but for real aircraft this is usually presented against indicated airspeed, with several lines representing different altitudes.

There are many different means of showing how to construct a VnV-n diagram in textbook and online, but most seem geared up to helping pilots understand them rather than aerospace engineers. For example - the following show good examples of how to formulate a VnV-n diagram, website one, website two.

For the following example, a representative jet trainer aircraft will be used, with structural limit loads of 3.0<nl<7.0-3.0<n_l<7.0 and ultimate loads of 5.0<nu<11.0-5.0<n_u<11.0.

For a range of flight speeds up to M=1.0M=1.0, the loads above can be plotted against EASEAS taking into account the nn values above only.

Stall limits on V-n diagram

At this stage, the VnV-n diagram isn’t particularly useful. You can hover over the graph and see that values of allowable load is constant with airspeed.

The graph will be adapted by inclusion of other limitations. The first of which is that at certain airspeeds, the aircraft will stall prior to nln_l being reached, and therefore this influences the speed at which manoeuvres can be attempted safely.

It can be shown from the definition of the load factor, and the definition of the lift coefficient that the load factor associated with stall is given by

nstall=CL,max12ρV2SWn_{stall}=\frac{C_{L,max}\,\tfrac{1}{2}\rho V^2\,S}{W}

If the aircraft is taken to have a wing area of 16m2^2, and a weight of 53kN, and a CL,maxC_{L,max} of 1.6 with a CL,minC_{L,min} of -1.0, then the variation of nstalln_{stall} can be overlaid on the previous graph

It can be appreciated that if the stall occurs for n<1n < 1, then steady flight is not possible at this speed.

Manoeuvre Speed

The intersection of the stall boundary and the limit load defines VAV_A, the Manoeuvre Speed. Sometimes this is called the corner speed.

At speeds below VAV_A, fore/aft motion of the stick cannot produce enough load for structural damage to occur as the flow will separate before reaching an incidence at which nln_l would occur. Hence at speeds below VAV_A, the aircraft is stall limited.

Hence VAV_A is the highest speed for safe application of maximum control deflection, whereas at speeds above VAV_A, the controls inputs must be limited to avoid overloading the airframe.

VAV_A can be defined in terms of the maxmimum lift coefficient and the aircraft parameters

VA=2nlWCLmaxρSV_A=\sqrt{\frac{2\cdot n_l\cdot W}{C_{L_{max}}\,\rho\,S}}
Manoeuvre Speed on a real aircraft

On a real aircraft, VAV_A may be defined at a lower speed allowing for a dynamic overshoot.

Furthermore, VAV_A is defined only for pure pitching motion. Combinations of pitch, roll, and yaw can lead to increased loads, as can dynamic motion which allows flow to stay attached for higher angles of attack than the steady state condition.

The latter proved fatal in 2001, when American Airlines Flight 587 took off from JFK. The Airbus A300-605R followed Japan Airlines Flight 47, and encountered wake turbulence.

The First Officer applied repeated opposite rudder inputs (which, as an aside, is a great way to set up a Dutch Roll oscillation), which increased the load on the vertical stabiliser until it ultimately sheared off - which occurred in <7s.

All 260 people aboard the aircraft, and 5 people on the ground were killed in the crash. This is America’s second-deadliest aviation accident.

There were other contributing factors to the accident, such as the light pedal forces on the aircraft misleading the pilots to the tail aerodynamic forces. But the salient point here is that the airframe was destroyed due to aerodynamic load at a velocity far lower than the manoeuvring speed.

The sobering story here is included to show the real world application, and limitations of the theory taught in this class. Manoeuvring speed is a useful tool to understand loads on an airframe, and we will use it to understand the limits in certain manoevures - whether they are stall limited or nn limited.

However, maoeuvring speed is is often poorly understood, and this can be disastrous.

High-Speed Limit

The graph above has been present with all speeds up the speed of sound at sea-level. Not all aircraft can achieve this speed, some due to powerplant limitations, and others due to dangerous aerodynamic phenomena including:

  • Torsional Divergence - where the wings twist up to the point that they snap off. Less common with wing aft-sweep, and impossible for certain planforms (e.g., the Spitfire had an imaginary divergence speed).
  • Flutter - a coupled bend/twist oscillation that is negatively aerodynamically damped which, again, can lead to the wings or other aerodynamic surface snapping off.
  • Control Reversal - where the moment provided by the control surface causes the entire wing to flex to the point that the incremental aerodynamic change due to control surface deflection is negated by the change to flexure of the wing. At speeds above the control reversal speed the controls are…well…reversed.
  • Buffeting - high frequency oscillatory aerodynamic phenomena that may occur at high speeds.

The mechanism of action is not of interest of this course - though if you took MMAE 304 Aerostructures, you will know how to determine divergence speed and control reversal speed.

Fundamentally, there will be a defined speed that the aircraft cannot fly above and this needs to be represented on the VnV-n diagram. This is called the dive speed, VDV_D.

For the jet trainer aircraft, consider VD=300V_D=300m/s. This can be represented on the figure as a vertical line on the right-hand side.

The area inside the graph above represents the envelope, and is the source of the phrase out of the envelope.

Loops

In the following, loops will be analysed using the definition of load factor and the basics of motion in a circle. If you can understand the forces during a loop, you can understand the forces during a general pull-up manoeuvre

If the pilot pulls back on the stick, provided there is sufficient thrust, the angle of attack will be increased and the lift will also increase. This increases the load factor.

Constant Radius Loop

A ‘perfect’ constant radius loop is:

  • Constant airspeed
  • Constant radius
Fig. 2: Constant Radius Loop

If the angular displacement is denoted by γ\gamma, with γ=0\gamma=0 being the bottom of the loop, increasing clockwise, then the equations of motion are, in the aircraft longitudinal direction

TDWsinγ=0T-D-W\sin\gamma=0

and

LWcosγ=WV2grL - W\cos\gamma=\frac{W\,V^2}{g\,r}

So, at different points in the circle the equations of motion give, with Fc=WV2grF_c = \frac{W\,V^2}{g\,r}:

Horizontal EquilibriumVertical Equilibrium
Point ATD=0T-D=0LW=FcL-W=F_c
Point BTDW=0T-D-W=0L=FcL=F_c
Point CTD=0T-D=0L+W=FcL+W=F_c
Point DTD+W=0T-D+W=0L=FcL=F_c

For the perfect loop, the normal acceleration is constant, and the load factor varies according to

n=cosγ+V2grn=\cos\gamma + \frac{V^2}{g\,r}

so at the different points of the loop above, the load factor is:

Load Factor
Point A1+V2gr1+\frac{V^2}{g\,r}
Point BV2gr\frac{V^2}{g\,r}
Point CV2gr1\frac{V^2}{g\,r}-1
Point DV2gr\frac{V^2}{g\,r}

Hence it is the load factor required to initiate the loop that will set the minimum turn radius, and hence the minimimum turn radius is given by

rmin=V2g(n1)r_{min}=\frac{V^2}{g\left(n-1\right)}

Hence the minimum value of the equation above will give the minimum radius for a constant radius loop. For the aircraft VnV-n diagram already constructed, the ratio of this can be plotted

For reasons that can be readily appreciated, the lift and the thrust must be continually varied to maintain constant FCF_C required for a constant radius loop.

For these reasons, your average loop looks something more like this - consider a loop with constant load factor.

Constant Load Factor Loop

The equation for the load factor in a constant radius loop can be rearranged for the radius

r=V2g(ncosγ)r=\frac{V^2}{g\,\left(n-\cos\gamma\right)}

radius, here might be a little confusing since in the above it refers to the radius of a theoretical circle that would be flown if, at any point in the loop, the instantaneous value of the pitching velocity were maintained.

Effectively, this means that the distance flow during a given section of the loop is inversely proportional to the load factor - so the aircraft flies further during the parts where cosγ=0\cos\gamma=0 since rr is larger there. This means that the shape flow is elongated vertically to give a tighter turn at the top and bottom of the loop.

Fig. 2: A loop with constant load factor and constant speed (representative, not to scale)

You should be able to answer qualitative questions about loops, and understand the radius at given points in a loop for different load factors, using the equations given.

The following code is included to help visualise the shape of a constant load factor loop. Run the code yourself - try and change the terms; maximium load factor, entry speed, and see the change to the loop shape.

Steady Turns

A steady turn is one for which there is no tangential acceleration:

  • Circular motion in a horizontal plane
  • Constant turn radius, RR,
  • Constant TAS
Fig. 3: General Steady Turn

The aircraft will have an angular velocity, ω\omega, which can be determined from the flightspeed and turn radius

V=RωV=R\,\omega

The centripetal (centre-seeking) acceleration required to maintain the turn is

ac=Rω2=RV2R2=V2R\begin{align} a_c&=R\,\omega^2\\ &=R\frac{V^2}{R^2}\\ &= \frac{V^2}{R} \end{align}

and hence from Newton’s second law, the associated force reqiured, FcF_c is

Fc=mac=mV2RF_c = m\cdot a_c = \frac{m\,V^2}{R}

This side force can be created in one of two ways:

Flat Turns

In a flat turn, the rudder is utilised to create the side force with wings held level. This is not preferred, because:

  • Very inefficient; small sideforce, and therefore large turn radius
  • Drag increase due to fuselage sideslip; increased fuel burn
  • Dynamic pressure imbalance on wings causes aileron adjustment to be required
  • Perceieved centrifugal force by occupants

Banked Turn

In a banked turn, the aircraft is rolled through angle ϕ\phi so part of the lift provides the sideforce:

Fig. 4: Banked Turn

In a banked turn, the lift is increased by the increment ΔL\Delta L, which can be related to the load factor

L+ΔL=nWL+\Delta L = n\cdot W Lcosϕ=WL\cos\phi = W     cosϕ=n1\implies \cos\phi = n^{-1}

This can be related to required force

Fc=mV2R=Lsinϕ=nWsinϕF_c = \frac{m\, V^2}{R} = L\cdot\sin\phi = n\cdot W\cdot\sin\phi WV2gR=nWsinϕ\frac{W\cdot V^2}{g\cdot R} = n\cdot W\cdot\sin\phi

The quantity V2/RV^2/R can be represented in three useful formats

V2R=gnsinϕ=gtanϕ=gn21\begin{align}\frac{V^2}{R}&=g\cdot n\cdot\sin\phi\\ &= g\cdot\tan\phi\\ &= g\sqrt{n^2-1}\end{align}

or the turn radius can be represented independently of the velocity by substituting LL for WW and using the definition of the lift coefficient

R=2WgρSCLsinϕR=\frac{2\,W}{g\,\rho\,S\,C_L\,\sin\phi}

Hence for a given speed:

  • For constant nn, turn radius is directly proportional to the square of flight speed
  • For constant VV, the turn radius is inversely proportional to the load factor

We can see the conditions that will minimise RR:

  • High density, therefore low altitude
  • Low wing loading (W/SW/S)
  • High lift coefficient
  • High bank angle

Turn Rate

The turn rate is the rate of change of heading, ψ\psi

ω=dψdt=VR\omega = \frac{\text{d}\psi}{\text{d}t} = \frac{V}{R} ω=gn21V=gρVCLsinψ2WS\omega = \frac{g\sqrt{n^2-1}}{V} = \frac{g\,\rho\,V\,C_L\,\sin\psi}{2\frac{W}{S}}