Battery Pack Grouping and Capacity Improvement for Electric Vehicles Based on a Genetic Algorithm

This paper proposes an optimal grouping method for battery packs of electric vehicles (EVs). Based on modeling the vehicle powertrain, analyzing the battery degradation performance and setting up the driving cycle of an EV, a genetic algorithm (GA) is applied to optimize the battery grouping topology with the objective of minimizing the total cost of ownership (TCO). The battery capacity and the serial and parallel amounts of the pack can thus be determined considering the influence of battery degradation. The results show that the optimized pack grouping can be solved by GA within around 9 min. Compared with the results of maximum discharge efficiency within a fixed lifetime, the proposed method can not only achieve a higher discharge efficiency, but also reduce the TCO by 2.29%. To enlarge the applications of the proposed method, the sensitivity to driving conditions is also analyzed to further prove the feasibility of the proposed method.


Introduction
Nowadays, consumption of fossil fuel and increment of greenhouse gas (GHG) emissions have become major challenges to environmental pollution.Vehicle electrification can mitigate these interactions by employing highly efficient motors to power vehicles [1,2].Currently, electric vehicles (EVs), hybrid electric vehicles (HEVs), plug-in HEVs (PHEVs) and fuel cell vehicles represent the developing directions of vehicle electrification.Among all the dominant solutions, EVs attract wide attention due to their relatively simple vehicle configuration, high operating efficiency, qualified driving performance and relatively mature energy storage technologies [3,4].In an EV, the energy storage system, usually a lithium-ion battery pack, becomes the essential and predominant component employed to drive the motor according to the driving demand.Therefore, it is considerably important and imperative to optimize the battery pack design with care to improve the vehicle economy.
For the battery pack in an EV, there exist hundreds of cells connected in series, parallel or even more complex topologies.These battery grouping technologies, including the series and parallel connected cell topologies, can greatly influence the whole pack specifications as well as the vehicle performance.Simply speaking, the series amount of all the cells can determine the vehicle rated voltage, which is related to the vehicle motor driving demands, and consequently, the parallel amount of the cells can influence the battery's overall capacity and the internal resistance.More importantly, the cost of the lithium-ion batteries occupies about 30-50% of the total cost of an EV.For instance, in the Nissan Leaf, the cost of the battery pack accounts for nearly 38% of the total vehicle cost [5].Given this, Energies 2017, 10, 439 2 of 15 the battery capacity and its grouping topology, i.e., the amount in series and in parallel, should be considered with care to provide an optimal grouping solution to meet the driving power demand, while simultaneously ensuring the vehicle economy [6].This consideration is the main research target in this paper.
Abundant research has been conducted on vehicle powertrain parameter optimization [7][8][9].The main research target usually involves battery and motor matching and their coordinated control.In [5], the whole process to determine the battery pack capacity is investigated to obtain the minimum battery cost for manufacturers.In view of the driving range, the pack size is designed to minimize the whole operation cost; however, the battery degradation influence should be considered with care to match different driving ranges.In [10], the battery size for PHEVs and extended range EVs (EREVs) with respect to different user types and different annual mileages is optimized considering energy consumption, total cost and greenhouse gas emission.The authors also analyze the battery size determination induced by policy influences in Germany.In addition, the battery life model is set up to be 12 years without considering the influence of degradation.In [11], a suitable battery size for prolonging the battery lifetime and decreasing the operating cost of a vehicle is evaluated; however, an important task is neglected in that the connection mode of the battery pack is not fully taken into account.In [12], an accurate battery life model is chosen for estimating the battery lifetime over different depths of discharge (DOD).It should be noted that there exists a tradeoff between DOD and battery size, as well as battery life.Furthermore, to acquire the minimum total cost of ownership (TCO), a global optimization algorithm is indirectly applied to make a closed-loop controller to reach the minimum operating cost.In [13], a multi-objective optimization, including energy management and battery capacity, is conducted based on hybrid fuel cell vehicles.After comparing this optimization with the optimized consequences of dynamic programming under different battery capacities, a soft-run energy management strategy is proposed.Meanwhile, the battery lifetime is also optimized considering battery size and hydrogen consumption.In [14], a tank-to-wheel method is applied to analyze the efficiency of energy recuperation and output considering the efficiency of converting electrical energy into mechanical energy.Based on an energy control algorithm, the energy efficiency is analyzed in detail when the battery size changes.In [15], a novel cost-optimal method is adopted to evaluate the charging, power management and battery degradation mitigation of PHEVs.The cost of electricity charged from the grid, fuel cost and battery life are considered and minimized based on convex programming (CP).Nonetheless, the pack grouping has not been analyzed yet due to its different serial-parallel construction focuses.
To summarize, it is crucial to consider the battery grouping topology and battery degradation to optimize the TCO.This criticality motivates the authors to conduct research into it.Although the TCO interactions between battery size and lifetime have been studied substantially, to the authors' best knowledge, researchers have seldom carried out studies taking the number of parallel modules and serial cells in a simplified connection pattern into account.Based on the above discussion, the overarching goal of this paper is herein proposed: to extend previous study by introducing a methodology for optimizing the lithium-ion battery pack grouping topology of EVs to satisfy customers' expectations with consideration for the TCO and battery degradation influence.
In order to optimize the TCO, it is imperative to analyze the vehicle powertrain, the battery degradation model, the daily operation cost, etc.The authors have considered all the related factors [16,17], as shown in Figure 1.The battery pack configuration can influence the cycle life, the discharging efficiency and the vehicle mass, and more directly, the battery cost, which is proportional to the battery capacity.Moreover, it can also determine the regeneration braking efficiency, and indirectly contribute to the operating cost.Therefore, these issues should be carefully taken into account to optimize the battery capacity and its grouping parameters.The proposed battery grouping topology is shown in Figure 2, in which the cell is a cylindrical lithium-ion battery.It can be observed that the number of parallel modules and the number of series modules should be determined.A commonly accepted design is that when the serial number increases, the voltage and internal resistance will consequently rise.Compared with the same power quantity, the battery current will decrease, thus improving the overall efficiency.In contrast, when the battery pack is short of voltage, the battery should release a larger current to satisfy the power demand, hence influencing the efficiency and cycle life [18].On the other hand, the amount of cells connected in parallel can affect the current and the whole pack resistance.More cells connected in parallel will certainly produce higher efficiency, while contributing a heavier weight and more expensive cost [6,19,20].Thus, a global optimization should be applied to determine these two topologies to minimize the battery lifetime cost and the total amount of all-life cost considering overall foregoing influences.To achieve the abovementioned target, the related vehicle powertrain model and the battery degradation model should be built first.Then the optimal number of parallel modules and serial cells can be modeled for minimizing the TCO under objective driving cycles.Since the target problem is non-convex and highly nonlinear with series discrete variables and constraints, genetic algorithm (GA) is herein elaborated to overcome the defect and to find the optimal solution [21][22][23][24], as seen in Section 3. Due to the fact that the comparison results are sensitive to the driving cycles, the corresponding analysis is performed.Finally, the results are validated to prove the feasibility of the proposed algorithm, followed by the conclusion marked in the end.The proposed battery grouping topology is shown in Figure 2, in which the cell is a cylindrical lithium-ion battery.It can be observed that the number of parallel modules and the number of series modules should be determined.A commonly accepted design is that when the serial number increases, the voltage and internal resistance will consequently rise.Compared with the same power quantity, the battery current will decrease, thus improving the overall efficiency.In contrast, when the battery pack is short of voltage, the battery should release a larger current to satisfy the power demand, hence influencing the efficiency and cycle life [18].On the other hand, the amount of cells connected in parallel can affect the current and the whole pack resistance.More cells connected in parallel will certainly produce higher efficiency, while contributing a heavier weight and more expensive cost [6,19,20].Thus, a global optimization should be applied to determine these two topologies to minimize the battery lifetime cost and the total amount of all-life cost considering overall foregoing influences.The proposed battery grouping topology is shown in Figure 2, in which the cell is a cylindrical lithium-ion battery.It can be observed that the number of parallel modules and the number of series modules should be determined.A commonly accepted design is that when the serial number increases, the voltage and internal resistance will consequently rise.Compared with the same power quantity, the battery current will decrease, thus improving the overall efficiency.In contrast, when the battery pack is short of voltage, the battery should release a larger current to satisfy the power demand, hence influencing the efficiency and cycle life [18].On the other hand, the amount of cells connected in parallel can affect the current and the whole pack resistance.More cells connected in parallel will certainly produce higher efficiency, while contributing a heavier weight and more expensive cost [6,19,20].Thus, a global optimization should be applied to determine these two topologies to minimize the battery lifetime cost and the total amount of all-life cost considering overall foregoing influences.To achieve the abovementioned target, the related vehicle powertrain model and the battery degradation model should be built first.Then the optimal number of parallel modules and serial cells can be modeled for minimizing the TCO under objective driving cycles.Since the target problem is non-convex and highly nonlinear with series discrete variables and constraints, genetic algorithm (GA) is herein elaborated to overcome the defect and to find the optimal solution [21][22][23][24], as seen in Section 3. Due to the fact that the comparison results are sensitive to the driving cycles, the corresponding analysis is performed.Finally, the results are validated to prove the feasibility of the proposed algorithm, followed by the conclusion marked in the end.To achieve the abovementioned target, the related vehicle powertrain model and the battery degradation model should be built first.Then the optimal number of parallel modules and serial cells can be modeled for minimizing the TCO under objective driving cycles.Since the target problem is non-convex and highly nonlinear with series discrete variables and constraints, genetic algorithm (GA) is herein elaborated to overcome the defect and to find the optimal solution [21][22][23][24], as seen in Section 3. Due to the fact that the comparison results are sensitive to the driving cycles, the corresponding analysis is performed.Finally, the results are validated to prove the feasibility of the proposed algorithm, followed by the conclusion marked in the end.

