Aerodynamic Moments (i.e., Torques) Robert Stengel, Aircraft Flight Dynamics, MAE 331, 2018
Learning Objectives
• Aerodynamic balance and moment • Aerodynamic center, center of pressure, neutral point, and static margin • Configuration and angle-of-attack effects on pitching moment and stability • Configuration and sideslip-angle effects on lateral-directional (i.e., rolling and yawing) aerodynamic moments Reading: • Tail design effects on airplane Flight Dynamics Aerodynamic Coefficients, 96-118 aerodynamics
Copyright 2018 by Robert Stengel. All rights reserved. For educational use only. http://www.princeton.edu/~stengel/MAE331.html 1 http://www.princeton.edu/~stengel/FlightDynamics.html
Review Questions § Why is induced drag proportional to angle of attack squared? § What spanwise lift distribution gives minimum induced drag? § Why can lift and drag coefficients be approximated by the Newtonian-flow assumption at very high angle of attack? § How does profile drag vary with Mach number? § What are some functions of secondary wing structures? § What is the primary function of leading edge extensions? § What is the “Area Rule”? 2
1 Handbook Approach to Aerodynamic Estimation • Build estimates from component effects – Technical reports, textbooks, ... – USAF Stability and Control DATCOM (download at http://www.pdas.com/datcomb.html) – USAF Digital DATCOM (see Wikipedia page) – ESDU Data Sheets (see Wikipedia page)
Interference Effects Tail Interference Aerodynamics Effects Fuselage Aerodynamics Wing Aerodynamics Interference Effects
3
Moments of the Airplane
4
2 Airplane Balance • Conventional aft-tail configuration – c.m. near wing's aerodynamic center [point at which wing's pitching moment coefficient is invariant with angle of attack ~25% mean aerodynamic chord (mac)] • Tailless airplane: c.m. ahead of the neutral point
Northrop N-9M
Douglas DC-3
5
Airplane Balance • Canard configuration: – Neutral point moved forward by canard surfaces – Center of mass may be behind the neutral point, requiring closed-loop stabilization • Fly-by-wire feedback control can expand envelope of allowable center-of-mass locations
Grumman X-29
McDonnell-Douglas X-36
6
3 Moment Produced By Force on a Particle
r! : Cross-product-equivalent matrix
Cross Product of Vectors
i j k x y z r × f = = (yfz − zfy )i + (zfx − xfz ) j + (xfy − yfx )k
fx fy fz
⎡ ⎤ ⎡ ⎤ yf − zf ⎡ ⎤ mx ⎢ ( z y ) ⎥ ⎡ 0 −z y ⎤ fx ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ m = my = zf − xf = r × f ! r"f = z 0 −x fy ⎢ ⎥ ⎢ ( x z ) ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ mz ⎢ ⎥ −y x 0 fz ⎣⎢ ⎦⎥ (xfy − yfx ) ⎣ ⎦ ⎣⎢ ⎦⎥ ⎣⎢ ⎦⎥ 7
Forces and Moments Acting on Entire Airplane Force Vector
! X $ # B & fB = # YB & # Z & " B %
Moment Vector
! L $ # B & mB = # M B & # N & " B % 8
4 Aerodynamic Force and Moment Vectors of the Airplane Force Vector ! f $ ! X $ # x & # B & f fB = ∫ # y &dx dydz = # YB & Surface# & # & f Z "# z %& " B %
Moment Vector " % " % (yfz − zfy ) L $ ' $ B ' $ ' mB = ∫ (zfx − xfz ) dx dydz =$ M B ' Surface$ ' $ ' $ ' N xfy − yfx # B & # ( ) & 9
Pitching Moment of the Airplane
10
5 Body-Axis Reference Frames • Reference frame origin is arbitrary; it is a fiducial point – x axis along centerline – Tip of nose: All values of x on airframe are negative, but nose shape could change – Forward-most bulkhead: Fixed for all manufacturing measurements – Center of mass: Rotational center, but changes with fuel use, payload, etc.
11
F-86 Nose Variations
F-86A F-86D
F-86E F-86D
12
6 Pitching Moment (moment about the y axis) Pressure and shear stress differentials x moment arms Integrate over the airplane surface to produce a net pitching moment
Body - Axis Pitching Moment = M B = − Δp x, y + Δs x, y x − x dx dy ∫∫ ⎣⎡ z ( ) z ( )⎦⎤( cm ) surface + Δp y,z + Δs y,z Δp z − z dydz ∫∫ ⎣⎡ x ( ) x ( )⎦⎤ x ( cm ) surface
13
Pitching Moment (moment about the y axis) Distributed effects can be aggregated to local centers of pressure indexed by i
I I M B ≈ −∑ Zi (xi − x cm ) + ∑ Xi (zi − zcm ) i=1 i=1 +Interference Effects + Pure Couples Net effect expressed as
M B = Cmq Sc 14
7 Pure Couple • Net force = 0 • Net moment ≠ 0
Rockets Cambered Lifting Surface
• Cross-sectional area, S(x) • x positive to the right Fuselage • At small α – Positive lift slope with dS(x)/dx > 0 – Negative lift slope with dS(x)/dx < 0 • Fuselage typically produces a destabilizing (positive) pitching
moment [ “Apparent mass” effect] 15
Lift Coefficient of a Cone
Munk’s airship theory (potential flow)
⎛ Length ⎞ Sbase C ! ±2 fcn α Lcone ⎝⎜ Diameter ⎠⎟ S ⎧ 0.4, Len Dia = 2 ⎫ ⎪ ⎪ S ! ±2 ⎨ 0.84, Len Dia = 5 ⎬ base α S ⎪ 0.94, Len Dia =10 ⎪ ⎩ ⎭ Sbase: cross-sectional area where flow separates
Munk, NACA-TR-184, 1924 16
8 Net Center of Pressure Local centers of pressure can be aggregated at a net center of pressure (or neutral point) along the body x axis
⎡ x C + x C + x C + ...⎤ ( cp N )wing ( cp N ) fuselage ( cp N )tail x = ⎣ ⎦ cpnet C Ntotal C = −C S = reference area N Z CA = −CX
Body Axes
Wind Axes (w.r.t, velocity vector)
17
Static Margin
• Static margin (SM) reflects the distance between the center of mass (cm) and the net center of pressure (cp) • Body axes • Normalized by mean aerodynamic chord • Does not reflect z position of center of pressure • Positive SM if cp is behind cm
100 x − x ( cm cpnet ) Static Margin ! SM = B , % c ≡ 100 h − h % ( cm cpnet )
xcm hcm ! 18 c
9 Static Margin
Static Margin = SM
100 xcm − xcp = ( net ), % c ≡ 100 h − h % ( cm cpnet )
19
Effect of Static Margin on Pitching Coefficient • Zero crossing determines trim angle of attack, i.e., sum of moments = 0 • Negative slope required for static stability
• Slope, ∂Cm/∂α, varies with static margin
C mo M = C + C α q Sc αTrim = − B ( mo mα ) C mα 20
10 Pitch-Moment Coefficient Sensitivity to Angle of Attack
For small angle of attack and no control deflection
Cm ≈ −CN hcm − hcp ≈ −CL hcm − hcp α αnet ( net ) αnet ( net ) x − x x x ⎛ cm cpwing ⎞ ⎛ cm − cpht ⎞ ≈ −CL − CL αwing αht ⎜ ⎟ ⎝⎜ c ⎠⎟ ⎝ c ⎠
referenced to wing area, S
21
Horizontal Tail Lift Sensitivity to Angle of Attack
2 ⎡ ⎤ ⎛ ∂ε ⎞ ⎛ Sht ⎞ ⎛ Vht ⎞ CL horizontal = CL 1− ηelas ⎢( αht ) ⎥ ( αht )ref S ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ tail = ht ⎝ ∂α ⎠ ⎝ S ⎠ V ⎣ ⎦ref = S ⎝ N ⎠
Vht : Airspeed at horizontal tail ε : Downwash angle due to wing at tail ∂ε ∂α : Sensitivity of downwash to angle of attack
ηelas : Aeroelastic effect
• Downwash effect on aft horizontal tail • Upwash effect on a canard (i.e., forward) surface 22
11 Aerodynamic Center and Center of Pressure of a Wing ∂C x = x for which m ≡ 0 ac ∂α
= xcp for a symmetric airfoil
≠ xcp for an asymmetric airfoil
NACA 0012
NACA 2412 Airfoil Tools http://airfoiltools.com NACA 4412
23
Effect of Static Margin on Pitching Moment For small angle of attack and no control deflection
M = C q Sc ≈ ⎡C − C h − h α ⎤q Sc B m ⎣ mo Nα ( cm cpnet ) ⎦
≈ ⎡C − C h − h α ⎤q Sc ⎣ mo Lα ( cm cpnet ) ⎦
24
12 Effect of Static Margin on Pitching Moment Sum of moments is zero in trimmed condition
M = C + C α q Sc B ( mo mα )
= 0 in trimmed (equilibrium) flight
Typically, static margin is positive and ∂Cm/∂α is negative for static pitch stability
25
Effect of Elevator Deflection on Pitching Coefficient Control deflection shifts curve up and down, affecting trim angle of attack
M = C + C α + C δE q Sc B ( mo mα mδE )
1 α = − C + C δ E Trim C ( mo mδ E ) mα Elevator deflection effectively changes C mo 26
13 Historical Factoids Aviation in The Great War • 1914-18: World War I changes the complexion of flying SPAD S.VII – Reconnaissance – Air superiority (dog fights) – Bombing – Personal transport • Wrights US monopoly broken by licensing for war effort • Aircraft Design – Biplanes, a few mono- and triplanes – Design for practical functions – Multiple engines, larger aircraft – Aft tails – Increased maneuverability, speed, g-loads, altitude – Improved piston engines – Tractor propellers 27
Maneuvering World War I Aircraft
• Maneuverable aircraft with idiosyncrasies – Rotary engine – Small tail surfaces DeHavilland DH-2 – Reliability issues • Maneuvering to stalls and spins • Snap roll: rudder and elevator • Barrel roll : aileron • Cross-control (e.g., right rudder, left stick) Fokker E.III – glide path control during landing – good view of landing point • Unintended snap rolls led to spins and accidents during takeoff or landing Sopwith Triplane
http://www..com/watch?v=6ETc1mNNQg8youtube
http://www.youtube.com/watch?v=OBH_Mb0Kj2s 28
14 Stability OR Control?
Stability AND Control • Need for better understanding of Flying (or Handling) Qualities – Stability and controllability characteristics as perceived by the pilot • Desired attributes: Stability of the S.E.-5 and controllability of the D.VII 29
Lateral-Directional Effects of Sideslip Angle
30
15 Rolling and Yawing Moments of the Airplane Distributed effects can be aggregated to local centers of pressure indexed by i
Rolling Moment I I LB ≈ ∑Zi (yi − y cm ) − ∑Yi (zi − zcm ) i=1 i=1 +Interference Effects + Pure Couples Yawing Moment
I I N B ≈ ∑Yi (xi − x cm ) − ∑ Xi (yi − ycm ) i=1 i=1
+Interference Effects + Pure Couples 31
Sideslip Angle Produces Side Force, Yawing Moment, and Rolling Moment
§ Sideslip usually a small angle ( 5 deg) § Side force generally not a significant effect § Yawing and rolling moments are principal effects
32
16 Side Force due to Sideslip Angle
∂C Y ≈ Y qS •β = C qS •β ∂β Yβ • Fuselage, vertical tail, nacelles, and wing are main contributors
CY ≈ CY + CY + CY + CY β ( β )Fuselage ( β )Vertical Tail ( β )Nacelles ( β )Wing
S = reference area
33
Side Force due to Sideslip Angle
⎛ ∂CY ⎞ ⎛ Svt ⎞ CY ≈ ηvt ⎜ ⎟ ( β )Vertical Tail ⎝⎜ ∂β ⎠⎟ ⎝ S ⎠ ref = Svt 2 SBase πdBase CY ≈ −2 ; SB = ( β )Fuselage S 4
2 CY ≈ −CD − kΓ ( β )Wing Parasite, Wing
ηvt = Vertical tail efficiency (p. 96, Flight Dynamics) π AR k = 1+ 1+ AR2
Γ = Wing dihedral angle, rad 34
17 Yawing Moment due to Sideslip Angle
∂C ⎛ ρV 2 ⎞ N ≈ n Sb • β = C qSb • β ∂β ⎝⎜ 2 ⎠⎟ nβ
§ Side force contributions times respective moment arms – Non-dimensional stability derivative
Cn ≈ Cn + Cn + Cn + Cn β ( β )Vertical Tail ( β )Fuselage ( β )Wing ( β )Propeller S = reference area 35
Yawing Moment due to Sideslip Angle Vertical tail contribution
Svtlvt Cn ≈ −CY ηvt −CY ηvt ( β )Vertical Tail βvt βvt VVT Sb
lvt Vertical tail length (+) = distance from center of mass to tail center of pressure = x − x [x is positive forward; both are negative numbers] cm cpvt %V 2 ( η =η 1+∂σ ' vt * vt elas ( ∂β) 2 &VN )
Svtlvt VVT = = Vertical Tail Volume Ratio Sb 36
18 Yawing Moment due to Sideslip Angle Fuselage contribution
−2K VolumeFuselage Cn = ( β )Fuselage Sb
1.3 " d % K = $1− max ' # Lengthfuselage & Wing (differential lift and induced drag) contribution
k1 2 Cn = 0.075CL Γ + fcn(Λ, AR,λ)CL ( β )Wing N N k C k C 2 (eq. 2.4-66, Flight Dynamics) ! 0 LN Γ + 1 LN AR
Seckel, from NACA TR-1098, 1950 37
Rolling Moment due to Sideslip Angle Dihedral effect L ≈ C qSb•β lβ
Cl ≈ Cl + Cl + Cl β ( β )Wing ( β )Wing− Fuselage ( β )Vertical Tail
Unequal lift on left and right wings induces rolling motion
38
19 Rolling Moment due to Sideslip Angle
• Wing vertical location effect: Crossflow produces up- and down-wash • Rolling effect depends on vertical location of the wing
• Vertical tail effect
39
Tail Design Effects
40
20 Tail Design Effects • Aerodynamics analogous to those of the wing • Longitudinal stability – Horizontal stabilizer – Short period natural frequency and damping • Directional stability – Vertical stabilizer (fin) • Ventral fins • Strakes • Leading-edge extensions • Multiple surfaces • Butterfly (V) tail – Dutch roll natural frequency and damping • Stall or spin prevention/ recovery • Avoid rudder lock (TBD) 41
Horizontal Tail Size and Location
Grumman XF5F North American F-86
§ 15-30% of wing area § ~ wing semi-span behind the c.m. § Must trim neutrally stable airplane at maximum lift in ground effect § Effect on short period mode § Horizontal Tail Volume: Typical value = 0.48
Sht lht VHT = S c 42
21 Tail Moment Sensitivity to Angle of Attack
2 ⎛ Vht ⎞ ⎛ ∂ε ⎞ ⎡⎛ Sht ⎞ ⎛ lht ⎞ ⎤ Cm = − CL 1− ηelas αht ( αht )ht ⎜ ⎟ ⎜ ⎟ ⎢⎜ ⎟ ⎜ ⎟ ⎥ ⎝ VN ⎠ ⎝ ∂α ⎠ ⎣⎝ S ⎠ ⎝ c ⎠ ⎦ 2 ⎛ Vht ⎞ ⎛ ∂ε ⎞
= − CL 1− ηelas ( αht )ht ⎜ ⎟ ⎜ ⎟ VHT ⎝ VN ⎠ ⎝ ∂α ⎠
S l ht ht Horizontal Tail Volume Ratio VHT = = Sc
43
Pitching Moment due to Elevator Deflection Normal force coefficient variation due to elevator deflection
∂CL Sht CL ! = τ htηht CL ≈ CN δ E ∂δ E ( α )ht S δ E τ ht = Carryover effect
ηht = Tail efficiency factor
ΔC = C δ E C = Horizontal tail lift-coefficient slope L N Nδ E ( α )ht
Sht = Horizontal tail reference area Pitching moment coefficient variation due to elevator deflection l ⎛ S l ⎞ C = C ht ≈ −τ η C ht ht m E N E ht ht L ⎜ ⎟ δ δ c ( α )ht ⎝ S c ⎠
= −τ htηht CL ( α )ht VHT 44
22 Downwash and Elasticity Also Effect Elevator Sensitivity
2 ⎡⎛ ∂CL ⎞ ⎤ ⎛ Vtail ⎞ ⎛ ∂ε ⎞ ⎛ Sht ⎞ = C = C 1− η ⎢⎜ ⎟ ⎥ ( Lδ E ) ( Lδ E ) ⎜ ⎟ ⎜ ⎟ elas ⎜ ⎟ ⎝ ⎠ ref = S ref = Sht ⎝ ⎠ ⎝ ⎠ ⎣ ∂δ E ht ⎦ref = S ⎝ VN ⎠ ∂α S
45
Vertical Tail Location and Size North American P-51B
• Analogous to horizontal tail volume • Effect on Dutch roll mode • Powerful rudder for spin recovery – Full-length rudder located behind the elevator – High horizontal tail so as not to block the flow over the rudder • Vertical Tail Volume: Typical value = 0.18 Curtiss SB2C
Svt lvt VVT = S b
Otto Koppen: “If they build more than one of these, they’re crazy!” http://en.wikipedia.org/wiki/Otto_C._Koppen
46
23 Lateral-Directional Control Surfaces
Elevons
Rudder
47
Yawing Moment due to Rudder Deflection Side force coefficient variation due to rudder deflection ⎛ ∂C ⎞ S C Y = ⎡ C ⎤ τ η vt Y R ! ⎜ ⎟ L vt vt ( δ )ref =S ⎣( α )vt ⎦ref =S ⎝ ∂δ R⎠ ref =S vt S ΔC = C δ R Y Yδ R
Yawing moment coefficient variation due to rudder deflection
l S l vt ⎡ ⎤ ⎛ vt vt ⎞ Cn = − CY ≈ − CL τ vtηvt ( δ R )ref =S ( δ R )ref =S ( α )vt ⎜ ⎟ b ⎣ ⎦ref =Svt ⎝ S b ⎠
⎡ ⎤ = −τ vtηvt CL ( α )vt VVT ⎣ ⎦ref =Svt
48
24 Rolling Moment due to Aileron Deflection
L ≈ Cl qSb •δ A Lδ A For a trapezoidal planform, subsonic flow
2 3 ⎛ C ⎞ CL Lδ ( a )3D ⎡1− k 1− k ⎤ Cl ! − (1− λ) ( δ A )3D ⎜ C ⎟ 1+ λ ⎢ 3 3 ⎥ ⎝ La ⎠ 2 D ⎣ ⎦ y k " , y = Inner edge of aileron, λ = Taper ratio b 2
Ayres Thrush Crop Duster 49
Next Time: Aircraft Performance
Reading: Flight Dynamics Aerodynamic Coefficients, 118-130
Learning Objectives Definitions of airspeed Performance parameters Steady cruising flight conditions Breguet range equations Optimize cruising flight for minimum thrust and power Flight envelope
50
25 Supplemental Material
51
Sopwith Camel
• Rotary engine induced gyroscopic coupling • Highly maneuverable • Aft fuel tank; when full, center of mass was too far aft for stability • Vertical tail too small, spin recovery not automatic with centering of controls
http://www.youtube.com/watch?v=3ApowyEXSXM 52
26 S.E.-5 vs. Fokker D.VII
• RAF S.E.-5: theoretical approach to design – Best WWI design from the Royal Aircraft Factory – Stationary engine – High dihedral – Stable spiral mode – High control forces – Poor maneuverability – Relatively safe and effective
• Fokker D.VII: empirical approach to design – Horn balances to reduce control forces – Stationary engine – Neutral-to-negative stability – Good maneuverability – Relatively dangerous
53
Planform Effect on Center of Pressure Variation with Mach Number
• Straight Wing – Subsonic center of pressure (c.p.) at ~1/4 mean aerodynamic chord (m.a.c.) – Transonic-supersonic c.p. at ~1/2 m.a.c. • Delta Wing – Subsonic-supersonic c.p. at ~2/3 m.a.c.
54
27 Subsonic Pitching Coefficient vs. Angle of Attack (0 < α < 90 )
55