[go: up one dir, main page]

CN113299951B - Method for observing cathode pressure and flow of proton exchange membrane fuel cell - Google Patents

Method for observing cathode pressure and flow of proton exchange membrane fuel cell Download PDF

Info

Publication number
CN113299951B
CN113299951B CN202110462076.6A CN202110462076A CN113299951B CN 113299951 B CN113299951 B CN 113299951B CN 202110462076 A CN202110462076 A CN 202110462076A CN 113299951 B CN113299951 B CN 113299951B
Authority
CN
China
Prior art keywords
observer
cathode
pressure
fuel cell
flow
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110462076.6A
Other languages
Chinese (zh)
Other versions
CN113299951A (en
Inventor
江维海
李丞
孙铂奇
朱仲文
李晓伟
聂振华
王森
周炼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Caac Yangzhou Automotive Engineering Research Institute Co ltd
China Automotive Technology and Research Center Co Ltd
CATARC Tianjin Automotive Engineering Research Institute Co Ltd
Original Assignee
Caac Yangzhou Automotive Engineering Research Institute Co ltd
China Automotive Technology and Research Center Co Ltd
CATARC Tianjin Automotive Engineering Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Caac Yangzhou Automotive Engineering Research Institute Co ltd, China Automotive Technology and Research Center Co Ltd, CATARC Tianjin Automotive Engineering Research Institute Co Ltd filed Critical Caac Yangzhou Automotive Engineering Research Institute Co ltd
Priority to CN202110462076.6A priority Critical patent/CN113299951B/en
Publication of CN113299951A publication Critical patent/CN113299951A/en
Application granted granted Critical
Publication of CN113299951B publication Critical patent/CN113299951B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M8/00Fuel cells; Manufacture thereof
    • H01M8/04Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
    • H01M8/04298Processes for controlling fuel cells or fuel cell systems
    • H01M8/04305Modeling, demonstration models of fuel cells, e.g. for training purposes
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M8/00Fuel cells; Manufacture thereof
    • H01M8/04Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
    • H01M8/04298Processes for controlling fuel cells or fuel cell systems
    • H01M8/04313Processes for controlling fuel cells or fuel cell systems characterised by the detection or assessment of variables; characterised by the detection or assessment of failure or abnormal function
    • H01M8/0438Pressure; Ambient pressure; Flow
    • H01M8/04395Pressure; Ambient pressure; Flow of cathode reactants at the inlet or inside the fuel cell
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M8/00Fuel cells; Manufacture thereof
    • H01M8/04Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
    • H01M8/04298Processes for controlling fuel cells or fuel cell systems
    • H01M8/04313Processes for controlling fuel cells or fuel cell systems characterised by the detection or assessment of variables; characterised by the detection or assessment of failure or abnormal function
    • H01M8/0438Pressure; Ambient pressure; Flow
    • H01M8/04425Pressure; Ambient pressure; Flow at auxiliary devices, e.g. reformers, compressors, burners
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01MPROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M8/00Fuel cells; Manufacture thereof
    • H01M8/04Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
    • H01M8/04298Processes for controlling fuel cells or fuel cell systems
    • H01M8/04992Processes for controlling fuel cells or fuel cell systems characterised by the implementation of mathematical or computational algorithms, e.g. feedback control loops, fuzzy logic, neural networks or artificial intelligence
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/30Hydrogen technology
    • Y02E60/50Fuel cells

Landscapes

  • Engineering & Computer Science (AREA)
  • General Chemical & Material Sciences (AREA)
  • Manufacturing & Machinery (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Fuel Cell (AREA)

Abstract

The invention provides a method for observing cathode pressure and flow of a proton exchange membrane fuel cell, which comprises the following steps: s1: establishing a lumped parameter model of a fuel cell air system; s2: based on the model, designing an adaptive observer to observe the pressure of a cathode cavity and the flow of a cathode inlet; s3: the observer design process considers the temperature difference between the inlet and the outlet of the manifold, which is caused by the intercooler, so as to ensure the convergence of the observer; s4: measurable signals of the fuel cell system are input into an observer, and real-time accurate observation of cathode flow and pressure is realized through real-time iterative operation. The method for observing the cathode pressure and the flow of the proton exchange membrane fuel cell has the advantages of convenient arrangement of a sensor for signal acquisition, strong realizability, effective increase of the observation precision and control effect of transient working conditions during high fuel cell power point migration.

Description

一种质子交换膜燃料电池阴极压力及流量观测方法A method for observing the cathode pressure and flow rate of a proton exchange membrane fuel cell

技术领域technical field

本发明创造属于燃料电池技术领域,尤其是涉及一种质子交换膜燃料电池阴极压力及流量观测方法。The invention belongs to the technical field of fuel cells, and in particular relates to a method for observing the cathode pressure and flow of a proton exchange membrane fuel cell.

背景技术Background technique

随着能源与环境问题的日益严峻,新能源汽车得到了人们的广泛关注,与纯电动汽车相比,燃料电池汽车具有加氢时间短,续驶里程长等优点。但是,目前燃料电池汽车的大规模生产与应用依然受到诸多技术难点的制约,对于汽车产品来说,可靠性与耐久性是评价汽车性能的重要指标,燃料电池在汽车动态工况下需要迅速响应功率需求的快速变化,以保证汽车的动力性要求,同时要求精确控制氢气与空气的流量,保证经济性的要求,空气供应不足会导致“氧饥饿”现象的出现,影响寿命,阴极压力过低会降低内部反应效率。因此,为了能够精确控制阴极入口空气流量以及阴极腔压力,需要对其进行反馈控制,然而上述物理量无法通过传感器直接测量,需要设计观测器对其进行实时估计。With the increasingly severe energy and environmental problems, new energy vehicles have received widespread attention. Compared with pure electric vehicles, fuel cell vehicles have the advantages of shorter hydrogenation time and longer driving range. However, the current large-scale production and application of fuel cell vehicles are still restricted by many technical difficulties. For automotive products, reliability and durability are important indicators for evaluating vehicle performance, and fuel cells need to respond quickly under dynamic vehicle conditions. Rapid changes in power demand to ensure the power requirements of the car, and at the same time require precise control of the flow of hydrogen and air to ensure economical requirements. Insufficient air supply will lead to the phenomenon of "oxygen starvation", affecting life, and the cathode pressure is too low It will reduce the internal reaction efficiency. Therefore, in order to accurately control the cathode inlet air flow and the cathode cavity pressure, feedback control is required. However, the above physical quantities cannot be directly measured by sensors, and an observer needs to be designed to estimate them in real time.

发明内容SUMMARY OF THE INVENTION

有鉴于此,本发明的目的在于通过设计观测器实现阴极腔压力和入堆口空气流量的实时精确估计,根据进气歧管入口流量、进出口温度、歧管内压力等可测量信号,间接估计出入堆口流量和阴极腔压力这两个无法通过传感器直接测量或测量困难的物理量。入堆口空气流量和阴极腔压力作为燃料电池的重要状态参数以及控制反馈量,在故障诊断、状态检测和空气系统解耦控制等关键技术领域至关重要。In view of this, the purpose of the present invention is to realize the real-time accurate estimation of the cathode cavity pressure and the air flow rate of the stack inlet by designing an observer, and indirectly estimate the pressure according to the measurable signals such as the inlet flow rate of the intake manifold, the temperature of the inlet and outlet, and the pressure in the manifold. The flow in and out of the stack port and the pressure in the cathode cavity cannot be directly measured or difficult to measure by sensors. As important state parameters and control feedback quantities of fuel cells, the air flow rate at the inlet port and the cathode cavity pressure are crucial in key technical fields such as fault diagnosis, state detection and air system decoupling control.

为达到上述目的,本发明创造的技术方案是这样实现的:In order to achieve the above object, the technical scheme created by the present invention is realized like this:

实现在行车过程中的实时观测,为燃料电池的状态检测及精确控制提供技术基础,具体步骤如下:Real-time observation during driving is realized, providing a technical basis for state detection and precise control of fuel cells. The specific steps are as follows:

1)建立燃料电池空气系统集总参数模型;1) Establish a lumped parameter model of the fuel cell air system;

