On the design of an electromagnetic aeroelastic energy harvester from nonlinear flutter

A new energy harvester by coupling the electromagnetic induction and the pitch vibration of a rigid wing is built up in this paper. It is aimed: (1) to harvest energy from the pitch limit cycle oscillation (LCO) of the wing due to the preloaded free-play nonlinearity; (2) to introduce a theoretical analysis scheme based on the equivalent linearized method into the design of this harvester. With the equivalent linearized method, the domains of the single stable LCO and double stable LCOs are respectively obtained. Combining the analytical and numerical solutions, the single stable LCO along with the stable limit cycle amplitude greater than its corresponding unstable one is recognized as the better mode for harvesting, since the larger limit cycle domain is induced and the more energy are yielded. Based on such chosen mode, analyses of varying parameters are conducted with respect to the plunge stiffness, pitch stiffness, distance of elastic axis from center of gravity, distance of geometric center from elastic axis, load resistance and magnetic flux density. Meanwhile, three indicators are applied to reveal their effects on the harvesting performances: (1) the size of limit cycle domain, (2) the onset velocity of LCO, and (3) the energy output.     

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.     

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.     

Terminal sliding mode control for aeroelastic systems

An aeroelastic system is a nonlinear system with two freedoms, i.e., the plunge displacement and the pitch angle, in a dynamic system model. A chaos effect or a limit cycle oscillation is presumably attributed to the nonlinear effect of the pitch angle mentioned above or the interaction between the aerodynamic behaviors. It is that a single trailing edge input in an aeroelastic system is employed as a way to suppress the limit cycle oscillation with an exclusive choice between the plunge displacement and the pitch angle for a control law design. Consequently, the remaining inevitably turns into an internal dynamics, whose stability is adversely affected by the flight speed and structure parameters, a problem improved by no means using a singe control input design. Toward this end, this work presents a controller design criterion with multiple input channels for both the leading and training edges to remove the uncertainty effect of internal dynamics, and render more room for the response design of the plunge displacement as well as the pitch angle. Mostly due to external disturbance and unknown uncertainty, there is a deviation between the intended and implemented system performances for a robust control design, a mainstream research issue in the modern control. As a consequence of a sliding mode control utilized here, the limit cycle oscillation suffered in an aeroelastic system is removed effectively by the use of a terminal sliding mode control, and the chattering phenomenon seen in the control signal is hence eliminated by his method. It is seen from simulations that the control system is stably assured to reach the target within a limited time frame with an addition of a saturation function to the control law.     

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.     

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.     

Limit cycle oscillation behavior of transonic control surface buzz considering free-play nonlinearity

The limit cycle oscillation (LCO) behaviors of control surface buzz in transonic flow are studied. Euler equations are employed to obtain the unsteady aerodynamic forces for Type B and Type C buzz analyses, and an all-movable control surface model, a wing/control surface model and a three-dimensional wing with a full-span control surface are adopted in the study. Aerodynamic and structural describing functions are used to deal with aerodynamic and structural nonlinearities, respectively. Then the buzz speed and buzz frequency are obtained by V-g method. The LCO behavior of the transonic control surface buzz system with linear structure exhibits subcritical or supercritical bifurcation at different Mach numbers. For nonlinear structural model with a free-play nonlinearity in the control surface deflection stiffness, the double LCO phenomenon is observed in certain range of flutter speed. The free-play nonlinearity changes the stability of LCOs at small amplitudes and turns the unstable LCO into a stable one. The LCO behavior is dominated by the aerodynamic nonlinearity for the case with large control surface oscillation amplitude but by the structural nonlinearity for the case with small amplitude. Good agreements between LCO behaviors obtained by the present method and available experimental data show that our study may help to explain the experimental observation in wind tunnel tests and to understand the physical mechanism of transonic control surface buzz.     

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.     

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.     

LCO flutter of cantilevered woven glass/epoxy laminate in subsonic flow

The paper presents a cantilevered composite wing, aeroelastic characteristics of idealized as a composite flat plate laminate. The composite laminate was made from woven glass fibers with epoxy matrix. The elastic and dynamic properties of the laminate were determined experimentally for aeroelastic calculations. Aeroelastic wind tunnel testing of the laminate was performed and the result showed that flutter, a dynamic instability occurred. The cantilevered laminate also displayed limit cycle amplitude, post-flutter oscillation. The experimental flutter velocity and frequency were verified by our computational analysis.     

