US20140316728A1 - System and method for soc estimation of a battery - Google Patents
System and method for soc estimation of a battery Download PDFInfo
- Publication number
- US20140316728A1 US20140316728A1 US14/192,867 US201414192867A US2014316728A1 US 20140316728 A1 US20140316728 A1 US 20140316728A1 US 201414192867 A US201414192867 A US 201414192867A US 2014316728 A1 US2014316728 A1 US 2014316728A1
- Authority
- US
- United States
- Prior art keywords
- battery pack
- soc
- battery
- soc estimation
- discharging
- 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
Images
Classifications
-
- G01R31/3651—
-
- 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
-
- G01R31/3606—
-
- 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/374—Arrangements 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
-
- 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/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
-
- 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/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
- G01R31/3828—Arrangements for monitoring battery or accumulator variables, e.g. SoC using current integration
-
- 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/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
- G01R31/3835—Arrangements for monitoring battery or accumulator variables, e.g. SoC involving only voltage measurements
Definitions
- the present invention is related to Li-ion battery management systems, and more particularly to a system and method for state of charge (SOC) estimation of Li-ion batteries.
- SOC state of charge
- Li-ion battery due to its advantages of high specific energy and green environmental protection, have been widely used as large capacity power supply in various fields such as electronic-powered automotive, aerospace, ship gradually.
- energy density of li-ion battery becomes higher and higher, quantity of battery units in a battery pack also becomes larger and larger.
- asymmetry developed among the batteries in the battery pack can cause one or more of the batteries overcharging or over-discharging, and subsequently lowers the performance of the battery pack in the whole, resulting in serious effect on the service life of the battery pack. Therefore, a battery management system for managing and monitoring the working state of the battery pack is indispensable.
- state of charge is an important reference parameter of the working state of the li-ion battery pack, and is usually employed to indicate remainder energy of the li-ion battery pack.
- SOC state of charge
- Accurate SOC estimation of the li-ion battery pack utilized in automobiles can not only tell drivers of correct estimated mileage of the automobiles, but also ensure charging/discharging optimization of the li-ion battery pack, which ensures safe utility of the li-ion battery pack.
- large currents may cause the battery pack to be overly discharged and subsequently destroy the battery pack. Therefore, real time collection of voltage, temperature, and charging/discharging current of each battery is important for accurate SOC estimation of the battery so as to prolong the life of the battery pack and increase performance of the automobile.
- the SOC can be estimated based on attribute parameters such as voltage, current, resistance, temperature of the battery.
- attribute parameters of the battery generally can change in accordance with the aging of the battery and other uncertain factors, such as random road conditions the automobile is going through.
- ampere-hour method which is also a relatively accurate method on SOC estimation.
- the ampere-hour method employs real time current integral to calculate ampere hour, and then revises temperature, self-discharging data and ageing parameters that can affect the SOC estimation, and eventually obtains a relatively accurate SOC value by using a revision function and said parameters.
- the above-mentioned method is still far from being sufficient for practical situations because there are many other factors that could practically affect SOC estimation of the battery, and because it is hard to achieve the revision function in practice. Therefore, to date the SOC value estimated by employing the ampere-hour method can be far from the real SOC value of the battery.
- FIG. 1 is a schematic diagram of a battery system device including a battery management system according to embodiments of the present disclosure.
- FIG. 2 is a schematic diagram of the battery management system according to a preferred embodiment.
- FIG. 3 is a flowchart of a method for obtaining parameters of a second order RC equivalent circuit simulating the battery unit in accordance with an exemplary disclosure of the present invention.
- FIG. 4A is an exemplary second order RC equivalent circuit simulating the battery unit.
- FIG. 4B shows the curve line of voltage changes corresponding to the discharging currents with time elapsing.
- FIG. 4C shows a fitting curve representing values of U ed of FIG. 4A after the discharge current being removed.
- FIG. 4D shows a fitted voltage loaded upon the series circuit composed of the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization of FIG. 4A .
- FIG. 5 is a flow chart of SOC estimation method of the battery pack in use according to a preferred embodiment of the present disclosure.
- FIG. 6 is an exemplary OCV-SOC mapping table.
- a battery system device 100 such as an electric bicycle, an electric vehicle or an integrated power storage system, generally comprises a battery pack 101 , a battery management system (BMS) 102 for managing the battery pack 101 , and a load 103 powered by the battery pack 101 .
- the battery pack 101 may comprise only one single battery, or is composed of many batteries serially connected one by one. In the condition that the battery pack 101 comprises just one battery, the battery pack 101 can also be called battery 101 .
- a single battery is marked as a battery unit, therefore, the battery pack 101 may comprise one or more battery units.
- the battery management system 102 is used to manage and maintain the battery pack 101 , comprising but not limited to providing over-voltage and/or over-current protection, state of charge (SOC) estimation of the battery pack 101 .
- the battery pack 101 and the battery management system 102 collectively form a power system 10 of the battery system device 100 .
- the load 103 may be any kind of power consumption device, such as motors employed by the electric bicycle or the electric vehicle.
- the battery pack 101 or the battery unit described hereinafter is lithium ion (Li-ion) typed.
- FIG. 2 is a detailed description of the battery management system 102 .
- the battery management system 102 comprises a controller 1020 , a storage unit 1021 , a measurement module 1022 , a balance module 1023 , a power supply 1024 , a communication module 1025 , and a protection module 1026 .
- the measurement module 1022 , the communication module 1025 and the protection module 1026 are electronically connected to the controller 1020 by way of photoelectric coupling isolation circuits (PCIC), respectively.
- PCIC photoelectric coupling isolation circuits
- the storage unit 1021 may be a memory integrated with the controller 1020 , such as a flash memory, a static random access memory (SRAM), an electrically erasable programmable read-only memory (EEPROM).
- the storage unit 1021 may be those storage apparatuses independent from but electrically connected to the controller 1021 , such as a solid state disk or a micro hard disk.
- the storage unit 1021 may be the combination of the memory and the storage apparatus.
- the storage unit 1021 stores program codes that can be executed by the controller 1021 to maintain the battery pack 101 , for example, estimating the SOC of the battery pack 101 .
- the storage unit 1021 is also used to store data generated during the SOC estimation in accordance with the preferred embodiment of the present invention.
- the controller 1020 is the PIC18F458 type chip produced by Microchip Technology Incorporation.
- the PIC18F458 type chip is an 8-bit micro-controller with 32 Kilobytes memory space for storing program codes, 1536 bytes-sized SRAM and 256 bytes-sized EEPROM.
- the PIC18F458 type chip further comprises a clock with frequency up to 32 MegaHz, which can operate as an external crystal oscillator circuit for, such as timing.
- Another outstanding feature of the PIC18F458 type chip is that the power consumption of the chip is relatively low. For example, under a sleeping mode, current consumed by the chip is only 0.2 ⁇ A.
- the measurement module 1022 is an important part of the battery management system 102 , and is used for measuring parameters such as current, voltage, and temperature of the battery pack 101 .
- the measured parameters are not only basis for the SOC estimation, but security guarantee of the battery pack 101 .
- the current measurement comprises measuring the current flowing through the battery pack 101 .
- the voltage measurement comprises measuring the voltage of the battery pack 101 , and the voltage of each battery unit in the battery pack 101 , providing the battery pack 101 comprises more than one battery units.
- the temperature measurement comprises measuring the temperature of the battery units in the battery pack 101 .
- the measured current, voltage and temperature data are stored in the storage unit 1021 .
- the balance module 1023 is used for balancing the charging/discharging current of the battery pack 101 .
- the balancing module 1023 comprises monitoring chip LTC6803 and other accompanied circuits provided by Linear Technology Corporation.
- the power supply 1024 is used for providing various reference voltages to different function modules, such as the controller 1021 , the monitoring chip LTC6803, etc.
- the power supply 1024 employs the LM2009 chip to provide fixed voltage, such as 15V, and uses other DC/DC regulator, such as LM7805 chip, to convert the fixed 15V voltage to various DC outputs, such as 10V, 5V, and provides the converted DC outputs to the different function modules.
- the communication module 1025 provides communication services for above-mentioned function modules.
- the communication services comprise a serial peripheral interface (SPI), a controller area network (CAN), and serial communication.
- SPI serial peripheral interface
- CAN controller area network
- serial communication serial communication
- the protection module 1026 provides protections under the conditions such as over-current, over-voltage, under-voltage, to ensure normal operation of the battery system device 100 , especially the battery pack 101 .
- electrochemical parameters of the battery pack 101 should be obtained.
- the electrochemical performance of the battery unit in the battery pack 101 generally presents non-linear characteristic.
- a second order RC equivalent circuit is employed to simulate the battery unit, in that the second order RC equivalent circuit owns the benefits comprising, for example, the order is relatively low, and parameters thereof are easily to be obtained.
- the second order RC equivalent circuit can be easily constructed and preferably simulates dynamic characteristics of the battery unit under working status.
- the dynamic characteristics comprise an ohmic resistance R 0 , an electrochemical polarization resistance R e , an electrochemical polarization capacitor C e , a concentration polarization resistance R d and a concentration polarization capacitor C d .
- FIG. 3 is a flowchart of a method for obtaining parameters of the second order RC equivalent circuit simulating the battery unit in accordance with an exemplary disclosure of the present invention.
- the second order RC equivalent circuit simulating the battery unit is constructed, as illustrated in FIG. 4A .
- the constructed second order RC equivalent circuit comprises the ohmic resistance R 0 , the electrochemical polarization resistance R c , the electrochemical polarization capacitor C c , the concentration polarization resistance R d and the concentration polarization capacitor C d .
- the electrochemical polarization resistance R e , the electrochemical polarization capacitor C e are connected in parallel to form a capacitance-resistance segment of electrochemical polarization, and the concentration polarization resistance R d and the concentration polarization capacitor C d are connected in parallel to form a capacitance-resistance segment of concentration polarization.
- the ohmic resistance R 0 , the capacitance-resistance segment of electrochemical polarization, and the capacitance-resistance segment of concentration polarization are connected in series.
- the experimental method comprises the following steps:
- Step I keeping the battery unit under a status without charging or discharging over one hour
- Step II discharging the battery unit continuously lasting for 900 seconds
- Step III keeping the battery unit under the status without charging or discharging for a relatively long time, such as one hour, up to the end of the experiment.
- FIG. 4B shows the curve line of voltage changes corresponding to the discharging currents with time elapsing.
- the exemplary embodiment firstly, placing one battery unit under a status without charging or discharging over one hour; then loading a discharge current pulse on a positive electrode and a negative electrode of the battery unit, for example, at the time point of 1800 s shown in FIG. 4B .
- a discharge current pulse is loaded, a responsive voltage upon between the positive electrode and the negative electrode of the battery unit gets a downward mutation.
- the discharge current pulse is removed, and the responsive voltage upon the positive electrode and the negative electrode of the battery unit gets a upward mutation.
- the status without charging or discharging of the battery could be defined as an idle status of the battery unit.
- the voltage upon the capacitance-resistance segment of concentration polarization composed of the concentration polarization resistance R d and the concentration polarization capacitor C d could not get a mutation, due to the exist of the concentration polarization capacitor C d .
- the voltage upon the capacitance-resistance segment of electrochemical polarization composed of the electrochemical polarization resistance R e and the electrochemical polarization capacitor C e also could not get a mutation. Therefore, the downward mutation and the upward mutation can only be caused by the ohmic resistance R 0 .
- the discharge current pulse is removed, the voltage upon the ohmic resistance R 0 would become 0 due to no currents flowing through the ohmic resistance R 0 .
- the upward mutation of the voltage upon the positive electrode and the negative electrode of the battery unit is caused by electronic discharging of the capacitors C d and C e , which are charged when the battery unit is discharging.
- the capacitors C d and C e discharge completely, the voltage upon the positive electrode and the negative electrode of the battery unit becomes stable at, for example, about 4.05 v as shown in FIG. 4B .
- the voltage upon the positive electrode and the negative electrode of the battery unit is close to electromotive force of the battery unit.
- the voltage response upon the positive electrode and the negative electrode of the battery unit at the time point t could be deemed as zero state response:
- U d and U e respectively stand for voltages loaded upon the capacitors C d and C e
- U ed represents the voltage loaded upon the series circuit composed of the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization.
- FIG. 4C shows a fitting curve representing values of U ed after the discharge current being removed.
- ⁇ e and ⁇ d can be achieved by using the exponent curve fitting module in the MATLAB software.
- the voltage upon the positive electrode and the negative electrode of the battery unit declines due to the ohmic resistance R 0 , the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization. Because the battery unit has been place at the idle status for a long time period, such as over one hour, the capacitor C e in the capacitance-resistance segment of electrochemical polarization and the capacitor C d in the capacitance-resistance segment of concentration polarization could be deemed as zero electron state.
- supposing the time point t of loading the discharge current pulse is zero time
- voltage response of the capacitor C e and C d to the discharge current pulse is a zero state response:
- the voltage U ed loaded upon the series circuit composed of the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization is fitted as shown in FIG. 4D .
- step S 307 putting the calculated ⁇ e and ⁇ d into the formula 1-2, and by using the exponent curve fitting method of the MATLAB software, the R d and R e can be achieved.
- the second order RC equivalent circuit is used to simulate the battery unit. Therefore, the identified parameters of the second order RC equivalent circuit are indeed dynamic characteristic paramters of the battery unit. According to above-described identifying processes, in one embodiment of the present disclosure, identified dynamic characteristic parameters of one exemplary battery unit are shown in following table:
- the parameters of the second order RC equivalent circuit simulating the battery unit which are identified by way of said-mentioned method are stored in the storage unit 1021 , and would be invoked in estimating processed when the battery unit is used in practice.
- FIG. 5 is a flow chart of SOC estimation method of the battery pack 101 in use according to a preferred embodiment of the present disclosure, which could be applied to the battery system device 100 of FIG. 1 .
- the SOC estimation method could be realized by the controller 1020 executing the programmed codes stored in the storage unit 1021 .
- the SOC estimation method is periodically executed by the controller 1020 to achieve a relatively accurate SOC of the battery pack 101 in different time points.
- a time cycle the method being executed periodically could be, for example, 10 millisecond (ms). In other embodiment, the time period could be set or configured as other value, such 15 ms.
- an once-through operation of the SOC estimation method is described as one SOC estimation process.
- the clock in the controller 1020 would record the time duration of the battery pack 101 being in the idle state in which there is no charging or discharging occurred to the battery pack 101 .
- the time duration of the battery pack 101 in the idle state is simplified as an idle time of the battery pack 101 .
- some program codes would be executed to determine whether the idle time of the battery pack 101 is longer than a predefined time period, such as one hour, at step S 501 .
- the controller 1020 reads a previous SOC value (here marked as SOC (k-1) ) generated during a previous SOC estimation process, and regards the SOC (k-1) as an initial SOC value of a current SOC estimation process.
- a previous SOC value here marked as SOC (k-1)
- an open circuit voltage (OCV) of the battery pack 101 is substantially equal to an electromotive force of the battery pack 101 . Therefore, if the idle time of the battery pack 101 is longer than the predefined time period, at step S 502 B, the measurement module 1022 measures a current voltage upon the positive electrode and the negative electrode of the battery pack 101 , subsequently, the controller 1020 queries an OCV-SOC mapping table to determine a SOC value corresponding to the current voltage, and records the determined SOC value as SOC (k-1) . The SOC (k-1) would act as an initial SOC value of a current SOC estimation process.
- step S 502 A and step S 502 B are alternative in one SOC estimation process according to a preferred embodiment of the present invention.
- the OCV-SOC mapping table is established by an experimental method.
- the experiment method for establishing the OCV-SOC mapping table comprises the steps of:
- step D) repeating said step C), under constant temperature of substantially 20 degree Celsius, until the battery pack 101 discharges completely;
- the programmable electronic load employs Chroma 6310A type programmable DC electronic load provided by Chroma Company.
- FIG. 6 is an exemplary OCV-SOC mapping table of an 18650 type li-ion battery, which is obtained by using the experiment method for establishing the OCV-SOC mapping table.
- the controller 1020 determines whether the battery pack 101 is in a working status.
- the working status is opposite to the idle status, that is, the working status means the battery pack 101 is discharging and/or charging.
- the controller 1020 determines the batter pack 101 is in the working status if the battery pack 101 is discharging and/or charging. If the battery pack 101 is not in the working status, or is in the idle status, the flow returns to step S 501 to determine whether an idle time of the battery pack 101 is longer than the predefined time period in another SOC estimation process.
- the measurement module 1022 measures current I (k) flowing through and voltage U (k) upon the two polarities of the battery pack 101 .
- a direction of the measured current I (k) shows the discharging or charging status of the battery pack 101 .
- the current I (k) is positive, and if the battery pack 101 is charging, the current I (k) is negative.
- an ampere hour method is employed to estimate SOC (k) in current SOC estimation process, based on the determined SOC (k-1) in step S 502 A or S 502 B.
- the ampere hour method is also called current integration method, which is a fundamental method to measure the SOC of batteries.
- current SOC estimation process after determining the SOC (k-1) in step S 502 A or S 502 B, the battery pack 101 is charging or charging from time point k-1 to k, such as 10 ms, a change of the SOC of the battery pack 101 could be represented as:
- ⁇ ⁇ ⁇ SOC 1 Q 0 ⁇ I ⁇ ( k ) ⁇ T ( 2 ⁇ - ⁇ 1 )
- Q 0 is a total quantity of electric discharge in fixed current of 0.1*C of the battery pack 101 , wherein C represents the nominal capacity of the battery pack 101 .
- I (k) represents the charging or discharging current of the battery pack 101 , wherein I (k) is positive if the battery pack 101 is discharged, otherwise the I (k) is negative.
- Peukert proposes an empirical formula to revise the SOC of the battery pack 101 working with changing currents, as shown in following formula 2-2:
- I discharging current
- t is a discharging time length
- n and K are constants that are determined by types and active materials of the battery pack 101 .
- the active material is lithium.
- k ⁇ is a revise value to the capacity volume of the battery pack 101 .
- SOC* (k) is an initial estimation of the SOC of the battery pack 101 by using the ampere-hour method, which would be further corrected in the following descriptions.
- a recursive least square (RLS) method is employed to, based on parameters R o(k-1) , R e(k-1) , C e(k-1) , R d(k-1) , C d(k-1) obtained in the previous SOC estimation process and the voltage U (k) and the current I (k) measured at step S 504 , revise and achieve another set of parameters R o(k) , R e(k) , C e(k) , R d(k) , C d(k) at the current time point k in the current SOC estimation process.
- RLS recursive least square
- the parameters R o(k-1) , R e(k-1) , C e(k-1) , R d(k-1) , C d(k-1) obtained in the previous SOC estimation process are stored in the storage unit 1021 . If it is a first SOC estimation process when the battery pack 101 is in use, the R o(k-1) , R e(k-1) , C e(k-1) , R d(k-1) , C d(k-1) would be those parameters identified according to the method described in FIG. 3 .
- G ⁇ ( s ) R 0 + R e 1 + ⁇ e ⁇ s + R d 1 + ⁇ d ⁇ s ( 2 ⁇ - ⁇ 6 )
- the battery transfer function 2-6 could be transformed to a discrete form as following expression 2-7 by way of the bilinear transformation method:
- G ⁇ ( z - 1 ) ⁇ 0 + ⁇ 1 ⁇ Z - 1 + ⁇ 2 ⁇ Z - 2 1 + ⁇ 1 ⁇ Z - 1 + ⁇ 2 ⁇ Z - 2 ( 2 ⁇ - ⁇ 7 )
- Expression 2-10 shows fundamental algorithm of the recursive least square (RLS) method:
- K (k) is a gain factor
- P (k-1) represents a covariance matrix upon the (k ⁇ 1) th measurement.
- the K (k) is set as 5
- P (k-1) is equal to ⁇ I n , wherein ⁇ is a giant number and is set as 10 5 in this disclosure
- I n is an identity matrix of size n, which is the n ⁇ n square matrix, with matrix elements being ones on the main diagonal and zeros elsewhere.
- the achieved parameters R o(k) , R e(k) , C e(k) , R d(k) , C d(k) are stored in the storage unit 1021 , and used as basis of a next SOC estimation process.
- SOC* (k) obtained at step S 505 would be revised by using the parameters R o(k) , R e(k) , C e(k) , R d(k) , C d(k) achieved at step S 506 in light of extended kalman filter (EKF) method.
- EKF extended kalman filter
- X k stands for a variable of state
- a k represents a gain matrix of the variable of state at the time point k
- B k represents a gain matrix of an input variable at the time pint k
- I k is the input variable at the time point k
- W k means process noises.
- Y k is an observer output vector
- C k is a transfer matrix of the variable of state at the time point k
- Z k means system observing value at the time point k.
- a state space equation employed by the EKF method is as following expression:
- SOC* (k) obtained at step S 505 would be served as an input of the expression 2-14.
- U e and U d are voltages respectively loaded upon the capacitors C e and C d , I (k) is the current obtained at the measurement step S 504 .
- the method for revising the SOC* (k) by using the EKF method comprises the following steps:
- Step 1 Determining the Linear Coefficient of the Variable of State
- a k [ 1 0 0 0 ⁇ - T ⁇ e ⁇ ( k ) 0 0 ⁇ - T ⁇ d ⁇ ( k ) ]
- B k [ - k ⁇ ⁇ T / Q 0 R e ( 1 - ⁇ - T ⁇ e ⁇ ( k ) ) R d ( 1 - ⁇ - T ⁇ d ⁇ ( k ) ) ]
- C k [ ⁇ U oc ⁇ SOC ⁇
- SOC S ⁇ O ⁇ ⁇ C ⁇ ( K ) - - 1 - 1 ] T
- D k - R 0 ⁇ ( k )
- Step 2 Initiating the Variable of State
- the initiation of the variable of state comprises initiating three components of X k-1 + with SOC*(k), 0, 0, respectively, and initiating P k-1 + with the varance of X k-1 + , i.e., var(X k-1 + ), as following:
- SOC*(k) is the result obtained at step S 505 .
- U k is the voltage obtained at the measurement step S 504
- P k is mean squared estimation error matrix at time point k
- L k is the Kalman system gain
- D w and D v stand respectively for system noises and measurement noises, which are determined by noises of practical systems.
- the system noise D w and measurement noise D v are both set as normal distribution noises with mean of 0 and squared error of 0.1.
Abstract
Description
- The present invention is related to Li-ion battery management systems, and more particularly to a system and method for state of charge (SOC) estimation of Li-ion batteries.
- As the energy-saving and environmental issues have become increasingly prominent, lithium ion (Li-ion) battery, due to its advantages of high specific energy and green environmental protection, have been widely used as large capacity power supply in various fields such as electronic-powered automotive, aerospace, ship gradually. With the development of li-ion battery technology, energy density of li-ion battery becomes higher and higher, quantity of battery units in a battery pack also becomes larger and larger. After long-time use of a battery pack, asymmetry developed among the batteries in the battery pack can cause one or more of the batteries overcharging or over-discharging, and subsequently lowers the performance of the battery pack in the whole, resulting in serious effect on the service life of the battery pack. Therefore, a battery management system for managing and monitoring the working state of the battery pack is indispensable.
- In practice, state of charge (SOC) is an important reference parameter of the working state of the li-ion battery pack, and is usually employed to indicate remainder energy of the li-ion battery pack. Accurate SOC estimation of the li-ion battery pack utilized in automobiles can not only tell drivers of correct estimated mileage of the automobiles, but also ensure charging/discharging optimization of the li-ion battery pack, which ensures safe utility of the li-ion battery pack. When the automobile is running, large currents may cause the battery pack to be overly discharged and subsequently destroy the battery pack. Therefore, real time collection of voltage, temperature, and charging/discharging current of each battery is important for accurate SOC estimation of the battery so as to prolong the life of the battery pack and increase performance of the automobile.
- The SOC can be estimated based on attribute parameters such as voltage, current, resistance, temperature of the battery. The attribute parameters of the battery generally can change in accordance with the aging of the battery and other uncertain factors, such as random road conditions the automobile is going through.
- Currently, the most popular method for SOC estimation of battery is the so-called ampere-hour method, which is also a relatively accurate method on SOC estimation. The ampere-hour method employs real time current integral to calculate ampere hour, and then revises temperature, self-discharging data and ageing parameters that can affect the SOC estimation, and eventually obtains a relatively accurate SOC value by using a revision function and said parameters. However, the above-mentioned method is still far from being sufficient for practical situations because there are many other factors that could practically affect SOC estimation of the battery, and because it is hard to achieve the revision function in practice. Therefore, to date the SOC value estimated by employing the ampere-hour method can be far from the real SOC value of the battery. Other existing methods for SOC estimation include constant current/voltage method, open circuit voltage method, specific density method, and so on. These methods each have more or less defects that would lead to inaccurate SOC value. Therefore, a novel method needs to be developed for accurate SOC estimation of a battery pack.
- The details as well as other features and advantages of this invention are set forth in the remainder of the specification and are shown in the accompanying drawings.
-
FIG. 1 is a schematic diagram of a battery system device including a battery management system according to embodiments of the present disclosure. -
FIG. 2 is a schematic diagram of the battery management system according to a preferred embodiment. -
FIG. 3 is a flowchart of a method for obtaining parameters of a second order RC equivalent circuit simulating the battery unit in accordance with an exemplary disclosure of the present invention. -
FIG. 4A is an exemplary second order RC equivalent circuit simulating the battery unit. -
FIG. 4B shows the curve line of voltage changes corresponding to the discharging currents with time elapsing. -
FIG. 4C shows a fitting curve representing values of Ued ofFIG. 4A after the discharge current being removed. -
FIG. 4D shows a fitted voltage loaded upon the series circuit composed of the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization ofFIG. 4A . -
FIG. 5 is a flow chart of SOC estimation method of the battery pack in use according to a preferred embodiment of the present disclosure. -
FIG. 6 is an exemplary OCV-SOC mapping table. - In the following description, for purposes of explanation and not limitation, specific details are set forth such as particular architectures, interfaces, techniques, etc. in order to provide a thorough understanding of the present invention. However, it will be apparent to those skilled in the art that the present invention may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.
- In the exemplary embodiment of the present invention, a
battery system device 100, such as an electric bicycle, an electric vehicle or an integrated power storage system, generally comprises abattery pack 101, a battery management system (BMS) 102 for managing thebattery pack 101, and aload 103 powered by thebattery pack 101. In the preferred embodiment of the present invention, thebattery pack 101 may comprise only one single battery, or is composed of many batteries serially connected one by one. In the condition that thebattery pack 101 comprises just one battery, thebattery pack 101 can also be calledbattery 101. For consistency, in this embodiment, a single battery is marked as a battery unit, therefore, thebattery pack 101 may comprise one or more battery units. Thebattery management system 102 is used to manage and maintain thebattery pack 101, comprising but not limited to providing over-voltage and/or over-current protection, state of charge (SOC) estimation of thebattery pack 101. In the exemplary embodiment, thebattery pack 101 and thebattery management system 102 collectively form apower system 10 of thebattery system device 100. Theload 103 may be any kind of power consumption device, such as motors employed by the electric bicycle or the electric vehicle. In the exemplary embodiment, thebattery pack 101 or the battery unit described hereinafter is lithium ion (Li-ion) typed. -
FIG. 2 is a detailed description of thebattery management system 102. In the preferred embodiment of the present disclosure, thebattery management system 102 comprises acontroller 1020, astorage unit 1021, ameasurement module 1022, abalance module 1023, apower supply 1024, acommunication module 1025, and aprotection module 1026. In a preferred embodiment, themeasurement module 1022, thecommunication module 1025 and theprotection module 1026 are electronically connected to thecontroller 1020 by way of photoelectric coupling isolation circuits (PCIC), respectively. - In the preferred embodiment, the
storage unit 1021 may be a memory integrated with thecontroller 1020, such as a flash memory, a static random access memory (SRAM), an electrically erasable programmable read-only memory (EEPROM). In other embodiments, thestorage unit 1021 may be those storage apparatuses independent from but electrically connected to thecontroller 1021, such as a solid state disk or a micro hard disk. In alternative embodiments, thestorage unit 1021 may be the combination of the memory and the storage apparatus. Thestorage unit 1021 stores program codes that can be executed by thecontroller 1021 to maintain thebattery pack 101, for example, estimating the SOC of thebattery pack 101. Thestorage unit 1021 is also used to store data generated during the SOC estimation in accordance with the preferred embodiment of the present invention. - In an exemplary embodiment, the
controller 1020 is the PIC18F458 type chip produced by Microchip Technology Incorporation. The PIC18F458 type chip is an 8-bit micro-controller with 32 Kilobytes memory space for storing program codes, 1536 bytes-sized SRAM and 256 bytes-sized EEPROM. The PIC18F458 type chip further comprises a clock with frequency up to 32 MegaHz, which can operate as an external crystal oscillator circuit for, such as timing. Another outstanding feature of the PIC18F458 type chip is that the power consumption of the chip is relatively low. For example, under a sleeping mode, current consumed by the chip is only 0.2 μA. - The
measurement module 1022 is an important part of thebattery management system 102, and is used for measuring parameters such as current, voltage, and temperature of thebattery pack 101. The measured parameters are not only basis for the SOC estimation, but security guarantee of thebattery pack 101. In this disclosure, the current measurement comprises measuring the current flowing through thebattery pack 101. The voltage measurement comprises measuring the voltage of thebattery pack 101, and the voltage of each battery unit in thebattery pack 101, providing thebattery pack 101 comprises more than one battery units. The temperature measurement comprises measuring the temperature of the battery units in thebattery pack 101. In this embodiment, the measured current, voltage and temperature data are stored in thestorage unit 1021. - The
balance module 1023 is used for balancing the charging/discharging current of thebattery pack 101. In an exemplary embodiment, thebalancing module 1023 comprises monitoring chip LTC6803 and other accompanied circuits provided by Linear Technology Corporation. - The
power supply 1024 is used for providing various reference voltages to different function modules, such as thecontroller 1021, the monitoring chip LTC6803, etc. In this embodiment, thepower supply 1024 employs the LM2009 chip to provide fixed voltage, such as 15V, and uses other DC/DC regulator, such as LM7805 chip, to convert the fixed 15V voltage to various DC outputs, such as 10V, 5V, and provides the converted DC outputs to the different function modules. - The
communication module 1025 provides communication services for above-mentioned function modules. The communication services comprise a serial peripheral interface (SPI), a controller area network (CAN), and serial communication. - The
protection module 1026 provides protections under the conditions such as over-current, over-voltage, under-voltage, to ensure normal operation of thebattery system device 100, especially thebattery pack 101. - In the preferred embodiment of the present invention, before the SOC estimation of the
battery pack 101, electrochemical parameters of thebattery pack 101 should be obtained. - 1. Obtaining Electrochemical Parameters of the Battery Pack
- The electrochemical performance of the battery unit in the
battery pack 101 generally presents non-linear characteristic. In the preferred disclosure of the present invention, a second order RC equivalent circuit is employed to simulate the battery unit, in that the second order RC equivalent circuit owns the benefits comprising, for example, the order is relatively low, and parameters thereof are easily to be obtained. Particularly, the second order RC equivalent circuit can be easily constructed and preferably simulates dynamic characteristics of the battery unit under working status. The dynamic characteristics comprise an ohmic resistance R0, an electrochemical polarization resistance Re, an electrochemical polarization capacitor Ce, a concentration polarization resistance Rd and a concentration polarization capacitor Cd. By using certain methods acted on the battery unit, the dynamic characteristics can be achieved to describe the electrochemical characteristics of the battery unit, subsequently to aid to estimate the state of charge of the battery unit in thebattery pack 101. -
FIG. 3 is a flowchart of a method for obtaining parameters of the second order RC equivalent circuit simulating the battery unit in accordance with an exemplary disclosure of the present invention. At step S301, the second order RC equivalent circuit simulating the battery unit is constructed, as illustrated inFIG. 4A . The constructed second order RC equivalent circuit comprises the ohmic resistance R0, the electrochemical polarization resistance Rc, the electrochemical polarization capacitor Cc, the concentration polarization resistance Rd and the concentration polarization capacitor Cd. In the exemplary embodiment, the electrochemical polarization resistance Re, the electrochemical polarization capacitor Ce are connected in parallel to form a capacitance-resistance segment of electrochemical polarization, and the concentration polarization resistance Rd and the concentration polarization capacitor Cd are connected in parallel to form a capacitance-resistance segment of concentration polarization. The ohmic resistance R0, the capacitance-resistance segment of electrochemical polarization, and the capacitance-resistance segment of concentration polarization are connected in series. In this embodiment, values of the ohmic resistance R0, the electrochemical polarization resistance Re, the electrochemical polarization capacitor Ce, the concentration polarization resistance Rd and the concentration polarization capacitor Cd would be obtained by an experimental method. In the exemplary embodiment, the experimental method comprises the following steps: - Step I: keeping the battery unit under a status without charging or discharging over one hour;
- Step II: discharging the battery unit continuously lasting for 900 seconds;
- Step III: keeping the battery unit under the status without charging or discharging for a relatively long time, such as one hour, up to the end of the experiment.
-
FIG. 4B shows the curve line of voltage changes corresponding to the discharging currents with time elapsing. - 1.1 Identifying R0
- In the exemplary embodiment, firstly, placing one battery unit under a status without charging or discharging over one hour; then loading a discharge current pulse on a positive electrode and a negative electrode of the battery unit, for example, at the time point of 1800 s shown in
FIG. 4B . When the discharge current pulse is loaded, a responsive voltage upon between the positive electrode and the negative electrode of the battery unit gets a downward mutation. After continuous discharging for 900 s, at the time point of the 2700 s shown inFIG. 4B , the discharge current pulse is removed, and the responsive voltage upon the positive electrode and the negative electrode of the battery unit gets a upward mutation. In the exemplary disclosure, the status without charging or discharging of the battery could be defined as an idle status of the battery unit. - Based on the second order RC equivalent circuit simulating the battery unit, the voltage upon the capacitance-resistance segment of concentration polarization composed of the concentration polarization resistance Rd and the concentration polarization capacitor Cd could not get a mutation, due to the exist of the concentration polarization capacitor Cd. In the same way, due to the exist of the electrochemical polarization capacitor Cc, the voltage upon the capacitance-resistance segment of electrochemical polarization composed of the electrochemical polarization resistance Re and the electrochemical polarization capacitor Ce also could not get a mutation. Therefore, the downward mutation and the upward mutation can only be caused by the ohmic resistance R0. According to the above-mentioned theory, at step S303, the resistance value R0 of the ohmic resistance R0 can be calculated by dividing the voltage change by the discharge current, that is, R0=ΔU/It.
- 1.2 Identifying τe and τd
- In the exemplary disclosure, time constants of the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization are respectively represented by τe and τd, wherein τe=Re*Ce, τd=Rd*Cd.
- As shown in
FIG. 4B , after the time point of 2700 s, the discharge current pulse is removed, the voltage upon the ohmic resistance R0 would become 0 due to no currents flowing through the ohmic resistance R0. The upward mutation of the voltage upon the positive electrode and the negative electrode of the battery unit, as shown inFIG. 4B , is caused by electronic discharging of the capacitors Cd and Ce, which are charged when the battery unit is discharging. When the capacitors Cd and Ce discharge completely, the voltage upon the positive electrode and the negative electrode of the battery unit becomes stable at, for example, about 4.05 v as shown inFIG. 4B . - After the discharge current pulse being removed, the voltage upon the positive electrode and the negative electrode of the battery unit is close to electromotive force of the battery unit. In an exemplary disclosure, supposing the time point t of removing the discharge current pulse is zero time, the voltage response upon the positive electrode and the negative electrode of the battery unit at the time point t could be deemed as zero state response:
-
- Here, Ud and Ue respectively stand for voltages loaded upon the capacitors Cd and Ce, Ued represents the voltage loaded upon the series circuit composed of the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization.
FIG. 4C shows a fitting curve representing values of Ued after the discharge current being removed. At step S305, with the MATLAB software provided by the MathWorks company, τe and τd can be achieved by using the exponent curve fitting module in the MATLAB software. - 1.3 Identifying Rd and Re
- According to the above-mentioned description, at the time point of 1800 s shown in
FIG. 4B , the voltage upon the positive electrode and the negative electrode of the battery unit declines due to the ohmic resistance R0, the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization. Because the battery unit has been place at the idle status for a long time period, such as over one hour, the capacitor Ce in the capacitance-resistance segment of electrochemical polarization and the capacitor Cd in the capacitance-resistance segment of concentration polarization could be deemed as zero electron state. In an exemplary disclosure, supposing the time point t of loading the discharge current pulse is zero time, voltage response of the capacitor Ce and Cd to the discharge current pulse is a zero state response: -
- The voltage Ued loaded upon the series circuit composed of the capacitance-resistance segment of electrochemical polarization and the capacitance-resistance segment of concentration polarization is fitted as shown in
FIG. 4D . At step S307, putting the calculated τe and τd into the formula 1-2, and by using the exponent curve fitting method of the MATLAB software, the Rd and Re can be achieved. - 1.4 Identifying Cd and Ce
- At step S309, based on the formula τ=R*C, it can be concluded that C is equal to τ/R. Because the time constants τe and τd, the resistance values Rd and Re of the polarization resistances Rd and Re are achieved, it is easy to achieve values Cd and Ce of the polarization capacitors Cd and Ce.
- In the preferred embodiment of the disclosure, the second order RC equivalent circuit is used to simulate the battery unit. Therefore, the identified parameters of the second order RC equivalent circuit are indeed dynamic characteristic paramters of the battery unit. According to above-described identifying processes, in one embodiment of the present disclosure, identified dynamic characteristic parameters of one exemplary battery unit are shown in following table:
-
Parameter Identified value R0 34.5 mΩ Rd 21.45 mΩ Re 10.6 mΩ Cd 4768.74 F Ce 3678.21 F - In the preferred embodiment, the parameters of the second order RC equivalent circuit simulating the battery unit which are identified by way of said-mentioned method are stored in the
storage unit 1021, and would be invoked in estimating processed when the battery unit is used in practice. - 2. SOC Estimation of the Battery Pack in Use
-
FIG. 5 is a flow chart of SOC estimation method of thebattery pack 101 in use according to a preferred embodiment of the present disclosure, which could be applied to thebattery system device 100 ofFIG. 1 . For example, the SOC estimation method could be realized by thecontroller 1020 executing the programmed codes stored in thestorage unit 1021. In the preferred embodiment, the SOC estimation method is periodically executed by thecontroller 1020 to achieve a relatively accurate SOC of thebattery pack 101 in different time points. A time cycle the method being executed periodically could be, for example, 10 millisecond (ms). In other embodiment, the time period could be set or configured as other value, such 15 ms. In the preferred embodiment, an once-through operation of the SOC estimation method is described as one SOC estimation process. - In the exemplary embodiment, the clock in the
controller 1020 would record the time duration of thebattery pack 101 being in the idle state in which there is no charging or discharging occurred to thebattery pack 101. For simplicity of description, the time duration of thebattery pack 101 in the idle state is simplified as an idle time of thebattery pack 101. At the beginning of each SOC estimation process, some program codes would be executed to determine whether the idle time of thebattery pack 101 is longer than a predefined time period, such as one hour, at step S501. - If the idle time of the
battery pack 101 is not longer than the predefined time period, at step S502A, thecontroller 1020 reads a previous SOC value (here marked as SOC(k-1)) generated during a previous SOC estimation process, and regards the SOC(k-1) as an initial SOC value of a current SOC estimation process. - In the preferred embodiment of the present invention, if the
battery pack 101 is in the idle state for a relatively long time, for example, over one hour, an open circuit voltage (OCV) of thebattery pack 101 is substantially equal to an electromotive force of thebattery pack 101. Therefore, if the idle time of thebattery pack 101 is longer than the predefined time period, at step S502B, themeasurement module 1022 measures a current voltage upon the positive electrode and the negative electrode of thebattery pack 101, subsequently, thecontroller 1020 queries an OCV-SOC mapping table to determine a SOC value corresponding to the current voltage, and records the determined SOC value as SOC(k-1). The SOC(k-1) would act as an initial SOC value of a current SOC estimation process. - It should be noted, step S502A and step S502B are alternative in one SOC estimation process according to a preferred embodiment of the present invention.
- In the preferred embodiment, the OCV-SOC mapping table is established by an experimental method. The experiment method for establishing the OCV-SOC mapping table comprises the steps of:
- A) charging the
battery pack 101, fully and completely; - B) placing the
battery pack 101 in the idle status, without charging or discharging, for over one hour; - C) discharging the
battery pack 101 to reduce 5 percent of the SOC of thebattery pack 101 by using a programmable electronic load, recording the open circuit voltage of thebattery pack 101 and subsequently placing thebattery pack 101 in the idle state for over one hour; - D) repeating said step C), under constant temperature of substantially 20 degree Celsius, until the
battery pack 101 discharges completely; - E) establishing the OCV-SOC mapping table based on the recorded open circuit voltage and corresponding SOC of the
battery pack 101. - In this exemplary disclosure, the programmable electronic load employs Chroma 6310A type programmable DC electronic load provided by Chroma Company.
FIG. 6 is an exemplary OCV-SOC mapping table of an 18650 type li-ion battery, which is obtained by using the experiment method for establishing the OCV-SOC mapping table. - At step S503, the
controller 1020 determines whether thebattery pack 101 is in a working status. In the exemplary disclosure, the working status is opposite to the idle status, that is, the working status means thebattery pack 101 is discharging and/or charging. In other words, thecontroller 1020 determines thebatter pack 101 is in the working status if thebattery pack 101 is discharging and/or charging. If thebattery pack 101 is not in the working status, or is in the idle status, the flow returns to step S501 to determine whether an idle time of thebattery pack 101 is longer than the predefined time period in another SOC estimation process. - If the
battery pack 101 is in the working status, at step S504, themeasurement module 1022 measures current I(k) flowing through and voltage U(k) upon the two polarities of thebattery pack 101. In the preferred embodiment of the present disclosure, when measuring the current I(k), a direction of the measured current I(k) shows the discharging or charging status of thebattery pack 101. Hereinafter, if thebattery pack 101 is discharging, the current I(k) is positive, and if thebattery pack 101 is charging, the current I(k) is negative. - At step S505, an ampere hour method is employed to estimate SOC(k) in current SOC estimation process, based on the determined SOC(k-1) in step S502A or S502B.
- 2.1 Ampere Hour Method
- The ampere hour method is also called current integration method, which is a fundamental method to measure the SOC of batteries. In current SOC estimation process, after determining the SOC(k-1) in step S502A or S502B, the
battery pack 101 is charging or charging from time point k-1 to k, such as 10 ms, a change of the SOC of thebattery pack 101 could be represented as: -
- Here, Q0 is a total quantity of electric discharge in fixed current of 0.1*C of the
battery pack 101, wherein C represents the nominal capacity of thebattery pack 101. I(k) represents the charging or discharging current of thebattery pack 101, wherein I(k) is positive if thebattery pack 101 is discharged, otherwise the I(k) is negative. - Peukert proposes an empirical formula to revise the SOC of the
battery pack 101 working with changing currents, as shown in following formula 2-2: -
In *t=K (2-2) - Here, I represents discharging current, t is a discharging time length, n and K are constants that are determined by types and active materials of the
battery pack 101. In the exemplary disclosure, the active material is lithium. - By multiplying two sides of the formula 2-2 by I1-n, a new expression is obtained:
-
Q=I*t=I l-n *K (2-3) - Here, if I=I0=0.1*C, Q is equal to Q0 mentioned above; if I=0.5*C, Q is marked as Q1. According to expression 2-3, Q0=I0 1-n*K, Q1=I1-n*K. Therefore, Q1/Q0=(I/I0)1-n. Assuming Q1/Q0=η, Q1=η*Q0. If kη=1/η, an initial SOC estimation value SOC*(k) considering charge/discharge rate is:
-
SOC*(k)=SOC(k-1) −k η *I(k)*T/Q 0 (2-4) - Here, kη is a revise value to the capacity volume of the
battery pack 101. SOC*(k) is an initial estimation of the SOC of thebattery pack 101 by using the ampere-hour method, which would be further corrected in the following descriptions. - 2.2 Recursive Least Square (RLS) Method
- At step S506, a recursive least square (RLS) method is employed to, based on parameters Ro(k-1), Re(k-1), Ce(k-1), Rd(k-1), Cd(k-1) obtained in the previous SOC estimation process and the voltage U(k) and the current I(k) measured at step S504, revise and achieve another set of parameters Ro(k), Re(k), Ce(k), Rd(k), Cd(k) at the current time point k in the current SOC estimation process. The parameters Ro(k-1), Re(k-1), Ce(k-1), Rd(k-1), Cd(k-1) obtained in the previous SOC estimation process are stored in the
storage unit 1021. If it is a first SOC estimation process when thebattery pack 101 is in use, the Ro(k-1), Re(k-1), Ce(k-1), Rd(k-1), Cd(k-1) would be those parameters identified according to the method described inFIG. 3 . - Based on the battery model shown in
FIG. 4A , and according to Kirchhoff laws and Laplace transform, it could be drawn upon changing time domain t to Laplace domain s: -
- Here, the most right part of the formula is defined as battery transfer function:
-
- The battery transfer function 2-6 could be transformed to a discrete form as following expression 2-7 by way of the bilinear transformation method:
-
- An difference equation of the expression 2-7 could be expressed as following:
-
U(k)=−α1 U(k−1)−α2 U(k−2)+β0 I(k)+β1 I(k−1)+β2 I(k−2) (2-8) - Here, given θ=[α1 α2 β0 β1 β2], hT(k)=[−U(k−1)−U(k−2) I(k) I(k−1) I(k−2)], it can be concluded:
-
U(k)=h T(k)θ+e(k) (2-9) - Expression 2-10 shows fundamental algorithm of the recursive least square (RLS) method:
-
- Here, K(k) is a gain factor, P(k-1) represents a covariance matrix upon the (k−1)th measurement. In the preferred embodiment of the present disclosure, the K(k) is set as 5, and P(k-1) is equal to αIn, wherein α is a giant number and is set as 105 in this disclosure, In is an identity matrix of size n, which is the n×n square matrix, with matrix elements being ones on the main diagonal and zeros elsewhere.
- According to the formula (2-10), {circumflex over (θ)}(k) could be calculated, which could be regarded as a current θ. Because θ=[α1 α2 β0 β1 β2], α1, α2, β0, β1, β2 could be subsequently obtained.
- Putting α1, α2, β0, β1, β2 into an inverse equation 2-11, the parameters Ro(k), Re(k), Ce(k), Rd(k), Cd(k) at the current time point k in the current SOC estimation process could be achieved. In the inverse equation 2-11, T is a sampling cycle, which is the time interval of the RLS method being executed.
-
- The achieved parameters Ro(k), Re(k), Ce(k), Rd(k), Cd(k) are stored in the
storage unit 1021, and used as basis of a next SOC estimation process. - At step S507, SOC*(k) obtained at step S505 would be revised by using the parameters Ro(k), Re(k), Ce(k), Rd(k), Cd(k) achieved at step S506 in light of extended kalman filter (EKF) method.
- 2.3 Extended Kalman Filter (EKF) Method
- Equation of state employed by the EKF method in the preferred embodiment could be expressed as following:
-
X k =A k X k-1 +B k I k +W k (2-12) - Here, Xk stands for a variable of state, Ak represents a gain matrix of the variable of state at the time point k, Bk represents a gain matrix of an input variable at the time pint k, Ik is the input variable at the time point k, Wk means process noises.
- An observer output equation of the EKF method is:
-
Y k =C k X k +Z k (2-13) - Here, Yk is an observer output vector, Ck is a transfer matrix of the variable of state at the time point k, Zk means system observing value at the time point k.
- A state space equation employed by the EKF method is as following expression:
-
- Here, SOC*(k) obtained at step S505 would be served as an input of the expression 2-14.
- In the expression 2-14, Ue and Ud are voltages respectively loaded upon the capacitors Ce and Cd, I(k) is the current obtained at the measurement step S504.
- In the exemplary embodiment of the present disclosure, the method for revising the SOC*(k) by using the EKF method comprises the following steps:
- Step 1: Determining the Linear Coefficient of the Variable of State
- Based on the parameters Ro(k), Re(k), Ce(k), Rd(k), Cd(k) achieved at step S506, one can compute τe and τd. Then, the following coefficient matrixes Ak, Bk and Dk can be calculated. According to the OCV-SOC curve of the
battery pack 101, coefficient matrixes Ck can be obtained. -
- Step 2: Initiating the Variable of State
- In this disclosure, the initiation of the variable of state comprises initiating three components of Xk-1 + with SOC*(k), 0, 0, respectively, and initiating Pk-1 + with the varance of Xk-1 +, i.e., var(Xk-1 +), as following:
- Xk-1 +=[SOC*(k) 0 0]T
- Pk-1 +=var(Xk-1 +)
- Here, SOC*(k) is the result obtained at step S505.
- Step 3: EKF Iteration
- The following expression 2-15 is an EKF iteration formula:
-
- In the expression 2-15, Uk is the voltage obtained at the measurement step S504, Pk is mean squared estimation error matrix at time point k, Lk is the Kalman system gain, Dw and Dv stand respectively for system noises and measurement noises, which are determined by noises of practical systems. In the exemplary disclosure, the system noise Dw and measurement noise Dv are both set as normal distribution noises with mean of 0 and squared error of 0.1.
- Iterating with the expression 2-15, a best variable of state would be achieved: Xk +=[SOC(k) Ue(k) Ud(k)]T. According to the best variable of state, a best state of charge, SOC(, can be retrieved.
- While the foregoing description and drawings represent the preferred embodiments of the present invention, it will be understood that various additions, modifications and substitutions may be made therein without departing from the spirit and scope of the present invention as defined in the accompanying claims. In particular, it will be clear to those skilled in the art that the present invention may be embodied in other specific forms, structures, arrangements, proportions, and with other elements, materials, and components, without departing from the spirit or essential characteristics thereof The presently disclosed embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims, and not limited to the foregoing description.
Claims (18)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/586,062 US10175302B2 (en) | 2013-06-20 | 2017-05-03 | Power system and state of charge estimation |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310245252 | 2013-06-20 | ||
CN2013102452526 | 2013-06-20 | ||
CN201310245252.6A CN104237791A (en) | 2013-06-20 | 2013-06-20 | Lithium battery charge state estimation method, battery management system and battery system |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/586,062 Division US10175302B2 (en) | 2013-06-20 | 2017-05-03 | Power system and state of charge estimation |
Publications (2)
Publication Number | Publication Date |
---|---|
US20140316728A1 true US20140316728A1 (en) | 2014-10-23 |
US9709635B2 US9709635B2 (en) | 2017-07-18 |
Family
ID=51729660
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/192,867 Expired - Fee Related US9709635B2 (en) | 2013-06-20 | 2014-02-27 | System and method for SOC estimation of a battery |
US15/586,062 Active US10175302B2 (en) | 2013-06-20 | 2017-05-03 | Power system and state of charge estimation |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/586,062 Active US10175302B2 (en) | 2013-06-20 | 2017-05-03 | Power system and state of charge estimation |
Country Status (2)
Country | Link |
---|---|
US (2) | US9709635B2 (en) |
CN (1) | CN104237791A (en) |
Cited By (48)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150046107A1 (en) * | 2012-10-26 | 2015-02-12 | Lg Chem, Ltd. | Apparatus and method for estimating state of charge of battery |
CN104502858A (en) * | 2014-12-31 | 2015-04-08 | 桂林电子科技大学 | Power battery SOC estimation method based on backward difference discrete model and system thereof |
US20150251556A1 (en) * | 2014-03-05 | 2015-09-10 | Ford Global Technologies, Llc | Battery model with robustness to cloud-specific communication issues |
CN104965179A (en) * | 2015-07-06 | 2015-10-07 | 首都师范大学 | Lithium ion storage battery temperature combinational circuit model and parameter identification method thereof |
CN105044440A (en) * | 2015-08-26 | 2015-11-11 | 苏州弗尔赛能源科技股份有限公司 | Fuel cell monolithic voltage inspection system based on LTC6803 |
CN105045971A (en) * | 2015-07-01 | 2015-11-11 | 西安交通大学 | Lithium battery fractional order discretization impedance model |
CN105093128A (en) * | 2015-08-31 | 2015-11-25 | 山东智洋电气股份有限公司 | Storage battery state of charge (SOC) estimation method based on extended Kalman filtering (EKF) |
CN105699907A (en) * | 2016-01-28 | 2016-06-22 | 广州市香港科大霍英东研究院 | A battery SOC estimation method and system based on dynamic impedance correction |
CN105738817A (en) * | 2016-01-29 | 2016-07-06 | 深圳市沃特玛电池有限公司 | Battery charge state estimation method based on AEKF and estimation system |
US9533597B2 (en) | 2014-03-05 | 2017-01-03 | Ford Global Technologies, Llc | Parameter identification offloading using cloud computing resources |
CN106970328A (en) * | 2017-01-17 | 2017-07-21 | 深圳市沛城电子科技有限公司 | A kind of SOC estimation method and device |
CN107390127A (en) * | 2017-07-11 | 2017-11-24 | 欣旺达电动汽车电池有限公司 | A kind of SOC estimation method |
CN107643488A (en) * | 2016-07-21 | 2018-01-30 | 神讯电脑(昆山)有限公司 | Corresponding to the metering method and its electronic installation of the battery electric quantity of temperature |
US20180147952A1 (en) * | 2016-11-25 | 2018-05-31 | Hyundai Motor Company | Method and system for controlling motors |
WO2018127296A1 (en) * | 2017-01-09 | 2018-07-12 | Volvo Truck Corporation | A method and arrangement for determining the state of charge of a battery pack |
CN108287316A (en) * | 2018-01-15 | 2018-07-17 | 厦门大学 | Accumulator method for estimating remaining capacity based on threshold spread Kalman Algorithm |
JP2018151175A (en) * | 2017-03-10 | 2018-09-27 | 株式会社デンソーテン | Estimation device, estimation method, and estimation program |
CN109031138A (en) * | 2018-06-29 | 2018-12-18 | 上海科列新能源技术有限公司 | A kind of safety evaluation method and device of power battery |
CN109061505A (en) * | 2018-08-28 | 2018-12-21 | 淮阴工学院 | A kind of detection method of lithium battery SOH |
CN109188293A (en) * | 2018-11-08 | 2019-01-11 | 武汉理工大学 | Based on new breath EKF lithium ion battery SOC estimation method of the covariance with fading factor |
US20190041467A1 (en) * | 2015-12-01 | 2019-02-07 | Omron Corporation | Remaining battery charge estimation system and remaining battery charge estimation method |
CN109726365A (en) * | 2018-12-05 | 2019-05-07 | 新奥数能科技有限公司 | A kind of method and apparatus of load forecast |
CN109901082A (en) * | 2017-12-08 | 2019-06-18 | 南京德朔实业有限公司 | Portable electric energy system and its measurement method |
JP2019124612A (en) * | 2018-01-18 | 2019-07-25 | 日立オートモティブシステムズ株式会社 | Secondary battery system |
CN110221219A (en) * | 2019-07-03 | 2019-09-10 | 中国民用航空飞行学院 | Airborne circumstance is got off the plane lithium battery SOC estimation method |
CN110245366A (en) * | 2018-03-08 | 2019-09-17 | 华为技术有限公司 | Dynamic power consumption estimation method, apparatus and system |
CN110275113A (en) * | 2019-06-25 | 2019-09-24 | 内蒙古工业大学 | A kind of lithium battery charge state estimation method |
CN110361652A (en) * | 2019-06-26 | 2019-10-22 | 河南理工大学 | A kind of Kalman filtering lithium battery SOC estimation method based on Model Parameter Optimization |
CN110414117A (en) * | 2019-07-23 | 2019-11-05 | 北京航空航天大学 | A kind of soft bag lithium ionic cell sealed reliable degree prediction technique |
CN110441702A (en) * | 2019-07-31 | 2019-11-12 | 湘潭大学 | A method of lithium ion battery charge capacity is estimated with Extended Kalman filter |
US10518640B2 (en) | 2015-11-17 | 2019-12-31 | Omron Corporation | Battery remaining capacity display device, battery system, and battery remaining capacity display method |
CN110837622A (en) * | 2019-11-26 | 2020-02-25 | 国网湖南省电力有限公司 | Lithium battery state of charge estimation method based on high-rate discharge |
CN111002834A (en) * | 2019-12-18 | 2020-04-14 | 中车株洲电力机车有限公司 | Method and system for predicting cruising mileage of residual electric quantity of traction battery of railway vehicle |
CN112119317A (en) * | 2019-01-23 | 2020-12-22 | 株式会社Lg化学 | Battery management device, battery management method and battery pack |
CN112255545A (en) * | 2019-07-05 | 2021-01-22 | 西南科技大学 | Lithium battery SOC estimation model based on square root extended Kalman filter |
CN112347415A (en) * | 2020-10-02 | 2021-02-09 | 广东电网有限责任公司广州供电局 | Prediction method based on short-circuit current zero crossing point prediction system |
EP3730957A4 (en) * | 2018-02-20 | 2021-03-17 | Lg Chem, Ltd. | Charge capacity calculation device and method for energy storage system |
CN112858929A (en) * | 2021-03-16 | 2021-05-28 | 上海理工大学 | Battery SOC estimation method based on fuzzy logic and extended Kalman filtering |
CN112881928A (en) * | 2021-03-24 | 2021-06-01 | 东风汽车集团股份有限公司 | Screening method for consistency of battery monomers |
CN113030555A (en) * | 2020-03-18 | 2021-06-25 | 深圳大学 | Energy storage open-circuit voltage estimation method and device, terminal equipment and storage medium |
US20210199724A1 (en) * | 2019-02-22 | 2021-07-01 | Lg Chem, Ltd. | Battery management system, battery management method, battery pack and electric vehicle |
US11069926B1 (en) * | 2019-02-14 | 2021-07-20 | Vcritonc Alpha, Inc. | Controlling ongoing battery system usage via parametric linear approximation |
CN113343443A (en) * | 2021-05-24 | 2021-09-03 | 国网综合能源服务集团有限公司 | Power distribution method for lithium battery prefabricated cabins with different SOC (state of charge) |
CN113341330A (en) * | 2021-05-25 | 2021-09-03 | 西南大学 | Lithium-sulfur power battery SOC estimation method based on OCV correction and Kalman filtering algorithm |
CN113900031A (en) * | 2021-09-14 | 2022-01-07 | 国网浙江省电力有限公司电力科学研究院 | SOC safety verification method after energy storage system is accessed to AGC |
CN114397581A (en) * | 2021-12-09 | 2022-04-26 | 国网天津市电力公司 | New energy automobile battery SOC anti-interference assessment method oriented to charging monitoring data of direct-current charging pile |
US11502530B2 (en) * | 2017-12-26 | 2022-11-15 | Panasonic Intellectual Property Management Co., Ltd. | Battery management device, battery system, and vehicle power supply system for managing battery state of charge level when in non-use state |
CN117811169A (en) * | 2024-02-29 | 2024-04-02 | 威海凯瑞电气股份有限公司 | Battery active equalization management system |
Families Citing this family (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015535993A (en) | 2012-09-25 | 2015-12-17 | スクート ネットワークス, インコーポレイテッドScoot Networks, Inc. | Vehicle access control system and method |
CN104535935B (en) * | 2014-12-31 | 2017-07-21 | 普天新能源车辆技术有限公司 | A kind of capacity check method and device of power battery pack |
JP2018509880A (en) * | 2015-01-13 | 2018-04-05 | ボルボ カー コーポレイション | Method and apparatus for determining the value of the energy state of a battery in an automobile |
CN104914380B (en) * | 2015-06-05 | 2018-11-06 | 上海科梁信息工程股份有限公司 | Recognize the method and system of SOC |
CN104868558B (en) * | 2015-06-10 | 2017-09-29 | 合肥联宝信息技术有限公司 | A kind of charge reminder apparatus and method of external battery for electronic equipment |
CN105044610B (en) * | 2015-07-10 | 2019-01-22 | 西安交通大学 | The high accuracy battery electricity evaluation method of current detecting is not necessarily to based on port voltage |
KR102468895B1 (en) * | 2015-07-21 | 2022-11-21 | 삼성전자주식회사 | Method and apparatus for estimating state of battery |
CN105445665A (en) * | 2015-11-12 | 2016-03-30 | 华晨汽车集团控股有限公司 | Method for estimating state of charge of battery through Kalman filtering |
CN105548896B (en) * | 2015-12-25 | 2019-04-09 | 南京航空航天大学 | Power battery SOC line closed loop estimation method based on N-2RC model |
CN105891721A (en) * | 2016-04-01 | 2016-08-24 | 深圳市清友能源技术有限公司 | SOC test method and SOC test device for battery management system |
CN105699910A (en) * | 2016-04-21 | 2016-06-22 | 中国计量大学 | Method for on-line estimating residual electric quantity of lithium battery |
CN106291375A (en) * | 2016-07-28 | 2017-01-04 | 河南许继仪表有限公司 | A kind of SOC estimation method based on cell degradation and device |
CN109564471B (en) * | 2016-08-12 | 2022-08-23 | 波士顿科学国际有限公司 | Distributed interactive medical visualization system with primary/secondary interaction features |
US10585552B2 (en) | 2016-08-12 | 2020-03-10 | Boston Scientific Scimed, Inc. | Distributed interactive medical visualization system with user interface features |
CN106443459A (en) * | 2016-09-06 | 2017-02-22 | 中国第汽车股份有限公司 | Evaluation method of state of charge of vehicle lithium ion power battery |
CN106707181A (en) * | 2016-12-05 | 2017-05-24 | 澳特卡新能源科技(上海)有限公司 | Cell parameter and charged state estimation method of lithium ion |
CN106816965B (en) * | 2017-01-18 | 2019-08-02 | 同济大学 | A kind of discrimination method of resonance type wireless charging system coil self-induction |
CN106951605A (en) * | 2017-03-02 | 2017-07-14 | 西南科技大学 | A kind of Li-ion batteries piles equivalent model construction method |
CN107390138B (en) * | 2017-09-13 | 2019-08-27 | 山东大学 | Power battery equivalent circuit model parameter iteration new method for identifying |
CN107643437B (en) * | 2017-11-10 | 2020-09-11 | 国家电网公司 | On-line intelligent conversion complex ratio current transformer |
CN107765187A (en) * | 2017-11-14 | 2018-03-06 | 佛山科学技术学院 | A kind of lithium battery charge state evaluation method |
CN109917298A (en) * | 2017-12-13 | 2019-06-21 | 北京创昱科技有限公司 | A kind of cell charge state prediction method and system |
US10630084B2 (en) * | 2017-12-21 | 2020-04-21 | International Business Machines Corporation | Battery management system for extending service life of a battery |
CN108400393B (en) * | 2018-01-17 | 2020-10-23 | 广州市香港科大霍英东研究院 | Battery management method and system suitable for echelon battery |
US11468503B2 (en) * | 2018-04-16 | 2022-10-11 | Bird Rides, Inc. | On-demand rental of electric vehicles |
WO2019204145A1 (en) | 2018-04-20 | 2019-10-24 | Bird Rides, Inc. | Remotely controlling use of an on-demand electric vehicle |
CN108802624B (en) * | 2018-06-19 | 2021-08-31 | 杭州电子科技大学 | Lithium battery SOC estimation method |
US11263690B2 (en) | 2018-08-20 | 2022-03-01 | Bird Rides, Inc. | On-demand rental of electric vehicles |
KR102259415B1 (en) * | 2018-08-29 | 2021-06-01 | 주식회사 엘지에너지솔루션 | Battery management apparatus, battery management method, battery pack and electric vehicle |
CN109581223B (en) * | 2018-11-29 | 2020-08-11 | 吉林大学 | Kalman filtering based core temperature estimation method of lithium ion battery pack |
CN109782188B (en) * | 2019-02-11 | 2022-03-04 | 北京东方计量测试研究所 | SOC characteristic calibration method of spacecraft battery pack simulator |
CN109669134A (en) * | 2019-02-27 | 2019-04-23 | 浙江科技学院 | A kind of evaluation method of the SOC based on Kalman filtering method |
CN110308396B (en) * | 2019-07-03 | 2020-09-25 | 华人运通(江苏)技术有限公司 | Battery state monitoring method, edge processor, system and storage medium |
CN110554320A (en) * | 2019-09-24 | 2019-12-10 | 东风航盛(武汉)汽车控制系统有限公司 | SOC estimation method of lithium ion battery |
CN110554321B (en) * | 2019-09-26 | 2021-05-11 | 长沙理工大学 | Method for detecting SOC (state of charge) of retired power battery in real time |
CN112213653B (en) * | 2019-10-30 | 2023-05-16 | 蜂巢能源科技有限公司 | Battery core state of charge estimation method of power battery and battery management system |
CN112904209B (en) * | 2019-11-19 | 2023-04-18 | 北京车和家信息技术有限公司 | Power battery power duration prediction method and device |
CN111025168A (en) * | 2019-11-29 | 2020-04-17 | 淮南师范学院 | Battery health state monitoring device and battery state of charge intelligent estimation method |
CN111864282B (en) * | 2020-07-28 | 2021-10-22 | 安徽江淮汽车集团股份有限公司 | Remaining power correction method, power automobile and readable storage medium |
CN111983479B (en) * | 2020-08-04 | 2021-05-04 | 珠海迈巨微电子有限责任公司 | Real-time establishing method and updating method of battery physical model and battery monitoring equipment |
CN112583066A (en) * | 2020-09-01 | 2021-03-30 | 骆驼集团新能源电池有限公司 | Forklift lithium iron phosphate battery charging method |
CN112265472A (en) * | 2020-10-29 | 2021-01-26 | 长城汽车股份有限公司 | Method and device for determining state of charge value of power battery of vehicle |
TWI790872B (en) * | 2021-12-23 | 2023-01-21 | 亞福儲能股份有限公司 | Battery management system and battery management method |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140244225A1 (en) * | 2013-02-24 | 2014-08-28 | The University Of Connecticut | Battery state of charge tracking, equivalent circuit selection and benchmarking |
-
2013
- 2013-06-20 CN CN201310245252.6A patent/CN104237791A/en active Pending
-
2014
- 2014-02-27 US US14/192,867 patent/US9709635B2/en not_active Expired - Fee Related
-
2017
- 2017-05-03 US US15/586,062 patent/US10175302B2/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140244225A1 (en) * | 2013-02-24 | 2014-08-28 | The University Of Connecticut | Battery state of charge tracking, equivalent circuit selection and benchmarking |
Cited By (60)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150046107A1 (en) * | 2012-10-26 | 2015-02-12 | Lg Chem, Ltd. | Apparatus and method for estimating state of charge of battery |
US9310441B2 (en) * | 2012-10-26 | 2016-04-12 | Lg Chem, Ltd. | Apparatus and method for estimating stage of charge of battery |
US9446678B2 (en) * | 2014-03-05 | 2016-09-20 | Ford Global Technologies, Llc | Battery model with robustness to cloud-specific communication issues |
US20150251556A1 (en) * | 2014-03-05 | 2015-09-10 | Ford Global Technologies, Llc | Battery model with robustness to cloud-specific communication issues |
US9533597B2 (en) | 2014-03-05 | 2017-01-03 | Ford Global Technologies, Llc | Parameter identification offloading using cloud computing resources |
CN104502858A (en) * | 2014-12-31 | 2015-04-08 | 桂林电子科技大学 | Power battery SOC estimation method based on backward difference discrete model and system thereof |
CN105045971A (en) * | 2015-07-01 | 2015-11-11 | 西安交通大学 | Lithium battery fractional order discretization impedance model |
CN104965179A (en) * | 2015-07-06 | 2015-10-07 | 首都师范大学 | Lithium ion storage battery temperature combinational circuit model and parameter identification method thereof |
CN105044440A (en) * | 2015-08-26 | 2015-11-11 | 苏州弗尔赛能源科技股份有限公司 | Fuel cell monolithic voltage inspection system based on LTC6803 |
CN105093128A (en) * | 2015-08-31 | 2015-11-25 | 山东智洋电气股份有限公司 | Storage battery state of charge (SOC) estimation method based on extended Kalman filtering (EKF) |
US10518640B2 (en) | 2015-11-17 | 2019-12-31 | Omron Corporation | Battery remaining capacity display device, battery system, and battery remaining capacity display method |
US10684328B2 (en) * | 2015-12-01 | 2020-06-16 | Omron Corporation | Remaining battery charge estimation system and remaining battery charge estimation method |
US20190041467A1 (en) * | 2015-12-01 | 2019-02-07 | Omron Corporation | Remaining battery charge estimation system and remaining battery charge estimation method |
CN105699907A (en) * | 2016-01-28 | 2016-06-22 | 广州市香港科大霍英东研究院 | A battery SOC estimation method and system based on dynamic impedance correction |
CN105738817A (en) * | 2016-01-29 | 2016-07-06 | 深圳市沃特玛电池有限公司 | Battery charge state estimation method based on AEKF and estimation system |
CN107643488A (en) * | 2016-07-21 | 2018-01-30 | 神讯电脑(昆山)有限公司 | Corresponding to the metering method and its electronic installation of the battery electric quantity of temperature |
US20180147952A1 (en) * | 2016-11-25 | 2018-05-31 | Hyundai Motor Company | Method and system for controlling motors |
US10189369B2 (en) * | 2016-11-25 | 2019-01-29 | Hyundai Motor Company | Method and system for controlling motors |
WO2018127296A1 (en) * | 2017-01-09 | 2018-07-12 | Volvo Truck Corporation | A method and arrangement for determining the state of charge of a battery pack |
CN110167783A (en) * | 2017-01-09 | 2019-08-23 | 沃尔沃卡车集团 | A kind of method and apparatus for determining the charged state of battery pack |
US11214150B2 (en) | 2017-01-09 | 2022-01-04 | Volvo Truck Corporation | Method and arrangement for determining the state of charge of a battery pack |
EP3565731B1 (en) | 2017-01-09 | 2022-10-19 | Volvo Truck Corporation | A method and arrangement for determining the state of charge of a battery pack |
CN106970328A (en) * | 2017-01-17 | 2017-07-21 | 深圳市沛城电子科技有限公司 | A kind of SOC estimation method and device |
JP2018151175A (en) * | 2017-03-10 | 2018-09-27 | 株式会社デンソーテン | Estimation device, estimation method, and estimation program |
JP7000030B2 (en) | 2017-03-10 | 2022-01-19 | 株式会社デンソーテン | Estimator, estimation method, and estimation program |
CN107390127A (en) * | 2017-07-11 | 2017-11-24 | 欣旺达电动汽车电池有限公司 | A kind of SOC estimation method |
CN109901082A (en) * | 2017-12-08 | 2019-06-18 | 南京德朔实业有限公司 | Portable electric energy system and its measurement method |
US11502530B2 (en) * | 2017-12-26 | 2022-11-15 | Panasonic Intellectual Property Management Co., Ltd. | Battery management device, battery system, and vehicle power supply system for managing battery state of charge level when in non-use state |
CN108287316A (en) * | 2018-01-15 | 2018-07-17 | 厦门大学 | Accumulator method for estimating remaining capacity based on threshold spread Kalman Algorithm |
WO2019142550A1 (en) * | 2018-01-18 | 2019-07-25 | 日立オートモティブシステムズ株式会社 | Secondary battery system |
JP2019124612A (en) * | 2018-01-18 | 2019-07-25 | 日立オートモティブシステムズ株式会社 | Secondary battery system |
JP7016704B2 (en) | 2018-01-18 | 2022-02-07 | ビークルエナジージャパン株式会社 | Rechargeable battery system |
EP3730957A4 (en) * | 2018-02-20 | 2021-03-17 | Lg Chem, Ltd. | Charge capacity calculation device and method for energy storage system |
US11467217B2 (en) | 2018-02-20 | 2022-10-11 | Lg Energy Solution, Ltd. | Charge capacity calculation device and method for energy storage system |
CN110245366A (en) * | 2018-03-08 | 2019-09-17 | 华为技术有限公司 | Dynamic power consumption estimation method, apparatus and system |
CN109031138A (en) * | 2018-06-29 | 2018-12-18 | 上海科列新能源技术有限公司 | A kind of safety evaluation method and device of power battery |
CN109061505A (en) * | 2018-08-28 | 2018-12-21 | 淮阴工学院 | A kind of detection method of lithium battery SOH |
CN109188293A (en) * | 2018-11-08 | 2019-01-11 | 武汉理工大学 | Based on new breath EKF lithium ion battery SOC estimation method of the covariance with fading factor |
CN109726365A (en) * | 2018-12-05 | 2019-05-07 | 新奥数能科技有限公司 | A kind of method and apparatus of load forecast |
CN112119317A (en) * | 2019-01-23 | 2020-12-22 | 株式会社Lg化学 | Battery management device, battery management method and battery pack |
US11069926B1 (en) * | 2019-02-14 | 2021-07-20 | Vcritonc Alpha, Inc. | Controlling ongoing battery system usage via parametric linear approximation |
US20210199724A1 (en) * | 2019-02-22 | 2021-07-01 | Lg Chem, Ltd. | Battery management system, battery management method, battery pack and electric vehicle |
US11567137B2 (en) * | 2019-02-22 | 2023-01-31 | Lg Energy Solution, Ltd. | Battery management system, battery management method, battery pack and electric vehicle |
CN110275113A (en) * | 2019-06-25 | 2019-09-24 | 内蒙古工业大学 | A kind of lithium battery charge state estimation method |
CN110361652A (en) * | 2019-06-26 | 2019-10-22 | 河南理工大学 | A kind of Kalman filtering lithium battery SOC estimation method based on Model Parameter Optimization |
CN110221219A (en) * | 2019-07-03 | 2019-09-10 | 中国民用航空飞行学院 | Airborne circumstance is got off the plane lithium battery SOC estimation method |
CN112255545A (en) * | 2019-07-05 | 2021-01-22 | 西南科技大学 | Lithium battery SOC estimation model based on square root extended Kalman filter |
CN110414117A (en) * | 2019-07-23 | 2019-11-05 | 北京航空航天大学 | A kind of soft bag lithium ionic cell sealed reliable degree prediction technique |
CN110441702A (en) * | 2019-07-31 | 2019-11-12 | 湘潭大学 | A method of lithium ion battery charge capacity is estimated with Extended Kalman filter |
CN110837622A (en) * | 2019-11-26 | 2020-02-25 | 国网湖南省电力有限公司 | Lithium battery state of charge estimation method based on high-rate discharge |
CN111002834A (en) * | 2019-12-18 | 2020-04-14 | 中车株洲电力机车有限公司 | Method and system for predicting cruising mileage of residual electric quantity of traction battery of railway vehicle |
CN113030555A (en) * | 2020-03-18 | 2021-06-25 | 深圳大学 | Energy storage open-circuit voltage estimation method and device, terminal equipment and storage medium |
CN112347415A (en) * | 2020-10-02 | 2021-02-09 | 广东电网有限责任公司广州供电局 | Prediction method based on short-circuit current zero crossing point prediction system |
CN112858929A (en) * | 2021-03-16 | 2021-05-28 | 上海理工大学 | Battery SOC estimation method based on fuzzy logic and extended Kalman filtering |
CN112881928A (en) * | 2021-03-24 | 2021-06-01 | 东风汽车集团股份有限公司 | Screening method for consistency of battery monomers |
CN113343443A (en) * | 2021-05-24 | 2021-09-03 | 国网综合能源服务集团有限公司 | Power distribution method for lithium battery prefabricated cabins with different SOC (state of charge) |
CN113341330A (en) * | 2021-05-25 | 2021-09-03 | 西南大学 | Lithium-sulfur power battery SOC estimation method based on OCV correction and Kalman filtering algorithm |
CN113900031A (en) * | 2021-09-14 | 2022-01-07 | 国网浙江省电力有限公司电力科学研究院 | SOC safety verification method after energy storage system is accessed to AGC |
CN114397581A (en) * | 2021-12-09 | 2022-04-26 | 国网天津市电力公司 | New energy automobile battery SOC anti-interference assessment method oriented to charging monitoring data of direct-current charging pile |
CN117811169A (en) * | 2024-02-29 | 2024-04-02 | 威海凯瑞电气股份有限公司 | Battery active equalization management system |
Also Published As
Publication number | Publication date |
---|---|
CN104237791A (en) | 2014-12-24 |
US20170234934A1 (en) | 2017-08-17 |
US10175302B2 (en) | 2019-01-08 |
US9709635B2 (en) | 2017-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10175302B2 (en) | Power system and state of charge estimation | |
CN110914696B (en) | Method and system for estimating battery open cell voltage, state of charge, and state of health during operation of a battery | |
JP4767558B2 (en) | Power supply state detection device, power supply device, and initial characteristic extraction device used for power supply device | |
He et al. | Online model-based estimation of state-of-charge and open-circuit voltage of lithium-ion batteries in electric vehicles | |
CN107368619B (en) | Extended Kalman filtering SOC estimation method | |
TWI384246B (en) | Apparatus and method for estimating resistance characteristics of battery based on open circuit voltage estimated by battery voltage variation | |
CN102662148B (en) | On-line feedback battery state of charge (SOC) predicting method | |
EP2321663B1 (en) | Apparatus and method for estimating state of health of battery based on battery voltage variation pattern | |
Shahriari et al. | Online state-of-health estimation of VRLA batteries using state of charge | |
EP3018753B1 (en) | Battery control method based on ageing-adaptive operation window | |
Chen et al. | Battery state of charge estimation based on a combined model of Extended Kalman Filter and neural networks | |
US8751086B2 (en) | Online battery capacity estimation | |
CN106814329A (en) | A kind of battery SOC On-line Estimation method based on double Kalman filtering algorithms | |
US20160363630A1 (en) | Systems and methods for estimating battery system parameters | |
JP6509725B2 (en) | Estimating the state of charge of the battery | |
JP7211420B2 (en) | Parameter estimation device, parameter estimation method and computer program | |
CN111929596A (en) | Method and device for acquiring battery capacity, storage medium and electronic equipment | |
EP4212896A1 (en) | Method for estimating state of charge of battery | |
KR20120046355A (en) | Apparatus and method for reporting exchange time of battery | |
Lu et al. | Modeling discharge characteristics for predicting battery remaining life | |
US20220196754A1 (en) | Method for detecting abnormal battery cell | |
Feng et al. | An improved lithium-ion battery model with temperature prediction considering entropy | |
CN111856286A (en) | DP-RC model-based battery power estimation method and device | |
WO2016054732A1 (en) | Method and system for estimating instantaneous state-of-charge of a lithium ion battery | |
CN109782177A (en) | A kind of acquisition methods of battery capacity, device and automobile |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: UNIVERSITY OF ELECTRONIC SCIENCE AND TECHNOLOGY OF Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ZHONG, QISHUI;LI, BAIHUA;LI, HUI;AND OTHERS;REEL/FRAME:032367/0480 Effective date: 20140124 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20210718 |