CA2585921C - State and parameter estimation for an electrochemical cell - Google Patents
State and parameter estimation for an electrochemical cell Download PDFInfo
- Publication number
- CA2585921C CA2585921C CA2585921A CA2585921A CA2585921C CA 2585921 C CA2585921 C CA 2585921C CA 2585921 A CA2585921 A CA 2585921A CA 2585921 A CA2585921 A CA 2585921A CA 2585921 C CA2585921 C CA 2585921C
- Authority
- CA
- Canada
- Prior art keywords
- prediction
- state
- parameter
- component configured
- internal
- 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.)
- Expired - Lifetime
Links
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/42—Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
- H01M10/48—Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/60—Methods or arrangements for performing computations using a digital non-denominational number representation, i.e. number representation without radix; Computing devices using combinations of denominational and non-denominational quantity representations, e.g. using difunction pulse trains, STEELE computers, phase computers
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- General Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Secondary Cells (AREA)
- Electrolytic Production Of Non-Metals, Compounds, Apparatuses Therefor (AREA)
- Fuel Cell (AREA)
Abstract
Methods and apparatus for estimating the state and parameters of an electrochemical cell. More particularly, for example, a method for estimating present states and present parameters of an electrochemical cell system comprising: estimating a state value of the electrochemical cell with a cell state filter to estimate the state value; estimating a parameter value of the electrochemical cell with a cell parameter filter to estimate the parameter value, and exchanging information between the cell state filter and the cell parameter filter. Also an apparatus configured to estimate present states and present parameters of an electrochemical cell comprising: a first component configured to estimate a cell state value; and a second component configured to estimate a cell parameter value. The first component and second component are in operable communication to exchange information there between.
Description
STATE AND PARAMETER ESTIMATION FOR AN ELECTROCHEMICAL CELL
Technical Field The present invention relates to methods and apparatus for estimation of battery pack system state and model parameters using digital filtering techniques. In particular, dual Kalman. filtering and dual extended Kalman filtering.
Background Art In the context of rechargeable battery pack technologies, it is desired in some applications to be able to estimate quantities that are descriptive of the present battery pack condition, but that may not be directly measured. Some of these quantities may change rapidly, such as the pack state-of-charge (SOC), which can traverse its entire range within minutes. Others may change very slowly, such as cell capacity, which might change as little as 20% in a decade or more of regular use. The quantities that tend to change quickly comprise the "state" of the system, and the quantities that tend to change slowly comprise the time varying "parameters" of the system.
In the context of the battery systems, particularly those that need to operate for long periods of time, as aggressively as possible without harming the battery life, for example, in Hybrid Electric Vehicles (HEVs), Battery Electric Vehicles (BEVs), laptop computer batteries, portable tool battery packs, and the like, it is desired that information regarding quickly varying parameters (e.g., SOC) be used to estimate how much battery energy is presently available to do work, and so forth. Further, it may be desirable to ascertain information regarding slowly varying parameters (e.g., total capacity) in order to keep the prior calculations precise over the lifetime of the pack, extending its useful service time, and help in determining the state-of-health (SOH) of the pack.
There are a number of existing methods for estimating the state of a cell, which are generally" concerned with estimating three quantities: SOC (a quickly varying quantity), power-fade, and capacity-fade (both slowly time varying).
Power fade may be calculated if the present and initial pack electrical resistances are known, and capacity fade may be calculated if present and initial pack total capacities are known, for example, although other methods may also be used.
Power- and capacity-fade are often lumped under the description "state-of-health" (SOH) . Some other information may be derived using the values of these variables, such as the maximum power available from the pack at any given time.
Additional state members or parameters may also be needed for specific applications, and individual algorithms would typically be required to find each one.
SOC is a value, typically reported in percent, that indicates the fraction of the cell capacity presently available to do work. A number of different approaches to estimating SOC have been- employed: a discharge test, ampere-hour counting (Coulomb counting), measuring the electrolyte, open-circuit voltage measurement, linear and nonlinear circuit modeling, impedance spectroscopy, measurement of internal resistance, coup de fouet, and some forms of Kalman filtering. The discharge test must completely discharge the cell in order to determine SOC. This test interrupts system
Technical Field The present invention relates to methods and apparatus for estimation of battery pack system state and model parameters using digital filtering techniques. In particular, dual Kalman. filtering and dual extended Kalman filtering.
Background Art In the context of rechargeable battery pack technologies, it is desired in some applications to be able to estimate quantities that are descriptive of the present battery pack condition, but that may not be directly measured. Some of these quantities may change rapidly, such as the pack state-of-charge (SOC), which can traverse its entire range within minutes. Others may change very slowly, such as cell capacity, which might change as little as 20% in a decade or more of regular use. The quantities that tend to change quickly comprise the "state" of the system, and the quantities that tend to change slowly comprise the time varying "parameters" of the system.
In the context of the battery systems, particularly those that need to operate for long periods of time, as aggressively as possible without harming the battery life, for example, in Hybrid Electric Vehicles (HEVs), Battery Electric Vehicles (BEVs), laptop computer batteries, portable tool battery packs, and the like, it is desired that information regarding quickly varying parameters (e.g., SOC) be used to estimate how much battery energy is presently available to do work, and so forth. Further, it may be desirable to ascertain information regarding slowly varying parameters (e.g., total capacity) in order to keep the prior calculations precise over the lifetime of the pack, extending its useful service time, and help in determining the state-of-health (SOH) of the pack.
There are a number of existing methods for estimating the state of a cell, which are generally" concerned with estimating three quantities: SOC (a quickly varying quantity), power-fade, and capacity-fade (both slowly time varying).
Power fade may be calculated if the present and initial pack electrical resistances are known, and capacity fade may be calculated if present and initial pack total capacities are known, for example, although other methods may also be used.
Power- and capacity-fade are often lumped under the description "state-of-health" (SOH) . Some other information may be derived using the values of these variables, such as the maximum power available from the pack at any given time.
Additional state members or parameters may also be needed for specific applications, and individual algorithms would typically be required to find each one.
SOC is a value, typically reported in percent, that indicates the fraction of the cell capacity presently available to do work. A number of different approaches to estimating SOC have been- employed: a discharge test, ampere-hour counting (Coulomb counting), measuring the electrolyte, open-circuit voltage measurement, linear and nonlinear circuit modeling, impedance spectroscopy, measurement of internal resistance, coup de fouet, and some forms of Kalman filtering. The discharge test must completely discharge the cell in order to determine SOC. This test interrupts system
2
3 PCT/KR2004/003103 function while the test is being performed and can be overly time consuming rendering it not useful for many applications.
Ampere-hour counting (Coulomb counting) is an "open loop"
methodology whose accuracy degrades over time by accumulated measurement error. Measuring the electrolyte is only feasible for vented lead-acid batteries, and therefore has limited applicability. Open-circuit voltage measurement may be performed only after extended periods of cell inactivity, and for cells with negligible hysteresis effect and does not work in a dynamic setting. Linear and nonlinear circuit modeling methods do not yield SOC directly; SOC must be inferred from the calculated values. Impedance spectroscopy requires making measurements not always available in a general application. Measurement of internal resistance is very sensitive to measurement error, and requires measurements not available in general applications. Coup de fouet works for lead-acid batteries only. Forms of Kalman filtering that do not use SOC as a filter state do not directly yield error bounds on the estimate. In another method, described in U.S.
Patent No. 6,534,954, a filter, preferably a Kalman filter is used to estimate SOC by employing a known mathematical model of cell dynamics and measurements of cell voltage, current, and temperature. This method directly estimates state values.
However, it does not address parameter values.
Not only is knowledge of SOC desired, but also knowledge of SOH. In this context, power fade refers to the phenomenon of increasing cell electrical resistance as the cell ages.
This increasing resistance causes the power that can be sourced/sunk by the cell to drop. Capacity fade refers to the phenomenon of decreasing cell total capacity as the cell ages. Both the cell's resistance and capacity are time-varying parameters. The prior art uses the following different approaches to estimate SOH: the discharge test, chemistry-dependent methods, Ohmic tests, and partial discharge. The discharge test completely discharges a fully charged cell in order to determine its total capacity. This test interrupts system function and wastes cell energy.
Chemistry-dependent methods include measuring the level of plate corrosion, electrolyte density, and "coup de fouet" for lead-acid batteries. Ohmic tests include resistance, conductance and impedance tests, perhaps combined with fuzzy-logic algorithms and/or neural networks. These methods require invasive measurements. Partial discharge and other methods compare cell-under-test to a good cell or model of a good cell.
There is a. need for a method to concurrently estimate the state and parameters of a cell. Furthermore, there is a need for tests that do not interrupt system function and do not waste energy, methods that are generally applicable (i.e., to different types of cell electrochemistries and to different applications), methods that do not require invasive measurements, and more rigorous approaches. There is a need for methods and apparatus for automatically estimating time-varying parameters, such as the cell's resistance and capacity. There is a need for a method that will work with different configurations of parallel and/or series cells in a battery pack.
Ampere-hour counting (Coulomb counting) is an "open loop"
methodology whose accuracy degrades over time by accumulated measurement error. Measuring the electrolyte is only feasible for vented lead-acid batteries, and therefore has limited applicability. Open-circuit voltage measurement may be performed only after extended periods of cell inactivity, and for cells with negligible hysteresis effect and does not work in a dynamic setting. Linear and nonlinear circuit modeling methods do not yield SOC directly; SOC must be inferred from the calculated values. Impedance spectroscopy requires making measurements not always available in a general application. Measurement of internal resistance is very sensitive to measurement error, and requires measurements not available in general applications. Coup de fouet works for lead-acid batteries only. Forms of Kalman filtering that do not use SOC as a filter state do not directly yield error bounds on the estimate. In another method, described in U.S.
Patent No. 6,534,954, a filter, preferably a Kalman filter is used to estimate SOC by employing a known mathematical model of cell dynamics and measurements of cell voltage, current, and temperature. This method directly estimates state values.
However, it does not address parameter values.
Not only is knowledge of SOC desired, but also knowledge of SOH. In this context, power fade refers to the phenomenon of increasing cell electrical resistance as the cell ages.
This increasing resistance causes the power that can be sourced/sunk by the cell to drop. Capacity fade refers to the phenomenon of decreasing cell total capacity as the cell ages. Both the cell's resistance and capacity are time-varying parameters. The prior art uses the following different approaches to estimate SOH: the discharge test, chemistry-dependent methods, Ohmic tests, and partial discharge. The discharge test completely discharges a fully charged cell in order to determine its total capacity. This test interrupts system function and wastes cell energy.
Chemistry-dependent methods include measuring the level of plate corrosion, electrolyte density, and "coup de fouet" for lead-acid batteries. Ohmic tests include resistance, conductance and impedance tests, perhaps combined with fuzzy-logic algorithms and/or neural networks. These methods require invasive measurements. Partial discharge and other methods compare cell-under-test to a good cell or model of a good cell.
There is a. need for a method to concurrently estimate the state and parameters of a cell. Furthermore, there is a need for tests that do not interrupt system function and do not waste energy, methods that are generally applicable (i.e., to different types of cell electrochemistries and to different applications), methods that do not require invasive measurements, and more rigorous approaches. There is a need for methods and apparatus for automatically estimating time-varying parameters, such as the cell's resistance and capacity. There is a need for a method that will work with different configurations of parallel and/or series cells in a battery pack.
4 Disclosure of the Invention Disclosed herein in one or more exemplary embodiments are methods and apparatus for estimating the state and parameters of an electrochemical cell. More particularly, for example, estimating state and parameter values of a cell.
A first aspect of the invention is a method for estimating present states and present parameters of an electrochemical cell system comprising: estimating a state value of the electrochemical cell with a cell state filter to estimate said state value; and estimating a parameter value of the electrochemical cell with a cell parameter filter to estimate said parameter value, wherein said estimating a state value comprises: making an internal state prediction of said cell; making an uncertainty prediction of said internal state prediction;
correcting said internal state prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal state prediction, said making an uncertainty prediction and said correcting to yield an ongoing estimation to said state and an ongoing uncertainty to said state estimating, and wherein said estimating a parameter value comprises: making an internal parameter prediction of said cell; making an uncertainty prediction of said internal parameter prediction; correcting said internal parameter prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal parameter prediction, said making an uncertainty prediction and
A first aspect of the invention is a method for estimating present states and present parameters of an electrochemical cell system comprising: estimating a state value of the electrochemical cell with a cell state filter to estimate said state value; and estimating a parameter value of the electrochemical cell with a cell parameter filter to estimate said parameter value, wherein said estimating a state value comprises: making an internal state prediction of said cell; making an uncertainty prediction of said internal state prediction;
correcting said internal state prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal state prediction, said making an uncertainty prediction and said correcting to yield an ongoing estimation to said state and an ongoing uncertainty to said state estimating, and wherein said estimating a parameter value comprises: making an internal parameter prediction of said cell; making an uncertainty prediction of said internal parameter prediction; correcting said internal parameter prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal parameter prediction, said making an uncertainty prediction and
5 said correcting to yield an ongoing estimation to said parameter and an ongoing uncertainty to said parameter estimation.
Another aspect is an apparatus configured to estimate present states and present parameters of an electrochemical cell comprising: a first component configured to estimate a cell state value; and a second component configured to estimate a cell parameter value, wherein said first component is configured to estimate a state value comprises: a component configured to make an internal state prediction of said cell; a component configured to make an uncertainty prediction of said internal state prediction; a component configured to correct said internal state prediction and said uncertainty prediction; and a component configured to apply an algorithm that iterates steps taken by said component configured to make an internal state prediction, said component .configured to make an uncertainty prediction and said component configured to correct to yield an ongoing estimation to said state and an ongoing uncertainty to said state estimating, and wherein said second component is configured to estimate a parameter value comprises: a component configured to make an internal parameter prediction of said cell; a component configured to make an uncertainty prediction of said internal parameter prediction; a component configured to correct said internal parameter prediction and said uncertainty prediction; and a component configured to apply an algorithm that iterates steps
Another aspect is an apparatus configured to estimate present states and present parameters of an electrochemical cell comprising: a first component configured to estimate a cell state value; and a second component configured to estimate a cell parameter value, wherein said first component is configured to estimate a state value comprises: a component configured to make an internal state prediction of said cell; a component configured to make an uncertainty prediction of said internal state prediction; a component configured to correct said internal state prediction and said uncertainty prediction; and a component configured to apply an algorithm that iterates steps taken by said component configured to make an internal state prediction, said component .configured to make an uncertainty prediction and said component configured to correct to yield an ongoing estimation to said state and an ongoing uncertainty to said state estimating, and wherein said second component is configured to estimate a parameter value comprises: a component configured to make an internal parameter prediction of said cell; a component configured to make an uncertainty prediction of said internal parameter prediction; a component configured to correct said internal parameter prediction and said uncertainty prediction; and a component configured to apply an algorithm that iterates steps
6 taken by said component configured to make an internal parameter prediction, said component configured to make an uncertainty prediction and said component configured to correct to yield an ongoing estimation to said parameter and an ongoing uncertainty to said parameter estimating.
Further, disclosed herein in another exemplary embodiment is a storage medium encoded with a machine-readable computer program code, wherein said storage medium includes instructions for causing a computer to implement a method for estimating present states and present parameters of an electrochemical cell comprising: estimating a state value of the electrochemical cell with a cell state filter to estimate said state value; and estimating a parameter value of the electrochemical cell with a cell parameter filter to estimate said parameter value, wherein said estimating a state value comprises: making an internal state prediction of said cell; making an uncertainty prediction of said internal state prediction; correcting said internal state prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal state prediction, said making an uncertainty prediction and said correcting to yield an ongoing estimation to said state and an ongoing uncertainty to said state estimating, and wherein said estimating a parameter value comprises: making an internal parameter prediction of said cell; making an uncertainty prediction of said internal parameter 6a prediction; correcting said internal parameter prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal parameter prediction, said making an uncertainty prediction and said correcting to yield an ongoing estimation to said parameter and an ongoing uncertainty to said parameter estimation.
Brief Description of the Drawings These and other features, aspects and advantages of the present invention will become better understood with regard to the following description, appended claims and accompanying drawings wherein like elements are numbered alike in the several Figures:
FIGURE 1 is a block diagram illustrating an exemplary system for state and parameter estimation in accordance with an exemplary embodiment of the invention;
and FIGURE 2 is a block diagram depicting a method of dual filtering, in accordance with an exemplary embodiment of the invention.
Best Mode For Carrying Out The Invention Disclosed herein are various embodiments and methods, systems and apparatus for the estimation of states and parameters of an electrochemical cell using dual filtering. Referring now to Figures 1 and 2, in the following description, numerous specific details are set forth in order to provide a more complete understanding of 6b the present invention. It will be appreciated that while the exemplary embodiments are described with reference to a battery cell, numerous electrochemical cells hereinafter referred to as a cell, may be employed, including, but not limited to, batteries, battery packs, ultracapacitors, capacitor banks, fuel cells, electrolysis cells, and the like, as well as combinations including at least one of the foregoing. Furthermore, it will be appreciated that a battery or battery pack may include a plurality of cells, where the exemplary embodiments disclosed herein are applied to one or more cells of the plurality.
One or more exemplary embodiments of the present invention estimate cell state and parameter values using dual filtering. One or more exemplary embodiments of the present invention estimate cell state and parameter values using dual Kalman filtering. Some embodiments of the present invention estimate cell state and parameter values using dual extended Kalman filtering. Some embodiments simultaneously estimate SOC, power- and/or capacity-fade, while others estimate additional cell state values and/or additional time-varying parameter values. It will further be appreciated that while the term filtering is employed for description and illustration of the exemplary embodiments, the terminology is intended to include methodologies of recursive prediction and correction commonly denoted as filtering, including but not limited to Kalman filtering and/or extended Kalman filtering.
FIG. 1 shows the components of the state and parameter estimator system 10 according an embodiment of the present invention. Electrochemical cell pack 20 comprising a plurality of cells 22, e.g., battery is connected to a load
Further, disclosed herein in another exemplary embodiment is a storage medium encoded with a machine-readable computer program code, wherein said storage medium includes instructions for causing a computer to implement a method for estimating present states and present parameters of an electrochemical cell comprising: estimating a state value of the electrochemical cell with a cell state filter to estimate said state value; and estimating a parameter value of the electrochemical cell with a cell parameter filter to estimate said parameter value, wherein said estimating a state value comprises: making an internal state prediction of said cell; making an uncertainty prediction of said internal state prediction; correcting said internal state prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal state prediction, said making an uncertainty prediction and said correcting to yield an ongoing estimation to said state and an ongoing uncertainty to said state estimating, and wherein said estimating a parameter value comprises: making an internal parameter prediction of said cell; making an uncertainty prediction of said internal parameter 6a prediction; correcting said internal parameter prediction and said uncertainty prediction; and applying an algorithm that iterates said making an internal parameter prediction, said making an uncertainty prediction and said correcting to yield an ongoing estimation to said parameter and an ongoing uncertainty to said parameter estimation.
Brief Description of the Drawings These and other features, aspects and advantages of the present invention will become better understood with regard to the following description, appended claims and accompanying drawings wherein like elements are numbered alike in the several Figures:
FIGURE 1 is a block diagram illustrating an exemplary system for state and parameter estimation in accordance with an exemplary embodiment of the invention;
and FIGURE 2 is a block diagram depicting a method of dual filtering, in accordance with an exemplary embodiment of the invention.
Best Mode For Carrying Out The Invention Disclosed herein are various embodiments and methods, systems and apparatus for the estimation of states and parameters of an electrochemical cell using dual filtering. Referring now to Figures 1 and 2, in the following description, numerous specific details are set forth in order to provide a more complete understanding of 6b the present invention. It will be appreciated that while the exemplary embodiments are described with reference to a battery cell, numerous electrochemical cells hereinafter referred to as a cell, may be employed, including, but not limited to, batteries, battery packs, ultracapacitors, capacitor banks, fuel cells, electrolysis cells, and the like, as well as combinations including at least one of the foregoing. Furthermore, it will be appreciated that a battery or battery pack may include a plurality of cells, where the exemplary embodiments disclosed herein are applied to one or more cells of the plurality.
One or more exemplary embodiments of the present invention estimate cell state and parameter values using dual filtering. One or more exemplary embodiments of the present invention estimate cell state and parameter values using dual Kalman filtering. Some embodiments of the present invention estimate cell state and parameter values using dual extended Kalman filtering. Some embodiments simultaneously estimate SOC, power- and/or capacity-fade, while others estimate additional cell state values and/or additional time-varying parameter values. It will further be appreciated that while the term filtering is employed for description and illustration of the exemplary embodiments, the terminology is intended to include methodologies of recursive prediction and correction commonly denoted as filtering, including but not limited to Kalman filtering and/or extended Kalman filtering.
FIG. 1 shows the components of the state and parameter estimator system 10 according an embodiment of the present invention. Electrochemical cell pack 20 comprising a plurality of cells 22, e.g., battery is connected to a load
7 circuit 30. For example, load circuit 30 could be a motor in an Electric Vehicle (EV) or a Hybrid Electric Vehicle (HEV).
An apparatus for measuring various cell characteristics and properties is provided as 40. The measurement apparatus 40 may include but not be limited to a device for measurement of cell terminal voltage such as a voltage sensor 42, e.g. a voltmeter and the like, while measurements of cell current are made with a current sensing device 44, e.g., an ammeter and the like. Optionally, measurements of cell temperature are made with a temperature sensor 46, e.g., a thermometer and the like. Pressure sensors and/or impedance sensors 48 are also possible and may be employed for selected types of cells. Various sensors may be employed as needed to evaluate the characteristics and properties of the cell(s). Voltage, current, and optionally temperature measurements are processed with an arithmetic circuit 50, e.g., processor or computer, which estimates the states and parameters of the cell(s). The system may also include a storage medium 52 comprising any computer usable storage medium known to one of ordinary skill in the art. The storage medium is in operable communication with arithmetic circuit 50 employing various means, including, but not limited to a propagated signal 54.
It should be appreciated that no instrument is required to take measurements from the internal chemical components of the cell 12 although such instrumentation may be used with this invention. Also note that all measurements may be non-invasive; that is, no signal must be injected into the system that might interfere with the proper operation of load circuit 30.
An apparatus for measuring various cell characteristics and properties is provided as 40. The measurement apparatus 40 may include but not be limited to a device for measurement of cell terminal voltage such as a voltage sensor 42, e.g. a voltmeter and the like, while measurements of cell current are made with a current sensing device 44, e.g., an ammeter and the like. Optionally, measurements of cell temperature are made with a temperature sensor 46, e.g., a thermometer and the like. Pressure sensors and/or impedance sensors 48 are also possible and may be employed for selected types of cells. Various sensors may be employed as needed to evaluate the characteristics and properties of the cell(s). Voltage, current, and optionally temperature measurements are processed with an arithmetic circuit 50, e.g., processor or computer, which estimates the states and parameters of the cell(s). The system may also include a storage medium 52 comprising any computer usable storage medium known to one of ordinary skill in the art. The storage medium is in operable communication with arithmetic circuit 50 employing various means, including, but not limited to a propagated signal 54.
It should be appreciated that no instrument is required to take measurements from the internal chemical components of the cell 12 although such instrumentation may be used with this invention. Also note that all measurements may be non-invasive; that is, no signal must be injected into the system that might interfere with the proper operation of load circuit 30.
8 In order to perform the prescribed functions and desired processing, as well as the computations therefore (e.g., the modeling, estimation of states and parameters prescribed herein, and the like), arithmetic circuit 50 may include, but not be limited to, a processor(s), gate array(s), custom logic, computer(s), memory, storage, register(s), timing, interrupt(s), communication interfaces, and input/output signal interfaces, as well as combinations comprising at least one of the foregoing. Arithmetic circuit 50 may also include inputs and input signal filtering and the like, to enable accurate sampling and conversion or acquisitions of signals from communications interfaces and inputs.
Additional features of arithmetic circuit 50 and certain processes therein are thoroughly discussed at a later point herein.
One or more embodiments of the invention may be implemented as new or updated firmware and software executed in arithmetic circuit 50 and/or other processing controllers.
Software functions include, but are not limited to firmware and may be implemented in hardware, software, or a combination thereof. Thus a distinct advantage of the present invention is that it may be implemented for use with existing and/or new processing systems for electrochemical cell charging and control.
In an exemplary embodiment, Arithmetic circuit 50 uses a mathematical model of the cell 12 that includes indicia of a dynamic system state, including, but not limited to, the SOC
as a model state. In one embodiment of the present invention, a discrete-time model is used. An exemplary model in a
Additional features of arithmetic circuit 50 and certain processes therein are thoroughly discussed at a later point herein.
One or more embodiments of the invention may be implemented as new or updated firmware and software executed in arithmetic circuit 50 and/or other processing controllers.
Software functions include, but are not limited to firmware and may be implemented in hardware, software, or a combination thereof. Thus a distinct advantage of the present invention is that it may be implemented for use with existing and/or new processing systems for electrochemical cell charging and control.
In an exemplary embodiment, Arithmetic circuit 50 uses a mathematical model of the cell 12 that includes indicia of a dynamic system state, including, but not limited to, the SOC
as a model state. In one embodiment of the present invention, a discrete-time model is used. An exemplary model in a
9 (possibly nonlinear) discrete-time state-space form has the form:
xk+l =f(xk,uk,ek) +Wk (1) Yk -g(xk,uk, k)+vk, where xk is the system state, dk is the set of time varying model parameters, ukis the exogenous input, ykis the system output, and wk and vkare "noise" inputs-all quantities may be scalars or vectors. f(-,=,=) and g(=,=;) are functions defined by the cell model being used. Non-time-varying numeric values required by the model may be embedded within f(=,,=) and g(=,=,=) , and are not included in k.
The system state includes, at least, a minimum amount of information, together with the present input and a mathematical model of the cell, needed to predict the present output. For a cell 12, the state might include: SOC, polarization voltage levels with respect to different time constants, and hysteresis levels, for example. The system exogenous input uk includes at minimum the present cell current ik , and may, optionally, include cell temperature (unless temperature change is itself modeled in the state) .
The system parameters k are the values that change only slowly with time, in such a way that they may not be directly determined with knowledge of the system measured input and output. These might include: cell capacity, resistance, polarization voltage time constant(s), polarization voltage blending factor(s), hysteresis blending factor(s), hysteresis rate constant(s), efficiency factor(s), and so forth. The model output Yk corresponds to physically measurable cell quantities or those directly computable from measured quantities at minimum for example, the cell voltage under load.
A mathematical model of parameter dynamics is also utilized. An exemplary model has the form:
k+1 = k + r k ( 2 dk - g(xkIukIek) +ek.
The first equation states that the parameters are essentially constant, but that they may change slowly over time, in this instance, modeled by a fictitious "noise"
process denoted, rk . The "output" of the optimum parameter dynamics is the cell output estimate plus some estimation error ek .
With models of the system state dynamics and parameter dynamics defined, in an exemplary embodiment, a procedure of dual filtering is applied. Once again, alternatively, dual Kalman filters may be employed, or dual extended Kalman filters. Furthermore, combinations of the abovementioned may also be employed. Table 1 is an exemplary implementation of the methodology and system utilizing dual extended Kalman filtering. The procedure is initialized by setting the parameter estimate to the best guess of the true parameters, e.g., 6 = E[B0] , and by setting the state estimate x to the best estimate of the cell state, e.g., x=E[x,] . The estimation-error covariance matrices EC and Et are also initialized. For example, an initialization of SOC might be estimated/based on a cell voltage in a look-up table, or information that was previously stored when a battery pack/cell was last powered down. Other examples might incorporate the length of time that the battery system had rested since powerdown and the like.
Table 1: Dual extended Kalman filter for state and weight update.
State-space models:
xk+1 f(xk,uk,Ok)+Wkand k+1 - k+1k vk g(xkIukIek)+ Vk dk _g(xkIukIek)+ek where wk , Vk , rk and ek are independent, zero-mean, Gaussian noise processes of covariance matrices E,, , Z, , E, and E, respectively.
Definitions:
_ af(xk-I,uk-1,lb x_ ag(xk,uk, k) dg(xk,u ,8) Ak-1 ax Ck ax Ck d k-1 xk_t_xk I k Xk_Xk 9=BA.
Initialization: For k=0, set o = E[90 ], EB O = E[(Bo - eo )( o - eo )T
xo = E[xo ], E C,o = E[(xo - $ )(xo - xo )T
Computation: For k =1,2,..., compute :
Time update for the weight filter [I corrected the following equation]
Bk - 0111--1 k 0,k-1 Time update for the state filter xk _flxk-l~uk-11ek) + T
EY k = Ak-1ZC k-1Ak-1 +Ew.
Measurement update for the state filter Lk = EX,k (Ck )T LCk Ex,k (Ck )T +E, xk -xk +Lk[Yk -g(xk,ukIek )]
Ex,k -(I-LkCk)Ex,k' Measurement update for the weight filter Lk - Ek (Ck) [Ck E/''g,k (Ck) +Ee Bk = /Bk + Lk [Yk - g(xk , uk , ek)}
Ek - (I LkCk )Ee In this example, several steps are performed in each measurement interval. First, the previous parameter estimate is propagated forward in time. The new parameter estimate is equal to the old parameter estimate, and the parameter error uncertainty is larger due to the passage of time (accommodated for in the model by the fictitious driving noise rk ). Various possibilities exist for updating the parameter uncertainty estimate-the table gives only one example. The state estimate and its uncertainty are propagated forward one step in time. A measurement of the cell output is made, and compared to the predicted output based on the state estimate, i and parameter estimate, 0 the difference is used to update the values of i and B. It may readily be appreciated that the steps outlined in the table may be performed in a variety of orders. While the table lists an exemplary ordering for the purposes of illustration, those skilled in the art will be able to identify many equivalent ordered sets of equations.
Turning now to FIG. 2 as well, an exemplary implementation of an exemplary embodiment of the invention is depicted. Two filters run in parallel are depicted. One filter 101/102 adapting the state estimate, x and one filter 103/104 adapting the parameter estimate, B . Additionally, information exchange between the filters permits updates from the parameter estimation filter 103/104 to the state estimation filter 101/102. Both filters 101/102, 103/104 have a time update or prediction aspect and a measurement update or correction aspect. State time update /prediction block 101 receives as input the previous exogenous input uk_, (which might include cell current and/or temperature, for example) along with the previously estimated system state value xk, and state uncertainty estimate Eck_1, and present predicted parameters Bk . The state time update/prediction block 101 provides predicted state 2k and predicted state uncertainty irk output to state measurement update/correction block 102, and to parameter measurement update/correction block 104. State measurement update/correction block 102 provides current system state estimate 241 and state uncertainty estimate Ck Parameter time update/prediction block 103 receives as input the previous time varying parameters estimate O and parameter uncertainty estimate Z+k_i . Parameter time update/prediction block 103 outputs predicted parameters k and predicted parameter uncertainty Z- to the parameter measurement update/correction block 104 and state time update block 101. Parameter measurement update block 104, which provides current parameter estimate ek and parameter uncertainty estimate Y-e<<. It will also be appreciated that a minus notation denotes that the vector is the result of the prediction components 101, 103 of the filters, while the plus notation denotes that the vector is the result of the correction component 102/104 of the filters.
Embodiments of this invention require a mathematical model of cell state and output dynamics for the particular application. In the examples above, this is accomplished by defining specific functions for f(;;) and g(=,=;) . An exemplary embodiment uses a cell model that includes effects due to one or more of the open-circuit-voltage (OCV) for the cell, internal resistance, voltage polarization time constants, and a hysteresis level. For the purpose of example, parameter values are fitted to this model structure to model the dynamics of high-power Lithium-Ion Polymer Battery (LiPB) cells, although the structure and methods presented here are general and apply to other electrochemistries.
In this example, SOC is captured by one state of the model. This equation is zk[m+1]=zk[m]-(r1At/Ck)Zk[m] (3) where At represents the inter-sample period (in seconds), Ck represents the capacity of cell number k in the pack (in ampere-seconds) , zk[m] is the SOC of cell k at time index m , ik , is the current out of cell k , and 77, is the Coulombic efficiency of a cell at current level ik.
In this example, the polarization voltage levels are captured by several filter states. If we let there be nf polarization voltage time constants, then Jk[in+1] = Affk[m]+Bfik[m]. (4) The matrix Af eRi."f'"f may be a diagonal matrix with real-valued polarization voltage time constants a, ===a"f . If so, the system is stable if all entries have magnitude less than one. The vector B f E R"f=X' may simply be set to n f "1"s . The entries of Bf are not critical as long as they are non-zero. The value of nfentries in the Afmatrix are chosen as part of the system identification procedure to best fit the model parameters to measured cell data. The Af and Bf matrices may vary with time and other factors pertinent to the present battery pack operating condition.
In this example, the hysteresis level is captured by a single state hk[m+1]=exp - 77;Zk[C]YAt) hk[rn]+ 1 -exp - 77rik[m]rAt sg(ik[m]), (5) where y is' the hysteresis rate constant, again found by system identification.
In this example, the overall model state is xk[m]= [fk[m]T hk[nz] zk[m]]T, where other orderings of states are possible. The state equation for the model is formed by combining all of the individual equations identified above.
In this example, the output equation that combines the state values to predict cell voltage is Vk [m] = OCV(zk 1M])+ G.fk [m] - Rik [m] + Mhk [m], (6) where Ge1i'xn, is a vector of polarization voltage blending factors g,===g,,, that blend the polarization voltage states together in the output, R is the cell resistance (different values may be used for discharge/charge), and M is the hysteresis blending factor. Note, G may be constrained such that the dc-gain from ikto Gfk is zero.
In this example, the parameters are T
0=[z7, C, a,...anj, 91...gnri, y, R, MI = (7) In any embodiment, the dual filters will adapt a state estimate and a parameter estimate so that a model input-output relationship matches the measured input-output data as closely as possible. This does not guarantee that the model state converges to physical state values. An exemplary embodiment takes extra steps to ensure that one model state converges to SOC. In yet another embodiment, the cell model used for dual filtering may be augmented with a secondary cell model that includes as outputs those states that must converge to their correct values. A specific example of such an augmented cell model (with output augmented with SOC) is g(xk,uk10 -[OCV(Zk)-Rik+hk+Gfkl (8) L Zk J
The augmented model output is compared to a measured output in the dual filter. In an exemplary embodiment, a measured value for SOC may be approximated using 2k derived as Yk OCV(Zk) - Rik OCV(Zk) ^ Yk +Rik (9) Zk = OCV ' (Vk +Rik ).
By measuring the voltage of a cell under load, the cell current, and having knowledge of R, (perhaps through B from a dual filter), and knowing the inverse OCV function for the cell chemistry, this example computes a noisy estimate of SOC, Zk .
In this example, a dual filter is run on this modified model, with the "measured" information in the measurement update being v k Zk Experimentation has shown that while the "noise" of 2k (short-term bias due to hysteresis effects and polarization voltages being ignored) prohibits it from being used as the primary estimator of SOC, its expected long-term behavior in a dynamic environment is accurate, and maintains the accuracy of the SOC state in the dual filter.
Another exemplary embodiment includes methods for estimating important aspects of SOH without employing a full dual filter. The full dual filter method may be computationally intensive. If precise values for the full set of cell model parameters are not necessary, then other methods potentially less complex or computationally intensive might be used. The exemplary methodologies determine cell capacity and resistance using filtering methods. The change in capacity and resistance from the nominal "new-cell" values give capacity fade and power fade, which are the most commonly employed indicators of cell SOH.
In this example, to estimate cell resistance using a filtering mechanism, we formulate a model:
Rk [m +1] = Rk [m] + rk [m] (10) Yk [m] = OCV(Zk [mD -1k [m]Rk [m] + ek [in]
where Rk[m] is the cell resistance and is modeled as a constant value with a fictitious noise process rk allowing adaptation. yk[m] is a crude estimate of the cell's voltage, ik is the cell current, and ekmodels estimation error. If we use an estimate of Zk from the state filter in a dual estimator, or from some other source, then we simply apply a filter to this model to estimate cell resistance. In the standard filter, we compare the model's prediction of yk with the true measured cell voltage, and use the difference to adapt Rk .
Note that the above model may be extended to handle different values of resistance for a variety of conditions of the cell. For example, differences based on charge and discharge, different SOCs, and different temperatures. The scalar Rk would be changed into a vector comprising all of the resistance values being modified, and the appropriate element from the vector would be used each time step of the filter during the calculations.
In this example, to estimate cell capacity using a filter, we again formulate a cell model:
Ck [m + 1] = Ck [m] + r/, [m]
0=zk [m] -z/c [m-1] + 17ii /c [m-1]At/Ck [m -1] +ek [m]. (11) Again, a filter is formulated using this model to produce a capacity estimate. As the filter runs, the computation in the second equation (right-hand-side) is compared to zero, and the difference is used to update the capacity estimate.
Note that good estimates of the present and previous states-of-charge are desired, possibly from a filter estimating SOC.
Estimated capacity may again be a function of temperature (and so forth), if desired, by employing a capacity vector, from which the appropriate element is used in each time step during calculations.
Industrial Applicability Thus, a method for simultaneous estimation of cell state and parameters has been described in conjunction with a number of specific embodiments. One or more embodiments use one or more Kalman filters. Some embodiments use one or more extended Kalman filters. Further, some embodiments include a mechanism to force convergence of state-of-charge. One or more embodiments include a simplified parameter filter to estimate resistance, while some embodiments include a simplified parameter filter to estimate total capacity. The present invention is applicable to a broad range of applications, and cell electrochemistries.
The disclosed method may be embodied in the form of computer-implemented processes and apparatuses for practicing those processes. The method can also be embodied in the form of computer program code containing instructions embodied in tangible media 52, such as floppy diskettes, CD-ROMs, hard drives, or any other computer-readable storage medium, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus capable of executing the method. The present method can also be embodied in the form of computer program code, for example, whether stored in a storage medium, loaded into and/or executed by a computer, or as data signal 54 transmitted whether a modulated carrier wave or not, over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus capable of executing the method. When implemented on a general-purpose microprocessor, the computer program code segments configure the microprocessor to create specific logic circuits.
It will be appreciated that the use of first and second or other similar nomenclature for denoting similar items is not intended to specify or imply any particular order unless otherwise stated. Furthermore, the use of the terminology "a" and "at least one of" shall each be associated with the meaning "one or more" unless specifically stated otherwise.
While the invention has been described with reference to an exemplary embodiment, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential scope thereof. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed as the best mode contemplated for carrying out this invention, but that the invention will include all embodiments falling within the scope of the appended claims.
xk+l =f(xk,uk,ek) +Wk (1) Yk -g(xk,uk, k)+vk, where xk is the system state, dk is the set of time varying model parameters, ukis the exogenous input, ykis the system output, and wk and vkare "noise" inputs-all quantities may be scalars or vectors. f(-,=,=) and g(=,=;) are functions defined by the cell model being used. Non-time-varying numeric values required by the model may be embedded within f(=,,=) and g(=,=,=) , and are not included in k.
The system state includes, at least, a minimum amount of information, together with the present input and a mathematical model of the cell, needed to predict the present output. For a cell 12, the state might include: SOC, polarization voltage levels with respect to different time constants, and hysteresis levels, for example. The system exogenous input uk includes at minimum the present cell current ik , and may, optionally, include cell temperature (unless temperature change is itself modeled in the state) .
The system parameters k are the values that change only slowly with time, in such a way that they may not be directly determined with knowledge of the system measured input and output. These might include: cell capacity, resistance, polarization voltage time constant(s), polarization voltage blending factor(s), hysteresis blending factor(s), hysteresis rate constant(s), efficiency factor(s), and so forth. The model output Yk corresponds to physically measurable cell quantities or those directly computable from measured quantities at minimum for example, the cell voltage under load.
A mathematical model of parameter dynamics is also utilized. An exemplary model has the form:
k+1 = k + r k ( 2 dk - g(xkIukIek) +ek.
The first equation states that the parameters are essentially constant, but that they may change slowly over time, in this instance, modeled by a fictitious "noise"
process denoted, rk . The "output" of the optimum parameter dynamics is the cell output estimate plus some estimation error ek .
With models of the system state dynamics and parameter dynamics defined, in an exemplary embodiment, a procedure of dual filtering is applied. Once again, alternatively, dual Kalman filters may be employed, or dual extended Kalman filters. Furthermore, combinations of the abovementioned may also be employed. Table 1 is an exemplary implementation of the methodology and system utilizing dual extended Kalman filtering. The procedure is initialized by setting the parameter estimate to the best guess of the true parameters, e.g., 6 = E[B0] , and by setting the state estimate x to the best estimate of the cell state, e.g., x=E[x,] . The estimation-error covariance matrices EC and Et are also initialized. For example, an initialization of SOC might be estimated/based on a cell voltage in a look-up table, or information that was previously stored when a battery pack/cell was last powered down. Other examples might incorporate the length of time that the battery system had rested since powerdown and the like.
Table 1: Dual extended Kalman filter for state and weight update.
State-space models:
xk+1 f(xk,uk,Ok)+Wkand k+1 - k+1k vk g(xkIukIek)+ Vk dk _g(xkIukIek)+ek where wk , Vk , rk and ek are independent, zero-mean, Gaussian noise processes of covariance matrices E,, , Z, , E, and E, respectively.
Definitions:
_ af(xk-I,uk-1,lb x_ ag(xk,uk, k) dg(xk,u ,8) Ak-1 ax Ck ax Ck d k-1 xk_t_xk I k Xk_Xk 9=BA.
Initialization: For k=0, set o = E[90 ], EB O = E[(Bo - eo )( o - eo )T
xo = E[xo ], E C,o = E[(xo - $ )(xo - xo )T
Computation: For k =1,2,..., compute :
Time update for the weight filter [I corrected the following equation]
Bk - 0111--1 k 0,k-1 Time update for the state filter xk _flxk-l~uk-11ek) + T
EY k = Ak-1ZC k-1Ak-1 +Ew.
Measurement update for the state filter Lk = EX,k (Ck )T LCk Ex,k (Ck )T +E, xk -xk +Lk[Yk -g(xk,ukIek )]
Ex,k -(I-LkCk)Ex,k' Measurement update for the weight filter Lk - Ek (Ck) [Ck E/''g,k (Ck) +Ee Bk = /Bk + Lk [Yk - g(xk , uk , ek)}
Ek - (I LkCk )Ee In this example, several steps are performed in each measurement interval. First, the previous parameter estimate is propagated forward in time. The new parameter estimate is equal to the old parameter estimate, and the parameter error uncertainty is larger due to the passage of time (accommodated for in the model by the fictitious driving noise rk ). Various possibilities exist for updating the parameter uncertainty estimate-the table gives only one example. The state estimate and its uncertainty are propagated forward one step in time. A measurement of the cell output is made, and compared to the predicted output based on the state estimate, i and parameter estimate, 0 the difference is used to update the values of i and B. It may readily be appreciated that the steps outlined in the table may be performed in a variety of orders. While the table lists an exemplary ordering for the purposes of illustration, those skilled in the art will be able to identify many equivalent ordered sets of equations.
Turning now to FIG. 2 as well, an exemplary implementation of an exemplary embodiment of the invention is depicted. Two filters run in parallel are depicted. One filter 101/102 adapting the state estimate, x and one filter 103/104 adapting the parameter estimate, B . Additionally, information exchange between the filters permits updates from the parameter estimation filter 103/104 to the state estimation filter 101/102. Both filters 101/102, 103/104 have a time update or prediction aspect and a measurement update or correction aspect. State time update /prediction block 101 receives as input the previous exogenous input uk_, (which might include cell current and/or temperature, for example) along with the previously estimated system state value xk, and state uncertainty estimate Eck_1, and present predicted parameters Bk . The state time update/prediction block 101 provides predicted state 2k and predicted state uncertainty irk output to state measurement update/correction block 102, and to parameter measurement update/correction block 104. State measurement update/correction block 102 provides current system state estimate 241 and state uncertainty estimate Ck Parameter time update/prediction block 103 receives as input the previous time varying parameters estimate O and parameter uncertainty estimate Z+k_i . Parameter time update/prediction block 103 outputs predicted parameters k and predicted parameter uncertainty Z- to the parameter measurement update/correction block 104 and state time update block 101. Parameter measurement update block 104, which provides current parameter estimate ek and parameter uncertainty estimate Y-e<<. It will also be appreciated that a minus notation denotes that the vector is the result of the prediction components 101, 103 of the filters, while the plus notation denotes that the vector is the result of the correction component 102/104 of the filters.
Embodiments of this invention require a mathematical model of cell state and output dynamics for the particular application. In the examples above, this is accomplished by defining specific functions for f(;;) and g(=,=;) . An exemplary embodiment uses a cell model that includes effects due to one or more of the open-circuit-voltage (OCV) for the cell, internal resistance, voltage polarization time constants, and a hysteresis level. For the purpose of example, parameter values are fitted to this model structure to model the dynamics of high-power Lithium-Ion Polymer Battery (LiPB) cells, although the structure and methods presented here are general and apply to other electrochemistries.
In this example, SOC is captured by one state of the model. This equation is zk[m+1]=zk[m]-(r1At/Ck)Zk[m] (3) where At represents the inter-sample period (in seconds), Ck represents the capacity of cell number k in the pack (in ampere-seconds) , zk[m] is the SOC of cell k at time index m , ik , is the current out of cell k , and 77, is the Coulombic efficiency of a cell at current level ik.
In this example, the polarization voltage levels are captured by several filter states. If we let there be nf polarization voltage time constants, then Jk[in+1] = Affk[m]+Bfik[m]. (4) The matrix Af eRi."f'"f may be a diagonal matrix with real-valued polarization voltage time constants a, ===a"f . If so, the system is stable if all entries have magnitude less than one. The vector B f E R"f=X' may simply be set to n f "1"s . The entries of Bf are not critical as long as they are non-zero. The value of nfentries in the Afmatrix are chosen as part of the system identification procedure to best fit the model parameters to measured cell data. The Af and Bf matrices may vary with time and other factors pertinent to the present battery pack operating condition.
In this example, the hysteresis level is captured by a single state hk[m+1]=exp - 77;Zk[C]YAt) hk[rn]+ 1 -exp - 77rik[m]rAt sg(ik[m]), (5) where y is' the hysteresis rate constant, again found by system identification.
In this example, the overall model state is xk[m]= [fk[m]T hk[nz] zk[m]]T, where other orderings of states are possible. The state equation for the model is formed by combining all of the individual equations identified above.
In this example, the output equation that combines the state values to predict cell voltage is Vk [m] = OCV(zk 1M])+ G.fk [m] - Rik [m] + Mhk [m], (6) where Ge1i'xn, is a vector of polarization voltage blending factors g,===g,,, that blend the polarization voltage states together in the output, R is the cell resistance (different values may be used for discharge/charge), and M is the hysteresis blending factor. Note, G may be constrained such that the dc-gain from ikto Gfk is zero.
In this example, the parameters are T
0=[z7, C, a,...anj, 91...gnri, y, R, MI = (7) In any embodiment, the dual filters will adapt a state estimate and a parameter estimate so that a model input-output relationship matches the measured input-output data as closely as possible. This does not guarantee that the model state converges to physical state values. An exemplary embodiment takes extra steps to ensure that one model state converges to SOC. In yet another embodiment, the cell model used for dual filtering may be augmented with a secondary cell model that includes as outputs those states that must converge to their correct values. A specific example of such an augmented cell model (with output augmented with SOC) is g(xk,uk10 -[OCV(Zk)-Rik+hk+Gfkl (8) L Zk J
The augmented model output is compared to a measured output in the dual filter. In an exemplary embodiment, a measured value for SOC may be approximated using 2k derived as Yk OCV(Zk) - Rik OCV(Zk) ^ Yk +Rik (9) Zk = OCV ' (Vk +Rik ).
By measuring the voltage of a cell under load, the cell current, and having knowledge of R, (perhaps through B from a dual filter), and knowing the inverse OCV function for the cell chemistry, this example computes a noisy estimate of SOC, Zk .
In this example, a dual filter is run on this modified model, with the "measured" information in the measurement update being v k Zk Experimentation has shown that while the "noise" of 2k (short-term bias due to hysteresis effects and polarization voltages being ignored) prohibits it from being used as the primary estimator of SOC, its expected long-term behavior in a dynamic environment is accurate, and maintains the accuracy of the SOC state in the dual filter.
Another exemplary embodiment includes methods for estimating important aspects of SOH without employing a full dual filter. The full dual filter method may be computationally intensive. If precise values for the full set of cell model parameters are not necessary, then other methods potentially less complex or computationally intensive might be used. The exemplary methodologies determine cell capacity and resistance using filtering methods. The change in capacity and resistance from the nominal "new-cell" values give capacity fade and power fade, which are the most commonly employed indicators of cell SOH.
In this example, to estimate cell resistance using a filtering mechanism, we formulate a model:
Rk [m +1] = Rk [m] + rk [m] (10) Yk [m] = OCV(Zk [mD -1k [m]Rk [m] + ek [in]
where Rk[m] is the cell resistance and is modeled as a constant value with a fictitious noise process rk allowing adaptation. yk[m] is a crude estimate of the cell's voltage, ik is the cell current, and ekmodels estimation error. If we use an estimate of Zk from the state filter in a dual estimator, or from some other source, then we simply apply a filter to this model to estimate cell resistance. In the standard filter, we compare the model's prediction of yk with the true measured cell voltage, and use the difference to adapt Rk .
Note that the above model may be extended to handle different values of resistance for a variety of conditions of the cell. For example, differences based on charge and discharge, different SOCs, and different temperatures. The scalar Rk would be changed into a vector comprising all of the resistance values being modified, and the appropriate element from the vector would be used each time step of the filter during the calculations.
In this example, to estimate cell capacity using a filter, we again formulate a cell model:
Ck [m + 1] = Ck [m] + r/, [m]
0=zk [m] -z/c [m-1] + 17ii /c [m-1]At/Ck [m -1] +ek [m]. (11) Again, a filter is formulated using this model to produce a capacity estimate. As the filter runs, the computation in the second equation (right-hand-side) is compared to zero, and the difference is used to update the capacity estimate.
Note that good estimates of the present and previous states-of-charge are desired, possibly from a filter estimating SOC.
Estimated capacity may again be a function of temperature (and so forth), if desired, by employing a capacity vector, from which the appropriate element is used in each time step during calculations.
Industrial Applicability Thus, a method for simultaneous estimation of cell state and parameters has been described in conjunction with a number of specific embodiments. One or more embodiments use one or more Kalman filters. Some embodiments use one or more extended Kalman filters. Further, some embodiments include a mechanism to force convergence of state-of-charge. One or more embodiments include a simplified parameter filter to estimate resistance, while some embodiments include a simplified parameter filter to estimate total capacity. The present invention is applicable to a broad range of applications, and cell electrochemistries.
The disclosed method may be embodied in the form of computer-implemented processes and apparatuses for practicing those processes. The method can also be embodied in the form of computer program code containing instructions embodied in tangible media 52, such as floppy diskettes, CD-ROMs, hard drives, or any other computer-readable storage medium, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus capable of executing the method. The present method can also be embodied in the form of computer program code, for example, whether stored in a storage medium, loaded into and/or executed by a computer, or as data signal 54 transmitted whether a modulated carrier wave or not, over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus capable of executing the method. When implemented on a general-purpose microprocessor, the computer program code segments configure the microprocessor to create specific logic circuits.
It will be appreciated that the use of first and second or other similar nomenclature for denoting similar items is not intended to specify or imply any particular order unless otherwise stated. Furthermore, the use of the terminology "a" and "at least one of" shall each be associated with the meaning "one or more" unless specifically stated otherwise.
While the invention has been described with reference to an exemplary embodiment, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential scope thereof. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed as the best mode contemplated for carrying out this invention, but that the invention will include all embodiments falling within the scope of the appended claims.
Claims (9)
1. A method for estimating present states and present parameters of an electrochemical cell system comprising:
estimating a state value of the electrochemical cell with a cell state filter to estimate said state value, said estimating the state value comprising:
making an internal state prediction of said cell comprising:
determining a parameter estimate;
determining a current measurement;
determining a voltage measurement;
using said parameter estimate, said current measurement and said voltage measurement in a parameter mathematical model to make said internal parameter prediction;
making a state uncertainty prediction of said internal state prediction comprising:
using said parameter estimate and said current measurement, said voltage measurement, and said temperature measurement in a state mathematical model to make said uncertainty prediction;
correcting said internal state prediction and said state uncertainty prediction; and applying a state algorithm that iterates said making an internal state prediction, said making a state uncertainty prediction and said correcting to yield an ongoing estimation to said internal state prediction and an ongoing uncertainty to said internal state prediction;
and estimating a parameter value of the electrochemical cell with a cell parameter filter to estimate said parameter value, said estimating the parameter value comprising:
making an internal state prediction of said cell;
making a parameter uncertainty prediction of said internal state prediction;
correcting said internal parameter prediction and said parameter uncertainty prediction; and applying a parameter algorithm that iterates said making an internal parameter prediction, said making a parameter uncertainty prediction and said correcting to yield an ongoing estimation to said internal parameter prediction and an ongoing uncertainty to said internal parameter prediction, wherein said internal parameter includes one or more of:
a polarization voltage time constant, a polarization voltage blending factor, a hysteresis blending factor, and a hysteresis rate constant.
estimating a state value of the electrochemical cell with a cell state filter to estimate said state value, said estimating the state value comprising:
making an internal state prediction of said cell comprising:
determining a parameter estimate;
determining a current measurement;
determining a voltage measurement;
using said parameter estimate, said current measurement and said voltage measurement in a parameter mathematical model to make said internal parameter prediction;
making a state uncertainty prediction of said internal state prediction comprising:
using said parameter estimate and said current measurement, said voltage measurement, and said temperature measurement in a state mathematical model to make said uncertainty prediction;
correcting said internal state prediction and said state uncertainty prediction; and applying a state algorithm that iterates said making an internal state prediction, said making a state uncertainty prediction and said correcting to yield an ongoing estimation to said internal state prediction and an ongoing uncertainty to said internal state prediction;
and estimating a parameter value of the electrochemical cell with a cell parameter filter to estimate said parameter value, said estimating the parameter value comprising:
making an internal state prediction of said cell;
making a parameter uncertainty prediction of said internal state prediction;
correcting said internal parameter prediction and said parameter uncertainty prediction; and applying a parameter algorithm that iterates said making an internal parameter prediction, said making a parameter uncertainty prediction and said correcting to yield an ongoing estimation to said internal parameter prediction and an ongoing uncertainty to said internal parameter prediction, wherein said internal parameter includes one or more of:
a polarization voltage time constant, a polarization voltage blending factor, a hysteresis blending factor, and a hysteresis rate constant.
2. The method of claim 1, wherein said correcting said internal state prediction comprises:
computing a gain factor;
computing a corrected internal state prediction using said gain factor, said voltage measurement and said internal state prediction; and computing a corrected state uncertainty prediction using said gain factor and said state uncertainty prediction.
computing a gain factor;
computing a corrected internal state prediction using said gain factor, said voltage measurement and said internal state prediction; and computing a corrected state uncertainty prediction using said gain factor and said state uncertainty prediction.
3. The method of claim 2, wherein said applying a state algorithm comprises using said corrected internal state prediction and said corrected state uncertainty prediction to obtain predictions for a next time step where said state algorithm repeats.
4. The method of claim 3, wherein said state algorithm is at least one of a Kalman Filter and an extended Kalman Filter.
5. The method of claim 4, wherein said internal state includes at least one of state-of-charge, voltage polarization levels, and hysteresis levels.
6. An apparatus configured to estimate present states and present parameters of an electrochemical cell comprising:
a first component configured to estimate a cell state value, said first component comprising:
a component configured to make an internal state prediction of said cell, said component comprising:
a component configured to determine a parameter estimate;
a component configured to determine a current measurement;
a component configured to determine a voltage measurement;
a component configured to use said parameter estimate and said current measurement and said voltage measurement in a mathematical model to make said internal state prediction;
a component configured to make a state uncertainty prediction of said internal state prediction;
a component configured to correct said internal state prediction and said state uncertainty prediction;
a component configured to apply a state algorithm that iterates steps taken by said component configured to make an internal state prediction, said component configured to make a state uncertainty prediction and said component configured to correct to yield an ongoing estimation to said internal state prediction and an ongoing uncertainty to said internal state prediction;
and a second component configured to estimate a cell parameter value, said second component comprising:
a component configured to make an internal parameter prediction of said cell;
a component configured to make a parameter uncertainty prediction of said internal parameter prediction;
a component configured to correct said internal parameter prediction and said parameter uncertainty prediction; and a component configured to apply a parameter algorithm that iterates steps taken by said component configured to make an internal parameter prediction, said component configured to make a parameter uncertainty prediction and said component configured to correct to yield an ongoing estimation to said internal parameter prediction and an ongoing uncertainty to said internal parameter prediction, wherein said internal parameter includes one or more of:
a polarization voltage time constant, a polarization voltage blending factor, a hysteresis blending factor, and a hysteresis rate constant.
a first component configured to estimate a cell state value, said first component comprising:
a component configured to make an internal state prediction of said cell, said component comprising:
a component configured to determine a parameter estimate;
a component configured to determine a current measurement;
a component configured to determine a voltage measurement;
a component configured to use said parameter estimate and said current measurement and said voltage measurement in a mathematical model to make said internal state prediction;
a component configured to make a state uncertainty prediction of said internal state prediction;
a component configured to correct said internal state prediction and said state uncertainty prediction;
a component configured to apply a state algorithm that iterates steps taken by said component configured to make an internal state prediction, said component configured to make a state uncertainty prediction and said component configured to correct to yield an ongoing estimation to said internal state prediction and an ongoing uncertainty to said internal state prediction;
and a second component configured to estimate a cell parameter value, said second component comprising:
a component configured to make an internal parameter prediction of said cell;
a component configured to make a parameter uncertainty prediction of said internal parameter prediction;
a component configured to correct said internal parameter prediction and said parameter uncertainty prediction; and a component configured to apply a parameter algorithm that iterates steps taken by said component configured to make an internal parameter prediction, said component configured to make a parameter uncertainty prediction and said component configured to correct to yield an ongoing estimation to said internal parameter prediction and an ongoing uncertainty to said internal parameter prediction, wherein said internal parameter includes one or more of:
a polarization voltage time constant, a polarization voltage blending factor, a hysteresis blending factor, and a hysteresis rate constant.
7. The apparatus of claim 6, wherein said component configured to make a state uncertainty prediction comprises a component configured to use said parameter estimate, a current measurement and a voltage measurement in a mathematical model to make said state uncertainty prediction.
8. The apparatus of claim 7, wherein said component configured to correct said internal state prediction comprises:
a component configured to compute a gain factor;
a component configured to compute a corrected internal state prediction using said gain factor, said voltage measurement and said internal state prediction;
and a component configured to compute a corrected state uncertainty prediction using said gain factor and said state uncertainty prediction.
a component configured to compute a gain factor;
a component configured to compute a corrected internal state prediction using said gain factor, said voltage measurement and said internal state prediction;
and a component configured to compute a corrected state uncertainty prediction using said gain factor and said state uncertainty prediction.
9. The apparatus of claim 8, wherein said component configured to apply a state algorithm comprises a component configured to use said corrected internal state prediction and said corrected state uncertainty prediction to obtain predictions for a next time step where said state algorithm repeats.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US10/985,617 | 2004-11-11 | ||
| US10/985,617 US8103485B2 (en) | 2004-11-11 | 2004-11-11 | State and parameter estimation for an electrochemical cell |
| PCT/KR2004/003103 WO2006052043A1 (en) | 2004-11-11 | 2004-11-29 | State and parameter estimation for an electrochemical cell |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2585921A1 CA2585921A1 (en) | 2006-05-18 |
| CA2585921C true CA2585921C (en) | 2012-03-20 |
Family
ID=36317422
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2585921A Expired - Lifetime CA2585921C (en) | 2004-11-11 | 2004-11-29 | State and parameter estimation for an electrochemical cell |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US8103485B2 (en) |
| JP (1) | JP4511600B2 (en) |
| KR (1) | KR100818520B1 (en) |
| BR (1) | BRPI0419118B1 (en) |
| CA (1) | CA2585921C (en) |
| RU (1) | RU2361333C2 (en) |
| TW (1) | TWI268637B (en) |
| WO (1) | WO2006052043A1 (en) |
Families Citing this family (113)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7321220B2 (en) * | 2003-11-20 | 2008-01-22 | Lg Chem, Ltd. | Method for calculating power capability of battery packs using advanced cell model predictive techniques |
| US7593821B2 (en) * | 2004-11-23 | 2009-09-22 | Lg Chem, Ltd. | Method and system for joint battery state and parameter estimation |
| US8446127B2 (en) * | 2005-08-03 | 2013-05-21 | California Institute Of Technology | Methods for thermodynamic evaluation of battery state of health |
| US7723957B2 (en) * | 2005-11-30 | 2010-05-25 | Lg Chem, Ltd. | System, method, and article of manufacture for determining an estimated battery parameter vector |
| US7400115B2 (en) | 2006-02-09 | 2008-07-15 | Lg Chem, Ltd. | System, method, and article of manufacture for determining an estimated combined battery state-parameter vector |
| KR100804698B1 (en) * | 2006-06-26 | 2008-02-18 | 삼성에스디아이 주식회사 | The method of assuming the state of charge of the battery, battery management system using the method and the driving method of the battery management system using the method |
| EP2160616A1 (en) * | 2007-07-02 | 2010-03-10 | ABB Research LTD | State of charge determination |
| US7830119B2 (en) | 2007-08-29 | 2010-11-09 | Gm Global Technology Operations, Inc. | Adaptive battery estimator and method |
| US8628872B2 (en) * | 2008-01-18 | 2014-01-14 | Lg Chem, Ltd. | Battery cell assembly and method for assembling the battery cell assembly |
| US7994755B2 (en) | 2008-01-30 | 2011-08-09 | Lg Chem, Ltd. | System, method, and article of manufacture for determining an estimated battery cell module state |
| US8067111B2 (en) * | 2008-06-30 | 2011-11-29 | Lg Chem, Ltd. | Battery module having battery cell assembly with heat exchanger |
| US8426050B2 (en) * | 2008-06-30 | 2013-04-23 | Lg Chem, Ltd. | Battery module having cooling manifold and method for cooling battery module |
| US7883793B2 (en) * | 2008-06-30 | 2011-02-08 | Lg Chem, Ltd. | Battery module having battery cell assemblies with alignment-coupling features |
| US9140501B2 (en) * | 2008-06-30 | 2015-09-22 | Lg Chem, Ltd. | Battery module having a rubber cooling manifold |
| US9759495B2 (en) * | 2008-06-30 | 2017-09-12 | Lg Chem, Ltd. | Battery cell assembly having heat exchanger with serpentine flow path |
| US8202645B2 (en) | 2008-10-06 | 2012-06-19 | Lg Chem, Ltd. | Battery cell assembly and method for assembling the battery cell assembly |
| US8116998B2 (en) | 2009-01-30 | 2012-02-14 | Bae Systems Controls, Inc. | Battery health assessment estimator |
| US9337456B2 (en) * | 2009-04-20 | 2016-05-10 | Lg Chem, Ltd. | Frame member, frame assembly and battery cell assembly made therefrom and methods of making the same |
| US8852778B2 (en) * | 2009-04-30 | 2014-10-07 | Lg Chem, Ltd. | Battery systems, battery modules, and method for cooling a battery module |
| US8403030B2 (en) * | 2009-04-30 | 2013-03-26 | Lg Chem, Ltd. | Cooling manifold |
| US8663828B2 (en) * | 2009-04-30 | 2014-03-04 | Lg Chem, Ltd. | Battery systems, battery module, and method for cooling the battery module |
| US8663829B2 (en) | 2009-04-30 | 2014-03-04 | Lg Chem, Ltd. | Battery systems, battery modules, and method for cooling a battery module |
| US8399118B2 (en) * | 2009-07-29 | 2013-03-19 | Lg Chem, Ltd. | Battery module and method for cooling the battery module |
| US8703318B2 (en) * | 2009-07-29 | 2014-04-22 | Lg Chem, Ltd. | Battery module and method for cooling the battery module |
| US8399119B2 (en) * | 2009-08-28 | 2013-03-19 | Lg Chem, Ltd. | Battery module and method for cooling the battery module |
| US8918299B2 (en) * | 2009-12-02 | 2014-12-23 | American Electric Vehicles, Inc. | System and method for maximizing a battery pack total energy metric |
| US8427105B2 (en) * | 2009-12-02 | 2013-04-23 | Gregory L. Plett | System and method for equalizing a battery pack during a battery pack charging process |
| US8041522B2 (en) * | 2009-12-02 | 2011-10-18 | American Electric Vehicles, Ind. | System and method for recursively estimating battery cell total capacity |
| TWI411796B (en) * | 2009-12-22 | 2013-10-11 | Ind Tech Res Inst | Apparatus for estimating battery's state of health |
| RU2491566C1 (en) * | 2010-02-18 | 2013-08-27 | Ниссан Мотор Ко., Лтд. | Device for estimation of battery state and method for estimation of battery state |
| JP5691592B2 (en) * | 2010-02-18 | 2015-04-01 | 日産自動車株式会社 | Battery state estimation device |
| US8341449B2 (en) | 2010-04-16 | 2012-12-25 | Lg Chem, Ltd. | Battery management system and method for transferring data within the battery management system |
| US9147916B2 (en) | 2010-04-17 | 2015-09-29 | Lg Chem, Ltd. | Battery cell assemblies |
| US8936394B2 (en) | 2010-05-25 | 2015-01-20 | GM Global Technology Operations LLC | Real-time estimation of cell core temperature during period of rest |
| US8942935B2 (en) * | 2010-06-14 | 2015-01-27 | Medtronic, Inc. | Charge level measurement |
| US9726732B2 (en) | 2010-06-22 | 2017-08-08 | GM Global Technology Operations LLC | Adaptive battery parameter extraction and SOC estimation for lithium-ion battery |
| AT508307B1 (en) * | 2010-06-24 | 2015-01-15 | Avl List Gmbh | METHOD FOR DETERMINING A STATE OF CELLS OF AN ENERGY STORAGE DEVICE |
| JP5842421B2 (en) * | 2010-07-20 | 2016-01-13 | 日産自動車株式会社 | Battery state estimation device |
| DE102010038646A1 (en) * | 2010-07-29 | 2012-02-02 | Robert Bosch Gmbh | Method and device for estimating the performance of at least one battery unit of a rechargeable battery |
| US8469404B2 (en) | 2010-08-23 | 2013-06-25 | Lg Chem, Ltd. | Connecting assembly |
| US8353315B2 (en) | 2010-08-23 | 2013-01-15 | Lg Chem, Ltd. | End cap |
| US8920956B2 (en) | 2010-08-23 | 2014-12-30 | Lg Chem, Ltd. | Battery system and manifold assembly having a manifold member and a connecting fitting |
| US9065292B2 (en) | 2010-08-23 | 2015-06-23 | California Institute Of Technology | Methods and systems for charging electrochemical cells |
| US8758922B2 (en) | 2010-08-23 | 2014-06-24 | Lg Chem, Ltd. | Battery system and manifold assembly with two manifold members removably coupled together |
| US9005799B2 (en) | 2010-08-25 | 2015-04-14 | Lg Chem, Ltd. | Battery module and methods for bonding cell terminals of battery cells together |
| US8662153B2 (en) | 2010-10-04 | 2014-03-04 | Lg Chem, Ltd. | Battery cell assembly, heat exchanger, and method for manufacturing the heat exchanger |
| TW201224485A (en) | 2010-12-02 | 2012-06-16 | Ind Tech Res Inst | State-of-charge estimation method and battery control unit |
| US20120179435A1 (en) * | 2011-01-10 | 2012-07-12 | Ford Global Technologies, Llc | Method For Determining A Power Capability For A Battery |
| US8698348B2 (en) | 2011-01-10 | 2014-04-15 | Ford Global Technologies, Llc | System and method for managing a power source in a vehicle |
| US8288031B1 (en) | 2011-03-28 | 2012-10-16 | Lg Chem, Ltd. | Battery disconnect unit and method of assembling the battery disconnect unit |
| US8449998B2 (en) | 2011-04-25 | 2013-05-28 | Lg Chem, Ltd. | Battery system and method for increasing an operational life of a battery cell |
| US9178192B2 (en) | 2011-05-13 | 2015-11-03 | Lg Chem, Ltd. | Battery module and method for manufacturing the battery module |
| CN102298118A (en) * | 2011-05-17 | 2011-12-28 | 杭州电子科技大学 | On-line synchronous estimating method for model parameters and remaining power of battery |
| US10234512B2 (en) * | 2011-06-11 | 2019-03-19 | Sendyne Corporation | Current-based cell modeling |
| US8974929B2 (en) | 2011-06-30 | 2015-03-10 | Lg Chem, Ltd. | Heating system for a battery module and method of heating the battery module |
| US8974928B2 (en) | 2011-06-30 | 2015-03-10 | Lg Chem, Ltd. | Heating system for a battery module and method of heating the battery module |
| US8993136B2 (en) | 2011-06-30 | 2015-03-31 | Lg Chem, Ltd. | Heating system for a battery module and method of heating the battery module |
| US8859119B2 (en) | 2011-06-30 | 2014-10-14 | Lg Chem, Ltd. | Heating system for a battery module and method of heating the battery module |
| US9496544B2 (en) | 2011-07-28 | 2016-11-15 | Lg Chem. Ltd. | Battery modules having interconnect members with vibration dampening portions |
| DE102011086803A1 (en) | 2011-11-22 | 2013-05-23 | Ford Global Technologies, Llc | Repair method of a cylinder surface by means of plasma spraying |
| DE102013200912B4 (en) | 2012-02-02 | 2018-05-30 | Ford Global Technologies, Llc | crankcase |
| US10556510B2 (en) | 2012-04-27 | 2020-02-11 | California Institute Of Technology | Accurate assessment of the state of charge of electrochemical cells |
| CN104471415B (en) | 2012-04-27 | 2017-09-01 | 加州理工学院 | embedded chip for battery applications |
| US9511467B2 (en) | 2013-06-10 | 2016-12-06 | Ford Global Technologies, Llc | Cylindrical surface profile cutting tool and process |
| US9079213B2 (en) | 2012-06-29 | 2015-07-14 | Ford Global Technologies, Llc | Method of determining coating uniformity of a coated surface |
| JP5733275B2 (en) * | 2012-07-13 | 2015-06-10 | トヨタ自動車株式会社 | Battery pack control device and battery pack reuse determination method |
| US8935043B2 (en) * | 2013-01-29 | 2015-01-13 | Ford Global Technologies, Llc | Temperature compensated battery parameter estimation |
| US10664562B2 (en) * | 2013-02-24 | 2020-05-26 | Fairchild Semiconductor Corporation and University of Connecticut | Battery state of charge tracking, equivalent circuit selection and benchmarking |
| CN104007390B (en) * | 2013-02-24 | 2018-03-20 | 快捷半导体(苏州)有限公司 | Battery state of charge tracking, equivalent circuit selection and reference test method and system |
| FR3009389B1 (en) * | 2013-08-02 | 2017-02-10 | Commissariat Energie Atomique | ENERGY MANAGEMENT IN A BATTERY |
| US9368841B2 (en) | 2013-08-30 | 2016-06-14 | Ford Global Technologies, Llc | Battery power capability estimation at vehicle start |
| KR20150029204A (en) * | 2013-09-09 | 2015-03-18 | 삼성에스디아이 주식회사 | Battery pack, apparatus including battery pack, and method of managing battery pack |
| FR3011084A1 (en) | 2013-09-25 | 2015-03-27 | St Microelectronics Grenoble 2 | METHOD FOR DETERMINING THE CHARGING STATE OF A BATTERY OF AN ELECTRONIC DEVICE |
| US9631595B2 (en) | 2013-09-26 | 2017-04-25 | Ford Global Technologies, Llc | Methods and systems for selective engine starting |
| JP2015078918A (en) * | 2013-10-17 | 2015-04-23 | ソニー株式会社 | Open-circuit voltage estimation device, power storage device and method of estimating open-circuit voltage |
| JP2015155859A (en) * | 2014-02-21 | 2015-08-27 | ソニー株式会社 | Battery remaining amount estimating device, battery pack, power storage device, electric vehicle, and battery remaining amount estimating method |
| US9382868B2 (en) | 2014-04-14 | 2016-07-05 | Ford Global Technologies, Llc | Cylinder bore surface profile and process |
| US9581988B2 (en) | 2014-06-05 | 2017-02-28 | Ford Global Technologies, Llc | Method and system for battery state of charge estimation |
| US10288691B2 (en) | 2014-06-05 | 2019-05-14 | Ford Global Technologies, Llc | Method and system for estimating battery model parameters to update battery models used for controls |
| EP2963434B1 (en) * | 2014-06-30 | 2021-08-11 | Foundation Of Soongsil University-Industry Cooperation | Battery state estimation method and system using dual extended kalman filter, and recording medium for performing the method |
| FR3025889B1 (en) * | 2014-09-12 | 2016-11-18 | Commissariat Energie Atomique | MANAGING THE RECHARGE OF THE BATTERY OF AN ELECTRIC VEHICLE |
| FR3029297B1 (en) * | 2014-11-28 | 2016-12-30 | Renault Sa | AUTOMATIC METHOD OF ESTIMATING THE CHARGING STATE OF A CELL OF A BATTERY |
| FR3029299B1 (en) * | 2014-11-28 | 2016-12-09 | Renault Sa | AUTOMATIC METHOD FOR DETERMINING THE CHARGING STATE OF A BATTERY |
| FR3029298B1 (en) * | 2014-11-28 | 2016-12-30 | Renault Sa | AUTOMATIC METHOD OF ESTIMATING THE CHARGING STATE OF A CELL OF A BATTERY |
| FR3029315B1 (en) * | 2014-11-28 | 2016-12-09 | Renault Sa | AUTOMATIC METHOD OF ESTIMATING THE CAPACITY OF A CELL OF A BATTERY |
| CN107533108B (en) * | 2015-03-27 | 2020-03-06 | 松下知识产权经营株式会社 | State estimation device and state estimation method for secondary battery |
| US10220453B2 (en) | 2015-10-30 | 2019-03-05 | Ford Motor Company | Milling tool with insert compensation |
| US10224579B2 (en) | 2015-12-31 | 2019-03-05 | Robert Bosch Gmbh | Evaluating capacity fade in dual insertion batteries using potential and temperature measurements |
| US10263447B2 (en) * | 2016-01-29 | 2019-04-16 | Robert Bosch Gmbh | Secondary battery management system |
| US10243385B2 (en) | 2016-01-29 | 2019-03-26 | Robert Bosch Gmbh | Secondary battery management system |
| US10686321B2 (en) | 2016-01-29 | 2020-06-16 | Robert Bosch Gmbh | Secondary battery management |
| US9960625B2 (en) | 2016-03-31 | 2018-05-01 | Robert Bosch Gmbh | Battery management system with multiple observers |
| US10447046B2 (en) | 2016-09-22 | 2019-10-15 | Robert Bosch Gmbh | Secondary battery management system with remote parameter estimation |
| CN106855612B (en) * | 2017-02-21 | 2019-09-24 | 山东大学 | The fractional order KiBaM battery model and parameter identification method of meter and non-linear capacity characteristic |
| JP6939057B2 (en) * | 2017-04-27 | 2021-09-22 | トヨタ自動車株式会社 | In-vehicle battery system and battery aging estimation method |
| RU197634U1 (en) * | 2017-10-30 | 2020-05-18 | Общество с ограниченной ответственностью "Бэттери Сервис Групп" | DISCHARGE-DIAGNOSTIC BATTERY DEVICE |
| CN108196195B (en) * | 2017-12-27 | 2020-12-22 | 银隆新能源股份有限公司 | Battery assembly method, device and equipment |
| CN108398642B (en) * | 2018-01-10 | 2020-06-02 | 中山大学 | An online calibration method for SOC of lithium-ion power battery |
| CN110687468B (en) * | 2018-06-19 | 2021-01-15 | 华为技术有限公司 | Method and device for estimating state of charge of battery |
| AT521643B1 (en) * | 2018-08-31 | 2020-09-15 | Avl List Gmbh | Method and battery management system for determining a state of health of a secondary battery |
| CN109061509A (en) * | 2018-09-12 | 2018-12-21 | 芜湖楚睿智能科技有限公司 | A kind of battery capacity remaining value evaluation method |
| US12103427B2 (en) * | 2019-10-17 | 2024-10-01 | Samsung Sdi Co., Ltd. | Method for determining the support of energy content and power of a battery |
| KR102680305B1 (en) | 2019-10-30 | 2024-06-28 | 주식회사 엘지에너지솔루션 | Apparatus for Controlling Power of Parallel Multi Battery Pack and Method thereof |
| KR102682497B1 (en) * | 2019-11-26 | 2024-07-08 | 현대자동차주식회사 | Vehicle and controlling method thereof |
| WO2022019664A1 (en) | 2020-07-21 | 2022-01-27 | 주식회사 엘지에너지솔루션 | Device and method for controlling output of parallel multi-pack module |
| CN111812538A (en) * | 2020-07-22 | 2020-10-23 | 兰州兰石恩力微电网有限公司 | Power battery evaluation system |
| KR102471890B1 (en) * | 2020-08-03 | 2022-11-29 | 삼성에스디아이 주식회사 | Method for simulating battery pack |
| CN112462282B (en) * | 2020-11-09 | 2022-03-18 | 西南大学 | Method for determining real-time state of charge of battery pack based on mechanism model |
| KR102790688B1 (en) * | 2021-11-08 | 2025-04-08 | 주식회사 플렉싱크 | A Battery SOH Forecastion Method and Battery Managing Method Using the Same |
| CN115508716A (en) * | 2022-10-27 | 2022-12-23 | 国网河南省电力公司电力科学研究院 | Method for estimating state of storage battery based on discharge voltage drop and alternating current impedance |
| US12248348B2 (en) * | 2023-06-28 | 2025-03-11 | TriSpace Technologies (OPC) Pvt. Ltd. | System and method for optimizing bill of material cost and power performance of platform system-on-chip for battery management system |
| US20250183387A1 (en) * | 2023-12-04 | 2025-06-05 | National Chung Shan Institute Of Science And Technology | Automatic pressure adjustment device for measuring pressurization demand of battery cells |
| US12436200B2 (en) * | 2024-02-21 | 2025-10-07 | Deere & Company | Battery load test in agricultural utility vehicle |
Family Cites Families (34)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3026690B2 (en) * | 1992-11-30 | 2000-03-27 | 株式会社リコー | Potential estimation device |
| US5825155A (en) * | 1993-08-09 | 1998-10-20 | Kabushiki Kaisha Toshiba | Battery set structure and charge/ discharge control apparatus for lithium-ion battery |
| US6064180A (en) | 1996-10-29 | 2000-05-16 | General Motors Corporation | Method and apparatus for determining battery state-of-charge using neural network architecture |
| US6018227A (en) | 1998-06-22 | 2000-01-25 | Stryker Corporation | Battery charger especially useful with sterilizable, rechargeable battery packs |
| EP1160953B1 (en) * | 2000-05-29 | 2009-12-02 | Panasonic Corporation | Method for charging battery |
| DE10056969A1 (en) * | 2000-11-17 | 2002-05-23 | Bosch Gmbh Robert | Determining battery charge involves computing charge in first range of operation on basis of model computation in which measured and computed battery voltages are equalized by feedback |
| US6359419B1 (en) * | 2000-12-27 | 2002-03-19 | General Motors Corporation | Quasi-adaptive method for determining a battery's state of charge |
| DE10106505A1 (en) * | 2001-02-13 | 2002-08-29 | Bosch Gmbh Robert | Method and device for condition detection of technical systems such as energy storage |
| DE10106508A1 (en) * | 2001-02-13 | 2002-08-29 | Bosch Gmbh Robert | Method and arrangement for determining the performance of a battery |
| JP4292721B2 (en) * | 2001-02-14 | 2009-07-08 | 株式会社日本自動車部品総合研究所 | Battery state control method for hybrid vehicle |
| US6441586B1 (en) * | 2001-03-23 | 2002-08-27 | General Motors Corporation | State of charge prediction method and apparatus for a battery |
| JP4193371B2 (en) | 2001-04-25 | 2008-12-10 | トヨタ自動車株式会社 | Battery capacity control device |
| US6618681B2 (en) * | 2001-05-02 | 2003-09-09 | Honeywell International Inc. | Method and apparatus for predicting the available energy of a battery |
| JP4523738B2 (en) * | 2001-06-07 | 2010-08-11 | パナソニック株式会社 | Secondary battery remaining capacity control method and apparatus |
| JP3934365B2 (en) | 2001-06-20 | 2007-06-20 | 松下電器産業株式会社 | Battery charge / discharge control method |
| US6876175B2 (en) | 2001-06-29 | 2005-04-05 | Robert Bosch Gmbh | Methods for determining the charge state and/or the power capacity of charge store |
| US6534954B1 (en) * | 2002-01-10 | 2003-03-18 | Compact Power Inc. | Method and apparatus for a battery state of charge estimator |
| JP3867581B2 (en) * | 2002-01-17 | 2007-01-10 | 松下電器産業株式会社 | Assembled battery system |
| DE10231700B4 (en) * | 2002-07-13 | 2006-06-14 | Vb Autobatterie Gmbh & Co. Kgaa | Method for determining the aging state of a storage battery with regard to the removable amount of charge and monitoring device |
| DE10240329B4 (en) * | 2002-08-31 | 2009-09-24 | Vb Autobatterie Gmbh & Co. Kgaa | Method for determining the charge quantity of a storage battery and monitoring device for a storage battery that can be taken from a fully charged storage battery |
| DE10252760B4 (en) * | 2002-11-13 | 2009-07-02 | Vb Autobatterie Gmbh & Co. Kgaa | Method for predicting the internal resistance of a storage battery and monitoring device for storage batteries |
| WO2004068157A1 (en) * | 2003-01-30 | 2004-08-12 | Robert Bosch Gmbh | State variable and parameter estimator comprising several partial models for an electrical energy storage device |
| US7317300B2 (en) * | 2003-06-23 | 2008-01-08 | Denso Corporation | Automotive battery state monitor apparatus |
| JP4045340B2 (en) * | 2003-08-13 | 2008-02-13 | 現代自動車株式会社 | Battery effective power calculation method and calculation system |
| US6927554B2 (en) * | 2003-08-28 | 2005-08-09 | General Motors Corporation | Simple optimal estimator for PbA state of charge |
| US7109685B2 (en) * | 2003-09-17 | 2006-09-19 | General Motors Corporation | Method for estimating states and parameters of an electrochemical cell |
| US7039534B1 (en) * | 2003-11-03 | 2006-05-02 | Ryno Ronald A | Charging monitoring systems |
| US7321220B2 (en) * | 2003-11-20 | 2008-01-22 | Lg Chem, Ltd. | Method for calculating power capability of battery packs using advanced cell model predictive techniques |
| CN100570388C (en) * | 2003-12-18 | 2009-12-16 | 株式会社Lg化学 | Apparatus and method for evaluating state of charge of a battery using a neural network |
| JP4583765B2 (en) * | 2004-01-14 | 2010-11-17 | 富士重工業株式会社 | Remaining capacity calculation device for power storage device |
| US7233128B2 (en) * | 2004-07-30 | 2007-06-19 | Ford Global Technologies, Llc | Calculation of state of charge offset using a closed integral method |
| US7197487B2 (en) * | 2005-03-16 | 2007-03-27 | Lg Chem, Ltd. | Apparatus and method for estimating battery state of charge |
| US7498772B2 (en) * | 2006-04-06 | 2009-03-03 | International Truck Intellectual Property Company, Llc | Method and system of modeling energy flow for vehicle battery diagnostic monitoring |
| KR100804698B1 (en) * | 2006-06-26 | 2008-02-18 | 삼성에스디아이 주식회사 | The method of assuming the state of charge of the battery, battery management system using the method and the driving method of the battery management system using the method |
-
2004
- 2004-11-11 US US10/985,617 patent/US8103485B2/en active Active
- 2004-11-29 JP JP2007541083A patent/JP4511600B2/en not_active Expired - Lifetime
- 2004-11-29 RU RU2007116586/09A patent/RU2361333C2/en active
- 2004-11-29 BR BRPI0419118A patent/BRPI0419118B1/en active IP Right Grant
- 2004-11-29 TW TW093136727A patent/TWI268637B/en not_active IP Right Cessation
- 2004-11-29 WO PCT/KR2004/003103 patent/WO2006052043A1/en not_active Ceased
- 2004-11-29 KR KR1020077010689A patent/KR100818520B1/en not_active Expired - Lifetime
- 2004-11-29 CA CA2585921A patent/CA2585921C/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| US20060100833A1 (en) | 2006-05-11 |
| KR100818520B1 (en) | 2008-03-31 |
| KR20070069197A (en) | 2007-07-02 |
| JP2008519977A (en) | 2008-06-12 |
| TWI268637B (en) | 2006-12-11 |
| RU2361333C2 (en) | 2009-07-10 |
| TW200616270A (en) | 2006-05-16 |
| BRPI0419118B1 (en) | 2019-08-27 |
| US8103485B2 (en) | 2012-01-24 |
| WO2006052043A1 (en) | 2006-05-18 |
| JP4511600B2 (en) | 2010-07-28 |
| RU2007116586A (en) | 2008-12-20 |
| CA2585921A1 (en) | 2006-05-18 |
| BRPI0419118A (en) | 2008-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2585921C (en) | State and parameter estimation for an electrochemical cell | |
| CA2588856C (en) | Method and system for battery state and parameter estimation | |
| US7593821B2 (en) | Method and system for joint battery state and parameter estimation | |
| US7315789B2 (en) | Method and system for battery parameter estimation | |
| US6534954B1 (en) | Method and apparatus for a battery state of charge estimator | |
| Wang et al. | Model‐based unscented Kalman filter observer design for lithium‐ion battery state of charge estimation | |
| US8035345B2 (en) | System, method, and article of manufacture for determining an estimated combined battery state-parameter vector | |
| Mastali et al. | Battery state of the charge estimation using Kalman filtering | |
| US7197487B2 (en) | Apparatus and method for estimating battery state of charge | |
| US7965059B2 (en) | System, method, and article of manufacture for determining an estimated battery parameter vector | |
| US20070236182A1 (en) | System and method for determining both an estimated battery state vector and an estimated battery parameter vector | |
| US20080162059A1 (en) | Secondary Battery Charge/Discharge Electricity Amount Estimation Method and Device, Secondary Battery Polarization Voltage Estimation Method and Device and Secondary Battery Remaining Capacity Estimation Method and Device | |
| EP4212896A1 (en) | Method for estimating state of charge of battery | |
| CA2588334C (en) | Method and system for joint battery state and parameter estimation | |
| KR100878123B1 (en) | Method and system for battery state and parameter estimation | |
| US12571845B2 (en) | Method and apparatus for calculating relative state-of-charge of battery | |
| KR100916510B1 (en) | Joint Battery Status and Parameter Estimation System and Method | |
| Liu et al. | Li-ion battery state of charge estimation of multi-type working conditions using a segmented multiple independent forgetting factors recursive least squares method | |
| US12618906B2 (en) | Method for estimating state of charge of battery |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request |