Energy Methods: Strain Energy in Members

Understanding how structures store energy as they deform is crucial for analyzing their strength, stiffness, and stability. Strain energy is the elastic energy stored within a material when it is deformed under load; it's formally defined as the work done by external forces in deforming a body, provided no energy is lost as heat. This concept forms the backbone of powerful energy methods, which allow engineers to solve complex deflection and load problems by considering the balance of energy in a system.

The Fundamental Concept: Strain Energy and Work

When you apply a force to a deformable body, you do work. If the material is linear elastic (obeying Hooke's Law) and the loading is applied gradually, this work is entirely stored internally as recoverable strain energy. For a simple axial member, like a spring or a rod, the force-deformation relationship is a straight line. The work done, and thus the strain energy stored ($U$), is the area under the force-deformation curve. This area is a triangle, leading to the fundamental relationship:

$$U=\frac{1}{2}P\delta$$

where $P$ is the applied force and $\delta$ is the corresponding deformation. This "one-half" factor is critical and only holds for linear elastic behavior. A more general approach considers energy stored per unit volume, known as strain energy density ($u$). For a general state of stress, it integrates the stress-strain curve. For a uniaxial, linear-elastic condition, it simplifies to $u=\frac{1}{2}\sigma ϵ$, where $\sigma$ is normal stress and $ϵ$ is normal strain. The total strain energy in a member is found by integrating this density over the entire volume: $U={\int }_{V}udV$.

Strain Energy for Axial Loading

Consider a prismatic bar of length $L$, cross-sectional area $A$, and modulus of elasticity $E$, subjected to a gradual axial load $P$. The normal stress is $\sigma =P/A$ (assuming constant), and the strain is $ϵ=\sigma /E=P/(AE)$. The deformation is $\delta =PL/(AE)$. Using the basic formula, $U=\frac{1}{2}P\delta =P^{2}L/(2AE)$. To express this in terms of the internal force (which is constant and equal to $P$ in this simple case), we write the general form for axial loading:

$$U={\int }_0^L\frac{N(x{)}^2}{2EA}dx$$

Here, $N(x)$ represents the internal axial force as a function of position along the member's length. You must use this integral form when the internal force or cross-sectional properties vary along the length. The derivation comes directly from integrating the strain energy density: $U={\int }_V\frac{{\sigma }^2}{2E}dV={\int }_0^L{\int }_A\frac{(N/A{)}^2}{2E}dAdx={\int }_0^L\frac{N^2}{2EA}dx$.

Strain Energy for Torsional Loading

For a circular shaft subjected to torque, the primary stress is shear stress. For a shaft of length $L$ and polar moment of inertia $J$, under a constant torque $T$, the angle of twist is $\varphi =TL/(GJ)$, where $G$ is the shear modulus. The work done by the gradually applied torque is $U=\frac{1}{2}T\varphi =T^{2}L/(2GJ)$. In its general integral form, accounting for varying internal torque $T(x)$, the expression is:

$$U={\int }_0^L\frac{T(x{)}^2}{2GJ}dx$$

This equation is analogous to the axial case, with internal torque $T$ replacing internal force $N$, shear modulus $G$ replacing elastic modulus $E$, and polar moment of inertia $J$ replacing area $A$. It's derived from the shear form of strain energy density, $u={\tau }^2/(2G)$, where $\tau$ is the shear stress, and integrating over the volume.

Strain Energy for Bending (Flexural Loading)

Bending members, like beams, store energy primarily due to flexural stress (normal stress from bending). For a beam segment in pure bending with moment $M$, the normal stress varies linearly from the neutral axis: $\sigma =-My/I$, where $I$ is the area moment of inertia and $y$ is the distance from the neutral axis. Substituting into the strain energy density formula and integrating over the cross-sectional area $A$ first yields the energy per unit length: $dU=\frac{M^2}{2EI}dx$. Therefore, the total strain energy for bending is:

$$U={\int }_0^L\frac{M(x{)}^2}{2EI}dx$$

Here, $M(x)$ is the internal bending moment as a function of position $x$ along the beam. This is often the dominant energy component in beam deflection problems. Crucially, this formula applies to bending about a principal axis and neglects the energy from shear deformation, which is typically very small for slender beams.

Superposition: Total Strain Energy from Combined Loading

Real structural members are often subject to multiple load types simultaneously—axial force, torque, and bending moment may all act on a single member, such as a crankshaft or a loaded frame. Because strain energy is a scalar quantity, the total strain energy stored in the member is simply the sum of the energies from each loading mode. You find this by integrating the contributions throughout the volume. The general expression becomes:

$$U_{total}={\int }_0^L\left(\frac{N(x{)}^2}{2EA}+\frac{T(x{)}^2}{2GJ}+\frac{M(x{)}^2}{2EI}\right)dx$$

This summation principle is powerful. For example, to find the total strain energy in a cantilever beam subjected to a vertical load at its free end, you would typically only need the bending term, as $N(x)$ and $T(x)$ are zero. However, if that same beam also had an axial load applied, you would add the integral for $N(x)$. You must evaluate each internal force function independently over the length of the member and sum the results.

Common Pitfalls

1. Forgetting the "1/2" and Squared Terms: The formulas $U=\frac{1}{2}P\delta$ and $U=\int \frac{F_{int}^2}{2k}dx$ are easy to misremember. A constant mistake is writing $U=P^{2}L/(EA)$ instead of $P^{2}L/(2EA)$. Always recall the triangular area under the linear force-deformation curve to remember the 1/2 factor.
2. Misapplying Formulas to Non-Prismatic or Variable-Load Members: Using the simple $U=F^{2}L/(2k)$ form when the internal force or cross-section varies is incorrect. You must use the integral form, breaking the member into segments where the internal force function or sectional property is constant, and then summing the contributions.
3. Incorrectly Superposing Loads Before Squaring: Strain energy is not linear with load because the internal forces are squared. You cannot calculate strain energy for each load separately using the total internal force function and then add. Instead, you must find the combined internal force function $N(x)$, $M(x)$, etc., from all loads acting together, square that function, and then integrate. Superposition of energy is only valid if the internal forces from different loads are uncoupled.
4. Neglecting to Use Consistent Units: This is a fundamental but frequent error in calculations. Ensure that force (N, lb), length (m, in), modulus (Pa, psi), and inertia ($m^4$, $i{n}^4$) are in a consistent system. Mixing units will lead to orders-of-magnitude errors in your energy result.

Summary

• Strain energy ($U$) is the elastic energy stored in a deformed body and equals the work done by gradually applied loads for linear-elastic materials, given by $U=\frac{1}{2}\times \text{Force}\times \text{Deformation}$.
• The internal strain energy for axial loading is $U=\int \frac{N(x{)}^2}{2EA}dx$, for torsional loading is $U=\int \frac{T(x{)}^2}{2GJ}dx$, and for bending is $U=\int \frac{M(x{)}^2}{2EI}dx$.
• These formulas derive from integrating the strain energy density ($u=\frac{1}{2}\sigma ϵ$ or $\frac{1}{2}\tau \gamma$) over the member's volume.
• For members under combined loading, the total strain energy is the sum of the contributions from each stress resultant: $U_{total}=\int \left(\frac{N^2}{2EA}+\frac{T^2}{2GJ}+\frac{M^2}{2EI}\right)dx$.
• Always use the integral form when internal forces or cross-sectional properties vary along the length, and be meticulous about the squared terms and the 1/2 coefficient in all calculations.