Multi-Objective Optimal Charging Method for Lithium-Ion Batteries

In order to optimize the charging of lithium-ion batteries, a multi-stage charging method that considers the charging time and energy loss as optimization targets has been proposed in this paper. First, a dynamic model based on a first-order circuit has been established, and the model parameters have been identified. Second, on the basis of the established model, we treat the objective function of the optimization problem as a weighted sum of charging time and energy loss. Finally, a dynamic programming algorithm (DP) has been used to calculate the charging current of the objective function. Simulation and experimental results show that the proposed charging method could effectively reduce the charging time and decrease the energy loss, compared with the constant-current constant-voltage charging method, under the premise of exerting little influence on the attenuation of battery capacity.


Introduction
As the key technology of battery management systems, the charging strategy influences the battery life and safe operation [1][2][3].Their long charging time is one of the main obstacles for the popularization of lithium batteries [4,5].Moreover, the battery loses energy in the form of heat dissipation during the charging process, and the literature [6] indicates that within a certain temperature range, the higher the battery temperature, the more serious the battery life decay.Therefore, effective reduction of the charging time and battery energy losses during the charging period has become one of the important directions for the research on lithium-ion battery charging strategies [7][8][9].
The constant-current constant-voltage (CCCV) charging method is one of the most widely used charging methods; further, it is simple and easy to control.However, different charging rates and cut-off voltages result in battery aging [10].Accordingly, a multi-stage constant current charging method for multi-objective optimization has been studied [11,12].The target optimization mainly aims to reduce the charging time, improve the charging efficiency, and prolong the battery life; however, it is necessary to calibrate numerous experimental data and therefore, such a method is difficult to use.Reference [13] proposed a fast charging method for a lithium battery close to complete discharge.However, such a method can result in an overheated battery at the expense of battery life and requires a better solid electrolyte interface or unique electrode structures to mitigate any potential lithium plating [14][15][16].References [17][18][19][20] involved the study of a pulse charging method, wherein the battery was stationary or discharged for a short period of time in order to slow down the polarization phenomenon during the charging process.However, the proposed method has higher requirements for the design of the charging power supply.
In order to improve the charging efficiency and prolong the battery life, a general charging voltage protocol based on battery state-of-health was proposed in [21], which aims to study the optimization of charging current.Moreover, the battery charging method has also been studied in optimization problems related to charging time, battery energy loss, and temperature rise in some studies.The effect of cut-off voltage, ambient temperature, and battery aging on the charging optimization results was also discussed in [22], whereas the influence of charge current on the Ohmic resistance, polarization resistor, and polarization capacitance of the battery has not been fully considered yet.Although battery charging current has been optimized from the perspective of charging time and battery energy losses in [23,24], the battery model that has been adopted is the voltage-resistor model, which ignores the influence of polarization phenomena on charging.
In summary, previous battery charge management studies mainly concentrated on the optimization of the charging current.In spite of the fact that some of them considered charging time and energy loss as characteristics of the charging method, the effect of the change of battery model parameters on the optimization of charging current was not considered during the optimization process.Based on the dynamic model of the equivalent circuit, a lithium cobalt-acid battery has been considered as the research objective in this study.Moreover, considering the influence of the model parameters on the state of charge (SOC) and charging current and considering the charging time and battery energy loss as the optimal targets, over-voltage is employed as a state variable, and a dynamic programming algorithm is also adopted to obtain the optimized charging current.Further, the influence of variable weighting factors on the optimization result is also discussed in this paper.

