Incremental Capacity Curve Peak Points-Based Regression Analysis for the State-of-Health Prediction of a Retired LiNiCoAlO 2 Series / Parallel Conﬁgured Battery Pack

: To recycle retired series / parallel battery packs, it is necessary to know their state-of-health (SOH) correctly. Unfortunately, voltage imbalances between the cells occur repeatedly during discharging / charging. The voltage ranges for the discharge / charge of a retired series / parallel battery pack are reduced owing to the voltage imbalances between the cells. To determine the accurate SOH of a retired series / parallel battery pack, it is necessary to calculate the total discharge capacity using fully discharging / charging tests. However, a fully discharging / charging test is impossible owing to the reduced voltage range. The SOH of a retired series / parallel battery pack with a voltage imbalance should be estimated within the reduced discharging / charging voltage range. This paper presents a regression analysis of the peak point in the incremental capacity (IC) curve from the fresh state to a 100-cycle aging state. Moreover, the SOH of the considered retired series / parallel battery pack was estimated using a regression analysis model. The error in the SOHs of the retired series / parallel battery pack and linear regression analysis model was within 1%, and hence a good accuracy is achieved.


Introduction
In recent years, with the increased use of electric vehicles (EVs), hybrid electric vehicles (HEVs), energy storage system (ESSs), and renewable energy storage system (RESSs), the development of a high-power and high-capacity lithium secondary battery have become active [1]. Because the output power of one cell is limited, multiple cells are connected in series and parallel according to the size of the voltage, current, and capacity of the application. When a battery pack consists of multiple cells, a voltage imbalance occurs between the cells. When such a battery pack is charged and discharged over many cycles, the above-mentioned cell-to-cell voltage imbalance occurs owing to the temperature variation, parasitic resistance, and cell aging. A voltage imbalance between each cell not only degrades the overall capacity utilization efficiency of the battery pack, but also accelerates the aging. Moreover, the full discharge/charge of a battery pack without balancing the voltage difference can lead to an overcharge/overdischarge. A cell-balancing circuit is used to prevent the voltage deviation between cells. However, because of the high manufacturing cost and power loss, it is inefficient.
In the case of electric transportation (ET), when the capacity of the battery pack is reduced to 80% of the fresh capacity, it should be replaced due to safety concerns [2]. Moreover, it is expected that battery's redox reaction, which is an important parameter to determine the capacity degradation of the battery [19,20]. With this method, the voltage plateaus in the discharging/charging curve and can be transformed into clearly identifiable peaks in the IC and differential voltage (DV) curve. The peaks in an IC curve represent the phase equilibria, whereas a DV curve represents the phase transitions [19,21]. Each peak in the curve has unique features, such as intensity and position, and it represents a specific electrochemical process occurring in the cell [19,22]. The extracted peak values and shape and position variation of the peaks are closely related to the battery capacity degradation, and therefore, can be used for monitoring battery aging. Furthermore, the changes of the IC curve peak point according to the battery aging is not affected by the various aging conditions, such as the depth of discharge (DOD), temperature, C-rate, and charging/discharging profile according to the application [22][23][24]. This IC-based method requires a few calculations, and they can be extracted in online applications [2]. This paper presents a regression analysis of the peak point in an IC curve from the fresh state to a 100-cycle aging state. And, the battery used in this research was a 12S28P retired series/parallel battery pack consisting of lithium nickel cobalt aluminum oxide (NCA) cells for a frequency regulation energy storage system (FRESS). First, the aging cycle test of the NCA cells was performed from the fresh state to the 100-cycle state. Second, the IC curve was extracted from the aging cycle test datasets. Third, to solve the measurement noise problems, the low-pass filter (LPF) smoothing method was used to obtain a smooth IC curve. Fourth, the peak point in the IC curve for the SOH estimation of the 12S28P retired NCA series/parallel battery pack was selected. Moreover, for comparison with the peak point in the IC curve method, analysis of the directly measured partial capacity in a short voltage range was performed. Fifth, the change in the IC curve peak point and directly measured partial capacity in the reduced voltage range caused by the aging cycle was modeled using linear regression analysis. Finally, using linear regression analysis, two aging models were experimentally verified on the 12S28P retired NCA series/parallel battery pack datasets.

