Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

Continuum mechanics elasticity review.

Posted by peeterjoot on April 23, 2012

[Click here for a PDF of this post with nicer formatting]


Review of key ideas and equations from the theory of elasticity portion of the class.

Strain Tensor

Identifying a point in a solid with coordinates x_i and the coordinates of that portion of the solid after displacement, we formed the difference as a measure of the displacement

\begin{aligned}u_i = x_i' - x_i.\end{aligned} \hspace{\stretch{1}}(2.1)

With du_i = {\partial {u_i}}/{\partial {x_j}} dx_j, we computed the difference in length (squared) for an element of the displaced solid and found

\begin{aligned}dx_k' dx_k' - dx_k dx_k = \left( \frac{\partial {u_j}}{\partial {x_i}} + \frac{\partial {u_i}}{\partial {x_j}} + \frac{\partial {u_k}}{\partial {x_i}} \frac{\partial {u_k}}{\partial {x_j}} \right) dx_i dx_j,\end{aligned} \hspace{\stretch{1}}(2.2)

or defining the \textit{strain tensor} e_{ij}, we have

\begin{aligned}(d\mathbf{x}')^2 - (d\mathbf{x})^2= 2 e_{ij} dx_i dx_j\end{aligned} \hspace{\stretch{1}}(2.3a)

\begin{aligned}e_{ij}=\frac{1}{{2}}\left( \frac{\partial {u_j}}{\partial {x_i}} + \frac{\partial {u_i}}{\partial {x_j}} + \frac{\partial {u_k}}{\partial {x_i}} \frac{\partial {u_k}}{\partial {x_j}} \right).\end{aligned} \hspace{\stretch{1}}(2.3b)

In this course we use only the linear terms and write

\begin{aligned}e_{ij}=\frac{1}{{2}}\left( \frac{\partial {u_j}}{\partial {x_i}} + \frac{\partial {u_i}}{\partial {x_j}} \right).\end{aligned} \hspace{\stretch{1}}(2.4)

Unresolved: Relating displacement and position by strain

In [1] it is pointed out that this strain tensor simply relates the displacement vector coordinates u_i to the coordinates at the point at which it is measured

\begin{aligned}u_i = e_{ij} x_j.\end{aligned} \hspace{\stretch{1}}(2.5)

When we get to fluid dynamics we perform a linear expansion of du_i and find something similar

\begin{aligned}dx_i' - dx_i = du_i = \frac{\partial {u_i}}{\partial {x_k}} dx_k = e_{ij} dx_k + \omega_{ij} dx_k\end{aligned} \hspace{\stretch{1}}(2.6)


\begin{aligned}\omega_{ij} = \frac{1}{{2}} \left( \frac{\partial {u_j}}{\partial {x_i}} +\frac{\partial {u_i}}{\partial {x_j}} \right).\end{aligned} \hspace{\stretch{1}}(2.7)

Except for the antisymmetric term, note the structural similarity of 2.5 and 2.6. Why is it that we neglect the vorticity tensor in statics?

Diagonal strain representation.

In a basis for which the strain tensor is diagonal, it was pointed out that we can write our difference in squared displacement as (for k = 1, 2, 3, no summation convention)

\begin{aligned}(dx_k')^2 - (dx_k)^2 = 2 e_{kk} dx_k dx_k\end{aligned} \hspace{\stretch{1}}(2.8)

from which we can rearrange, take roots, and apply a first order Taylor expansion to find (again no summation convention)

\begin{aligned}dx_k' \approx (1 + e_{kk}) dx_k.\end{aligned} \hspace{\stretch{1}}(2.9)

An approximation of the displaced volume was then found in terms of the strain tensor trace (summation convention back again)

\begin{aligned}dV' \approx (1 + e_{kk}) dV,\end{aligned} \hspace{\stretch{1}}(2.10)

allowing us to identify this trace as a relative difference in displaced volume

\begin{aligned}e_{kk} \approx \frac{dV' - dV}{dV}.\end{aligned} \hspace{\stretch{1}}(2.11)

Strain in cylindrical coordinates.

Useful in many practice problems are the cylindrical coordinate representation of the strain tensor

\begin{aligned}2 e_{rr} &= \frac{\partial {u_r}}{\partial {r}}  \\ 2 e_{\phi\phi} &= \frac{1}{{r}} \frac{\partial {u_\phi}}{\partial {\phi}} +\frac{1}{{r}} u_r  \\ 2 e_{zz} &= \frac{\partial {u_z}}{\partial {z}}  \\ 2 e_{zr} &= \frac{\partial {u_r}}{\partial {z}} + \frac{\partial {u_z}}{\partial {r}} \\ 2 e_{r\phi} &= \frac{\partial {u_\phi}}{\partial {r}} - \frac{1}{{r}} u_\phi + \frac{1}{{r}} \frac{\partial {u_r}}{\partial {\phi}} \\ 2 e_{\phi z} &= \frac{\partial {u_\phi}}{\partial {z}} +\frac{1}{{r}} \frac{\partial {u_z}}{\partial {\phi}}.\end{aligned} \hspace{\stretch{1}}(2.12)

This can be found in [2]. It was not derived there or in class, but is not too hard, even using the second order methods we used for the Cartesian form of the tensor.

An easier way to do this derivation (and understand what the coordinates represent) follows from the relation found in section 6 of [3]

\begin{aligned}2 \mathbf{e}_i e_{ij} n_j = 2 (\hat{\mathbf{n}} \cdot \boldsymbol{\nabla}) \mathbf{u} + \hat{\mathbf{n}} \times (\boldsymbol{\nabla} \times \mathbf{u}),\end{aligned} \hspace{\stretch{1}}(2.18)

where \hat{\mathbf{n}} is the normal to the surface at which we are measuring a force applied to the solid (our Cauchy tetrahedron).

The cylindrical tensor coordinates of 2.12 follow from
2.18 nicely taking \hat{\mathbf{n}} = \hat{\mathbf{r}}, \hat{\boldsymbol{\phi}}, \hat{\mathbf{z}} in turn.

Compatibility condition.

For a 2D strain tensor we found an interrelationship between the components of the strain tensor

\begin{aligned}2 \frac{\partial^2 e_{12}}{\partial x_1 \partial x_2} =\frac{\partial^2 {{e_{22}}}}{\partial {{x_1}}^2} +\frac{\partial^2 {{e_{11}}}}{\partial {{x_2}}^2},\end{aligned} \hspace{\stretch{1}}(2.19)

and called this the compatibility condition. It was claimed, but not demonstrated that this is what is required to ensure a deformation maintained a coherent solid geometry.

I wasn’t able to find any references to this compatibility condition in any of the texts I have, but found [4], [5], and [6]. It’s not terribly surprising to see Christoffel symbol and differential forms references on those pages, since one can imagine that we’d wish to look at the mappings of all the points in the object as it undergoes the transformation from the original to the deformed state.

Even with just three points in a plane, say \mathbf{a}, \mathbf{b}, \mathbf{c}, the general deformation of an object doesn’t seem like it’s the easiest thing to describe. We can imagine that these have trajectories in the deformation process \mathbf{a} = \mathbf{a}(\alpha, \mathbf{b} = \mathbf{b}(\beta), \mathbf{c} = \mathbf{c}(\gamma), with \mathbf{a}', \mathbf{b}', \mathbf{c}' at the end points of the trajectories. We’d want to look at displacement vectors \mathbf{u}_a, \mathbf{u}_b, \mathbf{u}_c along each of these trajectories, and then see how they must be related. Doing that carefully must result in this compatibility condition.

Stress tensor.

By sought and found a representation of the force per unit area acting on a body by expressing the components of that force as a set of divergence relations

\begin{aligned}f_i = \partial_k \sigma_{i k},\end{aligned} \hspace{\stretch{1}}(3.20)

and call the associated tensor \sigma_{ij} the \textit{stress}.

Unlike the strain, we don’t have any expectation that this tensor is symmetric, and identify the diagonal components (no sum) \sigma_{i i} as quantifying the amount of compressive or contractive force per unit area, whereas the cross terms of the stress tensor introduce shearing deformations in the solid.

With force balance arguments (the Cauchy tetrahedron) we found that the force per unit area on the solid, for a surface with unit normal pointing into the solid, was

\begin{aligned}\mathbf{t} = \mathbf{e}_i t_i = \mathbf{e}_i \sigma_{ij} n_j.\end{aligned} \hspace{\stretch{1}}(3.21)

Constitutive relation.

In the scope of this course we considered only Newtonian materials, those for which the stress and strain tensors are linearly related

\begin{aligned}\sigma_{ij} = c_{ijkl} e_{kl},\end{aligned} \hspace{\stretch{1}}(3.22)

and further restricted our attention to isotropic materials, which can be shown to have the form

\begin{aligned}\sigma_{ij} = \lambda e_{kk} \delta_{ij} + 2 \mu e_{ij},\end{aligned} \hspace{\stretch{1}}(3.23)

where \lambda and \mu are the Lame parameters and \mu is called the shear modulus (and viscosity in the context of fluids).

By computing the trace of the stress \sigma_{ii} we can invert this to find

\begin{aligned}2 \mu e_{ij} = \sigma_{ij} - \frac{\lambda}{3 \lambda + 2 \mu} \sigma_{kk} \delta_{ij}.\end{aligned} \hspace{\stretch{1}}(3.24)

Uniform hydrostatic compression.

With only normal components of the stress (no shear), and the stress having the same value in all directions, we find

\begin{aligned}\sigma_{ij} = ( 3 \lambda + 2 \mu ) e_{ij},\end{aligned} \hspace{\stretch{1}}(3.25)

and identify this combination -3 \lambda - 2 \mu as the pressure, linearly relating the stress and strain tensors

\begin{aligned}\sigma_{ij} = -p e_{ij}.\end{aligned} \hspace{\stretch{1}}(3.26)

With e_{ii} = (dV' - dV)/dV = \Delta V/V, we formed the Bulk modulus K with the value

\begin{aligned}K = \left( \lambda + \frac{2 \mu}{3} \right) = -\frac{p V}{\Delta V}.\end{aligned} \hspace{\stretch{1}}(3.27)

Uniaxial stress. Young’s modulus. Poisson’s ratio.

For the special case with only one non-zero stress component (we used \sigma_{11}) we were able to compute Young’s modulus E, the ratio between stress and strain in that direction

\begin{aligned}E = \frac{\sigma_{11}}{e_{11}} = \frac{\mu(3 \lambda + 2 \mu)}{\lambda + \mu }  = \frac{3 K \mu}{K + \mu/3}.\end{aligned} \hspace{\stretch{1}}(3.28)

Just because only one component of the stress is non-zero, does not mean that we have no deformation in any other directions. Introducing Poisson’s ratio \nu in terms of the ratio of the strains relative to the strain in the direction of the force we write and then subsequently found

\begin{aligned}\nu = -\frac{e_{22}}{e_{11}} = -\frac{e_{33}}{e_{11}} = \frac{\lambda}{2(\lambda + \mu)}.\end{aligned} \hspace{\stretch{1}}(3.29)

We were also able to find

We can also relate the Poisson’s ratio \nu to the shear modulus \mu

\begin{aligned}\mu = \frac{E}{2(1 + \nu)}\end{aligned} \hspace{\stretch{1}}(3.30)

\begin{aligned}\lambda = \frac{E \nu}{(1 - 2 \nu)(1 + \nu)}\end{aligned} \hspace{\stretch{1}}(3.31)

\begin{aligned}e_{11} &= \frac{1}{{E}}\left( \sigma_{11} - \nu(\sigma_{22} + \sigma_{33}) \right) \\ e_{22} &= \frac{1}{{E}}\left( \sigma_{22} - \nu(\sigma_{11} + \sigma_{33}) \right) \\ e_{33} &= \frac{1}{{E}}\left( \sigma_{33} - \nu(\sigma_{11} + \sigma_{22}) \right)\end{aligned} \hspace{\stretch{1}}(3.32)

Displacement propagation

It was argued that the equation relating the time evolution of a one of the vector displacement coordinates was given by

\begin{aligned}\rho \frac{\partial^2 {{u_i}}}{\partial {{t}}^2} = \frac{\partial {\sigma_{ij}}}{\partial {x_j}} + f_i,\end{aligned} \hspace{\stretch{1}}(4.35)

where the divergence term {\partial {\sigma_{ij}}}/{\partial {x_j}} is the internal force per unit volume on the object and f_i is the external force. Employing the constitutive relation we showed that this can be expanded as

\begin{aligned}\rho \frac{\partial^2 {{u_i}}}{\partial {{t}}^2} = (\lambda + \mu) \frac{\partial^2 u_k}{\partial x_i \partial x_k}+ \mu\frac{\partial^2 u_i}{\partial x_j^2},\end{aligned} \hspace{\stretch{1}}(4.36)

or in vector form

\begin{aligned}\rho \frac{\partial^2 {\mathbf{u}}}{\partial {{t}}^2} = (\lambda + \mu) \boldsymbol{\nabla} (\boldsymbol{\nabla} \cdot \mathbf{u}) + \mu \boldsymbol{\nabla}^2 \mathbf{u}.\end{aligned} \hspace{\stretch{1}}(4.37)


Operating on 4.37 with the divergence operator, and writing \Theta = \boldsymbol{\nabla} \cdot \mathbf{u}, a quantity that was our relative change in volume in the diagonal strain basis, we were able to find this divergence obeys a wave equation

\begin{aligned}\frac{\partial^2 {{\Theta}}}{\partial {{t}}^2} = \frac{\lambda + 2 \mu}{\rho} \boldsymbol{\nabla}^2 \Theta.\end{aligned} \hspace{\stretch{1}}(4.38)

We called these P-waves.


Similarly, operating on 4.37 with the curl operator, and writing \boldsymbol{\omega} = \boldsymbol{\nabla} \times \mathbf{u}, we were able to find this curl also obeys a wave equation

\begin{aligned}\rho \frac{\partial^2 {{\boldsymbol{\omega}}}}{\partial {{t}}^2} = \mu \boldsymbol{\nabla}^2 \boldsymbol{\omega}.\end{aligned} \hspace{\stretch{1}}(4.39)

These we called S-waves. We also noted that the (transverse) compression waves (P-waves) with speed C_T = \sqrt{\mu/\rho}, traveled faster than the (longitudinal) vorticity (S) waves with speed C_L = \sqrt{(\lambda + 2 \mu)/\rho} since \lambda > 0 and \mu > 0, and

\begin{aligned}\frac{C_L}{C_T} = \sqrt{\frac{ \lambda + 2 \mu}{\mu}} = \sqrt{ \frac{\lambda}{\mu} + 2}.\end{aligned} \hspace{\stretch{1}}(4.40)

Scalar and vector potential representation.

Assuming a vector displacement representation with gradient and curl components

\begin{aligned}\mathbf{u} = \boldsymbol{\nabla} \phi + \boldsymbol{\nabla} \times \mathbf{H},\end{aligned} \hspace{\stretch{1}}(4.41)

We found that the displacement time evolution equation split nicely into curl free and divergence free terms

\begin{aligned}\boldsymbol{\nabla}\left(\rho \frac{\partial^2 {{\phi}}}{\partial {{t}}^2} - (\lambda + 2\mu) \boldsymbol{\nabla}^2 \phi\right)+\boldsymbol{\nabla} \times\left(\rho \frac{\partial^2 {\mathbf{H}}}{\partial {{t}}^2} - \mu \boldsymbol{\nabla}^2 \mathbf{H}\right)= 0.\end{aligned} \hspace{\stretch{1}}(4.42)

When neglecting boundary value effects this could be written as a pair of independent equations

\begin{aligned}\rho \frac{\partial^2 {{\phi}}}{\partial {{t}}^2} - (\lambda + 2\mu) \boldsymbol{\nabla}^2 \phi = 0\end{aligned} \hspace{\stretch{1}}(4.43a)

\begin{aligned}\rho \frac{\partial^2 {\mathbf{H}}}{\partial {{t}}^2} - \mu \boldsymbol{\nabla}^2 \mathbf{H}= 0.\end{aligned} \hspace{\stretch{1}}(4.43b)

This are the irrotational (curl free) P-wave and solenoidal (divergence free) S-wave equations respectively.

Phasor description.

It was mentioned that we could assume a phasor representation for our potentials, writing

\begin{aligned}\phi = A \exp\left( i ( \mathbf{k} \cdot \mathbf{x} - \omega t) \right) \end{aligned} \hspace{\stretch{1}}(4.44a)

\begin{aligned}\mathbf{H} = \mathbf{B} \exp\left( i ( \mathbf{k} \cdot \mathbf{x} - \omega t) \right)\end{aligned} \hspace{\stretch{1}}(4.44b)


\begin{aligned}\mathbf{u} = i \mathbf{k} \phi + i \mathbf{k} \times \mathbf{H}.\end{aligned} \hspace{\stretch{1}}(4.45)

We did nothing with neither the potential nor the phasor theory for solid displacement time evolution, and presumably won’t on the exam either.

Some wave types

Some time was spent on non-qualitative descriptions and review of descriptions for solutions to the time evolution equations we did not attempt

  1. P-waves [7]. Irrotational, non volume preserving body wave.
  2. S-waves [8]. Divergence free body wave. Shearing forces are present and volume is preserved (slower than S-waves)
  3. Rayleigh wave [9]. A surface wave that propagates near the surface of a body without penetrating into it.
  4. Love wave [10]. A polarized shear surface wave with the shear displacements moving perpendicular to the direction of propagation.

For reasons that aren’t clear both the midterm and last years final ask us to spew this sort of stuff (instead of actually trying to do something analytic associated with them).


[1] R.P. Feynman, R.B. Leighton, and M.L. Sands. Feynman lectures on physics.[Lectures on physics], chapter Elastic Materials. Addison-Wesley Publishing Company. Reading, Massachusetts, 1963.

[2] L.D. Landau, EM Lifshitz, JB Sykes, WH Reid, and E.H. Dill. Theory of Elasticity: Vol. 7 of Course of Theoretical Physics. 1960.

[3] D.J. Acheson. Elementary fluid dynamics. Oxford University Press, USA, 1990.

[4] Wikipedia. Compatibility (mechanics) — wikipedia, the free encyclopedia [online]. 2011. [Online; accessed 23-April-2012].

[5] Wikipedia. Infinitesimal strain theory — wikipedia, the free encyclopedia [online]. 2012. [Online; accessed 23-April-2012].

[6] Wikipedia. Saint-venant’s compatibility condition — wikipedia, the free encyclopedia [online]. 2011. [Online; accessed 23-April-2012].\%27s_compatibility_condition&oldid=436103127.

[7] Wikipedia. P-wave — wikipedia, the free encyclopedia [online]. 2012. [Online; accessed 1-February-2012].

[8] Wikipedia. S-wave — wikipedia, the free encyclopedia [online]. 2011. [Online; accessed 1-February-2012].

[9] Wikipedia. Rayleigh wave — wikipedia, the free encyclopedia [online]. 2012. [Online; accessed 4-February-2012].

[10] Wikipedia. Love wave — wikipedia, the free encyclopedia [online]. 2012. [Online; accessed 4-February-2012].


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s

%d bloggers like this: