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