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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Finite element methods in fracture mechanics

Finite-element methodology specific to the analysis of fracture mechanics problems is reviewed. Primary emphasis is on the important algorithmic developments which have enhanced the numerical modeling of fracture processes. Methodologies to address elastostatic problems in two and three dimensions, elastodynamic problems, elastoplastic problems, special considerations for three-dimensional nonlinear problems, and the modeling of stable crack growth are reviewed. In addition, the future needs of the fracture community are discussed and open questions are identified.

Liebowitz, H.↗
NASA NTRS1989

Vector algorithms for geometrically nonlinear 3D finite element analysis

Algorithms for geometrically nonlinear finite element analysis are presented which exploit the vector processing capability of the VPS-32, which is closely related to the CYBER 205. By manipulating vectors (which are long lists of numbers) rather than individual numbers, very high processing speeds are obtained. Long vector lengths are obtained without extensive replication or reordering by storage of intermediate results in strategic patterns at all stages of the computations. Comparisons of execution times with those from programs using either scalar or other vector programming techniques indicate that the algorithms presented are quite efficient.

Whitcomb, John D.↗
NASA NTRS1991

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

Iannelli, G. S.↗
NASA NTRS1991

Fatigue crack growth reliability by probabilistic finite elements

Fusion of the probabilistic finite-element method and reliability analysis for probabilistic fatigue-crack growth is presented. A comprehensive method for determining the probability of fatigue failure for mixed-mode cyclic loading is also presented. The loading is mixed-mode with randomness in the initial and final crack lengths, initial crack angle and position, material properties, crack-growth law, crack-direction law, and loading. The methodology consists of calculating the reliability index via an optimization procedure which is used to calculate the probability of fatigue failure. Performance of the methodology presented is demonstrated on a classical mode-I fatigue problem.

Besterfield, Glen H.↗
NASA NTRS1993

Large-scale computation of incompressible viscous flow by least-squares finite element method

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to large-scale/three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations and results in symmetric, positive definite algebraic system which can be solved effectively by simple iterative methods. The first-order velocity-Bernoulli function-vorticity formulation for incompressible viscous flows is also tested. For three-dimensional cases, an additional compatibility equation, i.e., the divergence of the vorticity vector should be zero, is included to make the first-order system elliptic. The simple substitution of the Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. To show the validity of this scheme for large-scale computation, we give numerical results for 2D driven cavity problem at Re = 10000 with 408 x 400 bilinear elements. The flow in a 3D cavity is calculated at Re = 100, 400, and 1,000 with 50 x 50 x 50 trilinear elements. The Taylor-Goertler-like vortices are observed for Re = 1,000.

Jiang, Bo-Nan↗
NASA NTRS2007

Analysis of Fluid Gauge Sensor for Zero or Microgravity Conditions using Finite Element Method

In this paper the Finite Element Method (FEM) is presented for mass/volume gauging of a fluid in a tank subjected to zero or microgravity conditions. In this approach first mutual capacitances between electrodes embedded inside the tank are measured. Assuming the medium properties the mutual capacitances are also estimated using FEM approach. Using proper non-linear optimization the assumed properties are updated by minimizing the mean square error between estimated and measured capacitances values. Numerical results are presented to validate the present approach.

Deshpande, Manohar D.↗
NASA NTRS2007

An Approach for Assessing Delamination Propagation Capabilities in Commercial Finite Element Codes

An approach for assessing the delamination propagation capabilities in commercial finite element codes is presented and demonstrated for one code. For this investigation, the Double Cantilever Beam (DCB) specimen and the Single Leg Bending (SLB) specimen were chosen for full three-dimensional finite element simulations. First, benchmark results were created for both specimens. Second, starting from an initially straight front, the delamination was allowed to propagate. Good agreement between the load-displacement relationship obtained from the propagation analysis results and the benchmark results could be achieved by selecting the appropriate input parameters. Selecting the appropriate input parameters, however, was not straightforward and often required an iterative procedure. Qualitatively, the delamination front computed for the DCB specimen did not take the shape of a curved front as expected. However, the analysis of the SLB specimen yielded a curved front as may be expected from the distribution of the energy release rate and the failure index across the width of the specimen. Overall, the results are encouraging but further assessment on a structural level is required.

Krueger, Ronald↗
NASA NTRS2007

An Approach for Assessing Delamination Propagation Capabilities in Commercial Finite Element Codes

An approach to assessing the delamination propagation capabilities in commercial finite element codes is presented and demonstrated for one code. For this investigation, the Double Cantilever Beam (DCB) specimen and the Single Leg Bending (SLB) specimen were chosen for full three-dimensional finite element simulations. First, benchmark results were created for both specimens. Second, starting from an initially straight front, the delamination was allowed to propagate. Good agreement between the load-displacement relationship obtained from the propagation analysis results and the benchmark results could be achieved by selecting the appropriate input parameters. Selecting the appropriate input parameters, however, was not straightforward and often required an iterative procedure. Qualitatively, the delamination front computed for the DCB specimen did not take the shape of a curved front as expected. However, the analysis of the SLB specimen yielded a curved front as may be expected from the distribution of the energy release rate and the failure index across the width of the specimen. Overall, the results are encouraging but further assessment on a structural level is required.

Krueger, Ronald↗
NASA NTRS1994

Dynamical observer for a flexible beam via finite element approximations

The purpose of this view-graph presentation is a computational investigation of the closed-loop output feedback control of a Euler-Bernoulli beam based on finite element approximation. The observer is part of the classical observer plus state feedback control, but it is finite-dimensional. In the theoretical work on the subject it is assumed (and sometimes proved) that increasing the number of finite elements will improve accuracy of the control. In applications, this may be difficult to achieve because of numerical problems. The main difficulty in computing the observer and simulating its work is the presence of high frequency eigenvalues in the finite-element model and poor numerical conditioning of some of the system matrices (e.g. poor observability properties) when the dimension of the approximating system increases. This work dealt with some of these difficulties.

Manitius, Andre↗
NASA NTRS1992

Coupling finite element and integral equation solutions using decoupled boundary meshes

A method is outlined for calculating scattered fields from inhomogeneous penetrable objects using a coupled finite element-integral equation solution. The finite element equation can efficiently model fields in penetrable and inhomogeneous regions, while the integral equation exactly models fields on the finite element mesh boundary and in the exterior region. By decoupling the interior finite element and exterior integral equation meshes, considerable flexibility is found in both the number of field expansion points as well as their density. Only the nonmetal portions of the object need be modeled using a finite element expansion; exterior perfect conducting surfaces are modeled using an integral equation with a single unknown field since E(tan) is identically zero on these surfaces. Numerical convergence, accuracy, and stability at interior resonant frequencies are studied in detail.

Cwik, Tom↗
NASA NTRS2007

Finite Element Analysis of the LOLA Receiver Telescope Lens

This paper presents the finite element stress and distortion analysis completed on the Receiver Telescope lens of the Lunar Orbiter Laser Altimeter (LOLA). LOLA is one of six instruments on the Lunar Reconnaissance Orbiter (LRO), scheduled to launch in 2008. LOLA's main objective is to produce a high-resolution global lunar topographic model to aid in safe landings and enhance surface mobility in future exploration missions. The Receiver Telescope captures the laser pulses transmitted through a diffractive optical element (DOE) and reflected off the lunar surface. The largest lens of the Receiver Telescope, Lens 1, is a 150 mm diameter aspheric lens originally designed to be made of BK7 glass. The finite element model of the Receiver Telescope Lens 1 is comprised of solid elements and constrained in a manner consistent with the behavior of the mounting configuration of the Receiver Telescope tube. Twenty-one temperature load cases were mapped to the nodes based on thermal analysis completed by LOLA's lead thermal analyst, and loads were applied to simulate the preload applied from the ring flexure. The thermal environment of the baseline design (uncoated BK7 lens with no baffle) produces large radial and axial gradients in the lens. These large gradients create internal stresses that may lead to part failure, as well as significant bending that degrades optical performance. The high stresses and large distortions shown in the analysis precipitated a design change from BK7 glass to sapphire.

Matzinger, Elizabeth↗
NASA NTRS2006

SUPG Finite Element Simulations of Compressible Flows

The Streamline-Upwind Petrov-Galerkin (SUPG) finite element simulations of compressible flows is presented. The topics include: 1) Introduction; 2) SUPG Galerkin Finite Element Methods; 3) Applications; and 4) Bibliography.

Kirk, Brnjamin, S.↗
NASA NTRS1997

Error Estimation and h-Adaptivity for Optimal Finite Element Analysis

The objective of adaptive meshing and automatic error control in finite element analysis is to eliminate the need for the application engineer from re-meshing and re-running design simulations to verify numerical accuracy. The user should only need to enter the component geometry and a coarse finite element mesh. The software will then autonomously and adaptively refine this mesh where needed, reducing the error in the fields to a user prescribed value. The ideal end result of the simulation is a measurable quantity (e.g. scattered field, input impedance), calculated to a prescribed error, in less time and less machine memory than if the user applied typical uniform mesh refinement by hand. It would also allow for the simulation of larger objects since an optimal mesh is created.

Cwik, Tom↗
NASA NTRS2011

Energy Finite Element Analysis Developments for Vibration Analysis of Composite Aircraft Structures

The Energy Finite Element Analysis (EFEA) has been utilized successfully for modeling complex structural-acoustic systems with isotropic structural material properties. In this paper, a formulation for modeling structures made out of composite materials is presented. An approach based on spectral finite element analysis is utilized first for developing the equivalent material properties for the composite material. These equivalent properties are employed in the EFEA governing differential equations for representing the composite materials and deriving the element level matrices. The power transmission characteristics at connections between members made out of non-isotropic composite material are considered for deriving suitable power transmission coefficients at junctions of interconnected members. These coefficients are utilized for computing the joint matrix that is needed to assemble the global system of EFEA equations. The global system of EFEA equations is solved numerically and the vibration levels within the entire system can be computed. The new EFEA formulation for modeling composite laminate structures is validated through comparison to test data collected from a representative composite aircraft fuselage that is made out of a composite outer shell and composite frames and stiffeners. NASA Langley constructed the composite cylinder and conducted the test measurements utilized in this work.

Vlahopoulos, Nickolas↗
NASA NTRS1991

A finite element approach for the dynamic analysis of joint-dominated structures

A finite element method to model dynamic structural systems undergoing large rotations is presented. The dynamic systems are composed of rigid joint bodies and flexible beam elements. The configurations of these systems are subject to change due to the relative motion in the joints among interconnected elastic beams. A body fixed reference is defined for each joint body to describe the joint body's displacements. Using the finite element method and the kinematic relations between each flexible element and its corotational reference, the total displacement field of an element, which contains gross rigid as well as elastic effects, can be derived in terms of the translational and rotational displacements of the two end nodes. If one end of an element is hinged to a joint body, the joint body's displacements and the hinge degree of freedom at the end are used to represent the nodal displacements. This results in a highly coupled system of differential equations written in terms of hinge degrees of freedom as well as the rotational and translational displacements of joint bodies and element nodes.

Chang, Che-Wei↗
NASA NTRS2012

Development and Calibration of a System-Integrated Rotorcraft Finite Element Model for Impact Scenarios

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber (DEA) under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. The presence of this energy absorbing device reduced the peak impact acceleration levels by a factor of three. Accelerations and kinematic data collected from the crash tests were compared to a system-integrated finite element model of the test article developed in parallel with the test program. In preparation for the full-scale crash test, a series of sub-scale and MD-500 mass simulator tests were conducted to evaluate the impact performances of various components and subsystems, including new crush tubes and the DEA blocks. Parameters defined for the system-integrated finite element model were determined from these tests. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the full-scale crash test without the DEA. This combination of heuristic and quantitative methods identified modeling deficiencies, evaluated parameter importance, and proposed required model changes. The multidimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and copilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. One lesson learned was that this approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and pretest predictions. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, potentially reducing overall development costs.

Annett, Martin S.↗
NASA NTRS1972

Progress in nonlinear finite element analysis using asymptotic solution techniques.

The existence of standard forms of finite element equations for nonlinear analysis is observed. These are noted to provide a sufficiently special circumstance to admit the construction of asymptotic solution techniques. Although efforts in this direction have been made only recently, a number of distinct approaches to asymptotic solution techniques for finite element analysis have been put forward and are examined herein. Consideration is given to the differences and similarities of the approaches and the resulting methods of solution. Special characteristics regarding applicability, accuracy and efficiency are delineated.

Mallett, R. H.↗
NASA NTRS1985

Finite element thermal-structural analyses of a cable-stiffened orbiting antenna

Finite element thermal-structural analyses of a cable-stiffened orbiting antenna are presented. The determination of prestresses in the antenna is described first. Heating and thermal analyses for orbiting space structures are then discussed briefly. Structural deformations and stresses are presented for three finite element structural analysis approaches: (1) small deflections, (2) stress-stiffening, and (3) large deflections. The accuracy of the three analysis approaches is evaluated for the orbiting antenna at different prestress levels.

Thornton, E. A.↗