Control Strategy Optimization for Parallel Hybrid Electric Vehicles Using a Memetic Algorithm

Hybrid electric vehicle (HEV) control strategy is a management approach for generating, using, and saving energy. Therefore, the optimal control strategy is the sticking point to effectively manage hybrid electric vehicles. In order to realize the optimal control strategy, we use a robust evolutionary computation method called a “memetic algorithm (MA)” to optimize the control parameters in parallel HEVs. The “local search” mechanism implemented in the MA greatly enhances its search capabilities. In the implementation of the method, the fitness function combines with the ADvanced VehIcle SimulatOR (ADVISOR) and is set up according to an electric assist control strategy (EACS) to minimize the fuel consumption (FC) and emissions (HC, CO, and NOx) of the vehicle engine. At the same time, driving performance requirements are also considered in the method. Four different driving cycles, the new European driving cycle (NEDC), Federal Test Procedure (FTP), Economic Commission for Europe + Extra-Urban driving cycle (ECE + EUDC), and urban dynamometer driving schedule (UDDS) are carried out using the proposed method to find their respectively optimal control parameters. The results show that the proposed method effectively helps to reduce fuel consumption and emissions, as well as guarantee vehicle performance.


Introduction
The development of clean vehicles with high fuel economy and low emissions is gradually becoming mainstream in the automotive industry owing to the aggravation of the global energy crisis and environmental problems.In the field of vehicle engineering, conventional power systems driven by internal combustion engines (ICEs) have several disadvantages that adversely affect fuel economy and emissions.Furthermore, ICEs are generally over-designed approximately 10 times to meet the required vehicle driving performance that causes the cruising operating point to deviate away from the optimal operation point [1].Hybrid electric vehicles (HEVs) do not certainly require external battery charging and new infrastructure; therefore, many researchers have focused on HEV in the past few years.Additionally, their superior fuel economy and lower emissions with no compromise in dynamic performance make HEVs a viable solution for providing cleaner and more fuel-efficient vehicles.Tanoue et al. [2] pointed out that hybrid technology is of crucial importance for future automobiles.
An HEV includes at least two energy converters, i.e., ICE and an electric motor (EM), to generate the mechanical energy required to drive the vehicle and operate the on-board accessories.Therefore, energy flow management plays an important role in HEV efficiency.In order to make HEVs as efficient as possible, the HEV control strategy [3][4][5][6][7][8][9][10][11] is used to properly manage their energy components.By designing an appropriate control strategy, the HEV not only cuts down toxic exhaust emissions, but also maintains the performance of the road-driven vehicle while minimizing fuel consumption.A parameter-optimized control strategy can be used for solve this problem.
In early studies, rule-based control was widely used for parametric optimization.Dynamic programming (DP) can be used to decide the global optimal control parameters [12] when the driving cycle and vehicle performance are known.Pontryagin's minimum principle (PMP) may be used to determine the local optimal control parameters, and it can obtain an optimal trajectory in less computing time than DP [13].Many HEVs have been discussed in literature, Assanis et al. attempted to improve the fuel consumption (FC) while maintaining the driving performance within the standard limits and using the determined optimal sizes of the ICE, EM and battery pack [14].Sciarretta et al. proposed equivalent consumption minimization strategy (ECMS) a strategy based on the definition of the fuel equivalent of the electrical energy for the real-time load control of parallel HEVs [15].Johri and Filipi [16] developed a supervisory controller for HEVs based on the principles of reinforcement learning and neuro-dynamic programming to overcome the curse of dimensionality, and solving policy optimization for a system with very large design state space.Poursamad and Montazeri [17] presented a fuzzy logic controller that is tuned by a genetic algorithm for parallel HEVs to minimize fuel consumption and emissions.Panday and Bansal developed a fuel efficient energy management strategy which for power-split hybrid electric vehicle using modified state of charge estimation method.Their energy management strategy adapts GA to compute the optimal values of various governing parameters, and then uses Pontryagin's minimum principle to decide the threshold power at which engine is turned on [18].Several targets are usually considered simultaneously in HEV control strategy.To minimize vehicle fuel consumption and engine emissions, and maintain driving performance is the primary purpose of an HEV control strategy.Control strategies designed based on engineering intuition usually fail to achieve satisfactory overall system efficiency owing to the complex nature of HEVs.Several studies based on evolutionary computation have been proposed to determine the optimal parameters for the control strategy for HEVs.Montazeri and Poursamad [19] proposed a genetic algorithm (GA) to minimize the weighted sum of FC and emissions, and they also considered the Partnership for a New Generation of Vehicles (PNGV) performance requirements as constraints [20].Wu et al. [21] employed a particle swarm optimization (PSO) method to determine the optimal parameters of the powertrain and the control strategy to reduce the FC, emissions, and manufacturing costs of HEVs.Long and Nhan [22] used bees algorithm (BA) to minimize the weighted sum of FC and emissions, and considered the PNGV constraints for the vehicle performance.Hao et al. [23] employed the DIvided RECTangle (DIRECT) algorithm to optimize the extracted seven key parameters based on the optimization model of the key parameters from the perspective of fuel economy.More other control algorithms for HEVs can be referred to [24,25].
The characteristics of the powertrain system are highly non-linear and discontinuous, and may contain several local optima.However, many methods proposed in the literature lack local search capabilities that makes the optimal solution is difficult to obtain.In this study, a robust evolutionary computation method-memetic algorithm (MA) [26] is proposed to optimize the control parameters.The "local search" mechanism implemented in the MA greatly enhances its search capabilities.Therefore, it is particularly suitable for solving such problems.MAs were inspired by Dawkins' notion of a meme [27].They are similar, yet superior to the GA.MAs progress through a local search before becoming involved in the evolution process [28] to ensure that all chromosomes and offspring gain some experiences.To solve this problem, a single objective problem is used and a goal-attainment method is substituted for the original multi-objective optimization problem.The goal of optimization in this problem is to minimize the engine fuel consumption and emissions within the given criteria.In addition, vehicle performance requirements must also be maintained.The used fitness function is proposed according to an electric assist control strategy (EACS) [4].The proposed method is then performed for four different driving cycles, including new European driving cycle (NEDC), Federal Using MAs to solve the control strategy optimization problem associated with minimizing engine fuel consumption and emissions while maintaining driving performance saves a considerable amount of time and provides a global optimal solution.The MA proposed for the control strategy optimization correctly and quickly identifies the optimal parameters required for parallel HEVs.Our experimental results further indicate that MAs applied in the design of control strategy are effective in minimizing minimize fuel consumption and emissions while maintaining driving performance.

