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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Data-Driven Surrogate Modeling with Microstructure-Sensitivity of Viscoplastic Creep in Grade 91 Steel

Abstract To support the development of advanced steel alloys tailored to withstand extreme conditions, it is imperative to account for the mechanical performance of components, while considering the influence of local microstructure on the macroscopic response. To this end, this study focuses on the development of microstructure-sensitive constitutive models for the mechanical response of Grade 91 steel exposed to extreme thermo-mechanical environments. Polynomial chaos expansion (PCE) surrogates are used to emulate high-fidelity polycrystal simulations of the viscoplastic response of Grade 91 steel as a function of the microstructure fingerprint (e.g., dislocations and precipitates). To cover a wide temperature–stress domain, two separate PCE surrogates—one that captures softening and the other that captures hardening behavior—are combined using another (sparse) Gaussian process regression model. The resulting constitutive creep surrogate model is integrated within the MOOSE finite element framework to simulate the intricate effects of microstructure, in particular MX-phase precipitates, on a component with a graded microstructure. Surrogate sensitivity analysis is applied to quantify the relevant impact of spatially varying microstructure on the creep response in a test-case involving a Grade 91 alloy with a prototypical weld.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Engineering-scale Modeling of High-Temperature Creep and Creep Crack Growth in Alloy 316H

This document demonstrates completion of the goals described in the technical narrative of the Department of Energy’s Industry Funding Opportunity Announcement (iFOA) project entitled “Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing,” which relates to the development and demonstration of capabilities for conducting engineering-scale simulations of Alloy 316H components under high-temperature conditions. This work encompassed two major aspects: integrating and testing constitutive models for the creep response of 316H at high temperatures, and developing tools for modeling creep crack growth (CCG) in 316H. Two classes of constitutive models were used in this effort: phenomenological models based on behavior observed at the engineering scale, and reduced-order models (ROMs) that represent the nonlinear response of mesoscale models that capture the sensitivity to material microstructure and processing. Likewise, the approaches employed for CCG modeling considered both simplified engineering approaches and detailed simulations of creep and damage ahead of the crack tip. These developments, which were performed by utilizing the Grizzly code as well as the open-source libraries it depends on, strengthen Grizzly’s ability to support licensing and safety analyses of high-temperature reactor components.

316H↗
DOE OSTI2022

Bone: An Outstanding Composite Material

Bone is an outstanding, well-designed composite. It is constituted by a multi-level structure wherein its properties and behavior are dependent on its composition and structural organization at different length scales. The combination of unique mechanical properties with adaptive and self-healing abilities makes bone an innovative model for the future design of synthetic biomimetic composites with improved performance in bone repair and regeneration. However, the relation between structure and properties in bone is very complex. In this review article, we intend to describe the hierarchical organization of bone on progressively greater scales and present the basic concepts that are fundamental to understanding the arrangement-based mechanical properties at each length scale and their influence on bone’s overall structural behavior. The need for a better understanding of bone’s intricate composite structure is also highlighted.

Rosa, Natacha (ORCID:0000000199008851)↗
NASA NTRS1978

Interpolation and Extrapolation of Creep Rupture Data by the Minimum Commitment Method. Part 3: Analysis of Multiheats

The Minimum Commitment Method was applied to two sets of data for which multiple heat information was available. For one alloy, a 304 stainless steel studied in Japan, data on nine well characterized heats were used, while for a proprietary low alloy carbon steel studied in the United Kingdom data were available on seven heats - in many cases to very long rupture times. For this preliminary study no instability factors were used. It was discovered that heat-to-heat variations would be accounted for by introducing heat identifiers in the form A + B log sigma where sigma is the stress and the constants A and B depend only on the heat. With these identifiers all the data could be collapsed onto a single master curve, even though there was considerable scatter among heats. Using these identifiers together with the average behavior of all heats made possible the determination of an accurate constitutive equation for each individual heat. Two basic approaches are discussed for applying the results of the analysis.

Manson, S. S.↗
NASA NTRS1986

A viscoplastic constitutive theory for metal matrix composites at high temperature

A viscoplastic constitutive theory is presented for representing the high-temperature deformation behavior of metal matrix composites. The point of view taken is a continuum one where the composite is considered a material in its own right, with its own properties that can be determined for the composite as a whole. It is assumed that a single preferential (fiber) direction is identifiable at each material point (continuum element) admitting the idealization of local transverse isotropy. A key ingredient in this work is the specification of an experimental program for the complete determination of the material functions and parameters for characterizing a particular metal matrix composite. The parameters relating to the strength of anisotropy can be determined through tension/torsion tests on longitudinally and circumferentially reinforced thin-walled tubes. Fundamental aspects of the theory are explored through a geometric interpretation of some basic features analogous to those of the classical theory of plasticity.

Robinson, D. N.↗
NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗
NASA NTRS2007

GVIPS Models and Software

Two reports discuss, respectively, (1) the generalized viscoplasticity with potential structure (GVIPS) class of mathematical models and (2) the Constitutive Material Parameter Estimator (COMPARE) computer program. GVIPS models are constructed within a thermodynamics- and potential-based theoretical framework, wherein one uses internal state variables and derives constitutive equations for both the reversible (elastic) and the irreversible (viscoplastic) behaviors of materials. Because of the underlying potential structure, GVIPS models not only capture a variety of material behaviors but also are very computationally efficient. COMPARE comprises (1) an analysis core and (2) a C++-language subprogram that implements a Windows-based graphical user interface (GUI) for controlling the core. The GUI relieves the user of the sometimes tedious task of preparing data for the analysis core, freeing the user to concentrate on the task of fitting experimental data and ultimately obtaining a set of material parameters. The analysis core consists of three modules: one for GVIPS material models, an analysis module containing a specialized finite-element solution algorithm, and an optimization module. COMPARE solves the problem of finding GVIPS material parameters in the manner of a design-optimization problem in which the parameters are the design variables.

Arnold, Steven M.↗
NASA NTRS2010

Linear Stability of Binary Alloy Solidification for Unsteady Growth Rates

An extension of the Mullins and Sekerka (MS) linear stability analysis to the unsteady growth rate case is considered for dilute binary alloys. In particular, the stability of the planar interface during the initial solidification transient is studied in detail numerically. The rapid solidification case, when the system is traversing through the unstable region defined by the MS criterion, has also been treated. It has been observed that the onset of instability is quite accurately defined by the "quasi-stationary MS criterion", when the growth rate and other process parameters are taken as constants at a particular time of the growth process. A singular behavior of the governing equations for the perturbed quantities at the constitutional supercooling demarcation line has been observed. However, when the solidification process, during its transient, crosses this demarcation line, a planar interface is stable according to the linear analysis performed.

Mazuruk, K.↗
NASA NTRS2004

Nonlinear Thermoelastic Model for SMAs and SMA Hybrid Composites

A constitutive mathematical model has been developed that predicts the nonlinear thermomechanical behaviors of shape-memory-alloys (SMAs) and of shape-memory-alloy hybrid composite (SMAHC) structures, which are composite-material structures that contain embedded SMA actuators. SMAHC structures have been investigated for their potential utility in a variety of applications in which there are requirements for static or dynamic control of the shapes of structures, control of the thermoelastic responses of structures, or control of noise and vibrations. The present model overcomes deficiencies of prior, overly simplistic or qualitative models that have proven ineffective or intractable for engineering of SMAHC structures. The model is sophisticated enough to capture the essential features of the mechanics of SMAHC structures yet simple enough to accommodate input from fundamental engineering measurements and is in a form that is amenable to implementation in general-purpose structural analysis environments.

Turner, Travis L.↗
NASA NTRS2017

Grain Size-Dependent Crystal Plasticity Constitutive Model for Polycrystal Materials

A new method to introduce grain size-dependence in classical crystal plasticity constitutive model is developed by considering the resistance to dislocation motion in the grain boundary influence region as equivalent to that of a work hardening. A general framework for the size-dependent constitutive model is derived by implementing this method on a core and mantle model. The work hardening, equivalent to grain boundary effect, is realized by introducing a fictitious, pre-existing, plastic strain grain boundary influence region (mantle) following the principles of classical crystal plasticity. This fictitious plastic strain, in effect, increases the yield strength and decreases the initial hardening coefficient of the grain. With the thickness of grain boundary influence region and the distribution of introduced plastic strain remaining the same, the grain boundary effect increased as the grain size becomes smaller. A simplified model that considers the grain boundary effect on a grain average sense is also developed under this general framework. Implementation of this general framework to the specific cases of crystal plasticity constitutive models is demonstrated by considering the case of power law flow rule and hyperbolic-secant hardening rule. Finally, the grain size-dependent constitutive model is validated by comparing the predicted stress-strain behavior of polycrystal samples with different average grain sizes under uniaxial loading with the experimental results.

Grain Size↗
DOE OSTI2021

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗
NASA NTRS1988

Finite element (MARC) solution technologies for viscoplastic analyses

A need for development of realistic constitutive models for structural components operating at high temperatures, accompanied by appropriate solution technologies for stress/life analyses of these components is studied. Viscoplastic models provide a better description of inelastic behavior of materials, but their mathematical structure is very complex. The highly nonlinear and stiff nature of the constitutive equations makes analytical solutions difficult. Therefore, suitable solution, finite element or other numerical, technologies must be developed to make these models adaptable for better and rational designs of components. NASA-Lewis has developed several solution technologies and successfully applied them to the solution of a number of uniaxial and multiaxial problems. Some of these solution technologies are described along with the models and representative results. The solution technologies developed and presented encompass a wide range of models, such as, isotropic, anisotropic, metal matrix composites, and single crystal models.

Arya, V. K.↗
DOE OSTI2026

Molecular Dynamics Investigation of Oil Wetting on Synthetic Polymer Substrates