2)基于该模型设计自适应观测器对阴极腔压力和阴极入口流量进行观测;2) Based on this model, an adaptive observer is designed to observe the cathode cavity pressure and cathode inlet flow;

3)观测器设计过程考虑中冷器导致的歧管进出口温差,以提高观测精度,同时通过收敛性分析设计观测器参数,保证观测器收敛;3) In the observer design process, the temperature difference between the inlet and outlet of the manifold caused by the intercooler is considered to improve the observation accuracy, and at the same time, the observer parameters are designed through the convergence analysis to ensure the convergence of the observer;

4)将燃料电池系统的可测量信号(包括空压机出口温度和流量、阴极入口温度以及歧管内压力)输入到观测器中,通过实时迭代运算实现阴极流量和压力的实时精确观测。4) Input the measurable signals of the fuel cell system (including the air compressor outlet temperature and flow, the cathode inlet temperature and the pressure in the manifold) into the observer, and realize the real-time accurate observation of the cathode flow and pressure through real-time iterative operations.

所述步骤1)中,根据质量守恒定律、理想气体方程和假设条件,并考虑歧管进出口温差搭建面向控制的燃料电池空气系统集总参数模型。In the step 1), a control-oriented lumped parameter model of the fuel cell air system is built according to the law of conservation of mass, the ideal gas equation and assumptions, and considering the temperature difference between the inlet and outlet of the manifold.

所述模型假设条件具体包括:The model assumptions specifically include:

歧管内的空气视为理想气体;The air in the manifold is regarded as an ideal gas;

气体组分的压力、温度和浓度均匀分布;Uniform distribution of pressure, temperature and concentration of gas components;

阴极入口流量与压差成正比;The cathode inlet flow is proportional to the differential pressure;

不考虑热力学过程,视为绝热系统。The thermodynamic process is not considered, and it is regarded as an adiabatic system.

所述燃料电池空气系统集总参数模型的搭建需要考虑歧管进出口温差,其状态空间方程具体为:The construction of the fuel cell air system lumped parameter model needs to consider the temperature difference between the inlet and outlet of the manifold, and its state space equation is specifically:

x=Psm y=Psm x=P sm y=P sm

Figure BDA0003042659440000021
Figure BDA0003042659440000021

Figure BDA0003042659440000022
Figure BDA0003042659440000022

其中Psm为状态变量和输出量,u为输入量,包括进气管入口温度Tsm,in、进气管出口温度Tsm,out和空压机流量

Figure BDA0003042659440000023
R为气体常数,Vsm为进气管体积,Ma为空气摩尔质量,
Figure BDA0003042659440000031
为阴极入口流量。Among them, P sm is the state variable and output quantity, and u is the input quantity, including intake pipe inlet temperature T sm,in , intake pipe outlet temperature T sm,out and air compressor flow
Figure BDA0003042659440000023
R is the gas constant, V sm is the volume of the intake pipe, Ma is the air molar mass,
Figure BDA0003042659440000031
is the cathode inlet flow.

所述的步骤2)中,基于该模型设计自适应观测器对阴极腔压力和阴极入口流量进行观测,将阴极入口流量

Figure BDA0003042659440000032
作为时变参数,阴极腔压力Pca根据入口流量和进气管压力计算获得,具体观测器公式如下:In the described step 2), an adaptive observer is designed based on the model to observe the cathode cavity pressure and the cathode inlet flow rate, and the cathode inlet flow rate is
Figure BDA0003042659440000032
As a time-varying parameter, the cathode cavity pressure P ca is calculated from the inlet flow rate and the intake pipe pressure. The specific observer formula is as follows:

Figure BDA0003042659440000033
Figure BDA0003042659440000033

Figure BDA0003042659440000034
Figure BDA0003042659440000034

Figure BDA0003042659440000035
Figure BDA0003042659440000035

Figure BDA0003042659440000036
Figure BDA0003042659440000036

其中gca,in为入堆口流量系数,A、B、C和Φ为系数矩阵,Y为中间变量,L和Γ为观测器增益,Σ为任意正定对角矩阵。where g ca,in is the inlet flow coefficient, A, B, C and Φ are coefficient matrices, Y is an intermediate variable, L and Γ are observer gains, and Σ is an arbitrary positive definite diagonal matrix.

所述的步骤3)中,通过收敛性分析设计观测器参数,状态变量以及时变参数估计误差收敛性分析的具体步骤如下:In the described step 3), the observer parameters are designed by the convergence analysis, and the specific steps of the state variable and the time-varying parameter estimation error convergence analysis are as follows:

3)首先需要定义关于

Figure BDA0003042659440000037
Figure BDA0003042659440000038
的关系式:3) First need to define about
Figure BDA0003042659440000037
and
Figure BDA0003042659440000038
The relation of:

Figure BDA0003042659440000039
Figure BDA0003042659440000039

4)对上式求导得:4) Derive the above formula to get:

Figure BDA00030426594400000310
Figure BDA00030426594400000310

其中

Figure BDA00030426594400000311
求导后的共识可展开为如下形式:in
Figure BDA00030426594400000311
The consensus after derivation can be expanded into the following form:

Figure BDA00030426594400000312
Figure BDA00030426594400000312

根据状态空间方程以及观测器公式,可转化整理为如下形式:According to the state space equation and the observer formula, it can be transformed into the following form:

Figure BDA00030426594400000313
Figure BDA00030426594400000313

由于参数θ变化相对缓慢,因此可以将其导数近似为零,则上式可以简化为如下形式:Since the parameter θ changes relatively slowly, its derivative can be approximated to zero, and the above formula can be simplified to the following form:

Figure BDA00030426594400000314
Figure BDA00030426594400000314

若能够保证A+LC<0,即L<0,则η收敛于0。If it can be guaranteed that A+LC<0, that is, L<0, then η converges to 0.

3)由于

Figure BDA0003042659440000041
根据观测器公式,可推导出参数估计误差
Figure BDA0003042659440000042
的动态方程:3) Due to
Figure BDA0003042659440000041
According to the observer formula, the parameter estimation error can be deduced
Figure BDA0003042659440000042
The dynamic equation of :

Figure BDA0003042659440000043
Figure BDA0003042659440000043

由于∑为正定对角矩阵,若观测器增益Γ>0,则

Figure BDA0003042659440000044
收敛于0,同时Y为有界变量,能够保证
Figure BDA0003042659440000045
也收敛于0。Since ∑ is a positive definite diagonal matrix, if the observer gain Γ>0, then
Figure BDA0003042659440000044
Convergence to 0, and Y is a bounded variable, which can guarantee
Figure BDA0003042659440000045
also converges to 0.

相对于现有技术,本发明创造所述的一种质子交换膜燃料电池阴极压力及流量观测方法具有以下优势:Compared with the prior art, the method for observing the cathode pressure and flow rate of a proton exchange membrane fuel cell created by the present invention has the following advantages:

1)本发明通过引入自适应率,降低了对模型精度的要求,用于观测器设计的模型采用简化后的空气系统模型,有效提高运算速度,从而满足控制的实时性要求。1) The present invention reduces the requirement for model accuracy by introducing the self-adaptation rate, and the model used for the observer design adopts the simplified air system model, which effectively improves the operation speed and thus meets the real-time requirement of control.

2)本发明所设计的观测器通过收敛性分析提出了保证观测器收敛的参数范围,可通过参数设计保证观测过程在任意工况下的稳定性,满足燃料电池系统的安全性要求。2) The observer designed in the present invention proposes a parameter range to ensure the convergence of the observer through convergence analysis, and can ensure the stability of the observation process under any working conditions through parameter design, and meets the safety requirements of the fuel cell system.

3)本发明所设计的观测器输入量较少,用于信号采集的传感器布置方便,可实现性强。3) The observer designed in the present invention has less input, and the sensor used for signal acquisition is convenient to be arranged and highly implementable.

