Skip to main content

Plane stress/strain; 3D equilibrium equation

Short answer

Plane stress: thin plate with in-plane load, sigma_z=tau_xz=tau_yz=0. Plane strain: long uniform prism with constrained ends, eps_z=gamma_xz=gamma_yz=0. Three-dimensional equilibrium: dsigma_x/dx+dtau_xy/dy+dtau_xz/dz+X=0; dtau_xy/dx+dsigma_y/dy+dtau_yz/dz+Y=0; dtau_xz/dx+dtau_yz/dy+dsigma_z/dz+Z=0, or sigma_ij,j+bi=0.

Define plane stress and plane strain problems. Derive the differential equation of equilibrium for three-dimensional problems.

Plane stress and plane strain

Plane stress: a thin plate loaded only by forces in its own plane, so the stresses through the thickness are zero: σz=τxz=τyz=0\sigma_z=\tau_{xz}=\tau_{yz}=0. The stress state is two-dimensional, {σ}={σx,σy,τxy}T\{\sigma\}=\{\sigma_x,\sigma_y,\tau_{xy}\}^T, although the strain εz=−νE(σx+σy)≠0\varepsilon_z=-\dfrac{\nu}{E}(\sigma_x+\sigma_y)\neq0.

Plane strain: a long prism (length much larger than cross-section) with uniform cross-section and load that does not vary along its length and is perpendicular to it. Movement along the axis is prevented, so εz=γxz=γyz=0\varepsilon_z=\gamma_{xz}=\gamma_{yz}=0, but σz=ν(σx+σy)≠0\sigma_z=\nu(\sigma_x+\sigma_y)\neq0.

Examples: thin plate with a hole under in-plane load (plane stress); long gravity dam or retaining wall of uniform section (plane strain).

Differential equations of equilibrium in 3D

Consider an infinitesimal element dx dy dzdx\,dy\,dz in equilibrium with body forces per unit volume X,Y,ZX, Y, Z (e.g. self-weight). Stresses vary from face to face, e.g. the normal stress on the face at x+dxx+dx is σx+∂σx∂xdx\sigma_x+\dfrac{\partial\sigma_x}{\partial x}dx.

        y
        ^     sy + dsy
        |   +----------+
  txy   |   |          | sx + dsx
  <-----|   |  dx dy   |---->
   sx   +---+----------+---> x
            sy

Equilibrium of forces in the xx direction

Forces = stress ×\times area. Faces normal to xx have area dy dzdy\,dz, those normal to yy have area dx dzdx\,dz and those normal to zz have area dx dydx\,dy:

(∂σx∂xdx)dy dz+(∂τyx∂ydy)dx dz+(∂τzx∂zdz)dx dy+X dx dy dz=0\Big(\frac{\partial\sigma_x}{\partial x}dx\Big)dy\,dz+\Big(\frac{\partial\tau_{yx}}{\partial y}dy\Big)dx\,dz+\Big(\frac{\partial\tau_{zx}}{\partial z}dz\Big)dx\,dy+X\,dx\,dy\,dz=0

(the terms σx\sigma_x, τyx\tau_{yx}, τzx\tau_{zx} on opposite faces cancel). Dividing by dx dy dzdx\,dy\,dz:

∂σx∂x+∂τxy∂y+∂τxz∂z+X=0\frac{\partial\sigma_x}{\partial x}+\frac{\partial\tau_{xy}}{\partial y}+\frac{\partial\tau_{xz}}{\partial z}+X=0

using τyx=τxy\tau_{yx}=\tau_{xy} and τzx=τxz\tau_{zx}=\tau_{xz} (from moment equilibrium).

Similarly for yy and zz

∂τxy∂x+∂σy∂y+∂τyz∂z+Y=0∂τxz∂x+∂τyz∂y+∂σz∂z+Z=0\begin{aligned} \frac{\partial\tau_{xy}}{\partial x}+\frac{\partial\sigma_y}{\partial y}+\frac{\partial\tau_{yz}}{\partial z}+Y&=0\\ \frac{\partial\tau_{xz}}{\partial x}+\frac{\partial\tau_{yz}}{\partial y}+\frac{\partial\sigma_z}{\partial z}+Z&=0 \end{aligned}

Moment equilibrium about the three axes gives τxy=τyx\tau_{xy}=\tau_{yx}, τyz=τzy\tau_{yz}=\tau_{zy}, τzx=τxz\tau_{zx}=\tau_{xz} (symmetry of the stress tensor).

In compact form:  σij,j+bi=0\ \sigma_{ij,j}+b_i=0.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Computational Techniques in Civil Engineering chapter 3 questions →