Understanding the interaction of polymer surfaces with nonpolar, low surface tension liquids, or whether a substrate is oleophobic versus oleophilic, is critical for applications ranging from antifingerprint coatings to oil–water separation membranes and oil spill remediation. Despite its technological importance, the molecular mechanisms underpinning polymer–oil wetting are not well characterized. Here, we employ molecular dynamics simulations to investigate the behavior of n-hexadecane in contact with chemically distinct polymer surfaces, spanning eight constitutional unit chemistries as well as amorphous and crystalline morphologies. This permits a critical examination of both thermodynamic and dynamic descriptors of oleophobicity, including oil contact angle, dewetting free energy, interfacial diffusivity, and a proposed “ghost probe energy,” which does not require explicit simulation of oil–polymer interactions. We find that the oil contact angle does not reliably distinguish oleophobic behavior across polymer chemistries, whereas other metrics provide clearer and more consistent differentiation. Analysis of the results reveals that polymer–oil wetting behavior is primarily governed by interfacial van der Waals interactions and modulated by surface flexibility and morphology. Collectively, this work establishes a computational framework for characterizing oil wetting and provides additional insight into what molecular-level factors dictate trends in oleophobicity.

lipids↗
NASA NTRS1992

Thermo-elastoviscoplastic snapthrough behavior of cylindrical panels

The thermo-elastoviscoplastic snapthrough behavior of simply supported cylindrical panels is investigated. The analysis is based on nonlinear kinematic relations and nonlinear rate-dependent unified constitutive equations which include both Bodner-Partom's and Walker's material models. A finite element approach is employed to predict the inelastic buckling behavior. Numerical examples are given to demonstrate the effects of several parameters which include the temperature, thickness and flatness of the panel. Comparisons of buckling responses between Bodner-Partom's model and Walker's model are given. The creep buckling behavior, as an example of time-dependent inelastic deformation, is also presented.

Song, Y.↗
NASA NTRS1996

First Cycle Heterogeneous Deformation Behavior and Cyclic Shakedown Phenomena of Nitinol Near A(sub f) Temperatures

Experimental observations on the cyclic behavior of a NiTi alloy (Nitinol) at temperatures in the neighborhood of the A(sub f) (austenite finish) temperature are presented. The strongly heterogeneous nature of the deformation behavior of this material at temperatures within this regime during the first cycle is examined with emphasis placed on the difficulties that the existence of such phenomena pose on the formulation of realistic constitutive relations. It is further demonstrated that this heterogeneity of deformation persists on subsequent cycles with the result that the hysteretic cyclic behavior of these alloys can exhibit a point to point variation in an otherwise uniform geometry. The experimental observations on the deformation behavior of this alloy show that it is strongly dependent on temperature and prior deformation history of the sample, thus resulting in an almost intractable problem with respect to capturing an adequate constitutive description from either experiment or modeling.

Jones, H. N.↗
DOE OSTI2021

Elucidating the structural properties and reversible regional texture effect of GdB 6 under high pressure

As hard and refractory materials with high chemical resistance and mechanical strength, lanthanide hexaborides (LnB 6 ) have attracted much attention. Among the family of LnB 6 , gadolinium hexaboride (GdB 6 ) occupies a special position due to the half-filled 4f shell of Gd. Here, using in situ synchrotron radiation angle-dispersive X-ray diffraction in a diamond anvil cell at room temperature, the structural stability and compression behavior of GdB 6 are investigated extensively, GdB 6 is observed to be structurally stable up to 73 GPa, and the bulk modulus of 177 GPa is obtained under hydrostatic compression. In this paper, an interesting observation of pressure-induced spotty diffraction rings of GdB6 is reported, its formation mechanism can be well described by the reversible regional texture effect, which highlights the rearrangement of crystal grains under high pressure. The rearrangement mechanism has been well explained by investigating the pressure dependence of full width at half maximum, macro-differential stress and grain size. Collective grain rotation behavior motivated by stress difference is critical for the rearrangement process, the strong isotropy and strong stability of GdB 6 structure also provide necessary conditions for high-pressure grain rotation behavior. These results will help to promote the understanding of high-pressure structural properties of GdB 6 , and provide novel insights on the high-pressure grain behavior in hard materials with strong isotropy and three-dimensional skeleton constituted of strong covalent bonds.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Model calibration of the liquid mercury spallation target using evolutionary neural networks and sparse polynomial expansions

The mercury constitutive model predicting the strain and stress in the target vessel plays a central role in improving the lifetime prediction and future target designs of the mercury targets at the Spallation Neutron Source. We leverage the experiment strain data collected over multiple years to improve the mercury constitutive model through a combination of large scale simulations of the target behavior and the use of machine learning tools for parameter estimation. We present two interdisciplinary approaches for surrogate-based model calibration of expensive simulations using evolutionary neural networks and sparse polynomial expansions. The newly calibrated simulations achieve 7% average improvement on the prediction accuracy and 8% reduction in mean absolute error compared to previously reported reference parameters, with some individual sensors experiencing up to 30% improvement. The calibrated simulations can aid in fatigue analysis to estimate the mercury target lifetime, which reduces abrupt failure and saves tremendous amount of costs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