4)本发明考虑了空气系统空压机出口与入堆口的温差,在采集入堆口温度的基础上,增加了对空压机出口温度的采集,将二者同时作为观测器输入量,可消除基于单一温度的观测偏差,有效增加观测精度。4) The present invention considers the temperature difference between the air compressor outlet of the air system and the inlet of the stack, and on the basis of collecting the temperature of the inlet of the stack, the collection of the outlet temperature of the air compressor is added, and the two are simultaneously used as the input of the observer, It can eliminate the observation bias based on a single temperature and effectively increase the observation accuracy.

5)本发明通过观测器对阴极压力和阴极入堆口流量进行实时观测,为燃料电池状态检测以及故障诊断等功能的实现提供技术保障。5) The present invention conducts real-time observation of cathode pressure and cathode inlet flow rate through an observer, which provides technical support for the realization of functions such as fuel cell state detection and fault diagnosis.

6)本发明观测器所观测的阴极压力和入堆口流量可作为反馈控制量,从而有效提高燃料电池功率点迁移时瞬态工况的控制效果。6) The cathode pressure and the flow rate of the inlet port observed by the observer of the present invention can be used as feedback control quantities, thereby effectively improving the control effect of the transient working condition when the power point of the fuel cell is migrated.

附图说明Description of drawings

图1为本发明中观测器系统示意图;Fig. 1 is the schematic diagram of the observer system in the present invention;

图2为本发明中观测器内部算法结构示意图;2 is a schematic diagram of the internal algorithm structure of the observer in the present invention;

图3为本发明实例中FCU控制器硬件;Fig. 3 is FCU controller hardware in the example of the present invention;

图4为本发明中用于测试观测器的硬件在环系统架构;4 is a hardware-in-the-loop system architecture for testing an observer in the present invention;

图5为本发明中硬件在环实验平台控制界面;Fig. 5 is the hardware-in-the-loop experiment platform control interface in the present invention;

图6为本发明实例的测试工况;Fig. 6 is the test condition of the example of the present invention;

图7为本发明实例中入堆口流量观测结果;Fig. 7 is the observation result of inlet flow rate in the example of the present invention;

图8为本发明实例中阴极腔压力观测结果;Fig. 8 is the observation result of the cathode cavity pressure in the example of the present invention;

图9为本发明实例中基于观测器的过氧比控制结果。FIG. 9 shows the results of the observer-based peroxygen ratio control in the example of the present invention.

具体实施方式Detailed ways

除有定义外,以下实施例中所用的技术术语具有与本发明创造所属领域技术人员普遍理解的相同含义。以下实施例中所用的试验试剂,如无特殊说明,均为常规生化试剂;所述实验方法,如无特殊说明,均为常规方法。Unless otherwise defined, technical terms used in the following embodiments have the same meanings as commonly understood by those skilled in the art to which the present invention belongs. The test reagents used in the following examples are conventional biochemical reagents unless otherwise specified; the experimental methods are conventional methods unless otherwise specified.

下面结合实施例及附图来详细说明本发明创造。The present invention will be described in detail below with reference to the embodiments and the accompanying drawings.

图1为本发明实例所构建的质子交换膜燃料电池阴极压力及流量观测系统,该实例中的燃料电池空气系统主要包括空压机、中冷器、加湿器、冷凝器、背压阀以及各类传感器,该观测方法要求燃料电池空气系统中安装有温度传感器、压力传感器以及流量传感器,且流量传感器安装于空压机入口处用于测量进气歧管的入口空气流量,压力传感器安装于入堆口处用于测量进气歧管内的压力,两个温度传感器分别安装于空压机出口处和入堆口处用于测量空气经过中冷器后产生的温差,如图1所示,观测器通过采集各传感器的测量信号,通过内部运算实时计算得到入堆口空气流量以及阴极腔压力的估计值。Fig. 1 is a proton exchange membrane fuel cell cathode pressure and flow observation system constructed by an example of the present invention. The fuel cell air system in this example mainly includes an air compressor, an intercooler, a humidifier, a condenser, a back pressure valve and various This observation method requires that a temperature sensor, a pressure sensor and a flow sensor are installed in the fuel cell air system, and the flow sensor is installed at the inlet of the air compressor to measure the inlet air flow of the intake manifold, and the pressure sensor is installed at the inlet. The stack port is used to measure the pressure in the intake manifold, and two temperature sensors are installed at the air compressor outlet and the stack inlet to measure the temperature difference after the air passes through the intercooler, as shown in Figure 1. The device collects the measurement signals of each sensor, and calculates in real time through the internal operation to obtain the estimated value of the air flow rate of the inlet port and the pressure of the cathode cavity.

在本实例中,设计观测器之前需要先构建燃料电池空气系统集总参数模型,该模型为面向控制的模型,为了满足观测器的实时性要求,需要对模型进行简化处理,在简化过程中包括以下假设:In this example, the fuel cell air system lumped parameter model needs to be built before designing the observer. This model is a control-oriented model. In order to meet the real-time requirements of the observer, the model needs to be simplified. The simplified process includes: The following assumptions:

1)歧管内的空气视为理想气体;1) The air in the manifold is regarded as an ideal gas;

2)气体组分的压力、温度和浓度均匀分布;2) The pressure, temperature and concentration of gas components are uniformly distributed;

3)阴极入口流量与压差成正比;3) The cathode inlet flow is proportional to the pressure difference;

4)不考虑热力学过程,视为绝热系统。4) The thermodynamic process is not considered, and it is regarded as an adiabatic system.

