永磁同步电机多物理场仿真案例：电磁、谐响应与噪声分析，适合学习-编程实验室

永磁同步电机多物理场仿真案例，电磁—谐响应—噪声分析（NVH分析），该案例可以用于学习，具体参数见第一张图

咱们今天来整点硬核的——手把手拆解永磁同步电机的电磁振动噪声联合仿真。这玩意儿看着参数一大堆（具体数值请直接瞅图1），但说白了就是玩转三个场：电磁场抖机灵、结构场蹦迪、声学场唱歌。下面直接上干货！

一、电磁场先整活儿

电机不转起来的时候，磁场就像死宅一样瘫在铁芯里。用Python+FEMM搞个二维瞬态场仿真，核心代码长这样：

import femm femm.openfemm() femm.newdocument(0) femm.mi_probdef(0, 'millimeters', 'planar', 1e-8, 0, 30) femm.mi_drawrectangle(0,0,15,30) femm.mi_addmaterial('NdFeB', 1.05, 1.05, 0, 0, 0.5, 0, 0, 1, 0, 0, 0) femm.mi_addcircprop('PhaseA', 10, 1) # 10A电流走起 # 重点来了！设置运动边界 femm.mi_addboundprop('rotor_bdry', 0, 0, 0, 0, 0, 0, 0, 0, 5)

跑完仿真别急着关，用femm.mo_getcircuitvalues('PhaseA')抓取电磁力波形。这时候要注意看齿槽转矩的波动——这货就是后续振动的罪魁祸首！

二、结构场开始蹦跶

把电磁力加载到定子齿上，用ANSYS APDL玩谐响应分析时，记得这波骚操作：

! 导入电磁力数据 *DIM,force_data,TABLE,6,1000,,TIME *TREAD,force_data,force.txt ! 玩模态叠加法 ANTYPE,HARMIC HROPT,MSUP,,YES HARFRQ,0,4000 ! 扫到4kHz足够抓噪声 KBC,1 ! 阶梯加载不翻车 ! 关键设置：阻尼比给个玄学值0.03 DMPRAT,0.03

这里有个坑：网格要是太糙，共振峰会跑偏。建议在齿部用2mm的六面体网格，其他部位可以摸鱼用4mm。算完记得导出*.rst文件，后面声学场等着用呢。

三、声学场开演唱会

把结构振动速度映射到声学网格，用LMS Virtual.Lab干这事时，注意这个神操作：

vibData = importdata('vibration.csv'); mesh = ansysReadMesh('model.cdb'); % 玩声学逆向-NLH actool('set', 'MaxFreq', 4000); actool('assign', 'Velocity', vibData); actool('solve', 'BEM'); % 在驾驶员耳朵位置插值 fieldPoint = [1200, 500, 800]; spl = actool('getSPL', fieldPoint); plot(spl(:,1), spl(:,2), 'LineWidth',2);

重点看800Hz和2400Hz这两个频段，通常这里会出现啸叫峰值。要是声压级超过65dB(A)，赶紧回去检查电磁力的空间分布——八成是某个磁极的磁场跑偏了。