Introduction to The Principles of Virtual Work

CEE 421L. Matrix Structural Analysis

Department of Civil and Environmental Engineering
Duke University
Henri P. Gavin
Fall, 2014

Forces Fi , Displacements Di and D

Consider a linear elastic beam with N (real) forces, Fi , applied at points 1,2,
..., N .

F1

F2

F3

FN
x

(x,y) (x,y)
11111
00000
00000
11111
00000
11111
00000
11111

x=0

D1

D2

D*

D3

DN

11111
00000
00000
11111
00000
11111
00000
11111

x=L

All of these forces together produce deflections D1 , D2 , ... DN at the corresponding loading points. There are also deflections at other points on the
beam, for example the deflection D at point *. The sum of these forces produce stresses and strains  throughout the volume of the beam. With this
system of loading the external work and internal strain energy are:
W =

N
1X
1
Fi Di = {F }T {D}
2 i=1
2

(1)

and

1Z
U=
{}T {} dV
2 V
where the vectors of stresses {} and strains {} are written
{}T =

{}T =

xx yy zz
xx yy zz

(2)

xy xz yz

(3)

xy xz yz

(4)

CEE 421L. Matrix Structural Analysis Duke University Fall 2014, H.P. Gavin

i and D

Force F , Displacements D
Next, consider the same beam with a single point load, F , (This force could
be at one of the loading points, (1, , N ) but it doesnt have to be.)

F*
x

(x,y) (x,y)
11111
00000
00000
11111
00000
11111
00000
11111

x=0

D1

11111
00000
00000
11111
00000
11111
00000
11111

D2

D3

DN x=L

at point , as well as
The force F by itself produces a collocated deflection D
1, D
2 , ..., D
N , at the other locations on the beam.
deflections D
The force F by itself creates stresses
and strains  throughout the beam.
is
The external work of the force F passing through the displacement D
1
W = F D
(5)
,
2
and the total internal work associated with this single point load is
1Z
{
}T {} dV
(6)
U=
2 V
i + Di and D
+ D
Forces F + Fi , Displacements D
Finally, consider what happens when we apply F first, then apply all the other
(real) loads F1 , F2 , ..., FN , while holding F constant.

F1

F2

F3

FN

11111
00000
00000
11111
00000
11111
00000
11111

11111
00000
00000
11111
00000
11111
00000
11111

x=0 D+D
1 1
D+D
2

D+D x=L
N

D+D
3 3

D+D
*

*
CC BY-NC-ND H.P. Gavin

Introduction to The Principle Of Virtual Work

For the force F :

F 11111111111111111111111111111111111
00000000000000000000000000000000000
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
00000000000000000000000000000000000000
* 11111111111111111111111111111111111
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000

D
*

D +D
*
*

and the external work for force F is

1

W = F D
(7)
+ F D .
2
Note that F is held constant as the other loads, Fi , and deflections, Di , increase.
For the forces at the other locations, Fi , i = 1, . . . , N :

Fi

111111111111111111111111111111111111111
000000000000000000000000000000000000000
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111
000000000000000000000000000000000000000
111111111111111111111111111111111111111

Di

Di + D i
N
1X
W =
Fi Di
2 i=1

(8)

N
1
1X

+
F
D
+
Fi Di .
W = F D

2
2 i=1

(9)

The total external work is

CC BY-NC-ND H.P. Gavin

CEE 421L. Matrix Structural Analysis Duke University Fall 2014, H.P. Gavin

Now, lets find the internal work for the combination of the load F with all the
other loads Fi , i = 1, . . . , N . As before, we will say that the stresses and strains
caused by the force F are
and  and that the stresses and strains caused by
the set of forces Fi , i = 1, . . . , N are and . First we will look at the stresses
due to the force F . The stress
increases linearly until the strain  is attained.
After the force F is applied the stress
remains constant as the strains from
the forces Fi , i = 1, . . . , N are applied.

11111111111111111111111111111111111
00000000000000000000000000000000000
00000000000000000000000000000000000000
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000
00000000000000000000000000000000000000
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000
11111111111111111111111111111111111
00000000000000000000000000000000000000
0000000000000000000000000000000000011111111111111111111111111111111111111
11111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111

Next we will look at the stresses due to the forces F1 , F2 , ..., FN . As the load
F is applied, (i.e., before the forces Fi are applied), the stresses are all zero,
but the strains increase from 0 to . Once the force F has been applied, the
forces Fi are applied, and the strains increase linearly with the stress.

11111111111111111111111111111111111111
00000000000000000000000000000000000000
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111
00000000000000000000000000000000000000
11111111111111111111111111111111111111

The total internal work due to the combined actions of all the loads, applied
sequentially, is
U=

Z
1Z
1Z
{
}T {} dV + {
}T {} dV +
{}T {} dV.
V
V
V
2
2

(10)

CC BY-NC-ND H.P. Gavin

Introduction to The Principle Of Virtual Work

Equating the external and internal work (equations (9) and (10)), noting from
page 1 that
1Z
1
F D =
{
}T {} dV,
(11)
2
2 V
and from page 2 that
N
1X
1Z
Fi D i =
{}T {} dV,
2 i=1
2 V
we obtain the principle of virtual work:

F D =

Z
V

(12)

{
}T {} dV.

(13)

This expression holds for beams, bars, trusses, frames, plates, shells, bricks, etc.
etc. It is customary to call the left hand side of this expression the external
virtual work
= F D ,
W
(14)
and the right hand side the internal virtual work.
U =

Z
V

{
}T {} dV.

(15)

Recall the definitions of the various terms in the principle of virtual work (13).
D is the real deflection at some point in the structure (point *) caused by
the real forces F1 , F2 , ..., FN . Often the displacement D is the unknown
response to be found.
F is a virtual force in the direction and in the location of D .
Often we set this force equal to 1 unit of force.
 are the real strains associated with the system real forces Fi .
These strains are found by analyzing the structure with loads Fi .

are the virtual stresses caused by the virtual force F .

These stresses are found by analyzing the structure with load F .
In words, the external virtual work of a virtual force (F ) moving through a real
displacement (D ) equals the product of the virtual stresses ({
}) corresponding
to the virtual force (F ) with the real strains ({}) corresponding to the real
system of forces (Fi ), integrated over the volume of the solid.
CC BY-NC-ND H.P. Gavin

CEE 421L. Matrix Structural Analysis Duke University Fall 2014, H.P. Gavin

Equation (13) is sometimes called the principle of virtual forces, since the forces
and stresses correspond to the virtual loading system and the displacements
and strains correspond to the real system of loading. In a completely analogous
way we can write the the principle of virtual displacements:
N
X
i=1

i =
Fi D

Z
V

{}T {} dV.

(16)

in which,
Fi are a system of real forces at a set of coordinates (i = 1, , N ).
i are the virtual deflections collocated with the system of forces Fi , arisD
ing from some arbitrary system of virtual forces. Note that since these
virtual deflections are found from a system of virtual forces, the virtual
displacements are consistent with all the constraints (reaction points) on
the structure. In fact, the only requirement for the virtual displacements
is that they adhere to the constraints of the structure.
i.
 are the virtual strains associated with the virtual displacements D
are the real stresses caused by the system of real forces Fi .
In words, the external virtual work of a real forces (Fi ) moving through col i ) equals the product of real stresses ({})
located virtual displacements (D
associated with forces Fi and the virtual strains ({}) associated with virtual
i , integrated over the volume of the solid.
displacements D
The principle of virtual work applies to linear elastic and non-linear elastic
structures also.

CC BY-NC-ND H.P. Gavin

Introduction to The Principle Of Virtual Work

There are many ways in which the principle of virtual work is applied to problems in many fields of engineering and applied mathematics, including solid
mechanics, fluid mechanics, and electro-statics, to name a few. The principle
of virtual work is fundamental to the finite element method, which is used to
solve problems described by systems of partial differential equations in many
disciplines. Here is a typical example of how we can apply the principle of
virtual work to find the deflections at some point in an elastic solid. Consider
a beam carrying some loads
The problem is to find the displacement D at some point (*) on this structure.
Note that the principle of real work can not be applied to this problem.
The figure above shows the real displacements of the structure. This system is
the one we are really interested in and it has internal strains that are really
there.
To apply the principle of virtual work to this problem, we remove all of the externally applied loads in the figure above and apply a unit force in the direction
and location of the unknown displacement D . This unit virtual force (F ) will
cause bending and shear in this particular structure, which will have associated
real stresses. Knowing how to:
relate the real external loads to internal bending moments M (x) and shear
forces V (x);
Relate these internal moments and shears to strains, 
(x) and
Relate the unit virtual load to internal virtual bending moments M
shear forces V (x); and
relate those internal virtual moments and shears to virtual stresses,
,
equation (15) becomes
Z

1 D = {
}T {} dV.
(17)
V

The relationships between internal moments, shears, torques, and axial loads to
stresses and strains may be simplified by considering those cases individually.
Doing so simplifies the integral in equation (17) for these special cases.

CC BY-NC-ND H.P. Gavin