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⁻¹。两者相差可达一个数量级以上,这正是圆柱电芯"径向散热快、轴向散热慢"、方形电芯"大面散热快、厚度方向散热慢"的根源。

建模约定:本文统一采用圆柱坐标系,r 表示卷芯径向(面内),z 表示卷芯轴向(垂直极片方向)。方形电芯可类比为"厚度方向 = z,大面内 = r"。切换约定时必须显式说明,否则各向异性参数会张冠李戴。

二、参数测试方法

密度最容易获得:用电子天平测电芯质量 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 上。

电池热物性参数:物理角色 · 测试手段 · 换算关系 密度 ρ kg·m⁻³ 2200 ~ 2900 质量 / 体积 比热容 cp J·kg⁻¹·K⁻¹ 800 ~ 1200 DSC / 混合法 导热系数 k W·m⁻¹·K⁻¹ 面内 20~30 厚度 0.5~2 各向异性 kr, kz W·m⁻¹·K⁻¹ k_r / k_z ≈ 10~50 分方向辨识 天平 + 卡尺/排水 ρ = m / V DSC 热流法 蓝宝石参比 或 mcΔT = Pt 激光闪射 / 稳态热流 α = 0.1388 L²/t₁/₂ k = α ρ c_p 加热棒 + 多点热电偶 稳态温度场拟合 解耦 k_r 与 k_z 核心换算链 热扩散率 α × 密度 ρ × 比热容 c_p = 导热系数 k α 由激光闪射测得,ρ 与 c_p 独立测得,三者误差逐级传递至 k 误差传递提示 若 α、ρ、c_p 各含 5% 误差,k 的相对误差可达约 15%;各向异性辨识中 k_r 与 k_z 的耦合会进一步放大偏差 约定:r = 卷芯径向(面内),z = 卷芯轴向(垂直极片方向)

四、测试中的关键控制点

第一,样品代表性。DSC 样品只取几毫克,若只测正极涂层,得到的 cp 不能代表整芯;整芯建模应优先采用混合法或按质量分数加权。第二,温度依赖性。cp 与 k 随温度变化,25 ℃ 与 60 ℃ 下可相差 5%~15%,宽温域模型应测多个温度点并拟合温度系数。第三,接触热阻。稳态热流法中样品与夹具之间的接触热阻会直接叠加到测量结果上,需涂导热硅脂并施加恒定压力。第四,各向异性辨识的测点布置。测点应同时覆盖径向与轴向,且避免全部落在等温线上,否则 kr 与 kz 无法解耦。

常见误区:把供应商提供的"材料导热系数"直接当作电芯等效导热系数。极片涂层的 k 可能只有 1~3 W·m⁻¹·K⁻¹,但卷芯面内因集流体并联,等效值可达 20 W·m⁻¹·K⁻¹ 以上。二者不可混用。

五、参数汇总示例

下表给出一组用于教学示例的圆柱电芯参数(属于合理假设,非实测数据),可作为建模输入参考:

参数 符号 示例取值 单位 测试方法
等效密度ρ2550kg·m⁻³天平 + 排水法
等效比热容cp1050J·kg⁻¹·K⁻¹整芯混合法
径向热导率kr24W·m⁻¹·K⁻¹激光闪射 + 辨识
轴向热导率kz1.2W·m⁻¹·K⁻¹稳态热流 + 辨识
体积热容ρcp2.68×10⁶J·m⁻³·K⁻¹由 ρ 与 cp 相乘