Aerodynamic Modeling of an Unmanned Aerial Vehicle Using a Computational Fluid Dynamics Prediction Code

Total Page:16

File Type:pdf, Size:1020Kb

Aerodynamic Modeling of an Unmanned Aerial Vehicle Using a Computational Fluid Dynamics Prediction Code AERODYNAMIC MODELING OF AN UNMANNED AERIAL VEHICLE USING A COMPUTATIONAL FLUID DYNAMICS PREDICTION CODE A thesis presented to the faculty of the Russ College of Engineering and Technology of Ohio University In partial fulfillments of the requirements for the degree Master of Science Isaac D. Rose March 2009 © 2009 Isaac D. Rose. All Rights Reserved. This thesis titled AERODYNAMIC MODELING OF AN UNMANNED AERIAL VEHICLE USING A COMPUTATIONAL FLUID DYNAMICS PREDICTION CODE by ISAAC D. ROSE has been approved for the School of Electrical Engineering and Computer Science and the Russ College of Engineering and Technology by _____________________________________________________________ Douglas A. Lawrence Professor of Electrical Engineering and Computer Science _____________________________________________________________ Dennis Irwin Dean, Russ College of Engineering and Technology Abstract ROSE, ISAAC D., M.S., March 2009, Electrical Engineering AERODYNAMIC MODELING OF AN UNMANNED AERIAL VEHICLE USING A COMPUTATIONAL FLUID DYNAMICS PREDICTION CODE (192 pp.) Director of Thesis: Douglas A. Lawrence The process of creating a six degree-of-freedom model for an aerospace vehicle requires detailed knowledge of the aerodynamic characteristics. This thesis presents an implementation of a Computational Fluid Dynamics (CFD) prediction computer code to generate aerodynamic coefficients for the Brumby Mk. I Unmanned Aerial Vehicle (UAV). The aerodynamic coefficients include both the force and moment coefficients. These values are verified by creating a Matlab/Simulink six degree-of-freedom model. Approved: ______________________________________________________________ Douglas A. Lawrence Professor of Electrical Engineering and Computer Science Acknowledgments I would like to thank GOD through whom all things are possible. I would like to thank my wife and son for their understanding and encouragement. The hours spent working on this thesis were hours spent away from them. I would also like to thank Dr. Lawrence for the guidance and direction that he has given me over the years. The research presented in this thesis is a tribute to his resolve in the autonomous control of the Brumby Unmanned Arial Vehicle. Finally, I would like to thank the faculty of the department of Electrical Engineering and Computer Science for their help throughout the years. Table of Contents Abstract.................................................................................................................................3 Acknowledgments................................................................................................................4 Glossary of Variables..........................................................................................................14 Chapter 1: Introduction..................................................................................................19 1.1 Overview.............................................................................................................19 1.2 Motivation ..........................................................................................................20 1.3 Modeling Aerodynamic Forces and Moments....................................................21 1.4 Objectives............................................................................................................22 1.5 Thesis Organization.............................................................................................22 1.5.1 Missile DATCOM Input Parameters...........................................................22 1.5.2 Missile DATCOM Model of the Brumby UAV..........................................23 1.5.3 Equations of Motion....................................................................................23 1.5.4 Brumby UAV Model Simulation................................................................23 Chapter 2: Missile DATCOM Modeling Parameters Overview....................................24 2.1 Flight Conditions..................................................................................................26 2.2 Fuselage...............................................................................................................28 2.3 Primary Lifting Surface ......................................................................................29 2.4 Horizontal Stabilizer...........................................................................................33 2.5 Vertical Stabilizer................................................................................................35 2.6 Control Surfaces..................................................................................................36 2.7 Generating Additional Data................................................................................38 2.8 File Format and Content......................................................................................39 2.9 Missile DATCOM Example...............................................................................41 Chapter 3: Missile DATCOM Model of the Brumby Unmanned Aerial Vehicle .........72 3.1 Flight Conditions.................................................................................................72 3.2 Fuselage................................................................................................................75 3.3 Wing Planform....................................................................................................76 3.4 Vertical Stabilizer................................................................................................77 3.5 Control Surfaces .................................................................................................80 Chapter 4: Equations of Motion and Rigid Body Modeling .........................................86 4.1 Equations of Motion for A Rigid Body ..............................................................86 4.2 Aerodynamic Coefficients..................................................................................91 4.3 Six Degree-of-Freedom Aircraft Model............................................................101 Chapter 5: Simulation...................................................................................................103 5.1 Simulink Nonlinear Aircraft Model..................................................................103 5.2 Trimmed Aircraft Flight....................................................................................107 5.3 Linearized Aircraft Model.................................................................................108 5.4 Nonlinear Simulation Results...........................................................................112 5.4 Control Surface Doublet Simulation Results....................................................123 Chapter 6: Conclusions and Future Work....................................................................138 References........................................................................................................................140 Appendix A.1: for005.dat File.........................................................................................141 Appendix A.2 : Truncated for006.dat File.......................................................................145 Appendix A.3 : for003.dat File........................................................................................154 Appendix A.4 : for021.dat File........................................................................................155 Appendix B.1: Equations of Motion s-function...............................................................184 Appendix B.2: forces_moments.m..................................................................................189 Appendix B.3: datcomderive.m.......................................................................................192 List of Tables Table 2.1: Missile DATCOM Control Cards.....................................................................26 Table 2.2: Missile DATCOM Namelist FLTCON ..........................................................27 Table 2.3: Missile DATCOM Namelist REFQ..................................................................28 Table 2.4: Missile DATCOM Namelist AXIBOD............................................................29 Table 2.5: Missile DATCOM Namelist FINSET..............................................................30 Table 2.6: Missile DATCOM Namelist DEFLCT.............................................................38 Table 2.7: Missile DATCOM File Definitions..................................................................41 Table 3.1: Brumby UAV Flight Conditions (FLTCON) ..................................................73 Table 3.2: Brumby UAV Reference Values (REFQ)........................................................75 Table 3.3: Brumby UAV Body Definition (ASYM).........................................................76 Table 3.4: Brumby UAV Twin Vertical Tail Planform Definition (FINSET2)................79 Table 4.5: Brumby UAV Wing Planform Definition (FINSET1).....................................88 Table 5.1: Brumby UAV Mass Properties.......................................................................104 Table 5.2: S-function Functionality.................................................................................104 Table 5.3: DATCOMTableMex.dll Functionality...........................................................106 Table 5.4: Brumby UAV Control Input Trimmed Values (Case 1)................................107
Recommended publications
  • 6. Aerodynamic-Center Considerations of Wings I
    m 6. AERODYNAMIC-CENTER CONSIDERATIONS OF WINGS I AND WING-BODY COMBINATIONS By John E. Lsmar and William J. Afford, Jr. NASA Langley Research Center SUMMARY Aerodynamic-center variations with Mach number are considered for wings of different planform. The normalizing parameter used is the square root of the wing area, which provides a more meaningful basis for comparing the aerodynamic-center shifts than does the mean geometric chord. The theoretical methods used are shown to be adequate for predicting typical aerodynamlc-center shifts, and ways of minimizing the shifts for both fixed and variable-sweep wings are presented. INTRODUCTION In the design of supersonic aircraft, a detailed knowledge of the aerodynamlc-center position is important in order to minimize trim drag, maximize load-factor capability, and provide acceptable handling qualities. One of the principal contributions to the aerodynamic-center movement is the well-known change in load distributionwithMach number in going from subsonic to supersonic speeds. In addition, large aerodynamic-center var- iations are quite often associated with variable-geometry features such as variable wing sweep. The purpose of this paper is to review the choice of normalizing param- eters and the effects of Mach number on the aerodynamic-center movement of rigid wing-body combinations at low lift. For fixed wings the effects of both conventional and composite planforms on the aerodynamic-center shift are presented, and for variable-sweepwings the characteristic movements of aerodynamlc-center position wlth pivot location and with variable-geometry apex are discussed. Since systematic experimental investigations of the effects of planform on the aerodynamic-center movement with Mach number are still limited, the approach followed herein is to establish the validity of the computative processes by illustrative comparison wlth experiment and then to rely on theory to show the systematic variations.
    [Show full text]
  • CHAPTER TWO - Static Aeroelasticity – Unswept Wing Structural Loads and Performance 21 2.1 Background
    Static aeroelasticity – structural loads and performance CHAPTER TWO - Static Aeroelasticity – Unswept wing structural loads and performance 21 2.1 Background ........................................................................................................................... 21 2.1.2 Scope and purpose ....................................................................................................................... 21 2.1.2 The structures enterprise and its relation to aeroelasticity ............................................................ 22 2.1.3 The evolution of aircraft wing structures-form follows function ................................................ 24 2.2 Analytical modeling............................................................................................................... 30 2.2.1 The typical section, the flying door and Rayleigh-Ritz idealizations ................................................ 31 2.2.2 – Functional diagrams and operators – modeling the aeroelastic feedback process ....................... 33 2.3 Matrix structural analysis – stiffness matrices and strain energy .......................................... 34 2.4 An example - Construction of a structural stiffness matrix – the shear center concept ........ 38 2.5 Subsonic aerodynamics - fundamentals ................................................................................ 40 2.5.1 Reference points – the center of pressure..................................................................................... 44 2.5.2 A different
    [Show full text]
  • An Aerodynamic Comparative Analysis of Airfoils for Low-Speed Aircrafts
    Published by : International Journal of Engineering Research & Technology (IJERT) http://www.ijert.org ISSN: 2278-0181 Vol. 5 Issue 11, November-2016 An Aerodynamic Comparative Analysis of Airfoils for Low-Speed Aircrafts Sumit Sharma B.E. (Mechanical) MPCCET, Bhilai, C.G., India M.E. (Thermal) SSCET, SSGI, Bhilai, C.G., India Abstract - This paper presents investigation on airfoil S1223, NACA4412 at specified boundary conditions, A S819, S8037 and S1223 RTL at low speeds to find out the most comparison is done between results obtained from CFD suitable airfoil design to be used in the low speed aircrafts. simulation in ANSYS and results obtained from an The method used for analysis of airfoils is Computational experiment done in the wind tunnel. The comparative Fluid Dynamics (CFD). The numerical simulation of low result shows that, the result from experiment and from speed and high-lift airfoil has been done using ANSYS- FLUENT (version 16.2) to obtain drag coefficient, lift CFD simulation showing close agreement between these coefficient, coefficient of moment and Lift-to-Drag ratio over two methods. So, the CFD simulation using ANSYS- the airfoils for the comparative analysis of airfoils. The S1223 FLUENT can be used as an relative alternative to RTL airfoil has been chosen as the most suitable design for experimental method in determining drag and lift. the specified boundary conditions and the Mach number from 0.10 to 0.30. Jasminder Singh et al [2] (June- 2015), Analysis is done on NACA 4412 and Seling 1223 airfoils using Computational Key words: CFD, Lift, Drag, Pitching, Lift-to-Drag ratio.
    [Show full text]
  • Lifting Surface Lift and Pitch Considerations
    Stability and Control Some Characteristics of Lifting Surfaces, and Pitch-Moments The lifting surfaces of a vehicle generally include the wings, the horizontal and vertical tail, and other surfaces such as canards, winglets, etc. What they all have in common is that there purpose is to create an aerodynamic force and moment. Whether the surface is infinite (2-D) or finite (3-D), we can assign to it a force and moment. Usually we designate a 2-D lifting surface as an airfoil, and a 3-D lifting surface as a wing (or tail, or winglet, of whatever). Here all lifting surfaces will be assigned the name - wing. Generally the information concerning forces and moments on a wing are determined from wind tunnel tests or from computer codes. Unlike a force, a moment must be referred to a specific reference point. Often time the information given is not referenced to the point of your interest. Consequently we would like to be able to convert results referenced to one point to equivalent information with respect to a different reference point. We will later use these same ideas applied to a complete aircraft. Consider two identical wings mounted in a wind tunnel. One wing is supported at point A, and the other at point B. At each of these support points there is an electronic “balance” used to measure the forces and moments acting on the wing at those respective points. For simplicity the wings will be represented as flat plates. Si nce the wings are identical we can expect that the aerodynamic forces will be the same.
    [Show full text]
  • Preliminary Design and CFD Analysis of a Fire Surveillance Unmanned Aerial Vehicle Luis E
    TFAWS-08-1034 Preliminary Design and CFD Analysis of a Fire Surveillance Unmanned Aerial Vehicle Luis E. Casas1, Jon M. Hall1, Sean A. Montgomery1, Hiren G. Patel1, Sanjeev S. Samra1, Joe Si Tou1, Omar Quijano2, Nikos J. Mourtos3, Periklis P. Papadopoulos3 Aerospace Engineering San Jose State University One Washington Square San Jose, CA, 95192-0087 The Spartan Phoenix is an Unmanned Aerial Vehicle (UAV) designed for fire surveillance. It was inspired by the summer 2007 wildfires in Greece and California that caused billions of dollars in damage and claimed hundreds of lives. The Spartan Phoenix will map the perimeter of wildfires and provide real time feedback as to fire location and growth pattern. This paper presents the preliminary design of the UAV along with a CFD study of the horizontal and vertical stabilizers of the aircraft. Both surfaces use a NACA 0012 airfoil. First, the flow field over the airfoil was constructed using a CFD-GEOM grid. Next the grid was imported into CFD-FASTRAN and CFD-ACE to simulate a subsonic flow through the grid. The simulation was run using the Spartan Phoenix freestream conditions of 45 m/s airspeed at 00 angle-of-attack, and 20 m/s airspeed at 80 angle-of-attack, with pressure and temperature at 1 atm and 300 K respectively. The generated CFD solutions compared favorably to published data, solutions generated by the Institute of Computational Fluid Dynamics (iCFD), as well as results obtained from airfoil analysis using Sub2D, a potential flow simulation software. Nomenclature AR = aspect ratio = angle of attack CL = lift coefficient CL,min = lift coefficient for minimum drag CD = drag coefficient CDi = induced drag coefficient CD0 = zero lift drag coefficient CG = center of gravity c = chord length e = Oswald efficiency factor L/D = lift-to-drag ratio UAV = Unmanned Aerial Vehicle I.
    [Show full text]
  • 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]
  • Longitudinal Static Stability
    Longitudinal Static Stability Some definitions M pitching moment without dimensions C m 1 (so without influence of ρ, V and S) V2 Sc 2 it is a ‘shape’ parameter which varies with the angle of attack. Note the chord c in the denominator because of the unit Nm! L For the wing+aircraft we use the C L 1 surface area of the wing S! VS2 2 L For the tail we use the surface of the tail: SH ! C H LH 1 VS2 2 H Definition of aerodynamic center of a wing: The aerodynamic center (a.c.) is the point around which the moment does not change when the angle of attack changes. We can therefore use Cmac as a constant moment for all angles of attack. The aerodynamic center usually lies around a quarter chord from the leading edge. 1 Criterium for longitudinal static stability (see also Anderson § 7.5): We will look at the consequences of the position of the center of gravity, the wing and the tail for longitudinal static stability. For stability, we need a negative change of the pitching moment if there is a positive change of the angle of attack (and vice versa), so: CC 0 mm 00 0 Graphically this means Cm(α) has to be descending: For small changes we write: dC m 0 d We also write this as: C 0 m 2 When Cm(α) is descending, the Cm0 has to be positive to have a trim point where Cm = 0 and there is an equilibrium: So two conditions for stability: 1) Cm0 > 0; if lift = 0; pitching moment has to be positive (nose up) dC 2) m 0 ( or C 0 ); pitching moment has to become more negative when d m the angle of attack increases Condition 1 is easy to check.
    [Show full text]
  • Introduction to Aircraft Performance and Static Stability
    16.885J/ESD.35J Aircraft Systems Engineering Introduction to Aircraft Performance and Static Stability Prof. Earll Murman September 18, 2003 Today’s Topics • Specific fuel consumption and Breguet range equation • Transonic aerodynamic considerations • Aircraft Performance – Aircraft turning – Energy analysis – Operating envelope – Deep dive of other performance topics for jet transport aircraft in Lectures 6 and 7 • Aircraft longitudinal static stability Thrust Specific Fuel Consumption (TSFC) lb of fuel burned • Definition: TSFC (lb of thrust delivered)(hour) • Measure of jet engine effectiveness at converting fuel to useable thrust • Includes installation effects such as – bleed air for cabin, electric generator, etc.. – Inlet effects can be included (organizational dependent) • Typical numbers are in range of 0.3 to 0.9. Can be up to 1.5 • Terminology varies with time units used, and it is not all consistent. – TSFC uses hours – “c” is often used for TSFC lb of fuel burned – Another term used is c t (lb of thrust delivered)(sec) Breguet Range Equation • Change in aircraft weight = fuel burned dW ctTdt ct TSPC/3600 T thrust • Solve for dt and multiply by Vf to get ds VfdW VfW dW VfL dW ds Vfdt ctT ctT W ctD W • Set L/D, ct, Vf constant and integrate W 3600 L TO R Vf ln TSFC DW empty Insights from Breguet Range Equation W 3600 L TO R Vf ln TSFC DW empty 3600 represents propulsion effects. Lower TSFC is better TSFC V L represents aerodynamic effect. L/D is aerodynamic efficiency fD V La M L . a is constant above 36,000 ft.
    [Show full text]
  • Static Stability and Control
    CHAPTER 2 Static Stability and Control "lsn't it astonishing that all these secrets have been preserved for so many years just so that we could discover them!" Orville Wright, June 7, 1903 2.1 HISTORICAL PERSPECTIVE By the start of the 20th century, the aeronautical community had solved many of the technical problems necessary for achieving powered flight of a heavier-than-air aircraft. One problem still beyond the grasp of these early investigators was a lack of understanding of the relationship between stability and control as well as the influence of the pilot on the pilot-machine system. Most of the ideas regarding stability and control came from experiments with uncontrolled hand-launched gliders. Through such experiments, it was quickly discovered that for a successful flight the glider had to be inherently stable. Earlier aviation pioneers such as Albert Zahm in the United States, Alphonse Penaud in France, and Frederick Lanchester in England contributed to the notion of stability. Zahm, however, was the first to correctly outline the requirements for static stability in a paper he presented in 1893. In his paper, he analyzed the conditions necessary for obtaining a stable equilibrium for an airplane descending at a constant speed. Figure 2.1 shows a sketch of a glider from Zahm's paper. Zahm concluded that the center of gravity had to be in front of the aerodynamic force and the vehicle would require what he referred to as "longitudinal dihedral" to have a stable equilibrium point. In the terminology of today, he showed that, if the center of gravity was ahead of the wing aerodynamic center, then one would need a reflexed airfoil to be stable at a positive angle of attack.
    [Show full text]
  • Design and Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit
    Washington University in St. Louis Washington University Open Scholarship Engineering and Applied Science Theses & Dissertations McKelvey School of Engineering Spring 5-19-2017 Design and Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit Maria E. Ferguson Washington University in St Louis Follow this and additional works at: https://openscholarship.wustl.edu/eng_etds Part of the Aerodynamics and Fluid Mechanics Commons Recommended Citation Ferguson, Maria E., "Design and Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit" (2017). Engineering and Applied Science Theses & Dissertations. 227. https://openscholarship.wustl.edu/eng_etds/227 This Thesis is brought to you for free and open access by the McKelvey School of Engineering at Washington University Open Scholarship. It has been accepted for inclusion in Engineering and Applied Science Theses & Dissertations by an authorized administrator of Washington University Open Scholarship. For more information, please contact [email protected]. WASHINGTON UNIVERSITY IN ST. LOUIS School of Engineering and Applied Science Department of Mechanical Engineering and Materials Science Thesis Examination Committee: Dr. Ramesh Agarwal, Chair Dr. Qiulin Qu Dr. David Peters Design and Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit by Maria E. Ferguson A thesis presented to the School of Engineering and Applied Science of Washington University in St. Louis in partial fulfillment of the requirements for the degree of Master of Science May
    [Show full text]
  • And High Power, the Dynamic Pressure in the Shaded Area Can Be Much
    NAVWEPS 00-801-80 BASIC AERODYNAMICS and high power, the dynamic pressure in the net lift produced by the airfoil is difference shaded area can be much greater than the free between the lifts on the upper and lower sur- stream and this causes considerably greater faces. The point along the chord where the lift than at zero thrust. At high power con- distributed lift is effectively concentrated is ditions the induced flow also causes an effect termed the "center of pressure, c.p. - The similar to boundary layer control and increases center of pressure is essentially the ''center of the maximum lift angle of attack. The typical gravity" of the distributed lift pressure and four-engine propeller driven airplane may have the location of the c.p. is a function of camber 60 to 80 percent of the wing area affected by and section lift coefficient. the induced flow and power effects on stall Another aerodynamic reference point is the speeds may be considerable. Also, the lift of -aerodynamic center, a.c. - The aerodynamic the airplane at a given angle of attack and air- center is defined as the point along the chord speed will be greatly affected. Suppose the where all changes in lift effectively take place. airplane shown is in the process of landing To visualize the existence of such a point, flare from a power-on approach. If there is notice the change in pressure distribution with a sharp, sudden reduction of power, the air- angle of attack for the symmetrical airfoil plane may drop suddenly because of the reduced of figure 1.21.
    [Show full text]
  • Full-Scale Ornithopter
    The Flight Dynamics of a Full-Scale Ornithopter Tahir Rashid A thesis submitted in confomity with the requiiements for the degree of Master of Applied Science Graduate Department of Aerospace Science and Engineering in the University of Toronto O Copyright by Tahir Rashid, 1995 National Library Bibliothèque nationale du Canada Acquisitions and Acquisitions et Bibliographie Services services bibliographiques 395 Wdliion Street 395. nie Wellington OnewaûN K1AW Ottawa ON K1A ON4 Canaada Canada The author has granted a non- L'auteur a accordé une licence non exclusive licence allowing the exclusive permettant a la National Library of Canada to Bibliothèque nationale du Canada de reproduce, loan, distribute or sel1 reproduire, prêter, distribuer ou copies of this thesis in rnicroform, vendre des copies de cette thèse sous paper or electronic formats. la forme de microfiche/fiim, de reproduction sur papier ou sur format électronique. The author retauis ownership of the L'auteur conserve la propriéte du copyright in this thesis. Neither the droit d'auteur qui protège cette thèse. thesis nor substantial exûacts ffom it Ni la thèse ni des extraits substantiels may be printed or otherwise de celle-ci ne doivent être imprimés reproduced without the author's ou autrement reproduits sans son permission. autorisation. The Flight ûynamics of a FulCScale Ornithopter Master of Applied Science, 1995 Tahir Rashid Aerospace Science and Engineering University of Toronto This t hesis investigates the non-lineaif light dynamics of a full-scale flapping-wing aircraft (ornithopter). Using simplifying assumptions, the equations of motion were developed for a 2-wing-panel and 3-wing-panel rnodel.
    [Show full text]