4、热模型基础：电池产热机理与建模方法

做电池仿真这么多年，我越来越觉得热模型是个「绕不过去的坎」。你想想看，电池性能好不好，寿命长不长，安全不安全，温度几乎是个决定性因素。今天我们就来聊聊热模型的基础——电池到底怎么发热的，以及我们怎么把这些热算清楚。

4.1 电池产热的三大来源

电池发热，说白了就三个原因：焦耳热、反应热、极化热。我刚开始做热仿真时，总以为发热就是电流大了温度高，后来才发现事情没那么简单。

4.1.1 焦耳热（欧姆热）

这个最好理解。电流流过电池内部的电阻，就会产生热量。公式很简单：

Q_ohm = I² × R

其中 I 是电流，R 是电池内阻。内阻包括电子电阻（集流体、极耳、活性材料）和离子电阻（电解液、SEI膜）。

4.1.2 反应热（熵变热）

这个稍微绕一点。电池充放电时，正负极材料发生电化学反应，反应本身会吸热或放热。公式是：

Q_rxn = T × ΔS × (I / nF)

T 是温度，ΔS 是反应熵变，I 是电流，n 是电子转移数，F 是法拉第常数。

嗯，这里要注意：反应热可正可负。放电时通常是放热（正），但有些材料在特定 SOC 区间会吸热。我记得有次做 LFP 电池仿真，客户说「为什么放电到 50% SOC 时温度反而降了？」其实就是反应热在作怪。

4.1.3 极化热

极化热是「过电位」带来的热量。电池实际工作电压和开路电压之间的差值，就是过电位。这部分能量最终以热的形式耗散掉：

Q_pol = I × (V_oc - V_term)

极化又分活化极化和浓差极化。活化极化跟反应动力学有关，浓差极化跟离子扩散有关。大电流放电时，浓差极化会特别明显。

4.2 集总参数热模型

说白了，集总参数模型就是把整个电池当成一个「点」，不考虑内部温度分布。这个模型简单、计算快，适合系统级仿真。

核心公式就是能量守恒：

m × Cp × dT/dt = Q_gen - Q_diss

m 是电池质量，Cp 是比热容，Q_gen 是总产热，Q_diss 是散热。

散热通常用牛顿冷却定律：

Q_diss = h × A × (T - T_amb)

h 是对流换热系数，A 是散热面积，T_amb 是环境温度。

我个人的习惯是：做 BMS 热管理策略时，先用集总参数模型跑一遍，看看整体趋势。等方案定了，再用三维模型精算。

4.3 三维热模型

三维热模型就复杂多了。它要考虑电池内部各个位置的温度差异，比如电芯中心温度可能比表面高好几度。

控制方程是三维热传导方程：

ρ × Cp × ∂T/∂t = ∂/∂x(kx × ∂T/∂x) + ∂/∂y(ky × ∂T/∂y) + ∂/∂z(kz × ∂T/∂z) + Q_gen

ρ 是密度，Cp 是比热容，kx、ky、kz 是各向导热系数。

这里有个坑：电池的导热系数是各向异性的。沿极片方向导热好（约 20-30 W/m·K），垂直方向导热差（约 0.5-1 W/m·K）。我曾经见过有人用各向同性模型算方形电池，结果中心温度低估了 5°C 以上。

4.3.1 三维模型的边界条件

常见的边界条件有三种：

• 对流边界： 电池表面与空气或冷却液换热
• 辐射边界： 电池表面向环境辐射热量
• 绝热边界： 对称面或隔热面

实际项目中，对流边界最常见。但要注意：对流换热系数 h 不是常数，它跟风速、流体性质、表面状态都有关系。

4.3.2 三维模型的网格划分

网格划分是个技术活。网格太粗，算不准；网格太细，算不动。

我的建议是：

• 电芯内部：至少 3-5 层网格
• 极耳区域：加密处理
• 表面边界层：加密处理
• 总网格数：控制在 10 万-50 万之间

4.4 知识体系总览

下面这张图是我自己整理的，把热模型的知识结构串起来了。你可以把它当成一个「地图」，学完这章后对照着看看自己掌握了多少。

4.5 模型选择建议

到底用集总参数还是三维模型？我个人的经验是：

场景	推荐模型	理由
系统级仿真（整车热管理）	集总参数	计算快，参数少，够用
电芯内部温度分布	三维模型	需要捕捉热点
快充工况分析	三维模型	大电流下温差大
前期参数敏感性分析	集总参数	快速迭代
热失控仿真	三维模型	需要精确温度场

好了，热模型的基础就聊到这儿。记住：产热机理是根本，模型选择看需求。下次我们聊聊怎么用 Python 实现一个简单的集总参数热模型，到时候我会把代码一行行拆开讲。

---

公众号:蓝海资料掘金营，微信deep3321