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.
New method for calculating flow of air through and around helicopter rotor incorporated into General Rotorcraft Aeromechanical Stability Program (GRASP) (computer program for aeroelastic analysis). Flow about helicopter rotor represented by axisymmetric flow field in cylindrical region with actuator disk as source of flow.
This paper describes the calculation of flutter stability characteristics for a transonic forward swept fan configuration using a viscous aeroelastic analysis program. Unsteady Navier-Stokes equations are solved on a dynamically deforming, body fitted, grid to obtain the aeroelastic characteristics using the energy exchange method. The non-zero inter-blade phase angle is modeled using phase-lagged boundary conditions. Results obtained show good correlation with measurements. It is found that the location of shock and variation of shock strength strongly influenced stability. Also, outboard stations primarily contributed to stability characteristics. Results demonstrate that changes in blade shape impact the calculated aerodynamic damping, indicating importance of using accurate blade operating shape under centrifugal and steady aerodynamic loading for flutter prediction. It was found that the calculated aerodynamic damping was relatively insensitive to variation in natural frequency.
Flutter analysis of configurations with geometric structural nonlinearities typically is done with time-domain analysis. The results from this process are computationally expensive and can yield cumbersome results that may be difficult to manage/interpret. Compared to time-domain methods, frequency-domain flutter analysis can provide additional insight into the characteristics of a flutter stability problem. By linearizing the aeroelastic problem about the nonlinear equilibrium state, this work applies frequency-domain aeroelastic analysis to the Pazy wing, the subject of the Large Deformation Working Group in the Aeroelastic Prediction Workshop. Generalized aerodynamic forces (GAFs) are computed with both a doublet-lattice method and a computational fluid dynamics solver at a range of reduced frequencies as well as a range of dynamic pressures to account for the dependence of the mode shapes on the nonlinear equilibrium state. These GAFs are used in a p-k flutter solver, which is modified to handle the dependence of the stiffness matrix and GAFs on the dynamic pressure.
Aeroelastic analyses for advanced turbomachines are being developed for use at the NASA Glenn Research Center and industry. However, these analyses at present are used for turbomachinery design with uncertainties accounted for by using safety factors. This approach may lead to overly conservative designs, thereby reducing the potential of designing higher efficiency engines. An integration of the deterministic aeroelastic analysis methods with probabilistic analysis methods offers the potential to design efficient engines with fewer aeroelastic problems and to make a quantum leap toward designing safe reliable engines. In this research, probabilistic analysis is integrated with aeroelastic analysis: (1) to determine the parameters that most affect the aeroelastic characteristics (forced response and stability) of a turbomachine component such as a fan, compressor, or turbine and (2) to give the acceptable standard deviation on the design parameters for an aeroelastically stable system. The approach taken is to combine the aeroelastic analysis of the MISER (MIStuned Engine Response) code with the FPI (fast probability integration) code. The role of MISER is to provide the functional relationships that tie the structural and aerodynamic parameters (the primitive variables) to the forced response amplitudes and stability eigenvalues (the response properties). The role of FPI is to perform probabilistic analyses by utilizing the response properties generated by MISER. The results are a probability density function for the response properties. The probabilistic sensitivities of the response variables to uncertainty in primitive variables are obtained as a byproduct of the FPI technique. The combined analysis of aeroelastic and probabilistic analysis is applied to a 12-bladed cascade vibrating in bending and torsion. Out of the total 11 design parameters, 6 are considered as having probabilistic variation. The six parameters are space-to-chord ratio (SBYC), stagger angle (GAMA), elastic axis (ELAXS), Mach number (MACH), mass ratio (MASSR), and frequency ratio (WHWB). The cascade is considered to be in subsonic flow with Mach 0.7. The results of the probabilistic aeroelastic analysis are the probability density function of predicted aerodynamic damping and frequency for flutter and the response amplitudes for forced response.
Several aeroelastic stability analyses for a vertical fin similar to that of the National Aero-Space Plane are described. The objectives of the study were to design and obtain an experimental data base for a supersonic wind-tunnel model of the fin in order to examine the effects of thermal loading on the flutter characteristics. This paper describes the preliminary efforts to design the wind-tunnel model, including several of the geometric parameter variations that were analyzed. The dominant flutter mechanism involved a flap vibration mode and a fin bending mode. Variation of the thicknesses of flap and root flexures, used to attach the flap to the fin, and the fin to a support, significantly affected the flutter boundary. Uniform thermal loads, affecting only material properties, had little effect, as did the application of different uniform temperatures to each side of the fin. In contrast, the application of significant chord-wise thermal gradients induced stresses which reduced the flutter dynamic pressure by as much as 37 percent. For less extreme distributed loading, the low-aspect ratio fin was relatively unaffected.
As a collaborative effort among government aerospace research laboratories an advanced version of a widely used computational fluid dynamics code, OVERFLOW, was recently released. This latest version includes additions to model flexible rotating multiple blades. In this paper, the OVERFLOW code is applied to improve the accuracy of airload computations from the linear lifting line theory that uses displacements from beam model. Data transfers required at every revolution are managed through a Unix based script that runs jobs on large super-cluster computers. Results are demonstrated for the 4-bladed UH-60A helicopter. Deviations of computed data from flight data are evaluated. Fourier analysis post-processing that is suitable for aeroelastic stability computations are performed.
A time domain approach is used to determine the dynamic aeroelastic stability of a cascade of blades. The structural model for each blade is a typical section with two degrees of freedom. The aerodynamic model is the unsteady, two-dimensional, full-potential flow through the cascade of airfoils. The unsteady equations of motion for the structure and the fluid are integrated simultaneously in time starting with the steady flowfield and a small initial disturbance applied to the airfoils. The motion of each blade is analyzed to determine the aeroelastic stability of the cascade. The effect of interblade phase angle is included in the analysis by allowing each blade to have an independent motion and considering a number of blade passages. Calculations are made using an airfoil section and structural parameters that are representative of a propfan. The results are compared with those from a separate frequency domain analysis. Good agreement between the results is observed. With the time domain approach, it is possible to consider nonlinear structural models and nonlinear force-displacement relations. The method allows a realistic simulation of the motion of the fluid and the cascade blades for a better physical understanding and it also has the potential for saving computational time when compared to the frequency domain approach for the flutter analysis of cascades.
Flutter prediction is usually a knowledge-based analysis process that aims to reduce the cost of aeroelastic stability margin certification. However, early detection of flutter problems is beneficial in the development of unconventional aircraft. The recently developed automation tool ConceptFEA for structural sizing of aircraft concepts paves the way for rapid physics-based flutter prediction of aircraft concepts. A match-point iteration procedure using the p-k method is implemented for ConceptFEA with minimum user input requirements to generate flutter boundary points. A subsonic business jet concept and its high aspect-ratio wing variant are used to demonstrate how the newly developed flutter prediction capability can be used during aircraft conceptual design. Sized structures, flutter boundary curves, and flutter sensitivity analysis results are generated for these two concepts using ConceptFEA. The relevant equivalent plate theory is provided to show the quantitative relationships between a stiffened panel and its equivalent NASTRAN PSHELL panel. The rapid flutter prediction capability of ConceptFEA makes multidisciplinary collaborations between systems analysts and aeroelasticity experts feasible in practice.
This paper considers complex transcendental eigenvalue problems where one is interested in pairs of eigenvalues that are restricted to take real values only. Such eigenvalue problems arise in dynamic stability analysis of nonconservative physical systems, i.e., flutter analysis of aeroelastic systems. Some available solution methods are discussed and a new method is presented. Two computational approaches are described for analytical evaluation of the sensitivities of these eigenvalues when they are dependent on other parameters. The algorithms presented are illustrated through examples.
A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.
This paper considers complex transcendental eigenvalue problems where one is interested in pairs of eigenvalues that are restricted to take real values only. Such eigenvalue problems arise in dynamic stability analysis of nonconservative physical systems, i.e., flutter analysis of aeroelastic systems. Some available solution methods are discussed and a new method is presented. Two computational approaches are described for analytical evaluation of the sensitivities of these eigenvalues when they are dependent on other parameters. The algorithms presented are illustrated through examples.
An overview of the technical development of hingeless and bearingless rotors is given, with emphasis on aeroelastic and aeromechanical stability characteristics. Important considerations for theoretical analysis are discussed. Theoretical and experimental investigations of isolated blade flap-lag and flap-lag-torsion stability, and coupled rotor-body aeromechanical stability are described. Physical interpretation and important rotor system design parameters are emphasized. An overview of bearingless rotor dynamics is also included.
When designing and certifying aircraft, engineers must take into consideration aeroelastic effects such as flutter. Design and certification of a vehicle may require analysis of thousands of aeroelastic responses. Standard tools in the aerospace industry are based on linear aerodynamic models such as the doublet-lattice method, but these methods can be nonconservative in certain situations such as in the transonic regime. While computational fluid dynamics (CFD) is a higher fidelity alternative, the time-marching approach has a drastically increased computational cost compared to the linear aerodynamic methods. By taking advantage of the periodic nature of flutter, frequency-domain methods offer a more efficient alternative to time-marching CFD. In this work, a linearized frequency-domain method is implemented and verified in the stabilized finite-element solver in FUN3D. The linearized frequency-domain method is demonstrated and compared to other methods for traditional benchmark cases for computational aeroelasticity: the AGARD 445.6 wing, the Benchmark Supercritical Wing, and the Benchmark NACA 0012Wing.
This paper presents details of a flutter and stability analysis of aerospace structures such as hypersonic vehicles. Both structural and aerodynamic domains are discretized by the common finite element technique. A vibration analysis is first performed by the STARS code employing a block Lanczos solution scheme. This is followed by the generation of a linear aerodynamic grid for subsequent linear flutter analysis within subsonic and supersonic regimes of the flight envelope; the doublet lattice and constant pressure techniques are employed to generate the unsteady aerodynamic forces. Flutter analysis is then performed for several representative flight points. The nonlinear flutter solution is effected by first implementing a CFD solution of the entire vehicle. Thus, a 3-D unstructured grid for the entire flow domain is generated by a moving front technique. A finite element Euler solution is then implemented employing a quasi-implicit as well as an explicit solution scheme. A novel multidisciplinary analysis is next effected that employs modal and aerodynamic data to yield aerodynamic damping characteristics. Such analyses are performed for a number of flight points to yield a large set of pertinent data that define flight flutter characteristics of the vehicle. This paper outlines the finite-element-based integrated analysis procedures in detail, which is followed by the results of numerical analyses of flight flutter simulation.
A summary of the work performed during the progress period is presented. Analysis methods for predicting loads and instabilities of wind turbines were developed. Three new areas of research to aid the Advanced Turboprop Project (ATP) were initiated and developed. These three areas of research are aeroelastic analysis methods for cascades including blade and disk flexibility; stall flutter analysis; and computational aeroelasticity.
The lack of polar symmetry in two-bladed rotors leads to equations of motion with periodic coefficients in axial flight, which is contrary to three or more bladed rotors that result in constant coefficient equations. With periodic coefficients, the analysis becomes involved, as a result very few studies have been directed towards the analysis of two-bladed rotors. In this paper, the aeroelastic stability of two-bladed proprotor/pylon/wing combinations is examined in high speed axial flight. Several parametric studies are carried out to illustrate the special nature of two-bladed proprotors and to better understand the mechanism of whirl-flutter in such rotors. The wing beam bending mode for two-bladed rotors is found to be stable over the range of parameters examined, a behaviour very different from three-bladed rotors. Also, the wing torsion mode exhibits a new type of instability similar to a wing torsional divergence scouring at I/rev frequency. This type of behaviour is not seen in three and more bladed rotors. The interaction between wing chordwise bending and torsion modes is found to be much greater in the case of two-bladed rotors and, over the range of parameters considered, these two modes govern the stability of the system.
This investigation explores the effects of chordwise forces and deformations and steady-state deformation due to lift on the static and dynamic aeroelastic stability of a uniform cantilever wing. Results of this analysis are believed to have practical applications for high-performance sailplanes and certain RPV's. The airfoil cross section is assumed to be symmetric and camber bending is neglected. Motions in vertical bending, fore-and-aft bending, and torsion are considered. A differential equation model is developed, which included the nonlinear elastic bending-torsion coupling that accompanies even moderate deflections. A linearized expansion in small time-dependent deflections is made about a steady flight condition. The stability determinant of the linearized system then contains coefficients that depend on steady displacements. Loads derived from two-dimensional incompressible aerodynamic theory are used to obtain the majority of the results, but cases using three-dimensional subsonic compressible theory are also studied. The stability analysis is carried out in terms of the dynamically uncoupled natural modes of vibration of the uniform cantilever.
There is a potential for improving the performance and aeroelastic stability of tiltrotors through the use of elastically-coupled composite rotor blades. To study the characteristics of tiltrotors with these types of rotor blades it is necessary to formulate a new analysis which has the capabilities of modeling both a tiltrotor configuration and an anisotropic rotor blade. Background for these formulations is established in two preliminary investigations. In the first, the influence of several system design parameters on tiltrotor aeroelastic stability is examined for the high-speed axial flight mode using a newly-developed rigid-blade analysis with an elastic wing finite element model. The second preliminary investigation addresses the accuracy of using a one-dimensional beam analysis to predict frequencies of elastically-coupled highly-twisted rotor blades. Important aspects of the new aeroelastic formulations are the inclusion of a large steady pylon angle which controls tilt of the rotor system with respect to the airflow, the inclusion of elastic pitch-lag coupling terms related to rotor precone, the inclusion of hub-related degrees of freedom which enable modeling of a gimballed rotor system and engine drive-train dynamics, and additional elastic coupling terms which enable modeling of the anisotropic features for both the rotor blades and the tiltrotor wing. Accuracy of the new tiltrotor analysis is demonstrated by a comparison of the results produced for a baseline case with analytical and experimental results reported in the open literature. Two investigations of elastically tailored blades on a baseline tiltrotor are then conducted. One investigation shows that elastic bending-twist coupling of the rotor blade is a very effective means for increasing the flutter velocity of a tiltrotor, and the magnitude of coupling required does not have an adverse effect on performance or blade loads. The second investigation shows that passive blade twist control via elastic extension-twist coupling of the rotor blade has the capability of significantly improving tiltrotor aerodynamic performance. This concept, however, is shown to have, in general, a negative impact on stability characteristics.