Experimental Result of the Fresh Series/Parallel Battery Pack
The voltage imbalance of a battery pack mainly occurs between its series connections. As a battery pack repeatedly charges/discharges, this voltage imbalance is caused by the power loss due to the resistance between the series-connected terminals and difference in the temperature. Figure 1 shows the NCA series/parallel battery pack for a FRESS used for an experiment. This battery pack consisted of 336 NCA high-power cells with 12 series connections and 28 parallel connections. One NCA cell had a nominal capacity of 2.0 Ah, a cut-off voltage range of 4.2 V to 2.5 V, and a maximum discharge current rate of 15C-rate (30 A). In comparison, the 12S28P battery pack had a nominal capacity of 56 Ah, a cut-off voltage range of 49.5 V to 36 V, and maximum discharge current rate of 840 A. Figure 2 shows the aging current profile with the maximum current at a 4C-rate and response of the voltage in the experiment. Moreover, Figure 3 shows 12 discharging/charging voltages for each representative cell, measured using a voltage sensor, in the fresh tray consisting of 28 parallel-connected cells according to the aging current profile. It can be seen that the maximum voltage difference was 0.015 V, which is very small.

Experimental Result of the Retired Series/Parallel Battery Pack
The retired series/parallel battery pack was aged according to the aging profile presented in Figure 2. Moreover, there was no information about the aging process and current status datasets of the selected retired series/parallel battery pack. Therefore, a full discharge/charge test was conducted to determine the current status of the aged battery pack. Figure 4 shows 12 discharging/charging

Experimental Result of the Retired Series/Parallel Battery Pack
The retired series/parallel battery pack was aged according to the aging profile presented in Figure 2. Moreover, there was no information about the aging process and current status datasets of the selected retired series/parallel battery pack. Therefore, a full discharge/charge test was conducted to determine the current status of the aged battery pack. Figure 4 shows 12 discharging/charging voltages for each representative cell, measured using a voltage sensor, in the retired series/parallel battery pack consisting of 28 parallel-connected cells in a full charging and discharging 1C-rate experiment according to the cut-off voltage range of the fresh battery pack. It can be confirmed that the 7th, 10th, and 12th series stages with 28 parallel connections were overdischarged owing to the large voltage difference. Moreover, the 7th series stage was shorted to 0 V. Therefore, the retired series/parallel battery pack with a voltage imbalance should be used with consideration of the condition of each cell. Moreover, the discharge/charge range of the retired series/parallel battery pack in which the voltage imbalance occurred was narrowed to prevent the overdischarge/overcharge. It is necessary to calculate the discharge capacity using a fully discharging/charging test to determine the accurate SOH of the retired series/parallel battery pack. However, a fully discharging/charging test was impossible owing to the reduced voltage range caused by the voltage imbalance. Therefore, the SOH of the retired series/parallel battery pack within the voltage imbalance should be estimated using the information from the reduced discharge/charge voltage range. . Difference in 12 cell voltages of each representative cell, measured using a voltage sensor, in the fresh battery pack consisting 28 parallel-connected cells according to the discharging/charging current profile for aging (4C-rate).

Experimental Result of the Retired Series/Parallel Battery Pack
The retired series/parallel battery pack was aged according to the aging profile presented in Figure 2. Moreover, there was no information about the aging process and current status datasets of the selected retired series/parallel battery pack. Therefore, a full discharge/charge test was conducted to determine the current status of the aged battery pack. Figure 4 shows 12 discharging/charging voltages for each representative cell, measured using a voltage sensor, in the retired series/parallel battery pack consisting of 28 parallel-connected cells in a full charging and discharging 1C-rate experiment according to the cut-off voltage range of the fresh battery pack. It can be confirmed that the 7th, 10th, and 12th series stages with 28 parallel connections were overdischarged owing to the large voltage difference. Moreover, the 7th series stage was shorted to 0 V. Therefore, the retired series/parallel battery pack with a voltage imbalance should be used with consideration of the condition of each cell. Moreover, the discharge/charge range of the retired series/parallel battery pack in which the voltage imbalance occurred was narrowed to prevent the overdischarge/overcharge. It is necessary to calculate the discharge capacity using a fully discharging/charging test to determine the accurate SOH of the retired series/parallel battery pack. However, a fully discharging/charging test was impossible owing to the reduced voltage range caused by the voltage imbalance. Therefore, the SOH of the retired series/parallel battery pack within the voltage imbalance should be estimated using the information from the reduced discharge/charge voltage range.

7S Shorted to 0V
7S, 10S, 12S Over-discharged Figure 4. Difference in 12 cell voltages of each representative cell, measured using a voltage sensor, in the retired series/parallel battery pack consisting of 28 parallel-connected cells according to the discharging/charging current profile for aging (1C-rate). . Difference in 12 cell voltages of each representative cell, measured using a voltage sensor, in the retired series/parallel battery pack consisting of 28 parallel-connected cells according to the discharging/charging current profile for aging (1C-rate).

