1) An Equivalent Circuit Model Analysis
1) An Equivalent Circuit Model Analysis
1) An Equivalent Circuit Model Analysis
Jie Su , Maosong Lin, Shunli Wang, Jin Li, James Coffie-Ken and Fei Xie
Abstract
According to the demand of vehicle lithium-ion battery pack, the splice equivalent circuit model is constructed. First, a
joint experiment of intermittent discharge and hybrid power pulse characterization, basis of the requirements of para-
meter identification for the model, is designed to identify the parameters. Study shows that it can identify the para-
meters. Second, the splice equivalent circuit model of vehicle lithium-ion battery pack is simulated by MATLAB/Simulink,
which shows the model is feasible to describe the vehicle lithium-ion battery pack. The simulation error of ohmic inter-
nal resistance R0 should be less than 0.05 mO. Study suggests that the vehicle lithium-ion battery pack has a stable dis-
charge period within the state-of-charge range of [20%, 80%]. However, when stage of charge is below 20%, vehicle
lithium-ion battery pack is no longer stable and the parameters of the splice equivalent circuit model change
dramatically.
Keywords
Lithium-ion battery pack, pure electric vehicle, splice equivalent circuit model, hybrid power pulse characterization test,
parameter identification
Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License
(http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without
further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/
open-access-at-sage).
194 Measurement and Control 52(3-4)
nRC equivalent circuit model is proposed in Lai et al.26 resistance Rs is connected in parallel to characterize the
The high-order RC model is robust to the variations of self-discharge effect of the lithium-ion battery pack. Ro
model parameters and sensor errors. Wang et al.27 pro- is the benchmark ohmic internal resistance of the bat-
posed a splice equivalent circuit model (S-ECM) for the tery pack, which is used to describe the ohmic effect of
aerial lithium-ion battery pack, which can be used for the lithium-ion battery pack. Rd is a small variable of
the accurate estimation of state of charge (SoC). ohmic internal resistance caused by discharging. Rc is a
However, the difference between the monomers of small variable of ohmic internal resistance caused by
lithium-ion battery pack was not considered during the charging. They can reflect the dynamic characteristics
equivalent circuit modeling process in most of the of lithium-ion battery pack during discharging/charg-
above-mentioned researches. Disadvantages still exist ing. Rp is used to characterize the polarization resis-
in some of the simulated description of lithium-ion bat- tance of lithium-ion battery pack, and Cp is the
tery packs. In order to solve these issues, this paper polarization capacitance of lithium-ion battery pack. Rp
study the equivalent circuit model based on the first- and Cp constitute a first-order RC parallel circuit to
order RC model. First, the time-varying voltage source represent the polarization effect of lithium-ion battery
and the time-varying resistance are connected in series. pack and simulate the relaxation effect of lithium-ion
They are used to describe the effects of the equilibrium battery pack, thus realizing the transient response
state between monomers on open-circuit voltage (OCV) description of lithium-ion battery pack. Ud and Rd are
and ohmic internal resistance, respectively. Then, the used to describe the equilibrium state between internal
ohmic internal resistance is divided into benchmark interconnect monomers. They reflect the effects of the
ohmic internal resistance and charge/discharge differen- equilibrium state between monomers on OCV and
tial ohmic internal resistance to simulate the change ohmic internal resistance, respectively. UL is the termi-
characteristics of ohmic internal resistance of charging/ nal voltage of lithium-ion battery pack for pure electric
discharging. Meanwhile, the self-discharge internal vehicle, and I(t) are the current value at both ends.
resistance is used to simulate the self-discharge effect of
vehicle battery packs. After that, this paper discusses
one method to identify these parameters. Finally, this Parameter identification
model is applied to vehicle lithium-ion battery packs. In order to obtain model parameters, the intermittent
The rest of the paper is organized as follows. In sec-
discharge and hybrid power pulse characterization
tion ‘‘Theoretical analysis,’’ the construction method of
(HPPC) joint experiments of vehicle lithium-ion battery
equivalent circuit model and parameter identification
pace are designed. The process of designing HPPC is
method are discussed. In section ‘‘Experimental analy-
shown in Figure 2.
sis,’’ the proposed model is verified experimentally.
According to Figure 2, the parameters of lithium-ion
battery pack for pure electric vehicle can be obtained.
Theoretical analysis The counting method is as follows:
Construction of equivalent circuit model 1. The calculation of Uoc: When the lithium-ion bat-
According to the research needs of the pure electric tery pack is fully shelved, the internal electrochemi-
vehicle power battery packs, with existing research, cal reaction is stable, which is considered Uoc = UL
considering ohmic effect, polarization effect, and self- at this time. Therefore, in Figure 2, Uoc = U1.
discharge effect of the battery packs, the models-ECM 2. Calculate Ro: Ro describes the instantaneous vol-
is constructed. tage drop from U1 to U2 in Figure 2. Therefore
Figure 1 is the S-ECM model of vehicle lithium-ion
battery pack. Uoc is the OCV of the battery pack, which
represents the change of the terminal voltage OCV of U1 U2
R0 = ð1Þ
the vehicle lithium-ion battery pack, and a large I
Su et al. 195
,
Table 1. Experimental requirements.
8
8 Parameter Experimental requirements
8
89
8
Ro HPPC test28,29 in intermittent discharge
8 Rp HPPC test28,29 in intermittent discharge
8 Cp HPPC test28,29 in intermittent discharge
Rd Calculation of pulse discharge combined with
8
ohmic resistance
Rc Calculation of pulse charge combined with ohmic
W W W W resistance
WV
Ud Combining the evaluation results of equilibrium
Figure 2. HPPC testing process. states to get function
Rd Combining the evaluation results of equilibrium
states to get function
In formula (1), I is the current value of 1C5A pulse
discharge. SoC: state of charge; OCV: open-circuit voltage; HPPC: hybrid power
pulse characterization.
3. The calculation of Rd and Rc: The ohmic internal
resistance of vehicle lithium-ion batteries will cause
the terminal voltage to rise and fall instantaneously 6. The calculation of Ud and Rd: Ud and Rd are used
during pulse charging and discharging. However, to describe the equilibrium state between the mono-
there are differences in ohmic resistance to charg- mer of pure vehicle lithium-ion battery pack. Their
ing and discharging. The difference is expressed by calculation methods are described in formula (5)
Rd and Rc. Their formulas are as follows
8
Rd = U3 U 4
R0 >
> Ud = e UOC
I
U7 U8 ð2Þ >
> Rd = e R0
Rc = I R0 >
>
< n h
P i2
1 Uci EðUc Þ
e= n EðUc Þ ð5Þ
where U3, U4, U7, and U8 in formula (2) are shown in >
> i=1
Figure 2, and I is described in formula (1). >
> P
n
>
> ð 1
: E U cÞ = n Uci
i=1
4. The extraction of Rs: For the sake of obtaining Rs,
CC-CV is used to charge the pure vehicle lithium- In formula (5), n denotes the number of lithium-ion
ion battery packs to make SoC = 1. Then the bat- battery single monomer in the pure vehicle lithium-ion
tery is kept on hold for 1 h and the OCV of the bat- battery pack and Uci denotes the voltage value of ith
tery pack UOCRs is measured. The SoC value is lithium-ion battery pack single monomer only. e is the
SoCRs at this moment. Shelving for 30 days, the equilibrium state between the monomers, which is used
OCV of battery pack UOCRs2 is measured. The SoC to describe the degree of voltage discrepancy between
value is SoCRs2 at this moment. It is considered monomers.
that the whole self-discharge process is constant
current discharge and the current is Is. Assume that
the nominal capacity of the battery pack is Qn.
Experimental analysis
Accordingly, the calculation method of Rs can be Experimental requirements
obtained from formula (3)
There are nine parameters of the S-ECM model which
need to be identified. The identification requirements of
(
Rs = ðUOCRs2 UOCRs Þ these parameters are not identical exactly. The experi-
Is
DQ 3600ðSoCRs2 SoCRs ÞQn
ð3Þ mental requirements for identification are displayed in
Is = Dt = Dt Table 1.
5. Calculation of Rp and Cp: As shown in the HPPC
test shown in Figure 2, the formula for calculating Experimental design
the RC response parameter is as follows
To verify the correctness of the model and the accuracy
of parameter identification, the parameters of lithium
8
>
> Rp = U5 U
I
4 cobalt oxide battery commonly used in pure electric
< Cp = Rtp vehicles is identified by the power battery test system
ð4Þ (company: MACCOR; model: S4000; voltage accuracy:
> t = t3 t2
>
: U U
ln 1 5 4 0.02%; current accuracy: 0.02%). The parameters of
U1 U4
lithium-ion battery in this paper are shown in Table 2.
196 Measurement and Control 52(3-4)
Table 2. Battery parameters. Starting from SoC = 1, HPPC test is conducted every
5% SoC.
Project Parameter
,&$
8RF
8RF9
62&
R0
2.5
2.0
R0 (mΩ)
1.5
1.0
0.5
and discharging. By comparing the two curves, it is fluctuation is small. When SoC is below 20%, Cp
found that the values of Rd and Rc are identical basi- increases rapidly, far beyond 3000 F. After that, Ud and
cally at the same SoC value. When SoC is below 20%, Rd are analyzed in Table 3.
both Rd and Rc increase rapidly. From Table 3, it can be found that the seven bat-
Figure 7 shows that Rc and SoC are not significantly teries used in this paper have good consistency at the
correlated; for obtaining the relationship between Rp beginning of discharge. As deep discharge proceeds, the
and SoC, we can continue the least square fitting. battery pack becomes inconsistent slowly. When SoC is
When the SoC is above 20%, the polarization capaci- less than 15%, the inconsistency increases rapidly, Ud
tance Cp stays between 2000 and 3000 F, and the and Rd increase accordingly. Therefore, Ud and Rd can
198 Measurement and Control 52(3-4)
5G
5F
5P¡
100% 4.1292 4.1312 4.1313 4.1303 4.1283 4.1301 4.1337 0.0044 0.0000
95% 4.0864 4.0880 4.0877 4.0871 4.0862 4.0868 4.0904 0.0030 0.0000
90% 4.0462 4.0484 4.0472 4.0467 4.0467 4.0464 4.0508 0.0040 0.0000
85% 4.0053 4.0065 4.0048 4.0048 4.0056 4.0034 4.0080 0.0031 0.0000
80% 3.9630 3.9611 3.9578 3.9594 3.9626 3.9568 3.9616 0.0087 0.0000
75% 3.9245 3.9233 3.9201 3.9211 3.9242 3.9193 3.9237 0.0068 0.0000
70% 3.8884 3.8893 3.8870 3.8875 3.8891 3.8868 3.8899 0.0023 0.0000
65% 3.8619 3.8633 3.8618 3.8621 3.8632 3.8609 3.8642 0.0020 0.0000
60% 3.8402 3.8417 3.8405 3.8408 3.8416 3.8399 3.8429 0.0017 0.0000
55% 3.8222 3.8239 3.8224 3.8227 3.8235 3.8220 3.8246 0.0015 0.0000
50% 3.8083 3.8099 3.8085 3.8086 3.8097 3.8082 3.8104 0.0013 0.0000
45% 3.7945 3.7947 3.7932 3.7941 3.7953 3.7928 3.7948 0.0013 0.0000
40% 3.7797 3.7765 3.7737 3.7755 3.7783 3.7726 3.7761 0.0095 0.0000
35% 3.7600 3.7563 3.7537 3.7552 3.7583 3.7531 3.7558 0.0095 0.0000
30% 3.7465 3.7426 3.7405 3.7419 3.7446 3.7391 3.7426 0.0097 0.0000
25% 3.7272 3.7220 3.7197 3.7214 3.7245 3.7185 3.7218 0.0136 0.0000
20% 3.7053 3.6987 3.6967 3.6986 3.7015 3.6956 3.6986 0.0167 0.0000
15% 3.6913 3.6893 3.6870 3.6887 3.6902 3.6848 3.6896 0.0075 0.0000
10% 3.6257 3.5822 3.5639 3.5762 3.6035 3.5494 3.5813 1.0577 0.0005
5% 3.4016 3.2860 3.2266 3.2652 3.3453 3.1688 3.2838 10.5110 0.0062
describe the inconsistency between monomers. It also experimental data onto section ‘‘Experimental result’’ is
could quantitatively analyze their effects on OCV and imported into the simulation circuit and the ohmic
ohmic internal resistance. internal resistance R0 is taken as the observation. The
simulation results are shown in Figure 8.
It can be seen from Figure 8 that the simulation
Simulation and error analysis curve is highly consistent with the real value curve. For
The S-ECM shown in Figure 1 is simulated in the sake of quantitatively analyzing the simulation
MATLAB/Simulink to verify whether it can be used to results, the estimated error is calculated, as shown in
describe vehicle lithium-ion battery pack. The Figure 9.
Su et al. 199
Rp
20 Cp
19
8000
18
17
6000
Rp(mΩ)
Cp (F)
16
15
14 4000
13
12 2000
UHDO
VLPXODWLRQ
5P¡
HUURU
(UURUP¡
62&
packs based on the state of balance estimation. Appl battery pack based on the real-time detection and filter-
Energ 2016; 166: 44–58. ing. J Clean Prod 2018; 185: 187–197.
18. Mathew M, Janhunen S, Rashid M, et al. Comparative 25. Wu J, Xing XW, Liu XT, et al. Energy management strat-
analysis of lithium-ion battery resistance estimation tech- egy for grid-tied microgrids considering the energy storage
niques for battery management systems. Energies 2018; efficiency. IEEE T Ind Electron 2018; 65(12): 9539–9549.
11(6): 1490. 26. Lai X, Zheng YJ and Sun T. A comparative study of dif-
19. Chen J, Ouyang Q, Xu CF, et al. Neural network-based ferent equivalent circuit models for estimating state-of-
state of charge observer design for lithium-ion batteries. charge of lithium-ion batteries. Electrochim Acta 2018;
IEEE T Contr Syst T 2018; 26(1): 313–320. 259: 566–577.
20. Afshari HH, Attari M, Ahmed R, et al. Reliable state of 27. Wang SL, Fernandez C, Liu XH, et al. The parameter
charge and state of health estimation using the smooth identification method study of the Splice equivalent cir-
variable structure filter. Control Eng Pract 2018; 77: 1–14. cuit model for the aerial lithium-ion battery pack. Meas
21. Lai X, Qin C, Gao WK, et al. A state of charge estimator Control 2018; 51(5–6): 125–137.
based extended Kalman filter using an electrochemistry- 28. Yu QQ, Xiong R, Wang LY, et al. A comparative study
based equivalent circuit model for lithium-ion batteries. on open circuit voltage models for lithium-ion batteries.
Appl Sci 2018; 8(9): 1592. Chin J Mech Eng 2018; 31(1): 65.
22. Zhao LH, Ji GH and Liu ZY. Design and experiment of 29. INEE Laboratory. FreedomCAR battery test manual for
nonlinear observer with adaptive gains for battery state power-assist hybrid electric vehicles. U.S. Department of
of charge estimation. Energies 2017; 10(12): 2046. Energy, 2003, pp. 3–6, https://avt.inl.gov/sites/default/
23. Wei JW, Dong GZ, Chen ZH, et al. System state estima- files/pdf/battery/freedomcar_manual_04_15_03.pdf
tion and optimal energy control framework for multicell 30. Barai A, Uddin K, Widanage WD, et al. A study of the
lithium-ion battery system. Appl Energ 2017; 187: 37–49. influence of measurement timescale on internal resistance
24. Wang SL, Fernandez C, Chen MJ, et al. A novel safety characterisation methodologies for lithium-ion cells. Sci
anticipation estimation method for the aerial lithium-ion Rep 2018; 8: 21.