Physics:Michell solution

The Michell solution is a general solution to the elasticity equations in polar coordinates ([math]\displaystyle{ r, \theta \, }[/math]) developed by J. H. Michell. The solution is such that the stress components are in the form of a Fourier series in [math]\displaystyle{ \theta \, }[/math]. Michell showed that the general solution can be expressed in terms of an Airy stress function of the form

[math]\displaystyle{ \begin{align} \varphi(r,\theta) &= A_0~r^2 + B_0~r^2~\ln(r) + C_0~\ln(r) \\ & + \left(I_0~r^2 + I_1~r^2~\ln(r) + I_2~\ln(r) + I_3~\right) \theta \\ & + \left(A_1~r + B_1~r^{-1} + B_1'~r~\theta + C_1~r^3 + D_1~r~\ln(r)\right) \cos\theta \\ & + \left(E_1~r + F_1~r^{-1} + F_1'~r~\theta + G_1~r^3 + H_1~r~\ln(r)\right) \sin\theta \\ & + \sum_{n=2}^{\infty} \left(A_n~r^n + B_n~r^{-n} + C_n~r^{n+2} + D_n~r^{-n+2}\right)\cos(n\theta) \\ & + \sum_{n=2}^{\infty} \left(E_n~r^n + F_n~r^{-n} + G_n~r^{n+2} + H_n~r^{-n+2}\right)\sin(n\theta) \end{align} }[/math]

The terms [math]\displaystyle{ A_1~r~\cos\theta\, }[/math] and [math]\displaystyle{ E_1~r~\sin\theta\, }[/math] define a trivial null state of stress and are ignored.

Stress components

The stress components can be obtained by substituting the Michell solution into the equations for stress in terms of the Airy stress function (in cylindrical coordinates). A table of stress components is shown below.

[math]\displaystyle{ \varphi }[/math]	[math]\displaystyle{ \sigma_{rr}\, }[/math]	[math]\displaystyle{ \sigma_{r\theta}\, }[/math]	[math]\displaystyle{ \sigma_{\theta\theta}\, }[/math]
[math]\displaystyle{ r^2\, }[/math]	[math]\displaystyle{ 2 }[/math]	[math]\displaystyle{ 0 }[/math]	[math]\displaystyle{ 2 }[/math]
[math]\displaystyle{ r^2~\ln r }[/math]	[math]\displaystyle{ 2~\ln r + 1 }[/math]	[math]\displaystyle{ 0 }[/math]	[math]\displaystyle{ 2~\ln r + 3 }[/math]
[math]\displaystyle{ \ln r\, }[/math]	[math]\displaystyle{ r^{-2}\, }[/math]	[math]\displaystyle{ 0 }[/math]	[math]\displaystyle{ -r^{-2}\, }[/math]
[math]\displaystyle{ \theta\, }[/math]	[math]\displaystyle{ 0 }[/math]	[math]\displaystyle{ r^{-2}\, }[/math]	[math]\displaystyle{ 0 }[/math]
[math]\displaystyle{ r^3~\cos\theta \, }[/math]	[math]\displaystyle{ 2~r~\cos\theta \, }[/math]	[math]\displaystyle{ 2~r~\sin\theta \, }[/math]	[math]\displaystyle{ 6~r~\cos\theta \, }[/math]
[math]\displaystyle{ r\theta~\cos\theta \, }[/math]	[math]\displaystyle{ -2~r^{-1}~\sin\theta \, }[/math]	[math]\displaystyle{ 0 }[/math]	[math]\displaystyle{ 0 }[/math]
[math]\displaystyle{ r~\ln r~\cos\theta \, }[/math]	[math]\displaystyle{ r^{-1}~\cos\theta \, }[/math]	[math]\displaystyle{ r^{-1}~\sin\theta \, }[/math]	[math]\displaystyle{ r^{-1}~\cos\theta \, }[/math]
[math]\displaystyle{ r^{-1}~\cos\theta \, }[/math]	[math]\displaystyle{ -2~r^{-3}~\cos\theta \, }[/math]	[math]\displaystyle{ -2~r^{-3}~\sin\theta \, }[/math]	[math]\displaystyle{ 2~r^{-3}~\cos\theta \, }[/math]
[math]\displaystyle{ r^3~\sin\theta \, }[/math]	[math]\displaystyle{ 2~r~\sin\theta \, }[/math]	[math]\displaystyle{ -2~r~\cos\theta \, }[/math]	[math]\displaystyle{ 6~r~\sin\theta \, }[/math]
[math]\displaystyle{ r\theta~\sin\theta \, }[/math]	[math]\displaystyle{ 2~r^{-1}~\cos\theta \, }[/math]	[math]\displaystyle{ 0 }[/math]	[math]\displaystyle{ 0 }[/math]
[math]\displaystyle{ r~\ln r~\sin\theta \, }[/math]	[math]\displaystyle{ r^{-1}~\sin\theta \, }[/math]	[math]\displaystyle{ -r^{-1}~\cos\theta \, }[/math]	[math]\displaystyle{ r^{-1}~\sin\theta \, }[/math]
[math]\displaystyle{ r^{-1}~\sin\theta \, }[/math]	[math]\displaystyle{ -2~r^{-3}~\sin\theta \, }[/math]	[math]\displaystyle{ 2~r^{-3}~\cos\theta \, }[/math]	[math]\displaystyle{ 2~r^{-3}~\sin\theta \, }[/math]
[math]\displaystyle{ r^{n+2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -(n+1)(n-2)~r^n~\cos(n\theta) \, }[/math]	[math]\displaystyle{ n(n+1)~r^n~\sin(n\theta) \, }[/math]	[math]\displaystyle{ (n+1)(n+2)~r^n~\cos(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n+2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -(n+2)(n-1)~r^{-n}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -n(n-1)~r^{-n}~\sin(n\theta)\, }[/math]	[math]\displaystyle{ (n-1)(n-2)~r^{-n}~\cos(n\theta) }[/math]
[math]\displaystyle{ r^n~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -n(n-1)~r^{n-2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ n(n-1)~r^{n-2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ n(n-1)~r^{n-2}~\cos(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -n(n+1)~r^{-n-2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -n(n+1)~r^{-n-2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ n(n+1)~r^{-n-2}~\cos(n\theta) \, }[/math]
[math]\displaystyle{ r^{n+2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -(n+1)(n-2)~r^n~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -n(n+1)~r^n~\cos(n\theta) \, }[/math]	[math]\displaystyle{ (n+1)(n+2)~r^n~\sin(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n+2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -(n+2)(n-1)~r^{-n}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ n(n-1)~r^{-n}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ (n-1)(n-2)~r^{-n}~\sin(n\theta) \, }[/math]
[math]\displaystyle{ r^n~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -n(n-1)~r^{n-2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -n(n-1)~r^{n-2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ n(n-1)~r^{n-2}~\sin(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -n(n+1)~r^{-n-2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ n(n+1)~r^{-n-2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ n(n+1)~r^{-n-2}~\sin(n\theta) \, }[/math]

Displacement components

Displacements [math]\displaystyle{ (u_r, u_\theta) }[/math] can be obtained from the Michell solution by using the stress-strain and strain-displacement relations. A table of displacement components corresponding the terms in the Airy stress function for the Michell solution is given below. In this table

[math]\displaystyle{ \kappa = \begin{cases} 3 - 4~\nu & \rm{for~plane~strain} \\ \cfrac{3 - \nu}{1 + \nu} & \rm{for~plane~stress} \\ \end{cases} }[/math]

where [math]\displaystyle{ \nu }[/math] is the Poisson's ratio, and [math]\displaystyle{ \mu }[/math] is the shear modulus.

[math]\displaystyle{ \varphi }[/math]	[math]\displaystyle{ 2~\mu~u_r\, }[/math]	[math]\displaystyle{ 2~\mu~u_\theta\, }[/math]
[math]\displaystyle{ r^2\, }[/math]	[math]\displaystyle{ (\kappa-1)~r }[/math]	[math]\displaystyle{ 0 }[/math]
[math]\displaystyle{ r^2~\ln r }[/math]	[math]\displaystyle{ (\kappa-1)~r~\ln r - r }[/math]	[math]\displaystyle{ (\kappa + 1)~r~\theta }[/math]
[math]\displaystyle{ \ln r\, }[/math]	[math]\displaystyle{ -r^{-1}\, }[/math]	[math]\displaystyle{ 0 }[/math]
[math]\displaystyle{ \theta\, }[/math]	[math]\displaystyle{ 0 }[/math]	[math]\displaystyle{ -r^{-1}\, }[/math]
[math]\displaystyle{ r^3~\cos\theta \, }[/math]	[math]\displaystyle{ (\kappa-2)~r^2~\cos\theta \, }[/math]	[math]\displaystyle{ (\kappa+2)~r^2~\sin\theta \, }[/math]
[math]\displaystyle{ r\theta~\cos\theta \, }[/math]	[math]\displaystyle{ \frac{1}{2}[(\kappa-1) \theta~\cos\theta + \{1 - (\kappa+1) \ln r\} ~\sin\theta]\, }[/math]	[math]\displaystyle{ -\frac{1}{2}[(\kappa-1) \theta~\sin\theta + \{1 + (\kappa+1) \ln r\} ~\cos\theta]\, }[/math]
[math]\displaystyle{ r~\ln r~\cos\theta \, }[/math]	[math]\displaystyle{ \frac{1}{2}[(\kappa+1) \theta~\sin\theta - \{1 - (\kappa-1) \ln r\} ~\cos\theta] \, }[/math]	[math]\displaystyle{ \frac{1}{2}[(\kappa+1) \theta~\cos\theta - \{1 + (\kappa-1) \ln r\} ~\sin\theta] \, }[/math]
[math]\displaystyle{ r^{-1}~\cos\theta \, }[/math]	[math]\displaystyle{ r^{-2}~\cos\theta \, }[/math]	[math]\displaystyle{ r^{-2}~\sin\theta \, }[/math]
[math]\displaystyle{ r^3~\sin\theta \, }[/math]	[math]\displaystyle{ (\kappa-2)~r^2~\sin\theta \, }[/math]	[math]\displaystyle{ -(\kappa+2)~r^2~\cos\theta \, }[/math]
[math]\displaystyle{ r\theta~\sin\theta \, }[/math]	[math]\displaystyle{ \frac{1}{2}[(\kappa-1) \theta~\sin\theta - \{1 - (\kappa+1) \ln r\} ~\cos\theta]\, }[/math]	[math]\displaystyle{ \frac{1}{2}[(\kappa-1) \theta~\cos\theta - \{1 + (\kappa+1) \ln r\} ~\sin\theta]\, }[/math]
[math]\displaystyle{ r~\ln r~\sin\theta \, }[/math]	[math]\displaystyle{ -\frac{1}{2}[(\kappa+1) \theta~\cos\theta + \{1 - (\kappa-1) \ln r\} ~\sin\theta] \, }[/math]	[math]\displaystyle{ \frac{1}{2}[(\kappa+1) \theta~\sin\theta + \{1 + (\kappa-1) \ln r\} ~\cos\theta] \, }[/math]
[math]\displaystyle{ r^{-1}~\sin\theta \, }[/math]	[math]\displaystyle{ r^{-2}~\sin\theta \, }[/math]	[math]\displaystyle{ -r^{-2}~\cos\theta \, }[/math]
[math]\displaystyle{ r^{n+2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ (\kappa-n-1)~r^{n+1}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ (\kappa+n+1)~r^{n+1}~\sin(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n+2}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ (\kappa+n-1)~r^{-n+1}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -(\kappa-n+1)~r^{-n+1}~\sin(n\theta)\, }[/math]
[math]\displaystyle{ r^n~\cos(n\theta) \, }[/math]	[math]\displaystyle{ -n~r^{n-1}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ n~r^{n-1}~\sin(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ n~r^{-n-1}~\cos(n\theta) \, }[/math]	[math]\displaystyle{ n(~r^{-n-1}~\sin(n\theta) \, }[/math]
[math]\displaystyle{ r^{n+2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ (\kappa-n-1)~r^{n+1}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -(\kappa+n+1)~r^{n+1}~\cos(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n+2}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ (\kappa+n-1)~r^{-n+1}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ (\kappa-n+1)~r^{-n+1}~\cos(n\theta)\, }[/math]
[math]\displaystyle{ r^n~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -n~r^{n-1}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -n~r^{n-1}~\cos(n\theta) \, }[/math]
[math]\displaystyle{ r^{-n}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ n~r^{-n-1}~\sin(n\theta) \, }[/math]	[math]\displaystyle{ -n~r^{-n-1}~\cos(n\theta) \, }[/math]

Note that a rigid body displacement can be superposed on the Michell solution of the form

[math]\displaystyle{ \begin{align} u_r &= A~\cos\theta + B~\sin\theta \\ u_\theta &= -A~\sin\theta + B~\cos\theta + C~r\\ \end{align} }[/math]

to obtain an admissible displacement field.

See also

• Linear elasticity
• Flamant solution
• John Henry Michell