{'fraction_non_alphanumeric': 0.025656659268833988, 'fraction_numerical': 0.010413872155695924, 'mean_word_length': 4.703714007511476, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Research on the Optimal Charging Strategy for Li-Ion Batteries Based on Multi-Objective Optimization

Charging performance affects the commercial application of electric vehicles (EVs) significantly. This paper presents an optimal charging strategy for Li-ion batteries based on the voltage-based multistage constant current (VMCC) charging strategy. In order to satisfy the different charging demands of the EV users for charging time, charged capacity and energy loss, the multi-objective particle swarm optimization (MOPSO) algorithm is employed and the influences of charging stage number, charging cut-off voltage and weight factors of different charging goals are analyzed. Comparison experiments of the proposed charging strategy and the traditional normal and fast charging strategies are carried out. The experimental results demonstrate that the traditional normal and fast charging strategies can only satisfy a small range of EV users’ charging demand well while the proposed charging strategy can satisfy the whole range of the charging demand well. The relative increase in charging performance of the proposed charging strategy can reach more than 80% when compared to the normal and fast charging dependently.


Introduction
Electric vehicles have been the new development point of the motor industry because of the advantages of zero pollution and high energy efficiency [1][2][3]. Li-ion batteries with high energy density, good life cycle performance, and low self-discharge rate [4][5][6][7] are being more widely used in electric vehicles. However, the Li-ion battery charging strategy has become the bottleneck of their application, due to the slow charging speed and uncertain effects on battery life. Li-ion battery charging optimization technology has become one of the central difficult issues [8].
Many studies have researched the charging techniques and topologies, such as constant current (CC), constant voltage (CV), and constant current constant voltage (CC-CV). One of the most commonly used charging strategy is the CC-CV charging [9][10][11]. In a CC-CV charging process, the battery is first charged with constant current, and when the battery voltage reaches the cut-off value, the battery is charged with constant voltage. However, the constant voltage process extends the charging time. Much effort has been made to carry out a better charging algorithm. Gao [12] proposed a boundary-curve based optimum charging current strategy. The temperature rise and polarization voltage were used to obtain the boundary curves during the charging process. The charging time was decreased and battery capacity was increased. In [13], a constant polarization-based fuzzy-control charging pattern was proposed, while charging strategy was shortened, no obvious temperature rise was observed. Chen [14] proposed a duty-varied voltage charging strategy that can dynamically detect and track the suitable duty of the charging pulse.
Multistage constant current (MCC) charging is another widely researched charging method, which has the advantage of faster charging speed and better battery charging performance [15][16][17][18][19][20]. As illustrated in Figure 1, the entire charging process of MCC charging consists of a certain number of constant current charging processes. When the preset switching condition is reached, the charging process moves into the next charging stage until all the charging stages are completed. According to the switching condition, the MCC charging strategy can be divided into two categories: voltage-based MCC (VMCC) charging and state of charge (SOC)-based MCC (SMCC) charging. The research on MCC charging mostly focuses on the optimization of the charging current of each stage. In [17], a Taguchi-based algorithm is presented and orthogonal arrays are implemented to determine the optimal charging current pattern. Another paper [18] combined the Taguchi method and battery SOC estimation technology and proposed an SOC-based MCC charging strategy. Wang [19] regarded the problem as a combinatorial optimization problem, and proposed a PSO-based fuzzy-controlled method to search for the optimal pattern of the charging currents. In [20], an Ant Colony System algorithm was used to select the optimum charging current among five charging stages. In [21], the whole charging process was divided into ten parts according to SOC, and the ILP algorithm is adopted to solve the optimization problem. The charging time, charged capacity and energy loss are the most common concerns regarding performance in battery charging [22,23]. The traditional EV charging strategy can mainly be classified as normal charging or fast charging. For normal charging, the charging time can be hours and is usually used for an overnight charge, and fast charging is usually chosen when there is a demand for a short charging time. However, in daily uses, people have different charging requirements for overnight charging and fast charging. An optimal charging strategy should be able to provide different charging currents according to the different charging demands of the users. As a main contribution of this paper, an optimal charging strategy that can satisfy different charging demands of the EV users is proposed based on the MOPSO algorithm. Also, a simulation model of the VMCC charging algorithm is proposed to simulate the charging process of the VMCC charging strategy. As another contribution of this paper, the influence of charging parameters such as the charging cut-off voltage and the number of charging stages are discussed.
The rest of this paper is organized as follows. Section 2 establishes the Li-ion battery equivalent circuit model (ECM) and thermal model, studies the estimation of the battery model parameter and proposes the simulation model of the VMCC charging. Section 3 introduces the experimental setup and conducts the battery parameter estimation experiment. Section 4 formulates the multi-objective optimization problem. The MOPSO algorithm and the searching procedure are described in detail in Section 5. In Section 6, the influence of charging stage number, charging cut-off voltage and weight The charging time, charged capacity and energy loss are the most common concerns regarding performance in battery charging [22,23]. The traditional EV charging strategy can mainly be classified as normal charging or fast charging. For normal charging, the charging time can be hours and is usually used for an overnight charge, and fast charging is usually chosen when there is a demand for a short charging time. However, in daily uses, people have different charging requirements for overnight charging and fast charging. An optimal charging strategy should be able to provide different charging currents according to the different charging demands of the users. As a main contribution of this paper, an optimal charging strategy that can satisfy different charging demands of the EV users is proposed based on the MOPSO algorithm. Also, a simulation model of the VMCC charging algorithm is proposed to simulate the charging process of the VMCC charging strategy. As another contribution of this paper, the influence of charging parameters such as the charging cut-off voltage and the number of charging stages are discussed.
The rest of this paper is organized as follows. Section 2 establishes the Li-ion battery equivalent circuit model (ECM) and thermal model, studies the estimation of the battery model parameter and proposes the simulation model of the VMCC charging. Section 3 introduces the experimental setup and conducts the battery parameter estimation experiment. Section 4 formulates the multi-objective Energies 2017, 10, 709 3 of 15 optimization problem. The MOPSO algorithm and the searching procedure are described in detail in Section 5. In Section 6, the influence of charging stage number, charging cut-off voltage and weight factors of different charging goals are discussed and experiments are carried out to verify the performance of the proposed charging strategy. Finally, Section 7 concludes this paper.

Battery Equivalent Circuit Model
To optimize the charging currents of the VMCC charging procedure, an accurate and reliable battery model must be formulated when using model-based methods. The ECM (equivalent circuit model) approach uses the circuit elements to represent the electrochemical phenomena inside the battery to avoid the lengthy electrochemical process calculation while still capturing the dynamic performance of the battery [24,25]. In general, an ECM with one or two Resistance-Capacitance (RC) blocks without the parasitic branch is a common choice and for the batteries that have an obvious hysteresis phenomenon, a hysteresis model is also included [26]. In this paper, only the charging process is considered, so the one RC block model without the hysteresis model is chosen.
A one-RC network based ECM is shown in Figure 2. U oc is the battery open circuit voltage (OCV); U t is the battery terminal voltage; I L is the battery load current (assumed positive for discharging and negative for charging) and I p is the outflow current of C p . R r is the ohmic resistance which represents the electrical resistance in the battery during the accumulation and dissipation of charge in the electrical double-layer; U r is the voltage across R r . R p is the polarization resistance; and C p is the polarization capacity. The RC network consisted with R p and C p describes the mass transport effects and dynamic voltage performances of the battery. The polarization voltage on R r and C p is represented as U p . factors of different charging goals are discussed and experiments are carried out to verify the performance of the proposed charging strategy. Finally, Section 7 concludes this paper.

Battery Equivalent Circuit Model
To optimize the charging currents of the VMCC charging procedure, an accurate and reliable battery model must be formulated when using model-based methods. The ECM (equivalent circuit model) approach uses the circuit elements to represent the electrochemical phenomena inside the battery to avoid the lengthy electrochemical process calculation while still capturing the dynamic performance of the battery [24,25]. In general, an ECM with one or two Resistance-Capacitance (RC) blocks without the parasitic branch is a common choice and for the batteries that have an obvious hysteresis phenomenon, a hysteresis model is also included [26]. In this paper, only the charging process is considered, so the one RC block model without the hysteresis model is chosen.
A one-RC network based ECM is shown in Figure 2.
is the battery open circuit voltage (OCV); is the battery terminal voltage; is the battery load current (assumed positive for discharging and negative for charging) and is the outflow current of . is the ohmic resistance which represents the electrical resistance in the battery during the accumulation and dissipation of charge in the electrical double-layer; is the voltage across . is the polarization resistance; and is the polarization capacity. The RC network consisted with and describes the mass transport effects and dynamic voltage performances of the battery. The polarization voltage on and is represented as . According to the battery model, we have: The discrete format can be written as: where ( ) and ( ) is are values of and at the k-th step, respectively. = · is the time constant of the circuit and ∆ is the sampling time.
All the parameters of the battery model are the function of SOC, which is defined as: where denotes battery SOC at time t, 0 denotes the initial value of battery SOC (defined as 0 in this paper for the charging condition), is the rated capacity of the battery and η is the current efficiency, which is assumed to be 1 in this paper. According to the battery model, we have: The discrete format can be written as: where U p (k) and I L (k) is are values of U p and I L at the k-th step, respectively. τ = C p ·R p is the time constant of the circuit and ∆t is the sampling time.
All the parameters of the battery model are the function of SOC, which is defined as: where SOC t denotes battery SOC at time t, SOC 0 denotes the initial value of battery SOC (defined as 0 in this paper for the charging condition), C Q is the rated capacity of the battery and η is the current efficiency, which is assumed to be 1 in this paper. The energy loss of the battery can be calculated as:

Battery Parameter Estimation
The Least Square estimation method is used to estimate the battery parameter. Define the voltage difference between the terminal voltage and the OCV as: From Equation (1), U d can be calculated as: According to Equations (1) and (6), the system transform function can be calculated as: Setting the sampling time as 1 s and considering the slow variation of the battery parameter, the following relationship can be defined: The difference equation form of Equation (7) can be written as: where Define the following equation: Then, for n sets of sampled data, the following equations can be written: Then the least square method can be employed to estimate the battery parameter according to the following equation:

Battery Thermal Model
The value of the ECM components depends on SOC and battery temperature, so the battery thermal model is necessary for the prediction of battery voltage. The battery thermal model is consistent with thermal generation and thermal transfer. In this paper, the heat generation process of the battery is attributed to the ohmic heat and the inner battery temperature is assumed to be uniform. The thermal model of the battery is computed as:  (13) where Q is the heat generated from the ohmic resistance; C T = m batt ·C h,batt , m batt is the battery mass, C h,batt is the heat capacity of the battery; T is the battery temperature and T a is the ambient temperature; R E f f is the heat conduction coefficient between the battery and the ambience.
The discrete format of Equation (13) can be written as:

Simulation of the VMCC Charging Strategy
A simulation model of the VMCC charging strategy is established to simulate the charging process of the battery under the VMCC charging process. The charging time, charged capacity and energy loss of the VMCC charging process can be acquired directly from the simulation. The flow chart of the simulation is shown in Figure 3 and the pseudo code of the VMCC charging strategy are shown in Algorithm 1.

Algorithm 1.
Pseudo code of the VMCC charging strategy.

Simulation of the VMCC Charging Strategy
A simulation model of the VMCC charging strategy is established to simulate the charging process of the battery under the VMCC charging process. The charging time, charged capacity and energy loss of the VMCC charging process can be acquired directly from the simulation. The flow chart of the simulation is shown in Figure 3 and the pseudo code of the VMCC charging strategy are shown in Algorithm 1.

Experimental Setup
The experimental setup consists of (1) a LiFePO 4 test battery (Bak, Shenzhen, China); (2) a Chroma 72001 battery (Chroma ATE, Taoyuan, China) test system; (3) a host computer; and (4) a thermal chamber (ENVAN, Shanghai, China). The nominal voltage of the battery is 3.2 V, the recommended charging voltage is 3.65 V. The nominal capacity is 2.5 Ah. The Chroma 72001 battery test system is used to charge and discharge the battery with a maximum voltage of 5 V and a maximum current of 20 A with an accuracy of 1 mV voltage and 10 mA current. The thermal chamber is used to maintain the environmental temperature of the experiments. The experimental data such as current and voltage are measured by the Chroma 72001 and recorded by the host computer.

Battery Parameter Estimation
The pulse current charging profile is executed in 20 • C, 25 • C, 30 • C, 35 • C and 40 • C respectively to estimate the battery charging parameter under different temperature. The current pattern and corresponding battery terminal voltage under the temperature of 25 • C is shown in Figure 3a. Figure 3b shows the comparison of the voltage response between the experiment and simulation. The absolute error is also shown in Figure 3b. It can be seen that the simulated results coincide well with the experimental results. Figure 4 is the estimated battery parameter under the temperature of 25 • C.

Experimental Setup
The experimental setup consists of (1) a LiFePO4 test battery (Bak, Shenzhen, China); (2) a Chroma 72001 battery (Chroma ATE, Taoyuan, China) test system; (3) a host computer; and (4) a thermal chamber (ENVAN, Shanghai, China). The nominal voltage of the battery is 3.2 V, the recommended charging voltage is 3.65 V. The nominal capacity is 2.5 Ah. The Chroma 72001 battery test system is used to charge and discharge the battery with a maximum voltage of 5 V and a maximum current of 20 A with an accuracy of 1 mV voltage and 10 mA current. The thermal chamber is used to maintain the environmental temperature of the experiments. The experimental data such as current and voltage are measured by the Chroma 72001 and recorded by the host computer.

Battery Parameter Estimation
The pulse current charging profile is executed in 20 °C , 25 °C , 30 °C , 35 °C and 40 °C respectively to estimate the battery charging parameter under different temperature. The current pattern and corresponding battery terminal voltage under the temperature of 25 °C is shown in Figure 3a. Figure 3b shows the comparison of the voltage response between the experiment and simulation. The absolute error is also shown in Figure 3b. It can be seen that the simulated results coincide well with the experimental results. Figure 4 is the estimated battery parameter under the temperature of 25 °C .

Problem Formulation
The aim of this paper is to get a balanced charging performance between charging time, charged capacity and energy loss. These charging performances can be divided into two types: the smaller the better for charging time and energy loss, the larger the better for charged capacity. In order to unify

Problem Formulation
The aim of this paper is to get a balanced charging performance between charging time, charged capacity and energy loss. These charging performances can be divided into two types: the smaller the better for charging time and energy loss, the larger the better for charged capacity. In order to unify the two different types of problem, the concept of uncharged capacity (UCP) is adopted, which is defined as follows: The charged capacity is calculated as: where I L is the charging current of the battery. So, the multi-objective optimization problem can be formulated as:

Constraint Conditions
It should be noted that the maximum current for the battery used in this paper is 2 C (C rate, defined as a charge or discharge rate equal to the capacity of a battery divided by 1 h), so the maximum charging current is set as 2 C. Since the voltage limit of each stage is the same, if the charging current of the present stage is greater than that of the previous stage, the rise of battery polarization voltage and ohmic voltage will immediately make the battery terminal voltage higher than the cut-off voltage, so the charging current of the present stage should not exceed the previous stage. The current interval of 0.01 C is used in this paper and the minimum charging current is set as 0.1 C. Considering the overnight charging condition as the longest charging time, the charging time is constrained to be less than ten hours. The uncharged capacity is constrained to be less than 20% of the rated battery capacity considering the driving range and more than 10% of the rated battery capacity considering the battery cycling life [27,28].
According to the above statement, the constraint conditions can be expressed as:

Multi-Object Optimization
As for the problem of multi-objective optimization (MOO), there is not a global best solution that would be optimal for all the objectives simultaneously, so a compromised solution must be found. The MOPSO algorithm is proposed by Coello et al. [29] in 2004 to solve the multi-objective optimization problem and has been widely researched and applied [30][31][32]. Here, the MOPSO is employed because of its high searching ability and low time complexity [29]. In the MOPSO algorithm, every particle is a potential global optimum which has a fitness value evaluated by the object function.
According to the fitness value, the local best value P best and the global best value G best is used to guide the direction of the particles' movement. However, in the problem of MOO, the optimal value cannot be decided directly as different objects conflict with each other. To solve this problem, the Pareto Envelope and grid-making technique is introduced to the MOPSO algorithm and an external repository is used to keep a historical record of the nondominated vectors found along the search process. The detailed optimization procedure is summarized by the following steps.
Step 1 Particle initialization. An m-size population with random positions and zero velocities is generated. The position of the i-th particle is represented as . . x n i and the velocity of the i-th particle is represented as where n is the dimension of the search space.
Step 2 Evaluate the particles according to the simulation result of the MCC charging strategy.
Step 3 Store the nondominated particles in the repository (REP) and update the REP.
Step 4 Update the speed and velocity of each particle according to the following equations: where v n i (k + 1) and v n i (k) denote the current and previous velocities in the n-th dimension, respectively. x n i (k + 1) is the current position of the i-th particle in dimension n while x n i (k) is its previous position; ω is the inertia weight. Λ 1 is the cognitive learning factor, Λ 2 is the social learning factor. R 1 and R 2 are two random numbers in the range [0, 1]. P n best,i is the local best value and REP(h) is the global best value chosen from the repository.
To determine REP(h), the explored objectives are meshed with hypercubes, and a fitness value is assigned to the hypercubes containing more than one particle with the following function: where, n i is the number of particles located in the i-th hypercube.
Then REP(h) is selected by applying a roulette-wheel selection.
Step 5 Convergence determination. Compare the average fitness value of the current particle swarm with the previous particle swarm; if the difference is less than the threshold value α, the search procedure is terminated.

Step 6
Repeat the steps from Step 2 to Step 5 until the best point is found or a fixed number of iterations has been reached.

Parameter Setting
To make the searching procedure work effectively, proper parameters of the algorithm should be chosen. As for the number of the particles, the more particles used in the search procedure, the less steps of iterations are required. However, too many particles would increase the burden of computation. In this research, the number of particles is chosen as 100 to ensure good convergence and accuracy. The maximum iteration number is set as 200 and the inertia weight ω is set as 0.7; the constants Λ 1 and Λ 2 are set as 1.6 and 1.5 respectively and the grid number is chosen as 10.

Multi-Objective Decision Making
To get a balanced charging strategy, a compromised solution should be chosen from the resulted Pareto set. To address this issue, the weighted metric method is adopted in this paper. Since the three objects are in different units and scales, they are first normalized to be between 0 and 1 with the following equation: where f k (X) is the value of the k-th object, f k (X) is the normalized vale and f kmin (X) and f kmax (X) is the minimum and maximum value respectively. Then the weighted objective function can be written as: J ω = ω 1 T CT (I) + ω 2 Q EL (I) + ω 3 C UCP (I) (23) where ω 1 is the weight factor of charging time, ω 2 is the weight factor of uncharged capacity and ω 3 is the weight factor of energy loss.

Impact of Stage Numbers
Optimizations of three, five and seven charging stages are conducted to analyze the influence of charging stage numbers. The Pareto sets of the optimization result are shown in Figure 5. From the distribution of the Pareto set in Figure 6a, it is clear that there is not an optimal solution for all the three objectives. From Figure 5b,c, it can be observed that charging time is in conflict with uncharged capacity and energy loss. In Figure 5d, the Pareto set is a single point, which means the uncharged capacity and energy loss are both optimal in this point. Optimizations of three, five and seven charging stages are conducted to analyze the influence of charging stage numbers. The Pareto sets of the optimization result are shown in Figure 5. From the distribution of the Pareto set in Figure 6a, it is clear that there is not an optimal solution for all the three objectives. From Figure 5b,c, it can be observed that charging time is in conflict with uncharged capacity and energy loss. In Figure 5d, the Pareto set is a single point, which means the uncharged capacity and energy loss are both optimal in this point.  It can be observed from Figure 5b that as the stage number increases, the Pareto set moves towards the origin-both the optimized charging time and uncharged capacity decrease. The decrease is relatively large when the stage number increases from three to five while the difference between stage number of five and seven is not apparent. From Figure 5c, it can be seen that the optimization results are almost the same while stage number of three and five have a slightly better performance than stage number of seven. This is because the increase in the charging stage number augments the average charging current. This relationship can be further confirmed in Figure 5d. In Figure 5d, the optimization result is a single point and as the stage number increases, the energy loss decreases while the uncharged capacity stays the same. According to the former analysis, the charging stage number is chosen as five in this paper. It can be observed from Figure 5b that as the stage number increases, the Pareto set moves towards the origin-both the optimized charging time and uncharged capacity decrease. The decrease is relatively large when the stage number increases from three to five while the difference between stage number of five and seven is not apparent. From Figure 5c, it can be seen that the optimization results are almost the same while stage number of three and five have a slightly better performance than stage number of seven. This is because the increase in the charging stage number augments the average charging current. This relationship can be further confirmed in Figure 5d. In Figure 5d, the optimization result is a single point and as the stage number increases, the energy loss decreases while the uncharged capacity stays the same. According to the former analysis, the charging stage number is chosen as five in this paper.

Impact of Cut-Off Voltage
The effect of cut-off voltage on the charging performance is studied. To charge the battery to the target capacity, the cut-off voltage must be higher than the corresponding OCV, so the low-level cut-off voltage is set as 3.4 V. The high-level cut-off voltage is set as the manufacture's recommended voltage, which is 3.65 V. The medium-level cut-off voltage is set as 3.5 V. The resulted Pareto set is shown in Figure 6.
It can be observed in Figure 6a,c that a higher cut-off voltage leads to a smaller charging time and a higher energy loss. However, the difference of energy loss for the same charging time is

Impact of Cut-Off Voltage
The effect of cut-off voltage on the charging performance is studied. To charge the battery to the target capacity, the cut-off voltage must be higher than the corresponding OCV, so the low-level cut-off voltage is set as 3.4 V. The high-level cut-off voltage is set as the manufacture's recommended voltage, which is 3.65 V. The medium-level cut-off voltage is set as 3.5 V. The resulted Pareto set is shown in Figure 6.
It can be observed in Figure 6a,c that a higher cut-off voltage leads to a smaller charging time and a higher energy loss. However, the difference of energy loss for the same charging time is relatively small. From the optimization result in Figure 6b, it can be concluded that the cut-off voltage affects the charging time largely and a higher cut-off voltage also leads to a smaller uncharged capacity. In Figure 6d, the Pareto set of 3.65 V and 3.5 V coincide with each other and is better than the result of 3.4 V. In this paper, the manufacture's recommended cut-off voltage is used.

Impact of Weight Factors
According to the former analysis, a five stage VMCC charging strategy with the cut-off voltage of 3.65 V is chosen as the final charging strategy. The distribution of the resulted Pareto set is shown in Figure 7a. It can be found that the shortest charging time is 0.43 h (marked with "*"), where the corresponding charging current is I = [1.91 C, 1.73 C, 1.61 C, 1.55 C, 1. The resulted Pareto set provides all the possible optimal solutions for the five stage VMCC charging strategy. In actual usages, a compromised solution should be chosen from the Pareto set according to the specific charging demand. In this paper, the weight metric factor method is adopted to decide the optimal solution. To study the influence of weight factor, the charging performance under three different is compared, where 1 (the weight factor of charging time) is set as 0.4, 0.6 and 0.8, 2 (the weight factor of uncharged capacity) and 3 (the weight factor of energy loss) is set to be equal.
The charging current, SOC and energy loss performance under different weight factors are shown in Figure 7b-d. Figure 7b shows that high values of 1 result in high values of charging current in the first CC stage. In other words, as more emphasis is put on charging time, a larger charging current in the first CC stage is required to shorten the charging time. Figure 7c shows the performance of SOC. It is observed that by increasing 2 , more emphasis is put on uncharged capacity and a smaller uncharged capacity is acquired. It should also be noted that, as 2 increases further, the uncharged capacity does not decrease accordingly. This is caused by the constraint condition of the minimum uncharged capacity. Figure 7d is the performance of energy loss. As expected, high values of 3 result in lower energy loss as more emphasis is put on energy loss, which The resulted Pareto set provides all the possible optimal solutions for the five stage VMCC charging strategy. In actual usages, a compromised solution should be chosen from the Pareto set according to the specific charging demand. In this paper, the weight metric factor method is adopted to decide the optimal solution. To study the influence of weight factor, the charging performance under three different ω is compared, where ω 1 (the weight factor of charging time) is set as 0.4, 0.6 and 0.8, ω 2 (the weight factor of uncharged capacity) and ω 3 (the weight factor of energy loss) is set to be equal.
The charging current, SOC and energy loss performance under different weight factors are shown in Figure 7b-d. Figure 7b shows that high values of ω 1 result in high values of charging current in the first CC stage. In other words, as more emphasis is put on charging time, a larger charging current in the first CC stage is required to shorten the charging time. Figure 7c shows the performance of SOC. It is observed that by increasing ω 2 , more emphasis is put on uncharged capacity and a smaller uncharged capacity is acquired. It should also be noted that, as ω 2 increases further, the uncharged capacity does not decrease accordingly. This is caused by the constraint condition of the minimum uncharged capacity. Figure 7d is the performance of energy loss. As expected, high values of ω 3 result in lower energy loss as more emphasis is put on energy loss, which means a higher efficiency.

Experimental Comparison with the Traditional Charging Strategies
Experimental comparisons of the proposed charging strategy with the traditional charging strategies are carried out. Without loss of generality, the traditional EV charging strategy is classified as normal charging and fast charging. For normal charging, the charging strategy is defined as CC-CV charging with the CC stage current of 0.3 C and for fast charging, the charging strategy is defined as CC-CV charging with the CC stage current of 2.0 C. The concept of user satisfaction is adopted here to evaluate the charging performance of the optimized charging strategy, which is defined as: For the calculation of the weight factor ω, ω 1 is set to change from 0.1 to 0.9 with the interval of 0.1 while ω 2 is equal to ω 3 . According to the definition of η, it is easy to see that the larger the value of η, the better the charging performance satisfies the users' charging demand. Figure 8 shows the experimental comparisons of the proposed charging strategy with the traditional normal and fast charging strategies. It can be observed that, for the traditional normal charging strategy, the charging performance is good when the charging demand for charging time is low; however, as the demand for the charging time increases, the charging performance descends linearly. Conversely, for the traditional fast charging strategy, the charging performance is good when the charging demand for charging time is high, and as the demand for the charging time decreases, the charging performance decreases linearly. It can be concluded that the traditional normal and fast charging strategies only satisfy a small range of the charging demand. On the other hand, it can be seen that the proposed charging strategy always performs well no matter how the charging demand changes, which means the proposed charging strategy satisfies the whole range of the charging demand well. CC-CV charging with the CC stage current of 0.3 C and for fast charging, the charging strategy is defined as CC-CV charging with the CC stage current of 2.0 C. The concept of user satisfaction is adopted here to evaluate the charging performance of the optimized charging strategy, which is defined as: For the calculation of the weight factor ω, 1 is set to change from 0.1 to 0.9 with the interval of 0.1 while 2 is equal to 3 . According to the definition of , it is easy to see that the larger the value of , the better the charging performance satisfies the users' charging demand. Figure 8 shows the experimental comparisons of the proposed charging strategy with the traditional normal and fast charging strategies. It can be observed that, for the traditional normal charging strategy, the charging performance is good when the charging demand for charging time is low; however, as the demand for the charging time increases, the charging performance descends linearly. Conversely, for the traditional fast charging strategy, the charging performance is good when the charging demand for charging time is high, and as the demand for the charging time decreases, the charging performance decreases linearly. It can be concluded that the traditional normal and fast charging strategies only satisfy a small range of the charging demand. On the other hand, it can be seen that the proposed charging strategy always performs well no matter how the charging demand changes, which means the proposed charging strategy satisfies the whole range of the charging demand well. Figure 9 shows the relative increase in user satisfaction (represented with ) of the proposed charging strategy compared with the traditional normal and fast charging strategies. It can be seen that the increase can reach more than 80% when the proposed charging strategy is compared with the normal and fast charging strategies independently. In actual usage, when a short charging time is needed, the EV users may choose fast charging and otherwise the normal charging is used. In Figure 9, the "Switch point" is defined as the intersection of the "Fast charging" line and the "Normal charging" line. For ideal conditions, normal charging is used when 1 is smaller than the switch point and fast charging is used otherwise. In this case, the relative increase in charging performance of the proposed charging strategy can still reach almost 40%, which shows the advantage of the proposed charging strategy.   Figure 9 shows the relative increase in user satisfaction (represented with ε) of the proposed charging strategy compared with the traditional normal and fast charging strategies. It can be seen that the increase can reach more than 80% when the proposed charging strategy is compared with the normal and fast charging strategies independently. In actual usage, when a short charging time is needed, the EV users may choose fast charging and otherwise the normal charging is used. In Figure 9, the "Switch point" is defined as the intersection of the "Fast charging" line and the "Normal charging" line. For ideal conditions, normal charging is used when ω 1 is smaller than the switch point and fast charging is used otherwise. In this case, the relative increase in charging performance of the proposed charging strategy can still reach almost 40%, which shows the advantage of the proposed charging strategy.

Conclusions
An optimal battery charging strategy is proposed to satisfy the different charging demands of EV users for charging time, charge capacity and energy loss based on multi-objective optimization technology. The battery equivalent circuit model is adopted, and the least square method is used to estimate the battery parameters through a pulse current charging experiment. A VMCC simulation model is proposed and the MOPSO algorithm is used to solve the multi-objective optimization problem. The impacts of the charging stage number, charging cut-off voltage and weight factors of different charging goals are investigated in detail. We found that as the charging stage number increases, a better optimization result is achieved; however, when the stage number is greater than five, the increase in charging performance is not evident. Thus, a charging stage number of five is chosen. The analysis of charging cut-off voltage shows that the use of higher charging cut-off voltage benefits the charging performance and the manufacture's recommended charging cut-off voltage is used. The impacts of different weight factor combinations are also studied and a compromised optimal charging strategy is chosen based on the analysis result. Finally, an experimental comparison of the proposed charging strategy with the traditional EV charging strategies are performed. The comparison results show that the traditional normal and fast charging strategies can only satisfy a small range of charging demand while the proposed charging strategy satisfies the whole range of charging demand well. The relative increase in charging performance of the proposed charging strategy can reach more than 80% when compared to normal and fast charging.

Conclusions
An optimal battery charging strategy is proposed to satisfy the different charging demands of EV users for charging time, charge capacity and energy loss based on multi-objective optimization technology. The battery equivalent circuit model is adopted, and the least square method is used to estimate the battery parameters through a pulse current charging experiment. A VMCC simulation model is proposed and the MOPSO algorithm is used to solve the multi-objective optimization problem. The impacts of the charging stage number, charging cut-off voltage and weight factors of different charging goals are investigated in detail. We found that as the charging stage number increases, a better optimization result is achieved; however, when the stage number is greater than five, the increase in charging performance is not evident. Thus, a charging stage number of five is chosen. The analysis of charging cut-off voltage shows that the use of higher charging cut-off voltage benefits the charging performance and the manufacture's recommended charging cut-off voltage is used. The impacts of different weight factor combinations are also studied and a compromised optimal charging strategy is chosen based on the analysis result. Finally, an experimental comparison of the proposed charging strategy with the traditional EV charging strategies are performed. The comparison results show that the traditional normal and fast charging strategies can only satisfy a small range of charging demand while the proposed charging strategy satisfies the whole range of charging demand well. The relative increase in charging performance of the proposed charging strategy can reach more than 80% when compared to normal and fast charging.