Optimization Analysis of the Energy Management Strategy of the New Energy Hybrid 100% Low-Floor Tramcar Using a Genetic Algorithm

: Performance and economic efﬁciency of the fuel cell (FC)/battery/super capacitor (SC) hybrid 100% low-ﬂoor tramcar is mainly determined by its energy management strategy. In this paper, a train traction model was built to calculate the power output and energy consumption properties of the hybrid tramcar. With the purpose of reducing hydrogen consumption, the genetic algorithm was adopted to optimize the original energy management strategy. The results before and after the optimization show that the power requirement of the tramcar can be satisﬁed in both situations with the fuel cell (FC) module non-stopped. The maximum output power of the FC is reduced from 170 kW to 101.21 kW. As for the SC, a two-parallel connection module is used instead of the three-parallel one, and the power range changes from − 125~250 kW to − 67~153 kW. Under the original energy management strategy, the battery cannot be used efﬁciently with less exporting and absorbent power. Its utilization ratio is improved greatly after optimization. In sum, the equivalent total hydrogen consumption is reduced from 3.3469 kg to 2.8354 kg, dropping by more than 15%.


Introduction
At present, electric power is the main energy source of railway vehicles [1,2]. In China, more than 80% of the electric energy comes from thermal power generation [3]. Meanwhile, the "abundant coal-barren petroleum-lacking gas" energy structure compels the government to spend nearly 50% of coal on electric power generation annually [4,5]. As society develops, the increasing energy demand contributes to the energy crisis. Meanwhile, the large-scale application of primary energy promotes the CO 2 volume fraction in the atmosphere, which is the main cause of global warming. The emissions from the mobile pollution source or industry, such as HC, NO x , and particular matter (PM), greatly threaten human health [6,7]. In this environment, we believe the urban railway system is an effective way to ease the traffic pressure. In addition, the tramcar should have played a more important role in perfecting the city environment by getting rid of the power grid and adopting clean renewable energy as its power [8].
A fuel cell is a kind of power source which transforms chemical energy into electric energy without a combustion process. The proton exchange membrane fuel cell, which is fueled with hydrogen, has the potential for heavy-duty applications due to its advantages of a high conversion efficiency, better low-temperature adaptability, and non-emission [9][10][11]. However, the soft power output characteristic of the fuel cell cannot satisfy the great power demand of the railway vehicle, especially in transient operating conditions. Besides that, the FC cannot recover the braking energy due to the unidirectional energy flow route [12]. Therefore, some other energy modules should be introduced into the power chain, such as the super capacitor, battery, and unidirectional/bidirectional DC/DC. Based on this, a 100% low-floor tramcar is developed using the above introduced hybrid power system. Here, the FC is the primary energy source, while the battery is used for enlarging the energy storage quantity, absorbing the braking energy, and serving as the emergency traction power [13,14]. The super capacitor can also be used for braking energy recovery. Due to its high power density, it is capable of improving the tramcar acceleration response [15,16].
As for the new energy hybrid power tramcar with definite parameters, its performance and economic efficiency is mainly determined by its energy management strategy. Recently, the energy management of the hybrid system has become a research hot spot around the world. Li et al. [17] proposed a power management strategy based on adaptive droop control. By employing the equivalent consumption minimization strategy (ECMS) and multi-mode strategy, the hydrogen consumption of the hybrid tramway is optimized and the safe operating conditions are warranted to fulfill a better energy efficiency of the overall hybrid system and prolong the lifetime of each power source. Han et al. [18] proposed a multi-source coordination energy management strategy based on self-convergence droop control for a large-scale and high-power hybrid tramway. As introduced in their study, the self-convergence droop control approach is capable of limiting the circulating currents, balancing the SOC (state of charged) of each ESS (energy storage system), and modifying the power injections or absorption to fulfill the regulated voltages. Zhang et al. [19] developed a locomotive model using the Advanced Vehicle Simulator (ADVISOR). Based on the model, the power of the fuel cell (FC), battery, and motor are optimized using the bisection algorithm under certain constraints. Wang et al. [20] have presented the DC/DC power converters, including unidirectional boost and bidirectional types, by using time delay control to implement power management for hybrid power systems. Li et al. [21] have presented a power sharing strategy and a state machine strategy based on droop control for a power management system of a hybrid tramway. Fadil et al. [22] have developed a nonlinear controller on the basis of a nonlinear model of a hybrid energy storage system for electric vehicles by using Lyapunov stability design techniques. Wu et al. [23] have presented an adaptive control strategy based on parameter estimations to manage the power sharing in this hybrid power system and satisfy the system constraints. Li et al. [24] have presented a power sharing strategy and a state machine strategy for a proton exchange membrance fuel cell (PEMFC)/battery/super capacitor (SC) hybrid tramway.
These researches promote the performance and energy management strategy of the hybrid system. However, studies on the energy management of the FC/battery/SC hybrid system used for a tramcar, especially the optimization of the whole train simulation based on the railway line, have rarely been reported. In this paper, a whole train simulation model, which contains the energy storage module, power-train module, power bus module, and electric auxiliary system module, was built based on the independently developed FC/battery/SC 100% low-floor tramcar. Then, the genetic algorithm was adopted to optimize the original energy management strategy by selecting the Guanggu T1 uplink in Wuhan, China, for calculation. Therefore, the applicability and quantitative index of the economy of the optimized strategy are available with the help of the calculation model.

Parameters of the Prototype Train
There are two identical power systems in the FC/battery/SC hybrid 100% low-floor tramcar. The main parameters of the tramcar are shown in Table 1. In every power system, 3-in parallel SCs and 2-in parallel batteries are adopted, and the power source configurations are shown in Table 2.

Building of the Numerical Model
The structure of the numerical model of the new energy hybrid 100% low-floor tramcar is shown in Figure 1.

Building of the Numerical Model
The structure of the numerical model of the new energy hybrid 100% low-floor tramcar is shown in Figure 1.  The model starts from the demand of the driving cycle, and then deduces the torque, speed, and power values of every assembly unit. The emulational information flows along the vehicle dynamical model, wheel model, and drive system, and finally reaches the power assembly model. Then, the power demand is transmitted to the drive system according to the torque requirement and limitation. Next, it passes through the wheel module and finally reaches the vehicle module. The calculation error is considered in the reversed calculation process and then fed back to the vehicle module. Some key modules shown in Figure 1 will be introduced in the following parts.

Electric Auxiliary System Module
The electric auxiliary system module includes the air conditioner, illumination, and so on. Power dissipation of the equipment is considered in the module, and it is divided into the constant part and variable part separately, as shown in Equation (1).
where P a,cons is the constant part; and P(s) is a function of the train position, which is defined by the test data. The P(s) can also be replaced by a constant value, P cons , when the test data is unachievable. Then, Equation (1) can be replaced by Equation (2): According to the above analysis, part of the auxiliary system calculation model is shown in Figure 2. The model starts from the demand of the driving cycle, and then deduces the torque, speed, and power values of every assembly unit. The emulational information flows along the vehicle dynamical model, wheel model, and drive system, and finally reaches the power assembly model. Then, the power demand is transmitted to the drive system according to the torque requirement and limitation. Next, it passes through the wheel module and finally reaches the vehicle module. The calculation error is considered in the reversed calculation process and then fed back to the vehicle module. Some key modules shown in Figure 1 will be introduced in the following parts.

Electric Auxiliary System Module
The electric auxiliary system module includes the air conditioner, illumination, and so on. Power dissipation of the equipment is considered in the module, and it is divided into the constant part and variable part separately, as shown in Equation (1).
where , is the constant part; and ( ) is a function of the train position, which is defined by the test data. The ( ) can also be replaced by a constant value, , when the test data is unachievable. Then, Equation (1) can be replaced by Equation (2): According to the above analysis, part of the auxiliary system calculation model is shown in Figure 2.   All the input and output ports are included in the connector module. Here, the input port 1 is the power demand of the transmission system; input port 2 is the actual power output by the power producer which returns from the power bus; and input port 3 is the start switch. Output port 1 is the producer power requirement considering the power consumption of auxiliary systems transmitted to the power bus; and output port 2 is the actual output power returns to the drive system.
In the backward calculation routine, the system total power demand is the sum of the drive system power requirement and the power requirement of auxiliary systems. In the forward calculation routine, the exportable power of the drive system is the difference between actual output power returns from the power bus and the auxiliary system module power dissipation.

Power Bus Module
This sub-system transmits the demand and actual power of every power producer according to the load requirement and the energy management strategy. The structure of the power bus module is shown in Figure 3. All the input and output ports are included in the connector module. Here, the input port 1 is the power demand of the transmission system; input port 2 is the actual power output by the power producer which returns from the power bus; and input port 3 is the start switch. Output port 1 is the producer power requirement considering the power consumption of auxiliary systems transmitted to the power bus; and output port 2 is the actual output power returns to the drive system.
In the backward calculation routine, the system total power demand is the sum of the drive system power requirement and the power requirement of auxiliary systems. In the forward calculation routine, the exportable power of the drive system is the difference between actual output power returns from the power bus and the auxiliary system module power dissipation.

Power Bus Module
This sub-system transmits the demand and actual power of every power producer according to the load requirement and the energy management strategy. The structure of the power bus module is shown in Figure 3. Input port 1 is the power demand from the power producer; input port 2 is the actual power demand from the battery system, which is the difference between the actual output power and the total power demand of the train load; and input port 3 is the actual power output of the super capacitor.

