·技术交流·
燃料电池汽车空气系统的仿真与验证①
刘冬安,陈俊超,蒋永伟
(上海汽车集团股份有限公司商用车技术中心,上海 200438)
摘要:系统仿真作为现代科学研究方法的一个重要分支,对于系统的整体性能及其细节表现的研究
具有重要意义。燃料电池系统包含电堆、空压机、中冷器等,涉及到的零部件众多。通过系统仿真
手段对其进行前期的设计计算和后期的校核计算可以缩短燃料电池系统的研发周期,进而节约研
发成本,使得该手段在燃料电池汽车的研发方面具有巨大优势。利用Simcenter Amesim软件对
某款燃料电池汽车的燃料电池空气系统进行了系统匹配仿真分析,并且引入了机器学习中对于回
归分析的评价指标———平均绝对误差和均方误差对仿真结果进行评价,验证了仿真模型的准确性
与可行性。结果表明:平均绝对误差和均方误差均保持在4.40%以内,仿真结果与实测报文数据吻
合度较高。此仿真模型能较为真实地反映质子交换膜燃料电池的工作特性,能有效地用于燃料电
池系统性能的研究。
关键词:系统仿真;建模;Simcenter Amesim;燃料电池;空气系统
中图分类号:TM 911.4  文献标识码:A  文章编号:1008-7923(2020)03-0115-06
Simulation and Validation of Fuel Cell Vehicle Air Supply System
LIU Dong-an,CHEN Jun-chao,JIANG Yong-wei
(Shanghai Automotive Industry Corporation Commercial Vehicle Technical Center,Shanghai,200438,China.)
Abstract:The system simulation serves as a crucial branch of modern science research methodolo-
gy,which is of great significance to the performance of the whole system as well as its details.The
fuel cell system involves lots of components such as stack,air compressor,intercooler and so on.
By using the system simulation method can shorten the development period during the design
process,thus reducing the cost of the research,which shows great advantages in fuel cell vehicle s
research and development.An Simcenter Amesim model of the air supply system for one fuel cell
vehicle is developed and introduced mean absolute error and mean square error used in machine
learning to verify the model s accuracy and feasibility.The results indicate that both mean abso-
lute error and mean square error are less than 4.40%,which shows the simulation results are in
good agreement with the actual data.This model may reflect the proton exchange membrane fuel
cell s performance correctly,so it can be applied to the research of fuel cell systems effectively.
Keywords:System simulation;Modeling;Simcenter Amesim;Fuel cell;Air supply system
①作者简介:
刘冬安(1981-),男,河北人,高级工程师,博士。研究方向:燃料电池及储氢系统技术。联系地址:上海市杨浦区军工路2500号,邮编:200438。Email:LiuDongAn@saicmotor.com;Tel.:021-60569999-4301.
1 引言
日益严峻的环境污染和温室效应,对于能源的合理利用提出了全新的要求。作为全人类社会提供了一种全新的能源选择方案,燃料电池作为继水力、火力、核电之后的第四代发电技术,引起了世界各国政府和学者的大力重视。其中,氢燃料电池相对其他类型的燃料电池在移动出行领域方面,优势尤为明显。和传统汽车发动机能量转化效率普遍低于40%的情况相比,氢燃料电池的平均能量转化效率大概在47%左右甚至更高。此外,氢燃料电池的运行噪声较低,热辐射较低。相比二者的排放物,汽车发动机反应后产生的二氧化碳、氮氧化物、硫氧化物等均会对环境造成污染,而氢燃料电池的排放物仅仅为水,做到了真正意义上的零污染。鉴于上述优点,氢燃料电池汽车具有极其光明的应用前景[1-4]。而对比纯电动汽车的温度耐受面窄和续航里程焦虑等发展瓶颈,氢燃料电池系统的温度耐受范围
为-20℃~80℃,且续航里程基本上可以大于400km。很多专家学者都认为新能源汽车发展的最终目标是氢燃料电池汽车[5]。
燃料电池系统又名燃料电池发动机。燃料电池系统由空气系统、水热管理系统、电气系统和控制系统等组成,与内燃机的燃料供给系统、润滑系统、冷却系统、点火系统和起动系统有一定的相似之处。燃料电池的空气系统处于燃料电池堆的阴极侧。为了保证燃料电池堆电化学反应的顺利进行,需要对燃料电池堆的阴极源源不断地供应符合反应条件的空气流量。其中,燃料电池空气系统中的空压机及其控制器更是空气系统的核心零部件,对于燃料电池系统的正常运行起着极其重要的作用。空气系统的控制本质上属于负反馈控制,为此必须对空气的温度、压力和进入燃料电池堆的空气流量采用相应的传感器进行实时监测。传感器将读取到的信号传入电子控制器中,控制器采取相应的动作来控制被控元件[6]。对于车载工况,常采用(PEMFC:Proton Exchange Mem-brane Fuel Cell)质子交换膜燃料电池作为燃料电池汽车的动力源。
本文研究了燃料电池汽车空气系统的架构及其工作原理以及各子模块的数学模型,并且利用Sim-center Amesim系统仿真软件搭建了PEMFC燃料电池空气系统的仿真模型,并引入了回归分析中的平均绝对误差和均方误差两大评价指标对仿真结果进行评价,为燃料电池系统仿真的研究工作提供了一定的参考意义。2 燃料电池汽车空气系统的架构及其工作原理2.1 燃料电池汽车空气系统架构
空气系统主要由ACU(Air Compressor ControlUnit)空压机控制器、空压机、空滤、空气流量计、中冷器、膜增湿器、燃料电池电堆、冷凝器、蝶阀、温压一体传感器以及相应的空气管路等零部件组成。以上零部件共同工作,为燃料电池堆在适当的工况下提供相应流量的空气。燃料电池汽车空气系统的架构如图1所示。
图1 燃料电池汽车空气系统架构图
Fig.1 The block diagram of the fuel cell vehicle’s air supply system.
2.2 燃料电池车空气系统的工作原理
空气系统是由空压机控制器ACU(Air Com-pressor Control Unit)、整车控制器VCU(VehicleControl Unit)和燃料电池系统控制器FCU(FuelCell Control Unit)三者共同控制的。VCU根据整车的实际运行情况向FCU发出功率请求后,燃料电池系统的FCU向ACU发出相应的空气流量信号,ACU接收到空气流量信号后控制空压机的转速,从而控制进入电堆中空气的流量。FCU将电堆发出的实时功率监测值和VCU发出的请求功率值进行对比后,进行负反馈调节,空压机控制器ACU从而调节空压机转速,进而调节进入燃料电池电堆内部的空气的流量。空气系统尾部的背压是通过FCU对蝶阀的开度进
行控制来实现的,FCU与空压机控制器共同调节电堆内部空气的压力,保证燃料电池电化学反应的顺利进行。
3 空气系统的Amesim建模
燃料电池空气系统的主要零部件有空压机、中冷器、膜增湿器、冷凝器、蝶阀、管路、传感器和控制器等。仿真系统模型中还包括燃料电池电堆模型,目的是为了精确模拟空气系统对电堆性能的影响,仿真空气系统在各个工况条件下的表现。以上零部件共同构成了本文的Amesim模型。
3.1 空压机的Amesim建模
由于电堆模块的阴极侧的空气供给需要一定的空气压力,故需要采用空压机对空气进行压缩。本模
型采用双螺杆式空压机。空压机的Amesim模型如
图2所示。
图2 空压机的Amesim模型
Fig.2 The Amesim model of the air comp
ressor.Amesim中的空压机模型需要根据空压机特性
曲线输入两张对应的表格Table 1和Table 2,作为修正空气流量dmc和修正空压机转速wc的查询依据。空压机参数中需要输入的绝对压力比Pr、等熵效率ηis等,
均可由这两张表格来进行查询。对于绝对压力比[7]
,可以表示为:
Pr=Pdown
Pup
(1
)式中,Pdown———下游空气压力;Pup———上游空气压力。
对于等熵效率[8]
,可以表示为:
ηis=Pis
real
(2
)式中,Pis———等熵功率;Preal———实际功率。默认情况下,修正空气流量dmc的数学模型为:
dmc=dm·TupTst槡
·Pst
Pup
(3
)式中,dm—
——空压机中的实际空气流量;Tup—
—上游空气温度;Tst———标准空气温度;Pup———上游空气绝对压力;Pst———标准空气绝对压力。
默认情况下,修正空压机转速wc的数学模型为:
wc=w·
Tst
Tup槡
(4
)式中,w—
——空压机实际转速。当实际的修正空气流量dmc高于或低于默认的修正空气流量dmc的数值时,Amesim会根据
Table1表中的数值,
采用线性插值或多项式插值等方式,计算出实际的绝对压力比Pr和等熵效率ηis,进而协助Ame
sim软件完成整个空压机的数学模型的计算。
3.2 中冷器的Amesim建模
空气经过空压机的高速旋转压缩后,温度会上升至80℃以上。过高的空气温度会影响质子交换膜的寿命。故需要通过中冷器对压缩后的空气进行降温至大约60℃左右。中冷器的Amesim模型如图3所示。
图3 中冷器的Amesim模型
Fig
.3 The Amesim model of the inter-cooler.中冷器本质上属于换热器的一种,Amesim的中冷器采用效能-传热单元数法(ε-NTU)
来进行计算。效能ε的定义[9]
为:ε=(t′-t″
)max
t′1-t′
(5
)式(5)中,t′1———热流体的进口温度;t′
2—
——冷流体的进口温度;t
———热流体或冷流体的出口温度。该式中的分母为流体在换热器中可能发生的最大温差值,而分子则为冷流体或热流体在换热器中的实际温度差值中的大者。
传热单元数的定义[10]
为:
NTU=
kA
(qmc)min(6
)式中,k———平均传热系数;A—
——中冷器的有效换热面积;qm———流体的质量流量;c———流体的比热容。
对于空气系统的中冷器来说,属于两种流体各自均非混合的单流程叉流流动。德雷克提出了一个近
似关系[
11]:ε=1-expNTU0.22
[exp(-C*NTU0.78
)-1]{}
(7
式(7)中,C*称为热容比,
其定义为:C*=WminWmax=(qmc)min
(qmc)max
≤1
(8)3.3 膜增湿器的Amesim建模
质子交换膜需要有一定的含水量,PEMFC才会达到较大的反应速率。故需要利用膜增湿器对进堆前的空气进行增湿,使得空气湿度达到60%左右。膜增湿器的Amesim模型如图4所示。
图4 膜增湿器的Amesim模型Fig
.4 The Amesim model of the humidifier.膜增湿器在Am
esim中是预先设定好的一种超级零部件。所谓超级零部件,即许多子元件的组合。对于膜增湿器来说,本质上是由一个超级零部件。该超级零部件由湿度传感器、温度传感器、空气容腔、信号控制模型和信号传递模型组合而成。通过将一系列相关的子零件封装起来,隐藏内部具体的实现功能,
燃料电池电动汽车
仅仅通过超级原件对外展示必要的接口信息,是系统仿真常用的方法。膜增湿器超级零部件的Amesim模型如图5所示。
图5 膜增湿器超级零部件的Amesim子模型
Fig.5 The Amesim submodel of the humidifier super comp
onent.膜增湿器的数学模型[12,13]
如下式(9
)所示:lg10(Psat)=32-3 
152.2T
-7.3×lg10(
T)+2.4×1
0-9 T+1.8×10-6 T2
(9)膜增湿器中相对湿度定义为湿空气中水蒸气分
压力pvap
or与同一温度同样总压力的饱和湿空气中水蒸气分压力psat的比值。其计算公式[14]
如下式(10)所示:
φ=
pvapor
psat
(10)3.4 冷凝器的Amesim建模
空气经过燃料电池电堆反应后形成具有一定湿度的湿空气。湿空气为干空气和水蒸气的混合物。为了保
证蝶阀对空气背压调节的效率,需要反应后的空气在通入蝶阀之前将其进行冷凝处理,保证进入蝶阀的空气尽可能接近100%干空气。故需要采
用冷凝器将湿空气进行冷凝,湿空气经过冷凝器处理后变为干空气和100%冷凝水两部分,冷凝水通过尾排管排到车外。冷凝器的Amesim模型如图6所示。
图6 冷凝器的Amesim模型
Fig
.6 The Amesim model of the condenser.空气在冷凝器中发生对流换热过程,该物理过程
遵循牛顿冷却定律。冷凝器的数学模型[
15]
如下式(11
)所示:Q=hA(tf-
tw)(11)式中,Q———热流量,单位为W;h———对流换热系数,单位为W/(m2·K);A———有效换热面积,单位为m2;tf—
——流体的温度,单位为K;tw———冷凝器壁面的温度,单位为K。3.5 蝶阀的Amesim模型
空气进入燃料电池电堆反应以后,其中一部分反应后的空气重新返回膜增湿器对刚刚进入膜增湿器中的新空气进行增湿,另一部分反应后的空气则通过蝶阀排出至大气中。蝶阀通过控制其开度的大小,和
空压机一起控制进入电堆中空气的压力,保证膜电极两侧的空气和氢气的压力差保持在合理的范围之内。蝶阀的Amesim模型如图7所示。
蝶阀阀片旋转角度的变化范围为[0,90°]。在Amesim中输入蝶阀最大内径、
阀片的直径等参数以
图7 蝶阀的Amesim模型
Fig.7 The Amesim model of the throttle valve.
外,还需要外界对蝶阀进行相应的控制,将控制信号传递至蝶阀,蝶阀可以控制阀片的旋转角度,进而控制空气排出电堆的流量,调节整个空气系统的背压,从而保持电堆阴极侧空气压力在合理的范围之内。对于燃料电池系统的空气系统来说,可以根据电堆发出的电流需求来控制蝶阀的开度[16]。
3.6 整个空气系统的Amesim建模
燃料电池车的空气系统是由各个零部件通过一定的机械连接和电气连接相连而成。为了将更多关注点聚焦于空气系统本身,需采用一定的方式来简化整个燃料电池系统的仿真。采用恒压H
源来近似代替储氢气瓶,且电堆在不同工况下需要的散热量采用电堆在各个工作点的能量转换效率通过折算和插值得出。整个燃料电池车空气系统的Amesim模型如图8所示。
图8 整个空气系统的Amesim模型
Fig.8 The Amesim model of the air supply system.
整个空气系统仿真需要的主要参数如表1所示。
表1 整个空气系统仿真需要的主要参数
Table 1 The main required parameters of the whole air supply system
used for simulation.
参数值
空压机额定转速/rpm 8 800
阴极氧气初始体积分数/%21
阴极氮气初始体积分数/%78.999
阴极水蒸气初始体积分数/%0.0014 空气系统的仿真模型验证
系统仿真作为现代工程研发的一个重要手段,其应用常常表现在两个方面:在研发前期,对设计方案通过系统仿真进行模拟分析,评估方案的可行性,属于设计计算的内容;在研发后期,通过对设计好的系统进行模拟仿真,并与实际工况的测试数据进行对比,验证仿真模型的准确性,以便仿真模型和仿真经
验可以推广到未来的全新研发项目中,属于校核计算的内容。本文的内容属于研发后期的校核计算工作。4.1 仿真的评价指标
机器学习算法是人工智能中的一个重要分支。由于系统仿真本质上可以看作一种回归分析,故可以借鉴机器学习算法中对于回归分析的两个常用评价指标———平均绝对误差(MAE:Mean Absolute Error)和均方误差(MSE:Mean Square Error)来进行评价。MAE和MSE的公式详见式(12)和式(13)。
MAE=
ni=1|yi-yi
^|(12)MSE=
ni=1(yi-yi
^)2(13)4.2 空气系统的仿真模型验证
VCU根据整车的实际运行情况向FCU发出功率请求后,燃料电池系统的FCU向ACU发出相应的空气流量信号,ACU接收到空气流量信号后控制空压机的转速,从而控制进入电堆中空气的流量。FCU将电堆发出的实时功率监测值和VCU发出的请求功率值进行对比后,进行负反馈调节,空压机控制器ACU从而调节空压机转速,进而调节进入燃料电池电堆内部的空气的流量。空气系统尾部的背压是通过FCU对蝶阀的开度进行控制来实现的,FCU与空压机控制器共同调节电堆内部空气的压力,保证燃料电池电化学反应的顺利进行。
燃料电池系统的空气流量计采集反应所需的空气流量,从而将其传递至CAN网络上,可以通过实车报文采集方式将其采集。进而进行空气流量变化趋势的分析。图9为空气流量的仿真结果与实车报文的对比。如图所示,随着反应过程的逐渐进行,电堆的电流需求逐渐增大,空气流量也跟着相应地增大,最终保持在48g/s左右。仿真结果与实车报文的趋势基本一致。通过校核计算,得出MAE≤4.10%,MSE≤4.15%。
图10为空气进堆压力的仿真结果与实车报文的对比。如图所示,随着电堆电流需求的逐渐增大,空气进堆压力也跟着相应地增大,最终保持在180kPa