声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

123
返回列表 发新帖
楼主: yangxiaokang

[稳定性与分岔] 请教各位大侠转子-滑动轴承的碰摩问题

[复制链接]
发表于 2012-4-19 16:06 | 显示全部楼层
油膜力的程序是不对的,总有虚数部分泄露、
请参考
function oilforce=oilx(x,y,dx,dy,wi)
%R轴承半径,L轴承长度,wi激励频率,miu油膜黏度 dert为油膜平均厚度
R=0.025; L=0.012; miu=0.018;  dfai=1; dert=0.00011;
e=sqrt(x*x+y*y);
%delta为Sommerfeld修正系数
delta=miu*wi*R*L*(R/dert)*(R/dert)*(L/2.0/R)*(L/2.0/R);
%
ppp1=(y+2.0*dx)/(x-2.0*dy);
sign1=sign(ppp1);      
ppp2=y+2.0*dx;
sign2=sign(ppp2);
alpha=atan(ppp1)-pi/2.0*(sign1+sign2);
alphaa=atan((y*cos(alpha)-x*sin(alpha))/sqrt(abs(1.0-abs(x*x)-abs(y*y))));
fg=2.0*(pi/2.0+alphaa)/sqrt(abs(1.0-abs(x*x)-abs(y*y)));
fv=(2.0+(y*cos(alpha)-x*sin(alpha))*fg)/(1.0-abs(x*x)-abs(y*y));
fs=(x*cos(alpha)+y*sin(alpha))/(1.0-abs((x*cos(alpha)+y*sin(alpha))*(x*cos(alpha)+y*sin(alpha))));
f1=sqrt(abs(abs((x-2.0*dy)*(x-2.0*dy))+abs((y+2.0*dx)*(y+2.0*dx))))/(1.0-abs(x*x)-abs(y*y));
fx=-1.0*f1*(3.0*x*fv-sin(alpha)*fg-2.0*cos(alpha)*fs);
fy=-1.0*f1*(3.0*y*fv+cos(alpha)*fg-2.0*sin(alpha)*fs);
回复 支持 反对
分享到:

使用道具 举报

发表于 2012-8-22 18:44 | 显示全部楼层
得出的结果与原文的差距有点离谱,是不是应该换种积分方法?
发表于 2012-9-5 09:07 | 显示全部楼层
我也正做转子动力学呢   才列出来动力学方程 还没无量钢化呢啊
发表于 2012-9-11 15:48 | 显示全部楼层
发表于 2012-9-12 10:37 | 显示全部楼层
学习了,最近也在研究这方面
发表于 2012-9-17 14:30 | 显示全部楼层
这类转子油膜碰摩的文献里,大多数根据别人论文里的参数却做不出同样的结果来。
这种现象我遇见的太多了。真不知道怎么去解释好。同样的参数,按照别人同样的处理计算方法,同样的取点,就是做不出他们的那种结果来。可能是我自己有些问题,但有一些人的论文里就是有很明显错误的。
愿与广大同学共勉!

点评

应该把这种情况公布出来,让大家共同验证验证  发表于 2014-3-17 14:50
发表于 2012-9-17 14:46 | 显示全部楼层
可能是我自己有些问题,但有一些人的论文里就是有很明显错误的。
发表于 2013-6-26 16:54 | 显示全部楼层

21楼中定义的global w f   f=400:5:600和前面函数中f=0.1会不会有冲突?
发表于 2014-3-16 20:24 | 显示全部楼层
另外不知道你的分岔图怎么样了?现在对了吗?
“我自己认为是先做参数随频率的分岔图,然后在分岔出现的地方,取不同频率,再做轨迹图和映射图分析,不知道对否”
这个想法是对的,但是有时候不好选择参数的范围的时候,可以先尝试着计算一些点的轨迹图和映射图
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-23 20:11 , Processed in 0.089823 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表