# Design of Real-Time Control Based on DP and ECMS for PHEVs

Authors: Wei Wang; Zhenjiang Cai; Shaofei Liu
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

A real-time control is proposed for plug-in-hybrid electric vehicles (PHEVs) based on dynamic programming (DP) and equivalent fuel consumption minimization strategy (ECMS) in this study. Firstly, the resulting controls of mode selection and series mode are stored in tables through offline simulation of DP, and the parallel HEV mode uses ECMS-based real-time algorithm to reduce the application of maps and avoid manual adjustment of parameters. Secondly, the feedback energy management system (FMES) is built based on feedback from SoC, which takes into account the charge and discharge reaction (CDR) of the battery, and in order to make full use of the energy stored in the battery, the reference SoC is introduced. Finally, a comparative simulation on the proposed real-time controller is conducted against DP, the results show that the controller has a good performance, and the fuel consumption value of the real-time controller is close to the value using DP. The engine operating conditions are concentrated in the low fuel consumption area of the engine, and when the driving distance is known, the SoC can follow the reference SoC well to make full use of the energy stored in the battery.


---

## Body

## 1. Introduction

Hybrid electric vehicles use at least two power sources, usually driven by an internal-combustion engine associated with a motor, in order to minimize the fuel consumption and/or emissions. The energy management of a PHEV is often divided into two categories. The first concerns global optimization based on offline simulation. In this case, the vehicle speed is regulated to follow a speed cycle using a torque at the wheel controller. Examples of such methods include Pontryagin’s minimum principle [1, 2], dynamic programming (DP) [3–7], and genetic algorithm [8]. A second class of algorithms is real-time optimal control strategy that can be used to control a vehicle. Several algorithms have been proposed, some of which are based on rulers [9, 10] and Equivalent fuel Consumption Minimization Strategy (ECMS) [11–16], and others are approximate real-time control strategies based on DP [17–19]. ECMS has strong dynamic adaptability and can get similar results with DP in theory [20]; therefore, it has been extensively studied.In this paper, a real-time control for PHEV based on DP and ECMS is studied. Real-time implementation has remained a major challenge in the design of complex control systems. To address this hurdle, simple and efficient models and fast optimization algorithms are developed. The real-time controllers must be simple in order to be implemented with limited computation and memory resources. Moreover, manual tuning of control parameters should be avoided to reduce the calibration work efforts. DP can obtain global optimal solutions, and ECMS can realize real-time computing and can theoretically get similar results with DP. This study combines the advantages of both to establish a real-time controller.The contribution of the paper is to use the DP algorithm solving the optimal controls of driving cycle to establish the framework of FEMS. In order to fully utilize the potential of the battery, the charge and discharge reaction (CDR) of the battery is taken into account in the DP-based FEMS, and the reference SoC is introduced into the FEMS. The ECMS real-time algorithm is used for the parallel HEV mode to reduce the application of maps and avoid manual adjustment of parameters.


## 2. Hybrid Vehicle Modeling

For this study, two levels of modeling are considered. The first, called plant model (PM), shown in Figure1, is used to simulate the vehicle over speed cycles [21]. It only represents the longitudinal behavior and is designed for the energetic-consumption simulation. It includes the following:Dynamic response of engine torqueMotor model based on the characteristic map provided by the motor supplierDedicated hybrid transmission (DHT) model (including the shift strategy)Full dynamic vehicle modelHigh-voltage lithium battery model based on battery charge and discharge characteristicsFigure 1
Simulation model of PHEV in MATLAB-Simulink.An important part of PM is the fuel consumption model of engine. This is done only for fuel consumption using classical map and is validated according to real data results, as shown in Figure2.Figure 2
Fuel consumption map.Based on this PM, a simplified model, called Energy Consumption Model (ECM), has been derived. The purpose of this paper is not the vehicle modeling, but control law synthesis. So, only ECM is used to derive the optimization algorithm. PM is omitted here, but PM is used for the simulation results at the end of this paper. Figure1 is the simulation model of PHEV.


### 2.1. Energetic Consumption Modeling

The power flows of the PHEV and connections between components are shown in Figure3. The vehicle has three energy converters, an internal-combustion engine (ENG), a drive electric motor (DEM) connected through a dedicated hybrid transmission (DHT), and a generator electric motor (GEM) as a generator connected to the engine via DHT. Both electric machines can work in both motoring and generating modes. The main component parameters of the powertrain are listed in Table 1.Figure 3
Power flows of the hybrid vehicle powertrain. The arrows indicate the direction of the positive power flow.Table 1
Powertrain parameters.
ComponentsParametersValuesENGDisplacement (L)1.5Peak power (Kw)103Maximum speed (r·min−1)5500Peak torque (N·m)220GEMPeak power (Kw)9.5Peak torque (N·m)50DEMPeak power (Kw)17Peak torque (N·m)62DHTEngine gears number (—)6Motor gears number (—)2HvbatCapacity (Ah)20Energy (kWh)2Nominal voltage (V)46Maximum current (A)400As shown in Figure3, the powertrain allows the vehicle to be driven in the following four modes:Mode 1: one-motor pure electric mode: only the DEM is connected to DHT.Mode 2: two-motor pure electric mode: the DEM and GEM are connected to DHT.Mode 3: series HEV: only the DEM is connected to DHT. The ENG and GEM work as an auxiliary power unit (APU), producing electric power.Mode 4: parallel HEV: all energy converters are connected to the DHT.The following relations can be described as shown in Figure3:(1)ωwhk=ωekigbn⋅ired=ωgemkigem⋅ired=ωdemkigbj⋅ired,(2)Pmk=Pgemk+Pdemk,Tmk=igemigbjTgemk+Tdemk,(3)Twhk=iredηredηgbTekigbnk+Tmkigbjk,where n and j correspond to the engine transmission gear and the motor transmission gear, respectively.


## 2.1. Energetic Consumption Modeling

The power flows of the PHEV and connections between components are shown in Figure3. The vehicle has three energy converters, an internal-combustion engine (ENG), a drive electric motor (DEM) connected through a dedicated hybrid transmission (DHT), and a generator electric motor (GEM) as a generator connected to the engine via DHT. Both electric machines can work in both motoring and generating modes. The main component parameters of the powertrain are listed in Table 1.Figure 3
Power flows of the hybrid vehicle powertrain. The arrows indicate the direction of the positive power flow.Table 1
Powertrain parameters.
ComponentsParametersValuesENGDisplacement (L)1.5Peak power (Kw)103Maximum speed (r·min−1)5500Peak torque (N·m)220GEMPeak power (Kw)9.5Peak torque (N·m)50DEMPeak power (Kw)17Peak torque (N·m)62DHTEngine gears number (—)6Motor gears number (—)2HvbatCapacity (Ah)20Energy (kWh)2Nominal voltage (V)46Maximum current (A)400As shown in Figure3, the powertrain allows the vehicle to be driven in the following four modes:Mode 1: one-motor pure electric mode: only the DEM is connected to DHT.Mode 2: two-motor pure electric mode: the DEM and GEM are connected to DHT.Mode 3: series HEV: only the DEM is connected to DHT. The ENG and GEM work as an auxiliary power unit (APU), producing electric power.Mode 4: parallel HEV: all energy converters are connected to the DHT.The following relations can be described as shown in Figure3:(1)ωwhk=ωekigbn⋅ired=ωgemkigem⋅ired=ωdemkigbj⋅ired,(2)Pmk=Pgemk+Pdemk,Tmk=igemigbjTgemk+Tdemk,(3)Twhk=iredηredηgbTekigbnk+Tmkigbjk,where n and j correspond to the engine transmission gear and the motor transmission gear, respectively.


## 3. Optimal Control Problem

The objective in energy management for hybrid vehicles is to minimize the cumulative fuel consumption, which is equivalent to minimizing the power consumption of the engine.The battery is considered as a dynamical system, with the state of charge(4)xk+1=xk+PBTΔt,(5)xk+1=xk+ηBTPmTmk,ωdemkΔt.From (1) and (3), formula (5) can be expressed as follows:(6)xk+1=xk+ηBTPmTek,ωekΔt.The objective function is(7)J=∑k=0N−1m˙fTek,ωekΔt.The speeds and torques of both engine and motor are limited by the following mechanical constraints.Constraints on speeds:(8)ωm_min≤ωm≤ωm_max,(9)0≤ωe≤ωe_max.Constraints on torques:(10)Tm_min≤Tm≤Tm_max,(11)0≤Te≤Te_max.However, the constraints on state of charge are(12)xmin≤x≤xmax,(13)xN−x0=ΔSoC.WithΔSoC, the desired electric energy consumption over the speed cycle is called overall SoC variation.The relationships between the different torques and speeds, (2)–(4), allow writing the constraints (8) and (11) as(14)Te_min′k≤Tek≤Te_max′k,where(15)Te_min′=max0,Twhk/iredηredηgb−Tm_maxkigbjkigbnk,(16)Te_max′=minTe_max,Twhk/iredηredηgb−Tm_minkigbjkigbnk.For a given gear ratioigb, Te_min′ and Te_max′ define the interval of admissible values for engine torque. Several cases may happen, as follows:Te_min′=Te_max′=0: pure electric mode-engine speed is not high enough to close the clutchTe_min′=Te_max′≥0: the desired torque Twhk should be equal to the maximum torque of the powertrain determined by igbTe_min′≥Te_max′: the desired torque Twhk is greater than the powertrain torque capabilityTe_min′≤Te_max′: desired torque Twhk can be produced by both motor and engine


### 3.1. DP Formulation

Dynamic programming (DP) is a multistep decision process, which uses Bellman’s optimal principle for making hierarchical decisions and solving optimal controls [19].For an optimal decision, regardless of the initial state and initial decision (stage cost) dxk,xk−1,i, the remaining decisions (cost-to-go) Jk−1xk−1,i must be optimal for the first decision. That is, the second section of the optimal trajectory is also the optimal trajectory.The following is the equation of the multistep decision process:(17)Jkxk=minuk,idxk,xk−1,i+Jk−1xk−1,i,where Jkxk is the optimal value function of k-stage decision process starting state xk to the end state xf and uk,i is the control strategy at starting state xk of k-stage decision process so that the state is transferred to next state. In this paper, reverse solution is used.Figure4 shows the optimal path of WLTC using DP reverse solution (Figure 4(a)) and the cumulative fuel consumption of the corresponding optimal path (Figure 4(b)).Figure 4
The optimal results of WLTC using DP.
(a)(b)


### 3.2. ECMS Formulation

After dividing byηeqLHV, this results in the following objective function:(18)J=1ηeqLHV∑k=0N−1PeTek,ωekΔt.Introducing the Lagrangian parameterλk, the Hamiltonian function can be written as(19)H=1ηeqLHVPeTek,ωek+λkηBTTmkωdemk.In order to avoid exceeding the boundary value of the constraint condition, introducing an additional cost function, then, (19) can be rewritten as(20)H=1ηeqLHVPeTek,ωek+λk+γkηBTTmkωdemk,where(21)γk=0if constraints are not active,−Kif upper constraints are active,Kif upper constraints are active.In order to make the SoC meet the constraint condition (12), a penalty function is introduced:(22)pxk=1−xk−xxfkxmax−xmin/2a.Then, the Hamiltonian function can be rewritten as(23)H=1ηeqLHVPeTek,ωek+λk+γkηBTTmkωdemkpxk,where(24)sk=λk+γk.sk is the equivalent factor.According to Pontryagin’s minimum principle, the optimal controls are obtained by solving the minimum value of Hamiltonian function, shown as follows:(25)Tm∗=argminHsk,xk,Tmk,ωdemk.


#### 3.2.1. ECMS Algorithm

The following steps must be executed to implement ECMS, as also illustrated in Figure5:(1)
Identify the acceptable range of controlTgem,mink,Tgem,maxkand Tdem,mink,Tdem,maxkwhich satisfies the instantaneous torque constraints(2)
Discretize the intervalsTgem,mink,Tgem,maxkand Tdem,mink,Tdem,maxk into a finite number of controls Tdm,i and Tgm,i, where i=1,2,…,q and j=1,2,…,p, a total of q×p control candidates(3)
Calculate the equivalent fuel consumptionH corresponding to each control candidate(4)
Select the control valuesTgemk and Tdemk that minimize HFigure 5
ECMS algorithm flow.Steps 1 to 4 are computed at each instant of time over the entire duration of the driving cycle. This approach has been shown to closely approximate the global optimal solution.


## 3.1. DP Formulation

Dynamic programming (DP) is a multistep decision process, which uses Bellman’s optimal principle for making hierarchical decisions and solving optimal controls [19].For an optimal decision, regardless of the initial state and initial decision (stage cost) dxk,xk−1,i, the remaining decisions (cost-to-go) Jk−1xk−1,i must be optimal for the first decision. That is, the second section of the optimal trajectory is also the optimal trajectory.The following is the equation of the multistep decision process:(17)Jkxk=minuk,idxk,xk−1,i+Jk−1xk−1,i,where Jkxk is the optimal value function of k-stage decision process starting state xk to the end state xf and uk,i is the control strategy at starting state xk of k-stage decision process so that the state is transferred to next state. In this paper, reverse solution is used.Figure4 shows the optimal path of WLTC using DP reverse solution (Figure 4(a)) and the cumulative fuel consumption of the corresponding optimal path (Figure 4(b)).Figure 4
The optimal results of WLTC using DP.
(a)(b)


## 3.2. ECMS Formulation

After dividing byηeqLHV, this results in the following objective function:(18)J=1ηeqLHV∑k=0N−1PeTek,ωekΔt.Introducing the Lagrangian parameterλk, the Hamiltonian function can be written as(19)H=1ηeqLHVPeTek,ωek+λkηBTTmkωdemk.In order to avoid exceeding the boundary value of the constraint condition, introducing an additional cost function, then, (19) can be rewritten as(20)H=1ηeqLHVPeTek,ωek+λk+γkηBTTmkωdemk,where(21)γk=0if constraints are not active,−Kif upper constraints are active,Kif upper constraints are active.In order to make the SoC meet the constraint condition (12), a penalty function is introduced:(22)pxk=1−xk−xxfkxmax−xmin/2a.Then, the Hamiltonian function can be rewritten as(23)H=1ηeqLHVPeTek,ωek+λk+γkηBTTmkωdemkpxk,where(24)sk=λk+γk.sk is the equivalent factor.According to Pontryagin’s minimum principle, the optimal controls are obtained by solving the minimum value of Hamiltonian function, shown as follows:(25)Tm∗=argminHsk,xk,Tmk,ωdemk.


### 3.2.1. ECMS Algorithm

The following steps must be executed to implement ECMS, as also illustrated in Figure5:(1)
Identify the acceptable range of controlTgem,mink,Tgem,maxkand Tdem,mink,Tdem,maxkwhich satisfies the instantaneous torque constraints(2)
Discretize the intervalsTgem,mink,Tgem,maxkand Tdem,mink,Tdem,maxk into a finite number of controls Tdm,i and Tgm,i, where i=1,2,…,q and j=1,2,…,p, a total of q×p control candidates(3)
Calculate the equivalent fuel consumptionH corresponding to each control candidate(4)
Select the control valuesTgemk and Tdemk that minimize HFigure 5
ECMS algorithm flow.Steps 1 to 4 are computed at each instant of time over the entire duration of the driving cycle. This approach has been shown to closely approximate the global optimal solution.


## 3.2.1. ECMS Algorithm

The following steps must be executed to implement ECMS, as also illustrated in Figure5:(1)
Identify the acceptable range of controlTgem,mink,Tgem,maxkand Tdem,mink,Tdem,maxkwhich satisfies the instantaneous torque constraints(2)
Discretize the intervalsTgem,mink,Tgem,maxkand Tdem,mink,Tdem,maxk into a finite number of controls Tdm,i and Tgm,i, where i=1,2,…,q and j=1,2,…,p, a total of q×p control candidates(3)
Calculate the equivalent fuel consumptionH corresponding to each control candidate(4)
Select the control valuesTgemk and Tdemk that minimize HFigure 5
ECMS algorithm flow.Steps 1 to 4 are computed at each instant of time over the entire duration of the driving cycle. This approach has been shown to closely approximate the global optimal solution.


## 4. Control Design

In order to reduce the amount of memory use and improve the calculation speed, the offline simulation is used to calculate the fuel cost in series mode and mode selection for a given combinationTw,ωw,SoC [7, 19, 25]. Because the efficiency of the battery does not change greatly with the change of SoC in the desired operating region, the SoC is found to have minor effects on the optimal solution, so that effect is ignored.However, not only are all control variables stored in tables, but also some insights can be gained from the kinematic relations in (1)–(3) to reduce the amount of memory used:Mode 1 and mode 2:Tm can be directly calculated from Treq. Therefore, no tables are required.Mode 3: we only need to store the optimal line ofωAPU_opt (1-D) and Te_opt (1-D) as shown in Figure 6.Figure 6
The optimal line of APU.
(a)(b)The above maps are approximate estimates of the optimal controls of DP, which can be generated with the help of the Model-Based Calibration (MBC) toolbox of MathWorks.Mode 4 is implemented using ECMS algorithm, and the algorithm flow is shown in Figure4. There could be instances where an engine torque command produces the minimum cost but differs greatly from the previously selected engine torque. This can occur when higher engine torque and lower engine torque produce minimum costs that are close in value, which causes the Min function to alternate between higher and lower engine torque outputs. Therefore, the difference between the current engine power vector (Pe∗k) and the previously selected engine power (Pe∗k−1) is introduced into the Hamiltonian function and will help limit the rate at which the engine power (and torque) can change from time step to time step, and the Hamiltonian function (23) can be rewritten as(26)H=1ηeqLHVPeTek,ωek+skηBTTmkωdemkpxk+Pe∗k−Pe∗k−1.


### 4.1. Controller

The structure of the controller is shown in a block diagram in Figure7, which consists of three main subsystems. The first subsystem is the operation mode detection, combined with formulas (14)–(16) as the boundary condition of the mode selection; the second subsystem is operation mode management, which mainly realizes the transition of the four modes by the state machine; the third subsystem is torque distribution amnagement, which mainly realizes the torque distribution of pure electric mode (modes 1 and 2), series mode (mode 3), and parallel hybrid mode (mode 4).Figure 7
The structure of the controller.


## 4.1. Controller

The structure of the controller is shown in a block diagram in Figure7, which consists of three main subsystems. The first subsystem is the operation mode detection, combined with formulas (14)–(16) as the boundary condition of the mode selection; the second subsystem is operation mode management, which mainly realizes the transition of the four modes by the state machine; the third subsystem is torque distribution amnagement, which mainly realizes the torque distribution of pure electric mode (modes 1 and 2), series mode (mode 3), and parallel hybrid mode (mode 4).Figure 7
The structure of the controller.


## 5. Energy Management

Charge-depleting charge-sustaining strategy (CDCS) is to make use of all the stored electric energy in the battery. The PHEV is run as an electric vehicle until the SoC is under a certain limit and then operates as a hybrid in the charge-sustaining mode. It is guaranteed to make use of the stored electric energy, and it does not need information about the future driving mission, which is the main advantage of this strategy. Global optimal strategy based on DP is to mix usage of fuel and electricity throughout the driving cycle. Comparing the optimization-based strategies with the CDCS-based strategies, the optimization-based strategies may result in a lower fuel consumption than the CDCS-based strategies [24]. However, in order to use all the energy in the battery for the global optimal strategy, the distance of the driving cycle must be known.In order to make full use of the electric energy in the battery, in this paper a mix between global optimal strategy and CDCS strategies is implemented, and in order to reduce the application of maps and avoid manual adjustment of parameters, parallel HEV is implemented based on ECMS algorithm.


### 5.1. Charge-Discharge Reaction

In order to extend the life cycle of the battery, the charge and discharge reaction (CDR) of the battery is taken into account in the energy management strategy.The CDR of the battery is divided into 5 states: discharging, effective (Eff) discharging, normal, effective (Eff) charging, and critical (Crit) charging, as shown in Figure8. When the SoC is close to the maximum boundary value, the CDR is in discharging state. With the SoC gradual decrease, the CDR will be in the effective discharging state and the normal state and then in the effective charging state, and when the SoC is close to its minimum boundary value, the CDR will be in the critical charging state to avoid the voltage of the battery and the discharge depth of the battery into the nonlinear region [22, 23].Figure 8
Charge and discharge reaction.In Figure9, with SoC as the feedback variable, a feedback energy management system (FEMS) is established to maintain the SoC within an allowable interval, as shown in Figure 10(d). When the SoC decreases, the CDR also decreases accordingly; then, the FEMS will select charging maps, shown in Figure 11; when the SoC increases, the CDR also increases; then, the system will select the discharging maps. Each map is approximate estimates of the corresponding optimal trajectory of DP, which can be generated with the help of the Model-Based Calibration (MBC) toolbox of MathWorks.Figure 9
Feedback energy management system.Figure 10
The simulation results of WLTC without reference SoC.
(a)(b)(c)(d)Figure 11
The mode selection for effective charging.


### 5.2. Reference SoC

In order to make full use of all the energy stored in the battery, a blended strategy that the instantaneous optimal strategy based on ECMS is combined with CDCS strategy is implemented. In order to avoid SoC not reaching the final value of the reference SoC, when the end is reached, the strategy is to underestimate the approximate distance by 15% and use it as the horizon for the blended strategy and then switch to CS mode. This is achieved by setting a reference SoC [25], xrf, which is linear in the ratio of traveled distance versus expected distance according to equation (27). Minimum xrf is set to 0.3 in order to ensure that the final SoC is 0.3. The shape of xrf is shown in Figure 12.(27)Dt=Dreal0.85Dcycle,xrf=xf−x0Dt+x0,xmin≤xrf≤xmax,where xf is the minimum reference SoC.Figure 12
The outline of the reference SoC shape.In order to improve the robustness of the system, the PI controller is designed according to the following formula:(28)CDR=Kpxrf−x+Ki∫xrf−xdt.


### 5.3. Adaptive Optimal Supervisory Control

The adaptive optimal supervisory control is designed based on SoC feedback, which is to dynamically changesk (without using past driving information or trying to predict future driving behavior) to compare SoC changes and maintain its value near the reference value [26–28].An adaptation law based on the PI controller of the type:(29)s=s0+Kpxrf−x+Ki∫xrf−xdt.In (29), s0 represents the initial value of s at time t = 0, and Kp and Ki are the proportional and integral gains of the adaptation law. The initialization of this algorithm, i.e., the choice of s0, is arbitrary, and it can be done by averaging different optimal initial values obtained offline [28, 29].


## 5.1. Charge-Discharge Reaction

In order to extend the life cycle of the battery, the charge and discharge reaction (CDR) of the battery is taken into account in the energy management strategy.The CDR of the battery is divided into 5 states: discharging, effective (Eff) discharging, normal, effective (Eff) charging, and critical (Crit) charging, as shown in Figure8. When the SoC is close to the maximum boundary value, the CDR is in discharging state. With the SoC gradual decrease, the CDR will be in the effective discharging state and the normal state and then in the effective charging state, and when the SoC is close to its minimum boundary value, the CDR will be in the critical charging state to avoid the voltage of the battery and the discharge depth of the battery into the nonlinear region [22, 23].Figure 8
Charge and discharge reaction.In Figure9, with SoC as the feedback variable, a feedback energy management system (FEMS) is established to maintain the SoC within an allowable interval, as shown in Figure 10(d). When the SoC decreases, the CDR also decreases accordingly; then, the FEMS will select charging maps, shown in Figure 11; when the SoC increases, the CDR also increases; then, the system will select the discharging maps. Each map is approximate estimates of the corresponding optimal trajectory of DP, which can be generated with the help of the Model-Based Calibration (MBC) toolbox of MathWorks.Figure 9
Feedback energy management system.Figure 10
The simulation results of WLTC without reference SoC.
(a)(b)(c)(d)Figure 11
The mode selection for effective charging.


## 5.2. Reference SoC

In order to make full use of all the energy stored in the battery, a blended strategy that the instantaneous optimal strategy based on ECMS is combined with CDCS strategy is implemented. In order to avoid SoC not reaching the final value of the reference SoC, when the end is reached, the strategy is to underestimate the approximate distance by 15% and use it as the horizon for the blended strategy and then switch to CS mode. This is achieved by setting a reference SoC [25], xrf, which is linear in the ratio of traveled distance versus expected distance according to equation (27). Minimum xrf is set to 0.3 in order to ensure that the final SoC is 0.3. The shape of xrf is shown in Figure 12.(27)Dt=Dreal0.85Dcycle,xrf=xf−x0Dt+x0,xmin≤xrf≤xmax,where xf is the minimum reference SoC.Figure 12
The outline of the reference SoC shape.In order to improve the robustness of the system, the PI controller is designed according to the following formula:(28)CDR=Kpxrf−x+Ki∫xrf−xdt.


## 5.3. Adaptive Optimal Supervisory Control

The adaptive optimal supervisory control is designed based on SoC feedback, which is to dynamically changesk (without using past driving information or trying to predict future driving behavior) to compare SoC changes and maintain its value near the reference value [26–28].An adaptation law based on the PI controller of the type:(29)s=s0+Kpxrf−x+Ki∫xrf−xdt.In (29), s0 represents the initial value of s at time t = 0, and Kp and Ki are the proportional and integral gains of the adaptation law. The initialization of this algorithm, i.e., the choice of s0, is arbitrary, and it can be done by averaging different optimal initial values obtained offline [28, 29].


## 6. Simulation Result in MATLAB-Simulink

The controller is evaluated in a closed loop together with PM, and the simulation results are compared with the global optimal results of DP offline simulation.The offline simulation results of using DP reverse to solve WLTC are presented in Figure4. Figure 4(a) shows the optimal paths with different SoC initial values, and the cumulative fuel consumption of the corresponding optimal path is shown in Figure 4(b), and the SoC constraints are(30)25%≤x≤100%,xN=75%.In Figure4, the optimal paths with different SoC initial values converge to one path at 900 s, and the fluctuation range of SoC is in a larger interval [30, 95]; the average fuel consumption of all optimal paths after the WLTC is 820 g, corresponding to the one-hundred-kilometer fuel consumption which is 4.69 L.Figure13 shows that the engine operating points are concentrated in the low fuel consumption area of the engine and the speed is in the interval [1000 r/min, 3500 r/min]. It can also be seen from Figures 10(b) and 14(a) that the engine torque is mostly concentrated around 80 Nm, and the number of engine starts with the reference SoC (27 times) is lower than the number of engine starts without the reference SoC (31 times). Figure 10(c) is the trajectory of the equivalent factor, and the overall trend of the equivalent factor is stable with the peak upward. The larger the peak value, the greater the desire for engine power. Conversely, as shown in Figure 14(b), the equivalent fuel factor decreases with the decrease of the reference SoC, the peak is down, and the smaller the peak value, the greater the desire for motor power. The resulting SoC trajectories for the tested cycle are shown in Figures 10(d) and 14(c). In Figure 10(d), the SoC fluctuation range of the tested cycle is narrower than the DP offline simulation result in Figure 4(a), which is located in the interval [61, 65]. Compared with Figure 14(c), SoC can better follow the reference SoC, and the range of SoC variation is relatively large, indicating that the energy stored in the battery can be fully utilized. As shown in Figure 10(a), the measured vehicle speed can follow the target vehicle speed very well.Figure 13
Engine operating points’ distribution of WLTC.Figure 14
The simulation results of WLTC with reference SoC.
(a)(b)(c)In order to verify the adaptability of the controller to different tested cycles, in addition to the WLTC tested cycle, two tested cycles, China Urban Driving Cycle (CUDC) and NEDC, are also selected for simulation comparison. The results for the 3 tested cycles are shown in Table2.Table 2
Driving cycles’ simulation results.
Cycle infoWLTCCUDCNEDCWithout reference SoCmf4.813.534.71xN636159non31716With reference SoCmf4.213.234.23xN312931non27613Δmf0.610.20.48Dcycle23.165.910.95In fact, in WLTC testing, the final SoC may not reach exactly the target value (75%) of DP; therefore, in order to fairly compare fuel consumption results, a linear correlation between final SoC and fuel consumption is visible, which is easily approximated by the linear expression [30].(31)mf≈mf0+σΔSoC,where mf is the actual fuel consumption, mf0 is the value that would correspond to a zero SoC variation, and σ is a curve fitting coefficient that translates ΔSoC into a corresponding amount of fuel; here, σ≈s.In Table2, the fuel consumption of the WLTC without reference SoC is 4.81 L/100 km with the final SoC 63%. After correction, the fuel consumption is 4.83 L/100 km, which is 0.14 higher than the average fuel consumption of DP simulation with the final SoC value 75%. For the 3 test cycles, the fuel consumption without the reference SoC is higher than the fuel consumption with the reference SoC; the final value of the SoC without the reference SoC is close to the target value of 75%; the final value of the SoC with the reference SoC is close to 30%.


## 7. Conclusion

This study proposes a real-time control of PHEV based on DP and ECMS. In order to fully exploit the potential of the battery, combined with the CDR and CDCS, the FEMS was established, and the controller was evaluated by closed-loop simulation. The conclusion is as follows:(1)
This study proposes a real-time control of PHEV based on DP-ECMS, which is a suboptimal solution, and the results show that the real-time controller has good control ability and better robustness, and the fuel consumption value of the real-time controller is close to the offline simulation results of DP.(2)
The engine operating points are concentrated in the low fuel consumption area of the engine, and the engine starts and stops are evenly distributed. They effectively avoid alternating output between higher and lower engine torques.(3)
When the future driving distance is unknown, the controller can make the SoC within a admissible interval, but the SoC change range is relatively small, and the system cannot make full use of the energy stored in the battery. When the future driving distance is known, the system can make the SoC better follow the reference SoC, which can make full use of the energy stored in the battery; therefore, fuel economy is effectively improved.



---

*Source: 6667614-2021-02-03.xml*