A chaotic particle swarm optimization algorithm for solving optimal power system problem of electric vehicle

Transcription

1 Research Article A chaotic particle swarm optimization algorithm for solving optimal power system problem of electric vehicle Advances in Mechanical Engineering 2019, Vol. 11(3) 1 9 Ó The Author(s) 2019 DOI: / journals.sagepub.com/home/ade Tianjun Zhu 1,2, Hongyan Zheng 1 and Zonghao Ma 3 Abstract Transportation of electrification has become a hot issue in recent decades and the large-scale deployment of electric vehicles has yet to be actualized. This article proposes a powertrain parameter optimization design approach based on chaotic particle swarm optimization algorithm. To improve the driving and economy performance of pure electric vehicles, chaotic particle swarm optimization algorithm is adopted in this study to optimize principal parameters of vehicle power system. Vehicle dynamic performance simulations were carried out in the Cruise software, and the simulation results before and after optimization were compared. Simulation results show that optimized vehicles by chaotic particle swarm optimization can meet the expected dynamic performance and the driving range has been greatly improved. Meanwhile, it is also viable that the parameters of the optimal objective function can achieve the purpose of balancing the driving performance and economic performance, which provides a reference for the development of vehicle dynamic performance. Keywords Chaotic particle swarm optimization, power system optimization, driving range, electric vehicle Date received: 8 September 2018; accepted: 1 February 2019 Handling Editor: Dumitru Baleanu Introduction Energy crisis is an urgent issue nowadays, and more and more attention is paid to new energy resources. In automobile industry researchers and technicians focus on developing pure electric vehicles (EVs) of high efficiency, energy saving, and zero emissions. But the pure EV s core technology has not yet reached the top level; its cost and performance cannot fully meet the needs of the market. The difficulty in achieving the balance in dynamic and economy has become a bottleneck restricting the further development of pure EVs. In recent years, many scholars and research institutions have done a lot of research work in the optimization algorithm. Among them, various intelligent algorithms are the most prominent, such as genetic algorithm, evolutionary algorithm, direct algorithm, simulated annealing algorithm, multidisciplinary optimization algorithm, and various combinations of algorithms. 1 8 Fan 9 used the genetic algorithm to optimize the size and battery management logic of the hybrid powertrain and obtained reasonable component size and parameter matching solution. Wu et al. 10 adopted multi-objective evolutionary algorithm to optimize the 1 Department of Mechanical and Automotive Engineering, Zhaoqing University, Zhaoqing, China 2 College of Mechanical and Equipment Engineering, Hebei University of Engineering, Handan, China 3 China Automotive Technology & Research Center, Tianjin, China Corresponding author: Tianjun Zhu, College of Mechanical and Equipment Engineering, Hebei University of Engineering, 199 South Street, Handan City , China. happy.adam2012@hotmail.com Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).

2 2 Advances in Mechanical Engineering hybrid power, and optimized results are more conducive to battery life and improve power efficiency. S Chen 11 used particle swarm optimization (PSO) algorithm for pure EV transmission system; the results obtained proved that the vehicle economy and power have been effectively improved. Zhang et al. 12 adopted evolutionary algorithm to optimize the hybrid power system control strategy globally of hybrid electric city bus; the parameters of the control strategy can be determined accurately, which can improve the fuel economy and pollutant emission. Ant colony genetic algorithm is used to optimize the dynamic parameters of pure EVs. 13 The optimization results show that it is feasible to apply the ant colony algorithm to match the parameters of pure EV power system. Zhou et al. 14 proposed a hybrid electric vehicle (HEV) sizing method using Chaos-enhanced acceleration particle swarm optimization (CAPSO) algorithm to achieve the global optimization. The paper focuses on the component sizing optimization of HEV powertrain system, but does not deal with the optimization of HEV or EV economy and power performance. Jiang et al. 15 applied an improved chaotic particle swarm optimization (ICPSO) algorithm to optimize the key parameters of energy management strategy. The paper improves vehicle economic performance; the vehicle dynamic performance is not involved. Dong et al. 16 proposed a chaotic particle swarm optimization (CPSO) algorithm to solve the route optimization problem of parking lots. Pehlivanoğlu 17 applied a vibrational mutation PSO combining the idea of mutation strategy to solve path planning problem of unmanned aerial vehicle. From the literature review above, it is still an open issue to effectively make the powertrain components operate at high efficiency so that dynamic performance and the driving range have been greatly improved, especially in principal parameters of powertrain multiobjective optimization for pure EV. This article proposes a powertrain parameters design approach based on CPSO algorithm. The principal parameters of powertrain (including peak power of motor, battery s capacity and final drive ratio) are set as the optimized variables by CPSO, with the vehicle performance index (maximum speed, acceleration time, and maximum grade) being defined as the constraint conditions. The rest of this article is organized as follows. Section Establishment of vehicle model presents the vehicle models in detail. The CPSO algorithm for powertrain parameters is proposed in section CPSO algorithm. The optimization results and discussion are given in section Analysis of results. Finally, the conclusions and future work are given in section Conclusion. Establishment of vehicle model Whole vehicle model The whole vehicle model contains a wide range of modules such as battery, gearbox, motor, brake, tire, and MATLAB DLL module is related to the economy of EVs and include strategy of regenerative braking. It can establish a simple vehicle simulation model by selecting the appropriate parts from the model library and dragging them to the modeling panel. The parameterization of pure EV mainly includes that of whole vehicle, electromotor, and battery components. The EV vehicle model is shown in Figure 1. The EV vehicle module contains the basic information of the vehicle including vehicle curb weight, vehicle gross weight, wheel base, height of center of gravity, windward area, drag coefficient of air resistance, and so on. The important parameters of the vehicle module are shown in Table 1. Electromotor module Electric motor is one of the core components of EVs. The accuracy of the electric motor model has significant influence on the precision of EVs simulation results. The basic parameters of the electromotor module include the rated voltage, electromotor speed, electromotor inertia, maximum current, electromotor thermal characteristics, and so on. It is also necessary to input the external characteristic curve (Figure 2) of the electromotor and the map diagram of the electromotor efficiency (Figure 3). Battery module Battery, as energy storage device of EV, plays an important role in the development of EVs, while its electrochemical process is a complex nonlinear process, which is sensitive to environment. Lithium-ion battery, which has the advantages of environmental protection, good performance, high specific energy, and high specific power and is very popular in society, is used in this research. Taking the internal resistance into consideration, this research uses a simple battery model to calculate the current and finally get the battery charge and discharge power. The ideal battery power P batt, des(t) is expressed as follows P batt, des (t)=i 2 batt(t)r int + V batt, oc (t)i batt (t) ð1þ where R int is the battery internal resistance; V batt, oc (t) is open-circuit voltage. Equation (1) is valid for both charge and discharge modes.

3 Zhu et al. 3 Figure 1. EV whole vehicle model. Table 1. Vehicle module parameters. Parameters Value 0.94 motor efficiency map Vehicle curb weight (kg) 1150 Vehicle gross weight (kg) 1600 Length 3 width 3 height (mm) Wheel base (mm) 2600 Frontal area (m 2 ) 1.94 Drag coefficient 0.34 motor efficiency[%] Speed[r/min] Torque[Nm] Figure 3. The map diagram of electromotor efficiency. Figure 2. The external characteristic curve of electromotor. The important input parameters of the battery model are the rated capacity of the battery, the nominal voltage, the highest and lowest voltage, the series-parallel connection of voltages, and the relationship between the cell voltage and the State of Charge (SOC) at normal atmospheric temperature. The parameters in this curve are put into the battery model accurately, as shown in Figure 4.

