AUTONOMOUS FLIGHT CONTROL AND GUIDANCE SYSTEM OF ACCIDENT AIRCRAFT

Similar documents
Trajectory Optimization for Safe, Clean and Quiet Flight

A Study on Berth Maneuvering Using Ship Handling Simulator

Time Benefits of Free-Flight for a Commercial Aircraft

GUIDELINES FOR FLIGHT TIME MANAGEMENT AND SUSTAINABLE AIRCRAFT SEQUENCING

Validation of Integrated Safety-enhanced Intelligent flight control. Yoko Watanabe ONERA/DTIS-Toulouse

Mechanics of Frisbee Throwing

Advanced Flight Control System Failure States Airworthiness Requirements and Verification

Research on Fault Tolerant Controls within GARTEUR

AIRBUS FlyByWire How it really works

MPC FAULT-TOLERANT FLIGHT CONTROL CASE STUDY: FLIGHT Jan M. Maciejowski Colin N. Jones,1

Towards Autoomous ISR by a Team of Coopera;ng Gliders

Impact of Landing Fee Policy on Airlines Service Decisions, Financial Performance and Airport Congestion

Efficiency and Automation

ScienceDirect. Prediction of Commercial Aircraft Price using the COC & Aircraft Design Factors

Simulation of disturbances and modelling of expected train passenger delays

AiMT. Flight Trajectory Modelling to Increase General Aviation Safety. Advances in Military Technology Vol. 6, No. 1, June 2011

Research Article Study on Fleet Assignment Problem Model and Algorithm

Decision aid methodologies in transportation

An Analysis of Dynamic Actions on the Big Long River

ANALYSIS OF AIR TRAFFIC EFFICIENCY USING DYNAMIC PROGRAMMING TRAJECTORY OPTIMIZATION

Design Exploration of Low-Speed Wind Tunnel for Improved Flow Characteristics with STAR-CCM+ and Optimate+

PRAJWAL KHADGI Department of Industrial and Systems Engineering Northern Illinois University DeKalb, Illinois, USA

Optimization Model Integrated Flight Schedule and Maintenance Plans

Operational Evaluation of a Flight-deck Software Application

TWENTY-SECOND MEETING OF THE ASIA/PACIFIC AIR NAVIGATION PLANNING AND IMPLEMENTATION REGIONAL GROUP (APANPIRG/22)

Turboprop Propulsion System Malfunction Recog i n titi ion on an d R d Response

CFIT-Procedure Design Considerations. Use of VNAV on Conventional. Non-Precision Approach Procedures

HIGH SPEED FLIGHT DEMONSTRATION PROJECT

Fuel Conservation Reserve Fuel Optimization

Transportation Timetabling

Flight Trajectory Planning for Fixed Wing Aircraft in Loss Of Thrust Emergencies

UC Berkeley Working Papers

STEM FUTURES. Air Travel STEM. Works! Using Maths Tasks. About the Industry. About Your Task

CHAPTER 5 SIMULATION MODEL TO DETERMINE FREQUENCY OF A SINGLE BUS ROUTE WITH SINGLE AND MULTIPLE HEADWAYS

Research on Flight Operational Efficiency for Fuel and Noise

Evaluating the Robustness and Feasibility of Integer Programming and Dynamic Programming in Aircraft Sequencing Optimization

MULTIDISCIPLINARYMEETING REGARDING GLOBAL TRACKING

A Review of Airport Runway Scheduling

NOISE ABATEMENT PROCEDURES

Airport Simulation Technology in Airport Planning, Design and Operating Management

Special edition paper Development of a Crew Schedule Data Transfer System

Estimating the Risk of a New Launch Vehicle Using Historical Design Element Data

Gain-Scheduled Control of Blade Loads in a Wind Turbine-Generator System by Individual Blade Pitch Manipulation

PBN ROUTE SPACING AND CNS REQUIREMENTS (Presented by Secretariat)

REPUBLIC OF KENYA MINISTRY OF TRANSPORT AND INFRASTRUCTURE DEPARTMENT OF AIR ACCIDENT INVESTIGATION

Unmanned Aircraft System Loss of Link Procedure Evaluation Methodology

Index. Springer International Publishing AG 2018 I. Schagaev, B.R. Kirk, Active System Control, DOI /

FUZZY LOGIC BASED FLIGHT CONTROL SYSTEM FOR THE NIGERIAN AIR FORCE AIR BEETLE AIRCRAFT

FLIGHT PATH FOR THE FUTURE OF MOBILITY

Quiz 2 - Solution. Problem #1 (50 points) CEE 5614 Fall Date Due: Wednesday November 20, 2013 Instructor: Trani

Study on the assessment method for results of ship maneuvering training with the simulator

Unmanned Aircraft Systems Integration

LESSON PLAN Introduction (3 minutes)

Lessons Learnt From The EUROCONTROL Wake Impact Severity Assessment Flight Simulator Campaign

Optimized Profile Descents A.K.A. CDA A New Concept RTCA Airspace Working Group

Do Not Write Below Question Maximum Possible Points Score Total Points = 100

FILE AND FLY PROCEDURES AND TECHNIQUES FOR INTEGRATION OF UAVS IN CONTROLLED AIRSPACE

Air Traffic Complexity: An Input-Output Approach. Amy R Pritchett, Keumjin Lee and Eric JM Feron School of Aerospace Engineering Georgia Tech

Single Line Tethered Glider

EFFICIENT CONTROL OF ARRIVAL TIME AT A CONGESTED AIRPORT'S TERMINAL AREA

Safety Syllabus. VFR into IMC

THE IMPACTS OF AIRCRAFT INCIDENT ON THE UNIT OPERATING COSTS OF CIVIL AIRCRAFT

STUDY OVERVIEW MASTER PLAN GOALS AND OBJECTIVES

ANALYSIS OF THE CONTRIUBTION OF FLIGHTPLAN ROUTE SELECTION ON ENROUTE DELAYS USING RAMS

Congestion. Vikrant Vaze Prof. Cynthia Barnhart. Department of Civil and Environmental Engineering Massachusetts Institute of Technology

HOW TO IMPROVE HIGH-FREQUENCY BUS SERVICE RELIABILITY THROUGH SCHEDULING

Research on Pilots Development Planning

Bachelor of Science in Aviation Option in Aviation Management

SIMULATION MODELING AND ANALYSIS OF A NEW INTERNATIONAL TERMINAL

VISUALIZATION OF AIRSPACE COMPLEXITY BASED ON AIR TRAFFIC CONTROL DIFFICULTY

ATM Seminar 2015 OPTIMIZING INTEGRATED ARRIVAL, DEPARTURE AND SURFACE OPERATIONS UNDER UNCERTAINTY. Wednesday, June 24 nd 2015

All-Weather Operations Training Programme

Glass Cockpits in General Aviation Aircraft. Consequences for training and simulators. Fred Abbink

CASCADE OPERATIONAL FOCUS GROUP (OFG)

A RECURSION EVENT-DRIVEN MODEL TO SOLVE THE SINGLE AIRPORT GROUND-HOLDING PROBLEM

April 16, Erik Larson

EXPERIMENTAL ANALYSIS OF THE INTEGRATION OF MIXED SURVEILLANCE FREQUENCY INTO OCEANIC ATC OPERATIONS

H-II TRANSFER VEHICLE TRAJECTORY PLANNING AND FLIGHT OPERATION RESULTS

Airport Gate Assignment A Hybrid Model and Implementation

A Methodology for Integrated Conceptual Design of Aircraft Configuration and Operation to Reduce Environmental Impact

OPTIMAL PUSHBACK TIME WITH EXISTING UNCERTAINTIES AT BUSY AIRPORT

Technical Standard Order

Integrated Optimization of Arrival, Departure, and Surface Operations

Interval Management A Brief Overview of the Concept, Benefits, and Spacing Algorithms

USE OF 3D GIS IN ANALYSIS OF AIRSPACE OBSTRUCTIONS

Workshop. SESAR 2020 Concept. A Brief View of the Business Trajectory

Aircraft Arrival Sequencing: Creating order from disorder

Adaptive Trajectory Planning for Flight Management Systems

Instrument Proficiency Check Flight Record

Cross-sectional time-series analysis of airspace capacity in Europe

HEAD-UP DISPLAY (HUD), EQUIVALENT DISPLAYS AND VISION SYSTEMS

SYLLABUS INTRODUCTION TO ROTARY WING FLYING QUALITIES AND PERFORMANCE

Safety Analysis Tool for Automated Airspace Concepts (SafeATAC)

PBN and airspace concept

30 th Digital Avionics Systems Conference (DASC)

B.S. PROGRAM IN AVIATION TECHNOLOGY MANAGEMENT Course Descriptions

New Engine Option (A330neo) airplanes. These airplanes will have a novel or unusual design

SIMULATION OF ON-BOARD MODEL OF AIRLINER TO EVALUATE CAPABILITY OF TRAJECTORIES AND FLIGHT SAFETY

OPERATIONS CIRCULAR 01/2012. Subject: HEAD-UP DISPLAYS (HUD) AND ENHANCED VISION SYSTEMS (EVS)

ONLINE DELAY MANAGEMENT IN RAILWAYS - SIMULATION OF A TRAIN TIMETABLE

Transcription:

24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES AUTONOMOUS FLIGHT CONTROL AND GUIDANCE SYSTEM OF ACCIDENT AIRCRAFT Shinji Suzuki*, Fumihiro Kawamura**, Kazuya Masui** *The University of Tokyo, **Mitsubishi Heavy Industries Ltd, ***Japan Aerospace Exploration Agency (JAXA) Keywords: Safety, Control and Guidance, Optimization Abstract The purpose of this paper is to present a research program for autonomous control flight system of an aircraft in accidents or in trouble. The system will be expected to reduce the total loss number of accident aircrafts. Research subjects to be developed are automatic fault detection and identification, restructure algorithm of control system, and automatic guidance for the accident aircraft. System evaluations in computer simulation and flight test using an in-flight simulator are conducted. Three research groups from university, industry and national research laboratory are working together to integrate those technologies. This paper describes the outline of this research program, and the detail of the newly developed online flight path optimization method. 1 Introduction The safety of airliners is highly guaranteed throughout design, manufacturing, and operation processes. As aviation continues to grow, there is a concern that unless steps are taken to drastically reduce accident rates, increased flights will lead to more accidents. The National Civil Aviation Review Commission endorsed a goal to cut the fatal accident ratio by 80 percent by 2007[1]. While there are many reasons of fatal accidents of aircrafts, this paper considers the loss of control due to the failures of control surfaces. The failures include lock in place, hard over, and loss of effectiveness of control surfaces. These cases deteriorate the handling quality and increase the pilot workload, and then may lead to fatal crash. Autonomous flight control systems have been widely investigated to increase aircraft safety in civil and military aircrafts[2,3,4,5]. The autonomous flight control systems have the capability to detect and identify failures, to reconstruct control law, and to guide the damaged aircraft for safety flight or landing. The Society of Japanese Aerospace Companies has been promoting the research on autonomous flight control systems in corporation with industrial, academic and national research groups. The main feature of this research program is to demonstrate the total system capability including autonomous flight control and guidance systems for civil aircrafts throughout various simulation studies and flight tasting. This paper will outline the research program and focus on the online optimal flight trajectory design method and its computer simulations. The flight testing of the proposing method has been carried out from 2003 to 2004 by using in-flight simulator MuPAL- α and its results will be presented in another paper in this conference. 2 System Concept Present advanced aircrafts usually equip automatic flight control system that can greatly reduce pilot s workload and makes a significant contribution to air transportation safety. However, the automatic control system fundamentally assumes that every system is normal and has no trouble. Consequently, when aircrafts have some troubles or accidents, pilots have to manage to control their aircraft immediately. In order to control damaged or 1

SUZUKI S, KAWAMURA F, and MASUI K accident aircraft, firstly, pilots have to recognize or estimate damaged parts. If control systems such as engines or control surfaces have some troubles, pilots have to stabilize the aircraft by themselves. In some cases, this maneuver is very difficult since it is difficult to pinpoint damaged parts or the dynamic characteristics are greatly changed. Furthermore, if pilots have not experienced the troubles, the maneuver may be very difficult for them. Secondly, pilots have to guide the damaged aircraft to safety zones. If damaged aircrafts should be operated to lower altitude, or be landed to a nearest airport, pilots have to determine an adequate flight course in a short time. This research was intended to develop a new system that can cope with the above problem autonomously or can assist pilots by providing necessary information. Figure 1 indicates the proposing system concept that is divided into the following items to be developed. 2.1 Fault Detection and Identification System(FDIS): By comparing the normal states computed in the original system model with the real output data of an accident aircraft, the location and status are identified to minimize the difference between model output and flight data. 2.2 Restructurable Flight Control System(RFCS): The flight control law is restructured by using the information of damaged parts that are estimated in FDIS. The purpose of RFCS is to stabilize an accident aircraft by increasing the robustness of the original control law and restructuring the distribution law of control power. 2.3 Online Flight Trajectory Optimization System(OFTOS): Flight trajectories for emergency are computed by considering the system accidents estimated in FDIS. Real Time Algorithms are applied to optimize the flight path by onboard computers. 3 System Description Mathematical background of three systems Fig. 1. Future Autonomous Flight Control System 2

AUTONOMOUS FLIGHT CONTROL AND GUIDANCE SYSTEM OF ACCIDENT AIRCRAFT FDIS(Fault Detection and Identification System), RFCS(Resutructured Flight Control System), and OFTOS(Online Flight Trajectory Optimization System) will be explained briefly. 3.1 FDIS(Fault Detection and Identification System) FDIS is formulated as a system identification problem. The fundamental idea is to minimize the error between the real flight data measured from a set of sensors and the computed flight data from a mathematical flight model with system parameters to be estimated (as shown in Fig.2)[6]. While there is a wide range of faults or accidents, this study considers the damage of control surfaces. Therefore, the stability derivatives of each control surface must be identified by an online processing. The following three methods are developed for this purpose. MSLS(Modified Sequential Least Mean Square Method) A Sequential Least Mean Square Method is widely applied for linear system identification problems. Three linear momentum equations are considered for this problem. All stability derivatives including those of control surfaces and constant disturbances are estimated by using sensor output data, e.g., the time derivatives of pitch rate, roll rate and yaw rate signals. The mean square error is minimized sequentially with consideration of the constraints for each stability derivative. NN(Neural Network) The mean square error used in MSLS can be minimized by using the Backward Propagation method that is utilized in Neural Networks(NN). Although the linear model is used in this study, the NN has a potential for nonlinear estimation problems. EKF(Extended Kalman Filter) For nonlinear estimation problems, the Extended Kalman Filter is widely applied. This method linerlizes nonlinear dynamic equations and applies the Kalman Filter to obtain the system parameters that consist of control stability derivatives and disturbances. The Kalman Filter is designed to minimize the estimation error by considering the system noise and sensor noise. It is widely known that the adequate artificial inputs are necessary to obtain accurate estimations for any kind of estimation methods. It is noted that all of control surfaces are adequately moved during the Identification Process. Fig. 2. FDIS(Fault Detection and Identification System) 3.2 Restructurable Flight Control System(RFCS): The Flight Control System must be modified to compensate the deteriorated handling quality for damaged aircrafts. While a wide range of methods is proposed for this purpose, the robust control design[7] and the control mixing method are used in this research since the damage of control surface is considered[8]. RFCS is depicted in Fig.3 where FDIS sends information to the controller to restructure the control law. RFCS generates pseudo commands those are distributed to the deflection angle of each control surface. RFCS is divided into the longitudinal and lateral controllers as shown in Fig. 4. The pseudo inverse matrix is used to divide the pseudo commands into actual deflection angle of control surface. The information of estimated stability derivatives of each surface is utilized to construct the pseudo inverse matrix. 3.3 Online Flight Trajectory Optimization System(OFTOS) The guidance and the navigation of damaged aircrafts increase the pilot workload unusually. Consequently, the automatic 3

SUZUKI S, KAWAMURA F, and MASUI K generation of flight paths for the emergency flight is strongly recommended. Fig. 3. Restructurable Flight Control System(RFCS): Fig. 4. Distributor Robust Controller and Control The determination of flight trajectories satisfying several constraints is mathematically formulated by the optimal control theory that minimizes or maximizes a certain objective function with constraints. The optimal control theory is defined by the calculus of variations and some numerical computational methods have been developed[9]. Since these methods are time consuming and are not suitable for online computation, a new method is required for the on-line trajectory optimization. Several on-line methods based on heuristic search such as A* algorithm have been developed for path search problems[10]. Since these methods assume simple path search, they are not suitable for emergency flight that may require unusual flight patterns. This paper proposes the new method that combines the RTA* algorithm and the R-TABU search method[11]. The RTA* algorithm sequentially applies A* algorithm for on-line path search and the R-TABU search method is a kind of random search that controls search regions to improve computational efficiency. The new algorithm sequentially generates the new way point as shown in Fig 5 where the way point n and the flight path to the point n were determined in the previous step. The next way point n i and the flight path from the way point n to the point n i are computed before an airplane arrives at the way point n i. The new way point and the new flight path are determined to minimize the following objective function; J = g( n, ni) + h'( ni ) (1) where g(n,n i ) is an objective function calculated along the flight path from the way point n to n i, and h (n i ) is a heuristic function that estimates an objective function from the way point n i to the terminal point by using simple or heuristic methods. It should be noted that a set of constrains such as terminal conditions, and path constraints is included into objective functions by using a penalty function method. The R-TABU method is applied to minimize the objective function in Eq. (1). The outline of this method is illustrated in Fig. 6. This method uses sampling points selected at random around the initial solutions, and then the search range is shrunk gradually. The best solutions selected in this process are used for the next search. The parameters to be optimized are the coordinates at the way point n i and coefficients that determine the flight path from the way point n to the point n i. It is noted that the flight path is represented as a set of trigonometric functions. For example, the x coordinate is written as: N 0 + ( An cos t + Bn sin t) n= 1 x( t) = A ω ω (2) where the fundamental angular frequency is defined as ω = π /( tn i tn ). Since the flight path is represented as functions of time, the angle of attack, the bank angle, the load factor, and the thrust can be calculated by differentiating the flight path functions with considering a point mass model 4

AUTONOMOUS FLIGHT CONTROL AND GUIDANCE SYSTEM OF ACCIDENT AIRCRAFT of an airplane. The upper and lower limits of these values must be specified along the flight path to obtain the optimal flight trajectory. The way point is extended sequentially according the above process. When an airplane nears the terminal point within the distance R ref, the final flight path to the terminal point is computed. In that case, the terminal conditions such as the coordinates, the velocity vector are specified. 3 Numerical Simulation of Online Flight Trajectory Optimization System(OFTOS) The flight test using the experimental aircraft MuPAL- α [12] is conducted in this research program. Since the flight test results are reported in another paper[13], the numerical simulations carried out for the preparation of the flight test are summarized in this paper. MuPAL- α (Fig. 7) is based on Do-228 and has various capability and equipments for experiments. The initial and the terminal conditions are specified as shown in Table 1, when the emergency landing is considered. Table 1 Initial and Terminal conditions Initial Terminal Coordinate (10000, (0,0,100) (x,y,z)(m) -15000,1000) Velocity (m/s) 61.7 61.7 Path angle (deg) 0-1 Direction angle 170 90 (deg) Fig. 5. RTA* combined R-TAB method Fig. 6. R-TABU Search method Fig. 7. Multipurpose Experimental Airplane, MuPAL-α [12] Table 2 Constraints Minimum Maximum Velocity (m/s) 46.3 77.2 Angle of attack -5 15 (deg) Bank angle (deg) -20 20 Load factor 0.5 1.5 Thrust (N) -560~1040 10800~14800 Height (m) 0 - Table 2 shows the constraints specified in the trajectory optimization. It is noted that the maximum and minimum thrust values are assumed to be the linear functions of the flight velocity. The objective function to be minimized is the flight time from the initial to the terminal, since the emergency landing is considered in this simulation. Table 3 summarizes the parameters used in the optimization computation. Two different types of constraints and their optimized flight path are considered so as to check our method. Case 1 is a normal case and Case 2 is an emergency case that has different constrains, i.e., the deep left bank is limited. Two cases are compared in Table 4. 5

SUZUKI S, KAWAMURA F, and MASUI K Figure 8 illustrates the initial and terminal points and calculated flight paths. The normal aircraft (Case 1) firstly turn right and then turn left in order to minimize the flight time (337.3 sec). On the other hand, the emergency aircraft avoids a deep left turn, which leads a longer flight time (374.9 sec). Figure 9 that compares the flight velocity shows both aircraft try to maintain the maximum velocity to minimize the flight time, and reduce the velocity in the final stage in order to achieve the terminal velocity. Figures 10 and 11 illustrate the thrust and the angle of attack time histories. It is shown that the maximum thrust is not necessary since the paths are descent flights and the velocities reach the maximum values. However, the frequent thrust control is necessary in the final stage to accommodate the terminal conditions. Figure 11 shows the bank angle histories for two cases. It is confirmed that the strong constraint in the left bank angle is satisfied in the Case2. Table 3 Parameters in Optimization No. of samplings in one sub region 10 Shrink ratio of search sub region 3 Maximum no. of sampling in one search 100 Order of trigonometric function 5 Search time in RTA* (s) 40 Terminal search range in RTA*(m) 9000 Table 4 Two Cases Case1 Case2 Status Normal Accident Min left bank (deg) -20-7 Max right bank 20 20 (def) Flight time (s) 337.3 374.9 Fig. 8. Optimized Flight Paths (Case1: normal, Case 2: damaged aircraft) Fig. 9. Flight velocity (Case1:normal, Case2:Accident) Fig. 10. Thrust (Case1:normal, Case2:Accident) Fig. 11. Angle of Attack (Case1:normal, Case2:Accident) 6

AUTONOMOUS FLIGHT CONTROL AND GUIDANCE SYSTEM OF ACCIDENT AIRCRAFT Fig. 12. Bank Angle (Case1:normal, Case2:Accident) 4 Conclusion The research on the autonomous flight control and guidance systems of accident aircrafts have been conducted in Japan by the cooperative groups from industrial, academic, and governmental research institutes. The total system is basically divided into three subsystems, i.e., (1) Fault Detection and Identification System (FDIS), (2) Restructurable Flight Control System (RFCS), and (3) Online Flight Trajectory Optimization System (OFTOS). Wide range of simulation studies and the flight test using the experimental research airplane have been carried out. The present paper focuses on the simulation study on OFTOS (Online Flight Trajectory Optimization System), since FDIS and RFCS have been widely reported and the flight test will be presented in another paper[13]. On-line trajectory optimization techniques were applied to determine the emergency flight path with the consideration of aircraft accidents or failures. Since the unsteady flight dynamics should be considered in the emergency flight, the optimization problem should be solved in a strict manner by using onboard computers. The new method combining the RTA* algorithms with the R-TABU search makes it possible to provide the optimal flight path during the flight. The effectiveness of the presented method was confirmed from the simulation study. The flight test using an experimental airplane was conducted in our research, and the availability of the system was also confirmed, which will be presented in another paper[13]. The flight test of the total system including FDIS, RFCS, and OFTOS will be planned to develop the advanced autonomous flight control and guidance system. Finally it should be mentioned that this research program was funded by The Ministry of Economy, Trade and Industry (METI), and was organized by The Society of Japanese Aerospace Companies (SJAC). 5 References [1] NASA Aeronautical Blueprint, Toward a Bold New Era in Aviation, http://www.aerospace.nasa.gov/aboutus/tf/aero_bluep rint/index.html [2] Burken J J, and Burchman Jr. FW. Flight-Tests of Propulsion-Only Emergency Control System on MD- 11 Airplane, J. of Guidance, Control, and Dynamics, 20(5), pp. 980-987, 1997. [3] Ochi Y, and Knai K. Application of Restructurable Flight Control System to Large Transport Aircraft, J. of Guidance, Control, and Dynamics, 18(2), pp.365-370, 1995. [4] Brinker J S, and Wise K A. Flight Testing of Reconfigurable Control Law on the X-36 Tailless Aircraft, J. of Guidance, Control, and Dynamics, 24(5), pp. 903-909, 2001. [5] Bodson M. Reconfigurable Nonlinear Autopilot, J. of Guidance, Control, and Dynamics, 26(5), pp. 719-727, 2003. [6] Hamel P G, and Jategaonkar R V. Evaluation of Flight Vehicle System Identification, J. of Aircraft, 33(1), pp. 9-28, 1996. [7] McFarlane D C, Glocer K. Robust Controller Design Using Normalized Coprime Factor Plant Descriptions, Lecture Notes in Control and Information Series, Springer-Verlag, 1990. [8] Kameyama T, Senuma T, Mtaumoto Y, Ide M, Shimizu E, Hisano T, Yasui H, and Goto K. The Study on Self Repairing Flight Control System- Overview of Self Repairing Flight Control Law(in Japanese), Proceedings of the 36 th Aircraft Symposium, pp. 457-460, 1998. [9] Betts J T. Survey of Numerical Methods for Trajectory optimization, J. of Guidance, Control and Dynamics, 21(2), pp. 193-207, 1998. [10] Korf, R E. Real-Time Heuristic Search, Artificial Intelligence, 42(2-3), pp. 289-311, 1990. [11] Hu, N., Tabu Search Method with Random Moves for Globally Optimal Design, Int. J. of Numerical Methods in Engineering, 35, 1992, 1055-1070. [12] Masui K and Tsukano Y. Development of a new inflight simulator MuPAL-α. Proc AIAA Modeling and Simulation Technologies Conference and Exhibit, Denver, AIAA-2000-4574, pp 1-11, 2000. 7

SUZUKI S, KAWAMURA F, and MASUI K [13] Masui K, Tomita H, and Komatsu Y. Flight Experiment on a Flight Optimization for Aircraft in Trouble, ICAS2004 paper. 8