When the rod is thick Euler's theory fails to decribe properly the deformation, in these case we refer to Timoshenko beam theory. In this case, the governing equations are the following coupled system of ordinary differential equations:
$$
\frac{d^2}{d x^2}\left(EI\frac{d \varphi}{d x}\right) = q(x) \\
\frac{d w}{d x} = \varphi - \frac{1}{\kappa AG} \frac{d}{d x}\left(EI\frac{d \varphi}{d x}\right)
$$
where
- $L$ is the length of the beam.
- $A$ is the cross section area.
- $E$ is the elastic modulus.
- $G$ is the shear modulus.
- $I$ is the second moment of area.
- $\kappa$, called the Timoshenko shear coefficient, depends on the geometry ($\kappa = 5/6$ for a rectangular section).
- $q(x)$ is a distributed load (force per length).
Timoshenko beam theory for the static case is equivalent to the Euler–Bernoulli beam theory when the last term above is neglected, an approximation that is valid when
$\frac{3EI}{\kappa L^2 A G} \ll 1$, that is when the beam is sufficiently slender.
Combining the two equations gives, for a homogeneous beam of constant cross-section, the differential equation
$$ EI\frac{d^4 w}{d x^4} = q(x) - \frac{EI}{\kappa A G}\frac{d^2 q}{d x^2}$$
which leads to a different shape and a larger deflection with respect to the approximate Euler's theory.
(credits)