Fuzzy Logic-Based Aerodynamic Modeling with Continuous Differentiability

This paper presents a modeling method based on a fuzzy-logic algorithm to establish aerodynamic models by using the datasets from flight data recorder (FDR). The fuzzy-logic aerodynamic models are utilized to estimate more accurately the nonlinear unsteady aerodynamics for a transport aircraft, including the effects of atmospheric turbulence. The main objective in this paper is to present the model development and the resulting models with continuous differentiability. The uncertainty and correlation of the data points are estimated and improved bymonitoring amultivariable correlation coefficient in themodeling process.The latter is increased by applying a least square method to a set of data points to train a set of modeling coefficients. A commercial transport aircraft encountered severe atmospheric turbulence twice at transonic flight in descending phase is the study case in the present paper. The robustness and nonlinear interpolation capability of the fuzzy-logic algorithm are demonstrated in predicting the degradation in performance and stability characteristics of this transport in severe atmospheric turbulence with sudden plunging motion.


Introduction
During the operation of aircraft, it is very difficult to avoid the influence of adverse weather conditions.The effects of adverse atmospheric conditions usually are related to the dynamic aerodynamic effects that result from instantaneous changes of aircraft flight attitude.This situation, which cannot be easily flight-tested in the aircraft certification process, is the so-called abnormal flight condition [1].Such abnormal flight condition is very likely to cause the deployed aircraft out of control.In recent years, the transport aircraft response to the hazardous weathers due to atmospheric disturbance has been of great concern.The windshear, down burst, and atmospheric turbulence of hazardous weathers often lead to a high level of dynamic aerodynamic effects.These hazardous weathers usually have the characteristics of strong crosswind or vertical gust.It is not possible with the current technology to simulate the atmospheric turbulence effects by conventional flight simulators or wind tunnel on the aerodynamics, stability, and flying quality of aircraft [2].In addition, the aerodynamic derivative estimation from flight test data is always based on the concept of a small disturbance theory around an equilibrium state without having the influences of dynamic aerodynamic effects.
To provide the mitigation concepts and promote the understanding of aerodynamic responses of the transport aircraft under adverse weather conditions, a research project to evaluate the effects of weather hazards on performance and stability characteristics for transport aircraft is undertaken.
To effectively analyze the performance degradation and variations in static and dynamic stability of transport aircraft encountering hazardous weather, the nonlinear and dynamic (i.e., time-dependent) aerodynamic models based on flight data would be needed.
Traditional approaches of flight data analysis, such as the maximum likelihood method (MMLE) [3], the least-square or the regression method [4], apply the identification methods to data over an interval of time or other state variables to identify models with a quasisteady assumption over the interval.However, these conventional methods have not been shown to be capable of handling nonlinear and unsteady (i.e., dynamic) aerodynamic environments possibly exhibited in abnormal flight condition, in particular in hazardous weather, such as in-flight icing, windshear, and atmospheric turbulence.
In setting up the nonlinear and unsteady aerodynamic models to predict the required continuous derivatives of the function in aerodynamics for the analyses of stability and controllability, the traditional methods of using the tabulated data and the performance figures are very difficult to correlate complex functional relations among numerous parameters.
In earlier development of the fuzzy-logic algorithm, Zadeh [5] used the fuzzy sets to simulate physical parameters with membership functions.The disadvantage in this approach is that the predicted curves in the functional approximation tend to be piecewise continuous.Therefore, it was not suitable when derivatives of the function were needed in some applications of aviation technology.In 1985, Takagi and Sugeno [6] used the internal functions, instead of the fuzzy sets, in developing the fuzzy-logic algorithm.
It is well known that fuzzy models proposed by Takagi and Sugeno can be effectively represent unsteady and nonlinear system.Many important results in the applied areas, such as sensor fault detection [7], stabilization systems with timevarying delay [8], discrete-time state-delay systems [9], nonlinear itô stochastic systems with time-varying delay [10], a flexible air-breathing hypersonic vehicle [11], and fuzzy stochastic systems with time-varying delay [12] have great efforts in modern science and engineering, especially for some systems without having sufficient information as input databases for modeling [13,14].
In this approach, the main features are the internal functions, membership functions, and the output cells.The output curves in the prediction are smooth.Tan and Xie [15,16] applied the theory to simulate microelectronic processes with good accuracy.The present paper is based on the modeling technique suggested by Tan and Xie.This technique was first applied to aviation technology in 1997 and later to flight data [17][18][19][20][21][22][23][24].A new study to examine the effects of severe atmospheric turbulence on aircraft performance and stability in descending flight is presented in this paper.
Atmospheric turbulence has long been a difficult issue for the aviation community because it is not only the leading cause of serious injury in nonfatal accidents, but also most likely to cause the aircraft loss of control.Both the Federal Aviation Administration and National Aeronautics and Space Administration have research programs to deal with this problem.One of the objectives in NASA's Aviation Safety Program is to study the turbulence hazards through research, flight experiments, and data analysis.The hazard levels are quantified in terms of the peak values of the root mean square (RMS) of the normal loads with moving 5-sec intervals [25,26].The estimation was based on the assumption of von Karman's power spectral density [27].When von Karman's power spectral density is assumed, the Gaussian distribution is utilized to exhibit continuous turbulence.However, when localized severe turbulence or gust is included, the results will fall way off the correlation curve for the peak RMS normal loads to peak loads.Experience indicates that one main type of motion to cause flight injuries in atmospheric turbulence is the sudden plunging motion with the abrupt change in altitude.
A commercial transport aircraft encountered severe atmospheric turbulence twice in the descending phase.As a result, several passengers and cabin crews sustained injuries, because of which this event was classified as the aviation accident.Since one main type of motion to cause flight injuries is the sudden plunging motion with the abrupt change in altitude at transonic flight, so the main objective of this paper is to present the model development and the capability of resulting models to examine the static and dynamic stability characteristics in severe atmospheric turbulence with sudden plunging motion.