Time Average Geometric Moiré—Back to the Basics

Applicability of time average geometric moiré for elastic oscillating structures is analysed in this paper. Mathematical and numerical models describing the formation of time averaged fringes are carefully constructed without the assumption that dynamic deflections are described by a slowly varying function. Though time average geometric moiré is considered as a classical optical experimental technique, we show that well known relationship between the fringe order, amplitude of oscillation and pitch of the grating in state of equilibrium can be used only when the amplitude is small. Otherwise the inverse problem of fringe interpretation becomes much more complicated and is the object of analysis in this paper. We describe the interpretation of fringes produced by time average geometric moiré in detail and illustrate the complexity of the problem by numerical examples.     

Computational aeroelastic investigation of a transonic limit-cycle-oscillation experiment at a transport aircraft wing model

Aeroelastic measurements of a three-dimensional wing model, the so-called Aerostabil wing, were conducted in the Transonic Windtunnel Göttingen. This clean, backward-swept wing allowed the experimental investigation of limit cycle oscillations in a certain transonic parameter range. In this paper, a detailed insight into the observed physical phenomena, especially the measured limit cycle oscillations, is presented by means of CFD–CSM coupled simulations. These simulations on the basis of a detailed structural finite element model reveal the specific properties of the Aerostabil wing and furthermore allow investigating the unstable behavior of this windtunnel model for transonic flow settings. The aerodynamic characteristics include a two-shock system and large flow separation areas, further increasing the complexity of the aeroelastic problem. A structural single degree-of-freedom system is used for the prediction of the experimental stability range and the limit cycle oscillation investigations. Due to the good agreement of simulation and experiment the limit cycle oscillations can be explained by means of nonlinear aerodynamic effects.     

Nonlinear aeroelasticity of high-aspect-ratio wings excited by time-dependent thrust

Effects of engine placement on flutter characteristics of a very flexible high-aspect-ratio wing are investigated using the code NATASHA (Nonlinear Aeroelastic Trim And Stability of HALE Aircraft). Gravity for this class of wings plays an important role in flutter characteristics. In the absence of aerodynamic and gravitational forces and without an engine, the kinetic energy of the first two modes are calculated. Maximum and minimum flutter speed locations coincide with the area of minimum and maximum kinetic energy of the second bending and torsion modes. Time-dependent dynamic behavior of a turboshaft engine (JetCat SP5) is simulated with a transient engine model and the nonlinear aeroelastic response of the wing to the engine’s time-dependent thrust and dynamic excitation is presented. Below the flutter speed, at the wing tip and behind the elastic axis, the impulse engine excitation leads to a stable limit cycle oscillation; and for the ramp kind of excitation, beyond the flutter speed, at 75 % span, behind the elastic axis, it produces chaotic oscillation in the wing. Both the excitations above the flutter speed are stabilized, inboard of the wing.     

Experimental study of wings undergoing active root flapping and pitching

This paper presents the results of experiments carried out on mechanical wings undergoing active root flapping and pitching in the wind tunnel. The objective of the work is to investigate the effect of the pitch angle oscillations and wing profile on the aerodynamic forces generated by the wings. The experiments were repeated for a different reduced frequency, airspeed, flapping and pitching kinematics, geometric angle of attack and wing sections (one symmetric and two cambered airfoils). A specially designed mechanical flapper was used, modelled on large migrating birds. It is shown that, under pitch leading conditions, good thrust generation can be obtained at a wide range of Strouhal numbers if the pitch angle oscillation is adjusted accordingly. Consequently, high thrust was measured at both the lowest and highest tested Strouhal numbers. Furthermore, the work demonstrates that the aerodynamic forces can be sensitive to the Reynolds number, depending on the camber of the wings. Under pitch lagging conditions, where the effective angle of attack amplitude is highest, the symmetric wing was affected by the Reynolds number, generating less thrust at the lowest tested Reynolds value. In contrast, under pure flapping conditions, where the effective angle of attack amplitude was lower but still significant, it was the cambered wings that demonstrated Reynolds sensitivity.     

Aeroelastic coupling of geometrically nonlinear structures and linear unsteady aerodynamics: Two formulations

Two different time domain formulations of integrating commonly used frequency-domain unsteady aerodynamic models based on a modal approach with full order finite element models for structures with geometric nonlinearities are presented. Both approaches are tailored to flight vehicle configurations where geometric stiffness effects are important but where deformations are moderate, flow is attached, and linear unsteady aerodynamic modeling is adequate, such as low aspect ratio wings or joined-wing and strut-braced wings at small to moderate angles of attack. Results obtained using the two approaches are compared using both planar and non-planar wing configurations. Sub-critical and post-flutter speeds are considered. It is demonstrated that the two methods lead to the same steady solution for the sub-critical case after the transients subside. It is also shown that the two methods predict the amplitude and frequency of limit cycle oscillation (when present) with the same accuracy.     

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.     

The effect of phase angle and wing spacing on tandem flapping wings

In a tandem wing configuration, the hindwing often operates in the wake of the forewing and, hence, its performance is affected by the vortices shed by the forewing. Changes in the phase angle between the flapping motions of the fore and the hind wings, as well as the spacing between them, can affect the resulting vortex/wing and vortex/vortex interactions. This study uses 2D numerical simulations to investigate how these changes affect the leading dege vortexes (LEV) generated by the hindwing and the resulting effect on the lift and thrust coefficients as well as the efficiencies. The tandem wing configuration was simulated using an incompressible Navier-Stokes solver at a chord-based Reynolds number of 5 000. A harmonic single frequency sinusoidal oscillation consisting of a combined pitch and plunge motion was used for the flapping wing kinematics at a Strouhal number of 0.3. Four different spacings ranging from 0.1 chords to 1 chord were tested at three different phase angles, 0°, 90° and 180°. It was found that changes in the spacing and phase angle affected the timing of the interaction between the vortex shed from the forewing and the hindwing. Such an interaction affects the LEV formation on the hindwing and results in changes in aerodynamic force production and efficiencies of the hindwing. It is also observed that changing the phase angle has a similar effect as changing the spacing. The results further show that at different spacings the peak force generation occurs at different phase angles, as do the peak efficiencies.     

Scaling of flow separation on a pitching low aspect ratio plate

We use two different dye injection approaches, in two different water tunnels, to visualize the formation and subsequent evolution of leading-edge vortices and related separated structures, for a pitching low aspect ratio plate. The motion is a smoothed linear pitch ramp from 0° to 40° incidence, brief hold, and return to 0°, executed at reduced pitch rates ranging from 0.1 to 0.35 and about various pivot locations. All cases evince a leading edge vortex with pronounced axial flow, which leads to formation of large-scale, three-dimensional flow structures, culminating in a large vortical structure centered at the wing symmetry plane. Pitch is also compared to plunge, where the plunge-induced angle of attack is taken as the geometric pitch incidence angle, ignoring pitch-rate effects. At successively increasing values of convective time C/U, the three-dimensional patterns of the flow structure are remarkably similar for the pitching and plunging motions. The similarity of these patterns persists, though they are shifted in time, for variation of either the location of the pitching axis or the dimensionless pitch rate.     

THE FLOW-INDUCED VIBRATION OF A FLEXIBLE STRIP HANGING VERTICALLY IN A PARALLEL FLOW PART 1: TEMPORAL AEROELASTIC INSTABILITY

The effect of several parameters in a fluid–strip system are studied for linear/nonlinear models in detail. Such parameters are: the number of modes in the Galerkin discretization, the length of the strip, and the flow velocity. The present simulation clearly shows that when nonlinear forces are considered the response approaches a flutter-type limit cycle at supercritical flow speeds. The Reynolds number at which the strip begins to oscillate is about 10⁴–10⁵. With a further increase of the flow velocity the strip oscillates regularly. At higher flow velocities the oscillation becomes violent and irregular. The amplitude, frequency and drag coefficient at the limit cycle are presented as functions of the flow velocity for a given strip. The numerical predictions are in qualitative agreement with previous experimental data.     

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.