Adaptive Equivalent Consumption Minimization Strategy for Hybrid Heavy-Duty Truck Based on Driving Condition Recognition and Parameter Optimization

: The accurate determination and dynamic adjustment of key control parameters are challenges for equivalent consumption minimization strategy (ECMS) to be implemented in real-time control of hybrid electric vehicles. An adaptive real-time ECMS is proposed for hybrid heavy-duty truck in this paper. Three e ﬀ orts have been made in this study. First, six kinds of typical driving cycle for hybrid heavy-duty truck are obtained by hierarchical clustering algorithm, and a driving condition recognition (DCR) algorithm based on a neural network is put forward. Second, particle swarm optimization (PSO) is applied to optimize three key parameters of ECMS under a speciﬁed driving cycle, including equivalent factor, scale factor of penalty function


Literature Review
Heavy-duty trucks play an important role in the modernization of the national infrastructure.In the context of current heavy-duty truck market demand and strict fuel consumption and emission limits, energy saving and emissions reduction are of great significance to the development of heavy-duty trucks.In addition, the increasingly strict fuel consumption standard of the vehicle also makes it difficult to effectively resolve the contradiction between economic demands and power demands of traditional heavy-duty trucks [1,2].The development of an efficient hybrid power system is an effective way to meet the above challenges, and has become a common demand for the development of the world's heavy-duty truck industry.Establishing an appropriate energy management strategy (EMS) is the core difficulty of the design and development of hybrid heavy-duty trucks, and it is also the key to energy flow better.The motivation of this paper includes two main aspects.On one hand, for hybrid heavy-duty truck, inappropriate vehicle speed threshold for engine start-up and charge-sustainability of the battery will cause unnecessary fuel consumption, therefore, battery state of charge (SOC) and engine state should be considered.On the other hand, EF is affected by driving conditions, EF should vary in accordance with actual driving conditions, DCR algorithm should be studied.This paper tries to improve the control performance and real-time of ECMS by integrating the above two aspects, which makes this paper different from the previous contribution.

Contributions
The primary purpose of this paper is to propose a novel A-ECMS by combining DCR and optimization of three key factors for improving control performance of ECMS.Three main contributions have been made in this paper.(1) Based on driving cycles for heavy-duty vehicles applied in China and other countries, the typical sample database of driving cycles for heavy-duty truck is constructed by a hierarchical clustering algorithm, and a neural network is presented to recognize the above typical driving cycles, in preparation for incorporating driving conditions into ECMS.(2) Considering the effect of EF, battery SOC, and engine state on the control performance of ECMS, PSO is applied to optimize EF, scale factor of penalty function, and vehicle speed threshold for engine start-up simultaneously.
(3) A novel A-ECMS composed of key parameters optimization and DCR is proposed, the A-ECMS can achieve better fuel economy as well as keeping battery SOC fluctuate more smoothly.

Outline
This paper is organized as follows.The Section 2introduces the structure and system model of hybrid heavy-duty truck which is the basis of subsequent simulation.The Section 3 describes driving cycle generations based on hierarchical clustering algorithm and DCR derived from neural network.Section 4 conducts the optimization of three key parameters for ECMS by PSO.A novel A-ECMS combining above DCR and parameter optimization is presented in the Section 5. Conclusions are given in the sixth section.

Powertrain Structure and System Model
The heavy-duty truck investigated in this paper is a P2 hybrid powertrain configuration shown in Figure 1.The specifications of the hybrid heavy-duty truck are listed in Table 1.
Energies 2020, 13, x FOR PEER REVIEW 3 of 21 sustainability of the battery will cause unnecessary fuel consumption, therefore, battery state of charge (SOC) and engine state should be considered.On the other hand, EF is affected by driving conditions, EF should vary in accordance with actual driving conditions, DCR algorithm should be studied.This paper tries to improve the control performance and real-time of ECMS by integrating the above two aspects, which makes this paper different from the previous contribution.

Contributions
The primary purpose of this paper is to propose a novel A-ECMS by combining DCR and optimization of three key factors for improving control performance of ECMS.Three main contributions have been made in this paper.(1) Based on driving cycles for heavy-duty vehicles applied in China and other countries, the typical sample database of driving cycles for heavy-duty truck is constructed by a hierarchical clustering algorithm, and a neural network is presented to recognize the above typical driving cycles, in preparation for incorporating driving conditions into ECMS.(2) Considering the effect of EF, battery SOC, and engine state on the control performance of ECMS, PSO is applied to optimize EF, scale factor of penalty function, and vehicle speed threshold for engine start-up simultaneously.(3) A novel A-ECMS composed of key parameters optimization and DCR is proposed, the A-ECMS can achieve better fuel economy as well as keeping battery SOC fluctuate more smoothly.

Outline
This paper is organized as follows.The Section 2introduces the structure and system model of hybrid heavy-duty truck which is the basis of subsequent simulation.The Section 3 describes driving cycle generations based on hierarchical clustering algorithm and DCR derived from neural network.Section 4 conducts the optimization of three key parameters for ECMS by PSO.A novel A-ECMS combining above DCR and parameter optimization is presented in the Section 5. Conclusions are given in the sixth section.

Powertrain Structure and System Model
The heavy-duty truck investigated in this paper is a P2 hybrid powertrain configuration shown in Figure 1.The specifications of the hybrid heavy-duty truck are listed in Table 1.

