Practice Assignment 8
Finite Element Methods
ME G512
1. Determine the eigenvalues and eigenvectors for the stepped bar shown in Fig. 1. Also,
plot the first and second modes of vibration. (E = 3 × 107 psi, Density, ρ = 7.324 × 10-4
lb/in3).
2. Use two finite elements based on Euler-Bernoulli Beam theory and evaluate the natural
frequency for the beam shown in Fig. 2. Plot the first and second modes of vibration.
Material properties: E = 200GPa, ρ = 7840 kg/m3, I = 2000 mm4, A = 240 mm2.
Figure 1 Figure 2
Problem 1:
Global Stiffness Matrix :
KG =
3000000 -3000000 0
-3000000 6000000 -3000000
0 -3000000 3000000
Global Mass Matrix :
MG =
2.4413e-03 1.2207e-03 0
1.2207e-03 3.0517e-03 3.0517e-04
0 3.0517e-04 6.1033e-04
Practice Assignment 8
Finite Element Methods
ME G512
eigenvalues =
6.6838e+08
7.6101e+09
eigenvectors =
-1.4527e+01 -1.1572e+01
-1.7956e+01 3.7447e+01
Natural Frequencies in Hz are :
freq =
4.1146e+03
1.3884e+04
firstmodeshape =
0
-1.4527e+01
-1.7956e+01
Practice Assignment 8
Finite Element Methods
ME G512
secondmodeshape =
0
-1.1572e+01
3.7447e+01
Problem 2:
Global Stiffness Matrix :
KG =
1.7778e+05 2.6667e+04 -1.7778e+05 2.6667e+04 0 0
2.6667e+04 5.3333e+03 -2.6667e+04 2.6667e+03 0 0
-1.7778e+05 -2.6667e+04 3.5556e+05 0 -1.7778e+05 2.6667e+04
2.6667e+04 2.6667e+03 0 1.0667e+04 -2.6667e+04 2.6667e+03
0 0 -1.7778e+05 -2.6667e+04 1.7778e+05 -2.6667e+04
0 0 2.6667e+04 2.6667e+03 -2.6667e+04 5.3333e+03
Global Mass Matrix :
MG =
2.0966e-01 8.8704e-03 7.2576e-02 -5.2416e-03 0 0
8.8704e-03 4.8384e-04 5.2416e-03 -3.6288e-04 0 0
7.2576e-02 5.2416e-03 4.1933e-01 0 7.2576e-02 -5.2416e-03
-5.2416e-03 -3.6288e-04 0 9.6768e-04 5.2416e-03 -3.6288e-04
0 0 7.2576e-02 5.2416e-03 2.0966e-01 -8.8704e-03
Practice Assignment 8
Finite Element Methods
ME G512
0 0 -5.2416e-03 -3.6288e-04 -8.8704e-03 4.8384e-04
eigenvalues =
2.0298e+04
8.0998e+05
9.2655e+06
7.8053e+07
eigenvectors =
6.3970e-01 1.3720e+00 2.1505e-01 8.9887e-01
3.6522e+00 -1.3763e+00 -2.6945e+01 3.0793e+01
1.8841e+00 -1.9008e+00 2.1141e+00 3.5501e+00
4.3227e+00 -1.5253e+01 3.3981e+01 1.1437e+02
Natural Frequencies in Hz are :
freq =
2.2675e+01
1.4324e+02
4.8446e+02
1.4061e+03
firstmodeshape =
0
6.3970e-01
1.8841e+00
Practice Assignment 8
Finite Element Methods
ME G512
secondmodeshape =
0
1.3720e+00
-1.9008e+00