Jeffcott转子模型--方程求解

houzhengqi

　　我编写了一个Jeffcott转子动力学信号求解程序,发现总是不对，现将程序代码上传，请教大虾帮帮忙。

　　%%%%%%%%%%%%%主程序

　　clc;

　　clear;

　　global theta betta E kesai f f0 gama w0 w

　　theta = 0.16; %转静间隙

　　betta = 3.00; %刚度比

　　E = 0.04; %偏心矩

　　kesai = 0.12; %阻尼比

　　f = 0.12; %摩擦系数

　　f0 = 25; %转子固有频率

　　gama = 1.86; %频率比

　　w0 = 2*pi*f0;

　　w = w0*gama;

　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 调用函数

　　x0 = [0.086,0.23,0.149,0.18,0];

　　Tinterv = 0.1;

　　Tspan = 300;

　　[t,x] = ode45(@JeffcottFun,[0:Tinterv:Tspan],x0);N1=length(t)-1;

　　N2=round(0.5*N1)+1;

　　tt=t(N2:N1);

　　tt=tt-t(N2);

　　xx=x(N2:N1,:);

　　N=length(tt);

　　%%%%%%%时域波形图

　　subplot(2,2,1)

　　plot(tt,xx(:,1))

　　%%%%%%%%%%%%%调用程序

　　function FUN=JeffcottFun(t,x);

　　global theta betta E kesai f f0 gama w0 w

　　G=9800/(theta*w0^2);

　　e=sqrt(x(1)^2+x(3)^2);

　　if e>1

　　kr=-betta*(1-1/e);

　　else

　　kr=0;

　　end

　　x(2)

　　x(4)

　　z1 = -2*kesai*x(2) -x(1) +kr*(x(1)-f*x(3)) +E*gama^2*cos(gama*x(5)) ;

　　z2 = -2*kesai*x(4) -x(3) +kr*(f*x(1)+x(3)) +E*gama^2*sin(gama*x(5)) -G ;

　　FUN=[x(2);z1;x(4);z2;1];

　　哪位XDJM看过了，给点意见啊，谢谢

houzhengqi

补充说明：使用Matlab编的程序。

yejet

这是一个碰摩转子模型吧，把你的无量纲化方程给一下

houzhengqi

　　谢谢 yejet

　　无量纲方程在调用程序function FUN=JeffcottFun(t,x)里面

　　其中：

　　x(1)，x(2)分别代表无量纲横向和纵向 位移

　　x(3)，x(4)分别代表无量纲横向和纵向 速度

　　z1 , z2 分别代表无量纲横向和纵向 加速度

　　z1 = -2*kesai*x(2) -x(1) +kr*(x(1)-f*x(3)) +E*gama^2*cos(gama*x(5)) ;

　　z2 = -2*kesai*x(4) -x(3) +kr*(f*x(1)+x(3)) +E*gama^2*sin(gama*x(5)) -G ;

　　x(5)代表无量纲时间 wt

　　w=gama*w1,其中w1为转子固有频率，gama为频率比

　　调用ode45解方程，

yejet

这个我知道，只是感觉JeffcottFun写的和无量纲化的方程好像不对应

toes

就是就是，你给的方程和你的程序对不上，位移和速度搞混了。

jgwang

小弟刚刚接触转子动力学问题，用MATLAB仿真的效果不尽人意，特求教此问！
－－如果哪位XDJM知道，Jeffcott转子动力学的MATLAB仿真也请回答为盼！！
谢谢！！！

jgwang

1.MATLAB编程时，无量纲方程可以直接用来积分吗？
2.访真时发现，初值x0，积分步长，积分时间等都影响最终分叉结果，怎样取这些值呢？
谢谢诸位！！！