ME-114 Final Project Timothy Burkhard

ME-114 Final Project Timothy Burkhard - Phil Barat - John Gyurics 1

Background: ¢ Timothy Burhard and John Gyurics have recently finished designing a UAV model aircraft for the AUVSI competition. The design phase will soon be followed by building the composite aircraft, however the design will first be validated performing simulation analysis and also by building a much cheaper balsa version of it. Phil Barat is a Kite-Boarding enthusiast, who spends every waking moment possible out at the lake being tossed around by the forces of the wind, and is also very interested in the aero-dynamics involved in flight. For our Vibration and Controls project we decided to merge our goals and interests, and joined forces to develop a study of the longitudinal stability of an aircraft. 2

Objective: ¢ Perform a longitudinal Dynamic stability study on an aircraft, by working with the transfer functions that relates speed to the elevator in particular. Employ Matlab as our guiding tool to derive the Step. Response, Impulse Response, Root Locus, Bode Diagram, and the Sisotool. In “Airplane Flight Dynamics and Automatic Flight Controls”, Jan Roskam did a superb job at deriving the transfer functions for different style aircraft, and also demonstrating the correlation between a mass-spring-dampener system and a general aircraft dynamics system. Some of our work will utilize the resources made available by Jan, and further investigate and validate the work. 3

Along The Way: ¢ We started out by first analyzing a simple circuit, and generating the Root Locus Diagram, Body Diagram, Step Response Diagram and Impulse Response Diagram. We later turned our attention to Simulink, and performed an overall aircraft control loop analysis. Finally we tackled the Longitudinal Dynamic Stability Study on an aircraft, by performing the Step-Response, Impulse Response, Root Locus, Bode Diagram, and the development of the Sisotool. 4

Electronic Circuit 5

Bond Graph 6

MATLAB OUTPUT ¢ ¢ From the circuit, we enter the parameters of the two inductors, the resistor, and the capacitor in the campgsym. m file These lines need to be unsuppressed as well, to look like what is shown below: I 2 = 2 ; R 4 = 1 ; C 6 = 1 ; I 7 = 3 ; 7

MATLAB OUTPUT ¢ ¢ ¢ Then, the e 7 (the efforts link which represents voltage) bond graph link must be chosen to produce a transfer function, state-space form, and related figures. It will need to look like this: % e 7=P 2/I 2*R 4 -P 7/I 7*R 4 -Q 6/C 6 fprintf ('n Output # 1 is e 7' ) C(1, : ) = [-1/C 6, -1/I 7*R 4, 1/I 2*R 4]; D(1, : ) = [0]; 8

MATLAB OUTPUT ¢ ¢ ¢ ¢ ¢ State-space form: A MATRIX [ 0 1/3 [ [-1 -1/3 [ [ 0 1/3 C MATRIX [-1 -1/3 9 ¢ B MATRIX [0] [] [1] ¢ 0 ] ] 1/2 ] ] -1/2] ¢ ¢ ¢ 1/2] ¢ D MATRIX 0

MATLAB OUTPUT ¢ Characteristic Polynomial ¢ Transfer Functions Matrix H. . . ¢ 2 ¢ ¢ ¢ 3 2 ¢ s + 5/6 s + 1/3 s + 1/6 ¢ 10 s 3 ----------3 2 6 s +5 s +2 s+1

BODE Diagram and Root Locus 11

Step Response 12

Impulse Response 13

Impulse Response ¢ ¢ ¢ ¢ We can also use the SISOTOOL feature in MATLAB to design the control for this circuit. The following can be entered in the MATLAB command screen: num =[1 0 0] num = 1 0 0 >> den=[6 5 2 1] den = 6 5 2 1 >> g=tf(num, den) ¢ ¢ ¢ Transfer function: s^2 -----------6 s^3 + 5 s^2 + 2 s + 1 ¢ ¢ >> sisotool(g) 14

Impulse Response ¢ By moving the cursor over the plot, and following the directions at the bottom dialogue box, you can adjust the loop gain to design a desired system 15

Final Design Front View 16

Final Design Side View 17

Final Design Top View 18

Final Design Isometric View 19

20

Explanation of Static Stability ¢ Static stability: l l l A statically unstable body will tend to accelerate away from its equilibrium position when disturbed. A neutrally stable body will not accelerate toward or away from any position when disturbed; i. e. , there is no particular equilibrium position. A positively stable body will accelerate toward its equilibrium position when disturbed. 21

Explanation of Aircraft Axes ¢ Axes: l Movement about the lateral axis drawn from the nose to the tail of the aircraft is called roll. l Movement about the longitudinal axis drawn from wingtip to wingtip is called pitch. l Movement about the directional axis extending vertically through the fuselage is called yaw. 22

Lateral Stability ¢ Lateral (roll) stability is determined by: l l dihedral and wing placement. Center of Lift 3°. 75” 23

Longitudinal Stability ¢ Longitudinal (pitch) stability is determined by the longitudinal stability margin (LSM) LSM = 1” Center of Lift Center of Gravity 24

Aerodynamic Requirements: “Positive static stability about all three axes” Directional Stability ¢ Direction (yaw) stability is determined by the directional stability margin (DSM) DSM = 2” Center of Lateral Area Center of Gravity 25

Airframe Design: Wings The Concept of an Airfoil ¢ ¢ ¢ An airfoil profile is a 2 -D shape that is designed to behave in a particular way when subjected to external flow. Airfoils are usually designed to provide a desired coefficient of lift (CL) over a range of angles of attack (Ao. A) and Reynold’s numbers (Re) while minimizing drag. The characteristic behavior of an airfoil is empirically derived from wind tunnel testing and is plotted in Polar Curves. 26

Airframe Design: Wings Airfoil Selection ¢ ¢ ¢ At a cruising speed of 50 mph the UAV must generate a CL of. 202 in order to stay aloft. At a landing speed of 18 mph the UAV must generate a CL of 1. 26 in order to stay aloft. The NACA 3412 airfoil provides these lift coefficients at low drag. 27

Airframe Design: Fuselage Drag Reduction ¢ ¢ Top and side projections of the fuselage are NACA lowdrag bodies. A front projection of the fuselage is an ellipse. 28

General Stability Source: Airplane Flight Dynamics and Automated Flight Controls, Jan Roskam 29

Transfer Functions for all three axes 30 Source: Airplane Flight Dynamics and Automated Flight Controls, Jan Roskam

Overall Aircraft Control Loop 31

Airspeed Control by Elevator Deflection PID Loop 32

Roll Control by Aileron Deflection PI Loop 33

Altitude Step Response Airspeed Control = 20 knots 34

Altitude Step Response Airspeed Control = 25 knots 35

Dynamic Stability Criteria ¢ ¢ ¢ 36 A linear system is stable if and only if the real parts of the roots of the characteristic equation of the system are negative. A linear system is convergent (stable) if the roots of the characteristic equation of the system are real and negative. A linear system is divergent (unstable) if the roots of the characteristic equation of the system are real and positive. A linear system is oscillatory convergent (stable) if the real parts of the roots characteristic equation of the system are negative. A linear system is oscillatory divergent (unstable) if the real parts of the roots of the characteristic equation of the system are positive. A linear system is neutral stable if one of the roots of the characteristic equation of the system is zero or if the real parts of the roots of the characteristic equation of the system is zero.