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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1973

Analysis of axisymmetrical vibration of a partially liquid-filled elastic sphere by the method of Green's function

The longitudinal dynamic behavior of launch vehicles is largely determined by their huge liquid propellant masses which are spring supported by the elastic tank walls. Throughout much of the powered flight time, these masses constitute a high percentage of the entire vehicle mass and, therefore, may dominate the fundamental modes of the vehicle. In this report, a spherical container is considered. The analysis is based on a Galerkin approach, in the course of which a second-order differential equation must be solved. The solution has been obtained by the method of Green's function. This method is favorable because it displays the manner in which the analysis can be extended to partially liquid-filled general shells of revolution. The computer programs currently available for partially liquid-filled propellant tanks are based on the finite element methods and result in analytical models having as many as several hundred degrees of freedom. The method applied in this report results in a model having less than 10 degrees of freedom as can be shown by numerical evaluation. Therefore, it will be possible to analyze propellant tanks using much less computer time with comparable accuracy.

Glaser, R. F.↗

NASA NTRS1981

Analysis of cracks emanating from a circular hole in unidirectional fiber reinforced composites, part 2

An analytical method is developed for cracks emanating from a circular hole in an off-axis unidirectional fiber-reinforced composite. The method which is formulated by using conservation laws of elasticity and fundamental relationships in anisotropic fracture mechanics, provides a convenient and accurate means to examine the complicated crack behavior, when used in conjunction with a suitable numerical scheme such as the finite element method. The formulation is eventually reduced to a system of linear algebraic equations of mixed-mode stress intensity factors. Fracture parameters, describing crack-tip deformation and fracture in the composite, are obtained explicitly. Effects of material anisotropy and crack/hole geometry are examined also. Of particular interest are the energy release rates associated with crack extension; their values are evaluated for various cases. Results show that mixed-mode stress intensity factors and energy release rates associated with the cracks emanating from a hole change very appreciably with fiber orientation in the composite. K sub 1 and G increase monotonically with increasing theta; but K sub 2 reaches its maximum at theta = 45 deg, and then decreases gradually as theta increases further.

Wang, S. S.↗

NASA NTRS2008

Damage Diagnosis in Semiconductive Materials Using Electrical Impedance Measurements

Recent aerospace industry trends have resulted in an increased demand for real-time, effective techniques for in-flight structural health monitoring. A promising technique for damage diagnosis uses electrical impedance measurements of semiconductive materials. By applying a small electrical current into a material specimen and measuring the corresponding voltages at various locations on the specimen, changes in the electrical characteristics due to the presence of damage can be assessed. An artificial neural network uses these changes in electrical properties to provide an inverse solution that estimates the location and magnitude of the damage. The advantage of the electrical impedance method over other damage diagnosis techniques is that it uses the material as the sensor. Simple voltage measurements can be used instead of discrete sensors, resulting in a reduction in weight and system complexity. This research effort extends previous work by employing finite element method models to improve accuracy of complex models with anisotropic conductivities and by enhancing the computational efficiency of the inverse techniques. The paper demonstrates a proof of concept of a damage diagnosis approach using electrical impedance methods and a neural network as an effective tool for in-flight diagnosis of structural damage to aircraft components.

Ross, Richard W.↗

NASA NTRS2014

Design, Analysis, and On-Sun Evaluation of Reflective Strips Under Controlled Buckling

Solar concentrators are envisioned for use in a variety of space-based applications, including applications involving in situ resource utilization. Identifying solar concentrators that minimize mass and cost are of great interest, especially since launch cost is driven in part by the mass of the payload. Concentrators must also be able to survive the wide temperature excursions on the lunar surface. Identifying smart structures which compensate for changes in concentrator geometry brought about by temperature extremes are of interest. Some applications may benefit from the ability to change the concentrators focal pattern at will. This paper addresses a method of designing a single reflective strip to produce a desired focal pattern through the use of controlled buckling. Small variations in the cross section over the length of the reflective strip influence the distribution of light in the focal region. A finite element method of analysis is utilized here which calculates the curve produced for a given strip cross section and axial load. Varying axial force and strip cross section over the length of the reflective strip provide a means of optimizing ray convergence in the focal region. Careful selection of a tapered cross section yields a reflective strip that approximates a parabola. An array of reflective strips under controlled buckling produces a light weight concentrator and adjustments in the compression of individual strips provide a means of compensating for temperature excursions or changing the focal pattern at will.

solar collectors↗

NASA NTRS2014

Design, Analysis, and On-Sun Evaluation of Reflective Strips Under Controlled Buckling

Solar concentrators are envisioned for use in a variety of space-based applications, including applications involving in situ resource utilization. Identifying solar concentrators that minimize mass and cost are of great interest, especially since launch cost is driven in part by the mass of the payload. Concentrators must also be able to survive the wide temperature excursions on the lunar surface. Identifying smart structures which compensate for changes in concentrator geometry brought about by temperature extremes are of interest. Some applications may benefit from the ability to change the concentrator's focal pattern at will. This paper addresses a method of designing a single reflective strip to produce a desired focal pattern through the use of controlled buckling. Small variations in the cross section over the length of the reflective strip influence the distribution of light in the focal region. A finite element method of analysis is utilized here which calculates the curve produced for a given strip cross section and axial load. Varying axial force and strip cross section over the length of the reflective strip provide a means of optimizing ray convergence in the focal region. Careful selection of a tapered cross section yields a reflective strip that approximates a parabola. An array of reflective strips under controlled buckling produces a light weight concentrator and adjustments in the compression of individual strips provide a means of compensating for temperature excursions or changing the focal pattern at will.

buckling↗

NASA NTRS2017

Near Real-Time Probabilistic Damage Diagnosis Using Surrogate Modeling and High Performance Computing

This work investigates novel approaches to probabilistic damage diagnosis that utilize surrogate modeling and high performance computing (HPC) to achieve substantial computational speedup. Motivated by Digital Twin, a structural health management (SHM) paradigm that integrates vehicle-specific characteristics with continual in-situ damage diagnosis and prognosis, the methods studied herein yield near real-time damage assessments that could enable monitoring of a vehicle's health while it is operating (i.e. online SHM). High-fidelity modeling and uncertainty quantification (UQ), both critical to Digital Twin, are incorporated using finite element method simulations and Bayesian inference, respectively. The crux of the proposed Bayesian diagnosis methods, however, is the reformulation of the numerical sampling algorithms (e.g. Markov chain Monte Carlo) used to generate the resulting probabilistic damage estimates. To this end, three distinct methods are demonstrated for rapid sampling that utilize surrogate modeling and exploit various degrees of parallelism for leveraging HPC. The accuracy and computational efficiency of the methods are compared on the problem of strain-based crack identification in thin plates. While each approach has inherent problem-specific strengths and weaknesses, all approaches are shown to provide accurate probabilistic damage diagnoses and several orders of magnitude computational speedup relative to a baseline Bayesian diagnosis implementation.

Warner, James E.↗

NASA NTRS1984

Newton like: Minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗

NASA NTRS1985

Newton-like minimal residual methods applied to transonic flow calculations

Wong, Y. S.↗

NASA NTRS1987

Mechanical characteristics of filament-wound pressure vessel (burst pressure)

The finite element method is used to analyze the mechanical characteristics of a pressurized filament-wound (FW) pressure vessel, and to predict its burst pressure. The analysis takes into account the bending moment, the stretch-bend coupling effect, nonlinear stress-strain relations, and finite deflection. The analysis is based on two initial failure criteria for laminae, and two ultimate fracture criteria for laminated structures. The numerical results, obtained by applying the load incremental method to the isotensoid CFRP pressure vessel used in the launching of the Zikiken satellite, are in good agreement with the experimental burst pressure and fracture behaviors.

Iida, H.↗

NASA NTRS2016

An Immersed Boundary Method for Solving the Compressible Navier-Stokes Equations with Fluid Structure Interaction

An immersed boundary method for the compressible Navier-Stokes equation and the additional infrastructure that is needed to solve moving boundary problems and fully coupled fluid-structure interaction is described. All the methods described in this paper were implemented in NASA's LAVA solver framework. The underlying immersed boundary method is based on the locally stabilized immersed boundary method that was previously introduced by the authors. In the present paper this method is extended to account for all aspects that are involved for fluid structure interaction simulations, such as fast geometry queries and stencil computations, the treatment of freshly cleared cells, and the coupling of the computational fluid dynamics solver with a linear structural finite element method. The current approach is validated for moving boundary problems with prescribed body motion and fully coupled fluid structure interaction problems in 2D and 3D. As part of the validation procedure, results from the second AIAA aeroelastic prediction workshop are also presented. The current paper is regarded as a proof of concept study, while more advanced methods for fluid structure interaction are currently being investigated, such as geometric and material nonlinearities, and advanced coupling approaches.

Equations↗

NASA NTRS1994

Buckling and vibration analysis of laminated panels using VICONOPT

The analysis aspects of the 23,000-line FORTRAN program VIPASA with constraints and optimization (VICONOPT) are described. Overall stiffness matrices assembled from the earlier exact vibration and instability of plate assemblies including shear and anisotropy (VIPASA) flat plate stiffnesses are optionally coupled by Langrangian multipliers to find critical buckling loads, or natural frequencies of undamped vibration, of prismatic assemblies of anisotropic flat plates with arbitrarily located point supports or simple transverse supporting frames. The longitudunal continuity of typical wing and fuselage panels is closely approximated because the solutions are for the infinitely long structure obtained by repeating a bay and its supports longitudinally. Any longitudinally invariant in-plane plate stresses are permitted, and very rapid solutions are guaranteed by numerous refinements, including multilevel substructuring and a method for repetitive cross sections that is exact for regular polygons used to represent cylinders. Modal displacements and stresses in or between plies of laminated plates are calculated and plotted, with values being recovered at all nodes of substructures. Comparison with usual approximate finite-element methods confirms that, for comparably converged solutions, VICONOPT is typically between 100 and 10(exp 4) times faster.

Kennedy, David↗

NASA NTRS2013

ONERA-NASA Cooperative Effort on Liner Impedance Eduction

As part of a cooperation between ONERA and NASA, the liner impedance eduction methods developed by the two research centers are compared. The NASA technique relies on an objective function built on acoustic pressure measurements located on the wall opposite the test liner, and the propagation code solves the convected Helmholtz equation in uniform ow using a finite element method that implements a continuous Galerkin discretization. The ONERA method uses an objective function based either on wall acoustic pressure or on acoustic velocity acquired above the liner by Laser Doppler Anemometry, and the propagation code solves the linearized Euler equations by a discontinuous Galerkin discretization. Two acoustic liners are tested in both ONERA and NASA ow ducts and the measured data are treated with the corresponding impedance eduction method. The first liner is a wire mesh facesheet mounted onto a honeycomb core, designed to be linear with respect to incident sound pressure level and to grazing ow velocity. The second one is a conventional, nonlinear, perforate-over-honeycomb single layer liner. Configurations without and with ow are considered. For the nonlinear liner, the comparison of liner impedance educed by NASA and ONERA shows a sensitivity to the experimental conditions, namely to the nature of the source and to the sample width.

Primus, Julien↗

NASA NTRS1976

Disjoint design spaces in the optimization of harmonically excited structures

The paper deals with a method of optimizing harmonically or dynamically excited structures which removes the frequency constraint (that specifies that the structure's first natural frequency of vibration must be greater than the excitation frequency), replaces the displacement constraint by an inequality constraint on the allowable stress and applies finite element approximations to the continuous one-dimensional structures. When the problem is posed in this manner, the feasible design space is disconnected or 'disjoint'. An attempt is made to explain this disjointness by studying the optimal design of a thin-walled cantilevered rod subjected to steady, simple harmonic torsional excitation. The finite element method is applied to the optimization of a harmonically excited structure with an arbitrary number of design variables, with the disjointness of the design space, a complicating factor.

Johnson, E. H.↗

NASA NTRS2012

Nonlinear Analysis of Bonded Composite Tubular Lap Joints

The present study describes a semi-analytical solution method for predicting the geometrically nonlinear response of a bonded composite tubular single-lap joint subjected to general loading conditions. The transverse shear and normal stresses in the adhesive as well as membrane stress resultants and bending moments in the adherends are determined using this method. The method utilizes the principle of virtual work in conjunction with nonlinear thin-shell theory to model the adherends and a cylindrical shear lag model to represent the kinematics of the thin adhesive layer between the adherends. The kinematic boundary conditions are imposed by employing the Lagrange multiplier method. In the solution procedure, the displacement components for the tubular joint are approximated in terms of non-periodic and periodic B-Spline functions in the longitudinal and circumferential directions, respectively. The approach presented herein represents a rapid-solution alternative to the finite element method. The solution method was validated by comparison against a previously considered tubular single-lap joint. The steep variation of both peeling and shearing stresses near the adhesive edges was successfully captured. The applicability of the present method was also demonstrated by considering tubular bonded lap-joints subjected to pure bending and torsion.

E. Oterkus↗

NASA NTRS1987

Impact of solids on composite materials

The failure modes of composite materials as a result of low velocity impact were investigated by simulating the impact with a finite element analysis. An important facet of the project is the modeling of the impact of a solid onto cylindrical shells composed of composite materials. The model under development will simulate the delamination sustained when a composite material encounters impact from another rigid body. The computer equipment was installed, the computer network tested, and a finite element method model was developed to compare results with known experimental data. The model simulated the impact of a steel rod onto a rotating shaft. Pre-processing programs (GMESH and TANVEL) were developed to generate node and element data for the input into the three dimensional, dynamic finite element analysis code (DYNA3D). The finite element mesh was configured with a fine mesh near the impact zone and a coarser mesh for the impacting rod and the regions surrounding the impacting zone. For the computer simulation, five impacting loads were used to determine the time history of the stresses, the scribed surface areas, and the amount of ridging. The processing time of the computer codes amounted from 1 to 4 days. The calculated surface area were within 6-12 percent, relative error when compated to the actual scratch area.

Bronson, Arturo↗

NASA NTRS1990

The effect of coatings and liners on heat transfer in a dry shaft-bush tribosystem

The temperatures due to frictional heating within a solid lubricated or coated journal bearing were analyzed by using a finite element method. A solid model of the shaft-bush tribocontact was generated with an eight-node, three-dimensional, first-order isoparametric heat-transfer element and the Patran solid modeler software. The Patmar (Patran-Marc) translator was used to help develop the Marc-based finite element program for the system; this software was used on the Cray X-MP supercomputer to perform a finite element analysis of the contact. The analysis was performed for various liner materials, for thin, hard, wear-resistant coated bearings, and for different geometries and thermal cooling boundary conditions. The analyses indicated that thermal conductivity of the liner or coating material is the most vital thermal parameter that controls the interface temperature. In addition to design variations, the proximity of the cooling source to the heat-flux-generating interface is critically important to the temperature control in the system.

Ghosh, Mihir K.↗

NASA NTRS2000

An Integrated NDE and FEM Characterization of Composite Rotors

A structural assessment by integrating finite-element methods (FEM) and a nondestructive evaluation (NDE) of two flywheel rotor assemblies is presented. Composite rotor A is pancake like with a solid hub design, and composite rotor B is cylindrical with a hollow hub design. Detailed analyses under combined centrifugal and interference-fit loading are performed. Two- and three-dimensional stress analyses and two-dimensional fracture mechanics analyses are conducted. A comparison of the structural analysis results obtained with those extracted via NDE findings is reported. Contact effects due to press-fit conditions are evaluated. Stress results generated from the finite-element analyses were corroborated with the analytical solution. Cracks due to rotational loading up to 49 000 rpm for rotor A and 34 000 rpm for rotor B were successfully imaged with NDE and predicted with FEM and fracture mechanics analyses. A procedure that extends current structural analysis to a life prediction tool is also defined.

Abdul-Aziz, Ali↗

NASA NTRS1992

Method for improved prediction of bone fracture risk using bone mineral density in structural analysis

A non-invasive in-vivo method of analyzing a bone for fracture risk includes obtaining data from the bone such as by computed tomography or projection imaging which data represents a measure of bone material characteristics such as bone mineral density. The distribution of the bone material characteristics is used to generate a finite element method (FEM) mesh from which load capability of the bone can be determined. In determining load capability, the bone is mathematically compressed, and stress, strain force, force/area versus bone material characteristics are determined.

Cann, Christopher E.↗