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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
NASA NTRS2011

Shunted Piezoelectric Vibration Damping Analysis Including Centrifugal Loading Effects

Excessive vibration of turbomachinery blades causes high cycle fatigue problems which require damping treatments to mitigate vibration levels. One method is the use of piezoelectric materials as passive or active dampers. Based on the technical challenges and requirements learned from previous turbomachinery rotor blades research, an effort has been made to investigate the effectiveness of a shunted piezoelectric for the turbomachinery rotor blades vibration control, specifically for a condition with centrifugal rotation. While ample research has been performed on the use of a piezoelectric material with electric circuits to attempt to control the structural vibration damping, very little study has been done regarding rotational effects. The present study attempts to fill this void. Specifically, the objectives of this study are: (a) to create and analyze finite element models for harmonic forced response vibration analysis coupled with shunted piezoelectric circuits for engine blade operational conditions, (b) to validate the experimental test approaches with numerical results and vice versa, and (c) to establish a numerical modeling capability for vibration control using shunted piezoelectric circuits under rotation. Study has focused on a resonant damping control using shunted piezoelectric patches on plate specimens. Tests and analyses were performed for both non-spinning and spinning conditions. The finite element (FE) shunted piezoelectric circuit damping simulations were performed using the ANSYS Multiphysics code for the resistive and inductive circuit piezoelectric simulations of both conditions. The FE results showed a good correlation with experimental test results. Tests and analyses of shunted piezoelectric damping control, demonstrating with plate specimens, show a great potential to reduce blade vibrations under centrifugal loading.

Min, James B.↗
NASA NTRS2009

Aeroelastic-Acoustics Simulation of Flight Systems

This paper describes the details of a numerical finite element (FE) based analysis procedure and a resulting code for the simulation of the acoustics phenomenon arising from aeroelastic interactions. Both CFD and structural simulations are based on FE discretization employing unstructured grids. The sound pressure level (SPL) on structural surfaces is calculated from the root mean square (RMS) of the unsteady pressure and the acoustic wave frequencies are computed from a fast Fourier transform (FFT) of the unsteady pressure distribution as a function of time. The resulting tool proves to be unique as it is designed to analyze complex practical problems, involving large scale computations, in a routine fashion.

Gupta, kajal K.↗
NASA NTRS1995

Slot spiral antenna modeling using hybrid/mixed FE-BI technique

This report is concerned with the numerical simulation of the printed slot spiral antenna (SSA) using the finite element - boundary integral (FE-BI) analysis. It has been reported that the FE-BI technique is suited for modeling microstrip patch antennas of any shape, printed on layered planar structure or cylindrical platform and fed with a coaxial cable or a microstripline network underneath the radiating elements. However, direct application of the hybrid FE-BI technique to thin slot spirals requires excessive sampling rates to accurately simulate the geometry. To alleviate the meshing/modeling difficulties encountered with cavity-backed slot antennas, we describe a mixed finite element-boundary integral formulation. As in the past, the boundary integral is used to describe the radiation of the slot above the cavity. However the cavity is now modeled using a suitable mix of edge and node elements. The latter are used only at the aperture of the thin slot so that the nodes follow the center line of the slot. In this manner, regular size elements can be used regardless of the slot's width and any meshing restrictions are substantially relaxed. The proposed mixed element FE-BI formulation introduces three different computational regions and as expected this complicates the generation of the discrete system.

Gong, Jian↗
DOE OSTI2021

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Probabilistic Modeling of Intracranial Pressure Effects on Optic Nerve Biomechanics

Altered intracranial pressure (ICP) is involved/implicated in several ocular conditions: papilledema, glaucoma and Visual Impairment and Intracranial Pressure (VIIP) syndrome. The biomechanical effects of altered ICP on optic nerve head (ONH) tissues in these conditions are uncertain but likely important. We have quantified ICP-induced deformations of ONH tissues, using finite element (FE) and probabilistic modeling (Latin Hypercube Simulations (LHS)) to consider a range of tissue properties and relevant pressures.

biomechanics↗
NASA NTRSDate not supplied

Simulating High Energy Dynamic Impact of IM7/PEKK Continuous Fiber Laminated Thermoplastic Composite using Open Hole Coupon Experiments

The use of advanced thermoplastic composites (TPC) for structural applications in the aerospace and automotive industries has grown increasingly popular due to their high performance, high manufactured part output, and sustainability. Enhancing verification simulations of the progressive failure in TPCs at the coupon scale is needed to better validate advanced composite material models at larger length scales, such as element level panels. MAT_213, a tabulated composite material model, has been applied to simulate high velocity dynamic impact (HEDI) and is integrated in the LS-DYNA explicit finite element (FE) software. Previous efforts using MAT_213 have calibrated HEDI simulated failures with experimental data. Traditionally simulations have relied on a structured mesh to represent quasi-isotropic panels. For this study, a continuous fiber unidirectional tape, IM7/PEKK, was analyzed through a series of verification studies of notched laminate coupons, then validation analyses of quasi-isotropic panels under HEDI were compared with experimental data. Experiments using unidirectional TPC were performed to obtain tabulated stress vs strain curves using Digital Image Correlation (DIC) for input into MAT_213. Next, notched coupon experiments were used as an intermediate step to calibrate the damage and failure behavior in MAT_213. Lastly, investigations into differences between HEDI simulations using a fiber aligned mesh were compared with experimental failure patterns and damage sizes. The use of experimentally generated material behavior along with selective mesh layout improved HEDI predictions with minimal material parameter calibration.

Polymer matrix composites↗
NASA NTRS2017

Integrated Model of the Eye/Optic Nerve Head Biomechanical Environment

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a concern for long-duration space flight. Previously, it has been suggested that ocular changes observed in VIIP syndrome are related to the cephalad fluid shift that results in altered fluid pressures [1]. We are investigating the impact of changes in intracranial pressure (ICP) using a combination of numerical models, which simulate the effects of various environment conditions, including finite element (FE) models of the posterior eye. The specific interest is to understand how altered pressures due to gravitational changes affect the biomechanical environment of tissues of the posterior eye and optic nerve sheath. METHODS: Additional description of the numerical modeling is provided in the IWS abstract by Nelson et al. In brief, to simulate the effects of a cephalad fluid shift on the cardiovascular and ocular systems, we utilized a lumped-parameter compartment model of these systems. The outputs of this lumped-parameter model then inform boundary conditions (pressures) for a finite element model of the optic nerve head (Figure 1). As an example, we show here a simulation of postural change from supine to 15 degree head-down tilt (HDT), with primary outcomes being the predicted change in strains at the optic nerve head (ONH) region, specifically in the lamina cribrosa (LC), retrolaminar optic nerve, and prelaminar neural tissue (PLNT). The strain field can be decomposed into three orthogonal components, denoted as the first, second and third principal strains. We compare the peak tensile (first principal) and compressive (third principal) strains, since elevated strain alters cell phenotype and induces tissue remodeling. RESULTS AND CONCLUSIONS: Our lumped-parameter model predicted an IOP increase of c. 7 mmHg after 21 minutes of 15 degree HDT, which agreed with previous reports of IOP in HDT [1]. The corresponding FEM simulations predicted a relative increase in the magnitudes of the peak tensile and compressive strains in the lamina cribrosa of 42 and 43, respectively (Fig. 2). The corresponding changes in the optic nerve strains were 17 and 39, while in the PLNT they were 47 and 43. These magnitudes of relative elevations in peak strains may induce a phenotypic response in resident mechano-responsive resident cells [2]. This approach may be expanded to investigate other environmental changes (e.g. parabolic flight). Through our VIIP SCHOLAR project, we will validate and improve these integrated models by measuring patient-specific changes in optic nerve sheath geometry in patients with idiopathic intracranial hypertension before and after lumbar puncture and CSF removal.

biomechanics↗
NASA NTRSDate not supplied

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↗
DOE OSTI2022

Physics-informed machine learning assisted uncertainty quantification for the corrosion of dissimilar material joints

Jointing techniques like the Self-Piercing Riveting (SPR), Resistance Spot Welding (RSW) and Rivet-Weld (RW) joints are used for mass production of dissimilar material joints due to their high performance, short cycle time, and adaptability. However, the service life and safety usage of these joints can be largely impacted by the galvanic corrosion due to the difference in equilibrium potentials between the metals with the presence of electrolyte. Here, in this paper, we focus on Al-Fe galvanic corrosion and develop physics-informed machine learning based surrogate model for statistical corrosion analysis, which enables the reliability analysis of dissimilar material joints under corrosion environment. In this study, a physics-based finite element (FE) corrosion model has been developed to simulate the galvanic corrosion between a Fe cathode and an Al anode. Geometric and environmental factors including crevice gap, roughness of anode, conductivity, and the temperature of the electrolyte are investigated. Further, a thorough Uncertainty Quantification (UQ) analysis is conducted for the overall corrosion behavior of the Fe-Al joints. It is found that the electrolyte conductivity has the largest effects on the material loss and needs to be managed closely for better corrosion control. This will help in designing and manufacturing joints with improved corrosion performance.

42 ENGINEERING↗
NASA NTRS2023

Modeling and Simulation Credibility Assessments of Whole-Body Finite Element Computational Models for Use in NASA Extravehicular Activity Applications

Computational finite element (FE) models are used in suited astronaut injury risk assessments; however, these models’ verification, validation, and credibility (VV&C) procedures for simulating injuries in altered gravity environments are limited. Our study conducts VV&C assessments of THUMS and Elemance whole-body FE models for predicting suited astronaut injury biomechanics using eight credibility factors, as per NASA-STD-7009A. Credibility factor ordinal scores are assigned by reviewing existing documentation describing VV&C practices, and credibility sufficiency thresholds are assigned based on input from subject matter experts. Our results show the FE models are credible for suited astronaut injury investigation in specific ranges of kinematic and kinetic conditions correlating to highway and contact sports events. Nevertheless, these models are deficient when applied outside these ranges. Several credibility elevation strategies are prescribed to improve models’ credibility for the NASA-centric application domain.

Finite Element↗
NASA NTRS2023

Modeling and simulation credibility assessments of whole-body finite element computational models for use in NASA extravehicular activity applications

Computational finite element (FE) models are used in suited astronaut injury risk assessments; however, these models’ verification, validation, and credibility (VV&C) procedures for simulating injuries in altered gravity environments are limited. Our study conducts VV&C assessments of THUMS and Elemance whole-body FE models for predicting suited astronaut injury biomechanics using eight credibility factors, as per NASA-STD-7009A. Credibility factor ordinal scores are assigned by reviewing existing documentation describing VV&C practices, and credibility sufficiency thresholds are assigned based on input from subject matter experts. Our results show the FE models are credible for suited astronaut injury investigation in specific ranges of kinematic and kinetic conditions correlating to highway and contact sports events. Nevertheless, these models are deficient when applied outside these ranges. Several credibility elevation strategies are prescribed to improve models’ credibility for the NASA-centric application domain.

Finite Element↗
NASA NTRS2001

The Applicability of the Generalized Method of Cells for Analyzing Discontinuously Reinforced Composites

The paper begins with a short overview of the recent work done in the field of discontinuous reinforced composites, focusing on the different parameters which influence the material behavior of discontinuous reinforced composites, as well as the various analysis approaches undertaken. Based on this overview it became evident that in order to investigate the enumerated effects in an efficient and comprehensive manner, an alternative approach to the computationally intensive finite-element based micromechanics approach is required. Therefore, an investigation is conducted to demonstrate the utility of utilizing the generalized method of cells (GMC), a semi-analytical micromechanics-based approach, to simulate the elastic and elastoplastic material behavior of aligned short fiber composites. The results are compared with simulations using other micromechanical based mean field models and finite element (FE) unit cell models found in the literature given elastic material behavior, as well as finite element unit cell and a new semianalytical elastoplastic shear lag model in the inelastic range. GMC is shown to definitely have a window of applicability when simulating discontinuously reinforced composite material behavior.

Pahr, D. H.↗
NASA NTRS2001

The Applicability of the Generalized Method of Cells for Analyzing Discontinuously Reinforced Composites

The paper begins with a short overview of the recent work done in the field of discontinuous reinforced composites, focusing on the different parameters which influence the material behavior of discontinuous reinforced composites, as well as the various analysis approaches undertaken. Based on this overview it became evident, that in order to investigate the enumerated effects in an efficient and comprehensive manner, an alternative approach to the computationally intensive finite-element based micromechanics approach is required. Therefore, an investigation is conducted to demonstrate the utility of utilizing the generalized method of cells (GMC), a semi-analytical micromechanics-based approach, to simulate the elastic and elastoplastic material behavior of aligned short fiber composites. The results are compared with (1) simulations using other micromechanical based mean field models and finite element (FE) unit cell models found in the literature given elastic material behavior, as well as (2) finite element unit cell and a new semianalytical elastoplastic shear lag model in the inelastic range. GMC is shown to definitely have a window of applicability when simulating discontinuously reinforced composite material behavior.

Pahr, D. H.↗
DOE OSTI2022

Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea

We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.

depressurization↗
DOE OSTI2025

Finite Element Analysis of the TRUST Nonlinear Dynamics Testbed

This paper builds on prior work conducted within the Los Alamos National Laboratory (LANL) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) program. Specifically, it builds on finite element (FE) modeling efforts for the TRUST program’s Nonlinear Dynamics (ND) testbed. Historically, the FE model for the ND testbed has exclusively utilized a Lanczos eigensolver that linearly extracts the system’s natural frequencies. This paper investigates the Abaqus 2024’s explicit dynamic solver, which implements a central difference explicit solver. The central difference method used in Abaqus can capture nonlinear material responses in structural dynamic simulations making it suitable for the ND FE model.

42 ENGINEERING↗
DOE OSTI2021

Validation of a Pressure Dependent PBX 9501 Cookoff Model

A pressure dependent cookoff model for PBX 9501 was developed. This cookoff model was implemented in the finite element (FE) heat transfer code Aria and was used to simulate a set of cookoff experiments. This set is called the Large Scale Annular Cookoff (LSAC) experiments. Three dimensional and axisymmetric heat conduction FE models were implemented to simulate these experiments. The predictive ability of the PBX 9501 pressure dependent cookoff model was evaluated using Latin Hypercube sampling (LHS). Predictions of the thermal times to ignition and temperatures are compared with the experiments.

42 ENGINEERING↗