简介:这是一份《永磁同步电机Maxwell建模与仿真》的PDF实训资料,面向电力电子与传动控制专业本科生及电机设计初学者。内容围绕一台额定功率15千瓦、电压380伏、4极50赫兹的永磁同步电机,完整覆盖电机参数设置、RMxprt等效模型搭建、Maxwell 2D二维仿真,以及轴长与气隙长度等关键参数的分析优化,帮助读者掌握从参数设定到仿真结果解读的全流程。压缩包仅含1个PDF文件,大小约4.48MB,内含详细的模型图示、参数配置截图和仿真数据(如槽满率73.55%、D轴电抗6.11Ω、Q轴电抗31.60Ω等),并细致展示了RMxprt中定子、转子、绕组、槽形与永磁体材料的设定过程,以及Maxwell 2D的模型与仿真结果,可直接用于课程设计、毕业设计或工程参考。文档章节结构完整,包括参数、建模、仿真和参数分析四大部分,方便对照学习。目前已有1656人学习下载,适合需要快速上手Ansys Maxwell/RMxprt电机仿真,或希望参照实际案例理解永磁同步电机设计要点的读者。
1. 一台 RMxprt 效率 95.3% 的电机,为什么还要做 Maxwell 2D
拿到一份永磁同步电机设计参数,第一反应是直接丢进 RMxprt 把效率、转矩、功率因数算出来。RMxprt 确实快,一次满载求解不到一分钟,结果是效率 95.3016%、功率因数 0.992773、启动转矩 265.252N·m,怎么看都是个不错的四极 15kW 设计。但等效电路法有一个盲区:它算不出转子齿局部磁密到 3.0T 以上、定子轭部 1.62T 这种空间分布信息。Maxwell 2D 的有限元求解恰恰补的是这一块——用 1/4 模型验证气隙磁密波形、转矩脉动和局部饱和点,确认 RMxprt 的结果不是“平均意义上的正确”。这份材料适合电力电子与传动控制方向做电机本体验证、课程设计或毕业设计的同学,也适合手里有 RMxprt 模型但不确定有限元结果怎么对齐的工程师参考。
2. RMxprt 参数化建模:15kW 四极电机从铭牌到等效电路
2.1 先用一张参数表把电机定死
RMxprt 建模的前提是参数完整,缺一个槽形尺寸都会导致求解器报错。把原始参数整理成下面这张表,建模时逐项对照填,不遗漏。
| 参数类别 | 关键数值 |
|---|---|
| 额定功率 / 电压 / 频率 | 15kW / 380V / 50Hz |
| 极对数 / 同步转速 | 2(4 极)/ 1500rpm |
| 定子铁芯 | 外径 260mm,内径 170mm,轴长 190mm,叠压系数 0.95 |
| 定子槽 | 36 槽,单层星形绕组,节距 7,每槽 13 匝,1.16mm 铜线 5 股并绕 |
| 定子槽形(mm) | hs0=0.8,hs1=1.13,hs2=14.07,bs0=3.8,bs1=7.7,bs2=10.2 |
| 转子 | 外径 168.7mm,内径 60mm,内置 W 型永磁体,厚度 5.3mm,总宽 92mm |
| 转子槽形(mm) | hs0=0.8,hs1=1.27,hs2=13.73,bs0=2,bs1=6.4,bs2=5.5,rs=0 |
| 气隙 | 0.65mm |
| 材料 | 磁体 NTP264H(自定义),铁芯 D315_50,导条/端环铸铝 |
| 损耗设置 | 杂散损耗 2%,机械摩擦损耗 160W,风磨损耗 0W |
转速可以顺手验算:n = 60f/p = 60×50/2 = 1500rpm,频率、极对数、转速三者自洽,说明这是一个标准的工频直接并网或变频基频运行设计。
2.2 Machine 与 Stator:先把旋转电机的“壳”搭起来
RMxprt 的 Machine 页面选择 Three-Phase Induction Motor 旁边的 Synchronous Motor 类型时,注意选 Brushless Permanent-Magnet DC Motor 还是 Synchronous Motor 的差异。这份设计是内置式 W 型磁体、三相正弦供电,应该按 Synchronous Motor 处理,调整 Machine 里的定转子叠压系数为 0.95。
接着是 Stator 页,只需要填外径、内径、轴长和槽数。这里有一个容易忽略的点:轴长 lef=190mm 是铁芯净长,不是包含端部的总长。RMxprt 里 Stator Length 填 190mm 后,程序会自动乘以叠压系数 0.95 得到有效铁芯长度,这个值参与磁路计算,等效电路里所有跟磁通截面积相关的量都以它为准。如果填错成物理长度,磁密会整体偏低约 5%,后面 Maxwell 2D 对不齐。
2.3 Slot 与 Winding:槽形尺寸决定磁路截面积
Slot 参数页把定子槽形从圆形切换到自定义多边形时,需要把 hs0、hs1、hs2、bs0、bs1、bs2 全部填进去。这些尺寸不是随便给的:hs0、bs0 决定槽口宽度,直接影响槽漏抗;hs1、hs2 决定槽身深度,影响齿部磁路的宽度;bs1、bs2 决定槽身宽度,影响轭部磁路高度。原始参数里定子槽口 hs0=0.8mm、bs0=3.8mm,属于半闭口梨形槽,是常见的低压中小型电机槽型。
Winding 页面设置单层星形接法时,有四个关键参数:每槽导体数填 13,并联支路数默认 1,线圈节距填 7(跨 1 到 8 槽),并绕根数填 5。这里容易混淆的是“每槽匝数”和“每槽导体数”的换算。RMxprt 的 Winding 输入里每一槽的导体数为 13 匝乘以 5 股等于 65 根导体,但软件界面一般要求填每槽匝数,所以填 13。线径 1.16mm 乘以 5 股后,软件会自动计算等效截面积和槽满率。仿真结果里槽满率 73.5471%,这个数值偏高但不超限,工程上手工嵌线一般控制在 75% 以内,说明绕组设计接近工艺上限。
2.4 材料自定义:NTP264H 磁体的 B-H 属性
材料库没有 NTP264H 时,必须手动新建。打开材料管理器,Materials 列表中右键 New Material,先设相对磁导率为 1.05,再输入剩磁 Br 和矫顽力 Hc。原始资料里没有给出这个牌号的完整 B-H 数据,常见做法是按 Maxell 材料库中同类型钕铁硼磁体的参数推算:剩磁 Br 取 1.1~1.2T,矫顽力 Hc 取 800~900kA/m 之间,然后勾选 Nonlinear B-H Curve 输入去磁曲线第二象限的采样点。
有一个细节值得注意:RMxprt 材料页的永磁体默认方向是径向或平行充磁,但内置 W 型磁体是切向充磁为主。需要在 Pole 参数页的 Magnet 设置里把充磁方向改成 Alternate 或 Custom,软件会根据磁体形状自动计算磁化矢量。如果方向不对,空载气隙磁密会显著偏低——原始仿真结果气隙磁密 0.667467T,这个值对于内置式电机是合理的,说明磁化方向设置基本正确。
2.5 仿真结果:从电抗到磁密,哪些数据要重点看
RMxprt 求解后,Data 表里会列出稳态参数和满载数据,原始结果中几个关键值需要按工程经验校验:
| 仿真输出 | 数值 | 工程判断 |
|---|---|---|
| 槽满率 | 73.55% | 偏高但在手工嵌线可接受范围内 |
| D 轴电抗 | 6.11407Ω | 偏小,说明直轴磁路磁阻大 |
| Q 轴电抗 | 31.5992Ω | 远大于 D 轴电抗,典型的凸极率约 5.17 |
| 定子轭磁密 | 1.62085T | 接近 D315_50 硅钢片饱和拐点,需关注 |
| 额定转矩 | 95.4879N·m | 与 15kW@1500rpm 的理论值 95.5N·m 吻合 |
| 功率因数 | 0.992773 | 偏高,实际考虑谐波会略低 |
| 效率 | 95.3016% | 合理,杂散损耗 2% 已被计入 |
Q 轴电抗是 D 轴电抗的 5 倍以上,这是 W 型磁体布局的典型特征——直轴磁路经过磁体,磁阻大电抗小;交轴磁路大部分走铁芯,磁阻小电抗大。转矩角 68.04° 说明功角较大,电机工作在最大转矩附近,单位电流输出转矩的能力强。这些数据进入 Maxwell 2D 后会直接成为对照基准。
3. Maxwell 2D 有限元验证:边界条件、激励换算与转矩误差标定
3.1 为什么取 1/4 模型
Maxwell 2D 瞬态场计算的最大问题是时间和内存。完整 36 槽 4 极模型的网格量在默认设置下通常超过 30 万单元,而 1/4 模型只保留 9 槽 1 极,网格量减少到 8 万左右,单次求解时间从 40 分钟降到了 8 分钟以内,波形精度几乎不受影响。
取 1/4 模型的依据是电机结构的对称性:4 极电机在空间上每 90° 重复一次,定子 36 槽均匀分布,转子 W 型磁体在四个磁极下完全一致。只要三相绕组激励保持对称,磁场分布就天然满足半个极距的周期性。前提是绕组排布必须严格对称,原始参数里节距 7、每槽 13 匝、单层绕组满足这一条件。
3.2 边界条件:Master、Slave 与 Vector Potential 的配合
1/4 模型的边界处理方式,原始资料写得很清楚:曲线外边界用 Vector Potential(狄里克莱边界,A=0),横轴用 Master,纵轴用 Slave。三者的物理含义需要区分:
- Vector Potential 边界把模型外边界上的磁矢位固定在 0,等价于认为电机外部无漏磁,磁力线平行于边界切线方向。对于电机这种铁芯包在内部的几何体,外边界距离机壳足够远时这个近似完全成立。
- Master 和 Slave 是配合使用的周期边界对,Master 边界的磁场矢量会映射到 Slave 边界上,保证两个边界上的磁位连续。旋转机械里 Master/Slave 还分旋转周期性和平移周期性,这里 1/4 模型是旋转对称,需要在边界设置里把 Slave 的坐标旋转角度填 90°。
- 如果边界条件设反了(比如把对称边也设成 Vector Potential),求解不会报错,但磁力线会被强制垂直切开,转矩波形会出现周期性畸变,转矩平均值明显偏低。判断边界设对没设对的方法很简单:看磁密云图在 Master/Slave 交界处是否连续。
3.3 激励源设置:从线电压 380V 到 310.269V 的换算
Excitation 里接外部电路或者直接加电压源时,A 相激励表达式为 310.269sin(2pi50time+delta),B 相滞后 120°,C 相超前 120°。310.269 这个值不是拍脑袋定的,来源是线电压 380V 的有效值换算成相电压峰值:
import math # 线电压有效值 380V,星形接法 U_line_rms = 380.0 U_phase_rms = U_line_rms / math.sqrt(3) # 相电压有效值约 219.39V U_phase_peak = U_phase_rms * math.sqrt(2) # 相电压峰值约 310.27V print(f"相电压有效值: {U_phase_rms:.3f} V") print(f"相电压峰值: {U_phase_peak:.3f} V") print(f"激励幅值取整: 310.269 V")这段代码的逻辑是:星形接法下相电压等于线电压除以根号 3,得到有效值 219.39V,再乘以根号 2 得峰值 310.269V。delta 是转矩角,RMxprt 算出的转矩角 68.04° 换算成弧度后填入 delta。需要注意:derta 设置为 0 时电机输出转矩很小,必须填对;如果不清楚具体角度,可以用 Maxwell 的 Optimetrics 扫描 delta 从 0° 到 90°,找到最大平均转矩对应的角度当作初始工作点,但这会额外增加几倍求解时间。
原始报告里还提到电流源和电压源激励的对比:电压源激励下转矩平均值 260.0075N·m,电流源激励下波形更稳定但平均转矩略有不同。实际使用中,电压源更接近电机直接并网的真实工况,而电流源适合做控制策略仿真。如果只关心电磁设计验证,优先用电压源,因为电流源激励的幅值和相位需要先知道电流角,这个值在电磁方案未定型前是未知的。
3.4 转矩对比:RMxprt 与 Maxwell 2D 的误差标定
仿真输出的转矩数据是这组项目最重要的验证环节。RMxprt 启动工况转矩 265.252N·m,Maxwell 2D 电压源激励下转矩波形平均值为 260.0075N·m,误差只有 1.98%。这个误差水平说明 RMxprt 等效电路中的参数(尤其是饱和后的电感)和有限元耦合得比较好。
| 项目 | RMxprt | Maxwell 2D | 误差 | 说明 |
|---|---|---|---|---|
| 启动/额定转矩(N·m) | 265.252 | 260.0075 | 1.98% | 工程上 5% 以内可接受 |
| 气隙磁密(T) | 0.667467 | 0.6946 | 4.07% | 有限元考虑局部饱和,数值略高 |
| 定子齿磁密(T) | 1.40687 | 1.5 附近 | 约 6.6% | 齿部饱和区有限元更敏感 |
有两处需要注意:Maxwell 2D 求解的是 1/4 模型,转矩读取时要确认软件已经把模型周期数考虑进去——大多数版本会自动按模型比例放大转矩,但自定义边界时可能关掉这个设置,导致转矩变成实际的 1/4, 260N·m 读成 65N·m 这种低级错误不少见。另外一个问题是电磁转矩波形里包含明显的启动振荡,平均转矩的取法应该跳过前面几个周期,等波形稳定后再取平均值。
3.5 磁密分布与气隙磁密波形:有限元多算出的信息量
磁密分布图给出了 RMxprt 给不出的空间信息。定子齿和轭部磁密在 1.5T 左右,与 RMxprt 的 1.40687T 和 1.62085T 基本一致;但转子齿局部达到 3.0T 以上,远高于转子轭部平均值 1.07636T。这说明 W 型磁体之间的隔磁桥区域磁饱和极严重,磁力线被挤压在很窄的转子齿桥上。D315_50 硅钢片在 1.6T 以上就进入深度饱和区,3.0T 意味着该区域的相对磁导率已经降到接近空气的水平。
局部过饱和的直接代价是两个:一是该区域铁耗密度极高,热仿真时会出现局部热点;二是磁桥漏磁增加,实际气隙磁密会比理想磁路计算值低。对应做法是调整 Rib 尺寸(原始数据 8mm)和磁体间距 O1/O2,把磁桥加宽 1~2mm 或者在磁桥根部开辅助槽,让磁力线有更多的流通路径。
气隙磁密波形验证是这章最后一步。Maxwell 2D 提取气隙圆周上的径向磁密分量,得到接近正弦的波形,平均值 0.6946T 与 RMxprt 的 0.667467T 误差 4.07%。工程上这个偏差主要来自饱和效应的建模差异,有限元结果更可信。波形越接近正弦,反电动势谐波越小,转矩脉动越低——这也是后续参数扫描里要持续跟踪的指标。
4. 轴长与气隙的参数扫描:效率、转矩与饱和度的取舍
4.1 参数化扫描的正确打开方式
RMxprt 自带参数分析功能,不需要手动改完参数重新求解再记录结果。常见做法是在 Setup 里勾选 Parametric,把轴长 lef 或气隙长度定义为变量,设置起始值、终止值和步长,软件会自动跑完所有组合并输出效率、转矩、磁密等指标的曲线。
推荐的扫描设置如下:
| 参数 | 起始 | 终止 | 步长 | 关注指标 |
|---|---|---|---|---|
| 轴长 lef | 170mm | 210mm | 5mm | 效率、额定转矩、定子轭磁密 |
| 气隙长度 | 0.5mm | 0.8mm | 0.05mm | 气隙磁密、功率因数、启动转矩 |
扫描前要把之前在 2.3 节里填错的参数逐项确认,尤其是轴长,因为扫描结果里每一条曲线都会受到基础参数错误的影响。原始报告里轴长固定在 190mm,扫描以它为中心向两边扩展。
4.2 轴长 lef 的敏感性分析
轴长增加时,磁路截面积线性增大,磁密下降,定子铁耗减小但铜耗不变(电流不变),效率曲线呈先升后平的趋势。在 190mm 附近每增加 5mm,定子齿磁密大约下降 3%~4%,效率提升约 0.2 个百分点;但超过 200mm 后转矩提升趋于饱和,因为漏抗和端部漏磁也在同步增加。
轴长从 170mm 增加到 210mm,额定转矩从约 88N·m 升到 104N·m,接近线性。设计时如果空间允许,加大轴长是提升过载能力最直接的手段。但轴长还受机械振动和临界转速制约——轴越长,转子一阶弯曲临界转速越低,对于 1500rpm 的电机,190mm 已经是比较稳妥的上限,再长需要校核转轴挠度。
4.3 气隙长度对性能的双向影响
气隙从 0.5mm 加到 0.8mm,每增大 0.05mm,气隙磁密大约降低 1.5%~2%。原因很简单:气隙磁路的磁阻与气隙长度成正比,永磁体提供的磁动势不变,气隙长了磁压降就大,磁密自然掉。
但气隙并不是越小越好。小气隙带来高磁密的同时,齿槽转矩幅值和转矩脉动会明显增大,这是有限元扫描时最容易看出的趋势。0.5mm 气隙的齿槽转矩峰值可能比 0.65mm 大 50% 以上。另一个工程问题是装配公差:定子内径 170mm、转子外径 168.7mm,气隙 0.65mm,单边公差如果控制在 ±0.03mm,装配误差已经占气隙的 10%,气隙再小到 0.5mm 时,偏心导致的不平衡磁拉力会迅速上升,对轴承寿命不利。
原始参数选 0.65mm 是在磁密水平、齿槽转矩、工艺可行性之间取了一个折中。扫描结果如果显示 0.60mm 时效率提升不到 0.1% 但齿槽转矩增加 20%,就应该保留 0.65mm 不动。
4.4 扫描结果怎么落回设计
参数扫描结束后要做的不是选效率最高的点,而是找效率、转矩、磁密、工艺裕量四个维度的交集。比如轴长 210mm 效率最高,但转子轴承跨距变大,振动风险增加;气隙 0.55mm 磁密和转矩最好看,但装配良率下降。原始设计维持轴长 190mm、气隙 0.65mm 是合理的,扫描的价值在于确认这个组合不在性能悬崖边上——即相邻参数组合的性能变化率低于 1%,此时设计是鲁棒的。
扫描时刻关注的是转子齿最大磁密。原始数据已经出现 3.0T 的局部饱和,如果轴长缩短或气隙减小,这个值会冲得更高。当转子齿最大磁密超过 2.2T 时(即超过硅钢片饱和拐点约 40%),该方案就应该直接标记为不可用,优先通过增大 Rib 或调整磁桥结构来缓解,而不是继续优化其他参数。
5. 后处理技巧:用场计算器定位转子齿磁密饱和点
Maxwell 2D 求解完成后的第一件事不是读转矩,而是打开 Field Overlays 看磁密云图。云图上出现大片红色区域并不等于过饱和,必须先确认红色的绝对值是多少。默认色标是自动缩放的,最高值可能被某个极小区域的 3.2T 拉到上限,导致 1.5T 的区域显示成橙色——这会造成“处处饱和”的错觉。处理办法是手动把 Color Map 的上限固定到 2.0T,上限之外的区域用灰色覆盖,这样云图里超过 2.0T 的区域一眼就能锚定。
定位具体饱和点用场计算器而不是肉眼。在 Maxwell 2D 后处理里选 Geometry 菜单,创建一条穿过转子齿中线的 Polyline,然后在 Field Calculator 里构造表达式提取该路径上的 B 值。具体操作是用 Calculator 的 Named Expression 功能新建一个表达式:
Saturation_Margin = Mag(B) - 2.0[T]这个表达式的含义是磁密幅值减去 2.0T 参考线,得到的正数就是超过饱和阈值的裕量。把 Saturation_Margin 画到刚才的 Polyline 路径上,输出 B 随距离的曲线,曲线过零点的位置就是进入深度饱和的边界。如果整个转子齿路径上的 Saturation_Margin 都大于 0.3T,说明该齿全截面过饱和,仅靠参数微调已经不够,需要改磁桥拓扑。
对于 W 型磁体结构,重点扫三条路径:两个磁体之间的中间齿桥、磁体端部与转子外圆之间的边桥、以及转子内圆下方的连接桥。三个位置的磁密峰值通常不同,中间齿桥最容易过饱和,因为两侧磁体的磁通都汇聚到这里。用场计算器分别提取三条路径的峰值磁密,填入 Excel 后做一个小表格对比,就能量化哪条路径是漏磁的主要通道。
还有一个小技巧:将场计算器保存的表达式拖入项目树的 Output Variables 里,后续每次重新求解(比如扫描气隙长度)后,软件会自动按同一套表达式重新计算饱和裕量。不同设计方案的饱和情况就可以直接拉数据对比,不需要手工重复描点。这个习惯比每次重新打开云图人工读值要省时间得多。
本文还有配套的精品资源,点击获取