威尔逊theta法matlab程序及算例说明.docxVIP

下载本文档

24
0
约1.78千字
约 3页
2019-10-15 发布于上海
举报
版权申诉

威尔逊theta法matlab程序及算例说明.docx

1、本文档共3页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

单自由度： function [u,u1,u2]=wiltheta(m,k,c,d,e,f,p,theta,t) n=length(p); u=zeros(n,1);u1=zeros(n,1);u2=zeros(n,1);uu=zeros(n,1);pp=zeros(n,1); u(1)=d;u1(1)=e;u2(1)=f; l=6/(theta*t)^2;q=3/(theta*t);r=6/(theta*t);s=theta*t/2; kk=k+l*m+q*c; for i=1:n-1 pp(i)=p(i)+theta*(p(i+1)-p(i))+(l*u(i)+r*u1(i)+2*u2(i))*m+(q*u(i)+2*u1(i)+s*u2(i))*c; uu(i)=pp(i)/kk; u2(i+1)=u2(i)+l/theta*(uu(i)-u(i))-l*t*u1(i)-3/theta*u2(i); u1(i+1)=u1(i)+t/2*(u2(i+1)+u2(i)); u(i+1)=u(i)+t*u1(i)+t^2/6*(u2(i+1)+2*u2(i+1)); end end 实例：，初始时刻位移、速度、加速度均为零。以下为0.02-1时间段时程曲线。其中p=-1*[0.000508665 0.000906218 00000000.0040581 0000.0056473 0-0.000121305 -0-000000.001187971 00-0.000565749 -0-000-0-0000-0-0-0000.000469929 -0-02.34455E-05 0.0111213 0.000141692 -0.0116616 -0.0109857 00；多自由度 function [u u1 u2]=multi_wiltheta(M,K,C,d,e,f,P,theta,t) n=length(P(:,1));m=length(P(1,:)); u=zeros(n,m);u1=zeros(n,m);u2=zeros(n,m);uu=zeros(n,m);PP=zeros(n,m); u(:,1)=d;u1(:,1)=e;u2(:,1)=f; l=6/(theta*t)^2;q=3/(theta*t);r=6/(theta*t);s=theta*t/2; KK=K+l*M+q*C; for i=1:m-1 pp(:,i)=P(:,i)+theta*(P(:,i+1)-P(:,i))+M*(l*u(:,i)+r*u1(:,i)+2*u2(:,i))+C*(q*u(:,i)+2*u1(:,i)+s*u2(:,i)); uu(:,i)=KK\pp(:,i); u2(:,i+1)=u2(:,i)+l/theta*(uu(:,i)-u(:,i))-l*t*u1(:,i)-3/theta*u2(:,i); u1(:,i+1)=u1(:,i)+t/2*(u2(:,i+1)+u2(:,i)); u(:,i+1)=u(:,i)+t*u1(:,i)+t^2/6*(u2(:,i+1)+2*u2(:,i)); end end 实例：，假设初始时刻的位移、速度、加速度均为零。以下为0.1s-2s的时程曲线。