.. _stress_measures:

###############
Stress Measures
###############

*************
Cauchy Stress
*************

In this chapter we discuss the quantification of force intensity, or stress, within a body undergoing potentially large amounts of deformation.  We begin with the **Cauchy stress** tensor :math:`\mathbf{T}` and note that, provided we associate this object with the spatial configuration, this object can be interpreted exactly as in the infinitesimal case outlined in :numref:`linear_elastic_ivp_bvp`.  In the current notational framework, we interpret the components of :math:`\mathbf{T}`, denoted as :math:`T_{ij}`, which represent forces per unit areas in the spatial configuration at a given spatial point :math:`\mathbf{x} \in \varphi_t (\Omega)`.

It will be necessary in our description to consider related measures of stress defined in terms the reference and rotated configurations. To motivate this discussion, we reconsider the concept of traction discussed previously in the context of the infinitesimal elastic system. Recall that given a plane passing through the point of interest :math:`\mathbf{x}`, the traction, or force per unit area acting on this plane, is given by the formula

.. math::
   :label: stress_measure:eq:01

   t_i = T_{ij} n_j ,

where :math:`n_j` is the unit normal vector to the plane in question.

****************
Nanson's Formula
****************

We consider two differential vectors, :math:`\mathrm{d}\mathbf{r}_1` and :math:`\mathrm{d}\mathbf{r}_2`, as illustrated in :numref:`stress_measure-fig-nanson`.  We assume that :math:`\mathrm{d}\mathbf{r}_1` and :math:`\mathrm{d}\mathbf{r}_2` are linearly independent and that both have spatial point :math:`\mathbf{x}` as their base point.  We further assume that their orientations are such that the following relation from basic geometry holds:

.. math::
   :label: stress_measure:eq:02

   \mathrm{d}\mathbf{r}_1 \times \mathrm{d}\mathbf{r}_2 = \mathbf{n} \mathrm{d}a,

where :math:`\mathrm{d}a` is the (differential) area of the parallelogram encompassed by :math:`\mathrm{d}\mathbf{r}_1` and :math:`\mathrm{d}\mathbf{r}_2`.

.. _stress_measure-fig-nanson:

.. figure:: ../_static/figures/stress_measures-fig01.png
   :align: center
   :scale: 75 %

   Notation for derivation of Nanson's formula.

As discussed in :numref:`deformation_measures` (see :eq:`deformation_measure:eq:05`), we can think of the differential vectors :math:`\mathrm{d} \mathbf{r}_1` and :math:`\mathrm{d} \mathbf{r}_2` as the current positions of reference differential vectors :math:`\mathrm{d} \mathbf{R}_1` and :math:`\mathrm{d} \mathbf{R}_2`, which are based at :math:`\mathbf{X}=\varphi^{-1}_t (\mathbf{x})`.  In index notation, we can relate these two sets of differential vectors using the deformation gradient via
 
.. math::
   :label: stress_measure:eq:03

   (\mathrm{d}\mathbf{r}_1)_i = F_{iI} (\mathrm{d}\mathbf{R}_1)_I ,

and

.. math::
   :label: stress_measure:eq:04

   (\mathrm{d}\mathbf{r}_2)_i = F_{iI} (\mathrm{d}\mathbf{R}_2)_I .

We now seek to re-express :eq:`stress_measure:eq:02` in terms of reference quantities. Working in index notation, we write

.. math::
   :label: stress_measure:eq:05

   \begin{aligned}
   n_i \mathrm{d}a & = e_{ijk} F_{jJ} (\mathrm{d}\mathbf{R}_1)_J F_{kK} 
                       (\mathrm{d}\mathbf{R}_2)_K \\
                & = e_{ljk} \delta_{li} F_{jJ} (\mathrm{d}\mathbf{R}_1)_J F_{kK} 
                       (\mathrm{d}\mathbf{R}_2)_K \\
                & = e_{ljk} F_{lL} F^{-1}_{Li} F_{jJ} (\mathrm{d}\mathbf{R}_1)_J F_{kK} 
                       (\mathrm{d}\mathbf{R}_2)_K
   \end{aligned}

We extract and manipulate a particular product in the last line of :eq:`stress_measure:eq:05`, namely :math:`e_{ljk} F_{lL} F_{jJ} F_{kK}`. One can show by a case-by-case examination that the following relation holds:

.. math::
   :label: stress_measure:eq:06

   e_{ljk} F_{lL} F_{jJ} F_{kK} = e_{LJK} e_{ljk} F_{l1} F_{j2} F_{k3} .

Recall from :numref:`deformation_measures`, :eq:`deformation_measure:eq:11` that :math:`J=\text{ det}(\mathbf{F})` has the following representation in index notation:

.. math::
   :label: stress_measure:eq:07

   J = \text{ det}(\mathbf{F}) = e_{ljk} F_{l1} F_{j2} F_{k3}

Combination of :eq:`stress_measure:eq:05`, :eq:`stress_measure:eq:06`, and :eq:`stress_measure:eq:07` yields the following result:

.. math::
   :label: stress_measure:eq:08

   \begin{aligned}
   n_i \mathrm{d}a & = J e_{LJK} F^{-1}_{Li} (\mathrm{d}\mathbf{R}_1)_J 
                       (\mathrm{d}\mathbf{R}_2)_K \\
                   & = J F^{-1}_{Li} m_L \mathrm{d}A.
   \end{aligned}

In :eq:`stress_measure:eq:08`, :math:`\mathrm{d}A` is the differential reference area spanned by :math:`\mathrm{d}\mathbf{R}_1` and :math:`\mathrm{d}\mathbf{R}_2`, and :math:`\mathbf{m}` is the reference unit normal to this area.

In direct notation, we express this result as

.. math::
   :label: stress_measure:eq:09

   \mathbf{n}\mathrm{d}a = J \mathbf{F}^{-T} \mathbf{m} \mathrm{d}A .

:eq:`stress_measure:eq:09` is referred to as Nanson's formula and it is important, among other reasons, because it provides the appropriate change-of-variables formula for surface integrals in the reference and current configurations.

***************************************
First and Second Piola-Kirchhoff Stress
***************************************

We want to compute a differential force, which is the product of the traction acting on our plane at :math:`\mathbf{x}` and the differential area under consideration. Denoting this differential force by :math:`\mathrm{d}\mathbf{f}`, we write

.. math::
   :label: stress_measure:eq:10

   \mathrm{d}\mathbf{f} = \mathbf{t}\mathrm{d}a 
                     = \mathbf{T}\mathbf{n}\mathrm{d}a 
                     = J \mathbf{T}\mathbf{F}^{-T} \mathbf{m} \mathrm{d} A .

In examining :eq:`stress_measure:eq:10`, we find that the following definition is useful:

.. math::
   :label: stress_measure:eq:11

   \mathbf{P}(\mathbf{X}) = J(\mathbf{X}) \mathbf{T} \left( \varphi_t(\mathbf{X}) \right) 
                         \mathbf{F}^{-T} \left( \varphi_t(\mathbf{X}) \right) .

This allows us to write

.. math::
   :label: stress_measure:eq:12

   \mathrm{d}\mathbf{f} = \mathbf{P} \mathbf{m} \mathrm{d}A .

In :eq:`stress_measure:eq:12`, the product :math:`\mathbf{P}\mathbf{m}` represents a traction, the current force, :math:`\mathrm{d}\mathbf{f}`, divided by the reference area, :math:`\mathrm{d}A`.  The tensor :math:`\mathbf{P}` is called the (**First**) **Piola-Kirchhoff Stress** and :math:`\mathbf{P}\mathbf{m}` is called the Piola Traction.  Similar to the Piola Traction, the First Piola-Kirchhoff Stress measures stress in terms of forces in the current configuration and areas in the reference configuration.  The one-dimensional manifestation of this stress measure is the engineering stress, :math:`\sigma_E`, originally defined in :eq:`nonlinear_behavior:eq:03`.

It is worthy to note that :math:`\mathbf{P}` is neither a pure spatial nor a pure reference object. A reference object for stress can be constructed by performing a pull-back of the spatial Cauchy stress tensor :math:`\mathbf{T}` to the reference configuration:

.. math::
   :label: stress_measure:eq:13

   \begin{array} {l l}
   \mathbf{S}( \mathbf{X} ) & = J \mathbf{F}^{-1} \left( \varphi_t (\mathbf{X}) \right) 
                             \mathbf{T} \left( \varphi_t (\mathbf{X}) \right) 
                             \mathbf{F}^{-T} \left( \varphi_t (\mathbf{X}) \right) \\
                      {} & = \mathbf{F}^{-1} \left( \varphi_t (\mathbf{X}) \right) 
                             \mathbf{P} (\mathbf{X}) .
   \end{array}


:math:`\mathbf{S}` is called the **Second Piola-Kirchhoff Stress** and it is purely a reference object.  We note, in particular, that :math:`\mathbf{S}` is a symmetric tensor, while :math:`\mathbf{P}` is not symmetric in general.

This same concept of pull-back can be employed to define a stress tensor in the rotated configuration, which we shall denote as :math:`\pmb{\mathit{T}}`. This rotated tensor is defined as

.. math::
   :label: stress_measure:eq:14

   \pmb{\mathit{T}} = \mathbf{R}^{T} \left( \varphi_t (\mathbf{X}) \right) \mathbf{T} 
                   \left( \varphi_t (\mathbf{X}) \right) \mathbf{R} 
                   \left( \varphi_t (\mathbf{X}) \right) .


As was the case with the rotated configuration quantities introduced in :numref:`rates_of_deformation`, this definition will be of particular importance in the later examination of frame indifference.
