Professional Documents
Culture Documents
Simulation of UAV Systems
Simulation of UAV Systems
Symbols
m V & a a b F, Q, Y u, v, w p, q, r CD CL CLF cl cn cm r Ix, Iy, Iz c S b g T Airplane mass Velocity Time derivation of a Angle of attack Angle of drift Euler angles (roll, pitch, yaw) Components of velocity of the airplane mass centre relative to atmosphere Components of the angular velocity of the airplane Drag coefficient Lift coefficient Lift coefficient of the vertical tail surface Rolling moment coefficient Yawing moment coefficient Pitching moment coefficient Air density Moments of inertia about (x, y, z) axes Length of mean aerodynamic chord Wing area Wing span Standard gravitational acceleration Engine thrust
1 Introduction
Civil and military usage of low cost UAVs is becoming more needed. Possibly the most expensive design items are the control and navigation systems. Therefore, one of main questions that each system designer has to face is the selection of appropriate sensors for a specific autopilot system. Such sensors should satisfy the main requirements without contravening their boundaries. Higher sensor quality can lead to a significant rise in costs. In aircraft design this kind of consideration is especially important due to the safety requirements expressed in airworthiness standards. Therefore question is how to determine the optimal solution. This problem is mostly solved by the designers experience and by thorough testing. However, this can be very expensive, and involves with many risks in relation to flight safety. The problem can be resolved by using a suitable simulation method, for example in the Matlab Simulink program environment. This program can be considered as a facility fully competent for this task. An important factor is what is done to manipulate the functions of the program to achieve the autopilot design. A computer only solves logical problems. It cannot implement practical real world entities, and a computer simulation only simulates what in a sense the designer already knows. For the precise design solution it is necessary to have a mathematical model of the aircraft or at least the basic con-
Wing span
Overall length MTOW Max. speed Payload Endurance Engine
Fig. 1: UAV Sojka III
Czech Technical University Publishing House http://ctn.cvut.cz/ap/
109
Power
Pressure regulator Druck DPI 530
RS 232
PC
RS 232 RS 232
HPIB
Fig. 2: Automatic measurement system
straints in its movement. By using suitable simulations it is possible not only to evaluate the sensors, but also to optimise their filters and control algorithms.
2 Experiment
First phase of the project was to verify the sensor parameters declared by the manufacturer. In order to measure the sensor parameters (sensitivity, accuracy, stability, temperature dependence and hysteresis), it was necessary to adapt an existing automatic system by recording data into a file (see Fig. 2). The second phase in this project was to make a statistical evaluation of data obtained by the automatic measurement system. The validity of the measurements themselves was verified by an accuracy analysis of the measurement system and processing the statistical data. The most important quantity in the set of measured data was the pressure variation between two different altitudes, which could be measured very precisely. The entire evaluation of the measured data then helped to find sensor parameters and, consequently, to design a sensor model for the Matlab Simulink program. The designed model was a simplified version, because it reflected only parameters relevant for the specified UAV autopilot design. The sensor delay in this case could be ignored, because its value was negligible in comparison with the previously mentioned sensor parameters. The basic requirement for this project was to obtain data concerning the UAV design system. In this case, the UAV system was described in the following referential axes (Fig. 3) and by a set of differential equations (Equations 1, 2) [3, 4]:
1 rV 2Sc Ix - Iz , rp + cm 2 Iy Iy Iz - Ix 1 rV 2Sb pq + cn . 2 Iz Iz
(2)
All variables were calculated in a frame of differential equations. All aerodynamic parameters needed for these equations were obtained from wind tunnel experiments [4]. The airplane model in Simulink consists of partial blocks. These blocks represent basic mathematical operations or functions and tools that are necessary for modelling. For example, these could be memory, delay, signal sinks and sources (for an illustration, see Fig. 4). The mathematical model of the airplane was compiled as a continuous system, i.e., all calculations were not performed as a time defined sequence but the time interval changed due to the magnitude of the outcomes of the calculations. This solution guaranteed that all relevant quantities would be obtained. This was because the time between two calculations was in the order of milliseconds, which means more than two orders higher frequency than the highest own frequency of airplane movement. The model in Simulink was divided into subblocks representing the individual equations. The main initial conditions were defined in an external .m file, which had to be run before starting the main simulation in Simulink. The rest of the initial conditions and quantities could be set up directly in the Simulink scheme. These conditions were, for example, initial velocity, altitude and flight path. The input quantities to be modelled were: - elevator deflection - aileron deflection
Czech Technical University Publishing House http://ctn.cvut.cz/ap/
110
Fig. 4: Simulink scheme representing the calculation of the first equation from Equation 1
3 Results
The processed data from the automatic measurement system was used for designing a model corresponding to the basic parameters of the sensor. The results showed the significant temperature dependence. This dependence was easy to correct. However, taking into account the desired function of the sensor, this dependence could be ignored. Standard deviation was 0.031 By processing the 152 measured data items, the following results can be obtained:
-
All the calculated output quantities for visualisation could be shown in a graphical interface or saved. These quantities could be, for example, Euler angles and their derivations, translational and angular movements, speed, altitude, etc. Assuming that we had all necessary aerodynamic parameters, the design of a non-linear aircraft model was worked out. The model describes the aircraft behaviour in almost all standard phases of flight. Algorithms representing the autopilot control were also simulated. In order to make the simulation comprehensive, the model was extended by submodels of the wind and the actuators. A simplified diagram of the design simulation is shown in Fig. 5.
80 % of the results were within the interval 0.5m 60 % of the results were within the interval 0.3m 20 % of the results were within the interval 0.1m
Pitch autopilot
Actuators
Data visualisation
Roll autopilot
111
A decision on the suitability of using the sensor as an altimeter could have been made on the basis of the stated manufacturers data, from the supplied datasheets. Unfortunately this method was inapplicable in the case of using the pressure sensor as a vertical speedometer. The definition for calculating vertical speed is obvious from Fig. 6. and Eq. 3. Vertical speed is defined as change of altitude in time. Dh Vertical speed = (3) Dt It is evident from Graph 1. that the steps of discrete time derivation of altitude can cause undesirable step changing in measurement of vertical speed. These steps could be elimi-
Dt
Altitude Dh
Time
Fig. 6: Calculating the vertical speed
Graph 1
Graph 2
112
http://ctn.cvut.cz/ap/
nated by implementing a filter, but this would involve unacceptable time delays. Surprisingly, thanks to the robustness of the system, the system itself was resistant to sharp and large steps of indicated vertical speed. Graph 2 shows the resultant change of altitude. The control mechanism is made by keeping the vertical speed value (in this case 0) and keeping the set altitude value.
who designed the basic version of the automatic measurement system for pressure measurement.
References
[1] VTUL a PVO o. z. (Air Force Research Institute), Czech Republic, http://www.vtul.cz/. [2] http://www.mathworks.com/. [3] Etkin, B., Reid, L. D.: Dynamics of Flight Stability and rd Control, John Wiley & Sons, Inc. 1996, 3 ed. [4] Cech, Z.: Vyuit magnetickho pole Zem ke stabilizaci bezpilotnch prostedk (Usage of the Earths Magnetic Field for Stabilising UAV Flight)(in Czech), dissertation work, Military Academy Antonina Zapotockeho, Brno, CZ, 1989. [5] Proks, M.: Men modelu letounu Sojka-M4 v aerodynamickm tunelu. Technical report Z-3763/02, Aeronautical Research and Test Institute in Prague 2002. [6] http://www.pdiaero.com/downloads/download_files/MIL -F-8785C. pdf, section 3.7 Ing. Petr Kaovsk e-mail: [email protected] Department of Measurements Czech Technical University in Prague Faculty of Electrical Engineering Technick 2 166 27 Prague 6, Czech Republic Dr. Ladislav Smrcek Dr. Colin Goodchild Aerospace Engineering Department James Watt Building University of Glasgow Glasgow, G12 8QQ Scotland, U.K.
4 Conclusion
Simulated flight quantities were evaluated with comparing with real flight records and submitted by real Sojka operators. This correspondence certifies the correctness of the Simulink airplane model. The model of the UAV design system helps to create a powerful tool for a suitability testing of the autopilot and its sensors. This procedure can speed up the choice of sensors which makes a price reduction of their implementation into the UAV autopilot. Particular result from this study is a decision to use a low-cost altimeter for vertical speed measurements. This design method is also suitable for further utilisation in UAV design system simulation. Another possible utilisation of the method is in evaluating of sensor quality in ageing, for the purposes of flight safety.
5 Acknowledgment
The author thanks Mr. Z. Cech, from Delong Instruments Ltd. Brno, Czech Republic, who described the mathematical model of UAV Sojka, and Mr. V. Dvorak from CTU Prague,
Fig. 7
Fig. 8
Czech Technical University Publishing House http://ctn.cvut.cz/ap/
113