Flight Performance Uav
Flight Performance Uav
Flight Performance Uav
type of mission profile. This program can also be used to confirm the point performance
which is calculated by conventional methods and show the true ranges and endurances of
each segment. This method is also used for graphical representation of the actual mission
with a detailed breakdown of each segment. Unlike point performance, the user can check
the performance of the aircraft at any point during the flight and also fly it with different
speeds and altitude within one mission and check the feasibility of the mission. The example
of UAV in this paper is a fixed wing aircraft with the ROTAX 914UL turbofan engine. The
airfoil of the wing is NACA4412, and the typical aero data is generated for this airfoil. The
program includes the engine data, the propeller data, the aircraft mass data and the aero
data in particular. The propeller data is a complex function of velocity of aircraft and power
provided by the engine. The engine data depends on the propeller efficiency, height and
velocity of the aircraft. The aero data also depends on the velocity of the aircraft, thus the
whole program is a complex interdependent simulation and every parameter is carefully
modeled so that singularity in the simulation does not occur. The simulation is designed not
only for the particular UAV but also can be generalized for other UAVs as well by changing
the aircraft data. The features of the simulation include the user friendly technique to find
the performance of any aircraft by just changing the data. It helps in making unlimited
segments of the flight of the UAV to check whether a mission is possible before going on the
mission. It shows that the method in this paper is advantageous in flight performance
analysis in detail for UAV.
Nomenclature
DOF = Degree of freedom
UAV = Unmanned air vehicle
MATLAB = Matrix Laboratory
RPM = Revolutions per minute
Cf = Aerodynamic force coefficient
γ = Climb angle
•
ψ = Turn rate
ϕ = Bank angle
GUI = Graphical User Interface
η prop = Propeller efficiency
1
Post Graduate Student of Flight Vehicle Design, School of Aeronautics, Northwestern Polytechnical University,
Xi’an, P.R. China, AIAA Student Member.
2
Associate Professor in Manufacture Engineering of Aeronautical and Aerospace Engineering, School of
Aeronautics, Northwestern Polytechnical University, Xi’an, P.R. China.
1
American Institute of Aeronautics and Astronautics
Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.
η E −loss = Engine power loss
η Ins−loss = Propeller Installation loss
Rm = Radius of earth
Fengine = Thrust provided by the engine
D = Drag acting on the aircraft
L = Lift acting on the aircraft
Pengine = Engine power delivered
F = Force acting on aircraft
W = Weight of aircraft
m = Mass of aircraft
g = Gravity
a = Acceleration
h = Altitude of aircraft
S = Distance moved by aircraft
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
I. Introduction
T HIS work is an effort to develop a simulation to capture the flight performance of a fixed wing UAV including
mission profiling. The UAV model is a fully developed 6DOF program which captures the performance of the
UAV during all phases of flight. The mathematical model of the UAV is a complex program of interdependent
variables computed at every time step during the simulation. The goal of this exercise is to develop a fast and simple
performance tool which can show the flight performance of a UAV in much detail compared to the point
performance technique.
In this paper, a method is explained to calculate the performance of an aircraft with the help of analytical and
experimental data of the aircraft. The UAV has 6DOF and is fully functional in six-dimensions. In addition to
modeling the flight model of the UAV, special attention is given to fly the UAV with desired performance and
constantly checking the failure conditions of the UAV e.g. stalling and running out of fuel.
All the data of UAV including propeller performance, engine performance, structure, atmosphere and aero data
are fully utilized. An effort is made to couple all these data variables which are interdependent on each other to
make a flight program which can show the performance of the UAV in true sense. The whole flight simulation
program is simulated in a close loop. Careful attention is paid in the modeling of the aircraft so that singularity in the
closed loop does not occur. Commercial mathematical software MATLAB/SimuLink is used to model the trajectory
and natural flight of the aircraft. Simulation modeling of mathematical software is used to complete the process.
The flight dynamics parameters, angle of attack and flight path are also extracted at every small interval. These
intermediate properties are also saved for determining the performance evaluation of the UAV. The stalling of the
aircraft is checked at every instant and in case of stall, the mission of the aircraft is said to be failed. Manual mission
profiling is also available for the user so that the feasibility of a mission can be checked before actually going on the
mission. This technique is also helpful in evaluating the actual boundaries of the system at different intervals within
one mission.
The whole process is handled automatically and easily with the proposed technique as compared to the
convenctional methods of checking the performance of the UAV with the help of flight envelope or point
performance. Thus by using this technique, we can easily find out the borderline cases of the UAV performance
which in point performance may be hidden. This yields to understanding the UAV and its performance more
comprehensively and gives a more graphical understanding of the performance of UAV compared with conventional
methods.
The overall advantage of this technique is to develop a performance tool which can show the performance of
UAV in more detail compared to conventional methods and help in mission profiling of the UAV before the actual
flight. Also it can help the user in understanding the performance of the UAV much more easily with the help of
graphical representation of the performance of UAV.
2
American Institute of Aeronautics and Astronautics
II. Previous Approaches
Every new aircraft design starts with the preliminary system requirements including range, endurance and ceiling
of the aircraft. After the initial design is met, the system goes into detail designing which leads to detailed
aerodynamic analysis, structural analysis and control design. Performance analysis of UAV is done a lot in the past.
A good example of the performance calculations is done using point performance for teaching purposes1. Other
works on the flight performance include the 6DOF flight program of a helicopter2, testing and evaluation of aircraft
flight performance3 and sensitivity analysis of UAV4. Besides, the related work on flight performance of aircraft
includes the work by Sweeten5 and Ro6, but both are either limited to point performance or do not show the mission
profiling of the aircraft. This work will give additional information in terms of flight performance of the aircraft and
also the mission profiling with the ease in the hands of the user. A lot of work is found in the literature about the
preliminary designing of aircrafts7 and controller design of aircraft8. Modelling and simulation of UAVs are also
done in the past9, but the approach to the problem is entirely different. The approach is only limited to transition
phases, not including the whole flight of the aircraft.
Close work relating to this technique was done on the Meridian UAV by calculating its performance using
analytical and experimental data5. Flight performance of the aircraft is also calculated in this approach but is limited
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
to point performance of the aircraft, the data gathering technique in the above mentioned paper is similar to what is
done in this research. The data of ROTAX 914UL is taken in the same fashion as is taken in the above research on
Centurion 1.7. Meridian UAV also has sensitivity analysis4 done for its enhancement of flight performance. All of
this work is still unable to show the mission profiling of the aircraft.
Many aerodynamic and modelling equations for this paper are taken from the work done by Zhang10 and
Raymer11. The modelling equations are used to build a user friendly, changeable and improvable performance model
of UAV in case the model has to be improved.
A lot of help in simulation building and understanding of the MATLAB simulation tool can be seen in the work
done by Turevskiy12. Reference to this will also be used later in this paper. Simulink modelling of the aircraft in this
reference is a very good and effective way of modelling a new developing aircraft or UAV. However, extensive use
of SimuLink tools renders the performance program developer to make changes in the system required by developer
as the SimuLink packages are fixed. So the best way for the development of the simulation program for UAV is to
use the SimuLink tools which does not affect the system performance and which are useful, and make the rest of the
program by using the general equations of motions. This yields to changeable and better simulation program,
however this technique is more prone to human error but the advantage of mixing the Simulink blocks with user
created blocks is quite good.
For better understanding of the modelling of performance of the UAV, other work done on the UAVs was also
studied closely13. However, they did not affect the outcome of this research so their details are not mentioned.
1 0.15
Cd
CL
0.8 N-S
0.1 AVL
0.6
0.4
0.05
0.2
0 0
-5 0 5 10 15 20 -5 0 5 10 15 20
Alpha Alpha
Figure 2. Lift coefficient vs Angle of attack (deg). Figure 3. Drag coefficient vs Angle of attack (deg).
4
American Institute of Aeronautics and Astronautics
C. Engine performance and specifications
ROTAX 914UL engine was selected for this aircraft. The engine is shown in the Fig. 4 and the basic details of
the engine are given in Table 3. The engine loss is taken as 5% of the total Power.
The performance of the ROTAX 914 is taken from the operators manual of ROTAX 914UL15 and the effect of
altitude on power variance is ROTAX brochure14. The effect of SFC reduction with altitude for ROTAX 914UL was
not found specifically, a generalized SFC reduction was taken from the work done by Schulz16. The performance
curves used for performance evaluation of the system are given in the figures below. Fig. 5 shows the power change
with throttling and also the effect of altitude on Rotax 914UL power, Fig. 6 shows the SFC of the engine and also
the effect of altitude on SFC.
Engine RPM vs Max. Continous Power at different Altitude
80
Altitude = 0 km
Altitude = 2 km
70 Altitude = 4 km
Altitude = 6 km
Max. Continous Power (kW)
Altitude = 8 km
60
50
40
30
20
3000 3500 4000 4500 5000 5500 5800
Engine RPM
Figure 5. Engine RPM and Max. Continous Power
5
American Institute of Aeronautics and Astronautics
Engine RPM vs Fuel Consumption with different Altitude
28
Altitude = 0 km
26 Altitude = 2 km
Altitude = 4 km
Altitude = 6 km
24
Altitude = 8 km
20
18
16
14
12
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
75
70
65
60
55
50
45
40
35
50 100 150 200 250 300
Velocity of UAV (km/hr)
Eq. (7)
V
•
F = ma, m = ∫ m,V = ∫ a, S = ∫ V Eq. (8)
Pengine / SFC = fn(h, throttle ),η prop = fn(V , throttle ), gravity = fn(h ) Eq. (9)
B. Simulation Modeling
For simulating the performance of the UAV, MATLAB coupled with SimuLink is used. MATLAB is basically
used for storing the data and for variable definitions. MATLAB is also used for mission simulation which is done by
calling simulink’s simulations. Also the mission definition is done in MATLAB.
The inputs and outputs of the system are discussed in Table 4. The inputs of the system are defined by the user
before the start of simulation.
Inputs Outputs
Dead mass Range of each Segment / Total Range
Fuel mass Endurance of each Segment / Total Endurance
Lift & Drag coefficient wrt alpha Required angle of attack
Engine Power wrt height & RPM Throttle requirement for the aircraft
Fuel flow rate wrt RPM & Velocity Mass consumed during each segment
Propeller efficiency wrt RPM & Velocity 3D Mission profile
Initial Bank angle (For loiter only) Turn radius (For loiter only)
Climbing and Descending angle of attack Local flight path angle
Gear Ratio of engine Propeller RPM
Density of fuel Turn rate (For loiter only)
Table 4. Inputs and Outputs of the program.
Simulink is used for modeling of equations for UAV. SimuLink has advantage of graphical representation of
equations and also it gives more ease for the programmer to develop the programs in time domain which is the case
in the performance evaluation of UAV. Four different simulation programs are made for the ease of mission
definitions which include climb, cruise, loiter and descent. Different missions can be made with these programs by
calling these simulations again. Similar example of modeling the equations of motion of aircraft in Simulink is
found in the work done by MathWorks12. However, only few equations are selected in this work from the simulation
set for aerospace vehicle from SimuLink. Rest of the equations are modeled mathematically. Fig. 8 shows how
SimuLink in-build blocks were used for ease, Fig. 9 shows the user made blocks and Fig. 10 and Fig. 11 shows the
whole simulation blocks.
7
American Institute of Aeronautics and Astronautics
Figure 8. Atmospheric modeling with help of aerospace blockset.
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
V. Simulation Results
A user-made mission was defined using the simulation technique and mathematical model discussed above. The
mission definition and the details of the mission are given in the Fig. 12. The mission complete summary is given in
Table. 4. The mission is distributed in nine segments with segments in the following manner.
Climb→Cruise→Descent→Cruise before loiter→Loiter→Cruise after loiter→Climb→Cruise→Descent
9
American Institute of Aeronautics and Astronautics
Cruise Cruise
Cruise Cruise
Parameter Climb Descent before Loiter after Climb Descent
Out Back
loiter loiter
Total Time, hr 0.41 2.45 2.62 2.77 3.77 3.92 4.12 6.0 6.82
Segment
0.41 2.04 0.17 0.15 1.0 0.15 0.2 1.88 0.82
Duration, hr
Total Range,
54 504 532 562 562 592 622.5 1072.5 1184
km
Segment Range
54 450 28 30 562 30 30.5 450 111.5
, km
Initial Mass, kg 540 533 503.5 501 499 485.5 483.5 480 451.25
Final Mass, kg 533 503.5 501 499 485.5 483.5 480 451.25 439
Fuel
Consumption, 7 24.5 2.5 2 13.5 2 3.5 28.75 12.3
kg
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
Initial Height
0.01 7.0 7.0 5.0 5.0 5.0 5.0 8.0 8.0
(km)
Final Height
7.0 7.0 5.0 5.0 5.0 5.0 8.0 8.0 0.01
(km)
Cruise Velocity
- 220 - 200 200 200 - 240 -
(km/hr)
Total Fuel
101
Consumed, kg
Remaining Fuel 9 (8.9%)
Table 5. Mission summary of fixed wing UAV.
The simulation results of the fixed wing UAV are given in the figure plots below. In simulation, the blue color
denotes the climb, the red color shows cruise, the magenta color denotes descent and the loitering is shown by green
color. Fig. 13 shows the variation of velocity in the whole mission of the aircraft, the black line shows the stall
velocity and is for the confirmation of the mission that it did not stall during the course of flight. Fig. 14 shows the
three dimensional mission of the UAV. Fig. 15 shows the time vs height plot of the aircraft. Fig. 16 shows the thrust
variation with time. Fig. 17 shows the angle of attack variation of the aircraft. Fig. 18 shows the mass variation and
Fig. 19 shows the percentage throttle, Fig. 20 shows the variation of local flight path angle (pitch) with time. Fig. 21
shows the loitering pattern of the aircraft. Fig. 22 shows the variation of propeller efficiency and Fig. 23 shows the
propeller RPM variation with time.
Velocity vs Time plot
250
3-Dimensional mission profile
200
8
Velocity (km/hr)
6
Height (km)
150
4
2
100
0
1.5
Actual Velocity 1500
1
Stall Velocity 1000
50 0.5 500
0 1 2 3 4 5 6 7
Time(hr) Side Range (km) 0 0 Range (km)
Figure 13. Velocity vs Time plot of mission. Figure 14. Three dimensional plot of mission.
10
American Institute of Aeronautics and Astronautics
Height vs Time plot Thrust vs Time plot
8 1.6
7 1.4
6
1.2
5
Height (km)
Thrust (kN)
1
4
0.8
3
0.6
2
1 0.4
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
0 0.2
0 1 2 3 4 5 6 7 0 1 2 3 4 5 6 7
Time (hr) Time (hr)
Figure 15. Height vs Time plot of mission. Figure 16. Thrust vs time plot of mission.
Variation of Angle of Attack Mass change wrt Time
7 540
6 520
Angle of Attack (deg)
5
500
Mass (kg)
4
480
3
460
2
1 440
0 420
0 1 2 3 4 5 6 7 0 1 2 3 4 5 6 7
Time (hr) Time (hr)
Figure 17. Variation of angle of attack wrt Time. Figure 18. Mass variation of system wrt Time.
Throttling of UAV Flight path angle wrt Time
80 15
70
Local Flight Path Angle (deg)
10
Percentage Throttle (%)
60
5
50
0
40
-5
30
20 -10
0 1 2 3 4 5 6 7 0 1 2 3 4 5 6 7
Time (hr) Time (hr)
Figure 19. Throttle variation wrt Time. Figure 20. Flight path angle vs Time plot.
11
American Institute of Aeronautics and Astronautics
Loitering pattern of UAV Propeller efficiency wrt Time
1.4 85
1.2 80
0.8 70
0.6 65
0.4 60
0.2 55
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
0 50
561 561.5 562 562.5 0 1 2 3 4 5 6 7
Range (km) Time (hr)
Figure 21. Loiter pattern of UAV. Figure 22. Propeller efficiency wrt Time plot.
2100
2000
Propeller RPM (rpm)
1900
1800
1700
1600
1500
1400
0 1 2 3 4 5 6 7
Time (hr)
Figure 23. Propeller RPM variation wrt Time.
References
[1] Baughn, J. and Wolf, J., "Teaching Aircraft Flight Performance in Aeronautical Engineering Using PC-Based Flight
Simulation," AIAA Journal. SAE Technical Paper .1998.
[2] Subodh, B., and Richard, C., “Six DoF Dynamic Modeling and Flight Testing of a UAV Helicopter,” AIAA Modeling and
Simulation Technologies Conference and Exhibit. 15-18 Aug 2005.
[3] Burdun, I. and Mavris, D., "A Technique for Testing and Evaluation of Aircraft Flight Performance During Early Design
Phases," SAE Technical Paper 975541, 1997.
[4] David, A. R., Shahriar, K., Benjamin, C. S., and Valerie, J., “Modeling and Sensitivity Analysis of the Meridian Unmanned
Aircraft,” AIAA Infotech@Aerospace 2010 Conference. 20-22 Apr 2010.
[5] Benjamin, C. S., Dave, R., and Shahriar , K., “The Meridian UAV Flight Performance Analysis using Analytical and
Downloaded by PURDUE UNIVERSITY on January 20, 2015 | http://arc.aiaa.org | DOI: 10.2514/6.2012-2595
13
American Institute of Aeronautics and Astronautics