Configurations of HEVs
Before discussing the configurations of HEVs, we first mention the differences between a conventional HEV and a plug-in HEV.There are two fundamental differences between them, one is the main power source and the other is the whole energy efficiency of the two vehicle architectures.The main power source is a gasoline-powered ICE for conventional HEVs, while it is a battery-powered electric motor for plug-in HEVs.In conventional HEVs, the electric motor is used to complement the ICE and electricity is generated on board.Its energy savings are not substantial in comparison to plug-in HEVs.In plug-in HEVs, the ICE is used to complement the electric motor and grid-supplied electricity.Its energy savings are more substantial than conventional HEVs.Other types of prominent technology [29] are also used.Series and parallel type are currently two common configurations for HEVs.Another common type of HEV synthesizes the characteristics of both series and parallel types and is called dual-mode or multi-mode type [30].

Series HEV
The configuration of a series type HEV mainly consists of a fuel converter/ICE, a generator, an energy storage system/battery, and an electric motor, as shown in Figure 1.In the series HEV, the fuel converter does not directly drive the vehicle wheels.On the contrary, the series HEV uses a generator to convert mechanical power into electrical energy that is stored in the battery and to provide the electric motor power to achieve the torque required to drive the wheels.Using MAs to solve the control strategy optimization problem associated with minimizing engine fuel consumption and emissions while maintaining driving performance saves a considerable amount of time and provides a global optimal solution.The MA proposed for the control strategy optimization correctly and quickly identifies the optimal parameters required for parallel HEVs.Our experimental results further indicate that MAs applied in the design of control strategy are effective in minimizing minimize fuel consumption and emissions while maintaining driving performance.

Configurations of HEVs
Before discussing the configurations of HEVs, we first mention the differences between a conventional HEV and a plug-in HEV.There are two fundamental differences between them, one is the main power source and the other is the whole energy efficiency of the two vehicle architectures.The main power source is a gasoline-powered ICE for conventional HEVs, while it is a batterypowered electric motor for plug-in HEVs.In conventional HEVs, the electric motor is used to complement the ICE and electricity is generated on board.Its energy savings are not substantial in comparison to plug-in HEVs.In plug-in HEVs, the ICE is used to complement the electric motor and grid-supplied electricity.Its energy savings are more substantial than conventional HEVs.Other types of prominent technology [29] are also used.Series and parallel type are currently two common configurations for HEVs.Another common type of HEV synthesizes the characteristics of both series and parallel types and is called dual-mode or multi-mode type [30].

Series HEV
The configuration of a series type HEV mainly consists of a fuel converter/ICE, a generator, an energy storage system/battery, and an electric motor, as shown in Figure 1.In the series HEV, the fuel converter does not directly drive the vehicle wheels.On the contrary, the series HEV uses a generator to convert mechanical power into electrical energy that is stored in the battery and to provide the electric motor power to achieve the torque required to drive the wheels.

Parallel HEV
In the configuration of a parallel type HEV, the fuel converter and the electric motor both transmit power to the vehicle, as shown in Figure 2. The electric motor can also be used as a generator to charge the battery via the regenerative brake or by absorbing the excess power of the fuel converter.In comparisons, the parallel type is better than the series type in the configuration.The parallel type has a smaller fuel converter and a smaller electric motor than the series type while it provides the

Parallel HEV
In the configuration of a parallel type HEV, the fuel converter and the electric motor both transmit power to the vehicle, as shown in Figure 2. The electric motor can also be used as a generator to Energies 2017, 10, 305 4 of 21 charge the battery via the regenerative brake or by absorbing the excess power of the fuel converter.In comparisons, the parallel type is better than the series type in the configuration.The parallel type has a smaller fuel converter and a smaller electric motor than the series type while it provides the same performance.This characteristic makes the parallel type HEV more suitable for heavy-duty passenger vehicles than the series type.
Energies 2017, 10, 305 4 of 21 same performance.This characteristic makes the parallel type HEV more suitable for heavy-duty passenger vehicles than the series type.

Dual-Mode HEV
The configuration of the synthesized series-parallel type includes additional mechanical links as compared to the configuration of the series type, and it has an additional generator as compared to the configuration of the parallel type.That makes the synthesized series-parallel type relatively more complex and expensive than the series and parallel types.

Control Strategy for Parallel HEV
There are two power-drive units (i.e., fuel converter and electric motor) integrated into the parallel type; therefore, the purpose of the control strategy of the parallel HEV is to decide how to allocate the required torque between the fuel converter and the electric motor during driving.When positive torque is requested, the summation of torques for the fuel converter and electric motor is equal to the demand of the driver.On the contrary, when negative torque is requested, the summation of torques for the electric motor and brake is equal to the demand of the driver while the engine torque is zero.
Several control strategies have been employed for parallel HEVs.EACS [4] is the most commonly used control strategy; the fuel converter is the major energy provider, and the electric motor is the assistant component for fuel converter.The baseline control strategy (BCS) is used by the ADvanced VehIcle SimulatOR (ADVISOR) in a parallel HEV.Eight independent input parameters are defined to minimize the engine energy usage, and the variables are also usually defined for an EACS [19], as listed in Table 1.This EACS is a common method of hybrid control; examples of its application are the Toyota Prius [31] and Honda Insight [32].

Dual-Mode HEV
The configuration of the synthesized series-parallel type includes additional mechanical links as compared to the configuration of the series type, and it has an additional generator as compared to the configuration of the parallel type.That makes the synthesized series-parallel type relatively more complex and expensive than the series and parallel types.

Control Strategy for Parallel HEV
There are two power-drive units (i.e., fuel converter and electric motor) integrated into the parallel type; therefore, the purpose of the control strategy of the parallel HEV is to decide how to allocate the required torque between the fuel converter and the electric motor during driving.When positive torque is requested, the summation of torques for the fuel converter and electric motor is equal to the demand of the driver.On the contrary, when negative torque is requested, the summation of torques for the electric motor and brake is equal to the demand of the driver while the engine torque is zero.
Several control strategies have been employed for parallel HEVs.EACS [4] is the most commonly used control strategy; the fuel converter is the major energy provider, and the electric motor is the assistant component for fuel converter.The baseline control strategy (BCS) is used by the ADvanced VehIcle SimulatOR (ADVISOR) in a parallel HEV.Eight independent input parameters are defined to minimize the engine energy usage, and the variables are also usually defined for an EACS [19], as listed in Table 1.This EACS is a common method of hybrid control; examples of its application are the Toyota Prius [31] and Honda Insight [32].T ch cs_charge_trq An alternator-like torque loading on the engine to recharge the battery pack.

T min cs_min_trq_frac
When commanded at a lower torque, the engine will be manipulated at the threshold torque (minimum torque threshold = T min × T max ).Additionally, if SOC < SOC L , the electric motor serves as a generator.
T off cs_off_trq_frac When commanded at a lower torque and if SOC > SOC L , the engine will be shut down (minimum torque threshold = T off × T max ).
ELS L cs_electric_launch_spd_lo The lowest vehicle speed threshold.
ELS H cs_electric_launch_spd_hi The highest vehicle speed threshold.

D ch cs_charge_deplete_bool
To use charge deplete strategy or charge sustaining strategy.
The electric motor is used by the EACS in several specific ways: (1) When the battery state of charge (SOC) is more than SOC L , the engine will be shut down if the required speed is less than the ELS, which is also called the electric launch speed (as shown in Figure 4, Case 1).( 2) When the required torque is less than a minimum torque threshold (T off × T max ), the engine will also be shut down (as shown in Figure 4, Case 2).( 3) The electric motor will provide the total required torque when the engine is shut down (as shown in Figure 4, Cases 1 and 2, and Figure 3, Case 1).( 4) When the battery SOC is less than its SOC L , an alternator-like torque (T ch ) is provided from the engine to charge the battery (as shown in Figure 3, Case 2).This alternator-like charging torque is proportional to the difference between the SOC and the average of the SOC L and SOC H . (5) The engine charging torque is only applied when the engine is started up (as shown in Figure 3, Case 2).( 6) To avoid the engine working at an inefficient low torque status, the engine torque must be maintained at the minimum torque threshold (T min × T max ) (as shown in Figure 3, Case 3).
The fuel consumption, emissions, battery charge, and vehicle performance are greatly influenced by the control strategy parameters.Consequently, we proposed the robust evolutionary computation method MA to optimize the control parameters in this study.

Definition of Objective Function
In the optimization of the control strategy parameters, there are two different approaches to be considered [33,34], and they have already been applied to all types of evolutionary computationbased methods related to this problem [19,21,22].The two different approaches are to optimize a (6) To avoid the engine working at an inefficient low torque status, the engine torque must be maintained at the minimum torque threshold (Tmin × Tmax) (as shown in Figure 4, Case 3).
The fuel consumption, emissions, battery charge, and vehicle performance are greatly influenced by the control strategy parameters.Consequently, we proposed the robust evolutionary computation method MA to optimize the control parameters in this study.

