Internal state variable

K set of referential internal state variables I velocity gradient tensor... 

The theory is initially presented in the context of small deformations in Section 5.2. A set of internal state variables are introduced as primitive quantities, collectively represented by the symbol k. Qualitative concepts of inelastic deformation are rendered into precise mathematical statements regarding an elastic range bounded by an elastic limit surface, a stress-strain relation, and an evolution equation for the internal state variables. While these qualitative ideas lead in a natural way to the formulation of an elastic limit surface in strain space, an elastic limit surface in stress space arises as a consequence. An assumption that the external work done in small closed cycles of deformation should be nonnegative leads to the existence of an elastic potential and a normality condition. 

Specific applications of the theory are not considered in this chapter. Only one example, that of small deformation classical plasticity, is worked out in Section 5.3. The set of internal state variables k is taken to be comprised of... 

As with any constitutive theory, the particular forms of the constitutive functions must be constructed, and their parameters (material properties) must be evaluated for the particular materials whose response is to be predicted. In principle, they are to be evaluated from experimental data. Even when experimental data are available, it is often difficult to determine the functional forms of the constitutive functions, because data may be sparse or unavailable in important portions of the parameter space of interest. Micromechanical models of material deformation may be helpful in suggesting functional forms. Internal state variables are particularly useful in this regard, since they may often be connected directly to averages of micromechanical quantities. Often, forms of the constitutive functions are chosen for their mathematical or computational simplicity. When deformations are large, extrapolation of functions borrowed from small deformation theories can produce surprising and sometimes unfortunate results, due to the strong nonlinearities inherent in the kinematics of large deformations. The construction of adequate constitutive functions and their evaluation for particular... 

Motivated by the qualitative observations made above, a set of internal state variables deseribing the internal strueture of the material will be intro-dueed ab initio, denoted eolleetively by k. Their physieal meaning or preeise properties need not be established at this point, and they may inelude sealar, veetor, or tensor quantities. The following eonstitutive assumptions are now made ... 

Inelastic Loading. The strain lies on the elastic limit surface = 0, and the tangent to the strain history points in a direction outward from the elastic limit surface > 0. The material is said to be undergoing inelastic loading, and k is assumed to be a function of the strain s, the internal state variables k, and the strain rate k... 

This implies that, at constant k, the line integral of the differential form s de, parametrized by time t, taken over the closed curve h) zero. This is the integrability condition for the existence of a scalar function tj/ e) such that s = d j//de (see, e.g., Courant and John [13], Vol. 2, 1.10). This holds for an elastic closed cycle at any constant values of the internal state variables k. Therefore, in general, there exists a function ij/... 

The normality conditions (5.56) and (5.57) have essentially the same forms as those derived by Casey and Naghdi [1], [2], [3], but the interpretation is very different. In the present theory, it is clear that the inelastic strain rate e is always normal to the elastic limit surface in stress space. When applied to plasticity, e is the plastic strain rate, which may now be denoted e", and this is always normal to the elastic limit surface, which may now be called the yield surface. Naghdi et al. by contrast, took the internal state variables k to be comprised of the plastic strain e and a scalar hardening parameter k. In their theory, consequently, the plastic strain rate e , being contained in k in (5.57), is not itself normal to the yield surface. This confusion produces quite different results. 

The remainder of this section will be concerned with a particular case in which normality conditions hold. The constitutive equation for the internal state variables (5.11) involves the constitutive function a, and the normality conditions (5.56) and (5.57) involve an unknown scalar factor y. In some circumstances, a may be eliminated and y may be evaluated by using the consistency condition. These circumstances arise if b is nonsingular so that the normality condition in strain space (5.56j) may be solved for k... 

In this section, the general inelastic theory of Section 5.2 will be specialized to a simple phenomenological theory of plasticity. The inelastic strain rate tensor e may be identified with the plastic strain rate tensor e . In order to include isotropic and kinematic hardening, the set of internal state variables, denoted collectively by k in the previous theory, is reduced to the set (k, a) where k is a scalar representing isotropic hardening and a is a symmetric second-order tensor representing kinematic hardening. The elastic limit condition in stress space (5.25), now called a yield condition, becomes... 

Prager s rule of kinematic hardening is expressed by a = ce where c is a constant. Generalizing these concepts, the evolution equations for the internal state variables will be taken in the form... 

The choice (5.77) for the evolution equation for the plastic strain sets the evolution equations for the internal state variables (5.78) into the form (5.11) required for continuity. The consistency condition in the stress space description may be obtained by differentiating (5.73), or directly by expanding (5.29)... 

In order to consider the inelastic stress rate relation (5.111), some assumptions must be made about the properties of the set of internal state variables k. With the back stress discussed in Section 5.3 in mind, it will be assumed that k represents a single second-order tensor which is indifferent, i.e., it transforms under (A.50) like the Cauchy stress or the Almansi strain. Like the stress, k is not indifferent, but the Jaumann rate of k, defined in a manner analogous to (A.69), is. With these assumptions, precisely the same arguments... 

In Section 5.2 the set of internal state variables k was introduced. In the referential theory, a similar set of referential internal state variables K will be introduced in the same way without further physical identification at this stage. It will merely be assumed that each member of the set K is invariant under the coordinate transformation (A.50) representing a rigid rotation and translation of the coordinate frame. 

It may first be noted that the referential symmetric Piola-Kirchhoff stress tensor S and the spatial Cauchy stress tensor s are related by (A.39). Again with the back stress in mind, it will be assumed in this section that the set of internal state variables is comprised of a single second-order tensor whose referential and spatial forms are related by a similar equation, i.e., by... 

P.S. Follansbee and U.F. Kocks, A Constitutive Description of the Deformation of Copper Based on the Use of the Mechanical Threshold Stress as an Internal State Variable, Acta Metall. 36, 81-93 (1988). 

Everything in Swarm is an object with three main characteristics Name, Data and Rules. An object s Name consists of an ID that is used to send messages to the object, a type and a module name. An object s Data consists of whatever local data (i.e. internal state variables) the user wants an agent to possess. The Rules are functions to handle any messages that are sent to the object. The basic unit of Swarm is a swarm a collection of objects with a schedule of event over those objects. Swarm also supplies the user with an interface and analysis tools. 

The inelastic strain as an internal state variable is obtained from a linear evolution equation formulated with respect to the intermediate configuration ... 

We suppose that exists a non-negative kinetic energy Kt(/zt, v t), associated with each time-rate of change such that 0) = 0 and 0. In the absence of rp, the measure is sooner termed internal (state) variable and ruled by a first order evolution equation instead of a balance equation. 

In the next sections, the multiscale modeling methods are presented from the different disciplines perspectives. Clearly one could argue that overlaps occur, but the idea here is to present the multiscale methods from the paradigm from which they started. For example, the solid mechanics internal state variable theory includes mathematics, materials science, and numerical methods. However, it clearly started from a solid mechanics perspective and the starting point for mathematics, materials science, and numerical methods has led to other different multiscale methods. 

One effective hierarchical method for multiscale bridging is the use of thermodynamically constrained internal state variables (IS Vs) that can be physically based upon microstructure-property relations. It is a top-down approach, meaning the IS Vs exist at the macroscale but reach down to various subscales to receive pertinent information. The ISV theory owes much of its development to the state variable thermodynamics constructed by Helmholtz [4] and Maxwell [5]. The notion of ISV was introduced into thermodynamics by Onsager [6, 7] and was applied to continuum mechanics by Eckart [8, 9]. 

In terms of cyclic plasticity, Shenoy et al. [139, 140], Wang et al. [137], and McDowell [141] performed hierarchical multiscale modeling of Ni-based superalloys employing internal state variable theory. Fan et al. [142] performed a hierarchical multiscale modeling strategy for three length scales. 

Next to metals, probably the synthetic polymer-based composites have been modeled most by hierarchical multiscale methods. Different multiscale formulations have been approached top-down internal state variable approaches, self-consistent (or homogenization) theories, and nanoscale quantum-molecular scale methods. 

Allen [240], Costanzo et al. [241], and Krajcinovic [242-244] that thermodynamics of appropriate internal state variables for damage in composites under mechanical loads were addressed. Similar to the metal plasticity theoreticians, these researchers employed Coleman and Gurtin s [11] thermodynamic framework to determine the kinetic equations. However, unlike the metal internal state variable community that could quantify the evolution equations for dislocations and damage, these polymer-based composites theoreticians did not propose evolutionary rate equations, but just damage state equations. 

In order to start the multiscale modeling, internal state variables were adopted to reflect void/crack nucleation, void growth, and void coalescence from the casting microstructural features (porosity and particles) under different temperatures, strain rates, and deformation paths [115, 116, 221, 283]. Furthermore, internal state variables were used to reflect the dislocation density evolution that affects the work hardening rate and, thus, stress state under different temperatures and strain rates [25, 283-285]. In order to determine the pertinent effects of the microstructural features to be admitted into the internal state variable theory, several different length scale analyses were performed. Once the pertinent microstructural features were determined and included in the macroscale internal state variable model, notch tests [216, 286] and control arm tests were performed to validate the model s precision. After the validation process, optimization studies were performed to reduce the weight of the control arm [287-289]. 

In order to quantify the structure-property relations at each scale, a multiscale hierarchy of numerical simulations was performed, coupled with experiments, to determine the internal state variable equations of macroscale plasticity and damage... 

Now that the top-down internal state variable theory was established, the bottom-up simulations and experiments were required. At the atomic scale (nanometers), simulations were performed using Modified Embedded Atom Method, (MEAM) Baskes [176], potentials based upon interfacial atomistics of Baskes et al. [177] to determine the conditions when silicon fracture would occur versus silicon-interface debonding [156]. Atomistic simulations showed that a material with a pristine interface would incur interface debonding before silicon fracture. However, if a sufficient number of defects were present within the silicon, it would fracture before the interface would debond. Microstructural analysis of larger scale interrupted strain tests under tension revealed that both silicon fracture and debonding of the silicon-aluminum interface in the eutectic region would occur [290, 291]. 

---