c++ 提取傅里叶描述子_探讨一下一种仿真中电机谐波特性的提取方法。

芬兰读博太无聊了，这两天突然发现知乎是个交流一下学术的好地方，就把之前自己研究的一些东西拿出来分享一下，希望能够抛砖引玉。同时也相当于给自己做个阶段性的总结。

说明：内容大部分都是自己论文里面的东西，文章档次和知乎大佬比真的很低大家不要鄙视我，我很心虚的。没任何论文中的原图片，因为有点担心IEEE图片版权问题，不敢乱用。查了一下IEEE版权政策，似乎这种blog式的重复图片发表也要申请版权付钱，一张22欧的样子，我真的是付不起。目前使用的仿真软件是Altair Flux，这里给出的方法应该都适用于其它主流电机仿真软件如：Maxwell（确定可以），Jmag(也许吧，猜的)。

---

目的：为什么要提出这个问题或者说为什么有这个需求？因为，电机中的谐波会产生一系列有意思的现象，比如附加损耗，转矩波动，磁拉力，震动噪声等。但似乎很难分离出某一种谐波造成的影响具体有多大。大家在做仿真时候一直有这个疑问，因此，今天我们就来解决一下这个问题。

研究对象：为了说明这个解决上述问题，仅以一台实心转子异步电机为例，提取出气隙中的某一次谐波在转子中产生的损耗。同理，是该方法可以求出其它谐波特性的，比如转矩，磁拉力等。这里仅已转子谐波涡流损耗为例，其实仿真中我还提取了谐波转矩。该实心转子主要参数为 2 MW, 12000 r/m, 2p, 660 V，具体参数和详细设计见[1]-[2]。正常的瞬态场计算得到在转差率

时，总的转子涡流损耗为约13.3 kW（计及所有基波磁场和谐波磁场）。

谐波磁场的提取：感应电机气隙谐波磁场可以表示为

显然，该函数为一时空变换函数，简单点说，改函数随时间

和空间

同时改变。工程上，经常使用快速傅里叶分解(FFT)来提取谐波分量，但是传统1-D FFT仅适用于一元函数，而上述公式是有两个参变量。因此，采用2-D FFT来提取不同阶次以及频率下的磁密。有两种方法，第一种，按时间和空间分别做两次FFT，第二种，直接用Matlab中的“FFT2”命令（文章中用的就是“FFT2”），其结果见[2]。

二维傅里叶分解的必要性：因为传统1-D FFT局限性，表现为有谐波相互掺杂。比如，在采用传统空间傅里叶分解时候，5次相带谐波（阶次5， 频率

）和5次饱和谐波（阶次5， 频率

）是无法辨别的，因为阶次都是5。同理，基波（阶次1， 频率

）和定子一阶齿谐波（阶次35或者37， 频率

），是无法采用传统沿时间轴的傅里叶分解来区分的。

构建混合正弦磁场激励：上述过程已经成功提取气隙磁场磁密，并包含有气隙谐波磁场的幅值、阶次和频率。因此，理论上能够反向构建出该磁密。下图就是反向构建谐波磁场的模型。当然，构建的方法有很多种，这只是其中一种可能性。可以看出，定子区域完全被省略，只有一个很小的环形区域来提供正弦磁场。该区域其实是一块永磁体，磁导率很低（保证充磁不衰减），正弦充磁能够提供足够稳定的正弦磁场。

实心转子含正弦磁场激励（有点丑，不要怪我）

关键问题来了为何要混合充磁？什么叫混合充磁？这里的混合充磁是指，上图中的环形永磁体同时提供基波以及谐波磁场。原因如下，下面的式子为谐波引起的涡流在转子表面引起的透入深度，其中有一项很容易被忽略，就是该透入深度会随转子磁导率变换而发生变化。显然，电机空载和带载时候的转子材料是工作在

曲线上的不同区域的。因此，为了模拟带载时候的转子材料特性，充磁时候还必须提供基波。这样基波和谐波混合充磁，才能反映出该谐波的真实工作环境。

谐波损耗提取：上文已经描述出了所构建的模型，要求取谐波损耗只需要，构建一个只含基波的模型，求得转子涡流损耗

，再构建一个同时含有基波和谐波的模型，求得转子涡流损耗

。其差值

就是该谐波在这一工况下所产生的在转子表面所产生的谐波涡流损耗。

局限性：文章中还用了别的计算逻辑来相互证明了方法的可靠性。唯一的缺点就是，提出的模型彻底改变了磁路，定子磁路彻底消失了，导致定子饱和无法考虑进去。不过好在对于实心转子电机，转子通常设计在尽可能饱和或者过饱和状态，这样能提高功率因数[3]。论文中的样机也是，转子是饱和的，并且定子外径是设计得很大的以减少损耗。因此，对于该样机，大部分饱和磁场是转子产生。此外，定转子饱和所产生的谐波的频率和阶次是一模一样的， 因此2-D FFT也是无法区分开来的。

展望：该模型似乎还能用来做其它用途，比如，求取某谐波产生的转矩波动，不平衡磁拉力，负序磁场产生的损耗等。

补充：感谢楼下king的提问。关于谐波磁场，对转子饱和的影响。这个点提得非常好。其实在论文里面我还建立了一个模型，将所有谐波都构建进去。假设损耗

是含有所有谐波的转子损耗，

是含有所有谐波但缺少了一项谐波（就是你要分析的那个谐波）所产生的损耗。其差值

就是该谐波在转子表面引起的涡流损耗。两种思路在论文里都考虑了，研究发现谐波损耗差别很小。原因是，谐波幅值太小了，根本无法明显改变转子材料的工作点。也就是说谐波对转子饱和影响很小。修改于4.30.2019。

---

如有问题请大家批评指正！

参考文献：

1. Di, C., Petrov, I., & Pyrhönen, J. J. (2018). Modeling and mitigation of rotor eddy-current losses in high-speed solid-rotor induction machines by a virtual permanent magnet harmonic machine.IEEE Transactions on Magnetics,54(12), 1-12.
2. Di, C., Petrov, I., & Pyrhönen, J. J. (2019). Extraction of Rotor Eddy-Current Harmonic Losses in High-Speed Solid-Rotor Induction Machines by an Improved Virtual Permanent Magnet Harmonic Machine Model.IEEE Access,7, 27746-27755.
3. Pyrhonen, J., Nerg, J., Kurronen, P., & Lauber, U. (2010). High-speed high-output solid-rotor induction-motor technology for gas compression.IEEE Transactions on Industrial Electronics,57(1), 272-280.