Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Application of High Order Acoustic Finite Elements to Transmission Losses and Enclosure Problems

A family of acoustic finite elements was developed based on C continuity (acoustic pressure being the nodal variable) and the no-flow condition. The family include triangular, quadrilateral and hexahedral isoparametric elements with linear quadratic and cubic variation in modelling and distortion. Of greatest use in problems with irregular boundaries are the cubic isoparametric elements: the 32 node hexahedral element for three-dimensional systems; and the twelve node quadrilateral and ten node triangular elements for two-dimensional/axisymmetric applications. These elements were applied to problems involving cavity resonances, transmission loss in silencers and the study of end effects, using a Floating Point Systems 164 attached array processor accessed through an Amdahl 5860 mainframe. The elements are presently being used to study the end effects associated with duct terminations within finite enclosures. The transmission losses with various silencers and sidebranches in ducts is also being studied using the same elements.

Craggs, A.↗
DOE OSTI2019

Estimation of Errors in Stress Distributions Computed in Finite Element Simulations of Polycrystals

The accuracy of the stresses predicted from crystal plasticity-based finite element formulation depends on estimation and control of the errors associated with the discretization. In the current work, the errors in the stress distribution are estimated in virtual polycrystalline samples of α-phase titanium (hexagonal close-packed phase of Ti–6Al– 4 V). To estimate the error, the stress field, which does not possess inter-element continuity, is smoothed over a grain using an L 2 projection, thereby providing continuous stress distributions with inter-element continuity. The differences between the continuous (smooth) and discontinuous (raw) stress fields are calculated at individual Gauss quadrature points and used to estimate errors for corresponding elements and grains. Error estimations are performed for a Voronoi-tessellated microstructure, an equiaxed microstructure, and two microstructures with varying grain sizes for tensile loading extending into the fully plastic regime (≈ 5% extension). Magnitudes of the errors are found to depend on microstructural characteristics, particularly the shape and size of grains. Samples having variations in grain size or having less spherical grains exhibited higher errors than samples with uniformly sized, equiaxed grains, with the size variations having a more pronounced effect. In conclusion, errors correlate with proximity to grain boundaries at small (elastic) strains and with deformation-induced features (deformation bands) at large strains.

42 ENGINEERING↗
DOE OSTI2011

Finite element analysis of multilayer coextrusion.

Multilayer coextrusion has become a popular commercial process for producing complex polymeric products from soda bottles to reflective coatings. A numerical model of a multilayer coextrusion process is developed based on a finite element discretization and two different free-surface methods, an arbitrary-Lagrangian-Eulerian (ALE) moving mesh implementation and an Eulerian level set method, to understand the moving boundary problem associated with the polymer-polymer interface. The goal of this work is to have a numerical capability suitable for optimizing and troubleshooting the coextrusion process, circumventing flow instabilities such as ribbing and barring, and reducing variability in layer thickness. Though these instabilities can be both viscous and elastic in nature, for this work a generalized Newtonian description of the fluid is used. Models of varying degrees of complexity are investigated including stability analysis and direct three-dimensional finite element free surface approaches. The results of this work show how critical modeling can be to reduce build test cycles, improve material choices, and guide mold design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Efficient quadrature rules for finite element discretizations of nonlocal equations

Here, in this paper, we design efficient quadrature rules for finite element (FE) discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work, we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first show that the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient FE assembly.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 1: Model development and applications

