从RMxprt到Maxwell：交流绕组感应电动势仿真与绕组系数验证

简介：面向电机相关专业学生与电气工程技术人员，这是一份以ANSYS Maxwell为工具讲解交流绕组感应电动势和磁势计算的上机实践文档。文档从创建项目开始，逐步演示基于RMxprt模块建立同步电机模型，覆盖Machine、Stator、Rotor、Shaft等关键参数设置，并给出励磁电流7A、三相电流0A等典型激励条件，再通过求解与后处理观察磁力线、磁密云图及感应电动势随转速的变化；同时对比集中绕组与双层短距绕组的建模差异，有助于深入理解绕组结构与电磁性能的关系。资源为单一docx文件，压缩包大小约101.54MB，内容以图文步骤和操作说明为主，适合边看边操作。目前已有387人学习下载，可作为电机学课程上机练习、仿真实验或自学ANSYS Maxwell的参考资料。

1. 从 RMxprt 到 Maxwell：交流绕组感应电动势仿真先把绕组系数定住

电机学教材里，交流绕组感应电动势只有一行公式：E = 4.44fNΦk_w，k_w 是绕组系数，等于分布系数与节距系数的乘积。理论课讲清这个公式三分钟就够；可要在 ANSYS Maxwell 里把它算出来，得依次过两套软件：先在 RMxprt 里把同步电机的 Machine、Stator、Rotor、Shaft 四组参数填好，生成 Maxwell 模型后设置励磁电流 7A、三相电流 0A 的空载工况，再从磁力线、磁密云图和反电动势波形里把理论值验证回来。这套流程对两类人尤其合适：做电机课程设计或毕设的学生，需要把绕组系数落到具体模型上；电机厂的仿真工程师，要快速比较不同绕组方案的反电动势谐波。下面按实际操作顺序讲，参数背后的物理量比按钮位置重要。

2. RMxprt 同步电机参数设定：极槽配合、绕组节距与励磁绕组

新建项目并保存后，在 ANSYS Maxwell 里插入 RMxprt Design，进入电机模板选择窗口。文档里用的是电励磁同步电机模板，不同版本里名字略有差异，常见的是 Synchronous Machine Electrical Excitation 一类。这一章的四个子项必须按 Machine、Stator、Rotor、Shaft 的顺序填，RMxprt 的校验机制不允许跳项，任何一项不完整，后面的 Analyze 都不会通过。

2.1 Machine 项：极对数、额定转速与每极每相槽数

Machine 项里先定额定功率、额定电压、额定频率、极对数和额定转速。极对数 p 直接决定空载电动势频率：f = p·n/60，额定转速和频率必须自洽，否则后面生成 Maxwell 2D 模型时，转速和磁场频率对不上，反电动势波形的周期会和你预期差得很远。

极槽配合是另一个关键决策。设槽数为 Z、相数 m（这里取 3），每极每相槽数 q = Z/(2pm)。q 值决定分布系数 k_d = sin(qα/2)/(q·sin(α/2))，其中 α = 360°·p/Z 是槽距电角度。q 越大，绕组分布在越多槽里，电动势叠加后越接近正弦，但槽数多意味着加工成本上升。文档里的同步电机按整数槽绕组处理，q 取 2 到 4 都常见；如果槽数和极数组合导致 q 不是整数，RMxprt 会判定为分数槽并报警。分数槽绕组能做，但谐波特性和整数槽完全不同，第一次上机不建议碰。

2.2 Stator 项：槽形、线径与绕组系数

Stator 项填定子外径、内径、铁芯长度、槽形尺寸和材料。内径与转子外径共同决定气隙长度，气隙越小主磁路磁阻越小，同样励磁电流下气隙磁密越高，反电动势越大；但气隙太小会让齿槽转矩和谐波更突出。槽形参数经常被忽略，梨形槽需要填槽口宽、槽口高、槽深、槽底宽一组值。齿宽太窄，磁密容易超过硅钢片饱和点，反电动势波形削顶；槽面积太小，槽满率超过 75% 时 RMxprt 直接报警。材料选电工钢片库里的 DW 系列，仿真用的 B-H 曲线决定磁路是否饱和，不要为省事选线性材料。

绕组参数这一组直接对应绕组系数，核心项如下：绕组层数决定能否做双层短距，单层绕组两个线圈边跨整距，双层绕组上下层错开才能实现短距；线圈节距以槽数为单位，填小于极距的值就是短距；并联支路数影响每相串联匝数；每槽导体数决定匝数 N；线径影响电阻和槽满率。文档开头标注的励磁绕组线径选 1.25mm，在 7A 电流下折算电流密度约 5.7A/mm²，处于风冷电机常规范围。要特别说明：电流源激励的电磁仿真里，线径只影响铜损和槽满率校验，不改变反电动势波形，这一点新手常误解。

绕组系数不用手工按计算器，参数确认后可以用几行 Python 先算一遍，心里有数再进 Maxwell：

import math p = 2 # 极对数 Z = 48 # 定子槽数 m = 3 # 相数 tau = Z / (2 * p) # 极距, 槽数 = 12 q = Z / (2 * p * m) # 每极每相槽数 = 4 alpha = 360.0 * p / Z # 槽距电角度 = 15 度 kd = math.sin(q * math.radians(alpha) / 2) / (q * math.sin(math.radians(alpha) / 2)) beta = 10.0 / tau # 节距比, 10 槽节距 / 12 槽极距 = 5/6 kp1 = math.sin(beta * math.pi / 2) # 基波节距系数 kw1 = kd * kp1 # 基波绕组系数 print("q =", q, "kd =", round(kd, 4), "kp1 =", round(kp1, 4), "kw1 =", round(kw1, 4))

这段脚本里，alpha 换算成电角度后参与分布系数计算；beta 用节距除以极距得到，基波节距系数 kp1 = sin(βπ/2)。这里的 q、beta 都是无单位比值，脚本算出的 kw1 可以直接和后面仿真反电动势反推的绕组系数对账。如果 q 和 beta 填错，这里算出来的值和仿真结果会对不上，这是定位参数错误最快的方法。

2.3 Rotor 项与 Shaft 项：励磁绕组和转轴别漏填

Rotor 项填转子外径、凸极或隐极结构、极弧系数、励磁绕组匝数和线径。极弧系数决定气隙磁密沿圆周的分布形状，对反电动势波形的高次谐波有直接影响。励磁绕组数据最容易漏——RMxprt 需要匝数和线径填全才能算出励磁磁动势，生成 Maxwell 模型时才会出现励磁绕组端子，后面才能设置 7A 励磁电流。Shaft 项只填转轴直径和材料，不参与电磁计算，但它决定转子转动惯量，做瞬态启动或转速阶跃仿真时影响机械时间常数；只算稳态空载反电动势时，取结构钢默认值即可。

2.4 模型检查：RMxprt 在求解前拦下哪些错误

参数填完点 Analyze，RMxprt 先做内部校验。常见被拦下的问题：极对数和槽数组合导致每极每相槽数异常、励磁绕组没有闭合、槽满率超限、线径和电流密度组合不合理。遇到提示不要直接改数字硬凑，先回对应子项看单位——RMxprt 默认长度单位是毫米，但不同版本的下拉框可能切到英寸，定子内径差一个数量级是这类上机最常见的错误来源。校验通过后，右键 Design Settings 生成 Maxwell 2D 模型，进入激励和求解阶段。

3. 空载激励与瞬态求解：励磁电流 7A 下时间步长和收敛判据

3.1 激励源设置：为什么空载要把三相电流设为 0A

生成 Maxwell 模型后，设计树里会出现励磁绕组端子和 A、B、C 三相绕组端子。激励设置里最关键的理解点是“空载”：励磁电流 7A 建立主磁场，电枢三相电流取 0A，相当于把电枢绕组开路，只测励磁在电枢绕组里感应出的电动势，此时端电压就等于反电动势，这是后面所有对比的基准。很多人把空载误解为删掉电枢绕组，这是错的——绕组必须存在才能感应电压，电流为 0 和绕组不存在是两回事。

图形界面里，选中励磁绕组端子，右键 Assign Excitation → Current，填 7A；三个相绕组同样 Assign Excitation → Current，填 0A。用脚本批量设置时，常见写法如下：

import ScriptEnv ScriptEnv.Initialize("Ansoft.ElectronicsDesktop") oDesktop.RestoreWindow() oProject = oDesktop.GetActiveProject() oDesign = oProject.GetDesign("Maxwell2DDesign1") oModule = oDesign.GetModule("BoundarySetup") # 励磁绕组: 电流源 7A, 施加在线圈端子上 oModule.EditExcitation( "FieldWinding", "Type:=", "Current", "Current:=", "7A", "IsSolid:=", False) # 电枢三相: 空载开路, 电流 0A for ph in ["PhaseA", "PhaseB", "PhaseC"]: oModule.EditExcitation( ph, "Type:=", "Current", "Current:=", "0A", "IsSolid:=", False)

脚本里 EditExcitation 第一个参数是端子名，要和设计树里的实际命名一致，不一定是 PhaseA；Type 指定电流源激励；IsSolid 为 False 表示激励加在绕组端子上而不是实心导体截面。脚本的好处是反复修改参数做对比时不会漏改某个绕组，也比逐项点鼠标快。

注意：电流源激励下，绕组电阻和线径不会改变反电动势波形；真正影响它的是气隙长度、极弧系数、绕组分布和节距。

