A Novel Energy Management Strategy for Series Hybrid Electric Rescue Vehicle

The performance and fuel consumption of hybrid electric vehicle heavily depend on the EMS (energy management strategy).This paper presents a novel EMS for a series hybrid electric rescue vehicle. Firstly, considering the working characteristics of engine and battery, the EMS combining logic threshold and fuzzy control is proposed. Secondly, a fuzzy control optimization method based on IQGA (improved quantum genetic algorithm) is designed to achieve better fuel efficiency.Then, the modeling and simulation are completed by using MATLAB/Simulink; the results demonstrate that the fuel consumption can be decreased by 5.17% after IQGA optimization and that the optimization effect of IQGA is better than that of GA (genetic algorithm) and QGA (quantum genetic algorithm). Finally, the HILS (hardware in loop simulation) platform is constructedwith dSPACE; the HILS experiment shows that the proposed EMS can effectively improve the vehicle working efficiency, which can be applied to practical application.


Introduction
In recent years, natural disasters (such as earthquakes, debris flow, and heavy snowfall) have often occurred, causing inestimable loss to the people's life and property.As an important part of the emergency response system, rescue vehicles have a direct impact on the efficiency of disaster relief operations.The application of hybrid electric technology to rescue vehicles has aroused people's attention, by which the vehicle mobility and fuel efficiency can be significantly improved.At the same time, the generator can be used to provide power for lighting equipment, rescue equipment, and other electrical equipment, so the vehicle adaptability in harsh environment can be enhanced.However, the EMS plays a key role in the performance of hybrid electric system [1], which can be roughly classified into the following three categories.
(1) Logic Threshold Control Strategy.This strategy, which is widely used in many automobile enterprises, is simple and easy to develop, mainly including the "thermostat" control strategy and power following control strategy [2].By setting the logic threshold to establish control rules, each power component can work in the high efficient area as much as possible.However, the setting of threshold value depends on the experiences of experts, and the optimal control effect cannot be guaranteed [3,4].
(2) Optimization Based Control Strategy.By defining objective function, combining the constraint conditions, the control objective can be optimized by mathematical algorithm.This strategy can be subdivided into global optimization control strategy and instantaneous optimization control strategy.The former is designed based on given driving cycle, which is often used for vehicles with fixed routes, such as buses or commuter cars.Besides, it is also used as a criterion to evaluate other control strategies.Because the calculating quantity is large and complex, the optimization process cannot be carried out online.The representative strategies mainly include dynamic programming (DP) based strategy [5,6], Pontryagin's minimum principle (PMP) based strategy [7,8], Quadratic programming (QP) based strategy [9,10], divide rectangle (DIRECT) algorithm based strategy [11][12][13], and convex optimization based strategy [14].The latter is aimed at achieving optimal energy management in instantaneous state, which is not restricted by the specific driving cycle.The computation is relatively small and can be applied online, such as fuel consumption minimization strategy (ECMS) [15][16][17] and model predictive control strategy (MPC) [18][19][20].However, the instantaneous optimization belongs to local optimization, and the optimization effect depends on the precision of the mathematical model.
(3) Fuzzy Logic Control Strategy.With good real-time performance and robustness [21,22], the fuzzy logic control strategy does not rely on precise mathematical models, which is suitable for dealing with the complex, nonlinear and time-varying problems.In practical applications, the speed, SOC, required torque, and power are often used as input variables [23][24][25], and the power ratio of the engine and the battery is used as the output variable.However, the formulation of control logic is mainly based on engineering experience, and it is difficult to ensure the optimal power distribution.To achieve better performance, many strategies combining fuzzy control method with other control methods have been proposed by some scholars.In [26], a fuzzy controller was used to adjust the equivalent factors of ECMS, and the fuel economy of the hybrid electric heavy-duty vehicle was improved.In [27], the DP algorithm was used to optimize the performance of ISG hybrid electric vehicle, and combined with the optimization results, a more practical fuzzy controller was proposed.In [28][29][30], some intelligent algorithms, such as genetic algorithm and particle swarm optimization, were used to optimize the membership function or control rules of the fuzzy controller.In [31][32][33], the mixed-fuzzy control strategy was designed by using neural network, working condition identification, machine learning, and other intelligent techniques.However, the control effects of the above strategies still exist a certain improvement space.
Considering the advantages and disadvantages of the above three categories of methods, a novel EMS is proposed to improve the fuel economy of a series hybrid electric rescue vehicle in this paper.First, combined with threshold control and fuzzy control, the EMS is initially designed according the working characteristics of engine and battery.Then, an improved quantum genetic algorithm with better convergence performance is proposed and applied to optimal the fuzzy controller.The structure of the paper is as follows: in Section 2, the configuration of a series hybrid electric rescue vehicle is introduced; in Section 3, the EMS of the whole vehicle is studied; in Section 4, the method for optimizing fuzzy control based on IQGA is designed; in Section 5, the offline simulation and HILS experiment are described; finally, the paper is concluded in Section 6.

