转子-轴承油膜力程序问题，求指导
本帖最后由 伤痕累累 于 2012-4-19 17:23 编辑为什么时域图不是按照周期波动的，而是一直递减或者递增呢？
在文献中，作者给的油膜力公式就是这样的啊。按照公式来，是没有错的。倒搞不通为什么是这样的结果
大家有兴趣的一起来讨论下吧
function book_oilfilm_force
clear all
clc
y0=;
ff=2*pi;
fs=256*3;
tspan=;
options=odeset('rel',1e-6,'abs','1e-6');
=ode45(@oilfilm,tspan,y0,options);
plot(t,y(:,1));
title('水平方向时域图')
xlabel('时间/t')
ylabel('x')
figure
plot(y(:,1),y(:,3))
title('phase figure')
xlabel('x')
ylabel('dx')
function yy=oilfilm(t,y)
%油膜力表达式
%%%%%%%%%%%%%%%%
ppp1=(y(2)+2.0*y(3))/(y(1)-2.0*y(4));
sign1=sign(ppp1);
ppp2=y(2)+2.0*y(3);
sign2=sign(ppp2);
alpha=atan(ppp1)-pi/2.0*(sign1+sign2);
alphaa=atan((y(2)*cos(alpha)-y(1)*sin(alpha))/sqrt(abs(1.0-abs(y(1)*y(1))-abs(y(2)*y(2)))));
fg=2.0*(pi/2.0+alphaa)/sqrt(abs(1.0-abs(y(1)*y(1))-abs(y(2)*y(2))));
fv=(2.0+(y(2)*cos(alpha)-y(1)*sin(alpha))*fg)/(1.0-abs(y(1)*y(1))-abs(y(2)*y(2)));
fs=(y(1)*cos(alpha)+y(2)*sin(alpha))/(1.0-abs((y(1)*cos(alpha)+y(2)*sin(alpha))*(y(1)*cos(alpha)+y(2)*sin(alpha))));
f1=sqrt(abs(abs((y(1)-2.0*y(4))*(y(1)-2.0*y(4)))+abs((y(2)+2.0*y(3))*(y(2)+2.0*y(3)))))/(1.0-abs(y(1)*y(1))-abs(y(2)*y(2)));
fx=-1.0*f1*(3.0*y(1)*fv-sin(alpha)*fg-2.0*cos(alpha)*fs);
fy=-1.0*f1*(3.0*y(2)*fv+cos(alpha)*fg-2.0*sin(alpha)*fs);
%%%%%%%%%%%%%%%%%%%
%参数
segma=0.1;%文中的 σ
p=0.3;%小写的平，为 偏心率ρ
M=70.0;%文中的m
P=170/2.0;%转子重力，大写的P
Fx=fx/(segma*P);%无量纲化的油膜力
Fy=fy/(segma*P);
G=1/(segma*M); %无量纲化后的重力
yy=[y(3);
y(4);
p*sin(t)+Fx/M;
p*cos(t)-G+Fy/M;
];
本帖最后由 伤痕累累 于 2012-4-19 17:00 编辑
论文里的结果 本帖最后由 伤痕累累 于 2012-4-19 17:01 编辑
这个程序的结果 我想问下，哪篇论文呀 程序是你自己编写的吗？怎么没有注释呀 本帖最后由 伤痕累累 于 2012-4-19 17:16 编辑
回复 4 # paopaotaiqiang 的帖子
title: Chaotic motions of a rigid rotor in short journal bearings.
author :G.Adiletta,A.R.Guido.1996.Nonlinear Dynamics
如果下载不方便，我可以传递给你
回复 5 # paopaotaiqiang 的帖子
是我自己写的，好吧，我去注释一下 我想问下，你知道分叉图是怎么画吗？我写出来程序了，就是轴心轨迹的程序，但不知道分叉图要怎么做，关于轴承碰磨方面的 回复 8 # paopaotaiqiang 的帖子
分插图好做。 那要怎么做呢？给我说说吧，学习下 回复 10 # paopaotaiqiang 的帖子
这类程序在论坛里特别多。建议你到matlab版块去转转，会有收获的 这是我画的轴心轨迹 取了部分点做的图 回复 13 # paopaotaiqiang 的帖子
不收敛，而且近似圆。。。 你程序修改完啦，上传我看看哈，也许能帮你找到问题
页:
[1]
2