This study concerns the rotor dynamic characteristics of fluid-encompassed rotors, with special emphasis on shrouded pump impellers. The core of the study is a versatile and categorically new finite-element-based perturbation model, which is based on a rigorous flow analysis and what we have generically termed the 'virtually' deformable finite-element approach. The model is first applied to the case of a smooth annular seal for verification purposes. The rotor excitation components, in this sample problem, give rise to a purely cylindrical, purely conical, and a simultaneous cylindrical/conical rotor whirl around the housing centerline. In all cases, the computed results are compared to existing experimental and analytical data involving the same seal geometry and operating conditions. Next, two labyrinth-seal configurations, which share the same tooth-to-tooth chamber geometry but differ in the total number of chambers, were investigated. The results, in this case, are compared to experimental measurements for both seal configurations. The focus is finally shifted to the shrouded-impeller problem, where the stability effects of the leakage flow in the shroud-to-housing secondary passage are investigated. To this end, the computational model is applied to a typical shrouded-impeller pump stage, fabricated and rotor dynamically tested by Sulzer Bros., and the results compared to those of a simplified 'bulk-flow' analysis and Sulzer Bros.' test data. In addition to assessing the computed rotor dynamic coefficients, the shrouded-impeller study also covers a controversial topic, namely that of the leakage-passage inlet swirl, which was previously cited as the origin of highly unconventional (resonance-like) trends of the fluid-exerted forces. In order to validate this claim, a 'microscopic' study of the fluid/shroud interaction mechanism is conducted, with the focus being on the structure of the perturbed flow field associated with the impeller whirl. The conclusions of this study were solidified by the outcome of a numerical-certainty exercise, where the grid dependency of the numerical results is objectively examined. The final phase of the shrouded-impeller investigation involves the validation of a built-in assumption, in all other perturbation models, whereby single-harmonic tangential distributions of all the flow thermophysical properties are imposed. The last phase of the investigation course is aimed at verifying the fine details of the perturbed flow field in light of recent set of detailed LDA measurements in a smooth annular seal. Grid dependency of the fluid-induced forces is also investigated, and specific recommendations are made.

Baskharone, Erian A.↗
NASA NTRS1993

Helicopter trim analysis by shooting and finite element methods with optimally damped Newton iterations

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used, and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Achar, N. S.↗
NASA NTRS2000

A Strategy for Integrating a Large Finite Element Model: X-33 Lessons Learned

The X-33 vehicle is an advanced technology demonstrator sponsored by NASA. For the past three years the Structural Dynamics & Loads Group of NASA's Marshall Space Flight Center has had the task of integrating the X-33 vehicle structural finite element model. In that time, five versions of the integrated vehicle model have been produced and a strategy has evolved that would benefit anyone given the task of integrating structural finite element models that have been generated by various modelers and companies. The strategy that has been presented here consists of six decisions that need to be made. These six decisions are: purpose of model, units, common material list, model numbering, interface control, and archive format. This strategy has been proved and expanded from experience on the X-33 vehicle.

McGhee, David S.↗
NASA NTRS1980

Application of a data base management system to a finite element model

In today's software market, much effort is being expended on the development of data base management systems (DBMS). Most commercially available DBMS were designed for business use. However, the need for such systems within the engineering and scientific communities is becoming apparent. A potential DBMS application that appears attractive is the handling of data for finite element engineering models. The applications of a commercially available, business-oriented DBMS to a structural engineering, finite element model is explored. The model, DBMS, an approach to using the DBMS, advantages and disadvantages are described. Plans for research on a scientific and engineering DBMS are discussed.

Rogers, J. L., Jr.↗
NASA NTRS1981

Acoustic transmission in non-uniform ducts with mean flow. I - The method of weighted residuals. II - The finite element method

The problem of acoustic transmission through nonuniform ducts containing a high-speed subsonic flow is studied by means of the method of weighted residuals in the form of a modified Galerkin method and a Galerkin formulation of the finite element method. The method of weighted residuals is shown to employ the basis functions generated from eigenvalue calculations for the case of no flow, and is verified by comparison with exact eigenvalue calculations in the uniform duct case and numerical solutions of the one-dimensional form of the equations in the nonuniform duct case. The finite element scheme based on both the Galerkin method and the residual least squares method and employing eight-noded isoparametric elements is presented and used to investigate multimodal propagation by the coupling of the solution in the duct nonuniform section to modal expansions in uniform sections. Comparison of the results of the two methods reveals them to be in substantial agreement, and predicts the importance of multimodal interactions at high Mach numbers.