Vehicle Model Based on AMESim
Simulation plays an important role in the verification and adjustment of energy management strategies for HEVs.An accurate vehicle model is the basis for simulations.Based on AMESim, a quasi-static backward simulation model of hybrid heavy-duty vehicle is established using a combination of theoretical modeling and experimental modeling.The modeling process of key components is described below.

Engine Model
The energy management strategy of hybrid heavy-duty vehicle is concerned with engine fuel consumption, thus transient characteristics of the engine is not the research focus of this paper.A quasi-static engine model is established by a look-up table based on experimental data from bench test.The engine fuel consumption contour map and engine maximum torque curve are shown in Figure 2. The engine fuel consumption rate is usually described as a function of engine speed and torque, as shown in Equation (1).
b e = f (n e , T e ) where T e and n e denote the engine torque and engine rotational speed respectively, and b e is the brake-specific fuel consumption (BSFC) in g/kWh.

Vehicle Model Based on AMESim
Simulation plays an important role in the verification and adjustment of energy management strategies for HEVs.An accurate vehicle model is the basis for simulations.Based on AMESim, a quasi-static backward simulation model of hybrid heavy-duty vehicle is established using a combination of theoretical modeling and experimental modeling.The modeling process of key components is described below.

Engine Model
The energy management strategy of hybrid heavy-duty vehicle is concerned with engine fuel consumption, thus transient characteristics of the engine is not the research focus of this paper.A quasi-static engine model is established by a look-up table based on experimental data from bench test.The engine fuel consumption contour map and engine maximum torque curve are shown in Figure 2. The engine fuel consumption rate is usually described as a function of engine speed and torque, as shown in Equation (1).

( , )
e e e b f n T = where e T and e n denote the engine torque and engine rotational speed respectively, and e b is the brake-specific fuel consumption (BSFC) in g/kWh.

Motor Model
The electric motor is a permanent magnet synchronous motor (PMSM) in this paper.Similar to the engine model, the dynamic characteristics of the PMSM are also neglected here, and a quasi-static PMSM model is obtained based on the PMSM efficiency map derived from bench test data.The efficiency map and external characteristic curve of PMSM are presented in Figure 3.The PMSM efficiency can be expressed as Equation ( 2), which is a function of torque and speed.

Motor Model
The electric motor is a permanent magnet synchronous motor (PMSM) in this paper.Similar to the engine model, the dynamic characteristics of the PMSM are also neglected here, and a quasi-static PMSM Energies 2020, 13, 5407 5 of 20 model is obtained based on the PMSM efficiency map derived from bench test data.The efficiency map and external characteristic curve of PMSM are presented in Figure 3.The PMSM efficiency can be expressed as Equation (2), which is a function of torque and speed.
9550η em (t) T em (t)>0 where η em is the efficiency of the motor, ω em is the angular velocity of the motor in rad/s, n em is the rotation speed of the motor in r/min, T em is the desired torque of the motor, T em_max_dis (n em (t)) is the maximum output torque of the motor at the current speed when the motor operates as a traction motor, T em_max_char (n em (t)) is the maximum output torque of the motor at the current speed when the motor operates as a generator, and P b is the terminal battery power in kW.
Energies 2020, 13, x FOR PEER REVIEW 5 of 21 ( ) 0 1000 1000 60 9550   The battery is a complex nonlinear system which may be influenced by temperature, voltage, internal resistance and state of charge (SOC), thus it is difficult to establish a battery model with high accuracy.For the study of energy management strategy, complicated battery model is neither essential nor practical, a simplified control-oriented model with adequate accuracy is adopted in this paper.Rint model, RC model, Thevenin model, and PNGV model are four categories of simplified battery model that have been widely employed in the area of energy management.The Rint model is used in this paper due to its acceptable precision and simplicity.The Rint model is an internal resistance model based on experimental data, which consists of an ideal voltage source and an internal resistor consisted in series.The current and SOC of the battery can be formulated as follows:  The battery is a complex nonlinear system which may be influenced by temperature, voltage, internal resistance and state of charge (SOC), thus it is difficult to establish a battery model with high accuracy.For the study of energy management strategy, complicated battery model is neither essential nor practical, a simplified control-oriented model with adequate accuracy is adopted in this paper.Rint model, RC model, Thevenin model, and PNGV model are four categories of simplified battery model that have been widely employed in the area of energy management.The Rint model is used in this paper due to its acceptable precision and simplicity.The Rint model is an internal resistance model based on experimental data, which consists of an ideal voltage source and an internal resistor consisted in series.The current and SOC of the battery can be formulated as follows: Energies 2020, 13, 5407 6 of 20 where V oc is the battery open circuit voltage; R b is the battery internal resistance; I b is the battery current; P b denotes the battery power; SOC init and Q max are the battery initial SOC and battery capacity, respectively.

Vehicle Longitudinal Dynamic Model
Given that the fuel economy is mainly associated with longitudinal motion of the vehicle, only longitudinal dynamics are considered in this paper, without involving vertical dynamics and lateral dynamics of the vehicle.It is assumed that the vehicle operates on a horizontal road, the longitudinal dynamics model can be described as Equation (7).
where m is the mass; f is the rolling resistance coefficient; δ is the rotating mass coefficient; ρ air is the air density; C d is the wind resistance coefficient; A is the frontal area; a is the acceleration; v is the vehicle speed; r w is the wheel radius.
The relationship between wheel driving torque and coupling torque produced by the engine and motor can be expressed by Equation (8).
where g 0 is the total transmission ratio of the transmission and the final drive, η t is the transmission efficiency, and T w is the total required torque of the engine and motor.

Construction of Typical Sample Database of Driving Cycles
Constructing a typical sample database of driving cycles is the basis of improve robustness of ECMS for hybrid heavy-duty truck.This paper presents a novel design method of driving cycles.First, the driving cycles for heavy-duty truck applied by China and other countries in Table 2 are chosen as sample database of driving cycle, which covers different driving scenarios including urban, suburban and highway.Second, in order to decrease the types and difficulty of DCR, clustering analysis is utilized to classify above sample database of driving cycles by selecting appropriate characteristic parameters.Maximum speed, average speed, average operating speed, maximum acceleration, minimum acceleration, average acceleration of acceleration segment, average deceleration of deceleration segment, acceleration proportion, deceleration proportion, parking proportion, uniform proportion, and relative positive acceleration are chosen as characteristic parameters in this paper.The effectiveness of above 12 characteristic parameters are verified by a significance test called Kruskal-Wallis test [24].Before carrying out the Kruskal-Wallis test, enough samples are necessary for it.There are only ten driving cycles in sample database of driving cycles, the samples are too small when single driving cycle is selected as the sample.Therefore, driving samples are generated by truncating randomly from above standard driving cycle in Table 2 to enlarge sample database in this paper.Based on the driving cycles in Table 2, 200 new short driving cycles are generated randomly for each driving cycle in Table 2.The generation of driving samples can be stated as follow.Given a long driving cycle {C n } of length n.Firstly, a starting position n 0 is selected randomly.
Then, the truncated segment {c t } will be where t is driving period.When the driving period of HEV approaches or exceeds three minutes, it can reflect the characteristic of current driving condition [25].As a result, the length of each sample is defined as 200 s.200 samples under each standard driving cycle in Table 2 are generated by repeating the method for the Kruskal-Wallis test.Taking the average speed for an example, the significance test can be stated as follow.
Step 1: State the hypotheses.

The Null Hypothesis 1 (H1).
The distribution of the average speed samples coming from multiple driving cycles has no significant difference.
The Alternative Hypothesis (Ha).The distribution of the average speed samples coming from multiple driving cycles has significant difference.
Step 2: Calculate the test statistic.
where N is the total number of values in all samples; k is the number of sample groups; n i is the number of values contained in the ith sample; R i is the sum of ranks in ith sample.
Step 3: Calculate the result.The significance level in this paper is 0.05.The p-value corresponds to the test statistic, which turns out to be 1.1 × 10 −290 .Since the p-value is less than 0.05, the null hypothesis is rejected and the distribution of the average speed samples coming from multiple driving cycles has significant difference.
Also, the box plots of average speed for each sample is shown in Figure 4.The middle line of a box indicates the median of sample data under every driving cycle, as well as revealing significant differences of average speed under different driving cycle, while the upper edge and the lower edge of a box indicate the upper quartile and the lower quartile respectively.The significance test of other characteristic parameters can also be conducted by Kruskal-Wallis test like above average speed, and the p-value of other characteristic parameters is far less than 0.05, indicating that the characteristic parameters selected in this paper have significant differences under different driving cycles and can be used to represent the driving cycle.Misclassification will happen when a single driving cycle is identified together with a composite one, therefore it is necessary to divide composite driving cycles in sample database of driving cycles into a single driving cycle by construction method and components specified by China automotive test cycle-Part 2: Heavy-duty commercial vehicles (GB/T 38146.2-2019),and fuel consumption test methods for heavy-duty commercial vehicles (GB/T 27840-2011).The updated sample database of driving cycles is shown in Table 3.The hierarchical clustering analysis is used to cluster the sample database of driving cycles, and the process of hierarchical clustering analysis can be shown in Figure 5. First, the 12 characteristic parameters of each driving cycle are transformed to no-dimensional forms by Z-score standardized method.Then, the similarity among different driving cycles is characterized by angle cosine distance.As a result, the sample database of driving cycles is categorized into six types by angle cosine distance as presented in Table 4. Taking China automotive test cycle for heavy-duty commercial vehicles as priority principle, six driving cycles are chosen as representative of typical driving cycles to form each of the above six types, respectively.The curves of six typical driving cycles are shown in Figure 6.

Algorithm Design of DCR Base on Neural Network
Based on above clustering analysis result, back propagation (BP) neural network is adopted to design driving condition recognizer.The on-line recognition process of driving conditions is illustrated in Figure 7.
The neural network structure is composed of input layer with 12 neurons denoting 12 characteristics mentioned above, hidden layer with 10 neurons and output layer with 6 neurons denoting six typical driving cycles.The maximum number of iterations is 100.Further, seventy percent of sample data is used for learning and training, fifteen percent of sample data is used to validate, and remaining fifteen percent of sample data is used to test.And the training result is shown in Figure 8.After 88 iterations, the error is reduced to less than 0.0018 and the neural network can be used to identify the driving conditions online.

Algorithm Design of DCR Base on Neural Network
Based on above clustering analysis result, back propagation (BP) neural network is adopted to design driving condition recognizer.The on-line recognition process of driving conditions is illustrated in Figure 7.
characteristics mentioned above, hidden layer with 10 neurons and output layer with 6 neurons denoting six typical driving cycles.The maximum number of iterations is 100.Further, seventy percent of sample data is used for learning and training, fifteen percent of sample data is used to validate, and remaining fifteen percent of sample data is used to test.And the training result is shown in Figure 8.After 88 iterations, the error is reduced to less than 0.0018 and the neural network can be used to identify the driving conditions online.The neural network structure is composed of input layer with 12 neurons denoting 12 characteristics mentioned above, hidden layer with 10 neurons and output layer with 6 neurons denoting six typical driving cycles.The maximum number of iterations is 100.Further, seventy percent of sample data is used for learning and training, fifteen percent of sample data is used to validate, and remaining fifteen percent of sample data is used to test.And the training result is shown in Figure 8.After 88 iterations, the error is reduced to less than 0.0018 and the neural network can be used to identify the driving conditions online.Based on the assumption that driving conditions will not change suddenly in a short period of time, on-line recognition of driving conditions predicts driving condition for the foreseeable future by extracting information from driving condition in the past.The principle of on-line recognition is shown in Figure 9. Sampling period and prediction period are the two most important parameters for on-line recognition process.The recognition accuracy is determined by BP neural network error, the recognition accuracy is high when the error is low.According to the relation between sampling period and recognition accuracy as shown in Figure 10, the maximum recognition accuracy is 99.99% when sampling period is 275 s.However, the amount of computations will increase while the sampling period is set to larger.When the sampling period is 200 s, there is a turning point.Based on this, sampling period is set to 200 s and recognition accuracy is 99.7% in this paper.The prediction period can be set to 1 s, 3 s, 5 s, 10 s, 15 s, 20 s, and so on, the length of prediction period will influence on-line recognition accuracy.By a novel method of determining prediction period proposed by He [26], the relation between prediction period and on-line recognition accuracy can be found in Figure 11.According to Figure 11, the prediction period is set to 5 s considering both computation time and recognition accuracy.Based on the assumption that driving conditions will not change suddenly in a short period of time, on-line recognition of driving conditions predicts driving condition for the foreseeable future by extracting information from driving condition in the past.The principle of on-line recognition is shown in Figure 9. Sampling period and prediction period are the two most important parameters for on-line recognition process.The recognition accuracy is determined by BP neural network error, the recognition accuracy is high when the error is low.According to the relation between sampling period and recognition accuracy as shown in Figure 10, the maximum recognition accuracy is 99.99% when sampling period is 275 s.However, the amount of computations will increase while the sampling period is set to larger.When the sampling period is 200 s, there is a turning point.Based on this, sampling period is set to 200 s and recognition accuracy is 99.7% in this paper.The prediction period can be set to 1 s, 3 s, 5 s, 10 s, 15 s, 20 s, and so on, the length of prediction period will influence on-line recognition accuracy.By a novel method of determining prediction period proposed by He [26], the relation between prediction period and on-line recognition accuracy can be found in Figure 11.According to Figure 11, the prediction period is set to 5 s considering both computation time and recognition accuracy.
when sampling period is 275 s.However, the amount of computations will increase while the sampling period is set to larger.When the sampling period is 200 s, there is a turning point.Based on this, sampling period is set to 200 s and recognition accuracy is 99.7% in this paper.The prediction period can be set to 1 s, 3 s, 5 s, 10 s, 15 s, 20 s, and so on, the length of prediction period will influence on-line recognition accuracy.By a novel method of determining prediction period proposed by He [26], the relation between prediction period and on-line recognition accuracy can be found in Figure 11.According to Figure 11, the prediction period is set to 5 s considering both computation time and recognition accuracy.In this paper, the test driving condition is a composite driving cycle composed of 'CHTC_HT_CITY + CHTC_HT_SUB + CYC_WVUSUB + CYC_MANHATTAN+CHTC_C_SUB + CYC_HWFET', the recognition result of test driving condition with an recognition accuracy of 86.8% is shown in Figure 12 below.

Principles of ECMS
The key function of ECMS is given by Equation ( 12):  In this paper, the test driving condition is a composite driving cycle composed of 'CHTC_HT_CITY + CHTC_HT_SUB + CYC_WVUSUB + CYC_MANHATTAN+CHTC_C_SUB + CYC_HWFET', the recognition result of test driving condition with an recognition accuracy of 86.8% is shown in Figure 12 below.

Principles of ECMS
The key function of ECMS is given by Equation ( 12): In this paper, the test driving condition is a composite driving cycle composed of 'CHTC_HT_CITY + CHTC_HT_SUB + CYC_WVUSUB + CYC_MANHATTAN+CHTC_C_SUB + CYC_HWFET', the recognition result of test driving condition with an recognition accuracy of 86.8% is shown in Figure 12 below.In this paper, the test driving condition is a composite driving cycle composed of 'CHTC_HT_CITY + CHTC_HT_SUB + CYC_WVUSUB + CYC_MANHATTAN+CHTC_C_SUB + CYC_HWFET', the recognition result of test driving condition with an recognition accuracy of 86.8% is shown in Figure 12 below.

Principles of ECMS
The key function of ECMS is given by Equation ( 12):

Principles of ECMS
The key function of ECMS is given by Equation ( 12): .
m f_eqv (t) denotes total fuel consumption; . m f (t) denotes instantaneous fuel consumption of the engine; . m b (t) denotes equivalent fuel consumption of the battery as expressed in Equation (13).
where Q lhv denotes lower heating value of the fuel; P b denotes battery power; η dis and η chg are discharge efficiency and charge efficiency of the battery, respectively; s ef denotes EF.
A penalty function is used to keep battery charge-sustainability as shown in Equation (14).
where SOC target is target value of battery SOC, which is 50% in this paper; K p is scale factor of penalty function.As a result, the key function of ECMS can be determined as Equation (15). .
Considering normal operation range of the powertrain, constraint conditions are given by (16).
where T e_max and T e_min denote engine maximum torque and engine minimum torque, respectively; w e_max and w e_min denote engine maximum speed and engine minimum speed, respectively; T m_max and T m_min are motor maximum torque and motor minimum torque, respectively; w m_max and w m_min are motor maximum speed and motor minimum speed, respectively; SOC H and SOC L are battery maximum SOC and battery minimum SOC, respectively.According to the principle of ECMS shown in Figure 13 and key functions of ECMS, EF, scale factor of penalty function, and vehicle speed threshold for engine start-up are the three most parameters affecting control performance of ECMS, which should be optimized for further improvement of ECMS.

Parameter Optimization of ECMS Based on PSO
In this paper, PSO is introduced to optimize above three key parameters.The optimization process of PSO based on co-simulation platform combined with AMESim and Matlab/Simulink (R2016a, MathWorks, Natick, MA, USA) is shown in Figure 14.Considering the SOC balance and engine state, state variables and fitness function are regarded as Equation ( 17) during the optimization process.

Parameter Optimization of ECMS Based on PSO
In this paper, PSO is introduced to optimize above three key parameters.The optimization process of PSO based on co-simulation platform combined with AMESim and Matlab/Simulink (R2016a, MathWorks, Natick, MA, USA) is shown in Figure 14.

Parameter Optimization of ECMS Based on PSO
In this paper, PSO is introduced to optimize above three key parameters.The optimization process of PSO based on co-simulation platform combined with AMESim and Matlab/Simulink (R2016a, MathWorks, Natick, MA, USA) is shown in Figure 14.Considering the SOC balance and engine state, state variables and fitness function are regarded as Equation ( 17) during the optimization process.Considering the SOC balance and engine state, state variables and fitness function are regarded as Equation ( 17) during the optimization process.(17) where s, K p and v start respectively represent EF, scale factor of penalty function, and vehicle speed threshold for engine start-up; SOC final is terminated SOC; f(x) denotes fitness function; γ is the weighting factor, which decides the SOC penalty term.
Based on the co-simulation platform, the three key parameters of ECMS are optimized under CHTC_HT, C_WTVC, combined driving cycle of CHTC_HT and C_WTVC.The simulation results are illustrated in Figure 15.The specific optimization results for above three driving cycle can be found in Table 5.
Energies 2020, 13, x FOR PEER REVIEW 14 of 21      The flow chart of PSO-ECMS is shown in Figure 16, which can be divided into offline part and online part.Firstly, the PSO is used to optimize EF, scale factor of penalty function, and vehicle speed threshold for engine start-up offline under specific driving cycle.Secondly, parameter optimization results are applied to real-time ECMS.Finally, optimal torque distribution is solved by ECMS.

The Simulation and Discussion of PSO-ECMS
The flow chart of PSO-ECMS is shown in Figure 16, which can be divided into offline part and online part.Firstly, the PSO is used to optimize EF, scale factor of penalty function, and vehicle speed threshold for engine start-up offline under specific driving cycle.Secondly, parameter optimization results are applied to real-time ECMS.Finally, optimal torque distribution is solved by ECMS.In order to verify the effectiveness of proposed PSO-ECMS, the simulation analysis is carried out based on co-simulation platform under the CHTC_HT, C_WTVC and combined driving cycle respectively.The initial SOC is set to 50%, and the simulation results are illustrated in Figure 17 and Figure 18. Figure 17 shows that actual vehicle speed can keep track of target vehicle speed under combined driving cycle.Figure 18 show that PSO-ECMS can reasonably allocate engine torque and motor torque under various driving cycles, as well as ensuring that initial SOC and terminated SOC are equal within the limits of the vehicle's components.In order to verify the effectiveness of proposed PSO-ECMS, the simulation analysis is carried out based on co-simulation platform under the CHTC_HT, C_WTVC and combined driving cycle respectively.The initial SOC is set to 50%, and the simulation results are illustrated in Figures 17 and 18. Figure 17 shows that actual vehicle speed can keep track of target vehicle speed under combined driving cycle.Figure 18 show that PSO-ECMS can reasonably allocate engine torque and motor torque under various driving cycles, as well as ensuring that initial SOC and terminated SOC are equal within the limits of the vehicle's components.

The Simulation and Discussion of PSO-ECMS
The flow chart of PSO-ECMS is shown in Figure 16, which can be divided into offline part and online part.Firstly, the PSO is used to optimize EF, scale factor of penalty function, and vehicle speed threshold for engine start-up offline under specific driving cycle.Secondly, parameter optimization results are applied to real-time ECMS.Finally, optimal torque distribution is solved by ECMS.In order to verify the effectiveness of proposed PSO-ECMS, the simulation analysis is carried out based on co-simulation platform under the CHTC_HT, C_WTVC and combined driving cycle respectively.The initial SOC is set to 50%, and the simulation results are illustrated in Figure 17 and Figure 18. Figure 17 shows that actual vehicle speed can keep track of target vehicle speed under combined driving cycle.Figure 18 show that PSO-ECMS can reasonably allocate engine torque and motor torque under various driving cycles, as well as ensuring that initial SOC and terminated SOC are equal within the limits of the vehicle's components.

The Principle of A Novel A-ECMS
Combing above DCR and key parameter optimization of ECMS, a novel A-ECMS is presented in this paper.The principle of proposed A-ECMS is shown in Figure 19.Firstly, the PSO is used to solve optimal parameters of EF, scale factor of penalty function, and vehicle speed threshold for engine start-up under six typical driving cycles.The results are shown in Table 6 below.Secondly, actual driving conditions will be categorized as one of six typical driving cycles by applying a driving condition recognizer based on the BP neural network, and the identification results are one of six typical driving cycle.Then the corresponding three optimal parameters of ECMS are determined according to recognized driving condition.Finally, the optimal torque distribution of hybrid heavyduty truck can be obtained by ECMS under actual driving conditions.

The Principle of A Novel A-ECMS
Combing above DCR and key parameter optimization of ECMS, a novel A-ECMS is presented in this paper.The principle of proposed A-ECMS is shown in Figure 19.Firstly, the PSO is used to solve optimal parameters of EF, scale factor of penalty function, and vehicle speed threshold for engine under six typical driving cycles.The results are shown in Table 6 below.Secondly, actual driving conditions will be categorized as one of six typical driving cycles by applying a driving condition recognizer based on the BP neural network, and the identification results are one of six typical driving cycle.Then the corresponding three optimal parameters of ECMS are determined according to recognized driving condition.Finally, the optimal torque distribution of hybrid heavy-duty truck can be obtained by ECMS under actual driving conditions.

Comparison Analysis
The PSO-ECMS combines the PSO with the original ECMS introduced in the fourth section, which will determine the key parameters of ECMS.The A-ECMS combined the key parameters optimization and the DCR introduced in the fifth section, which will adapt to the real driving condition.Based on co-simulation platform, comparison analysis is performed between PSO-ECMS and A-ECMS under composite driving cycle composed of six typical driving cycles.The three key parameters of PSO-ECMS under the composite cycle can be attained according to Figure 16, and the equivalent factor, scale factor of penalty function, and vehicle speed threshold for engine start-up are 2.001, 2.738, and 20.096, respectively.The three key parameters of A-ECMS vary with the current driving conditions which can adapt to the composite cycle, as shown in Table 6, and the specific process can be attained in Figure 19.The simulation results are as follows.
Figure 20 demonstrates that A-ECMS can keep battery SOC stable around 50% without large fluctuations compared with PSO-ECMS.Moreover, due to the lack of historical information, the same SOC are obtained at the beginning 200 s of composite driving cycle for A-ECMS and PSO-ECMS.CHTC_HT_CITY is a kind of urban driving cycle and the vehicle speed is relatively low, thus battery SOC begins to slow down with motor-only drive mode at the begin.Battery SOC begins to rise when the truck runs under CHTC_HT_SUB.During CYC_MANHATTAN, the truck controlled by PSO-ECMS tends to work in motor-only drive mode, resulting in a large drop of battery SOC.However, A-ECMS has optimal key parameters of ECMS under the current driving cycle, so the battery SOC of the truck controlled by A-ECMS can be stable near target SOC of the battery and fluctuate more smoothly.

Comparison Analysis
The PSO-ECMS combines the PSO with the original ECMS introduced in the fourth section, which will determine the key parameters of ECMS.The A-ECMS combined the key parameters optimization and the DCR introduced in the fifth section, which will adapt to the real driving condition.Based on co-simulation platform, comparison analysis is performed between PSO-ECMS and A-ECMS under composite driving cycle composed of six typical driving cycles.The three key parameters of PSO-ECMS under the composite cycle can be attained according to Figure 16, and the equivalent factor, scale factor of penalty function, and vehicle speed threshold for engine start-up are 2.001, 2.738, and 20.096, respectively.The three key parameters of A-ECMS vary with the current driving conditions which can adapt to the composite cycle, as shown in Table 6, and the specific process can be attained in Figure 19.The simulation results are as follows.
Figure 20 demonstrates that A-ECMS can keep battery SOC stable around 50% without large fluctuations compared with PSO-ECMS.Moreover, due to the lack of historical information, the same SOC are obtained at the beginning 200 s of composite driving cycle for A-ECMS and PSO-ECMS.CHTC_HT_CITY is a kind of urban driving cycle and the vehicle speed is relatively low, thus battery SOC begins to slow down with motor-only drive mode at the begin.Battery SOC begins to rise when the truck runs under CHTC_HT_SUB.During CYC_MANHATTAN, the truck controlled by PSO-ECMS tends to work in motor-only drive mode, resulting in a large drop of battery SOC.However, A-ECMS has optimal key parameters of ECMS under the current driving cycle, so the battery SOC of the truck controlled by A-ECMS can be stable near target SOC of the battery and fluctuate more smoothly.Due to variable key parameters of A-ECMS, the engine state and motor state controlled by A-ECMS are different from that controlled by PSO-ECMS.Taking the time range of 3000 s-4000 s for an example, A-ECMS tends to make the truck work in engine drive mode while PSO-ECMS tends to make the truck work in motor-only drive mode, leading to different fuel consumption of the truck as illustrated in Figure 20.The comparison of fuel consumption for hybrid heavy-duty truck under ECMS, PSO-ECMS and A-ECMS are shown in Table 7, indicating the great potential of A-ECMS in terms of energy saving.

Discussion
In order to further verify the effectiveness and adaptivity of proposed A-ECMS, a verification driving cycle composed of 'C_WTVC + CYC_UDDSHDV + CYC_WVUINTER' is presented, which has not been trained by the BP neural network.For the verification driving cycle, the recognition result is shown in Figure 21, and the recognition accuracy of driving condition is 80.1%. Figure 21 shows that a majority of driving conditions can be recognized as one of six typical driving cycles, indicating the effectiveness of DCR for unknown driving conditions.The results of PSO-ECMS and A-ECMS under the verification driving cycle are shown in Figure 22.The fuel consumption of A-ECMS is 21.92 L/100 km with an improvement of 12.9%, compared to PSO-ECMS with the fuel consumption of 25.19 L/100 km.Due to variable key parameters of A-ECMS, the engine state and motor state controlled by A-ECMS are different from that controlled by PSO-ECMS.Taking the time range of 3000 s-4000 s for an example, A-ECMS tends to make the truck work in engine drive mode while PSO-ECMS tends to make the truck work in motor-only drive mode, leading to different fuel consumption of the truck as illustrated in Figure 20.The comparison of fuel consumption for hybrid heavy-duty truck under ECMS, PSO-ECMS and A-ECMS are shown in Table 7, indicating the great potential of A-ECMS in terms of energy saving.

Discussion
In order to further verify the effectiveness and adaptivity of proposed A-ECMS, a verification driving cycle composed of 'C_WTVC + CYC_UDDSHDV + CYC_WVUINTER' is presented, which has not been trained by the BP neural network.For the verification driving cycle, the recognition result is shown in Figure 21, and the recognition accuracy of driving condition is 80.1%. Figure 21 shows that a majority of driving conditions can be recognized as one of six typical driving cycles, indicating the effectiveness of DCR for unknown driving conditions.The results of PSO-ECMS and A-ECMS under the verification driving cycle are shown in Figure 22.The fuel consumption of A-ECMS is 21.92 L/100 km with an improvement of 12.9%, compared to PSO-ECMS with the fuel consumption of 25.19 L/100 km.Due to variable key parameters of A-ECMS, the engine state and motor state controlled by A-ECMS are different from that controlled by PSO-ECMS.Taking the time range of 3000 s-4000 s for an example, A-ECMS tends to make the truck work in engine drive mode while PSO-ECMS tends to make the truck work in motor-only drive mode, leading to different fuel consumption of the truck as illustrated in Figure 20.The comparison of fuel consumption for hybrid heavy-duty truck under ECMS, PSO-ECMS and A-ECMS are shown in Table 7, indicating the great potential of A-ECMS in terms of energy saving.

Discussion
In order to further verify the effectiveness and adaptivity of proposed A-ECMS, a verification driving cycle composed of 'C_WTVC + CYC_UDDSHDV + CYC_WVUINTER' is presented, which has not been trained by the BP neural network.For the verification driving cycle, the recognition result is shown in Figure 21, and the recognition accuracy of driving condition is 80.1%. Figure 21 shows that a majority of driving conditions can be recognized as one of six typical driving cycles, indicating the effectiveness of DCR for unknown driving conditions.The results of PSO-ECMS and A-ECMS under the verification driving cycle are shown in Figure 22.The fuel consumption of A-ECMS is 21.92 L/100 km with an improvement of 12.9%, compared to PSO-ECMS with the fuel consumption of 25.19 L/100 km.

Conclusions
In this paper, a novel A-ECMS for hybrid heavy-duty truck is proposed based on DCR and key parameter optimization of ECMS.The industrial and academic contributions of this paper can be summarized as follows: (1) Six kinds of typical driving cycle are obtained by hierarchical clustering analysis, and the neural network recognizer of driving condition has been developed on this basis, simulation results indicate that real driving condition can be accurately recognized with an accuracy of 86.8%.(2) PSO is applied to optimize three key parameters of ECMS under different driving cycles, including equivalent factor, scale factor of penalty function, and vehicle speed threshold for engine start-up, and the effectiveness of PSO is verified through simulation analysis under three different driving cycles based on co-simulation platform with AMESim and Matlab/Simulink.(3) A-ECMS combining driving condition recognition and PSO has been presented.The simulation results under composite driving cycle show that proposed A-ECMS and PSO-ECMS can achieve fuel economy improvement of 14.81% and 3.05%, respectively, compared to ECMS.The next step is to conduct hardware-in-the-loop test and real vehicle test validation of proposed strategy, as well as improving control performance further.The effect of driving style and driving behavior on A-ECMS will be considered in a future study.Furthermore, the recognition accuracy of driving conditions is no more than 90%, so improving the accuracy of DCR is another difficult task for us in the future.

Conclusions
In this paper, a novel A-ECMS for hybrid heavy-duty truck is proposed based on DCR and key parameter optimization of ECMS.The industrial and academic contributions of this paper can be summarized as follows: (1) Six kinds of typical driving cycle are obtained by hierarchical clustering analysis, and the neural network recognizer of driving condition has been developed on this basis, simulation results indicate that real driving condition can be accurately recognized with an accuracy of 86.8%.(2) PSO is applied to optimize three key parameters of ECMS under different driving cycles, including equivalent factor, scale factor of penalty function, and vehicle speed threshold for engine start-up, and the effectiveness of PSO is verified through simulation analysis under three different driving cycles based on co-simulation platform with AMESim and Matlab/Simulink.(3) A-ECMS combining driving condition recognition and PSO has been presented.The simulation results under composite driving cycle show that proposed A-ECMS and PSO-ECMS can achieve fuel economy improvement of 14.81% and 3.05%, respectively, compared to ECMS.The next step is to conduct hardware-in-the-loop test and real vehicle test validation of proposed strategy, as well as improving control performance further.The effect of driving style and driving behavior on A-ECMS will be considered in a future study.Furthermore, the recognition accuracy of driving conditions is no more than 90%, so improving the accuracy of DCR is another difficult task for us in the future.

Figure 2 .
Figure 2. Fuel consumption contour map of the engine.

Figure 2 .
Figure 2. Fuel consumption contour map of the engine.
the maximum output torque of the motor at the current speed when the motor operates as a generator, and b P is the terminal battery power in kW.

Figure 5 .
Figure 5.The process of the hierarchical clustering analysis.

Figure 5 .
Figure 5.The process of the hierarchical clustering analysis.

Figure 6 .
Figure 6.Curves of six typical driving cycle.

Figure 7 .
Figure 7. On-line recognition process of driving conditions.Figure 7. On-line recognition process of driving conditions.

Figure 7 .
Figure 7. On-line recognition process of driving conditions.Figure 7. On-line recognition process of driving conditions.

21 Figure 8 .
Figure 8. Training process of the BP neural network.

CrossentropyFigure 8 .
Figure 8. Training process of the BP neural network.

Figure 9 .
Figure 9.The principle of on-line recognition.

Figure 11 .
Figure 11.The relation between prediction period and recognition accuracy.

Figure 12 .
Figure 12.The recognition result of test driving condition.

Figure 10 . 21 Figure 10 .
Figure 10.The relation between sampling period and recognition accuracy.

Figure 11 .
Figure 11.The relation between prediction period and recognition accuracy.

Figure 12 .
Figure 12.The recognition result of test driving condition.

Figure 11 .
Figure 11.The relation between prediction period and recognition accuracy.

Energies 2020 , 21 Figure 10 .
Figure 10.The relation between sampling period and recognition accuracy.

Figure 11 .
Figure 11.The relation between prediction period and recognition accuracy.

Figure 12 .
Figure 12.The recognition result of test driving condition.

Figure 12 .
Figure 12.The recognition result of test driving condition.

Figure 14 .
Figure 14.The optimization process of PSO based on co-simulation platform.

Figure 14 .
Figure 14.The optimization process of PSO based on co-simulation platform.

Figure 14 .
Figure 14.The optimization process of PSO based on co-simulation platform.
where , p K and start v respectively represent EF, scale factor of penalty function, and vehicle speed threshold for engine start-up; final SOC is terminated SOC; f(x) denotes fitness function; γ is the weighting factor, which decides the SOC penalty term.Based on the co-simulation platform, the three key parameters of ECMS are optimized under CHTC_HT, C_WTVC, combined driving cycle of CHTC_HT and C_WTVC.The simulation results are illustrated in Figure15.The specific optimization results for above three driving cycle can be found in Table5 .

Figure 15 .
Figure 15.The parameter optimization results of ECMS under various driving cycles.

Figure 15 .
Figure 15.The parameter optimization results of ECMS under various driving cycles.

Figure 17 .
Figure 17.The actual speed and target speed under combined driving cycle.

Figure 17 .
Figure 17.The actual speed and target speed under combined driving cycle.Figure 17.The actual speed and target speed under combined driving cycle.

Figure 17 .
Figure 17.The actual speed and target speed under combined driving cycle.Figure 17.The actual speed and target speed under combined driving cycle.

Figure 18 .
Figure 18.State change of main components under various driving cycle.

Figure 18 .
Figure 18.State change of main components under various driving cycle.

Figure 19 .
Figure 19.The principle of proposed A-ECMS.

Figure 19 .
Figure 19.The principle of proposed A-ECMS.

Energies 2020 , 21 Figure 20 .
Figure 20.The simulation results of PSO-ECMS and A-ECMS under composite cycle.

Figure 21 .
Figure 21.The recognition result of the verification driving cycle.

Figure 20 .
Figure 20.The simulation results of PSO-ECMS and A-ECMS under composite cycle.

Figure 20 .
Figure 20.The simulation results of PSO-ECMS and A-ECMS under composite cycle.

Figure 21 .
Figure 21.The recognition result of the verification driving cycle.

Figure 21 .
Figure 21.The recognition result of the verification driving cycle.

Figure 22 .
Figure 22.The simulation results of PSO-ECMS and A-ECMS under verification driving cycle.

Figure 22 .
Figure 22.The simulation results of PSO-ECMS and A-ECMS under verification driving cycle.

Table 2 .
Sample database of driving cycles for hybrid heavy-duty truck.

Table 3 .
Updated sample database of driving cycles for hybrid heavy-duty truck.

Table 4 .
Six types for sample data of driving cycles.

Table 4 .
Six types for sample data of driving cycles.

Table 5 .
The optimal parameters of ECMS under three driving cycles.

Table 5 .
The optimal parameters of ECMS under three driving cycles.

Table 5 .
The optimal parameters of ECMS under three driving cycles.

Table 6 .
Three kinds of optimal parameters for ECMS under six typical driving cycles.

Table 6 .
Three kinds of optimal parameters for ECMS under six typical driving cycles.

Table 7 .
Fuel consumption of hybrid heavy-duty vehicle under ECMS, PSO-ECMS and A-ECMS.

Table 7 .
Fuel consumption of hybrid heavy-duty vehicle under ECMS, PSO-ECMS and A-ECMS.

Table 7 .
Fuel consumption of hybrid heavy-duty vehicle under ECMS, PSO-ECMS and A-ECMS.