Influence of friction and asymmetric freeplay on the limit cycle oscillation in aeroelastic system: An extended Hénon’s technique to temporal integration

Nonlinear effects such as friction and freeplay on the control surfaces can affect aeroelastic dynamics during flight. In particular, these nonlinearities can induce limit cycle oscillations (LCO), changing the system stability, and because of this it is essential to employ computational methods to predict this type of motion during the aircraft development cycle. In this context, the present article presents a matrix notation for describing the Hénon’s method used to reduce errors when considering piecewise linear nonlinearities in the numerical integration process. In addition, a new coordinate system is used to write the aeroelastic system of equations. The proposal defines a displacement vector with generalized and physical variables to simplify the computational implementation of the Hénon’s technique. Additionally, the article discusses the influence of asymmetric freeplay and friction on the LCO of an airfoil with control surface. The results show that the extended Hénon’s technique provides more accurate LCO predictions, that friction can change the frequency and amplitude of these motions, and the asymmetry of freeplay is important to determine the LCO behavior.     

A theoretical and experimental investigation of the effects of a steady angle of attack on the nonlinear flutter of a delta wing plate model

Limit cycle oscillations (LCO) of wings on certain modern high performance aircraft have been observed in flight and in wind tunnel experiments. Whether the physical mechanism that gives rise to this behavior is a fluid or structural nonlinearity or both is still uncertain. It has been shown that an aeroelastic theoretical model with only a structural nonlinearity can predict accurately the limit cycle behavior at low subsonic flow for a plate-like wing at zero angle of attack. Changes in the limit cycle and flutter behavior as the angle of attack is varied have also been observed in flight. It has been suggested that this sensitivity to angle of attack is due to a fluid nonlinearity. In this investigation, we study the flutter and limit cycle behavior of a wing in low subsonic flow at small steady angles of attack. Experimental results are compared to those predicted using an aeroelastic theoretical model with only a structural nonlinearity. Results from both experiment and theory show a change in flutter speed as the steady angle of attack is varied. Also the LCO magnitude increased at a given velocity as the angle of attack was increased for both the experiment and theory. While not proving that the observed sensitivity to angle of attack of LCO in aircraft is due to a structural nonlinearity, the results do show that a change in the aeroelastic behavior at angles of attack can be caused by a structural nonlinearity as well as a fluid nonlinearity. In this paper, only structural nonlinearities are considered, but an extension to include aerodynamic nonlinearities would be very worthwhile.     

Identification of Nonlinear Aeroelastic Systems Based on the Volterra Theory: Progress and Opportunities

The identification of nonlinear aeroelastic systems based on the Volterra theory of nonlinear systems is presented. Recent applications of the theory to problems in computational and experimental aeroelasticity are reviewed. Computational results include the development of computationally efficient reduced-order models (ROMs) using an Euler/Navier–Stokes flow solver and the analytical derivation of Volterra kernels for a nonlinear aeroelastic system. Experimental results include the identification of aerodynamic impulse responses, the application of higher-order spectra (HOS) to wind-tunnel flutter data, and the identification of nonlinear aeroelastic phenomena from flight flutter test data of the active aeroelastic wing (AAW) aircraft.     

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

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

Limit cycle oscillation suppression of 2-DOF airfoil using nonlinear energy sink

This paper presents a novel mechanical attachment, i.e., nonlinear energy sink （NES）, for suppressing the limit cycle oscillation （LCO） of an airfoil. The dynamic responses of a two-degree-of-freedom （2-DOF） airfoil coupled with an NES are studied with the harmonic balance method. Different structure parameters of the NES, i.e., mass ratio between the NES and airfoil, NES offset, NES damping, and nonlinear stiffness in the NES, are chosen for studying the effect of the LCO suppression on an aeroelastic system with a supercritical Hopf bifurcation or subcritical Hopf bifurcation, respectively. The results show that the structural parameters of the NES have different influence on the supercritical Hopf bifurcation system and the subcritical Hopf bifurcation system.     

Effects of damping on the speed of increase and amplitude of limit cycle for an aircraft braking system subjected to mode-coupling instability

A nonlinear model of an aircraft braking system is presented and used to investigate the effects of damping on the stability in Chevillot et al. (Arch Appl Mech 78(12):949–963, 2008). It has been shown that the addition of damping into the equations of motion does not lead systematically to the stabilization of the system. In the case of a mode-coupling instability, there is indeed an optimal ratio between the modal damping coefficients of the two modes in coalescence, that maximize the stable area. But the stable area is not a sufficient criterion. In dynamics, the amplitude of the vibrations and the transient behavior characterized by the speed of increase of the oscillations are best indicators. In this paper, the same nonlinear model of the aircraft braking system is used to compute time-history responses by integration of the full set of the nonlinear dynamic equations. The aim of the study is to evaluate the effects of damping on the nonlinear dynamics of the brake. It is shown that damping may be very efficient to significantly reduce and slow down the increase of the friction-induced vibrations. But, in the same way as for the stability area, there exists a value of the damping ratio that optimizes the effects of damping.     

Limit cycle oscillations of rectangular cantilever wings containing cubic nonlinearity in an incompressible flow

Limit cycle oscillations (LCO) as well as nonlinear aeroelastic analysis of rectangular cantilever wings with a cubic nonlinearity are investigated. Aeroelastic equations of a rectangular cantilever wing with two degrees of freedom in an incompressible potential flow are presented in the time domain. The harmonic balance method is modified to calculate the LCO frequency and amplitude for rectangular wings. In order to verify the derived formulation, flutter boundaries are obtained via a linear analysis of the derived system of equations for five different cases and compared with experimental data. Satisfactory results are gained through this comparison. The problem of finding the LCO frequency and amplitude is solved via applying the two methods discussed for two different cases with hardening cubic nonlinearities. The results from first-, third- and fifth-order harmonic balance methods are compared with the results of an exact numerical solution. A close agreement is obtained between these harmonic balance methods and the exact numerical solution of the governing aeroelastic equations. Finally, the nonlinear aeroelastic analysis of a rectangular cantilever wing with a softening nonlinearity is studied.     

An Investigation of Internal Resonance in Aeroelastic Systems

Although the study of internal resonance in mechanical systems has been given significant consideration, minimal attention has been given to internal resonance for systems which consider the presence of aerodynamic forces. Herein, the investigators examine the possible existence of internal resonances, and the related nonlinear pathologies that such responses may have, for an aeroelastic system which possesses nonlinear aerodynamic loads. Evidence of internal resonance is presented for specific classes of aeroelastic systems, and such adverse response indicates nonlinearities may lead to aeroelastic instabilities that are not predicted by traditional (linear) approaches.     

Coupled framework for limit-cycle oscillations modeling based on leading-edge vortex shedding

Current trends in the aircraft industry involve higher aspect-ratio wings made of lighter materials. These trends seek to reduce fuel emissions and increase flight efficiency by reducing drag to lift ratio and overall weight, respectively, of the aircraft. This results in reduced structural stiffness and coupling between the aeroelastic modes and flight dynamics. The flutter phenomenon is of particular interest for aeroelastic studies, and modeling post-flutter limit-cycle oscillations (LCO) is a challenging problem.Several studies have been developed to allow fast simulations of the highly non-linear aerodynamic situations, with leading-edge vortex modulation been a proved solution for modeling some forms of LCOs in airfoils. This article proposes a framework based on the 3D expansion of this method using strip theory and coupling with modal structural model for simulations of aerodynamic based non-linear phenomenon. A cantilevered flat plate is used for testing and validating the framework against wind-tunnel experiments and the industry standard approach. The results show that the proposed model is able to capture the main behavior of the LCO observed in the experiments and is directly comparable with the current approaches used at the industry. The framework allows for scalability and is also fast enough to provide time-based results in under two days for a desktop simulation, reducing the need of expensive cluster computations. Finally, since it is completely physics-based it allows for the engineer to get insights on the aerodynamic flow at a fraction of the cost of more detailed CFD models.     

LAUNCH VEHICLE BUFFETING WITH AEROELASTIC COUPLING EFFECTS

Launch vehicle structural responses can couple with transonic flow state transitions at the nose of payload fairings. This self-sustained coupling yields a nonlinear equation of motion that can be analyzed using the force–response relationship and the periodicity condition. The traditional analysis approach for this phenomenon, however, linearizes the equation of motion by converting the alternating flow forces into an aerodynamic damping term and defines a stability criterion as the response amplitude that yields zero net system damping. This work clarifies the relationship between the present and traditional methods, and compares results and conclusions. The feasibility of modifying a launch vehicle response analysis of buffeting (random pressure fluctuations caused by turbulent flow) to include aeroelastic coupling effects is also explored. The aerodynamic stiffness and damping terms formulated herein are consistent with trends observed in wind-tunnel test data. It is shown, however, that the modified buffet analysis can be inaccurate, particularly when the aeroelastic coupling contribution does not dominate the system response.     

Aeroelastic stability analysis considering a continuous flight envelope

This paper presents a new methodology to analyze aeroelastic stability in a continuous range of flight envelope with varying parameter of velocity and altitude. The focus of the paper is to demonstrate that linear matrix inequalities can be used to evaluate the aeroelastic stability in a region of flight envelope instead of a single point, like classical methods. The proposed methodology can also be used to study if a system remains stable during an arbitrary motion from one point to another in the flight envelope, i.e., when the problem becomes time-variant. The main idea is to represent the system as a polytopic differential inclusion system using rational function approximation to write the model in time domain. The theory is outlined and simulations are carried out on the benchmark AGARD 445.6 wing to demonstrate the method. The classical pk-method is used for comparing results and validating the approach. It is shown that this method is efficient to identify stability regions in the flight envelope.     

Stability analysis of whirl flutter in rotor-nacelle systems with freeplay nonlinearity

Nonlinear Dynamics - Tiltrotor aircraft are growing in prevalence due to the usefulness of their unique flight envelope. However, aeroelastic stability—particularly whirl flutter...     

飞行器非线性气动伺服弹性力学

现代飞行器日益呈现结构轻质化、控制系统宽通带和高权限的发展趋势. 因此, 非定常气动力、柔性结构和主动控制系统三者间的耦合力学成为重要的研究领域. 自20世纪80年代起, 航空界开始关注受控飞行器的气动弹性稳定性以及主动控制问题, 但对气动/结构的非线性效应、控制回路时滞对受控飞行器动力学行为的影响规律研究尚不充分. 研究这些影响规律不仅涉及非线性、高维数、多变参数和时滞效应等难题, 而且必须面对空气动力、飞行器结构、驱动机构、控制系统之间的强耦合问题. 其中的前沿难题是: 发展非线性气动伺服弹性动力学建模理论, 揭示上述因素诱发受控气动弹性振动的动力学机理, 开展气动伺服弹性控制风洞实验. 本文针对非线性气动伺服弹性力学所涉及的非线性非定常气动力建模、非线性结构动力学、气动伺服弹性控制律设计、气动伺服弹性实验, 总结相关研究现状和最新进展, 特别是近年来作者学术团队的研究成果, 并对进一步研究给出若干建议.      

Nonlinear airship aeroelasticity

The aeroelastic derivatives for today's aircraft are calculated in the concept phase using a standard procedure. This scheme has to be extended for large airships, due to various nonlinearities in structural and aerodynamic behaviour. In general, the structural model of an airship is physically as well as geometrically nonlinear. The main sources of nonlinearity are large deformations and the nonlinear material behaviour of membranes. The aerodynamic solution is also included in the nonlinear problem, because the deformed airship influences the surrounding flow. Due to these nonlinearities, the aeroelastic problem for airships can only be solved by an iterative procedure. As one possibility, the coupled aerodynamic and structural dynamic problem was handled using linked standard solvers. On the structural side, the Finite-Element program package ABAQUS was extended with an interface to the aerodynamic solver VSAERO. VSAERO is based on the aerodynamic panel method using potential flow theory. The equilibrium of the internal structural and the external aerodynamic forces leads to the structural response and a trimmed flight state for the specified flight conditions (e.g. speed, altitude). The application of small perturbations around a trimmed state produces reaction forces and moments. These constraint forces are then transferred into translational and rotational acceleration fields by performing an inertia relief analysis of the disturbed structural model. The change between the trimmed flight state and the disturbed one yields the respective aeroelastic derivatives. By including the calculated derivatives in the linearised equation of motion system, it is possible to judge the stability and controllability of the investigated airship.     

Studies of Store-Induced Limit-Cycle Oscillations Using a Model with Full System Nonlinearities

The authors investigate limit-cycle oscillations of a wing/store configuration. Unlike typical aeroelastic studies that are based upon a linearized form of the governing equations, herein full system nonlinearities are retained, and include transonic flow effects, coupled responses from the structure, and store-related kinematics and dynamics. Unsteady aerodynamic loads are modeled with the equations from transonic small disturbance theory. The structural dynamics for the cantilevered wing are modeled by the nonlinear equations of motion for a beam. The effects of general store-placement are modeled by the nonlinear equations of motion related to the position-induced nonlinear kinematics. Chordwise deformations of the wing surface, as well as pylon and store flexibility, are assumed negligible. Nonlinear responses are studied by examining bifurcation and related response characteristics using direct simulation. Particular attention is given to cases for which large-time, time-dependent behavior is dependent on initial conditions, as observed for some configurations in flight test. Comparisons of results in which selective nonlinearities are excluded indicate that the accurate prediction of nonlinear responses such as limit cycle oscillations (LCOs) may depend upon consideration of all nonlinearities related to the full system.     

An efficient method for nonlinear aeroelasticy of slender wings

This paper aims the nonlinear aeroelastic analysis of slender wings using a nonlinear structural model coupled with the linear unsteady aerodynamic model. High aspect ratio and flexibility are the specific characteristic of this type of wings. Wing flexibility, coupled with long wingspan can lead to large deflections during normal flight operation of an aircraft; therefore, a wing in vertical/forward-afterward/torsional motion using a third-order form of nonlinear general flexible Euler–Bernoulli beam equations is used for structural modeling. Unsteady linear aerodynamic strip theory based on the Wagner function is used for determination of aerodynamic loading on the wing. Combining these two types of formulation yields nonlinear integro-differentials aeroelastic equations. Using the Galerkin’s method and a mode summation technique, the governing equations will be solved by introducing a numerical method without the need to adding any aerodynamic state space variables and the corresponding equations related to these variables of the problem. The obtained equations are solved to predict the aeroelastic response of the problem. The obtained results for a test case are compared with those of some other works and show a good agreement between results.     

Assessment of low-altitude atmospheric turbulence models for aircraft aeroelasticity

We investigate the dynamic aeroelastic response of large but slow aircraft in low-altitude atmospheric turbulence. To this end, three turbulence models of increasing fidelity, namely, the one-dimensional von Kármán model, the two-dimensional Kaimal model and full three-dimensional wind fields extracted from large-eddy simulations (LES) are used to simulate ambient turbulence near the ground. Load calculations and flight trajectory predictions are conducted for a representative high-aspect-ratio wing aircraft, using a fully coupled nonlinear flight dynamics/aeroelastic model, when it operates in background atmospheric turbulence generated by the aforementioned models. Comparison of load envelopes and spectral content, on vehicles of varying flexibility, shows strong dependency between the selected turbulence model and aircraft aeroelastic response (e.g. 58% difference in the predicted magnitude of the wing root bending moment between LES and von Kármán models). This is mainly due to the presence of large flow structures at low altitudes that have comparable dimensions to the vehicle, and which despite the relatively small wind speeds within the Earth boundary layer, result in overall high load events for slow-moving vehicles. Results show that one-dimensional models that do not capture those effects provide fairly non-conservative load estimates and are unsuitable for very flexible airframe design.     

A time–frequency technique for the stability analysis of impulse responses from nonlinear aeroelastic systems

A time–frequency method is proposed for the analysis of response time histories from nonlinear aeroelastic systems. The approach is based on a time-varying curve-fit of the short time Fourier transform of the impulse response. It is shown that the method can be used in order to obtain a clear picture of the sub-critical stability of a number of aeroelastic systems with a variety of structural and aerodynamic nonlinearities. Additionally, frequency and amplitude information can be obtained for both the linear and nonlinear signatures of the response signals in the sub- and post-critical regions. Finally, it is shown that, given certain types of nonlinear functions, sub-critical damping trends can be extrapolated to predict bifurcation airspeeds.     

Response Statistics of Three-Degree-of-Freedom Nonlinear System to Narrow-Band Random Excitation

The principal resonance of a 3-DOF nonlinear system to narrow-band random external excitations is investigated. The method of multiple scales is used to derive the equations for modulation of amplitude and phase. The behavior, stability and bifurcation of steady-state responses are studied by means of qualitative analysis. The effects of damping, detuning, and excitation intensity on responses are analyzed. The theoretical analyses are verified by numerical results. Both theoretical analyses and numerical simulations show that when the intensity of the random excitation increases, the nontrivial steady state solution may change from a limit cycle to a diffused limit cycle. Under some conditions, co-existence of two kinds of stable steady-state solutions, saturation and jump phenomena may occur. The stationary probability density function of responses for the co-existence case is obtained approximately.     

Aircraft T-tail flutter predictions using computational fluid dynamics

The paper presents the application of computational aeroelasticity (CA) methods to the analysis of a T-tail stability in transonic regime. For this flow condition unsteady aerodynamics show a significant dependency from the aircraft equilibrium flight configuration, which rules both the position of shock waves in the flow field and the load distribution on the horizontal tail plane. Both these elements have an influence on the aerodynamic forces, and so on the aeroelastic stability of the system. The numerical procedure proposed allows to investigate flutter stability for a free-flying aircraft, iterating until convergence the following sequence of sub-problems: search for the trimmed condition for the deformable aircraft; linearize the system about the stated equilibrium point; predict the aeroelastic stability boundaries using the inferred linear model. An innovative approach based on sliding meshes allows to represent the changes of the computational fluid domain due to the motion of control surfaces used to trim the aircraft. To highlight the importance of keeping the linear model always aligned to the trim condition, and at the same time the capabilities of the computational fluid dynamics approach, the method is applied to a real aircraft with a T-tail configuration: the P180.