最近在搞8极48槽永磁同步电机的振动噪声分析,我用Motor-CAD把整套流程从头到尾跑通了一个demo——从电磁模型建立、额定点工况求解,到径向力波提取、结构模态匹配,再到瀑布图和声压级结果输出,一气呵成。这套组合拳在电驱动开发里用得特别多,尤其是新能源车主驱这类高转速、高功率密度场景,NVH性能直接影响整车主观感受。无论你是电机设计工程师、刚接触NVH仿真的学生,还是想从电磁角度切入振动噪声的新手,这篇文章都值得收藏,我会把操作步骤、参数设置、结果判读和踩坑记录都交代清楚。
1. 为什么我用8极48槽做这个demo:极槽配合的NVH底子
先说结论:8极48槽并不是随便选的,它是目前车用永磁同步电机里出现频率很高的拓扑之一。极对数p=4,槽数Z=48,算下来每极每相槽数q=48/(8×3)=2,属于整数槽绕组,双层绕组为主。这个配合在电磁性能上非常均衡——反电动势正弦性好、转矩脉动小、永磁体利用率高,同时对振动噪声来说,它的力波阶次分布也比较“友好”。
1.1 极槽搭配是如何影响振动的
电机运行时,气隙磁场相互作用会产生作用在定子齿和定子轭上的电磁力,但这个力不是静止的,它在空间上有一定波形,在时间上按一定频率波动。空间上的波形可以用“力波阶次”来描述,比如二阶力波就是圆周上呈现两个波峰两个波谷的分布,四阶就是四个波峰四个波谷。时间上的波动则表现为某个频率,通常和电频率的整数倍有关。
这个力波打在定子结构上,能不能引起明显振动,取决于它和定子本身固有模态的匹配程度。定子铁芯叠片、绕组、机壳装配在一起后,会有各种模态,常见的有0阶呼吸模态、2阶椭圆模态、3阶三角形模态、4阶方形模态等。如果电磁力波的空间阶次和结构模态阶次一致,同时时间频率又落在该模态固有频率附近,就会发生强烈的机电共振,噪声会瞬间飙起来。
8极48槽的优势就在于,它的基波磁场产生的径向力波,空间阶次是8阶,远高于电机结构最容易激发的低阶模态。低阶模态频率低、阻尼小、辐射噪声效率高,高阶模态往往频率高、响应弱,很难被激励起来。所以相对2极或4极电机来说,8极电机的NVH底子天生就好一些。当然,这并不意味着万事大吉,因为齿槽效应和磁路饱和还会带来大量谐波力波,仍可能在特定转速下踩中模态频率。
1.2 为什么说48槽是一个“甜点”槽数
槽数选多少,不是拍脑袋决定的。槽数太少,绕组磁动势谐波大,力波谐波丰富,转矩脉动也难看;槽数太多,齿槽细小,嵌线工艺变难,齿部刚度还可能会下降。48槽配合8极,每个齿距对应的机械角是360°/48=7.5°,折算到电角度是7.5°×4=30°,这个尺度对削弱齿槽转矩和齿谐波力非常有利。
另外,48槽的齿谐波力波阶次高。齿槽效应会在气隙磁导上附加空间谐波,这些谐波与永磁磁场相互作用后,会产生形如Z±2p、Z±4p这类高阶力波,也就是48-8=40阶、48+8=56阶等。阶次这么高的力波,对应的是波长很短、圆周上起伏很密的载荷分布,这样的载荷很难有效激发大型结构的低阶模态,所以在NVH分析里它们往往不是主要矛盾。真正要盯紧的,还是8阶、16阶、24阶这些相对低阶的力波,以及它们与模态频率的交叉点。
我在demo里设定的工况是额定转速6000rpm,对应电频率f_e=n×p/60=400Hz。径向电磁力基频为2倍电频率,也就是800Hz,这个频率在中小型永磁电机的定子2阶或3阶模态频率范围内很常见,所以这个转速点非常有分析价值。后文我会反复用这个数值举例,建议大家先记下来。
2. 先把原理讲明白:电磁径向力与结构模态是怎么耦合的
振动噪声分析如果只停留在“软件里点几个按钮出几张图”的层面,遇到异常结果会非常被动。我建议所有做仿真的,先把底层物理图像建立起来,后面调试demo才知道每一步在算什么。
2.1 径向电磁力的来源与计算
永磁同步电机气隙里的磁场,可以分解成径向分量与切向分量。定子齿面受到的电磁力,本质上是麦克斯韦应力张量在界面上的积分效果。工程上常用简化公式估算径向力密度:
f_r = (B_r² - B_t²) / (2μ₀)
其中B_r是气隙磁密的径向分量,B_t是切向分量,μ₀是真空磁导率。对大多数表面式或内置式永磁电机来说,B_r占主导,B_t相对较小,尤其在空载或轻载工况下,径向力密度正比于B_r²。这意味着,如果气隙磁密基波幅值是1T量级,力密度大约在几十kPa到上百kPa,这个力不算小,足以激励出明显振动。
关键是要理解,B_r不是纯净的正弦波。它包含基波极对数p=4的分量,也包含3次、5次、7次等高次谐波,还包含因定子开槽引起的齿槽谐波。力密度是磁密的平方,两个不同极对数的磁场分量相乘后,会得到空间阶次为两者之和或差的力波。比如基波磁场(4对极)与自身作用,产生8阶空间力波,时间频率上等于2f_e;基波与5次谐波(20对极)作用,产生16阶或24阶力波等。
Motor-CAD在做电磁瞬态求解时,会直接输出气隙中心线上的径向力密度波形,再通过傅里叶分解得到各阶力波的幅值和频率。这一步是整条NVH链路的数据源头,力波算不准,后边模态和声学结果全白搭。
2.2 模态匹配:NVH的“雷区”在哪
我把模态匹配比作一个粗糙的“搭积木”游戏。电磁力波是一堆不同空间阶次、不同频率的激励源,定子结构是一堆不同阶次、不同固有频率的振动系统。二者要在“空间阶次相同”和“时间频率接近”两个条件同时满足时,才会产生明显的共振。
空间阶次匹配这一点常被忽视。很多人只盯着频率谱,看到某个频率的力波幅值大就紧张,其实如果这个力波的空间阶次是40阶,而结构在对应频率处的40阶模态振型很“僵”,实际响应并不大。反过来,2阶或4阶这样的低阶力波,哪怕幅值并不算大,只要频率碰上2阶椭圆模态或0阶呼吸模态,振动响应会成倍放大,噪声瞬间就出来了。
所以我在看demo结果时,第一件事不是看声压级数字,而是看力波阶次表里哪些阶次和定子模态阶次“对得上号”。通常在8极48槽电机里,最值得关注的是2阶、4阶、0阶这几种结构模态,因为它们固有频率低、容易被激励。而对应8阶空间力波虽然幅值大,但由于对应的是结构8阶模态,频率通常很高,很多时候反而不是主要噪声源。这一步想清楚了,后边的优化方向也就清晰了。
3. Motor-CAD实操:从电磁模型到NVH结果一步步做出来
下面进入正题,我把这次demo在Motor-CAD里的操作流程完整捋一遍。我用的是Motor-CAD新版本,界面布局和旧版有些差异,但核心模块基本一致,大家可以参照自己的版本来对应。
3.1 几何参数与电磁模型配置
新建工程后先进入EMag模块,选择永磁同步电机拓扑,这里有两种常用路线:表贴式永磁(SPM)和内置式永磁(IPM)。如果做高速车用电驱demo,建议直接选IPM,V形或双层V形磁钢结构更贴近实际;如果只是想快速验证NVH流程,SPM更快,参数更少。我这次用的是内置V形磁钢结构,因为这类电机的磁阻转矩占比高,负载工况下的磁场饱和特性更复杂,分析价值更大。
关键几何参数我整理成一张表,大家在软件里照着填就能搭出一个能跑的基础模型:
- 定子外径:150mm
- 定子内径:92mm
- 气隙长度:0.8mm
- 转子外径:90.4mm
- 轴向长度:80mm
- 槽数:48
- 极数:8
- 绕组形式:双层整距叠绕组
- 每槽导体数:按反电动势目标反算
- 磁钢材料:钕铁硼N42SH
- 硅钢片:35JN250(或同类无取向硅钢)
这里气隙大小别拍脑袋设。气隙小了,磁密高,转矩密度大,但力波也会变大,而且装配工艺难度上升;气隙大了,NVH表现好,但永磁体用量和体积都上去了。0.8mm对于这个尺寸级别算是个折中方案,实际工程产品经常在0.5mm到1.0mm之间调。
绕组设置里要注意是“整距”还是“短距”,以及并联支路数。整数槽绕组常用短距来削弱5次和7次谐波,我在demo里选的是线圈跨距为9槽(总节距12槽),这样可以较好地抑制5、7次空间谐波,对力波频谱有直接改善效果。并联支路数设成1即可,省事,不会影响NVH结论。
电磁模型的激励方式,我建议用“电流激励”而不是“电压激励”。在Motor-CAD里可以直接给额定电流幅值、电流角,也可以选用自带的MTPA控制角计算。我习惯先用Id=0初算一遍,再看电流角对力波的影响。额定工况可以先设为峰值相电流250A(对应输出约100kW级),母线电压320V左右。这组参数不算精确匹配某个具体产品,但跑NVH趋势完全够用。
3.2 电磁场求解与径向力数据导出
模型搭好后,先跑一次空载,看反电动势波形和齿槽转矩;没问题后再跑负载瞬态。时间步长按电周期等分,我一般一个电周期取200到300步,太少的话力波高频分量会失真,太多则计算时间显著增加。对于8极、6000rpm的工况,电频率400Hz,一个电周期2.5ms,200步对应12.5μs一步,仿真0.2s就覆盖很多个电周期了。
求解完成后,关键后处理有两步。第一步是查看气隙磁密波形,确认基波幅值和波形畸变程度;第二步是提取气隙径向力密度波形。在Motor-CAD的后处理里,径向力密度曲线可以直接输出,软件内部已经完成了磁场平方和空间分布的处理。导出时要注意选择“相对于转子的角度”还是“相对于定子的角度”,如果暂态分析中网格动起来了,正确设置参考坐标系才能得到稳定的力波阶次。
力波的傅里叶分解结果建议整理成表格,重点看空间阶次0到12、时间频率在0到3000Hz范围内的分量。我这次算出来的主阶次大致是8阶@800Hz(基波磁场平方项)、16阶@1600Hz(基波与3次谐波互作用)、24阶@2400Hz,还有若干齿谐波相关的高阶分量。这些数据和理论预判基本吻合,说明模型参数没设错。
这里要特别提醒一点,电磁力波的时间采样率和FFT点数要足够。如果力波数据只包含少数几个电周期,频率分辨率会很低,和模态频率做匹配时会漏掉关键交点。我通常把FFT窗口设成至少16个电周期,频率分辨率能到25Hz左右,在工程预判阶段足够了。
3.3 NVH模块设置与结构参数输入
电磁结果算完,切换到NVH模块。Motor-CAD的NVH模块支持多种工作模式,最简单的做法是使用软件内置的2D结构模态求解器,对定子铁芯、绕组、机壳组成的截面做平面应力分析,得到各阶模态频率和振型。虽然2D模型不考虑轴向模态,但在电机轴向长度远大于截面尺寸时,2D分析能得到很不错的近似结果,特别适合demo阶段快速扫雷。
结构参数里,有几个地方特别影响模态计算结果:
- 定子铁芯材料属性:杨氏模量、密度、泊松比,注意叠片系数。铁芯是叠片结构,实际刚度比实体块要低,在软件里输入等效弹性模量时,一般要打9折甚至8折。
- 绕组质量:铜线、绝缘漆、槽楔合在一起,质量占定子总质量的比例相当大。NVH模块里可以设置槽内绕组填充系数,以及绕组端部质量如何处理。我是按填充系数0.45计算的,端部质量作为一个附加质量环加在铁芯两端。
- 机壳与螺栓连接:如果是自然冷却或水冷机壳,机壳和定子外圆的接触刚度很关键。过盈量大的话,两者近似刚性连接,模态频率会明显抬高;过盈量小或存在间隙的话,模态频率下降且阻尼上升。demo阶段可以分两种情况各算一次,看模态频率对边界条件的敏感性。
模态算完,NVH模块会自动把前面提取的电磁力波映射到结构网格表面,然后求解每个频率下的振动速度响应。这里要注意,力波的加载位置要选在定子齿面,而不是气隙中某根虚构的线。有的版本会默认把力加载到定子内圆上,这是对的,但如果模型里存在独立的气隙单元,要确认没有加载到气隙节点上。
完成振动响应计算后,软件再基于表面振动速度和声辐射效率,计算辐射声功率和指定测点的声压级。Motor-CAD里可以直接生成声压级的频率响应曲线,也可以选择在不同工况下组合成坎贝尔图。
4. 结果怎么看:阶次图、坎贝尔图与声学评价
跑完NVH计算,软件会吐出大量曲线和云图,新手很容易迷失。我讲一下我一般按什么顺序看结果、重点看什么。
4.1 阶次图和瀑布图里找共振点
阶次图的横轴通常是转速,纵轴是频率,颜色深浅代表振动幅值或声压级大小。电磁激励源有一个特点:它们的频率与转速严格成正比,所以在阶次图上会表现为从原点出发的一条条射线。比如2倍电频率激励表现为“2阶”的斜线,齿槽谐波表现为“12阶”甚至更高阶次的斜线(因为齿槽转矩周期数在8极48槽电机里是48,对应频率是48×f_r,也就是12×f_e)。这些斜线就是电磁力波的“轨迹”。
结构模态则表现为水平方向的亮带,因为模态频率基本不随转速变化。当斜线穿过水平亮带时,就是共振点。在瀑布图(三维频谱图)上,你会看到这些交叉点附近幅值骤然放大,形成类似“鲸鱼尾巴”或“山峰”的形状。这就是我们在仿真中要找的关键风险点。
我在这个demo里重点关注了800Hz附近——这是2倍电频率在6000rpm时的值。如果这个频率点落在结构2阶模态频率附近,那么在额定转速附近就会出现明显共振,噪声能量会集中在800Hz左右。解决思路无非是错开频率、降低力波幅值、增加阻尼三条路,后面我会具体展开。
4.2 声学指标怎么评价
声学结果包括声功率级(SWL)和声压级(SPL)。声压级要注意是否采用A计权,A计权模拟人耳对中低频不敏感、对2kHz到5kHz更敏感的特性,NVH性能评价中更常用A计权声压级。声压基准是20μPa,声功率基准是1×10⁻¹²W,许多新手容易在这两个基准上犯错。
对于电驱动总成,行业里常见的NVH目标是:在常用工况下,A计权声压级控制在某个限值内,同时避免在某一转速点出现尖锐的单频啸叫。单频成分哪怕总声压级不高,主观感受也会很差,这就是为什么除了看总声压级,还要看频谱里有没有明显的离散峰值。我在demo中特意做了6000rpm单点声压频谱,结果在800Hz处出现一个显著峰值,对应前面说的径向力基频与模态频率靠近的问题。
这里引入一个工程经验:判断是否需要优化,不能只看数值大小,还要看这个峰值离“听觉敏感区”和“共振区”有多近。如果峰值频率落在2kHz到5kHz的人耳敏感频段,哪怕声压级只有60dB(A)也可能被吐槽“啸叫”。反过来,如果峰值在200Hz以下,人耳不敏感,容忍度就高很多。
5. 调试过程中踩过的坑与解决思路
最后这部分是这次demo里最有价值的记录,很多问题不是看文档能提前预判的,只有真正动手跑一遍才会碰上。
5.1 常见报错与设置踩雷
先说说几个容易让人摸不着头脑的坑:
第一,力波频率范围不够导致模态匹配漏检。我第一次跑NVH时,瞬态仿真只算了20个电周期,FFT频率分辨率很粗,结果6000rpm附近的共振峰没有明显显示出来,我还以为是电机NVH表现不错。后来加长了瞬态仿真时间,峰值才暴露出来。结论是,做NVH用的电磁瞬态数据,时长一定要足够,频率分辨率至少做到25Hz以下,否则高频窄带共振会被平均掉。
第二,模态频率异常偏高或偏低。出现这种情况,先查材料弹性模量是否按各向同性输入。定子铁芯叠片在轴向和径向的模量差异很大,如果统一按实体钢的210GPa输入,模态频率能高出20%。另一个容易忽略的是绕组端部质量,漏掉端部后结构质量偏低,模态频率整体偏高,共振点判断就会偏移。我后来在模型里按照绕组总质量分摊到定子端部两侧,模态频率立刻降下来不少。
第三,2D模型与3D模型的偏差。2D模态分析不考虑轴向振型和端部效应,对轴向长度较短、机壳刚度较大的电机误差小;但如果你碰到的是轴向长、机壳薄的结构,2D结果只能用来做趋势判断。demo阶段可以先跑2D筛选出风险转速和风险阶次,真正出校核报告前再用3D模型验算。
第四,声压级结果出现负值或异常低值。检查测点距离和声辐射模型的设置,如果测点离电机太远,低频段声压会衰减到很低;如果边界条件里设置了吸声边界,也会影响结果。验算时可以把测点放在1m或0.5m这样的标准距离,方便和实测对标。
5.2 优化方向速览:从力波源头和结构路径两头堵
如果共振风险确认了,优化思路大体分两类:一类是降低激励源幅值,另一类是改变结构响应。
激励源侧最常见的工程手段是斜槽或斜极。8极48槽电机斜一个定子齿距(7.5°机械角,对应30°电角度),可以显著削弱齿槽转矩和齿谐波力波。注意斜槽后反电动势会略微下降,输出转矩会有一定损失,实际设计需要在NVH和性能之间取舍。分段斜极(比如分成三段,每段偏转一定角度)比连续斜槽更容易制造,效果也接近。另一个激励源侧手段是优化电流控制策略,比如调整电流角让某次谐波力波幅值最小,或者在PWM策略里避开某些开关频率边带。
结构侧优化包括增加定子轭厚度提高刚度、优化机壳加强筋布局、填充阻尼材料增加损耗因子等。这些措施改造成本高,通常放在产品后期调校阶段。
我在这次demo里用的一个快速验证方法非常直观:在NVH模块里分别关闭不同阶次的力波激励,看声压级贡献度排序。结果发现800Hz附近的峰值主要由8阶力波贡献,于是我把关注点放在了削弱基波磁场的谐波畸变上,比如调整磁钢厚度和极弧系数。仅调整极弧系数一项,800Hz处声压级就下降了约3dB(A),这在实际产品上已经是一个可感知的改善了。
5.3 结果验证与对标建议
仿真做完了,最好能用实测数据对标一次,不然永远不知道自己的设置是不是靠谱。最经济的验证方式是做一台简易样机,在台架上用加速度传感器测振动,用传声器测噪声,测点和仿真测点保持一致,记录相同工况下的频谱。
实测和仿真通常会有20%以内的频率偏差,这很正常,不用急着改模型。优先检查模态频率对么,比如锤击模态实验或阻抗分析,如果结构模态实测和仿真差5%以上,先修正结构参数;如果模态对准了,再对比响应峰值对应的转速点。一般来说,共振转速点对得准,就可以判断模型可信。
我个人有一个习惯,每次跑完demo都会把力波阶次表、定子模态频率、共振转速点三张表一起贴到项目记录里,方便后续调整参数后快速对比。这个方法看起来土,但在项目密集时非常提效。
最后再分享一个小技巧:做参数化扫描时,不要一开始就搞全因子设计。先用Motor-CAD的单变量扫描,把气隙、极弧系数、斜槽角度这几个影响最大的参数各拉一个趋势,找到敏感方向后再做多变量联合优化。我这次demo只扫了极弧系数和磁钢厚度两个参数,就基本锁定了优化方向,计算量比盲目组合少太多。
这套从电磁到振动再到声学的完整链路,现在在前期设计阶段就能跑通,对项目决策帮助巨大。你拿到一个电机方案,先在demo阶段摸清NVH底数,再决定要不要投样机,省下的时间和成本非常可观。