SEARCH · Search NASA

Results for “finite element problems”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Transient modeling/analysis of hyperbolic heat conduction problems employing mixed implicit-explicit alpha method

This paper describes the evaluation of mixed implicit-explicit finite element formulations for hyperbolic heat conduction problems involving non-Fourier effects. In particular, mixed implicit-explicit formulations employing the alpha method proposed by Hughes et al. (1987, 1990) are described for the numerical simulation of hyperbolic heat conduction models, which involves time-dependent relaxation effects. Existing analytical approaches for modeling/analysis of such models involve complex mathematical formulations for obtaining closed-form solutions, while in certain numerical formulations the difficulties include severe oscillatory solution behavior (which often disguises the true response) in the vicinity of the thermal disturbances, which propagate with finite velocities. In view of these factors, the alpha method is evaluated to assess the control of the amount of numerical dissipation for predicting the transient propagating thermal disturbances. Numerical test models are presented, and pertinent conclusions are drawn for the mixed-time integration simulation of hyperbolic heat conduction models involving non-Fourier effects.

Tamma, Kumar K.↗

NASA NTRS1979

The application of the least squares finite element method to Abel's integral equation

Abel's integral equation is the governing equation for certain problems in physics and engineering, such as radiation from distributed sources. The finite element method for the solution of this non-linear equation is presented for problems with cylindrical symmetry and the extension to more general integral equations is indicated. The technique was applied to an axisymmetric glow discharge problem and the results show excellent agreement with previously obtained solutions

Balasubramanian, R.↗

NASA NTRS1985

Probabilistic Structural Analysis Theory Development

The objective of the Probabilistic Structural Analysis Methods (PSAM) project is to develop analysis techniques and computer programs for predicting the probabilistic response of critical structural components for current and future space propulsion systems. This technology will play a central role in establishing system performance and durability. The first year's technical activity is concentrating on probabilistic finite element formulation strategy and code development. Work is also in progress to survey critical materials and space shuttle mian engine components. The probabilistic finite element computer program NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) is being developed. The final probabilistic code will have, in the general case, the capability of performing nonlinear dynamic of stochastic structures. It is the goal of the approximate methods effort to increase problem solving efficiency relative to finite element methods by using energy methods to generate trial solutions which satisfy the structural boundary conditions. These approximate methods will be less computer intensive relative to the finite element approach.

Burnside, O. H.↗

NASA NTRS1988

Coupled bending-torsion steady-state response of pretwisted, nonuniform rotating beams using a transfer-matrix method

Using the Newtonian method, the equations of motion are developed for the coupled bending-torsion steady-state response of beams rotating at constant angular velocity in a fixed plane. The resulting equations are valid to first order strain-displacement relationships for a long beam with all other nonlinear terms retained. In addition, the equations are valid for beams with the mass centroidal axis offset (eccentric) from the elastic axis, nonuniform mass and section properties, and variable twist. The solution of these coupled, nonlinear, nonhomogeneous, differential equations is obtained by modifying a Hunter linear second-order transfer-matrix solution procedure to solve the nonlinear differential equations and programming the solution for a desk-top personal computer. The modified transfer-matrix method was verified by comparing the solution for a rotating beam with a geometric, nonlinear, finite-element computer code solution; and for a simple rotating beam problem, the modified method demonstrated a significant advantage over the finite-element solution in accuracy, ease of solution, and actual computer processing time required to effect a solution.

Gray, Carl E., Jr.↗

NASA NTRS1993

Coupled structural, thermal, phase-change and electromagnetic analysis for superconductors, volume 1

This research program has dealt with the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements; (2) finite element modeling of the electromagnetic problem; (3) coupling of thermal and mechanical effects; and (4) computer implementation and solution of the superconductivity transition problem. The research was carried out over the period September 1988 through March 1993. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles; (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements; and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects; and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The grant has fully supported the thesis work of one doctoral student (James Schuler, who started on January 1989 and completed on January 1993), and partly supported another thesis (Carmelo Militello, who started graduate work on January 1988 completing on August 1991). Twenty-three publications have acknowledged full or part support from this grant, with 16 having appeared in archival journals and 3 in edited books or proceedings.

Felippa, Carlos A.↗

NASA NTRS1995

Adaptive Finite Element Methods for Continuum Damage Modeling

The paper presents an application of adaptive finite element methods to the modeling of low-cycle continuum damage and life prediction of high-temperature components. The major objective is to provide automated and accurate modeling of damaged zones through adaptive mesh refinement and adaptive time-stepping methods. The damage modeling methodology is implemented in an usual way by embedding damage evolution in the transient nonlinear solution of elasto-viscoplastic deformation problems. This nonlinear boundary-value problem is discretized by adaptive finite element methods. The automated h-adaptive mesh refinements are driven by error indicators, based on selected principal variables in the problem (stresses, non-elastic strains, damage, etc.). In the time domain, adaptive time-stepping is used, combined with a predictor-corrector time marching algorithm. The time selection is controlled by required time accuracy. In order to take into account strong temperature dependency of material parameters, the nonlinear structural solution a coupled with thermal analyses (one-way coupling). Several test examples illustrate the importance and benefits of adaptive mesh refinements in accurate prediction of damage levels and failure time.

Min, J. B.↗

NASA NTRS2019

Benchmarking Mixed Mode Failure in Progressive Damage and Failure Analysis Methods

The verification and validation of progressive-damage-analysis finite element methods are difficult but critical tasks to undertake during their development. Verification exercises assess whether a predictive analysis tool produces results that are consistent with the fundamental concepts and assumptions of the tool under evaluation. Ideally, closed-form analytical solutions can be derived for which method verification results can be compared. Problems selected for computational tool verification are often simple and isolate individual features of the tool. In the case of progressive damage finite element methods, verifications should be performed to evaluate the ability of the model to predict the initiation of damage and its growth through the finite element mesh under a variety of conditions. Mabson et al. proposed a test case of a unidirectional, fiber-reinforced plate with a center crack subjected to tensile loads to evaluate matrix crack propagation predictions. The problem was modeled using the Abaqus Hashin continuum damage mechanics (CDM) model for fiber-reinforced composites. Different combinations of matrix strengths and element sizes were used in the simulations, and the results were compared to a closed-form solution based on linear elastic fracture mechanics (LEFM). It was determined that the Abaqus CDM model could predict the LEFM solution of Mode I cracks only when the finite element mesh density met specific requirements based on the material properties. This paper presents closed-form LEFM solutions for a center notch mixed mode (CNMM) verification problem. Parametric finite element analyses were developed using progressive damage analysis methods of both the Discrete Damage Mechanics (DDM) and CDM classes. The progressive damage analysis methods applied in the analyses of the CNMM problem include CompDam and the Floating Node Method. Analyses were conducted with various mode mixities and element sizes to verify that the damage models were working as intended and to identify any limits of applicability.

Leone, Frank↗

NASA NTRS2007

Compatibility Condition in Theory of Solid Mechanics (Elasticity, Structures, and Design Optimization)

The strain formulation in elasticity and the compatibility condition in structural mechanics have neither been understood nor have they been utilized. This shortcoming prevented the formulation of a direct method to calculate stress. We have researched and understood the compatibility condition for linear problems in elasticity and in finite element analysis. This has lead to the completion of the method of force with stress (or stress resultant) as the primary unknown. The method in elasticity is referred to as the completed Beltrami-Michell formulation (CBMF), and it is the integrated force method (IFM) in structures. The dual integrated force method (IFMD) with displacement as the primary unknown has been formulated. IFM and IFMD produce identical responses. The variational derivation of the CBMF yielded the new boundary compatibility conditions. The CBMF can be used to solve stress, displacement, and mixed boundary value problems. The IFM in structures produced high-fidelity response even with a modest finite element model. The IFM has influenced structural design considerably. A fully utilized design method for strength and stiffness limitation has been developed. The singularity condition in optimization has been identified. The CBMF and IFM tensorial approaches are robust formulations because of simultaneous emphasis on the equilibrium equation and the compatibility condition.

Patnaik, Surya N.↗

NASA NTRS1997

Modeling of Structural-Acoustic Interaction Using Coupled FE/BE Method and Control of Interior Acoustic Pressure Using Piezoelectric Actuators

A coupled finite element (FE) and boundary element (BE) approach is presented to model full coupled structural/acoustic/piezoelectric systems. The dual reciprocity boundary element method is used so that the natural frequencies and mode shapes of the coupled system can be obtained, and to extend this approach to time dependent problems. The boundary element method is applied to interior acoustic domains, and the results are very accurate when compared with limited exact solutions. Structural-acoustic problems are then analyzed with the coupled finite element/boundary element method, where the finite element method models the structural domain and the boundary element method models the acoustic domain. Results for a system consisting of an isotropic panel and a cubic cavity are in good agreement with exact solutions and experiment data. The response of a composite panel backed cavity is then obtained. The results show that the mass and stiffness of piezoelectric layers have to be considered. The coupled finite element and boundary element equations are transformed into modal coordinates, which is more convenient for transient excitation. Several transient problems are solved based on this formulation. Two control designs, a linear quadratic regulator (LQR) and a feedforward controller, are applied to reduce the acoustic pressure inside the cavity based on the equations in modal coordinates. The results indicate that both controllers can reduce the interior acoustic pressure and the plate deflection.

Mei, Chuh↗

NASA NTRS1990

Error detection and control for nonlinear shell analysis

A problem-adaptive solution procedure for improving the reliability of finite element solutions to geometrically nonlinear shell-type problem is presented. The strategy incorporates automatic error detection and control and includes an iterative procedure which utilizes the solution at the same load step on a more refined model. Representative nonlinear shell problem are solved.

Mccleary, Susan L.↗

NASA NTRS1996

Free-Surface and Contact Line Motion of Liquid in Microgravity

This project involves fundamental studies of the role of nonlinearity in determining the motion of liquid masses under the principal influences of surface tension, viscosity and inertia. Issues to be explored are relevant to aspects of terrestrial processes, as well as being immediately applicable to fluid management in a low-gravity environment. Specific issues include: (1) the mechanic's of liquid masses in large-amplitude motions, (2) the influence of bounding surfaces on the motion, and (3) the ability of such surfaces to control liquid motion by wetting forces, especially when they are augmented by various surface treatments. Mathematical techniques include asymptotic analysis of the governing equations, for problem simplification, and numerical simulation, using both boundary-element and finite-difference methods. The flow problem is divided into an 'outer' or inviscid potential-flow region and one or more inner, or viscous dominated, regions. Relevant to one inner region, the vicinity of the contact line, we discuss time-dependent simulation of slow droplet motion, on a surface of variable wettability, using the lubrication approximation. The simulation uses a disjoining pressure model and reproduces realistic wetting-dewetting behavior.

Schwartz, Leonard W.↗

NASA NTRS1993

NASTRAN analysis for the Airmass Sunburst model 'C' Ultralight Aircraft

The purpose of this project was to create a three dimensional NASTRAN model of the Airmass Sunburst Ultralight comparable to one made for finite element analysis. A two dimensional sample problem will be calculated by hand and by NASTRAN to make sure that NASTRAN finds similar results. A three dimensional model, similar to the one analyzed by the finite element program, will be run on NASTRAN. A comparison will be done between the NASTRAN results and the finite element program results. This study will deal mainly with the aerodynamic loads on the wing and surrounding support structure at an attack angle of 10 degrees.

Verbestel, John↗

NASA NTRS1983

New discretization and solution techniques for incompressible viscous flow problems

Several topics arising in the finite element solution of the incompressible Navier-Stokes equations are considered. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. The role of artificial viscosity in viscous flow calculations is studied, emphasizing work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some modifications are mentioned.

Gunzburger, M. D.↗

NASA NTRS1986

Preconditioned conjugate residual methods for the solution of spectral equations

Conjugate residual methods for the solution of spectral equations are described. An inexact finite-difference operator is introduced as a preconditioner in the iterative procedures. Application of these techniques is limited to problems for which the symmetric part of the coefficient matrix is positive definite. Although the spectral equation is a very ill-conditioned and full matrix problem, the computational effort of the present iterative methods for solving such a system is comparable to that for the sparse matrix equations obtained from the application of either finite-difference or finite-element methods to the same problems. Numerical experiments are shown for a self-adjoint elliptic partial differential equation with Dirichlet boundary conditions, and comparison with other solution procedures for spectral equations is presented.

Wong, Y. S.↗

NASA NTRS1991

STARS: An integrated general-purpose finite element structural, aeroelastic, and aeroservoelastic analysis computer program

The details of an integrated general-purpose finite element structural analysis computer program which is also capable of solving complex multidisciplinary problems is presented. Thus, the SOLIDS module of the program possesses an extensive finite element library suitable for modeling most practical problems and is capable of solving statics, vibration, buckling, and dynamic response problems of complex structures, including spinning ones. The aerodynamic module, AERO, enables computation of unsteady aerodynamic forces for both subsonic and supersonic flow for subsequent flutter and divergence analysis of the structure. The associated aeroservoelastic analysis module, ASE, effects aero-structural-control stability analysis yielding frequency responses as well as damping characteristics of the structure. The program is written in standard FORTRAN to run on a wide variety of computers. Extensive graphics, preprocessing, and postprocessing routines are also available pertaining to a number of terminals.

Gupta, Kajal K.↗

NASA NTRS1988

Error indicators and accuracy improvements of finite element solutions

Practical and reliable estimators of the discretization errors in engineering problems are developed. Error indicators for identifying the regions or elements of the solution domain which are likely to have the largest discretization errors are presented, and a simple computational procedure for improving the accuracy of the finite element solutions for shell problems is given. The similarities between the proposed procedure and a preconditioned conjugate gradient (PCG) technique are identified and exploited to generate pointwise error indicators from the PCG technique. Numerical examples in the linear static analysis of shells are presented.

Noor, Ahmed K.↗

NASA NTRS1985

Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies

The particular problems investigated in the present study arise from nonlinear transient heat conduction. One of two types of nonlinearities considered is related to a material temperature dependence which is frequently needed to accurately model behavior over the range of temperature of engineering interest. The second nonlinearity is introduced by radiation boundary conditions. The finite element equations arising from the solution of nonlinear transient heat conduction problems are formulated. The finite element matrix equations are temporally discretized, and a nonlinear iterative solution algorithm is proposed. Algorithms for solving the linear problem are discussed, taking into account the form of the matrix equations, Gaussian elimination, cost, and iterative techniques. Attention is also given to approximate factorization, implementational aspects, and numerical results.

Winget, J. M.↗

NASA NTRS1993

Variational Trajectory Optimization Tool Set: Technical description and user's manual

The algorithms that comprise the Variational Trajectory Optimization Tool Set (VTOTS) package are briefly described. The VTOTS is a software package for solving nonlinear constrained optimal control problems from a wide range of engineering and scientific disciplines. The VTOTS package was specifically designed to minimize the amount of user programming; in fact, for problems that may be expressed in terms of analytical functions, the user needs only to define the problem in terms of symbolic variables. This version of the VTOTS does not support tabular data; thus, problems must be expressed in terms of analytical functions. The VTOTS package consists of two methods for solving nonlinear optimal control problems: a time-domain finite-element algorithm and a multiple shooting algorithm. These two algorithms, under the VTOTS package, may be run independently or jointly. The finite-element algorithm generates approximate solutions, whereas the shooting algorithm provides a more accurate solution to the optimization problem. A user's manual, some examples with results, and a brief description of the individual subroutines are included.

Bless, Robert R.↗