Vibration Analysis Exercise 4 6 Report
Vibration Analysis Exercise 4 6 Report
2. Methodology
To solve the system, the following steps were taken:
1. Construct the mass matrix M and stiffness matrix K.
2. Use MATLAB's generalized eigenvalue solver to compute the eigenvalues and
eigenvectors of the system.
3. The square roots of the eigenvalues give the natural frequencies in rad/s. Dividing by 2π
gives the frequencies in Hz.
4. The eigenvectors represent the mode shapes of the system.
5. Normalize the mode shapes for interpretation and plotting.
3. MATLAB Code
The following MATLAB code was used to perform the analysis:
m1 = 5;
m2 = 10;
k1 = 4000;
k2 = 6000;
M = [m1, 0;
0, m2];
omega_squared = diag(lambda);
omega_n = sqrt(omega_squared);
f_n = omega_n / (2 * pi);
figure;
subplot(1,2,1);
bar(mode_shapes(:,1));
title('Mode Shape 1');
xlabel('Mass Index'); ylabel('Amplitude');
grid on;
subplot(1,2,2);
bar(mode_shapes(:,2));
title('Mode Shape 2');
xlabel('Mass Index'); ylabel('Amplitude');
grid on;
4. Results
The natural frequencies and normalized mode shapes obtained from MATLAB are as
follows:
[ 0.81 -0.55 ]
[ 1.00 1.00 ]
5. Conclusion
The analysis successfully determined the two natural frequencies of the 2-DOF system. The
first mode represents both masses moving in-phase while the second mode corresponds to
out-of-phase motion. This report provides a MATLAB-based approach for modal analysis of
a simple mechanical system.