Charging Principle and Model Selection
In this study, a lithium cobalt acid battery is selected as the research object.Its details are listed in Table 1.A lithium-ion battery consists of three parts: the anode lithium compound, its intermediate electrolyte membrane, and cathode carbon.When it is charged, lithium ions are generated at the anode and embedded through the electrolyte membrane into the cathode material.As more lithium ions are embedded into the cathode material, the SOC level of the battery becomes higher.The chemical reaction during the lithium-ion battery charging process is given as follows: Moreover, during the battery charging process, since the lithium ion concentration fails to reach a balance within a short period of time, a voltage difference, i.e., the over-voltage, arises.In order to describe the charging process accurately, many battery models have been established in previous studies in order to reflect the mechanism changes in the charging process.Common lithium battery models include an equivalent circuit model and an electrochemical model; further, the equivalent circuit model is simple and can accurately reflect the running state of the battery [25,26].Therefore, this study refers to literature [27,28], and the first-order equivalent circuit model has been chosen as the research object, as shown in Figure 1.
Energies 2017, 10, 271 3 of 18 this study refers to literature [27,28], and the first-order equivalent circuit model has been chosen as the research object, as shown in Figure 1.As shown in Figure 1, the equivalent circuit model of the battery consists of the following parts: open-circuit voltage which is a controlled voltage source UOC, Ohmic resistance R0, polarization resistance R1, and polarization capacitance C1.According to [29], all the parameters in the model are related to the SOC and charge-discharge ratio of the battery under the condition of fixed ambient temperature.

Parameter Identification and Verification of Battery Model
In order to identify the parameters of the model for different SOC and charging ratios, the battery have been tested with several sets of specific charging ratios by adopting the flowchart shown in Figure 2.

Still 5min
The voltage is greater than the charge cut-off voltage

Constant current charging
The charge test is divided into N steps per change As evident from the test flow in Figure 2, in each test, the battery has been subjected to a charge process with a pulse of 10 s for every SOC of 5%.After the charging test is completed, the relationship between the battery terminal voltage and charging current corresponding to each pulse charging stage is fitted by using the least squares method mentioned in [30].The corresponding parameter identification results of the battery model are shown in Figure 3.As shown in Figure 1, the equivalent circuit model of the battery consists of the following parts: open-circuit voltage which is a controlled voltage source U OC , Ohmic resistance R 0 , polarization resistance R 1 , and polarization capacitance C 1 .According to [29], all the parameters in the model are related to the SOC and charge-discharge ratio of the battery under the condition of fixed ambient temperature.

Parameter Identification and Verification of Battery Model
In order to identify the parameters of the model for different SOC and charging ratios, the battery have been tested with several sets of specific charging ratios by adopting the flowchart shown in Figure 2. this study refers to literature [27,28], and the first-order equivalent circuit model has been chosen as the research object, as shown in Figure 1.As shown in Figure 1, the equivalent circuit model of the battery consists of the following parts: open-circuit voltage which is a controlled voltage source UOC, Ohmic resistance R0, polarization resistance R1, and polarization capacitance C1.According to [29], all the parameters in the model are related to the SOC and charge-discharge ratio of the battery under the condition of fixed ambient temperature.

Parameter Identification and Verification of Battery Model
In order to identify the parameters of the model for different SOC and charging ratios, the battery have been tested with several sets of specific charging ratios by adopting the flowchart shown in Figure 2.

Still 5min
The voltage is greater than the charge cut-off voltage

Constant current charging
The charge test is divided into N steps per change As evident from the test flow in Figure 2, in each test, the battery has been subjected to a charge process with a pulse of 10 s for every SOC of 5%.After the charging test is completed, the relationship between the battery terminal voltage and charging current corresponding to each pulse charging stage is fitted by using the least squares method mentioned in [30].The corresponding parameter identification results of the battery model are shown in Figure 3.As evident from the test flow in Figure 2, in each test, the battery has been subjected to a charge process with a pulse of 10 s for every SOC of 5%.After the charging test is completed, the relationship between the battery terminal voltage and charging current corresponding to each pulse charging stage is fitted by using the least squares method mentioned in [30].The corresponding parameter identification results of the battery model are shown in Figure 3.In order to verify the accuracy of the model parameters, this study selects charging current with different specific ratios to simulate the intermittent charging of the battery based on the targeted dynamic model.Firstly, the equivalent circuit model of a battery was established under Matlab/Simulink, and the charging simulation was carried out by using some specific current.Secondly the battery is tested with the same current as the simulation.Finally we compare the battery terminal voltage between simulation and test.Considering the simulation and actual test with charging ratio of 0.4C as an example, the current and voltage curves are shown in Figure 4. Figure 5 illustrates a comparison between the errors of the calculated model value and the actual test battery result with various charging ratios.The error is calculated as: where U t,r is the actual value of the battery terminal voltage at time t and U t,s is the simulation value.
Energies 2017, 10, 271 5 of 18 In order to verify the accuracy of the model parameters, this study selects charging current with different specific ratios to simulate the intermittent charging of the battery based on the targeted dynamic model.Firstly, the equivalent circuit model of a battery was established under Matlab/Simulink, and the charging simulation was carried out by using some specific current.Secondly the battery is tested with the same current as the simulation.Finally we compare the battery terminal voltage between simulation and test.Considering the simulation and actual test with charging ratio of 0.4C as an example, the current and voltage curves are shown in Figure 4. Figure 5 illustrates a comparison between the errors of the calculated model value and the actual test battery result with various charging ratios.The error is calculated as: where Ut,r is the actual value of the battery terminal voltage at time t and Ut,s is the simulation value.It is evident from Figure 4 that the battery terminal voltage of the targeted model is quite close to the actual battery test results when the charging ratio is 0.4C.It can be observed from Figure 5 that the average error of the battery terminal voltage calculated using the targeted model relative to the actual voltage is less than 0.1%, and the maximum error is less than or equal to 6.12% with different charging ratios.Therefore, it can be demonstrated that the targeted battery dynamic model can accurately represent the parameter change of the battery charging process.

