Energy Methods in Structural Analysis

Introduction

Structural analysis is commonly introduced through equilibrium equations, compatibility conditions, and force–deformation laws. Energy methods provide a complementary and often more efficient viewpoint: instead of determining every internal force first, the analyst relates applied work to stored strain energy and obtains selected displacements, rotations, reactions, or redundant forces directly.

These methods are especially useful for statically determinate and indeterminate beams, frames, and trusses. They are also fundamental in computational mechanics, vibration analysis, stability theory, and impact design. The principal methods developed in this chapter are:

• strain-energy calculations;
• Castigliano’s theorem;
• the unit-load or virtual-work method; and
• impact-loading analysis.

The methods apply most directly to linearly elastic members undergoing small deformation, although the underlying energy concepts extend to nonlinear systems with suitable modifications.

Learning Objectives

After studying this chapter, the student should be able to:

1. Define external work, complementary work, strain energy, and resilience.
2. Derive and use strain-energy expressions for axial, bending, torsional, and shear deformation.
3. Apply Castigliano’s theorem to calculate displacements, rotations, and redundant reactions.
4. Use virtual work and the unit-load method for beams, frames, and trusses.
5. Analyze gradually applied, suddenly applied, and impact loads.
6. Select an appropriate energy method and assess its assumptions and limitations.
7. Interpret results for practical design and structural safety.

1. Energy Concepts in Elastic Structures

When a load is applied gradually to an elastic structure, it does external work. Part of this work is stored internally as strain energy. If no energy is dissipated, the work done by the loads equals the strain energy stored in the structure.

For a single load P increasing gradually from zero to its final value, and a corresponding displacement Δ, the external work is

where δ is the instantaneous displacement. For a linear elastic system, P(δ) = Pδ/Δ, so

Here U is strain energy, with units of joules (J = N·m). The factor ½ applies only when the load is increased gradually and the load–displacement relationship is linear.

For several loads Pi and corresponding displacements Δi,

provided the structure is linear elastic and the displacements are measured in the directions of the loads. A more general expression is

where {P} is the load vector and {Δ} is the displacement vector.

1.1 Strain energy density

The strain energy per unit volume is called strain-energy density. For uniaxial stress,

where u is energy density (J/m³), σ is normal stress, and ε is normal strain. Under linear elasticity, σ = Eε, giving

where E is Young’s modulus. These expressions are equivalent only for a linearly elastic material.

2. Strain Energy of Structural Members

The total strain energy is obtained by integrating the energy density over the member volume. Different deformation modes produce different engineering formulas.

2.1 Axial deformation

For a prismatic bar of length L, constant area A, axial force P, and modulus E,

The axial displacement is

For a variable axial force or area,

where N(x) is the internal axial force.

2.2 Bending deformation

For a beam subjected to bending moment M(x), with flexural rigidity EI,

where I is the second moment of area about the bending axis. This is usually the dominant energy term for slender beams.

2.3 Torsional deformation

For a shaft carrying torque T(x), with shear modulus G and torsional constant J,

For a circular shaft, J is the polar second moment of area. For noncircular sections, J may be a torsional constant and is not always equal to the polar second moment.

2.4 Shear deformation

The shear contribution is

where V(x) is shear force, G is shear modulus, A is area, and κ is the shear-correction factor. For a rectangular section, κ is commonly taken as approximately 5/6 in elementary beam theory. Shear energy is often negligible in slender Euler–Bernoulli beams but may be important in deep beams, short beams, sandwich panels, and Timoshenko beam models.

2.5 Combined energy

When several independent deformation modes are present,

This additive form assumes small deformation and a linear constitutive model with no significant coupling between modes. In a frame, bending generally dominates, but axial and shear terms should not be omitted automatically.

3. Castigliano’s Theorem

Castigliano’s second theorem states that, for a linearly elastic structure, the partial derivative of the total strain energy with respect to an applied force gives the displacement in the direction of that force:

Similarly, the derivative with respect to an applied moment gives the corresponding rotation:

The theorem requires a conservative, elastic response and small displacements. The total energy must include all relevant deformation modes.

If internal forces depend on an applied load, the theorem can be written as

This general form is particularly useful for beams and frames.

3.1 Displacement of a simply supported beam

Consider a simply supported beam of span L carrying a central point load P. The bending moment diagram is symmetric. For the left half, 0 ≤ x ≤ L/2,

The total bending energy is

The central deflection is therefore

This result assumes a prismatic beam, constant E and I, negligible shear deformation, and a point load applied at midspan.

3.2 Rotation from an applied moment

