Piloted and batch simulations of the aeroservoelastic response are essential tools in the development of advanced ight control systems. In these simulations the number of differential equations must be suf ciently large to yield the required accuracy, yet small enough to enable real-time evaluationsof the aircraft ying qualities. The challenge of these con icting demands is reinforced by nonlinearities in the quasi-steady equations of motion and by the complex characteristics of the oscillatory forces. Our solution to the problem is based on a unique formulation that eliminates the need for auxiliary state variables to represent the unsteady aerodynamics. We also address transformations from the mean ight path axes to a body axes coordinate system and describe how the structural dynamic equations of motion are integrated with the quasi-steady, nonlinear, six-degree-of-freedom plant model. The uni ed model, which accurately preserves the roots of the dynamic aeroelastic system, includes provisions for control surface inputs and atmospheric turbulence.
Nomenclature
A, B= matrices in state-space formulation e = vector of generalized body axes coordinates F m = aerodynamic generalized force f, g, h = aircraft displacements in mean ight path system G NL = solution to nonlinear six-degree-of-freedom equations M m = inertial generalized force P = vector derived in p transform p, q, r = aircraft angular rates about body axes q = dynamic pressure S m = structural generalized force T = transformation matrix U = mean air ow velocity u, v, w = aircraft linear rates along body axes X, Y, Z = aircraft mean ight path axes system X ,Ȳ ,Z = aircraft body axes system z = state variable vector in body axes system a m,n = real part of aerodynamic approximation b m ,n = imaginary part of aerodynamic approximation c m, n = element of structural damping matrix D G = linear equations added to G NL d = vector of inputs " = vector of generalized mean ight path coordinateś = state variable vector in mean ight path system j m,n = element of structural stiffness matrix k = complex eigenvalue l m, n = element of mass matrix r = real part of eigenvalue U I = imaginary part of eigenvector U R = real part of eigenvector x = circular frequency