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

suifengsan

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

1. 
   
2. clc;clear all;
   
3. %永磁体内置式
   
4. num = 11; %11次谐波
   
5. g = 0.0025; %气隙宽度
   
6. taup = 0.06;%永磁体宽度
   
7. hm = 0.005; %永磁体径向厚度
   
8. taum = 0.044; %永磁体高度
   
9. hw = 0.0325; %绕组厚度
   
10. Rr = 0.078; %永磁体外径
    
11.                                     
    
12. miur = 1.05;%回复磁导率
    
13. Brem = 1.2;%剩磁         
    
14. Rm = 0.083; %永磁体内径                     
    
15. Rs = 0.118; %电枢内径
    
16. r = 0.08425;%气隙半径
    
17. alphap = taum/taup;
    
18. %利用卡特系数修正Rs
    
19. % taut=0.02;%齿距
    
20. % b0=0.0012;%槽宽
    
21. % g1=g+hm/miur;
    
22. % b=b0/(2*g1);
    
23. % v=(b0/(2*g1*tan(b))-0.5*log(1+b*2))*(4/pi);
    
24. % kc=taut/(taut-v*g1);%卡特系数
    
25. % ge=g+(kc-1)*g1;
    
26. % Rs=Rm+ge;
    
27. 
    
28. for n = 1:num
    
29. 
    
30.         
    
31.         Mn = (2*n-1)*pi/taup;
    
32.         Pn = (4/taup)*Brem*sin(2*n-1)*(pi/2)*alphap;
    
33.         a = 1e-3;
    
34.         fa = 0;
    
35.         fb = 0;
    
36.      for x = Mn*Rr:a:Mn*Rm
    
37.         f1 = besselk(1,x)/(besseli(1,x)*besselk(0,x)+besselk(1,x)*besseli(0,x));
    
38.         f2 = besseli(1,x)/(besseli(1,x)*besselk(0,x)+besselk(1,x)*besseli(0,x));
    
39.         fa = fa+f1;
    
40.         fb = fb+f2;
    
41.       end
    
42.         FAn = Pn/Mn*a*fa;
    
43.         FBn = Pn/Mn*a*fb;
    
44.         C1n = besseli(0,Mn*Rs);
    
45.         C2n=besselk(0,Mn*Rs);
    
46.         C3n = besseli(0,Mn*Rr);
    
47.         C4n=besselk(0,Mn*Rr);
    
48.         C5n = besseli(0,Mn*Rm);
    
49.         C6n=besselk(0,Mn*Rm);
    
50.         C7n = besseli(1,Mn*Rm);
    
51.         C8n=besselk(1,Mn*Rm);
    
52.         a1n = (FAn*C5n/C6n+FBn)/(miur*(C5n/C6n-C1n/C2n)-(C5n/C6n-C3n/C4n));
    
53.         a2n = (FAn*C7n/C8n-FBn)/((C7n/C8n+C1n/C2n)-(C7n/C8n+C3n/C4n));
    
54.         b1n = C1n/C2n*a1n;
    
55.         b2n = C3n/C4n*a2n;
    
56.         zn = 0;
    
57.         for z = 0:1e-3:2*taup
    
58.             zn = zn+1;
    
59.             B1z(n,zn) = (a1n*besseli(0,Mn*r)-b1n*besselk(0,Mn*r))*sin(Mn*z);
    
60.             B1r(n,zn) = -((a1n*besseli(1,Mn*r)+b1n*besselk(1,Mn*r))*cos(Mn*z));
    
61.         end
    
62.         B1zn=sum(B1z);B1rn=sum(B1r);
    
63. end
    
64. 
    
65. 
    
66. z = 0:1e-3:2*taup;
    
67. figure(1)
    
68. plot(z,B1zn,'g',z,B1rn,'b')
    
69. title('气隙处沿z方向与r方向磁密大小');
    
70. xlabel('Z');
    
71. ylabel('B');
    
72. legend('z','r');

复制代码

suifengsan

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

suifengsan

已经解决！！！！！！！！！！！！！！！！！！！