If a moment M₀ is applied to a structure, then

The result is a rotation in radians. Since radians are dimensionless, the units of ∂U/∂M₀ are J/(N·m), which is dimensionless, confirming unit consistency.

4. Castigliano’s Theorem for Statically Indeterminate Structures

Energy methods are particularly powerful for redundant structures. Suppose an unknown redundant reaction R is introduced. Compatibility may require the displacement at the redundant location to be zero. Castigliano’s theorem gives

If an imposed support settlement or prescribed displacement Δ₀ exists, the condition becomes

The sign convention must be maintained consistently. A redundant force is treated as an ordinary load while calculating internal-force functions.

Worked Example 1: Axial bar with a redundant support

A uniform bar of length L and area A is fixed at both ends and subjected to a tensile point load P at a location that divides the bar into lengths a and b, where a + b = L. Let the left and right support reactions be RA and RB. Equilibrium gives

The internal axial forces are constant in the two segments. Taking RB as the redundant and using compatibility, the displacement at the load point must be zero. The left segment elongation caused by its internal force and the right segment deformation must cancel. Equivalently, Castigliano’s theorem gives the force distribution. Because both segments have the same EA, the load divides in proportion to the opposite segment lengths:

Thus the force in the left segment is RA and the force in the right segment is RB. The result agrees with the physical interpretation that the shorter segment is stiffer and attracts the larger share of load. If A or E differs between segments, the load division is governed by axial stiffness EA rather than length alone.

5. Virtual Work and the Unit-Load Method

Virtual work determines a real displacement by applying a fictitious, unit load at the point and in the direction of the desired displacement. The structure’s real internal forces are denoted by N, M, T, and V, while the internal forces produced by the unit load are denoted by n, m, t, and v.

For a linearly elastic member,

For bending-dominated beams,

The unit load has magnitude 1 N when calculating a translation, or 1 N·m when calculating a rotation. The resulting displacement has units of metres or radians, respectively.

The method follows directly from the principle of virtual work: the external virtual work equals the internal virtual work for a compatible virtual displacement field. It is often more convenient than differentiating a lengthy energy expression.

5.1 Trusses

For a pin-jointed truss carrying axial forces only,

where Ni is the real member force, ni is the member force caused by a unit load, Li is member length, Ai is area, and Ei is Young’s modulus. Tension and compression are represented by positive and negative signs according to the chosen convention.

Worked Example 2: Deflection of a cantilever

A cantilever of length L carries a downward point load P at its free end. Determine the free-end deflection using virtual work. The real bending moment at a section x from the fixed end is

Apply a unit downward load at the free end. The corresponding virtual moment is

Therefore,

The negative signs multiply to give a positive downward displacement magnitude. This is the familiar Euler–Bernoulli result. If shear deformation is included, an additional term is PL/(κGA).

6. Maxwell–Betti Reciprocity

For a linear elastic structure, the displacement at point A caused by a load at point B equals the displacement at B caused by the same load applied at A. In notation,

This reciprocity theorem is useful for checking calculations, constructing flexibility matrices, and understanding structural influence coefficients. It depends on linear elasticity, small deformation, and conservative behavior. It does not generally hold in the same form for materially nonlinear, path-dependent, or dissipative systems.

7. Impact and Suddenly Applied Loads

Energy methods distinguish between gradual loading and dynamic loading. A gradually applied load allows the structure to respond quasi-statically. A suddenly applied load is imposed before the structure has time to deform, and an impact load transfers kinetic energy into structural deformation.

7.1 Suddenly applied load

Suppose a load P is applied instantaneously to a linear elastic spring of stiffness k. At maximum displacement Δmax, the external work is PΔmax and the stored energy is ½kΔmax². Since P = kΔst, where Δst is the static displacement,

so

Thus an abruptly applied load produces twice the static displacement and, in a linear elastic member, twice the static stress.

7.2 Impact of a falling mass

Let a mass W fall through height h onto a structure. Let Δst be the static deflection caused by W and Δmax the maximum impact deflection. Conservation of energy gives

Since k = W/Δst,

The impact factor relative to the static displacement is

For h = 0, IF = 2, corresponding to sudden application. As h increases, the impact factor increases. This idealization neglects damping, local contact deformation, plasticity, wave propagation, and mass distribution in the structure.

Worked Example 3: Impact on a bar

A 1,000 N weight falls 0.10 m onto a vertical elastic bar whose static deflection under the weight is 0.005 m. The maximum displacement is

The impact factor is approximately 7.40. The large value illustrates why drop loads can produce stresses many times greater than the corresponding static stress.

8. Choosing an Energy Method

The most suitable method depends on the desired result and structural model.

A good engineering workflow is to identify the deformation mode, draw the real internal-force diagrams, define the virtual system or redundant, perform the integration, and then check dimensions and limiting behavior.

9. Illustrative Case Study: Deflection of a Crane-Girder Bracket

Background

Consider an illustrative steel bracket supporting a small overhead crane. The bracket behaves approximately as a cantilever beam of length 1.2 m. A service load of 18 kN acts downward at the free end. The section has E = 200 GPa and second moment of area I = 8.0 × 10⁻⁶ m⁴. The design team must estimate the free-end deflection and determine whether shear deformation is significant.

Problem

Determine the bending deflection using the unit-load method. Assume the bracket is prismatic and initially neglect shear deformation.

Analysis

For a cantilever with an end load,

Substitution gives

Since 1.2³ = 1.728,

If the shear area is approximated as κA = 0.010 m² and G = 77 GPa, the shear contribution is

Thus the total estimated deflection is approximately 8.13 mm, and bending accounts for almost all of the deformation.

Outcome and learning points

The unit-load method separates bending and shear contributions and makes the governing deformation mechanism transparent. The bracket should nevertheless be checked for stress, local yielding, connection flexibility, fatigue, lateral stability, and serviceability limits. The calculation is not a substitute for a complete design standard check.

10. Assumptions, Limitations, and Professional Practice

Energy methods in this chapter generally assume:

• small displacements and rotations;
• linear elastic material behavior;
• known and constant material properties unless variation is included;
• idealized supports and connections;
• compatibility of deformation;
• negligible damping for basic impact calculations; and
• correct internal-force diagrams.

Common errors include using ½PΔ for a suddenly applied load, omitting axial energy in a frame, using a bending-only formula for a deep beam, confusing the polar moment with the torsional constant, and applying a unit load in the wrong direction.

For real structures, engineers must also consider geometric nonlinearity, residual stress, cracking, yielding, connection slip, vibration, fatigue, contact, and load uncertainty. Building codes and occupational safety regulations govern acceptable design practice. Energy calculations should be documented with sign conventions, units, assumptions, load combinations, and independent checks.

Modern finite-element programs use energy principles in their formulation. The principle of minimum potential energy underlies many stiffness methods, while virtual work is used in element derivation, error estimation, and adaptive mesh refinement. Software does not remove the need to verify boundary conditions, mesh adequacy, material models, and reaction equilibrium. Artificial intelligence may assist in model generation or anomaly detection, but responsibility for physical assumptions and safety remains with the engineer.

Chapter Summary

Strain energy is the work stored in an elastic structure. For linear elastic loading, it is one-half the product of load and corresponding displacement. Axial, bending, torsional, and shear energies are obtained by integrating the appropriate internal-force quantity divided by twice the relevant stiffness.

Castigliano’s theorem obtains a displacement by differentiating total strain energy with respect to an applied force, and a rotation by differentiating with respect to an applied moment. It is particularly effective for redundant reactions when compatibility conditions are imposed. The unit-load method, derived from virtual work, calculates a desired displacement using the product of real and virtual internal forces. For trusses, the summation form is especially convenient.

Impact analysis requires an energy balance rather than the gradual-loading relation ½PΔ. A suddenly applied load produces twice the static deflection in an ideal linear system, while a falling mass may generate much larger impact factors. All energy methods depend on appropriate assumptions, accurate internal-force diagrams, and careful interpretation of deformation modes.

Key Terms

Strain energy: Energy stored in a body because of elastic deformation.

External work: Work performed by applied loads as they move through displacements.

Strain-energy density: Strain energy stored per unit volume.

Castigliano’s theorem: The derivative of strain energy with respect to a load gives the corresponding displacement.

Unit-load method: A virtual-work method using a unit force or moment to calculate a displacement or rotation.

Virtual work: Work associated with a compatible imagined displacement system.

Flexural rigidity: The product EI governing bending stiffness.

Axial rigidity: The product EA governing axial stiffness.

Torsional rigidity: The product GJ governing torsional stiffness.

Shear-correction factor: A factor κ accounting for nonuniform shear stress distribution.

Redundant: An unknown force or moment that cannot be found from equilibrium alone.

Compatibility: A condition requiring connected parts of a structure to deform consistently.

Resilience: The capacity to store elastic strain energy.

Suddenly applied load: A load imposed instantaneously without initial deformation.

Impact factor: Ratio of maximum dynamic response to the corresponding static response.

Maxwell–Betti reciprocity: Equality of cross-displacements in a linear elastic structure.

Flexibility coefficient: Displacement produced at one coordinate by a unit load at another coordinate.

Conservative system: A system in which work is path independent and recoverable.

Euler–Bernoulli beam: A slender-beam model that neglects shear deformation.

Timoshenko beam: A beam model that includes shear deformation and rotational flexibility.

Review Questions

Short Answer Questions

1. Define strain energy and state its SI unit.
2. Under what conditions is U = ½PΔ valid?
3. State Castigliano’s theorem for displacement and rotation.
4. Why is a unit load applied in the virtual-work method?
5. What is the impact factor for a suddenly applied load in an ideal linear elastic system?

Long Answer Questions

1. Derive the bending strain-energy expression for a beam with variable bending moment and flexural rigidity.
2. Explain how Castigliano’s theorem can be used to determine a redundant reaction.
3. Derive the unit-load expression for the deflection of a beam, including bending and shear terms.
4. Compare the assumptions and applications of Euler–Bernoulli and Timoshenko beam models.
5. Derive the maximum-deflection equation for a falling weight using conservation of energy.

Analytical / Application Questions

1. A simply supported beam of span 4 m carries a 12 kN central load. Using the bending-energy method, derive its central deflection in terms of E and I, and evaluate it for E = 200 GPa and I = 6.0 × 10⁻⁶ m⁴.
2. A truss member carries a real axial force of 35 kN and has length 2.5 m, area 450 mm², and E = 200 GPa. A unit load at the required joint produces a member force of 0.62 kN per unit applied kilonewton. Determine that member’s contribution to joint displacement and explain the sign of the contribution if the member is in compression.
3. A 500 N mass falls 50 mm onto a spring-like structure with static deflection 2 mm under the mass. Calculate the impact factor and discuss why damping and local contact deformation may make the ideal result inaccurate.

Multiple Choice Questions with Answers

Question 1

For a linearly elastic member subjected to a gradually applied load P with final displacement Δ, the stored strain energy is:

A. PΔ

B. ½PΔ

C. 2PΔ

D. P/Δ

Correct Answer: B. ½PΔ

Question 2

Castigliano’s theorem gives a displacement by differentiating total strain energy with respect to:

A. Young’s modulus

B. cross-sectional area

C. the corresponding applied force

D. member length

Correct Answer: C. the corresponding applied force

Question 3

The bending strain energy of a beam is calculated from:

A. ∫M²/(2EI) dx

B. ∫M/(2EI) dx

C. ∫EI/M² dx

D. ∫V²/(2EA) dx

Correct Answer: A. ∫M²/(2EI) dx

Question 4

In the unit-load method, the virtual internal moment m is produced by:

A. the actual service load

B. a unit load applied at the desired displacement coordinate

C. the support reaction only

D. the material weight

Correct Answer: B. a unit load applied at the desired displacement coordinate

Question 5

For a pin-jointed truss with axial deformation only, the displacement contribution of member i is:

A. NᵢnᵢLᵢ/(AᵢEᵢ)

B. Nᵢ²/(2AᵢEᵢ)

C. MᵢmᵢLᵢ/(EI)

D. Vᵢ²Lᵢ/(κGA)

Correct Answer: A. NᵢnᵢLᵢ/(AᵢEᵢ)

Question 6

A suddenly applied load on an ideal linear elastic structure produces a maximum displacement equal to:

A. one-half the static displacement

B. the static displacement

C. twice the static displacement

D. four times the static displacement

Correct Answer: C. twice the static displacement

Question 7

For a falling weight W with drop height h and static deflection Δst, the impact factor is:

A. 1 − √(1 + 2h/Δst)

B. 1 + √(1 + 2h/Δst)

C. √(1 + h/Δst)

D. 2h/Δst

Correct Answer: B. 1 + √(1 + 2h/Δst)

Question 8

Which deformation mode is generally dominant in a slender Euler–Bernoulli beam?

A. Bending

B. Shear

C. Torsion

D. Thermal expansion

Correct Answer: A. Bending

Question 9

The compatibility equation for a redundant reaction R at a fixed support with no settlement is commonly written as:

A. ∂U/∂R = 1

B. ∂U/∂R = 0

C. U/R = 0

D. ∂R/∂U = 0

Correct Answer: B. ∂U/∂R = 0

Question 10

Maxwell–Betti reciprocity applies most directly to structures that are:

A. perfectly plastic and dissipative

B. linear elastic with small deformation

C. always dynamically unstable

D. independent of boundary conditions

Correct Answer: B. linear elastic with small deformation