Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element (FE) 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 37 records · Page 2

Development and Analysis of Energy Absorbing Subfloor Concepts to Improve eVTOL Crashworthiness

To help ensure safe transportation within electric Vertical Take-off and Landing (eVTOL) vehicles the National Aeronautics and Space Administration (NASA) has been developing novel energy absorbing (EA) design technologies to improve crashworthiness within the unique design constraints of these vehicles. As part of this effort, a series of lightweight energy absorbing subfloor concepts were developed for potential use within eVTOL vehicle design. The capability of the subfloor designs was first evaluated through finite element (FE) model simulation in both component and vehicle level impact conditions. Knowledge gained from these analyses were used to iterate upon the design prior to fabrication. Fabrication and testing of the subfloor designs has begun and will be used to verify predicted capability. Results from FE model analysis was used to down select to a final subfloor geometry for additional component level optimization and full-scale test validation.

eVTOL↗

Development of Energy Absorbing Subfloor Concepts to Improve Crashworthiness of eVTOL Aircraft

To help ensure safe transportation within electric Vertical Take-off and Landing (eVTOL) vehicles the National Aeronautics and Space Administration (NASA) has been developing novel energy absorbing (EA) design technologies to improve crashworthiness within the unique design constraints of these vehicles. As part of this effort, a series of lightweight energy absorbing subfloor concepts were developed for potential use within eVTOL vehicle design. The capability of the subfloor designs was first evaluated through finite element (FE) model simulation in both component and vehicle level impact conditions. Knowledge gained from these analyses were used to iterate upon the design prior to fabrication. Fabrication and testing of the subfloor designs has begun and will be used to verify predicted capability. Results from FE model analysis was used to down select to a final subfloor geometry for additional component level optimization and full-scale test validation.

Crashworthiness↗

Multi-scale characterization and simulation of impact welding between immiscible Mg/steel alloys

Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590, which are typically considered as un-weldable due to their high physical property disparities, low mutual solubility, and the lack of any intermetallic phases. Characterization results from scanning electron microscopy (SEM) and high-resolution transmission electron microscopy (HRTEM) of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles. The interlayer exhibits intact bonding with both DP590 and AZ31 substrates. To investigate the fundamental bond formation mechanisms at the interface, a finite element (FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition. Finally, taking the FE results at the interface as inputs, molecular dynamics (MD) simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling. The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer, leading to the metallurgical bond between Mg/steel alloys.

36 MATERIALS SCIENCE↗

Stress Relaxation Cracking of Alloys at Temperatures Higher Than 540°C

Type 347H stainless steel (347H SS), used in commercial concentrating solar power (CSP) thermal energy storage to store solar-salt at a temperature of 565°C, has been reported in the literature to be susceptible to stress-relaxation cracking (SRC). The welded heat-affected zone (HAZ) and fusion zone (FZ) of 347H SS, particularly in thick sections, are known to be susceptible to failure during post-weld heat treatment (reheat cracking). SRC could also occur after months or years under an elevated-temperature service environment. Two conditions must be present for failure to occur in the HAZ and/or FZ: 1) a modified or sensitized microstructure and 2) sufficiently high tensile residual stresses present at the elevated service temperature. The overarching goal of this project is to recommend SRC mitigation protocols to avoid susceptibility to fail through SRC at temperatures relevant for CSP. Post weld heat treatment conditions and alternative alloys are investigated as potential mitigation solutions to SRC. We used Gleeble thermomechanical simulation tests and finite element (FE) models to understand the susceptibility of 347H SS to SRC as a function of temperature, stress, and microstructure. We started with a literature review of weldability issues with 347H SS and techniques to mitigate SRC in 347H SS weldments. Next, we used Gleeble experiments to determine reheat cracking susceptibility in the simulated HAZ of 347SS weldments and an alternative alloy, 316L SS (NUCL 167 SPH) with boron added. We also performed Gleeble experiments to compare the reheat cracking susceptibility of 347H cross welded with two different fillers: E347, which is used in some commercial CSP TES tanks, and E16.8.2. We validated the experimental results with FE models of the residual stresses. We found that although the 316L (NUCL 167 SPH) is less susceptible to reheat cracking, it is slightly weaker than 347H and ASME BP&V codes limit its use to a service condition of 565°C. We also found that welds using E16.8.2 as the weld filler are less susceptible to failure than those using E347, likely due to the higher creep ductility of E16.8.2, and that it may be used as an alternative filler for repair welding of 347H welds or as the primary choice of filler for newly developed weld joints. We also found that post weld heat treatment could be a viable solution for mitigating stress in E347-347H SS thick, constrained welds, like those found in CSP tanks, and propose several options for mitigating SRC in existing and future TES tanks.

14 SOLAR ENERGY↗

The Impact of Ocular Pressures, Material Properties and Geometry on Optic Nerve Head Deformation

Alteration in intracranial pressure (ICP) has been associated with various diseases that cause visual impairment, including glaucoma, idiopathic intracranial hypertension and Visual Impairment and Intracranial Pressure (VIIP) syndrome. However, how changes in ICP lead to vision loss is unclear, although it is hypothesized to involve deformations of the tissues in the optic nerve head (ONH). Recently, understanding the effect of ICP alterations on ocular tissues has become a major concern for NASA, where 42 of astronauts that partake in long duration space missions suffer from VIIP syndrome. Astronauts with VIIP syndrome suffer from visual impairment and changes in ocular anatomy that persist after returning to earth (1). It is hypothesized that the cephalad fluid shift that occurs upon entering microgravity increases ICP, which leads to an altered biomechanical environment in the posterior globe and optic nerve sheath, and subsequently VIIP syndrome. Our goal was to develop a finite element (FE) model to simulate the acute effects of elevated ICP on the posterior eye. Here, we simulated how inter-individual differences affect the deformation of ONH tissues. Further, we examined how several different geometries influenced deformations when exposed to elevated ICP.

mathematical models↗

High Strain Rate Response of Adhesively Bonded Fiber-Reinforced Composite Joints – A Computational Study to Guide Experimental Design

Adhesively bonded carbon fiber-reinforced epoxy composite laminates are widely used in aerospace applications. During a high energy impact event, these laminates are often subjected to high strain rate loading. However, the influence of high strain rate loading on the response of these composite joints is not well understood. Computational finite element (FE) modeling and simulations are conducted to guide the design of high strain rate experiments. Two different experimental designs based on split Hopkinson bar were numerically modeled to simulate Mode I and Mode II types loading in the composite. In addition, the computational approach adopted in this study helps in understanding the high strain rate response of adhesively bonded composite joints subjected to nominally Mode I and Mode II loading. The modeling approach consists of a ply-level 3D FE model, a progressive damage constitutive model for the composite material behavior and a cohesive tie-break contact element for interlaminar delamination.

Ravindran, Suraj↗

Analysis of a Landing System for Planetary Payloads Utilizing Passive Energy Absorbing Composite Structure

Delivery of a payload from space to a planetary surface currently requires the development of an application specific landing system to protect the payload from forces imparted during impact with the planet surface. Often, active energy attenuating systems such as retro-rockets, deployable parachutes, and airbags are utilized within these landing systems to reduce landing impact energy. Unfortunately, these active systems come at a cost; active energy attenuating systems are susceptible to system faults which may limit or completely negate their energy attenuating capability. Additionally, components needing to be stowed such as fuel, parachutes, and airbags increase design complexity, cost, and weight. To overcome these limitations, this study examines the potential of passive energy attenuation through energy absorbing structural design and composite materials to mitigate landing loads for small payload planetary delivery. Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted extensive research into developing energy absorbing structures and components for the attenuation of impact energy under various loading conditions including aircraft crash and spacecraft impact. The current study leverages this research to design a lightweight planetary delivery system which utilizes unique outer mold line (OML) geometry and passive energy absorbing structural design to limit landing loads across potential planetary surface environments. The OML geometry is designed to control impact orientation and provide self-righting capabilities for slopped impact surfaces. The internal structure is composed of composite material structures arranged to provide energy absorption which is robust to impact angle and impact velocity. The developed planetary delivery design concept will be evaluated using finite element (FE) model analysis. Simulations of landing impacts with representative soil surface environments will be used to characterize the energy absorbing capabilities of the landing system. Sensitivity of predicted impact force to landing environment, impact angle, and impact velocity will be assessed to identify capabilities and limitations of the initial structural design. Results will be used to determine the feasibility of a lightweight composite structure to passively absorb landing energy for robust planetary payload delivery.

Crashworthiness↗

Orion Crew Member Injury Predictions during Land and Water Landings

A review of astronaut whole body impact tolerance is discussed for land or water landings of the next generation manned space capsule named Orion. LS-DYNA simulations of Orion capsule landings are performed to produce a low, moderate, and high probability of injury. The paper evaluates finite element (FE) seat and occupant simulations for assessing injury risk for the Orion crew and compares these simulations to whole body injury models commonly referred to as the Brinkley criteria. The FE seat and crash dummy models allow for varying the occupant restraint systems, cushion materials, side constraints, flailing of limbs, and detailed seat/occupant interactions to minimize landing injuries to the crew. The FE crash test dummies used in conjunction with the Brinkley criteria provides a useful set of tools for predicting potential crew injuries during vehicle landings.

Lawrence, Charles↗

The Use of a Vehicle Acceleration Exposure Limit Model and a Finite Element Crash Test Dummy Model to Evaluate the Risk of Injuries During Orion Crew Module Landings

A review of astronaut whole body impact tolerance is discussed for land or water landings of the next generation manned space capsule named Orion. LS-DYNA simulations of Orion capsule landings are performed to produce a low, moderate, and high probability of injury. The paper evaluates finite element (FE) seat and occupant simulations for assessing injury risk for the Orion crew and compares these simulations to whole body injury models commonly referred to as the Brinkley criteria. The FE seat and crash dummy models allow for varying the occupant restraint systems, cushion materials, side constraints, flailing of limbs, and detailed seat/occupant interactions to minimize landing injuries to the crew. The FE crash test dummies used in conjunction with the Brinkley criteria provides a useful set of tools for predicting potential crew injuries during vehicle landings.

Lawrence, Charles↗

Orion Crew Member Injury Predictions during Land and Water Landings

A review of astronaut whole body impact tolerance is discussed for land or water landings of the next generation manned space capsule named Orion. LS-DYNA simulations of Orion capsule landings are performed to produce a low, moderate, and high probability of injury. The paper evaluates finite element (FE) seat and occupant simulations for assessing injury risk for the Orion crew and compares these simulations to whole body injury models commonly referred to as the Brinkley criteria. The FE seat and crash dummy models allow for varying the occupant restraint systems, cushion materials, side constraints, flailing of limbs, and detailed seat/occupant interactions to minimize landing injuries to the crew. The FE crash test dummies used in conjunction with the Brinkley criteria provides a useful set of tools for predicting potential crew injuries during vehicle landings.

Lawrence, Charles↗

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗

Efficient finite element simulation of slot spirals, slot radomes and microwave structures

This progress report contains the following two documents: (1) 'Efficient Finite Element Simulation of Slot Antennas using Prismatic Elements' - A hybrid finite element-boundary integral (FE-BI) simulation technique is discussed to treat narrow slot antennas etched on a planar platform. Specifically, the prismatic elements are used to reduce the redundant sampling rates and ease the mesh generation process. Numerical results for an antenna slot and frequency selective surfaces are presented to demonstrate the validity and capability of the technique; and (2) 'Application and Design Guidelines of the PML Absorber for Finite Element Simulations of Microwave Packages' - The recently introduced perfectly matched layer (PML) uniaxial absorber for frequency domain finite element simulations has several advantages. In this paper we present the application of PML for microwave circuit simulations along with design guidelines to obtain a desired level of absorption. Different feeding techniques are also investigated for improved accuracy.

Gong, J.↗

Constitutive Modeling of the Facesheet to Core Interface in Honeycomb Sandwich Panels Subject to Mode I Delamination

A new cohesive zone traction-separation law, which includes the effects of fiber bridging, has been developed, implemented with a finite element (FE) model, and applied to simulate the delamination between the facesheet and core of a composite honeycomb sandwich panel. The proposed traction-separation law includes a standard initial cohesive component, which accounts for the initial interfacial stiffness and energy release rate, along with a new component to account for the fiber bridging contribution to the delamination process. Single cantilever beam tests on aluminum honeycomb sandwich panels with carbon fiber reinforced polymer facesheets were used to characterize and evaluate the new formulation and its finite element implementation. These tests, designed to evaluate the mode I toughness of the facesheet to core interface, exhibited significant fiber bridging and large crack process zones, giving rise to a concave downward concave upward pre-peak shape in the load-displacement curve. Unlike standard cohesive formulations, the proposed formulation captures this observed shape, and its results have been shown to be in excellent quantitative agreement with experimental load-displacement and apparent critical energy release rate results, representative of a payload fairing structure, as well as local strain fields measured with digital image correlation.

composites↗

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗

Coupling Finite Element and Finite Volume Simulation Within MOOSE

This work focuses on adding capability to couple finite element (FE) variables into finite volume (FV) physics within INL’s Multiphysics Object-Oriented Simulation Environment (MOOSE). This coupling can allow for improvement of multiphysics simulations where one set of physics is best suited for an FE discretization, while another set is best suited for an FV discretization. Electrohydrodynamics, which govern liquid metal reactor concepts and plasma dynamics, is a fitting example where the electromagnetic field equations are solved with FE and the fluid flow is solved with FV. The new FE to FV coupling method can be summarized as taking the element or face average of the FE variable value or gradient and applying that quantity directly in the FV equation objects.

97 MATHEMATICS AND COMPUTING↗

Micro-cantilever beam experiments and modeling in porous polycrystalline UO 2

Understanding the impact of microstructure on the thermo-mechanical behavior of oxide nuclear fuels is vital to predicting their performance through multiscale models. Evaluating the mechanical properties at the sub-grain length scale is key to developing these multiscale models. In this work, 3D finite element (FE) models were constructed to simulate the micrometer-scale bending of micro-cantilever beams fabricated using porous polycrystalline uranium dioxide (UO 2 ) and tested at room temperature. Here, the results showed that the porosity and elastic anisotropy of individual grains can play a significant role in determining the effective mechanical properties of the material deduced from the tests. Specifically, the porosity had a non-negligible effect, given that the pore size was of the same order of magnitude as the dimensions of the micro-beams. Correlations between load-deflection data, pore location, and elastic properties (effective Young's modulus) were investigated using UO 2 micro-beam FE models, where pore clusters were included and placed at different locations along the length of the beam. Results indicated that the presence of pore clusters near the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on the load-deflection behavior, with the porosity leading to a reduction of stiffness that is the largest for any location of the pore clusters. Furthermore, it was also found that pore clusters located towards the middle of the span and close to the end of the beam have a comparatively small effect on the load-deflection behavior. Therefore, it is concluded that accurate estimates of Young's modulus can be obtained from micro-cantilever experiments after accounting for porosity on the one third of the beam length close to the clamp. This, in turn, provides an avenue to improve microscale experiments and their analysis in porous, anisotropic elastic materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

MURR LEU structural and thermal hydraulics analyses: Part I – Preliminary irradiation thermo-mechanical behavior

The University of Missouri Research Reactor (MURR) is expected to be converted from highly enriched uranium (HEU, ≥ 20 wt% U-235) U-Al x dispersion fuel to low-enriched uranium (LEU, < 20 wt% U-235) with U-10Mo monolithic fuel. This work introduces high-fidelity irradiation thermo–mechanical (T-M) analysis of the MURR LEU focusing on changes in coolant channel gap thickness. Three-dimensional (3D) finite element (FE) models were developed to simulate the irradiation T-M behavior of the MURR LEU element with all 23 curved fuel plates, the two side plates, and the combs. It was shown that channel gap thickness changes were influenced not only by plate thickness variations due to fuel swelling and creep but also by the radial displacement of consecutive MURR LEU plates. Modeling the fuel element assembly captured side plate displacements, which were shown to reduce radial fuel plate displacements towards the convex side. The maximum local radial displacement in the element was predicted at the end of life (EOL) as 23.7 mil (602.0 µm) on the lateral centerline of plate 23 towards the convex side. The maximum stripe-averaged reduction in channel gap thickness, particularly relevant for thermal hydraulics (TH) safety analysis, was calculated as 15.9 mil (403.9 µm) in single-side heated channel 24 (the outermost channel). These results account for the thermal resistance from the oxide build-up on cladding surfaces which was shown to be up to 0.82 mil (20.8 µm) thick. It was demonstrated that accounting for oxide layer thermal resistance led to a 10 °C higher peak fuel temperature and a 4.4 mil (111.8 µm) greater maximum local radial displacement. The impact of the calculated channel gap thickness changes on the MURR LEU TH safety analysis is evaluated in Part II.

Irradiation thermo-mechanical analysis↗

Tardigrade-examples V0.1.0

Tardigrade-examples is a repository of computational workflows that exercise the Tardigrade software package. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling. These capabilities include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This repository investigates continuum upscaling of various direct numerical simulations (DNS) conducted in Abaqus finite element (FE), Ratel FE, and GEOS material point method (MPM) software. Verification of the upscaling workflow is first investigated by considering DNS of trivial stress states for homogeneous materials, results of which indicate that classical continuum behavior is recovered as expected. DNS of heterogeneous materials are then considered.

Allard, Thomas↗