THE UNIVERSITY OF GUYANA

FACULTY OF ENGINEERING AND TECHNOLOGY

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING

Course Code: CIV 4115

Matrix Stiffness Analysis _Rigid Frames
Example 3
SOLUTION

60 kN/m 120 kN
25 kN
1
EI 2 EI 3

5m
1.5EI
1.5EI
Y
2 4 5

0
X 6m 6m

Fig. 1 Given Rigid Frame

1 2 3
4 4 4

1 [1] [2]
2 3

[4]
[3]

Y
2 5
4

X
Discretisation and Numbering System -
DOFs, members [ ] and joints

DOFs = 4 (3 rotations + 1 side sway at joints 1, 2 and 3)

Member Properties

1
 Coordinates
Joints 12EI 6EI 4EI
(mm) Direction
L EI L3 2EI
Member Cosines L2 L
(m) (relative) L
l m (EI) (EI) (EI)
End 1 End 2 x1 x2 y1 y2
(EI)

(1) [1] [2] 0 6000 5000 5000 6 1 0 1.0 0.056 0.167 0.667 0.333

(2) [2] [3] 6000 12000 5000 5000 6 1 0 1.0 0.056 0.167 0.667 0.333

(3) [4] [2] 6000 6000 0 5000 5 0 1 1.5 0.144 0.36 1.2 0.6

(4) [5] [3] 12000 12000 0 5000 5 0 1 1.5 0.144 0.36 1.2 0.6

The member force-displacement relationship is in the form;

f   k d

 12EI 2 12EI 6EI 12EI 2 12EI 6EI 

 L3 m  lm  m  m lm  m
L3 L2 L3 L3 L2
 
  12EI lm 12EI 2
l
6EI
l
12EI
lm
12 EI
 3 l2
6EI 
l
 L3 L3 L2 L3 L L2 
 6EI 6EI 4EI 6EI 6EI 2EI 
  2 m l m  2 l
 k    L L2 L L2 L L 

12EI 12EI 6EI 12EI 2 12EI 6EI
 3 m2 lm m m  3 lm m 
 L L3 L2 L3 L L2 
 12EI 12EI 6EI 12EI 12EI 2 6EI 
 3
lm  3 l2  2 l  3 lm l  2 l
 L L L L L3 L 
 6EI 6EI 2EI 6EI 6EI 4EI 
  L2 m L2
l
L L2
m  2 l
L L 

Member 1

 f1x  F4  0 0 0 0 0 0   d1x  D4 
d1y= 0
f  0  
d2y= 0
 1y  0.0556 0.1667 0 0.0556 0.1667   d1y  0 
d1z= D1 d2z= D2  
 f1z  F1  0 0.1667 0.667 0 0.1667 0.333   d1z  D1 
   EI    
d1x= D4 (1) d2x= D4
f 2x  F4  0 0 0 0 0 0  d 2x  D4 
 f 2y  0 0.0556 0.1667 0 0.0556 0.1667   d 2y  0 
     
 f 2z  F2  0 0.1667 0.333 0 0.1667 0.667   d 2z  D2 

Giving, in condensed form

 F4  0 0 0 0  D4 
F    D 
 1   EI 0 0.667 0 0.333  1
 F4  0 0 0 0  D4 
     
 F2  0 0.333 0 0.667  D2 

2
 Member 2
 f1x  F4  0 0 0 0 0 0   d1x  D4 
f  0  
 1y  0.0556 0.1667 0 0.0556 0.1667   d1y  0 
 
 f1z  F2  0 0.1667 0.667 0 0.1667 0.333   d1z  D2  [2]
   EI    
f 2x  F4  0 0 0 0 0 0  d 2x  D4 
 f 2y  0 0.0556 0.1667 0 0.0556 0.1667   d 2y  0 
     
 f 2z  F3  0 0.1667 0.333 0 0.1667 0.667   d 2z  D3 

d1y= 0
d2y= 0
d1z= D2 d2z= D3

d1x= D4 (2) d2x= D4

Giving,

 F4  0 0 0 0   D4 
F  0 0.667 0 0.333 D 
 2   EI    2
 F4  0 0 0 0   D4 
     
 F3  0 0.333 0 0.667   D3 

Member 3

d2y= 0  f1x   0.144 0 0.36 0.144 0 0.36  d1x  0 

 f   0  
d2z= D2
 1y  0 0 0 0 0   d1y  0 
 
d2x= D4  f1z   0.36 0 1.2 0.36 0 0.6   d1z  0 
   EI    
(3) f 2x  F4   0.144 0 0.36 0.144 0 0.36  d 2x  D 4 
 f 2y   0 0 0 0 0 0   d 2y  0 
     
d1y= 0
 f 2z  F2   0.36 0 0.6 0.36 0 1.2   d 2z  D2 
d1z=0

d1x= 0

Giving,  F4  0.144 0.36  D4 

 F  =EI  0.36 1.2   D2 
 2 

3
 Member 4

d1y= 0

d1z= D3
 f1x   0.144 0 0.36 0.144 0 0.36  d1x  0 
 f   0  
0   d1y  0 
d1x= D4
 1y  0 0 0 0
 
(4)  f1z   0.36 0 1.2 0.36 0 0.6   d1z  0 
   EI    
d2y= 0 f 2x  F4   0.144 0 0.36 0.144 0 0.36  d 2x  D4 
d2z=0  f 2y   0 0 0 0 0 0   d 2y  0 
     
 f 2z  F3   0.36 0 0.6 0.36 0 1.2   d 2z  D3 
d2x= 0

Giving,  F4  0.144 0.36  D4 

 F  =EI  0.36 1.2   D3 
 3 

Assembling the global stiffness matrix [K]

0.667 0.333 0 0 
 0.333 2.534 0.333 0.36 
 K   EI  
 0 0.333 1.867 0.36 
 
 0 0.36 0.36 0.288

Assembling the global external force vector [F]:

-180 kNm 180 kNm -90 kNm 90 kNm

25 kN
3
180 kN (1) (2)
180 kN 60 kN 60 kN

(4)
(3)

(a) Free-body diagrams of members

4
Member (1)

wl 2 60(6) 2
FEM=    180 kNm
12 12

Member (2)
Pl 120(6)
FEM=    90 kNm
8 8
Global Force vector  F  becomes

F1   ( 180)  180 kNm

F2  (180  90)  90 kNm
F3  90 kNm
F4  25 kN
 180 
 90 
 F   
 90 
 
 25

The Structure force-displacement relationship is

 F =  K  D

180  0.667 0.333 0 0   D1 

 90   0.333 2.534 0.333 0.36  D 
103    EI    2
 90   0 0.333 1.867 0.36   D3 
     
 25  0 0.36 0.36 0.288 D4 

Solve for Displacements

1
 D1  0.667 0.333 0 0  180 
D   0.333 2.534 0.333 0.36   90 
 2   EI     103
 D3   0 0.333 1.867 0.36   90 
     
 D4   0 0.36 0.36 0.288  25 