Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Simulating the Impact Response of Three Full-Scale Crash Tests of Cessna 172 Aircraft

During the summer of 2015, a series of three full-scale crash tests were performed at the Landing and Impact Research Facility located at NASA Langley Research Center of Cessna 172 aircraft. The first test (Test 1) represented a flare-to-stall emergency or hard landing onto a rigid surface. The second test (Test 2) represented a controlled-flight- into-terrain (CFIT) with a nose down pitch attitude of the aircraft, which impacted onto soft soil. The third test (Test 3) also represented a CFIT with a nose up pitch attitude of the aircraft, which resulted in a tail strike condition. Test 3 was also conducted onto soft soil. These crash tests were performed for the purpose of evaluating the performance of Emergency Locator Transmitters and to generate impact test data for model calibration. Finite element models were generated and impact analyses were conducted to simulate the three impact conditions using the commercial nonlinear, transient dynamic finite element code, LS-DYNA®. The objective of this paper is to summarize test-analysis results for the three full-scale crash tests.

Jackson, Karen E.↗
DOE OSTI2016

BISON Fuel Performance Analysis of FeCrAl cladding with updated properties

In order to improve the accident tolerance of light water reactor (LWR) fuel, alternative cladding materials have been proposed to replace zirconium (Zr)-based alloys. Of these materials, there is a particular focus on iron-chromium-aluminum (FeCrAl) alloys due to much slower oxidation kinetics in high-temperature steam than Zr-alloys. This should decrease the energy release due to oxidation and allow the cladding to remain integral longer in the presence of high temperature steam, making accident mitigation more likely. As a continuation of the development for these alloys, suitability for normal operation must also be demonstrated. This research is focused on modeling the integral thermo-mechanical performance of FeCrAl-cladded fuel during normal reactor operation. Preliminary analysis has been performed to assess FeCrAl alloys (namely Alkrothal 720 and APMT) as a suitable fuel cladding replacement for Zr-alloys, using the MOOSE-based, finite-element fuel performance code BISON and the best available thermal-mechanical and irradiation-induced constitutive properties. These simulations identify the effects of the mechanical-stress and irradiation response of FeCrAl, and provide a comparison with Zr-alloys. In comparing these clad materials, fuel rods have been simulated for normal reactor operation and simple steady-state operation. Normal reactor operating conditions target the cladding performance over the rod lifetime (~4 cycles) for the highest-power rod in the highest-power fuel assembly under reactor power maneuvering. The power histories and axial temperature profiles input into BISON were generated from a neutronics study on full-core reactivity equivalence for FeCrAl using the 3D full core simulator NESTLE. Evolution of the FeCrAl cladding behavior over time is evaluated by using steady-state operating conditions such as a simple axial power profile, a constant cladding surface temperature, and a constant fuel power history. The fuel rod designs and operating conditions used are based off the Peach Bottom BWR and design consideration was given to minimize the neutronic penalty of the FeCrAl cladding by changing fuel enrichment and cladding thickness. As this study progressed, systematic parametric analysis of the fuel and cladding creep responses were also performed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Edge equilibrium code for tokamaks

The edge equilibrium code (EEC) described in this paper is developed for simulations of the near edge plasma using the finite element method. It solves the Grad-Shafranov equation in toroidal coordinate and uses adaptive grids aligned with magnetic field lines. Hermite finite elements are chosen for the numerical scheme. A fast Newton scheme which is the same as implemented in the equilibrium and stability code (ESC) is applied here to adjust the grids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Multidisciplinary tailoring of hot composite structures

A computational simulation procedure is described for multidisciplinary analysis and tailoring of layered multi-material hot composite engine structural components subjected to simultaneous multiple discipline-specific thermal, structural, vibration, and acoustic loads. The effect of aggressive environments is also simulated. The simulation is based on a three-dimensional finite element analysis technique in conjunction with structural mechanics codes, thermal/acoustic analysis methods, and tailoring procedures. The integrated multidisciplinary simulation procedure is general-purpose including the coupled effects of nonlinearities in structure geometry, material, loading, and environmental complexities. The composite material behavior is assessed at all composite scales, i.e., laminate/ply/constituents (fiber/matrix), via a nonlinear material characterization hygro-thermo-mechanical model. Sample tailoring cases exhibiting nonlinear material/loading/environmental behavior of aircraft engine fan blades, are presented. The various multidisciplinary loads lead to different tailored designs, even those competing with each other, as in the case of minimum material cost versus minimum structure weight and in the case of minimum vibration frequency versus minimum acoustic noise.

Singhal, Surendra N.↗
NASA NTRS1994

Comparison of Gap Elements and Contact Algorithm for 3D Contact Analysis of Spiral Bevel Gears

