J1-09 燃料电池系统一维建模与水/热耦合仿真

课程代码 J1-09　|　板块 J 仿真分析 CAE / J1 一维系统仿真 时长 约 4 学时　|　难度 L3 专家 适合对象 燃料电池系统仿真 / CAE 工程师（必修）；需要交付「可与台架、环境舱对标」的电堆系统模型的 FC 热管理正向设计与标定岗位（必修）；FC 控制与整车集成工程师了解模型接口、降阶档与可信域（选修） 前置 J1-01 1D 热管理仿真导论（KULI/AMESim/GT/Flowmaster） 1D 热管理仿真导论（KULI/AMESim/GT/Flowmaster）；J1-03 冷却液回路水力与热仿真建模 冷却液回路水力与热仿真建模；O1-01 燃料电池电堆产热与水/热管理耦合 燃料电池电堆产热与水/热管理耦合；O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角） PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径；O2-02 膜增湿器热湿管理与水平衡 膜增湿器热湿管理与水平衡 本讲义定位 只交建模工艺：电堆与 BOP 的物理结论怎么降成 1D 模型里的节点、源项与参数，水平衡怎么在求解器里闭合并逐步自检，冷启动怎么从「跑通一次」变成一个能与环境舱比的通过率，以及交付时必须随行的六项声明与可信域。它不重讲 FC 物理（O1-01 燃料电池电堆产热与水/热管理耦合 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角））、不重讲部件热设计（O2-01 空压机冷却与进气中冷 O2-02 膜增湿器热湿管理与水平衡 O2-03 氢气循环泵/引射器与供氢温控）、不讲三维方法（只给「何时必须上三维」的判据）、不讲通用数值内核与通用相关性方法学。

引言：这门课要防的，是一个「在额定点对得很好」的模型

先讲一个反复发生的失败，它几乎定义了本课的全部内容。

某 FCEV 系统的一维模型在额定点与台架对得很好——出堆温度差不到 1 K，散热器出口也对得上。它被拿去做冬季长坡持续大负荷的散热能力核算，算出来的出堆温度有余量，散热器就按这个结果定了尺寸。实车上路，反而是这个工况先撞出堆温度上限。

事后拆开看，根因是两处叠加，而且两处都不会让模型报错：

第一处：产热用的是从「LHV 口径效率数」反推出来的值，而模型里的热中性电压填的是 HHV 的那一个。基准混用，产热被系统性算小三成量级（换算的唯一归属是 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角），本课只负责把基准写死）。

第二处：模型里没有单独建阴极尾气带热这条支路，把整堆产热全额当成了冷却回路的负荷。这两处误差方向相反——第一处把产热算小，第二处把进入冷却液的份额算大——于是在额定点上，它们互相抵掉了一部分；剩下的那点差，标定时被一个「把堆侧换热能力调大」的动作补平了。

这个补偿参数在额定点附近是有效的。一旦换到低环温、尾气潜热占比变化的工况，两处误差的抵消关系被打破，补偿就失效了。而模型不会告诉你它失效了，它只会给出一个看起来正常的数。

⇒ 这就是本课要立的第一条判断：一个一维 FC 模型的可信度，不由它在标定点上的误差决定，由它的账本闭不闭合、以及它的每一项降阶有没有被写下来决定。

这门课的四条主线

主线一：粒度是被输出量倒推出来的，不是被物理细节推出来的。「电堆分几个节点」「增湿器建到第几档」「氢路能不能扔成边界条件」——这些问题没有一个通用答案，答案由「这一次的仿真要回答什么」决定，而且每换一次输出量就要重新过一遍。第 1、2、3 讲各在一个对象上把这条判据用一遍。

主线二：FC 的两个量级特征直接改写建模取舍。散热量与电功率同量级、可用温差只有内燃机的一半（承 O1-01 燃料电池电堆产热与水/热管理耦合 O3-01 FCEV 大散热量前端散热系统（低温差挑战））⇒ 回路走「大流量、小温差」⇒ 流量的相对误差几乎全额落在进出口温差上。这是 FC 模型标定与 BEV 最大的手感差异，第 2 讲量化它。

主线三：水平衡不是一条公式，是一本必须逐时间步闭合的账。这是本课与部件课最硬的分界线：O1-01 燃料电池电堆产热与水/热管理耦合 O2-02 膜增湿器热湿管理与水平衡 给的是水通量表达式，本课负责把它们接成一个在求解器里能闭合、能自检的账本。第 4 讲给出两条守恒残差、以及「残差不回零时按什么顺序查漏项」。

