4. 电机气隙磁场解析与仿真对比

气隙磁场是电磁振动与噪声的“源头信号”。定子电流在气隙中建立旋转磁动势，转子永磁体或励磁绕组再叠加自身磁动势，两者在气隙这一薄层空间内合成，形成随时间和空间同时变化的磁场。这个磁场沿气隙圆周的分布并不光滑，它包含基波与一系列空间谐波；同时由于齿槽、绕组分布、逆变器开关等因素，还叠加了时间谐波。径向电磁力密度近似正比于气隙磁密平方，因此磁密中哪怕幅值很小的谐波，也可能在特定阶次和频率上激起明显的定子振动。本节围绕空载与负载两种工况，说明磁场分布差异、谐波分解方法，以及解析法与有限元结果如何相互验证。

一、空载与负载磁场分布的本质差异

空载时，电枢电流近似为零，气隙磁场主要由转子磁动势建立。以表贴式永磁同步电机为例，若忽略齿槽，理想空载气隙磁密沿圆周接近矩形波，其空间谐波只含奇次项，极对数为 p 的基波幅值最高，3 次、5 次等谐波依次衰减。实际定子开槽后，齿槽磁导变化会对该矩形波产生“齿槽调制”，在空间谐波中引入与槽数相关的边带分量，典型阶次为 |ν ± kQ|（ν 为原谐波次数，Q 为槽数，k 为整数）。这些边带谐波阶次低、幅值虽小，却是齿槽转矩和低阶径向力的重要来源。

负载时，电枢绕组通入三相对称电流，产生与转子同步旋转的电枢磁动势。电枢反应磁场与转子磁场在气隙中矢量叠加：一方面基波合成磁密幅值和相位随负载电流变化，出现增磁或去磁；另一方面，电枢磁动势本身含有 5、7、11、13 等空间谐波，与转子谐波交叉调制后，会产生新的空间阶次。若采用 PWM 供电，电流中还含有开关频率及其边带的时间谐波，它们与基波磁场调制后形成频率为 f_sw ± kf₁ 的力波分量。因此，负载工况的磁密频谱通常比空载“更密”，低阶力波成分也更丰富。

二、谐波分解方法：从波形到阶次—频率谱

拿到气隙磁密沿圆周的分布后，第一步是空间谐波分解。在某一固定时刻 t₀，沿气隙圆周取 N 个等间距采样点，对 B(θ, t₀) 做离散傅里叶变换，得到各空间阶次 ν 的幅值和相位。工程上常取一个极距或整个圆周为一个周期，采样点数建议不少于 2 倍最高关心阶次，例如关心到 50 阶时，圆周采样不少于 256 点。分解结果中，ν = p 为基波，ν = 0 为平均分量，低阶 ν = 1、2 等往往对应不平衡或偏心，需要单独关注。

第二步是时间谐波分解。固定某一空间点（例如气隙中心线上某点），记录 B 随时间的变化，做 FFT 得到时间频率成分。空载时主要频率为 f₁ 及其倍数；负载 PWM 供电时，会出现 f_sw、f_sw ± f₁、f_sw ± 2f₁ 等边带。将空间阶次和时间频率配对，就得到二维谱 (ν, f)。径向力密度 P_r(θ,t) ≈ B²(θ,t)/(2μ₀)，两个空间阶次 ν₁、ν₂ 的磁密相乘，会产生阶次为 ν₁ ± ν₂ 的力波；两个时间频率 f₁、f₂ 相乘，会产生 f₁ ± f₂ 的力波频率。因此，低阶力波往往来自两个较高阶磁密的差频/差阶组合，这也是负载工况更容易激起低阶振动的原因。

三、解析法与有限元结果的相互验证

解析法和有限元法不是互相替代，而是互相校验。解析法基于子域法或等效磁路，把求解区域划分为气隙、槽、永磁体等子域，用分离变量法写出各子域磁位通解，再通过边界条件匹配系数。它的优势是能显式给出 B(θ,t) 的级数形式，便于快速扫描极弧系数、槽口宽度、电流相位等参数对特定谐波的影响。有限元法则把几何、材料非线性、外电路都纳入求解，结果更接近真实，但单次计算耗时较长。

验证时建议按以下顺序进行。第一，先对比空载线反电动势或空载气隙磁密基波幅值：解析法取理想矩形波傅里叶系数，有限元取气隙中线一周采样后 FFT，若两者基波幅值偏差在 5% 以内，说明几何与材料参数设置基本一致。第二，对比主要空间谐波阶次：例如 12 槽 10 极电机，空载磁密中应出现 ν = 5、7、17、19 等阶次，检查两种方法得到的阶次是否一致、幅值排序是否相同。第三，加载后对比电枢反应引起的基波幅值变化：解析法通过电枢磁动势与转子磁动势矢量合成计算，有限元通过瞬态场直接得到，若去磁/增磁趋势一致，说明电流相位和绕组系数设置正确。第四，对比力波频率：在定子齿顶取若干点，计算径向力密度时间波形并做 FFT，检查 f_sw ± f₁ 等特征频率是否同时出现。

四、一个可复现的对比示例

以一台 12 槽 10 极表贴式永磁电机为例，额定转速 3000 r/min，基波频率 f₁ = 250 Hz，PWM 开关频率 f_sw = 10 kHz。空载时，解析法取极弧系数 0.85，计算得基波磁密幅值约 0.82 T，5 次谐波约 0.09 T，7 次谐波约 0.05 T。有限元在气隙中线取 360 点，一个电周期内取 200 个时间步，FFT 后基波约 0.80 T，5 次约 0.085 T，7 次约 0.048 T，误差均在 5% 以内。负载时通入 q 轴电流 80 A，解析法预测基波幅值降至约 0.76 T，有限元得到约 0.74 T，去磁趋势一致。进一步对径向力密度做二维分解，发现 ν = 2、频率 500 Hz 附近的力波在负载时明显增强，对应电枢 5 次谐波与转子 7 次谐波的差阶组合，这为后续定子模态和声辐射分析提供了明确的激励源。

需要强调的是，解析法的精度高度依赖假设条件。若电机饱和严重、齿槽较深或端部效应显著，解析结果与有限元的偏差会增大，此时应以有限元为主、解析法为辅，用解析法解释趋势、用有限元确认量级。反过来，有限元网格过粗或时间步过大，也会丢失高阶谐波，因此网格在气隙处至少应保证 3 层以上，时间步应满足最高关心频率的采样要求。两者相互验证，才能让后续的电磁振动噪声计算建立在可靠的磁场数据之上。