The Train Traction Module
This module is used to calculate the traction of the whole train. The initial velocity used for iterative computations was not given to the module, but the average speed value between the present and last step was used instead. The power demand used for driving the wheel can be calculated by the product of the demand driving force in the current step and the average speed. Meanwhile, in the forward calculation route, the speed for the last step is derived from a speed estimation module.
The modeling method of this model is mainly based on the dynamic equation of the train, using the iterative method to calculate the acceleration of the vehicle and then gradually integrate the speed of the vehicle. The basic and accessional resistances are considered in the calculation. The Input port 1 is the power demand from the power producer; input port 2 is the actual power demand from the battery system, which is the difference between the actual output power and the total power demand of the train load; and input port 3 is the actual power output of the super capacitor.

The Train Traction Module
This module is used to calculate the traction of the whole train. The initial velocity used for iterative computations was not given to the module, but the average speed value between the present and last step was used instead. The power demand used for driving the wheel can be calculated by the product of the demand driving force in the current step and the average speed. Meanwhile, in the forward calculation route, the speed for the last step is derived from a speed estimation module.
The modeling method of this model is mainly based on the dynamic equation of the train, using the iterative method to calculate the acceleration of the vehicle and then gradually integrate the speed of the vehicle. The basic and accessional resistances are considered in the calculation. The structure of the train module is shown in Figure 4, which can be divided into the forward part and backward part. structure of the train module is shown in Figure 4, which can be divided into the forward part and backward part. (1) Backward calculation In the backward calculation routine, the vehicle demand speed tra_spd_r is the average value of the initial and finial speed in the current step. Also, the demanded traction is the sum of the start-basic resistance _ _ _ , the curve resistance tra_curve_resistance, and the slope resistance tra_slope_resistance.
(2) Forward calculation In the forward calculation routine, the sub-module train_speed calculates the initial speed of the next step according to the actual traction and speed from the train wheel module. This is the preparation step for the solution of the demand speed in the calculation path.
The principle of forward calculation is that the quadratic polynomial equation of the average speed in the current step is obtained according to the longitudinal dynamics equation of the train. The average speed of the current step v_aver can be obtained by solving this quadratic or one equation. Then, the final speed v_t is gained by the initial speed v_0 in the current step.
When the actual traction of the train tra_force_a is greater than or equal to zero, the average velocity of the current step v_aver is the smaller value between v_aver and tra_spd_a. Otherwise, the v_aver is taken by the maximum value between v_aver and tra_spd_a.

Energy Storage System Module
The energy storage system module in this study is used to simulate the on-board FC, battery, and SC modules. The module receives power demand instructions from the power bus, returns to the fuel cell, battery, and super capacitor of the actual output power, voltage, current, state of charged, and so on. The positive power value is defined as discharge.
(1) Fuel cell system module (1) Backward calculation In the backward calculation routine, the vehicle demand speed tra_spd_r is the average value of the initial and finial speed in the current step. Also, the demanded traction is the sum of the start-basic resistance tra_start_basic_resistance, the curve resistance tra_curve_resistance, and the slope resistance tra_slope_resistance.
(2) Forward calculation In the forward calculation routine, the sub-module train_speed calculates the initial speed of the next step according to the actual traction and speed from the train wheel module. This is the preparation step for the solution of the demand speed in the calculation path.
The principle of forward calculation is that the quadratic polynomial equation of the average speed in the current step is obtained according to the longitudinal dynamics equation of the train. The average speed of the current step v_aver can be obtained by solving this quadratic or one equation. Then, the final speed v_t is gained by the initial speed v_0 in the current step.
When the actual traction of the train tra_force_a is greater than or equal to zero, the average velocity of the current step v_aver is the smaller value between v_aver and tra_spd_a. Otherwise, the v_aver is taken by the maximum value between v_aver and tra_spd_a.

Energy Storage System Module
The energy storage system module in this study is used to simulate the on-board FC, battery, and SC modules. The module receives power demand instructions from the power bus, returns to the fuel cell, battery, and super capacitor of the actual output power, voltage, current, state of charged, and so on. The positive power value is defined as discharge.
(1) Fuel cell system module There are many kinds of fuel cells. In this study, the Ballard 150 kW PEMFC is adopted as the prototype for modeling. The operating process of the PEMFC is relatively complex. In order to deduce a simplified nonlinear fuel-cell module, some assumptions are made below: (1) Assuming that the working temperature of the FC stack is constant (the response time of the FC to temperature is about 102 s). (2) Assuming that both the anode and cathode of the fuel cell are well humidified, which means the humidity management is considered to be OK. (3) Assuming that the liquid water will not leave the stack, which means the water management is considered to be OK. (4) Assuming that the humidifier and temperature controller work in a reasonable range, which means temperature management is considered to be OK.
To simplify the dynamic model, assuming that the entry reactant mole fraction is constant, namely pure hydrogen (99.95%) is directly inhaled to the anode. Air mixing (according to Nitrogen and Oxygen 79:21) is sent into the cathode, assuming that all gases are ideal.
The equivalent circuit model of the proton exchange membrane fuel cell is shown in Figure 5. There are many kinds of fuel cells. In this study, the Ballard 150 kW PEMFC is adopted as the prototype for modeling. The operating process of the PEMFC is relatively complex. In order to deduce a simplified nonlinear fuel-cell module, some assumptions are made below: (1) Assuming that the working temperature of the FC stack is constant (the response time of the FC to temperature is about 102 s). (2) Assuming that both the anode and cathode of the fuel cell are well humidified, which means the humidity management is considered to be OK. (3) Assuming that the liquid water will not leave the stack, which means the water management is considered to be OK. (4) Assuming that the humidifier and temperature controller work in a reasonable range, which means temperature management is considered to be OK.
To simplify the dynamic model, assuming that the entry reactant mole fraction is constant, namely pure hydrogen (99.95%) is directly inhaled to the anode. Air mixing (according to Nitrogen and Oxygen 79:21) is sent into the cathode, assuming that all gases are ideal.
The equivalent circuit model of the proton exchange membrane fuel cell is shown in Figure 5. According to the steady-state voltage and current characteristics of the fuel cell, the output voltage of the fuel cell is represented by : where is the reversible voltage of fuel cell (open circuit voltage at steady state), V; is the activation voltage drop of the anode and cathode (also known as activation overpressure drop), V; is the voltage drop caused by reactant gas concentration drop or the mass transfer of oxygen and hydrogen (also known as the concentration overpressure drop), V; and is the ohm overvoltage caused by the proton passing through a solid electrolyte or the electron passing through a circuit, V.
The steady-state component of the fuel cell voltage is expressed as: According to the steady-state voltage and current characteristics of the fuel cell, the output voltage of the fuel cell is represented by V fc : where E is the reversible voltage of fuel cell (open circuit voltage at steady state), V; V a is the activation voltage drop of the anode and cathode (also known as activation overpressure drop), V; V c is the voltage drop caused by reactant gas concentration drop or the mass transfer of oxygen and hydrogen (also known as the concentration overpressure drop), V; and V ohm is the ohm overvoltage caused by the proton passing through a solid electrolyte or the electron passing through a circuit, V. The steady-state component of the fuel cell voltage V st is expressed as: The transient component of the fuel cell voltage V tr is: The calculation method of open circuit voltage is as follows: where K c is the voltage constant under rated operating condition; and E n is the Nernst voltage, V. The calculation method of activation voltage drop is as follows: where A is the slope of tafel, V; i 0 is the exchange current for fuel cells corresponding to activation voltage, A; and T d is the fuel cell activation response time, s.
(2) The battery module The equivalent circuit of the battery is shown in Figure 6.  = + The calculation method of open circuit voltage is as follows: where is the voltage constant under rated operating condition; and is the Nernst voltage, V.
The calculation method of activation voltage drop is as follows: where is the slope of tafel, V; is the exchange current for fuel cells corresponding to activation voltage, A; and is the fuel cell activation response time, s.
(2) The battery module The equivalent circuit of the battery is shown in Figure 6.
where Rint is the internal resistance. The basic method for SOC calculation is the Ann-time integral.
(3) SC system module The equivalent circuit of SC is shown in Figure 7.
where R int is the internal resistance. The basic method for SOC calculation is the Ann-time integral.
(3) SC system module The equivalent circuit of SC is shown in Figure 7. The open circuit voltage is calculated by Equation (14): where SOC is the capacity of the SC; and Vmax and Vmin are the maximum and minimum voltage, respectively. Generally, Vmax is the nominal voltage of the SC. If the SC can be emptied, the Vmin is zero. In the full charge state, namely the SOC = 1, the Voc is equal to Vmax. The equivalent series internal resistance of charge and discharge is a function of the internal monomer SC's temperature ess_tmp and current ess_i. The corresponding values of temperature, current, and internal resistance of charge and discharge are obtained by a discrete test, and then the continuous variation of the charge and discharge resistance with temperature and current is gained by interpolation. Compared to the battery model, there is only one set of internal resistance parameters in the SC module. The positive value of the current corresponds to the discharge resistance, while the negative value corresponds to the charge resistance. The internal resistance of the SC group is the result of single SC internal resistance multiplied by the series number and then divided by the number of parallel lines, that is: where r is the single super capacitor internal resistance; ns is the series number; and np is the parallel number.
According to the equivalent circuit: Equation (17) has two solutions. The voltage required to produce the same power is smaller for a large current situation. While the minimum working voltage of the battery is half of the open-circuit voltage Voc. Therefore, the smaller one in the solution is what we need, which means the current can be calculated by Equation (18): The open circuit voltage is calculated by Equation (14):

