Modeling Aircraft Performance and Stability On X-Plane
Modeling Aircraft Performance and Stability On X-Plane
Modeling Aircraft Performance and Stability On X-Plane
Christopher W. S. Thong1
Commercial flight simulator software has been long used by the public to hone their
piloting skills. X-Plane is the first commercially available flight simulator program
which directly outputs mathematical flight data without the need for programming
knowledge. This allows the public to readily simulate flight tests on a computer and also
acquire flight test data. This thesis models the F-15E in X-Plane and compares its
performance and handling characteristics to those predicted by theory and flight
manual data. These comparisons will be used to determine X-Plane’s usefulness as a
teaching and analytical tool in the demonstration of how forces and moments affect an
aircraft’s flight path in flight.
Nomenclature
Terms
CAM - Civil Aeronautics Manual
CFD - Computational Fluid Dynamics
CFT - Conformal Fuel Tanks
CG - Centre of Gravity
CPU - Central Processing Unit
DATCOM - DATa COMpendium
FADEC - Full Authority Digital Engine Control
FAR - Federal Aviation Regulation
FDM - Flight Dynamics Modules
FSX - Microsoft Flight Simulator X
GB - Gigabyte
Ghz - Gigahertz
GUI - Graphic User Interface
HUD - Heads-Up-Display
KIAS - Knots Indicated Air Speed
KTAS - Knots True Air Speed
NACA - National Advisory Committee for Aeronautics
NASA - National Aeronautics and Space Administration
RAM - Random Access Memory
RDS - Raymer Design Software
SFC - Specific Fuel Consumption
TR-824 - NACA Technical Report 824
USAF - United States Air Force
Variables
- Body frontal area
- Coefficient of Drag, consisting of skin friction, interference and pressure drag.
- Average aircraft wing chord
- Airfoil coefficient of lift
- Coefficient of drag on the body
- Maximum aircraft coefficient of lift
- Drag on body
- Wing Element Span
- Frequency
- Lift on wing element
1
Aeronautical Engineering Project, Thesis & Practical Experience ZEIT 4500
1
Final Thesis Report 2010, SEIT, UNSW@ADFA
- Pitching moment
- Dynamic pressure at flight condition
- Reference wing area
- Required Thrust
- Wing element velocity
- Stall speed
- Aircraft weight
-
- Angle of bank
ξ - Damping ratio
- Air density
- Natural frequency
I. Introduction
Computer flight simulation software has for decades been used extensively across the aerospace industry to
aid in the aircraft design. In this thesis, flight simulators will be taken to refer to computer software which
allows the user to pilot an aircraft with a visual representation of the cockpit and surroundings on the computer
monitor. It is also necessary to make the distinction between flight simulators and computer games which make
use of a similar cockpit GUI for the user to interact with. Thus, there is the additional constraint that flight
simulators must have a flight dynamics engine which sufficiently mimics real world physics. This is hard to
quantify so this thesis only considers software to be considered as a flight simulator if it focuses on the flight
dynamics and flying experience. While there are multiple benefits of flight simulation in relation to pilot
training (1 pp. 15-50), there are also numerous benefits for engineering applications. These include but are not
limited to:
time and cost savings as the aircraft can be simulated on the computer to determine its handling
characteristics and performance;
decreased personnel risk due to less actual test flights being flown ; and
the ability to test the effect of certain variables (such as the boundary layer effect or gravitational force)
on the aircraft which may be impossible to test in real life.
However, flight simulators built for engineering applications are usually expensive proprietary software,
built by major aircraft manufacturers to analyze aircraft before production. This thesis will use X-Plane, a
cheap, commercially available flight simulator, to model the F-15E aircraft.
Limitations of computer simulation are mainly based in the fundamental assumptions of the mathematical
model which drives internal calculations. Limitations may be exacerbated by the processing speed of the
computers used and accuracy required. Due to the usefulness of flight simulators, new software is constantly
being developed.
The thesis planning documents are contained in Appendices A, B and C whiles the essential and desirable
objectives of this thesis are listed below. Also, please note that due to the fact that numerous figures presented
in this thesis are in multiple colours, this report is best viewed in colour.
Essential Objectives
Collect dimensions, engine parameters and other required open source information on the F-15E
aircraft which will enable the creation of a model with maximum possible accuracy.
Build a model of the F-15E on X-Plane which can act as a teaching and building tool.
Obtain usable data relating to performance and handling qualities of the F-15E from public sources to
compare to X-Plane outputs.
Predict performance data and handling qualities of the F-15E by theoretical, non flight simulator means
to which the outputs obtained from X-Plane can be compared.
Compare the flight test outputs to predicted performance and stability data to assess the accuracy of the
simulator as well as the F-15E model.
Desirable Objectives
Develop as accurately as possible a working flight simulator model of the F-15E to assess the
possibility of making use of X-Plane to show the effect of changes on an aircraft.
2
Final Thesis Report 2010, SEIT, UNSW@ADFA
II. Methods of Predicting Aircraft Performance and Stability
There is a multitude of non-flight-simulator software which can provide estimates for aerodynamic
parameters, aerodynamic forces as well as stability and control characteristics of aircraft. These estimates will
be compared to certain X-Plane outputs so that the accuracy of using X-Plane can be established. The programs
elaborated on below are the more widely used software, each using different computational methods. It is also
possible to calculate certain parameters directly with the use of methods described by Raymer (2) and Roskam
(3) (4). Both these authors explain the method used to determine performance and stability characteristics and
provide necessary equations, graphs and historical data. They also explain the methodology as well as the
reasoning behind it and equation derivation in their books. However, their approaches are not covered in depth
in this section due to their relative simplicity.
3
Final Thesis Report 2010, SEIT, UNSW@ADFA
stated boundaries. CFD software does that by making use of
discretization methods such as the finite volume method, finite
element method, finite different method and the boundary
element method.
There are numerous ways of determining the fluid flow
within the boundaries. One could look at it as a two
dimensional situation or a three dimensional one. The user can
define a mesh which determines the number and position of grid
points at which the fluid flow must be determined. Figure 1: A Hybrid Grid around an F-15
FLUENT is extremely flexible and capable, with the ability (55)
to model both laminar and turbulent flows, eddy currents and
vortices with a variety of different user defined methods. It can calculate the force, moment and pressure
distributions on objects in a flow based on the test piece’s geometry. FLUENT also presents data in a variety of
pictorial forms, one of which is shown in Fig. 1.
While it is fairly straightforward to develop simple models and meshes for use in FLUENT, it is a much
more complex task to model an aircraft in dynamic flight and was thus impractical to attempt considering the
scope of this thesis. This method can be very accurate in determining the forces on an aircraft but it would be
extremely time consuming to construct the model and complete the simulation. Even making use of a small
section of the aircraft (such as half the vertical tail) would require much calculation and modeling time.
Furthermore, the use of this extremely exact prediction method means that input data such as dimensions and
airfoils used would have to be much more accurate than those obtained for this thesis.
D. Comments
Due to the unavailability of the necessary specifications required to model the F-15E in RDS and the
resource intensity of CFD, these programs were not used to predict the performance of handling qualities of the
aircraft. However, OpenDATCOM appears to be a feasible method of predicting the required outputs. Of
utmost importance is the ability to validate our results. With the above stated software, the computer completes
all calculations and returns the requested outputs. This means that users are unable to obtain intermediate
calculations unless one writes new software to extract this data. Thus, certain calculations as described by
Roskam will be completed to validate a sample of OpenDATCOM outputs. The reason why Roskam’s
textbooks will be used instead of Raymer’s is due to the fact that Roskam is much more comprehensive and
detailed.
A. FlightGear
FlightGear (15) is a free flight simulator maintained by Curt Olson which can be readily hacked to allow
users to tailor the program to their specific need. There are multiple different usable FDMs which make use of
different methods to model aerodynamic forces and moments to predict aircraft flight paths (16).
The most commonly used FDM is JSBSim (16) – an open source FDM in C++ which makes use of the
coefficient buildup method (17 pp. 5,16). JSBSim requires knowledge of the aerodynamic specifications such
as force and moment coefficients as well as stability and control derivatives (18 pp. 18 - 24). JSBSim also
requires data on graphs such as coefficient of lift versus angle of attack graphs. Due to its open source nature,
JSBSim is an extremely well developed FDM as an extensive range of users have contributed ideas for its
continual improvement (18 p. Preface).
Another notable FDM is YASim which, like X-Plane, takes primarily physical dimensions of the aircraft
into consideration for flight modeling (15 p. YASim Readme). The main difference between the two in terms of
data input is that YASim requires data such as cruise and approach information (15 p. YASim Readme) while
X-Plane takes airfoil data (19) instead. Thus, YASim, as with JSBSim, requires one to know to an extent the
abilities of the aircraft being modeled. There is a YASim helicopter FDM which makes use of blade element
theory (20) but throughout the process of literature review, no airplane FDMs were found which make use of
blade element theory.
4
Final Thesis Report 2010, SEIT, UNSW@ADFA
While calculations are performed to determine the aircraft’s flight path in simulation, there are no prebuilt
options in the FlightGear simulation GUI or command line prompts for the user to automatically extract any
data to a spreadsheet or text file (15), (16). Thus, the user would have to write a program to interface with
FlightGear if he wanted to use this software to determine certain data values such as maximum thrust or lift to
drag ratio. One could even take the simulator out of the loop and feed data straight into JSBSim (17) or YASim
and extract data from there. However, these possibilities would require programming knowledge beyond that of
this thesis’ author.
C. Comments
The discussed flight simulation programs have very different backgrounds but similar outcomes – a usable
flight simulator which can be used for pilot training but no built in method of easily extracting flight data from
the simulator. This is largely due to the fact that most commercially available flight simulators make use of
stability derivatives to determine an aircraft’s flight path. Eq. 1 (2 p. 470) shows the yaw stability derivative.
(1)
Not only is stability derivative information kept confidential by aircraft manufacturers and their clients, the
above equation indicates that stability derivatives change with aircraft weight, payload positioning and aircraft
damage. Thus, if one was to simulate an aircraft in one of the above flight simulators, he would usually have to
predict the stability derivative required to achieve certain handling qualities in flight. Thus, it is possible to
have a visual model of a Boeing-747 which actually flies like an F-15 due to the fact that the model’s flight
characteristics are determined by stability derivatives as opposed to the actual simulation model. As a result,
these flight simulators cannot predict aircraft motion like X-Plane does. However, it is possible for one to use
an FDM in FlightGear which does not make use of stability derivatives. However, information on this was
sparse and thus not elaborated on.
5
Final Thesis Report 2010, SEIT, UNSW@ADFA
Blade element theory is the method used by X-Plane to calculate an estimate of the forces acting on an
aircraft in flight. The aircraft’s lifting surfaces are divided into multiple longitudinal strips and the calculated
forces on each strip are then added together. This gives an estimate for the total force due to lifting surface such
as the wings and empennage on the aircraft. As an example, the lift equation shown in Eq. 2 (29) with dr, the
element span, being depicted in Fig. 2. The velocity used is the trigonometric sum of the longitudinal and
transverse component of flow over the wing. This allows X-Plane to consider the effect of sideslip on the forces
on the aircraft. The same method is used to determine the lift, drag and moment on all lifting surfaces made up
of an airfoil.
(2)
B. Plane Maker
Plane Maker (30) has numerous functions which allow the user to create nearly any possible aircraft
geometry. While Plane Maker creates the visual model used for simulation, it must be noted that unlike most
other commercial flight simulators, the dimensions of this model are directly taken into consideration for flight
modeling.
The modules contained within Plane Maker are used for designing the fuselage, wings, empennage,
undercarriage as well as miscellaneous bodies and wings such as wing pylons, speed brakes and CFTs. Due to
the fact that not all wings have wings of steady taper and airfoil, it will be necessary to model the aircraft with
multiple wing sections. Since the F-15E makes use of a NACA64A006.6 airfoil at the wing root and a
NACA64A203 airfoil at the wing tip, it is also necessary to model this in X-Plane with the right airfoil data. X-
Plane can interpolate between wing root and tip to obtain a smooth transition between the two airfoils (31)
which meets the above stated requirement.
In most of the modules, the user simply enters values into boxes for each stated parameter. However, for the
fuselage and miscellaneous bodies, it is possible to load a bitmap image onto the background of the GUI. A
bitmap image of the left side view and top view can be loaded. Plane Maker allows one to edit the fuselage
cross section at a maximum of 20 locations along the fuselage length, with up to nine nodes at each cross
section. This allows for much more accurate modeling of the fuselage than simply stating the radius and length,
as is the case in most software which can only model simple bodies of revolution.
When the designer is required to input the body coefficient of drag, the designer must include interference,
skin friction and pressure drag in this body coefficient of drag value (32). The drag on each body (such as the
fuselage) is then calculated with use of Eq. 3.
(3)
6
Final Thesis Report 2010, SEIT, UNSW@ADFA
Other Plane Maker inputs include engine parameters and locations, forward, default and aft centers of
gravity and coefficient of friction between the wheels and the ground. Plane Maker can also blend parts of the
aircraft together, leading to the ability to smoothly join the fuselage and wings together like on the F-16 aircraft.
If necessary, Plane Maker has a GUI to program artificial stability and autopilot systems into the aircraft.
A notable characteristic of X-Plane is that it is compatible with numerous other well-known software
platforms such as Blender, SketchUp and FSX. Blender and Sketchup are drawing programs which allow one
to create objects. Aircraft models can be loaded from Blender into Plane Maker while scenery objects from all
three of the mentioned programs can be uploaded into X-Plane (33).
C. Airfoil Maker
Airfoil Maker is used to set the characteristics of the airfoils used by an aircraft. General airfoil information
required includes the Reynolds number at which the airfoil data is entered, the thickness to chord ratio and drag
divergence mach number at zero lift. Since a wing’s coefficients of lift, drag and moment vary with Reynolds
number, it is possible for one to insert data for the wing at multiple different Reynolds numbers and X-Plane
will interpolate between the data to determine the most realistic coefficients at the flight Reynolds number. If
data at only one Reynolds
number is input, X-Plane
will approximate for other
Reynolds numbers.
Next, the coefficient to
angle of attack graphs
must be set. The graph
shown in Fig. 3 depicts
the plot for the NACA
64A006.6 airfoil. The
graph shows the
coefficients between an
angle of attack of -200 and
200. Notable points would
include the fact that
minimum drag appears to Figure 3: Coefficient of Lift versus Angle of Attack (57) for NACA 64A006.6. Each
occur at zero angle of graduation on the x and y-axes indicates 10 angle of attack and 0.1 coefficient value
attack, there exists a respectively.
maximum and minimum
lift coefficient which indicates the stall condition and that the
wing has a predominately pitch-down moment, as indicated by
the coefficient which is negative for most angles of attack.
Numerous controls are available which enable the user to
produce these graphs.
Airfoil Maker also presents the coefficient versus angle of
attack graphs from -1800 to 1800 for the users’ information. The
user thus modifies the graphs for the -200 and 200 region and
determines the shape of the graphs over the remaining angles of
attack. The user can then check his airfoil data by calculating
the airfoil lift on drag ratio by inserting the aspect ratio and
Oswald’s efficiency factor of the wing.
There are a number of airfoils already built in X-Plane.
These are mainly standard NACA airfoils. The F-15E uses the
NACA 64A006.6 and 64A203 airfoils (34) and the only NACA
64 airfoil available in Airfoil Maker is the NACA 64-208 airfoil.
Thus, appropriate airfoils need to be created for an accurate
F-15E model. Two sources will be consulted for this – the
NACA Technical Report 824 (TR-824) and Javafoil. TR-824
contains airfoil data for numerous different airfoils and the ones
most similar to the F-15E’s are the NACA 64-006 and NACA
64-206. Data includes mean line data, the coordinates for upper Figure 4: Graphs of Coefficient of Drag
and lower surfaces, predicted critical Mach numbers and against Angle of Attack at multiple
aerodynamic characteristics such as the coefficients of lift, drag Reynolds Numbers for NACA 64A006
and moment. These coefficients are necessary for modeling the (60 p. 177)
airfoil in Airfoil Maker but TR-824 only consists of data from -
7
Final Thesis Report 2010, SEIT, UNSW@ADFA
120 to 160, as shown in Fig. 4. There were even less data points in some other graphs. Javafoil (35), a free
airfoil analysis program was used to approximate the coefficients over the rest of the required -200 to 200 range.
Javafoil outputs are contained in Appendix I.
Javafoil results are highly agreeable (36) with University of Illinois Urbana-Champaign (UIUC) wind tunnel
tests which is indicative of Javafoil’s accuracy. The UIUC experiments’ airfoil data were published in three
volumes (37) (38) (39) between 1995 and 1998 but unfortunately do not contain any data relevant to the airfoils
used in the F-15E (40). TR-824 was completed in 1945 and thus, the accuracy of the UIUC experiments can be
expected to be superior. Furthermore, the X-Plane stock airfoil data was compared to both TR-824 and Javafoil
data and was found to be more convergent to Javafoil data. Thus, when there was a deviation between Javafoil
and TR-824 coefficient values, the Javafoil value was favoured.
8
Final Thesis Report 2010, SEIT, UNSW@ADFA
as described in Appendix D. Cycle dump also gives the resultant radius of gyration in each of the three axes for
the simulated aircraft and the aircraft’s static margin.
F. Limitations
Online forums indicate what users think of X-Plane and how it compares to other commercially available
flight simulator software. However, most forum members are interested in the simulation experience such as
memory usage and the smoothness of graphics. Other than the inability to model fuselage lift, technical
limitations discussed by forum members include the inaccuracy of X-Plane in modeling separated flow over
surfaces as well as the lack of vortex analysis and inaccuracy of transonic and supersonic modeling. Another
real world effect which is not simulated in X-Plane is the wing flex.
Of importance to this thesis is the issue of X-Plane’s transonic and supersonic modeling (26). While
compressible flow effects are considered using Prandtl-Glauert, transonic effects are only simulated by an
empirical Mach-divergent drag increase. This is the Mach number at which compressibility effects start to
become apparent. The only GUI in which the Mach-divergent drag can be input is Airfoil Maker. In
supersonic flight, the aircraft airfoils are assumed to have a diamond shape with the same thickness ratio as the
input airfoils (26). Pressures behind the shock waves are then calculated on each of the plates on the diamond
airfoil and summed to give the total pressures on each airfoil element. This would enable the determination of
lift and drag on each airfoil element. Additional drag on each component due to supersonic effects is calculated
by determining the compression shock on each component. This compression shock is based on the angular
difference between neighbouring components and the local flow.
One interesting limitation of X-Plane is that if display options are set too high for the computer being used,
the cycle rate for completing calculations will decrease and the aircraft will suffer from abnormal oscillations.
Thus, in flight tests, graphics options were set to a minimum.
There are also some technical limitations discussed by online magazines. One such discussed limitation is
the inaccuracy of secondary aileron effect modeling. In reality, if a step aileron input acts on an aircraft, the
aircraft first rolls about its longitudinal axis and yaws in the opposite direction. After the turn is established, the
aircraft should yaw in the direction of the turn and if the turn is continued, the nose should drop and the aircraft
should enter a spiral descent. Sim Pilot Magazine writer looked at how a Cessna 172 in X-Plane reacted to an
aileron input (41). It eventually rolled in the opposite direction due to an increase in airspeed and a pitch up
moment developed after the initial roll. Thus, the spiral dive never developed, emphasizing the differences
between X-Plane’s mathematical modeling and reality. This scenario was recreated in X-Plane a number of
times to confirm the article’s claim and it was validated. This indicates that the authors other claims are also
likely valid. Due to time constraints, no attempt was made to determine the exact deficiency in the X-Plane
modeling.
V. F-15E Overview
In order to accurately model the F-15E in X-Plane, accurate dimensions and engine parameters were
required. The F-15E makes use of the Pratt & Whitney F100-PW-220 or F100-PW-229 engines (42), (42).
Since the only performance charts obtainable were those of an F-15E equipped with two F100-PW-220 engines,
the F-15E was modeled as such. Information on the F-15E was limited to aircraft dimensions, some engine data
and CG locations.
A. Dimensions
Externally, the F-15A and F-15E base airframes are identical, with the largest visible differences being the
addition of CFTs and a tandem cockpit on the F-15E (43 p. 55). The most specific F-15 schematic found is
included in Appendix F and was used to base the F-15E model on. Actual photographs of the F-15E were used
to accurately model the canopy and CFTs.
Table 1: Data related to the F100-PW-220 (47)
B. Performance
Basic engine data can be seen in Table 1.
Performance charts and data included in the F-15E flight
manual includes stall speeds, one engine inoperable
engine rate of climb, all engines operable rate of climb,
fuel required to climb, combat ceiling, optimum long
range cruise, long range cruise, level flight envelope,
maximum speed at level flight, dive recovery, level flight acceleration and sustained level turns. Most charts
are available for a variety of different aircraft configurations.
Usage of X-Plane indicates that the most easily replicated results would be those for stall speeds, level flight
envelope, maximum speed at level flight as well as sustained level turns.
9
Final Thesis Report 2010, SEIT, UNSW@ADFA
An example of the data contained in the flight manual is shown in Appendix G (42 pp. A1-14) which
indicates the stall speed with gear and flaps up or down and provides a myriad of data which can be compared
with X-Plane outputs to test its ability to replicate real lift data.
C. Stability
Gross weight and CG locations for different payload combinations are contained in the F-15E flight manual
(42 pp. A1-8,9, A5-9). The aircraft has a 10 anhedral which has a slight destabilizing effect. This counters the
stabilizing effect of the high wing configuration, thus increasing the F-15E’s maneuverability.
VI. Creation of the X-Plane Model in Plane Maker and Airfoil Maker
In order to create a working model of the F-15E for use in X-Plane, it was necessary to study freely available
flight simulation models of the F-15E to fine-tune the method used to create the model. Unfortunately, all F-15
X-Plane models were found to be based on the model initially developed by Anthony Booher which meant that
there was only one original F-15 model available. This model was deconstructed piece by piece and numerous
errors, inefficiencies and points for improvement were found. One such inaccuracy with Anthony Booher’s F-
15E model is that it did not adhere to the official stall speeds.
Deconstruction and analysis of Anthony Booher’s F-15 model and other X-Plane fighter aircraft models
produced the following framework for producing the F-15E X-Plane Model. Numerous different components
were used to model the complex F-15E shape accurately.
A. Fuselage
The F-15E fuselage was relatively simple to
model. The cross section of the aircraft was
defined at 20 stations along the fuselage length
with the maximum of nine nodes selected at each
station. A bitmap image was loaded into Plane
Maker and each of the nine points can be dragged
into position. Accurate 3-view diagrams of the F-
15A were used to build the fuselage. Fig 6 shows
the top and side view of the model, with the top
view bitmap loaded. The body coefficient of drag
must also be defined in this section. This was
defined to be equal to 0.03, as indicated by
Hoerner’s estimations of drag numerous different Figure 6: The top and side view of the fuselage design GUI
(30), with a top view bitmap loaded to show how the
body shapes (44 pp. 3-9 to 3-22). Interference
definition points line up.
drag was neglected at this stage and would be
taken into consideration if total aircraft drag was lower than expected.
B. Airfoils
Six different airfoils were used on this aircraft, with two being
airfoils provided in the X-Plane package and four being created during
this thesis. The two stock airfoil files used were the flat_plate(very
thin).afl and NACA0009(symmetrical).afl which were used for the
horizontal and vertical tails as well as the connection of the
empennage to the fuselage. The specially created
NACA64A006.6R9.afl, NACA64A203R9.afl and MidptR9.afl airfoil
files were used on the wings and leading edge extension while
eagle_fuse.afl was used within the F-15E fuselage to mimic the F-
15E’s ability to create lift with its fuselage.
The horizontal and vertical tails were modeled as very thin flat
plates as indicated by observation of photographs. This airfoil did not
produce any moment between angle of attack of -200 and 200 and
produced no lift when the angle of attack was 00. In reality, however,
aircraft typically use an airfoil which generates a force opposing wing
lift (2 p. 471). The use of very thin plat plates was deduced to be a
reasonable assumption this is standard practice in the X-Plane
community. Many comparable aircraft provided in the X-Plane
package as well as those produced by the online community make use
of flat_plate(very thin).afl for the horizontal and vertical tails. Figure 7: Side view of the F-15E
Empennage
10
Final Thesis Report 2010, SEIT, UNSW@ADFA
The empennage connection to the fuselage and wings was
modeled by making use of the NACA0009(symmetrical).afl
airfoil file. Being a symmetrical airfoil, it did not produce lift at
00 angle of attack. Also, as seen in Fig. 7 and 8, the airfoil is a
very similar shape to the actual empennage connection.
The NACA64A006.6R9.afl and NACA64A203R9.afl airfoil
files were designed by making use of two sources – TR824 and
the Javafoil program. The relevant TR 824 and Javafoil charts
are included in Appendix H and I respectively which can be
compared to the X-Plane airfoil files shown in Appendix J.
MidptR9.afl was designed by setting each of its input parameters
to be the mean between those of NACA64A006.6R9.afl and
NACA64203R9.afl. This is because MidptR9.afl was used to
ensure that the airfoil at the midpoint between wing root and
wingtip was such that there would be a smooth transition in
airfoil shape from root to tip.
Eagle_fuse.afl was used to simulate the fuselage lift of the F- Figure 8: Bottom view of the F-15E
Empennage
15E. Due to the fact that X-Plane does not reduce the drag or
moment on overlapping components, the coefficient of drag and moment of this airfoil were set to be zero so
that the internal airfoil would not cause additional drag or moments.
11
Final Thesis Report 2010, SEIT, UNSW@ADFA
fuselage will be able to generate lift. This was determined by reducing the coefficient of lift of all airfoils on the
aircraft to zero. The way around this was to insert a dragless wing within the model’s fuselage and X-Plane will
calculate the lift on that wing while still adding the drag of the fuselage to the model.
D. Engines
The engines modeled in X-Plane were F100-PW-220s. Their location and shape was determined by the
three view diagrams. The NASA Glenn EngineSim Program (46) was used to determine additional engine
parameters which were unavailable in the F100-PW-220 information sheet (47). The engines were set to have
specific fuel consumptions of 0.01/hr so as to have minimal effect on aircraft weight as the test flights
progressed. This meant that special test flights would not need to be designed to factor in changing aircraft
weight. Engines are only simulated to have drag if the engines are idling or flamed out.
E. Control Surfaces
The size and shape of ailerons, flaps and the rudders were determined from the three view diagrams. The
entire horizontal tail plane deflected to act as an elevator. The deflection angles were determined by observing
the deflection angles used to create similar X-Plane models. The ailerons were set to deflect 25 0 upwards and
150 downwards, the elevators were set to deflect 300 upwards and 180 downwards while the rudders were set to
deflect 300 in each direction. The flaps were plain flaps and set to deflect 300 downwards at full deflection.
F. Systems
Plane Maker has the ability for a designer to select properties in the aircraft’s electrical, bus, pressurization,
hydraulics, starter, warning, radio altimeter, instrument arc limits, external light, internal light and spring
bouncer systems. Due to the fact that these systems did not alter the X-Plane models flight characteristics and
were not required for accurate flight modeling, this section of modeling was not considered.
G. Artificial Stability
The artificial stability GUI allows one to specify the artificial stability and autopilot characteristics of the
model. The artificial stability system aims to maintain the model’s angle of attack, sideslip and roll-rate.
Without the appropriate artificial stability system, the model suffers from unwanted oscillations. During test
flights to confirm the artificial stability settings, inappropriate settings occasionally caused the model to oscillate
at such a rate that the airframe exceeded its g-limit. Autopilot settings also needed to be entered and fine-tuned.
Again, inappropriate settings caused excessive oscillations which made accurate test flights impossible.
12
Final Thesis Report 2010, SEIT, UNSW@ADFA
K. Payload Table 2: Comparison of Payload Asymmetric Effects
The payloads used on the designed X-Plane
model were stock stores provided with the X-
Plane software. Table 2 provides information
on the stores used and their respective
asymmetric effects. Payload of differing
weights and location were tested to determine
if there was any measurable effect on the
model’s flight characteristics. Actual flight
tests or CFD would be useful to quantify the
effect of the tested payloads.
L. Cockpit Layout
Due to the fact that the cockpit layout of model is insignificant to the flight dynamics of the aircraft, the
cockpit layout of an existing model (48) which already had the essential autopilot switches and HUD was used.
This was to facilitate easier conduct of flight tests.
13
Final Thesis Report 2010, SEIT, UNSW@ADFA
VII. Prediction of F-15E Performance and Stability by Methods other than X-Plane
As stated above, this thesis aims to compare X-Plane flight test results to theoretical predictions as well as
flight manual data. This section will analyse the predictions and determine if they can be used as a baseline to
which flight test results can be compared.
A. OpenDATCOM
Unfortunately, through use of the program, it was discovered that there are some errors present in the
OpenDATCOM GUI. Thus, no outputs were able to be extracted from OpenDATCOM. However, it would be
possible to use the determined inputs in DATCOM+ or the original Digital Datcom program. This would only
be possible if the user has experience in using FORTRAN, which this thesis author does not. Thus, while not
used in comparison of X-plane outputs, DATCOM inputs are attached in Appendix M while explanations to the
terms are contained in (50) so that future work can be readily initiated with regards to the use of DATCOM to
further refine the F-15E model.
Table 3: Percentage Uncertainty in Flight
B. Flight Manual Data
Manual Stall Speed Data
The F-15E flight manual contains data and charts which
have aided in the design of the F-15 X-Plane model. The
charts with which the experimental data are compared to are
contained in Appendix G and N. The charts in Appendix G
indicate the F-15’s stall speed at 00, 150, 300, 450 and 600
angle of bank at a range of aircraft weights between 30000lbs
and 80000lbs while the charts in Appendix N shows the level
flight acceleration on the F-15 at set weights between
42800lbs and 66900lbs. The level flight acceleration graphs
taper off to the maximum speed for the thrust setting
(maximum or military) at 10000ft and 40000ft which is the
value we can compare the F-15 model’s maximum level flight
speed to. Maximum thrust is the aircraft’s maximum thrust
with afterburners while military thrust is the aircraft’s
maximum thrust without afterburners.
There is also an uncertainty associated Table 4: Percentage Uncertainty in Flight Manual Maximum
with the flight manual data. This is due to Thrust Data
the flight manual’s method of data
presentation. The data is all presented on
graphs and it is necessary to read the data
from the graphs. Thus, the uncertainty in the
flight manual data at each point would be
equal to half a division of the graph. Tables
3 and 4 show the percentage uncertainty
associated with each data point taken from
the flight manual. This analysis led to
suspicion that the flight manual data may
contain some anomalous results. This is
because the maximum level flight speed at
61500lbs did not conform to the trend of the
other data points. The maximum level flight
speed for all conditions was expected to decrease with increasing aircraft weight but the data at 61200lbs
appears to be an outlier. The data in the flight manual definitely has an associated uncertainty and this would
most likely be the cause of this discrepancy. Unfortunately, the uncertainty in flight manual data was not
included in the flight manual.
C. Theoretical Calculations
Since DATCOM could not be used to provide the bulk of prediction data, Roskam’s aircraft design manuals
were used to calculate the expected stall speed, maximum level flight speed as well as the phugoid and short
period characteristics.
The stall speed is simply obtained by making use of Eq. 4. The wing reference area was simply obtained
from (34), with air density being dependent on altitude while weight is the chosen variable. The maximum
aircraft coefficient of lift was obtained from flight manual stall speed data by making use of Eq. 4. The chosen
maximum aircraft coefficient of lift was the one which led to the least percentage deviation between
14
Final Thesis Report 2010, SEIT, UNSW@ADFA
theoretically predicted stall speeds and flight manual stall speeds. This method was used as the exact required
data was unobtainable in open source sources.
√ (4)
The theoretical and official clean stall speed values are as stated in Table 5. Due to the lack of open source
information on the maximum aircraft coefficient of lift, this value was determined from the official stall speed.
The differences between the official and theoretical stall speeds for clean and flapped configurations for the full
range of aircraft weights were at their minimums with a clean maximum aircraft lift coefficient of 2.36 and
flapped maximum aircraft lift coefficient of 2.46. Raymer indicates that the flapped maximum aircraft lift
coefficient should be 0.9 more than the unflapped value due to the use of plain flaps (2 p. 326). This is not the
case in these calculations which may indicate that the F-15E’s flaps are not as effective as those on other
aircraft.
The theoretical stall
Table 5: Comparison of Theoretical and Official Stall Speeds
speed in KTAS was
determined from Eq.4
and subsequently
changed to KIAS. The
maximum percentage
difference as seen is -
3.64%. As can be seen,
the theoretical stall
speed is higher than
official values at lower
weights but is lower
than official values at
higher weights. This
would indicate that the
assumption that the
aircraft has a single value for the maximum aircraft coefficient of lift is incorrect. The trend indicates that as the
F-15E’s weight is increase, its maximum aircraft coefficient of lift actually decreases. This is contrary to what
Raymer’s equations suggest in (2 pp. 316 - 327).
The maximum level flight speed can be determined with Raymer (2 p. 518). However, the Raymer
calculations require knowledge of the aircraft’s coefficient of drag. In order to determine the coefficient of drag
theoretically, it is necessary to factor in the skin friction drag, pressure drag, induced drag, interference drag as
well as wave drag. While not impossible to calculate, it will be extremely tedious to consider all required
factors.
Calculated engine power output would not be an appropriate baseline which X-Plane outputs could be
compared to. This is due to the fact that the F100-PW-220s on the F-15E make use of a FADEC and wastegate
to control the power output at different heights for maximum efficiency. Thus, maximum level flight speed
calculations will not be done.
The F-15E flight manual does not contain any numerical data or charts relating to phugoid or short period
motion. Thus, it was necessary to determine these characteristics through theoretical calculations. The
equations used in the calculations are contained in Appendix O. 31 input parameters are required, with 5 of
them either requiring excessive amounts of calculations (such as drag value) or unavailable data (such as )
and 4 parameters being dependent on flight test conditions. Those 5 parameters were thus determined through
flight tests. The rest of the input parameters were
Table 6: Uncertainty in Dynamic Longitudinal simply obtained from measurements of three-view
Stability Theoretical Calculations diagrams or available F-15E data. The 31 input
parameters used in this set of calculations are listed and
valued in Appendix R. The phugoid and short period
and ξ for the F-15E were thus calculated at altitudes of
10000ft and 40000ft and weights of 32000lbs, 55000lbs
and 80000lbs. The reason only six tests were conducted
for the phugoid and short period motion flight tests is
that much time is required for each test as well as its subsequent calculations.
Calculations were completed so as to determine the uncertainty associated with the short period and phugoid
and ξ due to measurements taken. The uncertainty is only due to measurements taken and does not include
15
Final Thesis Report 2010, SEIT, UNSW@ADFA
the uncertainty associated with and . This is because those two parameters were extremely hard to
measure accurately due to the fact that the aircraft suffers from slight oscillations. The oscillations can be seen
in all the X-Plane output graphs contained within the Appendices. The uncertainty associated with and
could not be measured due to this fact. Thus, a certain amount of trial and error was required to determine
and values which were generally correct. Further testing would enable more accurate predictions of
. This is because the is the only prediction which requires those two parameters as inputs.
Thus, since the uncertainty associated with and was not included in the measurements, the uncertainty
for ζ and is actually slightly more than as stated in Table 6.
The uncertainties in Table 6 also do not take into account the fact that Roskam’s equations are
approximations of real world effects. Assumptions are made so approximate an aircraft’s motion to the
equations. This is an avenue for uncertainty but cannot be measured.
One fact which can be deduced from Table 6 is that there is much more uncertainty associated with short
period motion characteristics than phugoid ones.
D. Summary of Results
The data from the official F-15E flight manual was originally expected to be perfect. However, the above
analysis indicated that certain data (such as the maximum speeds at 61500lbs) could be anomalous.
Furthermore, there would definitely be an uncertainty associated with the data which was not stated in the flight
manual. There were no anomalous data points in the stall speed charts. The uncertainty associated with reading
the flight manual charts was between 0.42% and 1.18% for maximum level flight speeds and between 0.56%
and 0.91% for stall speeds. Thus, the F-15E model would be expected to achieve a maximum level flight within
1.2% of flight manual data and stall speeds less than 1.0% different to flight manual data.
Theoretical stall speeds show close compliance with an uncertainty of only 3.7% to flight manual data. This
would indicate that it is unlikely that the flight manual stall speed charts would contain significant errors.
However, analysis indicates that the aircraft coefficient of lift would be a changing value. Due to the lack of
any definitive data on the aircraft coefficient of lift, this parameter will be taken as a constant.
Phugoid and short period and ξ cannot be compared to any real world data. Thus, the theoretical data
cannot be validated against anything. However, they can be expected to be fairly accurate in giving a general
estimate of the actual aircraft and ξ. This is because the and ξ align with expected values. As expected,
is very close to unity, while is closer to zero. The are also significantly
shorter that .
All considered, the predictions are expected to be fairly accurate baselines to which flight test data can be
compared.
A. Stall Speeds
The following method was considered for stall speed testing. It is the method stated by both CAM3 and
FAR Part 23 and summarized in (52 pp. 47-52). The aircraft should first be trimmed at approximately 1.5 times
the stall speed. Following that, the throttle should be brought back to idle and the deceleration rate to stall be
kept at 1 kn/s. The first two conditions could be met easily but there was no easy way of measuring and
controlling the deceleration. Thus, this method was altered slightly such that the test pilot simply keeps the
aircraft flying straight and level till the aircraft stalls.
16
Final Thesis Report 2010, SEIT, UNSW@ADFA
The stall is defined to be the point at which the aircraft pitches nose down uncontrollably in level flight.
This point can be confirmed by using X-Plane’s data output function. The aircraft’s elevator deflection, altitude
and calibrated airspeed can be displayed as a function of time.
The aircraft had no external payload, as stated in the official flight manual. The flight manual contains stall
graphs for the aircraft in the clean condition as well as stall graphs for the aircraft with flaps and gear down.
Thus, flight tests were conducted for the F-15E in both configurations. Appendix L shows the calculations
which indicate that the weight change in each flight test is negligible. While this is completely unrealistic in
mimicking real world physics, it is an advantage of making use of flight simulators and keeps the CG position
as well as the aircraft weight basically constant throughout the test flight.
In order to accurately measure the stall speed of an aircraft at a given altitude and weight, each test point was
tested three times. This is to ensure validity of data. Any anomalous data points would be replaced by data
from another test. The resultant stall speed is the average of the three data points. Due to the fact that official
stall speeds were determined to the closest whole number, flight test results shall also be presented to the closest
whole number.
Two main test groups were conducted – one for the F-15E in the clean configuration and one in the flaps and
gears down configuration. The model’s stall speed was determined at 32000lbs, 45000lbs, 55000lbs, 65000lbs,
75000lbs and 80000lbs as well as at angles of bank of 0 0, 150, 300, 450 and 600 for both F-15 configurations. All
stall tests were conducted at 10000ft.
B. Level Flight Acceleration and Maximum Level Flight Speed Flight Test Procedure
The actual flight test technique for the speed power test of a jet aircraft is quite complex and it involves
many mathematical calculations, including the use of non-dimensionalisation to ensure accurate results. This is
because as fuel is burnt, the aircraft becomes lighter. In a real world flight test, the aircraft would be flown
while varying test altitude in order to keep a constant (52 pp. 99, 100). This would require a large amount of
preflight planning as well as a great amount of pilot skill as it would be imperative to fly exactly the correct
flight path. However, as explained in Appendix L, the aircraft’s change in weight for these flight tests is
negligible. Thus, the method used would be to simply reach the required test altitude, stabilize at the test
altitude, engage the altitude hold autopilot, reach the initial test flight velocity then increase the throttle to its
appropriate maximum or military thrust setting and wait for the aircraft to reach its maximum speed.
The maximum level flight speed at maximum thrust was taken to be when the throttle was at maximum
thrust with afterburners engaged, altitude is constant and the aircraft was unable to accelerate anymore. The
maximum level flight speed at military thrust was the same but without afterburners.
The aircraft payload and weight were as stated for the stall speeds. The test flights were conducted with a
clean airframe. Due to the continuous nature of this flight test, only one test needs to be conducted at each
altitude and weight.
The tests were carried out at 10000ft and 40000ft, at 9 different aircraft weights between 42800lbs and
66900lbs. This is such that each available flight manual graph can be compared against. The tests at 10000ft
were conducted both at maximum thrust as well as at military thrust while those at 40000ft were conducted at
maximum thrust. This is due to limitations in flight manual data.
17
Final Thesis Report 2010, SEIT, UNSW@ADFA
Figure 13: Chart used to determine the
Phugoid Damping Ratio (52 p. 250)
(5)
√
Where (6)
The phugoid motion is a large amplitude variation of aircraft airspeed, pitch angle and altitude. From a side-
on perspective, the aircraft would appear to be travelling in a sinusoidal motion. The important values
determined by this flight test are and the . The of an aircraft is it’s resonant frequency, the
frequency at which oscillations will increase amplitude with increased oscillations. ξ is related to the ratio of an
oscillatory amplitude to the following amplitude. If the ξ is positive, it is a measure of how quickly an
oscillation converges. If the ξ is negative, it is a measure of how quickly an oscillation diverges.
The aircraft payload was as stated for the stall speeds. Due to the substantial time taken for each test and
their consequent required calculations, only three weights were considered for the phugoid motion tests. These
weights were 32000lbs, 50000lbs and 80000lbs. CG location and weight changes were not considered as
indicated by Appendix L.
Two flight tests were conducted at each flight condition so as to obtain results from both possible flight test
methods. The aircraft was allowed to undergo at least five full oscillations in each flight test so that an average
and can be obtained.
The flight tests were conducted at altitudes of 10000ft and 40000ft, meaning a total of twelve different flight
tests.
18
Final Thesis Report 2010, SEIT, UNSW@ADFA
The verification section will indicate the accuracy of X-Plane’s flight characteristic prediction. Verification
will be done in two ways - one will be to simply compare the obtained flight test results to the official and
theoretical values while the other will be to take into
account the change in tested flight characteristic with
the flight test variable. The second form of
verification is called the delta method as it simply
considers the delta, or change, of a variable. This
would enable any inaccuracies in modeling the F-15E
to be neglected. The verification section will also
determine the source of inconsistencies between flight
test, theoretical and official values and look at
possible ways to improve the F-15E model if it is
determined that further improvements are required.
It should be noted that in this section, any
reference to flight testing would be specifically to
flight tests conducted in X-Plane with the F-15E
model designed by this thesis author. Figure 15: Depiction of Lift and Weight Forces on a
Banked Aircraft (58)
A. Clean and Flapped Stall Speed Test Analysis
As an aircraft’s weight increases, so does its stall
speed, as indicated in Eq. 4. The stall speed also
increases with aircraft angle of bank. This is because
roll decreases the lift force working against aircraft
weight, as depicted in Fig. 15. These trends are visibly
present in the Fig. 16 and Fig. 17 plot. As angle of
bank and aircraft weight increases, so does stall speed.
It can also be seen that at two points in each graph, the
plots for stall speeds at 00 and 150 angle of bank
overlap. This can be attributed to pilot error.
Furthermore, the official clean stall speeds for
unbanked flight is114kts at 32000lbs and 174kts at Figure 16: Flight Test Clean Stall Speeds
75000lbs while those for 150 banked flight is 116kts at
32000lbs and 177kts at 75000lbs. These values are
extremely close and a slight deviation in data may lead
to the plots overlapping. Similarly, the flapped stall
speed values are very close together for 00 and 150
angle of bank. Thus, apart from the infrequent
overlapping, the trend behaves as expected. This would
indicate that the X-Plane software simulates stall speeds
correctly. The following section will quantitatively
consider the ability of X-Plane’s simulation software to
accurately model stall.
B. Clean and Flapped Stall Speed Verification Figure 17: Flight Test Flapped Stall Speeds
19
Final Thesis Report 2010, SEIT, UNSW@ADFA
X-Plane flapped stall speed test results and official values for the entire range of aircraft weight and angle of
bank is less than 3.4%, as indicated by Table 8.
The clean and flapped stall speeds indicate both a high degree of accuracy of the X-Plane F-15E model in
terms of its stall characteristics as well as a high degree of accuracy of the X-Plane software in modeling the
stall condition.
21
Final Thesis Report 2010, SEIT, UNSW@ADFA
designer only inputs one set of drag coefficient data without distinction between subsonic and supersonic
speeds.
Another difference is that the drag suddenly increases at Mach 1 instead of increasing from the critical Mach
number to a peak at around Mach 1. The drag coefficient obtained from flight tests then remains at nearly the
same level before further increasing and peaking at about Mach 1.1. This indicates that the mathematical
modeling involved in mimicking transonic and supersonic drag has significant room for improvement. The
imperfect transonic and supersonic drag modeling may be the reason for the F-15E aircraft to be unable to reach
its appropriate maximum velocities.
As indicated by the above analysis, the drag coefficient is accurately modeled at subsonic speeds. This
indicates that the drag coefficient of individual components is modeled correctly. At transonic and supersonic
speeds however, total drag coefficient is excessively high, more than twice the value of expected results. This
likely indicates that drag coefficient is overestimated by X-Plane. However, excessive acceleration in the
subsonic zone indicates that the X-Plane F-15E model has too much thrust at both tested altitudes.
Unfortunately, due to the lack of definitive thrust at 100% low compressor speed and the critical altitude of the
engines, accurate results would only be obtainable by trial and error. The correct thrust and drag settings will
improve and the method will be described in a later section.
22
Final Thesis Report 2010, SEIT, UNSW@ADFA
subsequent trough or vice versa. In order to increase accuracy, a minimum of five full oscillations were
conducted for each flight test. Multiple ratios of peak amplitude to trough amplitude were thus obtained and the
average value was determined for subsequent use. The peak to trough ratio was thus a fluctuating value
throughout the experiments. This led to a slight uncertainty in the obtained . Furthermore, the method used for
presenting the correlation between peak to trough ratio and was extremely coarse as shown in Fig 13. More
accurate graphs or data could be found so it would be expected that small changes in would not be detectable.
23
Final Thesis Report 2010, SEIT, UNSW@ADFA
be noted that this analysis does not indicate accuracy of the calculated theoretical predictions, merely the
possibility that they are correct based on these references.
1. Determination of the Optimum Aircraft Altitude and Weight for X-Plane Payload Flight Tests
Experiments were conducted to determine the normal spread of results between
Table 12: Percentage
identical flight tests. Nine flight tests were performed in each flight condition and the
Uncertainty in
standard deviation of the test results in each flight condition was determined as shown Standard Deviation
in Appendix U.
For a normal distribution, approximately 68% of measurements will lie within one
standard deviation of their mean (54 p. 34) and about 95% of measurements within two
standard deviations. Thus, the smaller the standard deviation, the lower the spread of
test results and the more accurate a comparison between different flight conditions.
Table 12 indicates the number of test flights required for a certain percentage
uncertainty in standard deviation. It was decided to conduct nine test flights such that a
low percentage uncertainty would be obtained while minimizing the amount of time
spent conducting test flights.
Once the optimum flight condition has been determined, the standard deviation will
be determined. This standard deviation would be such that there is a 95% probability that the flight parameter
would be within two standard deviations of the population mean and the uncertainty associated with this
standard deviation would be 25%.
The F-15E model was tested in the following conditions:
24
Final Thesis Report 2010, SEIT, UNSW@ADFA
Nine flights were conducted for the F-15E model in each of the above five flight conditions, with conditions
[1] and [4] being used as controls. The flight tests were conducted to determine the aircraft’s thrust, lift, drag,
side force, lift to drag ratio, total coefficient of lift, total coefficient of drag as well as the aileron, elevator and
rudder deflections. These parameters will be the only ones discussed as other parameters are either dependent
on the above parameters or are not meaningful to this analysis.
The standard deviation of the analyzed parameters for the above flight conditions were compared with the
control flights. A sample of the calculated standard deviations for the above set of flight tests is included in
Appendix U. Due to the large amount of data involved in each flight test, only a sample of standard deviations
for each flight condition is included in the Appendix. The data in Appendix U indicates that the standard
deviation of parameters was inversely proportionate to both altitude and aircraft weight. Thus, it would be
advantageous for the aircraft to have a lower altitude and weight. However, if the aircraft were flying too low,
X-Plane would take into consideration wing in ground effect which would increase the lift on the aircraft (19) ,
leading to errors in results. Ground effect affects aircraft when the aircraft is at a distance from the ground less
than half its wing span (13 p. 354). Thus, while excessive, the aircraft test altitude was chosen to be at 10000ft
such that it would be out of ground effect and would be at the same altitude as several other tests.
2. Determination of the Optimum Aircraft Speed Range for X-Plane Payload Flight Tests
The graphs of the above stated parameters against Table 13: Comparison of the Oscillation Amplitudes
true airspeed in flight conditions [1] and [4] are for Various Speed Ranges
contained in Appendix V and W respectively. These
graphs indicate that the F-15E model will suffer from
oscillations in lift, drag and side force which increase
with velocity. This trend is steady for the F-15E model
in both the clean and asymmetrical payload condition.
These oscillations are reflected in the graphs of aileron,
elevator and rudder deflections, also contained in
Appendices V and W. However, the control surface
deflection graphs for the clean aircraft indicate
abnormal aileron deflections between velocities of
260KTAS and 300KTAS as well as abnormal elevator
and rudder deflections between speeds of 340KTAS
and 370KTAS. The control surface deflection graphs
for the F-15E under payload shows similar trends.
Also, for the F-15E in a clean condition, the rudder
deflection graphs indicate that the rudder undergoes slight oscillations at velocities below 360KTAS.
As can be seen in Table 13, the speed range with the lowest oscillation amplitudes for the clean and
asymmetrical payload condition would be between 200KTAS and 250KTAS.
Thus, considering all the above factors, the optimum test speed range was determined to be 250KTAS to
200KTAS. This would minimize unwanted oscillations in aircraft forces and control surfaces, abnormal control
surface deflections as well as ensure the aircraft is well out of the transonic and supersonic zone.
As shown in Appendix V and W, the oscillations in aircraft lift follow a horizontal trend line. These
oscillation amplitudes can thus be determined easily. The lift was not exactly equal to weight because the
aircraft pitch varies with velocity and the resultant force which acts against aircraft weight would be a
trigonometric sum of aircraft thrust and total wing lift. The average total wing lift for all the clean flight tests at
10000ft was thus approximately 32030lbs, while that of the 2 x AIM-54 asymmetric payload tests was
34040lbs. Both these total wing lift values are slightly higher than the aircraft weight in those conditions. This
is because of two factors. Firstly, the actual force acting against gravity is a trigonometric sum of aircraft lift
and thrust, due to the fact that the aircraft is not flying at exactly 00 pitch. Secondly, the aircraft is undergoing
slight oscillations in pitch which causes the wing lift, total lift and resultant lift forces to vary.
25
Final Thesis Report 2010, SEIT, UNSW@ADFA
Table 14: Calculations for the Percentage Uncertainties for the Stated Flight Parameters.
parameters. Thus, other payloads would have standard deviations lower than those determined for the
aforementioned two aircraft configurations. The reason why the 2 x AIM-54 payload is considered an extreme
is not just that it’s weight exceeds that of a other payloads but also because it provides the greatest lateral
asymmetry (in ft-lbs).
One output parameter which needs explanation is the thrust. This is because the thrust is a negative value.
This is because the minimum possible aircraft weight was chosen – 32000lbs. This is only possible with nil
fuel. Thus, the engines are not operating so are producing drag, not thrust. This leads to the negative thrust
values. This is how X-Plane calculates drag on aircraft engines – it is factored into the thrust calculations.
2. Conclusion
The combination of all the above calculations indicates that the above parameter results from a single test
would be within 2.2% or population mean 95% of the time, with an associated uncertainty of 25%.
This leads to numerous possible deductions about the mathematical modeling involved in X-Plane. One
possible deduction is that there are very few possibilities or random errors occurring which may distort
experimental results. The likely reason for the uncertainty in the repeatability flight tests is that the autopilot
does not act in exactly the same manner in each flight test, leading to slight differences. Another conclusion
which can be made is that flight tests a very repeatable as long as the flight path and initial conditions are
readily controllable.
1. Expected Trends
It is expected that as lateral asymmetry increases, so
would aileron and rudder deflection. This would be
because increased lateral asymmetry would cause the
aircraft’s autopilot to increase aileron and rudder
deflection to as to maintain straight and level flight.
It is also expected that lift and drag would increase
with payload. This is because increased payload means
that total lift would have to increase to maintain straight
and level flight. However, total lift would be a Figure 23: Indication of Lateral Asymmetry Effects
trigonometric addition of wing and horizontal stabilizer
lift and thrust. Since the aircraft is in straight and level
flight, drag would be acting purely horizontally.
Fortunately, this total lift value is the value provided in
the X-Plane data output file so the total lift can be easily
charted.
26
Final Thesis Report 2010, SEIT, UNSW@ADFA
different from the rest in terms of rudder deflection. However, the 2xAIM-7 payload is the one with the largest
rudder deflection instead of the 2xAIM-54 payload. This is unexpected because the 2xAIM-54 payload has a
larger lateral asymmetry than the 2xAIM-7 payload. The reason for the aircrafts deviation from the expected
trend is unknown and retests indicate similar results. This could be due to a postulated complex relationship
between lift, drag, pitch and control surface deflection and would be able to be pinpointed by use of further
testing.
The side force on the F-15E did not show any sort of discernible trend. This is contrary to what was
expected. However, since the rudder deflection is so minimal, it is possible that the effect of payload on side
force cannot be detected.
Lift and drag showed mixed trends. As expected,
aircraft lift increases with payload weight. This can
be seen in Fig. 25. Aircraft drag also increases with
payload weight and size. The 2xAIM-54
configuration clearly has more drag than the 1xAIM-
54 configuration. This is to be expected as X-Plane
does drag calculations on the payloads.
As velocity increases, pitch decreases and thus
drag decreases due to the decrease in frontal area and
induced drag. However, when landing gear is
extended, aircraft drag increases with velocity. This Figure 25: Indication of Lateral Asymmetry Effects
is because when landing gear is extended, the
increase of drag on the landing gear far outweighs
the decrease in drag due to decrease pitch and is
visually represented in Fig. 26.
The GUI for creating payloads has a body
coefficient of lift input. As with other bodies, the
designer must include interference, skin friction and
pressure drag in this value before the total drag force
is calculated by Eq. 3. This is the only way X-Plane
determines the interference effect of the payload on
the lift over the wing directly above it.
It can be seen that X-Plane does simulate Figure 26: Indication of Lateral Asymmetry Effects
changes in the F-15E aircraft rather well in terms of
trending correctly. Unfortunately, quantitative analysis could not be completed for these flight tests in this
thesis due to the lack of data to which results could be compared.
D. Conclusion
The flight tests conducted for this section of the thesis have shed light on numerous aspects of X-Plane.
With the right flight conditions and number of test flights, the uncertainty in a result can be reduced to a
minimal value. The best flight condition to conduct flight tests in was determined to be a simple level flight
deceleration test between the velocities of 250KTAS and 200KTAS.
The desired aim of this thesis was to produce an F-15E X-Plane model which could show the effect of
changes to an aircraft. This change was modeled as increasing lateral asymmetry. The above analysis proves
that lateral asymmetry on the aircraft model does indeed produce expected trends. The effect of an alteration to
the F-15E airframe can also be tested. This would likely have to be conducted by adding a miscellaneous body
to the F-15E aircraft, taking note that an accurate coefficient of drag must be known beforehand. Only then can
accurate results be obtained.
One important point to note is that the X-Plane model suffers from slight oscillations at all flight conditions.
This is because of the time required to complete one cycle of calculations. Once one cycle is completed, the
flight parameters would have changed slightly, leading to the need for the autopilot to recorrect the aircraft’s
flight path. Even if a supercomputer were used, oscillations would still exist. The only difference would be that
the magnitude of oscillations would be tiny in comparison. The oscillations can be seen in the flight output data
charts contained in this thesis paper as well as in the appendices. In order to improve the precision of output
data, this issue can be addressed with the use of a more powerful computer. The computer used for these flight
tests has an Intel Core i3 CPU processor, clocked at 3.07Ghz with 3.87 GB of usable RAM.
27
Final Thesis Report 2010, SEIT, UNSW@ADFA
XI. Final Comments
The above described flight tests gave valuable insight as to the accuracy of the F-15E model as well as the
ability of the X-Plane software to correctly mimic real world physics. Any inaccuracies with regards to the
tested flight characteristics have been determined. These are summarized below and the results of further
improvements have been explained too. However, there are still many possibilities to improve on the designed
F-15E model. Numerous other flight tests can also be conducted so as to analyse the model’s characteristics in
other modes of flight, such as the spiral spin or sustained level flight turn.
D. Future Work
Further work can be conducted on the F-15E model to further improve the aircraft’s accuracy. This can be
facilitated by the conduct of additional flight tests which would indicate the necessary edits required on the F-
15E model. Additional flight tests to determine the F-15E’s takeoff distances, rotation speeds, single engine
28
Final Thesis Report 2010, SEIT, UNSW@ADFA
rates of climb, climb performance, descent performance and sustained level flight turn performance would be
ideal. While the flight manual contains data on other flight characteristics, the above stated characteristics are
the easiest to compare flight test results to.
Flight characteristics not mentioned in the flight manual can also be tested with the use of theoretical
predictions. The aircraft’s static longitudinal stability, dynamic lateral-directional stability and maneuvering
stability are a few characteristics which can be predicted by the use of (4). The directions for the correct
conduct of all the above mentioned flight tests are contained in (52).
Use of DATCOM will provide another avenue to which stability characteristics can be compared.
Knowledge of FORTRAN would enable one to use the DATCOM source code. Since slight alterations need to
be made to the aircraft before DATCOM can be used, such as the use of a single vertical tail as opposed to two,
an aircraft should be designed in X-Plane making use of the exact dimensions used in DATCOM. This would
ensure that the DATCOM outputs are valid for the F-15E.
With enough resources, CFD analysis of the F-15E could be conducted for more accurate determination of
forces and moments on the aircraft. This would provide a benchmark to which the payload comparison
deceleration test results can be compared such that those flight test results can be verified. The reason why it is
preferable to use X-Plane over CFD to determine the forces, moments and other characteristics of the F-15E in
flight is that X-Plane takes a much shorter time to do the essential calculations. CFD can be used to ensure that
the input parameters are correct and X-Plane can then be used to do the bulk of the calculations. However, if
CFD were to be used, it would be imperative to obtain much more accurate measurements of the F-15E. This
would be so as to achieve the best possible correlation between the X-Plane and CFD models and the real
aircraft.
Further study into possible flight conditions for precise measurements should be taken. This would enable
the F-15E model to be suffering from the least possible amount of oscillations. This would increase the
precision of data collected from X-Plane.
E. Conclusion
A large amount of research has been conducted thus far to ensure accuracy and validity of results. Future
work should focus on improvement of the F-15E model as well as verification of its abilities. The best possible
way to verify the F-15E model would be to make use of either CFD or actual F-15E flight test data. Even if the
X-Plane aircraft model is not extremely accurate, the delta method can be used to determine the effect of
changes to the aircraft. However, X-Plane has proven in this thesis to be a rather accurate tool in predicting the
flight characteristics of an aircraft if its input parameters are correct. This is seen in the very close correlation of
obtained flight test stall speeds to flight manual data, as well as the correlation of maximum level flight
acceleration and speeds to flight manual data once the aircraft was improved. This would lead to X-Plane being
a powerful prediction tool in determining the approximate result of a change to an aircraft, with more precise
analysis methods being used to focus on certain scenarios which may be been determined through use of X-
Plane.
Acknowledgements
This thesis would not have been possible without the help of my thesis supervisor – Andrew Neely.
Throughout this thesis, he has been a great help on this thesis and a source of many fantastic ideas for
improvement. I would also like to thank Anthony Booher, who not only provided me with privately designed
aircraft models which are unavailable online but also aided me in understanding X-Plane and Javafoil. Without
him, my aims would have been much harder to meet. Lastly, I would like to thank Austin Meyer. One would
not expect the creator of such a well selling piece of software to have the time to reply an endless stream of
emails about the mathematics behind X-Plane but Austin has proved to be an exception to this rule. To all
others who have helped me along the way – thank you.
29
Final Thesis Report 2010, SEIT, UNSW@ADFA
References
1. Williams, B. Microsoft Flight Simulator as a Training Aid. Newcastle, Western Australia : Aviation Supplies
& Academics, Inc., 2006.
2. Raymer, D. P. Aircraft Design: A Conceptual Approach (3rd Edition). Reston, Virginia : AIAA, 1999.
3. Roskam, J. Airplane Design Part VI: Preliminary Calculation of Aerodynamic, thrust and Power
Characteristics. Lawrence, Kansas : DAR Corporation, 2000.
4. —. Airplane Design Part VII: Determination of Stability, Control and Performance Characteristics; FAR and
Military Requirements. Lawrence, Kansas : DAR Corporation, 2000.
5. Hoak, D.E. USAF Stability and Control DATCOM, Volumes 1 to 4. Irvine, CA : Global Engineering
Documents, 1978.
6. Williams, J. E., Vukelich, S. R. The USAF Stability and Control Digital DATCOM. St. Louis, Missouri :
USAF, 1979. AFFDL-TR-79-3032.
7. Holy Cows Inc. DATCOM+. [Online] Mar 17, 2010. [Cited: Mar 28, 2010.]
http://www.holycows.net/datcom/.
8. Hoak, D.E. USAF Stability and Control DATCOM. Dayton, Ohio : Flight Control Division, Air Force Flight
Dynamics Laboratory, 1978.
9. Holy Cows Inc. DATCOM+, Software Package, Ver 2.6c. Holy Cows, Inc. [Online] Mar 17, 2010. [Cited:
Mar 20, 2010.] http://www.holycows.net/datcom/.
10. Matyushev, A. OpenDATCOM, Software Package, Ver 1.1. openAE. [Online] Jan 5, 2010. [Cited: Feb 8,
2010.] http://www.openae.org/software.
11. OpenDatcom. OpenDatcom. Open Aerospace Software Commmunity. [Online] January 05, 2010. [Cited:
March 25, 2010.] http://openae.org/software.
12. Raymer, D. P. RDS-Student Software. Reston, Virgina : AIAA, 2006.
13. —. RDS-Student Manual. Reston, Virginia : AIAA, 2006.
14. Wendt, J. F. Computational Fluid Dynamics: An Introduction (3rd Edition). Berlin, Germany : Springer,
2008.
15. Olson, C. L. FlightGear Ver 2.0, Software Package. FlightGear. [Online] 2010. [Cited: Mar 18, 2010.]
http://www.flightgear.org/.
16. —. The FlightGear Manual. FlightGear. [Online] Feb 25, 2010. [Cited: Mar 15, 2010.]
http://www.flightgear.org/Docs/getstart/getstart.html.
17. Berndt, J. S. JSBSim Reference Manual. JSBSim. [Online] Mar 20, 2010. [Cited: Mar 28, 2010.]
http://jsbsim.sourceforge.net/JSBSimReferenceManual.pdf.
18. —. JSBSim: An Open Source Flight Dynamics Model in C++. Providence, Rhode Island : AIAA Modeling
and Simulation Technologies Conference and Exhibit, Aug 16-19, 2004. Vols. AIAA 2004-4923, pp. 1-27.
19. Laminar Research. X-Plane 9 Operation Manual. s.l. : Graphsim Entertainment, 2008.
20. HHS. YASim helicopter model?? FlightGear Forum. [Online] Feb 26, 2008. [Cited: Apr 22, 2010.]
http://www.flightgear.org/forums/viewtopic.php?f=4&t=1049.
21. Away, Fly. Boost your Flight Simulator Experience with FSX Aircraft Add-ons. Fly Away. [Online] April
14, 2010. [Cited: May 2010, 13.] http://flyawaysimulation.com/article3983.html.
22. fsdownload. Flight Simulator X Downloads Directory for FSX FS2004 FS2002 FS2000 FS98 FS95.
www.fsdownload.com. [Online] Feb 23, 2008. [Cited: June 17, 2010.]
http://www.fsdownload.com/modules.php?name=News&file=article&sid=891.
23. Beckwith, J. Mudpond. Mudpond. [Online] 2010. [Cited: Apr 8, 2010.] http://www.mudpond.org/.
24. Sors, H. FS Airfile Editors and Design. Flight Dynamics and Navdata. [Online] 2010. [Cited: Apr 8, 2010.]
http://pagesperso-orange.fr/hsors/fsairfile.html.
25. Microsoft Flight Simulato. About the Aces Team. Microsoft Flight Simulator. [Online] Jan 26, 2009.
[Cited: Apr 10, 2010.]
http://www.microsoft.com/Products/Games/FSInsider/news/Pages/AMessageFromAces.aspx.
26. Laminar Research. Appendix A: How X-Plane Works. X-Plane Desktop Manual. [Online] 2010. [Cited:
Mar 20, 2010.] http://wiki.x-plane.com/Appendix_A:_How_X-Plane_Works.
27. —. X-Plane. [Online] 2010. [Cited: Apr 29, 2010.] http://www.x-plane.com/index_pro.html.
28. Taureau, N. Home of the X-Plane Flight Simulator Community. X-Plane.org. [Online] 2010. [Cited: Apr
10, 2010.] http://forums.x-plane.org/.
29. Ross, D. Structural Analysis of a Cessna 182 Wing During Manoeuvres using X-Plane Flight Simulator.
Canberra, Australia : UNSW@ADFA, 2008.
30. Laminar Research. Plane Maker, Software, v9.50rc.
31. —. Supplement: Airfoil-Maker. X-Plane Desktop Manual. [Online] 2010. [Cited: Apr 1, 2010.]
http://wiki.x-plane.com/Supplement:_Airfoil-Maker.
32. Meyer, A., HansVanKampen. Plane-Maker Menus. X-Plane On-Line Instruction Manual. [Online] 5 2002.
[Cited: 5 20, 2010.] http://roust.gotdns.com/instructions/Manual_Files/Plane-Maker%20Menus.html.
30
Final Thesis Report 2010, SEIT, UNSW@ADFA
33. Harris, J. X-Plane Scenery Tools. X-Plane Scenery. [Online] 2009. [Cited: Mar 14, 2010.]
http://marginal.org.uk/x-planescenery/tools.html.
34. Jane's all the World's Aircraft. Boeing F-15E Eagle. Jane's. [Online] Feb 26, 2010. [Cited: Mar 20, 2010.]
http://www.janes.com, search for "F-15E".
35. Hepperle, M. Javafoil, Software Package, Ver 2.13. Javafoil. [Online] Mar 1, 2007. [Cited: Apr 12, 2010.]
http://www.mh-aerotools.de/airfoils/jf_applet.htm.
36. —. Comparison between Theory and Experiment. Javafoil. [Online] Oct 15, 2007. [Cited: Aug 11, 2010.]
http://www.mh-aerotools.de/airfoils/methods.htm#ComparisonTheoryExperiment.
37. Selig, M.S., Guglielmo, J.J., Broeren, A.P., Giguere, P. Summary of Low-Speed Airfoil Data, Vol. 1.
Virginia Beach, VA : SoarTech Publications, 1995.
38. Selig, M.S., Lyon, C.A., Giguère, P., Ninham, C.N., Guglielmo, J.J. Summary of Low-Speed Airfoil Data,
Vol. 2. Virginia Beach, VA : SoarTech Publications, 1996.
39. Lyon, C.A., Broeren, A.P., Giguère, P., Gopalarathnam, A., Selig, M.S. Summary of Low-Speed Airfoil
Data, Vol. 3. Virginia Beach, VA : SoarTech Publications, 1998.
40. UIUC LSATs Airfoils Tested. UIUC Applied Aerodynamics Group. [Online] 2010. [Cited: August 11,
2010.] http://www.ae.illinois.edu/m-selig/uiuc_lsat_airfoilsTested.html.
41. Stock, C. Flight Dynamics: FSX and X-Plane Battle it out. Sim Pilot Magazine. [Online] 2010. [Cited: Apr
5, 2010.] http://www.simpilotnet.com/index.php?option=com_content&task=view&id=20&Itemid=9.
42. McDonnell Aircraft. USAF Series F-15E Aircraft Flight Manual. U.S.A. : McDonnell Aircraft, 1991. TO
1F-15E-1.
43. Jenkins, D. R. McDonnell Douglas F-15 Eagle. North Branch MN USA : Specialty Press, 1997.
44. Hoerner, S. F. Fluid-Dynamic Drag. Bakersfield, CA, USA : Hoerner Fluid Dynamics, 1992.
45. Cohen, E. Israel's Best Defense. New York : Orion Books, 1993.
46. National Aeronautics and Space Administration. EngineSim Version 1.7a. NASA. [Online] Oct 25, 2006.
[Cited: July 18, 2010.] http://www.grc.nasa.gov/WWW/K-12/airplane/ngnsim.html.
47. Jane's Aero-Engines. Pratt & Whitney F100. Jane's. [Online] Jul 10, 2009. [Cited: Mar 27, 2010.]
http://www.janes.com, search for "F100-PW-220".
48. Austin, D. , Booher,A. , Leger B. F-15E Strike Eagle, v9.1. X-Plane.org. [Online] Feb 3, 2009. [Cited: Feb
14, 2010.] http://forums.x-plane.org/index.php?autocom=downloads&showfile=7090.
49. Meyer, A. , Creator of X-Plane. Electronic Mail Correspondence. 2010.
50. Williams, J.E., Vukelich, S.R. The USAF Stability and Control Digital Datcom. Ohio : United States Air
Force, 1979. AFFDL-TR-79-3032, Vol I.
51. Smith, H. C. Understanding Performance flight Testing (2nd Edition). Hightstown, New Jersey : McGraw-
Hill, 2002.
52. Kimberlin, R. D. Flight Testing of Fixed Wing Aircraft. Reston, Virginia : AIAA, 2003.
53. Specification, U.S. Military. MIL-F-8785C. s.l. : U.S. Military, 1980.
54. Mendenhall, W. , Sincich, T. Statistics for Engineering and the Sciences. Upper Saddle River, New
Jersey : Pearson Prentice Hall, 2007.
55. Pointwise .Inc. F-15E Spin Analysis. Pointwise .Inc. [Online] 2010. [Cited: Mar 28, 2010.]
http://www.pointwise.com/apps/f15e.shtml.
56. Jane's Aero-Engines. General Electric F110. Jane's. [Online] May 29, 2009. [Cited: Mar 27, 2010.]
http://www.janes.com, use the "main search" to search for "F110-GE-129".
57. Laminar Research. Airfoil Maker, Software Package, Ver 9.50rc3. 2010.
58. Flickr.com. Flickr.com. [Online] [Cited: July 02, 2010.]
http://farm4.static.flickr.com/3428/3800915982_92a7db0bb0.jpg.
59. SRA D. M. CULLEN, USAF. DefenseImagery.Mil. [Online] Apr 30, 2001. [Cited: Feb 20, 2010.]
http://www.dodmedia.osd.mil/DVIC_View/Still_Details.cfm?SDAN=DFSD0405951&JPGPath=/Assets/Still/2
004/Air_Force/DF-SD-04-05951.JPG.
60. Abbot, I. H. , von Doenhoff,A. E. , Stivers Jr,L. S. Report No. 824 - Summary of Airfoil Data. Washington
D.C. : NACA, 1945.
31
Final Thesis Report 2010, SEIT, UNSW@ADFA