主线四：冷启动结论是统计量，不是单次结果。模型能跑通一次 −30 ℃ 启动，和这个平台的冷启动成功率是多少，是两件事。第 5 讲把 O1-04 低温冷启动（-30℃）与自加热策略 的失败判据做成模型侧通过率，并说明它只能与环境舱同格子的样本比分布。

全课的符号约定（含三处跨课口径提醒，⛔ 读之前先看）

本课与 O 板块、J1 板块的兄弟课共享大量符号。下面这张表把会撞车的那几个一次钉死。

符号	本课含义	单位	跨课口径
Q_stack	电堆在给定工作点放出的全部热功率（不含 BOP）	kW	与 O1-01 燃料电池电堆产热与水/热管理耦合 O1-02 电堆冷却回路与低电导率冷却液设计 完全一致
Q_coolant	真正进入冷却液的那一份热功率	kW	★ 与 O1-02 电堆冷却回路与低电导率冷却液设计 一致，本课沿用它的写法，⛔ 不另造 Q_cool
Q_exh	阴极尾气带走的显热 ＋ 水汽潜热	kW	本课首定；它是 Q_stack 与 Q_coolant 之间那本账里最大的一个扣项
Q_amb	电堆本体对环境的散失	kW	与 O1-05 停机吹扫、保温与整堆防冻设计 的 soak 口径同源
UA_s2c	电堆固体侧 → 冷却液的等效换热能力	kW/K	⛔⛔ 本课不用 UA_stack，见下方跨课提醒①
UA_amb	电堆模块对环境的总热导	W/K	＝ O1-05 停机吹扫、保温与整堆防冻设计 的 UA_stack，本课换名以避撞，见提醒①
E_tn	热中性电压，必须带基准：E_tn(HHV) / E_tn(LHV)	V/片	承 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角），数值与它一致
E_rev	可逆电压	V/片	承 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角）
V_cell	实际单片电压	V/片	承 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角）
N_cell	电堆片数	片	—
I	电堆电流	A	全课不写 i（i 留给电流密度，A/cm²）
ṁ	质量流量	kg/s	⛔ 流量不用 Q 表示（承 O1-02 电堆冷却回路与低电导率冷却液设计 的钉子）
V̇	体积流量	L/min 或 m³/s	同上
ΔT_coolant	出堆 − 进堆冷却液温差	K	承 O1-02 电堆冷却回路与低电导率冷却液设计
ω	湿度比（含湿量），kg 水 / kg 干气	—	承 O2-02 膜增湿器热湿管理与水平衡，⛔ 本课守恒一律用它、不用 RH
λ_air	空气化学计量比	无量纲	⚠ 见跨课提醒②
λ_H2,sys / λ_H2,stack-in	系统口径 / 堆入口口径的阳极化学计量比	无量纲	承 O2-03 氢气循环泵/引射器与供氢温控，两者⛔ 不可互代
ε_m	增湿器传质有效度	无量纲	承 O2-02 膜增湿器热湿管理与水平衡
ε_s2c	电堆固体 → 冷却液换热有效度	无量纲	本课首定
n_d	电渗拖曳系数	无量纲	承 O1-01 燃料电池电堆产热与水/热管理耦合
ΔT_mush	冰相变数值温区宽度	K	本课首定，⛔ 它是数值参数不是物理参数
R_w	以生成水归一化的水量闭合残差	无量纲	本课首定
R_e	以 Q_stack 归一化的能量闭合残差	无量纲	本课首定
P_succ	模型侧冷启动成功率	无量纲	本课首定

跨课提醒①（最要紧的一条）：UA_stack 在本课不用，因为它在 O1-05 停机吹扫、保温与整堆防冻设计 已经有另一个意思。 O1-05 停机吹扫、保温与整堆防冻设计 里 UA_stack ＝ 电堆模块对环境的总热导，量级是 3~25 W/K（做了保温且关阀多在 5~8 W/K）。而本课要的是电堆固体侧向冷却液的换热能力，同一台堆上它的量级是 kW/K 级——两者相差三个数量级，量纲却完全相同（W/K），互代之后模型照样能跑、照样收敛，只是把堆冷得慢了一千倍。 ⇒ 本课把这两个量分成两个符号：UA_s2c（→ 冷却液，kW/K）与 UA_amb（→ 环境，W/K，即 O1-05 停机吹扫、保温与整堆防冻设计 的 UA_stack）。 ⇒ 引用 O1-05 停机吹扫、保温与整堆防冻设计 的 soak 实测数时，进的是 UA_amb；⛔ 不得填进 UA_s2c。