Aging Cycle Test of the NCA Cell
It is necessary to determine the calendar-life data of the battery from the fresh state to the retired state to design the aging model of the retired series/parallel battery pack. However, there was no data of the battery pack during the aging period. Therefore, three fresh NCA cells were used for the aging cycle test. The experiments were performed on a three-channel MACCOR 4300K battery test system (Maccor, KR) with a current range of 0-15 A and voltage range of 0-5 V, as shown in Figure 5. The voltage and current accuracies were ±0.02% and ±0.05%, respectively. The data acquisition system had a time resolution of maximum 10 ms. The aging cycle test was performed in the temperature chamber within 25 • C. The 1C-rate (2 A) constant current (CC)-constant voltage (CV) full discharge/charge profile used to age the cells is presented in Figure 6a. The CC-CV full discharge/charge aging test was performed 100 times, and the capacity degradation for a particular slope according to the aging cycle test is presented in Figure 6b. The cells datasets of the aging cycle test were used for the capacity degradation modeling to estimate the SOH of the aged battery pack. voltage and current accuracies were ±0.02% and ±0.05%, respectively. The data acquisition system had a time resolution of maximum 10 ms. The aging cycle test was performed in the temperature chamber within 25 °C. The 1C-rate (2 A) constant current (CC)-constant voltage (CV) full discharge/charge profile used to age the cells is presented in Figure 6a. The CC-CV full discharge/charge aging test was performed 100 times, and the capacity degradation for a particular slope according to the aging cycle test is presented in Figure 6b. The cells datasets of the aging cycle test were used for the capacity degradation modeling to estimate the SOH of the aged battery pack.

Original IC Curve
An IC curve is widely used as a practical factor to analyze the capacity degradation of a battery and can be obtained from the discharge/charge state under the CC condition. Moreover, the IC is calculated using the capacity difference divided by the voltage difference for the same time, as shown in Equation (1). Figure 7 is an IC curve of cell no. 1 in the discharge state. Moreover, the sampling time of the IC curve dataset was 1 s. However, it was difficult to obtain the peaks based on the IC method because the parameter values had measurement noise.
in Equation (1). Figure 7 is an IC curve of cell no. 1 in the discharge state. Moreover, the sampling time of the IC curve dataset was 1 s. However, it was difficult to obtain the peaks based on the IC method because the parameter values had measurement noise.

Smoothing of the IC Curve
As is commonly done, a first-order LPF was used to smooth the measurement noise of the IC curve [19]. Equation (2) represents the discrete form of a first-order LPF, where yn−1 is the previous result of the LPF and τ is the value of the time constant, which determines the cut-off frequency (fcut) of the LPF, as presented in Equation (3). Here, xn is the input value of the LPF and ts is the sampling time of the input value. As τ increases, fcut becomes lower. Moreover, the IC curve becomes smoother, as presented in Figure 8. Because ts was 1 s, the fcut was lower than 1 Hz. The noise of the IC curve was properly removed at 100τ, when the τ was applied in the range from 10 to 150. Therefore, a firstorder LPF of 100τ was used to remove the measurement noise in the IC curve in this research.

Smoothing of the IC Curve
As is commonly done, a first-order LPF was used to smooth the measurement noise of the IC curve [19]. Equation (2) represents the discrete form of a first-order LPF, where y n−1 is the previous result of the LPF and τ is the value of the time constant, which determines the cut-off frequency (f cut ) of the LPF, as presented in Equation (3). Here, x n is the input value of the LPF and t s is the sampling time of the input value. As τ increases, f cut becomes lower. Moreover, the IC curve becomes smoother, as presented in Figure 8. Because t s was 1 s, the f cut was lower than 1 Hz. The noise of the IC curve was properly removed at 100τ, when the τ was applied in the range from 10 to 150. Therefore, a first-order LPF of 100τ was used to remove the measurement noise in the IC curve in this research.

Selection of the IC Curve Peak Point
The voltage operation range of the 12S28P NCA retired series/parallel battery pack was narrowed by the voltage imbalance to prevent the overdischarge/overcharge, as presented in Figure  9. The voltage operation range of the 12S28P NCA fresh battery pack was from 36-49.5 V, and for the retired series/parallel battery pack, it was from 40.8 V to 40 V. In the 12 series-connected stages, the normal discharge/charge range was from 3.6 V to 3.8 V. Therefore, the SOH of the retired series/parallel battery pack should be calculated using the information in the voltage operation range of 3.8-3.6 V. Figure 10 shows the IC curve peak point of a fresh cell no. 1 for 100 cycles according to

