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 PDFInfo
- 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
Links
- 239000000446 fuel Substances 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 17
- 239000012528 membrane Substances 0.000 title claims abstract description 10
- 230000003044 adaptive effect Effects 0.000 claims abstract description 6
- 238000012938 design process Methods 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 2
- 230000004907 flux Effects 0.000 claims 1
- 230000009466 transformation Effects 0.000 claims 1
- 230000001052 transient effect Effects 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 5
- 230000005012 migration Effects 0.000 abstract 1
- 238000013508 migration Methods 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 7
- 238000012360 testing method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 230000036284 oxygen consumption Effects 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 235000003642 hunger Nutrition 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000005984 hydrogenation reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011031 large-scale manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000037351 starvation Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04305—Modeling, demonstration models of fuel cells, e.g. for training purposes
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04313—Processes 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/0438—Pressure; Ambient pressure; Flow
- H01M8/04395—Pressure; Ambient pressure; Flow of cathode reactants at the inlet or inside the fuel cell
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04313—Processes 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/0438—Pressure; Ambient pressure; Flow
- H01M8/04425—Pressure; Ambient pressure; Flow at auxiliary devices, e.g. reformers, compressors, burners
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04992—Processes 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/50—Fuel 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
Description
技术领域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
其中Psm为状态变量和输出量,u为输入量,包括进气管入口温度Tsm,in、进气管出口温度Tsm,out和空压机流量R为气体常数,Vsm为进气管体积,Ma为空气摩尔质量,为阴极入口流量。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 R is the gas constant, V sm is the volume of the intake pipe, Ma is the air molar mass, is the cathode inlet flow.
所述的步骤2)中,基于该模型设计自适应观测器对阴极腔压力和阴极入口流量进行观测,将阴极入口流量作为时变参数,阴极腔压力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 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:
其中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)首先需要定义关于与的关系式:3) First need to define about and The relation of:
4)对上式求导得:4) Derive the above formula to get:
其中求导后的共识可展开为如下形式:in The consensus after derivation can be expanded into the following form:
根据状态空间方程以及观测器公式,可转化整理为如下形式:According to the state space equation and the observer formula, it can be transformed into the following form:
由于参数θ变化相对缓慢,因此可以将其导数近似为零,则上式可以简化为如下形式:Since the parameter θ changes relatively slowly, its derivative can be approximated to zero, and the above formula can be simplified to the following form:
若能够保证A+LC<0,即L<0,则η收敛于0。If it can be guaranteed that A+LC<0, that is, L<0, then η converges to 0.
3)由于根据观测器公式,可推导出参数估计误差的动态方程:3) Due to According to the observer formula, the parameter estimation error can be deduced The dynamic equation of :
由于∑为正定对角矩阵,若观测器增益Γ>0,则收敛于0,同时Y为有界变量,能够保证也收敛于0。Since ∑ is a positive definite diagonal matrix, if the observer gain Γ>0, then Convergence to 0, and Y is a bounded variable, which can guarantee 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:
其中为空压机出口流量(kg/s),fcp(·)代表空压机MAP,ω为空压机转速(rpm),k为压力比,k=Pa/Psm,Pa为大气压(1.013×105Pa),Psm为进气歧管压力(Pa)。为入堆口流量,由于入堆口两侧压差较小,非线性模型可以简化为关于压差的线型方程(3),gca,in为入堆口流量系数,可通过参数辨识确定,Pca为阴极腔压力。Vsm为进气歧管体积(m3),Ma为空气的摩尔质量(28.84g/mol),R为气体常数(8.314J/(K·mol)),Tsm为进气歧管温度(K)。本发明在建立进气歧管模型时考虑进出口温差,将方程(1)转化为如下形式:in 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). 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:
其中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:
其中为出堆口流量,gca,out为出堆口流量系数,Prm为排气歧管压力,为背压阀流量,W为修正系数,Trm为排气管温度,为背压阀开度(rad),a1和a2为开度系数,Vrm为排气歧管体积。in is the outlet flow rate, g ca,out is the outlet flow coefficient, P rm is the exhaust manifold pressure, is the back pressure valve flow, W is the correction coefficient, T rm is the exhaust pipe temperature, 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:
其中Tca为阴极腔温度,Vca为阴极腔体积,为每秒氧气消耗量,为氧气摩尔质量(32g/mol),n为膜电极数量,Ist为电堆电流(A),F为法拉第常数(96485C/mol)。where T ca is the cathode cavity temperature, V ca is the cathode cavity volume, is the oxygen consumption per second, 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:
其中状态变量x=Psm,输入量时变参数A=0,C=1。where the state variable x=P sm , the input quantity time-varying parameters A=0, C=1.
图2为观测器内部结构示意图,具体结构形式为:Figure 2 is a schematic diagram of the internal structure of the observer, and the specific structural form is:
其中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:
为了验证自适应观测器的收敛性,首先定义关于与的关系式:In order to verify the convergence of the adaptive observer, first define the and The relation of:
对上式求导得Derive the above formula
结合公式(16)和(17),上式可以写为如下形式:Combining formulas (16) and (17), the above formula can be written in the following form:
将公式(12)、(13)和(15)带入到公式(21)中可以得到Taking formulas (12), (13) and (15) into formula (21), we can get
由于参数变化相对缓慢,因此可以将其导数近似为零,则上式可以简化为如下形式:Since the parameter changes relatively slowly, its derivative can be approximated to zero, and the above formula can be simplified to the following form:
若能够保证A+LC<0,则η收敛于0,同时根据公式(14),参数估计误差的收敛性可以通过下式保证:If A+LC<0 can be guaranteed, then η converges to 0, and at the same time According to formula (14), the convergence of parameter estimation error can be guaranteed by the following formula:
由于∑为正定对角矩阵,若观测器增益Γ>0,则收敛于0,同时Y为有界变量,能够保证也收敛于0。Since ∑ is a positive definite diagonal matrix, if the observer gain Γ>0, then Convergence to 0, and Y is a bounded variable, which can guarantee 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
以上所述仅为本发明创造的较佳实施例而已,并不用以限制本发明创造,凡在本发明创造的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明创造的保护范围之内。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)
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)
| 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)
| 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)
| 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 |
-
2021
- 2021-04-27 CN CN202110462076.6A patent/CN113299951B/en active Active
Patent Citations (3)
| 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 |