Chapter 9
Beams
9.1 Stresses
Assume a rectangular beam in plane stress or plane strain subject to bending. (The “plane” direction is the \(y\)-direction.) The surface normal of the beam is oriented in \(z\)-direction and \(z=0\) is in the middle of the beam. Bending will give rise to internal stresses inside of the beam. We will require that these stresses comply with the external shear force \(Q(x)\) and bending moment \(M(x)\) in the weak or integral sense,
Because of the plane state, integration in \(y\)-direction will only yield a constant factor, the width \(t\) of the beam.
We can derive a condition equivalent to static equilibrium, \(\sigma _{ij,j}=0\), for the weak quantities defined in Eqs. \eqref{eq:beamweakforce} and \eqref{eq:beamweakmoment}. Taking the derivative of \(Q(x)\) yields
We call the quantity \(q(x)\equiv t \left [\sigma _{zz}(x, h/2) - \sigma _{zz}(x, -h/2)\right ]\) the line load of the beam. Next, we take the derivative of \(M(x)\),
The first term on the right hand side vanishes because the surfaces are traction free, \(\tau _{xz}=0\) at \(z=h/2\) and \(z=-h/2\). This yields
for the weak form of the equilibrium conditions.
Note: We have not yet talked about free surfaces, but of course every physical object will be bounded by surfaces. The surface can be loaded, i.e. there can a force acting on it. This force, normalized by area, is called a tractions. (We have already encountered a type of traction in form of the line load in the first chapters.) The traction is a vector field (defined over the surface of our object) with components \(t_x\), \(t_y\) and \(t_z\). Note that in the literature these are often denoted by \(P\) for the normal component, and \(Q_x\), \(Q_y\) for the components perpendicular to the surface. (Attention: Do not confuse this with the \(Q\) that we call the shear force here!)
These traction enter our equations as boundary condition. They fix three components of the stress tensor at the surface. For a surface with a normal in \(z\)-direction,
We now assume that the stress is a linear function of the position \(z\) perpendicular to the beam axis,
In what follow we show that this is a good assumption, i.e. we can fulfill force and moment equilibrium and the resulting strains fulfill the compatibility conditions.
Note: The theory derived in this chapter is commonly referred to as the Euler-Bernoulli beam theory. The starting point of this theory is typically not Eq. \eqref{eq:beamstress1}, but the assumption that each cross section will remain plane and undergo small rotations during deformation. These are sometimes called the Bernoulli assumptions. They imply that the strain \(\varepsilon _{xx}\propto z\) rather than the stress. It is often argued that \(\sigma _{xx}=E \varepsilon _{xx}\) but this of course ignores the other components of the strain tensor. As will be seen below, the Bernoulli assumptions are actually wrong but assuming a linear stress profile leads to the correct small strain expression for the deformation of a beam.
It is straightforward to compute the bending moment,
where \(A\) is the cross-section of the beam and
is the axial moment of inertia. For a rectangular beam of height \(h\) and width \(t\), \(I_y=h^3 t/12\). The moment of inertia is a geometric factor and depends on the shape of the cross-section of the beam.
We can rewrite Eq. \eqref{eq:beamstress1} as
Note that an additional longitudinal force \(L\) will simple be an additive contribution to Eq. \eqref{eq:beamstressxx}, \(\sigma _{xx}(x,z) = M(x)/I_y z + L/A\).
We can now use the condition for static equilibrium to compute the full stress tensor \(\t {\sigma }\). From \(\sigma _{xx,x} + \tau _{xz,z}=0\) we obtain
We can integrate this across the height \(h\) of the beam, keeping in mind that \(\tau _{xz}=0\) at a traction-free surface, to yield
Next, we use \(\sigma _{zz,z}+\tau _{xz,x}=0\) to obtain
We need to integrate this equation again, but now \(\sigma _{zz}\not =0\) at the surface since the beam is subject to a line load \(q(x)\). Rather, we need the condition that the loads on top and bottom surface of the beam balance, \(\sigma _{zz}(x,h/2)=-\sigma _{zz}(x,-h/2)\). This gives
where \(t\) is the width of the beam.
Note that in this derivation, we have required that the stress gives rise to a certain bending moment through Eq. \eqref{eq:beammoment}. Since we do not prescribe the specific stress state \(\sigma _{xx}(x, z)\) but only its integral, this is a weak condition. Similarly, integration of Eq. \eqref{eq:beamstressxz} in accordance with Eq. \eqref{eq:beamweakforce} gives
i.e. the force acting on the cross-section at position \(x\) along the beam. Again, this condition is fulfilled in the integral, i.e. in the weak sense.
9.2 Displacements
Now that we know the stress inside the beam, we can compute the displacement and thereby the deformation of the beam from Hooke’s law. Starting from Hooke’s law (in plain stress, i.e. for \(\sigma _{yy}=0\)),
and taking the derivative of Eq. \eqref{eq:beamhooke2} with respect to \(x\), we obtain
Combining these two equations yields
where we now insert Eqs. \eqref{eq:beamstressxx}, \eqref{eq:beamstressxz} and \eqref{eq:beamstresszz}. This gives
which at the surface of the beam \(w(x)=u_z(x,h/2)\) becomes
This equation is called the Euler-Bernoulli beam equation. By using \(M_{,xx}=-q(x)\) we can write this as
in terms of the line load \(q(x)\), or
if \(E I_y\) does not depend on position \(x\).