New methodology for the prediction of the

3 downloads 0 Views 608KB Size Report
compared with the XFoil results for different flight cases, expressed as various ..... The SVM-EGD approach was implemented in MATLAB. The results for the CL, ...
New methodology for the prediction of the aerodynamic coefficients of an ATR-42 scaled wing model Abdallah BEN MOSBAH, Ruxandra BOTEZ, Thien My DAO École de technologie supérieure (ÉTS), LARCASE, www.larcase.etsmtl.ca, Montreal, Quebec, H3C1K, Canada

l

Abstract A new approach for the prediction of lift, drag and moment coefficients is presented. This approach is based on the Support Vector Machines methodology, and on a new optimization algorithm, theExtended Great Deluge. The training and validation of this new combination method is realized using the aerodynamic coefficients of an ATR-42 wing model with Xfoil software. The results obtained with our approach are compared with the XFoil results for different flight cases, expressed as various combinations of angles of attack and Mach numbers. The main purpose of this methodology is to rapidly predict aircraft aerodynamic coefficients.

Introduction Supervised learning is a technique in which rules are automatically generated from a database. This learning database is characterized by a pair of inputs-outputs (xn,yn), where yn= f(xn). The objective of the supervised learning method is to determine a representation of the function f, called ''prediction function h''. This new function h provides an output y’=h(x’) for a new input x’. There are two types of problems that we seek to solve by means of supervised learning. The first is a ''regression problem'', in which the output associated with an input is a real number. The second type is a ''classification problem'', in which the output has a finite cardinal and a label should be assigned to a given input [1]. A regression problem is presented in this study, in which the desired outputs are real numbers. There are many supervised learning methods, such as neural networks (NN) and support vector machines (SVM), Page 1 of 8

which is used in this work. SVM, NN and fuzzy logic can be applied to design prediction or identification models. NN and fuzzy logic have been used extensively on control systems in aerospace. Neural networks have been used in multiple domains, including pattern classification, optimal control, and manufacturing [2–6]. In aerospace engineering, NNs can be applied to a large range of complex problems, as in [7]. Neural network have been used to resolve many problems in the aeronautic industry, such as: the detection and identification of structural damage [8], the modeling of aerodynamic characteristics from flight data [9, 10], the detection of unanticipated effects such as icing [11,12,13], and autopilot controllers and advanced control laws for applications such as carefree maneuvering [14, 15], as well as CL and CD aerodynamic coefficient prediction [16]. An experimental study on the use of smart sensing and neural networks to strain loads for different airflow cases are presented by Lunia et al. [17], in which the authors used a fiber-optic sensor to train and verify the neural network performance. They used a multilayer network to map the inputs and the outputs of a nonlinear system, thereby creating a three-layered neural network, with three neurons in the first layer, five in the second, and one neuron in the third. An hyperbolic tangent sigmoid function is used in the first and the second layer; and a linear function in the third layer. The network was developed through supervised learning, in which parameters are adjusted to achieve the target outputs for given inputs. Scott [18] developed an adaptive neural network-based control system that integrated three control systems that were developed and tested. One system used flutter suppression control laws, a second system employed a predictive NN control scheme, and the third system used an NN in an inverse model control scheme. Suresh et al. [19] used recurrent neural networks for the prediction of lift coefficients at high angle of attack. In their

