Switched Energy Management Strategy for Fuel Cell Hybrid Vehicle Based on Switch Network

: Environmentally friendly and pollution-free fuel cell / lithium battery hybrid vehicles have received the attention of the community in recent years. It is imperative for fuel cell / lithium battery hybrid vehicles to use the energy management strategy (EMS) to distribute the output power of each power source to improve fuel economy and system life. In practical application, inconsistency of battery pack will lead to security hazard and capacity degradation. However, few EMS take the inconsistency of battery pack into account. Also, the current battery equalization strategy rarely discusses how to perform the equilibrium process while meeting the power demand of vehicle. To solve these issues, a novel equalization energy management strategy (EEMS) based on the switch network is proposed at ﬁrst. Then, a switched energy management strategy (SEMS) that switches between the EEMS and the equivalent consumption minimization strategy (ECMS) is proposed and implemented in the fuel cell / lithium battery hybrid system to validate its e ﬀ ectiveness. The results show that the proposed SEMS can ameliorate the inconsistency of series lithium battery pack while meeting the power demand of vehicle’s normal operation. It can improve the safety and durability of the system and reduce the equalization time. Besides, it has good expansibility and no energy waste.


Introduction
In recent years, the world's environmental pollution and energy shortages have become increasingly serious. Not only in China, but also worldwide, the carbon emissions of metallurgical industry, shipping, aerospace and other industries have attracted the attention of the community [1][2][3][4][5]. Many regions and industries have released relevant regulations to limit carbon emissions. So, the development of traditional fuel vehicles has been limited [6][7][8]. To reduce exhaust emissions and energy consumption, it has become an irresistible trend for the automotive industry to replace fossil fuels by using electricity and other clean energy [9][10][11]. Batteries have high energy density, fast dynamic response, and no pollution. They are widely used in electric vehicles. However, batteries also have some shortcomings, such as low power density, long charging time, and short life [12,13]. Fuel cells with high efficiency and no pollution are ideal fossil energy alternatives. However, fuel cell stacks require flow supply, diffusion and catalysis, which lead to slow dynamic response. Hybrid vehicles can combine the advantages of fuel cells and lithium batteries, and have high development prospects [14,15].
EMS is the basis for the good performance of hybrid vehicles. The key issue is the energy distribution in hybrid systems under various operating conditions [16,17]. EMS mainly includes rule-based control strategies and optimization-based control strategies [18][19][20].
SOC of battery pack can be controlled to a specific value. The real-time simulation of a typical vehicle driving cycle verifies that the SEMS can balance the series lithium battery pack during the normal operation of the vehicle by controlling the switch state to improve the safety and durability of the system. Besides, it has good expansibility and no energy waste.
The rest of the paper is organized as follows. Section 2 includes three parts: Part 2.1 describes the modeling of the proton exchange membrane fuel cell. Part 2.2 describes the methods and results of the battery test, and the modeling of series lithium battery pack is described. Part 2.3 describes the modeling of the DC/DC converter. Section 3 includes two parts: Part 3.1 first introduces the circuit structure of the switch network, and then describes the specific rules of the EEMS. Part 3.2 describes the algorithm of the ECMS. Section 4 gives the simulation results and analysis. The conclusion is drawn in Section 5.

Hybrid System Modeling
Fuel cell/lithium battery hybrid system can make full use of their respective advantages. The main power is provided by the fuel cell, and the lithium battery pack is used as an auxiliary power source to provide partial power and recycle braking energy during the operation of the vehicle. The structure of the hybrid system is shown in Figure 1. The fuel cell is connected to the DC bus through a unidirectional DC/DC converter, and the lithium battery pack is connected to the bus through a bidirectional DC/DC converter to provide auxiliary power.
Energies 2020, 13, x FOR PEER REVIEW 3 of 22 lithium battery pack can be improved under the premise of ensuring the vehicle's normal operation, and the SOC of battery pack can be controlled to a specific value. The real-time simulation of a typical vehicle driving cycle verifies that the SEMS can balance the series lithium battery pack during the normal operation of the vehicle by controlling the switch state to improve the safety and durability of the system. Besides, it has good expansibility and no energy waste. The rest of the paper is organized as follows. Section 2 includes three parts: Part 2.1 describes the modeling of the proton exchange membrane fuel cell. Part 2.2 describes the methods and results of the battery test, and the modeling of series lithium battery pack is described. Part 2.3 describes the modeling of the DC/DC converter. Section 3 includes two parts: Part 3.1 first introduces the circuit structure of the switch network, and then describes the specific rules of the EEMS. Part 3.2 describes the algorithm of the ECMS. Section 4 gives the simulation results and analysis. The conclusion is drawn in Section 5.

Hybrid System Modeling
Fuel cell/lithium battery hybrid system can make full use of their respective advantages. The main power is provided by the fuel cell, and the lithium battery pack is used as an auxiliary power source to provide partial power and recycle braking energy during the operation of the vehicle. The structure of the hybrid system is shown in Figure 1. The fuel cell is connected to the DC bus through a unidirectional DC/DC converter, and the lithium battery pack is connected to the bus through a bidirectional DC/DC converter to provide auxiliary power.

Proton Exchange Membrane Fuel Cell Modeling
The single cell voltage of proton exchange membrane fuel cell can be calculated by open circuit voltage, activation loss, ohmic loss and concentration loss [38]: The open circuit voltage is calculated according to the change of Gibbs free energy in the chemical reaction: where, is the fuel cell stack temperature. 2 and 2 are hydrogen partial pressure and oxygen partial pressure, respectively. They are expressed in atm.
The activation loss is a result of the need to cause electron transfer and to break and form chemical bonds in the anode and cathode [39]. Part of the energy is lost in driving the chemical reaction that transfers the electrons to and from the electrodes [40]. It can be described as follows:

Proton Exchange Membrane Fuel Cell Modeling
The single cell voltage of proton exchange membrane fuel cell can be calculated by open circuit voltage, activation loss, ohmic loss and concentration loss [38]: The open circuit voltage is calculated according to the change of Gibbs free energy in the chemical reaction: where, T f c is the fuel cell stack temperature. p H 2 and p O 2 are hydrogen partial pressure and oxygen partial pressure, respectively. They are expressed in atm. The activation loss is a result of the need to cause electron transfer and to break and form chemical bonds in the anode and cathode [39]. Part of the energy is lost in driving the chemical reaction that transfers the electrons to and from the electrodes [40]. It can be described as follows: where, c 1 is a constant derived from experience, and the values of v 0 and v a can be calculated according to the stack temperature, cathode pressure, and oxygen partial pressure: where, p ca , p ca , and p ca are expressed in bar. i is the stack current density, which can be obtained by stack current I f c and battery activation area A: the saturation pressure p sat is calculated as follows: the ohmic loss is caused by the resistance of the electrode material and the polymer membrane, which can be expressed as follows: where, r ohm is the internal electrical resistance which has the unit of Ω·cm 2 . t m is the thickness of the membrane and σ m is the membrane conductivity, which is a function of the membrane water content λ m and the stack temperature: where, b 11 , b 12 and b 2 are constants, which are generally obtained by experience. The concentration loss is caused by the change in the reactant concentration on the electrode surface and can be described as follows: where, i max , c 2 and c 3 are constants that depend on the temperature and the reactant partial pressure and can be determined empirically.
The fuel cell stack consists of many single cells in series, so the total voltage output is as follows:   Figure 2 shows the curve of the voltage and power of the fuel cell changing with the current density.
Energies 2020, 13, x FOR PEER REVIEW 5 of 22 Figure 2 shows the curve of the voltage and power of the fuel cell changing with the current density. .

Lithium Battery Pack Modeling
Six lithium batteries with different aging degrees were tested by the Arbin battery test bench (Arbin BT-5HC-5V/100A) as shown in Figure 3. The test bench had 3 scales (100 A/10 A/1 A) and its operating voltage was 0-5 V. The resolution of the current and voltage was ±0.5‰. The lithium battery was placed in a temperature chamber at the temperature of 298.15 K. The computer is used to load experiment procedures and collect data. Through the open circuit voltage test and hybrid pulse power characteristic (HPPC) test, the capacity test and parameter identification, the data of lithium battery can be obtained from Guo and Xiang's experiment in Ref. [41][42][43][44]. The discharge rate of the lithium battery capacity test is 1 C, and the capacity of test batteries is shown in Table 2.

Lithium Battery Pack Modeling
Six lithium batteries with different aging degrees were tested by the Arbin battery test bench (Arbin BT-5HC-5V/100A) as shown in Figure 3. The test bench had 3 scales (100 A/10 A/1 A) and its operating voltage was 0-5 V. The resolution of the current and voltage was ±0.5% . The lithium battery was placed in a temperature chamber at the temperature of 298.15 K. The computer is used to load experiment procedures and collect data. Through the open circuit voltage test and hybrid pulse power characteristic (HPPC) test, the capacity test and parameter identification, the data of lithium battery can be obtained from Guo and Xiang's experiment in Ref. [41][42][43][44]. . Figure 2. The voltage and power of fuel cell stack.