Configuration of the Series Hybrid Electric Rescue Vehicle
In series hybrid electric system, there is no mechanical connection between the engine and the driving wheels, the speed, and torque of the engine can be controlled to any operating point on its speed-torque map.So by optimizing control, the fuel consumption and emission of the engine can be improved.As shown in Figure 1, the four-wheel drive structure with two motors is adopted to the rescue vehicle in this paper, which can not only increase the vehicle power performance, but also provide better regenerative braking effect.The electric drive system mainly includes general controller, internal-combustion-engine-generator power unit (IGPU) and its controller, battery pack and its management system (BMS), motors and motor controllers, rectifier, and DC/DC converter.The general controller is the core component, which is used to coordinate the power distribution of IGPU and battery pack according to the driver instruction, vehicle running state and other feedback information.As an auxiliary energy source, the battery pack provides additional power when the IGPU output power is insufficient, such as the situation of vehicle uphill and acceleration.The two motors can be used not only to provide driving force, but also to work as a generator to charge the battery when the vehicle is braking.The specific parameters of the vehicle are shown in Table 1.

Energy Management Strategy Design
Under the condition of satisfying vehicle power performance, how to coordinate the energy distribution between IGPU and battery to realize the optimization of vehicle efficiency is a primary problem that needed to be resolved for EMS.
The following principles are taken into consideration in the formulation of EMS: (1) In order to improve engine efficiency, the minimum fuel consumption curve control method is adopted.(2) In order to make the battery not only able to provide power in time, but also able to recover the regenerative braking energy effectively, the battery SOC is controlled in the high efficiency area, and a certain margin is reserved for charging.(3) On the premise of braking stability, in order to make maximum use of the two motors to recover braking energy, the braking force distribution method based on I curve is adopted.
Through the experiments, the IGPU fuel consumption characteristics and battery internal resistance characteristics are obtained as shown in Figures 2 and 3, respectively.It is known that, when the IGPU output power is less than 50KW, the fuel economy is bad, and when SOC is in the range of 0.5∼0.9, the battery working efficiency is high.Therefore, the IGPU output power is controlled within the range of 50∼ 260KW, and the battery SOC is controlled within the range of 0.5∼0.8(SOC low =0.5, SOC high =0.8).According to the driving mode switch, battery SOC and required power, the vehicle working state is divided into four operating modes: pure electric traction mode, IGPU traction mode, hybrid traction mode, and braking mode.Each mode switches by the set threshold, and the flowchart of the EMS is shown in Figure 4.   is the required power of two motors,    is the target power of IGPU,    is the target power of battery, and    is the target power of brake resistance.Define  max is the maximum charging power of battery, and  max is the maximum discharge power of battery.Assume that the charging current and power of the battery are negative and that the discharge current and power of the battery are positive; the energy distribution of different modes is as follows.
3.1.Pure Electric Traction Mode.In this mode, the battery alone supplies its power to meet the power demand.And this mode includes two cases (module 1 and module 2 in Figure 4): when the pure electric switch is opened, the vehicle is in pure electric drive mode A; when 0 ≤   ≤ 50KW and SOC > SOC low , the vehicle is in pure electric drive mode B. Both of these two cases are satisfied with the following formula:

IGPU Traction Mode.
In this mode, IGPU not only provides power to propel the vehicle, but also provides power to charge the battery.There are also two cases included in this mode (module 3 and module 4 in Figure 4).When 0 ≤   ≤ 50KW and SOC ≤ SOC low , the vehicle is in IGPU traction mode A, and the battery is charged with the maximum charging current: When   >50KW and SOC≤ SOC low , the vehicle is in IGPU traction mode B:

Hybrid Traction
Mode.This mode is also divided into two cases (module 5 and module 6 in Figure 4).When   >50KW and SOC≥ SOC high , the vehicle is in hybrid traction mode A, at this time, the required power is mainly supplied by the battery, and the additional power is supplemented by IGPU; that is, When   >50KW and SOC low < SOC < SOC high , the vehicle is in hybrid traction mode B. Considering the efficiency of battery and the fuel economy of engine, the power distribution between the battery and IGPU is carried out by a Mamdani fuzzy logic controller (FLC) with double inputs and single output.The IGPU load coefficient  ( =   /260) and battery SOC are selected as input variables, power distribution coefficient   is selected as output variable.The target power of the battery and IGPU can be expressed as The relative membership function diagram is shown in Figure 5. , SOC, and   all have five fuzzy subsets: VS, MS, M, MB, VB, and the basic domain of , SOC, and   , respectively, are [0.1, 1.2], [0.5, 0.8], and [0.7, 1.5].In the design process of FLC, the following principles should be noted: When  is small,   should be raised so as to make the engine operate in its optimal operation region.When  is large, k g should be reduced so that the battery can supplement the additional power.When SOC is large or small,   should be reduced or raised so that the battery can be discharged or charged in time to keep SOC in the area with high efficiency.By using the rule form of "if  is A, and SOC is B, then   is C", 25 fuzzy rules are established as shown in Table 2.

Braking Mode.
Under driving mode,   can be calculated according to the angular displacement of accelerator pedal.However, for braking mode,   cannot be calculated by using the angular displacement of brake pedal, because the electrical regenerative braking coexists with mechanical braking.In this paper, the braking force of the motors is  calculated by using of I curve, and then   is obtained combining with the vehicle speed.I curve is an ideal braking force distribution curve, as shown in Figure 6, where j is the deceleration of the vehicle during braking.Under this condition, the front axle braking force  1 and rear axle braking force  2 are satisfied with the following formula: where G is the vehicle gravity, ℎ  is the height of vehicle centroid, b is the distance between vehicle centroid and rear axle, and L is the wheelbase.
Besides, the sum of  1 and  2 is equal to the overall braking force; that is, where  is conversion coefficient of brake pedal angle displacement ( ∈ [−1, 0]) and   max is the maximum I curve j = 1 .0 g j = 0 .9 g j = 0 .8 g j = 0 .7 g j = 0 .6 g j = 0 .5 g j = 0 .4 g j = 0 .3 g j = 0 . 2 g j = 0 . 1 g Rear braking force ratio Front braking force ratio (F 1 /G) braking force of the vehicle.Combining formula (6) and formula (7),  1 and  2 can be obtained. 1 and  2 are usually composed of mechanical braking force and regenerative braking force.For more energy recapture, the regenerative braking should be given priority.If  1 ( 2 ) is smaller than the maximum regenerative braking force of motor, the total braking force is provided by the motor and the mechanical brake is not used, as shown by point a in Figure 6.If  1 ( 2 ) is greater than the maximum regenerative braking force of motor, the motor is controlled to produce the maximum regenerative braking force and the remaining force  1 ( 2 ) is provided by mechanical brake as shown by point b in Figure 6; that is, where  1 ,  2 are the target regenerative braking force on the front axle and rear axle, respectively. max is the maximum regenerative force of motor, and its value is limited by the motor speed N. In addition, it should be noted that the regenerative braking is affected by braking strength Z, vehicle speed V, and battery SOC.
When 0<Z<0.7, the motor can be controlled for regenerative braking.When Z≥0.7, the vehicle is in an emergency braking state.Considering braking safety, the motor should not work, and the mechanical brake has to produce the total braking force as required.The braking strength correction factor  1 is introduced as When V≤12km/h, it is hard for the motors to produce brake torque due to the low speed.Therefore, with the vehicle speed drops to 12 km/h, the regenerative braking force should Mathematical Problems in Engineering be gently reduced to 0. The speed correction factor  2 is introduced as When the value of SOC is very high (SOC>0.99),the battery will be damaged if it is still charged.Therefore, with the battery SOC up to 0.99, the regenerative braking force should also be gently reduced to 0. The SOC correction factor  3 is introduced as The actual regenerative braking force on the front axle and rear axle   1 and   2 can be expressed as Therefore, the   under braking state can be expressed as where r is wheel radius and  0 is final drive ratio.The braking state can be divided into two cases according to SOC value (module 7 and module 8 in Figure 4).When SOC< SOC low , the vehicle is in braking mode A. The battery should be charged with maximum current, and the brake resistance does not work; that is, where   is the maximum regenerative braking coefficient.When SOC≥ SOC low , the vehicle is in braking mode B. The battery is charged normally and IPGU does not work.If the recovery energy is too much, the excess energy can be consumed in the form of thermal energy by using the brake resistance; that is,

Implementation of the IQGA for FLC
Because the FLC design is mainly based on the engineering experiences of experts, there is room for improvement.In this section, in order to improve the fuel economy of hybrid power system, the IQGA is proposed to optimize the fuzzy control strategy.

4.1.
Improved Quantum Genetic Algorithm.The quantum genetic algorithm (QGA) is a probability search algorithm based on the theory of quantum computing, which is proposed by Han [34].In QGA, quantum state vector is used for genetic encoding, and qubit rotation gate is used to realize chromosome evolution.Therefore, the algorithm has strong parallelism, which can achieve better calculation results than conventional genetic algorithm (GA).However, its convergence performance leaves room to improve.The IQGA is proposed in this paper based on the improvement of QGA.The dynamic adjusting rotation gate is adopted to improve the convergence speed, and quantum catastrophic operation is introduced to overcome premature convergence and improve global optimization ability.

Qubit Encode.
A qubit can be interpreted not only "0" state or "1" state, but also any superposition of them.It can be expressed as where  and  are complex numbers, satisfying || 2 + || 2 = 1.|| 2 and || 2 represent the probability of qubit in 0 state and 1 state, respectively.The multiqubits can be used to represent the chromosome with multigenes, as follows: where    is the j-th individual chromosome in t-th generation, k is the number of qubits used for each gene encoding, and m is the number of genes in the chromosome.
4.1.2.Measure.The measuring method for the t-th generation population () = (  1 ,   2 , . . .,    ) is as follows: For each qubit, generate a random number N rand between 0 and 1; if N rand ≤ || 2 , the qubit is 0; otherwise 1.Then the binary Then the decimal measurement value of the th chromosome    = ( 1 ,  2 , ⋅ ⋅ ⋅,   ) can be calculated.By the same method, the decimal value of the t-th generation population The updating operation of qubit is as follows: where (  ,   )  and (   ,    )  are the probability amplitude of th qubit before and after updating, respectively. is the rotation angle, in general,  ∈ [0.001, 0.05], and it can be expressed as where Δ  and (  ,   ) are the value and direction of the rotation angle, and they can be obtained by Table 3.In Table 3,   and   , respectively, are the th bit of the current measuring individual and optimal individual; (  ) and (), respectively, are the fitness values of the current measuring individual and optimal individual. is a dynamic value, which is used to control the convergence speed of the algorithm; it can be expressed as where t is the evolution generation,  is the largest evolution generation, and  max and  min are the maximum rotation angle and minimum rotation angle, respectively. is decreased with the increase of evolution generation, so as to realize the dynamic adaptive adjustment of the qubit rotation gate.

Population Catastrophe.
After several generations of evolution, the algorithm may fall into the local optimal solution.The population catastrophe can be used to solve this problem.For each evolution, the average fitness of the population  V can be calculated as follows: where   represents the fitness value of the th individual and n is the number of individuals in the population.The fitness values which is greater than  V are averaged to get   V and define the catastrophe factor Δ as If Δ is very small, it indicates that the diversity of the population is poor, and the catastrophe should be carried out at this time; that is, only the best individual is retained and the others should be regenerated.

FLC Optimization.
The membership function is the key factor affecting the performance of FLC, so the decisive parameters of the membership function are optimized in this paper.Take  as an example; as shown in Figure 7, the location of the vertexes of the red triangles A, B, and C determines the shape of the membership function.Therefore, the distance between the line segments AD, AE, and EC (i.e., x1, x2, and x3) can be used as the variables to be optimized.The number of variables can be reduced effectively by this method, which greatly reduces the computational complexity.The same coding method is used to optimize the membership function of SOC and   .Thus there are 9 variables to be optimized; that is, there are 9 genes in the chromosome.
During the driving cycle, both the battery and the enginegenerator set need to provide energy for the vehicle.To eliminate the effect of battery energy on fuel consumption,  the battery SOC of final state before and after optimization should be controlled at the same level.The ideal situation is that the final state SOC difference ΔSOC is controlled to 0, but in this case, the amount of calculation will become very large.In order to improve the speed of simulation, a certain tolerance of ΔSOC can be set.In this paper, the tolerance is 0.5% [35,36]; that is, the constraint on SOC is defined as Under the condition of satisfying vehicle power performance, aiming at the goal of low fuel consumption in urban and high speed driving cycle, the objective function is obtained as follows [37]:  = 1 0.55/   + 0.45/   (26) where    and    are the engine fuel consumption (L) under the urban and high speed driving cycle, respectively.
The optimization process of the fuzzy controller based on IQGA is shown in Figure 8, and the specific steps are as follows: (1) Initialize the population ( 0 ); randomly generate n chromosomes encoded with qubits.
(3) Assign ( 0 ) to FLC, evaluate ( 0 ), and record the optimal individual   0  and its fitness   0  .(4) Decide whether the calculation process can be finished; if it meets the termination condition, end, or go on.

Simulation and HILS Experiment
In this section, the effectiveness of the proposed EMS for series hybrid electric rescue vehicle is validated.Firstly, the simulation model of the vehicle is built with Simulink.Then, the offline optimization simulation is carried out.Finally,  (1) Vehicle Dynamics Model.The vehicle dynamics model is built based on longitudinal dynamics characteristics of the vehicle; the torque of wheel   can be expressed as where  1 and  2 are the transmission efficiency of the front and rear motor, respectively. 1 and  2 are the output torque of the front and rear motor, respectively;  ℎ1 and  ℎ2 are the mechanical braking torque of the front and rear axle, respectively.  can also be expressed as where m is the vehicle mass (kg), g is the gravity acceleration, C r is the rolling resistance coefficient,  is the road-grade angle,   is the air drag coefficient,  is the air density, A is the frontal area of the vehicle (m 2 ),   is the correction coefficient of rotating mass, and u is the vehicle speed (m/s).The front and rear motor speed N 1 and N 2 (r/min) have the following relationship with u: (2) Engine-Generator Model.Because there is velocity coupling between the engine and generator, they are modeled as a whole.The equivalent circuit of engine-generator set is shown in Figure 10 [39].Ignoring the internal resistance and torque loss, the DC bus voltage U dc and the generator electromagnetic torque T g can be expressed as where K e and K x are the electromotive force coefficient and equivalent impedance coefficient of the generator, respectively, and they can be obtained from the bench test data. g is the rotation angular velocity of generator; I dc is the DC bus current.The engine and the generator are connected by a speed increase box, and the speed increasing ratio is i eg .According to the torque balance, the following equations can be obtained: where T eng is the engine torque and J e and J g are the moment of inertia of engine and generator, respectively.n eng and n g are the rotational speed of engine and generator, Mathematical Problems in Engineering respectively.In addition, the engine fuel consumption V e can be calculated from formula (32): where b e is the fuel consumption rate (b e =f (T eng , n eng )),   is the fuel density, and   is the engine output efficiency.
(3) Battery Model.The battery model is built by using the typical Rint model [40], and the equivalent circuit is shown in Figure 11.U b , V 0 , I b , R ch , and R dch are the terminal voltage, open-circuit voltage, internal current (positive when discharge), charge resistance, and discharge resistance of the battery, respectively.Among them, the values of V 0 , R ch , and R dch are all related to the battery SOC, so the following equations can be obtained: The ampere hour method [41] is applied to estimate the battery SOC.Suppose the maximum capacity of the battery is  max , and the initial capacity of the battery is Q 0 , then the battery SOC at time t can be expressed as 5.2.Simulation.The driving cycle used in the simulation, as shown in Figure 12, is mainly composed of UDDS urban driving cycle and HWFET high speed driving cycle.In order to compare the optimization effect of different algorithms, GA, QGA, and IQGA are used to optimize the FLC, respectively.Supposing that the genetic algebra is 200, that the population size is 50, and that the binary length of each variable is 10, the simulation is carried out by using MATLAB/Simulink.The  changing curves of global fitness with different strategies are shown in Figure 13.It can be seen that, under the condition of same population size and generation, the fitness of IQGA reaches the minimum value 3.8 after 80 generations, while the fitness of QGA reaches the minimum value 3.85 after 114 generations and the fitness of GA reaches the minimum value 4.02 after the 104 generations.It shows that IQGA has faster convergence speed and better optimization result than QGA and GA, which is more suitable for solving the optimization problem of EMS.The fuzzy control membership function optimized by IQGA is shown in Figure 14.
Figures 15 and 16, respectively, reveal the distribution of IGPU operating points before and after IQGA optimization.It can be seen that, compared with the original fuzzy control strategy, the optimized IGPU operating points are more distributed in the high efficiency region.The changing curves of battery SOC before and after optimization are shown in Figure 17.It shows that the initial SOC value is 0.7; after IQGA optimization, the final SOC value changed from 0.677 to 0.674.It is observed that Δ= -0.44%, which satisfies the SOC constraints set in (25).Figure 18 shows the fuel consumption before and after optimization.It can be seen that the fuel consumption is 8.71L before optimization, while it drops by 5.17% to 8.26L after optimization.Therefore, the results indicate that, with IQGA optimized control strategy,  the fuel consumption can be effectively reduced and the IGPU efficiency can be significantly improved.

HILS Experiment.
In offline simulation, the driver's operation habits and the real-time operation effect of control strategy are not considered.In order to validate the proposed EMS in real-world implementation, the HILS experiment based on driver and controller is carried out with dSPACE.Using dSPACE, the Simulink models could be converted into C codes and run reliably [42,43].The schematic diagram  of HILS platform is shown in Figure 19.The driver controls the acceleration or brake pedal to produce analog signals which are transmitted to the general controller through the A/D interface.After calculation and analysis, the control signal produced by general controller is input to dSPACE through the CAN bus.According to the control instruction, the engine-generator model, battery model, motor model, and vehicle dynamics model are run in dSPACE.Meanwhile, the operating parameters are displayed on the monitor.The host computer is used to modify, compile, and download the internal models of dSPACE.The hardware of HILS platform is shown in Figure 20.
According to the driver's driving habits, the HILS experiment is carried out under a typical 70 seconds driving condition, and the results are shown in Figure 21. Figure 21(a) shows the driver's operation signal, and Figure 21(b) shows the vehicle speed.As we can see, in 1-11s, the vehicle experienced a slight acceleration, and the speed increased from 0 to 11 km/h.In 12-17s, the vehicle experienced a slight braking, and the speed fell to 2 km/h.In 18-33s, the vehicle experienced a medium strength acceleration, and the speed rose to 47 km/h.In 34-41s, the vehicle experienced a medium strength braking, and the speed fell to 14 km/h.In 42-62s, the vehicle experienced two types of quick acceleration, and the speed rose to 88 km/h.In 63-65s and 69-70s, the vehicle was in the state of gliding, and the speed decreased for the driving resistance.In 66-68s, the vehicle experienced an emergency braking, and the speed fell to 3 km/h.It can be seen that the vehicle speed curve is smooth and stable, which can follow the driver's acceleration or braking intention very well.
Figure 21(c) shows the motor demand power   , IGPU output power   , and battery output power   .It can be seen that the battery plays a good role in reducing peak and filling valley: when the demand power is negative, the energy can be recovered effectively; when the demand power is high, the compensation power can be output in time.Meanwhile, the engine works in the efficient area above 50KW, and when the demand power is small, it can be properly adjusted to output larger power to charge the battery.
In Figure 21(d), the torque distribution of the front and rear axle is shown.When the vehicle is driven, the driving torque on the front axle is the same with that of rear axle.When the vehicle is braking, the braking torque on the front and rear axle can be distributed according to the I curve.
Figures 21(e) and 21(f) show the distribution of mechanical torque and motor torque on front and rear axle, respectively.It can be seen that when the vehicle is driven, the motor torque can basically meet the demand torque, but when the vehicle speed is high, the motor torque is not sufficient enough for accelerating.When the vehicle is braking, it needs to be analyzed in different situations.In 12-17s, because the speed is lower than 12 km/h, only mechanical braking is carried out.In 34-38s, the braking strength is less than 0.7, and the require torque is smaller than the maximum braking torque of motor, so only regenerative braking is carried out.In 39-41s, the speed is reduced from 24 km/h to 14 km/h, so the regenerative braking torque gradually becomes smaller and the mechanical braking torque gradually becomes larger.In 66-68s, the vehicle is in an emergency braking state, and the braking strength is more than 0.7, so the braking torque is totally provided by the mechanical brake.

Conclusions
In this paper, an energy management strategy combining fuzzy control and threshold control was proposed for a series hybrid electric rescue vehicle, and a FLC optimization method based on IQGA was introduced to further improve the fuel economy.The effectiveness and real-time performance of the proposed EMS were validated by simulation and HILS experiment.The simulation results show that, compared with QGA and GA, the IQGA optimization algorithm has faster convergence speed and better optimization result.With IQGA optimization EMS, the engine efficiency is obviously improved, and the fuel consumption can be reduced by 5.17% compared with that before optimization.The HILS experiment shows that, with the proposed EMS, the real-time requirements of the vehicle can be well satisfied, the working characteristics of IGPU and battery can be effectively exploited, and the braking energy can be fully absorbed by the front and rear motor.In conclusion, the EMS proposed in this paper has achieved good control effect, which can be used in real-time control of the series hybrid electric rescue vehicle.

Figure 5 :
Figure 5: Membership function of the FLC.

Figure 9 :
Figure 9: Simulation model of series hybrid electric rescue vehicle.

Figure 10 :
Figure 10: The equivalent circuit of engine-generator set.

Figure 11 :
Figure 11: The equivalent circuit of battery.

Figure 13 :
Figure 13: Optimization process of different algorithms.

Figure 21 :
Figure 21: The results of HILS experiment.

Table 1 :
Specific parameters of the hybrid rescue vehicle.
represents the decimal value corresponding to the m-th gene   = (   −1 ⋅ ⋅ ⋅  2  1 ), which is in the range between  max and  min  ;   can be converted to   as follows: