Search NASA⌕ Search

SEARCH · Search NASA

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

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.↗

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↗

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.↗

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↗

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↗

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.↗

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.↗

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.↗

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.↗

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↗

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↗

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)↗

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↗

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.↗

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.↗

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.↗

Utah FORGE: Southwestern Utah Magnetotelluric (MT) Data

This comprehensive magnetotellurics (MT) dataset, which covers southwestern Utah, integrates 600 sites from various surveys, including those from the Utah FORGE, SubTER, and Play Fairway projects, all of which are linked below. The core of this dataset is the use of a 3D finite element (FE) algorithm for inversion imaging, allowing for detailed analysis of the MT data. Provided here are all .edi files, the 3D finite element mesh (with a readme file), graphic slices, node graphics, resistivity graphics, and reports on the acquisition areas and inversion.

15 GEOTHERMAL ENERGY↗

Finite element thermal-structural modeling of orbiting truss structures

A description of an integrated finite element (FE) thermal-structural approach for accurate and efficient modeling of large space structures is presented. A geometric model with a common discretization for all analyses is employed. It uses improved thermal elements and the results from the thermal analysis directly in the structural analysis without any intervening data processing. The differences between the conventional FE approach as implemented in large programs and an integrated FE approach currently under development are described. Considerations for thermal modeling of truss members is discussed and three thermal truss finite elements are presented. The performance of these elements was evaluated for typical truss members neglecting joint effects. A simple truss with metallic joints and composite members was studied to evaluate the effectiveness of the approach for realistic truss designs. A study of the effects of aluminum joints on the thermal deformations of a simple, plane truss with composite members showed that joint effects may be significant. Further study is needed to assess the role of joint effects on the deformation of large trusses.

Thornton, E. A.↗