麦克斯韦张量法（2）

麦克斯韦张量法（2）

上一篇文章分享了麦克斯韦张量法的原理，今天给大家介绍如何基于matlab实现！本文以电机负载一个时刻下沿圆周方向进行分析。

将气隙磁密谐波的径向和切向磁密的FFT后的幅值和角度导出后，接着导入matlab中，保存好变量工作区后，进行代码编写调试。

%%麦克斯韦张量法计算电机转矩
clear all
clc
load torque_240716.mat
R=0.037; %气隙半径（m）                 
lef=0.125;%电机轴向长度（m）                    
u0=4*pi*1e-7;%真空磁导率              
n=40;%气隙谐波总阶数                           
m=1;                          
Tk=zeros(n,m);             
Tdeg=zeros(n,m);

j=1;
   for i=1:n    
    Tdeg(i,j)=-1*cos(Brdeg(i,j)-Btdeg(i,j));
   %各阶次径向和切向磁密角度余弦值
    Tk(i,j)=pi*R*R*lef/u0*Brmag(i,j)*Btmag(i,j)*Tdeg(i,j);
   end
   %各阶次径向和切向磁密相乘得到力密度
Te=zeros(40,1);
for i=1:40
    Te(i)=sum(Tk(i,:));
 end
bar(Te);

以下是运行结果。以12槽10极表贴式电机为例，X坐标表示气隙谐波阶次，但值得注意的是从0阶次开始，也就是说X=6，即5阶气隙磁密谐波所贡献的转矩。以这台电机为例，该电机工作谐波为5阶，贡献转矩约为16.37Nm。这也印证了，SPM的转矩贡献由单一阶次谐波产生有效转矩的结论！

但麦克斯韦张量法的应用潜力并非局限，它更应该使用在多工作谐波电机的分析中，如永磁游标电机、磁阻电机等的！当然，也欢迎广大电机电磁科研民工基于此代码继续完善！

如有永磁电机问题可以关注下我的B站账号！！！B站搜：鹤晓枫 https://www.bilibili.com/video/BV1vd4y1W74g/