approach, the lift coefficient was obtained from wind tunnel tests. Fei et al. [20] evaluated the air speed, the angle of attack and the angle of sideslip fundamental parameters in the control of flying bodies. They have proposed a new experimental methodology by which the flight parameters are inferred from multiple hot-film flow speed sensors mounted on the surface of the wing of a Micro Air Vehicle (MAV). In order to obtain a good mathematical relationship between the sensor readings and the flight parameters, they propose using micro hot-film flow speed sensor arrays and a back-propagation neural network to determine three flight parameters: air speed, angle of attack and angle of sideslip [20]. Peyada [21] proposes using Feed Forward Neural Networks to estimate aircraft parameters from flight data. That method uses Feed Forward Neural Networks to establish a neural model that can be used th to predict the time histories of motion variables at the (k + 1) instant, given that the measured initial conditions correspond to th the k instant [21]. A neural network based on a flush air data sensing system and demonstrated on a mini air vehicle was presented by Samy et al. [22]. Xuan et al. [23] presented a fuzzy neural network controller, which had the advantages of both fuzzy control and neural network, for the control of uncertain parameters for nonlinear time-varying systems. Fuzzy logic can be used to model highly non-linear, multidimensional systems, including those with parameter variations, or where the sensors’ signals are not accurate enough for other models [24]. De Jesus-Mota et Botez [25] proposed a new technique for helicopter model identification from flight data tests based on neural networks. The dynamic behavior of the helicopter was identified with a recurrence method and an optimization procedure based on neural network theory and tuning of the initial conditions [25]. A aeroservoelastic model was presente by o ly et al. , ], in which neural network and fuzzy logic algorithms identify the multi-input and the multi-output system of the F/A-18 aircraft. A flight parameter control system based on neural networks has been proposed by Mosbah et al. [28], in which the proposed NNs are optimized using a metaheuristic algorithm, called the extended great deluge (EDG). Their approach predicts pressure distributions and aerodynamic coefficients from the known parameters (angle of attack, Mach number, etc.). Kouba et al. [29,30] proposed an identification model, based on fuzzy logic methods, to identify the nonlinear aircraft models for many flight-test cases; their proposed model used an F/A-18 aircraft. A new method for the realization of two neuro-fuzzy controllers for a morphing wing design application is presented by Grigorie et al. [31]. The propose controllers’ main function is to correlate each set of pressure differences, calculated between the optimized and the reference airfoil, with each of the airfoil deformations produce by the actuators’ system [31]. Several other uses of fuzzy logic are presented by the same authors in identification and control areas in [3237].

Support Vector Machines (SVM) The SVM method is used to solve regression and classification problems. The idea is to estimate a prediction function h from the original function f given in the learning phase. y= f(x) Page 2 of 8

y’=h(x’) SVM method

This estimation is made from a training set of N samples. A tube of wi th ε is efine aroun the esire outputs, so that all the predicted values should be inside this tube. However, Smola et al. [38] selected their prediction function to be flat and smooth so that it could be approximated with a linear function.

(1) where ω is weight of the inputs space an b is a billet Є R In order to satisfy these two conditions, the input space is converted to a larger size space, in which the transfer function for this conversion is enote by Φ, as seen in equation (1). As described by Vapnik [39] , the optimization of the prediction function h(x) in the larger space is characterized by the resolution of the following system [1]:

min (,  ) 

N 1 2   C  ( i   i* ) 2 i 1

(2)

subject to:

((.( xi ))  b)  yii     i*

 i ,  i*  0

i  1...N

where:  C is a regularization parameter of the compromise between the complexity of the model, for example the planarity of the function f, and the degree of tolerance for the deviation of the samples against ε; an 

 i ,  i*

are variables of samples which are outside of the

ε-tube (Figure 1).

Figure 1. Illustration of the principle of slack variables ζi and ζi* [49]

An optimization problem has a dual form in which the objective function and the constraints are ''strictly convex''. If the conditions are satisfied, the ''dual problem'' becomes equivalent to the ''original problem''. For this reason, we use a Lagrangian function, where the Lagrangian is defined as the sum of the objective function and a linear combination of constraints whose coefficients (αi ≥ 0) are calle Lagrange multipliers [1]. As explained by Smola et al. [38], by introducing Lagrange multipliers, our optimization problem defined in equation (2)

becomes a dual form, and the shape of the ''regression function'' is defined as follows: ∑

(3)

where:  αi are the Lagrange multipliers;  M is the subset of samples corresponding to the nonzero Lagrange multiplier; and  K is the kernel function representing a scalar product in the re-description space. The kernel function is used to construct the decision surface (hyper-plan) in the input space. Some examples of kernel functions are presented in [1]: -Linear kernel function: -Polynomial kernel function: -Gaussian kernel function:





)

Optimization of the SVM parameters To obtain a good regression it is necessary to use the appropriate optimized SVM parameters, such as C, ε and K. In fact the parameters C and ε are optimized. It is also important to use the correct type of kernel function K and to optimize its parameters, such as its degree d, for example, in the case of its polynomial form. To optimize these three parameters, many techniques have been proposed. For example, Keerthi and al.[40] used a technique to minimize the margin. Cherkasky and al.[41] proposed an analytical method to optimize the parameters C and ε when a Gaussian kernel function wasselected. The meta-heuristic methods were also used to optimize the SVM parameters, while Ping and al. [42] used a simulated annealing algorithm to fin the values of C an ε parameters. Since the quality of results necessarily depends on the quality of the SVM parameters, we propose an original hybridization of the SVM method with the ''Extended Great Deluge'' (EGD) algorithm to optimize these parameters. The EGD algorithm is described.

