CN109061508A - A kind of estimation method of electric automobile lithium battery SOH - Google Patents
A kind of estimation method of electric automobile lithium battery SOH Download PDFInfo
- Publication number
- CN109061508A CN109061508A CN201811056160.2A CN201811056160A CN109061508A CN 109061508 A CN109061508 A CN 109061508A CN 201811056160 A CN201811056160 A CN 201811056160A CN 109061508 A CN109061508 A CN 109061508A
- Authority
- CN
- China
- Prior art keywords
- soc
- lithium battery
- particle
- internal resistance
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Secondary Cells (AREA)
Abstract
The present invention relates to the estimation methods of electric automobile lithium battery SOH a kind of, which comprises the following steps: 1) the Order RC equivalent model of construction force lithium battery；2) parameter identification of Order RC equivalent model, including ohmic internal resistance R0And Order RC link parameter；3) ohmic internal resistance R is calculated using improvement UPF method0；4) SOH value of dynamic lithium battery is calculated.Compared with prior art, the present invention has many advantages, such as that tracking prediction ability is strong, estimated accuracy is high.
Description
Technical field
The present invention relates to dynamic lithium battery technical fields, more particularly, to the estimation side of electric automobile lithium battery SOH a kind of
Method.
Background technique
Global Oil resource is increasingly depleted, and the popularize caused problem of environmental pollution and energy crisis of fuel-engined vehicle cause
The highest attention of people.Electric car promotes countries in the world to throw as a kind of environmentally friendly, resource-conserving vehicles
Enter the research and development that a large amount of manpower and material resources are engaged in electric car.Core component of the power battery as electric car, ecology potential must
Key effect is played to the performance of electric car.Currently, application prospect of the lithium ion battery on electric car is best, to its property
The research of energy also becomes hot spot in recent years, and the state parameter that characterize battery performance mainly has state-of-charge (State of
Charge, SOC) and health status (State of Health, SOH).State-of-charge SOC characterizes the energy storage capacity of battery, shadow
Performances, the health status SOH such as the course continuation mileage of electric car is rung then to characterize the service life cycle of battery, directly affect electricity
The reliability and safety of electrical automobile.Research compared to SOC, domestic and foreign scholars are less to the research of SOH, and SOH is that battery is ground
The weak link studied carefully, accurately monitors SOH and estimation helps to improve battery management system (Battery
Management System, BMS), SOC estimated accuracy is improved, overcharge, overdischarge is avoided, facilitates fault diagnosis and safety
Early warning is of great significance to electric car general safety stable operation.
Currently, the main approaches in terms of lithium battery SOH estimation have: defining method (" Intelligent
prognostics for battery health monitoring based on sample entropy."Expert
Systems with Applications ", " State-of-Health estimation of Lithium-ion battery
Based on interval capacity "), also known as volumetric method, direct electric discharge, be substantially by the capacity of battery, electricity,
The proportionate relationship of the Feature changes such as power carries out SOH assessment.Internal resistance method (" grind by power battery pack health status evaluation method
Study carefully ", " State-of-health monitoring of lithium-ion batteries in electric vehicles
By on-board internal resistance estimation "), mainly by establishing the pass between internal resistance and SOH
System is to estimate SOH.Neural network (Neural Network, NN) (" A Neural Network Based State-of-
Health Estimation of Lithium-ion Battery in Electric Vehicles ", " grey neural network
Model estimation on line lithium ion battery SOH "), SOH is obtained by obtaining a large amount of experimental data come training pattern.Kalman's filter
Wave method (Kalman Filter, KF) (" the 18650 lithium ion battery health estimations based on UKF ", " based on the lithium of AUKF from
Sub- cell health state estimation "), battery-based state-space model passes through the SOH of recursion iterative estimate battery.Particle filter
(" a kind of lithium ion battery health forecast algorithm based on particle filter " " is based on SVR- to wave method (Particle Filter, PF)
The unmanned plane cell health state diagnosis research of PF ", " State-of-health estimation of lithium-ion
battery packs in electric vehicles based on genetic resampling particle
Filter "), basic thought is with one group of sample come the Posterior probability distribution of approximate representation system, then approximate using this
Indicate the state to estimate nonlinear system.Definition method is the effective ways for the evaluation cell health state that industry is generally acknowledged, but needs
By battery off-line measurement, difficulty is realized for vehicle mounted dynamic battery group, and special test load is needed to carry out battery
Repeated charge, experimental period is long, and practical value is not high.It is related to its internal resistance that internal resistance method establishes cell degradation degree
Property, its health status is estimated by accurately measuring internal resistance of cell value.Internal resistance method can be realized power battery health status
On-line Estimation, but the internal resistance of cell belongs to small signal, it is relatively difficult to accurately measure, so this method obtains not yet at present
To a large amount of application.
The fundamental characteristics of neural network is non-linear, and has very strong generalization ability and parallel processing capability, is suitble to mould
Quasi- electric automobile power battery working condition, versatility are stronger.But it is to need to carry out battery the shortcomings that neural network
The precision that a large amount of experiment is trained model with obtaining comprehensive sample data, and estimates is largely by training number
According to influence and be difficult to realize on-line training.Kalman filtering algorithm is essentially recursive feedback, be substantially carried out the time update and
Measuring state updates.Since Kalman filtering is only applicable to linear system, at present using it is more be all its innovatory algorithm, such as expand
Open up Kalman filtering (Extended Kalman Filter, EKF) and Unscented kalman filtering (Unscented Kalman
Filter, UKF) etc., it is aided with other improvement means, to adapt to the nonlinear system of battery.But EKF and UKF are only applicable to
Non-linear Gaussian system effectively overcomes though UKF has abandoned the traditional method that EKF linearizes nonlinear function
The defect that EKF estimated accuracy is low, stability is poor.But UKF is still the posterior probability density with Gaussian Profile come approximation system state
Function, if error is larger in the case where non-gaussian.
Particle filter is able to solve nonlinear and non-Gaussian problem, and nonlinear and non-Gaussian system as research battery is come
It says, effective method of can yet be regarded as, and application of the PF algorithm in terms of lithium battery SOH estimation be not also extensive, has larger
Research space.In recent years, " a kind of lithium ion battery health forecast algorithm based on particle filter " is proposed with lumped parameter mould
Particle filter algorithm based on type predicts battery SOH." the unmanned plane cell health state diagnosis research based on SVR-PF "
Propose it is a kind of based on support vector regression particle filter (Support Vector Machine Particle Filter,
SVR-PF cell health state diagnostic method).In order to solve the degenerate problem of particle filter algorithm, " State-of-health
estimation of lithium-ion battery packs in electric vehicles based on genetic
Resampling particle filter " genetic algorithm is combined with particle filter algorithm, most using forgetting factor recursion
Small square law carries out parameter identification, so that the estimation for carrying out battery pack internal resistance is tested to evaluate SOH.Generally speaking, particle filter
Algorithm has certain superiority for the solution of nonlinear and non-Gaussian system, but often will appear particle during realization
The problems such as degeneration, particle Loss of diversity, and calculation amount is larger, needs to make improvements.
Summary of the invention
It is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide a kind of electric car lithium electricity
The estimation method of pond SOH.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of estimation method of electric automobile lithium battery SOH, comprising the following steps:
1) the Order RC equivalent model of construction force lithium battery；
2) parameter identification of Order RC equivalent model, including ohmic internal resistance R0And Order RC link parameter；
3) ohmic internal resistance R is calculated using improvement UPF method0；
4) SOH value of dynamic lithium battery is calculated.
In the step 1), Order RC equivalent model is that polarization resistance in parallel is added based on Thevenin model
RpiWith polarization capacity CpiOrder RC equivalent-circuit model is constructed, and by the U of two RC linksp1And Up2, ohmic internal resistance R0With
SOC constructs observational equation as observed quantity as state variable x, the end voltage U of battery.
The circuit equation of the Order RC equivalent model are as follows:
U=UOCV(SOC)-IR0-Up1-Up2
Wherein, R0For ohmic internal resistance, UOCVIt (SOC) is the open-circuit voltage under SOC, I is the charging and discharging currents of battery, and U is end
Voltage, Up1And Up2For the voltage of two RC links, UpiFor polarizing voltage.
The observational equation are as follows:
Wherein, T is sampling time, QNFor battery rated capacity, IkFor charging and discharging currents, w1,k、w2,k、w3,k、w4,kTo be
System state-noise, νkFor systematic observation noise, Cp1、Cp2The capacitor of respectively two RC links, Rp1、Rp2Respectively two RC rings
The resistance of section, subscript k indicate sampling instant.
In the step 2), the parameter in second-order model is recognized using HPPC pulse charge-discharge test, specifically:
21) it is tested by volume test, demarcates the rated capacity Q of lithium batteryN；
22) relational expression between low current charge-discharge test fitting acquisition SOC and open-circuit voltage OCV is carried out；
23) parameter in Order RC equivalent-circuit model is recognized by HPPC pulse charge-discharge test, obtains charge and discharge
The ohmic internal resistance of each SOC point on direction, and be fitted to obtain the relationship expression of SOC and ohmic internal resistance on charge and discharge direction
Formula.
The step 3) specifically includes the following steps:
31) sampling instant k=0 is enabled, total number of particles N, and init state variable x are set0, covariance matrix P and just
Beginning weightI=1,2 ..., N；
32) sampling instant k=k+1 is enabled, particle is predicted using UKF algorithm, complete the estimation of weighting particle and is adopted
Sample improves UPF algorithm by UKF and generates the importance density function of PF algorithm to carry out forecast updating to particle；
33) weight of all N number of particles is calculated separatelyAnd it is normalized to obtain normalization weight
34) number of effective particle is countedGiven threshold Nth=2N/3, if Neff> Nth,
Step 35) is then carried out, if Neff≤Nth, then to particle collectionCarry out resampling；
35) optimal estimation and the corresponding variance matrix of each particle of k sample moment state are obtained；
36) terminate when sampling number reaches setting value, otherwise return step 32).
In the step 33), i-th of particle weightCalculating formula are as follows:
Wherein, S is the standard deviation that observational equation introduces noise, Upred_ukfFor the mean value of observation prediction.
The resampling of the step 34) specifically:
341) to all normalization weightsIt is ranked up from small to large, and will be with normalization weightCorresponding grain
Sub- estimated valueWith true valueDifferenceIt is ranked up；
342) weight will be normalizedWith differenceQuotientIt is ranked up from small to large, if existing in quotient
Two groups of equal data then delete more new sort after lesser one group of weight of normalization；
343) setting survival of the fittest threshold value ωthFor quotientThe 2/5 of average value, ifThen retain correspondence
Particle weight and particle estimated value otherwise delete corresponding particle weight and particle estimated value, and more new particle collection.
In the step 4), the calculating formula of the SOH value of dynamic lithium battery are as follows:
Wherein, REOLIt (SOC) is the internal resistance value at the end of battery life, RNOWIt (SOC) is the current internal resistance value of battery, i.e., it is optimal
Particle estimated valueIn ohmic internal resistance, RNEW(SOC) internal resistance value when dispatching from the factory for battery.
Compared with prior art, the invention has the following advantages that
The present invention has initially set up dynamic lithium battery Order RC equivalent-circuit model, experimental data is then based on, using most
Small two, which multiply curve-fitting method, picks out parameter in model.It proposes according to the R in charge and discharge process at difference SOC point0Value is estimated
SOH is counted, and using improvement UPF algorithm come tracking prediction R0.For the accuracy of the mentioned model and method of the verifying present invention, using changing
R into UPF algorithm to battery in constant current charge-discharge0Carry out real-time tracking prediction, and the reality obtained with second-order model
Value compares, the results showed that improves UPF algorithm to R0Tracking prediction ability be better than UKF algorithm, by comparing same loop
The estimated accuracy of SOH, again demonstrates side proposed by the invention under number difference SOC point and identical SOC point difference cycle-index
The superiority of method.
Detailed description of the invention
Fig. 1 is dynamic lithium battery Order RC equivalent-circuit model.
Relational graph of the Fig. 2 between SOC and open-circuit voltage.
Fig. 3 is that SOC is 30%~40% charging pulse experimental result.
Fig. 4 is that SOC is 40%~30% discharge pulse experimental result.
Fig. 5 is HPPC charging experiment SOC-R0Curve matching.
Fig. 6 is HPPC discharge test SOC-R0Curve matching.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
In order to improve the calculating simplicity and estimated accuracy of electric automobile lithium battery (dynamic lithium battery) SOH, root of the present invention
According to ohmic internal resistance R in battery charge and discharge process0SOH is estimated with the situation of change by stages of remaining capacity SOC, and using improvement
Carry out tracking prediction R without mark particle filter (Unscented Particle Filter, UPF) algorithm0, i.e., on the basis of UPF algorithm
On, discard tradition method for resampling, carries out resampling using the particle selection algorithm of the weight sequence survival of the fittest.Through the invention
The sample degeneracy problem that UPF algorithm can be improved, can effectively improve the estimated accuracy of dynamic lithium battery SOH.
Step of the invention is as follows:
Step 1: establish dynamic lithium battery Order RC equivalent model:
Equivalent-circuit model common at present mainly has RINT model, PNGV model, Thevenin model etc..Wherein,
RINT model belongs to ideal model, and there is no consider influence of the extraneous factor to battery；The internal resistance parameter and battery of PNGV model
DC internal resistance is more accurate, but the performance under different temperatures environment is unsatisfactory；Thevenin equivalent-circuit model energy ratio
The static and dynamic performance of preferable reaction cell, but original Thevenin model is limited to the responding ability of dynamic current.
The present invention is constructed on the basis of considering Thevenin model by increasing polarization resistance and polarization capacity in parallel
Order RC equivalent-circuit model, as shown in Figure 1.
Step 2: the parameter identification of Order RC battery model:
The constant current charge-discharge for carrying out at least 3 times to completely new lithium battery first is tested, if discharge capacity difference is 2% every time
Within, using this discharge capacity as battery rated capacity.
Then carry out 0.05C low current charge-discharge test fitting SOC and open-circuit voltage (Open Circuit Voltage,
OCV the relationship between) chooses 5%~95% SOC point (being spaced 5% between each SOC point) and its correspondence in experimentation
Charge and discharge stage voltage value, the charging/discharging voltage of each SOC point is averaged as the open-circuit voltage at each SOC point
Value.
It is distinguished followed by HPPC (Hybrid Pulse Power Characterization) pulse charge-discharge test
Know the parameter in second-order model, is 5%, 10%, 20% ... 90% in SOC, carries out HPPC charging respectively on 95% 11 points
Pulse test and the experiment of HPPC discharge pulse.Charging pulse is 0.2C electric current charging 30min, stands 90min；Discharge pulse is
0.2C current discharge 30min stands 90min, to the parameter identification of R0 and RC link, needs to combine in pulse charge-discharge test
Battery each SOC point response process, it is hereby achieved that HPPC charge and discharge Process Model Identification parameter.
Step 3: it is calculated based on the R0 for improving UPF method:
Calculating process is specific as follows:
(1) enabling k=0, population N is 50, init state variable x0=[Up1Up2R0SOC]TWith covariance matrix P and set
Initial weightI=1,2 ..., N, wherein Up1And Up2For the polarizing voltage of polarizing voltage RC link, ohmic internal resistance
R0In addition state variable of the SOC as system.
(2) weight prediction, the sampling of particle: k=k+1 carries out forecast updating to particle using UKF algorithm and calculates
Sigma point,To suggest that distribution function isExtract particle
(3) N number of particle weight is calculatedAnd carry out weight normalization
(4) judge whether resampling: calculating effective particle numberGiven threshold Nth=2N/3, if Neff
> Nth, (5) are gone to, otherwise to particle collectionResampling.
First particle weights during resamplingAccording to sequence sequence from small to large, while handle corresponds
Particle estimated valueWith true valueDifferenceIt is ranked up.ThenDivided byIt obtainsIt willAccording to sequence sequence from small to large, if the result being divided byIn have two groups of equal data, then retain grain
Sub- weightBiggish one group.?The 2/5 of average value is used as threshold value ωthIt selects the superior and eliminates the inferior, if a certain momentGreater than ωth, then corresponding particle weightsWith particle estimated valueIt remains.Grain after resampling
Subset isWherein
(5) optimal estimation of k moment state is calculatedAnd the corresponding variance matrix of each particle
(6) if k reaches the number of iterations of setting, algorithm terminates, and otherwise goes to (2).
The present invention realizes SOC and R using UPF algorithm is improved0State estimation, because of R0Change with the variation of SOC, estimates
Internal resistance R can be constantly corrected while counting SOC0Value improves UPF algorithm to R to further improve0Tracking prediction effect.
Estimation is participated in by improving the small particle of UPF algorithm picks difference comparsion, improves sample degeneracy and Loss of diversity problem, from
And make internal resistance R0Estimation it is more accurate.
Step 4: dynamic lithium battery SOH is calculated:
The health status of lithium battery refers to the current capacity of lithium battery, i.e., under certain condition, lithium battery can fill
Enter or release the percentage of electricity Yu battery nominal capacity.For pure electric automobile, when power battery SOH drops to 80%
When, it is believed that the end-of-life of battery.Since internal resistance can increase as battery SOH is reduced in the use process of battery, so
The health status of battery can be defined from the angle that the internal resistance of cell changes, it may be assumed that
RNOWIt (SOC) is the current internal resistance value of battery, RNEW(SOC) be battery factory when internal resistance value, REOLIt (SOC) is battery
Internal resistance value when end-of-life, RNOW(SOC) estimated to obtain by step 3), that is, optimal particle estimated value
In R0。RNEW(SOC) and REOL(SOC) it is solved and is obtained by linear equation in two unknowns group, because after carrying out cycle charge-discharge experiment
The practical charge/discharge capacity of battery can be obtained, so that practical SOH is obtained, the practical R during cycle charge-discharge0It again can be from
SOC and R0It is obtained in the matched curve of relationship, according to the calculating formula of SOH value, with RNEW(SOC) and REOL(SOC) unknown for two
Number, only it is to be understood that practical SOH and practical R under different charge and discharge number0Value (corresponding RNOW(SOC)), so that it may which it is primary that column write binary
Equation group obtains RNEW(SOC) and REOL(SOC), the current internal resistance value R of battery is finally estimated that by step 3)NOW(SOC),
And then cell health state (SOH) can be estimated according to formula (1).
Embodiment:
The embodiment that the present embodiment provides summary of the invention is divided into four steps.
1, model is established
As shown in Figure 1, according to Kirchoff s voltage, current law, the circuit equation of model are as follows:
U=UOCV(SOC)-IR0-Up1-Up2 (2)
Wherein R0For ohmic internal resistance, UOCVIt (SOC) is open-circuit voltage, there are non-linear relation, R with SOCpi(i=1,2) it is
Polarization resistance, Cpi(i=1,2) is polarization capacity, Upi(i=1,2) is polarizing voltage, and I is the charging and discharging currents of battery, and U is end
Voltage.
Choose the voltage U of two RC linksp1And Up2, ohmic internal resistance R0In addition state variable of the SOC as system chooses electricity
The end voltage in pond is as observed quantity, shown in structural regime equation such as formula (4), shown in observational equation such as formula (5).
Wherein, T is the sampling time；QNFor battery rated capacity；IkFor charging and discharging currents；wi,k~(0, Qi) it is system mode
Noise, wherein i=1,2,3,4；vk~(0, R) is systematic observation noise.
2, parameter identification
It is tested first by volume test, the rated capacity for demarcating lithium battery is 3.2Ah.Secondly low current charge and discharge are carried out
Relationship between experimental fit SOC and open-circuit voltage (Open Circuit Voltage, OCV), as shown in Figure 2.
The polynomial fitting of SOC and OCV relationship can be obtained by testing above are as follows:
Then HPPC (Hybrid Pulse Power Characterization) pulse charge-discharge test is carried out to recognize
Parameter in second-order model.To R0And the parameter identification of RC link, it needs in conjunction with battery in pulse charge-discharge test each
The response process of SOC point.
It is individually below 30% charging 30 minutes to SOC with SOC is that be 40% electric discharge 30 minutes be 40%, SOC to SOC
For 30%, internal resistance R is sketched0In the identification process in HPPC charge and discharge stage.
Then HPPC (Hybrid Pulse Power Characterization) pulse charge-discharge test is carried out to recognize
Parameter in second-order model.To R0And the parameter identification of RC link, it needs in conjunction with battery in pulse charge-discharge test each
The response process of SOC point.
It is individually below 30% charging 30 minutes to SOC with SOC is that be 40% electric discharge 30 minutes be 40%, SOC to SOC
For 30%, internal resistance R is sketched0In the identification process in HPPC charge and discharge stage.
As shown in figure 3, the voltage that this stage A → B, which is SOC, to be mutated when being 30%, C → D and E → F are SOC when being 40%
The voltage of mutation, B → C are 30 minutes constant-current charge virtual voltage curves, and the C moment starts to stand, and D → E is to stand 90 minutes in fact
Border voltage curve.
As shown in figure 4, the voltage that this stage A → B, which is SOC, to be mutated when being 40%, C → D and E → F are SOC when being 30%
The voltage of mutation, B → C are 30 minutes constant-current discharge virtual voltage curves, and the C moment starts to stand, and D → E is to stand 90 minutes in fact
Border voltage curve.
The section charge and discharge stage B → C is zero state response:
The section charge and discharge stage D → E is zero input response:
It takes battery charging and discharging to terminate moment and stands the average value of end transient voltage mutation to calculate R0, calculation formula
It is as follows:
According to the above method, successively identification obtains the ohmic internal resistance R of each SOC point in charge and discharge direction0, in the R of identification0Base
On plinth, SOC-R is obtained by polynomial fitting method0Relation curve, Fig. 5 are HPPC charging experiment SOC-R0Matched curve, Fig. 6 are
HPPC discharge test SOC-R0Matched curve.
HPPC charging experiment SOC-R0Six rank multinomial fitting formulas are as follows:
R0=0.0571 × SOC6-0.3695×SOC5+1.3369×SOC4-2.1067×SOC3
+1.5349×SOC2-0.5171×SOC+0.1244 (10)
HPPC discharge test SOC-R0Six rank multinomial fitting formulas are as follows:
R0=6.2941 × SOC6-21.1941×SOC5+28.2674×SOC4-18.9771×SOC3
+6.7523×SOC2-1.2194×SOC+0.1467 (11)
3, it is calculated based on the R0 for improving UPF method
Summary improves UPF algorithm calculating process by taking the 20th constant-current charge process as an example:
(1) enabling k=0, population N is 50, init state variable x0=[0 0 0.1262 0.04]TAnd covariance matrix
P=10-4×[1 0 0 0；0 1 00；0 0 10；000 1], if initial weight isI=1,2 ..., N.It adopts
Process noise matrix is Q=[10-9 10-10 10-11 10-8]T, observation noise matrix is R=0.001, sampling period T=
1s。
(2) weight prediction, the sampling of particle: k=k+1 carries out forecast updating to particle using UKF algorithm and calculates 9
Sigma point set and its weight coefficient.
(3) according to formulaN number of particle weight is calculated, and carries out weight normalizationS is the standard deviation that observational equation introduces noise.
(4) judge whether resampling: calculating effective particle numberGiven threshold Nth=2N/3, if Neff
> Nth, (5) are gone to, otherwise to particle collectionResampling.
First particle weights during resamplingAccording to sequence sequence from small to large, while handle corresponds
Particle estimated valueWith true valueDifferenceIt is ranked up.ThenDivided byIt obtainsIt willAccording to sequence sequence from small to large, if the result being divided byIn have two groups of equal data, then retain grain
Sub- weightBiggish one group.?The 2/5 of average value is used as threshold value ωthIt selects the superior and eliminates the inferior, if a certain momentGreater than ωth, then corresponding particle weightsWith particle estimated valueIt remains.Grain after resampling
Subset isWherein
(5) optimal estimation of k moment state is calculatedAnd the corresponding variance matrix of each particle
(6) if k reaches the sampling number 9400 of setting, algorithm terminates, and otherwise goes to (2).
Finally obtain state variable matrixThe practical SOH difference of 20th cycle charging and the 40th cycle charging experiment
For 0.9635 and 0.9549, when SOC is 0.5, the practical R of the 20th cycle charging and the experiment of the 40th cycle charging0Respectively
For 0.058 Ω and 0.0622 Ω, column write linear equation in two unknowns group:
R can be obtainedNEWIt (0.5) is 0.0402 Ω, REOL(0.5) it is 0.528 Ω, when SOC is 0.5, improves the of UPF estimation
The R of 20 cycle chargings0Value is 0.0589 Ω, to can estimate the SOH value of the 20th cycle charging according to formula (1).
4, dynamic lithium battery SOH is calculated
SOH=(0.528-0.0589)/(0.528-0.0402)=0.9617 can be obtained according to formula (1).
Illustrate effect of the invention to further choose, has chosen the different SOC points of the 20th constant current charge-discharge process
To compare the method for the present invention (improving UPF method) and not improve UKF method for the estimated accuracy of SOH, as shown in table 1.
The corresponding SOH value of the 20th constant current charge-discharge process difference SOC point of table 1
As shown in Table 1, either charging process still during discharge, this hair is estimated for the SOH of different SOC points
The method of bright proposition is superior to unmodified UPF method, to illustrate the superiority of improvement UPF algorithm.
Claims (9)
1. a kind of estimation method of electric automobile lithium battery SOH, which comprises the following steps:
1) the Order RC equivalent model of construction force lithium battery；
2) parameter identification of Order RC equivalent model, including ohmic internal resistance R0And Order RC link parameter；
3) ohmic internal resistance R is calculated using improvement UPF method0；
4) SOH value of dynamic lithium battery is calculated.
2. the estimation method of electric automobile lithium battery SOH according to claim 1 a kind of, which is characterized in that the step
It is rapid 1) in, Order RC equivalent model is that polarization resistance R in parallel is added based on Thevenin modelpiWith polarization capacity Cpi
Order RC equivalent-circuit model is constructed, and by the U of two RC linksp1And Up2, ohmic internal resistance R0With SOC as state variable x,
The end voltage U of battery constructs observational equation as observed quantity.
3. the estimation method of electric automobile lithium battery SOH according to claim 2 a kind of, which is characterized in that described two
The circuit equation of rank RC equivalent model are as follows:
U=UOCV(SOC)-IR0-Up1-Up2
Wherein, R0For ohmic internal resistance, UOCVIt (SOC) is the open-circuit voltage under SOC, I is the charging and discharging currents of battery, and U is end voltage,
Up1And Up2For the voltage of two RC links, UpiFor polarizing voltage.
4. the estimation method of electric automobile lithium battery SOH according to claim 3 a kind of, which is characterized in that the sight
Survey equation are as follows:
Wherein, T is sampling time, QNFor battery rated capacity, IkFor charging and discharging currents, w1,k、w2,k、w3,k、w4,kFor system mode
Noise, νkFor systematic observation noise, Cp1、Cp2The capacitor of respectively two RC links, Rp1、Rp2The electricity of respectively two RC links
Resistance, subscript k indicate sampling instant.
5. the estimation method of electric automobile lithium battery SOH according to claim 1 a kind of, which is characterized in that the step
It is rapid 2) in, the parameter in second-order model is recognized using HPPC pulse charge-discharge test, specifically:
21) it is tested by volume test, demarcates the rated capacity Q of lithium batteryN；
22) relational expression between low current charge-discharge test fitting acquisition SOC and open-circuit voltage OCV is carried out；
23) parameter in Order RC equivalent-circuit model is recognized by HPPC pulse charge-discharge test, obtains charge and discharge direction
The ohmic internal resistance of upper each SOC point, and be fitted to obtain the relational expression of SOC and ohmic internal resistance on charge and discharge direction.
6. the estimation method of electric automobile lithium battery SOH according to claim 4 a kind of, which is characterized in that the step
It is rapid 3) specifically includes the following steps:
31) sampling instant k=0 is enabled, total number of particles N, and init state variable x are set0, covariance matrix P and initial weight
32) sampling instant k=k+1 is enabled, particle is predicted using UKF algorithm, completes the estimation and sampling of weighting particle；
33) weight of all N number of particles is calculated separatelyAnd it is normalized to obtain normalization weight
34) number of effective particle is countedGiven threshold Nth=2N/3, if Neff> Nth, then carry out
Step 35), if Neff≤Nth, then to particle collectionCarry out resampling；
35) optimal estimation and the corresponding variance matrix of each particle of k sample moment state are obtained；
36) terminate when sampling number reaches setting value, otherwise return step 32).
7. the estimation method of electric automobile lithium battery SOH according to claim 6 a kind of, which is characterized in that the step
It is rapid 33) in, i-th of particle weightCalculating formula are as follows:
Wherein, S is the standard deviation that observational equation introduces noise, Upred_ukfFor the mean value of observation prediction.
8. the estimation method of electric automobile lithium battery SOH according to claim 7 a kind of, which is characterized in that the step
Rapid resampling 34) specifically:
341) to all normalization weightsIt is ranked up from small to large, and will be with normalization weightCorresponding particle is estimated
EvaluationWith true valueDifferenceIt is ranked up；
342) weight will be normalizedWith differenceQuotientIt is ranked up from small to large, if there are two groups in quotient
Equal data then delete more new sort after lesser one group of weight of normalization；
343) setting survival of the fittest threshold value ωthFor quotientThe 2/5 of average value, ifThen retain corresponding grain
Otherwise sub- weight and particle estimated value delete corresponding particle weight and particle estimated value, and more new particle collection.
9. the estimation method of electric automobile lithium battery SOH according to claim 1 a kind of, which is characterized in that the step
It is rapid 4) in, the calculating formula of the SOH value of dynamic lithium battery are as follows:
Wherein, REOLIt (SOC) is the internal resistance value at the end of battery life, RNOWIt (SOC) is the current internal resistance value of battery, i.e. optimal particle
Estimated valueIn ohmic internal resistance, RNEW(SOC) internal resistance value when dispatching from the factory for battery.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811056160.2A CN109061508A (en) | 2018-09-11 | 2018-09-11 | A kind of estimation method of electric automobile lithium battery SOH |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811056160.2A CN109061508A (en) | 2018-09-11 | 2018-09-11 | A kind of estimation method of electric automobile lithium battery SOH |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109061508A true CN109061508A (en) | 2018-12-21 |
Family
ID=64760201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811056160.2A Pending CN109061508A (en) | 2018-09-11 | 2018-09-11 | A kind of estimation method of electric automobile lithium battery SOH |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109061508A (en) |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109856561A (en) * | 2019-01-30 | 2019-06-07 | 北京长城华冠汽车科技股份有限公司 | A kind of health state evaluation method and apparatus of Vehicular dynamic battery group |
| CN110058160A (en) * | 2019-04-29 | 2019-07-26 | 西安理工大学 | The prediction technique of lithium battery health status based on SREKF |
| CN110133532A (en) * | 2019-05-29 | 2019-08-16 | 东北电力大学 | A kind of health status calculation method based on battery open circuit voltage |
| CN110320472A (en) * | 2019-05-17 | 2019-10-11 | 枣庄学院 | A kind of self-correction SOC estimation method for mining lithium battery |
| CN110398698A (en) * | 2019-07-30 | 2019-11-01 | 惠州市科达星辰技术有限公司 | A kind of method of battery management system SOH estimation |
| CN110456279A (en) * | 2019-08-15 | 2019-11-15 | 长安大学 | A kind of power battery cloud management system based on data-driven model |
| CN110456274A (en) * | 2019-08-29 | 2019-11-15 | 清华大学 | Battery impulse heats temperature rise rate estimation method |
| CN110806541A (en) * | 2019-11-12 | 2020-02-18 | 湖南大学 | AD-BAS-based lithium battery model parameter identification method |
| CN111090047A (en) * | 2019-12-09 | 2020-05-01 | 泉州装备制造研究所 | Lithium battery health state estimation method based on multi-model fusion |
| CN112595987A (en) * | 2020-11-28 | 2021-04-02 | 国网河南省电力公司电力科学研究院 | Lithium battery life estimation method based on mixed pulse voltage change |
| CN113075560A (en) * | 2021-04-19 | 2021-07-06 | 南京邮电大学 | Online estimation method for health state of power lithium ion battery |
| FR3106415A1 (en) | 2020-01-21 | 2021-07-23 | IFP Energies Nouvelles | Rapid and offline diagnostic procedures for accumulators and associated devices |
| CN113391212A (en) * | 2021-06-23 | 2021-09-14 | 山东大学 | Lithium ion battery equivalent circuit parameter online identification method and system |
| CN113552489A (en) * | 2021-06-29 | 2021-10-26 | 广东工业大学 | Power station standby power supply state of charge estimation method based on grey neural network |
| CN113557439A (en) * | 2019-03-12 | 2021-10-26 | 三星Sdi株式会社 | Method for estimating state of health of battery |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106354964A (en) * | 2016-09-08 | 2017-01-25 | 厦门理工学院 | Charge state estimation method of lithium-ion capacitor for electric automobile |
| CN106405433A (en) * | 2016-11-04 | 2017-02-15 | 首都师范大学 | Extended Kalman particle filtering based SOC (State Of Charge) estimation method and system |
| CN107367698A (en) * | 2017-08-25 | 2017-11-21 | 北京交通大学 | The health status Forecasting Methodology of electric automobile lithium battery group |
| CN107390127A (en) * | 2017-07-11 | 2017-11-24 | 欣旺达电动汽车电池有限公司 | A kind of SOC estimation method |
| US20180143257A1 (en) * | 2016-11-21 | 2018-05-24 | Battelle Energy Alliance, Llc | Systems and methods for estimation and prediction of battery health and performance |
| CN108445401A (en) * | 2018-02-09 | 2018-08-24 | 深圳市鹏诚新能源科技有限公司 | On-line Estimation method, electronic device and the storage medium of battery charge state SOC |
-
2018
- 2018-09-11 CN CN201811056160.2A patent/CN109061508A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106354964A (en) * | 2016-09-08 | 2017-01-25 | 厦门理工学院 | Charge state estimation method of lithium-ion capacitor for electric automobile |
| CN106405433A (en) * | 2016-11-04 | 2017-02-15 | 首都师范大学 | Extended Kalman particle filtering based SOC (State Of Charge) estimation method and system |
| US20180143257A1 (en) * | 2016-11-21 | 2018-05-24 | Battelle Energy Alliance, Llc | Systems and methods for estimation and prediction of battery health and performance |
| CN107390127A (en) * | 2017-07-11 | 2017-11-24 | 欣旺达电动汽车电池有限公司 | A kind of SOC estimation method |
| CN107367698A (en) * | 2017-08-25 | 2017-11-21 | 北京交通大学 | The health status Forecasting Methodology of electric automobile lithium battery group |
| CN108445401A (en) * | 2018-02-09 | 2018-08-24 | 深圳市鹏诚新能源科技有限公司 | On-line Estimation method, electronic device and the storage medium of battery charge state SOC |
Non-Patent Citations (3)
| Title |
|---|
| 刘力舟: "锂离子电池测试系统设计与实现", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 夏飞 等: "改进粒子滤波和小波包在汽轮机振动诊断中的应用", 《系统仿真学报》 * |
| 邓涛 等: "电动汽车动力电池 OH估计方法探讨", 《现代制造工程》 * |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109856561A (en) * | 2019-01-30 | 2019-06-07 | 北京长城华冠汽车科技股份有限公司 | A kind of health state evaluation method and apparatus of Vehicular dynamic battery group |
| JP7292404B2 (en) | 2019-03-12 | 2023-06-16 | 三星エスディアイ株式会社 | Method of estimating battery health |
| CN113557439A (en) * | 2019-03-12 | 2021-10-26 | 三星Sdi株式会社 | Method for estimating state of health of battery |
| JP2022524827A (en) * | 2019-03-12 | 2022-05-10 | 三星エスディアイ株式会社 | Battery health estimation method |
| EP3940405A4 (en) * | 2019-03-12 | 2022-12-21 | Samsung SDI Co., Ltd. | Method for estimating state of health of battery |
| CN110058160A (en) * | 2019-04-29 | 2019-07-26 | 西安理工大学 | The prediction technique of lithium battery health status based on SREKF |
| CN110320472A (en) * | 2019-05-17 | 2019-10-11 | 枣庄学院 | A kind of self-correction SOC estimation method for mining lithium battery |
| CN110320472B (en) * | 2019-05-17 | 2021-06-01 | 枣庄学院 | Self-correction SOC estimation method for mining lithium battery |
| CN110133532A (en) * | 2019-05-29 | 2019-08-16 | 东北电力大学 | A kind of health status calculation method based on battery open circuit voltage |
| CN110398698B (en) * | 2019-07-30 | 2021-05-07 | 惠州市科达星辰技术有限公司 | SOH estimation method for battery management system |
| CN110398698A (en) * | 2019-07-30 | 2019-11-01 | 惠州市科达星辰技术有限公司 | A kind of method of battery management system SOH estimation |
| CN110456279A (en) * | 2019-08-15 | 2019-11-15 | 长安大学 | A kind of power battery cloud management system based on data-driven model |
| CN110456274A (en) * | 2019-08-29 | 2019-11-15 | 清华大学 | Battery impulse heats temperature rise rate estimation method |
| CN110806541A (en) * | 2019-11-12 | 2020-02-18 | 湖南大学 | AD-BAS-based lithium battery model parameter identification method |
| CN111090047B (en) * | 2019-12-09 | 2022-01-28 | 泉州装备制造研究所 | Lithium battery health state estimation method based on multi-model fusion |
| CN111090047A (en) * | 2019-12-09 | 2020-05-01 | 泉州装备制造研究所 | Lithium battery health state estimation method based on multi-model fusion |
| WO2021148324A1 (en) | 2020-01-21 | 2021-07-29 | IFP Energies Nouvelles | Method for rapid offline diagnosis of accumulators, and related devices |
| FR3106415A1 (en) | 2020-01-21 | 2021-07-23 | IFP Energies Nouvelles | Rapid and offline diagnostic procedures for accumulators and associated devices |
| CN112595987A (en) * | 2020-11-28 | 2021-04-02 | 国网河南省电力公司电力科学研究院 | Lithium battery life estimation method based on mixed pulse voltage change |
| CN113075560A (en) * | 2021-04-19 | 2021-07-06 | 南京邮电大学 | Online estimation method for health state of power lithium ion battery |
| CN113075560B (en) * | 2021-04-19 | 2022-11-18 | 南京邮电大学 | Online estimation method for health state of power lithium ion battery |
| CN113391212A (en) * | 2021-06-23 | 2021-09-14 | 山东大学 | Lithium ion battery equivalent circuit parameter online identification method and system |
| CN113391212B (en) * | 2021-06-23 | 2022-05-17 | 山东大学 | Lithium ion battery equivalent circuit parameter online identification method and system |
| CN113552489A (en) * | 2021-06-29 | 2021-10-26 | 广东工业大学 | Power station standby power supply state of charge estimation method based on grey neural network |
| CN113552489B (en) * | 2021-06-29 | 2022-03-11 | 广东工业大学 | Power station standby power supply state of charge estimation method based on grey neural network |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109061508A (en) | A kind of estimation method of electric automobile lithium battery SOH | |
| Chen et al. | A review of lithium-ion battery for electric vehicle applications and beyond | |
| CN106443474B (en) | A kind of electrokinetic cell system service life Decline traits quickly know method for distinguishing | |
| Zhou et al. | A low-complexity state of charge estimation method for series-connected lithium-ion battery pack used in electric vehicles | |
| CN103792495B (en) | Method for evaluating battery performance based on Delphi method and grey relation theory | |
| CN105807226B (en) | Lithium ion battery SOC Forecasting Methodologies and device based on equivalent-circuit model | |
| Zhang et al. | Aging performance characterization and state-of-health assessment of retired lithium-ion battery modules | |
| CN109358293B (en) | Lithium ion battery SOC estimation method based on IPF | |
| Shen et al. | Neural network-based residual capacity indicator for nickel-metal hydride batteries in electric vehicles | |
| Li et al. | A novel state estimation approach based on adaptive unscented Kalman filter for electric vehicles | |
| Elmahdi et al. | Fitting the OCV-SOC relationship of a battery lithium-ion using genetic algorithm method | |
| CN111308356A (en) | SOC estimation method with weighted ampere-hour integration | |
| CN110221221A (en) | Charge states of lithium ion battery and health status combined estimation method | |
| Pei et al. | An equivalent circuit model for lithium battery of electric vehicle considering self-healing characteristic | |
| CN106772104A (en) | A kind of electrokinetic cell SOC value evaluation method | |
| CN111064253A (en) | Battery health degree rapid evaluation method based on average discrete Frechet distance | |
| Ye et al. | Model-based state-of-charge estimation approach of the Lithium-ion battery using an improved adaptive particle filter | |
| Li et al. | The open-circuit voltage characteristic and state of charge estimation for lithium-ion batteries based on an improved estimation algorithm | |
| CN103683427A (en) | Improved storage battery pack charging system and SOC (State of Charge) estimation method thereof | |
| CN104300184A (en) | Self-adaptive-fuzzy reasoning system-based method for achieving charging and discharging of battery | |
| Wei et al. | State-of-charge estimation for lithium-ion batteries based on temperature-based fractional-order model and dual fractional-order kalman filter | |
| Singh et al. | State of charge and health estimation of batteries for electric vehicles applications: key issues and challenges | |
| Song et al. | Series-connected lithium-ion battery pack health modeling with cell inconsistency evaluation | |
| CN106004481A (en) | SOH value estimation method for battery pack of hybrid electric vehicle | |
| Deng et al. | An improved algorithm of SOC testing based on open-circuit voltage-ampere hour method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181221 |