4 4 Advances in Mechanical Engineering small braking energy. When the speed is greater than or equal to 120 km/h, the energy recovery is not set because of the fast speed and safety. When the vehicle speed is from 10 to 120 km/h and SOC is less than 95%, the motor torque and power are always set as the maximum output to maximize energy recovery. Figure 4. The relationship between charging and discharging battery voltage and SOC. MATLAB DLL control module The MATLAB DLL control module contains the regenerative braking control strategy, which is closely related to the vehicle economy. The control strategy of the module is shown in Figure 5. As shown in Figure 5, the regenerative braking control strategy of the module receives data about brake pedal and vehicle speed. And then the motor load, the front brake pedal pressure and the rear brake pedal pressure are calculated through a series of algorithms. Through a series of algorithms, the motor load, front brake pedal pressure, and rear brake pedal pressure are obtained. When the speed is less than or equal to 10 km/h, the energy recovery is not set because of the CPSO algorithm The optimization of EV power system aims to improve the economy and reduce vehicle cost, while ensuring the dynamic performance. In fact, it basically is a multiobjective optimization issue with constraints. The traditional PSO algorithm brings about many drawbacks in dealing with multi-objective optimization problems with constraints. In this research, it combines chaos algorithm with particle warm optimization algorithm to develop a new optimization algorithm CPSO. In order to prevent some particles from stagnating in the iterative process and fall into the local non-optimal region, a chaotic variable S is generated at the same time as the optimization variable X, and the optimization variable is transformed into a chaotic variable by the carrier, so that the optimization variable is different from the previous state, the results are calculated in the anti-chaotic way to obtain the optimal results, thus solving the basic particle swarm algorithm premature convergence and other issues. The specific steps of CPSO algorithm are as follows: 1. Determine the particle size, the learning factor, the number of iterations and the weighting Figure 5. Regenerative braking control strategy.

5 Zhu et al. 5 factor, and initialize the velocity and position of each particle randomly. 2. The velocity and position of the particles in equation (2) are updated according to the basic particle swarm algorithm Vi k +1 = wv k i c 1 r 1 (x k pbest, i xk i ) + c 2 r 2 (x k gbest xk i ) x k i + 1 = x k i + n k i + 1 ð2þ 3. Chaotization in the optimal location (x gbest )of particle swarm. Map x gbest to the domain [0, 1] of the logistic equation by equation (3) y k 1 = xk gbest Rk min R k max Rk min ð3þ Get chaotic sequence (y k =(y k 1, yk 2,..., yk M )) after M times of iteration for y k 1 by logistic equation (y k n + 1 = myk n (1 yk n )). Map the chaotic sequences back to the original solution space using the following equation x k gbest, m =Rk min +(Rk max Rk min )yk m m=1, 2,..., M It produces a feasible sequence of chaotic variables: x k gbest =(xk gbest, 1, xk gbest, 2,..., xk gbest, M ) Calculate the fitness value of vector in each feasible solution and preserve feasible solution vector, written for Xg k, in accordance with the global best fitness value. 4. Select a particle from the current particle swarm randomly and replace the position vector of the selected particle with the position vector of x k g. 5. Jump to Step 2, until the calculation is up to the maximum iterations number or satisfactory solution. Because in the traditional PSO algorithm, the initial and evolutionary stages of particles are random, which leads to aimless updating of x gbest and x pbest, and accelerate the convergence speed of particles in the evolutionary stage. The CPSO algorithm based on chaotic optimization can effectively avoid local extremum and improve the convergence accuracy and speed of the algorithm. The flow chart of the chaotic particle swarm algorithm is shown in Figure 6. As shown in Figure 6, the key idea of CPSO algorithm is shown in two aspects. First, the initial position and velocity of particles are reset by chaotic sequence. After chaotic processing, the diversity of particle swarm traversal is optimized, while the original stochastic Velocity and position of initial particles Calculate the fitness of each particle If the current value of a particle is greater than the global extreme value, the current value is the individual extreme value Termination conditions judgment Yes The end characteristics of traditional PSO algorithm are preserved. On the basis of the initial population, the initial population is selected from among them. Second, a chaotic sequence is formed based on the optimal position of all the current population to the end of the search, and the best position individual in the chaotic sequence is used to replace the position of one individual in the current population. Because there are many adjacent points in the local optimal solution in each rotation, the search algorithm of chaotic sequence is introduced. Thus, the inert individual can get rid of the local minimum and find the optimal solution quickly. Optimization goals Acceleration demand rather than maximum speed and climbing capability determines the peak power of the motor, so the peak power of the motor is directly related to the vehicle acceleration performance. Power battery capacity directly determines the size of the tram economy. The selection of transmission ratio of main reducer is also closely related to the vehicle performance of EVs. Therefore, the above three parameters are chosen as the optimization variables, expressed as x 3 i = ½p e max, C, i 0 Š Next Generation Particles Chaotic optimization of the optimal position g best is carried out to obtain the optimal * solution g best, and the particle position in the current population is updated with the optimal solution. Update Position Update speed Figure 6. The flow chart of the chaotic particle swarm algorithm. No ð4þ where p e max is the peak power of motor, C is battery s capacity, and i 0 is final drive ratio.

