Search NASA⌕ Search

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.

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

At least 235 records · Page 13

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.↗

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.↗

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↗

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.↗

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.↗

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.↗

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↗

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↗

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.↗

Water impact analysis of space shuttle solid rocket motor by the finite element method

Preliminary analysis showed that the doubly curved triangular shell elements were too stiff for these shell structures. The doubly curved quadrilateral shell elements were found to give much improved results. A total of six load cases were analyzed in this study. The load cases were either those resulting from a static test using reaction straps to simulate the drop conditions or under assumed hydrodynamic conditions resulting from a drop test. The latter hydrodynamic conditions were obtained through an emperical fit of available data. Results obtained from a linear analysis were found to be consistent with results obtained elsewhere with NASTRAN and BOSOR. The nonlinear analysis showed that the originally assumed loads would result in failure of the shell structures. The nonlinear analysis also showed that it was useful to apply internal pressure as a stabilizing influence on collapse. A final analysis with an updated estimate of load conditions resulted in linear behavior up to full load.

Buyukozturk, O.↗

A finite element method for potential aerodynamics around complex configurations

A general formulation for steady and oscillatory, subsonic and supersonic, potential linearized aerodynamic flow around complex configurations is presented. A linear integral equation relating the unknown potential on the surface of the body to the known downwash is used. The formulation is applied to the analysis of the flow field around wings and wing-body combinations. The surface is divided into small quadrilateral elements which are approximated with a hyperboloidal surface. The potential is assumed to be constant within each element. This yields a set of linear algebraic equations. The coefficients are evaluated analytically. Numerical results for steady and oscillatory, subsonic and supersonic flows indicate that the method, intrinsically general and flexible, is also fast, accurate and in excellent agreement with existing results.

Chen, L.-T.↗

A finite-element method for lifting surfaces in steady incompressible subsonic flow

The problem of potential steady subsonic flow for lifting surfaces is considered. This problem requires the solution of an integral equation relating the value of the potential discontinuity on the lifting surface and its wake to the values of the normal derivative of the potential which are known from the boundary conditions. The lifting surface is divided into small (quadrilateral hyperboloidal) surface elements, which are described in terms of the Cartesian components of the four corner points. The values of the potential discontinuity and the normal derivative of the potential are assumed to be constant within each element and equal to their values at the centroids of the elements. This yields a set of linear algebraic equations. Numerical results are in good agreement with existing ones.

Morino, L.↗

Some inconsistencies of the finite element method as applied to inelastic response

The inadequacy of a two noded beam-column element with a linear axial and a cubic transverse displacement field for inelastic analysis is demonstrated. For complete equilibrium satisfaction in the linear elastic range a three noded beam-column element is shown to be consistent. Next, the sensitivity of the inelastic response to numerical solutions of the inelastic response of a cantilever beam resulting from approximate integration of strain energy are brought out and finally, consequences of this on the nonlinear transient response of structures are considered.

Kamat, M. P.↗

Fillet Weld Stress Using Finite Element Methods

Average elastic Von Mises equivalent stresses were calculated along the throat of a single lap fillet weld. The average elastic stresses were compared to initial yield and to plastic instability conditions to modify conventional design formulas is presented. The factor is a linear function of the thicknesses of the parent plates attached by the fillet weld.

Lehnhoff, T. F.↗

Analysis of hourglass instabilities and control in underintegrated finite element methods

Belytschko et al. (1981, 1984) has developed stabilization methods for the treatment of underintegrated FEM problems; these methods involve the computation of an underintegrated stiffness matrix, which is rank-deficient, and the addition of a stabilization matrix which effectively eliminates the spurious modes. An attempt is presently made to give this a priori stabilization method a mathematical means of support. Attention is also given to an a posteriority stabilization method for hourglass control, in which an approximate solution of the underintegrated system is obtained and then subjected to a special projection in order to eliminate the hourglass modes. A proof is obtained for the convergence of this stabilized underintegrated approximation to the exact solution of a model problem at almost the same rate (as the mesh is refined) as the fully integrated solutions.

Jacquotte, O.-P.↗