Unit 8 Strain Energy: Structure

Download as pdf or txt
Download as pdf or txt
You are on page 1of 24

Strain Energy

UNIT 8 STRAIN ENERGY


Structure
8.1 Introduction
Objectives
8.2 Potential Energy
8.3 Concept of Strain Energy
8.3.1 Strain Energy in Axial Loading
8.3.2 Strain Energy in Shear
8.3.3 Strain Energy in Bending
8.3.4 Strain Energy in Torsion

8.4 Suddenly Applied Loads


8.5 Impact Loads
8.5.1 Stresses Due to Axial Loads
8.5.2 Stresses Due to Bending Loads

8.6 Castigliano’s Theorem


8.6.1 Application of Castigliano’s Theorem to Angular Movements
8.6.2 Application to Truss Problems
8.6.3 Application to Beam Problems

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.

8.2 POTENTIAL ENERGY


Energy is defined as capacity to do work. It may exist in many forms, e.g. mechanical
(potential or kinetic), thermal, nuclear, chemical, etc. It is also well known that energy is
conserved through transformation from one form to another. The potential energy of a
body is the form of energy which is stored by virtue of the work which was previously
been done on that body, e.g. in lifting it to some height above a datum. In this case, work
is done on the body by displacing it against the force of gravity. This potential energy is
129
Forces and Stresses converted to an equivalent kinetic energy if the body is allowed to fall from the height by
in Beams developing velocity of fall. Let us now consider the case of a deformable body subjected
to certain external forces or what we call as loads. Because of the property of
deformability of the body, each of the points of application of these forces may undergo
certain displacement. Thus, the force moves through a certain distance, which is
equivalent to doing work on the body, i.e. a certain energy imparted to the body. One
may now ask a question, ‘where did all this energy go’? because it had to be conserved.
This question, we shall answer in the next section.

8.3 CONCEPT OF STRAIN ENERGY


From some of the earlier units, we know that when a deformable body is subjected to
external forces the body internally develops stresses and strains. For example, if we apply
a compressive force on a spring, the spring deforms developing internally shear stresses
and shear strains. At the same time we observe that the point of application of the
compressive force moves through a distance of the spring deflection (amount of
contraction). We also know that this compressed spring has a capacity to do work when
released. This means that a certain amount of energy is stored in the spring in its strained
state. Earlier we noted that work was done on the spring by the compressive force. Now,
when we link the above two statements we realise that the work done on the spring is
transformed into energy stored in it by virtue of the strain in the spring. This form of
energy emanating out of the internal strain in the body is known as strain energy. Thus,
strain energy is a particular form of potential energy which is stored within materials
which have been subjected to strain as a result of work done on the material.
Also, we have seen that strain energy,
U = external work done.
Thus, for a gradually applied load the work done in straining the body is given by the
shaded area under the load-extension graph of Figure 8.1.
1
U = PΔ
2
Load, P

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 :

(Stress) 2 Young's modulus


= = × (Strain) 2
2 × Young's modulus 2

8.3.2 Strain Energy in Shear


Consider now an element of length ds of the bar subjected to a shear load Q at one end
causing deformation through an angle γ (the shear strain) and a shear deflection δ, as
shown in Figure 8.3.

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

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

where dx is the undeformed element length.


If M is constant throughout the length L, then

M2 L
U =
2 EI

8.3.4 Strain Energy in Torsion


The element ds is now considered subject to a torque T as shown in Figure 8.5 producing
an angle of twist dθ.

ds

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

The required diameter of the aluminium shaft is 65.8 mm.


Example 8.2
Two shafts are of the same material, length and weight. One is solid and of
100 mm diameter, the outer is hollow. If the hollow shaft is to be store 25% more
energy than the solid shaft when transmitting torque, what must be its internal and
external diameters? Assume the same maximum shear stress in both the shafts.
Solution
Let A be the solid shaft of diameter Da and B be the hollow shaft with external
diameter Db and internal diameter db. Since both the shafts are to be the same
material and of equal weight, their volumes should be equal.

π Db2 π ( Db2 − db2 )


Therefore, ×L= ×L
4 4

or Da2 = Db2 − d b2 = (100) 2 mm 2

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

From the relationship, this reduces to


2
⎛ Db ⎞ 4 J b 4 ( Db4 − db4 )
⎜ ⎟ = × = ×
⎝ Da ⎠ 5 Ja 5 Da4

4 ( Db4 − db4 ) 4 ( Db2 + db2 ) ( Db2 − db2 )


Then, ∴ Db2 = × = ×
5 Da4 5 Da2

4
or Db2 = × ( Db4 + db4 ) because Db2 − db2 = Da2 = 104
5
134
4 2 Strain Energy
or ∴ Db2 = ( Db + Db2 − 102 )
5

3Db2 = 4 × 104 or Db = 115.47 mm

Therefore, ∴ db2 = (115.47) 2 − 104 or db = 57.74 mm

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.

8.4 SUDDENLY APPLIED LOADS


If a load P is applied gradually to a bar to produce an extension δ, the load-extension
graph will be as shown in Figure 8.1. However, if a load P′ is applied suddenly to the bar
(by suddenly we mean that the value of P′ is reached instantaneously but not gradually as
is the case with P) so as to produce the same extension δ, the load extension graph now
looks as in Figure 8.6.

Work Done (Suddenly Applied Load) W ′ = P ′ δ

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.

8.5 IMPACT LOADS


If a body that possesses kinetic energy by virtue of its motion strikes an elastic member
and is brought to rest, it can be conservatively assumed that all the kinetic energy of the
body is converted into strain energy in the member. Equally, if the body falls from a
certain height, the potential energy that the body has lost or the kinetic energy gained by
the body will be converted into strain energy in the member. Such a loading where there
is a transfer of kinetic energy into strain energy is termed as impact loading.

8.5.1 Stresses Due to Axial Loads


Consider the problem of a vertical rod of strength L and cross-sectional area A shown in
Figure 8.7. A load W, concentric with the rod and free to slide vertically, is dropped from
a height h on the end stop which may be rigid collar.

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 ⎦

The instantaneous deflection δ then can be obtained from


σL
δ=
E
If the load is not dropped but suddenly applied from effectively zero height (h = 0 in the
above equation), the stress developed would be
σ = 2σ st

It is a result which we have already predicted in the earlier section.


Thus, it can be seen that the stresses resulting from impact loads will be even higher and
any design work where impact loading is possible should always employ a safety factor
well in excess of two.

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.

8.5.2 Stresses Due to Bending Loads


Consider a simply supported beam as shown in Figure 8.8 subjected to an impact load
produced by a weight W falling from a height h. it will be assumed that the weight of the
beam can be neglected in comparison with W and that the work done by W is completely
transformed into strain energy of bending.

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

Substituting this in the above work equation, we get


W δ
× × δ = W ( h + δ)
2 δ st

On simplifying, we get

δ2 − 2δ st δ − 2δ st h = 0

Solution of this quadratic equation gives

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

What will be the instantaneous stress and elongation of a 25 mm diameter bar,


2.6 m long, suspended vertically, if a mass of 10 kg falls through a height of
300 mm onto a collar which is rigidly attached to the bottom end of the bar? Take
E = 200 GPa.
Solution

W ⎛ 2 EhA ⎞
We have σ= ⎜⎜ 1 + 1 + ⎟
A ⎝ LW ⎟⎠

10 × 9.81 × 4 ⎛ 2 × 200 × 109 × 0.3 × π × (0.025) 2 ⎞


= ⎜1 + 1 + ⎟
π (0.025) 2 ⎜⎝ 4 × 10 × 9.81 × 2.6 ⎟

= 96.24 MPa
σL 96.24 × 106 × 2.6
Now, δ= = = 0.00125 m = 1.25 mm
E 200 × 109
Thus, the maximums tress and elongation are 96.24 MPa and 1.25 mm,
respectively.
This elongation of 1.25 mm is very small compared to the height of fall namely
300 mm. Let us find out the stress if this elongation were to be neglected in the
calculation. When it is neglected the work strain energy equation becomes.
σ2
Wh = × AL
2E