Lithium Battery Pack Modeling
Six lithium batteries with different aging degrees were tested by the Arbin battery test bench (Arbin BT-5HC-5V/100A) as shown in Figure 3. The test bench had 3 scales (100 A/10 A/1 A) and its operating voltage was 0-5 V. The resolution of the current and voltage was ±0.5‰. The lithium battery was placed in a temperature chamber at the temperature of 298.15 K. The computer is used to load experiment procedures and collect data. Through the open circuit voltage test and hybrid pulse power characteristic (HPPC) test, the capacity test and parameter identification, the data of lithium battery can be obtained from Guo and Xiang's experiment in Ref. [41][42][43][44]. The discharge rate of the lithium battery capacity test is 1 C, and the capacity of test batteries is shown in Table 2.   The discharge rate of the lithium battery capacity test is 1 C, and the capacity of test batteries is shown in Table 2. The open circuit voltage test is based on No.1 battery. Figure 4 shows the comparison between the measured data and fitted curve.  Figure 4 shows the comparison between the measured data and fitted curve.
the root means square error of the fitted curve and the measured data is 0.006, indicating that the fitted curve has high accuracy.
HPPC test is intended to determine dynamic power capacity, resistance over the battery's useable voltage range using a test profile that incorporates both discharge and regen pulses. The test profile can be seen in Ref. [45] of Jeffrey. The charge and discharge resistance of the lithium battery was measured by HPPC test. The resistance of No.1 battery is shown in Figure 5:  Here, the fitted curve formula is: equation/formula the root means square error of the fitted curve and the measured data is 0.006, indicating that the fitted curve has high accuracy. HPPC test is intended to determine dynamic power capacity, resistance over the battery's useable voltage range using a test profile that incorporates both discharge and regen pulses. The test profile can be seen in Ref. [45] of Jeffrey. The charge and discharge resistance of the lithium battery was measured by HPPC test. The resistance of No.1 battery is shown in Figure 5: Energies 2020, 13, 247 7 of 23 the root means square error of the fitted curve and the measured data is 0.006, indicating that the fitted curve has high accuracy.
HPPC test is intended to determine dynamic power capacity, resistance over the battery's useable voltage range using a test profile that incorporates both discharge and regen pulses. The test profile can be seen in Ref. [45] of Jeffrey. The charge and discharge resistance of the lithium battery was measured by HPPC test. The resistance of No.1 battery is shown in Figure 5: According the switch network state in the Section 3, the battery in the circuit may be in two states: connected or disconnected. Therefore, the voltage and resistance of each battery also have two states (ignoring the resistance of the wire): According the switch network state in the Section 3, the battery in the circuit may be in two states: connected or disconnected. Therefore, the voltage and resistance of each battery also have two states (ignoring the resistance of the wire): The equivalent model of the lithium battery pack is the RINT model. According to Ohm's law, the output voltage V pack of the series lithium battery pack is as follows: where, I bat is the current of the batteries. V oc_pack is the sum of the open circuit voltage of all connected batteries. R pack is the sum of the resistance of all connected batteries. The relationship between the current and the power can be derived from Equation (16): where, P pack is the charge or discharge power of the series lithium battery pack. For a single lithium battery, SOC can be calculated by current integration method: where, SOC t 0 is the initial value of SOC and Q is the capacity of the battery. Table 3 shows the parameter of single lithium battery used in the simulation:

DC/DC Converter Modeling
The voltages of the fuel cell and the lithium battery pack are matched to the bus voltage by unidirectional DC/DC converter and bidirectional DC/DC converter. When the system is in steady state, the function of DC/DC converter can be described by the following equations: where, V in and I in are the input voltage and current of the converter respectively; V out and I out are the output voltage and current of the converter respectively; η dc is the efficiency of the converter; the output voltage can be controlled by the duty cycle D.

Switched Energy Management Strategy
The block diagram of the SEMS is shown in Figure 6. The SEMS includes two energy management strategies. The SEMS performs real-time switching between EEMS and ECMS by comparing the standard deviation of the lithium battery SOC with a threshold.  The flow of the SEMS is: According to the power demand and the SOC status of the current lithium battery pack, the SEMS switches to EEMS mode or ECMS mode to obtain the reference output power of the lithium battery pack _ and the switch network state. The series lithium battery model calculates the reference output current of the single lithium battery  (22). The fuel cell system power is calculated based on the actual output power of the lithium battery pack and the power demand according to Equation (20) and the constraint Equation (21).

