SEARCH · Engineering Papers

Results for “Finite element model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Study for the optimization of a transport aircraft wing for maximum fuel efficiency. Volume 1: Methodology, criteria, aeroelastic model definition and results

Work performed in the design of a transport aircraft wing for maximum fuel efficiency is documented with emphasis on design criteria, design methodology, and three design configurations. The design database includes complete finite element model description, sizing data, geometry data, loads data, and inertial data. A design process which satisfies the economics and practical aspects of a real design is illustrated. The cooperative study relationship between the contractor and NASA during the course of the contract is also discussed.

Radovcich, N. A.↗

NASA NTRS1998

Numerical Modeling of Solidification in Space With MEPHISTO-4

A pre-flight analysis of the directional solidification of BiSn with MEPHISTO-4 is presented. Simplified Bridgman growth under microgravity conditions is simulated using a two dimensional finite element model. This numerical model is a single domain, pseudo-steady state model, and includes the effects of both thermal and solutal convection. The results show that for all orientations of the applied steady state gravity vector, of magnitude 1 micro-g, the directional solidification process remains diffusion controlled. The maximum convective velocity was found to be 4.424 x 10(exp -5) cm/s for the horizontal Bridgman growth configuration. This value is an order of magnitude lower than the growth velocity. The maximum and minimum values or solute concentration in the liquid at the crystal-melt interface were 13.867 at.% and 13.722 at.%, respectively. This gives a radial segregation value of xi = 1.046% at the interface. A secondary objective of this work was to compare the results obtained to those that consider thermal convection only (no solutal convection). It was found that the convective flow patterns in simulations which included solutal convection were significantly different from those which ignored solutal convection. The level of radial segregation predicted by the current simulations is an order of magnitude lower than that found in simulations which ignore solutal convection. The final aim was to investigate the effect of g-jitter on the crystal growth process. A simulation was performed to calculate the system response to a 1 second, 100 micro-g gravity impulse acting normal to the direction of growth. This pulse is consistent with that induced by Orbiter thruster firings. The results obtained indicate that such a gravity pulse causes an increase in the level of radial solute segregation at the interface from the steady state values. The maximum value of solute concentration in the liquid was found to be 13.888 at.%, the minimum value calculated was 13.706 at.%, yielding a radial segregation value of xi = 1.31% at the interface. These values occurred 126 seconds after the pulse terminated. Thus it is anticipated that the process will remain diffusion controlled even when subjected to such g-jitter.

Simpson, James E.↗

NASA NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗

DOE OSTI2021

A Combined Experimental and Modeling Approach to Investigate the Performance of Joint Between AZ31 Magnesium and Uncoated DP590 Steel Using Friction Stir-Assisted Scribe Technique

In this study, friction stir assisted scribe technique is used to achieve a viable joint between AZ31 magnesium alloy and uncoated DP590 steel. Here, the mechanical properties and performance of the joint are investigated using finite element modeling and electron microscopy. The joint strength can be ascribed to two factors, the mechanical interlocking and the metallurgical bonding of two materials at their interface. The contribution of both factors is intended, and the metallurgical bonding is studied to understand the basic joining phenomenon.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Simulation of Creep Deformation and Failure in Graded AM Microstructures

This report describes modeling tools and techniques developed to simulate the long-term material performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF). A physics-based Crystal Plasticity Finite Element model is used to simulate creep in microstructures and to study the roles of grain morphology, porosity, and texture. We describe our modeling methodology, including an orientation-mapping technique to capture the spatially varying crystallographic orientation that results from the build conditions. Our study of microstructural features shows that AM microstructures produced by LPBF tend to creep faster in the build direction, while texture and grain boundaries strengthen the transverse directions. However, when grain-boundary porosity and the consequent cavity growth are included in the model, the transverse directions begin to creep faster. In examining texture, the results indicate that spatially varying orientation arising from the build conditions increases anisotropy in the material, making it critical to account for orientation gradients in the material to accurately model its mechanical behavior. We also describe a material-model calibration campaign in which we calibrated the constitutive model specifically for LPBF 316H stainless steel at 725℃ for both solution-annealed and as-built conditions. Finally, these tools and techniques are used to model creep in microstructures representing different regions of an LPBF material with graded microstructure, owing to intentional variation in processing conditions. The creep simulation results show good agreement with experimental data across all three microstructures, with future work planned to study rupture in the material.

36 MATERIALS SCIENCE↗

NASA NTRS1981

SPAR analysis of LDEF vibration characteristics

