2. 动力电池基础
锂离子电池是当前电动汽车与储能系统最核心的电化学储能载体。理解其工作原理、电化学特性、产热机理与热失控概念，是开展热管理系统建模与仿真的前提。热管理工程师若只掌握传热学与流体力学，而不理解电池内部的电化学反应与产热源项，所建立的模型往往缺乏物理一致性，难以准确预测极端工况下的温度分布与热安全边界。本章从电化学基础出发，逐步推导至产热模型与热失控判据，为后续章节的CFD仿真与集总参数模型奠定物理基础。
2.1 锂离子电池工作原理与电化学特性
锂离子电池的本质是依靠锂离子在正负极之间的浓度差驱动嵌入与脱嵌反应，从而实现化学能与电能的相互转换。以典型的石墨负极与三元正极（NCM）体系为例，充电时锂离子从正极晶格中脱出，经电解液迁移穿过隔膜，嵌入负极石墨层间；放电过程则完全反向。这一过程可用如下简化的电极反应表示：
正极：Li(NiCoMn)O₂ ⇌ Li₁₋ₓ(NiCoMn)O₂ + xLi⁺ + xe⁻
负极：C₆ + xLi⁺ + xe⁻ ⇌ LiₓC₆
总反应：Li(NiCoMn)O₂ + C₆ ⇌ Li₁₋ₓ(NiCoMn)O₂ + LiₓC₆从热管理建模的角度，需要重点关注以下电化学特性参数：
- 开路电压（OCV）：OCV是SOC与温度的非线性函数，其温度导数∂OCV/∂T直接决定反应熵热的大小。不同材料体系的熵热系数差异显著，磷酸铁锂在低SOC区间呈现明显的熵热峰，而三元材料相对平缓。
- 内阻组成：欧姆内阻（R_ohm）来自电子导电与离子导电路径，极化内阻（R_ct+R_diff）来自电荷转移与固相/液相扩散。内阻随SOC、温度、老化状态变化，是产热计算的核心输入。
- 容量与倍率特性：可用容量受放电倍率影响（Peukert效应），高倍率下有效容量下降，同时产热率急剧上升。
2.2 产热机理与热源项量化
锂离子电池在充放电过程中的产热可分解为五个物理来源，其中前三项为常规工况下的主要产热项，后两项在滥用条件下才显著贡献。准确量化各热源项是建立热模型的基础。
| 产热类型 | 物理来源 | 数学表达式 | 占比（常规工况） | 关键影响因素 |
|---|---|---|---|---|
| 不可逆焦耳热 | 欧姆内阻消耗电能 | Q₁ = I²·R_ohm | 30%~50% | 电流、温度、SOC |
| 极化热 | 电荷转移与扩散过电位 | Q₂ = I·(η_ct + η_diff) | 20%~40% | 倍率、电解液电导率 |
| 反应熵热 | 锂离子嵌入/脱嵌熵变 | Q₃ = I·T·(∂OCV/∂T) | 10%~30% | SOC、材料体系、温度 |
| 混合热 | 浓度梯度松弛释放 | Q₄ = ∫Δc·(∂μ/∂c)dV | <5% | 倍率变化、静置时间 |
| 分解热 | SEI/电解液/正极分解 | Q₅ = m·H_r·k(T)·exp(-Ea/RT) | 滥用条件下主导 | 温度、材料热稳定性 |
在工程仿真中，通常将前四项合并为可测的等效产热率，通过等温量热仪（IBC）或加速量热仪（ARC）标定。对于集总参数模型，可采用Bernardi方程近似：
Q_gen ≈ I·(OCV - V_terminal) + I·T·(∂OCV/∂T)
= I²·R_total + I·T·(∂OCV/∂T)其中第一项为不可逆热（包含欧姆与极化），第二项为可逆熵热。该公式在中等倍率（≤2C）下误差小于10%，是工程建模的首选简化形式。
2.3 热失控概念与演变过程
热失控是指电池内部产热速率远大于散热速率，导致温度持续上升，进而加速副反应放热，形成正反馈直至电池起火或爆炸的现象。理解热失控的触发条件与演变阶段，对于设计热管理系统的安全策略至关重要。
热失控的演变可分为四个阶段：
- SEI膜分解阶段（80~120℃）：当电池局部温度超过SEI膜的热稳定温度，SEI膜开始分解，释放热量并暴露负极活性表面。此阶段产热速率较慢，但为后续反应提供了触发条件。
- 负极-电解液反应阶段（120~200℃）：SEI膜分解后，电解液与嵌入锂的石墨直接接触，发生剧烈放热反应，产生可燃气体（H₂、C₂H₄等）。此阶段温升速率明显加快。
- 正极分解阶段（200~300℃）：三元正极材料发生晶格氧释放，与电解液发生燃烧反应，释放大量热量与氧气。此阶段温升速率可达10℃/s以上。
- 电解液燃烧与热失控蔓延（>300℃）：电解液在高温下燃烧，热量向相邻电芯传递，引发模组级热失控蔓延。
2.4 案例拆解：某三元锂电池快充产热仿真与验证
为将上述理论应用于实际工程，以下给出一个完整的案例闭环。某车企开发一款搭载80kWh三元锂电池包的纯电动车型，电池采用NCM811体系，额定容量120Ah，标称电压3.65V。在4C快充（480A）工况下，电池包温升过快，需要建立产热模型并优化热管理策略。
初始现象与目标：实测4C快充10分钟后，电芯中心温度达到52℃，表面温度44℃，温差8℃。目标是将电芯最高温度控制在45℃以内，温差控制在5℃以内。
关键输入数据：
- HPPC测试获得不同SOC（10%~90%）与温度（-20℃~50℃）下的内阻数据，4C倍率下总内阻约0.8mΩ（25℃，50%SOC）
- 等温量热仪测得熵热系数∂OCV/∂T在SOC=30%时达到峰值-0.4mV/K
- 电芯热物性：k_r=28W/(m·K)（面内），k_z=0.8W/(m·K)（厚度方向），Cp=950J/(kg·K)，ρ=2400kg/m³
分析判断过程：采用Bernardi方程计算产热率。在4C、50%SOC、25℃条件下：
不可逆热：Q_irr = I²·R_total = 480² × 0.0008 = 184.3W
可逆熵热：Q_rev = I·T·(∂OCV/∂T) = 480 × 298 × (-0.0004) = -57.2W
总产热：Q_gen = 184.3 - 57.2 = 127.1W可见在50%SOC附近，熵热为吸热（负值），对温升有抑制作用。但在SOC=80%以上，熵热系数转为正值，产热率将显著增加。进一步建立三维CFD模型，将产热率作为体积热源加载至电芯几何模型中，考虑冷却板与导热胶的接触热阻。
实际操作步骤：
- 在STAR-CCM+中建立电芯-导热胶-冷却板-流道的共轭传热模型，网格量约800万
- 将HPPC内阻数据与熵热系数表导入用户自定义函数（UDF），实时计算产热率
- 设置冷却液入口温度25℃，流量8L/min，模拟4C快充全过程
- 对比仿真与实测温度曲线，校准接触热阻与冷却板换热系数
结果对比：仿真预测电芯中心最高温度50.8℃，与实测52℃偏差1.2℃，误差2.3%，满足工程精度要求。通过优化冷却流道设计（增加扰流片）与降低入口温度至22℃，仿真显示最高温度降至44.5℃，温差4.2℃，达成目标。
风险与代价：降低冷却液入口温度会增加压缩机功耗，实测整车能耗增加约3%。此外，低温下大倍率充电可能引发析锂风险，需在BMS中增加低温加热策略。
复盘经验：产热模型精度主要受内阻数据准确性影响。建议在建模前对电芯进行充分的老化测试，因为老化后内阻增加20%~30%，产热率相应上升。同时，熵热系数在SOC两端变化剧烈，模型必须覆盖全SOC范围。