位移图、放大系数以及自由振动、阻尼振动和强制阻尼振动的传递率（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客     

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

SDoF 系统的振动分析，位移图、放大系数以及自由振动、阻尼振动和强制阻尼振动的传递率。 

📚2 运行结果

部分代码：

%Let frequency of forcing element be w'
%let the ratio of forcing frequency to natural frequency be W, i.e. W = (w'/wn)

W = linspace(0,3);

zeta_1 = 0;
tr_1 = sqrt(1.+(W.*(2*zeta_1)).^2) ./ (sqrt(((1-(W.^2)).^2) + ((W.*(2*zeta_1)).^2)));
plot(W,tr_1)
hold on

zeta_2 = 0.2;
tr_2 = sqrt(1.+(W.*(2*zeta_2)).^2) ./ (sqrt(((1-(W.^2)).^2) + ((W.*(2*zeta_2)).^2)));
plot(W,tr_2)

zeta_3 = 0.5;
tr_3 = sqrt(1.+(W.*(2*zeta_3)).^2) ./ (sqrt(((1-(W.^2)).^2) + ((W.*(2*zeta_3)).^2)));
plot(W,tr_3)

zeta_4 = 0.6;
tr_4 = sqrt(1.+(W.*(2*zeta_4)).^2) ./ (sqrt(((1-(W.^2)).^2) + ((W.*(2*zeta_4)).^2)));
plot(W,tr_4)

xticks([0 0.5 1 sqrt(2) 2 2.5 3])
xticklabels({'0','0.5','1','√2','2','2.5', '3'})
plot(sqrt(2),0,'rx')
hold off
legend(["ζ=0","ζ=0.2","ζ=0.5","ζ=0.6"])
grid
xline(sqrt(2),"--")
xlim([0 3])
ylim([0 4])
xlabel('w/wn')
ylabel('Transmissibility')

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]R. Chittin Varunav (2023). Vibration analysis of SDoF systems

🌈4 Matlab代码实现