Eversman, W.↗
NASA NTRS1975

A triangular thin shell finite element: Nonlinear analysis

Aspects of the formulation of a triangular thin shell finite element which pertain to geometrically nonlinear (small strain, finite displacement) behavior are described. The procedure for solution of the resulting nonlinear algebraic equations combines a one-step incremental (tangent stiffness) approach with one iteration in the Newton-Raphson mode. A method is presented which permits a rational estimation of step size in this procedure. Limit points are calculated by means of a superposition scheme coupled to the incremental side of the solution procedure while bifurcation points are calculated through a process of interpolation of the determinants of the tangent-stiffness matrix. Numerical results are obtained for a flat plate and two curved shell problems and are compared with alternative solutions.

Thomas, G. R.↗
NASA NTRS1983

Finite element computation with parallel VLSI

This paper describes a parallel processing computer consisting of a 16-bit microcomputer as a master processor which controls and coordinates the activities of 8086/8087 VLSI chip set slave processors working in parallel. The hardware is inexpensive and can be flexibly configured and programmed to perform various functions. This makes it a useful research tool for the development of, and experimentation with parallel mathematical algorithms. Application of the hardware to computational tasks involved in the finite element analysis method is demonstrated by the generation and assembly of beam finite element stiffness matrices. A number of possible schemes for the implementation of N-elements on N- or n-processors (N is greater than n) are described, and the speedup factors of their time consumption are determined as a function of the number of available parallel processors.

Mcgregor, J.↗
NASA NTRS2003

Full-Scale Crash Test and Finite Element Simulation of a Composite Prototype Helicopter

A full-scale crash test of a prototype composite helicopter was performed at the Impact Dynamics Research Facility at NASA Langley Research Center in 1999 to obtain data for validation of a finite element crash simulation. The helicopter was the flight test article built by Sikorsky Aircraft during the Advanced Composite Airframe Program (ACAP). The composite helicopter was designed to meet the stringent Military Standard (MIL-STD-1290A) crashworthiness criteria and was outfitted with two crew and two troop seats and four anthropomorphic dummies. The test was performed at 38-ft/s vertical and 32.5-ft/s horizontal velocity onto a rigid surface. An existing modal-vibration model of the Sikorsky ACAP helicopter was converted into a model suitable for crash simulation. A two-stage modeling approach was implemented and an external user-defined subroutine was developed to represent the complex landing gear response. The crash simulation was executed with a nonlinear, explicit transient dynamic finite element code. Predictions of structural deformation and failure, the sequence of events, and the dynamic response of the airframe structure were generated and the numerical results were correlated with the experimental data to validate the simulation. The test results, the model development, and the test-analysis correlation are described.

Jackson, Karen E.↗
NASA NTRS1974

A finite element analysis of fatigue crack closure

Experiments have shown that fatigue cracks close at positive loads during constant-amplitude load cycling. The crack-closure phenomenon is caused by residual plastic deformations remaining in the wake of an advancing crack tip. The present paper is concerned with the application of a two-dimensional, nonlinear, finite-element analysis for predicting crack-closure and crack-opening stresses during cyclic loading. A two-dimensional finite-element computer program, which accounts for both elastic-plastic material behavior and changing boundary conditions associated with crack extension and intermittent contact of the crack surfaces under cyclic loading, has been developed. An efficient technique to account for changing boundary conditions was also incorporated into the nonlinear analysis program. This program was subsequently used to study crack extension and crack closure under constant-amplitude and two-level block loading. The calculated crack-closure and crack-opening stresses were qualitatively consistent with experimental observations.

Newman, J. C., Jr.↗
NASA NTRS2018

Modeling Space Radiation Induced Bone Changes in Rat Femurs through Finite Element Analysis

