The goal of any problem in solid mechanics is to determine the distribution of displacements, strains, and stresses throughout a body subjected to external forces. This requires a set of governing equations that are based on three core physical principles: the balance of forces (equilibrium), the geometry of deformation (kinematics), and the material’s response (constitutive law).

For a three-dimensional elastic body, we must solve for a total of 15 unknown field quantities at every point \((x, y, z)\) within the body.

The 15 Unknowns

The 15 unknowns can be grouped into three categories:

  1. Displacement Vector (3 unknowns): These describe how each point in the body moves.
    • \(u(x,y,z)\), \(v(x,y,z)\), \(w(x,y,z)\)
  2. Strain Tensor (6 independent unknowns): These describe the deformation (stretching and shearing) of the material. The strain tensor is symmetric (\(\epsilon_{ij} = \epsilon_{ji}\)), so it has 6 unique components.
    • Normal Strains: \(\epsilon_{xx}, \epsilon_{yy}, \epsilon_{zz}\)
    • Shear Strains: \(\epsilon_{xy}, \epsilon_{yz}, \epsilon_{xz}\)
  3. Stress Tensor (6 independent unknowns): These describe the internal forces acting on infinitesimal surfaces within the material. Due to the balance of moments, the stress tensor is also symmetric (\(\sigma_{ij} = \sigma_{ji}\)), giving it 6 unique components.
    • Normal Stresses: \(\sigma_{xx}, \sigma_{yy}, \sigma_{zz}\)
    • Shear Stresses: \(\sigma_{xy}, \sigma_{yz}, \sigma_{xz}\)

Total Unknowns = 3 (Displacements) + 6 (Strains) + 6 (Stresses) = 15 Quantities.

To solve for these 15 unknowns, we need an equal number of independent equations, which are provided by the fundamental laws of continuum mechanics.

The 15 Governing Equations

The 15 equations are derived from three fundamental principles:

1. Equilibrium Equations (3 Equations)

Considering infinitesimal element within the body, it follows from Newton’s second law (\(\Sigma \mathbf{F} = m\mathbf{a}\)) that the stress components satisfy the following equations of motion: \[ \frac{\partial \sigma_{xx}}{\partial x} + \frac{\partial \sigma_{yx}}{\partial y} + \frac{\partial \sigma_{zx}}{\partial z}+\rho b_x = \rho \frac{\partial ^2 u}{\partial t^2} \] \[ \frac{\partial \sigma_{xy}}{\partial x} + \frac{\partial \sigma_{yy}}{\partial y} + \frac{\partial \sigma_{zy}}{\partial z}+\rho b_y = \rho \frac{\partial ^2 v}{\partial t^2}\] \[ \frac{\partial \sigma_{xz}}{\partial x} + \frac{\partial \sigma_{yz}}{\partial y} + \frac{\partial \sigma_{zz}}{\partial z}+\rho b_z= \rho \frac{\partial ^2 w}{\partial t^2} \] In the above equations, \(\mathbf{u}=\begin{bmatrix} u & v & w\end{bmatrix}\) is the displacement field and \(\mathbf{b}\) is the body force per unit mass.

The equilibrium equations are often written in compact index notation as \[ \sigma_{ji,j}+\rho b_i=\rho \ddot{u}_i,\qquad (i=1, 2, 3) \] where \(_{,j}\) means differentiation with respect to \(x_j\), summation over the repeated index \(j\) is implied (Einstein notation), and a double dot is used to denote second derivative with respect to time.

Another form of writing the above equation is as follows \[ \nabla\cdot\boldsymbol{\sigma}+\rho\mathbf{b}=\rho \ddot{\mathbf{u}}. \]

2. Kinematic (Strain-Displacement) Equations (6 Equations) These are geometric relationships that define the components of the strain tensor in terms of the derivatives of the displacement vector. They are valid for the assumption of small deformations. \[ \epsilon_{xx} = \frac{\partial u}{\partial x} \] \[ \epsilon_{yy} = \frac{\partial v}{\partial y} \] \[ \epsilon_{zz} = \frac{\partial w}{\partial z} \] \[ \epsilon_{xy} = \frac{1}{2} \left( \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} \right) \] \[ \epsilon_{yz} = \frac{1}{2} \left( \frac{\partial v}{\partial z} + \frac{\partial w}{\partial y} \right) \] \[ \epsilon_{xz} = \frac{1}{2} \left( \frac{\partial u}{\partial z} + \frac{\partial w}{\partial x} \right) \] In a compact form: \[ \epsilon_{ij}=\frac{1}{2}\left(u_{i,j}+u_{j,i}\right), \] or \[ \boldsymbol{\epsilon}=\frac{1}{2}\left[\nabla \boldsymbol{\epsilon}+\nabla (\boldsymbol{\epsilon})^T\right]. \] 3. Constitutive Equations / Generalized Hooke’s Law (6 Equations) These equations describe the material’s intrinsic behavior by relating stress to strain. In general: \[ \sigma_{ij}=C_{ijkl}\epsilon_{kl} \] * Notice that summation over the repeated indices \(k\) and \(l\) is implied. That is, \[C_{ijkl}\epsilon_{kl}=\sum_{k=1}^3\sum_{l=1}^3 C_{ijkl}\epsilon_{kl}.\] For a linear, isotropic elastic material, this relationship is defined by two material constants, typically Young’s Modulus (\(E\)) and Poisson’s Ratio (\(\nu\)). \[ \sigma_{xx} = \frac{E}{(1+\nu)(1-2\nu)} [(1-\nu)\epsilon_{xx} + \nu(\epsilon_{yy} + \epsilon_{zz})] \] \[ \sigma_{yy} = \frac{E}{(1+\nu)(1-2\nu)} [(1-\nu)\epsilon_{yy} + \nu(\epsilon_{xx} + \epsilon_{zz})] \] \[ \sigma_{zz} = \frac{E}{(1+\nu)(1-2\nu)} [(1-\nu)\epsilon_{zz} + \nu(\epsilon_{xx} + \epsilon_{yy})] \] \[ \sigma_{xy} = \frac{E}{1+\nu} \epsilon_{xy} \] \[ \sigma_{yz} = \frac{E}{1+\nu} \epsilon_{yz} \] \[ \sigma_{xz} = \frac{E}{1+\nu} \epsilon_{xz} \] or in a compact way \[ \sigma_{ij}=\lambda\epsilon_{kk}\delta_{ij}+2\mu \epsilon_{ij}, \] or \[ \boldsymbol{\sigma}=\lambda(\operatorname{tr}\boldsymbol{\epsilon})\boldsymbol{I}+2\mu \boldsymbol{\epsilon}. \] #### Summary

The theory of linearized elasticity is a closed and complete mathematical framework. We have established a system of 15 unknowns and a corresponding set of 15 independent equations:

This equality ensures that the problem is mathematically well-posed. By applying appropriate boundary conditions (i.e., specifying the forces or displacements on the surface of the body), a unique solution for the stress, strain, and displacement fields can be determined.