Unit 8 Strain Energy: Structure
Unit 8 Strain Energy: Structure
Unit 8 Strain Energy: Structure
8.7 Summary
8.8 Answers to SAQs
8.1 INTRODUCTION
This unit introduces you to the concept of strain energy, how it is related to stress and
strain and the method of estimating this energy in different modes of application of forces
or loads on the body. It also explains to you how to obtain the stresses in the case of
sudden and impact loads using the strain energy concept. Further, it deals with
Castigliano’s theorem which can be employed to estimate the defections of deformable
bodies under loads.
Objectives
After studying this unit, you should be able to
• conceptualise the potential energy,
• estimate the strain energy under different types of loads,
• obtain the stresses in a body when subjected to suddenly applied or impact
loads, and
• determine the deflections of beams through the application of Castigliano’s
theorem.
Extension, Δ
Figure 8.1
The strain energy per unit volumes is after referred to as resilience. The resilience at the
yield point or proof stress is termed as proof resilience. The strain energy stored upto the
fracture of the material is generally known as toughness.
In the earlier units, it is seen that loads may be applied on a body in different modes, viz.,
as axial loads, as shear loads, as bending loads or as torsional loads.
In the following sections, we shall derive expressions for the strain energy stored under
each of the above forms of loading.
8.3.1 Strain Energy in Axial Loading
Consider the bar of length L and cross-sectional area A shown in Figure 8.2 to which an
axial load P is applied. Under this load, the bar may extend or contract (extend if P is
130
tension and contract if P is compression) by an amount Δ where the P – Δ relationship is Strain Energy
a shown in Figure 8.1.
Cross-sectional Area A
P P
L Δ
Figure 8.2
The work done by the force P is equal to the strain energy stored in the bar.
1
Thus, U= PΔ
2
But, we know that the deformation, Δ = strain × L
Stress P L
Δ= ×L= ×
E A E
1 PL P 2 L
∴ U = P× =
2 AE 2 AE
⎛P⎞
If we rewrite the above expression in terms of the stress ⎜ ⎟ and the strain in the bar,
⎝ A⎠
1 P P 1
U = × × × AL = × Stress × Strain × Volume
2 A AE 2
1
The strain energy per unit volume of the bar = × Stress × Strain
2
1 σ2 E ε2
= ×σ×ε= =
2 2E 2
Hence, strain energy per unit volume of the material is as follows :
ds
γ
Q
Figure 8.3
131
Forces and Stresses 1 1
in Beams Strain energy, dU, of the element = work done = Q δ = Q γ ds
2 2
Shear Stress, τ Q
But, shear strain, γ = =
Shear Modulus, G GA
1 Q Q 2 ds
∴ dU = Q× ds =
2 GA 2GA
∴ Total strain energy in the bar, resulting from shear is as follows :
L
Q 2 ds Q 2 L
U = ∫ =
2 AG 2 AG
0
π2
= × Volume of the bar
2G
8.3.3 Strain Energy in Bending
Let the element now be subjected to a constant bending moment M causing it to bend into
an arc of radius R and subtending an angle dθ at the centre, as in Figure 8.4. The beam
will also have moved through an angle dθ.
Strain energy of the element, dU = work done
1
= × moment × angle of rotation
2
M
M dθ
dx
dx
R
dθ
Figure 8.4
1
= M ×dθ
2
1 ds
= M× (since, ds = R d θ)
2 R
1 M ⎛ M E⎞
= M× × ds ⎜ since, = ⎟
2 EI ⎝ I R⎠
M 2 ds
=
2 EI
L
M 2 ds
Total strain energy in bending, U = ∫
0
2 EI
132
L Strain Energy
M 2 dx
For small curvatures this reduces to U = ∫
0
2 EI
M2 L
U =
2 EI
ds
dθ
T
Figure 8.5
1
Strain energy = work done = T dθ
2
T dθ T ds
But, from the theory of simple torsion, =G× or d θ =
J ds GJ
Where J is the polar moment of inertia and G is the shear modulus of the material.
L
T 2 ds
Therefore, total strain energy resulting from torsion, U = ∫ 2 GJ
0
T2 L
For constant T, we get U = .
2 GJ
Example 8.1
Determine the diameter of an aluninium shaft which is designed to store the same
amount of strain energy per unit volume as a 50 mm diameter steel shaft of the
same length. Both shafts are subjected to equal compressive axial loads. What will
be the ratio of the stresses set up in the two shafts?
Take Esteel = Es = 201 GN/m2 and Ealuminium = Ea = 67 GN/m2.
Solution
Since the strain energy per unit volume of the two shafts is equal, we get
σ2a σ2 σ2 E 67 1
= s or a2 = a = =
2 Ea 2 Es σs Es 201 3
σs
∴ σ s = 3 σa or = 3
σa
load
But we know stress, σ = ,
area of cross-section
133
Forces and Stresses ⎡ ⎧ ⎫⎤
in Beams ⎡⎛ ⎞⎤ ⎢ ⎪ ⎪⎥
⎢⎜ P ⎟ ⎥ ⎢ ⎪ P ⎪⎥
Thus, ⎢⎜ ⎟⎥ = ⎢ 3 ⎨ ⎬⎥
⎢⎜ π Ds ⎟ ⎥ ⎢ ⎪ ⎛ π Da ⎞ ⎪ ⎥
2 2
⎢⎜ ⎟⎥ ⎪ ⎜⎜ 4 ⎟⎟ ⎪ ⎥
⎣⎝ 4 ⎠ ⎦ ⎢ ⎩⎝ ⎠ ⎭⎦
⎣
Da4
i.e. =3
Ds4
or ∴ Da4 = 3 × (50) 4
1
Da = (3 × (50)4 ) 4 = 65.8 mm
Since the maximum shear stress is the same for both the shafts, we have,
τ a = τb
Ta Da Tb Db ⎛ T × r T × D⎞
= ⎜ since τ = = ⎟
2Ja 2J B ⎝ J 2J ⎠
Ta Db J a
Therefore, = ×
Tb Da J b
But it is also given that the strain energy of B = 1.25 × strain energy of A.
2
Tb2 L 5 Ta2 L ⎛T ⎞ 4 J
Thus, = × or ⎜ a ⎟ = × a
2GJ b 4 2GJ B ⎝ Tb ⎠ 5 Jb
4
or Db2 = × ( Db4 + db4 ) because Db2 − db2 = Da2 = 104
5
134
4 2 Strain Energy
or ∴ Db2 = ( Db + Db2 − 102 )
5
Thus, the internal and external diameters of the hollow tube are 57.74 mm and
115.47 mm, respectively.
SAQ 1
(a) Calculate the strain energy stored in a bar of circular cross-section of
diameter and length of 0.2 m and 2 m, respectively, when it is
(i) subjected to a tensile load of 25 kN,
(ii) subjected to a torque of 25 kN m, and
(iii) subjected to a uniform bending moment of 25 kN m.
Take E = 208 GN/m2, G = 80 GN/m2 for the bar material.
(b) Compare the strain energies of two bars of the same material and length and
carrying the same gradually applied compressive load if one is 25 mm
diameter throughout and the other is turned to 20 mm diameter over half its
length, the remainder being 25 mm diameter.
If both bars are subjected to pure torsion only, compare the torsional strain
energies stored if the shear stress in both bars is limited to 75 N/mm2.
(c) A 50 mm square alloy-steel bar 1 m long is a part of a machine and must
resist an axial load of 100 N m. What must be the proportional limit of the
steel to safely resist the energy load elastically with a factor of safety of 4? E
= 200 GPa.
P Work Done, W = ½ P δ
1 (Gradually Applied Load)
P
Extension
δ
Figure 8.6
135
Forces and Stresses Since the bar is strained to the same extent in both the cases, the strain energy should be
in Beams the same for the two cases. However, strain energy is equal to the work done by the load.
∴ W = W′
1 P
or × P δ = P′ δ or P′ =
2 2
This means that the suddenly applied load which is required to produce a certain value of
instantaneous strain is half the equivalent value of static load (gradually applied load)
required to produce the same strain. Vice-versa, one can say that the effect of a suddenly
applied load is two times that achieved when the load is applied gradually. Thus, greater
care must be exercised while designing structural elements that are susceptible to sudden
loading conditions.
Load, W
Rigid Collar δ
Figure 8.7
The rod is suspended from a rigid support. When the load is dropped it will produce a
maximum instantaneous extension δ in the bar, at the point of impact. Then, the load
would have done work equal to W (h + δ). This then must be equal to the strain energy
stored in the bar. If σ is the instantaneous stress produced in the bar.
σ2
The strain energy, U = × AL
2E
σ2
Hence, W ( h + δ) = AL
2E
σL
but δ=
2E
136
Strain Energy
⎛ σL ⎞ σ 2
Therefore, W ⎜h + ⎟= × AL
⎝ E ⎠ 2E
2W 2WEh
or σ2 − ×σ− =0
A L
2W
Denoting as σ st , which would have been the stress in the bar if W were to have been
A
applied gradually, the above equation can be written as
2 Eh
σ2 − 2σ σ st − σ st = 0
L
The solution of this equation results as follows :
⎡ 2 Eh ⎤ W ⎡ 2 EhA ⎤
σ = σ st ⎢1 + 1 + ⎥= ⎢1 + 1 + ⎥
⎣⎢ L × σ st ⎥⎦ A ⎣ LW ⎦
SAQ 2
An axial pull of 50 kN is suddenly applied to a steel rod 2 m long and 1000 mm2 in
cross-section. Calculate the strain energy that can be absorbed, if E = 200 GPa.
Figure 8.8
137
Forces and Stresses Then, the strain energy, U = Work done by W = W (h + δ)
in Beams
where δ is the instantaneous maximum deflection under the load W.
If We is an equivalent static load which, when applied gradually, produces the same
deflection as the impact load W, the work done by We will be equal to the work done
by W.
1
i.e. We δ = W ( h + δ)
2
Let δst be the deflection of the beam if W were applied gradually. Then, by the linear
relationship of load to deflection, it can be written
W We δ
= or We = W ×
δ st δ δ st
On simplifying, we get
δ2 − 2δ st δ − 2δ st h = 0
2h
δ2 = δ st ± δ st 1+
δ st
⎛ 2h ⎞
δ = δ st ⎜ 1 ± 1 + ⎟
⎜ δ st ⎟⎠
⎝
This then will be the instantaneous deflection under the impact load W.
Once δ is known one can obtain We and the maximum bending moment, Mmax produced
in the beam by the application of We. Then, the maximum stress σmax occurring in the
beam can be obtained from the bending equation as
M max C
σmax = ,
I
where I is the moment of inertia for the cross-section and C is the extreme fibre distance
from the neutral axis.
SAQ 3
(a) A bar 3 m long and 5 cm diameter hangs vertically and has a collar securely
attached at the lower end. Find the maximum stress induced when,
(i) a weight 2.5 kN falls from a height of 12 cm on the collar, and
(ii) a weight of 25 kN falls from a height of 1 cm on the collar.
(b) An unknown weight falls through 10 mm onto a collar rigidly attached to
the lower end of a vertical bar 3 m long and 600 mm2 in section. If the
maximum instantaneous extension is known to be 2 mm, what is the
corresponding stress and the value of the unknown weight.
Take E = 200 kN/mm2.
138
Example 8.3 Strain Energy
W ⎛ 2 EhA ⎞
We have σ= ⎜⎜ 1 + 1 + ⎟
A ⎝ LW ⎟⎠
This shows that the error in the stress calculation is insignificant even if the
elongation is neglected.
Example 8.4
When a concentrated force of 1 kN is applied at the midspan point of a simply
supported beam a static deflection of 5 mm is produced. The same load produces a
maximum stress of 158 MN/m2. Determine the magnitude of the instantaneous
stress produced when a weight of 10 kg is allowed to fall through a height of
12 mm onto the beam at midspan. What will be the instantaneous deflection?
Solution
Let We be the equivalent gradually applied load which produces the same static
deflection δmax as that of the impact load.
We 1 kN
Then, = = 200 kN/m
δmax 0.005 m
Now, equating the work done by the impact load to that done by the equivalent
static load, we get
We
δmax = W (h + δmax )
2
We We 10 × 9.81 ⎛ W ⎞
i.e. × = ⎜ 0.012 + e ⎟
2 200 1000 ⎝ 200 ⎠
where δm is the displacement of the point of application of the load Pi in the direction
of Pm.
140
Let us prove the theorem now. Consider the beam in Figure 8.9 with forces PA, PB, PC, B
Strain Energy
etc. acting at points A, B, C, etc.
PA PB PC
(+ Δ PA)
A B C
Unloaded Beam Position Beam
loaded with PA, PB, PC
δA δB δC
ΔδA ΔδB
Figure 8.9
Let δA, δB, δC, etc. be the deflections in the direction of loads at the points A, B, C, . . .,
B
respectively.
Then, the strain energy of the system is equal to the work done.
1 1 1
U = PA δ A + PB δ B + PC δC + . . .
2 2 2
If one of the loads PA is now increased by an amount Δ PA, the changes in the deflection
will be Δ δA, Δ δB, Δ δC, etc. as shown in Figure 8.9. In increasing the load from PA to
B
(PA + Δ PA), extra work is done and it is as follows (refer Figure 8.10 also).
⎛ 1 ⎞
Extra work done at A = ⎜ PA + Δ PA ⎟ × Δ δ A
⎝ 2 ⎠
Extra work done at B = PB × Δ δ B
Extra work done at C = PC × Δ δC
Similarly, in the case of other load also. Then, the total extra work one,
1
= PA Δ δ A + Δ PB Δ δ B + PB Δ δ B + PC Δ δC + . . .
2
This should be equal to the increase in strain energy Δ U.
Neglecting the product of small quantities (Δ PA Δ δA) since they are very small as
compared to the other terms.
Δ U = PA Δ δ A + PB Δ δ B + PC Δ δC + . . .
Load At A Load At B
PA + Δ PA PA
Extra Work Extra Work
PA = (PA + ½ Δ PA) ΔδA = PB +Δ δB
Deflection at A Deflection at B
δA δA + Δ δ A δB δB + ΔδB
Figure 8.10
141
Forces and Stresses Let us now consider the case where we apply the loads (PA + Δ PA), PB, PC, etc. gradually
B
in Beams to the full value. As the final loads in the case are exactly the same as the earlier case and
as all the loads are applied statically, the final deflections should also be the same.
Then, the strain energy in the case, U′ (neglecting product of small terms Δ PA Δ δA)
would be as follows :
1 1 1
U′ = ( PA + Δ PA ) (δ A + Δ δ A ) + PB (δ B + Δ δ B ) + PC (δC + Δ δC ) + . . .
2 2 2
1 1 1 1 1
= PA δ A + PA Δ δ A + Δ PA δ A + PB δ B + PB Δ δ B
2 2 2 2 2
1 1
+ PC δC + PC Δ δC + . . .
2 2
1 1 1 1 1
= PA δ A + PB δ B + PC δC + . . . + PA Δ δ A + PB Δ δ B
2 2 2 2 2
1 1
+ PC Δ δC + . . . + Δ PA δ A
2 2
1 1 1 1
=U + PA Δ δ A + PB Δ δ B + PC Δ δC + . . . + PA δ A
2 2 2 2
However, since loads as well as deflections are the same in both the cases, the total strain
energy should also be the same. Hence,
U ′ = U + ΔU
1 1 1 1
or ΔU = PA Δ δ A + PB Δ δ B + PC Δ δC + . . . + Δ PA δ A
2 2 2 2
∴ 2Δ U = PA Δ δ A + PB Δ δ B + PC Δ δC + . . . + Δ PA δ A
Earlier we have obtained
Δ U = PA Δ δ A + PB Δ δ B + PC Δ δC + . . .
Now, from the last two expressions, we get
ΔU
Δ U = Δ PA δ A or = δA
Δ PA
∂U
In the limit, δ A =
∂PA
∂U ∂U
Similarly, it can be shown that = δB , = δC , etc.
∂PB ∂PC
∂U
Thus, in general, = δM .
∂PM
As a simple example of the application of Castigliano’s theorem, let us consider a bar of
cross-sectional area A and length L under axial tension P. From earlier section we have
known that the strain energy U in this case is given as
P2 L
U =
2 AE
If we differentiate U with P, we get,
∂U PL
=
∂P AE
which we recognise as the extension of the bar under the load P, i.e. the displacement of
the point of application of the load P.
142
8.6.1 Application of Castigliano’s Theorem to Angular Movements Strain Energy
The terms load and displacement in Castigliano’s theorem have a more general meaning
than is indicated in the previous section. The load may be a couple and the displacement
a rotation produced by the couple. Thus, as per Castigliano’s theorem
“If the total strain energy, expressed in terms of the external moments, be partially
differentiated with respect to one of the moments, the result is the angular deflection (in
radians) of the point of application of that moment and in its direction.”
This form of the statement will be useful in determining rotations in either bending or
torsion problems.
8.6.2 Application to Truss Problems
In this section, we shall apply Castigliano’s theorem to determine the deflections in a
truss.
Let a truss with m members be subjected to a system of loads P1, P2, . . . , etc. The forces
in the members of the truss can be obtained by standard methods. Let Fi be the force in
the ith member due to the above system of forces. Let Li and Ai be respectively the length
and area of cross-section of the member. Then the strain energy stored in this member
will be
Fi 2 Li
2 Ai E
where E is the Young’s modulus of the material.
Then, the total strain energy stored in the truss can be written as follows :
m Fi 2 Li
U = ∑
i =1 2 Ai E
Now, as per Castigliano’s theorem the deflection δn under the load Pn in its direction will
be
∂U ∂ m Fi 2 Li m ⎛ F L ⎞ ⎛ ∂F ⎞
δn = = ∑ = ∑ ⎜ i i⎟⎜ i ⎟
∂Pn ∂Pn i = 1 2 Ai E i = 1 ⎝ Ai E ⎠ ⎝ ∂Pn ⎠
Let us observe this equation closely.
∂Fi
is the rate at which Fi varies with Pn with all other loads remaining constant or in
∂Pn
other words it is the force in the member when a unit load is applied instead of Pn, while
all other loads are zero.
The above concept can still be applied to find the deflection in any direction at any joint
of the truss even when no load is acting in that direction or at that point. In this case, we
introduce a fictitious force Q at the pint in the direction in which we need the deflection.
Then we determine the member forces for the load system which now includes Q also.
m Fi 2 Li
Then, U = ∑ (Fi now includes the effect of Q also)
i =1 2 Ai E
∂U m F L ∂F
Therefore, δQ = = ∑ i i i
∂Q i = 1 Ai E ∂Q
However, in actual problem there is no Q or its value is zero. So we put Q = 0 in the
above expression and obtain δQ.
The concept of unit load can be brought in even here. Now when we include Q in the
load system Fi can be written as
Fi = FiP + FiQ (principle of superposition)
143
Forces and Stresses where FiP and FiQ are the member forces only due to the Ps and Q, respectively.
in Beams
( FiP + FiQ )2 Li
Then, U =∑
2 Ai Ei
m ( FiP + FiQ ) Li ∂
and δQ = ∑ × ( FiP + FiQ )
i =1 Ai E ∂Q
∂FiP
(since FiP does not involve Q and, hence, = 0 ).
∂Q
When we put Q = 0 in the above expression it reduces to
m ⎛F L ⎞ ⎛ ∂FiQ ⎞
δQ = ∑ ⎜ iP i ⎟⎜ ⎟
i = 1 ⎝ Ai E ⎠ ⎝ ∂Q ⎠
∂FiQ
is nothing but the member force due to a unit load in the direction of Q at the point
∂Q
of application of Q while all other loads are zero.
Example 8.5
In the truss shown in Figure 8.11, the cross-sectional area of its top and bottom
chords is equal to 2.4 × 103 mm2, while the cross-sectional area of its diagonal and
vertical members is equal to 1.6 × 103 mm2. The modulus of elasticity of the
material form which the members are made is E = 200 GPa. Compute the vertical
and horizontal components of the deflection at joint 7.
Solution
Vertical Deflection at Joint 7
Since there is a 40 kN acting vertically at joint 7, this is a case of straight
forward application of Castigliano’s theorem.
As a first step, the forces in the members of the truss (F) are determined
under the action of the two 40 KN loads that are acting on the truss. These
FL
forces are shown in column 2 of Table 8.1. Then, the respective are
AE
calculated. These are presented in column 3.
2 9 4 10 6
6 8
4
5 7 3.2 k N
1 2 3
1 3 5 7
40 k N 40 k N
3 @ 2.4 m = 7.2 m
Figure 8.11
144
Table 8.1 Strain Energy
Membe F* A L FL ∂F FL ∂F
r (kN) (102 (mm) .
AE ∂P AE ∂P
mm2)
1 − 60 2.4 2400 − 0.30 − 1.5 0.450
2 − 60 2.4 2400 − 0.30 − 1.5 0.450
3 − 30 2.4 2400 − 0.15 − 0.75 0.1125
4 − 100 1.6 4000 − 1.25 − 1.25 1.5625
5 40 1.6 3200 0.40 0 0.0000
6 50 1.6 4000 0.625 1.25 0.78125
7 − 40 1.6 3200 − 0.40 − 1.00 0.400
8 50 1.6 4000 0.625 1.25 0.78125
9 120 2.4 2400 0.60 2.25 1.350
10 30 2.4 2400 0.15 0.75 0.1125
⎛ FL ⎞ ⎛ ∂F ⎞ + 6.000**
δ7 =∑⎜ ⎟⎜ ⎟ mm
(vertical) ⎝ AE ⎠ ⎝ ∂P ⎠
Next, all the external loads are replaced by a single unit load acting
vertically down at joint 7 and the member forces are obtained. The result
∂F
will be equal to . These values are presented in column 6 of the
∂P
Table 8.1. Then, the rest of the calculations are completed to obtain
FL ∂F
∑ . which according to Castigliano’s theorem will be equal to the
AE ∂P
vertical deflection at joint 7. We can see this deflection is 6 mm downwards.
Horizontal Deflection at Joint 7
In this case, there is no horizontal load at joint 7. So to apply Castigliano’s
theorem, we have to first introduce a fictitious horizontal load Q at joint 7,
then find the members forces and differentiate these forces with respect to
FL ∂F
Q. Finally, ∑ . is determined. Then, Q unit is made equal to zero in
AE ∂Q
this expression to get δ7 (horizontal). However, here we use the unit load
concept and solve this problem. But you are advised to check the result by
the above procedure.
Here also, the first step is to determine the member forces (F) under the
external load system. Thus, column 2 of the table remains the same. Now
the external load system is replaced by a single horizontal unit load (because
we require horizontal deflection). Let us assume that this unit load acts left
to right. Now, the member forces are determined again under the action of
⎛ ∂F ⎞
this unit load alone. The result will be equivalent to ⎜ ⎟
⎝ ∂Q ⎠Q = 1
at 7. These values are shown in column 6 of Table 8.2.
⎛ FL ∂F ⎞
Now the horizontal deflection at joint 7 is given by ⎜ . ⎟ . The results
⎝ AE ∂Q ⎠
are presented in Table 8.2. From this table, it can be seen that horizontal
deflection of joint 7 is 0.75 mm towards the left.
Table 8.2 145
Forces and Stresses Membe F* A L FL ∂F FL ∂F
in Beams .
r (kN) (102 (mm)
AE ∂Q AE ∂Q
mm2)
1 − 60 2.4 2400 − 0.30 1 − 0.30
2 − 60 2.4 2400 − 0.30 1 − 0.30
3 − 30 2.4 2400 − 0.15 1 − 0.15
4 − 100 1.6 4000 − 1.25 0 0.0
5 40 1.6 3200 0.40 0 0.0
6 50 1.6 4000 0.625 0 0.0
7 − 40 1.6 3200 − 0.40 0 0.0
8 50 1.6 4000 0.625 0 0.0
9 120 2.4 2400 0.60 0 0.0
10 30 2.4 2400 0.15 0 0.0
⎛ FL ⎞ ⎛ ∂F ⎞ − 0.75* mm
δ7 =∑⎜ ⎟⎜ ⎟
(horizontal) ⎝ AE ⎠ ⎝ ∂Q ⎠
* Negative sign indicates that the deflection is opposite to the direction of unit load.
SAQ 5
Determine the vertical and horizontal deflection of point A in the pin-joined frame
shown in Figure 8.12. All members have EA = 20 MN.
Given is AB = AC = BC = 1 m and DB = 0.5 m.
B
A
C
50 kN
Figure 8.12
The bending moment will be a function of the loads. Then, by Castigliano’s theorem the
deflection under one of the loads, P, in its direction will be
L L
∂U ∂ M 2 dx M ∂M
δ= =
∂P ∂P ∫ 2 EI
= ∫ .
EI ∂P
dx
0 0
146
If we have to find the deflection at a point where there is no load, a fictitious load Q has Strain Energy
to be introduced at that point in the direction in which the deflection is needed and U is
∂U
computed. Then, is computed. If in this expression Q is made equal to zero, the
∂Q
result will be the deflection at the pint of application of Q.
L
M 2 dx
U = ∫ 2 EI
0
⎡ ∂U ⎤ ⎡ ∂ L
M 2 dx ⎤ ⎡ L M ∂M ⎤
δ=⎢ ⎥ =⎢
⎣ ∂Q ⎦ Q = 0 ⎢⎣ ∂Q
∫ 2 EI ⎥⎦
⎥ = ⎢∫ .
⎢⎣ 0 EI ∂Q
dx ⎥
⎥⎦ Q = 0
0 Q=0
∂M
From the similarities with the truss problems, it can easily be recognised that will
∂Q
be the moment at a section due to a unit load at the point of application of Q with no
other loads present on the beam. We shall now apply these concepts to some practical
problems.
Example 8.6
Determine the vertical deflection at the tip of the cantilever shown in Figure 8.13.
P
E, I
Figure 8.13
Solution
Measuring the distance x from the free end, the bending moment at any section
distant x from the tip is given as
M = − Px
L L
M 2 ds P2 x2 P 2 L3
Then, U = ∫ 2 EI
= ∫ 2 EI
dx =
6 EI
0 0
∂U PL3
P= =
∂P 3EI
The positive sign shows that the deflection is in the direction of P, i.e. downwards.
∂U
It may be seen that can be written as follows :
∂P
L
⎛ − Px ⎞
∫ ⎜
⎝ EI ⎠
⎟ (− x) dx
0
Px M ∂M
where − = and − x =
EI EI ∂P
147
Forces and Stresses Example 8.7
in Beams
Determine the deflection (vertical) at a distance ‘a’ from the tip in Example 8.6.
Solution
Since in this case there is no load at the point where the deflection is needed, we
introduce a fictitious load Q as shown in Figure 8.14.
Q P
Figure 8.14
∂M
Then, M = − Px and = 0 for 0 ≤ x ≤ a
∂Q
∂M
and M = − Px − Q ( x − a) and = − ( x − a) for a ≤ x ≤ L
∂Q
L a L
M ∂M ⎛ Px ⎞ ⎛ [− Q ( x − a ) − Px] [− ( x − a)] ⎞
∴ ∫ EI ∂Q
dx = ∫ ⎜− ⎟ (0) dx +
⎝ EI ⎠
∫ ⎜
⎝ EI
⎟ dx
⎠
0 0 0
L
⎛ Px ( x − a) + Q ( x − a) 2 ⎞
= ∫ ⎜⎜
EI
⎟⎟ dx
0 ⎝ ⎠
⎡ L M ∂M ⎤ L
Px ( x − a)
δa = ⎢ ∫ . dx ⎥ = ∫ dx
⎢⎣ 0 EI ∂Q ⎥⎦ Q = 0 0
EI
⎡⎛ 2 ⎞L ⎤
x 2 ⎥ P ⎡ L3 a3 L2 a ⎤
L
P ⎢⎜ ( x − a ) x ⎟ −
=
EI ⎢⎜⎝ 2 ⎟⎠ ∫ 2
dx = ⎢ +
⎥ EI ⎣⎢ 3 6
− ⎥
2 ⎦⎥
⎣ a 0 ⎦
Example 8.8
Determine the slope of the cantilever given in Example 8.6 at its free end.
Solution
To determine rotation in a beam we should have a couple acting at that point.
Since at the free end, there is no couple, we introduce a fictitious couple M′ as
shown in Figure 8.15.
P
M′
A
θA
L
Figure 8.15
148
∂M Strain Energy
Now M = − Px − M ′ and = −1
∂M ′
⎡ L M ∂M ⎤ ⎡ L ( Px + M ′) × 1 ⎤
∴ θA = ⎢∫ . ⎥ = ⎢∫ dx ⎥
⎢⎣ 0 EI ∂M ′ ⎥⎦ M ′ = 0 ⎢⎣ 0 EI ⎥⎦ M ′ = 0
PL2
θA = (Positive sign indicates rotation in the sense of M′, i.e. clockwise).
2 EI
Thus, we see that the tip slope of a cantilever loaded with a concentrated vertical
PL2
force P at its free end is .
2 EI
SAQ 6
(a) Using Castigliano’s theorem, determine the deflection at the centre of the
beam loaded as shown in Figure 8.16. The EI is constant.
L/2 L/2
Figure 8.16
(b) A steel rod, of flexural rigidity EI, forms a cantilever ABC lying in a vertical
plane as shown in Figure 8.17. A horizontal load P acts at C.
Calculate
(i) the horizontal deflection of C,
(ii) the vertical deflection of C, and
(iii) the slope at B.
Consider the strain energy resulting from bending only.
C P
Figure 8.17
149
Forces and Stresses
in Beams 8.7 SUMMARY
In this unit, we have
(a) introduced the concept of potential energy,
(b) introduced the concept of strain energy,
(c) shown the relationship between flexural work and strain energy,
(d) derived the expressions for strain energy in different modes of loading,
(e) seen the differences between gradually applied loads, suddenly applied
loads and impact loads,
(f) derived the expressions for instantaneous stresses and deflection for axial
and bending impact loads,
(g) stated and proved Castigliano’s theorem, and
(h) applied Castigliano’s theorems to determine the deflections in trusses and
beams.
150
Strain Energy
FURTHER READING
Prasad, Jainti, Strength of Materials, CBS Publishers and Distributors.
Timoshenko, Stephan, Strength of Materials – Part I and II, CBS Publishing and
Distributors.
Schaum’s Outline Series, (1989), Strength of Materials, Second Edition, McGraw-Hill
Book Company.
Popov, E. P., (1993), Mechanics of Materials, Second Edition, Prentice-Hall of India
Private Limited.
Ryder, G. H., (1993), Strength of Materials, Educational Low Priced Books Scheme.
Case, John, Chilver, L., Ross, Carl, T. F., (1993), Strength of Materials and Structures,
Third Edition, Educational Low Priced Books Scheme.
151
Forces and Stresses
in Beams FORCES AND STRESSES IN BEAMS
Computation of appropriate kinds of stresses in a structural member is a pre-requisite for
designing the structural member. These stresses are developed in the member due to
various forces acting on it. In this block, an attempt has been made to familiarise you
with various types of forces and stresses developed in the structural members.
This block consists of five units.
Unit 4 deals with the most important parameters used in designing, namely shear force
and bending moment. Unit starts with the classification of supports, beams and loadings.
It covers all types of possible combinations of loads and beams for shear force and
bending moment computation through a large number of examples.
Unit 5 covers the details of bending stresses developed across the cross-section of a
beam. It also discusses the theory of simple or pure bending and concept of moment of
resistance. Unit concludes with some practical applications of these concepts.
Unit 6 deals with stress distribution in the beams. In the beginning, unit covers the
concept of section modulus and gives the bending stress distribution diagrams for all
important geometric shapes of beam sections. It also covers the flitched beam and shear
stress distribution pattern of different cross-sections of the beam.
Unit 7 discusses the effect of combination of direct and bending stresses on the members
and gives the distribution pattern of these stresses across the cross-section. It also evolves
the condition for no tension in a section.
Unit 8 presents some concepts of strain energy useful in axial loading, shear, bending and
torsion. It also covers the Castigliaon’s theorem applied to angular movement, truss and
beam problems.
152