📚 课程目录第4章电池产热速率模型

4. 电池产热速率模型

电池产热速率模型是热模型建模的“源项”。无论后续用一维热阻网络、三维有限元还是降阶模型，热源算错，温度场就失去意义。本节围绕三类常用产热模型展开：Bernardi产热模型、简化产热模型和电化学-热耦合产热模型，并说明不同工况下产热特性的差异与判读方法。

一、Bernardi产热模型：从能量守恒到可计算形式

Bernardi模型的核心思想是把电池视为一个均匀产热的控制体，用能量守恒把产热速率拆成不可逆热和可逆熵变热两部分。其常用形式为：

q = I·(U_ocv − U) + I·T·(∂U_ocv/∂T)

其中 q 为产热速率（W），I 为电流（A，规定放电为正、充电为负，全文统一此约定），U_ocv 为开路电压（V），U 为端电压（V），T 为温度（K），∂U_ocv/∂T 为熵变系数（V/K）。第一项 I·(U_ocv − U) 是不可逆热，来自欧姆内阻、极化内阻和传质阻力；第二项 I·T·(∂U_ocv/∂T) 是可逆热，来自锂嵌入/脱出过程的熵变，放电时可能吸热，充电时可能放热。

实际使用时，U_ocv − U 可由内阻近似：U_ocv − U ≈ I·R_total，因此 q ≈ I²·R_total + I·T·(∂U_ocv/∂T)。R_total 包含欧姆内阻 R_ohm、电荷转移内阻 R_ct 和扩散内阻 R_diff，通常随 SOC 和温度变化。熵变系数 ∂U_ocv/∂T 需通过不同温度下的 OCV 曲线测得，典型值在 −0.5～0.5 mV/K 之间，符号随 SOC 变化。

教学示例：某 3.7 V/50 Ah 三元锂电池，在 25 ℃、1C 放电时，R_total 取 1.2 mΩ，∂U_ocv/∂T 取 −0.1 mV/K。则不可逆热约为 50²×0.0012 = 3.0 W，可逆热约为 50×298×(−0.0001) = −1.49 W，总产热约 1.51 W。该数值为合理假设示例，用于说明量级与符号关系，并非实测结果。

二、简化产热模型：工程快速估算的取舍

当缺少熵变数据或只关心温升趋势时，可忽略可逆热，得到简化产热模型：

q ≈ I²·R_total

进一步，若把 R_total 视为常数，则产热速率与电流平方成正比，这是分析倍率工况最常用的近似。简化模型适合：初步热设计、散热方案对比、BMS 产热估算。但它会低估充电初期的吸热效应，也会忽略低温下内阻急剧上升带来的产热增加。

为兼顾精度与计算量，工程上常把 R_total 写成 SOC 和温度的函数，用二维查表实现：

R_total = interp2(SOC_grid, T_grid, R_table, SOC, T)
q = I² · R_total + I · T · dUocv_dT(SOC)

这样既保留 Bernardi 的两项结构，又避免在线求解电化学方程，是热管理控制器中较常见的折中方案。

三、电化学-热耦合产热模型：从“总热”到“分层热”

Bernardi 模型给出的是电池整体产热，无法回答“热在哪里产生”。电化学-热耦合模型把产热分解到正极、负极、隔膜和电解液等区域，并与 P2D（准二维）或单粒子模型耦合。其产热项通常包括：

欧姆热：j_ohm = σ_eff·(∂φ_s/∂x)² + κ_eff·(∂φ_e/∂x)²，分布在固相和液相；
反应热：j_rxn = a_s·i_loc·(φ_s − φ_e − U)，集中在电极颗粒表面；
极化热：由浓差极化、电荷转移极化贡献；
可逆熵热：j_rev = a_s·i_loc·T·(∂U/∂T)。

耦合方式为：电化学模型输出各区域产热速率 → 输入热模型 → 热模型更新温度 → 温度反馈修正电化学参数（如扩散系数、反应速率常数）。该模型精度高，但计算量大，常用于电芯设计、极耳优化和热失控前期分析。

四、不同工况下的产热特性分析

产热速率不是常数，它随倍率、温度、SOC 和充放电方向显著变化。理解这些趋势，才能正确设置热模型边界和验证工况。

倍率影响：不可逆热与 I² 成正比，因此 2C 放电的欧姆热约为 1C 的 4 倍。可逆热与 I 成正比，倍率越高，不可逆热占比越大，总产热越接近 I²R 规律。小倍率下可逆热甚至可能使总产热为负（吸热），表现为放电初期温度短暂下降。

温度影响：低温下电解液黏度上升、锂离子扩散变慢，R_diff 和 R_ct 显著增大，导致产热增加；高温下内阻下降，但副反应和 SEI 分解可能带来额外产热。因此低温充电是热管理重点工况。

SOC 影响：石墨负极在低 SOC 区熵变系数可能为正，高 SOC 区为负；正极材料（如 NCM）的熵变曲线形状不同。两者叠加后，全电池 ∂U_ocv/∂T 随 SOC 呈非单调变化，导致可逆热在充放电过程中可能变号。

充放电方向：放电时 I 为正，充电时 I 为负。若 ∂U_ocv/∂T 为负，放电时 I·T·(∂U_ocv/∂T) 为负（吸热），充电时该项为正（放热）。这就是为什么充电末期温升往往比放电更明显。

注意：Bernardi 模型假设电池内部温度均匀、产热均匀，这在电芯尺寸较大或倍率较高时误差明显。此时应改用分层产热模型或三维电化学-热耦合模型，至少把极耳和卷芯中心作为独立节点处理。

五、模型选择与验证建议

实际项目中可按以下思路选择：若只需估算整包温升趋势，用简化模型 q = I²·R_total(SOC,T) 即可；若需分析充放电温升差异和低温加热策略，用 Bernardi 模型并保留熵热项；若需优化极耳位置、隔膜热稳定性或热失控传播，则必须用电化学-热耦合模型。

验证时，建议至少设置三组工况：1C 放电、2C 放电、−10 ℃ 充电。用实测表面温度和端电压反推产热速率，与模型计算值对比。若偏差超过 15%，优先检查 R_total 的 SOC/温度依赖是否准确，以及熵变系数符号是否与约定一致。产热速率模型是热模型的输入，输入不准，后续温度场标定就只是“拟合误差”，而不是真正的参数辨识。