3、电磁力波计算：基于有限元法的电磁力波提取方法，气隙磁密的傅里叶分解，力波阶数与谐波次数的关系

好，咱们进入正题。这一节讲的是电磁力波怎么算。说实话，很多刚入行的工程师一听到“力波计算”就头大，觉得公式多、概念绕。其实没那么玄乎。我个人的习惯是，先搞清楚力波从哪来，再谈怎么算。

电磁力波，说白了就是电机定转子之间那个气隙磁场相互作用产生的力。这个力不是均匀的，它会随着空间位置和时间变化。你想想看，电机转起来的时候，气隙里的磁密在变，力自然也在变。这些变化的力作用在定子齿上，就成了振动和噪声的源头。

3.1 基于有限元法的电磁力波提取方法

怎么把力波提取出来？最直接的办法就是有限元法。我当年刚做电机设计时，还用手算过力波，那叫一个痛苦。后来用了有限元，效率高多了。

有限元法提取力波，核心就三步：

1. 建好模型——把电机的几何结构、材料属性、绕组设置都弄对。这一步马虎不得，模型不准后面全白搭。
2. 算磁场——求解麦克斯韦方程组，得到每个时刻、每个位置的气隙磁密。一般用瞬态场求解器。
3. 算力——用麦克斯韦应力张量法，把磁密转换成力。公式长这样：

σ_r = (B_r² - B_θ²) / (2μ₀)
σ_θ = (B_r * B_θ) / μ₀

其中σ_r是径向力密度，σ_θ是切向力密度。B_r和B_θ分别是径向和切向气隙磁密。μ₀是真空磁导率。

我在项目中遇到过一个问题：有限元网格太粗，导致力波计算结果波动很大。后来我把气隙区域的网格加密到0.1mm级别，结果就稳定了。嗯，这里要注意，气隙网格质量直接影响力波精度。

3.2 气隙磁密的傅里叶分解

有限元算出来的磁密是时域信号，一堆离散点。但我们需要知道它包含哪些频率成分、哪些空间阶次。这时候就要用傅里叶分解了。

傅里叶分解分两步：

• 空间傅里叶分解——沿圆周方向，把磁密分解成不同空间阶次（极对数）的分量。
• 时间傅里叶分解——随时间变化，把磁密分解成不同时间频率的分量。

实际工程中，我们通常做二维傅里叶变换，同时得到空间阶次和时间频率。代码实现大概是这样：

% 假设 B 是二维数组，行代表空间位置，列代表时间步
B_fft = fft2(B);  % 二维FFT
% 提取空间阶次和时间频率
spatial_order = (0:size(B,1)-1) - floor(size(B,1)/2);
time_order = (0:size(B,2)-1) - floor(size(B,2)/2);

分解完之后，你会看到磁密由一系列谐波组成。每个谐波都有三个特征：幅值、空间阶次、时间频率。这三个参数决定了它会产生什么样的力波。

3.3 力波阶数与谐波次数的关系

这是最核心的部分。力波不是凭空产生的，它是两个磁密谐波相互作用的结果。具体来说：

假设气隙中有两个磁密谐波：

• 谐波A：空间阶次ν₁，时间频率f₁
• 谐波B：空间阶次ν₂，时间频率f₂

它们相互作用产生的力波，其阶数和频率为：

力波空间阶次 r = |ν₁ ± ν₂|
力波时间频率 f = |f₁ ± f₂|

为什么会这样？因为力是磁密的平方，而平方运算在频域里就是卷积。两个谐波卷积，就会产生和频与差频分量。

举个例子：

• 基波（ν=1, f=50Hz）与5次谐波（ν=5, f=250Hz）相互作用
• 会产生 r = |1-5| = 4 阶力波，频率 f = |50-250| = 200Hz
• 也会产生 r = |1+5| = 6 阶力波，频率 f = |50+250| = 300Hz

你看，一个简单的相互作用，就能产生两个力波分量。

为了让你更直观地理解，我画了一张图，展示力波阶数与谐波次数的关系：

最后，我把常见的力波阶次与谐波次数的对应关系整理成了一张表，方便你查阅：

谐波A阶次ν₁	谐波B阶次ν₂	力波阶次r = |ν₁-ν₂|	力波阶次r = ν₁+ν₂	典型来源
1（基波）	1（基波）	0	2	基波自作用
1（基波）	5（5次谐波）	4	6	定子绕组谐波
1（基波）	7（7次谐波）	6	8	定子绕组谐波
1（基波）	Z±1（齿谐波）	Z 或 Z-2	Z+2	定子开槽
p（转子极对数）	p（转子极对数）	0	2p	转子磁场

这张表你最好记在心里。做电机设计时，看到某个谐波，马上就能想到它会产生哪些力波。我在项目里排查振动问题时，经常拿着这张表对照FFT结果，很快就能锁定噪声源。

好了，这一节的内容就到这。力波计算是振动噪声分析的基础，搞懂了它，后面分析模态、预测噪声就顺了。