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 595 records · Page 33

Design of optical mirror structures

The structural requirements for large optical telescope mirrors was studied with a particular emphasis placed on the three-meter Large Space Telescope primary mirror. Analysis approaches through finite element methods were evaluated with the testing and verification of a number of element types suitable for particular mirror loadings and configurations. The environmental conditions that a mirror will experience were defined and a candidate list of suitable mirror materials with their properties compiled. The relation of the mirror mechanical behavior to the optical performance is discussed and a number of suitable design criteria are proposed and implemented. A general outline of a systematic method to obtain the best structure for the three-meter diffraction-limited system is outlined. Finite element programs, using the STRUDL 2 analysis system, were written for specific mirror structures encompassing all types of active and passive mirror designs. Parametric studies on support locations, effects of shear deformation, diameter to thickness ratios, lightweight and sandwich mirror configurations, and thin shell active mirror needs were performed.

Soosaar, K.↗

Quantification of Saint-Venant's principle for a general prismatic member

A method to quantify Saint-Venant's principle for general prismatic members is described. Finite elements are used to model the member cross-section to allow the analysis of arbitrary geometries, material properties, and nonhomogeneities. The variation in the structural response along the member's length, away from the loaded end, is modelled with an exponential decay which has been established in the literature to be the correct form. Many analytic results exist in the literature for special cases. Because of its ability to handle arbitrary cross sectional geometries, this semi-analytic finite element method can provide a firm footing for the utilization of Saint-Venant's principle for whole new classes of problems. Results for several typical problems are given.

Goetschel, D. B.↗

ENSAERO - A multidisciplinary program for fluid/structural interaction studies of aerospace vehicles

A computational procedure is developed that uses a moving zonal grid concept to model complex flexible aerospace vehicles. The Euler/Navier-Stokes equations are used to model the flow, and computations are made using efficient methods based on both central and upwind schemes. The structure is represented by a finite element method which can model general aerospace vehicles. Provisions are made to accommodate other disciplines such as controls and thermal loads. The code is capable of computing unsteady flows on flexible wings with vortical flows. Adaptation of this procedure for parallel processing and validation for complete aerospace configurations is in progress.

Guruswamy, G. P.↗

Computational Design of Materials: Planetary Entry to Electric Aircraft and Beyond

NASA's projects and missions push the bounds of what is possible. To support the agency's work, materials development must stay on the cutting edge in order to keep pace. Today, researchers at NASA Ames Research Center perform multiscale modeling to aid the development of new materials and provide insight into existing ones. Multiscale modeling enables researchers to determine micro- and macroscale properties by connecting computational methods ranging from the atomic level (density functional theory, molecular dynamics) to the macroscale (finite element method). The output of one level is passed on as input to the next level, creating a powerful predictive model.

Materials Design↗

A method for the geometrically nonlinear analysis of compressively loaded prismatic composite structures

A method was developed for the geometrically nonlinear analysis of the static response of thin-walled stiffened composite structures loaded in uniaxial or biaxial compression. The method is applicable to arbitrary prismatic configurations composed of linked plate strips, such as stiffened panels and thin-walled columns. The longitudinal ends of the structure are assumed to be simply supported, and geometric shape imperfections can be modeled. The method can predict the nonlinear phenomena of postbuckling strength and imperfection sensitivity which are exhibited by some buckling-dominated structures. The method is computer-based and is semi-analytic in nature, making it computationally economical in comparison to finite element methods. The method uses a perturbation approach based on the use of a series of buckling mode shapes to represent displacement contributions associated with nonlinear response. Displacement contributions which are of second order in the model amplitudes are incorported in addition to the buckling mode shapes. The principle of virtual work is applied using a finite basis of buckling modes, and terms through the third order in the model amplitudes are retained. A set of cubic nonlinear algebraic equations are obtained, from which approximate equilibrium solutions are determined. Buckling mode shapes for the general class of structure are obtained using the VIPASA analysis code within the PASCO stiffened-panel design code. Thus, subject to some additional restrictions in loading and plate anisotropy, structures which can be modeled with respect to buckling behavior by VIPASA can be analyzed with respect to nonlinear response using the new method. Results obtained using the method are compared with both experimental and analytical results in the literature. The configurations investigated include several different unstiffened and blade-stiffening panel configurations, featuring both homogeneous, isotropic materials, and laminated composite material.

Stoll, Frederick↗

Finite element discrete ordinates method for radiative transfer in non-rotationally invariant scattering media - Application to the leaf canopy problem

A finite element discrete ordinates method for solving the radiative transfer equation in nonrotationally invariant scattering media has been applied to the lead-canopy problem, and results are presented on the cross sections and the reflection functions. The method is based on a unique implementation of the Galerkin integral law formulation of the transport equation. For both near-normal and grazing incidences, the transfer functions of leaf canopies are found to be strongly anisotropic, with relatively more scattered flux in the vertical directions. It is suggested that the assumption of isotropic scattering in leaf canopies is not valid.

