Strain Energy Method Tutorial

Energy stored equals work done — and the 1/2 factor that everyone forgets.

Mechanics of MaterialsMechanical Engineering Year 2Free preview
⏱️ About 16 min

Slowly pull a spring and it pushes back; release it and it snaps, returning every joule you put in. Where did that energy hide while you were pulling? It was stored inside the material itself, as <em>strain energy</em>. A loaded bar, a twisted shaft, and a bent beam are all, in this view, springs — they absorb work as they deform and can hand it back. Looking at structures through the lens of energy is a quiet revolution: instead of solving force balance at every section, you equate the work the load does to the energy the member stores, and out pop deflections, impact factors, and whole classes of problems that resist a direct attack. This final module begins with the single idea on which all of it rests.

💡
The big idea: As a load is applied gradually from zero to its full value <em>P</em>, the work it does is <strong>W = ½ P δ</strong> — the area of the triangle under the load-deflection line. For a linear-elastic member that work is stored as <strong>strain energy</strong> U, so W = U. For a bar in axial tension the strain energy is U = P²L/(2AE); for a shaft in torsion U = T²L/(2GJ); for a beam in bending U = ∫ M²/(2EI) dx. The factor of ½ appears everywhere because the load builds up from zero, and it is the key to the work-energy method for finding deflections under a single load.
🎯 By the end, you'll be able to
  • Explain why strain energy carries a 1/2 factor (load ramps from zero)
  • Write the strain-energy formula for axial, torsion, and bending members
  • Apply the work-energy principle W = U to find a deflection
  • Relate strain-energy density u = sigma^2/(2E) to the total energy U

Work done by a gradually applied load

Suppose a load grows smoothly from zero up to some final value P, deflecting the member by δ along the way. The work done is not P × δ. At the instant the deflection was only δ/2, the load was only P/2; at δ/4 it was P/4. Averaged over the whole motion, only half the final load was acting, so the work is

W = ½ P δ

Geometrically, W is the area of the triangle under the straight load-deflection line — which is exactly half of the bounding rectangle P × δ. This is the 1/2 that students drop most often. Get it right here and every energy formula in the module falls into place.

🔑 Strain energy is NOT force x displacement

The single most common error in energy methods is writing U = P δ. It is half of that. The load does not act at its full value throughout the deflection — it ramps up from zero, so the work (and the stored energy) carries a factor of ½.

  • Gradually applied load: W = ½ P δ. Energy = triangular area.
  • If you ever see U = P δ (no ½), that is the work of a load held constant at P while the displacement happens — the suddenly-applied-then-held case, not the gradual one. We meet it next lesson, where it produces the famous impact factor of 2.
A load-deflection plot. A straight line rises linearly from the origin to the working point at deflection delta and load P. The strain energy U is the shaded triangle under the line, equal to one-half times P times delta. A faint rectangle behind it represents the full force times displacement product: the strain energy is exactly half of that rectangle, because the load ramps up from zero rather than acting at full value throughout the deflection. force x displacement (whole rectangle) delta deflection P load U = (1/2) P delta strain energy = shaded area

A load-deflection plot with a straight line from the origin to the working point (delta, P). The shaded triangle under the line is the strain energy U, equal to one-half P delta. A faint rectangle behind it shows the full force times displacement product; the energy is exactly half of it.

Strain energy is the area under the load-deflection line — a triangle of height P and base delta. The energy is half of the force-times-displacement rectangle.

Strain-energy density

Zoom in to a single particle of material. Its share of the energy per unit volume — the strain-energy density u — is the area under its stress-strain curve. For a linear-elastic material (σ = Eε) that area is a triangle too:

\[ u = \tfrac{1}{2}\,\sigma\varepsilon = \frac{\sigma^{2}}{2E} = \tfrac{1}{2}E\varepsilon^{2} \]
Strain-energy density (energy per unit volume) for a linear-elastic material. The total strain energy is U = integral of u over the volume. The 1/2 is the same factor as the work triangle.

The three loading forms

Substituting the relevant stress into the density and integrating over the volume gives three workhorse formulas — one for each simple loading. Memorise the pattern, not the derivation: each is (load squared) x (length) divided by (2 x stiffness).

\[ U_{\text{axial}} = \frac{P^{2}L}{2AE}, \qquad U_{\text{torsion}} = \frac{T^{2}L}{2GJ}, \qquad U_{\text{bending}} = \int_{0}^{L}\frac{M^{2}}{2EI}\,dx \]
Strain energy for a uniform bar in axial tension, a circular shaft in torsion, and a beam in bending. Each is quadratic in the load and inversely proportional to a stiffness (AE, GJ, or EI); the 1/2 is present in all three.
✨ Energy is quadratic in the load

Notice that every strain-energy formula is proportional to the square of the load. Double the load and you store four times the energy; halve it and the energy drops to a quarter. This steep dependence is the root cause of impact being so destructive: a modest increase in load dumps a great deal more energy into a member, and that energy has to go somewhere — usually into larger stress and deflection. The interactive simulator below lets you feel the P-squared relationship directly.

