Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive behavior”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Shear Band Formation in Plane Strain Experiments of Sand

A series of biaxial (plane strain) experiments were conducted on three sands under low (15 kPa) and high (100 kPa) confining pressure conditions to investigate the effects of specimen density, confining pressure, and sand grain size and shape on the constitutive and stability behavior of granular materials. The three sands used in the experiments were fine-, medium-, and coarse-grained uniform silica sands with rounded, subangular, and angular grains, respectively. Specimen deformation was readily monitored and analyzed with the help of a grid pattern imprinted on the latex membrane. The overall stress-strain behavior is strongly dependent on the specimen density, confining pressure, sand grain texture, and the resulting failure mode(s). That became evident in different degrees of softening responses at various axial strains. The relationship between the constitutive behavior and the specimens' modes of instability is presented. The failure in all specimens was characterized by two distinct and opposite shear bands. It was found that the measured dilatancy angles increase as the sand grains' angularities and sizes increase. The measured shear band inclination angles are also presented and compared with classical Coulomb and Roscoe solutions.

Alshibli, Khalid A.↗
NASA NTRS1983

Development of constitutive models for cyclic plasticity and creep behavior of super alloys at high temperature

An uncoupled constitutive model for predicting the transient response of thermal and rate dependent, inelastic material behavior was developed. The uncoupled model assumes that there is a temperature below which the total strain consists essentially of elastic and rate insensitive inelastic strains only. Above this temperature, the rate dependent inelastic strain (creep) dominates. The rate insensitive inelastic strain component is modelled in an incremental form with a yield function, blow rule and hardening law. Revisions to the hardening rule permit the model to predict temperature-dependent kinematic-isotropic hardening behavior, cyclic saturation, asymmetric stress-strain response upon stress reversal, and variable Bauschinger effect. The rate dependent inelastic strain component is modelled using a rate equation in terms of back stress, drag stress and exponent n as functions of temperature and strain. A sequence of hysteresis loops and relaxation tests are utilized to define the rate dependent inelastic strain rate. Evaluation of the model has been performed by comparison with experiments involving various thermal and mechanical load histories on 5086 aluminum alloy, 304 stainless steel and Hastelloy X.

Haisler, W. E.↗
NASA NTRS1988

Unified constitutive model for single crystal deformation behavior with applications

Single crystal materials are being used in gas turbine airfoils and are candidates for other hot section components because of their increased temperature capabilities and resistance to thermal fatigue. Development of a constitutive model which assesses the inelastic behavior of these materials has been studied in 2 NASA programs: Life Prediction and Constitutive Models for Engine Hot Section Anisotropic Materials and Biaxial Constitutive Equation Development for Single Crystals. The model has been fit to a large body of constitutive data for single crystal PWA 1480 material. The model uses a unified approach for computing total inelastic strains (creep plus plasticity) on crystallographic slip systems reproducing observed directional and strain rate effects as a natural consequence of the summed slip system quantities. The model includes several of the effects that have been reported to influence deformation in single crystal materials, such as shear stress, latent hardening, and cross slip. The model is operational in a commercial Finite Element code and is being installed in a Boundary Element Method code.

Walker, K. P.↗
DOE OSTI2026

A physically based mechanical model for Mullins effect in thermoplastic polyurethanes

Despite decades of research, connecting the chemical and physical structure of thermoplastic polyurethanes to their mechanical properties remains highly challenging. Of particular note are their large-deformation and rate-dependent behaviors, which vary greatly with molecular chemistry, including the type and relative content of soft and hard segments. In this work, we develop a physically motivated mechanical theory for predicting the behavior of thermoplastic polyurethanes. The theory incorporates a representation of microstructural evolution during mechanical deformation, which captures the signatures of stress softening over cyclic loading (commonly referred to as the Mullins effect). There are only eight physically motivated fitting parameters, including a direct dependence on the hard segment fraction. The model predicts that increasing the hard segment fraction leads to higher stiffness and greater energy dissipation, in quantitative agreement with published experimental data. Furthermore, we provide a comprehensive analysis of the model and validate its predictions across several independent datasets focused on mechanical characterization. Direct comparisons to experimental data demonstrate its predictive capability on the effect of loading rate, cyclic deformations, and applied tension or compression. Altogether, this work establishes a predictive framework that connects polymer chemistry and microstructure to emergent mechanical behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanical behavior and texture evolution of WE43 magnesium-rare earth alloy in Split-Hopkinson Pressure Bar and Taylor Impact Cylinder Testing

Mechanical behavior and texture evolution of Mg rare-earth alloy WE43 is investigated for strain-rates 10 -3 /s to upwards of 10 5 /s for the two material conditions - as-cast (AC) and T6 age hardened, rolled plate (RT6). The high strain-rate behavior is tested using both Taylor cylinder impact tests (TC) and split Hopkinson pressure bar tests (SHB) and bulk textures are obtained using neutron diffraction. Unlike the quasi-static strained material, AC and RT6 SHB retained high hardening rates throughout the test, even up to 30% true strain. Moreover, the high strain-rate data revealed that the RT6 material has a much higher strength than the AC material, but similar hardening rates despite significantly different initial texture. The flow stress near yield increased up to I 0% for RT6 and up to 30% for AC as the strain-rate increased six orders of magnitude from quasi-static rates 10 -3 /s to 10 3 /s. Neither material exhibited significant plastic anisotropy over the broad range of strain rates, despite the fact that the RT6 material had a moderately strong initial texture. In the TC tests, the geometric cross-sectional changes and texture along the cylinder from the cylindrical sample foot to head are measured and from the neutron diffraction texture analysis, upper-bound estimates of twin volume fraction are obtained as well as dislocation density from analyzing diffraction peak broadening. Recorded geometrical changes along several loading directions show that the material has deformed remarkably homogeneously under impact. Analysis of deformed textures indicates that {101¯2} extension deformation twinning occurred in the RT6 condition over the range of strain rates, with an upper bound estimate of 40% twin volume fraction for approximately 0.10-0.25 true strain. The peak texture components after the impact have their c-axes closely aligned with the impact direction. As a result, these observations are presented and rationalized in the paper.

42 ENGINEERING↗
NASA NTRS1985

Constitutive modeling of superalloy single crystals with verification testing

The goal is the development of constitutive equations to describe the elevated temperature stress-strain behavior of single crystal turbine blade alloys. The program includes both the development of a suitable model and verification of the model through elevated temperature-torsion testing. A constitutive model is derived from postulated constitutive behavior on individual crystallographic slip systems. The behavior of the entire single crystal is then arrived at by summing up the slip on all the operative crystallographic slip systems. This type of formulation has a number of important advantages, including the prediction orientation dependence and the ability to directly represent the constitutive behavior in terms which metallurgists use in describing the micromechanisms. Here, the model is briefly described, followed by the experimental set-up and some experimental findings to date.

Jordan, Eric↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2021

Recurrent neural network-based multiaxial plasticity model with regularization for physics-informed constraints

We report a recurrent neural network (RNN) based model is developed as a surrogate to predict nonlinear plastic response under multiaxial loading. The RNN-based model is trained and tested on stress versus strain curves generated using a numerical solution based on the classical radial return method. Besides simply learning the basic constitutive relationship, a novel approach is taken to enforce certain physical conditions. Specifically, regularization is employed to maintain non-negative plastic power density throughout the loading history thereby ensuring monotonically increasing plastic work and thermodynamic consistency. Enforcing physics in this manner permits coupling of the data-driven RNN approach with physics-based knowledge and laws. This has the effect of reducing the necessary amount of data and ensuring known physical laws are not violated. Since, once trained, the model need not perform the expensive task of solving nonlinear equations, its efficiency is orders of magnitude greater than its numerical counterpart. The RNN-based model has been trained on varied sets of data and the accuracy on test datasets validated. The developed model is general and robust and has widespread application such as in the simulation of metal forming, large scale plasticity, and part life prediction.

42 ENGINEERING↗
DOE OSTI2022

Analytic model of dislocation density evolution in fcc polycrystals accounting for dislocation generation, storage, and dynamic recovery mechanisms

Here, an analytic model of the evolution of dislocation density in fcc polycrystals is described. The evolution equations approximately account for most known dislocation storage, dynamic recovery, and dislocation generation mechanisms in fcc polycrystals. Specifically, the model incorporates network (forest) and grain boundary storage, mobile-network and mobile–mobile annihilation, screw–screw annihilation via athermal and thermal single cross-slip, generation by double cross-slip (Koehler mechanism, including dipole formation), Frank-Read sources, grain boundary nucleation, and mobile–immobile dislocation nucleation due to shock loading. Single cross-slip is assumed to proceed through the Friedel–Escaig (FE) mechanism; the corresponding activation energy is calculated using a modified FE model. The activation energy for double cross-slip is calculated for the first time by extending the FE model. The exact evolution equations are integro-differential equations, and as such are difficult to implement in a code; hence, the evolution equations are simplified by making several approximations. Preliminary results on copper are presented, including comparisons to experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling the non-Schmid crystallographic slip in MAX phases

We present a crystal plasticity constitutive relation for the description of experimentally observed non-Schmid crystallographic slip in a class of ternary carbides and nitrides commonly referred to as MAX phases. In the constitutive relation, we assume that the evolution of the slip system strength in MAX phases has two components – a classical component that depends on the Taylor cumulative shear strain and a non-Schmid component that depends on the stress normal to the slip plane. The non-Schmid crystal plasticity constitutive relation is then used to carry out finite element simulations of micropillar compression of single crystals of two MAX phases, Ti 2 AlC and Ti 3 AlC 2 . The finite element simulations not only quantitatively predict the stress – strain response of a wide range of crystallographic orientations of the micropillars but also rationalize the non-uniform deformation and the deformed shape of the micropillars observed in the experiments for the two materials. As a result, parametric studies are also carried out to quantify the role of the non-Schmid effect and understand the effects of key experimental parameters on the stress – strain response of the micropillars of the two MAX phases.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deep material network via a quilting strategy: visualization for explainability and recursive training for improved accuracy

Recent developments integrating micromechanics and neural networks offer promising paths for rapid predictions of the response of heterogeneous materials with similar accuracy as direct numerical simulations. The deep material network is one such approaches, featuring a multi-layer network and micromechanics building blocks trained on anisotropic linear elastic properties. Once trained, the network acts as a reduced-order model, which can extrapolate the material’s behavior to more general constitutive laws, including nonlinear behaviors, without the need to be retrained. However, current training methods initialize network parameters randomly, incurring inevitable training and calibration errors. Here, we introduce a way to visualize the network parameters as an analogous unit cell and use this visualization to “quilt” patches of shallower networks to initialize deeper networks for a recursive training strategy. The result is an improvement in the accuracy and calibration performance of the network and an intuitive visual representation of the network for better explainability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗
NASA NTRS1983

Strainrange partitioning: A total strain range version

Procedures are presented for expressing the Strainrange Partitioning (SRP) method for creep fatigue life prediction in terms of total strain range. Inelastic and elastic strain-range - life relations are summed to give total strain-range - life relations. The life components due to inelastic strains are dealt with using conventional SRP procedures while the life components due to elastic strains are expressed as families of time-dependent terms for each type of SRP cycle. Cyclic constitutive material behavior plays an important role in establishing the elastic strain-range - life relations as well as the partitioning of the inelastic strains. To apply the approach, however, it is not necessary to have to determine the magnitude of the inelastic strain range. The total strain SRP approach is evaluated and verified using two nickel base superalloys, AF2-1DA and Rene 95. Excellent agreement is demonstrated between observed and predicted cyclic lifetimes with 70 to 80 percent of the predicted lives falling within factors of two of the observed lives. The total strain-range SRP approach should be of considerable practical value to designers who are faced with creep-fatigue problems for which the inelastic strains cannot be calculated with sufficient accuracy to make reliable life predictions by the conventional inelastic strain range SRP approach.

Halford, G. R.↗
NASA NTRS1985

Finite element analysis of notch behavior using a state variable constitutive equation

The state variable constitutive equation of Bodner and Partom was used to calculate the load-strain response of Inconel 718 at 649 C in the root of a notch. The constitutive equation was used with the Bodner-Partom evolution equation and with a second evolution equation that was derived from a potential function of the stress and state variable. Data used in determining constants for the constitutive models was from one-dimensional smooth bar tests. The response was calculated for a plane stress condition at the root of the notch with a finite element code using constant strain triangular elements. Results from both evolution equations compared favorably with the observed experimental response. The accuracy and efficiency of the finite element calculations also compared favorably to existing methods.

Dame, L. T.↗
NASA NTRS1983

Strainrange partitioning - A total strain range version

Procedures are presented for expressing the Strainrange Partitioning (SRP) method for creep fatigue life prediction in terms of total strain range. Inelastic and elastic strain-range - life relations are summed to give total strain-range - life relations. The life components due to inelastic strains are dealt with using conventional SRP procedures while the life components due to elastic strains are expressed as families of time-dependent terms for each type of SRP cycle. Cyclic constitutive material behavior plays an important role in establishing the elastic strain-range life relations as well as the partitioning of the inelastic strains. To apply the approach, however, it is not necessary to have to determine the magnitude of the inelastic strain range. The total strain SRP approach is evaluated and verified using two nickel base superalloys, AF2-1DA and Rene 95. Excellent agreement is demonstrated between observed and predicted cyclic lifetimes with 70 to 80 percent of the predicted lives falling within factors of two of the observed lives. The total strain-range SRP approach should be of considerable practical value to designers who are faced with creep-fatigue problems for which the inelastic strains cannot be calculated with sufficient accuracy to make reliable life predictions by the conventional inelastic strain range SRP approach.

Halford, G. R.↗