This paper presents a structural dynamic modeling of the Long Duration Exposure Facility (LDEF), which is a Space Shuttle payload of passive scientific experiments contained in trays mounted on a large cylindrically shaped structure. Special detailed finite element modeling, using the SPAR system of computer programs was required to obtain good agreement between analytical and test vibrations modes. Experimental trays contributed significantly to overall LDEF stiffnesses, and these contributions were realistically represented for each tray by the stiffness matrix of an equivalent orthotropic panel in the overall LDEF SPAR model. Orthotropic stiffnesses for this panel were obtained from finely detailed statically loaded tray SPAR models in which stiffness coupling was accounted for along with partial relative sliding allowed by the tray clamping attachments. Joint boundary conditions were also significant in the structural dynamic modeling of LDEF, and static data proved valuable in assessing modeling of local end fittings.

Edighoffer, H. H.↗

NASA NTRS1988

Improved method for stress and compatibility analysis of multicomponent rotating systems

An improved method of analyzing multicomponent rotating assemblies for the determination of operating stresses and component compatibility has been developed. In this method, a single finite element model is developed which contains all of the separate components in the rotating assembly. This is made possible by using gap elements to simulate the contact surfaces between components. The MARC finite element computer program is then used to perform the analysis. This improved method is less time consuming and more reliable than the conventional method of analyzing such systems. Results are presented for two different stress-compatibility analyses of a six-component axial flow compressor rotor. The results for the previously used flexibility analysis method are compared with those for the improved analysis method. The stresses predicted by each method compare quite well with each other. The predictions of the component compatibility, as well as the magnitude of the forces at the contact surfaces, also compare well for these two analysis procedures. It is therefore recommended that the improved analysis method be used to determine the stress-compatibility characterstics of multicomponent rotating systems.

Carek, Gerald A.↗

NASA NTRS2004

Numerical and Experimental Dynamic Characteristics of Thin-Film Membranes

Presented is a total-Lagrangian displacement-based non-linear finite-element model of thin-film membranes for static and dynamic large-displacement analyses. The membrane theory fully accounts for geometric non-linearities. Fully non-linear static analysis followed by linear modal analysis is performed for an inflated circular cylindrical Kapton membrane tube under different pressures, and for a rectangular membrane under different tension loads at four comers. Finite element results show that shell modes dominate the dynamics of the inflated tube when the inflation pressure is low, and that vibration modes localized along four edges dominate the dynamics of the rectangular membrane. Numerical dynamic characteristics of the two membrane structures were experimentally verified using a Polytec PI PSV-200 scanning laser vibrometer and an EAGLE-500 8-camera motion analysis system.

Young, Leyland G.↗

NASA NTRS1991

A rheologically layered three-dimensional model of the San Andreas fault in central and southern California

The effects of rheological parameters and the fault slip distribution on the horizontal and vertical deformation in the vicinity of the fault are investigated using 3D kinematic finite element models of the San Andreas fault in central and southern California. It is shown that fault models with different rheological stratification schemes and slip distributions predict characteristic deformation patterns. Models that do not include aseismic slip below the fault locking depth predict deformation patterns that are strongly dependent on time since the last earthquake, while models that incorporate the aseismic slip below the locking depth depend on time to a significantly lesser degree.

Williams, Charles A.↗

DOE OSTI2022

Computational Modeling of Multi-Pass Rolling Parameters Effect on Resulting Fuel Foil Shape

A focus of the U.S. Department of Energy is to improve production yield and reduce the cost of Low Enriched Uranium (LEU)-molybdenum alloy (U-10Mo) monolithic fuel plates that will be replacing High Enriched Uranium (HEU) oxide dispersion fuels used currently in the United States High Performance Research Reactors (USHPRR). One area of improvement is lowering the high transverse waviness and longitudinal waviness currently present within rolled foils prior to cladding to produce fuel plates. Traditional rolling manufacturing techniques for other metal foils use winders to pull and straighten the foil as it is rolled back and forth to the final thickness. This approach cannot be used to roll thin U-10Mo foils (0.008-0.025” thick) because only small castings can be rolled due to nuclear criticality safety concerns. As a result, the fuel foils are too short (1 m in length) to use traditional winders. Therefore, it is crucial to identify other rolling parameters (i.e., roller friction, axial tension load, roller diameter, and roll pass reduction percent) that might reduce transverse waviness and longitudinal waviness in the rolled fuel foil and develop a high-yield, low-cost multi-pass rolling manufacturing process. This report documents a systematic finite element modeling study to investigate the effects of numerous rolling parameters to reduce resulting transverse waviness and longitudinal waviness in the fuel foil during multi-pass rolling of U-10Mo foils. The rolling of a U-10Mo plate with initial dimensions of 1”x1”x 0.048” is modeled using Abaqus CAE. This rolling is modeled to undergo eight 20% reduction roll passes to a final fuel foil thickness of 0.01”. The elastic-plastic constitutive model of the U-10Mo alloy was input to the fuel foil rolling model. The rollers were modeled as rigid bodies. A comparison of rolling friction coefficients of 0.3 and 0.7 over a wide range of applied axial tension loads were investigated in order to evaluate the effect of using a lubricant during rolling. The effect of roller diameter on the resulting transverse waviness and longitudinal waviness of the fuel foil over a wide range of axial tension loads were also investigated by modeling rollers 7/8” and 3.75” in diameter. The results of this systematic finite element method study will aid manufacturers in producing low transverse waviness and reduced longitudinal waviness in U-10Mo fuel foils.

U-10Mo, FEA, Rolling, Residual Stress, Fuel Foil↗

NASA NTRS1980

Design, fabrication and test of graphite/polyimide composite joints and attachments for advanced aerospace vehicles

Bonded and bolted designs are presented for each of four major attachment types. Prepreg processing problems are discussed and quality control data are given for lots 2W4604, 2W4632 and 2W4643. Preliminary design allowables test results for tension tests and compression tests of laminates are included. The final small specimen test matrix is defined and the configuration of symmetric step-lap joint specimens are shown. Finite element modeling studies of a double lap joint were performed to evaluate the number of elements required through the adhesive thickness to assess effects of various joint parameters on stress distributions. Results of finite element analyses assessing the effect of an adhesive fillet on the stress distribution in a double lap joint are examined.

Skoumal, D. E.↗

NASA NTRS1987

Design modeling for shape optimization

Some important aspects of design modeling for shape optimization are discussed for both stamped sheet metal components and cast solid components. For stamped components the basis for the modeling approach is a boundary design function. Design parameters control the shape of 2-D regions. For more complex, folded plate components, the 2-D regions can be assembled using translation and rotation operations. The analysis model is automatically created using a mesh generation procedure requiring only boundary data. For less complex solid components, it was found that this approach is not suitable. For these structures, the finite element models are typically created using very sophisticated graphical modeling systems. A new approach which overlays a parameterized surface design model on an existing analysis model is described. To summarize, the future needs for solid shape design is described in terms of an extension of the previously described 2-D capability.

Botkin, M. E.↗

NASA NTRS1996

Crack Development in Cross-Ply Laminates Under Uniaxial Tension

This study addresses matrix-dominated failures in carbon fiber/polymer matrix composite laminates in a cross-ply lay-up. The events of interest are interlaminar fracture in the form of transverse cracks in the 90' plies and longitudinal splitting in the 0 deg plies and interlaminar fracture in the form of 0/90 delamination. These events were observed using various nondestructive evaluation (NDE) techniques during static tensile tests. Acoustic emission (AE) x radiography, and edge view microscopy were the principal ones utilized in a real-time environment. A comparison of the NDE results with an analytical model based on the classical linear fracture mechanics concept of strain energy release rate as a criterion for crack growth was performed. The virtual crack closure theory was incorporated with a finite element model to generate strain energy release rate curves for the analytical case. Celion carbon fiber/polyimide matrix (G30-500/PMR-15) was the material tested with cross-ply lay-ups of (0(2)/90(6))s and (0(4)/90(4))s. The test specimens contained thermally induced cracks caused by the high-temperature processing. The analytical model was updated to compensate for the initial damage and to study further accumulation by taking into account the crack interactions. By correlating the experimental and analytical data, the critical energy release rates were found for the observable events of interest.

Gyekenyesi, Andrew L.↗

NASA NTRS2008

High-Temperature Strain Sensing for Aerospace Applications

Thermal protection systems (TPS) and hot structures are utilizing advanced materials that operate at temperatures that exceed abilities to measure structural performance. Robust strain sensors that operate accurately and reliably beyond 1800 F are needed but do not exist. These shortcomings hinder the ability to validate analysis and modeling techniques and hinders the ability to optimize structural designs. This presentation examines high-temperature strain sensing for aerospace applications and, more specifically, seeks to provide strain data for validating finite element models and thermal-structural analyses. Efforts have been made to develop sensor attachment techniques for relevant structural materials at the small test specimen level and to perform laboratory tests to characterize sensor and generate corrections to apply to indicated strains. Areas highlighted in this presentation include sensors, sensor attachment techniques, laboratory evaluation/characterization of strain measurement, and sensor use in large-scale structures.

Piazza, Anthony↗

NASA NTRS2018

Automation of Structural Sizing of Aircraft Concepts Under Static Aeroelastic Constraints

