Peroxide ratio actual control method and system based on C15-20 air compressor numerical modeling

Abstract

The invention discloses a peroxide ratio actual control method and system based on C15-20 air compressor numerical modeling. The method comprises the steps that 1, numerical extraction is conducted on the air mass flow and the compression ratio at each rotating speed in a C15-20 air compressor MAP; 2, fitting a tool box by utilizing a curve of MATLAB (matrix laboratory) to obtain a binary quartic polynomial fitting result after regularization of the air mass flow with respect to the compression ratio and the rotating speed center of the air compressor; 3, based on the fitting result, the actual compression ratio and rotating speed of the air compressor and the voltage value of the driving motor when the peroxide ratio is equal to 2 are calculated with the load current step as external disturbance; and step 4, performing PID control on the voltage of the driving motor by taking the peroxide ratio = 2 as a control target to obtain a control result of the voltage of the driving motor under the load current step. According to the method, the numerical modeling work of the ROTREX C15 series air compressor is perfected, the control result of the C15-20 air compressor is quantified, and the application range of modeling simulation and control of the ROTREX C15 series air compressor is enlarged.

Description

Technical Field

[0001] The present invention belongs to the technical field of fuel cells, and in particular relates to a method and system for actually controlling an excess oxygen ratio based on numerical modeling of a C15-20 air compressor. Background Art

[0002] As a leading field of hydrogen energy, hydrogen fuel cells have received extensive attention. PEMFC (Proton Exchange Membrane Fuel Cell, PEMFC) is currently the most widely used and most promising automotive hydrogen fuel cell due to its advantages such as short startup time and low operating temperature.

[0003] The PEMFC engine system is the core component of the fuel cell and determines the overall working condition of the battery. Establishing an accurate PEMFC engine system simulation model is a prerequisite for the commercialization of PEMFC. The calculation accuracy of the air mass flow delivered by the air compressor determines the calculation accuracy of the oxygen pressure, nitrogen pressure and output voltage after the air enters the battery stack, which is the most important prerequisite for the accuracy of PEMFC engine system modeling and control.

[0004] In recent years, with the popularization of new energy vehicles, the power demand for fuel cells has become increasingly higher. Figure 1 As shown in the figure, the C15 series air compressor produced by Denmark's ROTREX company can be used in PEMFC with an output power range of 60~175 kW, and is a commonly used air compressor in high-power fuel cell companies. The C15 series air compressor is mainly divided into three models: C15-16, C15-60 and C15-20. Among them, the numerical modeling of the C15-16 and C15-60 model air compressors has been completed by relevant personnel.

[0005] In the existing technology, researchers only determined the air mass flow rate of the C15-20 air compressor under a certain working condition, and the MAP (Manifold Absolute Pressure Sensor, MAP, intake pressure sensor) diagram of the entire air compressor is not universal, which limits its modeling accuracy at other temperatures or pressures. The modeling method of simulating and calculating the geometric model of the air compressor through fluid mechanics has many parameters, complex models, large calculations, low efficiency, and cannot be coupled with the intake control strategy, making it difficult to apply to actual fuel cell systems. Summary of the invention

[0006] The purpose of the present invention is to provide a method and system for actually controlling the excess oxygen ratio based on numerical modeling of a C15-20 air compressor, which can truly reflect the working process of the C15-20 air compressor and realize accurate modeling and actual control application of all models of C15 series air compressors, thereby providing data support for the application of this series of air compressors in fuel cell systems.

[0007] The present invention provides a method for actually controlling the excess oxygen ratio based on numerical modeling of a C15-20 air compressor, comprising the following steps:

[0008] Step 1, based on the image digitization tool of Origin, numerically extract the air mass flow and compression ratio at each speed in the MAP diagram of the C15-20 air compressor;

