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.
The evaluation of the dynamic behavior of a rotating system is possible by means of modal parameters (Eigenvalues and Eigenvectors). A mixed analytical and experimetal approach is used to identify the modal parameters of a specially designed test rig. The modal identification is done both for nonrotating as well as rotating systems. These modal parameters are used to validate a developed Finite Element Model.
This paper presents an overview of the development and qualification test campaign for the primary structure of the European Service Module of ORION, the NASA spacecraft which will serve the future human exploration missions to the Moon, Mars and beyond. Under an agreement between NASA and ESA, the ORION will be powered by a European Service Module (ESM), providing also water and oxygen for astronauts' life sustainability. The development and qualification of the European Service Module (ESM) is under ESA responsibility with Airbus Defense and Space as the prime contractor. Thales Alenia Space Italia is responsible for design development, manufacturing, assembly and qualification of the Structure subsystem. The European Service Module, installed onto the launch adapter, shall support the crew module with its adapter and a launch abort system. It shall sustain: - A combination of global and local launch loads during lift off and ascent phases, - On orbit loads induced by engine firing for orbital transfers and attitude control. The ESM structure is based on a core made of Composite Fiber Reinforced Polymer (CFRP) sandwich panels complemented by aluminum alloy platforms, longerons and secondary structures. A development campaign has been implemented in order to define and validate composite parts' strength allowable values for design: coupon tests at material level, test at component level up to breadboards tests performed on main structural components (composite to metallic joints, and at panels' discontinuities). An incremental approach as defined in [1] has been followed. A qualification static test campaign at primary structure assembly level has been implemented in order to validate the design against static stiffness and ultimate strength as well as to correlate the structural Finite Element Model (FEM) used for sizing and confirm the margins of safety. The tests have been performed successfully by Thales Alenia Space Italia (TAS-I) on two flight representative structural models (STA1, STA2), in Turin facilities (Italy) between August 2015 and March 2017, with engineering support of technical representatives from Airbus, ESA, NASA and LMCO. The main development and qualification test activities and associated results are presented and discussed in the paper
A unique single degree-of-freedom approximation technique has been developed to enable rapid application of a temporally-defined multi-spectral semi-narrow-band loading for generation of realistic stress/cycle values compared to a resonant analysis. The technique uses the harmonic analysis at resonance of a high-fidelity finite element model to produce a transfer function, which is then used to calibrate the response of the SDOF model. A standard numerical ordinary differential equation solver is then used to obtain the temporal response, and its histogram is used in a fatigue/fracture model. This technique is related to other SDOF methods used widely in industry, such as Miles' Equation and the Shock Response Spectra, but it is unique in that it produces a realistic time history of the response. The most obvious error in the process, which is the effect of closely-spaced modes, was also assessed using the parallel application of several SDOF models, and the error is shown to be small. The application of this unique and tractable reduced-order methodology has enabled the SLS program to avoid substantial cost and schedule penalties if a redesign or change of material were required. It has also enabled quick analysis of a number of other structures undergoing the same or similar excitation fields, and quick assessment when the excitation and structural configuration has been altered due to design changes in the system.
The advanced materials team investigated the use of additively manufactured metallic lattice structures for mitigating impact response in a Davis gun earth penetrator impact experiment. High-fidelity finite element models were developed and validated with quasistatic experiments. These models were then used to simulate the response of such lattices when subjected to the acceleration loads expected in the Davis gun experiment. Results reveal how the impact mitigation performance of lattices can change drastically at a certain relative density. Based on these observations, an experiment deck was designed to probe the response of lattices with different relative densities during the Davis gun phase 2 shots. The expected performance of these lattices is predicted before testing based on simulation results. The results of the Davis gun phase 2 shots are expected to provide data which will be used to assess the predictive capability of the finite element simulations in such a complex impact environment.
This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.
Ever-tightening regulations on fuel economy, and the likely future regulation of carbon emissions, demand persistent innovation in vehicle design to reduce vehicle mass. Classical methods for computational mass reduction include sizing, shape and topology optimization. One of the few remaining options for weight reduction can be found in materials engineering and material design optimization. Apart from considering different types of materials, by adding material diversity and composite materials, an appealing option in automotive design is to engineer steel alloys for the purpose of reducing plate thickness while retaining sufficient strength and ductility required for durability and safety. A project to develop computational material models for advanced high strength steel is currently being executed under the auspices of the United States Automotive Materials Partnership (USAMP) funded by the US Department of Energy. Under this program, new Third Generation Advanced High Strength Steel (i.e., 3GAHSS) are being designed, tested and integrated with the remaining design variables of a benchmark vehicle Finite Element model. The objectives of the project are to integrate atomistic, microstructural, forming and performance models to create an integrated computational materials engineering (ICME) toolkit for 3GAHSS. The mechanical properties of Advanced High Strength Steels (AHSS) are controlled by many factors, including phase composition and distribution in the overall microstructure, volume fraction, size and morphology of phase constituents as well as stability of the metastable retained austenite phase. The complex phase transformation and deformation mechanisms in these steels make the well-established traditional techniques obsolete, and a multi-scale microstructure-based modeling approach following the ICME [0]strategy was therefore chosen in this project. Multi-scale modeling as a major area of research and development is an outgrowth of the Comprehensive Test Ban Treaty of 1996 which banned surface testing of nuclear devices [1]. This had the effect that experimental work was reduced from large scale tests to multiscale experiments to provide material models with validation at different length scales. In the subsequent years industry realized that multi-scale modeling and simulation-based design were transferable to the design optimization of any structural system. Horstemeyer [1] lists a number of advantages of the use of multiscale modeling. Among these are: the reduction of product development time by alleviating costly trial-and-error iterations as well as the reduction of product costs through innovations in material, product and process designs. Multi-scale modeling can reduce the number of costly large scale experiments and can increase product quality by providing more accurate predictions. Research tends to be focussed on each particular length scale, which enhances accuracy in the long term. This paper serves as an introduction to the LS-OPT and LS-DYNA methodology for multi-scale modeling. It mainly focuses on an approach to integrate material identification using material models of different length scales. As an example, a multi-scale material identification strategy, consisting of a Crystal Plasticity (CP) material model and a homogenized State Variable (SV) model, is discussed and the parameter identification of the individual material models of different length scales is demonstrated. The paper concludes with thoughts on integrating the multi-scale methodology into the overall vehicle design.
The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms with variable time steps, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, were selected and finite element models of the configurations are described. These problems include a space shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system, and a model of the space shuttle orbiter wing. Results generally indicate a preference for implicit oer explicit algorithms for solution of transient structural heat transfer problems when the governing equations are stiff.
Small variations in gravity anomalies and topographic elevation observed in areas that have undergone highly variable amounts of upper crustal thinning can be satisfactorily explained by ductile flow of lower crustal material under the proper conditions. The boundary between the unextended Colorado Plateau and a strongly extended domain in the Basin and Range Province in the Lake Mead (Nevada) region is examined. Finite element modeling of Newtonian flow and power law creep shows that flow over the length scale of the eastern Basin and Range (500) km or more) corresponding to upper crustal extension by a factor of 1.4-3 over 10 million years requires effective viscosities less than 10 to the 18th - 10 to the 20th Pa s for ductile channels 10-25 km thick. Modeling suggests that these effective viscosities may be sustained by lower crustal material deforming at laboratory-derived power law creep rates. The longer-scale flow may require elevated crustal temperatures (more than 700 C), depending on the composition and material properties assumed. Under the boundary conditions assumed in this study the linear viscous flow models yield a satisfactory approximation to deformation by power law creep. This work suggests that flow in the lower crust may be a viable mechanism for producing small variations in total crustal thickness between strongly extended and less extended regions, and thereby explaining the relative uniformity in gravity and topography between such regions.
This paper discusses the effect of pressure on the dynamics of pre-stiffened structures such as the Advanced Solid Rocket Motor (ASRM). Previous work in which the stiffness terms resulting from constant pressure were derived has been extended to enable modeling of nonconstant pressure applied over nonenclosed volumes. These conditions will result in nonsymmetric terms in the global stiffness matrix which will not cancel out. Three new pressure stiffness elements incorporating these nonsymmetric terms have been implemented as dummy elements in COSMIC NASTRAN and have been tested on various simple examples as well as an existing ASRM NASTRAN finite element model. The results indicate that for all load cases of practical interest to the ASRM program, the nonsymmetric terms have very little effect on the dynamic characteristics. In addition, the pressure stiffness elements developed in the previous work which assumed constant pressure gave virtually the same results as the new elements even for problems in which the pressures are not constant. The original elements appear to work well as long as the pressure gradient across any individual element is no larger than about 0.75 psi/inch. The new elements are therefore most useful for determining the conditions under which the original pressure stiffness elements can be used.
The James Webb Space Telescope (JWST) Integrated Science Instrument Module (ISIM) Structure is a precision optical metering structure for the JWST science instruments. Optomechanical performance requirements place stringent limits on the allowable thermal distortion of the metering structure. A significant effort was completed to develop capabilities to predict and metrologize cryogenic thermal distortion of the ISIM Structure. This paper focuses on thermal distortion finite element modeling, analysis, and model validation. Extensive thermal distortion analysis was completed during the design phase for the ISIM Structure to demonstrate that thermal distortion requirements were achieved. Comparison of measurements from recently completed cryogenic testing and model predictions demonstrate the adequacy of thermal distortion modeling uncertainty factors adopted during the design phase, and provide bounds on the accuracy of the model predictions. This paper will provide an overview of the test configurations, describe the thermal distortion models of the tests, and provide a comparison of test results and analytical predictions from the models.
How does experimental mechanics contribute to aircraft aeroelastic airworthiness? The aviation community is always improving analysis and testing techniques to realize quicker, cheaper, reliable solutions. A ground vibration test (GVT) is conducted to identify structural mode shapes, frequencies, and damping values to validate analytical models used for flutter analysis, which shows whether a structure has acceptable aeroelastic flutter margins for airworthiness. The FLL recently conducted a GVT of the Passive Aeroelastic Tailored (PAT) Wing using an experimental modal technique called Fixed Base Correction (FBC). The GVT objective was to obtain the PAT wing modal characteristics to compare test results with finite element model (FEM) results for the tow-steered wingbox.
This paper presents an aeroservoelastic flight dynamic modeling approach for NASA Lift+Cruise eVTOL aircraft in forward flight. A conceptual aerodynamic analysis for forward flight is performed in CHARM, a rotorcraft aeromechanics unsteady vortex-lattice flow solver. A trim analysis is performed to obtain the steady state aerodynamic information at a trim point to be used for the aeroelastic analysis. A NASTRAN finite element model for the Lift+Cruise vehicle is available. The structural dynamic mode shapes are extracted from the NASTRAN for the aeroelastic analysis. An aeroservoelastic model is developed using the Theodorsen’s theory. A stability and control derivative analysis is performed to formulate a flight dynamic model of the Lift+Cruise vehicle. An aeroservoelastic flight dynamic model is formulated by coupling these two models together. Flight dynamic and flutter analyses are performed to determine vehicle stability and flutter-free operation.
Utilizing inverse uncertainty quantification techniques, structural health monitoring can be integrated with damage progression models to form probabilistic predictions of a structure's remaining useful life. However, damage evolution in realistic structures is physically complex. Accurately representing this behavior requires high-fidelity models which are typically computationally prohibitive. In the present work, a high-fidelity finite element model is represented by a surrogate model, reducing computation times. The new approach is used with damage diagnosis data to form a probabilistic prediction of remaining useful life for a test specimen under mixed-mode conditions.
In this work, a method is developed to calculate the length into a sample to which a Frequency Domain Thermoreflectance (FDTR) measurement is sensitive. Sensing depth and sensing radius are defined as limiting cases for the spherically spreading FDTR measurement. A finite element model for FDTR measurements is developed in COMSOL multiphysics and used to calculate sensing depth and sensing radius for silicon and silicon dioxide samples for a variety of frequencies and laser spot sizes. The model is compared to experimental FDTR measurements. Design recommendations for sample thickness are made for experiments where semi-infinite sample depth is desirable. For measurements using a metal transducer layer, the recommended sample thickness is three thermal penetration depths, as calculated from the lowest measurement frequency.
The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms with variable time steps, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite-element models of the configurations are described. These problems include a Space Shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system, and a model of the Space Shuttle Orbiter wing. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are 'stiff' (typical of many practical problems such as insulated metal structures).
This work investigates the response of thermoplastic composites under impact loading. Panels fabricated with the material system under consideration – carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) – were subject to impact with an Aluminum projectile, at low and high velocities. At lower velocities, the projectile rebounded from the panel after impact and caused damage e.g. delaminations. Finite element models were developed to model this test using LS-DYNA®, with material model MAT299, a continuum damage model based on the LaRC04 criteria, to capture the dynamic response and damage modes in the T700/LMPAEK panels. Interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA. Previous work by the authors showed that delaminations were consistently underpredicted in size, relative to the test. It was surmised that this was due to effects of strain rate on intralaminar and interlaminar response of the T700/LMPAEK material system. Therefore, the current work is focused on incorporating strain rate effects in modeling the interlaminar response. Rate dependence was incorporated using ‘one-way tie-break contact’ in LS-DYNA, wherein, the Mode-II fracture toughness of the LMPAEK resin, G IIc , was assumed to be an exponential function of the strain rate. In particular, G IIc , was modeled to exponentially decrease at higher strain rates. Models with this form of rate dependence favorably correlated with experimental data regarding delamination extent and rebound velocity, relative to prior models without strain rate effects.
The topics discussed in this volume include lunar endogenic rocks and processes; lunar regoliths and breccias; the terrestrial planets; shergottites; primitive materials, exposure, and atmospheres; and impacts and crater tectonics. Papers are presented on the petrology and geochemistry of alkali gabbronorites from lunar breccia 67975; the formation of Apollo 17 orange and black glass beads; mixing levels, the Apennine Front soil component, and compositional trends in the Apollo 15 soils; the meteorite component of Apollo 16 noritic impact melt breccias; and constraints on the lithospheric structure of Venus from mechanical models and tectonic surface features. Consideration is also given to a fractal interpretation of topography and geoid spectra on the earth, moon, Venus, and Mars; rare earth patterns in shergottite phosphates and residues; nuclide production by primary cosmic-ray protons; gas chromatographic instrumentation for the analysis of aerosols and gases in Titan's atmosphere; and finite-element models of non-Newtonian crater relaxation.
A lumped parameter model of a single-stage rotor system is developed which includes a new mathematical model for precision rolling element bearings to clarify this issue quantitatively and qualitatively. A comprehensive bearing stiffness matrix of dimension six is proposed which clearly demonstrates the coupling between the rotating shaft and casing motion. An eigensolution and the forced harmonic solution are obtained numerically. The natural frequencies predicted are found to be in excellent agreement with a finite element model.