Scribd
Upload
70 views6 pages

Assignment 2 Modal Analysis of Mechanical Systems: Submitted By, Pranav R ME15M081

The document describes performing a modal analysis of a mechanical system using MATLAB. It defines mass and stiffness matrices for a 6 degree of freedom system and calculates the natural frequencies and mode shapes. It then plots the frequency response functions between each degree of freedom considering 6 modes and 3 modes. Finally, it plots the time response of each degree of freedom to visualize the system response.

Uploaded by

Pranav Raveendran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
70 views6 pages

Assignment 2 Modal Analysis of Mechanical Systems: Submitted By, Pranav R ME15M081

The document describes performing a modal analysis of a mechanical system using MATLAB. It defines mass and stiffness matrices for a 6 degree of freedom system and calculates the natural frequencies and mode shapes. It then plots the frequency response functions between each degree of freedom considering 6 modes and 3 modes. Finally, it plots the time response of each degree of freedom to visualize the system response.

Uploaded by

Pranav Raveendran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 6

ASSIGNMENT 2

MODAL ANALYSIS OF MECHANICAL SYSTEMS

Submitted by,

PRANAV R

ME15M081
MATLAB CODE
clear all
clc
m1=5.4e5;
m2=3.6e5;
m3=1.8e5;
m4=1.2e5;
m5=1e5;
m6=0.8e5;
M= [5.4e5 0 0 0 0 0; 0 3.6e5 0 0 0 0; 0 0 1.8e5 0 0 0; 0 0 0 1.2e5 0 0; 0 0 0
0 1e5 0; 0 0 0 0 0 0.8e5];
k1=430e6;
k2=290e6;
k3=200e6;
k4=160e6;
k5=120e6;
k6=80e6;
K= [k1+k2 -k2 0 0 0 0;-k2 k2+k3 -k3 0 0 0;0 -k3 k3+k4 -k4 0 0;0 0 -k4 k4+k5 -
k5 0;0 0 0 -k5 k5+k6 -k6; 0 0 0 0 -k6 k6];
xeta= [0.01 0.015 0.02 0.025 0.03 0.035];
F= [20000 0 0 0 0 0];
syms f;
[si, eigenvalue]=eig(K,M);
vm = [0 0 0 0 0 0;si];
m=sqrt(eigenvalue);
naturalfreq=diag(m)
phi= [si(:,1) si(:,2) si(:,3) si(:,4) si(:,5) si(:,6)];
for i=1:6
damped_natural_freq(i,1)=naturalfreq(i,1)*sqrt(1-xeta(1,i)^2);
end
damped_natural_freq
for i=1:6
x=vm(:,i);
y=[0;1;2;3;4;5;6];
figure(1)
subplot(1,6,i);
plot(x,y,'b','Marker','x','linewidth',1.5);
heading= ['Mode',num2str(i)];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
grid on
ylabel('DOF Number','fontsize',10)
xlabel('Amp1itude','fontsize',10)
end
for j=1:3
for k=1:3
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3))+phi(k,4)*phi(j,4)./(naturalfreq(4)^2-
w.^2+1i*2*xeta(4)*w*naturalfreq(4))+phi(k,5)*phi(j,5)./(naturalfreq(5)^2-
w.^2+1i*2*xeta(5)*w*naturalfreq(5))+phi(k,6)*phi(j,6)./(naturalfreq(6)^2-
w.^2+1i*2*xeta(6)*w*naturalfreq(6));
figure(2)
subplot(3,3,j+3*(k-1));
plot(f,20*log10(abs(alpha)));
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('freq Hz');
ylabel('\alpha dB');
hold on;
end
end
for j=1:3
for k=1:3
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3));
figure(2)
subplot(3,3,j+3*(k-1));
plot(f,20*log10(abs(alpha)),'r');
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('Freq Hz');
ylabel('\alpha dB');
hold on;
end
end
legend ('considering 6 modes', 'considering 3 modes');
for j=1:3
j=j+3;
for k=1:3
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3))+phi(k,4)*phi(j,4)./(naturalfreq(4)^2-
w.^2+1i*2*xeta(4)*w*naturalfreq(4))+phi(k,5)*phi(j,5)./(naturalfreq(5)^2-
w.^2+1i*2*xeta(5)*w*naturalfreq(5))+phi(k,6)*phi(j,6)./(naturalfreq(6)^2-
w.^2+1i*2*xeta(6)*w*naturalfreq(6));
figure(3)
subplot(3,3,j-3+3*(k-1));
plot(f,20*log10(abs(alpha)));
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('freq Hz');
ylabel('\alpha dB');
hold on;
end
end
for j=1:3
j=j+3;
for k=1:3
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3));
figure(3)
subplot(3,3,j-3+3*(k-1));
plot(f,20*log10(abs(alpha)),'r');
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('freq Hz');
ylabel('\alpha dB');
hold on;
end
end
legend('considering 6 modes','considering 3 modes');
for j=1:3
for k=1:3
k=k+3;
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3))+phi(k,4)*phi(j,4)./(naturalfreq(4)^2-
w.^2+1i*2*xeta(4)*w*naturalfreq(4))+phi(k,5)*phi(j,5)./(naturalfreq(5)^2-
w.^2+1i*2*xeta(5)*w*naturalfreq(5))+phi(k,6)*phi(j,6)./(naturalfreq(6)^2-
w.^2+1i*2*xeta(6)*w*naturalfreq(6));
figure(4)
subplot(3,3,j+3*(k-3-1));
plot(f,20*log10(abs(alpha)));
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('freq Hz');
ylabel('\alpha dB');
hold on;
end
end
for j=1:3
for k=1:3
k=k+3;
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3));
figure(4)
subplot(3,3,j+3*(k-3-1));
plot(f,20*log10(abs(alpha)),'r');
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'fontsize',10);
xlabel('freq Hz');
ylabel('\alpha dB');
hold on;
end
end
legend('considering 6 modes', 'considering 3 modes');
for j=1:3
j=j+3;
for k=1:3
k=k+3;
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3))+phi(k,4)*phi(j,4)./(naturalfreq(4)^2-
w.^2+1i*2*xeta(4)*w*naturalfreq(4))+phi(k,5)*phi(j,5)./(naturalfreq(5)^2-
w.^2+1i*2*xeta(5)*w*naturalfreq(5))+phi(k,6)*phi(j,6)./(naturalfreq(6)^2-
w.^2+1i*2*xeta(6)*w*naturalfreq(6));
figure(5)
subplot(3,3,j-3+3*(k-3-1));
plot(f,20*log10(abs(alpha)));
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('freq Hz');
ylabel('\alpha dB');
hold on;
end
end
for j=1:3
j=j+3;
for k=1:3
k=k+3;
syms w;
w=0:0.01:70;
f=w/(2*pi);
alpha=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3));
figure(5)
subplot(3,3,j-3+3*(k-3-1));
plot(f,20*log10(abs(alpha)),'r');
heading=['\alpha','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('freq Hz');
ylabel('\alpha dB');
hold on;
end
end
legend('considering 6 modes', 'considering 3 moddes');
for j=1:6
k=1;
syms w;
w=0:0.01:70;
f=w/(2*pi);
response=phi(k,1)*phi(j,1)./(naturalfreq(1)^2-
w.^2+1i*2*xeta(1)*w*naturalfreq(1))+phi(k,2)*phi(j,2)./(naturalfreq(2)^2-
w.^2+1i*2*xeta(2)*w*naturalfreq(2))+phi(k,3)*phi(j,3)./(naturalfreq(3)^2-
w.^2+1i*2*xeta(3)*w*naturalfreq(3))+phi(k,4)*phi(j,4)./(naturalfreq(4)^2-
w.^2+1i*2*xeta(4)*w*naturalfreq(4))+phi(k,5)*phi(j,5)./(naturalfreq(5)^2-
w.^2+1i*2*xeta(5)*w*naturalfreq(5))+phi(k,6)*phi(j,6)./(naturalfreq(6)^2-
w.^2+1i*2*xeta(6)*w*naturalfreq(6));
figure(6)
subplot(1,6,j);
plot(f,abs(response)*10e3);
heading=['response','(',num2str(j),',',num2str(k),') vs freq'];
title(heading,'FontSize',14,'FontWeight','bold','FontName','Arial');
xlabel('freq Hz');
ylabel('response mm');
hold on;
end

RESULTS

NATURAL FREQUENCIES

1. 10.7332 Hz
2. 21.6154 Hz
3. 33.7464 Hz
4. 43.1530 Hz
5. 50.9517 Hz
6. 62.0344 Hz

DAMPED NATURAL FREQUENCIES

1. 10.7327 Hz
2. 21.6130 Hz
3. 33.7396 Hz
4. 43.1396 Hz
5. 50.9288 Hz
6. 61.9964 Hz

You might also like