[0009] Step 2, based on the extracted values, using the curve fitting toolbox of MATLAB, obtain the bivariate quartic polynomial fitting result of the air mass flow rate after regularization with respect to the compression ratio and the center of the speed of the air compressor;

[0010] Step 3, based on the obtained fitting results, taking the load current step as the external disturbance, calculate the actual air compressor compression ratio, speed and driving motor voltage value when the oxygen ratio is 2;

[0011] Step 4, taking the overoxygen ratio = 2 as the control target, PID control is performed on the drive motor voltage to obtain the control result of the drive motor voltage under the load current step.

[0012] Furthermore, the equation of the two-variable quartic polynomial in step 2 is as follows:

[0013] (1)

[0014] In the formula, is the air mass flow rate of the air compressor; for The compression ratio of the air compressor after center regularization; for The air compressor speed after center regularization, unit: rpm; i and j are and The number of items; are the polynomial coefficients.

[0015] Furthermore, the expression of the center regularization in step 2 is as follows:

[0016] (2)

[0017] is the compression ratio; is the air compressor speed.

[0018] Furthermore, the step 2 further comprises:

[0019] The residual sum of squares, R square and root mean square were used as error indicators to perform polynomial fitting error analysis and obtain the analysis results.

[0020] Furthermore, the step 3 specifically includes:

[0021] The dynamic equation of the air compressor when working is as follows:

[0022] (3)

[0023] In the formula, is the air compressor input torque, unit: N·m; is the actual torque of the air compressor, unit: N·m; is the specific heat capacity of air at constant pressure, which is 1004 J·kg -1 ·K -1 ; and is the motor constant driving the air compressor, in units of V·rad -1 ·s and N·m·A -1 ; and are the motor mechanical efficiency and the air compressor efficiency, respectively, which are dimensionless variables; is the driving motor voltage, unit: V; is the air mass flow rate of the air compressor, unit: kg·s -1 ; is the resistance of the motor armature, unit: Ω; The temperature of air flowing out of the air compressor, unit: K; is the temperature of air entering the air compressor, which is 298.15 K; is the motor angular velocity, unit: rad·s -1 ; is the angular velocity of the air compressor impeller, unit: rad·s -1 ; is the air compressor rotation inertia, unit: kg·m 2 ; is the air mass flow rate under actual working conditions, unit: kg·s -1 , calculated according to the following formula:

[0024] (4)

[0025] In the formula, For standard atmospheric pressure, take 101325 Pa;

[0026] Setting the actual excess oxygen ratio Limit the air mass flow rate involved in the reaction:

[0027] (5)

[0028] In the formula, is the number of single cells, take 762; is the load current, unit: A; is the relative molecular mass of oxygen, taken as 0.032 kg·mol -1 ; F is Faraday constant, taken as 96485 C·mol -1 .

[0029] Assuming the air compressor inlet pressure is 1 bar, the compression ratio change rate is calculated as follows:

[0030] (6)

[0031] In the formula, is the thermodynamic constant of air, which is 286.9 J·kg -1 ·K -1 ; is the volume of the gas supply pipeline, unit: m 3 ; is the air mass flow rate of the air supply pipeline, unit: kg·s -1 , calculated according to the following formula:

[0032] (7)

[0033] In the formula, is the cathode flow field pressure, unit: bar; is the air supply resistance coefficient, which is 0.3629 kg·s -1 bar -1 ; is the outlet pipe pressure of the air compressor, unit: bar;

[0034] set up , calculated ; Then set follow The value of the change, according to formula (6), The rate of change is 0, at this time ,because , in order to determine ; Determined according to MAP diagram or Any one of the ; When in formula (3) When the rate of change of is 0, we can obtain , and finally obtain .

[0035] Furthermore, the step 4 specifically includes:

[0036] For the target peroxygen ratio OER ref Compared with the actual peroxygen ratio The difference e(t) is controlled using the PID algorithm , as shown in formulas (11) to (12):

[0037] (11)

[0038] (12)

[0039] Pick =50, =5 and =0, so follow change.

[0040] The present invention also provides an actual control system for the excess oxygen ratio of a C15-20 air compressor based on numerical modeling, comprising a control module, wherein the control module executes the actual control method for the excess oxygen ratio based on numerical modeling of a C15-20 air compressor.

[0041] The present invention also provides a non-transitory computer-readable storage medium, which stores computer instructions. When the computer instructions are executed by a processor, the actual control method of the excess oxygen ratio based on numerical modeling of a C15-20 air compressor is implemented.

[0042] The present invention also provides an electronic device, comprising:

[0043] A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the actual control method of the excess oxygen ratio based on numerical modeling of a C15-20 air compressor by executing the computer instructions.

[0044] Through the above scheme, through the actual control method and system of the excess oxygen ratio based on the numerical modeling of the C15-20 air compressor, the numerical relationship between the speed, compression ratio and air mass flow rate of the C15-20 air compressor can be accurately reflected, the numerical modeling work of the ROTREX C15 series air compressors is improved, the control results of the C15-20 air compressors are quantified, and the application scope of the modeling, simulation and control of the ROTREX C15 series air compressors is expanded.

[0045] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention and implement it according to the contents of the specification, the following is a detailed description of the preferred embodiments of the present invention in conjunction with the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] Figure 1 This is the structural diagram of the C15 series air compressor produced by ROTREX of Denmark;

[0047] Figure 2 It is a flow chart of the actual control method of the excess oxygen ratio based on numerical modeling of a C15-20 air compressor of the present invention;

[0048] Figure 3 It is a schematic diagram of numerical modeling of an air compressor and a practical example based on the modeling in an embodiment of the present invention;

[0049] Figure 4 is a MAP diagram of a C15-20 air compressor in one embodiment of the present invention;

[0050] Figure 5 is a fitting result diagram in one embodiment of the present invention;

[0051] Figure 6 is a graph showing the speed control result of an air compressor in one embodiment of the present invention;

[0052] Figure 7 is a diagram of the voltage control result of the driving motor in one embodiment of the present invention;

[0053] Figure 8 It is a structural schematic diagram of an electronic device of the present invention. DETAILED DESCRIPTION

[0054] The specific implementation of the present invention is further described in detail below in conjunction with the accompanying drawings and examples. The following examples are used to illustrate the present invention, but are not intended to limit the scope of the present invention.

[0055] Ginseng Figure 2 As shown, this embodiment provides an actual control method of excess oxygen ratio (ORE) based on numerical modeling of a C15-20 air compressor, comprising the following steps:

[0056] Step S1, based on the image digitization tool of Origin, numerically extract the air mass flow and compression ratio at each speed in the MAP map of the C15-20 air compressor. According to the intake pressure sensor (MAP) map of the C15-20 air compressor, the air mass flow and compression ratio at each speed in the MAP map are numerically extracted using the image digitization function of Origin.

[0057] Step S2, based on the extracted values, using the curve fitting toolbox of MATLAB, obtain the binary quartic polynomial fitting result of the air mass flow rate with respect to the air compressor compression ratio and the center regularization of the speed. Using the curve fitting toolbox of MATLAB, a binary quartic polynomial equation of the air mass flow rate with respect to the air compressor speed and compression ratio is established, and the obtained equation fitting accuracy is as high as 99.3%.

[0058] Step S3, based on the obtained fitting result, taking the load current step as the external disturbance, calculate the actual air compressor compression ratio, speed and driving motor voltage value when the oxygen ratio is 2.

[0059] Step S4, taking the overoxygen ratio = 2 as the control target, performing PID control on the drive motor voltage to obtain the control result of the drive motor voltage under the load current step.

[0060] The present invention is described in further detail below.

