CN111830418A - SOC estimation method considering external environment influence of soft package battery - Google Patents
SOC estimation method considering external environment influence of soft package battery Download PDFInfo
- Publication number
- CN111830418A CN111830418A CN202010655623.8A CN202010655623A CN111830418A CN 111830418 A CN111830418 A CN 111830418A CN 202010655623 A CN202010655623 A CN 202010655623A CN 111830418 A CN111830418 A CN 111830418A
- Authority
- CN
- China
- Prior art keywords
- battery
- voltage
- equation
- soft package
- package battery
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000002474 experimental method Methods 0.000 claims abstract description 50
- 238000012937 correction Methods 0.000 claims abstract description 13
- 230000003044 adaptive effect Effects 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000005070 sampling Methods 0.000 claims description 13
- 230000010287 polarization Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 238000012546 transfer Methods 0.000 claims description 12
- 238000007599 discharging Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 229910000831 Steel Inorganic materials 0.000 claims description 8
- 239000010959 steel Substances 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 229910052751 metal Inorganic materials 0.000 claims description 5
- 239000002184 metal Substances 0.000 claims description 5
- 230000000052 comparative effect Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 239000004033 plastic Substances 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 3
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 5
- 229910052744 lithium Inorganic materials 0.000 description 5
- 229910052782 aluminium Inorganic materials 0.000 description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000000178 monomer Substances 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000005022 packaging material Substances 0.000 description 1
- 239000002985 plastic film Substances 0.000 description 1
- 229920006255 plastic film Polymers 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/385—Arrangements for measuring battery or accumulator variables
- G01R31/387—Determining ampere-hour charge capacity or SoC
- G01R31/388—Determining ampere-hour charge capacity or SoC involving voltage measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/389—Measuring internal impedance, internal conductance or related variables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/396—Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Secondary Cells (AREA)
Abstract
The invention discloses an SOC estimation method considering the external environment influence of a soft package battery, which comprises the following steps: establishing a state equation and acquiring parameters; carrying out a battery capacity calibration experiment, a constant current pulse experiment and an urban road running condition experiment; obtaining the battery capacity under a single influence factor by using data recorded in a battery capacity calibration experiment, and obtaining a function equation by using battery capacity calibration experiment data fitting; obtaining open-circuit voltage values under different conditions by using constant-current pulse experimental data, fitting power functions under different conditions, and further fitting the coefficient of the same order to influence factors; and substituting the obtained correction of the battery capacity and the battery open-circuit voltage based on the external environment into the adaptive unscented Kalman filtering to obtain an estimation result. The invention increases the influence factor of external mechanical load on the basis of temperature and discharge rate, can be closer to the real scene of the soft package battery in real vehicle application, and simultaneously enables the charge state estimation result of the battery to be more accurate.
Description
Technical Field
The invention belongs to the technical field of power battery management, and particularly relates to an SOC estimation method considering the external environment influence of a soft package battery.
Background
With the rapid development of technology and the strategic national demand for energy structure transformation, new energy automobiles have attracted extensive attention as an important means for alleviating energy crisis and environmental pollution. The most mature products of new energy automobiles are mainly electric automobiles, and the power source of the new energy automobiles is that battery monomers are combined into battery packs in a series-parallel connection mode, and then a plurality of battery packs are combined with hardware required by various systems to form a battery pack. In order to ensure that the power battery supplies power to the whole vehicle in a good state, a Battery Management System (BMS) needs to monitor the environment of each power battery and measure required data by using various sensors, including: charge-discharge current, terminal voltage, temperature, etc.; the possibility of occurrence of an emergency of the battery is greatly reduced, and obtaining an accurate SOC value of the battery is a prerequisite for ensuring that the BMS is in a good operation state.
The two most common methods for current SOC estimation are the ampere-hour integration method and the model-based estimation method, respectively. The ampere-hour integration method is based on an SOC (state of charge) definition formula, and under the condition that the total capacity and the charging and discharging coulombic efficiency of the battery are known, the SOC value of the battery at the current moment can be calculated only by acquiring the charging and discharging current value of the battery. However, this method can only obtain accurate results when the battery is in a constant working environment, the environment of the battery is constantly changed when the battery is in real-vehicle operation, and as the working time of the battery increases, the aging of the battery is inevitable, and various parameters are changed, so this method is not suitable for being used in a real-vehicle BMS system. The model-based estimation method is to finally acquire the SOC value by establishing an accurate battery model and combining different algorithms, and the algorithms can filter data noise and have certain robustness, so the method is more suitable for a dynamic working environment. Batteries can be classified into: cylinder lithium cell, square battery and laminate polymer battery. At present, many electric automobiles represented by Tesla MODEL S adopt a cylindrical lithium battery as a power source, and the main reasons are as follows: good consistency, low cost and good mechanical property of the monomer. However, the battery still has many problems, such as safety of the battery and high complexity of the battery system after the battery is assembled. The soft package battery has the advantages that the soft package battery has different military prominences in recent years, the market share is continuously improved, the soft package battery exceeds other two types of batteries, the packaging material and the structure of the soft package battery enable the soft package battery to have the advantages which are not possessed by a plurality of metal shell batteries, for example, the safety performance is good, the soft package battery is packaged by an aluminum plastic film structurally, and the soft package battery generally expands air and cracks when a safety problem occurs, and the soft package battery is not exploded like a steel shell battery or an aluminum shell battery; the weight of the soft package battery is 40% lighter than that of a steel shell lithium battery with the same capacity, and is 20% lighter than that of an aluminum shell lithium battery; the internal resistance is small, the internal resistance of the soft package battery is smaller than that of a lithium battery, and the self power consumption of the battery can be greatly reduced; the cycle performance is good, the cycle life of the soft package battery is longer, and the cycle attenuation of 100 times is 4-7% less than that of an aluminum shell; the design is flexible, the appearance can be changed into any shape, the thickness can be thinner, and the novel battery cell model can be developed according to the customization of the requirements of customers.
Compared with a metal shell battery, the soft package battery is carried in a real vehicle and is influenced by more factors, the SOC estimation of the battery is indirectly influenced except that the ambient temperature and the discharge rate influence the battery capacity and the open-circuit voltage, and the influence caused by external mechanical load after the soft package battery is grouped is not negligible. Therefore, the invention provides the SOC estimation method for correcting the environment where the soft package battery is located according to the point, so that a more accurate SOC value of the soft package battery in a dynamic environment can be obtained.
Disclosure of Invention
In view of the defects of the prior art, the invention aims to provide an SOC estimation method considering the external environment influence of a pouch battery, so as to solve the problems of large estimation error, poor robustness and the like caused by insufficient consideration factors for the state estimation of the pouch battery in the prior art.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
the invention discloses an SOC estimation method considering the external environment influence of a soft package battery, which comprises the following steps of:
step S1: establishing a state equation based on a second-order RC model, and acquiring parameters in the current second-order RC model by using an online parameter identification method in combination with the measured battery current and terminal voltage data;
step S2: selecting a soft package battery to be tested, selecting the test ranges of external mechanical load, temperature and discharge rate, and changing a single variable under a set standard environment to perform a corresponding battery capacity calibration experiment, a constant current pulse experiment and a general urban road driving condition experiment;
step S3: obtaining the battery capacity under a single influence factor by using data recorded in a battery capacity calibration experiment, introducing a capacity correction factor after normalization processing, defining the capacity correction factor as a function related to the three influence factors, and fitting the battery capacity calibration experiment data to obtain a function equation;
step S4: obtaining open-circuit voltage values under different conditions by using the constant-current pulse experimental data in the step S2, defining a fifth-order power function OCV-SOC by combining corresponding SOC values, fitting the power functions under different conditions, and further fitting the coefficient of the same-order term to influence factors to obtain corresponding results;
step S5: under the running conditions of urban roads under different external environments, the correction of the battery capacity and the battery open-circuit voltage based on the external environments, which is obtained by the experimental data, is substituted into the adaptive unscented Kalman filtering, so that the final estimation result is obtained.
Further, the laplace equation of the second-order RC model in step S1 is described as:
in the formula of UOCIs the open circuit voltage of the battery, R0Is the ohmic internal resistance of the battery,R1For electrochemical polarization resistance, R2For concentration polarization internal resistance, C1For electrochemical polarization of capacitance, C2Is a concentration polarization resistance, UtIs the terminal voltage of the battery, also called load voltage, I is the current flowing through the battery, i.e. load current;
the corresponding transfer equation is:
wherein, tau1=R1C1,τ2=R2C2;
The transfer equation is mapped from the S plane to the Z plane by a bilinear transformation method, the transformation equation isThe equation after conversion is summarized as follows:
then it is converted into a difference equation:
y(k)=Uoc(k)-Ut(k)=a1y(k-1)+a2y(k-2)+a3I(k)+a4I(k-1)+a5I(k-2)
in the formula, k, k-1, k-2 are respectively corresponding sampling time, and the open-circuit voltage in the differential equation is simplified as follows:
Uoc(k)=Uoc(k-1)=Uoc(k-2)
the difference equation can be finally collated as:
Ut(k)=(1-a1-a2)Uoc(k)+a1Ut(k-1)+a2Ut(k-2)-a3I(k)-a4I(k-1)-a5I(k-2)
defining a correlation vector:
θ(k)=[(1-a1-a2)Uoc(k) a1a2-a3-a4-a5]T
the final simplified difference equation is expressed as:
further, the online parameter identification method in step S1 adopts a least square method with forgetting factor, and its formula is as follows:
wherein,the parameter vector estimated at the previous moment, y (k) is the actual observed value of the system at the current moment, K (k) is a gain matrix, P (k) is a covariance matrix of theta (k), and lambda is a forgetting factor;
calculating to obtain theta (k) of each sampling moment by using a least square method with forgetting factors, and further solving to obtain final battery model parameters; order:
bringing the above formula into G (z)-1) The following equation is obtained:
and (3) the transfer function obtained by reverse deduction and the initial transfer function formula are equal according to the same term coefficient to obtain:
further derivation yields:
and the parameter value of the second-order RC model can be obtained in real time by using an online parameter identification method in combination with the current and the battery terminal voltage measured by the sensor when the soft package battery actually works.
Further, the step S2 specifically includes:
step S21: the method includes the steps that the influence of external mechanical load on the soft-package battery is researched, an experimental device is designed, steel pressure plates are arranged at the upper portion and the lower portion of the experimental device, the soft-package battery is placed in the middle of the pressure plates, and a layer of plastic pad is arranged between the soft-package battery and the steel pressure plates, so that the battery is prevented from being in direct contact with the pressure plates; the four corners of the two pressure plates are fastened by screws and nuts; selecting an external mechanical load range of 0-1500N, and selecting a group of experimental points every 300N, namely setting 6 groups of comparative experiments; before the experiment is started, a metal resistance wire of a pressure sensor is bonded on the surface of a soft package battery and protected by a covering layer, and a lead-out wire is connected to a display instrument, so that the monitoring of mechanical load applied to the outside of the soft package battery is realized; the soft package battery is charged to the upper cut-off voltage through constant current and constant voltage before the test, then stands for a period of time, and is discharged to the lower cut-off voltage through constant current and constant voltage under the set mechanical load to carry out the capacity calibration test; after each group of experiments is finished, the mechanical load is reset to zero, then the battery is recharged and stands still to prepare for next capacity calibration, and 0N and 25 ℃ are used as standard environments;
step S22: researching the influence of temperature on a soft package battery, selecting a corresponding temperature interval of-10-40 ℃, taking 25 ℃ as a reference and 5 ℃ as a temperature interval, respectively performing a capacity calibration experiment and a constant current pulse experiment, fully charging and standing the battery for a period of time in a standard environment before each group of experiments, then performing constant current and constant voltage discharge to a lower cut-off voltage at the current set temperature to realize capacity calibration, recovering the temperature to 25 ℃ and standing the battery for a period of time, fully charging and standing the battery again, then continuously discharging for 6 minutes at a discharge rate of 1C, standing the battery for 30 minutes, measuring the terminal voltage of the battery at the terminal moment of standing, taking the terminal voltage at the moment as a corresponding open-circuit voltage, and repeating the operation for 10 times, namely, finally emptying the battery; at this moment, capacity calibration and constant current pulse experiments at the same temperature are completed, and then the experiment temperature is changed through a constant temperature box to repeat the operations;
step S23: the method is characterized by researching the influence of discharge rate on the soft package battery, selecting corresponding discharge intervals of 0-5C, performing capacity experiments with 0.5C as an interval between 0C and 1C as an interval between 1C and 5C, fully charging the battery in a standard environment before each experiment, standing for a period of time, and then discharging to a lower cut-off voltage through constant current and constant voltage at a current set temperature to realize capacity calibration.
Further, the step S3 specifically includes: further processing is performed on each set of capacity values obtained in step S2, with the battery capacity at 0N, 25 ℃ and 1C discharge rate as denominator and the capacities measured under other different conditions as numerators, to perform normalization processing, and the obtained result is the capacity correction factor α described in step S3, which is defined as a function of three influencing factors, namely:
α=f(F T R)
in the formula, F is the external mechanical load of the battery, T is the working environment temperature of the battery, R is the discharge rate, and the specific form of the function is obtained by a fitting method.
Further, the step S4 specifically includes: processing the constant current pulse experimental data at the temperature changed in the step S22, and obtaining a power function formula by using the open-circuit voltage and the corresponding SOC value recorded at each group of temperatures in a fitting manner, wherein the specific form is as follows:
UOC=k5z5+k4z4+k3z3+k2z2+k1z+k0
in the formula, z is an SOC value, each term coefficient in the above formula is defined as a function of temperature T, and fitting is performed again by using the same term coefficient values obtained by fitting at different temperatures, that is:
ki=ai,4T4+ai,3T3+ai,2T2+ai,1T+ai,0
therefore, the accurate OCV-SOC relational expression can be obtained under different external environments.
Further, the specific steps of the Adaptive Unscented Kalman Filter (AUKF) in step S5 are as follows:
s51: obtaining 2n +1 sigma points and corresponding weight omega through one-time UT conversion, wherein n is the dimension of the state vector X;
in the formula, subscript i is the ith column of the square root of the matrix;is the state vector mean, P is the corresponding covariance matrix;
in the formula, subscript m is a mean value, c is a covariance, and superscript is a corresponding sampling point; λ ═ α2(n + k) -n is a scaling parameter, alpha is selected to control the distribution state of sampling points, k and beta are both parameters to be selected, and beta is more than or equal to 0;
establishing a state equation of the soft package battery based on a second-order RC model:
S52: a set of sample points is obtained by the UT transform formula of S51:
substituting the prediction into a state equation of S51 to obtain a one-step prediction of each point, wherein i is 1-2 n + 1;
X(i)(k+1|k)=f[k,X(i)(k|k)]
s53: calculating a one-step prediction and covariance matrix of the state vectors (U1, U2, z), and obtaining by weighted summation of prediction values of the sigma point set:
s54: and on the basis of the predicted value in the step S53, carrying out UT transformation again to generate a new sigma point set:
s55: substituting the point set into an observation equation to obtain a predicted observed quantity:
Z(i)(k+1|k)=h[X(i)(k+1|k)];
s56: and obtaining the predicted observed quantity of the sigma point set by the above steps, and obtaining the mean value and covariance of system prediction by weighted summation:
s57: calculating a gain matrix:
s58: updating the state vector x and the corresponding covariance using the gain matrix and the one-step predictor of the state vector:
s59: process noise Q by voltage innovationkAnd measuring the noise RkUpdating:
in the formula,iis the voltage innovation of the cell model at time k, i.e. the difference between the voltage measurement and the estimate, L is the window size, HkA covariance approximation of the voltage innovation at time k;
the last term z of the state vector x of step S58 is the resulting SOC value.
The invention has the beneficial effects that:
compared with a cylindrical battery, the soft package battery needs to consider more external environmental factors to influence various parameters of the battery, the external mechanical load influencing factor is increased on the basis of temperature and discharge rate, the soft package battery can be closer to a real scene when the soft package battery is applied in a real vehicle, and meanwhile, the estimation result of the charge state of the battery is more accurate.
Drawings
FIG. 1 is a diagram of a second order RC model.
FIG. 2 is a schematic diagram of an experimental apparatus.
Fig. 3 is an overall flow chart of the present invention.
Detailed Description
In order to facilitate understanding of those skilled in the art, the present invention will be further described with reference to the following examples and drawings, which are not intended to limit the present invention.
Referring to fig. 3, the SOC estimation method considering the external environment influence of the pouch battery according to the present invention includes the following steps:
step S1: establishing a state equation based on a second-order RC model, and acquiring parameters in the current second-order RC model by using an online parameter identification method in combination with the measured battery current and terminal voltage data;
referring to FIG. 1, the Laplace equation for the second order RC model is described as:
in the formula of UOCIs the open circuit voltage of the battery, R0Is the ohmic internal resistance, R, of the battery1For electrochemical polarization resistance, R2For concentration polarization internal resistance, C1For electrochemical polarization of capacitance, C2Is a concentration polarization resistance, UtIs the terminal voltage of the battery, also called load voltage, I is the current flowing through the battery, i.e. load current;
the corresponding transfer equation is:
wherein, tau1=R1C1,τ2=R2C2;
The transfer equation is mapped from the S plane to the Z plane by a bilinear transformation method, the transformation equation isThe equation after conversion is summarized as follows:
then it is converted into a difference equation:
y(k)=Uoc(k)-Ut(k)=a1y(k-1)+a2y(k-2)+a3I(k)+a4I(k-1)+a5I(k-2)
in the formula, k, k-1, k-2 are respectively corresponding sampling time, and the open-circuit voltage in the differential equation is simplified as follows:
Uoc(k)=Uoc(k-1)=Uoc(k-2)
the difference equation can be finally collated as:
Ut(k)=(1-a1-a2)Uoc(k)+a1Ut(k-1)+a2Ut(k-2)-a3I(k)-a4I(k-1)-a5I(k-2)
defining a correlation vector:
θ(k)=[(1-a1-a2)Uoc(k) a1a2-a3-a4-a5]T
the final simplified difference equation is expressed as:
the online parameter identification method adopts a least square method with forgetting factors, and the formula is as follows:
wherein,the parameter vector estimated at the previous moment, y (k) is the actual observed value of the system at the current moment, K (k) is a gain matrix, P (k) is a covariance matrix of theta (k), and lambda is a forgetting factor;
calculating to obtain theta (k) of each sampling moment by using a least square method with forgetting factors, and further solving to obtain final battery model parameters; order:
bringing the above formula into G (z)-1) The following equation is obtained:
and (3) the transfer function obtained by reverse deduction and the initial transfer function formula are equal according to the same term coefficient to obtain:
further derivation yields:
and the parameter value of the second-order RC model can be obtained in real time by using an online parameter identification method in combination with the current and the battery terminal voltage measured by the sensor when the soft package battery actually works.
Step S2: selecting a soft package battery to be tested, selecting the test ranges of external mechanical load, temperature and discharge rate, and changing a single variable under a set standard environment to perform a corresponding battery capacity calibration experiment, a constant current pulse experiment and a general urban road driving condition experiment;
step S21: the influence of external mechanical load on the soft-package battery is researched, an experimental device shown in figure 2 is designed, steel pressure plates are arranged above and below the experimental device, the soft-package battery is placed in the middle of the pressure plates, and a layer of plastic pad is arranged between the soft-package battery and the steel pressure plates, so that the battery is prevented from being in direct contact with the pressure plates; the four corners of the two pressure plates are fastened by screws and nuts; selecting an external mechanical load range of 0-1500N, and selecting a group of experimental points every 300N, namely setting 6 groups of comparative experiments; before the experiment is started, a metal resistance wire of a pressure sensor is bonded on the surface of a soft package battery and protected by a covering layer, and a lead-out wire is connected to a display instrument, so that the monitoring of mechanical load applied to the outside of the soft package battery is realized; the soft package battery is charged to the upper cut-off voltage through constant current and constant voltage before the test, then stands for a period of time, and is discharged to the lower cut-off voltage through constant current and constant voltage under the set mechanical load to carry out the capacity calibration test; after each group of experiments is finished, the mechanical load is reset to zero, then the battery is recharged and stands still to prepare for next capacity calibration, and 0N and 25 ℃ are used as standard environments;
step S22: researching the influence of temperature on a soft package battery, selecting a corresponding temperature interval of-10-40 ℃, taking 25 ℃ as a reference and 5 ℃ as a temperature interval, respectively performing a capacity calibration experiment and a constant current pulse experiment, fully charging and standing the battery for a period of time in a standard environment before each group of experiments, then discharging to a lower cut-off voltage through constant current and constant voltage at the current set temperature to realize capacity calibration, recovering the temperature to 25 ℃ and standing the battery for a period of time, fully charging and standing the battery again, then continuously discharging for 6 minutes (namely discharging for 10% each time) at a discharge rate of 1C, standing the battery for 30 minutes, measuring the terminal voltage of the battery at the terminal moment of standing, taking the terminal voltage as a corresponding open-circuit voltage, and repeating the operation for 10 times, namely obtaining the final emptying of the battery; at this moment, capacity calibration and constant current pulse experiments at the same temperature are completed, and then the experiment temperature is changed through a constant temperature box to repeat the operations;
step S23: the method is characterized by researching the influence of discharge rate on the soft package battery, selecting corresponding discharge intervals of 0-5C, performing capacity experiments with 0.5C as an interval between 0C and 1C as an interval between 1C and 5C, fully charging the battery in a standard environment before each experiment, standing for a period of time, and then discharging to a lower cut-off voltage through constant current and constant voltage at a current set temperature to realize capacity calibration.
Step S3: obtaining the battery capacity under a single influence factor by using data recorded in a battery capacity calibration experiment, introducing a capacity correction factor after normalization processing, defining the capacity correction factor as a function related to the three influence factors, and obtaining a function equation by using battery capacity calibration experiment data fitting;
further processing is performed on each set of capacity values obtained in step S2, with the battery capacity at 0N, 25 ℃ and 1C discharge rate as denominator and the capacities measured under other different conditions as numerators, to perform normalization processing, and the obtained result is the capacity correction factor α described in step S3, which is defined as a function of three influencing factors, namely:
α=f(F T R)
in the formula, F is the external mechanical load of the battery, T is the working environment temperature of the battery, R is the discharge rate, and the specific form of the function is obtained by a fitting method.
Step S4: obtaining open-circuit voltage values under different conditions by using the constant-current pulse experimental data in the step S2, defining a fifth-order power function OCV-SOC by combining corresponding SOC values, fitting the power functions under different conditions, and further fitting the coefficient of the same-order term to influence factors to obtain corresponding results;
the step S4 specifically includes: processing the constant current pulse experimental data at the temperature changed in the step S22, and obtaining a power function formula by using the open-circuit voltage and the corresponding SOC value recorded at each group of temperatures in a fitting manner, wherein the specific form is as follows:
UOC=k5z5+k4z4+k3z3+k2z2+k1z+k0
in the formula, z is an SOC value, each term coefficient in the above formula is defined as a function of temperature T, and fitting is performed again by using the same term coefficient values obtained by fitting at different temperatures, that is:
ki=ai,4T4+ai,3T3+ai,2T2+ai,1T+ai,0
therefore, the accurate OCV-SOC relational expression can be obtained under different external environments.
Step S5: under the urban road running working conditions under different external environments, the correction of the battery capacity and the battery Open Circuit Voltage (OCV) based on the external environment, which is obtained by the experimental data, is substituted into an Adaptive Unscented Kalman Filter (AUKF), so that a final estimation result is obtained;
the specific steps of the Adaptive Unscented Kalman Filtering (AUKF) are as follows:
s51: obtaining 2n +1 sigma points and corresponding weight omega through one-time UT conversion, wherein n is the dimension of the state vector X;
in the formula, subscript i is the ith column of the square root of the matrix;is the state vector mean, P is the corresponding covariance matrix;
in the formula, subscript m is a mean value, c is a covariance, and superscript is a corresponding sampling point; λ ═ α2(n + k) -n is a scaling parameter, alpha is selected to control the distribution state of sampling points, k and beta are both parameters to be selected, and beta is more than or equal to 0;
establishing a state equation of the soft package battery based on a second-order RC model:
S52: a set of sampling points (i.e., a sigma point set) is obtained by the UT transform formula of S51:
substituting the prediction into a state equation of S51 to obtain a one-step prediction of each point, wherein i is 1-2 n + 1;
X(i)(k+1|k)=f[k,X(i)(k|k)]
s53: calculating a one-step prediction and covariance matrix of the state vectors (U1, U2, z), and obtaining by weighted summation of prediction values of the sigma point set:
s54: and on the basis of the predicted value in the step S53, carrying out UT transformation again to generate a new sigma point set:
s55: substituting the point set into an observation equation to obtain a predicted observed quantity:
Z(i)(k+1|k)=h[X(i)(k+1|k)];
s56: and obtaining the predicted observed quantity of the sigma point set by the above steps, and obtaining the mean value and covariance of system prediction by weighted summation:
s57: calculating a gain matrix:
s58: updating the state vector x and the corresponding covariance using the gain matrix and the one-step predictor of the state vector:
s59: process noise Q by voltage innovationkAnd measuring the noise RkUpdating:
in the formula,iis the voltage innovation of the cell model at time k, i.e. the difference between the voltage measurement and the estimate, L is the window size, HkA covariance approximation of the voltage innovation at time k;
the last term z of the state vector x in step S58 is the finally obtained SOC value.
While the invention has been described in terms of its preferred embodiments, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention.
Claims (7)
1. An SOC estimation method considering external environment influence of a soft package battery is characterized by comprising the following steps:
step S1: establishing a state equation based on a second-order RC model, and acquiring parameters in the current second-order RC model by using an online parameter identification method in combination with the measured battery current and terminal voltage data;
step S2: selecting a soft package battery to be tested, selecting the test ranges of external mechanical load, temperature and discharge rate, and changing a single variable under a set standard environment to perform a corresponding battery capacity calibration experiment, a constant current pulse experiment and a general urban road driving condition experiment;
step S3: obtaining the battery capacity under a single influence factor by using data recorded in a battery capacity calibration experiment, introducing a capacity correction factor after normalization processing, defining the capacity correction factor as a function related to the three influence factors, and fitting the battery capacity calibration experiment data to obtain a function equation;
step S4: obtaining open-circuit voltage values under different conditions by using the constant-current pulse experimental data in the step S2, defining a fifth-order power function OCV-SOC by combining corresponding SOC values, fitting the power functions under different conditions, and further fitting the coefficient of the same-order term to influence factors to obtain corresponding results;
step S5: under the running conditions of urban roads under different external environments, the correction of the battery capacity and the battery open-circuit voltage based on the external environments, which is obtained by the experimental data, is substituted into the adaptive unscented Kalman filtering, so that the final estimation result is obtained.
2. The SOC estimation method considering the external environment influence of the soft-package battery according to claim 1, wherein the Laplace equation of the second-order RC model in the step S1 is described as follows:
in the formula of UOCIs the open circuit voltage of the battery, R0Is the ohmic internal resistance, R, of the battery1For electrochemical polarization resistance, R2For concentration polarization internal resistance, C1For electrochemical polarization of capacitance, C2Is a concentration polarization resistance, UtIs the terminal voltage of the battery, also called load voltage, I is the current flowing through the battery, i.e. load current;
the corresponding transfer equation is:
wherein, tau1=R1C1,τ2=R2C2;
The transfer equation is mapped from the S plane to the Z plane by a bilinear transformation method, the transformation equation isThe equation after conversion is summarized as follows:
then it is converted into a difference equation:
y(k)=Uoc(k)-Ut(k)=a1y(k-1)+a2y(k-2)+a3I(k)+a4I(k-1)+a5I(k-2)
in the formula, k, k-1, k-2 are respectively corresponding sampling time, and the open-circuit voltage in the differential equation is simplified as follows:
Uoc(k)=Uoc(k-1)=Uoc(k-2)
the difference equation can be finally collated as:
Ut(k)=(1-a1-a2)Uoc(k)+a1Ut(k-1)+a2Ut(k-2)-a3I(k)-a4I(k-1)-a5I(k-2)
defining a correlation vector:
θ(k)=[(1-a1-a2)Uoc(k) a1a2-a3-a4-a5]T
the final simplified difference equation is expressed as:
3. the SOC estimation method considering the external environment influence of the pouch battery as claimed in claim 1, wherein the online parameter identification method in step S1 adopts a least square method with forgetting factor, and the formula is as follows:
wherein,the parameter vector estimated at the previous moment, y (k) is the actual observed value of the system at the current moment, K (k) is a gain matrix, P (k) is a covariance matrix of theta (k), and lambda is a forgetting factor;
calculating to obtain theta (k) of each sampling moment by using a least square method with forgetting factors, and further solving to obtain final battery model parameters; order:
bringing the above formula into G (z)-1) The following equation is obtained:
and (3) the transfer function obtained by reverse deduction and the initial transfer function formula are equal according to the same term coefficient to obtain:
further derivation yields:
and the parameter value of the second-order RC model can be obtained in real time by using an online parameter identification method in combination with the current and the battery terminal voltage measured by the sensor when the soft package battery actually works.
4. The SOC estimation method considering the external environment influence of the soft package battery according to claim 1, wherein the step S2 specifically comprises:
step S21: the method includes the steps that the influence of external mechanical load on the soft-package battery is researched, an experimental device is designed, steel pressure plates are arranged at the upper portion and the lower portion of the experimental device, the soft-package battery is placed in the middle of the pressure plates, and a layer of plastic pad is arranged between the soft-package battery and the steel pressure plates, so that the battery is prevented from being in direct contact with the pressure plates; the four corners of the two pressure plates are fastened by screws and nuts; selecting an external mechanical load range of 0-1500N, and selecting a group of experimental points every 300N, namely setting 6 groups of comparative experiments; before the experiment is started, a metal resistance wire of a pressure sensor is bonded on the surface of a soft package battery and protected by a covering layer, and a lead-out wire is connected to a display instrument, so that the monitoring of mechanical load applied to the outside of the soft package battery is realized; the soft package battery is charged to the upper cut-off voltage through constant current and constant voltage before the test, then stands for a period of time, and is discharged to the lower cut-off voltage through constant current and constant voltage under the set mechanical load to carry out the capacity calibration test; after each group of experiments is finished, the mechanical load is reset to zero, then the battery is recharged and stands still to prepare for next capacity calibration, and 0N and 25 ℃ are used as standard environments;
step S22: researching the influence of temperature on a soft package battery, selecting a corresponding temperature interval of-10-40 ℃, taking 25 ℃ as a reference and 5 ℃ as a temperature interval, respectively performing a capacity calibration experiment and a constant current pulse experiment, fully charging and standing the battery for a period of time in a standard environment before each group of experiments, then performing constant current and constant voltage discharge to a lower cut-off voltage at the current set temperature to realize capacity calibration, recovering the temperature to 25 ℃ and standing the battery for a period of time, fully charging and standing the battery again, then continuously discharging for 6 minutes at a discharge rate of 1C, standing the battery for 30 minutes, measuring the terminal voltage of the battery at the terminal moment of standing, taking the terminal voltage at the moment as a corresponding open-circuit voltage, and repeating the operation for 10 times, namely, finally emptying the battery; at this moment, capacity calibration and constant current pulse experiments at the same temperature are completed, and then the experiment temperature is changed through a constant temperature box to repeat the operations;
step S23: the method is characterized by researching the influence of discharge rate on the soft package battery, selecting corresponding discharge intervals of 0-5C, performing capacity experiments with 0.5C as an interval between 0C and 1C as an interval between 1C and 5C, fully charging the battery in a standard environment before each experiment, standing for a period of time, and then discharging to a lower cut-off voltage through constant current and constant voltage at a current set temperature to realize capacity calibration.
5. The SOC estimation method considering the external environment influence of the soft package battery according to claim 4, wherein the step S3 specifically comprises: further processing is performed on each set of capacity values obtained in step S2, with the battery capacity at 0N, 25 ℃ and 1C discharge rate as denominator and the capacities measured under other different conditions as numerators, to perform normalization processing, and the obtained result is the capacity correction factor α described in step S3, which is defined as a function of three influencing factors, namely:
α=f(F T R)
in the formula, F is the external mechanical load of the battery, T is the working environment temperature of the battery, R is the discharge rate, and the specific form of the function is obtained by a fitting method.
6. The SOC estimation method considering the external environment influence of the soft package battery according to claim 4, wherein the step S4 specifically comprises: processing the constant current pulse experimental data at the temperature changed in the step S22, and obtaining a power function formula by using the open-circuit voltage and the corresponding SOC value recorded at each group of temperatures in a fitting manner, wherein the specific form is as follows:
UOC=k5z5+k4z4+k3z3+k2z2+k1z+k0
in the formula, z is an SOC value, each term coefficient in the above formula is defined as a function of temperature T, and fitting is performed again by using the same term coefficient values obtained by fitting at different temperatures, that is:
ki=ai,4T4+ai,3T3+ai,2T2+ai,1T+ai,0
therefore, the accurate OCV-SOC relational expression can be obtained under different external environments.
7. The SOC estimation method considering the external environment influence of the soft package battery according to claim 1, wherein the adaptive unscented Kalman filtering in the step S5 specifically comprises the following steps:
s51: obtaining 2n +1 sigma points and corresponding weight omega through one-time UT conversion, wherein n is the dimension of the state vector X;
in the formula, subscript i is the ith column of the square root of the matrix;is the state vector mean, P is the corresponding covariance matrix;
in the formula, subscript m is a mean value, c is a covariance, and superscript is a corresponding sampling point; λ ═ α2(n + k) -n is a scaling parameter, alpha is selected to control the distribution state of sampling points, k and beta are both parameters to be selected, and beta is more than or equal to 0;
establishing a state equation of the soft package battery based on a second-order RC model:
S52: a set of sample points is obtained by the UT transform formula of S51:
substituting the prediction into a state equation of S51 to obtain a one-step prediction of each point, wherein i is 1-2 n + 1;
X(i)(k+1|k)=f[k,X(i)(k|k)];
s53: calculating a one-step prediction and covariance matrix of the state vectors (U1, U2, z), and obtaining by weighted summation of prediction values of the sigma point set:
s54: and on the basis of the predicted value in the step S53, carrying out UT transformation again to generate a new sigma point set:
s55: substituting the point set into an observation equation to obtain a predicted observed quantity:
Z(i)(k+1|k)=h[X(i)(k+1|k)];
s56: and obtaining the predicted observed quantity of the sigma point set by the above steps, and obtaining the mean value and covariance of system prediction by weighted summation:
s57: calculating a gain matrix:
s58: updating the state vector x and the corresponding covariance using the gain matrix and the one-step predictor of the state vector:
s59: process noise Q by voltage innovationkAnd measuring the noise RkUpdating:
in the formula,iis the voltage innovation of the cell model at time k, i.e. the difference between the voltage measurement and the estimate, L is the window size, HkA covariance approximation of the voltage innovation at time k;
the last term z of the state vector x in step S58 is the finally obtained SOC value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010655623.8A CN111830418B (en) | 2020-07-09 | 2020-07-09 | SOC estimation method considering external environment influence of soft package battery |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010655623.8A CN111830418B (en) | 2020-07-09 | 2020-07-09 | SOC estimation method considering external environment influence of soft package battery |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111830418A true CN111830418A (en) | 2020-10-27 |
CN111830418B CN111830418B (en) | 2021-05-11 |
Family
ID=72900318
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010655623.8A Active CN111830418B (en) | 2020-07-09 | 2020-07-09 | SOC estimation method considering external environment influence of soft package battery |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111830418B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112051504A (en) * | 2020-08-13 | 2020-12-08 | 联合汽车电子有限公司 | Method and device for predicting battery capacity, terminal and computer-readable storage medium |
CN112630672A (en) * | 2020-12-31 | 2021-04-09 | 湖北亿纬动力有限公司 | Method for correcting capacity of lithium battery and application thereof |
CN114236401A (en) * | 2021-12-20 | 2022-03-25 | 上海正泰电源系统有限公司 | Battery state estimation method based on adaptive particle swarm optimization |
CN114487024A (en) * | 2021-12-31 | 2022-05-13 | 河南省日立信股份有限公司 | Calibration fitting method of palladium alloy hydrogen sensor based on power function |
CN116224099A (en) * | 2023-05-06 | 2023-06-06 | 力高(山东)新能源技术股份有限公司 | Method for dynamically and adaptively estimating battery SOC |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150226810A1 (en) * | 2013-11-12 | 2015-08-13 | Infineon Technologies Ag | Sensor arrangement, battery cell and energy system |
CN105759216A (en) * | 2016-02-26 | 2016-07-13 | 同济大学 | Method for estimating state of charge of soft package lithium-ion battery |
CN106054084A (en) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | Power battery SOC estimation method |
CN106716158A (en) * | 2014-06-11 | 2017-05-24 | 北京交通大学 | Method and device for estimating state of charge of battery |
CN106909716A (en) * | 2017-01-19 | 2017-06-30 | 东北电力大学 | The ferric phosphate lithium cell modeling of meter and capacity loss and SOC methods of estimation |
CN108872868A (en) * | 2018-06-12 | 2018-11-23 | 无锡智安能科技有限公司 | Detect integrating device, the system and method for battery pack state-of-charge and health status |
CN109560245A (en) * | 2018-10-18 | 2019-04-02 | 欣旺达电子股份有限公司 | Symmetrical cells and preparation method thereof |
CN109633472A (en) * | 2019-01-23 | 2019-04-16 | 刘平 | A kind of state-of-charge algorithm for estimating of lithium battery monomer |
CN111308351A (en) * | 2019-10-18 | 2020-06-19 | 南京航空航天大学 | Low-temperature environment power battery SOC estimation method, storage medium and equipment |
-
2020
- 2020-07-09 CN CN202010655623.8A patent/CN111830418B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150226810A1 (en) * | 2013-11-12 | 2015-08-13 | Infineon Technologies Ag | Sensor arrangement, battery cell and energy system |
CN106716158A (en) * | 2014-06-11 | 2017-05-24 | 北京交通大学 | Method and device for estimating state of charge of battery |
CN105759216A (en) * | 2016-02-26 | 2016-07-13 | 同济大学 | Method for estimating state of charge of soft package lithium-ion battery |
CN106054084A (en) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | Power battery SOC estimation method |
CN106909716A (en) * | 2017-01-19 | 2017-06-30 | 东北电力大学 | The ferric phosphate lithium cell modeling of meter and capacity loss and SOC methods of estimation |
CN108872868A (en) * | 2018-06-12 | 2018-11-23 | 无锡智安能科技有限公司 | Detect integrating device, the system and method for battery pack state-of-charge and health status |
CN109560245A (en) * | 2018-10-18 | 2019-04-02 | 欣旺达电子股份有限公司 | Symmetrical cells and preparation method thereof |
CN109633472A (en) * | 2019-01-23 | 2019-04-16 | 刘平 | A kind of state-of-charge algorithm for estimating of lithium battery monomer |
CN111308351A (en) * | 2019-10-18 | 2020-06-19 | 南京航空航天大学 | Low-temperature environment power battery SOC estimation method, storage medium and equipment |
Non-Patent Citations (3)
Title |
---|
ZHAO WANZHONG: ""Multi-objective Optimization of Differential Steering System of Electric Vehicle with Motorized Wheels"", 《TRANSACTIONS OF NANJING UNIVERSITY OF AERONAUTICS AND ASTRONAUTICS》 * |
赵万忠: ""基于无迹卡尔曼滤波的汽车状态参数估计"", 《华南理工大学学报》 * |
陈嵩: ""考虑迟滞效应的电池极化电压误差修正方法"", 《电网技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112051504A (en) * | 2020-08-13 | 2020-12-08 | 联合汽车电子有限公司 | Method and device for predicting battery capacity, terminal and computer-readable storage medium |
CN112630672A (en) * | 2020-12-31 | 2021-04-09 | 湖北亿纬动力有限公司 | Method for correcting capacity of lithium battery and application thereof |
CN112630672B (en) * | 2020-12-31 | 2023-10-03 | 湖北亿纬动力有限公司 | Lithium battery capacity correction method and application thereof |
CN114236401A (en) * | 2021-12-20 | 2022-03-25 | 上海正泰电源系统有限公司 | Battery state estimation method based on adaptive particle swarm optimization |
CN114236401B (en) * | 2021-12-20 | 2023-11-28 | 上海正泰电源系统有限公司 | Battery state estimation method based on self-adaptive particle swarm algorithm |
CN114487024A (en) * | 2021-12-31 | 2022-05-13 | 河南省日立信股份有限公司 | Calibration fitting method of palladium alloy hydrogen sensor based on power function |
CN114487024B (en) * | 2021-12-31 | 2023-11-03 | 河南省日立信股份有限公司 | Calibration fitting method of palladium alloy hydrogen sensor based on power function |
CN116224099A (en) * | 2023-05-06 | 2023-06-06 | 力高(山东)新能源技术股份有限公司 | Method for dynamically and adaptively estimating battery SOC |
Also Published As
Publication number | Publication date |
---|---|
CN111830418B (en) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111830418B (en) | SOC estimation method considering external environment influence of soft package battery | |
CN110501652B (en) | Rapid assessment method and device for available capacity of retired lithium battery | |
Wang et al. | Online dynamic equalization adjustment of high-power lithium-ion battery packs based on the state of balance estimation | |
CN102662148B (en) | On-line feedback battery state of charge (SOC) predicting method | |
CN107748336B (en) | On-line estimation method and system for state of charge of lithium ion battery | |
CN110196393B (en) | Combined on-line estimation method for lithium battery charge state, energy state and power state | |
CN105425154B (en) | A kind of method of the state-of-charge for the power battery pack for estimating electric automobile | |
CN107247235A (en) | A kind of battery capacity evaluation method for considering batteries in parallel connection difference | |
Ren et al. | Battery remaining discharge energy estimation based on prediction of future operating conditions | |
CN111352032A (en) | Lithium battery dynamic peak power prediction method | |
CN110109028A (en) | A kind of power battery remaining life indirect predictions method | |
CN103529399A (en) | Simulation method of lead acid battery-based improved PNGV model | |
CN110554321B (en) | Method for detecting SOC (state of charge) of retired power battery in real time | |
CN109239602B (en) | Method for estimating ohmic internal resistance of power battery | |
Li et al. | A new parameter estimation algorithm for an electrical analogue battery model | |
CN113253140A (en) | Battery health state online estimation method | |
CN105223512A (en) | Based on the method for the dynamic rectification dump energy of battery behavior | |
CN107015152A (en) | The discrimination method of battery pack inconsistency | |
CN112098849B (en) | Lithium battery residual electric quantity estimation method based on integral Kalman filtering | |
CN114740386A (en) | Lithium ion battery state-of-charge estimation method based on health state | |
CN111707953A (en) | Lithium battery SOC online estimation method based on backward smoothing filtering framework | |
CN115494400B (en) | Lithium battery lithium separation state online monitoring method based on ensemble learning | |
CN114200321A (en) | Lithium ion battery variable-order equivalent circuit model modeling method | |
CN103513188B (en) | The electricity computing method of battery cell in a kind of electric system energy storage station | |
CN117471324A (en) | Method and device for evaluating consistency of single capacity in battery module |
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 |