Extended Great Deluge algorithm The local search procedure called "Extended Great Deluge" (EGD) was introduced by Dueck [48] in 1993. It considers a local search algorithm in which bad solutions might be obtained, whose values are smaller than a certain limit B. This limit B decreases monotonically (in the case of minimization problems) during the search. The initial value of B is equal to the "objective function", and for each iteration its value ecreases by a fixe Δ in minimization problems, an increases by the same value of Δ for maximization problems. The step of Δ represents an ''input parameter'' in this approach. During the search, B is the limit between a feasible and a not-feasible area of research, it serves to orientate the solution of the problem towards the ''feasible area''. In other words, the neighborhood of the solution s* is cut by the limit B and the research is only to be conducted on one side, below or above the limit B, depending on the objective function Page 3 of 8

(minimization or maximization of the objective function). The increase or ecrease of with Δ can be consi ere as a control process giving the desired solution. In the beginning stages of research, the actual solution has the ability to move in both directions, and can also be found inside the feasible portion limited by B. Otherwise, there is a great chance of accepting poorsolutions, because the limit B is located at a long distance from the chosen solution s* and a small part of its neighborhood may be cut off. During the search, the limit B moves closer to the value of the current solution, the search space becomes smaller and the possibility of improving the solution becomes lower, leading to the end of the process research [7]. The first application realized with this approach was the optimization of the exam timetable problem treated by Burke et al. [43]. The results proved the effectiveness of this algorithm [43]. Several of these results have been improved using approaches such as ''taboos search''. Figure 2 shows the steps of the EGD algorithm [44]. The second application of the EGD algorithm was the optimization of the group scheduling problem by Mosbah [45, 46, 47, 48] (first author of this paper), in which the EGD algorithm gave better results than genetic algorithm (GA) and simulated annealing (SA). The advantage of the EGD is that only one parameter needs to be a juste , Δ , the rate of increase or decrease of B. Burke et al. [43] have shown that the convergence time of the algorithm epen e on the Δ value. In ee , increase of the Δ value will result in a ecrease of the convergence time, but the quality of solutions can also degrade, which explains the importance of choosing the best Δ to obtain a goo compromise between quality of results and calculation time. The EDG algorithm is a metaheuristic type in which the optimum solution is not guaranteed, because the search process is based on a randomly selected initial solution, which is the initial boundary B and the value of ΔB [48]. The first step of the algorithm, as shown in Figure 2, consists of ran omly initializing Δ an the initial value of the SVM parameters (S). Next, the efficiency α (the Error) of S is calculated, and its value is assigned to B. The neighborhood * N(S) of S is then defined; a neighboring solution S belonging * to N(S) is randomly selected; and the new solution N(S ) is compared with the old solution, N(S) and B. If two conditions * (α(S )≤ α(S) and B) are not satisfied, then a new neighboring solution S is selected; if one of these conditions is true, then * the solution S is accepted and B=B-ΔB is recalculated. Finally, the stopping criterion is tested, and if the number of iterations has been reached, the algorithm ends. A new neighboring solution S is selected and the algorithm is tested again [47].

New proposed SVM-EGD algorithm In this work, a hybrid SVM- EGD method is proposed to control the lift, drag and moment coefficients for different combinations of angles of attack and Mach number values. Figure 3 shows the steps used to obtain the optimal SVM parameters so that the error between the calculated coefficients using XFoil and the calculated coefficients using SVM-EGD remains as small as possible.

In our new algorithm, we use a qualitative performance measure describing the learning abilities of a given trained SVM method, that is, the training error is expressed as the mean sum of the squared residuals (errors) in the training data [40]:



(4)

where t is the time when the y is calculated, y(xk) is the th estimated output of the SVM method for the k input xk and N is the number of data points used in the training set. The performance and results obtained with our proposed algorithm are further evaluated to determine the lift, drag and moment coefficients for varying angles of attack and Mach numbers using Xfoil values.

