SEARCH · Search NASA

Results for “Boundary finite element”

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 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS2012

TPSAS-NF1676L-13539-DND

The performance of variable stiffness composite panels is characterized using experiments and structural analyses. Fabricated using a fiber placement system, the fiber angle in these panels varies continuously within each ply. Measured panel shapes and thermal responses are correlated with finite element analyses. Large geometric imperfections require unique fixtures to both straighten the panel edges and provide structural test boundary conditions. Panels show linear prebuckling and nonlinear postbuckling under applied end shortening, which occur at higher loads than a baseline panel. Panel responses to in-plane loading are computed using geometric and materially nonlinear analyses that include mechanical prestresses and measured imperfections. Analytical strains and displacements also correlate well with measured results.

Chauncey Wu↗

NASA NTRS2023

Acoustics of a Supersonic Mach 1.4 Axisymmetric Spike Inlet

An engine inlet for a supersonic commercial aircraft is analyzed for noise radiation in low-speed configurations for which auxiliary inlet doors are open. Duct modes generated by the fan are estimated by post-processing unsteady Reynold averaged Navier-Stokes (URANS) simulations of the fan and are used as boundary conditions to the inlet domain. Far-field inlet radiated tone levels are computed using commercial acoustics finite element software. The results indicated that the inlet tone levels at cutback and sideline conditions are higher by more than 20 dB than the approach levels.

Supersonics↗

NASA NTRS1987

Thermal finite-element analysis of space shuttle main engine turbine blade

Finite-element, transient heat transfer analyses were performed for the first-stage blades of the space shuttle main engine (SSME) high-pressure fuel turbopump. The analyses were based on test engine data provided by Rocketdyne. Heat transfer coefficients were predicted by performing a boundary-layer analysis at steady-state conditions with the STAN5 boundary-layer code. Two different peak-temperature overshoots were evaluated for the startup transient. Cutoff transient conditions were also analyzed. A reduced gas temperature profile based on actual thermocouple data was also considered. Transient heat transfer analyses were conducted with the MARC finite-element computer code.

Abdul-Aziz, Ali↗

NASA NTRS1983

An investigation of the accuracy of finite difference methods in the solution of linear elasticity problems

The accuracy of the finite difference method in the solution of linear elasticity problems that involve either a stress discontinuity or a stress singularity is considered. Solutions to three elasticity problems are discussed in detail: a semi-infinite plane subjected to a uniform load over a portion of its boundary; a bimetallic plate under uniform tensile stress; and a long, midplane symmetric, fiber reinforced laminate subjected to uniform axial strain. Finite difference solutions to the three problems are compared with finite element solutions to corresponding problems. For the first problem a comparison with the exact solution is also made. The finite difference formulations for the three problems are based on second order finite difference formulas that provide for variable spacings in two perpendicular directions. Forward and backward difference formulas are used near boundaries where their use eliminates the need for fictitious grid points.

Bauld, N. R., Jr.↗

NASA NTRS1998

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

NASA NTRS1973

Finite Element Method for Thermal Analysis

A two- and three-dimensional, finite-element thermal-analysis program which handles conduction with internal heat generation, convection, radiation, specified flux, and specified temperature boundary conditions is presented. Elements used in the program are the triangle and tetrahedron for two- and three-dimensional analysis, respectively. The theory used in the program is developed, and several sample problems demonstrating the capability and reliability of the program are presented. A guide to using the program, description of the input cards, and program listing are included.

Heuser, J.↗

NASA NTRS1989

The coupling effects of kinematics and flexibility on the Lagrangian dynamic formulation of open chain deformable links

A nonlinear Lagrangian formulation for the spatial kinematic and dynamic analysis of open chain deformable links consisting of cylindrical joints that connect pairs of flexible links is developed. The special cases of revolute or prismatic joint can also be obtained from the kinematic equations. The kinematic equations are described using a 4x4 matrix method. The configuration of each deformable link in the open loop kinematic chain is identified using a coupled set of relative joint variables, constant geometric parameters, and elastic coordinates. The elastic coordinates define the link deformation with respect to a selected joint coordinate system that is consistent with the kinematic constraints on the boundary of the deformable link. These coordinates can be introduced using approximation techniques such as Rayleigh-Ritz method, finite element technique or any other desired approach. The large relative motion between two neighboring links are defined by a set of joint coordinates which describes the large relative translational and rotational motion between two neighboring joint coordinate systems. The origin of these coordinate systems are rigidly attached to the neighboring links at the joint definition points along the axis of motion.

