Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Modelling”

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,099 records · Page 61

Synthesis, Microstructure and Properties of Nickel Aluminide Foams

Two Ph.D. students were involved in the project: Mr. Christopher Schuh (part-time, graduated in Spring 2001) and Ms. Andrea Hodge (full-time, graduated Summer 2002). One post-doctoral fellow, Dr. Heeman Choe, worked full-time on the project from July to December 2002. A new process to aluminize and chromize nickel foams was created. A kinetic aluminization model was developed. Creep testing was conducted on the foams. A finite-element model and a simplified analytical model for foam creep were produced. Four articles were written: one is published, two are accepted for publication, and one is in preparation. Ms. Hodge spent four months at NASA Glenn Research Center (9-12/2001 and 2-3/2002) under the supervision of Dr. Nathal. She conducted research on NiAl foam fabrication, mechanical testing and numerical modeling. She gave a talk at the ASM annual conference in November 2001 and presented her results at NASA in December 2001.

Dunand, David C.↗

Structural Dynamics Modeling of HIRENASD in Support of the Aeroelastic Prediction Workshop

An Aeroelastic Prediction Workshop (AePW) was held in April 2012 using three aeroelasticity case study wind tunnel tests for assessing the capabilities of various codes in making aeroelasticity predictions. One of these case studies was known as the HIRENASD model that was tested in the European Transonic Wind Tunnel (ETW). This paper summarizes the development of a standardized enhanced analytical HIRENASD structural model for use in the AePW effort. The modifications to the HIRENASD finite element model were validated by comparing modal frequencies, evaluating modal assurance criteria, comparing leading edge, trailing edge and twist of the wing with experiment and by performing steady and unsteady CFD analyses for one of the test conditions on the same grid, and identical processing of results.

Wieseman, Carol↗

Lumped Parameter Models of the Central Nervous System for VIIP Research

INTRODUCTION: Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit, such as to Mars and asteroids, expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome [1]. It has been hypothesized that the headward shift of cerebral spinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn induces VIIP syndrome through biomechanical pathways [1, 2]. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the realted IWS abstracts submitted by Nelson et al., Feola et al. and Ethier et al. METHODS: We have developed a nine compartment CNS model (Figure 1) capable of both time-dependent and steady state fluid transport simulations, based on the works of Stevens et al. [3]. The breakdown of compartments within the model includes: vascular (3), CSF (2), brain (1) and extracranial (3). The boundary pressure in the Central Arteries [A] node is prescribed using an oscillating pressure function PA(t) simulating the carotid pulsatile pressure wave as developed by Linninger et al. [4]. For each time step, pressures are integrated through time using an adaptive-timestep 4th and 5th order Runga-Kutta solver. Once pressures are found, constitutive equations are used to solve for flowrates (Q) between each compartment. In addition to fluid flow between the different compartments, compliance (C) interactions between neighboring compartments are represented. We are also developing a second CNS model based on the works of Linninger et al. [4] which takes a more granular approach to represent the interactions of the intracranial and spinal compartments with the inclusion of arteries, arterioles, capillaries, venules, veins, venous sinus, and ventricles. The flow through the arteries, veins and CSF compartments are governed by continuity, momentum and distensibility balance equations. Furthermore, unlike the Stevens et al. approach, the Monro-Kellie doctrine of constant cranial volume and the bi-phasic nature of the brain parenchyma are implemented. These features appear to be more consistent with the physiologic and anatomical behavior of the CNS, and follow a modeling philosophy similar to the lumped parameter eye model that is intended to be integrated with the CNS model. However, Linninger’s approach has never been implemented to include hydrostatic gradient and microgravity simulation capabilities. Therefore, we aim at implement this modeling approach for spaceflight simulations and assess its overall applicability to VIIP research. OBJECTIVES: We will present verification and validation test results for both models, as well as head-to-head comparison to explore their strengths and limitations with respect to mathematical implementation and physiological significance for VIIP research. In doing so, we hope to provide some guidance to the HRP research community on how to appropriately leverage lumped parameter models for space biomedical research.

