哪位做过王佳斌的无槽直线电机推力仿真？？？
附上论文，还有自己写的径向充磁永磁体内置时的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'); 论文忘记发了 题目是：A General Framework for the Analysis and Design of Tubular Linear Permanent Magnet Machines 比较经典的一篇论文 内置式的没有考虑硅钢片磁化曲线？ 论文里应该有考虑这个，我就想按照论文里的公式复现一遍，发现现在高词谐波的幅值越来越大 ，正常高次谐波应该幅值很小才对，看了很多天没有找到原因，拜托大佬帮我看看再。
可有偿帮忙！！
页:
[1]