Virtual Work 1011

Structural Analysis III
Dr. C. Caprani 1

Virtual Work
3rd Year
Structural Engineering

2010/11

Dr. Colin Caprani
Dr. C. Caprani 2
Contents
1. Introduction ......................................................................................................... 4
1.1 General .............................................................................................................. 4
1.2 Background ....................................................................................................... 5
2. The Principle of Virtual Work ......................................................................... 14
2.1 Definition ........................................................................................................ 14
2.2 Virtual Displacements..................................................................................... 17
2.3 Virtual Forces ................................................................................................. 18
2.4 Simple Proof using Virtual Displacements .................................................... 19
2.5 Internal and External Virtual Work ................................................................ 20
3. Application of Virtual Displacements ............................................................. 22
3.1 Rigid Bodies ................................................................................................... 22
3.2 Deformable Bodies ......................................................................................... 29
3.3 Problems ......................................................................................................... 40
4. Application of Virtual Forces ........................................................................... 42
4.1 Basis ................................................................................................................ 42
4.2 Deflection of Trusses ...................................................................................... 43
4.3 Deflection of Beams and Frames .................................................................... 50
4.4 Integration of Bending Moments .................................................................... 57
4.5 Problems ......................................................................................................... 60
5. Virtual Work for Indeterminate Structures ................................................... 68
5.1 General Approach ........................................................................................... 68
5.2 Using Virtual Work to Find the Multiplier ..................................................... 70
5.3 Basic Example Propped Cantilever ............................................................. 73
5.4 Indeterminate Trusses ..................................................................................... 74
5.5 Indeterminate Frames ..................................................................................... 78
5.6 Continuous Beams .......................................................................................... 86
5.7 Problems ......................................................................................................... 95
Dr. C. Caprani 3
6. Virtual Work for Self-Stressed Structures ..................................................... 98
6.1 Background ..................................................................................................... 98
6.2 Trusses .......................................................................................................... 104
6.3 Beams ............................................................................................................ 110
6.4 Frames ........................................................................................................... 113
6.5 Problems ....................................................................................................... 118
7. Past Exam Questions ....................................................................................... 120
7.1 Summer 1997 ................................................................................................ 120
7.2 Summer 1998 ................................................................................................ 121
7.3 Summer 1999 ................................................................................................ 122
7.4 Summer 2000 ................................................................................................ 123
7.5 Summer 2001 ................................................................................................ 124
7.6 Summer 2002 ................................................................................................ 125
7.7 Summer 2004 ................................................................................................ 126
7.8 Summer 2007 ................................................................................................ 127
7.9 Semester 2, 2008 ........................................................................................... 128
7.10 Semester 2, 2009 ........................................................................................ 129
7.11 Semester 2, 2009 ........................................................................................ 130
7.12 Semester 2, 2010 ........................................................................................ 131
7.13 Semester 2, 2010 ........................................................................................ 132
8. References ........................................................................................................ 133
9. Appendix Volume Integrals ......................................................................... 134

Dr. C. Caprani 4
1. Introduction
1.1 General
Virtual Work is a fundamental theory in the mechanics of bodies. So fundamental in
fact, that Newtons 3 equations of equilibrium can be derived from it. Virtual Work
provides a basis upon which vectorial mechanics (i.e. Newtons laws) can be linked
to the energy methods (i.e. Lagrangian methods) which are the basis for finite
element analysis and advanced mechanics of materials.

Virtual Work allows us to solve determinate and indeterminate structures and to
calculate their deflections. That is, it can achieve everything that all the other
methods together can achieve.

Before starting into Virtual Work there are some background concepts and theories
that need to be covered.

Dr. C. Caprani 5
1.2 Background
Strain Energy and Work Done
Strain energy is the amount of energy stored in a structural member due to
deformation caused by an external load. For example, consider this simple spring:

We can see that as it is loaded by a gradually increasing force, F, it elongates. We can
graph this as:

The line OA does not have to be straight, that is, the constitutive law of the springs
material does not have to be linear.

Dr. C. Caprani 6
An increase in the force of a small amount, F o results in a small increase in
deflection, y o . The work done during this movement (force displacement) is the
average force during the course of the movement, times the displacement undergone.
This is the same as the hatched trapezoidal area above. Thus, the increase in work
associated with this movement is:

( )
2
2
F F F
U y
F y
F y
F y
o
o o
o o
o
o
+ +
=
= +
~
(1.1)

where we can neglect second-order quantities. As 0 y o , we get:

dU F dy =

The total work done when a load is gradually applied from 0 up to a force of F is the
summation of all such small increases in work, i.e.:

0
y
U F dy =
}
(1.2)

This is the dotted area underneath the load-deflection curve of earlier and represents
the work done during the elongation of the spring. This work (or energy, as they are
the same thing) is stored in the spring and is called strain energy and denoted U.

If the load-displacement curve is that of a linearly-elastic material then F ky = where
k is the constant of proportionality (or the spring stiffness). In this case, the dotted
area under the load-deflection curve is a triangle.
Dr. C. Caprani 7

As we know that the work done is the area under this curve, then the work done by
the load in moving through the displacement the External Work Done,
e
W - is given
by:

1
2
e
W Fy = (1.3)

We can also calculate the strain energy, or Internal Work Done,
I
W , by:

0
0
2
1
2
y
y
U F dy
ky dy
ky
=
=
=
}
}

Also, since F ky = , we then have:

( )
1 1
2 2
I
W U ky y Fy = = =

Dr. C. Caprani 8
But this is the external work done,
e
W . Hence we have:

e I
W W = (1.4)

Which we may have expected from the Law of Conservation of Energy. Thus:

The external work done by external forces moving through external displacements is
equal to the strain energy stored in the material.

Dr. C. Caprani 9
Law of Conservation of Energy
For structural analysis this can be stated as:

Consider a structural system that is isolated such it neither gives nor receives
energy; the total energy of this system remains constant.

The isolation of the structure is key: we can consider a structure isolated once we
have identified and accounted for all sources of restraint and loading. For example, to
neglect the self-weight of a beam is problematic as the beam receives gravitational
energy not accounted for, possibly leading to collapse.

In the spring and force example, we have accounted for all restraints and loading (for
example we have ignored gravity by having no mass). Thus the total potential energy
of the system, H, is constant both before and after the deformation.

In structural analysis the relevant forms of energy are the potential energy of the load
and the strain energy of the material. We usually ignore heat and other energies.

Potential Energy of the Load
Since after the deformation the spring has gained strain energy, the load must have
lost potential energy, V. Hence, after deformation we have for the total potential
energy:

2
1
2
U V
ky Fy
H = +
=
(1.5)

In which the negative sign indicates a loss of potential energy for the load.
Dr. C. Caprani 10
Principle of Minimum Total Potential Energy
If we plot the total potential energy against y, equation (1.5), we get a quadratic curve
similar to:

Consider a numerical example, with the following parameters, 10 kN/m k = and
10 kN F = giving the equilibrium deflection as 1 m y F k = = . We can plot the
following quantities:

- Internal Strain Energy, or Internal Work:
2 2 2
1 1
10 5
2 2
I
U W ky y y = = = =
- Potential Energy: 10 V Fy y = =
- Total Potential Energy:
2
5 10 U V y y H= + =
- External Work:
1 1
10 5
2 2
e
W Py y y = = =

and we get the following plots (split into two for clarity):

Dr. C. Caprani 11
-10
-5
0
5
10
15
20
25
30
35
-0.5 0.0 0.5 1.0 1.5 2.0 2.5
Deflection, y
E
n
e
r
g
y
Total Potential Energy
External Work
Internal Work
Equilibrium

-40
-30
-20
-10
0
10
20
30
40
-0.5 0.0 0.5 1.0 1.5 2.0 2.5
Deflection, y
E
n
e
r
g
y
Total Potential Energy
Strain Energy
Potential Energy
Equilibrium

From these graphs we can see that because
I
W increases quadratically with y, while
the
e
W increases only linearly,
I
W always catches up with
e
W , and there will always
be a non-zero equilibrium point where
e I
W W = .

Dr. C. Caprani 12
Admittedly, these plots are mathematical: the deflection of the spring will not take up
any value; it takes that value which achieves equilibrium. At this point we consider a
small variation in the total potential energy of the system. Considering F and k to be
constant, we can only alter y. The effect of this small variation in y is:

( ) ( ) ( ) ( )
( ) ( )
( ) ( )
2
2
2
2
1 1
2 2
1 1
2
2 2
1
2
y y y k y y F y y ky Fy
k y y F y k y
ky F y k y
o o o
o o o
o o
H + H = + + +
= +
= +
(1.6)

Similarly to a first derivate, for H to be an extreme (either maximum or minimum),
the first variation must vanish:

( )
( )
1
0 ky F y o o H= = (1.7)

Therefore:

0 ky F = (1.8)

Which we recognize to be the 0
x
F =
. Thus equilibrium occurs when H is an

extreme.

Before introducing more complicating maths, an example of the above variation in
equilibrium position is the following. Think of a shopkeeper testing an old type of
scales for balance she slightly lifts one side, and if it returns to position, and no
large rotations occur, she concludes the scales is in balance. She has imposed a
Dr. C. Caprani 13
variation in displacement, and finds that since no further displacement occurs, the
structure was originally in equilibrium.

Examining the second variation (similar to a second derivate):

( )
( )
2
2
1
0
2
k y o o H = > (1.9)

We can see it is always positive and hence the extreme found was a minimum. This is
a particular proof of the general principle that structures take up deformations that
minimize the total potential energy to achieve equilibrium. In short, nature is lazy!

To summarize our findings:
- Every isolated structure has a total potential energy;
- Equilibrium occurs when structures can minimise this energy;
- A small variation of the total potential energy vanishes when the structure is in
equilibrium.

These concepts are brought together in the Principle of Virtual Work.

Dr. C. Caprani 14
2. The Principle of Virtual Work
2.1 Definition
Based upon the Principle of Minimum Total Potential Energy, we can see that any
small variation about equilibrium must do no work. Thus, the Principle of Virtual
Work states that:

A body is in equilibrium if, and only if, the virtual work of all forces acting on
the body is zero.

In this context, the word virtual means having the effect of, but not the actual form
of, what is specified. Thus we can imagine ways in which to impose virtual work,
without worrying about how it might be achieved in the physical world.

We need to remind ourselves of equilibrium between internal and external forces, as
well as compatibility of displacement, between internal and external displacements
for a very general structure:

