houzhengqi 发表于 2006-5-8 15:43

Jeffcott转子模型--方程求解

本帖最后由 VibInfo 于 2016-4-22 15:23 编辑

  我编写了一个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 = ;

  Tinterv = 0.1;

  Tspan = 300;

   = ode45(@JeffcottFun,,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=;

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

houzhengqi 发表于 2006-5-8 15:45

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

yejet 发表于 2006-5-8 15:59

回复:(houzhengqi)Jeffcott转子模型--方程求解

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

houzhengqi 发表于 2006-5-8 17:01

回复:(houzhengqi)Jeffcott转子模型--方程求解

本帖最后由 VibInfo 于 2016-4-22 15:23 编辑

  谢谢 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 发表于 2006-5-8 17:13

回复:(houzhengqi)Jeffcott转子模型--方程求解

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

toes 发表于 2006-5-8 17:19

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

jgwang 发表于 2006-10-8 20:56

请教楼主:问题解决了吗?

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

jgwang 发表于 2006-12-3 10:50

还是没有弄明白,哪位高人指点一下?

1.MATLAB编程时,无量纲方程可以直接用来积分吗?
2.访真时发现,初值x0,积分步长,积分时间等都影响最终分叉结果,怎样取这些值呢?
谢谢诸位!!!
页: [1]
查看完整版本: Jeffcott转子模型--方程求解