Simply supported beam—Castigliano's method - All this
And now it’s all this
I just said what I said and it was wrong
Or was taken wrong
Next post<br>Previous post
Simply supported beam—Castigliano's method
June 1, 2026 at 8:06 AM by Dr. Drang
Continuing our series on the many ways to get the center deflection of a uniformly loaded beam, we come to another energy-based technique: Castigliano’s method.
Castigliano’s second theorem provides a relationship between displacements, forces, and the strain energy in a linearly elastic system. The strain energy, U, is the potential energy of the system due to its deformation. For a generic elastic body, like this,
Castigliano’s second theorem can be as:
If an elastic system is mounted so that rigid-body displacements of<br>the entire system are impossible and certain external point forces<br>P1, P2, … act on the system, in addition to distributed loads and thermal<br>strains, the displacement component δi of the point of application of<br>force Pi in the direction of force Pi is determined by the equation
δi=∂U∂Pi
I took this from Henry Langhaar’s Energy Methods in Applied Mechanics, a favorite text of mine. I’ve altered the variable names to avoid conflicting with some variable names we’ve used in previous posts in this series.
Looking this over and thinking about it in terms of our problem,
two concerns come to mind:
Our problem doesn’t have a point force at the center of the beam, where we want to determine the deflection.
Even if we did have a point force at the center of the beam, how do we write an expression for the strain energy in terms of that force?
Let’s tackle the second concern first. Recall from our two Rayleigh-Ritz posts (Fourier and polynomial) that the potential energy associated with beam bending is
∫0L12EI(y″)2dx
This is the strain energy, U. Unfortunately, it’s written in terms of deflection, not force, but we can fix that. We saw early on that the curvature, y″, and the bending moment, M, are proportional:
M=−EIy″ory″=−MEI
We substitute the second of these into the expression for U to get
U=∫0LM22EIdx
Because the bending moment can be written in terms of the applied loads, we’ve solved the second concern.
We solve the first concern with a trick. We can pretend there’s a load, P, at the center of the beam, and write out an expression for the bending moment, M, that includes P.
We then do the integration to get U, take its derivative with respect to P, and—here’s the trick—evaluate the derivative for P=0. This may seem like cheating, but it’s perfectly legitimate. Imagine P is very very small—it won’t contribute much to the deflection. So why not just make it zero?
Let’s go ahead and see what happens. The upward support reaction force at end of the beam above is (wL+P)/2, so we can draw a free-body diagram for the left portion of the beam like this:
where xL/2. The bending moment in the left half of the beam is therefore
M=12[(wL+P)x−wx2]
and its square is
M2=14[(wL+P)2x2−2w(wL+P)x3+w2x4]
This formula applies only in the left half of the beam, but because of symmetry we know that
U=∫0LM22EIdx=2∫0L/2M22EIdx=∫0L/2M2EIdx
So
U=14EI∫0L/2[(wL+P)2x2−2w(wL+P)x3+w2x4]dx=14EI[13(wL+P)2x3−12w(wL+P)x4+15w2x5]0L/2=14EI[L324(wL+P)2−L432w(wL+P)+L5160w2]
Taking the derivative with respect to P gives us
∂U∂P=14EI[L312(wL+P)−L432w]
Therefore the downward deflection at the center of the beam is
δ=∂U∂P|P=0=wL448EI−wL4128EI=5wL4384EI
as we’ve now seen many times.
Since this is based on Castigliano’s second theorem, you may be wondering what kinds of problems we use his first theorem to solve. It’s a good question, but one I can’t answer. In the 45 years since I first learned of Castigliano’s theorems, I’ve never used his first one.
Next post<br>Previous post
Site search
-->
Meta
drdrang at leancrew
Blog archive
RSS feed
JSON feed
Twitter stream<br>--><br>Mastodon
GitHub repositories
Recent posts
Credits
This work is licensed under a<br>Creative Commons Attribution-Share Alike 3.0 Unported License.
© 2005–2026, Dr. Drang