跨课提醒②：λ_air 在 FC 语境里是空气化学计量比，⛔ 不是空气导热系数。 O1-01 燃料电池电堆产热与水/热管理耦合 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角） O2-01 空压机冷却与进气中冷 等 FC 课里 λ_air ＝ 空气化学计量比（无量纲）；而在 C5-03 低温起雾机理与主动防雾控制 这类玻璃与遮挡课里，λ_air ＝ 静止空气的导热系数（≈0.026 W/(m·K)）。两者不在同一个板块里出现，但同一个仿真工程师可能同时读这两门课。⇒ 本课凡出现 λ_air 一律指化学计量比；导热系数在本课只在物性表里出现且一律写全称。

跨课提醒③：Q_coolant 的分账口径必须随数字一起交付。 O1-02 电堆冷却回路与低电导率冷却液设计 已经把这本账钉死：HHV 口径的 Q_stack 是 Q_coolant 的上界（一项都没扣，且成立还需要「回路侧加项之和 ＜ 扣项之和」这个前提）；LHV 口径那个工程近似既不是上界也不是下界。本课的建模动作是把那本账做成模型里真实存在的支路，而不是在后处理里减一个数——原因见第 2 讲。

数字上分三类，全课一致

• 物性与可复算的常数（F、M_H2O、L_f、E_tn 的显示值、乙二醇水溶液比热区间）给准确值并写明来源级别；
• 行业典型区间只给区间并标「典型」；
• 平台相关量（N_cell、V_cell、UA_s2c、(m·c_p)_stack、ΔT_coolant、ε_m、残留水量）本课给出典型参考带 ＋ 常见值锚点 ＋ 拿到它的路，用途是让你判断「自己手上那个数落不落在合理量级」。⚠ 它们是典型值、⛔ 不是验收依据，必须以本平台实测确认；
• 凡带代入值的算例，输入一律是教学假设值，⛔ 不对应任何平台、⛔ 禁止照抄取用。

所引标准本课未逐字核对原文，故只写标准号与它管什么对象，⛔ 不写条款内容与限值；版次以现行有效目录为准，引用前须核。

---

第 1 讲　电堆准稳态电化学—热子模型：降阶粒度怎么定、参数从哪来

1.1 先问模型要回答什么，再定分几个节点

「电堆在 1D 模型里分几个热节点」——这个问题被问的次数远多于它被正确回答的次数。常见的错误答法有两种：一种是「越细越准，分片级最好」；另一种是「1D 嘛，一个节点就够了」。两种都是在回答一个没被问清楚的问题。

正确的问法是：这一次仿真的输出量里，有没有出现只有更细的节点才能算出来的量？

按这条判据，粒度分成三档，逐档说清它给不了什么：

第 1 档：整堆集总单热节点。整堆一个热容、一个温度 T_stack、一个向冷却液的换热能力 UA_s2c。适用输出量：回路热负荷 Q_coolant、出堆冷却液温度、整堆升温 / 降温曲线、散热器与前端能力核算的输入。这一档在整车联合仿真里用得最多，也够用。

第 2 档：沿叠层方向分节点。把 N_cell 片按位置分成若干段（典型 3~7 段，端板段单独成段），每段一个热节点、一个电压。适用输出量：片间温差、最低片电压、端板段的升温滞后、冷启动时哪一段最先「掉队」。

第 3 档：一维根本给不了，必须转三维。只要输出量里出现这四类中的任何一类——片内脊-沟局部温度梯度、集管进出口的片间流量分配不均、端板边缘的二维冷点、局部液态水堵塞的位置——一维即失效。本课只给这条判据，方法本体见 J2-01 CFD 基础、湍流模型与网格划分（Fluent/STAR-CCM+）（网格与湍流）与 J2-04 共轭传热（CHT）仿真方法（固体—流体共轭传热），电堆内部固体热阻链见 O1-06 电堆内部传热路径与冷却流场板热设计。

