CN113609792A - Multidisciplinary modeling method based on power flow - Google Patents

Multidisciplinary modeling method based on power flow Download PDF

Info

Publication number
CN113609792A
CN113609792A CN202110575533.2A CN202110575533A CN113609792A CN 113609792 A CN113609792 A CN 113609792A CN 202110575533 A CN202110575533 A CN 202110575533A CN 113609792 A CN113609792 A CN 113609792A
Authority
CN
China
Prior art keywords
model
fuel cell
power system
hydrogen fuel
resistance
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.)
Granted
Application number
CN202110575533.2A
Other languages
Chinese (zh)
Other versions
CN113609792B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202110575533.2A priority Critical patent/CN113609792B/en
Publication of CN113609792A publication Critical patent/CN113609792A/en
Application granted granted Critical
Publication of CN113609792B publication Critical patent/CN113609792B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Fuel Cell (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

The invention discloses a multidisciplinary modeling method based on power flow. The method comprises the following steps of taking a hydrogen fuel cell automobile power system as a specific research object, starting from a power system energy supply device, and establishing a hydrogen fuel cell electrochemical mechanism model; the internal electrochemical reaction kinetics and the electrochemical thermodynamic mechanism of the fuel cell are researched by constructing a voltage output model, the internal association mechanism of electrode back pressure and effective partial pressure of gas is analyzed by constructing a back pressure regulation model, the model which directly takes gas flow or effective partial pressure as output in the prior art is replaced, and the dynamic accurate regulation of the output performance of the fuel cell is realized; aiming at the energy consumption of the load of the hydrogen fuel cell automobile power system, resistance models of all parts of the load of the power system are established and fused into a dynamic load model of the power system. The invention takes power as a link, and realizes the cooperative analysis of the driving force and resistance of the hydrogen fuel cell automobile and the energy supply and consumption of a power system.

Description

Multidisciplinary modeling method based on power flow
Technical Field
The invention belongs to the field of multidisciplinary modeling of complex products, and particularly relates to a multidisciplinary modeling method based on power flow.
Background
In 2020, the quantity of motor vehicles in China is up to 3.72 hundred million, which is increased by 3.56% in 2019, wherein the new energy automobile part is increased by 29.18% in 2019, and the quantity of the motor vehicles is increased to 492 thousands. At present, related matching technologies of pure electric vehicles are mature, but the electric energy of the pure electric vehicles completely depends on the power grid supply, so the environmental protection degree of the pure electric vehicles depends on the environmental protection property of a power generation mode, but at present, China still mainly uses traditional thermal power generation, has large pollution and is limited by a lithium battery technology, and the development of the pure electric vehicles faces a bottleneck; the hydrogen fuel cell has high energy density and larger lifting space, and the hydrogen fuel cell automobile has more environmental-friendly technical advantages and wider development prospect along with the development of vehicle-mounted hydrogen supply, liquid hydrogen storage and transportation and other supporting technologies.
The insufficient performance of the hydrogen fuel cell power system seriously hinders the progress of China in manufacturing hydrogen fuel cell automobiles with practical commercial values, and the hydrogen fuel cell automobile power system has a complex structure, is wide in subject related field and has many technical difficulties, and needs to be researched by means of a digital prototype technology to shorten the research and development period and reduce the research and development cost. Through research on a modeling method of a hydrogen fuel cell automobile power system, a power system model capable of accurately and efficiently solving various characteristics of the power system can be established.
Disclosure of Invention
In order to solve the problems in the background art, the invention provides a power flow-based modeling method for a hydrogen fuel cell automobile power system. The energy supply of the power system needs to be adjusted according to energy consumption, starting from two aspects of energy supply and energy consumption, and establishing a hydrogen fuel cell model and a power system load model at the same time, so that the synergistic solution of the energy supply and the energy consumption of the power system is realized, and the energy output characteristic of the power system is truly and comprehensively reflected.
Establishing a voltage output model by researching the electrochemical reaction mechanism of the hydrogen fuel cell; a backpressure regulation model is established by researching a backpressure regulation mechanism of the fuel cell, and the two models are fused according to the internal logic of the electrochemistry of the fuel cell to obtain a hydrogen fuel cell electrochemistry mechanism model facing to the voltage output characteristic. And aiming at the energy consumption of the load part of the power system, establishing a dynamic load model of the power system according to the automobile dynamics and the automobile kinematics.
A hydrogen fuel cell automobile power system modeling method based on power flow adopts the following technical scheme:
the hydrogen fuel cell automobile power system model comprises a hydrogen fuel cell electrochemical mechanism model and a power system load model; the electrochemical mechanism model of the hydrogen fuel cell mainly comprises an electrochemical electromotive force model, an activation polarization overvoltage model, a concentration polarization overvoltage model, an ohm polarization overvoltage model and a back pressure regulation model; electrochemical mechanism model V of hydrogen fuel celloutThe calculation formula of (a) is as follows:
Vout＝ENernst-Vact-Vconc-Vohm
wherein ,ENernstAs an electrochemical electromotive force model, VactTo activate the polarization overvoltage model, VconcIs a concentration polarization overvoltage model, VohmAn ohmic polarization overvoltage model;
the power system load model mainly comprises a road resistance model, an air resistance model, an acceleration resistance model and a transmission efficiency model; power system load model FtThe calculation formula of (a) is as follows:
Ft＝Fr+Fw+Fa+Fη
wherein ,FrAs a model of road resistance, FwAs an air resistance model, FaFor the acceleration resistance model composition, FηIs a transmission efficiency module.
The electrochemical electromotive force model ENernstThe calculation formula of (a) is specifically as follows:
Figure BDA0003084220410000021
wherein Δ G is the amount of change in gibbs free energy; f is a Faraday constant; Δ S is the variation of the standard molar entropy; r is a universal gas constant;
Figure BDA0003084220410000022
is an effective partial pressure of hydrogen;
Figure BDA0003084220410000023
is the effective partial pressure of oxygen; n is the number of mobile electrons; t isstIs the fuel cell operating temperature; t isrefIs the reference temperature.
The activated polarization overvoltage model VactThe calculation formula of (a) is specifically as follows:
Figure BDA0003084220410000024
wherein ,TstIs the fuel cell operating temperature; t isrefIs a reference temperature; pcBack pressure for the cathode of the cell;
Figure BDA0003084220410000025
is the saturated water vapor pressure; i is the external circuit current density; x is the number of1,x2…x14All are electrochemical correlation coefficients.
Effective partial pressure of hydrogen
Figure BDA0003084220410000026
And effective partial pressure of oxygen
Figure BDA0003084220410000027
The method comprises the following steps of constructing a backpressure regulation model:
Figure BDA0003084220410000028
Figure BDA0003084220410000031
wherein ,PaBack pressure of the battery anode; m is1,m2…m5All are electrochemical correlation coefficients.
The concentration polarization overvoltage model VconcThe calculation formula of (a) is specifically as follows:
Figure BDA0003084220410000032
wherein R is a universal gas constant; t isstIs the fuel cell operating temperature; f is a Faraday constant; i is the current magnitude in the circuit; i islimIs the limiting current of the fuel cell.
The ohmic polarization overvoltage model VohmThe calculation formula of (a) is specifically as follows:
Figure BDA0003084220410000033
wherein I is the current magnitude in the circuit; s is the effective activation area of the proton exchange membrane; rcIs the equivalent impedance of the fuel cell external circuit; l is the proton exchange membrane thickness; gamma is the resistivity correlation coefficient of the proton exchange membrane;
wherein v is the water content of the proton exchange membrane, and the calculation formula is as follows:
Figure BDA0003084220410000034
Pcfor cell cathode back pressure, lambdacK is the water content correlation coefficient of the proton exchange membrane for the stoichiometric ratio of the cathode reactant of the battery.
The output parameter of the electrochemical mechanism model of the hydrogen fuel cell is the output power P of the hydrogen fuel celleout：
Peout＝Vout×I
The concrete construction processes of a road resistance model, an air resistance model, an acceleration resistance model and a transmission efficiency model in the power system load model are as follows:
1) road resistance model FrThe construction of (1):
1.1) carrying out moment calculation by taking the middle points of the front and rear wheels of the automobile and the road contact surface as centers, and calculating the normal reaction force of each wheel in the running process of the automobile, wherein the specific formula is as follows:
Figure BDA0003084220410000035
wherein G is the total gravity of the automobile; b1Distance from the front wheel to the center of gravity of the vehicle; b2The distance from the rear wheel to the center of gravity of the vehicle; l is the distance between the front and rear wheels; alpha is the slope angle of the road; h isgThe distance between the gravity center of the vehicle and the wheel axle; u is the vehicle speed; fzw1 and Fzw2Respectively, an aerodynamic lift force acting on the vehicle body and located above the grounding points of the left front and left rear wheels, Fzw3 and Fzw4Respectively acting on the vehicle body and positioned above the grounding points of the right front wheel and the right rear wheel; t isf1 and Tf2Rolling resistance couple moments, T, on the left front and left rear wheels, respectivelyf3 and Tf4Respectively rolling resistance couple moments on the right front wheel and the right rear wheel; t isjw1 and Tjw2The inertia resistance couple moments, T, on the left front and left rear wheels respectivelyjw3 and Tjw4The moment of inertia resistance couple on the right front wheel and the right rear wheel respectively; fz1 and Fz2Normal reaction forces of a left front wheel and a left rear wheel of the automobile are respectively; fz3、Fz4Normal reaction forces of the right front wheel and the right rear wheel are respectively; m is the total mass of the automobile;
wherein ,
Figure BDA0003084220410000041
Tf(1,2,3,4)＝Grsfcosα，
Figure BDA0003084220410000042
Figure BDA0003084220410000043
wherein ,rsIs the wheel radius; clf and CrfThe air lift coefficients of the front wheel and the rear wheel are respectively; i iswThe moment of inertia of each wheel; a is the frontal area; rhoaIs the air density; f is a rolling resistance coefficient of the rolling element,
Figure BDA0003084220410000044
f0,f1,f4is an inherent coefficient that affects the rolling resistance of a tire;
1.2) F obtained according to step 1.1)z1,Fz2,Fz3,Fz4Calculating the sum F of normal reaction forces borne by each tire of the automobilezThe concrete formula is as follows:
Fz＝Fz1+Fz2+Fz3+Fz4
＝Gcosα-ClfAρau2-CrfAρau2
1.3) calculating by the following formula to obtain a road resistance model Fr：
Figure BDA0003084220410000045
2) Air resistance model FwThe construction of (1):
Figure BDA0003084220410000046
wherein ,FwIs the air resistance; cDIs the air resistance coefficient.
3) Acceleration resistance model FaThe construction of (1):
Figure BDA0003084220410000051
wherein ,FaIs acceleration resistance; i iswThe moment of inertia of each wheel; r issIs the wheel radius;
4) transmission efficiency model FηThe construction of (1):
Figure BDA0003084220410000052
wherein T is motor torque; i.e. i0The transmission ratio of the main speed reducer is set; etaTThe transmission efficiency of the power system is improved.
5) Obtaining a load model F of the power system according to a calculation formula of each load resistance in the power systemt：
Figure BDA0003084220410000053
The power system load model FtThe output parameters are acceleration a, speed u and power consumption power of the power system
Figure BDA0003084220410000054
Wherein the acceleration a and the speed u are output parameters of the acceleration resistance model
1) The calculation formula of the acceleration a is as follows:
Figure BDA0003084220410000055
wherein a is the acceleration generated by the power system; etaTThe transmission efficiency of the power system is improved; peIs the driving power of the power system, i.e. is equal to the output power, P, of the hydrogen fuel celle＝Peout＝Vout×I；
2) The velocity u is calculated as follows:
Figure BDA0003084220410000056
wherein ,t0Is the initial time; u. of0At an initial speed, i.e. t0The vehicle running speed at the moment.
3) Power consumption of power system
Figure BDA0003084220410000057
The calculation formula of (a) is as follows:
Figure BDA0003084220410000058
the invention has the beneficial effects that:
1. compared with the traditional model in which the reaction gas flow or the effective partial pressure is used as input, the electrochemical mechanism model part of the hydrogen fuel cell uses the back pressure of the anode and the cathode as input, so that the model is more in line with the actual situation in the production and application of the fuel cell and obtains a more accurate solving result.
2. Through the established electrochemical mechanism model of the hydrogen fuel cell and the dynamic load model of the power system, the collaborative analysis of the driving force and resistance of the hydrogen fuel cell automobile and the energy supply and consumption of the power system by taking the power as a link is realized.
Drawings
FIG. 1 is a schematic diagram of a model of the electrochemical mechanism of a hydrogen fuel cell;
FIG. 2 is a model diagram of the electrochemical Enernst electromotive force portion.
FIG. 3 is a schematic diagram of a portion of the activated polarization overvoltage.
FIG. 4 is a schematic diagram of a concentration polarization overvoltage portion.
FIG. 5 is a model diagram of the ohmic polarization overvoltage portion.
FIG. 6 is a graph of a voltage output model.
Fig. 7 is a model of hydrogen fuel cell back pressure regulation.
FIG. 8 is a model of dynamic loading of the powertrain.
FIG. 9 is a model diagram of a road resistance segment.
FIG. 10 is a schematic view of a model of an air resistance section.
FIG. 11 is a model diagram of an acceleration resistance portion.
Detailed Description
The invention is further illustrated by the following figures and examples.
The invention comprises a hydrogen fuel cell electrochemical mechanism model and a dynamic system dynamic load model. The electrochemical mechanism model of the hydrogen fuel cell is used for modeling a power system energy supply device, and the dynamic load model of the power system is used for modeling the energy consumption of the load of the power system. By the aid of the hydrogen fuel cell electrochemical mechanism model and the power system dynamic load model, the driving force and resistance of the hydrogen fuel cell automobile and the energy supply and consumption of the power system are cooperatively analyzed by taking power as a link.
As shown in fig. 1, the electrochemical mechanism model of the hydrogen fuel cell includes two parts, a voltage output model and a back pressure regulation model. The voltage output model takes the operation temperature, the current, the effective partial pressure of hydrogen, the effective partial pressure of oxygen and the cathode back pressure of the fuel cell as input parameters and takes the output voltage of the fuel cell as an output parameter. Compared with the existing hydrogen fuel cell related model, the hydrogen fuel cell voltage output model disclosed by the invention has the advantages that the factors such as the effective partial pressure of the reaction gas, the cell operation temperature, the chemical metering ratio of the reaction gas, the humidity in the cell and the like are calculated in detail, the influence factors of the fuel cell are considered more comprehensively, and the solving result of the output voltage is more accurate.
As shown in FIG. 6, the voltage output model includes an electrochemical electromotive force model Enernst, an activated polarization overvoltage model VactConcentration polarization overvoltage model VconcAnd ohm polarization overvoltage model Vohm：
1) The electrochemical electromotive force model Enernst is shown in fig. 2:
Figure BDA0003084220410000071
the content of Func1 in FIG. 2 is
Figure BDA0003084220410000072
2) Activated polarization overvoltage model VactAs shown in fig. 3:
Figure BDA0003084220410000073
in FIG. 3, Func1 isx2(Tst-Tref) (ii) a Func2 has the content of
Figure BDA0003084220410000074
Func3 has the content of
Figure BDA0003084220410000075
Func4 has the content of
Figure BDA0003084220410000076
Func5 has the content of
Figure BDA0003084220410000077
Func6 has the content of
Figure BDA0003084220410000078
Func7 has the content of
Figure BDA0003084220410000079
Func8 has the content of
Figure BDA00030842204100000710
wherein ,x1,x2…x14Has an actual value of x1＝0.279；x2＝0.00085；x3＝0.00004308；x4＝ 1.01325；x5＝0.5；x6＝8.63811；x7＝0.5736；x8＝0.00058；x9＝0.00018；x10＝ 0.166；x11＝0.1173；x12＝0.01618；x13＝0.00001618；x14＝10。
3) Concentration polarization overvoltage model VconcAs shown in fig. 4:
Figure BDA00030842204100000711
in FIG. 4, Func1 is
Figure BDA00030842204100000712
4) Ohm polarization overvoltage model VohmAs shown in fig. 5:
Figure BDA00030842204100000713
wherein v is the water content of the proton exchange membrane, and the calculation formula is as follows:
Figure BDA00030842204100000714
wherein, the actual value of k is 6.3348.
In FIG. 5, Func1 is
Figure BDA00030842204100000715
Func2 has the content of
Figure BDA00030842204100000716
Func3 has the content of
Figure BDA00030842204100000717
Func4 has the content of
Figure BDA00030842204100000718
Func5 has the content of
Figure BDA00030842204100000719
Func6 has the content of
Figure BDA0003084220410000081
Func7 has the content of
Figure BDA0003084220410000082
Func8 has the content of
Figure BDA0003084220410000083
Func9 has the content of
Figure BDA0003084220410000084
Func10 has the content of
Figure BDA0003084220410000085
The backpressure regulation model takes the backpressure of the anode and the cathode, the current and the operation temperature of the fuel cell as input, takes the effective partial pressure of hydrogen and oxygen as output, analyzes the internal association mechanism of the electrode backpressure and the effective partial pressure of gas, replaces the prior model which directly takes the gas flow or the effective partial pressure as output, better accords with the actual application scene of the hydrogen fuel cell, and realizes more accurate regulation of the output performance of the fuel cell. Fig. 7 shows a back pressure regulation model:
Figure BDA0003084220410000086
Figure BDA0003084220410000087
in FIG. 7, Func1 is
Figure BDA0003084220410000088
Func2 has the content of
Figure BDA0003084220410000089
Func3 has the content of
Figure BDA00030842204100000810
Func4 has the content of
Figure BDA00030842204100000811
Func5 has the content of
Figure BDA00030842204100000812
wherein ,m1,m2…m5The values of (A) are as follows:
m1＝3.7619；m2＝0.291；m3＝0.832；m4＝1.635；m5＝1.334。
as shown in FIG. 8, the dynamic load model of the power system includes three parts, namely a road resistance model, an air resistance model and an acceleration resistance model. The driving power, the initial speed, the road parameters and the like of the power system are used as input, and the acceleration, the current speed and the resistance power are used as output. And (3) independently and completely modeling each part of resistance of the hydrogen fuel cell automobile power system load through the relation between the driving power of the power system and the power consumed by each part of resistance.
Func1 in FIG. 8 is included as a transmission efficiency model
Figure BDA00030842204100000813
And (4) partial calculation.
The road resistance model is shown in fig. 9:
Figure BDA0003084220410000091
in FIG. 9, Func1 is Gcos alpha-CLfAρau2-CrfAρau2(ii) a Func2 has the content of
Figure BDA0003084220410000092
wherein ,f0,f1,f4The actual values of (a) are as follows: f. of0＝0.0076；f1＝0.000056；f4＝0
Air resistance model FwAs shown in fig. 10:
Figure BDA0003084220410000093
acceleration resistance model FaAs shown in fig. 11:
Figure BDA0003084220410000094
acceleration resistorThe output parameter of the force model is Fa、a、u。
The invention provides a multi-model integrated modeling method of a power system, which integrates a power system model according to the logic relation between a hydrogen fuel cell electrochemical mechanism model and a power system dynamic load model and the flow direction of mutual feedback information, and describes a hydrogen fuel cell automobile simulation calculation process. When the simulation of the hydrogen fuel cell automobile power system is carried out, the output power can be obtained by multiplying the output voltage of the hydrogen fuel cell electrochemical mechanism model and the circuit current. In a power system independently powered by a hydrogen fuel cell, the output power of the hydrogen fuel cell is used as an input parameter of a dynamic load model of the power system, and simulation analysis of the power system can be performed. The powertrain dynamic load model may calculate the current power consumption and power demand of the powertrain. The electrochemical mechanism model of the hydrogen fuel cell needs to adjust the output power based on the above, transmit the output power of the hydrogen fuel cell to a load model of a power system, and calculate the current acceleration and speed of the automobile by combining the driving power and the resistance power, thereby forming iterative operation.

