Elasticity/Energy methods example 3

Example 3 : Torsion of a cylinder
Suppose that the cross-section of the cylinder is bounded by the curves:
 * $$\begin{align}

x_1 & = c \\ x_1 & = d \\ x_2 & = g(x_1) \\ x_2 & = -g(x_1) \end{align}$$

A statically admissible Prandtl stress function for this cross-section is

\phi = [x_2^2 - g^2(x_1)] f(x_1) $$ with the restrictions that $$f(x_1)\,$$ is twice continuously differentiable and $$f(c) = f(d) = 0\,$$. We seek to derive the best approximate Prandtl stress function of this form by minimizing the complementary energy.

You can show that the complementary energy per unit length of the cylinder can be expressed as

\Pi^c[\phi(x_1,x_2)] = \int_{\mathcal S} \left[\frac{1}{2\mu}(\phi_{1}^2   + \phi_{,2}^2) - 2\alpha\phi\right]~dA $$ Plugging in the given form of $$\phi$$ and after some algebra, we get

\Pi^c[f(x_1)] = \frac{15}{16\mu} \int_c^d \left[\frac{1}{2} g^5(f^{'})^2 + \frac{5}{2}g^4 g^{'} f f^{'} + \frac{15}{4}g^3 (g^{'} f)^2 + \frac{5}{4}g^3 f^2 + \frac{5}{2}\mu\alpha g^3 f\right]~dx_1 $$

Taking the variation of $$\Pi$$, and after considerable manipulation, we get

\delta\Pi^c = \frac{15}{16\mu} \int_c^d \frac{g^5}{4}\left[f^{} - \frac{5}{g^2}(g^2 f)^{} + \frac{10}{g^2} f + \frac{10\mu\alpha}{g^2}\right]\delta f~dx_1 $$

Now, if we consider the cross-section to be rectangular, then we have $$c = -a$$, $$d = a$$, and $$g = b$$. Therefore, the above equation reduces to

\delta\Pi^c = \frac{15}{16\mu} \int_{-a}^a \frac{b^5}{4}\left[f^{} - 5f^{} + \frac{10}{b^2} f + \frac{10\mu\alpha}{b^2}\right]\delta f~dx_1 $$ Therefore, the function $$f$$ that minimizes $$\Pi^c$$ satisfies the equation

- 4f^{''} + \frac{10}{b^2} f + \frac{10\mu\alpha}{b^2} = 0 $$ or,

f^{''} - \frac{5}{2b^2} f - \frac{5\mu\alpha}{2b^2} = 0 $$ with the static admissibility conditions $$f(-a) = f(a) = 0\,$$.

The general solution of the above equation is

f = A\sinh\left(\sqrt{2.5}~\frac{x_1}{b}\right) + B\cosh\left(\sqrt{2.5}~\frac{x_1}{b}\right) - \mu\alpha $$ Therefore, from the BCs,

A = 0; B = \frac{\mu\alpha}{\cosh{\sqrt{2.5}~a/b}} $$

We can substitute back into $$\phi$$ to get the approximate Prandtl stress function for this problem. The error between the exact and this approximate solution is generally less than 1%.