Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscoplasticity”

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.

117 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A reinvestigation into Munson’s model for room closure in bedded rock salt

Accurate predictions of room closure are important for hazardous waste repositories in rock salt formations, such as the Waste Isolation Pilot Plant (WIPP). When Munson and co-workers simulated several room closure experiments conducted at the WIPP during the 1980’s and 1990’s, their simulated closure curves closely agreed with the closure measurements. A careful review of their work, however, raised concerns and prompted the reinvestigation in this paper. Here, to begin the reinvestigation, Munson’s legacy Room D closure simulation was reasonably recreated in a current-day finite element code. Next, special care was taken to obtain numerically converged results, re-introduce the anhydrite strata intermittently ignored by Munson, and calibrate the Munson–Dawson (M–D) constitutive model for salt as much as possible from laboratory test measurements. When this new model was used to simulate Room D’s closure, it under-predicted the horizontal and vertical closure rates by 2.34x and 3.10x, respectively, at 5.7 years after room excavation. As a result, the M–D model was extended to capture the newly established creep behavior at low equivalent stresses (<8 MPa) and replace the Tresca with the Hosford equivalent stress. Simulations using the new M–D model over-predicted the horizontal closure rate by 1.15x and under-predicted the vertical closure rate by 1.08x at 5.7 years, averaged over three room closure experiments. Although further improvements could be made, the new model has a stronger scientific foundation than Munson’s legacy model and appears ready for careful engineering use.

58 GEOSCIENCES↗
DOE OSTI2026

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗
DOE OSTI2024

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗
DOE OSTI2022

Study of the interplay between lower-order and higher-order energetic strain-gradient effects in polycrystal plasticity

in this report strain-gradient (SG) plasticity refers to a class of non-local theories in which gradients of plastic slip determine the storage of geometrically necessary dislocations, introducing a length-scale dependence in the mechanical behavior of crystalline materials, which is otherwise lacking in local theories. In this work, we incorporate lower-order (LO) and higher-order energetic (HOE) strain-gradient effects into a crystal plasticity fast Fourier transform (FFT)-based formulation to investigate the interplay of the length scale that each strain-gradient term introduces at the microscale, and the mechanical properties that result at the macroscale. For an applicable range of length scales, we consider two systems: a 1-D two-phase face centered cubic (FCC) laminate and a 3-D FCC polycrystal, and two uniaxial deformation modes: monotonic tension and cyclic tension–compression. We show that increases in the individual LO and HOE length scales increase the hardening rate and strength of the material, respectively. When combined, the strong LO hardening is less pronounced than the effect alone due to the lowering of the gradients due to the HOE microstress. We demonstrate that the LO and HOE hardening manifest as “isotropic” (yield surface expansion) and “kinematic” (yield surface shift) effects, respectively, consistent with their theoretical origins. We show that in cyclic loading, the Bauschinger effect emerges in both local and non-local calculations and link its origins and severity to the behavior in the strain field, slip-system rates, and the HOE microforce.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Decoupling indentation size and strain rate effects during nanoindentation: A case study in tungsten

Materials indented at small scales may simultaneously exhibit indentation size and strain rate effects which complicate the identification of the mechanisms that control deformation and strength. Here, in this work, we explore the possibility that indentation size and rate effects in some materials can be decoupled in a simple way. Nanoindentation tests with various load-time histories were carried out to measure the hardness of a tungsten single crystal over a wide range of indentation depths (~500–3600 nm) and indentation strain rates (~5 ×10 –5 –2 ×10 –1 s –1 ). Under these conditions, this material exhibits significant indentation size and rate effects, but the size effect is, to a good approximation, independent of strain rate. It is shown that this behavior can be understood by the Nix-Gao model for the indentation size effect modified to include the effects of a strain rate dependent friction stress. As a consequence, the size and rate dependencies of the hardness can be expressed as the sum of two independent terms: H($\dot{ε}$ i , h c )=H f ($\dot{ε}$ i ) + (H 0 – H f ) $\sqrt{1+ \frac{h^*}{h_c}}$, where H($\dot{ε}$ i , h c ) is the hardness at given indentation strain rate ($\dot{ε}$ i ) and contact depth (h c ), H f ($\dot{ε}$ i ) is the hardness contributed by the rate dependent friction stress, and (H 0 – H f ) and h* are size and rate independent constants that follow from the Nix-Gao analysis. This formula, together with an expression for the rate dependence of H f , was successfully applied to decouple the indentation size and rate effects observed in tungsten. In addition, the physics underlying the rate independence of indentation size effect is discussed, which provides guidance for application of the proposed approach to other materials.

