Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element methods”

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 685 records · Page 38

Probability of failure and risk assessment of propulsion structural components

The probabilistic structural analysis method (PSAM) was developed to analyze the effects of fluctuating loads, variable material properties, and uncertain analytical models especially for high performance structures such as the Space Shuttle Main Engine turbopump blades. Risk is calculated after expensive service experience. However, probabilistic structural analysis provides a rational alternative method to quantify uncertainties in the structural performance and durability. NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) was developed as a probabilistic structural analysis computer code which integrates finite element methods and reliability algorithms, capable to predicting the probability distributions of structural response variables such as stress, displacement, natural frequencies, and buckling loads. This computer code is detailed.

Shiao, Michael C.↗

Cyclic-stress analysis of notches for supersonic transport conditions

The feasibility of using the finite element method to account for the effects of cyclic load and temperature on local stresses and strains at a notch was demonstrated. The behavior of a notched titanium panel was studied under variable loads and temperatures representative of flight conditions for the lower wing surface of a Supersonic Transport (SST). The analysis was performed with the use of the BOPACE finite-element computer program which provides capability to determine high temperature and large viscoplastic effects caused by cyclic thermal and mechanical loads. The analysis involves the development of the finite-element model as well as determination of the structural behavior of the notched panel. Results are presented for twelve SST flights comprised of five different load-temperature cycles. The results show the approach is feasible, but material response to cyclic loads, temperatures, and hold times requires improved understanding to allow proper modeling of the material.

Source record↗

The truss-beam and space technology

A method is presented herein to determine the dynamic characteristics of truss-beam systems for both planar and space structures. Results from this approach are compared with those of the finite element method and test results. The truss-beam model is based on the assumptions of negligible web contribution, i.e., very large flexible structures. Future experimental evaluation of the plane truss-system is planned at North Carolina Agricultural and Technical State University and this experimental program is described.

Abu-Saba, Elias G.↗

Finite elements and fluid dynamics

Difficulties concerning a use of the finite element method in the solution of the nonlinear equations of fluid dynamics are partly related to various 'hidden' instabilities which often arise in fluid calculations. The instabilities are typically due to boundary effects or nonlinearities. It is shown that in certain cases these instabilities can be avoided if certain conservation laws are satisfied, and that the latter are often intimately related to finite elements.

Fix, G.↗

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗

Application of the extended boundary condition method to Monte Carlo simulations of scattering of waves by two-dimensional random rough surfaces

The extended boundary condition method is applied to Monte Carlo simulations of two-dimensional random rough surface scattering. The numerical results are compared with one-dimensional random rough surfaces obtained from the finite-element method. It is found that the mean scattered intensity from two-dimensional rough surfaces differs from that of one dimension for rough surfaces with large slopes.

Tsang, L.↗

An AWE Implementation for Electromagnetic Analysis

Although full wave electromagnetic systems are large and cumbersome to solve, typically only a few parameters, such as input impedance, S parameters, and far field pattern, are needed by the designer or analyst. A reduced order modeling of these parameters is therefore an important consideration in minimizing the the CPU requirements. The Asymptotic Waveform Evaluation (AWE) method is one approach to construct a reduced order model of the input impedance or other useful electromagnetic parameters. We demonstrate its application and validity when used in conjunction with the finite element method to simulate full wave electromagnetic problems.

Gong, Jian↗

Singularity Preserving Numerical Methods for Boundary Integral Equations

In the past twelve months (May 8, 1995 - May 8, 1996), under the cooperative agreement with Division of Multidisciplinary Optimization at NASA Langley, we have accomplished the following five projects: a note on the finite element method with singular basis functions; numerical quadrature for weakly singular integrals; superconvergence of degenerate kernel method; superconvergence of the iterated collocation method for Hammersteion equations; and singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel. This final report consists of five papers describing these projects. Each project is preceeded by a brief abstract.

Kaneko, Hideaki↗

Vibrations of cantilevered shallow cylindrical shells of rectangular planform

A cantilevered, shallow shell of circular cylindrical curvature and rectangular planform exhibits free vibration behavior which differs considerably from that of a cantilevered beam or of a flat plate. Some numerical results can be found for the problem in the previously published literature, mainly obtained by using various finite element methods. The present paper is the first definitive study of the problem, presenting accurate non-dimensional frequency parameters for wide ranges of aspect ratio, shallowness ratio and thickness ratio. The analysis is based upon shallow shell theory. Numerical results are obtained by using the Ritz method, with algebraic polynomial trial functions for the displacements. Convergence is investigated, with attention being given both to the number of terms taken for each co-ordinate direction and for each of the three components of displacement. Accuracy of the results is also established by comparison with finite element results for shallow shells and with other accurate flat plate solutions.

Leissa, A. W.↗

Calculation of exact vibration modes for plane grillages by the dynamic stiffness method

A dynamic stiffness method is developed for the calculation of the exact modal parameters for plane grillages which consist of straight and uniform beams with coincident elastic and inertial axes. Elementary bending-torsion beam theory is utilized, and bending translation is restricted to one direction. The exact bending-torsion dynamic stiffness matrix is obtained for a straight and uniform beam element with coincident elastic and inertial axes. The element stiffness matrices are assembled using the standard procedure of the static stiffness method to form the dynamic stiffness matrix of the complete grillage. The exact natural frequencies, mode shapes, and generalized masses of the grillage are then calculated by solving a nonlinear eigenvalue problem based on the dynamic stiffness matrix. The exact modal solutions for an example grillage are calculated and compared with the approximate solutions obtained by using the finite element method.

