Next Article in Journal
Smart Blade Count Selection to Align Modal Propagation Angle with Stator Stagger Angle for Low-Noise Ducted Fan Designs
Next Article in Special Issue
Flight Simulation of Fire-Fighting Aircraft Based on Multi-Factor Coupling Modeling of Forest Fire
Previous Article in Journal
Parametric Design Method and Lift/Drag Characteristics Analysis for a Wide-Range, Wing-Morphing Glide Vehicle
Previous Article in Special Issue
Gliding Footprint Calculation Method for Aircraft with Thrust Failure Based on Six-Degree-of-Freedom Flight Envelope and Back-Propagation Artificial Neural Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fuzzy Modeling Framework Using Sector Non-Linearity Techniques for Fixed-Wing Aircrafts

by
Pablo Brusola
1,
Sergio Garcia-Nieto
1,*,
Jose Vicente Salcedo
1,
Miguel Martinez
1 and
Robert H. Bishop
2
1
Instituto de Automática e Informática Industrial, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
2
College of Engineering, University of South Florida, 4202 East Fowler Avenue ENB118, Tampa, FL 33620-5350, USA
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(4), 258; https://doi.org/10.3390/aerospace11040258
Submission received: 30 January 2024 / Revised: 13 March 2024 / Accepted: 20 March 2024 / Published: 26 March 2024
(This article belongs to the Special Issue Advanced Aircraft Technology)

Abstract

:
This paper presents a mathematical modeling approach utilizing a fuzzy modeling framework for fixed-wing aircraft systems with the goal of creating a highly desirable mathematical representation for model-based control design applications. The starting point is a mathematical model comprising fifteen non-linear ordinary differential equations representing the dynamic and kinematic behavior applicable to a wide range of fixed-wing aircraft systems. Here, the proposed mathematical modeling framework is applied to the AIRBUS A310 model developed by ONERA. The proposed fuzzy modeling framework takes advantage of sector non-linearity red techniques to recast all the non-linear terms from the original model to a set of combined fuzzy rules. The result of this fuzzification is a more suitable mathematical description from the control system design point of view. Therefore, the combination of this fuzzy model and the wide range of control techniques available in the literature for such kind of models, like parallel and non-parallel distributed compensation control laws using linear matrix inequality optimization, enables the development of control algorithms that guarantee stability conditions for a wide range of operations points, avoiding the classical gain scheduling schemes, where the stability issues can be extremely challenging.

1. Introduction

Piloting a modern aircraft, including of crewed, remotely piloted, and fully autonomous types, relies on automatic control systems to insure robust stability and high performance under a wide variety of conditions. Robust stability, often measured by phase and gain margin, is a closed-loop concept that quantifies the ability of a feedback control system to remain stable in the presence of uncertainties, both in the mathematical model and in cases of unexpected changes in the external environmental conditions. The level of performance required is application-specific and often measured in terms of rise time, overshoots, settling time, and other time-domain performances. Many phases of flight, such as final approach with autonomous landing, demand that many variables representing the aircraft attitude state (e.g., aircraft body reference frame relative to a ground reference frame) and translational state (e.g., position and velocity relative to the ground) are controlled simultaneously. The level of complexity of any model-based control system design process coupled with the ability to guarantee robust stability and the requisite closed-loop performance depends on an accurate mathematical aircraft model describing both the translational motion and the orientation. In general, the more closely the mathematical model reflects the specific aircraft and the environment it flies in (that is, the real world), the better the actively controlled closed loop system performance. Improvements in control system design methodologies stemming from the availability of accurate yet tractable mathematical models when combined with effective navigation systems directly contributes to safe, high-performance flight.
The classical mathematical model of the aircraft translation and orientation comprises fifteen non-linear differential equations [1]. Many control system design methodologies employ linearized mathematical models based on Taylor series approximation of the fifteen non-linear differential equations about an equilibrium point representing a flight regime of interest. For example, it is a well known and often studied problem of providing a smooth flight of a passenger aircraft in the presence of wind turbulence. The mathematical model can be linearized using wing level; constant velocity flight condition and the feedback control system design based on linear mathematical models strive to keep the flight smooth and near that equilibrium condition. Large deviations from the equilibrium condition, say, due to an unexpected large wind shear, can move the aircraft out of the linear region and potentially lead to performance and stability issues. In a single flight, there are many flight regimes from take-off to landing represented by different equilibrium points. One of the extended solutions for this issue is the combination of different controllers with the same structure (i.e., gain scheduling) [2].
This paper presents dynamic models for aircraft control design, which implies accurate models with suitable structures for optimizing the control design stage. One of the aims of this paper is to present a general modeling framework based on fuzzy logic for control design rather than linearization around various equilibrium points. This work explores ways to obtain a general fuzzy model of a fixed-wing aircraft using the sector non-linearity technique due to its proven ability to describe dynamic behavior in the field of non-linear systems (page 6 of [3,4,5,6]). The significant advantage of obtaining a fuzzy representation of the original aircraft dynamic model is the possibility to design controllers like parallel (PDC) and non-parallel (non-PDC) distributed compensation control laws [5,7,8] using linear matrix inequality (LMI) optimization [9]. These control algorithms guarantee stability conditions for a wide range of operation points and flight conditions using the same controller, avoiding the classical gain scheduling schemes where stability issues can be highly challenging.
The application of fuzzy logic in the field of aerospace applications has been inherent to the emergence of the first contributions in this area, as shown in [10,11,12]. These earlier works explored the use of fuzzy logic as alternatives in control algorithms, state estimation or fault detection. Afterwards, several contributions have been presented where experimental data and clustering techniques have been combined with fuzzy logic for the identification of dynamic models (black-box type) and/or controllers. Some examples are [13,14,15].
Since the use of fuzzy logic has provided satisfactory results almost from its first applications, it has been widely used in many different areas of the aerospace sector such as dynamic control in hyper-sonic aircraft surfaces [16], health monitoring for aircraft engines [17], autopilots for helicopters [18] or soft sensors for angle-of-attack (AoA) estimation [19,20]. Furthermore, the application of fuzzy logic in novel areas of the aerospace sector, such as Unmanned Aerial Vehicles (UAVs) or Autonomous Underwater Vehicles (AUVs), is also currently an area of interest, as shown in [21,22,23,24,25,26].
The approach presented in this work is to define a global non-linear fuzzy model representing the global aircraft dynamics. In general, there are two approaches for constructing fuzzy models: (i) identification using input–output data (as mentioned before [13,14,15,27]), and (ii) derivation from given non-linear system equations. The proposed fuzzy modeling framework is focused on the second approach and uses the idea of sector non-linearity. This method was introduced by Tanaka and Wang in (page 10 of [3]), and it enables the formulation of an equivalent fuzzy model from an original non-linear system using linear subsystems. This approach has been addressed previously in the aerospace field, as shown, for example, in [26,28,29].
The sector non-linearity technique is based on representing any static non-linear term k ( t ) of a model with a fuzzy quantity, where the jth rule is of the form
IF  z n 1 ( t ) is M z n 1 i n 1 and … and z n s ( t ) is M z n s i n s
THEN k ( t ) = A k j ( t ) j = 1 , , 2 s j = i n 1 + i n 2 · 2 1 + + i n s · 2 s 1 , i n h { 1 , 2 } h = 1 , , s ,
where z n h ( t ) are the premise variables, M z n h i n h are the membership functions and A k j ( t ) are constant values (zero-order model) or linear functions depending on state variables x ( t ) (first-order model) of non-linear models. The premise variables are known functions which may depend on the state variables and/or time [30].
In this paper, the main objective is to present a new version of the conventional fixed-wing aircraft dynamic model that uses fuzzy logic methods, in particular a way to obtain an equivalent and complete fuzzy model of such non-linear model. This is the main contribution of this work, which distinguishes it from other contributions such as [28,29], where the dynamics of the aircraft are only partially described (longitudinal model), the rotation matrices between the different axes are not included and simplifications are performed through the application of Taylor series to make the problem more solvable. For all these reasons, this article proposes a much more general framework that allows to obtain a global and complete model with a higher accuracy.
Finally, in future research, a Takagi–Sugeno fuzzy model (page 6 of [3]) will be developed starting from the fuzzy model obtained in this manuscript. This Takagi–Sugeno model will be used to design non-linear TS fuzzy controllers. On the one hand, for a better understanding, the validation and examples of the proposed fuzzy modeling are performed using parameters of an Airbus A310 aircraft [31].
This paper is organized as follows. In Section 2, the classical aircraft dynamic model is presented, where the attitude dynamics are represented by quaternions. Section 3 introduces the sector non-linearity technique for fuzzy modeling. Section 4 provides the construction of the fuzzy model by selecting model limits and applying them in the dynamic model. Section 5 presents simulation results and validation with the Onera non-linear model benchmark of A310 [31]. Section 6 offers conclusions.

2. Classic Aircraft Model with Quaternions

2.1. Quaternions in Aircraft Dynamic Models

It is common to use Euler angles to represent aircraft orientation. They have a very physical interpretation, more readily enabling visualization of orientation by humans. However, when utilizing Euler angles, it is possible to experience the so-called gimbal lock. For example, if using a roll–pitch–yaw description of the orientation, gimbal lock is a numerical condition that occurs when the pitch angle is ±90°. A potential solution, therefore, is to represent the orientation in some other way. Quaternions are widely used as attitude representation parameters for rigid bodies such as aircrafts [32,33]. For the dynamic aircraft models in this paper, quaternions are used to represent the orientation between the various reference frames due to their computational efficiency and excellent numerical properties. Computing Euler angles and associated coordinate transformation matrices to and from quaternions is readily accomplished. The three main reference frames utilized here are the body axis reference frame, the Earth-linked reference frame, and the wind axis reference frame described below. Quaternions are used to represent the orientation of the various Earth-linked vertical reference frames to the aircraft body axis reference frames, as well as the orientation between the aircraft wind axes reference frame and the aircraft body axes reference frame (Figure 1).
  • Body axis reference frame [ x b , y b , z b ] : is the origin at the center of gravity fixed to the aircraft structure. The x-axis is the longitudinal aircraft axis in which the positioning direction is forward, the y-axis is perpendicular to the aircraft plane of symmetry-positive direction out the right wing, and the z-axis is perpendicular to the other two and oriented downward.
  • Earth-linked axis reference frame [ x v , y v , z v ] : moves with the aircraft and its origin; it is also the center of gravity. The x-axis is oriented to the North, the y-axis to the East and z-axis perpendicular to the other two and oriented downward. The x-y plane is parallel to the Earth’s surface.
  • Wind axis reference frame [ x w , y w , z w ] : is a particular body axis frame because the x-axis is aligned with the airspeed velocity vector which is always tangent to the trajectory. The y-axis is perpendicular to the x-axis in the same plane as x b and y b , and the z-axis is perpendicular to the other two and oriented downward.
To increase the numerical robustness of the aircraft attitude, that is to say, to avoid the so-called gimbal lock, a set of parameters is formulated using the same reference frame in each equation of the aircraft dynamic model. For that reason, it is important to define the quaternion transformation coordinate matrix. The transformation matrix from body axis reference frame to the Earth-linked frame [32,33] is represented by 
R b v ( Q ( t ) ) = q 0 2 + q 1 2 q 2 2 q 3 2 2 ( q 1 q 2 q 3 q 0 ) 2 ( q 1 q 3 + q 2 q 0 ) 2 ( q 1 q 2 + q 3 q 0 ) q 0 2 q 1 2 + q 2 2 q 3 2 2 ( q 2 q 3 q 1 q 0 ) 2 ( q 1 q 3 q 2 q 0 ) 2 ( q 2 q 3 + q 1 q 0 ) q 0 2 q 1 2 q 2 2 + q 3 2
where Q ( t ) = [ q 0 , q 1 , q 2 , q 3 ] T is the attitude quaternion that satisfies | | Q ( t ) | | = 1 for all t. Matrix R b v ( Q ( t ) ) is an orthonormal matrix, hence R v b ( Q ( t ) ) = R b v ( Q ( t ) ) T . To visualize the aircraft attitude, it is useful to extract the Euler angles from the quaternion. For a roll–pitch–yaw Euler angle sequence, we have [32,33]
Φ = ϕ ( t ) θ ( t ) ψ ( t ) = tan 1 ( 2 ( q 2 q 3 + q 0 q 1 ) q 3 2 q 2 2 q 1 2 + q 0 2 ) sin 1 ( 2 ( q 1 q 3 q 0 q 2 ) ) tan 1 ( 2 ( q 2 q 1 + q 0 q 3 ) q 3 2 q 2 2 + q 1 2 + q 0 2 )
where ϕ , θ and ψ are the roll, pitch, yaw angles formed by the rotations of the x b , y b and z b axis, respectively.
The aircraft dynamic model is given by force, moment, and kinematic equations [1]. The aircraft attitude kinematics are given by
Q ˙ ( t ) = 1 2 T ( Q ( t ) ) Ω ( t )
where
T ( Q ( t ) ) = q 1 q 2 q 3 q 0 q 3 q 2 q 3 q 0 q 1 q 2 q 1 q 0
and Ω ( t ) = [ p ( t ) , q ( t ) , r ( t ) ] T is the aircraft angular velocity vector describing the angular velocity of the body reference frame with respect to the Earth-linked reference frame and represented in the body reference frame.
Here, we present useful relationships between the various forms of attitude representation (Euler angles, quaternions and transformation matrices) and the three key reference frames (body axis reference frame, the Earth-linked reference frame and the wind axis reference frame). We first consider the attitude representation of the wind axis reference frame to the body axis reference frame. The angle of attack, denoted by α ( t ) , represents rotations around the y-axis and the sideslip angle, denoted by β ( t ) , representing rotations around the z-axis. The transformation matrix from the wind axis reference frame to the body axis reference frame is given in [32]
R w b α ( t ) , β ( t ) = cos α ( t ) cos β ( t ) cos α ( t ) sin β ( t ) sin α ( t ) sin β ( t ) cos β ( t ) 0 sin α ( t ) cos β ( t ) sin α ( t ) sin β ( t ) cos α ( t )
The attitude quaternion that represents the orientation of the wind axis reference frame relative to the body axis reference frame is given by
q w b ( α ^ ( t ) , β ^ ( t ) ) = cos α ^ ( t ) cos β ^ ( t ) sin α ^ ( t ) sin β ^ ( t ) cos β ^ ( t ) sin α ^ ( t ) cos α ^ ( t ) sin β ^ ( t )
where we define half-angles α ^ ( t ) : = α ( t ) / 2 and β ^ ( t ) : = β ( t ) / 2 . We let q w b = q 0 q 1 q 2 q 3 T R 4 . Then, the transformation matrix from the wind axis reference frame to the body axis reference frame is given by
R w b ( q w b ) = 1 2 ( q 0 2 + q 1 2 q 2 2 q 3 2 ) q 1 q 2 q 0 q 3 q 1 q 2 + q 0 q 3 1 2 ( q 0 2 q 1 2 + q 2 2 q 3 2 ) q 0 q 2 q 1 q 3 2 q 2 q 3 q 0 q 2 q 1 q 3 0 1 2 ( q 0 2 q 1 2 q 2 2 + q 3 2 )

2.2. Aircraft Dynamic Model

The mathematical model describing the translation and orientation of a fixed-wing aircraft may be summarised as follows [1]:      
V ˙ ( t ) = 1 m F ( t ) Ω ( t ) × V ( t ) Ω ˙ ( t ) = I 1 ( M ( t ) Ω ( t ) × I Ω ( t ) ) Q ˙ ( t ) = 1 2 T ( Q ( t ) ) Ω ( t ) X ˙ ( t ) = R b v ( Q ( t ) ) V ( t )
where m is the total mass of the system, V ( t ) = [ u ( t ) , v ( t ) , w ( t ) ] T is the translational velocity in the body reference frame, Ω ( t ) = ( p ( t ) , q ( t ) , r ( t ) ) is the angular velocity represented in the body axis reference frame, X ( t ) is the translational position of the aircraft in the Earth-linked reference frame, I is the inertia matrix where I x x , I y y , I z z are the moments of inertia and I x z is a product of inertia. The products of inertia, I x y and I y z , related to longitudinal plane ( Y b = 0 ) are both neglected because of the aircraft’s symmetry with respect to this plane. Vectors F and M are the aircraft forces and moments, respectively.

2.2.1. Forces and Moments