基于上述假设,建立燃料电池空气系统模型,包括进气歧管模型、排气歧管模型和阴极腔模型,其中进气歧管模型又包含空压机子模型、入堆口子模型以及歧管内压力动态模型,其状态空间方程可描述为:Based on the above assumptions, the fuel cell air system model is established, including the intake manifold model, the exhaust manifold model and the cathode cavity model. The intake manifold model also includes the air compressor sub-model, the stack inlet sub-model and the pressure in the manifold. The dynamic model, its state space equation can be described as:

Figure BDA0003042659440000061
Figure BDA0003042659440000061

Figure BDA0003042659440000062
Figure BDA0003042659440000062

Figure BDA0003042659440000063
Figure BDA0003042659440000063

其中

Figure BDA0003042659440000064
为空压机出口流量(kg/s),fcp(·)代表空压机MAP,ω为空压机转速(rpm),k为压力比,k=Pa/Psm,Pa为大气压(1.013×105Pa),Psm为进气歧管压力(Pa)。
Figure BDA0003042659440000065
为入堆口流量,由于入堆口两侧压差较小,非线性模型可以简化为关于压差的线型方程(3),gca,in为入堆口流量系数,可通过参数辨识确定,Pca为阴极腔压力。Vsm为进气歧管体积(m3),Ma为空气的摩尔质量(28.84g/mol),R为气体常数(8.314J/(K·mol)),Tsm为进气歧管温度(K)。本发明在建立进气歧管模型时考虑进出口温差,将方程(1)转化为如下形式:in
Figure BDA0003042659440000064
is the air compressor outlet flow (kg/s), f cp (·) represents the air compressor MAP, ω is the air compressor rotational speed (rpm), k is the pressure ratio, k=P a /P sm , and P a is the atmospheric pressure (1.013×10 5 Pa), and P sm is the intake manifold pressure (Pa).
Figure BDA0003042659440000065
is the flow rate of the reactor inlet, because the pressure difference between the two sides of the reactor inlet is small, the nonlinear model can be simplified to the linear equation (3) about the pressure difference, g ca,in is the inlet flow coefficient, which can be determined by parameter identification , P ca is the cathode cavity pressure. V sm is the intake manifold volume (m 3 ), Ma is the molar mass of air (28.84 g /mol), R is the gas constant (8.314 J/(K·mol)), and T sm is the intake manifold temperature (K). The present invention considers the inlet and outlet temperature difference when establishing the intake manifold model, and transforms equation (1) into the following form:

Figure BDA0003042659440000066
Figure BDA0003042659440000066

其中Tsm,in为进气歧管入口温度、Tsm,out为进气歧管出口温度。where T sm,in is the intake manifold inlet temperature, and T sm,out is the intake manifold outlet temperature.

排气歧管模型可分为出堆口子模型、背压阀模型以及歧管内压力动态子模型。其中出堆口模型和入堆口模型类似,采用关于阴极腔和排气歧管压差的线性方程描述;背压阀的结构形式为蝶阀,其模型可采用有变效截面积的喷嘴模型来描述。其状态空间方程可描述为:The exhaust manifold model can be divided into a stack outlet sub-model, a back pressure valve model and a pressure dynamic sub-model in the manifold. The reactor outlet model is similar to the reactor inlet model, and is described by a linear equation about the pressure difference between the cathode cavity and the exhaust manifold; the structure of the back pressure valve is a butterfly valve, and its model can be a nozzle model with a variable effective cross-sectional area. describe. Its state space equation can be described as:

Figure BDA0003042659440000071
Figure BDA0003042659440000071

Figure BDA0003042659440000072
Figure BDA0003042659440000072

Figure BDA0003042659440000073
Figure BDA0003042659440000073

Figure BDA0003042659440000074
Figure BDA0003042659440000074

Figure BDA0003042659440000075
Figure BDA0003042659440000075

其中

Figure BDA0003042659440000076
为出堆口流量,gca,out为出堆口流量系数,Prm为排气歧管压力,
Figure BDA0003042659440000077
为背压阀流量,W为修正系数,Trm为排气管温度,
Figure BDA0003042659440000078
为背压阀开度(rad),a1和a2为开度系数,Vrm为排气歧管体积。in
Figure BDA0003042659440000076
is the outlet flow rate, g ca,out is the outlet flow coefficient, P rm is the exhaust manifold pressure,
Figure BDA0003042659440000077
is the back pressure valve flow, W is the correction coefficient, T rm is the exhaust pipe temperature,
Figure BDA0003042659440000078
is the back pressure valve opening (rad), a 1 and a 2 are the opening coefficients, and V rm is the exhaust manifold volume.

阴极腔模型可根据入堆口子模型、出堆口子模型以及氧气消耗子模型来描述阴极腔压力的动态变化,其状态空间方程可描述为:The cathode cavity model can describe the dynamic change of the cathode cavity pressure according to the reactor inlet submodel, the reactor outlet submodel and the oxygen consumption submodel. Its state space equation can be described as:

Figure BDA0003042659440000079
Figure BDA0003042659440000079

Figure BDA00030426594400000710
Figure BDA00030426594400000710

其中Tca为阴极腔温度,Vca为阴极腔体积,

