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

heyong2002

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

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

FtpAdmin

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

yanyongju

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

21172485

这个是我以前写的，没多少注释，不过不难，希望可以帮你
% newmark方法的程序实现
clc
clear
m=[6/12,1/12;1/12,2/12];
k=[2*3,-1*2;-1*2,1*2];
f1=[0;1];
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=[0;0];
d1=[0;0];

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

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

njhebin

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

iewoug

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

欧阳中华

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

iewoug

原帖由 欧阳中华 于 2006-12-4 21:54 发表
.
    M、K和C时变时，威尔逊cita法或纽马克法仍然适用的根据那里有？

参加克拉夫－结构动力学 非线性结构反应分析部分，逐步积分法本就是按照K,C随时变来给出的。
平时我们所用大多为非时变的特殊简单情况