Energy Management of Parallel-Connected Cells in Electric Vehicles Based on Fuzzy Logic Control

Inconsistencies that are associated with parallel-connected cells used in electric vehicles induce varied states of charge (SOCs) in each cell. Thus, loop current in the battery pack is inevitable, and this reduces overall capacity, energy utilization rate, and pack lifetime. However, no method is available to address loop current. To reduce loop current and the resulting battery inconsistency, a parallel-connected cell pack (PCCP) model that considers thermal effects is established, and a novel Simscape model that is based on PCCP is successfully constructed. Furthermore, the strategy of parallel-connected cell energy management (PCCEM) is proposed to utilize fuzzy logic control (FLC) strategy, which automatically adjusts the number of cells in a circuit in accordance with the load demand, and turns on the first N switches in the corresponding SOC order. The New European Driving Cycle (NEDC) driving cycle simulation shows that the PCCEM strategy considerably reduces loop current and improves the consistency of battery performance and the utilization rate of battery power.


Introduction
The lithium battery is associated with numerous concerns because of its high-energy density, large capacity, and no memory effect [1][2][3].Nonetheless, the lithium battery has been considered one of the major viable energy storage media in electric cars [4][5][6][7].To meet the voltage and power requirements of cars, if we use small-capacity and low-power lithium cells, as in Tesla's Model S, then the cells should be connected in parallel or/and series fashion [8][9][10][11].A parallel-connected cell pack (PCCP) can be regarded as a larger single cell by omitting the differences among cells.However, given the differences in internal resistance, capacity, aging property, and temperature of cells, the discrepancy in the state of charge (SOC) and/or voltage among the parallel batteries is inevitable [12][13][14].As a result, circulating current among batteries with different SOC and/or voltage is unavoidable.Connecting cells in parallel, which entails a number of engineering complications that must be resolved, is not the first choice.Still, when we apply cells into a large-capacity demand area, such as an electric vehicle, we have no choice but to connect cells in parallel.However, the overall performance of the Li-ion battery deteriorates over time because the system is repeatedly cycled through iterations of charge and discharge.Although the cells used in an electric vehicle at a state of health (SOH) = 100% would have almost the same capacities and SOC, the battery would have deteriorated capacities after being cells would be left with less than 80% of their original capacities, and some cells would have lower SOC than others.As shown in Figure 1, several cells have high capacity or high SOC, whereas others have relatively low capacity or low SOC.The gradually increasing inconsistencies in the cells will increase the difference in working current, which results in uneven temperature distribution.Figure 2a shows that loop current will occur between parallel cells in standby mode when all the switches are closed.Consequently, various problems, including increased battery temperature, energy loss and reduced service time, will occur simultaneously.Moreover, the uneven distribution of the current and temperature will produce different cell aging rates.Many researchers have discussed cell inconsistency from different perspectives.Bruen et al. [15] evaluated the variation in energy among cells and investigated its probable effect on the design of highly parallelized battery systems.Yang et al. [1] introduced the influence of uneven temperature distribution on PCCP discharge characteristics and service life.They pointed out that the cell with higher temperature will have a higher discharge current during the early stage of discharge and a lower current during the later stages.Moo et al. studied the current characteristics of PCCP charge and discharge.They reported that the battery pack releases more energy and increases energy efficiency through controlling the PCCP's discharge sequence.Hsu et al. [16] introduced an energy management method that increases uninterrupted power supply (UPS) energy output by alternately driving the load.Barre and Maher et al. [17,18] pointed out that the current of each cell cannot be measured by traditional battery management system (BMS) and that the differences among currents cannot be detected.Offer et al. [19] reported that cells closer to the main circuit have higher currents than those of cells that are further from the main circuit.This condition causes capacity difference in the long term, which indicates that parallel cells have different SOC values despite having the same terminal voltage.Gong et al. [4] studied the characteristics of the parallel four-cell pack and reported the adverse effects of inconsistency on the pack, such as different aging degrees and accelerated aging rate, which in turn increases the inconsistency of the pack.Bruen et al. [20] established a model and reported that significant differences in current occur and affect cell aging rate and temperature distribution.Energies 2017, 10, 404 2 of 13 cells would be left with less than 80% of their original capacities, and some cells would have lower SOC than others.As shown in Figure 1, several cells have high capacity or high SOC, whereas others have relatively low capacity or low SOC.The gradually increasing inconsistencies in the cells will increase the difference in working current, which results in uneven temperature distribution.Figure 2a shows that loop current will occur between parallel cells in standby mode when all the switches are closed.Consequently, various problems, including increased battery temperature, energy loss and reduced service time, will occur simultaneously.Moreover, the uneven distribution of the current and temperature will produce different cell aging rates.Many researchers have discussed cell inconsistency from different perspectives.Bruen et al. [15] evaluated the variation in energy among cells and investigated its probable effect on the design of highly parallelized battery systems.Yang et al. [1] introduced the influence of uneven temperature distribution on PCCP discharge characteristics and service life.They pointed out that the cell with higher temperature will have a higher discharge current during the early stage of discharge and a lower current during the later stages.Moo et al. studied the current characteristics of PCCP charge and discharge.They reported that the battery pack releases more energy and increases energy efficiency through controlling the PCCP's discharge sequence.Hsu et al. [16] introduced an energy management method that increases uninterrupted power supply (UPS) energy output by alternately driving the load.Barre and Maher et al. [17,18] pointed out that the current of each cell cannot be measured by traditional battery management system (BMS) and that the differences among currents cannot be detected.Offer et al. [19] reported that cells closer to the main circuit have higher currents than those of cells that are further from the main circuit.This condition causes capacity difference in the long term, which indicates that parallel cells have different SOC values despite having the same terminal voltage.Gong et al. [4] studied the characteristics of the parallel four-cell pack and reported the adverse effects of inconsistency on the pack, such as different aging degrees and accelerated aging rate, which in turn increases the inconsistency of the pack.Bruen et al. [20] established a model and reported that significant differences in current occur and affect cell aging rate and temperature distribution.Many researchers have discussed cell inconsistency from different perspectives.Bruen et al. [15] evaluated the variation in energy among cells and investigated its probable effect on the design of highly parallelized battery systems.Yang et al. [1] introduced the influence of uneven temperature distribution on PCCP discharge characteristics and service life.They pointed out that the cell with higher temperature will have a higher discharge current during the early stage of discharge and a lower current during the later stages.Moo et al. studied the current characteristics of PCCP charge and discharge.They reported that the battery pack releases more energy and increases energy efficiency through controlling the PCCP's discharge sequence.Hsu et al. [16] introduced an energy management method that increases uninterrupted power supply (UPS) energy output by alternately driving the load.Barre and Maher et al. [17,18] pointed out that the current of each cell cannot be measured by traditional battery management system (BMS) and that the differences among currents cannot be detected.Offer et al. [19] reported that cells closer to the main circuit have higher currents than those of cells that are further from the main circuit.This condition causes capacity difference in the long term, which indicates that parallel cells have different SOC values despite having the same terminal voltage.Gong et al. [4] studied the characteristics of the parallel four-cell pack and reported the adverse effects Energies 2017, 10, 404 3 of 13 of inconsistency on the pack, such as different aging degrees and accelerated aging rate, which in turn increases the inconsistency of the pack.Bruen et al. [20] established a model and reported that significant differences in current occur and affect cell aging rate and temperature distribution.
The above analysis indicates that several studies have provided valuable results.However, prior studies have mainly focused on phenomenon analysis, modeling, and simulation.No rational method has been proposed to handle circulating currents in parallel-connected cells.Therefore, the present study investigates the energy management of parallel-connected cells.Figure 2 shows the operating principle.To minimize the circulating current, we can control the switch S11-SMNon the basis of power demand, which ensures that the cell with a higher SOC has a higher discharge priority.If the discharge current is low, then only part of the battery will be involved in discharging.Cells with lower SOCs will be disconnected from the circuit.Then, the cell with a lower state of charge (SOH) is isolated from the pack when its performance decays to a certain value by turning off the corresponding switch.Therefore, this study attempts to reduce the circulating current by appropriately controlling each cell's discharge process.
First, this study builds a parallel-connected cell pack (PCCP) model that considers thermal effect.The PCCP model is constructed on the foundation of a previous work [21].Then, the PCCEM strategy is proposed.The strategy changes the working mode of each cell in accordance with real conditions and load on the basis of fuzzy logic control (FLC).Finally, the PCCEM algorithm is applied to the PCCP model to evaluate its influence on a highly parallel battery system.The NEDC driving cycle simulation and hardware-in-the-loop experiments show that this strategy improves the consistency of battery performance and the utilization rate of battery power.
The experimental method and model in this article serve as a guide for the energy management of parallel-connected cells to optimize the energy utilization rate.Section 2 provides the single cell model and Simulink model of parallel-connected cells.Section 3 presents the PCCEM strategy.Section 4 discusses the simulation results.Section 5 presents the conclusions and prospects.

