Chapter 05a
Chapter 05a
Learning Objectives
• To derive the two-dimensional arbitrarily oriented
beam element stiffness matrix
• To demonstrate solutions of rigid plane frames by
the direct stiffness method
• To describe how to handle inclined or skewed
supports
The local axes and are located along the beam element and
transverse to the beam element, respectively, and the global
axes x’ and y’ are located to be convenient for the total
structure.
S C 0 0 0
0
0 0 1 0 0 0
T
0 0 0 S C 0
0 0 0 0 0 1
u1 v1 1 u2 v2 2
12S 2 12SC 6LS 12S 2 12SC 6LS
12SC 12C 2 6LC 12SC 12C 2 6LC
EI 6LS 6LC 4L2 6LS 6LC 2L2
k 3
L 12S 2 12SC 6LS 12S 2 12SC 6LS
12SC 12C 2 6LC 12SC 12C 2 6LC
6LS 6LC 2L2 6LS 6LC 4L2
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 5/52
f1x C 1 0 0 C 1 0 0 u 1
f 0 12C 2 6 LC 2 0 12C 2 6 LC 2 v 1
1y
m 1 0 6 LC 2 4C 2 L2 0 6 LC 2 2C 2 L2 1
f2x C 1 0 0 C1 0 0 u 2
f2y 0 12C 2 6 LC 2 0 12C 2 6 LC 2 v 2
m 2 0 6 LC 2 2C 2 L2 0 6 LC 2 4C 2 L2 2
AE EI
C1 C2
L L3
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 6/52
u1 C S 0 0 0 0 u1
v S C 0 0 0 0 v
1 1
1 0 0 1 0 0 0 1 d Td
u2 0 0 0 C S 0 u2
v 2 0 0 0 S C 0 v 2
2 0 0 0 0 0 1 2
L L L L L L
12I 6I 12I 12I 2 6I
AS 2 C 2
2
C A 2 CS AS 2 C
2
C
L L L L L
6I 6I
4I C 2I
k L
E S
L L
12I 12I 6I
AC 2 2 S 2 A CS S
L L2 L
12I 6I
symmetric AS 2 2 C 2 C
L L
4I
Let E = 30 x 106 psi and A = 10 in2 for all elements, and let I =
200 in4 for elements 1 and 3, and I = 100 in4 for element 2.
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 10/52
C 0 S 1
E 30 106
250,000 lb 3
L 120 in
u1 v1 1 u2 v2 2
0.167 0 10 0.167 0 10
0 10 0 0 10 0
10 0 800 10 0 400 lb
k (1) 250,000
0.167 0 10 0.167 0 10 in
0 10 0 0 10 0
10 0 400 10 0 800
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 11/52
C 1 S0
E 30 106
250,000 lb 3
L 120 in
u2 v2 2 u3 v3 3
10 0 0 10 0 0
0 0.0835 5 0 0.0835 5
0 5 400 0 5 200 lb
k (2) 250,000
10 0 0 10 0 0 in
0 0.0835 5 0 0.0835 5
0 5 200 0 5 400
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 12/52
C 0 S 1
E 30 106
250,000 lb 3
L 120 in
u3 v3 3 u4 v4 4
0.167 0 10 0.167 0 10
0 10 0 0 10 0
10 0 800 10 0 400 lb
k (3) 250,000
0.167 0 10 0.167 0 10 in
0 10 0 0 10 0
10 0 400 10 0 800
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 13/52
x
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 14/52
0 1 0 0 0 0 u1 0 0
1 0 0 0
0 0 v1 0
0
0 0 1 0 0 0 1 0 0
Td u 0.211 in
0 0 0 0 1 0 2 0.00148 in
0 0 0 1 0 0 v 2 0.00148 in 0.211 in
0 0 0 0 0 1 2 0.00153 rad 0.00153 rad
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 15/52
C1 0 0 C1 0 0 C1
2
AE 10 in 30 10 psi
6
2.5 106 lb in
0 6LC2
12C2 6LC2 0 12C2
L 120 in
0 6LC2 4C2L2 0 6LC2 2C2L2
k
C1 0 0 0 0
C1
EI 200 in 30 10 psi
4 6
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 2C2L2 0 6LC2
4C2L2 C2 3,472.22 lb in
L3 120 in
3
f1x 3,700 lb
f
1y 4,990 lb
m1 376 k in
f2x 3,700 lb
f2y 4,990 lb
m2 223 k in
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 16/52
1 0 0 0 0 0 u2 0.211 in 0.211 in
0
1 0 0 0 0 v 2 0.00148 in 0.00148 in
0 0 1 0 0 0 2 0.00153 rad 0.00153 rad
Td
0 0 0 1 0 0 u3 0.209 in 0.209 in
0 0 0 0 1 0 v 3 0.00148 in 0.00148 in
0 0 0 0 0 1 3 0.00149 rad 0.00149 rad
C1 0 0 C1 0 0 C1
2
AE 10 in 30 10 psi
6
2.5 106 lb in
0 6LC2
12C2 6LC2 0 12C2
L 120 in
0 6LC2 4C2L2 0 6LC2 2C2L2
k
C1 0 0 0 0
C1
EI 100 in 30 10 psi
4 6
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 2C2L 2
0 6LC2
4C2L2 C2 1,736.11 lb in
L3 120 in
3
f2x 5,010 lb
f 3,700 lb
2y
m2 223 k in
f3x 5,010 lb
f3y 3,700 lb
m3 221 k in
0 1 0 0 0 u3 0.209 in
0 0.00148 in
1 0 0 0 0 v 3 0.00148 in 0.209 in
0
0 0 1 0 0 3 0.00149 rad 0.00149 rad
0
Td
0 0 0 0 1 0 u4 0 0
0 0 0 1 0 0 v 4 0 0
0 0 0 0 0 1 4 0 0
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 18/52
C1 0 0 C1 0 0 C1
2
AE 10 in 30 10 psi
6
2.5 106 lb in
0 6LC2
12C2 6LC2 0 12C2
L 120 in
0 6LC2 4C2L2 0 6LC2 2C2L2
k
C1 0 0 0 0
C1
EI 200 in 30 10 psi
4 6
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 2C2L2 0 6LC2
4C2L2 C2 3,472.22 lb in
L3 120 in
3
f3x 3,700 lb
f 5,010 lb
3y
m3 226 k in
f4x 3,700 lb
f4y 5,010 lb
m4 375 k in
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 19/52
Let E = 30 x 106 psi and A = 100 in2 for all elements, and
let I = 1,000 in4 for all elements.
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 20/52
wL (1,000 lb / ft )40 ft
f2 y f3 y 20k
2 2
wL2 (1,000 lb / ft )(40 ft )2
m2 m3 133,333 lb ft
12 12
1,600 k in
C 0.707 S 0.707
E 30 106
58.93 k 3
L 12 30 2 in
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 21/52
u2 v2 2
2,948 2,945 491
k (1)
2,945 2,948 491 k
in
491 491 235,700
C 1 S 0
E 30 106
52.5 k 3
L 12 40 in
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 22/52
u2 v2 2
6,250 0 0
k (2)
0 3.25 781.25 k
in
0 781.25 250,000
u2 0.0033 in
v 2 0.0097 in
0.0033 rad
2
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 23/52
u2 0.0033 in
v 2 0.0097 in
0.0033 rad
2
x
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 24/52
Element 1 Element 2
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 27/52
C 0.447 S 0.895
12I 12(800)
0.0334 in 2
12 44.7
2 2
L
u v
6I 6(800) 4 4 4
8.95 in 3 90.0 178 448
L 12 44.7
k (1)
178 359 244 k in
E 30 10 6
448 244 179,000
55.9 k 3
L 12 44.7 in
C 0.447 S 0.895
12I 12(800)
0.0334 in 2
12 44.7
2 2
L
u v
6I 6(800) 4 4 4
8.95 in 3 90.0 178 448
L 12 44.7
k (2) 178 359 244 k in
E 30 10 6
448 244 179,000
55.9 k 3
L 12 44.7 in
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 32/52
C 1 S0
12I 12(800)
0.0267 in 2
12 50
2 2
L
u v
6I 6(800) 4 4 4
8.0 in 3 400 0 0
L 12 50
k (3)
0 1.334 400 k in
E 30 10 6
0 400 160,000
50.0 k 3
L 12 50 in
u4 0.0103 in
v 4 0.000956 in
0.00172 rad
4
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 33/52
1.67 k k(1) Td
0.88 k
158 k in
kd
f(1)
1.67 k
0.88 k
311 k in
2.44 k k(2) Td
0.877 k
158 k in
kd
f(2)
2.44 k
0.877 k
312 k in
Therefore, the above equations are the final local nodal forces
2.44 k
0.877 k
158 k in
kd
f(2)
2.44 k
0.877 k
312 k in
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 37/52
1 0 0 0 0 0 4
u0.0103 in 0.0103 in
0
1 0 0 0 0 v 0.000956
4 in 0.000956 in
0 0 1 0 0 0 0.00172
4 rad 0.00172 rad
Td
0 0 0 1 0 0 u0
3 0
0 0 0 0 1 0 v0
3
0
0 0 0 0 0 1 30 0
4.12 k k(3) Td
0.687 k
275 k in
kd
f(3)
4.12 k
0.687 k
137 k in
Therefore, the above equations are the final local nodal forces
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 39/52
For the bar, A = 1 x 10-3 m2, for the beam, A = 2 x 10-3 m2,
I = 5 x 10-5 m4, and L = 3 m. Let E = 210 GPa for both
elements.
Bar
Beam
C 1 S 0
12I 12(5 10 5 )
6.67 105 m 2
L2 (3)2
6I 6(5 10 5 ) u1 v1 1
104 m3
L 3
2 0 0
k (1)
70 10 0
3
0.067 0.10 kN
E 210 106 m
70 106 kN / m 3
L 3 0 0.10 0.20
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 40/52
C 0.707 S 0.707
u1 v1
k (2)
103 m 2 210 106 kN m2 0.5 0.5
0.5 0.5 m
kN
4.24 m
u1 v1
0.354 0.354 kN
k (2) 70 103 m
0.354 0.354
2.354 0.354 0
K 70 10 0.354 0.421 0.10 kN m
3
0 0.10 0.20
u1 0.00388 m
v1 0.0225 m
0.0113 rad
1
AE
f1x Cu1 Sv1 670 kN
L
AE
f3x Cu1 Sv1 670 kN
L
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 42/52
1 0 0 0 0 1
0 u0.00388 m 0.00388 m
0
1 0 0 0 0 v0.0225
1 m 0.0225 m
0 0 1 0 0 0 0.0113
1 rad 0.0113 rad
d Td
0 0 0 1 0 0 u02 0
0 0 0 0 1 0 v 02
0
0 0 0 0 0 1 20 0
C1 0 0 C1 0 0
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 4C2L2 0 6LC2 2C2L2
k
1C 0 0 C1 0 0
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 2C2L2 0 6LC2 4C2L2
AE (0.002)210 106
C1 140 103 kN
m
L 3
EI 210 106 (5 10 5 )
C2 388.89 kN m
3
3
L3
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 43/52
Assume A = 2 x 10-2 m2, I = 2 x 10-4 m4, and E = 210 GPa for all
elements.
0.015 0 0.03
The parts of k associated 5
k 5.25 10 0
(1)
2 0 kN m
with node 2 are:
0.03 0 0.08
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 45/52
2 0 0
The parts of k associated (2) 5
k 4.2 10 0 0.0096 0.024 kN m
with node 2 are:
0 0.024 0.08
0.8480 0 0.0158
K 10 06
1.0500 0.0101 kN m
0.0158 0.0101 0.0756
u2 4.95 10 m
6
6
v 2 2.56 10 m
2.66 10 4 rad
2
0 1 0 0 0 0 u01 0
1 0 0 0 0 0 v01 0
0 0 1 0 0 0 01 0
d Td 6 6
0 0 0 0 1 0 u4.95
2 10 m 2.56 10 m
0 0 0 1 0 0 v2.56
2
10 6 m 4.95 10 6 m
0 0 0 0 0 1 2.66
2 10 4 rad 2.66 10 4 rad
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 47/52
C1 0 0 C1 0 0
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 4C2L2 0 6LC2 2C2L2
k
1C 0 0 C1 0 0
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 2C2L2 0 6LC2 4C2L2
EI 210 106 (5 10 5 )
C2 388.89 kN m
3
3
L3
1 0 0 0 0 0 u24.95
10 6 m 4.95 10 6 m
0 1 0 0 0 0 v22.56
6
10 m 2.56 10 m
6
0 0 1 0 0 0 22.66
10 4 rad 2.66 10 4 rad
d Td u
0 0 0 1 0 0 30 0
0 0 0 0 1 0 v 30 0
0 0 0 0 0 1 30 0
C1 0 0 C1 0 0
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 4C2L2 0 6LC2 2C2L2
k
1C 0 0 C1 0 0
0 12C2 6LC2 0 12C2 6LC2
0 6LC2 2C2L2 0 6LC2 4C2L2
AE (2 10 2 )210 106
C1 0.84 106 kN
m
L 5
EI 210 106 (2 10 4 )
C2 336 kN m
5
3
L3
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 49/52
cos sin 0
d '3 [t3 ]d3 t3 sin cos 0
0 0 1
CIVL 7/8117 Chapter 5 - Plane Frame and Grid Equations - Part 1 51/52
Both the identity matrix [I] and the matrix [t3] are 3 x 3 matrices.
F1x u1
F v
1y 1
M1 1
F2 x u2
F2 y [Ti ][K ][Ti ]
T
v2
M
2 2
F '3 x u '3 Elemental
F ' v ' coordinates
3y 3
M3 3
End of Chapter 5a