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.
A method of mounting a cryogenically cooled, lightweight, double arch, class mirror for infrared, astronomical telescopes was developed. A 50 cm, fused silica mirror was modified for use in a new mount configuration. The flexures and the finite element analysis of the mirror stresses are reported.
Composite micromechanics are used to derive equations for predicting the elastic and thermal properties of unidirectional intraply hybrid composites. The results predicted using these equations are compared with those predicted using approximate equations based on the rule of mixtures, linear laminate theory, finite element analysis and limited experimental data. The comparisons for three different intraply hybrids indicate that all four methods predict approximately the same elastic properties and are in good agreement with measured data. The micromechanics equations and linear laminate theory predict about the same values for thermal expansion coefficients. The micromechanics equations predict through-the-thickness properties which are in good agreement with the finite element results.
The objective is to develop micromechanical models for predicting the stiffness and strength properties of textile composite materials. Two models are presented to predict the homogeneous elastic constants and coefficients of thermal expansion of a textile composite. The first model is based on rigorous finite element analysis of the textile composite unit-cell. Periodic boundary conditions are enforced between opposite faces of the unit-cell to simulate deformations accurately. The second model implements the selective averaging method (SAM), which is based on a judicious combination of stiffness and compliance averaging. For thin textile composites, both models can predict the plate stiffness coefficients and plate thermal coefficients. The finite element procedure is extended to compute the thermal residual microstresses, and to estimate the initial failure envelope for textile composites.
Nearly five decades have passed since G. Lundberg and A. Palmgren published their life theory in 1947 and 1952 and it was adopted as an ANSI/ABMA and ISO standard in 1950 and 1953. Subsequently, many variations and deviations from their life theory have been proposed, the most recent being that of E. Ioannides and T.A. Harris in 1985. This paper presents a critical analysis comparing the results of different life theories and discussing their implications in the design and analysis of rolling-element bearings. Variations in the stress-life relation and in the critical stress related to bearing life are discussed using stress fields obtained from three-dimensional, finite-element analysis of a ball in a nonconforming race under varying load. The results showed that for a ninth power stress-life exponent the Lundberg-Palmgren theory best predicts life as exhibited by most air-melted bearing steels. For a 12th power relation reflected by modern bearing steels, a Zaretsky-modified Weibull equation is superior. The assumption of a fatigue-limiting stress distorts the stress-life exponent and overpredicts life.
The RWDMES is a tool for modeling the disturbances imparted on spacecraft by spinning reaction wheels. Reaction wheels are usually the largest disturbance source on a precision pointing spacecraft, and can be the dominating source of pointing error. Accurate knowledge of the disturbance environment is critical to accurate prediction of the pointing performance. In the past, it has been difficult to extract an accurate wheel disturbance model since the forcing mechanisms are difficult to model physically, and the forcing amplitudes are filtered by the dynamics of the reaction wheel. RWDMES captures the wheel-induced disturbances using a hybrid physical/empirical model that is extracted directly from measured forcing data. The empirical models capture the tonal forces that occur at harmonics of the spin rate, and the broadband forces that arise from random effects. The empirical forcing functions are filtered by a physical model of the wheel structure that includes spin-rate-dependent moments (gyroscopic terms). The resulting hybrid model creates a highly accurate prediction of wheel-induced forces. It accounts for variation in disturbance frequency, as well as the shifts in structural amplification by the whirl modes, as the spin rate changes. This software provides a point-and-click environment for producing accurate models with minimal user effort. Where conventional approaches may take weeks to produce a model of variable quality, RWDMES can create a demonstrably high accuracy model in two hours. The software consists of a graphical user interface (GUI) that enables the user to specify all analysis parameters, to evaluate analysis results and to iteratively refine the model. Underlying algorithms automatically extract disturbance harmonics, initialize and tune harmonic models, and initialize and tune broadband noise models. The component steps are described in the RWDMES user s guide and include: converting time domain data to waterfall PSDs (power spectral densities); converting PSDs to order analysis data; extracting harmonics; initializing and simultaneously tuning a harmonic model and a wheel structural model; initializing and tuning a broadband model; and verifying the harmonic/broadband/structural model against the measurement data. Functional operation is through a MATLAB GUI that loads test data, performs the various analyses, plots evaluation data for assessment and refinement of analysis parameters, and exports the data to documentation or downstream analysis code. The harmonic models are defined as specified functions of frequency, typically speed-squared. The reaction wheel structural model is realized as mass, damping, and stiffness matrices (typically from a finite element analysis package) with the addition of a gyroscopic forcing matrix. The broadband noise model is realized as a set of speed-dependent filters. The tuning of the combined model is performed using nonlinear least squares techniques. RWDMES is implemented as a MATLAB toolbox comprising the Fit Manager for performing the model extraction, Data Manager for managing input data and output models, the Gyro Manager for modifying wheel structural models, and the Harmonic Editor for evaluating and tuning harmonic models. This software was validated using data from Goodrich E wheels, and from GSFC Lunar Reconnaissance Orbiter (LRO) wheels. The validation testing proved that RWDMES has the capability to extract accurate disturbance models from flight reaction wheels with minimal user effort.
Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.
A three-dimensional finite-element analysis of a middle-crack tension specimen subjected to mode I loading was performed to study the stress singularity along the crack front. The specimen was modeled using 20-node isoparametric elements with collapsed nonsingular elements at the crack front. The displacements and stresses from the analysis were used to estimate the power of singularities, by a log-log regression analysis, along the crack front. Analyses showed that finite-sized cracked bodies have two singular stress fields. Because of two singular stress fields near the free surface and the classical square root singularity elsewhere, the strain energy release rate appears to be an appropriate parameter all along the crack front.
A numerical and experimental investigation of the bending behavior of six eight-ply graphite-epoxy circular cylinders is presented. Bending is induced by applying a known end-rotation to each end of the cylinders, analogous to a beam in bending. The cylinders have a nominal radius of 6 inches, a length-to-radius ratio of 2 and 5, and a radius-to-thickness ratio of approximately 160. A (+/- 45/0/90)S quasi-isotropic layup and two orthotropic layups, (+/- 45/0 sub 2)S and (+/- 45/90 sub 2)S, are studied. A geometrically nonlinear special-purpose analysis, based on Donnell's nonlinear shell equations, is developed to study the prebuckling responses and gain insight into the effects of non-ideal boundary conditions and initial geometric imperfections. A geometrically nonlinear finite element analysis is utilized to compare with the prebuckling solutions of the special-purpose analysis and to study the buckling and post buckling responses of both geometrically perfect and imperfect cylinders. The imperfect cylinder geometries are represented by an analytical approximation of the measured shape imperfections. Extensive experimental data are obtained from quasi-static tests of the cylinders using a test fixture specifically designed for the present investigation. A description of the test fixture is included. The experimental data are compared to predictions for both perfect and imperfect cylinder geometries. Prebuckling results are presented in the form of displacement and strain profiles. Buckling end-rotations, moments, and strains are reported, and predicted mode shapes are presented. Observed and predicted moment vs. end-rotation relations, deflection patterns, and strain profiles are illustrated for the post buckling responses. It is found that a geometrically nonlinear boundary layer behavior characterizes the prebuckling responses. The boundary layer behavior is sensitive to laminate orthotropy, cylinder geometry, initial geometric imperfections, applied end-rotation, and non-ideal boundary conditions. Buckling end-rotations, strains, and moments are influenced by laminate orthotropy and initial geometric imperfections. Measured buckling results correlate well with predictions for the geometrically imperfect specimens. The postbuckling analyses predict equilibrium paths with a number of scallop-shaped branches that correspond to unique deflection patterns. The observed postbuckling deflection patterns and measured strain profiles show striking similarities to the predictions in some cases. Ultimate failure of the cylinders is attributed to an interlaminar shear failure mode along the nodal lines of the postbuckling deflection patterns.
This paper presents the results of a preliminary structural sizing and alternate material trade study for NASA s Crew Exploration Vehicle (CEV) Crew Module (CM). This critical CEV component will house the astronauts during ascent, docking with the International Space Station, reentry, and landing. The alternate material design study considers three materials beyond the standard metallic (aluminum alloy) design that resulted from an earlier NASA Smart Buyer Team analysis. These materials are graphite/epoxy composite laminates, discontinuously reinforced SiC/Al (DRA) composites, and a novel integrated panel material/concept known as WebCore. Using the HyperSizer (Collier Research and Development Corporation) structural sizing software and NASTRAN finite element analysis code, a comparison is made among these materials for the three composite CM concepts considered by the 2006 NASA Engineering and Safety Center Composite Crew Module project.
We have investigated laser-patterning and devitrification as a novel method to realize emergent properties in modified regions at the ribbon surface consisting of periodic localized laser spots thermally treated under rapid heating and cooling conditions. Fe-based amorphous ribbons were annealed for systematically varied laser power. Extremely high heating and cooling rates on the order of 10(exp 8)-10(exp 10)K/s have been estimated by finite element analysis simulations. Observations suggest surface melting followed by rapid solidification fast enough to quench and form an amorphous structure upon cooling. Diffusion of heat occurs from the laser irradiated surface region so that a finite volume of material surrounding this region rises in temperature above the crystallization temperature relevant for conventional isothermal annealing experiments. The underlying mechanism of laser annealing and consequences of heat transients involved are discussed in terms of impacts upon micro/nanostructure, residual stresses, and magnetic domain structure surrounding the laser irradiated region. The study illustrates the potential to exploit spatially optimized phase transformations in a scalable manufacturing process of amorphous and nanocrystalline alloys to locally access otherwise inaccessible extreme heating and cooling rates.
Design and analysis results are reported for a panel that is a 16th arc-segment of a full 33-ft diameter cylindrical barrel section of a payload fairing structure. Six such panels could be used to construct the fairing barrel, and, as such, compression buckling testing of a 16th arc-segment panel would serve as a validation test of the buckling analyses used to design the fairing panels. In this report, linear and nonlinear buckling analyses have been performed using finite element software for 16th arc-segment panels composed of aluminum honeycomb core with graphiteepoxy composite facesheets and an alternative fiber reinforced foam (FRF) composite sandwich design. The cross sections of both concepts were sized to represent realistic Space Launch Systems (SLS) Payload Fairing panels. Based on shell-based linear buckling analyses, smaller, more manageable buckling test panel dimensions were determined such that the panel would still be expected to buckle with a circumferential (as opposed to column-like) mode with significant separation between the first and second buckling modes. More detailed nonlinear buckling analyses were then conducted for honeycomb panels of various sizes using both Abaqus and ANSYS finite element codes, and for the smaller size panel, a solid-based finite element analysis was conducted. Finally, for the smaller size FRF panel, nonlinear buckling analysis was performed wherein geometric imperfections measured from an actual manufactured FRF were included. It was found that the measured imperfection did not significantly affect the panel's predicted buckling response
A parallel algorithm for the solution of the generalized eigenproblem in linear elastic finite element analysis, (K)(phi)=(M)(phi)(omega), where (K) and (M) are of order N, and (omega) is of order q is presented. The parallel algorithm is based on a completely connected parallel architecture in which each processor is allowed to communicate with all other processors. The algorithm has been successfully implemented on a tightly coupled multiple-instruction-multiple-data (MIMD) parallel processing computer, Cray X-MP. A finite element model is divided into m domains each of which is assumed to process n elements. Each domain is then assigned to a processor, or to a logical processor (task) if the number of domains exceeds the number of physical processors. The macro-tasking library routines are used in mapping each domain to a user task. Computational speed-up and efficiency are used to determine the effectiveness of the algorithm. The effect of the number of domains, the number of degrees-of-freedom located along the global fronts and the dimension of the subspace on the performance of the algorithm are investigated. For a 64-element rectangular plate, speed-ups of 1.86, 3.13, 3.18 and 3.61 are achieved on two, four, six and eight processors, respectively.
A computer program was developed for designing a low vibration gearbox. The code is based on a finite element shell analysis, a modal analysis, and a structural optimization method. In the finite element analysis, a triangular shell element with 18 degrees-of-freedom is used. In the optimization method, the overall vibration energy of the gearbox is used as the objective function and is minimized at the exciting frequency by varying the finite element thickness. Modal analysis is used to derive the sensitivity of the vibration energy with respect to the design variable. The sensitivity is representative of both eigenvalues and eigenvectors. The optimum value is computed by the gradient projection method and a unidimensional search procedure under the constraint condition of constant weight. The computer code is applied to a design problem derived from an experimental gearbox in use at the NASA Lewis Research Center. The top plate and two side plates of the gearbox are redesigned and the contribution of each surface to the total vibration is determined. Results show that optimization of the top plate alone is effective in reducing total gearbox vibration.
The shape and stability of two-dimensional finite-amplitude cellular interfaces arising during directional solidification are compared for several solidification models that account differently for latent heat released at the interface, unequal thermal conductivities of melt and solid, and solute diffusivity in the solid. Finite-element analysis and computer-implemented perturbation methods are used to analyze the families of steadily growing cellular forms that evolve from the planar state. In all models a secondary bifurcation between different families of finite-amplitude cells exists that halves the spatial wavelength of the stable interface. The quantitative location of this transition is very dependent on the details of the model. Large amounts of solute diffusion in the solid retard the growth of large-amplitude cells.
The papers presented in this volume cover a variety of topics, including sound intensity, structural intensity, modal analysis and synthesis, statistical energy analysis and energy methods, passive and active damping, and boundary element methods. Attention is also given to diagnostics and condition monitoring, material characterization and nondestructive evaluation, active noise and vibration control, sound radiation and scattering, and finite element analysis.
Metallic pressure tanks used in space missions are inherently vulnerable to hypervelocity impacts from micrometeoroids and orbital debris; thereby knowledge of impact damage and its effect on the tank integrity is crucial to a spacecraft risk assessment. This paper describes tests that have been performed to assess the effects of hypervelocity impact (HVI) damage on Titanium alloy (Ti-6Al-4V) pressure vessels burst pressure and characteristics. The tests consisted of a pair of HVI impact tests on water-filled Ti-6Al-4V tanks (water being used as a surrogate to the actual propellant) and subsequent burst tests as well as a burst test on an undamaged control tank. The tanks were placed behind Aluminum (Al) shields and then each was impacted with a 7 km/s projectile. The resulting impact debris plumes partially penetrated the Ti-6Al-4V tank surfaces resulting in a distribution of craters. During the burst tests, the tank that failed at a lower burst pressure did appear to have the failure initiating at a crater site with observed spall cracks. A fracture mechanics analysis showed that the tanks failure at the impact location may have been due to a spall crack that formed upon impact of a fragmentation on the Titanium surface. This result was corroborated with a finite element analysis from calculated Von-Mises and hoop stresses.
The flight mirror assemblies (FMA) for X-ray telescopes similar to that of the International X-ray Observatory (IXO) concept consist of several thousands of individual mirror segments. The size, shape, and location of these mirrors affect many characteristics of the telescope design. Mission requirements among other factors in turn restrict mirror segment parameters such as thickness, axial- length, azimuthal span, and mass density. This paper provides an overview of the critical relationships relating to mirror segment size and configuration throughout the design and analysis of an X-ray mirror assembly. A computational analysis is presented in the form of ray tracing pairs of thin X-ray mirror segments of varying sizes aligned in gravity and supported using kinematic constraints with corresponding self weight distortions calculated using finite element analysis (FEA). The work in this paper may be used as a starting point for determining mirror segment sizes for X-ray missions like that of IXO and beyond.
A finite element technique using overlaid plane stress elements is presented for modeling a pinned composite joint. The technique allows for discrete modeling of the pin and web regions of the joint in a two-dimensional (2-D) finite element analysis. Thus, predictive capability is substantially increased without the added complexity of a three-dimensional (3-D) analysis. The overlay technique requires a thorough evaluation of the model constraints between the joint components. Bearing panels are used to account for a reduced compressive modulus under bearing loads. Application of the plane stress overlay technique to the Filament Wound Case (FWC) Program is discussed. The experimental approach to determine bearing panel moduli using a double lap shear test is detailed, and finite element results are correlated with full-scale FWC joint test results.