Login
OR
Don't worry. We'll fix this for you!
Energy methods provide powerful tools for analyzing problems in Mechanics of Solids. Unlike classical force-balance methods, which rely heavily on equilibrium equations and can become cumbersome for complex structures, energy methods exploit the fundamental principle of conservation of energy to determine deflections, stresses, and strains efficiently.
At their core, these methods involve two key ideas:
Using these concepts, energy methods relate forces and displacements through energy expressions, enabling elegant solutions to problems involving beams, bars, shafts, and trusses.
Discovering deflections by energy methods is particularly advantageous when dealing with statically indeterminate structures, where classical equilibrium alone is insufficient.
To begin, let's understand the nature of work and energy in deformable solids.
When a force \( P \) causes a displacement \( \delta \) at its point of application, the work done by this force is
Work done \( W = P \times \delta \)
This work is stored as strain energy inside the material as it deforms elastically.
Strain energy is the energy stored per unit volume in the material due to deformation under loading. It represents the area under the stress-strain curve up to the current deformation.
This energy can be categorized according to the type of loading:
1. Axial Loading: Consider a bar of length \( L \), cross-sectional area \( A \), subjected to axial load \( P \). The strain energy stored is
\[ U = \frac{P^2 L}{2 A E} \]where \( E \) is the modulus of elasticity of the material.2. Bending: For beams under bending moment \( M(x) \) varying along the length, strain energy is
\[ U = \int_0^L \frac{M(x)^2}{2 E I} \, dx \]where \( I \) is the moment of inertia of the beam's cross-section.3. Torsion: A circular shaft of length \( L \) under torque \( T \) stores strain energy
\[ U = \frac{T^2 L}{2 G J} \]where \( G \) is the shear modulus and \( J \) is the polar moment of inertia.Understanding these expressions is crucial as they form the basis of energy methods used to evaluate deflections and rotations.
Castigliano's theorems provide a direct way to find deflections and rotations in elastic structures by differentiating strain energy with respect to loads or moments.
This theorem states:
If the strain energy \( U \) of a linear elastic structure can be expressed as a function of applied forces \( F_i \), then the partial derivative of \( U \) with respect to any force \( F_i \) gives the displacement \( \delta_i \) in the direction of that force.
Mathematically,
\[\delta_i = \frac{\partial U}{\partial F_i}\]Similarly, if the forces are constant and moments \( M_i \) act, then partial derivatives of \( U \) with respect to moments give angular rotations.
For rotation angle \( \theta_i \),
\[\theta_i = \frac{\partial U}{\partial M_i}\]These theorems allow transformations from energy stored due to loads directly into deflections - a very powerful approach in structural analysis.
In the above diagram, the free end deflection \( \delta \) due to load \( P \) can be found by first expressing strain energy \( U \) in terms of \( P \), then differentiating.
The Unit Load Method is a practical application of Castigliano's theorem used particularly for beam deflections. It is well-suited for both statically determinate and certain indeterminate structures.
The method follows these key steps:
graph TD A[Apply Real Loads] --> B[Determine Internal Forces (M, V, N)] B --> C[Apply Unit Load at Point of Desired Deflection] C --> D[Compute Internal Forces due to Unit Load (M₁, V₁, N₁)] D --> E[Calculate Deflection using Energy Expression and Integration] E --> F[Apply Boundary Conditions and Solve]
Step-by-step, you determine the internal force distributions due to the actual load and the artificial unit load separately and use their products integrated over the structure length to find deflections.
This avoids direct displacement calculations and leverages energy expressions for simplified analysis.
Step 1: Write expression for bending moment \( M(x) \) at distance \( x \) from fixed end (0 ≤ x ≤ L):
\( M(x) = -P (L - x) \)
(Negative sign indicates sagging moment)
Step 2: Write strain energy in bending:
\( U = \int_0^L \frac{M(x)^2}{2 E I} dx = \int_0^L \frac{P^2 (L - x)^2}{2 E I} dx \)
Step 3: Evaluate integral:
\[ U = \frac{P^2}{2 E I} \int_0^L (L - x)^2 dx = \frac{P^2}{2 E I} \left[ \frac{(L - x)^3}{3} \right]_0^L = \frac{P^2}{2 E I} \frac{L^3}{3} \]
Step 4: Apply Castigliano's first theorem:
\[ \delta = \frac{\partial U}{\partial P} = \frac{\partial}{\partial P} \left( \frac{P^2 L^3}{6 E I} \right) = \frac{2 P L^3}{6 E I} = \frac{P L^3}{3 E I} \]
Step 5: Substitute values:
\( \delta = \frac{1000 \times (2)^3}{3 \times 2 \times 10^{11} \times 4 \times 10^{-6}} = \frac{1000 \times 8}{3 \times 2 \times 10^{11} \times 4 \times 10^{-6}} \)
\( = \frac{8000}{2.4 \times 10^{6}} = 0.00333\,m = 3.33\,mm \)
Answer: Deflection at free end is approximately 3.33 mm.
Step 1: Calculate bending moment due to real load \( w \):
\( M(x) = \frac{w x}{2}(L - x) \), for \( 0 \leq x \leq L \)
Step 2: Apply unit load at midpoint, i.e. at \( x = L/2 = 2\, m \). Bending moment due to unit load:
\( M_1(x) = \left\{ \begin{array}{ll} \frac{x}{2} & 0 \leq x \leq 2 \\ \frac{L - x}{2} & 2 \leq x \leq 4 \end{array} \right. \)
Step 3: Calculate deflection using:
\[ \delta = \frac{1}{E I} \int_0^L M(x) M_1(x) dx \]
Split integral at \( x=2 \):
\[ \delta = \frac{1}{E I} \left[ \int_0^2 M(x) M_1(x) dx + \int_2^4 M(x) M_1(x) dx \right] \]
Step 4: Integrate (detailed algebra omitted for brevity):
Substituting \( w = 2000\, N/m \), \( L = 4 \), and integrating yields:
\[ \delta = \frac{23.99}{E I} \]
Step 5: Calculate numeric value:
\[ \delta = \frac{23.99}{2 \times 10^{11} \times 5 \times 10^{-6}} = 2.399 \times 10^{-5} \, m = 0.024 \, mm \]
Answer: Deflection at beam midpoint is approximately 0.024 mm.
Step 1: Calculate polar moment of inertia \( J \) for shaft:
\[ J = \frac{\pi d^4}{32} = \frac{\pi (0.05)^4}{32} = 3.07 \times 10^{-7} \, m^4 \]
Step 2: Calculate strain energy stored:
\[ U = \frac{T^2 L}{2 G J} = \frac{(500)^2 \times 1.5}{2 \times 8 \times 10^{10} \times 3.07 \times 10^{-7}} = 76.3 \, J \]
Step 3: Calculate angle of twist:
\[ \theta = \frac{T L}{G J} = \frac{500 \times 1.5}{8 \times 10^{10} \times 3.07 \times 10^{-7}} = 0.0306 \, rad = 1.75^\circ \]
Answer: Strain energy stored = 76.3 J, angle of twist = 1.75°.
Step 1: Calculate cross-sectional area:
\[ A = \frac{\pi d^2}{4} = \frac{\pi (0.02)^2}{4} = 3.14 \times 10^{-4} \, m^2 \]
Step 2: Calculate strain energy:
\[ U = \frac{P^2 L}{2 A E} = \frac{(10,000)^2 \times 3}{2 \times 3.14 \times 10^{-4} \times 2 \times 10^{11}} = 2.4\, J \]
Answer: The strain energy stored is 2.4 joules.
Step 1: Calculate strain energy due to axial load:
\[ U_{axial} = \frac{P^2 L}{2 A E} = \frac{(20,000)^2 \times 3}{2 \times 6 \times 10^{-4} \times 2 \times 10^{11}} = 10\, J \]
Step 2: Calculate strain energy due to bending:
\[ U_{bending} = \int_0^3 \frac{M(x)^2}{2 E I} dx = \int_0^3 \frac{(200 x)^2}{2 \times 2 \times 10^{11} \times 8 \times 10^{-6}} dx \]
\[ = \frac{1}{3.2 \times 10^6} \int_0^3 40,000 x^2 dx = \frac{40,000}{3.2 \times 10^6} \times \left[ \frac{x^3}{3} \right]_0^3 = \frac{40,000 \times 27}{3.2 \times 10^6 \times 3} = 0.1125\, J \]
Step 3: Total strain energy:
\[ U_{total} = U_{axial} + U_{bending} = 10 + 0.1125 = 10.1125 \, J \]
Answer: Total strain energy stored is approximately 10.11 J.
When to use: Applying Castigliano's theorems or unit load method.
When to use: Solving deflections in symmetric loading or beam geometry.
When to use: Exam preparations and solving energy method problems.
When to use: Beams or shafts under axial, bending, and torsional loads simultaneously.
When to use: When direct integration is complex or time is limited.
Progress tracking is paywalled — subscribe to mark subtopics as understood and save your streak.
Go to practice →
Your Score
Your Rank
| Rank | Name | Score |
|---|
