Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 811 records · Page 45

Analysis of aircraft tires via semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell.

Noor, Ahmed K.↗

Study of the available finite element software packages at KSC

The interaction among the three finite element software packages, SDRCI/I-DEAS, MSC/NASTRAN, and I/FEM, used at NASA, Kennedy Space Center is addressed. The procedures for using more than one of these application software packages to model and analyze a structure design are discussed. Design and stress analysis of a solid rocket booster fixture is illustrated by using four different combinations of the three software packages. Their results are compared and show small yet acceptable differences.

Lu, Chu-Ho↗

Structural behavior of scientific balloons - Finite element simulation and verification

An off-the-shelf nonlinear finite element code was used to analyze fully inflated scientific balloons. The thin balloon film was modeled by shell bending elements. Numerical difficulties caused by insignificant bending stiffness terms were overcome by introducing some artificial bending stiffness. This approximation is justified by the fact that in thin shells with nonzero Gaussian curvature the membrane solution component is essentially independent of the bending solution component. Perturbation of the coveraged solution by increasing the bending stiffness by a full decade verified this assertion. This analytical approach was experimentally verified. As a result of this verification process it was discovered that the generally accepted linearly visco-elastic model for polyethylene film is inappropriate for a significant planar (as opposed to uniaxial) stress state. A linear elastic model presents a good approximation for planar stress states.

Schur, Willi W.↗

Adaptive finite element strategies for shell structures

The present paper extends existing finite element adaptive refinement (AR) techniques to shell structures, which have heretofore been neglected in the AR literature. Specific challenges in applying AR to shell structures include: (1) physical discontinuities (e.g., stiffener intersections); (2) boundary layers; (3) sensitivity to geometric imperfections; (4) the sensitivity of most shell elements to mesh distortion, constraint definition and/or thinness; and (5) intrinsic geometric nonlinearity. All of these challenges but (5) are addressed here.

Stanley, G.↗

Adaptive explicit and implicit finite element methods for transient thermal analysis

The application of adaptive finite element methods to the solution of transient heat conduction problems in two dimensions is investigated. The computational domain is represented by an unstructured assembly of linear triangular elements and the mesh adaptation is achieved by local regeneration of the grid, using an error estimation procedure coupled to an automatic triangular mesh generator. Two alternative solution procedures are considered. In the first procedure, the solution is advanced by explicit timestepping, with domain decomposition being used to improve the computational efficiency of the method. In the second procedure, an algorithm for constructing continuous lines which pass only once through each node of the mesh is employed. The lines are used as the basis of a fully implicit method, in which the equation system is solved by line relaxation using a block tridiagonal equation solver. The numerical performance of the two procedures is compared for the analysis of a problem involving a moving heat source applied to a convectively cooled cylindrical leading edge.

Probert, E. J.↗

Finite-element time evolution operator for the anharmonic oscillator

The finite-element approach to lattice field theory is both highly accurate (relative errors approximately 1/N(exp 2), where N is the number of lattice points) and exactly unitary (in the sense that canonical commutation relations are exactly preserved at the lattice sites). In this talk I construct matrix elements for dynamical variables and for the time evolution operator for the anharmonic oscillator, for which the continuum Hamiltonian is H = p(exp 2)/2 + lambda q(exp 4)/4. Construction of such matrix elements does not require solving the implicit equations of motion. Low order approximations turn out to be extremely accurate. For example, the matrix element of the time evolution operator in the harmonic oscillator ground state gives a results for the anharmonic oscillator ground state energy accurate to better than 1 percent, while a two-state approximation reduces the error to less than 0.1 percent.

Milton, Kimball A.↗

Finite element analysis of the Space Shuttle 2.5-inch frangible nut

Finite element analysis of the Space Shuttle 2.5-inch frangible nut was conducted to improve understanding of the current design and proposed design changes to this explosively-actuated nut. The 2.5-inch frangible nut is used in two places to attach the aft end of the Space Shuttle Orbiter to the External Tank. Both 2.5-inch frangible nuts must function to complete safe separation. The 2.5-inch frangible nut contains two explosive boosters containing RDX explosive each capable of splitting the nut in half, on command from the Orbiter computers. To ensure separation, the boosters are designed to be redundant. The detonation of one booster is sufficient to split the nut in half. However, beginning in 1987 some production lots of 2.5-inch frangible nuts have demonstrated an inability to separate using only a single booster. The cause of the failure has been attributed to differences in the material properties and response of the Inconel 718 from which the 2.5-inch frangible nut is manufactured. Subsequent tests have resulted in design modifications of the boosters and frangible nut. Model development and initial analysis was conducted by Sandia National Laboratories (SNL) under funding from NASA Lyndon B. Johnson Space Center (NASA-JSC) starting in 1992. Modeling codes previously developed by SNL were transferred to NASA-JSC for further analysis on this and other devices. An explosive bolt with NASA Standard Detonator (NSD) charge, a 3/4-inch frangible nut, and the Super*Zip linear separation system are being modeled by NASA-JSC.

Darin N. McKinnis↗

A Least-Squares Finite Element Method for Electromagnetic Scattering Problems

The least-squares finite element method (LSFEM) is applied to electromagnetic scattering and radar cross section (RCS) calculations. In contrast to most existing numerical approaches, in which divergence-free constraints are omitted, the LSFF-M directly incorporates two divergence equations in the discretization process. The importance of including the divergence equations is demonstrated by showing that otherwise spurious solutions with large divergence occur near the scatterers. The LSFEM is based on unstructured grids and possesses full flexibility in handling complex geometry and local refinement Moreover, the LSFEM does not require any special handling, such as upwinding, staggered grids, artificial dissipation, flux-differencing, etc. Implicit time discretization is used and the scheme is unconditionally stable. By using a matrix-free iterative method, the computational cost and memory requirement for the present scheme is competitive with other approaches. The accuracy of the LSFEM is verified by several benchmark test problems.

Wu, Jie↗

Finite Element Modeling of Transient Thermography Inspection of Composite Materials

Several finite element models of defects such as debond and void have been developed for composite panels subjected to transient thermography inspection. Since the exact nature of the heat generated from the flash lamps is unknown, direct comparison between FEA and experimental results is not possible. However, some similarity of the results has been observed. The shape of the time curve that simulates the heat flux from the flash lamps has minimal effect on the temperature profiles. Double the number of flash lamps could increase the contrast of thermal image and define the shape of defect better.

Chu, Tsuchin Philip↗

Extension and Validation of a Hybrid Particle-Finite Element Method for Hypervelocity Impact Simulation

The hybrid particle-finite element method of Fahrenthold and Horban, developed for the simulation of hypervelocity impact problems, has been extended to include new formulations of the particle-element kinematics, additional constitutive models, and an improved numerical implementation. The extended formulation has been validated in three dimensional simulations of published impact experiments. The test cases demonstrate good agreement with experiment, good parallel speedup, and numerical convergence of the simulation results.

Fahrenthold, Eric P.↗

Coupling Finite Element and Integral Equation Representations To Efficiently Model Large Three-Dimensional Scattering Objects

The usefulness of finite element modeling follows from the ability to accurately simulate the gemetry and three-dimensional fields on the scale of a fraction of a wavelength. To make this meodelling practical for engineering design,it is necessary to integrate the stages of geometry modeling and mesh generation, numbeical solution and display of fields.

Integral Equations Finite Element↗

Notch Sensitivity of Woven Ceramic Matrix Composites Under Tensile Loading: An Experimental, Analytical, and Finite Element Study

