Output Feedback Modular Adaptive Control of a Nonlinear Prototypical Wing Section

The paper presents nonlinear adaptive control systems for the control of limit cycle oscillations of a prototypical wing section with structural nonlinearities using only output feedback. The chosen model describes the plunge and pitch motion of a wing. The model includes plunge and pitch nonlinearities, and has a single control surface for the purpose of control. Using a canonical representation of the aeroelastic system, a modular output feedback adaptive control system consisting of an input-to-state stabilizing controller and a passive identifier (an observer and adaptation law) is derived. In the closed-loop system, asymptotic stabilization of the pitch and plunge motion is accomplished. Simulation results show that the control system is effective in regulating the state vector to the origin in spite of large parameter uncertainties.     

Limit Cycle Oscillation and Orbital Stability in Aeroelastic Systems with Torsional Nonlinearity

The paper treats the question of the existence of limit cycleoscillations of prototypical aeroelastic wing sections with structuralnonlinearity using the describing function method. The chosen dynamicmodel describes the nonlinear plunge and pitch motion of a wing. Themodel includes an asymmetric structural nonlinearity in the pitchdegree-of-freedom. The dual-input describing functions of thenonlinearity are derived for the limit cycle analysis. Analyticalexpressions for the average value, and the amplitude and frequency ofoscillation of pitch and plunge responses are obtained. Based on ananalytical approach as well as the Nyquist criterion, stability of thelimit cycles is examined. Numerical results are presented for a set ofvalues of the flow velocities and the locations of the elastic axiswhich show that the predicted limit cycle oscillation amplitude andfrequency as well as the mean value are quite close to the actualvalues. Furthermore, for the chosen model with linear aerodynamics, itis seen that the amplitude of the pitch limit cycle oscillation does notalways increase with the flow velocity for certain elastic axislocations.     

Effects of multiple structural nonlinearities on limit cycle oscillation of missile control fin

Aeroelastic analyses are performed for a 2-D typical section model with multiple nonlinearities. The differences between a system with multiple nonlinearities in its pitch and plunge spring and a system with a single nonlinearity in its pitch are thoroughly investigated. The unsteady supersonic aerodynamic forces are calculated by the doublet point method (DPM). The iterative V-g method is used for a multiple-nonlinear aeroelastic analysis in the frequency domain and the freeplay nonlinearity is linearized using a describing function method. In the time domain, the DPM unsteady aerodynamic forces, which are based on a function of the reduced frequency, are approximated by the minimum state approximation method. Consequently, multiple structural nonlinearities in the 2-D typical wing section model are influenced by the pitch to plunge frequency ratio. This result is important in that it demonstrates that the flutter speed is closely connected with the frequency ratio, considering that both pitch and plunge nonlinearities result in a higher flutter speed boundary than a conventional aeroelastic system with only one pitch nonlinearity. Furthermore, the gap size of the freeplay affects the amplitude of the limit cycle oscillation (LCO) to gap size ratio.     

Nonlinear state feedback control design to eliminate subcritical limit cycle oscillations in aeroelastic systems

A strictly nonlinear state feedback control law is designed for an aeroelastic system to eliminate subcritical limit cycle oscillations. Numerical continuation techniques and harmonic balance methods are employed to generate analytical estimates of limit cycle oscillation commencement velocity and its sensitivity with respect to the introduced control parameters. The obtained estimates are used in a multiobjective optimization framework to generate optimal control parameters which maximize the limit cycle oscillation commencement velocity while minimizing the control cost. Numerical simulations are used to show that the assumed nonlinear state feedback law with the optimal control parameters successfully eliminates any existing subcritical limit cycle oscillations by converting it to supercritical limit cycle oscillations, thereby guaranteeing safe operation of the system in its flight envelope.     

Spectrographic analysis for the modal testing of nonlinear aeroelastic systems

The spectrograph is a signal-processing tool often used for the frequency domain analysis of time-varying signals. When the signal to be analyzed is a function of time, the spectrograph represents the frequency content of the signal as a sequence of power spectra that change with time. In this paper, the usefulness of the technique is demonstrated in its application to the analysis of the time history response of a nonlinear aeroelastic system. The aeroelastic system is modelled analytically as a two-dimensional, rigid airfoil section free to move in both the bending and pitching directions and possessing a rigid flap. The airfoil is mounted by torsional and translational springs attached at the elastic axis, and the flap is used to provide the forcing input to the system. The nonlinear system is obtained by introducing a freeplay type of nonlinearity in the pitch degree-of-freedom restoring moment. The airfoil is immersed in an aerodynamic flow environment, modelled using incompressible thin airfoil theory for unsteady oscillatory motion. The equations of motion are solved using a fourth-order Runge–Kutta numerical integration technique to provide time-history solutions of the response of the airfoil in the pitch and plunge directions. Time-histories are obtained for the nonlinear responses of the linear and nonlinear aeroelastic systems to a sine-sweep input. The time-histories are analyzed using the spectrographic technique, and the frequency content of the response is plotted directly as a function of the input frequency. Results show that the combination of the sine-sweep input with the spectrographic analysis permits a unique insight into the behavior of the nonlinear system with a minimum of testing. It is shown that the frequency of the nonlinear system response is a function of the input frequency and one other characteristic frequency that can be associated with the limit cycle oscillations of the same nonlinear system subject to a transient input.     

Low speed flutter and limit cycle oscillations of a two-degree-of-freedom flat plate in a wind tunnel

This paper explores the dynamical response of a two-degree-of-freedom flat plate undergoing classical coupled-mode flutter in a wind tunnel. Tests are performed at low Reynolds number (Re~2.5×10⁴), using an aeroelastic set-up that enables high amplitude pitch–plunge motion. Starting from rest and increasing the flow velocity, an unstable behaviour is first observed at the merging of frequencies: after a transient growth period the system enters a low amplitude limit-cycle oscillation regime with slowly varying amplitude. For higher velocity the system transitions to higher-amplitude and stable limit cycle oscillations (LCO) with amplitude increasing with the flow velocity. Decreasing the velocity from this upper LCO branch the system remains in stable self-sustained oscillations down to 85% of the critical velocity. Starting from rest, the system can also move toward a stable LCO regime if a significant perturbation is imposed. Those results show that both the flutter boundary and post-critical behaviour are affected by nonlinear mechanisms. They also suggest that nonlinear aerodynamic effects play a significant role.     

FLUTTER OF AN AIRFOIL WITH A CUBIC RESTORING FORCE

In this paper, the effect of a cubic structural restoring force on the flutter characteristics of a two-dimensional airfoil placed in an incompressible flow is investigated. The aeroelastic equations of motion are written as a system of eight first-order ordinary differential equations. Given the initial values of plunge and pitch displacements and their velocities, the system of equations is integrated numerically using a fourth order Runge-Kutta scheme. Results for soft and hard springs are presented for a pitch degree-of-freedom nonlinearity. The study shows the dependence of the divergence flutter boundary on initial conditions for a soft spring. For a hard spring, the nonlinear flutter boundary is independent of initial conditions for the spring constants considered. The flutter speed is identical to that for a linear spring. Divergent flutter is not encountered, but instead limit-cycle oscillation occurs for velocities greater than the flutter speed. The behaviour of the airfoil is also analysed using analytical techniques developed for nonlinear dynamical systems. The Hopf bifurcation point is determined analytically and the amplitude of the limit-cycle oscillation in post-Hopf bifurcation for a hard spring is predicted using an asymptotic theory. The frequency of the limit-cycle oscillation is estimated from an approximate method. Comparisons with numerical simulations are carried out and the accuracy of the approximate method is discussed. The analysis can readily be extended to study limit-cycle oscillation of airfoils with nonlinear polynomial spring forces in both plunge and pitch degrees of freedom.     

Reduced order nonlinear aeroelasticity of swept composite wings using compressible indicial unsteady aerodynamics

Nonlinear dynamic aeroelasticity of composite wings in compressible flows is investigated. To provide a reasonable model for the problem, the composite wing is modeled as a thin walled beam (TWB) with circumferentially asymmetric stiffness layup configuration. The structural model considers nonlinear strain displacement relations and a number of non-classical effects, such as transverse shear and warping inhibition. Geometrically nonlinear terms of up to third order are retained in the formulation. Unsteady aerodynamic loads are calculated according to a compressible model, described by indicial function approximations in the time domain. The aeroelastic system of equations is augmented by the differential equations governing the aerodynamics lag states to derive the final explicit form of the coupled fluid-structure equations of motion. The final nonlinear governing aeroelastic system of equations is solved using the eigenvectors of the linear structural equations of motion to approximate the spatial variation of the corresponding degrees of freedom in the Ritz solution method. Direct time integrations of the nonlinear equations of motion representing the full aeroelastic system are conducted using the well-known Runge–Kutta method. A comprehensive insight is provided over the effect of parameters such as the lamination fiber angle and the sweep angle on the stability margins and the limit cycle oscillation behavior of the system. Integration of the interpolation method employed for the evaluation of compressible indicial functions at any Mach number in the subsonic compressible range to the derivation process of the third order nonlinear aeroelastic system of equations based on TWB theory is done for the first time. Results show that flutter speeds obtained by the incompressible unsteady aerodynamics are not conservative and as the backward sweep angle of the wing is increased, post-flutter aeroelastic response of the wing becomes more well-behaved.     

跨音速极限环型颤振的高效数值分析方法

事先建立一个低阶的非线性、非定常气动力模型是开展非线性流场中气动弹性问题研究的一个捷径. 基于CFD方法, 通过计算结构在流场中自激振动的响应来获得系统的训练数据. 采用带输出反馈的循环RBF神经网络, 建立时域非线性气动力降阶模型.耦合结构运动方程和非线性气动力降阶模型, 采用杂交的线性多步方法计算结构在不同速度(动压)下的响应历程, 从而获得模型极限环随速度(动压)变化的特性. 两个典型的跨音速极限环型颤振算例表明, 基于气动力降阶模型方法的计算结果与直接CFD仿真结果吻合很好, 与后者相比其将计算效率提高了1~2个数量级.      

An analytical and experimental investigation into limit-cycle oscillations of an aeroelastic system

We perform an analytical and experimental investigation into the dynamics of an aeroelastic system consisting of a plunging and pitching rigid airfoil supported by a linear spring in the plunge degree of freedom and a nonlinear spring in the pitch degree of freedom. The experimental results show that the onset of flutter takes place at a speed smaller than the one predicted by a quasi-steady aerodynamic approximation. On the other hand, the unsteady representation of the aerodynamic loads accurately predicts the experimental value. The linear analysis details the difference in both formulation and provides an explanation for this difference. Nonlinear analysis is then performed to identify the nonlinear coefficients of the pitch spring. The normal form of the Hopf bifurcation is then derived to characterize the type of instability. It is demonstrated that the instability of the considered aeroelastic system is supercritical as observed in the experiments.     

Dimensionless modeling and nonlinear analysis of a coupled pitch–plunge–leadlag airfoil-based piezoaeroelastic energy harvesting system

In this paper, an airfoil-based piezoaeroelastic energy harvesting system is proposed with an additional supporting device to harvest the mechanical energy from the leadlag motion. A dimensionless dynamic model is built considering the large-effective-angle-of-attack vibrations causing (1) the nonlinear coupling between the pitch–plunge–leadlag motions, (2) the inertia nonlinearity, and (3) the aerodynamic nonlinearity modeled by the ONERA dynamic stall model. Cubic hardening stiffness is introduced in the pitch degree of freedom for persistent vibrations with acceptable amplitude beyond the flutter boundary. The nonlinear aeroelastic response and the average power output are numerically studied. Limit cycle oscillations are observed and, as the flow velocity exceeds a secondary critical speed, the system experiences complex vibrations. The power output from the leadlag motion is smaller than that from the plunge motion, whereas the gap is narrowed with increasing flow velocity. Case studies are performed toward the effects of several dimensionless system parameters, including the nonlinear torsional stiffness, airfoil mass eccentricity, airfoil radius of gyration, mass of the supporting devices, and load resistances in the external circuits. The optimal parameter values for the power outputs from the plunge and leadlag motions are, respectively, obtained. Beyond the secondary critical speed, it is shown that the variations of the power outputs with those parameters become irregular with fluctuations and multiple local maximums. The bifurcation analysis shows the mutual transitions between the limit cycle oscillations, multi-periodic vibrations, and possible chaos. The influences of these complex vibrations on the power outputs are discussed.     

Assessment of added mass effects on flutter boundaries using the Leishman–Beddoes dynamic stall model

We consider the dynamics of a typical airfoil section both in forced and free oscillations and investigate the importance of the added mass terms, i.e. the second derivatives in time of the pitch angle and plunge displacement. The structural behaviour is modelled by linear springs in pitch and plunge and the aerodynamic loading represented by our interpretation of the state-space version of the Leishman–Beddoes semi-empirical model. The added mass terms are often neglected since this leads to an explicit system of ODEs amenable for solution using standard ODE solvers. We analyse the effect of neglecting the added mass terms in forced oscillations about a set of mean angles of incidence by comparing the solutions obtained with the explicit and implicit systems of ODEs and conclude that their differences amount to a time lag that increases at a constant rate with increases of the reduced frequency. To determine the effect of the added mass terms in free oscillations, we introduce a spring offset angle to obtain static equilibrium positions at various degrees of incidence. We analyse the stability of the explicit and implicit aeroelastic systems about those positions and compare the locations of the respective flutter points calculated as Hopf bifurcation points. For low values of the spring offset angle, added mass effects are significant for low values of the mass ratio, or the ratio of natural frequencies, of the aeroelastic system. For high values of the spring offset angle, corresponding to stall flutter, we observe that their effect is greater for large values of the mass ratio.     

CHARACTERIZATION OF THE BEHAVIOUR OF A SIMPLE AEROSERVOELASTIC SYSTEM WITH CONTROL NONLINEARITIES

The characterization of the behaviour of nonlinear aeroelastic systems has become a very important research topic in the Aerospace Industry. However, most work carried to-date has concentrated upon systems containing structural or aerodynamic nonlinearities. The purpose of this paper is to study the stability of a simple aeroservoelastic system with nonlinearities in the control system and power control unit. The work considers both structural and control law nonlinearities and assesses the stability of the system response using bifurcation diagrams. It is shown that simple feedback systems designed to increase the stability of the linearized system also stabilize the nonlinear system, although their effects can be less pronounced. Additionally, a nonlinear control law designed to limit the control surface pitch response was found to increase the flutter speed considerably by forcing the system to undergo limit cycle oscillations instead of fluttering. Finally, friction was found to affect the damping of the system but not its stability, as long as the amplitude of the frictional force is low enough not to cause stoppages in the motion.     

Analysis of a piecewise linear aeroelastic system with and without tuned vibration absorber

Nonlinear Dynamics - The dynamics of a two-degrees-of-freedom (pitch–plunge) aeroelastic system is investigated. The aerodynamic force is modeled as a piecewise linear function of the...     

Nonlinear Analysis of Store-Induced Limit Cycle Oscillations

Flight tests of modern high-performance fighter aircraft reveal the presence of limit cycle oscillation (LCO) responses for aircraft with certain external store configurations. Conventional linear aeroelastic analysis predicts flutter for conditions well beyond the operational envelope, yet these store-induced LCO responses occur at flight conditions within the flight envelope. Several nonlinear sources may be present, including aerodynamic effects such as flow separation and shock-boundary layer interaction and structural effects such as stiffening, damping, and system kinematics. No complete theory has been forwarded to accurately explain the mechanisms responsible. This research examines a two degree-of-freedom aeroelastic system which possesses kinematic nonlinearities and a strong nonlinearity in pitch stiffness. Nonlinear analysis techniques are used to gain insight into the characteristics of the behavior of the system. Numerical simulation is used to verify and validate the analysis. It is found that when system damping is low, the system clearly exhibits nonlinear interaction between aeroelastic modes. It is also shown that although certain applied forcing conditions may appear negligible, these same forces produce large amplitude LCOs under specific realizable circumstances.     

Nonlinear aeroelastic analysis of a two-dimensional wing with control surface in supersonic flow

Based on the piston theory of supersonic flow and the energy method, a two dimensional wing with a control surface in supersonic flow is theoretically modeled, in which the cubic stiffness in the torsional direction of the control surface is considered. An approximate method of the cha- otic response analysis of the nonlinear aeroelastic system is studied, the main idea of which is that under the condi- tion of stable limit cycle flutter of the aeroelastic system, the vibrations in the plunging and pitching of the wing can approximately be considered to be simple harmonic excita- tion to the control surface. The motion of the control surface can approximately be modeled by a nonlinear oscillation of one-degree-of-freedom. The range of the chaotic response of the aeroelastic system is approximately determined by means of the chaotic response of the nonlinear oscillator. The rich dynamic behaviors of the control surface are represented as bifurcation diagrams, phase-plane portraits and PS diagrams. The theoretical analysis is verified by the numerical results.     

Frequency lock-in during nonlinear vibration of an airfoil coupled with van der Pol Oscillator

The frequency lock-in during the nonlinear vibration of a turbomachinery blade is modeled using a spring-mounted airfoil coupled with a van der Pol Oscillator (VDP) oscillator. The proposed reduced-order model uses the nonlinear VDP oscillator to represent the oscillatory nature of wake dynamics caused by the vortex shedding. The damping term in the VDP oscillator is assumed to be nonlinear. The coupled equations governing the pitch and plunge motion of an airfoil are used to approximate the vibration of a turbomachinery blade. Springs having cubic-order nonlinearity for their stiffnesses are used to mount the airfoil. The unsteady lift acting on the blade is modeled using a self-excited nonlinear wake oscillator. The model for wake dynamics takes into account the influence of blade inertia. The nonlinear coupled three degrees of freedom (dof) aeroelastic system is studied for instability resulting in the frequency lock-in phenomenon. The equations are transformed into non-dimensional form, and then the frequencies of the coupled system are plotted to demonstrate the frequency lock-in. Further, the method of multiple scales is used to derive modulation equations which represent the amplitude and phase of the oscillation. The results obtained using the method of multiple scales are compared with direct numerical solutions to verify the present modeling method. The steady-state amplitudes of the response are plotted against the detuning parameter, which represents the frequency response curve. Further, the sensitivity of non-dimensional parameters such as coupling coefficients, mass ratio, reduced velocity, static unbalance, structural damping coefficient and the ratio of uncoupled pitch and plunge natural frequencies on the frequency response is investigated. The study revealed that parameters such as mass ratio, reduced velocity, structural damping coefficient, and coupling coefficients have a stronger influence in suppressing the amplitude of vibration. Meanwhile, parameters such as the frequency ratio, static unbalance, reduced velocity, and mass ratio significantly affect the range of frequency in which the lock-in phenomenon happens. Further, linear perturbation analysis is done to understand the qualitative effect of the system parameters such as coupling coefficients, mass ratio, frequency ratio, and static unbalance on the range of lock-in.     

Response analysis of a pitch–plunge airfoil with structural and aerodynamic nonlinearities subjected to randomly fluctuating flows

This study focuses on numerically investigating the response dynamics of a pitch–plunge airfoil with structural nonlinearity under dynamic stall conditions. The aeroelastic responses are investigated for both deterministic and randomly time varying flow conditions. To that end, a pitch–plunge airfoil under dynamic stall condition is considered and the nonlinear aerodynamic loads are computed using a Leishman–Beddoes formulation. It is shown that the presence of structural nonlinearities can give rise to a variety of dynamical responses in the pre-flutter regime. Next, a response analysis under the presence of a randomly fluctuating wind is carried out. It is demonstrated that the route to flutter occurs via a regime of pre-flutter oscillations called intermittency. Finally, the manifestation of these stochastic responses is characterized by invoking stochastic bifurcation concepts. The route to flutter via intermittency is presented in terms of topological changes occurring in the joint-probability density function of the state variables.     

Aeroelastic inverse: Estimation of aerodynamic loads during large amplitude limit cycle oscillations

This paper presents an algorithm to compute the aerodynamic forces and moments of an aeroelastic wing undergoing large amplitude heave and pitch limit cycle oscillations. The technique is based on inverting the equations of motion to solve for the lift and moment experienced by the wing. Bayesian inferencing is used to estimate the structural parameters of the system and generate credible intervals on the lift and moment calculations. The inversion technique is applied to study the affect of mass coupling on limit cycle oscillation amplitude. Examining the force, power, and energy of the system, the reasons for amplitude growth with wind speed can be determined. The results demonstrate that the influence of mass coupling on the pitch–heave difference is the driving factor in amplitude variation. The pitch–heave phase difference not only controls how much aerodynamic energy is transferred into the system but also how the aerodynamic energy is distributed between the degrees of freedom.     

An analytical control law of length rate for tethered satellite system

Based on the nonlinear dynamic equations of a tethered satellite system with three-dimensional attitude motion, an analytical tether length rate control law for deployment is derived from the equilibrium positions of the system and the scheme of the value range of the expected in-plane pitch angle. The proposed control law can guarantee that the tensional force acting on the end of the tether remains positive. The oscillation of the out-of-plane roll motion in conjunction with the in-plane pitch motion is effectively suppressed during deployment control. The analytical control law is still applicable, even if the system runs on a Keplerian elliptical orbit with a large eccentricity. The local stability of the non-autonomous system during deployment control is analyzed using the Floquet theory, and the global behavior is numerically verified using simple cell mapping. The numerical simulations in the paper demonstrate the proposed analytical control law.