单自由度非线性体系线加速度法程序编写_第1页
单自由度非线性体系线加速度法程序编写_第2页
单自由度非线性体系线加速度法程序编写_第3页
单自由度非线性体系线加速度法程序编写_第4页
单自由度非线性体系线加速度法程序编写_第5页
全文预览已结束

下载本文档

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

文档简介

1、单自由度非线性体系线加速度法程序编写本程序主要应用线加速度法数值求解单自由度非线性体系的位移曲线。整个程序通过Matlab编写。1.线加速法原理1.1增量方程t时刻：f(t)+f(t)+f(t)=P(t)IDst+AtI时刻：f(t+At)+f(t+At)+f(t+At)=P(t+At)IDsAf(t)+Af(t)+Af(t)=AP(t)其中：Af(t)=f(t+At)-f(t)=my(t+At)-my(t)=m-Ay(t)IIIAf(t)Qk(t)-Ay(t)AfD(t)-c(t)-Ay(t)AP(t)=P(t+At)-P(t)贝U可得增量方程mAy(t)+c(t)Ay(t)+k(t)Ay(

2、t)=AP(t)1.2线加速度法Ay(t)二y(t)At+罟(At)2+Ay(t)二y(Ay(t)二y(t)At+罟(At)2+y(4)(t)雄Ay(t)-*(t)-3y(t)Ay(t)=y(t)At+雄Ay(t)-*(t)-3y(t)Ay(t)=3Ay(t)-3y(t)-Ay(t)=3Ay(t)-3y(t)-Aty(t)At2Aty(t+At)=y(t)+Ay(t)y(t+At)=y(t)+Ay(t)y(t+At)=y(t)+Ay(t)2.程序编写流程图Ay=AP/TOC o 1-5 h z6At63其中AP(t)=AP(t)+m丁y(t)+3y(t)+c(t)3y(t)+片y(t)k(t)

3、=k(t)m+c(t)2(At)2Aty(t)=y(t)+Ay(t) HYPERLINK l bookmark4 o Current Document ii-1i-1y(t)=y(t)+Ay(t)ii-1i-1y(t)=-P(t)-f(t)-f(t)imiDisi3.具体应用算例P(t)zb3kNt(s)0+050.80.13.具体应用算例P(t)zb3kNt(s)0+050.80.1IW=15kN、y(t)例：求位移时程曲线，一恢复力时程曲线，p(t)最大位移，最大恢复力，(亦尹始时静止g卜7通过所编写的程序，得到以下结果：当At=0.1s时，ymax=0.1423,fsmax=3y-t曲线fs-t曲线当At=0.05s时，ymax=0.1563,fsmax=3当At=0.025s时，ymax=0.1583,fsmax=3从以上图像中可以看到，当t越小时，所得到的结构精度越高，但可以看到当

人人文库> 全部分类> 图纸下载 > 毕业设计

温馨提示

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

最新文档

评论

0/150

提交评论