Multi-Objective Optimal Charging Method
During the battery charging, the charging time and battery energy loss are important indexes to evaluate the charging method.When analyzing the charging process using the first-order equivalent circuit battery model, the battery's energy loss originates from the Ohmic resistance and polarized internal resistance.However, the charging time is usually inversely proportional to the charging ratio [31].The key problem of the charging method is to balance the charging time and battery energy loss.

Determination of the Boundary Condition of Charging Current
Zhang et al. [32] observed that the maximum charge current of a lithium-ion battery is usually related to the current SOC.Accordingly, this study performs many charging tests, and subsequently, the test results are fitted; finally, the maximum charge current corresponding to the SOC is obtained.By using the test flow shown in Figure 6, the maximum charge capacity corresponding to the charge current ratio in the range of 0.1C-1C can be obtained.The ratio of maximum charge capacity to the rated capacity, i.e.,  , or in other words, the percentage of maximum charge capacity, serves as an independent variable, and the charging current acts as a dependent variable to obtain a fitting.Thus, the fitted curve shown in Figure 7 is obtained, and the results are given as follows: According to [33,34], when the charging current is greater than the security threshold corresponding to the current SOC, the battery aging accelerates or heat might arise out of the control.In order to avoid damage to the battery, the maximum charge current corresponding to the aforementioned charge capacity  is considered as the boundary condition of the charging current selection for different SOCs, i.e.It is evident from Figure 4 that the battery terminal voltage of the targeted model is quite close to the actual battery test results when the charging ratio is 0.4C.It can be observed from Figure 5 that the average error of the battery terminal voltage calculated using the targeted model relative to the actual voltage is less than 0.1%, and the maximum error is less than or equal to 6.12% with different charging ratios.Therefore, it can be demonstrated that the targeted battery dynamic model can accurately represent the parameter change of the battery charging process.

Multi-Objective Optimal Charging Method
During the battery charging, the charging time and battery energy loss are important indexes to evaluate the charging method.When analyzing the charging process using the first-order equivalent circuit battery model, the battery's energy loss originates from the Ohmic resistance and polarized internal resistance.However, the charging time is usually inversely proportional to the charging ratio [31].The key problem of the charging method is to balance the charging time and battery energy loss.

Determination of the Boundary Condition of Charging Current
Zhang et al. [32] observed that the maximum charge current of a lithium-ion battery is usually related to the current SOC.Accordingly, this study performs many charging tests, and subsequently, the test results are fitted; finally, the maximum charge current corresponding to the SOC is obtained.By using the test flow shown in Figure 6, the maximum charge capacity corresponding to the charge current ratio in the range of 0.1C-1C can be obtained.The ratio of maximum charge capacity to the rated capacity, i.e., ξ, or in other words, the percentage of maximum charge capacity, serves as an independent variable, and the charging current acts as a dependent variable to obtain a fitting.Thus, the fitted curve shown in Figure 7 is obtained, and the results are given as follows: According to [33,34], when the charging current is greater than the security threshold corresponding to the current SOC, the battery aging accelerates or heat might arise out of the control.In order to avoid damage to the battery, the maximum charge current corresponding to the aforementioned charge capacity ξ is considered as the boundary condition of the charging current selection for different SOCs, i.e., I max (SOC) = The maximum charging current(A) The percentage of maximum charge capacity(%) Figure 7. Relationship between the maximum charge capacity and its corresponding maximum charge current.