6 6 Advances in Mechanical Engineering Optimization goals Driving range. According to international standards, the driving range under the New European Driving Cycle (NEDC) condition and constant 60 km/h condition are important indicators of vehicle economy. The driving range in the NEDC is calculated by integrating the velocity function at each moment F 1 (X )=S j = Xn i = 1 t j ð + 1 t j V j (t)dt, j = 1, 2,..., N ð5þ where N is the number of segments; V j (t) is the mathematical expression of the velocity, and n is the number of the driving cycles. Energy consumption. Energy consumption is evaluated through the simulations of two driving conditions: constant 60 km/h condition and NEDC condition. As the NEDC condition is not a constant, instead it is composed of four urban driving conditions and one suburb driving condition P = v mgf + C DAv h t 21:15 ð6þ S = Cv ð7þ P where S is the driving range; P is power of constant speed driving; h t is mechanical efficiency; m is vehicle curb weight; g is gravitational acceleration, here g = 9.8 m/s 2 ; f is coefficient of rolling resistance. C D is drag coefficient; A is frontal area; v is constant cruising speed of 60 km/h. In equation (6), the portion in the brackets indicates energy consumption of vehicle resistance. The first part is the energy consumption of rolling resistance and the second part is the energy consumption of air resistance ð ð w N j = F t v(t)dt = t t T tq i 0 h t r wheel v(t)dt ð8þ where F t is driving force; v(t) is function of velocity varying with time; T tq is torque of motor; r wheel is wheel radius; w N j energy consumption under the j working conditions. The energy consumption of NEDC is as follows w N = Xn j = 1 s j w N j ð9þ where s j is energy consumption correction coefficient. Considering the influencing factors of opening/closing air conditioning of vehicles under high and low temperature conditions, we use energy consumption correction coefficient to update the real energy consumption. The influence of air conditioning on vehicle energy consumption is at least 5% 15%. Establishment of constraints In the establishment of the constraints, it is necessary to take the dynamic performance and the ratio of the final drive into account, and also consider the battery energy requirements decided by driving mileage under constant 60 km/h condition and the battery discharge power under the 55% SOC. 1. Constraint conditions under the maximum velocity g 1 (x)= 3600h tp max mgf + C DAv 2 max ø 0 v max 21:15 ð10þ 2. Constraint conditions under the acceleration time g 2 (x)=t g m ð v b " ð v b 0 1 P e max n b mgf + C DAv 2 21:15 # 1 P e max n b mgf + C DAv 2 21:15 dv ø 0 ð11þ where t g is acceleration time;v b is velocity corresponding to the base speed of the motor. 3. Constraint conditions under the maximum climbing degree g 3 (X )= 30P emaxi 0 h t pn b mg cos a max +mg sin a max + C DAv 2 ø 0 21:15 ð12þ where a max is maximum slope angle of road; n b is the base speed of the motor. i 0 is the final drive ratio. 4. Constraint conditions under the driving range of single full charge g 4 (X)= Ch DODh mgf + C DAv 2 21: ø 0 ð13þ where h DOD is discharge efficiency. h is efficiency of motor. 5. Constraint conditions under the power of the electromotor and the maximum power demand of the electromotor

