Cantilever beam
Click here to go back to Solid Mechanics
On this page we conduct numerical studies of the cantilever beam, a common numerical benchmark in elastostatics.
Exact solution
Consider a beam od dimensions $L \times D$ with the origin of the coordinate system at $(x,y) = (0,D/2)$. The left end of the beam is fixed and the right end of the beam is loaded with a vertical traction: \begin{equation} t_y(y) = \frac{P}{2I}\left(\frac{D^2}{4}-y^2\right). \end{equation} The exact solution of this problem is available in (Timoshenko and Goodier, 1977) and can be summarized as follows. The displacements in $x$- and $y-$ direction are given by \begin{align*} u_x &= -\frac{P}{6EI}\left(\left(6L-3x\right)x + \left(2+\nu\right)\left(y^2-\frac{D^2}{4}\right)\right) \\ u_y &= \frac{P}{6EI}\left(3\nu y^2(L-x) + (4+5\nu)\frac{D^2 x}{4} + (3L-x)x^2\right) \end{align*} where $E$ is Youngs modulus, $\nu$ is the Poisson ratio and $I$ is the moment of inertia for a beam with rectangular cross section and unit thickness: \[ I = \frac{D^3}{12}. \]
Numerical solution
For the numerical solution we first choose the following parameters:
- Loading: $P = -1000$ N
- Young's modulus: $E = 3 \times 10^7$ N/m2
