FEM Assignment 1
FEM Assignment 1
1. Solve the deflection and stresses in the given truss manually by FEM, if both angles are 60 o and
the length of each member is 4m. The area is A=20cm2 and E=2.1x1011pa. A vertical load of 80kN
and horizontal force 60kN are applied as shown and the supports at the left and right are hinged
and roller as shown.
3 2
o o
60 60
1
1 2
f1y, v1 f2y, v2
f1x, u1 f2x, u2
𝑐2 𝑠𝑐 −𝑐 2 − 𝑠𝑐
2
[𝐾]𝑒 = 𝑘 [ 𝑠𝑐2 𝑠2 −𝑠𝑐 −𝑠 ]
−𝑐 −𝑠𝑐 𝑐 2 𝑠𝑐
2
−𝑠𝑐 −𝑠 2 𝑠𝑐 𝑠
U1 V1 U2 V2
1 0 −1 0 U1
0 0 0 0 V1
[𝐾](1) = 𝑘[ ] ………………. (1)
−1 0 1 0 U2
0 0 0 0 V2
For element (2) =120°
f3y, v3
edof(5,6) f3x, u3
Ө =120°
edof(3,4) f2x, u2
f2y, v2
cos2 120 sin 120 cos 120 −cos2 120 −sin 120 cos 120
[𝐾]𝑒 = 𝑘 [ sin 120 2cos 120 sin2 120 − sin 120 cos 120 − sin2 120 ]
−cos 120 −sin 120 cos 120 cos2 12 0 sin 120 cos 120
−sin 120 cos 120 − sin2 120 sin 120 cos 120 sin2 120
U2V2U3V3
F3y, v3
F3x, u3
f1y, v1
f1x, u1
U1V1U3V3
1 0 −10 0 0 U1
0 0 00 0 0 V1
−1 0 10 0 0 U2
[𝐾](1) =𝑘
0 0 00 0 0 V2
0 0 00 0 0 0
[ 0 0 00 0 0] 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0.25 −0.433 −0.25 0.433 U2
[𝐾](2) =𝑘
0 0 −0.433 0.75 0.433 −0.75 V2
0 0 −0.25 0.433 0.25 −0.433 U3
[0 0 0.433 −0.75 −0.433 0.75 ] V3
U1V1U2V2U3V3
𝑓3𝑥 60,000
= −0.25𝑈2 + 0.5𝑈3 + 0, f3x=60 kN10.5∗107 = 0.00057143 = −0.25𝑈2 +
𝑘
0.5𝑈3 … … … … … … … … (𝑏)
𝑓3𝑦
= 0.433𝑈2 + 0𝑈3 + 1.5𝑉3, f3y=80 kN
𝑘
−80,000
= −0.000761904 = 0.433𝑈2 +0+1.5𝑉3…………...(c)
10.5∗107
0=𝑈2 − 0.000505651
𝑈2 = 0.000505651𝑚
𝑈2 = 0.505651𝑚𝑚
= 0.505651𝑚𝑚…………………………………………………..………Answer
=0.001395685m=1.395685mm…………………………………………….. Answer
=0.00065391m =-0.65391mm…………………………………………Answer
𝑈3=1.395685mm
𝑉3=-0.65391mm
For element(1)
∈𝐴 ∈𝐴
𝑓2 = (𝑢2− 𝑢1 ) = (𝑈2− 𝑈1 ) , 𝑈2 = 𝑢2
𝐿 𝐿
𝑈1 = 𝑢1 = 0
𝑓2 ∈
𝜎1 = = (𝑈 𝑈 )
𝐴 𝐿 2− 1
∈𝑈2 1011 𝑁 0.000505651𝑚 𝑁
𝜎1 = =2.1* ∗ =26546677.5
𝐿 𝑚2 4𝑚 𝑚2
=26.5466775*106 𝑁/𝑚2 ………………………………………. Answer
Forelement(2)
𝑓3 = 𝑘(𝑢2 −𝑢3 )
∈𝐴
𝑓3 = (−𝑈3 cos 60° + 𝑉3 sin 60° +𝑈2 cos 60° )
𝐿
𝑓3 𝐸
𝜎2 = = (cos 60°(𝑈2 − 𝑈3 ) +
𝐴 𝐿
2.1*1011 N/m2
𝑉3 sin 60° )= [cos 60°(0.000505651-0.001395685)+(-0.00065391) sin 60°]
4m
= -53094142.5N/𝑚2
𝑓3 −𝑓1 −∈
𝜎3 = = = (−𝑈3 cos 60° − 𝑉3 sin 60° )
𝐴 𝐴 𝐿
−2.1∗1011 𝑁/𝑚2
= [− 0.001395685 cos 60° − (−0.00065391) sin 60°]𝑚
4𝑚
= 6905840.981N/𝑚2