Energy Optimal Control Strategy of PHEV Based on PMP Algorithm

Under the global voice of “energy saving” and the current boom in the development of energy storage technology at home and abroad, energy optimal control of the whole hybrid electric vehicle power system, as one of the core technologies of electric vehicles, is bound to become a hot target of “clean energy” vehicle development and research. This paper considers the constraints to the performance of energy storage system in Parallel Hybrid Electric Vehicle (PHEV), fromwhich lithium-ion battery frequently charges/discharges, PHEV largely consumes energy of fuel, and their are difficulty in energy recovery and other issues in a single cycle; the research uses lithium-ion battery combined with super-capacitor (SC), which is hybrid energy storage system (Li-SC HESS), working together with internal combustion engine (ICE) to drive PHEV. Combined with PSO-PI controller and Li-SC HESS internal power limited management approach, the research proposes the PHEV energy optimal control strategy. It is based on revised Pontryagin’s minimum principle (PMP) algorithm, which establishes the PHEV vehicle simulation model through ADVISOR software and verifies the effectiveness and feasibility. Finally, the results show that the energy optimization control strategy can improve the instantaneity of tracking PHEVminimum fuel consumption track, implement energy saving, and prolong the life of lithium-ion batteries and thereby can improve hybrid energy storage system performance.


Introduction
In 2015, the incident that the illegal use of emissions control software of Volkswagen further triggers a strong concern on cars "energy saving" problem.So it is essential to use energy storage technologies to alleviate energy waste.Considering the key factor that restricts battery performance and service life in lithium-ion battery, this paper applies Li-SC HESS to PHEV.However, how to coordinate each of storage units charging/discharging and optimize energy management of vehicle power system in real time is the key to implement the performance optimization of Li-SC HESS and minimize fuel consumption in PHEV.Santucci et al. [1] proposes a new dynamic optimization method by combining model predictive control with DP algorithm and considering the state of charge (SOC) each of SC and lithiumion battery, together with simplified battery aging models and others.Masih-Tehrani et al. [2] considers factors like fuel consumption, periodic battery replacement, and others and proposes a DP algorithm based on PHEV management costs, which considerably improves energy economy.But DP algorithm has a larger calculation and is more timeconsuming, so it is difficult for practical application, while the energy optimal control strategy based on PMP algorithm is becoming research focus in recent years with the advantages of its fast speed and the smaller amount of calculation than DP algorithm.
In this paper, PMP global optimization control algorithm has been adopted, which combines with Li-SC HESS internal power limited management strategy, thus taking an energy optimal control to PHEV, having an optimal management of lithium-ion battery charging/discharging states, improving the Li-SC HESS performance, and meanwhile ensuring that the vehicle fuel consumption during running can track the fuel consumption minimum trajectory in real time.S-PHEV, wherein the SHEV control system is not considered because of its relatively simple control and its large energy conversion, which easily poses a threat to lithium-ion batteries life.PHEV can make energy more efficient and fuel economy relatively higher, so in order to reduce fuel consumption as much as possible this paper selects PHEV powertrain structure, with its power-type composition architecture being shown in Figure 1, and the corresponding power flows of PHEV power system are shown in Figure 2. Its main components include gearbox, clutch, mechanical coupling, reducer, alternator, ICE, fuel tank, hybrid energy storage system, and other modules.The vehicle is equipped with a 5-speed manual gearbox; the AC alternator is connected to ICE which owns a fixed speed ratio, to reduce the weight and volume of mechanical transmission and simplify the structure.

PHEV Driven Model
From Figure 1,   ,  ice , and  alt are, respectively, the wheels, ICE, and alternator torque, Nm;   ,  ice , and  alt are the wheels, ICE, and alternator speed, rpm;  gb is the gearbox efficiency;  is the number of gearbox's gears that corresponds to its gear ratio ();  ℎV is fuel energy density, kJ/kg.
The vehicle power flow of PHEV power system is shown in Figure 2, the bidirection arrow indicates energy flow being bidirectional.Under the power-train operating mode in PHEV drive cycle, one part of power comes from energy conversion of ICE; the other part is derived from Li-SC HESS power driving.Considering the constraints of traffic environment and road situation, the electric vehicles are often in the state of standstill, started, acceleration or smooth driving, and deceleration.If lithium-ion batteries and SC output power can be allocated reasonably, Li-SC HESS can effectively reduce fuel consumption and improve PHEV economic performance.When vehicle starts, in order to reduce energy consumption of ICE, SC provides high instantaneous electric power by electric supercharging device firstly, reducing lithium-ion lifetime impact because of its high current discharging.When PHEV drives in steady, lithium-ion batteries, as the main power of Li-SC HESS, supply the portion of the power to alternator.When the vehicle accelerates or drives uphill, in order to meet the demand of high power portion, flat lithium-ion battery charging/discharging process, SC fills up the vacancy of alternator power gaps peak change preferentially.When the vehicle brakes or drives downhill, SC and lithium-ion batteries recycle braking energy at the same time, and SC takes priority to recycle the energy of high current power at the beginning.

Li-SC HESS Model.
As PHEV energy storage device, Li-SC HESS owns its control module which is also one  part of the whole vehicle control.Selecting Li-SC HESS appropriate equivalent circuit model, which includes its topology, the equivalent circuit model of each storage element, and control variables, is the precondition for the study on PHEV power system energy optimization control strategy.When the vehicle is in state of accelerating or moving uphill, this structure can quickly respond to the demand of alternator's high instantaneous power and it can also recover energy when the vehicle is in state of decelerating and braking.Besides, with respect to the topology structure between lithium-ion batteries and SC, each Li-SC HESS topology will make a difference in its component parts, different practical applications, and different connections, like the 4 kinds as follows: SC and lithium batteries parallel directly, SC and lithium-ion batteries parallel through each independent DC/DC converter, lithium-ion batteries (/SC) connects to DC/DC, and then parallel with SC (/lithium-ion batteries).Considering the advantages of HESS combination type, energy conversion efficiency, the complexity of system structure, and others, SC connects to DC/DC converter and then parallels with lithium-ion batteries (Figure 3), which is more suitable for the requirements of current dynamics and the low cost of PHEV.Compared with the independently configured DC/DC converters in parallel, its control method is relatively simpler, and the efficiency is relatively higher, which can better meet the needs of reducing vehicle fuel consumption in this paper.

Lithium-Ion Battery Equivalent
Model.This paper selects first-order RC parallel circuit as lithium-ion battery equivalent model (Figure 4), due to its low cost and less calculation; now the model has been accepted by some lithium-ion battery businesses, such as Sony and Panasonic.
Under different road conditions, the current will change dramatically with different power demands of the alternator; environmental factors such as battery temperature will also be affected, thus forming a threat to the lithium-ion battery life and safety.From Figure 4, lithium-ion battery state of charge SOC bat () is represented as follows:

𝑚 ⋅ 𝑐
From ( 1),  bat.0 is battery capacity; ( bat ()) and ( bat ()) are battery impact factor of charging/discharging rate and thermal effect, which, respectively, represent the battery current  bat () and temperature  bat (); SOC bat,0 is the initial value of battery state of charge (SOC);  stands for the battery quality;   represents battery's specific heat capacity;  bat is the terminal voltage of batteries' equivalent model; ℎ  means the convective heat transfer coefficient;  is the equivalent surface area;   is battery's ambient temperature.
From (1), SOC bat () is related to  bat () and  bat () in the current time, as  bat () is a function of  bat ().So SOC bat () can be eventually represented by  bat (), as with formula (2) as follows: Similarly, the instantaneous power  bat () and the corresponding current  bat (SOC bat ()) are shown in formula (3), wherein  bat is the battery equivalent resistance. (3)

SC Equivalent Model.
Similarly, considering the structural flexibility and lower computing costs, choosing the equivalent first-order RC parallel circuit structure as SC equivalent model, like Figure 5, the physical significance of each parameter is clear enough.Under the PHEV power systems complex conditions, this model can well describe large current and voltage fluctuations, which also has a relatively high degree of practical engineering.
From Figure 5, the voltage  sc () of  sc is similar to  bat () above; the expression SOC sc () of SC, instantaneous power  sc (), and the current  sc (SOC sc ()) are shown in formula (4) as follows: Where  sc is SC equivalent resistance, which is composed of the equivalent circuit of  ep sc paralleling with  sc and then connecting with  es sc . sc,max is the maximum voltage of  sc . sc () is the instantaneous power of SC.SOC sc,0 is the initial value of SOC.

Problem Description of PHEV Energy Management Optimization
where  is the end time in calculation,  fuel ( ice (),  ice ()) is the instantaneous fuel power at time , and the corresponding speed and torque of ICE are  ice () and  ice ().

Objective Function
From (6), () is the number of vehicle gear, the driving cycle is usually defined by   () and ().When   () and () are already known, it is pretty easy to deduce the kinetics equation (5) of ICE fuel consumption function at the desired wheel torque   ().

ICE and Alternator Physical Model Control Constraints.
ICE is a complex system; many of its physical phenomena are not easy to be modeled, like the burning process.In this case, the research ignores the temperature dependence of ICE and its dynamic characteristics, and then get the distribution graph of instantaneous fuel consumption function under the action of  ice () and  ice () with the static look-up table (LUT) showing as Figure 6.Similarly, as given  alt and  alt , combined with the relevant alternator LUT method, efficiency function of electric motor and the corresponding maximum current graph can be derived from Figures 7 and 8.
As it can be seen from Figures 6-8, when ICE and alternator speed are given, the corresponding torques are   constrained by their maximum available torque.So formula (7) is defined as follows: Formula (6) shows that the spindle torque  ps () =  ice ()+ alt (), considering the alternator torque constraints, at any time , both of the ICE minimum  ice,min ( ice ()) and maximum  ice,max ( ice ()) satisfies formula (8).Since lithium-ion battery and SC are all energy buffer devices, in the continued charging state, assessing fuel economy of the energy optimization control, PMP algorithm should meet the need of the end constraints conditions of ΔSOC bat ≈ 0 and ΔSOC sc ≈ 0, that is, formula (10).
From (11), In order to connect the dynamic characteristics of both lithium-ion battery and SC, this research considers SC characteristics that it has the priority to respond large current changes quickly and its protection to overcharge and overdischarge, and then it brings in a dynamic buffer variable Φ(SOC sc ) as a penalty function to restrain Li-SC HESS dynamic processes, which are described in the formula (12).

Solution of Extreme Value of Li-SC HESS Output Coefficient. According to Hamiltonian function equation (
where Ω = { ice ∈ ( ice,min ( ice (),  ice,max ( ice ())));  DC, ∈ ( DC, ,  DC, )} is the capacities of control variable  ice and  DC, .From (13)-( 15) it shows that solving the optimal control problem is transformed into solving the initial conditions SOC bat,0 and SOC sc,0 of Li-SC HESS state of charge and costate variable initial value  0 = ( 10 ,  20 ,  0 ).As SOC bat,0 , SOC sc,0 can be given directly; then it is further simplified into solving the initial output coefficients  10 and  20 of each energy storage element and the penalty intensity factor  0 under the constraints of the boundary condition in the vehicle driving cycle.Besides, the initial value of the costate  0 requires the minimum value of the control variables  ice and  DC, within the allowable range Ω.So define  1 ≜ − 1 / bat (SOC bat ) bat,0 ;  2 ≜ − 2 /SOC sc  sc , the Hamiltonian mathematical model of the system can be expressed as formula (16).
(16) From ( 16),  fuel ( ice ,  ice ),  bat, (SOC bat ), and  sc, (SOC sc ,  DC, ), respectively, correspond to ICE fuel consumption power, internal lithium-ion battery power, and SC power at current time, so it is the same with  1 ,  2 , and   as the weighting factor of Li-SC HESS.It is apparent that the practical significance of Hamiltonian function is described as the equivalent fuel power function, it is also to be the sum of PHEV weighted power within a certain drive period, which is consistent with the law of conservation of energy.So it verifies the feasibility of PMP algorithm in the application of the actual object.So PMP algorithm can be transformed into an online "-control" method under the optimal solution of minimum fuel consumption.

"𝜆-Control" Based PSO-PI Real-Time Optimization.
Although the costate variable initial value  0 is a constant in off-line, HESS output power differs from different cycle driving conditions in drive cycle.So using the PMP algorithm is difficult to ensure the real-time characteristic of "-control."In order to respond to the online noncausal optimal control strategy, this research uses PSO (Particle Swarm Optimization) algorithm to optimize the parameters of PI closed-loop controller (PSO-PI controller), which is to improve the characteristics of flexibility and adaptability of feedback closed-loop controller, besides its robustness and fast convergence speed, easily implementing and high computational efficiency.
Considering SOC() of each energy storage unit of Li-SC HESS, the research assumes that the reference value of Li-SC HESS SOC() is SOC ref , under the computer control system environment; set the sampling period  and  =  and introduce the PI feedback closed-loop equation ( 17), corresponding to the PSO-PI control block diagram (Figure 9).(17) From ( 17), there are two parameters   and   of PI controller to be optimized.According to the ITAE (integral of time multiplied by the absolute value of error) indicators, considering the steady-state error, the performance of settling time, small overshoot, and oversmooth, it uses criterion ITAE of PSO algorithm to calculate the objective function.In addition, every potential optimal solution of optimization problems to be optimized in PSO algorithm represents a particle in one of the solvable space, such as particle  which corresponds to the fitness value of -particle fitness function.This research introduces the particle current position   = ( 1 ,  2 , . . .,   ),  = 1, 2, . . ., , current speed ]  = (] 1 , ] 2 , . . ., ]  ), all particles' best flying position trajectory   = ( 1 ,  2 , . . .,   ), monomer extreme value  best, = ( best,1 ,  best,2 , . . .,  best, ), group extreme value  gbest, = ( gbest,1 ,  gbest,2 , . . .,  gbest, ), and inertia weight ℎ and then updates and iterates the particle according to formula (18).
From ( 18),  = 1, 2, . . ., ; ℎ is the inertia weight;  1 and  2 are the random numbers between (0, 1);  1 and  2 are the nonnegative constant evolution factor;  , and ] , are the updated position and speed of particle  at the th iteration among -dimensional space; ℎ initial is the initial inertia weight;  max is the maximum number of iterations; ℎ end is the inertia weight when  max .Taking ℎ initial = 0.9 and ℎ end = 0.4 to ensure a strong initial global search capability, the latter part of the algorithm facilitates local search.
The main steps of PSO-PI controller parameter optimization are as follows: (1) Assuming that the particle  has parameters   ,   , group scale, current iteration number  and the iteration maximum number  max , inertia weight, learning factor, monomer extreme value  best, , and group extreme value  gbest, , and so on, then initialize   and ]  of particle  randomly.
(2) Update the  , and ] , of particle  according to formula (18) and then calculate its fitness value   .
(3) Compare   with the corresponding  best, ; if   >  best, , update the position of  best, instead of the current position of   .
(4) Similarly, compare   with the corresponding  gbest, ; if   >  gbest, , update the position of  gbest, instead of the current position of   .
(6) Output the optimized parameter values   and   .
Since PI controller is not adaptive itself, add PSO algorithm to adjust the parameters of the controller; the selfadaptability is improved to achieve the purpose of fast tracking and controlling the covariable in PMP algorithm.When group scale is set to 30, the maximum calculation period to is set to 100, and both  1 and  2 are set to 1.50425, the convergence curve of the best individual fitness function is shown in Figure 10.

Management of Li-SC HESS Limited Power.
Above, the research takes the vehicle dynamic character and minimal fuel consumption as the main analysis object and initially establishes energy optimization management method of PHEV power system.Although the allocated processing factors can ensure the coordinated allocation of power between lithium-ion battery and SC, SOC constraints during controlling just to prevent Li-SC HESS overcharge and overdischarge, which are the basic conditions.In order to further improve Li-SC HESS performance, it needs to make a real-time management of the power of Li-SC HESS each energy storage unit during the charging/discharging state.The vehicle energy optimization control flow chart is shown in Figure 11.
According to the lithium-ion battery's characteristics of low power density, strong energy density, and limited life, the research adopts the preresponse principle of SC when the demanded power of the alternator changes and dictates that when SOC sc of SC reaches the limiting value of serious charge/overcharge, Li-SC HESS prohibits its charge/discharge.When SOC sc has not reached serious critical limits, it will be divided into the normal working section (SOC low , SOC high ) and power limitation management section (SOC high , SOC max ) ∪ (SOC min , SOC low ).That is as follows.
Δ bat and Δ sc are the amended power of lithium-ion battery and SC, respectively, and there is Δ bat = −Δ sc .During the normal operation, SOC sc ∈ (SOC low , SOC high ), Δ sc = 0, and the power of each storage device does not change.When the discharging power exceeds the limit, the rule of power correction is shown in formula (19).Similarly, there is formula (20) when the charging power exceeds the limit.(20)

Results and Analysis
PHEV with Li-SC HESS in this paper is obtained by the secondary development of PHEV model based on ADVISOR software (Figure 12).Taking into account of the majority of domestic small car users daily and mainly using in the city, the research adopts the urban road cycling conditions (CYC-UDDS).The realtime curve of driving cycle speed is shown in Figure 13(a).The gear position curve of the driving cycle is shown in Figure 13(b).Integrating Figures 13(a) and 13(b), it can be seen that the vehicle speed of the driving cycle is welltracked with gear position changes, which basically meets the evaluation requirements of vehicles driving cycle actually.It also verifies the feasibility of energy-optimized controlling electric vehicles by PMP algorithm.
(1) The comparison results of PHEV power system before and after: when working normally in the vehicle driving cycle, Li-SC HESS can provide or absorb some of the energy through alternator, reducing fuel consumption ICE.Figures 14(a) and 14(b), respectively, represent the alternator torque curve of vehicle power system, controlled by PMP energy optimization algorithm.
Comparing Figure 14(a) with Figure 14(b), it is easy to know that, after using PMP algorithm, the alternator torque curve fluctuation changes more dramatically than before and the alternator power requirements increase significantly.According to the conservation of energy, it indicates that the part of ICE fuel consumption absorbed by alternator increases clearly.
(2) The result of real-time optimization of output coefficients by PSO algorithm: in order to demonstrate the characteristic of real-time tracking by PSO-PI controller, the research obtains the output coefficient () curve (Figure 15) of Hamiltonian function under the action of single lithiumion battery.
From Figure 15, this control strategy can also adjust ICE to moving along its track of minimum fuel energy consumption through the alternator real-timely, and it also further validates the rationality and feasibility of PMP algorithm.The simulation situation between the single lithium-ion battery energy storage and Li-SC HESS.It can be known from Figures 16(a) and 16(b) that the state of charge of lithium battery or SC is consistent with the end constraints of Pontryagin's minimum principle.Besides, as is shown in Figures 17(a The charge and discharge currents of lithium-ion battery are clearly smaller than the single lithium-ion batteries, which not only well smooth the charging/discharging process of batteries, but also obviously reduce the corresponding lithiumion batteries' charging/discharging cycle times.Finally, it reflects the good coordinate ability between Li-SC HESS each energy storage element, which can extend battery's life.And it also validates the validity and effectiveness of this energyoptimized control method for designed PHEV with Li-SC HESS.

Summary
The research designs a kind of PHEV that introduces Li-SC HESS.Compared with traditional vehicles, This PHEV has the advantages of both ICE and Li-SC HESS.For example, the internal ICE can use existing gas station resources and reduce the overall investment costs.Meanwhile, it can alleviate the difficulty more effectively than pure electric vehicles when solving the problems brought from defrosting, air conditioners, and other pieces of large energy consumption equipment.Li-SC HESS can help extend battery life and extend the driving ranges of cars.In particular, the embedment of SC makes Li-SC HESS well suited to start the vehicle, the speed change, and energy recovery during braking.Mainly, PHEV energy optimization control strategy can effectively reduce vehicle exhaust emissions, benefit for the urban environment, which has a high research value.

Figure 2 :
Figure 2: Power flow schematic of PHEV power system.

Figure 4 :
Figure 4: Lithium-ion battery equivalent model based on first-order RC parallel circuit.

Figure 5 :
Figure 5: SC equivalent model based on first-order RC parallel circuit.
Figure 6: Distribution graph of instantaneous fuel consumption function under the action of  ice (),  ice ().
. c u r r e n t (A ) 0 7500 15000 A lt. sp ee d (r ps )

Figure 7 :
Figure 7: Graph of Alt.efficiency function curve under the action of  alt ,  alt .

Figure 8 :
Figure 8: Graph of Alt.maximum current curve under the action of  alt ,  alt .

Figure 11 :
Figure 11: The overall flow chart of PHEV energy optimization control.
) and 17(b) together with Figure 16(b), since Li-SC HESS is embedded with SC in the vehicle driving cycle, it can significantly reduce the output of lithium-ion battery.
Journal of Control Science and Engineering 2.1.Vehicle Structure of PHEV.According to vehicle powertrain structure, HEV can be divided into SHEV, PHEV, and 2 Control Variable Constraints.PHEV energy optimal management requires considering the restrictions of the objective function from a global point of view, mainly from vehicle structure constraints; both physical model constraints of ICE and alternator and Li-SC HESS state constraints.