J1-06 电池热管理一维建模与快充温升仿真
课程代码 J1-06 · 板块 J 仿真分析 CAE / J1 一维系统仿真 时长 约 4.0 小时(5 讲 + 1 次「Bernardi 生热 + 液冷冷板 + 冷却回路 → 4C 快充温升 → 流量扫描定最小流量」实操) 适合对象 电池热管理 CAE 工程师(主修)· 快充冷却设计工程师(必修)· 电池系统与充电策略工程师(选修)· BMS 热保护标定工程师(选修) 前置 D1-01 动力电池产热机理与产热速率建模 动力电池产热机理与产热速率建模(Bernardi 两项、熵热系数、内阻的脉宽口径)· J1-03 冷却液回路水力与热仿真建模 冷却液回路水力与热仿真建模(ε-NTU 降阶、UA 是工况量、稳态能量守恒口径)· A1-06 热阻网络法与集总参数建模入门 热阻网络法与集总参数建模入门 · D1-07 电芯热学物性的实测与索取:c_p、各向异性 k 与内部等效热阻的取值来源、测法与验收 电芯热学物性的实测与索取——⚠ 该前置在建,读者点过去目前是空的;本课不代讲它的内容,凡涉及有效比热 c_p、等效导热与内部等效热阻的取值,本课一律按待定输入处理,只写「它进哪个式子、给错了会怎样、去哪儿拿」 本讲义定位 讲师授课蓝本 / 学员自学讲义。本课是 J1 板块里唯一把电化学产热源接进热回路的一门——J1-03 教你把回路解出来,本课教你把一个会随温度和电流反过来变的热源挂上去,然后在 10 分钟量级的窗口里把温度算准。⛔ 本课不做电芯选型、不做充电协议设计、不做冷板流道 CFD,只交出「这版冷却能撑住多大电流、峰值落在哪、温差从哪儿来」这三个数与它们的口径
按任务导航
| 你手上的事 | 直接看 | 还要配 |
|---|---|---|
| 手上有电芯规格书,要搭第一版电—热耦合模型 | 第 1 讲 1.2~1.5:ECM 接 Bernardi 的两项与符号约定 | 第 5 讲 5.2:内阻按什么脉宽取 |
| 模型跑出来温升明显偏低,要查是哪一环 | 第 2 讲 2.4:把 UA_cp 直接当散热系数用是最常见的一处高估 | 第 3 讲 3.4:四条曲线的偏差方向对照表 |
| 要回答「这版冷却最多支持多大持续充电电流」 | 第 4 讲 4.2:Q̇(I) 抛物线与散热能力直线求交,本课算到 4.05C | 第 4 讲 4.3:流量加到无穷大也只到 4.57C,因为固体侧封顶 |
| 评审问「6C 能不能上」 | 第 4 讲 4.4:稳态无解、瞬态可行——两个答案都对,问的不是同一件事 | 图14:两条判据线画在同一张图上 |
| 包内温差超标,要判该改冷却拓扑还是改 TIM | 第 3 讲 3.6:三项预算拆账,本课算例里介质侧独占 8.14 K | D1-02 电芯温度不均与温差控制目标 温差口径四要素;D2-05 冷板与电芯的接触热阻与 TIM 应用 TIM 那一级 |
| 要给 BMS 标热保护门限提供仿真依据 | 第 4 讲 4.6:降功率一旦介入,稳态守恒式就不成立了 | 第 4 讲 4.5:预热到可充温度的能量与时间账 |
| 要判「析锂风险」在仿真里怎么交 | 第 4 讲 4.4:交的是 (T_min, I) 轨迹去查包络,⛔ 不是交一个 ΔT_pack | D4-04 高倍率工况下的温差与析锂风险控制 析锂判据本体;本课只做仿真侧接口 |
| 模型和实测对不上,要判是内阻错还是熵热系数错 | 第 5 讲 5.3:双倍率量热把两个未知量分开 | 第 5 讲 5.4:单点回标为什么必然把误差全记到内阻头上 |
| 要把这个电池模型接进整车联合仿真 | 第 5 讲 5.6:四件交付物与接口契约 | ⛔ 联合仿真本体 → J1-04 整车热管理系统级联合仿真 |
| 冷板各支路流量不均,要判它进不进这个模型 | 第 2 讲 2.6:集总模型里没有这个自由度 | ⛔ 流量均匀性本体 → J3-02 电池包液冷流道 CFD 与流量均匀性 |
| 充电结束停在原地,模型说温度几乎不掉 | 第 2 讲 2.7:UA_amb 被当成零了;三态取值不同 | ⛔ 分路径 UA 拆账 → D3-03 电池包保温设计与静置热损 |
| 要写模型的可用边界声明 | 第 5 讲 5.5:外推到什么倍率、SOH、环温之外就不保证 | 第 5 讲 5.7:温度包线与试验相关性 → J7-01 仿真—试验相关性验证(Correlation) |
⛔ 本课不覆盖 · 栏一(每样带去向):产热八类项的完整全集、三支极化按时间尺度的分离、熵热系数的电化学来源 → D1-01 动力电池产热机理与产热速率建模;包内温差的四条来源链、温差指标的口径四要素与验收写法 → D1-02 电芯温度不均与温差控制目标;电芯有效比热、各向异性导热与内部等效热阻的测法与验收 → D1-07 电芯热学物性的实测与索取:c_p、各向异性 k 与内部等效热阻的取值来源、测法与验收(⚠ 在建);TIM 那一级的分解式、BLT 公差链与选型树 → D2-05 冷板与电芯的接触热阻与 TIM 应用;冷板流道设计、板内流量均匀性与压降 → D2-01 电池液冷冷板设计:流道、压降与温差均匀性 / J3-02 电池包液冷流道 CFD 与流量均匀性;直冷系统的膨胀阀控制与过热度策略本体 → D2-02 电池直冷(制冷剂直冷)系统设计与控制;充电协议本身怎么定、充电桩侧的功率协商 → D4-01 快充产热与温升预测建模 / D4-03 充电过程温控策略与充电功率协同;析锂机理、析锂包络怎么测、串并联方向的自限差别 → D4-04 高倍率工况下的温差与析锂风险控制;加热手段选型与脉冲自加热 → D3-01 电池加热策略:PTC / 热泵 / 液热 / D3-02 脉冲自加热(交流自加热)低温快速升温;保温叠层、分路径 UA 拆账与热桥清单 → D3-03 电池包保温设计与静置热损;车底外流场里对流换热系数怎么提取 → J2-03 整车外流场与气动—热耦合仿真;热失控、热扩散与安全边界 → D5-01 热失控热扩散抑制(热管理与结构协同);整车级联合仿真的接口与工况链 → J1-04 整车热管理系统级联合仿真;温度包线与试验相关性判据 → J7-01 仿真—试验相关性验证(Correlation);水力网络求解、泵曲线与定压点 → J1-03 冷却液回路水力与热仿真建模;制冷剂侧循环建模与物性调用 → J1-02 制冷循环一维建模与制冷剂物性。
⛔ 本课不覆盖 · 栏二(真正不给的数):① 任何电芯的真实内阻、熵热系数与有效比热——本课全部用教学假设值,⛔ 一律不得取用,取源路径见 5.2;② 析锂包络(不同 SOC / 温度下的最大允许充电电流)——它只能来自电芯供应商的三电极或参考电极试验,⛔ 本课不给任何一条线,只给「拿它来干什么、仿真侧要交什么去查它」;③ 包内温差的验收限值——业内常用的工程经验目标是把包内最大温差压在 5 K 量级以内,这是经验值不是法规值(承 D1-02 电芯温度不均与温差控制目标),须用自家电芯的老化与析锂数据校核,⛔ 不得当验收线直接抄;④ 降功率门限与恢复迟滞的具体温度值——它是整车标定量,本课只给「门限一旦介入,模型里哪条恒等式就失效」;⑤ 标准条款的具体内容与限值——本课引用的标准一律只写它管什么与去向,⚠ 版次、实施日期与代替关系以现行目录为准,引用前须核对原文,⛔ 正文任何一处都不把版次单值化成一个确定的年份。
⚠ 栏二之外的平台相关量(整包有效热容、冷板等效 UA、包体漏热 UA_amb 三态、冷却液流量、温度上限与降功率门限、充电台阶划分)本课照给典型参考区间与常见值锚点,逐条标「⚠ 典型值、⛔ 不是验收依据,须实测确认」与拿到它的路。汇总表在 1.8 与 2.8。
引言:一份跑得很顺、峰值温度却低了 3 K 的模型
一位工程师接到的任务很清楚:给一个 800 V 平台的电池包做 4C 快充温升仿真,回答「20%→80% SOC 这一程,包最热的那颗芯到几度」。他手上有电芯规格书(容量、标称电压、25 ℃ 下的直流内阻)、有冷板的样本 UA、有冷却回路的流量。工具是现成的一维系统仿真平台,电池库里有 Bernardi 生热模块,冷板库里有换热器模块,连线、给边界、按运行。
模型三分钟就跑完了,曲线很漂亮:温度从 25 ℃ 平滑升到 30.5 ℃,没有振荡,没有不收敛,峰值出现在充电结束那一刻。他把这条曲线交了出去,结论是「峰值 30.5 ℃,距离 45 ℃ 上限余量充足,4C 冷却没有问题,可以考虑往 5C 探一探」。
三个月后台架实测:峰值 34.9 ℃,而且峰值不在充电结束,出现在 SOC 70% 附近——也就是电流已经从 4C 降到 3C 之后的那一段。差了 4.4 K,位置也不对。更麻烦的是「往 5C 探」这条建议:实测把电流台阶抬上去之后,包温直接撞上降功率门限,充电时间不降反升。
复盘发现的不是一个错,是四个方向不同的错叠在一起,而它们互相抵消了一部分,所以模型曲线看起来不但不荒唐,还相当「合理」:
- 冷板的样本 UA 被直接当成了散热系数——他写的是 Q̇_cool = UA_cp·(T_cell − T_进液),而冷却液流过冷板是会被加热的,真正该进这个式子的是 ε·C_cool,比 UA_cp 小。⇒ 高估散热,低估温升(本课算例里这一项值 0.98 K)。
- 内阻取的是 1 s 脉宽的 DCIR——因为规格书上只有这一个数。而 4C 台阶要持续 270 s,该用的是长脉宽的内阻。⇒ 低估产热,低估温升(本课算例里 2.50 K)。
- 内阻按 25 ℃ 定死没随温度修正——包温实际升到了 31 ℃,那里的内阻更低。⇒ 高估产热,高估温升(1.22 K,方向与前两项相反)。
- 熵热项被漏掉了,只算了 I²R——这颗电芯在充电段的熵热是吸热的。⇒ 高估产热,高估温升(1.68 K,方向也相反)。
四项净下来只剩 −1.5 K 左右,加上边界条件与热容的其余偏差,最终落成了 4.4 K。这就是本课要处理的那件事:电—热耦合模型的偏差极少是「差了一点」,它是一组方向已知、各自可算的系统性偏差在互相抵消。互相抵消的后果比单向偏差更坏——单向偏差你换一个工况就会发现,抵消掉的偏差在标定工况上对得很好,换一个倍率就散架,而它散架的时候你已经把「可以往 5C 探」写进了报告。
峰值位置的那件事也不是巧合。峰值出现在 SOC 70%(3C 段的末尾)而不是充电结束,是因为热时间常数比充电窗口还长:这个包的 τ = m·c_p/UA_eff ≈ 12.7 min,而整程只有 11.5 min。温度根本来不及跟上电流台阶,它一直在追一个不断下移的稳态目标,于是峰值落在「热量还在往里灌、而稳态目标已经开始下移」的那个交点上。⇒ 快充温升本质上是一个瞬态问题,用「算一个 UA 够不够」的稳态思路去做,不是精度差一点,是问错了问题(这条判据承 D4-01 快充产热与温升预测建模,本课把它变成一个可以手算的数)。
本课的两根钉子,先在这里立起来。
钉子①:电—热耦合模型的每一处简化,都对应一个方向可判、量级可算的偏差;简化是允许的,但方向与量级必须写进交付说明。这句话的重量在后半句。本课不会告诉你「必须建熵热」「必须做温度修正」——那要看你的倍率和精度要求。本课要求的是:每断开一处耦合,你都得能当场说出「这让温度偏高还是偏低、在这个算例里值几 K」。说不出方向的简化,等于把一个未知量伪装成了零。全课用同一组教学假设参数,把七处简化的偏差逐个算成一个 K 数,汇总表在 3.4。
钉子②:冷却对充电功率的约束有两条不同的天花板,而它们不是同一个数。第一条是稳态天花板——产热与散热持平时的电流,它由 UA_eff·(T_上限 − T_进液) 定;第二条是瞬态天花板——在给定时长内热容能吞下多少热而不越限,它由 m·c_p 与时长定。本课算例里稳态天花板是 4.05C,而 6C 的稳态方程根本无解(需要的 UA_eff 是固体侧物理上限的 1.78 倍);可 6C 充 180 s 的瞬态峰值只有 36.2 ℃,安全得很。⇒ 「6C 能不能上」这个问题,答案取决于问的是哪条天花板,而两条天花板对应完全不同的设计动作:稳态不够要动界面与冷源,瞬态不够要动热容与预调温。⛔ 把两条混成一条,是超充热管理里最贵的一类误判。
★ 还要先点破一条本课最容易被读反的说法:「加大冷却液流量就能压住快充温升」——这句话在本课算例里成立,但成立的原因和绝大多数人以为的不是同一条。多数人以为加流量是把冷板的 UA_cp 加上去了;本课算例里,流量从 15.0 加到 45.0 L/min(三倍),UA_cp 只从 686.1 涨到 733.5 W/K(+6.9%),因为控制热阻压倒性地在固体侧(电芯内部 + TIM + 界面 + 板壁占 93.33%)。真正被流量改掉的是另外两样:冷却液自己的温升 ΔT_coolant = Q̇/(ṁ·c_p)(从 16.28 K 降到 5.43 K)和由它派生的效能 ε。⇒ 这条误读会导出一个具体的错误动作:他会去换一台压力更高的水泵,而不是去改流道拓扑或者收 TIM 的公差——因为他以为收益在 UA 上,而 UA 那一侧的天花板是 771.4 W/K,已经近在眼前。第 2 讲 2.5 与第 3 讲 3.5 各兑现一次。
第 1 讲 电池产热怎么建:ECM、Bernardi 两项,与它们随 SOC / 温度 / 倍率的漂移
1.1 这一讲要立的东西
一维快充仿真的入口是产热源项。它写错,后面所有的冷却设计、流量扫描、功率协同全部建在流沙上——而产热源项写错是不会报错的,它只会给你一条形状对、幅值错的温升曲线。这一讲要立四件:
- 符号约定——本课沿用 D1-01 动力电池产热机理与产热速率建模 的「放电为正」约定,并给出充电段的书写法;
- 两项的身份——Q_irr 从 ECM 出,Q_rev 从熵热系数表出,两者数据来源完全不同;
- 门槛判据 I* = T·|∂U_oc/∂T|/R——它一句话回答「这个倍率下熵热能不能忽略」;
- 三个自变量的漂移——内阻随 SOC、温度、脉宽变,取定值的偏差方向可判。
1.2 符号约定:本课全程「放电为正」,充电段用 I_ch 写
本课与 D1-01 动力电池产热机理与产热速率建模 共用同一套约定,⛔ 不另起一套:
- 电流 I:放电为正。⇒ 快充段 I < 0。
- 为免全篇写负号,本课引入 I_ch ≡ −I > 0,读作「充电电流幅值」,单位 A。凡本课写 4C、400.0 A,指的都是 I_ch。
- 熵热系数 ∂U_oc/∂T:单位 V/K 或 mV/K,本课教学取值 −0.100 mV/K(与 D1-01 动力电池产热机理与产热速率建模 教学算例同值同口径,⚠ 教学假设值,⛔ 不对应任何电芯)。
⚠ 这一条必须在模型文档第一页写死,⛔ 不许口头约定。为什么:换成「充电为正」,Bernardi 熵热项的吸放热标签整行互换,而模型照常收敛、曲线照常平滑。约定错的模型和约定对的模型,在一条 4C 曲线上的差别是 1.68 K——足够被读成「网格不够细」或者「热容取偏了」,而不会被读成「符号错了」。
易错点①:模型里用「充电为正」而参数表里的熵热系数按「放电为正」测出,两处不一致。为什么会错:熵热系数的符号约定通常写在试验报告的方法一节而不是数据表旁边,参数录入的人只看到一个 −0.100 mV/K。怎么当场发现:跑一次 0.5C 恒流充电、绝热边界,看包温是升还是降——本课算例里这个工况净产热为 −31.8 W(净吸热,温度应当缓慢下降);若模型给出温升,符号约定必然有一处反了。
1.3 Bernardi 两项:一项从 ECM 出,一项 ECM 给不了
Bernardi 生热的完整式子是:
Q̇ = I·(U_oc − V) − I·T·(∂U_oc/∂T)
两项的身份完全不同,这一点决定了它们的数据来源和出错方式都不一样。
第一项 · 不可逆热(极化 + 欧姆)。I·(U_oc − V) 就是「实际端电压偏离开路电压的那部分乘以电流」。它恒为正——放电时 V < U_oc,充电时 V > U_oc,两种情况下这一项都放热。在一维模型里,它通常被写成 ECM 里三个电阻上的耗散之和;若只保留欧姆项并把极化折算进去,就退化成工程上最常用的那个写法:
Q̇_irr ≈ I²·R_int
这里的 R_int 是一个把欧姆与被建模到的那几支极化打包在一起的等效内阻,它的值取决于你按多长的脉宽去测(1.6 展开)。
第二项 · 可逆热(熵热)。−I·T·(∂U_oc/∂T) 是电极反应熵变的直接后果,⛔ 不是一个拟合参数(这一条承 D1-01 动力电池产热机理与产热速率建模,本课不重推)。它的符号由两个因子共同决定:电流方向,以及熵热系数本身的符号。按本课约定(放电为正、I_ch = −I):
Q̇_rev = +I_ch·T·(∂U_oc/∂T)
取本课教学值 ∂U_oc/∂T = −1.00×10⁻⁴ V/K,⇒ 充电段 Q̇_rev < 0,吸热。
⚠⚠ 这里有一个整类的陷阱,必须点破:ECM 无论做到几阶,都算不出 Q_rev。等效电路模型的每一个元件都在描述「电压相对于 U_oc 的偏离」,而熵热来自 U_oc 本身随温度的移动——它在 ECM 的建模对象之外。⇒ 你把二阶 ECM 换成三阶、把 RC 参数辨识得再准,Q_rev 都仍然是零。而模型不会报错:它会给你一条更平滑、更「精细」的温升曲线,而那条曲线在低倍率段的符号都是错的。这是「没有被建进模型的现象,模型不会报错,它会如实地报告这个现象不存在」的一个具体实例。
算例 1-A(教学算例,全部为教学假设值)。取本课贯穿全篇的教学包 TP-1(参数表见 1.8):192S1P、单体 100.0 Ah、单体内阻 R_cell = 0.550 mΩ(25 ℃ / 50% SOC / 10 s 脉宽)⇒ 包内阻 R_pack = 192 × 0.550 mΩ = 105.6 mΩ;熵热系数 −1.00×10⁻⁴ V/K;温度取 35.0 ℃ = 308.15 K。
4C 充电(I_ch = 400.0 A):
| 项 | 代入 | 结果 |
|---|---|---|
| 不可逆热 Q̇_irr | 400.0² × 0.1056 Ω | +16 896 W = +16.90 kW |
| 可逆热 Q̇_rev | 192 × 400.0 × 308.15 × (−1.00×10⁻⁴) V/K | −2 367 W = −2.37 kW |
| 净产热 Q̇_gen | 16 896 − 2 367 | +14 529 W = 14.53 kW |
⇒ 熵热在 4C 下把产热削掉了 14.0%(分母是 16 896 W 这一项)。结论怎么用:漏掉熵热的模型在 4C 段高估产热 16.3%(分母换成净值 14 529 W),折成峰值温度是高估 1.68 K(3.4 逐项算)。这是一个偏保守的偏差——它不会让你把一个不安全的方案判成安全,但它会让你多买冷却能力。
算例 1-B(同一组参数,换 0.5C)。I_ch = 50.0 A:
| 项 | 代入 | 结果 |
|---|---|---|
| Q̇_irr | 50.0² × 0.1056 | +264.0 W |
| Q̇_rev | 192 × 50.0 × 308.15 × (−1.00×10⁻⁴) | −295.8 W |
| 净产热 | 264.0 − 295.8 | −31.8 W(净吸热) |
⇒ 同一颗电芯、同一组参数,0.5C 充电净吸热,4C 充电净放热 14.53 kW。漏掉熵热的模型在 0.5C 段给出 +264.0 W,连符号都是错的。这就是为什么「熵热在快充里可以忽略」这句话必须连条件一起说——它在 4C 成立(值 14.0%),在慢充段完全不成立。
1.4 门槛判据 I*:一句话回答「这个倍率下熵热能不能扔」
把两项的比值写出来(承 D1-01 动力电池产热机理与产热速率建模 的归纳,本课直接取用):
|Q̇_rev| / Q̇_irr = T·|∂U_oc/∂T| / (I·R_cell) = I* / I,其中 I* = T·|∂U_oc/∂T| / R_cell
vars:T 电芯绝对温度(K)· ∂U_oc/∂T 熵热系数(V/K)· R_cell 单体等效内阻(Ω)· I 电流幅值(A)· I* 门槛电流(A),即两项恰好相等的那个电流。 适用条件:恒流、该段内 R_cell 与 ∂U_oc/∂T 近似不变、单体与包同比例(串联包两侧都乘 N,比值不变)。 失效边界:① ∂U_oc/∂T 在某些 SOC 点会变号,穿零附近 I* → 0,此式给出的「熵热可忽略」在那一小段是自动成立的、⛔ 不能外推到整个 SOC 区间;② R_cell 随温度变化大(低温段尤甚),I* 随之漂移,⛔ 不是一个常数;③ 该式只比较两项幅值,⛔ 不告诉你符号,扔与不扔的后果还要看 1.3 那张符号表。
算例 1-C。I* = 308.15 × 1.00×10⁻⁴ / 5.50×10⁻⁴ = 0.030815 / 0.00055 = 56.03 A = 0.56C。
⇒ 判据一句话:工作电流在 0.56C 以下时熵热是主项,以上时不可逆热是主项,比值按 I*/I 反比下降。4C 处比值 = 56.03/400.0 = 0.140(与算例 1-A 的 14.0% 逐位对上);6C 处 = 0.093。
易错点②:把「快充里熵热可忽略」当成无条件全称,于是整个模型库里都不建熵热。为什么会错:这个结论是在高倍率段量出来的,而同一份模型会被拿去跑慢充、跑标定、跑热管理待机段,那些工况全部落在 I* 左边。怎么当场发现:拿模型跑一次 0.3C 绝热充电,若温度单调上升且斜率与 I² 成比例,熵热必然没建。
易错点③:把 I* 当成一个电芯常数写进内部规范。为什么会错:I* 的分子含 T 与 ∂U_oc/∂T(后者随 SOC 变号)、分母含 R_cell(随温度强烈变化),三个都在动。怎么当场发现:在同一份规范里同时出现「I* = 某个固定安培数」与「内阻随温度修正表」,两者必然矛盾——把修正表的两端代进去,I* 会给出两个差很远的值。
1.5 产热随 SOC / 温度 / 倍率变:三个自变量,三种错法
R_int 不是一个数,它是一张至少三维的表:R_int(SOC, T, 脉宽),而工程上常常还要加第四维 SOH。⚠ 快充仿真里最常见的做法是取一个 25 ℃ / 50% SOC / 10 s 的单值,这个做法可以,但三处偏差的方向必须写出来:
| 简化 | 偏差方向 | 本课算例量级(4C,峰值温度) | 为什么是这个方向 |
|---|---|---|---|
| 不随温度修正(固定 25 ℃ 值) | 高估产热、高估温升 | +1.22 K | 内阻随温度下降;包温升到 31.4 ℃ 时内阻已比 25 ℃ 低约 11.5%(教学取 −1.8%/K) |
| 不随 SOC 修正(固定 50% 值) | 方向取决于 SOC 区间,⛔ 不可写成单向 | 本课不给单值 | 内阻在两端(低 SOC 与高 SOC)通常抬头,中段平坦;20%→80% 这一程主要走中段,偏差小;若程序段延伸到 90% 以上,偏差转为低估 |
| 脉宽取短(1 s 而非 10 s) | 低估产热、低估温升 | −2.50 K | 短脉宽只抓到欧姆与快极化,慢极化那一支还没建立起来 |
算例 1-D(温度修正的量级)。教学取内阻温度系数 −1.8%/K(⚠ 教学假设值,⛔ 不对应任何电芯,真实值须由该电芯的 DCIR-温度表给出)。包温从 25.0 ℃ 升到 31.4 ℃(本课基线工况的峰值,见 3.3),ΔT = 6.4 K ⇒ R_cell = 0.550 × (1 − 0.018 × 6.4) = 0.4866 mΩ,降 11.5%。据此重算 4C 段净产热:Q̇_irr = 400.0² × 192 × 4.866×10⁻⁴ = 14 950 W,Q̇_rev 不变 −2 367 W ⇒ Q̇_gen = 12 583 W,比固定 25 ℃ 的 14 529 W 低 13.4%。折成峰值温度:30.22 ℃ vs 31.44 ℃,固定内阻高估 1.22 K。
⇒ 这一条怎么用:温度修正是偏保守方向的简化,快充定尺寸阶段可以不做(它给你余量),但做相关性对标时必须做——否则你会把模型偏高的那 1.22 K 归因到冷却系数上,然后去调一个本来就对的 UA。
算例 1-E(脉宽的量级)。教学取同一电芯三档脉宽内阻:R(1 s) = 0.420 mΩ、R(10 s) = 0.550 mΩ、R(30 s) = 0.630 mΩ(⚠ 教学假设值)。4C 段用 1 s 值:Q̇_irr = 400.0² × 192 × 4.20×10⁻⁴ = 12 902 W ⇒ Q̇_gen = 10 536 W,比 10 s 值低 27.5%。折成峰值温度:28.94 ℃ vs 31.44 ℃,低估 2.50 K。
⚠ 脉宽该取多长,判据是「你这一段恒流持续多久」,⛔ 不是「规格书上有哪个」。本课 4C 台阶持续 270 s ⇒ 10 s 值已经偏短,30 s 值更接近;再长的脉宽会开始把 SOC 漂移与温升混进去,那要靠 5.3 的量热回标而不是脉冲。⇒ 这是本课的一条判据:脉冲宽度的选择就是在选「这个 R 代表哪几支极化」(承 D1-01 动力电池产热机理与产热速率建模)。
1.6 单体 → 模组 → 包:三档热网络,各自能回答什么问题
产热建好之后,下一个决定是把热挂在几个节点上。这不是精度问题,是自由度问题——每一档降阶都会让某一整类问题在模型里「不存在」(承 A1-06 热阻网络法与集总参数建模入门 的热阻网络方法本体,本课只做快充场景的取舍)。
| 档 | 节点数(本课 TP-1) | 能回答 | 在模型里不存在的问题 |
|---|---|---|---|
| 集总一节点 | 1 | 包平均温升、峰值出现时刻、冷却能力够不够 | 包内温差恒为零;「哪颗芯最热」得到的答案是「没有这回事」 |
| 分区多节点 | 6~12(按冷板流路分区) | 沿程温升带来的区间温差、进出口两端的差别 | 区内芯间差;同一区内 TIM 公差分散看不到 |
| 每芯一节点 + 芯内分层 | 192 × 2~3 | 芯间温差、芯内表面—芯部梯度 | 电芯内部真实三维场;极耳局部热点(那是三维的事) |
⚠ 这张表最要紧的是最右列。集总一节点的模型跑快充,出水温度会算得好好的,而包内温差恒为零——不是算得不准,是模型里没有一个状态量在记录它。⇒ 若你的交付要回答「ΔT_pack 是多少」,集总模型无论怎么调参数都给不出来,必须升到分区档,或者用 3.6 那条预算式在集总模型外面另算一遍。
本课的取法:主线用集总一节点(因为快充的第一问是「峰值到几度」),温差用 3.6 的预算式外挂。这是一个刻意的降阶,它的损失就写在上表最右列,⛔ 不许拿集总模型的输出去回答温差问题。
1.7 常见错法与它们的自查动作
易错点④:把电芯外的欧姆热(汇流排、连接片、极柱、高压接插件)摊进 I²R_int。为什么会错:R_int 的口径是「电芯两个极柱之间」,而汇流排热是集中源,它落在结构件上、走的散热路径与电芯本体不同;摊进去等于把一份本该出现在结构件节点上的热挪到了电芯节点上,直接抬高电芯温度。怎么当场发现:把包内阻与「192 × 单体内阻」比一比——若模型里的包内阻明显大于这个乘积,多出来的那一段就是被摊进来的电芯外欧姆热,应当单独建节点。本课 TP-1 的 R_pack = 105.6 mΩ 严格等于 192 × 0.550 mΩ,电芯外欧姆热另建(教学取 8.0 mΩ 挂在结构节点上,4C 时 400.0² × 0.008 = 1 280 W)。
易错点⑤:用容量与电压算出的「充电功率 × (1 − 效率)」当产热。为什么会错:那个「效率」是一个把电芯、汇流排、接插件、充电机通讯损耗全打包的整机数,它的分母与边界都不是电芯;而且它是一个标称值,不随倍率变。怎么当场发现:用它算 0.5C 与 4C 两点,产热会严格按电流一次方变化——而真实的不可逆热按二次方变,熵热按一次方变,两者比例随电流变。⇒ 一次方的产热曲线是这个错法的指纹。
易错点⑥:把 SOC 用电压查表得到,而不用安时积分。为什么会错:快充段极化压降很大,端电压与 OCV 差得远(本课 4C 时单体压降 = 400.0 × 0.550 mΩ = 0.220 V),用电压查 SOC 会得到一个明显偏高的 SOC,于是充电台阶提前切换,产热曲线整体前移。怎么当场发现:把模型的 SOC 轨迹与「初始 SOC + ∫I dt/C_n」对一遍,快充段两者应当逐点重合,⛔ 差一个百分点以上就是查表混进来了。
易错点⑦:SOH 的影响用「内阻乘一个系数」一刀切。为什么会错:老化同时改写内阻和可用容量和熵热系数曲线,三者的变化率不同;只乘内阻系数会让同一 SOC 区间的充入电量算错,进而算错时长。怎么当场发现:看模型在 SOH 80% 档下 20%→80% 的充电时长——若它与新电芯完全相同,容量那一维必然没跟着改。⛔ SOH 分档建模超出本课射程,去向 D1-04 电池老化、寿命与温度的耦合关系。
1.8 教学包 TP-1 参数表(⚠ 全部为教学假设值)
⚠⚠ 下表是一整套为了让本课每一个算例都能在一张纸上手算复核而构造的自洽参数。它不对应任何电芯、任何电池包、任何平台,⛔ 禁止照抄取用,⛔ 不得据它反算任何真实产品的能力。凡本课出现的数字,只要能追到这张表,就是教学假设值。
| 量 | 符号 | 本课取值 | 单位 | 分类 | 拿到它的路 |
|---|---|---|---|---|---|
| 串并联 | — | 192S1P | — | 教学假设 | — |
| 单体标称容量 | C_n | 100.0 | Ah | 教学假设 | 电芯规格书 |
| 包标称电压 | U_pack | 710.4(192 × 3.70) | V | 教学假设 | — |
| 包能量 | E_pack | 71.0 | kWh | 教学假设 | — |
| 单体内阻(25 ℃ / 50% SOC / 10 s) | R_cell | 0.550 | mΩ | 教学假设 | DCIR 脉冲试验,见 5.2 |
| 包内阻 | R_pack | 105.6 | mΩ | 教学假设 | = 192 × R_cell |
| 熵热系数 | ∂U_oc/∂T | −0.100(放电为正) | mV/K | 教学假设 | 电位—温度弛豫试验,见 5.2 |
| 内阻温度系数 | — | −1.8 | %/K | 教学假设 | DCIR-温度表 |
| 整包有效热容 | m·c_p | 450.0 | kJ/K | 平台相关量,⚠ 典型值 | ⚠ 属主 D1-07 电芯热学物性的实测与索取:c_p、各向异性 k 与内部等效热阻的取值来源、测法与验收(在建);过渡做法见 5.2 |
| 单体冷却面积 | A_cell | 0.0045 | m² | 教学假设 | 冷板贴合面几何 |
| 冷却面总面积 | A_tot | 0.864(= 192 × A_cell) | m² | 教学假设 | — |
| 面热阻链合计(30.0 L/min) | R″_tot | 12.0×10⁻⁴ | K·m²/W | 教学假设 | 拆账见 2.3 |
| 冷板等效 UA(30.0 L/min) | UA_cp | 720.0 | W/K | 平台相关量,⚠ 典型值 | = A_tot/R″_tot;件级台架反算见 L1-10 冷板/液冷总成件级热性能台架试验(模拟热源法) |
| UA_cp 上限(流量 → ∞) | — | 771.4 | W/K | 教学假设 | = A_tot/R″_固体侧 |
| 冷却液密度 / 比热(50/50 乙二醇水) | ρ / c_p | 1050 / 3400 | kg/m³ / J·kg⁻¹·K⁻¹ | 行业典型区间 | 冷却液物性表,⚠ 随浓度与温度变 |
| 基线冷却液流量 | V̇ | 30.0 | L/min | 平台相关量,⚠ 典型值 | 回路水力解,见 J1-03 冷却液回路水力与热仿真建模 |
| 包体漏热(静置 / 充电停车 / 行驶) | UA_amb | 12.0 / 18.0 / 35.0 | W/K | 平台相关量,⚠ 典型值 | 分态反标,属主 D3-03 电池包保温设计与静置热损 |
| 电芯温度上限 | T_limit | 45.0 | ℃ | 平台相关量,⚠ 典型值 | 电芯规格书 + 整车标定 |
| 降功率门限 | T_derate | 42.0 | ℃ | 平台相关量,⚠ 典型值 | BMS 标定,见 D4-03 充电过程温控策略与充电功率协同 |
⚠ 标「平台相关量」的六项,本课给的是典型参考值,⛔ 不是验收依据,须由本项目实测确认。它们的常见区间:整包有效热容随包结构与冷却介质持液量在数百 kJ/K 量级(本课锚点 450.0);冷板等效 UA 随冷却面积与界面质量在数百到一千余 W/K 量级(锚点 720.0);包体漏热 UA_amb 在十几到数十 W/K 量级且三态不同(锚点 12.0 / 18.0 / 35.0);快充回路冷却液流量在十几到数十 L/min 量级(锚点 30.0);电芯温度上限常见在四十几摄氏度量级(锚点 45.0),降功率门限设在其下若干开尔文(锚点 42.0)。
关于 m·c_p 的特别声明。它的属主是 D1-07 电芯热学物性的实测与索取:c_p、各向异性 k 与内部等效热阻的取值来源、测法与验收(电芯热学物性的实测与索取),⚠ 该课在建。本课不代讲它的测法与验收,只说清三件与仿真直接相关的:① 它进的是集总温升方程的分母,⇒ 给错 10% 温升斜率就错 10%,而曲线形状完全不变、模型不会报错;② 它的口径必须写明算不算壳体、模组结构件、冷板本体与板内持液量——本课 450.0 kJ/K 是「电芯 + 模组结构件 + 冷板本体」的口径,⛔ 不含回路里冷板以外的冷却液;③ 在 D1-07 电芯热学物性的实测与索取:c_p、各向异性 k 与内部等效热阻的取值来源、测法与验收 落盘之前,过渡做法是把它当成一个待标定参数,用 5.5 的降温段反标(那一段的斜率只由 m·c_p 与 UA 两个量决定)。
本讲小结:① 符号约定「放电为正」必须写在模型文档第一页,充电段用 I_ch = −I 书写;② Bernardi 两项数据来源不同——Q_irr 从 ECM 出,Q_rev 从熵热系数表出,而 ECM 无论几阶都给不出 Q_rev;③ 门槛判据 I* = T·|∂U_oc/∂T|/R_cell 一句话回答「熵热能不能扔」,本课算例 I* = 56.03 A = 0.56C,4C 处熵热占 14.0%、0.5C 处净产热为 −31.8 W(连符号都会错);④ 内阻的三个自变量给出三个方向不同的偏差(温度 +1.22 K、脉宽 −2.50 K、SOC 方向取决于区间),其中两处方向相反会互相抵消;⑤ 热网络三档降阶,右列「在模型里不存在的问题」才是选档判据——集总模型的包内温差恒为零;⑥ 七条易错点各配一个当场自查动作,其中三条(符号约定、SOC 口径、电芯外欧姆热)只要跑一个特定工况就能露馅。⇒ 产热源项立住之后,下一步是把冷却接上去——而那里藏着本课最贵的一处高估:把冷板的 UA 直接当成散热系数用。那是第 2 讲。
后面还有 4 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做