Vehicle Modeling and Analysis
Based on the aforementioned statement [5,16], the target function for minimizing the TCO should be related to the battery grouping topology and the driving conditions, as well as the circumstance conditions, such as the electricity price, battery lifetime, etc.Therefore, the target function can be formulated: where N parallel equals the number of parallel modules and N series states the number of series cells in a pack module.Y driving_cycle , depicting the driving behavior under different driving cycles, is modeled in relation to the daily trip length D. Z represents the underlying scenario assumptions, including the electricity price, the assumed battery lifetime, the production cost of batteries, motor cost, etc.
In order to determine the vehicle battery pack parameters, the vehicle model should be built first to provide the fundamental analysis of the demanded energy and power.

Vehicle Model
As shown in Figure 3, the powertrain structure of the proposed research target mainly comprises of a battery pack, a traction motor, a fixed gear, as well as a differential case [25,26].In this paper, a permanent-magnet synchronous electric motor (PMSM) is harnessed to power the vehicle.The main parameters of the motor are shown in Table 1.The power density is supposed to be 2.5 kW/kg [27].It can be found that the motor's maximum power is 80 kW, and the mass is 32 kg.In this paper, a high power lithium-ion cylindrical battery cell AHR32113 (A123 Systems, Livonia, USA) [28,29] is selected to construct a battery pack, whose detailed parameters are listed in Table 2, from which we can see that the nominal capacity is 4.4 Ah, the rated cell voltage is 3.3 V, and the operating temperature ranges from −30

Vehicle Modeling and Analysis
Based on the aforementioned statement [5,16], the target function for minimizing the TCO should be related to the battery grouping topology and the driving conditions, as well as the circumstance conditions, such as the electricity price, battery lifetime, etc.Therefore, the target function can be formulated: In order to determine the vehicle battery pack parameters, the vehicle model should be built first to provide the fundamental analysis of the demanded energy and power.

Vehicle Model
As shown in Figure 3, the powertrain structure of the proposed research target mainly comprises of a battery pack, a traction motor, a fixed gear, as well as a differential case [25,26].In this paper, a permanent-magnet synchronous electric motor (PMSM) is harnessed to power the vehicle.The main parameters of the motor are shown in Table 1.The power density is supposed to be 2.5 kW/kg [27].It can be found that the motor's maximum power is 80 kW, and the mass is 32 kg.In this paper, a high power lithium-ion cylindrical battery cell AHR32113 (A123 Systems, Livonia, USA) [28,29] is selected to construct a battery pack, whose detailed parameters are listed in Table 2, from which we can see that the nominal capacity is 4.4 Ah, the rated cell voltage is 3.3 V, and the operating temperature ranges from −30 °C to 55 °C.With the battery cell specification, the weight of battery pack m pack can be formulated as: where m cell is the weight of the battery cells, N bat is the number of battery cells, and ϕ is the packaging coefficient, which is specified as 1.25 in this paper [30].Now, the total vehicle mass m total can be calculated: where m veh is the vehicle chassis mass, m mot is the electric motor mass, and m pc is the mass of the passengers and cargo and is specified to be 150 kg on average.Based on the kinematic functions, the EV's driving dynamics model can be formulated: where P(t) means the driving power, P w (t) equals the power of the aerodynamic drag force, and P f (t) represents the power of rolling resistance.P i (t) represents the gradient resistance power.P j (t) is the acceleration (or deceleration) power demand.P acc , equaling 269 W, is the accessory power [31].g represents the gravitational acceleration, and equals 9.81 m/s 2 .α namely the slope of the EV, is near zero most of the time.J total denotes the total energy consumption in each driving range in kWh.
The remaining vehicle parameters are shown in Table 3.For the EV, the regenerative braking system can automatically recover the braking energy.To optimize the battery capacity, the braking performance should be properly considered.To simplify the problem, the regenerative braking function can be indicated as: where W represents the regenerative braking energy in kWh, and a k and v k are the acceleration/deceleration in m/s 2 and velocity in m/s at moment k, respectively.µ denotes the factor of average energy recovery and equals 0.323, as discussed in [32].t initial and t terminal state the initial time and terminal time of the braking operation.

Battery Model
In order to study the interactions between the amount of series and parallel cells, the discharge efficiency influences induced by the resistance of the pack should be analyzed in advance.Given the amount of series and parallel units [33], the connected mode among each unit can be simplified and indexed, as shown in Figure 4.Each cell within a parallel unit connects with R i , which is equal to the resistance of the copper conductors, and the contact resistor R C expresses the aluminum clamp resistance.Their values are specified to be 0.1 milliohm and 3 milliohm, respectively [34], and the cell resistor is set to be 10 milliohm.In this paper, the inconsistency of the pack is not taken into account for simplification.Therefore, the resistance function R parallel_unit of the parallel unit can be furnished via the equivalent transformation between Y and ∆ connections, which is only determined by the parallel amount of cells after ensuring the aforementioned pre-set values: Hence, the overall resistance R pack of the pack can be given as: The value of current I and discharge efficiency η discharge can be described as: where U is the open-circuit voltage of the battery pack and equals 3.3 × N series in this paper.
= ×  1000 3600 where W represents the regenerative braking energy in kWh, and k a and k v are the acceleration/deceleration in 2 m/s and velocity in m/s at moment k , respectively.μ denotes the factor of average energy recovery and equals 0.323, as discussed in [32].initial t and terminal t state the initial time and terminal time of the braking operation.

Battery Model
In order to study the interactions between the amount of series and parallel cells, the discharge efficiency influences induced by the resistance of the pack should be analyzed in advance.Given the amount of series and parallel units [33], the connected mode among each unit can be simplified and indexed, as shown in Figure 4.Each cell within a parallel unit connects with i R , which is equal to the resistance of the copper conductors, and the contact resistor C R expresses the aluminum clamp resistance.Their values are specified to be 0.1 milliohm and 3 milliohm, respectively [34], and the cell resistor is set to be 10 milliohm.In this paper, the inconsistency of the pack is not taken into account for simplification.Therefore, the resistance function of the parallel unit can be furnished via the equivalent transformation between Y and ∆ connections, which is only determined by the parallel amount of cells after ensuring the aforementioned pre-set values: Hence, the overall resistance pack R of the pack can be given as: The value of current I and discharge efficiency discharge η can be described as: where U is the open-circuit voltage of the battery pack and equals 3.3 × series N in this paper.Traditionally, the widely adopted battery degradation model is calculated as an empirical function of the charge/discharge rate, DOD and battery capacity [35][36][37].Indeed, the battery degradation model is rather complex and it is quite difficult to build a uniform relationship between it and the usage operation.This is not the research focus of this paper, and a simplified method, i.e., the Peterson model, Energies 2017, 10, 439 7 of 15 is herein employed [38].Considering one charge process per day, the battery energy consumption can be described as: , discharge where w charge and w discharge are the energy processed for the daily charging and discharging respectively.
η discharge represents the discharge efficiency.L day denotes the driving range per day.E dis means the electric energy consumption per kilometer in Wh/km under the identified driving cycle, and is proportional to the vehicle mass.Consequently, by inputting ( 6) through ( 8) into ( 9), the battery lifetime T bat can be deduced: where E cell is a single battery cell's energy capacity in kWh/cell.α discharge and α charge are the battery discharging and charging degradation factors and equal 3.46 × 10 −5 and 1.72 × 10 −5 , respectively [39].Note that D drive denotes the driving days per year and equals 300 days in this paper, and r EOL indexes the window function of DOD.In this paper, we suppose that when the battery capacity degrades to 80% rated values, the battery should be replaced [40].

Vehicle Cost Model
Based on the above discussion, the TCO consists of four parts: the vehicle body cost, powertrain cost, pack cost, as well as the operation cost.In this paper, an A00 class vehicle chassis is selected as the target platform, of which the mass is set to be 800 kg and the average cost is specified to be 2686 euros.The average price of a lithium-ion battery is specified to be 268.6 euros/kWh [27].Now, the total cost of the battery cells C pack can be formulated: The average annual ownership cost C year for the EV can be calculated: where T veh is the life cycle of both the vehicle and battery pack in years, C basic_vehicle is the vehicle basic cost, C p is the powertrain cost and P e , equaling 0.07 euros/kWh, represents the average electricity price based on investigation [41].
Based on the described model, GA is then applied to optimize the battery grouping configuration and to find the optimal solutions.

Optimization Framework
In this paper, GA is employed to optimize the battery pack grouping topology for EVs.GA is qualified to control the highly nonlinear system, and is strongly skilled in finding the optimal solutions.Typically, the process of GA includes encoding, selection, crossover and mutation operations [23].During each generation, the individuals with a higher fitness value are elected to produce new individuals for the next offspring.The optimal populations are unchangeable and are directly considered as elitists for the next offspring.By selecting two individuals as parents to be hybridized, new offspring, namely children individuals, can be generated.The mutation process is integrated to reduce the possibility of converging into a local optimum by choosing a random element of a chromosome.After these steps, the whole process continuously iterates until meeting the terminal Energies 2017, 10, 439 8 of 15 conditions [42,43].Generally, the termination conditions include: (1) finding a satisfactory solution, (2) reaching the set fitness values, and (3) exceeding the budget time or number of generations.
By optimization, the number of generations, population, crossover rate and mutation rate are set to be 60, 40, 0.8 and 0.2, respectively.The whole application is depicted in Figure 5. First, the initialization is set up, followed by calculating the TCO.Then, GA is applied to minimize the TCO until the termination conditions are achieved.Finally, the whole calculation process ends and the battery grouping information and TCO are determined.
Typically, the process of GA includes encoding, selection, crossover and mutation operations [23].During each generation, the individuals with a higher fitness value are elected to produce new individuals for the next offspring.The optimal populations are unchangeable and are directly considered as elitists for the next offspring.By selecting two individuals as parents to be hybridized, new offspring, namely children individuals, can be generated.The mutation process is integrated to reduce the possibility of converging into a local optimum by choosing a random element of a chromosome.After these steps, the whole process continuously iterates until meeting the terminal conditions [42,43].Generally, the termination conditions include: (1) finding a satisfactory solution, (2) reaching the set fitness values, and (3) exceeding the budget time or number of generations.
By optimization, the number of generations, population, crossover rate and mutation rate are set to be 60, 40, 0.8 and 0.2, respectively.The whole application is depicted in Figure 5. First, the initialization is set up, followed by calculating the TCO.Then, GA is applied to minimize the TCO until the termination conditions are achieved.Finally, the whole calculation process ends and the battery grouping information and TCO are determined.From the foregoing elaboration, the optimization can be processed with a series of optimal constraints.Due to the regulation that the most basic quality standard of battery life is 8 years [44], the limitation with respect to battery life is enforced to be more than 8 years considering its routine employment by a user.It is necessary to mention that due to the powering demand and voltage ratings of the motor and its driving system, the amount of series cells series N is set up to range from 80 to 120 to ensure the driving efficiency.In addition, the highest value of parallel N is limited to 100, since when the parallel amount is 100 and the series amount is 80, the total battery available capacity is approximately equal to 117 kWh, which should be enough to power the vehicle.Moreover, in view of the interaction between DOD and all-life cycles, the limited DOD ranges from 0 to 0.8.The other constraints relative to the vehicle components, e.g., maximum motor power, are regarded as common settings.To sum up, all the constraints can be formulated as: In order to optimize the battery size, the driving cycle should be determined with attention.In this paper, the driving pattern of users derived from Germany is adopted to simulate the daily driving distribution for automobiles, as shown in Figure 6  From the foregoing elaboration, the optimization can be processed with a series of optimal constraints.Due to the regulation that the most basic quality standard of battery life is 8 years [44], the limitation with respect to battery life is enforced to be more than 8 years considering its routine employment by a user.It is necessary to mention that due to the powering demand and voltage ratings of the motor and its driving system, the amount of series cells N series is set up to range from 80 to 120 to ensure the driving efficiency.In addition, the highest value of N parallel is limited to 100, since when the parallel amount is 100 and the series amount is 80, the total battery available capacity is approximately equal to 117 kWh, which should be enough to power the vehicle.Moreover, in view of the interaction between DOD and all-life cycles, the limited DOD ranges from 0 to 0.8.The other constraints relative to the vehicle components, e.g., maximum motor power, are regarded as common settings.To sum up, all the constraints can be formulated as: In order to optimize the battery size, the driving cycle should be determined with attention.In this paper, the driving pattern of users derived from Germany is adopted to simulate the daily driving distribution for automobiles, as shown in Figure 6 [45].Here, a probability density function H(s) can be fitted to calculate the average daily driving length, which will not interact with the investigative accuracy.The daily driving length S daily can be determined via integral operation: Energies 2017, 10, 439 9 of 15 where s denotes the daily driving distance.Based on Figure 6, the daily and annual driving lengths can be calculated, which are 41.22 km and 12366 km, respectively.
Energies 2017, 10, 439 9 of 15 ( ) H s can be fitted to calculate the average daily driving length, which will not interact with the investigative accuracy.The daily driving length daily S can be determined via integral operation: where s denotes the daily driving distance.Based on Figure 6, the daily and annual driving lengths can be calculated, which are 41.22 km and 12366 km, respectively.In this paper, a mixed driving cycle is employed to construct the daily driving conditions, which is composed of 40% urban cycle and 60% extra-urban cycle [10].In [10], the authors calculated the weights of 40% and 60% through the average speed and average trip distance when driving 365 days per year; however, the situation of driving every day in a year rarely happens in the real world.Therefore, we propose that the share is still set to be 40% urban cycle and 60% extra-urban cycle but with 300 driving days per year.Here, based on the New European Driving Cycle (NEDC), exactly combined with the above two general kinds of portions, a mixed driving cycle is selected to simulate a typical cycle to estimate a real-world cycle, as indicated in Figure 7 [46].Its average speed is 32.21 km/h, slightly lower than the 33.60 km/h under the NEDC cycle.Commonly, the basic driving performance requirement is the most pivotal portion of the whole cycle, and even the most attractive to the customers or manufacturers to some extent.In this paper, the basic EV performance requirements derived from China Battery Electric Passenger Cars Specifications (GB/T 28382-2012) [47] are implemented, as listed in Table 4.The acceleration time from 0 to 50 km/h is limited to 10 s, and the maximum grade ability is specified to be more than 20%.In addition, to satisfy the daily driving demand, the all-electric range (AER) should be more than 80 km.In this paper, a mixed driving cycle is employed to construct the daily driving conditions, which is composed of 40% urban cycle and 60% extra-urban cycle [10].In [10], the authors calculated the weights of 40% and 60% through the average speed and average trip distance when driving 365 days per year; however, the situation of driving every day in a year rarely happens in the real world.Therefore, we propose that the share is still set to be 40% urban cycle and 60% extra-urban cycle but with 300 driving days per year.Here, based on the New European Driving Cycle (NEDC), exactly combined with the above two general kinds of portions, a mixed driving cycle is selected to simulate a typical cycle to estimate a real-world cycle, as indicated in Figure 7 [46].Its average speed is 32.21 km/h, slightly lower than the 33.60 km/h under the NEDC cycle.
Energies 2017, 10, 439 9 of 15 ( ) H s can be fitted to calculate the average daily driving length, which will not interact with the investigative accuracy.The daily driving length daily S can be determined via integral operation: =  ( ) daily s S Hssds (14) where s denotes the daily driving distance.Based on Figure 6, the daily and annual driving lengths can be calculated, which are 41.22 km and 12366 km, respectively.In this paper, a mixed driving cycle is employed to construct the daily driving conditions, which is composed of 40% urban cycle and 60% extra-urban cycle [10].In [10], the authors calculated the weights of 40% and 60% through the average speed and average trip distance when driving 365 days per year; however, the situation of driving every day in a year rarely happens in the real world.Therefore, we propose that the share is still set to be 40% urban cycle and 60% extra-urban cycle but with 300 driving days per year.Here, based on the New European Driving Cycle (NEDC), exactly combined with the above two general kinds of portions, a mixed driving cycle is selected to simulate a typical cycle to estimate a real-world cycle, as indicated in Figure 7 [46].Its average speed is 32.21 km/h, slightly lower than the 33.60 km/h under the NEDC cycle.Commonly, the basic driving performance requirement is the most pivotal portion of the whole cycle, and even the most attractive to the customers or manufacturers to some extent.In this paper, the basic EV performance requirements derived from China Battery Electric Passenger Cars Specifications (GB/T 28382-2012) [47] are implemented, as listed in Table 4.The acceleration time from 0 to 50 km/h is limited to 10 s, and the maximum grade ability is specified to be more than 20%.In addition, to satisfy the daily driving demand, the all-electric range (AER) should be more than 80 km.Commonly, the basic driving performance requirement is the most pivotal portion of the whole cycle, and even the most attractive to the customers or manufacturers to some extent.In this paper, the basic EV performance requirements derived from China Battery Electric Passenger Cars Specifications (GB/T 28382-2012) [47] are implemented, as listed in Table 4.The acceleration time from 0 to 50 km/h is limited to 10 s, and the maximum grade ability is specified to be more than 20%.In addition, to satisfy the daily driving demand, the all-electric range (AER) should be more than 80 km.
The next step is to perform the corresponding analysis and results validation to prove the feasibility of the proposed method.

Results and Discussion
Due to the fact that the battery pack is a complicated system, an assumption is herein made that the temperature of the battery pack is consistent at 25 • C.Moreover, in order to ensure the accuracy of the modeling, the powertrain of the EV is taken from Autonomie, which is an effective and well-known simulation platform for EVs and PHEVs [48].The whole validation is shown in Figure 8, yielding that the setting parameters are capable of obtaining a good result within 60 generations.The total calculation duration is within 9 min with a desktop computer using a Core i7 CPU and 8 gigabyte memory storage.The average TCO per year and the related results are exhibited in Table 5.After optimization, the parallel number of cells of the battery pack is 40, and the serial number is 113.The battery cost is 1687.93euros, which occupies an overriding portion of 66.13% in the TCO and is apparently inverse to the sum of the operation cost and based vehicle cost (33.87%).In addition, the acceleration times to go from 0 to km/h and from 0 to 80 km/h are 6.7 s and 10.6 s, respectively.The grade ability performance and maximum speed requirement can also be satisfied.In addition, the driving range is 298.96km, reaching the design requirement.The next step is to perform the corresponding analysis and results validation to prove the feasibility of the proposed method.

Results and Discussion
Due to the fact that the battery pack is a complicated system, an assumption is herein made that the temperature of the battery pack is consistent at 25 °C.Moreover, in order to ensure the accuracy of the modeling, the powertrain of the EV is taken from Autonomie, which is an effective and wellknown simulation platform for EVs and PHEVs [48].The whole validation is shown in Figure 8, yielding that the setting parameters are capable of obtaining a good result within 60 generations.The total calculation duration is within 9 min with a desktop computer using a Core i7 CPU and 8 gigabyte memory storage.The average TCO per year and the related results are exhibited in Table 5.After optimization, the parallel number of cells of the battery pack is 40, and the serial number is 113.The battery cost is 1687.93euros, which occupies an overriding portion of 66.13% in the TCO and is apparently inverse to the sum of the operation cost and based vehicle cost (33.87%).In addition, the acceleration times to go from 0 to 50 km/h and from 0 to 80 km/h are 6.7 s and 10.6 s, respectively.The grade ability performance and maximum speed requirement can also be satisfied.In addition, the driving range is 298.96km, reaching the design requirement.Figure 9 shows the TCO variations when the parallel and serial amounts of cells are supposed separately to be constant.It can be observed that when the parallel amount is set to be 40, the TCO reaches the minimum point when the serial amount is 113; similarly, when the serial amount is set to be 113, the TCO can achieve the minimum value of a parallel amount of 40.Thus, this observation proves that the grouping amounts of 113 and 40 reach the optimal solution and therefore justifies the feasibility of the proposed algorithm.Essentially, the main target of the proposed optimization method is to achieve the highest efficiency of the battery operation.
Energies 2017, 10, 439 11 of 15 Figure 9 shows the TCO variations when the parallel and serial amounts of cells are supposed separately to be constant.It can be observed that when the parallel amount is set to be 40, the TCO reaches the minimum point when the serial amount is 113; similarly, when the serial amount is set to be 113, the TCO can achieve the minimum value of a parallel amount of 40.Thus, this observation proves that the grouping amounts of 113 and 40 reach the optimal solution and therefore justifies the feasibility of the proposed algorithm.Essentially, the main target of the proposed optimization method is to achieve the highest efficiency of the battery operation.Based on the discussion in Section 2 above, the efficiency and lifetime of the battery pack are two key variables directly affecting the remaining variables.To further extend the analysis, these two variables are herein discussed.To easily compare these influences, the battery grouping topology is determined with regard to optimizing the highest efficiency with a fixed lifetime.The corresponding result is shown in Figure 10.It can be observed that the proposed solution can reduce the TCO by a maximum value of 2.29%, compared with those based on the constraint of fixed time.Now it can be proven that the determined topology can reach the highest efficiency by way of the proposed algorithm.To analyze the influences between the solutions and AER, an AER with increasing operation lifetime is taken into account, of which the groups are also selected by the aforementioned method, i.e., highest efficiency and constant pre-set lifetime, as shown in Figure 11.Apparently, the AER drops linearly, thus the function can be fitted with a linear function.In addition, it is necessary to note that the slopes of attenuation with rising operation lifetime are almost consistent and the differences between two lines is about 35 km.Based on the discussion in Section 2 above, the efficiency and lifetime of the battery pack are two key variables directly affecting the remaining variables.To further extend the analysis, these two variables are herein discussed.To easily compare these influences, the battery grouping topology is determined with regard to optimizing the highest efficiency with a fixed lifetime.The corresponding result is shown in Figure 10.It can be observed that the proposed solution can reduce the TCO by a maximum value of 2.29%, compared with those based on the constraint of fixed time.Now it can be proven that the determined topology can reach the highest efficiency by way of the proposed algorithm.To analyze the influences between the solutions and AER, an AER with increasing operation lifetime is taken into account, of which the groups are also selected by the aforementioned method, i.e., highest efficiency and constant pre-set lifetime, as shown in Figure 11.Apparently, the AER drops linearly, thus the function can be fitted with a linear function.In addition, it is necessary to note that the slopes of attenuation with rising operation lifetime are almost consistent and the differences between two lines is about 35 km.
Energies 2017, 10, 439 11 of 15 Figure 9 shows the TCO variations when the parallel and serial amounts of cells are supposed separately to be constant.It can be observed that when the parallel amount is set to be 40, the TCO reaches the minimum point when the serial amount is 113; similarly, when the serial amount is set to be 113, the TCO can achieve the minimum value of a parallel amount of 40.Thus, this observation proves that the grouping amounts of 113 and 40 reach the optimal solution and therefore justifies the feasibility of the proposed algorithm.Essentially, the main target of the proposed optimization method is to achieve the highest efficiency of the battery operation.Based on the discussion in Section 2 above, the efficiency and lifetime of the battery pack are two key variables directly affecting the remaining variables.To further extend the analysis, these two variables are herein discussed.To easily compare these influences, the battery grouping topology is determined with regard to optimizing the highest efficiency with a fixed lifetime.The corresponding result is shown in Figure 10.It can be observed that the proposed solution can reduce the TCO by a maximum value of 2.29%, compared with those based on the constraint of fixed time.Now it can be proven that the determined topology can reach the highest efficiency by way of the proposed algorithm.To analyze the influences between the solutions and AER, an AER with increasing operation lifetime is taken into account, of which the groups are also selected by the aforementioned method, i.e., highest efficiency and constant pre-set lifetime, as shown in Figure 11.Apparently, the AER drops linearly, thus the function can be fitted with a linear function.In addition, it is necessary to note that the slopes of attenuation with rising operation lifetime are almost consistent and the differences between two lines is about 35 km.

Conclusions
This paper develops an optimization methodology for battery grouping configuration with the objective to minimize the TCO based on the driving patterns of users.Based on modeling the vehicle powertrain and the battery degradation, GA is applied to acquire the expected number of parallel modules and serial cells with the target of optimizing the TCO.The proposed approach is proven to be effective in optimizing the battery economy.In addition, aiming to implement this approach in practice, different cycles are employed to extend the applications so that manufacturers can adaptively regulate the pack grouping topologies according to the actual demand.
Currently, the driving cycles used in this paper are based on knowing the driving patterns.To enlarge the applications of the proposed method, an actual vehicle test and hardware-in-the-loop test need to be carried out as the next step.Additionally, further investigation of traveling habits and related factors, such as influences of the weather/climate and traffic information, should also be included in future studies.

Figure 1 .
Figure 1.Structure of the battery combination optimization problem for BEVs.

Figure 2 .
Figure 2. The battery pack in the parallel-series manner.

Figure 1 .
Figure 1.Structure of the battery combination optimization problem for BEVs.

Figure 2 .
Figure 2. The battery pack in the parallel-series manner.

Figure 2 .
Figure 2. The battery pack in the parallel-series manner.
driving behavior under different driving cycles, is modeled in relation to the daily trip length D .Z represents the underlying scenario assumptions, including the electricity price, the assumed battery lifetime, the production cost of batteries, motor cost, etc.

Figure 3 .
Figure 3.The EV powertrain topology of a central drive powertrain.

Figure 3 .
Figure 3.The EV powertrain topology of a central drive powertrain.

Figure 4 .
Figure 4. Schematic of series connection among each parallel unit.

Figure 4 .
Figure 4. Schematic of series connection among each parallel unit.

Figure 5 .
Figure 5.The optimization procedure using GA.

Figure 5 .
Figure 5.The optimization procedure using GA.

Figure 6 .
Figure 6.Distribution of daily driving distance.

Figure 7 .
Figure 7.The mixture cycle based on NEDC.

Figure 6 .
Figure 6.Distribution of daily driving distance.

Figure 6 .
Figure 6.Distribution of daily driving distance.

Figure 7 .
Figure 7.The mixture cycle based on NEDC.

Figure 7 .
Figure 7.The mixture cycle based on NEDC.

Figure 8 .
Figure 8.The GA method calculation process.

Figure 8 .
Figure 8.The GA method calculation process.

Figure 10 .
Figure 10.The comparison of TCO with the highest discharge efficiency and varying fixed lifetime.(a) The comparison of TCO; (b) The varied discharge efficiency.

Figure 10 .
Figure 10.The comparison of TCO with the highest discharge efficiency and varying fixed lifetime.(a) The comparison of TCO; (b) The varied discharge efficiency.Figure 10.The comparison of TCO with the highest discharge efficiency and varying fixed lifetime.(a) The comparison of TCO; (b) The varied discharge efficiency.

Figure 10 .
Figure 10.The comparison of TCO with the highest discharge efficiency and varying fixed lifetime.(a) The comparison of TCO; (b) The varied discharge efficiency.Figure 10.The comparison of TCO with the highest discharge efficiency and varying fixed lifetime.(a) The comparison of TCO; (b) The varied discharge efficiency.

Table 1 .
The specifications of the baseline electric motor.

Table 1 .
The specifications of the baseline electric motor.

Table 4 .
Vehicle performance requirements.

Table 4 .
Vehicle performance requirements.

Table 5 .
The optimal solutions by genetic algorithm.

Table 5 .
The optimal solutions by genetic algorithm.

Table 6 .
The results under a mixed cycle combined with UDDS and HWFET cycles.