Hostname: page-component-586b7cd67f-2plfb Total loading time: 0 Render date: 2024-11-25T06:19:46.517Z Has data issue: false hasContentIssue false

An efficient reduced-order framework for active/passive hybrid flutter suppression

Published online by Cambridge University Press:  02 May 2022

D. F. Li*
Affiliation:
College of Water Resources and Architectural Engineering, Cold and Arid Regions Water Engineering Safety Research Center, Northwest A & F University, Yangling, Shaanxi, PR China
Z. Z. Wang
Affiliation:
College of Water Resources and Architectural Engineering, Cold and Arid Regions Water Engineering Safety Research Center, Northwest A & F University, Yangling, Shaanxi, PR China
A. Da Ronch
Affiliation:
Engineering and the Environment, University of Southampton, Southampton, UK
G. Chen
Affiliation:
State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace, Xi’an Jiaotong University, Xi’an, Shaanxi, PR China
*
*Corresponding author email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Flutter suppression is an important measure to improve fatigue life and enhance the performance of aircraft in modern aircraft design. In order to design more effective controllers for flutter suppression with high efficiency, an efficient reduced-order framework for active/passive hybrid flutter suppression is proposed. The traditional CFD-based ROMs have been successfully applied to active flutter suppression with high accuracy and efficiency. But, when a structure modification is made such as in aeroelastic tailoring and aeroelastic structural optimisation, the structural model should be updated, and the expensive, time-consuming CFD-based ROMs have to be reconstructed; such a process is impractical for passive flutter suppression. To overcome the realistic challenge, an efficient reduced-order framework for active/passive hybrid flutter suppression is proposed by extending an efficient aeroelastic CFD-based POD/ROM which we have developed. The proposed framework is demonstrated and evaluated using an improved AGARD 445.6 wing model. The results show that the proposed framework can accurately predict the aeroelastic response for active/passive hybrid flutter suppression with high efficiency. It provides a powerful tool for active/passive hybrid flutter suppression, and therefore, is ideally suited to design more effective controllers, and may have the potential to reduce the overall cost of aircraft design.

Type
Research Article
Copyright
© The Author(s), 2022. Published by Cambridge University Press on behalf of Royal Aeronautical Society

Nomenclature

F

nonlinear numerical flux

w

vector of conservative flow variables

u

structural displacement vector

u c

control surface deflection

A

diagonal matrix of cell volumes

H

gradient of the flux function with respect to the vector of fluid variables

E

gradient of the cell volumes with respect to the generalised coordinates

G

gradient of the flux function with respect to the generalised coordinates

C

gradient of the flux function with respect to the velocities

L

gradient of the flux function with respect to control surface deflection

Q

semi-positive definite weighting symmetric matrix

R

positive definite weighting symmetric matrix

K

stiffness matrix

M

mass matrix

$ \boldsymbol{\Omega} $

subspaces in the proper orthogonal decomposition

$ \textbf{x}^{n\times 1} $

state vector of the linearised time domain equations

$ \textbf{X} $

snapshots matrix

$\boldsymbol{\Psi}_{r}$

r-order proper orthogonal decomposition basis

$\xi_i $

the contribution of the i-th snapshot to the original system

$\boldsymbol{\Theta} $

physical displacement

${\lambda ^i}$

i-th eigenvalue

$\boldsymbol{\phi}^{i}$

i-th eigenvector

E

Young’s modulus of the wing structure

$ \rho $

density of the wing structure

Ma

Mach number

AOA

angle of attack

$ V_{\infty} $

freestream speed

1.0 Introduction

Flutter is a self-excited oscillation phenomenon resulting from the dynamic instability of a coupling between unsteady aerodynamic loads, inertial and elastic forces of structure. When flutter occurs, the vibration amplitudes of the wing become dynamically unstable, significantly affecting the fatigue life of aircraft, and may even lead to an unexpected catastrophe. Flutter suppression is one of important objective of many investigations in the aeroelastic field, not only to ensure the flight stability and safety of aircraft and avoid catastrophic failures, but also to enhance aircraft performance and expand the flutter boundary. Different strategies have been proposed for flutter suppression, including passive control and active control [Reference Kuzmina, Amiryants, Schweiger, Cooper, Amprikidis and Sensberg1Reference Liu3]. Passive control is the earliest technology that has been widely applied in flutter suppression. It takes a mass balance and change in structural stiffness as main measures for flutter suppression. However, passive control may cause weight increase, which has adverse effects on the flight performance. Active control is a relatively new technology when employed for flutter suppression. Such technologies change the aerodynamic load through the deflection of the control surface to suppress flutter. The benefits of active control usually have a higher control authority and versatility as compared to passive control, while the benefits of the latter are that passive control is inherently stable and much less difficult and expensive to implement.

In order to design more effective controllers for flutter suppression, an active/passive hybrid control method could be beneficial [Reference Reich4, Reference Gregory, Mart and Crawley5]. Leão et al. [Reference Le, de Lima, Donadon and Cunha-Filho6] pointed out that a better understanding and improvements on the aeroelastic behaviour of wing structures using active/passive hybrid flutter suppression strategies are nowadays key issues in designing advanced aerospace structures. To design a more effective controller for flutter suppression, some researchers have proposed and developed several active/passive hybrid methods. Barzegari et al. [Reference Barzegari, Dardel and Fathi7] demonstrated an active/passive hybrid method for aeroelastic control, shape memory alloy wires are employed for passive control of the wings by increasing the wing’s stiffness, and the active control system design is based on the state-dependent Riccati equation method. Roeland et al. [Reference de Breuker, Binder and Wildschek8] developed a dedicated aeroservoelastic model that can quickly assess gust loads in the presence of flaps and bending-twisting coupling, and demonstrated that a combination of combined active and passive control method outperforms separate flutter control and suppression methods. Silva et al. [Reference Silva, Silvestre, Donadon, Santos, Guimar, Neto, da Silva, Versiani, Gonzalez and Bertolin2] described a flutter suppression testbed based on active/passive hybrid control. The active control is used to flutter suppression by changing aerodynamic load through deflection of control surface. The passive control employs a shape-memory alloy to provide additional torsional stiffness no full stop-torsional stiffness and the results show that an active/passive hybrid strategy can significantly suppress flutter at a relatively low cost. Schweiger et al. [Reference Schweiger, Krammer and Coetzee9] presented that the combination of passive control with active control for flutter suppression, the aircraft performance can be improved at considerable weight savings. Weisshaar et al. [Reference Weisshaar and Duke10] combined wing aeroelastic stiffness tailoring with active control surface design to create a “control-friendly structure” that reduces the flutter excitation. It uses a full-span deforming control surfaces to provide effective subsonic induced drag control with precise control of the spanwise lift distribution to achieve the purpose of flutter suppression. Handojo et al. [Reference Handojo, Lancelot and de Breuker11] presented an active/passive hybrid loads alleviation method on a generic mid-range aircraft configuration. The passive control is carried out by changing the wing box stiffness properties using composite material and variable lamination parameters. For the active control, the manoeuver loads alleviation is implemented by defining a symmetric aileron deflection during pull-up and push-down manoeuvers. All the above research calls for the development of an efficient framework for active/passive hybrid flutter suppression with high efficiency, which has the ability to design more effective controllers.

High-performance modern aircraft often cruise in transonic regime to achieve high performance. Therefore, flutter suppression is a particularly important aspect for modern aircraft design. For flutter suppression, a high-fidelity aeroelastic analysis model is highly needed. Linear potential flow theory provides a low-cost and high efficiency way for aeroelastic analysis. Unfortunately, linear aerodynamic theories usually fail due to the inherently nonlinear characteristics of the transonic flow. Computational fluid dynamics (CFD) is a feasible alternative method to model these flow nonlinearities. However, these full order models (FOMs) including finite element analysis (FEA) and CFD require large computer memory and computing time, which is impractical to be routinely used in flutter suppression.