Figure BDA00030426594400000711
为每秒氧气消耗量,
Figure BDA00030426594400000712
为氧气摩尔质量(32g/mol),n为膜电极数量,Ist为电堆电流(A),F为法拉第常数(96485C/mol)。where T ca is the cathode cavity temperature, V ca is the cathode cavity volume,
Figure BDA00030426594400000711
is the oxygen consumption per second,
Figure BDA00030426594400000712
is the oxygen molar mass (32 g/mol), n is the number of membrane electrodes, I st is the stack current (A), and F is the Faraday constant (96485 C/mol).

根据本实例中建立的燃料电池空气系统模型设计状态观测器,为了便于描述,首先将进气歧管模型写成状态空间方程的形式:The state observer is designed according to the fuel cell air system model established in this example. For the convenience of description, the intake manifold model is first written in the form of a state space equation:

Figure BDA0003042659440000081
Figure BDA0003042659440000081

其中状态变量x=Psm,输入量

Figure BDA0003042659440000082
时变参数
Figure BDA0003042659440000083
A=0,
Figure BDA0003042659440000084
C=1。where the state variable x=P sm , the input quantity
Figure BDA0003042659440000082
time-varying parameters
Figure BDA0003042659440000083
A=0,
Figure BDA0003042659440000084
C=1.

图2为观测器内部结构示意图,具体结构形式为:Figure 2 is a schematic diagram of the internal structure of the observer, and the specific structural form is:

Figure BDA0003042659440000085
Figure BDA0003042659440000085

Figure BDA0003042659440000086
Figure BDA0003042659440000086

Figure BDA0003042659440000087
Figure BDA0003042659440000087

Figure BDA0003042659440000088
Figure BDA0003042659440000088

其中L和Γ为观测器增益,∑为正定对角矩阵,定义状态变量估计误差与参数估计误差分别为如下形式:where L and Γ are the observer gains, Σ is a positive definite diagonal matrix, and the state variable estimation error and parameter estimation error are defined as follows:

Figure BDA0003042659440000089
Figure BDA0003042659440000089

Figure BDA00030426594400000810
Figure BDA00030426594400000810

为了验证自适应观测器的收敛性,首先定义关于

Figure BDA00030426594400000811
Figure BDA00030426594400000812
的关系式:In order to verify the convergence of the adaptive observer, first define the
Figure BDA00030426594400000811
and
Figure BDA00030426594400000812
The relation of:

Figure BDA00030426594400000813
Figure BDA00030426594400000813

对上式求导得Derive the above formula

Figure BDA00030426594400000814
Figure BDA00030426594400000814

结合公式(16)和(17),上式可以写为如下形式:Combining formulas (16) and (17), the above formula can be written in the following form:

Figure BDA00030426594400000815
Figure BDA00030426594400000815

将公式(12)、(13)和(15)带入到公式(21)中可以得到Taking formulas (12), (13) and (15) into formula (21), we can get

Figure BDA00030426594400000816
Figure BDA00030426594400000816

由于参数变化相对缓慢,因此可以将其导数近似为零,则上式可以简化为如下形式:Since the parameter changes relatively slowly, its derivative can be approximated to zero, and the above formula can be simplified to the following form:

Figure BDA0003042659440000091
Figure BDA0003042659440000091

若能够保证A+LC<0,则η收敛于0,同时

Figure BDA0003042659440000092
根据公式(14),参数估计误差的收敛性可以通过下式保证:If A+LC<0 can be guaranteed, then η converges to 0, and at the same time
Figure BDA0003042659440000092
According to formula (14), the convergence of parameter estimation error can be guaranteed by the following formula:

Figure BDA0003042659440000093
Figure BDA0003042659440000093

由于∑为正定对角矩阵,若观测器增益Γ>0,则

Figure BDA0003042659440000094
收敛于0,同时Y为有界变量,能够保证
Figure BDA0003042659440000095
也收敛于0。Since ∑ is a positive definite diagonal matrix, if the observer gain Γ>0, then
Figure BDA0003042659440000094
Convergence to 0, and Y is a bounded variable, which can guarantee
Figure BDA0003042659440000095
also converges to 0.

