SEARCH · Engineering Papers
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.
Numerical-simulation method for aeroelasticity analysis of helicopter rotor blade combines established techniques for analysis of aerodynamics and vibrations of blade. Application of method clearly shows elasticity of blade modifies flow and, consequently, aerodynamic loads on blade.
Predictive methods/computer codes for the computational simulation of acoustic fatigue resistance of hot composite structures subjected to acoustic excitation emanating from an adjacent vibrating component are discussed. Select codes developed over the past two decades at the NASA Lewis Research Center are used. The codes include computation of acoustic noise generated from a vibrating component, degradation in material properties of a composite laminate at use temperature, dynamic response of acoustically excited hot multilayered composite structure, degradation in the first ply strength of the excited structure due to acoustic loading, and acoustic fatigue resistance of the excited structure, including the propulsion environment. Effects of the laminate lay-up and environment on the acoustic fatigue life are evaluated. The results show that, by keeping the angled plies on the outer surface of the laminate, a substantial increase in the acoustic fatigue life is obtained. The effect of environment (temperature and moisture) is to relieve the residual stresses leading to an increase in the acoustic fatigue life of the excited panel.
This paper presents predictive methods/codes for computational simulation of acoustic fatigue resistance of hot composite structures subjected to acoustic excitation emanating from an adjacent vibrating component. Select codes developed over the past two decades at the NASA Lewis Research Center are used. The codes include computation of (1) acoustic noise generated from a vibrating component, (2) degradation in material properties of the composite laminate at use temperature, (3) dynamic response of acoustically excited hot multilayered composite structure, (4) degradation in the first-ply strength of the excited structure due to acoustic loading, and (5) acoustic fatigue resistance of the excited structure, including propulsion environment. Effects of the laminate lay-up and environment on the acoustic fatigue life are evaluated. The results show that, by keeping the angled plies on the outer surface of the laminate, a substantial increase in the acoustic fatigue life is obtained. The effect of environment (temperature and moisure) is to relieve the residual stresses leading to an increase in the acoustic fatigue life of the excited panel.
A summary of three-dimensional simulations on the hypersonic rarefied flows in an effort to understand the highly nonequilibrium flows about space vehicles entering the Earth's atmosphere for a realistic estimation of the aerothermal loads is presented. Calculations are performed using the direct simulation Monte Carlo method with a five-species reacting gas model, which accounts for rotational and vibrational internal energies. Results are obtained for the external flows about various bodies in the transitional flow regime. For the cases considered, convective heating, flowfield structure and overall aerodynamic coefficients are presented and comparisons are made with the available experimental data. The agreement between the calculated and measured results are very good.
Due to the unavailability and, later, prohibitive cost of the computational power required, many phenomena in nonlinear dynamic systems have in the past been addressed in terms of linear systems. Linear systems respond to periodic inputs with periodic outputs, and may be characterized in the time domain or in the frequency domain as convenient. Reduction to the frequency domain is frequently desireable to reduce the amount of computation required for solution. Nonlinear systems are only soluble in the time domain, and may exhibit a time history which is extremely sensitive to initial conditions. Such systems are termed chaotic. Dynamic buckling, aeroelasticity, fatigue analysis, control systems and electromechanical actuators are among the areas where chaotic vibrations have been observed. Direct transient analysis over a long time period presents a ready means of simulating the behavior of self-excited or externally excited nonlinear systems for a range of experimental parameters, either to characterize chaotic behavior for development of load spectra, or to define its envelope and preclude its occurrence.
The Advanced Land Observing Satellite (ALOS) was launched on January 24 2006 and has been operated successfully since then. This satellite has the attitude dynamics characterized by three large flexible structures, four large moving components, and stringent attitude/pointing stability requirements. In particular, it has one of the largest solar array paddles. Presented in this paper are flight data analyses and modeling of spacecraft attitude motion induced by the large solar array paddle. On orbit attitude dynamics was first characterized and summarized. These characteristic motions associated with the solar array paddle were identified and assessed. These motions are thermally induced motion, the pitch excitation by the paddle drive, and the role excitation. The thermally induced motion and the pitch excitation by the paddle drive were modeled and simulated to verify the mechanics of the motions. The control law updates implemented to mitigate the attitude vibrations are also reported.
Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.
A new, experimental method of measuring the material damping of test models is presented. The method involves measuring the decay of free vibrations as the model is lofted into free-fall in a vacuum. Vibration signals are transmitted through the vacuum by a miniature telemetry system. The resulting data are recorded and later analyzed using a fast Fourier transform technique to determine the percent critical damping. The experimental apparatus offers a unique means to reset the experiment without breaking the vacuum, which reduces the time between runs.
The paper presents analytical and experimental comparison of two control laws for a laboratory structure designed to simulate large space structures. The first control law is the standard linear quadratic law, which is optimal but requires model reduction for practical implementation. The second control law is a new simple direct feedback control law designed to minimize control forces while guaranteeing stability. The optimal control law was found to be only slightly better than the direct feedback law even in terms of the quadratic performance index. Moreover, the optimal control law provided almost no margin of stability for the unmodeled modes while the direct feedback law provided significant stability margins to all modes. The above results were verified experimentally using a digital implementation of the control laws. Excellent agreement between the analytical prediction and experimental measurements was observed.
A rolling joint is provided for obtaining slewing maneuvers for various apparatus including space structures, space vehicles, robotic manipulators, and simulators. Two noncircular cylinders, namely a drive and a driven cylinder, are provided in driving contact with one another. This contact is maintained by two pairs of generally S-shaped bands, each pair forming a generally 8-shaped coupling tightly about the circumferential periphery of the noncircular drive and driven cylinders. A stationarily fixed arm extends between and is rotatably journalled with a drive axle and a spindle axle respectively extending through selected rotational points of the drive cylinder and of the driven cylinder. The noncircular cylinders are profiled to obtain the desired varying gear ratio. The novelty of the present invention resides in using specifically profiled noncircular cylinders to obtain a desired varying gear ratio.
The effects of vibratory angular acceleration on detection thresholds for constant angular acceleration in a dynamic flight simulator are reported in three experiments. Detection thresholds were determined for 10 pilots and four nonpilots using a random, double-staircase procedure while the subjects sat erect in a device which rotated about an earth-vertical axis. Constant angular acceleration were presented for 0.5 and 1.0 s with concurrent, vibratory angular acceleration at 1 and 5 Hz, and thresholds with no vibratory angular acceleration were established. The thresholds were obtained while the subjects observed a visual reference in the enclosed cockpit in two experiments and in total darkness in a third. The results confirmed earlier experiments showing an inverse relationship between the duration of constant angular acceleration and detection threshold and showed that the detection thresholds in darkness were higher than with a visual reference present. Two analyses of variance revealed no significant differences in thresholds across the three vibration conditions. These results indicate that vibratory angular acceleration of fairly high levels can be present in a dynamic flight simulator without masking the pilot's ability to detect either maneuver or disturbance motions.
The quick suppression of the structural vibrations excited by bang-bang (BB) type time-optional slew maneuvers via modal-dashpot design of velocity output feedback control was investigated. Simulation studies were conducted, and modal dashpots were designed for the SCOLE flexible body dynamics. A two-stage approach was proposed for rapid slewing and precision pointing/retargeting of large, flexible space systems: (1) slew the whole system like a rigid body in a minimum time under specified limits on the control moments and forces, and (2) damp out the excited structural vibrations afterwards. This approach was found promising. High-power modal/dashpots can suppress very large vibrations, and can add a desirable amount of active damping to modeled modes. Unmodeled modes can also receive some concomitant active damping, as a benefit of spillover. Results also show that not all BB type rapid pointing maneuvers will excite large structural vibrations. When properly selected small forces (e.g., vernier thrusters) are used to complete the specified slew maneuver in the shortest time, even BB-type maneuvers will excite only small vibrations (e.g., 0.3 ft peak deflection for a 130 ft beam).
The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.
The ability of the ITD identification algorithm in identifying a complete set of structural modal parameters using a large number of free-response time histories simultaneously in one analysis, assuming a math model with a high number of degrees-of-freedom, has been studied. Identification results using simulated free responses of a uniform rectangular plate, with 225 measurement stations, and experimental responses from a ground vibration test of the Long Duration Exposure Facility (LDEF) Space Shuttle payload, with 142 measurement stations, are presented. As many as 300 degrees-of-freedom were allowed in analyzing these data. In general, the use of a significantly oversized math model in the identification process was found to maintain or increase identification accuracy and to identify modes of low response level that are not identified with smaller math model sizes. The concept of a Mode Shape Correlation Constant is introduced for use when more than one identification analysis of the same structure are conducted. This constant quantifies the degree of correlation between any two sets of complex mode shapes identified using different excitation conditions, different user-selectable algorithm constants, or overlapping sets of measurements.
The purpose of this project is to test the equivalence of inertial and passive gravitational mass in an Earth-orbiting satellite. A ground-based experiment is now well developed. It consists of comparing the motions of two cylindrical test masses suspended in precision superconducting magnetic bearings and free to move along the horizontal (axis) direction. The masses are made of niobium and lead-plated aluminum. A position detector based on a SQUID magnetometer measures the differential motion between the masses. The periods of the masses are matched by adjustment of the position detector until the system is insensitive to common mode signals, and so that the experiment is less sensitive to seismic vibration. The apparatus is contained in a twelve inch helium dewar suspended in a vibration isolation stand. The stand achieves 30 db isolation from horizontal motions between 0.1 and 60 Hz, by simulating the motion of a 200 meter long pendulum with an air bearing. With this attenuation of seismic noise and a common mode rejection ratio of 10 to the 5th power in the differential mode, the ground based apparatus should have a sensitivity to equivalence principle violations of one part in 10 to the 13th power; the satellite version might have a sensitivity of one part in 10 to the 17th power.
The ability of the ITD identification algorithm in identifying a complete set of structural modal parameters using a large number of free-response time histories simultaneously in one analysis, assuming a math model with a high number of degrees-of-freedom, has been studied. Identification results using simulated free responses of a uniform rectangular plate, with 225 measurement stations, and experimental responses from a ground vibration test of the Long Duration Exposure Facility (LDEF) Space Shuttle payload, with 142 measurement stations, are presented. As many as 300 degrees-of-freedom were allowed in analyzing these data. In general, the use of a significantly oversized math model in the identification process was found to maintain or increase identification accuracy and to identify modes of low response level that are not identified with smaller math model sizes. The concept of a Mode Shape Correlation Constant is introduced for use when more than one identification analysis of the same structure are conducted. This constant quantifies the degree of correlation between any two sets of complex mode shapes identified using different excitation conditions, different user-selectable algorithm constants, or overlapping sets of measurements.
The topic of validating structural optimization methods by use of experimental results is addressed. The need for validating the methods as a way of effecting a greater and an accelerated acceptance of formal optimization methods by practicing engineering designers is described. The range of validation strategies is defined which includes comparison of optimization results with more traditional design approaches, establishing the accuracy of analyses used, and finally experimental validation of the optimization results. Examples of the use of experimental results to validate optimization techniques are described. The examples include experimental validation of the following: optimum design of a trussed beam; combined control-structure design of a cable-supported beam simulating an actively controlled space structure; minimum weight design of a beam with frequency constraints; minimization of the vibration response of helicopter rotor blade; minimum weight design of a turbine blade disk; aeroelastic optimization of an aircraft vertical fin; airfoil shape optimization for drag minimization; optimization of the shape of a hole in a plate for stress minimization; optimization to minimize beam dynamic response; and structural optimization of a low vibration helicopter rotor.
The quality of reentry simulation for Shuttle, HERMES, Sanger, and NASP systematically suffers from the strong non-equilibrium of rotational and vibrational temperature due to the rapid acceleration of the test gas in the nozzle. Therefore the determination of temperatures is necessary and, if possible, preferable by a non-intrusive technique. The specific interests of this review are optical techniques such as electron beam fluorescence, laser-induced fluorescence, and coherent anti-Stokes Raman scattering. The capabilities available for local measurements with temporal resolution and quantitative accuracy are discussed for velocity, temperature, density, species concentrations, and fluctuations due to turbulence. The applicability of these methods of measurement is presented and discussed for the coming topic in aerothermodynamics: experimental techniques of hot gases in high enthalpy flows.