Timoshenko-Ehrenfest beam theory


The Timoshenko-Ehrenfest beam theory was developed by Stephen Timoshenko and Paul Ehrenfest early in the 20th century. The model takes into account shear deformation and rotational bending effects, making it suitable for describing the behaviour of thick beams, sandwich composite beams, or beams subject to high-frequency excitation when the wavelength approaches the thickness of the beam. The resulting equation is of 4th order but, unlike Euler–Bernoulli beam theory, there is also a second-order partial derivative present. Physically, taking into account the added mechanisms of deformation effectively lowers the stiffness of the beam, while the result is a larger deflection under a static load and lower predicted eigenfrequencies for a given set of boundary conditions. The latter effect is more noticeable for higher frequencies as the wavelength becomes shorter, and thus the distance between opposing shear forces decreases.
Rotary inertia effect was introduced by Bresse and Rayleigh.
If the shear modulus of the beam material approaches infinity—and thus the beam becomes rigid in shear—and if rotational inertia effects are neglected, Timoshenko beam theory converges towards ordinary beam theory.

Quasistatic Timoshenko beam

In static Timoshenko beam theory without axial effects, the displacements of the beam are assumed to be given by
where are the coordinates of a point in the beam, are the components of the displacement vector in the three coordinate directions, is the angle of rotation of the normal to the mid-surface of the beam, and is the displacement of the mid-surface in the -direction.
The governing equations are the following coupled system of ordinary differential equations:
The Timoshenko beam theory for the static case is equivalent to the Euler-Bernoulli theory when the last term above is neglected, an approximation that is valid when
where
Combining the two equations gives, for a homogeneous beam of constant cross-section,
The bending moment and the shear force in the beam are related to the displacement and the rotation. These relations, for a linear elastic Timoshenko beam, are:
Then, from the strain-displacement relations for small strains, the non-zero strains based on the Timoshenko assumptions are
Since the actual shear strain in the beam is not constant over the cross section we introduce a correction factor such that
The variation in the internal energy of the beam is
Define
Then
Integration by parts, and noting that because of the boundary conditions the variations are zero at the ends of the beam, leads to
The variation in the external work done on the beam by a transverse load per unit length is
Then, for a quasistatic beam, the principle of virtual work gives
The governing equations for the beam are, from the fundamental theorem of variational calculus,
For a linear elastic beam
Therefore the governing equations for the beam may be expressed as
Combining the two equations together gives

Boundary conditions

The two equations that describe the deformation of a Timoshenko beam have to be augmented with boundary conditions if they are to be solved. Four boundary conditions are needed for the problem to be well-posed. Typical boundary conditions are:
For a cantilever beam, one boundary is clamped while the other is free. Let us use a right handed coordinate system where the direction is positive towards right and the direction is positive upward. Following normal convention, we assume that positive forces act in the positive directions of the and axes and positive moments act in the clockwise direction. We also assume that the sign convention of the stress resultants is such that positive bending moments compress the material at the bottom of the beam and positive shear forces rotate the beam in a counterclockwise direction.
Let us assume that the clamped end is at and the free end is at. If a point load is applied to the free end in the positive direction, a free body diagram of the beam gives us
and
Therefore, from the expressions for the bending moment and shear force, we have
Integration of the first equation, and application of the boundary condition at, leads to
The second equation can then be written as
Integration and application of the boundary condition at gives
The axial stress is given by

Dynamic Timoshenko beam

In Timoshenko beam theory without axial effects, the displacements of the beam are assumed to be given by
where are the coordinates of a point in the beam, are the components of the displacement vector in the three coordinate directions, is the angle of rotation of the normal to the mid-surface of the beam, and is the displacement of the mid-surface in the -direction.
Starting from the above assumption, the Timoshenko beam theory, allowing for vibrations, may be described with the coupled linear partial differential equations:
where the dependent variables are, the translational displacement of the beam, and, the angular displacement. Note that unlike the Euler-Bernoulli theory, the angular deflection is another variable and not approximated by the slope of the deflection. Also,
These parameters are not necessarily constants.
For a linear elastic, isotropic, homogeneous beam of constant cross-section these two equations can be combined to give
From equation, assuming appropriate smoothness, we have
Differentiating equation gives
Substituting equation,, into equation and rearrange, we get
The Timoshenko equation predicts a critical frequency
For normal modes the Timoshenko equation can be solved. Being a fourth order equation, there are four independent solutions, two oscillatory and two evanescent for frequencies below.
For frequencies larger than all solutions are oscillatory and, as consequence, a second spectrum appears.

Axial effects

If the displacements of the beam are given by
where is an additional displacement in the -direction, then the governing equations of a Timoshenko beam take the form
where and is an externally applied axial force. Any external axial force is balanced by the stress resultant
where is the axial stress and the thickness of the beam has been assumed to be.
The combined beam equation with axial force effects included is

Damping

If, in addition to axial forces, we assume a damping force that is proportional to the velocity with the form
the coupled governing equations for a Timoshenko beam take the form
and the combined equation becomes
A caveat to this Ansatz damping force is that, whereas viscosity leads to a frequency-dependent and amplitude-independent damping rate of beam oscillations, the empirically measured damping rates are frequency-insensitive, but depend on the amplitude of beam deflection.

Shear coefficient

Determining the shear coefficient is not straightforward ; generally it must satisfy:
The shear coefficient depends on Poisson's ratio. The attempts to provide precise expressions were made by many scientists, including Stephen Timoshenko, Raymond D. Mindlin, G. R. Cowper, N. G. Stephen, J. R. Hutchinson etc.. In engineering practice, the expressions by Stephen Timoshenko are sufficient in most cases. In 1975 Kaneko published an excellent review of studies of the shear coefficient. More recently new experimental data show that the shear coefficient is underestimated.
According to Cowper for solid rectangular cross-sections,
and for solid circular cross-sections,
where is Poisson's ratio.