The forces applied to the aircraft can be decomposed into three components: engine thrust, gravity and aerodynamic forces. We have F ( t ) = F e n g ( t ) + F g ( t ) + F a ( t ) , all expressed in the body axis reference frame. It is assumed that the thrust is aligned with the longitudinal axis body x-axis and given by
F e n g ( t ) = T ( t ) 0 0 ,
where T is the total engine thrust. The gravity force, naturally expressed in the Earth-linked axis frame and transformed to the body axis reference frame, is given by
F g ( t ) = R v b ( Q ( t ) ) 0 0 m g = m g 2 ( q 1 q 3 q 2 q 0 ) 2 ( q 2 q 3 + q 1 q 0 ) q 0 2 q 1 2 q 2 2 + q 3 2
where g is the gravity near the Earth’s surface. The aerodynamic forces are naturally expressed in the wind axis reference frame and transformed to the body axis as
F a ( t ) = q d S R w b ( q w b ) C = q d S R w b C X C Y C Z .
Aerodynamic coefficients C X , C Y and C Z are functions of the system variables, q d is the dynamic pressure and S is the aerodynamic reference area.
The moment about the center of gravity of the aircraft results from the engine and aerodynamic forces. So we have M ( t ) = M e n g + M a . The model is based on the A310 turbofan aircraft. For this reason, it is assumed that there are two engines delivering the same thrust and the distance of the midpoint of the motors to the center of gravity is null along the y-axis and the x-axis. It is assumed that the engines are located below the center of gravity at z e n g . The moment resulting from the thrust is
M e n g ( t ) = G E × F e n g = 0 z e n g · F e n g x ( t ) 0
where G E is the distance between the center of gravity “G” and the center of the gravity of engines “E”. The aerodynamic moment is composed of two main terms. The first is directly proportional to moment coefficients C l , C m and C n that correspond to roll, pitch and yaw moments, respectively. The second term is the moment resulting from aerodynamic forces F a ( t ) applied at the aerodynamic center, denoted by A, which may differ from the center of gravity. The total moment equation is given by
M a ( t ) = q d S c ¯ C l C m C n + G A × F a
where c ¯ is the aerodynamic mean chord and GA is the distance vector between the center of the gravity and the aerodynamic center, where it is assumed that the distance component is only about the x-axis, denoted by x a .

2.2.2. Aircraft Aerodynamic Model

In the Civilian Aircraft Landing Challenge [31], the aerodynamic coefficients, C L (lift coefficient), C Y (lateral coefficient), and  C D (drag coefficient), are given by
C L = C L 0 + C L α α ( t ) + c ¯ V a ( t ) C L q q ( t ) + C L δ e δ e + + C L H e λ L H L G ( t ) C Y = C Y β β ( t ) + C Y r δ r C D = C D 0 + C D α α ( t ) + C D α 2 α ( t ) 2
Similarly, the moment coefficients about the x-, y- and z-axes are given by
C l = C l β β ( t ) + c ¯ V a ( t ) ( C l p p ( t ) + ( C l r 0 + C l r α α ( t ) ) r ( t ) ) + + C l δ a δ a + C δ r δ r C m = C m 0 + C m α α ( t ) + c ¯ V a ( t ) C m q q ( t ) + C m δ e δ e + + ( C m H 0 + C m H α α ( t ) ) e λ m H L G ( t ) C n = ( C n β 0 + C n β α α ( t ) ) β ( t ) + c ¯ V a ( t ) ( C n r r ( t ) + + ( C n p 0 + C n p α α ( t ) ) p ( t ) ) + + C n δ a δ a + C n δ r δ r
Note that for C L and C m equations, the term that depends on the height of the main landing gear above the runway, denoted by H L G ( t ) , describes the ground effect. Terms δ e , δ a and δ r represent the deflections of the control surface of the elevators, ailerons and rudder, respectively. These equations are functions of angular rate variables, p ( t ) , q ( t ) , r ( t ) , the sideslip angle, β ( t ) , the angle-of-attack denoted by α ( t ) , and the true airspeed denoted by V a ( t ) . The airspeed is affected by wind, denoted as W = W x , W y W z T expressed in the Earth-linked axis reference frame. Airspeed V a = V a x ( t ) , V a y ( t ) V a z ( t ) T in the body axis reference frame is given by
V a = V ( t ) R v b ( Q ) W
Angle of attack α ( t ) and sideslip β ( t ) are computed via
α ( t ) = tan 1 ( V a z ( t ) V a x ( t ) ) a n d β ( t ) = sin 1 ( V a y ( t ) V a ( t ) )
Finally, the relationships between C L , C D and C X , C Z are given by
C X = C D , C Z = C L .

3. Sector Non-Linearity

Sector non-linearity in fuzzy model construction appears in page 10 of [3], and is based on the following idea. We consider a non-linear system with a static non-linear term defined by function f ( x ( t ) ) . The aim is to find the global sector such that f ( x ( t ) ) [ a 1 a 2 ] x ( t ) + [ b 1 b 2 ] . Figure 2 illustrates this approach. Sometimes it is difficult to find global sectors for general non-linear functions. In this case, we can consider a local sector non-linearity where two vertical lines become this local sector where d < x ( t ) < d , and where [ d d ] are the limits of variable x ( t ) as illustrated in Figure 3.
Therefore, it is necessary to establish numerical limits [ d d ] of the variables to be approximated by the local sector non-linearity technique. In the case of a fixed-wing aircraft, those ranges include variables such as angle of attack, sideslip and airspeed. The information detailed in operational flight manuals or other technical information provided by aircraft designers [34,35] details the valid operational ranges as valuable sources to establish the particular value of limit [ d d ] for each variable.
This approach is applied to obtain the fuzzy representation of static non-linear terms which appear in the model of fixed-wing aircraft according to Equation (1).
Example 1. 
The following non-linear function is the expression of angle-of-attack based on classic aircraft modeling:
α ( t ) = tan 1 x 1 ( t ) x 2 ( t ) = tan 1 ( s ( t ) )
s ( t ) = x 1 ( t ) x 2 ( t )
We obtain a fuzzy representation in two steps. In the first step, we deal with the equation between α and s. Figure 4 shows the graphical representation of α ( t ) where two bounds of the global sector might be [ 1 , 0 ] s ( t ) . Therefore, based on these two simple sectors, α ( t ) can be represented using Equation (1) with two rules:
α ( t ) = M α 1 s ( t ) + M α 2 0
where
M α 1 + M α 2 = 1
The M α 1 and M α 2 consequents are s ( t ) and 0; then,
A α 1 ( t ) = s ( t ) , A α 2 ( t ) = 0
Moreover, α ( t ) is the reference variable of the fuzzy model. Note that the membership function is dependent on α ( t ) ; then, it is important to remark that in this function, s ( t ) = tan ( α ( t ) ) . Therefore, the membership functions can be calculated as
M α 1 = α ( t ) t a n ( α ( t ) ) , M α 2 = 1 M α 1 = t a n ( α ( t ) ) α ( t ) t a n ( α ( t ) )
It should be noted that, in general, s 0 , but if this happens, we set M α 1 = 1 and M α 2 = 0 . For the rest of the paper, these membership functions M k i are termed Type I.
In the second step, we consider s ( t ) . In the angle-of-attack case, s = x 1 x 2 where x 1 ( t ) = V a z ( t ) and x 2 ( t ) = V a x ( t ) . In this non-linear case, the global sector cannot be defined, so a local sector isused instead, where the quotient is bounded by x 1 x 2 m i n , x 1 x 2 m a x .
The reason for these limits is the same as that mentioned at the beginning of Example 1. The outcome of the local sector non-linearity is this fuzzy quantity with two rules:
s ( t ) = x 1 ( t ) x 2 ( t ) = N α 1 x 1 x 2 m a x + N α 2 x 1 x 2 m i n
where membership functions are
N α 1 = x 1 ( t ) x 2 ( t ) x 1 x 2 m i n x 1 x 2 m a x x 1 x 2 m i n N α 2 = x 1 x 2 m a x x 1 ( t ) x 2 ( t ) x 1 x 2 m a x x 1 x 2 m i n
In this case, the  N α 1 and N α 2 consequents are x 1 x 2 m a x and x 1 x 2 m i n ; then,
A s 1 ( t ) = x 1 x 2 m a x , A s 2 ( t ) = x 1 x 2 m i n
All membership functions with this structure are defined as type IV and labelled as N k i .
For representing the α ( t ) fuzzy model, the chosen value limits were x 1 x 2 = V a z ( t ) V a x ( t ) [ 40 , 40 ] . In fuzzy models, it is usual to name the membership functions as “Positive”, “Negative”, “Zero” and “Not Zero” (see Figure 5 and Figure 6). Finally, the fuzzy model for α ( t ) model has four rules after combining fuzzy Euantities (23) and (27), and it is presented in Table 1, where
A α 1 = s | s = V a z V a x m a x =40
A α 2 = s | s = V a z V a x m i n = 40
A α 3 =0
A α 4 =0
The defuzzification is carried out as
α ( t ) = j = 1 4 h j A α j
where h 1 = M α 1 N α 1 , h 2 = M α 1 N α 2 , h 3 = M α 2 N α 1 , h 4 = M α 2 N α 2 . This fuzzy model is of the zero order because all the consequents are constants.
For further information regarding this example, see Example 3 in page 14 of [3].