为了验证本发明观测器的效果,采用MATLAB/Simulink搭建燃料电池模型,并通过硬件在环(HIL)实验进行验证。测试环境基于NI的Veristand进行搭建,控制器硬件如图3所示,采用TC277D作为运算核心,并集成了AD、CAN等数据通讯模块,具备PWM输入输出、高低边驱动、高压采集以及LIN通讯等功能,硬件在环测试系统架构如图4所示,硬件在环实验平台控制界面如图5所示。燃料电池空气系统模型及观测器参数由表1给出,实验工况所给出的目标电流和阴极腔压力如图6所示,目标电流在100-150A之间阶跃变化,目标压力在2-2.2bar之间阶跃变化,相应的入堆口流量的观测结果如图7所示,阴极腔压力的观测结果如图8所示,流量与压力的传感器测量值采用硬件在环模型输出值近似表示,通过与观测器估计值的对比可以看到,无论在稳态工况还是在瞬态工况下都具有较高的观测精度,且相比于基于单一温度输入的观测器来说具有更高的瞬态精度。在目标电流从140A到130A阶跃变化以及目标压力从2.1bar到2.05bar阶跃变化的瞬态工况下,从局部放大图可以看到观测值具有较好的跟随效果。为了进一步验证观测器对阴极空气系统控制的有益效果,图9对比了有无观测器情况下过氧比的控制效果,在瞬态工况下,基于观测器的控制方法具有更好的瞬态控制精度。In order to verify the effect of the observer of the present invention, a fuel cell model is built by using MATLAB/Simulink, and verified by hardware-in-the-loop (HIL) experiments. The test environment is built based on NI's Veristand. The controller hardware is shown in Figure 3, using TC277D as the computing core, and integrating AD, CAN and other data communication modules, with PWM input and output, high and low side drive, high voltage acquisition and LIN communication, etc. Function, the hardware-in-the-loop test system architecture is shown in Figure 4, and the hardware-in-the-loop experimental platform control interface is shown in Figure 5. The fuel cell air system model and observer parameters are given in Table 1. The target current and cathode cavity pressure given by the experimental conditions are shown in Figure 6. The target current changes stepwise between 100-150A, and the target pressure is 2 The step change between -2.2bar, the corresponding observation results of the inlet flow are shown in Figure 7, and the observation results of the cathode cavity pressure are shown in Figure 8. The measured values of the flow and pressure sensors adopt the output value of the hardware-in-the-loop model Approximate representation, through the comparison with the estimated value of the observer, it can be seen that it has higher observation accuracy in both steady state and transient conditions, and has a higher observation accuracy than the observer based on a single temperature input. Higher transient accuracy. Under the transient conditions of the step change of target current from 140A to 130A and the step change of target pressure from 2.1bar to 2.05bar, it can be seen from the partial enlarged view that the observed value has a good follow-up effect. In order to further verify the beneficial effect of the observer on the control of the cathode air system, Figure 9 compares the control effect of the peroxygen ratio with and without the observer. Under the transient condition, the observer-based control method has better transient control precision.

表1模型及观测器参数Table 1 Model and observer parameters

Figure BDA0003042659440000096
Figure BDA0003042659440000096

Figure BDA0003042659440000101
Figure BDA0003042659440000101

以上所述仅为本发明创造的较佳实施例而已,并不用以限制本发明创造,凡在本发明创造的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明创造的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the within the scope of protection of the present invention.

Claims (3)

1. A method for observing the cathode pressure and flow of a proton exchange membrane fuel cell is characterized in that: the method comprises the following steps:
s1: establishing a lumped parameter model of the fuel cell air system;
s2: based on the model, designing an adaptive observer to observe the pressure of a cathode cavity and the flow of a cathode inlet;
s3: the observer design process considers the temperature difference between an inlet and an outlet of the manifold, which is caused by the intercooler, so that the convergence of the observer is ensured;
s4: inputting a measurable signal of the fuel cell system into an observer, and realizing real-time accurate observation of cathode flow and pressure through real-time iterative operation;
the construction of the fuel cell air system lumped parameter model needs to consider the temperature difference between an inlet and an outlet of a manifold, and the state space equation is as follows:
x=P sm y=P sm
Figure FDA0003811752240000011
Figure FDA0003811752240000012
wherein P is sm Is the state variable and the output quantity, u is the input quantity,
T sm,in is the inlet temperature of the air inlet pipe,
T sm,out is the temperature at the outlet of the air inlet pipe,
Figure FDA0003811752240000013
is the flow rate of the air compressor,
r is the gas constant, and R is the gas constant,
V sm is the volume of the air inlet pipe,
M a the molar mass of the air is the mass of the air,
Figure FDA0003811752240000014
is the cathode inlet flow rate;
in S2, the adaptive observer is designed based on the model to observe the cathode cavity pressure and the cathode inlet flow rate, and the cathode inlet flow rate is measured
Figure FDA0003811752240000021
Cathode cavity pressure P as a time varying parameter ca The observer is obtained by calculation according to the inlet flow and the pressure of the air inlet pipe, and the specific observer formula is as follows:
Figure FDA0003811752240000022
Figure FDA0003811752240000023
Figure FDA0003811752240000024
Figure FDA0003811752240000025
wherein g is ca,in For the flux coefficient of the reactor inlet, A, B, C and phi are coefficient matrixes, Y is an intermediate variable, L and gamma are observer gains, and sigma is an arbitrary positive definite diagonal matrix;
in S3, the convergence analysis of the observer parameters, state variables, and time-varying parameter estimation errors is designed through convergence analysis, which includes the following specific steps:
1) first, it is necessary to define
Figure FDA0003811752240000026
And
Figure FDA0003811752240000027
the relation of (c):
Figure FDA0003811752240000028
2) the above formula is derived:
Figure FDA0003811752240000029
wherein
Figure FDA00038117522400000210
The derived consensus can be expanded into the following form:
Figure FDA00038117522400000211
according to the state space equation and the observer formula, the transformation and the adjustment can be in the following form:
Figure FDA00038117522400000212
since the parameter θ changes relatively slowly, its derivative can be approximated to zero, and the above equation can be simplified to the following form:
Figure FDA00038117522400000213
if the A + LC is less than 0, namely L is less than 0, then eta converges to 0;
due to the fact that
Figure FDA00038117522400000214
Deducing parameter estimation error according to observer formula
Figure FDA00038117522400000215
The dynamic equation of (c):
Figure FDA00038117522400000216
since sigma is a positive definite diagonal matrix, if observer gain gamma > 0, then
Figure FDA0003811752240000031
Convergence to 0, with Y being a bounded variable, can be guaranteed
Figure FDA0003811752240000032
Also converging to 0.
2. The method for observing the cathode pressure and flow of a proton exchange membrane fuel cell according to claim 1, wherein: in S1, the total parameter model builds a control-oriented fuel cell air system lumped parameter model according to the mass conservation law, the ideal gas equation and the assumed conditions, and according to the manifold inlet-outlet temperature difference, where the model assumed conditions specifically include: the air in the manifold is considered an ideal gas; the pressure, temperature and concentration of the gas component are uniformly distributed; the cathode inlet flow is proportional to the pressure differential; irrespective of the thermodynamic process, an adiabatic system is considered.
3. The method for observing the cathode pressure and flow of a proton exchange membrane fuel cell according to claim 1, wherein: the measurable signals in S4 include air compressor outlet temperature and flow, cathode inlet temperature, and manifold pressure.
CN202110462076.6A 2021-04-27 2021-04-27 Method for observing cathode pressure and flow of proton exchange membrane fuel cell Active CN113299951B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110462076.6A CN113299951B (en) 2021-04-27 2021-04-27 Method for observing cathode pressure and flow of proton exchange membrane fuel cell

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110462076.6A CN113299951B (en) 2021-04-27 2021-04-27 Method for observing cathode pressure and flow of proton exchange membrane fuel cell

