ANSYS 电机受力仿真：8 极 IPM 电机径向电磁力密度波形与空间阶次频谱

项目过程

本项目为某驱动电机 NVH 团队的电磁力评估项目，需求是量化 8 极 IPM 永磁同步电机在额定工况下的气隙径向电磁力密度 pr，并做空间阶次傅里叶分解，定位激励定子齿部 NVH 噪声的主谐波阶次。客户给定：8 极 48 槽 IPM 永磁电机、额定转速 4000 rpm、额定转矩 80 N·m、定子外径 200 mm、气隙长度 0.8 mm、磁钢 N35SH。

本项目采用 ANSYS Maxwell 2D 静磁场 + 自研 Python 后处理。Maxwell 模型取一对极（机械 45°），5 万单元自适应网格（局部加密到气隙 0.05 mm），激励为三相正弦电流（峰值 200 A，电气角 30°）。后处理：从气隙中线 360 个等角度点提取径向磁密 Br 与切向磁密 Bt，按 Maxwell 应力张量 pr=B²r/2μ0 计算径向力密度，再做 FFT 空间阶次分解（每对极 4 阶单位，0~24 次谐波）。

案例分析

径向力密度 pr 沿气隙一周呈现 8 个明显波瓣（峰谷交替），峰幅 pr_max≈410 kPa，谷值接近 0。8 个波瓣恰好对应 8 极磁场的”8 阶空间基波”，位置与磁极中心对齐，与理论预期完全吻合。波形并非完美正弦——峰顶部出现”双驼峰”，这是磁钢饱和 + 定子开槽引起的非正弦畸变。

FFT 阶次分解显示主谐波为 n=8（200 kPa），即 8 极磁场基波；n=16 谐波约 70 kPa（齿槽谐波），n=24 谐波约 22 kPa（磁极斜槽或偏心引起），更高次谐波 <10 kPa。n=0 平均力 185 kPa 是稳态径向拉力，不引起振动。n=8 主谐波的振动激励频率 = 转速 × 极对数 = 4000/60 × 4 = 266.7 Hz，与客户 NVH 试验中识别的 268 Hz 电磁噪声一致。

总结教训

本项目最容易踩的坑是 Maxwell 应力张量的”虚功位移”误用。客户最初用 Maxwell 自动虚功法计算电磁力，结果 N·m 转矩数值偏差 8%，后处理 pr 偏差更大；本项目严格按 pr=B²r/2μ0 公式提取，与经验一致。

另一坑是 FFT 窗口选择。空间阶次 048 是必须的，但客户最初只取 016，丢失了 n=24+ 高次谐波；这些高次虽然幅值小，但可能与定子高频模态（>5000 Hz）共振。

项目经验

经验一：n=8 主谐波幅值 200 kPa 是定子齿部电磁噪声的核心指标。本电机定子 48 槽，对应齿部”齿响应函数”在 8 阶最大放大倍数 1.6，实际定子齿振动位移约 0.8 μm（NVH 限值 1.0 μm 内）。如果客户需要降噪 20%，可考虑磁极斜槽 1 个槽距或转子开辅助槽。

经验二：n=16、n=24 高次谐波幅值较小（70、22 kPa），但频率更高——n=24 对应 f=800 Hz，可能与机壳 1 阶模态共振，需要单独评估。

经验三：客户最初考虑用 12 极 72 槽替代方案，理论上 n=12 阶波幅值会下降，但本项目证明齿数增加也带来齿槽谐波，反而得不偿失。

有限元仿真