Learning topic
Castigliano's Theorem and Mohr's Integral
Calculate beam and frame deflections with Castigliano’s theorem, the unit-load Mohr integral, and Vereshchagin’s diagram multiplication rule.
Castigliano's theorem and the Mohr integral make it possible to determine individual linear and angular displacements of elastic member systems using strain energy or products of internal force resultants from two loading states.
Castigliano's theorem
For a linear-elastic system, the displacement in the direction of a generalized force is obtained by differentiating the strain energy with respect to that force:
$$\delta_i=\frac{\partial U}{\partial P_i}.$$
Similarly, the rotation at a generalized applied moment is:
$$\varphi_i=\frac{\partial U}{\partial M_i}.$$
- $U$ — strain energy;
- $P_i$ — generalized force;
- $\delta_i$ — corresponding linear displacement;
- $M_i$ — generalized moment;
- $\varphi_i$ — corresponding rotation.
If the required load is absent from the actual system, an auxiliary force or moment may be introduced in the required direction, $U$ written as a function of that parameter, differentiated, and the auxiliary parameter then set to zero.
For a beam in which bending strain energy dominates, $U=\int M^2/(2EI)\,dx$. If $M$ depends on force $P$, differentiation gives:
$$\delta_P=\frac{\partial U}{\partial P}=\int\frac{M}{EI}\frac{\partial M}{\partial P}\,dx.$$
For a linear system, $\partial M/\partial P$ corresponds to the bending-moment diagram produced by a unit force in the direction of $P$, showing the connection between Castigliano's theorem and the unit-load method.
Auxiliary force or moment
If no real force acts at the point of the required displacement, introduce an auxiliary parameter $X$ in the required direction. Write the internal force resultants as functions of $X$, evaluate $\partial U/\partial X$, and set $X=0$ after differentiation.
Mohr integral
The displacement at a specified point and direction can be determined by the unit-load method. For bending of a beam or frame:
$$\delta=\int_0^L\frac{M(x)\,\bar M(x)}{EI}\,dx,$$
where $M(x)$ is the bending moment from the real loading and $\bar M(x)$ is the moment from a unit force applied at the point and in the direction of the required linear displacement.
For a required rotation, apply a unit moment.
For a general member, applicable contributions may be summed:
$$\delta=\int\frac{N\bar N}{EA}\,dx+\int\frac{M\bar M}{EI}\,dx+\int\frac{T\bar T}{GJ_p}\,dx+\int\frac{Q\bar Q}{\kappa GA}\,dx,$$
- $N$, $M$, $T$, $Q$ — internal force resultants from the real loading;
- $\bar N$, $\bar M$, $\bar T$, $\bar Q$ — corresponding resultants from the unit-load state;
- $E$, $G$, $A$, $I$, $J_p$, $\kappa$ — stiffness and section parameters.
The sign of $\delta$ indicates whether the actual displacement agrees with the direction of the introduced unit load.
For a linear displacement, apply a unit force in the direction of the required displacement. For a rotation, apply a unit moment. The sign of the result indicates whether the actual displacement agrees with the direction of the unit load.
Unit-load procedure
- Analyze the system under the real loading and construct the required internal-force diagrams.
- Remove the real loading and apply a unit force or moment at the point and in the direction of the required displacement.
- Construct the unit-load diagrams.
- Write the Mohr integral over all segments and members.
- Use the actual $EA$, $EI$, $GJ_p$ and, where necessary, shear stiffness.
- Evaluate the integrals and sum all contributions.
Vereshchagin rule
For a segment with constant $EI$, the Mohr integral $\int M\bar M/(EI)\,dx$ can be evaluated graphically when one of the two multiplied diagrams is linear over that segment.
- Construct the real-load bending-moment diagram $M$ and the unit-load diagram $\bar M$ for the required displacement.
- Divide the beam or frame into segments where $EI$ is constant and the diagrams have convenient geometric shapes.
- For each simple part of one diagram, determine its algebraic area $\Omega$ and the coordinate of its centroid.
- On the other diagram, which must be linear over that segment, determine the ordinate $\bar M_c$ beneath the centroid of the first diagram area.
- Calculate the contribution $\Delta=\Omega\bar M_c/(EI)$, retaining the signs of the diagrams.
- Sum the contributions from all segments. A positive result acts in the direction of the introduced unit load; a negative result acts in the opposite direction.
If neither diagram is linear on a segment, the simple area-times-centroid-ordinate rule is generally not applicable without additional subdivision or direct integration.
The Vereshchagin rule is especially efficient for beams and frames with constant-$EI$ segments and simple diagrams. It replaces part of the analytical integration with operations involving diagram areas and ordinates.
Reciprocity
For a linearly elastic system, the work of one set of forces through the displacements caused by another set equals the reciprocal work:
$$\sum_i P_i^{(1)}\delta_i^{(2)}=\sum_i P_i^{(2)}\delta_i^{(1)}.$$
As a special case for two unit forces, Maxwell's reciprocal-displacement theorem gives:
$$\delta_{ij}=\delta_{ji}.$$
- $P_i^{(1)}$, $P_i^{(2)}$ — generalized forces in loading states 1 and 2;
- $\delta_i^{(1)}$, $\delta_i^{(2)}$ — corresponding displacements;
- $\delta_{ij}$ — displacement at coordinate $i$ caused by a unit force at coordinate $j$.
Thus the displacement at point $i$ in the direction of force $i$ caused by a unit force at point $j$ equals the corresponding displacement at point $j$ caused by a unit force at point $i$, provided the system is linearly elastic and reciprocity conditions are satisfied.
Reciprocity theorems provide a useful check on unit-load calculations and can sometimes suggest a simpler reciprocal loading state.
Common errors
Typical mistakes include choosing the wrong direction for the unit load, omitting members or segments, mixing diagram signs, using one $EI$ where stiffness changes, and mechanically applying the Vereshchagin rule to two nonlinear diagrams.
About this topic
Mohr's integral evaluates structural displacements from arbitrary load configurations by applying dummy unit forces or unit moments. This section covers unit load method integration and Vereshchagin's visual diagram multiplication method for beams and frames using centroidal ordinates.