Hallauer, W. L., Jr.↗

Steady, Oscillatory and Unsteady, Subsonic and Supersonic Aerodynamics (SOUSSA) for complex aircraft configurations

The Green's function method and the computer program SOUSSA (Steady Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed. The Green's function method is applied to the fully unsteady potential equation yielding an integro-differential-delay equation. This equation is approximated by a set of differential-delay equations in time using the finite element method. The Laplace transform is used to yield a matrix relating the velocity potential to the normal wash. The matrix of the generalized aerodynamic forces is obtained by premultiplying and postmultiplying the matrices relating generalized forces to the potential and the normal wash by the generalized coordinates. The program SOUSSA is compared with existing numerical results. Results indicate that the program is not only general, flexible, and easy to use, but also accurate and fast.

Morino, L.↗

Theoretical performance of foil journal bearings

A modified forward iteration approach for the coupled solution of foil bearings is presented. The method is used to predict the steady state theoretical performance of a journal type gas bearing constructed from an inextensible shell supported by an elastic foundation. Bending effects are treated as negligible. Finite element methods are used to predict both the foil deflections and the pressure distribution in the gas film.

Carpino, M.↗

Realistic Thermographic Simulation of Impact Damage with Quadrupole Method

Flash thermography has been shown to be an effective method for detection of delaminations in carbon fiber reinforced polymer (CFRP) composite structures. Improved understanding of thelimitations of the technique can be obtained by simulating the inspection process. Time domainfinite difference and finite element methods are well suited for such simulations but can be computationally intensive. An alternate method is to solve the Laplace transform of the heat equation, then invert the Laplace transform to produce a time domain response. Often this is referred to as the thermal quadrupole method and is well suited for simulating flash thermography. The quadrupole method has been used extensively for simulating one-dimensional heat flow in multilayer systems with contact resistances at the interfaces. It is also applicable for three dimensional configurations, in particular for simulation of changes in thermographic response due to delaminations in composites. This presentation discusses three-dimensional simulations with realistic shapes and contact resistances at the interfaces between layers. The results of these simulations are compared to thermography measurements on composites specimens with impact induced delaminations.

Thermography↗

A space radiation transport method development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest-order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard finite element method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 ms and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re-configurable computing and could be utilized in the final design as verification of the deterministic method optimized design. Published by Elsevier Ltd on behalf of COSPAR.

NASA Center JSC↗

Finite-element analysis of dynamic fracture

Applications of the finite element method to the two dimensional elastodynamics of cracked structures are presented. Stress intensity factors are computed for two problems involving stationary cracks. The first serves as a vehicle for discussing lumped-mass and consistent-mass characterizations of inertia. In the second problem, the behavior of a photoelastic dynamic tear test specimen is determined for the time prior to crack propagation. Some results of a finite element simulation of rapid crack propagation in an infinite body are discussed.

Aberson, J. A.↗

Moment method with isoparametric elements for three-dimensional anisotropic scatterers

A novel method for computing the frequency-domain electromagnetic fields scattered from, and penetrating into, arbitrarily shaped, three-dimensional, lossy, inhomogeneous anisotropic scatters is presented. The method is based on a general volume integrodifferential formulation of the scattering problem and consists of the numerical solution of the coupled integral equations by the moment method and point matching. The numerical model of the scatterer is obtained by parametric volume elements, and the basis functions used to represent the field within each element are the same used in the finite-element method. Element integration problems due to the singular kernel of the integral equations are treated in some detail. Numerical results for both the isotropic and the anisotropic spherical scatterer are presented, including comparisons with results obtained by different numerical methods for the isotropic cases considered. The capability of the numerical code presented to deal with cases where the material parameters of the scatterer are given by singular matrices is discussed for two particular examples.

Graglia, Roberto D.↗

Realistic Thermographic Simulation of Impact Damage with Quadrupole Method

Flash thermography has been shown to be an effective method for detection of delaminations in carbon fiber reinforced polymer (CFRP) composite structures. Improved understanding of the limitations of the technique can be obtained by simulating the inspection process. Time domain finite difference and finite element methods are well suited for such simulations but can be computationally intensive. An alternate method is to solve the Laplace transform of the heat equation, then invert the Laplace transform to produce a time domain response. Often this is referred to as the thermal quadrupole method and is well suited for simulating flash thermography. The quadrupole method has been used extensively for simulating one-dimensional heat flow in multilayer systems with contact resistances at the interfaces. It is also applicable for three dimensional configurations, in particular for simulation of changes in thermographic response due to delaminations in composites. This presentation discusses three-dimensional simulations with realistic shapes and contact resistances at the interfaces between layers. The results of these simulations are compared to thermography measurements on composites specimens with impact induced delaminations.

thermography↗

Reduction methods for nonlinear steady-state thermal analysis

Hybrid analysis techniques and a problem adaptive computational algorithm are presented for predicting the nonlinear steady state temperature distribution in structures and solids. In the proposed techniques, the structure is discretized by using the finite element method. The vector of nodal temperatures is then expressed as a linear combination of a small number of global temperature modes (or basis vectors), and the Bubnov Galerkin technique is used to compute the amplitudes of the global modes. The global temperature modes (or basis vectors) are chosen to include the various order derivatives of the nodal temperature vector with respect to preselected path parameter(s). The potential of the proposed reduction methods for solution of large scale, nonlinear thermal problems is discussed and the effectiveness of the methods is demonstrated by means of numerical examples, including steady state conduction, convection, and radiation modes of heat transfer.

Noor, A. K.↗