Next Article in Journal
Air Quality Index and Air Pollutant Concentration Prediction Based on Machine Learning Algorithms
Next Article in Special Issue
Performance Analysis of a Hybrid Electric Vehicle with Multiple Converter Configuration
Previous Article in Journal
The Ultrafast Laser Ablation of Li(Ni0.6Mn0.2Co0.2)O2 Electrodes with High Mass Loading
Previous Article in Special Issue
Optimal Charging of Plug-In Electric Vehicle: Considering Travel Behavior Uncertainties and Battery Degradation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Optimization for Plug-In 4WD Hybrid Electric Vehicle Powertrain

1
School of Traffic & Transportation Engineering, Changsha University of Science & Technology, Changsha 410114, China
2
Jiangxi Province Key Laboratory of Precision Drive & Control, Nanchang Institute of Technology, Nanchang 330099, China
3
National Demonstrating Center for Experimental Civil Engineering Education, Hunan City University, Yiyang 410003, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(19), 4068; https://doi.org/10.3390/app9194068
Submission received: 23 August 2019 / Revised: 22 September 2019 / Accepted: 24 September 2019 / Published: 29 September 2019
(This article belongs to the Special Issue Hybrid Vehicle Technologies for a Sustainable Future Mobility)

Abstract

:
This paper focuses on the parameter optimization for the CVT (a continuously variable transmission) based plug-in 4WD (4-wheel drive) hybrid electric vehicle powertrain. First, the plug-in 4WD hybrid electric vehicle (plug-in 4WD HEV)’s energy management strategy based on the CD (charge depleting) and CS (charge sustain) mode is developed. Then, the multi-objective optimization’s mathematical model, which aims at minimizing the electric energy consumption under the CD stage, the fuel consumption under the CS stage and the acceleration time from 0–120 km/h, is established. Finally, the multi-objective parameter optimization problem is solved using an evolutionary based non-dominated sorting genetic algorithms-II (NSGA-II) approach. Some of the results are compared with the original scheme and the classical weight approach. Compared with the original scheme, the best compromise solution (i.e., electric energy consumption, fuel consumption and acceleration time) obtained using the NSGA-II approach are reduced by 1.21%, 6.18% and 5.49%, respectively. Compared with the weight approach, the Pareto optimal solutions obtained using NSGA-II approach are better distributed over the entire Pareto optimal front, as well as the best compromise solution is also better.

1. Introduction

Owing to regulations on fuel economy and emissions become more and more stringent, the development of electrified vehicles in recent years have been a surging trend [1,2]. Plug-in hybrid electric vehicles (PHEVs), which achieve a longer all-electric range compared to conventional hybrid electric vehicles and have no driver range anxiety compared to pure battery electric vehicles, become an important research direction in the field of electric vehicles [3,4].
The parameters of PHEV’s powertrain have a significant impact on the electric energy consumption, fuel consumption and dynamic performance of the vehicle. Therefore, parameter optimization of the powertrain is the basis of the vehicle development. The research on PHEV’s parameter optimization has gone through the following process. In the initial stage, scholars did not realize the problem that powertrain parameters and control strategy parameters coupled and affected vehicles’ performance together. Therefore, the parameter optimization during this stage is mainly focused on independent optimization of powertrain parameters, instead of joint optimizing the parameters of powertrain and control strategy [5,6]. With the in-depth research, attention has been paid to the coupling problem of powertrain parameters and control strategy parameters. Then, simultaneously optimization of powertrain parameters and control strategy parameters begins to appear in some literature [7,8,9]. Actually, literatures [5,6,7,8,9] take parameter optimization as a single objective optimization problem, the parameters of PHEV are optimized to only improve the economy of the vehicle. However, the objectives of vehicle economy, power system cost, emission and power performance usually conflict with each other. The optimization only considering the vehicle economy will result in achieving optimal economic performance at the expense of other objectives that conflict with it. Therefore, in recent years, multi-objective optimization has been applied to optimize PHEV’s parameters. In reference [10], energy storage system costs and fuel economy are considered as the objective functions for hybrid electric vehicle’s battery size optimization. Reference [11] presents a multi-objective optimization methodology for vehicles design considering the parameters for design and macro level operating strategy. Reference [12] applies multi-objective algorithms (minimization of the couples cost and fuel, cost and LCA (life cycle impact) CO2eq, fuel and LCA CO2eq) to perform the powertrain components optimization. Reference [13] takes powertrain cost, fuel consumption and emission as multi-objectives to optimize the parameters of powertrain and control parameters simultaneously. In [14], the usage cost, acceleration performance and mode discrimination of the whole vehicle under a new European driving cycle (NEDC) are considered as the optimization objectives to optimize the transmission system of a new hybrid power system. Reference [15] aims at minimizing the cost of the power supply and the energy flow of batteries, and adopts the convex optimization algorithm to optimize the power supply’s parameters and energy management strategy.
Multi-objective optimization and its applications have been an important area of research for over two decades now [16]. The key of solving the multi-objective optimization problem is how to get the Pareto solution set. Traditional algorithms for solving multi-objective optimization problems include the linear weighting method [17], constraint method [18], mini–max method [19], goal programming method [20] and goal satisfaction method [21]. In order to get the Pareto solution set, these methods need to run many times, which reduces the efficiency of the algorithm. Evolutionary algorithms are a group-based global optimization algorithm that simulate the evolution process of natural organisms [22]. The evolutionary algorithm has achieved great success in solving power system’s multi-objective optimization problems [16,23,24,25], e.g., a strength Pareto evolutionary algorithm is used to solve the multi-objective reactive power price clearing problem [23]; a novel-efficient evolutionary-based multi-objective optimization approach is proposed to solve multi-objective optimal power flow problems [16,24] and the multi-objective strength Pareto evolutionary algorithm 2+ has been employed to solve the congestion management problem [25]. Therefore, it has become an important method for solving multi-objective optimization problems.
Parameter optimization of PHEVs is related to the energy management strategy, and energy management strategies need to be developed before parameter optimization. At present, the research on the energy management strategy for PHEVs mainly focuses on the development of the advanced optimization algorithm, such as the algorithm based on the minimum equivalent fuel consumption [26,27,28,29], the dynamic programming algorithm [30,31,32], stochastic dynamic programming [33], the algorithm based on convex optimization [2,34] and the model predictive control algorithm [35,36,37,38]. Although the above-mentioned optimization algorithm can obtain the local or global optimum, it is difficult to apply to real vehicle control for hardly knowing the driving cycles beforehand or the large amount of calculation. However, the rule control strategy based on the charge depleting–charge sustain (CD–CS) mode does not need to know the driving cycles beforehand, and the calculation is small, so it is widely used in the real vehicle control of PHEVs.
As a summary of the entire literature review, in order to complete the plug-in 4WD HEV’s parameter optimization well, the simultaneous optimization for the main parameters of powertrain and control strategy is necessary, multi-objective optimization should be taken into account and the rule control strategy based on the CD–CS mode for plug-in 4WD HEVs should be developed. However, there is still a large shortage for plug-in 4WD HEV’s parameter optimization. Firstly, it is still challenging to select reasonably multi-objective functions, which can evaluate the optimal performance of plug-in 4WD HEV well. Secondly, although the evolutionary algorithm has achieved great success in solving power system’s multi-objective optimization problems, it rarely applies to solve plug-in 4WD HEV’s multi-objective parameter optimization problems. Finally, the rule control strategy based on the CD–CS mode for plug-in 2WD hybrid electric vehicles is mature, but there is little literature on the rule control strategy based on the CD–CS mode for plug-in 4WD HEVs.
To address the challenges summarized above, in this study, the energy management strategy based on the CD–CS mode for plug-in 4WD HEV is developed. Then, the reasonably multi-objective functions, which are composed of electric energy consumption under the CD stage, fuel consumption under the CS stage and acceleration time from 0–120 km/h, are established. Finally, the evolutionary based NSGA-II (non-dominated sorting genetic algorithms-II) approach is selected to simultaneously optimize the parameters of the powertrain and control strategy.
The outline of this paper is as follows. The structure and dynamic model of the powertrain are provided in Section 2. The energy management strategy based on the CD–CS mode is developed in Section 3. Mathematical model of multi-objective optimization is built in Section 4. Optimization algorithm is proposed in Section 5. Optimization results are discussed in Section 6. Finally, conclusions are summarized in Section 7.

2. The Structure and Dynamic Model of the Powertrain

2.1. Structure of the Plug-In 4WD Hybrid Electric Vehicle

This study focuses on a plug-in 4WD HEV. Its structure layout is shown in Figure 1. Its powertrain is mainly composed of gasoline engine, an integrated starter and generator motor (ISG motor), rear-drive motor, CVT, one-way clutch C1, the front final drive and the rear final drive. The ISG motor mounted on the gearbox input (referred to as P2 [39]), and the rear-drive motor mounted on the rear driving axle (referred to as P4 [39]).

2.2. The Dynamic Model of the Powertrain

The powertrain dynamic model is mainly applied to the development of the energy management strategy and the evaluation of economy and dynamic performance. Therefore, the model was deduced by quasi-static modeling technology [40,41] in this study. The total driving force at the wheel is described as:
F t ( t ) = m g f r cos ( α ( t ) ) + m g sin ( α ( t ) ) + c d A 21.15 v ( t ) 2 + m d v ( t ) d t
where m , c d , A and f r are the vehicle mass, the coefficient of air resistance, the windward area and the rolling resistance coefficient, respectively. v ( t ) and α ( t ) are the vehicle speed and road gradient, respectively.
During deceleration, F t ( t ) is negative, and the front and the rear axle braking force’s distribution follows the fixed factor β . Assuming the front and the rear axle braking force are F f b ( t ) and F r b ( t ) , respectively. The braking force are calculated by:
{ F t ( t ) = F f b ( t ) + F r b ( t ) F f b ( t ) = β F t ( t ) F r b ( t ) = ( 1 β ) F t ( t )
Then, the torque and power demand at the wheel is described as:
{ T r e q ( t ) = F t ( t ) R P r e q ( t ) = F t ( t ) v ( t )
where R is the wheel’s radius.
The torque and angular speed relationship of powertrain elements can be described as:
{ ( T e ( t ) + T i s g ( t ) ) i c v t i f o + T m ( t ) i r o = T r e q ( t ) ω e ( t ) = ω i s g ( t ) = ω w i c v t i f o ω m = ω w i r o
where T e , T i s g and T m are the output torque of the engine, ISG motor and the rear-drive motor, respectively. ω e , ω i s g , ω m and ω w are the angular speed of the engine, ISG motor, the rear-drive motor and the wheel, respectively. i c v t , i f o and i r o are the gear ratio of CVT, the front final drive and the rear final drive, respectively.
The electric power demanded or generated by the ISG motor and the rear-drive motor is described as:
{ P b ( t ) = P b i ( t ) + P b m ( t ) P b i ( t ) = { P i s g ( t ) / η i s g    P i s g ( t ) 0 P i s g ( t ) η i s g    P i s g ( t ) < 0   P b m ( t ) = { P m ( t ) / η m    P m ( t ) 0 P m ( t ) η m    P m ( t ) < 0  
where P b ( t ) , P b i ( t ) , P b m ( t ) , P i s g ( t ) and P m ( t ) are the battery terminal power, the electric power demanded or generated by the ISG motor, the electric power demanded or generated by the rear-drive motor, ISG motor’s output power and the rear-drive motor’s output power, respectively. η i s g and η m are the efficiency of the ISG motor and the rear-drive motor, respectively. These efficiencies are related to the motors’ working point and obtained by looking up tables, and the base motor/generator for the ISG motor is shown in Figure 2a with a maximum torque of 115 Nm and maximum speed of 6000 rpm, the base motor/generator for the rear-drive motor is shown in Figure 2b, with a maximum torque of 145 Nm and maximum speed of 9200 rpm.
Then, the battery’s state of charge (SOC) is calculated by:
{ I ( t ) = V o c + V o c 2 4 R int P b ( t ) 2 R int S O C ( t + 1 ) = S O C ( t ) I ( t ) Δ t Q 0
where V o c and R int are the battery’s terminal voltage and the battery’s internal resistance, respectively. They are related to SOC and obtained by looking up tables according to SOC, as shown in Figure 3a Q 0 is the nominal battery capacity and Δ t is the sample time.
The fuel consumption can be solved by looking up the table, and the base engine is shown in Figure 3b.

3. The Energy Management Strategy Based on the CD–CS Mode

The energy management strategy based on the CD–CS mode for a plug-in 4WD HEV has two operating modes (CD mode and CS mode). When SOC ≥ SOCcd (SOCcd is the switch threshold between the CD mode and the CS mode) is satisfied, the vehicle operates in the CD mode. However, when the SOC drops to the lower limit, the vehicle runs in the CS mode to sustain the battery SOC. The CD and CS operating modes switching strategy is shown in Figure 4, where SOCo is the margin of switching back to the CD mode from the CS mode. In this paper, the value of SOCcd and SOCo were set to 0.3 and 0.05, respectively.

3.1. CD Mode

When SoC ≥ SOCcd, the vehicle operates in the CD mode, its corresponding control strategy flow is shown in Figure 5. V, Treq, Tmmax and Tisgmax are the vehicle’s speed, the demand torque at the wheel, the rear-drive motor’s maximum output torque and the ISG motor’s maximum output torque, respectively. ifo, iro and icvt are the speed ratio of the front final drive, the rear final drive and CVT, respectively. Their definitions are as follows: i f o = ω c v t _ o u t / ω w , i r o = ω m / ω w and i c v t = ω i s g / ω c v t _ o u t , ω i s g , ω m and ω w are the same as Equation (4) and ω c v t _ o u t is the angular speed of the CVT output shaft. Tmmax_w is the equivalent wheel torque of rear-drive motor’s maximum output torque, it is expressed as Tmmax_w = Tmmax × iro. Tisgmax_w is the equivalent wheel torque of ISG motor’s maximum output torque, it is expressed as Tisgmax_w = Tisgmax × ifo × icvt × ηcvt. There are three drive modes (rear motor drive mode, dual motor drive mode and 4WD hybrid mode). In the rear motor drive mode, the vehicle is only driven by the rear-drive motor. In the dual motor drive mode, the vehicle is jointly driven by the rear-drive motor and the ISG motor. The rear-drive motor, the ISG motor and engine work jointly to meet the driving demand in the 4WD hybrid mode. If 0 ≤ Treq ≤ Tmmax_w is satisfied, the vehicle drives in the rear motor drive mode. When the vehicle runs at a low speed, and if Treq > Tmmax_w is satisfied, the vehicle drives in the dual motor drive mode. When the vehicle speed is higher than the predetermined speed Vo, and if Tmmax_w ≤ Treq ≤ (Tmmax_w + Tisgmax_w) is satisfied, the vehicle drives in the dual motor drive mode, otherwise the vehicle drives in the 4WD hybrid mode. Regenerative braking can be employed to recover energy in the braking mode, the braking strategy will be introduced in Section 3.3.

3.2. CS Mode

When SoC ≤ SOCcd, vehicle operates in the CS mode. The control strategy of the CS mode is shown in Figure 6. SOCl, Temax, Tel and Teh are the lower limit of SOC in the CS mode, the engine’s maximum output torque, engine’s minimum output torque and engine output torque’s upper limit under the pure engine drive mode, respectively. Temax_w, Tel_w and Teh_w are the equivalent wheel torque of engine’s maximum output torque, engine’s minimum output torque and engine’s upper limit output torque, respectively. They are expressed as Temax_w = Temax × ifo × icvt × ηcvt, Tel_w = Tel × ifo × icvt × ηcvt and Teh_w = Teh × ifo × icvt × ηcvt, respectively. Tel and Teh are obtained by:
{ T e h = T e o p t + k u p ( T e max T e o p t ) T e l = k l o w T e o p t
where T e o p t is the output torque corresponding to the engine’s minimum fuel consumption rate at a certain speed and k u p and k l o w are the control strategy parameters that need to be optimized. In addition to satisfying Equation (7), T e h and T e l also need to satisfy these inequality constraints: T e o p t < T e h T e max and 0 < T e l T e o p t , therefore, these parameters k u p and k l o w are less than 1, and their ranges are usually set to [0.1,1]. The initial values of these parameters are generated in the population initialization process of the evolutionary algorithm.
There are four drive modes (rear motor drive mode, pure engine drive mode, 2WD hybrid mode and charging mode). When the SOC is less than SOCl, it needs to charge the battery packs, if the demand torque is less than Temax_w, the vehicle drives in the charging mode, the torque distribution in this mode is carried out by:
{ T e = T e l , T m = 0 , T i s g = max ( T i s g _ min , ( T r e q / ( i f o i c v t η c v t ) T e ) )   if   0 T r e q T e l _ w T e = T e o p t , T m = 0 , T i s g = max ( T i s g _ min , ( T r e q / ( i f o i c v t η c v t ) T e ) )   if   T e l _ w < T r e q T e o p t i f o i c v t η c v t T e = T e max , T m = 0 , T i s g = max ( T i s g _ min , ( T r e q / ( i f o i c v t η c v t ) T e ) )   if   T e o p t i f o i c v t η c v t < T r e q T e max _ w  
When the SOC is higher than SOCl, and if Treq ≤ Tel_w is satisfied, the vehicle drives in the rear motor drive mode, if Tel_w ≤ Treq ≤ Teh_w is satisfied, the pure engine drive mode is adopted, otherwise the vehicle drives in the 2WD hybrid mode.

3.3. The Braking Strategy

The braking strategy of the above CD and CS mode is shown in Figure 7. In the figure, Z is the braking strength and Zm is the braking strength, which can be provided by the rear-drive motor’s peak torque for the whole vehicle, its value is 0.09. Freq is the total braking force on the wheel. Ffb is the braking force distributed to the front axle. Frb is the braking force distributed to the rear axle. β is the braking force’s distribution factor. Fm and Fisg are the braking force distributed to the rear-drive motor and ISG motor, respectively. Fmmax and Fisgmax are the maximum braking force, which can be distributed to rear-drive motor and ISG motor, respectively. Fμ1 and Fμ2 are the braking force distributed to the front wheel hydraulic brake and the rear wheel hydraulic brake, respectively.
The rear axle has no transmission, so in the low braking strength case (0 < Z ≤ Zm), only the rear-drive motor braking is used to ensure high energy recovery efficiency. When the braking strength is higher than Zm, the braking force is distributed to the front and rear axles. If the braking strength is between Zm and 0.7, the motor brake and hydraulic brake work together, the front axle braking force is distributed to the ISG motor brake and hydraulic brake according to the value of the Ffb and Fisgmax and the rear axle braking force is distributed to the rear-drive motor brake and hydraulic brake according to the value of the Frb and Fmmax. If SOC > 0.85 or the braking strength is higher than 0.7, regenerative braking is not allowed, only hydraulic braking works.

3.4. Strategy Validation

To verify the feasibility of the proposed strategy for plug-in 4WD HEV, a hardware-in-the-loop (HIL) test system was established, as shown in Figure 8a. A real-world vehicle controller was selected for the verification of the proposed strategy. The schematic diagrams of the HIL tests are shown in Figure 8b. It mainly includes three parts: The vehicle controller, test cabinet and host computer. The proposed strategy runs in the vehicle controller. The test cabinet integrates the power management module, fault injection module, signal conditioning module, signal acquisition and control module, system connection and conversion device, industrial computer, etc. It is mainly used to provide signal for the vehicle controller and process the output signal from the vehicle controller. The vehicle’s real-time simulation model and fuel and battery consumption were solved by the host computer.
A CVT-based plug-in 4WD HEV was selected as the study object. The main parameters of this vehicle are shown in Table 1. Eight repeated NEDCs (new European driving cycles) were chosen for the simulation. The initial SOC and SOCcd were 0.95 and 0.3, respectively, and the allowable SOC range was 0.25–0.95.
The HIL test results are shown in Figure 9. The SOC’s change curve, as shown in Figure 9a, indicated that the vehicle ran in the CD mode at the beginning, then, the vehicle ran in the CS mode when SOC was below 0.3. The torque’s change curve is shown in Figure 9b. In the CD mode: Treq was less than Tmmax_w, so only the rear-drive motor provided the demand torque in the driving mode. In the CS mode: The engine provided the main drive torque, the ISG motor gave the assisted drive torque and the rear-drive motor worked only when required drive torque was less than Tel_w. In the brake mode, in most cases, both the rear-drive motor and ISG motor provided the braking torque. Summing up, the proposed strategy could distribute the demand torque between the engine, rear-drive motor and ISG motor according to the already established rules. Therefore, the proposed strategy for the plug-in 4WD HEV was feasible.

4. Mathematical Model of the Multi-Objective Optimization

When establishing the objective function, scholars mainly take fuel consumption, equivalent fuel consumption or use the cost under certain driving cycles as the objective function of the economy index. In fact, PHEVs usually run in the CD and CS stage. In the CD stage, the vehicle is basically driven by the battery, the emission and fuel consumption in this stage are approximately zero, so it is reasonable to select electric energy consumption as the objective function of economy index in this stage; in the CS stage, energy consumption mainly comes from the engine, so fuel consumption should be chosen as the objective function of the economy index in this stage.
The objectives of vehicle economy and power performance usually conflict with each other. Therefore, in the optimization process, if the power performance is not considered, the vehicle economy is only taken as the optimization objective, the power performance will be greatly sacrificed to obtain the optimal economy. The power performance evaluation indices include the vehicle’s maximum speed, acceleration time and the maximum climbing gradient. When dealing with the PHEVs parameter optimization problem, the acceleration time is generally chosen as the evaluation index of the power performance. In this study, the acceleration time of 0–120 km/h is selected as the objective function of the power performance index.
The parameters of powertrain (maximum power of engine, rear-drive motor and ISG motor, speed ratio of the front final drive and the rear final drive) and control strategy parameters (kup and klow) have significant effects on the electric energy consumption of the CD stage, fuel consumption of the CS stage and acceleration time from 0 to 120 km/h. Therefore, this paper chose the above parameters to optimize, as shown in Table 2.
The multi-objective parameter optimization problem is a highly non-linear and non-convex optimization problem due to a large number of variables and limit constraints. This optimization problem minimizes objective functions subjected to a set of inequality constraints. Mathematically, this optimization problem is defined as,
{ M i n i m i z e   J ( X ) = [ f e l e ( X ) , f f u e l ( X ) , f a c c ( X ) ] T S u b j e c t e d   t o   { a j X j b j N e min N e ( X ) N e max T e min T e ( X ) T e max N m min N m ( X ) N m max T m min T m ( X ) T m max N i s g min N i s g ( X ) N i s g max T i s g min T i s g ( X ) T i s g max X = [ P e max , P m max , P i s g max , i f o , i r o , k u p , k l o w ] T
where the object function   J ( X ) is a vector object function assembled by three scalar object functions, which are the 100 km fuel consumption f f u e l ( X ) , 100 km electric energy consumption f e l e ( X ) and acceleration time from 0 to 120 km/h f a c c ( X ) , as illustrated in Equations (10)–(12). X is the decision vector, which includes the maximum power of the engine P e max , maximum power of the rear-drive motor P m max , maximum power of the ISG motor P i s g max , speed ratio of the front final drive i f o , speed ratio of the rear final drive i r o and control strategy parameters ( k u p and k l o w ), X j is the decision variable, a and b are the lower and upper limits of each decision variable, respectively. The value of a and b is given in the “Range” column of Table 2. N i max and N i min are the upper and lower limits of the components’ speed, respectively. T i max and T i min are the upper and lower limits of the components’ output torque, respectively.
In order to obtain the above objective functions’ values, different simulations are needed.
When the electric energy consumption function is simulated, the vehicle runs in the CD stage through the SOC constraint. The 100 km electric energy consumption under the CD stage is calculated by:
{ f e l e ( X ) = 100 3600 S 1 k = 0 N e l e P b ( X , k ) Vehicle   is   running   in   the   CD   stage ,   and   the   SOC   constraint   is   described   as { S O C ( 0 ) = 0.95 0.25 S O C ( k ) 0.95 0.3 S O C ( N e l e ) 0.95
where k represents the kth simulation step, one simulation step was set to one second in this paper. N e l e is the total step for the electric energy consumption simulation, so k = 1 , 2 , , N e l e 1 . P b is the battery power, its unit is kw and S 1 is the driving mileage in the CD stage, its unit is km. S O C ( 0 ) is the initial SOC value. S O C ( k ) is the SOC value in the kth simulation step. S O C ( N e l e ) is the SOC value at the end of the simulation.
When the fuel consumption function is simulated, the vehicle runs in the CS stage through constraining the final SOC to be equal to the initial SOC. The 100 km fuel consumption under the CS stage is calculated by:
{ f f u e l ( X ) = 100 ρ S 2 k = 0 N f u e l m ˙ f ( X , k ) Vehicle   is   running   in   the   CS   stage ,   and   the   SOC   constraint   is   described   as { S O C ( 0 ) = 0.3 0.25 S O C ( k ) 0.95 0.25 S O C ( N f u e l ) 0.35
where m ˙ f is the instantaneous fuel consumption, it can be obtained by looking up the table by the engine speed and engine torque, its unit is g/s. S 2 is the driving mileage in the CS stage. ρ is the density of gasoline. N f u e l is the total steps for the fuel consumption simulation.
When the electric energy consumption and fuel consumption functions are simulated, the required torque and speed at the wheel can be obtained from the given driving cycle, while the acceleration time function is simulated, the required torque at the wheel is determined by the outputs of the engine, rear-drive motor and ISG motor. The speed at the wheel is determined after vehicle’s velocity discretization at the interval of ΔV. Then, the acceleration time from 0 to 120 km/h is regarded as the time accumulation between two velocity steps. The acceleration time function is described as:
f a c c ( X ) = k = 0 N a c c Δ V a ( X , k )
where Δ V is the velocity interval selected as 1 km/h, N a c c are the total velocity steps, its value is 121. a ( X , k ) is the acceleration at velocity step k , it can be calculated by:
{ a ( X , k ) = T w ( X , k ) R F d r a g ( X , k ) δ m T w ( X , k ) = ( T e ( k ) + T i s g ( k ) ) i c v t i f o + T m ( k ) i r o F d r a g ( X , k ) = m g sin α + m g f cos α + C d A V ( k ) 2 21.15
where T w ( X , k ) is the driving torque at the wheel, R is the radius of the wheel and F d r a g ( X , k ) is the sum of rolling resistance, gradient resistance and air resistance.

5. Optimization Algorithm

Compared with single-objective optimization problems, the solution of multi-objective optimization problems is more difficult. In single-objective optimization problems, the optimal solution is unique, while in multi-objective optimization problems, the optimal solution is not unique, there is a set of compromise solutions, called the Pareto solution set. In this paper, the evolutionary based NSGA-II approach was selected to obtain the Pareto solution set of the multi-objective parameter optimization problems.
The flow chart of multi-objective parameter optimization based on the NSGA-II algorithm is shown in Figure 10. Firstly, the population was initialized, and then the vehicle ran in the CD stage, CS stage and 0–120 km/h acceleration stage by setting different operating conditions, and the objective function values in different operating stages were calculated. On the basis of obtaining f e l e , f f u e l and f a c c , the hierarchy and crowding degree of each individual could be obtained by fast non-dominated ranking and crowding degree calculation. Individuals with small hierarchy and high crowding degree were selected to enter the mating pool. After crossover and mutation, the offspring were generated, the offspring and the parent were merged into a unified population, and the elite strategy was used to generate a new generation of population. Then the above steps were repeated until the termination conditions were met.

6. Optimization Results and Analysis

The powertrain’s dynamic model, energy management strategy and calculation model of the objective function were built by using MATLAB/Simulink simulation software, in which the driving cycles were selected under the urban road environment, as shown in Figure 11, the CD stage ran under one FUDS (federal urban driving schedule) driving cycle and the CS stage ran under three repetitive FUDS driving cycles. The parameters of the dynamic model were as follows: The wind resistance coefficient was 0.32 and the windward area was 2.28. The rolling resistance coefficient was 0.0135, the wheel radius was 0.307 and the conversion coefficient of rotary inertia was 1.02. The vehicle mass was calculated by:
m = 1100 + 3 . 2 P e max + 63 P m max / 30 + 58 P i s g max / 30
The parameters of the NSGA-II algorithm were as follows: Population size was 100, the Pareto fraction was 0.2, the maximum generation was 500, the termination algebra was 500, the deviation of fitness function was 1 × 10–100, the crossover rate was 0.8 and the mutation rate was 0.1. After optimization, we got the three-dimensional Pareto solution set as shown in Figure 12a. There are 20 non-dominant solutions. Almost all the optimal solutions distribute on a narrow space surface.
The term Pareto front graph between the acceleration time and electric energy consumption is shown in Figure 12b. From the graph, the achievable extreme electric energy consumption was obtained at non-dominated solution 3 (NDS-3, as shown in Figure 12), its value was 11.63 kwh/100 km. The achievable extreme acceleration time was obtained at non-dominated solution 1 (NDS-1, as shown in Figure 12), its value was 6.5 s. It can be seen that the 100 km electric energy consumption was in conflict with the acceleration time from 0–120 km/h, so it was impossible to achieve the extreme electric energy consumption and the extreme acceleration capacity at the same time. The main reason was that these two goals have different requirements for powertrain parameters. Electric energy consumption is related to the parameters of the maximum power of the rear-drive motor and ISG motor, speed ratio of the front final drive and the rear final drive. In order to obtain the extreme electric energy consumption, the optimal above-mentioned parameters need to be achieved to make the rear-drive motor and ISG motor’s efficiencies highest. However, extreme acceleration capacity was achieved by increasing the maximum power of the engine, the rear-drive motor and ISG motor as much as possible on the condition that all kinds of constraints were satisfied. The comparison of the rear-drive motor and ISG motor’s efficiencies between NDS-1 and NDS-3 is provided in Table 3. The rear-drive motor and ISG motor’s operating efficiencies of NDS-3 was obviously higher than that of NDS-1.
The term Pareto front graph between the acceleration time and fuel consumption is shown in Figure 12c. From the graph, the achievable extreme fuel consumption was obtained at non-dominated solution 4 (NDS-4, as shown in Figure 12), its value was 4.34 L/100 km. It can be seen that the 100 km fuel consumption was also in conflict with the acceleration time from 0–120 km/h. The maximum power of the engine played a major role. The effects of the engine’s maximum power on the fuel consumption and acceleration time are shown in Figure 13. The points in Figure 13 corresponded to the points in the term Pareto front in Figure 12. From Figure 13, fuel consumption increased with the increase of the maximum power of engine (the reason was that the load rate of the engine decreased with the increase of the maximum power of engine). However, the acceleration time decreased with the increase of engine maximum power.
From Figure 12, in order to achieve one extreme goal, solutions of NDS-1, NDS-3 and NDS-4 sacrificed a lot of the other goals. Therefore, a best compromise solution should be found after obtaining the Pareto optimal front. The fuzzy min–max approach in [14] can be applied to find the best compromise solution. The pth objective function is represented by a membership function hp, defined as:
h p = { 1      f p f p min f p max f p f p max f p min f p min < f p < f p max 0      f p f p max
where p = 1 , 2 , , N o b j , N o b j is the number of the objective functions for the member of the Pareto optimal front and f p min and f p max are the minimum and maximum values of the pth objective function of Pareto optimal set. Each member of the Pareto optimal front is evaluated using h, defined as:
h = p = 1 N o b j h p
The best compromise solution is that having the maximum value of h. The solution of NDS-2 is the best compromise solution after using the fuzzy min–max approach. Therefore, the solution of NDS-2 was chosen as the final optimal solution in this paper. The results for NDS-1, NDS-2, NDS-3, NDS-4 and the original scheme (OS) are shown in Table 4. The comparison between the original scheme (OS) and the final optimization scheme (NDS-2) is also shown in Table 4. The original scheme’s powertrain parameters and control strategy parameters were obtained from the prototype vehicle manufactured by an automobile company. The original scheme’s 100 km electric energy consumption, 100 km fuel consumption and acceleration time were achieved from a simulation using the uniform simulation model built in this study. According to this table, compared with the original scheme (OS), the final optimized scheme (NDS-2)’s 100 km electric energy consumption, 100 km fuel consumption and acceleration time were reduced by 1.21%, 6.18% and 5.49%, respectively.
In this study, the Pareto optimal front obtained from the NSGA-II approach was compared with the weight approach. The weight approach is described as follows:
{ J = w 1 J e l e + w 2 J f u e l + w 3 J a c c i = 1 3 w i = 1 w i [ 0 : 0.25 : 1 ] ,   i = 1 , 2 , 3
The weight approach transforms the multi-objective optimization into single-objective optimization by setting weight factors w 1 , w 2 and w 3 . By varying these weight factors, different portions of the Pareto optimal front can be generated [14]. In this study, the varying step of all weight factors was 0.25, then, there were 15 combinations of weight factors for plug-in HEV parameter optimization. These weight factor combinations were computed one by one to obtain a group of optimal solutions. These optimal solutions were the Pareto optimal solutions searching by the weight approach.
The distribution of optimal solutions in the Pareto optimal front obtained using the weight approach is shown in Figure 14. Comparing Figure 14 with Figure 12, it can be seen that the Pareto optimal solutions obtained by the NSGA-II algorithm are well distributed over the entire Pareto optimal front. However, the optimal solutions obtained by the weight approach were non-uniformly distributed over the entire Pareto optimal front. Further, all the non-dominated solutions could not be obtained and some of the solutions obtained were inferior.
After obtaining the Pareto optimal front, the fuzzy min–max approach in [14] was applied to find the best compromise solution. The best compromise solutions obtained by using the NSGA-II algorithm and the weight approach are shown in Table 5. The best compromise solutions obtained by using the NSGA-II algorithm and the weight approach were (12.24 kwh/100 km, 5.01 L/100 km and 8.6 s) and (11.82 kwh/100 km, 5.27 L/100 km and 8.8 s), respectively. This shows that the best compromise solution obtained by using the NSGA-II algorithm was better than the weight approach.

7. Conclusions

The plug-in 4WD HEV’s energy management control strategy based on the CD–CS mode was developed. Then, the multi-objective optimization’s mathematical model, which aims at minimizing the electric energy consumption under the CD stage, the fuel consumption under the CS stage and the acceleration time from 0–120 km/h, was established. Finally, the multi-objective parameter optimization problem was solved using the evolutionary based NSGA-II approach. Some of the results were compared with the original scheme and the classical weight approach.
From the term Pareto front graph obtained using the NSGA-II approach, it could be seen that the two objectives of 100 km electric energy consumption and 100 km fuel consumption were in conflict with the acceleration time from 0–120 km/h, respectively. It was impossible for them to achieve the optimum at the same time. Therefore, the best compromise solution (NDS-2) was chosen as the final optimized scheme. Compared with the original scheme, the final optimized scheme’s 100 km electric energy consumption, 100 km fuel consumption and acceleration time were reduced by 1.21%, 6.18% and 5.49% respectively.
Compared with the weight approach, the Pareto optimal solutions obtained using the NSGA-II approach were better distributed over the entire Pareto optimal front, as well as the best compromise solution was also better.
NSGA-II is an evolutionary based multi-objective optimization algorithm, the main drawback of this approach is computationally burdensome. Therefore, pursuing a more efficient algorithm to solve plug-in 4WD HEV’s multi-objective parameter optimization problems is future research work.

Author Contributions

Y.C. wrote the paper, built the mathematical model of the multi-objective optimization and analyzed the optimization results; Y.Z. built the dynamic model of the powertrain and developed control strategy based on CD-CS mode; J.Y. provided suggestions on mathematical model building and paper writing; Z.W. provided optimization algorithms and made revisions to the manuscript.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (Grant No. 51678075; 51665020), Major Special Science and Technology Project of Hunan Province (Grant No. 2019SK2171), Hunan Province Key Laboratory of Smart Roadway and Cooperative Vehicle-Infrastructure Systems (Grant No. 2017TP1016) and Key Research and Development Project of Jiangxi Province (Grant No. 20192BBEL50015).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations were used in this manuscript:
CDcharge depleting
CScharge sustain
CVTa continuously variable transmission
4WD4-wheel drive
PHEVsplug-in hybrid electric vehicles
NSGA-IINon-dominated Sorting Genetic Algorithms-II
Athe windward area, its unit is m2
a ( X , k ) the acceleration at velocity step k, its unit is m/s2
α ( t ) road gradient
c d the coefficient of air resistance
f r the rolling resistance coefficient
f e l e ( X ) the 100 km electric energy consumption under CD stage, its unit is kwh/100 km
f f u e l ( X ) the 100 km fuel consumption under CS stage, its unit is L/100 km
f a c c ( X ) the acceleration time from 0 to 120 km/h, its unit is s
i gear ratio
k u p control strategy parameter
k l o w control strategy parameter
m ˙ f the instantaneous fuel consumption, its unit is g/s
N e l e the total steps for electric energy consumption simulation
N f u e l the total steps for fuel consumption simulation
N a c c the total velocity steps
P b the battery terminal power, its unit is kw
P e max the maximum power of engine, its unit is kw
P m max the maximum power of rear-drive motor, its unit is kw
P i s g max the maximum power of ISG motor, its unit is kw
Q 0 the nominal battery capacity, its unit is Ah
R the wheel’s radius, its unit is m
R int the battery’s internal resistance, its unit is Ω
S 1 the driving mileage in CD stage, its unit is km
S 2 the driving mileage in CS stage, its unit is km
T e the output torque of the engine, its unit is Nm
T i s g the output torque of ISG motor, its unit is Nm
T m the output torque of the rear-drive motor, its unit is Nm
v ( t ) the vehicle speed, its unit is km/h
V o c the battery’s terminal voltage, its unit is V
ω the angular speed, its unit is rad/s
η i s g the efficiency of ISG motor
η m the efficiency of the rear-drive motor

References

  1. Hu, X.; Martinez, C.M.; Yang, Y. Charging, power management, and battery degradation mitigation in plug-in hybrid electric vehicles: A unified cost-optimal approach. Mech. Syst. Signal Process. 2017, 87, 4–16. [Google Scholar] [CrossRef]
  2. Yao, M.; Qin, D.; Zhou, X.; Zhan, S.; Zeng, Y. Integrated optimal control of transmission ratio and power split ratio for a CVT-based plug-in hybrid electric vehicle. Mech. Mach. Theory 2019, 136, 52–71. [Google Scholar] [CrossRef]
  3. Xie, S.; Li, H.; Xin, Z.; Liu, T.; Wei, L. A pontryagin minimum principle-based adaptive equivalent consumption minimum strategy for a plug-in hybrid electric bus on a fixed route. Energies 2017, 10, 1379. [Google Scholar] [CrossRef]
  4. Liu, Y.; Li, J.; Ye, M.; Qin, D.; Zhang, Y.; Lei, Z. Optimal energy management strategy for a plug-in hybrid electric vehicle based on road grade information. Energies 2017, 10, 412. [Google Scholar] [CrossRef]
  5. Roy, H.K.; McGordon, A.; Jennings, P.A. A generalized powertrain design optimization methodology to reduce fuel economy variability in hybrid electric vehicles. IEEE Trans. Veh. Technol. 2013, 63, 1055–1070. [Google Scholar] [CrossRef]
  6. Wang, Q.; Duan, B.; Wang, P.; Gong, Y.; Zhu, Q. Optimization of powertrain transmission parameters of plug-in hybrid electric vehicle. J. Jilin Univ. Technol. Ed. 2017, 47, 1–7. [Google Scholar]
  7. Zan, S.; Qin, D.; Zeng, Y. A study on the parameter optimization of hybrid electric vehicle based on multiple driving cycles. Automot. Eng. 2016, 38, 922–928. [Google Scholar]
  8. Lei, Z.; Cheng, D.; Liu, Y.; Qin, D.; Zhang, Y.; Xie, Q. A dynamic control strategy for hybrid electric vehicles based on parameter optimization for multiple driving cycles and driving pattern recognition. Energies 2017, 10, 54. [Google Scholar] [CrossRef]
  9. Liu, Y.; Li, J.; Qin, D.; Lei, Z.; Xie, Q.Y. Parameter optimization of hybrid electric vehicle based on multi-cycle optimization algorithm. J. Mech. Eng. 2017, 53, 61–69. [Google Scholar] [CrossRef]
  10. Shahverdi, M.; Mazzola, M.S.; Grice, Q.; Doude, M. Pareto front of energy storage size and series HEV fuel economy using bandwidth-based control strategy. IEEE Trans. Transp. Electrif. 2016, 2, 36–51. [Google Scholar] [CrossRef]
  11. Dimitrova, Z.; Marechal, F. Techno-economic design of hybrid electric vehicles and possibilities of the multi-objective optimization structure. Appl. Energy 2016, 161, 746–759. [Google Scholar] [CrossRef]
  12. Ribau, J.P.; Silva, C.M.; Sousa, J.M.C. Efficiency, cost and life cycle CO2 optimization of fuel cell hybrid and plug-in hybrid urban buses. Appl. Energy 2014, 129, 320–335. [Google Scholar] [CrossRef]
  13. Zeng, Y.; Qin, D.; Su, L.; Yao, M. Cost, fuel consumption and emission multi-objective parameter optimization for the powertrain of a Plug-in HEV. Automot. Eng. 2016, 38, 397–402. [Google Scholar]
  14. Huang, K.; Wang, Q.; Qiu, M.; Zhang, Y.; Liu, Z.; Geng, L. Parameter optimization of multi-mode hybrid vehicle considering mode-switching frequency. J. Xi’an Jiaotong Univ. 2019, 53, 1–8. [Google Scholar]
  15. Song, C.; Zhou, F.; Xiao, F.; Chang, C.; Shao, Y. Parameter matching of on-board hybrid energy storage system based on convex optimization method. J. Mech. Eng. 2017, 53, 44–51. [Google Scholar] [CrossRef]
  16. Reddy, S.S.; Bijwe, P.R. Multi-objective optimal power flow using efficient evolutionary algorithm. Int. J. Emerg. Electr. Power Syst. 2017, 18, 1–21. [Google Scholar] [CrossRef]
  17. Roghanian, E.; Sadjadi, S.J.; Aryanezhad, M.B. A probabilistic bi-level linear multi-objective programming problem to supply chain planning. Appl. Math. Comput. 2007, 188, 786–800. [Google Scholar] [CrossRef]
  18. Hsu, G.J.; Leung, P.S.; Ching, C.T. Energy planning in taiwan: An alternative approach using a multiobjective programming and input-output model. Energy J. 1988, 9, 53–72. [Google Scholar] [CrossRef]
  19. Osyczka, A. An approach to multicriterion optimization problems for engineering design. Comput. Methods Appl. Mech. Eng. 1978, 15, 309–333. [Google Scholar] [CrossRef]
  20. Steuer, R.E. Multiple criteria optimization: Theory, computation and application. J. Oper. Res. Soc. 1988, 39, 879. [Google Scholar]
  21. He, Y.; Mei, J.; Zang, J.; Xie, S.; Zhang, F. Multicriteria optimization design for end effector mounting bracket of a high speed and heavy load palletizing robot. Math. Probl. Eng. 2018, 2018, 6049635. [Google Scholar] [CrossRef]
  22. Zhang, X.; Song, X.; Qiu, W.; You, Y.; Deng, N. Multi-objective optimization of tension leg platform using evolutionary algorithm based on surrogate model. Ocean Eng. 2018, 148, 612–631. [Google Scholar] [CrossRef]
  23. Reddy, S.S.; Abhyankar, A.R.; Bijwe, P.R. Reactive power price clearing using multi-objective optimization. Energy 2011, 36, 3579–3589. [Google Scholar] [CrossRef]
  24. Reddy, S.S.; Bijwe, P.R. Differential evolution-based efficient multi-objective optimal power flow. Neural Comput. Appl. 2019, 31, 509–522. [Google Scholar] [CrossRef]
  25. Reddy, S.S. Multi-objective based congestion management using generation rescheduling and load shedding. IEEE Trans. Power Syst. 2017, 32, 852–863. [Google Scholar]
  26. Zeng, Y.; Cai, Y.; Kou, G.; Gao, W.; Qin, D. Energy management for plug-in hybrid electric vehicle based on adaptive Simplified-ECMS. Sustainability 2018, 10, 2060. [Google Scholar] [CrossRef]
  27. Zeng, Y.; Sheng, J.; Li, M. Adaptive real-time energy management strategy for plug-in hybrid electric vehicle based on Simplified-ECMS and a novel driving pattern recognition method. Math. Probl. Eng. 2018, 2018, 5816861. [Google Scholar] [CrossRef]
  28. Xie, S.; Hu, X.; Lang, K.; Qi, S.; Liu, T. Powering mode-integrated energy management strategy for a plug-in hybrid electric truck with an automatic mechanical transmission based on pontryagin’s minimum principle. Sustainability 2018, 10, 3758. [Google Scholar] [CrossRef]
  29. Zeng, Y.; Cai, Y.; Chu, C.; Kou, G.; Gao, W. Integrated Energy and Catalyst Thermal Management for Plug-In Hybrid Electric Vehicles. Energies 2018, 11, 1761. [Google Scholar] [CrossRef]
  30. Liu, J.; Chen, Y.; Zhan, J.; Shang, F. Heuristic Dynamic Programming Based Online Energy Management Strategy for Plug-In Hybrid Electric Vehicles. IEEE Trans. Veh. Technol. 2019, 68, 4479–4493. [Google Scholar] [CrossRef]
  31. Xiao, R.; Liu, B.; Shen, J.; Guo, N.; Yan, W.; Chen, Z. Comparisons of energy management methods for a parallel plug-in hybrid electric vehicle between the convex optimization and dynamic programming. Appl. Sci. 2018, 8, 218. [Google Scholar] [CrossRef]
  32. Moura, S.J.; Fathy, H.K.; Callaway, D.S.; Stein, J.L. A stochastic optimal control approach for power management in plug-In hybrid electric vehicles. IEEE Trans. Control Syst. Technol. 2011, 19, 545–555. [Google Scholar] [CrossRef]
  33. Fletcher, T.; Thring, R.; Watkinson, M. An Energy Management Strategy to concurrently optimise fuel consumption & PEM fuel cell lifetime in a hybrid vehicle. Int. J. Hydrogen Energy 2016, 41, 21503–21515. [Google Scholar] [Green Version]
  34. Wu, X.; Hu, X.; Yin, X.; Li, L.; Zeng, X.; Volker, P. Convex programming energy management and components sizing of a plug-in fuel cell urban logistics vehicle. J. Power Sources 2019, 423, 358–366. [Google Scholar] [CrossRef]
  35. Liu, Y.; Li, J.; Chen, Z.; Qin, D.; Zhang, Y. Research on a multi-objective hierarchical prediction energy management strategy for range extended fuel cell vehicles. J. Power Sources 2019, 429, 55–66. [Google Scholar] [CrossRef]
  36. Wang, H.; Huang, Y.; Khajepour, A.; Song, Q. Model predictive control-based energy management strategy for a series hybrid electric tracked vehicle. Appl. Energy 2016, 182, 105–114. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, H.; Li, X.; Wang, W.; Han, L.; Xiang, X. Markov velocity predictor and radial basis function neural network-based real-time energy management strategy for plug-in hybrid electric vehicles. Energy 2018, 152, 427–444. [Google Scholar] [CrossRef]
  38. Zeng, X.; Wang, J. A parallel hybrid electric vehicle energy management strategy using stochastic model predictive control with road grade preview. IEEE Trans. Control Syst. Technol. 2015, 23, 1–15. [Google Scholar] [CrossRef]
  39. Yang, Y.; Hu, X.; Pei, H.; Peng, Y. Comparison of power-split and parallel hybrid powertrain architectures with a single electric machine: Dynamic programming approach. Appl. Energy 2016, 168, 683–690. [Google Scholar] [CrossRef]
  40. Cipek, M.; Pavkovi, D.; Petri, J. A control-oriented simulation model of a power-split hybrid electric vehicle. Appl. Energy 2013, 101, 121–133. [Google Scholar] [CrossRef]
  41. Zhou, X.; Qin, D.; Hu, J. Multi-objective optimization design and performance evaluation for plug-in hybrid electric vehicle powertrains. Appl. Energy 2017, 208, 1608–1625. [Google Scholar] [CrossRef]
Figure 1. The plug-in 4-wheel drive (4WD) hybrid electric vehicle powertrain.
Figure 1. The plug-in 4-wheel drive (4WD) hybrid electric vehicle powertrain.
Applsci 09 04068 g001
Figure 2. (a) The base integrated starter and generator (ISG) motor efficiency map and (b) the base rear-drive motor efficiency map.
Figure 2. (a) The base integrated starter and generator (ISG) motor efficiency map and (b) the base rear-drive motor efficiency map.
Applsci 09 04068 g002
Figure 3. (a) The battery’s terminal voltage and internal resistance and (b) the engine’s fuel consumption.
Figure 3. (a) The battery’s terminal voltage and internal resistance and (b) the engine’s fuel consumption.
Applsci 09 04068 g003
Figure 4. Operating mode switching strategy.
Figure 4. Operating mode switching strategy.
Applsci 09 04068 g004
Figure 5. Control strategy in charge depleting (CD) mode.
Figure 5. Control strategy in charge depleting (CD) mode.
Applsci 09 04068 g005
Figure 6. Control strategy in charge sustain (CS) mode.
Figure 6. Control strategy in charge sustain (CS) mode.
Applsci 09 04068 g006
Figure 7. The braking strategy.
Figure 7. The braking strategy.
Applsci 09 04068 g007
Figure 8. (a) The hardware-in-the-loop test system. (b) Schematic diagrams of the hardware-in-the-loop test system.
Figure 8. (a) The hardware-in-the-loop test system. (b) Schematic diagrams of the hardware-in-the-loop test system.
Applsci 09 04068 g008
Figure 9. The hardware-in-the-loop (HIL) test results. (a) The state of charge’s (SOC’s) change curve. (b) The torque’s change curve.
Figure 9. The hardware-in-the-loop (HIL) test results. (a) The state of charge’s (SOC’s) change curve. (b) The torque’s change curve.
Applsci 09 04068 g009
Figure 10. The flow chart of multi-objective parameter optimization based on the non-dominated sorting genetic algorithms-II (NSGA-II) algorithm.
Figure 10. The flow chart of multi-objective parameter optimization based on the non-dominated sorting genetic algorithms-II (NSGA-II) algorithm.
Applsci 09 04068 g010
Figure 11. FUDS (federal urban driving schedule) driving cycle.
Figure 11. FUDS (federal urban driving schedule) driving cycle.
Applsci 09 04068 g011
Figure 12. The optimal solutions. (a) The three-dimensional Pareto solution set; (b) the term Pareto front graph between the acceleration time and electric energy consumption and (c) the term Pareto front graph between the acceleration time and fuel economy.
Figure 12. The optimal solutions. (a) The three-dimensional Pareto solution set; (b) the term Pareto front graph between the acceleration time and electric energy consumption and (c) the term Pareto front graph between the acceleration time and fuel economy.
Applsci 09 04068 g012
Figure 13. Effects of the engine’s maximum power on the fuel consumption and acceleration time.
Figure 13. Effects of the engine’s maximum power on the fuel consumption and acceleration time.
Applsci 09 04068 g013
Figure 14. The Pareto optimal front using the weight approach. (a) The term Pareto front graph between the acceleration time and electric energy consumption and (b) the term Pareto front graph between the acceleration time and fuel economy.
Figure 14. The Pareto optimal front using the weight approach. (a) The term Pareto front graph between the acceleration time and electric energy consumption and (b) the term Pareto front graph between the acceleration time and fuel economy.
Applsci 09 04068 g014
Table 1. Main parameters of the plug-in 4WD hybrid electric vehicle (HEV).
Table 1. Main parameters of the plug-in 4WD hybrid electric vehicle (HEV).
NameValue
Vehicle mass1590 kg
Windward area2.265 m2
Tire radius0.307 m
Coefficient of air resistance0.0135
Max power of the engine93 kW
Max power of the rear-drive motor55 kW
Max power of the ISG motor30 kW
Rated capacity of the battery30 Ah
Table 2. The powertrain and control strategy’s parameters.
Table 2. The powertrain and control strategy’s parameters.
ParametersDescriptionsRange
Pe max/kwEngine’s maximum power[63, 120]
Pm max/kwRear-drive motor’s maximum power[65, 110]
Pisg max/kwISG motor’s maximum power[20, 60]
ifoSpeed ratio of the front final drive[3.81, 6.92]
iroSpeed ratio of the rear final drive[4.4, 8.68]
kupControl strategy parameter[0.1, 1]
klowControl strategy parameter[0.1, 1]
Table 3. Comparison of the rear-drive motor and ISG motor’s efficiencies between NDS-1 and NDS-3.
Table 3. Comparison of the rear-drive motor and ISG motor’s efficiencies between NDS-1 and NDS-3.
The Rear-Drive MotorISG Motor
0.7 ≤ ηmot < 0.8ηmot ≥ 0.80.7 ≤ ηisg < 0.8ηisg ≥ 0.8
NDS-10.740.260.540.46
NDS-30.630.370.370.63
Table 4. The results for NDS-1, NDS-2, NDS-3, NDS-4 and the original scheme (OS).
Table 4. The results for NDS-1, NDS-2, NDS-3, NDS-4 and the original scheme (OS).
PemaxPmmaxPisg maxifoirokupklowfeleffuelfacc
NDS-1118102496.327.520.690.3313.986.946.5
NDS-36666226.97.240.70.9711.635.439.1
NDS-46567213.875.950.120.7912.174.3310.2
ContrastOS9355355.186.680.540.7312.395.349.1
NDS-27273305.267.420.710.7812.245.018.6
Reduce by--------------−1.21%−6.18%−5.49%
Table 5. The best compromise solution obtained by using the NSGA-II algorithm and the weight approach.
Table 5. The best compromise solution obtained by using the NSGA-II algorithm and the weight approach.
PemaxPmmaxPisg maxifoirokupklowfeleffuelfacc
NSGA-II approach7273305.267.420.710.7812.245.018.6
Weight approach8760255.897.060.510.5311.825.278.8

Share and Cite

MDPI and ACS Style

Wang, Z.; Cai, Y.; Zeng, Y.; Yu, J. Multi-Objective Optimization for Plug-In 4WD Hybrid Electric Vehicle Powertrain. Appl. Sci. 2019, 9, 4068. https://doi.org/10.3390/app9194068

AMA Style

Wang Z, Cai Y, Zeng Y, Yu J. Multi-Objective Optimization for Plug-In 4WD Hybrid Electric Vehicle Powertrain. Applied Sciences. 2019; 9(19):4068. https://doi.org/10.3390/app9194068

Chicago/Turabian Style

Wang, Zhengwu, Yang Cai, Yuping Zeng, and Jie Yu. 2019. "Multi-Objective Optimization for Plug-In 4WD Hybrid Electric Vehicle Powertrain" Applied Sciences 9, no. 19: 4068. https://doi.org/10.3390/app9194068

APA Style

Wang, Z., Cai, Y., Zeng, Y., & Yu, J. (2019). Multi-Objective Optimization for Plug-In 4WD Hybrid Electric Vehicle Powertrain. Applied Sciences, 9(19), 4068. https://doi.org/10.3390/app9194068

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop