Search NASA⌕ Search

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.

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

At least 307 records · Page 17

On local total strain redistribution using a simplified cyclic inelastic analysis based on an elastic solution

Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction purposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure has been found to predict stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load cycled problems. This study derived and incorporated Neuber type corrections in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was exercised on a mechanically load cycled benchmark notched plate problem. Excellent agreement was found between the predicted material response and nonlinear finite element solutions for the problem. The simplified analysis computer program used 0.3 percent of the CPU time required for a nonlinear finite element analysis.

Hwang, S. Y.↗

On local total strain redistribution using a simplified cyclic inelastic analysis based on an elastic solution

Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction purposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure has been found to predict stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load cycled problems. This study derived and incorporated Neuber type corrections in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was exercised on a mechanically load cycled benchmark notched plate problem. Excellent agreement was found between the predicted material response and nonlinear finite element solutions for the problem. The simplified analysis computer program used 0.3 percent of the CPU time required for a nonlinear finite element analysis.

Hwang, S. Y.↗

Solution-adaptive finite element method in computational fracture mechanics

Some recent results obtained using solution-adaptive finite element method in linear elastic two-dimensional fracture mechanics problems are presented. The focus is on the basic issue of adaptive finite element method for validating the applications of new methodology to fracture mechanics problems by computing demonstration problems and comparing the stress intensity factors to analytical results.

Min, J. B.↗

On the performance of explicit and implicit algorithms for transient thermal analysis

The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite element models of the configurations are discribed. These problems include a space shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system and a model of the space shuttle orbiter wing. Calculations were carried out using the SPAR finite element program, the MITAS lumped parameter program and a special purpose finite element program incorporating the GEAR algorithms. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are stiff. Careful attention to modeling detail such as avoiding thin or short high-conducting elements can sometimes reduce the stiffness to the extent that explicit methods become advantageous.

Adelman, H. M.↗

Micromechanical predictions of crack initiation, propagation and crack growth resistance in boron/aluminum composites

An elastoplastic, axisymmetric finite element model has been used to predict the initiation and propagation of a crack in a composite model consisting of a single broken boron fiber embedded in an annular sheath of aluminum matrix. The accuracy of the axisymmetric finite element model for crack problems has been established by solving the classical problem of a penny-shaped crack in a thick cylindrical rod under axial tension. Also, the stress intensity factors predicted by the present numerical model are compared with continuum results. A constant displacement boundary condition applied during an increment of crack growth permits a substantial amount of stable crack growth in the matrix material. The concept of Crack Growth Resistance Curves (KR-curves) has been used to determine the point of crack instability

Mahishi, J. M.↗

Propagating plane harmonic waves through finite length plates of variable thickness using finite element techniques

An analysis is given using finite element techniques which addresses the propagaton of a uniform incident pressure wave through a finite diameter axisymmetric tapered plate immersed in a fluid. The approach utilized in developing a finite element solution to this problem is based upon a technique for axisymmetric fluid structure interaction problems. The problem addressed is that of a 10 inch diameter axisymmetric fixed plate totally immersed in a fluid. The plate increases in thickness from approximately 0.01 inches thick at the center to 0.421 inches thick at a radius of 5 inches. Against each face of the tapered plate a cylindrical fluid volume was represented extending five wavelengths off the plate in the axial direction. The outer boundary of the fluid and plate regions were represented as a rigid encasement cylinder as was nearly the case in the physical problem. The primary objective of the analysis is to determine the form of the transmitted pressure distribution on the downstream side of the plate.

Clark, J. H.↗

Programming substructure computations for elliptic problems on the CHiP system

A number of studies have been conducted with the aim to apply parallel computation to problems associated with solving finite element equations arising in structural mechanics and fluid dynamics. These studies have provided many important results. The present investigation is concerned with a set of experiments designed to test two ideas, including configurability and substructuring. The considered algorithms and tests are intended for implementation on the Configurable, Highly Parallel (CHiP) family of architecture described by Snyder (1982). The ChiP computer is composed of homogeneous processing elements (PEs) placed at regular intervals in a lattice of programmable switches. Two examples of the role of configurability and substructuring for simple iterative algorithms are considered, giving attention to conjugate gradient iterations, and tridiagonal systems of equations.

