Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element expansion”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Residual Stresses Modeled in Thermal Barrier Coatings

Thermal barrier coating (TBC) applications continue to increase as the need for greater engine efficiency in aircraft and land-based gas turbines increases. However, durability and reliability issues limit the benefits that can be derived from TBC's. A thorough understanding of the mechanisms that cause TBC failure is a key to increasing, as well as predicting, TBC durability. Oxidation of the bond coat has been repeatedly identified as one of the major factors affecting the durability of the ceramic top coat during service. However, the mechanisms by which oxidation facilitates TBC failure are poorly understood and require further characterization. In addition, researchers have suspected that other bond coat and top coat factors might influence TBC thermal fatigue life, both separately and through interactions with the mechanism of oxidation. These other factors include the bond coat coefficient of thermal expansion, the bond coat roughness, and the creep behavior of both the ceramic and bond coat layers. Although it is difficult to design an experiment to examine these factors unambiguously, it is possible to design a computer modeling "experiment" to examine the action and interaction of these factors, as well as to determine failure drivers for TBC's. Previous computer models have examined some of these factors separately to determine their effect on coating residual stresses, but none have examined all the factors concurrently. The purpose of this research, which was performed at DCT, Inc., in contract with the NASA Lewis Research Center, was to develop an inclusive finite element model to characterize the effects of oxidation on the residual stresses within the TBC system during thermal cycling as well as to examine the interaction of oxidation with the other factors affecting TBC life. The plasma sprayed, two-layer thermal barrier coating that was modeled incorporated a superalloy substrate, a NiCrAlY bond coat, and a ZrO2-8 wt % Y2O3 ceramic top coat. We examined the effect on stress during burner rig thermal cycling of the following independent variables: creep in the bond coat and top coat, oxidation, bond coat coefficient of thermal expansion, number of thermal cycles, and interfacial roughness. All these factors were suspected of influencing TBC failure. The model showed that all the material properties studied had a significant effect on the coating's residual stresses if the interface of the bond coat was rough. Bond coat expansion, bond coat oxidation, and bond coat creep had the highest effect on coating stresses, and these were highly interactive. The model also showed that the mechanism of stress generation during thermal cycling changed with the number of thermal cycles. Bond coat and top coat creep dominated stress generation during early thermal cycles, greatly increasing delamination stresses at the peaks of the bond coat. Therefore, creep is the prime driver for delamination cracking early in life, but cracking is limited to the bond coat peak region. Oxidation of the bond coat, on the other hand, tended to dominate stress generation during later cycles by greatly increasing delamination stresses over bond coat valleys. These results indicate that oxidation is the driver for the continued cracking necessary to cause ceramic layer spallation.

Freborg, A. M.↗
NASA NTRS1987

Dynamic characteristics of a vibrating beam with periodic variation in bending stiffness

A detailed dynamic analysis is performed of a vibrating beam with bending stiffness periodic in the spatial coordinate. Using a perturbation expansion technique the free vibration solution is obtained in a closed-form, and the effects of system parameters on beam response are explored. It is found that periodic stiffness acts to modulate the modal displacements from the characteristic shape of a simple sine wave. The results are verified by a finite element solution and through experimental testing.

Townsend, John S.↗
NASA NTRS1984

Three-dimensional Navier-Stokes calculations of multiple interacting vortex rings

Results from a finite-difference Navier-Stokes code for three-dimensional, unsteady, vortical flows in unbounded domains are presented and analyzed in this paper. The vortical flows presented are representative of vortex rings and other closed vortical tubes or structures in fluid mechanics. Such structures are important elements in fluid flows such as jets, atmospheric turbulence, and the far-field wakes of aircraft, and studies of their interaction may aid in an understanding of complex fluid flows. The paper demonstrates that computational methods can be used as a viable alternative or supplement to experimental techniques for studying the physics of vortex flows. The separate visualization of vortex stretching, convection, and diffusion is presented in this paper for a single elliptical vortex ring.The calculations employ a truncated series expansion technique to simulate the unbounded nature of the fluid flow with a finite computational domain, which is a more accurate technique than the conventional freestream boundary specification. The numerical divergence of the three-dimensional vorticity field is considered as a useful estimate of truncation error, and the use of a kinetic energy decay law as a calculation check is demonstrated. Results from the Navier-Stokes code are presented for the unsteady motion of two and four vortex rings along parallel axes, and the results agree qualitatively with experimental flow visualization.

Chamberlain, J. P.↗
NASA NTRS1991

Residual thermal strains and stresses in nickel aluminide matrix composites

Thermally induced residual strains and stresses developed during postfabrication cooling in Saphikon/NiAl and tungsten/NiAl high-temperature composites are investigated through three-dimensional elastoplastic finite-element analyses. Average axial and transverse strains in the matrix are found to be tensile and compressive, respectively, and similar for both Saphikon and W-fiber-reinforced NiAl composites. It is suggested that the residual matrix stresses and strains are controlled more by the low-matrix yield stress than by the fiber/matrix expansion mismatch. Residual thermal strains in the matrix of these composites are measured by using a neutron-diffraction technique; the measured axial and transverse strains in the matrix are found to be in agreement with the computed values.

Saigal, A.↗
NASA NTRS2018

Advanced Mirror Technology Development (AMTD) II Modal Test of A 1.5 m Glass Slumped Mirror

The Advanced Mirror Technology Development (AMTD) project is in Phase 2 of a multiyear effort initiated in Fiscal Year 2012 to mature toward the next technology readiness level critical technologies required to enable 4-m-or-larger monolithic or segmented ultraviolet, optical, and infrared (UVOIR) space telescope primary-mirror assemblies for general astrophysics and ultra-high-contrast observations of exoplanets. As part of AMTD II, a free-free modal test was performed of a light weighted slumped 1.5 m mirror made of Corning Ultra Low Expansion (ULE®) material. The test article and support structure were suspended via bungee to simulate a free-free environment. Modes were excited by roaming an instrumented modal test hammer and responses were measured. Predicted and measured frequencies are presented as well as Modal Assurance Criteria (MAC) results to compare the mode shapes. The finite element mirror model used for pre-test predictions and posttest comparisons was provided by the mirror vendor, Harris Corporation. The mirror FEM included deformations of the ribs that were a result of the slumping process. Modal test frequencies matched predictions within the 5% target with the exception of one mode and that pair differed by 5.2%. Of the seven modes measured and predicted, four had MAC values meeting the target of ≥ 0.90, one was just under and two were notably below the target.

Knight, J. Brent↗
NASA NTRS1984

The solution of non-linear hyperbolic equation systems by the finite element method

A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.

Loehner, R.↗
NASA NTRS1984

Investigation of Toughened Neat Resins and Their Relations to Advanced Composite Mechanical Properties

Four epoxies (Hercules 3502, Hercules 2220-1 and 2220-3, and Ciba-Geigy Fibredux 914 epoxy) were evaluated to judge their merits within the ACEE program and also to provide material property input for a micromechanics finite element program. A casting procedure was developed for the four matrix systems mentioned and mechanical properties were studied at three test temperatures and two moisture contents. A scanning electron microscope study was performed to catalogue fracture surfaces for the unreinforced (neat) epoxies. In addition, tension testing, shear testing, a limited amount of fracture toughness testing, and thermal expansion coefficient and moisture expansion coefficient testing were carried out. Mechanical material properties were generated for input into a micromechanics program and to characterize these four matrix systems.

R. S. Zimmerman↗
NASA NTRS2000

Finite-Element Modeling of 3C-SiC Membranes

Finite-element modeling (FEM) of 3C-SiC thin-film membranes on Si substrates was used to determine the residual stress and center deflection with applied pressure. The anisotropic, three-dimensional model includes the entire 3C-SiC membrane and Si substrate with appropriate material properties and boundary conditions. Residual stress due to the thermal-expansion-coefficient mismatch between the3C-SiC film and Si substrate was included in the model. Both before-and after-etching, residual stresses were calculated. In-plane membrane stress and normal deflection with applied pressure were also calculated. FEM results predict a tensile residual stress fo 259 MPa in the 3C-SiC membrane before etching. This decreases to 247 MPa after etching the substrate below the membrane. The residual stress experimentally measured on sample made at Case Western Reserve University was 280 MPa on post-etched membranes. This is excellent agreement when an additional 30-40 MPa of residual stress to account for lattice mismatch is added to the FEM results.

DeAnna, R. G.↗
NASA NTRS1983

Parallel solution of finite element equations

The paper examines several parallel processing solution algorithms for finite element equations arising in linear equilibrium problems. Two basic groups of algorithms, direct and iterative, are investigated with respect to a number of parallel computer architectures and associated selection criteria. The direct algorithms include: LR-Gauss, Crout, Cholesky, Cyclic Reduction and WZ-factorization. The iterative methods examined are: Accelerated Gauss-Seidel, Surrogate Stiffness, Jacobi, Series Expansion, and Energy Monte Carlo. For real-time applications, where the object is to minimize the execution time, Cyclic Reduction appears to be best suited. This assumes a computer with an unlimited number of parallel processors. However, for computers with a limited number of parallel processors that must be used efficiently, both Gauss factorization and Jacobi-like iterative methods rank favorably.

Salama, M.↗
NASA NTRS2018

Simulation of Stochastic Mud-Crack Damage Formation in an Environmental Barrier Coating

The integrated Finite Element Analysis–Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button–sized disk–shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300oC to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.

residual tensile stress↗
NASA NTRS1991

Scattering of Lamb waves in a composite plate

A combined analytical and finite element technique is developed to gain a better understanding of the scattering of elastic waves by defects. This hybrid method is capable of predicting scattered displacements from arbitrary shaped defects as well as inclusions of different material. The continuity of traction and displacements at the boundaries of the two areas provided the necessary equations to find the nodal displacements and expansion coefficients. Results clearly illustrate the influence of increasing crack depth on the scattered signal.

Bratton, Robert↗
NASA NTRS1999

Modeling of SiC Lateral Resonant Devices Over a Broad Temperature Range

Finite-element analysis (FEA) modal results of 3C-SiC lateral resonant devices anchored to a Si substrate are presented as resonant frequency versus temperature. The suspended elements are etched from a 2 micron, 3C-SiC film grown at 1600 K on a 500 micron-thick, Si substrate. The analysis includes, temperature-dependent properties, shape change due to volume expansion with temperature, and thermal stress caused by differential thermal expansion of different materials. Two designs are considered: type I has anchor locations close to the geometric centroid and a small shuttle; type 11 has a large shuttle with anchors far from the centroid, The resonant frequency decreases approximately 3.5% over a 1000 K temperature increase for the type-I device, and behaves according to theory. The resonant frequency of the type-11 device decreases by 2% over the first 400 K, then rises slightly over the remaining 600 K. This device deviates from theory because of the high thermal stress induced in the beams. The thermal stress is caused by the differential thermal expansion of the suspended element relative to the substrate. The results show that the device geometry must be properly chosen if the resonant frequency of that device will be used to calculate the temperature coefficient of Young's modulus. These results apply only to resonators of one material on a substrate of a different material.

DeAnna, Russell G.↗
DOE OSTI2025

Boosting efficiency and reducing graph reliance: Basis adaptation integration in Bayesian multi-fidelity networks

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

42 ENGINEERING↗
NASA NTRS1983

An analytical investigation of shape control of large space structures by applied temperatures

An analytical procedure for the static shape control of flexible space structures subjected to thermal distortions is developed which is based on prescribing temperatures in control elements having much higher coefficients of thermal expansion than the main structure. The temperatures at the control elements are defined so as to minimize the overall thermal distortion of the structure from its ideal shape, and a matrix equation is obtained which can be solved for the set of optimum control temperatures. A formulation of the procedure for continuous structures governed by differential equations and a formulation for discrete (finite element modeled) structures governed by matrix equations are presented. The equations from the continuous formulation are employed for the shape control of a simple beam distorted by nonuniform heating, and the discrete formulation is applied in a general purpose finite-element structural analysis computer program for the shape control of a 750 m radiometer antenna reflector dish subjected to orbital heating. A reduction in thermal distortion by a factor of nearly 50 was obtained with the use of only seven control elements. Results for four different sets of control locations for the antenna are presented in which reductions in distortion of up to a factor of four were obtained.

Haftka, R. T.↗
NASA NTRS2020

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗
NASA NTRS1990

Micromechanical thermal analysis of interphase region in a titanium aluminide MMC

The high reactivity between the fiber and matrix in silicon carbide/titanium aluminide MMCs leads to the formation of brittle reaction products at the fiber/matrix (F/M) interface. Also, the high thermal expansion coefficient mismatch between the fiber and matrix leads to high tensile residual stresses at the F/M interface, and this can lead to premature cracking during cooldown. One solution to these problems is the use of a metallic fiber coating like Ta which acts as an F/M diffusion barrier and reacts with the matrix to form a beta stabilized compliant layer. A finite element micromechanics analysis was performed to study the effects of Ta and beta interphase layers on the thermal residual stresses during consolidation. A 5-micron-thick beta layer reduced cool-down stresses by 8 percent compared to a 2 percent reduction computed for a Ta layer of the same thickness. Plastic yielding in the Ta was not effective in reducing cool-down stresses. Compliant alpha-2 particles next to the stiffer gamma particles reduced stresses in the gamma particles by less than 2 percent. A simple closed form analysis was developed to calculate thermal residual stresses in a fiber/interphase/matrix system. A 2-micron-thick Ag interphase layer was found to reduce residual stresses by about 11 percent.

Naik, R. A.↗
NASA NTRS1978

Evaluation and expansion of an analytical model for fatigue of notched composite laminates

The analytical and experimental study performed to expand the existing static and fatigue failure analysis is described. The analytical effort extended the analysis to include interlaminar effects, while the experimental effort developed methods to obtain basic experimental data required as input to the analysis. The static failure analysis for notched laminates was modified to include interlaminar effects near the notch. Three dimensional elastic and two dimensional elasticplastic finite element analysis were performed for some notched laminates.

Ramkumar, R. L.↗
NASA NTRS2001

GRCop-84 Developed for Rocket Engines

GRCop-84 (Cu-8 at.% Cr-4 at.% Nb) was developed at the NASA Glenn Research Center for use in regeneratively cooled rocket engines. This alloy possesses an excellent combination of conductivity, thermal expansion, strength, creep resistance, ductility, and low-cycle fatigue (LCF) life. In comparison to the currently used alloy, NARloy-Z (Cu-3 wt % Ag-0.5 wt % Zr), GRCop-84's significantly better properties give it the potential for significant gains in engine performance and reliability. For GRCop-84 to be safely used in a rocket engine, it is critical to develop a detailed data base of the important thermophysical and mechanical properties. Work has focused on five major properties: thermal conductivity, thermal expansion, tensile strength, creep resistance, and LCF life. The analysis went beyond simply measuring the data and reporting averages. A detailed statistical analysis was conducted that allowed regression of the data over the entire temperature range tested and establishment of design minimums. The design values were expressed as simple mathematical formulas that are highly amenable to use in computer codes such as in finite element analysis and related computer modeling work.

Ellis, David L.↗