Beam Displacements

David Roylance

Department of Materials Science and Engineering

Massachusetts Institute of Technology

Cambridge, MA 02139

November 30, 2000

Introduction
We want to be able to predict the deection of beams in bending, because many applications
have limitations on the amount of deection that can be tolerated. Another common need for
deection analysis arises from materials testing, in which the transverse deection induced by a
bending load is measured. If we know the relation expected between the load and the deection,
we can back out the material properties (specically the modulus) from the measurement. We
will show, for instance, that the deection at the midpoint of a beam subjected to three-point
bending (beam loaded at its center and simply supported at its edges) is
3
48
where the length and the moment of inertia are geometrical parameters. If the ratio of
to is measured experimentally, the modulus can be determined. A stiness measured this
way is called the exural modulus.
There are a number of approaches to the beam deection problem, and many texts spend
a good deal of print on this subject. The following treatment outlines only a few of the more
straightforward methods, more with a goal of understanding the general concepts than with
developing a lot of facility for doing them manually. In practice, design engineers will usually
consult handbook tabulations of deection formulas as needed, so even before the computer age
many of these methods were a bit academic.
=

Multiple integration
In Module 12, we saw how two integrations of the loading function () produces rst the shear
function () and then the moment function ():
=
=

() + 1

(1)

() + 2

(2)

where the constants of integration 1 and 2 are evaluated from suitable boundary conditions on
and . (If singularity functions are used, the boundary conditions are included explicitly and
the integration constants 1 and 2 are identically zero.) From Eqn. 6 in Module 13, the curvature
1

,xx () is just the moment divided by the section modulus . Another two integrations then
give
1

,x () =

() + 3

(3)

() =

,x () + 4

(4)

where 3 and 4 are determined from boundary conditions on slope or deection.

Example 1

Figure 1: Three-point bending.

As an illustration of this process, consider the case of three-point bending shown in Fig. 1. This
geometry is often used in materials testing, as it avoids the need to clamp the specimen to the testing
apparatus. If the load is applied at the midpoint, the reaction forces at and are equal to half the
applied load. The loading function is then
() =

hi1 h i1
2
2

Integrating according to the above scheme:

() =
() =

hi0 + h i0
2
2

hi1 h i1
2
2

(5)

hi2 h i2 + 3
4
2
2
From symmetry, the beam has zero slope at the midpoint. Hence ,x = 0 @ = L/2, so 3 can be found
to be 2 /16. Integrating again:
,x () =

2
hi3 h i3
+ 4
12
6
2
16
The deection is zero at the left end, so 4 = 0. Rearranging, the beam deection is given by
() =




3
3
2
=
4 3 8h i
48
2

(6)

The maximum deection occurs at = L/2, which we can evaluate just before the singularity term
activates. Then
max =

3
48

(7)

This expression is much used in exural testing, and is the example used to begin this module.

Before the loading function () can be written, the reaction forces at the beam supports
must be determined. If the beam is statically determinate, as in the above example, this can
be done by invoking the equations of static equilibrium. Static determinacy means only two
reaction forces or moments can be present, since we have only a force balance in the direction
transverse to the beam axis and one moment equation available. A simply supported beam (one
resting on only two supports) or a simply cantilevered beam are examples of such determinate
beams; in the former case there is one reaction force at each support, and in the latter case there
is one transverse force and one moment at the clamped end.
Of course, there is no stringent engineering reason to limit the number of beam supports
to those sucient for static equilibrium. Adding extra supports will limit deformations and
stresses, and this will often be worthwhile in spite of the extra construction expense. But the
analysis is now a bit more complicated, since not all of the unknown reactions can be found from
the equations of static equilibrium. In these statically indeterminate cases it will be necessary
to invoke geometrical constraints to develop enough equations to solve the problem.
This is done by writing the slope and deection equations, carrying the unknown reaction
forces and moments as undetermined parameters. The slopes and deections are then set to
their known values at the supports, and the resulting equations solved for the unknowns. If
for instance a beam is resting on three supports, there will be three unknown reaction forces,
and we will need a total of ve equations: three for the unknown forces and two more for the
constants of integration that arise when the slope and deection equations are written. Two
of these equations are given by static equilibrium, and three more are obtained by setting the
deections at the supports to zero. The following example illustrates the procedure, which is
straightforward although tedious if done manually.
Example 2
Consider a triply-supported beam of length = 15 as shown in Fig. 2, carrying a constant uniform load
of = 10. There are not sucient equilibrium equations to determine the reaction forces , , and
, so these are left as unknowns while multiple integration is used to develop a deection equation:
() = hi1 + h 7.5i1 + h 15i1 10hi0
Z

() =
Z
() =
Z
0 () =

() = hi0 h 7.5i0 h 15i0 + 10hi1

