CN111929581B - Method for predicting internal and external temperatures of power lithium battery - Google Patents

Method for predicting internal and external temperatures of power lithium battery Download PDF

Info

Publication number
CN111929581B
CN111929581B CN202010505249.3A CN202010505249A CN111929581B CN 111929581 B CN111929581 B CN 111929581B CN 202010505249 A CN202010505249 A CN 202010505249A CN 111929581 B CN111929581 B CN 111929581B
Authority
CN
China
Prior art keywords
model
amb
equation
follows
lithium ion
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.)
Active
Application number
CN202010505249.3A
Other languages
Chinese (zh)
Other versions
CN111929581A (en
Inventor
庞辉
郭龙
武龙星
牟联晶
晋佳敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Technology
Original Assignee
Xian University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian University of Technology filed Critical Xian University of Technology
Priority to CN202010505249.3A priority Critical patent/CN111929581B/en
Publication of CN111929581A publication Critical patent/CN111929581A/en
Application granted granted Critical
Publication of CN111929581B publication Critical patent/CN111929581B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/374Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC] with means for correcting the measurement for temperature or ageing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Secondary Cells (AREA)

Abstract

The invention provides a method for predicting internal and external temperatures of a power lithium battery, which comprises the following steps of firstly, constructing a lithium ion battery electric-thermal coupling model considering environmental temperature; identifying the parameters of the model by a forgetting factor least square method to respectively obtain model parameters under different environments; then, carrying out simulation layer verification on the lithium ion battery electrical-thermal coupling model established by the parameter identification result by combining with dynamic working condition test data; and finally, constructing a state observer considering the environment temperature to carry out application level verification on the model. The invention solves the problems of large state estimation error and the like caused by inaccurate models due to different environmental temperatures, and plays an important role in the utilization efficiency, the service life and the whole vehicle performance of the power battery pack.

Description

Method for predicting internal and external temperatures of power lithium battery
Technical Field
The invention belongs to the technical field of lithium ion batteries, and particularly relates to a method for predicting internal and external temperatures of a power lithium battery.
Background
In recent years, lithium ion batteries have been widely used in the field of electric vehicles due to their high specific energy, high specific power, no memory effect, and long cycle life. The lithium ion power battery is a main energy supplier of a New energy automobile (NEV), has a plurality of advantages of light weight, low discharge rate, high energy density and the like, and is widely applied to the NEV. While a Battery Management System (BMS) is critical to the safe and efficient operation of lithium batteries (at a suitable voltage, operating temperature range) in electric vehicle applications, the System is usually connected to hundreds of Battery cells in a series/parallel configuration to meet the high power and high voltage requirements of the vehicle. In addition, thermal management is a key part of the BMS, and improper control of the temperature of the battery is a real threat to the traveling safety of the new energy automobile. Therefore, accurate modeling of lithium batteries is one of the most challenging tasks and difficulties to develop efficient battery management systems.
For this reason, many researchers at home and abroad have conducted a great deal of research on the mathematical model of the lithium ion battery. The method mainly comprises an electrochemical model and an equivalent circuit model. The equivalent circuit model has the characteristics of few parameters, convenient modeling and simple calculation in the aspect of calculating the external characteristics of the battery, so the equivalent circuit model is widely applied. As the internal and external temperatures of the lithium battery are difficult to directly measure, the electric-thermal coupling model is proposed by scholars for estimating the internal and external temperatures of the lithium battery and estimating the voltage, and the electric heating model has more parameters and needs to consider the influence of the environmental temperature on modeling.
It is noted that this section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
Disclosure of Invention
The invention aims to provide a method for predicting the internal and external temperatures of a power lithium battery, which solves the problems of large state estimation error and the like caused by inaccurate models due to different environmental temperatures in the prior art, and the provided modeling method can be controlled within a small error range in the estimation of the output internal and external temperatures.
In order to achieve the purpose, the invention adopts the following technical scheme:
the method for predicting the internal and external temperatures of the power lithium battery comprises the following steps:
s1: constructing an electric-thermal coupling model of the lithium ion battery considering the environmental temperature;
s2: identifying parameters of the lithium ion battery electrical thermal coupling model in the step S1 by a forgetting factor least square method to respectively obtain model parameters under different environments;
s3: carrying out simulation layer verification on the lithium ion battery electrical-thermal coupling model established by the parameter identification result in combination with dynamic working condition test data;
s4: and comparing the measured dynamic working condition experimental data with terminal voltage simulation data obtained by model simulation, performing application level verification on the model established in the step S1, and obtaining prediction of the internal and external temperatures of the lithium battery based on an Extended Kalman Filtering (EKF) method.
Further, the mathematical relationship of the electrical characteristics of the lithium ion battery electrical thermal coupling model in step S1 is as follows:
Figure GDA0003829982130000021
wherein, U t Is the battery terminal voltage; u shape OC Represents an open circuit voltage; r is ohmic internal resistance; r 1 (Tamb) and C 1 (T amb ) Electrochemical polarization resistance and electrochemical polarization fractional order capacitance, respectively, and is the ambient temperature T amb A function of (a); r is 2 (T amb ) And C 2 (T amb ) Concentration polarization resistance and concentration polarization fractional order capacitance respectively and is ambient temperature T amb A function of (a); i represents a load current; u shape 1 And U 2 Electrochemical polarization voltage and concentration polarization voltage are respectively represented.
Further, the mathematical relationship of the thermal characteristics of the lithium ion battery electrical and thermal coupling model in step S1 is as follows:
Figure GDA0003829982130000031
Figure GDA0003829982130000032
Q=(U OC -U t )I (4)
definition of
T is =T in -T amb (5)
T Ss =T S -T amb (6)
The above equations (2) and (3) are modified as follows:
Figure GDA0003829982130000033
Figure GDA0003829982130000034
the above equations (7) and (8) are modified as follows:
Figure GDA0003829982130000035
Figure GDA0003829982130000036
wherein, C C Is the cell core heat capacity; c S Is the thermal capacity of the housing; t is a unit of S Is the cell surface temperature; t is a unit of in Is the battery core temperature; q represents a calorific value; r is i And R o Respectively representing equivalent internal and external thermal resistances;
further, the application level verification in step S4 specifically includes:
discretization of equation (1) can yield:
Figure GDA0003829982130000041
discretization can be done for equation (9), where Δ t represents the sampling time:
Figure GDA0003829982130000042
discretization of equation (4) can yield:
Q(k)=I(k)(U OC (k)-U t (k)) (13)
formula (13) may be substituted for formula (12):
Figure GDA0003829982130000043
discretization of equation (10) can yield:
Figure GDA0003829982130000044
namely that
Figure GDA0003829982130000045
Definition of X (k) = [ SOC (k), U) 1 (k),U 2 (k),T is (k),T Ss (k)] T As a state variable of the system, Y (k) is an output voltage of the system, i.e., a terminal voltage; the state space equation of the system is set as follows:
Figure GDA0003829982130000046
wherein
g(X(k))=f OCV (SOC(k))-U 1 (k)-U 2 (k)-I(k)f R (T amb )(k) (18)
From the formulae (11), (14) and (16), the matrices A and B are:
Figure GDA0003829982130000051
definition of
X(k+1)=F(X(k),I(k),ω) (19)
Combining the model with an EKF algorithm, the method comprises the following specific steps:
step 1, initializing;
state initialization X 0 =E[X 0 ]
Covariance initialization P 0 =E[(X 0 -E(X 0 ))(X 0 -E(X 0 )) T ]
Noise initialization Q and R (omega and v are both Gaussian noise)
Step 2, predicting the state;
the state estimation value at the moment k can be calculated from the state value at the moment k-1:
Figure GDA0003829982130000052
in the formula (I), the compound is shown in the specification,
Figure GDA0003829982130000053
represents the prediction result at the time k;
the estimated value of the terminal voltage at the moment k is as follows:
Y k =g(X k )+R k (21)
the covariance estimate at time k is:
P k/k-1 =A k-1 P k-1 A k-1 T +Q k-1 (22)
step 3, updating and correcting;
according to the actual observed voltage y (k), updating the state value and covariance of the system, specifically as follows:
(1) Computing matrix H k I.e. jacobian matrix
Figure GDA0003829982130000061
Wherein the content of the first and second substances,
Figure GDA0003829982130000062
the derivative of the second order Taylor expansion of equation (4) with respect to SOC;
(2) Kalman gain calculation
K k =P k/k-1 H k T (H k P k/k-1 H k T +R k ) -1 (24)
(3) Measurement update
The state correction equation at the moment k is as follows:
Figure GDA0003829982130000063
the error covariance update equation at time k is:
P k/k =P k/k-1 -K k H k P k/k-1 (26)
the invention has the beneficial effects that:
1) The modeling method provided by the invention can be controlled within a smaller error range in the estimation of the output internal temperature and the output external temperature, so that the effectiveness and the accuracy of the parameter identification method are verified, the application of the electric vehicle in a wide temperature range is improved, and the problems of larger state estimation error and the like caused by inaccurate models due to different environmental temperatures are solved;
2) The method has great significance for state estimation and energy management applied to the battery management system of the electric vehicle; the method plays an important role in the utilization efficiency, the service life and the performance of the whole power battery pack.
Drawings
FIG. 1 is a schematic diagram of the present invention;
FIG. 2 is a schematic diagram of an electro-thermal coupling model according to the present invention;
FIG. 3 is a schematic diagram of a thermal model in an electrothermal coupling model according to the present invention;
FIG. 4 is a diagram illustrating the relationship between OCV-SOC-Tamb according to the present invention at different ambient temperatures;
FIG. 5 is a schematic diagram showing the relationship between the discharge capacity and the ambient temperature at different ambient temperatures according to the present invention;
FIG. 6 is a schematic diagram of the variation of the input current of the battery under the DST & US06 condition of the present invention;
FIG. 7 is a schematic diagram of the voltage variation of the battery terminal under the DST & US06 working condition of the present invention;
FIG. 8 is a schematic diagram of errors between a terminal voltage test value and a simulation value under a US06 working condition at an ambient temperature of-10 ℃ according to the present invention;
FIG. 9 is a schematic diagram showing a comparison between a terminal voltage test value and a simulation value under a US06 working condition at an ambient temperature of 50 ℃ in accordance with the present invention;
FIG. 10 is a schematic diagram of errors between a terminal voltage test value and a simulation value under US06 condition at an ambient temperature of 50 deg.C;
FIG. 11 is a schematic flow chart of the forgetting factor least square method in the present invention;
FIG. 12 is a schematic diagram showing the comparison between the external temperature test value and the simulation value under the working condition of US06 at the ambient temperature of-10 ℃ according to the present invention;
FIG. 13 is a schematic diagram showing the error of the comparison between the external temperature test value and the simulation value under the working condition of US06 at the ambient temperature of-10 deg.C;
FIG. 14 is a schematic diagram showing comparison between the external temperature test value and the simulation value under the working condition of US06 at the ambient temperature of 50 ℃ in accordance with the present invention;
FIG. 15 is a schematic diagram showing the error of the comparison between the external temperature test value and the simulation value under the working condition of US06 at the ambient temperature of 50 ℃ in accordance with the present invention;
FIG. 16 is a schematic diagram showing the comparison between the internal temperature test value and the simulation value under the working condition of US06 at the ambient temperature of-10 ℃ according to the present invention;
FIG. 17 is a schematic diagram showing the comparison error between the internal temperature test value and the simulation value under the working condition of US06 at the ambient temperature of-10 deg.C;
FIG. 18 is a schematic diagram showing the comparison between the internal temperature test value and the simulation value under the working condition of US06 at the ambient temperature of 50 ℃ in accordance with the present invention;
FIG. 19 is a schematic diagram showing the comparison error between the internal temperature test value and the simulation value under the working condition of US06 at the ambient temperature of 50 ℃.
Detailed Description
Example embodiments will now be described more fully with reference to the accompanying drawings. Example embodiments may, however, be embodied in many different forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of example embodiments to those skilled in the art. The described features or characteristics may be combined in any suitable manner in one or more embodiments.
As shown in FIG. 1, the invention discloses a method for predicting the internal and external temperatures of a power lithium battery, which specifically comprises the following steps:
step 1, constructing an electric-thermal coupling model of the lithium ion battery considering the environmental temperature;
step 2, identifying parameters of the battery electro-thermal coupling model; specifically, a least Forgetting factor method (FFLS) is adopted to identify model parameters, and model parameters under charging and discharging conditions at different environmental temperatures are obtained respectively;
step 3, carrying out simulation layer verification on the model established in the step 1 by combining with dynamic working condition test data;
and 4, comparing the measured dynamic working condition experimental data with terminal voltage simulation data obtained by model simulation, verifying the application level of the model established in the step 1, and predicting the internal and external temperatures of the lithium battery based on an EKF method.
The battery electro-thermal coupling model constructed in the step 1, wherein the mathematical relation of the electrical characteristic part is as follows:
Figure GDA0003829982130000081
wherein, U t Is the battery terminal voltage; u shape OC Represents an open circuit voltage; r is ohmic internal resistance; r 1 (Tamb) and C 1 (T amb ) Electrochemical polarization resistance and electrochemical polarization fractional order capacitance, respectively, and is the ambient temperature T amb Function of (2);R 2 (T amb ) And C 2 (T amb ) Concentration polarization resistance and concentration polarization fractional order capacitance respectively and is ambient temperature T amb A function of (a); i represents a load current; u shape 1 And U 2 Electrochemical polarization voltage and concentration polarization voltage are respectively represented.
Wherein, the mathematical relation of the thermal characteristic part is as follows: :
Figure GDA0003829982130000091
Figure GDA0003829982130000092
Q=(U OC -U t )I (4)
definition of
T is =T in -T amb (5)
T Ss =T S -T amb (6)
The above equations (2) and (3) are modified as follows:
Figure GDA0003829982130000093
Figure GDA0003829982130000094
the above equations (7) and (8) are modified as follows:
Figure GDA0003829982130000095
Figure GDA0003829982130000096
wherein, C C Is the cell core heat capacity; c S Is the thermal capacity of the housing; t is S Is the cell surface temperature; t is in Is the battery core temperature; q represents a calorific value; r i And R o Respectively, representing equivalent internal and external thermal resistances.
The application layer verification in the step S4 specifically includes:
discretization of equation (1) can yield:
Figure GDA0003829982130000101
discretization can be done for equation (9), where Δ t represents the sampling time:
Figure GDA0003829982130000102
discretization of equation (4) can yield:
Q(k)=I(k)(U OC (k)-U t (k)) (13)
formula (13) may be substituted for formula (12):
Figure GDA0003829982130000103
discretization of equation (10) can yield:
Figure GDA0003829982130000104
namely, it is
Figure GDA0003829982130000105
Definition of X (k) = [ SOC (k), U) 1 (k),U 2 (k),T is (k),T Ss (k)] T As a state variable of the system, Y (k) is an output voltage of the system, i.e., a terminal voltage; the state space equation of the system is set as follows:
Figure GDA0003829982130000106
wherein
g(X(k))=f OCV (SOC(k))-U 1 (k)-U 2 (k)-I(k)f R (T amb )(k) (18)
From the formulae (11), (14) and (16), the matrices A and B are:
Figure GDA0003829982130000111
definition of
X(k+1)=F(X(k),I(k),ω) (19)
Combining the model with an EKF algorithm, the method comprises the following specific steps:
step 1, initializing;
state initialization X 0 =E[X 0 ]
Covariance initialization P 0 =E[(X 0 -E(X 0 ))(X 0 -E(X 0 )) T ]
Noise initialization Q and R (both omega and nu are Gaussian noise)
Step 2, predicting the state;
the state estimation value at the moment k can be calculated from the state value at the moment k-1:
Figure GDA0003829982130000112
in the formula (I), the compound is shown in the specification,
Figure GDA0003829982130000113
represents the prediction result at the time k;
the estimated value of the terminal voltage at the moment k is as follows:
Y k =g(X k )+R k (21)
the covariance estimate at time k is:
P k/k-1 =A k-1 P k-1 A k-1 T +Q k-1 (22)
step 3, updating and correcting;
according to the actual observed voltage y (k), updating the state value and covariance of the system, specifically as follows:
(1) Computing matrix H k I.e. jacobi matrix
Figure GDA0003829982130000121
Wherein, the first and the second end of the pipe are connected with each other,
Figure GDA0003829982130000122
the derivative of the second order Taylor expansion of equation (4) with respect to SOC;
(2) Kalman gain calculation
K k =P k/k-1 H k T (H k P k/k-1 H k T +R k ) -1 (24)
(3) Measurement update
The state correction equation at the moment k is as follows:
Figure GDA0003829982130000123
the error covariance update equation at time k is:
P k/k =P k/k-1 -K k H k P k/k-1 (26)
in order to verify the accuracy of the estimation of the internal and external temperatures of the battery electrical-thermal coupling model established by the parameter identification result, firstly, an electrothermal coupling equivalent circuit model considering the influence of the environmental temperature is established in the MATLAB/Simulink environment. Secondly, identifying model parameters through an FFLS algorithm, verifying a simulation layer of the battery electrical and thermal coupling model established by the parameter identification result through a US06 working condition (US 06) and comparing US06 dynamic working condition experimental data obtained by experimental measurement with terminal voltage simulation data obtained by model simulation, wherein the results are respectively shown in a figure 7, a figure 8, a figure 10 and a figure 12, and errors are shown in a figure 9, a figure 11 and a figure 13. And finally, obtaining the estimation of the internal and external temperatures of the lithium battery based on an EKF method, and comparing the estimation with experimental data.
Through the experiment, the estimated root mean square errors of the terminal voltage under the working conditions of different environmental temperatures US06 can be obtained, wherein the estimated root mean square errors are 0.096 at minus 10 ℃, 0.128 at 20 ℃ and 0.067 at 50 ℃, and in addition, the estimated root mean square errors of the internal temperature and the external temperature are respectively 0.127,0.079 at minus 10 ℃; the internal and external temperature estimation errors at 20 ℃ were 0.121,0.075, respectively; the error in the estimation of the internal and external temperatures at 50 c was 0.063,0.039, respectively.
The data of simulation and application level show that the modeling method provided by the invention can be controlled in a smaller error range in the output internal and external temperature estimation, the effectiveness and the accuracy of the parameter identification method of the invention are verified, the application of the electric vehicle in a wide temperature range is improved, the problems of larger state estimation error and the like caused by inaccurate model due to different environmental temperatures are solved, and the method has great significance for the state estimation and the energy management of a battery management system of the electric vehicle; the device plays an important role in the utilization efficiency, the service life and the performance of the whole vehicle of the power battery pack.
Other embodiments of the invention will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the invention being indicated by the following claims.

Claims (1)

1. A method for predicting the internal and external temperatures of a power lithium battery is characterized by comprising the following steps:
s1: constructing an electric-thermal coupling model of the lithium ion battery considering the environmental temperature;
s2: identifying parameters of the lithium ion battery electrical thermal coupling model in the step S1 by a forgetting factor least square method to respectively obtain model parameters under different environments;
s3: carrying out simulation layer verification on the lithium ion battery electrical-thermal coupling model established by the parameter identification result in combination with dynamic working condition test data;
s4: comparing the measured dynamic working condition experimental data with terminal voltage simulation data obtained by model simulation, performing application level verification on the model established in the step S1, and obtaining prediction of the internal and external temperatures of the lithium battery based on an extended Kalman filtering method;
the mathematical relation of the electrical characteristic part of the lithium ion battery electrical thermal coupling model in the step S1 is as follows:
Figure FDA0003829982120000011
wherein, U t Is the battery terminal voltage; u shape OC Represents an open circuit voltage; r is ohmic internal resistance; r 1 (Tamb) and C 1 (T amb ) Electrochemical polarization resistance and electrochemical polarization fractional order capacitance, respectively, and is the ambient temperature T amb A function of (a); r is 2 (T amb ) And C 2 (T amb ) Concentration polarization resistance and concentration polarization fractional order capacitance respectively and is ambient temperature T amb A function of (a); i represents a load current; u shape 1 And U 2 Respectively representing electrochemical polarization voltage and concentration polarization voltage;
the mathematical relation of the thermal characteristic part of the lithium ion battery electrical thermal coupling model in the step S1 is as follows:
Figure FDA0003829982120000012
Figure FDA0003829982120000013
Q=(U OC -U t )I (4)
definition of
T is =T in -T amb (5)
T Ss =T S -T amb (6)
The above equations (2) and (3) are modified as follows:
Figure FDA0003829982120000021
Figure FDA0003829982120000022
the above equations (7) and (8) are modified as follows:
Figure FDA0003829982120000023
Figure FDA0003829982120000024
wherein, C C Is the cell core heat capacity; c S Is the thermal capacity of the housing; t is a unit of S Is the cell surface temperature; t is in Is the battery core temperature; q represents a calorific value; r i And R o Respectively representing equivalent internal and external thermal resistances;
the application level verification in the step S4 specifically includes:
discretization of equation (1) can yield:
Figure FDA0003829982120000025
discretization can be done for equation (9), where Δ t represents the sampling time:
Figure FDA0003829982120000026
discretization of equation (4) can yield:
Q(k)=I(k)(U OC (k)-U t (k)) (13)
formula (13) may be substituted for formula (12):
Figure FDA0003829982120000031
discretization of equation (10) can yield:
Figure FDA0003829982120000032
namely, it is
Figure FDA0003829982120000033
Define X (k) = [ SOC (k), U) 1 (k),U 2 (k),T is (k),T Ss (k)] T As a state variable of the system, Y (k) is an output voltage of the system, i.e., a terminal voltage; the state space equation of the system is set as follows:
Figure FDA0003829982120000034
wherein
g(X(k))=f OCV (SOC(k))-U 1 (k)-U 2 (k)-I(k)f R (T amb )(k) (18)
From the formulae (11), (14) and (16), the matrices a and B are:
Figure FDA0003829982120000035
definition of
X(k+1)=F(X(k),I(k),ω) (19)
Combining the model with an EKF algorithm, the method comprises the following specific steps:
step 1, initializing;
state initialization X 0 =E[X 0 ]
Covariance initialization P 0 =E[(X 0 -E(X 0 ))(X 0 -E(X 0 )) T ]
Noise initialization, wherein Q and R, omega and v are Gaussian noise
Step 2, predicting the state;
the state estimated value at the time k can be calculated from the state value at the time k-1:
Figure FDA0003829982120000041
in the formula (I), the compound is shown in the specification,
Figure FDA0003829982120000042
represents the prediction result at the time k;
the estimated value of the terminal voltage at the moment k is as follows:
Y k =g(X k )+R k (21)
the covariance estimate at time k is:
P k/k-1 =A k-1 P k-1 A k-1 T +Q k-1 (22)
step 3, updating and correcting;
according to the actual observed voltage y (k), updating the state value and covariance of the system, specifically as follows:
(1) Calculating the matrix H k I.e. jacobian matrix
Figure FDA0003829982120000043
Wherein the content of the first and second substances,
Figure FDA0003829982120000044
a derivative of the second order Taylor expansion of equation (4) with respect to SOC;
(2) Kalman gain calculation
K k =P k/k-1 H k T (H k P k/k-1 H k T +R k ) -1 (24)
(3) Measurement update
The state correction equation at the moment k is as follows:
Figure FDA0003829982120000045
the error covariance update equation at time k is:
P k/k =P k/k-1 -K k H k P k/k-1 (26)。
CN202010505249.3A 2020-06-05 2020-06-05 Method for predicting internal and external temperatures of power lithium battery Active CN111929581B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010505249.3A CN111929581B (en) 2020-06-05 2020-06-05 Method for predicting internal and external temperatures of power lithium battery

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010505249.3A CN111929581B (en) 2020-06-05 2020-06-05 Method for predicting internal and external temperatures of power lithium battery

Publications (2)

Publication Number Publication Date
CN111929581A CN111929581A (en) 2020-11-13
CN111929581B true CN111929581B (en) 2022-10-21

Family

ID=73316769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010505249.3A Active CN111929581B (en) 2020-06-05 2020-06-05 Method for predicting internal and external temperatures of power lithium battery

Country Status (1)

Country Link
CN (1) CN111929581B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112798971B (en) * 2020-12-30 2022-08-02 浙大城市学院 Soft-package type lithium ion battery coupling electric thermal model
CN113281655B (en) * 2021-05-20 2022-03-04 中南大学 Predictive control method and device for internal heating of power battery in low-temperature environment
CN113420471B (en) * 2021-06-01 2022-02-18 上海交通大学 Power lithium battery thermal model construction and establishment method and system based on electrochemical mechanism
CN113722926B (en) * 2021-09-07 2023-09-19 重庆大学 Square lithium battery electric heating coupling modeling error source analysis method
CN113884901A (en) * 2021-10-11 2022-01-04 西安交通大学 Battery surface temperature distribution estimation method and system
CN114184958B (en) * 2021-12-10 2024-04-12 中国华能集团清洁能源技术研究院有限公司 Lithium ion battery modeling method considering thermal characteristics
CN114757026B (en) * 2022-04-08 2024-06-07 昆明理工大学 Electrochemical coupling modeling method for full-working-condition multi-scale power lithium battery
CN116819344B (en) * 2023-08-25 2024-01-23 宁德时代新能源科技股份有限公司 Lithium battery nucleation overpotential prediction method, device, vehicle and medium

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015188610A1 (en) * 2014-06-11 2015-12-17 北京交通大学 Method and device for estimating state of charge of battery
WO2016134496A1 (en) * 2015-02-28 2016-09-01 北京交通大学 Method and apparatus for estimating state of charge of lithium ion battery
CN107741569A (en) * 2017-11-16 2018-02-27 温州大学 A kind of evaluation method of the lithium battery charge state based on segment extension Kalman filtering
CN108490356A (en) * 2018-03-06 2018-09-04 福建工程学院 A kind of lithium battery SOC estimation method for improving EKF algorithms
CN109669134A (en) * 2019-02-27 2019-04-23 浙江科技学院 A kind of evaluation method of the SOC based on Kalman filtering method
CN110361652A (en) * 2019-06-26 2019-10-22 河南理工大学 A kind of Kalman filtering lithium battery SOC estimation method based on Model Parameter Optimization
CN110412483A (en) * 2019-05-06 2019-11-05 西安理工大学 A kind of vehicle lithium battery state-of-charge estimation method based on interconnection observer
CN110515011A (en) * 2019-09-16 2019-11-29 华南农业大学 A kind of Accurate Estimation Method of lithium-ion-power cell SOC
CN110795851A (en) * 2019-10-30 2020-02-14 西安理工大学 Lithium ion battery modeling method considering environmental temperature influence
CN110954831A (en) * 2019-12-06 2020-04-03 重庆大学 Multi-time scale square lithium battery SOC and SOT joint estimation method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995464B (en) * 2014-05-26 2016-04-13 北京理工大学 A kind ofly estimate the parameter of the power system of electric vehicle and the method for state
EP2963434B1 (en) * 2014-06-30 2021-08-11 Foundation Of Soongsil University-Industry Cooperation Battery state estimation method and system using dual extended kalman filter, and recording medium for performing the method
GB2537406B (en) * 2015-04-16 2017-10-18 Oxis Energy Ltd Method and apparatus for determining the state of health and state of charge of lithium sulfur batteries
US10367235B2 (en) * 2016-02-19 2019-07-30 Cps Technology Holdings Llc Systems and methods for real-time parameter estimation of a rechargeable battery
FR3084170B1 (en) * 2018-07-17 2020-07-17 Zodiac Aero Electric METHOD FOR DETERMINING A STATE NOISE COVARIANCE MATRIX FOR THE ADJUSTMENT OF A BATTERY CHARGE OBSERVER AND DEVICE THEREOF

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015188610A1 (en) * 2014-06-11 2015-12-17 北京交通大学 Method and device for estimating state of charge of battery
WO2016134496A1 (en) * 2015-02-28 2016-09-01 北京交通大学 Method and apparatus for estimating state of charge of lithium ion battery
CN107741569A (en) * 2017-11-16 2018-02-27 温州大学 A kind of evaluation method of the lithium battery charge state based on segment extension Kalman filtering
CN108490356A (en) * 2018-03-06 2018-09-04 福建工程学院 A kind of lithium battery SOC estimation method for improving EKF algorithms
CN109669134A (en) * 2019-02-27 2019-04-23 浙江科技学院 A kind of evaluation method of the SOC based on Kalman filtering method
CN110412483A (en) * 2019-05-06 2019-11-05 西安理工大学 A kind of vehicle lithium battery state-of-charge estimation method based on interconnection observer
CN110361652A (en) * 2019-06-26 2019-10-22 河南理工大学 A kind of Kalman filtering lithium battery SOC estimation method based on Model Parameter Optimization
CN110515011A (en) * 2019-09-16 2019-11-29 华南农业大学 A kind of Accurate Estimation Method of lithium-ion-power cell SOC
CN110795851A (en) * 2019-10-30 2020-02-14 西安理工大学 Lithium ion battery modeling method considering environmental temperature influence
CN110954831A (en) * 2019-12-06 2020-04-03 重庆大学 Multi-time scale square lithium battery SOC and SOT joint estimation method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Parameter identification and state-of-charge estimation approach for enhanced lithium–ion battery equivalent circuit model considering influence of ambient temperatures;Hui Pang 等;《Chinese Physics B》;20191015;全文 *
Parameter identification and systematic validation of an enhanced single-particle model with aging degradation physics for Li-ion batteries;Hui pang 等;《Electrochimica Acta》;20190402;全文 *
一种基于简化电化学模型的锂电池互联状态观测器;庞辉 等;《物理学报》;20181106;全文 *
在线参数辨识和扩展卡尔曼算法的锂离子电池SOC估算研究;李博文等;《自动化仪表》;20200320(第03期);全文 *

Also Published As

Publication number Publication date
CN111929581A (en) 2020-11-13

Similar Documents

Publication Publication Date Title
CN111929581B (en) Method for predicting internal and external temperatures of power lithium battery
Wang et al. Model based insulation fault diagnosis for lithium-ion battery pack in electric vehicles
US9377512B2 (en) Battery state estimator combining electrochemical solid-state concentration model with empirical equivalent-circuit model
JP5683175B2 (en) An improved method for estimating the unmeasurable properties of electrochemical systems
CN102981125B (en) A kind of electrokinetic cell SOC method of estimation based on RC equivalent model
CN113049962B (en) LSTM-based energy storage device operation situation deduction method
Li et al. A combination state of charge estimation method for ternary polymer lithium battery considering temperature influence
CN111914503B (en) Method for establishing electric heating coupling model for power input of lithium ion battery
CN109900937B (en) Lithium battery charge state estimation method with temperature compensation function
Shi et al. A novel lumped thermal characteristic modeling strategy for the online adaptive temperature and parameter co-estimation of vehicle lithium-ion batteries
CN103675683A (en) Lithium battery state of charge (SOC) estimation method
CN110795851A (en) Lithium ion battery modeling method considering environmental temperature influence
Wang et al. Estimation of battery open-circuit voltage and state of charge based on dynamic matrix control-extended Kalman filter algorithm
CN104535934A (en) Online feed-forward compensating power battery charge state estimating method and system
CN112710955B (en) Algorithm for improving battery capacity estimation precision
Xia et al. Joint estimation of the state-of-energy and state-of-charge of lithium-ion batteries under a wide temperature range based on the fusion modeling and online parameter prediction
Luan et al. Estimation of state of charge for hybrid unmanned aerial vehicle Li-ion power battery for considering rapid temperature change
CN113093016A (en) Power battery management system fault diagnosis method based on uncertain noise filtering
CN115877232A (en) Lithium ion battery internal temperature estimation method based on Kalman filtering
Qin et al. Joint SOC–SOP estimation method for lithium-ion batteries based on electro-thermal model and multi-parameter constraints
Tian et al. Parallel-connected battery module modeling based on physical characteristics in multiple domains and heterogeneous characteristic analysis
CN112926186B (en) Online estimation method for internal temperature of single super capacitor
Liu et al. A novel Kalman-filter-based battery internal temperature estimation method based on an enhanced electro-thermal coupling model
CN106772075B (en) A kind of online battery impedance model optimization method considering temperature gradient
Jia et al. Temperature enhanced early detection of internal short circuits in lithium-ion batteries using an extended Kalman filter

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
GR01 Patent grant
GR01 Patent grant