Engineering PapersSearch

SEARCH · Engineering Papers

Results for “FEM”

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 55 records · Page 3

Algebraic Nonoverlapping Domain Decomposition Methods for Stabilized FEM and FV Discretizations

We consider preconditioning methods for convection dominated fluid flow problems based on a nonoverlapping Schur complement domain decomposition procedure for arbitrary triangulated domains. The triangulation is first partitioned into a number of subdomains and interfaces which induce a natural 2 x 2 partitioning of the p.d.e. discretization matrix. We view the Schur complement induced by this partitioning as an algebraically derived coarse space approximation. This avoids the known difficulties associated with the direct formation of an effective coarse discretization for advection dominated equations. By considering various approximations of the block factorization of the 2 x 2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 using MPI message passing protocol. A number of 2-D CFD calculations will be presented for both scalar advection-diffusion equations and the Euler equations discretized using stabilized finite element and finite volume methods. These results show very good scalability of the preconditioner for various discretizations as the number of processors is increased while the number of degrees of freedom per processor is fixed.

Barth, Timothy J.

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy

Eigenvalues of Rectangular Waveguide Using FEM With Hybrid Elements

A finite element analysis using hybrid triangular-rectangular elements is developed to estimate eigenvalues of a rectangular waveguide. Use of rectangular vector-edge finite elements in the vicinity of the PEC boundary and triangular elements in the interior region more accurately models the physical nature of the electromagnetic field, and consequently quicken the convergence.

Deshpande, Manohar D.

Inverse FEM for Full-Field Reconstruction of Elastic Deformations in Shear Deformable Plates and Shells

The inverse problem of real-time reconstruction of full-field structural displacements is addressed through the application of a new variational formulation leading to versatile, robust and computationally efficient inverse shell finite element analysis. Utilizing surface strain measurements from strain sensors mounted on the load-carrying structural components, the methodology enables accurate computations of the three-dimensional displacement field. This high fidelity computational technology is essential for providing feedback to the actuation and control systems of the next generation of aerospace vehicles.

Tessler, Alexander

Flow Separation Side Loads Excitation of Rocket Nozzle FEM

Modern rocket nozzles are designed to operate over a wide range of altitudes, and are also built with large aspect ratios to enable high efficiencies. Nozzles designed to operate over specific regions of a trajectory are being replaced in modern launch vehicles by those that are designed to operate from earth to orbit. This is happening in parallel with modern manufacturing and wall cooling techniques allowing for larger aspect ratio nozzles to be produced. Such nozzles, though operating over a large range of altitudes and ambient pressures, are typically designed for one specific altitude. Above that altitude the nozzle flow is 'underexpanded' and below that altitude, the nozzle flow is 'overexpanded'. In both conditions the nozzle produces less than the maximum possible thrust at that altitude. Usually the nozzle design altitude is well above sea level, leaving the nozzle flow in an overexpanded state for its start up as well as for its ground testing where, if it is a reusable nozzle such as the Space Shuttle Main Engine (SSME), the nozzle will operate for the majority of its life. Overexpansion in a rocket nozzle presents the critical, and sometimes design driving, problem of flow separation induced side loads. To increase their understanding of nozzle side loads, engineers at MSFC began an investigation in 2000 into the phenomenon through a task entitled "Characterization and Accurate Modeling of Rocket Engine Nozzle Side Loads", led by A. Brown. The stated objective of this study was to develop a methodology to accurately predict the character and magnitude of nozzle side loads. The study included further hot-fire testing of the MC-l engine, cold flow testing of subscale nozzles, CFD analyses of both hot-fire and cold flow nozzle testing, and finite element (fe.) analysis of the MC-1 engine and cold flow tested nozzles. A follow on task included an effort to formulate a simplified methodology for modeling a side load during a two nodal diameter fluid/structure interaction for a single moment in time.

Smalley, Kurt B.

Modal Correlation of Complex Aerospace Joints Using Automated Variable Substitution

A critical task involved with being able to predict flight loads accurately in aerospace finite element models (FEMs) is the prior verification of the FEMs by conducting modal survey testing (MST). Experience comparing dynamic response of initial FEMs to MST data tends to demonstrate that FEMs can have unacceptable accuracy even when best modeling practices are followed. One inherent source of inaccuracy in linear dynamic FEMs is the modeling of nonlinear joints with mechanisms such as spherical bearings. These joints are usually designed to freely translate or rotate under the high levels of loading experienced in flight. Engineers who create linear FEMs conventionally model these joints without any stiffness in the mechanism degrees of freedom to meet this design intent. However, inaccuracy is observed during test validation of these FEMs, which usually relies on low-level force excitation orders of magnitude below flight load levels. This low-level modal test rarely overcomes the joint friction that is present, and thus the mechanism joints are able to react loads. This divide between the test results and the FEM creates a significant challenge to the engineer who is performing the correlation, in that the engineer has no basis for what stiffness value should be used to make the FEM match the test results. A compounding challenge is that complicated built-up aerospace structures commonly have multiple joints through a load path where each joint will “stick” and “slip” at different levels of force input. Explicitly matching the dynamics of a system containing these nonlinear mechanisms would require a nonlinear FEM, which is prohibitively costly for dynamic simulations of most aerospace systems. The objective of this paper is to present a workflow that can efficiently cycle through many iterations of a FEM, allowing a Monte Carlo style examination of the design space to identify candidate stiffness values for nonlinear mechanism joints. The outlined approach is specific to MSC Nastran and utilizes MSC Nastran’s symbolic substitution capabilities, coupled with the IMAT™ and Attune™ software packages developed by ATA Engineering, Inc. (ATA). The workflow is demonstrated with a case study from the correlation effort for The Boeing Company’s Crew Space Transportation (CST)-100 Starliner FEM. Keywords: Correlation, MSC Nastran, IMAT™, Attune™, Finite Element Model, Modal Testing

Correlation

Finite elements: Theory and application

Recent advances in FEM techniques and applications are discussed in reviews and reports presented at the ICASE/LaRC workshop held in Hampton, VA in July 1986. Topics addressed include FEM approaches for partial differential equations, mixed FEMs, singular FEMs, FEMs for hyperbolic systems, iterative methods for elliptic finite-element equations on general meshes, mathematical aspects of FEMS for incompressible viscous flows, and gradient weighted moving finite elements in two dimensions. Consideration is given to adaptive flux-corrected FEM transport techniques for CFD, mixed and singular finite elements and the field BEM, p and h-p versions of the FEM, transient analysis methods in computational dynamics, and FEMs for integrated flow/thermal/structural analysis.

Dwoyer, D. L.

TPSAS-NF1676L-33512-DND

INTRODUCTION: In an effort to quantify the predictive capability of the Livermore Software Technology Corporation (LSTC)-provided Hybrid III (small female, mid-size male, and large male) anthropomorphic test device (ATD) finite element models (FEMs) within aerospace loading environments, the National Aeronautics and Space Administration (NASA) performed a series of sled tests with the isolated head-neck complex of each ATD. These tests and subsequent modeling efforts revealed significant deficiencies in the large male head neck FEM. Differences in both the shape and material properties between the FEM and the physical ATD were identified. In order to improve the predictive response of the FEM, a model development effort was undertaken to update the head-neck complex to better represent the geometric and material properties of the physical ATD.METHODS: To better match the physical geometry, a new neck component model was meshed and implemented into the Hybrid III large male head-neck complex. The updated geometry and weight were verified against measurements of the physical ATD to insure accuracy. The primary source of material discrepancy between the physical ATD and FEM was found to be in the material definition of the rubber neck pucks. Using LS-Opt the material properties of this part were calibrated to match physical response of the head-neck complex in the isolated head-neck testing performed. RESULTS AND DISCUSSION: The updated geometry and calibrated material properties were shown to significantly improve the FEM predictive accuracy under lateral, rear, and frontal impacts with combined horizontal and vertical loading. The results of this study demonstrate an effective means for improving ATD FEM response through isolated component level testing and material parameter optimization. The improved and quantified accuracy of the Hybrid III large male FEM head-neck complex stemming from with work lends to confidence in its future use for occupant protection evaluation of large male occupants

J Putnam

Development of Test-Analysis Models (TAM) for correlation of dynamic test and analysis results

The primary objective of structural analysis of aerospace applications is to obtain a verified finite element model (FEM). The verified FEM can be used for loads analysis, evaluate structural modifications, or design control systems. Verification of the FEM is generally obtained as the result of correlating test and FEM models. A test analysis model (TAM) is very useful in the correlation process. A TAM is essentially a FEM reduced to the size of the test model, which attempts to preserve the dynamic characteristics of the original FEM in the analysis range of interest. Numerous methods for generating TAMs have been developed in the literature. The major emphasis of this paper is a description of the procedures necessary for creation of the TAM and the correlation of the reduced models with the FEM or the test results. Herein, three methods are discussed, namely Guyan, Improved Reduced System (IRS), and Hybrid. Also included are the procedures for performing these analyses using MSC/NASTRAN. Finally, application of the TAM process is demonstrated with an experimental test configuration of a ten bay cantilevered truss structure.

Angelucci, Filippo

Test and Simulation of a Fokker F28 Crash Landing

Test and Simulation of a Fokker F28 Crash Landing – Jacob B. Putnam 1, Justin D. Littell 1, Karen E. Jackson 2; 1 NASA Langley Research Center, Hampton, VA, 2 National Institute of Aerospace, Hampton, VA. In June of 2019, the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) conducted a full-scale crash test of a Fokker F28 MK1000 aircraft. Aircraft crash-worthiness is typically evaluated through component level tests (i.e. vertical drops of fuselage subsections or isolated seat tests). Finite Element Models (FEMs) are also used to bridge the gap between component testing and full-scale crash prediction. The full-scale crash test performed in conjunction with previous subsection testing of a Fokker F28 fuselage provided the opportunity to evaluate differences in crash-worthiness predictions between full- and sub-scale testing and to quantify the predictive capability of FEMs in the aerospace crash environment. In this study vehicle crash-worthiness was quantified through anthropometric test devices (ATDs) included in both the full-scale crash test as well as fuselage section drop tests previously performed at NASA LaRC. A FEM of the Fokker F28 aircraft was developed and simulated within the tested environment. Vehicle and ATD response predictions were compared between test and simulation. The International Organization for Standardization ISO/TR 16250 curve comparison methodology was used to provide a quantitative assessment of predictive accuracy for both the vehicle and ATD FEMs. Determination of aircraft crash-worthiness made between the three evaluations methodologies (full-vehicle crash test, component test, and FEM simulation) were then compared. Structural differences as well as the more complex loading environment achieved within the full-vehicle test resulted in increased injury risk compared to that predicted within the component level tests. The FEM simulations were found to produce a more realistic prediction of injury risk within the full-vehicle crash environment than the component level testing. Limitations of both component testing and FEM simulation within aerospace crash environment were identified.

impact testing

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification

Structural Sizing of a Composite Transonic Truss-Braced Wing

Accurate finite element modeling (FEM) is a vital part of the modern aircraft design process. As aircraft become increasingly complex, the time-consuming nature of detailed FEM approaches comes at a significant cost to program timeline and budget. In an effort produce modeling efforts that are sufficiently accurate and minimally costly, a proposed FEM approach and optimization scheme for a composite aircraft is explored, including a discussion of the manufacturing constraints of a highly tailored composite panel design. A high-fidelity structural model of a transonic truss-braced wing (TTBW) is generated, and the components are sized by structural optimization to satisfy buckling and strength constraints while subjected to critical maneuver loads. The structural modeling approaches and sizing of a TTBW are discussed, including details for FEM approaches, verification of an approximated FEM approach, a sizing optimization using the optimization software LS-OPT, and a manufacturing trial of integrally stiffened composite panels conducted to explore the validity of highly tailored composites as a design consideration. The results of study discussed herein indicate that the proposed FEM approach is suitable for modeling composite-construction aircraft and for use in sizing optimization. Further efforts regarding integrally stiffened composite panels will indicate the suitability of this method for the integration of highly tailored composite panels into the design and optimization process, given the manufacturability of such panels as shown here.

Finite Element

Structural Sizing of a Transonic Truss-Braced Wing

Accurate finite element modeling (FEM) is a vital part of the modern aircraft design process. As aircraft become increasingly complex, the time-consuming nature of detailed FEM approaches comes at a significant cost to program timeline and budget. In an effort produce modeling efforts that are sufficiently accurate and minimally costly, a proposed FEM approach and optimization scheme for a composite aircraft is explored, including a discussion of the manufacturing constraints of a highly tailored composite panel design. A high-fidelity structural model of a transonic truss-braced wing (TTBW) is generated, and the components are sized by structural optimization to satisfy buckling and strength constraints while subjected to critical maneuver loads. The structural modeling approaches and sizing of a TTBW are discussed, including details for FEM approaches, verification of an approximated FEM approach, a sizing optimization using the optimization software LS-OPT, and a manufacturing trial of integrally stiffened composite panels conducted to explore the validity of highly tailored composites as a design consideration. The results of study discussed herein indicate that the proposed FEM approach is suitable for modeling composite-construction aircraft and for use in sizing optimization. Further efforts regarding integrally stiffened composite panels will indicate the suitability of this method for the integration of highly tailored composite panels into the design and optimization process, given the manufacturability of such panels as shown here.

Finite Element

Solving sparse finite element problems on neuromorphic hardware

The finite element method (FEM) is one of the most important and ubiquitous numerical methods for solving partial differential equations (PDEs) on computers for scientific and engineering discovery. Applying the FEM to larger and more detailed scientific models has driven advances in high-performance computing for decades. Here we demonstrate that scalable spiking neuromorphic hardware can directly implement the FEM by constructing a spiking neural network that solves the large, sparse, linear systems of equations at the core of the FEM. We show that for the Poisson equation, a fundamental PDE in science and engineering, our neural circuit achieves meaningful levels of numerical accuracy and close to ideal scaling on modern, inherently parallel and energy-efficient neuromorphic hardware, specifically Intel’s Loihi 2 neuromorphic platform. We illustrate extensions to irregular mesh geometries in both two and three dimensions as well as other PDEs such as linear elasticity. Our spiking neural network is constructed from a recurrent network model of the brain’s motor cortex and, in contrast to black-box deep artificial neural network-based methods for PDEs, directly translates the well-understood and trusted mathematics of the FEM to a natively spiking neuromorphic algorithm.

Applied mathematics