Theoretical Development
Since the nonlinear unsteady aerodynamic models are established by using flight data, modeling technique is important and need to be carefully considered.Factors that affect the modeling procedures include the mathematical tool to set up system model and the method to identify parameters of model structure.Modeling procedures start from separating the input data into many groups, and nonlinear relations are set up between each input-output data space.In order to obtain all the corresponding output, the present paper uses internal functions instead of fuzzy sets [5,28,29] to generate the output of the model.
The general idea of the present FLM technique is to set up the relations between its input and output variables of the whole system.The internal functions, membership functions, and outputs are three basic elements for the FLM approach.Two main tasks are involved in the FLM process.One is the identification of the coefficients of the internal functions, which is called parameter identification.The other one is structure identification to identify the optimal structure of fuzzy cells of the model.Details of the FLM technique are described in the followings.

Internal Functions.
The fuzzy-logic model uses many internal functions to cover the defined ranges of the influencing parameters.Although the form of internal functions is very simple, it can represent a highly nonlinear relationship between the input and output for the whole system.These internal functions are assumed to be linear functions of input variables [15,16] as follows: where    ,  = 0, 1, 2, . . ., , are the coefficients of internal functions   , and  is number of input variables.
The recorded data in QAR, such as flight altitude (ℎ), calibrated airspeed (CAS), and angle of attack (), is chosen as the input variables for a specific fuzzy model.In the present paper,   is denoted as an estimated aerodynamic coefficient of force or moment, and   are the variables of the input data.The numbers of the internal functions (i.e., cell's numbers) are quantified by the membership functions.

Membership Functions.
The values of each fuzzy variable, such as the angle of attack, are divided into several ranges, each of which represents a membership function with (  ) as its membership grade.One membership function from each variable constitutes a fuzzy cell.For the th cell, the corresponding membership grades are represented by    (  ),  = 1, 2, . . ., .In other words, the membership functions allow the membership grades of the internal functions for a given set of input variables to be calculated.For a given system with input variables  1 ,  2 , . . .,   , . . .,   , the recorded values of each input variables are normalized by using (  −  , min )/( , max −  , min ) to transform them into the ranges of [0, 1].The membership grading also ranges from 0 to 1; "0" meaning no effect from the corresponding internal function, and "1" meaning a full effect.Generally, overlapped straight lines, triangles, or trapezoids are frequently the shapes used to represent the functions.Because overlapped triangular membership function is simple and involves less computing time, it represents the grades of membership functions in the present FLM technique.
The total number of fuzzy cells is  =  1 ×  2 × ⋅ ⋅ ⋅ ×   × ⋅ ⋅ ⋅ ×   .For a variable   , the number of membership function is   .Each fuzzy cell is in a different combination from others formed by taking one membership function from each input variable.
The membership functions of triangular shapes are illustrated in Figure 1.In Figure 1, although the membership functions are continuous functions, there are discontinuities in slopes at some points.However, differentiation of membership functions is not performed in estimating derivatives.In the present application, aerodynamic derivatives are all estimated with a central difference scheme, which will be presented later.

Fuzzy Rule Inference.
A fuzzy cell is formed by taking one membership function from each variable.The total number of cells is the number of possible combinations by taking one membership function from each input variable.For every cell, it has a fuzzy rule to guide the input and output relations.The rule of the th cell [15,16] is stated as , and  2 is   2 ( 2 ) , and . . ., and   is    (  ) , then the cell output is where  = 1, 2, . . .,  is the index of the cells,  is the total number of cells of the model;   ( 1 ,  2 , . . .,   , . . .,   ) is the internal function with parameters   0 ,   1 , . . .,    , . . .,    to be determined, and    (  ) denotes the membership function for   .Each function covers a certain range of input variables.

Defuzzification.
In each fuzzy cell, the contribution to the outcome (i.e., the cell output) is based on the internal function, (5).The final prediction of the outcome is the weighted average of all cell outputs after the process of reasoning algorithm (i.e., fuzzy rule inference).Because of this weighting among many factors over large ranges of possibilities, the word "fuzzy" is derived to describe the method.However, its prediction after defuzzification is in a deterministic value (a "crisp" one).The output estimated by the fuzzy-logic algorithm corresponding to the th input ( 1, ,  2, , . . .,  , , . . .,  , ) is as follows: In ( 6) product [  ( 1, ), . . .,   ( , ), . . .,   ( , )] is the weighted value of the th cell and the index  of the data set, where  = 1, 2, . . ., , and  is the total number of the data record and the "product" stands product operator of its elements in this paper.

Parameter Identification.
Given a set of membership functions for each input variable, the unknown coefficients of the internal functions are determined by using the gradientdescent method.The gradient descent is also known as steepest descent.The accuracy of the established aerodynamic model through the fuzzy-logic algorithm is estimated by the sum of squared errors (SSEs) and the multiple correlation coefficients ( 2 ): In ( 7) and ( 8), where ŷ , the output of the fuzzy-logic model at point , is estimated by (6);   is the data point used for the model training at point ;  is the mean of the sample data, and  is the total number of data points.The model training is to determine the unknown coefficients of the internal functions,    , by maximizing the value of  2 .These coefficients are determined by the following iterative formula to minimize the sum of squared error (see (7)): where   is the convergence factor or the step size in the gradient method; subscript index  denotes the iteration sequence.
After simplification, (9) becomes the following.For  = 0, The present FLM technique was utilized to conduct NASA research projects [17,18].The unsteady aerodynamic models for F-16XL were based on dynamic wind-tunnel data in [17], and those were from flight test data in [18].The longitudinal and lateral-directional derivatives would be useful in aircraft stability analysis and control system design.Although the data source or application was not the same, the iteration during the search sequence stops of modeling should be the same with a principle for judgment.
The iteration during the search sequence stops when one of the following three criteria is satisfied: In the previos criteria, Cost = SSE  is the sum of squared errors (SSE) in current iteration to be denoted by "Cost" and RER = (cost current − cost previous)/cost current to be denoted by "RER" (i.e., the relative error) for simplicity in descriptions;  1 and  2 are the required precision criteria; and  max is a specified maximum iteration numbers.Given membership functions and the training data, this parameter identification procedure can be applied to establish a fuzzy-logic model.Although the procedure is well understood, to obtain a good model with the appropriate membership functions is quite challenging and requires numerical experimentation.

Model Structure Identification.
In the fuzzy-logic model, the model structure is indicated by the number of membership functions for each variable.For a fuzzy-logic model with multiple variables, the structure is the combination of the numbers and forms of the membership functions assigned to all input variables.Since the sequence defines the one-toone relationship between the numbers and the forms for each variable, the structure can be uniquely described by numbers of the membership functions.
The model structure is determined by maximizing the correlation coefficient, (8).A search forward algorithm has been employed for the identification.At each search stage, there may be many fuzzy-logic models with different structure combinations.The search stage numbers are denoted by   .Out of all the possible intermediate fuzzy-logic models at each search stage, for an efficient search, only some structures are developed and evaluated.Two selection criteria, to be given below, are used to choose these structures.With the incremental sequence and the selection criteria, the search forward algorithm is summarized as follows.
( max corresponding to all picking values is the optimal structure within a sensible   . The previous process is illustrated in Figure 2. The structure identification is tied up with parameter identification.In the structure identification, parameter identification to determine the -parameters according to ( 9) is also needed.The value of   is gradually changed in each iteration until the composite of   (i.e., optimal model structure) being fund, while the searching stage reaches  2 max , for the model structure identification; and the composite of   is fixed, the number of iteration is unfixed until both the values of  2 and RER reach the requirements in the final the parameter identification.

Application to Aircraft Aerodynamic Modeling
3.1.Flight Data.In the present paper, a commercial transport aircraft encountered sever atmospheric turbulence twice during descending phase.The first and the second atmospheric ( The necessary data in the FDR to determine the aerodynamics for this aircraft is time (), CAS, pressure altitude (ℎ), Euler angles (, , and ), the longitudinal, lateral, and vertical accelerations (  ,   ,   ), angle of attack (), aileron deflection (  ), elevator   , rudder   , stabilizer   , engine EPR, outside air temperature, wind speed (  ), wind direction, and fuel flow rate.Since only the normal acceleration is recorded in 8-Hz resolution (i.e., 8 points per second), all other parameters are interpolated with a monotone cubic spline to the same sampling rate.

Compatibility Analysis.
Typically, the longitudinal, lateral, and vertical accelerations (  ,   ,   ) along the (, , )body axes of aircraft, angle of attack , and the Euler angles (, , and ), as well as all control deflections are available and recorded in the FDR of all transport aircraft.Since the recorded flight data may contain errors (or called biases), compatibility analysis is performed to remove them by satisfying the following kinematic equations [19]: where  is the gravitational acceleration and  is the flight speed,  is sideslip angle,  is roll rate,  is pitch rate, and  is yaw rate in (13).Bias analysis is based on the least-square concept.The biases are:    ,    ,    ,   ,   ,   ,   ,   ,   ,   ,   ,   , respectively, for   ,   ,   , and so forth, representing random noise and systematic errors, and are estimated by minimizing the sum of squares of differences between the two sides of the kinematic equations through the following algorithm: where super fix " " stands for the mean value and the subscript "" indicates the measured or recorded values.The cost function is defined as where  is a weighting diagonal matrix with elements being 1.0 except the one for the slowly varying flight speed being 10.0 and ⇀  is calculated with a central difference scheme with  ⇀   , which is the measured value of  ⇀  .The steepest descent optimization method is adopted to minimize the cost function.As a result of the analysis, time derivatives, such as α , φ , and θ , are all estimated by using the central difference method; variables not present in the FDR, such as , , , , can be determined by cubic spline interpolation.
Note that since all flight variables recorded are based on the body axes, it is more convenient to estimate the force and moment coefficients for aircraft on the same axes system.Therefore the latter are obtained from the following flight dynamic equations [30]: where  is the aircraft mass;  is the dynamic pressure;  is the wing reference area;   ,   , and   are the longitudinal aerodynamic force and moment coefficients;   ,   , and   are the lateral-directional aerodynamic force and moment coefficients; and   ,   , and   are the moments of inertia about -, -, and -axes, respectively.The products of inertia,   ,   , and   , are assumed zero in the present case; but they are included in the equations because nonzero values may be available in other applications.The terms,   and   , represent the thrust contributions to the force in the direction of -axes, and to the pitching moment, respectively.

Equivalent Harmonic Motion.
The reduced frequency is a parameter to indicate the degree of unsteadiness in unsteady aerodynamics and is estimated in this paper by fitting the local trajectory with a harmonic motion.In the static case, the reduced frequency approaches 0. Large values of the reduced frequency imply the importance of unsteady aerodynamic effect.For longitudinal aerodynamics, the equivalent harmonic motion is the one based on the angle of attack variation (i.e., Theodorsen's theory).For lateral-directional aerodynamics, it is based on the time variation of roll angle [20].
For the longitudinal motion, the time history of the angle of attack () and time rate of angle of attack (d/d, or α ) are fitted with one of a harmonic motion at any instant as follows: where those terms on the left hand side of ( 18) are given and the unknowns are the local mean angle of attack (), the local amplitude of the harmonic motion (), the phase lag (), and the estimated harmonic frequency ( 1 ).These unknowns are calculated through an optimization method by minimizing the following cost function (least squares): In (19), where   is the measured value at point  and  is the number of the data points used in the optimization.For the case in the present study,  = 20 is found to be the best value.The 20 points preceding and including the current time are employed in (19).The least-square method is found to converge well and gives reasonably accurate results.The lateral-directional equivalent reduced frequency is computed in the same manner.
The local equivalent reduced frequency in the longitudinal motion is defined as where  is the mean chord length of wing airfoil section;  1 is the estimated harmonic frequency with time history of angle of attack.The lateral-directional equivalent reduced frequency is defined as where  is the wing span;  2 is the estimated harmonic frequency with time history of roll angle.

Fuzzy-Logic Aerodynamic Models.
Modeling means to establish the numerical relationship among certain variables of interest.In the fuzzy-logic model, more complete necessary influencing flight variables can be included to capture all possible effects on aircraft response to atmospheric disturbances [19,20].For longitudinal aerodynamics, the models are assumed to be of the form   ,   ,   =  (, α , ,  1 , ,   , , ,   , ) , where the left hand side represents the coefficients of axial force (  ), normal force (  ), and pitching moment (  ), respectively.All variables on the right hand side of ( 22) have been defined in the previous section.It should be noted that the stabilizer angle (  ) is included here, because it varies, though slowly, in flight to provide pitch trim (i.e., reducing the total static pitching moment to 0.0).The roll rate is included here because it is known that an aircraft under high aerodynamic loads at transonic speeds may have its longitudinal stability derivatives affected when additional disturbance due to roll rate is imposed.For the lateral-directional aerodynamics,   ,   ,   =  (, , , , ,  2 ,   ,   , , α , β ) , (23) where the left hand side represents the coefficients of side force (  ), rolling moment (  ), and yawing moment (  ), respectively.

Numerical Results and Discussions
The commercial transport aircraft encountered sever atmospheric turbulences in revenue flights in descending phase.
To examine the static and dynamic stability characteristics, it is imperative to understand the flight environment in detail.The dynamic aerodynamic effects can be expected to be very significant under the circumstances of instantaneous changes of , ℎ, and  in transonic flight.Since  reaches the value about 5.0 deg. at transonic flight, compressibility effect is important.It should be noted that the turbulent vertical wind field was not measured or estimated in the FDR; but it is included in the total .The pitch angle () does not vary as much as the angle of attack (), reaching only 2.5 deg.The variations of  are large, especially during the second turbulence encounter.The magnitudes of pitch rate () and roll rate () are shown in Figure 4(c); the variation of  is larger than that of .The yaw rate is not shown because it is small throughout.

Pilot's Responses in Operations.
The flight data recorder indicated that the autopilot was deactivated from  = 7489 sec to  = 7493 sec, for 4 seconds.It is of interest to examine the pilot's response in operations after the autopilot being tripped off.The time history of   and   is shown in Figure 4(d).The magnitude of   becomes positive to keep pitch down in descending.The variation of   in the second turbulence encounter is larger than that of the first one.While the roll angle () is changing fast in the negative direction right before  = 7643 sec,   is increasingly more positive, reaching +13 deg at  = 7645 sec, as shown in Figure 4(d).On the other hand, opposite changes occur around  = 7646 sec.The   is not shown because it is small throughout.
Both the pitch and roll rates are significant in severe atmospheric turbulence, especially the latter.The yaw rate is not shown in Figure 4, because it is small throughout.In normal flight conditions, pitch and roll rates are typically small, in particular at transonic speeds, to avoid nonlinear Figure 5: The main aerodynamic coefficients predicted by the nonlinear and unsteady aerodynamic models.
aerodynamic effects.The autopilot was tripped off in severe atmospheric turbulence.The pitch angle and roll angle are the only cues to the pilots in making the control decision.

Analysis of Model Predictions.
In the present study, the accuracy of the established unsteady aerodynamic models with six aerodynamic coefficients by using FLM technique is estimated by the sum of squared errors (SSE) and the square of multiple correlation coefficients ( 2 ).All the aerodynamic derivatives in the study of stability and controllability are calculated with these aerodynamic models of aerodynamic coefficients.
Figure 5 presents the aerodynamic coefficients of normal force   , pitching moment   , rolling moment   , and yawing moment   predicted by the unsteady aerodynamic models.The predicted data by the final models has a good match with the flight data.
The final main aerodynamic models of aerodynamic coefficients consist of many fuzzy rules for each coefficient as described in Tables 1 and 2. In Tables 1 and 2, the numbers below each input variable represents the number of membership function.The total number of fuzzy cells () in each model is the product of each number which presented in column 3. The last column shows the final multiple correlation coefficients ( 2 ).The accuracy of the established aerodynamic model through the fuzzy-logic algorithm can be judged by the multiple correlation coefficients ( 2 ).

Analysis of Stability Characteristics.
Since the sudden plunging motion with the abrupt change in altitude affects the flight safety the most in severe atmospheric turbulence, so the time period of 7482 ∼ 7494 sec is emphasized in examining the characteristics of stability.In order to examine the variations in characteristics, the units of all aerodynamic derivatives are converted to rad −1 .It should be noted that these derivatives are evaluated at the instantaneous conditions, instead of about the trim conditions as have been traditionally done.The time history of main longitudinal and lateral-directional static stability derivatives along the flight path is presented in Figure 6.From the point of view in static stability, initially, the configuration has longitudinal stability (  > 0 and   < 0) in Figure 6(a), stable longitudinal damping (  < 0) in Figure 6(b), lateral stability (  < 0) and directional stability (  > 0) in Figure 6(c), roll damping (  < 0), and directional damping (  < 0) in Figure 6(d).During the plunging motion,   > 0 and   > 0 in the period of  = 7483.0∼ 7484.5 sec, and the static stability becomes unstable with insufficient roll damping (  ≈ 0) and directional damping (  > 0).The aerodynamic instability is most likely caused by the motion that produces a timedependent pressure distribution on the aircraft surface involving compressibility effects.
Figure 7 presents the time history of main longitudinal and lateral-directional oscillatory derivatives along the flight path to associate with α and β -derivatives.Note that in Figure 7(a), the oscillatory derivatives are defined as In Figure 7(c), the oscillatory derivatives are defined as The values of oscillatory derivatives are equivalent to the combinations of steady damping and dynamic derivatives in the previous equations ( 24)∼ (25).The use of oscillatory derivative instead of steady damping ones is more consistent with the actual case of the aerodynamic damping in the analysis of stability characteristics.To be stable, (  ) osc > 0, (  ) osc < 0, (  ) osc < 0, and (  ) osc < 0. Physically, if it is unstable, the motion will be divergent in oscillatory motions.The values in the period of plunging motion have some differences between oscillatory and steady damping derivatives in Figures 7(a) and 7(c) due to the effects of the dynamic derivatives (i.e., α -and β -derivatives).Figures 7(b) and 7(d) show dynamic derivatives of stability.To be stable,   α > 0,   α < 0,   β < 0, and   β > 0. The magnitudes of   α and   α have significant variations and   α < 0 in the period of  = 7482.5∼ 7491 sec in Figure 7(b).It should be noted that   α represents the virtual mass effect and is particularly large in transonic flow to affect the plunging motion [31].The (  ) osc and (  ) osc are inadequate in oscillatory damping in the periods of  = 7484 ∼ 7486 sec and  = 7490 ∼ 7494 sec, as shown in Figure 7(a).The effect of α -derivative on (  ) osc is to improve the stability in pitch in the periods of  = 7482.5∼ 7484 sec and 7486 ∼ 7488.5 sec.The magnitudes of   β in Figure 7 ) ) In essence, the effects of α -derivative on (  ) osc and β derivative on (  ) osc are small.However, the effect of αderivative on (  ) osc is to improve the stability in pitch; while the effects of β -derivative are to cause more directional instability.These results indicate that the turbulent crosswind has some adverse effects on directional stability and damping.Although the dynamic derivatives tend to be small for the present configuration, this is much helpful to understand the unknown factors of instability characteristics.

Concluding Remarks
The main objective in this paper was to present the development of nonlinear and unsteady aerodynamic models based on FLM method and the resulting models having the capability to create the required continuous derivatives for the investigative analysis.The static and dynamic stability characteristics in severe atmospheric turbulence with sudden plunging motion in descending flight were examined by using those continuous derivatives.The numerical results and discussions of application to aerodynamic modeling were concluded as follows.
(1) The uncertainty and correlation of the data points were estimated and improved by monitoring a multivariable correlation coefficient in the modeling process.
(2) The predicted results showed that the final models could not only produce relatively accurate aerodynamic coefficients but also generate smooth coefficient curves with the required independent variables to provide reasonable local aerodynamic derivatives.
(3) The fuzzy logic-based aerodynamic modeling method was shown to be capable of handling nonlinear and unsteady aerodynamic environments possible exhibited in abnormal flight condition.
The analytical results could provide the mitigation concepts and promote the understanding of aerodynamic responses of the commercial transport aircraft in severe atmospheric turbulence.

Figure 3 :
Figure 3: The time history of flight variables in severe atmospheric turbulence during the descending phase in transonic flight.

4. 1 .
Aerodynamic and Flight Environments.The corresponding flight data for this aircraft in severe atmospheric turbulence during descending flight is presented in Figure3.The dataset of time span from  = 7480 ∼ 7739 sec used for the modeling is extracted from the FDR.This aircraft encounters severe atmospheric turbulence twice during this time span.In Figure3(a), the variations of normal acceleration (  ) show the highest   being 2.05 g around  = 7483 sec and the lowest being −1.05 g around  = 7484 sec in the first turbulence encounter, the highest   being 1.91 g around  = 7682 sec and the lowest being −0.16 g around  = 7684 sec in the second one.Figure3(b) shows that variation of  is approximately  (deg)

Figure 4 :
Figure 4: The variations of dynamic characteristics and control variables in severe atmospheric turbulence during the descending phase in transonic flight.
Figure 4 presents the variations of dynamic characteristics and control variables.The variation ranges of  are 4 deg.to −1.8 deg. in the first turbulence encounter and 5.5 deg. to −2 deg. in the second one; the time history of  is about −2 deg.with the magnitude of small fluctuation, as indicated in Figure 4(a).The time history of pitch angle () and roll angle () is shown in Figure 4(b).

Figure 7 :
Figure 7: The time history of main longitudinal and lateral-directional oscillatory derivatives along the flight path.

Table 1 :
The final main aerodynamic models of longitudinal aerodynamics.

Table 2 :
The final main aerodynamic models of lateral-directional aerodynamics.
(  ):Membership function for input variable     ,   ,   : Acceleration components along the (, , )-body axes of aircraft, respectively,  : W i n gs p a n ,m CAS: Calibrated airspeed, knows, m/sec   ,   ,   : Longitudinal aerodynamic force and moment coefficients   ,   ,   : Lateral-directional aerodynamic force and moment coefficients : M e a na e r o d y n a m i cc h o r d ,m ℎ: A l t i t u d e ,m   ,   ,   : Moments of inertia about the -, -, and axes, respectively, kg⋅m 2   ,   ,   : Products of inertia, kg⋅m 2  1 ,  2 : Sideslip angle, deg and time rate of sideslip angle, deg/sec., respectively   ,   ,   : Control deflection angles of aileron, elevator, and rudder, respectively, deg   : Stabilizer angle, deg , , : Euler angles in roll, pitch, and yaw, respectively, deg  1 ,  2 : Estimated harmonic frequencies with time histories of angle of attack and roll angle, respectively, sec −1 .