Control Strategy for PHEB Based on Actual Driving Cycle with Driving Style Characteristic

. To exert fully the energy economy performance of plug-in hybrid electric buses (PHEBs) and enhance the adaptability to different drivers and driving cycles, a control strategy for PHEB based on actual driving cycle with driving style characteristicis proposed in this paper. Through the actual city bus driving data, collected in real time, 6 actual driving cycles with driving style characteristic are fitted by using Principal Component Analysis (PCA) and Cluster Analysis (CA). Based on the 6 driving cycles, the key parameters of rule-based control strategy are optimized and established by a combinatorial optimization algorithm in Isight. Then, an identification model to recognize the current condition based on the Learning Vector Quantization (LVQ) neural network has been built and trained offline, which is integrated in the control strategy for PHEB to invoke the corresponding optimized key control parameters in real time. A hardware-in-the-loop (HIL) test is conducted, and the result shows that the proposed strategy could improve the energy consumption by 4.94%, compared with the original rule-based control strategy, and its validity and practicability are fully verified.


Introduction
Internal combustion engine (ICE) and electric machine (EM) are collectively used as the power sources of PHEB by electromechanical coupling system, which makes PHEB acquire more superior economy and emission performance compared with hybrid electric bus (HEB), and wider driving range and better environment adaptability compared with electric bus (EB).But, on account of the diversity of power sources, the complexity of electromechanical coupling system, and the variety of road driving cycles, there are some challenges for PHEB energy meticulous management [1,2].
At present, the PHEB strategies mainly include rulebased and optimization-based strategies, where the former include the strategies based on deterministic rules and fuzzy rules, and the latter include global optimization strategies based on dynamic programming (DP), etc., and real-time optimization strategies based on model predictive control (MPC), etc. [3,4].In [5], the rule-based strategy was used to allocate the output between the hybrid energy storage system (HESS) and the assistance power unit (APU), and the MPC was introduced to regulate the output between the battery and the ultracapacitor.Reference [6] predicted the electric power demand by the state transition probability matrices of the electrical power with different speeds, and the result was that the overall loss incurred by the whole HESS is reduced.Although the above control strategies have obtained some certain positive results, they rely on the standard driving cycles too much, resulting in a certain deviation with actual driving scenes.What is more, there are some studies which focus on the drivers' driving behaviors when developing the strategies.For example, the demanding torque correction factor was adjusted to meet different driving intentions in [7].In [8], the driver behaviors were regarded as stochastic disturbance inputs, and the stochastic MPC was modified with the equivalent consumption minimization strategy (ECMS).However, these control strategies do not Collecting Data take the driving cycles and driving styles into simultaneous consideration, and the economy of PHEB can be further improved.So it is of great significance to fit the actual driving cycles with driving style characteristic and then optimize the control strategy for PHEB to exert fully its comprehensive performance.
In this paper, a control strategy for PHEB based on actual driving cycle with driving style characteristic is proposed, as shown in Figure 1, mainly including the progress to fit out actual driving cycles with driving style characteristic and the progress to optimize the control strategy based on the above driving cycles.The structure of this paper is as follows.A large amount of city bus road driving data is collected in real time and divided into different types of driving cycles, and driving styles through different characteristic parameters by using PCA and CA, and then 6 actual driving cycles with driving style characteristic are fitted finally in Section 2. In Section 3, the key parameters of the control strategy are optimized by Isight automatic optimization platform by using a combinatorial optimization algorithm with the above actual driving cycles, and a control strategy for PHEB is built to invoke the corresponding optimized key control parameters in real time with a driving cycle identification model, which has been trained through PCA and LVQ neural network.A HIL test is designed to verify the validity and practicability of the control strategy for PHEB based on actual driving cycle with driving style characteristic in Section 4. Finally, Section 5 concludes the paper.

Actual Driving Cycles with Driving Style Characteristic
Combined with the traffic characteristics of city bus, the kinematic fragment analysis is usually served as the driving cycle development method and then through PCA and CA, and then some kinematics fragments with driving cycle and driving style characteristics will be selected to fit out the actual driving cycles with driving style characteristic [9,10].

Data Acquisition and Preprocessing.
The city bus driving data are acquired through China Automotive Test Cycles (CATC) information system, and the schematic of data acquisition is expounded in Figure 2.
Taking into consideration that the bus lines are fixed and the buses always run on different lines with different speeds, the bus lines can reflect the city bus driving characteristic to some extent, and the chosen bus lines are expounded in Figure 3(a), where the 21 bus lines covered the main road, the subtrunk road, the branch road, the loop line, and the outside area.And the bus type and the equipment, used in the data acquisition progress, are shown in Figures 3(b) and 3(c), respectively.
Through the vehicle terminal equipment, the city bus operating data is collected in real time.The preprocessing progress is the kinematic fragment analysis progress; that is, several characteristic parameters are extracted to represent the kinematic fragments.In this paper, a total of 20182 kinematic segments are used, and there are 12 characteristic parameters, as shown in Table 1, which are introduced to describe the characteristics of each kinematics fragment accurately and comprehensively.So the total actual driving  sample data can be represented by a 20182 (sample size)×12 (eigenvalue) matrix, denoted as X, as shown in Table 2.

PCA and CA.
Although the 12 characteristic parameters can describe the kinematics fragment, there are some considerable differences in the influence extent.Therefore, the PCA is introduced to analyze influence extent of each characteristic parameter to simplify the following CA.Before PCA, to eliminate the different dimensions influence, the above matrix X needs to be standardized firstly.The matrix, after being standardized, is denoted as Y, and the calculation method is where x ij is the jth characteristic parameter in ith kinematics fragment,   is the mean of jth characteristic parameter, sd i is the standard deviation of jth characteristic parameter, and n is the total number of kinematics fragments.Then the PCA is conducted according to [11,12]; if the cumulative contribution rate of the top several principal components can get to 80%, it can be taken for granted that the top several principal components are retained.The result of PCA is shown in Table 3, and the cumulative contribution rate of the top 4 principal components has reached 83.767%, which can be used to characterize the kinematic fragment.And the principal component score of each kinematic fragment, which can characterize the pertinence between principal components and kinematic fragment, is illustrated in Table 4.
To achieve ideal clustering effect, the optimal clustering number needs to be determined primarily.The optimal clustering number is determined by Calinski-Harabasz (CH) index [13], which is as follows: where n is the clustering number, c is the category of current driving cycle, trD(k) is the trace of scatter matrix within the same driving cycle class, and trA(k) is the trace of scatter matrix between different driving cycle classes.From Equation (2), a larger CH index represents a closer relation within the same class and a more dispersed relation between different classes.The CH index result is shown in Figure 4; the optimal clustering number is 3, which are divided as congested, urban, and suburban, and the cluster centers are shown in Table 5.The driving style characteristic is embodied in the fact that the energy consumption is often quite different under different driving styles when the same vehicles are in the same route and similar driving conditions.The average acceleration and the acceleration standard deviations are introduced to represent driving style characteristic, which can embody the driver's demand for power and the dispersion degree of the acceleration in the driving process [7].Consider where n is the sampling frequency, a i is the acceleration value at time i corresponding to the fragment, and a is the acceleration sampling value. 2 kinds of driving styles, economic and dynamic, are clustered by combining  and sd(a) after the above driving cycle's CA, and the cluster centers are shown in Table 6.

Driving Cycle with Driving Style Characteristic.
According to the Constant Volume Sampling (CVS) and the experience of fitting driving cycle, 6 driving cycles with driving style characteristic, of which the duration is about 1200s, are fitted out [14].When selecting kinematics fragment, the Euclidean distance between the kinematics fragment and the cluster center is introduced; these kinematics fragments will be selected, which have the minimum sum of driving cycle Euclidean distance ( - ) and driving cycle Euclidean distance ( - ) between the Ith kinematics fragment and the cluster center of the Jth column that the kinematics fragment belongs to, where S I,N is the Nth principal component score of the Ith kinematics fragment,    is the Nth principal component score of the driving cycle cluster center, N=4 in this paper, a I and sd(a) I are the average acceleration and the acceleration standard deviations of the Ith kinematics fragment, and a J and sd(a) J are the driving style cluster center.What is more, in order to facilitate the vehicle to reproduce and follow the driving cycles well, it is necessary to smooth driving cycles that have been fitted.And Figure 5 shows the different kinds of city bus driving cycles under different driving styles after smooth filtering.
Considering that a certain identification time period in the following control strategy development is needed, it is where t i is the duration of the cluster i in comprehensive driving cycle, T overall is the duration of all kinematics fragments, T i is the duration of all kinematics fragments in cluster i, and t dc is the duration of the comprehensive driving cycle, about 1200s.The comprehensive driving cycle is shown in Figure 6, whose duration is 1174s, including 273s congested, 413s urban, and 488s suburban.

PHEB Control Strategy Based on Driving Cycle with Driving Style Characteristic
As is shown in Figure 7, the control strategy based on driving cycle with driving style characteristic is to add an identification module to recognize the current run condition and invoke the corresponding optimal control parameters to realize the adaptive effect of running in different conditions and driving styles.A particular point that needs to be explained is that the control strategy cannot identify the current condition in the first period; it will call the optimal control parameters of comprehensive driving cycle.

The Rule-Based Control Strategy.
A 12m single axis parallel-series PHEB is selected as the research object; the PHEB system configuration and main technical parameters are shown in Figure 8 and Table 7, respectively.The rule-based control strategy is designed for the PHEB system, which is shown in Table 8.The basic operation mode is divided into pure electric driving (PED), hybrid driving charge depleting (HDCD), and hybrid driving charge sustaining (HDCS) based on the battery SOC profile [15].

LVQ Neural Network Identification Model.
The identification model is built to identify the driving cycles with driving style characteristic during the driving process.Although it is crucial to use more characteristic parameters for accurate identification, its complexity will be increased because these parameters are related to each other and have different dimensions.So an identification method is introduced by combining PCA with LVQ neural network according to [16]; the flow chart is shown in Figure 9.
Considering the impact of the identification period on identification effect, the identification accuracy has been nearly 90% within 120s period in the actual driving condition, as shown in Figure 10, so that 120s is selected as the identification period.
In order to prove the effectiveness and universality of the identification model, a comprehensive operating condition of the reacquisition of one month's bus driving condition data is constructed as the testing driving cycle, which lasts 1324s; the recognition result is shown in Figure 11.
From Figure 11, '1, ' '2, ' '3, ' '4, ' '5, ' and '6' represent congested-economic, congested-dynamic, urban-economic, urban-dynamic, suburban-economic, and suburbandynamic, respectively.In the first 120s period, the driving cycle series is identified as '0'; that is to say, there are no identification result, because the period cannot satisfy the dual identification model.The rest period identification results are in line with the trend of velocity change and can reflect the driving cycle and driving style relatively more accurately and stably.

Key Control
where the F(x) is the objective function of the optimization question, F fuel , F electrical are, respectively, fuel and electrical consumption per 100km,   (x)≥0 is the dynamic constraint, n is the number of constraints, x j are the key control strategy parameters, x i j , x h j are, respectively, the range value set of each variable, and m is the number of variable.The key control parameters and the constraints are introduced in Table 9.
The PHEB model built in AVL-Cruise and the control strategy built in MATLAB/Simulink are integrated in the Isight automatic optimization platform, illustrated in Figure 12.The key control parameters are optimized by using a compositional optimization algorithm, Multi-island Genetic Algorithm (MIGA), and Sequential Quadratic Programming (SQP) Algorithm.
The optimal key control parameters values of the above driving cycles with driving style characteristic are shown in Table 10.From Table 10, there is a large distinction in different driving cycles, which will improve the vehicle's robustness and adaptability to complicated work conditions.

Hardware-in-the-Loop Test
To verify the validity and practicability of the control strategy proposed in this paper, a HIL test is introduced, as is shown in Figure 13.The HIL test platform consists of a hardware system and a software system, where the hardware system  Input Layer Competitive Layer Linear Output Layer      full play, and the battery SOC changes down more slowly.That is, the control strategy proposed in this paper can function better than the electric energy utilization efficiency of the hybrid power system.As shown in Figure 16, (a) (c) are, respectively, the torque distribution diagram of the control strategy proposed in this paper and the rule-based control strategy, and (b) (d) are, respectively, the partially enlarged diagram from 340s to 440s of (a) (c).From (b) (d), it can be seen that the working time and frequency of ICE are significantly reduced using the control strategy proposed in this paper, which makes the EM more involved to drive the whole vehicle, when it is in the inefficient working area of the ICE, and thus reduces the comprehensive energy consumption.
In the test progress, the fuel and electricity consumption are obtained by the integral at different moments and converted into the consumption of 100 km.In addition, the electricity consumption is converted into fuel consumption by Equation (6), and the sum of fuel and electricity is as the comprehensive energy consumption.The result is shown in Table 11; the proposed control strategy could improve the energy consumption by 4.94% compared with the original rule-based control strategy, and its validity and practicability are verified.

Conclusions
Using the city bus data collected in real time, through PCA and CA methods, the optimal number of driving cycles to cluster is determined by CH index, and the 6 typical bus driving cycles with the driving style characteristic.Through the Isight automatic optimization platform, the key control parameters are optimized by taking advantage of the 6 driving cycles with driving style characteristic.And an identification model for driving cycle with driving style characteristic is established to invoke the corresponding optimized key control parameters in real time.
Through the HIL test, the energy economy of the control strategy in this paper is improved by 4.94%, compared with the rule-based control strategy, which proves that the control strategy can better improve the torque distribution between ICE and EM.
) = F fuel + F electrical /3.02 s.t.g i (x) ≥ 0, i=1, 2, 3. ..n x l j ≤ x j ≤ x ℎ j , j = 1, 2, 3. ..m Col Col Col Col Col Col Col Col Col Col Col Col Col Col Col Col Col Col Col Col Collec lec lec lec lec lec lec lec lec lec lec lec lec lec lec lec lec lec lec lec lectin tin tin tin tin tin tin tin tin tin tin tin tin tin tin tin tin tin tin tin ting D g ata ata ata ata ata ata ata a a a ata ata ata ata ata ata ata ata ata Dynamic constraint e acceleration time 0∼18km/h ≤ 5s e acceleration time 0∼50km/h ≤ 20s e absolute value of Speed tracking error ≤3km/h e maximum climbing degree (15km/h)≥

Figure 1 :
Figure 1: Structure schematic of this research.

Figure 2 :
Figure 2: Schematic diagram of data acquisition process.
(a) Bus lines for data acquisition (b) A bus for data acquisition OBD interface GPS antenna Working indicator light GPRS antenna (c) Data acquisition equipment

Figure 3 :
Figure 3: Routes, bus, and equipment for data acquisition.

Figure 7 :Figure 8 :
Figure 7: Control strategy based on dual identification.

Figure 11 :
Figure 11: Recognition result of dual identification.

Figure 16 :
Figure 16: Demand torque and ICE torque of two strategies.

Table 3 :
Contribution rates of each principal component.

Table 4 :
Principal component score matrix.

Table 5 :
Principal component cluster centers.

Table 6 :
Driving style cluster centers in 3 driving cycles.

Table 7 :
Basic parameters of whole vehicle.

Table 8 :
The whole vehicle torque distribution method.   are the ICE minimum and maximum output torque at the current speed, respectively;   ,   and   are the torque of ICE, EM, and ISG at the current speed, respectively.
is the battery SOC limit value when entering HDCD mode;   is the battery SOC limit value when HDCS mode;   is the supercapacitor SOC; MC is the state of mode clutch (0 means the clutch is broken, 1 means clutch is engaged);   is the speed of the motor at the current moment;   is the speed limit value with clutch engagement;   is the required torque at the current moment;    ,

Table 9 :
Optimization parameters and conditions.
where   0∼18 is the acceleration time (0∼18km/h);   0∼50 is the acceleration time (0∼50km/h);  V is the absolute value of speed tracking error;  15 is the maximum gradeability (15km/h); e h, e l are the correction coefficients of ICE working upper limit and lower limit, respectively; disc is the discharge coefficient of the supercapacitor; cs c is the joint speed of mode clutch at CS stage.

Table 10 :
The change of vehicle control strategy parameters before and after the optimization of various conditions.

Table 11 :
Energy comparison of hardware in the loop test.