Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS CALCULATION”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Thermal stress analysis of composites in the space environment

A finite element micromechanics approach was utilized to investigate the thermally induced stress fields in continuous fiber reinforced polymer matrix composites at temperatures typical of spacecraft operating environments. The influence of laminate orientation was investigated with a simple global/local formulation. Thermal stress calculations were used to predict probable damage initiation locations, and the results were compared to experimentally observed damage in several epoxy matrix composites. The influence of an interphase region on the interfacial stress states was investigated.

Bowles, David E.↗
NASA NTRS1974

Experimental evaluation of stresses in cylindrically hollow (drilled) balls

An experimental stress analysis was undertaken to evaluate stresses within cylindrically hollow (drilled) bearing balls proportioned for 40, 50, and 60 percent mass reductions. Strain gage rosettes were used to determine principal strains and stresses in the steel ball models statically loaded in various orientations. Results are reported for 127 mm OD balls under 44,500 N loads. Similitude considerations permit these results to be applied to calculate stresses in actual size drilled bearing balls proportioned to these mass reductions.

Nypan, L. J.↗
NASA NTRS1980

Effect of thermal cycling on ZrO2-Y2O3 thermal barrier coatings

The paper studies the comparative life of plasma-sprayed ZrO2-Y2O3 thermal barrier coatings on NiCrAlY bond coats on Rene 41 in short (4 min) and long (57 min) thermal cycles at 1040 C in a 0.3-Mach flame. Attention is given to determining the effect of short- and long-duration cycles on ZrO2-Y2O3 coatings, the cause of any cycle frequency effects, and methods to improve tolerance to thermal stress. Short cycles greatly reduced the life of the ceramic coating in terms of time at temperatures as compared to longer cycles, the failed coating indicating compressive failure. The experiments and stress calculations show that repeatedly subjecting a ceramic coating to high rates of initial heating has a more destructive influence on the coating than sustained operation at temperature. The effect of such thermal compressive stresses might be minimized through coating deposition and thickness control and by turbine cycle measurement to keep starting heating rates below critical values.

Mcdonald, G.↗
DOE OSTI2024

Optimal local truncation error method for 3-D elasticity interface problems

The paper deals with a new effective numerical technique on unfitted Cartesian meshes for simulations of heterogeneous elastic materials. Here, we develop the optimal local truncation error method (OLTEM) with 27- point stencils (similar to those for linear finite elements) for the 3-D time-independent elasticity equations with irregular interfaces. Only displacement unknowns at each internal Cartesian grid point are used. The interface conditions are added to the expression for the local truncation error and do not change the width of the stencils. The unknown stencil coefficients are calculated by the minimization of the local truncation error of the stencil equations and yield the optimal second order of accuracy for OLTEM with the 27-point stencils on unfitted Cartesian meshes. A new post-processing procedure for accurate stress calculations has been developed. Similar to basic computations it uses OLTEM with the 27-point stencils and the elasticity equations. The post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing numerical techniques (e.g., with finite elements). Numerical experiments show that at an accuracy of 0.1% for stresses, OLTEM with the new post-processing procedure significantly (by 10 5 -10 9 times) reduces the number of degrees of freedom compared to linear finite elements. OLTEM with the 27-point stencils yields even more accurate results than high-order finite elements with wider stencils.

42 ENGINEERING↗
NASA NTRS1994

Observations, models, and mechanisms of failure of surface rocks surrounding planetary surface loads

Geophysical models of flexural stresses in an elastic lithosphere due to an axisymmetric surface load typically predict a transition with increased distance from the center of the load of radial thrust faults to strike-slip faults to concentric normal faults. These model predictions are in conflict with the absence of annular zones of strike-slip faults around prominent loads such as lunar maria, Martian volcanoes, and the Martian Tharsis rise. We suggest that this paradox arises from difficulties in relating failure criteria for brittle rocks to the stress models. Indications that model stresses are inappropriate for use in fault-type prediction include (1) tensile principal stresses larger than realistic values of rock tensile strength, and/or (2) stress differences significantly larger than those allowed by rock-strength criteria. Predictions of surface faulting that are consistent with observations can be obtained instead by using tensile and shear failure criteria, along with calculated stress differences and trajectories, with model stress states not greatly in excess of the maximum allowed by rock fracture criteria.

Schultz, R. A.↗
DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗
DOE OSTI2017

BISON Fuel Performance Analysis of IFA-796 Rod 3 & 4 and Investigation of the Impact of Fuel Creep

In order to improve the accident tolerance of light water reactor (LWR) fuel, alternative cladding materials have been proposed to replace the currently used zirconium (Zr)-based alloys. Of these materials, there is a particular focus on iron-chromiumaluminum (FeCrAl) alloys because they exhibit slower oxidation kinetics in high-temperature steam than Zr-alloys. This should decrease the energy release due to oxidation and slow cladding consumption in the presence of high temperature steam. These alloys should also exhibit increased “coping time” in the event of an accident scenario by improving the mechanical performance at high temperatures, allowing greater flexibility to achieve core cooling. As a continuation of the development of these alloys, in-reactor irradiation testing of FeCrAl cladded fuel rods has started. In order to provide insight on the possible behavior of these fuel rods as they undergo irradiation in the Halden Boiling Water Reactor, engineering analysis has been performed using FeCrAl material models implemented into the BISON fuel performance code. This milestone report provides an update on the ongoing development of modeling capability to predict FeCrAl cladding fuel performance and to provide an early look at the possible behavior of planned in-reactor FeCrAl cladding experiments. In particular, this report consists of two separate analyses. The first analysis consists of fuel performance simulations of IFA-796 rod 4 and two segments of rod 3. These simulations utilize previously implemented material models for the C35M FeCrAl alloy and UO2 to provide a bounding behavior analysis corresponding to variation of the initial fuel cladding gap thickness within the fuel rod. The second analysis is an assessment of the fuel and cladding stress states after modification of the fuel creep model that is currently implemented in the BISON fuel performance code. Effects from modifying the fuel creep model were identified for the BISON simulations of the IFA-796 rod 4 experiment, but show that varying the creep model (within the range investigated here) only provide a minimal increase in the fuel radius and maximum cladding hoop stress. Continued investigation of fuel behavioral models will include benchmarking the modified fuel creep model against available experimental data, as well as an investigation of the role that fuel cracking will play in the compliance of the fuel. Correctly calculating stress evolution in the fuel is key to assessing fuel behavior up to gap closure and the subsequent deformation of the cladding due to PCMI. The inclusion of frictional contact should also be investigated to determine the axial elongation of the fuel rods for comparison with data from this experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2001

Implementing a Nitrogen-Based Model for Autotrophic Respiration Using Satellite and Field Observations

The rate of carbon accumulation by terrestrial plant communities in a process-level, mechanistic modeling is the difference of the rate of gross photosynthesis by a canopy (A(sub g)) and autotrophic respiration (R) of the stand. Observations for different biomes often show that R to be a large and variable fraction of A(sub g), ca. 35% to 75%, although other studies suggest the ratio of R and A(sub g) to be less variable. Here, R has been calculated according to the two compartment model as being the sum of maintenance and growth components. The maintenance respiration of foliage and living fine roots for different biomes has been determined objectively from observed nitrogen content of these organs. The sapwood maintenance respiration is based on pipe theory, and checked against an independently derived equation considering sapwood biomass and its maintenance coefficient. The growth respiration has been calculated from the difference of A(sub g) and maintenance respiration. The A(sub g) is obtained as the product of biome-specific radiation use efficiency for gross photosynthesis under unstressed conditions and intercepted photosynthetically active radiation, and adjusted for stress. Calculations have been done using satellite and ground observations for 36 consecutive months (1987-1989) over large contiguous areas (ca. 10(exp 5) sq km) of boreal forests, crop land, temperate deciduous forest, temperate grassland, tropical deciduous forest, tropical evergreen forest, tropical savanna, and tundra. The ratio of annual respiration and gross photosynthesis, (R/A(sub g)), is found to be 0.5-0.6 for temperate and cold adopted biome areas, but somewhat higher for tropical biome areas (0.6-0.7). Interannual variation of the fluxes is found to be generally less than 15%. Calculated fluxes are compared with observations and several previous estimates. Results of sensitivity analysis are presented for uncertainties in parameterization and input data. It is found that uncertainty in determining maintenance respiration for tropical biomes is such that R/A(sub g) for these biomes could be similar to that for temperate biomes.

Choudhury, Bhaskar J.↗
NASA NTRS1991

Stress intensity factors for large aspect ratio surface and corner cracks at a semi-circular notch in a tension specimen

Stress intensity factor solutions for semielliptic surface and quarter-elliptic corner cracks emanating from a semicircular notch in a tension specimen are presented. A three-dimensional finite-element analysis in conjunction with the equivalent domain integral was used to calculate stress intensity factors (SIF). SIF solutions for surface or corner crack (crack length to depth ratio of 2) at a notch are presented for a wide range of crack sizes and notch radii. Results showed that the SIF are larger for larger crack lengths and for larger notch radii. The SIF are nearly constant all along the crack front for deep surface cracks and for all corner cracks analysed.

Shivakumar, K. N.↗
DOE OSTI2017

Einstein gravity 3-point functions from conformal field theory

We report on stress tensor correlation functions in four-dimensional conformal field theories with large N and a sparse spectrum. Theories in this class are expected to have local holographic duals, so effective field theory in anti-de Sitter implies that the stress tensor sector should exhibit universal, gravity-like behavior. At the linearized level, the hallmark of locality in the emergent geometry is that stress tensor three-point functions $\langle$T T T $\rangle$, normally specified by three constants, should approach a universal structure controlled by a single parameter as the gap to higher spin operators is increased. We demonstrate this phenomenon by a direct CFT calculation. Stress tensor exchange, by itself, violates causality and unitarity unless the three-point functions are carefully tuned, and the unique consistent choice exactly matches the prediction of Einstein gravity. Under some assumptions about the other potential contributions, we conclude that this structure is universal, and in particular, that the anomaly coefficients satisfy a ≈ c as conjectured by Camanho et al. The argument is based on causality of a four-point function, with kinematics designed to probe bulk locality, and invokes the chaos bound of Maldacena, Shenker, and Stanford.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

A boundary element alternating method for two-dimensional mixed-mode fracture problems

A boundary element alternating method, denoted herein as BEAM, is presented for two dimensional fracture problems. This is an iterative method which alternates between two solutions. An analytical solution for arbitrary polynomial normal and tangential pressure distributions applied to the crack faces of an embedded crack in an infinite plate is used as the fundamental solution in the alternating method. A boundary element method for an uncracked finite plate is the second solution. For problems of edge cracks a technique of utilizing finite elements with BEAM is presented to overcome the inherent singularity in boundary element stress calculation near the boundaries. Several computational aspects that make the algorithm efficient are presented. Finally, the BEAM is applied to a variety of two dimensional crack problems with different configurations and loadings to assess the validity of the method. The method gives accurate stress intensity factors with minimal computing effort.

Raju, I. S.↗
DOE OSTI2022

Residual Stress Measurements in Extreme Environments for Hazardous, Layered Specimens

In nuclear fuel plates of low-enriched U-10Mo (LEU) clad with aluminum by hot isostatic pressing (HIP), post-irradiation stresses arising during reactor shutdown are a major concern for safe reactor operations. Measurement of those residual stresses has not previously been possible because the high radioactivity of the plates requires handling only by remote manipulation in a hot cell. The incremental slitting method for measuring through-thickness stress profiles was modified, and a system for automated, remote operation was built and tested. Experimental modifications consisted of replacing electric-discharge machining (EDM) with a small end mill and strain-gauge measurements with cantilever displacement measurements. The inverse method used to calculate stresses was the pulse-regularization method modified to allow discontinuities across material interfaces. The new system was validated by comparing with conventional slitting on a depleted U-10Mo (DU) fuel plate. The new system was applied to two measurements each on six as-fabricated (pre-irradiation) LEU miniature fuel plates. Variations between the measurements at two locations in the same plate were strongly correlated with measured geometrical heterogeneity in the plate—a tilt in the fuel foil. Compressive stresses in the U-10Mo were shown to increase from 20 to 250 MPa as the ratio of aluminum thickness to U-10Mo thickness increased causing increased constraint during cooling. Faster cooling rates during processing also increased stress magnitudes. The measurements trends agreed with data in the literature from similar plates made with DU, which further validates the method. Because other methods are impractical in a hot cell, the modified slitting method is now poised for the first measurements of post-irradiation stresses.

36 MATERIALS SCIENCE↗
NASA NTRS2001

C-Coupon Studies of CMCS: Fracture Behavior and Microstructural Characterization

A curved beam 'C-coupon' was used to assess fracture behavior in a Sylramic(tm)/melt infiltration (MI) SiC matrix composite. Failure stresses and fracture mechanisms, as determined by optical and scanning electron microstructural analysis, are compared with finite element stress calculations to analyze failure modes. Material microstructure was found to have a strong influence on mechanical behavior. Fracture occurs in interlaminar tension (ILT), provided that the ratio of ILT to tensile strength for the material is less than the ratio of radial to hoop stresses for the C-coupon geometry. Utilization of 3D architectures to improve interlaminar strength requires significant development efforts to incorporate through thickness fibers in regions with high curvatures while maintaining uniform thickness, radius, and microstructure.

Hurwitz, Frances I.↗
DOE OSTI2023

Methods for Estimating Hydrogen Fuel Tank Characteristics