Vera, J.↗

VIIP: Central Nervous System (CNS) Modeling

Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome. It has been hypothesized that the headward shift of cerebrospinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn may then induce VIIP syndrome through interaction with various biomechanical pathways. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the related IWS works submitted by Nelson et al., Feola et al. and Ethier et al.

intraocular pressure↗

Influence of Population Variation of Physiological Parameters in Computational Models of Space Physiology

The earliest manifestations of Visual Impairment and Intracranial Pressure (VIIP) syndrome become evident after months of spaceflight and include a variety of ophthalmic changes, including posterior globe flattening and distension of the optic nerve sheath. Prevailing evidence links the occurrence of VIIP to the cephalic fluid shift induced by microgravity and the subsequent pressure changes around the optic nerve and eye. Deducing the etiology of VIIP is challenging due to the wide range of physiological parameters that may be influenced by spaceflight and are required to address a realistic spectrum of physiological responses. Here, we report on the application of an efficient approach to interrogating physiological parameter space through computational modeling. Specifically, we assess the influence of uncertainty in input parameters for two models of VIIP syndrome: a lumped-parameter model (LPM) of the cardiovascular and central nervous systems, and a finite-element model (FEM) of the posterior eye, optic nerve head (ONH) and optic nerve sheath. Methods: To investigate the parameter space in each model, we employed Latin hypercube sampling partial rank correlation coefficient (LHSPRCC) strategies. LHS techniques outperform Monte Carlo approaches by enforcing efficient sampling across the entire range of all parameters. The PRCC method estimates the sensitivity of model outputs to these parameters while adjusting for the linear effects of all other inputs. The LPM analysis addressed uncertainties in 42 physiological parameters, such as initial compartmental volume and nominal compartment percentage of total cardiac output in the supine state, while the FEM evaluated the effects on biomechanical strain from uncertainties in 23 material and pressure parameters for the ocular anatomy. Results and Conclusion: The LPM analysis identified several key factors including high sensitivity to the initial fluid distribution. The FEM study found that intraocular pressure and intracranial pressure had dominant impact on the peak strains in the ONH and retro-laminar optic nerve, respectively; optic nerve and lamina cribrosa stiffness were also important. This investigation illustrates the ability of LHSPRCC to identify the most influential physiological parameters, which must therefore be well-characterized to produce the most accurate numerical results.

biomechanics↗

Comparison of AGR-3/4 Fission Product Transport Model to Measurements and Extraction of Diffusivities via Analytical Fits

A one-dimensional finite-element model based on the MOOSE framework has been constructed, and overall transport estimates are predicted based on available diffusivities and sorption isotherms in the extant literature. The predicted transport profiles are compared to experimental results, where such results are available, including total inventory analysis of gamma emitters, liquid scintillation results for Sr-90 concentrations, radial deconsolidation leach-burn-leach within compacts, tomographic reconstructions based on detailed gamma scanning and destructive physical sampling of several rings. The magnitude of expected transport to the sink ring is generally predicted by the model for silver and cesium, however, modeled detailed concentration profiles are not consistent with experimental results. Strontium and europium transport is not well characterized by model results.

07 ISOTOPE AND RADIATION SOURCES↗

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗

Second order tensor finite element

The results of a research and software development effort are presented for the finite element modeling of the static and dynamic behavior of anisotropic materials, with emphasis on single crystal alloys. Various versions of two dimensional and three dimensional hybrid finite elements were implemented and compared with displacement-based elements. Both static and dynamic cases are considered. The hybrid elements developed in the project were incorporated into the SPAR finite element code. In an extension of the first phase of the project, optimization of experimental tests for anisotropic materials was addressed. In particular, the problem of calculating material properties from tensile tests and of calculating stresses from strain measurements were considered. For both cases, numerical procedures and software for the optimization of strain gauge and material axes orientation were developed.

