newmark迭代的MATLAB程序实现

21172485

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

复制代码

bd040590

没说明 叫人怎么用

double29

这是什么啊？不知道哦！

lemon1983

建议搂主给大家补上说明

cao

这个程序其实较好理解，先搞懂newmark方法，一般的动力学有关的书上会有newmark方法的编程思路，结合着看一下，应该十分钟内您就会明白．
第一段是一些初值设置；第二段就是迭代算法！
批评指正，相互学习！

lemon1983

f1(2)=0;什么意思？

xinyuxf

楼主能不能说得明白一些啊？谢谢，想学习一下

zjlygwt

原帖由 lemon1983 于 2006-8-30 20:11 发表
f1(2)=0;什么意思？

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

baiyang

请问版主有没有newmark分析分叉混沌的程序，我们共同探讨一下，qq254785836

suffer

原帖由 baiyang 于 2006-11-12 17:40 发表
请问版主有没有newmark分析分叉混沌的程序，我们共同探讨一下，qq254785836

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

当然对于非线性问题，一般认为增量形式的newmark法效果会比较好

suffer

原帖由 zjlygwt 于 2006-9-11 13:56 发表

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

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

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

[ 本帖最后由 suffer 于 2006-11-13 10:18 编辑 ]

daihao213

最好补充一下注释

realyyy

见讨论：http://forum.vibunion.com/forum/ ... ighlight=%2Brealyyy

sally1982

请问有没有wilson法的程序，俺在毕业设计里面用到的，可是在f(t+θ*dt)里面不知道怎么叠带，还请帮忙

aishuren

d00那样求对吗？与书上讲的相差太大，得出的结果一定是发散的