The pressure vessels needed to store hydrogen for next-generation hydrogen fuel cell vehicles are expected to be a substantial portion of the total system mass, volume, and cost. Gravimetric capacity, volumetric capacity, and cost per kilogram of usable hydrogen are key performance metrics that the U.S. Department of Energy (DOE) uses to determine the viability of hydrogen fuel cell systems. Research and development related to hydrogen storage systems covers a wide range of potential operating conditions, from cryogenic temperatures to high temperatures (above ambient) and low pressure to high pressure. Researchers at PNNL have developed methods for estimating these key pressure vessel characteristics to support on-board hydrogen storage system design and performance evaluation and to support decision-making about DOE hydrogen storage system research investments. This article describes the pressure tank estimation methodology that has been used as a stand-alone calculation and has been incorporated into larger system evaluation tools. The methodology estimates the geometry, mass, and material cost of type I, type III, and type IV pressure vessels based on operating pressure and material strength at the system's operating temperature, using classical thin-wall and thick-wall pressure vessel stress calculations. The geometry, mass, and material cost requirements of the pressure vessel have significant impacts on the total system performance. For example, hydrogen storage materials that can separately achieve a very high hydrogen density can be deemed impractical for use in fuel cell vehicle hydrogen storage systems because the pressure tank containing them is too large, heavy, or expensive. This article describes the design philosophy and analytical process of the tank characteristic estimation methodology, which has been implemented in spreadsheet calculation tools and system-level analysis tools used by DOE researchers. Each of the three tank types (type I, type III, and type IV) uses a different analysis methodology with some common elements. This article also provides examples of implementing the methodology to perform parametric studies of all three pressure vessel types. The goal of this article is to present the methodology in sufficient detail so it can be implemented in other hydrogen fuel cell vehicle design and analysis tools.

42 ENGINEERING↗
DOE OSTI2020

GPU acceleration of all-electron electronic structure theory using localized numeric atom-centered basis functions

We present an implementation of all-electron density-functional theory for massively parallel GPU-based platforms, using localized atom-centered basis functions and real-space integration grids. Special attention is paid to domain decomposition of the problem on non-uniform grids, which enables compute- and memory-parallel execution across thousands of nodes for real-space operations, e.g. the update of the electron density, the integration of the real-space Hamiltonian matrix, and calculation of Pulay forces. To assess the performance of our GPU implementation, we performed benchmarks on three different architectures using a 103-material test set. We find that operations which rely on dense serial linear algebra show dramatic speedups from GPU acceleration: in particular, SCF iterations including force and stress calculations exhibit speedups ranging from 4.5 to 6.6. For the architectures and problem types investigated here, this translates to an expected overall speedup between 3–4 for the entire calculation (including non-GPU accelerated parts), for problems featuring several tens to hundreds of atoms. Additional calculations for a 375-atom Bi2Se3 bilayer show that the present GPU strategy scales for large-scale distributed-parallel simulations.

42 ENGINEERING↗
DOE OSTI2016

Investigation of wellbore microannulus permeability under stress via experimental wellbore mock-up and finite element modeling

This research aims to describe the microannulus region of the cement sheath-steel casing interface in terms of its compressibility and permeability. Here, a wellbore system mock-up was used for lab-scale testing, and was subjected to confining and casing pressures in a pressure vessel while measuring gas flow along the specimen’s axis. The flow was interpreted as the hydraulic aperture of the microannuli. Numerical joint models were used to calculate stress and displacement conditions of the microannulus region, where the mechanical stiffness and hydraulic aperture were altered in response to the imposed stress state and displacement across the joint interface.

58 GEOSCIENCES↗
DOE OSTI2023

FINITE ELEMENT MODEL MESH REFINEMENT EFFECTS ON QUALIFICATION OF NUCLEAR GRADE GRAPHITE COMPONENTS

The American Society of Mechanical Engineers (ASME) provides the full and simplified design-by-analysis probabilistic assessments for determining acceptance of nuclear grade graphite core components. The assessments can be characterized by three parts: (1) a component stress distribution, often determined by a finite element (FE) model; (2) a Weibull probability density function (pdf) that characterizes the experimental tensile strength distribution; and (3) the post-processor, which combines the FE model and the Weibull strength distribution in accordance with the full and simplified assessments to determine component acceptance. It is known that the level of mesh refinement in FE models can affect the modeled component’s calculated stress distribution. Depending on the component geometry, the stress distribution may converge with sufficient refinement. It was previously unknown whether the acceptance decision resulting from the full and simplified assessments might change even with sufficient mesh refinement. This study explores that question using experimental strength results for a dog-bone geometry for two graphite grades, IG-110 and PCEA. The simplified assessment has two criteria that must be met, the first limits the combined membrane stress by the allowable stress and the second limits the peak equivalent stress by the allowable stress scaled by the ratio of flexural to tensile strength. In the application of the simplified assessment, convergence of the peak equivalent stress required extreme mesh refinement, however, the acceptance decision was not affected. It is hypothesized that more complex geometries with stress concentrations may present mesh refinement effects on the simplified assessment acceptance decision. Mesh refinement did affect the acceptance decision in the full assessment for the applied pressure loadings in this study. This work suggests component stress distribution convergence is not a sufficient criteria for POF convergence in the full assessment and that mesh refinement should continue until the POF has converged, especially where the resulting POF is bordering the SRC acceptable POF limit.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