不平衡磁拉力表达式

% 不平衡磁拉力表达式 Fx  Fy  %
  Rr = *;         %转子半径 m
  Lr = *;        %转子长 m

  kj = *;              %气隙基波磁动势系数 无单位
  Ij = *;            %励磁电流

  mu_0 = *;     %空气磁导系数 H/m  

  delta_0 = *;   %均匀气隙大小 m

  Fj = kj^2 * Ij^2;              %励磁电流的基波磁动势幅值

 

u（1）=x1;   (x1,y1为转子外圆几何中心坐标)

u(3)=y1; 

% 不平衡磁拉力表达式 Fx  Fy  %
  Fj = kj * Ij;              %励磁电流的基波磁动势幅值
  e = sqrt(u(1).^2 + u(3).^2);                           %转子轴心轨迹
  F_coeffi = R * L *pi * mu_0 *Fj^2 /delta_0.^2;
  Fx_ump = F_coeffi * (1/2 * e + 5/8 * e.^3) * u(1)./e;
  Fy_ump = F_coeffi * (1/2 * e + 5/8 * e.^3) * u(3)./e;       %不平衡磁拉力