7 Zhu et al. 7 Table 2. Comparison of vehicle parameters before and after optimization. Table 3. Comparison of vehicle performance parameters before and after optimization. Parameter Before optimization After optimization Performance parameter Before optimization After optimization Maximum power of electromotor (kw) Powertrain ratio Battery capacity (Ah) g 5 (X)=k U 2 r soc P emax ø 0 ð14þ where r soc is the internal resistance of MnNiCo lithium battery (0.94 mo). U is battery voltage. 6. Constraint conditions under the powertrain parameters g 6 (x)=i 0 0:377 v maxr wheel n max ø 0 ð15þ Maximum velocity (km/h) Time of full load acceleration km/h (t) Time of full load acceleration km/h (t) Maximum 27% 29% climbing degree Mileage in NEDC (km) Rate of energy consumption in NEDC (kwh/100 km) Mileage in 60 km/h (km) Rate of energy consumption in 60 km/h (kwh/100 km) NEDC: new european driving cycle. where r wheel is rolling radius;n max is peak speed of motor. Generally, the penalty function method is used to handle constraints in optimization problems, it can convert constrained optimization problem to unconstrained optimization problems. Therefore, the new objective function is expressed in equation (16) f i (X)= X j Xn (m) d j F j (X)+ r k = 1 min (0, g k (x)); i=1, 2,..., p ð16þ where n is number of constraint condition; r (m) is penalty factor, and r (m).1; with the iteration of the optimization procedure, r (m + 1) = Cr (m), C. 1. d j is weight coefficient. A relative optimal solution of f i (x) can be figured out when an appropriate m is chosen. If the penalty factor is too big or too small, the accurate value of the new objective function will be affected. Therefore, it is necessary to select the appropriate penalty factor through many times of adjustment and verification. Analysis of results In this article, the scale of the particle swarm is 40, the learning factor is 2.05, and number of iterations is about The optimal variables are peak power of motor, power battery capacity, and final drive ratio in turn. The objective of optimization is to minimize energy consumption and better driving performance. Figure 7. Driving mileage under NEDC condition before the optimization. CPSO algorithm is developed in MATLAB. After a series of calculations, the optimized parameters are compared with the original ones, as shown in Table 2. The optimized parameters were put into the CRUISE for simulation. The simulation results after optimization are compared with the original ones as shown in Table 3. The driving mileage under NEDC condition before and after the optimization are shown in Figure 7 and Figure 8 respectively. The comparison results under NEDC condition show that the driving mileage of EV after the optimization is significantly improved from to km. The driving mileage under constant 60 km/h condition before and after the optimization are shown in Figure 9 and Figure 10 respectively.

