J6-03 电机电磁—热耦合仿真
课程代码 J6-03 · 板块 J 仿真分析 CAE / J6 多物理场与电子热 时长 约 4.0 小时(6 讲 + 1 次电磁—热耦合仿真实操) 适合对象 电驱/电机方向的仿真工程师(必修);参与热仿真环节的电机设计工程师 前置 E1-04 Hairpin 扁线电机的散热特点 Hairpin 扁线电机的散热特点(交流损耗认知);E1-05 电机温升限值、绝缘等级与降额 电机温升限值、绝缘等级与降额;A1-06 热阻网络法与集总参数建模入门 热阻网络法与集总参数建模入门(节点划分与 R/C 反标定的方法论宿主) 本讲义定位 讲师授课蓝本 / 学员自学讲义
按任务导航
| 你手上的事 | 直接看 | 还要配 |
|---|---|---|
| 接手一份别人交来的电机温升仿真结果,要判它可不可信 | 引言:三句质询 · 第 2 讲 2.5(映射守恒核对) | 第 1 讲 1.5(铁耗系数用的是哪套标定条件)· 第 5 讲 5.2(反标定顺序对不对) |
| 定这次要不要做双向耦合 | 第 3 讲 3.2 三档判据 + 3.3 耦合频度 | 第 1 讲 1.4(三类损耗的温度依赖方向不一致,这是双向耦合值不值得做的根) |
| 把电磁 FEA 的损耗结果交给做热的人 | 第 2 讲:分区输出 → 栅格扫描 → 量纲换算 → 守恒核对 | 第 4 讲 4.2(节点清单:热源到底该注到哪个节点) |
| 搭一张电机 LPTN,要给支路取参 | 第 4 讲 4.3 三类关键支路参数与它们各自的属主课 | 第 5 讲 5.1~5.3(哪些参数只能靠温升试验反标定) |
| 判短时过载能撑多久、退磁裕度够不够 | 第 6 讲 6.2(绝缘侧)· 6.4(退磁侧,看拐点⛔ 不看剩磁) | 第 3 讲 3.4(这类判断该用 LPTN 还是 3D) |
| 把离线模型降阶交给上车的在线热模型 | 第 5 讲 5.5 交接口径与冻结清单 | E1-07 电机不可测温度的在线估计与热保护标定:绕组热点、磁钢温度与退磁保护(那门课是接收端)· J7-02 模型降阶(ROM)与控制导向建模(降阶方法本体) |
⛔ 本课不覆盖 · 栏一(不讲哪些题目):电磁 FEA 本身的建模、网格、时步与磁路设计——本课只消费它的损耗输出;多物理场耦合的通用架构与跨工具协同 ⇒ J6-02 电—热—流多物理场耦合仿真;共轭传热(CHT)的求解设置与网格判据 ⇒ J2-04 共轭传热(CHT)仿真方法;热阻网络法本体、节点划分方法论与 R/C 反标定的方法学 ⇒ A1-06 热阻网络法与集总参数建模入门(本课只做电机这一个对象上的落地);机壳水道与油冷结构的换热强化设计 ⇒ E1-01 电机水冷(机壳水道)设计与换热强化、E1-03 油水混合冷却与转子冷却设计;气隙对流与端腔换热的关联式、风摩损的机理与标度 ⇒ E1-06 干式电机的转子与端部热路:气隙对流、端部空腔与风摩损;交流损耗的机理、逐层阶梯与 Dowell 一族公式 ⇒ E1-04 Hairpin 扁线电机的散热特点;绝缘热分级与温升限值的等级号与限值 ⇒ E1-05 电机温升限值、绝缘等级与降额;磁体牌号、耐温档与磁导系数的材料侧 ⇒ E1-03 油水混合冷却与转子冷却设计、F5-02 轮毂电机散热与簧下热管理;上车实时的温度估计、退磁保护与降额执行 ⇒ E1-07 电机不可测温度的在线估计与热保护标定:绕组热点、磁钢温度与退磁保护;模型降阶与定点化的方法本体 ⇒ J7-02 模型降阶(ROM)与控制导向建模、J7-04 整车热管理数字孪生与实时模型;仿真—试验相关性的验收判据与合格线 ⇒ J7-01 仿真—试验相关性验证(Correlation);台架台位、埋点与试验规程 ⇒ L1-08 电驱总成台架热性能试验:多温限同步判读与峰值包络标定。
⛔ 栏二(真正不给的数),三类,每样带去向:① 未取得原文的标准条款与限值——任何绝缘等级对应的极限热点温度、任何温升试验的判读条款与考核点,本课只写标准号、它管什么与去向(E1-05 电机温升限值、绝缘等级与降额 与 IEC 60085、IEC 60034-1 / GB/T 755;⚠ 版次以现行目录为准,引用前须核);② 本项目的验收阈值——过载时长限值、退磁裕度合格线、相关性合格线、任何降额门限;③ 某一台电机、某一个牌号的实际参数值——损耗 map 的具体数值、接触热阻与槽内等效导热的本平台取值、气隙换热系数的本机取值、磁体的最高工作温度与退磁曲线拐点,这些由你自家的磁路、牌号与温升试验给出。
⚠ 与栏二相对:平台相关量本课照给典型参考值(区间两端 + 常见值锚点 + 免责句 + 拿到它的路),⛔ 不在「不给的数」之列——但一律⛔ 不是验收依据,你的平台须实测确认。
引言:温度不是热仿真一个人的事
一份电机温升仿真报告放在你面前。绕组热点 T_hotspot 的曲线画得很干净,磁体温度 T_m 在过载 30 秒末尾停在一个「还有裕度」的位置,云图配色专业,网格量写着七位数。会上没人质疑它。
而它可能整条都是错的,原因跟热仿真本身没有任何关系:喂进去的那份损耗,是在一个跟这次工况不一样的温度上算出来的。
这就是这门课存在的理由。电机温升是损耗与散热两侧的博弈,而这两侧分属两个学科、常常分属两个人、有时候还分属两个部门:
- 只做热仿真的人,拿到的是一份「损耗数」。损耗从哪来、分布在哪、为什么是这个形状——他猜不出来,也无从判断这份数对不对。
- 只做电磁仿真的人,算损耗时必须先给定一个温度(铜的电阻率要温度、磁体的剩磁要温度)。那个温度是他假设的,而真实温度恰恰是热仿真的输出。
⇒ 两侧各自都能跑出「看起来合理」的结果,而把它们接起来的那条缝,没有任何一个求解器会替你检查。
三句质询:拿到任何一份电机温升结果,先问这三句
① 这份损耗是在什么温度下算出来的?——是恒定 20 ℃ 的默认值,还是某个假设的热态温度,还是与热模型迭代收敛过的?答不出来,后面的温度曲线就是一条自洽但没有根的曲线。
② 损耗从电磁模型搬到热模型时,总量对上了吗?——电磁 FEA 给的是损耗密度(W/m³),热网络节点要的是功率(W)。中间隔着一次按单元体积的积分与一次分区映射。逐区域核过总量没有?偏差多少?答不出来,那份热量可能已经悄悄少了一块,或者多了一块。
③ 热点是节点温度还是平均温度?——退磁看的是磁体最热那一小块,绝缘看的是绕组最热那一点。拿一个把整块磁体抹平的节点温度去比退磁判据,风险会被系统性漏判。
三句答不齐的结果,⛔ 不许拿去定型。这不是流程洁癖:这三句问的分别是损耗侧的正确性、接口的守恒性、判据侧的口径,它们各自都能让最终结论反号,而三者中没有一个会让求解器报错。
本课的符号与口径,一次钉死
全课通用,后面不再打断推导。⚠ 与本课程体系已有课程共享的量,这里逐条写明是沿用还是本课首定——这一层没有任何工具替你检查,读者互代出错是系统性的。
| 符号 | 量与单位 | 口径与出处 |
|---|---|---|
| P_cu,dc / P_cu,ac | 直流铜耗 / 交流铜耗,W | 全绕组总量。P_cu,ac = P_cu,dc·k_ac |
| k_ac | 交流因子,无量纲 | 沿用 E1-04 Hairpin 扁线电机的散热特点,但本课用的是全绕组等效值;⚠ E1-04 Hairpin 扁线电机的散热特点 另有逐层的 k_ac,j,两者不是同一个数(顶层可比平均高数倍),⛔ 不得互代 |
| P_fe | 铁耗,W | 含磁滞与涡流两项;本课凡写 P_fe 均指定子铁心 + 转子铁心的合计,分区时写 P_fe,s / P_fe,r |
| P_pm | 永磁体涡流损耗,W | 本课首定;⛔ 它不属于 P_fe,注入的节点也不同(磁钢节点,⛔ 不是定子轭节点) |
| P_mech | 机械损耗,W | 轴承摩擦 + 风摩;风摩的机理与标度沿用 E1-06 干式电机的转子与端部热路:气隙对流、端部空腔与风摩损 |
| P_loss | 总损耗,W | = P_cu,ac + P_fe + P_pm + P_mech |
| q''' | 损耗密度,W/m³ | 电磁 FEA 的原始输出形式;⛔ 它不能直接当热网络的源项 |
| Q_i | 节点 i 的热源功率,W | = ∫ q''' dV,积分域=映射到该节点的那些单元 |
| ε_map | 映射守恒偏差,无量纲 | 本课首定:ε_map = (ΣQ_i − P_loss)/P_loss,逐区域各算一次 + 全机算一次 |
| C_i / R_ij / T_i | 节点热容 J/K / 节点间热阻 K/W / 节点温度 ℃ | 沿用 A1-06 热阻网络法与集总参数建模入门 与 E1-07 电机不可测温度的在线估计与热保护标定:绕组热点、磁钢温度与退磁保护 的 LPTN 记法 |
| R_th,* | 一切热阻,K/W | 沿用 E1-04 Hairpin 扁线电机的散热特点 的禁令:⛔ 全课不裸写 R。电阻另写 R_s(定子相电阻,Ω) |
| T_hotspot / T_m | 绕组绝缘热点温度 / 磁体温度,℃ | 沿用 E1-07 电机不可测温度的在线估计与热保护标定:绕组热点、磁钢温度与退磁保护 |
| T_coolant,in | 冷却液(油冷机取油)入口温度,℃ | 沿用 E1-07 电机不可测温度的在线估计与热保护标定:绕组热点、磁钢温度与退磁保护:全课热阻链的参考零点取它,⛔ 不取出口温度、⛔ 不取任何 NTC |
| λ_∥ / λ_⊥ | 槽内等效导热,沿槽向 / 垂直槽向,W/(m·K) | 沿用 E1-04 Hairpin 扁线电机的散热特点:两个方向服从两条不同的混合律,⛔ 不许合并成一个 λ |
| h_gap / h_end | 气隙 / 端腔换热系数,W/(m²·K) | 沿用 E1-06 干式电机的转子与端部热路:气隙对流、端部空腔与风摩损;同单位不同对象,⛔ 不得互相代入 |
| α_Br / α_Hcj | 剩磁 / 内禀矫顽力的温度系数,1/K | 沿用 E1-07 电机不可测温度的在线估计与热保护标定:绕组热点、磁钢温度与退磁保护:一律取正数,负号写在式子里,⛔ 不把负号折进系数 |
| B_r(T) / H_cj(T) / B_knee(T) | 剩磁 / 内禀矫顽力 / 退磁曲线拐点磁密 | 退磁判据看 B_knee,⛔ 不看 B_r(第 6 讲 6.4) |
★ 一条跨课的口径提醒,必须现在说清:本课判「哪里是热点」时沿用 E1-04 Hairpin 扁线电机的散热特点 的基准——热点因子 = 归到该点的损耗(W)× 该点到冷却液的热阻(K/W),⛔ 不写成「局部损耗密度 × 局部热阻」。两种写法在量纲上都自洽,但后者会让人用一个体积平均的密度去乘一个面积平均的热阻,得到一个既不是峰值也不是均值的数。
第 1 讲 · 电机损耗全貌与温度依赖
这一讲要交两样东西:一是四类损耗各自的机理与标度——不知道它随转速怎么长,你就不知道该在栅格的哪一片加密;二是它们各自的温度依赖方向——这是判「要不要做双向耦合」的唯一依据,而且它的结论会让不少人意外。
1.1 四类损耗:先把账目分清
① 铜耗,两部分。直流部分 P_cu,dc = 3·I²·R_s,只跟电流有效值与相电阻有关;交流附加来自集肤效应与邻近效应,写成 P_cu,ac = P_cu,dc·k_ac(f, T)。⚠ 这里的 f 是电频率而不是转速——两者差一个极对数 p,f = p·n/60。机理、逐层阶梯与 k_ac 的算法整套归 E1-04 Hairpin 扁线电机的散热特点,本课只做一件事:把 k_ac 当成一个随频率与温度变的量输进来,⛔ 不许把它冻成常数(这是本讲最后要重点说的一条)。
② 铁耗,两项。磁滞损耗随频率一次方长、涡流损耗随频率平方长,合起来就是那条人人会写的分离式(见「关键公式详解」式①)。⚠ 车用驱动电机的实际铁耗不止这两项——Bertotti 的三项式还有一个异常损耗(附加损耗)项,在高速、薄硅钢的场合它不可忽略。本课在式①的适用边界里把这一条单列。
③ 永磁体涡流损耗 P_pm。定子开槽与电流谐波在气隙里造成的磁场脉动,穿进转子表层与磁体内部,在导电的磁体里感生涡流。⚠ 这一项在整机损耗账上通常不大,但它注入的位置最糟糕:磁体在转子上,散热要么跨气隙、要么走轴与转子铁芯,是全机热阻最大的一条路径之一 ⇒ 单位瓦特造成的温升,它比铜耗高得多。⛔ 因此不能按「占比小就并进转子铁耗」处理——第 4 讲 4.2 会把这条做成一条硬规矩。
④ 机械损耗 P_mech:轴承摩擦 + 风摩(windage)。风摩的阻力矩约按 ω² 标度、功率约按 ω³ 标度(机理与关联式归 E1-06 干式电机的转子与端部热路:气隙对流、端部空腔与风摩损,本课不重讲)。量级上:主流 15,000~20,000 rpm 的乘用车驱动电机上,机械损耗通常只占总损耗的百分之一量级(数十到数百 W);要到 80,000~120,000 rpm 级的燃料电池空压机与 eVTOL 电推进电机上,它才翻身成为主项损耗。⚠ 典型区间与量级判断,来源等级=行业惯例,⛔ 不是验收依据;你的平台须按自己的转速域与轴承/密封型式实测确认。
⚠⚠ 量级小 ≠ 可以省。在 LPTN 里,风摩往往是转子/气隙节点唯一的本地源项——把它并进「机械损耗总账」挂到机壳上,转子那一路就没有热源了,转子与磁体温度会被系统性算低。⇒ 按注入位置分账,⛔ 不按占比分账。
1.2 为什么必须按转速—扭矩栅格给,而不能给几个点
一台电机的损耗不是一个数,是转速—扭矩平面上的一个曲面。原因在四类损耗的自变量不同:铜耗主要跟电流(≈ 扭矩)走,铁耗与磁体涡流耗主要跟频率(≈ 转速)走,机械损耗只跟转速走。⇒ 在 (n, T_e) 平面上,它们的等值线族方向完全不同:铜耗的等值线大致横着(沿转速轴),铁耗的大致竖着(沿扭矩轴)。
⚠ 这里的 T_e 是电磁扭矩(N·m),沿用 E1-07 电机不可测温度的在线估计与热保护标定:绕组热点、磁钢温度与退磁保护 的约定:本课温度一律写 T 带下标且下标不为 e,T_e 唯一表示扭矩,⛔ 全课不裸写 T。
后果:只在三五个「典型工况点」上标定损耗,然后线性插值铺满整个平面,会在两处系统性出错——① 弱磁区(高转速、低扭矩),那里铁耗因为深度饱和与谐波而非线性上翘,线性插值一定低估;② 短时峰值区(高扭矩),那里电流大、k_ac 与饱和都在变,而这个区往往恰恰是标定点最稀的地方(因为台架上不好久驻)。这正是本课要防的第一个陷阱:标定点稀的地方,正是热风险最高的地方。
1.3 高速与低速,谁更热?——这个问题问得不对
学员常问「电机是低速大扭矩更热还是高速更热」。这个问题没有跨平台的答案,因为它取决于两件本平台的事:① 冷却结构对哪一路更有效(水套治定子外圈,对转子磁体几乎无能为力);② 工况谱在哪一片驻留得久。
⇒ 正确的提法是:在你自己的损耗 map 上,把每一格的损耗按注入节点拆开,各自乘上该节点到冷却液的热阻,看哪一格的乘积最大。这正是引言里那条沿用 E1-04 Hairpin 扁线电机的散热特点 的热点基准。这样问出来的答案通常是两个热点、两套工况:绕组热点落在低速大扭矩长驻留(铜耗大 + 风冷弱),磁体热点落在高速弱磁(P_pm + 铁耗大 + 跨气隙散热差)。⛔ 只校核其中一套,另一套就是裸奔。
1.4 四类损耗的温度依赖方向:三个方向,不是一个
这是本讲的核心结论,也是第 3 讲耦合架构选型的唯一依据。
| 损耗 | 温度升高时 | 机理 | 强弱 |
|---|---|---|---|
| 铜耗 P_cu | 变大 ↑ | 铜的电阻率是正温度系数 | 强,是主导项 |
| 铁耗 P_fe | 略微变小 ↓ | 硅钢的电阻率随温度升高 ⇒ 涡流损耗下降;磁滞项也弱降 | 弱 |
| 永磁体涡流耗 P_pm | 略微变小 ↓ | 剩磁 B_r 随温度下降 ⇒ 气隙磁场脉动幅值下降;⚠ 磁体电阻率同时上升,也压涡流 | 弱 |
| 机械损耗 P_mech | 基本不变 — | 与温度无直接关系;⚠ 油冷机的甩油/搅油损耗例外,它随油温降黏而下降 | — |
工程量级:铜的电阻温度系数 α_cu 在 20 ℃ 参考点约 0.00393 /K(纯铜;来源等级=材料物性手册值,绕组用铜的实测值随牌号与加工态略有出入)。⇒ 绕组从 20 ℃ 升到 140 ℃(ΔT = 120 K),相电阻 R_s 上升约 47%,直流铜耗在同电流下同比例上升。⚠ 这是一个可复算的物性推论,⛔ 不是任何平台的验收依据;实际电机上电流本身会因控制策略而变,⛔ 不能直接把 47% 当成铜耗的增量。
⇒ 三个方向不一致,这就是「温度高损耗就大」这句话错在哪里。它对铜耗是对的,对铁耗与磁体涡流耗是反的。净效果取决于这一格里铜耗占比多少:低速大扭矩区铜耗占大头 ⇒ 净效果是正反馈,温度与损耗互相推高;高速弱磁区铁耗占大头 ⇒ 净效果反而是微弱的负反馈,温升会自己收敛得比预期快一点。
⚠ 但⛔ 别把这条推成「高速区不用做双向耦合」:高速区真正的风险不在总损耗的净变化上,而在磁体温度上——它走的是退磁那条判据,看的是 T_m 本身,跟总损耗涨没涨没关系。
1.5 铁耗系数的适用边界:这是本课最容易被整段跳过的一节
式①里的 k_h、k_e(或 Bertotti 的三系数)从哪来?手册与商用软件材料库里的那一套,是按正弦交变磁通、单一频率,在爱泼斯坦方圈或环样上拟合出来的。而车用驱动电机的实际工况,在三条上都与这个标定条件不同:
① PWM 供电的高频谐波。逆变器输出的不是正弦电压,磁通里叠着开关频率附近的纹波。这些高频分量幅值虽小,但涡流损耗按 f² 走 ⇒ 它们对铁耗的贡献被频率平方放大。⚠ 偏差随开关频率上升而上升。
② 齿部的旋转磁通。爱泼斯坦方圈上磁通是交变的(沿一个方向来回),而电机齿根与轭部交界处的磁通轨迹是一个椭圆——它在旋转。旋转磁化的损耗与交变磁化不是同一条曲线,尤其在高磁密下差别明显。
③ 弱磁区的深度饱和。局部磁密进入饱和段后,B-H 关系强非线性,按小信号拟合的系数不再成立。
⇒ 三条的方向是一致的:都让实际铁耗高于模型值。这就是「照搬软件默认材料库、按正弦条件算 PWM 电机」会系统性低估铁耗的机理,而低估的后果是温升裕度与退磁裕度双双偏乐观——两个判据一起朝着危险方向错。
修正路径,三条,按代价从低到高:
| 做法 | 怎么做 | 代价 | 治哪一条 |
|---|---|---|---|
| 谐波分解逐次叠加 | 把磁通波形做 FFT,逐次谐波各自代入式①再相加 | 低(后处理即可) | 主要治 ① |
| 非正弦修正(iGSE 一类) | 用改进的通用 Steinmetz 一族方法,直接吃波形的 dB/dt | 中 | 治 ①,部分治 ② |
| 厂商按 PWM 条件实测的损耗曲面 | 向硅钢供应商索取在接近你的开关频率与直流偏置条件下实测的损耗数据 | 高(要谈、要等) | ①②③ 都覆盖得最好 |
⛔⛔ 无论走哪条,最终一律以自家电机的试验标定系数为准,⛔ 不用软件默认材料库值。这一条没有折中:默认值的标定条件与你的工况不重合,而它永远不会报错。
1.6 本讲小结与交付物
本讲交出的是损耗侧的账,三样:
- 四类损耗的分账,并且写明各自该注入哪个节点(⛔ 不按占比合并,按位置分账);
- 各自的温度依赖方向(↑ / ↓ / ↓ / —),它是第 3 讲判耦合架构的输入;
- 铁耗系数的来源声明:这套系数是哪来的、按什么条件标定的、做没做非正弦修正。⚠ 这条声明必须写进交付物,⛔ 不能只留在做的人脑子里——下一个接手的人无从反查。
后面还有 5 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做