Selection of the IC Curve Peak Point
The voltage operation range of the 12S28P NCA retired series/parallel battery pack was narrowed by the voltage imbalance to prevent the overdischarge/overcharge, as presented in Figure 9. The voltage operation range of the 12S28P NCA fresh battery pack was from 36-49.5 V, and for the retired series/parallel battery pack, it was from 40.8 V to 40 V. In the 12 series-connected stages, the normal discharge/charge range was from 3.6 V to 3.8 V. Therefore, the SOH of the retired series/parallel battery pack should be calculated using the information in the voltage operation range of 3.8-3.6 V. Figure 10 shows the IC curve peak point of a fresh cell no. 1 for 100 cycles according to the aging cycle test. It was confirmed that a total of three peak points and two valley points occurred in the full discharge/charge range. However, all five points cannot be used to estimation of SOH due to voltage imbalance in the battery pack. As the maximum voltage difference increases, the number of available IC curve peak points decreases. The maximum chargeable voltage and the maximum dischargeable voltage decrease by the magnitude of the maximum voltage difference. So, the retired battery pack used in this paper has a 0.4V of maximum voltage difference and voltage operation range of 3.6-3.8V, as shown in Figure 9. And, the available IC curve peak point is only valley point 1. Furthermore, the IC curve peak point changed with a decreasing capacity, depending on the type of battery. However, the storage capability information of most lithium-ion batteries were determined by the degree of the redox reaction at the negative electrode of the battery [24,25]. The same applied to the NCA series batteries used in this paper. With the five points shown in Figure 10, the storage capability range at the cathode of the battery was peak point 1 and valley point 1 [24,25]. So, the valley point 1 in the IC curve occurring between 3.6 V and 3.7 V satisfied the common discharge/charge voltage operation range of the retired series/parallel battery pack and available factor of capacity fading, as shown in Figure 10. In this study, the IC curve valley point 1 of fresh cell no. 1 was used as a factor for the capacity-degrading model of the retired series/parallel battery pack. Figure 11a is the IC curve valley point 1 of fresh cell no. 1 according to the aging cycle test with 100 cycles. It was confirmed that the IC curve valley point 1 increased with the number of aging cycle tests.

Analysis of the Partial Capacity in the Short Voltage Range
The analysis of the tendency of the directly measured capacity value was performed in the same voltage range to compare the tendency of an IC curve peak point with that found in Section 4.1. Figure 12 shows the tendency of the directly measured capacity value of cell no. 1 within the voltage range from 3.6 V to 3.7 V. It was confirmed that the partial capacity value of cell no. 1 in this short voltage range decreased with the number of aging cycle tests. This tendency indicated that the partial capacity value in the short voltage range decreased as the capacity value in the full operating voltage range decreased. Therefore, the directly measured partial capacity value of cell no. 1 in the short voltage range was also suitable for use as a factor of the capacity fading model. The analysis of the tendency of the directly measured capacity value was performed in the same voltage range to compare the tendency of an IC curve peak point with that found in Section 4.1. Figure 12 shows the tendency of the directly measured capacity value of cell no. 1 within the voltage range from 3.6 V to 3.7 V. It was confirmed that the partial capacity value of cell no. 1 in this short voltage range decreased with the number of aging cycle tests. This tendency indicated that the partial capacity value in the short voltage range decreased as the capacity value in the full operating voltage range decreased. Therefore, the directly measured partial capacity value of cell no. 1 in the short voltage range was also suitable for use as a factor of the capacity fading model.

Linear Regression Method
Linear regression analysis is a statistical analysis method for identifying the relationship between independent and dependent variables [26,27]. A linear regression model with one independent variable is called a simple linear regression model, and with two or more independent variables, it is called a multiple linear regression model. In this study, a simple linear regression model was used to model the capacity degradation of the retired series/parallel battery pack. To determine the applicability of linear regression analysis, correlation analysis should be performed to confirm the relationship between the independent and dependent variables. Correlation analysis is a method of determining the linear relationship between independent and dependent variables. Moreover, the correlation coefficient, r, is a statistic that shows the linear correlation between two variables, as expressed in Equation (4) [26]. The x and y denote independent and dependent variables, respectively, ̅ and are the average values of the independent and dependent variables, respectively, and n is the number of data points.