() = hi1 + h 7.5i1 + h 15i1

() =

10 2
hi
2

10
hi2 +
h 7.5i2 +
h 15i2 hi3 + 1
2
2
2
6

Figure 2: Uniformly loaded beam resting on three supports.

Z
() =

0 () =

10
hi3 +
h 7.5i3 +
h 15i3 hi4 + 1 + 2
6
24
6
6

These equations have 5 unknowns: , , , 1 , and 2 . These must be obtained from the two
equilibrium equations
X
= 0 = + +
X

2
2
and the three known zero displacements at the supports
= 0 =

(0) = (L/2) = () = 0
Although the process is straightforward, there is a lot of algebra to wade through. Statically indeterminate
beams tend to generate tedious mathematics, but fortunately this can be reduced greatly by modern
software. Follow how easily this example is handled by the Maple V package (some of the Maple responses
removed for brevity):

>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>

# read the library containing the Heaviside function

readlib(Heaviside);

# use the Heaviside function to define singularity functions;

#
sfn(x,a,n) is same is <x-a>^n

sfn := proc(x,a,n) (x-a)^n * Heaviside(x-a) end;

# define the deflection function:

y := (x)-> (Ra/6)*sfn(x,0,3)+(Rb/6)*sfn(x,7.5,3)+(Rc/6)*sfn(x,15,3)

-(10/24)*sfn(x,0,4)+c1*x+c2;

# Now define the five constraint equations; first vertical equilibrium:

eq1 := 0=Ra+Rb+Rc-(10*15);

# rotational equilibrium:

eq2 := 0=(10*15*7.5)-Rb*7.5-Rc*15;

# Now the three zero displacements at the supports:

eq3 := y(0)=0;

eq4 := y(7.5)=0;

eq5 := y(15)=0;

# set precision; 4 digits is enough:

> Digits:=4;

> # solve the 5 equations for the 5 unknowns:

> solve({eq1,eq2,eq3,eq4,eq5},{Ra,Rb,Rc,c1,c2});

{c2 = 0, c1 = -87.82, Rb = 93.78, Ra = 28.11, Rc = 28.11}

> # assign the known values for plotting purposes:
> c1:=-87.82;c2:=0;Ra:=28.11;Rb:=93.78;Rc:=28.11;
> # the equation of the deflection curve is:
> y(x);
3
3
4.686 x Heaviside(x) + 15.63 (x - 7.5) Heaviside(x - 7.5)

>
>
>
>

3
4
+ 4.686 (x - 15) Heaviside(x - 15) - 5/12 x Heaviside(x) - 87.82 x
# plot the deflection curve:
plot(y(x),x=0..15);
# The maximum deflection occurs at the quarter points:
y(15/4);
-164.7

The plot of the deection curve is shown in Fig. 3.

Figure 3: Deection curve () for uniformly loaded triply-supported beam (Note dierence
in horizontal and vertical scales).

Energy method
The strain energy in bending as given by Eqn. 8 of Module 13 can be used to nd deections,
and this may be more convenient than successive integration if the deection at only a single
point is desired. Castiglianos Theorem gives the deection congruent to a load as
U

=
=
P P
5

2
2

It is usually more convenient to do the dierentiation before the integration, since this lowers
the order of the expression in the integrand:
Z

where here and are assumed not to vary with .

The shear contribution to bending can be obtained similarly. Knowing the shear stress
= / (omitting the subscript on for now), the strain energy due to shear can be
written
=

2
=
2

2
2

2

2

The integral over the cross-sectional area is a purely geometrical factor, and we can write
Z

(8)

where the is a dimensionless form factor for shear dened as

= 2

(9)

Figure 4: Rectangular beam section.

Evaluating for rectangular sections for illustration (see Fig. 4), we have in that case
3
12

 

(h/2)

= +

2
2
= bh,

()
=
(3 /12)2

h/2

h/2

1
6
=
2

Hence is the same for all rectangular sections, regardless of their particular dimensions.
Similarly, it can be shown (see Prob. 3) that for solid circular sections = 10/9 and for hollow
circular sections = 2.
6

Example 3
If for instance we are seeking the deection under the load in the three-point bending example done
earlier, we can dierentiate the moment given in Eqn. 5 to obtain
M
1

