TMS BOOK · ACADEMY 课题
动力电池产热机理与产热速率建模
电池热管理的一切设计都始于一个数：这颗电芯此刻发多少热。本课把产热从「I²R」这个粗糙印象拆回 Bernardi 方程，讲清焦耳热与熵热各自的脾气、参数怎么实测、以及怎么从电芯放大到整包并交给 1D/3D 仿真。
D1-01 · 电池热管理 BTMS / 电池产热与热失控基础
待认领listed
→
众筹中recruiting
→
已开课open · Course 售卖
一课题可并存多讲师版本
课程大纲
学完能做什么
- 用 Bernardi 方程分解出可逆熵热与不可逆极化热，并判断各自量级
- 从 HPPC 脉冲数据算出随温度/SOC/倍率变化的内阻表
- 选出与问题匹配的产热模型颗粒度（集总/ECM/MSMD/电化学）
- 把电芯级产热率换算成模组与包的体积热源，接进 CFD 或 1D 回路
- 用量热试验验证产热模型，判断误差来自参数还是结构
内容大纲
- 电芯里的热从哪来
- 不可逆热：欧姆极化 + 活化极化 + 浓差极化，三者随倍率与温度的占比迁移
- 可逆热：熵变热，随 SOC 变号——充电吸热/放电放热并非恒定
- 副反应热与接触电阻热：正常工况可忽略，滥用与老化末期不可忽略
- 极耳与汇流排：局部电流密度集中，是包内的隐藏热点
- 判据：小倍率下熵热可与焦耳热同量级，大倍率下极化热主导
- Bernardi 方程与产热项拆解
- Q = I(U_oc − U) − I·T·(∂U_oc/∂T) 的物理来源：电能未转成化学能的部分 + 熵变
- 第一项恒为正（耗散），第二项可正可负（熵热符号随 SOC 平台翻转）
- 为什么用 (U_oc − U) 而不是 I·R：前者是实测量，后者依赖内阻定义口径
- 从单体功率 W 到体积产热率 W/m³：卷芯体积 vs 电芯外廓体积的口径差
- 简化到 Q ≈ I²R 的适用边界：中高倍率、非平台区、只求温升量级
- 参数怎么来：内阻、熵热系数与热学物性数据包
- HPPC 脉冲法：不同脉冲宽度得到不同 R，欧姆内阻/极化内阻要分开取
- 内阻的三维表 R(T, SOC, I)：低温段与低 SOC 段是强非线性区，采点要加密
- 熵热系数 ∂U_oc/∂T：恒温静置 + 阶梯变温测 OCV 斜率，量级在 mV/K 以内
- EIS 交流阻抗与 DCIR 的区别：不同频段对应不同物理过程，别混用
- 老化会同时抬高 R 与改变熵热特性——热设计参数应取 EOL 而非 BOL
- 热学参数（c_p 与各向异性 k∥/k⊥、内部等效热阻）的测法、取样适用性与向电芯厂索取时的验收口径，见 D1-07；本课只用其数值、不测它
- 从电芯到模组到包：尺度与模型选型
- 集总参数（单节点）：算包级平均温升与热平衡，秒级出结果，看不到温差
- ECM + 热网络：控制导向与 1D 回路联仿的主力，可跑全工况
- MSMD / NTGK 类分布式模型：需要包内温度场与电流分布时用
- P2D 电化学模型：研究析锂、快充边界时用，整包算不动，用于标定与降阶
- 各向异性导热：卷绕/叠片面内与厚度方向导热系数常差一个量级，必须分别赋值
- 建模落地与验证
- 等温量热仪测总产热：给模型一个不含结构误差的标尺
- ARC 绝热温升反推产热：适合大倍率与近失控边界，注意热惯量因子修正
- 整包温升试验对标：先对总能量，再对温度场分布，顺序反了会用结构误差补参数误差
- 误差归因：总产热对不上 → 参数问题；总量对得上而分布错 → 热阻链或流场问题
- 交付：一张 R(T,SOC,I) 表 + 一条 ∂U_oc/∂T(SOC) 曲线 + 一份验证报告 + 一份电芯热参数数据包（ρ、c_p、k∥/k⊥、内部等效热阻及各自测试状态声明）
关键公式
Q_gen = I(U_oc − U) − I·T·(∂U_oc/∂T)
Bernardi 产热方程：第一项极化耗散热（恒正），第二项可逆熵热（随 SOC 可正可负）
Q_irr = I·(U_oc − U) ≈ I²·R_int
不可逆极化热；仅当内阻口径与工况脉冲宽度匹配时才可用 I²R 近似
Q_rev = −I·T·(∂U_oc/∂T)
可逆熵热；低倍率长时工况（如慢充、静置放电）下不可忽略
R_DC = (U_oc − U(t_pulse))/I
HPPC 直流内阻；t_pulse 取值不同得到的 R 不同，须与目标工况时间尺度一致
ρ·c_p·∂T/∂t = ∇·(k∇T) + q̇_v
电芯内部能量方程；q̇_v 即 Bernardi 产热折算的体积热源，k 需按各向异性赋值
q̇_v = Q_gen / V_cell
体积产热率；交给 CFD/1D 的标准接口量，注意 V_cell 用卷芯体积还是外廓体积要声明
关键概念
Bernardi 方程熵热系数欧姆极化活化极化浓差极化DCIRHPPCOCV体积产热率各向异性导热ECM 等效电路MSMD
推荐工具与标准
GT-SUITE / GT-AutoLion（电-热-老化耦合） Simcenter Amesim 电池库 Fluent MSMD / STAR-CCM+（包内温度场） COMSOL 电池模块 / PyBaMM（P2D 与降阶） Arbin / Neware 充放电台 + 等温量热仪
GB/T 31486-2024（电动汽车用动力蓄电池电性能要求及试验方法；2024-09-29 发布、2025-04-01 实施，推荐性国标）——射程仅到电性能试验口径（容量、能量、功率、直流内阻的测试条件与判定），可用来对齐本课内阻实测的取值口径；产热率、熵热系数与量热标定不在其射程内，本课的产热参数包不能声称「按该标准测得」 IEC 62660-1（电动道路车辆用锂离子电池性能试验） USABC / INL 的 HPPC 测试规程 （本课不再引 SAE J1798：其继任件 J1798/1 只管铅酸与镍氢模块、J1798/2 虽覆盖锂离子但全线是电性能评级——无量热、无产热率测定、不给任何热学判据，给不出本课需要的输入量。本课的对口依据即同列的 USABC/INL HPPC 规程、IEC 62660-1 与 GB/T 31486-2024；熵热系数与产热量本身目前无对口的公开标准方法，等温量热仪与 ARC 反推属实验室方法，按「方法 + 测试状态声明」写，不挂标准号充数。）
工程案例
某纯电 SUV 做 0℃ 环境起充、电芯未预热的 4C 快充温升预测：沿用 25℃ 常温内阻建模，把整条温升曲线连同末段峰值一起算低了一大截。误差其实发生在充电前段的低温区（该段内阻较常温有数倍抬升），但温度是产热对时间的积分，所以最终体现在末段峰值上——这一点不补明起充温度就读不通：若按默认常温起充，4C 下电芯越充越热，用常温内阻反而会高估末段产热。补测低温段内阻后幅值才对上；再加入熵热项，平台区的曲线形状也对上了（该倍率下熵热约为 I²R 的一成到两成，只改形状、改不动那「一大截」的幅值）。结论落在低温快充窗口上：预热到多少温度才允许起 4C、加热与冷却的切换门限怎么定，而不是一路上调冷板换热能力。
动手做
交付物 · 拿电芯 HPPC 数据：① 按 1 s / 10 s 提取 R，画 R(T,SOC) 曲面标出强非线性区；② 用 Bernardi 算 2C 放电段 Q_irr 与 Q_rev 占比；③ 折算体积产热率接进 1D 热网络对标实测。交付：内阻表 + 产热率曲线 + 归因说明。
常见误区
- 只算 I²R 漏掉熵热：小倍率或电压平台附近，熵热可与焦耳热同量级且可能为负（吸热），漏掉会把慢充温升算高、把某些 SOC 段的走势算反
- 用常温内阻算低温工况：内阻对温度强非线性，低温段抬升可达数倍，低温快充产热会被严重低估——内阻必须建成温度的函数
- 用稳态/长脉冲内阻算瞬态工况：HPPC 取 10 s 与取 1 s 的 R 不是一个数，脉冲宽度须与目标工况时间尺度对齐
- 把电芯当各向同性材料：卷绕/叠片电芯面内与垂直方向导热系数常差一个量级，各向同性假设会把电芯内部温差算平、把极耳热点抹掉
- 用 1C 标定的模型外推 4C：极化项高度非线性，低倍率标定的参数外推到超充区必失真，须在目标倍率区间内标定
- 把本课的三类极化直接搬去燃料电池：这里的欧姆/活化/浓差极化是锂离子电池口径（载流子是 Li⁺ 在电极与电解液之间嵌入/脱出），PEMFC 同名的三类极化在载流子、反应物供给方式与失效后果上完全不同，转 FCEV 方向不可直接代入——对照见 O1-08
相关课题