SEARCH · Search NASA

Results for “Finite element method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1985

A high-precision finite element method for shock-tube calculations

A two-pass explicit scheme is developed in order to exploit some of the capabilities of finite difference modeling (FDM) for finite element modeling (FEM), which offers the opportunity to account for any type of geometry in fluid flow modeling. Features of the first-order upwind and the Lax-Wendroff high precision explicit finite difference algorithms are reviewed. A flux limiter is developed for FEM to serve as an analog for the single limiter function which has been defined for the various FDMs. It is shown that an antidiffusive limiter must be introduced into the weighting function which normally multiplies the time-derivative term in the variational equation. The two-pass scheme which results is demonstrated to be the equivalent of FDMs with five-point support. However, the present scheme is valid only for one-dimensional calculations and linear shape functions for shock tube flow phenomena. Further work is required for its use with nonlinear hyperbolic systems.

Hughes, T. J. R.↗

NASA NTRS1994

A coupled/uncoupled deformation and fatigue damage algorithm utilizing the finite element method

A fatigue damage computational algorithm utilizing a multiaxial, isothermal, continuum based fatigue damage model for unidirectional metal matrix composites has been implemented into the commercial finite element code MARC using MARC user subroutines. Damage is introduced into the finite element solution through the concept of effective stress which fully couples the fatigue damage calculations with the finite element deformation solution. An axisymmetric stress analysis was performed on a circumferentially reinforced ring, wherein both the matrix cladding and the composite core were assumed to behave elastic-perfectly plastic. The composite core behavior was represented using Hill's anisotropic continuum based plasticity model, and similarly, the matrix cladding was represented by an isotropic plasticity model. Results are presented in the form of S-N curves and damage distribution plots.

Wilt, Thomas E.↗

NASA NTRS1992

Material nonlinear analysis via mixed-iterative finite element method

The performance of elastic-plastic mixed-iterative analysis is examined through a set of convergence studies. Membrane and bending behaviors are tested using 4-node quadrilateral finite elements. The membrane result is excellent, which indicates the implementation of elastic-plastic mixed-iterative analysis is appropriate. On the other hand, further research to improve bending performance of the method seems to be warranted.

Sutjahjo, Edhi↗

NASA NTRS1990

An error analysis of least-squares finite element method of velocity-pressure-vorticity formulation for Stokes problem

A theoretical proof of the optimal rate of convergence for the least-squares method is developed for the Stokes problem based on the velocity-pressure-vorticity formula. The 2D Stokes problem is analyzed to define the product space and its inner product, and the a priori estimates are derived to give the finite-element approximation. The least-squares method is found to converge at the optimal rate for equal-order interpolation.

Chang, Ching L.↗

NASA NTRS1991

Finite element methods for integrated aerodynamic heating analysis

This report gives a description of the work which has been undertaken during the second year of a three year research program. The objectives of the program are to produce finite element based procedures for the solution of the large scale practical problems which are of interest to the Aerothermal Loads Branch (ALB) at NASA Langley Research Establishment. The problems of interest range from Euler simulations of full three dimensional vehicle configurations to local analyses of three dimensional viscous laminar flow. Adaptive meshes produced for both steady state and transient problems are to be considered. An important feature of the work is the provision of specialized techniques which can be used at ALB for the development of an integrated fluid/thermal/structural modeling capability.

Morgan, K.↗

NASA NTRS1977

Hamilton's principle - Finite-element methods and flexible body dynamics

A variational formulation is given for the equations of motion for an unconstrained elastic body, and Hamilton's principle is used to derive the equations of motion and deformation of the body. Finite-element approximations are developed for these dynamical equations with respect to a body axis system satisfying the mean axis conditions. The free-body influence matrix for the body then is developed in terms of the finite-element model parameters.

Cavin, R. K., III↗

NASA NTRS1991

A locally refined rectangular grid finite element method - Application to computational fluid dynamics and computational physics

The present FEM technique addresses both linear and nonlinear boundary value problems encountered in computational physics by handling general three-dimensional regions, boundary conditions, and material properties. The box finite elements used are defined by a Cartesian grid independent of the boundary definition, and local refinements proceed by dividing a given box element into eight subelements. Discretization employs trilinear approximations on the box elements; special element stiffness matrices are included for boxes cut by any boundary surface. Illustrative results are presented for representative aerodynamics problems involving up to 400,000 elements.

Young, David P.↗

NASA NTRS1989

A finite element method for nonlinear panel flutter

