DOCSLIB.ORG

05 Longitudinal Stability Derivatives

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
05 Longitudinal Stability Derivatives Flight Dynamics and Control Lecture 5: Longitudinal stability Derivatives G. Dimitriadis University of Liege 1 Previously on AERO0003-1 • We developed linearized equations of motion Longitudinal direction " % "m 0 0 0 0% " v % −Yv −(Yp + mWe ) −(Yr − mUe ) −mgcosθe −mgsinθe " v % " Yξ Yς % $ ' $ ' $ ' $ ' $ ' 0 I x −I xz 0 0 p $ −L −L −L 0 0 ' p Lξ Lς $ ' $ ' v p r $ ' $ ' "ξ% $ ' $ ' $ ' $ ' $ ' 0 −I xz Iz 0 0 r + −N −N −N 0 0 r = Nξ Nς $ ' $ ' $ ' $ v p r ' $ ' $ ' #ς & $ 0 0 0 1 0' $ϕ ' $ 0 −1 0 0 0 ' $ϕ ' $ 0 0 ' $ ' #$ 0 0 0 0 1&' #$ψ &' − #$ψ &' $ 0 0 ' # 0 Lateral0 direction1 0 0 & # & " % " − − − − θ % m −Xw 0 0 "u % Xu Xw (Xq mWe ) mgcos e "u % " X X % $ ' $ ' η τ $ ' $ ' $ ' $ 0 m − Z 0 0' w $ − − − + θ ' w Zη Zτ "η% ( w ) $ ' + Zu Zw (Zq mUe ) mgsin e $ ' = $ ' $ ' $ ' $ ' $ ' − $ q ' $ q ' Mη Mτ #τ & $ 0 M w I y 0' $−M u −M w −M q 0 ' $ ' $ ' $ ' $ 0 0 0 1' #θ & $ − ' #θ & # 0 0 & 2 # & # 0 0 1 0 & Longitudinal stability derivatives • It has already been stated that the best way to obtain the values of the stability derivatives is to measure them. • However, it is still useful to discuss simplified methods of estimating these coefficients. • Such estimates can be used, for example, in the preliminary design of aircraft. • This lecture will treat longitudinal stability derivatives. 3 Simple example • We keep the quasi-steady aerodynamic assumption. • Assume that the lift of an aircraft lies entirely in the z direction: 1 2 Z = ρU SCL 2 • where CL is the lift coefficient, assumed to be constant in this simple example. • Now consider that there is a perturbation in the airspeed such that U=Ue+u, Ue being the equilibrium airspeed 4 Simple example (2) • The total aerodynamic force in the z direction becomes: 1 2 1 2 2 Z = ρ (Ue + u) SCL = ρ (U e + 2Ueu + u )SCL 2 2 • We linearize such that u2=0: 1 2 Z = ρ (U e + 2Ueu)SCL 2 • Now we can calculate the derivative of Z with respect to u: ! ∂Z Zu = = ρUeSCL 5 ∂u Lift and Drag • This simple example considered only one aerodynamic force and one perturbation. Now we need to consider all aerodynamic loads and all perturbations. • Lift and drag are defined in directions perpendicular and parallel to the direction of motion: U Lift e αe θe V0 γe Xa W Drag e mg Za Pitching Moment, Ma 6 Longitudinal geometry and loads S l = T T L VT T Lw Sc ac cg MT M0 l T mg hc h denotes the cg c position h0 denotes the h c From Aircraft Design Lectures 0 aerodynamic centre 7 Za, Xa and Ma • Clearly, the aerodynamic forces Za and Xa and the aerodynamic moment Ma can be written as: Za = −Lift × cosα − Drag × sinα X a = Lift × sinα − Drag × cosα Ma = M0 + Lift(h − h0)c − LiftT lT • The last equation was obtained during the Aircraft Design lectures • We have assumed that the wing lift is 8 approximately equal to the total lift. In coefficient form • The lift and drag can be written, as usual, in terms of the lift and drag coefficients, so that 1 2 Za = − ρQ∞S(CL cosα +CD sinα) 2 1 2 Xa = ρQ∞S(CL sinα −CD cosα) 2 = 1 ρ 2 + − − M a Q∞Sc (Cm0 CL (h h0 ) VTCLT ) 2 • Using the notation of Aircraft Design, C = C + aα L L0 a=lift curve slope 2 CL C = C + e=Oswald factor 9 D D0 eπAR Angle of attack • Consider an aircraft in equilibrium at angle of attack αe, perturbed by small airspeeds u, w and small pitch rate q. Ue Ue+u −1 W + w −1 W −1 w α = tan e ≈ tan e + tan αe α Ue + u Ue Ue or V0 We w α = αe + Ue We+w 10 Effective camber • A pitching wing can be seen to feature an additional camber distribution due to the pitching motion. Local flow velocity due to pitch rate Additional contribution to C.G. angle of attack q α = − = α eff (h1 h) h0 c /4 Ue h V0 h1 = 3c /4 11 Wing-tail flow geometry • For a static aircraft, the downwash effect of the wing deflects the free stream flow seen by the tailplane by an angle α. • Total angle of attack of tail: αT=α-ε+ηT Tailplane Wing ηT α αT α V0 η ε Elevator V0 βη Trim tab • The total lift on the tailplane is given by: = + α + η + β 12 CLT a0 a1 T a2 a3 η Unsteady wing downwash • Now, consider a dynamic aircraft. A change in downwash at the wing because of a change in pitch angle at time t will at the tail Δt seconds later. • The time difference Δt is the time it takes for air leaving the wing to reach the tail and can be approximated as Δt=lt/Ue. • Then: ∂ε ∂ε ε (t) = α − α!Δt ∂α ∂α • So that: ' ∂ε ) ∂ε l C = aη + a 1− α + a t q + a η LT 1 T 1( * 1 2 ∂α ∂α Ue 13 Lift, drag and moment • In terms of perturbations around a static equilibrium condition, the lift, drag and moment can be written as: = + a + − CL CLe (w (h1 h)q) Ue 2aC = + Le + − CD CDe (w (h1 h)q) Ueeπ AR " " ∂ε %% ∂ε = + − − a1 − a + − − lt − η Cm Cme $h h0 VT $1 '' (w (h1 h)q) VT a1 q VT a2 # a # ∂α &&Ue ∂α Ue Second order perturbation terms in the drag expression have been ignored 14 Added mass terms • The acceleration stability derivatives, e.g. M w , X w , Z w result from the added mass effect of the aircrafts motion. • This motion carries along a certain mass of air around the aircraft. • As a consequence, the aircraft has more inertia than it would have had if it was moving in a vacuum • This is known as the added mass or apparent mass effect 15 Added mass terms (2) • Added mass can be seen as the inertia of cylinder of diameter c and length b that is moving with the wing. • Some added mass terms: c Lift = ρπ S(w + U q) added 4 e c % c ( c % 3c ( π Moment = ρπ S' h − * w − ρπ S' − h* U q − ρU c 2Sq added 4 & 2) 4 & 4 ) e 16 e 16 Perturbation • Consider that, as the flight is perturbed, the total airspeed is: 2 2 2 2 Q∞ = (Ue + u) +(We + w) = V0 + 2Ueu + 2Wew • The cosine and sine terms of the angle of attack become: # w & w cosα = cos% αe + ( = cosα e − sinα e $ Ue ' Ue # w & w sinα = sin% α e + ( = sinαe + cosαe $ Ue ' Ue 17 Force Za • The total force Za becomes: # & 1 2 1 2 a 2aCL = − ρ α + α − ρ % + − α + e + − α( Za Q∞S(CLe cos CDe sin ) Q∞S (w (h1 h)q)cos (w (h1 h)q)sin 2 2 $Ue Ueeπ AR ' c − ρπ S(w! +Ueq)cosα 4 • After linearization: = − ρ + α + α − 1 ρ 2 − α + α Za Zae (Ueu Wew)S(CLe cos e CDe sin e ) V0 S( CLe sin e CDe cos e )w 2Ue " % 1 2 a 2aCLe c − ρV0 S$ (w +(h1 − h)q)cosαe + (w +(h1 − h)q)sinαe '− ρπ S(w! +Ueq)cosαe 2 #Ue Ueeπ AR & 4 18 Z stability derivatives • Differentiating Za with respect to w, u, q etc yields the Z stability derivatives as: ! = −ρ α + α Zu UeS(CLe cos e CDe sin e ) ! = −ρ α + α − 1 ρ 2 − α + α Zw WeS(CLe cos e CDe sin e ) V0 S( CLe sin e CDe cos e ) 2Ue " % 1 2 a 2aCLe − ρV0 S$ cosαe + sinαe ' 2 #Ue Ueeπ AR & " 2aC % ! 1 2 a Le c Zq = − ρV0 S$ (h1 − h)cosαe + (h1 − h)sinαe '− ρπ SUe cosαe 2 #Ue Ueeπ AR & 4 ! c Zw" = −ρπ S cosαe 4 19 Force Xa • The total force Xa becomes: # & 1 2 1 2 a 2aCL = ρ α − α + ρ % + − α − e + − α( Xa Q∞S(CLe sin CDe cos ) Q∞S (w (h1 h)q)sin (w (h1 h)q)cos 2 2 $Ue Ueeπ AR ' c + ρπ S(w! +Ueq)sinα 4 • After linearization: = + ρ + α − α + 1 ρ 2 α + α Xa Xae (Ueu Wew)S(CLe sin e CDe cos e ) V S(CLe cos e CDe sin e )w 2Ue " % 1 2 a 2aCLe c + ρV0 S$ (w +(h1 − h)q)sinαe − (w +(h1 − h)q)cosαe '+ρπ S(w! +Ueq)sinαe 2 #U U eπ AR & 4 e e 20 X stability derivatives • Differentiating Xa with respect to w, u, q etc yields the X stability derivatives as: ! = ρ α − α Xu UeS(CLe sin e CDe cos e ) ! = ρ α − α + 1 ρ 2 α + α Xw WeS(CLe sin e CDe cos e ) V S(CLe cos e CDe sin e ) 2Ue " % 1 2 a 2aCLe + ρV0 S$ sinαe − cosαe ' 2 #Ue Ueeπ AR & " 2aC % ! 1 2 a Le c Xq = ρV0 S$ (h1 − h)sinαe − (h1 − h)cosαe '+ ρπ SUe sinαe 2 #Ue Ueeπ AR & 4 ! c Xw" = ρπ Ssinαe 4 21 Moment Ma • The total moment Ma becomes: # # ∂ε && = 1 ρ 2 + 1 ρ 2 − − a1 − a + − M a Q∞ScCme Q∞Sc %h h0 VT %1 (( (w (h1 h)q) 2 2 $ a $ ∂α ''Ue 1 2 ∂ε lt 1 2 − ρQ∞ScVT a1 q − ρV ScVT a2η 2 ∂α Ue 2 c # c & c #3c & π 2 + ρπ S%h − (w! − ρπ S% − h(Ueq − ρUec Sq 4 $ 2 ' 4 $ 4 ' 16 • After linearization: " " ∂ε %% = + ρ + + 1 ρ 2 − − a1 − a + − M a M ae (Ueu Wew)ScCme V0 Sc $h h0 VT $1 '' (w (h1 h)q) 2 # a # ∂α &&Ue 1 2 ∂ε lt 1 2 − ρV0 ScVT a1 q − ρV0 ScVT a2η 2 ∂α Ue 2 c " c % c "3c % π 2 + ρπ S$h − 'w! − ρπ S$ − h'Ueq − ρUec Sq 22 4 # 2 & 4 # 4 & 16 M stability derivatives • Differentiating Ma with respect to w, u, q etc yields the M stability derivatives as: ! = ρ Mu UeScCme " " ∂ε %% ! = ρ + 1 ρ 2 − − a1 − a M w WeScCme V0 Sc $h h0 VT $1 '' 2 # a # ∂α &&Ue " " %% ! 1 2 a1 ∂ε a M q = ρV0 Sc $h − h0 −VT $1− '' (h1 − h) 2 # a # ∂α &&Ue " % 1 2 ∂ε lt c 3c π 2 − ρV0 ScVT a1 − ρπ S$ − h'Ue − ρUec S 2 ∂α Ue 4 # 4 & 16 ! c " c % M w" = ρπ S$h − ' 4 # 2 & 23 ! 1 2 Mη = − ρV0 ScVT a2 2 Discussion of longitudinal derivatives • Notice that this methodology for the stability derivatives is by no means exact.
Load more

