弱磁控制的电磁场仿真
弱磁控制电磁场的理论基础
什么是弱磁控制
弱磁是指故意削弱磁石的力量吗？ 不是很浪费吗？
乍看起来浪费，但基础速度以上要旋转的话，会撞上逆变器的电压限制。通过在d轴注入负向电流来抵消磁石的磁通，从而抑制反电动势，在逆变器电压限制范围内实现高速旋转。这对电动汽车的高速公路行驶是必需的技术。
电动汽车难道不是用齿轮箱改变速度吗？
很多电动汽车没有变速器，也就是所谓的单档。在内燃机汽车中，会换到5档或6档进行高速行驶，而电动汽车的电动机则通过弱磁控制来实现这一点，能够达到定额转速的3到4倍。换句话说，弱磁控制是让电动机进行"虚拟变速"的技术。
明白了，就是用改变电流的方法代替齿轮！ 但如果削弱磁石的力量，扭矩不会下降吗？
完全正确。在弱磁控制区域，扭矩随转速反比下降。也就是说，输出（=扭矩×转速）基本保持常数的"恒功率区域"。在高速公路100 km/h巡航时所需扭矩很小，这样的性能就足够了。
弱磁控制（Field Weakening Control）是永久磁石同步电动机（IPMSM/SPMSM）在超过基础速度的高速动作中实现的电流控制技术。通过在d轴注入负向电流（$i_d < 0$）来电气消除永磁石磁通 $\psi_m$，抑制反电动势，使得在逆变器输出电压约束范围内实现高速旋转。
dq轴电压方程和约束条件
具体来说，控制的约束用什么数学公式表示呢？
从IPMSM的dq轴电压方程出发。假定稳态（$di/dt = 0$），得到这样的式子：
其中 $R_s$ 是绕组电阻，$L_d, L_q$ 是dq轴电感，$\omega_e$ 是电气角速度，$\psi_m$ 是永磁体磁通链数。在高速区域 $R_s$ 的电压降较小，可以忽略，则**电压限制条件**可以整理为：
同时，逆变器的电流额定值施加**电流限制条件**：
也就是说，电压和电流都有上限，不能超出这个范围运行，对吧？
是的。这里很重要的一点是，在 $i_d$-$i_q$ 平面上看电压限制式，它是一个椭圆。转速越高，$\omega_e$ 越大，椭圆就越缩小。所以在高速区域，可工作的电流范围会越来越窄。
电压限制圆和电流限制圆
将电压限制条件用 $i_d$-$i_q$ 平面重新写出，得到：
这是中心 $(-\psi_m/L_d, \, 0)$、长短轴分别为 $V_{\max}/(\omega_e L_d)$ 和 $V_{\max}/(\omega_e L_q)$ 的椭圆。在IPMSM中由于 $L_d < L_q$，所以d轴方向是短轴。随着转速上升，椭圆缩小，与电流限制圆 $i_d^2 + i_q^2 = I_{\max}^2$ 的交叉区域会发生变化。
为什么椭圆的中心不在原点，而是在 $(-\psi_m/L_d, 0)$ 呢？
因为永磁体的磁通 $\psi_m$ 始终存在于d轴方向，所以电压会被"消耗"。即使d轴电流为零，磁石磁通产生的反电动势也会出现。所以椭圆的中心向 $i_d = -\psi_m/L_d$ 偏移。弱磁控制本质上就是把工作点向这个椭圆中心移动。
MTPA控制和MTPV控制
经常听到MTPA和MTPV这些词，有什么区别吗？
简单地说，是按速度区域来应用的两种最优电流策略。
MTPA（单位电流最大扭矩）： 在电流限制圆内部，选择给定电流幅值下产生最大扭矩的电流矢量角 $\beta$ 的策略。IPMSM的扭矩式：
为了活用磁阻扭矩项 $(L_d - L_q) i_d i_q$，在 $i_d < 0$ 的区域动作。MTPA轨迹在 $i_d$-$i_q$ 平面上呈双曲线形的曲线。
MTPV（单位电压最大扭矩）： 在电压限制椭圆上选择产生最大扭矩的工作点的策略。高速区域在达到电压限制后使用，沿着椭圆边界移动工作点。MTPV轨迹的条件式是：
那么，随着速度增加，工作点会沿着MTPA轨迹→电压限制椭圆上→MTPV轨迹这样移动吗？
完全正确。低速区域沿MTPA轨迹实现最高效运行，基础速度附近达到电压限制后开始沿椭圆移动，更高速区域过渡到MTPV轨迹。灵活地切换这三个区域是控制的关键，在CAE中需要对各种速度-扭矩组合计算最优电流点，进而构建"效率图"。
减磁风险的物理
听说如果d轴电流太大，磁石会坏掉，是真的吗？
与其说是坏掉，不如说是"减磁"。钕磁石有去磁曲线（B-H曲线的第二象限），当工作点超过膝点（曲线的折点）时，磁力会不可逆地下降。弱磁控制中由于要注入大的负向d轴电流，磁石内部的磁通密度有可能接近膝点，存在减磁的危险。
温度也有关系吧？ 夏天的时候好像很危险...
正是这里最危险。钕磁石的保磁力随温度上升急剧降低。20°C时有2000 kA/m的保磁力，到150°C时可能下降到一半以下。所以在CAE中，一定要在最恶劣条件（最高温度×最大d轴电流×过渡电流超调）下确认磁石工作点。
减磁判定通过FEM求出的磁石内部磁通密度分布与对应温度的去磁曲线进行对照而进行。磁石角部（边角）由于反磁场集中，局部磁通密度容易降低，需要特别注意。
特斯拉Model 3的电动机设计是"弱磁控制的教科书"的原因
特斯拉Model 3的后电动机（IPMSM类型）定额约5000转，最高转速约18000转。也就是说，基础速度的3.6倍在弱磁控制下运转。这之所以可能，是因为通过磁石嵌入式结构（V字配置）活用了磁阻扭矩，以及SiC逆变器实现了高电压利用率。在设计阶段，FEM对弱磁时的减磁余裕进行了大量分析。日本汽车OEM也采用了类似的设计流程，JMAG和Ansys Maxwell的FEM-回路耦合分析发挥了中心作用。
弱磁控制电磁场的数值计算方法
FEM电动机建模基础
用FEM仿真弱磁控制，具体进行什么计算呢？
大致分为两种方法。一种是进行2D截面FEM磁场解析，计算各旋转角度位置的扭矩、磁通链数、损耗。另一种是把FEM和驱动回路（逆变器）耦合，直接仿真弱磁控制算法本身。
电动机电磁场分析中，使用矢量势 $\mathbf{A}$ 求解准静态磁场方程：
其中 $\nu$ 是磁阻率（$\nu = 1/\mu$）、$\mathbf{J}_s$ 是线圈电流密度、$\sigma \partial \mathbf{A}/\partial t$ 是涡流项、$\mathbf{M}_r$ 是永磁体的剩余磁化矢量。
2D分析就够了吗？ 电动机是三维物体啊。
很敏锐。用硅钢片堆叠的电动机，轴向电流占主导，所以2D截面分析通常能得到足够精度。通过补正系数来处理斜切（捩转）和端部线圈的影响。但要精确评估3D特有的端部效应，比如斜切对扭矩纹波的降低作用或端部线圈漏磁感，需要3D分析。
回路耦合分析
FEM和回路的耦合具体怎么做呢？
从FEM的各线圈区域提取磁通链数 $\psi$ 和反电动势，然后和外部回路的基尔霍夫定律联立求解。线圈电流影响磁场，磁场影响线圈反电动势，二者相互影响，所以需要每个时间步长反复迭代求解。
FEM-回路耦合的基本方程组如下：
其中 $K_{ff}$ 是FEM的磁场刚度矩阵，$K_{cc}$ 是回路的阻抗矩阵，$K_{fc}$, $K_{cf}$ 是FEM与回路的耦合项，$A$ 是节点矢量势，$I$ 是线圈电流，$V_s$ 是逆变器输出电压。
逆变器的PWM开关也要一起算吗？ 计算量好像会很大...
取决于目标。只想看弱磁区域的平均扭矩和效率的话，用基波电压源就可以。如果要分析扭矩纹波或EMC，就需要以PWM载波时间刻度（数微秒）进行求解。实务中通常先用基波分析制作效率图，然后只对关心的工作点进行PWM耦合详细检讨。
损耗计算（铁损、铜损、磁石涡流损）
要正确评估弱磁区域的效率，需要精确计算以下损耗成分：
| 损耗成分 | 计算式 | 弱磁时的变化 |
|---|---|---|
| 铜损 $P_{cu}$ | $\frac{3}{2}R_s(i_d^2 + i_q^2)$ | $|i_d|$ 增加而增大（最大损耗因子） |
| 滞后铁损 $P_h$ | $k_h f B_m^{\alpha}$ | 频率 $f$ 增加而增大，$B_m$ 降低而减小 |
| 涡流铁损 $P_e$ | $k_e f^2 B_m^2$ | 与 $f^2$ 成比例急剧增加（主要损耗） |
| 磁石涡流损 | FEM时间域分析直接计算 | 高次谐波磁通增加而增大 |
| 机械损 | 轴承摩擦 + 风损 | 与转速的2～3次方成比例 |
铁损计算式中的 $B_m$ 用什么呢？ 与磁石的 $B_r$ 不同吧？
$B_m$ 是各单元处的磁通密度振幅。进行1电气周期的瞬态分析，获取各有限元的磁通密度波形。从中用FFT分离基波和谐波成分，分别计算各频率的铁损。再对全部单元积分，这就是FEM型铁损计算。
效率图的构建方法
效率图怎么制作呢？ 要对所有速度×扭矩组合进行FEM解析吗？
那样的话计算量太多了。实务中先用FEM"参数化扫描dq轴电流"进行磁场分析。比如 $i_d$ 从-200A到0A、$i_q$ 从0A到300A各扫20点，就是400个工况的静磁场分析。每个工况得到扭矩、磁通链数、铁损，建成查找表。然后用这个查找表，为任意速度-扭矩指令找到最优电流。
效率图构建的步骤如下：
- FEM参数分析： 在 $(i_d, i_q)$ 网格点进行静磁场解析 → 生成扭矩 $T(i_d, i_q)$、$\psi_d(i_d, i_q)$、$\psi_q(i_d, i_q)$、铁损 $P_{fe}(i_d, i_q)$ 的图表
- 最优电流搜索： 对各速度-扭矩指令，在满足电压限制、电流限制的约束下，搜索损耗最小的 $(i_d^*, i_q^*)$
- 效率计算： 绘制 $\eta = P_{out} / (P_{out} + P_{cu} + P_{fe} + P_{mech})$ 的等效率线（如95%、90%、85%）
非线性收敛和磁饱和处理
硅钢片的B-H曲线是非线性的。FEM中怎么处理呢？
用牛顿-拉夫逊法进行非线性迭代。每次迭代时，从当前磁通密度查表B-H曲线，更新磁导率，重新构建刚度矩阵并求解。弱磁区域中，d轴电流导致磁石附近的铁心局部脱饱和，收敛有时会改善。而q轴大电流会使齿部进入深饱和，收敛可能恶化。
非线性收敛的判定标准通常是残差范数的相对值，建议 $||R||/||R_0|| < 10^{-4}$。磁饱和深度大的情况下，可配合直线搜索或自适应松弛提高鲁棒性。
弱磁控制电磁场的实务应用
弱磁解析的工作流程
实际进行弱磁仿真时，按什么步骤进行呢？
一个典型的工作流程我来详细说明。
第1步：几何建模
- 绘制定子、转子、永磁体、空气域（包括气隙）的截面形状
- 利用对称性：8极48槽可用1/8模型（1极）进行分析
- 定义旋转界面用于滑动网格
第2步：材料设置
- 硅钢片：非线性B-H曲线 + 铁损系数（Steinmetz公式或Bertotti公式）
- 永磁体：剩余磁化 $B_r$、矫顽力 $H_{cJ}$、温度系数、去磁曲线（第2象限）
- 导体：铜的导电率（温度相关）
第3步：参数励磁分析
- 定义d轴电流 $i_d$ 和q轴电流 $i_q$ 的网格点（如20×20 = 400点）
- 每点执行1电气周期的瞬态分析（旋转角度步长：电气角1～2°）
- 提取扭矩平均值、磁通链数（$\psi_d$、$\psi_q$）、铁损
第4步：效率图构建
- 从查找表数据进行最优电流搜索（速度×扭矩各点损耗最小化）
- MTPA/弱磁/MTPV区域的自动判别
- 绘制等效率线（如95%、90%、85%）的等高线图
第5步：减磁风险评估
- 在最大d轴电流×最高温度条件下进行FEM分析
- 磁石各部位的磁通密度分布与去磁曲线对比
- 减磁余裕（到膝点的距离）的定量评估
电动机电磁场分析的网格策略
网格中最需要注意的地方在哪里？
电动机特有的注意点有几个。
- 空气气隙： 配置最少3～5层单元。气隙长0.5 mm需要0.1 mm以下的单元尺寸。这是对扭矩精度影响最直接的关键区域
- 磁石边缘： 为了进行局部减磁判定，磁石角部需要细化网格
- 槽口部分： 磁通密度变动较大的区域。影响齿槽扭矩精度
- 转子桥接： 磁饱和深的薄壁部分。至少需要2～3层单元
- 旋转界面： 滑动网格的分割数与旋转角度步长的一致性需要注意
| 区域 | 推荐单元尺寸 | 单元类型 | 原因 |
|---|---|---|---|
| 空气气隙 | 0.05～0.15 mm | 三角形1阶 | 扭矩精度·麦克斯韦应力张量 |
| 磁石内部 | 0.3～1.0 mm | 三角形1阶 | 减磁评估所需分辨率 |
| 齿部·轭 | 0.5～2.0 mm | 三角形1阶 | 饱和·铁损计算所需 |
| 外部空气域 | 5～20 mm | 三角形1阶 | 漏磁通衰减的表现 |
边界条件和对称性的应用
8极48槽电动机使用1/8模型时，边界条件怎么设置呢？
周向两端应用**周期边界条件**（反周期边界）。相邻极的磁极相反，所以用 $A(0°) = -A(45°)$ 这样的反周期条件。外周边界用 $A = 0$（迪利克雷边界）防止磁通泄漏。外部空气域在定子外径的1.5～2倍就足够了。
减磁余裕检查的实务步骤
减磁检查的实务步骤请教一下。判定基准是什么？
步骤是这样的。首先定义最恶劣条件：
- 磁石温度：最高使用温度（如150°C～180°C）
- d轴电流：最大（包括短路故障时）
- q轴电流：零（最严格地评估d轴电流的影响）
在这个条件下运行FEM分析，确认磁石内部的最小磁通密度 $B_{\min}$。然后与对应温度去磁曲线的膝点 $B_{knee}$ 比较：
- $B_{\min} > B_{knee} + \Delta B_{margin}$ 则OK
- 余裕 $\Delta B_{margin}$ 通常确保 0.05～0.1T
- 若 $B_{\min}$ 低于 $B_{knee}$，通过改变磁石等级（如N改为SH）、限制d轴电流上限、改变磁石形状（增加磁通障碍）等方式对策
"弱磁减磁"——电动汽车开发现场的实际失败案例
某国产电动汽车开发项目因弱磁控制电流限制值设定错误，过大的d轴电流导致钕磁石部分减磁的事例。减磁表现为扭矩下降、效率恶化，产品回收前才原因特定。根本原因是"弱磁区域的磁石工作点未包含在仿真中"。现在的标准流程是在FEM上全温度范围确认弱磁时磁石工作点，并作为设计规范管理与去磁曲线的余裕（demagnetization margin）。
弱磁控制电磁场的软件比较
弱磁解析对应工具比较
弱磁仿真可以用什么软件？
从弱磁解析的观点对电动机电磁场分析的主要工具进行比较。
| 功能 | JMAG | Ansys Maxwell | Motor-CAD | COMSOL |
|---|---|---|---|---|
| 2D FEM瞬态分析 | ◎ | ◎ | -（LPM） | ○ |
| 3D FEM分析 | ◎ | ◎ | - | ○ |
| FEM-回路耦合 | ◎ | ◎（Twin Builder） | - | ○ |
| 效率图自动生成 | ◎（标准功能） | ○（需脚本） | ◎（主功能） | △ |
| 减磁解析 | ◎（温度耦合） | ○ | △（简易） | ○ |
| MTPA/MTPV自动计算 | ◎ | ○ | ◎ | △ |
| 磁石温度相关 | ◎ | ○ | ◎ | ○ |
| 计算速度（效率图） | 中（数小时～） | 中（数小时～） | ◎（数分钟） | 低 |
| 日语支持 | ◎（国产） | ○ | ○ | ○ |
JMAG的弱磁解析
据说日本汽车厂商多用JMAG，弱磁解析功能充实吗？
JMAG针对电动机设计特别下功夫，弱磁相关的功能非常充实。用JMAG-Designer的"效率图"功能，从dq轴电流参数分析到效率图生成基本可以半自动执行。特别是减磁分析，能温度耦合进行，可以考虑磁石温度分布后评估各部位的减磁余裕。
JMAG弱磁解析的工作流：
- 模板模型创建： 用电动机模板功能定义几何形状
- dq轴参数分析： 用"特性解析"模块自动扫描 $(i_d, i_q)$ 网格
- 效率图生成： 用损耗计算结果进行最优电流搜索和效率等高线制作
- 减磁检查： 用"减磁耐力分析"自动评估磁石工作点与去磁曲线的余裕
Ansys Maxwell的弱磁解析
Ansys Maxwell怎么样呢？
Maxwell的强点是RMxprt（分析设计）→ Maxwell 2D/3D（FEM）→ Twin Builder（系统仿真）→ Ansys Motor-CAD（效率图）的联动。RMxprt设计的模型自动导入Maxwell进行FEM验证，Twin Builder可进行包含逆变器PWM控制的时间域仿真。不过效率图的自动生成没有JMAG那么简单，实务中多数需要用Python脚本（PyAEDT）自动化。
Motor-CAD的系统级设计
Motor-CAD不是FEM吧？ 也能做弱磁分析吗？
Motor-CAD是集中参数模型（LPM）基础，但在效率图高速计算上有压倒性优势。一个速度×扭矩点的计算以秒计，所以可以在改变冷却条件下即时更新效率图，或筛选100个转子形状。这在概念设计阶段是无敌的武器。但磁石角部局部减磁评估不用FEM精度达不到，所以最终验证还是要交给JMAG或Maxwell。
JMAG vs Motor-CAD——"精密"和"速度"的二刀流
在弱磁区域的电动机解析工具中，日本国内JMAG（JSOL制，FEM）广为好评。弱磁时的铁损、铜损、减磁余裕的温度依赖计算工作流充实，被汽车OEM和一级供应商广泛采用。Motor-CAD（Ansys制）强在热-电-机械的快速耦合分析，系统级效率图计算速度压倒性更快。详细设计用JMAG，概念设计用Motor-CAD的分工比较普遍。近来出现的混合工作流将JMAG分析结果导入Motor-CAD，组合了FEM精度和LPM速度，也开始普及。
弱磁控制电磁场的故障排除
弱磁仿真的典型问题
老师，弱磁仿真容易出现什么问题？
好问题。弱磁特有的陷阱我来总结一下。
| 症状 | 原因 | 对策 |
|---|---|---|
| 高速域效率图扭矩不变为零 | 电压限制条件实现错误。未考虑SVPWM的电压利用率 $1/\sqrt{3}$ | 正确设置 $V_{\max} = V_{dc}/\sqrt{3}$ |
| FEM和等效回路扭矩不一致 | FEM中 $L_d$, $L_q$ 随电流变化，等效回路却用常数 | 将FEM得到的 $L_d(i_d, i_q)$、$L_q(i_d, i_q)$ 图表反映到等效回路 |
| 减磁余裕过于严格 | 用单一最高温度评估磁石（实际有温度分布） | 用热分析结果的温度分布。磁石中心部比周边更冷 |
| 弱磁区铁损被低估 | 仅计算基波铁损，忽视谐波铁损 | 用FFT分离各次谐波，逐次计算铁损再积分 |
| 非线性收敛在弱磁条件恶化 | 转子桥深饱和与磁石脱饱和同时出现，磁导率空间变化大 | 启用牛顿-拉夫逊直线搜索，初始解用低电流工况 |
| MTPV区域效率图缺失 | $\psi_m/L_d > I_{\max}$（椭圆中心在电流限制圆外），此设计不存在MTPV轨迹 | 重新审视设计参数，使 $\psi_m/L_d \leq I_{\max}$，实现宽定出力域 |
扭矩纹波和谐波问题
听说弱磁域扭矩纹波会增大，为什么？
弱磁控制中d轴电流增加，定子齿部的饱和模式改变。通常运转时的谐波成分与之不同，导致齿槽扭矩和扭矩纹波的波形变化。特别是IPMSM由于磁阻扭矩的空间谐波成分会改变，基础速度以下纹波很小，但弱磁域会急剧增大的例子不少。
这能从仿真预测吗？
当然可以。在弱磁条件（$i_d$ 增大）下进行FEM瞬态分析，计算扭矩波形，纹波直接算出。但要出精度需要充分降低旋转角度刻度（电气角0.5～1°以下）和充分细化空气间隙网格。考虑用斜切最优化对策时需要3D分析。
调试检查清单
仿真结果不对时，有系统的检查步骤吗？
弱磁仿真的系统化调试步骤教你。
- 单位系统检查： $\psi_m$ [Wb]、$L_d$/$L_q$ [H]、$i_d$/$i_q$ [A]、$\omega_e$ [rad/s] 是否整合一致
- 空载反电动势验证： FEM空载分析的反电动势是否与 $E = \omega_e \psi_m$ 一致。不一致表明磁石材料定义有问题
- 电感合理性： FEM计算的 $L_d$、$L_q$ 是否在教科书范围（$L_q/L_d$ = 2～8）内
- 扭矩理论值比对： FEM扭矩与 $T = \frac{3}{2}p[\psi_m i_q + (L_d - L_q)i_d i_q]$ 比较。大偏差说明饱和影响或定式问题
- 基础速度电压确认： 基础速度时的电压是否到达 $V_{\max}$。未达表明 $V_{\max}$ 定义或巻线电阻有问题
- 减磁检查温度： 去磁曲线是否用了正确温度。是不是用了20°C数据而评估150°C条件
按这个步骤逐项排查，大部分问题都能解决了。谢谢！
是的。弱磁解析是电气、磁场、热3个物理相互作用的复杂问题，一个个分离原因是铁则。"全部一起改"会陷入死循环。
帮助
更详细
错误