Gannon, D.↗

Supercomputer implementation of finite element algorithms for high speed compressible flows

Prediction of compressible flow phenomena using the finite element method is of recent origin and considerable interest. Two shock capturing finite element formulations for high speed compressible flows are described. A Taylor-Galerkin formulation uses a Taylor series expansion in time coupled with a Galerkin weighted residual statement. The Taylor-Galerkin algorithms use explicit artificial dissipation, and the performance of three dissipation models are compared. A Petrov-Galerkin algorithm has as its basis the concepts of streamline upwinding. Vectorization strategies are developed to implement the finite element formulations on the NASA Langley VPS-32. The vectorization scheme results in finite element programs that use vectors of length of the order of the number of nodes or elements. The use of the vectorization procedure speeds up processing rates by over two orders of magnitude. The Taylor-Galerkin and Petrov-Galerkin algorithms are evaluated for 2D inviscid flows on criteria such as solution accuracy, shock resolution, computational speed and storage requirements. The convergence rates for both algorithms are enhanced by local time-stepping schemes. Extension of the vectorization procedure for predicting 2D viscous and 3D inviscid flows are demonstrated. Conclusions are drawn regarding the applicability of the finite element procedures for realistic problems that require hundreds of thousands of nodes.

Thornton, E. A.↗

A feasibility study on the numerical modeling of interior noise fields

New developments in electronics and computer control have made active methods of noise suppression feasible. To date, however, active methods of noise control have been applied principally in geometrically uniform spaces with localized sound sources. The purpose of this paper is to demonstrate how finite element techniques may be used to model sound fields and active control inside geometrically complex spaces with non-uniform acoustic boundaries and multiple spatially separated sources. A closed-form solution is derived for the case of a one-dimensional tube with a vibrating piston at each end. Results are compared with finite element solutions and the problem is then generalized to two-dimensions with non-uniform geometry and non-uniform acoustic boundaries.

Abrahamson, A. L.↗

A hybrid nonlinear programming method for design optimization

Solutions to engineering design problems formulated as nonlinear programming (NLP) problems usually require the use of more than one optimization technique. Moreover, the interaction between the user (analysis/synthesis) program and the NLP system can lead to interface, scaling, or convergence problems. An NLP solution system is presented that seeks to solve these problems by providing a programming system to ease the user-system interface. A simple set of rules is used to select an optimization technique or to switch from one technique to another in an attempt to detect, diagnose, and solve some potential problems. Numerical examples involving finite element based optimal design of space trusses and rotor bearing systems are used to illustrate the applicability of the proposed methodology.

Rajan, S. D.↗

A spectral overlay method for dissimilar materials

The spectral overlay method is applied to examine the stress distribution along the interface between the two different materials. The essential feature of the spectral overlay method is that the high resolution of localized steep gradients can be achieved by overlaying a spectral interpolant on a standard finite element mesh. In order to evaluate this method and compare it with conventional finite element method, a test problem with high gradients in Poisson equation has been devised. The results show that this method is very powerful in capturing the shape of the gradient field and its peak value. Solutions for two problems with interfaces between dissimilar materials are then given to illustrate the effectiveness of this method.

Belytschko, T.↗

A hybrid formulation of a component mode synthesis method

Component mode synthesis is a substructuring technique frequently employed in structural dynamics. In this method, a given structure is subdivided into components or substructures, each of which is analyzed independently for natural frequencies and for mode shapes. The substructure mode shapes are then assembled to give displacement shapes or load patterns of the original structure. An analytical justification of the basic concept is presented using spectral decompositions, and a variant substructuring approach where intersubstructure continuity is enforced in a weak form is derived. This leads to a hybrid formulation of the basic method which is particularly suitable for assembling heterogeneous substructures and analyzing nonconforming and incompatible finite element substructure models. For problems where both the basic and hybrid methods are applicable, the hybrid variant can be computationally more advantageous.

Farhat, Charbel↗

Discretized energy minimization in a wave guide with point sources

An anti-noise problem on a finite time interval is solved by minimization of a quadratic functional on the Hilbert space of square integrable controls. To this end, the one-dimensional wave equation with point sources and pointwise reflecting boundary conditions is decomposed into a system for the two propagating components of waves. Wellposedness of this system is proved for a class of data that includes piecewise linear initial conditions and piecewise constant forcing functions. It is shown that for such data the optimal piecewise constant control is the solution of a sparse linear system. Methods for its computational treatment are presented as well as examples of their applicability. The convergence of discrete approximations to the general optimization problem is demonstrated by finite element methods.

Propst, G.↗

An aspect ration bound for triangulating a d-grid cut by a hyperplane

We consider the problem of triangulating a d-dimensional uniform grid of d-cubes that is cut by a k-dimensional affine subspace. The goal is to obtain a triangulation with bounded aspect ratio. To achieve this goal, we allow some of the box faces near the affine subspace to be displaced. This problem has applications to finite element mesh generation. For general d and k, the bound on aspect ratio that we attain is double-exponential in d. For the important special case of d = 3, the aspect ratio bound is small enough that the technique is useful in practice.

Mitchell, Scott A.↗

Plasticity Tool for Predicting Shear Nonlinearity of Unidirectional Laminates Under Multiaxial Loading

This study implements a plasticity tool to predict the nonlinear shear behavior of unidirectional composite laminates under multiaxial loadings, with an intent to further develop the tool for use in composite progressive damage analysis. The steps for developing the plasticity tool include establishing a general quadratic yield function, deriving the incremental elasto-plastic stress-strain relations using the yield function with associated flow rule, and integrating the elasto-plastic stress-strain relations with a modified Euler method and a substepping scheme. Micromechanics analyses are performed to obtain normal and shear stress-strain curves that are used in determining the plasticity parameters of the yield function. By analyzing a micromechanics model, a virtual testing approach is used to replace costly experimental tests for obtaining stress-strain responses of composites under various loadings. The predicted elastic moduli and Poisson's ratios are in good agreement with experimental data. The substepping scheme for integrating the elasto-plastic stress-strain relations is suitable for working with displacement-based finite element codes. An illustration problem is solved to show that the plasticity tool can predict the nonlinear shear behavior for a unidirectional laminate subjected to multiaxial loadings.

Wang, John T.↗

Modeling holes and voids in three dimensions using a single element within the Extended Finite Element framework

This paper introduces a highly efficient computational approach using the extended finite element method to model three-dimensional mechanics of holes and voids. The defect is encapsulated within a single element with a tailored exponential enrichment function that captures the behavior near the discontinuity boundary and can account for the pressure acting on the void’s surface through an equivalent nodal force vector. The method’s results were compared against analytical solutions and conventional finite element method results for problems involving single and multiple holes. The comparison demonstrated that the method provides accurate stress concentration estimates while adhering to computational constraints.

Adnan Shahriar↗

Finite elements of nonlinear continua.

The finite element method is extended to a broad class of practical nonlinear problems, treating both theory and applications from a general and unifying point of view. The thermomechanical principles of continuous media and the properties of the finite element method are outlined, and are brought together to produce discrete physical models of nonlinear continua. The mathematical properties of the models are analyzed, and the numerical solution of the equations governing the discrete models is examined. The application of the models to nonlinear problems in finite elasticity, viscoelasticity, heat conduction, and thermoviscoelasticity is discussed. Other specific topics include the topological properties of finite element models, applications to linear and nonlinear boundary value problems, convergence, continuum thermodynamics, finite elasticity, solutions to nonlinear partial differential equations, and discrete models of the nonlinear thermomechanical behavior of dissipative media.

Oden, J. T.↗