1. 永磁同步电机NVH仿真概述
永磁同步电机（Permanent Magnet Synchronous Motor, PMSM）作为新能源汽车、工业伺服等领域的核心动力部件，其噪声、振动与声振粗糙度（Noise, Vibration and Harshness, NVH）性能直接影响产品品质。多物理场联合仿真技术通过耦合电磁场、结构场和声学场的相互作用，能够准确预测电机的振动噪声特性。
这个案例完整呈现了从电磁力计算到噪声预测的全流程，主要解决三个核心问题：
- 如何准确计算时空分布的电磁力密度
- 如何将电磁力转化为结构振动响应
- 如何将结构振动映射为可听噪声
提示：本案例采用Python+FEMM进行电磁场计算，ANSYS APDL处理结构谐响应，LMS Virtual.Lab完成声学仿真，形成完整的工具链。
2. 电磁场仿真关键实现
2.1 二维瞬态场建模要点
使用FEMM（Finite Element Method Magnetics）进行二维电磁场仿真时，需特别注意以下参数设置：
python复制femm.mi_probdef(0, 'millimeters', 'planar', 1e-8, 0, 30)
- 单位制选择毫米（'millimeters'）与机械设计保持一致
- 平面问题类型（'planar'）适合轴向长度远大于截面的电机
- 1e-8的精度要求确保齿槽转矩计算准确
- 30°的对称周期可减少计算量
2.2 材料定义与边界条件
永磁体材料参数设置示例：
python复制femm.mi_addmaterial('NdFeB', 1.05, 1.05, 0, 0, 0.5, 0, 0, 1, 0, 0, 0)
- 1.05的相对磁导率对应典型钕铁硼材料
- 0.5T的矫顽力需根据实际磁钢牌号调整
- 运动边界
rotor_bdry的5mm/s速度需与后续谐响应频率对应
2.3 电磁力提取技巧
通过电路属性获取相电流电磁力：
python复制femm.mo_getcircuitvalues('PhaseA')
- 需保存各时间步的齿部径向/切向力密度
- 建议采样率至少为最高分析频率的10倍（本例4kHz对应40kHz采样）
- 力密度分布应导出为CSV格式供结构场使用
3. 结构谐响应分析实战
3.1 电磁力加载方法
ANSYS APDL中采用表数组方式加载时变电磁力：
apdl复制*DIM,force_data,TABLE,6,1000,,TIME
*TREAD,force_data,force.txt
- 6个自由度分量需与电磁仿真坐标系一致
- 1000个时间点需覆盖至少3个电周期
- 阶梯加载（KBC,1）避免引入虚假高频成分
3.2 模态叠加法设置
谐响应分析关键参数：
apdl复制HROPT,MSUP,,YES
HARFRQ,0,4000
DMPRAT,0.03
- 模态叠加法（MSUP）计算效率高于直接法
- 0-4000Hz范围覆盖人耳敏感频段
- 0.03阻尼比适用于多数硅钢片结构
3.3 网格划分策略
不同区域的网格尺寸建议：
| 部件 | 推荐尺寸 | 单元类型 | 说明 |
|---|---|---|---|
| 定子齿部 | 2mm | SOLID185 | 力密度变化剧烈区域 |
| 定子轭部 | 4mm | SOLID185 | 采用扫略网格提高质量 |
| 机壳 | 6mm | SHELL181 | 薄壁结构适用壳单元 |
注意：电磁与结构网格无需一致，但需保证力映射时节点位置对应准确
4. 声学仿真关键技术
4.1 振动速度映射
在LMS Virtual.Lab中进行数据转换：
matlab复制vibData = importdata('vibration.csv');
mesh = ansysReadMesh('model.cdb');
- 需确保振动数据与声学网格时间/空间同步
- 建议采用节点直接映射法保证精度
- 对于大型模型可使用降阶技术加速计算
4.2 边界元法设置
声学求解器关键参数：
matlab复制actool('set', 'MaxFreq', 4000);
actool('assign', 'Velocity', vibData);
actool('solve', 'BEM');
- 边界元法（BEM）适合开放场噪声问题
- 4000Hz上限需与结构分析保持一致
- 建议采用自适应积分提高高频精度
4.3 声压级评价
典型场点声压级计算：
matlab复制fieldPoint = [1200, 500, 800];
spl = actool('getSPL', fieldPoint);
- 驾驶员耳位通常距电机300-500mm
- A计权结果更接近人耳感知特性
- 重点关注800-2500Hz的啸叫频段
5. 工程经验与故障排除
5.1 典型问题解决方案
| 现象 | 可能原因 | 解决措施 |
|---|---|---|
| 共振峰频率偏移 | 网格尺寸过大 | 齿部加密至1.5mm以下 |
| 声压级高频成分异常 | 力密度采样不足 | 提高电磁仿真时间分辨率 |
| 计算不收敛 | 阻尼设置不合理 | 采用瑞利阻尼替代固定阻尼比 |
| 噪声频谱出现谐波失真 | 斜极角度未正确建模 | 电磁场采用3D斜极模型 |
5.2 性能优化技巧
-
并行计算配置：
- FEMM支持多线程求解
- ANSYS使用DSDB模式分布式计算
- 声学仿真采用FMM加速算法
-
模型简化建议：
- 对称周期可缩减至30°模型
- 非承载部件采用质量点替代
- 远场声学使用无限元边界
-
后处理加速方法：
- 电磁力采用FFT频域映射
- 结构响应使用模态参与因子筛选
- 声学结果预存场点数据库
5.3 斜极设计影响
永磁体斜极可显著改变噪声特性：
- 单段斜极：降低齿谐波幅值
- 多段斜极：分散噪声能量谱
- 最优斜极角度通常为1-2个齿距
实际项目中，我通常会先完成直极模型仿真确定主要噪声源，再引入斜极进行优化设计。这种方法能大幅减少计算迭代次数。