Three dimensional stress analysis of spiral bevel gears in mesh using the finite element method is presented. A finite element model is generated by solving equations that identify tooth surface coordinates. Contact is simulated by the automatic generation of nonpenetration constraints. This method is compared to a finite element contact analysis conducted with gap elements.

Bibel, G. D.↗
NASA NTRS2020

Overview of Coupon Testing of an IM7/8552 Composite Required to Characterize High-Energy Impact Dynamic Material Models

As part of NASA’s Advanced Composites Project (ACP), a significant research effort was established to develop and mature analysis methods for composite materials particularly applicable to high-energy dynamic impact applications. An overall goal was to develop and advance analysis methods that can be characterized using fundamental coupon level tests. Based on this fundamental characterization, simulations at the higher levels of a building block approach can be purely predictive. The tests described in this report were performed on a unidirectional carbon fiber polymer matrix composite, IM7/8552 tape, as this material system is representative of a modern aerospace grade carbon fiber tape system. These tests were chosen to provide fundamental material characterization data for three representative progressive damage and failure analysis material models available within the LS-DYNA commercial transient dynamic finite element code for use in the simulation of the dynamic and impact response of composite materials. Two of the material models, MAT162 and MAT261, are well established material models within LS-DYNA. The third material model, MAT213, is under active development. Testing consisted primarily of ASTM-type standard tests, and the tests that were performed were consistent with the needs of all three models. Additionally, coupon level tests were performed that can be used to characterize a delamination model (e.g., cohesive zone or tie-break contact). While, the tests described in this report provide fundamental material characterization data for all three of the material models mentioned above, the data is applicable to additional material models not discussed.

Justusson, Brian P.↗
DOE OSTI2013

Development of a graphite probe calorimeter for absolute clinical dosimetry

The aim of this work is to present the numerical design optimization, construction, and experimental proof of concept of a graphite probe calorimeter (GPC) conceived for dose measurement in the clinical environment (U.S. provisional patent 61/652,540). A finite element method (FEM) based numerical heat transfer study was conducted using a commercial software package to explore the feasibility of the GPC and to optimize the shape, dimensions, and materials used in its design. A functioning prototype was constructed inhouse and used to perform dose to water measurements under a 6 MV photon beam at 400 and 1000 MU/min, in a thermally insulated water phantom. Heat loss correction factors were determined using FEM analysis while the radiation field perturbation and the graphite to water absorbed dose conversion factors were calculated using Monte Carlo simulations. The difference in the average measured dose to water for the 400 and 1000 MU/min runs using the TG-51 protocol and the GPC was 0.2% and 1.2%, respectively. Heat loss correction factors ranged from 1.001 to 1.002, while the product of the perturbation and dose conversion factors was calculated to be 1.130. The combined relative uncertainty was estimated to be 1.4%, with the largest contributors being the specific heat capacity of the graphite (type B, 0.8%) and the reproducibility, defined as the standard deviation of the mean measured dose (type A, 0.6%). By establishing the feasibility of using the GPC as a practical clinical absolute photon dosimeter, this work lays the foundation for further device enhancements, including the development of an isothermal mode of operation and an overall miniaturization, making it potentially suitable for use in small and composite radiation fields. It is anticipated that, through the incorporation of isothermal stabilization provided by temperature controllers, a subpercent overall uncertainty will be achieved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2010

Modification of a Macromechanical Finite-Element Based Model for Impact Analysis of Triaxially-Braided Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. For the current analytical approach, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The constitutive model requires stiffness and strength properties of an equivalent unidirectional composite. Simplified micromechanics methods are used to determine the equivalent stiffness properties, and results from coupon level tests on the braided composite are utilized to back out the required strength properties. Simulations of quasi-static coupon tests of several representative braided composites are conducted to demonstrate the correlation of the model. Impact simulations of a represented braided composites are conducted to demonstrate the capability of the model to predict the penetration velocity and damage patterns obtained experimentally.

Goldberg, Robert K.↗
NASA NTRSDate not supplied

Design and Analysis of Shape Memory Spring Tires for Martian and Lunar Rover Vehicles

