Next Article in Journal
Decarbonizing the Chilean Electric Power System: A Prospective Analysis of Alternative Carbon Emissions Policies
Previous Article in Journal
Assessment on the Application of Facilitated Transport Membranes in Cement Plants for CO2 Capture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Fragile H Nonlinear Observer for State of Charge Estimation of Lithium-Ion Battery Based on a Fractional-Order Model

1
Institute of Energy Storage and New Materials Technology, Dongfang Electric Corporation Science and Technology Research Institute Co., Ltd., 18 Xixin Avenue High-Tech Zone West Park, Chengdu 611731, China
2
School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Energies 2021, 14(16), 4771; https://doi.org/10.3390/en14164771
Submission received: 15 June 2021 / Revised: 2 July 2021 / Accepted: 21 July 2021 / Published: 5 August 2021

Abstract

:
This paper deals with the state of charge (SOC) estimation of lithium-ion battery (LIB) in electric vehicles (EVs). In order to accurately describe the dynamic behavior of the battery, a fractional 2nd-order RC model of the battery pack is established. The factional-order battery state equations are characterized by the continuous frequency distributed model. Then, in order to ensure the effective function of nonlinear function, Lipschitz condition and unilateral Lipschitz condition are proposed to solve the problem of nonlinear output equation in the process of observer design. Next, the linear matrix (LMIS) inequality based on Lyapunov’s stability theory and H method is presented as a description of the design criteria for non-fragile observer. Compared with the existing literature that adopts observers, the proposed method takes the advantages of fractional-order systems in modeling accuracy, the robustness of H method in restricting the unknown variables, and the non-fragile property for tolerating slow drifts on observer gain. Finally, The LiCoO 2 LIB module is utilized to verify the effectiveness of the proposed observer method in different operation conditions. Experimental results show that the maximum estimation accuracy of the proposed non-fragile observer under three different dynamic conditions is less than 2%.

1. Introduction

Due to the huge fossil energy consumption since the industrial revolution, electric vehicles (EVs) have been greatly developed [1]. The key energy source of electric vehicles, lithium-ion batteries (LIBs), is considered to be an important way of energy storage because of its high energy storage per unit volume, stable performance, and the potential for lower production costs [2].The state of charge (SOC) of LIB, defined as the ratio of the remaining capacity to the nominal capacity of the battery is crucial to prevent battery from over-charging or over-discharging, fire, and even explosion. Meanwhile, SOC stands for the most basic parameter in a battery management system (BMS) in EVs. The electrochemistry complexity of the battery’s internal reaction usually causes the high nonlinearities of battery systems. Moreover, in practice, SOC cannot be measured directly by sensors.
In the existing studies, a series of equivalent circuit models (ECMs) is carried out to capture nonlinearities of LIB comprehensively, such as Rint Model [3], Thevenin model [4], 2nd-order RC (resistor-capacitor) model [5], etc. More recently, many researchers have found that nonlinear systems have fractional properties [6]. The fractional-order (FO) battery models, which consist of fractional elements, including Warburg element [7], constant phase element (CPE) [7,8], and fractional capacitor [9,10,11,12], have been introduced to model the nonlinear characteristics of battery. The three electrical components, indeed, can be unified by the same equation. Compared with integer-order (IO) models, fractional-order models can model nonlinearities of battery more accurately to achieve better approximation of the battery system [7,8,11,12]. Among them, the Kalman filter-based approaches [7,9] often suppose that the model and measurement noises are Gaussian, and assume that both the covariance matrix of measurement noise and model error are known. On the other hand, the observer-based approaches [8,11,12] do not demand knowledge on noise distributions which are considered to be more practical, of which the observer in [8] with constant gain was slightly weak in tracking the nonlinear dynamic process. In fact, only for integer-order battery models, the observers with dynamic gain were proposed to adapt the nonlinear dynamics [13,14]. Moreover, in order to improve the role of the observer to suppress the measurement noise and modeling error, the H observers for integer-order battery model were introduced in [15,16]. Finally, note that the fractional-order calculus, an extension of the classical integer-order calculus, is developed in various fields such as control engineering [17,18], signal processing [19,20], and system modeling, including supercapacitors [21] and LIBs [7].
In practical applications, many SOC estimation methods are implemented digitally on microcontrollers [22], and the FOC is also actualized by numerical techniques [7,8,9,10,11,12]. Meanwhile, some perturbations such as external disturbances and slow drifts on the observer gain may cause serious deterioration of the observer performance [20,23,24]. This brings up a fragility problem that has been of great interest in the control theory area [20,24,25,26,27]. As the observer gains are usually obtained from LMIs, which are calculated off-line [8,13,14,15,16], a non-fragile observer that tolerates some perturbations or gain fluctuations [25] and allows for on-line tuning of model parameters to retrench the cost of implementation is essential. The authors of [26] focused on a class of nonlinear fractional-order uncertain systems with bounded perturbation on the observer gain. In [20], the authors dealt with the problem of non-fragile observer design for a class of Lipschitz nonlinear fractional-order systems. In [24,25,27], non-fragile observers based on H performance have been investigated.
For the purpose of capturing nonlinearities of LIB in modeling and considering the advantages of the H observer and its practical applications, a non-fragile H observer method with dynamic gain based on fractional 2nd-order RC model is proposed in this paper. Compared with the existing literature, the main contributions of this paper are as follows:
  • Kalman filter-based approaches [7,9] usually assume that the model and measurement noises are Gaussian, which is not the case in actual situations. In order to restrict the effects of the non-Gaussian model and measurement noises [28] for fractional-order battery models, a H -based nonlinear observer is proposed. Furthermore, the proposed observer with dynamic gain is also considered to outperform the constant gain observer [8]. The dynamic gain means that the proposed observer can estimate the SOC accurately in various dynamic operation conditions.
  • The fragility problem involves a trade-off between implementation accuracy and performance deterioration of the observer implementation. It has caused extensive concern in [20,24,25,26,27,29], but it has not be taken seriously in research areas of SOC estimation. Therefore, to further improve the accuracy of SOC estimation of observer-based methods, we put forward the non-fragile observer. For SOC estimation of fractional-order battery model, the non-fragile observer design criterion for the gain of proposed observer is crucial to solve the fragility problem in practical and real-time applications. Considering the additive perturbation on observer gain with known bound make it implement more properly than the off-line calculation of the observer gain in [8]. Moreover, it is vital to improve the reliability and robustness of all observer-based methods for SOC estimation.