Establishment of Optimization Objective Function
This study considers the battery energy loss and charging time as the optimization objective function.During the process of formulating this function, the charging process is first divided into N constant current charging stages, and the parameters of the model are constant in each charging stage.During the entire charging process, the battery's energy loss can be expressed as:  The maximum charging current(A) The percentage of maximum charge capacity(%) Figure 7. Relationship between the maximum charge capacity and its corresponding maximum charge current.

Establishment of Optimization Objective Function
This study considers the battery energy loss and charging time as the optimization objective function.During the process of formulating this function, the charging process is first divided into N constant current charging stages, and the parameters of the model are constant in each charging

Establishment of Optimization Objective Function
This study considers the battery energy loss and charging time as the optimization objective function.During the process of formulating this function, the charging process is first divided into N constant current charging stages, and the parameters of the model are constant in each charging stage.During the entire charging process, the battery's energy loss can be expressed as: where I j represents the charge current in the jth charging phase, and R 0,j (I j ) represents the corresponding Ohmic resistance of the jth phase, which is affected by the current SOC and charging current.Furthermore, W 1,j is the energy loss of the polarization resistance in the jth phase, which is related to the charging current I j , polarization capacitance C 1,j (I j ), and polarization resistance R 1,j (I j ) of the current stage; further, ∆ j indicates the time duration of the jth charging phase.
The SOC changes in each charging phase are expressed as follows: where Cap represents the battery rated capacity with the unit of ampere hour and η is the Coulomb efficiency.In this paper, 1 has been taken [30], and when the charging process is divided by 1% of ∆SOC, the charging time of each charging phase is: When the sampling time of the charging system is ∆t, the number of samples in each charging phase is: The energy loss of the polarization resistance at each charging stage separated by the sampling time ∆t is given as: where I 1,j,k is the current of the polarization resistance R 1,j (I j ) at the K-sampling moment in the jth charging phase.According to the full response equation of the first-order circuit, the following can be obtained: where U 1,j represents the polarization voltage at the starting point of the jth charging phase.By substituting (10) into (9), the following can be obtained: Substituting ( 11) into (5), the total energy consumption of the battery during multi-stage constant current charging can be obtained as: During the entire charging process, the charging time of the battery is: By setting the sampling time of the system as ∆t = 1s, the objective function of the optimized charging method presented in this paper can be expressed as: The following can be obtained after further simplification: where α is the obtained weighting factor, which could be adjusted between 0 and 1 to achieve the optimization of different focus on charging time and energy loss.M serves as the compensating coefficient, which allows the charging time and energy loss to share the same order of magnitude.The control variable is I j , i.e., the charging current in the charging stage, and the state variable is the voltage U 1,j at the ends of the polarization resistance R 1,j (I j ) at the beginning of each charging stage.
Owing to the parallel connection of polarization resistance and polarization capacitor in the equivalent circuit model, the voltage drop in the polarization resistance cannot be changed suddenly, and hence, the initial over-voltage of each charging stage is equal to the termination over-voltage of the corresponding previous charging stage.Accordingly, the initial over-voltage U 1,j of each charging stage is chosen as the state variable, and the state transfer equation is: The boundary conditions of the charge currents indicate that, when the SOC is large, the charging current is small, which leads to the slow charging of the last few stages.At the end of the charge, it is believed that the current in the polarization resistance reaches the steady state value I N .Thus, according to (16), the initial over-voltage during the final charging phase can be obtained, i.e., the state variables U 1,N = R 1,N (I N ) × I N At the initial stage of charging, the state variable of such a phase is U 1,1 = 0, because the voltage of polarization resistance is zero.