Claims (9)

1. A hydrogen fuel cell automobile power system modeling method based on power flow is characterized in that a hydrogen fuel cell automobile power system model comprises a hydrogen fuel cell electrochemical mechanism model and a power system load model;
the electrochemical mechanism model of the hydrogen fuel cell mainly comprises an electrochemical electromotive force model, an activation polarization overvoltage model, a concentration polarization overvoltage model, an ohm polarization overvoltage model and a back pressure regulation model; electrochemical mechanism model V of hydrogen fuel celloutThe calculation formula of (a) is as follows:
Vout＝ENernst-Vact-Vconc-Vohm
wherein ,ENernstAs an electrochemical electromotive force model, VactTo activate the polarization overvoltage model, VconcIs a concentration polarization overvoltage model, VohmAn ohmic polarization overvoltage model;
power systemThe load model mainly comprises a road resistance model, an air resistance model, an acceleration resistance model and a transmission efficiency model; power system load model FtThe calculation formula of (a) is as follows:
Ft＝Fr+Fw+Fa+Fη
wherein ,FrAs a model of road resistance, FwAs an air resistance model, FaFor the acceleration resistance model composition, FηIs a transmission efficiency module.
2. The modeling method for a power-flow-based hydrogen fuel cell automotive power system of claim 1, wherein the electrochemical electromotive force model ENernstThe calculation formula of (a) is specifically as follows:
Figure FDA0003084220400000011
wherein Δ G is the amount of change in gibbs free energy; f is a Faraday constant; Δ S is the variation of the standard molar entropy; r is a universal gas constant;
Figure FDA0003084220400000012
is an effective partial pressure of hydrogen;
Figure FDA0003084220400000013
is the effective partial pressure of oxygen; n is the number of mobile electrons; t isstIs the fuel cell operating temperature; t isrefIs the reference temperature.
3. The power-flow-based hydrogen fuel cell automotive power system modeling method of claim 1, wherein the activated polarization overvoltage model VactThe calculation formula of (a) is specifically as follows:
Figure FDA0003084220400000014
wherein ,TstIs the fuel cell operating temperature; t isrefIs a reference temperature; pcBack pressure for the cathode of the cell;
Figure FDA0003084220400000015
is the saturated water vapor pressure; i is the external circuit current density; x is the number of1,x2…x14All are electrochemical correlation coefficients.
4. A power flow based hydrogen fuel cell automotive power system modeling method as claimed in claim 2 or 3, characterized by effective partial pressure of hydrogen gas
Figure FDA0003084220400000021
And effective partial pressure of oxygen
Figure FDA0003084220400000022
The method comprises the following steps of constructing a backpressure regulation model:
Figure FDA0003084220400000023
Figure FDA0003084220400000024
wherein ,PaBack pressure of the battery anode; m is1,m2…m5All are electrochemical correlation coefficients.
5. The power flow-based hydrogen fuel cell automotive power system modeling method of claim 1, wherein the concentration polarization overvoltage model VconcThe calculation formula of (a) is specifically as follows:
Figure FDA0003084220400000025
wherein R is a universal gas constant; t isstIs the fuel cell operating temperature; f is a Faraday constant; i is the current magnitude in the circuit; i islimIs the limiting current of the fuel cell.
6. The modeling method for a power flow based hydrogen fuel cell automotive power system as claimed in claim 1, wherein said ohmic polarization overvoltage model VohmThe calculation formula of (a) is specifically as follows:
Figure FDA0003084220400000026
wherein I is the current magnitude in the circuit; s is the effective activation area of the proton exchange membrane; rcIs the fuel cell external circuit equivalent impedance; l is the proton exchange membrane thickness; gamma is the resistivity correlation coefficient of the proton exchange membrane;
wherein v is the water content of the proton exchange membrane, and the calculation formula is as follows:
Figure FDA0003084220400000027
Pcfor cell cathode back pressure, lambdacK is the water content correlation coefficient of the proton exchange membrane.
7. The modeling method for hydrogen fuel cell automobile power system based on power flow as claimed in claim 1, wherein the output parameter of the electrochemical mechanism model of the hydrogen fuel cell is the output power P of the hydrogen fuel celleout：
Peout＝Vout×I
8. The modeling method for the power system of the hydrogen fuel cell automobile based on power flow as claimed in claim 1, wherein the road resistance model, the air resistance model, the acceleration resistance model and the transmission efficiency model in the power system load model are constructed by the following steps:
1) road resistance model FrThe construction of (1):
1.1) carrying out moment calculation by taking the middle points of the front and rear wheels of the automobile and the road contact surface as centers, and calculating the normal reaction force of each wheel in the running process of the automobile, wherein the specific formula is as follows:
Figure FDA0003084220400000031
wherein G is the total gravity of the automobile; b1Distance from the front wheel to the center of gravity of the vehicle; b2The distance from the rear wheel to the center of gravity of the vehicle; l is the distance between the front and rear wheels; alpha is the slope angle of the road; h isgThe distance between the gravity center of the vehicle and the wheel axle; u is the vehicle speed; fzw1 and Fzw2Respectively, an aerodynamic lift force acting on the vehicle body and located above the grounding points of the left front wheel and the left rear wheel, Fzw3 and Fzw4Respectively acting on the vehicle body and positioned above the grounding points of the right front wheel and the right rear wheel; t isf1 and Tf2Rolling resistance couple moments, T, on the left front and left rear wheels, respectivelyf3 and Tf4Respectively rolling resistance couple moments on the right front wheel and the right rear wheel; t isjw1 and Tjw2The inertia resistance couple moments, T, on the left front and left rear wheels respectivelyjw3 and Tjw4The inertia resistance couple moments on the right front wheel and the right rear wheel respectively; fz1 and Fz2Normal reaction forces of a left front wheel and a left rear wheel of the automobile are respectively; fz3、Fz4Normal reaction forces of the right front wheel and the right rear wheel are respectively; m is the total mass of the automobile;
wherein ,
Figure FDA0003084220400000032
Tf(1,2,3,4)＝Grsfcosα，
Figure FDA0003084220400000033
Figure FDA0003084220400000034
wherein ,rsIs the wheel radius; clf and CrfThe air lift coefficients of the front wheel and the rear wheel are respectively; i iswThe moment of inertia of each wheel; a is the frontal area; rhoaIs the air density; f is a rolling resistance coefficient of the rolling element,
Figure FDA0003084220400000035
f0,f1,f4is an inherent coefficient that affects the rolling resistance of a tire;
1.2) F obtained according to step 1.1)z1,Fz2,Fz3,Fz4Calculating the sum F of normal reaction forces borne by each tire of the automobilezThe concrete formula is as follows:
Fz＝Fz1+Fz2+Fz3+Fz4
＝Gcosα-ClfAρau2-CrfAρau2
1.3) calculating by the following formula to obtain a road resistance model Fr：
Figure FDA0003084220400000041
2) Air resistance model FwThe construction of (1):
Figure FDA0003084220400000042
wherein ,FwIs the air resistance; cDIs the air resistance coefficient.
3) Acceleration resistance model FaThe construction of (1):
Figure FDA0003084220400000043
wherein ,FaIs acceleration resistance; i iswThe moment of inertia of each wheel; r issIs the wheel radius;
4) transmission efficiency model FηThe construction of (1):
Figure FDA0003084220400000044
wherein T is motor torque; i.e. i0The transmission ratio of the main speed reducer is set; etaTThe transmission efficiency of the power system is improved.
5) Obtaining a load model F of the power system according to a calculation formula of each load resistance in the power systemt：
Figure FDA0003084220400000045
9. The power flow-based hydrogen fuel cell automotive power system modeling method of claim 8, characterized in that the power system load model FtThe output parameters of are acceleration a, speed u and power consumption of the power system
Figure FDA0003084220400000046
Wherein the acceleration a and the speed u are output parameters of the acceleration resistance model:
1) the calculation formula of the acceleration a is as follows:
Figure FDA0003084220400000047
wherein a is the acceleration generated by the power system; etaTThe transmission efficiency of the power system is improved; peFor driving power of the power system, i.e. equal to hydrogen fuelOutput power of the battery, Pe＝Peout＝Vout×I；
2) The velocity u is calculated as follows:
Figure FDA0003084220400000051
wherein ,t0Is the initial time; u. of0At an initial speed, i.e. t0The vehicle running speed at the moment.
3) Power consumption of power system
Figure FDA0003084220400000052
The calculation formula of (a) is as follows:
Figure FDA0003084220400000053
CN202110575533.2A 2021-05-26 2021-05-26 Multidisciplinary modeling method based on power flow Active CN113609792B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110575533.2A CN113609792B (en) 2021-05-26 2021-05-26 Multidisciplinary modeling method based on power flow

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110575533.2A CN113609792B (en) 2021-05-26 2021-05-26 Multidisciplinary modeling method based on power flow