First-Order RC Equivalent Circuit Model
PCCP energy management is based on accurate and real-time SOC estimation.Thus, the battery SOC must be accurate under all conditions [22].Hence, constructing a battery model to estimate the SOC is advisable.A number of researchers [17,18,23,24] mention several battery models and describe cell charge and discharge characteristics.The first-order equivalent circuit model is very common in simulating cell voltage change because the model is relatively simple, can easily obtain parameters, and runs in real time.The core goal of the equivalent circuit model is to simulate the actual battery voltage response of the current input [20].Huria et al. [21] reported that the first-order resistance capacitance (RC) circuit model is precise for most industrial applications.Increasing the order will significantly elevate the amount of calculation, but will not increase accuracy greatly [25].The first-order equivalent circuit model is advocated by Thomas and Huria et al. for its accuracy and feasibility, and many researchers [26][27][28][29] have adapted this model in their studies.The equivalent circuit model, which contains a voltage source, series resistor, and a pair of RC, accurately describes cell characteristics.In practice, Coulomb counting method suffers from drift.According to the suggestion of Seaman et al. [27], we made the resistors and capacitance dependent on SOC and temperature to increase the fidelity of the battery model.Figure 3 shows the classic first-order RC equivalent circuit model.
In Figure 3, E m represents the open circuit voltage; R 1 and C 1 represent the resistance and capacitance of polarization effect, respectively; and R 0 represents the ohmic resistance of battery.E m , R 1 , C 1 , and R 0 are functions of SOC and temperature.Specifically, these four elements can be obtained from the two-dimensional look-up tables: simulating cell voltage change because the model is relatively simple, can easily obtain parameters, and runs in real time.The core goal of the equivalent circuit model is to simulate the actual battery voltage response of the current input [20].Huria et al. [21] reported that the first-order resistance capacitance (RC) circuit model is precise for most industrial applications.Increasing the order will significantly elevate the amount of calculation, but will not increase accuracy greatly [25].The first-order equivalent circuit model is advocated by Thomas and Huria et al. for its accuracy and feasibility, and many researchers [26][27][28][29] have adapted this model in their studies.The equivalent circuit model, which contains a voltage source, series resistor, and a pair of RC, accurately describes cell characteristics.In practice, Coulomb counting method suffers from drift.According to the suggestion of Seaman et al. [27], we made the resistors and capacitance dependent on SOC and temperature to increase the fidelity of the battery model.Figure 3 shows the classic first-order RC equivalent circuit model.Those look-up tables were obtained by using a parameter estimation tool in Simulink Design Optimization with series pulse discharge data under different temperatures.Additionally, these look-up tables were selected based on seven different points of SOC.The pulse discharge curve for each temperature was run individually through estimation.This will produce a set of one-dimensional look-up tables versus SOC for the four parameters at each temperature.Simulink Design Optimization iteratively simulated the discharge profile in Simscape while comparing the simulation results with experimental data.Then, using nonlinear least squares algorithm, the error gradient across each of the 28 parameters (four tables × seven breakpoints) to minimize the sum of squared error.
The battery model in the reference [21], which is a fidelity electrical model with thermal dependence, was used for characterization and simulation of high-power lithium battery cells.The cell temperature is computed by using the thermal model in [21].The inner cell temperature is assumed to be uniform and taken as the average temperature inside the cell.SOC was calculated through Coulomb counting by using Formula (5).SOC (t 0 ) is calculated based on each measured open circuit voltage (OCV) reading using linear interpolation.SOC (t 0 ) = f (OCV, T), and only when the dv/dt < 4 uV/s condition is satisfied are OCV readings taken; SOH is calculated using Formulas (6) and (7) according to [30].In Formulas (6), SOC 1 and SOC 2 are the SOCs of the battery before and after discharging, respectively, and they are calculated only when the OCV is accurate enough.
where SOC(t 0 ) is the SOC at the beginning of the discharge or charge; I(t) is the working current; C M is the current capacity of the battery; C N is the brand new capacity of the battery; Passed Charge is the total capacity flow out the battery; and SOC 1 and SOC 2 are the SOCs of the battery before and after discharging, respectively.Many battery models in the current literature do not consider thermal influence.However, Huria et al. [21] added a thermal model to the equivalent circuit model and creatively progressed the Simscape model.As shown in Figure 4, the Simscape model contains three input terminals (environment temperature, positive, and negative) and three output terminals (battery temperature, SOC, and energy consumption power).
Energies 2017, 10, 404 5 of 13 These outputs can be used to analyze the working state of the cell and manage energy flow in PCCPs.The Simscape model is suitable for Matlab simulation and analysis.The parameters of the equivalent circuit model can be estimated through the least square method with experimental data, which can be obtained by pulse discharge process under different temperatures and SOCs.Then, the parameter results are incorporated into lookup tables that are related to SOC and temperature (T).The computer determines the parameters of the model via the look-up tables.This method is very practical.Moreover, the Simscape model has beentested by Wang and Huria et al. [14,21].Their results indicated that the model responds accurately in a given battery voltage fluctuation.

