# State storage

 Many of the constitutive models in Abaqus require tensors to be stored to define the state at a material calculation point. Such “material state tensors” are stored as their components in a local, orthonormal, system at the material calculation point. The orientation of that system with respect to the global $(X,Y,Z)$ spatial system is stored as a rotation from the global axis system. The purpose of this section is to define the manner in which such tensors are stored and updated.
 Related Topics In Other Guides Conventions

ProductsAbaqus/StandardAbaqus/Explicit

Three types of local basis are used in Abaqus for material calculations. For isotropic materials in continuum elements the global, spatial, $(X,Y,Z)$ system is used—the material basis is fixed in time. For isotropic materials in structural surface elements (shells and membranes) the local system is defined by the standard Abaqus convention described in Conventions, and for beams and trusses it is defined with the 1-direction along the axis of the member and the 2- and 3-directions in material directions in the cross-section. Thus, with isotropic materials the material basis is always the same as the element basis, although for structural elements the material basis changes with time. For anisotropic materials the material basis must be defined by the user and rotates with the average rigid body spin of the material. In this case the material basis and the element basis are not the same.

We refer to this local material basis at time t as $eiM|t$, where the superscript $M$ indicates that the basis is associated with material calculations and $|t$ means that the basis is taken at time t. In this section Latin subscripts (like the $i$ above) take the range 1–3, while Greek subscripts will take the range 1–2.

Any tensor associated with the material's state, $a$, say (such as the stress tensor $σ$), is stored in terms of its components along the material basis:

$a=ai⁢j⁢eiM⁢ejM.$

The increment from time t to time $t+Δ⁢t$ of local motion at the material calculation point is defined by the incremental deformation gradient,

$Δ⁢F⁢=def⁢∂⁡xt+Δ⁢t∂⁡xt.$

This matrix is calculated from the gradient interpolator of the finite element and the coordinates of the element's nodes at times t and $t+Δ⁢t$.

The polar decomposition of $Δ⁢F$ is

$Δ⁢F=Δ⁢V⋅Δ⁢R,$

where $Δ⁢R$ is the average rigid body rotation at the material point and $Δ⁢V$ is a pure stretch matrix:

$Δ⁢V=∑I=13Δ⁢λI⁢nI⁢nI$

(here $Δ⁢λI$ is a principal stretch ratio and $nI$ is a principal stretch direction).

During an increment any material state tensor changes according to

$a|t+Δ⁢t=Δ⁢a∇+Δ⁢R⋅a|t⋅Δ⁢RT,$

where $Δ⁢a∇$ is the change in $a$ caused by constitutive behavior and $Δ⁢R$ is the average incremental rigid body rotation of the material. Since material tensors are written in terms of their components in the material basis system, this update is computed as

$(ai⁢j)|t+Δ⁢t=Δ⁢ai⁢j∇+Δ⁢Ri⁢k⁢(ak⁢l)|t⁢Δ⁢Rj⁢l.$

It is, therefore, necessary to project $Δ⁢R$ onto the material basis systems at the start and end of the increment to define the update of material tensor components:

$Δ⁢Ri⁢j=eiM|t+Δ⁢t⋅Δ⁢R⋅ejM|t.$

For isotropic materials the $eiM$ have been chosen for geometric convenience only, so the $Δ⁢Ri⁢j$ are quite general.

For anisotropic materials the material basis system, $eiM$, rotates with the average rigid body rotation of the material, $Δ⁢R$, and so is updated by

$eiM|t+Δ⁢t=Δ⁢R⋅eiM|t.$

In this case we see that

$Δ⁢Ri⁢j=δi⁢j,$

and so the update of a material tensor's components simplifies to

$(ai⁢j)|t+Δ⁢t=Δ⁢ai⁢j∇+(ai⁢j)|t.$

However, since in this case the material basis system is not the same as the element basis system, $eiE$ (which is chosen for geometric convenience for element calculations), transformations must be done to change basis system. Specifically, at the start of the material calculation routines, the strain increments are rotated from the element basis into the material basis:

$Δ⁢εi⁢jM=Δ⁢εk⁢lE⁢[eiM⋅ekE]⁢[ejM⋅elE].$

Likewise, at the end of the material calculation routines, the stress increments are rotated back to the element basis for integration into the discretized approximation to the equilibrium equations:

$Δ⁢σi⁢jE=Δ⁢σk⁢lM⁢[eiE⋅ekM]⁢[ejE⋅elM].$

In addition, the material stiffness matrix,

$Di⁢j⁢k⁢l⁢=def⁢∂⁡σi⁢j∂⁡εk⁢l$

must also be rotated from the material basis to the element basis:

$Di⁢j⁢k⁢lE=Dm⁢n⁢p⁢qM⁢[eiE⋅emM]⁢[ejE⋅enM]⁢[ekE⋅epM]⁢[elE⋅eqM].$

For a shell or membrane only two-dimensional rotations are required—for example,

$Δ⁢εα⁢βM=Δ⁢εγ⁢δE⁢[eαM⋅eγE]⁢[eβM⋅eδE],$

since $e3M≡e3E$ because both are unit vectors along the normal to the surface.