Shape memory alloys (SMAs) have played an important role in various innovative engineering and medical applications, such as aerospace actuators, vibration damping devices, and coronary stents. In applications, shape memory alloys are commonly utilized in two fundamentally different ways: (i) making use of the superelasticity/ pseudoelasticity (SE/PE) phenomena, as in applications in biomedical engineering, and (ii) taking advantages of the shape memory effect (SME), as is used for actuators. Their ability to act in such vastly different capacities is mainly due to their unique capability to recover large amounts of deformation produced by either applied stresses or temperature changes. One recent emerging application in use of SMAs has been in the area of non-pneumatic tire designs for Martian or Lunar roving vehicles. These vehicles require tires that are capable of traversing rugged terrain while withstanding extreme temperatures and atmospheric conditions. Inspired by the flexible wire mesh tires used on three Lunar Roving Vehicle (LRV) missions to the Moon on Apollo 15, 16, and 17, a new compliant tire technology was developed by the NASA Glenn Research Center (GRC) and Goodyear Tire & Rubber, known as the Spring Tire. The Spring Tire consists of several hundred coiled springs woven into a flexible mesh and formed into the shape of a tire. Like the LRV wire mesh tires, the original Spring Tires were made from spring steel and were prone to permanent deformation when undergoing high localized loads. Later, a new iteration of this technology was invented, known as the ‘Superelastic Tire’. This new technology incorporated the use of superelastic SMA springs, which could effectively undergo approximately 30 times more reversible deformation than the steel spring. It also provided even greater durability and allowed for more flexibility in design, such as the use of other structural forms to reduce mass or increase load carrying capacity. Because of the unique nature of both the SMA material and the complex interactions between the springs, designing Spring Tires for a specific application requires extensive effort. Historically, design decisions have relied on full-scale empirical testing; however, this is very expensive and time consuming, especially when multiple iterations are needed. Therefore, developing a large-scale, robust, and predictive numerical model entailing complex spring interactions and the shape memory material behavior within a tire construct is the first essential step towards a successful design program. The current work focuses on implementation of the user-defined Shape Memory Alloy (SMA) model, otherwise known as SMA-GVIPs, in the Finite Element analysis (FEA) program ABAQUS for large-scale simulations of the GRC-developed Spring Tires made of SMA. The novelty of this work lies in the thorough, detail-oriented, and computationally efficient finite element analysis of full-scale SMA tires. A well-thought material characterization plan followed by model validation and a systemic sensitivity study on SMA tires has never been reported in the previous literature. The main objective of this study is to help the team improve and optimize the structural design of the SMA tires through in-depth numerical analysis and sensitivity studies. Various design variables (wire diameter, coil diameter, pitch, bead angle, and number of springs) were varied to study their influence on the global load-displacement response of the tire construct. A detailed investigation of the three-dimensional stress states was also carried out to enhance our understanding of the local changes as the tire goes through global deformation. It was concluded that a robust numerical model with a good predictive capability, together with a thoughtfully crafted sensitivity study can result in improved design iterations required to reach a desired tire performance while, significantly reducing manufacturing, labor and testing expenses. A summary of the Finite Element (FE) model construction will be presented together with a description of the user-defined SMA model, characterization process, experimental results, model validation, and numerical sensitivity study results.

shape memory alloys↗
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 NTRS1996

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Yu, Sheng-Tao↗
DOE OSTI2015

Liquid crystal free energy relaxation by a theoretically informed Monte Carlo method using a finite element quadrature approach

A theoretically informed Monte Carlo method is proposed for Monte Carlo simulation of liquid crystals on the basis of theoretical representations in terms of coarse-grained free energy functionals. The free energy functional is described in the framework of the Landau-de Gennes formalism. A piecewise finite element discretization is used to approximate the alignment field, thereby providing an excellent geometrical representation of curved interfaces and accurate integration of the free energy. The method is suitable for situations where the free energy functional includes highly non-linear terms, including chirality or high-order deformation modes. Finally, the validity of the method is established by comparing the results of Monte Carlo simulations to traditional Ginzburg-Landau minimizations of the free energy using a finite difference scheme, and its usefulness is demonstrated in the context of simulations of chiral liquid crystal droplets with and without nanoparticle inclusions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Geometric Sensitivity of Residual Fields in Metal Additive Manufacturing

The sensitivity of geometrical parameters on the residual stress and distortion as computed by a metal additive manufacturing simulation is presented. A nozzle-like geometry was chosen with its length, wall thickness, and outer radius varied. A Latin hypercube sampling of 300 specimens produced the combinations of the geometrical parameters. Full thermomechanical simulations predicted the residual fields. Statistical analysis including regression curves shows that there is a strong correlation between these residual fields and the geometry, hinting that it may be possible to predictably influence the residual displacements and residual stresses by varying the overall geometry.

Multiscale Modeling↗
DOE OSTI2019

Finite Element Analysis and In-Situ Measurement of Out-of-Plane Distortion in Thin Plate TIG Welding