The stress concentrations associated with circular notches and subjected to uniform tensile loading in woven ceramic matrix composites (CMCs) have been investigated for high-efficient turbine engine applications. The CMC's were composed of Nicalon silicon carbide woven fabric in SiNC matrix manufactured through polymer impregnation process (PIP). Several combinations of hole diameter/plate width ratios and ply orientations were considered in this study. In the first part, the stress concentrations were calculated measuring strain distributions surrounding the hole using strain gages at different locations of the specimens during the initial portion of the stress-strain curve before any microdamage developed. The stress concentration was also calculated analytically using Lekhnitskii's solution for orthotropic plates. A finite-width correction factor for anisotropic and orthotropic composite plate was considered. The stress distributions surrounding the circular hole of a CMC's plate were further studied using finite element analysis. Both solid and shell elements were considered. The experimental results were compared with both the analytical and finite element solutions. Extensive optical and scanning electron microscopic examinations were carried out for identifying the fracture behavior and failure mechanisms of both the notched and notched specimens. The stress concentration factors (SCF) determined by analytical method overpredicted the experimental results. But the numerical solution underpredicted the experimental SCF. Stress concentration factors are shown to increase with enlarged hole size and the effects of ply orientations on stress concentration factors are observed to be negligible. In all the cases, the crack initiated at the notch edge and propagated along the width towards the edge of the specimens.

Haque, A.↗

Three-dimensional elastic-plastic finite-element analysis of fatigue crack propagation

Fatigue cracks are a major problem in designing structures subjected to cyclic loading. Cracks frequently occur in structures such as aircraft and spacecraft. The inspection intervals of many aircraft structures are based on crack-propagation lives. Therefore, improved prediction of propagation lives under flight-load conditions (variable-amplitude loading) are needed to provide more realistic design criteria for these structures. The main thrust was to develop a three-dimensional, nonlinear, elastic-plastic, finite element program capable of extending a crack and changing boundary conditions for the model under consideration. The finite-element model is composed of 8-noded (linear-strain) isoparametric elements. In the analysis, the material is assumed to be elastic-perfectly plastic. The cycle stress-strain curve for the material is shown Zienkiewicz's initial-stress method, von Mises's yield criterion, and Drucker's normality condition under small-strain assumptions are used to account for plasticity. The three-dimensional analysis is capable of extending the crack and changing boundary conditions under cyclic loading.

Goglia, G. L.↗

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗

Cyclic creep analysis from elastic finite-element solutions

A uniaxial approach was developed for calculating cyclic creep and stress relaxation at the critical location of a structure subjected to cyclic thermomechanical loading. This approach was incorporated into a simplified analytical procedure for predicting the stress-strain history at a crack initiation site for life prediction purposes. An elastic finite-element solution for the problem was used as input for the simplified procedure. The creep analysis includes a self-adaptive time incrementing scheme. Cumulative creep is the sum of the initial creep, the recovery from the stress relaxation and the incremental creep. The simplified analysis was exercised for four cases involving a benchmark notched plate problem. Comparisons were made with elastic-plastic-creep solutions for these cases using the MARC nonlinear finite-element computer code.

Kaufman, A.↗

A finite-element perturbation approach to fluid/rotor interaction in turbomachinery elements. I - Theory. II - Application

The vibrational characteristics of a rotor that is in contact with a fluid in an annular clearance gap are investigated. The disturbance under consideration involves the axis of rotation, and includes a virtual lateral eccentricity, together with a whirling motion around the housing centerline. The fluid reaction components arise from infinitesimally small deformations with varied magnitudes which are experienced by an assembly of finite elements in the rotor-to-housing gap. A perturbation model is presented in which the perturbation equations emerge from the flow-governing equations in their discrete finite-element form. It is concluded that restrictions on the rotor-to-housing gap geometry, or the manner in which the rotor virtual eccentricity occurs are practically nonexisting. This model is used to compute the rotordynamic coefficients of an annular seal. The rotordynamic behavior of a hydraulic seal with a clearance gap depth/length ratio of 0.01 is analyzed under a cylindrical type of rotor whirl and several running speeds.

Baskharone, E. A.↗

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms

A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.

Prosser, W. H.↗