Oden, J. Tinsley↗

Applications of Displacement Transfer Functions to Deformed Shape Predictions of the G-III Swept-Wing Structure

In support of the Adaptive Compliant Trailing Edge [ACTE] project at the NASA Armstrong Flight Research Center, displacement transfer functions were applied to the swept wing of a Gulfstream G-III airplane (Gulfstream Aerospace Corporation, Savannah, Georgia) to obtain deformed shape predictions. Four strainsensing lines (two on the lower surface, two on the upper surface) were used to calculate the deformed shape of the G III wing under bending and torsion. There being an insufficient number of surface strain sensors, the existing G III wing box finite element model was used to generate simulated surface strains for input to the displacement transfer functions. The resulting predicted deflections have good correlation with the finite-element generated deflections as well as the measured deflections from the ground load calibration test. The convergence study showed that the displacement prediction error at the G III wing tip can be reduced by increasing the number of strain stations (for each strain-sensing line) down to a minimum error of l.6 percent at 17 strain stations; using more than 17 strain stations yielded no benefit because the error slightly increased to 1.9% when 32 strain stations were used.

displacement theory↗

THERMAL STRESS ANALYSIS OF A SPENT NUCLEAR FUEL CANISTER

The potential for chloride induced stress corrosion cracking (CISCC) in spent nuclear fuel dry storage canisters is a current topic of research by the US Department of Energy Spent Fuel and Waste Science and Technology program. One of the important prerequisites for CISCC is a tensile residual stress state. This study utilizes computational models to provide an initial analysis of thermal stresses in a generic spent nuclear fuel canister. A STAR-CCM+ thermal fluid model provides a temperature profile analysis of the canister based on four different ambient temperature conditions. An ANSYS APDL finite element model incorporates the temperature profiles to analyze the thermal stresses in the canister. The calculated thermal stress magnitudes are in the range of 10 MPa to 80 MPa, which could be significant for crack propagation through the canister wall. The maximum thermal stress is not a concern for the structural failure of the canister, but it is high enough that it is a significant contribution to the total stress state of the canister wall, which also includes residual stress from fabrication, internal pressure from helium cover gas, and potential transient mechanical loads, such as earthquakes. This paper describes the initial finite element thermal stress analysis that was completed in 2020, reports the initial findings, and identifies areas where model refinement is still needed to complete this work in the future.

Jensen, Ben J.↗

Moderate Velocity Ball Impact of a Mock High-Explosive

Modeling of thermal and mechanical events in high-explosive materials is complicated by the composite nature of the material, which experiences viscoelastic and plastic deformations and sustains damage in the form of microcracks that can dominate its overall behavior. A mechanical event of interest is projectile interaction with the material, which leads to extreme local deformation and adiabatic heating, which can potentially lead to adverse outcomes in an energetic material. Simulations of such an event predicted large local temperature rises near the path of a spherical projectile, but these were experimentally unconfirmed and hence potentially non-physical. This work concerns the experimental verification of local temperatures both at the surface and in the wake of a spherical projectile penetrating a mock (unreactive) high-explosive at {approx}700 m/s. Fast response thermocouples were embedded radially in a mid-plane of a cylindrical target, which was bonded around the thermocouples with epoxy and recorded by an oscilloscope through a low-pass filter with a bandwidth of 500 Hz. A peak temperature rise of 70 K was measured both at the equator of the projectile and in its wake, in good agreement with the temperature predicted in the minimally distorted elements at those locations by a finite element model in ABAQUS employing the ViscoSCRAM constitutive model. Further work is needed to elucidate the extreme temperature rises in material undergoing crushing or fragmentation, which is difficult to predict with meshed finite element methods due to element distortion, and also challenging to quantify experimentally.

36 MATERIALS SCIENCE↗

Numerical Manufacturing And Design Tool (NuMAD v2.0) for Wind Turbine Blades: User’s Guide