4. Aircraft Fuzzy Model

As shown in Section 2.2, there is a large number of non-linear terms associated with aircraft dynamics. As a consequence, these terms are the main objective for applying the sector non-linearity methodology. Here, we provide a brief definition of all the types of membership functions and their consequent use following a similar procedure to Example 1 with membership Types I and IV. The different types on membership functions are chosen on the basis of the different structures of the most repeated non-linear terms.

4.1. Types of Membership Functions

Example 1 introduces membership Types I and IV, and following the same structure, this section describes Types II, III, V, VI and VII.

4.1.1. Type II: β = sin 1 ( x ( t ) )

As it can be observed in Figure 7, the function can be delimited by the global sector as β = sin 1 ( x ( t ) ) [ π 2 , 1 ] x ( t ) . The global sector non-linearity method results in the following membership functions of Type II as illustrated in Figure 8:
E β 1 = β ( t ) s i n ( β ( t ) ) s i n ( β ( t ) ) ( π 2 1 ) a n d E β 2 = π 2 s i n ( β ( t ) ) β ( t ) s i n ( β ( t ) ) ( π 2 1 )
where π 2 x ( t ) and x ( t ) are their consequents respectively. Consequently, this fuzzy model is of the first order.

4.1.2. Type III: z 1 ( t ) = x ( t )

This non-linear term appears in Equations (14) and (15). In this particular case, a local sector non-linearity is applied, where bounds are defined as [ x m i n , x m a x ] . Due to the non-negative value for the function, x m i n must be forced to 0. The membership functions can be defined as follows:
F z 1 1 = x ( t ) x m a x a n d F z 1 2 = x m a x x ( t ) x m a x
where F z 1 1 and F z 1 2 consequents are x m a x and x m i n = 0 respectively.

4.1.3. Type V: z 2 ( t ) = e x ( t )

To understand the fuzzy model from this type of equation, we consider that the exponential functions are used to model the aircraft ground effect. For this reason, variable x ( t ) that expresses the distance of the aircraft from the ground will never be negative. Therefore, the minimum value will be x ( t ) = 0 , offering the maximum of function z 2 ( t ) = 1 . In consequence, the first local sector is given by the x ( t ) + 1 line equation. The other global sector is defined by its asymptote when x ( t ) reaches high values, and when it is placed at z 2 = 0 . This global sector non-linearity is illustrated in Figure 9. The membership functions of Type V are described in the following Equations:
G z 2 1 = e x ( t ) 1 x ( t ) , G z 2 2 = 1 x ( t ) e x ( t ) 1 x ( t )
where G z 2 1 and G z 2 2 consequents are 1 x ( t ) and 0, respectively. The membership functions are illustrated in Figure 10.

4.1.4. Particularities of Type IV

Example 1 describes the membership function of quotient input terms as x 1 ( t ) x 2 ( t ) . These functions are denoted as N k i and are used in further membership structures. First, in quadratic terms ( x ( t ) 2 ), where the membership function structure is equivalent to type IV, their consequents are x m a x and x m i n , respectively. The linear equation terms also use the N k i form, and its consequents are a x m a x + b and a x m i n + b , respectively. For better understanding, in the summary Table 2 and Table 3, the membership function for quadratic terms is H k i (Type VI) and for linear functions, it is J k i (Type VII).

4.2. Fuzzy Application to the Aaircraft Model

Once the different fuzzy model types are defined, the aircraft dynamic model introduced in Equations (8)–(18) can be reformulated in terms of fuzzy logic. In this new structure, the non-linear terms (except the polynomial forms) can be found in two sets of equations: the wind effects and aerodynamic coefficients.

4.2.1. Wind Effects

The variables included in non-linear terms and related to wind effects are α ( t ) , β ( t ) and V a ( t ) . The  fuzzy model of the angle of attack, α ( t ) , is described in Table 1, the  fuzzy model of the sideslip angle, β ( t ) , in Table 2, and the fuzzy model of the true airspeed, V a ( t ) , in Table 3.

4.2.2. Aerodynamic Coefficients

The non-linear terms of the aerodynamic coefficients are reformulated here using the sector non-linearity method. In particular, Equation (28) introduces non-linear terms C L i , C D i , C l i , C m i and C n i where i refers the specific coefficient sub-term from the C L , C D , C l , C m and C n aerodynamic coefficients, where      
C L 1 = c ¯ V a ( t ) C L q q ( t ) C L 2 = C L H e λ L H L G ( t ) C D 2 = C D α 2 α ( t ) 2 C l 1 = c ¯ V a ( t ) C l p p ( t ) C l 2 = c ¯ V a ( t ) r ( t ) ( C l r 0 + C l r α α ( t ) ) C m 1 = c ¯ V a ( t ) C m q q ( t ) C m 2 = ( C m H 0 + C m H α α ( t ) ) e λ m H L G ( t ) C n 1 = c ¯ V a ( t ) C n r r ( t ) C n 2 = c ¯ V a ( t ) p ( t ) ( C n p 0 + C n p α α ( t ) ) C n 3 = ( C n β 0 + C n β α α ( t ) ) β ( t )
The aerodynamic coefficients in (32) are a combination of different non-linear functions of the types previously described. The fuzzy models of the ten aerodynamic coefficients are described in Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12 and Table 13. The maximums and minimums are left as unknown variables and the user should define them. This is because there are different bound conditions for each aircraft type and the objective of this paper is to propose a general model, that is to say, a full flight dynamics model for later application introducing all the aircraft constants.
Figure 11 presents the procedure to combine all the fuzzy models presented in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12 and Table 13 for obtaining a full fuzzy model for the aircraft.

5. Results and Validation

Due to the complexity of the formulation presented in this work and the large number of fuzzy terms involved, the Fuzzy Modelling Toolbox for Aircraft Systems (FMA Toolbox) is developed in Matlab. The FMA Toolbox reformulates fixed-wing aircraft models in terms of fuzzy logic, and it can be downloaded at GitHub [36]. Furthermore, to clarify the set of functions and scripts developed in this Toolbox, the detailed user manual can be consulted in [37]. Also, the main implementation of the fuzzy non-linear model for fixed-wing aircraft can be found in [38], where Equations (2)–(18) (introduced in Section 2) are coded in terms of fuzzy logic. This toolbox uses Algorithm 1.
Algorithm 1 Aircraft fuzzy model
1:
Set t = 0 . Set initial conditions of the model.
2:
while  t t m a x  do
3:
    Read the values of the inputs at t
4:
    Compute non-linear terms at t using fuzzy models from Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12 and Table 13
5:
    Substitute non-linear terms in aircraft model Equations (8).
6:
    Integrate these equations and obtain all the states and outputs at t + δ t
7:
    Set t = t + δ t
8:
end while

5.1. Inputs and Initial Conditions

Once the classical aircraft model is completely reformulated in terms of fuzzy logic, it is necessary to test and compare the behavior of these two models using real data. The aircraft dynamic model used here is based on Onera’s benchmark of the landing challenge [31]. Thus, it is logical to put the initial values of this model into the fuzzy model. Table 14 shows all the parameters of the Airbus A310 aircraft used here. It should be emphasized that the results are simulated in a short time frame and that is why the aircraft mass is considered constant, as well as the inertia matrix. The initial conditions of the simulation are calculated as classic Taylor series equilibrium points searching for the actuator trim values. For a more realistic model, these actuators are also modeled and their dynamics are approximated by magnitude and rate limit first-order filters whose characteristics are summarized in Table 15. To obtain significant simulation results, the actuators must vary with different combinations as shown in Figure 12. Finally, the fuzzy model needs maximum and minimum values to be defined. Table 16 defines very conservative limits; for example, α ( t ) and β ( t ) are bounded by a very large interval. Therefore, it is important to remark that some of these value limits can be selected more restrictive, if desired.

5.2. Simulation Results

To compare the non-linear dynamic model of Civilian Aircraft Landing Challenge [31] and the fuzzy logic model, both models are simulated under equivalent conditions. Particularly, simulations are carried out using the Toolbox developed in [36]. These simulations are executed on a PC with Intel Core i5 CPU at 1.6 GHz with 16 GB of DDR3 RAM and using MATLAB 2020B. A computation time of 72 s is necessary to obtain the dynamic evolution of the simulation model for 100 s.
The original non-linear dynamical model developed in Equation (8) uses a first-order algorithm with a fixed time step of 0.05 seconds, usually known as ODE1 [39], to solve the differential equations. For that reason, to validate and compare the fuzzy model proposed with the original one in equivalent conditions, the Toolbox developed uses the same first-order approximation (ODE1) as a default. However, higher-order solvers, such as ODE2 or ODE3 [39], are available in the Toolbox and can used to simulate the aircraft fuzzy model, always using a fixed time step. Inn addition, the aircraft fuzzy model is simulated by applying Algorithm 1.
On the other hand, the Variance Accounted For (VAF) index is used to provide an objective value of the percent of similitude between the two output models [12]. The VAF is governed by the following equation:
V A F i = 1 v a r ( y i y ^ i ) v a r ( y i ) 100 % ,
where y i is the measured output (classic model) and y ^ i is the estimated output (fuzzy model) for the ith component. Figure 13 and Figure 14 show the response of the state variables in both models. Note that the fuzzy modeling method is an exact method, and so the VAFs have high values. Nevertheless, the index is not 100% because the non-linear model of [31] uses the integration of the Euler angles while this integration is based on quaternions in the fuzzy model. In this simulation, the difference is insignificant due the short period of simulation. Therefore, the error increases as simulation time increases.

6. Conclusions and Future Works

A general framework for modelling fixed-wing aircraft using fuzzy structures and sector non-linearity techniques is proposed. The main contribution the obtention of a more general and accurate non-linear model than other alternatives proposed in the existing literature, such as [28,29]. Moreover, the approach presented allows for the expression of aerodynamic coefficients, moments, forces and wind effects, all of them expressed in terms of fuzzy logic and the body frame. This is possible thanks to the capability of the global and local sector non-linearity technique to approximate non-linear terms with fuzzy logic rules and with the same accuracy as the original non-linear model. In addition, the proposed quaternions form contributes to the good performance of the fuzzy logic model by removing discontinuities and more complex trigonometric functions. However, the good mathematical performance of the quaternions has the negative effect of losing the physical concept of comparing the Euler angles. Nevertheless, it is possible to find an equivalent Euler representation from quaternions.
The main limitation of obtaining the equivalent fuzzy model is the requirement to limit the maximum and minimum values for some of the system variables. This is due to the application of the local sector non-linearity technique, as detailed in Section 3 and Section 4. However, this limitation is not extremely constraining in the case of a fixed-wing aircraft, as there is often a significant amount of information available (theoretical and experimental) from manufacturers and the research community on the limits of these variables.
A Matlab toolbox is implemented [36] to simplify the application of this new approach for any standard fixed-wing aircraft defined by classical parameters such as geometry, mass, coefficients and bound conditions, where the outcome of the toolbox is a newly reformulated fuzzy model in terms of quaternions.
On the other hand, in Section 5, using this Matlab toolbox, a fuzzy equivalent model for the Onera non-linear model benchmark of the A310 [31] is obtained. Also, a comparison between the AIRBUS A310 model developed by ONERA and the obtained fuzzy model is performed. As a consequence, it is concluded that both models are equivalent.
The following research step consists of defining in terms of the Takagi–Sugeno structure (page 6 of [3]) the fuzzy model obtained in this work. This representation is suitable for designing PDC and non-PDC fuzzy controllers using well-known methodologies based on linear matrix inequalities [5,6,40].

Author Contributions

Conceptualization, P.B., S.G.-N. and J.V.S.; methodology, S.G.-N. and J.V.S.; software, P.B. and S.G.-N.; validation, P.B., S.G.-N., J.V.S., M.M. and R.H.B.; formal analysis, investigation, writing, review and editing, P.B., S.G.-N., J.V.S., M.M. and R.H.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Spain government via MCIN/AEI/10.13039/501100 011033 [project PID2020-119468RA-I00].

Data Availability Statement

The software and data presented in this work are available at https://github.com/sergarro/FuzzyAircraftModellingToolbox/tree/master.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

A k j rule consequent
E k i type II membership function
F k i type III membership function
G k i type V membership function
H k i type VI membership function
J k i type VII membership function
M k i type I membership function
N k i type IV membership function
inumber of membership function
jnumber of rule
kvariable reference of fuzzy model
α angle of attack, rad
β aerodynamic sideslip angle, rad
C D , C Y , C L drag, lateral, and lift force coefficients
C l , C m , C n roll, pitch and yaw moment coefficients
c ¯ aerodynamic mean chord, m
δ a , δ e , δ r ailerons, elevators, and rudder deflections, rad
EPRexhaust pressure ratio
Fresulting force vector acting on aircraft body, N
F g , F a , F e n g gravity, aerodynamic, and engine forces, N
F e n g x x component of the of the engines forces, N
ggravitational field intensity
near the Earth’s surface, m/s2
H L G landing gear height, m
Iaircraft inertia matrix, kg/m2
h j membership function product
I { x x , y y , z z } moments of inertia on { x b , y b , z b } axes, kg/m2
I { x y , x z , y z } products of inertia on { x b , y b , z b } axes, kg/m2
λ L , λ m lift and pitch ground effect coefficients
Mresulting moment vector acting on aircraft body, N.m
maircraft total mass, kg
Ω aircraft angular velocity vector, rad/s
p, q, raircraft angular X, Y, and Z velocity components with
respect to ground
and expressed in body axes, rad/s
Φ euler angles matrix, rad
ϕ, θ, ψaircraft Euler angle components, rad
Qdefinition of quaternion matrix
q w b attitude represented from wind axis to body axis
q d dynamic pressure, kg/(m.s2)
q 0 , q 1 , q 2 , q 3 quaternion components based on Euler angles
q 0 , q 1 , q 2 , q 3 quaternion components based on angle of attack and
sideslip angle
R b v reference system rotation matrix from body axis
to Earth-linked vertical frame
R w b reference system rotation matrix from wind axis
to body axis
Sprincipal wing surface, m2
Vaircraft velocity vector, m/s
V a true airspeed, m/s
VAFvariance accounted for index
u, v, waircraft X, Y, and Z velocity components with respect
to ground and expressed in body axes, m/s
x A X component of aerodynamic aircraft center position, m
x b , y b , z b components of body axis frame, m
x v , y v , z v components of earth-linked axis frame, m
x w , y w , z w components of wind axis frame, m
Wwind velocity vector, m/s
z e n g Z component of engine gravity center in body axis, m