Transient distortion of thin plate in the welding process usually has a complicated mode and large magnitude. Quantitative measurement and prediction of full-field distortion are challenging and rarely reported up to now. In this study, the out-of-plane distortion of a thin plate during the Tungsten Inert Gas (TIG) welding process was measured using the digital image correlation (DIC) method. A simulation model based on thermal elastic–plastic finite element method (FEM) and DIC measured geometric imperfection were developed for accurate prediction of the transient welding distortion. The numerical results and experimental data agreed very well in both out-of-plane deformation modes and magnitudes of the plate at different stages of welding. The maximum out-of-plane distortion was larger than 4 mm during welding which can cause instability of arc length and heat input. The distance change between welding torch and plate surface was investigated under different initial deflections of the plate before welding. The plate with flat geometry shows the minimum transient and final gap change. In addition, the relationship between heat input and welding distortion was clarified through a series of numerical analyses. As a result, optimization of welding heat input can be performed based on numerical results to avoid excessive welding distortion.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Investigation of uncertainty in CO 2 reservoir models: A sensitivity analysis of relative permeability parameter values

Numerical reservoir models of CO 2 injection in saline formations rely on parameterization of laboratory-measured pore-scale processes. Here, we have performed a parameter sensitivity study and Monte Carlo simulations to determine the normalized change in total CO 2 injected using the finite element heat and mass-transfer code (FEHM) numerical reservoir simulator. Experimentally measured relative permeability parameter values were used to generate distribution functions for parameter sampling. The parameter sensitivity study analyzed five different levels for each of the relative permeability model parameters. All but one of the parameters changed the CO 2 injectivity by <10%, less than the geostatistical uncertainty that applies to all large subsurface systems due to natural geophysical variability and inherently small sample sizes. The exception was the end-point CO 2 relative permeability, k$0\atop{r}$ CO2 , the maximum attainable effective CO 2 permeability during CO 2 invasion, which changed CO2 injectivity by as much as 80%. Similarly, Monte Carlo simulation using 1000 realizations of relative permeability parameters showed no relationship between CO 2 injectivity and any of the parameters but k$0\atop{r}$ CO2 , which had a very strong (R 2 = 0.9685) power law relationship with total CO 2 injected. Model sensitivity to k$0\atop{r}$ CO2 points to the importance of accurate core flood and wettability measurements.

58 GEOSCIENCES↗
DOE OSTI2018

Impacts of hydrogeological characteristics on groundwater-level changes induced by earthquakes

Changes in groundwater level during earthquakes have been reported worldwide. In this study, field observations of co-seismic groundwater-level changes in wells under different aquifer conditions and sampling intervals due to near-field earthquake events in Taiwan are presented. Sustained changes, usually observed immediately after earthquakes, are found in the confined aquifer. Oscillatory changes due to the dynamic strain triggered by passing earthquake waves can only be recorded by a high-frequency data logger. While co-seismic changes recover rapidly in an unconfined aquifer, they can sustain for months or longer in a confined aquifer. Three monitoring wells with long-term groundwater-level data were examined to understand the association of co-seismic changes with local hydrogeological conditions. The finite element software ABAQUS is used to simulate the pore-pressure changes induced by the displacements due to fault rupture. The calculated co-seismic change in pore pressure is related to the compressibility of the formation. The recovery rate of the change is rapid in the unconfined aquifer due to the hydrostatic condition at the water table, but slow in the confined aquifer due to the less permeable confining layer. Fracturing of the confining layer during earthquakes may enhance the dissipation of pore pressure and induce the discharge of the confined aquifer. The study results indicated that aquifer characteristics play an important role in determining groundwater-level changes during and after earthquakes.

58 GEOSCIENCES↗
DOE OSTI2020

Non-destructive simulation of node defects in additively manufactured lattice structures

Additive Manufacturing (AM), commonly referred to as 3D printing, offers the ability to not only fabricate geometrically complex lattice structures but parts in which lattice topologies in-fill volumes bounded by complex surface geometries. However, current AM processes produce defects on the strut and node elements which make up the lattice structure. This creates an inherent difference between the as-designed and as-fabricated geometries, which negatively affects predictions (via numerical simulation) of the lattice’s mechanical performance. Although experimental and numerical analysis of an AM lattice’s bulk structure, unit cell and struts have been performed, there exists almost no research data on the mechanical response of the individual as-manufactured lattice node elements. Here we propose a methodology that, for the first time, allows non-destructive quantification of the mechanical response of node elements within an as-manufactured lattice structure. A custom-developed tool is used to extract and classify each individual node geometry from micro-computed tomography scans of an AM fabricated lattice. Voxel-based finite element meshes are generated for numerical simulation and the mechanical response distribution is compared to that of the idealised computer-aided design model. The method demonstrates compatibility with Uncertainty Quantification methods that provide opportunities for efficient prediction of a population of nodal responses from sampled data. Overall, the non-destructive and automated nature of the node extraction and response evaluation is promising for its application in qualification and certification of additively manufactured lattice structures.

36 MATERIALS SCIENCE↗