Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element models”

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 217 records · Page 12

Finite element modeling of stress in the Nazca plate - Driving forces and plate boundary earthquakes

The state of stress within the Nazca plate due to plate driving forces and large plate boundary earthquakes has been analyzed by applying a finite element method using the wave front solution technique to models of the intraplate stress field in a single plate using a refined grid. Although only static elastic models have been explicitly calculated, certain limiting cases of an elastic plate over a viscous asthenosphere were also treated. A state of nearly east-west compression inferred from the source mechanism of thrust earthquakes in the interior of the plate requires ridge pushing forces. The net pulling force on the oceanic plate by the subducted slab has a maximum value comparable to pushing forces. The estimated horizontal deviatoric stress in intraplate regions, based on potential forces associated with the ridge, is on the order of a few hundred bars. The intraplate stress field in the region of the 1960 earthquake may change by a few tens of bars at most once the asthenosphere has relaxed, with changes on the order of one bar occurring at greater distances into the plate. The changes in the intraplate stress field are probably not noticeable unless the lithosphere is near failure.

Richardson, R. M.↗

A finite element model with nonviscous damping

A constitutive law by which structural damping is modeled as a relationship between stress, strain, and strain rate in a material is used in conjunction with the finite element method to develop general integral expressions for viscous and nonviscous damping matrices. To solve the set of nonlinear equations resulting from the presence of nonviscous damping, a solution technique is developed by modifying the Newmark method to accommodate an iterative solution and treat the nonviscous damping as a pseudo-force. The technique is then checked for accuracy and convergence in single- and multi-degree-of-freedom problems, and is found to be accurate and efficient for initial-condition problems with small nonviscous damping.

Roussos, L. A.↗

Material Property Estimation in Thin Battery Components Using Guided Wave Measurement, Experimental Dispersion Curve Extraction and Finite Element Modeling

At NASA, we are investigating nondestructive evaluation (NDE) and structural health monitoring (SHM) techniques to detect precursors of thermal runaway failure in lithium metal based lithium ion batteries (LIB). The approach is centered on computational simulation models to guide inspection and aid in interpretation of results. Since lithium metal LIBs have combinations of solid and fluid-filled porous materials, obtaining accurate material properties is both challenging and critical for successful simulation of battery inspection. To this end, we investigated a multi-verification approach for material characterization of thin battery components. First, a laser Doppler vibrometer (LDV) was used to measure guided wave fields in thin (microns-thick) battery components subject to broadband excitation. The time-space wavefield data was converted to frequency-wave number data to extract guided wave dispersion curves. A data visualization and post processing graphics user interface (GUI) was developed at NASA to aid the data exploration and analysis. Due to the thinness of the samples, low frequency-thickness-product plate wave approximations allowed for the calculation of closed-form solutions for material elastic property estimation. These approximations were then verified by calculating the Lamb wave solutions using the previously obtained material properties. Finally, the estimated elastic material properties were implemented in a COMSOL simulation model, and dispersion curves were extracted from simulation results. The dispersion curves and derived material properties were then compared to the analysis results from the experimental data. These comparisons informed on the accuracy of the measured material properties and helped demonstrate the accuracy of the finite element analysis (FEA) computational models. This assessment will prove vital when we start simulating more complex multi-layer components and poroelastic media. This paper gives a brief background of the problem space, outlines the workflow for data analysis and verification, shows results from the workflow, and gives an overview of future plans for simulation of lithium metal LIB inspection.

Peter Juarez↗

Finite element modelling of buried structures

In many structures the final stress states are dependent on the sequence of construction or the stress states at various stages of construction are of interest. Such problems can be analyzed using finite element programs that have the capability of adding (birthing) elements to simulate the progress of construction. However, the usual procedure of assembling elements may lead to numerical instabilities or stress states that are unrealistic. Both problems are demonstrated in the analysis of a structure using the program ADINA. A technique which combines application of a preload with element birthing to overcome these problems is described and illustrated.

Playdon, D. K.↗

Finite element models of non-Newtonian crater relaxation

