FEM Assignment 4 & 5
FEM Assignment 4 & 5
𝐹1𝑦 12 6𝐿 − 12 6𝐿 𝑑1𝑦
𝑚 2 2 Ф1
𝐸𝐼 6𝐿 4𝐿 − 6𝐿 2𝐿
[ 1 ] = 𝐿3 [ ]
𝐹𝑦 −12 − 6𝐿 12 − 6𝐿 𝑑2𝑦
𝑚2 6𝐿 2𝐿2 − 6𝐿 4𝐿2 { Ф2 }
12 3𝐿 − 12 3𝐿 𝑑2𝑦 𝐹2𝑦
2 2
8 𝐸 𝐼 3𝐿 𝐿 − 3𝐿 𝐿 /2 Ф2 𝑚
[ ] = [ 2 ] -------- 2
𝐿3 −12 − 3𝐿 12 − 3𝐿 𝑑3𝑦 𝐹3𝑦
3𝐿 𝐿2 /2 − 3𝐿 𝐿2 { Ф3 } 𝑚3
Delete 1st row and 1st column, 2nd row and 2nd column, 5th row and 5th column in the
stiffness matrix,
24 0 3𝐿
𝐿2 𝑑2𝑦 −𝑃
2
8𝐸𝐼 0 2𝐿 { Ф2 } = { 0 }
2
𝐿3
𝐿2 Ф3 0
[ 3𝐿 𝐿2 ]
2
Solving,
− 𝟕 𝑷 𝑳𝟑
The mid span deflection, 𝒅𝟐𝒚 =
𝟕𝟔𝟖 𝑬 𝑰
𝟏
3. Evaluate ∫−𝟏(𝒙𝟒 + 𝒙𝟐 ) 𝒅𝒙 by applying 3-point Gaussian quadrature.
Solution:
1
Given: Integral, I = ∫−1(𝑥 4 + 𝑥 2 ) 𝑑𝑥
f(x) = 𝑥 4 + 𝑥 2
3
𝑥1 = √5 = 0.774596669;
𝑥2 = 0;
3
𝑥3 = - √5 = - 0.774596669
5
𝑤1 = = 0.555555;
9
8
𝑤2 = = 0.888888;
9
5
𝑤3 = 9 = 0.555555
f (𝑥1 ) = 𝑥1 4 + 𝑥1 2
f (𝑥2 ) = 𝑥2 4 + 𝑥2 2
= (0)4 + (0)2 = 0
f (𝑥3 ) = 𝑥3 4 + 𝑥3 2
Solution:
Given: Natural co-ordinates of point P,
ε = 0.6, η = 0.8
Cartesian co-ordinates of point 1, 2, 3, and 4
𝑥1 = 3; 𝑦1 = 2; 𝑥2 = 9; 𝑦2 = 4; 𝑥3 = 6; 𝑦3 = 8; 𝑥4 = 4; 𝑦4 = 5
The tetrahedron element, with three translational degrees of freedom per node, is shown
in the global xyz coordinate system in Figure. (the global coordinates are denoted as
x, y, z instead of X, Y, Z, for simplicity). Since there are 12 nodal degrees of freedom
Q3i-2, Q3i-1, Q3i, Q3j-2 ...,Q3l and three displacement components u, v, and w, we choose
the displacement variation to be linear as
and similar expressions for v(x, y, z) and w(x, y, z). Thus, the displacement field can be
expressed in matrix form as
and similar expressions for v(x, y, z) and w(x, y, z). Thus, the displacement field can be
expressed in matrix form as
Since the matrices [B] and [D] are independent of x, y, and z, the stiffness matrix can be
obtained by carrying out matrix multiplications as
6. Derive the Shape functions for Constant Strain Triangular (CST) element.
𝟏
7. Evaluate the Jacobian matrix [J] at ε = η = for the linear quadrilateral element
𝟐
shown in figure.
Solution:
1 1
Given: Natural coordinates at point, P coordinates ε = 2 = 0.5 and η = = 0.5
2
Cartesian co-ordinates of point 1, 2, 3, and 4
𝑥1 = 4; 𝑦1 = 4; 𝑥2 = 7; 𝑦2 = 5; 𝑥3 = 8; 𝑦3 = 10; 𝑥4 = 3; 𝑦4 = 8
2.25 0.875
Result: Jacobian Matrix, [ J ] = [ ] Ans.
0.25 2.375
8. Calculate the temperature distribution in a one-dimension fin with physical
properties given in figure. The fin is rectangular in shape and is 120 mm long, 40
mm wide and 10 mm thick. Assume that convection heat loss occurs from the end
of the fin. Use two elements. Take k = 0.3 W/mm 0C; h = 1 x 𝟏𝟎−𝟑 W/mm2 0C; 𝑻𝜶
= 200 C
9. Explain the importance of Eigen value and Eigen vectors.
Answer: Importance of Eigenvalues and Eigenvectors:
Eigenvalues and eigenvectors play a crucial role in Finite Element Methods, especially
in problems involving dynamic analysis, stability, and vibration analysis:
1. Vibration Analysis:
o In structural dynamics, eigenvalues represent the natural frequencies of a
system, which are critical for understanding its response to vibrations.
o Eigenvectors correspond to the mode shapes that describe how the
structure deforms at each natural frequency.
2. Dynamic Analysis:
o They help determine the critical frequencies to avoid resonance, ensuring
the safety and reliability of structures and mechanical systems.
3. Stability Analysis:
o Eigenvalues are used to analyze buckling problems in structures. The
smallest eigenvalue indicates the critical load at which buckling occurs.
4. Simplifying Complex Systems:
o Eigenvectors form a basis to decouple complex systems of equations into
independent equations, simplifying computations in dynamic and stability
analysis.
5. Energy Methods:
o The eigenvalues relate to the energy levels in the system and aid in
understanding the potential and kinetic energy distribution.
Thus, eigenvalues and eigenvectors are indispensable in ensuring that FEM solutions
accurately predict the behavior of engineering systems under various conditions.
10. Derive the stiffness matrix for one-dimensional Heat conduction element.
11. Explain the following:
(i) Consistent mass matrix, (ii)Lumped mass matrix
Answer:
i) Consistent mass matrix for bar element:
mass matrix, [m] = ∫ 𝜌 [𝑁]𝑇 [N] dv
𝑙 𝑁1
= ρ A ∫0 { } [𝑁1 𝑁2 ] dx
𝑁2
𝑥 𝑥
Substituting, 𝑁1 = 1 − ; 𝑎𝑛𝑑 𝑁2 = 𝑙
𝑙
𝑥
𝑙 1− 𝑥 𝑥
𝑙
[m] = ρ A ∫0 { 𝑥 } [(1 − 𝑙 ) ( 𝑙 )] dx
𝑙
𝑥 𝑥 𝑥2
𝑙 (1 − )2 −
𝑙 𝑙 𝑙2
= ρ A ∫0 [ 𝑥 𝑥2 𝑥2
] dx
−
𝑙 𝑙2 𝑙2