Publications (2)

Publication Number Publication Date
CN113609792A true CN113609792A (en) 2021-11-05
CN113609792B CN113609792B (en) 2023-08-04

Family

ID=78336471

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110575533.2A Active CN113609792B (en) 2021-05-26 2021-05-26 Multidisciplinary modeling method based on power flow

Country Status (1)

Country Link
CN (1) CN113609792B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117374327A (en) * 2023-12-07 2024-01-09 北京理工大学 Fuel cell system safety domain modeling method, system and equipment

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR19980017054A (en) * 1996-08-30 1998-06-05 박병재 Prediction and extension device of electric vehicle and its method
CN108303875A (en) * 2017-12-31 2018-07-20 湖南沃森电气科技有限公司 A kind of control method of electric power load for testing simulator and its system
CN109657348A (en) * 2018-12-18 2019-04-19 安徽江淮汽车集团股份有限公司 One proton exchanging film fuel battery modeling method
US20190181465A1 (en) * 2017-12-12 2019-06-13 Hyundai Motor Company Fuel cell control method and control system
CN112349936A (en) * 2020-10-27 2021-02-09 中船动力有限公司 Marine fuel cell modeling and generalized predictive control method based on system identification
CN112820906A (en) * 2021-01-15 2021-05-18 湖南理工学院 Comprehensive evaluation method for thermodynamic performance of vehicle fuel cell under dynamic working condition

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR19980017054A (en) * 1996-08-30 1998-06-05 박병재 Prediction and extension device of electric vehicle and its method
US20190181465A1 (en) * 2017-12-12 2019-06-13 Hyundai Motor Company Fuel cell control method and control system
CN108303875A (en) * 2017-12-31 2018-07-20 湖南沃森电气科技有限公司 A kind of control method of electric power load for testing simulator and its system
CN109657348A (en) * 2018-12-18 2019-04-19 安徽江淮汽车集团股份有限公司 One proton exchanging film fuel battery modeling method
CN112349936A (en) * 2020-10-27 2021-02-09 中船动力有限公司 Marine fuel cell modeling and generalized predictive control method based on system identification
CN112820906A (en) * 2021-01-15 2021-05-18 湖南理工学院 Comprehensive evaluation method for thermodynamic performance of vehicle fuel cell under dynamic working condition

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周鹏: "燃料电池系统建模与仿真研究", 万方中国学位论文数据库 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117374327A (en) * 2023-12-07 2024-01-09 北京理工大学 Fuel cell system safety domain modeling method, system and equipment
CN117374327B (en) * 2023-12-07 2024-02-23 北京理工大学 Fuel cell system safety domain modeling method, system and equipment

Also Published As

Publication number Publication date
CN113609792B (en) 2023-08-04

Similar Documents

Publication Publication Date Title
Eren et al. A fuzzy logic based supervisory controller for an FC/UC hybrid vehicular power system
Basma et al. Fuel cell electric tractor-trailers: Technology overview and fuel economy
Turner et al. Modeling a PEM fuel cell for use in a hybrid electric vehicle
Robledo et al. Hydrogen fuel cell scooter with plug-out features for combined transport and residential power generation
CN113609792B (en) Multidisciplinary modeling method based on power flow
CN114347866B (en) VMD processing-based fuel cell energy management method
Karaoğlan et al. Investigation of the effects of battery types and power management algorithms on drive cycle simulation for a range-extended electric vehicle powertrain
CN113815437B (en) Predictive energy management method for fuel cell hybrid electric vehicle
CN111753402B (en) Modeling method for auxiliary power source lithium battery of fuel cell automobile
Rizzo et al. Converting conventional cars in mild hybrid solar vehicles
Zhang et al. Optimal sizing of fuel cell electric vehicle powertrain considering multiple objectives
Song et al. Study on the fuel economy of fuel cell electric vehicle based on rule-based energy management strategies
BOYACIOĞLU et al. Modeling and Control of a PEM Fuel Cell Hybrid Energy System Used in a Vehicle with Fuzzy Logic Method
Mourad A proposed fuel cell vehicle for reducing CO2 emissions and its contribution to reducing greenhouse gas emissions
Xiao et al. An optimized energy management strategy for fuel cell hybrid vehicles
Wang et al. Prediction and Analysis of Energy Saving Potential for Fuel Cell Vehicles
Bai et al. Research on Key Technologies of New Energy Fuel Cell Electric Vehicle
Parreira et al. Hydrogen as the fuel of the future and its application in electric mobility
Bernabé et al. Development of a Hybrid Vehicular Prototype with H2/air PEMFC Stack and Li-Ion Battery Module
Sun et al. Modeling and Simulation of PEM Fuel Cell/Battery Hybrid Vehicle
Andriyashin et al. Fuel cells as a future of automotive industry: A review of current developments
Alam et al. Realization of hybrid-electric powertrain system for a three wheeler auto taxi
Ommi et al. Evaluation of Model and Performance of Fuel Cell Hybrid Electric Vehicle in Different Drive Cycles
Zheng et al. Strategies and evaluation of fuel economy in fuel cell hybrid vehicles
Jee et al. Modeling and Simulation for PEMFC/Solar Panel Hybrid Vehicle With Solar Water Electrolysis System

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