Recommended publications
  • Introduction to Aircraft Stability and Control Course Notes for M&AE 5070
    Introduction to Aircraft Stability and Control Course Notes for M&AE 5070 David A. Caughey Sibley School of Mechanical & Aerospace Engineering Cornell University Ithaca, New York 14853-7501 2011 2 Contents 1 Introduction to Flight Dynamics 1 1.1 Introduction....................................... 1 1.2 Nomenclature........................................ 3 1.2.1 Implications of Vehicle Symmetry . 4 1.2.2 AerodynamicControls .............................. 5 1.2.3 Force and Moment Coefficients . 5 1.2.4 Atmospheric Properties . 6 2 Aerodynamic Background 11 2.1 Introduction....................................... 11 2.2 Lifting surface geometry and nomenclature . 12 2.2.1 Geometric properties of trapezoidal wings . 13 2.3 Aerodynamic properties of airfoils . ..... 14 2.4 Aerodynamic properties of finite wings . 17 2.5 Fuselage contribution to pitch stiffness . 19 2.6 Wing-tail interference . 20 2.7 ControlSurfaces ..................................... 20 3 Static Longitudinal Stability and Control 25 3.1 ControlFixedStability.............................. ..... 25 v vi CONTENTS 3.2 Static Longitudinal Control . 28 3.2.1 Longitudinal Maneuvers – the Pull-up . 29 3.3 Control Surface Hinge Moments . 33 3.3.1 Control Surface Hinge Moments . 33 3.3.2 Control free Neutral Point . 35 3.3.3 TrimTabs...................................... 36 3.3.4 ControlForceforTrim. 37 3.3.5 Control-force for Maneuver . 39 3.4 Forward and Aft Limits of C.G. Position . ......... 41 4 Dynamical Equations for Flight Vehicles 45 4.1 BasicEquationsofMotion. ..... 45 4.1.1 ForceEquations .................................. 46 4.1.2 MomentEquations................................. 49 4.2 Linearized Equations of Motion . 50 4.3 Representation of Aerodynamic Forces and Moments . 52 4.3.1 Longitudinal Stability Derivatives . 54 4.3.2 Lateral/Directional Stability Derivatives .
    [Show full text]
  • Determination of Static and Dynamic Stability Derivatives Using Beggar
    Air Force Institute of Technology AFIT Scholar Theses and Dissertations Student Graduate Works 3-12-2008 Determination of Static and Dynamic Stability Derivatives Using Beggar Michael E. Bartowitz Follow this and additional works at: https://scholar.afit.edu/etd Part of the Aerodynamics and Fluid Mechanics Commons Recommended Citation Bartowitz, Michael E., "Determination of Static and Dynamic Stability Derivatives Using Beggar" (2008). Theses and Dissertations. 2673. https://scholar.afit.edu/etd/2673 This Thesis is brought to you for free and open access by the Student Graduate Works at AFIT Scholar. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of AFIT Scholar. For more information, please contact [email protected]. Determination of Static and Dynamic Stability Coefficients Using Beggar THESIS Michael E. Bartowitz, Second Lieutenant, USAF AFIT/GAE/ENY/08-M02 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE INSTITUTE OF TECHNOLOGY Wright-Patterson Air Force Base, Ohio APPROVED FOR PUBLIC RELEASE; DISTRIBUTION UNLIMITED. The views expressed in this thesis are those of the author and do not reflect the official policy or position of the United States Air Force, Department of Defense, or the United States Government. AFIT/GAE/ENY/08-M02 Determination of Static and Dynamic Stability Coefficients Using Beggar THESIS Presented to the Faculty Department of Aeronautics and Astronautics Graduate School of Engineering and Management Air Force Institute of Technology Air University Air Education and Training Command In Partial Fulfillment of the Requirements for the Degree of Master of Science in Aeronautical Engineering Michael E. Bartowitz, B.S.A.E.
    [Show full text]
  • Mmittee for Aeronautics, Langley Field, Va., June 21, 1954
    METHODS OF PREDICTING HESLICOPTER STABILITY By Robert J. Tapscott aad F. B. Gustafson Langley Aeronautical Laboratory Langley Field, Va. CLASSIFICATION CHAhGED 10 h BY AUTHORITY OF S 3 w material contains iniormatlon affe United States witMn the me- of the espionage laws, Title 18, U sion or revelatloo of which in arqr mammr to an unautbrized person i I f MMITTEE Q: FOR AERONAUTICS NACA RM L34GO3 NATIONAL ADVISORY COWLITTEE FOR AERONAUTICS METHODS OF PREDICTING IIELICOPTER STABILITY By Robert J. Tapscott and F. B. Gustafson Some of the methods of predicting rotor stability derivatives have been reviewed. The methods by which these rotor derivatives are employed to estimate helicopter stability characteristics have been summarized. It is concluded that, although these methods are not always feasible for predicting absolute values of the stability of the helicopter, the effects on stability of changes in individual derivatives can generally be estimated satisfactorily. INTRODUCTION In order to predict helicopter stability - for example, to esti- mate theoretically whether a prospective helicopter will meet the flying-qualities requirements - both the applicable equations of motion and the necessary stability derivatives must be determined. The processes for handling equations of motion have been well established in conjunction with airplanes, and the modification of these procedures for helicopter use has been found a secondary problem in comparison with the provision of values of stability derivatives. The prediction of stability derivatives, in general, requires a knowledge of the contributions of both the rotor and the fuselage. Fuselage characteristics are not open to as specific an analysis as rotor characteristics, and preliminary estimates can be handled on the basis of data from previous designs and from wind-tunnel model tests.
    [Show full text]
  • Longitudinal Stability and Control Derivatives of a Jet Fighter Airplane Extracted from Flight Test Data by Utilizing Maximum Likelihood Estimation
    https://ntrs.nasa.gov/search.jsp?R=19720010363 2020-03-23T14:20:50+00:00Z View metadata, citation and similar papers at core.ac.uk brought to you by CORE 1 provided by NASA Technical Reports Server N A S A T E C H N I C A L N O T E QP^QP^ NASA TN D-6532 c^1 ^'^&-A/ S ^t-^^ ^s~ m Q t-^ g 5 yy.- RETl/S -S s ^^ " s ^^r^'I. vi. LONGITUDINAL STABILITY AND CONTROL DERIVATIVES OF A JET FIGHTER AIRPLANE EXTRACTED FROM FLIGHT TEST DATA BY UTILIZING MAXIMUM LIKELIHOOD ESTIMATION by George G. Steinmetz, Russell V. Parrish, and Roland L. Bowles / '*>'.'.' ^'''.. / '''... ^' . Langley Research Center \ '. ' '/'' / Hampton, Va. 23365 r/ NATIONAL AERONAUTICS AND SPACE ADMINISTRATION WASHINGTON, D. C. MARCH 1972 f TECH LIBRAHY KAFB, NM D1333b^ 1. Report No. 2. Government Accession No. 3. Recipient's Catalog No. NASA TN D-6532 4. Title and Subtitle 5. Report Date LONGITUDINAL STABILITY AND CONTROL DERIVATIVES March 1972__________ OF A JET FIGHTER AIRPLANE EXTRACTED FROM g performing Organization Code------- FLIGHT TEST DATA BY UTILIZING MAXIMUM LIKELIHOOD ESTIMATION 7. Author(s) 8. Performing Organization Report No. George G. Steinmetz, Russell V. Parrish, L-7882 and Roland L. Bowles 10. Work Unit No. 9. Performing Organization Name and Address 136-62-01-03 .NASA Langley Research Center n. contract Grant No. Hampton, Va. 23365 13. Type of Report and Period Covered 12. Sponsoring Agency Name and Address Technical Note National Aeronautics and Space Administration 14. sponsoring Agency code Washington, D.C. 20546 15. Supplementary Notes 16.
    [Show full text]
  • Determination of Stability and Control Derivatives for a Modern Light Composite Twin Engine Airplane
    Theses - Daytona Beach Dissertations and Theses Fall 2009 Determination of Stability and Control Derivatives for a Modern Light Composite Twin Engine Airplane Monica M. Londono Embry-Riddle Aeronautical University - Daytona Beach Follow this and additional works at: https://commons.erau.edu/db-theses Part of the Aerodynamics and Fluid Mechanics Commons, and the Aviation Commons Scholarly Commons Citation Londono, Monica M., "Determination of Stability and Control Derivatives for a Modern Light Composite Twin Engine Airplane" (2009). Theses - Daytona Beach. 124. https://commons.erau.edu/db-theses/124 This thesis is brought to you for free and open access by Embry-Riddle Aeronautical University – Daytona Beach at ERAU Scholarly Commons. It has been accepted for inclusion in the Theses - Daytona Beach collection by an authorized administrator of ERAU Scholarly Commons. For more information, please contact [email protected]. DETERMINATION OF STABILITY AND CONTROL DERIVATIVES FOR A MODERN LIGHT COMPOSITE TWIN ENGINE AIRPLANE by Monica M. Londono A Thesis Submitted to the Graduate Studies Office In Partial Fulfillment of the Requirements for the Degree of Master of Science in Aerospace Engineering Embry-Riddle Aeronautical University Daytona Beach, Florida Fall 2009 UMI Number: EP31998 INFORMATION TO USERS The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleed-through, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion.
    [Show full text]
  • X-29A Lateral–Directional Stability and Control Derivatives Extracted from High-Angle-Of- Attack Flight Data
    NASA Technical Paper 3664 December 1996 X-29A Lateral–Directional Stability and Control Derivatives Extracted From High-Angle-of- Attack Flight Data Kenneth W. Iliff Kon-Sheng Charles Wang NASA Technical Paper 3664 1996 X-29A Lateral–Directional Stability and Control Derivatives Extracted From High-Angle-of- Attack Flight Data Kenneth W. Iliff Dryden Flight Research Center Edwards, California Kon-Sheng Charles Wang SPARTA, Inc. Lancaster, California National Aeronautics and Space Administration Office of Management Scientific and Technical Information Program CONTENTS Page ABSTRACT . 1 NOMENCLATURE . 1 Acronyms . 1 Symbols . 1 INTRODUCTION . 3 FLIGHT PROGRAM OVERVIEW . 4 VEHICLE DESCRIPTION . 5 INSTRUMENTATION AND DATA ACQUISITION . 7 PREDICTIONS AND ENVELOPE EXPANSION METHODS . 8 PARAMETER IDENTIFICATION METHODOLOGY . 9 RESULTS AND DISCUSSION. 13 Sideslip Derivatives . 14 Aileron Derivatives. 15 Rudder Derivatives . 16 Rotary Derivatives . 16 Aerodynamic and Instrumentation Biases . 17 CONCLUDING REMARKS . 18 REFERENCES . 19 FIGURES 1. X-29A aircraft, number 2 . 22 2. Three-view drawing of the X-29A showing major dimensions . 23 3. The maximum likelihood estimation concept with state and measurement noise . 24 4. Time-history data for a typical high-AOA, lateral–directional, PID maneuver . 25 5. Sideslip derivatives as functions of AOA . 28 6. Aileron derivatives as functions of AOA . 29 7. Rudder derivatives as functions of AOA . 31 8. Rotary derivatives as functions of AOA . 32 9. Aerodynamic and instrumentation biases as functions of AOA . 34 iii ABSTRACT The lateral–directional stability and control derivatives of the X-29A number 2 are extracted from flight data over an angle-of-attack range of 4° to 53° using a parameter identification algorithm.
    [Show full text]
  • Computation of Aircraft Stability Derivatives Using an Automatic Differentiation Adjoint Approach
    AIAA JOURNAL Vol. 49, No. 12, December 2011 Computation of Aircraft Stability Derivatives Using an Automatic Differentiation Adjoint Approach Charles A. Mader∗ University of Toronto, Toronto, Ontario M3H 5T6, Canada and Joaquim R. R. A. Martins† University of Michigan, Ann Arbor, Michigan 48109 DOI: 10.2514/1.J051147 This paper presents a method for the computation of the static and dynamic stability derivatives of arbitrary aircraft configurations. Three-dimensional computational fluid dynamics are used in this method to simulate the flow characteristics around the configuration, and a moving-grid formulation is included in the flow solver to handle the rotational physics necessary for the computation of the dynamic derivatives. To obtain the stability derivatives, the computational fluid dynamics code is differentiated using the automatic differentiation adjoint (ADjoint) approach. This approach enables the efficient and accurate computation of derivatives for a wide variety of variables, including the dynamic model states that are typical of the stability derivatives. To demonstrate the effectiveness of this approach, stability derivatives are computed for a NACA 0012 airfoil and an ONERA M6 wing. Nomenclature v = flow velocity with respect to the moving grid b = span V = aircraft speed w = velocity of the moving grid CD = aircraft drag coefficient fi x = grid coordinates (Sec. III.A), design variables (Sec. III.C) CL = aircraft lift coef cient xref = center of moment and rotation with respect to wing Cl = airfoil lift coefficient (Sec. III.B), aircraft roll moment coefficient (Sec. IV) root/airfoil leading edge fi = angle of attack Cm = airfoil pitch moment coef cient (Sec. III.B), aircraft pitch fi = sideslip angle moment coef cient (Sec.
    [Show full text]
  • A COMPREHENSIVE REVIEW of VERTICAL TAIL DESIGN 1 Introduction
    A COMPREHENSIVE REVIEW OF VERTICAL TAIL DESIGN Fabrizio Nicolosi Danilo Ciliberti Pierluigi Della Vecchia Associate Professor Post-Doc Assistant Professor University of Naples University of Naples University of Naples “Federico II” – Dept. of “Federico II” – Dept. of Industrial “Federico II” – Dept. of Industrial Industrial Engineering (DII) Engineering (DII) Engineering (DII) Via Claudio 21, 80125, Via Claudio 21, 80125, Via Claudio 21, 80125, Naples – ITALY Naples – ITALY Naples – ITALY [email protected] [email protected] [email protected] Salvatore Corcione Vincenzo Cusati Post-Doc Consultant University of Naples “Federico II” – University of Naples “Federico II” – Dept. of Industrial Engineering (DII) Dept. of Industrial Engineering (DII) Via Claudio 21, 80125, Via Claudio 21, 80125, Naples – ITALY Naples – ITALY [email protected] [email protected] Abstract. This work deals with a comprehensive review of vertical tail design methods for aircraft directional stability and vertical tail sizing. The focus on aircraft directional stability is due to the significant discrepancies that classical semi-empirical methods, as USAF DATCOM and ESDU, provide for some configurations, since they are based on NACA wind tunnel tests about models not representative of an actual transport airplane. The authors performed RANS CFD simulations to calculate the aerodynamic interference among aircraft parts for hundreds configurations of a generic regional turboprop aircraft, providing useful results that have been collected in a new vertical tail preliminary design method, named VeDSC. Semi-empirical methods have been put in comparison on a regional turboprop aircraft, where the VeDSC method shows a strong agreement with numerical results. A wind tunnel investigation involving more than 180 configurations has validated the numerical approach.
    [Show full text]
  • Flight Dynamics Summary
    Flight Dynamics Summary 1. Introduction In this summary we examine the flight dynamics of aircraft. But before we do that, we must examine some basic ideas necessary to explore the secrets of flight dynamics. 1.1 Basic concepts 1.1.1 Controlling an airplane To control an aircraft, control surfaces are generally used. Examples are elevators, flaps and spoilers. When dealing with control surfaces, we can make a distinction between primary and secondary flight control surfaces. When primary control surfaces fail, the whole aircraft becomes uncontrollable. (Examples are elevators, ailerons and rudders.) However, when secondary control surfaces fail, the aircraft is just a bit harder to control. (Examples are flaps and trim tabs.) The whole system that is necessary to control the aircraft is called the control system. When a control system provides direct feedback to the pilot, it is called a reversible system. (For example, when using a mechanical control system, the pilot feels forces on his stick.) If there is no direct feedback, then we have an irreversible system. (An example is a fly-by-wire system.) 1.1.2 Making assumptions In this summary, we want to describe the flight dynamics with equations. This is, however, very difficult. To simplify it a bit, we have to make some simplifying assumptions. We assume that . • There is a flat Earth. (The Earth’s curvature is zero.) • There is a non-rotating Earth. (No Coriolis accelerations and such are present.) • The aircraft has constant mass. • The aircraft is a rigid body. • The aircraft is symmetric. • There are no rotating masses, like turbines.
    [Show full text]
  • Calculation of Lateral-Directional Stability Derivatives for Wing-Body Combinations with and Without Jet-Interaction Effects
    https://ntrs.nasa.gov/search.jsp?R=19770023140 2020-03-22T07:24:50+00:00Z View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by NASA Technical Reports Server N 7 7 - 3 G 0 8 4 NASA Contractor Report 145251 CALCULATION OF LATERAL-DIRECTIONAL STABILITY DERIVATIVES FOR WING-BODY COMBINATIONS WITH AND WITHOUT JET-INTERACTION EFFECTS C. EDWARD LAN NASA GRANT NSG-1139 THE UNIVERSITY OF KANSAS CENTER FOR RESEARCH, INC, LAWRENCE, KS 66045 AUGUST 1977 NASA National Aeronautics and Space Administration Langley Research Center Hampton, Virginia 23665 Calculation of Lateral-Directional Stability Derivatives for Wing-Body Combinations with and without Jet-Interaction Effects by C. Edward Lan Technical Report CRINC-FRL-281-1 August 1977 The Flight Research Laboratory The University of Kansas Center for Research, Inc. Lawrence, Kansas 66045 Prepared under NASA Grant NSG 1139 for Langley Research Center National Aeronautics and Space Administration Table of Contents Page Summary 1 1. List of Symbols 2 2. Introduction 6 3. Theoretical Formulation 3.1 Fuselage Effect without Jet Interaction 7 3.2 Fuselage Effect with Jet Interaction 15 3.3 Forces and Moments 16 4. Numerical Results and Discussions 4.1 Fuselage Alone 29 4.2 Wing or Wing-Body Characteristics without Jet Interaction 29 4.3 Effect of Jet Interaction 36 5. Concluding Remarks 43 6. References 44 Appendix A Influence Coefficient Matrices for Wing-Body- Jet Interaction 48 Appendix B Derivation of fn(x£^) in the formulation for the Fuselage effect 55 Appendix C Calculation of Fuselage Aerodynamic Characteristics 59 Appendix D A Trigonometric Interpolation Formula for the Calculation of Streamwise Vortex Density 53 Appendix E Calculation of Tip Suction 68 Summary A theoretical method is presented for predicting the lateral- directional stability derivatives of wing-body combinations with or without the blowing jet effect.
    [Show full text]
  • Methodology to Analyze Attitude Stability Of
    Methodology to analyze attitude stability of satellites subjected to aerodynamic torques David Mostaza-Prieto1 and Peter C. E. Roberts2 The University of Manchester; Manchester; EnglandM139P L; UnitedKingdom: A method to analyze the stability of satellites flying in low Earth orbits is presented. Starting from the classical rotational equations of motion including gravity gradient, the influence of aerodynamic torques is studied both in the pitch and roll-yaw axes. Using concepts used to assess aircraft stability, such as stability derivatives, several ideas and recommendations are presented. The work aims to link geometrical prop- erties with attitude response, providing a way to design satellites that actively take into account environmental torques, either by minimizing them or by utilizing them to meet mission goals. 1 Researcher, School of Mechanical, Aerospace and Civil Engineering. 2 Lecturer, School of Mechanical, Aerospace and Civil Engineering. 1 Nomenclature 2 Ix,Iy,Iz = principal moments of inertia [kg:m ] !0 = angular orbital velocity in a circular orbit [rad=s] φ = roll angle [rad] θ = pitch angle [rad] = yaw angle [rad] M = disturbance torque [N:m] Va = velocity w.r.t wind [m=s] ρ = density [kg=m3] 2 Sref = aerodynamic reference area [m ] Lref = aerodynamic reference length [m] Cl = roll moment coefficient Cm = pitch moment coefficient Cn = yaw moment coefficient C~F = aerodynamic force coefficients in body axes α = angle of attack [rad] β = angle of sideslip [rad] δ1, δ2 = angles defining the orientation of aerodynamic velocity in orbital frame [rad] L = vehicle length [m] D = diameter of the vehicle base area [m2] xcog = distance from the leading face to the center of mass [m] CD = drag coefficient Cmα = stability derivative, slope of the pitch coefficient w.r.t angle of attack Torb = orbital period I.
    [Show full text]
  • STABILITY DERIVATIVES of CONESAT by Murray Tobak
    I STABILITY DERIVATIVES OF CONESAT SUPERSONIC SPEEDS By Murray Tobak and William R. Wehre.nd Ames Aeronautical Laboratory MoHett Field, Calif. Washington September 1956 \ —. --- —. =..-.. .> ..= . .. -- - IECHUBMY KAFB,NM lY NATIONAL ADVISORY COMMITHW FOR AERONAUT‘ Iulllill!lulllllulllll DQLL703 mcm%ccwmm 3788 . STABILITYDERIVATIVES OF CONES AT SUPERSONIC SPEEDS By Murray Tobak and WillMm R. Webrend , The aerodynamic stability derivatives due to pitching velocity and vertical acceleration are calculated by use of potential theory for cir- cular cones traveling at supersonic speeds. The analysis is based on two theoretical techniques used successfdly previously in application to the case of uniform axial and inclined flow. ~ the ftist, potential solu- tions for sxial flow snclcrossfl.oware deriveclfrm the first-order wave equation but in application to caltiations for the forces no approxima- # tions sre made either to the tangency condition or to the isentropic pressure relatim. The secona methoa consists in conibmg the first- . oraer crossflow potential with an axial-flow potential correct to second or&?r. Closed-form solutions by both methods are found for a cone,”and numerical results for the stability derivatives are pres~teii as a func- tion of Mach number for cones having semivertex angles of 10° anti20°. lh addition, expressions for the forces, maments, and stability derivatives of arbitrary bodies of revolution are obtainea using Newtonian impact theory. Numerical results for cones compsre well with those obtainea from the conibinedfirst- ana second-order potential theory at the highest Mach nuniberfor which the latter theory is applicable. INTRODUCTION The hportance of the body as a lift-producing congmnent of aircraft flying at supersonic speea has occasioned a great aea,lof theoretical work frm which it is now possible to calculate the bodyts static aerodynamic properties with good accuracy (see, e.g., refs.
    [Show full text]