2WhE 2 × 10 × 9.81 × 0.3 × 200 × 109 × 4


or σ= = = 96.04 MPa
AL π × (0.025) 2 × 2.6

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 ⎠

∴ We2 − 0.9162 We − 0.47088 = 0


139
Forces and Stresses
0.1962 + (0.1962)2 + 4 × 0.47088
in Beams
Hence, We =
2
= 0.981 + 0.6932 = 0.791 kN
We 0.791
But δmax = = = 0.003955 m = 3.955 mm ; 4 mm
200 200
It is given that the load of 1 kN produces a maximum stress of 158 MN/m2.
Therefore, a load of We should produce a maximum stress, σmax equal to
(158 We) MN/m2.

∴ σmax = 158 × 0.791 = 124.978 MN/m 2 ; 125 MN/m 2


Thus, the maximum stress and maximum deflection are 125 MN/m2 and 4 mm,
respectively.
SAQ 4
(a) When a load of 20 kN is gradually applied at a certain point on a beam, it
produces a deflection of 13 mm and a maximum bending stress of
75 MN/m2. From what height can a load of 5 kN fall onto the beam at this
point if the maximum bending stress is to be 150 MN/m2.
(b) A beam uniform cross-section with centroid at mid-depth and length 7 m, is
simply supported at its ends and carries a point load of 5 kN at 3 m from one
end. If the maximum bending stress is not to exceed 90 MN /m2 and the
beam is 150 mm deep,
(i) working from the first principles, find the deflection under the load,
and
(ii) what load dropped from a height of 75 mm onto the beam at 3 m from
one end would produce a stress of 150 MN/m2 at the point of
application of the load?
Take E = 200 GN/m2.

8.6 CASTIGLIANO’S THEOREM


In the previous sections we have seen that the strain energy stored in a elastic body is
equal to the work done by the external forces on the body during the process of
deformation. This suggests that there should be some way by which one can determine
the displacements of the body given the strain energy stored in the body. A general
theory for determining displacements of points in an elastic body was first stated by
Castigliano. Castigliano stated that
“If the total strain energy of a body or a framework is expressed in terms of the external
loads and is partially differentiate with respect to one of the loads, the result is the
displacement of the point of application of the load in the direction of the load.”
Mathematically speaking if the total strain energy, U, is a function of the loads
Pi (i = 1 to n), then
∂U
δm =
∂Pm

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

Beam loaded with PA, PB, PC ,


ΔδC etc. Plus Extra Load Δ PA

Δδ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

m ( FiP + FiQ ) Li ⎛ ∂F ∂FiQ ⎞


= ∑ × ⎜ iP + ⎟
i =1 Ai E ⎝ ∂Q ∂Q ⎠
m ( FiP + FiQ ) ∂FiQ
δQ = ∑ × Li ×
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 ⎠

* Sign convention : − ve compression and + ve tension.


** Positive sign indicates that the deflection is in the direction of the load. i.e. downwards.

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

8.6.3 Application to Beam Problems


In beams, the strain energy stored is mostly due to bending with shear contributing less
and this contribution is generally neglected. Thus, the strain energy U can be written as
L L
M 2 ds M 2 dx
U =∫ =∫ (for very small curvatures)
0
2 EI 0
2 EI

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

Hence, the vertical deflection below

∂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.

W Per Unit Length

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.

8.8 ANSWERS TO SAQs


SAQ 1
(a) 0.096 N m, 49.7 N m, 38.2 N m
(b) 0.78, 0.22
(c) 235 MPa
SAQ 2
50 kN m
SAQ 3
(a) 142.83 N/mm2, 143.6 N/mm2
(b) 133.3 N/mm2, 6666.7 N
SAQ 4
(a) 78 mm
(b) 24 mm, 1.45 kN
SAQ 5
9.2 mm, 1.45 kN
SAQ 6
5WL4
(a)
786 EI
Pa 2 Pab 2 Pab
(b) ( a + 3b), ,
3EI 2 EI EI

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

You might also like