本研究报告介绍了一项关于质子交换膜（Proton Exchange Membrane， PEM）燃料电池动态建模与验证的原创性研究。该研究由Caisheng Wang（学生会员）、M. Hashem Nehrir（高级会员）和Steven R. Shaw（会员）共同完成，他们均来自美国蒙大拿州立大学电气与计算机工程系。该研究成果以论文形式《Dynamic Models and Model Validation for PEM Fuel Cells Using Electrical Circuits》发表于电气电子工程师学会（IEEE）旗下的期刊《IEEE Transactions on Energy Conversion》2005年6月刊，第20卷第2期。
一、 学术背景与研究目标
本研究隶属于燃料电池系统建模与控制领域，具体聚焦于质子交换膜燃料电池的动态特性。PEM燃料电池作为一种分布式发电（Distributed Generation， DG）技术，因其不受地理限制、排放低、功率密度高等优点，在分布式发电和电动汽车领域展现出巨大潜力。然而，与风电、光伏等其他DG技术相比，PEM燃料电池对电负载瞬态变化的响应速度较慢，这主要源于其内部缓慢的电化学与热力学响应过程。这种动态特性限制了其在应对电机启动、电动汽车加速或系统故障等瞬态工况时的性能。因此，深入理解并准确预测PEM燃料电池的动态行为，对于设计有效的控制系统、优化系统集成和提升整体性能至关重要。
在本文发表之前，已有文献报道了燃料电池的稳态模型和一些动态模型。这些模型大多基于经验公式或电化学反应原理。近年来，一些研究开始关注燃料电池电学端口特性的动态建模。文献[17]详细解释了燃料电池的电化学特性，并提出了一个包含双电层充电效应电容的简单等效电路。然而，一个能够全面集成电化学、物质传输、热力学动态过程，并便于与电力电子设备联合仿真的综合性动态模型，仍然是当时研究的重点和难点。
基于此背景，本研究旨在开发一种基于电路原理的PEM燃料电池动态模型。其核心目标包括：1) 通过机理分析，建立包含活化损失、欧姆电压降、浓度电压降的等效内部电压源和等效电阻；2) 在模型中整合双电层充电效应和燃料电池内部的热力学特性；3) 在MATLAB/Simulink和PSpice两种仿真环境中分别实现该动态模型；4) 利用真实的500W PEM燃料电池堆（Avista Labs SR-12）的实验数据，对所建模型的稳态和瞬态响应进行验证；5) 最终提供一个可用于燃料电池实时控制及多源电力系统应用研究的有效工具。
二、 详细研究流程
本研究的工作流程可分为以下几个关键步骤：
第一步：模型机理分析与数学建模。 研究首先基于一系列合理的假设（如一维处理、理想气体、流道压力恒定、工作温度低于100℃等），从燃料电池内部物理化学过程出发，建立了完整的数学模型。这个流程包含了多个子流程： 1. 气体扩散过程建模： 基于Stefan-Maxwell公式，推导了阳极和阴极多孔电极内氢气（H₂）和氧气（O₂）的有效分压计算公式（公式(1)-(13)）。这涉及将复杂的扩散过程简化为一维传输，并应用法拉第定律关联电流与反应物消耗速率，最终得到有效分压与电流、流道压力、扩散层厚度等参数的函数关系。 2. 物质守恒方程建立： 利用理想气体状态方程，建立了氢气和氧气有效分压随时间变化的动态方程（公式(14)-(15)）。在稳态下，分压变化率为零；在瞬态下，考虑了燃料和氧化剂流量变化相对于负载电流变化的延迟效应，并用一阶延迟模型来描述（公式(18)）。 3. 输出电压计算模型构建： 这是模型的核心。首先，利用能斯特方程（Nernst equation，公式(20)）计算燃料电池的可逆电势（内部电势），该电势是温度、氢气分压和氧气分压的函数。同时，将燃料/氧化剂流延迟的影响建模为一个从内部电势中减去的电压项V_d（公式(22)-(23)）。然后，从内部电势中依次减去三种主要的电压损失：a) 活化电压损失(V_act)：采用Tafel方程（公式(27)）和一个经验修正项（公式(28)）进行计算，并将其等效为一个与电流和温度相关的电阻R_act（公式(29)）。b) 欧姆电压损失(V_ohm)：由膜电阻、接触电阻等引起，表示为电流和温度的函数（公式(30)-(31)），等效为电阻R_ohm。c) 浓度电压损失(V_conc)：由高电流密度下反应物传输限制引起，根据菲克定律和法拉第定律推导得出（公式(33)），等效为电阻R_conc（公式(34)）。最终，单电池输出电压V_out = E - V_act - V_ohm - V_conc，而电堆电压为单电池电压乘以电堆串联电池数（公式(25)-(26)）。 4. 双电层充电效应集成： 研究考虑了PEM燃料电池中多孔电极与膜界面形成的双电层（Double-Layer）的电容效应。该效应被建模为一个跨接在活化电阻和浓度电阻两端的大电容C（数量级可达数法拉），如图2所示。这使得流经活化/浓度支路的电流分为两部分：一部分用于电化学反应，另一部分用于对双电层电容充电。因此，计算输出电压时，用于产生活化/浓度损失的电压是基于电容两端电压V_C（公式(35)）而非电流I，最终输出电压公式更新为公式(36)。 5. 热力学能量平衡集成： 建立了燃料电池堆的温度动态模型。通过能量平衡方程（公式(37)），计算化学反应产生的净热量。该净热量等于可用化学能（基于吉布斯自由能变化，公式(38)-(39)）减去电输出功率（公式(40)），再减去反应物/产物带来的显热和潜热（公式(41)）以及通过对流散失到环境的热量（公式(42)）。净热量导致电堆温度变化，其动态关系由公式(43)描述，其中包含了电堆总质量和总比热容。
第二步：模型在MATLAB/Simulink中的实现。 基于上述数学模型，研究团队在MATLAB/Simulink环境中构建了动态模型。如图3所示，该模型以阳极/阴极压力、初始温度、环境温度为输入，负载电流和时间作为关键变量。模型内部通过一系列相互关联的功能模块（子模块）来计算气体有效分压、内部电势、各项电压损失（考虑双电层电容动态）、并求解能量平衡方程以获得实时温度。温度与电流作为反馈信号，参与内部电势和各项损失的计算，从而形成一个闭环的动态系统，最终输出燃料电池的端电压。
第三步：模型在PSpice中的等效电路实现。 为了便于燃料电池与电力电子变换器等外部电气设备进行联合仿真，研究团队进一步在PSpice中开发了等效电路模型。该模型的总体框图如图4所示，其创新之处在于将所有的物理和化学过程都用等效的电气元件和受控源来表征。 1. 内部电势电路： 如图5所示，使用电流和温度控制的电压源来模拟能斯特方程及燃料延迟电压V_d。 2. 活化损失电路： 如图6所示，将活化电压损失V_act分解为与温度相关的常数部分V_a0和与电流、温度相关的电阻压降部分。后者通过一个由常数电阻、电流控制电阻和温度控制电阻串联组成的等效电阻R_act来模拟。PSpice中的多项式电流控制电压源（HPOLY）和模拟行为建模（ABM）功能被用来实现这些受控电阻。 3. 欧姆与浓度损失电路： 采用与活化损失电路类似的方法，分别构建了由常数、电流相关、温度相关电阻部分串联而成的等效电阻R_ohm（图7）和R_conc（图8）。 4. 热力学模块电路： 这是本模型的一个亮点。研究利用热电类比（Thermal-Electrical Analogy）原理，将热学量（热流、温差、热容、热阻）映射为电学量（电流、电压、电容、电阻），如表I所示。基于此，将能量平衡方程（公式(37)-(43)）转换成了如图9所示的RC电路。其中，电压源代表环境温度，电容C_th上的电压代表电堆温度，电阻R_conv代表对流热阻，而电流源代表内部产生的净热功率（即电损失功率I*(V_act + V_ohm + V_conc)）。
第四步：模型参数确定与实验验证。 研究使用Avista Labs SR-12型500W PEM燃料电池堆的规格参数（表II）作为模型输入。部分模型参数（如双电层电容C、对流换热系数h等）通过文献或实验测量确定（表III）。为了验证模型，搭建了如图10所示的实验平台：使用可编程电子负载（Chroma 63112）施加负载，采用LEM公司的电流和电压传感器采集信号，使用K型热电偶测量温度，并通过数据采集卡将数据录入PC。验证实验分为两部分： 1. 稳态特性验证： 将负载电流从1.1A缓慢阶梯增加至20.5A（每步0.2A，持续40秒），测量并记录稳态的电压-电流（V-I）曲线和功率-电流（P-I）曲线，同时记录温度变化。 2. 瞬态特性验证： 对负载进行快速阶跃变化，测量燃料电池输出电压的瞬态响应波形，以考察模型对短时间常数（双电层效应）和长时间常数（流量延迟、热力学响应）动态过程的预测能力。
三、 主要研究结果
在稳态特性验证方面： 图11展示了SR-12电堆实验测得的平均V-I曲线与Simulink、PSpice模型仿真结果的对比。两者吻合度非常高。曲线清晰地显示出三个区域：起始段的快速电压下降主要由活化损失导致；中间段的近似线性下降主要由欧姆损失主导；末端在高电流密度下的急剧电压下降则归因于浓度损失。图12的P-I曲线对比也显示模型能准确预测最大功率点及其后的功率下降趋势。图13的温度响应对比显示，模型预测的温度变化趋势与实验数据一致，在最差情况下的偏差在3.5K以内，证明了集成热力学模型的有效性。
在瞬态特性验证方面： 模型成功捕捉了PEM燃料电池在不同时间尺度上的动态特性。图14展示了负载电流快速阶跃时，模型在短时间范围（秒级）内的响应。输出电压在阶跃瞬间因欧姆电阻压降而骤降，随后由于双电层电容的充放电效应，电压呈指数形式缓慢衰减至新的稳态值。这一过程的时间常数τ = R_act * C，由于R_act较小（本研究中小于0.05欧姆），即使电容C很大（数法拉），时间常数仍小于1秒。图15则展示了在长时间范围（数十秒至数分钟）内的瞬态响应。当负载阶跃上升后，输出电压会先降至一个较低值，然后随着燃料/氧化剂流量的缓慢调整以及电堆温度的逐渐变化，电压缓慢恢复至稳态值。这一慢动态主要由公式(18)描述的流量延迟和公式(43)描述的热惯性所支配。Simulink和PSpice模型的仿真结果与实验测量数据均吻合良好，证明了所建模型能够全面、准确地模拟PEM燃料电池从毫秒级到分钟级的动态行为。
这些结果逻辑严密地支撑了研究结论： 首先，稳态V-I/P-I曲线的精确匹配，证明了模型在静态工作点参数计算（如分压、能斯特电势、各项损失）方面的正确性。其次，短时瞬态响应的吻合，验证了双电层电容效应建模的准确性。最后，长时瞬态响应的成功预测，综合验证了物质传输延迟模型和热力学动态模型的正确性。这三方面的验证结果共同表明，所提出的基于电路的动态模型是一个高度可靠的仿真工具。
四、 研究结论与价值
本研究成功开发并验证了用于PEM燃料电池的动态电路模型。主要结论是：所提出的模型能够精确预测PEM燃料电池堆在稳态和各种瞬态工况下的电气响应，并能有效模拟电堆的温度动态变化。模型在MATLAB/Simulink和PSpice两种平台上的实现，为不同应用场景提供了便利。
本研究的价值体现在以下几个方面： 1. 科学价值： 它提供了一个从电气工程师视角出发的、统一的建模框架，将燃料电池内部复杂的电化学、传质、传热过程巧妙地映射为易于理解和分析的电路元件（电压源、电阻、电容、受控源）和热电类比电路。这深化了对PEM燃料电池多物理场耦合动态行为的机理认识。 2. 应用价值： 模型具有重要的工程实用价值。MATLAB/Simulink模型适用于控制系统设计和算法开发；而PSpice等效电路模型则便于将燃料电池模型直接嵌入到包含电力电子变换器、电机负载等在内的完整电气系统电路中进行联合仿真，这对于评估燃料电池在混合动力系统、分布式发电系统中的动态性能至关重要。 3. 工具价值： 验证后的模型可作为后续相关研究的基础工具，例如用于燃料电池的实时控制器设计、故障诊断、寿命预测以及多能源系统（如燃料电池-蓄电池混合系统）的优化运行研究。
五、 研究亮点
六、 其他有价值的内容
研究中详细推导的气体扩散方程和物质守恒方程，为理解反应物在电极内的传输过程及其对输出电压的影响提供了清晰的数学描述。对活化、欧姆、浓度三种电压损失的机理分析与等效电阻公式，有助于定量分析不同运行条件（电流、温度）下各种损失的比例变化。此外，论文中提供的完整符号表（Nomenclature）和详细的模型参数列表（表III），为其他研究者复现和应用该模型提供了极大的便利。最后，研究中对实验配置（图10）的说明，也为从事燃料电池实验验证的研究人员提供了有价值的参考。