Myneni, Ranga B.↗

Curvilinear Displacement Transfer Functions for Deformed Shape Predictions of Curved Structures Using Distributed Surface Strains

Curvilinear Displacement Transfer Functions were formulated for deformed shape predictions of different curved structures using surface strains. The embedded curved beam (depth-wise cross section of a curved structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. The discrete approach enabled piecewise integrations of a curvature-strain differential equation for the embedded curved beam to yield closed-form Curvilinear Displacement Transfer Functions, which are written in terms of embedded curved-beam geometrical parameters and surface strains. By inputting the surface strain data, the Curvilinear Displacement Transfer Functions can transform surface strains into deflections along each embedded curved beam for mapping out the overall structural deformed shapes. The finite-element method was used to analytically generate the surface strains of the curved beams. The deformed shape prediction accuracies were then determined by comparing the theoretical deflections with the finite-element-generated deflections, which were used as yardsticks. By introducing the correction factors in simple mathematical forms, the Curvilinear Displacement Transfer Functions can be quite accurate for shape predictions of different curved-beam structures ranging from limit case of straight beam up to semicircular curved beam.

Ko, William L.↗

Simulating Fatigue Crack Growth in Spiral Bevel Pinion

This project investigates computational modeling of fatigue crack growth in spiral bevel gears. Current work is a continuation of the previous efforts made to use the Boundary Element Method (BEM) to simulate tooth-bending fatigue failure in spiral bevel gears. This report summarizes new results predicting crack trajectory and fatigue life for a spiral bevel pinion using the Finite Element Method (FEM). Predicting crack trajectories is important in determining the failure mode of a gear. Cracks propagating through the rim may result in catastrophic failure, whereas the gear may remain intact if one tooth fails and this may allow for early detection of failure. Being able to predict crack trajectories is insightful for the designer. However, predicting growth of three-dimensional arbitrary cracks is complicated due to the difficulty of creating three-dimensional models, the computing power required, and absence of closed- form solutions of the problem. Another focus of this project was performing three-dimensional contact analysis of a spiral bevel gear set incorporating cracks. These analyses were significant in determining the influence of change of tooth flexibility due to crack growth on the magnitude and location of contact loads. This is an important concern since change in contact loads might lead to differences in SIFs and therefore result in alteration of the crack trajectory. Contact analyses performed in this report showed the expected trend of decreasing tooth loads carried by the cracked tooth with increasing crack length. Decrease in tooth loads lead to differences between SIFs extracted from finite element contact analysis and finite element analysis with Hertz contact loads. This effect became more pronounced as the crack grew.

Ural, Ani↗

ACCESS 3. Approximation concepts code for efficient structural synthesis: User's guide

A user's guide is presented for ACCESS-3, a research oriented program which combines dual methods and a collection of approximation concepts to achieve excellent efficiency in structural synthesis. The finite element method is used for structural analysis and dual algorithms of mathematical programming are applied in the design optimization procedure. This program retains all of the ACCESS-2 capabilities and the data preparation formats are fully compatible. Four distinct optimizer options were added: interior point penalty function method (NEWSUMT); second order primal projection method (PRIMAL2); second order Newton-type dual method (DUAL2); and first order gradient projection-type dual method (DUAL1). A pure discrete and mixed continuous-discrete design variable capability, and zero order approximation of the stress constraints are also included.

Fleury, C.↗

Parallel iterative methods for sparse linear and nonlinear equations

As three-dimensional models are gaining importance, iterative methods will become almost mandatory. Among these, preconditioned Krylov subspace methods have been viewed as the most efficient and reliable, when solving linear as well as nonlinear systems of equations. There has been several different approaches taken to adapt iterative methods for supercomputers. Some of these approaches are discussed and the methods that deal more specifically with general unstructured sparse matrices, such as those arising from finite element methods, are emphasized.

Saad, Youcef↗

Simultaneous Aerodynamic and Structural Design Optimization (SASDO) for a 3-D Wing

The formulation and implementation of an optimization method called Simultaneous Aerodynamic and Structural Design Optimization (SASDO) is shown as an extension of the Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) method. It is extended by the inclusion of structure element sizing parameters as design variables and Finite Element Method (FEM) analysis responses as constraints. The method aims to reduce the computational expense. incurred in performing shape and sizing optimization using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis, FEM structural analysis and sensitivity analysis tools. SASDO is applied to a simple. isolated, 3-D wing in inviscid flow. Results show that the method finds the saine local optimum as a conventional optimization method with some reduction in the computational cost and without significant modifications; to the analysis tools.

Gumbert, Clyde R.↗

Space Shuttle Main Engine structural analysis and data reduction/evaluation. Volume 6: Primary nozzle diffuser analysis

The primary nozzle diffuser routes fuel from the main fuel valve on the Space Shuttle Main Engine (SSME) to the nozzle coolant inlet mainfold, main combustion chamber coolant inlet mainfold, chamber coolant valve, and the augmented spark igniters. The diffuser also includes the fuel system purge check valve connection. A static stress analysis was performed on the diffuser because no detailed analysis was done on this part in the past. Structural concerns were in the area of the welds because approximately 10 percent are in areas inaccessible by X-ray testing devices. Flow dynamics and thermodynamics were not included in the analysis load case. Constant internal pressure at maximum SSME power was used instead. A three-dimensional, finite element method was generated using ANSYS version 4.3A on the Lockheed VAX 11/785 computer to perform the stress computations. IDEAS Supertab on a Sun 3/60 computer was used to create the finite element model. Rocketdyne drawing number RS009156 was used for the model interpretation. The flight diffuser is denoted as -101. A description of the model, boundary conditions/load case, material properties, structural analysis/results, and a summary are included for documentation.

Foley, Michael J.↗

Geometrically nonlinear analysis of laminated elastic structures

Laminated composite plates and shells that can be used to model automobile bodies, aircraft wings and fuselages, and pressure vessels among many other were analyzed. The finite element method, a numerical technique for engineering analysis of structures, is used to model the geometry and approximate the solution. Various alternative formulations for analyzing laminated plates and shells are developed and their finite element models are tested for accuracy and economy in computation. These include the shear deformation laminate theory and degenerated 3-D elasticity theory for laminates.

Reddy, J. N.↗

Combined, nonlinear aerodynamic and structural method for the aeroelastic design of a three-dimensional wing in supersonic flow

An iterative procedure for the static aeroelastic design of a flexible wing at supersonic speeds has been developed. The procedure combines a nonlinear, full-potential solver (NCOREL) with an equivalent plate structural analysis method. The NCOREL method yields significantly improved aerodynamic estimates compared to linear theory. The equivalent plate structural analysis method demonstrates an order of magnitude reduction in computer memory and execution time compared to finite-element methods. A highly swept wing is analyzed at high lift using this aeroelastic procedure. The results indicate that the wing deforms favorably due to aerodynamic loading and, consequently, that the inviscid drag levels do not vary at the required lift coefficient although the angle of attack varies significantly. A sensitivity analysis of the type required for optimization studies was also performed with the aeroelastic design procedure.

Pittman, J. L.↗

A finite element conjugate gradient FFT method for scattering

An extension of a two dimensional formulation is presented for a three dimensional body of revolution. With the introduction of a Fourier expansion of the vector electric and magnetic fields, a coupled two dimensional system is generated and solved via the finite element method. An exact boundary condition is employed to terminate the mesh and the fast fourier transformation (FFT) is used to evaluate the boundary integrals for low O(n) memory demand when an iterative solution algorithm is used. By virtue of the finite element method, the algorithm is applicable to structures of arbitrary material composition. Several improvements to the two dimensional algorithm are also described. These include: (1) modifications for terminating the mesh at circular boundaries without distorting the convolutionality of the boundary integrals; (2) the development of nonproprietary mesh generation routines for two dimensional applications; (3) the development of preprocessors for interfacing SDRC IDEAS with the main algorithm; and (4) the development of post-processing algorithms based on the public domain package GRAFIC to generate two and three dimensional gray level and color field maps.

Collins, Jeffery D.↗

Implicit Extrapolation Methods for Variable Coefficient Problems

Implicit extrapolation methods for the solution of partial differential equations are based on applying the extrapolation principle indirectly. Multigrid tau-extrapolation is a special case of this idea. In the context of multilevel finite element methods, an algorithm of this type can be used to raise the approximation order, even when the meshes are nonuniform or locally refined. Here previous results are generalized to the variable coefficient case and thus become applicable for nonlinear problems. The implicit extrapolation multigrid algorithm converges to the solution of a higher order finite element system. This is obtained without explicitly constructing higher order stiffness matrices but by applying extrapolation in a natural form within the algorithm. The algorithm requires only a small change of a basic low order multigrid method.

Jung, M.↗

Use of finite element techniques in the determination of the acoustic properties of turbofan inlets

This paper describes the application of the finite-element method in combination with Galerkin's method in the determination of the acoustic properties of turbofan inlets containing high-subsonic-Mach-number flows. An approximate solution to the steady inviscid flow field is obtained using an integral method for calculating the potential flow field in the inlet with a correction to account for compressibility effects. The accuracy of the finite-element technique in predicting the acoustic properties of annular ducts has been checked by comparison with available analytical solutions for the problems of planeand spinning-wave propagation through a hard-walled annular duct with a constant mean flow. Results are presented comparing low-frequency plane-wave propagation through a hard-walled turbofan inlet containing a one-dimensional flow with the same inlet containing a fully two-dimensional axisymmetric mean flow. It is shown that when one-dimensional mean flow is assumed to exist in the duct, the plane wave propagates with relatively little distortion. However, propagation of a plane wave through the fully two-dimensional flow field in the inlet produces severe distortion due to the excitation of higher-order modes.

Sigman, R. K.↗