IPM电机(内嵌式磁石型)
IPM电机(内嵌式磁石型)的理论基础
概述
老师,我经常听到IPM电机这个名字,但它和SPM到底有什么根本区别?
IPM(Interior Permanent Magnet)电机是一种将永久磁石埋入转子铁芯内部的结构。与表面贴磁石的SPM不同,IPM可以利用铁芯的凸极性。也就是说,d轴和q轴的电感会产生差异。
电感有差异的话有什么好处?
可以使用磁石转矩加上磁阻转矩。这是IPM最大的优势。EV驱动用的主力电机之所以采用IPM,是因为它可以同时实现高转矩密度和宽定功率运行范围。
支配方程
IPM电机的转矩用数式怎么表达?
在dq坐标系中的转矩式如下。
其中 $p$ 是极对数,$\psi_m$ 是磁石磁通链数,$L_d$、$L_q$ 是d轴和q轴电感,$i_d$、$i_q$ 是d轴和q轴电流。
第一项是磁石转矩,第二项是磁阻转矩,对吧?$L_d - L_q$ 越大磁阻转矩就越大?
完全正确。在IPM中,磁石会阻挡d轴磁路,所以 $L_d < L_q$,凸极比 $\xi = L_q / L_d$ 通常在1.5~3之间。凸极比越大,磁阻转矩的贡献就越大。
另外,电压方程也很重要。
从这个电压方程可以看出,高速域的电压限制决定了弱磁控制的必要性。
我明白了。用JMAG进行IPM解析时,就是从FEM中提取这个式子的参数,对吧?
没错。JMAG和Ansys Maxwell都可以按照旋转角计算磁通链数,然后求出dq轴电感。由于存在非线性性,需要注意电流大小会导致电感变化。
电磁界解析的基础方程
FEM实际求解的方程是什么?
采用二维磁标量势 $A_z$ 的扩散方程。
其中 $\nu$ 是磁阻率(透磁率的倒数),$\mathbf{J}_0$ 是外部电流密度,$\mathbf{M}$ 是磁化矢量,$\sigma$ 是电导率。
磁石部分用 $\nabla \times \mathbf{M}$ 表示,涡电流用 $\sigma \partial A / \partial t$ 项表示,对吗?
完全正确。铁芯的非线性B-H特性作为 $\nu(B)$ 被纳入,Newton-Raphson法执行非线性迭代。在IPM中,磁饱和很强,这种非线性处理是解析精度的关键。
实务上的注意点
IPM电机解析需要特别注意的是什么?
让我整理一下重要要点。
- 桥接部分的网格:保持磁石的薄桥(0.5~1mm)会发生严重磁饱和。在这里至少放置3层单元
- 非线性B-H曲线:电磁钢板的种类(35H300、20HIM等)会大大改变特性
- 磁石的温度依存性:NdFeB磁石的剩余磁通密度随温度系数 $\alpha_B \approx -0.12\%/°C$ 下降
- 电流相位角扫描:要找到MTPA(最大转矩/电流)控制点,需要在 $\beta$ 角从0°~90°之间扫描
用JMAG的话,磁石桥接部分可以用自动网格处理吗?
JMAG-Designer有薄板部分的自动识别功能,但不要过度依赖。一定要目视确认网格质量,检查桥接部分的纵横比是否在5以下。
"奖励"磁阻转矩让IPM成为王者
IPM电机战胜表面磁石型(SPM)的最大原因是"磁阻转矩"这一物理现象。除了磁石的转矩外,d轴和q轴电感的差异还产生了这个"奖励转矩",占总转矩的20~30%。丰田在设计普锐斯时发现(或者说巧妙地利用了)这一特性,用同样的磁石数量能产生更大的转矩。"磁石工作,铁也工作"的二刀流就是IPM的本质。
IPM电机(内嵌式磁石型)的数值计算手法
数值解法的详细
IPM电机的解析,具体的计算流程是怎样的?
IPM电机的典型电磁界解析流程如下。
1. 静磁界解析:为各个电流相位角创建dq轴磁通链数地图
2. 瞬态解析:随转子旋转计算转矩波形和感应电压波形
先用静磁界取参数,然后进行瞬态解析,是这样吗?
完全正确。JMAG的"dq轴参数化解析"功能可以自动为 $i_d$-$i_q$ 平面上的每个格点计算磁通链数。这是效率地图制作的基础数据。
dq轴电感的提取
电感是怎样提取的?
从有限元法计算得到的磁通链数 $\psi_d$、$\psi_q$ 计算如下。
这是依赖于电流大小的非线性电感,对吧?
完全正确。由于磁饱和,电感会随着电流大小以非线性的方式变化。特别是q轴方向穿过铁芯,容易饱和,大电流时 $L_q$ 会显著下降。忽视这一点会在弱磁领域的特性预测中产生巨大误差。
转子运动的耦合
旋转的转子和定子之间的气隙怎样处理?
主要有三种方法。
| 方法 | 概述 | 适用工具 |
|---|---|---|
| 滑动网格法 | 按转动重新连接气隙面网格 | JMAG、Maxwell |
| 带网格法 | 在气隙中设置旋转带进行重新网格划分 | JMAG |
| 变形法 | 通过变形气隙单元 | COMSOL |
JMAG的滑动网格,连接面的磁通连续性没问题吗?
JMAG通过插值处理来保持连续性,但气隙中心的周向网格分割数很重要。最少需要每个极距30~50分割。分割过粗的话齿槽转矩的精度会恶化。
转矩计算手法
FEM计算得到的转矩值可信吗?
转矩计算有多种方法。
- Maxwell应力张量法:积分气隙面上的 $B_r \cdot B_\theta / \mu_0$
- 虚功法:从微小旋转的能量变化计算
- Arkkio法:气隙体积内的体积积分
JMAG默认使用Arkkio法,网格依赖性较小。Maxwell主要使用虚功法。
用多种方法计算,结果一致的话信度就高了,是吧?
完全同意。特别是在评估像齿槽转矩这样的微小量时,用2种以上的方法比较是最佳实践。
为什么IPM解析"2D比3D精度高"
IPM电机的FEM解析,90%以上采用2D解析。你可能会想"3D应该更准确吧?",但其实不一定。因为IPM转子是由叠层钢板组成,轴向磁通变化很小,2D模型效果非常好。只要用补正系数近似端部效应(线圈端部的磁通泄漏),2D的计算速度更快,重复计算的精度也更稳定。3D解析留给"想要详细查看斜度和端部效应"时作为最终验证手段。这是现场的常识。
IPM电机(内嵌式磁石型)的实务应用
实践指南
实际设计IPM电机时,应该按什么步骤进行?
实务设计流程如下所示。
第1步:规格确定
- 额定转矩、转速、输出
- 外径、轴长的约束
- 直流电压、最大电流
第2步:初始设计(解析方法)
- 极数、槽数的选择(例:8极48槽)
- 磁石形状(V形、I形、三角形)的选择
- 用MotorCAD解析尺寸模板确定初始尺寸
第3步:FEM详细解析
第4步:优化
- 磁石形状参数的敏感度分析
- 用JMAG的优化功能或ModeFRONTIER耦合进行多目标优化
MotorCAD和JMAG结合使用是标准做法吗?
MotorCAD在热、电磁的简易计算中表现优异,初期检讨很快。JMAG基于FEM可以进行精密的电磁界解析。将两者耦合使用是高效的设计流程。
网格设计的实践
IPM电机特有的网格注意事项请说一下。
按各部位逐一说明。
| 部位 | 推荐网格尺寸 | 原因 |
|---|---|---|
| 气隙 | 0.2~0.5mm(3层以上) | 与转矩精度直接相关 |
| 磁石桥接 | 厚度的1/3以下 | 准确捕捉磁饱和 |
| 槽开口部 | 0.5mm以下 | 漏磁通的精度 |
| 背面轭 | 2~5mm | 降低计算成本 |
气隙和桥接部要特别细致,背面轭可以粗一点,是这样吗?
完全正确。关键是合理管理整体节点数。2D模型在5万~20万单元,3D模型在50万~200万单元是实用范围。
常见失败与对策
初学者容易掉入的坑有哪些?
| 失败 | 原因 | 对策 |
|---|---|---|
| 转矩与理论值相差很大 | 磁石的B-H曲线不正确 | 从厂家目录中准确输入数值。也设置温度依存性 |
| 齿槽转矩的解析值比实测大 | 网格过粗 | 气隙周向分割改为极距50分割以上 |
| 弱磁域效率不吻合 | $L_d$ 的饱和特性不足 | 在宽电流范围内建立电感地图 |
| 铁损被低估 | 未考虑高调波铁损 | 减小时间步,输入PWM波形 |
磁石的B-H曲线可以从厂家获取吗?
可以从信越化学、TDK等磁石厂家获取数据表。JMAG内置了主要磁石等级的材料库,从中选择最可靠。
磁石V形配置"无论试多少次都好用"的设计方法
IPM转子的磁石配置经常看到V形,这不是巧合,而是数十年试错的结晶。将V形角度优化可以使磁阻转矩最大,进一步调整磁通屏障的形状还能抑制高调波成分。某厂商的开发者曾说"V形顶角只改变5°,转矩就会变化几个百分点。没有FEM已经无法判断"。实践中这种"V形+屏障形状"的参数化研究成为了设计的关键。
IPM电机(内嵌式磁石型)的软件比较
商用工具比较
IPM电机的设计、解析可以使用什么工具?
按用途整理一下主要工具。
| 工具 | 开发商 | 强项 | 主要用途 |
|---|---|---|---|
| JMAG-Designer | JSOL | 电机专用模板、材料数据库 | 电磁界FEM解析 |
| Ansys Maxwell | Ansys | Workbench集成、多物理耦合 | 电磁界FEM解析 |
| Motor-CAD | Ansys | 高速初期设计、热/电磁耦合 | 概念设计、热设计 |
| PSIM | Powersim | 电路-FEM协作仿真 | 控制系统设计 |
| COMSOL | COMSOL AB | 自定义方程自由度高 | 研究、特殊解析 |
| SPEED | CDadapco/Siemens | 基于解析方法的高速设计 | 初期检讨 |
JMAG-Designer
JMAG在电机解析中为什么这么强?
JMAG有针对电机设计的专用功能。
- 电机模板:提供了IPM、SPM、异步机、SRM等参数化模型
- 材料数据库:内置约500种电磁钢板、约200种磁石
- dq轴映射:自动在 $i_d$-$i_q$ 平面上生成磁通链数地图
- 效率地图生成:配合损耗计算,自动生成转矩-转速平面的效率地图
有了模板的话初学者也容易着手,是吧?
Ansys Maxwell + Motor-CAD
ANSYS这一侧的组合怎样分工?
典型工作流程是:用Motor-CAD做初期设计和热设计,用Maxwell做精密FEM解析。ANSYS在Workbench平台上,很容易进行结构解析(Mechanical)和热流体解析(Fluent)的耦合。这是ANSYS的强项。
NVH解析也能在Workbench上完成吗?
用Maxwell计算电磁力,转送给Mechanical进行结构模态解析、加载响应解析。从噪声预测一直到全流程贯通,这是ANSYS的优势。
工具选定的指南
最后,应该选哪一个?
来整理一下判断标准。
- 电机专业、成本优先 → JMAG(日语支持也很完善)
- 需要多物理耦合 → Ansys Maxwell + Workbench
- 想与控制系统同时设计 → PSIM + FEM耦合
- 想快速进行热设计 → Motor-CAD
- 研究目的的自定义定式 → COMSOL
汽车厂商的实际情况怎样?
国内JMAG的市场占有率很高。丰田、日产、本田等EV开发部门都广泛使用JMAG。国外也有不少采用Maxwell+Motor-CAD组合的。最近还有一些公司同时使用JMAG和Motor-CAD。
"免费工具+商用工具"的分工是聪明的选择
IPM设计的现场,已经广泛采用"用FEMM(免费)快速验证概念,再用Ansys Maxwell(商用)精密解析"的两阶段方式。FEMM是2D静磁界专用的极简UI工具,学习成本低,几分钟就能建模。商用工具在非线性过渡解析和最优化引擎耦合上很强大,但许可证费用很高。"概念验证用免费,最后解析用商用"的角色分工能大幅改善IPM设计团队的成本效率。
IPM电机(内嵌式磁石型)的先进研究
先进话题
IPM电机领域现在的研究热点是什么?
有几个重要方向。
无重稀土磁石
传统的NdFeB磁石需要添加Dy(镝)以确保耐热性,但稀有性和成本都是问题。最近,不含Dy的磁石等级和使用铁氧体磁石的IPM设计开始引起关注。
铁氧体磁石的剩余磁通密度大约是一半吧。用那样的磁石IPM能成立吗?
通过扩大凸极比,提高磁阻转矩的比例来应对。比如凸极比3以上的设计,用铁氧体磁石也能达到实用级的转矩密度。用JMAG的拓扑优化功能来进行磁路设计的研究正在进展。
拓扑优化
拓扑优化与结构的优化有区别吗?
与参数优化(改变尺寸)不同,拓扑优化是改变材料分布本身。转子内部的铁/空气/磁石配置可以从零开始探索。
这里 $\rho$ 是各单元的材料密度变量(0:空气,1:铁)。
JMAG可以进行拓扑优化吗?
JMAG v22以后搭载了拓扑优化功能。ANSYS通过与OptiSlang耦合也可以实现。不过在实务中,制造性约束(冲裁可能的形状吗)的过滤很重要。
AI代理模型
FEM解析一个案例需要数十分钟,所以探索很多设计点会花很多时间。最近,用学过FEM结果的神经网络(代理模型)快速探索设计空间的方法在扩大。
JMAG的数据能直接传给Python搭建AI模型吗?
JMAG有Python API,可以CSV/NumPy格式输出解析结果。拿这个用scikit-learn或PyTorch学习,就能建立电感地图和效率地图的高速预测模型。JMAG的"AI活用最优化"功能也能利用。
800V系应对
EV的800V化在进行中,对IPM电机有什么影响?
电压上升导致绕组绝缘设计更严格。缠绕间可能出现局部放电(PD: Partial Discharge)的风险。高dv/dt还会导致轴承电蚀。用Ansys Maxwell和Q3D评估浮性容量,用Workbench的结构解析定量评估电蚀风险的方法开始受关注。
多材料磁石突破IPM极限的研究
传统IPM都是"烧结钕磁铁"单一选择,但先进研究中"拓扑优化"能自由设计磁石形状和材料。比如将磁石部分做成铁氧体与钕的混合体,在减少稀有金属使用量的同时可以维持性能。更进一步,AMRF(先进磁石研究设施)还在尝试用3D打印机任意造型磁石,IPM的"磁石必须是长方形"这个常识正在崩溃。
IPM电机(内嵌式磁石型)的故障处理
故障排除
IPM电机解析中常见的故障有哪些?
来整理一下实务中遇到的代表性问题。
1. 转矩在FEM和控制仿真中不匹配
JMAG计算的转矩与PSIM电路仿真结果差20%,这是为什么?
原因和对策:
- FEM的电感地图过粗 → 增加 $i_d$-$i_q$ 格点(最少10×10,推荐20×20)
- 电路模型一侧没考虑铁损 → 将铁损等效电阻并联接到d轴
- 忽视磁饱和的交叉项 → 使用同时包含 $\psi_d(i_d, i_q)$ 和 $\psi_q(i_d, i_q)$ 双变量依赖的地图
2. 齿槽转矩的解析值与实测不符
齿槽转矩的解析值变成了实测的两倍。
可能原因:
- 网格不足:气隙周向分割过粗 → 改为极距最少60分割
- 转动步长过粗 → 采用机械角0.5°以下的步长计算
- 3D效应:2D解析无法考虑斜度效果 → 采用3D或多分片法
反过来实测比解析值小的情况,是因为制造误差(磁石着磁不均、组装偏心)在实机中叠加。
3. 非线性收敛很慢
Newton-Raphson迭代超过100次还不收敛。
对策:
- 扩展B-H曲线的高磁场部分:如果材料数据在2T左右就截断,在饱和区会收敛不好
- 改进初值:用前一个时间步的解作为初值
- 引入缓和系数:在Newton法的更新量上乘以0.5~0.8的系数
- JMAG特有:启用"非线性迭代的自动缓和"选项
4. 3D解析太慢
3D模型每个时间步要30分钟。全转子计算要花很多天。
加速方法:
| 方法 | 效果 | 适用条件 |
|---|---|---|
| 利用周期对称性 | 计算量减少1/极对数 | 结构有对称性 |
| 2D+端部效果补正 | 100倍以上加速 | 轴向磁通可忽略 |
| 多分片法 | 用2D近似3D效果 | 斜度评估 |
| 并行计算(MPI) | 与核心数成正比 | HPC环境 |
先用2D把设计定细,最后只验证几个3D案例,这是比较现实的吗?
完全正确。用2D跑100多个优化案例,最后只有2~3个候选方案用3D验证,这是实务的定律。JMAG和Maxwell都是按这个流程来设计功能的。
5. 减磁判定
在最大电流条件下,磁石是否会减磁,应该怎样检查?
步骤:在最大负d轴电流+最高温度条件下(比如150°C)进行静磁界解析,检查磁石内部的减磁场 $H$ 是否超过保磁力 $H_{cj}(T)$。JMAG有"减磁率地图"显示功能,可以用等高线来可视化各磁石单元的减磁率。
FEM没有问题实机却出现的"磁通屏障裂纹"
IPM转子的薄桥(磁通屏障端)在高速旋转时容易因遗心力裂开。应力解析用FEM(结构)来做,但问题是"同时还有磁吸力"。电磁和结构分别解析的话,磁压力荷载会被遗漏。许多实机故障就是从这个"耦合解析遗漏"产生的。最近电磁-结构耦合解析成为了标准,但仍然有"电磁OK、强度OK,实机却裂了"的情况。
相关话题
价值
更多
错误