Acknowledgment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 256 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 256

Abstract The present paper is devoted to our recent results on an aircraft wing model in a subsonic airflow. The model is governed by a system of two coupled integro-differential equations and a twoparameter family of boundary conditions modeling the action of the self-straining actuators. The system of equations of motion is equivalent to a single operator evolution-convolution equation in the energy space. The Laplace transform of the solution involves the so-called generalized resolvent operator, which is a finite meromorphic operator-valued function of the spectral parameter defined on the complex plane with a branchcut along the negative real semiaxis. The main findings in this study are (a) asymptotic distribution of the poles (called aeroelastic modes) of the generalized resolvent and (b) the Riesz basis property of the mode shapes.