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看过了，给点意见啊，谢谢
补充说明：使用Matlab编的程序。
回复：（houzhengqi）Jeffcott转子模型--方程求解
这是一个碰摩转子模型吧，把你的无量纲化方程给一下回复：（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解方程，
回复：（houzhengqi）Jeffcott转子模型--方程求解
这个我知道，只是感觉JeffcottFun写的和无量纲化的方程好像不对应 就是就是，你给的方程和你的程序对不上，位移和速度搞混了。请教楼主：问题解决了吗？
小弟刚刚接触转子动力学问题，用MATLAB仿真的效果不尽人意，特求教此问！－－如果哪位XDJM知道，Jeffcott转子动力学的MATLAB仿真也请回答为盼！！
谢谢！！！
还是没有弄明白，哪位高人指点一下？
1.MATLAB编程时，无量纲方程可以直接用来积分吗？2.访真时发现，初值x0，积分步长，积分时间等都影响最终分叉结果，怎样取这些值呢？
谢谢诸位！！！
页:
[1]