Power Management Strategy of Hybrid Electric Vehicles Based on Quadratic Performance Index

Abstract: An energy management strategy (EMS) considering both optimality and real-time performance has become a challenge for the development of hybrid electric vehicles (HEVs) in recent years. Previous EMSes based on the optimal control theory minimize the fuel consumption, but cannot be directly implemented in real-time because of the requirement for a prior knowledge of the entire driving cycle. This paper presents an innovative design concept and method to obtain a power management strategy for HEVs, which is independent of future driving conditions. A quadratic performance index is designed to ensure the vehicle drivability, maintain the battery energy sustainability and average and smooth the engine power and motor power to indirectly reduce fuel consumption. To further improve the fuel economy, two rules are adopted to avoid the inefficient engine operation by switching control modes between the electric and hybrid modes according to the required driving power. The derived power of the engine and motor are related to current vehicle velocity and battery residual energy, as well as their desired values. The simulation results over different driving cycles in Advanced Vehicle Simulator (ADVISOR) show that the proposed strategy can significantly improve the fuel economy, which is very close to the optimal strategy based on Pontryagin’s minimum principle.


Introduction
As increasingly concerning on the deterioration in air quality and decrease in petroleum resources, a great interest is shown in the development of safe, clean, and high-efficient transportation.It has been well recognized that the electric vehicle (EV), hybrid electric vehicle (HEV), and fuel cell electric vehicle (FCEV) are the most promising solution to the problem of land transportation in the future [1].Since showing improvement in fuel consumption with minimum extra cost, HEVs are considered to offer the best promise in the short to mid-term [2].In HEVs, the internal combustion engine (ICE) provides driving power together with one or more reversible energy storage systems (ESSes), such as a bank of batteries and ultra-capacitors, which are generally used as an energy buffering unit to recycle braking energy and change engine operating points as well as provide an extra degree of freedom for energy management strategies (EMSes).Undeniably, the introduction of ESSes makes driving modes more flexible but EMSes more complicated.Therefore, it is especially important to design an excellent EMS for HEV development and application [3,4].
As early as 1997, Jalil used a set of predefined rules based on the battery state of charge (SOC) and power demand to assign the power to the engine, battery, or a combination of both, for a series HEV to ensure the high efficiency of the engine and battery operation [5].Such rule-based strategies developed from heuristic ideas are widely used in HEVs [6][7][8] because they can be implemented Energies 2015, 8,[12458][12459][12460][12461][12462][12463][12464][12465][12466][12467][12468][12469][12470][12471][12472][12473] in real-time, but making rules commonly depends on engineering experience, known mathematical models, large numbers of experimental results, etc., having limited benefits for fuel economy.In order to improve the fuel economy of rule-based methods, Schouten and Duan directly adopted fuzzy rules instead of deterministic ones to improve the operation efficiency of vehicle system in 2002 and 2003 [9,10].Zhao added a fuzzy algorithm to modify the rules in 2013 [11].However, the fuel-saving potential of HEVs cannot be fully exploited because the membership functions and rules of a fuzzy controller are designed based on human expertise and heuristics.To further reduce fuel consumption, fuzzy controllers were modified by offline optimizing membership functions and rules through a particle swarm optimization algorithm [12], a genetic algorithm [13], or a machine learning algorithm [14].Alternatively, a learning vector quantization neural network [15] or a fuzzy neural network [16] is used in a fuzzy EMS to identify the driving cycle style periodically.
By contrast with above heuristic-based strategies, EMSes based on optimal control theory, such as dynamic programming (DP) and Pontryagin's minimum principle (PMP) have been investigated quite intensively in recent years.For a prior known driving cycle, DP discretizes continuous states and control values into finite grids and decomposes the overall dynamic optimization problem into a sequence of sub-problems.The cost function of each sub-problem is the fuel consumption from current step to last step.By calculating backwards along the horizon based on Bellman's Principle of Optimality, DP generated an optimal EMS of HEVs, but with a large computational load that exponentially increases as the state variables increase in number [17][18][19][20].Theoretically, if the whole driving cycle is known in advance and the performance index is defined as an integral of fuel consumption rate, the EMS obtained from DP will minimize the total fuel consumption.However, such an optimal solution is only suitable for the known driving cycle rather than other ones.The requirement for future driving demands to be known in advance leads to a real-time problem, so DP always acts as a benchmark to assess other EMSs [21,22].To apply the optimal results of DP in real-time, many attempts have been made, such as extracting rules from optimal results [23,24], modeling the power demand as a random Markov chain [25,26] or predicting future driving conditions [27][28][29].
As another popular theory-based method, PMP formulates the optimal control problem of HEVs as a two-point boundary value problem of nonlinear differential equations.In this method, the whole driving cycle still needs to be known in advance but the computation load is much less than DP [30].A study on the comparison between PMP and DP demonstrated that the optimal results generated by PMP are very close to DP [31], so PMP can also be a benchmark.Furthermore, Kim, Cha, and Peng had proved in 2011 that under the assumption that the battery SOC varies within a small range, the Lagrange multiplier λ, which is also called as a co-state, is a constant [32].On the other hand, λ can be interpreted as an equivalent factor to equate the electrical usage of a battery to virtual fuel consumption.When λ is a constant, the dynamic optimization problem is converted into a static one, and the equivalent consumption minimization strategy (ECMS) is developed [33,34].However, λ is still very sensitive to driving cycles.Much work has been reported to solve the above problem, such as developing a function to estimate equivalent factors through observing a number of optimal results calculated from DP and PMP [35], designing a driving cycle recognizer using a neural adaptive network [36], or adding an online algorithm to ECMS framework to periodically refresh equivalent factors combined with the past and predicted vehicle speed and GPS data [37].However, if these modified theory-based methods are applied to real vehicles, they are sub-optimal, complex, and time-consuming.
In summary, the rule-based strategies are suitable for real-time applications but with limited fuel economy while the optimal theory-based strategies have the real-time problem caused by two main reasons.One is that the future driving cycle (or vehicle speed commands) should be known prior to deciding the control parameters (for example, the co-state of PMP).Another is that their calculation is relatively large.In this paper, we apply linear quadratic optimal control theory to solve the power management problem of HEVs for the first time, overcoming the shortcomings of Energies 2015, 8, 12458-12473 existing theory-based strategies with little loss in fuel economy.The proposed power management strategy is termed as the quadratic performance index strategy (QPIS), whose engine power and motor power are related to current vehicle velocity and battery residual energy as well as their desired values, and independent of future driving conditions.The fuel economy of QPIS is significantly improved from two aspects: one is to average and smooth the engine and motor power to indirectly reduce fuel consumption through a designed quadratic performance index; another is to avoid the inefficient engine operation by switching control modes based on requested driving power.The simulation results over various driving cycles show that with the same weight coefficients of quadratic performance index, the QPIS has excellent control performance on vehicle drivability, SOC sustainability, especially the fuel consumptions of QPIS are very close to that of PMP.
The remainder of this paper is organized as follows.Section 2 introduces a nonlinear vehicle model as the controlled plant and linearizes this nonlinear model for deriving QPIS.The main innovation of this paper is presented in Section 3, including constructing a quadratic performance index combing with rules and deriving the control law.In Section 4, comparative simulations in Advanced Vehicle Simulator (ADVISOR) over different driving cycles, road slopes, and vehicle total masses are performed and the results confirm the good performance of QPIS.Finally, conclusions are summarized in Section 5.

Vehicle Model
In general, the electromechanical coupling systems of HEVs are classified into three categories: torque coupling, speed coupling, and power coupling systems.The QPIS proposed in this paper is suitable for power-coupling HEVs, whose configurations are depicted in Figure 1, satisfying the assumptions as follows.(a) A continuously-variable transmission (CVT) adopted in Figure 1a or a planetary gear mechanism adopted in Figure 1b (such as the Toyota Prius and Ford Escape Hybrid) makes it possible for the engine to always operate on an optimal operating line as the bold solid one plotted in Figure 2. Every engine operating point of optimal operating line is confined to a specific output torque and speed and has the minimum fuel consumption [32].In other words, for a given engine power P ice , we can determine an optimal engine operating point " T * ice , ω * ice ‰ based on this optimal operating line.Thus control variables of energy management for HEVs are reduced from the torque and speed to only the power.
(b) Through reasonably choosing battery capacity, the battery SOC mainly varies in a narrow region, for example, 0.6-0.8,so the charge-discharge characteristics of battery are almost invariable.In other words, the open-circuit voltage and equivalent internal resistance of battery can be regarded as constants [32].Every engine operating point of optimal operating line is confined to a specific output torque and speed and has the minimum fuel consumption [32].In other words, for a given engine power ice P , we can determine an optimal engine operating point * * ice ice , T      based on this optimal operating line.Thus control variables of energy management for HEVs are reduced from the torque and speed to only the power.(c) The motor driving system has sufficient capability of short-time overload (overtorque or overpower), adequate field-weaking range, and wide high-efficiency area.The motor efficiency is not sensitive to engine operating points.
(d) Because the dynamic response time of engine or motor is much shorter than that of vehicle accelerating or decelerating and battery charging or discharging, the dynamic process of engine and motor can be neglected and only static efficiency models need to be considered [22].
In this study, a pedal position is interpreted as a velocity demand v ˚.The proposed algorithm can calculate the engine power P * ice and motor power P * ess based on current velocity, desired velocity v ˚, current SOC, and desired SOC.For the vehicle in Figure 1a, the optimal operation point " T * ice , ω * ice ‰ can be determined by P * ice , then we can jointly regulate the engine throttle and transmission ratio of CVT to make the engine run on the optimal operating line and satisfy T * ice ˆω* ice " P * ice .Meanwhile, the motor speed ω * ess is determined by the CVT ratio and P * ess " T * ess ˆω* ess can be satisfied by regulating the motor torque T * ess .For the vehicle in Figure 1b, the ring gear is connected to the final drive, so current vehicle velocity dictates the speed of the motor and ring gear.By jointly regulating the generator speed and engine throttle, the engine can run on the optimal operating line and satisfy P ice " P * ice .At the same time, we can regulate the motor torque to ensure the sum of the motor and generator power is equal to P * ess .In the following, we take the hybrid system in Figure 1a as an example to derive the QPIS.The main parameters of the vehicle originated from ADVISOR are listed in Table 1.

Nonlinear HEV Model for Simulations
The MAP of Engine FC_SI41emis and its optimal operating line are shown in Figure 2. The corresponding optimal fuel consumption line is plotted in Figure 3, which shows the relation between P ice and the fuel consumption rate .m.The fuel consumption of control algorithms in simulations is the integral of .m over the whole driving cycle.

Nonlinear HEV Model for Simulations
The MAP of Engine FC_SI41emis and its optimal operating line are shown in Figure 2. The corresponding optimal fuel consumption line is plotted in Figure 3, which shows the relation between ice P and the fuel consumption rate m  .The fuel consumption of control algorithms in simulations is the integral of m  over the whole driving cycle.If the vehicle is running at the velocity v , the driving force F can be calculated by: where  is the rotating mass coefficient; m is the vehicle total mass, including the passengers and cargo; g is the gravitational acceleration constant; r f is the rolling resistance coefficient;  is the road slope;  is the air density; D C is the aerodynamic drag coefficient; f A is the vehicle frontal area, and t  is the transmission efficiency, which is the function of vehicle velocity, load torque and CVT ratio, req req 1, 0 1, 0 and req P is the requested power, which satisfies: req ice ess where ice P is the engine power and ess P is the motor power.
Multiplying Equation ( 1) by v , we will get the vehicle dynamic model: The energy storage system is a 6Ah Ni-MH battery.The battery power bat P satisfies: where E is the total energy of battery; SOC E  is the residual energy of battery; V and R are the open-circuit voltage and equivalent internal resistance, respectively [32].If the vehicle is running at the velocity v, the driving force F can be calculated by: where δ is the rotating mass coefficient; m is the vehicle total mass, including the passengers and cargo; g is the gravitational acceleration constant; f r is the rolling resistance coefficient; θ is the road slope; ρ is the air density; C D is the aerodynamic drag coefficient; A f is the vehicle frontal area, and η t is the transmission efficiency, which is the function of vehicle velocity, load torque and CVT ratio, p " # 1, P req ą 0 ´1, P req ď 0 and P req is the requested power, which satisfies: where P ice is the engine power and P ess is the motor power.
Multiplying Equation ( 1) by v, we will get the vehicle dynamic model: The energy storage system is a 6Ah Ni-MH battery.The battery power P bat satisfies: where E is the total energy of battery; E ¨SOC is the residual energy of battery; V and R are the open-circuit voltage and equivalent internal resistance, respectively [32].
Energies 2015, 8, 12458-12473 As assumptions stated above, a permanent magnet motor is selected, whose efficiency map is shown in Figure 4.The motor power P ess can be expressed as: where η m is the motor efficiency, which is the function of motor torque and motor speed, Energies 2015, 8 7 As assumptions stated above, a permanent magnet motor is selected, whose efficiency map is shown in Figure 4.The motor power ess P can be expressed as:

Linear HEV Model for QPIS
To utilize linear quadratic optimal control theory to derive control functions of QPIS, we should establish a linear model of HEV.The dot line shown in Figure 5 depicts the relationship between v and the power to overcome resistance x v  as the state variable:

Linear HEV Model for QPIS
To utilize linear quadratic optimal control theory to derive control functions of QPIS, we should establish a linear model of HEV.The dot line shown in Figure 5 depicts the relationship between v and the power to overcome resistance P r " mg f r vcosθ `mgvsinθ `1 2 ρC D A f v 3 , when θ " 0. The correlation between v and P r should be fitted as a parabola, f v 2 (the solid line, which is available in the involved power range), to obtain the linear model, so we can replace P r with f v 2 ( f " 15.57).In addition, η t should be replaced by average efficiency, η t " 0.862.Then the vehicle dynamic model can be expressed as Equation ( 6) with x 1 " v 2 as the state variable:   For the motor, the efficiency map is symmetric about the horizontal axis as shown in Figure 4, based on which the average efficiency of motor is η m " 0.8.Thus, the linear relation between the state variable x 2 " E ¨SOC and P ess is: In the above linearization process, the linear vehicle dynamic model ( 6) is obtained when θ " 0. However, in fact, vehicles usually run on slopes, which degrades the vehicle drivability and SOC sustainability.To overcome the influences caused by road slopes, two integral actions, x 3 " ´pv ˚q2 ´v2 ¯dt and x 4 " E ¨pSOC ˚´SOCqdt, are added to the HEV model as two extended state variables.As known from linear quadratic optimal control theory, the control law is the feedback of system states, so the control variables of QPIS contain not only the feedback of current states (v and E ¨SOC) but also the integrals of the deviations of actual states from commands (v ˚and E ¨SOC * ).Such integral actions can eliminate the influences of various driving conditions on the control performance of QPIS, especially overcome the influences of road slopes effectively.6) and ( 7), we can establish the linear HEV system as: .
x " Ax `Bu u `Br z ˚ (8) where A " , and B r " are the coefficient matrixes of the system; z ˚ff are the commands; v ˚is the desired vehicle velocity decided by the pedal position, and SOC ˚is the desired battery SOC (a constant that the SOC changes around to efficiently use and protect the battery).It should be noted that the aim of linearizing the original nonlinear vehicle model is just to utilize linear quadratic optimal control theory to obtain the QPIS.The HEV model to be controlled by the strategies involved in the simulations is the same nonlinear one introduced in Section 2.1.

Power Management Strategy Based on Quadratic Performance Index
To compare with the proposed strategy in this paper, PMP is briefly introduced at first.In fact, applying PMP to solve the minimum fuel consumption problem of HEVs is to search for the motor power P ess ptq to minimize the fuel consumption under a specific driving cycle.Since the driving cycle is known previously, the requested power P req ptq can be calculated based on vehicle parameters when the vehicle is running along the vehicle velocity line of driving cycle.In this way, for the given P req ptq, Energies 2015, 8, 12458-12473 which satisfies the Equation (2), the above minimum fuel consumption control problem is to calculate P ess ptq to minimize the integral performance index as: where .
m pP ice q is the fuel consumption rate of the engine when its output power is P ice .The relationship between .m pP ice q and P ice is plotted in Figure 3 based on the assumption (a) in Section 2. Meanwhile, P ess ptq satisfies Equations ( 4) and ( 5), the terminal constraint condition: and the maximum and minimum constraints of P ice and P ess .
According to PMP, the necessary condition that the solution of above optimal control problem should satisfy is to minimize the Hamiltonian: for each sampling instant [32].In Equation ( 11), λ ă 0 is called co-state and satisfies the co-state equation .
λ " ´BH B pE ¨SOCq m and are not the functions of SOC, (i.e., the assumption (b) in Section 2), λ is a constant [32].Known from the Introduction, λ can be interpreted as an equivalent factor to equate the electrical usage of a battery to the virtual fuel consumption.Therefore, an empirical value of equivalent factor can be chosen as the initial co-state.For example, 30 kWh of battery energy corresponds to 10 L gasoline; then we can calculate the empirical value that equals ´6.935 ˆ10 ´5 as the initial co-state.For a specific driving cycle and a set co-state, P ess ptq can be calculated by minimizing Equation (11) in its feasible range for each sampling instant t.If SOC ´t f ōf the battery controlled by P ess ptq satisfies ˇˇSOC ´t f ¯´SOC pt 0 q ˇˇď 0.05, P ess ptq is the result that we desire in this paper.If SOC ´t f ¯´SOC pt 0 q ą 0.05 (or SOC ´t f ¯´SOC pt 0 q ă ´0.05), a new co-state, whose absolute value is smaller (or greater) than the absolute value of empirical co-state, is chosen to repeat the above calculation process.Obviously, to achieve optimal fuel economy, λ are different for different driving cycles.Moreover, .m and dpE¨SOCq dt in Equation ( 11) are both nonlinear; thus, it is a relatively large calculation to minimize the Equation (11) in the feasible range of P ess ptq for each sampling instant t.In the following, an EMS is obtained based on the quadratic performance index to overcome the disadvantages of PMP with little loss in fuel economy.

Power-Split Strategy Based on Quadratic Performance Index
For the system (8), we try to find a linear state feedback control u " Lx to ensure the system ˚ff is set to be constant, after the dynamic regulation process is completed, the vehicle velocity is ?x 1 " v ˚, the battery residual energy is x 2 " E ¨SOC ˚, the engine power is u 1 " f pv ˚q2 {η p t , the motor power is u 2 " 0, and the corresponding outputs of the integral regulator, x 3 and x 4 , are constant and satisfies u ss " Lx ss .We define these steady states and control inputs as x ss and u ss , the deviation of actual states from steady ones as r x " x ´xss , and the Energies 2015, 8, 12458-12473 deviation of actual control inputs from steady ones as r u " u ´uss .Since t Ñ 8 , r x Ñ 0 , and r u Ñ 0 , the performance index about r x and r u: should be finite, and the solution to minimize the above quadratic performance index can be obtained based on the regulator theory of optimal control theory.The control law has the form of state feedback as: where K ptq " K ptq T is the solution of differential Riccati equation; and L ptq " ´R´1 B T u K ptq is the state feedback matrix.Known from the quadratic optimal control theory, if t f is large enough, K ptq convergences to its terminal value only when t approaches t f , and in the most time of " t 0 , t f ı , K ptq is constant and satisfies the algebraic Riccati equation: In this way, r u " Lr x and u ss " Lx ss , where L " ´R´1 B T u K, the control variable to minimize r J is a form of state feedback: u " r u `uss " Lr x `Lx ss " Lx Additionally, the system stability can be ensured because A `Bu L is a stable matrix based on linear quadratic optimal control theory.
For the system (8), the changes of B u in value caused by the energy flow direction (indicated by p and k) lead to the variation of K.For the simplification of this algorithm, B u is rewritten as the product of two matrices: Setting fl and u 1 " ηu is the engine and motor power excluding the loss power of transmission, motor, and battery, the system (8) can be rewritten as: .
where B 1 u is constant so that the variation of K caused by B u can be avoided.Accordingly, the performance index is rewritten as: where needs to ensure vehicle drivability, prevent large fluctuations of SOC, cut peaks and fill valleys of engine power, and smooth the engine power to indirectly reduce fuel consumption.
For the optimal control problem of the system (18) and performance index (19), the solution has the form of Equation ( 16): where L 1 " ´R1 ´1 B 1 T u K 1 is the matrix of state feedback and K 1 satisfies: Then, the engine power and motor power can be calculated by: where η "

Two Control Modes
To further improve the fuel economy, two rules are designed to switch control modes (electric mode and hybrid mode) based on P req for avoiding inefficient engine operation.The selection of a switch point between the two control modes is combined with engine characteristics, vehicle parameters, and the principle of benefiting fuel economy.In this paper, the switch point is P req " P 0 , which is decided by comparing the energy conversion efficiency of a vehicle propelled separately by an engine and a motor.That is, when P req " P 0 , the efficiency of vehicle propelled solely by an engine is equal to the product of the efficiency of a vehicle propelled solely by a motor and the efficiency of the battery charged by an engine [38].The two control modes are realized by choosing different weight coefficients of the quadratic performance index as follows: (1) If P req ď P 0 : the battery should provide the total driving power to avoid inefficient engine operation or recover the braking energy that will be stored in the battery.Thus, we should set γ 5 Ñ 8 to make P ice approach to zero, and γ 2 " 0, γ 4 " 0 to temporarily not consider the constraint of SOC unless SOC reaches the maximum or minimum value.
(2) If P req ą P 0 : the engine and battery should provide requested power together.The battery shares the driving energy to restrain fluctuations of P ice .Now the constraint of SOC should be necessarily involved in the performance index for battery energy sustainability.
In general, the weight coefficients of two control modes can be calibrated by a test, where the desired SOC is 0.7 and the desired vehicle velocity line is shown in Figure 6.In the calibration, to tune the weight coefficients of Q is similar to design a proportional integral (PI) controller: γ 1 and γ 2 correspond to proportional gains, γ 3 and γ 4 correspond to integral gains.In other words, increasing γ 1 and γ 3 results in faster response and smaller tracking errors of v, but paying the price of larger fluctuations of SOC and more power provided by the engine and motor that may exceed their feasible bounds.Similarly, increasing γ 2 and γ 4 quickens the response of SOC and maintains SOC nearer to 0.7 with a little influence on trajectories of v, P ice and P ess , but weakens the buffering effect of the battery.Additionally, the fluctuations ranges of P ice (or P ess ) can be restricted within the feasible bounds by appropriately increasing γ 5 (or γ 6 ) with a slight increment of tracking errors of v.According to the above changing laws about the control effect over weight coefficients, we can tune the weight coefficients of two control modes by compromising between the tracking errors of states and fluctuation ranges of control variables until obtaining a relatively better result of vehicle velocity Energies 2015, 8, 12458-12473 error, final SOC, and fuel economy.In this paper, the calibrated weight coefficients and corresponding K 1 and L 1 are given in Table 2.The test results with above weight coefficients are shown in Figure 6." 80 0 0 80 Actually, the more the detailed operational conditions are divided, the better the control effect of QPIS is obtained, but much more complex algorithms are needed, which will weaken the adaptability of QPIS.Eventually, the concrete form of QPIS is: The control system of QPIS is shown in Figure 7.It can be observed that the control variables of QPIS are simple linear functions of current system states x without future driving conditions.This interesting feature makes QPIS simply structured and particularly suitable for real-time optimization control.Furthermore, u 1 , and u 2 are switched based on P req , which is the sum of P ice and P ess calculated by QPIS.To avoid frequent switches of control functions at switch points, we necessarily add a hysteresis loop when P req " P 0 , of which the switch on and off points are P on " 8 kW and P o f f " 2 kW, respectively, and the engine power and motor power of u 1 and u 2 are limited in their respective feasible ranges.Actually, the more the detailed operational conditions are divided, the better the control effect of QPIS is obtained, but much more complex algorithms are needed, which will weaken the adaptability of QPIS.Eventually, the concrete form of QPIS is:

Simulation Results and Discussion
To verify the control performance of proposed QPIS, under different driving cycles, road slopes, and vehicle total masses are conducted in ADVISOR, where the units are redesigned to make ADVISOR, which is a simulator combined with forward and backward approaches, suitable for forward simulation.

Various Driving Cycle
During the simulations, the weight coefficients of QPIS are the same ones given in Table 2.The initial and desired SOC are set to 0.7. Figure 8 depicts the trajectories of vehicle velocity, battery SOC, engine power and motor power obtained from QPIS under the NEDC cycle.It is found that the desired velocity trajectory is properly tracked with a small average error e v (0.064 m/s).The SOC trajectory stays near 0.7 with the deviation within 0.06.When req P is increasing, but less than on P , the control function is 1 u .If req P continues to increase and becomes greater than on P , the control function is switched to 2 u .In reverse, if req P decreases, the points to switch control functions are off P .
Throughout the driving cycle, the average value of ess P is about zero and the battery as an energy buffer effectively recycles the braking energy and restrains large fluctuations of ice P .

Simulation Results and Discussion
To verify the control performance of proposed QPIS, simulations under different driving cycles, road slopes, and vehicle total masses are conducted in ADVISOR, where the units are redesigned to make ADVISOR, which is a simulator combined with forward and backward approaches, suitable for forward simulation.

Various Driving Cycle
During the simulations, the weight coefficients of QPIS are the same ones given in Table 2.The initial and desired SOC are set to 0.7. Figure 8 depicts the trajectories of vehicle velocity, battery SOC, engine power and motor power obtained from QPIS under the NEDC cycle.It is found that the desired velocity trajectory is properly tracked with a small average error (0.064 m/s).The SOC trajectory stays near 0.7 with the deviation within 0.06.When P req is increasing, but less than P on , the control function is u 1 .If P req continues to increase and becomes greater than P on , the control function is switched to u 2 .In reverse, if P req decreases, the points to switch control functions are P off .Throughout the driving cycle, the average value of P ess is about zero and the battery as an energy buffer effectively recycles the braking energy and restrains large fluctuations of P ice .
To confirm that QPIS is independent of driving cycles, more simulations are performed under five other driving cycles.The simulation results of QPIS over six driving cycles, compared with the conventional vehicles, the default energy management strategy (DEMS) in ADVISOR and PMP are summarized in Table 3, where FC cv is the fuel consumption of conventional cars.It is observed that SOC ´t f ¯of DEMS and QPIS are unequal to SOC pt 0 q.To compare the DEMS and QPIS with PMP, the deviation of SOC ´t f ¯from SOC pt 0 q is converted to virtual fuel consumption by an equivalent factor s, which can be calculated by: where E is the total energy of battery; D f is the fuel density; Q f is the low heating value; η e is the average efficiency of engine used to charge the battery; η m is the average efficiency of motor; and η cd is the average of the battery charging and discharging efficiency.The results of fuel consumption (FC) in where E is the total energy of battery; f D is the fuel density; f Q is the low heating value; e  is the average efficiency of engine used to charge the battery; m  is the average efficiency of motor; and cd  is the average of the battery charging and discharging efficiency.The results of fuel consumption (FC) in Table 3 contain the converted fuel consumption.The DEMS, which is the electric assistant control strategy based on predefined rules in ADVISOR, determines the control mode (including motor  It is amazingly found in Table 3 that even if driving cycles are changed, the control effect of QPIS with the same group of weight coefficients remains as excellent as the results under NEDC.Compared to conventional cars, the three EMSes significantly reduce the fuel consumption.Especially, the QPIS and PMP achieve higher fuel economy than DEMS.More importantly, the fuel consumption improvement of QPIS over six driving cycles is quite close (just slightly lower, i.e., less than 4.8% at worst and 0.4% at best) to that of PMP.It is demonstrated that QPIS, which inherits the advantages of real-time performance of DEMS and the improvement in fuel economy of PMP, can be applied to the power assignment for HEVs even if the future driving cycle is unknown.

Various Road Slopes and Vehicle Total Masses
The influences of θ and m on the control effect of QPIS are involved in this section because vehicles usually run on slopes and the vehicle total mass may change in practice.In ADVISOR, θ can be easily added to a driving cycle as a function of distance and m can be easily changed by setting the vehicle parameter.Tables 4 and 5 show the results of DEMS, QPIS, and PMP when the HEV runs

Figure 2 .
Figure 2. Optimal operating line of the engine.

Figure 2 .
Figure 2. Optimal operating line of the engine.
bat ą 0 ´1, P bat ď 0 (P bat ą 0 indicates that the battery is discharging and P bat ď 0 indicates that the battery is charging).

P
m  is the motor efficiency, which is the function of motor torque and motor speed,  indicates that the battery is discharging and bat 0 P  indicates that the battery is charging).

Figure 4 .
Figure 4. Efficiency map of the motor.
correlation between v and r P should be fitted as a parabola, 2 fv (the solid line, which is available in the involved power range), to obtain the linear model, so we can replace r P with 2 fv ( 15.57 f  ).In addition, t  should be replaced by average efficiency, 0.862 t   .Then the vehicle dynamic model can be expressed as Equation (6) with 2 1

Figure 4 .
Figure 4. Efficiency map of the motor.

2 fvFigure 5 .
Figure 5. Relationship between r P and 2 fv .For the battery, when the SOC changes within the interval [0.6,0.8] and bat P varies in the interval [−15 kW,0], the battery charging efficiency calculated by

Figure 5 .
Figure 5. Relationship between P r and f v 2 .

Figure 8 .
Figure 8. Simulation results of QPIS over NEDC cycle in ADVISOR To confirm that QPIS is independent of driving cycles, more simulations are performed under five other driving cycles.The simulation results of QPIS over six driving cycles, compared with the conventional vehicles, the default energy management strategy (DEMS) in ADVISOR and PMP are summarized in Table 3, where FCcv is the fuel consumption of conventional cars.It is observed that   SOC f t of DEMS and QPIS are unequal to

Figure 8 .
Figure 8. Simulation results of QPIS over NEDC cycle in ADVISOR.

Table 2 .
Weight coefficients and state feedback matrixes of two rules.
Table 3 contain the converted fuel consumption.The DEMS, which is the electric assistant control strategy based on predefined rules in ADVISOR, determines the control mode (including motor regenerating, motor only, engine only, and hybrid mode) based on predefined threshold values, Energies 2015, 8, 12458-12473 as well as actual states and driving commands.It is suitable for real-time control, but with limited fuel economy.Energies 2015, 8

Table 3 .
Simulation results of DEMS, PMP, and QPIS over various driving cycles in ADVISOR.