3.2 求解选项：停止时间、时间步长和非线性残差怎么配

求解类型必须选瞬态场 Transient。静磁场只能算磁力线和磁密，算不出随时间变化的电动势；涡流场是频域求解，给定频率算稳态，也拿不到完整波形。只有瞬态场能把磁链随转子位置变化的过程逐时间步解出来。进入 Analysis Setup 后，核心参数是停止时间、时间步长、非线性残差三个。

停止时间至少覆盖一两个电周期。电周期 T = 60/(p·n)，极对数 p=2、转速 1500rpm 时 f=50Hz、T=20ms，停止时间取 40ms 够看两个周期。时间步长决定谐波分析能分辨到多少次，采样频率 fs = 1/Δt，要看 20 次以内谐波，步长取 T/200 到 T/500，也就是 0.1ms 到 0.04ms。步长太大，反电动势波形带锯齿，FFT 出来的高频成分是假的。非线性残差控制每步迭代里磁导率更新的收敛精度，默认 1e-6 对硅钢片电机足够，改到 1e-7 只会让计算时间翻倍。求解设置的常见脚本：

oModule = oDesign.GetModule("AnalysisSetup") oModule.InsertSetup("Transient", ["NAME:Setup1", "Enabled:=", True, "StopTime:=", "40ms", # 按电周期换算, 至少两个周期 "TimeStep:=", "0.1ms", # 小于电周期的 1/200 "SaveFieldsType:=", "Every N Steps", "N Steps:=", "1", "NonlinearResidual:=", 1e-6, # 磁导率迭代收敛阈值 "NonlinearMaxIter:=", 10]) # 单步最大迭代次数

这里 StopTime 和 TimeStep 按 3.2 的表对应转速换算，NonlinearMaxIter 如果频繁触发，先回头查网格而不是调大这个数。网格方面，气隙处至少要三层单元，否则气隙磁导算不准，反电动势波形会带空间谐波毛刺。自动网格在常规极槽配合下够用，但手动细化一次气隙网格，是后处理波形干净与否的分水岭。

不同转速下的时间参数换算如下：

3.3 模型检查与求解：报错信息怎么定位

设置完点 Validate，Maxwell 会检查几何、边界、激励、运动设置四类内容。最容易出现的两类报错：一是运动 Band 没有完整包住转子，报运动体相关错误，回 3D Modeler 检查 Band 半径和长度；二是绕组端子没有构成闭合线圈，报激励相关错误，说明绕组正负端子没配对。改成电流源 0A 的空载工况后如果求解器报不收敛，通常不是激励问题，而是网格太粗或铁芯 B-H 曲线数据不完整，先查材料再查网格。正常求解时每时间步 4 到 6 次迭代收敛，如果某一步迭代到上限还没收敛，波形会在该时刻出现孤立尖峰，后处理时要能识别这是数值问题而不是电机本身的特性。

4. 后处理判读：磁力线、磁密云图与反电动势随转速的变化

4.1 磁力线：极对数、对称性与漏磁的一次性检查

后处理调出 Fields 面板，选磁矢位 A 的等值线图看磁力线，重点检查三件事：极对数是否符合预期、磁路是否对称、漏磁是否明显。p 对极电机应该看到 2p 个磁力线汇聚中心，空载下相邻磁极极性交替。如果磁力线分布左右不对称，多半是运动 Band 或边界设置不对称造成的假象。漏磁主要看槽口和磁极边缘的小闭环磁力线，它不参与绕组交链，对反电动势大小无贡献，但影响绕组漏感，做短路工况时才显性化。用 Animate 功能让磁力线随时间刷新，可以看到磁极随转子转动、磁场旋转的过程，和理论课的旋转磁场概念直接对应。

4.2 磁密云图：饱和程度决定反电动势波形会不会削顶

磁密云图看 B 幅值分布，重点看齿部和轭部颜色。冷轧硅钢片在 1.5T 到 1.7T 以下线性度良好，超过 1.8T 后磁导率快速下降进入饱和区。齿部磁密大面积超过 1.8T，气隙磁密波形会被削顶，反电动势出现平顶、谐波增大，这是磁路设计问题，回 RMxprt 里加齿宽或加铁芯长度缓解，不是绕组问题。比云图更能说明问题的是气隙磁密沿圆周的径向分量曲线：在气隙中间画一条圆弧，把 B 的径向分量画成 1D 曲线。磁密波形接近正弦，反电动势才接近正弦；磁密有平顶或马鞍形，反电动势就有对应次数的谐波，这个因果链是后处理里最值得养成的判断习惯。

4.3 感应电动势波形：导出数据、FFT 与理论值对账

建 Transient Report，纵轴取 A 相绕组感应电压。空载下 A 相端电压就是反电动势，理论值满足 E = 4.44fNΦk_w。为了定量核对，把波形导出成 CSV，用 Python 做 FFT 看各次谐波幅值，顺便算 THD，常见做法如下：

import numpy as np data = np.loadtxt("emf_a.csv", delimiter=",", skiprows=1) t = data[:, 0] * 1e3 # 时间, ms emf = data[:, 1] # A 相电压, V dt = (t[1] - t[0]) * 1e-3 # 采样间隔, s fs = 1.0 / dt N = len(emf) win = np.hanning(N) # 抑制频谱泄漏 spec = np.fft.rfft(emf * win) # 加窗后做 FFT freq = np.fft.rfftfreq(N, dt) amp = 2.0 * np.abs(spec) / np.sum(win) # 幅度归一化 f1 = 50.0 # p=2, n=1500rpm 时的基波频率 for k in [1, 5, 7, 11, 13]: idx = np.argmin(np.abs(freq - k * f1)) print(f"{k:2d} 次: {amp[idx]:6.3f} V")

这段代码里，hanning 窗用来处理非整周期截断的频谱泄漏；归一化除以窗函数求和而不是 N，是因为加窗后信号能量被加权，直接用 2/N 会把幅值算小。基波频率 f1 由极对数和转速决定，改成其他转速要同步修改。FFT 结果里基波幅值应接近理论值 4.44fNΦk_w，5、7 次谐波越小说明绕组和磁路设计越干净。如果报告里取的是线电压而不是相电压，基波幅值会大 √3 倍，对账前先统一这个口径。

提示：FFT 前先确认波形已进入稳态，只取最后两个周期的数据段，否则启动暂态段会往频谱里注入虚假谐波。

4.4 转速改变后幅值和频率怎么变

改变转速是文档里单列的一条 Tips，背后是两个公式：频率 f = p·n/60，反电动势 E ∝ fΦ ∝ nΦ。励磁电流不变、磁路不饱和时 Φ 基本不变，转速升高到两倍，频率和电动势幅值同时翻倍。在 Maxwell 里改 Motion Setup 的机械转速（注意是机械转速 rpm，不是角速度），改完重新求解，记录转速、频率、基波幅值三组数据验证线性关系：

如果测出来不是线性关系，优先怀疑两条：磁路饱和——转速高、频率高，铁耗影响变大，波形幅值偏低；或者时间步长没跟着电周期调整，高频下采样点不足，FFT 结果偏小。把步长按 3.2 的规则重设后一般就对了。

5. 集中绕组与双层短距绕组：用节距系数解释反电动势谐波差异

5.1 集中绕组：q=1 时分布系数等于 1 意味着什么

集中绕组指每极每相槽数 q=1 的绕组，线圈绕在单个齿上，端部短、下线简单。q=1 时分布系数 k_d = 1，绕组边在空间上没有分散，反电动势波形直接继承气隙磁密形状，5、7 次谐波只能靠磁路去压，比如调极弧系数或降低磁密。集中绕组的线圈跨距小，节距系数也不等于 1，谐波谱和整数槽分布绕组明显不同。拿它当对比基准，可以把分布效应和短距效应分开看：先看 q=1 的谐波底子，再看双层短距能把哪几次谐波压下去。

5.2 双层短距绕组：节距系数把 5/7 次谐波压到四分之一

双层短距绕组把每个槽分成上下两层，上层线圈边属于一个线圈，下层属于相邻线圈，线圈节距可以比极距短几个槽，得到 β < 1。第 ν 次谐波的节距系数为 k_pν = sin(νβπ/2)，以 β=5/6 代入：

基波只降到 0.966，损失 3.4%，而 5 次、7 次被压到 0.259，衰减近四分之三——这两个次数正是三相电机反电动势里最主要的低次谐波，所以短距用很小代价换来波形质量提升。气隙磁势谐波同理，短距对电枢磁势的削弱和它对空载电动势的削弱来自同一套绕组系数，这就是文档标题把“电动势和磁势”并列的原因。

RMxprt 里做短距，绕组层数必须选 Double Layer，线圈节距填小于极距的槽数。以 48 槽 4 极电机为例，极距 τ=12 槽，整距节距填 12，β=5/6 对应节距填 10。验证时复制一份模型，只把绕组层数改成 Double Layer、节距填 10，其余几何、材料、励磁都不动，重新生成 Maxwell 模型求解，然后对比两个模型的 A 相反电动势 FFT 谱。预期结果是双层短距的基波略低、5/7 次谐波显著下降、波形 THD 明显变好。如果 5 次谐波没有下降，先回 RMxprt 检查层数是否选成了 Single Layer、节距是否填的是极距——这是这个实验里最常见的两个原因。

本文还有配套的精品资源，点击获取