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 127 records · Page 7

Non-uniform Stress-free Strains in a Spherically Symmetrical Nano-sized Particle and Its Applications to Lithium-ion Batteries

The stress-free strain originated from local chemical composition and phase transformation can significantly alter the microstructures of materials; and then affect their properties. In this paper, we developed an analytical method to calculate stress-strain field due to the non-uniform stress-free strain in a spherically symmetrical particle. Applying the method to a lithium ion (Li-ion) battery electrode, the evolution of Li-ion concentration and strain field during the lithiation process is studied. Our studies reveal that the maximum strain in the electrode generally occurs on surface of sample, and is mainly dependent on the difference of Li-ion concentration of surface and of center in sample. Decreasing the difference of Li-ion concentration can efficiently decrease the maximum strain so that cracks of electrodes can been prevented. Our analytical results provide a useful guidance for practical applications of energy storage materials.

36 MATERIALS SCIENCE↗

Stress Mapping in Glass-to-Metal Seals using Indentation Crack Lengths.

Predicting the residual stress which develops during fabrication of a glass-to-metal compression seal requires material models that can accurately predict the effects of processing on the sealing glass. Validation of the predictions requires measurements on representative test geometries to accurately capture the interaction between the seal materials during a processing cycle required to form the seal, which consists of a temperature excursion through the glass transition temperature of the sealing glass. To this end, a concentric seal test geometry, referred to as a short cylinder seal, consisting of a stainless steel shell enveloping a commercial sealing glass disk has been designed, fabricated, and characterized as a model validation test geometry. To obtain data to test/validate finite element (FE) stress model predictions of this geometry, spatially-resolved residual stress was calculated from the measured lengths of the cracks emanating from radially positioned Vickers indents in the glass disk portion of the seal. The indentation crack length method is described, and the spatially-resolved residual stress determined experimentally are compared to FE stress predictions made using a nonlinear viscoelastic material model adapted to inorganic sealing glasses and an updated rate dependent material model for 304L stainless steel. The measurement method is a first to achieve a degree of success for measuring spatially resolved residual stress in a glass-bearing geometry and a favorable comparison between measurements and simulation was observed.

36 MATERIALS SCIENCE↗

Stress Mapping in Glass-to-Metal Seals using Indentation Crack Lengths

Predicting the residual stress which develops during fabrication of a glass-to-metal compression seal requires material models that can accurately predict the effects of processing on the sealing glass. Validation of the predictions requires measurements on representative test geometries to accurately capture the interaction between the seal materials during a processing cycle required to form the seal, which consists of a temperature excursion through the glass transition temperature of the sealing glass. To this end, a concentric seal test geometry, referred to as a short cylinder seal, consisting of a stainless steel shell enveloping a commercial sealing glass disk has been designed, fabricated, and characterized as a model validation test geometry. To obtain data to test/validate finite element (FE) stress model predictions of this geometry, spatially-resolved residual stress was calculated from the measured lengths of the cracks emanating from radially positioned Vickers indents in the glass disk portion of the seal. The indentation crack length method is described, and the spatially-resolved residual stress determined experimentally are compared to FE stress predictions made using a nonlinear viscoelastic material model adapted to inorganic sealing glasses and an updated rate dependent material model for 304L stainless steel. The measurement method is a first to achieve a degree of success for measuring spatially resolved residual stress in a glass-bearing geometry and a favorable comparison between measurements and simulation was observed.

36 MATERIALS SCIENCE↗

Comparative energetics of carbon storage molecules in green algae

Several members of the green algae possess the ability to produce lipids and/or high value compounds in significant quantities. While for several of these green algal species induction of increased lipid production has been shown, and cultivation of species for high value molecules occurs at production scale, the molecular mechanisms governing over-accumulation of molecules synthesized from isoprenoid precursors, carotenoids, for example, have received far less attention. Here, we present a calculation of the required ATP equivalencies per carbon atom and reducing power equivalencies as NADH/NADPH (NAD(P)H) per carbon atom for the isoprenoid molecules ..beta..-carotene (C40), astaxanthin (C40), and squalene (C30). We compared energetic requirements of carbohydrates, triacylglycerol, and isoprenoid molecules under a gradient of conditions of cellular stress. Our calculations revealed slightly less ATP and NAD(P)H equivalency per carbon atom between triacylglycerol and the three isoprenoid molecules. Based on our results, we propose that the driving force for differences in accumulation patterns of carotenoids vs. triacylglycerols in algal cells under stress is largely dependent on the presence and regulation of bypass mechanisms at metabolic junction bottlenecks, like pyruvate dehydrogenase (PDH), within particular species. We provide a discussion of several molecular mechanisms that may influence carbon partitioning within different groups of green algae, including metabolic inhibition through accumulation of specific substrates related to ATP and reducing equivalent production (NAD(P)H) as well as cellular compartmentalization. This work contributes to the ongoing discussion of cellular homeostatic regulation during stress, as well as the potential mechanisms driving long-term carbon storage as it relates to energy and redox states within the algal cell.