Publications (2)

Publication Number Publication Date
CN113299951A CN113299951A (en) 2021-08-24
CN113299951B true CN113299951B (en) 2022-09-27

Family

ID=77320422

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110462076.6A Active CN113299951B (en) 2021-04-27 2021-04-27 Method for observing cathode pressure and flow of proton exchange membrane fuel cell

Country Status (1)

Country Link
CN (1) CN113299951B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102022103159A1 (en) * 2022-02-10 2023-08-10 Bayerische Motoren Werke Aktiengesellschaft Method and device for determining the oxidizing agent mass flow in an energy converter
CN114595531B (en) * 2022-03-11 2022-10-14 武汉雄韬氢雄燃料电池科技有限公司 Air inlet manifold design method based on double-stack flow distribution consistency
CN119481170A (en) * 2024-12-24 2025-02-18 国创氢能科技有限公司 A control method for fuel cell oxygen supply system based on BP neural network and active disturbance rejection control
CN119890356A (en) * 2025-03-27 2025-04-25 合肥工业大学 Shutdown purging method for fuel cell system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007165082A (en) * 2005-12-13 2007-06-28 Nissan Motor Co Ltd Fuel cell system
CN108241128A (en) * 2018-01-09 2018-07-03 西南交通大学 A State Estimation Method for Proton Exchange Membrane Fuel Cell System
CN212011145U (en) * 2020-03-12 2020-11-24 吉林大学 A fuel cell with decoupled control

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8323841B2 (en) * 2008-05-06 2012-12-04 GM Global Technology Operations LLC Anode loop observer for fuel cell systems
CN111707945A (en) * 2020-05-11 2020-09-25 同济大学 A method for observing the internal state of a proton exchange membrane fuel cell

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007165082A (en) * 2005-12-13 2007-06-28 Nissan Motor Co Ltd Fuel cell system
CN108241128A (en) * 2018-01-09 2018-07-03 西南交通大学 A State Estimation Method for Proton Exchange Membrane Fuel Cell System
CN212011145U (en) * 2020-03-12 2020-11-24 吉林大学 A fuel cell with decoupled control

Also Published As

Publication number Publication date
CN113299951A (en) 2021-08-24

Similar Documents

Publication Publication Date Title
CN113299951B (en) Method for observing cathode pressure and flow of proton exchange membrane fuel cell
CN109524693B (en) Model Predictive Control Method for Fuel Cell Air Supply System
CN112072142B (en) Fuel cell control method and system based on model predictive control
CN111342086B (en) A fuel cell air-peroxygen ratio and flow pressure coordinated control method and system
CN111261903B (en) Model-based proton exchange membrane fuel cell impedance online estimation method
CN103016184B (en) Electronic controlled diesel egr system NOx discharge monitoring method and device
CN108241128B (en) A State Estimation Method for Proton Exchange Membrane Fuel Cell System
CN106405431A (en) Fuel cell automation test system
CN113675444B (en) Fuel cell air system decoupling control method and device and storage medium
CN112635802A (en) Hydrogen control method for vehicle proton exchange membrane fuel cell system
CN102663219A (en) Fuel cell output prediction method and system based on mixing model
CN108091909A (en) It is a kind of based on optimal peroxide than fuel battery air flow control methods
CN111063919B (en) Fuel cell anode state monitoring method
CN111707945A (en) A method for observing the internal state of a proton exchange membrane fuel cell
CN110611111A (en) Method for measuring instantaneous hydrogen flow of fuel cell hydrogen system
CN111196184A (en) Vehicle, and control method and control device for vehicle fuel cell air intake system
CN111952646A (en) A decoupling control method and system for a fuel cell air system
CN109684704B (en) An online calibration method of engine intake air flow based on velocity density model
CN115188999A (en) Control method of air supply loop of fuel cell system
CN116093385A (en) A Method for Estimating Nitrogen Concentration in Fuel Cell Anode Based on Multipoint Voltage
CN115692783B (en) Control method and system for vehicle fuel cell air system
CN113644301A (en) Cathode flow and pressure control method and system for vehicle fuel cell
CN120015877A (en) Bypass air flow estimation method, device and stack air flow estimation method
US9564648B2 (en) Anode leak location detection
CN117832546B (en) A control-oriented mechanism-data modeling approach for fuel cell air systems

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant