Flutter and forced response of mistuned rotors using standing wave analysis
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Dugundji, J..
Explore the source record for details and available documents.
Two of the more common procedures for analyzing the stability and forced response of equations with periodic coefficients are reviewed: the use of Floquet methods, and the use of multiblade coordinate and harmonic balance methods. The analysis procedures of these periodic coefficient systems are compared with those of the more familiar constant coefficient systems. Previously announced in STAR as N82-23702
A standing wave approach is applied to the analysis of the flutter and forced response of tuned and mistuned rotors. The traditional traveling wave cascade airforces are recast into standing wave arbitrary motion form using Pade approximants, and the resulting equations of motion are written in the matrix form. Applications for vibration modes, flutter, and forced response are discussed. It is noted that the standing wave methods may prove to be more versatile for dealing with certain applications, such as coupling flutter with forced response and dynamic shaft problems, transient impulses on the rotor, low-order engine excitation, bearing motions, and mistuning effects in rotors.
Horizontal axis wind turbine dynamics were studied. The following findings are summarized: (1) review of the MOSTAS computer programs for dynamic analysis of horizontal axis wind turbines; (2) review of various analysis methods for rotating systems with periodic coefficients; (3) review of structural dynamics analysis tools for large wind turbine; (4) experiments for yaw characteristics of a rotating rotor; (5) development of a finite element model for rotors; (6) development of simple models for aeroelastics; and (7) development of simple models for stability and response of wind turbines on flexible towers.
A horizontal axis wind turbine was constructed to study the effect of coning angle on the yawing moments produced. Coning angles of 0 deg, +10 deg and -10 deg were studied in the upwind and downwind cases. Moment and rotational frequency of the blades at each yaw angle setting were taken. It was found that as the coning angle increased from -10 deg to +10 deg in either the upwind or downwind case the stability decreased. The downwind case was slightly more stable for all coning angles than was the upwind case. It is found that all the previous cases were stable for high rotation speeds, but at lower rotation speeds, they were all unstable and could not self start unless held in the wind.
Two of the more common procedures for analyzing the stability and forced response of equations with periodic coefficients are reviewed: the use of Floquet methods, and the use of multiblade coordinate and harmonic balance methods. The analysis procedures of these periodic coefficient systems are compared with those of the more familiar constant coefficient systems.
The aeroelastic behavior of a transonic compressor rotor operated in the MIT Blowdown Compressor Facility has been examined by means of piezoelectric motion sensors at the base of each of the 23 blades. Excitation has been observed due to rotating stall, due to an incipient flutter, and due to the facility startup transient. A method has been found for determining the aerodynamic damping force by modal analysis of the blade motion. Application of this technique to the example of excitation by rotating stall has led to the conclusions that the blade loading decreases in the stall cell, and that the damping force on the blades in the clean flow is in phase with blade velocity but opposite it in sign, leading to a logarithmic decrement of 0.2. This method of force derivation has quite general applicability as it requires only blade motion data such as are routinely acquired with strain gages. It is argued that models are needed for aerodynamic damping which focus on the effects of near neighbors of a given blade, since flutter often results in large response of isolated blades or small groups of blades.
The paper outlines the nonlinear dynamic analysis of an isolated three-degree flap-lag-feather wind turbine blade under a gravity field and with shear flow. Lagrangian equations are used to derive the nonlinear equations of motion of blade for arbitrarily large angular deflections. The limit cycle analysis for forced oscillations and the determination of the principal parametric resonance of the blade due to periodic forces from the gravity field and wind shear are performed using the harmonic balance method. Results are obtained first for a two-degree flap-lag blade, then the effect of the third degree of freedom (feather) is studied. The self-excited flutter solutions are obtained for a uniform wind and with gravity forces neglected. The effects of several parameters on the blade stability are examined, including coning angle, structural damping, Lock number, and feather frequency. The limit cycle flutter solution of a typical configuration shows a substantial nonlinear softening spring behavior.
A preliminary study conducted by Kerrebrock et al. (1976) has shown that the torsional rigidity of untwisted thin blades of a transonic compressor can be reduced significantly by transient thermal stresses. The aerodynamic loads have various effects on blade vibration. One effect is that gas bending loads may result in a bending-torsion coupling which may change the characteristics of the torsion and bending vibration of the blade. For a general study of transient-temperature distribution within a rotor stage, a finite-element heat-conduction analysis was developed. The blade and shroud are divided into annular elements. With a temperature distribution obtained from the heat-conduction analysis and a prescribed gas bending load distribution along the blade span, the static deformation and moment distributions of the blade can be solved iteratively using the finite-element method. The reduction of the torsional rigidity of pretwisted blades caused by the thermal stress effect is then computed. The dynamic behavior of the blade is studied by a modified Galerkin's method.
Coupled differential equations of motion are given for application to a rotating, pretwisted and heated beam under the effects of thermal stresses and gas bending loads. The circumferential modes of the multi-blade vibration of a bladed-disk rotor stage were studied. A finite element method was developed for the dynamic and static deformation analysis of the blade. The deformations of a bladed disk and a shrouded-bladed disk were studied by introducing a special bladed-disk element and a special shrouded-blade element. Some features of the vibration of part-span-shrouded, bladed-disk rotor stages are discussed. The static deformation, thermal stress and gas bending effects on the blade vibration were presented previously.
An aeroelastic investigation of horizontal axis wind turbines is described. The study is divided into two simpler areas; (1) the aeroelastic stability of a single blade on a rigid tower; and (2) the mechanical vibrations of the rotor system on a flexible tower. Some resulting instabilities and forced vibration behavior are described.
The basic whirl stability of a rotating windmill on a flexible tower is reviewed. Effects of unbalance, gravity force, gyroscopic moments, and aerodynamics are discussed. Some experimental results on a small model windmill are given.
Tests were conducted in a small wind tunnel to study the torsional stall flutter behavior of a two-dimensional flat-plate wing pivoted about the midchord. The nonlinear static divergence equilibrium properties of the wing were well predicted from the measured static moment characteristics. Large amplitude limit cycles ranging from plus or minus 11 degrees to plus or minus 100 degrees were observed. Stall flutter occurred above a critical value of a reduced frequency of about 2. Self-excitation occurred for initial angles of attack between 0 and 8 degrees. Nondimensional harmonic coefficients were extracted from the free transient vibration tests for amplitudes up to 80 degrees.
An experimental investigation is made of the purely torsional stall flutter of a two-dimensional wing pivoted about the midchord, and also of the bending-torsion stall flutter of a two-dimensional wing pivoted about the quarterchord. For the purely torsional flutter case, large amplitude limit cycles ranging from + or - 11 to + or - 160 degrees were observed. Nondimensional harmonic coefficients were extracted from the free transient vibration tests for amplitudes up to 80 degrees. Reasonable nondimensional correlation was obtained for several wing configurations. For the bending-torsion flutter case, large amplitude coupled limit cycles were observed with torsional amplitudes as large as + or - 40 degrees. The torsion amplitudes first increased, then decreased with increasing velocity. Additionally, a small amplitude, predominantly torsional flutter was observed when the static equilibrium angle was near the stall angle.
A systematic way of applying both perturbation methods and harmonic balance methods to nonlinear panel flutter problems is developed here. Results obtained by both these methods for two-dimensional simply supported and three-dimensional clamped-clamped plates with six modes agree well with those obtained by the straightforward direct integration method, yet require less computer time and provide better insight into the solutions. Effects of viscoelastic structural damping on the flutter stability boundary are generally found to be destabilizing and the postflutter behavior becomes more explosive. The methods developed here may be of interest in related vibration problems.
Nonlinear panel flutter analyses by perturbation and harmonic balance methods