В задаче есть две различные компоненты напряжения: радиальная (сила давления действует в направлении радиус-вектора) и тангенциальная (в плоскости, перпендикулярной радиус-вектору). Уравнения для них записываются как:
\[\left\{\begin{array}lE\varepsilon_r=\sigma_r-2\mu\sigma_\tau\\E\varepsilon_\tau=\sigma_\tau-\mu\sigma_r-\mu\sigma_\tau\end{array}\right.\implies\sigma_\tau=\frac1{1-\mu}\left(E\varepsilon_\tau+\mu\sigma_r\right)\implies\\\implies E\varepsilon_r=\sigma_r-\frac{2\mu}{1-\mu}\left(E\varepsilon_\tau+\mu\sigma_r\right)=\frac{(1-2\mu)(1+\mu)}{1-\mu}\sigma_r-\frac{2\mu}{1-\mu}E\varepsilon_\tau.\]Нам осталось разобраться с $\sigma_r$ и $\varepsilon_{r,\tau}$. Для $\varepsilon_{r,\tau}$ <<очевидно, что>>:\[\varepsilon_r=\frac{\mathrm du}{\mathrm dr},\quad\varepsilon_\tau=\frac ur.\]С $\sigma_r$ всё сложнее. Рассмотрим, какие силы действуют на участок $[r,r+\mathrm dr]$, видимый под телесным углом $\mathrm d\alpha\times\mathrm d\beta$. Помимо объёмной силы тяжести $\left(\frac{4\pi\rho G}3r\cdot\rho\,r\mathrm d\alpha\,r\mathrm d\beta\,\mathrm dr\right)$, в условие механического равновесия надо включить силу, действующую на <<боковую поверхность>> этого участка $\left(-2\sigma_\tau r\mathrm dr\mathrm d\alpha\mathrm d\beta\right)$ и силы, действующие на <<крышку>> $(-\sigma_r(r+\mathrm dr)(r+\mathrm dr)^2\mathrm d\alpha\mathrm d\beta)$ и <<дно>> $(-\sigma_r(r)r^2\mathrm d\alpha\mathrm d\beta)$:\[-\frac{4\pi\rho^2G}3r^3+\frac{\mathrm d}{\mathrm dr}\left[r^2\sigma_r\right]=2r\sigma_\tau=\frac{2r}{1-\mu}\left(E\varepsilon_\tau+\mu\sigma_r\right).\]Глядя на полученные уравнения, хочется найти $u$ и $\sigma_r$ в виде каких-то не очень сложных многочленов. Учитывая, что $\sigma_r(R)=0$, $u(0)=0$, напрашивается поиск решения в виде $u(r)=-r(A-Br^2)$, $\sigma_r=-C(R^2-r^2)$. Отсюда получаем:\[\left\{\begin{array}{l}
E \frac{\mathrm d}{\mathrm d r}\left[r\left(A-B r^2\right)\right]=C \frac{(1-2 \mu)(1+\mu)}{1-\mu}\left(R^2-r^2\right)-\frac{2 \mu}{1-\mu} E\left(A-B r^2\right) \\
\frac{4 \pi g^2 G}{3} r^3+C \frac{\mathrm d}{\mathrm d r}\left[r^2\left(R^2-r^2\right)\right]=\frac{2 r}{1-\mu}\left[E\left(A-B r^2\right)+\mu C\left(R^2-r^2\right)\right]
\end{array}\right.\implies\\\implies\left\{\begin{array}{lc}A=&\frac{2\pi\rho^2R^2(1-2 \mu)(3-\mu)}{15E(1-\mu)}\\B=&\frac{2\pi\rho^2(1-2 \mu)(1+\mu)}{15E(1-\mu)}\\C=&\frac{2\pi\rho^2G(3-\mu)}{15(1-\mu)}\end{array}\right.\implies u(r)=-\frac{2\pi\rho^2R^2(1-2 \mu)(1+\mu)}{15E(1-\mu)} r\left[\frac{3-\mu}{1+\mu}-\frac{r^2}{R^2}\right]\]