6 DOF Simulation

8/20/2012

DHRUBA

Summer Training Report |

Dhruba Paul
Dual(Ae & Av)/09/08 A4717209008

6 DOF Simulation

Table of Contents 1. Development of the 6DOF Nonlinear Model in MATLAB/Simulink 1.2Coordinate System 1.2.1 Euler Method 1.3 The 6DOF Equations of Motion System 1.4 Aerodynamics Module

1.5 Thrust Force Module 1.6 Gravity Module

List of Figures Figure 1 Figure 2 Dhruba/Coefficients Figure 2.1 dhruba/Coefficients/Body rate damping Figure 2.1.1. dhruba/Coefficients/Body rate damping/Coefficient w. r. t r Figure2.2 dhruba/Coefficients/Coefficient w r t actuators Figure 2.2.1 dhruba/Coefficients/Coefficient w r t actuators/Rudder/Angle Conversion Figure 2.2.2 dhruba/Coefficients/Coefficient w r t actuators/Coeff w r t elaevator Figure 2.2.3 dhruba/Coefficients/Coefficient w r t actuators/Coefficient w rt aileron Figure 2.3 dhruba/Coefficients/Coefficient wrt to alpha beta Figure: 3 Dhruba/gravity model Figure:4 Dhruba/Flight Parameters Figure: 5 Dhruba/Servo Command

2|Page

6 DOF Simulation

Acknowledgement I express my sincere gratitude to, my mentor and guide, Prof M. S. Prasad for his able guidance, continuous support and cooperation throughout my research, without which the present work would not have been possible. I would also like to thank all my friends for the constant support and help in the successful completion of my project. And above all, to the Almighty God, whose never ceaseing love and for the continued guidance and protection. Thank you to all mentioned above.

SIGNATURE:

(STUDENT.)

3|Page

6 DOF Simulation

1 Development of the 6DOF Nonlinear Model in MATLAB/Simulink

In this chapter, the development of a generic 6DOF nonlinear model is presented. This 6DOF nonlinear model is constructed using the stability derivatives of aircraft. This 6DOF model is built based on the MATLAB/Simulink platform. The structure of the 6DOF model can be broken in several modules as illustrated in Figure 1. Discussions in this chapter are associated with each individual module. Appendix A shows the details of the 6DOF models construction layer by layer.
Figure 1

1.1 Coordinate System

The Earth-Centered Earth-Fixed (ECEF) coordinate system is used in constructing the 6DOF nonlinear model. The origin of the ECEF coordinate system is defined at the center of the Earth. This ECEF frame is rotating about the Earth-Centered Inertial (ECI) reference frame. The representation of the ECEF frame from the ECI frame is simplified by considering only a constant rotation about the Z-axis. The body of interest, the aircraft, is assumed to be rigid and is represented in the ECEF frame. The ECEF coordinate system is implemented in the 6DOF ECEF (Euler) block from the Aerospace Blockset in MATLAB/Simulink as shown in Appendix A.

1.2Aircraft Attitude Representations

In aircraft simulations, there is a need for a mathematical technique to describe the position and orientation of the body of interest in both inertial and non-inertial coordinate systems. Various methods have been developed to determine the aircrafts attitude representation. These include the well-known Euler angles the Euler-Axis rotation parameters, the direction cosines
4|Page

6 DOF Simulation

1.2.1 Euler Method

The Euler method is a well-known approach used for aircraft attitude representations and computations. The Euler method is defined based on the definition of Euler angle representations which can be found in many texts such as References. Using the Euler method, a vector transformation from the Earth-fixed axis to the Body-fixed axis and vice versa can be obtained through the transformation matrix shown in Eq.1 . The notations used in this equation are defined as follows: Sx = Sin(x) and Cx = Cos(x). The subscript b refers to the Body-fixed frames, and the subscript f refers to the Earth-fixed frame.

Equation 1

Since the Euler angles are functions of time, they need to be updated at each time step throughout the simulation process. The Euler kinematic equations are introduced for this purpose. Using the numerical integration method, the equations can be solved to compute the rate of change of the Euler angles in each given time step. The Euler kinematic equations are shhown below.

Equation 2 Note that when the pitch angle () is equal to 90 degrees, the kinematic equations cannot be defined and thus cannot be solved. This is known as the gimbal lock singularity problem.

5|Page

6 DOF Simulation

1.3 The 6DOF Equations of Motion System

The core of the 6DOF nonlinear model are the EOM used in this simulation. In this 6DOF model, the standard six degrees of freedom nonlinear differential equations for a conventional fixed wing aircraft are used. These 6DOF nonlinear differential equations are written as

Forces Equations:
Equation 4

Equation 4

Moment Equations:

Equation 5

From a mathematical viewpoint, Eq.4 and 5 and form a set of six nonlinear differential equations with six unknown variables (u, v, w, p, q, r). Each variable is presented in different equations and interacts with the others. Therefore, the equations cannot be solved individually. The total solution for the system can be obtained only by applying numerical integration to all equations for each given time step.
6|Page

6 DOF Simulation

For the force equations, the inputs for the system are the three major forces applied to the aircraft. These are the aerodynamic forces, thrust forces, and gravity forces. These forces are nonlinear and time variant. The components of each force are broken down in detail and are discussed in Section , Section and Section The inputs for the moment equation are the three moment terms applied to the three aircraft Body-fixed axes. These moments are generated from the aerodynamic and thrust forces with respect to the aircrafts center of gravity (C.G.). As the gravity forces are defined at the C.G. location, no moment is introduced by the gravity forces. The assumption is made that the thrust line passes through the C.G. Therefore, the terms that contribute to the moment equations are only associated with the aerodynamic forces. These aerodynamic moment terms are broken down one-by-one and are presented in Section. The variables P, Q, R from the 6DOF equations are applied to the Euler kinematic equations Eq. 6-8 to update the euler terms. The updated euler terms are then used for vector coordinate transformations. The 6DOF nonlinear equations are implemented in the 6DOF ECEF (Euler) block from the Aerospace Blockset in MATLAB/Simulink as shown in Appendix A.

1.4 Aerodynamics Module

As previously mentioned, the aerodynamics forces are one of the major forces applied to the aircraft, and these forces create aerodynamic moments that contribute to the moment equations. The main purpose of this aerodynamics module is to estimate the values of the aerodynamic forces and moments in the aircraft body frame. A high level block diagram for this module can be simplified as seen in Figure Figure 2

7|Page

6 DOF Simulation

The aerodynamic forces are composed of three forces, which are the lift, drag, and sideforce. These forces generate moments with respect to the center of gravity about the X, Y, Z-axis and are described as the rolling moment, pitching moment, and yawing moment. The component build-up method is used to generate the forces and moments. The total forces and moments
8|Page

6 DOF Simulation

that act on the aircraft are simply the summation of the forces and moments contributed by each component. The aerodynamic forces and moments are first described as dimensionless coefficients, which are associated with the stability and control derivatives given in the modeling results. Equations above are used to implement the aerodynamic force and moment coefficients. Note that when modeling the aircraft in this study, only one set of derivative values at the trim speed. The six aerodynamic force and moment equations are implemented in the subsystem Aerodynamics Coefficients as shown on Appendix . This subsystem is broken down into three different subsystems They are: 1) Coefficient w.r.t. Body Rate, 2) Coefficient w.r.t. Alpha and Beta Angles, and 3) Coefficient w.r.t. Deflection Angles for each Control Surface, which creates the inputs for equations shown below. Starting with the first vector, they are the drag coefficient, sideforce coefficient, lift coefficient, rolling moment coefficient, pitching moment coefficient, and yawing moment coefficient equations, respectively. They are illustrated on pages through in Appendix . The final dimensional values of the forces and moments are then computed from the dimensionless coefficients using the following equations.

9|Page

6 DOF Simulation

The force components shown in Eqs, ,are defined in the stability axis. To make them useable within Eq., which is defined in the Body-fixed axis, a coordinate frame transformation is needed to transfer the stability axis forces to Body-fixed axis forces using the equations expressedbelow.

The implementation of force and moment equations is within the Aerodynamic Forces and Moments block from the Aerospace Blockset in MATLAB/Simulink and is shown on page in Appendix A .

1.5 Thrust Force Module

Thrust force is assumed to be constant at trim condition. The implementation of Thrust force in MATLAB/Simulink and is shown on page in Appendix A .

1.6 Gravity Module

Gravity is the last major force that is applied to the aircraft system. In most simple simulation models, gravity is always assumed to be constant regardless of the aircraft position and altitude with respect to the Earth. The WGS84 Gravity Model block implements the mathematical representation of the geocentric equipotential ellipsoid of the World Geodetic System (WGS84). The block output is the Earth's gravity at a specific location. The WGS84 gravity calculations are based on the assumption of a geocentric equipotential ellipsoid of revolution. Since the gravity potential is assumed to be the same everywhere on the ellipsoid, there must be a specific theoretical gravity potential that can be uniquely determined from the four independent constants defining the ellipsoid. Use of the WGS84 Taylor Series model should be limited to low geodetic heights. It is sufficient near the surface when submicrogal precision is not necessary. At medium and high geodetic heights, it is less accurate. Use of the WGS84 Close Approximation model should be limited to a geodetic height of 20,000.0 m (approximately 65,620.0 feet). Below this height, it gives results with submicrogal precision.

10 | P a g e

6 DOF Simulation

1.7 Atmosphere Module

During the simulation process, updated atmospheric data are frequently required to give feedback to some of the modules in order to provide the necessary computational information. For this reason, an atmosphere module is required to provide the latest atmospheric data for the current altitude. The ISA Atmosphere Model block implements the mathematical representation of the international standard atmosphere values for ambient temperature, pressure, density, and speed of sound for the input geopotential altitude. The use of this ISA Atmosphere model is shown on page in Appendix A.

11 | P a g e

6 DOF Non Linear Model developed using Matlab/Simulink

Scope3 Scope1 Scope

Scope2

25000 V (ft/s) thrust

e

X e (ft) Coef
s tab xyz

F F body Add

(lbf)

Body Euler Angles (rad) DCM

Flight Parameter Bus

bar

q <signal3> CG [3x1] coefficient out Constant V

b

be

PLANT BUS CP M body

V (ft/s)
b

Actuator Coefficients [6; 0; 0] Constant1 Aerodynamic Forces and Moments

<signal5>

Fixed Mass M
xyz

(rad/s) (lbf-ft) d /dt A (ft/s 2 )

b

6DoF (Euler Angles)

Aerodynamic angles PLANT BUS Gravity Force

Total Air Speed

Dynamic Pressure1 gravity model PLANT BUS Mach

Temperature

Pressure

Flight Parameters

APPENDIX A

Pilot Command Pilot command

Actuator

6 DOF Simulation

Pilot command

12 | P a g e

Servo Command

Figure 1

6 DOF Simulation

Dhruba/Coefficients
Figure 2

13 | P a g e

6 DOF Simulation

Body rate damping

Figure 2.1 dhruba/Coefficients/Body rate damping

Product

-KGain

Coefficient w. r. t p

-KProduct1 Gain1

1 Coefficient w.r t q Body Rate Damping Divide

-KProduct2 Gain2

Add

Coefficient w. r. t r 1 Plant Bus

<signal6>

2 Flight Parameters

<signal2>

14 | P a g e

6 DOF Simulation

Coefficient w. r. t r
Figure 2.1.1. dhruba/Coefficients/Body rate damping/Coefficient w. r. t r

Coefficient w.r t q
Figure 2.1.2. dhruba/Coefficients/Body rate damping/Coefficient w.r t q

15 | P a g e

6 DOF Simulation

Coefficient w r t actuators
Figure2.2 dhruba/Coefficients/Coefficient w r t actuators

Angle Conversion
Figure 2.2.1 dhruba/Coefficients/Coefficient w r t actuators/Rudder/Angle Conversion

16 | P a g e

6 DOF Simulation

Coeff w r t elaevator
Figure 2.2.2 dhruba/Coefficients/Coefficient w r t actuators/Coeff w r t elaevator

Coefficient w rt aileron
Figure 2.2.3 dhruba/Coefficients/Coefficient w r t actuators/Coefficient w rt aileron

17 | P a g e

6 DOF Simulation

18 | P a g e

6 DOF Simulation

Coefficient wrt to alpha beta

Figure 2.3 dhruba/Coefficients/Coefficient wrt to alpha beta

19 | P a g e

6 DOF Simulation

Coefficient alpha
Figure 2.3.1 dhruba/Coefficients/Coefficient wrt to alpha beta/Coefficient alpha

Coefficient w r t zero alpha

Figure 2.3.2. dhruba/Coefficients/Coefficient wrt to alpha beta/coefficient w r t zero alpha

20 | P a g e

6 DOF Simulation

Coefficient w.r t alphadot

Figure 2.3.3 dhruba/Coefficients/Coefficient wrt to alpha beta/Coefficient w.r t alphadot

21 | P a g e

6 DOF Simulation

Coefficient w.r t beta

Figure 2.3.4. dhruba/Coefficients/Coefficient wrt to alpha beta/Coefficient w.r t beta

22 | P a g e

6 DOF Simulation

Dhruba/gravity model
Figure: 3

23 | P a g e

6 DOF Simulation

Dhruba/Flight Parameters Figure:4

Dhruba/Servo Command
Figure: 5

24 | P a g e

6 DOF Simulation

Appendix B Result for simulation:

For our simulation we have considered boeing 747. We require stability derivatives ,cg location ,cp location aerodynanmic center , weight, inertia tensor, chord, span ,aspect ratio. We got these values from r.c nelsons flight stability and control. We consider trim condition, where the flight is at an altitude of 20000 feet and a constant thrust 0f 13000 lbs. We are giving pilot command of elevator deflection of 2 degree for 2 seconds and analyse the response of aircraft which is shown by the graphs below.

25 | P a g e

6 DOF Simulation

Future Work
The USAF digital DATCOM program should be used as another modeling method to potentially improve the models used in this research. The 6DOF model can be interfaced with other existing autopilot systems to develop a hardware-in-the-loop simulation platform.

References
1. Center of Remote Sensing of Ice Sheets Research Institute, Lawrence, KS, URL: http://www.cresis.ku.edu [cited 19 September 2008]. 2. Bernstein, L., Bosch, P., et. al., Climate Change 2007: Synthesis Report. Summary for Policy Markers, Intergovernmental Panel on Climate Change [online publication], URL: http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr_spm.pdf [cited 20 August 2008]. 3. Donovan, W. R., The Design of an Uninhabited Air Vehicle for Remote Sensing in the Cryosphere, Masters Thesis, Department of Aerospace Engineering, The University of Kansas, Lawrence, KS, 18 December 2007. 4. Donovan, W. R., CReSIS UAV Critical Design Review: The Meridian, Technical Report 123, Center for Remote Sensing of Ice Sheets, The University of Kansas, Lawrence, KS, 25 June 2007. 5. Underwood, S., Performance and Emission Characteristics of an Aircraft Turbo diesel Engine using JET-A Fuel, Masters Thesis, Department of Aerospace Engineering, The University of Kansas, Lawrence, KS, 05 May 2008.

6.Anderson, J. D., and Wendt, J. F., Computational Fluid Dynamics: An Introduction, 2nd ed., Springer-Velag, New York, 1996. 7. Anderson, J. D., Computational Fluid Dynamics: The Basics with Applications, McGraw Hill Inc., 1995.

8. Phillips, W. F., Mechanics of Flight, John Wiley & Sons, Hoboken, NJ, 2004. 9. Pamadi, B. N., Performance, Stability, Dynamics, and Control of Airplanes, 2nd edition, AIAA Education Series, Reston, VA, 2004. 10. Roskam, J., Aircraft Flight Dynamics and Automatic Flight Controls (Part I), DAR Corporation, Lawrence, KS, 2003.

11,MATLAB/Simulink, Software Package, Version R2012-a, The MathWorks Inc., Natick, MA.

26 | P a g e