ANALYSIS OF BEAMS USING MATRIX STIFFNESS METHOD
Determine the joint displacements, member end forces, and support reactions for the three-span
continuous beam as shown below. Draw the analytical model. Adapt solutions of sample problem.
Joint Details
Degrees of
Joint x (m) y (m) Freedom
1 0 0 5 6
2 6 0 1 2
3 10 0 7 3
4 20 0 8 4
Member Details
Beginning End Length
Member Joint Joint Xe Xb Ye Yb (m)
1 1 2 6 0 0 0 6
2 2 3 10 6 0 0 4
3 3 4 20 10 0 0 10
ANALYSIS OF BEAMS USING MATRIX STIFFNESS METHOD
Member 1 Member Stiffness Matrix [k]
E 28 Gpa 5 6 1 2
I 8700000000 mm^4 13533.3 40600 -13533.3 40600 5
L 6 m 40600 162400 -40600 81200 6
EI/L^3 1127.777778 kN/m -13533.3 -40600 13533.33 -40600 1
40600 81200 -40600 162400 2
Fixed-end Force Vector {Qf}
FSb 63 Kn 5
FMb 54 Kn-m 6
FSe 27 kN 1
FMe -36 kN-m 2
Member 2 Member Stiffness Matrix [k]
E 28 Gpa 1 2 7 3
I 5800000000 mm^4 30450 60900 -30450 60900 1
L 4 m 60900 162400 -60900 81200 2
EI/L^3 2537.5 kN/m -30450 -60900 30450 -60900 7
60900 81200 -60900 162400 3
Fixed-end Force Vector {Qf}
FSb 0 Kn 1
FMb 0 Kn-m 2
Fse 0 kN 7
Fme 0 kN-m 3
ANALYSIS OF BEAMS USING MATRIX STIFFNESS METHOD
Member 3 Member Stiffness Matrix [k]
E 28 Gpa 7 3 8 4
I 5800000000 mm^2 1948.8 9744 -1948.8 9744 7
L 10 m 9744 64960 -9744 32480 3
EI/L^3 162.4 kN/m -1948.8 -9744 1948.8 -9744 8
9744 32480 -9744 64960 4
Fixed-end Force Vector {Qf}
FSb 75 Kn 7
FMb 187.5 Kn-m 3
Fse 75 kN 8
Fme -187.5 kN-m 4
Structure
Stiffness Matrix
1 2 3 4
43983.3333 20300 60900 0 1
20300 324800 81200 0 2
S=
60900 81200 227360 32480 3
0 0 32480 64960 4
Joint
Load
Vector
P
-200 1 kN
0 2 kN-m
-90 3 kN-m
0 4 kN-m
Structure Fixed
Load Vector
Pf
27 1 kN
-36 2 kN
187.5 3 kN
-187.5 4 kN
ANALYSIS OF BEAMS USING MATRIX STIFFNESS METHOD
Solving for the displacement
matrix:
-200 27 43983 20300 60900 0 d1
0 -36 20300 324800 81200 0 d2
- = x
-90 187.5 60900 81200 227360 32480 d3
0 -187.5 0 0 32480 64960 d4
-227 43983 20300 60900 0 d1
36 20300 324800 81200 0 d2
= x
-277.5 60900 81200 227360 32480 d3
187.5 0 0 32480 64960 d4
d1 -0.0044729 m
d2 0.00056143 rad
d3 -0.0006842 rad
d4 0.00322847 rad
Member 1
{u} Member End Displacement {u}
u1 0 5
u2 0 6
u3 -0.004472863 1
u4 0.000561433 2
{Q} - {Qf} = [k]{u} Local Member End Forces {Q} - {Qf}
5 6 1 2
Q1-Qf1 83.32690996 kN 13533.33 40600 -13533.3 40600 5 0
Q2-Qf2 227.1865621 kN-m 40600 162400 -40600 81200 6 0
Q3-Qf3 -83.32690996 kN -13533.3 -40600 13533.33 -40600 1 -0.00447
Q4-Qf4 272.7748977 kN-m 40600 81200 -40600 162400 2 0.000561
{Q} = [k]{u} + {Qf} Local Member End Forces {Q}
Q1 146.32691 kN 5
Q2 281.1865621 kN-m 6
Q3 -56.32690996 kN 1
Q4 236.7748977 kN-m 2
ANALYSIS OF BEAMS USING MATRIX STIFFNESS METHOD
Member 2
{u} Member End Displacement {u}
u1 -0.004472863 1
u2 0.000561433 2
u3 0 7
u4 -0.000684165 3
{Q} - {Qf} = [k]{u} Local Member End Forces {Q} - {Qf}
1 2 7 3
Q1-Qf1 -143.67309 kN 30450 60900 -30450 60900 1 -0.00447
Q2-Qf2 -236.7748977 kN-m 60900 162400 -60900 81200 2 0.000561
Q3-Qf3 143.67309 kN -30450 -60900 30450 -60900 7 0
Q4-Qf4 -337.9174625 kN-m 60900 81200 -60900 162400 3 -0.00068
{Q} = [k]{u} + {Qf} Local Member End Forces {Q}
Q1 -143.67309 kN 1
Q2 -236.7748977 kN-m 2
Q3 143.67309 kN 7
Q4 -337.9174625 kN-m 3
Member 3
{u} Member End Displacement {u}
u1 0 7
u2 -0.000684165 3
u3 0 8
u4 0.003228474 4
{Q} - {Qf} = [k]{u} Local Member End Forces {Q} - {Qf}
7 3 8 4
Q1-Qf1 24.79174625 kN 1948.8 9744 -1948.8 9744 7 0
Q2-Qf2 60.41746248 kN-m 9744 64960 -9744 32480 3 -0.00068
Q3-Qf3 -24.79174625 kN -1948.8 -9744 1948.8 -9744 8 0
Q4-Qf4 187.5 kN-m 9744 32480 -9744 64960 4 0.003228
{Q} = [k]{u} + {Qf} Local Member End Forces {Q}
Q1 99.79174625 kN 7
Q2 247.9174625 kN-m 3
Q3 50.20825375 kN 8
Q4 0 kN-m 4
ANALYSIS OF BEAMS USING MATRIX STIFFNESS METHOD
Reaction
Vector
Reactions
146.3269 5 kN
281.1866 6 kN-m
243.4648 7 kN
50.20825 8 kN