⚠ 这条判据要按每一次的输出量重新过一遍，⛔ 不是给这个平台定死一档。同一台堆，做整车热平衡时第 1 档够，做冷启动成功率时必须第 2 档（原因见第 5 讲：P_succ 的判据里含 min_i V_cell,i，单节点模型根本没有这个量）。

1.2 产热源项：一律按电压尺算，⛔ 不许按电功率乘系数

是什么。电堆热源项只有一个正确写法：

Q_stack = N_cell · (E_tn − V_cell) · I

这个式子的含义在 O1-01 燃料电池电堆产热与水/热管理耦合 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角） 已经推透：把每一片的电压画成一根尺，E_tn 在上、E_rev 居中、V_cell 在下；(E_tn − E_rev) 那一截是可逆熵热，(E_rev − V_cell) 那一截是不可逆极化热，两截加起来正是 (E_tn − V_cell)。本课不重推，只负责把它正确落进模型。

⛔ 禁止的写法：Q_stack ≈ P_el × 某个系数。这种写法在 BEV 的电驱建模里可以接受（那里损耗确实近似正比于功率），在 FC 上会出两个问题：① 那个「系数」本身随工作点大幅变化（低负荷高电压点与额定点差得很远）；② 它把 HHV / LHV 的基准问题藏起来了——你填进去的系数到底是哪个基准算出来的，事后无从查证。

⚠ O1-01 燃料电池电堆产热与水/热管理耦合 给的 Q_stack ≈ P_stack·(1/η_stack − 1) 是估算式，不是建模式。它用来在方案阶段快速判量级（「散热量与电功率同量级」这个结论就是从它来的），⛔ 不用来做模型里的源项——因为 η 的口径与 E_tn 的口径必须成对，而估算式把这对关系藏进了一个数里。

HHV / LHV 基准怎么在模型里写死。三个动作，缺一不可：

1. 在模型参数表里把 E_tn 写成带基准后缀的两个独立参数：E_tn_HHV = 1.48、E_tn_LHV = 1.25（V/片，25 ℃ 显示值，承 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角）），并用一个显式开关选，⛔ 不许直接填一个裸的 1.48；
2. 凡是从供应商那里拿来的效率数、比功率数、热负荷数，入模前逐个标出它的基准；标不出的，退回去问，⛔ 不许猜；
3. 交付页第①项写明本模型用的是哪一端。

算例 1-A：同一个工作点，两套基准的产热差多少（教学假设值，⛔ 不对应任何平台）。

假设 N_cell ＝ 300 片，额定点 V_cell ＝ 0.65 V/片，I ＝ 400 A。

• HHV 尺：Q_stack(HHV) ＝ 300 × (1.48 − 0.65) × 400 ＝ 300 × 0.83 × 400 ＝ 99.6 kW
• LHV 尺：Q_stack(LHV) ＝ 300 × (1.25 − 0.65) × 400 ＝ 300 × 0.60 × 400 ＝ 72.0 kW
• 同点电功率：P_el ＝ 300 × 0.65 × 400 ＝ 78.0 kW

读出三件事： ① 两套基准差 (99.6 − 72.0)/99.6 ≈ 27.7%——这就是「差三成量级」的来处； ② Q_stack 与 P_el 同量级（99.6 vs 78.0）——这是 FCEV 散热难的第一性来源（承 O1-01 燃料电池电堆产热与水/热管理耦合 O3-01 FCEV 大散热量前端散热系统（低温差挑战）），也是本课第 2 讲所有取舍的前提； ③ 两套基准之差 99.6 − 72.0 ＝ 27.6 kW，恰是「生成水从液态记账改成气态记账」那一笔潜热。它不会凭空消失——第 2 讲会看到，它正是 Q_exh 里潜热那一项的来处。

算例 1-B：如果基准混用会算成什么样（教学假设值）。

某供应商给的规格书写「额定点系统效率 50%（LHV）」。工程师据此反推产热：P_el ＝ 78 kW ⇒ 输入化学能 ＝ 78/0.50 ＝ 156 kW ⇒ 产热 ＝ 156 − 78 ＝ 78 kW。然后他把这个 78 kW 填进模型，同时模型里 E_tn 用的是 HHV 的 1.48。

后果：模型标定时会发现算出的出堆温度偏低（因为真值 99.6 kW 被填成了 78 kW，少了 21.6 kW，约 22%）。为了对上台架实测的出堆温度，标定工程师唯一能动的旋钮是把 UA_s2c 调小——把堆冷得慢一点，温度就上去了。于是一个被算小 22% 的热源，和一个被调小的换热能力，在额定点上互相抵消，模型看起来完全正常。