Solution of Charging Current using a Dynamic Programming Algorithm
Dynamic programming theory, proposed by the U.S. mathematician R. Bellman, is an overall optimization method aiming to transform complex decision problems into a series of sub-stage decision problems, and it has been widely used in offline global optimization of nonlinear systems.In this study, the charging method is based on the battery equivalent dynamic model, and the parameters of the model are nonlinear with the changes of SOC and charging current.Therefore, the dynamic programming algorithm is suitable for solving the current of the proposed charging method.The dynamic programming model for obtaining the charging current is shown in Table 2.
Table 2. Mathematical model of dynamic programming algorithm.

Phase variables
The charge stage is divided according to ∆SOC; each charging stage is denoted by j = 1, 2, ..., N.
State variables U 1,j is the over-voltage at the initial moment of the jth charging phase.
Decision variable I j is the charging current during the jth charging phase.
Enable decision collection D j (I j ) = {0 ≤ I j ≤ I max (∆SOC) × j} State transfer equation The inverse solution method is used to solve the problem.The N + 1 charging stage is added as the boundary condition and the objective value function at this stage is zero, i.e., f N+1 = 0. First, the extreme value of f N is obtained; further, the expression of optimal charging current in the N stage is obtained, which is I N,opt = M N (U 1,N ).Second, if U 1,N and I N,opt are used for obtaining f N , the f N expression about U 1,N is obtained, which is P N (U 1,N ), According to the state transfer equation, P N (U 1,N ) can be further expressed as P N (φ(U 1,N-1 ,I N-1 )).After storing I N,opt and f N , the same method is used to calculate f N-1 , and the expressions of I N-1,opt and f N-1 are represented as M N-1 (U 1,N-1,opt ) and P N-1 (φ(U 1,N-2 ,I N-2 )), respectively.After storing I N-1,opt and f N-1 , f N-2 is calculated.Thus, this process is performed step by step, obtaining each stage of the charging current I j,opt in the expression for U 1,j , i.e., Mj(U1,j).According to the initial condition U 1,1 = 0 and the state transition equation U 1,j = φ(U 1,j-1 ,I j-1 ), the optimal charging current I j at each stage is calculated.The solution flow is shown in Figure 8.

Optimization Results for Different Weighting Factors
When the weighting factor changes from 0 to 1 in the optimization target, the optimization result of the proposed charging method is shown in Figure 9.The charging range of the battery corresponding to each weighting factor is in the range of 0 to 95% of SOC.

Optimization Results for Different Weighting Factors
When the weighting factor changes from 0 to 1 in the optimization target, the optimization result of the proposed charging method is shown in Figure 9.The charging range of the battery corresponding to each weighting factor is in the range of 0 to 95% of SOC.

Optimization Results for Different Weighting Factors
When the weighting factor changes from 0 to 1 in the optimization target, the optimization result of the proposed charging method is shown in Figure 9.The charging range of the battery corresponding to each weighting factor is in the range of 0 to 95% of SOC.
(a) Figure 9a illustrates an optimized charge current curve for different weighting factors, and Figure 9b illustrates the time and energy loss for different weighting factors corresponding to the proposed charging method.It can be observed from the graph that, when the charging time decreases, the battery energy loss increases, which indicates that the two problems are contradictory with respect to the optimization of battery charging.

Analysis of the Weighting Factors
In order to verify the effectiveness of the proposed charging method and to provide reference for selecting the optimum weighting factor, this study has compared the proposed method and the CCCV charging method.The weighting factor the proposed method is in the range of 0.1-0.9,whereas the charging current ratio of the CCCV charging method is 0.5C.The charging range is in the SOC range of 0 to 95%.The charging time and battery energy loss of both these methods are demonstrated in Table 3.The comparison between the two charging methods assumes the following approach: Figure 9a illustrates an optimized charge current curve for different weighting factors, and Figure 9b illustrates the time and energy loss for different weighting factors corresponding to the proposed charging method.It can be observed from the graph that, when the charging time decreases, the battery energy loss increases, which indicates that the two problems are contradictory with respect to the optimization of battery charging.

Analysis of the Weighting Factors
In order to verify the effectiveness of the proposed charging method and to provide reference for selecting the optimum weighting factor, this study has compared the proposed method and the CCCV charging method.The weighting factor of the proposed method is in the range of 0.1-0.9,whereas the charging current ratio of the CCCV charging method is 0.5C.The charging range is in the SOC range of 0 to 95%.The charging time and battery energy loss of both these methods are demonstrated in Table 3.The comparison between the two charging methods assumes the following approach: The aforementioned equations present the percentage time and energy consumption, respectively, of the optimal charging method when compared to the CCCV charging method; t α and E α , respectively, represent the time and energy loss of the proposed charging method corresponding to the weighting factor α; t CCCV and E CCCV , respectively, represent the time and energy loss of the CCCV charging method of 0.5C.The comparison results are shown in Figure 10.When 'Energy save' or 'Time reduce' is positive, it means that the proposed charging method is superior to the constant current constant voltage charging method in reducing the energy loss or reducing the charging time, and the greater the absolute value of 'Energy save' or 'Time reduce', the better the optimization effect and vice versa.
As evident from Figure 10, compared with the CCCV charging method of 0.5C, when the weighting factor is 0.1-0.4,the charging method can reduce the charging time, but the energy loss is increased.However, when the weighting factor is 0.6-0.9, the charging method can reduce the energy loss, but the charging time increases.It is only at the weighting factor of 0.5 that the energy loss decreases by 5.45% and the charging time decreases by 7.66%.This results in a trade-off in optimizing charging time and energy loss in order to change the charging current according to the changes of the internal parameters of the battery along with the changes of the SOC When the internal resistance of the SOC is relatively large, the charging current of the optimized calculation can be relatively reduced in order to reduce the energy consumption of the battery.When the resistance value corresponding to the SOC is small, the optimized calculation of the charging current can be appropriately improved in order to reduce the charging time.
The aforementioned equations present the percentage time and energy consumption, respectively, of the optimal charging method when compared to the CCCV charging method; tα and Eα, respectively, represent the time and energy loss of the proposed charging method corresponding to the weighting factor α; tCCCV and ECCCV, respectively, represent the time and energy loss of the CCCV charging method of 0.5C.The comparison results are shown in Figure 10.When 'Energy save' or 'Time reduce' is positive, it means that the proposed charging method is superior to the constant current constant voltage charging method in reducing the energy loss or reducing the charging time, and the greater the absolute value of 'Energy save' or 'Time reduce', the better the optimization effect and vice versa.As evident from Figure 10, compared with the CCCV charging method of 0.5C, when the weighting factor is 0.1-0.4,the charging method can reduce the charging time, but the energy loss is increased.However, when the weighting factor is 0.6-0.9, the charging method can reduce the energy loss, but the charging time increases.It is only at the weighting factor of 0.5 that the energy loss decreases by 5.45% and the charging time decreases by 7.66%.This results in a trade-off in optimizing charging time and energy loss in order to change the charging current according to the changes of the internal parameters of the battery along with the changes of the SOC When the internal resistance of the SOC is relatively large, the charging current of the optimized calculation can be relatively reduced in order to reduce the energy consumption of the battery.When the resistance value corresponding to the SOC is small, the optimized calculation of the charging current can be appropriately improved in order to reduce the charging time.This study uses the surface temperature of the battery to reflect the energy loss in the charging process indirectly.Figure 11 illustrates the comparison between the proposed charging method with weight factor of 0.5 and CCCV charging method of 0.5C in the charging.The charging range is the SOC of 0 to 95%.This study uses the surface temperature of the battery to reflect the energy loss in the charging process indirectly.Figure 11 illustrates the comparison between the proposed charging method with weight factor of 0.5 and CCCV charging method of 0.5C in the charging.The charging range is the SOC of 0 to 95%.As evident from Figure 11, the maximum temperature rise, average temperature rise, and charging time of the proposed charging method for the weighting factor of 0.5 are reduced, as compared to the CCCV charging method of 0.5°C.The effectiveness of the charging method proposed in this paper and the rationality behind choosing the charging method with the weighting factor of 0.5 to optimize the charging time and energy loss are explained.

Attenuation Comparison of Battery Capacity
In order to evaluate the charging method objectively, this study performs a comparative analysis between the charging method and CCCV charging method from the perspective of battery capacity attenuation.The empirical formula describing the decay of battery life can be obtained from literature [35]: where Caploss is the percentage of battery capacity attenuation, and K and δ represent the parameters related to battery charge-discharge ratio and charge-discharge cut-off voltage and temperature; further, x represents the battery charge-discharge cycles or the throughput capacity of the battery.In literature [35], the relationship between the parameters and the charge-discharge ratio was calibrated at 25 °C with the lithium-ion battery which considers the cathode material as LiCoO2 and the anode material as graphite.Specifically, x represents charge-discharge cycles and the calibration results are shown in Table 4.In order to apply the model of battery decay and calibration parameters of [35] to the research in this study, the aging model has been further deduced, and (19) can be changed to: where Ah represents the throughput capacity of the battery.Further, (20) can be modified to: As evident from Figure 11, the maximum temperature rise, average temperature rise, and charging time of the proposed charging method for the weighting factor of 0.5 are reduced, as compared to the CCCV charging method of 0.5 • C. The effectiveness of the charging method proposed in this paper and the rationality behind choosing the charging method with the weighting factor of 0.5 to optimize the charging time and energy loss are explained.

Attenuation Comparison of Battery Capacity
In order to evaluate the charging method objectively, this study performs a comparative analysis between the charging method and CCCV charging method from the perspective of battery capacity attenuation.The empirical formula describing the decay of battery life can be obtained from literature [35]: where Cap loss is the percentage of battery capacity attenuation, and K and δ represent the parameters related to battery charge-discharge ratio and charge-discharge cut-off voltage and temperature; further, x represents the battery charge-discharge cycles or the throughput capacity of the battery.In literature [35], the relationship between the parameters and the charge-discharge ratio was calibrated at 25 • C with the lithium-ion battery which considers the cathode material as LiCoO 2 and the anode material as graphite.Specifically, x represents charge-discharge cycles and the calibration results are shown in Table 4.In order to apply the model of battery decay and calibration parameters of [35] to the research in this study, the aging model has been further deduced, and (19) can be changed to: where Ah represents the throughput capacity of the battery.Further, (20) can be modified to: By using (20) to achieve the derivation of Ah, we obtain: Substituting ( 21) into (22) in order to achieve difference, we obtain: where Cap loss,p+1 and Cap loss,p , respectively, represent the cumulative loss of the battery at p and p + 1 moments, and ∆Ah p represents the throughput capacity of the battery during the p to p + 1 interval, i.e., According to (20), we obtain: In the battery attenuation model adopted in this study, the charge-discharge current of the battery has a significant influence on the battery attenuation.In order to evaluate the advantages of the charging method proposed in this paper, the 0.5C current ratio will be unified used when discharging during the battery charge and discharge testing process.After multiple charge and discharge cycles, compared with the optimized charging method corresponding to the weighting factor of 0.5 and the CCCV charging method, the decay of battery capacity attenuation is shown in Figure 12.Further, the relative capacity decay rate µ is calculated as: where Cap loss,opt and Cap loss,cccv represent the percentage of the battery capacity attenuation of the proposed charging method and CCCV voltage charging method of 0.5C, respectively, Energies 2017, 10, 271 15 of 18 By using (20) to achieve the derivation of Ah, we obtain: Substituting ( 21) into (22) in order to achieve difference, we obtain: where Caploss,p+1 and Caploss,p, respectively, represent the cumulative loss of the battery at p and p+1 moments, and ΔAhp represents the throughput capacity of the battery during the p to p+1 interval, i.e., According to (20), we obtain: In the battery attenuation model adopted in this study, the charge-discharge current of the battery has a significant influence on the battery attenuation.In order to evaluate the advantages of the charging method proposed in this paper, the 0.5C current ratio will be unified used when discharging during the battery charge and discharge testing process.After multiple charge and discharge cycles, compared with the optimized charging method corresponding to the weighting factor of 0.5 and the CCCV charging method, the decay of battery capacity attenuation is shown in Figure 12.Further, the relative capacity decay rate μ is calculated as:  It can be observed in Figure 12 that, compared with the traditional CCCV charging method, a similarity can be observed in the effect of the optimized charging method on the declining battery life.Reference [36] indicates that lithium cobalt acid batteries are mainly applied in portable mobile electronics.If the electronic products are only used for one year and the daily charge is indispensable, compared to the CCCV charging method, the proposed method reduces the charging time and energy consumption by 7.66% and 5.45%, respectively, by adopting the charging method proposed in this paper, whereas the capacity attenuation increases by merely 1.25%.
Therefore, considering the influence of charging time, battery energy consumption, and capacity attenuation, the charging method proposed in this paper can effectively reduce the charging time and battery energy loss without accelerating the attenuation of the battery capacity.

Conclusions
By establishing a first-order dynamic equivalent circuit model of the battery and realizing parameter identification, a charging method that considers the battery charging time and battery energy loss as the optimized objectives has been proposed in this paper, and the optimal current value has been obtained by using dynamic programming.According to the simulation and experimental results, during charging, reducing the time and reducing the energy loss are contradictory problems; therefore, the charging curve can be obtained with different emphasis on charging time and energy loss by changing the weighting factor.The results of discussion on the weighting factors in this paper are as follows.Compared with the CCCV charging method, when the weighting factor is 0.5, the energy loss and time are reduced by 5.45% and 7.66%, respectively, whereas the effect of the two charging methods on the attenuation of the battery capacity is almost the same.Therefore, the multi-objective optimal charging method can effectively solve the comprehensive optimization problem of the charging time and energy loss.

Figure 1 .
Figure 1.First-order circuit model of battery.

Figure 2 .
Figure 2. Test flow of battery model parameters.

Figure 1 .
Figure 1.First-order circuit model of battery.

Figure 1 .
Figure 1.First-order circuit model of battery.

Figure 2 .
Figure 2. Test flow of battery model parameters.

Figure 2 .
Figure 2. Test flow of battery model parameters.

Figure 3 .
Figure 3. Identification results of battery model parameters.(a) Open circuit voltage identification results; (b) Ohmic resistance identification results; (c) Polarization resistance identification results; (d) Polarization capacitance identification results.

Figure 3 .
Figure 3. Identification results of battery model parameters.(a) Open circuit voltage identification results; (b) Ohmic resistance identification results; (c) Polarization resistance identification results; (d) Polarization capacitance identification results.

Figure 4 .
Figure 4. Simulation and test results with charging ratio of 0. 4C.(a) Charging current curve; (b) Voltage curves for simulation and testing.Figure 4. Simulation and test results with charging ratio of 0. 4C.(a) Charging current curve; (b) Voltage curves for simulation and testing.

Figure 4 .Figure 5 .
Figure 4. Simulation and test results with charging ratio of 0. 4C.(a) Charging current curve; (b) Voltage curves for simulation and testing.Figure 4. Simulation and test results with charging ratio of 0. 4C.(a) Charging current curve; (b) Voltage curves for simulation and testing.

Figure 5 .
Figure 5.Comparison between the errors of the calculated value and the test result of the equivalent circuit model.

Figure 6 .
Figure 6.Test flow for the maximum charge capacity with different charging currents.

Figure 7 .
Figure 7. Relationship between the maximum charge capacity and its corresponding maximum charge current.

Figure 8 .
Figure 8. Flowchart of dynamic programming algorithm adopted to obtain current.

Figure 8 .
Figure 8. Flowchart of dynamic programming algorithm adopted to obtain current.

Figure 8 .
Figure 8. Flowchart of dynamic programming algorithm adopted to obtain current.

Figure 9 .
Figure 9.Comparison between charging current curves for different weights and charging time and energy loss for different weights.(a) Charging current curve; (b) Charging time and energy loss.

Figure 10 .
Figure 10.Comparison between charging method and CCCV charging method for different weight factors.

Figure 10 .
Figure 10.Comparison between charging method and CCCV charging method for different weight factors.

Figure 11 .
Figure 11.Surface temperature rise of the battery during the charging process.

Figure 11 .
Figure 11.Surface temperature rise of the battery during the charging process.

Figure 12 .
Figure 12.Effect of the two charging methods on the attenuation of battery capacity.

Figure 12 .
Figure 12.Effect of the two charging methods on the attenuation of battery capacity.

Table 1 .
Nominal specifications of the battery.

Table 3 .
Time and energy loss of the proposed charging method and CCCV charging method.

Table 3 .
Time and energy loss of the proposed charging method and CCCV charging method.

Table 4 .
Calibration results for parameters k and  .

Table 4 .
Calibration results for parameters k and δ.