The Simulink Model for Parallel-Cell Pack
Given the electrical connection of traditional PCCP cells, controlling the current distribution among the cells is difficult.Moreover, monitoring the current of every cell is impractical.Thus, some of the cellsare likely to exceed the current threshold, although the average current is in a reasonable range.Therefore, it is necessary to establish a PCCP model to investigate PCCP management.Some researchers [15,31] proposed a parallel-connected cell pack model on the assumption that all the batteries are the same, which is obviously impractical.Currently, most PCCP models donot consider thermal influence, thus making most models imprecise.This study utilized the first-order RC equivalent model to build a model of eight parallel-connected cell packs.Cell parameterswere estimated from eight different aging conditions of LiNixMnyCozO2 (NMC) cells.Thermal influence was considered in the estimation process [21].
As shown in Figure 5, the PCCP model contains eight cells that are connected to the switch and direct current (DC) bus.There is no loop current among the cells with different SOCs when the switches are turned off.MATLAB S-function is used to calculate the number of cells to be connected in the circuit on the basis of power demand.Then, the cells are sorted by SOC value from high to low.The first N cells are then connected to the circuit.All the cells will be in the disconnected state under no-load conditions to completely avoid loop current.Hsu et al. [16] also recommended that all cells should be disconnected.The automobile dynamics model calculates the required current at theworking condition.The fuzzy control algorithm subsystem and MATLAB S-function sort the SOCs and select which switches to close.With the PCCP model, we can conduct the vehicle dynamic driving cycle test, simulate the depth discharge of cells, and analyze the resistance of the isolated cell, as well as its influence on the pack.

The Simulink Model for Parallel-Cell Pack
Given the electrical connection of traditional PCCP cells, controlling the current distribution among the cells is difficult.Moreover, monitoring the current of every cell is impractical.Thus, some of the cellsare likely to exceed the current threshold, although the average current is in a reasonable range.Therefore, it is necessary to establish a PCCP model to investigate PCCP management.Some researchers [15,31] proposed a parallel-connected cell pack model on the assumption that all the batteries are the same, which is obviously impractical.Currently, most PCCP models donot consider thermal influence, thus making most models imprecise.This study utilized the first-order RC equivalent model to build a model of eight parallel-connected cell packs.Cell parameterswere estimated from eight different aging conditions of LiNi x Mn y Co z O 2 (NMC) cells.Thermal influence was considered in the estimation process [21].
As shown in Figure 5, the PCCP model contains eight cells that are connected to the switch and direct current (DC) bus.There is no loop current among the cells with different SOCs when the switches are turned off.MATLAB S-function is used to calculate the number of cells to be connected in the circuit on the basis of power demand.Then, the cells are sorted by SOC value from high to low.The first N cells are then connected to the circuit.All the cells will be in the disconnected state under no-load conditions to completely avoid loop current.Hsu et al. [16] also recommended that all cells should be disconnected.The automobile dynamics model calculates the required current at theworking condition.The fuzzy control algorithm subsystem and MATLAB S-function sort the SOCs and select which switches to close.With the PCCP model, we can conduct the vehicle dynamic driving cycle test, simulate the depth discharge of cells, and analyze the resistance of the isolated cell, as well as its influence on the pack.

Principle and Process of PCCEM Strategy
The basic principles of selecting cells for battery packs are to reduce the circulating current; only the cell with the highest SOC should be selected to discharge when power demand is small; when the demand increases, the cell is successively connected with the highest SOC among the rest cells.The PCCEM flow chart in Figure 6 introduces methods to control the current discharge of each cell.First, the BMS measures the voltage, current, and temperature of batteries.SOCs are then calculated and sorted, and power demand is classified into different levels.The PCCEM algorithm determines the number (N) of cells to be connected into the circuit and ΔSOC in accordance with cell load level and SOC.After this step, the algorithm will turn on the first N cell and turn off others to drive the load and loop infinitely.When the load is small, only cells with higher SOC are connected into the circuit.At this time, cells with lower SOC have no electrical connection with cells with higher SOC.Hence, no loop current exists between the cells and the BMS estimates the cell SOC and SOH accurately under open loop conditions.As the load increases, more cells gradually access the circuit in SOC sequence.When the power demand is relatively high, the discharge current will also increase.We should connect all the cells in the circuit to reduce energy loss because energy dissipation is proportional to the square of the current (Q = I 2 × R), and connecting all the cells into the circuit will reduce overall internal resistance.When the pack SOC (generally, the minimum value of cell SOC is considered as the pack SOC) is low, we should appropriately increase the number of cells in the circuit even in light load to reduce internal resistance and prevent voltage drop to reach the cut-off voltage early.Conversely, the ideal situation in the charging mode is to charge cells with lower SOC first, followed by cells with higher SOC, and finally, cells with the highest SOC.This order will increase the energy utilization rate.Thus, the cell's status and load determines the accessibility of the cell to the circuit.The PCCEM system also can isolate poorly performing cells by turning off the corresponding switches.Then, the system will send a message to the owner.
As mentioned above, the control switches of several cells will quickly turn on or off to control the discharging current of corresponding cells when cell load and SOC change.However, a problem arises when the SOC order changes during the discharge process: any tiny changes in the SOC order will make the control frequently unit turn on or off the switches that have corresponding cells with similar SOCs.

Principle and Process of PCCEM Strategy
The basic principles of selecting cells for battery packs are to reduce the circulating current; only the cell with the highest SOC should be selected to discharge when power demand is small; when the demand increases, the cell is successively connected with the highest SOC among the rest cells.The PCCEM flow chart in Figure 6 introduces methods to control the current discharge of each cell.First, the BMS measures the voltage, current, and temperature of batteries.SOCs are then calculated and sorted, and power demand is classified into different levels.The PCCEM algorithm determines the number (N) of cells to be connected into the circuit and ∆SOC in accordance with cell load level and SOC.After this step, the algorithm will turn on the first N cell and turn off others to drive the load and loop infinitely.When the load is small, only cells with higher SOC are connected into the circuit.At this time, cells with lower SOC have no electrical connection with cells with higher SOC.Hence, no loop current exists between the cells and the BMS estimates the cell SOC and SOH accurately under open loop conditions.As the load increases, more cells gradually access the circuit in SOC sequence.When the power demand is relatively high, the discharge current will also increase.We should connect all the cells in the circuit to reduce energy loss because energy dissipation is proportional to the square of the current (Q = I 2 × R), and connecting all the cells into the circuit will reduce overall internal resistance.When the pack SOC (generally, the minimum value of cell SOC is considered as the pack SOC) is low, we should appropriately increase the number of cells in the circuit even in light load to reduce internal resistance and prevent voltage drop to reach the cut-off voltage early.Conversely, the ideal situation in the charging mode is to charge cells with lower SOC first, followed by cells with higher SOC, and finally, cells with the highest SOC.This order will increase the energy utilization rate.Thus, the cell's status and load determines the accessibility of the cell to the circuit.The PCCEM system also can isolate poorly performing cells by turning off the corresponding switches.Then, the system will send a message to the owner.
As mentioned above, the control switches of several cells will quickly turn on or off to control the discharging current of corresponding cells when cell load and SOC change.However, a problem arises when the SOC order changes during the discharge process: any tiny changes in the SOC order Energies 2017, 10, 404 7 of 13 will make the control frequently unit turn on or off the switches that have corresponding cells with similar SOCs.To avoid this issue, the frequency of SOC sorting should be controlled within a reasonable range.A variable ΔSOC is introduced, which is a threshold value.Only when the SOC of any of the cell changes beyond the ΔSOC will the PCCEM reorder the SOCs.As the discharge current increases, the SOC changes more quickly.Thus, increasing ΔSOC is advantageous because this will reduce the on-off frequency of the switches, specifically the ΔSOC changes with the discharge current.Under heavy load, all cells will be connected to the bus.Controlling the current distribution is unnecessary at this time.

Fuzzy Logic Control Strategy
Traditional control theory is powerful in a linear system, but is difficult to control in a lithium battery, which is a nonlinear system [32,33].FLC can simulate human knowledge and reasoning and is robust even when measurements are imprecise or when input characteristics change [34][35][36].To explore the influence of different factors, such as power demand and discharge current, this study uses fuzzy control strategy in the PCCEM strategy.FLC calculates the optimal number of cells connected in the circuit and the value of ΔSOC.
The two FLC inputs are cell SOC and demand current.These inputs can be quantified into different grades and expressed with linguistic variables, such as small and medium.The membership functions of discharge current and battery SOC, which are based on extensive experimental data and theoretical analysis, are indicated in Figure 7. Figure 7a shows the membership function of demand current.The fuzzy subset L, M, H signifies that the current demand is low, medium, or high, respectively.Figure 7b shows the membership functions of pack SOC, whereas the fuzzy subset L, M, and H shows the battery pack SOC as insufficient, medium, and enough, respectively.To avoid this issue, the frequency of SOC sorting should be controlled within a reasonable range.A variable ∆SOC is introduced, which is a threshold value.Only when the SOC of any of the cell changes beyond the ∆SOC will the PCCEM reorder the SOCs.As the discharge current increases, the SOC changes more quickly.Thus, increasing ∆SOC is advantageous because this will reduce the on-off frequency of the switches, specifically the ∆SOC changes with the discharge current.Under heavy load, all cells will be connected to the bus.Controlling the current distribution is unnecessary at this time.

Fuzzy Logic Control Strategy
Traditional control theory is powerful in a linear system, but is difficult to control in a lithium battery, which is a nonlinear system [32,33].FLC can simulate human knowledge and reasoning and is robust even when measurements are imprecise or when input characteristics change [34][35][36].To explore the influence of different factors, such as power demand and discharge current, this study uses fuzzy control strategy in the PCCEM strategy.FLC calculates the optimal number of cells connected in the circuit and the value of ∆SOC.
The two FLC inputs are cell SOC and demand current.These inputs can be quantified into different grades and expressed with linguistic variables, such as small and medium.The membership functions of discharge current and battery SOC, which are based on extensive experimental data and theoretical analysis, are indicated in Figure 7. Figure 7a shows the membership function of demand current.The fuzzy subset L, M, H signifies that the current demand is low, medium, or high, respectively.Figure 7b shows the membership functions of pack SOC, whereas the fuzzy subset L, M, and H shows the battery pack SOC as insufficient, medium, and enough, respectively.To avoid this issue, the frequency of SOC sorting should be controlled within a reasonable range.A variable ΔSOC is introduced, which is a threshold value.Only when the SOC of any of the cell changes beyond the ΔSOC will the PCCEM reorder the SOCs.As the discharge current increases, the SOC changes more quickly.Thus, increasing ΔSOC is advantageous because this will reduce the on-off frequency of the switches, specifically the ΔSOC changes with the discharge current.Under heavy load, all cells will be connected to the bus.Controlling the current distribution is unnecessary at this time.

Fuzzy Logic Control Strategy
Traditional control theory is powerful in a linear system, but is difficult to control in a lithium battery, which is a nonlinear system [32,33].FLC can simulate human knowledge and reasoning and is robust even when measurements are imprecise or when input characteristics change [34][35][36].To explore the influence of different factors, such as power demand and discharge current, this study uses fuzzy control strategy in the PCCEM strategy.FLC calculates the optimal number of cells connected in the circuit and the value of ΔSOC.
The two FLC inputs are cell SOC and demand current.These inputs can be quantified into different grades and expressed with linguistic variables, such as small and medium.The membership functions of discharge current and battery SOC, which are based on extensive experimental data and theoretical analysis, are indicated in Figure 7. Figure 7a shows the membership function of demand current.The fuzzy subset L, M, H signifies that the current demand is low, medium, or high, respectively.Figure 7b shows the membership functions of pack SOC, whereas the fuzzy subset L, M, and H shows the battery pack SOC as insufficient, medium, and enough, respectively.The objective of managing the discharge current is to make every cell work with high efficiency and safety.The rules for management, which are based on experience and experimental results, are summarized in Figure 8. Figure 8a shows that the number (N) of cells connected to the circuit increases with increasing discharge current, but decreases with increasing SOC. Figure 8b shows that ∆SOC, which is the interval to reorder SOCs, increases with increasing discharge current and with increasing SOC.
Energies 2017, 10, 404 8 of 13 The objective of managing the discharge current is to make every cell work with high efficiency and safety.The rules for management, which are based on experience and experimental results, are summarized in Figure 8. Figure 8a shows that the number (N) of cells connected to the circuit increases with increasing discharge current, but decreases with increasing SOC. Figure 8b shows that ΔSOC, which is the interval to reorder SOCs, increases with increasing discharge current and with increasing SOC.

Discussion of Simulation Results
To verify the control effect of PCCEM, a NEDC simulation experiment is conducted with the PCCP model presented in Figure 5.The parameters of the batteries in the PCCP were obtained from eight different aged batteries (power-oriented, max voltage 4.2 V, 31 Ah lithium ion cells of the LiNixMnyCozO2 chemistry) using the parameter identification method in [21].The initial distributions in SOH and capacity among cells in the PCCP are listed as Table 1.The load of the NEDC cycle is a vehicle simulation model, in which mass M = 3000 kg, road dragging coefficient Ψ = 0.0076 + 0.000056 × uv, correction coefficient of rotating mass δ = 1.2, coefficient of air resistance CD = 0.3, and windward area A = 2.1 m 2 .
Figure 9a shows the profiles of vehicle speed and pack discharge current.The idle condition during 0-1000 s and 2117-3000 s highlights the performance of PCCEM.

Discussion of Simulation Results
To verify the control effect of PCCEM, a NEDC simulation experiment is conducted with the PCCP model presented in Figure 5.The parameters of the batteries in the PCCP were obtained from eight different aged batteries (power-oriented, max voltage 4.2 V, 31 Ah lithium ion cells of the LiNi x Mn y Co z O 2 chemistry) using the parameter identification method in [21].The initial distributions in SOH and capacity among cells in the PCCP are listed as Table 1.The load of the NEDC cycle is a vehicle simulation model, in which mass M = 3000 kg, road dragging coefficient Ψ = 0.0076 + 0.000056 × u v , correction coefficient of rotating mass δ = 1.2, coefficient of air resistance C D = 0.3, and windward area A = 2.1 m 2 .
Figure 9a shows the profiles of vehicle speed and pack discharge current.The idle condition during 0-1000 s and 2117-3000 s highlights the performance of PCCEM.
Energies 2017, 10, 404 8 of 13 The objective of managing the discharge current is to make every cell work with high efficiency and safety.The rules for management, which are based on experience and experimental results, are summarized in Figure 8. Figure 8a shows that the number (N) of cells connected to the circuit increases with increasing discharge current, but decreases with increasing SOC. Figure 8b shows that ΔSOC, which is the interval to reorder SOCs, increases with increasing discharge current and with increasing SOC.

Discussion of Simulation Results
To verify the control effect of PCCEM, a NEDC simulation experiment is conducted with the PCCP model presented in Figure 5.The parameters of the batteries in the PCCP were obtained from eight different aged batteries (power-oriented, max voltage 4.2 V, 31 Ah lithium ion cells of the LiNixMnyCozO2 chemistry) using the parameter identification method in [21].The initial distributions in SOH and capacity among cells in the PCCP are listed as Table 1.The load of the NEDC cycle is a vehicle simulation model, in which mass M = 3000 kg, road dragging coefficient Ψ = 0.0076 + 0.000056 × uv, correction coefficient of rotating mass δ = 1.2, coefficient of air resistance CD = 0.3, and windward area A = 2.1 m 2 .
Figure 9a shows the profiles of vehicle speed and pack discharge current.The idle condition during 0-1000 s and 2117-3000 s highlights the performance of PCCEM.  Figure 10a shows the SOCs without PCCEM strategy.Figure 10b shows the SOCs with PCCEM.During the driving cycles, the SOC of the eight cells changed in different ways.During 0-1000 s and 2117-3000 s, cells with higher SOC (Cells 1-4) are discharged, whereas the others (Cells 5-8) are charged.Loop current is shown in Figure 11a.All SOCs remain unchanged while the PCCEM strategy is on.Loop current is close to zero in Figure 11b.During 1000-2117 s, cells with a high SOC simultaneously supply the current for the load and charge cells with a low SOC, as shown in Figure 10a.In Figure 10b, the cells discharge in the SOC order and only cells with a higher SOC supply power when in light load.Then, cells with a low SOC are disconnected from the circuit.Therefore, loop current decreases significantly.
Energies 2017, 10, 404 9 of 13 Figure 10a shows the SOCs without PCCEM strategy.Figure 10b shows the SOCs with PCCEM.During the driving cycles, the SOC of the eight cells changed in different ways.During 0-1000 s and 2117-3000 s, cells with higher SOC (Cells 1-4) are discharged, whereas the others (Cells 5-8) are charged.Loop current is shown in Figure 11a.All SOCs remain unchanged while the PCCEM strategy is on.Loop current is close to zero in Figure 11b.During 1000-2117 s, cells with a high SOC simultaneously supply the current for the load and charge cells with a low SOC, as shown in Figure 10a.In Figure 10b, the cells discharge in the SOC order and only cells with a higher SOC supply power when in light load.Then, cells with a low SOC are disconnected from the circuit.Therefore, loop current decreases significantly.Figure 11a shows the cell current without PCCEM strategy.Figure 11b shows the cell current with PCCEM.During 0-1000 s and 2117-3000 s, the demand power is zero and the pack does not need to discharge, as shown in Figure 9b.However, Figure 11a shows that loop current exists if a discrepancy in the SOC occurs in the parallel cells.Figure 11b clearly shows that with the PCCEM strategy, all the cells are disconnected from the bus when the power demand is zero.Hence, the current of each cell is zero, and no loop current will occur.Figure 12a shows the cell temperature without PCCEM strategy, and Figure 12b shows the cell temperature with PCCEM.The power dissipated inside the cell (Ps) can be calculated as: where I and I1 are the current through the resistors R and R1, respectively.The meaning of R0 and R1 are presented in Figure 3. Cell temperature is computed by solving the heat equation of a homogeneous body that is exchanging heat with its environment [19].Comparing Figure 12a,b reveals that the temperature of each cell is different from each other even in idle mode, given loop current.When PCCEM is on at Figure 11a shows the cell current without PCCEM strategy.Figure 11b shows the cell current with PCCEM.During 0-1000 s and 2117-3000 s, the demand power is zero and the pack does not need to discharge, as shown in Figure 9b.However, Figure 11a shows that loop current exists if a discrepancy in the SOC occurs in the parallel cells.Figure 11b clearly shows that with the PCCEM strategy, all the cells are disconnected from the bus when the power demand is zero.Hence, the current of each cell is zero, and no loop current will occur.
Energies 2017, 10, 404 9 of 13 Figure 10a shows the SOCs without PCCEM strategy.Figure 10b shows the SOCs with PCCEM.During the driving cycles, the SOC of the eight cells changed in different ways.During 0-1000 s and 2117-3000 s, cells with higher SOC (Cells 1-4) are discharged, whereas the others (Cells 5-8) are charged.Loop current is shown in Figure 11a.All SOCs remain unchanged while the PCCEM strategy is on.Loop current is close to zero in Figure 11b.During 1000-2117 s, cells with a high SOC simultaneously supply the current for the load and charge cells with a low SOC, as shown in Figure 10a.In Figure 10b, the cells discharge in the SOC order and only cells with a higher SOC supply power when in light load.Then, cells with a low SOC are disconnected from the circuit.Therefore, loop current decreases significantly.Figure 11a shows the cell current without PCCEM strategy.Figure 11b shows the cell current with PCCEM.During 0-1000 s and 2117-3000 s, the demand power is zero and the pack does not need to discharge, as shown in Figure 9b.However, Figure 11a shows that loop current exists if a discrepancy in the SOC occurs in the parallel cells.Figure 11b clearly shows that with the PCCEM strategy, all the cells are disconnected from the bus when the power demand is zero.Hence, the current of each cell is zero, and no loop current will occur.Figure 12a shows the cell temperature without PCCEM strategy, and Figure 12b shows the cell temperature with PCCEM.The power dissipated inside the cell (Ps) can be calculated as: where I and I1 are the current through the resistors R and R1, respectively.The meaning of R0 and R1 are presented in Figure 3.
Cell temperature is computed by solving the heat equation of a homogeneous body that is exchanging heat with its environment [19].Comparing Figure 12a,b reveals that the temperature of each cell is different from each other even in idle mode, given loop current.When PCCEM is on at   Figure 12a shows the cell temperature without PCCEM strategy, and Figure 12b shows the cell temperature with PCCEM.The power dissipated inside the cell (P s ) can be calculated as: Energies 2017, 10, 404 10 of 13 where I and I 1 are the current through the resistors R and R 1 , respectively.The meaning of R 0 and R 1 are presented in Figure 3. Cell temperature is computed by solving the heat equation of a homogeneous body that is exchanging heat with its environment [19].Comparing Figure 12a,b reveals that the temperature of each cell is different from each other even in idle mode, given loop current.When PCCEM is on at the end of NEDC, the highest temperature is 20.47 • C, which is lower than 20.59 • C in Figure 12a.Thus, the temperature difference decreases.
The consumed energy during the NEDC cycle can be calculated by where u m (t) is the voltage of the cell and i m (t) is the current flow through the cell.
The dissipated energy per cell can be calculated via the time integration of P s .Then, adding the values yields dissipated energy.
Our calculations show that the energy consumed by the NEDC cycle is 24.75 kJ.As shown in Figure 13, the dissipated energy is 2.145 kJ, which accounts for 8.66% of the consumed energy.The distinction of dissipated energy is very clear at the beginning of the discharge circle.When PCCEM is on, the circuit current almost disappears, especially at 0-1000 s, where the dissipated energy is zero and the PCCEM of the experiment is approximately 1.171 kJ.At the end of this working condition, the PCCEM strategy reduces the dissipated energy from 2.145 to 0.571 kJ, which indicates that 73.3% is cut off.Namely, during the NEDC cycle, 6.36% of the consumed energy is avoided by the PCCEM.
Energies 2017, 10, 404 10 of 13 the end of NEDC, the highest temperature is 20.47 °C, which is lower than 20.59 °C in Figure 12a.Thus, the temperature difference decreases.The consumed energy during the NEDC cycle can be calculated by where um(t) is the voltage of the cell and im(t) is the current flow through the cell.
The dissipated energy per cell can be calculated via the time integration of Ps.Then, adding the values yields dissipated energy.
Our calculations show that the energy consumed by the NEDC cycle is 24.75 kJ.As shown in Figure 13, the dissipated energy is 2.145 kJ, which accounts for 8.66% of the consumed energy.The distinction of dissipated energy is very clear at the beginning of the discharge circle.When PCCEM is on, the circuit current almost disappears, especially at 0-1000 s, where the dissipated energy is zero and the PCCEM of the experiment is approximately 1.171 kJ.At the end of this working condition, the PCCEM strategy reduces the dissipated energy from 2.145 to 0.571 kJ, which indicates that 73.3% is cut off.Namely, during the NEDC cycle, 6.36% of the consumed energy is avoided by the PCCEM.In the simulation, the closed resistance of each switch is 0.01 Ohm, and the open conductance is 10 -8 1/Ohm.The switching frequency was restricted to very small scales; thus, we calculated for the power consumed by the switches closed resistance and the current flow through it.The formula is: The result is shown in Figure 14.When t = 3000 s, Es = 0.597 kJ when the power consumed by the  the end of NEDC, the highest temperature is 20.47 °C, which is lower than 20.59 °C in Figure 12a.Thus, the temperature difference decreases.The consumed energy during the NEDC cycle can be calculated by where um(t) is the voltage of the cell and im(t) is the current flow through the cell.
The dissipated energy per cell can be calculated via the time integration of Ps.Then, adding the values yields dissipated energy.
Our calculations show that the energy consumed by the NEDC cycle is 24.75 kJ.As shown in Figure 13, the dissipated energy is 2.145 kJ, which accounts for 8.66% of the consumed energy.The distinction of dissipated energy is very clear at the beginning of the discharge circle.When PCCEM is on, the circuit current almost disappears, especially at 0-1000 s, where the dissipated energy is zero and the PCCEM of the experiment is approximately 1.171 kJ.At the end of this working condition, the PCCEM strategy reduces the dissipated energy from 2.145 to 0.571 kJ, which indicates that 73.3% is cut off.Namely, during the NEDC cycle, 6.36% of the consumed energy is avoided by the PCCEM.In the simulation, the closed resistance of each switch is 0.01 Ohm, and the open conductance is 10 -8 1/Ohm.The switching frequency was restricted to very small scales; thus, we calculated for the power consumed by the switches closed resistance and the current flow through it.The formula is: The result is shown in Figure 14.When t = 3000 s, Es = 0.597 kJ when the power consumed by the In the simulation, the closed resistance of each switch is 0.01 Ohm, and the open conductance is 10 −8 1/Ohm.The switching frequency was restricted to very small scales; thus, we calculated for the power consumed by the switches closed resistance and the current flow through it.The formula is: The result is shown in Figure 14.When t = 3000 s, E s = 0.597 kJ when the power consumed by the switches is taken into account.The PCCEM strategy reduces the dissipated energy from 2.145 to 1.168 kJ, which indicates that 45.5% is the cut off, rather than 73.3%.Therefore, the power consumed by the switches is not negligible in the experiment.
The experimental results verify that the PCCEM algorithm can isolate poorly performing cells and arrange the discharge order of cells to effectively reduce loop current.Hence, the algorithm simultaneously reduces the pack-dissipated energy and improves the utilization rate of battery power.
Energies 2017, 10, 404 11 of 13 kJ, which indicates that 45.5% is the cut off, rather than 73.3%.Therefore, the power consumed by the switches is not negligible in the experiment.
The experimental results verify that the PCCEM algorithm can isolate poorly performing cells and arrange the discharge order of cells to effectively reduce loop current.Hence, the algorithm simultaneously reduces the pack-dissipated energy and improves the utilization rate of battery power.