start

Initialize :∆ , S

Calculate the efficiency of S (s) B(S)

Define the neighborhood N (S)

Select randomly a neighboring solution Figure 2. The Extended Great Deluge [44] ok

no (S*) ≤(S) or (S*) ≤ B

Accept S*and put S=S*

B= B-∆ Stopping criterion

ok END

no

Page 4 of 8

Price-Paidoussis wind tunnel

Start

The validation of this new approach will be achieved by experimental tests using Price-Paidoussis the wind tunnel of LARCASE,shown in Figure 4. This wind tunnel has two test 2 chambers; a first chamber with a section equal to 0.3 x 0.6 m that provides a speed up to 60 m/s, and a second chamber test 2 with a section of 0.6 x 0.9 m that provides a speed up to 40 m/s [50].

Initialize SVM parameters, Δ , number of iteration and error_int

learning vectors :inputs x and outputs y B

error_int

Define the neighborhood N

Select randomly a neighboring solution S*

Learning

Output vector y' Figure 4.

Price-Paidoussis wind tunnel

Error (y,y')

no

Implementation of the SVM-EGD algorithm and preliminary results Error(y,y')≤error(y,y'-1) Or

The objective is to determine the aerodynamic lift (Cl), drag (Cd) and moment (Cm) coefficients for different values of Mach numbers and angles of attack on the ATR-42 wing.

Error(y') < B

ok

For the ''learning'' of the SVM, a database obtained using Xfoil software was used. A total of 101 values of lift (Cl), drag

Accept solution S=S*

Iteration number?

no

ok

SVM-parameters optimized and selected

Predicted outputs y_pred with SVM_EGD

New input data

Figure 3.

Page 5 of 8

Proposed algorithm

B=B- Δ

(Cd) and moment (Cm) coefficients for combinations of angles o o o of attack α between -5 to 5 (0.1 per step) and a Mach -1 number = 0.11 (40 m s ) were used. These test cases were selected so that they could be validated using the PraicePaidoussis wind tunnel. The validation data set was composed of 11 random vectors and the test data set had had 11 random vectors. The optimal values of the SVM parameters (kernel, C, degree d and Ɛ) were obtained using the EGD algorithm. The best results were obtained using the Gaussian kernel with the parameters C= 8036518, d=1 and Ɛ=9.29 10-6. The SVM-EGD approach was implemented in MATLAB. The results for the CL, CD and CM are values are presented in Tables 1 to 3 and Figures 5 to 7.

Table 1. Comparison for lift coefficients CL (original vs prediction) for different airflow cases

0,11757

Original Lift coefficient -0,2676

Predicted Lift coefficient -0,2708196

0,11768

-0,1562

-0,15527248

0,594

-2,7

0,11777

-0,0654

-0,06557326

-0,265 -7,931

α

Mach number

-4,7 -3,6

Error (%) -1,203

-2

0,11784

0,0078

0,00841861

-1,3

0,11791

0,0814

0,07942546

2,426

-0,1

0,11803

0,198

0,19817848

-0,090

0,6

0,1181

0,2648

0,26317644

0,613

1,6

0,1182

0,364

0,36395321

0,013

2,6

0,1183

0,4684

0,46793551

0,099

3,9

0,11843

0,6176

0,61415461

0,558

4,7

0,11851

0,7219

0,72523071

-0,461

Figure 6.

Drag coefficient CD versus angle of attack α

Table 3. Comparison of CM coefficients (original vs prediction) for different airflow cases

Figure 5. Lift coefficient CL versus angle of attack α Table 2. Comparison of CD coefficients (original vs prediction) for different airflow cases α

Mach number

Original Drag coefficient

-4,7

0,11757

0,01101

Predicted Drag coefficient 0,01103101

-3,6

0,11768

0,01082

0,01083617

-0,149

-2,7

0,11777

0,01082

0,01083611

-0,149

-2

0,11784

0,01065

0,010641

0,084

-1,3

0,11791

0,01041

0,01045531

-0,435

-0,1

0,11803

0,00925

0,00922795

0,238

0,6

0,1181

0,00877

0,00875986

0,116

1,6

0,1182

0,00863

0,00865606

-0,302

2,6

0,1183

0,00905

0,00904081

0,102

3,9

0,11843

0,00973

0,00970529

0,254

4,7

0,11851

0,01018

0,01020146

-0,211

0,11757

Original Moment coefficient -0,0342

Predicted Moment coefficient -0,03427513

0,11768

-0,031

-0,03098787

0,039

-2,7

0,11777

-0,0286

-0,02856293

0,130

-2

0,11784

-0,0271

-0,02708922

0,040

-1,3

0,11791

-0,0258

-0,02585505

-0,213

-0,1

0,11803

-0,0221

-0,02206207

0,172

0,6

0,1181

-0,0193

-0,0192996

0,002

1,6

0,1182

-0,0155

-0,0155388

-0,250

2,6

0,1183

-0,0128

-0,01278234

0,138

3,9

0,11843

-0,0128

-0,01249511

2,382

4,7

0,11851

-0,0156

-0,01555607

0,282

α

Mach number

-4,7 -3,6

-0,220

Error (%) -0,191

Figure 7.

Page 6 of 8

Error %

Moment coefficient CM versus angle of attack α

Conclusions Anew algorithm using the SVM-EGD approach was used to optimize the values of SVM parameters. The obtained results were acceptable and can be further improved using experimental wind tunnel test validation. The average percentage error did not exceed a maximum of 1.3 % for the lift coefficient, and most of the error values were lower than 0.2% for the drag and 0.35% for the moment coefficients. Based on these results, the SVM-EGD approach was found to be robust and accurate. Future work will include verifying the performance of the proposed approach with experimental wind tunnel tests.

References El Asli, N., ‘‘Approche hybri e basée sur les machines à vecteurs de support et les algorithmes génétiques pour l’estimation es coûts e fabrication,’’ Master thesis, École de technologie supérieure, Quebec University, 2008. 2. Wrong, B. K., Vincent, S. L., and Lam J., "A bibliography of neural network business applications research: 1994– 1998," Computers and Operation Research, 27:1045– 1076, 2002. 3. Hunt, K. J., Sbarbaro, D., Zbikowski, R., and Gawthrop, P. J., "Neural networks for control systems – a survey," Automatica, 28:1083–1112,1992. 4. Udo, G. J., "Neural networks application in manufacturing process," Computers and Industrial Engineering, 23 (1– 4):97–100, 1992. 5. Wong, B. K., Bodnovich, T. A., and Selvi, Y., "Neural network application in business: a review and analysis of the literature (1988–1995)," Decision Support Systems , 19:301–320, 1997. 6. Chen, D., and Burrel, P., "On the optimal structure design of multilayer feedforward neural networks for pattern recognition," Int. J. Pattern Recognition and Artif. Intell. 6 (4):375–398, 2002. 7. Faller, W. E., and Schreck, S. J., "Neural networks: Applications and opportunities in aeronautics," Progress in Aerospace Sciences, 32:433–456,1996. 8. Rauch, H. E., Kline-Schoder, R. J., Adams, J. C., and Youssef, H. M., "Fault detection, isolation, and reconfiguration for aircraft using neural networks," AIAA Paper:1993-3870, 1993. 9. Linse, D. J. , and Stengel, R. F., "Identification of aerodynamic coefficients using computational neural networks," Journal of Guidance, Control and Dynamics, 16(6):1018–1025, 1993. 10. Amin, S. M., Gerhart, V., and Rodin, E.Y., "System identification via artificial neural networks: applications to online aircraft parameter estimation," AIAA Paper:975612, 1997. 11. Johnson, M. D. , and Rokhsak, K. "Using artificial neural network and self-organizing maps for detection of airframe icing,"Journal of Aircraft, 38(2):224–230, 2001. 12. Aykan, R., ''Kalman filter and neural network-based icing identification applied to A340 aircraft dynamics,'' Aircraft Engineering and Aerospace Technology: An International Journal , 77 (1): 23–33, 2005. 1.

Page 7 of 8

13. Johnson, M.D., and Rokhsaz, K., “Using artificial Neural networks and self organizing maps for detection of airframe icing,” 000 Atmospheric Flight Mechanics Conference, AIAA-2000-4099., 2000. 14. Napolitano, M. R., and Kincheloe, M., "On-line learning neural network controllers for autopilot systems," Journal of Guidance, Control and Dynamics, 33(6):1008–1015, 1995. 15. Yavrucuk, I., Prasad, J.V. R., Calise, A., "Adaptive limit detection and avoidance for carefree manoeuvring," AIAA Paper: 2001-4003, 2001. 16. Wallach, R., De Mattos, B. S., and Da Mota Girardi, R., '' Aerodynamic coefficient prediction of a general transport aircraft using neural network,'' 25th International Congress of the Aeronautical Sciences ICAS 2007. 17. Lunia, A., Isaac, K. M., Chandrashekhara, K., and Watkins, S.E., "Aerodynamic testing of a smart composite wing using fiber-optic strain sensing and neural networks,” Smart MaterIals and Structures, 9:767–773, 2000. 18. Scott, R. C., ‘‘Active control of win -tunnel model aeroelastic response using neural networks’’, Journal of Guidance, Control, and Dynamics, 23(6):1100-1108, 2000. 19. Suresh, S., Omkar, S.N., Mani, V., and Guru Prakash, T. N., "Lift coefficient prediction at high angle of attack using recurrent neural network," Aerospace Science and Technology, 7:595-602, 2003. 20. Fei, H., Zhu, R., Zhou, Z., and Wang, J., ''Aircraft flight parameter detection based on a neural network using multiple hot-film flow speed sensors'', Smart Materials and Structures, 16:1239–1245, 2007. 21. Peyada, N. K., and Ghosh, A. K., ''Aircraft parameter estimation using a new filtering technique based upon a neural network and Gauss-Newton method,'' The Aeronautical Journal, 113(1142): 243-252, April 2009. 22. Samy, I., Postlethwaite, I., and Gu, D.-W., ''NeuralNetwork-Based Flush Air Data Sensing System Demonstrated on a Mini Air Vehicle,'' Journal of Aircraft, 47(1), January–February 2010. 23. Xuan, C.-Z., Chen, Z., Wu, P., Zhang, Y., and Guo, W., “Stu y of Fuzzy Neural Network on Win Velocity control of Low-Speed Wind Tunnel,” International Conference on Electrical and Control Engineering, 2010. 24. Sivanandam, S. N., Sumathi, S., and Deepa, S. N., "Introduction to Fuzzy Logic using MATLAB', Springer, Berlin Heidelberg, 2007. 25. De Jesus Mota, S., and Botez, R., "New identification method based on neural network for helicopters from flight test data," AIAA Atmospheric Flight Mechanics Conference , Chicago, IL, USA, Aug. 10-13, 2009. 26. o ly, N. , and Botez, R.M., ''New Approach for the Identification and Validation of a Nonlinear F/A-18 Model by Use of Neural Networks'', Neural Networks, IEEE Transactions, 21(11):1759-1765, 2010. 27. o ly, N., Botez, R., and Kouba, G., "Identification of a non-linear F/A-18 model by the use of fuzzy logic and neural network methods,"Proceedings of the Institution of Mechanical Engineers, Part G : Journal of Aerospace Engineering, 225(5):559-574, 2011. 28. Mosbah, A. B., Botez, R., and Dao, T. M., "New methodology for calculating flight parameters with neural network - EGD method," AIAA Modeling and Simulation Technologies (MST) Conference, Boston, MA, USA, Aug 19-22, 2013. 29. Kouba, G., Botez, R., and Bo ly, N., "Identification of F/A18 model from flight tests using the fuzzy logic method,"

47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, USA, Jan. 5-8, 2009. 30. Kouba, G., Botez, R., and Boely, N.,"Fuzzy logic method use in F/A-18 aircraft model identification," Journal of Aircraft, 47(1):10-17, 2010. 31. Grigorie, T. L., Botez, R. , and Popov, A. V.,"Adaptive neuro-fuzzy controllers for an open-loop morphing wing system," Proceedings of the Institution of Mechanical Engineers, Part G : Journal of Aerospace Engineering, 223(7):965-975, 2009. 32. Grigorie, T. L., and Botez, R.," Adaptive neuro-fuzzy inference system-based controllers for smart material actuator modeling ", Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, 223(6):655-668, June 1, 2009 33. Grigorie, T. L., and Botez, R.," Applications of fuzzy logic in the design and control of a morphing wing using smart material actuators," Fuzzy controllers, theory and applications:253-296, USA : INTECH, 2011 34. Grigorie, T. L., Botez, R., Popov, A. V., Mamou, M., and Mébarki, Y.,"Application of fuzzy logic in the design and control of a morphing wing using smart material actuators ». In 58th Aeronautics Conference and AGM (AERO 2011) (Montreal, Canada, Apr. 26-28, 2011) Kanata, Canada, Canadian Aeronautics and Space Institute. 35. Grigorie, T. L., Botez, R., Popov, A. V., Mamou, M., and Mébarki, Y., ''An Intelligent Controller based Fuzzy Logic Techniques for a Morphing Wing Actuation System using Shape Memory Alloy'', 52nd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials ,Conference 19th, 4 - 7 April 2011, Denver, Colorado. 36. Grigorie, T. L., Botez, R., Popov, A. V., Mamou, M., and Mébarki, Y.," A hybrid fuzzy logic proportional integralderivative and conventional on-off controller for morphing wing actuation using shape memory alloy Part 1: Morphing system mechanisms and controller architecture design ," Aeronautical Journal, 116(1179):433-449, 2012. 37. Grigorie, T. L., Botez, R., Popov, A. V., Mamou, M., and Mébarki, Y. ,"A hybrid fuzzy logic proportionalintegralderivative and conventional on-off controller for morphing wing actuation using shape memory alloy Part 2: Controller implementation and validation," Aeronautical Journal, 116(1179):451-465, 2012. 38. Smola, A.J. , and Scholkopf , B., ''A tutorial on support vector regression, Statistic and Computing,'' 14:199-222, 2004. 39. Vapnik, V., ''Tree remarks on the support vector method of function estimation,'' Advances in Kernel methods support vector learning, 25-42, 1999. 40. Keerthi, S., ''Efficient tuning of SVM hyperparameters using radius/margin bound and iterative algorithms,'' IEEE Trans. on Neural Networks, 13(5):1225-1229, 2002. 41. Cherkasky V, and Ma, Y., "Practical selection of SVM parameters and noise estimation for SVM regression," neural networks, 17:113-26, 2004. 42. Ping F.P., and Hong W.C., "Support vector machines with simulated annealing algorithms in electricity load forecasting," Energy Conversion and Management, 46(17):2669-2688, 2005. 43. Burke, E., Bykov, Y., Newell, J., and Petrovic, S., "A timepredefined local search approach to exam timetabling problems," IIE Transactions, 36(6):509–528, 2004. 44. Mosbah, A. B., and Dao, T. M., ''Optimimization of Group Scheduling Problem Using the Hybrid Meta-heuristic Page 8 of 8

Extended Great Deluge (EGD) Approach: A Case Study,'' Journal of Management and Engineering Integration, 4(2):1-13, 2011. 45. Mosbah, A., B., ‘‘Optimisation de l'ordonnancement cellulaire avec métaheuristiques,’’ Master thesis, École de technologie supérieure, Quebec University, 2011. 46. Mosbah, A.B. and Dao, T. M., "Optimization of group

scheduling using simulation with the meta-heuristic Extended Great Deluge (EGD) approach," Industrial Engineering and Engineering Management (IEEM), IEEE International Conference, 275-280, 2010. 47. Mosbah, A. B., and Dao, T. M., ''Optimization of Manufacturing Cell Formation with Extended Great Deluge Meta-heuristic Approach,'' International Journal of Services Operations and Informatics, 7(4):280293,2013. 48. Dueck, G., ''New optimization heuristics. The great deluge algorithm and the record-to-record travel,'' Journal of Computational Physics, 104(1):86–92, 1993. 49. Kecman, V., ''Learning and soft computing-support vector machines, Neural networks and fuzzy logic models,'' The MIT Press, 2001. 50. Mosbah, A. B., Flores Salinas, M., Botez, R., and Dao, T. M., "New Methodology for Wind Tunnel Calibration Using Neural Networks - EGD Approach," SAE International Journal of Aerospace, 6(2):761-766, 2013.

Acknowledgments The first author would like to thank Dr Ruxandra Botez for the opportunity to work on this project, for her constant encouragement and support. Thanks are also due to Dr Thien My Dao for his generosity and for his well needed advices. Many thank you to LARCASE team for their cooperation.