Linear Regression Method
Linear regression analysis is a statistical analysis method for identifying the relationship between independent and dependent variables [26,27]. A linear regression model with one independent variable is called a simple linear regression model, and with two or more independent variables, it is called a multiple linear regression model. In this study, a simple linear regression model was used to model the capacity degradation of the retired series/parallel battery pack. To determine the applicability of linear regression analysis, correlation analysis should be performed to confirm the relationship between the independent and dependent variables. Correlation analysis is a method of determining the linear relationship between independent and dependent variables. Moreover, the correlation coefficient, r, is a statistic that shows the linear correlation between two variables, as expressed in Equation (4) [26]. The x and y denote independent and dependent variables, respectively, x and y are the average values of the independent and dependent variables, respectively, and n is the number of data points.
A correlation coefficient, r, value close to the absolute value of 1 indicates that there is a strong relationship between the independent and dependent variable. Equation (5) represents a simple linear regression model [26]. The above equation expresses the effect of an independent variable, x, on a dependent variable, y. The influence of the independent variable, x, on the dependent variable, y, is expressed as an equation. The effect of x is represented by the magnitude and sign of the regression coefficient, β. The regression coefficient, α, implies that the change in y does not affect the change in x. Moreover, the error term ε includes all of the influences except the influence of x on y.
Because the regression curve given in Equation (5) cannot represent the data of the entire population, the fitted regression equation given in Equation (6) is used [26]. The fitted regression curve is the linear regression line of the measurement from the specimen. The values of the errors, e i , and sum of squared errors (SSE) represent the difference between the fitted regression line,ŷ i , and actual data, y i , as expressed in Equations (5) and (6) [26]. Calculating a and b to minimize the SSE leads to the linear regression equation that best represents the relationship between x and y. The least-squares method is used to minimize the SSE, as presented in Equations (7) and (8) [26]. The regression coefficients, α and β, are calculated to minimize the SSE using the partial differentiation method. Following the regression coefficient calculation, the coefficient of determination, R 2 , is calculated, as given in Equations (6), (10), and (11) [26]. The sum of squares (SSR) represents the unexplained variance from the linear regression equation of the total variance. The coefficient of determination is a variable that can be explained using the estimated regression equation of the total variance of the variables. The coefficient of determination has a value between 0 and 1 A coefficient of determination, R 2 , close to the value of 1 indicates a strong relationship between the variables and estimated regression equation.
The results of the linear regression modeling according to the IC curve peak points and partial capacity of cell no. 1 are shown in Figures 13 and 14. Figures 13a and 14a shows the IC curve peak points and partial capacity of cell no. 1, respectively, according to the SOH changes. In Sections 4.1 and 4.2, it was confirmed that the IC curve valley point 1 and partial capacity of cell no. 1 increased and decreased as the capacity decreased according to the aging cycle test, as shown in Figures 11b  and 12b. The capacity decrease was indicative of the SOH decrease. The SOH can be calculated using Equation (12), where Capacity current is the discharge capacity of the battery in the current state and Capacity fresh is the discharge capacity of the battery in the fresh state.
Therefore, the changes in the IC curve valley point 1 and the partial SOC of cell no. 1 according to the SOH decreases presented in Figures 11b and 12b and Equation (11) are shown in Figures 13a  and 14a. The datasets in Figures 13a and 14a were used to estimate the linear regression equation. First, the correlation coefficient, r, was calculated using Equation (4), and the results were 0.9796 and 0.9864. It was confirmed that the partial capacity of cell no.1 had a slightly higher correlation with SOH than the IC curve valley point 1 of cell no. 1 had with SOH. However, the resulting values of both the parameters were large. Therefore, the independent variables of the IC curve valley point 1 and partial capacity of cell no. 1 and dependent variable SOH had a strong correlation, as seen through the linear regression analysis. Moreover, Figures 13b and 14b show a linear equation for these datasets. It was confirmed that the IC curve valley point 1 and the partial capacity of cell No.1 had negative and positive regression coefficients, respectively, as the SOH decreased.

