1. Introduction
The introduction of a power battery can make up for the shortcomings of fuel cell hybrid electric vehicles (FCHEV), such as the inability to recover braking energy, slow start speed and soft output characteristics. The dual power source (fuel cell and battery pack) can make the fuel cell hybrid electric vehicles (FCHEVs) produce a better power performance, but how to make the power source power distribution more reasonable and better improve the economy is a research difficulty. Based on previous experience, researchers developed rule-based energy management algorithms, such as thermostatic control strategy (TCS) [
1] and a power following control strategy (PFCS) [
2,
3]. Fuzzy control strategy (FCS) [
4,
5,
6] and fuzzy control strategy optimized by other algorithms [
7] can adapt to the requirements of vehicle nonlinear control and effectively distribute the power between the power sources of fuel cell hybrid vehicles. However, due to the lack of road condition information, they are difficult to further improve the working efficiency and the economy of power sources in complex working conditions. Another control strategy based on optimization, such as dynamic programming (DP) [
8,
9,
10], are widely used in hybrid electric vehicle energy management strategy because they can achieve global optimization. However, those methods will increase the computational burden and make it difficult to realize the online application. In order to simplify the calculation, some strategies, such as equivalent consumption minimization strategy (ECMS) [
11,
12,
13], Pontryagin minimum principle strategy (PMPS) [
14,
15] and stochastic dynamic programming (SDP) [
16], further improve the energy management performance on the basis of effectively reducing the calculation amount. For some intelligent algorithms, such as particle swarm optimization (PSO) [
17] and genetic algorithm (GA) [
18], the fuel economy can also be improved by optimizing some relevant parameters based on the rule-based control strategy.
Working conditions have a profound impact on the economy and power source performance of FCHEVs. Ahmadi et al. [
19] investigated the influence of driving patterns, and they found that various driving patterns under different conditions could affect the degradation of a fuel cell, and then affect the economy of the fuel cell vehicles. Raykin et al. [
20] investigated the influence of driving patterns under different working conditions and an electric power supply on the well-to-wheel energy use and greenhouse gases of a plug-in hybrid electric vehicle (PHEV). When formulating the FCHEVs’ energy control strategy, some references mentioned that they took single working condition into account, and there were certain limitations in improving the economy under different working conditions. Moreover, they did not consider the efficient working area of a fuel cell (FC) and battery pack to give full play to their respective advantages. Under the condition that working conditions can be identified, the energy management strategy of FCHEV should be adjusted according to the actual situation to achieve efficient and reliable power distribution among power sources, improve economy and extend the service life of power sources.
A lot of scholars have studied working condition identification. References [
21,
22,
23,
24] based on a fuzzy control recognizer, realized the identification of driving conditions. However, membership functions and rules of the fuzzy controller were selected and formulated based on personal experience, and the ideal effect could be achieved after multiple debugging. Clustering methods also play a role in the field of driving conditions recognition [
25,
26]. In [
25], working conditions were divided into five typical working conditions by way of a clustering analysis method, then working conditions were identified by a Euclid approach degree. Yu et al. [
26] identified high impact factors affecting pattern characteristics from static and quasi-static environment and traffic information, then proposed a trip/route division algorithm based on data clustering method. However, the selection of initial clustering center affected the clustering analysis results. Recently, machine learning has been further applied. Neural networks, such as back-propagation (BP) neural network [
27] and learning vector quantization (LVQ) neural network [
28,
29], involve first, characteristic parameters that have an important influence on driving conditions being selected as the input, then, the identification period of the working condition samples are classified. After training the samples, the prediction of future working conditions can be realized. However, the accuracy of neural network depends on its structure. Chen [
30] et al. proposed an improved hierarchical clustering algorithm to divide the driving cycle data into four groups, and then applied a support vector machine (SVM) to predict driving conditions based on the clustering results.
The least square support vector machines (LSSVM) based on support vector machines (SVM), compared with SVM, can complete a prediction in a shorter time and has a great generalization ability. Moreover, LSSVM is not subject to the set of algorithm structures and has good robustness in handling regression and classification problems.
In order to improve the performance of FCHEV, this paper proposes a driving condition recognizer. By extracting feature parameters and segmenting recognition segments from driving conditions information, LSSVM optimized by CV is used to realize working condition recognition. Energy management controllers based on a fuzzy control under different working conditions are established and optimized. Combined with the driving conditions identification, the energy management controller adopts corresponding fuzzy control strategy according to driving conditions to improve the performance of FCHEV.
2. Vehicle Structure and Parameters
The FCHEV was a front-drive vehicle with the structure shown in
Figure 1. The fuel cell system was connected to the Controller Area Network (CAN) bus through a one-way DC/DC converter, while the battery pack was directly connected to the CAN bus. The motor drives the vehicle through the final drive and differential. The complete vehicle parameters of a fuel cell hybrid electric vehicle are shown in
Table 1.
In this paper, the vehicle model of FCHEV was established in AVL Cruise, as shown in
Figure 2, and the control strategy model was established in Matlab/Simulink, shown in
Figure 3. In
Figure 2, the overall simulation model includes driver module, fuel cell system, power battery pack, motor and controller, one-way DC/DC converter, final drive, and energy management module. The blue line and red line represent mechanical connection electrical connection, respectively.
2.1. Fuel Cell Module
The fuel cells in this paper were proton exchange membrane fuel cells (PEMFC), and they were built out of membrane electrode assemblies (MEA), which included the electrodes, electrolyte, anode catalyst layer, cathode catalyst layer (CCL), and gas diffusion layer (GDL). The detailed modeling process is found in references [
31,
32]. In the fuel cell component, in addition to the fuel cell, there was a simple compressor model, and its properties are shown in
Table 2. The compressor delivered hydrogen continuously to the fuel cell stack, which generated electricity to drive the motor.
The voltage of the fuel cell electrochemical model is calculated as follows:
where
Ufc is the output voltage,
Uoc is the ideal open circuit voltage,
η0 is the cathode voltage loss,
Vact is the activation over potential,
VCCL is the voltage loss caused by the oxygen transmission loss in the cathode catalyst layer (CCL),
VGDL is the voltage loss caused by the oxygen transmission loss in the anode catalyst layer,
j0 and
Ist are the electric flow density and current of the stack, while
Aarea is the effective area of the fuel cell,
R is the ohmic internal resistance of the fuel cell. The activation loss can be defined as follows.
where
bTf is the Tafel slope which describes the speed of the chemical reaction, and
ccc is the oxygen concentration in the cathode channel, while
cci is the oxygen concentration at the channel inlet. Moreover,
ja and
j* can be defined as
where
i* is the volumetric exchange current density, and
Spc is the CCL proton conductivity, in addition,
lCCL is the thickness of the CCL.
The voltage loss
VCCL can be defined as
where
F is the Faraday constant,
DCCL is the oxygen diffusion coefficient in the CCL.
j*
l and
B can be defined as
where
DGDL is the oxygen diffusion coefficient in the GDL, while
lGDL is the thickness of GDL.
The voltage loss
VGDL can be defined as
Assuming that the fuel cell stack consists of
n fuel cell cells, the output power of the fuel cell stack is
The efficiency of fuel cell stack can be expressed as follows:
The single fuel cell properties are shown in
Table 3.
2.2. Power Battery Pack
The lithium battery selected in this paper had a capacity of 24 Ah and a rated voltage of 3.3 V, and its specific parameters are shown in
Table 4. Its equivalent circuit model adopted the Rint model, as shown in
Figure 4a. The voltage of the battery output to the CAN bus is:
where
UOCV is the open circuit voltage of lithium battery,
UOut is the output voltage,
Ib and
R0 are the current and ohmic internal resistance of lithium battery respectively.
SOC, an important parameter of a lithium battery, is expressed by the following equation:
where
η is the coulomb efficiency, in this paper,
η = 1, SOC
0 was the initial value, sampling time Δ
t = 1 s, and
CP was the actual capacity of the battery. Through experiments, the parameters relationships of the battery are shown in
Figure 4b.
3. Typical Driving Conditions
Working conditions of a vehicle have an important impact on economy and power distribution. Therefore, a more efficient energy management strategy can be developed by predicting the future working conditions.
In this paper, three typical driving conditions were selected, namely UDDS (Urban Dynamometer Driving Schedule), EUDC (Extra Urban Driving Cycle) and US06 (Highway Driving Schedule), as shown in
Figure 5, which corresponded to an urban condition, suburban condition and highway condition, respectively. In an urban working condition, the vehicle speed is low and frequent parking occurs. The average vehicle speed is less than 35 km·h
−1, moreover, the vehicle is in a state of low power output. The speed is fast in highway conditions, and the average speed is about 70 km·h
−1, in addition, the output power of the car is relatively large. The suburban working condition is in the middle of the two, with an average speed of about 60 km·h
−1.
3.1. Selection of Working Condition Characteristic Parameters
The selection of characteristic parameters of working conditions is the key to accurately identifying future working conditions. In principle, more characteristic parameters is more helpful for prediction, but that requires high computational power. In contrast, too few characteristic parameters cannot cover the information of working conditions, which may lead to a large prediction deviation. Many scholars have studied the selection of characteristic parameters of driving conditions [
25,
30,
33,
34,
35]. Based on some research and the importance of each parameter in driving conditions identification, six common characteristic parameters were selected, i.e., acceleration time/total time (
rc), deceleration time/total time (
rdc), time of uniform speed/total time (
ru), average speed (
va), average acceleration (
ac) and average deceleration (
adc). Six characteristic parameters of three working conditions are shown in
Table 5.
3.2. Dividing of Working Condition Samples
The time length of the working conditions samples, namely, the identification cycle and update of identification cycle, will also have an impact on the working condition recognition. The specific segmentation of the working condition recognition samples is shown in
Figure 6. △
T is the identification period, therefore, six characteristic parameters in this period of time can be calculated to identify the working conditions of this sample. While △
s is the update of the period, that is, the time difference between the beginning of the previous cycle segmentation and the beginning of the current cycle segmentation. If △
T is too long, although it contains more information, it will increase useless information and calculation burden, which will reduce the effect of recognition. If △
T is too short, it will not accurately reflect the real situation of working conditions. Similarly, a too small △
s leads to frequent cycle switching, which will cause a burden on the processor, while a too large △
s is not conducive to the timely switching of working conditions. References [
35,
36] studied in detail the effect of △
T and △
s on the accuracy of working conditions identification. Based on considering the accuracy and calculation cost, △
T = 100 s, △
s = 3 s.
4. Working Condition Identification Model Based on LSSVM
4.1. Least Squares Support Vector Machine
LSSVM is able to classify samples by mapping them into high-latitude feature Spaces. LSSVM replaces the inequality constraints of problems in SVM with a set of linear equality constraints, thus simplifying the solution of Lagrange multipliers. A training set is considered with
n data samples to be (
Xi,
yi), where input data
Xi ∈ R
n, output data
yi ∈ R. A linear function in the high-level feature space will be used to fit the samples.
where
φ(
X) is a nonlinear mapping function,
ω is the weight vector in the feature space, and
b is the bias term.
According to the principle of structural risk minimization and taking into account the complexity of function and fitting error, the optimization problem of LSSVM can be expressed as:
where
ξi is the error variable and
C is the penalty factor.
Converting Equation (16) to unconstrained functions by building Lagrange functions and solving this Lagrange function, the classification prediction model of LSSVM can be obtained, as shown in Equation (8), and its structure is shown in
Figure 7. Combined with
Section 3, six characteristic parameters are taken as the input of the LSSVM, and the output is the working condition categories:
where the radial basis function (RBF) is selected as the kernel function, namely
,
α is the Lagrange multiplier.
4.2. The Influence of Key Parameters on the Accuracy of LSSVM
If σ→0, then K(X, Xi)→0, which means that all the mapped points have the same distance from each other, that is, there is no clustering phenomenon. However, If σ→∞, then K(X, Xi)→1, which means that all sample points will be divided into the same class and cannot be distinguished. As for the penalty factor C, if C is too large, ξi→0, the tolerance of samples between boundaries is very low, and there are less misclassifications, which means the fitting of samples is good, however, the prediction effect is not always good; on the other hand, if the value of C is too small, there are more samples between two boundaries, resulting in a greater possibility of misclassification, and the fitting of samples decreases.
The accuracy of LSSVM’s model depends on the kernel parameter σ and the penalty factor C. A too large σ will reduce the model’s accuracy, but a too small σ will lead to overfitting. The penalty factor C will affect the error and complexity of the model. Therefore, in this study, the cross-validation method was used to obtain the optimal parameters.
4.3. The K-Fold Cross-Validation for Optimizing LSSVM
Cross-validation has been widely used to estimate prediction errors. In this work,
K-fold cross-validation combined with grid search was applied to optimize LSSVM, which could overcome the limitations of the holdout validation [
37]. The steps to optimize LSSVM were as follows:
- (1)
Establish grid coordinates. Let a = [–10, 10], b = [–10, 10], and the step size is 0.5, then the mesh points of the model parameters are σ = 2a and C = 2b respectively. In this work, the exponential function was selected to divide the grid, which would ensure that the parameter value was not negative.
- (2)
Divide the sample data and calculate the test error. The training data are divided into K subsets (K = 10, which means that the CV is 10-fold cross-validation method). For each group (σ, C) in the grid, a 10-fold cross validation method was applied to iterate the training data 10 times, and the mean value of the mean square error (MSE) of the test results under this group of parameters could be obtained.
- (3)
Get the optimal combination of parameters. Repeat (2) to replace the parameter σ and C, and calculate the mean square deviation of the training model under all the parameter combinations in the grid in turn. After comparing one by one, the parameter combination corresponding to the minimum mean square deviation is the optimal parameter combination in the grid interval.
In order to present equidistant grid search results more clearly, grid coordinates (σ, C) are converted to logarithmic coordinates (log2σ, log2C).
5. Fuzzy Energy Management Strategy Based on Working Condition Identification
Fuzzy control based on the theory of fuzzy mathematics, fuzzes the actual input and output, and formulates rules through experience. These kinds of simulation of a human’s approximate reasoning and comprehensive decision-making process has good robustness and adaptability. Fuzzy energy management strategies [
4,
5,
6] developed by some researchers were aimed at a single working condition. In addition, fuzzy control rules based on personal experience are difficult to deal with complex multi-working conditions. Therefore, on the basis of condition identification, three fuzzy energy management strategies were formulated to deal with urban, suburban and expressway conditions, respectively. Besides, PSO is used to optimize the fuzzy control under various working conditions with total equivalent energy consumption as the objective function, and the adaptive switching effect is achieved through the identification of working conditions. It should be noted that the following fuzzy controller and optimization take the urban working condition as an example.
5.1. Fuzzy Controller Design
- (1)
Selection of input and output variables of fuzzy controller.
The SOC of the battery pack and the total power demand
Pr of FCHEV were selected as the input of the fuzzy controller, while the output is the output power
Pfc of the fuel cell. The power demand relationship is as follows:
where
Pb is the output power of the battery, and
Pr includes the power of the drive motor and the power consumed by accessories.
- (2)
Fuzzy distribution of input and output variables.
The range of FCHEV’s total power demand Pr is [0, 60] (kw), and its fuzzy subsets are very small, small, medium, large and very large, i.e., {VS, S, M, L, VL}; the SOC range of power battery is [0, 1], and the fuzzy subsets are {VL, L, M, H, VH}, representing very low, low, medium, high and very high; the range of fuel cell’s output power Pfc is [0, 50] (kw), hence its fuzzy subsets {VL, L, M, H, VH} represent very low, low, medium, high and very high.
- (3)
Fuzzy control rules.
The fuzzy control rules of FCHEV are formulated according to the following principles:
- ①
When the SOC of the power battery is too low, the output power of the fuel cell should not only meet the requirements of driving the vehicle, but also charge the battery to make the SOC of the power battery rise to a reasonable range (SOC = 40–80%).
- ②
When the SOC of the power battery and the demand power are both medium level, the fuel cell acts as the active power source and changes with the demand power. SOC of battery fluctuates in a reasonable range, which is beneficial to prolonging battery life.
- ③
When SOC is too high, the power battery acts as the main power source, and the output power of the fuel cell is as small as possible to reduce the SOC to a reasonable range.
- ④
When the demand power is too large, the power battery and fuel cell provide output power together.
5.2. Fuzzy Controller Optimization Based on PSO
As an optimization algorithm, PSO is a solution to reducing the influence of making fuzzy control strategy based on personal experience. In this paper, the membership functions of the input and output of the fuzzy controller were selected as the parameters to be optimized, and the objective function was total equivalent energy consumption (TEEC) of the power sources, i.e.,
where
Efc(
x) and
Eb(
x) are the equivalent electric energy consumption of the fuel cell and electric energy consumption of battery, respectively, while
Gi(
x) is the constraint condition of the vehicle, such as the time of acceleration and SOC fluctuation range of the battery pack.
The distributions of control rules under urban working condition before and after optimization are shown in
Figure 8.
5.3. Fuzzy Energy Management Based on Condition Identification
After identifying working conditions by LSSVM, the corresponding fuzzy control rules are selected by the fuzzy controller according to the working conditions. The flow chart of the energy management strategy based on working conditions identification is shown in
Figure 9. Firstly, the characteristic parameters were extracted from the working condition information and sample segmentations were determined, and then working condition identifications were carried out by LSSVM. Fuzzy control strategies were optimized under three working conditions, and corresponding fuzzy control rule was selected under a specific working condition to realize the adaptive switching of the control strategy under complex working conditions.