% 1.無阻尼體系的簡諧荷載反應
clc;clear;
t=0:0.05:30;%輸入反應時長
m=input;%輸入體系質量，例：m=20；
k=input;%輸入體系剛度，例：k=40；
cwy=2;%輸入體系初位移；
csd=1;輸入體系初速度；
A=cwy;
wn=sqrt(k/m);
w=5;
beta=w/wn;%頻率比
P0=10;
B=csd/wn-P0/k*(beta/(1-beta*beta));
q1=A*cos(wn*t)+B*sin(wn*t)+P0/k*(1/((1-beta*beta)))*sin(wn*t);
subplot(1,2,1);
plot(t,q1,'r')
?xlabel('t')
?ylabel('位移')
?title(['無阻尼體系的簡諧荷載反應'])
?
?% 2.有阻尼體系的簡諧荷載反應
kexi=0.05;%阻尼比
wd=wn*sqrt(1-kexi*kexi);
ust=0.5;%等效靜位移
C=ust.*(1-beta*beta)/((1-beta*beta)^2+(2*kexi*beta)^2);
D=-ust.*(2*kexi*beta)/((1-beta*beta)^2+(2*kexi*beta)^2);
q2=exp(-kexi*wn*t).*(A*cos(wd*t)+B*sin(wd*t))+C*sin(w*t)+D*cos(w*t);
subplot(1,2,2);
plot(t,q2,'r')
?xlabel('t')
?ylabel('位移')
?title(['有阻尼體系的簡諧荷載反應'])

?

總結

以上是生活随笔為你收集整理的单自由度体系对简谐荷载的反应的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。