Panel flutter is the self excited oscillations of a plate in supersonic flow. Linear theory gives no information about the panel's deflections and stresses. Hence, the service life of the panel cannot be predicted by linear methods. Nonlinear structural theory determines the limit cycle oscillating frequency, and also panel deflections and stresses. Panel fatigue life, therefore, can be predicted. For a more thorough understanding of panel flutter behavior, the geometric nonlinearity effects due to large deflections must be considered in the formulation. The Galerkin's method was used in the investigation of limit cycle oscillations of panels in supersonic flow in the spatial domain, and the panel deflections. All nonlinear flutter methods and results presently available were limited to a linear or linearized aerodynamic theory. Extension of the finite element flutter formulation to include nonlinear hypersonic aerodynamic loading and finite element nonlinear panel flutter results are presented.

Mei, Chuh↗

NASA NTRS1987

Isoparametric finite element methods to estimate discontinuous coefficients in two-dimensional elliptic equations

A summary of findings is presented on the problem of estimating discontinuous coefficients, including the locations of interfaces at which discontinuities occur, for elliptic equations on a two-dimensional domain. Finite-element approximations are constructed that are based on an underlying grid that adaptively moves to follow locations of interfaces as they are iteratively updated in the numerical estimation process.

Lamm, Patricia K.↗

NASA NTRS1979

Multi-level adaptive finite element methods. 1: Variation problems

A general numerical strategy for solving partial differential equations and other functional problems by cycling between coarser and finer levels of discretization is described. Optimal discretization schemes are provided together with very fast general solvers. It is described in terms of finite element discretizations of general nonlinear minimization problems. The basic processes (relaxation sweeps, fine-grid-to-coarse-grid transfers of residuals, coarse-to-fine interpolations of corrections) are directly and naturally determined by the objective functional and the sequence of approximation spaces. The natural processes, however, are not always optimal. Concrete examples are given and some new techniques are reviewed. Including the local truncation extrapolation and a multilevel procedure for inexpensively solving chains of many boundary value problems, such as those arising in the solution of time-dependent problems.

Brandt, A.↗

NASA NTRS1986

The use of finite element methods for determining the cross-sectional temperature distribution in heat pipes

A model that is currently used to predict the priming and performance limitations of a monogroove heat pipe is expanded to include the boiling limitation and the cross-sectional temperature distribution as determined from a multidimensional finite element analysis technique. The improved model is verified experimentally and shown to accurately predict the cross-sectional temperature distribution when the heat flux distribution is known. The model provides a way to estimate the level at which nucleate boiling and the associated dryout of the capillary wick occurs.

Peterson, G. P.↗

NASA NTRS1989

A finite-element method for large-amplitude, two-dimensional panel flutter at hypersonic speeds

The nonlinear flutter behavior of a two-dimensional panel in hypersonic flow is investigated analytically. An FEM formulation based unsteady third-order piston theory (Ashley and Zartarian, 1956; McIntosh, 1970) and taking nonlinear structural and aerodynamic phenomena into account is derived; the solution procedure is outlined; and typical results are presented in extensive tables and graphs. A 12-element finite-element solution obtained using an alternative method for linearizing the assumed limit-cycle time function is shown to give predictions in good agreement with classical analytical results for large-amplitude vibration in a vacuum and large-amplitude panel flutter, using linear aerodynamics.

Mei, Chuh↗

NASA NTRS1984

Finite element method - A companion in experimental mechanics

The hybrid experimental-numerical procedure for structural analysis is described by its applications in fracture mechanics. The procedure was first verified by the excellent agreements between the dynamic stress intensity factors obtained directly by dynamic photoelasticity and those generated by the hybrid procedure where a dynamic finite element code was executed in its generation mode. The hybrid procedure was then used to determine the dynamic fracture toughness of reaction bonded silicon nitride.

Kobayashi, A. S.↗

NASA NTRS1986

An adaptive characteristic Petrov-Galerkin finite element method for convection-dominated linear and nonlinear parabolic problems in one space variable

The present adaptive FEM technique for convection-dominated problems is based on a Petrov-Galerkin scheme for spatial approximation, whose typical time-step employs test functions chosen to yield an approximate solution coinciding with the exact solutions at the finite element grid nodes. The derivation of truly local a posteriori error estimates is made possible by this procedure, which is also shown to be a very effective solver by the numerical examples presented.

Demkowicz, L.↗

NASA NTRS1992

Structural optimization of thin shells using finite element method

The objective of the present work was the structural optimization of thin shell structures that are subjected to stress and displacement constraints. In order to accomplish this, the structural optimization computer program DESAP1 was modified and improved. In the static analysis part of the DESAP1 computer program the torsional spring elements, which are used to analyze thin, shallow shell structures, were eliminated by modifying the membrane stiffness matrix of the triangular elements in the local coordinate system and adding a fictitious rotational stiffness matrix. This simplified the DESAP1 program input, improved the accuracy of the analysis, and saved computation time. In the optimization part of the DESAP1 program the stress ratio formula, which redesigns the thickness of each finite element of the structure, was solved by an analytical method. This scheme replaced the iterative solution that was previously used in the DESAP1 program, thus increasing the accuracy and speed of the design. The modified program was used to design a thin, cylindrical shell structure with optimum weight, and the results are reported in this paper.

Gotsis, Pascal K.↗

NASA NTRSDate not supplied

Modeling Uprighting of the Orion Crew Module using a Finite Element Method with Wave Coupling and Structural Deformability

The Crew Module Uprighting System (CMUS) for the Orion Crew Module (CM) includes five bags that inflate after splashdown to rotate and upright the CM as necessary. Wave conditions can be severe enough to damage the bags and tethers, which could make the recovery operation hazardous. The CMUS therefore has design limits that ultimately constrain where the CM can land based upon sea conditions. This manuscript describes a full fluid-structure interaction computational model being used for the CMUS. It is a finite element model that uses explicit time integration to directly simulate uprighting, a computational approach usually reserved for short duration, highly transient phenomena and novelly applied here. The goal of this current effort is to determine the feasibility of this approach for CMUS uprighting, and to estimate the loads the CMUS bags experience in different wave conditions. The results from this work will be used to supplement other experimental and computational data in support of NASA’s Artemis program.

John Puryear↗

NASA NTRS1996

Applications of Parallel Computation in Micro-Mechanics and Finite Element Method

This project discusses the application of parallel computations related with respect to material analyses. Briefly speaking, we analyze some kind of material by elements computations. We call an element a cell here. A cell is divided into a number of subelements called subcells and all subcells in a cell have the identical structure. The detailed structure will be given later in this paper. It is obvious that the problem is "well-structured". SIMD machine would be a better choice. In this paper we try to look into the potentials of SIMD machine in dealing with finite element computation by developing appropriate algorithms on MasPar, a SIMD parallel machine. In section 2, the architecture of MasPar will be discussed. A brief review of the parallel programming language MPL also is given in that section. In section 3, some general parallel algorithms which might be useful to the project will be proposed. And, combining with the algorithms, some features of MPL will be discussed in more detail. In section 4, the computational structure of cell/subcell model will be given. The idea of designing the parallel algorithm for the model will be demonstrated. Finally in section 5, a summary will be given.

Tan, Hui-Qian↗

NASA NTRS1998

Deterministic and Probabilistic Creep and Creep Rupture Enhancement to CARES/Creep: Multiaxial Creep Life Prediction of Ceramic Structures Using Continuum Damage Mechanics and the Finite Element Method

High temperature and long duration applications of monolithic ceramics can place their failure mode in the creep rupture regime. A previous model advanced by the authors described a methodology by which the creep rupture life of a loaded component can be predicted. That model was based on the life fraction damage accumulation rule in association with the modified Monkman-Grant creep ripture criterion However, that model did not take into account the deteriorating state of the material due to creep damage (e.g., cavitation) as time elapsed. In addition, the material creep parameters used in that life prediction methodology, were based on uniaxial creep curves displaying primary and secondary creep behavior, with no tertiary regime. The objective of this paper is to present a creep life prediction methodology based on a modified form of the Kachanov-Rabotnov continuum damage mechanics (CDM) theory. In this theory, the uniaxial creep rate is described in terms of stress, temperature, time, and the current state of material damage. This scalar damage state parameter is basically an abstract measure of the current state of material damage due to creep deformation. The damage rate is assumed to vary with stress, temperature, time, and the current state of damage itself. Multiaxial creep and creep rupture formulations of the CDM approach are presented in this paper. Parameter estimation methodologies based on nonlinear regression analysis are also described for both, isothermal constant stress states and anisothermal variable stress conditions This creep life prediction methodology was preliminarily added to the integrated design code CARES/Creep (Ceramics Analysis and Reliability Evaluation of Structures/Creep), which is a postprocessor program to commercially available finite element analysis (FEA) packages. Two examples, showing comparisons between experimental and predicted creep lives of ceramic specimens, are used to demonstrate the viability of this methodology and the CARES/Creep program.

Jadaan, Osama M.↗