Strain Invariants and Their Engineering Applications

In the field of continuum mechanics and solid mechanics, the strain tensor $\boldsymbol{\varepsilon}$ serves as the fundamental descriptor of a body's deformation state. However, a significant challenge arises when working with tensor components: their numerical values are inherently dependent on the chosen coordinate system. If an observer rotates their frame of reference, the individual components of the strain tensor will change, even though the physical deformation of the material remains identical.

To overcome this coordinate dependency, we utilize the concept of strain invariants. Invariants are scalar-valued polynomial functions of the strain tensor components that remain constant under any orthogonal transformation (rotation). By focusing on these invariants, engineers and physicists can describe the state of deformation through objective and intrinsic properties that are independent of the observer's orientation.

Mathematical Definitions and Physical Significance

For a symmetric third-order strain tensor $\boldsymbol{\varepsilon}$, there are three primary invariants, typically denoted as $I_1, I_2,$ and $I_3$. These are closely related to the principal strains ($\varepsilon_1, \varepsilon_2, \varepsilon_3$), which are the eigenvalues of the tensor.

1. The First Invariant ($I_1$)

The first invariant is defined as the trace of the strain tensor, representing the sum of its diagonal components or the sum of the principal strains:
$$ I_1 = \text{tr}(\boldsymbol{\varepsilon}) = \varepsilon_{11} + \varepsilon_{22} + \varepsilon_{33} = \varepsilon_1 + \varepsilon_2 + \varepsilon_3 $$
Physical Significance: $I_1$ is directly associated with the volumetric strain. It quantifies the change in volume of a material element. In the study of incompressible materials (such as many rubbers or liquids), the volume remains constant during deformation, which mathematically implies that $I_1 = 0$.

2. The Second Invariant ($I_2$)

The second invariant is the sum of the principal strains taken two at a time, or the sum of the second-order principal minors of the tensor:
$$ I_2 = \varepsilon_1\varepsilon_2 + \varepsilon_2\varepsilon_3 + \varepsilon_3\varepsilon_1 $$
Physical Significance: $I_2$ is intimately linked to the shape change or distortion of the material. While $I_1$ tracks volume, $I_2$ helps characterize how the material is being stretched or sheared without necessarily changing its total volume. It is a critical component in constructing energy density functions for deviatoric (shape-changing) behavior.

3. The Third Invariant ($I_3$)

The third invariant is the determinant of the strain tensor, which is the product of the principal strains:
$$ I_3 = \det(\boldsymbol{\varepsilon}) = \varepsilon_1 \varepsilon_2 \varepsilon_3 $$
Physical Significance: $I_3$ provides information regarding the ratio of the deformed volume to the original volume. It serves as a higher-order measure of the volumetric state of the material.

Role in Constitutive Modeling

One of the most vital applications of strain invariants is in the formulation of constitutive equations—the mathematical models that describe how a specific material responds to applied loads.

The Principle of Objectivity

To ensure that a physical model is realistic, it must satisfy the principle of material frame indifference (or objectivity). This principle dictates that the constitutive response of a material should not change simply because the coordinate system is rotated. By expressing constitutive laws as functions of strain invariants rather than raw tensor components, engineers automatically satisfy this requirement.

Hyperelasticity and Energy Functions

For isotropic hyperelastic materials, the strain energy density function $W$ (the energy stored per unit volume) is expressed solely as a function of the invariants:
$$ W = W(I_1, I_2, I_3) $$
This approach allows for the elegant modeling of complex, non-linear materials. For instance, the widely used Mooney-Rivlin or Ogden models for elastomers are built upon these invariant relationships. Once the energy function is defined, the resulting stress tensor $\boldsymbol{\sigma}$ can be derived through differentiation:
$$ \boldsymbol{\sigma} = \frac{\partial W}{\partial \boldsymbol{\varepsilon}} = \sum_{i=1}^{3} \frac{\partial W}{\partial I_i} \frac{\partial I_i}{\partial \boldsymbol{\varepsilon}} $$
This mathematical framework simplifies the transition from theoretical energy states to practical stress analysis in complex engineering components.

Engineering Applications: Yielding and Plasticity

In the realm of metal plasticity and structural integrity, strain invariants are indispensable for predicting when a material will transition from elastic to plastic behavior.

1. Distortion Energy and the von Mises Criterion

In many metals, plastic deformation is driven not by changes in volume, but by changes in shape. This is captured by the concept of distortion energy density ($U_d$). For an isotropic linear elastic material, the distortion energy is related to the invariants as follows:
$$ U_d = \frac{1}{2} \mu \left( I_2 - \frac{1}{3}I_1^2 \right) $$
where $\mu$ is the shear modulus. This concept forms the theoretical backbone of the von Mises yield criterion. The criterion posits that a material begins to yield when the distortion energy reaches a critical threshold. This is the standard method used in Finite Element Analysis (FEA) software to predict the onset of permanent deformation in structural steel and aluminum components.

2. Principal Strain Analysis

While some theories, such as the Maximum Shear Strain Theory, are often discussed in terms of the differences between principal strains, these principal strains are themselves the roots of the characteristic equation:
$$ \det(\boldsymbol{\varepsilon} - \lambda \mathbf{I}) = 0 $$
The coefficients of this characteristic equation are precisely the strain invariants $I_1, I_2,$ and $I_3$. Thus, even when working with principal values, the underlying mathematical structure is governed by the invariants.

Computational Example

To illustrate the practical calculation of these values, consider a point in a material where the strain tensor in a Cartesian coordinate system is given by:
$$ \boldsymbol{\varepsilon} = \begin{bmatrix} 0.01 & 0.005 & 0 \ 0.005 & 0.02 & 0 \ 0 & 0 & 0.015 \end{bmatrix} $$

We can compute the invariants as follows:

• First Invariant ($I_1$):
  $$ I_1 = 0.01 + 0.02 + 0.015 = \mathbf{0.045} $$

• Third Invariant ($I_3$):
  $$ I_3 = \det(\boldsymbol{\varepsilon}) = 0.01(0.02 \times 0.015 - 0) - 0.005(0.005 \times 0.015 - 0) + 0 $$
  $$ I_3 = 0.000003 - 0.000000375 = \mathbf{0.000002625} $$

• Second Invariant ($I_2$):
  Using the sum of principal minors:
  $$ I_2 = \begin{vmatrix} 0.01 & 0.005 \ 0.005 & 0.02 \end{vmatrix} + \begin{vmatrix} 0.01 & 0 \ 0 & 0.015 \end{vmatrix} + \begin{vmatrix} 0.02 & 0 \ 0 & 0.015 \end{vmatrix} $$
  $$ I_2 = (0.0002 - 0.000025) + (0.00015) + (0.0003) = \mathbf{0.000625} $$

These invariant values provide a unique "fingerprint" of the deformation at that point. Regardless of how the coordinate axes are subsequently rotated, these three numbers will remain identical, ensuring consistent results in any numerical simulation.

Conclusion

Strain invariants serve as the essential bridge between abstract tensor algebra and practical engineering application. By providing a coordinate-independent description of deformation, they allow for the creation of robust, objective constitutive models. From the complex hyperelasticity of polymers to the critical yield predictions in metal forming, the mastery of strain invariants is fundamental to understanding and predicting how materials behave under real-world loading conditions.