Equalization Energy Management Strategy Based on Switch Network
The switch network [37] shown in Figure 7 can control the number of working batteries by controlling the switches of branch and bypass. For example, when the switch 11 is closed and 12 is opened at the same time, the battery 1 can be connected to the circuit. On the contrary, the battery 1 can be disconnected. According to the SOC state of each lithium battery and the current power demand, the EEMS performs power distribution, and the switch controller performs real-time control When the standard deviation is greater than the threshold, it indicates that the lithium batteries are not consistent. The SEMS works in EEMS mode. In EEMS mode, the switch network controller can control whether the lithium battery is connected to the circuit by controlling the state of the switch network. When the standard deviation is less than the threshold, it indicates that the lithium batteries are consistent. The SEMS works in ECMS mode. In ECMS mode, all lithium batteries are connected to the circuit.
The flow of the SEMS is: According to the power demand and the SOC status of the current lithium battery pack, the SEMS switches to EEMS mode or ECMS mode to obtain the reference output power of the lithium battery pack P pack_re f and the switch network state. The series lithium battery model calculates the reference output current of the single lithium battery I bat_re f according to Formula (16) by reference output power P pack_re f , switch network state, resistance and open circuit voltage of each single lithium battery. Each single lithium battery model updates the SOC status and output Energies 2020, 13, 247 9 of 23 power according to the reference output current I bat_re f and the constraint Formula (22). The fuel cell system power is calculated based on the actual output power of the lithium battery pack P pack and the power demand according to Equation (20) and the constraint Equation (21).

Equalization Energy Management Strategy Based on Switch Network
The switch network [37] shown in Figure 7 can control the number of working batteries by controlling the switches of branch and bypass. For example, when the switch S 11 is closed and S 12 is opened at the same time, the battery B 1 can be connected to the circuit. On the contrary, the battery B 1 can be disconnected. According to the SOC state of each lithium battery and the current power demand, the EEMS performs power distribution, and the switch controller performs real-time control on the switch network.
The flow of the SEMS is: According to the power demand and the SOC status of the current lithium battery pack, the SEMS switches to EEMS mode or ECMS mode to obtain the reference output power of the lithium battery pack _ and the switch network state. The series lithium battery model calculates the reference output current of the single lithium battery  (22). The fuel cell system power is calculated based on the actual output power of the lithium battery pack and the power demand according to Equation (20) and the constraint Equation (21).

Equalization Energy Management Strategy Based on Switch Network
The switch network [37] shown in Figure 7 can control the number of working batteries by controlling the switches of branch and bypass. For example, when the switch 11 is closed and 12 is opened at the same time, the battery 1 can be connected to the circuit. On the contrary, the battery  The main goal of EEMS is to improve the inconsistency of series lithium batteries as soon as possible while meeting the power demand. Therefore, lithium battery can be divided into two categories: high SOC and low SOC according to whether the SOC of lithium battery is higher than the average SOC of the series lithium battery. Count the number of two types of batteries in real time and make decisions according to the number of two types of batteries and the power demand. When the number of high SOC batteries (N high ) is greater than or equal to the number of low SOC batteries (N low ), the power distribution and switching rules are determined as shown in Table 4. Table 4. Rules of equalization energy management strategy (EEMS) (N high ≥ N low ).

State
Power Demand P pack_ref Switch Network State (Disconnect) Table 4, P l is the current load power demand. The fuel cell maximum power P f cmax and minimum power P f cmin are 300 W and 0 W respectively, which is a 25 times reduction in fuel cell simulation results. P b_target is the No.1 battery maximum reference power when the SOC of No.1 battery is equal to the SOC at equilibrium point in ECMS mode. Briefly speaking, in this paper, P b_target is the power when the SOC of No.1 battery is 0.6 and the discharge current is 4 A. P batmax is the charging or discharging power of each battery when the charging or discharging current is I max . N is the number of batteries in the subscript state, h and l are higher than and lower than respectively, max, min and mean are the maximum, minimum and average value of SOC respectively. For example, N h_mean is the number of batteries whose SOC is higher than the mean.
For state 1: Power demand is lower than the minimum power of fuel cell, which indicates that the vehicle is in the state of braking energy recovery. Therefore, the switch network control strategy is: disconnect the batteries whose SOC are equal to max, and charge other batteries.
For state 2: Power demand is higher than the minimum power and lower than the maximum power of the fuel cell. Since the high SOC batteries are the majority, in order to achieve the goal of rapid equalization, the fuel cell needs to charge the low SOC batteries while meeting the power demand. The switch network control strategy is: disconnect the batteries whose SOC are higher than mean, and charge other batteries.
For state 3: Power demand is higher than the maximum power of the fuel cell and lower than the sum of the maximum power of the fuel cell and the reference power of high SOC batteries. In order to meet the power demand, the high SOC batteries outputs the maximum power, and the rest power is supplemented by the fuel cell. The switch network control strategy is: disconnect the batteries whose SOC are lower than mean.
For state 4: The batteries whose SOC are higher than min output the maximum power, and the rest power is supplemented by the fuel cell. The switch network control strategy is: disconnect the batteries whose SOC are equal to min.
For state 5: The power demand exceeds the limit output power of the hybrid system, so all batteries are connected to the circuit for maximum power output.
When the number of high SOC batteries (N high ) is less than the number of low SOC batteries (N low ), only states 2 and states 3 control strategies are different. The power distribution and switching rules are determined as shown in Table 5. Table 5. Rules of EEMS (N high < N low ).

State
Power Demand P pack_ref Switch Network State (Disconnect) 6 P f cmin ≤ P l < N h_mean ·P b_target P l − P f cmin SOC ≥ mean 7 N h_mean ·P b_target ≤ P l < P f cmax + N h_mean ·P b_target P batmax SOC < mean For state 6: Power demand is higher than the minimum power and lower than the reference power of high SOC batteries. Since the low SOC batteries are the majority, in order to achieve the goal of rapid equalization, the high SOC batteries need to supply the power. The switch network control strategy is: disconnect the batteries whose SOC are lower than mean.
For state 7: In order to meet the power demand, the high SOC batteries outputs the maximum power, and the rest power is supplemented by the fuel cell. The switch network control strategy is: disconnect the batteries whose SOC are lower than mean.
In EEMS mode, after obtaining the battery pack reference power P pack_re f by the above rules, the fuel cell reference power is calculated by the following equation: and in any mode, the fuel cell and lithium battery meet the following constraints: The actual power of the fuel cell meets the following constraint: Energies 2020, 13, 247

of 23
The current of batteries meets the following constraint: since the battery life and safety are closely related to the depth of discharge (DOD), the SOC of batteries meets the following constraint: where, SOC l and SOC h are the lower and upper limit of SOC, respectively. If the SOC of a battery is lower than SOC l , it can only be connected to the circuit when it is charged. If the SOC of a battery is higher than SOC h , it can only be connected to the circuit when it is discharged.

Equivalent Consumption Minimization Strategy
When the standard deviation is less than the threshold, it represents the end of the equalization process. The consistency of lithium battery pack is good, which can be regarded as the SOC of each battery is equal. However, EEMS can only make SOC tend to be consistent while the demand power is satisfied and cannot control the SOC value at the end of equilibrium. Therefore, it is necessary to combine other EMS to make SOC tend to target value, and through this kind of combination, fuel economy and system durability can be considered.
ECMS belongs to the real-time optimization strategy. By transforming the power of the battery into the equivalent hydrogen consumption of the fuel cell, the power distribution among the power sources is solved to obtain the best fuel economy. It has been proved to have good fuel economy in Xu, Garcia and Motapon's references [27,46,47].
The fuel economy optimization problem of the fuel cell/battery hybrid system can be described as an optimization problem of controlling the power of the fuel cell and the battery to minimize the total hydrogen consumption of the hybrid system under the constraints of the system, as shown in Equation (25): where, s l and s h are the lower and upper value of the battery pack SOC control range. V pack_min and V pack_max are the lower and upper limit of battery pack voltage, respectively. k is linear coefficient, which can be expressed by: where, µ is SOC balance coefficient. It determines the speed of SOC reaching the equilibrium point. The value ranges from 0 to 1. C f c is the instantaneous hydrogen consumption of the fuel cell, which can usually be expressed as a linear function of the fuel cell output power, as shown in Equation (27): C pack is the equivalent instantaneous hydrogen consumption of the battery pack, which can be expressed as Equation (28): where, P pack is the power of the battery pack. When discharging, P pack is positive; when charging, P pack is negative. η chg_avg and η dis_avg are the mean efficiency of charging and discharging, respectively. P f c_avg and C f c_avg are the mean power and the mean hydrogen consumption. η chg and η dis are the charging and discharging efficiency of the battery pack, which is a function of the current SOC of the battery pack, and can be expressed as Equation (29): where, R dis and R chg are the battery pack discharging and charging resistance, respectively. It can be obtained by the experiment in Section 2. To simplify the presentation, the symbols are defined as follows: substituting Equations (26), (28)- (30) into Equation (24), the optimized power of the lithium battery pack can be obtained as follows: (31) it can be seen from Equation (31) that the optimized power of the battery pack is a function of the current battery pack SOC and the balance coefficient µ. Figure 8 shows the relationship between the optimized power of the battery pack and the current SOC when the balance coefficient µ changes from 0 to 1.
it can be seen from Equation (31) that the optimized power of the battery pack is a function of the current battery pack SOC and the balance coefficient . Figure 8 shows the relationship between the optimized power of the battery pack and the current SOC when the balance coefficient changes from 0 to 1.

Verification and Discussion
To verify that the SEMS can not only meet the power demand for vehicle to work, but also ameliorate the inconsistency of series lithium battery pack, a real-time simulation under a new city European driving cycle (NEDC) condition is conducted. In order to match the power of the 6 series

Verification and Discussion
To verify that the SEMS can not only meet the power demand for vehicle to work, but also ameliorate the inconsistency of series lithium battery pack, a real-time simulation under a new city European driving cycle (NEDC) condition is conducted. In order to match the power of the 6 series batteries, it was also reduced 25 times in the simulation. Figure 9 shows the power demand in the simulation. It can be seen that the entire process lasts about 1370 s, and the peak power is 450 W. When the power value is positive, it denotes that at this moment, the hybrid power system provides the energy to the vehicles. When the power value is below zero, it denotes that the vehicle is recycling the energy during braking. In this simulation, two different initial SOC distribution conditions of battery pack are set as: ℎ ℎ ≥ and ℎ ℎ < . The threshold of the switch is set as 0.001. The sampling time is set as 1 s. The and ℎ in Equation (24) are set as 0.1 and 0.9, respectively. is set as 0.6 in Equation (26). and ℎ are set as 0.4 and 0.8, respectively. Therefore, in ECMs mode, SOC will eventually tend to be around 0.6.
Under the first case ( ℎ ℎ ≥ ), the initial values of SOC from the No.1 battery to the No.6 battery are respectively set as: 0.35, 0.45, 0.6, 0.62, 0.65, 0.68. Figure 10 shows the results of the power demand and the output power of each energy source. It can be seen that the entire process lasts about 1370 s, and the peak power is 450 W. When the power value is positive, it denotes that at this moment, the hybrid power system provides the energy to the vehicles. When the power value is below zero, it denotes that the vehicle is recycling the energy during braking. In this simulation, two different initial SOC distribution conditions of battery pack are set as: N high ≥ N low and N high < N low . The threshold of the switch is set as 0.001. The sampling time is set as 1 s. The SOC l and SOC h in Equation (24) are set as 0.1 and 0.9, respectively. µ is set as 0.6 in Equation (26). s l and s h are set as 0.4 and 0.8, respectively. Therefore, in ECMs mode, SOC will eventually tend to be around 0.6.
Under the first case (N high ≥ N low ), the initial values of SOC from the No.1 battery to the No.6 battery are respectively set as: 0.35, 0.45, 0.6, 0.62, 0.65, 0.68. Figure 10 shows the results of the power demand and the output power of each energy source.
Energies 2020, 13, 247 14 of 23 time is set as 1 s. The and ℎ in Equation (24) are set as 0.1 and 0.9, respectively. is set as 0.6 in Equation (26). and ℎ are set as 0.4 and 0.8, respectively. Therefore, in ECMs mode, SOC will eventually tend to be around 0.6.
Under the first case ( ℎ ℎ ≥ ), the initial values of SOC from the No.1 battery to the No.6 battery are respectively set as: 0.35, 0.45, 0.6, 0.62, 0.65, 0.68. Figure 10 shows the results of the power demand and the output power of each energy source.  In Figure 10, the green dotted line represents the deviation between the power demand and the power that the hybrid system provides. It can be seen during 800 s to 1200 s, there are some periods when power demand exceeds the limit output power of hybrid power system. Thus, both the fuel cell and the batteries output the maximum power. During other periods, the EMS meets the power demands. Figure 11 is the SOC changes of each battery and the standard deviation σ SOC of the battery pack. Figures 12 and 13 are the current of working batteries and the state change of the SEMS. In the first 600 s, σ SOC is higher than the threshold, the series lithium battery pack shows poor consistency. Besides, there are only the No.1 battery and the No.2 battery keeping at the low SOC state. On most occasions, the power demand is lower than the maximum power of the fuel cell, so the control strategy usually works at state 2 where fuel cells supply power to the vehicle and the batteries. At state 2, the switch network control strategy can make the batteries whose SOC is higher than mean disconnected, and with time on, the No.1 battery and the No.2 battery are charged. At 240 s, after SOC of the No.3 battery arrives at the mean, then it fluctuates around the mean with control cycle. And at 280 s, the SOC of the No.2 battery reaches the average value. As the procedure goes on, the SOC of each battery gets close to the mean. At last, when it is 620 s, σ SOC is lower than the threshold, which denotes the SOC of the lithium battery pack is similar and the equilibrium process finished. Then the SEMS turns into the ECMS mode (state 8), where all the batteries are connected into the circuit. Due to the different aging conditions of each battery in the lithium battery pack, their capacity and internal resistance are in different states. Thus, in ECMS mode, even though each battery is charged or discharged with the same current, their SOC changes and output power are still different. In the first case, it can be seen in Figure 11, the SOC of all batteries is 0.66 when the equilibrium process finished. It is very close to the equilibrium point 0.6, the reference output power and the current of the battery pack under ECMS mode are small. Therefore, the SOC of each battery changes slowly, and the state of SEMS does not change very frequently in Figure 13. At 1040 s, the power demand exceeds the limit output power of the hybrid system, and both the fuel cell and the lithium battery are discharged with the maximum power. Since the capacity of the No.6 battery is the smallest, its SOC varies the fastest. Finally, when σ SOC is higher than the threshold, which denotes that the inconsistency of the battery pack is higher than its limit value, then the SEMS turns into the EEMS mode. After 1180 s, the second equilibrium process end, the SEMS turns into ECMS mode.
slowly, and the state of SEMS does not change very frequently in Figure 13. At 1040 s, the power demand exceeds the limit output power of the hybrid system, and both the fuel cell and the lithium battery are discharged with the maximum power. Since the capacity of the No.6 battery is the smallest, its SOC varies the fastest. Finally, when is higher than the threshold, which denotes that the inconsistency of the battery pack is higher than its limit value, then the SEMS turns into the EEMS mode. After 1180 s, the second equilibrium process end, the SEMS turns into ECMS mode.     For the second case, the simulation is also carried out, while the initial value of SOC is set as: 0.45, 0.35, 0.15, 0.1, 0.08, 0.05. Figure 14 shows the results of the power demand and the output power of each energy source. For the second case, the simulation is also carried out, while the initial value of SOC is set as: 0.45, 0.35, 0.15, 0.1, 0.08, 0.05. Figure 14 shows the results of the power demand and the output power of each energy source. It can be seen by comparing the power distribution results with that under the first case, the lithium battery pack is more likely at discharging state. From the deviation curve, although the distribution method is different, the SEMS can satisfy the power demand when the maximum output power of the hybrid power system is above the power demand. Figures 15-17 are also the SOC change of each battery and the standard deviation of the battery pack, the current of working batteries and the state change of the SEMS under the second case. During the first 780 s, is higher than the threshold, the inconsistency of the series lithium battery pack is severe, and No.1 battery and No.2 battery are at high SOC state. The SEMS mostly works at state 6 and state 7. Meanwhile, the switch network control strategy disconnect the batteries whose SOC is lower than the . When the power demand is not high, the power is supplied by the batteries in the high SOC state. When the power demand is high, the high SOC batteries outputs the maximum power, and the fuel cell supply power when necessary. With the No.1 battery and No.2 It can be seen by comparing the power distribution results with that under the first case, the lithium battery pack is more likely at discharging state. From the deviation curve, although the distribution method is different, the SEMS can satisfy the power demand when the maximum output power of the hybrid power system is above the power demand. Figures 15-17 are also the SOC change of each battery and the standard deviation σ SOC of the battery pack, the current of working batteries and the state change of the SEMS under the second case. During the first 780 s, σ SOC is higher than the threshold, the inconsistency of the series lithium battery pack is severe, and No.1 battery and No.2 battery are at high SOC state. The SEMS mostly works at state 6 and state 7. Meanwhile, the switch network control strategy disconnect the batteries whose SOC is lower than the mean. When the power demand is not high, the power is supplied by the batteries in the high SOC state. When the power demand is high, the high SOC batteries outputs the maximum power, and the fuel cell supply power when necessary. With the No.1 battery and No.2 battery discharging, mean continues to decline. When it is 320 s, the SOC of the No.3 battery reaches the mean. And then it fluctuates around the mean with control cycle. As the procedure goes on, all batteries SOC are tend to be consistent at 760 s, the equilibrium process finished. The SEMS turns into the ECMS mode. However, unlike the first case, the SOC of all batteries is 0.14 when the equilibrium process finished in the second case. It is far away from the equilibrium point 0.6. Therefore, it can be seen in Figures 8 and 16, the SEMS charge all batteries with high current in ECMS mode. It causes the SOC of batteries to change rapidly until σ SOC is higher than the threshold, and the SEMS turns into the EEMS mode again. Therefore, the SEMS switches back and forth between EEMS and ECMS mode during 760 s to 1040 s. The power demand exceeds the limit output power of the hybrid system during 1040 s to 1180 s, and all batteries are discharged with the maximum power. After 1180 s, The SEMS returned to charging batteries in ECMS mode and equalizing in EEMS mode at most of the time. And the SOC of batteries will be controlled around 0.6 at last. The enlarged image in the upper right corner in Figure 15 shows the over-discharge protection function in Equation (24) of the SEMS. In Figure 17, the SEMS is in state 4 during 142 s to 144 s. According to the state 4 of in EEMS mode, the SOC of No.5 battery is higher than the min (the SOC of No.6 battery), so it should be discharged with the maximum current. However, its SOC during this period is 0.082, which is lower than the SOC l in Equation (24). In order to prevent the DOD from exceeding the limit, the No.5 battery is disconnected from the circuit and is only connected to the circuit during charging. It can also be seen that the SOC of No.4 battery is higher than the SOC l , so the No.4 battery is discharged, resulting in the SOC of No.4 battery falling to 0.1.     In order to verify the effect of the threshold on the equalization quality, some data are listed in Table 6.
is the maximum value of the from the time when the battery pack reaches equilibrium for the first time to the end of the driving cycle (except the time when the demand power exceeds the limit of hybrid system). The switching times is the sum of the switching times of all battery switches in the entire driving cycle. The equalization time is the time when the battery pack reaches equilibrium for the first time. It can be seen from Table 6 that as the threshold decreases, the battery inconsistency becomes better. This means less safety risk for the battery pack. Alas, reducing the threshold means that more In order to verify the effect of the threshold on the equalization quality, some data are listed in Table 6. σ max is the maximum value of the σ SOC from the time when the battery pack reaches equilibrium for the first time to the end of the driving cycle (except the time when the demand power exceeds the limit of hybrid system). The switching times is the sum of the switching times of all battery switches in the entire driving cycle. The equalization time is the time when the battery pack reaches equilibrium for the first time. It can be seen from Table 6 that as the threshold decreases, the battery inconsistency becomes better. This means less safety risk for the battery pack. Alas, reducing the threshold means that more switching times are needed, and the life of the switch network will be shortened. As the threshold decreases, the equalization time increases, but the increase is not obvious. In actual system design, priority should be given to battery inconsistency and the life of the switch network.

Conclusions
Inconsistency of battery will lead to security hazard and capacity degradation. And the key issue of EMS is to the energy distribution in hybrid systems. However, few EMS take the battery inconsistency into account, which will lead to energy distribution error, battery damage and even security risks. The current battery equalization strategy rarely discusses how to perform the equilibrium process while meeting the power demand of vehicle. To solve these problems, a method is proposed in this paper.
Firstly, a hybrid power system of fuel cell and lithium battery is established based on the simulation of fuel cell and the test results of lithium batteries. Six single battery models with different parameters are used to represent the batteries with different aging degrees in the modeling of lithium battery. Then, a switch control strategy based on the circuit structure of Manenti in Ref. [37] is designed to control the connected state of batteries. After that, an EEMS based on the switch control strategy is designed to improve the inconsistency of the battery pack. ECMS is used to control the range of SOC and ensure the efficiency of the system. The whole EMS (SEMS) switching between EEMS and ECMS in real time based on the comparison of σ SOC with a threshold. In the Matlab/Simulink environment, NEDC working condition was selected for simulation verification. The results show that the SEMS can ameliorate the inconsistency of series lithium battery pack by distribute the power of fuel cell and the series battery pack under the premise of ensuring the vehicle's normal operation. In case 1 and case 2, the inconsistency of lithium battery pack is reduced to below the threshold value for the first time, which costs about 700 s on average. After that time, when the power demand does not exceed the limit of the hybrid system, the inconsistency of the battery pack is effectively controlled. In the whole driving cycle, the number of switching times is about 800 and the selection of threshold affects its equalization quality and switch network life. Therefore, it can be seen that the SEMS can improve the safety and durability of the system, save equalization time. Besides, it has good expansibility and no energy waste. However, there is a disadvantage: with the increase of the number of batteries in series, the number of switches will increase rapidly due to the characteristics of the switch network.
There are some limitations in this article: there is no model for quantitative description of battery life improvement. The battery degradation model and the effect of working conditions on the battery life will be considered for further verification. The future work will optimize the fuel economy of EEMS according to the resistance distribution of lithium battery. In addition, the response speed and durability of fuel cell will be considered in the power distribution method.

Conflicts of Interest:
The authors declare no conflict of interest.