[0061] Ginseng Figure 3 As shown, the numerical modeling is based on the MAP (intake pressure sensor) diagram of the centrifugal air compressor in the engine system, and a binary quartic polynomial of the air mass flow rate on the compression ratio and speed of the air compressor is established and an error analysis is performed; the actual control example application is to use the load current as the external disturbance, based on the calculation principle of the compression ratio and speed of the air compressor, calculate the actual compression ratio and speed when OER=2, and then obtain the actual drive motor voltage, and obtain the control result of the drive motor voltage through the PID algorithm. The specific calculation method is as follows:

[0062] 1. Numerical modeling of air compressor based on MAP diagram

[0063] like Figure 4 As shown, Figure 4 The relationship between compression ratio and air mass flow is described in the paper. The compression ratio (ratio of outlet pressure to inlet pressure) and speed of the air compressor determine its mass flow. Therefore, the specific steps of numerical modeling are as follows: Based on the MAP diagram of the air compressor, the following compression ratio is constructed by using a binary quartic polynomial fitting: , air compressor speed Compressor air mass flow The relationship is:

[0064] (1)

[0065] In the formula, for The compression ratio of the air compressor after center regularization; for The air compressor speed after center regularization, unit: rpm; i and j are and The number of items; is the polynomial coefficient. The specific expression of center regularization is as follows:

[0066] (2)

[0067] The polynomial coefficients obtained based on formula (4) are shown in Table 1. The specific fitting results are as follows: Figure 5 shown.

[0068] Table 1 Fitting polynomial coefficients Polynomial coefficient symbols Polynomial coefficient values <![CDATA[p 00 ]]> 0.0765 <![CDATA[p 10 ]]> 0.2082 <![CDATA[p 01 ]]> -0.167 <![CDATA[p 20 ]]> -0.2627 <![CDATA[p 11 ]]> 0.5192 <![CDATA[p 02 ]]> -0.2122 <![CDATA[p 30 ]]> 0.1395 <![CDATA[p 21 ]]> -0.4844 <![CDATA[p 12 ]]> 0.441 <![CDATA[p 03 ]]> -0.1072 <![CDATA[p 40 ]]> -0.01839 <![CDATA[p 31 ]]> 0.1176 <![CDATA[p 22 ]]> -0.172 <![CDATA[p 13 ]]> 0.06038

[0069] Use the three error analysis indicators provided by the fitting tool - residual sum of squares (SSE), R square (R 2 ) and root mean square error (RMSE) were used to perform polynomial fitting error analysis. The analysis results are shown in Table 2.

[0070] Table 2 Results of fitting error analysis for C15-20 air compressor Analysis indicator name / symbol Analyze indicator values Residual Sum of Squares / SSE 0.00008 <![CDATA[R squared / R 2 > 0.9931 Root Mean Square Error / RMSE 0.0033

[0071] As can be seen from Table 2, compared with the R 2 0.9805, R of C15-20 air compressor 2 It is higher, indicating that the fitting results have an excellent agreement with the MAP diagram, so this modeling method can be considered accurate.

[0072] 2. Application of actual control examples based on numerical modeling

[0073] The actual example first analyzes the calculation principle of the air compressor compression ratio and speed, and then calculates the air compressor compression ratio and speed when OER=2 based on the Ballard HD6-V2 fuel cell system parameters, and then calculates the drive motor voltage value, and obtains its control result through the PID algorithm.

[0074] Calculation principle of the example:

[0075] The dynamic equation of the air compressor at work is as follows:

[0076] (3)

[0077] In the formula, is the air compressor input torque, unit: N·m; is the actual torque of the air compressor, unit: N·m; is the specific heat capacity of air at constant pressure, which is 1004 J·kg -1 ·K -1 ; and is the motor constant driving the air compressor, in units of V·rad -1 ·s and N·m·A -1 ; and are the motor mechanical efficiency and the air compressor efficiency, respectively, which are dimensionless variables; is the driving motor voltage, unit: V; is the air mass flow rate of the air compressor, unit: kg·s -1 ; is the resistance of the motor armature, unit: Ω; The temperature of air flowing out of the air compressor, unit: K; is the temperature of air entering the air compressor, which is 298.15 K; is the motor angular velocity, unit: rad·s -1 ; is the angular velocity of the air compressor impeller, in rad·s -1 ; is the air compressor rotation inertia, unit: kg·m 2 ; To obtain the air mass flow rate under actual working conditions, unit: kg·s -1 , calculated according to the following formula:

[0078] (4)

[0079] In the formula, For standard atmospheric pressure, take 101325 Pa.

[0080] In order to ensure sufficient oxygen, it is necessary to set a suitable actual excess oxygen ratio. To limit the air mass flow rate involved in the reaction:

[0081] (5)

[0082] In the formula, is the number of single cells, take 762; is the load current, unit: A; is the relative molecular mass of oxygen, taken as 0.032 kg·mol -1 .

[0083] Assuming the air compressor inlet pressure is 1 bar, the compression ratio change rate is calculated as follows:

[0084] (6)

[0085] In the formula, is the thermodynamic constant of air, which is 286.9 J·kg -1 ·K-1 ; is the volume of the gas supply pipeline, unit: m 3 ; is the air mass flow rate of the air supply pipeline, unit: kg·s -1 , calculated according to the following formula:

[0086] (7)

[0087] In the formula, is the cathode flow field pressure, unit: bar; is the air supply resistance coefficient, which is 0.3629 kg·s -1 bar -1 ; It is the outlet pipe pressure of the air compressor, unit: bar.

[0088] Combining the above formulas, we can derive the parameter calculation principle: First, set , calculated ; Then set follow The value of the change, according to formula (6), The rate of change is 0, at this time ,because , it can be determined ; According to the MAP diagram, only or Any one of the items can determine the other item, that is, determine ; Finally, when in formula (3) When the rate of change of is 0, we can get , and finally obtain .

[0089] Calculation process of the example:

[0090] Based on the calculation principle of the example, all the original parameters in the formula can be calculated , and The original parameters involved in the example are shown in Table 3.

[0091] Table 3 Original parameters of the example Parameter / Unit Numeric <![CDATA[Rotational inertia J of air compressor / kg·m 2 > <![CDATA[5×10 -5 ]]> <![CDATA[Motor constant K v / V·rad -1 ·s]]> 0.026 <![CDATA[Motor constant K t / N·m·A -1 > 0.036 <![CDATA[Drive motor efficiency η cm > 0.97 <![CDATA[Air compressor efficiency η cp > 0.74 <![CDATA[The resistance R of the motor armature cm / Ω]]> 0.01 <![CDATA[Volume V of the air supply pipeline sm / m 3 > 0.025

[0092] Based on the above original parameters, first set [0,100,300,500,700] (unit: s) for a total of 5 Step moment, The step value matrix is ​​[160, 170, 180, 190, 200] (unit: A), and we get The numerical matrix of is:

[0093] (8)

[0094] Then make Follow The specific numerical matrix is ​​[1.6, 1.7, 1.8, 1.9, 2.0] (unit: bar), then we get The numerical calculation results are:

[0095] (9)

[0096] After querying the MAP diagram, we get The result is [13347,14139,14869,15551,16200] (unit: / rad·s -1 ). Convert equation (3) to After the equation form is obtained, substituting all the values ​​into The resulting matrix is:

[0097] (10)

[0098] , and The numerical results are shown in Table 4.

[0099] Table 4 Calculation results of example parameters Load current / A <![CDATA[Air mass flow rate / (kg·s -1 )]]> Air compressor compression ratio <![CDATA[Air compressor speed / (rad·s -1 )]]> Driving motor voltage / V 160 0.087 1.84 13347 27.52 170 0.092 1.95 14139 29.16 180 0.098 2.07 14869 30.68 190 0.103 2.18 15551 32.10 200 0.109 2.30 16200 33.45

[0100] In this embodiment, the difference between the target OER and the actual OER is controlled by using a PID algorithm. , as shown in formulas (11) to (12):

[0101] (11)

[0102] (12)

[0103] Example =50, =5 and =0, then follow The control results of the air compressor speed and the drive motor voltage are shown in Figure 6 and Figure 7 .

[0104] This embodiment uses Origin and MATLAB fitting toolbox to complete the numerical modeling of the C15-20 air compressor based on the MAP diagram, and then completes the actual control example based on numerical modeling. First, based on Origin's image digitization tool, the air compressor compression ratio, air compressor speed and air mass flow values ​​are identified; then, based on the above values, the MATLAB fitting toolbox is used to give the binary quartic polynomial fitting results and error analysis of the air compressor compression ratio and air compressor speed center regularization. Based on the obtained MAP diagram fitting results, the calculation principle of the air compressor compression ratio and speed is analyzed, and then based on the Ballard HD6-V2 fuel cell system parameters, the load current is used as the external disturbance, and the air compressor compression ratio and speed when OER=2 are calculated. The value of the drive motor voltage is then obtained; finally, the control result of the drive motor voltage is obtained through the PID algorithm.

[0105] The present invention can accurately reflect the numerical relationship between the rotation speed, compression ratio and air mass flow rate of the C15-20 air compressor, improve the numerical modeling of the ROTREX C15 series air compressors, quantify the control results of the C15-20 air compressors, and expand the application scope of modeling, simulation and control of the ROTREX C15 series air compressors.

[0106] This embodiment also provides an actual control system of the excess oxygen ratio based on numerical modeling of a C15-20 air compressor, including a control module, which executes the actual control method of the excess oxygen ratio based on numerical modeling of a C15-20 air compressor.

[0107] This embodiment also provides a non-transitory computer-readable storage medium, which stores computer instructions. When the computer instructions are executed by a processor, the actual control method of the excess oxygen ratio based on numerical modeling of the C15-20 air compressor is implemented.

[0108] Ginseng Figure 8 As shown, this embodiment also provides an electronic device, including:

[0109] A memory 201 and a processor 202, wherein the memory 201 and the processor 202 are communicatively connected to each other, the memory 201 stores computer instructions, and the processor 202 executes the computer instructions to execute the actual control method of the excess oxygen ratio based on numerical modeling of a C15-20 air compressor.

[0110] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. It should be pointed out that a person skilled in the art can make several improvements and modifications without departing from the technical principles of the present invention, and these improvements and modifications should also be regarded as within the scope of protection of the present invention.

Claims

1. A practical control method of excess oxygen ratio based on numerical modeling of C15-20 air compressor, characterized in that: The steps include: Step 1, based on the image digitization tool of Origin, numerically extract the air mass flow and compression ratio at each speed in the MAP diagram of the C15-20 air compressor; Step 2, based on the extracted values, using the curve fitting toolbox of MATLAB, obtain the bivariate quartic polynomial fitting result of the air mass flow rate after regularization with respect to the compression ratio and the center of the speed of the air compressor; Step 3, based on the obtained fitting results, taking the load current step as the external disturbance, calculate the actual air compressor compression ratio, speed and driving motor voltage value when the oxygen ratio is 2; Step 4, taking the overoxygen ratio = 2 as the control target, PID control is performed on the drive motor voltage to obtain the control result of the drive motor voltage under the load current step.

2. The actual control method of excess oxygen ratio based on numerical modeling of C15-20 air compressor according to claim 1 is characterized in that: The equation of the two-variable quartic polynomial in step 2 is as follows: (1) In the formula, is the air mass flow rate of the air compressor; for The compression ratio of the air compressor after center regularization; for The air compressor speed after center regularization, unit: rpm; i and j are and The number of items; are the polynomial coefficients.

3. The actual control method of excess oxygen ratio based on numerical modeling of C15-20 air compressor according to claim 2 is characterized in that: The expression of the center regularization in step 2 is as follows: (2) In the formula, is the compression ratio; is the air compressor speed.

4. The actual control method of excess oxygen ratio based on numerical modeling of C15-20 air compressor according to claim 3 is characterized in that: The step 2 also includes: The residual sum of squares, R square and root mean square were used as error indicators to perform polynomial fitting error analysis and obtain the analysis results.

5. The actual control method of excess oxygen ratio based on numerical modeling of C15-20 air compressor according to claim 4 is characterized in that: The step 3 specifically includes: The dynamic equation of the air compressor when working is as follows: (3) In the formula, is the air compressor input torque, unit: N·m; is the actual torque of the air compressor, unit: N·m; is the specific heat capacity of air at constant pressure, which is 1004 J·kg -1 ·K -1 ; and is the motor constant driving the air compressor, in units of V·rad -1 ·s and N·m·A -1 ; and are the motor mechanical efficiency and the air compressor efficiency, respectively, which are dimensionless variables; is the driving motor voltage, unit: V; is the air mass flow rate of the air compressor, unit: kg·s -1 ; is the resistance of the motor armature, unit: Ω; The temperature of air flowing out of the air compressor, unit: K; is the temperature of air entering the air compressor, which is 298.15 K; is the motor angular velocity, unit: rad·s -1 ; is the angular velocity of the air compressor impeller, unit: rad·s -1 ; is the air compressor rotation inertia, unit: kg·m 2 ; is the air mass flow rate under actual working conditions, unit: kg·s -1 , calculated according to the following formula: (4) In the formula, For standard atmospheric pressure, take 101325 Pa; Setting the actual excess oxygen ratio Limit the air mass flow rate involved in the reaction: (5) In the formula, is the number of single cells, take 762; is the load current, unit: A; is the relative molecular mass of oxygen, taken as 0.032 kg·mol -1 ; F is Faraday constant, taken as 96485 C·mol -1 ; Assuming the air compressor inlet pressure is 1 bar, the compression ratio change rate is calculated as follows: (6) In the formula, is the thermodynamic constant of air, which is 286.9 J·kg -1 ·K -1 ; is the volume of the gas supply pipeline, unit: m 3 ; is the air mass flow rate of the air supply pipeline, unit: kg·s -1 , calculated according to the following formula: (7) In the formula, is the cathode flow field pressure, unit: bar; is the air supply resistance coefficient, which is 0.3629 kg·s -1 bar -1 ; is the outlet pipe pressure of the air compressor, unit: bar; set up , calculated ; Then set follow The value of the change, according to formula (6), The rate of change is 0, at this time ,because , in order to determine ; Determined according to MAP diagram or Any one of the ; When in formula (3) When the rate of change of is 0, we can obtain , and finally obtain .

6. The actual control method of excess oxygen ratio based on numerical modeling of C15-20 air compressor according to claim 5 is characterized in that: The step 4 specifically includes: For the target peroxygen ratio OER ref Compared with the actual peroxygen ratio The difference e(t) is controlled using the PID algorithm , as shown in formulas (11) to (12): (11) (12) Pick =50, =5 and =0, so follow change.

7. A practical control system of excess oxygen ratio based on numerical modeling of C15-20 air compressor, characterized in that: It includes a control module, which executes an actual control method of excess oxygen ratio based on numerical modeling of a C15-20 air compressor as described in any one of claims 1 to 6.

8. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium stores computer instructions, which, when executed by a processor, implement an actual control method for the excess oxygen ratio based on numerical modeling of a C15-20 air compressor as described in any one of claims 1 to 6.

9. An electronic device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the actual control method of the excess oxygen ratio based on numerical modeling of a C15-20 air compressor as described in any one of claims 1 to 6 by executing the computer instructions.

Images