The content of this article is arranged as follows. In Section 2, the fractional-order 2 RC equivalent circuit model of the battery pack is explained. Section 3 consists of experiments and a least-square method mixed with physical behavior based identification method for model parameters. In Section 4, the nonlinear observer design criterion is presented to design the SOC estimator. Experimental validations of the SOC estimation are discussed in Section 5. Conclusions are given in Section 6.

2. Fractional-Order Modeling for LiCoO 2 Lithium-Ion Battery (LIB) Module

2.1. Constant Phase Element (CPE)

The electrochemical impedance spectroscopy (EIS) method is based on the measurement of the electrical voltage response of a LIB by applying a wide frequency range of a.c. signal. CPE is an equivalent electrical circuit component to describe the investigated LIB impedance, that is, an imperfect capacitor [30]. Meanwhile, it was introduced in [31] that the majority of capacitors show fractional properties in practical dynamic conditions. In consequence, the electrical impedance of the CPE could be represented as the following form [9]:
Z C P E ( j ω ) = 1 C C P E ( j ω ) α , 0 < α < 1
or
Z C P E ( s ) = 1 C C P E s α , s = j ω , 0 < α < 1
where Z C P E ( j ω ) = u C P E ( j ω ) / i C P E ( j ω ) is the electrical impedance of the CPE when the a.c. signal i C P E ( j ω ) is applied with frequency ω ; i C P E and u C P E denote the current and voltage of the CPE, respectively; C C P E is the capacitance of the CPE, j is the imaginary unit; α is the fractional-order; and s denotes the Laplace variable.
It follows that
i C P E ( t ) = C C P E d α u C P E ( t ) d t α

2.2. Fractional 2nd-Order RC Battery Model

The Grünwald–Letnikov (GL) definition for generalized fractional integro-differential system is commonly used.
Definition 1.
The α-th order GL fractional derivative of x ( t ) is defined by [9,32]
D t α x ( t ) = lim h 0 h α j = 0 [ t / h ] ω j α x ( t j h ) ,
where D t α = d α d t α , α > 0 1 , α = 0 0 t ( d τ ) α . α < 0 is the continuous integro-differential operator, the factor ω j α = ( 1 ) j α j , ω 0 α = 1 , α j is the binomial coefficient, thus ω j α = ( 1 α + 1 j ) ω j 1 α , and [t/h] is the integer part of t / h .
According to Definition 1, the following GL fractional-order integral form when α < 0 can be obtained:
D α x ( t ) = 0 t x ( τ ) ( d τ ) α = lim h 0 h α j = 0 [ t / h ] ω j α x ( t j h ) .
As is mentioned in the introduction, various kinds of ECMs have been widely used in battery modeling, most of which comprise one or more capacitors. Generally, the capacitors of the ECMs are assumed to be integer-order. As is mentioned in [31], the fractional capacitor model can reproduce and predict the capacitor’s behavior much better than any other theory. Considering the fractional properties of capacitor, in this paper, the fractional 2nd-order RC battery model is introduced by adopting two CPEs to replace the two integer-order capacitors in the commonly used integer 2nd-order RC ECM [5,33], which had been studied in [7,8]. Another reason to employ this model is that we find it can perform better in modeling the change of terminal voltage than the integer 2nd-order RC model. As such, the fractional 2nd-order RC ECM is employed here, and its schematic diagram is shown in Figure 1.
In Figure 1, the notation U OC denotes the OCV related to SOC; I T is the operating current, which is positive in the discharge process and negative in the charge process; U T indicates the terminal voltage; and R 0 is Ohmic resistance. The notations R 1 and C 1 are the electrochemical polarization resistance and constant phase element CPE 1 , respectively. R 2 and C 2 are the concentration polarization resistance and constant phase element CPE 2 , respectively. In addition, U 1 and U 2 denote the voltages of CPE 1 and CPE 2 , respectively.
Next, we will build a battery model based on Kirchhoff Voltage laws with relationship Equation (3).
State Equation:
d α 1 U 1 d t α 1 = 1 R 1 C 1 U 1 + 1 C 1 I T , d α 2 U 2 d t α 2 = 1 R 2 C 2 U 2 + 1 C 2 I T , d S O C d t = 1 Q n I T .
Output equation:
U T = U O C ( S O C ) R 0 I T U 1 U 2 ,
where Q n is the nominal capacity of the battery, and 0 < α 1 , α 2 < 1 are the fractional orders of the CPE 1 and CPE 2 , respectively.

3. Experiments and Identification of Model Parameters

In this section, we want to identify the unknown parameters α 1 , α 2 , R 0 , R 1 , R 2 , C 1 , C 2 and nonlinear function U O C ( S O C ) in battery Equations (6) and (7) by experiments. First, let us introduce our test bench.

3.1. Test Bench

The test bench is illustrated in Figure 2; two sets of devices are used to measure the battery characteristics at the same time and record the battery voltage and current data. Device 1 is NEWARE BTS-4000. Because of its high measurement accuracy and low measurement noise, the measured value of this device is used as a reference value. Meanwhile, using a set of home-made equipment to measure battery voltage and current data, the devices have measurement noise and errors during the data collection process, which can be used to simulate the real measurement environment and verify the effectiveness of the proposed SOC estimation method. Finally, the cell type used in the experiments is NCR18650BE, a LIB with cathode material lithium cobalt oxide and anode material graphite made by Panasonic Corp., Kadoma, Japan. Its nominal voltage is 3.6 V, and the nominal capacity is 3.2 Ah. Moreover, the test object is a battery pack with ten batteries connected in parallel. So the nominal capacity of the battery pack Q n equals to 32 Ah.

3.2. Parameters and Nonlinear Function Identification

Step 1. Identify the parameters R 0 , R 1 , R 2 , C 1 , C 2 , and U O C ( S O C ) . Here, both the orders α 1 and α 2 are assumed to be unity. Using the method employed in [5] with current pulses, these parameters R 0 ( m Ω ) , R 1 ( m Ω ) , R 2 ( m Ω ) , C 1 ( k F ) , and C 2 ( k F ) are identified as the blue points in Figure 3.
The mean values of the SOC interval [ 0.2 , 1 ] in Figure 3 are set to the constant identified results. The eighth-order polynomial is employed to fit the nonlinear function of OCV and SOC and it is given by
U O C ( S O C ) = p 1 S O C 8 + p 2 S O C 7 + + p 9 ,
where [ p 1 , p 2 , , p 9 ] = [ 775.25 , 3273.15 , 5691.59 , 5264.44 , 2795.32 , 860.49 , 148.66 , 14.13 , 2.78 ] . The measurement data and fitted curve of OCV vs. SOC are shown in Figure 4. As a result, the parameters are shown in Table 1.
Step 2. Identify the fractional orders α 1 and α 2 . Based on the identified parameters in Step 1, the fractional orders α 1 and α 2 are further employed to improve modeling accuracy. To this end, the least-square (LS) method is used to minimize the modeling error, that is,
min α 1 , α 2 j = 1 N U T * ( j ) U T ( j ) 2
where N is the number of the samples, U T * ( j ) is the measured terminal voltage, and U T ( j ) is the predictive voltage using the identified model. Notice that for given fractional orders α 1 and α 2 , the predictive voltage U T ( j ) is given by
U T = U O C ( S O C ) R 0 I T h α 1 j = 0 k ω j α 1 R 1 I ( t j h ) U 1 ( t j h ) R 1 C 1 h α 2 j = 0 k ω j α 2 R 2 I ( t j h ) U 2 ( t j h ) R 2 C 2 ,
where h is the sample time. The identification results of parameters α 1 and α 2 are also shown in Table 1.

4. Non-Fragile H Nonlinear Observer for SOC Estimation Based on Fractional-Order Model

4.1. Battery State-Space Model

The battery Equations (6) and (7) can be rewritten as
D α x = A x + B u , x ( 0 ) = x 0 , y = h ( x ) + C x R 0 u ,
where α = [ α 1 , α 2 , 1 ] , x = [ U 1 , U 2 , S O C ] T , y = U T , u = I T , h ( x ) = U O C ( S O C ) , x 0 is the initial state and
A = 1 τ 1 0 0 0 1 τ 2 0 0 0 0 , B = 1 C 1 1 C 2 1 Q n , C T = 1 1 0 .
As the battery charging and discharging process involves complicated physical and chemical reactions, the battery Equation (18) are further rewritten as
D α x = A x + B u + d x , x ( 0 ) = x 0 , y = h ( x ) + C x R 0 u + d y ,
where d x represents the state disturbance caused by the modeling error of current and measurement noise, and d y represents the output disturbance caused by measurement noise of current and terminal voltage. In the observer-based SOC estimation method, assuming that only d x and d y bounded, the range is d x < and d y < .

4.2. Discussions for SOC-OCV Function U O C ( S O C )

The estimated values of the state x and S O C are x ^ and S O C ^ , respectively, and h ˜ h ( x ) h ( x ^ ) . The S O C values range from [0,1] and the function of U O C ( S O C ) is a monotonically increasing. Then, we have
β min U ˙ O C ( S O C ) β max ,
with β min , β max > 0 . In addition, the battery used in the experiment during the test, according to the SOC-OCV curve Equation (8), is satisfied by the the above Equation (13) with
β min = 0.4551 , β max = 4.4708 .
In order to design a nonlinear observer for the above system Equation (12), we first need to discuss the properties of the nonlinear function h ( x ) .
There exists matrix Q, W, for any x , x ^ with S O C , S O C ^ [ 0 , 1 ] , such that the nonlinear functions h ( x ) satisfies the following one-sided Lipschitz condition Equation (15) and Lipschitz condition Equation (16):
x ˜ T h x x = x ^ T h ˜ x ˜ T Q x ˜ ,
h ˜ T h ˜ x ˜ T W 2 x ˜ ,
where Q = 0 0 0 0 0 0 0 0 β min 2 , W = 0 0 0 0 0 0 0 0 β max 2 .
Proof. 
Please see the authors’ previous work [5]. □

4.3. Non-Fragile Observer Designed for SOC Estimation

The following two lemmas are crucial for designing the observer based on fractional-order model and establish non-fragile observer design criterion, respectively.
Lemma 1.
The factional-order system [17,18,19]
D t α x ( t ) = f ( x ( t ) ) ,
can be characterized by the continuous frequency distributed model
Z x ( ω , t ) t = ω Z x ( ω , t ) + f ( x ( t ) ) , x ( t ) = 0 μ ( ω ) Z x ( ω , t ) d ω ,
where ω is the elementary frequency, z x ( ω , t ) is the infinite dimension distributed state variable, and μ ( ω ) = s i n ( α π ) π ω α .
Lemma 2.
Ref. [20] Let X R n , Y R n , and ε > 0 . Then we have X T Y + Y T X ε Y T Y + 1 ε X T X .
Considering the application of the above mentioned observer-based methods [5,8,11,12,13,14], the non-fragile observer design criterion is introduced here. These observers would be implemented digitally on microcontrollers as many other SOC estimation methods [22]. Meanwhile, as is mentioned in [8], the real application of FOC should be actualized on microcontrollers by numerical techniques. However, the electrical signal loss or interference on the D/A converter (DAC) and Hall sensor is unavoidable in the signals’ conversion and transmission process. Besides, there exist rounding error when the microcontroller unit (MCU) processes floating point numbers, thermal noise of electric conductor, stray electromagnetic field interference, etc. All these bring perturbations to the observer gain lead to a fragility problem, and the fragility problem indeed requires a trade-off of observer implementation accuracy and performance deterioration [23]. Thus, if the fragility problem is not considered, the estimation accuracy of observer-based methods will not be ideal enough as it is designed by H performance. The non-fragile observer provides a optimization algorithm to improve practicability of SOC estimation.
Therefore, we further develop the nonlinear observer based on integer-order model in [5,13,14] to a non-fragile observer based on fractional-order model that tolerates some perturbations as follows:
D α x ^ = A x ^ + B u + ( L + Δ L ) H ( x ^ ) ( y y ^ ) , x ^ ( 0 ) = x ^ 0 , y ^ = h ( x ^ ) + C x ^ R 0 u ,
where y ^ is the estimate of the real terminal voltage y, L is a matrix that will be designed later, but the vector ( L + Δ L ) h x x = x ^ T is the actual observer gain that is dynamic. Let h x x = x ^ T H ( x ^ ) for simplification. Δ L H ( x ^ ) is additive perturbation on the error gain with the bound satisfying Δ L δ and δ is a positive constant. Note that the upper bound δ will be calculated later. It follows that the error dynamics is given by
D α x ˜ = ( A ( L + Δ L ) H ( x ^ ) C ) x ˜ ( L + Δ L ) H ( x ^ ) h ˜ + ( E ( L + Δ L ) H ( x ^ ) F ) d
where d = d x d y is the synthetic disturbance, E = [ I 0 , 0 ] , F = [ 0 , 1 ] , and I 0 denotes the appropriate dimensional identity matrix.
The main contribution of this paper, the non-fragile H design criteria with H performance, is then established as shown in the Figure 5.
Theorem 1.
Regarding to the non-fragile observer Equation (19), it has a stable observation under system disturbance Equation (12), for given positive scalars ϵ < 2 β min 2 , and γ, if there exist positive real number ε 1 , ε 2 , ε 3 , matrix P = P T > 0 and vector S, while the observer gain L is the solution of the following LMIs:
Ξ < 0 R 2 S < β m a x R 2 S > β m i n
where Ξ = Ξ 11 P E R 1 S F P P P Ξ 22 0 0 0 1 ε 1 I 0 0 0 * 1 ε 2 I 0 0 1 ε 3 I 0 and Ξ 11 = P A + A T P R 1 S C C T S T R 1 T + ϵ I 0 2 Q + δ 1 C T C + δ 2 W , Ξ 22 = δ 3 F T F γ 2 I 0 , δ 1 = δ 2 ε 1 , δ 2 = δ 2 ε 2 , δ 3 = δ 2 ε 3 , then the battery Equation (12) and the non-fragile nonlinear observer Equation (19) L = P 1 satisfy the H performance with the given attenuation γ > 0 , that is
t = 0 Γ x ˜ 2 d t < λ m a x ( P ) x ˜ ( 0 ) 2 + γ 2 ϵ 0 Γ d ( t ) 2 d t
for any Γ > 0 .
Proof. 
From Lemma 1, the observer error dynamic system Equation (20) can characterized by
Z x ˜ ( ω , t ) t = ω Z x ˜ ( ω , t ) + D α x ˜ , x ˜ = 0 μ ( ω ) Z x ˜ ( ω , t ) d ω ,
where μ ( ω ) = diag { μ 1 ( ω ) , μ 2 ( ω ) , μ 3 ( ω ) } = diag { s i n ( α 1 π ) π ω α 1 , s i n ( α 2 π ) π ω α 2 , s i n π π ω 1 } . Note that for 0 < α i < 1 , i = 1 , 2 , ω > 0 , μ i ( ω ) = s i n ( α i π ) π ω α i and s i n π π ω 1 are always positive. Therefore, let the candidate of the Lyapunov function be
V ( x ) = 0 μ ( ω ) Z x ˜ T ( ω , t ) P Z x ˜ ( ω , t ) d ω ,
and the observer gain L = P 1 .
Taking the derivative of V ( x ) along the error System (20) causes
V ˙ ( x ) = 2 0 μ ( ω ) Z x ˜ T P ω Z x ˜ ( ω , t ) + D α x ˜ d ω .
V ˙ ( x ) = 2 0 μ ( ω ) Z x ˜ T P ω Z x ˜ ( ω , t ) + D α x ˜ d ω , 2 x ˜ T P D α x ˜ d ω 2 x ˜ T [ ( P A R 1 S C Q ) x ˜ + ( P E R 1 S F ) d ] 2 x ˜ T P Δ L H ( x ^ ) C x ˜ 2 x ˜ T P Δ L H ( x ^ ) h ˜ 2 x ˜ T P Δ L H ( x ^ ) F d ,
V ˙ ( x ) 2 x ˜ T [ ( P A R 1 S C Q ) x ˜ + ( P E R 1 S F ) d ] + ε 1 x ˜ T P P T x ˜ + ε 2 x ˜ T P P T x ˜ + ε 3 x ˜ T P P T x ˜ + 1 ε 1 x ˜ T C T H T ( x ^ ) Δ L T Δ L H ( x ^ ) C x ˜ + 1 ε 2 h ˜ T H T ( x ^ ) Δ L T Δ L H ( x ^ ) h ˜ + 1 ε 3 d T F T H T ( x ^ ) Δ L T Δ L H ( x ^ ) F d , 2 x ˜ T [ ( P A R 1 S C Q + 1 2 ( ε 1 + ε 2 + ε 3 ) P P T + δ 2 2 ε 1 C T C + δ 2 2 ε 2 W ) x ˜ + ( P E R 1 S F ) d ] + d T δ 2 ε 3 F T F d .
By applying L = P 1 and the oneside Lipschitz condition Equation (15), and substituting the term x ˜ ( t ) in Equation (23), the above Equation (25) can be rewritten as Equation (26) at the top of the next page.
Whereas Δ L H ( x ^ ) δ , applying Lemma 2 and the Lipschitz condition Equation (16) to Equation (26) concludes Equation (27) at the top of the next page.
To establish then H performance, we have
V ˙ ( x ) + ϵ x ˜ T x ˜ γ 2 d T d < 0 ,
where ϵ < 2 β min 2 is used to guarantee that the matrix Ξ 11 could be negative definite. Applying the Schur complement lemma to Equation (28), it is clear to get:
x ˜ d T Ξ x ˜ d < 0 ,
where Ξ is defined by criterion Equation (21). Therefore, we see that
Ξ < 0 , R 2 S < β m a x , R 2 S > β m i n
which yield to Equation (21). Note that 0 Γ V ˙ ( x ˜ ) d t = V ( Γ ) V ( 0 ) V ( 0 ) λ max ( P )
x ˜ ( 0 ) 2 . Therefore, we have 0 Γ x ˜ T x ˜ d t λ max ( P ) ϵ x ˜ ( 0 ) 2 + γ 2 ϵ 0 Γ d T d d t for any Γ > 0 . Therefore, the criterion Equations (21) and (22) could be directly proved. This ends the proof. □

5. Experimental Validations of the SOC Estimation

5.1. Calculate the Gain L and Bound δ

In the subsection, the LMIs (21) are utilized to calculate the gain L. However, the upper bound δ of uncertainty Δ L is unknown in practice. As such, the gain L is obtained by solving the following optimization problem (21):
max L δ s . t . LMIs .
Here, the purpose is to obtain the greatest robustness for the unknown uncertainty Δ L .
Using the MATLAB LMI Toolbox [34] to solve the convex optimization problem (31), the feasible solutions of design critiron Equation (21) with ϵ = 0.1 in Theorem 1 can be obtained and some of them are list by γ = 5.36 , δ = 1.864 × 10 2 and ε 1 = 1.35 × 10 7 , ε 2 = 7.46 × 10 7 , ε 3 = 1.03 × 10 9 ,
L = 1.21 × 10 5 8.72 × 10 8 4.36 × 10 7 8.72 × 10 8 2.32 × 10 4 3.27 × 10 4 4.36 × 10 7 3.27 × 10 4 0.0055 .
In addition, the observer in [8] is used to highlight the advantages of the nonlinear observer proposed in this paper. Note that in the following figures, the lines named as Real SOC are directly obtained by the standard measurement data based on the ampere-hour counting, but the other lines are calculated by the home-made measurement data.

5.2. Numerical Implementation and Issue for Fractional-Order Observer

There are several numerical methods for fractional-order systems implementations [6,8]. The most commonly used among them is the aforementioned approximate GL definition, where the step size of h is assumed to be very small. By GL definition, the numerical implementation for proposed fractional-order observer Equation (19) is obtained as
x ^ ( k + 1 ) = ( A h α + diag ( α ) ) x ^ ( k ) i = 2 k + 1 ( 1 ) i α j x ^ ( k + 1 i ) + B h α u ( k ) + h α L h x x = x ^ ( k ) T ( y ( k ) y ^ ( k ) ) ,
for k 1 , and x ^ ( 1 ) = ( A h α + diag ( α ) ) x ^ ( 0 ) + B h α u ( 0 ) + h α L h x x = x ^ ( 0 ) T ( y ( 0 ) y ^ ( 0 ) ) , x ^ ( 0 ) = x ^ 0 ,
y ( k ) = h ( x ^ ( k ) ) + C x ^ ( k ) R 0 u ( k ) ,
where h α = diag { h α 1 h α 2 h } , α j = diag { α 1 j α 2 j 1 j } . Note that L is substituted by L + Δ L for proposed non-fragile observer (19). Note that the short memory principle [6,35] is adopted to overcome the increasing amount of calculation from the second term on the right side of Equation (33).

5.3. Experimental Results

We adopt three operation condition: the urban dynamometer driving schedule (UDDS), the new European driving cycle (NEDC), and the highway fuel economy test (HWFET) cycle to evaluate the SOC estimation accuracy of the proposed observer methods. In order to show accurately the SOC estimation results of proposed non-fragile observer for IO model (IONFO) and FO model (FONFO) and the observer for FO model in [8], the mean SOC estimation errors of the above observers are demonstrated in Table 2.

5.3.1. UDDS

The current and terminal voltage profiles are shown in Figure 6. The SOC estimation profiles and the SOC estimation error profiles are shown in Figure 7 and Figure 8, respectively.

5.3.2. NEDC

The current and terminal voltage profiles are shown in Figure 9. The SOC estimation profiles and the SOC estimation error profiles are shown in Figure 10 and Figure 11, respectively.

5.3.3. HWFET

The current and terminal voltage profiles are shown in Figure 12. The SOC estimation profiles and the SOC estimation error profiles are shown in Figure 13 and Figure 14, respectively.

5.4. SOC Estimation Results and Evaluation

In order to better discuss the effectiveness of the proposed method under various dynamic conditions, we will display the mean/maximum errors after overcoming the initial SOC error in the following table.
After that, by observing Figure 8, Figure 11 and Figure 14 and Table 2, the following conclusions can be obtained.
  • As is shown in Figure 8, Figure 11 and Figure 14 and Table 2 with the initial condition is given as [ U 1 U 2 S O C ] = [ 0 0 80 % ] , we draw the conclusion that the SOC estimation accuracy of the proposed observer for the fractional-order battery model is higher than the integer-order battery model. This implies that smaller modeling error brings higher SOC estimation accuracy. This fact is in accordance with the former studies [7,9] in fractional modeling for this battery. Meanwhile, the convergence rate of proposed observer for the battery model is slightly lower than integer-order battery model, which is the main drawback for fractional-order battery with the identified parameters of orders 0 < α 1 , α 2 < 1 .
  • In Figure 8, Figure 11 and Figure 14 and Table 2, we can see that the SOC estimation accuracy of proposed observer is superior to the observer in [8]. This is principally because of the Lipschitz constant γ mentioned in [8]. The value γ is close to 1 in [8] and the SOC estimation error bound is very sensitive to this value as is claimed. Moreover, according to eighth-order fitted curve Equation (8) of the SOC-OCV relationship, the coefficient of linear term of the polynomial is 14.13, which is larger than it in [8]. Moreover, the Lipschitz constant of U OC ( S O C ) function without the linear term is calculated by γ = 9.66 , which is much more larger than γ in [8] as well. In fact, the fifth-order polynomial used for curve fitting in [8] demonstrates poor accuracy for SOC-OCV relationship of the LiCoO 2 LIB module. Because of the high nonlinearity of SOC-OCV relationship, the observer method in [8] shows lower accuracy below expectation for the battery used in our study. On the contrary, the H performance of the proposed nonlinear observer can deal with the different Lipschitz constant in modeling for various types of batteries, which presents advanced universal property than the observer in [8].
  • It is shown in Figure 8, Figure 11 and Figure 14 and Table 2 that the SOC estimation of non-fragile observer Equation (19) shows more accuracy compared to H observer for integer-order model in [5]. This means that, if the gain drifts are taken into account in observer design procedure, the proposed non-fragile observer Equation (19) outperforms the observer in [5]. Therefore, clearly, non-fragile observer functions in tolerating additive perturbation or gain fluctuations on the observer gain. Thus, the proposed non-fragile H observer design criterion can not only restrict the effect of the modeling error and measurement noises on the state estimation, but also limit effect of the gain drifts on observer.
Based on the above analysis, we can see that proposed observer for the fractional-order battery model is more effective than the integer-order battery model. The proposed method can be generally used for SOC estimation of various types of batteries. The non-fragile observer design criterion is robust under the effects of gain drifts.

6. Conclusions

This paper presents a non-fragile nonlinear observer based on a time-varying gain algorithm which can estimate the SOC of the LIBs in electric vehicles accurately. The fractional 2nd-order RC model compared with integer 2nd-order RC model is introduced to model the charging/discharging behavior of the LIB. Moreover, a least-square method mixed with physical behavior based identification method is proposed for model parameters identification. The OCV–SOC relationship is fitted using the eighth-order polynomial for better consideration of the test battery’ property. The principle of H performance design criterion is to reduce the effect of the non-Gaussian system and measurement noises. Furthermore, a non-fragile observer design criterion is proposed to tolerate observer gain drifts to achieve a robust SOC estimation algorithm. Three operation conditions are applied to evaluate the performance of the proposed method by comparing with one former study. Experimental SOC estimation results and evaluation reveal the superiority of the proposed method.
In our study, by comparing the SOC estimation accuracy of fractional-order and integer-order, we reinforce impression that the influence of modeling error is much larger than the measurement noise. As is shown in the parameters identification section, the parameters changes rapidly in low SOC range due to the serious polarization effects. Therefore, despite of the better considering of fractional-order modeling, our future research is towards modeling and online parameters identification in low SOC range.

Author Contributions

Data curation, formal analysis, funding acquisition, Z.Z.; writing—review & editing, D.Z.; writing—original draft, methodology, N.X.; aupervision, project administration, Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 51775450) and Sichuan Science and Technology Program (2019JDRC0025).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ehsani, M.; Gao, Y.; Gay, S.E.; Emadi, A. Modern Electric, Hybrid Electric, and Fuel Cell Vehicles; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  2. Li, X.; Yuan, C.; Li, X.; Wang, Z. State of health estimation for Li-Ion battery using incremental capacity analysis and Gaussian process regression. Energy 2020, 190, 116467. [Google Scholar] [CrossRef]
  3. Tian, Y.; Xia, B.; Sun, W.; Xu, Z.; Zheng, W. A modified model based state of charge estimation of power lithium-ion batteries using unscented Kalman filter. J. Power Sources 2014, 270, 619–626. [Google Scholar] [CrossRef]
  4. Li, Z.; Xiong, R.; Mu, H.; He, H.; Wang, C. A novel parameter and state-of-charge determining method of lithium- ion battery for electric vehicles. Appl. Energy 2017, 207, 361–371. [Google Scholar] [CrossRef]
  5. Zhu, Q.; Xiong, N.; Yang, M.; Huang, R.; Hu, G. State of Charge Estimation for Lithium-Ion Battery Based on Nonlinear Observer: An H Method. Energies 2017, 10, 679. [Google Scholar] [CrossRef] [Green Version]
  6. Monje, C.A.; Gao, Y.; Gay, S.E.; Emadi, A. Fractional-Order Systems and Controls: Fundamentals and Applications; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  7. Xu, J.; Mi, C.C.; Cao, B.; Cao, J. A new method to estimate the state of charge of lithium-ion batteries based on the battery impedance model. J. Power Sources 2013, 233, 277–284. [Google Scholar] [CrossRef]
  8. Wang, B.; Liu, Z.; Li, S.E.; Moura, S.J.; Peng, H. State-of-Charge Estimation for Lithium-Ion Batteries Based on a Nonlinear Fractional Model. IEEE Trans. Control. Syst. Technol. 2016, 25, 3–11. [Google Scholar] [CrossRef]
  9. Liu, C.; Liu, W.; Wang, L.; Hu, G.; Ma, L.; Ren, B. A new method of modeling and state of charge estimation of the battery. J. Power Sources 2016, 320, 1–12. [Google Scholar] [CrossRef]
  10. Xiao, R.; Shen, J.; Li, X.; Yan, W.; Pan, E.; Chen, Z. Comparisons of modeling and state of charge estimation for lithium-ion battery based on fractional order and integral order methods. Energies 2016, 9, 184. [Google Scholar] [CrossRef]
  11. Zhong, F.; Li, H.; Zhong, S.; Zhong, Q.; Yin, C. An SOC estimation approach based on adaptive sliding mode observer and fractional order equivalent circuit model for lithium-ion batteries. Commun. Nonlinear Sci. Numer. Simul. 2014, 24, 127–144. [Google Scholar] [CrossRef]
  12. Zhong, Q.; Zhong, F.; Cheng, J.; Li, H.; Zhong, S. State of charge estimation of lithium-ion batteries using fractional order sliding mode observer. ISA Trans. 2016, 66, 448–459. [Google Scholar] [CrossRef]
  13. Xia, B.; Chen, C.; Tian, Y.; Sun, W.; Xu, Z.; Zheng, W. A novel method for state of charge estimation of lithium-ion batteries using a nonlinear observer. J. Power Sources 2014, 270, 359–366. [Google Scholar] [CrossRef]
  14. Liu, L.; Wang, L.Y.; Chen, Z.; Wang, C.; Lin, F.; Wang, H. Integrated system identification and state-of-charge estimation of battery systems. IEEE Trans. Energy Convers. 2013, 28, 12–23. [Google Scholar] [CrossRef]
  15. Lin, C.; Mu, H.; Xiong, R.; Shen, W. A novel multi-model probability battery state of charge estimation approach for electric vehicles using H-infinity algorithm. Appl. Energy 2016, 166, 76–83. [Google Scholar] [CrossRef]
  16. Zhang, F.; Liu, G.; Fang, L.; Wang, H. Estimation of battery state of charge with H observer: Applied to a robot for inspecting power transmission lines. IEEE Trans. Ind. Electron. 2012, 59, 1086–1095. [Google Scholar] [CrossRef]
  17. Ibrahima, N.; Taous-Meriem, L.K.; Mohamed, D. H Adaptive observer for nonlinear fractional-order systems. Int. J. Adapt. Control. Signal Process. 2017, 31, 314–331. [Google Scholar] [CrossRef]
  18. Heleschewitz, D.; Matignon, D. Diffusive Realisations of Fractional Integrodifferential Operators: Structural Analysis under Approximation. In IFAC Proceedings Volumes; Elsevier: Amsterdam, The Netherlands, 1998; Volume 31, pp. 227–232. [Google Scholar]
  19. Trigeassou, J.C.; Maamri, N.; Sabatier, J.; Oustaloup, A. A Lyapunov approach to the stability of fractional differential equations. Signal Process. 2010, 91, 437–445. [Google Scholar] [CrossRef]
  20. Boroujeni, E.A.; Momeni, H.R. Non-fragile nonlinear fractional order observer design for a class of nonlinear fractional order systems. Signal Process. 2012, 92, 2365–2370. [Google Scholar] [CrossRef]
  21. Zhang, L.; Hu, X.; Wang, Z.; Sun, F.; Dorrell, D.G. Fractional-order modeling and State-of-Charge estimation for ultracapacitors. J. Power Sources 2016, 314, 28–34. [Google Scholar] [CrossRef]
  22. Waag, W.; Fleischer, C.; Uwe, D. Critical review of the methods for monitoring of lithium-ion batteries in electric and hybrid vehicles. J. Power Sources 2014, 258, 321–339. [Google Scholar] [CrossRef]
  23. Dorato, P. Non-fragile controller design: An overview. In Proceedings of the American Control Conference, Philadelphia, PA, USA, 26 June 1998; pp. 2829–2831. [Google Scholar] [CrossRef]
  24. Boroujeni, E.A.; Momeni, H.R. An iterative method to design optimal non-fragile H observer for Lipschitz nonlinear fractional-order systems. Nonlinear Dyn. 2015, 84, 1801–1810. [Google Scholar] [CrossRef]
  25. Sakthivel, R.; Anbuvithya, R.; Mathiyalagan, K.; Prakash, P. Combined H and passivity state estimation of memristive neural networks with random gain fluctuations. Neurocomputing 2015, 168, 1111–1120. [Google Scholar] [CrossRef]
  26. Wang, J.; Lu, J.; Li, C. Observer-based robust stabilisation of a class of non-linear fractional-order uncertain systems: An linear matrix inequalitie approach. IET Control. Theory Appl. 2012, 2757–2764. [Google Scholar] [CrossRef]
  27. Mathiyalagan, K.; Park, J.H.; Jung, H.Y.; Sakthivel, R. Non-fragile Observer-Based H Control for Discrete-Time Systems Using Passivity Theory. Circuits Syst. Signal Process. 2015, 34, 2499–2516. [Google Scholar] [CrossRef]
  28. Zhou, K.; Doyle, J.C.; Keith, G. Robust and Optimal Control; Prentice-Hall: Hoboken, NJ, USA, 1996. [Google Scholar]
  29. Gao, L.; Jiang, X.; Wang, D. Observer-based robust finite time H sliding mode control for Markovian switching systems with mode-dependent time-varying delay and incomplete transition rate. ISA Trans. 2015, 61, 29–48. [Google Scholar] [CrossRef] [PubMed]
  30. Kochowski, S.; Nitsch, K. Description of the frequency behaviour of metal-SiO2-GaAs structure characteristics by electrical equivalent circuit with constant phase element. Thin Solid Films 2002, 415, 133–137. [Google Scholar] [CrossRef]
  31. Westerlund, S.; Ekstam, L. Capacitor theory. IEEE Trans. Dielectr. Electr. Insul. 1994, 1, 826–839. [Google Scholar] [CrossRef]
  32. Sabatier, J.; Agrawal, O.P.; Tenreiro Machado, J.A. Advances in Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  33. He, H.; Zhang, X.; Xiong, R.; Xu, Y.; Guo, H. Online model-based estimation of state-of-charge and open-circuit voltage of lithium-ion batteries in electric vehicles. Energy 2012, 39, 310–318. [Google Scholar] [CrossRef]
  34. Boyd, S.; Ghaoui, L.E.; Feron, E.; Balakrishnan, V. Linear Matrix Inequalities in System and Control Theory; SIAM: Philadelphia, PA, USA, 1994. [Google Scholar]
  35. Mu, H.; Xiong, R.; Zheng, H.; Chang, Y.; Chen, Z. A novel fractional order model based state-of-charge estimation method for lithium-ion battery. Appl. Energy 2017, 384–393. [Google Scholar] [CrossRef]