09 BIOMASS FUELS↗

Combined shape and topology optimization for minimization of maximal von Mises stress

Here, this work shows that a combined shape and topology optimization method can produce optimal 2D designs with minimal stress subject to a volume constraint. The method represents the surface explicitly and discretizes the domain into a simplicial complex which adapts both structural shape and topology. By performing repeated topology and shape optimizations and adaptive mesh updates, we can minimize the maximum von Mises stress using the p-norm stress measure with p-values as high as 30, provided that the stress is calculated with sufficient accuracy.

42 ENGINEERING↗

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗

Mitigating cutting-induced plasticity in the contour method, Part 2: Numerical analysis

Cutting-induced plasticity can have a significant effect on the measurement accuracy of the contour method. The present study examines the benefit of a double-embedded cutting configuration that relies on self-restraint of the specimen, relative to conventional edge-crack cutting configurations. A series of finite element analyses are used to simulate the planar sectioning performed during double-embedded and conventional edge-crack contour cutting configurations. The results of numerical analyses are first compared to measured results to validate the cutting simulations. The simulations are then used to compare the efficacy of different cutting configurations by predicting the deviation of the residual stress profile from an original (pre-cutting) reference stress field, and the extent of cutting-induced plasticity. Comparisons reveal that while the double-embedded cutting configuration produces the most accurate residual stress measurements, the highest levels of plastic flow are generated in this process. As a result, this cutting-induced plastic deformation is, however, largely confined to small ligaments formed as a consequence of the sample sectioning process, and as such it does not significantly affect the back-calculated residual stress field.

36 MATERIALS SCIENCE↗

Computational Modeling of Heterogeneity of Stress, Charge, and Cyclic Damage in Composite Electrodes of Li-Ion Batteries

Charge heterogeneity is a prevalent feature in many electrochemical systems. In a commercial cathode of Li-ion batteries, the composite is hierarchically structured across multiple length scales including the sub-micron single-crystal primary-particle domains up to the macroscopic particle ensembles. The redox kinetics of charge transfer and mass transport strongly couples with mechanical stresses. This interplay catalyzes substantial heterogeneity in the charge (re)distribution, stresses, and mechanical damage in the composite electrode during charging and discharging. We assess the heterogeneous electrochemistry and mechanics in a LiNi x Mn y Co z O 2 (NMC) cathode using a fully coupled electro-chemo-mechanics model at the cell level. A microstructure-resolved model is constructed based on the synchrotron X-ray tomography data. We calculate the stress field in the composite and then quantitatively evaluate the kinetics of surface charge transfer and Li transport biased by mechanical stresses. We further model the cyclic behavior of the cell. The repetitive deformation of the active particles and the weakening of the interfacial strength cause gradual increase of the interfacial debonding. The mechanical damage impedes electron transfer, incurs more charge heterogeneity, and results in the capacity degradation in batteries over cycles.

25 ENERGY STORAGE↗

Porosity evolution under increasing tension in wire-arc additively manufactured aluminum using in-situ micro-computed tomography and convolutional neural network

Internal defects such as porosities are often formed in additively manufactured metal components. The pores nucleate, grow, and coalesce to form cracks under loads, leading to eventual catastrophic failure. Here, in this paper, the full-field porosity evolution, including pore growth and coalescence in a wire-arc additively manufactured (WAAM) aluminum alloy cylinder under tension is observed with in-situ X-ray micro-computed tomography (μCT). The pore size distribution, density, and tensile stress are calculated from the volumetric images analyzed by a convolutional neural network (CNN) algorithm, which provides rapid analysis of 12,950 slice images from μCT volumetric images at the reference state and 13 tensile strains. The results show the quantitative evolution of the growth and coalescence of macropores under tension. A strong correlation is found between the local pore volume fraction and the true tensile stress when the tensile strain is larger than 5%.

36 MATERIALS SCIENCE↗

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Precipitation Behavior and Hot Deformation Characteristics of Electron Beam Smelted Inconel 740 Superalloy

