Bahir Dar University
Bahir Dar Institute of Technology (BiT)
Faculty of Civil and Water Resource Engineering
Post Graduate Program in Structural Engineering
Finite Element Method in Structures
Assignment-4
Prepared by: Natinael Wolde Azale
ID: BDU1022155
Submitted to: Ins. Seyfe Nigussie
Submission Date: July 14, 2020
Finite Element Method in Structures Assignment -4
#Q-1. Determine the nodal displacement, element stresses and support reactions of
the axially loaded bar as shown in Fig.1 Take E = 200 GPa and P = 30 kN
Figure-1
Solution
✓ Step-1: Assign the nodal point and the naming member (element) for each structure,
Number of node =Number of element +1, and there is only one degree of freedom at each
node. In this problem there are three element and four node.
Node-Element connectivity by table:
Element Node-1 Node-2
1 1 2
2 2 3
3 3 4
✓ Step-2: Compute element stiffens matrix for each element and assemble Global stiffness
matrix
Element-1 Element-2 Element-3
𝑬𝟏 𝑨𝟏 𝟏 −𝟏 𝐸2 𝐴2 1 −1 𝐸3 𝐴3 1 −1
[𝒌]𝟏 = [ ] [𝑘]2 = [ ] [𝑘]3 = [ ]
𝑳𝟏 −𝟏 𝟏 𝐿2 −1 1 𝐿3 −1 1
1 2 2 3 3 4
333.33 −333.33 𝟐 266.67 −266.67 𝟑
[𝒌]𝟏 = 𝟏𝟎𝟑 [ 𝟑𝟑𝟑. 𝟑𝟑 −𝟑𝟑𝟑. 𝟑𝟑 𝟏 3
[𝒌]𝟐 = 10 [ 3
] ] [𝒌]𝟑 = 10 [ ]
−𝟑𝟑𝟑. 𝟑𝟑 𝟑𝟑𝟑. 𝟑𝟑 𝟐 −333.33 333.33 𝟑 −266.67 266.67 𝟒
Master of Structural Engineering Page 1
Finite Element Method in Structures Assignment -4
✓ Assemble Global stiffens matrix
[𝑲] = [𝒌]𝟏 + [𝒌]𝟐 + [𝒌]𝟑
𝟏 𝟐 𝟑 𝟒
𝟑𝟑𝟑. 𝟑𝟑𝟑 −𝟑𝟑𝟑. 𝟑𝟑𝟑 𝟎 𝟎 𝟏
−𝟑𝟑𝟑. 𝟑𝟑 𝟔𝟔𝟔. 𝟔𝟔𝟔 −𝟑𝟑𝟑. 𝟑𝟑 𝟎 𝟐
[𝑲] = 𝟏𝟎𝟑 [ ]
𝟎 −𝟑𝟑𝟑. 𝟑𝟑 𝟔𝟎𝟎. 𝟎𝟎𝟎 −𝟐𝟔𝟔. 𝟔𝟕 𝟑
𝟎 𝟎 −𝟐𝟔𝟔. 𝟔𝟕 𝟐𝟔𝟔. 𝟔𝟕 𝟒
✓ Step-3: Determine the nodal displacement, element stresses and support reactions by
applying finite element formulation
I. Nodal Displacement
[𝐾]{𝐷} = {𝐹}
333.333 −333.333 0 0 𝑢1 𝐹1
−333.33 666.666 −333.33 0 𝑢2 𝐹2
𝟏𝟎𝟑 [ ]{ } = { }
0 −333.33 600.000 −266.67 𝑢3 𝐹3
0 0 −266.67 266.67 𝑢4 𝐹4
By applying boundary conditions 𝑢1 = 𝑢4 = 0
𝑢1 0
𝑢2 𝑢2
{𝑢} = {𝐷} = {𝑢 } = { }
3 𝑢3
𝑢4 0
Nodal Force Vector are
𝐹1 0 0
𝐹
{𝐹} = { 2 } = {𝑃} = 𝟏𝟎𝟑 {30}
𝐹3 0 0
𝐹4 0 0
333.333 −333.333 0 0 0 0
−333.33 𝟔𝟔𝟔. 𝟔𝟔𝟔 −𝟑𝟑𝟑. 𝟑𝟑 0 𝑢2 𝟑𝟎
𝟏𝟎𝟑 [ ] { } = 𝟏𝟎𝟑 { }
0 −𝟑𝟑𝟑. 𝟑𝟑 𝟔𝟎𝟎. 𝟎𝟎𝟎 −266.67 𝑢3 𝟎
0 0 −266.67 266.67 0 0
Reduced matrix equation
666.666 −333.33 𝑢1 30
[ ] {𝑢 } = { }
−333.33 600.00 2 0
By using Cramer’s Rule
30 −333.33
| | 18,000
𝒖𝟏 = 0 600.00 = = 𝟎. 𝟎𝟔𝟐𝟑𝟎𝟖 𝒎𝒎
666.666 −333.33
| | 288,888.89
−333.333 600.00
Master of Structural Engineering Page 2
Finite Element Method in Structures Assignment -4
666.666 30
| | 10,000
𝑢2 = −333.33 0 = = 𝟎. 𝟎𝟑𝟒𝟔𝟏𝟓 𝒎𝒎
666.666 −333.33
| | 288,888.89
−333.333 600.00
𝑢1 0
𝑢2 0.062308
∴ 𝑁𝑜𝑑𝑎𝑙 𝑑𝑖𝑠𝑝𝑙𝑎𝑐𝑒𝑚𝑒𝑛𝑡 𝑖𝑠 {𝑢} = {𝐷} = {𝑢 } = { }
3 0.034615
𝑢4 0
II. Determine reaction
𝑢1 0
𝑢2 3 0.062308
❖ 𝑅1 = [𝑘11 𝑘12 𝑘13 𝑘14 ] {𝑢 } = 10 [333.33 −333.33 0 0] { }
3 0.034615
𝑢4 0
∴ 𝑅1 = 103 [−333.333 ∗ 0.062308] = −𝟐𝟎. 𝟕𝟕 𝑲𝑵
𝑢1 0
𝑢2
❖ 𝑅2 = [𝑘41 𝑘42 𝑘43 𝑘44 ] {𝑢 } = 103 [0 0 −266.67 266.67] {0.062308}
3 0.034615
𝑢4 0
∴ 𝑅2 = 103 [−266.67 ∗ 0.034615] = −𝟗. 𝟐𝟑𝑲𝑵
III. Determine Stress Element
1. Stress on Element-1
𝑑𝑢 𝐸1 𝑢1 2𝑋105 0
(𝝈)𝟏 = 𝐸1 = [−1 1] {𝑢 } = [−1 1] { } = 𝟖𝟑. 𝟎𝟖 𝑵/𝒎𝒎𝟐
𝑑𝑥 𝐿 2 150 0.062308
2. Stress on Element-2
𝑑𝑢 𝐸2 𝑢2 2𝑋105
(𝝈)𝟐 = 𝐸2 = [−1 1] {𝑢 } = [−1 1] {0.062308} = −𝟑𝟔. 𝟗𝟐𝑵/𝒎𝒎𝟐
𝑑𝑥 𝐿 3 150 0.034615
3. Stress on Element-3
𝑑𝑢 𝐸3 𝑢3 2𝑋105
(𝝈)𝟑 = 𝐸3 = [−1 1] {𝑢 } = [−1 1] {0.034615} = −𝟐𝟑. 𝟎𝟖 𝑵/𝒎𝒎𝟐
𝑑𝑥 𝐿 4 150 0
Master of Structural Engineering Page 3
Finite Element Method in Structures Assignment -4
#Q-2: Determine the forces in the members of the truss shown in Fig.2.Take E
= 200GPa, A = 2000 mm2.
Figure-2
Solution
✓ Step-1: Assign the nodal point and the naming member (element) for each structure,
Number of node =Number of element +1, and there is two degree of freedom at each node.
In this problem there are five element and four node.Hence total degree of freedom in this
problem is =Number of nodes X degree of freedom at each node = 4x2=8.
Master of Structural Engineering Page 4
Finite Element Method in Structures Assignment -4
Node-Element connectivity by table:
Element Node-1 Node-2 Angle(θ) l(cosθ) m(sinθ) L(m) EA/L(103 N/mm)
1 1 2 0 1.0 0.0 4.0 100.00
2 2 3 90 0.0 1.0 3.0 133.33
3 3 4 180 -1.0 0.0 4.0 100.00
4 1 4 90 0.0 1.0 3.0 133.33
5 1 3 36.68 0.8 0.6 5.0 80.00
✓ Step-2: Compute element stiffens matrix for each element and assemble Global stiffness
matrix
Compute element stiffens matrix
Formula Element-1 Element-2
1 2 3 4 3 4 5 6
𝒍𝟐 𝒍𝒎 −𝒍𝟐 −𝒍𝒎 1 3
𝑬𝒆 𝑨𝒆 𝒍𝒎 𝟐
100 0 −100 0 0 0 0 0
[𝒌𝒆 ] = [ 𝟐 𝒎 −𝒍𝒎 −𝒎𝟐 ]
[𝑘1 ]=103 [ 0 0 0 0 2 0 133.33 0 −333.33 4
𝒍𝒆 −𝒍 −𝒍𝒎 𝒍𝟐 𝒍𝒎 ] [𝑘2 ]=103 [ ]
−𝒍𝒎 −𝒎𝟐 𝒍𝒎 𝒎𝟐 −100 0 100 0 3 0 0 0 0 5
0 0 0 0 4 0 −333.33 0 333.33 6
Element-3 Element-4 Element-5
5 6 7 8 1 2 7 8 1 2 5 6
1 1 1
100 0 −100 0 0 0 0 0 51.2 38.4 −51.2 −38.4
2 2 2
[𝑘3 ]=10 [ 0
3 0 0 0
] [𝑘4 ]=10 [0
3 133.33 0 −133.33
] [𝑘5 ]=10 [ 38.4
3 28.8 −38.4 −28.8
]
−100 0 100 0 3 0 0 0 0 7 −51.2 −38.4 512 38.4 5
0 0 0 0 0 −133.33 0 133.33 −38.4 −28.8 38.4 28.8
4 8 6
Master of Structural Engineering Page 5
Finite Element Method in Structures Assignment -4
Assemble Global stiffens matrix
[𝑲] = [𝒌]𝟏 + [𝒌]𝟐 + [𝒌]𝟑 + [𝒌]𝟒 + [𝒌]𝟓 + [𝒌]𝟔 + [𝒌]𝟕 + [𝒌]𝟖
𝟏 𝟐 𝟑 𝟒 𝟓 𝟔 𝟕 𝟖
151.2 38.40 −100 0 −51.20 −38.4 0 0 𝟏
38.40 162.13 0 0 −38.40 −28.8 0 −133.33 𝟐
−100 0 100 0 0 0 0 0 𝟑
[𝑲] = 𝟏𝟎𝟑 0 0 0 133.33 0 −133.33 0 0 𝟒
−51.2 −38.4 0 0 151.20 38.40 −100 0 𝟓
−38.40 −28.8 0 −133.33 38.40 162.13 0 0 𝟔
0 0 0 0 −100 0 100 0 𝟕
[ 0 −133.33 0 0 0 0 0 −133.33] 𝟖
✓ Step-3: Determine the nodal displacement and determine the forces in the members of the
truss by applying finite element formulation
a) Nodal Displacement
[𝑲]{𝑫} = {𝑭}
151.2 38.40 −100 0 −51.20 −38.4 0 0 𝑢1 𝐹1
38.40 162.13 0 0 −38.40 −28.8 0 −133.33 𝑣1 𝐹2
−100 0 100 0 0 0 0 0 𝑢2 𝐹3
0 0 0 133.33 0 −133.33 0 0 𝑣2 𝐹4
𝟏𝟎𝟑 𝑢3 = 𝐹5
−51.2 −38.4 0 0 151.20 38.40 −100 0
−38.40 −28.8 0 −133.33 38.40 162.13 0 0 𝑣3 𝐹6
0 0 0 0 −100 0 100 0 𝑢 4 𝐹7
[ 0 −133.33 0 0 0 0 0 −133.33 ] { 𝑣4 } {𝐹8 }
By applying boundary conditions 𝑢1 = 𝑣1 = 𝑣2 = 𝑢4 = 𝑣4 = 0
Nodal Force Vector are
𝐹1 0 0
𝐹2 0 0
𝐹3 0 0
𝐹4 0 0
{𝐹} = = = 𝟏𝟎𝟑 𝑁
𝐹5 0 0
𝐹6 𝑃 −30
𝐹7 0 0
{𝐹8 } { 0 } { 0 }
Master of Structural Engineering Page 6
Finite Element Method in Structures Assignment -4
151.2 38.40 −100 0 −51.20 −38.4 0 0 0 0
38.40 162.13 0 0 −38.40 −28.8 0 −133.33 0 0
−100 0 𝟏𝟎𝟎 0 𝟎 𝟎 0 0 𝒖𝟐 0
𝟑 0 0 0 133.33 0 −133.33 0 0 0 𝟑 0
𝟏𝟎 = 𝟏𝟎
−51.2 −38.4 𝟎 0 𝟏𝟓𝟏. 𝟐𝟎 𝟑𝟖. 𝟒𝟎 −100 0 𝒖𝟑 0
−38.40 −28.8 𝟎 −133.33 𝟑𝟖. 𝟒𝟎 𝟏𝟔𝟐. 𝟏𝟑 0 0 𝒗𝟑 −30
0 0 0 0 −100 0 100 0 0 0
[ 0 −133.33 0 0 0 0 0 −133.33] { 0 } { 0 }
Reduced matrix equation
100 0 0 𝑢2 0
10 [ 03 𝑢 3
151.20 38.40 ] { 3 } = 10 { 0 }
0 38.40 162.13 𝑣3 −30
By using Cramer’s Rule
0 0 0
| 0 151.20 38.40 |
−30 38.40 162.13 = 0
∴ [𝑢2 ] = = 𝟎 𝒎𝒎
100 0 0 2,304,000
| 0 151.20 38.40 |
0 38.40 162.13
100 0 0
| 0 0 38.40 |
∴ [𝑢3 ] = 0 −30 162.13 = 115,200 = 𝟎. 𝟎𝟓 𝒎𝒎
100 0 0 2,304,000
| 0 151.20 38.40 |
0 38.40 162.13
100 0 0
| 0 151.20 0 |
∴ [𝑣3 ] = 0 38.40 −30 = −453,600 = −𝟎. 𝟏𝟗𝟔𝟖𝟕𝟓 𝒎𝒎
100 0 0 2,304,000
| 0 151.20 38.40 |
0 38.40 162.13
𝒖𝟏 𝟎
𝒗𝟏 𝟎
𝒖𝟐 𝟎
𝒖 𝒗𝟐 𝟎
∴ The nodal displacement is [ ] = 𝒖 = 𝒎𝒎
𝒗 𝟑 𝟎. 𝟎𝟓
𝒗𝟑 −𝟎. 𝟏𝟗𝟔𝟖𝟕𝟓
𝒖𝟒 𝟎
{𝒗𝟒 } { 𝟎 }
Master of Structural Engineering Page 7
Finite Element Method in Structures Assignment -4
b) Determine the forces in the members of the truss
1) Forces in the members on element-1
𝒖𝟏 𝟎
𝐸1 𝐴1 𝒗𝟏
𝐹1 = [−𝑙 −𝑚 𝑙 𝑚] { } = 100𝑥103 [−1 0 1 0] {𝟎} = 𝟎 𝒌𝑵
𝐿1 𝒖𝟐 𝟎
𝒗𝟐 𝟎
2) Forces in the members on element-2
𝑢2 0
𝐸2 𝐴2 𝑣2 0
𝐹2 = [−𝑙 −𝑚 𝑙 𝑚] {𝑢 } = 133.33𝑥103 [0 −1 0 1] { } = −𝟐𝟔. 𝟐𝟓 𝒌𝑵
𝐿2 3 0.05
𝑣3 −0.196875
3) Forces in the members on element-3
𝑢3 0.05
𝐸3 𝐴3 𝑣 −0.196875
𝐹3 = [−𝑙 −𝑚 𝑙 𝑚] {𝑢3 } = 100𝑥103 [1 0 −1 0] { } = 5 𝑘𝑁
𝐿3 4 0
𝑣4 0
4) Forces in the members on element-4
𝒖𝟒 0
𝐸4 𝐴4 𝒗𝟒 0
𝐹4 = [−𝑙 −𝑚 𝑙 𝑚] {𝒖 } = 133.33𝑥103 [0 −1 0 1] { } = 0 𝒌𝑵
𝐿4 𝟓 0
𝒗𝟓 0
5) Forces in the members on element-5
𝒖𝟓 0
𝐸5 𝐴5 𝒗𝟓 3 0
𝐹5 = [−𝑙 −𝑚 𝑙 𝑚] {𝒖 } = 80𝑥10 [−0.8 −0.6 0.8 0.6] { } = −𝟔. 𝟐𝟓 𝒌𝑵
𝐿5 𝟔 0.05
𝒗𝟔 −0.196875
c) Determine the reaction of the truss
Master of Structural Engineering Page 8
Finite Element Method in Structures Assignment -4
0
0
0
0
❖ 3[
𝑅1 = 10 151.2 38.40 −100 0 −51.20 −38.40 0 0] = 5 𝑘𝑁
0.05
−0.196875
0
{ 0 }
0
0
0
0
❖ 3
R 2 = 10 [38.40 162.133 0 0 −38.4 −28.8 0 −133.33] = 3.75 kN
0.05
−0.196875
0
{ 0 }
0
0
0
0
❖ R 3 = 103 [0 0 0 −133.33 0 −133.33 0 0] = 26.50 kN
0.05
−0.196875
0
{ 0 }
0
0
0
0
❖ R 4 = 103 [0 0 0 0 −100 0 100 0] = −5 kN
0.05
−0.196875
0
{ 0 }
0
0
0
0
❖ R 5 = 103 [0 −133.33 0 0 0 0 0 −133.33] = 0 kN
0.05
−0.196875
0
{ 0 }
Figure 2.2: Member force and reaction on the truss
Master of Structural Engineering Page 9
Finite Element Method in Structures Assignment -4
#Q-3: Determine the deflection, rotations, and reactions force of beam with spring
system shown in Figure-3 Given E = 210 GPa, k= 200 kN/m, L =3m, I= 2x10-4 m4
and P = 50 kN
Figure-3
Solution:
✓ Step-1: Assign the nodal point and the naming member (element) for each structure, the
beam has a roller (hinge) support at node 2 and a spring support at node 3. To solve this
problem, we use two beam element and one spring element .it has four nodes and neglect
axial deformation. Node-Element connectivity by table:
Element Node-1 Node-2
1 1 2
2 2 3
3 3 4
Master of Structural Engineering Page 10
Finite Element Method in Structures Assignment -4
✓ Step-2: Compute element stiffens matrix for each element and assemble Global stiffness
matrix
Element(member) -1 Element-3/spring/
V1 1 V2 2 V3 V4
12 6𝐿 −12 6𝐿 V1 Let
𝒌𝑳𝟑
[ 𝑘𝑠 ] = [ 𝑘 −𝑘
𝐸𝐼 6𝐿 4𝐿2 −6𝐿 2𝐿2 1 V3
[𝑘1 ] = 3 [ ] ] 𝒌′ =
𝐿 −12 −6𝐿 12 −6𝐿 V2 𝑬𝑰
6𝐿 2𝐿2 −6𝐿 4𝐿2 2
−𝑘 𝑘 V4
Element(member) -2 𝑬𝑰 𝑘′ −𝑘′]
[𝑘𝑠 ] = [
𝑳𝟑 −𝑘′ 𝑘′
V2 2 V3 3
12 6𝐿 −12 6𝐿 V2
𝐸𝐼 6𝐿 4𝐿2
−6𝐿 2𝐿2 2 𝒌𝑳𝟑 𝟐𝟎𝟎∗(𝟑𝟎𝟎𝟎)𝟑
[𝑘2 ] = 3 [ ] 𝒌′ = = 𝟐𝟏𝟎∗𝟏𝟎𝟑𝑿𝟐∗𝟏𝟎𝟖 = 𝟎. 𝟏𝟐𝟖𝟓𝟕𝟏𝟒𝟐𝟗
𝐿 −12 −6𝐿 12 −6𝐿 V3 𝑬𝑰
6𝐿 2𝐿2 −6𝐿 4𝐿2 3
Assemble Global stiffens matrix
[𝑲] = [𝒌]𝟏 + [𝒌]𝟐 + [𝒌]𝟑
𝑽𝟏 𝜽𝟏 𝑽𝟐 𝜽𝟐 𝑽𝟑 𝜽𝟑 𝑽𝟒
12 6L −12 6𝐿 0 0 0 𝑽𝟏
6𝐿 4𝐿2 −6𝐿 2𝐿2 0 0 0 𝜽𝟏
𝑬𝑰 −12 −6𝐿 24 0 −12 6𝐿 0 𝑽𝟐
[𝑲] = 𝟑 6𝐿 2𝐿2 0 8𝐿2 −6𝐿 2𝐿2 0 𝜽𝟐
𝑳 0 0 −12 −6𝐿 12 + 𝒌′ −6𝐿 −𝒌′ 𝑽𝟑
0 0 6𝐿 2𝐿2 −6𝐿 4𝐿2 0 𝜽𝟑
[ 0 0 0 0 −𝒌′ 0 −𝒌′] 𝑽𝟒
✓ Step-3: Determine the nodal displacement(deflection and rotation) and support reactions by
applying finite element formulation
i. Nodal Displacement(deflection and rotation)
[𝐾]{𝐷} = {𝐹}
Master of Structural Engineering Page 11
Finite Element Method in Structures Assignment -4
12 6L −12 6𝐿 0 0 0 𝑽𝟏 𝑭𝟏𝒚
6𝐿 4𝐿2 −6𝐿 2𝐿2 0 0 0 𝜽𝟏 𝑴𝟏
𝑽𝟐 𝑭
𝑬𝑰 −12 −6𝐿 24 0 −12 6𝐿 0 𝟐𝒚
6𝐿 2𝐿2 0 8𝐿2 −6𝐿 2𝐿2 0 𝜽𝟐 = 𝑴𝟐
𝑳𝟑 𝑭𝟑𝒚
0 0 −12 −6𝐿 12 + 𝒌′ −6𝐿 −𝒌′ 𝑽𝟑
0 0 6𝐿 2𝐿2 −6𝐿 4𝐿2 0 𝜽𝟑 𝑴𝟑
[ 0 0 0 0 −𝒌′ 0 −𝒌′] {𝑽𝟒 } {𝑭𝟒𝒚 }
We now by applying boundary conditions 𝑽𝟏 = 𝜽𝟏 = 𝑽𝟐 = 𝑽𝟒 = 𝟎
Nodal Force Vector are
𝑭𝟏𝒚 0 0
𝑴𝟏 0 0
𝑭𝟐𝒚 0 0
{𝐹} = 𝑴𝟐 = 0 = 0 𝑘𝑁
𝑭𝟑𝒚 −𝑃 −50
𝑴𝟑 0 0
{𝑭𝟒𝒚 } { 0 } { 0 }
12 6L −12 6𝐿 0 0 0 0 0
6𝐿 4𝐿2 −6𝐿 2𝐿2 0 0 0 0 0
𝑬𝑰 −12 −6𝐿 24 0 −12 6𝐿 0 0 0
6𝐿 2𝐿2 0 𝟖𝑳𝟐 −𝟔𝑳 𝟐𝑳𝟐 0 𝜽𝟐 = 𝟎
𝑳𝟑
0 0 −12 −𝟔𝑳 𝟏𝟐 + 𝒌′ −𝟔𝑳 −𝒌′ 𝑽𝟑 −𝟓𝟎
0 0 6𝐿 𝟐𝑳𝟐 −𝟔𝑳 𝟒𝑳𝟐 0 𝜽𝟑 𝟎
[ 0 0 0 0 −𝒌′ 0 −𝒌′ ] { 0 } { 0 }
✓ Deleting the first three and seventh equation (rows and columns), we have the following
reduced equation,
𝟐
𝟐𝑳𝟐 𝜽𝟐 0
𝑬𝑰 𝟖𝑳 −𝟔𝑳
𝟑
[ −𝟔𝑳 𝟏𝟐 + 𝒌′ −𝟔𝑳 ] { 𝑉 3 } = 𝟏𝟎 { −50 } 𝑁
𝑳𝟑
𝟐𝑳𝟐
−𝟔𝑳 𝟒𝑳𝟐 𝜽𝟑 0
𝟐𝟏 ∗ 𝟏𝟎𝟒 ∗ 𝟐 ∗ 𝟏𝟎𝟖 8 ∗ 3000
2
−6 ∗ 3000 2 ∗ 30002 𝜽𝟐 0
𝑉 𝟑
[−6 ∗ 3000 12 + 0.128571429 −6 ∗ 3000 ] { 3 } = 𝟏𝟎 { −50} 𝑁
𝟑𝟎𝟎𝟎𝟑
2 ∗ 30002 −6 ∗ 3000 4 ∗ 3000 2 𝜽𝟑 0
7.2 ∗ 107 −1.8 ∗ 104 1.8 ∗ 107 𝜃2 0
1555.556 [−1.8 ∗ 104 4 ] { 𝑉 } = 103 {−50} 𝑁
12.128571429 −1.8 ∗ 10 3
1.8 ∗ 107 −1.8 ∗ 104 3.6 ∗ 107 𝜃3 0
𝟏. 𝟏𝟐 ∗ 𝟏𝟎𝟖 −𝟐. 𝟖 ∗ 𝟏𝟎𝟒 𝟐. 𝟖 ∗ 𝟏𝟎𝟕 𝜽𝟐 0
[−𝟐. 𝟖 ∗ 𝟏𝟎𝟒 𝟒 ] { 𝑉 } = {−50} 𝑁
𝟏. 𝟖𝟖𝟔𝟔𝟕 −𝟐. 𝟖 ∗ 𝟏𝟎 3
𝟐. 𝟖 ∗ 𝟏𝟎𝟕 −𝟐. 𝟖 ∗ 𝟏𝟎𝟒 𝟓. 𝟔 ∗ 𝟏𝟎 𝟕 𝜽 𝟑 0
Master of Structural Engineering Page 12
Finite Element Method in Structures Assignment -4
By using Cramer’s Rule
0 −2.8 ∗ 104 2.8 ∗ 107
|−50 1.88667 −2.8 ∗ 104 |
0 −2.8 ∗ 104 5.6 ∗ 107 −3.92 ∗ 1013
∴ [𝜽𝟐 ] = = = −𝟎. 𝟎𝟎𝟐𝟒𝟗 𝒓𝒂𝒅
1.12 ∗ 108 −2.8 ∗ 104 2.8 ∗ 107 1.57323 ∗ 1016
|−2.8 ∗ 104 1.88667 −2.8 ∗ 104 |
2.8 ∗ 107 −2.8 ∗ 104 5.6 ∗ 107
1.12 ∗ 108 0 2.8 ∗ 107
|−2.8 ∗ 10 −50 −2.8 ∗ 104 |
4
2.8 ∗ 107 0 5.6 ∗ 107 −2.744 ∗ 1017
∴ [𝑉3 ] = = = −17.442 𝑚𝑚 = −𝟎. 𝟎𝟏𝟕𝟒𝟒𝟐 𝒎
1.12 ∗ 108 −2.8 ∗ 104 2.8 ∗ 107 1.5732 ∗ 1016
|−2.8 ∗ 104 1.88667 −2.8 ∗ 104 |
2.8 ∗ 10 7 −2.8 ∗ 104 5.6 ∗ 107
1.12 ∗ 108 −2.8 ∗ 104 0
|−2.8 ∗ 104 1.88667 −50|
2.8 ∗ 107
−2.8 ∗ 104
0 −1.176 ∗ 1014
[𝜽 ]
∴ 𝟑 = = = −𝟎. 𝟎𝟎𝟕𝟒𝟕 𝒓𝒂𝒅
1.12 ∗ 108 −2.8 ∗ 104 2.8 ∗ 107 1.57323 ∗ 1016
|−2.8 ∗ 104 1.88667 −2.8 ∗ 104 |
2.8 ∗ 107 −2.8 ∗ 104 5.6 ∗ 107
𝑽𝟏 𝟎 𝒎𝒎 𝟎 𝒎
𝜽𝟏 𝟎 𝒓𝒂𝒅 𝟎 𝒓𝒂𝒅
𝑽𝟐 𝟎 𝒎𝒎 𝟎 𝒎
𝑽
∴ 𝐓𝐡𝐞 𝐧𝐨𝐝𝐚𝐥 𝐝𝐢𝐬𝐩𝐥𝐚𝐜𝐞𝐦𝐞𝐧𝐭 𝐢𝐬 = [𝑫] = [ ] = 𝜽𝟐 = −𝟎. 𝟎𝟎𝟐𝟒𝟗 𝒓𝒂𝒅 = −𝟎. 𝟎𝟎𝟐𝟒𝟗 𝒓𝒂𝒅
𝜽 𝑽𝟑 𝒎𝒎 𝒎
−𝟏𝟕. 𝟒𝟒𝟐 −𝟎. 𝟎𝟏𝟕𝟒𝟒𝟐
𝜽𝟑 −𝟎. 𝟎𝟎𝟕𝟒𝟕 𝒓𝒂𝒅 −𝟎. 𝟎𝟎𝟕𝟒𝟕 𝒓𝒂𝒅
{𝑽𝟒 } { 𝟎 } 𝒎𝒎 { 𝟎 } 𝒎
ii. Reactions Forces
From the global FE equation, we obtain the nodal reaction forces as,
[𝐾]{𝐷} = {𝐹}
12 6L −12 6𝐿 0 0 0 𝟎 𝑭𝟏𝒚
6𝐿 4𝐿2 −6𝐿 2𝐿2 0 0 0 𝟎 𝑴𝟏
𝑭𝟐𝒚
𝑬𝑰 −12 −6𝐿 24 0 −12 6𝐿 0 𝟎
6𝐿 2𝐿2 0 8𝐿2 −6𝐿 2𝐿2 0 −𝟎. 𝟎𝟎𝟐𝟒𝟗 = 𝑴𝟐
𝑳𝟑 𝑭𝟑𝒚
0 0 −12 −6𝐿 12 + 𝒌′ −6𝐿 −𝒌′ −𝟏𝟕. 𝟒𝟒𝟐
0 0 6𝐿 2𝐿2 −6𝐿 4𝐿2 0 −𝟎. 𝟎𝟎𝟕𝟒𝟕 𝑴𝟑
[ 0 0 0 0 −𝒌′ 0 −𝒌′ ] { 𝟎 } {𝑭𝟒𝒚 }
From this we find out rection forces
Master of Structural Engineering Page 13
Finite Element Method in Structures Assignment -4
0
0
0
EI 𝟐𝟏∗𝟏𝟎𝟒 ∗𝟐∗𝟏𝟎𝟖
❖ 𝑭𝟏𝒚 = 3 [12 6L −12 6L 0 0 0] −0.00249 = ∗ 6 ∗ 3000 ∗ −0.00249 = −𝟔𝟗. 𝟗𝒌𝑵
L 𝟑𝟎𝟎𝟎𝟑
−17.442
−0.00747
{ 0 }
0
0
0
EI 𝟐𝟏∗𝟏𝟎𝟒 ∗𝟐∗𝟏𝟎𝟖
❖ M1 = 3 [6L 4L2 −6L 2L2 0 0 0] −0.00249 = ∗ 2 ∗ 30002 ∗ −0.00249 = −𝟔𝟗. 𝟕 𝐤𝐍𝐦
L 𝟑𝟎𝟎𝟎𝟑
−17.442
−0.00747
{ 0 }
0
0
0
𝑬𝑰 𝟐𝟏∗𝟏𝟎𝟒 ∗𝟐∗𝟏𝟎𝟖
❖ F2y =
𝑳𝟑
[−12 −6L 24 0 −12 6L 0] −0.00249 = 𝟑𝟎𝟎𝟎𝟑
(−12 ∗ −17.442 + 6 ∗ 3000 ∗ −0.00747) = 𝟏𝟏𝟔. 𝟒 𝐤𝐍
−17.442
−0.00747
{ 0 }
0
0
0
𝑬𝑰 𝟐𝟏∗𝟏𝟎𝟒 ∗𝟐∗𝟏𝟎𝟖
❖ M2 = [6L 2L2 0 8L2 −6L 2L2 0] −0.00249 = (8 ∗ 30002 ∗ −0.00249 − 6 ∗ 3000 ∗
𝑳𝟑 𝟑𝟎𝟎𝟎𝟑
−17.442
−0.00747
{ 0 }
−17.442 + 2 ∗ 30002 ∗ −0.00747) = 𝟎𝒌𝑵𝒎
0
0
0
𝑬𝑰 𝟐𝟏∗𝟏𝟎𝟒 ∗𝟐∗𝟏𝟎𝟖
❖ 𝑭𝟑𝒚 = [0 0 −12 −6L 12 + k′ −6L k′] −0.00249 = (−6 ∗ 3000 ∗ −0.00249 + 12.13 ∗
𝑳𝟑 𝟑𝟎𝟎𝟎𝟑
−17.442
−0.00747
{ 0 }
−17.442 − 6 ∗ 3000 ∗ −0.00747) = −𝟓𝟎𝒌𝑵
0
0
0
𝑬𝑰 𝟐𝟏∗𝟏𝟎𝟒 ∗𝟐∗𝟏𝟎𝟖
❖ M3 = [𝟎 𝟎 𝟔𝑳 2L2 −𝟔𝑳 4L2 𝟎] −0.00249 = (2 ∗ 30002 ∗ −0.00249 − 6 ∗ 3000 ∗
𝑳𝟑 𝟑𝟎𝟎𝟎𝟑
−17.442
−0.00747
{ 0 }
−17.442 + 4 ∗ 30002 ∗ −0.00747) = 𝟎 𝒌𝑵𝒎
0
0
0
𝑬𝑰 21∗104 ∗2∗108
❖ 𝑭𝟒𝒚 = [0 0 0 0 −𝑘′ 0 −𝑘′] −0.00249 = (−0.128𝟔 ∗ −17.442) = 𝟑. 𝟓𝟎𝒌𝑵
𝑳𝟑 30003
−17.442
−0.00747
{ 0 }
Master of Structural Engineering Page 14
Finite Element Method in Structures Assignment -4
𝐅𝟏𝐲
−𝟔𝟗. 𝟗 𝐤𝐍
𝐌𝟏 −𝟔𝟗. 𝟕 𝐤𝐍𝐦
𝐅𝟐𝐲 𝟏𝟏𝟔. 𝟒 𝐤𝐍
∴ 𝐓𝐡𝐞 𝐬𝐮𝐩𝐩𝐨𝐫𝐭 𝐫𝐞𝐚𝐜𝐭𝐢𝐨𝐧 𝐢𝐬 𝐌𝟐 = 𝟎 𝐤𝐍𝐦
𝐅𝟑𝐲 −𝟓𝟎 𝐤𝐍
𝐌𝟑 𝟎 𝐤𝐍𝐦
𝐅 { 𝟑. 𝟓𝟎 } 𝐤𝐍
{ 𝟒𝐲 }
Master of Structural Engineering Page 15
