Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite elements”

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 145 records · Page 8

Probabilistic finite elements for transient analysis in nonlinear continua

The probabilistic finite element method (PFEM), which is a combination of finite element methods and second-moment analysis, is formulated for linear and nonlinear continua with inhomogeneous random fields. Analogous to the discretization of the displacement field in finite element methods, the random field is also discretized. The formulation is simplified by transforming the correlated variables to a set of uncorrelated variables through an eigenvalue orthogonalization. Furthermore, it is shown that a reduced set of the uncorrelated variables is sufficient for the second-moment analysis. Based on the linear formulation of the PFEM, the method is then extended to transient analysis in nonlinear continua. The accuracy and efficiency of the method is demonstrated by application to a one-dimensional, elastic/plastic wave propagation problem. The moments calculated compare favorably with those obtained by Monte Carlo simulation. Also, the procedure is amenable to implementation in deterministic FEM based computer programs.

Liu, W. K.↗

Ablative Thermal Response Analysis Using the Finite Element Method

A review of the classic techniques used to solve ablative thermal response problems is presented. The advantages and disadvantages of both the finite element and finite difference methods are described. As a first step in developing a three dimensional finite element based ablative thermal response capability, a one dimensional computer tool has been developed. The finite element method is used to discretize the governing differential equations and Galerkin's method of weighted residuals is used to derive the element equations. A code to code comparison between the current 1-D tool and the 1-D Fully Implicit Ablation and Thermal Response Program (FIAT) has been performed.

Dec John A.↗

A least-squares finite element method for incompressible Navier-Stokes problem

Most finite element schemes for solving the Navier-Stokes equations can be categorized into the Galerkin mixed method and the penalty method. The mixed method leads to a saddle-point problem. In order to guarantee the existence of a solution, the combination of velocity and pressure interpolations requires satisfaction of the Ladyzhenskaya Babuska Brezzi (LBB) consistency condition which precludes the use of equal order interpolations and many seemingly natural pairs of velocity and pressure elements. In a previous paper a least-squares finite element method based on the first order velocity-pressure-vorticity formulation for the Stokes problem was proposed. This method leads to a minimization problem. The choice of combination of elements is thus not subject to the LBB condition. The numerical experiments exhibit the optimal rate of convergence for all variables with equal order interpolations. A theoretical error analysis supports the numerical results. In this paper the least-squares finite element method is extended to solving the incompressible Navier-Stokes problem.

Jiang, Bonan↗

Finite element solution algorithm for incompressible fluid dynamics

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the transient motion of a viscous incompressible fluid, i.e., hydrodynamics. Dependent variable transformation renders the differential equation description uniformly elliptic. The finite element algorithm is established using the Galerkin criterion on a local basis within the Method of Weighted Residuals. It is unconstrained with respect to system linearity, computational mesh uniformity or solution domain closure regularity. The finite element matrices are established using a linear 'natural coordinate function' description. Computational solutions using the COMOC computer program illustrate the various features of the algorithm including recirculating flows.

Baker, A. J.↗

Finite-Element Methods for Real-Time Simulation of Surgery

Two finite-element methods have been developed for mathematical modeling of the time-dependent behaviors of deformable objects and, more specifically, the mechanical responses of soft tissues and organs in contact with surgical tools. These methods may afford the computational efficiency needed to satisfy the requirement to obtain computational results in real time for simulating surgical procedures as described in Simulation System for Training in Laparoscopic Surgery (NPO-21192) on page 31 in this issue of NASA Tech Briefs. Simulation of the behavior of soft tissue in real time is a challenging problem because of the complexity of soft-tissue mechanics. The responses of soft tissues are characterized by nonlinearities and by spatial inhomogeneities and rate and time dependences of material properties. Finite-element methods seem promising for integrating these characteristics of tissues into computational models of organs, but they demand much central-processing-unit (CPU) time and memory, and the demand increases with the number of nodes and degrees of freedom in a given finite-element model. Hence, as finite-element models become more realistic, it becomes more difficult to compute solutions in real time. In both of the present methods, one uses approximate mathematical models trading some accuracy for computational efficiency and thereby increasing the feasibility of attaining real-time up36 NASA Tech Briefs, October 2003 date rates. The first of these methods is based on modal analysis. In this method, one reduces the number of differential equations by selecting only the most significant vibration modes of an object (typically, a suitable number of the lowest-frequency modes) for computing deformations of the object in response to applied forces.

Basdogan, Cagatay↗

Evaluation of an improved finite-element thermal stress calculation technique

A procedure for generating accurate thermal stresses with coarse finite element grids (Ojalvo's method) is described. The procedure is based on the observation that for linear thermoelastic problems, the thermal stresses may be envisioned as being composed of two contributions; the first due to the strains in the structure which depend on the integral of the temperature distribution over the finite element and the second due to the local variation of the temperature in the element. The first contribution can be accurately predicted with a coarse finite-element mesh. The resulting strain distribution can then be combined via the constitutive relations with detailed temperatures from a separate thermal analysis. The result is accurate thermal stresses from coarse finite element structural models even where the temperature distributions have sharp variations. The range of applicability of the method for various classes of thermostructural problems such as in-plane or bending type problems and the effect of the nature of the temperature distribution and edge constraints are addressed. Ojalvo's method is used in conjunction with the SPAR finite element program. Results are obtained for rods, membranes, a box beam and a stiffened panel.

Camarda, C. J.↗

Analysis of Large Quasistatic Deformations of Inelastic Solids by a New Stress Based Finite Element Method

A new hybrid stress finite element algorithm suitable for analyses of large quasistatic deformation of inelastic solids is presented. Principal variables in the formulation are the nominal stress rate and spin. The finite element equations which result are discrete versions of the equations of compatibility and angular momentum balance. Consistent reformulation of the constitutive equation and accurate and stable time integration of the stress are discussed at length. Examples which bring out the feasibility and performance of the algorithm conclude the work.

Reed, Kenneth W.↗

A new parallel-vector finite element analysis software on distributed-memory computers

A new parallel-vector finite element analysis software package MPFEA (Massively Parallel-vector Finite Element Analysis) is developed for large-scale structural analysis on massively parallel computers with distributed-memory. MPFEA is designed for parallel generation and assembly of the global finite element stiffness matrices as well as parallel solution of the simultaneous linear equations, since these are often the major time-consuming parts of a finite element analysis. Block-skyline storage scheme along with vector-unrolling techniques are used to enhance the vector performance. Communications among processors are carried out concurrently with arithmetic operations to reduce the total execution time. Numerical results on the Intel iPSC/860 computers (such as the Intel Gamma with 128 processors and the Intel Touchstone Delta with 512 processors) are presented, including an aircraft structure and some very large truss structures, to demonstrate the efficiency and accuracy of MPFEA.

Qin, Jiangning↗

Nonlinear crack analysis with finite elements

The application of finite element techniques to the analytic representation of the nonlinear behavior of arbitrary two-dimensional bodies containing cracks is discussed. Specific methods are proposed using which it should be possible to obtain information concerning: the description of the maximum, minimum, and residual near-tip stress and strain fields; the effects of crack closure on the near-tip behavior of stress and strain fields during cyclic loading into the plastic range; the stress-strain and displacement field behavior associated with a nonstationary crack; and the effects of large rotation near the crack tip.

Armen, H., Jr.↗

A numerically efficient finite element hydroelastic analysis

A finite element hydroelastic analysis formulation is developed on the basis of Toupin's complementary variational principle. Emphasis is placed on the special case of an incompressible fluid model which is applicable to propellant tank hydroelastic analysis. A concise fluid inertia representation results from the assumption of incompressibility and the hydroelastic equations reduce to a simplified form associated with the structure alone. The efficiency of the incompressible hydroelastic formulation is enhanced for both fluid and structure by introduction of harmonic reduction as an alternative to Guyan reduction. The theoretical developments are implemented in the NASTRAN Program and the technique is verified and demonstrated as an efficient and accurate approach with a series of illustrative problems including the 1/8-scale Space Shuttle external tank.

Coppolino, R.↗

Progressive failure of notched composite laminates using finite elements

An incremental, nonlinear, finite element program capable of tracing the damage that occurs in a notched laminated composite, before ultimate failure takes place, is described. Computer predictions for the ultimate failure load and failure modes are compared to experimental results for two laminates containing holes and loaded in uniaxial tension. The comparison is quite favorable.

Nuismer, R. J.↗

A numerically efficient finite element hydroelastic analysis

A finite element hydroelastic analysis formulation is developed on the basis of Toupin's complementary variational principle. Emphasis is placed on the special case of an incompressible fluid model which is applicable to propellant tank hydroelastic analysis. A concise fluid inertia representation results from the assumption of incompressibility and the hydroelastic equations reduce to a simplified form associated with the structure alone. The efficiency of the incompressible hydroelastic formulation in unhanced for both fluid and structure by introduction of harmonic reduction as an alternative to Guyan reduction. The theoretical developments are implemented in the NASTRAN Program and the technique is verified and demonstrated as an efficient and accurate approach with a series of illustrative problems including the 1/8 scale space shuttle external tank.

Coppolino, R. N.↗

Origin of multi-ring basin ridge systems - An upper limit to elastic deformation based on a finite-element model

A finite-element model of a 500-km diameter lunar basin flooded with 10 km of basalt was used to represent the central portion of a multi-ring basin undergoing elastic deformation as a result of basin subsidence. Gravitational forces due to the weight of the basalt were set as the initial forces in the model; elastic constants computed from seismic velocities of lunar samples provided strength specifications. In addition, the crust was constrained by permitting no displacements at a depth of 60 km. Because of the strength of the lithosphere, surface subsidence amounted to only 20 m. This magnitude of elastic deformation does not equal that found on the lunar surface, but the orientation of principal stress and the position of circumferential compressional stress are compatible with the evolution of inner-ring and radial ridge systems by gravitational-induced downdropping.

Maxwell, T. A.↗

Grid generation for two-dimensional finite element flowfield computation

The finite element method for fluid dynamics was used to develop a two dimensional mesh generation scheme. The method consists of shearing and conformal maps with upper and lower surfaces handled independently to allow sharp leading edges. The method also generates meshes of triangular or quadrilateral elements.

Tatum, K. E.↗

The solution of non-linear hyperbolic equation systems by the finite element method

A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.

Loehner, R.↗

An adaptive remeshing method for finite element thermal analysis

A finite element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