The Inconel 740 superalloy was prepared by the electron beam smelting (EBS) technology, the precipitation behavior and strengthening mechanism were studied, and the hot deformation characteristics of EBS 740 superalloy were investigated. The results indicate that the EBS 740 superalloy is mainly strengthened by the mechanism of weakly coupled dislocation shearing, and the resulting critical shear stress is calculated to be 234.6 MPa. The deformation parameters show a great influence on the flow behavior of EBS 740 superalloy. The strain rate sensitivity exponent increases with the increasing of deformation temperature, and the strain hardening exponent shows a decreasing trend with the increasing of strain. The activation energy of EBS 740 above 800 °C is measured to be 408.43 kJ/mol, which is higher than the 740H superalloy. A hyperbolic-sine-type relationship can be observed between the peak stress and Zener–Hollomon parameter. Nevertheless, the influence of deformation parameters is found to be considerably different at temperatures below and above 800 °C. The size of dynamic recrystallization (DRX) grains decreases with the increasing of strain rate when the strain rate is lower than 1/s, and reverse law can be found at higher strain rate. As a result, a piecewise function is established between the DRX grain size and hot working parameters.

36 MATERIALS SCIENCE↗

Cold pulse and rotation reversals with turbulence spreading and residual stress

Transport modeling based on inclusion of turbulence spreading and residual stresses shows internal rotation reversals and polarity reversal of cold pulses, with a clear indication of nonlocal transport effects due to fast spreading in the turbulence intensity field. The effects of turbulence spreading and residual stress are calculated from the gradient of the turbulence intensity. In the model presented in this paper, the flux is carried by the turbulence intensity field, which in itself is subject to radial transport effects. The pulse polarity inversion and the rotation profile reversal positions are close to the radial location of the stable/unstable transition. Both effects have no direct explanation within the framework of classical transport modeling, where the fluxes are related directly to the linear growth rates, the turbulence intensity profile is not considered and the corresponding residual stress is absent. Our simulations are in qualitative agreement with measurements from ohmically heated plasmas. Rotation reversal at a finite radius is found in situations not displaying saturated confinement, which we identify as situations where the plasma is nearly everywhere unstable. As an additional and new effect, the model predicts a perturbation of the velocity profile following a cold pulse from the edge. This allows direct experimental confirmation of both the existence of residual stress caused by turbulence intensity profiles and fundamental ideas of transport modeling presented here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian estimations of orientation distribution functions from small-angle scattering enable direct prediction of mechanical stress in anisotropic materials

Properties of soft materials are influenced by their anisotropic structuring under nonequilibrium fields. Although anisotropic structure-property relationships have been extensively explored theoretically, comparison to experiments requires determination of the microstructural orientation probability distribution function (OPDF) of microstructural elements. Small angle scattering (SAS) measurements encode information about the OPDF, but tools to navigate this connection are incomplete. Here, we develop and validate an explicit framework to link arbitrary OPDFs to SAS measurements. Specifically, we propose, validate, and apply a method, maximum a posteriori scattering inference (MAPSI), whereby the OPDF may be obtained from SAS measurements using a Bayesian estimation method. Using this method, we obtain estimates of the full 3D OPDF for two model semidilute fd-virus (rodlike) dispersions at concentrations that are approximately equal to and twice the overlap concentration. From the OPDF, we calculate its second and fourth moments and compare these to predictions for a dilute suspension of rigid rods and to a recent theory for semidilute suspensions. Finally, we use both the theoretical and measured moments to calculate the stress, both for dilute and semidilute suspensions. These predictions are not only compared to each other, but also to measured values of the shear stress, and point to new insights into the behavior of suspensions of highly elongated particles in the transition between dilute and semidilute behavior. We also use this new framework to provide perspective on the connection between scalar parameterizations of scattering and the OPDF that have frequently been used in the past. Here, the new tools developed in this work provide an unprecedented path toward experimental validation of dynamical theories of rodlike colloids and polymers, and for measurement of nonequilibrium structures and stresses of other complex fluids and soft materials with SAS.

36 MATERIALS SCIENCE↗

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗

In situ characterization of residual stress evolution during heat treatment of SiC/SiC ceramic matrix composites using high-energy X-ray diffraction

Volumetric strains were measured in silicon carbide/silicon carbide melt-infiltrated ceramic matrix composites (CMCs) at ambient and high temperatures using high-energy synchrotron X-ray diffraction (XRD). Both silicon and silicon carbide constituents were interrogated utilizing a broad spectrum of diffracting planes that would be largely inaccessible to common laboratory XRD equipment. Residual room-temperature principal strains in the melt-infiltrated silicon phase were found to be approximately 1100 mu epsilon in compression, corresponding to stresses of approximately 300 MPa using simplifying constitutive assumptions. Residual room-temperature principal strains in silicon carbide particles found throughout the matrix were approximately 500 mu epsilon in tension, corresponding to approximately 300 MPa. Residual strains were found to decrease considerably as temperatures increased from ambient temperature to 1250 degrees C. Additionally, residual strains returned to approximately preheat treatment values after cool-down to ambient temperature. Strain measurements in the silicon phase were found to be significantly affected by dissolved boron dopant levels causing contraction of the silicon lattice. This contraction must be accounted for in high-temperature experiments for accurate calculation of stresses in the silicon phase.

36 MATERIALS SCIENCE↗

Residual lattice strain in quartzites as a potential palaeo-piezometer

SUMMARY If a crystal lattice is subjected to a stress, it becomes distorted and no longer represents the ideal crystal symmetry, and if the stress introduces defects such as dislocations, some of this distortion is preserved after the applied stress is removed. In this study, we investigate lattice distortion in quartz at the micron scale with synchrotron X-ray Laue diffraction. From Laue images the local deviatoric strain tensor is derived and corresponding stresses are calculated based on elastic properties. The method is applied to metasedimentary quartzites from the Bergell Alps that were deformed at conditions of greenschist facies metamorphism. The residual palaeostrain is represented in maps of the deviatoric strain tensor components and with deviatoric strain axis pole figures. Data suggest overall shortening perpendicular to the schistosity plane but with considerable asymmetry relative to foliation and lineation, probably attributed to simple shear. Crystallographic pole figures from Laue diffraction agree with neutron diffraction and EBSD measurements and display quartz c-axes girdle distributions with maxima also perpendicular to schistosity. The method shows promise to be used as a palaeo-piezometer to unravel the stress field during tectonic deformation.

58 GEOSCIENCES↗

Enhancements to BISON U-Zr Metallic Fuel X447 Example Problem

As development of a metallic fuel modeling capability in BISON has progressed, the need for an example problem used as a comparison basis was observed. Collaborative work between researchers at Los Alamos National Laboratory (LANL) and Idaho National Laboratory (INL) then proceeded to determine a viable rod to use as the basis and create a BISON input deck utilizing as many metallic fuel models as feasible. The basis chosen was what would be considered a generic rod from subassembly X447, an assembly irradiated in EBR-II towards the end of its operating life, heavily based on reported data for fuel pin DP11. Thus, the approach was adopted to use flow characteristics from subassembly X447 as a basis for the convective heat transfer solution, power history and axial power profiles that are representative of rod DP11 from subassembly X447. The rod simulated is a U-10Zr wt% (U-22.5Zr at%) composition. A 2D-RZ mesh would be used to capture axial thermal hydraulic effects, axial swelling and stress-strain calculations over the full length of the rod. After initial work was invested, a refinement of the various models and input parameters was conducted to ensure consistency between operator-declared conditions, model input requirements and those represented in the example problem. This report serves as a synopsis of the enhancements and refinements to the example problem conducted throughout the 2016 fiscal year.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-Energy Density and Superhard Nitrogen-Rich B-N Compounds

The pressure-induced transformation of diatomic nitrogen into nonmolecular polymeric phases may produce potentially useful high-energy-density materials. Here, we combine first-principles calculations with structure searching to predict a new class of nitrogen-rich boron nitrides with a stoichiometry of B 3 N 5 that are stable or metastable relative to solid N 2 and h -BN at ambient pressure. The most stable phase at ambient pressure has a layered structure ( h - B 3 N 5 ) containing hexagonal B 3 N 3 layers sandwiched with intercalated freely rotating N 2 molecules. At 15 GPa, a three-dimensional C 2 2 2 1 structure with single N-N bonds becomes the most stable. This pressure is much lower than that required for triple-to-single bond transformation in pure solid nitrogen (110 GPa). More importantly, C 2 2 2 1 - B 3 N 5 is metastable, and can be recovered under ambient conditions. Its energy density of ~ 3.44 kJ / g makes it a possible high-energy-density material. In addition, stress-strain calculations estimate a Vicker’s hardness of ~ 4 4 GPa . Structure searching reveals a new clathrate sodalitelike BN structure that is metastable under ambient conditions.

36 MATERIALS SCIENCE↗