Newmark-Beta算例.doc_第1页
Newmark-Beta算例.doc_第2页
Newmark-Beta算例.doc_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

某四层框架结构在顶部受一个简谐荷载的作用，力的作用时间t1=5.0s，F0=100N。试计算响应的时间为20s（计算分2000步完成）。阻尼矩阵由Rayleigh阻尼构造。结构基本计算简图计算Newmark-beta法的matlab程序m=1,2,3,4;m=diag(m);k=800-80000;-8002400-16000;0-16004800-3200;00-32008000;c=0.05*m+0.02*k;f0=100;t1=5;nt=2000;dt=0.01;gama=0.5;beta=0.25;a0=1/ beta/dt/dt;a1=gama/ beta /dt;a2=1/beta/dt;a3=1/2/ beta -1;a4=gama/beta-1;a5=dt/2*(gama/ beta-2);a6=dt*(1-gama);a7=dt*gama;d=zeros(4,nt+1);v=zeros(4,nt+1);a=zeros(4,nt+1);tt=zeros(1,nt+1);fori=2:(nt+1)tt(1,1)=1*dt;tt(1,i)=i*dt;t=(i-1)*dt;if(tt1),f=f0*sin(4*pi*t/t1);0;0;0;else f=0;0;0;0;endke=k+a0*m+a1*c;fe=f+m*(a0*d(:,i-1)+a2*v(:,i-1)+a3*a(:,i-1)+c*(a1*d(:,i-1)+a4*v(:,i-1)+a5*a(:,i-1);d(:,i)=inv(ke)*fe;a(:,i)=a0*(d(:,i)-d(:,i-1)-a2*v(:,i-1)-a3*a(:,i-1);v(:,i)=v(:,i-1)+a6*a(:,i-1)+a7*a(:,i);endfigure; plot(tt(1,:),d(1,:)ylabel(位移/m,fontsize,20); grid on;xlabel(时间/s,f
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论