Conclusions and Prospects
This study built a Simscape PCCP model that considers thermal effect.The PCCEM strategy is proposed on the basis of the PCCP model and fuzzy control logic.The PCCEM strategy automatically adjusts the number of cells in the circuit in accordance with load demand and turns on the first N switches in the corresponding SOC order.This method can be easily utilized in a single parallel string.
For a pack with many parallel strings, for example, 64 cells, PCCEM can treat every eight cells as a group.Thus, eight groups were obtained.If we parallel eight groups and manage these groups with PCCEM, then we are managing a big group with 64 cells using PCCEM.The PCCP model is designed especially for aging cells, and the reuse of Li-ion EV batteries for energy storage systems in stationary settings is a promising technology to support improved management of the demand and supply of electricity [37].Reused batteries from electric vehicles, hybrid electric vehicles, and plug-in hybrid electric vehicles present an excellent, cost-effective option for energy storage applications that can help build 'smart grid' technologies, such as computer-based remote control, automation, and information management, to improve the reliability, efficiency and economics of electricity generation and distribution (U.S.Department of Energy, 2014).The PCCP model can be used to control the discharge current of reused batteries to reduce loop current among them.However, the model has many limitations; for example, each cell needs a control switch and a current sensor, which will increase the cost.In addition, control in large-scale applications will become extremely complicated.This study is only a preliminary study.We believe that more work can be done with the PCCP model to manage the measurement of PCCP energy in the future.

Conclusions and Prospects
This study built a Simscape PCCP model that considers thermal effect.The PCCEM strategy is proposed on the basis of the PCCP model and fuzzy control logic.The PCCEM strategy automatically adjusts the number of cells in the circuit in accordance with load demand and turns on the first N switches in the corresponding SOC order.This method can be easily utilized in a single parallel string.
For a pack with many parallel strings, for example, 64 cells, PCCEM can treat every eight cells as a group.Thus, eight groups were obtained.If we parallel eight groups and manage these groups with PCCEM, then we are managing a big group with 64 cells using PCCEM.The PCCP model is designed especially for aging cells, and the reuse of Li-ion EV batteries for energy storage systems in stationary settings is a promising technology to support improved management of the demand and supply of electricity [37].Reused batteries from electric vehicles, hybrid electric vehicles, and plug-in hybrid electric vehicles present an excellent, cost-effective option for energy storage applications that can help build 'smart grid' technologies, such as computer-based remote control, automation, and information management, to improve the reliability, efficiency and economics of electricity generation and distribution (U.S.Department of Energy, 2014).The PCCP model can be used to control the discharge current of reused batteries to reduce loop current among them.However, the model has many limitations; for example, each cell needs a control switch and a current sensor, which will increase the cost.In addition, control in large-scale applications will become extremely complicated.This study is only a preliminary study.We believe that more work can be done with the PCCP model to manage the measurement of PCCP energy in the future.

Figure 1 .
Figure 1.Inconsistency in cells' capacities and states of change (SOCs).

Figure 1 .
Figure 1.Inconsistency in cells' capacities and states of change (SOCs).

Figure 3 .
Figure 3. First-order RC equivalent circuit model.Figure 3. First-order RC equivalent circuit model.

Figure 3 .
Figure 3. First-order RC equivalent circuit model.Figure 3. First-order RC equivalent circuit model.
, the Simscape model has beentested by Wang and Huria et al.[14,21].Their results indicated that the model responds accurately in a given battery voltage fluctuation.

Figure 4 .
Figure 4. Simscape model of the first-order RC equivalent circuit.

Figure 4 .
Figure 4. Simscape model of the first-order RC equivalent circuit.

Figure 5 .
Figure 5. Simulink model of the eight parallel-cell pack.

Figure 5 .
Figure 5. Simulink model of the eight parallel-cell pack.

Figure 7 .Figure 6 .
Figure 7. Membership functions of inputs.(a) Membership of discharge current; and (b) membership of battery SOC.

Figure 7 .Figure 7 .
Figure 7. Membership functions of inputs.(a) Membership of discharge current; and (b) membership of battery SOC.

Figure 8 .
Figure 8. Decision surface of fuzzy logic controller.(a) The membership of N; and (b) the membership of SOC.

Figure 9 .
Figure 9.The profiles of speed and discharge current.(a) Speed profile for the experiment; and (b) discharge current profile of the packs.

Figure 8 .
Figure 8. Decision surface of fuzzy logic controller.(a) The membership of N; and (b) the membership of SOC.

Figure 8 .
Figure 8. Decision surface of fuzzy logic controller.(a) The membership of N; and (b) the membership of SOC.

Figure 9 .
Figure 9.The profiles of speed and discharge current.(a) Speed profile for the experiment; and (b) discharge current profile of the packs.

Figure 9 .
Figure 9.The profiles of speed and discharge current.(a) Speed profile for the experiment; and (b) discharge current profile of the packs.

Figure 10 .
Figure 10.Comparison of the SOCs.(a) SOC of each cell without PCCEM; and (b) SOC of each cell with PCCEM.

Figure 11 .
Figure 11.Comparison of loop currents.(a) Current of each cell without PCCEM; and (b) current of each cell with PCCEM.

Figure 10 .
Figure 10.Comparison of the SOCs.(a) SOC of each cell without PCCEM; and (b) SOC of each cell with PCCEM.

Figure 10 .
Figure 10.Comparison of the SOCs.(a) SOC of each cell without PCCEM; and (b) SOC of each cell with PCCEM.

Figure 11 .
Figure 11.Comparison of loop currents.(a) Current of each cell without PCCEM; and (b) current of each cell with PCCEM.

Figure 11 .
Figure 11.Comparison of loop currents.(a) Current of each cell without PCCEM; and (b) current of each cell with PCCEM.

Figure 12 .
Figure 12.Comparison of cell temperatures.(a) Temperature of each cell without PCCEM; and (b) temperature of each cell with PCCEM.

Figure 12 .
Figure 12.Comparison of cell temperatures.(a) Temperature of each cell without PCCEM; and (b) temperature of each cell with PCCEM.

Figure 12 .
Figure 12.Comparison of cell temperatures.(a) Temperature of each cell without PCCEM; and (b) temperature of each cell with PCCEM.

Table 1 .
Initial distributions in SOH and capacity among cells.

Table 1 .
Initial distributions in SOH and capacity among cells.

Table 1 .
Initial distributions in SOH and capacity among cells.
Comparison of dissipated energy.