Figure 1. The schematic diagram of the fractional 2nd-order RC model.
Figure 1. The schematic diagram of the fractional 2nd-order RC model.
Energies 14 04771 g001
Figure 2. Block diagram of the test bench.
Figure 2. Block diagram of the test bench.
Energies 14 04771 g002
Figure 3. Identification results and mean value in [0.2, 1] of parameters R 0 , R 1 , R 2 , C 1 , C 2 .
Figure 3. Identification results and mean value in [0.2, 1] of parameters R 0 , R 1 , R 2 , C 1 , C 2 .
Energies 14 04771 g003
Figure 4. Measurement data and fitted curve of OCV vs. SOC.
Figure 4. Measurement data and fitted curve of OCV vs. SOC.
Energies 14 04771 g004
Figure 5. The implement flowchart of non-fragile H observer.
Figure 5. The implement flowchart of non-fragile H observer.
Energies 14 04771 g005
Figure 6. Current and terminal voltage profiles of UDDS.
Figure 6. Current and terminal voltage profiles of UDDS.
Energies 14 04771 g006
Figure 7. SOC estimation profiles for UDDS.
Figure 7. SOC estimation profiles for UDDS.
Energies 14 04771 g007
Figure 8. SOC estimation error profiles for UDDS.
Figure 8. SOC estimation error profiles for UDDS.
Energies 14 04771 g008
Figure 9. Current and terminal voltage profiles of NEDC.
Figure 9. Current and terminal voltage profiles of NEDC.
Energies 14 04771 g009
Figure 10. SOC estimation profiles for NEDC.
Figure 10. SOC estimation profiles for NEDC.
Energies 14 04771 g010
Figure 11. SOC estimation error profiles for NEDC.
Figure 11. SOC estimation error profiles for NEDC.
Energies 14 04771 g011
Figure 12. Current and terminal voltage profiles of HWFET.
Figure 12. Current and terminal voltage profiles of HWFET.
Energies 14 04771 g012
Figure 13. SOC estimation profiles for HWFET.
Figure 13. SOC estimation profiles for HWFET.
Energies 14 04771 g013
Figure 14. SOC estimation error profiles for HWFET.
Figure 14. SOC estimation error profiles for HWFET.
Energies 14 04771 g014
Table 1. Identification results of the parameters.
Table 1. Identification results of the parameters.
α 1 α 2 R 0 R 1 R 2 C 1 C 2
0.910.836.9742.6048.20110.259125.688
Table 2. Mean/maximum SOC estimation errors of observers after convergence.
Table 2. Mean/maximum SOC estimation errors of observers after convergence.
Operation ConditionUDDSNEDCHWFET
IONFO0.85%/2.62%0.71%/1.64%0.51%/2.72%
FONFO0.73%/1.64%0.53%/1.24%%0.41%/1.83%
Observer in Ref. [8]0.83%/3.44%0.94%/2.23%0.55%/3.58%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Z.; Zhou, D.; Xiong, N.; Zhu, Q. Non-Fragile H Nonlinear Observer for State of Charge Estimation of Lithium-Ion Battery Based on a Fractional-Order Model. Energies 2021, 14, 4771. https://doi.org/10.3390/en14164771

AMA Style

Zhang Z, Zhou D, Xiong N, Zhu Q. Non-Fragile H Nonlinear Observer for State of Charge Estimation of Lithium-Ion Battery Based on a Fractional-Order Model. Energies. 2021; 14(16):4771. https://doi.org/10.3390/en14164771

Chicago/Turabian Style

Zhang, Zhongwei, Dan Zhou, Neng Xiong, and Qiao Zhu. 2021. "Non-Fragile H Nonlinear Observer for State of Charge Estimation of Lithium-Ion Battery Based on a Fractional-Order Model" Energies 14, no. 16: 4771. https://doi.org/10.3390/en14164771

APA Style

Zhang, Z., Zhou, D., Xiong, N., & Zhu, Q. (2021). Non-Fragile H Nonlinear Observer for State of Charge Estimation of Lithium-Ion Battery Based on a Fractional-Order Model. Energies, 14(16), 4771. https://doi.org/10.3390/en14164771

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop