The Deformation Gradient Tensor

The deformation gradient $\mathbf{F}$ is the most important quantity in large deformation mechanics. It captures how infinitesimal material elements stretch, rotate, and shear at each point.

Definition

The deformation gradient is the gradient of the deformation mapping:

$$\mathbf{F} = \frac{\partial \boldsymbol{\varphi}}{\partial \mathbf{X}} = \nabla_{\mathbf{X}} \boldsymbol{\varphi}$$

Physical Interpretation

Consider an infinitesimal line element $d\mathbf{X}$ in the reference configuration. After deformation, it becomes:

$$d\mathbf{x} = \mathbf{F} \cdot d\mathbf{X}$$

So $\mathbf{F}$ tells us how small material fibers transform: they can stretch, rotate, and change direction.

Example: Circle to Ellipse

An infinitesimal circle in the reference configuration becomes an ellipse in the current configuration. The principal stretches (eigenvalues of $\sqrt{\mathbf{F}^T\mathbf{F}}$) determine the ellipse's axes.

Circles transform into ellipses. The shape of each ellipse shows the local F tensor.

Computing F from Deformation Mapping

Given the deformation mapping from the previous lesson:

$$\varphi_1 = (1 - 0.25t) X_1 + t(X_2-2)^2 + 2.5t$$

$$\varphi_2 = (1 + 0.25t) X_2 + 0.25t$$

The deformation gradient is:

$$\mathbf{F} = \begin{bmatrix} \frac{\partial \varphi_1}{\partial X_1} & \frac{\partial \varphi_1}{\partial X_2} \\ \frac{\partial \varphi_2}{\partial X_1} & \frac{\partial \varphi_2}{\partial X_2} \end{bmatrix} = \begin{bmatrix} 1 - 0.25t & 2t(X_2-2) \\ 0 & 1 + 0.25t \end{bmatrix}$$

Notice that $\mathbf{F}$ varies with position — it's a tensor field, not a single value.

Python Code

import numpy as np

def deformation_mapping(X, t):
    """x = phi(X, t)"""
    x1 = (1 - 0.25*t) * X[0] + t * (X[1] - 2)**2 + 2.5*t
    x2 = (1 + 0.25*t) * X[1] + 0.25*t
    return np.array([x1, x2])

def deformation_gradient(X, t):
    """F = d(phi)/dX computed analytically"""
    F = np.array([
        [1 - 0.25*t, 2*t*(X[1] - 2)],
        [0, 1 + 0.25*t]
    ])
    return F

# Example: compute F at a point
X = np.array([1.0, 3.0])
t = 0.5
F = deformation_gradient(X, t)
print(f"F = \n{F}")
print(f"det(F) = {np.linalg.det(F):.3f}")

Key Properties of F

1. Determinant (Volume Change)

$$J = \det(\mathbf{F})$$

$J > 1$: Volume increases (expansion)
$J < 1$: Volume decreases (compression)
$J = 1$: Incompressible (rubber, fluids)
$J \leq 0$: Invalid deformation (element inversion)

2. Polar Decomposition

$$\mathbf{F} = \mathbf{R} \cdot \mathbf{U} = \mathbf{V} \cdot \mathbf{R}$$

$\mathbf{R}$: Rotation tensor (orthogonal)
$\mathbf{U}$: Right stretch tensor (symmetric positive definite)
$\mathbf{V}$: Left stretch tensor (symmetric positive definite)

This separates rotation from stretching — useful for material modeling.

Automotive Application: Tire Rubber

During cornering, tire sidewall rubber experiences:

Large stretch in circumferential direction
Compression in radial direction
Shear due to lateral forces

Computing $\mathbf{F}$ at each point reveals:

Where material is most stretched (fatigue hotspots)
Where det(F) < 1 (compression, potential buckling)
Principal stretch directions (cord ply loading)

Key Takeaways

F transforms infinitesimal elements: $d\mathbf{x} = \mathbf{F} \cdot d\mathbf{X}$
F is computed as $\partial\varphi_i / \partial X_J$
det(F) = volume ratio (must be positive)
Polar decomposition separates rotation from stretch
F varies with position — it's a tensor field

What's Next

In the next lesson, we'll explore different stress measures and how they relate to the deformation gradient. Understanding which stress measure to use is essential for interpreting FEA results.