原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
11/8/1512:45PMD:\study\TONGJI\StructuralEngineer...\interpolationmethod.m1of1
%%InterpolationMethod
%%Givendata
m=100;%%mass
k=100;%%stiffness
dr=0.02;%%dampingratio
v(1)=0;%%intialvelocityattimet=0
x(1)=0;%%initialdiscementattimet=0
t(1)=0;
dt=0.02;%%(timeintervalofElcentrodata)
p=-ElCentro*m;%%Inputearthquakeforce
wn=sqrt(k/m);
wd=wn*sqrt(1-dr^2);
dc=2*dr*sqrt(k*m);%%dampingcoefficient
acc(1)=(p(1)-dc*v(1)-k*x(1))/m;
%%Coefficientsofrecurrenceformula
A=exp(-dr*wn*dt)*(dr/sqrt(1-dr^2)*sin(wd*dt)+cos(wd*dt));
B=exp(-dr*wn*dt)*(1/wd*sin(wd*dt));
C=1/k*(2*dr/wn/dt+exp(-dr*wn*dt)*(((1-2*dr^2)/wd/dt-dr/sqrt(1-dr^2))*sin(wd*dt)-(1
+2*dr/wn/dt)*cos(wd*dt)));
D=1/k*(1-2*dr/wn/dt+exp(-dr*wn*dt)*((2*dr^2-1)/dt*sin(wd*dt)+2*dr/wn/dt*cos(wd*dt)));
a=-exp(dr*wn*dt)*(wn/sqrt(1-dr^2)*sin(wd*dt));
b=exp(-dr*wn*dt)*(cos(wd*dt)-dr/sqrt(1-dr^2)*sin(wd*dt));
c=1/k*(-1/dt+exp(-dr*wn*dt)*((wn/sqrt(1-dr^2)+dr/dt/sqrt(1-dr^2))*sin(wd*dt)+1/dt*cos
(wd*dt)));
d=1/k/dt*(1-exp(-dr*wn*dt)*(dr/sqrt(1-dr^2)*sin(wd*dt)+cos(wd*dt)));
%%Recurrenceformula
fori=1:1:15728
x(i+1)=A*x(i)+B*v(i)+C*p(i)+D*p(i+1);
v(i+1)=a*x(i)+b*v(i)+c*p(i)+d*p(i+1);
acc(i+1)=(p(i+1)-dc*v(i+1)-k*x(i+1))/m;
t(i+1)=t(i)+.02;
end
%%Plotofresponse
plot(t,acc,b-);
xlabel(Time(s));
ylabel(acceleration(m/s^2));
title(Responsespectruma-t)
gridon
11/8/1512:45PMD:\study\TONGJI\StructuralEngineer...\插值法.m1of1
%%插值法
%%给定数据m=100；%%质量k=100;%%刚度dr=0.02;%%阻
尼比v(1)=0;%%时刻t=0x(1)=0的初始速度；%%时刻t=0t(1)=0时
的初始位移；dt=0.02;%%（Elcentrodata的时间间隔）p‑
ElCentro*m;%%输入力wnsqrt(k/m);wdwn*sqrt(1‑dr
^2);
文档评论（0）