= hi1 h( i1
P
2
2
Then
=

Z 

hi1 h i1
2
2



1 1
hi (
h i1
2
2

Expanding this and adjusting the limits of integration to account for singularity functions that have not
been activated:
)
(Z
"

 
2 #
Z
2

+
=

4
2
2
0
L/2
=

3
48

as before.
The contribution of shear to the deection can be found by using = P/2 in the equation for strain
energy. For the case of a rectangular beam with = 6/5 we have:
(P/2)2 (6/5)

2
6
U
=
=
P
20

=
P,

The shear contribution can be compared with the bending contribution by replacing with 12I/2 (since
= and = 3 /12). Then the ratio of the shear to bending contributions is
32
2 /40
=
3 /24
52
Hence the importance of the shear term scales as (h/L)2 , i.e. quadratically as the span-to-depth ratio.

The energy method is often convenient for systems having complicated geometries and com
bined loading. For slender shafts transmitting axial, torsional, bending and shearing loads the
strain energy is
Z

2
2
2
2
+
+
+
2 2
2
2

(10)

Example 4
Consider a cantilevered circular beam as shown in Fig. 5 that tapers from radius 1 to 2 over the length
. We wish to determine the deection caused by a force applied to the free end of the beam, at an
angle from the horizontal. Turning to Maple to avoid the algebraic tedium, the dimensional parameters
needed in Eqn. 10 are dened as:

Figure 5: Tapered circular beam.

>
>
>
>

r
A
Iz
Jp

:=
:=
:=
:=

proc
proc
proc
proc

(x)
(r)
(r)
(r)

r1 + (r2-r1)*(x/L) end;
Pi*(r(x))^2 end;
Pi*(r(x))^4 /4 end;
Pi*(r(x))^4 /2 end;

where r(x) is the radius, A(r) is the section area, Iz is the rectangular moment of inertia, and Jp is the
polar moment of inertia. The axial, bending, and shear loads are given in terms of as
> P := F* cos(theta);

> V := F* sin(theta);

> M := proc (x) -F* sin(theta) * x end;

The strain energies corresponding to tension, bending and shear are

>
>
>
>

U1
U2
U3
U

:=
:=
:=
:=

P^2/(2*E*A(r));
(M(x))^2/(2*E*Iz(r));
V^2*(10/9)/(2*G*A(r));
int( U1+U2+U3, x=0..L);

Finally, the deection congruent to the load is obtained by dierentiating the total strain energy:
> dF := diff(U,F);
The result of these manipulations yields


12 2 122 cos2 + 922 cos2 + 10 22 10 22 cos2
=
9 1 23 E
This displacement is in the direction of the applied force ; the horizontal and vertical deections of the
end of the beam are then
= cos
= sin

Superposition
In practice, many beams will be loaded in a complicated manner consisting of several concentrated or distributed loads acting at various locations along the beam. Although these multipleload cases can be solved from scratch using the methods described above, it is often easier to
solve the problem by superposing solutions of simpler problems whose solutions are tabulated.
Fig. 6 gives an abbreviated collection of deection formulas1 that will suce for many problems.
The superposition approach is valid since the governing equations are linear; hence the response
to a combination of loads is the sum of the responses that would be generated by each separate
load acting alone.

Figure 6: Deections for cantilevered and simply-supported beams, under concentrated and
distributed loading.

Example 5
We wish to nd the equation of the deection curve for a simply-supported beam loaded in symmetric
four-point bending as shown in Fig. 7. From Fig. 6, the deection of a beam with a single load at a
1

A more exhaustive listing is available in W.C. Young, Roarks Formulas for Stress and Strain, McGraw-Hill,
New York, 1989.

Figure 7: Four-point bending.





3
3
2
2
. Our present problem is just
distance from the left end is () = 6
h i +
two such loads acting simultaneously, so we have




( )

h i3 3 + 2 ( )2
() =
6




h ( )i3 3 + 2 2
+
6

In some cases the designer may not need the entire deection curve, and superposition of
tabulated results for maximum deection and slope is equally valid.

Problems
1. (a)(h) Write expressions for the slope and deection curves of the beams shown here.

Prob. 1
2. (a)(h) Use MapleV (or other) software to plot the slope and deection curves for the beams
in Prob. 1, using the values (as needed) = 25 in, a = 15 in, w = 10 lb/in, P = 150 lb.
3. Show that the shape factor for shear for a circular cross section is

10


= 2

10
=
2

4. (a)(b) Determine the deection curves for the beams shown here. Plot these curves for
the the values (as needed) = 25 in, a = 5 in, w = 10 lb/in, P = 150 lb.

Prob. 4
5. (a) Determine the deection of a coil spring under the inuence of an axial force ,
including the contribution of bending, direct shear, and torsional shear eects. Using
= 1 mm and = 10 mm, compute the relative magnitudes of the three contributions.
(b) Repeat the solution in (a), but take the axial load to be placed at the outer radius of
the coil.

Prob. 5
6. (a)(c) Use the method of superposition to write expressions for the deection curve ()
for the cases shown here.

11

Prob. 6

12

MIT OpenCourseWare
http://ocw.mit.edu

3.11 Mechanics of Materials

Fall 1999

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.