Module

4
Analysis of Statically
Indeterminate
Structures by the Direct
Stiffness Method
Version 2 CE IIT, Kharagpur
 Lesson
25
The Direct Stiffness
Method: Truss Analysis
(Continued)
Version 2 CE IIT, Kharagpur
Instructional Objectives
After reading this chapter the student will be able to
1. Transform member stiffness matrix from local to global co-ordinate system.
2. Assemble member stiffness matrices to obtain the global stiffness matrix.
3. Analyse plane truss by the direct stiffness matrix.
4. Analyse plane truss supported on inclined roller supports.

25.1 Introduction
In the previous lesson, the direct stiffness method as applied to trusses was
discussed. The transformation of force and displacement from local co-ordinate
system to global co-ordinate system were accomplished by single transformation
matrix. Also assembly of the member stiffness matrices was discussed. In this
lesson few plane trusses are analysed using the direct stiffness method. Also the
problem of inclined support will be discussed.

Example 25.1
Analyse the truss shown in Fig. 25.1a and evaluate reactions. Assume EA to be
constant for all the members.

Version 2 CE IIT, Kharagpur

The numbering of joints and members are shown in Fig. 25.1b. Also, the possible
displacements (degrees of freedom) at each node are indicated. Here lower
numbers are used to indicate unconstrained degrees of freedom and higher
numbers are used for constrained degrees of freedom. Thus displacements 6,7
and 8 are zero due to boundary conditions.

First write down stiffness matrix of each member in global co-ordinate system
and assemble them to obtain global stiffness matrix.

Element 1: = 60, L = 4.619 m. Nodal points 4-1

0.25 0.433 0.25 0.433

0.433 0.433 0.75
[ ]
k =
1 EA 0.75
4.619 0.25 0.433 0.25 0.433
(1)

0.433 0.75 0.433 0.75

Element 2: = 90, L = 4.00 m. Nodal points 2-1

Version 2 CE IIT, Kharagpur

 0 0 0 0
0 1
[k ] 2
=
EA 0 1
4.0 0 0 0 0
(2)

0 1 0 1

Element 3: = 120, L = 4.619 m. Nodal points 3-1

0.25 0.433 0.25 0.433

0.433 0.75
[k ]
3
=
EA
4.619 0.25
0.75
0.433
0.433
0.25 0.433
(3)

0.433 0.75 0.433 0.75

Element 4: = 0, L = 2.31 m. Nodal points 4-2

1 0 1 0
0
[k ]
4
=
EA 0
2.31 1
0 0
0 1 0
(4)

0 0 0 0

Element 5: = 0, L = 2.31 m. Nodal points 2-3

1 0 1 0
0
[k ]
5
=
EA 0
2.31 1
0 0
0 1 0
(5)

0 0 0 0

The assembled global stiffness matrix of the truss is of the order 8 8 . Now
assemble the global stiffness matrix. Note that the element k111 of the member
stiffness matrix of truss member 1 goes to location (7,7 ) of global stiffness matrix.
On the member stiffness matrix the corresponding global degrees of freedom are
indicated to facilitate assembling. Thus,

Version 2 CE IIT, Kharagpur

 1 2 3 4 5 6 7 8
0.108 0 0 0 0.054 0.094 0.054 0.094 1
0 0.575 0 0.25 0.094 0.162 0.094 0.162 2

0 0 0.866 0 0.433 0 0.433 0 3

0 0.25 0 0.25 0 0 0 0 4
[K ] = EA
0.054 0.094 0.433 0 0.487 0.094 0 0 5

0.094 0.162 0 0 0.094 0.162 0 0 6
0.054 0.094 0.433 0 0 0 0.487 0.0934 7

0.094 0.162 0 0 0 0 0.0934 0.162 8

(6)

Writing the load-displacement relation for the truss, yields

p1 0.108 0 0 0 0.054 0.094 0.054 0.094 u1

p 0 0.575 0 0.25 0.094 0.162 0.094 0.162 u 2
2
p3 0 0 0.866 0 0.433 0 0.433 0 u 3

p4 0 0.25 0 0.25 0 0 0 0 u 4
= EA
p5 0.054 0.094 0.433 0 0.487 0.094 0 0 u 5

p6 0.094 0.162 0 0 0.094 0.162 0 0 u 6
0.054 0.094 0.433
p7 0 0 0 0.487 0.0934 u 7

p 0.094 0.162 0.0934 0.162 u 8
8 0 0 0 0
(7)

The displacements u1 to u 5 are unknown. The displacements u 6 = u 7 = u8 = 0 .

Also p1 = p 2 = p3 = p5 = 0 . But p4 = 10 kN .

0 0.108 0 0 0 0.054 u1
0

0
0.575 0 0.25 0.094 u 2

0 = EA 0 0 0.866 0 0.433 u 3 (8)
10
0 0.25 0 0.25 0 u 4

0 0.054 0.094 0.433 0 0.487 u 5

Solving which, the unknown displacements are evaluated. Thus,

Version 2 CE IIT, Kharagpur

 6.668 34.64 6.668 74.642 13.334
u1 = ; u2 = ; u3 = ; u4 = ; u5 = (9)
AE AE AE AE AE
Now reactions are evaluated from equation,

6.668
34.64
p6 0.094 0.162 0 0 0.094
1
p 7 = EA 0.054 0.094 0.433 0 0 6.668 (10)
p EA
8 0.094 0.162 0 0 0 74.642

13.334

Thus,

p6 = 5.00 kN ; p7 = 0 ; p8 = 5.00 kN . (11)

Now calculate individual member forces.

Member 1: l = 0.50 ; m = 0.866 ; L = 4.619 m .

u 7
u

{p'1 } = AE
[l m l m ] 8
4.619 u1
u 2

AE 1 6.667
{ p '1} = [ 0.5 0.866] = 5.77 kN (12)
4.619 AE 34.64

Member 2: l = 0 ; m = 1.0 ; L = 4.0 m .

u 3
u

{p'1 } = AE
[l m l m ] 4
4 .0 u1
u 2

AE 1 74.642
{ p '1} = [1 1] = 10.0 kN (13)
4.619 AE 34.64

Member 3: l = 0.50 ; m = 0.866 ; L = 4.619 m .

Version 2 CE IIT, Kharagpur

 u 5
u

{p'1 } = AE
[l m l ]
m 6
4.619 u1
u 2

13.334
AE 1
{ p '1} = [ 0.5 0.5 0.866] 6.667 = 5.77 kN (14)
4.619 AE
34.64

Member 4: l = 1.0 ; m = 0 ; L = 2.31.0 m .

u 7
u

{p'1 } = AE
[l m l m ] 8
2.31 u 3
u 4

AE 1 0
{ p '1} = [1 1] = 2.88 kN (15)
2.31 AE 6.667

Member 5: l = 1.0 ; m = 0 ; L = 2.31.0 m .

u 3
u

{p'1 } = AE
[l m l m ] 4
2.31 u 5
u 6

AE 1 6.667
{ p '1} = [1 1] = 2.88 kN (16)
2.31 AE 13.334

Example 25.2
Determine the forces in the truss shown in Fig. 25.2a by the direct stiffness
method. Assume that all members have the same axial rigidity.

Version 2 CE IIT, Kharagpur

Version 2 CE IIT, Kharagpur
The joint and member numbers are indicated in Fig. 25.2b. The possible degree
of freedom are also shown in Fig. 25.2b. In the given problem u1 ,u 2 and u 3
represent unconstrained degrees of freedom and u 4 = u 5 = u 6 = u 7 = u8 = 0 due
to boundary condition. First let us generate stiffness matrix for each of the six
members in global co-ordinate system.

Element 1: = 0, L = 5.00 m. Nodal points 2-1

3 4 1 2
1 0 1 0 3
0
[ ]
k1 =
EA 0
5.0 1
0 0
0 1 0
4
1
(1)

0 0 0 0 2

Element 2: = 90, L = 5.00 m. Nodal points 4-1

Version 2 CE IIT, Kharagpur

 7 8 1 2
0 0 0 0 7
0 1 8
[k ]
2
=
EA 0 1
5.0 0 0 0 0 1
(2)

0 1 0 1 2

Element 3: = 0, L = 5.00 m. Nodal points 3-4

5 6 7 8
1 0 1 0 5
0
[k ]
3
=
EA 0
5.0 1
0 0
0 1 0
6
7
(3)

0 0 0 0 8

Element 4: = 90, L = 5.00 m. Nodal points 3-2

5 6 3 4
0 0 0 0 5
0 1 6
[k ]
4
=
EA 0 1
5.0 0 0 0 0 3
(4)

0 1 0 1 4

Element 5: = 45, L = 7.07 m. Nodal points 3-1

5 6 1 2
0.5 0.5 0.5 0.5 5
0.5 0.5 0.5 0.5 6
[k ]
5
=
EA
7.07 0.5 0.5 0.5 0.5 1
(5)

0.5 0.5 0.5 0.5 2

Element 6: = 135, L = 7.07 m. Nodal points 4-2

Version 2 CE IIT, Kharagpur

 7 8 3 4
0.5 0.5 0.5 0.5 7
0.5 0.5 8
EA 0.5 0.5
[k ]
6
=
7.07 0.5 0.5 0.5 0.5 3
(6)

0.5 0.5 0.5 0.5 4

There are eight possible global degrees of freedom for the truss shown in the
figure. Hence the global stiffness matrix is of the order ( 8 8 ). On the member
stiffness matrix, the corresponding global degrees of freedom are indicated to
facilitate assembly. Thus the global stiffness matrix is,

0.271 0.071 0.2 0 0.071 0.071 0 0

0.071 0.271 0 0 0.071 0.071 0 0.20

0.20 0 0.271 0.071 0 0 0.071 0.071

0 0 0.071 0.271 0.071 0.20 0.071 0.071
[K ] = AE (7)
0.071 0.071 0 0 0.271 0.071 0.20 0

0.071 0.071 0 0.20 0.071 0.271 0 0
0 0 0.071 0.071 0.20 0 0.271 0.071

0 0.20 0.071 0.071 0 0 0.071 0.271

The force-displacement relation for the truss is,

p1 0.271 0.071 0.2 0 0.071 0.071 0 0 u1

p 0.071 0.271 0 0 0.071 0.071 0 0.20 u 2
2
p3 0.20 0 0.271 0.071 0 0 0.071 0.071 u 3

p4 0 0 0.071 0.271 0.071 0.20 0.071 0.071 u 4
= EA
p5 0.071 0.071 0 0 0.271 0.071 0.20 0 u 5

p6 0.071 0.071 0 0.20 0.071 0.271 0 0 u 6
0
p7 0 0.071 0.071 0.20 0 0.271 0.071 u 7

p 0.20 0.071 0.071 0.071 0.271 u 8
8 0 0 0
(8)

u1 ,u 2 and u 3 are unknowns.

The displacements
Here, p1 = 5 kN ; p2 = 10 ; p3 = 0 and u 4 = u 5 = u 6 = u 7 = u 8 = 0 .

Version 2 CE IIT, Kharagpur

 5 0.271 0.071 0.2 0 0.071 0.071 0 0 u1
10 0.071 0.271 0 0 0.071 0.071 0 0.20 u 2

0 0.20 0 0.271 0.071 0 0 0.071 0.071 u 3

p4 0 0 0.071 0.271 0.071 0.20 0.071 0.071 0
= EA
0.071 0.071 0.071 0.20 0 0
p5
0 0 0.271

p6 0.071 0.071 0 0.20 0.071 0.271 0 0 0
0
p7 0 0.071 0.071 0.20 0 0.271 0.071 0

p 0.20 0.071 0.071 0.071 0.271 0
8 0 0 0
(9)

Thus,

5 0.271 0.071 0.20 u1

10 = 0.071 0.271 0 u 2 (10)
0 0.20 0.271 u 3
0

Solving which, yields

72.855 55.97 53.825

u1 = ; u2 = ; u3 =
AE AE AE

Now reactions are evaluated from the equation,

p4 0 0 0.071
p 0.071 0.071 0 u1
5 (11)
p 6 = 0.071 0.071 0 u 2
p 0 0

0.071 u 3
7
p8 0 0.20 0.071

p4 = 3.80 kN ; p5 = 1.19 kN ; p6 = 1.19 kN ; p7 = 3.8 0 kN ; p8 = 15.00 kN

In the next step evaluate forces in members.

Element 1: = 0, L = 5.00 m. Nodal points 2-1

Version 2 CE IIT, Kharagpur

 u 3
u

{p'1 } = AE
[l m l ]
m 4
5.0 u1
u 2

AE 1 53.825
{ p '1} = [1 1] = 3.80 kN (12)
5.0 AE 72.855

Element 2: = 90, L = 5.00 m. Nodal points 4-1

u 7
u

{p'1 } = AE
[l m l m ] 8
5 u1
u 2

AE 1 0
{ p '1} = [1 1] = 11.19kN (13)
5 AE 55.97

Element 3: = 0, L = 5.00 m. Nodal points 3-4

{p'1 } = AE
[0 ]{0} = 0 (14)
5

Element 4: = 90, L = 5.00 m. Nodal points 3-2

u 5
u

{p'1 } = AE
[l m l m ] 6
5 u 3
u 4
AE 1
{ p '1} = [ 0] {53.825} = 0 (15)
5 AE

Element 5: = 45, L = 7.07 m. Nodal points 3-1

u 5
u

{p'1 } = AE
[l m l m ] 6
7.07 u1
u 2

Version 2 CE IIT, Kharagpur

 AE 1 72.855
{ p '1} = [ 0.707 0.707] = 1.688 kN (16)
7.07 AE 55.97

Element 6: = 135, L = 7.07 m. Nodal points 4-2

u 7
u

{p'1 } = AE
[l m l m ] 8
7.07 u 3
u 4

AE 1
{ p '1} = [ 0.707] {53.825} = 5.38 kN (17)
7.07 AE

25.2 Inclined supports

Sometimes the truss is supported on a roller placed on an oblique plane (vide
Fig. 25.3a). At a roller support, the displacement perpendicular to roller support is
zero. i.e. displacement along y" is zero in the present case.

Version 2 CE IIT, Kharagpur

If the stiffness matrix of the entire truss is formulated in global co-ordinate system
then the displacements along y are not zero at the oblique support. So, a special
procedure has to be adopted for incorporating the inclined support in the analysis
of truss just described. One way to handle inclined support is to replace the
inclined support by a member having large cross sectional area as shown in Fig.
25.3b but having the length comparable with other members meeting at that joint.
The inclined member is so placed that its centroidal axis is perpendicular to the
inclined plane. Since the area of cross section of this new member is very high, it
does not allow any displacement along its centroidal axis of the joint A . Another
method of incorporating inclined support in the analysis is to suitably modify the
member stiffness matrix of all the members meeting at the inclined support.

Version 2 CE IIT, Kharagpur

Consider a truss member as shown in Fig. 25.4. The nodes are numbered as 1
and 2. At 2, it is connected to a inclined support. Let x' y ' be the local co-ordinate
axes of the member. At node 1, the global co-ordinate system xy is also shown.
At node 2, consider nodal co-ordinate system as x" y" , where y" is perpendicular
to oblique support. Let u '1 and u ' 2 be the displacements of nodes 1 and 2 in the
local co-ordinate system. Let u1 , v1 be the nodal displacements of node 1 in
global co-ordinate system xy . Let u "2 , v "2 be the nodal displacements along x" -
and y" - are in the local co-ordinate system x" y" at node 2. Then from Fig. 25.4,

u '1 = u1 cos x + v1 sin x

u ' 2 = u"2 cos x" + v"2 sin x" (25.1)

This may be written as

u1
u '1 cos x sin x 0 0 v1
=
u ' 2 0 0 cos x" sin x" u"2 (25.2)
v"2
Denoting l = cos x ; m = sin x ; l" = cos x" ; m" = sin x"

u1

u '1 l m 0 0 v1
= (25.3a)
u ' 2 0 0 l" m" u"2
v"2

or {u'} = [T ']{u}
where [T '] is the displacement transformation matrix.

Version 2 CE IIT, Kharagpur

Similarly referring to Fig. 25.5, the force p '1 has components along x and
y axes. Hence

p1 = p '1 cos x (25.4a)

p 2 = p'1 sin x (25.4b)

Version 2 CE IIT, Kharagpur

Similarly, at node 2, the force p ' 2 has components along x" and y" axes.
p "3 = p '2 cos x (25.5a)

p "4 = p '2 sin x (25.5b)

The relation between forces in the global and local co-ordinate system may be
written as,

p1 cos x 0
p sin 0 p '1
2
=
x

p "3 0 cos x p '2 (25.6)
p "4 0
sin x

{p} = [T ']T {p '} (25.7)

Using displacement and force transformation matrices, the stiffness matrix for
member having inclined support is obtained.

[k ] = [T ']T [k '][T ']

l 0
m 0 AE 1 1 l m 0 0
[k ] =
0 l" L 1 1 0 0 l" m" (25.8)

0 m"

Simplifying,

l2 lm ll" lm"

m2 ml" mm"
[k ] =
EA lm
L ll" ml" l"2 l" m" (25.9)

lm" mm" l" m" m"2

If we use this stiffness matrix, then it is easy to incorporate the condition of zero
displacement perpendicular to the inclined support in the stiffness matrix. This is
shown by a simple example.

Version 2 CE IIT, Kharagpur

Example 25.3
Analyse the truss shown in Fig. 25.6a by stiffness method. Assume axial rigidity
EA to be constant for all members.

Version 2 CE IIT, Kharagpur

The nodes and members are numbered in Fig. 25.6b. The global co-ordinate
axes are shown at node 3. At node 2, roller is supported on inclined support.
Hence it is required to use nodal co-ordinates x" y" at node 2 so that u 4 could be
set to zero. All the possible displacement degrees of freedom are also shown in
the figure. In the first step calculate member stiffness matrix.

Member 1: x = 143.13, x" = 6.87, L = 5.00 m. Nodal points 1-2

l = 0.80 ; m = 0.6 ; l" = 0.993 ; m" = 0.12 .

1 2 3 4
0.64 0.48 0.794 0.096 1
0.48 0.596 0.072
[ ]
k1 =
EA 0.36
5.0 0.794 0.596 0.986 0.119
2
3
(1)

0.096 0.072 0.119 0.014 4

Version 2 CE IIT, Kharagpur

Member 2: x = 0, x" = 30, L = 4.00 m. Nodal points 2-3
l = 1 ; m = 0 ; l" = 0.866 ; m" = 0.50 .

Version 2 CE IIT, Kharagpur

 5 6 3 4
0.64 0.48 0.794 0.096 5
0.48 0.596 0.072
[k ]
2
=
EA 0.36
4.0 0.794 0.596 0.986 0.119
6
3
(2)

0.096 0.072 0.119 0.014 4

Member 3: x = 90, L = 3.00 m. , l = 0 ; m = 1 Nodal points 3-1

.
5 6 1 2
0 0 0 0 5
0 1 6
[k ]
3
=
EA 0 1
3.0 0 0 0 0 1
(3)

0 1 0 1 2

For the present problem, the global stiffness matrix is of the order (6 6 ) . The
global stiffness matrix for the entire truss is.

1 2 3 4 5 6
0.128 0.096 0.159 0.019 0 0 1
0.096 0.405 0.119 0.014 0 0.333 2

0.159 0.119 0.385 0.132 0.217 0 3
[k ] = EA
(4)
0.019 .0.014 0.132 0.065 0.125 0 4
0 0 0.217 0.125 0.25 0 5

0 0.333 0 0 0 0.333 6

Writing load-displacement equation for the truss for unconstrained degrees of

freedom,

5 0.128 0.096 0.159 u1

5 = 0.096 0.405 0.119 u 2 (5)
0 0.159 0.119 0.385 u
3

Solving ,

77.408 3.728 33.12

u1 = ; u2 = ; u3 = (6)
AE AE AE

Version 2 CE IIT, Kharagpur

Now reactions are evaluated from the equation

p4 0.019 .0.014 0.132 77.40

1
p 5 = AE 0 0 0.217 3.728 (6)
p AE
6 0 0.333 0 33.12

p4 = 2.85 kN ; p5 = 7.19 kN ; p6 = 1.24 kN

Summary
Sometimes the truss is supported on a roller placed on an oblique plane. In such
situations, the direct stiffness method as discussed in the previous lesson needs
to be properly modified to make the displacement perpendicular to the roller
support as zero. In the present approach, the inclined support is handled in the
analysis by suitably modifying the member stiffness matrices of all members
meeting at the inclined support. A few problems are solved to illustrate the
procedure.

Version 2 CE IIT, Kharagpur