Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous finite element”

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 163 records · Page 9

Benchmark notch test for life prediction

The laser Interferometric Strain Displacement Gage (ISDG) was used to measure local strains in notched Inconel 718 test bars subjected to six different load histories at 649 C (1200 F) and including effects of tensile and compressive hold periods. The measurements were compared to simplified Neuber notch analysis predictions of notch root stress and strain. The actual strains incurred at the root of a discontinuity in cyclically loaded test samples subjected to inelastic deformation at high temperature where creep deformations readily occur were determined. The steady state cyclic, stress-strain response at the root of the discontinuity was analyzed. Flat, double notched uniaxially loaded fatigue specimens manufactured from the nickel base, superalloy Inconel 718 were used. The ISDG was used to obtain cycle by cycle recordings of notch root strain during continuous and hold time cycling at 649 C. Comparisons to Neuber and finite element model analyses were made. The results obtained provide a benchmark data set in high technology design where notch fatigue life is the predominant component service life limitation.

Domas, P. A.↗

Elastic/viscoplastic behavior of fiber-reinforced thermoplastic composites

An elastic/viscoplastic constitutive model was used to characterize the nonlinear and rate dependent behavior of a continuous fiber-reinforced thermoplastic composite. This model was incorporated into a finite element program for the analysis of laminated plates and shells. Details on the finite element formulation with the proposed constitutive model were presented. The numerical results were compared with experimental data for uniaxial tension and three-point bending tests of (+ or - 45 deg)3s APC-2 laminates.

Wang, C.↗

LS-DYNA Analysis of a Full-Scale Helicopter Crash Test

A full-scale crash test of an MD-500 helicopter was conducted in December 2009 at NASA Langley's Landing and Impact Research facility (LandIR). The MD-500 helicopter was fitted with a composite honeycomb Deployable Energy Absorber (DEA) and tested under vertical and horizontal impact velocities of 26 ft/sec and 40 ft/sec, respectively. The objectives of the test were to evaluate the performance of the DEA concept under realistic crash conditions and to generate test data for validation of a system integrated LS-DYNA finite element model. In preparation for the full-scale crash test, a series of sub-scale and MD-500 mass simulator tests was conducted to evaluate the impact performances of various components, including a new crush tube and the DEA blocks. Parameters defined within the system integrated finite element model were determined from these tests. The objective of this paper is to summarize the finite element models developed and analyses performed, beginning with pre-test and continuing through post test validation.

Annett, Martin S.↗

Improving Kinetic Consistency Across Elastic Model Transitions With Quadratic Inequality Constrained Weighted Least Squares (WLSQI)

Flexible body modeling presents a large challenge in the development of simulations to aid in design of flight control systems for launch and landing vehicles. Typically, the flexible body model is not a single continuous model but are rather discrete sets of Linear Time Invariant (LTI) Finite Element Models (FEM) incremented by propellant levels. This introduces the problem of smoothly transitioning modal and physical states of the vehicle when switching from one FEM to the next. This paper introduces a new approach to optimally transition flexible body states with weighted least squares, building off previous methods.

Convex Optimization↗

Improving Kinetic Consistency Across Elastic Model Transitions With Quadratic Inequality Constrained Weighted Least Squares (WLSQI)

Flexible body modeling presents a large challenge in the development of simulations to aid in design of flight control systems for launch and landing vehicles. Typically, the flexible body model is not a single continuous model but are rather discrete sets of Linear Time Invariant (LTI) Finite Element Models (FEM) incremented by propellant levels. This introduces the problem of smoothly transitioning modal and physical states of the vehicle when switching from one FEM to the next. This paper introduces a new approach to optimally transition flexible body states with weighted least squares, building off previous methods.

GN&C↗

R-Adaptivity to Enable Compression of Elementary Computations in Extreme-Scale Finite Element Simulators

Modern computing systems are capable of exascale calculations, which are revolutionizing the development and application of high-fidelity numerical models in computational science and engineering. While these systems continue to grow in processing power, the available system memory has not increased commensurately, and electrical power consumption continues to grow. A predominant approach to limit the memory usage in large-scale applications is to exploit the abundant processing power and continually recompute many low-level simulation quantities, rather than storing them. However, this approach can adversely impact the throughput of the simulation and diminish the benefits of modern computing architectures. We present three novel contributions to reduce the memory burden while maintaining, and sometimes improving, performance in simulations based on finite element discretizations. The first contribution develops dictionary-based data compression schemes that detect and exploit the structure of the discretization, due to redundancies across the finite element mesh. While these schemes are shown to reduce memory requirements by more than 99% on meshes with large numbers of identical mesh cells, there are applications where this structure does not exist. The second contribution leverages a recently developed augmented Lagrangian optimization algorithm to enable r-adaptivity for meshes with the goal of enhancing the redundancies in the mesh. The third contribution extends these methods to patch-based linear solvers and preconditioners by compressing local matrices. Numerical results demonstrate the effectiveness of the proposed methods to detect, enhance and exploit mesh structure on a suite of examples inspired by large-scale applications.