⇒ 这就是引言里那个案例的第一处根因。它不可能被残差自检抓到（第 4 讲的能量守恒检的是账本内部闭不闭合，而这里账本是闭合的，只是整本账的量级错了），只能靠交付页第①项的声明 ＋ 第 6 讲的独立验证点（低负荷高电压点）暴露。

1.3 极化曲线怎么入模：三种做法与各自的失效边界

模型里必须有一个 V_cell ＝ f(I, T, p, λ_air) 的关系，否则 Q_stack 算不出来。它有三种实现，失效边界完全不同：

做法一：查表插值。直接把台架实测的 i-V 点做二维（或多维）插值。最稳、最少假设，因为它不引入任何函数形式。失效边界：外推区完全无效——插值算法在数据区外要么外推出物理上不可能的值，要么被钳在边界值上，两种都不报错。⇒ 用它的前提是工况矩阵完全落在数据区内，且必须在模型里加一个「输入超出表范围」的告警，⛔ 不许静默钳位。

做法二：分段半经验拟合（全工况模型的默认选择）。把曲线拆成三段，每段的参数有物理含义：活化段（对数形式）、欧姆线性段（斜率即面比电阻）、传质塌陷项（在极限电流附近陡降）。三段的机理与读法承 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角），本课只用它的结果。 它的好处是参数可以随温度、压力、λ_air 修正——而这正是热耦合模型必须的：温度一变，欧姆段斜率变（膜电导率随含水量与温度变）、传质段的极限电流变。查表做法要覆盖这些维度，表会大到不可标定。

做法三：纯多项式拟合。⛔ 禁止用于外推。高阶多项式在数据区内可以拟合得极好（R² 高到没话说），在数据区外会剧烈反向弯曲。⚠ 危险的地方在于：它在数据区内与半经验拟合几乎无法区分，你不会在标定阶段发现问题，只会在跑一个新工况时得到一个「电压比开路电压还高」或「电压变成负数」的结果。

算例 1-C：外推会错多少（教学假设值，⛔ 不对应任何平台）。

假设台架数据覆盖 0.2~1.2 A/cm²，单片有效面积 300 cm²，N_cell ＝ 300，E_tn(HHV) ＝ 1.48 V/片。三阶多项式是在 0.2~1.2 这段数据上最小二乘拟合出来的——它在区内与半经验几乎重合（差 < 0.005 V/片，远小于台架测量分散度）。

第一档：外推到 1.45 A/cm²（超出上界 21%），I ＝ 435 A。

• 半经验（含传质塌陷项）：V_cell ≈ 0.55 V/片 ⇒ Q_stack ＝ 300 × (1.48 − 0.55) × 435 ≈ 121.4 kW
• 三阶多项式：V_cell ≈ 0.60 V/片 ⇒ Q_stack ＝ 300 × (1.48 − 0.60) × 435 ≈ 114.8 kW ⇒ 差 6.6 kW ＝ 5.4%

第二档：外推到 1.80 A/cm²（超出上界 50%），I ＝ 540 A。

• 半经验：V_cell ≈ 0.33 V/片 ⇒ Q_stack ＝ 300 × (1.48 − 0.33) × 540 ≈ 186.3 kW
• 三阶多项式：V_cell ≈ 0.49 V/片 ⇒ Q_stack ＝ 300 × (1.48 − 0.49) × 540 ≈ 160.4 kW ⇒ 差 25.9 kW ＝ 13.9%

读出两件事： ① 误差随外推距离迅速放大（21% 外推差 5.4%，50% 外推差 13.9%），而且方向一律是把峰值产热算小——峰值点恰恰是散热能力核算最关心的那一点； ② ⚠ 标定阶段看不出来。两条曲线在数据区内的差别小于台架测量分散度，任何拟合优度指标都不会报警。⇒ 唯一的防线是「工况矩阵会不会跑出数据区」这个事前判断，⛔ 不是事后的拟合检查。

⚠ 三种做法共同的一条：曲线一律是 BOL（初期）口径。EOL 的电压下移会让同一电流下的产热抬升，这一档的定尺寸口径归 O1-07 电堆衰减与散热余量：从 BOL 到 EOL 的产热抬升与定尺寸。⚠ 该课目前是登记稿、大纲未结构化，读者点过去是空的；在其落盘前，本课的动作只有一个：在交付页第①项旁边注明「本模型 i-V 为 BOL 口径，未含衰减抬升」，⛔ 不代讲 EOL 的取值方法。

1.4 为什么电化学取准稳态，以及这个假设什么时候失效

本课的默认假设：在热耦合的时间尺度上，把 i-V 关系当代数式用。

依据是时间常数分离：双电层充放电在毫秒级、膜含水量重分布在秒级，而电堆的热时间常数（(m·c_p)_stack / 换热能力）在几十秒到分钟级。⇒ 当输出量是温度时，电化学侧的动态在热的尺度上已经完全松弛，用代数关系不会引入可观误差。

这个假设有一个明确的失效条件：输出量本身变成了「电压跌落」。

两种情形会撞上它：

• 冷启动初期——膜含水量、催化层冰堵、电导率都在剧烈变化，电压是被这些动态主宰的；
• 秒级强变载——空压机跟不上，λ_air 瞬时跌破目标，传质段提前塌陷（承 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角）），电压跌落的主因在气路而不在热。

⇒ 处置见第 5 讲。这里只记住判据：输出量是温度 ⇒ 准稳态可用；输出量是电压跌落 ⇒ 准稳态必须放弃。

⚠ 一个常见的误操作：为了「保险」，一上来就建全动态电化学模型。代价不只是算得慢——全动态模型引入的那批时间常数参数（双电层电容、膜水扩散时间常数）在你的实测数据里没有对应的观测量，标定时它们只能被优化器随便分配。⇒ 这是「分得太细却没有对应观测量」的典型，模型在标定工况上照样能拟合，误差只是在参数之间乱窜。

1.5 热容必须分列，⛔ 不许并成一个数

是什么。模型里的 (m·c_p)_stack 不能是一个数，至少要拆成两块：

分块	包含	为什么单列
片间热容	MEA ＋ 双极板 ＋ 片内滞留冷却液	与冷却液直接换热，升温最快
端部热容	端板 ＋ 集流板 ＋ 拉杆 ＋ 绝缘板	升温最慢的一段，且与冷却液的换热路径长

为什么这条是硬要求。端板是整堆升温最慢的部分（端部偏冷的机理承 O1-01 燃料电池电堆产热与水/热管理耦合 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角））。如果把端板热容并进整堆平均热容，模型算出来的是一条整体均匀升温的曲线——它会让冷启动时间算短、成功率算高。

⇒ 误差方向偏危险。这一点值得单独强调：本课绝大多数建模简化的误差方向是偏保守（算得更热、更难），只有这一条相反。偏保守的简化在设计上只是浪费，偏危险的简化会让一个实际过不了的平台在模型上通过。

算例 1-D：并成一个数把冷启动时间算短多少（教学假设值，⛔ 不对应任何平台）。

假设：片间有效热容 25 kJ/K、端部有效热容 8 kJ/K；自加热注入功率 12 kW，其中约 80% 落在片间、20% 传到端部（简化的两节点分配）。从 −30 ℃ 升到 0 ℃，需升温 30 K。

并成一个数的算法：总热容 ＝ 25 + 8 ＝ 33 kJ/K；升温时间 t ＝ 33 × 30 / 12 ＝ 82.5 s。

分列的算法：

• 片间：t ＝ 25 × 30 / (12 × 0.8) ＝ 25 × 30 / 9.6 ＝ 78.1 s
• 端部：t ＝ 8 × 30 / (12 × 0.2) ＝ 8 × 30 / 2.4 ＝ 100 s

⇒ 整堆真正「全部越过 0 ℃」的时间是 100 s，而并成一个数给的是 82.5 s，算短了 17.5%。若门限 t_lim 恰好卡在 90 s，这个模型会把一次实际失败的启动判成成功。

⚠ 上面那个 80/20 的分配比是教学假设，真实分配取决于自加热的注入方式（承 O1-04 低温冷启动（-30℃）与自加热策略）与堆的结构。⇒ 拿到它的路：在环境舱里做一次带端板测点的自加热试验，用两个温升斜率反解。

1.6 参数从哪来，以及「拿不到数据时怎么办」

本讲需要的参数只有五个，逐个给来源与典型带：

