Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Chapter 2 Lecture: Stress and Strain

4/7 Quiz QR Code

Purpose

To describe how the Earth deforms and how forces are transmitted through it, we need a framework for:

These two quantities are linked through constitutive relationships, which ultimately control:

In this chapter, we build that framework.


Learning Objectives

By the end of this lecture, you should be able to:


1. The Stress Tensor

1.1 Traction

Imagine cutting an infinitesimal plane inside a material. Forces act across that plane.

The force per unit area is called the traction:

t(n^)=(tx,ty,tz)\mathbf{t}(\hat{\mathbf{n}}) = (t_x, t_y, t_z)

where n^\hat{\mathbf{n}} is the unit normal to the plane.

We can always decompose traction into:

Key symmetry:

t(n^)=t(n^)\mathbf{t}(-\hat{\mathbf{n}}) = -\mathbf{t}(\hat{\mathbf{n}})

👉 This just says: flip the plane, and the force reverses.

Special case (fluids):

Fluids cannot support shear stress, so:

t=Pn^\mathbf{t} = -P \hat{\mathbf{n}}

Only normal forces (pressure) exist.


1.2 Stress Tensor

Instead of describing traction separately for every possible plane, we define a single object:

The stress tensor:

τ=[τxxτxyτxzτyxτyyτyzτzxτzyτzz]\boldsymbol{\tau} = \begin{bmatrix} \tau_{xx} & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & \tau_{yy} & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & \tau_{zz} \end{bmatrix}

Interpretation:

The key idea is:

t(n^)=τn^\mathbf{t}(\hat{\mathbf{n}}) = \boldsymbol{\tau} \hat{\mathbf{n}}

👉 The stress tensor maps plane orientation → traction.


1.3 Why is Stress Symmetric?

Physics requires that a tiny volume cannot spontaneously rotate.

This leads to:

τij=τji\tau_{ij} = \tau_{ji}

So:

👉 This is a big simplification.


1.4 Principal Stresses

There exist special directions where traction has no shear component:

τn^=λn^\boldsymbol{\tau} \hat{\mathbf{n}} = \lambda \hat{\mathbf{n}}

This is an eigenvalue problem:

(τλI)n^=0(\boldsymbol{\tau} - \lambda I)\hat{\mathbf{n}} = 0

👉 In these directions, stress is purely normal.


1.5 Principal Coordinate System

If we rotate into the principal directions:

τR=[τ1000τ2000τ3]\boldsymbol{\tau}_R = \begin{bmatrix} \tau_1 & 0 & 0 \\ 0 & \tau_2 & 0 \\ 0 & 0 & \tau_3 \end{bmatrix}

👉 No shear stresses remain.

This is often the most physically meaningful coordinate system.


1.6 Maximum Shear Stress

The largest shear stress occurs on planes at 45° to the principal axes:

τmax=τ1τ32\tau_{\max} = \frac{\tau_1 - \tau_3}{2}

👉 This is directly related to failure and faulting.


1.7 Hydrostatic vs Deviatoric Stress

We can split stress into:

Hydrostatic (mean) stress:

τm=τ1+τ2+τ33\tau_m = \frac{\tau_1 + \tau_2 + \tau_3}{3}

Deviatoric stress:

τD=ττmI\boldsymbol{\tau}_D = \boldsymbol{\tau} - \tau_m I

Key takeaway:

Earthquakes are driven by deviatoric stress, not hydrostatic pressure.


2. The Strain Tensor

2.1 Displacement Field

We describe motion using displacement:

u(r0,t)=rr0\mathbf{u}(\mathbf{r}_0, t) = \mathbf{r} - \mathbf{r}_0

2.2 From Displacement → Strain

Strain describes how displacement varies in space:

J=uJ = \nabla \mathbf{u}

This is the displacement gradient tensor.


2.3 Strain Tensor

We split deformation into:

J=e+ωJ = \mathbf{e} + \boldsymbol{\omega}

Strain is:

eij=12(uixj+ujxi)e_{ij} = \frac{1}{2} \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right)

👉 Only the symmetric part changes shape.


2.4 Physical Meaning

👉 Strain tells you how the material deforms, not where it moves.


2.5 Volume Change (Dilatation)

Δ=u=tr(e)\Delta = \nabla \cdot \mathbf{u} = \text{tr}(\mathbf{e})

2.6 Principal Strains

Strain behaves just like stress:

en^=λn^\mathbf{e} \hat{\mathbf{n}} = \lambda \hat{\mathbf{n}}

👉 There are directions with pure extension/compression and no shear.


2.7 Connection to Seismic Waves


2.8 Typical Strain Magnitudes

106\sim 10^{-6}

👉 Very small → justifies linear elasticity.


3. Stress–Strain Relationship (Elasticity)

To connect forces and deformation, we need a constitutive law.


3.1 General Linear Form

τij=cijklekl\tau_{ij} = c_{ijkl} e_{kl}

3.2 Isotropic Case

For most Earth materials (first approximation):

τij=λδijekk+2μeij\tau_{ij} = \lambda \delta_{ij} e_{kk} + 2\mu e_{ij}

Only two parameters:


3.3 Physical Interpretation

👉 μ\mu is especially important—it controls whether S-waves exist.


3.4 Seismic Velocities

Elasticity directly determines wave speeds:

α=λ+2μρ,β=μρ\alpha = \sqrt{\frac{\lambda + 2\mu}{\rho}}, \quad \beta = \sqrt{\frac{\mu}{\rho}}

Summary

This framework underpins everything we do in seismology: wave propagation, earthquake mechanics, and Earth structure.