Transonic and Supersonic Flutter Trend Investigation of a Variable-sweep Wing
Wind tunnel analysis of flutter trend for variable sweep aircraft at transonic and supersonic speeds
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Wind tunnel analysis of flutter trend for variable sweep aircraft at transonic and supersonic speeds
The flutter characteristics of several wings with an aspect-ratio of 4.0, a taper ratio of 0.2, and a quarter-chord sweepback of 45 deg. have been investigated analytically for Mach numbers up to 2.0. The calculations were based on the modified-strip-analysis method, the subsonic-kernel-function method, piston theory, and quasi-steady second-order theory. Results of t h e analysis and comparisons with experiment indicated that: (1) Flutter speeds were accurately predicted by the modified strip analysis, although accuracy at t h e highest Mach numbers required the use of nonlinear aerodynamic theory (which accounts for effects of wing thickness) for the calculation of the aerodynamic parameters. (2) An abrupt increase of flutter-speed coefficient with increasing Mach number, observed experimentally in the transonic range, was also indicated by the modified strip analysis. (3) In the low supersonic range for some densities, a discontinuous variation of flutter frequency with Mach number was indicated by the modified strip analysis. An abrupt change of frequency appeared experimentally in the transonic range. (4) Differences in flutter-speed-coefficient levels obtained from tests at low supersonic Mach numbers in two wind tunnels were also predicted by the modified strip analysis and were shown to be caused primarily by differences in mass ratio. (5) Flutter speeds calculated by the subsonic-kernel-function method were in good agreement with experiment and with the results of the modified strip analysis. (6) Flutter speed obtained from piston theory and from quasi-steady second-order theory were higher than experimental values by at least 38 percent.
A modal flutter analysis for biaxially loaded, orthotropic panels, using linear piston-theory aerodynamics, was extended in order to include the effects of inplane shear loading. Flutter boundaries for shear loads up to buckling are calculated for simply supported, isotropic panels of various length-width ratios and for a square, isotropic panel with elastic boundary conditions along the leading and trailing edges. These flutter boundaries are used to define conservative design curves. Sample calculations made using these design curves indicate that practical panels, which have otherwise been adequately designed, could become flutter critical if the inplane shear loads approach the buckling value.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The utilization and operation of a set of six computer programs for the prediction of panel flutter at supersonic speeds by finite element methods are described. The programs run individually to determine the flutter behavior of a flat panel where the finite elements which model the panel each have four degrees of freedom (DOF), a curved panel where the finite elements each have four DOF, and a curved panel where the finite elements each have six DOF. The panels are assumed to be of infinite aspect ratio and are subjected to either simply-supported or clamped boundary conditions. The aerodynamics used by these programs are based on piston theory. Application of the program is illustrated by sample cases where the number of beam finite elements equals four, the in-plane tension parameter is 0.0, the maximum camber to panel length ratio for a curved panel case is 0.05, and the Mach number is 2.0. This memorandum provides a user's guide for these programs, describes the parameters that are used, and contains sample output from each of the programs.
Flutter analyses were developed to predict the onset of supersonic unstalled flutter of a cascade of two-dimensional airfoils. The first of these analyzes the onset of supersonic flutter at low levels of aerodynamic loading (i.e., backpressure), while the second examines the occurrence of supersonic flutter at moderate levels of aerodynamic loading. Both of these analyses are based on the linearized unsteady inviscid equations of gas dynamics to model the flow field surrounding the cascade. These analyses are utilized in a parametric study to show the effects of cascade geometry, inlet Mach number, and backpressure on the onset of single and multi degree of freedom unstalled supersonic flutter. Several of the results are correlated against experimental qualitative observation to validate the models.
Recently two flutter analyses have been developed at NASA Lewis Research Center to predict the onset of supersonic unstalled flutter of a cascade of two-dimensional airfoils. The first of these analyzes the onset of supersonic flutter at low levels of aerodynamic loading (i.e., backpressure), while the second examines the occurrence of supersonic flutter at moderate levels of aerodynamic loading. Both of these analyses are based on the linearized unsteady inviscid equations of gas dynamics to model the flow field surrounding the cascade. The details of the development of the solution to each of these models have been published. The objective of the present paper is to utilize these analyses in a parametric study to show the effects of cascade geometry, inlet Mach number, and backpressure on the onset of single and multi degree of freedom unstalled supersonic flutter. Several of the results from this study are correlated against experimental qualitative observation to validate the models.
An analytical model is developed for predicting the onset of supersonic stall bending flutter in axial flow compressors. The analysis is based on a modified two dimensional, compressible, unsteady actuator disk theory. It is applied to a rotor blade row by considering a cascade of airfoils whose geometry and dynamic response coincide with those of a rotor blade element at 85 percent of the span height (measured from the hub). The rotor blades are assumed to be unshrouded (i.e., free standing) and to vibrate in their first flexural mode. The effects of shock waves and flow separation are included in the model through quasi-steady, empirical, rotor total-pressure-loss and deviation-angle correlations. The actuator disk model predicts the unsteady aerodynamic force acting on the cascade blading as a function of the steady flow field entering the cascade and the geometry and dynamic response of the cascade. Calculations show that the present model predicts the existence of a bending flutter mode at supersonic inlet Mach numbers. This flutter mode is suppressed by increasing the reduced frequency of the system or by reducing the steady state aerodynamic loading on the cascade. The validity of the model for predicting flutter is demonstrated by correlating the measured flutter boundary of a high speed fan stage with its predicted boundary. This correlation uses a level of damping for the blade row (i.e., the log decrement of the rotor system) that is estimated from the experimental flutter data. The predicted flutter boundary is shown to be in good agreement with the measured boundary.
Supersonic asymmetric flutter and divergence characteristics of truncated conical shells with ring-supported edges under three types of flow-induced instabilities
Lifting surfaces supersonic-hypersonic flutter at angle of attack determined by shock expansion, Newtonian flow and local flow piston theory
An analytical model is proposed for predicting the onset of supersonic stall bending flutter in high-speed rotors. The analysis is based on a modified two-dimensional, compressible, unsteady actuator disk theory. The stability boundary predicted by the analysis is shown to be in good agreement with the measured boundary of a high speed fan. The prediction that the flutter mode would be a forward traveling wave sensitive to wheel speed and aerodynamic loading is confirmed by experimental measurements. In addition, the analysis shows that reduced frequency and dynamic head also play a significant role in establishing the supersonic stall bending flutter boundary of an unshrouded fan.
Subsonic and supersonic flutter calculations on swept and unswept wings to evaluate effects of variations in aerodynamic parameters and density
An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.
Experiment and theory of supersonic flutter in panels - length/width and stress ratios to prevent flutter of various thicknesses at mach 3
A numerical technique is presented that determines the minimum-weight thickness distribution for a one-dimensional, simply supported panel subjected to a parallel supersonic flow on one side. An aeroelastic eigenvalue, which characterizes flutter speed, is held constant. The governing differential equations are approximated by sets of difference equations adjoined to the weight function via a penalty function. A conjugate gradient method then solves the resulting sequence of unconstrained minimization problems. Numerical results are obtained without imposing a minimum thickness constraint; the effect of constant inplane tensile and compressive stresses is also investigated.
Flutter tests have been made on flat panels having a 1/4 inch-thick plastic-foam core covered with thin fiber-glass laminates. The testing was done in the Langley Unitary Plan wind tunnel at Mach numbers from 1.76 t o 2.87. The flutter boundary for these panels was found to be near the flutter boundary of thin metal panels when compared on the basis of an equivalent panel stiffness. The results also demonstrated that the depth of the cavity behind the panel has a pronounced influence on flutter. Changing the cavity depth from 1 1/2 inches to 1/2 inch reduced the dynamic pressure at start of flutter by 40 percent. No flutter was obtained when the spacers on the back of the panel were against the bottom of the cavity.
The flutter characteristics of a generic arrow-wing supersonic transport configuration are studied. The wing configuration has a 3 percent biconvex airfoil and a leading-edge sweep of 73 deg out to a cranked tip with a 60 deg leading-edge sweep. The ground vibration tests and flutter test procedure are described. The effects of flutter on engine nacelles, fuel loading, wing-mounted vertical fin, wing angle-of-attack, and wing tip mass and stiffness distributions are analyzed. The data reveal that engine nacelles reduce the transonic flutter dynamic pressure by 25-30 percent; fuel loadings decrease dynamic pressures by 25 percent; 4-6 deg wing angles-of-attack cause steep transonic boundaries; and 5-10 percent changes in flutter dynamic pressures are the result of the wing-mounted vertical fin and wing-tip mass and stiffness distributions.