Parameter Sensitivity Analysis for Fractional-Order Modeling of Lithium-Ion Batteries

This paper presents a novel-fractional-order lithium-ion battery model that is suitable for use in embedded applications. The proposed model uses fractional calculus with an improved Oustaloup approximation method to describe all the internal battery dynamic behaviors. The fractional-order model parameters, such as equivalent circuit component coefficients and fractional-order values, are identified by a genetic algorithm. A modeling parameters sensitivity study using the statistical Multi-Parameter Sensitivity Analysis (MPSA) method is then performed and discussed in detail. Through the analysis, the dynamic effects of parameters on the model output performance are obtained. It has been found out from the analysis that the fractional-order values and their corresponding internal dynamics have different degrees of impact on model outputs. Thus, they are considered as crucial parameters to accurately describe a battery’s dynamic voltage responses. To experimentally verify the accuracy of developed fractional-order model and evaluate its performance, the experimental tests are conducted with a hybrid pulse test and a dynamic stress test (DST) on two different types of lithium-ion batteries. The results demonstrate the accuracy and usefulness of the proposed fractional-order model on battery dynamic behavior prediction.


Introduction
Nowadays, electric vehicles (EVs), including hybrid electric vehicles (HEVs) and fuel-cell-based HEVs (FCHEV), are considered to be an effective way to achieve significant fuel consumption and carbon emission reductions [1]. Compared with fuel cells [2,3], lithium-ion batteries are considered the principal candidate for short-and medium-range EVs due to their advantages of high energy density, longer cycle life, easier management, safer chemistry and broad operation temperature range. To ensure reliable and safe operation, a battery management system (BMS) is necessary for the lithium-ion battery packs in EVs [4,5]. An important requirement for the BMS is that it should be able to provide an accurate estimation of state of charge (SOC) and state of health (SOH) for the battery. In general, designing an accurate BMS estimation algorithm for SOC and SOH requires a battery model that can accurately describe the battery terminal voltage behaviors in both steady-state and dynamic transient conditions.
Generally speaking, battery models presented in literature can be divided in two categories: (1) equivalent-circuit models (ECMs), and (2) physics-based electrochemical models. ECMs use equivalent voltage sources, resistors and capacitors to characterize the battery output voltage feature. Due to their mathematical simplicity and acceptable precision, ECMs are widely used in BMS for model, large computations and memories are needed due to its complex electrochemical PDEs, thus it cannot be directly used in BMS for real-time applications. In addition, most fractional calculus approaches used in these fractional models are complex and non-recursive, which leads to a heavy computational burden. Furthermore, a common negligence of these fractional models in the literature is that, the influence of the parameters, especially the effect of the values of fractional calculus on the battery model performance, has never been analyzed or discussed for the battery model development. In fact, by definition, the fractional-order model can show various voltage dynamic characteristics by changing the values of fractional-order parameters. Thus, the lack of knowledge about the sensitivity and dynamic effect of each fractional-order parameter can lead to an inaccurate fractional-order model output result.
In this paper, a novel-fractional-order model for lithium-ion batteries is presented at first. The proposed fractional-order model can accurately describe all the lithium-ion battery internal electrochemical behaviors, including battery Ohmic losses, lithium-ion migration through the solid electrolyte interface (SEI) layer, charge transfer kinetics, double layer effects, and lithium-ion diffusion processes using free fractional-orders elements. The proposed fractional-order model is developed and solved using an improved Oustaloup approximation algorithm [22,23], which can recursively approximate free fractional-orders transfer function from a set of integer-order filters for lithium-ion battery models. For the parameter identification process of the proposed model, a genetic algorithm is used. Secondly, a full analysis of modeling parameter sensitivity, including circuit components coefficients and fractional-order values, is performed. The experimental validation of the proposed model is then performed on two different types of lithium-ion batteries. In order to highlight the usefulness of the developed model, a comparison between the prediction results of the proposed fractional-order model and a classic second-RC model is also presented. Experimental results show that, the proposed fractional-order battery model can accurately describe the battery terminal voltage with dynamic transient behaviors. In addition, the parameter sensitivities analyses provide insights into the influence of fractional-order parameters, which is important for accurate description of dynamic behaviors of battery model. The novel contributions of this paper can be summarized as follows:

‚
All the battery internal dynamic behaviors, including SEI layer, charge transfer kinetics, double layer effects, and lithium-ion diffusion processes are accurately captured by simple free fractional-orders elements in the proposed model. With the proposed model, the battery voltage dynamics of charging and discharging can be precisely reproduced without using complex electrochemical equations. Thus, the proposed model can accurately predict dynamic behaviors inside the battery by using much less computational resources compared to classic electrochemical equations-based models. This feature makes the proposed model particularly suitable for implementation in real time embedded applications.

‚
The presented new fractional calculus method solves the boundary fitting problems using the innovative improved Oustaloup recursive algorithm, which, to the best of our knowledge, has never been used for lithium-ion battery models in the literature. This novel approach combines a boundary fitting optimization algorithm with a traditional Oustaloup recursive approximation method, thus the calculation accuracy and efficiency are both satisfied. In this case, a satisfactory fitting result can be efficiently achieved in the whole pre-specified frequency range. Compared to other fractional calculus approximation algorithms, the proposed efficient approach in this paper provides a more precise approximation of free fractional-orders derivative by a series of integer-order filters. From our point of view, this approach is a novel contribution in approximating lithium-ion battery free fractional-order models.

‚
The modeling parameter sensitivities, especially the effect of the values of fractional calculus on the battery model performance, are analyzed in this paper. This analysis provides insights into the influence of fractional-order parameters, and further shows which internal dynamic behaviors have more significant effects on battery terminal voltage.

Fractional-Order Models
In order to have a general understanding of fractional-order system, more definitions of fractional calculus are given in the Appendix A. In this section, the proposed novel fractional-order lithium-ion battery model is presented based on the battery electrochemical impedances approach [19]. The mathematical modelling equations and the fractional transfer functions are described and solved using an improved fractional-order system approximation method [22,23]. Finally, the parameters identification process is shown as well.

An Improved Oustaloup Recursive Approximation Algorithm
The numerical integer-order approach of derivation for fractional transfer functions is critical in fractional-order systems. As an efficient approach, the Oustaloup recursive approximation algorithm [22] is widely used to approximate fractional-order transfer functions, and hence selected in this paper. The Oustaloup algorithm uses a set of integer-order filters for approximation of free fractional-order derivatives in a certain frequency range. However, it has been indicated that, in the Oustaloup algorithm, the numerical fitting accuracy is not satisfactory in the regions near the lower bound ω a and upper bound ω b .
To overcome this boundary fitting drawback of the Oustaloup algorithm, an improved filter design is proposed in [23], which combines a boundary fitting optimization algorithm with traditional Oustaloup recursive approximation method, thus maintain simultaneously the calculation accuracy and efficiency.
In this paper, the notation s in Laplace domain is used to describe fractional differentiation: Then the general single-input single-output (SISO) system can be expressed in a fractional-order transfer function form: where µ 0 , µ 1 , . . . , µ m and ν 0 , ν 1 , . . . , ν n are system parameters, α 0 ă α 1 ă . . . ă α m and β 1 ă β 2 ă¨¨¨ă β n are the fractional orders, which are all real numbers. Based on the first-order approximation of Taylor series expansion of Equation (A23) in Appendix A, the improved filter is thus given by: A more detailed derivation can be found in Appendix A. In this paper, this improved Oustaloup approximation is applied to the fractional-order transfer functions of the proposed lithium-ion battery model.

Structure of Equivalent Fractional-Order Electrochemical Impedance Model
The fractional-order electrochemical impedance modeling approach of a lithium-ion battery consists of different dynamic phenomena inside the battery: lithium-ion migration through the solid electrolyte interface (SEI) layer, activation kinetics in both negative and positive electrodes, double layer effects at the interfaces of electrolytes, and lithium-ion diffusion processes in the active material of the electrodes. In order to characterize these phenomena, the electrochemical impedance spectroscopy (EIS) method is generally used. The lithium-ion battery internal dynamics phenomena, the corresponding EIS Nyquist diagram and the proposed equivalent fractional-order impedance circuit diagram are shown in Figure 1.  It can be seen from the Figure 1 that, the EIS Nyquist diagram can be divided into four regions [24][25][26]:


Region A can be represented by a single Ohmic resistance , this area corresponds to battery Ohmic losses, including the resistances of current collectors, active material, electrolyte and separator.  High frequency region B, which can be represented by a solid electrolyte interface (SEI) resistance , a SEI layer capacitance , and a fractional element δ , corresponds to the lithium-ion migration through the SEI film layer.  Medium frequency region C, which can be represented by a charge transfer resistance , a double layer capacitance , and a fractional element δ , corresponds to the activation kinetics and double layer effects at the interface of electrolytes;  Low frequency region D, which can be represented by a Warburg element and a fractional element δ , corresponds to lithium-ion diffusion processes in the active material of the electrodes.
Based on the analysis above, each region in the battery EIS Nyquist diagram can be represented by an equivalent fractional-order circuit unit, and each circuit unit contributes to the Nyquist curve shape as indicated by the arrows. It can be seen from the Figure 1 that, the EIS Nyquist diagram can be divided into four regions [24][25][26]:

‚
Region A can be represented by a single Ohmic resistance R o , this area corresponds to battery Ohmic losses, including the resistances of current collectors, active material, electrolyte and separator.
‚ High frequency region B, which can be represented by a solid electrolyte interface (SEI) resistance R sei , a SEI layer capacitance C sei , and a fractional element δ α1 , corresponds to the lithium-ion migration through the SEI film layer.

‚
Medium frequency region C, which can be represented by a charge transfer resistance R ct , a double layer capacitance C dl , and a fractional element δ α2 , corresponds to the activation kinetics and double layer effects at the interface of electrolytes; ‚ Low frequency region D, which can be represented by a Warburg element Z w and a fractional element δ α3 , corresponds to lithium-ion diffusion processes in the active material of the electrodes.
Based on the analysis above, each region in the battery EIS Nyquist diagram can be represented by an equivalent fractional-order circuit unit, and each circuit unit contributes to the Nyquist curve shape as indicated by the arrows.

Mathematical Description of Fractional-Order Model
The open circuit voltage U ocv of the battery, which is a function of battery state-of-charge (SOCq [14,15], is given by: where K 0 , K 1 , K 2 , M are coefficients to be identified, and s represents the sign of current described as follow: where ε is a small positive number.
The SOC represented as the ratio of the remaining capacity to the fully charged nominal capacity of the battery, is given by: where SOCp0q represents the initial state of charge, C n represents the battery cell nominal capacity in Ampere hour, η is the coulomb efficiency, I L is the battery current (positive for discharging, negative for charging). Thus, from Figure 1, the corresponding fractional-order differential equation for each equivalent electrical component can be expressed by: where D α denotes α order of fractional derivative (0 ď α ď 1), U sei , U ct and U w represents the voltage across equivalent resistance R sei , R ct and Warburg element Z w respectively. α 1 , α 2 and α 3 are values of fractional derivative, W d is the coefficient of Warburg element, and U T is the terminal voltage. Thus the relationship between U T , U ocv can be expressed by: where U F is the voltage across the components including R o , R sei and R ct . Thus, the total battery electrochemical impedance can be expressed as follow: From Equation (9), the fractional-order transfer function in Laplace domain can be given by Equation (10). The integer-order approximation approach of fractional calculus (3) can be then applied to the obtained fractional-order transfer function (10).

Model Parameters Identification Method
As a commonly used strategy of parameter identification, the recursive least squares method is well-suited for general ECM [14]. For the proposed fractional-order model identification, genetic Energies 2016, 9, 123 7 of 26 algorithm (GA) is particularly useful for such multi-parametric and nonlinear system estimation [27].
The purpose of GA is to find optimal solution for objective function f´χ g i¯a s follows: , .
- (11) whereχ g t is the estimated parameter values at time t for the generation g, Y t is the measured output value at time t, andŶ´χ g t¯i s the predicted output value. The parameters identification framework design for the proposed model is shown in Figure 2.

Model Parameters Identification Method
As a commonly used strategy of parameter identification, the recursive least squares method is well-suited for general ECM [14]. For the proposed fractional-order model identification, genetic algorithm (GA) is particularly useful for such multi-parametric and nonlinear system estimation [27]. The purpose of GA is to find optimal solution for objective function ̂ as follows: where ̂ is the estimated parameter values at time for the generation , is the measured output value at time , and ̂ is the predicted output value The parameters identification framework design for the proposed model is shown in Figure 2.
The detailed steps of Algorithm (1)-(4) are expressed as follows:  Initialize the original generation parameters θ within the given parameters variation range, which are achieved using recursive least squares [28]. Then the lumped parameters μ , μ , μ , μ , μ , μ , μ and ν ν , ν , ν , ν , ν , ν are generated for the fractional-order transfer function (10).  Fractional order value α α , α , α is used for calculating the fractional derivative using improved Oustaloup approximation Algorithm (3), where 9 and 10 [23].  The detailed expression of fractional transfer function can be obtained by μ , ν and s . Meanwhile, battery current is defined as the input of fractional transfer function, to calculate the output voltage of the transfer function.  The predicted terminal voltage ̂ can be obtained using (8), then the fitness value can be obtained from estimation error ̂ . Based on the fitness value, GA will take a series of actions, including elitism selection, crossover and mutation, and output the new individual parameters for next estimation iteration.

Experimental Validation of the Fractional-Order Model
The fractional-order battery modelling technique and the GA-based identification method are applied to real applications, and the experimental results and analyses are performed in this section. The parameters need to be estimated are given by: The detailed steps of Algorithm (1)-(4) are expressed as follows: ‚ Initialize the original generation parameters θ within the given parameters variation range, which are achieved using recursive least squares [28]. Then the lumped parameters µ " rµ 0 , µ 1 , µ 2 , µ 3 , µ 4 , µ 5 , µ 6 s and ν " rν 0 , ν 1 , ν 2 , ν 3 , ν 4 , ν 5 s are generated for the fractional-order transfer function (10).

‚
The detailed expression of fractional transfer function can be obtained by µ, ν and s α . Meanwhile, battery current I L is defined as the input of fractional transfer function, to calculate the output voltage U F of the transfer function.

‚
The predicted terminal voltageŶ´χ g t¯c an be obtained using (8), then the fitness value can be obtained from estimation error Y t´Ŷ´χ g t¯. Based on the fitness value, GA will take a series of actions, including elitism selection, crossover and mutation, and output the new individual parameters for next estimation iteration.

Experimental Validation of the Fractional-Order Model
The fractional-order battery modelling technique and the GA-based identification method are applied to real applications, and the experimental results and analyses are performed in this section.

Experiment Setup
In order to experimentally validate the developed fraction-order lithium-ion battery model, two types of lithium-ion cells were used in our experimental tests. The first cell is a rare earth lithium yttrium battery (WB-LYP40AHA) manufactured by the Winston Company (Shenzhen, China), the second one is a lithium-ion phosphate battery (LiFePO 4 ) manufactured by the Tenergy Company (Fremont, CA, USA). The experimental platform is developed as shown in Figure 3, and the main parameters of the two batteries are listed in Table 1. In the experimental setup, the battery cell is only connected to a programmable DC power supply and an electronic load. The power supply and electronic load are used to perform the battery charging and discharging cycles at the desired rates, respectively. Both charging and discharging cycles signals are generated by a PC using the LabVIEW software (National Instrument, Austin, TX, USA) and sent to power components via digital interfaces (GPIB and RS232). All the battery terminal voltage and current measurements from the different sensors were centralized, treated and stored in the supervising PC.

Experiment Setup
In order to experimentally validate the developed fraction-order lithium-ion battery model, two types of lithium-ion cells were used in our experimental tests. The first cell is a rare earth lithium yttrium battery (WB-LYP40AHA) manufactured by the Winston Company (Shenzhen, China), the second one is a lithium-ion phosphate battery (LiFePO4) manufactured by the Tenergy Company (Fremont, CA, USA). The experimental platform is developed as shown in Figure 3, and the main parameters of the two batteries are listed in Table 1. In the experimental setup, the battery cell is only connected to a programmable DC power supply and an electronic load. The power supply and electronic load are used to perform the battery charging and discharging cycles at the desired rates, respectively. Both charging and discharging cycles signals are generated by a PC using the LabVIEW software (National Instrument, Austin, TX, USA) and sent to power components via digital interfaces (GPIB and RS232). All the battery terminal voltage and current measurements from the different sensors were centralized, treated and stored in the supervising PC.

Battery
Nominal capacity (Ah) In order to verify the validity of the proposed model and to evaluate its performance, experiments are conducted with a hybrid pulse test and a dynamic stress test (DST) [29]. As shown in Figure 4, the hybrid pulse test is a sequence of pulse cycles including standard Hybrid Pulse Power  In order to verify the validity of the proposed model and to evaluate its performance, experiments are conducted with a hybrid pulse test and a dynamic stress test (DST) [29]. As shown in Figure 4, the hybrid pulse test is a sequence of pulse cycles including standard Hybrid Pulse Power Characterization (HPPC) profile and a custom-defined discharging/charging pulse profile [29]. The purpose of hybrid pulse test is to provide battery dynamic datasets for all SOC range for parameters identification. On Characterization (HPPC) profile and a custom-defined discharging/charging pulse profile [29]. The purpose of hybrid pulse test is to provide battery dynamic datasets for all SOC range for parameters identification. On the other hand, DST simulates the driving cycle-based conditions, and aims to evaluate the proposed model performance as shown in Figure 5.  Characterization (HPPC) profile and a custom-defined discharging/charging pulse profile [29]. The purpose of hybrid pulse test is to provide battery dynamic datasets for all SOC range for parameters identification. On the other hand, DST simulates the driving cycle-based conditions, and aims to evaluate the proposed model performance as shown in Figure 5.

Identification Results
The parameters estimation process ( Figure 2) have been performed with the recorded datasets of the hybrid pulse test (Figure 4). A number of generations 200 is selected to achieve the optimal solution for objective function ̂ . The pre-computed coefficients , , , are used in Equation (4) based on a previous work presented in [14,15]. The identified model parameters are given in Table 2. In order to validate the identification results with the proposed model, the measured voltages of hybrid pulse test on two type batteries are compared with the fractional-order model output voltages in Figure 6.

Identification Results
The parameters estimation process ( Figure 2) have been performed with the recorded datasets of the hybrid pulse test (Figure 4). A number of generations g " 200 is selected to achieve the optimal solution for objective function f´χ g i¯. The pre-computed coefficients K 0 , K 1 , K 2 , M are used in Equation (4) based on a previous work presented in [14,15]. The identified model parameters are given in Table 2. In order to validate the identification results with the proposed model, the measured voltages of hybrid pulse test on two type batteries are compared with the fractional-order model output voltages in Figure 6.

Identification Results
The parameters estimation process ( Figure 2) have been performed with the recorded datasets of the hybrid pulse test (Figure 4)

. A number of generations
200 is selected to achieve the optimal solution for objective function ̂ . The pre-computed coefficients , , , are used in Equation (4) based on a previous work presented in [14,15]. The identified model parameters are given in Table 2. In order to validate the identification results with the proposed model, the measured voltages of hybrid pulse test on two type batteries are compared with the fractional-order model output voltages in Figure 6.  From Figure 6, it can be seen that the model prediction is very close to the experimental measurement over the entire SOC range.

Accuracy of the Proposed Fractional-Order Model
To verify the accuracy of proposed model and evaluate its performance compared to other conventional modeling approaches, a classical equivalent circuit second-order RC (2-RC) model [30][31][32][33] is adopted here to use as the reference.
As shown in the Figure  From the zoom-in comparisons shown in Figures 8 and 9, it can be concluded that the proposed fractional-order model results show a very good accuracy compared with the experimental measurements over the entire SOC range. More specifically, the voltage dynamics of charging and discharging are correctly and precisely reproduced by the fractional-order model, which cannot be accurately predicted by the classical 2-RC model. From Figure 6, it can be seen that the model prediction is very close to the experimental measurement over the entire SOC range.

Accuracy of the Proposed Fractional-Order Model
To verify the accuracy of proposed model and evaluate its performance compared to other conventional modeling approaches, a classical equivalent circuit second-order RC (2-RC) model [30][31][32][33] is adopted here to use as the reference.
As shown in the Figure 7, the battery terminal voltage V can be estimated based on the model applied current I and battery state of charge (SOC), which can be calculated using the Coulomb counting method. The open circuit voltage V OC is a function of SOC. The series resistance R s and two parallel RC circuit including R t,s , R t,l , C t,l , and C t,l describe how the battery terminal voltage V dynamically changes under a variable applied current I. From Figure 6, it can be seen that the model prediction is very close to the experimental measurement over the entire SOC range.

Accuracy of the Proposed Fractional-Order Model
To verify the accuracy of proposed model and evaluate its performance compared to other conventional modeling approaches, a classical equivalent circuit second-order RC (2-RC) model [30][31][32][33] is adopted here to use as the reference.
As shown in the Figure  From the zoom-in comparisons shown in Figures 8 and 9, it can be concluded that the proposed fractional-order model results show a very good accuracy compared with the experimental measurements over the entire SOC range. More specifically, the voltage dynamics of charging and discharging are correctly and precisely reproduced by the fractional-order model, which cannot be accurately predicted by the classical 2-RC model. From the zoom-in comparisons shown in Figures 8 and 9 it can be concluded that the proposed fractional-order model results show a very good accuracy compared with the experimental measurements over the entire SOC range. More specifically, the voltage dynamics of charging and discharging are correctly and precisely reproduced by the fractional-order model, which cannot be accurately predicted by the classical 2-RC model.   To further show the performance and advantages of the proposed fractional-order model, the 2-RC model output voltages and fractional-order model output voltages are compared with measured voltages for the one single cycle in DST, as shown in Figure 10. To further show the performance and advantages of the proposed fractional-order model, the 2-RC model output voltages and fractional-order model output voltages are compared with measured voltages for the one single cycle in DST, as shown in Figure 10. In addition, the model prediction error numerical results during the entire test including rootmean-square error (RMSE) and maximum estimation error (MAE) are shown in Table 3. It can be seen that the proposed fractional-order model can achieve better accuracy compared to the 2-RC model.

Model Parameters Sensitivity Analysis
The fractional-order battery model can achieve a better dynamic modelling accuracy compared to the integer-order model, as shown and discussed in the previous section. However, it should be noted that each fractional-order parameter has its own influence with a different degree of sensitivity on the model output accuracy. Thus a parameter sensitivity study for a fractional-order model represents an important step for the model development, because such a study can indicate the most and least influencing parameter and leads to a better understanding of the model output dynamic, In addition, the model prediction error numerical results during the entire test including root-mean-square error (RMSE) and maximum estimation error (MAE) are shown in Table 3. It can be seen that the proposed fractional-order model can achieve better accuracy compared to the 2-RC model.

Model Parameters Sensitivity Analysis
The fractional-order battery model can achieve a better dynamic modelling accuracy compared to the integer-order model, as shown and discussed in the previous section. However, it should be noted that each fractional-order parameter has its own influence with a different degree of sensitivity on the model output accuracy. Thus a parameter sensitivity study for a fractional-order model represents an important step for the model development, because such a study can indicate the most and least influencing parameter and leads to a better understanding of the model output dynamic, which is mainly affected by fractional-order parameters. In this section, a modelling parameter sensitivity analysis is presented and discussed in detail, based on the previously developed fractional-order lithium-ion battery model.

Sensitivity Analysis Method
To analyze the influence of different modelling parameters on the fractional-order model output voltage, the parameters sensitivity study is performed and discussed using statistical Multi-Parametric Sensitivity Analyze (MPSA) method [34,35]. The different steps of the MPSA method can be expressed as follows:

‚
Select the fractional-order model parameters need to be analyzed, set an appropriate variation range (numerical boundary) of each parameter.
‚ Generate a series of 700 independent random numbers with a uniform distribution within in the selected parameter design range.
‚ Run the fractional-order model using selected typical parameter values in order to get the model typical output voltage V typical,i , where i is the operating current (positive for discharging, negative for charging); Run the fractional-order model using generated sets of 700 independent numbers for one parameter to get the model distribution output voltages V distribution,i,n , where n is the generated numbers. Then the relative sensitivity criteria at different operating current values i of the selected parameter can be computed using the following equation: ‚ In order to calculate the overall parameter sensitivity over the whole battery operating current range, the parameter sensitivity index value Si can be defined as the sum of the relative sensitivity criteria at different battery operating current values:

MPSA Method Results and Discussion
In the proposed fractional-order model, eight individual parameters have been selected for the sensitivity analysis. These parameters are: circuit component coefficients (R sei , R ct , C sei , C dl , W d ) and fractional-order parameters pα 1 , α 2 , α 3 q. The typical values of these parameters have been chosen based on the GA-identified results, and variation range is set to˘20% of their typical value. By applying the previously described MPSA method, the sensitivity analysis results of fractional-order output voltage for each parameter over the experimental battery operating current rate range (from´1C to 1C) are summarized in Table 4.
It can be clearly seen from the Table 4 that, different parameters have different influence to the fractional-order model characteristics. Generally, they can be divided into three categories: highly sensitive parameters R sei , α 2 (the index value is greater than 1.0ˆ10´1), sensitive parameters C dl , α 1 , C sei (the index value ranges from 1.0ˆ10´4 to 1.0ˆ10´1) and insensitive parameters R ct , α 3 , W d (the index value is less than 1.0ˆ10´4). The results indicate that, in the process of GA identification, the fractional-order model error between the output voltage and the experimentation could be effectively corrected by changing the values of highly sensitive and sensitive parameters. Thus, in order to improve the accuracy of fractional-order model and the efficiency of the identification process, the circuit element parameters pR sei , C sei , C dl q and fractional-order parameters (α 1 , α 2 q should be carefully chosen. It should also be noted that, each parameters correspond different electrochemical dynamic behaviors as mentioned in Section 2.2. The parameter sensitivity analysis results further show that the lithium-ion migration through the SEI film layer (associated parameter R sei ), and the double layer effects at the interface of electrolytes (associated parameters C dl , α 2 ) have an high impact on the battery terminal voltage. The capacitance of the SEI layer (associated parameters C sei , α 1 ) have an medium impact on the battery terminal voltage. The activation kinetics in both negative and positive electrodes (associated parameter R ct ), and the lithium-ion diffusion processes in the active material of the electrodes (associated parameters W d , α 3 ) have an low impact on the battery terminal voltage. This meaningful information reveals that the battery internal dynamics have different degree influences on battery performance.
In addition to the analysis of index value for each parameter, the individual relative sensitivity criteria values of three fractional-order parameters at different operating current values are also presented. The relative sensitivity index value of fractional-order factors are presented in Figure 11. It can be seen from the Figure 11 that, fractional-order factor α 1 influences the battery output voltage over the entire current rate range, especially around˘0.4C point and at high current rate (˘1C).  The results indicate that, in the process of GA identification, the fractional-order model error between the output voltage and the experimentation could be effectively corrected by changing the values of highly sensitive and sensitive parameters. Thus, in order to improve the accuracy of fractional-order model and the efficiency of the identification process, the circuit element parameters , , and fractional-order parameters (α , α should be carefully chosen. It should also be noted that, each parameters correspond different electrochemical dynamic behaviors as mentioned in Section 2.2. The parameter sensitivity analysis results further show that the lithium-ion migration through the SEI film layer (associated parameter ), and the double layer effects at the interface of electrolytes (associated parameters , α ) have an high impact on the battery terminal voltage. The capacitance of the SEI layer (associated parameters , α ) have an medium impact on the battery terminal voltage. The activation kinetics in both negative and positive electrodes (associated parameter ), and the lithium-ion diffusion processes in the active material of the electrodes (associated parameters , α ) have an low impact on the battery terminal voltage. This meaningful information reveals that the battery internal dynamics have different degree influences on battery performance.
In addition to the analysis of index value for each parameter, the individual relative sensitivity criteria values of three fractional-order parameters at different operating current values are also presented. The relative sensitivity index value of fractional-order factors are presented in Figure 11. It can be seen from the Figure 11 that, fractional-order factor α influences the battery output voltage over the entire current rate range, especially around 0.4C point and at high current rate ( 1C). The sensitivity criteria analysis results of fractional-order factors α and α are given in Figure 12 and Figure 13, respectively. It can be seen from the two figures that, the highly sensitive factor α is the most sensitive parameter among the three fractional-order parameters. Furthermore, its value has a higher impact on the model output voltage at high current rate range. α is an insensitive parameter to the fractional-order model, an increase of the current rate increases the relative sensitivity of the parameter α also. The sensitivity criteria analysis results of fractional-order factors α 2 and α 3 are given in Figures 12  and 13 respectively. It can be seen from the two figures that, the highly sensitive factor α 2 is the most sensitive parameter among the three fractional-order parameters. Furthermore, its value has a higher  In order to give a comprehensive analysis of index value for each fractional-order parameter over the entire state of charge (SOC) range, the individual relative sensitivity criteria values of three fractional-order parameters at constant 1C to 5C discharge current are presented in Figures 14-16.   In order to give a comprehensive analysis of index value for each fractional-order parameter over the entire state of charge (SOC) range, the individual relative sensitivity criteria values of three fractional-order parameters at constant 1C to 5C discharge current are presented in Figures 14-16.  In order to give a comprehensive analysis of index value for each fractional-order parameter over the entire state of charge (SOC) range, the individual relative sensitivity criteria values of three fractional-order parameters at constant 1C to 5C discharge current are presented in Figures 14-16.  In order to give a comprehensive analysis of index value for each fractional-order parameter over the entire state of charge (SOC) range, the individual relative sensitivity criteria values of three fractional-order parameters at constant 1C to 5C discharge current are presented in Figures 14-16.    It can be seen again from the three figures that the highly sensitive factor α is the most sensitive parameter among the three fractional-order parameters. All the three fractional-order parameters relative sensitivity increase with the current rates. α has the lowest impact on the model output voltage at 75% SOC. A decrease of the SOC increases the relative sensitivity of the parameters α and α .
The above sensitivity analysis results of these parameters can thus provide useful information during the parameter fitting process of the model.

Dynamic Effect of Fractional-Order Parameters
According to the analysis results in the previous section, the proposed model output voltage is more sensitive to fractional-order parameters at medium and high current rate. For this reason, one single cycle in the hybrid pulse test (from figure 4a), which includes various discharging /charging current rate pulses, is selected as the input for the proposed model to further investigate the influence of fractional-order parameters (α , α , α , on the model output dynamic behavior, as shown in Figure 17.    It can be seen again from the three figures that the highly sensitive factor α is the most sensitive parameter among the three fractional-order parameters. All the three fractional-order parameters relative sensitivity increase with the current rates. α has the lowest impact on the model output voltage at 75% SOC. A decrease of the SOC increases the relative sensitivity of the parameters α and α .
The above sensitivity analysis results of these parameters can thus provide useful information during the parameter fitting process of the model.

Dynamic Effect of Fractional-Order Parameters
According to the analysis results in the previous section, the proposed model output voltage is more sensitive to fractional-order parameters at medium and high current rate. For this reason, one single cycle in the hybrid pulse test (from figure 4a), which includes various discharging /charging current rate pulses, is selected as the input for the proposed model to further investigate the influence of fractional-order parameters (α , α , α , on the model output dynamic behavior, as shown in Figure 17.   It can be seen again from the three figures that the highly sensitive factor α 2 is the most sensitive parameter among the three fractional-order parameters. All the three fractional-order parameters relative sensitivity increase with the current rates. α 1 has the lowest impact on the model output voltage at 75% SOC. A decrease of the SOC increases the relative sensitivity of the parameters α 2 and α 3 .
The above sensitivity analysis results of these parameters can thus provide useful information during the parameter fitting process of the model.

Dynamic Effect of Fractional-Order Parameters
According to the analysis results in the previous section, the proposed model output voltage is more sensitive to fractional-order parameters at medium and high current rate. For this reason, one single cycle in the hybrid pulse test (from Figure 4a), which includes various discharging /charging current rate pulses, is selected as the input for the proposed model to further investigate the influence of fractional-order parameters (α 1 , α 2 , α 3 q, on the model output dynamic behavior, as shown in Figure 17. In this study, each fractional-order parameter changes its value from fractional values (0.1-0.9) to integer value of 1. The corresponding fractional-order model output voltages are then presented in Figure 18-20.
From Figure 18, it is easy to observe that the sensitive fractional-order parameter α have an important influence on the dynamic behavior of model output. Compared with classic integer-order models, the fractional-order model can effectively better describe the battery dynamic voltage responses with additional modeling degree of freedom. The dynamic effect of the fractional-order for highly sensitive fractional-order parameter α is shown in Figure 19. The variation of α has indeed a stronger influence than α , especially at high current. These results are in agreement with the previous analysis presented in Section 4.2.
In addition, it can be clearly seen from the Figure 19 that, as the parameter α increase from 0.1 to 1, the voltage dynamic transient time [36] gradually increases. In this study, each fractional-order parameter changes its value from fractional values (0.1-0.9) to integer value of 1. The corresponding fractional-order model output voltages are then presented in  From Figure 18, it is easy to observe that the sensitive fractional-order parameter α 1 have an important influence on the dynamic behavior of model output. Compared with classic integer-order models, the fractional-order model can effectively better describe the battery dynamic voltage responses with additional modeling degree of freedom. In this study, each fractional-order parameter changes its value from fractional values (0.1-0.9) to integer value of 1. The corresponding fractional-order model output voltages are then presented in Figure 18-20.
From Figure 18, it is easy to observe that the sensitive fractional-order parameter α have an important influence on the dynamic behavior of model output. Compared with classic integer-order models, the fractional-order model can effectively better describe the battery dynamic voltage responses with additional modeling degree of freedom. The dynamic effect of the fractional-order for highly sensitive fractional-order parameter α is shown in Figure 19. The variation of α has indeed a stronger influence than α , especially at high current. These results are in agreement with the previous analysis presented in Section 4.2. In addition, it can be clearly seen from the Figure 19 that, as the parameter α increase from 0.1 to 1, the voltage dynamic transient time [36] gradually increases. The dynamic effect of the fractional-order for highly sensitive fractional-order parameter α 2 is shown in Figure 19. The variation of α 2 has indeed a stronger influence than α 1 , especially at high current. These results are in agreement with the previous analysis presented in Section 4.2. In addition, it can be clearly seen from the Figure 19 that, as the parameter α 2 increase from 0.1 to 1, the voltage dynamic transient time [36] gradually increases.  Figure 20 shows the different results of fractional-order model output voltage by changing value of parameter α . In order to give a clearer show, the zoom-in at about 371 s and 1581 s for these ten curves are given. It can be seen clearly from figure that, the insensitive parameter α has very small impact on the model output voltage. This results are also in agreement with the previous analysis presented in Section 4.2. Summarizing the above analysis, it can be concluded that, the fractional-order parameters are more objective and original to reveal the nature of the dynamic physical phenomena.

Conclusions
In this paper, a novel fractional-order of lithium-ion battery has been developed in the first part. The fractional-order modeling approach, as well as the mathematical modeling equations, are introduced and explained based on an improved Oustaloup fractional-order approximation method. In this case, all the internal lithium-ion battery dynamic behaviors can be accurately captured by the proposed free fractional-orders elements. Compared to the classic integer-order equivalent circuit   Figure 20 shows the different results of fractional-order model output voltage by changing value of parameter α 3 . In order to give a clearer show, the zoom-in at about 371 s and 1581 s for these ten curves are given. It can be seen clearly from figure that, the insensitive parameter α 3 has very small impact on the model output voltage. This results are also in agreement with the previous analysis presented in Section 4.2.  Figure 20 shows the different results of fractional-order model output voltage by changing value of parameter α . In order to give a clearer show, the zoom-in at about 371 s and 1581 s for these ten curves are given. It can be seen clearly from figure that, the insensitive parameter α has very small impact on the model output voltage. This results are also in agreement with the previous analysis presented in Section 4.2. Summarizing the above analysis, it can be concluded that, the fractional-order parameters are more objective and original to reveal the nature of the dynamic physical phenomena.

Conclusions
In this paper, a novel fractional-order of lithium-ion battery has been developed in the first part. The fractional-order modeling approach, as well as the mathematical modeling equations, are introduced and explained based on an improved Oustaloup fractional-order approximation method. In this case, all the internal lithium-ion battery dynamic behaviors can be accurately captured by the proposed free fractional-orders elements. Compared to the classic integer-order equivalent circuit  Summarizing the above analysis, it can be concluded that, the fractional-order parameters are more objective and original to reveal the nature of the dynamic physical phenomena.

Conclusions
In this paper, a novel fractional-order of lithium-ion battery has been developed in the first part. The fractional-order modeling approach, as well as the mathematical modeling equations, are introduced and explained based on an improved Oustaloup fractional-order approximation method. In this case, all the internal lithium-ion battery dynamic behaviors can be accurately captured by the proposed free fractional-orders elements. Compared to the classic integer-order equivalent circuit models in the literature, a fractional-order model has the advantage of higher modeling accuracy especially for the battery dynamic voltage behavior under transient operation conditions. The parameter identification of the proposed model is conducted using genetic algorithm, and the model has been validated experimentally with two different types of lithium-ion batteries. In order to evaluate its performance, one classical equivalent circuit second-order RC model is used for comparison purpose. The results show that the proposed model can achieve better accuracy over the entire battery SOC range, and the dynamic of charging/discharging voltages can be accurately reproduced by the fractional-order model.
Although parameter determination and identification are very important in the lithium-ion battery model development, especially for the fractional-order modelling which includes additional fractional-order parameters, very little information has been published so far in the literature on the analysis of sensitivity and dynamic effect of fractional-order parameters. This paper investigates in particular in the second part the sensitivities and dynamic effects of the parameters in the proposed fractional-order model. The parameter sensitivity analysis are performed using a statistical multi-parametric sensitivity analysis method. The effect of each parameter value deviation (˘20%) on the output voltages at different operating current values are presented and discussed. From the analysis results, the modelling parameters are divided into three categories based on their sensitivity index value: highly sensitive parameters (R 1 and α 2 ), sensitive parameters (C dl,1 , C dl,2 , α 1 ) and insensitive parameters (R 2 , α 3 , W d ). These analysis results also reveal that the internal battery dynamics have different degrees of influence on battery performance. In addition, the influence of three fractional-order parameters value (α 1 , α 2 , α 3 ) in proposed model on the output voltage are also shown and discussed. The analysis results can provide important modelling information for dynamic model development and parameter identification.