Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Element Simulation”

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.

At least 163 records · Page 9

Multibody instantly centered moving Lagrangian observer schemes. I - Formulation. II - Application to vehicular simulations

A finite element modeling method employing an instantly centered moving Lagrangian observer is proposed to solve multibody problems involving several rotating components, each with its own rotational history. The present technique allows the steady-state behavior of the problem simulation to be rendered stationary and time independent. In the second part, the method is extended to model the steady and transient response of ground-based automotive-type vehicular systems, including the modeling of vehicular obstruction rollover events.

Padovan, Joseph↗

Finite Element Analysis of Adaptive-Stiffening and Shape-Control SMA Hybrid Composites

Shape memory alloy hybrid composites with adaptive-stiffening or morphing functions are simulated using finite element analysis. The composite structure is a laminated fiber-polymer composite beam with embedded SMA ribbons at various positions with respect to the neutral axis of the beam. Adaptive stiffening or morphing is activated via selective resistance heating of the SMA ribbons or uniform thermal loads on the beam. The thermomechanical behavior of these composites was simulated in ABAQUS using user-defined SMA elements. The examples demonstrate the usefulness of the methods for the design and simulation of SMA hybrid composites. Keywords: shape memory alloys, Nitinol, ABAQUS, finite element analysis, post-buckling control, shape control, deflection control, adaptive stiffening, morphing, constitutive modeling, user element

Gao, Xiujie↗

Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms

This paper presents a final report on Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms. The focus of this research is to derive and implement: 1) Wavelet based methodologies for the compression, transmission, decoding, and visualization of three dimensional finite element geometry and simulation data in a network environment; 2) methodologies for interactive algorithm monitoring and tracking in computational mechanics; and 3) Methodologies for interactive algorithm steering for the acceleration of large scale finite element simulations. Also included in this report are appendices describing the derivation of wavelet based Particle Image Velocity algorithms and reduced order input-output models for nonlinear systems by utilizing wavelet approximations.

Kurdila, Andrew J.↗

Modeling and Simulation Credibility Assessment of Finite Element Models for Fall From Heights Injury Scenarios

Finite element (FE) computational models of the whole human body are being used to investigate the injury modalities and mechanisms arising from suited astronauts’ ExtraVehicular Activities (EVAs). These FE whole body models – simplified Elemance and THUMS pedestrian versions – have undergone a series of verification and validation (V&V) assessments for applications in automotive, domestic, and contact sports arenas. However, the FE models’ Modeling and Simulation (M&S) credibility assessments for simulating suited astronaut injuries due to fall from heights in altered gravity environments are limited. The current study evaluates the M&S credibility of simplified Elemance and THUMS pedestrian FE models based on eight credibility factors – data pedigree, input pedigree, code verification, solution verification, conceptual validation, referent validation, results uncertainty, and results robustness – as per NASA-STD-7009A [1]. These credibility factors were assigned ordinal scores ranging between 0 to 4, with a score of 0 in the case of insufficient evidence and 4 when all data for the model is known, used appropriately, and traceable to the Real-World System (RWS) with acceptable accuracy, precision, and uncertainty. The credibility factor sufficiency thresholds were assigned upon guidance from NASA EVA subject matter experts. The factor scores were ascertained through the evaluation of M&S credibility assessment procedures described in the published literature and FE model reports. The M&S credibility scores for the Elemance and THUMS FE models vary between 1 and 3 for data pedigree, input pedigree, and referent validation credibility factors, and between 0 and 1 for code verification, solution verification, conceptual validation, results uncertainty and results robustness. The results show that these FE models are credible for simulating fall from heights injury scenarios in specific kinematic and kinetic ranges that correlate with highway and contact sports-related external conditions. Multiple credibility elevation strategies are prescribed to improve the Elemance and THUMS models’ credibility for use in the EVA fall from height injury scenarios.

R A Perkins↗

Effects of Finite Element Resolution in the Simulation of Magnetospheric Particle Motion

This document describes research done in conjunction with a degree program. The purpose of the research was to compare particle trajectories in a specified set of global electric and magnetic fields; to study the effect of mesh spacing, resulting in an evaluation of adequate spacing resolution; and to study time-dependent fields in the context of substorm dipolarizations of the magnetospheric tail.

Hansen, Richard↗

An Enriched Shell Finite Element for Progressive Damage Simulation in Composite Laminates

A formulation is presented for an enriched shell nite element capable of progressive damage simulation in composite laminates. The element uses a discrete adaptive splitting approach for damage representation that allows for a straightforward model creation procedure based on an initially low delity mesh. The enriched element is veri ed for Mode I, Mode II, and mixed Mode I/II delamination simulation using numerical benchmark data. Experimental validation is performed using test data from a delamination-migration experiment. Good correlation was found between the enriched shell element model results and the numerical and experimental data sets. The work presented in this paper is meant to serve as a rst milestone in the enriched element's development with an ultimate goal of simulating three-dimensional progressive damage processes in multidirectional laminates.

McElroy, Mark W.↗

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗

FE Simulation of SMA Seal for Mars Sample Return

Several NASA rovers and landers have been on Mars and performed successful in-situ exploration. Returning Martian samples to Earth for extensive analysis is of great interest to the planetary science community. Current Mars sample return architecture would require leaving the acquired samples on Mars for years before being retrieved by subsequent mission. Each sample would be sealed securely to keep its integrity. A reliable seal technique that does not affect the integrity of the samples and uses a simple low-mass tool is required. The shape memory alloy (SMA) seal technique is a promising candidate. A study of the thermal performances of several primary designs of a SMA seal for sample tubes by finite element (FE) simulation are presented in this paper. The results show sealing the sample tube by SMA plugs and controlling the sample temperature below the allowed temperature level are feasible.

seal↗

Comparison of Nonlinear Random Response Using Equivalent Linearization and Numerical Simulation

A recently developed finite-element-based equivalent linearization approach for the analysis of random vibrations of geometrically nonlinear multiple degree-of-freedom structures is validated. The validation is based on comparisons with results from a finite element based numerical simulation analysis using a numerical integration technique in physical coordinates. In particular, results for the case of a clamped-clamped beam are considered for an extensive load range to establish the limits of validity of the equivalent linearization approach.

Rizzi, Stephen A.↗

Optical laboratory solution and error model simulation of a linear time-varying finite element equation

The use of simplified error models to accurately simulate and evaluate the performance of an optical linear-algebra processor is described. The optical architecture used to perform banded matrix-vector products is reviewed, along with a linear dynamic finite-element case study. The laboratory hardware and ac-modulation technique used are presented. The individual processor error-source models and their simulator implementation are detailed. Several significant simplifications are introduced to ease the computational requirements and complexity of the simulations. The error models are verified with a laboratory implementation of the processor, and are used to evaluate its potential performance.

Taylor, B. K.↗

Engine dynamic analysis with general nonlinear finite-element codes. I Overall approach and development of bearing damper element

NASA-sponsored research on engine dynamic simulation using general finite element nonlinear time transient computer codes available on the open market is reviewed. The approach taken was to develop software packages to model engine components which are not typically found on dynamical structures and are therefore not already computer codes. The software package developed for squeeze-film bearing dampers is outlined, and the results of a parametric study of damper pressure for a variety of specified circular orbits are presented for both long-bearing and short-bearing solutions. The data from a four-degree-of-freedom rotor-damper-stator model under conditions of small rotor unbalance through large rotor unbalance are also given.

Adams, M. L.↗

Residual Strength Analyses of Monolithic Structures

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.

Forth, Scott↗

Verification of a three-dimensional resin transfer molding process simulation model

