Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element (FE)”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Astronaut Bone Medical Standards Derived from Finite Element (FE) Models of QCT Scans from Population Studies

This work was accomplished in support of the Finite Element [FE] Strength Task Group, NASA Johnson Space Center [JSC], Houston, TX. This group was charged with the task of developing rules for using finite-element [FE] bone-strength measures to construct operating bands for bone health that are relevant to astronauts following exposure to spaceflight. FE modeling is a computational tool used by engineers to estimate the failure loads of complex structures. Recently, some engineers have used this tool to characterize the failure loads of the hip in population studies that also monitored fracture outcomes. A Directed Research Task was authorized in July, 2012 to investigate FE data from these population studies to derive these proposed standards of bone health as a function of age and gender. The proposed standards make use of an FE-based index that integrates multiple contributors to bone strength, an expanded evaluation that is critical after an astronaut is exposed to spaceflight. The current index of bone health used by NASA is the measurement of areal BMD. There was a concern voiced by a research and clinical advisory panel that the sole use of areal BMD would be insufficient to fully evaluate the effects of spaceflight on the hip. Hence, NASA may not have a full understanding of fracture risk, both during and after a mission, and may be poorly estimating in-flight countermeasure efficacy. The FE Strength Task Group - composed of principal investigators of the aforementioned population studies and of FE modelers -donated some of its population QCT data to estimate of hip bone strength by FE modeling for this specific purpose. Consequently, Human Health Countermeasures [HHC] has compiled a dataset of FE hip strengths, generated by a single FE modeling approach, from human subjects (approx.1060) with ages covering the age range of the astronauts. The dataset has been analyzed to generate a set of FE strength cutoffs for the following scenarios: a) Qualify an applicant for astronaut candidacy, b) Qualify an astronaut for a long-duration (LD) mission, c) Qualify a veteran LD astronaut for a second LD mission, and d) Establish a non-permissible, minimum hip strength following a given mission architecture. This abstract will present the FE-based standards accepted by the FE Strength Task Group for its recommendation to HHC in January 2015.

Sibonga, J. D.↗
DOE OSTI2023

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗
NASA NTRS2004

Eccentric loading of microtensile specimens

Ceramic materials have a lower density than most metals and are capable of performing at extremely high temperatures. The utility of these materials is obvious; however, the fracture strength of brittle materials is not easily predicted and often varies greatly. Characteristically, brittle materials lack ductility and do not yield as other materials. Ceramics materials are naturally populated with microscopic cracks due to fabrication techniques. Upon application of a load, stress concentration occurs at the root of these cracks and fracture will eventually occur at some not easily predicted strength. In order to use ceramics in any application some design methodology must exist from which a component can be placed into service. This design methodology is CARES/LIFE (Ceramics Analysis and Reliability Evaluation of Structures) which has been developed and refined at NASA over the last several decades. The CARES/LIFE computer program predicts the probability of failure of a ceramic component over its service life. CARES combines finite element results from a commercial FE (finite element) package such as ANSYS and experimental results to compute the abovementioned probability of failure. Over the course of several tests CARES has had great success in predicting the life of various ceramic components and has been used throughout industry. The latest challenge is to verify that CARES is valid for MEMS (Micro-Electro Mechanical Systems). To investigate a series of microtensile specimens were fractured in the laboratory. From this data, material parameters were determined and used to predict a distribution of strength for other specimens that exhibit a known stress concentration. If the prediction matches the experimental results then these parameters can be applied to a desired component outside of the laboratory. During testing nearly half of the tensile Specimens fractured at a location that was not expected and hence not captured in the FE model. It has been my duty to investigate the nature of this phenomenon in hopes of finding a better correlation between theory and empirical results. To investigate I built complete FE models of all of the tensile specimens using ANSYS. It is suspected that some misalignment naturally occurs during testing and thus additional bending stresses are present in the specimens. I modeled this eccentric loading and ran several FE trials using ANSYS/PDS (a probabilistic design system in ANSYS). My objective this summer has been familiarize myself with the CARES/LIFE program in hopes of using it in conjunction with ANSYS to help verify that CARES is applicable to MEMS-scale (greater that 1 micron, less than 1 millimeter) components.

Trapp, Mark A.↗
DOE OSTI2021

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

Design of Rigid-Flex PCB Robotics Leveraging Validated Finite Element Simulations

The use of rigid-flex printed circuit board (PCB) as primary structure has the potential to reduce the weight and volume of robotic systems. In the case of robotics for interplanetary exploration, these systems can leverage origami-inspired folding for increased mobility options and reduced storage volume. Folding rigid-flex PCB robotics can be constructed with rigid PCB connected by short Nomex fabric hinges coupled with flex PCB ribbon cables that permits enhanced system flexibility and energy dissipation to promote impact survivability. This paper presents a design methodology of rigid-flex PCB systems with an emphasis on impact resistance. The design process considers solder joint adequacy, panel bending, and fracture using a finite element (FE) model. The proposed design methodology is developed using a case study with NASA JPL's Pop-Up Folding Flat Explorer Robot (PUFFER). First, the finite-element (FE) modeling methodology is presented with consideration to both frequency and time-domain modeling applications, which include operational self-contact analysis and high impact scenarios. The time-domain impact modeling methodology utilizes hyperelastic material properties for the Nomex hinges. This modeling method is validated using image correlation of PUFFER drop tests. A flowchart is presented to guide users through a validated Abaqus modeling procedure for highly flexible rigid-flex systems. Next, a case study is presented in which PUFFER is subject to drop heights representative of falls into Lunar pits and then the design is refined for a more optimum impact performance. Finally, the results of the case study are used to inform a generalized design methodology for rigid-flex PCB robotics subject to high-impact loads with the considerations presented.

John Bell↗
NASA NTRS2022

Design of Rigid-Flex PCB Robotics Leveraging Validated Finite Element Simulations

The use of rigid-flex printed circuit board (PCB) as primary structure has the potential to reduce the weight and volume of robotic systems. In the case of robotics for interplanetary exploration, these systems can leverage origamiinspired folding for increased mobility options and reduced storage volume. Folding rigid-flex PCB robotics can be constructed with rigid PCB connected by short Nomex fabric hinges coupled with flex PCB ribbon cables that permits enhanced system flexibility and energy dissipation to promote impact survivability. This paper presents a design methodology of rigidflex PCB systems with an emphasis on impact resistance. The design process considers solder joint adequacy, panel bending, and fracture using a finite element (FE) model. The proposed design methodology is developed using a case study with NASA JPL’s Pop-Up Folding Flat Explorer Robot (PUFFER). First, the finite-element (FE) modeling methodology is presented with consideration to both frequency and time-domain modeling applications, which include operational self-contact analysis and high impact scenarios. The time-domain impact modeling methodology utilizes hyperelastic material properties for the Nomex hinges. This modeling method is validated using image correlation of PUFFER drop tests. A flowchart is presented to guide users through a validated Abaqus modeling procedure for highly flexible rigid-flex systems. Next, a case study is presented in which PUFFER is subject to drop heights representative of falls into Lunar pits and then the design is refined for a more optimum impact performance. Finally, the results of the case study are used to inform a generalized design methodology for rigid-flex PCB robotics subject to high-impact loads with the considerations presented.

de la Croix, Jean-Pierre↗
NASA NTRS2002

Coupling Finite Element and Meshless Local Petrov-Galerkin Methods for Two-Dimensional Potential Problems

A coupled finite element (FE) method and meshless local Petrov-Galerkin (MLPG) method for analyzing two-dimensional potential problems is presented in this paper. The analysis domain is subdivided into two regions, a finite element (FE) region and a meshless (MM) region. A single weighted residual form is written for the entire domain. Independent trial and test functions are assumed in the FE and MM regions. A transition region is created between the two regions. The transition region blends the trial and test functions of the FE and MM regions. The trial function blending is achieved using a technique similar to the 'Coons patch' method that is widely used in computer-aided geometric design. The test function blending is achieved by using either FE or MM test functions on the nodes in the transition element. The technique was evaluated by applying the coupled method to two potential problems governed by the Poisson equation. The coupled method passed all the patch test problems and gave accurate solutions for the problems studied.

Chen, T.↗
NASA NTRS1992

Coupling finite and boundary element methods for two-dimensional potential problems

A finite element-boundary element (FE-BE) coupling method based on a weighted residual variational method is presented for potential problems, governed by either the Laplace or the Poisson equations. In this method, a portion of the domain of interest is modeled by finite elements (FE) and the remainder of the region by boundary elements (BE). Because the BE fundamental solutions are valid for infinite domains, a procedure that limits the effects of the BE fundamental solution to a small region adjacent to the FE region, called the transition region (TR), is developed. This procedure involves a judicious choice of functions called the transition (T) functions that have unit values on the BE-TR interface and zero values on the FE-TR interface. The present FE-BE coupling algorithm is shown to be independent of the extent of the transition region and the choice of the transition functions. Therefore, transition regions that extend to only one layer of elements between FE and BE regions and the use of simple linear transition functions work well.

Krishnamurthy, T.↗
NASA NTRS2000

Effect of Crystal Orientation on Fatigue Failure of Single Crystal Nickel Base Turbine Blade Superalloys

High Cycle Fatigue (HCF) induced failures in aircraft gas-turbine engines is a pervasive problem affecting a wide range of components and materials. HCF is currently the primary cause of component failures in gas turbine aircraft engines. Turbine blades in high performance aircraft and rocket engines are increasingly being made of single crystal nickel superalloys. Single-crystal Nickel-base superalloys were developed to provide superior creep, stress rupture, melt resistance and thermomechanical fatigue capabilities over polycrystalline alloys previously used in the production of turbine blades and vanes. Currently the most widely used single crystal turbine blade superalloys are PWA 1480/1493 and PWA 1484. These alloys play an important role in commercial, military and space propulsion systems. PWA1493, identical to PWA1480, but with tighter chemical constituent control, is used in the NASA SSME (Space Shuttle Main Engine) alternate turbopump, a liquid hydrogen fueled rocket engine. Objectives for this paper are motivated by the need for developing failure criteria and fatigue life evaluation procedures for high temperature single crystal components, using available fatigue data and finite element modeling of turbine blades. Using the FE (finite element) stress analysis results and the fatigue life relations developed, the effect of variation of primary and secondary crystal orientations on life is determined, at critical blade locations. The most advantageous crystal orientation for a given blade design is determined. Results presented demonstrates that control of secondary and primary crystallographic orientation has the potential to optimize blade design by increasing its resistance to fatigue crack growth without adding additional weight or cost.

Arakere, Nagaraj K.↗
DOE OSTI2022

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗
NASA NTRS1999

User's Guide for ENSAERO_FE Parallel Finite Element Solver

A high fidelity parallel static structural analysis capability is created and interfaced to the multidisciplinary analysis package ENSAERO-MPI of Ames Research Center. This new module replaces ENSAERO's lower fidelity simple finite element and modal modules. Full aircraft structures may be more accurately modeled using the new finite element capability. Parallel computation is performed by breaking the full structure into multiple substructures. This approach is conceptually similar to ENSAERO's multizonal fluid analysis capability. The new substructure code is used to solve the structural finite element equations for each substructure in parallel. NASTRANKOSMIC is utilized as a front end for this code. Its full library of elements can be used to create an accurate and realistic aircraft model. It is used to create the stiffness matrices for each substructure. The new parallel code then uses an iterative preconditioned conjugate gradient method to solve the global structural equations for the substructure boundary nodes.