References

  1. Klein, V.; Morelli, E.A. Aircraft System Identification: Theory and Practice; American Institute of Aeronautics & Astronautics: Reston, VI, USA, 2006; Chapter 3; pp. 40–42. [Google Scholar]
  2. Apkarian, P. Advandec Gain-Scheduling Techniques for Uncertain Systems. IEEE Trans. Control. Syst. Technol. 1998, 6, 209–228. [Google Scholar] [CrossRef]
  3. Tanaka, K.; Wang, H. Fuzzy Control Systems Design and Analysis: A Linear Matrix Inequality Approach; Wiley-Blackwell: Hoboken, NJ, USA, 2001. [Google Scholar]
  4. Angelov, P.P.; Filev, D.P. An approach to online identification of Takagi-Sugeno fuzzy models. IEEE Trans. Syst. Man Cybern. Part B 2004, 34, 484–498. [Google Scholar] [CrossRef]
  5. Wang, H.O.; Li, J.; Niemann, D.; Tanaka, K. TS fuzzy model with linear rule consequence and PDC controller: A universal framework for nonlinear control systems. In Proceedings of the Ninth IEEE International Conference on Fuzzy Systems. FUZZ-IEEE 2000 (Cat. No. 00CH37063), Padua, Italy, 18–23 July 2022; Volume 2, pp. 549–554. [Google Scholar]
  6. Wang, H.O.; Tanaka, K.; Griffin, M.F. An approach to fuzzy control of nonlinear systems: Stability and design issues. IEEE Trans. Fuzzy Syst. 1996, 4, 14–23. [Google Scholar] [CrossRef]
  7. Li, J.; Wang, H.O.; Niemann, D.; Tanaka, K. Dynamic parallel distributed compensation for Takagi–Sugeno fuzzy systems: An LMI approach. Inf. Sci. 2000, 123, 201–221. [Google Scholar] [CrossRef]
  8. Chang, W.J.; Ku, C.C.; Chang, C.H. PDC and Non-PDC fuzzy control with relaxed stability conditions for contintuous-time multiplicative noised fuzzy systems. J. Frankl. Inst. 2012, 349, 2664–2686. [Google Scholar] [CrossRef]
  9. Boyd, S.; El Ghaoui, L.; Feron, E.; Balakrishnan, V. Linear Matrix Inequalities in System and Control Theory; SIAM: Philadelphia, PA, USA, 1994. [Google Scholar]
  10. Larkin, L.I. A fuzzy logic controller for aircraft flight control. In Proceedings of the 23rd IEEE Conference on Decision and Control, Las Vegas, NE, USA, 12–14 December 1984; pp. 894–897. [Google Scholar]
  11. Chiu, S.; Chand, S. Fuzzy controller design and stability analysis for an aircraft model. In Proceedings of the 1991 American Control Conference, Brighton, UK, 11–13 December 1991; pp. 821–826. [Google Scholar]
  12. Oosterom, M.; Babuska, R. Virtual sensor for fault detection and isolation in flight control systems-fuzzy modeling approach. In Proceedings of the 39th IEEE conference on decision and control (Cat. No. 00CH37187), Sydney, NSW, Australia, 12–15 December 2000; Volume 3, pp. 2645–2650. [Google Scholar]
  13. Pan, M.; Wang, H.; Huang, J. T–S Fuzzy Modeling for Aircraft Engines: The Clustering and Identification Approach. Energies 2019, 12, 3284. [Google Scholar] [CrossRef]
  14. Weinstein, R.; Hubbard, J.E.; Cunningham, M. Fuzzy modeling and parallel distributed compensation for aircraft flight control from simulated flight data. In Proceedings of the 2018 Atmospheric Flight Mechanics Conference, Atlanta, GA, USA, 25–29 June 2018; p. 3313. [Google Scholar]
  15. Singh, D.J.; Verma, N.K.; Ghosh, A.; Sanwale, J.; Malagaudanavar, A. Aerodynamic Parameter Modeling Using TS Fuzzy Systems from Flight Data. In Computational Intelligence: Theories, Applications and Future Directions-Volume I: ICCI-2017; Springer: Singapore, 2019; pp. 425–437. [Google Scholar]
  16. Gao, D.; Sun, Z.; Du, T. Dynamic Surface Control of Hypersonic Aircraft Using Fuzzy Logic System. In Proceedings of the 2007 IEEE International Conference on Automation and Logistics, Jinan, China, 18–21 August 2007. [Google Scholar]
  17. Demirci, S.; Hajiyev, C.; Schwenke, A. Fuzzy logic-based automated engine health monitoring for commercial aircraft. Aircraft Eng. Aerosp. Technol. 2000, 80, 516–525. [Google Scholar] [CrossRef]
  18. Phillips, C.; Karr, C.L.; Walker, G. Helicopter Flight Control with Fuzzy Logic and Genetic Algorithms. Genet. Algorithms Fuzzy Log. Syst. Soft Comput. Perspect. 1996, 9, 175–184. [Google Scholar] [CrossRef]
  19. Singh, D.J.; Dev, R.; Verma, N.K.; Ghosh, A.; Malagaudanavar, A. Estimating angle of attack of an ATAS aircraft using TS fuzzy model. In Proceedings of the 2017 International Conference on Sensing, Diagnostics, Prognostics, and Control (SDPC), Shanghai, China, 16–18 August 2017; pp. 451–455. [Google Scholar]
  20. Singh, D.J.; Agrawal, P.; Verma, N.K.; Ghosh, A.; Malagaudanavar, A. Interval type-2 TS fuzzy model for angle of attack sensor of the aircraft. J. Intell. Fuzzy Syst. 2018, 34, 3891–3901. [Google Scholar] [CrossRef]
  21. Kadmiry, B.; Driankov, D. A fuzzy flight controller combining linguistic and model-based fuzzy control. Fuzzy Sets Syst. 2004, 146, 313–347. [Google Scholar] [CrossRef]
  22. Petritoli, E.; Leccese, F.; Cagnetti, M. Takagi-Sugeno Discrete Fuzzy Modeling: An IoT Controlled ABS for UAV. In Proceedings of the 2019 II Workshop on Metrology for Industry 4.0 and IoT (MetroInd4.0 IoT), Naples, Italy, 4–6 June 2019; pp. 191–195. [Google Scholar] [CrossRef]
  23. Mao, X.; Zhang, H.; Wang, Y. Flocking of quad-rotor UAVs with fuzzy control. ISA Trans. 2018, 74, 185–193. [Google Scholar] [CrossRef] [PubMed]
  24. Kayacan, E.; Maslim, R. Type-2 fuzzy logic trajectory tracking control of quadrotor VTOL aircraft with elliptic membership functions. IEEE/ASME Trans. Mechatronics 2016, 22, 339–348. [Google Scholar] [CrossRef]
  25. Rashid, M.; Sheikh, N.; Raza, A.; Raziq, A.; Rehman, A.R.A.; Baber, J.; Basit, A. Development of Dynamic Modeling and Fuzzy Logic System by Classical and Modern Strategies for the Control of Quadcopter. Am. Sci. Res. J. Eng. Technol. Sci. (ASRJETS) 2019, 62, 108–114. [Google Scholar]
  26. Leccese, F.; Cagnetti, M.; Giarnetti, S.; Petritoli, E.; Luisetto, I.; Tuti, S.; Durovic-Pejcev, R.; Dordevic, T.; Tomašević, A.; Bursić, V.; et al. A Simple Takagi-Sugeno Fuzzy Modelling Case Study for an Underwater Glider Control System. In Proceedings of the 2018 IEEE International Workshop on Metrology for the Sea, Bari, Italy, 8–10 October 2018; pp. 262–267. [Google Scholar] [CrossRef]
  27. Takagi, T.; Sugeno, M. Fuzzy Identification of Systems and Its Applications to Modeling and Control. IEEE Trans. Syst. Man, Cybern. 1985; SMC-15, 116–132. [Google Scholar]
  28. Hušek, P.; Narenathreyas, K. Aircraft longitudinal motion control based on Takagi–Sugeno fuzzy model. Appl. Soft Comput. 2016, 49, 269–278. [Google Scholar] [CrossRef]
  29. Tan, B.; Mu, J. Fuzzy Modeling and Control of Longitudinal Channel of F16 Aircraft at a Fixed Speed. In Proceedings of the 2018 2nd IEEE Advanced Information Management, Communicates, Electronic and Automation Control Conference (IMCEC), Xi’an, China, 25–27 May 2018; pp. 631–636. [Google Scholar]
  30. Tognetti, E.S.; Linhares, T.M. Dynamic output feedback controller design for uncertain Takagi–Sugeno fuzzy systems: A premise variable selection approach. IEEE Trans. Fuzzy Syst. 2020, 29, 1590–1600. [Google Scholar] [CrossRef]
  31. Biannic, J.M.; Boada-Bauxell, J. A Civilian Aircraft Landing Challenge. Onera Benchmark 2016. Available online: https://w3.onera.fr/smac/sites/w3.onera.fr.smac/files/CALC_v2.pdf (accessed on 30 January 2024).
  32. Diebel, J. Representing Attitude: Euler Angles, Unit Quaternions, and Rotation Vectors. Matrix 2006, 58, 1–35. [Google Scholar]
  33. Groß ekatthöfer, K.; Yoon, Z. Introduction into Quaternions for Spacecraft Attitude Representation. Ph.D. Thesis, Technical University of Berlin, Berlin, Germany, 2012. [Google Scholar]
  34. Lufthansa-A330. Flight Crew Operating Manual. 2024. Available online: https://www.avialogs.com/reader.php?jid=991#p=3 (accessed on 30 January 2024).
  35. Airbus-A310. A310 Flight Deck and Systems Briefing for Pilots. 2024. Available online: https://www.avialogs.com/aircraft-a/airbus/item/989-a310-flight-deck-and-systems-briefing-for-pilots (accessed on 30 January 2024).
  36. Brusola, P.; Garcia-Nieto, S. Fuzzy Modelling Toolbox for Aircraft Systems-Code. 2024. Available online: https://github.com/sergarro/FuzzyAircraftModellingToolbox.git (accessed on 30 January 2024).
  37. Brusola, P.; Garcia-Nieto, S. Fuzzy Modelling Toolbox for Aircraft Systems-User Manual. 2024. Available online: https://github.com/sergarro/FuzzyAircraftModellingToolbox/blob/master/User%20Manual.pdf (accessed on 30 January 2024).
  38. Brusola, P.; Garcia-Nieto, S. Fuzzy Modelling Toolbox for Aircraft Systems-Non-Linear TS Aircraft Model. 2024. Available online: https://github.com/sergarro/fuzzyAircraftModelling/blob/master/lib/dFuzzyModeldt.m (accessed on 30 January 2024).
  39. Jelonek, P.Z. Numerical Techniques in MATLAB: Differential Equations and Non-Linear Dynamics. 2018. Available online: https://warwick.ac.uk/fac/soc/economics/current/modules/rm/notes1/research_methods_matlab_3.pdf (accessed on 23 March 2024).
  40. Sala, A.; Ariño, C. Asymptotically necessary and sufficient conditions for stability and performance in fuzzy control: Applications of Polya’s theorem. Fuzzy Sets Syst. 2007, 158, 2671–2686. [Google Scholar] [CrossRef]