This paper presents an automation process for structural sizing of subsonic and supersonic aircraft concepts under static aeroelastic constraints. The automation process starts with an OpenVSP geometry and ends with a PATRAN plot of a NASTRAN solution for static aeroelastic analysis or optimization. ModelCenter is used to integrate all analysis codes with easy-to-use interfaces. Automation tools are developed to streamline the setup process and avoid user errors. Fuel is distributed by solving an optimization problem to match the center of gravity of aircraft at a specified flight condition. Fuel weights are also automatically attached to the structural model as point masses. All other weights used in FLOPS mission analysis (excluding fuselage and wing structural weights) are automatically attached to or smeared on the structural model. For any given OpenVSP geometry and FLOPS analysis data, a static aeroelastic sizing model for NASTRAN analysis can be generated in a couple of hours. The empirical fuselage and wing structural weights from FLOPS are replaced by structural panel weights from the sized finite-element model. Three supersonic and two subsonic aircraft concepts are used to demonstrate the automation process as a physics-based weight estimation tool for aircraft conceptual design.

Li, Wu↗

NASA NTRSDate not supplied

SLS Integrated Modal Test Uncertainty Quantification using the Hybrid Parametric Variation Method

Uncertainty in structural loading during launch is a significant concern in the development of spacecraft and launch vehicles. Small variations in launch vehicle and payload mode shapes and their interaction can result in significant variation in system loads. In many cases involving large aerospace systems it is difficult, not economical, or impossible to perform a system modal test. However, it is still vital to obtain test results that can be compared with analytical predictions to validate models. Instead, the “Building Block Approach” is used in which system components are tested individually. Component models are correlated and updated to agree as best they can with test results. The Space Launch System consists of a number of components that are assembled into a launch vehicle. Finite element models of the components are developed, reduced to Hurty/Craig-Bampton models and assembled to represent different phases of flight. The only opportunity to obtain modal test data from an assembled Space Launch System will be during the Integrated Modal Test. There is always uncertainty in every model, which flows into uncertainty in predicted system results. Uncertainty Quantification is used to determine statistical bounds on prediction accuracy based on model uncertainty. For the Space Launch System, model uncertainty is at the Hurty/Craig-Bampton component level. Uncertainty in the Hurty/Craig-Bampton components is quantified using the hybrid parametric variation approach that combines parametric and nonparametric uncertainty. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. This type of uncertainty cannot be represented by variations infinite element model input parameters and thus cannot be included in a parametric approach. However, model-form uncertainty can be modeled using a nonparametric approach based on random matrix theory. The hybrid parametric variation method requires the selection of dispersion values for the Hurty/Craig-Bampton fixed-interface eigenvalues, and the Hurty/Craig-Bampton stiffness matrices. Component test/analysis frequency error is used to identify the fixed-interface eigenvalue dispersions, while test/analysis cross-orthogonality is used to identify stiffness dispersion values. The hybrid parametric variation uncertainty quantification approach is applied to the Space Launch System Integrated Modal Test configuration. Monte Carlo analysis is performed, and statistics are determined for modal correlation metrics, frequency response from Integrated Modal Test shakers to selected accelerometers, as well as other metrics for determining how well target modes are excited and identified. If the predicted uncertainty envelopes future Integrated Modal Test results, then there will be increased confidence in the utility of the component-based hybrid parametric variation uncertainty quantification approach.

Uncertainty Quantification↗

NASA NTRS2010

Modeling Geometry and Progressive Failure of Material Interfaces in Plain Weave Composites

A procedure combining a geometrically nonlinear, explicit-dynamics contact analysis, computer aided design techniques, and elasticity-based mesh adjustment is proposed to efficiently generate realistic finite element models for meso-mechanical analysis of progressive failure in textile composites. In the procedure, the geometry of fiber tows is obtained by imposing a fictitious expansion on the tows. Meshes resulting from the procedure are conformal with the computed tow-tow and tow-matrix interfaces but are incongruent at the interfaces. The mesh interfaces are treated as cohesive contact surfaces not only to resolve the incongruence but also to simulate progressive failure. The method is employed to simulate debonding at the material interfaces in a ceramic-matrix plain weave composite with matrix porosity and in a polymeric matrix plain weave composite without matrix porosity, both subject to uniaxial cyclic loading. The numerical results indicate progression of the interfacial damage during every loading and reverse loading event in a constant strain amplitude cyclic process. However, the composites show different patterns of damage advancement.

Hsu, Su-Yuen↗

NASA NTRS1976

A finite element-analytical method for modeling a structure in an infinite fluid

A method is described from which the interaction of an elastic structure with an infinite acoustic fluid is determined. The displacements of the structure and the pressure field of the immediate surrounding fluid are modeled by finite elements, and the remaining pressure field of the infinite fluid region is given by an analytical expression. This method yields a frequency dependent boundary condition for the outer fluid boundary when applied to the frequency response of an elastic beam in contact with an acoustic fluid. The frequency response of the beam is determined using NASTRAN, and compares favorably to the exact solution which is also presented. The effect of the fluid on the response of the structure at low and high frequencies is due to added mass and damping characteristics, respectively.

Zarda, P. R.↗