Eldred, Lloyd B.↗
NASA NTRS2013

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

Thermal Fault Tolerance Analysis of Carbon Fiber Rope Barrier Systems for Use in the Reusable Solid Rocket Motor ( RSRM) Nozzle Joints

Carbon Fiber Rope (CFR) thermal barrier systems are being considered for use in several RSRM (Reusable Solid Rocket Motor) nozzle joints as a replacement for the current assembly gap close-out process/design. This study provides for development and test verification of analysis methods used for flow-thermal modeling of a CFR thermal barrier subject to fault conditions such as rope combustion gas blow-by and CFR splice failure. Global model development is based on a 1-D (one dimensional) transient volume filling approach where the flow conditions are calculated as a function of internal 'pipe' and porous media 'Darcy' flow correlations. Combustion gas flow rates are calculated for the CFR on a per-linear inch basis and solved simultaneously with a detailed thermal-gas dynamic model of a local region of gas blow by (or splice fault). Effects of gas compressibility, friction and heat transfer are accounted for the model. Computational Fluid Dynamic (CFD) solutions of the fault regions are used to characterize the local flow field, quantify the amount of free jet spreading and assist in the determination of impingement film coefficients on the nozzle housings. Gas to wall heat transfer is simulated by a large thermal finite element grid of the local structure. The employed numerical technique loosely couples the FE (Finite Element) solution with the gas dynamics solution of the faulted region. All free constants that appear in the governing equations are calibrated by hot fire sub-scale test. The calibrated model is used to make flight predictions using motor aft end environments and timelines. Model results indicate that CFR barrier systems provide a near 'vented joint' style of pressurization. Hypothetical fault conditions considered in this study (blow by, splice defect) are relatively benign in terms of overall heating to nozzle metal housing structures.

Clayton, J. Louie↗
NASA NTRS2017

Optimal Reference Strain Structure for Studying Dynamic Responses of Flexible Rockets

In the proposed paper, the optimal design of reference strain structures (RSS) will be performed targeting for the accurate observation of the dynamic bending and torsion deformation of a flexible rocket. It will provide the detailed description of the finite-element (FE) model of a notional flexible rocket created in MSC.Patran. The RSS will be attached longitudinally along the side of the rocket and to track the deformation of the thin-walled structure under external loads. An integrated surrogate-based multi-objective optimization approach will be developed to find the optimal design of the RSS using the FE model. The Kriging method will be used to construct the surrogate model. For the data sampling and the performance evaluation, static/transient analyses will be performed with MSC.Natran/Patran. The multi-objective optimization will be solved with NSGA-II to minimize the difference between the strains of the launch vehicle and RSS. Finally, the performance of the optimal RSS will be evaluated by checking its strain-tracking capability in different numerical simulations of the flexible rocket.

finite-element (FE)↗
DOE OSTI2024

A Study on Particle Trajectory Error in Finite-Element Particle-in-Cell Algorithms

Particle-in-cell (PIC) algorithms are widely used for the simulation of kinetic plasmas. PIC algorithms account for the interaction between charged particles in a plasma and the electromagnetic field in ambient space, including self-field interactions. The objective of this article is to study the error in charged particle trajectories present in finite-element (FE)-based PIC algorithms on unstructured meshes. We study how the trajectory error behaves according to the FE mesh resolution and the matrix solver employed in the FE algorithm. The study is performed by considering a trajectory established by a parabolic electric potential and an axial magnetic force acting on the charged particle. Under a proper choice combination of electric and magnetic field parameters, the 2-D particle trajectories comprise closed orbits. Numerical errors cause small perturbations on the trajectories, with cumulative effects. As a result, the resulting orbital trajectories exhibit distortions including spurious apsidal precession. These distortions provide a clear imprint of the numerical error. Here, we also study the numerical error in a quantitative fashion by computing the distance norm between the trajectories obtained by the exact fields and numerical fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.↗