suifengsan 发表于 2019-5-31 16:24

做王家斌的无槽气隙磁密仿真,磁密仿出来很乱,帮忙看一下吧

附上论文,还有自己写的径向充磁永磁体内置时的matlab程序,仿真出来结果在前两次还正常,当谐波次数高了之后幅值很大,大致出来的磁密凌乱,同样的方法径向充磁永磁体外置的B就仿真出来了(这里只是求的气隙处磁密),很奇怪,帮忙看一下吧!!!

clc;clear all;
%永磁体内置式
num = 11; %11次谐波
g = 0.0025; %气隙宽度
taup = 0.06;%永磁体宽度
hm = 0.005; %永磁体径向厚度
taum = 0.044; %永磁体高度
hw = 0.0325; %绕组厚度
Rr = 0.078; %永磁体外径
                                    
miur = 1.05;%回复磁导率
Brem = 1.2;%剩磁         
Rm = 0.083; %永磁体内径                     
Rs = 0.118; %电枢内径
r = 0.08425;%气隙半径
alphap = taum/taup;
%利用卡特系数修正Rs
% taut=0.02;%齿距
% b0=0.0012;%槽宽
% g1=g+hm/miur;
% b=b0/(2*g1);
% v=(b0/(2*g1*tan(b))-0.5*log(1+b*2))*(4/pi);
% kc=taut/(taut-v*g1);%卡特系数
% ge=g+(kc-1)*g1;
% Rs=Rm+ge;

for n = 1:num

      
      Mn = (2*n-1)*pi/taup;
      Pn = (4/taup)*Brem*sin(2*n-1)*(pi/2)*alphap;
      a = 1e-3;
      fa = 0;
      fb = 0;
   for x = Mn*Rr:a:Mn*Rm
      f1 = besselk(1,x)/(besseli(1,x)*besselk(0,x)+besselk(1,x)*besseli(0,x));
      f2 = besseli(1,x)/(besseli(1,x)*besselk(0,x)+besselk(1,x)*besseli(0,x));
      fa = fa+f1;
      fb = fb+f2;
      end
      FAn = Pn/Mn*a*fa;
      FBn = Pn/Mn*a*fb;
      C1n = besseli(0,Mn*Rs);
      C2n=besselk(0,Mn*Rs);
      C3n = besseli(0,Mn*Rr);
      C4n=besselk(0,Mn*Rr);
      C5n = besseli(0,Mn*Rm);
      C6n=besselk(0,Mn*Rm);
      C7n = besseli(1,Mn*Rm);
      C8n=besselk(1,Mn*Rm);
      a1n = (FAn*C5n/C6n+FBn)/(miur*(C5n/C6n-C1n/C2n)-(C5n/C6n-C3n/C4n));
      a2n = (FAn*C7n/C8n-FBn)/((C7n/C8n+C1n/C2n)-(C7n/C8n+C3n/C4n));
      b1n = C1n/C2n*a1n;
      b2n = C3n/C4n*a2n;
      zn = 0;
      for z = 0:1e-3:2*taup
            zn = zn+1;
            B1z(n,zn) = (a1n*besseli(0,Mn*r)-b1n*besselk(0,Mn*r))*sin(Mn*z);
            B1r(n,zn) = -((a1n*besseli(1,Mn*r)+b1n*besselk(1,Mn*r))*cos(Mn*z));
      end
      B1zn=sum(B1z);B1rn=sum(B1r);
end


z = 0:1e-3:2*taup;
figure(1)
plot(z,B1zn,'g',z,B1rn,'b')
title('气隙处沿z方向与r方向磁密大小');
xlabel('Z');
ylabel('B');
legend('z','r');

suifengsan 发表于 2019-5-31 16:27

这篇论文的附件 其中有公式, 图片是仿真出来的永磁体内置时气隙磁密的结果

suifengsan 发表于 2019-12-18 14:33

已经解决!!!!!!!!!!!!!!!!!!!
页: [1]
查看完整版本: 做王家斌的无槽气隙磁密仿真,磁密仿出来很乱,帮忙看一下吧