参数	符号	典型参考带	常见值锚点	拿到它的路
片数	N_cell	200~450 片	乘用车 80~120 kW 级多在 300 片附近	电堆规格书（确定值）
单片额定电压	V_cell	0.60~0.72 V/片	额定点多在 0.65 V/片附近	i-V 实测曲线
片间有效热容	(m·c_p)_cell,tot	15~45 kJ/K	80~120 kW 级多在 25 kJ/K附近	① 按 BOM 逐件累加；② 由 soak 实测反解
端部有效热容	(m·c_p)_end	4~15 kJ/K	多在 8 kJ/K附近	按端板 ＋ 集流板 ＋ 拉杆的 BOM 累加
堆侧换热能力	UA_s2c	3~15 kW/K	额定流量下多在 8 kW/K附近	台架上扫流量，用 ε_s2c 反解（见第 2 讲）

⚠ 以上五项中后三项是平台相关量：区间与锚点是典型值，⛔ 不是验收依据，必须以本平台实测确认。它们的用途只有一个——让你判断自己标出来的数落不落在合理量级。若你标出来的 UA_s2c 是 40 kW/K 或 0.5 kW/K，那不是「这个平台特殊」，那是标定过程里有别的错误被这个参数吸收了。

⚠ 另注意 (m·c_p) 的口径：O1-05 停机吹扫、保温与整堆防冻设计 给的 (m·c_p)_stack ＝ 20~60 kJ/K 是含堆内滞留冷却液的整堆值。本课把它拆成片间 ＋ 端部两项，两项之和应落在 O1-05 停机吹扫、保温与整堆防冻设计 那个带内——这是一个免费的量级自检，⛔ 不要跳过。

拿不到逐片数据时怎么办：唯一正确的动作是降档并声明。

如果供应商只给整堆的 i-V 与整堆热容，不给逐片分布，那么你只能建整堆单节点模型，并在交付页第②项写明「本模型不评片间温差与最低片电压」。

⛔ 不许用假设的分布凑出多节点模型。比如「假设片间温差呈抛物线分布，端部两片各低 5 K」——这种做法会造出一个看起来能给出片间温差的模型，而那个温差完全是假设的产物。后果在第 5 讲会兑现：P_succ 的判据里含 min_i V_cell,i，用假设分布算出来的最低片电压会得到一个看不出错的假成功率。

算例 1-E：整堆热容的量级自检（教学假设值）。

按上表锚点：片间 25 kJ/K ＋ 端部 8 kJ/K ＝ 33 kJ/K。落在 O1-05 停机吹扫、保温与整堆防冻设计 的 20~60 kJ/K 带内、且接近它给的 30 kJ/K 锚点 ⇒ 量级自检通过。

再用 τ 做第二次自检：若这台堆在环境舱里做纯 soak（不加热、不通气、关阀）实测出的时间常数 τ ≈ 5500 s，则 UA_amb ＝ 33000 / 5500 ＝ 6.0 W/K——落在 O1-05 停机吹扫、保温与整堆防冻设计 给的「做了保温且关阀多在 5~8 W/K」带内 ⇒ 通过。 ⚠ 注意这里算出来的是 UA_amb（W/K 级），⛔ 不是 UA_s2c（kW/K 级）——它们差三个数量级，见引言的跨课提醒①。

1.7 本讲不重讲什么

• 三段极化的机理与曲线读法——前置 O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角），本讲只讲它怎么入模；
• 可逆熵热与不可逆极化热的拆分、膜内水传递三路径——前置 O1-01 燃料电池电堆产热与水/热管理耦合，本讲只讲它们怎么降成源项与状态量；
• 催化层—GDL—双极板的固体侧热阻链与冷却腔构型——O1-06 电堆内部传热路径与冷却流场板热设计，本讲把整条链压成一个 UA_s2c；
• HHV / LHV 两套口径的换算本身——O1-08 PEMFC 工作原理与电堆构造导论：极化曲线、Nernst 方程与两套效率口径（热管理视角） 是唯一归属，本讲只负责声明基准。

1.8 本讲小结

粒度由输出量倒推，三档各写清「给不了什么」；产热一律走电压尺、⛔ 不走电功率乘系数，且 E_tn 必须带基准；极化曲线三种入模做法差别全在外推区，全工况默认选分段半经验；电化学取准稳态的前提是「输出量是温度」，一旦输出量变成电压跌落就必须放弃；热容分列是硬要求，因为并成一个数的误差方向偏危险；拿不到逐片数据就老实降档并在交付页声明，⛔ 不许用假设的分布凑节点。