前阵子接了个预研需求:用 Motor CAD 对一台 8 极 48 槽永磁同步电机做振动噪声分析。说实话,刚开始我以为这只是个常规的“电磁方案 + 简单模态”的组合,但真正动手之后才发现,要想把电磁、结构、声学三个环节串成一条可靠的计算链路,中间有不少细节值得专门记录下来。这篇文章就是这次 demo 的完整复盘,从极槽配合关系讲到声功率级输出,把每一步的思考过程和踩坑经验都写出来。
如果你正在做电机本体设计、NVH 预研,或者准备用 Motor CAD 评估电机的声品质,这篇文章应该对你有用。它不会教你所有按钮的位置,但会让你知道为什么这样做、有哪些隐藏风险,以及怎么快速判断结果可不可信。
1. 先搞清楚“8级48槽”意味着什么
1.1 整数槽绕组的底子和极槽配合
拿到“8 级 48 槽”这个组合时,多数电机工程师心里会快速算一个数:每极每相槽数
[ q = \frac{Z}{2 p m} = \frac{48}{2 \times 4 \times 3} = 2 ]
整数槽,而且是双层分布绕组的典型方案。这个“整数槽”身份太重要了,它直接决定了电磁力波的底子怎么样。
对比一下常见的分数槽集中绕组,比如 10 极 12 槽,(q \approx 0.4),绕组端部短、齿槽转矩小,但代价是气隙磁场里会出现大量低阶空间力波,尤其容易激发出 0 阶呼吸模态和 2 阶椭圆模态,NVH 天然吃亏。而 8 极 48 槽这种整数槽方案,绕组对称性好、反电动势波形更干净,低阶力波相对少,理论上振动噪声的底子是更好的。
但“更好”不代表“没问题”。整数槽电机里最典型的 NVH 敌人是齿谐波力波和 2p 阶力波。8 极电机的极对数是 4,基波磁场相互作用后会产生空间阶次 (2p=8) 的径向力波。8 阶模态在定子结构上属于较高阶模态,固有频率通常不低,但一旦转速范围拉得宽,或者变频器谐波丰富,仍然可能撞上。
1.2 电磁力波的阶次与频率是振动噪声的激励源
永磁同步电机的振动噪声本质上是“电磁力波 + 结构模态 + 声学辐射”三者的耦合问题。激励源是作用在定子齿面上的径向电磁力密度,工程上最常用麦克斯韦应力张量法算:
[ F_r \approx \frac{B_r^2}{2 \mu_0} ]
也就是气隙径向磁密的平方除以两倍真空磁导率。这一步在 Motor CAD 里是自动算的,但你要理解它的物理含义,否则后处理时很容易被海量频谱图淹没。
力波有两个维度必须同时看:一个是空间阶次,一个是时间频率。空间阶次决定它会“压弯”定子的哪一种模态,时间频率决定你会在哪个转速、哪个声频点上看到峰值。典型的整数槽永磁同步电机,径向电磁力的主要时间频率是电频率的偶数倍,比如 (2f_e)、(6f_e)、(12f_e),再加上电流谐波产生的旁频。空间阶次方面,除了 0 阶和 2p 阶,还有齿谐波相互作用产生的高阶力波,比如 44 阶、52 阶这些,它们虽然空间阶次高,但对应的时间频率往往落在人耳敏感的中高频范围,会产生明显的电磁噪声。
所以做 demo 的第一步绝不是急着点“分析”,而是先在电磁设计里把气隙磁密波形、电磁力频谱和空间阶次导出来看一眼。如果激励源的频率峰值刚好能和定子某一阶模态的固有频率重合,那后面所有的结构优化都会事倍功半。反过来,如果激励源本身就很干净,结构设计只要不“主动作死”,NVH 结果大概率不会差。
2. Motor CAD 振动噪声分析的模型链路
2.1 从电磁模型到结构模型的搭建
Motor CAD 最舒服的地方是它把电磁、热、结构、声学塞在一个模型里,不需要像传统流程那样在 Maxwell 和 Mechanical 之间反复导入导出。我这次的 demo 流程是先在电磁模块把电机性能算到“基本满意”,再切到振动噪声模块。
新建一个 PMSM 工程时,注意选择正确的拓扑类型。8 极 48 槽如果做成表贴式,转子结构选 Surface PM 就行;如果做成内置式,要选 V 型或一字型,并检查磁桥尺寸。这一步不要偷懒,磁桥的饱和程度直接影响气隙磁密谐波,进而影响电磁力。
绕组设置里,48 槽 8 极双层绕组的默认节距一般是 5 槽,整距是 6 槽,实际设计常采用短距来削弱 5、7 次谐波。在 Motor CAD 里可以直接设定每槽导体数、并联支路数、绕组层数,软件会自动计算绕组系数。如果你用的是“自动绕组编辑器”,务必看一眼生成的绕组排布对不对,尤其是跨槽方向,错了的话后面一切结果都是废的。
材料分配方面,定子铁心要用非线性硅钢片 B-H 曲线,磁钢根据温度选牌号,比如 N42SH 或者 N35UH。这里有个很容易被忽视的点:dashboard 上的电磁计算用峰值电流还是 RMS 电流,一定要统一。我见过有人在电磁模块里用 RMS 电流算出了正确的扭矩,但切到振动噪声模块时默认变成了峰值电流,结果电磁力大了整整 (\sqrt{2}) 倍,谐响应峰高得离谱。
2.2 材料、约束与网格设置
切换到振动噪声模块后,Motor CAD 会根据电磁几何自动生成结构模型。但“自动生成”不等于“什么都不用管”,以下几项必须手调。
第一是材料参数。定子铁心在电磁计算里只关心导磁性能,但结构计算里需要密度、杨氏模量和泊松比。叠片铁心是各向异性的,径向和轴向的等效弹性模量不一样,不过工程简化时常用各向同性近似,密度取 7600~7800 kg/m³,弹性模量取 180~200 GPa。绕组问题更麻烦:它由铜线、绝缘漆、槽楔混合而成,结构上不是连续实体,常见做法是把它折算成等效质量块压在定子齿槽里,密度按“铜重 + 绝缘重量 / 槽体积”估算,弹性模量取较低的值,甚至设置为一种“配重材料”。如果你把绕组当成实心铜块建模,模态频率会明显偏高,因为刚度贡献被算多了。
第二是约束设置。电机在实际安装中有卧式底脚安装、端盖法兰安装、自由悬挂等状态,不同约束下的固有频率可能相差 20% 以上。做 demo 时通常先按“自由边界”算一版,看结构的本征特性;再对比一版“底部固定”,看安装刚度带来多大变化。注意约束并不是越多越好,约束过强会把共振峰频率拉得过偏,后期的改进方向也会被误导。
第三是网格。结构网格不需要像电磁 FEA 那么细,但定子齿顶面、齿根圆角、机壳散热筋这些力传递路径上的位置要局部加密。Motor CAD 的默认网格往往是为了“能算”,不是为了“算准”,我第一次跑出来的声功率级曲线带有明显的网格波纹,加密后就平滑了。建议先用粗网格跑通流程,再加密验证关键频段的结果。
第四是阻尼。结构模态阻尼比通常取 0.5%~2%,具体取决于是否考虑铁心叠片间的摩擦和总装后的连接阻尼。如果阻尼设成 0,谐响应在某些共振频率点会冲出无穷大的尖峰,数值上表现为一个窄到不真实的超高峰,这在实际结构中是不存在的。做趋势对比时,统一阻尼比比精确的阻尼值更重要。
3. 实操一遍:从电磁力到声功率级
3.1 电磁部分快速出“合格”模型
这个 demo 用的电机基本参数大致是这样:定子外径 160 mm、内径 100 mm、铁心长度 80 mm、气隙 0.8 mm,表贴式磁钢。额定转速我设为 4000 rpm,对应电频率
[ f_e = \frac{4000}{60} \times 4 = 266.7 \text{ Hz} ]
电磁力主要分量自然落在约 533 Hz((2f_e))、1600 Hz((6f_e))等位置。如果你用变频器供电,还要把载波频率附近的边频带纳入关注范围,比如 8 kHz 载波减 2 kHz 左右的边带。
电磁分析时我通常先跑磁路法快算,再切到 FEA 精确求解。Motor CAD 的 FEA 求解器在时间步长足够小的情况下,能直接给出每个时刻的气隙磁密分布。为了后续振动噪声分析,时间步长要设置得足够密,至少要覆盖最高关注频率的 1/2 或 1/3,否则高频分量会被混叠。我当时设置的步长对应最高分析频率 10 kHz,奈奎斯特频率足够,但如果你只关心 2 kHz 以下,步长可以放宽,计算时间会快很多。
这里要特别说一句:电磁网格和时间步是振动噪声结果的“上游精度”。如果你在电磁侧把谐波砍了,结构侧必然看不到对应频率的振动峰。所以做 NVH 的电磁模型,不能只验证扭矩对不对,还要单独看一下空载和负载下的气隙磁密 FFT,确认主要谐波阶次都在。
3.2 振动噪声部分:模态、谐响应、声学
电磁模型算完后,切到振动噪声模块。我的操作习惯是分三步走:
第一步是模态分析。让 Motor CAD 计算 0~10 kHz 范围内的固有频率和振型。重点关注 0 阶呼吸模态、2 阶椭圆、3 阶三角形、4 阶四边形这些低频弯曲模态。8 极 48 槽整数槽电机由于 2p=8 阶力波存在,也要留意 8 阶模态的位置,但那个频率通常偏高。把模态频率列出来,和电磁激励频谱对照,能很快找到潜在共振点。
第二步是谐响应分析。把电磁模块算出的径向和切向电磁力映射到结构网格上,施加在定子齿面,扫频方式计算定子表面的振动速度和加速度。频率步长建议不大于 5 Hz,共振峰较窄时可以局部加密到 1 Hz。别小看这一步,步长太大可能直接把共振峰“跨过去”,导致你误判这个转速没问题。
第三步是声学计算。Motor CAD 的噪声模块可以在定子周围定义声学计算边界,计算远场声压级(SPL)或声功率级(SWL)。通常建议同时看线性 dB 和 A 计权 dB,因为人耳对中频敏感,A 计权后某些低频峰会被压下去,但它们在实际测试中可能依然影响振动。
结果出来后,第一个要看的图是“振动速度 vs 频率”和“声功率级 vs 频率”的曲线。如果某个频率点的声功率级明显突出,顺着频率往回找:这个频率是不是 (2f_e) 的倍频?对应的空间力波阶次是多少?定子在这个频率附近有没有模态?三步一串联,问题定位就很快。对这台 8 极 48 槽电机,我仿真发现在 6f_e 附近有一个明显峰值,对应 1600 Hz,正好接近定子 2 阶模态的固有频率,属于典型的力波-模态耦合。
对于不同转速,建议扫一个转速范围而不是只看额定点。把转速从 1000 rpm 扫到 5000 rpm,得到“坎贝尔图”,可以直观看到共振转速带。如果某一转速区间内电磁激励频率线和模态频率线有交叉,那就在这个转速附近细化谐响应计算,确认峰值量级。
3.3 用参数扫描做方案对比
demo 阶段最大的价值不是“算出一个绝对声功率级”,而是“快速对比方案趋势”。我在这个项目里用 Motor CAD 的自动扫描功能扫了三个参数:极弧系数、磁钢偏心距和是否斜极。
以斜极为例:不做斜极时,8 极 48 槽整数槽电机的齿谐波比较强,对应高频噪声;把磁钢分段斜极或者定子斜槽后,齿谐波磁场被削弱,声功率级峰值下降明显。仿真结果里,斜极方案的总声功率级比不斜极低 3~5 dB,这在人耳听感上是能清晰分辨的改善。参数扫描时要注意固定不变项,比如额定扭矩、反电动势、效率允许范围,不能只看 NVH 一项指标盲目优化。
Motor CAD 的扫描效率比传统 3D FEA 高很多,一个转速点的电磁 + 振动 + 噪声计算几分钟就能完成,适合在概念阶段筛掉明显有 NVH 缺陷的方案。粗筛出两三个方案后,再拿详细的多物理场工具做精度复核,这是目前比较实用的工程路径。
4. 实操中的常见问题与避坑记录
4.1 力波映射和网格不一致造成的假峰值
我这次 demo 遇到最迷惑的问题是谐响应结果里出现了一排“梳齿状”的尖峰,频率间隔非常规律,像是什么周期性噪声,而不是单个共振峰。排查了很久,发现是电磁网格和结构网格不一致导致的插值误差。
具体来说,电磁 FEA 网格在气隙处把磁密算在单元中心或节点上,映射到结构网格时,如果两套网格的密度、位置相差太大,力密度分布会出现高频“锯齿”。这个锯齿被结构模型响应放大后,表现为一串虚假的窄峰。解决方法很简单:把结构网格的定子齿顶区域刻意加密到电磁网格的同等量级,或者提高电磁模型的时间分辨力,避免夸大高频分量。
做这种多物理场耦合分析,我建议先跑一个非常稀疏的网格版本,确认整体趋势没问题后,再加密网格看结果变化。如果加密后峰值位置移动超过 1% 或幅值变化超过 2~3 dB,说明结果还没有网格无关性,需要继续加密直到收敛。
4.2 模态仿真与测试对不上
另一个经典问题是仿真出的固有频率和样机锤击模态测试结果相差 10%~20%。这里面原因很典型:叠片铁心的等效弹性模量估算不准、绕组浸漆后对刚度的贡献被忽略、机壳螺栓连接刚度被高估或低估。
我自己的经验是:绕组不能只当“死重”。中等功率电机浸漆后,绕组和定子铁心之间的接触刚度相当可观,尤其槽楔灌漆后,整个绕组会参与结构刚度。处理办法是把绕组等效弹性模量从完全不加硬的“配重”逐渐提高到接近绝缘浸漆复合材料的真实水平,再多比对一两台样机标定。
另外,测试时的边界条件也很关键。自由悬挂的电机和底脚锁死的电机,模态频率完全两码事。仿真时加约束的位置、面积要参考实际安装面,不要随手加一个“固定约束”。固定约束过多会把刚体模态错误消除,导致低阶弹性模态频率虚高。
4.3 单位、频率范围和后处理的坑
振动噪声模块涉及结构和声学,单位制混用是重灾区。电机的定子尺寸你用 mm,弹性模量用 GPa,密度用 kg/m³,但计算结果里的力是 N,压力是 MPa,声压是 Pa。不同模块可能有自己的默认单位体系,切换后如果不检查单位一致性,谐响应幅值能差出几个数量级。我的习惯是流程开始前先统一单位显示,重大操作后重新核对一遍,特别看力密度和声压的量纲。
频率范围方面,只算到 2 kHz 是很常见的失误。永磁同步电机的高频电磁噪声往往落在 5~10 kHz 区间,尤其是 PWM 供电时,载波边带会产生明显的哨叫声。如果你用谐响应分析只看低频,会漏掉真正让人烦躁的噪声。把分析范围放到 10 kHz,虽然计算时间增加,但对 NVH 判断来说是必要的。
后处理的时候要注意声功率级和声压级的区别。声功率级是声源固有的总辐射能力,不随距离变化;声压级是某个测点处的实际声压,随距离衰减。很多初学者看到“某个距离处的声压级”突然低了,以为是改进有效,其实只是距离远了。我习惯用声功率级做方案对比,用声压级辅助判断实际人耳位置效果。
最后说一句关于 demo 版本本身的事。如果你用的也是试用授权,通常模型数量、网格规模、导出格式会有一定限制,结果图也难免带水印。这不影响趋势判断和原理学习,但如果你要拿数据出正式报告,还是建议用正规授权。网上不少人不清楚 demo license 怎么安装、找不到入口,我建议直接找软件原厂或正规代理商申请试用序列号,不要依赖来路不明的“绿色版”,仿真软件许可机制复杂,版本错乱、计算出错的风险很大,最后算出来的东西没人敢用。
4.4 坎贝尔图和瀑布图的读法
除了单点频率响应,坎贝尔图是判断共振转速最直观的工具。横轴是转速,纵轴是频率,颜色深浅代表振动噪声幅值。电磁激励频率会随着转速线性上升,比如 (2f_e)、(6f_e) 这些斜线;结构固有频率则近似水平。斜线和水平线相交的地方,就是潜在共振点。
读坎贝尔图时要小心颜色标尺。默认情况下颜色映射可能会把极个别大峰压艳,导致其它中等峰完全看不清。我会手动归一化色标,比如把最大值截到某个水平,再看中间频段有没有“隐性”峰值。另外,转速扫描步长不能太粗,否则共振带之间的窄带会被忽略,一般转速步长取 50~100 rpm 比较合适。
5. 这个 demo 做完后,我的一点体会
Motor CAD 的振动噪声模块非常适合设计早期做“横向对比”。它不需要精细的三维 CAD 模型,能够在几分钟内给出电磁力、模态和声辐射的整体趋势,这在概念设计阶段是极强的优势。但你也别指望它替代详细结构有限元,毕竟等效刚度和边界条件都做了简化,绝对数值只能当参考。
我个人在实际项目中的做法是:先用 Motor CAD 把 8 极 48 槽这类整数槽方案的极槽配合、斜极角度、极弧系数快速扫描一轮,筛出 NVH 底线合规的两三套方案;然后再对这几套方案做详细电磁仿真和结构模态复核,最后才出样机。这个流程可以把 80% 的打样风险在电脑前消灭掉,样机阶段集中验证制造和材料偏差,比“做了再测、测了再改”省钱得多。
这次 8 极 48 槽台电机,最终通过调整极弧系数和增加磁钢分段斜极,把主要阶次的电磁力幅值明显压了下去。当然,不同尺寸的电机表现差异很大,这个数值不能直接套到你自己的设计上,但分析思路是可复用的:先找激励源,再对比模态,最后优化解决。做完这个 demo,我最大的感触是,NVH 仿真难不在软件操作,难在你能不能把“力的频率和阶次”“结构的固有模态”“声辐射的效率”这三张图在脑子里合到一块。这三张图想明白了,剩下的只是点按钮的时间。