Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms
Abstract
:1. Introduction
1.1. Motivation
1.2. Contribution
1.3. Manuscript Organization
2. Nonlinear 6-DOF Dynamic Model of a Research Aircraft
The Aircraft
3. Trim Condition
4. Linearization
5. Parameters
6. Main Result: Linearization by Means of the Genetic Algorithm (GA)
6.1. Genetic Algorithm (GA)
6.2. Linearization via GA
6.2.1. Transition Matrix
6.2.2. Fitness Function
- Line 1: The name of the function is , and the argument is the vector of the variables to be found by the GA in order to minimize this function.
- Line 2: Global variables , , , and t are used in this function, where is the number of samples considered during the linearization, is a matrix containing the values of the nonlinear states , from the numerical simulation at desired instants, represents the vector of initial conditions, and is a vector with the corresponding time instants at which matrix has been obtained.
- Line 3: Matrix is constructed from elements of vector x from a previous iteration.
- Line 4: Matrix is initialized as a matrix of zeros with rows (number of samples) and four columns (states ).
- Lines 5–7: A loop iterates over the elements of t, using the exponential matrix operation to compute the states of at each one of the time instants included in vector t.
- Lines 8–11: Squared errors for each state, with .
- Line 12: The total error f is computed as the square root of the sum of squared errors.
- Line 13: The function returns f as the result of the fitness evaluation.
6.3. Numerical Simulations
6.3.1. Simulation 1: GA Linearization Using 31 Samples
Algorithm 1 Function |
|
6.3.2. Simulation 2: GA Linearization Using 66 Samples
6.4. Linearization via Least Squares (LS) Method
7. Conclusions
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
Notation and Definitions
A | Matrix for LTI longitudinal dynamics |
Jacobian matrix of nonlinear equations with respect to the state variables | |
B | Input Matrix for LTI system |
Jacobian matrix of nonlinear equations with respect to inputs | |
Cosine function | |
C | Matrix for output definition in state form |
E | Jacobian matrix of nonlinear equations with respect to the first-order state variables |
Force vector in body frame | |
Jacobian matrix of nonlinear equations with respect to the first-order state variables | |
h | Step size for numerical derivative, |
Aircraft flying height | |
H | Weight diagonal matrix for optimization criteria |
Transformation matrix of angular rates from body axes to the inertial system | |
Inertial tensor of body center of gravity | |
m | Aircraft mass |
Moment vector in body axes | |
Moment due to pitch rate | |
Moment due to horizontal velocity in body axes | |
Moment due to vertical speed in body axis | |
p | Roll rate in x body axis |
Angular acceleration in the x-axis in body frame | |
Aircraft displacement to the east direction (navigation frame) | |
Aircraft displacement to north direction (navigation frame) | |
q | Pitch rate in y body axis |
Angular acceleration in the y-axis in body frame | |
r | Yaw rate in body axes |
Angular acceleration on the z-axis in body frame | |
Sine function | |
T | Step size for genetic algorithm |
u | Flight velocity in x-axes in body frame from the equilibrium state (cruise flight) |
U | Input vector control |
Flight acceleration in x-axes in the body frame and from the equilibrium state (cruise flight) | |
Input vector control for equilibrium condition | |
Input for equilibrium condition in the linear system | |
v | Lateral component for flight speed from equilibrium condition, in body frame |
Aircraft velocity from the center of mass to Earth axis, super-index b indicates body frame | |
Flight velocity in trim condition | |
Lateral component for acceleration from equilibrium point, in body frame |
w | Vertical component for flight velocity from equilibrium condition, in body frame |
Vertical component for acceleration from equilibrium condition, in body frame | |
X | Twelve dimensional state vector |
x | x displacement on body frame |
Initial condition for non linear model | |
Vector state in trim condition | |
Forces in x body axis due to pitch angle | |
State values in equilibrium point for the linear system definition | |
First-order state values in equilibrium point for the linear system definition | |
First-order vector state | |
Matrix containing the values of the nonlinear states | |
Forces in x body axis due to pitch rate | |
Forces in the x-body axis due to velocity in the x-axis of the body frame | |
Forces in the x-body axis due to velocity in the z-axis of the body frame | |
y | Lateral displacement in body axis |
z | Vertical displacement in body axis |
Forces in z-body axis due to pitch angle | |
Forces in z-body axis due to pitch rate | |
Forces in the z-body axis due to velocity in the x-axis of the body frame | |
Forces in the z-body axis due to velocity in the z-axis of the body frame | |
Longitudinal flight path angle | |
Pitch angle | |
Pitch angle for trim condition | |
Pitch rate | |
Vector of Euler angle rates, angular rates of aircraft respect to the inertial system | |
Roll angle | |
Roll rate | |
Transition matrix | |
Transition matrix with initial condition | |
Yaw angle | |
Yaw rate | |
Angular rate vector in body axis | |
Angular acceleration vector in body axis |
References
- Hanks, H.T.; Tran, H.T. Mathematical and Experimental Modeling of Physical and Biological Processes; Chapman and Hall/CRC: Boca Raton, FL, USA, 2009. [Google Scholar]
- Nelles, O. Nonlinear System Identification: From Classical Approaches to Neural Networks and Fuzzy Models; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2001. [Google Scholar]
- Martinez, O.A.; Cardona, M. State of the Art and Future Trends on Unmanned Aerial Vehicle. In Proceedings of the 2018 International Conference on Research in Intelligent and Computing in Engineering (RICE), San Salvador, El Salvador, 22–24 August 2018; pp. 1–6. [Google Scholar] [CrossRef]
- Brunton, S.L.; Proctor, J.L.; Kutz, N.J. Discovering Governing Equations From Data by Sparse Identification of Nonlinear Dynamical Systems. Proc. Natl. Acad. Sci. USA 2016, 113, 3932–3937. [Google Scholar] [CrossRef] [PubMed]
- Grauer, J.A.; Morelli, E.A. Introduction to the Advances in Aircraft System Identification from Flight Test Data Virtual Collections. J. Aircr. 2023, 60, 1329. [Google Scholar] [CrossRef]
- Hosseini, B.; Steinert, A.; Hofmann, R.; Fang, X.; Steffensen, R.; Holzapfe, F. Advancements in the Theory and Practice of Flight Vehicle System Identification. J. Aircr. 2023, 6, 1419–1436. [Google Scholar] [CrossRef]
- Jategaonkar, R.V. Flight Vehicle System Identification: A Time-Domain Methodology; American Institute of Aeronautics and Astronautics, Inc.: Reston, VA, USA, 2015. [Google Scholar]
- Deiler, C.; Mönnich, W.; Seher-Wei, S.; Wartmann, J. Retrospective and Recent Examples of Aircraft and Rotorcraft System Identification at DLR. J. Aircr. 1996, 33, 1–26. [Google Scholar] [CrossRef]
- Tischler, M.B.; Remple, R. System Identification. Engineering Methods with Flight Test Examples; AIIA Education: Richmond, VA, USA, 2006. [Google Scholar]
- Cedillo, S.G.T.; Bonello, P. Empirical identification of the inverse model of a squeeze-film damper bearing using neural networks and its application to a nonlinear inverse problem. J. Vib. Control 2018, 24, 357–378. [Google Scholar] [CrossRef]
- Kirkpatrick, K.; May, J.; Valasekz, J. Aircraft System Identification Using Artficial Neural Networks. In Proceedings of the 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Grapevine, TX, USA, 7–10 January 2013; pp. 1–12. [Google Scholar]
- Ulinowicz, M.; Narkiewicz, J. Aircraft parameter Identification Using Genetic Algorithm. In Proceedings of the 29th Congress of the International Council of the Aeronautical Sciences, St. Petersburg, Russia, 7–12 September 2014; pp. 317–334. [Google Scholar]
- Bruce, P.; Kellett, M. Maximum likelihood identification of linear aircraft dynamics using a hybrid genetic algorithm. In Proceedings of the 36th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 12–15 January 1998; pp. 357–378. [Google Scholar]
- Roy, A.G.; Peyada, N.K. Longitudinal Aircraft Parameter Estimation Using Neuro-Fuzzy and Genetic Algorithm Based Method. In Proceedings of the AIAA Atmospheric Flight Mechanics Conference, Grapevine, TX, USA, 9–13 January 2017; pp. 1–12. [Google Scholar]
- Anderson, M.B. Genetic Algorithms in Aerospace Design: Substantial Progress, Tremendous Potential; Technical Report; Sverdrup Technology Inc.: Allen Park, MI, USA; TEAS Group: Manchester, UK, 2002. [Google Scholar]
- Brunton, S.L.; Kutz, J.N. Data Driven Science and Engineering Machine Learning, Dynamical Systems, and Control; Department of Applied Mathematics, University of Washington: Philadelphia, PA, USA, 2017. [Google Scholar]
- Miri, N.; Mohammadzaheri, M.; Chen, L. An enhanced physics-based model to estimate the displacement of piezoelectric actuators. J. Intell. Mater. Syst. Struct. 2015, 26, 1442–1451. [Google Scholar] [CrossRef]
- Kailath, T.; Sayed, A.H.; Hassibi, B. Linear Estimation; Prentice Hall: Englewood Cliffs, NJ, USA, 2000. [Google Scholar]
- Ogata, K. Modern Control Engineering; Pearson Education: Saddle River, NJ, USA, 2010. [Google Scholar]
- Lambrechts, P.; Bennani, S.; Looye, G.; Moormann, D. The RCAM design challenge problem description. In Proceedings of the Robust Flight Control, Toulouse, France, 7 March 1997; Magni, J.F., Bennani, S., Terlouw, J., Eds.; Springer: Berlin/Heidelberg, Germany, 1997; pp. 149–179. [Google Scholar]
- Lambrechts, P.F. Technical Publication TP-088-3, GARTEUR-FM(AG08); Technical Report; Group for Aeronautical Research and Technology in Europe (GARTEUR): Madrid, Spain, 1997. [Google Scholar]
- Stevens, B.L.; Lewis, F.L.; Johnson, E.N. Aircraft Control and Simulation: Dynamics, Controls Design, and Autonomous Systems, 3rd ed.; Wiley: Hoboken, NJ, USA, 2016. [Google Scholar]
- Hong-Gi, L. Linearization of Nonlinear Control Systems; Springer: Singapore, 2022. [Google Scholar]
- Kiusalaas, J. Numerical Methods in Engineering with MATLAB®, 3rd ed.; Cambridge University Press: Cambridge, UK, 2015. [Google Scholar] [CrossRef]
- Miller, G. Numerical Analysis for Engineers and Scientists; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar] [CrossRef]
- Roskam, J. Airplane Flight Dynamics and Automatic Flight Controls; DAR Corporation: Lawrence, KS, USA, 2003. [Google Scholar]
- Blakelock, J.H. Automatic Control of Aircraft and Missiles; Wiley: New York, NY, USA, 1991. [Google Scholar]
- McLean, D. Automatic Flight Control Systems; Prentice Hall: New York, NY, USA, 1990. [Google Scholar]
- Holland, J.H. Adaptation in Natural and Artificial Systems; University of Michigan Press: Cambridge, MA, USA, 1975. [Google Scholar]
- Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning; Addison-Wesley: Boston, MA, USA, 1989. [Google Scholar]
- Mitchell, M. An Introduction to Genetic Algorithms; MIT Press: Cambridge, MA, USA, 1998. [Google Scholar]
- Haupt, R.L.; Haupt, S.E. Practical Genetic Algorithms; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
- López-González, A.; Meda-Campaña, J.A.; Hernández-Martínez, E.G.; Paniagua-Contro, P. Multi robot distance based formation using Parallel Genetic Algorithm. Appl. Soft Comput. 2020, 86, 105929. [Google Scholar] [CrossRef]
- Liu, S.; Acar, P. Parameter Space Exploration of Cellular Mechanical Metamaterials Using Genetic Algorithms. AIAA J. 2023, 61, 3633–3643. [Google Scholar] [CrossRef]
- Isakhani, H.; Yue, S.; Xiong, C.; Chen, W. Aerodynamic Analysis and Optimization of Gliding Locust Wing Using Nash Genetic Algorithm. AIAA J. 2021, 59, 4002–4013. [Google Scholar] [CrossRef]
- Shu, J.I.; Wang, Y.; Brown, A.; Kaminsky, A. Genetic-Algorithm-Guided Development of Parametric Aeroelastic Reduced-Order Models with State-Consistence Enforcement. AIAA J. 2023, 61, 3976–3994. [Google Scholar] [CrossRef]
- MathWorks. Genetic Algorithm Options (Global Optimization Toolbox); Mathworks: Natick, MA, USA, 2023. [Google Scholar]
Linearization via | MSE | Execution Time |
---|---|---|
Section 4 | —– | |
GA (31 samples) | ||
GA (66 samples) | ||
LS (31 samples) | —– | |
LS (full dataset) | —– |
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. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Peña-García, R.; Velázquez-Sánchez, R.D.; Gómez-Daza-Argumedo, C.; Escobedo-Alva, J.O.; Tapia-Herrera, R.; Meda-Campaña, J.A. Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms. Aerospace 2024, 11, 142. https://doi.org/10.3390/aerospace11020142
Peña-García R, Velázquez-Sánchez RD, Gómez-Daza-Argumedo C, Escobedo-Alva JO, Tapia-Herrera R, Meda-Campaña JA. Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms. Aerospace. 2024; 11(2):142. https://doi.org/10.3390/aerospace11020142
Chicago/Turabian StylePeña-García, Raymundo, Rodolfo Daniel Velázquez-Sánchez, Cristian Gómez-Daza-Argumedo, Jonathan Omega Escobedo-Alva, Ricardo Tapia-Herrera, and Jesús Alberto Meda-Campaña. 2024. "Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms" Aerospace 11, no. 2: 142. https://doi.org/10.3390/aerospace11020142
APA StylePeña-García, R., Velázquez-Sánchez, R. D., Gómez-Daza-Argumedo, C., Escobedo-Alva, J. O., Tapia-Herrera, R., & Meda-Campaña, J. A. (2024). Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms. Aerospace, 11(2), 142. https://doi.org/10.3390/aerospace11020142