转子碰摩分岔，附程序，请教各位指错

xiaoshihanlan · 发表于 2012-4-23 20:41

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

请教各位前辈老师，刚接触分岔这块，做一个Jeffcott转子系统碰摩的分岔图，跟着书上方程编写的matlab程序，就是出不来正确的图，请教各位帮忙看看是哪里错了？
方程：function dx=pengmo_1(t,x,flag,gamma)
beta=3.0;
alfa=0.5;
epsilon=0.12;
f=0.12;
deta=0.16;
u=0.136;
w0=50*pi;
g=9.8;
w=w0*gamma;
dx=[x(2);
-2*epsilon*x(2)-x(1)-beta*(1-1/sqrt(x(1)^2+x(3)^2))*(x(1)-f*x(3))-alfa*(x(1)^2+x(3)^2)*x(1)+u*(gamma)^2*cos(w0*gamma*t);
x(4);
-2*epsilon*x(4)-x(3)-beta*(1-1/sqrt(x(1)^2+x(3)^2))*(f*x(1)+x(3))-alfa*(x(1)^2+x(3)^2)*x(3)+u*(gamma)^2*sin(w0*gamma*t)-g/((w^2)*deta)];
主程序：w0=50*pi;
gamma=1.10:0.01:1.90;
x0=[0.1;0;0.1;0];
for h=1:length(gamma)
disp(gamma(h));
T=2*pi/(w0*gamma(h));
[t x]=ode45('pengmo_1',[0:T/100:200*T],x0,[],gamma(h));
plot(gamma(h),x(19500:200:end,1),'k.');hold on
end

xiaoshihanlan · 发表于 2012-4-23 20:49

hsfy919 · 发表于 2012-4-28 11:43

回复 2 # xiaoshihanlan 的帖子

1、既然是碰磨问题，那么判断是否接触的程序呢？
2、主程序中应该是用闪频法画的庞加莱截面，那么周期取点也不对。
3、你是在那本书上看到的？

xiaoshihanlan · 发表于 2012-4-29 09:51

回复 3 # hsfy919 的帖子

书上只给了微分方程，其它是自己写的，还要先判断是否发生碰摩？不懂。

伤痕累累 · 发表于 2012-4-29 10:31

回复 4 # xiaoshihanlan 的帖子

你这是哪本书上的？碰摩的程序根本就不对

xiaoshihanlan · 发表于 2012-4-29 11:59

回复 5 # 伤痕累累的帖子

在一本转子故障实验分析书上看到的，刚开始接触就跟着人家的方程做了,毕设设计做转子碰摩，现在就是程序搞不懂啊，请前辈帮帮忙

伤痕累累 · 发表于 2012-4-29 18:17

回复 7 # xiaoshihanlan 的帖子

你是本科是吧？

xiaoshihanlan · 发表于 2012-4-29 20:42

回复 8 # 伤痕累累的帖子

是的，在做毕业设计

账号		自动登录	找回密码
密码			我要加入

[稳定性与分岔] 转子碰摩分岔，附程序，请教各位指错

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。