Introduction to Finite Strain Theory Under Large Deformations

In classical linear elasticity, we often rely on the assumption of infinitesimal deformations. This simplification assumes that the displacement gradients are so small that the difference between the original and deformed shapes is negligible. However, when dealing with soft biological tissues, highly elastic polymers, or structural components undergoing extreme loading, this assumption collapses.

In such regimes, we encounter geometric nonlinearity (large changes in shape and orientation) and material nonlinearity (complex stress-strain relationships). To accurately capture these phenomena, we must transition from small-strain theory to Finite Strain Theory. This framework provides a mathematically rigorous way to describe deformation by accounting for the distinction between the initial state of a body and its current, deformed state.

1. Kinematics: The Geometry of Deformation

To describe large deformations, we must first define the spatial context in which the deformation occurs. We distinguish between two primary states:

• Reference Configuration ($\mathcal{B}_0$): The undeformed state of the body, where every material point is identified by its initial position vector $\mathbf{X}$.
• Current Configuration ($\mathcal{B}$): The deformed state of the body, where points are identified by their spatial coordinates $\mathbf{x}$.

The transition from the reference to the current state is governed by a continuous mapping function, known as the motion, denoted as $\boldsymbol{\chi}(\mathbf{X})$. Consequently, the position of a point in the current configuration is $\mathbf{x} = \boldsymbol{\chi}(\mathbf{X})$.

The Deformation Gradient ($\mathbf{F}$)

The fundamental quantity in finite strain kinematics is the deformation gradient, $\mathbf{F}$. It is defined as the gradient of the motion with respect to the reference coordinates:
$$\mathbf{F} = \frac{\partial \mathbf{x}}{\partial \mathbf{X}} = \nabla_{\mathbf{X}}\boldsymbol{\chi}$$
$\mathbf{F}$ is a second-order tensor that encapsulates all information regarding the local deformation, including stretching, shearing, and rigid-body rotation.

A critical scalar derived from $\mathbf{F}$ is the Jacobian, $J = \det \mathbf{F}$. The Jacobian represents the ratio of the deformed volume to the initial volume. In the case of incompressible materials (such as many rubbers and biological cells), the volume remains constant during deformation, leading to the constraint:
$$J = 1$$

2. Strain Tensors: Quantifying Deformation

Because $\mathbf{F}$ contains both rotation and stretching, it is not a pure measure of "strain." To isolate the actual shape change from rigid-body rotations, we utilize various strain tensors. These are generally categorized by whether they are defined relative to the reference or the current configuration.

Key Insight: While many measures exist, the Green–Lagrange strain is particularly favored in computational mechanics because it vanishes under rigid-body motion and simplifies the derivation of energy-based constitutive laws.

3. Stress Measures: The Force-Area Relationship

In small-strain theory, we use a single stress tensor. In finite strain theory, we must distinguish between forces acting on the original (undeformed) area and forces acting on the actual (deformed) area.

• Cauchy Stress ($\boldsymbol{\sigma}$): Often called the "true stress," this tensor represents the force per unit deformed area. It is the most physically intuitive measure for observing the actual state of internal forces within a material.
• First Piola–Kirchhoff Stress ($\mathbf{P}$): This is a "two-point" tensor that relates the force in the current configuration to the area in the reference configuration. It is frequently used in the weak form of Finite Element Method (FEM) equations.
• Second Piola–Kirchhoff Stress ($\mathbf{S}$): This tensor is purely defined in the reference configuration. It is symmetric and is mathematically "paired" with the Green–Lagrange strain, making it the standard choice for constructing hyperelastic potential functions.

The relationship between these measures is given by:
$$\mathbf{P} = J \boldsymbol{\sigma} \mathbf{F}^{-\mathrm{T}}, \quad \mathbf{S} = \mathbf{F}^{-1} \mathbf{P}$$

4. Constitutive Modeling: Hyperelasticity and Plasticity

A constitutive model defines how a specific material responds to a given strain. In the finite strain regime, we often use Hyperelasticity, where the stress is derived from a Strain Energy Density Function $W$.

4.1 Hyperelastic Models

For a hyperelastic material, the Second Piola–Kirchhoff stress is derived directly from the energy potential:
$$\mathbf{S} = 2 \frac{\partial W}{\partial \mathbf{C}}$$

Common models include:

• Neo-Hookean Model: A simple extension of Hooke's law for large deformations, suitable for many elastomers.
  $$W = \frac{\mu}{2}(\text{tr},\mathbf{C} - 3) - \mu \ln J + \frac{\lambda}{2}(\ln J)^2$$
• Mooney–Rivlin Model: A more sophisticated model that accounts for additional non-linearities by using invariants of the strain tensor, providing a better fit for rubber-like materials.

4.2 Large Deformation Plasticity

When materials undergo permanent deformation, we often employ the $J_2$ flow theory. Here, the plastic flow is typically modeled using the Hencky (logarithmic) strain, as it provides a more additive and physically consistent description of cumulative plastic deformation.

5. Numerical Illustration: Uniaxial Tension

To bridge theory and practice, let us consider a simple case: a Neo-Hookean elastic rod undergoing uniaxial stretching.

Given Parameters:

• Initial length $L_0 = 1$ m.
• Final length $L = 1.2$ m (a 20% stretch).
• Material constants: $\mu = 0.8$ MPa, $\lambda = 1.2$ MPa.

Step 1: Kinematics
The stretch ratio $\lambda_z = L/L_0 = 1.2$. The deformation gradient $\mathbf{F}$ is:
$$\mathbf{F} = \text{diag}(1.2, 1, 1)$$
The Jacobian $J = \det \mathbf{F} = 1.2$.

Step 2: Calculating Stress
Using the Neo-Hookean relation for the First Piola–Kirchhoff stress component $P_{zz}$:
$$P_{zz} = \mu(\lambda_z - \lambda_z^{-1}) + \lambda \ln J (\lambda_z^{-1})$$
Substituting the values:
$$P_{zz} = 0.8(1.2 - 0.8333) + 1.2 \ln(1.2) \times 0.8333 \approx 0.475 \text{ MPa}$$

Step 3: Finding True Stress
The Cauchy stress $\sigma_{zz}$ is:
$$\sigma_{zz} = \frac{P_{zz}}{J} = \frac{0.475}{1.2} \approx 0.396 \text{ MPa}$$

Observation: If we had used linear elasticity ($\sigma = E\varepsilon$), the result would likely be significantly higher. The finite strain approach correctly captures the geometric softening effect inherent in large-scale stretching.

6. Implementation Best Practices

When implementing these theories in numerical solvers (like custom FEM codes), several pitfalls must be avoided:

• Maintain Tensor Symmetry: Ensure that $\mathbf{C}$ and $\mathbf{E}$ remain symmetric during numerical integration to prevent non-physical energy generation.
• Numerical Stability: As $J \to 0$ (extreme compression), the terms involving $\mathbf{F}^{-1}$ or $\ln J$ can explode. Use regularization techniques (e.g., adding a small $\epsilon$ to the determinant) to maintain stability.
• The Consistent Tangent: For Newton-Raphson iterations to converge efficiently, you must provide the consistent tangent modulus (the derivative of the stress with respect to the deformation gradient). Using a simplified linear tangent will lead to poor convergence in highly non-linear problems.

7. Summary

Finite Strain Theory is an essential framework for modern mechanics. By utilizing the deformation gradient $\mathbf{F}$ as the cornerstone, it allows us to navigate the complex relationship between reference and current configurations. Through the careful selection of strain measures (like Green–Lagrange) and stress measures (like Cauchy or Piola–Kirchhoff), and by employing robust hyperelastic constitutive models, we can accurately simulate everything from the stretching of a rubber band to the complex deformation of human arteries.