Experimental evidence was obtained to complete the verification of the parameters needed for input to a three-dimensional finite element model simulating the resin flow and cure through an orthotropic fabric preform. The material characterizations completed include resin kinetics and viscosity models, as well as preform permeability and compaction models. The steady-state and advancing front permeability measurement methods are compared. The results indicate that both methods yield similar permeabilities for a plain weave, bi-axial fiberglass fabric. Also, a method to determine principal directions and permeabilities is discussed and results are shown for a multi-axial warp knit preform. The flow of resin through a blade-stiffened preform was modeled and experiments were completed to verify the results. The predicted inlet pressure was approximately 65% of the measured value. A parametric study was performed to explain differences in measured and predicted flow front advancement and inlet pressures. Furthermore, PR-500 epoxy resin/IM7 8HS carbon fabric flat panels were fabricated by the Resin Transfer Molding process. Tests were completed utilizing both perimeter injection and center-port injection as resin inlet boundary conditions. The mold was instrumented with FDEMS sensors, pressure transducers, and thermocouples to monitor the process conditions. Results include a comparison of predicted and measured inlet pressures and flow front position. For the perimeter injection case, the measured inlet pressure and flow front results compared well to the predicted results. The results of the center-port injection case showed that the predicted inlet pressure was approximately 50% of the measured inlet pressure. Also, measured flow front position data did not agree well with the predicted results. Possible reasons for error include fiber deformation at the resin inlet and a lag in FDEMS sensor wet-out due to low mold pressures.

Fingerson, John C.↗

Residual Strength Characterization of a Curved Integrally-Stiffened Panel

Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. 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 intension 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 curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. 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 (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).

Seshadri, B. R.↗

Material Variability Effects on Damage Development Within Composite Adhesively Bonded Joints

Use of composite adhesively bonded joints (ABJ) is of critical importance to the adoption of composite materials in the automotive industry, as ABJ enable lower stress concentrations as compared to conventional mechanically fastened joints. ABJ are better suited for joining composite materials as compared to fastened joints because fastened joints require drilling of holes which locally affect composite material structure. Composite materials and adhesives are subject to unavoidable stochastic local material variations which make different failure scenarios possible. An experimentally tested ABJ configuration is simulated using finite element analysis (FEA). Experimentally, under tension, the joints failed by three major failure modes, with peak loads ranging from 13.0-16.1 kips. Progressive failure analysis tools are used to simulate damage development within each material within the joint. The simulation agreed well with the average experimental peak load. Stochastically occurring adhesive porosity and matrix-fiber micro-disbonding were numerically simulated. The simulations revealed a similar trend as observed experimentally: joints which failed at higher peak loads had lower levels of damage within the face-sheets of the composite panels which were adhesively bonded; these joints which failed at higher peak loads also had greater damage in the doubler of the experimentally tested double lap joint configuration.

Richard Larson↗

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↗

Comparisons of Kinematics and Dynamics Simulation Software Tools

Kinematic and dynamic analyses for moving bodies are essential to system engineers and designers in the process of design and validations. 3D visualization and motion simulation plus finite element analysis (FEA) give engineers a better way to present ideas and results. Marshall Space Flight Center (MSFC) system engineering researchers are currently using IGRIP from DELMIA Inc. as a kinematic simulation tool for discrete bodies motion simulations. Although IGRIP is an excellent tool for kinematic simulation with some dynamic analysis capabilities in robotic control, explorations of other alternatives with more powerful dynamic analysis and FEA capabilities are necessary. Kinematics analysis will only examine the displacement, velocity, and acceleration of the mechanism without considering effects from masses of components. With dynamic analysis and FEA, effects such as the forces or torques at the joint due to mass and inertia of components can be identified. With keen market competition, ALGOR Mechanical Event Simulation (MES), MSC visualNastran 4D, Unigraphics Motion+, and Pro/MECHANICA were chosen for explorations. In this study, comparisons between software tools were presented in terms of following categories: graphical user interface (GUI), import capability, tutorial availability, ease of use, kinematic simulation capability, dynamic simulation capability, FEA capability, graphical output, technical support, and cost. Propulsion Test Article (PTA) with Fastrac engine model exported from IGRIP and an office chair mechanism were used as examples for simulations.

Yeu-Sheng Paul Shiue↗