声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 15577|回复: 15

[共享资源] newmark迭代的MATLAB程序实现

[复制链接]
发表于 2006-4-23 21:54 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
  1. % newmark方法的程序实现
  2. clc
  3. clear
  4. m=[6/12,1/12;1/12,2/12];
  5. k=[2*3,-1*2;-1*2,1*2];
  6. f1=[0;1];
  7. d2=inv(m)*f1;
  8. b=0.25;
  9. r=0.5;
  10. dt=0.12;
  11. a0=1/(b*(dt)^2);
  12. a1=r/(b*dt);
  13. a2=1/(b*dt);
  14. a3=1/(2*b)-1;
  15. a4=r/b-1;
  16. a5=0.5*dt*((r/b)-2);
  17. a6=dt*(1-r);
  18. a7=r*dt;
  19. pk=k+a0*m;
  20. pk=inv(pk);
  21. d0=[0;0];
  22. d1=[0;0];

  23. for i=1:10
  24. t=0.12*i;
  25. if t>0.5
  26. f1(2)=0;
  27. end
  28. ff=0.0;
  29. d00=0.0;
  30. ff=ff+m*(a0*d0+a2*d1+a3*d2);
  31. pf1=f1+ff;
  32. d00=d00+pk*pf1;
  33. d22=a0*(d00-d0)-a2*d1-a3*d2;
  34. d11=d1+a6*d2+a7*d22;
  35. d0=d00;
  36. d1=d11;
  37. d2=d22;
  38. a(i)=d00(1);
  39. b(i)=d00(2);
  40. end
  41. a
  42. b
复制代码

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2006-6-14 10:46 | 显示全部楼层
没说明 叫人怎么用
发表于 2006-8-25 15:23 | 显示全部楼层
这是什么啊?不知道哦!
发表于 2006-8-28 21:31 | 显示全部楼层
建议搂主给大家补上说明
发表于 2006-8-29 09:23 | 显示全部楼层
这个程序其实较好理解,先搞懂newmark方法,一般的动力学有关的书上会有newmark方法的编程思路,结合着看一下,应该十分钟内您就会明白.
第一段是一些初值设置;第二段就是迭代算法!
批评指正,相互学习!
发表于 2006-8-30 20:11 | 显示全部楼层
f1(2)=0;什么意思?
发表于 2006-8-31 11:35 | 显示全部楼层
楼主能不能说得明白一些啊?谢谢,想学习一下
发表于 2006-9-11 13:56 | 显示全部楼层
原帖由 lemon1983 于 2006-8-30 20:11 发表
f1(2)=0;什么意思?

f是外菏载向量???
下面迭代里面的f怎么又对这个进行说明??
发表于 2006-11-12 17:40 | 显示全部楼层

太牛了

请问版主有没有newmark分析分叉混沌的程序,我们共同探讨一下,qq254785836
发表于 2006-11-13 10:15 | 显示全部楼层
原帖由 baiyang 于 2006-11-12 17:40 发表
请问版主有没有newmark分析分叉混沌的程序,我们共同探讨一下,qq254785836


newmark法适用于分析微分方程的直接积分法
无所谓分析什么问题,只要你能把问题抽象成微分方程的形式,基本可以用该算法进行分析

当然对于非线性问题,一般认为增量形式的newmark法效果会比较好
发表于 2006-11-13 10:16 | 显示全部楼层
原帖由 zjlygwt 于 2006-9-11 13:56 发表

f是外菏载向量???
下面迭代里面的f怎么又对这个进行说明??


f就是标准微分方程的右端项
后边是否要变动者完全取决于右端项的形式

比如在正弦激励下,显然他是和时间有关的

[ 本帖最后由 suffer 于 2006-11-13 10:18 编辑 ]
发表于 2006-11-27 19:02 | 显示全部楼层
最好补充一下注释
发表于 2006-12-7 10:20 | 显示全部楼层
发表于 2006-12-12 11:32 | 显示全部楼层
请问有没有wilson法的程序,俺在毕业设计里面用到的,可是在f(t+θ*dt)里面不知道怎么叠带,还请帮忙
发表于 2008-9-24 10:09 | 显示全部楼层

d00那样求对吗?与书上讲的相差太大,得出的结果一定是发散的
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-10 16:10 , Processed in 0.064588 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表