Definition of Objective Function
In the optimization of the control strategy parameters, there are two different approaches to be considered [33,34], and they have already been applied to all types of evolutionary computation-based methods related to this problem [19,21,22].The two different approaches are to optimize a single objective function that is a weighted aggregation of some of the goals and to define the other goals as constraints.Several simultaneous goals including the FC, emissions, charge requirement, and driving performance are managed in the HEV control strategy.In the other words, the HEV control strategy must minimize the weighted sum of the FC and emissions (HC, CO and NO x ) while meeting the charge sustaining requirement and driving performance.We thus define FC, HC, CO, and NO x as a single objective function that is a weighted aggregation as shown in Equation (1): where x is a vector that consists of eight parameters of the control strategy as shown in above-mentioned Table 1.Chromosome encoding in Section 5.The design of MA for optimize control strategy of parallel HEVs in details); w fc , w hc , w co , and w nox are defined as weighting factors employed to respectively investigate the influences of the different objectives of FC, HC, CO, and NO x on the optimization results.Furthermore, we use the PNGV [20] as the dynamic performance requirements to guarantee that the vehicle performance is still maintained during optimization.Table 2 describes the seven PNGV dynamic performance constraints considered.In order to apply the MA directly to the constrained optimization problem, the seven PNGV constraints are integrated into the objective function [35] as shown in Equation ( 2): where x is a vector of the parameters of the control strategy as mentioned above and it is a solution within the solution space S; Objective(x) is the objective function; and c i (x) is the PNGV dynamic performance constraints.

The Design of MA for Optimized Control Strategy of Parallel HEVs
MAs were inspired by Dawkins' conception [27,36].The adjective 'memetic' is derived from the term 'meme', which denotes an analogue to a gene in the context of cultural evolution [26].Similar to GAs, the evolutionary computations involved, such as selection, crossover and mutation are effective in implementing optimal solutions for the control strategy of parallel HEVs.However, in contrast to GAs, MAs progress through a local refinement of individuals before becoming involved in the evolution process and thus assure that all individuals and offspring gain some experiences.After each run, individuals in an MA share information with each other, and superior solutions based on a fitness rule are refined from generation to generation [37].
The flowchart of the proposed MA for the optimize control strategy is shown in Figure 5.There are eight dissociated procedures included: (1) chromosome encoding; (2) population initialization; (3) local search; (4) fitness evaluation; (5) ADVISOR evaluation; (6) judgment on termination conditions; (7) selection, crossover, and mutation processes; and (8) replacement process.These eight procedures are described below.
Energies 2017, 10, 305 7 of 21 In order to apply the MA directly to the constrained optimization problem, the seven PNGV constraints are integrated into the objective function [35] as shown in Equation ( 2): where x is a vector of the parameters of the control strategy as mentioned above and it is a solution within the solution space S; Objective( ) is the objective function; and ( ) is the PNGV dynamic performance constraints.

The Design of MA for Optimized Control Strategy of Parallel HEVs
MAs were inspired by Dawkins' conception [27,36].The adjective 'memetic' is derived from the term 'meme', which denotes an analogue to a gene in the context of cultural evolution [26].Similar to GAs, the evolutionary computations involved, such as selection, crossover and mutation are effective in implementing optimal solutions for the control strategy of parallel HEVs.However, in contrast to GAs, MAs progress through a local refinement of individuals before becoming involved in the evolution process and thus assure that all individuals and offspring gain some experiences.After each run, individuals in an MA share information with each other, and superior solutions based on a fitness rule are refined from generation to generation [37].
The flowchart of the proposed MA for the optimize control strategy is shown in Figure 5.There are eight dissociated procedures included: (1) chromosome encoding; (2) population initialization; (3) local search; (4) fitness evaluation; (5) ADVISOR evaluation; (6) judgment on termination conditions; (7) selection, crossover, and mutation processes; and (8) replacement process.These eight procedures are described below.(1) Chromosome encoding In order to use the MA to optimize the control strategy parameters in parallel HEVs, a chromosome encoding x must first be defined.In this study, we use eight independent control strategy parameters of an EACS as shown in Table 1 as the chromosome encoding for this optimization problem.The chromosome encoding x is a vector with eight dimensions as shown in Equation (3): where x represents a vector that encodes eight parameters of the control strategy.(1) Chromosome encoding In order to use the MA to optimize the control strategy parameters in parallel HEVs, a chromosome encoding x must first be defined.In this study, we use eight independent control strategy parameters of an EACS as shown in Table 1 as the chromosome encoding for this optimization problem.The chromosome encoding x is a vector with eight dimensions as shown in Equation ( 3): where x represents a vector that encodes eight parameters of the control strategy.
Energies 2017, 10, 305 8 of 21 (2) Population initialization When performing the MA algorithm for the control strategy optimization, at the start, dozens of chromosomes x = (SOC L , SOC H , T ch , T min , T off , ELS L , ELS H , D ch ) are randomly generated for an initial population without duplicates.The values of SOC L , SOC H , T ch , T min , T off , ELS L , ELS H and D ch are randomly generated between default bounds given by the EACS constraint.Their default bounds are presented in Table 3. (

3) Local search
The local search is the highly capable of identifying superior individuals (i.e., chromosomes with good fitness) from amongst the neighbors of an original individual.The experience of the original individual will be improved by getting experience from the neighbors, and thus a local optimum solution can be acquired.At the beginning of the algorithm, the local search process is performed by all individuals in order to acquire the local optimality of each individual.Furthermore, the local search process is also performed by new generated individuals in each iteration in order for the local optimum to always be retained.Finally, a global optimum is determined.We applied the pseudo-code of the local search to the proposed method which is as follows.(4) Fitness evaluation The fitness function is a function used to evaluate the fitness of the chromosomes in the individuals.To make use of the MA in the simultaneous optimization of parallel HEVs in the study, the fitness function adapts the inverse of the objective function Objective(x) in Equation ( 1) and the seven PNGV dynamic performance constraints are considered as penalty functions to maximize the fitness value as shown in Equation (4): Energies 2017, 10, 305 9 of 21 where y i are penalty factor selected by trial and error [22] as indicated in Table 2, and p i (x) are penalty functions that are related to the i-th constraint c i (x) in Table 2.The seven penalty functions p i (x), i = 1, 2, 3, . . ., 6, 7 for the PNGV dynamic performance constraints are respectively shown in the following Equations ( 5)- (11).

•
Maximum speed • Maximum acceleration p 5 (x) = 0, if max acceleration satisfied the PNGV max acceleration PNGV max acceleration − max acceleration, if max acceleration < PNGV max acceleration (9) • Distance in 5 s p 6 (x) = 0, if distance in 5 s satisfied the distance in 5 s PNGV distance in 5 s − distance in 5 s, if distance in 5 s < PNGV distance in 5 s (10) • Gradeability (5) ADVISOR evaluation ADVISOR is an abbreviation of "ADvanced VehIcle SimulatOR" which is written in the MATLAB/Simulink environment and was developed by the National Renewable Energy Laboratory in November 1994.It provides the vehicle engineering community with an easy-to-use, flexible, yet robust and supported analysis package for advanced vehicle modeling [38].In order to effectively measure the FC, emissions (HC, CO, and NO x ), and driving performance for the used control strategy in parallel HEVs, we have combined the ADVISOR with the proposed MA to evaluate the fitness in chromosomes.When the fitness function is called, the control strategy variables will be passed on to the ADVISOR to calculate the FC, emissions (HC, CO, and NO x ), and driving performance and return their values.These values will then be used to evaluate the fitness of the control strategy.
Energies 2017, 10, 305 10 of 21 (6) Judgment on termination conditions The termination conditions provide the evolutionary computation method to stop meaningless and endless computations.Here, two termination conditions are employed in the proposed method.One is that the maximum number of iterations is reached, and the other is that the fitness value is no longer updated in several following iterations (i.e., the computation has converged).( 7) Selection, crossover, and mutation processes The election, crossover, and mutation processes are the fundamental operations for the evolutionary computation of an MA.The proposed MA only performs either the crossover or mutation process.In the selection process, two individuals are randomly selected from the population for the following crossover or mutation operations.In the crossover process, uniform crossover is applied to the proposed MA.When the crossover probability is satisfied, the selected two individuals will be processed by the crossover operation.Figure 6a illustrates the flowchart of the crossover process.In the mutation process, one-point mutation is implemented in the proposed MA.When the mutation probability is satisfied, the selected two individuals will be processed by the mutation operation.Figure 6b illustrates the flowchart of the mutation process.
Energies 2017, 10, 305 10 of 21 (6) Judgment on termination conditions The termination conditions provide the evolutionary computation method to stop meaningless and endless computations.Here, two termination conditions are employed in the proposed method.One is that the maximum number of iterations is reached, and the other is that the fitness value is no longer updated in several following iterations (i.e., the computation has converged).( 7) Selection, crossover, and mutation processes The election, crossover, and mutation processes are the fundamental operations for the evolutionary computation of an MA.The proposed MA only performs either the crossover or mutation process.In the selection process, two individuals are randomly selected from the population for the following crossover or mutation operations.In the crossover process, uniform crossover is applied to the proposed MA.When the crossover probability is satisfied, the selected two individuals will be processed by the crossover operation.Figure 6a illustrates the flowchart of the crossover process.In the mutation process, one-point mutation is implemented in the proposed MA.When the mutation probability is satisfied, the selected two individuals will be processed by the mutation operation.Figure 6b illustrates the flowchart of the mutation process.In order to update the individuals of the population from iteration to iteration, the replacement process is an essential and crucial operation.In the proposed MA, the worst individuals in the population will be replaced with the new individuals.This process is repeated in each next generation until one of the termination conditions is met.(8) Replacement process In order to update the individuals of the population from iteration to iteration, the replacement process is an essential and crucial operation.In the proposed MA, the worst individuals in the population will be replaced with the new individuals.This process is repeated in each next generation until one of the termination conditions is met.

Driving Cycles for Optimized Control Strategy of Parallel HEVs
Four different driving cycles, including NEDC, FTP, ECE + EUDC, and UDDS are carried out by the proposed method to find out the optimal control parameters for each driving cycle.The parameters used are listed in Table 4.
The NEDC eliminates the idling period, which causes the engine to start at 0 s, and the emission sampling begins at the same time as the entire cycle [39].The FTP is used for regulatory emission testing of heavy-duty on-road engines in the United States (40 CFR 86.1333).It is based on the UDDS chassis dynamometer driving cycle and includes "motoring" segments, and therefore, requires a DC or AC electric dynamometer capable of both absorbing and supplying power [40].The ECE + EUDC is performed on a chassis dynamometer and was used for EU type approval testing of emissions and fuel consumption in light-duty vehicles [EEC Directive 90/C81/01].This cycle is used for the emission certification of light-duty vehicles in Europe and is also known as the MVEG-A cycle.The entire cycle includes four ECE segments repeated without interruption, followed by one EUDC segment.Before the test, the vehicle is allowed to soak for at least 6 h at a test temperature of 20-30

Parameter Settings
Four parameters are set in the proposed method, i.e., the number of generations, the population size, the probability of crossover, and the probability of mutation.Their values were set to 50, 20, 1.0, and 1.0, respectively.The parameters are chosen by "trial and error."Furthermore, our experience in algorithmic design has also helped us to select suitable parameters for the problem.The published papers for the algorithms used by the authors can also be referred to [36,37,[42][43][44][45][46][47][48][49].Furthermore, the components for the parallel HEV model are set as shown in Table 5, and the parameters for the used parallel HEV model are set as shown in Table 6.The versions and types are the parameters used in ADVISOR.The variables are the parameters used in ADVISOR.

Results for the Used Driving Cycles
Table 7 shows the initial values and optimal values obtained, and Figure 7 shows the histogram for fitness values obtained using the proposed MA method with the parameters of the control strategy for four driving cycles: NEDC, FTP, ECE + EUDC, and UDDS.The result for the FC and emissions is presented in Table 8 and Figure 8. Further, the result for the dynamic performance is presented in Table 9 and Figure 9. On comparing the results presented in Table 8, Figure 8, Table 9, and Figure 9, we observe that the results obtained using the initial and optimal values of the parameters of the control strategy are very different.The optimal values of the parameters of the control strategy found using the proposed MA method indeed help to reduce the integral FC, HC, CO and NO x emissions.Furthermore, the dynamic performance not only simultaneously meets the PNGV constraints but also exceeds its requirements.From the above results, the ability of the proposed MA method to optimize the control parameters in parallel HEVs is confirmed.The results used in this study are freely accessible and provided online for those interested [50].6.5 6.5 -6.5 6.5 6.5 -6.5    6.5 6.5 -6.5 6.5 6.5 -6.5

Improvement in Overall Effectiveness
First, for the observation of the fitness values as shown in Table 7 and Figure 7, the fitness values of the optimal results are 0.011, 0.013, 0.008, and 0.012, which are better than those at the initial conditions in the driving cycles of NEDC, FTP, ECE + EUDC, and UDDS, respectively (as shown in Figure 10).This shows that the proposed MA method is quite useful for selecting a feasible control strategy from the given bounds in the EACS constraints (see Table 3).The improvement in the fitness value means the improvement in overall effectiveness, including the reduction of the FC and emissions, and the improvement or maintain of dynamic performance.Different driving cycles show distinct degrees of improvement owing to the different parameters of each driving cycle (see Table 4).

Improvement in Overall Effectiveness
First, for the observation of the fitness values as shown in Table 7 and Figure 7, the fitness values of the optimal results are 0.011, 0.013, 0.008, and 0.012, which are better than those at the initial conditions in the driving cycles of NEDC, FTP, ECE + EUDC, and UDDS, respectively (as shown in Figure 10).This shows that the proposed MA method is quite useful for selecting a feasible control strategy from the given bounds in the EACS constraints (see Table 3).The improvement in the fitness value means the improvement in overall effectiveness, including the reduction of the FC and emissions, and the improvement or maintain of dynamic performance.Different driving cycles show distinct degrees of improvement owing to the different parameters of each driving cycle (see Table 4).

Improvement in FC and Emissions
For the FC and emissions (HC, CO, and NOx) as shown in Table 8 and Figure 8, we observe that the FC is improved for the driving cycles of the NEDC and FTP.The FC is improved by 0.72, 5.49, −0.36, and −0.05 mi/gal for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.
The HC emission is reduced for all the driving cycles except UDDS.The HC emission is reduced by 0.10, 0.03, 0.02, and −0.04 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The CO emission is reduced for all the driving cycles.The CO emission is substantially reduced by 5.12, 7.62, 0.16, and 1.92 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The NOx emission is improved for the driving cycles of the NEDC and ECE + EUDC.The NOx emission is reduced by 0.16, −0.15, 0.15, and −0.03 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The improvement results for the FC and emissions are shown in Figure 11.Although some values of the parameters are not greatly improved, their improved values are reflected in the other parameters or the following resultant dynamic performance.
. The histogram of the improvement in overall effectiveness between the results using initial and optimal values of the parameters of the control strategy.

Improvement in FC and Emissions
For the FC and emissions (HC, CO, and NO x ) as shown in Table 8 and Figure 8, we observe that the FC is improved for the driving cycles of the NEDC and FTP.The FC is improved by 0.72, 5.49, −0.36, and −0.05 mi/gal for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.
The HC emission is reduced for all the driving cycles except UDDS.The HC emission is reduced by 0.10, 0.03, 0.02, and −0.04 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The CO emission is reduced for all the driving cycles.The CO emission is substantially reduced by 5.12, 7.62, 0.16, and 1.92 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The NO x emission is improved for the driving cycles of the NEDC and ECE + EUDC.The NO x emission is reduced by 0.16, −0.15, 0.15, and −0.03 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The improvement results for the FC and emissions are shown in Figure 11.Although some values of the parameters are not greatly improved, their improved values are reflected in the other parameters or the following resultant dynamic performance.

Improvement in FC and Emissions
For the FC and emissions (HC, CO, and NOx) as shown in Table 8 and Figure 8, we observe that the FC is improved for the driving cycles of the NEDC and FTP.The FC is improved by 0.72, 5.49, −0.36, and −0.05 mi/gal for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.
The HC emission is reduced for all the driving cycles except UDDS.The HC emission is reduced by 0.10, 0.03, 0.02, and −0.04 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The CO emission is reduced for all the driving cycles.The CO emission is substantially reduced by 5.12, 7.62, 0.16, and 1.92 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The NOx emission is improved for the driving cycles of the NEDC and ECE + EUDC.The NOx emission is reduced by 0.16, −0.15, 0.15, and −0.03 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The improvement results for the FC and emissions are shown in Figure 11.Although some values of the parameters are not greatly improved, their improved values are reflected in the other parameters or the following resultant dynamic performance.With respect to the dynamic performance, all the optimal values of the control strategy satisfy the seven PNGV dynamic performance constraints and the dynamic performance shows improvements from the initial conditions as shown in Table 2. Table 9 and Figure 9 show the improved results.The improvement in the dynamic performance is shown in Figure 10.For the NEDC, the acceleration time is reduced by 4.98, 2.96, and 13.29 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased by 0.003 mi/s, the maximum acceleration is maintained at the original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.005 mi, and the gradeability is maintained at 6.5%.For the FTP, the acceleration time is reduced by 3.12, 1.94, and 8.56 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased by 0.002 mi/s, the maximum acceleration is maintained at the original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.003 mi, and the gradeability increased from 0% to 6.5%.For the ECE + EUDC, the acceleration time is reduced by 3.79, 2.35, and 10.31 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased by 0.003 mi/s, the acceleration is maintained at its original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.004 mi, and the gradeability is maintained at the original value of 6.5%.Finally, for the UDDS, the acceleration time is reduced by 5.61, 3.27, and 15.08 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased in 0.003 mi/s, the maximum acceleration is maintained at its original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.006 mi, and the gradeability improves from 0% to 6.5%.There improvements in the dynamic performance are shown in Figure 12.
Energies 2017, 10, 305 The histogram of the improvements in FC and emissions between the results using initial and optimal values of the parameters of control strategy.

Improvement in Dynamic Performance
With respect to the dynamic performance, all the optimal values of the control strategy satisfy the seven PNGV dynamic performance constraints and the dynamic performance shows improvements from the initial conditions as shown in Table 2. Table 9 and Figure 9 show the improved results.The improvement in the dynamic performance is shown in Figure 10.For the NEDC, the acceleration time is reduced by 4.98, 2.96, and 13.29 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased by 0.003 mi/s, the maximum acceleration is maintained at the original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.005 mi, and the gradeability is maintained at 6.5%.For the FTP, the acceleration time is reduced by 3.12, 1.94, and 8.56 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased by 0.002 mi/s, the maximum acceleration is maintained at the original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.003 mi, and the gradeability increased from 0% to 6.5%.For the ECE + EUDC, the acceleration time is reduced by 3.79, 2.35, and 10.31 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased by 0.003 mi/s, the maximum acceleration is maintained at its original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.004 mi, and the gradeability is maintained at the original value of 6.5%.Finally, for the UDDS, the acceleration time is reduced by 5.61, 3.27, and 15.08 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively.The maximum speed is increased in 0.003 mi/s, the maximum acceleration is maintained at its original value of 0.003 mi/s 2 , the distance traveled in 5 s is increased by 0.006 mi, and the gradeability improves from 0% to 6.5%.There improvements in the dynamic performance are shown in Figure 12.

Conclusions
This study proposes a robust evolutionary computation method called a "memetic algorithm (MA)" to optimize the control parameters in parallel HEVs.Its "local search" mechanism greatly enhances the search capabilities for obtaining optimal values of the parameters of the control strategy.

Conclusions
This study proposes a robust evolutionary computation method called a "memetic algorithm (MA)" to optimize the control parameters in parallel HEVs.Its "local search" mechanism greatly enhances the search capabilities for obtaining optimal values of the parameters of the control strategy.We combine the ADVISOR with the fitness function in the proposed MA method and use it to effectively evaluate the FC and emissions (HC, CO, and NO x ) of the vehicle engine based on the EACS.Furthermore, we consider that the driving performance requirements, which are referred to as the PNGV constraints, must be satisfied simultaneously.In order to verify the proposed method, we have employed four different driving cycles-NEDC, FTP, ECE + EUDC, and UDDS; the results for these driving cycles are verified.The fitness values of the optimal results are 0.011, 0.013, 0.008, and 0.012, which are better than those of the initial conditions in the driving cycles of the NEDC, FTP, ECE + EUDC, and UDDS, respectively.For the FC and emissions (HC, CO, and NO x ), the FC is improved by 0.72 and 5.49 mi/gal for the NEDC and FTP, respectively.The HC emission is reduced by 0.10, 0.03, and 0.02 g/mi for the NEDC, FTP, and ECE + EUDC, respectively.The CO emission is substantially reduced by 5.12, 7.62, 0.16, and 1.92 g/mi for the NEDC, FTP, ECE + EUDC, and UDDS, respectively.The NO x emission is a reduced by 0.16 and 0.15 g/mi for the NEDC and ECE + EUDC, respectively.With regard to the dynamic performance, all the optimal values of the control strategy satisfy the seven PNGV dynamic performance constraints.The acceleration time is reduced by 4.98, 2.96, and 13.29 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively; the maximum speed is increased by 0.003 mi/s; and the distance traveled in 5 s is increased by 0.005 mi for the NEDC.The acceleration time is respectively reduced in 3.12, 1.94, and 8.56 s for 0-60 mph, 40-60 mph, and 0-85 mph; the maximum speed is increased in 0.002 mi/s; the distance traveled in 5 s is increased in 0.003 mi for the FTP.The acceleration time is reduced by 3.79, 2.35, and 10.31 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively; the maximum speed is increased by 0.003 mi/s; and the distance traveled in 5 s is increased by 0.004 mi for the ECE + EUDC.The acceleration time is reduced by 5.61, 3.27, and 15.08 s for 0-60 mph, 40-60 mph, and 0-85 mph, respectively; the maximum speed is increased by 0.003 mi/s; and the distance traveled in 5 s is increased by 0.006 mi for the UDDS.Furthermore, the maximum acceleration is maintained at 0.003 mi/s 2 and the gradeability reached a value of 6.5% in all the considered driving cycles.The results indicate that the proposed MA method is capable of determining the optimal parameters of the control strategy for parallel HEVs.It helps to improve the fuel consumption and reduce the emissions, as well as guarantee vehicle performance.

Figure 1 .
Figure 1.Series configuration of a hybrid electric vehicle (HEV).

Figure 2 .
Figure 2. Parallel configuration of an HEV.

Figure 3 .
Figure 3. Electric assist control strategy (EACS) for high state of charge (SOC).Figure 4. Electric assist control strategy (EACS) for high state of charge (SOC).

Figure 4 .
Figure 3. Electric assist control strategy (EACS) for high state of charge (SOC).Figure 4. Electric assist control strategy (EACS) for high state of charge (SOC).

Figure 5 .
Figure 5. Flowchart of memetic algorithm (MA) for optimize control strategy.

Figure 5 .
Figure 5. Flowchart of memetic algorithm (MA) for optimize control strategy.

Figure 6 .
Figure 6.Flowchart for the (a) crossover process; and (b) mutation process.

Figure 6 .
Figure 6.Flowchart for the (a) crossover process; and (b) mutation process.

Figure 7 .
Figure 7.The histogram of fitness values for initial and optimal values of the parameters of the control strategy.

Figure 8 .
Figure 8.The histogram for FC and emissions using initial and optimal values of the parameters of control strategy.(A) FC and emissions based on NEDC cycle; (B) FC and emissions based on FTP cycle; (C) and emissions based on ECE + EUDC cycle, and (D) FC and emissions based on UDDS cycle.

Figure 8 .
Figure 8.The histogram for FC and emissions using initial and optimal values of the parameters of control strategy.(A) FC and emissions based on NEDC cycle; (B) FC and emissions based on FTP cycle; (C) and emissions based on ECE + EUDC cycle; and (D) FC and emissions based on UDDS cycle.

Figure 8 .
Figure 8.The histogram for FC and emissions using initial and optimal values of the parameters of control strategy.(A) FC and emissions based on NEDC cycle; (B) FC and emissions based on FTP cycle; (C) and emissions based on ECE + EUDC cycle, and (D) FC and emissions based on UDDS cycle.

Figure 9 .
Figure 9.The histogram for dynamic performance using initial and optimal values of the parameters of control strategy.(A) dynamic performance based on NEDC cycle; (B) dynamic performance based on FTP cycle; (C) dynamic performance based on ECE + EUDC cycle, and (D) dynamic performance based on UDDS cycle.

Figure 9 .
Figure 9.The histogram for dynamic performance using initial and optimal values of the parameters of control strategy.(A) dynamic performance based on NEDC cycle; (B) dynamic performance based on FTP cycle; (C) dynamic performance based on ECE + EUDC cycle; and (D) dynamic performance based on UDDS cycle.

Figure 10 .
Figure10.The histogram of the improvement in overall effectiveness between the results using initial and optimal values of the parameters of the control strategy.

Energies 2017, 10 , 305 16 of 21 Figure 10 .
Figure10.The histogram of the improvement in overall effectiveness between the results using initial and optimal values of the parameters of the control strategy.

Figure 12 .
Figure 12.The histogram of the improvements in the dynamic performance between the results using initial and optimal values of the parameters of the control strategy.

Figure 12 .
Figure 12.The histogram of the improvements in the dynamic performance between the results using initial and optimal values of the parameters of the control strategy.

Table 1 .
Eight independent input parameters/variables for an electric assist control strategy (EACS).

Table 3 .
Range of the parameters of the EACS constraints for the proposed method.

Local Search Pseudo-Code for the Proposed Method
• C. It is then started and allowed to idle for 40 s [39].The UDDS is also called FTP-72 (Federal Test Procedure) or LA-4 cycle.It is used for light-duty vehicle testing.The cycle simulates an urban route of 7.5 miles with frequent stops.The maximum speed is 56.7 mph, and the average speed is 19.6 mph.Two phases are included in the cycle: (1) 505 s (3.59 miles at 25.6 mph average speed); and (2) 867 s [41].
Data source: ADVISOR, United States Environmental Protection Agency (EPA), and DieselNet.

Table 5 .
Components for the parallel HEV model.

Table 6 .
Parameters for the used parallel HEV model.

Table 7 .
The initial values and optimal values of the parameters of the control strategy.

Table 7 .
The initial values and optimal values of the parameters of the control strategy.The histogram of fitness values for initial and optimal values of the parameters of the control strategy.

Table 8 .
The result for fuel consumption (FC) and emissions using initial and optimal values of the parameters of the control strategy.

Table 8 .
The result for fuel consumption (FC) and emissions using initial and optimal values of the parameters of the control strategy.

Table 9 .
The result for the dynamic performance using initial and optimal values of the parameters of the control strategy.

Table 9 .
The result for the dynamic performance using initial and optimal values of the parameters of the control strategy.

Table 9 .
The result for the dynamic performance using initial and optimal values of the parameters of the control strategy.
The histogram of the improvements in FC and emissions between the results using initial and optimal values of the parameters of control7.3.3.Improvement in Dynamic Performance