HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Problem 1: A plane frame is fixed at the
C p0 D
ends as shown in Figure 1. Assume that:
E = 2x1011 Pa, ν 0.3 , L = 1 m,
p0 15kN / m M = p0L2
A 100 cm2 ,J 1000 cm 4 L/2
Find:
y
a) Displacements and rotation at C. B L
b) Reaction forces at A and D.
c) Plot the bending moment diagram. x
A α 450
Figure 1.
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Solution
Discretize the model as shown in the figure below.
q8 q5
The Boolean matrix:
Node i Node j q7 q4
3
Element 1 2 3 4 5 6 q9 2
2 q6
1 1 2 3 7 8 9
2 7 8 9 4 5 6
1
q2
1 q1
q3
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Element stiffness matrix:
- Element 1: 1 2 3 7 8 9
A B' A B' B 'L A B' B ' A B 'L 1
2 2 2 2 2 2 2 2 2 2 2 2
A B' B 'L B ' A A B' B 'L 2
2 2 2 2 2 2 2 2 2 2
4J B 'L B 'L 2J 3
E
K 1 L 2 2 2
A B'
2 2
A B'
2
B 'L
7
2 2 2 2 2 2
A B' B 'L
8
2 2 2 2
4J 6J
sym 9 where: B '
2 L2
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
- Element 2:
7 8 9 4 5 6
A 0 0 A 0 0 7
BL BL
B 0 B 8
2 2
BL
E 4J 0 2J 9
K 2 L 2
A 0 0 4
B
BL
5
2
sym 4J 6
12J
where: B
L2
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Global stiffness matrix:
A B' A B' B 'L A B' B ' A B 'L
0 0 0
2 2 2 2 2 2 2 2 2 2 2 2
A B' B 'L B ' A A B' B 'L
0 0 0
2 2 2 2 2 2 2 2 2 2
4J B 'L B 'L 2J
0 0 0
2 2 2 2 2 2
A 0 0 A 0 0
BL BL
E B B
L
K 2
0
2
BL
4J 0 2J
2
A B' A B' B 'L
A
2 2 2 2 2 2
A B' B 'L B L
B
2 2 2 2 2
4J
sym 4J
2
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Load vector:
- Load vector of Element 1:
3M 1
2 4L
0 0 0
2 3M 3M 2
2 2 2L 4L
0 0 0
2 M M 3
0 1 0 0 4 4
P 1
2 3M 3M
0 0 0 7
2 2 2L 4L
3M
2 M
0 0 0 8
2 4 4L
M
0 0 0 1 9
4
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
- Load vector of Element 2:
p0 L
1 0 2 7
0 0 p L 0 8
0
0 0 2 0 9
P 2
0 1 p0 L p0 L 4
0 0 2 2
0 5
0 0
0 6
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
- Nodal load vector:
X1
Y2
Y1
M1
3 2 X2
X2
2 M2 P n Y2
M
2
0
1 0
0
Y1
X1
1
M1
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
- Global load vector:
3M 3M
X
1 4L
4L
3 M Y 3 M
4L X1 1
4L
M Y1 M M
4 M1 1
4
p0 L p0 L
X2 X2
2 2
P 0 Y2
Y2
M
0 2
M2
0
3M p L 3M p L
0
0 0
4L 2 4L 2
3 M 0 3 M
4 L 4L
M M
4 4
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
We obtain the system of equations:
A B' A B' B 'L A B' B ' A
B 'L
0 0 0 3M
2 2 2 2 2 2 2 2 2 2 2 2
X
4L
1
A B' B 'L B ' A A B' B 'L
0 0 0 3M
2 2 2 2 2 2 2 2 2 2 q1 Y1 4L
4J B 'L B 'L 2J
0 0 0 q 2 M M
2 2 2 2 2 2 q 1
4
A 3
A 0 0 0 0
q4 X p0 L
E B
BL
0 B
BL 2 2
2 2 5
q
L Y2
BL
4J 0 2J 6
q M2
2
q
3M p L
A B' A B' B 'L 7
0
A q 4L 2
2 2 2 2 2 2 8
3M
A B' B 'L B L q9
B 4L
2 2 2 2 2 M
4J
sym 4J 4
2
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Impose boundary conditions q1 q2 q4 q5 0 (m ) and q3 q6 0 (rad) , the
equation becomes:
A B' A B' B 'L 3M p0 L
4L 2
A
2 2 2 2 2 2
q7
E A B' A B ' BL B 'L 3M
B q8
L 2 2 2 2 2 2 2 4L
B 'L q
M
BL B 'L 9
(4 2 2 )J
2 2 2 2 2 4
a) Solve the equation and substitute the values of the parameters, we get:
q7 0.0149 10 3 m
3
q8 0.0326 10 m
0.2885 10 3 rad
9
q
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Displacement results:
q 0
T
3 3 3
0 0 0 0 0 0.0149 10 0.0326 10 0.2885 10
b) Reaction forces:
E A B' E B ' A EB ' 3M
X1 q7 q8 q9 22328 (N )
L 2 2 L 2 2 2 2 4L
E B ' A E A B' EB ' 3M
Y1 q7 q8 q9 2679 (N )
L 2 2 L 2 2 2 2 4L
EB ' EB ' 2EJ M
M1 q7 q8 q9 2732(Nm )
2 2 2 2 L 4
EA p0 L
X2 q7 37328(N )
L 2
EB EB EB 2EJ
Y2 q8 q9 2679(N ) M2 q8 q9 762(Nm )
L 2 2 L
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
c) Bending moment:
Element 1:
0
3 2L 0
2L2
3 2L 3 2L 3 2L
2 2L2 1017
EJ 2 2 2 2 0
M1 [S ']1 {q }1 3 3
(Nm )
L 3 2L 3 2L 3 2L 3 2L 0.0149 10 1833
2L2 2 2L2
0.0326 10
3
2 2 2 2
3
0.2885 10
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Element 2:
0.0149 10 3
3
0.0326 10
EJ 0 6L 4L2 0 6L 2L2 0.2885 10 3 1916
M 2 [S ']2 {q } 2 3 2 (Nm )
L 0 6L 2L2 0 6L 4L 0 762
0
0
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Bending moment diagram:
1833 762
7500
+
1916 3750
3750
7500
1017
TP. Hồ Chí Minh, 10/29/2023
HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY
FACULTY OF APPLIED SCIENCE
DEPARTMENT OF ENGINEERING MECHANICS
Bending moment diagram:
1833 762
1916 1917
2733
TP. Hồ Chí Minh, 10/29/2023