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 1,441 records · Page 80

Finite element methods for integrated aerodynamic heating analysis

This report gives a description of the work which has been undertaken during the second year of a three year research program. The objectives of the program are to produce finite element based procedures for the solution of the large scale practical problems which are of interest to the Aerothermal Loads Branch (ALB) at NASA Langley Research Establishment. The problems of interest range from Euler simulations of full three dimensional vehicle configurations to local analyses of three dimensional viscous laminar flow. Adaptive meshes produced for both steady state and transient problems are to be considered. An important feature of the work is the provision of specialized techniques which can be used at ALB for the development of an integrated fluid/thermal/structural modeling capability.

Morgan, K.↗

Fluid-structural interactions using Navier-Stokes flow equations coupled with shell finite element structures

A computational procedure is presented to study fluid-structural interaction problems for three-dimensional aerospace structures. The flow is modeled using the three-dimensional unsteady Euler/Navier-Stokes equations and solved using the finite-difference approach. The three dimensional structure is modeled using shell/plate finite-element formulation. The two disciplines are coupled using a domain decomposition approach. Accurate procedures both in time and space are developed to combine the solutions from the flow equations with those of the structural equations. Time accuracy is maintained using aeroelastic configuration-adaptive moving grids that are computed every time step. The work done by aerodynamic forces due to structural deformations is preserved using consistent loads. The present procedure is validated by computing the aeroelastic response of a wing and comparing with experiment. Results are illustrated for a typical wing-body configuration.

Guruswamy, Guru P.↗

MGGHAT: Elliptic PDE software with adaptive refinement, multigrid and high order finite elements

MGGHAT (MultiGrid Galerkin Hierarchical Adaptive Triangles) is a program for the solution of linear second order elliptic partial differential equations in two dimensional polygonal domains. This program is now available for public use. It is a finite element method with linear, quadratic or cubic elements over triangles. The adaptive refinement via newest vertex bisection and the multigrid iteration are both based on a hierarchical basis formulation. Visualization is available at run time through an X Window display, and a posteriori through output files that can be used as GNUPLOT input. In this paper, we describe the methods used by MGGHAT, define the problem domain for which it is appropriate, illustrate use of the program, show numerical and graphical examples, and explain how to obtain the software.

Mitchell, William F.↗

Finite Element Analysis of Magnetic Damping Effects on G-Jitter Induced Fluid Flow

This paper reports some interim results on numerical modeling and analyses of magnetic damping of g-jitter driven fluid flow in microgravity. A finite element model is developed to represent the fluid flow, thermal and solute transport phenomena in a 2-D cavity under g-jitter conditions with and without an applied magnetic field. The numerical model is checked by comparing with analytical solutions obtained for a simple parallel plate channel flow driven by g-jitter in a transverse magnetic field. The model is then applied to study the effect of steady state g-jitter induced oscillation and on the solute redistribution in the liquid that bears direct relevance to the Bridgman-Stockbarger single crystal growth processes. A selection of computed results is presented and the results indicate that an applied magnetic field can effectively damp the velocity caused by g-jitter and help to reduce the time variation of solute redistribution.

Pan, Bo↗

Integration of MAC/GMC into CalculiX, an Open Source Finite Element Code

An analysis framework is presented that makes available multiscale analysis of composite structures using the open-source FEA solver package CalculiX CrunchiX (CCX). At the center of this framework is the coupling of the Finite Element Analysis - Micromechanics Analysis Code (FEAMAC) library from NASA's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) to allow micromechanics analysis. The results show that the proposed coupling can be used with appropriate care for multiscale FEA simulations of composite materials. The largest error reported in this validation was in a four-point bend test specimen with an error of less than 1% difference in the maximum deflection of the beam.

Yapor Genao, Francisco A.↗

Finite element modeling of stress in the Nazca plate - Driving forces and plate boundary earthquakes

The state of stress within the Nazca plate due to plate driving forces and large plate boundary earthquakes has been analyzed by applying a finite element method using the wave front solution technique to models of the intraplate stress field in a single plate using a refined grid. Although only static elastic models have been explicitly calculated, certain limiting cases of an elastic plate over a viscous asthenosphere were also treated. A state of nearly east-west compression inferred from the source mechanism of thrust earthquakes in the interior of the plate requires ridge pushing forces. The net pulling force on the oceanic plate by the subducted slab has a maximum value comparable to pushing forces. The estimated horizontal deviatoric stress in intraplate regions, based on potential forces associated with the ridge, is on the order of a few hundred bars. The intraplate stress field in the region of the 1960 earthquake may change by a few tens of bars at most once the asthenosphere has relaxed, with changes on the order of one bar occurring at greater distances into the plate. The changes in the intraplate stress field are probably not noticeable unless the lithosphere is near failure.

Richardson, R. M.↗

Finite element analysis of electromagnetic propagation in an absorbing wave guide

Wave guides play a significant role in microwave space communication systems. The attenuation per unit length of the guide depends on its construction and design frequency range. A finite element Galerkin formulation has been developed to study TM electromagnetic propagation in complex two-dimensional absorbing wave guides. The analysis models the electromagnetic absorptive characteristics of a general wave guide which could be used to determine wall losses or simulate resistive terminations fitted into the ends of a guide. It is believed that the general conclusions drawn by using this simpler two-dimensional geometry will be fundamentally the same for other geometries.

Baumeister, Kenneth J.↗

An independent refinement and integration procedure in multiregion finite element analysis

An independent refinement and integration procedure is developed to couple together independently modeled (global and local) regions in a single analysis. The models can have different levels of refinement and along the interface between them the finite element nodes need not coincide with one another. In the local model all the nodes except the nodes at the interface are statically condensed and the reduced stiffness matrix is obtained. For this static condensation a modified frontal solution technique is employed. A spline interpolation function that satisfies the linear isotropic plate bending differential equation is used to relate the local model interface nodal displacements to the global model interface displacements. The proposed independent refinement and integration procedure is evaluated by applying it to two- and three-dimensional cases involving inplane and out-of-plane deformation. The procedure yielded very accurate results for all the examples studied.

Krishnamurthy, T.↗

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.↗

Transpose of finite-element data

TRANSPOSE computer program examines single point of structural-analysis model under many loads. Program helps in data reduction and analysis, saves output for subsequent postprocessing, and reduces time required for structure analysis. Essentially, program transposes finite-element data from one loading condition for all element and node data to one-node-point data for all loading conditions.

Furuike, T.↗

Finite element solutions for crack-tip behavior in small-scale yielding

The subject considered is the stress and deformation fields in a cracked elastic-plastic power law hardening material under plane strain tensile loading. An incremental plasticity finite element formulation is developed for accurate analysis of the complete field problem including the extensively deformed near tip region, the elastic-plastic region, and the remote elastic region. The formulation has general applicability and was used to solve the small scale yielding problem for a set of material hardening exponents. Distributions of stress, strain, and crack opening displacement at the crack tip and through the elastic-plastic zone are presented as a function of the elastic stress intensity factor and material properties.

Tracey, D. M.↗

Application of the finite element method to rotary-wing aeroelasticity

Recent research in rotary-wing aeroelasticity has indicated that all fundamental problems in this area are inherently nonlinear. The non-linearities in this problem are due to the inclusion of finite slopes, due to moderate deflections, in the structural, inertia and aerodynamic operators associated with this aeroelastic problem. In this paper the equations of motion, which are both time and space dependent, for the aeroelastic problem are first formulated in P.D.E. form. Next the equations are linearized about a suitable equilibrium position. The spatial dependence in these equations is discretized using a local Galerkin method of weighted residuals resulting in a finite element formulation of the aeroelastic problem. As an illustration the method is applied to the coupled flap-lag problem of a helicopter rotor blade in hover. Comparison of the solutions with previously published solutions establishes the convergence properties of the method. It is concluded that this formulation is a practical tool for solving rotary-wing aeroelastic stability or response problems.

Friedmann, P.↗

Finite element analysis of moisture effects in graphite-epoxy composites

A numerical analysis procedure is developed for the micro-mechanical investigation of moisture effects on graphite-epoxy composites. The procedure employs a nonlinear constitutive relationship in a finite element formulation. A quasi-Newton method is used to solve the resulting nonlinear algebraic equations. The results presented indicate that the stresses at the graphite-epoxy interface have a strong dependence on moisture content, fiber spacing, and applied load.

Tobing, F. E. L.↗

Validation of finite element and boundary element methods for predicting structural vibration and radiated noise

Analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise radiated from the box. The FEM was used to predict the vibration, and the surface vibration was used as input to the BEM to predict the sound intensity and sound power. Vibration predicted by the FEM model was validated by experimental modal analysis. Noise predicted by the BEM was validated by sound intensity measurements. Three types of results are presented for the total radiated sound power: (1) sound power predicted by the BEM modeling using vibration data measured on the surface of the box; (2) sound power predicted by the FEM/BEM model; and (3) sound power measured by a sound intensity scan. The sound power predicted from the BEM model using measured vibration data yields an excellent prediction of radiated noise. The sound power predicted by the combined FEM/BEM model also gives a good prediction of radiated noise except for a shift of the natural frequencies that are due to limitations in the FEM model.

Seybert, A. F.↗

Static Aeroelastic and Longitudinal Trim Model of Flexible Wing Aircraft Using Finite-Element Vortex-Lattice Coupled Solution

This paper presents a static aeroelastic model and longitudinal trim model for the analysis of a flexible wing transport aircraft. The static aeroelastic model is built using a structural model based on finite-element modeling and coupled to an aerodynamic model that uses vortex-lattice solution. An automatic geometry generation tool is used to close the loop between the structural and aerodynamic models. The aeroelastic model is extended for the development of a three degree-of-freedom longitudinal trim model for an aircraft with flexible wings. The resulting flexible aircraft longitudinal trim model is used to simultaneously compute the static aeroelastic shape for the aircraft model and the longitudinal state inputs to maintain an aircraft trim state. The framework is applied to an aircraft model based on the NASA Generic Transport Model (GTM) with wing structures allowed to flexibly deformed referred to as the Elastically Shaped Aircraft Concept (ESAC). The ESAC wing mass and stiffness properties are based on a baseline "stiff" values representative of current generation transport aircraft.

Aeroelasticity↗

A split finite element algorithm for the compressible Navier-Stokes equations

An accurate and efficient numerical solution algorithm is established for solution of the high Reynolds number limit of the Navier-Stokes equations governing the multidimensional flow of a compressible essentially inviscid fluid. Finite element interpolation theory is used within a dissipative formulation established using Galerkin criteria within the Method of Weighted Residuals. An implicit iterative solution algorithm is developed, employing tensor product bases within a fractional steps integration procedure, that significantly enhances solution economy concurrent with sharply reduced computer hardware demands. The algorithm is evaluated for resolution of steep field gradients and coarse grid accuracy using both linear and quadratic tensor product interpolation bases. Numerical solutions for linear and nonlinear, one, two and three dimensional examples confirm and extend the linearized theoretical analyses, and results are compared to competitive finite difference derived algorithms.

Baker, A. J.↗

Finite Element Analysis System Workflow Tools

A collection of MATLAB functions and class definitions called System Workflow Tools (SWFT) are available to semi-automate steps in the simulation process. Some of these steps are often simple and routine for smaller finite element models, but if done directly by an analyst can quickly become labor intensive, cumbersome, and error prone for larger, system level models. Some of SWFT’s capabilities demonstrated in this report includes writing Sierra input decks and processing Quantities of Interest (QOI) from results files. SWFT also writes scripts in order to utilize other software programs such as Cubit (separating system level CAD into subassemblies and components, creating nodesets and sidesets), DAKOTA (ensemble management), and ParaView (contour plots and animations). Detailed commands and workflows from mesh generation to report generation are provided as examples for analysts to utilize SWFT capabilities.

97 MATHEMATICS AND COMPUTING↗