Changizi, Koorosh↗

NASA NTRS1993

Compression buckling response of tailored rectangular composite plates

Buckling resistance is often a controlling criterion in the design of structures composed of plate elements. Design concepts that lead to increased buckling loads (or strains) of these plate elements can directly lower the structural cost and/or weight by a number of means. This study quantifies the improvements that can be achieved in compression buckling loads of rectangular composite plates by using a simple stiffness-tailoring concept. The approach is to position the unidirectional lamina through the thickness and over the planform of the plate so that the buckling load is increased with no loss in in-plane stiffness or increase in weight. Finite element analyses have been used to determine the effects of tailoring on the buckling loads of plates with various boundary conditions, aspect ratios, thicknesses, and membrane stiffnesses. Increases in buckling loads (or strains) of 200 percent or more compared to the uniform plate-buckling loads are shown possible with this tailoring concept.

Biggers, Sherrill B.↗

NASA NTRS1991

Improved compression buckling for rectangular composite plates by stiffness tailoring

Buckling resistance is often a controlling criterion in the design of structural elements. Design concepts that lead to increased buckling loads (or strains) can directly lower the structural cost and/or weight by a number of means. This study quantifies the improvements that can be achieved in compression buckling loads of rectangular composite plates by using a simple stiffness tailoring concept. The approach is to position the unidirectional lamina through the thickness and over the planform of the plate so that the buckling load is increased with no loss in in-plane stiffness or increase in weight. Finite element analyses have been used to determine the effects of tailoring on the buckling load of plates with various boundary conditions, aspect ratios, thicknesses, and membrane stiffnesses. Increases in buckling loads (or strains) of nearly 200 percent over the uniform plate buckling loads are shown possible with this tailoring concept.

Biggers, Sherrill B.↗

NASA NTRS1994

A Method for Thermal Analysis of Spiral Bevel Gears

A modeling method for analyzing the three-dimensional thermal behavior of spiral bevel gears has been developed. The model surfaces are generated through application of differential geometry to the manufacturing process for face-milled spiral bevel gears. Contact on the gear surface is found by combining tooth contact analysis with three-dimensional Hertzian theory. The tooth contact analysis provides the principle curvatures and orientations of the two surfaces. This information is then used directly in the Hertzian analysis to find the contact size and maximum pressure. Heat generation during meshing is determined as a function of the applied load, sliding velocity, and coefficient of friction. Each of these factors change as the point of contact changes during meshing. A nonlinear finite element program was used to conduct the heat transfer analysis. This program permitted the time- and position-varying boundary conditions, found in operation, to be applied to a one-tooth model. An example model and analytical results are presented.

Handschuh, Robert F.↗

NASA NTRS1992

An improved plate theory of order (1,2) for thick composite laminates

A new (1,2)-order theory is proposed for the linear elasto-static analysis of laminated composite plates. The basic assumptions are those concerning the distribution through the laminate thickness of the displacements, transverse shear strains and the transverse normal stress, with these quantities regarded as some weighted averages of their exact elasticity theory representations. The displacement expansions are linear for the inplane components and quadratic for the transverse component, whereas the transverse shear strains and transverse normal stress are respectively quadratic and cubic through the thickness. The main distinguishing feature of the theory is that all strain and stress components are expressed in terms of the assumed displacements prior to the application of a variational principle. This is accomplished by an a priori least-square compatibility requirement for the transverse strains and by requiring exact stress boundary conditions at the top and bottom plate surfaces. Equations of equilibrium and associated Poisson boundary conditions are derived from the virtual work principle. It is shown that the theory is particularly suited for finite element discretization as it requires simple C(sup 0)- and C(sup -1)-continuous displacement interpolation fields. Analytic solutions for the problem of cylindrical bending are derived and compared with the exact elasticity solutions and those of our earlier (1,2)-order theory based on the assumed displacements and transverse strains.

Tessler, A.↗

NASA NTRS1988

The effect of elastic boundary conditions on the dynamic response of rectangular plates

Natural frequencies and forced steady-state harmonic response for the vibration of uniform rectangular plates with edges elastically restrained against rotation and transverse translation are addressed. A single mode Rayleigh-Ritz solution is derived using functions that describe the normal modes of vibration of a beam whose ends are elastically restrained. The finite element solution is obtained for comparison. MACSYMA symbolic manipulation system is implemented as an aid to the mathematical rigor of the Ritz approach, and NASTRAN finite element code is used to model the mechanical system. Comparisons are made to published results and the solutions of this study are found to give lower frequencies for some values of boundary restraint. Steady-state harmonic amplitudes of displacement and acceleration are found to agree favorably for the two solutions. Low predictions of steady-state strain from NASTRAN result in some cases when compared to the Ritz values. Finally, a subjective assessment is made about the merit of using MACSYMA and NASTRAN.

Brewer, Terry K.↗

NASA NTRS1992

A probabilistic approach to the evaluation of fatigue damage in a space propulsion system injector element

The fatigue damage of a space propulsion system component is computed using probabilistic structural analysis methods. The analysis takes into account the variations in static and dynamic loads, the uncertainty in structural damping, and the scatter in material fatigue resistance. The key elements of the probabilistic approach include: (1) a numerical engine model for describing the global component interface loads consistent with engine balance, (2) models for computing the local mode boundary conditions on the component, (3) a static structural analysis model that captures the strains at a damage critical location as a function of engine performance variables, (4) a finite element model for assessment of the random amplitude stress response due to random base and pressure excitations with uncertain power and correlation length, and (5) advanced first order reliability methods for computing the probabilities associated with the fatigue damage.

Rajagopal, K. R.↗

NASA NTRS2003

Nonlinear Aeroacoustics Computations by the Space-Time CE/SE Method

The Space-Time Conservation Element and Solution Element Method, or CE/SE Method for short, is a recently developed numerical method for conservation laws. Despite its second order accuracy in space and time, it possesses low dispersion errors and low dissipation. The method is robust enough to cover a wide range of compressible flows: from weak linear acoustic waves to strong discontinuous waves (shocks). An outstanding feature of the CE/SE scheme is its truly multi-dimensional, simple but effective non-reflecting boundary condition (NRBC), which is particularly valuable for computational aeroacoustics (CAA). In nature, the method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its careful treatment of the surface fluxes and geometry, it is different from the existing schemes. Currently, the CE/SE scheme has been developed to a matured stage that a 3-D unstructured CE/SE Navier-Stokes solver is already available. However, in the present review paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen and sketched in section 2. Then applications of the 2-D and 3-D CE/SE schemes to linear, and in particular, nonlinear aeroacoustics are depicted in sections 3, 4, and 5 to demonstrate its robustness and capability.

Loh, Ching Y.↗

NASA NTRS1979

Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods

The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.

Periaux, J.↗

NASA NTRS1985

Algorithmic and theoretical results on computation of incompressible viscous flows by finite element methods

Primitive variable as well as streamfunction-vorticity and pure streamfunction formulations are discussed. For the primitive variable case alternative choices of the viscous stress term are shown to produce natural boundary conditions which are well suited for matching to various far field conditions. For the other cases recent analytical results, including error estimates are described, and an optical algorithm for pressure recovery as well as treatment for multiply connected domains are given.

Gunzburger, M. D.↗

NASA NTRS1974

Elastic-plastic analysis of a propagating crack under cyclic loading

Development and application of a two-dimensional finite-element analysis to predict crack-closure and crack-opening stresses during specified histories of cyclic loading. An existing finite-element computer program which accounts for elastic-plastic material behavior under cyclic loading was modified to account for changing boundary conditions - crack growth and intermittent contact of crack surfaces. This program was subsequently used to study the crack-closure behavior under constant-amplitude and simple block-program loading.

Newman, J. C., Jr.↗

NASA NTRS1984

Spectral methods for exterior elliptic problems

Spectral approximations for exterior elliptic problems in two dimensions are discussed. As in the conventional finite difference or finite element methods, the accuracy of the numerical solutions is limited by the order of the numerical farfield conditions. A spectral boundary treatment is introduced at infinity which is compatible with the infinite order interior spectral scheme. Computational results are presented to demonstrate the spectral accuracy attainable. Although a simple Laplace problem is examined, the analysis covers more complex and general cases.

Canuto, C.↗