8 8 Advances in Mechanical Engineering Figure 8. Driving mileage under NEDC condition after the optimization. Figure 10. Driving mileage under constant 60 km/h condition after the optimization. Figure 9. Driving mileage under constant 60 km/h condition before the optimization. The comparison of simulation results before and after optimization under constant 60 km/h shows that the driving mileage of EV after the optimization is significantly improved from km to km. It is similar to the simulation results under NEDC condition, the effect of economic performance optimization is obvious. To summarize, according to above optimization results, the maximum velocity has decreased but acceleration time reduced slightly, and maximum climbing degree has increased. The driving mileage with NEDC condition of EV increases to km, and the mileage with constant 60 km/h condition of EV increases to km. With the increase of battery capacity, there is a rise in the vehicle weight and then followed by a small climb in the rate of energy consumption. It is feasible to optimize the parameters of power system for pure EV by CPSO, which extends the application range of particle swarm algorithm. Conclusion In this article, the parameter optimization of power system for pure EV was investigated using CPSO method in the MATLAB and CRUISE software. Through establishing the objective function of the longitudinal vehicle dynamics, the CPSO algorithm was developed to optimize the parameter of power system. Some performance simulations were carried out and the simulation results before and after optimization were compared. It can be concluded that it is feasible to apply the chaotic particle swarm to the optimization of the vehicle longitudinal dynamics target. Meanwhile, it is also viable that the parameters of the optimal objective function can achieve the purpose of balancing the power performance and economic performance, which provides a reference for the development of vehicle power performance. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors would like to acknowledge the financial support of Collaborative Innovation and Platform Environment Construction Project of Guangdong Province (2015A ) and Hebei provincial High Level Talents Foundation of Hebei Province (A ) and Hebei Province Key R & D Project ( D). The authors also would like to acknowledge the financial support of Natural Science Foundation of Characteristic Innovation project of Guangdong Ordinary University and Science and

9 Zhu et al. 9 Technology Research Project of Hebei Higher Education (ZD ). The authors also would like to acknowledge the financial support of subsidized projects for overseas students (CL201705). ORCID id Tianjun Zhu References 1. Saini V, Singh S, Nv S, et al. Genetic algorithm based gear shift optimization for electric vehicles. SAE Int J Alternat Powertr 2016; 5: Chakraborty D, Vaz W and Nandi A. Optimal driving during electric vehicle acceleration using evolutionary algorithms. Appl Soft Comput 2015; 34: Hao J, Yu Z, Zhi Z, et al. Optimization of key parameters of energy management strategy for hybrid electric vehicle using direct algorithm. Energies 2016; 9: Li S and Yang G. Evolutionary algorithm and hybrid power system optimization. Beijing, China: Machinery Industry Press, Chen X. Parameter matching and optimization of hybrid SUV. PhD Dissertation, Wuhan University of Technology, Wuhan, China, Liu L. Matching of powertrain parameters of a pure electric vehicle. J HeFei Univ Technol 2007; 5: Yang C. Dynamic modeling and simulation of electric vehicle. PhD Dissertation, Wuhan University of Technology, Wuhan, China, Chen Z and Xiong R. Particle swarm optimization-based optimal power management of plug-in hybrid electric vehicles considering uncertain driving conditions. J Energy 2016; 96: Fan BS. Multidisciplinary optimization of hybrid electric vehicle: component sizing and power management logic. PhD Dissertation, University of Waterloo, Waterloo, ON, Canada, Wu L, Wang Y, Yuan X, et al. Multi-objective optimization of HEV fuel economy and emission using the selfadaptive differential evolution algorithm. IEEE T Veh Technol 2011; 60: Chen S. Optimal energy management of a hybrid electric powertrain system using improved particle swarm optimization. J Appl Energ 2015; 160: Zhang X, Song J, Tian Y, et al. Control strategy optimization of hybrid electric vehicle based on multi-objective genetic algorithm. J Mech Eng 2009; 45: Liu X. Study on dynamic parameter matching and simulation of pure electric vehicle. PhD Dissertation, Chongqing Jiaotong University, Chongqing, China, Zhou Q, Zhang W, Cash S, et al. Intelligent sizing of a series hybrid electric power-train system based on Chaosenhanced accelerated particle swarm optimization. Appl Energ 2017; 189: Jiang YC, Lin C and Cao WK. Parameters optimization for extended-range electric vehicle based on improved chaotic particle swarm optimization. Grid Distrib Comput 2016; 9: Dong N, Fang X and Wu A. A novel chaotic particle swarm optimization algorithm for parking space guidance. Math Probl Eng 2016; 2016: Pehlivanoğlu Y. A new particle swarm optimization method for the path planning of UAV in 3D environment. Aeronaut Space Technol 2012; 5: 1 14.