42 ENGINEERING↗
DOE OSTI2022

Elevated temperature contact creep and friction of nickel-based superalloys using machine learning assisted finite element analysis

Nickel-based superalloys with superior thermochemical, mechanical, and tribological properties are highly utilized for critical components in several high temperature applications such as gas turbines and nuclear reactors. Inconel 617, in particular, is considered as one of the main candidate superalloys for tribo-components in very-high-temperature gas-cooled nuclear reactors. Recent findings indicate that this alloy grows unique surface oxide especially in a high-temperature helium environment with distinctive wear, friction, and contact properties. Here, this study investigates the high temperature contact area evolution and frictional behavior of Inconel 617 using finite element simulation and provides predictive models for the contact and friction performance at different normal loads, dwell times, and temperatures. High temperature helium-aged Inconel 617 top surface properties (up to 600 °C) are utilized along with a single asperity-based deformable elastic-plastic contact model under combined normal and tangential loading. Machine learning is used to assist the finite element results and to predict friction coefficient as well as contact area evolution. While a small difference is observed in the instantaneous friction coefficient (no dwell time) for all temperatures, friction coefficient increases considerably with dwell time. This shows that the effect of contact creep for longer dwell times significantly dominants the effect of high temperature variation in basic mechanical parameters such as modulus and yield strength. It is found that increasing temperature and dwell times lead to the friction coefficient increase, yet the dominance of dwell time effects decreases at higher temperatures and loads. While the analysis is presented for Inconel 617, the methodology is easy to be generalized and can be applied to other HT alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shockwave dissipation by interface-dominated porous structures

The advent of additive manufacturing (AM) has enabled topological control of structures at the micrometer scale, transforming the properties of polymers for a variety of applications. Examples include tailored mechanical responses, acoustic properties, and thermal properties. Porous polymer materials are a class of materials used for shock and blast mitigation, yet they frequently possess a lack of structural order and are largely developed and evaluated via trial-and-error. Here, we demonstrate control of shockwave dissipation through interface-dominated structures prepared by AM using 2-photon polymerization. A fractal structure with voids, or free surfaces, arranged less than 100 μm apart, allows for rarefaction interactions on the timescale of the shockwave loading. Simulations and dynamic x-ray phase contrast imaging experiments show that fractal structures with interfaces assembled within a “critical” volume reduce shockwave stress and wave velocity by over an order of magnitude within the first unit cell.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Viscoelastic Material Characterization and Modeling of Photovoltaic Module Packaging Materials for Direct Finite-Element Method Input

Numerical tools, such as the finite-element method, are increasingly used to design and evaluate the photovoltaic (PV) modules, providing for the reduction of development time and improved performance and reliability. However, high-fidelity material models are necessary to accurately model the complex structural behavior of the involved packaging materials. A common simplification used in recent years is to model the polymer materials (i.e., encapsulant and backsheet) as linear elastic, which will lead to inaccurate results. Therefore, in this work, we present a thorough characterization of the time- and temperature-dependent mechanical response of predominant PV module encapsulant and backsheet materials. Based on this material characterization, we developed and experimentally validated generalized Maxwell models to describe each material's viscoelastic response. In addition, we included measurements of the coefficient of thermal expansion and presented all material models in such a fashion for direct input into commercial finite-element method modeling software.

FEM↗