Next Article in Journal
The Microstructure Evolution and Mechanical Properties of Rotary Friction Welded Duplex Stainless Steel Pipe
Next Article in Special Issue
Effects of La Addition on Microstructure Evolution and Thermal Stability of Cu-2.35Ni-0.59Si Sheet
Previous Article in Journal
Encapsulation of Silicon Nano Powders via Electrospinning as Lithium Ion Battery Anode Materials
Previous Article in Special Issue
Energy Storage and Dissipation in Consecutive Tensile Load-Unload Cycles of Gum Metal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Review: Inelastic Constitutive Modeling: Polycrystalline Materials †

1
Department of Civil and Environmental Engineering, Washkewicz College of Engineering, Cleveland State University, Cleveland, OH 44115, USA
2
Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
*
Author to whom correspondence should be addressed.
This manuscript has been authored by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan, accessed on 24 April 2023).
Materials 2023, 16(9), 3564; https://doi.org/10.3390/ma16093564
Submission received: 25 February 2023 / Revised: 13 April 2023 / Accepted: 25 April 2023 / Published: 6 May 2023

Abstract

:
This article provides a literature review that details the development of inelastic constitutive modeling as it relates to polycrystalline materials. This review distinguishes between inelastic constitutive models that account for nonlinear behavior at the microstructural level, time-independent classic plasticity models, and time-dependent unified models. Particular emphasis is placed on understanding the underlying theoretical framework for unified viscoplasticity models where creep and classical plasticity behavior are considered the result of applied boundary conditions instead of separable rates representing distinct physical mechanisms. This article establishes a clear understanding of the advantages of the unified approach to improve material modeling. This review also discusses recent topics in constitutive modeling that offer new techniques that bridge the gap between the microstructure and the continuum.

Graphical Abstract

1. Introduction

Nonlinear hereditary inelastic deformation behavior can occur in many materials utilized at elevated service temperatures. This behavior can include creep, rate sensitivity, and plasticity. Accurate assessments of nonlinear stress and deformation behavior are important in predicting the operational life and overall performance of critical engineering systems. Inelastic constitutive models have been developed and deployed to meet these assessment needs. The models must predict nonlinear behavior under complex thermomechanical load paths. This includes capturing phenomena such as Bauschinger’s effect, cyclic softening and hardening, stress relaxation, and ratcheting when present in high-temperature applications.
In this article, focus is given to theoretical concepts that historically supported the development of unified inelastic constitutive models. The field is rich, and because of this, there will be some oversights here, and apologies are offered a priori. The discussion begins by mentioning the work of Andrade [1,2] with soft metals (e.g., lead). Andrade modeled transient creep strain by assuming it is proportional to the cubic root of time. Expressions for different stages of creep (primary, steady state, but not tertiary) have been offered by Norton [3] and Bailey [4]. The generalization of Norton’s power law to multiaxial states of stress using a potential function is highlighted later in the viscoplasticity section of this paper. Odqvist [5] published a multiaxial formulation for steady-state (secondary) creep prior to the modern-day work on unified viscoplasticity models, and he later extended the model to include primary creep [6]. However, as indicated later, experimental evidence supports the use of unified models in modeling the first two stages of creep as opposed to a piecemeal approach. In addition, the interaction between creep and plasticity must be addressed in a unified fashion. Nonetheless, some attention has been given here to certain material science models that explain nonlinear behavior at the microstructural level. An understanding of material microstructure is always necessary to develop multiaxial continuum-level constitutive models. Reaction-rate and power-law formulations for microstructural constitutive relationships are discussed. The work in developing rheological viscoelastic constitutive models by Schapery [7] and Findley et al. [8] is briefly presented since they point the way to the endochronic viscoplastic models that appear later. The objective of this article is to highlight the genesis of unified time-dependent viscoplastic models.
Early work by Besseling [9] suggested separate strain rates for plasticity and creep that were based on partitioning a total strain rate. Specifically, he proposed a one-dimensional model where the total strain rate could be partitioned as follows:
d ε T d t = d ε E d t + d ε C d t + d ε P d t
where
εTTotal strain
εEElastic strain
εCCreep strain
εPPlastic strain
The inelastic components of this expression are as follows:
d ε I d t =   Inelastic   Strain   Rate   = d ε C d t + d ε P d t = ε ˙ I ε C , ε P , t
where the two rates representing creep and plasticity are identified as separate and distinct deformation mechanisms.
The equations above are associated with microstructural creep mechanisms that are thermally activated. For materials utilized in fabricating nuclear reactor components that function as moderators, e.g., graphite, additional creep mechanisms arise that are activated by irradiation. As Onimus et al. [10] point out, Equation (2) should be modified as follows:
d ε I d t = I n e l a s t i c   S t r a i n   R a t e = d ε P d t + d ε thermal - creep   d t + d ε irradiation - creep   d t + d ε swelling / growth   d t = ε ˙ I ε P , ε thermal - creep   , ε irradiation - creep   , ε swelling / growh   , t
Expressions for creep strain rates are identified as εthermal-creep, εirradiation-creep, and εswelling/growth in publications cited above or based on models established through material science concepts. The forms for creep strain rates based on material science are dependent on various facets of a material microstructure. For example, dislocation mechanics include models for dislocation glide, dislocation creep, and diffusional flow. These mechanisms and the attending creep strain models are discussed in the following section.
For the plasticity component of Equation (2), incremental theories are deployed. Classic plasticity theory poses a constitutive relationship that is independent of the rate of loading. However, plasticity is path-dependent, and load sequences are often represented as a process evolving in time. Hence, classical plasticity theories are posed on an incremental basis, and Equation (2) would be better stated in a differential format:
d ε I = d ε P + d ε thermal - creep   + d ε irradiation - creep + d ε swelling / growth
Equation (4) can be easily integrated with respect to time provided that the load path considered can be expressed as increments of time. Details associated with incremental plasticity models are presented in a later section.
The modeling of plasticity and creep can be accomplished separately, or represented as a single rate, i.e., as a unified model. Nonunified models allow creep strain to be rate-dependent, but plastic strains on the other hand are rate-independent. Alternatively, in unified constitutive models, creep and plastic strains are represented by a single inelastic strain measure that is considered inherently rate-dependent. Nonunified models lack the ability to predict interactions between creep and plasticity, e.g., ratcheting (cyclic creep) and stress relaxation. Experimental evidence such as accumulated creep strain affecting the size of the threshold (yield) surface and limited creep strain recovery during the unloaded periods accompanied by significant reductions in the resistance to creep upon subsequent reloading have been observed that cannot be captured by nonunified models. Jaske et al. [11] as well as Pugh and Robinson [12] have pointed out that the behavior of metallic alloys under cyclic load at elevated service temperatures is indicative of an interaction going on between the creep and plasticity components of Equation (2). Many authors, including Robinson [13,14,15,16,17,18,19,20,21] have pointed out that separating creep and plasticity into distinct strain rates as is performed in nonunified models leads to an inability to capture the interactive behavior exhibited experimentally. Phenomenologically based unified creep–plasticity interaction models were developed from theoretical frameworks based on a single inelastic rate ε ˙ I . The unified manner of representing a constitutive relationship enables a comparison of various aspects of the model to nonlinear behaviors. The details of a number of these unified models are presented later.
The dichotomy between classical mechanisms (plasticity and material science creep) associated with Equation (2) and unified modeling approaches is presented in Figure 1. In that overview figure, inelastic constitutive models are first segregated into three sub-categories: models that explain nonlinear behavior at the microstructural level, time independent classic plasticity models, and time-dependent unified models. The models where the time-dependent inelastic strain is derived from a single potential or a single integral function are classified as unified models.
Plasticity, as well as thermally activated material science creep models, and creep models based on mechanism activated by irradiation are discussed separately in the following two sections for historical completeness. The remainder of the article primarily focuses on unified viscoplasticity models where creep and classical plasticity behavior are considered the result of applied boundary conditions instead of separable rates representing distinct physical mechanisms. The existing models for a unified strain rate (e.g., Robinson et al. [13,14,15,16,17,18,19,20,21], Chaboche et al. [22,23,24], Bodner et al. [25,26,27,28,29,30,31,32,33,34,35], Miller et al. [36,37,38,39,40,41,42,43,44,45,46,47,48] and Walker et al. [49,50,51]) are based on constructs involving internal state variables that capture several inelastic hardening behaviors. Chan et al. [52], as well as Allen and Harris [53], provide overviews of the viscoplastic models just mentioned.
There are micromechanics models for polymer composite material systems such as the models proposed by Chamis [54]. Here, each component of the microstructure is treated as a continuum. These models are discussed elsewhere, and the rest of the article focuses on models that capture the macro-level response of materials exhibiting inelastic behavior.

2. Early Efforts in Modeling Creep Behavior

The Norton–Bailey [3,4] creep model mentioned earlier is a power law formulation that is not associated with a specific microstructural mechanism. In this effort, the Norton–Bailey [3,4] model is categorized as a material science model and is included with other microstructural creep strain models since the model is deployed both at the microstructural level and the continuum level (although the model is one-dimensional). The model was derived predicated on the assumption that creep strain is functionally dependent on stress, time, and temperature. The functional form for the Norton–Bailey [3,4] creep model is based on the following separable formulation
ε C = F 1 σ F 2 t F 3 T
Several forms for F1 can be utilized, including Norton’s [3] model
F 1 σ = A σ m
an exponential format
F 1 σ = C e x p [ σ σ 0 ]
and a hyperbolic sine format
F 1 σ = s i n h m [ σ σ 0 ]
The hyperbolic function is convenient since at the limits for low-stress values, the function approaches values given by the power law, and at high-stress values, the hyperbolic takes on values nearly equal to the exponential function.
The time function F2(t) is also formulated as a power law in time
F 2 t = a t n
Equation (9) is attributable to Bailey [4]. Finally, the temperature function is stipulated as
F 3 T = F 3 [ exp ( Q R T ) ]
In this expression,
Q = activation energy
R = Boltzmann’s constant
Functions F2 and F3 are combined such that
F 4 t , T = α [ t exp ( Q R T ) ]
and
ε C = A σ m [ t exp ( Q R T ) ] n
Under isothermal conditions (T = constant),
ε C = A σ m t n
Under the assumption that the length of exposure time to temperature and stress level are the major factors affecting the creep strain rate,
ε C = F 1 σ F 2 t = A σ m α t n = α A σ m t n
Taking a time derivative of Equation (14) leads to a rate formulation of the Norton creep law
d ε C d t = A ~ σ m t n ~
Note that the constants
  A ~ = α A n
and
n ~ = n 1
are assumed to be temperature-dependent.

3. Creep Strain Models Based on Microstructural Mechanisms

When considering continuum-level deformation for polycrystalline materials, it is important to be cognizant of the microstructural mechanisms driving the deformation processes. There are two primary mechanisms that drive creep strain rates. The first mechanism is diffusional creep (Coble creep [55], as well as Nabarro [56] and Herring [57] creep) where vacancies are driven through the crystal lattice of the material. The second process involves the translation of dislocations through the atomic matrix of the material (see Weertman [58]). For either mechanism, the overall stiffness of the material decreases with time. Stress–strain curves that are initially linear in the elastic range bend over as total strains accumulate. The decrease in stress represents the hardening of the material when subjected to monotonic loads. The decrease in strain rate to a steady-state value in a material subject to a constant load represents two material mechanisms competing and eventually balance a hardening response with a recovery response.
Figure 2 represents a generic deformation mechanism map. In this figure, dominant microstructural mechanisms are mapped to various load and temperature regimes. At low homologous temperatures, i.e., <(0.5) T/Tm, where T is the service temperature and Tm is the material melting temperature, deformations are primarily elastic until high-stress regimes are attained. At moderate temperatures and high stress, dislocations tend to glide on slip planes. Contrasted with other diffusional creep mechanisms, Coble creep [55] is similar to Nabarro–Herring creep [56,57] in that both are dominant at lower stress levels and higher temperatures. Dislocation glide mechanisms occur at higher strength levels at all temperatures. Coble creep [55] occurs through the diffusion of atoms in a polycrystalline material along grain boundaries. At elevated temperatures, dislocation mechanisms are aided by vacancy diffusion. Vacancy diffusion also causes deformations without any dislocation movement involved.
Models exist that attempt to bridge the gap between microstructural-level modeling and continuum-level modeling. These include the rate-dependent deformation models of Ponter and Leckie [59]. Here, a simple constitutive framework incorporates hardening and recovery mechanisms associated with dislocation motion. Similarly, Rice [60] models dislocation motion using the normality structure which has served many times as a fundamental principle in continuum-level macroscopic constitutive laws. In Rice’s [60] model the local microstructural rearrangements proceed at a rate governed by a thermodynamic force.

4. Irradiation-Induced Creep

Another form of plastic deformation is irradiation-induced creep, which occurs in materials exposed to radiation under applied stress. Irradiation-induced creep occurs at temperatures well below 50% of the melting temperature of a material because the collisions of the radiation with the lattice atoms produce the defects (self-interstitials and lattice vacancies) that allow creep to occur rather than relying on thermal lattice vibrations to produce the creep. The mechanisms of diffusional creep (i.e., Coble creep [55] and Nabarro–Herring creep [56,57]) are not enhanced by irradiation [61,62]. Instead, the primary mechanisms of irradiation creep are the motion of dislocations, i.e., either enhanced climb or climb and glide of dislocations, and the preferential formation of loops. Stress-induced preferential absorption (SIPA) occurs when the climb of dislocations is accelerated by preferential absorption generally of interstitials. Stress-induced climb and glide (SICG) occurs when irradiation defects allow dislocations to climb over pinning points and the subsequent glide is the source of the creep strain. When the applied stress affects the climb velocity (i.e., SIPA with SICG), the resulting behavior is referred to as the preferred absorption glide (PAG). Stress-induced preferential nucleation (SIPN) occurs when dislocation loops nucleate in the crystal structure with a preferred orientation.
Creep of alpha-uranium was observed in the 1950s [63]. The mechanism proposed was that internal stresses due to anisotropic growth of the crystals, when combined with externally applied stress, can overcome the yield strength, thereby causing plastic deformation. However, Hesketh [64] instead proposed a mechanism that depends on the anisotropic growth of crystals, the SIPN of dislocation loops, and the climb of the SIPN loops. In austenitic stainless steels and nickel-based alloys, primary creep is mainly due to dislocation slip (SICG) that is reduced as the material hardens (increased pinning point density as fluence increases), while steady-state creep is thought to be primarily due to SICG and SIPA, but the mechanism can change as experimental parameters change [62]. Ferritic–Martensitic steels undergo creep, but SIPA and PAG are not sufficient to describe the total creep, so it is accepted that SIPN is also an active mechanism [62]. In zirconium alloys, irradiation creep is thought to be due to both thermal dislocation slip with irradiation-enhanced climb and SIPA. Cubic silicon carbide (3C-SiC) irradiation creep has a stress exponent near one, indicating a SIPA mechanism, but microscopic studies have also observed an increase in the number of loops in the {111} planes, indicating some influence of SIPN [65].
Irradiation creep in graphite has been studied quite heavily, yet there is still only limited agreement on the microstructural mechanisms. Primary creep in graphite is saturated around one elastic strain unit, defined as applied stress divided by Young’s modulus and is primarily thought to be driven by dislocation bowing [66,67,68] and interactions with pinning points. The secondary creep (i.e., steady state) mechanism was originally proposed to be driven by anisotropic grain growth, but that was shown to not be the case when boron additions changed the rate of grain growth but did not change the creep rate [69]. The long-accepted mechanism is the pinning–unpinning mechanism proposed by Kelly and Foreman [70]. Strain in this mechanism is driven by dislocation glide, and the rate is controlled by the creation and destruction of 2-4 atom clusters that act as the dislocation pinning points. With this being a glide-based mechanism (SICG), it should have a stress dependence of two. Newer published irradiation creep results [66,71,72,73] suggest a SIPA mechanism. Results from Kennedy in the 1960s [66] showed that changing stress during irradiation did not affect the primary creep strain, which provides additional results that disagree with the pinning–unpinning mechanism and instead support a SIPA model.

5. Time-Independent Models Based on Continuum Principles: Plasticity

Historically, the classic plasticity models with flow theories and hardening laws have been proposed to model material time-independent inelastic behavior, i.e., plasticity. The theoretical basis for yielding under complex states of stress at a continuum level stretches back to Coulomb’s [74] effort in defining failure using his nascent concepts in soil mechanics. Yielding in the classic plasticity sense had its origins in the nineteenth century when Tresca [75] conducted a series of punching and extrusion experiments that led him to postulate that yielding occurred when the maximum shear stress at a point reached a critical value.
Von Mises [76] developed a mathematically based yield criterion, which was later explained by Henry [77] as plastic yield when the elastic shear strain energy at a point reached a critical value. Von Mises [76] separately published equations for rigid-perfectly plastic materials, which resembled Levy’s [78] work. In 1913, Von Mises [76] proposed his widely accepted J2 yield criterion that stipulated that yielding took place when critical shear values were attained on octahedral planes (see Nadai’s work [79]). Huber [80] had earlier published essentially the same criterion. In 1948, Hill [81] proposed the first anisotropic yield criterion. Duwez [82] established a stress–strain relationship for plasticity in single-crystal materials at the microstructural level. In that work, the existence of a “secondary structure” in crystals was adopted, where the material is assumed to be composed of N parallel sub-elements. This model was later modified by Besseling [9].
In Figure 1, time-independent plasticity theories are loosely grouped into models based on flow theories and endochronic theories. Plastic flow theories [83,84,85] laid the foundation for the formulation of constitutive equations for plasticity modeling. The 1940s saw the advent of the classical theory; Hill [81], Koiter [83], Prager [84], and Drucker [85], among others, brought together many fundamental aspects of the theory into a single focus. The endochronic theory of plasticity, an intrinsic time theory proposed by Valanis [86,87,88,89] without a yield surface, was able to capture the phenomena of cyclic hardening, initial strain problems, and cross-hardening, unlike classical plasticity. Here, the same constitutive equation is applied for loading and unloading processes; Bazant and Bhat [90] showed that endochronic theories can capture inelasticity and failure in concrete.
Models formulated at a continuum level in a multiaxial setting are useful to engineers designing components in high-temperature applications (nuclear reactors, combustor sections of jet engines, etc.). In continuum theories, stress, strain, and other field quantities are posed at a mathematical point. In applying constitutive models to structural components with internal microstructures, a continuum point must be reinterpreted as a continuum element of finite dimension. The continuum element must be small enough to be nearly homogenous in the presence of stress and strain, as well as temperature gradients, but large enough relative to the microstructure being homogenized. Based on these concepts and well-chosen phenomenological experiments, a system of constitutive relationships can be formulated in support of conducting structural analyses for engineering design. Thus, in this section and succeeding sections, constitutive models have been reviewed that help design components, not the material. Microstructural constitutive models focus more on designing the material. At the continuum level, one designs with the material. At the microstructural level, one designs the material.
Deformation analyses proceed on an incremental basis in classic theories of time-independent plasticity. One calculates increments in plastic strain (Pij). There are three essential elements to a work hardening-based incremental inelastic constitutive law. First is the existence of a threshold function delineating elastic states of stress—a yield criterion. The threshold function typically serves as a potential function. Consider the yield criterion attributed to Tresca [75], von Mises [76], and Hill [81], all of which are stipulated at the continuum level. The second element is a hardening rule, also referred to as an evolutionary law. A hardening rule mathematically describes the evolution of a potential function (i.e., how a material “hardens”) to accommodate the effects of inelastic deformation. The third element is the flow rule. Mendelsohn [91] as well as Chen and Han [92] describe the link between incremental plastic strain and incremental changes in the state of stress captured by the flow rule. Inelastic modeling is primarily based on two types of flow rules. The first is referred to as the associated flow rule, where the threshold function behaves like a potential function. Historically constitutive models for ductile metals are modeled based on this associated flow rule (see Equation (18) below).
The second type is known as a non-associated flow rule, where the inelastic strain rate vector is not normal to the threshold surface. Associated flow rules are more prevalent, so this type of constitutive model will be focused on.
Since classical time-independent incremental plasticity models assume that increments in plastic strains (Pij) are normal to a yield surface. Flow rules, often referred to as normality rules, take the following form:
d ε P i j = λ f σ i j        λ > 0
Thus, an increment in plastic strain is coincident with the gradient of the scalar-valued flow potential function, f. Here, f represents the yield function that is dependent on the stress tensor σij such that
f = f (σij, Hα)
and it incorporates a state variable vector Hα. This vector can contain a single state variable, i.e., K, for isotropic hardening; six state variables associated with kinematic hardening and represented by the symmetric tensor αij; or seven if both isotropic hardening and kinematic hardening are present, i.e.,
Hα= 1, ⋯, 6 = H (K, αij)
A loading rule must be established before inelastic (plastic) strains can be quantified. The loading rule determines whether inelastic (plastic) strains occur along an incremental load path. Since the yield (threshold) function identifies elastic states of stress, stress states outside of the surface of the threshold function are mathematically inaccessible. States of stress within the threshold surface represent elastic states of stress. The inaccessible states of stress outside the boundary of the yield function can be subsequently absorbed within the surface by evolving the boundary of the threshold function as plastic strains accumulate. This evolution process eventually migrates to the functional boundary sufficiently so that stress states beyond the boundary are subsumed. The threshold-potential function therefore must be dependent on stress, as well as the appropriate internal state variables.
A change in the stress state that does not change the inelastic state variable vector Hα corresponds to unloading into or around (tangent to) the elastic stress region. As unloading takes place,
df < 0
and the value of the components of the inelastic state variable vector scalar does not change. When a material is loaded such that a state of stress is on the inelastic threshold surface (f = 0) and an increment in stress, ij, is applied such that
f σ i j d σ i j > 0
then the increment in stress gives rise to an increment in inelastic strain and a corresponding change in the inelastic state of the material, i.e.,
dHα ≠ 0
This change in inelastic state corresponds to the material hardening and in turn impacts the stress–strain curve. The presence of inelastic strains can be detected through the nonlinear behavior of the stress–strain curve or by unloading the material and noting the permanent strains. In general, the change in an inelastic state variable follows the following tensor expression
dHα = Gijα   Pij
The incremental stress ij evolves the threshold function. The inelastic behavior of a material is best described mathematically by
d ε i j p 0 d H α 0 f σ i j , H α = 0 a n d f σ i j d σ i j > 0
or
d ε i j p = 0 d H α = 0 f σ i j + d σ i j , H α < 0 o r f σ i j , H α = 0 a n d f σ i j d σ i j < 0
Specific forms for Pij and dHα depend on the potential function adopted and the hardening laws assumed for the material as inelastic deformations take place. Note that the form of the inelastic potential function can be utilized to extract different types of anisotropic behavior from a plasticity model, as well as different behavior in tension and compression (see Green and Mkrtichian [93]). The theoretical framework for classic plasticity provided above serves as a frame of reference for the time-dependent viscoplastic models discussed next.

6. Time-Dependent Isotropic Unified Models: Viscoplasticity

Predicting time-dependent stress depends not only on the current strain state but also on the previously accrued strains, i.e., the material strain history. Two approaches can be taken in characterizing nonlinear hereditary stress–strain behavior. They involve solving either a system of differential equations or a system of integrals. The first, discussed in detail here, is referred to as internal-state variable representation. Horstemeyer and Bamman [94] provide an excellent historical overview of the internal state variable approach to modeling constitutive behavior. The state variable method develops a system of differential equations that quantify an inelastic strain rate as well as differential equations that govern the evolution of the state of the material during deformation. The viscoplastic model based on a system of differential equations requires a flow rule, an evolutionary law to allow for changes in the state variables adopted, and thermodynamic stability constructs. Adopting a state variable representation is a convenient and flexible thermodynamic format for the description of nonlinear hereditary behavior. The second approach develops a system of integrals to account for stress–strain history. The integral approach, or endochronic theory, represents the inelastic behavior of materials using the notion of intrinsic time. The formulation of a flow rule in endochronic models is not predicated on the use of flow potentials. A third approach is a hybrid approach based on characterizing specific changes in a material’s microstructure during the deformation process. Finite element analysis is employed to develop a representative volume element (RVE) at the microstructural level. Models based on these RVEs can then be used to analyze components at the continuum level once the RVE model is constructed. However, the microstructural models tend to ignore recovery mechanisms that appear explicitly in continuum-level modeling. With the RVE approach, an Arrhenius function and reaction rate theory is applied to a slip process of choice to quantify a strain rate. This third approach will be addressed later in the discussion on how to bridge microstructural concepts and continuum concepts.
As Ponter and Leckie [59] point out, the internal state variable method mirrors many elements of classical plasticity presented in the previous section. They also indicate that internal state variable models are a marriage of classical plasticity concepts with the Norton [3] creep law outlined earlier. A flow law is required to quantify the inelastic strain rates. Rice [60] outlined thermodynamic justification for the existence of a flow potential from which inelastic strain rates are obtained. Ponter and Leckie [59] generalized Norton’s [3] creep law to viscoplasticity for multiaxial states as
d ε i j I d t = σ i j ϕ n + 1 n + 1
where the bracketed term represents Rice’s [60] scalar flow potential function, denoted as Ω in the literature. An evolutionary law (see below) is needed to characterize the change in the internal state variables used to characterize the material microstructure. A third element of an internal state variable model is a criterion stipulating when inelastic deformations will take place. Here, it will be referred to as a threshold function. The threshold function delineates regions of elastic behavior and inelastic behavior in the six-dimensional stress space. For viscoplastic constitutive models based on the internal-state variable method, the threshold function can serve as a flow potential function in Equation (27) above. Threshold functions are key in deriving growth laws that capture various types of hardening mechanisms (e.g., isotropic, kinematic, irradiation) using the internal state variables adopted for a particular model.
The governing differential equations of an internal state variable formulation for viscoplasticity are associated with the normality structure of the threshold potential function Ω. Consider the following the dependence for the threshold potential function suggested by Ponter and Leckie [59]:
Ω = Ω (F, G)
where F is dependent on both the stress state (σij) and the inelastic state variable (αij), i.e.,
F = F (σij, αij)
The second function G
G = G (αij)
is only dependent on inelastic state variables, which are represented by a second-order tensor. Now
Ω = Ω (σij, αij)
and the flow law according to Rice [60] is
ε . i j I = Ω σ i j
Ponter and Leckie [59] stipulated
α i j . = 1 h α i j Ω α i j
where h is a scalar-valued function of the internal state variables only. Equations (32) and (33) stipulate the normality structure of the inelastic strain rate and the rate of change in the inelastic state variable with the threshold potential function Ω.
Based on the dislocation mechanics work of Orowan [95], creep is a competitive process between work hardening, where dislocations in the microstructure interact and local dislocation densities increase and dynamic recovery mechanisms reduce the density of dislocations. Dynamic recovery is a thermally activated mechanism in which immobile dislocations that have piled up are released and the stored energy in the microstructure is diminished. This release corresponds to an obvious decrease in dislocation density. When these two mechanisms balance, the result is a steady state condition. Mitra and McLean [96] suggested how to verify Orowan’s [95] competing process concept experimentally. From a uniaxial stress perspective consider
σ = σ (t, ε)
Taking the differential of this relationship leads to
            d σ = σ t d t + h σ ε d ε = R d t + H d ε
Here, R is the magnitude of the recovery rate and H is a work-hardening coefficient. Note that stress is constant when the material undergoes steady state creep. Thus, for steady state creep conditions,
= 0
As a result, the following relationship can be obtained from Equation (35)
0 = R d t + H d ε ε t s t e a d y s t a t e c r e e p = R H
This relationship was originally proposed in Orowan [95]. It took several decades of experimental work to produce the data that support the mathematical framework in Equation (37). The reader is referred to the seminal efforts of Mitra and McLean [96] as well as others (identified in Gan’s [97] overview) for the experimental data supporting Equation (37). In general, the experimental work carefully measured creep strain rates, the rate of recovery (R), and the hardening coefficient (H) over a range of temperatures and applied stress.
Observations are made in Walker [50] concerning the growth and recovery of the internal state variables due to hardening in the presence and absence of inelastic deformation, respectively. As the growth of the internal state variable due to hardening is dependent on the presence of inelastic deformation, only recovery and annealing changes in the internal state variable are possible during the elastic unloading phase of a thermomechanical cycle. Freed and Robinson [98], Freed and Chaboche [99], and Arnold and Saleeb [100] provide thermodynamic rationales as to the forms of the evolutionary laws for internal state variables. The evolutionary models, along with the flow equations for inelastic strain rate can be utilized to obtain accurate theoretical solutions for a wide variety of time-dependent boundary value problems posed in elevated temperature environments. A general multiaxial formulation for an internal state variable evolutionary law utilized in unified viscoplastic theories takes the following form:
d α i j d t = h α i j d ε i j I d t r α i j a i j
Given the relationships suggested by Rice [60] in Equation (32), and in Equation (33) suggested by Ponter and Leckie [59], obtaining specific formulations for Equation (38) is strictly dependent on the expression used to characterize the flow potential Ω. For each of the viscoplastic models presented in the subsequent sections, the expression for Ω is identified. Viscoplastic models constructed based on flow potentials allow for the convenient use of tensorial invariant theory to extend the models and capture a variety of mechanical responses, e.g., anisotropic behavior and different behavior in tension and compression.
Just as in classical plasticity, the unified viscoplastic models reported on here satisfy Drucker’s [85] postulates. In addition, viscoplastic models can satisfy Ponter’s inequalities [101,102], which are stipulated for time-dependent constitutive models. The postulates and inequalities relate to thermodynamic admissibility restrictions outlined in Arnold and Saleeb [100].

7. Unified Viscoplasticity Models Based on Potential Functions

Fundamentally, the unified approach treats all aspects of inelastic deformation (plasticity, creep, and stress relaxation) with a consistent set of flow equations and evolutionary equations that track the time-dependent behavior of internal state variables. Viscoplasticity models have been used in the nuclear industry to design reactor components. In the aerospace industry, viscoplastic models are used to design combustor section components and leading-edge technologies in hypersonic applications. This approach is robust and has found widespread use.
Allen and Harris [53] point out that several of the viscoplastic constitutive theories they surveyed had similar elements. The models utilize a set of internal state variables that provide locally averaged representations of microstructure such as dislocation rearrangement and grain boundary sliding. The following multiaxial models make the continuum assumption. Unified constitutive equations can be characterized as mathematically “stiff”. The system of coupled partial differential equations used to solve boundary value problems with viscoplasticity models depends on variables that are susceptible to large changes over small time increments of the independent variables. This “stiff” behavior occurs usually with the onset of a significant amount of inelastic strain in the load step and is due primarily to the nonlinear nature of the functional forms.

7.1. Robinson’s Model

Robinson’s model [12,13,14,15,16,17,18,19,20,21] is a unified viscoplastic model based on potential function. This model has well-intended similarities to the structure of the classical plasticity model discussed earlier. The Robinson [15] model possesses a yield criterion (flow potential), and it was originally proposed in a full multiaxial formulation. Under isothermal conditions, the flow potential is dependent on stress, and if kinematic hardening is assumed, the state variable is taken as a deviatoric second-order tensor, i.e., aij. For this model,
Ω = K 2 1 2 μ f F d F + R H g G d G
where μ, R, H, and K are material constants. As the work of Bridgman [103] and others indicates, inelastic deformation is essentially unaffected by the hydrostatic component of the stress state. Robinson [15] originally took stress dependence in terms of the deviatoric components of the applied stress
S i j = σ i j σ k k δ i j 3
and similarly, the second order deviatoric state variable tensor mentioned is
a i j = α i j α k k δ i j 3
Robinson further identified an effective stress as
Σij = Sij − aij
The dependence upon the deviatoric effective stress Sij and the deviatoric internal state variable αij are introduced through the scalar function
F = Fij )
     = F (J2)
where
J 2 = 1 2 i j j i
Similarly,
G = G a i j    = G J 2
where
J 2 = 1 2 a i j a j i
Finally, for the isotropic version of Robinson’s model, the function F takes the following form
F = J 2 K 2 1
and
G = J 2 K 2
F serves as a threshold function. Inelastic strains occur only for states of stress where F is greater than zero. The threshold stress K is generally a scalar state variable and accounts for isotropic hardening (or softening). Note that the concept of a threshold function was introduced by Bingham [104] and later generalized by Hohenemser and Prager [105]. The threshold function is also referred to as the Bingham–Prager threshold function.
Taking the partial derivative of Ω with respect to σij as indicated in Equation (32) with the integral definition of Ω identified in Equation (39) leads to
2 μ d ε i j I d t = f F i j
Robinson [15] then specialized f(F) as a power law function
f (F) = Fn
where n is an additional material constant. Thus,
2 μ d ε i j I d t = F n i j
Similarly, taking the partial derivative of Ω with respect to the internal state variable aij as indicated in Equation (49), the following evolutionary law is obtained for the change in state
d a i j d t = h G d ε i j I d t R H h G g G a i j
Robinson [15] then specialized h(G) as
h G = H G β
and
g (G) = Gm
Thus
d a i j d t = H G β d ε i j I d t R H H G β G m a i j
Equation (55) is consistent with the physically based and well-accepted Bailey–Orowan theory [95,106]. The equation suggests that two competing mechanisms are present that control the evolution of the inelastic state. The first term represents a work-hardening mechanism that proceeds with inelastic strain. The second term represents a recovery process causing a softening that competes with the hardening term. When these two terms balance, the material attains a steady inelastic state.
Equations (51) and (55) comprise the multiaxial statement of the isotropic version of the Robinson viscoplastic constitutive model. The form of the flow potential function above is specialized to a J2 material. Several modifications, generalizations, and improvements of the original Robinson approach have been developed over the years, including the works by Duffy [107] and Saleeb et al. [108].

7.2. Chaboche’s Model

The Chaboche [22,23,24] model is a unified viscoplastic model also based on a potential function. This model includes isotropic and kinematic hardening state variables to capture the Bauschinger effect and cyclic hardening. Here, there are three internal state variables: a second-order tensorial back stress αij to account for kinematic hardening and the other two, yield strength Y and drag stress K, to account for isotropic hardening.
f = 3 2 S i j α i j S i j α i j R K
Note that K represents the yield stress, R is an isotropic hardening state variable, and aij is a back stress (the kinematic state variable). As Chaboche [23] points out, there is a choice in the formulation of the viscoplastic potential Ω through this function f. Equation (56) is a J2 for f. Other tensorial invariants can be used to elicit different responses out of the model (see Duffy [107] and Saleeb et al. [108]). In general, the Chaboche [23] model is presented in a Norton [3] power law formulation as follows:
Ω = K n + 1 f K n + 1
Along with the deviatoric Cauchy stress introduced in the previous section, Chaboche [109] utilizes a deviatoric back stress
a i j = α i j α k k δ i j 3
The inelastic strain rate is obtained from a viscoplastic potential function as follows:
d ε i j I d t = Ω σ i j                = Ω f f σ i j           = λ f σ i j
where the partial derivative of f with respect to σij is a gradient to the inelastic flow surface defined by f. The last form is recognized as a Prandtl–Reuss equation [110,111], where the multiplier λ is the magnitude of the inelastic strain rate vector. The direction of the gradient vector and the associated magnitude of the gradient vector are directly dependent on the choice of the formulation for Ω and f. Given the Norton power law formulation for Ω and a J2 yield function, the flow law takes the form
d ε i j I d t = K 2 2 μ 3 2 S i j α i j S i j α i j R K K n 3 2 S i j α i j 2 3 3 2 S i j α i j 3 2 S i j α i j
All three inelastic state variables evolve through competitive processes that include work hardening, deformation-induced dynamic recovery, and thermally induced state recovery. The evolution of the internal state can also include terms that vary linearly with the external variable rates. For isotropic hardening, the evolutionary law for the single-state variable (K) is
d K d t = b p . Q R
where
p . = 2 3 d ε i j I d t d ε i j I d t
where Q and b are, respectively, the hardening saturation state and the rate of convergence to this state. These are parameters used to define the application of Chaboche’s model to strain range partitioning in fatigue.
For non-isothermal kinematic hardening,
d a i j d t = 2 3 C d ε i j I d t γ a i j p . γ τ T a i j a i j M m 1 a i j
For isothermal kinematic hardening,
d a i j d t = 2 3 C d ε i j I d t γ a i j p .
This is a formulation of an evolutionary law for a tensorial state variable that is similar to Robinson’s [15] given in Equation (52).

8. Viscoplastic Constitutive Models Not Based on Potential Functions

For models described in this subsection, the original formulations for the model were not presented as a derivation from a viscoplastic flow potential function. Subsequent presentations have included flow potentials for some of the models (see, for example, Kim and Oden [112]). Although functions are used to delineate elastic regions of the stress space from inelastic regions, the functions were not utilized to express a potential normality feature.

8.1. Bodner’s Model

Starting with Bodner’s early work [25], a major emphasis in the 1980s focused on combining the materials science flow rules for creep (specifically Norton’s [3]) with flow rules from plasticity. This effort forced different modeling approaches together, and internal state variable theory would join the concepts. Internal state variable theory gained influence as researchers embraced unified-creep-plasticity theories.
Bodner and Partom [26] proposed a flow rule with roots in classical plasticity and influenced by phenomenological observations, but which does not involve the specification of a yield surface. This theory is offered as a means of characterizing the isotropic hardening of certain materials. As in Chaboche’s model [109], where an invariant of the inelastic strain rate is utilized, Bodner et al. [32] adopt a Prandtl–Reuss [110,111] type of flow law:
d ε i j I d t = λ S i j
Squaring both sides of Equation (29) leads to
d ε i j I d t d ε i j I d t = D 2 I                    = λ 2 J 2
or the second invariant of inelastic strain rate is equal to the second invariant of deviatoric stress multiplied by the constant λ2, i.e.,
D 2 I = λ 2 J 2                     = F J 2 , T , Z k
Here, T is temperature and Zk is a vector of inelastic state variables. One of the earliest forms for the flow law was based on a single state variable, Z, and was intended to capture isotropic hardening behavior. Specifically,
D 2 I = D 0 e x p Z 2 3 J 2 n n + 1 n
where D0 is a limiting inelastic strain rate in shear and n impacts strain rate sensitivity, as it does in all models rooted in Norton’s [3] creep law.
Bodner and Partom [27] postulated that the rate of change in a scalar-state variable was a function of inelastic work. For steady-state creep, there must be two competing state variables. Bodner’s model [113] is one of the earliest unified models that captures isotropic hardening through a nonrecoverable isotropic scalar state variable ZI and a directional second-order state variable for kinematic (directional) hardening ZD. These state variables combine linearly:
Z = ZI + ZD
The evolutionary law for the isotropic state variable is
d Z I d t = m 1 Z 1 Z I t d W . I d t A 1 Z 1 Z I t Z 2 Z 1 r 1
Here, m1 is a hardening rate, and
d W I d t = σ i j d ε i j I d t
is the inelastic work rate and is a measure of work hardening. The negative part of Equation (69), representing the first state variable, can be interpreted as dynamic recovery. It is a thermal or “static” recovery of hardening where Z2 is the stable (minimum) value of ZI at a given temperature and A1 and r1 are temperature-dependent material constants.
Kinematic (directional) hardening is represented by a second-order symmetric tensor βij and a second order tensor of direction cosines from the current state of stress, i.e., ZD takes the form
ZD (t) = βij (t) dij (t)

8.2. Walker’s Model

The Walker model [50] was originally developed in an integral form by modifying the constitutive relation for a three-parameter viscoelastic solid. Walker [50] points out that later a differential form of the model was proposed, and this differential format is presented here. Two state variables ωij and K were introduced into the viscoelastic model to account for kinematic (directional) and isotropic hardening. The back stress ωij accounts for kinematic hardening and Bauschinger effects. The isotropic state variable K affects cyclic hardening or softening in a material undergoing repeated loads.
The Walker model [50] is composed of four rate equations and an exponential function. The formulation for the inelastic strain rate depends on the applied deviatoric stress tensor, Sij, and the two state variables ωij and K in the following manner
d ε i j I d t = 2 3 3 S i j 2 ω i j 3 S i j 2 ω i j K n 3 S i j 2 ω i j 2 3 3 S i j 2 ω i j 3 S i j 2 ω i j
The evolutionary law for the back stress state variable is stipulated as
d ω i j d t = n 1 + n 2 + n 1 θ θ t d ε i j I d t ω i j ω ~ i j n 1 ε i j I d G d t 1 n 2 n 2 θ θ t
where n1, n2, and ω i j ~ are material constants. Here, Θ represents temperature, and G is a recovery function. Although not given as a rate equation, the current value of the state variable K depends on the norm of the inelastic strain rate
d R d t = 2 3 d ε i j I d t d ε i j I d t
as follows
K = K1K2 exp(−n3R)
In Equation (74), the time rate of the change in the recovery function G is
d G d t = n 4 + n 5 e x p n 6 R d R d t + n 7 2 ω i j ω i j 3 m 1 2
The model parameters n3, n4, n5, n6, n7, K1, K2, and m are material constants. All material constants in this model are considered temperature-dependent. See Cassenti [114] for details regarding temperature effects impacting the constants in the Walker model [50], the Miller model [36], and the Krieg, Swearengen, and Rhode model [115].

8.3. Miller’s Model

The Miller [36,37,38,39,40,41,42,43,44,45,46,47,48] viscoplastic model captures a wide range of physical phenomena, e.g., cyclic hardening, softening, and saturation behavior. Walker [50], Allen and Harris [53], and Cassenti [114] point out that in the original Miller model, the inelastic strain rate was formulated as
d ε i j I d t = B θ s i n h 3 J 2 K 3 2 n 3 S i j 2 ω i j 3 J 2
and the rate of change for the kinematic state variable is
d ω i j d t = H 1 d ε i j I d t H 1 B θ 1 A 1 2 ω i j ω i j 3 n ω i j 2 ω i j ω i j 3
In these two expressions, B, A1, and H1 are material constants. The rate of change in the scalar state variable is
d K d t = H 2 d R d t C 2 + 2 ω i j ω i j 3 A 2 A 1 K n H 2 C 2 B θ s i n h A 2 K 3 n
In this last expression, A2, C2, and H2, are material constants, and
θ = e x p Q * k T f o r T 0.6 T m θ = e x p Q * k T l n 1 + 0.6 T m T f o r T < 0.6 T m
just as in the Walker [50] model,
d R d t = 2 3 d ε i j I d t d ε i j I d t
The hyperbolic sine function was adopted in this model to capture the creep responses more accurately in a material. Hence, the hyperbolic sine formulation is maintained in the derivation of the thermal-recovery variables of back-stress and drag-stress evolution laws. This particular unified model possesses the ability to predict the Bauschinger effect, history dependence, temperature dependence, anelasticity, and multiaxial response phenomena.

8.4. Krieg, Swearengen, and Rhode’s Model

The model developed by Krieg et al. [115] is one of the first to incorporate both a drag-stress and a back-stress term in a viscoplastic model. The inelastic strain rate is identical to the expression in Walker’s model [50], i.e.,
d ε i j I d t = 2 3 3 S i j 2 ω i j 3 S i j 2 ω i j K n 3 S i j 2 ω i j 2 3 3 S i j 2 ω i j 3 S i j 2 ω i j
However, the evolutionary law for the kinematic state variable is different:
d ω i j d t = A 1 d ε i j I d t A 2 ω i j 2 ω p q ω p q 3 e x p 2 A 3 ω p q ω p q 3 1
The current value of the isotropic state variable K is given by the expression
K = A 4 d R d t A 5 K K 0 n
The value for K depends on the norm of the inelastic strain rate as follows:
d R d t = 2 3 d ε i j I d t d ε i j I d t
No explicit provision exists to model cyclic hardening in either of the state variables. The constants A1 through A5 appearing in the growth laws for the inelastic strain rate and the evolutionary laws for the internal state variables do not explicitly depend on the cumulative inelastic deformation. The linear hardening terms and the recovery term in the state variable evolutionary equation yield stress–strain curves and hysteresis loops that exhibit the same tri-linear character as Miller’s theory.

8.5. Hart’s Model

The Hart [116,117,118] viscoplastic model is one of the early unified creep-plasticity models based on internal state variables. Hart [118] took a microstructural viewpoint in developing his continuum model. He was intent on modeling creep behavior based on experimental evidence of the existence of barriers to dislocations and suggested the existence of a barrier where dislocations “piled up” and a resistance or friction opposing dislocation movement within barriers. Here, the relation between the applied stress, the internal stress, and the glide friction stress was derived as the internal stress was shown to be linearly proportional to a stored anelastic strain. Hart’s model [118] is unique because the isotropic state variable, typically represented as K in other viscoplastic models, is constant. Since this parameter is a constant in the Hart [118] model, it cannot be considered an internal state variable.
The inelastic strain rate is similar to the expression in Walker’s model [50], i.e.,
d ε i j I d t = a * 2 3 3 S i j 2 ω i j 3 S i j 2 ω i j K n 3 S i j 2 ω i j 2 3 3 S i j 2 ω i j 3 S i j 2 ω i j
where Sij is the applied deviatoric stress tensor and ωij is a kinematic internal state variable. The parameter a* is a material constant. The evolutionary law for the kinematic internal state variable is
d ω i j d t = K d ε i j I d t f σ * m ω i j 2 ω p q ω p q 3 l n σ * 2 ω p q ω p q 3 1 β
The variable σ* could be considered a secondary scalar state variable since it serves only to modify the equilibrium stress state variable ωij. The rate of change of this scalar variable is given by the following expression:
d σ * d t = C f σ * m + 1 σ * , ω i j l n σ * ω p q ω p q 2 3 1 β
Note that
σ * , ω i j = γ σ *
although Delph [119] offered other forms for Equation (94). In the above expressions, λ, μ, a*, K, n, f, σ0*, β, C, γ, and δ are material constants that depend on temperature. In addition, σ0* is the initial value of σ*.

9. Integral-Based Viscoplastic Models

The previous viscoplastic models are examples of representing nonlinear constitutive models using a differential equations format. Whereas differential representation is presented in the stress space, the integral representation of material behavior is conducted in the strain/deformation space. As Walker [50] points out, several previous viscoplastic models have both differential and integral formulations. The endochronic models focus on the inelastic response of materials through the use of memory integrals known as memory kernels. Early on, Bouc [120,121] investigated the mathematics of hysteresis, which led to the study of a class of applicable functional operators. In structural engineering, the Bouc–Wen [121,122] model of hysteresis is typically employed to describe nonlinear/inelastic hysteretic systems. In the context of engineering mechanics and structural dynamics, the Bouc–Wen [121,122] model is a phenomenological model that captures both the linear–elastic and elasto-plastic restoring forces in systems that exhibit hysteretic phenomena.
The Bouc–Wen [121,122] model is widely employed for modeling the cyclic behavior of structures in seismic engineering. Specifically, jointed connections dissipate energy through friction over localized regions near the connection and are known to exhibit history-dependent or hysteretic behavior. The Bouc–Wen [121,122] model and subsequent variations are rate-independent hysteresis models. The hysteresis model proposed by Valanis [87] is a rate-dependent hysteresis model. Both the Bouc–Wen [121,122] and Valanis [87] models are grounded in the mechanical modeling concepts found in viscoelasticity. Viscous damping can also be formulated by describing the damping force as a function of the cumulative history of the system, i.e., memory effects. Endochronic models proposed by Valanis [87] and the Bouc–Wen [121,122] are two important examples of hereditary models in endochronic theory that are used in modeling nonlinear material behavior. The Valanis [87] model and its extension to viscoplasticity behavior is focused on here.

Valanis Model

The main principle of a hysteretic model is known as intrinsic time, which relates the deformation history to a deformation memory component. Initially, the endochronic theory introduced by Valanis [87] described a nondecreasing function dependent on the strain tensor εij or the stress tensor σij. The original Valanis [87] endochronic model is built on a single integrated framework and is independent of a yield surface. Resembling the viscoelastic theory, the endochronic theory replaces real time with an auxiliary time variable, referred to as intrinsic time. The constitutive theory developed on the endochronic theory characterizes the hysteresis and strain-hardening behavior of some metals independent of yield conditions, flow rule, or any hardening rules. The intrinsic time-dependent stress evolution rule for the endochronic theory is derived by a convolution integral between the strain tensor εij and the scalar function called a memory kernel. According to endochronic theory, the time-dependent constitutive relationship is defined as
σ i j t = δ i j λ + 2 μ 3 ε k k t + 0 t G z t z ξ ε i j ξ δ i j 3 ε k k ξ d
When the memory kernel G is an exponential function, an incremental form of endochronic flow rules exists. The exponential formulation is commonly used in standard deformation analyses and applications. In general,
G z t z ξ = G 1 e x p a 1 z t z ξ + G 2 e x p a 2 z t z ξ
The function z(t) is defined by the following expression
z t = 1 β l n 1 β R t
where
R t = 0 t θ ξ f θ ξ d ξ
and
θ ξ = 2 3 c i j ξ c i j ξ
Finally,
c i j ξ = i j λ ε k k ξ + 2 μ ε i j ξ k σ i j ξ                  0 < k < 1
And in the expression above, λ, μ, G1, G2, a1, a2, β, and k are material constants.
Several authors, including Bazant and Bhat [90], developed methods of extending endochronic theory to concrete, clay, and sand using the Valanis [87] model as a point of departure. This was carried out by introducing intrinsic time to the hydrostatic component of the stress state and incorporating inelastic dilatancy by way of the shear strain. The Valanis [86,87,88,89] approach to modeling is useful in predicting the nonlinear material behavior of metal-forming crucibles obtained from a powder-metal-forming process. Locating a threshold surface for a high-temperature powder metal can be difficult if the powder metal material is not relatively close to being fully dense. The endochronic theories [123] can be useful from the standpoint that a yield surface is not required to perform an analysis. A subsequent version of the endochronic theory was developed by Valanis [124], where the intrinsic time was redefined as the path length in the inelastic strain space. The new endochronic model could predict a temporal stress response for several deformation processes. Wu and Ho [125] introduced another functional form for the dependency upon the intrinsic time scale to the equivalent deviatoric plastic strain rate and applied endochronic theory to investigate the transient creep of material.

10. Design Scale

In the preceding discussions, modeling inelastic deformation behavior has been presented primarily at a continuum level with references to the physics at the microstructural/mesoscale. Historically, components have been designed at the continuum level, and that design viewpoint still dominates today. The continuum-level design approach is taught to all undergraduate design engineers. However, the continuum-level models must have a basis in physics and a connection to the mesoscale models that have been presented. If that connection to physics is not there, then the continuum-level model is a curve fit. Moreover, if the goal is to improve material performance, then one designs the material at the microstructural level.
The different microstructures within a material (grain morphology, size distributions, anisotropy, and crystallographic orientation; the presence of flaws and porosity; the physical and chemical properties of the intergranular interfaces, etc.) directly influence the performance properties of a material. A link between microstructure and continuum can provide valuable insight into the design of components fabricated from high-performance materials. It should be noted that the design perspective at which a material/component/system is modeled strongly depends on the experimental information available. Internal state variable (ISV) concepts have been formalized as a continuum-level parameter that accounts for dislocation mechanics. This provides a link to the prediction of the movement of dislocations through a microstructure. System-level parameters such as flexibility and stiffness coefficients are based on continuum-level parameters and geometry.
Relative to time-dependent behavior, constitutive models are available at both levels with convenient bridges identified by researchers over the years. In the next section, a discussion is provided on those links to different design scales.

11. The Bridge from Microstructure to Continuum: Time-Dependent Behavior

ISVs account for local unified creep and plasticity mechanisms driving time-dependent deformation behavior at the microstructural level. Coleman and Gurtin [126] connected the behavior of a material under applied load at the microstructural level with the incorporation of the internal state variables in the material model. Lubliner [127,128,129] also showed that a material with a property of fading memory can be described with internal state variables. The fading memory property refers to the hypothesis that the values of stress are dependent on recent deformations when compared to the deformations incurred in the distant past.
To capture the damage in materials, Horstemeyer et al. [130] postulated an ISV theory combining the hardening equations and porosity evolution equations. They pointed out how a finite element model that explicitly incorporates ISVs can be used to capture the structure–property relations in the realm of large deformations. A finite element simulation accommodates different grain sizes, particle sizes, pore sizes, and volume fractions within each element so that a material with a heterogeneous microstructure can be represented within the whole mesh. Based on the plasticity/damage approach proposed by Bamman et al. [131], Horstemeyer et al. [130] made improvements to the theory by adding ISV equations that introduce grain particle size, particle volume fraction, grain size, pore volume fraction, and nearest neighbor distances of the particles and pores. This allowed the use of heterogeneous distributions of microstructures throughout a finite element mesh. This multiscale modeling approach demonstrated that the standard assumption of homogeneous distributions of microstructural features such as porosity (i.e., the classical approach) can lead to erroneous results and conclusions. This is the benefit of using microstructural property-based ISV theories in conjunction with finite element analysis.
The deformation process in a material generally involves nonuniform fields present at the microstructural level. As noted earlier, Rice [132] examined the relationship between inelastic strain and slip displacements. He concluded that the slip within individual grains may be considered a deformation in the usual continuum sense and that each point of the inelastic strain can be broken down into a finite set of simple shearing strains of the respective slip planes and directions.
Thus, at the microstructural level, the material can be considered an ensemble of distinct entities, e.g., individual grains in polycrystalline metals. Many of the early continuum-level viscoplastic models adopted the Rice [132] interpretation to make a connection between the continuum level and the microstructure.
One must be cognizant of the material evolution caused by dislocation motion taking place in the material microstructure, and this can be achieved by studying well-accepted material science models. Blum and Eisenlohr [133] point out that insight obtained from material science models relative to how dislocation structures and microstructural morphology impact deformation behavior should be reflected in continuum-level models. However, material science models tend to not paint a complete picture in the sense that primary creep (the transient response) is captured in some material science models and ignored in other models that focus on steady-state response. Material science models are typically not posed relative to multiaxial stress states. Several schools of thought exist on the topic of representing material behavior in a fully three-dimensional setting with microstructural-level information explicitly incorporated. Consider an approach that employs finite element analysis to develop a representative volume element (RVE) that captures microstructural level features. Hill [134] first suggested how to construct models that can be constructed based on these RVEs that analyze components at the continuum level in a fully three-dimensional setting. The RVE approach has its merits in that the physics of dislocation mechanics is captured. Modeling effort using the RVE approach requires an interpretation of boundary conditions, as well as knowledge of initial dislocation densities.
The ability to observe, characterize, simulate, and then design components at the mesoscale level requires advances in techniques. For example, to observe how the microstructure of the test specimen evolves, High Energy X-ray Diffraction Microscopy (HEDM) is the frontrunner because of its ability to look at the mesoscale structure. This ability to see with microscale resolution helps in characterizing the polycrystal orientation and strain states under thermal and mechanicals loads. There is a need to characterize the initial dislocation field in situ to model from the microstructural level up to the continuum level. The challenge is in defining the right boundary conditions; in the case of a single crystal with free surfaces, free boundaries (dislocations are free to vanish when they reach the surface) are defined or closed boundaries (the dislocations can penetrate the surface) are defined to capture the behavior of a polycrystal. At the microstructural level, the dislocations used to define the initial dislocation density are randomly distributed with respect to their length and number. These definitions, integrated over the simulated region, can be used to derive the global quantities, e.g., the stress–strain curves, dislocation densities, and local qualities calculated at every step of the simulation. In some regions, these stress states must be allowed to vary, where the trajectories of the dislocation path and the glide plane are tracked.
There is a need to bridge the two scales and master the art of assembling structural and functional microstructural models into larger continuum models that allow for the design of complex structural components. Polycrystalline materials are among the most important due to their ability to distribute load via plastic deformation. The two main mechanisms through which plastic strains are expressed are the flow of dislocations and twinning, both of which are characterized as in-line defects in crystals. These mechanisms are understood individually, but their interactions that affect the behavior at the mesoscale level remain a challenge. For example, elastic moduli are calculated using first principles but the computation for stress–strain curves, as required by engineers, is a challenge because of the complicated dislocation mechanics at the atomistic level.
The phenomenological models discussed in earlier sections are continuum-level deformation models. The phenomenological constitutive models based on dislocation theories incorporate a temporal relationship between continuum level stress and strain. Currently derived phenomenological models do not explicitly account for the interaction between material microstructure and the continuum. However, well-thought-out phenomenological models should in some fashion mirror the evolution of the material at the microstructural level mathematically. Phenomenological models may not explicitly capture in detail the material behavior within the microstructure. However, the question should be asked as to what level should a component be designed at. Should one start the design of components at the atomic/molecular level, or should one appeal to nanomechanics/micromechanics, or the mesomechanics continuum level? Historically, engineering design has been cognizant of material behavior at the microstructural level, but component designs are executed at a continuum level. To begin a discussion on bridging these schools of thought, a brief overview is provided regarding the details supporting representative volume elements (RVE).
There are always gaps between the modeling efforts at one level relative to the next level. A prevalent philosophy has been that the design engineer designs components using phenomenological models and continuum mechanics, while material scientists design the material with microstructural models. The modeling at both levels should have a symbiotic relationship. Bodner and Partom [27] attempted to bridge the gap between the material microstructure and the continuum by including certain physical concepts supporting dislocation mechanics. Others followed by similarly pointing to physical mechanisms that their continuum models capture. Here, a transitional approach is discussed where the physics of the microstructure can be captured by a representative volume element. Instead of pointing to microstructural mechanisms and attempting to mirror the mechanics, the intent is to explicitly incorporate the mechanisms via a discretized finite element mesh. This approach evolved with the intent to develop a collaboration between material science and engineering mechanics. Mesomechanics allows the adoption of a heterogeneous medium with a distinct microstructure and micromechanics that evolve the microstructure as opposed to a continuum with averaged properties. The limitations associated with mesomechanics include the need for a precise description of the microstructure and the inability to capture the complexities in the evolution of complex material structures.
Non-linear deformations and time-dependent material behavior based on the concept of representative volume element were first adopted by Hill [134]. Since finite element analysis is based on the continuum, mechanics and small deformations are typically used under the assumption of linear deformations. To circumvent the limitations of heterogeneity and linear processes, a few methods have been proposed in the form of asymptomatic homogenization theory by Terada and Kikuchi [135] and Ghosh et al. [136,137] and volume averaging by Smit et al. [138] and Feyel and Chaboche [139]. In contrast to assigning homogenized constitutive behavior at integration points, the average behavior associated with the microstructure is assigned at each integration point of the macromesh of the representative volume element (RVE).
To achieve precise detail at the microstructural level and to better understand the link between microstructure and continuum, a polycrystal plasticity finite element model was proposed by Holm et al. [140]. The authors addressed the limitations of classical constitutive models by allowing an atoms-up approach that also produces a continuum down the path. Their model includes geometry for pores as well as the grain structure. The authors use a paradigm where a combined computational and experimental assessment of the microstructure leads to the homogenization of the material. The homogenized properties are then utilized in continuum level engineering design. Digital image correlation coupled with electron-backscatter diffraction microscopy provide strain maps that can validate microstructure. Characteristic features and their distributions and trends are revealed that validate simulation results, e.g., plastic strain bands that appear oriented at 45° in components subjected to uniaxial loads.
A discrete element method proposed by Horner et al. [141], along with Iwashita and Oda [142] in collaboration with Tordesillas [143], suggests the adoption of the smallest possible representative volume element (RVE) to ensure high resolution at the microstructural level to capture the dynamics and kinematics of the particles to link to the evolution of the microstructural mechanics.
The micromorphic continuum theory is outlined in Mindlin [144], where the material is observed as a collection of deformable “points” used to model material parameters and generate models at the nanoscale. This was accomplished by using numerous bridging scale methods in conjunction with continuum, quantum, and molecular mechanics as demonstrated in Liu et al. [144,145]. To further the application of these methods, a new theory called the multiresolution continuum theory was proposed by Liu et al. [144,145] for heterogeneous materials. Gao et al. [146] expanded the micromorphic theory by incorporating geometric relevant material defects.
Finally, by incorporating damage mechanics at the microstructural level in the classical continuum models also helps bridge the gap and capture the material behavior at several sub-levels, as was proposed in the work of Holm [147], who continued the development of future multiscale modeling theories by adopting various approaches. The upscaling method and resolved-scale method are reviewed in Fish et al. [148]. Holm et al. [148] state that understanding various processes of microstructural scale damage where the interaction between the grain boundary with dislocations and twins is important and critical in modeling damage mechanics. Continuum models fail to predict fracture. The goal is that the microscale damage models can bridge this gap. Methods have been developed to bridge the gap across length and time scales, as discussed in Fish et al. [149].

12. Concluding Remarks

For over a century, the materials community has been proposing methods to characterize the nonlinear stress–strain behavior exhibited by materials. In this review article, a collection of works is presented that depicts the development of inelastic constitutive models that have applications at the continuum level, a homogenized polycrystalline material. The discussion begins with classical plasticity theories and that discussion carries through to time-dependent viscoplastic models.
Focus has been given to material science models that physically explain nonlinear behavior at the microstructural level. An understanding of material microstructure is always necessary in developing accurate multiaxial continuum-level constitutive models that characterize the responses of engineering components modeled with continuum-level perspectives.
For further reading, attention is drawn to the work of Vladmir Buljak and Gianluca Ranzi [149], which captures the distinction between various inelastic constitutive models from a different perspective.
A bridge between microstructural concepts and continuum concepts is convenient. This bridge is discussed in the context of time-dependent behavior constitutive models.
Many forms of inelastic constitutive models are presented that include differential formulations as well as integral forms. Models are available that aid in the design of components that function in high-temperature environments. This would be accomplished by incorporating three-dimensional continuum-level inelastic models in finite element analysis software. A discussion regarding that effort is available elsewhere.

Funding

This research was sponsored by the Advanced Reactor Technologies program of the U.S. Department of Energy. Oak Ridge National Laboratory is managed by UT-Battelle LLC under contract number DE-AC05-00OR22725 for the U.S. Department of Energy.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Andrade, E.N.D.C. On the Viscous Flow in Metals, and Allied Phenomena. Proc. R. Soc. Sect. A–Math. Phys. Sci. 1910, 84, 1–12. [Google Scholar]
  2. Andrade, E.N.D.C. The Flow in Metals under Large Constant Stresses. Proc. R. Soc. Sect. A–Math. Phys. Sci. 1914, 90, 329–342. [Google Scholar]
  3. Norton, F.H. The Creep of Steel at High Temperatures; McGraw-Hill: London, UK, 1929. [Google Scholar]
  4. Bailey, R.W. Creep of Steel Under Simple and Compound Stress. Engineering 1930, 121, 129–265. [Google Scholar]
  5. Odqvist, F.K.G. International Congress of Applied Mechanics, Cambridge, Proceedings; Cambridge University Press: Cambridge, UK, 1935; pp. 228–229. [Google Scholar]
  6. Odqvist, F.K.G. Historical Survey of the Development of Creep Mechanics from its Beginnings in the Last Century to 1970. In Creep in Structures, IUTAM, 3rd Symposium; Springer: Leicester, UK, 1980; pp. 1–12. [Google Scholar]
  7. Schapery, R.A. Nonlinear Viscoelastic and Viscoplastic Constitutive Equations Based on Thermodynamics. Mech. Time-Depend. Mater. 1997, 1, 209–240. [Google Scholar] [CrossRef]
  8. Findley, W.N.; Lai, J.S.; Onaran, K. Creep and Relaxation of Nonlinear Viscoelastic Materials; Dover Publications Inc.: New York, NY, USA, 1976; ISBN 978-0486660165. [Google Scholar]
  9. Besseling, J.F. Theory of elastic, plastic, and creep deformations of an initially isotropic material showing anisotropic strain-hardening, creep recovery, and secondary creep. J. Appl. Mech. 1958, 25, 529. [Google Scholar] [CrossRef]
  10. Onimus, F.; Jourdan, T.; Xu, C.; Campbell, A.A.; Griffiths, M. Irradiation Creep in Materials. In Comprehensive Nuclear Materials, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2021; pp. 301–366. [Google Scholar]
  11. Jaske, C.E.; Leis, B.N.; Pugh, C.E. Monotonic and Cyclic Stress-Strain Response of Annealed 21/4Cr-1Mo Steel. In Symposium on Structural Materials for Service at Elevated Temperatures in Nuclear Power Generation; Schaefer, A.O., Ed.; ASME Winter Annual Meeting: New York, NY, USA, 1975. [Google Scholar]
  12. Pugh, C.E.; Robinson, D.N. Some Trends in Constitutive Equation Model Development for High-Temperature Behavior of Fast-Reactor Structural Alloys. Nucl. Eng. Des. 1978, 48, 269. [Google Scholar] [CrossRef]
  13. Robinson, D.N. ORNL TM-5110; Candidate Creep-Recovery Model for 2 1/4 Cr-1 Mo Steel and Its Experimental Implementation. Energy Research and Development Administration: Oak Ridge, TN, USA, 1975.
  14. Robinson, D.N. ORNL-5136; Developments Toward Refined Constitutive Laws for Reactor System Metals. Energy Research and Development Administration: Oak Ridge, TN, USA, 1975; pp. 15–23.
  15. Robinson, D.N. ORNL TM-5571; On the Concept of a Flow Potential and the Stress-Strain Relations of Reactor System Metals. Energy Research and Development Administration: Oak Ridge, TN, USA, 1976.
  16. Robinson, D.N. ORNL TM-5235; Tests for Examining the Concept of a Flow Potential in the Stress-Strain Relations of Reactor System Metals. Energy Research and Development Administration: Oak Ridge, TN, USA, 1976.
  17. Robinson, D.N. ORNL TM-5339; Developments toward Refined Constitutive Laws. Energy Research and Development Administration: Oak Ridge, TN, USA, 1977.
  18. Robinson, D.N. ORNL TM-5969; A Unified Creep-Plasticity Model for Structural Metals at High Temperature. Energy Research and Development Administration: Oak Ridge, TN, USA, 1978.
  19. Robinson, D.N.; Swindeman, R.W. ORNL TM-8444; Unified Creep-plasticity Constitutive Equations for 2-1/4 Cr-1Mo Steel at Elevated Temperature. Energy Research and Development Administration: Oak Ridge, TN, USA, 1982.
  20. Robinson, D. N Constitutive Relationships for Anisotropic High Temperature Alloys. Nucl. Eng. Des. 1984, 83, 389–396. [Google Scholar] [CrossRef]
  21. Robinson, D.N.; Bartolotta, P.A. NASA CR-174836; Viscoplastic Constitutive Relationships with Dependence on Thermomechanical History. Energy Research and Development Administration: Oak Ridge, TN, USA, 1985.
  22. Chaboche, J.L. Bulletin de L’Academie des Sciences. Série Des Sci. Tech. 1977, 1, 33. [Google Scholar]
  23. Chaboche, J.L. A Review of Some Plasticity and Viscoplasticity Constitutive Theories. Int. J. Plast. 2008, 1, 33. [Google Scholar] [CrossRef]
  24. Chaboche, J.L.; Nouailhas, D. On Various Nonlinear Kinematic Hardening Rules in Cyclic Plasticity and Viscoplasticity. In Proceedings of the Computational Mechanics ’88: Theory and Applications: Proceedings of the International Conference on Computational Engineering Science, Atlanta, GA, USA, 10–14 April 1988; pp. 19.i.1–19.i.4. [Google Scholar]
  25. Bodner, S.R. Constitutive Equations for Dynamic Material Behavior. In Mechanical Behavior of Materials Under Dynamic Loading; U.S. Lindholm, Ed.; Springer: New York, NY, USA, 1968; pp. 176–190. [Google Scholar]
  26. Bodner, S.R.; Partom, Y. A Large Deformation Elastic-Viscoplastic Analysis of a Thick-Walled Spherical Shell. ASME J. Appl. Mech. 1972, 39, 751–757. [Google Scholar] [CrossRef]
  27. Bodner, S.R.; Partom, Y. Constitutive Equations for Elastic-Viscoplastic Strain Hardening Materials. ASME J. Appl. Mech. 1978, 42, 385–389. [Google Scholar] [CrossRef]
  28. Bodner, S.R.; Partom, I.; Partom, Y. Uniaxial Cyclic Loading of Elastic-viscoplastic Materials. ASME J. Appl. Mech. 1979, 46, 805–810. [Google Scholar] [CrossRef]
  29. Bodner, S.R. A Procedure for Including Damage in Constitutive Equations for Elastic- viscoplastic Work-hardening Materials. In Proceedings of the IUTAM Symposium on Physical Nonlinearities in Structural Analysis, Senlis, France, 27–30 May 1981; pp. 21–28. [Google Scholar]
  30. Stouffer, D.C.; Bodner, S.R. A Relationship Between Theory and Experiment for a State Variable Constitutive Equation. In Mechanical Testing for Deformation Model Development; ASTM STP: West Conshohocken, PN, USA, 1982; Volume 765, pp. 239–250. [Google Scholar]
  31. Bodner, S.R. Evolution Equations for Anisotropic Hardening and Damage of Elastic-Viscoplastic Materials. In Proceedings of the Conference on Plasticity Today, Udine, Italy, 13 July 1983. [Google Scholar]
  32. Bodner, S.R.; Stouffer, D.C. Comments on Anisotropic Plastic Flow and Incompressibility. Int. J. Eng. Sci. 1983, 21, 211–215. [Google Scholar] [CrossRef]
  33. Lindholm, U.S.; Chan, K.S.; Bodner, S.R.; Weber, R.M.; Walker, K.P.; Cassenti, B.N. Constitutive Modeling for Isotropic Materials (HOST). In First Annual Contract Report NASA CR-174718; NASA: Washington, DC, USA, 1984. [Google Scholar]
  34. Lindholm, U.S.; Chan, K.S.; Bodner, S.R.; Weber, R.M.; Walker, K.P.; Cassenti, B.N. Constitutive Modelling for Isotropic Materials (HOST). In 2nd Annual Contract Report, NASA CR-174980; NASA: Washington, DC, USA, 1985. [Google Scholar]
  35. Chan, K.S.; Bodner, S.R.; Lindholm, U.S. Phenomenological Modeling of Hardening and Thermal Recovery in Metals. J. Eng. Mater. Technol. Am. Soc. Mech. Eng. 1988, 110, 1–8. [Google Scholar] [CrossRef]
  36. Miller, A.K. An Inelastic Constitutive Model for Monotonic, Cyclic and Creep Deformation: Part 1–Equation Development and Analytical Procedures. J. Eng. Mater. Technol. 1976, 99, 97–105. [Google Scholar] [CrossRef]
  37. Miller, A.K. Modelling of Cyclic Plasticity: Improvements in Simulating Normal and Anomalous Bauschinger Effects. J. Eng. Mater. Technol. 1980, 102, 215–220. [Google Scholar] [CrossRef]
  38. Schmidt, C.G.; Miller, A.K. A Unified Phenomenological Model for Non-elastic Deformation of Type 316 Stainless Steel-Part I: Development of the Model and Calculation of the Material Constants. Res. Mech. 1981, 3, 109–129. [Google Scholar]
  39. Schmidt, C.G.; Miller, A.K. A Unified Phenomenological Model for Non-elastic Deformation of Type 316 Stainless Steel-Part II: Fitting and Predictive Capabilities. Res. Mech. 1981, 3, 175–193. [Google Scholar]
  40. Ruano, O.A.; Miller, A.K.; Sherby, O.D. The Influence of Pipe Diffusion on the Creep of Fine-grained Materials. Mater. Sci. Eng. 1981, 51, 9–16. [Google Scholar] [CrossRef]
  41. Kassner, M.E.; Rubin, K.A.; Miller, A.K. Verification of a Microstructurally based Equation for Elevated-temperature Transient Isotropic Hardening. In Strength of Metals and Alloys; Pergamon Press: Oxford, UK, 1982; Volume 2, pp. 581–587. [Google Scholar]
  42. Miller, A.K.; Ziaai-Moayyed, A.A. Some Critical Experimental Tests of the MATMOD Constitutive Equations with Respect to Directional Hardening and Cyclic Deformation. In Mechanical Testing for Deformation Model Development, STP 765; ASTM: West Conshohocken, PN, USA, 1982; pp. 202–222. [Google Scholar]
  43. Miller, A.K.; Kassner; Sherby, O.D. The Separate Roles of Subgrains and Forest Dislocations in the Isotropic Hardening of Type 304 Stainless Steel. Met. Trans. 1982, A, 13A. [Google Scholar]
  44. Schmidt, C.G.; Miller, A.K. The Effect of Solutes on the Strength and Strain Hardening Behavior of Alloys. Acta Met. 1982, 30, 615–625. [Google Scholar] [CrossRef]
  45. Lowe, T.C.; Miller, A.K. Improved Constitutive Equations for Modelling Strain Softening-Part 1: Conceptual Development and Part 2: Predictions for Aluminum. J. Eng. Mater. Technol. 1984, 106, 337–348. [Google Scholar] [CrossRef]
  46. Hedling, D.E.; Miller, A.K. The Incorporation of Yield Surface Distortion into a Unified Constitutive Model, Part I: Equation Development. Acta Mech. 1987, 69, 9–23. [Google Scholar] [CrossRef]
  47. Miller, A.K. Unified Constitutive Equations for Creep and Plasticity; Elsevier: London, UK, 1987. [Google Scholar]
  48. Miller, A.K.; Obalueki, A.O.; Lee, C.W.; Tanaka, T.G.; Lee, S.B. A Unified Model for Fatigue Crack Initiation and Growth, with Emphasis on Short-Crack Behavior, Crack Closure Effects Variable–Temperature Fatigue and Creep-Fatigue Interaction. Mat. Sci. Eng. 1988, A103, 71–93. [Google Scholar] [CrossRef]
  49. Walker, K.P. Representation of Hastelloy-X Behavior at Elevated Temperature with a Functional Theory of Viscoplasticity. In Proceedings of the ASME Pressure Vessels Conference, San Francisco, CA, USA, 12–15 August 1980. [Google Scholar]
  50. Walker, K.P. NASA CR-165533; Research and Development Program for Non-linear Structural Modelling with Advanced Time-temperature Dependent Constitutive Relationships. NASA: Washington, DC, USA, 1981.
  51. Freed, A.D.; Walker, K.P. NASA TM 102338; Refinements in a Viscoplastic Model. NASA: Washington, DC, USA, 1989.
  52. Chan, K.S.; Lindholm, U.S.; Bodner, S.R.; Walker, K.P. High Temperature Deformation Under Uniaxial Loading: Theory and Experiment. J. Eng. Mater. Technol. ASME 1989, 111, 345–353. [Google Scholar] [CrossRef]
  53. Allen, D.H.; Harris, C.E. NASA TM 102727; A Review of Models for Nonlinear Constitutive Models for Metals. NASA: Washington, DC, USA, 1990.
  54. Chamis, C.C. NASA TM 83320; Simplified Composite Micromechanics Equations for Hygral, Thermal and Mechanical Properties. NASA: Washington, DC, USA, 1983.
  55. Coble, R.L. A Model for Boundary Diffusion Controlled Creep in Polycrystalline Materials. J. Appl. Phys. 1963, 34, 1679–1682. [Google Scholar] [CrossRef]
  56. Nabarro, F.R.N. Report of Conference on Strength of Solids. Phys. Soc. 1948, 2, 75–90. [Google Scholar]
  57. Herring, C. Diffusional Viscosity of a Polycrystalline Solid. J. Appl. Phys. 1950, 21, 437. [Google Scholar] [CrossRef]
  58. Weertman, J.R. Theory of steady-state creep based on dislocation climb. J. Appl. Phys. 1955, 26, 1213–1217. [Google Scholar] [CrossRef]
  59. Ponter, A.R.S.; Leckie, F.A. Constitutive Relationships for the Time Dependent Deformation of Metals. ASME J. Eng. Mater. Technol. 1976, 98, 47–51. [Google Scholar] [CrossRef]
  60. Rice, J.R. On the Structure of Stress-Strain Relations for Time Dependent Plastic Deformation in Metals. J. Appl. Mech. 1970, 37, 728–737. [Google Scholar] [CrossRef]
  61. Was, G.S. Fundamentals of Radiation Materials Science-Metals and Alloys, 2nd ed.; Springer: New York, NY, USA, 2017. [Google Scholar]
  62. Onimus, F.; Jourdan, T.; Xu, C.; Campbell, A.A.; Griffiths, M. Chapter 1.10: Irradiation Creep in Materials. In Comprehensive Nuclear Materials, 2nd ed.; Konings, R.J.M., Stoller, R.E., Eds.; Elsevier: Oxford, UK, 2020; pp. 310–366. [Google Scholar]
  63. Roberts, A.C.; Cottrell, A.H. LXXII. Creep of Alpha Uranium During Irradiation with Neutrons. A J. Theor. Exp. Appl. Phys. 1956, 1, 711–717. [Google Scholar] [CrossRef]
  64. Hesketh, R.V. A Possible Mechanism of Irradiation Creep and Its Reference to Uranium. Philos. Mag. A J. Theor. Exp. Appl. Phys. 1962, 7, 1417–1420. [Google Scholar] [CrossRef]
  65. Koyanagi, T.; Katoh, Y.; Ozawa, K.; Shimoda, K.; Hinoki, T.; Snead, L.L. Neutron-Irradiation Creep of Silicon Carbide Materials Beyond the Initial Transient. J. Nucl. Mater. 2016, 478, 97–111. [Google Scholar] [CrossRef]
  66. Campbell, A.A. Historical Experiment to Measure Irradiation-Induced Creep of Graphite. Carbon 2018, 139, 279–288. [Google Scholar] [CrossRef]
  67. Hesketh, R.V. The Mechanisms of Irradiation Creep in Graphite. Philos. Mag. 1965, 11, 917–927. [Google Scholar] [CrossRef]
  68. Kelly, B.T.; Brocklehurst, J.E. ORNL/NPR-92-58; Analysis of Irradiation Creep in Reactor Graphite. Energy Research and Development Administration: Oak Ridge, TN, USA, 1993.
  69. Kelly, B.T.; Brocklehurst, J.E. UKAEA Reactor Group Studies of Irradiation-Induced Creep in Graphite. J. Nucl. Mater. 1977, 65, 79–85. [Google Scholar] [CrossRef]
  70. Kelly, B.T.; Foreman, A.J.E. The Theory of Irradiation Creep in Reactor Graphite–The Dislocation Pinning-Unpinning Model. Carbon 1974, 12, 151–158. [Google Scholar] [CrossRef]
  71. Campbell, A.A.; Was, G.S. Proton Irradiation-Induced Creep of Ultra-Fine Grain Graphite. Carbon 2014, 77, 993–1010. [Google Scholar] [CrossRef]
  72. Sarkar, A.; Eapen, J.; Raj, A.; Murty, K.L.; Burchell, T.D. Modeling Irradiation Creep of Graphite Using Rate Theory. J. Nucl. Mater. 2016, 473, 197–205. [Google Scholar] [CrossRef]
  73. Heggie, M.I.; Suarez-Martinez, I.; Davidson, C.; Haffenden, G. Buckle, Ruck and Tuck: A Proposed New Model for the Response of Graphite to Neutron Irradiation. J. Nucl. Mater. 2011, 413, 150–155. [Google Scholar] [CrossRef]
  74. Coulomb, C.A. Essai sur une application des regles des maximis et minimis a quelquels problemesde statique relatifs, a la architecture. Mem. Acad. Roy. Div. Sav. 1776, 7, 343–387. [Google Scholar]
  75. Tresca, H.E. Mémoire sur l’écoulement des corps solides soumis à de fortes pressions. Comptes Rendus De L’academie Des Sci. 1864, 59, 754–759. [Google Scholar]
  76. Von Mises, R. Mechanik der festen Körper im plastisch deformablen Zustand. Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen. Math.-Phys. Kl. 1913, 582–592. [Google Scholar]
  77. Hencky, H. Zur Theorie plasticher Deformationen und der hierdurch im Material hervorgerufenen Nachspannungen. In Proceedings of the Conference of Applied Mechanics, Delft, The Netherlands, 29 February 1924. [Google Scholar]
  78. Levy, M. Extrait du mémoire sur les equations générales des mouvements intérieures des corps solides ductiles au dela des limites ou l’élasticité pourrait les ramener à leur premier état. J. Math. Pures. Appl. 1871, 16, 369–372. [Google Scholar]
  79. Nadai, A. Theory of Flow and Fracture of Solids, 2nd ed.; McGraw-Hill Book Co., Inc.: New York, NY, USA, 1950; Volume 1. [Google Scholar]
  80. Huber, M.T. Specific Strain Work as a Measure of Material Effort (in Polish: Właściwa Praca Odkształcenia Jako Miara Wytężenia Materyału; Czasopismo Techniczne: Lwów, Ukraine, 1904; pp. 34–40, 49–50, 61–62, 80–81. [Google Scholar]
  81. Hill, R. The Mathematical Theory of Plasticity. Bull. Am. Math. Soc. 1952, 58, 507–512. [Google Scholar]
  82. Duwez, P. On the Plasticity of Crystals. Phys. Rev. 1935, 47, 494–501. [Google Scholar] [CrossRef]
  83. Koiters, W.T. Stress-Strain Relations, Uniqueness and Variational Theorems for Elastic-Plastic Materials with a Singular Yield Surface. Q. Appl. Math. 1953, 11, 350–354. [Google Scholar] [CrossRef]
  84. Prager, W. The Theory of Plasticity: A Survey of Recent Achievements. Proc. Inst. Mech. Eng. 1955, 169, 41–57. [Google Scholar] [CrossRef]
  85. Drucker, D.C. A Definition of Stable Inelastic Material. J. Appl. Mech. 1959, 26, 101–106. [Google Scholar] [CrossRef]
  86. Valanis, K.C. Unified Theory of Thermomechanical Behavior of Viscoelastic Materials. In Mechanical Behavior of Materials under Dynamic Loads; Springer: Berlin, Germany, 1968; p. 343. [Google Scholar]
  87. Valanis, K.C. A Theory of Viscoplasticity Without a Yield Surface Part I. General Theory. Arch. Mech. 1971, 23, 517–533. [Google Scholar]
  88. Valanis, K.C. A Theory of Viscoplasticity Without a Yield Surface Part II. Application to Mechanical Behavior of Metals. Arch. Mech. 1971, 23, 535–551. [Google Scholar]
  89. Valanis, K.C. Report G-224/DME-78-01; Fundamental Consequences of a New Intrinsic Time Measure: Plasticity as a Limit of the Endochronic Theory. Division of Materials Engineering-The University of Iowa (NSF Supported): Iowa City, IA, USA, 1978.
  90. Bazant, Z.P.; Bhat, P.D. Endochronic theory of inelasticity and failure of concrete. ASCE J. Eng. Mech. Div. 1976, 102, 701–722. [Google Scholar] [CrossRef]
  91. Mendelsohn, A. Plasticity: Theory and Applications; MacMillan: New York, NY, USA, 1968. [Google Scholar]
  92. Chen, W.F.; Han, D.J. Plasticity for Structural Engineers; Gau Lih Book Co. Ltd.: Taipei, Taiwan, 1995. [Google Scholar]
  93. Green, A.E.; Mkrtichian, J.Z. Elastic Solids with Different Moduli in Tension and Compression. J. Elast. 1977, 7, 369–386. [Google Scholar] [CrossRef]
  94. Horstemeyer, M.F.; Bamman, D.J. Historical Review of Internal State Variable Theory for Inelasticity. Int. J. Plast. 2010, 26, 1310–1334. [Google Scholar] [CrossRef]
  95. Orowan, E. The Creep of Metals. J. West. Scotl. Iron Steel Inst. 1946, 54, 45–96. [Google Scholar]
  96. Mitra, S.R.; McLean, D. Work Hardening and Recovery in Creep. In Proceedings of the Royal Society, London, UK, 24 October 1961. [Google Scholar]
  97. Gan, D. A study of the Bailey-Orowan Equation of Creep. J. Mater. Sci. 1982, 17, 89–99. [Google Scholar] [CrossRef]
  98. Freed, A.D.; Robinson, D.N. A Theory of Viscoplasticity Accounting for Internal Damage. In Proceedings of the Pressure Vessels and Piping Conference, PVP 123, San Diego, CA, USA, 23 June 1987; pp. 119–136. [Google Scholar]
  99. Freed, A.D.; Chaboche, J.L. NASA TM-101288; Viscoplasticity: A Thermodynamic Formulation. NASA: Washington, DC, USA, 1989.
  100. Arnold, S.; Saleeb, A. On the Thermodynamic Framework of Generalized Coupled Thermoelastic–Viscoplastic Damage Modeling. Int. J. Plast. 1994, 10, 263–278. [Google Scholar] [CrossRef]
  101. Ponter, A.R.S. Dynamic Behavior of Components Composed of Strain and Work Hardening Viscoplastic Materials. Int. J. Solids Struct. 1980, 16, 793–806. [Google Scholar] [CrossRef]
  102. Ponter, A.R.S. Convexity and associated continuum properties of a class of constitutive relationships. J. De Mec. 1976, 15, 527–542. [Google Scholar]
  103. Bridgman, P.W. Flow and Fracture; Metals Plasticity, Technical Publication; Harvard University Press: Cambridge, MA, USA, 1944; p. 1782. [Google Scholar]
  104. Bingham, E.C. Fluidity and Plasticity; McGraw-Hill: New York, NY, USA, 1922. [Google Scholar]
  105. Hohenemser, K.; Prager, W. Ueber die Ansaetze der Mechanik isotroper Kontinua. Zeit Fuer Angew. Math. Und Mech. 1932, 12, 216–226. [Google Scholar] [CrossRef]
  106. Bailey, R.W. Note on the Softening of Strain Hardened Metals and its Relation to Creep. J. Inst. Met. 1926, 35, 27. [Google Scholar]
  107. Duffy, S.F. A Viscoplastic Constitutive Theory for Transversely Isotropic Metal Alloys. Ph.D. Thesis, University of Akron, Akron, OH, USA, 1987. [Google Scholar]
  108. Saleeb, A.; Arnold, S.; Castelli, M.; Wilt, T.; Graf, W. A General Hereditary Multi- Mechanism Based Deformation Model with Application to the Visco-elastoplastic Response of Titanium Alloys. Int. J. Plast. 2001, 17, 1305–1350. [Google Scholar] [CrossRef]
  109. Chaboche, J.-L.; Gaubert, A.; Kanouté, P.; Longuet, A.; Azzouz, F.; Mazière, M. Viscoplastic Constitutive Equations of Combustion Chamber Materials Including Cyclic Hardening and Dynamic Strain Aging. Int. J. Plast. 2013, 46, 1–22. [Google Scholar] [CrossRef]
  110. Prandtl, L. Spannungsverteilung in Plastischen Körpern. In Proceedings of the First International Congress on Applied Mechanics, Delft, The Netherlands, 22–26 April 1924; Technische Boekhandel en Drukkerij. Waltman, J., Jr., Ed.; Waltman: Delft, The Netherlands, 1925; pp. 43–45. [Google Scholar]
  111. Reuss, E. Berücksichtigung der Elastischen Formänderungen in der Plastizitätstheorie. Z. Für Angew. Math. Und Mech. 1930, 10, 266–274. [Google Scholar] [CrossRef]
  112. Kim, S.J.; Oden, J.T. Generalized Flow Potentials in Finite Elastoplasticity–II Examples. Int. J. Eng. Sci. 1985, 23, 515–530. [Google Scholar] [CrossRef]
  113. Bodner, S.R. Review of a Unified Elastic—Viscoplastic Theory. In Unified Constitutive Equations for Creep and Plasticity; Miller, A.K., Ed.; Springer: Dordrecht, The Netherlands, 1987; ISBN 978-94-010-8039-2. [Google Scholar]
  114. Cassenti, B.N. NASA CR-168191; Research and Development Program for the Development of Advanced Time-Temperature Dependent Constitutive Relationships, Volume I-Theoretical Discussion. NASA: Washington, DC, USA, 1983.
  115. Krieg, R.D.; Swearengen, J.C.; Rhode, R.B. A Physically-based Internal Variable Model for Rate-dependent Plasticity. In Proceedings of the ASHE/CSME PVP Conference, Montreal, QC, Canada, 25–29 June 1978; pp. 15–27. [Google Scholar]
  116. Hart, E.W. A Phenomenological Theory for Plastic Deformation of Polycrystalline Metals. Acta Metall. 1970, 18, 599–610. [Google Scholar] [CrossRef]
  117. Hart, E.W.; Solomon, H.D. Load Relaxation Studies of Polycrystalline High Purity Aluminum. Acta Metall. 1973, 21, 295–307. [Google Scholar] [CrossRef]
  118. Hart, E.W. Constitutive Relations for the Non-elastic Deformation of Metals. J. Eng. Mater. Technol. 1976, 98, 193. [Google Scholar] [CrossRef]
  119. Delph, T.J. A Comparative Study of Two State-Variable Constitutive Theories. ASME J. Mater. Technol. 1980, 102, 327–336. [Google Scholar] [CrossRef]
  120. Bouc, R. Forced Vibration of Mechanical Systems with Hysteresis. In Proceedings of the Fourth Conference on Nonlinear Oscillation, Prague Czechoslovakia, Czech Republic, 30 May 1967; p. 315. [Google Scholar]
  121. Bouc, R. Modèle Mathématique d’Hystérésis: Application aux Systèmes à un Degré de Liberté. Acustica 1971, 24, 16–25. [Google Scholar]
  122. Wen, Y.K. Method for Random Vibration of Hysteretic Systems. ASCE J. Eng. Mech. 1976, 102, 249–263. [Google Scholar] [CrossRef]
  123. Khoei, A.R. Chapter 7 Advanced Plasticity Models – Endochronic Plasticity Theory. In Computational Plasticity in Powder Forming Processes; Elsevier: Amsterdam, The Netherlands, 2005; ISBN 9780080529707. [Google Scholar]
  124. Valanis, K.C. Fundamental Consequences of a New Intrinsic Time Measure-Plasticity as a Limit of the Endochronic Theory. Arch. Mech. Stossowanej 1980, 32, 171–191. [Google Scholar]
  125. Wu, H.C.; Ho, C.C. An Investigation of Transient Creep by Means of Endochronic Viscoplasticity and Experiment. ASME J. Mater. Technol. 1995, 117, 260–268. [Google Scholar] [CrossRef]
  126. Coleman, B.D.; Gurtin, M.E. Thermodynamics with Internal State Variables. J. Chem. Phys. 1967, 47, 597–613. [Google Scholar] [CrossRef]
  127. Lubliner, J. On the Thermodynamic Foundations of Non-linear Solid Mechanics. Int. J. Non-Linear Mech. 1972, 7, 237–254. [Google Scholar] [CrossRef]
  128. Lubliner, J. On the Structure of the Rate Equations of Materials with Internal Variables. Acta Mech. 1973, 17, 109–119. [Google Scholar] [CrossRef]
  129. Lubliner, J. A Simple Theory of Plasticity. Int. J. Solids Struct. 1974, 10, 313–319. [Google Scholar] [CrossRef]
  130. Horstemeyer, M.F.; Gokhale, A.M. Modeling Stress State Dependent Damage Evolution in a Cast Al–Si–Mg aluminum Alloy. Theor. Appl. Fract. Mech. 2000, 33, 31–47. [Google Scholar] [CrossRef]
  131. Bammann, D.J.; Horstemeyer, M.F. Failure in Ductile Materials Using Finite Element Methods. In Structural Crashworthiness and Failure; CRC Press: Boca Raton, FL, USA, 1993; pp. 1–54. [Google Scholar]
  132. Rice, J.R. Inelastic Constitutive Relations for Solids: An Internal Variable Theory and Its Application to Metal Plasticity. J. Mech. Phys. Solids 1971, 19, 433–455. [Google Scholar] [CrossRef]
  133. Blum, W.; Eisenlohr, P. Dislocation Mechanics of Creep. Mater. Sci. Eng. 2009, 510–511, 7–13. [Google Scholar] [CrossRef]
  134. Hill, R. Elastic Properties of Reinforced Solids: Some Theoretical Principles. J. Mech. Phys. Solids 1963, 11, 357–372. [Google Scholar] [CrossRef]
  135. Terada, K.; Kikuchi, N. Nonlinear Homogenization Method for Practical Applications. Am. Soc. Mech. Eng. Appl. Mech. Div. 1995, 212, 1–16. [Google Scholar]
  136. Ghosh, S.; Lee, K.; Moorthy, S. Multiple Scale Analysis of Heterogeneous Elastic Structures Using Homogenization Theory and Voronoi Cell Finite Element Method. Int. J. Solids Struct. 1995, 32, 27–62. [Google Scholar] [CrossRef]
  137. Ghosh, S.; Lee, K.; Moorthy, S. Two Scale Analysis of Heterogeneous Elastic-Plastic Materials with Asymptotic Homogenization and Voronoi cell Finite Element Model. Computer. Methods Appl. Mech. Eng. 1996, 132, 63–116. [Google Scholar] [CrossRef]
  138. Smit, R.J.; Brekelmans, W.M.; Meijer, H.E. Prediction of the Mechanical Behavior of Nonlinear Heterogeneous Systems by Multi-Level Finite Element Modeling. Comput. Methods Appl. Mech. Eng. 1998, 155, 181–192. [Google Scholar] [CrossRef]
  139. Feyel, F.; Chaboche, J.L. FE2 Multiscale Approach for Modelling the Elastoviscoplastic Behaviour of Long Fibre SiC/Ti Composite Materials. Comput. Methods Appl. Mech. Eng. 2000, 183, 309–330. [Google Scholar] [CrossRef]
  140. Holm, E.A.; Battaile, C.C.; Fang, H.E.; E Buchheit, T.; Wellman, G.W. Making the Connection between Microstructure and Mechanics; U.S. Department of Energy Office of Scientific and Technical Information: Oak Ridge, TN, USA, 2003.
  141. Horner, D.A.; Peters, J.F.; Carrillo, A. Large Scale Discrete Element Modeling of Vehicle-Soil Interaction. J. Eng. Mech. 2001, 127, 1027–1032. [Google Scholar] [CrossRef]
  142. Iwashita, K.; Oda, M. Micro-Deformation Mechanism of Shear Banding Process Based on Modified Distinct Element Method. Powder Technol. 2000, 109, 192–205. [Google Scholar] [CrossRef]
  143. Tordesillas, A.; Walsh, S.D.C.; Gardiner, B.S. Bridging the Length Scales: Micromechanics of Granular Media. BIT Numer. Math. 2004, 44, 539–556. [Google Scholar] [CrossRef]
  144. Mindlin, R.D. Microstructure in Linear Elasticity; Columbia University New York Department of Civil Engineering and Engineering Mechanics: New York, NY, USA, 1963. [Google Scholar]
  145. Liu, W.K.; Qian, D.; Gonella, S.; Li, S.; Chen, W.; Chirputkar, S. Multiscale Methods for Mechanical Science of Complex Materials: Bridging from Quantum to Stochastic Multiresolution Continuum. Int. J. Numer. Methods Eng. 2010, 83, 1039–1080. [Google Scholar] [CrossRef]
  146. Gao, H.; Huang, Y.; Nix, W.D.; Hutchinson, J.W. Mechanism-Based Strain Gradient Plasticity—I. Theory. J. Mech. Phys. Solids 1999, 47, 1239–1263. [Google Scholar] [CrossRef]
  147. Holm, E. A Overview: Incorporating Microstructural Scale Damage into Continuum Models for Performance of Structural Materials; Computational Materials Science and Engineering Sandia National Laboratories: Albuquerque, NM, USA, 2012. [Google Scholar]
  148. Fish, J.; Wagner, J.; Keten, S. Mesoscopic and Multiscale Modelling in Materials. Nat. Mater. 2021, 20, 774–786. [Google Scholar] [CrossRef] [PubMed]
  149. Vladmir Buljak Gianluca Ranzi Constitutive Modeling of Engineering Material; Academic Press: Cambridge, MA, USA, 2021.
Figure 1. Classification of Inelastic Constitutive Models.
Figure 1. Classification of Inelastic Constitutive Models.
Materials 16 03564 g001
Figure 2. Generic deformation mechanism map.
Figure 2. Generic deformation mechanism map.
Materials 16 03564 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Baig, M.; Owusu-Danquah, J.; Campbell, A.A.; Duffy, S.F. Review: Inelastic Constitutive Modeling: Polycrystalline Materials. Materials 2023, 16, 3564. https://doi.org/10.3390/ma16093564

AMA Style

Baig M, Owusu-Danquah J, Campbell AA, Duffy SF. Review: Inelastic Constitutive Modeling: Polycrystalline Materials. Materials. 2023; 16(9):3564. https://doi.org/10.3390/ma16093564

Chicago/Turabian Style

Baig, Mirza, Josiah Owusu-Danquah, Anne A. Campbell, and Stephen F. Duffy. 2023. "Review: Inelastic Constitutive Modeling: Polycrystalline Materials" Materials 16, no. 9: 3564. https://doi.org/10.3390/ma16093564

APA Style

Baig, M., Owusu-Danquah, J., Campbell, A. A., & Duffy, S. F. (2023). Review: Inelastic Constitutive Modeling: Polycrystalline Materials. Materials, 16(9), 3564. https://doi.org/10.3390/ma16093564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop