Roymech engineering encyclopedia

Strain Energy Methods and Castiglianos Theorem


Strain energy methods provide an efficient means of calculating specific deflections and rotations in beams and structures. Castigliano's theorem relates displacement to the derivative of the total strain energy and can often avoid the need to derive the complete elastic curve of a member.

Need help choosing a beam calculation? See the Beam Design and Analysis Guide for guidance on beam types, loading, bending, shear, deflection and selecting the appropriate RoyMech calculation method.

This page introduces the application of Castigliano's theorem and includes worked examples for a simply supported beam, a cantilever and a simple frame. For the strain-energy expressions used here, see Strain Energy.

Contents

Energy Methods for Beam Deflection

When external forces are applied to an elastic beam or structure, the structure deforms and the external forces perform work. For a conservative linear-elastic system loaded gradually from zero, this work is stored as elastic strain energy in the structure.

Under these assumptions:

W = U

where:

  • W = external work done
  • U = elastic strain energy stored in the structure

When the loads are removed within the elastic range, the deformation is recovered and the stored elastic strain energy is released.

This energy relationship provides an alternative route for determining deflections and rotations. Rather than deriving the complete deflected shape of a beam, an energy method can be used to calculate the displacement at a particular point and in a particular direction.

Castigliano's Theorem

Castigliano's theorem provides a convenient method for finding a particular displacement or rotation in an elastic structure. Unlike direct integration of the beam differential equation, it does not normally provide an expression for displacement at every point along the member.

For a linearly elastic structure undergoing small displacements, the displacement corresponding to an applied force is obtained by differentiating the total strain energy with respect to that force.

Linear Displacement

Castigliano theorem equation for linear displacement

δi = ∂ U ∂ Fi

where δi is the displacement at the point of application of force Fi, measured in the direction of that force.

Angular Displacement

The corresponding relationship for rotation is:

Castigliano theorem equation for angular displacement

θi = ∂ U ∂ Mi

where θi is the rotation corresponding to the applied moment Mi.

Displacement Where No Load Exists

Castigliano's theorem is also useful when the required displacement is at a point where no real force or moment acts.

A fictitious force is introduced at the required point and in the direction of the required displacement. The total strain energy is written including this fictitious load, differentiated with respect to it, and the fictitious load is then set equal to zero.

Similarly, a fictitious moment can be introduced when a rotation is required at a point where no real moment acts.

Procedure for Applying Castigliano's Theorem

Linear Deflection

To determine a linear displacement δf in the direction of a real or fictitious force Ff:

  1. Determine the internal force and moment expressions throughout the structure, including the real loads and the fictitious force where required.
  2. Form an expression for the total elastic strain energy U. Depending on the problem, this may include strain energy due to axial force, bending, torsion and shear.
  3. Differentiate the total strain energy with respect to the force:
δf = ∂ U ∂ Ff
  1. If Ff is fictitious, set Ff = 0 after differentiation and evaluate the resulting expression.

Angular Deflection

To determine a rotation θf in the direction of a real or fictitious moment Mf, use the same procedure, but differentiate the total strain energy with respect to the moment:

θf = ∂ U ∂ Mf

If the moment is fictitious, set Mf = 0 after differentiation.

Example 1 – Simply Supported Beam

Consider a simply supported rectangular beam carrying a central force F. The vertical deflection at the load point is required.

Span, l 2 m
Width, b 0.1 m
Depth, h 0.05 m
Force, F 10,000 N
Young's modulus, E 206 GPa
Shear modulus, G 78.610 GPa
Second moment of area, I 4.17 × 10-6 m4

Simply supported beam carrying a central load for Castigliano deflection example

The beam is rectangular and the strain energy due to both bending and transverse shear is included. Because the beam and loading are symmetrical, the total strain energy can be obtained by integrating over half the beam and doubling the result.

For the left-hand half of the beam:

M = F2 x V = F2

Strain Energy

Using the bending and transverse-shear strain-energy relationships developed on the Strain Energy page:

Total bending and transverse shear strain energy for simply supported beam

The image is retained so that the original derivation can be checked before it is removed or modified.

Deflection from Castigliano's Theorem

The deflection at the point of application of F is:

δF = ∂ U ∂ F

Castigliano calculation of central beam deflection including bending and transverse shear

For the dimensions and elastic properties used in this example, the shear component of deflection is small compared with the bending component. The original calculation gives a transverse-shear contribution of approximately 0.8% of the bending deflection.

Transverse-shear deflection is often small for slender homogeneous beams, but it should not be assumed negligible for short or deep beams, materials with relatively low shear modulus, or structural forms in which shear deformation is significant.

Example 2 – Cantilever with a Fictitious Load

Consider the cantilever shown below. A real force F acts at a distance b from the free end. The deflection at the free end is required.

Because no real load acts at the free end in the required direction, a fictitious force Ff is introduced there.

Cantilever beam with applied load and fictitious free-end load for Castigliano analysis

For this example, only bending strain energy is included.

Bending-Moment Expressions

The bending moment must be expressed for the appropriate regions of the cantilever, including both the real and fictitious forces.

Bending moment expressions for cantilever with real and fictitious loads

Strain Energy

Using the bending strain-energy relationship:

U = ∫ 0 L M2 2EI dx

the strain energy for the separate regions of the cantilever is:

Bending strain energy expression for cantilever Castigliano example

Differentiate with Respect to the Fictitious Force

From Castigliano's theorem:

δf = ∂ U ∂ Ff

Differentiation of cantilever strain energy with respect to fictitious load

Set the Fictitious Force to Zero

After differentiation, set:

Ff = 0

and integrate over the appropriate beam lengths.

Final free-end deflection calculation for cantilever using Castigliano theorem

Example 3 – Frame Deflection

Consider the frame shown below. One leg is free to move horizontally and the other is pinned. A horizontal force F is applied to the free leg, and the horizontal displacement at that point is required.

The original example neglects strain energy due to transverse shear.

Frame subjected to horizontal force for Castigliano displacement calculation

Internal Forces

For the members identified in the diagram, the original example uses the following internal-force relationships:

MCD = Fy MBC = F L1 MAB = Fy NBC = F

The force in member BC is axial, so its axial strain energy must be included in addition to the bending strain energy of the frame members.

Total Strain Energy

Total bending and axial strain energy for frame Castigliano example

Horizontal Deflection

The displacement in the direction of the applied force is:

δ = ∂ U ∂ F

Castigliano calculation of horizontal frame displacement

Application of Energy Methods

Castigliano's theorem is particularly useful when the displacement or rotation is required at a specific point rather than as a complete equation for the deflected shape. It can be applied to structures in which the total strain energy contains contributions from bending, axial loading, torsion and shear.

The energy terms included in a calculation should reflect the deformation mechanisms that are significant for the structure being analysed. For many slender beams bending dominates, while shear, axial or torsional deformation may become important in other geometries.

The theorem in the form used on this page assumes elastic behaviour and the conditions appropriate to the strain-energy formulation. For nonlinear material behaviour, large deformation or other non-linear effects, the applicable energy relationship must be considered separately.

Related Beam References