计算磁悬浮转子的固有频率问题

lnly

不知道大家有没有搞磁悬浮转子，我想问一下是不是在转子静止不工作的时候计算的频率才叫固有频率，而当转子工作的时候考虑外力的影响计算的就不是固有频率？
我是研究磁悬浮轴承－转子系统耦合，我想考虑磁悬浮轴承之间的耦合对转子固有频率的影响，但耦合力是在系统工作的时候由于外力的作用导致转子偏心才产生的，所以我不太清楚是不是这时计算的频率就不是固有频率了？
要是有谁是研究磁悬浮轴承－转子系统的希望大家指教指教啊！谢谢了！

咕噜噜

系统的固有频率和外界力无关，是线性系统的本征特征，外力的作用不会对系统固有频率产生影响，除非外力改变了系统的刚度或刚度矩阵
你的磁悬浮轴承－转子应该是一个系统，而且是不接触电磁作用力，在新的系统中你应该考虑整个系统可能存在哪几个频率，而不单单是转子。
不知道你的具体机构，所以无法肯定您的系统是几个自由度，要分清楚多自由度系统的频率和单自由度系统频率

lnly

我研究的整个系统就是一个磁悬浮轴承－转子系统，转子是靠两个径向磁力轴承悬浮起来的，还有一个轴向轴承控制转子的轴向位移．转子与轴承之间没有接触．系统是5自由度的，分别是x,y,z及绕x,y方向的转角.由于不考虑径向和轴向的耦合，可以把径向和轴向分开计算，所以可按４自由度计算．加耦合电磁力后，系统的刚度矩阵改变了．那就说耦合的电磁力影响转子的固有频率了，是吗？
＂你的磁悬浮轴承－转子应该是一个系统，而且是不接触电磁作用力，在新的系统中你应该考虑整个系统可能存在哪几个频率，而不单单是转子。＂这句话是什么意思啊，不就是转子在旋转吗，那不就是考虑转子的频率，还什么有频率啊？没明白？能说的再明白点吗？嘿嘿，我比较苯！什么都不懂！
下面是一个计算四自由度转子固有频率的程序，但总觉得好象不对却不知道错在哪里，能麻烦你帮忙看一下吗？谢谢了！
clc
clear
%syms w m l a b j1 j2 k1 k2 k3 k4 k5 k6
w=36000;m=4.94;l=1.981;a=0.98;b=1.1;j1=0.0066;j2=0.0133;k1=3.15;k2=0.018;k3=243.24;k4=204;k5=1.22e+6;k6=1.02e+6;

M=[m*b/l,m*a/l,0,0;0,0,m*b/l,m*a/l;-j2/l,j2/l,0,0;0,0,-j2/l,j2/l];
C=[k3*k2,k4*k2,0,0;0,0,k3*k2,k4*k2;-a*k3*k2,b*k4*k2,w*j1/l,-w*j1/l;-w*j1/l,w*j1/l,-a*k3*k2,b*k4*k2];
K=[(k5+k3*k1),(k6+k4*k1),0,0;0,0,(k5+k3*k1),(k6+k4*k1);-a*(k3*k1+k5),b*(k6+k4*k1),0,0;0,0,-a*(k3*k1+k5),b*(k6+k4*k1)];
I=eye(4);
mm=zeros(4);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%转换为标准特征值问题
G=[mm,M;-inv(M)*K,-inv(M)*C];
cri=eig(G);

pp=imag((cri));
wn=(abs(imag(cri)));
firstwn=wn(６)     %%%第一阶固有频率
kkk=min(abs((real(cri))))
secondwn=wn(３)   %%%第二阶固有频率
kkk=min(abs((real(cri))));
thirdwn=wn(１)    %%%第三阶固有频率通过程序算出的特征值，有正实数也有负实数，还有共轭复数，固有频率到底应取哪个值？