Equilibrium Compatibility
Dr. C. Caprani 15
Virtual Work
There are two principles that arise from consideration of virtual work, and we use
either as suited to the unknown quantity (force or displacement) of the problem under
analysis.

Principle of Virtual Displacements:
Virtual work is the work done by the actual forces acting on the body moving through
a virtual displacement.

This means we solve an equilibrium problem through geometry, as shown:

(Equilibrium)
Statics
(Compatibility)
Geometry
Principle of Virtual Displacements
External Virtual Work Internal Virtual Work
Real Force
(e.g. Point load)

Virtual Displacement
(e.g. vertical deflection)
Real Force
(e.g. bending moment)

Virtual Displacement
(e.g. curvature)
equals
Either of these is the
unknown quantity of interest
External compatible
displacements
Virtual Displacements:
Geometry
Internal compatible
displacements

Dr. C. Caprani 16

Principle of Virtual Forces:
Virtual work is the work done by a virtual force acting on the body moving through
the actual displacements.

This means we solve a geometry problem through equilibrium as shown below:

(Equilibrium)
Statics
(Compatibility)
Geometry
Principle of Virtual Forces
External Virtual Work Internal Virtual Work
Virtual Force
(e.g. Point load)

Real Displacement
(e.g. vertical deflection)
Virtual Force
(e.g. bending moment)

Real Displacement
(e.g. curvature)

equals
The unknown quantity of interest
External real
forces
Internal Real Displacements:
Statics
Internal real
forces
Constitutive
relations
Internal real
displacements

Dr. C. Caprani 17
2.2 Virtual Displacements
A virtual displacement is a displacement that is only imagined to occur. This is
exactly what we did when we considered the vanishing of the first variation of H; we
found equilibrium. Thus the application of a virtual displacement is a means to find
this first variation of H.

So given any real force, F, acting on a body to which we apply a virtual
displacement. If the virtual displacement at the location of and in the direction of F is
y o , then the force F does virtual work W F y o o = .

There are requirements on what is permissible as a virtual displacement. For
example, in the simple proof of the Principle of Virtual Work (to follow) it can be
seen that it is assumed that the directions of the forces applied to P remain
unchanged. Thus:
- virtual displacements must be small enough such that the force directions are
maintained.

The other very important requirement is that of compatibility:
- virtual displacements within a body must be geometrically compatible with
the original structure. That is, geometrical constraints (i.e. supports) and
member continuity must be maintained.

In summary, virtual displacements are not real, they can be physically impossible but
they must be compatible with the geometry of the original structure and they must be
small enough so that the original geometry is not significantly altered.

As the deflections usually encountered in structures do not change the overall
geometry of the structure, this requirement does not cause problems.
Dr. C. Caprani 18
2.3 Virtual Forces
So far we have only considered small virtual displacements and real forces. The
virtual displacements are arbitrary: they have no relation to the forces in the system,
or its actual deformations. Therefore virtual work applies to any set of forces in
equilibrium and to any set of compatible displacements and we are not restricted to
considering only real force systems and virtual displacements. Hence, we can use a
virtual force system and real displacements. Obviously, in structural design it is these
real displacements that are of interest and so virtual forces are used often.

A virtual force is a force imagined to be applied and is then moved through the actual
deformations of the body, thus causing virtual work.

So if at a particular location of a structure, we have a deflection, y, and impose a
virtual force at the same location and in the same direction of F o we then have the
virtual work W y F o o = .

Virtual forces must form an equilibrium set of their own. For example, if a virtual
force is applied to the end of a spring there will be virtual stresses in the spring as
well as a virtual reaction.
Dr. C. Caprani 19
2.4 Simple Proof using Virtual Displacements
We can prove the Principle of Virtual Work quite simply, as follows. Consider a
particle P under the influence of a number of forces
1
, ,
n
F F which have a resultant
force,
R
F . Apply a virtual displacement of y o to P, moving it to P, as shown:

The virtual work done by each of the forces is:

1 1 n n R R
W F y F y F y o o o o = + + =

Where
1
y o is the virtual displacement along the line of action of
1
F and so on. Now
if the particle P is in equilibrium, then the forces
1
, ,
n
F F have no resultant. That is,
there is no net force. Hence we have:

1 1
0 0
R n n
W y F y F y o o o o = = + + = (2.1)

Proving that when a particle is in equilibrium the virtual work of all the forces acting
on it sum to zero. Conversely, a particle is only in equilibrium if the virtual work
done during a virtual displacement is zero.
Dr. C. Caprani 20
2.5 Internal and External Virtual Work
Consider the spring we started with, as shown again below. Firstly it is unloaded.
Then a load, F, is applied, causing a deflection y. It is now in its equilibrium position
and we apply a virtual displacement, y o to the end of the spring, as shown:

A free body diagram of the end of the spring is:

Thus the virtual work done by the two forces acting on the end of the spring is:

W F y ky y o o o =

If the spring is to be in equilibrium we must then have:

0
0
W
F y ky y
F y ky y
F ky
o
o o
o o
=
=
=
=

Dr. C. Caprani 21
That is, the force in the spring must equal the applied force, as we already know.
However, if we define the following:
- External virtual work,
E
W F y o o = ;
- Internal virtual work,
I
W ky y o o = ;

We then have:

0
0
I E
W
W W
o
o o
=
=

Thus:

E I
W W o o = (2.2)

which states that the external virtual work must equal the internal virtual work for a
structure to be in equilibrium.

It is in this form that the Principle of Virtual Work finds most use.

Of significance also in the equating of internal and external virtual work, is that there
are no requirements for the material to have any particular behaviour. That is, virtual
work applies to all bodies, whether linearly-elastic, elastic, elasto-plastic, plastic etc.
Thus the principle has more general application than most other methods of analysis.

Internal and external virtual work can arise from either virtual displacements or
virtual forces.
Dr. C. Caprani 22
3. Application of Virtual Displacements
3.1 Rigid Bodies
Basis
Rigid bodies do not deform and so there is no internal virtual work done. Thus:

0
0
E I
i i
W
W W
F y
o
o o
o
=
=
=
(3.1)

A simple application is to find the reactions of statically determinate structures.
However, to do so, we need to make use of the following principle:

Principle of Substitution of Constraints
Having to keep the constraints in place is a limitation of virtual work. However, we
can substitute the restraining force (i.e. the reaction) in place of the restraint itself.
That is, we are turning a geometric constraint into a force constraint. This is the
Principle of Substitution of Constraints. We can use this principle to calculate
unknown reactions:

1. Replace the reaction with its appropriate force in the same direction (or sense);
2. Impose a virtual displacement on the structure;
3. Calculate the reaction, knowing 0 W o = .

Dr. C. Caprani 23
Reactions of Determinate and Indeterminate Structures
For statically determinate structures, removing a restraint will always render a
mechanism (or rigid body) and so the reactions of statically determinate structures are
easily obtained using virtual work. For indeterminate structures, removing a restraint
does not leave a mechanism and hence the virtual displacements are harder to
establish since the body is not rigid.

Dr. C. Caprani 24
Example 1
Determine the reactions for the following beam:

Following the Principle of Substitution of Constraints, we replace the geometric
constraints (i.e. supports), by their force counterparts (i.e. reactions) to get the free-
body-diagram of the whole beam:

Next, we impose a virtual displacement on the beam. Note that the displacement is
completely arbitrary, and the beam remains rigid throughout:
Dr. C. Caprani 25

In the above figure, we have imposed a virtual displacement of
A
y o at A and then
imposed a virtual rotation of
A
ou about A. The equation of virtual work is:

0
0
E I
i i
W
W W
F y
o
o o
o
=
=
=

Hence:

0
A A C B B
V y P y V y o o o + =

Relating the virtual movements of points B and C to the virtual displacements gives:

( ) ( )
0
A A A A B A A
V y P y a V y L o o ou o ou + + + =

And rearranging gives:

( ) ( )
0
A B A B A
V V P y V L Pa o ou + + =
Dr. C. Caprani 26
And here is the power of virtual work: since we are free to choose any value we want
for the virtual displacements (i.e. they are completely arbitrary), we can choose
0
A
ou = and 0
A
y o = , which gives the following two equations:

( )
0
0
A B A
A B
V V P y
V V P
o + =
+ =

( )
0
0
B A
B
V L Pa
V L Pa
ou =
=

But the first equation is just 0
y
F =
whilst the second is the same as

about 0 M A=
. Thus equilibrium equations occur naturally within the virtual

work framework. Note also that the two choices made for the virtual displacements
correspond to the following virtual displaced configurations: Draw them

Dr. C. Caprani 27
Example 2
For the following truss, calculate the reaction
C
V :

Firstly, set up the free-body-diagram of the whole truss:

Dr. C. Caprani 28
Next we release the constraint corresponding to reaction
C
V and replace it by the
unknown force
C
V and we apply a virtual displacement to the truss to get:

Hence the virtual work done is:

0
10 0
2
5 kN
E I
C
C
W
W W
y
V y
V
o
o o
o
o
=
=
+ =
=

Note that the reaction is an external force to the structure, and that no internal virtual
work is done since the members do not undergo virtual deformation. The truss rotates
as a rigid body about the support A.

Problem:
Find the horizontal reactions of the truss.
Dr. C. Caprani 29
3.2 Deformable Bodies
Basis
For a virtual displacement we have:

0
E I
i i i i
W
W W
F y P e
o
o o
o o
=
=
=
}

In which, for the external virtual work,
i
F represents an externally applied force (or
moment) and
i
y o its virtual displacement. And for the internal virtual work,
i
P
represents the internal force in member i and
i
e o its virtual deformation. Different
stress resultants have different forms of internal work, and we will examine these.
Lastly, the summations reflect the fact that all work done must be accounted for.

Remember in the above, each the displacements must be compatible and the forces
must be in equilibrium, summarized as:

i i i i
F y P e o o =

Set of forces in
equilibrium
Set of compatible
displacements
Dr. C. Caprani 30
These displacements are completely arbitrary (i.e. we can choose them to suit our
purpose) and bear no relation to the forces above. For example, the simple truss
shown has a set of forces in equilibrium along with its actual deformed shape. Also
shown is an arbitrary set of permissible compatible displacements:

Loading and dimensions Equilibrium forces and actual deformation

Compatible set of displacements

Dr. C. Caprani 31
Internal Virtual Work by Axial Force
Members subject to axial force may have the following:
- real force by a virtual displacement:

I
W P e o o =

- virtual force by a real displacement:

I
W e P o o =

We have avoided integrals over the length of the member since we will only consider
prismatic members.
Dr. C. Caprani 32
Internal Virtual Work in Bending
The internal virtual work done in bending is one of:
- real moment by a virtual curvature:

I
W M o ok =

- virtual moment by a real curvature:

I
W M o k o =

The above expressions are valid at a single position in a beam.

When virtual rotations are required along the length of the beam, the easiest way to
do this is by applying virtual loads that in turn cause virtual moments and hence
virtual curvatures. We must sum all of the real moments by virtual curvatures along
the length of the beam, hence we have:

0 0
0
L L
I x x x x
L
x
x
x
W M dx M dx
M
M dx
EI
o ok k o
o
= =
=
} }
}

Dr. C. Caprani 33
Internal Virtual Work in Shear
At a single point in a beam, the shear strain, , is given by:

v
V
GA
=

where V is the applied shear force, G is the shear modulus and
v
A is the cross-section
area effective in shear which is explained below. The internal virtual work done in
shear is thus:
- real shear force by a virtual shear strain:

I
v
V
W V V
GA
o
o o = =

- virtual shear force by a real shear strain:

I
v
V
W V V
GA
o o o = =

In these expressions, the area of cross section effective in shear comes about because
the shear stress (and hence strain) is not constant over a cross section. The stress
v
V A is an average stress, equivalent in work terms to the actual uneven stress over
the full cross section, A. We say
v
A A k = where k is the shear factor for the cross
section. Some values of k are 1.2 for rectangular sections, 1.1 for circular sections,
and 2.0 for thin-walled circular sections.

The above expressions are valid at a single position in a beam and must be integrated
along the length of the member as was done for moments and rotations.
Dr. C. Caprani 34
Internal Virtual Work in Torsion
At a single point in a member, the twist, | , is given by:

T
GJ
| =

where T is the applied torque, G is the shear modulus, J is the polar moment of
inertia. The internal virtual work done in torsion is thus:
- real torque by a virtual twist:

I
T
W T T
GJ
o
o o| = =

- virtual torque by a real twist:

I
T
W T T
GJ
o o | o = =

Once again, the above expressions are valid at a single position in a beam and must
be integrated along the length of the member as was done for moments and rotations.

Note the similarity between the expressions for the four internal virtual works.

Dr. C. Caprani 35
Example 3
For the beam of Example 1 (shown again), find the bending moment at C.

To solve this, we want to impose a virtual displacement configuration that only
allows the unknown of interest, i.e.
C
M , to do any work. Thus choose the following:

Since portions AC and CB remain straight (or unbent) no internal virtual work is done
in these sections. Thus the only internal work is done at C by the beam moving
through the rotation at C. Thus:
Dr. C. Caprani 36

0
E I
C C C
W
W W
P y M
o
o o
o ou
=
=
=

Bu the rotation at C is made up as:

C CA CB
C C
C
y y
a b
a b
y
ab
ou ou ou
o o
o
= +
= +
+
| |
=
|
\ .

But a b L + = , hence:

C C C
C
L
P y M y
ab
Pab
M
L
o o
| |
=
|
\ .
=

We can verify this from the reactions found previously: ( )
C B
M V b Pa L b = = .

Dr. C. Caprani 37
Sign Convention for Rotations
When imposing a virtual rotation, if the side that is already in tension due to the real
force system elongates, then positive virtual work is done.

Dr. C. Caprani 38
Example 4
Calculate the force
1
F in the truss shown:

To do this, we introduce a virtual displacement along the length of member 1. We do
this so that member 2 does not change length during this virtual displacement and so
does no virtual work. Note also that compatibility of the truss is maintained. For
example, the members still meet at the same joint, and the support conditions are met.

Dr. C. Caprani 39
The virtual work done is then:

1
1
0
10
2
10
5 2 kN
2
E I
W
W W
y
F y
F
o
o o
o
o
=
=
=
= =

Note some points on the signs used:
1. Negative external work is done because the 10 kN load moves upwards; i.e. the
reverse direction to its action.
2. We assumed member 1 to be in compression but then applied a virtual
elongation which lengthened the member thus reducing its internal virtual
work. Hence negative internal work is done.
3. We initially assumed
1
F to be in compression and we obtained a positive
answer confirming our assumption.

Problem:
Investigate the vertical and horizontal equilibrium of the loaded joint by considering
vertical and horizontal virtual displacements separately.

Dr. C. Caprani 40
3.3 Problems
1. For the truss shown, calculate the vertical reaction at C and the forces in the
members, using virtual work.

2. For the truss shown, find the forces in the members, using virtual work:

Dr. C. Caprani 41

3. Using virtual work, calculate the reactions for the beams shown, and the
bending moments at salient points.

Beam 1

Beam 2

Beam 3

Dr. C. Caprani 42
4. Application of Virtual Forces
4.1 Basis
When virtual forces are applied, we have:

0
E I
i i i i
W
W W
y F e P
o
o o
o o
=
=
=

And again note that we have an equilibrium set of forces and a compatible set of
displacements:

In this case the displacements are the real displacements that occur when the structure
is in equilibrium and the virtual forces are any set of arbitrary forces that are in
equilibrium.

i i i i
y F e P o o =

Set of compatible
displacements
Set of forces in
equilibrium

Dr. C. Caprani 43
4.2 Deflection of Trusses
Illustrative Example
For the truss shown below, we have the actual displacements shown with two
possible sets of virtual forces.

Actual Loading and (Compatible) Displacements
Dr. C. Caprani 44

Virtual Force (Equilibrium) Systems

In this truss, we know we know:
1. The forces in the members (got from virtual displacements or statics);
2. Thus we can calculate the member extensions,
i
e as:

i
i
PL
e
EA
| |
=
|
\ .

3. Also, as we can choose what our virtual force F o is (usually unity), we have:

0
E I
i i i i
i
i
W
W W
y F e P
PL
y F P
EA
o
o o
o o
o o
=
=
=
| |
=
|
\ .

4. Since in this equation, y is the only unknown, we can calculate the deflection
of the truss.

Dr. C. Caprani 45
Example 5
Given that E = 200 kN/mm
2
and A = 100 mm
2
for each member, for the truss shown
below, calculate the vertical and horizontal deflection of joint B.

In these problems we will always choose 1 F o = . Hence we apply a unit virtual force
to joint B. We apply the force in the direction of the deflection required. In this way
no work is done along deflections that are not required. Hence we have:

For horizontal deflection For vertical deflection

Dr. C. Caprani 46
The forces and elongations of the truss members are:
- Member AB:
50 kN
50 5000
200 100
12.5 mm
AB
AB
P
e
= +
+
| |
=
|
\ .
= +

- Member BC:
30 kN
30 3000
200 100
4.5 mm
BC
BC
P
e
=

| |
=
|
\ .
=

Note that by taking tension to be positive, elongations are positive and contractions
are negative.

Horizontal Deflection:
( ) ( )
0
1 12.5 0 4.5 1
4.5 mm
E I
i i i i
AB BC
W
W W
y F e P
y
y
o
o o
o o
=
=
=
= + +
=

Vertical Deflection:
5 3
1 12.5 4.5
4 4
18.4 mm
i i i i
AB BC
y F e P
y
y
o o =
| | | |
= +
| |
\ . \ .
= +

Note that the sign of the result indicates whether the deflection occurs in the same
direction as the applied force. Hence, joint B moves 4.5 mm to the left.

Dr. C. Caprani 47
Example 6
Determine the vertical and horizontal deflection of joint D of the truss shown. Take E
= 200 kN/mm
2
and member areas, A = 1000 mm
2
for all members except AE and BD
where A = 1000\2 mm
2
.

The elements of the virtual work equation are:
- Compatible deformations: The actual displacements that the truss undergoes;
- Equilibrium set: The external virtual force applied at the location of the
required deflection and the resulting internal member virtual forces.

Firstly we analyse the truss to determine the member forces in order to calculate the
actual deformations of each member:

Dr. C. Caprani 48
Next, we apply a unit virtual force in the vertical direction at joint D. However, by
linear superposition, we know that the internal forces due to a unit load will be 1/150
times those of the 150 kN load.

For the horizontal deflection at D, we apply a unit horizontal virtual force as shown:

Equations of Virtual Work
0
1
0
2
0
1
1
E I
i i i i
DV i
i
DH i
i
W
W W
y F e P
P L
y P
EA
P L
y P
EA
o
o o
o o
o
o
=
=
=
| |
=
|
\ .
| |
=
|
\ .

In which:
-
0
P are the forces due to the 150 kN load;
-
1
P o are the virtual forces due to the unit virtual force applied in the vertical
direction at D:

0
1
150
P
P o =
-
2
P o are the virtual forces due to the unit virtual force in the horizontal
direction at D.

Dr. C. Caprani 49
Using a table is easiest because of the larger number of members:

Member
L A
0
P
0
1
150
P
P o =
2
P o
0
1
P L
P
A
o
| |
|
\ .

0
2
P L
P
A
o
| |
|
\ .

(mm) (mm
2
) (kN) (kN) (kN) (kN/mm)kN (kN/mm)kN
AB
AE
AF
BC
BD
BE
CD
DE
EF
2000
20002
2000
2000
20002
2000
2000
2000
2000
1000
10002
1000
1000
10002
1000
1000
1000
1000
+150
+1502
0
0
+1502
-150
0
-150
-300
+1
+12
0
0
+12
-1
0
-1
-2
0
0
0
0
0
0
0
+1
+1
+300
+600
0
0
+600
+300
0
+300
+1200
0
0
0
0
0
0
0
-300
-600
=

3300 -900

E is left out because it is common. Returning to the equations, we now have:
0
1
1
1
3300
16.5 mm
200
DV i
i
DV
P L
y P
E A
y
o
| |
=
|
\ .
+
= = +

Which indicates a downwards deflection and for the horizontal deflection:
0
2
1
1
900
4.5 mm
200
DH i
i
DH
P L
y P
E A
y
o
| |
=
|
\ .
= =

The sign indicates that it is deflecting to the left.
Dr. C. Caprani 50
4.3 Deflection of Beams and Frames
Example 7
Using virtual work, calculate the deflection at the centre of the beam shown, given
that EI is constant.

To calculate the deflection at C, we will be using virtual forces. Therefore the two
relevant sets are:
- Compatibility set: the actual deflection at C and the rotations that occur along
the length of the beam;
- Equilibrium set: a unit virtual force applied at C which is in equilibrium with
the internal virtual moments it causes.

Compatibility Set:
The external deflection at C is what is of interest to us. To calculate the rotations
along the length of the beam we have:

0
L
x
x
M
dx
EI
k =
}

Hence we need to establish the bending moments along the beam:
Dr. C. Caprani 51

For AC the bending moment is given by (and similarly for B to C):

2
x
P
M x =

Equilibrium Set:
As we choose the value for 1 F o = , we are only left to calculate the virtual moments:

For AC the internal virtual moments are given by:

1
2
x
M x o =
Dr. C. Caprani 52
Virtual Work Equation

0
E I
i i i i
W
W W
y F M
o
o o
o k o
=
=
=

Substitute in the values we have for the real rotations and the virtual moments, and
use the fact that the bending moment diagrams are symmetrical:

2
0
2
0
2
2
0
2
3
0
3
3
1 2
2 1
2 2
2
4
2 3
6 8
48
L
x
x
L
L
L
M
y M dx
EI
P
y x x dx
EI
P
x dx
EI
P x
EI
P L
EI
PL
EI
o
(
=
(

| | | |
=
| |
\ . \ .
=
(
=
(

=
=
}
}
}

Which is a result we expected.

Dr. C. Caprani 53
Example 8
Find the vertical and horizontal displacement of the prismatic curved bar cantilever
shown:

Even though it is curved, from the statics, the bending moment at any point of course
remains force by distance, so from the following diagram:

Dr. C. Caprani 54
At any angle | , we therefore have:

( )
( )
cos
1 cos
M P R R
PR
|
|
|
=
=

To find the displacements, we follow our usual procedure and place a unit load at the
location of, and in the direction of, the required displacement.

Vertical Displacement
The virtual bending moment induced by the vertical unit load shown, is related to that
for P and is thus:

( )
1 cos M R
|
o | =

Thus our virtual work equations are:

Dr. C. Caprani 55

2
0
0
1
E I
BV
W
W W
M ds
M
M Rd
EI
t
|
|
o
o o
o k o
o |
=
=
=
=
}
}

In which we have used the relation ds R d| = to change the integration from along
the length of the member to around the angle | . Next we introduce our equations for
the real and virtual bending moments:

( )
( )
( )
( )
2
0
2
3
2
0
3 3
1 cos
1 1 cos
1 cos
3 4 5.42
BV
PR
R Rd
EI
PR
d
EI
PR PR
EI EI
t
t
|
o | |
| |
t
=
=
= ~
}
}

Horizontal Displacement
In this case, the virtual bending moment is:

sin M R
|
o | =

Dr. C. Caprani 56
Thus the virtual work equations give:

( )
( )
2
0
2
3
0
2
3
0
3 3
1 cos
1 sin
sin sin cos
sin2
sin
2
1
2 2
BH
PR
R Rd
EI
PR
d
EI
PR
d
EI
PR PR
EI EI
t
t
t
|
o | |
| | | |
|
| |
=
=
| |
=
|
\ .
| |
= =
|
\ .
}
}
}

Dr. C. Caprani 57
4.4 Integration of Bending Moments
We are often faced with the integration of being moment diagrams when using virtual
work to calculate the deflections of bending members. And as bending moment
diagrams only have a limited number of shapes, a table of volume integrals is used:

This table is at the back page of these notes for ease of reference.

Example 9
Using the table of volume integrals, verify the answer to Example 7.

In this case, the virtual work equation becomes:

| | | |
2
0
3
1 2
2
shape shape
2 1
3 4 4 2
48
L
x
x
x x
M
y M dx
EI
y M M
EI
PL L L
EI
PL
EI
o
o
(
=
(

=
(
=
(

=
}

In which the formula
1
3
jkl is used from the table.
Dr. C. Caprani 58
Example 10 Summer 07 Part (a)
For the frame shown, determine the horizontal deflection of joint C. Neglect axial
effects in the members and take
3 2
36 10 kNm EI = .

FIG. Q3(a)
A
B
C
2
0

k
N
/
m
100 kN
D

Firstly we establish the real bending moment diagram:

Next, as usual, we place a unit load at the location of, and in the direction of, the
required displacement:
Dr. C. Caprani 59

Now we have the following for the virtual work equation:

0
1
E I
CH
W
W W
M ds
M
M ds
EI
o
o o
o k o
o
=
=
=
=
}
}

Next, using the table of volume integrals, we have:

( )( )
( )
( )( )( )
1 1 1
440 2 4 2 2 160 2 440 4
3 6
1659.3 1386.7
3046
AB BC
M
M ds
EI EI
EI EI
EI
o

( (
= + +
`
( (

)
= +
=
}

Hence:

3
3
3046 3046
1 10 84.6 mm
36 10
CH
EI
o = = =

Dr. C. Caprani 60
4.5 Problems
1. Determine the vertical deflection of joint E of the truss shown. Take E = 200
kN/mm
2
2
for all members except AC where A =
1000\2 mm
2
. (Ans. 12.2 mm).

FIG. Q3

60 kN
A
B C
D
E

2. Determine the vertical and horizontal deflection of joint C of the truss shown.
Take E = 10 kN/mm
2
2
for all members except
AC where A = 1000\2 mm
2
and CE where A = 2500 mm
2
.
(Ans. 32/7 mm horizontal, 74/7 mm vertical).

FIG. Q3(a)
A
B C
D
E

50 kN
3
4

Dr. C. Caprani 61
3. Determine the horizontal deflection of joint A and the vertical deflection of joint
B of the truss shown. Take E = 200 kN/mm
2
2

for all members except BD where A = 1000\2 mm2 and AB where A = 2500
mm
2
. (Ans. 15.3 mm; 0 mm)
A
B E
C D
90 kN
FIG. Q2(a)

4. Determine the vertical and horizontal deflection of joint C of the truss shown.
Take E = 200 kN/mm
2
2
AC where A = 150\2 mm2 and CD where A = 250 mm
2
. (Ans. 8.78 mm; 0.2 mm)

FIG. Q3
A
B
C
D
100 kN

Dr. C. Caprani 62
5. Determine the vertical deflection of joint D of the truss shown. Take E = 200
kN/mm
2
2
for all members except BC where A =
1000\2 mm
2
. (Ans. 24 mm)

C
B
A
D
E
FIG. Q3

90 kN

6. Verify that the deflection at the centre of a simply-supported beam under a
uniformly distributed load is given by:

4
5
384
C
wL
EI
o =

7. Show that the flexural deflection at the tip of a cantilever, when it is subjected to
a point load on the end is:

3
3
B
PL
EI
o =

8. Show that the rotation that occurs at the supports of a simply supported beam
with a point load in the middle is:

2
16
C
PL
EI
u =

Dr. C. Caprani 63
9. Show that the vertical and horizontal deflections at B of the following frame are:

3 3
;
2 4
By Bx
PR PR
EI EI
o o = =

10. For the frame of Example 10, using virtual work, verify the following
displacements in which the following directions are positive: vertical upwards;
horizontal to the right, and; clockwise:
- Rotation at A is
1176.3
EI
;
- Vertical deflection at B is
3046
EI
;
- Horizontal deflection at D is
8758.7
EI
;
- Rotation at D is
1481.5
EI
.

Dr. C. Caprani 64
11. For the cantilever beam below, show that the vertical deflection of B, allowing
for both flexural and shear deformation, is
3
3
v
PL PL
EI GA
+ .

12. For the beam below, show that the vertical deflection of D is zero and that that of
C is
1850 3
EI
:

Dr. C. Caprani 65
Problems 13 to 17 are considered Genius Level!

13. Show that the vertical deflection at the centre of a simply supported beam subject
to UDL, allowing for both flexural and shear deformation, is
4 2
5
384 8
v
wL wL
EI GA
+ .

14. For the frame shown, show that the vertical deflection of point D is
4594
EI
+.
Neglect axial deformation and take
3 2
120 10 kNm EI = .

Dr. C. Caprani 66
15. For the frame shown, determine the vertical deflection of C and the rotation at D
about the global x axis. Take the section and material properties as follows:
2 2 2
6 4 6 4
6600 mm 205 kN/mm 100 kN/mm
36 10 mm 74 10 mm .
A E G
I J
= = =
= =

(Ans. 12.23 mm|; 0.0176 rads ACW)

16. Determine the horizontal and vertical deflections at B in the following frame.
Neglect shear and axial effects.

(Ans.
4
2
Bx
wR
EI
o = ;
4
3
2
By
wR
EI
t
o = |)
Dr. C. Caprani 67

17. For the frame shown, neglecting shear effects, determine the vertical deflection of
points B and A.

(Ans.
( )
3
3
By
P wa b
EI
o
+
= +;
( )
3
3 4 2
1
3 8 3 2
Ay
P wa b
Pa wa wa ab
Pa
EI GJ
o
( + | |
= + + + + +
( |
\ .

)

Dr. C. Caprani 68
5. Virtual Work for Indeterminate Structures
5.1 General Approach
Using the concept of compatibility of displacement, any indeterminate structure can
be split up into a primary and reactant structures. In the case of a one-degree
indeterminate structure, we have:

Final = Primary + Reactant

At this point the primary structure can be analysed. However, we further break up the
reactant structure, using linear superposition:

Reactant = Multiplier Unit Reactant

We summarize this process as:

0 1
M M M o = +
- M is the force system in the original structure (in this case moments);
-
0
M is the primary structure force system;
-
1
M is the unit reactant structure force system.
Dr. C. Caprani 69
For a truss, the procedure is the same:

Final = Primary + Reactant

Reactant = Multiplier Unit Reactant

The final system forces are:

0 1
P P P o = +

The primary structure can be analysed, as can the unit reactant structure. Therefore,
the only unknown is the multiplier, o .

We use virtual work to calculate the multiplier o .

Dr. C. Caprani 70
5.2 Using Virtual Work to Find the Multiplier
We must identify the two sets for use:
- Displacement set: We use the actual displacements that occur in the real
structure;
- Equilibrium set: We use the unit reactant structures set of forces as the
equilibrium set. We do this, as the unit reactant is always a determinate structure
and has a configuration similar to that of the displacement set.

A typical example of these sets for an indeterminate structure is shown:

Compatibility Set Equilibrium Set
Internal and external real
displacements linked through
constitutive relations
Internal (bending moments) and external
virtual forces are linked through
equilibrium

Note that the Compatibility Set has been chosen as the real displacements of the
structure. The virtual force (equilibrium) set chosen is a statically determine sub-
structure of the original structure. This makes it easy to use statics to determine the
internal virtual forces from the external virtual force. Lastly, note that the external
virtual force does no virtual work since the external real displacement at B is zero.
This makes the right hand side of the equation zero, allowing us to solve for the
multiplier.

Dr. C. Caprani 71
The virtual work equation (written for trusses) gives:

1
0
0 1
E I
i i i i
i
i
W
W W
y F e P
PL
P
EA
o
o o
o o
o
=
=
=
| |
=
|
\ .

There is zero external virtual work. This is because the only virtual force applied is
internal; no external virtual force applied. Also note that the real deformations that
occur in the members are in terms of P, the unknown final forces. Hence, substituting
0 1
P P P o o = + (where o is now used to indicate virtual nature):

( )
( )
0 1
1
0 1
1 1
2
1
0 1
0
0
i
i
i i
i i
i i
i i
i i
P P L
P
EA
P L P L
P P
EA EA
P L
P P L
EA EA
o o
o
o
o o o
o
o
o
| |
+
| =
|
\ .
| | | |
= +
| |
\ . \ .

= +

For beams and frames, the same development is:

0
0 1
E I
i i i i
i
i
W
W W
y F M
M
M
EI
o
o o
o k o
o
=
=
=
=
}
}

Where again there is no external displacement of the virtual force.
Dr. C. Caprani 72

Note that we use the contour integral symbol | |
}
simply to indicate that we integrate
around the structure, accounting for all members in the beam/frame (i.e. integrate
along the length of each member separately, and then sum the results).

Also, substitute
0 1
M M M o o = + to get:

( )
( )
0 1
1
2
1
0 1
0
0
i
i
i
i
i i
M M
M dx
EI
M
M M
dx dx
EI EI
o o
o
o
o
o
+
=
= +
}
} }

Thus in both bases we have a single equation with only one unknown, o . We can
establish values for the other two terms and then solve for o and the structure as a
whole.
Dr. C. Caprani 73
5.3 Basic Example Propped Cantilever
Example 11
Calculate the bending moment diagram for the following prismatic propped
cantilever. Find also the deflection under the point load.

Solution to be done in class.
Dr. C. Caprani 74
5.4 Indeterminate Trusses
Example 12
Calculate the forces in the truss shown and find the horizontal deflection at C. Take
EA to be 1010
4
kN for all members.

Solve for the Truss Member Forces
Choose member AC as the redundant:

Next analyse for the
0
P and
1
P force systems.

Dr. C. Caprani 75

P
0
System P
1
System

Using a table:

Member
L
0
P
1
P o
0 1
P P L o ( )
2
1
P L o
(mm) (kN) (kN) 10
4
10
4

AB
BC
CD
AC
BD
4000
3000
4000
5000
5000
+40
0
0
0
-50
- 4/5
- 3/5
- 4/5
1
1
-12.8
0
0
0
-25
0.265
0.108
0.256
0.5
0.5
=

-37.8 1.62

Hence:
( )
2
1
0 1
4 4
0
37.8 10 1.62 10
i i
i i
i i
P L
P P L
EA EA
EA EA
o
o
o
o

= +

= +

And so
Dr. C. Caprani 76
37.8
23.33
1.62
o
= =

The remaining forces are obtained from the compatibility equation:

Member
0
P
1
P o
0 1
P P P o o = +
(kN) (kN) (kN)
AB
BC
CD
AC
BD
+40
0
0
0
-50
- 4/5
- 3/5
- 4/5
1
1
21.36
-14
-18.67
23.33
-26.67

Note that the redundant always has a force the same as the multiplier.

Calculate the Horizontal Deflection at C
To calculate the horizontal deflection at C, using virtual work, the two relevant sets
are:
- Compatibility set: the actual deflection at C and the real deformations that
occur in the actual structure;
- Equilibrium set: a horizontal unit virtual force applied at C to a portion of the
actual structure, yet ensuring equilibrium.

We do not have to apply the virtual force to the full structure. Remembering that the
only requirement on the virtual force system is that it is in equilibrium; choose the
force systems as follows:

Dr. C. Caprani 77

Thus we have:

4 4
0
1
23.33 5000 5 18.67 4000 4
10 10 3 10 10 3
2.94 mm
E I
i i i i
CH i
i
CH
W
W W
y F e P
PL
y P
EA
y
o
o o
o o
o
=
=
=
| |
=
|
\ .

| | | |
= +
| |

\ . \ .
=

Because we have chosen only two members for our virtual force system, only these
members do work and the calculation is greatly simplified.
Dr. C. Caprani 78
5.5 Indeterminate Frames
Example 13
For the frame shown, calculate the reactions and draw the bending moment diagram.
Determine the horizontal deflection at joint C. Take EI = 4010
3
kNm
2
, constant for
all members.

Break the frame up into its reactant and primary structures:

M System = M
0
System + o M
1
System

Establish the
0
M and
1
M force systems:

Dr. C. Caprani 79

M
0
System M
1
System

Apply the virtual work equation:
( )
2
1
0 1
0
i
i
i i
M
M M
dx dx
EI EI
o
o
o
= +
} }

We will be using the table of volume integrals to quicken calculations. Therefore we
can only consider lengths of members for which the correct shape of bending moment
diagram is available. Also, we must choose sign convention: we consider tension on
the outside of the frame to be positive.

As each term has several components we consider them separately:

Term 1 -
0 1
i
i
M M
dx
EI
o
}
:
- AD: we have graphically, and from the volume integral table:

( ) ( )( )
1 2
1 1
6 600 360 4 11520
2 2
j k k l + = + =
Dr. C. Caprani 80
- DB: Similarly:

( )( ) 6 360 4 8640 jkl = =
- BC:

( )( )
1 1
6 360 6 3240
4 4
jkl = =
Hence:
0 1
23400
i
i
M M
dx
EI EI
o
=
}

Term 2:
( )
| |
( )( ) ( )( )
2
1
1 1
3
1 1
6 6 8 6 6 6
3
360
i
AB
BC i
AB
BC
M
dx jkl jkl
EI EI
EI
EI
o

(
= +
`
(

)

(
= + (
`

(

)
=
}

Note that Term 2 is always easier to calculate as it is only ever made up of straight
line bending moment diagrams.

Thus we have:

Dr. C. Caprani 81
( )
2
1
0 1
0
23400 360
i
i
i i
M
M M
dx dx
EI EI
EI EI
o
o
o
o
= +
= +
} }

And so:

23400
65.0
360
o = =

Thus the vertical reaction at C is 65.0 kN upwards. With this information we can
solve for the moments using
0 1
M M M o = + (or by just using statics) and the
remaining reactions to get:

Dr. C. Caprani 82
To calculate the horizontal deflection at C using virtual work, the two relevant sets
are:
- Compatibility set: the actual deflection at C and the real deformations
(rotations) that occur in the actual structure;
- Equilibrium set: a horizontal unit virtual force applied at C to a determinate
portion of the actual structure.

Choose the following force system as it is easily solved:

Thus we have:

0
1
E I
i i i i
x
CH x
W
W W
y F M
M
y M dx
EI
o
o o
o u o
o
=
=
=
(
=
(

}
}

Dr. C. Caprani 83
The real bending moment diagram, M, is awkward to use with the integral table.
Remembering that
0 1
M M M o = + simplifies the calculation by using:

To give:

And so using the table formulae, we have:

Dr. C. Caprani 84

( )( ) ( )( )
( )( ) ( )( )
( )( ) ( )( )
1 1
4 6 4 2 360 600 4 4
2 6
1 1
4 8 6 4 4 8 360 4
2 2
1 1
8 6 6 8 360 6
3 4
x
x
AD
DB
BC
M
M dx
EI
o o
o
o
( (
= + + +
( (

(
+ +
(

(
+
(

}

Which gives us:

( )
1
288 16480
2240
0.056
56 mm
x
Cx x
M
M dx
EI
EI
EI
o o
o
(
=
(

= + (

=
=
=
}

The answer is positive, indicating that the structure moves to the right at C: the same
direction in which the unit virtual force was applied.

We could have chosen any other statically determinate sub-structure and obtained the
same result. Some sub-structures will make the analysis easier to perform.

Exercise:
Verify that the same deflection is obtained by using a sub-structure obtained by
removing the vertical support at C and applying the virtual force. Does this sub-
structure lead to an easier analysis for the deflection?
Dr. C. Caprani 85
Example 14
For the frame shown, calculate the reactions and draw the bending moment diagram.
Determine the vertical deflection at joint A. Take EI = 2010
3
kNm
2
, constant for all
members.

To be done in class.

Dr. C. Caprani 86
5.6 Continuous Beams
Basis
In this section we will introduce structures that are more than 1 degree statically
indeterminate. We do so to show that virtual work is easily extensible to multiply-
indeterminate structures, and also to give a method for such beams that is easily
worked out, and put into a spreadsheet.

Consider the example 3-span beam. It is 2 degrees indeterminate, and so we introduce
2 hinges (i.e. moment releases) at the support locations and unit reactant moments in
their place, as shown:

Dr. C. Caprani 87
Alongside these systems, we have their bending moment diagrams:

Using the idea of the multiplier and superposition again, we can see that:

0 1 2
1 2
M M M M o o o o = + +

The virtual work equation is:

0
0 1
E I
W
W W
M
o
o o
k o
=
=
=
}

Dr. C. Caprani 88
There is no external virtual work done since the unit moments are applied internally.
Since we have two virtual force systems in equilibrium and one real compatible
system, we have two equations:

1
0
M
M
EI
o =
}
and
2
0
M
M
EI
o =
}

For the first equation, expanding the expression for the real moment system, M:

( )
0 1 2
1 2
1
0 1 1 1 2 1
1 2
0
0
M M M
M
EI
M M M M M M
EI EI EI
o o o o
o
o o o o o
o o
+ +
=

+ + =
}
} } }

In which weve dropped the contour integral it being understood that we sum for all
members. Similarly for the second virtual moments, we have:

0 2 1 2 2 2
1 2
0
M M M M M M
EI EI EI
o o o o o
o o

+ + =
} } }

Thus we have two equations and so we can solve for
1
o and
2
o . Usually we write
this as a matrix equation:

0 1 1 1 2 1
1
0 2 1 2 2 2
2
0
0
M M M M M M
EI EI EI
M M M M M M
EI EI EI
o o o o o
o
o
o o o o o
(
(

( + =
` ` `
( ) )

(

)
} } }
} } }

Each of the integral terms is easily found using the integral tables, and the equation
solved.
Dr. C. Caprani 89

Similarly the (simpler) equation for a 2-span beam is:

0 1 1 1
1
0
M M M M
EI EI
o o o
o

+ =
} }

And the equation for a 4-span beam is:

0 1 1 1 2 1
1
0 2 1 2 2 2 3 2
2
0 3 2 3 3 3 3
0
0
0
0
0
M M M M M M
EI EI EI
M M M M M M M M
EI EI EI EI
M M M M M M
EI EI EI
o o o o o
o
o o o o o o o
o
o
o o o o o
(
(
(

(
+ =
` ` `
(

(
) )

(

(
)
} } }
} } } }
} } }

The diagram shows why there are no terms involving
3 1
M M o o , and why it is only
adjacent spans that have non-zero integrals:

Dr. C. Caprani 90
Example 15
Using virtual work, find the bending moment diagram for the following beam:

Proceeding as described above, we introduce releases (hinges) at the support points,
apply the unit virtual moments, and find the corresponding bending moment
diagrams:

Next we need to evaluate each term in the matrix virtual work equation. Well take
the two hard ones first:
Dr. C. Caprani 91
-
( )( )( ) ( )( )( )
( )
0 1
1 1 1 1
50 1 4 2 45 1 6
6 2 3
1 95
50 45
AB BC
M M
EI EI EI
EI EI
o
( (
= + +
( (

= =
}

Note that since
1
0 M o = for span CD, there is no term for it above. Similarly, for the
following evaluation, there will be no term for span AB. Note also the 2EI term for
member BC this could be easily overlooked.

-
( )( )( ) ( )( )( )
( )
0 2
1 1 1 1
45 1 6 93.75 1 5
2 3 3
1 201.25
45 156.25
BC CD
M M
EI EI EI
EI EI
o
( (
= +
( (

= =
}

The following integrals are more straightforward since they are all triangles:

-
( )( )( ) ( )( )( )
1 1
1 1 1 1 2.333
1 1 4 1 1 6
3 2 3
AB BC
M M
EI EI EI EI
o o
( (
= + =
( (

}

-
( )( )( )
2 1
1 1 0.5
1 1 6
2 6
BC
M M
EI EI EI
o o
(
= =
(

}

-
1 2
0.5 M M
EI EI
o o
=
}
, since it is equal to
2 1
M M
EI
o o
}
by the commutative
property of multiplication.

- ( )( )( ) ( )( )( )
2 2
1 1 1 1 2.667
1 1 6 1 1 5
2 3 3
BC CD
M M
EI EI EI EI
o o
( (
= + =
( (

}

With all the terms evaluated, enter them into the matrix equation:
Dr. C. Caprani 92

0 1 1 1 2 1
1
0 2 1 2 2 2
2
0
0
M M M M M M
EI EI EI
M M M M M M
EI EI EI
o o o o o
o
o
o o o o o
(
(

( + =
` ` `
( ) )

(

)
} } }
} } }

To give:

1
2
95 2.333 0.5 0
1 1
201.25 0.5 2.667 0 EI EI
o
o
(
+ =
` ` `
(
) ) )

And solve, as follows:

( )
1
2
1
2
2.333 0.5 95
0.5 2.667 201.25
2.667 0.5 95
1
0.5 2.333 201.25 2.333 2.667 0.5 0.5
25.57
70.67
o
o
o
o
(
=
` `
(
) )
(
=
` `
(

) )

=
`
)

Now using our superposition equation for moments,
0 1 2
1 2
M M M M o o o o = + + ,
we can show that the multipliers are just the hogging support moments:

0 25.57 1 70.67 0 25.57 kNm
0 25.57 0 70.67 1 70.67 kNm
B
C
M
M
= + + =
= + + =

From these we get the final BMD:

Dr. C. Caprani 93

Dr. C. Caprani 94
Spreadsheet Solution
A simple spreadsheet for a 3-span beam with centre-span point load and UDL capabilities, showing Example 1, is:

Span 1 Span 2 Span 3 VW Equations
Length m 4 6 5
EI x 1 2 1 -95.000 + 2.333 0.500 * alpha1 = 0
PL kN 50 0 0 -201.250 0.500 2.667 alpha2 0
UDL kN/m 0 10 30
BMDs PL 50 0 0 2.333 0.500 * alpha1 = 95.000
UDL 0 45 93.75 0.500 2.667 alpha2 201.250
alpha1 = 0.446512 -0.083721 * 95.000
M0 x M1 Span 1 Span 2 Span 3 Totals alpha2 -0.083721 0.390698 201.250
PL -50.000 0.000 0.000
UDL 0.000 -90.000 0.000 alpha1 = 25.57
Total -50.000 -45.000 0.000 -95.000 alpha2 70.67
M1 x M1 Span 1 Span 2 Span 3
Total 1.333 1.000 0.000 2.333
Mid-Span and Support Moments
M2 x M1 Span 1 Span 2 Span 3 Mab 37.2 kNm
Total 0.000 0.500 0.000 0.500 Mb -25.6 kNm
Mbc -3.1 kNm
M0 x M2 Span 1 Span 2 Span 3 Mc -70.7 kNm
PL 0.000 0.000 0.000 Mcd 58.4 kNm
UDL 0.000 -90.000 -156.250
Total 0.000 -45.000 -156.250 -201.250
Total 0.000 0.500 0.000 0.500
Total 0.000 1.000 1.667 2.667

Dr. C. Caprani 95
5.7 Problems
1. For the truss shown, calculate the vertical deflection at C when member EB is
not present and when it is present. Members AB, DC and DE have EA = 60
10
3
kN; members EB, BC and AD have EA = 100 10
3
kN, and; member BD
has EA = 80 10
3
kN. (Ans. 15 mm and 11.28 mm)

2. Autumn 2007 For the truss shown, determine the force in each member and
determine the horizontal deflection of joint B. Take
3
200 10 kN EA= for all
members. (Ans. Choosing BD: 50 2 o = ; 2 mm)

FIG. Q3
A
C
B
D
100 kN

Dr. C. Caprani 96
3. For the frame shown, calculate the bending moment diagram, and verify the
following deflections:
856.2
Bx
EI
o = and
327.6
Cy
EI
o = +. (Ans.
30.1 kNm
B
M = )

4. Summer 2007, Part (b): For the frame shown, draw the bending moment
diagram and determine the horizontal deflection of joint C. Neglect axial
effects in the members and take
3 2
36 10 kNm EI = . (Ans. 75.8 kN
D
H = ;
19.0 mm)
FIG. Q3(b)
A
B
C
2
0

k
N
/
m
100 kN
D

Dr. C. Caprani 97
5. For the beam of Example 1, show that the deflections at the centre of each
span, taking downwards as positive, are:
- Span AB:
41.1
EI
;
- Span BC:
23.9
EI
;
- Span CD:
133.72
EI
.

6. Using virtual work, analyse the following prismatic beam to show that the
support moments are 51.3 kNm
B
M = and 31.1 kNm
C
M = , and draw the
bending moment diagram:

Dr. C. Caprani 98
6. Virtual Work for Self-Stressed Structures
6.1 Background
Introduction
Self-stressed structures are structures that have stresses induced, not only by external
loading, but also by any of the following:
- Temperature change of some of the members (e.g. solar gain);
- Lack of fit of members from fabrication:
o Error in the length of the member;
o Ends not square and so a rotational lack of fit;
- Incorrect support location from imperfect construction;
- Non-rigid (i.e. spring) supports due to imperfect construction.

Since any form of fabrication or construction is never perfect, it is very important for
us to know the effect (in terms of bending moment, shear forces etc.) that such errors,
even when small, can have on the structure.

Here we introduce these sources, and examine their effect on the virtual work
equation. Note that many of these sources of error can exist concurrently. In such
cases we add together the effects from each source.

Dr. C. Caprani 99
Temperature Change
The source of self-stressing in this case is that the temperature change causes a
member to elongate:

( )
T
L L T o A = A

where o is the coefficient of linear thermal expansion (change in length, per unit
length, per degree Celsius), L is the original member length and T A is the
temperature change.

Since temperature changes change the length of a member, the internal virtual work
is affected. Assuming a truss member is being analysed, we now have changes in
length due to force and temperature, so the total change in length of the member is:

( )
PL
e L T
EA
o = + A

Hence the internal virtual work for this member is:

( )
I
W e P
PL
L T P
EA
o o
o o
=
| |
= + A
|
\ .

Dr. C. Caprani 100
Linear Lack of Fit
For a linear lack of fit, the member needs to be artificially elongated or shortened to
fit it into place, thus introducing additional stresses. This is denoted:

L

Considering a truss member subject to external loading, the total change in length
will be the deformation due to loading and the linear lack of fit:

L
PL
e
EA
= +

Hence the internal virtual work for this member is:

I
L
W e P
PL
P
EA
o o
o
=
| |
= +
|
\ .

Dr. C. Caprani 101
Rotational Lack of Fit
A rotational lack of fit, which applies to frames only, occurs when the end of a
member is not square. Thus the member needs to be artificially rotated to get it into
place, as shown below. This is denoted as:

u

Considering a frame member which has a lack of fit,
u
and a virtual moment M o at
the same point, then the internal virtual work done at this point is:

I
W M
u
o o =

This must be added to the other forms of internal virtual work. Not also that the signs
must be carefully chosen so that the virtual moment closes the gap we will see this
more clearly in an example.

Dr. C. Caprani 102
Errors in Support Location
The support can be misplaced horizontally and/or vertically. It is denoted:

S

A misplaced support affects the external movements of a structure, and so contributes
to the external virtual work. Denoting the virtual reaction at the support, in the
direction of the misplacement as R o , then we have:

e S
W R o o =

Dr. C. Caprani 103
Spring Supports
For spring supports we will know the spring constant for the support, denoted:

S
k

Since movements of a support are external, spring support movements affect the
external virtual work. The real displacement
S
A that occurs is:

S S
R k A =

In which R is the real support reaction in the direction of the spring. Further, since R
will be known in terms of the multiplier and virtual reaction, R o , we have:

0
R R R o o = +

Hence:

( )
0
S S
R R k o o A = +

And so the external work done is:

( )
2
0
e S
S S
W R
R
R
R
k k
o o
o
o o
= A
= +

The only unknown here is o which is solved for from the virtual work equation.

Dr. C. Caprani 104
6.2 Trusses
Example 16
Here we take the truss of Example 6 and examine the effects of:
- Member ED was found to be 5 mm too long upon arrival at site;
- Member AB is subject to a temperature increase of +100 C.
For this truss, E = 200 kN/mm
2;
member areas, A = 1000 mm
2
AE and BD where A = 1000\2 mm
2
, and; coefficient of expansion is
6 -1
10 10 C o
= .

The change in length due to the temperature change is:

( )
( )( )( )
6
10 10 2000 100
2 mm
T
L L T o
A = A
= +
=

Vertical Displacement of Joint D
The virtual work equation is now:

Dr. C. Caprani 105
0
1 1 1
0
1 5 2
E I
i i i i
DV i ED AB
i
W
W W
y F e P
P L
y P P P
EA
o
o o
o o
o o o
=
=
=
| |
= + +
|
\ .

In Example 6 we established various values in this equation:

-
0
1
16.5
i
i
P L
P
EA
o
| |
=
|
\ .

-
1
1
ED
P o =
-
1
1
AB
P o = +

Hence we have:

( ) ( )
1 16.5 5 1 2 1
13.5 mm
DV
DV
y
y
= + + +
=

Horizontal Displacement of Joint D:
Similarly, copying values from Example 6, we have:

( ) ( )
0
2 2 2
0
1 5 2
4.5 5 1 2 0
0.5 mm to the right
E I
i i i i
DH i ED AB
i
DH
W
W W
y F e P
P L
y P P P
EA
y
o
o o
o o
o o o
=
=
=
| |
= + +
|
\ .
= + + +
= +

Dr. C. Caprani 106
Example 17
Here we use the truss of Example 1 and examine, separately, the effects of:
- Member AC was found to be 3.6 mm too long upon arrival on site;
- Member BC is subject to a temperature reduction of -50 C;
- Support D is surveyed and found to sit 5 mm too far to the right.
Take EA to be 1010
4
kN for all members and the coefficient of expansion to be
6 -1
10 10 C o
= .

Error in Length:
To find the new multiplier, we include this effect in the virtual work equation:

1 1
0
0
E I
i i i i
i L AC
i
W
W W
y F e P
PL
P P
EA
o
o o
o o
o o
=
=
=
| |
= +
|
\ .

But
0 1
P P P o o = + , hence:
Dr. C. Caprani 107
( )
( )
0 1
1 1
0 1
1 1 1
2
1
0 1
1
0
0
i L AC
i
i i L AC
i i
i i
i i
L AC
i i
P P L
P P
EA
P L P L
P P P
EA EA
P L
P P L
P
EA EA
o o
o o
o
o o o o
o
o
o o
| |
+
| = +
|
\ .
| | | |
= + +
| |
\ . \ .

= + +

As can be seen, the
1
L AC
P o term is simply added to the usual VW equation. From
before, we have the various values of the summations, and so have:

( )( )
4 4
37.8 10 1.62 10
0 3.6 1
EA EA
o

= + + +
4
4
37.8 10 3.6
1.11
1.62 10
EA
o

= = +

Thus member AC is 1.1 kN in tension. Note the change: without the error in fit it was
23.3 kN in tension and so the error in length has reduced the tension by 22.2 kN.

Temperature Change:
The same derivation from the VW equation gives us:

( )
2
1
0 1
1
0
i i
i i
T BC
i i
P L
P P L
L P
EA EA
o
o
o o

= + + A

The change in length due to the temperature change is:

( )
( )( )( )
6
10 10 3000 50 1.5 mm
T
L L T o
A = A
= =

Dr. C. Caprani 108

Noting that the virtual force in member BC is
1
0.6
BC
P o = , we have:

( )( )
4 4
37.8 10 1.62 10
0 1.5 0.6
EA EA
o

= + +

Giving:

4
4
37.8 10 0.9
17.7
1.62 10
EA
o

= = +

Thus member AC is 17.7 kN in tension, a decrease of 5.6 kN from 23.3 kN.

Error in Support Location:
Modifying the VW equation gives:

( )
2
1
0 1
i i
i i
S D
i i
P L
P P L
H
EA EA
o
o
o o

= +

The value of the virtual horizontal reaction is found from the virtual force system to
be 0.6 kN to the right. Hence:

( )( )
4 4
37.8 10 1.62 10
5 0.6
EA EA
o

+ = +
Note the sign on the support displacement: since the real movement is along the same
direction as the virtual force, it does positive virtual work. Solving:

4
4
3 37.8 10
4.8
1.62 10
EA
o
+
= = +

Dr. C. Caprani 109

Thus member AC is 4.8 kN in tension; a reduction of 18.5 kN.

All Effects Together:
In this case, the virtual work equation is:

( )
2
1
0 1
1 1
i i
i i
S D L AC T BC
i i
P L
P P L
H P L P
EA EA
o
o
o o o o

= + + + A

Substituting the various values in gives:

( )( ) ( )( ) ( )( )
4 4
37.8 10 1.62 10
0.5 0.6 3.6 1 1.5 0.6
EA EA
o

+ = + + + +

And solving:

( )
4
4
3.6 0.9 3 37.8 10
22.9
1.62 10
EA
o
+
= =

And so member AC is 22.9 kN in compression. Thus should be the same as the
original force of 23.3 kN plus all the changes induced by the errors:

23.3 22.2 5.6 18.5 23.0 kN =

Exercise:
Find the remaining forces in the truss and compare to the forces without the presence
of self-stresses.
Dr. C. Caprani 110
6.3 Beams
Example 18
Consider the following 2-span beam with central spring support. Determine an
expression for the central support reaction in terms of its spring stiffness.

To do this, we will consider the central support as the redundant:

The virtual work equation, accounting for the spring, is:

Dr. C. Caprani 111
( )
2
1
0 1
1
i
i
S
i i
M
M M
dx dx
EI EI
o
o
o
A = +
} }

The deflection at the central support will be:

B
S
V
k k
o
A = =

Since the reaction at B is 1 o . Noting that the deflection of the spring will be
opposite to the unit load, and using the volume integrals, we have:

( ) ( )
2
4 3
3 4
2 5 2 1
12 2 2 3 2 2
5
24 6
5
6 24
l wl l l
l l
k EI EI
wl l
k EI EI
EI l wl
k
o
o
o
o
o
(
| | (
| | | || |
= +
( | | | |
(
\ . \ .\ .
\ .

= +
| |
=
|
\ .

And finally:

4
3
3
5
24
4
5
24 4
wl
EI
l
k
wl
EI
kl
o =
+
=
+

It is important we understand the implications of this result:

Dr. C. Caprani 112
- For no support present, 0 k = and so 24EI k meaning 0 o and there is no
support reaction (as we might expect);
- For k = we have the perfectly rigid (usual) roller support and so 24 0 EI k
giving us
5
4
wl o - a result we established previously;
- For intermediate relative stiffnesses, the value of the reaction is between these
extremes. Therefore, we as the designer have the ability to choose the reaction
most favourable to us.

A plot of the reaction and relative stiffnesses is:

0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.E-10 1.E-08 1.E-06 1.E-04 1.E-02 1.E+00 1.E+02 1.E+04 1.E+06 1.E+08 1.E+10
Relative Stiffness (EI/kL
3
)
C
e
n
t
r
a
l

R
e
a
c
t
i
o
n

(
o
/
w
L
)

Dr. C. Caprani 113
6.4 Frames
Example 19
The following frame, in addition to its loading, is subject to:
- Support A is located 10 mm too far to the left;
- End C of member BC is
3
1.2 10 rads
out of square, as shown;

- Member CD is 12 mm too short.
Determine the bending moment diagram. Take EI = 36 10
3
kNm
2
for each member.

We will choose the horizontal reaction at A as the redundant. Since we are dealing
with a linear lack of fit in member CD, we need to allow for the virtual work done by
the axial forces in this member and so we solve for the axial force diagrams also.

For the primary structure:

0 0
0: 90 0 90 kN
x D D
F H H = = =

Dr. C. Caprani 114
2
0 0
about 0:
6
24 90 4 90 2 6 0 162 kN
2
D D
M A
V V
=
+ + = =

0 0
0: 24 6 162 0 18 kN
y A A
F V V = + = =

Also, 90 6 540 kNm
C
M = = giving:

For the redundant structure:

Dr. C. Caprani 115
1 1
1 1
1
about 0: 1 2 6 0 kN
3
1 1
0: 0 kN
3 3
D D
y A A
M A V V
F V V
= + = = +
= + = = |

The virtual work equation, accounting for the relevant effects is:

( )
2
1
0 1
1 1 1
0
E I
i
i
S A C L CD
i i
W
W W
M
M M
H dx dx M P
EI EI
u
o
o o
o
o
o o o
=
=
= + + +
} }

We take each term in turn:
(a)
1
S A
H : The applied unit load is in the same direction as the error in the
support location, and keeping all units in metres:

( )
1 3 3
10 10 1 10 10 kN m
S A
H

= + =

(b)
0 1
i
i
M M
dx
EI
o
}
: Using the volume integrals:

( )( ) ( )( )
0 1
1 1 1
4 3 6 540 6 6 540 6
12 3
12420
i
BC CD i
M M
dx
EI EI
EI
o
( (
= + +
`
( (

)
=
}

(c)
( )
2
1
i
i
M
dx
EI
o
}
: Again using the integrals table:

Dr. C. Caprani 116
( )
( )( ) ( ) ( ) ( )
( )( )
2
1
1 1 1 1
4 4 4 4 2 4 6 6 4 2 6 6
3 6
1 1
6 6 6
3
i
AB BC i
CD
M
dx
EI EI EI
EI
o
( (
= + + + +
( (

(
+
(

}

( )
2
1
245.33 i
i
M
dx
EI EI
o
=
}

(d)
1
C
M
u
o : The virtual bending moment at C is 6 kNm with tension on the
inside of the frame. This is in the opposite direction to that needed to close the
lack of fit and so the sign of this term is negative, as shown:

( ) ( )
1 3 3
1.2 10 6 7.2 10 kN m
C
M
u
o

= =

(e)
1
L CD
P o : Since the member is too short, the lack of fit is negative, whilst the
virtual force is in tension and so positive:

( ) ( )
1 3 3
1
12 10 4 10 kN m
3
L CD
P o

| |
= + =
|
\ .

Substituting these values into the equation, along with
3 2
36 10 kNm EI = , gives:
Dr. C. Caprani 117

3 3 3
3 3
12420 245.33
10 10 7.2 10 4 10
36 10 36 10
o

= +

And so we can solve for o :

10 345 6.815 7.2 4
10 345 7.2 4
6.815
53.73
o
o
= +
+ + +
=
=

And so the horizontal reaction at A is:

0 1
0 53.73 1 53.73 kN
A A A
H H H o = + = + =

and similarly for the other reactions. Also, using the superposition equation for
moments,
0 1
M M M o = + , we have:

Dr. C. Caprani 118
6.5 Problems
1. For the truss of Example 1, show that the force in member AC is 22.2 kN in
compression when there is no external loads present, and only the lack of fit of
3.6 mm of member AC.

2. For the truss of Example 1, determine the horizontal deflection of joint C due
to each of the errors separately, and then combined.

3. For the beam of Example 1, show that the stiffness of the spring support that
optimizes the bending moments is
3
89.5EI l , i.e. makes the sagging and
hogging moments equal in magnitude.

4. For the following beam, show that the vertical deflection at C can be given by:

3
2 1
3
Cy
L
P
EI k
o
| |
= +
|
\ .

Dr. C. Caprani 119
5. For the frame of Example 1, show that the horizontal deflection of joint C due
to the applied loads only is 77.5 mm. Find the deflection of joint C due to both
the loads and errors given.

6. For the following frame, the support at D was found to yield horizontally by
0.08 mm/kN. Also, the end of member BC is not square, as shown. Draw the
bending moment diagram and determine the horizontal deflection at D. Take
EI = 10010
3
kNm
2
. (Ans. 5.43 mm to the right)

Dr. C. Caprani 120
7. Past Exam Questions
7.1 Summer 1997

Answers:
(a) 30.0 mm
Cx
o = ;
(b) 12.9 mm
Cx
o = ;
(c) 21.1 mm
BD
= too long.
Dr. C. Caprani 121
7.2 Summer 1998

Answers:
(a) 45.0 mm
Cx
o = ;
(b) 26.2 mm
Cx
o = .
Dr. C. Caprani 122
7.3 Summer 1999

Answers:
(a) 42.5 kN
E
V = |;
(b) 28.6 mm
Cy
o = +.
Dr. C. Caprani 123
7.4 Summer 2000

Answers:
(a) 2.78 kN
A
V = +;
(b) 0.47 mm
C
o = +.
Dr. C. Caprani 124
7.5 Summer 2001

Answers:
(a) 44.4 kN
BD
P = + ;
(b) 300 kN
BD
P = + .

Dr. C. Caprani 125
7.6 Summer 2002

Answers:
(a) 1.35 mm
L
= shorter;
(b) 160 kN
C
V = |.
Dr. C. Caprani 126
7.7 Summer 2004

Answers:
(a) 2.19 mm
Ay
o = +;
(b) 5 mm
L
= shorter.

Dr. C. Caprani 127
7.8 Summer 2007
Using Virtual Work:
(i) For the frame in Fig. Q3(a), determine the horizontal deflection of joint C;
(ii) Draw the bending moment diagram for the frame in Fig. Q3(b);
(iii) Determine the horizontal deflection of joint C for the frame in Fig. Q3(b).

Note:
You may neglect axial effects in the members.
Take
3 2
36 10 kNm EI = for all members.

FIG. Q3(a)
A
B
C
2
0

k
N
/
m
100 kN
D
FIG. Q3(b)
A
B
C
2
0

k
N
/
m
100 kN
D

Answers:
(a) 84.6 mm
Cx
o = ;
(b) 75.8 kN
D
H = ;
(c) 19.02 mm
Cx
o = .

Dr. C. Caprani 128
7.9 Semester 2, 2008
For the frame shown in Fig. Q2, using Virtual Work:
(a) Determine the following:
(i) reactions;
(ii) bending moment diagram.

(b) After construction, it is found that the support at C is not rigid but a spring support of stiffness 2000 3 kN/m k = .
Taking this into account, determine the following:
(i) vertical displacement of joint C;
(ii) bending moment diagram.
(40 marks)

Note:
You may neglect axial effects in the members.
Take
3 2
36 10 kNm EI = for all members.

FIG. Q2
B
A
10 kN/m
C
20 kN

Answers:
(a) 35.63 kN
C
V = |;
(b) 45 mm
Cy
o = +.
Dr. C. Caprani 129
7.10 Semester 2, 2009
QUESTION 2

Using Virtual Work, for the frame shown in Fig. Q2, do the following:

(i) Determine the reactions;

(ii) Draw the bending moment diagram;

(iii) Draw the deflected shape of the structure;

(iv) Determine the horizontal deflection of joint D.

(25 marks)

Note:
- Take
2
200 kN/mm E = and
7 4
45 10 mm I = for all members.

A
150 kN

FIG. Q2
B C
D

Answers:
(a) 225 kN
B
V = |;
(b) 81.25 mm
Dx
o = .
Dr. C. Caprani 130
7.11 Semester 2, 2009
QUESTION 3

For the structure shown in Fig. Q2, due to an error on site, it is required to make the moment at support A zero. This is
to be done by the introduction of a spring support at B, in place of the roller support, as shown in Fig. Q3. Using Virtual
Work:

(i) Determine the spring stiffness required so that the moment at A is zero.

(ii) Drawn the revised bending moment diagram;

(iii) Draw the revised deflected shape of the structure;

(iv) Determine the vertical deflection of the structure at B.

(25 marks)

Note:
- Take
2
200 kN/mm E = and
7 4
45 10 mm I = for all members.
- You may use any relevant results from your workings for Q2, but in doing so acknowledge their source.

A
150 kN

FIG. Q3
B C
D
k

Answers:
(a)
3
20 10 kN/m k = ;
(b) 7.5 mm
By
o = +.
Dr. C. Caprani 131
7.12 Semester 2, 2010
QUESTION 2

Using Virtual Work, for the continuous beam shown in Fig. Q2, do the following:

(i) Draw the bending moment diagram, noting all important values;

(ii) Draw the deflected shape of the structure;

(iii) Determine the deflection at E.

(25 marks)

Note:
- Take
2
200 kN/mm E = and
7 4
20 10 mm I = for all members;
- The following expression may be of assistance:

0 1 1 1 2 1
1
0 2 1 2 2 2
2
0
0
M M M M M M
EI EI EI
M M M M M M
EI EI EI
o o o o o
o
o
o o o o o
(
(

( + =
` ` `
(
) )

(

)
} } }
} } }

A
B

FIG. Q2
D
C
60 kN/m
100 kN
40 kN/m
E

Answers:
(a) 174.5 kNm; 97.6 kNm
B C
M M = = ;
(b) 16.9 mm
Ey
o = +.
Dr. C. Caprani 132
7.13 Semester 2, 2010
QUESTION 3

Using Virtual Work, for the truss shown in Fig. Q3, do the following:

(i) Determine the member forces;

(ii) Draw the member forces using conventional arrow notation;

(v) Draw the deflected shape of the structure;

(iii) Determine the vertical deflection of joint E.

(25 marks)

Note:
- Take
2
200 kN/mm E = ;
-
2
1500 mm A= for members BC, BF, CE, and EF;
-
2
1500 2 mm A= for members BE and CF;
-
2
2000 mm A= for members AF and DE;
-
2
2500 mm A= for members AB and CD.

FIG. Q3

60 kN
F
B C
D
E
60 kN
A

Answers:
(a) 15 2 kN
BE
F = compression;
(b) 4.25 mm
Ey
o = +.
Dr. C. Caprani 133
8. References
- Charlton, T.M., Analysis of Statically Indeterminate Frameworks, Longmans,
1961.
- Charlton, T.M., Energy Principles in Theory of Structures, Oxford University
Press, 1973.
- Coates, R.C., Coutie, M.G. and Kong, F.K., Structural Analysis, Chapman and
Hall, 1987.
- Davies, G.A.O., Virtual Work in Structural Analysis, John Wiley & Sons, 1982.
- Dym, C.L., Structural Modeling and Analysis, Cambridge University Press, 2005.
- Gregory, M.S., Introduction to Extremum Principles, Butterworths, London, 1969.
- Guarracino, F. and Walker, A., Energy Methods in Structural Mechanics, Thomas
Telford, 1999.
- Heyman, J., Beams and Framed Structures, 2nd Edn., Pergamon Press, 1974.
- Kong, F.K., Prentis, J.M. and Charlton, T.M., Principle of virtual work for a
general deformable body a simple proof, The Structural Engineer, Vol. 61A,
No. 6, 1983.
- Neal, B.G., Structural Theorems and their Applications, Pergamon Press, 1964.
- Thompson, F. and Haywood, G.G., Structural Analysis Using Virtual Work,
Chapman and Hall, 1986.
- Rajasekaran, S. and Sankarasubramanian, G., Computational Structural
Mechanics, Prentis Hall of India, New Delhi, 2001.
- Richards, T.H., Energy methods in Stress Analysis, Ellis Horwood, 1977.
- Wallerstein, D.V. , A Variational Approach to Structural Analysis, Wiley, 2001.
Dr. C. Caprani 134
9. Appendix Volume Integrals

l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
i
k
l
k
l
k
l
k
l
k

l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
i
k
l
k
l
k
l
k
l
k

l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
i
k
l
k
l
k
l
k
l
k

l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
i
k
l
k
l
k
l
k
l
k

l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

1
3
jkl
1
6
jkl ( )
1 2
1
2
6
j j kl +
1
2
jkl
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

1
6
jkl
1
3
jkl ( )
1 2
1
2
6
j j kl +
1
2
jkl
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

( )
1 2
1
2
6
j k k l + ( )
1 2
1
2
6
j k k l +
( )
( )
1 1 2
2 1 2
1
2
6
2
j k k
j k k l
+ +
+ (

( )
1 2
1
2
j k k l +
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

1
2
jkl
1
2
jkl ( )
1 2
1
2
j j kl + jkl
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

( )
1
6
jk l a + ( )
1
6
jk l b +
( )
( )
1
2
1
6
j l b
j l a k
+ +
+ (

1
2
jkl
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

5
12
jkl
1
4
jkl ( )
1 2
1
3 5
12
j j kl +
2
3
jkl
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

1
4
jkl
5
12
jkl ( )
1 2
1
5 3
12
j j kl +
2
3
jkl
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

1
4
jkl
1
12
jkl ( )
1 2
1
3
12
j j kl +
1
3
jkl
l
j
l
j
l
j j
1 2
l
j
a b
k
l
k
l
k
l
k k
1 2
l
k
l
k
l
k
l
k
l
k

1
12
jkl
1
4
jkl ( )
1 2
1
3
12
j j kl +
1
3
jkl
k
l

1
3
jkl
1
3
jkl ( )
1 2
1
3
j j kl +
2
3
jkl

Virtual Work 1011

Uploaded by

Copyright:

Available Formats

Virtual Work 1011

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Virtual Work 1011

Uploaded by

Copyright:

Available Formats

What is the principle of virtual work and how is it used in structural analysis?

How are virtual displacements and forces applied to rigid and deformable bodies?

Structural Analysis III

. Thus equilibrium occurs when H is an

whilst the second is the same as

. Thus equilibrium equations occur naturally within the virtual

out of square, as shown;

You might also like