有没有哪位高人做过Muszynska的密封力模型
本帖最后由 fanhuasijin 于 2011-5-5 09:20 编辑我现在做的是密封力作用下的转子非线性研究，转子是单圆盘Jeffcott转子模型，密封力采用的是
Muszynska模型。用matlab仿真结果总是不收敛，有没有人遇到过这种问题，最后怎么解决的。请大家多指教。下面就是该模型的M文件 function dxdt=Muszynska(t,x)
global w;
Ke=1.0e+008;De=500; %轴的刚度、转子外阻尼
u=2.0e-005; %流体粘度l=0.05; %密封长度 c=0.001; %密封间隙 R=0.5; %转子半径 tao1=0.49; b=0.15;m=3000; %转子质量 m0=-0.25;n0=0.066;v=50; % 轴向速度p=200; % 流体密度 z=0.1; %水力损失 n=2.5; r=0.2;%转子的相对偏心距 %未标注的为具体的密封参数
Ra=p*v*c/u; %轴向雷诺数
Rv=p*R*c*w/u; %周向雷诺数
ac=Ra/Rv;
lam=n0*Ra^(m0)*(1+(1/(2*ac))^2)^((1+m0)/2);
o=lam*l/c;
T=1/v;
B1=1/(1+(2*ac)^2);
ap=(1+z+2*o)*p*v^2/2; %密封前后压差
B=1+4*ac^2*B1*(1+m0);
E=(1+z)/(2*(1+z+2*o));
u0=2*o^2*E*(1-m0)/(1+z+2*o);
u1=2*o^2*(E/o+B*(1/6+E)/2)/(1+z+2*o);
u2=o*(1/6+E)/(1+z+2*o);
u3=pi*ap*R/lam;
mf=u2*u3*T^2;
M=m+mf;
K0=u0*u3;
D0=u1*u3*T;
tao=tao1*(1-r)^b;
K=K0*(1-r)^(-n);
D=D0*(1-r)^(-n);
K1=(Ke+K-tao^2*w^2*mf)/(M*w^2);
K2=D*tao/(M*w);
D1=(De+D)/(M*w);
D2=2*mf*tao/M;
G=9.8*m/(M*c*w^2);
dxdt=做过的高手进来帮帮忙吧，急呀
以上所有参数值均来自一篇已发表的博士论文
下面是该模型的无量纲模型，但是后面缺了转子质量偏心力，在M文件里有
页:
[1]