Figure 1. Aircraft axis reference [ x b , y b , z b ] , [ x v , y v , z v ] and [ x w , y w , z w ] , where α and β are the angle of attack and of sideslip, respectively, and V a i r is the airspeed.
Figure 1. Aircraft axis reference [ x b , y b , z b ] , [ x v , y v , z v ] and [ x w , y w , z w ] , where α and β are the angle of attack and of sideslip, respectively, and V a i r is the airspeed.
Aerospace 11 00258 g001
Figure 2. Global sector non-linearity and its membership functions.
Figure 2. Global sector non-linearity and its membership functions.
Aerospace 11 00258 g002
Figure 3. Local sector non-linearity.
Figure 3. Local sector non-linearity.
Aerospace 11 00258 g003
Figure 4. α ( t ) = tan 1 ( s ) global sector non-linearity.
Figure 4. α ( t ) = tan 1 ( s ) global sector non-linearity.
Aerospace 11 00258 g004
Figure 5. Membership functions M 1 and M 2 , Type I.
Figure 5. Membership functions M 1 and M 2 , Type I.
Aerospace 11 00258 g005
Figure 6. Membership functions N 1 and N 2 , Type IV.
Figure 6. Membership functions N 1 and N 2 , Type IV.
Aerospace 11 00258 g006
Figure 7. sin 1 x ( t ) global sector non-linearity.
Figure 7. sin 1 x ( t ) global sector non-linearity.
Aerospace 11 00258 g007
Figure 8. sin 1 x ( t ) membership functions.
Figure 8. sin 1 x ( t ) membership functions.
Aerospace 11 00258 g008
Figure 9. e x ( t ) global sector non-linearity.
Figure 9. e x ( t ) global sector non-linearity.
Aerospace 11 00258 g009
Figure 10. G z 2 1 and G z 2 2 membership functions.
Figure 10. G z 2 1 and G z 2 2 membership functions.
Aerospace 11 00258 g010
Figure 11. Flowchart which depicts the modeling process step by step.
Figure 11. Flowchart which depicts the modeling process step by step.
Aerospace 11 00258 g011
Figure 12. Simulated actuator values.
Figure 12. Simulated actuator values.
Aerospace 11 00258 g012
Figure 13. Validation of Ω and V.
Figure 13. Validation of Ω and V.
Aerospace 11 00258 g013
Figure 14. Validation of Φ and X.
Figure 14. Validation of Φ and X.
Aerospace 11 00258 g014
Table 1. α Fuzzy Model.
Table 1. α Fuzzy Model.
Model Rule 1Model Rule 2
IF  α ( t ) is “Zero” andIF  α ( t ) is “Zero” and
V a z ( t ) V a x ( t ) is “Positive” V a z ( t ) V a x ( t ) is “Negative”
THEN  α ( t ) = A α 1 THEN  α ( t ) = A α 2
Model Rule 3Model Rule 4
IF  α ( t ) is “Not zero” andIF  α ( t ) is “Not Zero” and
V a z ( t ) V a x ( t ) is “Positive” V a z ( t ) V a x ( t ) is “Negative”
THEN  α ( t ) = A α 3 THEN  α ( t ) = A α 4
Table 2. β ( t ) Fuzzy Model.
Table 2. β ( t ) Fuzzy Model.
Function:  a r c s i n ( V a y ( t ) V a ( t ) ) ,Membership function types: II and IV
Max and Min:  ( V a y V a ) m a x = 1 and ( V a y V a ) m i n = 1
Model Rule 1Model Rule 2
IF  β ( t ) is E β 1 and V a y ( t ) V a ( t ) is N β 1 IF  β ( t ) is E β 1 and V a y ( t ) V a ( t ) is N β 2
THEN  β ( t ) = A β 1 = π 2 THEN  β ( t ) = A β 2 = π 2
Model Rule 3Model Rule 4
IF  β ( t ) is E β 2 and V a y ( t ) V a ( t ) is N β 1 IF  β ( t ) is E β 2 and V a y ( t ) V a ( t ) is N β 2
THEN  β ( t ) = A β 3 = 1 THEN  β ( t ) = A β 4 = 1
Table 3. V a ( t ) Fuzzy Model.
Table 3. V a ( t ) Fuzzy Model.
Function:  V a x ( t ) 2 + V a y ( t ) 2 + V a z ( t ) 2 ,Membership function types: III
Max and Min:  ( V a x 2 + V a y 2 + V a z 2 ) m a x and ( V a x 2 + V a y 2 + V a z 2 ) m i n = 0
Model Rule 1Model Rule 2
IF  V a x ( t ) 2 + V a y ( t ) 2 + V a z ( t ) 2 is F V a 1 IF  V a x ( t ) 2 + V a y ( t ) 2 + V a z ( t ) 2 is F V a 2
THEN  V a ( t ) = A V a 1 = ( V a x 2 + V a y 2 + V a z 2 ) m a x THEN  V a ( t ) = A V a 2 = 0
Table 4. C L 1 Fuzzy Model.
Table 4. C L 1 Fuzzy Model.
Function:  C L 1 c ¯ V a ( t ) C L q q ( t ) ,Membership function types: IV
Max and Min:  ( q V a ) m a x and ( q V a ) m i n
Model Rule 1Model Rule 2
IF  q ( t ) V a ( t ) is N C L 1 1 IF  q ( t ) V a ( t ) is N C L 1 2
THEN  C L 1 = A C L 1 1 = c ¯ C L q ( q V a ) m a x THEN  C L 1 = A C L 1 2 = c ¯ C L q ( q V a ) m i n
Table 5. C L 2 Fuzzy Model.
Table 5. C L 2 Fuzzy Model.
Function:  C L 2 = C L H e λ L H L G ( t ) ,Membership function types: V
Max and Min:  λ L ( H L G ) m a x and λ L ( H L G ) m i n = 0
Model Rule 1Model Rule 2
IF  λ L H L G ( t ) is G C L 2 1 IF  λ L H L G ( t ) is G C L 2 2
THEN  C L 2 = A C L 2 1 = C L H ( λ L H L G ( t ) + 1 ) THEN  C L 2 = A C L 2 2 = 0
Table 6. C D 2 Fuzzy Model.
Table 6. C D 2 Fuzzy Model.
Function:  C D 2 = C D α 2 α ( t ) 2 ,Membership function types: VI
Max and Min:  α m a x and α m i n
Model Rule 1Model Rule 2
IF  α ( t ) is H C D 2 1 IF  α ( t ) is H C D 2 2
THEN  C D 2 = A C D 2 1 = C D α 2 α m a x α ( t ) THEN  C L 2 = A C L 2 2 = C D α 2 α m i n α ( t )
Table 7. C l 1 Fuzzy Model.
Table 7. C l 1 Fuzzy Model.
Function:  C l 1 = c ¯ V a ( t ) C l p p ( t ) ,Membership function types: IV
Max and Min:  ( p V a ) m a x and ( p V a ) m i n
Model Rule 1Model Rule 2
IF  p ( t ) V a ( t ) is N C l 1 1 IF  p ( t ) V a ( t ) is N C l 1 2
THEN  C l 1 = A C l 1 1 = c ¯ C l p ( p V a ) m a x THEN  C l 1 = A C L 1 2 = c ¯ C l p ( p V a ) m i n
Table 8. C l 2 Fuzzy Model.
Table 8. C l 2 Fuzzy Model.
Function:  C l 2 = c ¯ V a ( t ) r ( t ) ( C l r 0 + C l r α α ( t ) )
Membership function types: IV and VII
Max and Min:  ( r V a ) m a x and ( r V a ) m i n | α m a x and α m i n
Model Rule 1
IF  r ( t ) V a ( t ) is N C l 2 1 and α ( t ) is J C l 2 1
THEN  C l 2 = A C l 2 1 = ( r V a ) m a x ( C l r 0 + C l r α α m a x )
Model Rule 2
IF  r ( t ) V a ( t ) is N C l 2 1 and α ( t ) is J C l 2 2
THEN  C l 2 = A C l 2 2 = ( r V a ) m a x ( C l r 0 + C l r α α m i n )
Model Rule 3
IF  r ( t ) V a ( t ) is N C l 2 2 and α ( t ) is J C l 2 1
THEN  C l 2 = A C l 2 3 = ( r V a ) m i n ( C l r 0 + C l r α α m a x )
Model Rule 4
IF  r ( t ) V a ( t ) is N C l 2 2 and α ( t ) is J C l 2 2
THEN  C l 2 = A C l 2 4 = ( r V a ) m i n ( C l r 0 + C l r α α m i n )
Table 9. C m 1 Fuzzy Model.
Table 9. C m 1 Fuzzy Model.
Function:  C m 1 = c ¯ V a ( t ) C m q q ( t ) ,Membership function types: IV
Max and Min:  ( q V a ) m a x and ( q V a ) m i n
Model Rule 1Model Rule 2
IF  q ( t ) V a ( t ) is N C m 1 1 IF  q ( t ) V a ( t ) is N C m 1 2
THEN  C m 1 = A C m 1 1 = c ¯ C m q ( q V a ) m a x THEN  C m 1 = A C m 1 2 = c ¯ C m q ( q V a ) m i n
Table 10. C m 2 Fuzzy Model.
Table 10. C m 2 Fuzzy Model.
Function:  C m 2 = ( C m H 0 + C m H α α ( t ) ) e λ m H L G ( t )
Membership function types: V and VII
Max and Min:  λ m H L G , m a x and λ m H L G , m a x | α m a x and α m i n
Model Rule 1
IF  λ m H L G ( t ) is G C m 2 1 and α ( t ) is J C m 2 1
THEN  C m 2 = A C m 2 1 = ( λ m H L G ( t ) + 1 ) ( C m H 0 + C m H α α m a x )
Model Rule 2
IF  λ m H L G ( t ) is G C m 2 1 and α ( t ) is J C m 2 2
THEN  C m 2 = A C m 2 2 = ( λ m H L G ( t ) + 1 ) ( C m H 0 + C m H α α m i n )
Model Rule 3
IF  λ m H L G ( t ) is G C m 2 2 and α ( t ) is J C m 2 1
THEN  C m 2 = A C m 2 3 = 0
Model Rule 4
IF  λ m H L G ( t ) is G C m 2 2 and α ( t ) is J C m 2 2
THEN  C m 2 = A C m 2 4 = 0
Table 11. C n 1 Fuzzy Model.
Table 11. C n 1 Fuzzy Model.
Function:  c ¯ V a ( t ) C n r r ( t ) ,Membership function types: IV
Max and Min:  ( r V a ) m a x and ( r V a ) m i n
Model Rule 1Model Rule 2
IF  r ( t ) V a ( t ) is N C n 1 1 IF  r ( t ) V a ( t ) is N C n 1 2
THEN  C n 1 = A C n 1 1 = c ¯ C n r ( r V a ) m a x THEN  C n 1 = A C n 1 2 = c ¯ C n r ( r V a ) m i n
Table 12. C n 2 Fuzzy Model.
Table 12. C n 2 Fuzzy Model.
Function:  C n 2 = c ¯ V a ( t ) p ( t ) ( C n p 0 + C n p α α ( t ) )
Membership function types: IV and VII
Max and Min:  ( p V a ) m a x and ( p V a ) m i n | α m a x and α m i n
Model Rule 1
IF  p ( t ) V a ( t ) is N C n 2 1 and α ( t ) is J C n 2 1
THEN  C n 2 = A C n 2 1 = ( p V a ) m a x ( C n p 0 + C n p α α m a x )
Model Rule 2
IF  p ( t ) V a ( t ) is N C n 2 1 and α ( t ) is J C n 2 2
THEN  C n 2 = A C n 2 2 = ( p V a ) m a x ( C n p 0 + C n p α α m i n )
Model Rule 3
IF  p ( t ) V a ( t ) is N C n 2 2 and α ( t ) is J C n 2 1
THEN  C n 2 = A C n 2 3 = ( p V a ) m i n ( C n p 0 + C n p α α m a x )
Model Rule 4
IF  p ( t ) V a ( t ) is N C n 2 2 and α ( t ) is J C n 2 2
THEN  C n 2 = A C l 2 4 = ( p V a ) m i n ( C n p 0 + C n p α α m i n )
Table 13. C n 3 Fuzzy Model.
Table 13. C n 3 Fuzzy Model.
Function:  C n 3 = ( C n β 0 + C n β α α ( t ) ) β ( t )
Membership function types: IV and VII
Max and Min:  β m a x and β m i n | α m a x and α m i n
Model Rule 1
IF  β ( t ) is H C n 3 1 and α ( t ) is J C n 3 1
THEN  C n 3 = A C n 3 1 = β m a x ( C n β 0 + C n β α α m a x )
Model Rule 2
IF  β ( t ) is H C n 3 1 and α ( t ) is J C n 3 2
THEN  C n 3 = A C n 3 2 = β m a x ( C n β 0 + C n β α α m i n )
Model Rule 3
IF  β ( t ) is H C n 3 2 and α ( t ) is J C n 3 1
THEN  C n 3 = A C n 3 3 = β m i n ( C n β 0 + C n β α α m a x )
Model Rule 4
IF  β ( t ) is H C n 3 2 and α ( t ) is J C n 3 2
THEN  C n 3 = A C n 3 4 = β m i n ( C n β 0 + C n β α α m i n )
Table 14. Airbus A310 parameters.
Table 14. Airbus A310 parameters.
Mass and Geometry
S = 360 m2 c ¯  =  7.500 m z e n g  = 2 m m 0  =  150 , 000 kg I x x  =  10 7 kg/m2
I y y  =  1.600 · 10 7 kg/m2 I z z  =  2.400 · 10 7 kg/m2 I x z  =  10 6 kg/m2 x C G  = 21 m
Aerodynamic Coefficients
C L 0 = 0.900 C L α = 5.500 C L q = 3.300 C L δ e = 0.320 C L H = 0.200
λ L = 0.120 C D 0 = 0.065 C D α = 0.400 C D α 2 = 1.550 C Y β = 0.700
C Y δ r = 0.250 C l β = 3 C l p = 15 C l r 0 = 5 C l r α = 35
C δ a = 0.700 C l δ r = 0.200 C m 0 = 0.300 C m α = 1.500 C m q = 12
C m δ e = 1.200 C m h 0 = 0.090 C m h α = 0.900 λ m = 0.150 C n β 0 = 0.850
C n β α = 1.950 C n p 0 = 3 C n p α = 35 C n r = 7 C n δ a = 0.04
C n δ r = 1.250
Table 15. Engine and actuator characteristics.
Table 15. Engine and actuator characteristics.
ParameterTime-ConstantLower-BoundUpper-BoundRate-Limit
Engines (EPR)2 s0.9501.6000.100
Ailerons ( δ A )0.060 s−55 deg55 deg60 deg/s
Elevators ( δ E )0.070 s−25 deg25 deg20 deg/s
Rudder ( δ R )0.200 s−30 deg30 deg30 deg/s
Table 16. Fuzzy model limits.
Table 16. Fuzzy model limits.
ParameterLower-BoundUpper-BoundUnit
α π 2 π 2 rad
β π 2 π 2 rad
p V a 1 1rad/m
q V a 1 1rad/m
r V a 1 1rad/m
V a 2 1500m2/s2
V a z V a x 40 40(-)
V a y V a 1 1(-)
H L G 040m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Brusola, P.; Garcia-Nieto, S.; Salcedo, J.V.; Martinez, M.; Bishop, R.H. Fuzzy Modeling Framework Using Sector Non-Linearity Techniques for Fixed-Wing Aircrafts. Aerospace 2024, 11, 258. https://doi.org/10.3390/aerospace11040258

AMA Style

Brusola P, Garcia-Nieto S, Salcedo JV, Martinez M, Bishop RH. Fuzzy Modeling Framework Using Sector Non-Linearity Techniques for Fixed-Wing Aircrafts. Aerospace. 2024; 11(4):258. https://doi.org/10.3390/aerospace11040258

Chicago/Turabian Style

Brusola, Pablo, Sergio Garcia-Nieto, Jose Vicente Salcedo, Miguel Martinez, and Robert H. Bishop. 2024. "Fuzzy Modeling Framework Using Sector Non-Linearity Techniques for Fixed-Wing Aircrafts" Aerospace 11, no. 4: 258. https://doi.org/10.3390/aerospace11040258

APA Style

Brusola, P., Garcia-Nieto, S., Salcedo, J. V., Martinez, M., & Bishop, R. H. (2024). Fuzzy Modeling Framework Using Sector Non-Linearity Techniques for Fixed-Wing Aircrafts. Aerospace, 11(4), 258. https://doi.org/10.3390/aerospace11040258

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop