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.
Wake model and computer program to compute geometries, flows and velocity influence coefficients for helicopter blade load calculations
This User's Guide applies to the three dimensional viscous flow forward marching analysis, PEPSIG, as used for the calculation of the helicopter tip vortex flow field. The guide presents a discussion of the program flow and subroutines, as well as a list of sample input and output.
A numerical method for solving the Euler equations for multiblade rotors has been developed and some preliminary results reported. The numerical scheme is a combination of several recent methods and algorithm improvements, adapted to the particular requirements of rotor-body interactions. A cylindrical basic grid has been used to study conventional multiblade helicopter rotors. Test calculations have been made for two- and six-blade rotors in hover and for a two-blade rotor in forward flight, under transonic tip conditions but without lift. The results show good agreement with experimental data.
Recently, curved vortex elements have been developed for use in helicopter free wake calculations. As previously documented, these elements allow more efficient and accurate calculation of the wake-on-wake interactions than do traditional straight vortex elements. Curved elements have been applied here to the prediction of the rotor wake geometry and its associated inflow distribution in low speed forward flight. By coupling the wake model to a blade dynamic analysis, a simulation capable of predicting blade motion and hub moments in forward flight has been developed. The resulting code has been validated by successful predictions of the lateral flapping of articulated rotors in low speed level flight, and additional simulators of climb and descent have been undertaken. Successful correlation of other experimental data, including moments due to cyclic on hingeless rotors and transient thrust response in hover, has given encouraging evidence of the capabilities of this code as a tool for analysis of the low speed regime.
A numerical method for solving the Euler equations for multiblade rotors has been developed and some preliminary results reported. The numerical scheme is a combination of several recent methods and algorithm improvements, adapted to the particular requirements of rotor-body interactions. A cylindrical basic grid has been used to study conventional multiblade helicopter rotors. Test calculations have been made for two- and six-blade rotors in hover and for a two-blade rotor in forward flight, under transonic tip conditions but without lift. The results show good agreement with experimental data.
The empirical relation between the induced velocity, thrust, and rate of vertical descent of a helicopter rotor was calculated from wind tunnel force tests on four model rotors by the application of blade-element theory to the measured values of the thrust, torque, blade angle, and equivalent free-stream rate of descent. The model tests covered the useful range of C(sub t)/sigma(sub e) (where C(sub t) is the thrust coefficient and sigma(sub e) is the effective solidity) and the range of vertical descent from hovering to descent velocities slightly greater than those for autorotation. The three bladed models, each of which had an effective solidity of 0.05 and NACA 0015 blade airfoil sections, were as follows: (1) constant-chord, untwisted blades of 3-ft radius; (2) untwisted blades of 3-ft radius having a 3/1 taper; (3) constant-chord blades of 3-ft radius having a linear twist of 12 degrees (washout) from axis of rotation to tip; and (4) constant-chord, untwisted blades of 2-ft radius. Because of the incorporation of a correction for blade dynamic twist and the use of a method of measuring the approximate equivalent free-stream velocity, it is believed that the data obtained from this program are more applicable to free-flight calculations than the data from previous model tests.
Mathematical models based on the finite element method of structural analysis, as embodied in the NASTRAN computer code, are routinely used by the helicopter industry to calculate airframe static internal loads used for sizing structural members. Historically, less reliance has been placed on the vibration predictions based on these models. Beginning in the early 1980's NASA's Langley Research Center initiated an industry wide program with the objective of engendering the needed trust in vibration predictions using these models and establishing a body of modeling guides which would enable confident future prediction of airframe vibration as part of the regular design process. Emphasis in this paper is placed on the successful modeling of the Army/Boeing CH-47D which showed reasonable correlation with test data. A principal finding indicates that improved dynamic analysis requires greater attention to detail and perhaps a finer mesh, especially the mass distribution, than the usual stress model. Post program modeling efforts show improved correlation placing key modal frequencies in the b/rev range with 4 percent of the test frequencies.
The empirical relation between the induced velocity, thrust, and rate of vertical descent of a helicopter rotor was calculated from wind tunnel force tests on four model rotors by the application of blade-element theory to the measured values of the thrust, torque, blade angle, and equivalent free-stream rate of descent. The model tests covered the useful range of C(sub t)/sigma(sub e) (where C(sub t) is the thrust coefficient and sigma(sub e) is the effective solidity) and the range of vertical descent from hovering to descent velocities slightly greater than those for autorotation. The three bladed models, each of which had an effective solidity of 0.05 and NACA 0015 blade airfoil sections, were as follows: (1) constant-chord, untwisted blades of 3-ft radius; (2) untwisted blades of 3-ft radius having a 3/1 taper; (3) constant-chord blades of 3-ft radius having a linear twist of 12 degrees (washout) from axis of rotation to tip; and (4) constant-chord, untwisted blades of 2-ft radius. Because of the incorporation of a correction for blade dynamic twist and the use of a method of measuring the approximate equivalent free-stream velocity, it is believed that the data obtained from this program are more applicable to free-flight calculations than the data from previous model tests.
Boeing Helicopters, together with other U.S. Helicopter manufacturers, participated in a finite element applications program to give the United States a superior capability to utilize finite element analysis models in support of helicopter airframe design. The program was sponsored by the NASA Langley Research Center. Under this program, an activity was sponsored to evaluate existing analysis methods applicable to calculate coupled rotor-airframe vibrations. The helicopter used in this evaluation was the AH-1G helicopter. The results of the Boeing Helicopters efforts are summarized. The planned analytical procedure is reviewed. Changes to the planned procedure are discussed, and results of the correlation study are presented.
Analytical investigation of three dimensional turbulent boundary layers on helicopter rotors including cross flow derivatives and effects of centrifugal and Coriolis forces
Helicopter airframe vibration is examined using calculations and measurements for the SA349/2 research helicopter. The hub loads, which transmit excitations to the fuselage, are predicted using a comprehensive rotorcraft analysis and correlated with measuring hub loads. The predicted and measured hub loads are then coupled with finite element models representing the SA349/2 fuselage. The resulting vertical acceleration at the pilot seat is examined. Adjustments are made to the airframe structural models to examine the sensitivity of predicted vertical acceleration to the model. Changes of a few percent to the damping and frequency of specific models lead to large reductions in predicted vibration, and to major improvements in the correlations with measured pilot-seat vertical acceleration.
Helicopter airframe vibration is examined using calculation and measurements for the SA349/2 research helicopter. The hub loads, which transmit excitation to the fuselage, are predicted using a comprehensive rotorcraft analysis and correlated with measured hub loads. The predicted and measured hub loads are then coupled with finite element models representing the SA349/2 fuselage. The resulting vertical acceleration at the pilot seat is examined. Adjustments are made to the airframe structural models to examine the sensitivity of predicted vertical acceleration to the model. Changes of a few percent to the damping and frequency of specific modes lead to large reductions in predicted vibration and to major improvements in the correlations with measured pilot seat vertical acceleration.
The use of computational techniques in predicting lift coefficients and pressure distributions of two dimenstional airfoil sections was studied. The computer code FL06/IBL was used to solve the compressible, two dimensional flow about four different airfoil sections. The lift coefficients of the airfoils were calculated at various angles of attack at subsonic Mach numbers and compared with experimental data.
Computer programs were developed to implement the computational scheme arising from Van Holten's asymptotic method for calculating airloads on a helicopter rotor blade in forward flight, and a similar technique which is based on a discretized version of the method. The basic outlines of the two programs are presented, followed by separate descriptions of the input requirements and output format. Two examples illustrating job entry with appropriate input data and corresponding output are included. Appendices contain a sample table of lift coefficient data for the NACA 0012 air foil and listings of the two programs.
A method for the calculation of the flow field of a helicopter rotor blade in hover is presented. The approach uses a finite volume solution of the three dimensional Euler equations for the blade near field. In the Euler solver the velocity field is decomposed into two parts. One is the induced velocity of the vortex wake extending below the blade, found from a free wake calculation procedure. The other part is the unknown additional velocity field of the rotor blade. This approach eliminates numerical diffusion of the rolled up wake vorticity due to truncation error and artificial viscosity. Also, the effects of the far wake are included in the limited computational domain. Solutions are presented for an isolated wing and a model helicopter rotor and compared to experiment.
Excessive vibration is one of the most prevalent technical obstacles encountered in the development of new rotorcraft. The inability to predict these vibrations is primarily due to deficiencies in analysis and simulation tools. The Langley Rotorcraft Structural Dynamics Program was instituted in 1984 to meet long term industry needs in the area of rotorcraft vibration prediction. As a part of the Langley program, this research endeavors to develop an efficient means of coupling the rotor to the airframe for preliminary design analysis of helicopter airframe vibrations. The main effort was to modify the existing computer program for modeling the dynamic and aerodynamic behavior of rotorcraft called DYSCO (DYnamic System COupler) to calculate the rotor impedance. DYSCO was recently developed for the U.S. Army and has proven to be adaptable for the inclusion of new solution methods. The solution procedure developed to use DYSCO for the calculation of rotor impedance is presented. Verification of the procedure by comparison with a known solution for a simple wind turbine model is about 75 percent completed, and initial results are encouraging. After the wind turbine impedance is confirmed, the verification effort will continue by comparison to solutions of a more sophisticated rotorcraft model. Future work includes determination of the sensitivity of the rotorcraft airframe vibrations to helicopter flight conditions and rotor modeling assumptions. When completed, this research will ascertain the feasibility and efficiency of the impedance matching method of rotor-airframe coupling for use in the analysis of airframe vibrations during the preliminary rotorcraft design process.
A finite difference code for predicting the high speed flow over the advancing helicopter rotor is presented. The code solves the low frequency, transonic small disturbance equation and is suitable for modeling the effects of advancing blade unsteadiness on blades of nearly arbitrary planform. The method employs a quasi-conservative mixed differencing scheme and solves the resulting difference equations by an alternating direction scheme. Computed results showed good agreement with experimental blade pressure data and illustrate some of the effects of varying the rotor planform. The flow unsteadiness is shown to be an indispensible part of a transonic solution. Close to the tip at high advance ratio, cross flow effects can significantly affect the solution.
A theoretical study was conducted to develop an analytical prediction method for helicopter main rotor noise due to the ingestion of atmospheric turbulence. This study incorporates an atmospheric turbulence model, a rotor mean flow contraction model and a rapid distortion turbulence model which together determine the statistics of the non-isotropic turbulence at the rotor plane. Inputs to the combined mean inflow and turbulence models are controlled by atmospheric wind characteristics and helicopter operating conditions. A generalized acoustic source model was used to predict the far field noise generated by the non-isotropic flow incident on the rotor. Absolute levels for acoustic spectra and directivity patterns were calculated for full scale helicopters, without the use of empirical or adjustable constants. Comparisons between isotropic and non-isotropic turbulence at the rotor face demonstrated pronounced differences in acoustic spectra. Turning and contraction of the flow for hover and low speed vertical ascent cases result in a 3 dB increase in the acoustic spectrum energy and a 10 dB increase in tone levels. Compared to trailing edge noise, turbulence ingestion noise is the dominant noise mechanism below approximately 30 rotor harmonics, while above 100 harmonics, trailing edge noise levels exceed turbulence ingestion noise by 25 dB.