原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
单自由度体系杜哈梅积分
function y=kst(t0,t1,t2,ts,m,b0,b1,w0,c)
t0=input(请输入起始时间:t0= );t1=input(请输入荷载消失时间:t1= );t2=input(请输入想要的时间:t2= ); ts=input(请输入时间步长:ts= );
m=input(请输入质量:m= ) ;b0=input(请输入荷载截距:b0= );b1=input(荷载消失时的荷载:b1= );k=input(请输入刚度:k= ) ;
c=input(请输入阻尼比:c= );
w0=sqrt(k/m);w1=w0*sqrt(1-c^2);
t=t0:ts:t2;
for i=1:(length(t))
x=linspace(t(1),t(length(t)))
p=interp1([t0 t1],[b0 b1],t);
p(find(isnan(p)==1)) = 0;
px=linspace(p(1),p(length(t)));
a=px.*exp(c*w0*x).*cos(w1*x);
A=trapz(x,a);
b=px.*exp(c*w0*x).*sin(w1*x);
B=trapz(x,b);
y=exp(-c*w0*t).*(A.*sin(w1*t)-B.*cos(w1*t))./(m*w1)
v=diff(y)
a0=diff(y,2)
end
ymax=max(y)
figure
plot(t,y);
此程序为复合梯形法计算冲击荷载作用下的杜哈梅积分。
以P(t)=-1250000*(t+0.08)的冲击荷载为例，质量：m=6.4；阻尼比c=0.05；刚度：k=34847.77 N/m.将参数输入程序得到以下结果：
冲击荷载作用下的体系的速度图
冲击荷载作用下的体系加速度图
文档评论（0）