The effects of a non-Newtonian rheology on the profiles of relaxing craters (such as those seen on the surfaces of the icy Galilean and Saturnian satellites) were studied. Two-dimensional finite element simulations of non-Newtonian viscous flow were performed, and the results were compared with those associated with Newtonian rheology. Viscous relaxation of craters in a non-Newtonian medium was significantly different from that in a Newtonian medium. Crater rims are observed to relax at a more rapid rate in a non-Newtonian region as a result of the movement of the low viscosity region to underneath the crater rim after the initial relaxation of the bowl. Significant differences are also found when central depth is plotted as a function of time. For a Newtonian medium, crater relaxation is exponential in form. In contrast, non-Newtonian crater relaxation is initially rapid, in response to the large initial stresses and small viscosities; however, as stresses decrease, this relaxation becomes extremely gradual.

Thomas, Paul J.↗

Interlaminar fracture toughness: Three-dimensional finite element modeling for end-notch and mixed-mode flexure

A computational procedure is described for evaluating End-Notch-Flexure (ENF) and Mixed-Mode-Flexure (MMF) interlaminar fracture toughness in unidirectional fiber composites. The procedure consists of a three-dimensional finite element analysis in conjunction with the strain energy release rate concept and with composite micromechanics. The procedure is used to analyze select cases of ENF and MMF. The strain energy release rate predicted by this procedure is in good agreement with limited experimental data. The procedure is used to identify significant parameters associated with interlaminar fracture toughness. It is also used to determine the critical strain energy release rate and its attendant crack length in ENF and/or MMF. This computational procedure has considerable versatility/generality and provides extensive information about interlaminar fracture toughness in fiber composites.

Murthy, P. L. N.↗

Determination of efficiencies, loss mechanisms, and performance degradation factors in chopper controlled dc vehical motors. Section 2: The time dependent finite element modeling of the electromagnetic field in electrical machines: Methods and applications

The time dependent solution of the magnetic field is introduced as a method for accounting for the variation, in time, of the machine parameters in predicting and analyzing the performance of the electrical machines. The method of time dependent finite element was used in combination with an also time dependent construction of a grid for the air gap region. The Maxwell stress tensor was used to calculate the airgap torque from the magnetic vector potential distribution. Incremental inductances were defined and calculated as functions of time, depending on eddy currents and saturation. The currents in all the machine circuits were calculated in the time domain based on these inductances, which were continuously updated. The method was applied to a chopper controlled DC series motor used for electric vehicle drive, and to a salient pole sychronous motor with damper bars. Simulation results were compared to experimentally obtained ones.

Hamilton, H. B.↗

Finite element modeling of truss structures with frequency-dependent material damping

A physically motivated modelling technique for structural dynamic analysis that accommodates frequency dependent material damping was developed. Key features of the technique are the introduction of augmenting thermodynamic fields (AFT) to interact with the usual mechanical displacement field, and the treatment of the resulting coupled governing equations using finite element analysis methods. The AFT method is fully compatible with current structural finite element analysis techniques. The method is demonstrated in the dynamic analysis of a 10-bay planar truss structure, a structure representative of those contemplated for use in future space systems.

Lesieutre, George A.↗

Finite Element Modeling Techniques for Analysis of VIIP

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a major health concern for long-duration space missions. Currently, it is thought that a cephalad fluid shift in microgravity causes elevated intracranial pressure (ICP) that is transmitted along the optic nerve sheath (ONS). We hypothesize that this in turn leads to alteration and remodeling of connective tissue in the posterior eye which impacts vision. Finite element (FE) analysis is a powerful tool for examining the effects of mechanical loads in complex geometries. Our goal is to build a FE analysis framework to understand the response of the lamina cribrosa and optic nerve head to elevations in ICP in VIIP.

eye (anatomy)↗

Finite element modeling of electromagnetic fields and waves using NASTRAN

The various formulations of Maxwell's equations are reviewed with emphasis on those formulations which most readily form analogies with Navier's equations. Analogies involving scalar and vector potentials and electric and magnetic field components are presented. Formulations allowing for media with dielectric and conducting properties are emphasized. It is demonstrated that many problems in electromagnetism can be solved using the NASTRAN finite element code. Several fundamental problems involving time harmonic solutions of Maxwell's equations with known analytic solutions are solved using NASTRAN to demonstrate convergence and mesh requirements. Mesh requirements are studied as a function of frequency, conductivity, and dielectric properties. Applications in both low frequency and high frequency are highlighted. The low frequency problems demonstrate the ability to solve problems involving media inhomogeneity and unbounded domains. The high frequency applications demonstrate the ability to handle problems with large boundary to wavelength ratios.

