Stress rates

We define the rate of Kirchhoff stress for use in constitutive equations as an appropriate stress measure for stress-sensitive materials.

The following topics are discussed:

Related Topics
In Other Guides


Many of the materials we wish to model with Abaqus are history dependent, and it is common for the constitutive equations to appear in rate form. In Stress measures it was suggested that an appropriate stress measure for stress-sensitive materials (such as yielding materials) is the Kirchhoff stress. We, therefore, need to define the rate of Kirchhoff stress for use in the constitutive equations. This definition is not simply the material time rate of Kirchhoff stress, because the Kirchhoff stress components are associated with spatial directions in the current configuration (recall that the Kirchhoff stress is Jσ, where J is the volume change from the reference configuration and σ is the Cauchy stress, defined by t=σn, where t and n are vectors in the current configuration).

To illustrate the issue, consider a uniaxial tension specimen under constant axial force P, lying along the x-axis at time t1 and rotated—with the axial force held constant—to lie along the y-axis at time t2 (see Figure 1).

Figure 1. Rotated specimen.

Write the stress components on the global (1,2,3) rectangular Cartesian basis. At time t1, σ11=P/A, and all other σij=0, while at time t2, σ22=P/A, and all other σij=0. Obviously during t1t2, dσ110 and dσ220, but equally clearly this rate of change of stress has nothing to do with the constitutive response of the material making up the bar. (A materially based stress, such as the second Piola-Kirchhoff stress, would stay constant during the above rotation, because its components are associated with a material basis.) The problem, then, is that the components of σ or τ are associated with current directions in space and, therefore, dτ and dσ will be nonzero if there is pure rigid body rotation, even though from a constitutive point of view the material is unchanged. Thus, we must divide the increment of σ or τ into two parts—one attributable to rigid body motion only and a remainder that is then, presumably, associated with the rate form of the stress-strain law.

We can derive a simple result for this purpose for any matrix whose components are associated with spatial directions. At some time t imagine attaching to a material point a set of base vectors, eα, α=1,2,3. These vectors cannot stretch but are defined to spin with the same spin as the material. Recall that the spatial gradient of the material particle velocity at a point, v/x, was decomposed into a rate of deformation and a spin,

D=12(vx+[vx]T) and W=12(vx-[vx]T).

One of the concepts of the motion of the base vectors eα in Abaqus is that


Another concept of the motion of the base vectors eα used in Abaqus is


where Ω=R˙RT. Here R is the rigid body rotation in the polar decomposition of the deformation gradient F. The differences between these two concepts are significant only if finite rotation of a material point is accompanied by finite shear.

Now consider any matrix T based on the current configuration: we can write it in terms of its components in the eα directions:


Taking the time derivative then gives


The second and third terms are the rate of T caused by the rigid body spin, so the first term is that part of T caused by other effects (in the case of stress, the rate associated with the constitutive response), called the corotational rate of T. From the definitions of eα as rigid base vectors that can be considered to spin with either W or Ω, we can write two corotational rates of T as




where TJ and TG are called the Jaumann and Green-Naghdi rates, respectively.

We, thus, have the total rate of any matrix associated with spatial directions in the current configuration as the sum of the corotational rate of the matrix and a rate caused purely by the local spin or rigid body rotation. For example, the Jaumann rate of change of Kirchhoff stress can be written as


We are assuming that the constitutive theory will define dJdt(Jσ), the corotational stress rate per reference volume, in terms of the rate of deformation and past history, so this equation provides a convenient link between that material model and the overall change in “true” (Cauchy) stress (which is the stress measure defined directly from the equilibrium equations). In Mechanical Constitutive Theories where the constitutive models in Abaqus are discussed, “stress rate” per reference volume will mean ddt(Jσ), the corotational rate of Kirchhoff stress, which is the stress measure work conjugate to the rate of deformation.