Sandia National Laboratories has an on-going effort to reduce the cost of energy and improve reliability for wind systems through improved blade design and manufacture. As part of this effort, a software tool named NuMAD (Numerical Manufacturing And Design) has been developed to greatly simplify the process of creating a three-dimensional finite element model for a modern wind turbine blade. NuMAD manages all blade information including databases of airfoils, materials, and material placement to enable efficient creation of models. NuMAD is a stand-alone, user-friendly, graphical pre-processor for the ANSYS® commercial finite element package. The blade information contained in the NuMAD database is also used to manage capabilities such as output for CFD mesh, computation of blade cross section properties, and aeroelastic instability analysis of the blade. This user’s manual describes the capabilities and usage of NuMAD.

Berg, Jonathan Charles↗

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗

Structural Model Tuning Capability in an Object-Oriented Multidisciplinary Design, Analysis, and Optimization Tool

Updating the finite element model using measured data is a challenging problem in the area of structural dynamics. The model updating process requires not only satisfactory correlations between analytical and experimental results, but also the retention of dynamic properties of structures. Accurate rigid body dynamics are important for flight control system design and aeroelastic trim analysis. Minimizing the difference between analytical and experimental results is a type of optimization problem. In this research, a multidisciplinary design, analysis, and optimization [MDAO] tool is introduced to optimize the objective function and constraints such that the mass properties, the natural frequencies, and the mode shapes are matched to the target data as well as the mass matrix being orthogonalized.

Lung, Shun-fat↗

Structural Model Tuning Capability in an Object-Oriented Multidisciplinary Design, Analysis, and Optimization Tool

Updating the finite element model using measured data is a challenging problem in the area of structural dynamics. The model updating process requires not only satisfactory correlations between analytical and experimental results, but also the retention of dynamic properties of structures. Accurate rigid body dynamics are important for flight control system design and aeroelastic trim analysis. Minimizing the difference between analytical and experimental results is a type of optimization problem. In this research, a multidisciplinary design, analysis, and optimization (MDAO) tool is introduced to optimize the objective function and constraints such that the mass properties, the natural frequencies, and the mode shapes are matched to the target data as well as the mass matrix being orthogonalized.

Lung, Shun-fat↗

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗

Finite Rotation Analysis of Highly Thin and Flexible Structures

Deployable space structures such as sunshields and solar sails are extremely thin and highly flexible with limited bending rigidity. For analytical investigation of their responses during deployment and operation in space, these structures can be modeled as thin shells. The present work examines the applicability of the solid shell element formulation to modeling of deployable space structures. The solid shell element formulation that models a shell as a three-dimensional solid is convenient in that no rotational parameters are needed for the description of kinematics of deformation. However, shell elements may suffer from element locking as the thickness becomes smaller unless special care is taken. It is shown that, when combined with the assumed strain formulation, the solid shell element formulation results in finite element models that are free of locking even for extremely thin structures. Accordingly, they can be used for analysis of highly flexible space structures undergoing geometrically nonlinear finite rotations.

Clarke, Greg V.↗

A Modeling Approach for Time-Dependent Geometry Applied to Transient Heat Transfer of Aluminum Electrolysis Cells

The thermal balance of aluminum electrolysis cells (AEC) have to be rigorously controlled in order to improve the efficiency and sustainability of this industrial process. A new modeling strategy is developed to consider the displacements of solid bodies and moving boundaries in finite element models. The transient thermal-electric modeling of the AEC demonstrates the effect of an increase in operating voltage on both the anode cover and the side ledge. With higher heat generation, the anode cover deteriorates and the side ledge thickness decreases. Since the anode cover is characterized by irreversible transformations, the top heat dissipation remains higher even when the operating voltage comes back to its typical value. For the first time, the transient temperature and electric fields throughout the anode life are simulated and validated by industrial measurements. The modeling predictions have been validated from instrumented anodes and manual measurements, all performed on operating AEC.

30 DIRECT ENERGY CONVERSION↗