Approximate Floquet Analysis of Parametrically Excited Multi-Degree-of-Freedom Systems With Application to Wind Turbines

2018 ◽  
Vol 141 (1) ◽  
Author(s):  
Gizem D. Acar ◽  
Brian F. Feeny

General responses of multi-degrees-of-freedom (MDOF) systems with parametric stiffness are studied. A Floquet-type solution, which is a product between an exponential part and a periodic part, is assumed, and applying harmonic balance, an eigenvalue problem is found. Solving the eigenvalue problem, frequency content of the solution and response to arbitrary initial conditions are determined. Using the eigenvalues and the eigenvectors, the system response is written in terms of “Floquet modes,” which are nonsynchronous, contrary to linear modes. Studying the eigenvalues (i.e., characteristic exponents), stability of the solution is investigated. The approach is applied to MDOF systems, including an example of a three-blade wind turbine, where the equations of motion have parametric stiffness terms due to gravity. The analytical solutions are also compared to numerical simulations for verification.

Author(s):  
Theodore S. Brockett ◽  
Jerzy T. Sawicki

A six-degree-of-freedom non-linear model is developed using Lagrange’s equation. The model is used to estimate transient fan-stage dynamic response during a fan-blade-out event in a turbo fan engine. The coupled degrees of freedom in the model include the fan whirl in the fan plane, the torsional response of the fan and low-pressure turbines (LPTs) about the engine centerline, the radial position of the released blade fragment, and the angular rotation of the trailing blade from its free state due to acceleration of the released blade. The released blade is assumed to slide radially outward along the trailing blade without friction. The external loading applied to the system includes fan imbalance, the remaining fan blades machining away the rub strip, rubbing of the blades with the fan case, and slowly-varying torques on the low pressure (LP) spool as engine performance degrades. The machining of the abradable imparts tangential loading on the fan blades as momentum is transferred to the liberated rub strip material. After application of the initial conditions including angular positions, angular velocities, released blade fragment position, and torsional wind-up, the governing equations are integrated forward in time from the instant the blade fragment is released. A reasonable match to test data is shown. Parameters affecting the fan-system response are varied to study the impact on fan peak lateral whirl amplitude, peak LP shaft torque, and peak loading on the trailing blade. It is found that the rub strip and mass eccentricity have the strongest influence on the LP shaft torsional loading. It is found that mass eccentricity has the largest influence on peak fan whirl. It is also found that released blade mass and attachment stiffness have the largest influence on the trailing blade loading.


Author(s):  
Renan F. Corrêa ◽  
Flávio D. Marques

Abstract Aeroelastic systems have nonlinearities that provide a wide variety of complex dynamic behaviors. Nonlinear effects can be avoided in practical applications, as in instability suppression or desired, for instance, in the energy harvesting design. In the technical literature, there are surveys on nonlinear aeroelastic systems and the different manners they manifest. More recently, the bistable spring effect has been studied as an acceptable nonlinear behavior applied to mechanical vibration problems. The application of the bistable spring effect to aeroelastic problems is still not explored thoroughly. This paper contributes to analyzing the nonlinear dynamics of a typical airfoil section mounted on bistable spring support at plunging motion. The equations of motion are based on the typical aeroelastic section model with three degrees-of-freedom. Moreover, a hardening nonlinearity in pitch is also considered. A preliminary analysis of the bistable spring geometry’s influence in its restoring force and the elastic potential energy is performed. The response of the system is investigated for a set of geometrical configurations. It is possible to identify post-flutter motion regions, the so-called intrawell, and interwell. Results reveal that the transition between intrawell to interwell regions occurs smoothly, depending on the initial conditions. The bistable effect on the aeroelastic system can be advantageous in energy extraction problems due to the jump in oscillation amplitudes. Furthermore, the hardening effect in pitching motion reduces the limit cycle oscillation amplitudes and also delays the occurrence of the snap-through.


Universe ◽  
2018 ◽  
Vol 4 (9) ◽  
pp. 95 ◽  
Author(s):  
Gianluca Calcagni

We present a method to solve the nonlinear dynamical equations of motion in gravitational theories with fundamental nonlocalities of a certain type. For these specific form factors, which appear in some renormalizable theories, the number of field degrees of freedom and of initial conditions is finite.


Author(s):  
Daniele Dessi ◽  
Franco Mastroddi

A three degrees of freedom aeroelastic typical section with control surface is theoretically modeled including nonlinear springs and augmented states for linear unsteady aerodynamic description. The system response is determined by time marching of the governing equations by using a standard Runge-Kutta algorithm in conjunction with a ‘shooting method’ to find out stable and unstable limit cycles along with stability reversal in the neighborhood of the Hopf bifurcation. Furthermore, the equations of motion are analyzed by a singular perturbation technique, specifically, by using a normal form method. Approximate analytical expressions for amplitudes and frequencies of limit cycles are obtained and the terms which are responsible of the nonlinear system behavior are identified.


2018 ◽  
Vol 18 (06) ◽  
pp. 1850085 ◽  
Author(s):  
Vahid Jahangiri ◽  
Mir Mohammad Ettefagh

Stability and dynamic modeling of the floating wind turbine (FWT) is a crucial challenge in designing of the type of structures. In this paper, the tension leg platform (TLP) type FWT is modeled as a multibody system considering the flexibility between the nacelle and tower. The flexibility of the FWT is modeled as a torsional spring and damper. It has 6 degrees of freedom (DOFs) related to the large-amplitude translation and rotation of the tower and 4 DOFs related to the relative rotation between the rotor-nacelle assembly and the tower. First, the nonlinear equations of motion are derived by the theory of momentum cloud based on the conservation of momentum. Then, the equations of motion are solved and the system is simulated in MATLAB. Moreover, the effect of flexibility between the nacelle and tower is investigated via the dynamic response. The stability of the system in three different environmental conditions is studied. Finally, the spring and damping coefficients for the system response to get near to instability are determined, by which the critical region is defined. The simulation results demonstrate the importance of the flexibility between the nacelle and tower on the overall behavior of the system and its stability.


Author(s):  
André de Souza Mendes ◽  
Agenor de Toledo Fleury ◽  
Marko Ackermann ◽  
Fabrizio Leonardi ◽  
Roberto Bortolussi

This article addresses the yaw stability of articulated vehicles by assessing the influence of the road-tire friction coefficient on the convergence region of a particular equilibrium condition. In addition, the boundaries of this region are compared to the boundaries of the non-jackknife and non-rollover regions to distinguish the instability phenomenon, jackknife or roll-over, responsible for this delimitation. The vehicle configuration considered in this analysis is composed by one tractor unit and one towed unit connected through an articulation point, for instance, a tractor-semitrailer combination. A nonlinear articulated bicycle model with four degrees of freedom is used together with a nonlinear lateral force tire model. To estimate the convergence region, the phase trajectory method is used. The equations of motion of the mathematical model are numerically integrated for different initial conditions in the phase plane, and the state orbits are monitored in order to verify the convergence point and the occurrence of instability events. In all cases, the longitudinal force on each tire, such as traction and braking, is not considered. The results show the existence of convergence regions delimited only by jackknife events, for low values of the friction coefficient, and only by rollover events, for high values of the friction coefficient. Moreover, the transition between these two conditions as the friction coefficient is changed is graphically presented. The main contributions of this article are the identification of the abrupt reduction of the convergence region as the value of the friction coefficient increases and the distinction of the instability events, jackknife or rollover, that define the boundaries of the convergence region.


2020 ◽  
Vol 25 (4) ◽  
pp. 116-129
Author(s):  
O.S. Lanets ◽  
V.T. Dmytriv ◽  
V.M. Borovets ◽  
I.A. Derevenko ◽  
I.M. Horodetskyy

AbstractThe article deals with atwo-mass above resonant oscillatory system of an eccentric-pendulum type vibrating table. Based on the model of a vibrating oscillatory system with three masses, the system of differential equations of motion of oscillating masses with five degrees of freedom is compiled using generalized Lagrange equations of the second kind. For given values of mechanical parameters of the oscillatory system and initial conditions, the autonomous system of differential equations of motion of oscillating masses is solved by the numerical Rosenbrock method. The results of analytical modelling are verified by experimental studies. The two-mass vibration system with eccentric-pendulum drive in resonant oscillation mode is characterized by an instantaneous start and stop of the drive without prolonged transient modes. Parasitic oscillations of the working body, as a body with distributed mass, are minimal at the frequency of forced oscillations.


1966 ◽  
Vol 25 ◽  
pp. 373
Author(s):  
Y. Kozai

The motion of an artificial satellite around the Moon is much more complicated than that around the Earth, since the shape of the Moon is a triaxial ellipsoid and the effect of the Earth on the motion is very important even for a very close satellite.The differential equations of motion of the satellite are written in canonical form of three degrees of freedom with time depending Hamiltonian. By eliminating short-periodic terms depending on the mean longitude of the satellite and by assuming that the Earth is moving on the lunar equator, however, the equations are reduced to those of two degrees of freedom with an energy integral.Since the mean motion of the Earth around the Moon is more rapid than the secular motion of the argument of pericentre of the satellite by a factor of one order, the terms depending on the longitude of the Earth can be eliminated, and the degree of freedom is reduced to one.Then the motion can be discussed by drawing equi-energy curves in two-dimensional space. According to these figures satellites with high inclination have large possibilities of falling down to the lunar surface even if the initial eccentricities are very small.The principal properties of the motion are not changed even if plausible values ofJ3andJ4of the Moon are included.This paper has been published in Publ. astr. Soc.Japan15, 301, 1963.


2020 ◽  
Vol 1 (1) ◽  
pp. 93-102
Author(s):  
Carsten Strzalka ◽  
◽  
Manfred Zehn ◽  

For the analysis of structural components, the finite element method (FEM) has become the most widely applied tool for numerical stress- and subsequent durability analyses. In industrial application advanced FE-models result in high numbers of degrees of freedom, making dynamic analyses time-consuming and expensive. As detailed finite element models are necessary for accurate stress results, the resulting data and connected numerical effort from dynamic stress analysis can be high. For the reduction of that effort, sophisticated methods have been developed to limit numerical calculations and processing of data to only small fractions of the global model. Therefore, detailed knowledge of the position of a component’s highly stressed areas is of great advantage for any present or subsequent analysis steps. In this paper an efficient method for the a priori detection of highly stressed areas of force-excited components is presented, based on modal stress superposition. As the component’s dynamic response and corresponding stress is always a function of its excitation, special attention is paid to the influence of the loading position. Based on the frequency domain solution of the modally decoupled equations of motion, a coefficient for a priori weighted superposition of modal von Mises stress fields is developed and validated on a simply supported cantilever beam structure with variable loading positions. The proposed approach is then applied to a simplified industrial model of a twist beam rear axle.


2017 ◽  
Vol 139 (4) ◽  
Author(s):  
Samuel F. Asokanthan ◽  
Soroush Arghavan ◽  
Mohamed Bognash

Effect of stochastic fluctuations in angular velocity on the stability of two degrees-of-freedom ring-type microelectromechanical systems (MEMS) gyroscopes is investigated. The governing stochastic differential equations (SDEs) are discretized using the higher-order Milstein scheme in order to numerically predict the system response assuming the fluctuations to be white noise. Simulations via Euler scheme as well as a measure of largest Lyapunov exponents (LLEs) are employed for validation purposes due to lack of similar analytical or experimental data. The response of the gyroscope under different noise fluctuation magnitudes has been computed to ascertain the stability behavior of the system. External noise that affect the gyroscope dynamic behavior typically results from environment factors and the nature of the system operation can be exerted on the system at any frequency range depending on the source. Hence, a parametric study is performed to assess the noise intensity stability threshold for a number of damping ratio values. The stability investigation predicts the form of threshold fluctuation intensity dependence on damping ratio. Under typical gyroscope operating conditions, nominal input angular velocity magnitude and mass mismatch appear to have minimal influence on system stability.


Sign in / Sign up

Export Citation Format

Share Document