heyong2002 发表于 2006-11-21 20:52

求助多自由度威尔逊cita法或纽马克法的程序

质量和阻尼及刚度矩阵都随时间时变,求助多自由度威尔逊cita法或纽马克法的程序

哪位有多自由度威尔逊cita法或纽马克法的程序,matlab或c 的都行
本人需求一4自由度动力学方程响应,质量和阻尼及刚度矩阵都随时间时变,需求出最后的数值解.希望高手给指点下编程方法和技巧

FtpAdmin 发表于 2006-11-23 08:53

搜索论坛，论坛由不少这两种算法的程序

yanyongju 发表于 2006-11-23 17:13

www.ansys.edu8.com上有两个这样的程序,可以去看看

21172485 发表于 2006-11-24 08:19

这个是我以前写的，没多少注释，不过不难，希望可以帮你
% newmark方法的程序实现
clc
clear
m=;
k=;
f1=;
d2=inv(m)*f1;
b=0.25;
r=0.5;
dt=0.12;
a0=1/(b*(dt)^2);
a1=r/(b*dt);
a2=1/(b*dt);
a3=1/(2*b)-1;
a4=r/b-1;
a5=0.5*dt*((r/b)-2);
a6=dt*(1-r);
a7=r*dt;
pk=k+a0*m;
pk=inv(pk);
d0=;
d1=;

for i=1:10
t=0.12*i;
if t>0.5
f1(2)=0;
end
ff=0.0;
d00=0.0;
ff=ff+m*(a0*d0+a2*d1+a3*d2);
pf1=f1+ff;
d00=d00+pk*pf1;
d22=a0*(d00-d0)-a2*d1-a3*d2;
d11=d1+a6*d2+a7*d22;
d0=d00;
d1=d11;
d2=d22;
a(i)=d00(1);
b(i)=d00(2);
end
a
b

heyong2002 发表于 2006-11-24 14:39

我的M K C都与时间有关是时变的,不过先谢谢两位的帮助:handshake 了.

njhebin 发表于 2006-11-24 16:36

这两种方法对时变与时不变没有什么差异，用起来就让M K C随着时间步长的迭代改变就可以了！

iewoug 发表于 2006-12-4 19:56

离散开的话，M K C设置为多时间维的三维形式就可以了

欧阳中华 发表于 2006-12-4 21:54

.
    M、K和C时变时，威尔逊cita法或纽马克法仍然适用的根据那里有？

iewoug 发表于 2006-12-5 13:15

原帖由 欧阳中华 于 2006-12-4 21:54 发表
.
    M、K和C时变时，威尔逊cita法或纽马克法仍然适用的根据那里有？
参加克拉夫－结构动力学 非线性结构反应分析部分，逐步积分法本就是按照K,C随时变来给出的。
平时我们所用大多为非时变的特殊简单情况
页: [1]
查看完整版本: 求助多自由度威尔逊cita法或纽马克法的程序