开关磁阻电动机（SRM）的电磁场仿真
开关磁阻电动机（SRM）电磁场理论基础
SRM的工作原理
开关磁阻电动机不使用永磁体吗？不需要稀土？
完全正确。转子只是简单的铁制凸极结构，没有磁石也没有线圈。稀土无关，成本便宜、高温耐性强是优点。不过转矩脉动和噪声较大，这是需要解决的问题。现在Nidec的洗衣机用SRM通过AI控制已实现噪声降低60%的案例。
没有磁石怎样产生转矩呢？完全想象不出…
简单来说，利用"铁被磁石吸引的力"。在定子线圈中通电产生电磁体，转子铁制凸极就被吸向那里。这就是磁阻转矩。磁路电感最大的位置，转子倾向于运动到那里。
那么按时序切换线圈电流，就能连续旋转了？
正是这样。典型的8/6结构（定子8极·转子6极）中，依次通电A相→B相→C相→D相，转子凸极逐个被定子凸极吸引而旋转。但相切换的时序和电流波形至关重要，弄错会导致没有转矩甚至反转。
SRM基本结构特征总结如下。
| 项目 | SRM | PMSM（参考） |
|---|---|---|
| 转子结构 | 凸极铁心仅有（无磁石·线圈） | 永磁体嵌入/贴表面 |
| 稀土使用 | 无 | 钕磁体使用 |
| 耐热性 | 250°C以上也能工作 | 受磁石消磁温度限制（约150°C） |
| 转矩脉动 | 大（20〜40%） | 小（3〜8%） |
| 电磁噪声 | 大（问题） | 小 |
| 成本 | 低 | 磁体成本高 |
| 主要应用 | 家电、工业、EV（研发阶段） | EV、伺服 |
转矩产生机制与控制方程
SRM的转矩公式是什么？和PMSM完全不同？
根本不同。PMSM转矩是"电流×磁石磁通"，但SRM转矩来自电感变化率。线性区域（无磁饱和）的瞬时转矩是这样的：
其中 $I$ 是相电流 [A]，$L$ 是电感 [H]，$\theta$ 是转子角度 [rad]。这个公式有3个重点。
- 转矩与电流平方成正比 — 与电流方向无关，可用单向电流（单极）驱动。PMSM的正弦交流不必要。
- 仅当 $dL/d\theta > 0$ 时产生正转矩 — 需在电感增加区间通电，在减少前断电。
- 忽略磁饱和的近似式 — 实际上铁心饱和，这个线性式不够。
饱和会怎样？公式变吗？
会变。饱和时电感同时依赖角度和电流，变成 $L(\theta, I)$。此时需从共能（余能）推导精确转矩：
$\psi(\theta, i)$ 是链交磁通 [Wb]。此函数无法解析求得，实务标准方法是用FEM创建磁通-电流-角度的3D映射 $\psi(\theta, I)$，再数值微分求转矩。
SRM各相的电路方程如下：
右侧第2项是自感项（相当于 $L \cdot dI/dt$），第3项是反电动势项（$e = \omega \cdot d\psi/d\theta$，与转速成正比）。
电感特性曲线 L(θ,I)
$L(\theta, I)$ 的映射是什么样？想象不出…
低电流时接近梯形波。转子凸极与定子凸极对齐（对齐位置）时最大 $L_{\max}$，错开时（非对齐）最小 $L_{\min}$。
但电流增大后，对齐位置铁心饱和，$L_{\max}$ 大幅下降。而 $L_{\min}$ 由气隙主导变化不大。也就是说饱和时 $dL/d\theta$ 减小，转矩下降。这是SRM设计最大的困扰。
| 参数 | 典型值（8/6 SRM, 1kW级） | 说明 |
|---|---|---|
| $L_{\max}$（低电流） | 40〜80 mH | 对齐位置、线性区 |
| $L_{\max}$（额定电流） | 15〜30 mH | 饱和大幅下降 |
| $L_{\min}$ | 5〜10 mH | 非对齐位置、电流依存性小 |
| $L_{\max}/L_{\min}$比 | 3〜8（线性）/ 2〜4（饱和） | 比值越大高转矩 |
磁饱和与非线性B-H曲线
磁饱和具体在多少磁通密度发生？材料不同？
因使用的电磁钢板而异。SRM常用材料比如：
- 35A300（日本薄板）：饱和磁通密度 $B_s \approx 1.7$ T。鉄损少但饱和较快
- 50A470：$B_s \approx 1.8$ T。成本便宜但鉄损略大
- 10JNEX900（JFE的6.5% Si钢）：$B_s \approx 1.5$ T，高频鉄损极小。高速SRM适用
实务中1.5〜1.8 T开始B-H曲线变缓，相对磁导率从数千骤降到数十。FEM中忠实输入非线性B-H曲线是精度的关键。
转矩脉动与NVH（噪声振动）
SRM噪声具体是什么机制？真的比PMSM吵这么多？
吵。会发出"嗞"或"嘎"的金属声。原因很明确：径向电磁力使定子变形成椭圆。
SRM相通电ON/OFF时径向力急剧变化，激励定子铁心的固有振动模式（特别是椭圆模式=呼吸模式）。加振力频率与定子固有频率一致就共振，产生大振动和噪声。FEM的任务是：
- 电磁场解析求径向力时间波形
- FFT分析加振力频谱
- 结构模态解析求定子固有振动数
- 避免加振频率与固有频率重合的设计
也就是需要电磁场和结构的耦合解析？
是的。所以SRM的NVH解析流程是"电磁场FEM → 径向力提取 → 结构FEM（模态解析）→ 声学解析"这样的多阶段。JMAG和Maxwell有结构解析力数据导出功能，通常就用这个。
SRM作为"无磁石的EV电动机"备受关注
SRM完全不用永磁体，对中国稀土出口管制风险免疫。转子只是层叠钢板冲压重叠，制造成本仅为PMSM的1/2到1/3。Nidec在2023年开始量产家用洗衣机用SRM，通过AI最优电流波形实现传统SRM相比噪声降低60%。但Jaguar I-PACE的初期原型SRM虽然考虑过，却因NVH要求没能满足而回到PMSM。"安静的SRM"仍然是电动机业界的热门课题，电磁场FEM和控制仿真的精度是商品化的关键。
开关磁阻电动机（SRM）电磁场数值计算方法
电磁场FEM的定式化
SRM电磁场解析就是用FEM解麦克斯韦方程，对？但具体怎样定式化？
SRM解析频率（〜数十kHz）可忽略位移电流，用准静态近似（涡电流方程）。2D解析中只有矢量势 $A_z$ 的一个分量：
其中 $\nu = 1/\mu$ 是磁阻系数（透磁率倒数），$\sigma$ 是电导率，$\mathbf{J}_s$ 是施加电流密度。因 $\nu$ 依赖 $B$（也就是 $A$），方程是非线性的。
2D就够了？不需要3D？
SRM的基本转矩·电感计算2D-FEM精度足够，因为轴向结构均匀。但这些情况需要3D：
- 转子斜边结构 — 为降低转矩脉动而将转子拧转的设计
- 端部效应评估 — 线圈端的漏电感和端部涡电流
- 轴向磁通分量不能忽视 — 轴向磁通型SRM
FEM的弱形式用Galerkin法，用三角形单元（2D）或四面体/棱柱单元（3D）进行空间离散化。SRM电磁场FEM的关键要素技术如下。
| 要素技术 | 用途 | SRM重要度 |
|---|---|---|
| 节点单元（标量 $A_z$） | 2D截面解析 | 最基本。计算成本低 |
| 边单元（Nedelec单元） | 3D解析（矢量 $\mathbf{A}$） | 避免虚假模式必需 |
| 2阶单元（6节点三角形） | 气隙磁通精度提高 | 直接影响转矩精度。推荐 |
| 滑动网格 | 处理转子旋转 | 必需（下述） |
非线性牛顿-拉夫逊法
处理非线性B-H曲线，需在每个时间步进行牛顿-拉夫逊迭代？收敛没问题？
是的。在每个时间步进行牛顿-拉夫逊反复迭代。残差方程这样：
其中 $[S(\mathbf{A})]$ 是非线性刚度矩阵（依赖 $\nu(B)$），$[M]$ 是质量矩阵（涡电流项），$\{\mathbf{F}\}$ 是电流源矢量。雅可比矩阵是：
收敛的诀窍有3个：
- B-H曲线样条插值要光滑 — 折线近似会使 $d\nu/dB^2$ 不连续，雅可比振动
- 初值用前一步的解 — 从零开始收敛慢
- 引入阻尼（衰减系数） — $\Delta \mathbf{A}^{(k+1)} = \alpha \cdot [J]^{-1}\mathbf{R}$，取 $\alpha=0.5\sim0.8$
通常3〜8次反复达到 $\|R\|/\|R_0\| < 10^{-4}$ 收敛。超过这个次数说明B-H数据或网格有问题。
时间步进与旋转运动处理
转子旋转时网格怎么办？每步重新生成？
好问题。实务通常用两种方法：
1. 滑动网格法（Band法）
在气隙内建立圆周方向的"滑动带"，转子侧网格和定子侧网格独立生成。旋转时仅更新带上的节点对应关系。JMAG和Ansys Maxwell的标准方法。不需要网格重生成，但带上单元尺寸和旋转角度步长需协调，否则插值误差影响转矩。
2. 气隙重网格法
每步重新网格化气隙区域。比滑动带方法计算成本高，但更通用。COMSOL常用。
时间步长 $\Delta t$ 选择依据转速和通电角度如下：
$\theta_{\text{step}}$ 是一相通电角度 [rad]，$\omega$ 是角速度 [rad/s]。比如6000 rpm的8/6 SRM，$\theta_{\text{step}} = 15°$，$\omega = 200\pi$ rad/s，所以 $\Delta t \leq 4.2 \times 10^{-6}$ s（约4 μs）。
电路耦合解析
实际SRM由逆变器驱动吧？光用FEM不能得到电流波形啊？
很敏锐。要认真做SRM电磁场FEM，与驱动电路（逆变器）的耦合是绕不过去的。典型的SRM非对称半桥电路中，上下开关管和续流二极管状态决定了端子电压切换为 $+V_{dc}$、$0$、$-V_{dc}$。
耦合方法有两种：
- 直接耦合：FEM求解器内含电路方程。JMAG·Maxwell标准支持。电压源驱动FEM计算电流。精度最高
- 表格基础耦合：FEM预先计算 $\psi(\theta, I)$ 映射，嵌入Simulink控制模型。计算快得多但忽略了相互电感和过渡饱和效果
初期设计用表格法，详细设计·NVH用直接耦合，这是现场普遍做法。
开关磁阻电动机（SRM）电磁场实务应用
SRM解析总体流程
SRM电磁场解析从头到尾怎么流程？想看整体情况。
SRM设计解析大体是这样流程：
- 形状参数确定 — 定子/转子极数（例：8/6, 12/8）、凸极宽、气隙长、堆叠长、绕组匝数
- 2D断面模型创建 — CAD或FEM工具参数化造型
- 材料设定 — B-H曲线（非线性）、铁损系数（Steinmetz常数）、铜电阻率
- 网格生成 — 气隙和凸极先端优先细分
- 静磁场解析（先做简单的） — 各转子角度计算 $L(\theta, I)$ 映射
- 过渡解析 — 电路耦合计算实际电流波形·转矩波形·铁损
- 后处理 — 评估转矩脉动率、效率、径向力谱
- 参数扫描/优化 — 凸极宽·匝数·控制角的最优化
不能直接从过渡解析开始？
不能。先静磁场解析算出 $L(\theta, I)$ 映射，用手工计算或MATLAB验证"大概的性能"是必须的。直接过渡解析要花好多个小时，错误发现太晚。现场一般是"静磁场2天验证 → 过渡解析3天 → 参数扫描1周"这样的时间感。
网格策略（气隙·凸极）
SRM的网格特别要注意的地方是哪里？
SRM网格失误最痛的是气隙。这里的磁通密度分布直接决定转矩计算精度。参考如下：
| 区域 | 推荐单元尺寸 | 理由 |
|---|---|---|
| 气隙（径向） | 气隙长的1/3〜1/5（例：0.5mm气隙则0.1〜0.17mm） | 麦克斯韦应力张量精度确保 |
| 气隙（周向） | 凸极先端圆弧40分割以上 | 转子旋转时转矩波形平滑度 |
| 凸极先端（定子/转子） | 0.2〜0.5 mm | 局部饱和的精确捕捉 |
| 轭部 | 1〜3 mm | 磁通密度梯度缓，可粗 |
| 线圈区 | 1〜2 mm | 电流密度均匀，可粗 |
| 外部空气区 | 5〜20 mm（逐级增大） | 远场精度要求低 |
外部空气区要扩到多远？无限远怎么办…
实务上定子外径的2〜3倍半径到外周边界，施加零Dirichlet条件 $A_z = 0$ 就通用。更小会产生"看不见的墙"磁通反射导致转矩异常。更优雅的做法是用Kelvin变换（无限单元），JMAG和Maxwell都有这功能。
边界条件与对称性利用
8/6 SRM全周模型太浪费，能用对称性？
能用。8/6 SRM中定子极数和转子极数的最大公约数 = GCD(8,6) = 2，所以1/2模型（180°分）是最小对称单位。但要注意：
- 周期边界条件（反对称/同型对称）要正确设。SRM一般用反对称条件 $A_z(\theta) = -A_z(\theta + 180°)$
- 看相互电感时1/2模型够，考虑转子偏心和制造偏差则需全模型
- 12/8 SRM则 GCD(12,8) = 4，可用1/4模型（90°分）
后处理与转矩计算验证
FEM中转矩计算有几种方法？选哪个？
主要有3种。为保证可信度至少用2种方法验证一致是铁律：
| 方法 | 原理 | 精度 | SRM注意事项 |
|---|---|---|---|
| 麦克斯韦应力张量法 | 气隙中的面积分 | 网格依存性高 | 积分面设在气隙中心。网格粗时误差大 |
| 虚功原理 | 共能的角度微分 | 高 | 需要两个转子角度的解。扰动量 $\delta\theta$ 选择要注意 |
| Erying-Tunks法 | 单元级局部力计算 | 极高 | JMAG标准。网格依存性小，推荐 |
麦克斯韦应力张量法网格敏感让我很困扰…
是的。所以网格收敛性检查必须做。转矩值随网格细化而不再变化（收敛）为止，至少3个网格密度阶段计算。特别是气隙单元层数从2层→4层→6层变化时转矩变化在1%以内就合格。
SRM噪声问题 — "刺耳音"怎样消除
SRM实际运转会产生切换时"嗞嗞""啪啪"的刺耳电磁噪声，这一直是SRM搭载EV的最大障碍。噪声根源是相切换时径向电磁力导致定子椭圆变形（呼吸模式）。对策包括：(1) 增加极数，提高加振次数（12/8→16/12等）；(2) 电流波形从方波改正弦风格，实现软开关；(3) 定子外周加肋条，移动固有频率等。现场需要结合多个对策，通过FEM预测加振力频谱，同时优化结构和控制，找到最优解。
开关磁阻电动机（SRM）电磁场软件对比
SRM设计对应工具对比
SRM解析用什么软件？费用多少？
主要有4个。SRM解析中"电路耦合"和"非线性收敛稳定性"是选工具的关键：
| 工具 | 开发方 | SRM强项 | 弱点 | 许可证 |
|---|---|---|---|---|
| JMAG-Designer | JSOL（日本） | 有SRM模板，非线性收敛稳定，铁损模型丰富，NVH联动充实 | 海外文献多用Maxwell | 节点锁定/浮动（年300万〜） |
| Ansys Maxwell | Ansys（美国） | Simplorer/Twin Builder制控耦合强大，自适应网格自动生成 | 初始网格粗时转矩脉动有误差 | Ansys套件内（年500万〜） |
| COMSOL Multiphysics | COMSOL AB（瑞典） | 电磁-结构-声学完全耦合在一工具。自定制性高 | 无电动机专用模板。电路耦合略复杂 | 模块制（需AC/DC模块，年200万〜） |
| Motor-CAD | Ansys/Motor Design Ltd（英国） | 包括SRM的各种电动机设计专化。高速参数扫描 | FEM内装但网格细控不可。NVH无 | Ansys套件或单体（年150万〜） |
开源方案没有？大学经费吃不消…
有的。FEMM（Finite Element Method Magnetics）是2D磁场解析专用免费工具，SRM的静磁场解析和电感计算用它完全可以。支持Lua/Python脚本自动化。但无过渡解析和电路耦合，需自己用MATLAB/Python写外部循环。还有Elmer FEM（芬兰CSC开发）是开源3D电磁过渡解析工具，但SRM的滑动网格实现还不成熟。
用途别工具选择指南
| 用途 | 推荐工具 | 理由 |
|---|---|---|
| 大学研究（论文级） | FEMM + Python | 免费，脚本驱动再现性好 |
| 初期设计·参数扫描（产业） | Motor-CAD | 数分钟运100个案例，速度快 |
| 详细设计（转矩精度优先） | JMAG / Maxwell | 非线性FEM+电路耦合可信 |
| NVH解析（噪声对策） | JMAG → Nastran/Actran联动 | 径向力→结构模态→声学联动 |
| 多物理耦合（含发热） | COMSOL | 电磁-热-结构在一工具完成 |
| 控制最优化 | Maxwell + Simplorer / JMAG-RT + Simulink | 控制闭环与FEM耦合 |
"Simulink + FEM耦合"何以成为SRM控制优化的标配
SRM设计不能只靠电磁场FEM，需加入电力变换器（逆变器）和控制算法的系统级仿真。原因是相电流波形最优化同时影响噪声和效率。JMAG的"JMAG-RT"和Ansys Maxwell的"Simplorer"都提供了Matlab/Simulink接口，将FEM计算出的$\psi(\theta, I)$表格嵌入Simulink控制模型这一工作流已成标准。这种"电磁场FEM + 控制仿真耦合"是当前SRM控制优化的必需品。
开关磁阻电动机（SRM）电磁场故障排除
非线性收敛失败
老师，SRM过渡解析跑牛顿-拉夫逊没收敛就停了…怎么办？
SRM非线性收敛失败的原因95%是以下4个：
| 原因 | 症状 | 对策 |
|---|---|---|
| B-H曲线不连续·非单调 | 残差振动不收敛 | 用3阶样条光滑补间B-H数据。全域验证 $dB/dH > 0$ |
| 时间步长太大 | 电流急变时发散 | 将 $\Delta t$ 减半重算。通电ON/OFF时启用自动步长控制 |
| 气隙网格太粗 | 转矩振荡，出现负值 | 气隙径向最少3层单元 |
| 初始磁化状态不协调 | 最初几步发散 | 无通电（$I=0$）状态下逐步升电流的"升压"过程 |
B-H数据直接用厂商规格书可以？
规格书点数少（约10点）直接用会变折线，微分不连续导致雅可比振动。最少细化到30〜50点光滑曲线。高磁场侧（$B > 2.0$ T）数据常缺，用 $\mu_0$ 外推是标准做法。JMAG有B-H自动外推·平滑化功能，用这个很安心。
转矩波形异常振荡
计算收敛了，但转矩波形锯齿状，一点不光滑。这是对的结果？
SRM本来转矩脉动大，但要分"物理脉动"和"数值噪声"。区别方法简单：网格细化到2倍，再算一次。锯齿减少就是网格问题，不变就是物理脉动。
数值噪声主因是：
- 麦克斯韦应力张量法的网格依存性 — 改用虚功法或Erying-Tunks法
- 滑动带单元不协调 — 带上单元尺寸和旋转步长要同步。旋转角等于一个单元最理想
- 时间步长过粗 — 电气角度0.5°以下细度
初学者常见错误
老师，你见过的SRM解析新手最常犯的错误是什么？
排名前5的新手错误是这样的：
| # | 错误内容 | 结果 | 正确做法 |
|---|---|---|---|
| 1 | B-H曲线用线性（固定 $\mu_r$） | 转矩变成实测的2〜3倍 | 必须用非线性B-H曲线 |
| 2 | 空气区太小（定子外径的边界） | 磁通外周反射，转矩异常 | 外径的2〜3倍范围确保空气 |
| 3 | 通电角（turn-on/turn-off角）设错 | 产生负转矩，平均转矩猛跌 | 仅在 $dL/d\theta > 0$ 的区间通电 |
| 4 | 对称边界条件弄反（同型↔反对称） | 磁通分布物理上不可能 | 根据通电模式选对对称/反对称 |
| 5 | 忽略铁损，效率估计过高 | 实机比仿真效率低10〜20% | 用Steinmetz或FEM铁损模型（JMAG的iron loss求解器等） |
SRM结构简单的分，仿真"设置的正确性"直接决定精度。先用FEMM的静磁场跑一遍8/6 SRM试试。手工计算验证理解会深很多。
很好
更详细
报告