Investigation of Limit Cycle Oscillations Using Aeroelastic-Harmonic Balance Method

This work investigates limit cycle oscillations in the transonic regime. A novel approach to predict Limit Cycle Oscillations using high fidelity analysis is exploited to accelerate calculations. The method used is an Aeroeasltic Harmonic Balance approach, which has been proven to be efficient and able to predict periodic phenomena. The behaviour of limit cycle oscillations is analysed using uncertainty quantification tools based on polynomial chaos expansions. To improve the efficiency of the sampling process for the polynomial-chaos expansions an adaptive sampling procedure is used. These methods are exercised using two problems: a pitch/plunge aerofoil and a delta-wing. Results indicate that Mach n. variability is determinant to the amplitude of the LCO for the 2D test case, whereas for the wing case analysed here, variability in the Mach n. has an almost negligible influence in amplitude variation and the LCO frequency variability has an almost linear relation with Mach number. Further test cases are required to understand the generality of these results.

Prediction of limit cycle oscillations under uncertainty using a Harmonic Balance method

The Harmonic Balance method is an attractive solution for computing periodic responses and can be an alternative to time domain methods, at a reduced computational cost. The current paper investigates using a Harmonic Balance method for simulating limit cycle oscillations under uncertainty. The Harmonic Balance method is used in conjunction with a non-intrusive polynomial-chaos approach to propagate variability and is validated against Monte Carlo analysis. Results show the potential of the approach for a range of nonlinear dynamical systems, including a full wing configuration exhibiting supercritical and subcritical bifurcations, at a fraction of the cost of performing time domain simulations.

Prediction of Transonic Limit-Cycle Oscillations Using an Aeroelastic Harmonic Balance Method

This work proposes a novel approach to compute transonic limit-cycle oscillations using high-fidelity analysis. Computational-Fluid-Dynamics based harmonic balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling harmonic balance methods to accurately capture limit-cycle oscillations; this is achieved by defining a frequency-updating procedure based on a coupled computational-fluid-dynamics/computational-structural-dynamics harmonic balance formulation to find the limit-cycle oscillation condition. A pitch/plunge airfoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both limit-cycle oscillation amplitude and frequency while producing at least a one-order-of-magnitude reduction in computational time.

Reducing Parametric Uncertainty in Limit Cycle Oscillations

The assimilation of discrete higher fidelity data points with model predictions can be used to achieve a reduction in the uncertainty of the model input parameters which generate accurate predictions. The problem investigated here involves the prediction of limit-cycle oscillations using a High-Dimensional Harmonic Balance method (HDHB). The efficiency of the HDHB method is exploited to enable calibration of structural input parameters using a Bayesian inference technique. Markov-chain Monte Carlo is employed to sample the posterior distributions. Parameter estimation is carried out on both a pitch/plunge aerofoil and Goland wing configuration. In both cases significant refinement was achieved in the distribution of possible structural parameters allowing better predictions of their
true deterministic values.

Die aeroelastische Stabilitätsanalyse – ein praxisnaher Ansatz zur intervalltheoretischen Betrachtung von Modellierungsunsicherheiten am Flugzeug

Die Untersuchung des dynamischen aeroelastischen Stabilitätsverhaltens von Flugzeugen erfordert sehr komplexe Rechenmodelle, welche die wesentlichen elastomechanischen und instationären aerodynamischen Eigenschaften der Konstruktion wiedergeben sollen. Bei der Modellbildung müssen einerseits Vereinfachungen und Idealisierungen im Rahmen der Anwendung der Finite Elemente Methode und der aerodynamischen Theorie vorgenommen werden, deren Auswirkungen auf das Simulationsergebnis zu bewerten sind. Andererseits können die strukturdynamischen Kenngrößen durch den Standschwingungsversuch identifiziert werden, wobei die Ergebnisse Messungenauigkeiten enthalten. Für eine robuste Flatteruntersuchung müssen die identifizierten Unwägbarkeiten in allen Prozessschritten über die Festlegung von unteren und oberen Schranken konservativ ermittelt werden, um für alle Flugzustände eine ausreichende Flatterstabilität sicherzustellen. Zu diesem Zweck wird in der vorliegenden Arbeit ein Rechenverfahren entwickelt, welches die klassische Flatteranalyse mit den Methoden der Fuzzy- und Intervallarithmetik verbindet. Dabei werden die Flatterbewegungsgleichungen als parameterabhängiges nichtlineares Eigenwertproblem formuliert. Die Änderung der komplexen Eigenlösung infolge eines veränderlichen Einflussparameters wird mit der Methode der numerischen Fortsetzung ausgehend von der nominalen Startlösung verfolgt. Ein modifizierter Newton-Iterations-Algorithmus kommt zur Anwendung. Als Ergebnis liegen die berechneten aeroelastischen Dämpfungs- und Frequenzverläufe in Abhängigkeit von der Fluggeschwindigkeit mit Unschärfebändern vor.

The use of Gramian matrices for aeroelastic stability analysis

Most of the established procedures for analysis of aeroelastic flutter in the development of aircraft are based on frequency domain methods. Proposing new methodologies in this field is always a challenge, because the new methods need to be validated by many experimental procedures. With the interest for new flight control systems and nonlinear behavior of aeroelastic structures, other strategies may be necessary to complete the analysis of such systems. If the aeroelastic model can be written in time domain, using state-space formulation, for instance, then many of the tools used in stability analysis of dynamic systems may be used to help providing an insight into the aeroelastic phenomenon. In this respect, this paper presents a discussion on the use of Gramian matrices to determine conditions of aeroelastic flutter. The main goal of this work is to introduce how observability gramian matrix can be used to identify the system instability. To explain the approach, the theory is outlined and simulations are carried out on two benchmark problems. Results are compared with classical methods to validate the approach and a reduction of computational time is obtained for the second example. © 2013 Douglas Domingues Bueno et al.

Influence of structural nonlinearities in stall-induced aeroelastic response of pitching airfoils

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Characterization of grazing bifurcation in airfoils with control surface freeplay nonlinearity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Representation and analysis of control surface freeplay nonlinearity

Different representations for a control surface freeplay nonlinearity in a three degree of freedom aeroelastic system are assessed. These are the discontinuous, polynomial and hyperbolic tangent representations. The Duhamel formulation is used to model the aerodynamic loads. Assessment of the validity of these representations is performed through comparison with previous experimental observations. The results show that the instability and nonlinear response characteristics are accurately predicted when using the discontinuous and hyperbolic tangent representations. On the other hand, the polynomial representation fails to predict chaotic motions observed in the experiments. (c) 2012 Elsevier Ltd. All rights reserved.

Modeling and identification of freeplay nonlinearity

A nonlinear analysis is performed for the purpose of identification of the pitch freeplay nonlinearity and its effect on the type of bifurcation of a two degree-of-freedom aeroelastic system. The databases for the identification are generated from experimental investigations of a pitch-plunge rigid airfoil supported by a nonlinear torsional spring. Experimental data and linear analysis are performed to validate the parameters of the linearized equations. Based on the periodic responses of the experimental data which included the flutter frequency and its third harmonics, the freeplay nonlinearity is approximated by a polynomial expansion up to the third order. This representation allows us to use the normal form of the Hopf bifurcation to characterize the type of instability. Based on numerical integrations, the coefficients of the polynomial expansion representing the freeplay nonlinearity are identified. Published by Elsevier Ltd.

Bifurcation analysis of an aeroelastic system with concentrated nonlinearities

Analytical and numerical analyses of the nonlinear response of a three-degree-of-freedom nonlinear aeroelastic system are performed. Particularly, the effects of concentrated structural nonlinearities on the different motions are determined. The concentrated nonlinearities are introduced in the pitch, plunge, and flap springs by adding cubic stiffness in each of them. Quasi-steady approximation and the Duhamel formulation are used to model the aerodynamic loads. Using the quasi-steady approach, we derive the normal form of the Hopf bifurcation associated with the system's instability. Using the nonlinear form, three configurations including supercritical and subcritical aeroelastic systems are defined and analyzed numerically. The characteristics of these different configurations in terms of stability and motions are evaluated. The usefulness of the two aerodynamic formulations in the prediction of the different motions beyond the bifurcation is discussed.

Simulazione numerica di fenomeni aeroelastici con applicazione alle strutture civili

Il presente lavoro tratta lo studio dei fenomeni aeroelastici di interazione fra fluido e struttura, con il fine di provare a simularli mediante l’ausilio di un codice agli elementi finiti. Nel primo capitolo sono fornite alcune nozioni di fluidodinamica, in modo da rendere chiari i passaggi teorici fondamentali che portano alle equazioni di Navier-Stokes governanti il moto dei fluidi viscosi. Inoltre è illustrato il fenomeno della formazione di vortici a valle dei corpi tozzi dovuto alla separazione dello strato limite laminare, con descrizione anche di alcuni risultati ottenuti dalle simulazioni numeriche. Nel secondo capitolo vengono presi in rassegna i principali fenomeni di interazione fra fluido e struttura, cercando di metterne in luce le fondamenta della trattazione analitica e le ipotesi sotto le quali tale trattazione è valida. Chiaramente si tratta solo di una panoramica che non entra in merito degli sviluppi della ricerca più recente ma fornisce le basi per affrontare i vari problemi di instabilità strutturale dovuti a un particolare fenomeno di interazione con il vento. Il terzo capitolo contiene una trattazione più approfondita del fenomeno di instabilità per flutter. Tra tutti i fenomeni di instabilità aeroelastica delle strutture il flutter risulta il più temibile, soprattutto per i ponti di grande luce. Per questo si è ritenuto opportuno dedicargli un capitolo, in modo da illustrare i vari procedimenti con cui si riesce a determinare analiticamente la velocità critica di flutter di un impalcato da ponte, a partire dalle funzioni sperimentali denominate derivate di flutter. Al termine del capitolo è illustrato il procedimento con cui si ricavano sperimentalmente le derivate di flutter di un impalcato da ponte. Nel quarto capitolo è presentato l’esempio di studio dell’impalcato del ponte Tsing Ma ad Hong Kong. Sono riportati i risultati analitici dei calcoli della velocità di flutter e di divergenza torsionale dell’impalcato e i risultati delle simulazioni numeriche effettuate per stimare i coefficienti aerodinamici statici e il comportamento dinamico della struttura soggetta all’azione del vento. Considerazioni e commenti sui risultati ottenuti e sui metodi di modellazione numerica adottati completano l’elaborato.

On the identification of bridge decks aeroelastic coefficients: the covariance block hankel matrix method

Il seguente elaborato si concentra sull'identifi�cazione strutturale di sistemi soggetti a sollecitazioni aeroelastiche e nello speci�fico l'attenzione viene rivolta ad impalcati da ponte. Si analizzano i concetti principali caratterizzanti il campo dell'aeroelasticità indagando i fattori dominanti che entrano in gioco sul piano teorico. In seguito, si considera il metodo di identifi�cazione strutturale chiamato Covariance Block Hankel Matrix (CBHM) utilizzato come strumento di derivazione dei coeffi�cienti aeroelastici propri della struttura. Infi�ne, si indaga il comportamento di questo metodo di identi�ficazione al variare di una serie di parametri chiave e all'interno di diversi scenari, visionando risultati ottenuti tramite una serie di test eff�ettuati per provare l'a�dattabilità del metodo stesso al variare delle condizioni che caratterizzano il sistema.

MODEL UPDATING AND STRUCTURAL HEALTH MONITORING OF HORIZONTAL AXIS WIND TURBINES VIA ADVANCED SPINNING FINITE ELEMENTS AND STOCHASTIC SUBSPACE IDENTIFICATION METHODS

Wind energy has been one of the most growing sectors of the nation’s renewable energy portfolio for the past decade, and the same tendency is being projected for the upcoming years given the aggressive governmental policies for the reduction of fossil fuel dependency. Great technological expectation and outstanding commercial penetration has shown the so called Horizontal Axis Wind Turbines (HAWT) technologies. Given its great acceptance, size evolution of wind turbines over time has increased exponentially. However, safety and economical concerns have emerged as a result of the newly design tendencies for massive scale wind turbine structures presenting high slenderness ratios and complex shapes, typically located in remote areas (e.g. offshore wind farms). In this regard, safety operation requires not only having first-hand information regarding actual structural dynamic conditions under aerodynamic action, but also a deep understanding of the environmental factors in which these multibody rotating structures operate. Given the cyclo-stochastic patterns of the wind loading exerting pressure on a HAWT, a probabilistic framework is appropriate to characterize the risk of failure in terms of resistance and serviceability conditions, at any given time. Furthermore, sources of uncertainty such as material imperfections, buffeting and flutter, aeroelastic damping, gyroscopic effects, turbulence, among others, have pleaded for the use of a more sophisticated mathematical framework that could properly handle all these sources of indetermination. The attainable modeling complexity that arises as a result of these characterizations demands a data-driven experimental validation methodology to calibrate and corroborate the model. For this aim, System Identification (SI) techniques offer a spectrum of well-established numerical methods appropriated for stationary, deterministic, and data-driven numerical schemes, capable of predicting actual dynamic states (eigenrealizations) of traditional time-invariant dynamic systems. As a consequence, it is proposed a modified data-driven SI metric based on the so called Subspace Realization Theory, now adapted for stochastic non-stationary and timevarying systems, as is the case of HAWT’s complex aerodynamics. Simultaneously, this investigation explores the characterization of the turbine loading and response envelopes for critical failure modes of the structural components the wind turbine is made of. In the long run, both aerodynamic framework (theoretical model) and system identification (experimental model) will be merged in a numerical engine formulated as a search algorithm for model updating, also known as Adaptive Simulated Annealing (ASA) process. This iterative engine is based on a set of function minimizations computed by a metric called Modal Assurance Criterion (MAC). In summary, the Thesis is composed of four major parts: (1) development of an analytical aerodynamic framework that predicts interacted wind-structure stochastic loads on wind turbine components; (2) development of a novel tapered-swept-corved Spinning Finite Element (SFE) that includes dampedgyroscopic effects and axial-flexural-torsional coupling; (3) a novel data-driven structural health monitoring (SHM) algorithm via stochastic subspace identification methods; and (4) a numerical search (optimization) engine based on ASA and MAC capable of updating the SFE aerodynamic model.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.