Moyer, E. Thomas, Jr.↗

Explicit Finite Element Modeling of Multilayer Composite Fabric for Gas Turbine Engine Containment Systems: Ballistic Impact Testing - Part 2

Under the Federal Aviation Administration's Airworthiness Assurance Center of Excellence and the Aircraft Catastrophic Failure Prevention Program, National Aeronautics and Space Administration Glenn Research Center collaborated with Arizona State University, Honeywell Engines, Systems and Services, and SRI International to develop improved computational models for designing fabric-based engine containment systems. In the study described in this report, ballistic impact tests were conducted on layered dry fabric rings to provide impact response data for calibrating and verifying the improved numerical models. This report provides data on projectile velocity, impact and residual energy, and fabric deformation for a number of different test conditions.

ZYLON↗

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗

Development of a Full-Scale Finite Element Model of the Fokker F28 Fellowship Aircraft and Crash Simulation Predictions

In June 2019, a full-scale crash test of a Fokker F28 Fellowship aircraft was conducted as part of a joint National Aeronautics and Space Administration/Federal Aviation Administration (NASA/FAA) project to investigate the performance of transport aircraft under realistic crash conditions. The test objectives were to provide data for assessment of transport aircraft crashworthiness under realistic impact conditions and to generate test data for model validation. The test article was loaded with transport aircraft seats in a 3+2 configuration. A total of 24 instrumented ATDs were placed in the seats and restrained. The test article weighed 33,306-lb. (15,107.3-kg) and, during the crash test, impacted a 2-ft. (0.61-m) high soil bed at 65.3-ft/s (19.9-m/s) forward and 31.8-ft/s (9.7-m/s) vertical velocity. The full-scale aircraft crash test was simulated using the commercial nonlinear explicit transient dynamic finite element code, LS-DYNA. This paper will provide a description of the test article and the crash test conditions, document the F28 full-scale model development, and present test-analysis comparisons. The test-analysis comparisons include inertial properties, kinematic responses, structural acceleration responses, and structural deformation and failure.

Karen E Jackson↗

Neotectonics of Asia: Thin-shell finite-element models with faults

As India pushed into and beneath the south margin of Asia in Cenozoic time, it added a great volume of crust, which may have been (1) emplaced locally beneath Tibet, (2) distributed as regional crustal thickening of Asia, (3) converted to mantle eclogite by high-pressure metamorphism, or (4) extruded eastward to increase the area of Asia. The amount of eastward extrusion is especially controversial: plane-stress computer models of finite strain in a continuum lithosphere show minimal escape, while laboratory and theoretical plane-strain models of finite strain in a faulted lithosphere show escape as the dominant mode. We suggest computing the present (or neo)tectonics by use of the known fault network and available data on fault activity, geodesy, and stress to select the best model. We apply a new thin-shell method which can represent a faulted lithosphere of realistic rheology on a sphere, and provided predictions of present velocities, fault slip rates, and stresses for various trial rheologies and boundary conditions. To minimize artificial boundaries, the models include all of Asia east of 40 deg E and span 100 deg on the globe. The primary unknowns are the friction coefficient of faults within Asia and the amounts of shear traction applied to Asia in the Himalayan and oceanic subduction zones at its margins. Data on Quaternary fault activity prove to be most useful in rating the models. Best results are obtained with a very low fault friction of 0.085. This major heterogeneity shows that unfaulted continum models cannot be expected to give accurate simulations of the orogeny. But, even with such weak faults, only a fraction of the internal deformation is expressed as fault slip; this means that rigid microplate models cannot represent the kinematics either. A universal feature of the better models is that eastern China and southeast Asia flow rapidly eastward with respect to Siberia. The rate of escape is very sensitive to the level of shear traction in the Pacific subduction zones, which is below 6 MPa. Because this flow occurs across a wide range of latitudes, the net eastward escape is greater than the rate of crustal addition in the Himalaya. The crustal budget is balanced by extension and thinning, primarily within the Tibetan plateau and the Baikal rift. The low level of deviation stresses in the best models suggests that topographic stress plays a major role in the orogeny; thus, we have to expect that different topography in the past may have been linked with fundamentally different modes of continental collision.

Kong, Xianghong↗