Scribd
Upload
18 views1 page

Coding

This document contains equations and code to model the transmission of force through a coupled mechanical system with a foundation, screen, and driving frequency. It calculates and plots the force transmitted, foundation displacement, and screen displacement as functions of the driving frequency for both decoupled and fully coupled models. The key outputs are plots showing how these values change with frequency, with indications of any peaks in the displacement values.

Uploaded by

Syafiq Shebli
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
18 views1 page

Coding

This document contains equations and code to model the transmission of force through a coupled mechanical system with a foundation, screen, and driving frequency. It calculates and plots the force transmitted, foundation displacement, and screen displacement as functions of the driving frequency for both decoupled and fully coupled models. The key outputs are plots showing how these values change with frequency, with indications of any peaks in the displacement values.

Uploaded by

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

%%Transmission of Force

k_f=299680000;k_eq=579620;x_f=0.000009815;
m_f=66735.285;m=3500;f_0=222151;w_n=12.8688;w=75.3982;f_t=6664.53;
w=[0:150];
f_tD=abs((w_n.^2)./(w.^2-w_n.^2)).*f_0;
TR1=k_eq.*k_f-k_eq.*m_f.*w.^2;
TR2=((k_eq-m.*w.^2).*(k_eq+k_f-m_f.*w.^2))-k_eq.^2;
TR=(TR1./TR2);
f_tC=abs(TR.*f_0);
figure (1)
[Peak, PeakIdx] = findpeaks(f_tC);
plot(w,f_tD,'b--');hold on;
plot(w,f_tC,'r');hold on;
text(w(PeakIdx), Peak, sprintf('Peak=%6.3f', Peak));
ylabel('Force Tansmitted,F_T(m)');
xlabel('Frequency,(rad/s)');
title('Force Transmitted Vs Driving Frequency');
legend('De-Coupled','Fully Coupled');

%%Displacement of Foundation
XD=abs(f_t./(k_f-m_f.*w.^2));
TR1=k_eq;
TR2=((k_eq-m.*w.^2).*(k_eq+k_f-m_f.*w.^2))-k_eq.^2;
figure(2);
XC=abs((TR1./TR2).*f_0);
plot(w,XD,'b--');hold on;
plot(w,XC,'r');hold on;
ylabel('Foundation Displacement, X_F (m)');
xlabel('Frequency,(rad/s)');
title('Foundation Displacement Vs Driving Frequency');
legend('De-coupled','Fully Coupled');

%%Displacement of Screen
XD=abs(f_0./(k_eq-m.*w.^2));
TR1=(k_eq+k_f)-(k_eq.*m_f.*w.^2);
TR2=((k_eq-m.*w.^2).*(k_eq+k_f-m_f.*w.^2))-k_eq.^2;
figure (3);
XC=abs((TR1./TR2).*f_0);
[Peak, PeakIdx]=findpeaks(XD);
plot(w,XD,'r--');hold on;
plot(w,XC,'b');hold on;
text(w(PeakIdx), Peak, sprintf('Peak=%6.3f',Peak));
ylabel('Screen Displacement, X(m)');
xlabel('Frequency, (rad/s)');
title('Screen Displacement Vs Driving Frequency');
legend('De-Coupled','Fully Coupled');

You might also like