# Elasticity/Disk with hole
## Disk with a central hole

* An elastic disk with a central circular hole
Under general loading, for the stresses and displacements to be single-valued and continuous, they must be periodic in $\theta$, e.g., $\sigma _{11}(r,\theta )=\sigma _{11}(r,\theta +2m\pi )$.
An Airy stress function appropriate from this situation is
    ${\text{(83)}}\qquad \varphi =\sum _{n=0}^{\infty }f_{n}(r)\cos(n\theta )+\sum _{n=0}^{\infty }g_{n}(r)\sin(n\theta )$
In the absence of body forces,
    ${\text{(84)}}\qquad \nabla ^{4}{\varphi }=\nabla ^{2}{(\nabla ^{2}{\varphi })}=0~;~~\nabla ^{2}=\left({\cfrac {\partial ^{2}}{\partial r^{2}}}+{\cfrac {1}{r}}{\cfrac {\partial }{\partial r}}+{\cfrac {1}{r^{2}}}{\cfrac {\partial ^{2}}{\partial \theta ^{2}}}\right)$
Plug in $\varphi$.
    ${\begin{aligned}\nabla ^{2}{\varphi }=&\sum _{n=0}^{\infty }\left[f_{n}^{''}(r)\cos(n\theta )+{\cfrac {1}{r}}f_{n}^{'}(r)\cos(n\theta )-{\cfrac {n^{2}}{r^{2}}}f_{n}(r)\cos(n\theta )\right]+\\&\sum _{n=0}^{\infty }\left[g_{n}^{''}(r)\sin(n\theta )+{\cfrac {1}{r}}g_{n}^{'}(r)\sin(n\theta )-{\cfrac {n^{2}}{r^{2}}}g_{n}(r)\sin(n\theta )\right]\qquad {\text{(85)}}\end{aligned}}$
or,
    ${\text{(86)}}\qquad \nabla ^{2}{\varphi }=\sum _{n=0}^{\infty }F_{n}(r)\cos(n\theta )+\sum _{n=0}^{\infty }G_{n}(r)\sin(n\theta )$
Therefore,
    ${\begin{aligned}\nabla ^{4}{\varphi }=&\sum _{n=0}^{\infty }\left[F_{n}^{''}(r)\cos(n\theta )+{\cfrac {1}{r}}F_{n}^{'}(r)\cos(n\theta )-{\cfrac {n^{2}}{r^{2}}}F_{n}(r)\cos(n\theta )\right]+\\&\sum _{n=0}^{\infty }\left[G_{n}^{''}(r)\sin(n\theta )+{\cfrac {1}{r}}G_{n}^{'}(r)\sin(n\theta )-{\cfrac {n^{2}}{r^{2}}}G_{n}(r)\sin(n\theta )\right]\qquad {\text{(87)}}\end{aligned}}$
To satisfy the compatibility condition $\nabla ^{4}{\varphi }=0$, we need
    ${\begin{aligned}{\text{(88)}}\qquad F_{n}^{''}(r)+{\cfrac {1}{r}}F_{n}^{'}(r)-{\cfrac {n^{2}}{r^{2}}}F_{n}(r)&=0\\{\text{(89)}}\qquad G_{n}^{''}(r)+{\cfrac {1}{r}}G_{n}^{'}(r)-{\cfrac {n^{2}}{r^{2}}}G_{n}(r)&=0\end{aligned}}$
The general solution of these Euler-Cauchy type equations is
    ${\begin{aligned}{\text{(90)}}\qquad F_{n}(r)&=A_{1}r^{n}+B_{1}r^{-n}\\{\text{(91)}}\qquad G_{n}(r)&=C_{1}r^{n}+D_{1}r^{-n}\end{aligned}}$
We can use either to determine $f_{n}(r)$. Thus,
    ${\text{(92)}}\qquad f_{n}^{''}(r)+{\cfrac {1}{r}}f_{n}^{'}(r)-{\cfrac {n^{2}}{r^{2}}}f_{n}(r)=A_{1}r^{n}+B_{1}r^{-n}$
or,
    ${\text{(93)}}\qquad r^{2}f_{n}^{''}(r)+rf_{n}^{'}(r)-n^{2}f_{n}(r)=A_{1}r^{n+2}+B_{1}r^{-n+2}$
The homogeneous and particular solutions of this equation are
    ${\begin{aligned}{\text{(94)}}\qquad f_{n}^{h}(r)&=A_{2}r^{n}+B_{2}r^{-n}\\{\text{(95)}}\qquad f_{n}^{p}(r)&=A_{1}r^{n+2}+B_{1}r^{-n+2}\end{aligned}}$
Hence, the general solution is
    ${\text{(96)}}\qquad f_{n}(r)=A_{1}r^{n+2}+B_{1}r^{-n+2}+A_{2}r^{n}+B_{2}r^{-n}$
This form is valid for $n>1$. If $n=0,1$, alternative forms are obtained. Thus,
    ${\begin{aligned}{\text{(97)}}\qquad f_{0}(r)&=A_{O}r^{2}+B_{0}r^{2}\ln r+C_{0}+D_{0}\ln r\\{\text{(98)}}\qquad f_{1}(r)&=A_{1}r^{3}+B_{1}r+C_{1}r\ln r+D_{1}r^{-1}\\{\text{(99)}}\qquad f_{n}(r)&=A_{n}r^{n+2}+B_{n}r^{n}+C_{n}r^{-n+2}+D_{n}r^{-n}~,~~n>1\end{aligned}}$
Terms in $f_{n}$ are chosen according to the specific problem of interest.

### Traction BCs

    at $r=a$
    ${\text{(100)}}\qquad \sigma _{rr}=T_{1}(\theta )~,~~\sigma _{r\theta }=T_{2}(\theta )$
    at $r=b$
    ${\text{(101)}}\qquad \sigma _{rr}=T_{3}(\theta )~,~~\sigma _{r\theta }=T_{4}(\theta )$
Express $T_{i}(\theta )$ in Fourier series form.
    ${\text{(102)}}\qquad T_{i}(\theta )=\sum _{n=0}^{\infty }A_{ni}\cos(n\theta )+\sum _{n=0}^{\infty }B_{ni}\sin(n\theta )~,~~i=1,2,3,4$
Terms in $T_{i}$ are chosen according to the specific problem of interest.