3. 电池热物性参数及其测试方法
电池热模型能否给出可信的温度场,取决于四个基础热物性参数是否准确:密度 ρ、比热容 cp、导热系数 k,以及卷芯层面的各向异性热导率 kr、kz。它们分别回答四个问题:热惯性有多大、储热能力有多强、热量沿哪个方向传得快、以及温度变化有多"迟钝"。本节先厘清物理意义,再给出可操作的测量原理与数据处理方式。
一、四个参数的物理意义与量级
密度 ρ(kg·m⁻³)决定单位体积的质量,是热容量的乘数因子。锂离子电芯的等效密度通常在 2200~2800 kg·m⁻³ 之间,磷酸铁锂体系因正极材料密度较低,常落在 2200~2500 kg·m⁻³,三元体系可接近 2700~2900 kg·m⁻³。注意这里说的是"等效密度":电芯由活性物质、集流体、隔膜、电解液、壳体共同组成,建模时把整体质量除以整体体积即可,无需逐层区分。
比热容 cp(J·kg⁻¹·K⁻¹)表示单位质量升高 1 K 所需的热量。它直接决定温升速率:在产热功率相同的情况下,cp 越大,温升越慢。锂离子电芯的等效比热容一般在 800~1200 J·kg⁻¹·K⁻¹,常见取值约 1000 J·kg⁻¹·K⁻¹。它与密度相乘得到体积热容 ρcp(J·m⁻³·K⁻¹),这才是热模型中真正出现的系数。
导热系数 k(W·m⁻¹·K⁻¹)描述材料传导热量的能力。电芯内部是典型的各向异性结构:沿卷绕/叠片平面方向(径向 r 或面内方向),热量可沿金属集流体(铝箔 k≈237、铜箔 k≈400 W·m⁻¹·K⁻¹)快速传导,等效导热系数通常为 20~30 W·m⁻¹·K⁻¹;而垂直于极片方向(轴向 z 或厚度方向),热量必须穿过导热性差的隔膜与活性物质层,等效导热系数往往只有 0.5~2 W·m⁻¹·K⁻¹。两者相差可达一个数量级以上,这正是圆柱电芯"径向散热快、轴向散热慢"、方形电芯"大面散热快、厚度方向散热慢"的根源。
二、参数测试方法
密度最容易获得:用电子天平测电芯质量 m,用游标卡尺或排水法测外形体积 V,ρ = m/V。若需分层密度,可拆解后分别称重各组分,但工程建模一般只需整体等效值。
比热容常用两种方法。其一是差示扫描量热法(DSC):取少量电极或隔膜样品,以已知比热容的蓝宝石为参比,在恒定升温速率下测量热流差,由热流比直接得到 cp。其二是整芯混合法:把电芯放入绝热或准绝热环境,用已知功率 P 加热时间 t,测量温升 ΔT,由能量守恒 mcpΔT = Pt 反推 cp。后者更贴近整芯等效值,但要求散热损失足够小,通常需做散热修正或采用真空绝热腔。
导热系数的测量分方向进行。面内方向可用激光闪射法:对薄片样品施加激光脉冲,记录背面温升曲线,由半升温时间 t1/2 计算热扩散率 α = 0.1388·L²/t1/2(L 为样品厚度),再由 k = αρcp 得到导热系数。厚度方向则常用稳态热流法:在样品两侧建立稳定温差,测量通过的热流密度 q,由 k = qL/ΔT 计算。对整芯各向异性,也可用3ω 法或热流计法配合不同方向的温度测点反演。
各向异性热导率的工程做法是"分方向测 + 参数辨识"。以圆柱电芯为例:在表面不同高度布置热电偶,中心插入加热棒施加已知功率,记录稳态温度分布,再用二维热传导模型拟合 kr 与 kz。由于 kr 与 kz 对温度场的影响模式不同——kr 主要影响径向温差,kz 主要影响轴向温度梯度——只要测点位置合理,二者可以被解耦辨识。
三、参数关系与测试链路
下图给出四个参数的物理角色、测试手段与相互换算关系。可以看到,导热系数并非直接测量量,而是由热扩散率、密度、比热容三者相乘得到,因此 α、ρ、cp 的误差会逐级传递到 k 上。
四、测试中的关键控制点
第一,样品代表性。DSC 样品只取几毫克,若只测正极涂层,得到的 cp 不能代表整芯;整芯建模应优先采用混合法或按质量分数加权。第二,温度依赖性。cp 与 k 随温度变化,25 ℃ 与 60 ℃ 下可相差 5%~15%,宽温域模型应测多个温度点并拟合温度系数。第三,接触热阻。稳态热流法中样品与夹具之间的接触热阻会直接叠加到测量结果上,需涂导热硅脂并施加恒定压力。第四,各向异性辨识的测点布置。测点应同时覆盖径向与轴向,且避免全部落在等温线上,否则 kr 与 kz 无法解耦。
五、参数汇总示例
下表给出一组用于教学示例的圆柱电芯参数(属于合理假设,非实测数据),可作为建模输入参考:
| 参数 | 符号 | 示例取值 | 单位 | 测试方法 |
|---|---|---|---|---|
| 等效密度 | ρ | 2550 | kg·m⁻³ | 天平 + 排水法 |
| 等效比热容 | cp | 1050 | J·kg⁻¹·K⁻¹ | 整芯混合法 |
| 径向热导率 | kr | 24 | W·m⁻¹·K⁻¹ | 激光闪射 + 辨识 |
| 轴向热导率 | kz | 1.2 | W·m⁻¹·K⁻¹ | 稳态热流 + 辨识 |
| 体积热容 | ρcp | 2.68×10⁶ | J·m⁻³·K⁻¹ | 由 ρ 与 cp 相乘 |