As the duration of manned missions outside of the Earth’s protective shielding increase, astronauts are at risk for exposure to space radiation. Various organ systems may be damaged due to exposure. This study investigates the bonestrength changes using finite element modeling of Long Evans rats (n=85) subjected to graded, head-only proton (0, 10, 25, and 100 cGy, 150 MeV/n) and 28silicon (0, 10, 25, and 50 cGy, 300 MeV/n) radiation. The strength of the femoral neck will be examined due its clinical relevance to hip fractures. It has been shown in previous studies that bone mineral density was not reduced at the site of fracture. These findings question whether measurements of bone mineral density may be used to assess risk of hip fracture. The mechanisms leading to the irregular relationship between bone density and strength are still uncertain within literature and investigated to greater extent in clinical applications. Finite element analysis within this study simulated physiological loading of the femoral neck. No significant changes in femoral neck strength were found across doses of proton or 28silicon head-only radiation. Future work includes performing mechanical testing of the bone samples. Moving from mouse to larger animal models may also provide the increased lifespan for assessing the long-term outcomes of radiation exposure.

Dale Johnson↗
NASA NTRS1999

A Strategy for Integrating a Large Finite Element Model Using MSC NASTRAN/PATRAN: X-33 Lessons Learned

The X-33 vehicle is an advanced technology demonstrator sponsored by NASA. For the past 3 years the Structural Dynamics and Loads Branch of NASA's Marshall Space Flight Center has had the task of integrating the X-33 vehicle structural finite element model. In that time, five versions of the integrated vehicle model have been produced and a strategy has evolved that would benefit anyone given the task of integrating structural finite element models that have been generated by various modelers and companies. The strategy that has been presented here consists of six decisions that need to be made: purpose of models, units, common materials list, model numbering, interface control, and archive format. This strategy has been proven and expanded from experience on the X-33 vehicle.

McGhee, D. S.↗
DOE OSTI2019

Development of an Advanced Finite Element Model and Parametric Study to Evaluate Cement Sheath Barrier

Cement failure is known as one of the major causes for loss of well control events. Cement design is considered as one of the top technological knowledge gaps in high-pressure high-temperature oil and gas exploration. The primary objective of this paper is to perform a parametric analysis and identify critical parameters affecting the mechanical integrity of the set cement sheath. To achieve the objective, three-dimensional finite element models consisting of concentric casings and annular cement sheath were created. The finite element model was validated by analytical calculations. Performance of cement sheath was assessed by analyzing radial, hoop, and maximum shear stresses at different loading conditions. A parametric study was conducted by individually varying influencing factors such as cement material properties, sheath dimensions, and wellbore pressure loads. Values of all parameters were normalized and represented on the same plot against mechanical stresses. Such response curves can be used to estimate whether cement will structurally fail because of various operational loads or material aging. The plot can also be utilized to rank various factors in terms of influence on cement’s performance. Sensitivity response reveals that wellbore pressure, cement material properties, and annulus pressure are major parameters influencing mechanical stresses in neat class G cement. The order of importance depends on the type of stress. Results indicate interfacial bond failure and radial cracking to be the more likely modes of failure for class G cement. Cement response curves can help design engineers and regulators alike in quickly evaluating short-term or long-term fitness-for-service of cement sheath from the perspective of structural integrity. Industry standards and guidelines can be improved by adding performance curves for standard cement recipes.

Energy & Fuels↗
NASA NTRS2003

An Integrated Magnetic Circuit Model and Finite Element Model Approach to Magnetic Bearing Design

A code for designing magnetic bearings is described. The code generates curves from magnetic circuit equations relating important bearing performance parameters. Bearing parameters selected from the curves by a designer to meet the requirements of a particular application are input directly by the code into a three-dimensional finite element analysis preprocessor. This means that a three-dimensional computer model of the bearing being developed is immediately available for viewing. The finite element model solution can be used to show areas of magnetic saturation and make more accurate predictions of the bearing load capacity, current stiffness, position stiffness, and inductance than the magnetic circuit equations did at the start of the design process. In summary, the code combines one-dimensional and three-dimensional modeling methods for designing magnetic bearings.

Provenza, Andrew J.↗