The Original Energy Management Strategy
where SOC is the capacity of the SC; and V max and V min are the maximum and minimum voltage, respectively. Generally, V max is the nominal voltage of the SC. If the SC can be emptied, the V min is zero. In the full charge state, namely the SOC = 1, the V oc is equal to V max . The equivalent series internal resistance of charge and discharge is a function of the internal monomer SC's temperature ess_tmp and current ess_i. The corresponding values of temperature, current, and internal resistance of charge and discharge are obtained by a discrete test, and then the continuous variation of the charge and discharge resistance with temperature and current is gained by interpolation. Compared to the battery model, there is only one set of internal resistance parameters in the SC module. The positive value of the current corresponds to the discharge resistance, while the negative value corresponds to the charge resistance. The internal resistance of the SC group is the result of single SC internal resistance multiplied by the series number and then divided by the number of parallel lines, that is: where r is the single super capacitor internal resistance; ns is the series number; and np is the parallel number. According to the equivalent circuit: Equation (17) has two solutions. The voltage required to produce the same power is smaller for a large current situation. While the minimum working voltage of the battery is half of the open-circuit voltage V oc . Therefore, the smaller one in the solution is what we need, which means the current can be calculated by Equation (18):

The Original Energy Management Strategy
The original energy management strategy of the FC/battery/SC hybrid 100% low-floor tramcar is a power following principle. The first step of the strategy is to calculate the power restriction of the power source under current conditions. Define the charge/discharge cut-off SOC of the battery and the charge/discharge cut-off voltage of the capacitor. If the mass of the remaining hydrogen is greater than the lowest available hydrogen of the fuel cell, the maximum output power of the fuel cell is defined as the product of the minimum voltage and the maximum current, and the minimum output power of the fuel cell is the product of the maximum voltage and the minimum current. If the mass of the remaining hydrogen is smaller than the lowest available hydrogen of the fuel cell, the maximum and minimum output power is 0. The DC bus side maximum power of the fuel cell is the maximum output power multiplied by DC efficiency. The DC bus minimum power is the minimum output power multiplied by DC efficiency.
If the capacitor voltage is between the charge/discharge cut-off voltages, the maximum output power of the capacitor is the capacitor voltage multiplied by the maximum limiting current of the capacitor. If capacitor SOC is not in this range, the maximum output power of the capacitor is 0. Capacitor minimum output power (i.e., charging power, negative) is the terminal voltage of the capacitor multiplied by the minimum current of the capacitor (negative). The DC bus side maximum power of the capacitor is the maximum output power of the capacitor multiplied by DC efficiency. The DC bus minimum power of the capacitor is the minimum output power of the capacitor multiplied by DC efficiency.
If the battery voltage is between the charge/discharge cut-off voltages, the maximum output power is the battery voltage multiplied by the maximum limiting current. The minimum output power (i.e., charging power, negative) is the current terminal voltage multiplied by the minimum current (negative). If the SOC of the capacitor is not in this range, the maximum and minimum output power of the battery is 0. The DC bus side maximum power of the battery is the maximum output power of capacitor multiplied by DC efficiency. The DC bus minimum power of the battery is the minimum output power multiplied by DC efficiency.
The maximum output power is the sum of the maximum powers of the fuel cell, battery, and SC. Then, the calculation of power distribution is as follows: First, the power on the FC bus side is output based on the minimum allowable power. Then, the power of the FC minus the bus power will be left over. This part of the power will be given priority to the charging capacitor, and then the battery will be charged (considering the maximum charging power of the battery and capacitor).
If the demand power of the bus is equal to or greater than the minimum power of the FC DC side and the demand power of the bus is smaller than the minimum power of the FC DC side, then the output power of the FC is the demand power of the bus with the batteries and SCs do not work. If the demand power of the bus is greater than the maximum power of the FC DC side, fuel cells are powered by the maximum allowable capacity and the remaining part of the demand power is output by the capacitor.
The final actual output power of the FC is the bus demand power of the fuel DC, divided by DC efficiency.

Simulation Results and Analysis
Considering the previously introduced train parameters, energy storage device configuration, and energy allocation strategy, the circuit simulation was carried out based on the T1 line of Wuhan Guanggu. The T1 line has 23 platforms and 22 main parking spaces, as shown in Table 3. All intersections are considered in the calculation to ensure that the supplied energy storage device is equipped with sufficient energy.
Due to the line speed limitation, the maximum train speed was controlled to be no more than 50 km/h. The initial SOC of the battery is 90%, and the initial voltage of the SC is 510 V. In this report, both the maximum discharge and charge currents of the battery are set as 4C, according to the experimental results. These two parameters were not changed in the optimization process to reduce the variables. The maximum charging power is about one third of the maximum discharge power. In order to speed up the optimization of the genetic algorithm, we assume that the maximum battery absorption power is one third of the maximum output power in the energy distribution strategy. _ch_bt = (1/3)·P_disch_bt (19) In the same way, due to the SC maximum output current being 700 A and the maximum absorption current being 300 A, the energy distribution strategy assumes that the maximum absorption power of the SC is three sevenths of the maximum output power.

P_ch_uc = (3/7)·P_disch_uc
(20) Figure 8 shows the train position-speed and time-power diagrams. The dashed line represents the expected value, and the solid line represents the actual running curve of the train. It seems that the train can travel at the desired speed and power with the configuration of the FC/battery/SC hybrid power 100% low-floor tramcar. Table 3. Station and intersections data information of T1. 1  164  64  0  2  300  0  30  3  805  64  30  4  1231  64  30  5  1295  0  30  6  1715  0  30  7  2070  0  30  8  1796  64  30  9  2249  64  30  10  2320  0  30  11  2768  0  30  12  2948  64  30  13  3438  0  30  14  3638  64  30  15  3910  0  30  16  4110  64  30  17  4853  0  30  18  5053  64  30  19  5315  0  30  20  5515  64  30  21  6217  64  30  22  7014  64  30  23  7870  64  30  24  8966  64  30  25  9316  0  30  26  10,100  0  30  27 10,215 64 30 28 10,500 0 30 29 10,820 0 30 30 10,970 64 30   Figure 9 is the output power curves of FC, SC, and power battery versus time. We can see that the FC is not shut down in the whole working condition, and the maximum output power is 170 kW. In the practical engineering use condition, the fuel cell starts slowly with a slow dynamic response. Since the service life of the FC can be reduced when working in a high frequency start-stop condition, the FC does not shut down when the tramcar operates on certain rail tracks. The output power of the SC fluctuates between −125 kW and 250 kW. This is the result of the maximum charge current of the SC is three parts of the maximum discharge current. Under the existing strategy, the battery has a low utilization rate, which means it neither exports nor absorbs much energy most of the time. Figure 10 shows the time-power state curve. The solid line is the battery SOC, and the dotted line is the SC SOC. It seems that the SOC of SC varies between 18% and 100%, while the battery SOC is limited by the energy management system (the battery SOC allows a range of 40~90%), and the whole operating process is basically maintained at about 90%.

Station Number Position (m) Platform Length (m) Stopping Time (s)
Appl. Sci. 2018, 8, x 13 of 22 Figure 9 is the output power curves of FC, SC, and power battery versus time. We can see that the FC is not shut down in the whole working condition, and the maximum output power is 170 kW. In the practical engineering use condition, the fuel cell starts slowly with a slow dynamic response. Since the service life of the FC can be reduced when working in a high frequency start-stop condition, the FC does not shut down when the tramcar operates on certain rail tracks. The output power of the SC fluctuates between −125 kW and 250 kW. This is the result of the maximum charge current of the SC is three parts of the maximum discharge current. Under the existing strategy, the battery has a low utilization rate, which means it neither exports nor absorbs much energy most of the time. Figure 10 shows the time-power state curve. The solid line is the battery SOC, and the dotted line is the SC SOC. It seems that the SOC of SC varies between 18% and 100%, while the battery SOC is limited by the energy management system (the battery SOC allows a range of 40~90%), and the whole operating process is basically maintained at about 90%.  (2) The battery has not been used effectively during the traction stage of the train; (3) The regenerative braking energy loss is great in the braking stage.
If the optimization of the train energy management strategy is carried out, it needs to reduce the energy output of FC and SC, improve the utilization of the power battery in the traction stage of the train, and reduce the loss of regenerative braking energy.   By calculation, the hydrogen consumption of the train is 3.3469 kg, and the running time is 2178 s. The simulation results show that: (1) The energy configuration and energy management strategy of the new energy hybrid tramcar can meet the operational requirements of the experimental line; (2) The battery has not been used effectively during the traction stage of the train; (3) The regenerative braking energy loss is great in the braking stage.
If the optimization of the train energy management strategy is carried out, it needs to reduce the energy output of FC and SC, improve the utilization of the power battery in the traction stage of the train, and reduce the loss of regenerative braking energy. If the optimization of the train energy management strategy is carried out, it needs to reduce the energy output of FC and SC, improve the utilization of the power battery in the traction stage of the train, and reduce the loss of regenerative braking energy.

Selection of Optimization Method
At present, energy management strategies based on optimization methods can be divided into two categories. One type is based on static data tables or historical data, a global optimal energy management method for energy optimization control under specific operating conditions. The other type is based on the real-time status of the vehicle or the current parameters of the recurrence control, which can usually guarantee a local or instantaneous optimal energy allocation strategy [25,26].
Since the instantaneous optimal energy allocation strategy is generally aimed at automotive operating conditions, and depending on the operating characteristics of the tram, its operating conditions can be determined. Therefore, the global optimal energy management method is selected here. Based on the energy management strategy of the global optimal, the most representative are the dynamic programming control method [27], Pontryagin minimum method [28], genetic algorithm [29], and other intelligent control methods. Among them, the genetic algorithm pays no attention to the system model. It just needs to determine how the evaluation function can be evolved through optimization. It is also advantageous in providing more convenience in changing the energy allocation strategy in the original tram and it is easy to implement in the project. So, the genetic algorithm was adopted in this report to optimize the energy allocation strategy in the original tramcar [30].

Target Fitness Function
The optimization goal of this study is to reduce the fuel consumption of the fuel cell without reducing the performance of the tram. Therefore, the target fitness value function used to evaluate the tram energy allocation strategy is: Among them, X is the set of variables related to calculating the fitness value including the tram, and f (X) is the fitness value function. S is the search space where X is located. m H 2 (X) is the mass of hydrogen consumed by trams after a complete operating condition. m chg,bat (X) and m chg,sc (X) are the electrical energy consumed by the grid to charge the battery and capacitor to their initial state after a complete operating condition, respectively, and the electrical energy is converted into fuel quality: m chg,bat (X) = E chg,bat (SOC init,bat − SOC end,bat )·η f uel,ele ζ dc (21) m chg,sc (X) = E chg,sc (SOC init,sc − SOC end,sc )·η f uel,ele Among them: ζ dc is the efficiency of the energy passing through the DC/DC converter when the grid is charging the battery and the super capacitor; E chg,bat () and E chg,sc () are the required energy calculation functions (kWh), and are related to the state of the charge of the battery/super capacitor at the beginning and the end of the journey, respectively; SOC init,bat and SOC end,bat are the battery states at the beginning and the end of the journey, respectively; SOC init,sc and SOC end,sc are the SOC of the super capacitor at the beginning and the end of the journey, respectively; η f uel,ele is the conversion coefficient of electric energy into fuel quality, which is based on previous research, and according to a more conservative estimation, the lower price of hydrogen fuel is ¥20 per kilogram, and the higher price of industrial electricity is ¥0.9 per kilowatt-hour, so it is converted here. The coefficient is taken as 0.9/20; ζ dc is the efficiency of the DC/DC that the energy passes when charging the battery/super capacitor by grid.
X is subject to the following constraints: Among them, SOC bt_min and SOC bt_max are the minimum and maximum values of the allowable range of the battery state, which are taken as 40% and 90% in this study, respectively; and SOC uc_min and SOC uc_min are the minimum and maximum values of the allowable capacitor state, which are taken as 18% and 96% in this study, respectively (100% in simulation).

Process of Energy Allocation Strategy Optimization
According to the description of Section 2.1, the structure of the energy allocation strategy for tram is shown in Figure 11.
From the foregoing, we can see that the optimization object of the genetic algorithm is the power target in the energy allocation strategy: P_disch_fc, P_disch_bt, and P_disch_uc, which represent the maximum output power of the fuel cell, the maximum output power of the power battery, and the maximum output power of the super capacitor. After determining the optimization objectives and the optimization target fitness value function, the genetic algorithm can be optimized for the energy distribution strategy of the tram. The energy allocation strategy genetic algorithm optimization process is shown in Figure 12. is subject to the following constraints: Among them, _ and _ are the minimum and maximum values of the allowable range of the battery state, which are taken as 40% and 90% in this study, respectively; and _ and _ are the minimum and maximum values of the allowable capacitor state, which are taken as 18% and 96% in this study, respectively (100% in simulation).

Process of Energy Allocation Strategy Optimization
According to the description of Section 2.1, the structure of the energy allocation strategy for tram is shown in Figure 11.
From the foregoing, we can see that the optimization object of the genetic algorithm is the power target in the energy allocation strategy: P_disch_fc, P_disch_bt, and P_disch_uc, which represent the maximum output power of the fuel cell, the maximum output power of the power battery, and the maximum output power of the super capacitor. After determining the optimization objectives and the optimization target fitness value function, the genetic algorithm can be optimized for the energy distribution strategy of the tram. The energy allocation strategy genetic algorithm optimization process is shown in Figure 12.

Simulation Calculation and Analysis of Optimization Strategy
According to the simulation method in Section 3.2, the genetic algorithm optimization simulation of the energy distribution strategy of trams is carried out. In the simulation, the number of individuals for each generation of population is 50, and a total of 100 generations have been evolved. The simulation results are as follows: P_disch_bt = 91.57 kW; P_disch_uc = 153.74 kW; P_disch_fc = 101.21 kW Figure 13 shows the evolutionary process diagram of the genetic algorithm simulation. It can be seen that at the 36th generation, the fitness level reached the maximum, which corresponds to the minimum fuel consumption of the tram. From the 36th generation to the 100th generation, the optimal solution of the genetic algorithm remains unchanged, that is, the algorithm converges. In the end, the entire journey fuel consumption is equivalent to 2.8354 kg.

Simulation Calculation and Analysis of Optimization Strategy
According to the simulation method in Section 3.2, the genetic algorithm optimization simulation of the energy distribution strategy of trams is carried out. In the simulation, the number of individuals for each generation of population is 50, and a total of 100 generations have been evolved. The simulation results are as follows: P_disch_bt = 91.57 kW; P_disch_uc = 153.74 kW; P_disch_fc = 101.21 kW Figure 13 shows the evolutionary process diagram of the genetic algorithm simulation. It can be seen that at the 36th generation, the fitness level reached the maximum, which corresponds to the minimum fuel consumption of the tram. From the 36th generation to the 100th generation, the optimal solution of the genetic algorithm remains unchanged, that is, the algorithm converges. In the end, the entire journey fuel consumption is equivalent to 2.8354 kg.      Figure 15 is the power output characteristics of each energy storage system after optimization. We can see that the FC is not shut down in the whole working condition, and the maximum output power is 101.21 kW. The output power of the SC fluctuates between −125 kw and 250 kW. This is because the maximum charge current of the SC is three sevenths of the maximum discharge current. As for the battery, the utilization ratio is improved efficiently by exporting and absorbing more energy, which is the main reason for using the optimized energy management strategy to reduce the fuel consumption.     Figure 15 is the power output characteristics of each energy storage system after optimization. We can see that the FC is not shut down in the whole working condition, and the maximum output power is 101.21 kW. The output power of the SC fluctuates between −125 kw and 250 kW. This is because the maximum charge current of the SC is three sevenths of the maximum discharge current. As for the battery, the utilization ratio is improved efficiently by exporting and absorbing more energy, which is the main reason for using the optimized energy management strategy to reduce the fuel consumption.  Figure 15 is the power output characteristics of each energy storage system after optimization. We can see that the FC is not shut down in the whole working condition, and the maximum output power is 101.21 kW. The output power of the SC fluctuates between −125 kw and 250 kW. This is because the maximum charge current of the SC is three sevenths of the maximum discharge current. As for the battery, the utilization ratio is improved efficiently by exporting and absorbing more energy, which is the main reason for using the optimized energy management strategy to reduce the fuel consumption.  Figure 16 shows the tram running time-power status curve. It can be seen that the SOC of the super capacitor changes from 18% to 100%, and the SOC of the power battery drops significantly. This further proves that the optimized power distribution strategy has improved the utilization rate of power batteries. It should be noted that this figure only shows the status information of the battery and the super capacitor during the running of the tram, and does not record the process of charging the battery and capacitor to the initial state after the tram arrives at the station. This charging process has been considered in the optimization process.  Figure 16 shows the tram running time-power status curve. It can be seen that the SOC of the super capacitor changes from 18% to 100%, and the SOC of the power battery drops significantly. This further proves that the optimized power distribution strategy has improved the utilization rate of power batteries. It should be noted that this figure only shows the status information of the battery and the super capacitor during the running of the tram, and does not record the process of charging the battery and capacitor to the initial state after the tram arrives at the station. This charging process has been considered in the optimization process.

Conclusions
(1) A simulation model was built based on the FC/battery/SC hybrid 100% low-floor tramcar, which can be used to predict the performance and fuel consumption of the new energy hybrid system. (2) Under the original energy management strategy, the fuel consumption of the energy storage system is relatively high with a low capacity utilization ratio. The FC module never stops during the whole process, with a maximum output power of 170 kW. The power range of the super capacitor is −125 kW to 250 kW, which is the result of the charge current being three-sevenths of the discharge current. The battery was not used efficiently for its tiny exporting and absorbent energy. The equivalent energy consumption is 3.3469 kg for the whole line simulation. (3) The genetic algorithm can be used to realize the optimization solution of the energy management strategy. The optimal objects are set as the maximum output power of the FC P_disch_fc, battery P_disch_bt, and SC P_disch_uc. The maximum fitness value corresponding to the minimum fuel consumption and the algorithm converged when the optimal solution remained unchanged. (4) The optimized energy management strategy can improve the economic efficiency of the new energy hybrid tramcar. The maximum output power of the FC dropped to 101.2 kW and never stopped during the whole running process. The SC module is reduced from 3-in parallel to 2-in parallel, and its output power drops to −67 kW to 153 kW. As for the battery, the use ratio is improved greatly by exporting and absorbing more energy. The equivalent energy

Conclusions
(1) A simulation model was built based on the FC/battery/SC hybrid 100% low-floor tramcar, which can be used to predict the performance and fuel consumption of the new energy hybrid system. (2) Under the original energy management strategy, the fuel consumption of the energy storage system is relatively high with a low capacity utilization ratio. The FC module never stops during the whole process, with a maximum output power of 170 kW. The power range of the super capacitor is −125 kW to 250 kW, which is the result of the charge current being three-sevenths of the discharge current. The battery was not used efficiently for its tiny exporting and absorbent energy. The equivalent energy consumption is 3.3469 kg for the whole line simulation. (3) The genetic algorithm can be used to realize the optimization solution of the energy management strategy. The optimal objects are set as the maximum output power of the FC P_disch_fc, battery P_disch_bt, and SC P_disch_uc. The maximum fitness value corresponding to the minimum fuel consumption and the algorithm converged when the optimal solution remained unchanged. (4) The optimized energy management strategy can improve the economic efficiency of the new energy hybrid tramcar. The maximum output power of the FC dropped to 101.2 kW and never stopped during the whole running process. The SC module is reduced from 3-in parallel to 2-in parallel, and its output power drops to −67 kW to 153 kW. As for the battery, the use ratio is improved greatly by exporting and absorbing more energy. The equivalent energy consumption is 2.8354 kg for the whole line simulation, which is a reduction of 15% compared to the original strategy. Funding: This research received no external funding.

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