To overcome the computational cost limitation of CFD-based aeroelastic analysis, some researchers have proposed and developed several CFD-based reduced order models (ROMs). These models extract key data of the fluid system to generate a low-dimensional system, it may significantly reduce the computational cost while maintain similar accuracy to the FOMs [Reference Chen, Zhao and Huang12, Reference Kou and Zhang13]. System identification [Reference Zhang, Lv, Diwu and Zhong14] and Proper Orthogonal Decomposition (POD) [Reference Lu, Jin, Chen, Yang, Hou, Zhang, Li and Fu15] are the two most commonly used ROMs in CFD-based aeroelastic analysis. The CFD-based ROMs have been also exercised for transonic [Reference Zhou, Chen and Li16] and hypersonic [Reference Cao, Nie, Pan, Cai and Qu17] aeroelastic analysis, limit-cycle oscillation (LCO) [Reference Dowell, Thomas and Hall18], gust response analysis [Reference Zhou, Chen, da Ronch and Li19], aerodynamic shape optimisation [Reference Wu, Zhang, Peng and Wang20, Reference Li, Li, Zhang, Lu and Yuan21], and blade aerodynamic force [Reference Zhou Li, Yang, Luo, Zhang and Ni Yuan22, Reference Sajadmanesh, Mojaddam, Mohseni and Nikparto23]. In particular, the last decade has witnessed various studies on the developments of the CFD-based ROMs for active flutter suppression. Chen et al. [Reference Gang, Jian and Yueming24] proposed an active flutter control law design method based on an improved POD-BT/ROM and demonstrated by the Goland wing+/store system, and balanced truncation method was applied to further reduce the dimension of the time-domain POD/ROM. Song et al. [Reference Song, Qian, Wang, Pant, Chin and Brenner25] presented a holistic methodology and framework that uniquely integrates the aerodynamics ROM, structural dynamics model, and active structural control into a single platform to enable integrated, ultra-fast ROMs for efficient aeroelastic and aeroservoelastic analysis, and controller design. Nie et al. [Reference Nie, Yang and Zhang26] used a CFD-based ROM coupled with computational structural dynamics (CSD) to predict the flutter. The closed-loop control systems designed by the sliding mode control (SMC) and linear quadratic Gaussian (LQG) control were constructed with ROM/CSD to suppress the AGARD 445.6 wing flutter. The results indicate that the CFD-based ROM can replace the CFD calculation accurately and the designed control laws for the coupled ROM/CSD system can effectively suppress flutter. Liu et al. [Reference Liu, Huang, Zhao and Hu27] designed a ROM for an elastic wing with multiple control surfaces. The CFD-based ROM enables one to predict the unsteady aerodynamic loads induced by the deflections of the wing and leading-edge and trailing-edge control surfaces with high accuracy and efficiency. Hoseini et al. [Reference Hoseini and Hodges28] using POD/ROM to design a flutter suppression system for a joined 3D/1D finite element model of damaged High Altitude-Long Endurance (HALE) aircraft wings. Patterson et al. [Reference Patterson and Friedmann29] describe a surrogate-based ROM for simulating the unsteady aerodynamic effects of pulsed-jet flow control on a two-dimensional aerofoil. The ROM is constructed by extracting time-domain sample data for the unsteady lift, moment and drag from CFD simulations of the actuation. All of above-mentioned studies show that CFD-based ROMs have shown acceptable accuracy and efficiency for active flutter suppression.

Unfortunately, traditional CFD-based ROMs are all formulated for a fixed flight condition for a frozen model configuration, such as structural parameters (e.g. mass, stiffness). Even worse is that the influence of structural mode shapes change will also affect the solution of fluid model. In order to ensure accuracy, both the structural model and CFD-based ROM have to be reconstructed as the exact model, which destroys the computational efficiency of the CFD-based ROM. Therefore, traditional CFD-based ROMs cannot be directly applied for passive flutter suppression. In recent years, researchers have been focusing and trying to make up the shortcoming of traditional CFD-based ROMs. Fenwick et al. [Reference Fenwick, Jones and Gaitonde30] used a linear interpolation on a set of available ROMs to obtain a new ROM without the need for a rebuild, for mass redistribution (e.g. fuel load redistribution) of a wing. Zhang et al. [Reference Zhang, Chen and Ye31] demonstrated a method that can replace the CFD solver used in the process of existing unsteady CFD-based auto regressive with eXogenous input (ARX) ROM based on radial basis function (RBF) interpolation functions. It can quickly produce the aerodynamic responses corresponding to the mode excitations for arbitrary mode shapes associated with changes of the root boundary condition. Winter et al. [Reference Winter, Heckmeier and Breitsamter32] presented two novel CFD-based ROMs methodologies robust to variations in the structural mode shapes due to the additional lumped mass. Li et al. [Reference Li, Kou and Zhang33] presented an unsteady aerodynamic model based on LSTM (long short-term memory) network from deep learning theory, and the model can accurately capture the dynamic characteristics of aerodynamic and aeroelastic systems for varying flow and structural parameters. We have been developed and implemented an efficient aeroelastic CFD-based POD/ROM [Reference Chen, Li, Zhou, da Ronch and Li34], which avoids the burden of reconstructing CFD-based POD/ROMs associated the exact model. The efficient aeroelastic CFD-based POD/ROM has been successfully applied to passive flutter suppression with high accuracy and efficiency, such as aeroelastic tailoring [Reference Li, Gong, da Ronch, Chen and Li35] and aeroelastic structural optimisation [Reference Li, da Ronch, Chen and Li36]. However, most of above-mentioned studies only considered the CFD-based ROMs for active or passive flutter suppression alone.

With the previous paragraphs as background, in order to design more effective controllers for flutter suppression with high efficiency, this paper proposes an efficient reduced-order framework for active/passive hybrid flutter suppression by extending the efficient aeroelastic CFD-based POD/ROM. The proposed efficient reduced-order framework would be convenient to design more effective controllers with high efficiency, and provide a potential efficient tool for flutter suppression, which will allow the reduction of overall cost of the aircraft design. This paper is constructed as follows: Section 2 gives the brief introduction to the standard CFD-based POD/ROM construction procedure for active flutter suppression; section 3 the proposed efficient reduced-order framework for active/passive hybrid flutter suppression is described in detail; In section 4, the proposed framework was demonstrated and evaluated by an improved AGARD 445.6 aeroelastic wing model. Section 5 is the conclusion and discussion of the proposed reduced-order framework.

2.0 Overview of active flutter suppression

To understand the fundamental principles underlying of the proposed efficient reduced-order framework for active/passive hybrid flutter suppression, it is essential to first understand the standard CFD-based POD/ROM construction procedure for active flutter suppression.

2.1. Model reduction for aeroelastic problem

POD is a method to rebuild the behaviour of the overall system with a small number of degrees of freedom as in the case of unsteady flows for aeroelastic analysis [Reference Hall, Thomas and Dowell37, Reference Lieu, Farhat and Lesoinne38]. Here, POD is employed to construct a CFD-based ROM for Euler equations in active flutter suppression. Compared with the traditional aeroelastic linearised flow equations, the aerodynamic force is jointly determined by the structural deformation and the deflection of control surface. First, the unsteady Euler equations are linearised around a mean flow solution (equilibrium), with the flow solution satisfying the relation $\textbf{F}({\textbf{w}_0},{\rm{ }}{\textbf{u}_0},{\rm{ }}{\textbf{u}_{c0}})$ , where ${\textbf{w}_0}\,=\,{\boldsymbol{0}}$ , $\dot{\textbf{w}}_0\,=\,{\boldsymbol{0}}$ , $\textbf{u}_0\,=\,{\boldsymbol{0}}$ , $\dot{\textbf{u}}_0\,=\,{\boldsymbol{0}}$ , $\textbf{u}_{c0}\,=\,{\boldsymbol{0}}$ , $\dot{\textbf{u}}_{c0}\,=\,{\boldsymbol{0}}$ . Here, F indicates the nonlinear numerical flux resulting from the spatial discretisation, w is the vector of conservative flow variables, u is the structural displacement vector, $\textbf{u}_c$ denotes the control surface deflection. Supposing that $\partial \textbf{w},{\rm{ }}\partial \textbf{u},{\rm{ }}\partial \textbf{u}_c$ are small perturbations around the steady state variables. Then, one obtains the linearised flow equations about the steady-state solution $(\textbf{w}_0,{\rm{ }}\textbf{u}_0,{\rm{ }}\textbf{u}_{c0})$ through a first order Taylor’s expansion:

(1) \begin{align}\textbf{A}_{\boldsymbol{0}}\dot {\textbf w} + \textbf{Hw} + (\textbf{E} + \textbf{C})\dot {\textbf u} + \textbf{Gu} + \textbf{Lu}_c = \boldsymbol{0}\end{align}

where

\begin{align*}\textbf H = \frac{{\partial \textbf F}}{{\partial \textbf w}}({\textbf w_0},{\textbf u_0},{\textbf u_{c0}})\quad \textbf G = \frac{{\partial \textbf F}}{{\partial \textbf u}}({\textbf w_0},{\textbf u_0},{\textbf u_{c0}})\quad \textbf E = \frac{{\partial \textbf A}}{{\partial \textbf u}}{\textbf w_0}\quad \textbf C = \frac{{\partial \textbf F}}{{\partial \dot {\textbf u}}}({\textbf w_0},{\textbf u_0},{\textbf u_{c0}})\quad \textbf L = \frac{{\partial \textbf F}}{{\partial {\textbf u_c}}}({\textbf w_0},{\textbf u_0},{\textbf u_{c0}})\end{align*}

Here, A denotes the diagonal matrix of cell volumes. The matrices H, E, G, C and L are the first-order terms of a Taylor expansion of the non-dimensional numerical flux function around the equilibrium ( ${\textbf w_0},{\textbf u_0},{\textbf u_{c0}}$ ). The matrix H is the gradient of the numerical flux function with respect to the vector of fluid variables. The matrix E indicates the gradient of the cell volumes with respect to the generalised coordinates. Matrices G and C are the gradients of the flux function with respect to the generalised coordinates and their velocities, respectively. Finally, the matrix L indicates the gradient of the numerical flux function with respect to control surface deflection. Note that the matrices G, E and C need to be recomputed for modified structure.

In active flutter suppression system, the elastic influence of the control surface can be ignored, so the linearisation of the structural dynamic equation around an equilibrium state can be written as:

(2) \begin{align}\bar {\textbf M}\ddot {\textbf u} + {\bar {\textbf D}_0}\dot {\textbf u} + {\bar {\textbf K}_s}\textbf u = {\textbf P_0} \textbf w\end{align}

where

\begin{align*}{\bar {\textbf K}_0} = \frac{{\partial {\textbf f^{int}}}}{{\partial \textbf u}}({\textbf u_0},{\dot {\textbf u}_0})\quad {\bar {\textbf K}_s} = {\bar {\textbf K}_0} - \frac{{\partial {\textbf f^{ext}}}}{{\partial \textbf u}}({\textbf w_0},{\textbf u_0})\quad {\bar {\textbf D}_0} = \frac{{\partial {\textbf f^{int}}}}{{\partial \dot {\textbf u}}}({\textbf u_0},{\dot {\textbf u}_0})\quad {\textbf P_0} = \frac{{\partial {\textbf f^{ext}}}}{{\partial \textbf w}}({\textbf w_0},{\textbf u_0})\end{align*}

For the analysis of the system stability, the terms $\dfrac{{\partial {\textbf f^{ext}}}}{{\partial \textbf u}}$ and ${\bar {\textbf D}_0}$ can be neglected, which leads to the structural dynamic equation:

(3) \begin{align}\bar {\textbf M}\ddot {\textbf u} + {\bar {\textbf K}_0} \textbf u = {\textbf P_0} \textbf w\end{align}

The linearised CFD-based aeroelastic equations of Eq. (1) and Eq. (3) are too large for aeroelastic system. Therefore, the POD is used to reduce the full order models (FOMs). Denote $\{ {\textbf x^k}\} $ , $k = 1,2,3 \ldots m$ , a set of data, with ${\textbf x^k}$ is the n-dimensional space, and m is the number of snapshots. The POD method searches an m-dimensional proper orthogonal subspace, ${\boldsymbol \Psi} \in {\textbf R^{n \times m}}$ , to minimise the mapping errors from ${\boldsymbol \Psi}$

(4) \begin{align} \textbf G = \mathop {\min }\limits_{\boldsymbol \Phi} \sum\limits_{k = 1}^m {\left\| {{\textbf x^k} - {\boldsymbol \Omega} {\boldsymbol \Omega} {^T}{\textbf x^k}} \right\|} = \sum\limits_{k = 1}^m {\left\| {{\textbf x^k} - {\boldsymbol \Psi} {{\boldsymbol \Psi} ^T}{\textbf x^k}} \right\|} ,\quad {{\boldsymbol \Omega} ^H}{\boldsymbol \Omega} = {\textbf I} \end{align}

The minimisation problem is equivalent to [Reference Lieu39]:

(5) \begin{align} \textbf H = \mathop {\max }\limits_{\boldsymbol \Phi} \sum\limits_{k = 1}^m {\frac{{\left\langle {{{\left( {{\textbf x^k},{\boldsymbol \Omega} } \right)}^2}} \right\rangle }}{{{{\left\| {\boldsymbol \Omega} \right\|}^2}}}} = \sum\limits_{k = 1}^m {\frac{{\left\langle {{{\left( {{\textbf x^k},{\boldsymbol \Psi} } \right)}^2}} \right\rangle }}{{{{\left\| {\boldsymbol \Psi} \right\|}^2}}}} ,\quad {{\boldsymbol \Omega} ^H}{\boldsymbol \Omega} = \textbf I\end{align}

The constraint optimisation problem in Eq. (5) is transformed into the following Lagrange equation

(6) \begin{align} J\left( {\boldsymbol \Omega} \right) = \sum\limits_{k = 1}^m {{{\left( {{\textbf x^k},{\boldsymbol \Omega} } \right)}^2}} - \lambda \left( {{{\left\| {\boldsymbol \Omega} \right\|}^2} - 1} \right)\end{align}

Solving the partial derivative of the objective function $J(\boldsymbol \Omega )$ with respect to $\boldsymbol \Omega $ gives

(7) \begin{align} \frac{d}{{d{\boldsymbol \Omega} }}J\left( {\boldsymbol\Omega} \right) = 2{\textbf X}{\textbf X^H}{\boldsymbol \Omega} - 2\lambda {\boldsymbol \Omega} \end{align}

where $\textbf X = \{ {\textbf x_1},{\textbf x_2}, \cdots {\textbf x_m}\} \in {\textbf R^{n \times m}}$ is a matrix containing m snapshots as columns. By setting Eq. (7) to zero; thus, the following equation is obtained:

(8) \begin{align} \left( {\textbf X{\textbf X^T} - \xi \textbf I} \right){\boldsymbol \Psi} = 0\end{align}

The problem is transformed into solving the eigenvalue problem of the POD kernel, $\textbf X{\textbf X^H}$ . The eigenvalue problem has very large size, as $\textbf X{\textbf X^H} \in {\textbf R^{n \times n}}$ . Because $\textbf X{\textbf X^H}$ and $\textbf X{\textbf X^H}$ have the same eigenvalues, so we can obtain ${\boldsymbol \Psi} $ as:

(9) \begin{align} \left\{ \begin{array}{l}\textbf X{\textbf X^T}\textbf V = \textbf V{\boldsymbol\Lambda} \\[5pt] {\boldsymbol \Psi} = {\textbf {XV}}{\boldsymbol \Lambda ^{ - 1/2}}\end{array} \right.\end{align}

where ${\boldsymbol \Psi} = [{{\boldsymbol \psi} _1},{{\boldsymbol \psi} _2}, \ldots ,{{\boldsymbol \psi} _m}]$ , ${\boldsymbol \Lambda} = diag({\xi _1},{\xi _2} \cdots {\xi _m})$ , ${\xi _1} \ge {\xi _2} \ge ,\ldots, \ge {\xi _m}$ . The value of ${\xi _i}$ represents the contribution of the i-th snapshot to the original system. To build an aeroelastic ROM, it is possible to retain the first r-order POD modes ${{\boldsymbol \Psi} _r} = [{{\boldsymbol \psi} _1},{{\boldsymbol \psi} _2},\ldots,{{\boldsymbol \psi} _r}]$ while retaining most of the energy of the original system. By projecting the full-order series ${{\boldsymbol x}^{n \times 1}}$ on the r-order POD modes ${{\boldsymbol \Psi} _r} = [{{\boldsymbol \psi} _1},{{\boldsymbol \psi} _2},\ldots,{{\boldsymbol \psi} _r}]$ , we can reduce the full order system to a reduced r-order system

(10) \begin{align} \left\{ \begin{array}{l}{{\dot {\textbf x}}_r} = {\boldsymbol \Psi} _r^T{\textbf A_a}{{\boldsymbol \Psi} _r}{\textbf x_r} + {\boldsymbol \Psi} _r^T{\textbf B_a}\textbf y{\rm{ + }}{\boldsymbol \Psi} _r^T{\textbf L_a}{\textbf u_c}\\[5pt] {\textbf f^{ext}} = {\textbf C_a}{\textbf x_r}\end{array} \right.\end{align}

where ${\textbf A_a} = -\textbf A_0^{ - 1}\textbf H$ , ${\textbf B_a} = - \textbf A_0^{ - 1}\left[ {\textbf E + {\textbf {C G}}} \right]$ , ${\textbf L_a} = - \textbf A_0^{ - 1}\textbf L$ , ${\textbf C_a} = \textbf C{{\boldsymbol \Psi} _r}$ , $\textbf y = {\left[ {\dot {\textbf u}}\;{\textbf u} \right]^T}$ .

The above steps outline the process of generating an unsteady aerodynamics ROM in active flutter suppression. The aerodynamic force is jointly determined by the structural deformation and the deflection of control surface. The resulting aeroelastic model is obtained by coupling the structural dynamic equations, Eq. (3), with Eq. (10) for the ROM of the fluid.

2.2. LQR control design

A Linear Quadratic Regulator (LQR) regulator is used in current study for the active flutter control. The LQR assumes the availability of the entire state vector. In this work practice, some elements of the state vector may not be measurable, so sub-optimal control is generally used. The measured inputs are the displacement and velocity at wing tip where the sensors are located. First, the LQR state feedback design method states that

(11) \begin{align} {\textbf u_c} = - {\textbf K^*}{\textbf x_{ase}}\end{align}

where ${\textbf x_{ase}} = {\left[ {{\textbf x^r},{\rm{ }}\textbf y} \right]^T}$ . The suboptimal control design method assumes that only a few linear combinations of system states can be directly measured from the sensors. Let $\tilde {\textbf K} = {\textbf K^*}\textbf C_{ase}^T/({\textbf C_{ase}}\textbf C_{ase}^T)$ , ${\textbf C_{ase}} = \left[ {\textbf 0,{\rm{ }}{\textbf C_a}} \right]$ . The output vector feedback control law is formulated as

(12) \begin{align} {\textbf u_c} = - \tilde {\textbf K}{\textbf C_{ase}}{\textbf x_{ase}} = - \tilde {\textbf {Ky}}\end{align}

For a conventional LQR method, aiming to minimise the cost function

(13) \begin{align} J = \sum\limits_{k = 0}^\infty {\left\{ {{\textbf x^T}\textbf {Qx} + \textbf u_c^T\textbf R{\textbf u_c}} \right\}} \end{align}

where Q is a semi-positive definite weighting symmetric matrix, and R is a positive definite weighting symmetric matrix. The optimal problem of minimising the objective function J can be converted to the solution of:

(14) \begin{align} \begin{array}{l}\tilde {\textbf K} = {\textbf R^{ - 1}}\textbf {BP}\\[5pt] \textbf{PA} + {\textbf A^T}\textbf P - \textbf {PB}{\textbf R^{ - 1}}{\textbf B^T}\textbf P + \textbf Q = 0\end{array}\end{align}

where P is the solution of the Riccati equation, which can be calculated by many iterative algorithms. Like the conventional LQR design method, the optimal state feedback gains were obtained in this work using MATLAB /SIMULINK.

In the standard CFD-based POD/ROM construction procedure for active flutter suppression, when the structure displacement is obtained, the control law which designed by LQR is used to control the deflection of control surface at each time step, and add the influence of the control surface deflection to the aerodynamic force into the aeroelastic system to suppress flutter.

3.0 Efficient reduced-order framework

3.1. Structural dynamic reanalysis method

It is an important aspect to recalculate the eigenvalues of structure in some structural optimisation problems. For large-scale eigenvalue problem, one main obstacle is the high computational cost, and they require large computer memories and computing time. To overcome the obstacle, a number of methods including direct methods and approximate methods, have been proposed to ease the eigenvalue reanalysis [Reference Song, Chen and Sun40]. Direct methods are commonly based on the Sherman-Morrison-Woodbury (SMW) formula [Reference Akgün, Garcelon and Haftka41, Reference Huang, Wang and Li42]. But its application is limited by the rank of the modification, and also suffer from high computational costs. Different from direct methods, approximate methods aim at obtaining the modal characteristics of modified structures [Reference Beliveau, Cogan, Lallement and Ayer43Reference Kirsch, Bogomolni and Sheinman45], in which the mode shapes are approximated as a linear combination of basic mode shapes. The advantage in doing so is a significant reduction of computational costs, and the applicability is extended for high-rank modifications of structures. The extended Kirsch combined method [Reference Chen and Yang44, Reference Chen, Yang and Lian46] is an efficient method for the case of large modifications of the structural parameters. More details relevant to this work are given next.

Consider an initial structure with stiffness matrix ${\textbf K_0}$ and mass matrix ${\textbf M_0}$ . The corresponding mode shapes ${\boldsymbol{\phi}} _0^i$ and modal frequencies $\lambda _0^i$ are calculated by solving the equations:

(15) \begin{align} {\textbf K_0}{\boldsymbol{\phi}} _0^i = \lambda _0^i{\textbf M_0}{\boldsymbol \phi} _0^i\end{align}

${\textbf K_0}$ and ${\textbf M_0}$ are perturbed into the form ${\textbf K_0} + \Delta \textbf K$ and ${\textbf M_0} + \Delta \textbf M$ , $\textbf K = {\textbf K_0} + \Delta \textbf K$ , $\textbf M = {\textbf M_0} + \Delta \textbf M$ , in which $\Delta \textbf K$ and $\Delta \textbf M$ are the perturbations in the stiffness matrix and mass matrix, respectively. The eigenvalue problem for the modified structure becomes:

(16) \begin{align} {\textbf K}{{\boldsymbol{\phi}} ^i} = {\lambda ^i}{\textbf M}{{\boldsymbol{\phi}} ^i}\end{align}

where ${\lambda ^i}$ and ${{\boldsymbol{\phi}} ^i}$ are the i-th eigenvalue and eigenvector of the modified structure, respectively.

Extended Kirsch combined method use the second-order eigenvector terms [Reference Chen, Wu and Yang47] as the basis vectors in the following modal shape reduced basis:

(17) \begin{align} {{\boldsymbol{\phi}} ^i} = {\boldsymbol{\phi}} _B^i{\textbf z^i},\quad {\textbf z^i} = {(z_0^i,{\rm{ }}z_1^i,{\rm{ }}z_2^i)^T} \in {{\textbf R}^{3 \times 1}}\end{align}

where

(18) \begin{align} {\boldsymbol{\phi}} _B^i = \left[ {{\boldsymbol{\phi}} _0^i,{\rm{ }}{\boldsymbol{\phi}} _1^i,{\rm{ }}{\boldsymbol{\phi}} _2^i} \right]\end{align}
(19) \begin{align} \lambda _1^i = {({\boldsymbol{\phi}} _0^i)^T}\Delta {{\textbf K}{\boldsymbol{\phi}}} _0^i - \lambda _0^i{({\boldsymbol{\phi}} _0^i)^T}\Delta {\textbf M}{\boldsymbol{\phi}} _0^i{\rm{ }}\end{align}
(20) \begin{align} {\boldsymbol{\phi}} _1^i = \sum\limits_{s = 1,s \ne i}^n {\frac{1}{{\lambda _0^i - \lambda _0^s}}} \left[{({\boldsymbol{\phi}} _0^s)^T}(\Delta {\textbf K} - \lambda _0^i\Delta {\textbf M}){\boldsymbol{\phi}} _0^i\right]{\boldsymbol{\phi}} _0^s - \frac{1}{2}\left[{({\boldsymbol{\phi}} _0^i)^T}\Delta {\textbf M}{\boldsymbol{\phi}} _0^i\right]{\boldsymbol{\phi}} _0^i = {{\boldsymbol{\phi}} _0}{\textbf Z}_1^i\end{align}
(21) \begin{align} \lambda _2^i & = {({\boldsymbol{\phi}} _0^i)^ T}\Delta {\textbf K}{\boldsymbol{\phi}} _1^i - \lambda _0^i{({\boldsymbol{\phi}} _0^i)^T}\Delta {\textbf M}{\boldsymbol{\phi}} _1^i - \lambda _1^i{({\boldsymbol{\phi}} _0^i)^T}{{\textbf M}_0}{\boldsymbol{\phi}} _1^i - \lambda _1^i{({\boldsymbol{\phi}} _0^i)^T}\Delta {\textbf M}{\boldsymbol{\phi}} _0^i{\rm{ }} \nonumber \\ i,s & = 1,2,3 \ldots ..,{\rm{ }}\;number{\rm{ }}\;of{\rm{ }}\;modes\end{align}
(22) \begin{align} {\boldsymbol{\phi}} _2^i & = \sum\limits_{s = 1,s \ne i}^n {\frac{1}{{\lambda _0^i - \lambda _0^s}}} \left[{({\boldsymbol{\phi}} _0^s)^T}(\Delta {\textbf K} - \lambda _0^i\Delta {\textbf M}){\boldsymbol{\phi}} _1^i - \lambda _1^i{({\boldsymbol{\phi}} _0^s)^T}({{\textbf M}_0}{\boldsymbol{\phi}} _1^i + \Delta {\textbf M}{\boldsymbol{\phi}} _0^i)\right]{\boldsymbol{\phi}} _0^s \nonumber \\ & \quad - \frac{1}{2}\left[{({\boldsymbol{\phi}} _0^i)^T}\Delta {\textbf M}{\boldsymbol{\phi}} _1^i + {({\boldsymbol{\phi}} _1^i)^T}({\textbf M_0}{\boldsymbol{\phi}} _1^i + \Delta {\textbf M}{\boldsymbol{\phi}} _0^i)\right]{\boldsymbol{\phi}} _0^i = {{\boldsymbol{\phi}} _0}{\textbf Z}_2^i\end{align}

where $\lambda _1^i$ and $\lambda _2^i$ are the i-th first-order and second-order eigenvalue of the modified structure, $\boldsymbol{\phi} _1^i$ and $\boldsymbol{\phi} _2^i$ are the i-th first-order and second-order eigenvector of the modified structure, respectively. The coefficient vector, ${\textbf z^i}$ , contains three unknowns (for a second-order perturbation). Substituting Eq. (20) and Eq. (22) into Eq. (17), ${\boldsymbol \Phi} $ can be written as

(23) \begin{align} {\boldsymbol \Phi} = \left[ {\begin{array}{*{20}{c}}{{\boldsymbol{\phi}} _0^{1}\ldots{\boldsymbol{\phi}} _0^i}\\{\boldsymbol{0}}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{{\boldsymbol{\phi}} _0^1\ldots{\boldsymbol{\phi}} _0^i}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{\boldsymbol{0}}\\{{\boldsymbol{\phi}} _0^1\ldots{\boldsymbol{\phi}} _0^i}\end{array}} \right]{\left[ {\begin{array}{{c}}{{\rm {\textbf I}}_0^1\ldots{\rm {\textbf I}}_0^i}\\{\boldsymbol{0}}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{{\textbf Z}_1^1 \ldots {\textbf Z}_1^i}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{\boldsymbol{0}}\\{{\textbf Z}_2^1 \ldots {\textbf Z}_2^i}\end{array}} \right]^T}\left[ {\begin{array}{{c}}{{{\textbf z}^1}}\\{\boldsymbol{0}}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{{{\textbf z}^2}}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{\boldsymbol{0}}\\{{{\textbf z}^3}}\end{array}} \right] = {\Phi _0}{\textbf Z}\end{align}

where

(24) \begin{align} {\textbf Z} = {\left[ {\begin{array}{*{20}{c}}{{\textbf I}_0^1 \ldots {\textbf I}_0^i}\\{\boldsymbol{0}}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{{\textbf Z}_1^1 \ldots {\textbf Z}_1^i}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{\boldsymbol{0}}\\{{\textbf Z}_2^1 \ldots {\textbf Z}_2^i}\end{array}} \right]^T}\left[ {\begin{array}{{c}}{{{\textbf z}^1}}\\{\boldsymbol{0}}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{{{\textbf z}^2}}\\{\boldsymbol{0}}\end{array}\begin{array}{{c}}{\boldsymbol{0}}\\{\boldsymbol{0}}\\{{{\textbf z}^3}}\end{array}} \right]\end{align}

Substituting Eq. (17) into the modified analysis equations Eq. (16), and premultiplying by ${({\boldsymbol{\phi}} _B^i)^T}$ , one obtains:

(25) \begin{align} {({\boldsymbol{\phi}} _B^i)^T}({\textbf K_0} + \Delta {\textbf K}){\boldsymbol{\phi}} _B^i{{\textbf z}^i} = {\lambda ^i}{({\boldsymbol{\phi}} _B^i)^T}({\textbf M_0} + \Delta \textbf M){\boldsymbol{\phi}} _B^i{z^i}\end{align}

Introducing the notation

(26) \begin{align} {\textbf K}_R^i = {({\boldsymbol{\phi}} _B^i)^T}({{\textbf K}_0} + \Delta {\textbf K}{\boldsymbol{\phi}})\end{align}
(27) \begin{align} {\textbf M}_R^i = {({\boldsymbol{\phi}} _B^i)^T}({{\textbf M}_0} + \Delta {\textbf M}){\boldsymbol{\phi}} _B^i\end{align}

and substituting Eq. (26) and (27) into Eq. (25), we can obtain a set of ( $3 \times 3$ ) matrix equation

(28) \begin{align} {\textbf K}_R^i{{\textbf z}^i} = {\lambda ^i}{\textbf M}_R^i{{\textbf z}^i}\end{align}

Thus, the coefficient vector ${{\textbf z}^i}$ is evaluated from Eq. (28). The i-th eigenvector of the modified structure is obtained by substituting ${{\textbf z}^i}$ into Eq. (17) and Z is obtained by substituting ${{\textbf z}^i}$ into Eq. (24). It needs to be noted that the first-order eigenvector of Eq. (28) is only used [Reference Chen and Yang44].

Finally, the i-th eigenvalue of the modified structure, $\lambda _K^i$ , is computed using the Rayleigh quotient:

(29) \begin{align} \lambda _K^i = \frac{{{{({{\boldsymbol{\phi}} ^i})}^T}({{\textbf K}_0} + \Delta {\textbf K}){{\boldsymbol{\phi}} ^i}}}{{{{({{\boldsymbol{\phi}} ^i})}^T}{{\textbf M}_0}{{\boldsymbol{\phi}} ^i}}}\end{align}

3.2. Efficient reduced-order framework

In section 2.1, the unsteady aerodynamics ROM should be regenerated for each structural modification. Therefore, it cannot be used for active/passive hybrid flutter suppression. In the proposed efficient reduced-order framework, the initial structural mode shapes ${{\boldsymbol \Phi} _0}$ is taken as the basic mode shapes for basic POD/ROM construction. For the modified structure, the physical displacement ${\boldsymbol \Theta} $ of the wing can be written as

(30) \begin{align} {\boldsymbol \Theta} = {\boldsymbol \Phi} {\textbf u}\end{align}

Substituting Eq. (23) into Eq. (30), the physical displacement of the wing can also be written as

(31) \begin{align} {\boldsymbol \Theta} = {{\boldsymbol \Phi} _0}(\textbf {Zu})\end{align}

and ${{\textbf u}_b} = {\textbf {Zu}}$ , ${{\textbf u}_b}$ is artificially defined as basic generalised displacements.

For the modified structure in the efficient reduced-order framework, the matrices G, E and C of the linearised flow solver, Eq. (1), substituting the relation ${{\textbf u}_b} = {\textbf {Zu}}$ , the matrices may be rewritten in terms of the vector of basic generalised displacements:

\begin{align*}{\textbf G} = \frac{{\partial {\textbf F}}}{{\partial {\textbf u}}}({{\textbf w}_0},{{\textbf u}_0},{\dot {\textbf u}_0}) = \frac{{\partial {\textbf F}}}{{{{\textbf Z}^{ - 1}}\partial {{\textbf u}_b}}}({{\textbf w}_0},{{\textbf u}_0},{\dot {\textbf u}_0}) = {\textbf Z}{{\textbf G}_b}\end{align*}
(32) \begin{align} {\textbf E} = \frac{{\partial {\textbf A}}}{{\partial {\textbf u}}}{{\textbf w}_0} = \frac{{\partial {\textbf A}}}{{{{\textbf Z}^{ - 1}}\partial {{\textbf u}_b}}}{{\textbf w}_0} = {\textbf Z}{{\textbf E}_b}\end{align}
\begin{align*}{\textbf C} = \frac{{\partial {\textbf F}}}{{\partial \dot {\textbf u}}}({{\textbf w}_0},{{\textbf u}_0},{\dot {\textbf u}_0}) = \frac{{\partial {\textbf F}}}{{{{\textbf Z}^{ - 1}}\partial {{\dot {\textbf u}}_b}}}({{\textbf w}_0},{{\textbf u}_0},{\dot {\textbf u}_0}) = {\textbf Z}{{\textbf C}_b}\end{align*}

where $\textbf{G}_{b}$ , $\textbf{E}_{b}$ , $\textbf{C}_{b}$ are the first order terms in a Taylor series expansion of the basis reduced r-order aeroelastic model. In flutter suppression, since the elastic influence of the control surface can be ignored, the gradient of the numerical flux function with respect to control surface deflection $\textbf{L}_{a}$ is not changed. Now, the ROM in the proposed efficient reduced-order framework for active/passive hybrid flutter suppression is written as

(33) \begin{align} \left\{ \begin{array}{l}{{\dot {\textbf x}}_r} = {\boldsymbol \Psi} _r^T{\textbf A}{{\boldsymbol \Psi} _r}{{\textbf x}_r} + {\boldsymbol \Psi} _r^T{\textbf Z}{{\textbf B}_b}{\textbf y} + {\boldsymbol \Psi} _r^T{{\textbf L}_a}{{\textbf u}_c}\\[5pt] {{\textbf f}^{ext}} = {\textbf Z}{{\textbf C}_b}{{\boldsymbol \Psi} _r}{{\textbf x}_r}\end{array} \right.\end{align}

where ${\textbf A} = - {\textbf A}_0^{ - 1}{\textbf H}$ , ${{\textbf B}_b} = - {\textbf A}_0^{ - 1}\left[ {{{\textbf E}_b} + {{\textbf C}_b}\; {{\textbf G}_b}} \right]$ , ${\textbf y} = {\left[ {\dot {\textbf u}}\;{\textbf u} \right]^T}$ , ${{\textbf C}_b} = {{\textbf P}_b}$ .

The structural dynamic equations of modified structure in the efficient framework are:

(34) \begin{align} \bar {\textbf M}\ddot {\textbf u} + \bar {\textbf {Ku}} = {{\textbf f}^{ext}}\end{align}

here, $\bar {\textbf M} = {{\textbf Z}^T}{\boldsymbol \Phi} _0^T({{\textbf M}_0} + \Delta {\textbf M}){\Phi _0}{\textbf Z}$ , $\bar {\textbf K} = {{\textbf Z}^T}{\boldsymbol \Phi} _0^T({{\textbf K}_0} + \Delta {\textbf K}){{\boldsymbol \Phi} _0}{\textbf Z}$ . In addition, when $\Delta {\textbf K} = {\boldsymbol{0}}$ and $\Delta {\textbf M} = {\boldsymbol{0}}$ , Z is the unit matrix, Eq. (33) and Eq. (34) are equivalent to Eq. (10) and Eq. (3), respectively. For active/passive hybrid flutter suppression, adding a mass balance and changing the structural stiffness were made by varying the structural parameters, both the structural model and the time-consuming CFD-based ROM have to be reconstructed to ensure accuracy. When the aeroelastic response is obtained, the control surface rotation is updated at each time step by control law which is designed by LQR and then add the influence of deflection of control surface to the aerodynamic force into the aeroelastic system to achieve the purpose of flutter suppression. The aeroelastic response of the modified structure is then obtained by the proposed efficient reduced-order framework, without the expensive time-consuming reconstruction procedure of the new POD basis. Therefore, the proposed framework can significantly reduce the computational cost. The proposed framework provides a potential powerful tool for active/passive hybrid flutter suppression that would be convenient to design more effective controllers to improve aircraft fatigue life and enhance aircraft performance.

4. Numerical results and discussion

4.1. POD/ROM solver validation

The accuracy of POD/ROM solver is validated by the AGARD 445.6 aeroelastic wing model [Reference Yates48, Reference Silva, Chwalowski and Perry49]. The aerofoil section is a NACA 65A004. The density of wing material is 381.98 kg/m3. The elastic modulus in the spanwise direction and the chordwise direction is 3.151Gpa and 0.416GPa, respectively. The shear modulus is 0.4392GPa, and the Poisson’s ratio is 0.31 [Reference Zhong and Xu50]. The structural model based on shell elements, shown in Fig. 1(a), consists of 231 nodes and 200 elements, and the thickness distribution of the wing was governed by aereofoil shape. For the fluid model, a multi-block structured mesh was employed, as shown in Fig. 1(b). The spatial convergence of the CFD mesh was analysed in Zhou et al. [Reference Qiang, Chen, Li and da Ronch51], which reported a good agreement of the results from both the medium and fine grids. In this paper, the total number of grid points on the medium grid, herein used, is 223,146 ( $99 \times 49 \times 46$ ).

Figure 1. AGARD 445.6 wing: (a) structural model and trailing-edge control surface, and (b) surface CFD mesh.

Figure 2 shows the flutter boundary computed by the full-order CFD/CSD aeroelastic model and POD/ROM solver, experimental data [Reference Yates48] are also included. The agreement between the POD/ROM, CFD/CSD aeroelastic model and experimental data is good for all Mach numbers considered (0.499 to 1.141), including the well-known transonic dip of the flutter speed. The accuracy of full-order CFD/CSD aeroelastic model and POD/ROM solver have been evaluated over the years in a number of aeroelastic studies [Reference Zhou, Chen and Li16, Reference Chen, Li, Zhou, da Ronch and Li34, Reference Li, da Ronch, Chen and Li36, Reference Zhou, Li, da Ronch, Chen and Li52, Reference Chen, Sun and Li53].The good agreement indicates that the POD/ROM solver has near the same accuracy as the nonlinear CFD model at small or moderate perturbation.

Figure 2. AGARD 445.6 wing flutter boundary.

4.2. Accuracy evaluation of the efficient aeroelastic CFD-based POD/ROM

The accuracy and efficiency of the efficient CFD-based POD/ROM has been successfully verified in [Reference Chen, Li, Zhou, da Ronch and Li34]. The AGARD 445.6 wing structural model is divided into four sections along spanwise direction, as shown in Fig. 1(a). Each section of the identical material properties, but material properties at different sections may be different. The material properties are assumed to vary as

  1. Section 1, ${E_1} = (1 + 3\varepsilon ){E_0}$ , ${\rho _1} = (1 + 3\varepsilon ){\rho _0}$

  2. Section 2, ${E_2} = (1 + 2\varepsilon ){E_0}$ , ${\rho _2} = (1 + 2\varepsilon ){\rho _0}$

  3. Section 3, ${E_3} = (1 + \varepsilon ){E_0}$ , ${\rho _3} = (1 + \varepsilon ){\rho _0}$

  4. Section 4, ${E_4} = {E_0}$ , ${\rho _4} = {\rho _0}$

where ${E_0}$ and ${\rho _0}$ are the Young’s modulus and density of the original wing structure, respectively. It is worth observing that the choice above is used for demonstration purposes and shall not be taken here as a limiting case of the present method.

Three cases for the improved AGARD 445.6 wing are considered: $\varepsilon$ = 1/12, 1/6, 1/3. For conciseness, the flutter speed at Ma = 0.960 (air density is 0.06341 kg/m3) and angle-of-attack AOA = 0deg (the flutter speed is 294m/s for standard AGARD 445.6 wing model) using the efficient aeroelastic CFD-based POD/ROM is presented in Table 1. Reference data is taken from the exact model in which both the structural model and the CFD-based POD/ROM have to be reconstructed for each modified structure. Although the flutter speed error increase for increasing level of the structural modifications, it is worth observing that the maximum difference is limited to about 2%. It indicates that the efficient CFD-based POD/ROM has good accuracy for flutter speeds prediction, even for very large range ( $\varepsilon$ = 1/3) of structural parameter variation.

Table 1. Flutter speed obtained by exact model and the efficient aeroelastic CFD-based POD/ROM

4.3. Accuracy evaluation of the efficient reduced-order framework

After evaluating the accuracy of POD/ROM solver and the efficient aeroelastic CFD-based POD/ROM, the whole proposed efficient reduced-order framework for active/passive hybrid flutter suppression will be evaluated using an improved AGARD 445.6 wing model in this section.

The structural model of the improved AGARD 445.6 wing were modified to accommodate a trailing-edge control surface at the tip of the wing, with dimensions equal to 20% of the wing span and 30% of the wing chord, is shown in Fig. 1(a). Two sensors were added at the wing tip to measure the displacements and velocities of the wing. The motion of the control surface is treated as an additional mode shape. The mode shape for the control surface deflection on CFD surface, is shown in Fig. 3. It can be seen that the control surface deflection on the CFD surface keeps the form of control surface deflection of structural model. It needs to be noted that rotations of the control surface are small, so modeling gaps can be ignored [Reference Zhou, Li, da Ronch, Chen and Li52].

Figure 3. The mode shape for the control surface deflection on the CFD surface (in red).

For simplifying the description, only the aeroelastic responses at Ma = 0.960 (air density is 0.06341Kg/m3) and AOA = 0deg are presented in this paper. For the three cases ( $\varepsilon$ = 1/12, 1/6, 1/3), the freestream speed is set to 340m/s, 360m/s and 395m/s, respectively, which are about 9% higher than the flutter speed. The weight matrices, Q 1 = 0.1I 60 × 60 and R = 10, were used for the LQR synthesis. Figures 46 show the aeroelastic responses obtain by proposed efficient reduced-order framework. Reference data are representative of the exact system, where the aeroelastic CFD-based POD/ROM has to be reconstructed for each modified structure. For all cases, although the discrepancy between the two methods increases for increasing level of the structural modifications, the good agreement between the exact system and the proposed framework are founded. The good agreement indicate that the proposed efficient reduced-order framework has good predictive capability for aeroelastic response. As can be seen from Figs. 46, at same freestream speed condition, the aeroelastic response of passive control alone continuously diverges that the system is unstable, and the aeroelastic responses of the hybrid control is convergent that the system is stable. It indicates that the proposed framework is ideally suitable for active/passive hybrid flutter suppression to design more effective controllers with good accuracy. In addition, Figs. 4-6 also indicate that the efficient aeroelastic CFD-based POD/ROM has good accuracy for aeroelastic responses prediction.

Figure 4. Comparison of the aeroelastic response for $\varepsilon$ = 1/12 (V =340m/s).

Figure 5. Comparison of the aeroelastic response for $\varepsilon$ = 1/6 (V =365m/s).

Figure 6. Comparison of the aeroelastic response for $\varepsilon$ = 1/3 (V =395m/s).

To further demonstrate the predictive capability of the proposed efficient reduced-order framework, the fifth modal coordinate related to the rotation of the trailing-edge control surface is presented in Fig. 7. The reference data are taken from the exact system. It can be seen that in all cases, the rotation of the trailing-edge control surface is consistent. It indicates again that the proposed framework has good accuracy for predict aeroelastic response, without reconstructing a new set of CFD-based POD/ROM for modified structure. All the above comparisons indicate that the proposed efficient reduced-order framework can capture the aeroelastic response with good accuracy, and ideally suitable for active/passive hybrid flutter suppression to design more effective controllers with high accuracy. It should be noted that the proposed framework would work well for small to moderate deformations but not for large variations.

Figure 7. Trailing-edge control surface rotation.

Finally, a quantification of flutter speed errors for the proposed efficient reduced-order framework and exact system is presented in Table 2. It is not unexpected that the accuracy of the proposed framework degrades with the increasing structural modifications. However, it is worth observing that the maximum difference is limited to about 2.5%. For all cases, compared with passive control alone, the flutter speed increased by 36.3%, 36.0% and 37.4%, respectively. Compared with active control alone [Reference Zhou, Li, da Ronch, Chen and Li52], the flutter speed increased by 7.3%, 13.9% and 27.3%, respectively. It is indicate that compared with passive flutter suppression alone [Reference Chen, Li, Zhou, da Ronch and Li34] and active flutter suppression alone [Reference Zhou, Li, da Ronch, Chen and Li52], the proposed framework can greatly expand the flutter boundary. All the above comparisons indicate that the proposed efficient reduced-order framework is very suitable for active/passive hybrid flutter suppression to design more effective controllers with good accuracy.

Table 2. Flutter speed obtained by exact system and the efficient reduced-order framework

4.4. Efficiency evaluation of the efficient reduced-order framework

The computational efficiency is one of the most important criteria of the proposed efficient reduced-order framework. All analyses were performed on a Windows 10 system PC with Intel® Core(TM) i7-9700K CPU (3.60 GHz, 8 cores, but only one core used) and 32 GB RAM.

For the exact system, each modified structure need regenerate a set of POD modes, requiring about 16h per configuration. For the three cases, required about 48h. In contract, the set of POD modes is generated only once for the proposed efficient reduced-order framework. For the proposed efficient reduced-order framework, it only requires about 16h. It is reasonable to consider about 100 cases and 20 values of the freestream dynamic pressure to assess the aeroelastic stability. With information reported in Table 3, the exact system would require over 1,600 CPU hours (this consists of 16h × 100 and 2.69s × 20 × 100, totaling 1,609.49h), whereas the proposed framework just require over 17 CPU hours (16h × 1 and 2.92s × 20 × 100, totaling 17.62h). It is obvious that the expected speed-up is proportional to the number of modified structural models, the computational advantage of the efficient reduced order framework becomes more significant the more the number of modified structural model is considered. All the demonstration cases indicated that the proposed efficient reduced-order framework has high computational efficiency and can be applied for active/passive hybrid flutter suppression to design more effective controllers.

Table 3. Computational cost of the exact system and the efficient reduced order framework

5.0 Conclusions

In modern aircraft design, in order to design more effective controllers for flutter suppression to improve fatigue life and enhance performance, an active/passive hybrid control method could be beneficial. The CFD-based ROMs have been successfully applied to active flutter suppression with high accuracy and efficiency. We have developed and implemented an efficient aeroelastic CFD-based POD/ROM, which has been successfully applied to passive flutter suppression, such as aeroelastic tailoring and aeroelastic structural optimisation. In order to design more effective controllers with high efficiency, an efficient reduced-order framework for active/passive hybrid flutter suppression is proposed by extending the efficient aeroelastic CFD-based POD/ROM in this research effort. It creates a potential efficient tool for active/passive hybrid flutter suppression, which will allow the reduction of overall cost of the aircraft design.

The accuracy and efficiency of the proposed efficient reduced-order framework for active/passive hybrid flutter suppression was demonstrated and evaluated by an improved AGARD 445.6 wing model. Firstly, the accuracy of the POD/ROM solver and efficient aeroelastic CFD-based POD/ROM were validated. Then, the accuracy of the whole proposed framework was evaluated by comparing the aeroelastic responses including the generalised displacements and the rotation of the trailing-edge control surface. The good agreements of the numerical results show that the proposed efficient reduced-order framework can capture the aeroelastic response with high accuracy. The flutter speeds were also compared and although the flutter speed errors increase as the level of structural modifications, the max error is still less than 2.5%. The computational efficiency of the proposed framework, which generates only once the set of POD modes, the computational cost is reduced obviously. It is also obvious that the efficiency improvement is proportional to the number of modified structural models.

The greatest advantage of the proposed efficient reduced-order framework is not only has good predictive capability for aeroelastic response, but also offers a great reduction of the computational cost. All above research indicate that the proposed framework has the ability to be applied to active/passive hybrid flutter suppression. The proposed efficient reduced-order framework provides a potential powerful tool for active/passive hybrid flutter suppression that would be convenient to design more effective controllers with high efficiency and may have the potential to reduce of overall cost of aircraft design.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (No. 51775518, 11672225, 11511130053 and 11272005), the National Program on Key Research Projects (No. MJ-2015-F-010), the Shannxi Province Natural Science Foundation (No. 2016JM1007), and the funds for the Central Universities (2014XJJ0126).

Conflict of interest statement

There is no conflict of interest.

References

Kuzmina, S., Amiryants, G., Schweiger, J., Cooper, J., Amprikidis, M. and Sensberg, O. Review and outlook on active and passive aeroelastic design concepts for future aircraft. ICAS 2002 Congress, 2002, pp 813.Google Scholar
Silva, G.C., Silvestre, F.V.J., Donadon, M.C.V., Santos, O.S., Guimar, ES, Neto, A.N.B., da Silva, R.G.A., Versiani, T.D.S.S., Gonzalez, P.J. and Bertolin, R.M. Active and passive control for acceleration reduction of an aeroelastic typical wing section. J. Vib. Control., 2018, 24, pp 26732687.CrossRefGoogle Scholar
Liu, Y. Comparison of the passive and active control gust alleviation of a flying-wing aircraft. J. Eng., 2019, 2019, pp 421423.CrossRefGoogle Scholar
Reich, G.W. Design and modeling of an active aeroelastic wing, Massachusetts Institute of Technology, 1994, Cambridge, MA.Google Scholar
Gregory, R.H., Mart, L.C. and Crawley, E. An active aeroelastic wing model for vibration and flutter suppression. 36th Structures, Structural Dynamics and Materials Conference 1995.Google Scholar
Le, O.L., de Lima, A., Donadon, M. and Cunha-Filho, A. Dynamic and aeroelastic behavior of composite plates with multimode resonant shunted piezoceramics in series. Compos. Struct., 2016, 153, pp 815824.CrossRefGoogle Scholar
Barzegari, M.M., Dardel, M. and Fathi, A. Control of aeroelastic characteristics of cantilever wing with smart materials using state-dependent Riccati equation method. J. Intell. Mater. Syst. Struct., 2015, 26, pp 9881005.CrossRefGoogle Scholar
de Breuker, R., Binder, S. and Wildschek, A. Combined Active and Passive Loads Alleviation through Aeroelastic Tailoring and Control Surface/Control System Optimization. 2018 AIAA Aerospace Sciences Meeting 2018, 0764.Google Scholar
Schweiger, J., Krammer, J. and Coetzee, E. MDO application for active flexible aircraft design. 7th AIAA/USAF/NASA/ISSMO Symposium on Multidisciplinary Analysis and Optimization 1998, 4835.CrossRefGoogle Scholar
Weisshaar, T.A. and Duke, D.K. Induced drag reduction using aeroelastic tailoring with adaptive control surfaces. J. Aircr., 2006, 43, pp 157164.CrossRefGoogle Scholar
Handojo, V., Lancelot, P. and de Breuker, R. Implementation of Active and Passive Load Alleviation Methods on a Generic mid-Range Aircraft Configuration. 2018 Multidisciplinary Analysis and Optimization Conference 2018, 3573.Google Scholar
Chen, Z., Zhao, Y. and Huang, R. Parametric reduced-order modeling of unsteady aerodynamics for hypersonic vehicles. Aerosp. Sci. Technol., 2019, 87, pp 114.CrossRefGoogle Scholar
Kou, J. and Zhang, W. A hybrid reduced-order framework for complex aeroelastic simulations. Aerosp. Sci. Technol., 2019, 84, pp 880–94.CrossRefGoogle Scholar
Zhang, W., Lv, Z., Diwu, Q. and Zhong, H. A flutter prediction method with low cost and low risk from test data. Aerosp. Sci. Technol., 2019, 86, pp 542557.CrossRefGoogle Scholar
Lu, K., Jin, Y., Chen, Y., Yang, Y., Hou, L., Zhang, Z., Li, Z. and Fu, C. Review for order reduction based on proper orthogonal decomposition and outlooks of applications in mechanical systems. Mech. Syst. Signal Process., 2019, 123, pp 264297.CrossRefGoogle Scholar
Zhou, Q., Chen, G. and Li, Y. A reduced order model based on block Arnoldi method for aeroelastic system. Int. J. Appl. Mech., 2014, 6, pp 1450069.CrossRefGoogle Scholar
Cao, C., Nie, C., Pan, S., Cai, J. and Qu, K. A constrained reduced-order method for fast prediction of steady hypersonic flows. Aerosp. Sci. Technol., 2019, 91, pp 679690.CrossRefGoogle Scholar
Dowell, E.H., Thomas, J.P. and Hall, K.C. Transonic limit cycle oscillation analysis using reduced order aerodynamic models. J. Fluids Struct., 2004, 19, pp 1727.CrossRefGoogle Scholar
Zhou, Q., Chen, G., da Ronch, A. and Li, Y. Reduced order unsteady aerodynamic model of a rigid aerofoil in gust encounters. Aerosp. Sci. Technol., 2017, 63, pp 203213.CrossRefGoogle Scholar
Wu, X., Zhang, W., Peng, X. and Wang, Z. Benchmark aerodynamic shape optimization with the POD-based CST airfoil parametric method. Aerosp. Sci. Technol., 2019, 84, pp 632640.CrossRefGoogle Scholar
Li, L.Z., Li, J.J., Zhang, J., Lu, K. and Yuan, M.N. Aerodynamic shape optimization by continually moving ROM. Aerosp. Sci. Technol., 2020, 99, p 105729.CrossRefGoogle Scholar
Zhou Li, L., Yang, M., Luo, X., Zhang, J. and Ni Yuan, M. A spline ROM of blade aerodynamic force to upstream wake. Aerosp. Sci. Technol., 2019, 84, pp 650660.Google Scholar
Sajadmanesh, S.M., Mojaddam, M., Mohseni, A. and Nikparto, A. Numerical identification of separation bubble in an ultra-high-lift turbine cascade using URANS simulation and proper orthogonal decomposition. Aerosp. Sci. Technol., 2019, 93, p 105329.CrossRefGoogle Scholar
Gang, C., Jian, S. and Yueming, L. Active flutter suppression control law design method based on balanced proper orthogonal decomposition reduced order model. Nonlinear Dyn., 2012, 70, pp 112.CrossRefGoogle Scholar
Song, H., Qian, J., Wang, Y., Pant, K., Chin, A.W., Brenner, M.J. Development of Aeroelastic and Aeroservoelastic Reduced Order Models for Active Structural Control. 56th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 2015, 2055.CrossRefGoogle Scholar
Nie, X., Yang, G. and Zhang, M. Investigation on transonic flutter active suppression with CFD-Based ROMs. Sci. China Phys. Mech. Astron., 2015, 58, pp 110.CrossRefGoogle Scholar
Liu, H., Huang, R., Zhao, Y. and Hu, H. Reduced-order modeling of unsteady aerodynamics for an elastic wing with control surfaces. J. Aerosp. Eng., 2017, 30, p 04016083.CrossRefGoogle Scholar
Hoseini, H. and Hodges, D.H. Flutter suppression for finite element modeling of damaged hale aircraft wings. 2018 AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 2018, 1205.Google Scholar
Patterson, R.P. and Friedmann, P.P. Application of a CFD-Based Surrogate Approach for Active Flow Control Modeling. AIAA Scitech 2019 Forum 2019, 1706.CrossRefGoogle Scholar
Fenwick, C., Jones, D. and Gaitonde, A. Consideration of ROM interpolation for aeroelastic design using structural modification. 18th AIAA Computational Fluid Dynamics Conference, Miami, FL 2007.CrossRefGoogle Scholar
Zhang, W., Chen, K. and Ye, Z. Unsteady aerodynamic reduced-order modeling of an aeroelastic wing using arbitrary mode shapes. J. Fluids Struct., 2015, 58, pp 254270.CrossRefGoogle Scholar
Winter, M., Heckmeier, F.M. and Breitsamter, C. CFD-based aeroelastic reduced-order modeling robust to structural parameter variations. Aerosp. Sci. Technol., 2017, 67, pp 1330.CrossRefGoogle Scholar
Li, K., Kou, J. and Zhang, W. Deep neural network for unsteady aerodynamic and aeroelastic modeling across multiple Mach numbers. Nonlinear Dyn., 2019, 96, pp 21572177.CrossRefGoogle Scholar
Chen, G., Li, D., Zhou, Q., da Ronch, A. and Li, Y. Efficient aeroelastic reduced order model with global structural modifications. Aerosp. Sci. Technol., 2018, 76, pp 113.CrossRefGoogle Scholar
Li, D., Gong, C., da Ronch, A., Chen, G. and Li, Y. An efficient implementation of aeroelastic tailoring based on efficient computational fluid dynamics-based reduced order model. J. Fluids Struct., 2019, 84, pp 182198.CrossRefGoogle Scholar
Li, D., da Ronch, A., Chen, G. and Li, Y. Aeroelastic global structural optimization using an efficient CFD-based reduced order model. Aerosp. Sci. Technol., 2019, 94, p 105354.CrossRefGoogle Scholar
Hall, K.C., Thomas, J.P. and Dowell, E.H. Proper orthogonal decomposition technique for transonic unsteady aerodynamic flows. AIAA J., 2000, 38, pp 18531862.CrossRefGoogle Scholar
Lieu, T., Farhat, C., Lesoinne, M. POD-based aeroelastic analysis of a complete F-16 configuration: ROM adaptation and demonstration. 46th AIAA/ASME/ASCE/AHS/ASC structures, structural dynamics and materials conference 2005, 2295.CrossRefGoogle Scholar
Lieu, T. Adaptation of reduced order models for applications in aeroelasticity, University of Colorado at Boulder, 2004, Boulder, CO.Google Scholar
Song, Q., Chen, P. and Sun, S. An exact reanalysis algorithm for local non-topological high-rank structural modifications in finite element analysis. Comput. Struct., 2014, 143, pp 6072.CrossRefGoogle Scholar
Akgün, M.A., Garcelon, J.H. and Haftka, R.T. Fast exact linear and non-linear structural reanalysis and the Sherman–Morrison–Woodbury formulas. Int. J. Numer. Methods Eng., 2001, 50, pp 15871606.CrossRefGoogle Scholar
Huang, G., Wang, H. and Li, G. An exact reanalysis method for structures with local modifications. Struct. Multidiscipl. Optim., 2016, 54, pp 499509.CrossRefGoogle Scholar
Beliveau, J.-G., Cogan, S., Lallement, G. and Ayer, F. Iterative least-squares calculation for modal eigenvector sensitivity. AIAA J., 1996, 34, pp 385391.CrossRefGoogle Scholar
Chen, S.H. and Yang, X.W. Extended Kirsch combined method for eigenvalue reanalysis. AIAA J., 2000, 38, pp 927930.CrossRefGoogle Scholar
Kirsch, U., Bogomolni, M. and Sheinman, I. Efficient dynamic reanalysis of structures. J. Struct. Eng., 2007, 133, pp 440448.CrossRefGoogle Scholar
Chen, S., Yang, X. and Lian, H. Comparison of several eigenvalue reanalysis methods for modified structures. Struct. Multidiscipl. Optim., 2000, 20, pp 253259.CrossRefGoogle Scholar
Chen, S.H., Wu, X.M. and Yang, Z.J. Eigensolution reanalysis of modified structures using epsilon©\algorithm. Int. J. Numer. Methods Eng., 2006, 66, pp 21152130.CrossRefGoogle Scholar
Yates, E.C. Jr AGARD standard aeroelastic configurations for dynamic response. Candidate configuration I.-wing 445.6. 1987.Google Scholar
Silva, W.A., Chwalowski, P. and Perry, B. III Evaluation of linear, inviscid, viscous, and reduced-order modelling aeroelastic solutions of the AGARD 445.6 wing using root locus analysis. Int. J. Comput. Fluid Dyn., 2014, 28, pp 122139.CrossRefGoogle Scholar
Zhong, J. and Xu, Z. Coupled fluid structure analysis for wing 445.6 flutter using a fast dynamic mesh technology. Int. J. Comput. Fluid Dyn., 2016, 30, pp 531542.CrossRefGoogle Scholar
Qiang, Z., Chen, G., Li, Y. and da Ronch, A. Aeroelastic moving gust responses and alleviation based on CFD. AIAA Modeling and Simulation Technologies Conference 2016, 3837.CrossRefGoogle Scholar
Zhou, Q., Li, D.-F., da Ronch, A., Chen, G. and Li, Y.-M. Computational fluid dynamics-based transonic flutter suppression with control delay. J. Fluids Struct., 2016, 66, pp 183206.CrossRefGoogle Scholar
Chen, G., Sun, J. and Li, Y.-M. Adaptive reduced-order-model-based control-law design for active flutter suppression. J. Aircr., 2012, 49, pp 973980.CrossRefGoogle Scholar
Figure 0

Figure 1. AGARD 445.6 wing: (a) structural model and trailing-edge control surface, and (b) surface CFD mesh.

Figure 1

Figure 2. AGARD 445.6 wing flutter boundary.

Figure 2

Table 1. Flutter speed obtained by exact model and the efficient aeroelastic CFD-based POD/ROM

Figure 3

Figure 3. The mode shape for the control surface deflection on the CFD surface (in red).

Figure 4

Figure 4. Comparison of the aeroelastic response for $\varepsilon$ = 1/12 (V=340m/s).

Figure 5

Figure 5. Comparison of the aeroelastic response for $\varepsilon$ = 1/6 (V=365m/s).

Figure 6

Figure 6. Comparison of the aeroelastic response for $\varepsilon$ = 1/3 (V=395m/s).

Figure 7

Figure 7. Trailing-edge control surface rotation.

Figure 8

Table 2. Flutter speed obtained by exact system and the efficient reduced-order framework

Figure 9

Table 3. Computational cost of the exact system and the efficient reduced order framework