97 MATHEMATICS AND COMPUTING↗

A variational method for finite element stress recovery and error estimation

A variational method for obtaining smoothed stresses from a finite element derived nonsmooth stress field is presented. The method is based on minimizing a functional involving discrete least-squares error plus a penalty constraint that ensures smoothness of the stress field. An equivalent accuracy criterion is developed for the smoothing analysis which results in a C sup 1-continuous smoothed stress field possessing the same order of accuracy as that found at the superconvergent optimal stress points of the original finite element analysis. Application of the smoothing analysis to residual error estimation is also demonstrated.

Tessler, A.↗

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗

Approaches to solar cell design for pulsed laser power receivers

Using a laser to beam power from Earth to a photovoltaic receiver in space could be a technology with applications to many space missions. Extremely high average-power lasers would be required in a wavelength range of 700-1000 nm. However, high-power lasers inherently operate in a pulsed format. Existing solar cells are not well designed to respond to pulsed incident power. To better understand cell response to pulsed illumination at high intensity, the PC-1D finite-element computer model was used to analyze the response of solar cells to continuous and pulsed laser illumination. Over 50 percent efficiency was calculated for both InP and GaAs cells under steady-state illumination near the optimum wavelength. The time-dependent response of a high-efficiency GaAs concentrator cell to a laser pulse was modeled, and the effect of laser intensity, wavelength, and bias point was studied. Three main effects decrease the efficiency of a solar cell under pulsed laser illumination: series resistance, L-C 'ringing' with the output circuit, and current limiting due to the output inductance. The problems can be solved either by changing the pulse shape or designing a solar cell to accept the pulsed input. Cell design possibilities discussed are a high-efficiency, light-trapping silicon cell, and a monolithic, low-inductance GaAs cell.

Jain, Raj K.↗

Thermal stress analysis of composites in the space environment

A finite element micromechanics approach was utilized to investigate the thermally induced stress fields in continuous fiber reinforced polymer matrix composites at temperatures typical of spacecraft operating environments. The influence of laminate orientation was investigated with a simple global/local formulation. Thermal stress calculations were used to predict probable damage initiation locations, and the results were compared to experimentally observed damage in several epoxy matrix composites. The influence of an interphase region on the interfacial stress states was investigated.

Bowles, David E.↗

An improved plate theory of order (1,2) for thick composite laminates

A new (1,2)-order theory is proposed for the linear elasto-static analysis of laminated composite plates. The basic assumptions are those concerning the distribution through the laminate thickness of the displacements, transverse shear strains and the transverse normal stress, with these quantities regarded as some weighted averages of their exact elasticity theory representations. The displacement expansions are linear for the inplane components and quadratic for the transverse component, whereas the transverse shear strains and transverse normal stress are respectively quadratic and cubic through the thickness. The main distinguishing feature of the theory is that all strain and stress components are expressed in terms of the assumed displacements prior to the application of a variational principle. This is accomplished by an a priori least-square compatibility requirement for the transverse strains and by requiring exact stress boundary conditions at the top and bottom plate surfaces. Equations of equilibrium and associated Poisson boundary conditions are derived from the virtual work principle. It is shown that the theory is particularly suited for finite element discretization as it requires simple C(sup 0)- and C(sup -1)-continuous displacement interpolation fields. Analytic solutions for the problem of cylindrical bending are derived and compared with the exact elasticity solutions and those of our earlier (1,2)-order theory based on the assumed displacements and transverse strains.

Tessler, A.↗

Trabecular Bone Mechanical Properties and Fractal Dimension

Countermeasures for reducing bone loss and muscle atrophy due to extended exposure to the microgravity environment of space are continuing to be developed and improved. An important component of this effort is finite element modeling of the lower extremity and spinal column. These models will permit analysis and evaluation specific to each individual and thereby provide more efficient and effective exercise protocols. Inflight countermeasures and post-flight rehabilitation can then be customized and targeted on a case-by-case basis. Recent Summer Faculty Fellowship participants have focused upon finite element mesh generation, muscle force estimation, and fractal calculations of trabecular bone microstructure. Methods have been developed for generating the three-dimensional geometry of the femur from serial section magnetic resonance images (MRI). The use of MRI as an imaging modality avoids excessive exposure to radiation associated with X-ray based methods. These images can also detect trabecular bone microstructure and architecture. The goal of the current research is to determine the degree to which the fractal dimension of trabecular architecture can be used to predict the mechanical properties of trabecular bone tissue. The elastic modulus and the ultimate strength (or strain) can then be estimated from non-invasive, non-radiating imaging and incorporated into the finite element models to more accurately represent the bone tissue of each individual of interest. Trabecular bone specimens from the proximal tibia are being studied in this first phase of the work. Detailed protocols and procedures have been developed for carrying test specimens through all of the steps of a multi-faceted test program. The test program begins with MRI and X-ray imaging of the whole bones before excising a smaller workpiece from the proximal tibia region. High resolution MRI scans are then made and the piece further cut into slabs (roughly 1 cm thick). The slabs are X-rayed again and also scanned using dual-energy X-ray absorptiometry (DEXA). Cube specimens are then cut from the slabs and tested mechanically in compression. Correlations between mechanical properties and fractal dimension will then be examined to assess and quantify the predictive capability of the fractal calculations.

Hogan, Harry A.↗

System-Integrated Finite Element Analysis of a Full-Scale Helicopter Crash Test with Deployable Energy Absorbers

A full-scale crash test of an MD-500 helicopter was conducted in December 2009 at NASA Langley's Landing and Impact Research facility (LandIR). The MD-500 helicopter was fitted with a composite honeycomb Deployable Energy Absorber (DEA) and tested under vertical and horizontal impact velocities of 26-ft/sec and 40-ft/sec, respectively. The objectives of the test were to evaluate the performance of the DEA concept under realistic crash conditions and to generate test data for validation of a system integrated finite element model. In preparation for the full-scale crash test, a series of sub-scale and MD-500 mass simulator tests was conducted to evaluate the impact performances of various components, including a new crush tube and the DEA blocks. Parameters defined within the system integrated finite element model were determined from these tests. The objective of this paper is to summarize the finite element models developed and analyses performed, beginning with pre-test predictions and continuing through post-test validation.

Annett, Martin S.↗

TPSAS-NF1676L-13539-DND

The performance of variable stiffness composite panels is characterized using experiments and structural analyses. Fabricated using a fiber placement system, the fiber angle in these panels varies continuously within each ply. Measured panel shapes and thermal responses are correlated with finite element analyses. Large geometric imperfections require unique fixtures to both straighten the panel edges and provide structural test boundary conditions. Panels show linear prebuckling and nonlinear postbuckling under applied end shortening, which occur at higher loads than a baseline panel. Panel responses to in-plane loading are computed using geometric and materially nonlinear analyses that include mechanical prestresses and measured imperfections. Analytical strains and displacements also correlate well with measured results.

Chauncey Wu↗

On Practical Aspects of Variational Consistency in Contact Dynamics

Usage of contact mechanics methodologies is a pervasive modeling requirement in dynamic simulations. While for some trivial problems, solutions taken from analytical geometry are available, use of a finite element framework is common to achieve formulation generality. This work explores two dynamic contact formulations: one based on the traditional node-to-segment (NTS) approach, and a variationally consistent segment-to-segment (STS) mortar formulation. The NTS formulation employed here enforces the constraints kinematically (i.e., the interpenetration is enforced to the solver tolerance), whereas the mortar approach uses Lagrange multipliers to enforce the contact constraints. Both approaches are implemented in the open-source finite element framework Multiphysics Object-Oriented Simulation Environment (MOOSE). The results highlight two relevant contact-interface-related dynamic phenomena in finite element simulations. First, stabilization of contact constraints is discussed, taking into account the evolution of the total energy in a benchmark problem. Second, the influence of finite element discretization on both of the aforementioned contact formulations is analyzed by exercising a large-deformation example with continuous relative sliding. Variationally consistent contact approaches such as the mortar formulation lead to improved energy preservation and avoid spurious excitation of the system's frequencies. This is especially relevant in settings where inertia and vibrations are of importance.

42 ENGINEERING↗

Scale-bridging with the extended/generalized finite element method for linear elastodynamics

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. Furthermore, this introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

36 MATERIALS SCIENCE↗

Vibration of thick composite laminates - Analytic theory and finite element approximations

Two versions of a (1,2) higher-order laminate plate theory are extended to elastodynamics, and their performance on free vibration problems for homogeneous and laminated plates is examined. Both versions are based on continuously varying through-the-thickness linear and parabolic distributions for the inplane and transverse displacements, respectively. They also use continuous parabolic distributions for the transverse shear strains. The theory is capable of providing an improved accuracy as compared to other approximate theories and accurately predicting the lowest thickness-stretch frequencies which cannot be obtained with shear-deformable theories. The theory is considered to be particularly well-suited for generating simple and efficient variationally-based displacement finite elements leading to C sup 0 elements with only five degrees of freedom per node.

Tessler, A.↗