电机热管理分析 — 电磁-热耦合的理论和实务
电机热管理的理论基础
为什么需要电磁-热耦合
老师，电机的温度分析需要与电磁分析分开进行吗？我以为电磁和热是可以分离计算的…
很好的问题。结论是需要进行耦合分析。之所以这样做，是因为电磁和热之间会双向影响。
具体来说，电机的发热源有三种——铜损、铁损、机械损。这些损失导致温度上升，然后会发生什么呢…
- 铜的电阻率增加 → 铜损进一步增加（正反馈）
- 永磁体的剩余磁通密度 $B_r$ 下降 → 转矩降低，为了输出相同的功率需要更大电流，进一步发热
- 铁损特性也会改变 → 硅钢片的B-H曲线随温度变化
也就是说，"电磁分析计算损失 → 热分析计算温度 → 用该温度更新材料物性再做电磁分析"这样的循环是必需的。
哇，这么相互影响啊！温度上升时铜的电阻增加我能理解，但磁石也受影响吗…
对。特别是钕铁硼（NdFeB）磁石很厉害。温度系数约为 $-0.12\%/°\text{C}$，也就是说常温时 $B_r = 1.3\,\text{T}$ 的磁石，到150°C时就会降到约 $1.14\,\text{T}$。而且超过150°C时不可逆失磁风险急剧上升。一旦不可逆失磁发生了，冷却后也恢复不了。
不可逆…很吓人呢。EV的电机比如说在高速公路爬坡时连续高负荷，那不会出问题吗？
正是这种工况是设计的最大课题。所以通过电磁-热耦合分析来"最坏条件下的磁石温度"能得到准确预测，这是电机设计最重要的工作之一。
损失的分类和发热机制
铜损、铁损、机械损，各自占的比例大概是多少呢？
比例是随电机种类和运行条件变化的，但典型的EV用IPMSM情况大概是这样。
| 损失类型 | 发生位置 | 占全损失比例（参考） | 温度依存性 |
|---|---|---|---|
| 铜损 $P_{cu}$ | 定子绕组 | 50～70% | 温度上升时增大 |
| 铁损 $P_{fe}$ | 定子铁芯、转子铁芯 | 15～30% | 温度上升时略减小 |
| 机械损 $P_{mech}$ | 轴承、气隙（风损） | 5～15% | 转速相关 |
| 磁石涡流损 | 永磁体 | 1～5% | 高次谐波相关 |
铜损占一半以上啊！铜损的公式是什么样的？
基本是 $P_{cu} = I^2 R$，但考虑温度依存性就是这样：
其中 $R_0$ 是基准温度 $T_0$（通常20°C）的电阻，$\alpha_{Cu} \approx 0.00393\,/°\text{C}$ 是铜的电阻温度系数。
比如，某电机在20°C时铜损为10W，温度达到150°C时：
$P_{cu}(150) = 10 \times [1 + 0.00393 \times (150 - 20)] = 10 \times 1.511 \approx 15.1\,\text{W}$
也就是说铜损增加了约50%。这个在没有耦合分析时就会被遗漏。
增加50%这么大啊… 那铁损怎么计算？
铁损历来用Steinmetz公式，改进版的iGSE（改进的广义Steinmetz方程）在业界最常用。
第1项：滞后损（$k_h$：滞后系数，$\beta \approx 1.6 \sim 2.0$）
第2项：经典涡流损（$k_e$：涡流系数）
第3项：异常涡流损（$k_{ex}$：磁壁移动相关）
系数有3个，这些怎么求啊？实验吗？
对，从钢板厂商的数据（Epstein试验、环形芯试验的结果）反演。比如新日铁住金的35H300硅钢，50Hz、1.5T时铁损为 $3.0\,\text{W/kg}$ 这样的数据，就从这个反演出系数。
支配方程
损失分清楚了，那温度分布怎么计算？
基本就是三维非定常热传导方程。电磁分析得到的损失密度 $q_{loss}$ 作为发热项（源项）输入。
$\rho$: 密度 [kg/m³]、$c_p$: 比热 [J/(kg·K)]、$k$: 热导率 [W/(m·K)]、$q_{loss}$: 体积发热密度 [W/m³]
关键是 $q_{loss}$ 本身是温度 $T$ 的函数。温度升高，铜损变大，铁损特性也变。所以 $q_{loss}(\mathbf{x}, T)$ 这样写的意思就是这个。
啊，发热量本身是温度的函数，所以不能一次解出来啊。
其实。加上电机各部件之间的热阻也很重要。特别是这些界面会大大影响温度分布：
| 界面 | 热阻的原因 | 典型值 |
|---|---|---|
| 绕组-槽壁 | 浸渍漆 + 空气层 | $R_{th} \approx 0.01 \sim 0.05\,\text{K·m²/W}$ |
| 铁芯-机壳 | 压入接触热阻 | $R_{th} \approx 10^{-4} \sim 10^{-3}\,\text{K·m²/W}$ |
| 磁石-转子铁芯 | 粘合剂层 | $R_{th} \approx 0.005 \sim 0.02\,\text{K·m²/W}$ |
| 气隙 | 空气低热导率 + 对流 | $h \approx 50 \sim 300\,\text{W/(m²·K)}$ |
永磁体的失磁风险
刚才不可逆失磁的话，具体是温度多少就危险啊？
钕铁硼磁石的等级不同，危险温度也不同。总结一下：
| 磁石等级 | $B_r$ (20°C) | 最大使用温度 | $\alpha_{B_r}$ [%/°C] |
|---|---|---|---|
| N52（高$B_r$、低耐温） | 1.43 T | 80°C | -0.12 |
| N42SH（中等$B_r$、高耐温） | 1.30 T | 150°C | -0.11 |
| N38UH（低$B_r$、超高耐温） | 1.23 T | 180°C | -0.10 |
| SmCo（钐钴） | 1.05 T | 300°C | -0.035 |
不可逆失磁是否发生，不只看温度，还看反磁场强度。d轴电流大的弱磁运行中高温，B-H曲线上的工作点会掉到拐点以下，导致不可逆失磁。所以"最大电流 × 最高温度"这个最坏组合来检验失磁裕度是铁律。
弱磁运行中危险啊。高速公路全力加速时磁石坏了…想想就恐怖。
绝缘等级和许可温度
磁石以外，绕组绝缘也有温度限制吧？
对。绕线的绝缘被膜有IEC 60085定义的耐热等级，超过就绝缘寿命急剧下降。
| 绝缘等级 | 最高许可温度 | 代表绝缘材料 | 用途 |
|---|---|---|---|
| B级 | 130°C | 聚酯系 | 产业用通用电机 |
| F级 | 155°C | 环氧系 | 一般EV电机 |
| H级 | 180°C | 硅酮系 | 高功率EV、航空 |
| N级 | 200°C | 聚酰亚胺系 | 特殊高温环境 |
实务中常见的是"F级绝缘、B级温升"——也就是155°C耐热的材料，但设计上实际温升只控制在130°C以内。这样就有25°C的安全裕度。
Steinmetz的经验式 — 100多年还在用
1890年代，Charles Steinmetz从大量实验数据中总结出铁损与磁通密度、频率的关系式。当时没有电脑，全是手摇计算和实验反复。惊人的是，这个"Steinmetz式"到了现在，JMAG、Maxwell、Motor-CAD这些最前沿的电机设计工具里还在用。理论跟不上磁区行为的复杂性，但实验数据库的力量跨越了130多年还活着。前人的执着真是让人敬佩。
电机热管理的数值计算方法
LPTN（集中参数热路网）
电机温度计算都用有限元吗？有限元耗时啊…
好眼光。其实电机热分析有两大方法：
- LPTN（集中参数热路网） — 集中定数的热路网模型
- 有限元法（FEM） — 3D分布定数模型
设计初期的高速参数化研究用LPTN，详细设计的温度分布确认用FEM。业界就这样分工的。
LPTN…电路一样的东西？
完全是这样！电路的类比用在热问题上：
| 电气回路 | 热回路 | 单位 |
|---|---|---|
| 电压 $V$ | 温度 $T$ | °C 或 K |
| 电流 $I$ | 热流 $\dot{Q}$ | W |
| 电阻 $R$ | 热阻 $R_{th}$ | K/W |
| 电容 $C$ | 热容 $C_{th}$ | J/K |
| 电流源 | 发热源（损失） | W |
电机的各部件（绕组、齿、轭、磁石、轴、机壳等）做成节点，部件间用热阻连接。一个典型电机约20～50个节点的网络。
$[C_{th}]$: 热容矩阵，$[G_{th}]$: 热导纳矩阵（$G = 1/R_{th}$），$\{P_{loss}\}$: 各节点发热量向量
速度怎么样？
LPTN一个运行点毫秒级解出来。效率特性图全域（比如100×100 = 10,000个点）数十秒就算完。有限元同样的计算要耗数日。所以Motor-CAD这样的工具是LPTN基础的。
有限元耦合分析
有限元连成是怎么做的？
有限元的电磁-热耦合通常这样流程：
- 电磁有限元做1转周期（或数个电气周期）的过渡分析 → 各单元的时间平均损失密度算出
- 损失密度图映射到热有限元的源项
- 热有限元计算温度分布
- 从温度分布更新材料物性（铜电阻率、磁石$B_r$、铁损系数）
- 用更新物性重做电磁有限元
- 温度收敛为止（通常3～5次迭代就收敛）重复步骤2～5
电磁和热别的网格用吗？
很好的问题。电磁和热的要求很不一样：
- 电磁网格: 气隙要超细（0.1mm以下要素）才能算出磁通密度准。2D断面数万要素。
- 热网格: 气隙用对流边界条件代替，只有固体部分网格。3D数十万～数百万要素。
的同网格使、"映射"処理必要。JMAGMaxwell的機能内蔵。
冷却回路建模
冷却系怎么建模？水冷套管、油冷喷射之类的。
冷却回路建模分3个等级，精度和复杂度逐级升高：
等级1：固定热传递系数
冷却面设定一个常数 $h$。最简单但精度低。
水冷套管大约 $h \approx 3000 \sim 5000\,\text{W/(m²·K)}$。
等级2：1D冷媒流路模型
冷媒温度的上游→下游变化作为1D流路计算。从Nusselt数相关式动态算出 $h$。
其中 $P_{wet}$ 是润湿周长，$\dot{m}$ 是质量流量。Motor-CAD采用这种方式。
等级3：3D CFD耦合
冷却流路也用Fluent/Star-CCM+做3D网格流体分析。精度最高，但计算代价也最大。套管形状优化或油冷喷射设计时用。
实务上用哪个多？
设计阶段圧倒的多数等级2。等级3冷却系的几何固最終检証段阶使事多。EV业界最近増的ATF（自动变速箱油）油冷、线圈端直接喷射。这等级2精度出、CFD必要场合多。
耦合策略：弱耦合 vs 强耦合
弱耦合和强耦合什么区别？
| 项目 | 弱耦合（单向/序列） | 强耦合（双向/完全耦合） |
|---|---|---|
| 数据交换 | 电磁→热单向，或少数反复 | 各時间步双向交换 |
| 计算代价 | 低（各求解器1～5次运行） | 高（每步都收敛反复） |
| 精度 | 定常緩过渡充分 | 急激过渡现象必要 |
| 工具例 | JMAG→Fluent、Maxwell→Icepak | COMSOL（内蔵耦合） |
| 适用场景 | 定常额定运行，效率特性图 | 锁定转子，过负荷瞬间 |
实务的90%以上弱耦合充分。强耦合本当必要的、回転（拘束通电）数秒温度数百度変极端案例。
油冷 vs 水冷 — 实务迷热设计的分岐点
EV用的冷却方式议论的「油冷（ATF喷射）vs 水冷（）」的选択。水冷热抵抗小制御反面、直接触热经路长。油冷直接ATF喷射的局所冷却强力、油的撹拌损增。実务连続出力重视水冷、出力重视油冷选案例多。最近的EV両方并用冷却。
电机热管理的实务应用
分析流程
老师，实际做电机电磁-热耦合分析时，从哪里开始？
典型的工作流程是这样的：
- 电机规格整理: 额定/最大电流、转速范围、冷却方式、绝缘等级
- 电磁有限元（2D过渡分析）: 代表运行点（额定/最大转矩/最高转速）下的损失计算
- LPTN或3D热有限元: 损失映射后的温度分布获取
- 材料物性的温度修正: 铜电阻、磁石$B_r$更新后重做电磁有限元
- 收敛确认: 温度变化1°C以内为止反复
- 失磁检验: 最高温度时磁石工作点是否超过拐点
- 工况过渡分析: WLTP、US06等行驶工况的温度历史确认
网格划分的要点
电机热分析的网格有什么特别注意的？
电机特有的网格注意事项总结一下：
- 绕组建模: 数千条铜线逐个网格划分不现实。用等效截面建模，占积率反映在等效物性中。
- 硅钢叠层: 逐片建模不可行。沿叠层方向设各向异性热导率的均质体近似。
- 线圈端: 2D电磁分析中不存在，但铜损的20～30%在线圈端发生。3D热分析中必须建模。
- 气隙: 电磁网格中超微细（0.1mm以下要素），热网格中用对流边界条件代替
- 接触热阻面: 铁芯-机壳界面、磁石-转子界面设接触热阻（TCR）
线圈端漏掉的话，绕组温度会低估2～3成吧？
完全是。线圈端是热上最不利的场所。槽内绕组靠近铁芯，热容易散，但线圈端凸出到端部空间的空气中，冷却受限。最高温度几乎100%在线圈端发生。这里疏漏了，设计就崩了。
边界条件和热传递系数的设置
热传递系数 $h$ 的值怎么定？有标准值吗？
参考值是有的，但实际上流速、温度、形状依赖性很强，从Nusselt数相关式算出来才正确。
| 冷却面 | $h$ 范围 [W/(m²·K)] | 说明 |
|---|---|---|
| 自然对流（外表面） | 5～25 | 密闭型电机外表面 |
| 强制空冷（风扇） | 30～100 | 全闭外扇式（TEFC） |
| 水冷套管（层流） | 500～2,000 | 流量不足时 |
| 水冷套管（湍流） | 3,000～8,000 | 流量充分时 |
| ATF油冷（喷射） | 500～3,000 | 喷射速度、角度相关 |
| 气隙（转子表面） | 50～300 | 转速相关。Taylor-Couette流 |
特气隙的热传达初心者苦労。回转的间的Taylor-Couette流、Taylor数 $Ta$ 遷移决：
$\omega$: 角速度、$r_m$: 平均半径、$\delta$: 气隙长度、$\nu$: 動粘性係数。$Ta > 1700$ Taylor涡発生、热传达急増。
验证和有效性确认
分析结果正、確認？
電機热分析的验证3段階行：
- 保全: 入力損失的合計 = 冷却系的放热量 + 蓄热量（過渡的場合）。合模型。
- 感度分析: 不確実性的大参数（接触热抵抗、绕组等効热导率） ±30% 振温度的影響確認。
- 実験的照合: 热电対Pt100绕組温度、赤外線表面温度、冷媒的入出口温度差（法）比较。許容差 ±10°C 一般的目标。
±10°C结構厳。実際合的？
接触热抵抗绕组等効物性実験、±5°C以内合可能。逆言、的参数「教科書的値」的使20～30°C。実験数据分析精度的鍵。
电机热管理的软件比较
主要工具比较
电机电磁-热耦合的分析工具、选択肢？
大分「専用工具」「汎用多物理场工具」。
| 工具 | 开发方 | 热分析方法 | 电磁耦合 | 优势 |
|---|---|---|---|---|
| Motor-CAD | Ansys（旧MDL） | LPTN + CFD | Maxwell联动 | 设计初期的高速探索 |
| JMAG-Designer | JSOL | 3D FEM热 | 内蔵耦合 | 電磁-熱-结构的一体分析 |
| Maxwell + Icepak | Ansys | 3D FEM/CFD | Workbench联动 | 詳細3D CFD冷却分析 |
| COMSOL | COMSOL AB | 3D FEM | 完全内蔵耦合 | 任意的物理柔軟連成 |
| Flux + Portunus | Altair | LPTN + FEM | 系统連成 | 制御系的統合 |
| MotorWizard + FEMM | OSS + 商用 | 外部連携 | 手動 | 低 |
Motor-CAD最近聞、FEM的精度大丈夫？
疑問。Motor-CAD的LPTN50年以上的学術研究基热回路網、適切定常温度±5°C以内的精度出。局所的温度分布（例内的绕组温度的不均一性）見。FEM的出番。
实务中的使分：
- 设计探索（100案以上的） → Motor-CAD（LPTN）
- 詳細设计的验证（候補2～3案的最終確認） → JMAG / Maxwell+Icepak（FEM）
- 冷却回路的优化 → Star-CCM+ / Fluent（3D CFD）
选择指南
予算限中小企業？
予算目的応段階的揃的：
- 予算: FEMM（2D電磁、無料）+ Excel/Python手動LPTN。学術研究十分。
- 年間100～300万円: Motor-CAD単体。设计探索十分。
- 年間500～1000万円: JMAG + 3D熱分析。日本的的標準的構成。
- 年間1000万円以上: Ansys一式（Maxwell + Motor-CAD + Fluent + Icepak）。大手OEM的統合環境。
Motor-CAD業界標準理由
Motor-CAD（英Motor Design Ltd社、2022年Ansys买収）自動車航空业界广最大的理由、電磁分析热分析同工具内耦合点。従来「電磁FEA損失分布出、別的热FEM入力」2必要、数据変换的模型的不整合頻発。Motor-CADLPTN的热分析電磁分析自動、设计反復大幅短縮。50程度的LPTN複雑流体冷却的詳細表現——的限界理解。
电机热管理的故障排查
常见错误和对策
先生電機热分析「、…」？
何度（笑）。故障症状別总结。
| 现象 | 原因 | 对策 |
|---|---|---|
| 绕组温度比实测高50°C以上 | 绕组等效热导率设定过低（浸渍不足仮定） | 占积率浸渍率実態合等効物性再计算 |
| 磁石温度比实测低很多 | 磁石的渦電流損無視 | 磁石的分割導電率考慮渦電流損计算追加 |
| 连成反复振动不收敛 | 銅損的温度依存性强、正的发散 | 缓和係数（$\alpha = 0.3 \sim 0.7$）更新量制限 |
| 冷媒出口温度不符 | 冷媒流量设定実機異 / 無視 | 実測流量数据使用。流確認 |
| 机壳表面温度偏高 | 外表面的自然对流 $h$ 低 | 放热的効果正模型化。実測 $h$ 同定 |
| 定常分析非现实的高温 | 热的逃道（断热境界的设定） | 全的（入力電力 = 出力 + 全損失 = 全放热）確認 |
缓和係数意味？
連成反復温度 $T$ 更新、前回的値今回的计算値的加重平均取：
$\alpha = 1$ 计算値的使（振動）。$\alpha = 0.5$ 半分更新（安定収束遅）。経験的 $\alpha = 0.5 \sim 0.7$ 良。
调试检查清单
「分析合！」的、系统的方法教！
的检查清单上順、大抵的問題特定：
- 保全: 全損失 [W] ≒ 冷媒的吸热量 $(\dot{m} c_p \Delta T_{coolant})$ + 外表面放热量 ？
- 损失的妥当性: 電磁FEM的損失値効率整合？ 効率90%的入力10kW全損失約1kW。
- 材料物性的異方性: 硅鋼板的叠層方向的热伝導率面内同値？（）
- 接触热抵抗: 铁芯-機壳界面的 $R_{th}$ 设定？ （完全接触）的非現実的低温度。
- 線圈端: 3D模型線圈端模型化？ 省略最高温度見落。
- 冷媒温度的上流-下流勾配: 冷媒温度一定仮定？ 実際下流行温度上。
- 単位系: 電磁FEM的损失密度 [W/m³] 热FEM的入力単位整合？
的、実践的。保存！ 電機的热分析、電磁熱流体的3的知識全部必要大変、。
的通。電磁-熱耦合CAE的中最「多物理场」領域的一。大変、EV時代引張。頑張！