📝 Worked example: A steel rod (E = 200 GPa) has a solid circular cross-section of diameter d = 25 mm and length L = 1.5 m. It carries an axial tensile load P = 40 kN applied gradually. Find (a) the strain energy stored, and (b) use the work-energy principle W = U to recover the rod's elongation.
  1. Cross-sectional area: A = pi d^2 / 4 = pi (25)^2 / 4 = 490.9 mm^2.
  2. Use N and mm: E = 200 000 N/mm^2, L = 1500 mm, P = 40 000 N.
  3. (a) U = P^2 L / (2 A E) = (40 000)^2 (1500) / (2 x 490.9 x 200 000).
  4. Numerator = 1.6 x 10^9 x 1500 = 2.40 x 10^12. Denominator = 1.964 x 10^8.
  5. U = 1.222 x 10^4 N.mm = 12.2 J.
  6. (b) Work-energy: W = U, and W = (1/2) P delta, so delta = 2 U / P.
  7. delta = 2 x 12 222 / 40 000 = 0.611 mm.
  8. Check against the direct formula delta = P L / (A E) = 40 000 x 1500 / (490.9 x 200 000) = 0.611 mm. The two agree, as they must.
✓ (a) U = 12.2 J. (b) From W = U, delta = 2U/P = 0.611 mm, matching the direct deflection formula.
🎮 Interactive: strain energy vs. load LIVE
Predict first: Hold length L = 1.5 m and diameter d = 25 mm fixed, then double the load P from 20 to 40 kN. The strain energy QUADRUPLES — that is the P-squared dependence at work. Now shrink d from 25 to 15 mm at fixed load: a thinner rod stores MORE energy (it stretches more for the same force), so U rises.

An interactive calculator with sliders for axial load P, rod length L, and diameter d, and a magnitude bar showing the strain energy U in joules.

Strain energy U = P^2 L / (2 A E) for a steel rod. Drag the load to feel the quadratic dependence; thinner rods store more energy for the same load.
✏️ Practice: A solid circular steel shaft (G = 80 GPa) has diameter d = 40 mm and length L = 1.2 m. It carries a torque T = 500 N·m. Using U = T²L/(2GJ) with J = πd⁴/32, find the strain energy stored, in J. (Express T in N·mm: 500 N·m = 500 000 N·mm.)
J
Solution
  1. J = pi d^4 / 32 = pi (40)^4 / 32 = pi x 2.560 x 10^6 / 32 = 251 327 mm^4.
  2. T = 500 000 N.mm, so T^2 = 2.50 x 10^11.
  3. U = T^2 L / (2 G J) = (2.50 x 10^11)(1200) / (2 x 80 000 x 251 327).
  4. Numerator = 3.00 x 10^14. Denominator = 4.021 x 10^10.
  5. U = 7460 N.mm = 7.46 J.
✏️ Practice: A steel cantilever beam (E = 200 GPa) of length L = 1.0 m carries an end load P = 2 kN. Its rectangular cross-section is b = 20 mm wide and d = 40 mm deep. The bending strain energy is U = ∫ M²/(2EI) dx = P²L³/(6EI) (with M = Px measured from the tip and I = bd³/12). Find U, in J.
J
Solution
  1. I = b d^3 / 12 = 20 x (40)^3 / 12 = 20 x 64 000 / 12 = 106 667 mm^4.
  2. P = 2 000 N, L = 1000 mm. P^2 = 4.0 x 10^6, L^3 = 1.0 x 10^9.
  3. U = P^2 L^3 / (6 E I) = (4.0 x 10^6)(1.0 x 10^9) / (6 x 200 000 x 106 667).
  4. Numerator = 4.0 x 10^15. Denominator = 1.280 x 10^11.
  5. U = 31 250 N.mm = 31.2 J.

Check your understanding

1. The work done by a load applied gradually from zero to P, producing deflection delta, is:
Because the load ramps up from zero, its average value over the deflection is P/2, giving W = (1/2) P delta — the area of the triangle under the load-deflection line.
2. The axial strain energy of a bar is proportional to:
U = P^2 L / (2 A E) is quadratic in P: doubling the load quadruples the stored energy.
3. The work-energy principle W = U for a single gradually applied load lets you find:
Setting the external work (1/2) P delta equal to the internal strain energy U gives one equation in the single unknown deflection delta.
✅ Key takeaways
  • A gradually applied load does work W = (1/2) P delta — half of force times displacement, because the load ramps from zero.
  • That work is stored as strain energy U; for linear-elastic members W = U.
  • Axial U = P^2 L / (2 A E); torsion U = T^2 L / (2 G J); bending U = integral M^2 / (2 E I) dx.
  • Strain-energy density u = sigma^2 / (2 E); total U integrates u over the volume.
  • Setting W = U recovers a deflection under a single load — the work-energy method.
➡️ We derived everything for a load applied gradually. But many real loads arrive suddenly — a weight dropped onto a beam, a hammer blow, a vehicle hitting a buffer. Energy balance still holds, yet the result is dramatically larger. The next lesson derives the impact factor and shows why a load that is simply released — not even dropped from a height — already doubles the static deflection.
Want to test yourself on this? Try the Mechanical Aptitude test →
🎓 Go deeper: external courses & trusted references