SOH current =
Capacity current Capacity f resh The coefficient of determination, R 2 , was calculated using Equation (12), and the results were 0.9774 and 0.9849. It was confirmed that the partial capacity of cell no. 1 had a slightly larger coefficient of determination, R 2 , than that of the IC curve valley point 1 of cell no. 1. These values of the coefficient of determination, R 2 , indicating that the relationship between the variables and estimated regression. The linearity of x and y is high according to the R 2 is close to 1. Therefore, it was confirmed that the derived linear regression equations of the two methods were suitable for use as battery aging models without further statistical verification. Therefore, the changes in the IC curve valley point 1 and the partial SOC of cell no. 1 according to the SOH decreases presented in Figure 11b and Figure 12b and Equation (11) are shown in Figure  13a and Figure 14a. The datasets in Figure 13a and Figure 14a were used to estimate the linear regression equation. First, the correlation coefficient, r, was calculated using Equation (4), and the results were 0.9796 and 0.9864. It was confirmed that the partial capacity of cell no.1 had a slightly higher correlation with SOH than the IC curve valley point 1 of cell no. 1 had with SOH. However, the resulting values of both the parameters were large. Therefore, the independent variables of the IC curve valley point 1 and partial capacity of cell no. 1 and dependent variable SOH had a strong correlation, as seen through the linear regression analysis. Moreover, Figure 13b and Figure 14b show a linear equation for these datasets. It was confirmed that the IC curve valley point 1 and the partial capacity of cell No.1 had negative and positive regression coefficients, respectively, as the SOH decreased.
The coefficient of determination, R 2 , was calculated using Equation (12), and the results were 0.9774 and 0.9849. It was confirmed that the partial capacity of cell no. 1 had a slightly larger coefficient of determination, R 2 , than that of the IC curve valley point 1 of cell no. 1. These values of the coefficient of determination, R 2 , indicating that the relationship between the variables and estimated regression. The linearity of x and y is high according to the R 2 is close to 1. Therefore, it was confirmed that the derived linear regression equations of the two methods were suitable for use as battery aging models without further statistical verification. Figure 15 exhibits the calculation step of the SOH using a linear regression model. The IC curve was calculated in the discharge interval, and the IC curve peak value was extracted in the terminal  Figure 15 exhibits the calculation step of the SOH using a linear regression model. The IC curve was calculated in the discharge interval, and the IC curve peak value was extracted in the terminal voltage range between 3.6 V and 3.7 V. Moreover, the extracted IC curve peak point was used as the input value of the linear regression model to calculate the SOH, as shown in Figure 12b. The directly measured partial capacity was also extracted in the terminal voltage range between 3.6 V and 3.7 V. The comparison of the linear regression model of cell no. 1 and the real SOH according to the IC curve peak point of cell nos. 2 and 3 as a reference value is shown in Figures 16 and 17. Furthermore, the numerical comparison results of the SOH obtained using the linear regression model of cell no. 1 and a real SOH as a reference value are listed in Tables 1 and 2. The error of the SOH was calculated using Equation (13). SOH real was calculated using the real capacity value obtained using an off-line experiment. SOH model was calculated using the linear regression model. Figure 18 shows the error in the SOH using the IC curve peak point and partial capacity according to the aging cycle. It was confirmed that the maximum errors in the SOH with the IC curve peak point and partial SOC were 0.9771% and 0.8653%. All the accuracy is within 1%. Furthermore, it was confirmed that the average error in the estimated SOH was also within 1%, as shown in Figure 18. Therefore, it was confirmed that both the designed aging models with the IC curve peak point and partial capacity exhibited a good accuracy in the cell-based verification.       Figure 19 shows the IC curve valley point 1 of the 12 series-connected stages of the retired series/parallel battery pack using the datasets in Figure 4. The 7S series stage of the short circuit was excluded. These datasets in Figure 19 were used to estimate the SOH of the retired series/parallel battery pack based on the capacity fading model shown in Figure 13b. The value of the off-line current capacity was needed to verify the estimated SOH value of the retired series/parallel battery pack. However, it was impossible to perform a full discharge/charge of a retired series/parallel battery pack with a voltage imbalance. Therefore, this retired pack was disassembled, as shown in Figure 20. Following the disassembly of the retired series/parallel battery pack, a full discharge/charge test was conducted to obtain the value of the current capacity. These current capacities and SOHs of the 12 series-connected stages of the disassembled retired series/parallel battery pack are presented in Figure 21. It was confirmed that an imbalance in the SOH was caused by the imbalance in the voltage in the 12 series-connected stages.  Figure 19 shows the IC curve valley point 1 of the 12 series-connected stages of the retired series/parallel battery pack using the datasets in Figure 4. The 7S series stage of the short circuit was excluded. These datasets in Figure 19 were used to estimate the SOH of the retired series/parallel battery pack based on the capacity fading model shown in Figure 13b. The value of the off-line current capacity was needed to verify the estimated SOH value of the retired series/parallel battery pack. However, it was impossible to perform a full discharge/charge of a retired series/parallel battery pack with a voltage imbalance. Therefore, this retired pack was disassembled, as shown in Figure 20. Following the disassembly of the retired series/parallel battery pack, a full discharge/charge test was conducted to obtain the value of the current capacity. These current capacities and SOHs of the 12 series-connected stages of the disassembled retired series/parallel battery pack are presented in Figure 21. It was confirmed that an imbalance in the SOH was caused by the imbalance in the voltage in the 12 series-connected stages.   Figure 19 shows the IC curve valley point 1 of the 12 series-connected stages of the retired series/parallel battery pack using the datasets in Figure 4. The 7S series stage of the short circuit was excluded. These datasets in Figure 19 were used to estimate the SOH of the retired series/parallel battery pack based on the capacity fading model shown in Figure 13b. The value of the off-line current capacity was needed to verify the estimated SOH value of the retired series/parallel battery pack. However, it was impossible to perform a full discharge/charge of a retired series/parallel battery pack with a voltage imbalance. Therefore, this retired pack was disassembled, as shown in Figure 20. Following the disassembly of the retired series/parallel battery pack, a full discharge/charge test was conducted to obtain the value of the current capacity. These current capacities and SOHs of the 12 series-connected stages of the disassembled retired series/parallel battery pack are presented in Figure 21. It was confirmed that an imbalance in the SOH was caused by the imbalance in the voltage in the 12 series-connected stages.    The maximum and minimum values of the SOH were 0.96 and 0.88, respectively, and the maximum difference in the SOH was 0.08. By converting this value of the SOH error to a percentage, it was confirmed that there was a very large SOH difference of 8%. If this retired series/parallel battery pack continues to be used, the speed of the SOH imbalance and voltage imbalance will increase. Therefore, it is important to estimate the current SOH of the series-connected stage of the retired series/parallel battery pack. Figures 22 and 23 show the comparison of the linear regression model of cell no. 1 and datasets of the IC curve peak point and partial capacity according to the real SOH value of the 12 seriesconnected stages of the retired series/parallel battery pack. Moreover, Tables 3 and 4 show the numerical results of Figures 22 and 23. It was confirmed that the maximum error in the SOH from the IC curve peak point method and partial capacity method was 0.7867% and 4.3288%, respectively, and the maximum error in the SOH from the partial capacity method was extremely large. Figure 24 shows the error in the SOH obtained with the IC curve peak point and partial capacity according to the number of the 12 series-connected stages of the retired series/parallel battery pack. It was confirmed that the average errors in the estimated SOH using the IC curve peak point method and partial capacity method were 0.4863% and 1.6309%, respectively. The SOH error value with the IC curve peak point method was within 1%, but the result for the partial capacity method had a large value of over 1% of the target of the retried battery pack. These results suggest that the directly measured partial capacity method was not suitable for the different aging conditions of the battery. The aging profiles of the unit fresh NCA cell and fresh NCA series/parallel battery pack were different, as shown in Figures 2 and 5. Moreover, the capacity fading model for the SOH estimation was built with the unit cell information, but the retired series/parallel battery pack had a different aging history than of the unit cell. Also, the retired series/parallel battery pack had a difference in the depth of discharge between the series terminals owing to the voltage imbalance, as shown in Figure  9. It is necessary that the aging model is designed with the data under the same aging condition for the application to estimate the SOH with the partial capacity method, in which the battery is used. However, the IC curve peak point method closely matched with the phase transition during the Li + intercalation and de-intercalation process [18]. This process is directly related to the capacity fading of the battery. Therefore, the IC curve peak point method was suitable for the various aging conditions of the battery. Thus, the SOH estimation method using the capacity fading model via the linear regression analysis according to the IC curve peak point and SOH in the narrowed discharge/charge operation voltage range was sufficient. The maximum and minimum values of the SOH were 0.96 and 0.88, respectively, and the maximum difference in the SOH was 0.08. By converting this value of the SOH error to a percentage, it was confirmed that there was a very large SOH difference of 8%. If this retired series/parallel battery pack continues to be used, the speed of the SOH imbalance and voltage imbalance will increase. Therefore, it is important to estimate the current SOH of the series-connected stage of the retired series/parallel battery pack. Figures 22 and 23 show the comparison of the linear regression model of cell no. 1 and datasets of the IC curve peak point and partial capacity according to the real SOH value of the 12 series-connected stages of the retired series/parallel battery pack. Moreover, Tables 3 and 4 show the numerical results of Figures 22 and 23. It was confirmed that the maximum error in the SOH from the IC curve peak point method and partial capacity method was 0.7867% and 4.3288%, respectively, and the maximum error in the SOH from the partial capacity method was extremely large. Figure 24 shows the error in the SOH obtained with the IC curve peak point and partial capacity according to the number of the 12 series-connected stages of the retired series/parallel battery pack. It was confirmed that the average errors in the estimated SOH using the IC curve peak point method and partial capacity method were 0.4863% and 1.6309%, respectively. The SOH error value with the IC curve peak point method was within 1%, but the result for the partial capacity method had a large value of over 1% of the target of the retried battery pack. These results suggest that the directly measured partial capacity method was not suitable for the different aging conditions of the battery. The aging profiles of the unit fresh NCA cell and fresh NCA series/parallel battery pack were different, as shown in Figures 2 and 5. Moreover, the capacity fading model for the SOH estimation was built with the unit cell information, but the retired series/parallel battery pack had a different aging history than of the unit cell. Also, the retired series/parallel battery pack had a difference in the depth of discharge between the series terminals owing to the voltage imbalance, as shown in Figure 9. It is necessary that the aging model is designed with the data under the same aging condition for the application to estimate the SOH with the partial capacity method, in which the battery is used. However, the IC curve peak point method closely matched with the phase transition during the Li + intercalation and de-intercalation process [18]. This process is directly related to the capacity fading of the battery. Therefore, the IC curve peak point method was suitable for the various aging conditions of the battery. Thus, the SOH estimation method using the capacity fading model via the linear regression analysis according to the IC curve peak point and SOH in the narrowed discharge/charge operation voltage range was sufficient.

Results of the Retired Cell Datasets
In Sections 5.1 and 5.2, the datasets used for verification came from batteries with the aging period before 80% SOH. Batteries with the aging period before 80% SOH are not retired batteries. The battery pack used in Section 5.2 was a short-circuited battery due to the large cell-to-cell voltage imbalance, but was not sufficient to verify the battery SOH values after 0.8. The additional verification was performed with the 600 CC-CV cycle test using retired cell datasets to verify the aging period after 80% SOH of the proposed SOH estimation model, as shown in Figure 25. Furthermore, Figures 26 and 27 show the comparison between the linear regression model of cell no. 1 and datasets of the IC curve peak point and partial capacity according to the real SOH value of the retired cell. Furthermore, the numerical comparison results of the SOH obtained using the linear regression model of cell no. 1 and a real SOH as a reference value are listed in Tables 5 and 6. It was confirmed that the maximum error in the SOH from the IC curve peak point method and partial capacity method were 1.303744% and 3.625032%, respectively, and the maximum error in the SOH from the partial capacity method was extremely large. Figure 18 shows the error in the SOH obtained with the IC curve peak point and partial capacity according to the number of the aging cycle. It was confirmed that the average errors in the estimated SOH using the IC curve peak point method and partial capacity method were 0.4967% and 2.0158%, respectively. The SOH error value with the IC curve peak point method was within 0.5%, but the partial capacity method has a large error value of over 2% of the target of the retired battery pack. These results suggest that the directly measured partial capacity method was not suitable for the retired battery below a SOH of 0.8. Furthermore, the IC curve peak point method was suitable for the retired battery of the aging period after 80% SOH. Thus, the SOH estimation method using the capacity fading model via the linear regression analysis according to the IC curve peak point and SOH for recycling of retired battery was effective.

Conclusions
This paper proposed a SOH estimation method for a retired series/parallel battery pack with a voltage imbalance. The voltage imbalance and narrowed voltage operation range were confirmed using the discharge/charge tests of the 12S28P NCA retired series/parallel battery pack. Moreover, an IC curve peak point method was proposed to estimate the SOH using the information of the narrowed voltage operation range. The linear regression analysis method was performed with the IC curve peak point and SOH value based on the capacity information of the fresh cell according to a 100-cycle aging test to build the capacity fading model. First, to verify the capacity fading model using linear regression analysis, the datasets of other cells according to the 100-cycle aging test were used. All the errors in the SOH were within 1%, and it was confirmed via the linear regression model that the capacity fading model had a good accuracy. The directly measured partial capacity analysis in a short voltage range was performed to be compared with the result of the IC curve peak point. Moreover, all the errors in the SOH were within 1%, and it was confirmed by the linear regression model that the partial capacity model had a good accuracy using cell-based verification. Second, to verify the capacity fading model of cell no. 1 via linear regression, datasets of the 12 series-connected stages of the retired series/parallel battery pack were used. Moreover, the result of the error in the SOH was within 1%. However, the result of the error in SOH with the partial capacity method was over 1%. Therefore, it was confirmed by the linear regression analysis that the capacity fading model with the IC curve peak point was suitable for SOH estimation of a retired series/parallel battery pack with a voltage imbalance. If the SOH was estimated in the units of a retired series/parallel battery pack with voltage imbalances without disassembling it, then the time and cost of the capacity-based screening for recycling the retired series/parallel battery pack were effectively reduced. Moreover, it was easy to determine whether the retired series/parallel battery pack could be used. Finally, to verify the capacity fading model by linear regression analysis after the aging period of SOH 80%, the datasets of other cells according to the 600-cycle aging test were used. The result of the error in SOH with the IC curve peak point method was within 0.5%. However, the result of the error in SOH with the partial capacity method was over 2%. Therefore, it was confirmed via the linear regression analysis that the capacity fading model with the IC curve peak point was suitable for SOH estimation of a retired cell after the aging period of 80% SOH. The future works will require verification in the aging intervals after the 80% SOH value of the retired battery pack. This future work will for the SOH estimation of recycled retired series/parallel battery packs.