CA2547012C - Method for calculating power capability of battery packs using advanced cell model predictive techniques - Google Patents
Method for calculating power capability of battery packs using advanced cell model predictive techniques Download PDFInfo
- Publication number
- CA2547012C CA2547012C CA002547012A CA2547012A CA2547012C CA 2547012 C CA2547012 C CA 2547012C CA 002547012 A CA002547012 A CA 002547012A CA 2547012 A CA2547012 A CA 2547012A CA 2547012 C CA2547012 C CA 2547012C
- Authority
- CA
- Canada
- Prior art keywords
- battery
- charge
- cell
- state
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/3644—Constructional arrangements
- G01R31/3648—Constructional arrangements comprising digital calculation means, e.g. for performing an algorithm
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J7/00—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
- G01R31/3828—Arrangements for monitoring battery or accumulator variables, e.g. SoC using current integration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/374—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC] with means for correcting the measurement for temperature or ageing
Abstract
The present invention relates to a method and an apparatus for estimating discharge and charge power of battery applications, including battery packs used in Hybrid Electric Vehicles (HEV) and Electric Vehicles (EV). One charge/discharge power estimating method incorporates voltage, state-of-charge (SOC), power, and current design constraints and works for a user-specified prediction time horizon of ¥Ät. At least two cell models are used in calculating maximum charge/discharge power based on voltage limits. The first is a simple cell model that uses a Taylor-series expansion to linearize the equation involved. The second is a more complex and accurate model that models cell dynamics in discrete-time state-space form. The cell model can incorporate a inputs such as temperature, resistance, capacity, etc. One advantage of using model-based approach is that the same model may be used in both Kalman-filtering to produce the SOC and the estimation of maximum charge/discharge current based on voltage limits.
Description
METHOD FOR CALCULATING POWER CAPABILITY OF BATTERY PACKS
USING ADVANCED CELL MODEL PREDICTIVE TECHNIQUES
Technical Field The present invention relates to the implementation of a method and apparatus for estimating battery charge power and discharge power.
Background Art A number of high-performance battery applications require precise real-time estimates of the power available to be sourced by the battery pack. For example, in Hybrid Electric Vehicles (HEVs) and Battery Electric Vehicles (BEVs), the vehicle controller requires continuous up-to-date information from the Battery Management System (BMS) regarding the power that may be supplied to the electric motor from the battery pack, and power that may be supplied to the pack via regenerative braking or by active recharging via the motor. One current technique in the art, called the HPPC (Hybrid Pulse Power Characterization) method, performs this task of estimation by using the voltage limits to calculate the maximum charge and discharge limits. As described in the PNGV (Partnership for New Generation Vehicles) Battery Test Manual, Revision 3, February 2001, published by the Idaho National Engineering and Environment Laboratory of the U.S. Department of Energy, the HPPC method estimates maximum cell power by considering only operational design limits on voltage. It does not consider design limits on current, power, or the battery state-of-charge (SOC). Also the method produces a crude prediction for horizon fit. Each cell in the battery pack is modeled by the approximate relationship vk(t) = OCV (zk(t)) - R x ik(t), ( 1 ) where OCV(zk(t)) is the open-circuit-voltage of cell k at its present state-of-charge (zk(t)) and R is a constant representing a cell's internal resistance. Different values of R may be used for charge and discharge currents, if desired, and are denoted as R~''g and Rdls, respectively.
Since the design limits vmin -_< vk(t) < vI"a" must be enforced, the maximum discharge current may be calculated as constrained by voltage, as shown below dis,volt ~C'V (zk (t)) - 'Umin Zmax,k - Rdis The maximum magnitude charge current may be similarly calculated based on voltage. Note, however, that charge current is assumed negative in sign by convention employed in the present invention (although the opposite convention may be used with minor modifications to the method), so that maximum-magnitude current is a minimum in the signed sense. It is chg,volt OCV (zk(t)) - vmax Zmin,k ' Rchg ~ ( Pack power is then calculated as chg ( chg,volt) Pmin - nsnp Trig JC 'UTnaxZmin,k Pdis - n n ITlkI1 (v Zdis,wlt) .
max s p min m~,k This prior art charge calculation method is limited in several respects. First, as noted above, the method does not use operational design limits in SOC, maximum current, or maximum power in the computation. More importantly, the cell model used is too primitive to give precise results.
Overly optimistic or pessimistic values could be generated, either posing a safety of battery-health hazard or causing inefficient battery use.
What is desired is a new method and apparatus for battery charge estimation based on a battery cell model.
Such a cell model would be combined with a maximum-power algorithm that uses the cell model to give better power prediction. The new method would also take in operational design limits such as SOC, current, and power.
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 where:
Fig. 1A is a flow chart that outlines the maximum discharge estimation according to an embodiment of the present invention;
Fig. 1B is a flow chart that outlines the minimum charge estimation according to an embodiment of the present invention;
Fig. 2 is a schematic block diagram showing the sensor components of a power estimating embodiment of the present invention;
USING ADVANCED CELL MODEL PREDICTIVE TECHNIQUES
Technical Field The present invention relates to the implementation of a method and apparatus for estimating battery charge power and discharge power.
Background Art A number of high-performance battery applications require precise real-time estimates of the power available to be sourced by the battery pack. For example, in Hybrid Electric Vehicles (HEVs) and Battery Electric Vehicles (BEVs), the vehicle controller requires continuous up-to-date information from the Battery Management System (BMS) regarding the power that may be supplied to the electric motor from the battery pack, and power that may be supplied to the pack via regenerative braking or by active recharging via the motor. One current technique in the art, called the HPPC (Hybrid Pulse Power Characterization) method, performs this task of estimation by using the voltage limits to calculate the maximum charge and discharge limits. As described in the PNGV (Partnership for New Generation Vehicles) Battery Test Manual, Revision 3, February 2001, published by the Idaho National Engineering and Environment Laboratory of the U.S. Department of Energy, the HPPC method estimates maximum cell power by considering only operational design limits on voltage. It does not consider design limits on current, power, or the battery state-of-charge (SOC). Also the method produces a crude prediction for horizon fit. Each cell in the battery pack is modeled by the approximate relationship vk(t) = OCV (zk(t)) - R x ik(t), ( 1 ) where OCV(zk(t)) is the open-circuit-voltage of cell k at its present state-of-charge (zk(t)) and R is a constant representing a cell's internal resistance. Different values of R may be used for charge and discharge currents, if desired, and are denoted as R~''g and Rdls, respectively.
Since the design limits vmin -_< vk(t) < vI"a" must be enforced, the maximum discharge current may be calculated as constrained by voltage, as shown below dis,volt ~C'V (zk (t)) - 'Umin Zmax,k - Rdis The maximum magnitude charge current may be similarly calculated based on voltage. Note, however, that charge current is assumed negative in sign by convention employed in the present invention (although the opposite convention may be used with minor modifications to the method), so that maximum-magnitude current is a minimum in the signed sense. It is chg,volt OCV (zk(t)) - vmax Zmin,k ' Rchg ~ ( Pack power is then calculated as chg ( chg,volt) Pmin - nsnp Trig JC 'UTnaxZmin,k Pdis - n n ITlkI1 (v Zdis,wlt) .
max s p min m~,k This prior art charge calculation method is limited in several respects. First, as noted above, the method does not use operational design limits in SOC, maximum current, or maximum power in the computation. More importantly, the cell model used is too primitive to give precise results.
Overly optimistic or pessimistic values could be generated, either posing a safety of battery-health hazard or causing inefficient battery use.
What is desired is a new method and apparatus for battery charge estimation based on a battery cell model.
Such a cell model would be combined with a maximum-power algorithm that uses the cell model to give better power prediction. The new method would also take in operational design limits such as SOC, current, and power.
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 where:
Fig. 1A is a flow chart that outlines the maximum discharge estimation according to an embodiment of the present invention;
Fig. 1B is a flow chart that outlines the minimum charge estimation according to an embodiment of the present invention;
Fig. 2 is a schematic block diagram showing the sensor components of a power estimating embodiment of the present invention;
Fig. 3 is an example plot of open-circuit-voltage (OCV) as a function of state-of-charge for a particular cell electrochemistry;
Fig. 4 is an example plot showing the derivative of OCV
as a function of state-of-charge for a particular cell electrochemistry;
Fig. 5 is a plot showing the voltage prediction using the cell model of the present invention;
Fig 6 is a zoom-in of the plot of voltage prediction for one UDDS cycle at around 50o state-of-charge;
Fig. 7 is a state-of-charge trace for cell test;
Fig. 8 is a plot comparing static maximum power calculations as functions of SOC for the PNGV HPPC method and Method I of the present invention;
Fig. 9 is a plot showing that discharge power capability estimates for cell cycle test comprising sixteen UDDS
cycles over an SOC range of 90~ down to 100;
Fig. 10 is zoomed-in plot of Fig. 9, showing about one UDDS cycle;
Fig. 11 is a plot showing charging power capability estimates for cell cycle test comprising sixteen UDDS
cycles over an SOC range of 90o down to 10%; and Fig. 12 is zoomed-in plot of Fig. 11, showing about one UDDS cycle.
Fig. 4 is an example plot showing the derivative of OCV
as a function of state-of-charge for a particular cell electrochemistry;
Fig. 5 is a plot showing the voltage prediction using the cell model of the present invention;
Fig 6 is a zoom-in of the plot of voltage prediction for one UDDS cycle at around 50o state-of-charge;
Fig. 7 is a state-of-charge trace for cell test;
Fig. 8 is a plot comparing static maximum power calculations as functions of SOC for the PNGV HPPC method and Method I of the present invention;
Fig. 9 is a plot showing that discharge power capability estimates for cell cycle test comprising sixteen UDDS
cycles over an SOC range of 90~ down to 100;
Fig. 10 is zoomed-in plot of Fig. 9, showing about one UDDS cycle;
Fig. 11 is a plot showing charging power capability estimates for cell cycle test comprising sixteen UDDS
cycles over an SOC range of 90o down to 10%; and Fig. 12 is zoomed-in plot of Fig. 11, showing about one UDDS cycle.
Disclosure of the Invention The present invention relates to a method and an apparatus for estimating discharge and charge power of battery applications, including battery packs used in Hybrid Electric Vehicles (HEVs) and Electric Vehicles (EVs).
One embodiment is a charge prediction method that incorporates voltage, state-of-charge, power, and current design constraints, works for a user-specified prediction horizon Ot, and is more robust and precise than the state of the art. The embodiment has the option of allowing different modeling parameters during battery operation to accommodate highly dynamic batteries used in Hybrid Electric Vehicles (HEVs) and Electric Vehicles (EVs) where such previous implementations were difficult.
An embodiment of the present invention calculates maximum charge/discharge power by calculating the maximum charge/discharge current using any combination of four primary limits:
1. state-of-charge (SOC) limits 2. voltage limits 3. current limits 4. power limits In one embodiment, the minimum absolute charge/discharge current value from the calculations using state-of-charge (SOC), voltage, and current limits is then chosen to obtain the maximum absolute charge/discharge power. In one embodiment, the maximum absolute charge/discharge power is checked to ensure it is within the power limits. In one embodiment, the maximum absolute charge/discharge power is calculated in a way as to not violate any combination of the limits that may be used.
Prior methods do not use SOC limits in their estimation of maximum charge/discharge power. The present invention incorporates the SOC of the battery cell or battery pack to estimate the maximum charge/discharge current. The estimation explicitly includes a user-defined time horizon Ot. In one embodiment, the SOC is obtained by using a Kalman filter. The SOC that is produced by Kalman filtering also yields an estimate of the uncertainty value, which can be used in the maximum charge/discharge calculation to yield a confidence level of the maximum charge/discharge current estimate.
Methods of the present inventions improve upon prior art estimation of power based on voltage limits. In the present invention, voltage limits are used to calculate the maximum charge/discharge current in a way that includes a user-defined time horizon 0t. Two primary cell model embodiments are in the present invention for the calculation of maximum charge/discharge power based on voltage limits. The first is a simple cell model that uses a Taylor-series expansion to linearize the equation involved. The second is a more complex and accurate cell model that models cell dynamics in discrete-time state-space form. The cell model can incorporate a variety of inputs such as temperature, resistance, capacity, etc. One advantage of using model-based approach is that the same cell model may be used in both Kalman filtering to produce the SOC and the estimation of maximum charge/discharge current based of voltage limits.
Embodiments of the present invention also include methods of charge estimation based on any combination of the voltage, current, power, or SOC limits described above.
For example, charge estimation can be based on voltage limits only, or combined with current limits, SOC limits and/or power limits.
Embodiments of the present invention are directed to a power apparatus that takes in data measurements from the battery such as current, voltage, temperature, and feeding such measurements to an arithmetic circuit, which includes calculation means that performs the calculation methods disclosed in the present invention to estimate the absolute maximum charge or discharge power.
Best Mode for Carrying Out the Invention Embodiments of the present invention relates to battery charge estimation for any battery-powered application. In one embodiment, the estimator method and apparatus find the maximum absolute battery charge and/or discharge power (based on present battery pack conditions) that may be maintained for Ot seconds without violating pre-set limits on cell voltage, state-of-charge, power, or current.
Figs. 1A and 1B illustrates an overview of the embodiments of the present invention. Fig. 1A shows a method for finding the maximum discharge power for a user-defined time horizon 0t, i.e. how much power may be drawn from the battery continuously for use for the next 0t time period. In vehicle applications, accurate estimation of maximum discharge power can help prevent the hazardous occurrence of over-drawing the battery.
In step 10, the maximum discharge current is calculated based on pre-set limits on state-of-charge. The estimation explicitly includes a user-defined time horizon Ot. In one embodiment, the SOC is obtained by using a kalman filtering method. The SOC that is produced by Kalman filtering also yields an estimate of the uncertainty value, which can be .
5' used in the maximum charge/discharge calculation to yield a confidence level of the maximum charge/discharge current estimation. In another embodiment, a simple state-of-charge is used. Step 10 is further described in the section titled "Calculation Based on State-of-Charge (SOC) Limits."
The maximum discharge current is calculated based on pre-set limits on voltage in step 12. The present invention has two main model embodiments for calculation of maximum charge/discharge power based on voltage limits, although it is understood that other cell models could be used. Both overcome the limitation of prior art discharge estimation methods of giving a crude prediction of time horizon fit.
The first is a simple cell model that uses a Taylor-series expansion to linearize the equation involved. The second is a more complex and accurate cell model that models cell dynamics in discrete-time state-space form. The cell model can incorporate a variety of inputs such as temperature, resistance, capacity, etc. The two cell models are further described in the section titled "Calculation Based on Voltage Limits."
Then In step 14, the maximum discharge current is calculated based on pre-set limits on current. In step 16, the minimum of the three calculated current values from steps 10, 12, and 14 is chosen. It is understood that the execution order of steps 10, 12, 14 is interchangeable. It is further understood that any combination of steps 10, 12, and 14 may be omitted, if desired, in .an implementation.
Using the chosen discharge current value, step 18 calculates the maximum discharge power. The calculated pack power may be further refined in order to not violate individual cell or battery pack power design limits.
Fig. 1B shows a method for finding the maximum absolute charge power for a user-defined time horizon 0t, i.e. how much power can be put back into the battery continuously for the next 0t time period. The details and progression of the method mirror that of Fig. 1A. Since charge current is considered to have a negative sign, the maximum absolute current is the minimum current in a signed sense. In step 20, the minimum charge current is calculated based on pre-set limits on state-of-charge. Again the SOC can be a simple one or one obtained using the Kalman filtering method. Then the minimum charge current is calculated based on pre-set limits on voltage in step 22 in accordance with a cell model, such as one of the two cell models described . in the present disclosure. Then in step 24, the minimum charge current is calculated based on pre-set limits on current. Then, in step 26, the maximum of the three calculated current values from steps 20, 22, 24 is chosen.
Note again that the execution order of steps 20, 22, 24 is interchangeable. It is further understood that any combination of steps 20, 22, 24 may be used and any of the steps 20, 22, 24 may be omitted, if desires, in an implementation. Using the chosen charge current value, step 28 calculates the minimum charge power. The calculated pack power may be further refined in order to not violate individual cell or battery pack power design limits.
One embodiment is a charge prediction method that incorporates voltage, state-of-charge, power, and current design constraints, works for a user-specified prediction horizon Ot, and is more robust and precise than the state of the art. The embodiment has the option of allowing different modeling parameters during battery operation to accommodate highly dynamic batteries used in Hybrid Electric Vehicles (HEVs) and Electric Vehicles (EVs) where such previous implementations were difficult.
An embodiment of the present invention calculates maximum charge/discharge power by calculating the maximum charge/discharge current using any combination of four primary limits:
1. state-of-charge (SOC) limits 2. voltage limits 3. current limits 4. power limits In one embodiment, the minimum absolute charge/discharge current value from the calculations using state-of-charge (SOC), voltage, and current limits is then chosen to obtain the maximum absolute charge/discharge power. In one embodiment, the maximum absolute charge/discharge power is checked to ensure it is within the power limits. In one embodiment, the maximum absolute charge/discharge power is calculated in a way as to not violate any combination of the limits that may be used.
Prior methods do not use SOC limits in their estimation of maximum charge/discharge power. The present invention incorporates the SOC of the battery cell or battery pack to estimate the maximum charge/discharge current. The estimation explicitly includes a user-defined time horizon Ot. In one embodiment, the SOC is obtained by using a Kalman filter. The SOC that is produced by Kalman filtering also yields an estimate of the uncertainty value, which can be used in the maximum charge/discharge calculation to yield a confidence level of the maximum charge/discharge current estimate.
Methods of the present inventions improve upon prior art estimation of power based on voltage limits. In the present invention, voltage limits are used to calculate the maximum charge/discharge current in a way that includes a user-defined time horizon 0t. Two primary cell model embodiments are in the present invention for the calculation of maximum charge/discharge power based on voltage limits. The first is a simple cell model that uses a Taylor-series expansion to linearize the equation involved. The second is a more complex and accurate cell model that models cell dynamics in discrete-time state-space form. The cell model can incorporate a variety of inputs such as temperature, resistance, capacity, etc. One advantage of using model-based approach is that the same cell model may be used in both Kalman filtering to produce the SOC and the estimation of maximum charge/discharge current based of voltage limits.
Embodiments of the present invention also include methods of charge estimation based on any combination of the voltage, current, power, or SOC limits described above.
For example, charge estimation can be based on voltage limits only, or combined with current limits, SOC limits and/or power limits.
Embodiments of the present invention are directed to a power apparatus that takes in data measurements from the battery such as current, voltage, temperature, and feeding such measurements to an arithmetic circuit, which includes calculation means that performs the calculation methods disclosed in the present invention to estimate the absolute maximum charge or discharge power.
Best Mode for Carrying Out the Invention Embodiments of the present invention relates to battery charge estimation for any battery-powered application. In one embodiment, the estimator method and apparatus find the maximum absolute battery charge and/or discharge power (based on present battery pack conditions) that may be maintained for Ot seconds without violating pre-set limits on cell voltage, state-of-charge, power, or current.
Figs. 1A and 1B illustrates an overview of the embodiments of the present invention. Fig. 1A shows a method for finding the maximum discharge power for a user-defined time horizon 0t, i.e. how much power may be drawn from the battery continuously for use for the next 0t time period. In vehicle applications, accurate estimation of maximum discharge power can help prevent the hazardous occurrence of over-drawing the battery.
In step 10, the maximum discharge current is calculated based on pre-set limits on state-of-charge. The estimation explicitly includes a user-defined time horizon Ot. In one embodiment, the SOC is obtained by using a kalman filtering method. The SOC that is produced by Kalman filtering also yields an estimate of the uncertainty value, which can be .
5' used in the maximum charge/discharge calculation to yield a confidence level of the maximum charge/discharge current estimation. In another embodiment, a simple state-of-charge is used. Step 10 is further described in the section titled "Calculation Based on State-of-Charge (SOC) Limits."
The maximum discharge current is calculated based on pre-set limits on voltage in step 12. The present invention has two main model embodiments for calculation of maximum charge/discharge power based on voltage limits, although it is understood that other cell models could be used. Both overcome the limitation of prior art discharge estimation methods of giving a crude prediction of time horizon fit.
The first is a simple cell model that uses a Taylor-series expansion to linearize the equation involved. The second is a more complex and accurate cell model that models cell dynamics in discrete-time state-space form. The cell model can incorporate a variety of inputs such as temperature, resistance, capacity, etc. The two cell models are further described in the section titled "Calculation Based on Voltage Limits."
Then In step 14, the maximum discharge current is calculated based on pre-set limits on current. In step 16, the minimum of the three calculated current values from steps 10, 12, and 14 is chosen. It is understood that the execution order of steps 10, 12, 14 is interchangeable. It is further understood that any combination of steps 10, 12, and 14 may be omitted, if desired, in .an implementation.
Using the chosen discharge current value, step 18 calculates the maximum discharge power. The calculated pack power may be further refined in order to not violate individual cell or battery pack power design limits.
Fig. 1B shows a method for finding the maximum absolute charge power for a user-defined time horizon 0t, i.e. how much power can be put back into the battery continuously for the next 0t time period. The details and progression of the method mirror that of Fig. 1A. Since charge current is considered to have a negative sign, the maximum absolute current is the minimum current in a signed sense. In step 20, the minimum charge current is calculated based on pre-set limits on state-of-charge. Again the SOC can be a simple one or one obtained using the Kalman filtering method. Then the minimum charge current is calculated based on pre-set limits on voltage in step 22 in accordance with a cell model, such as one of the two cell models described . in the present disclosure. Then in step 24, the minimum charge current is calculated based on pre-set limits on current. Then, in step 26, the maximum of the three calculated current values from steps 20, 22, 24 is chosen.
Note again that the execution order of steps 20, 22, 24 is interchangeable. It is further understood that any combination of steps 20, 22, 24 may be used and any of the steps 20, 22, 24 may be omitted, if desires, in an implementation. Using the chosen charge current value, step 28 calculates the minimum charge power. The calculated pack power may be further refined in order to not violate individual cell or battery pack power design limits.
It is noted that modifications may be made to the method embodiments as shown in Figs. 1A and 1B. For example, any or all of the current calculation steps based on state-of-charge and voltage limits may be removed. Also, the present invention discloses several methods of calculating maximum absolute charge and discharge current based on state-of-charge, voltage limits, and current limits.
One embodiment of the present invention estimates the maximum absolute charge and/or discharge power of a battery pack. The battery pack may be, for example, a battery pack used in a hybrid electric vehicle or an electric vehicle.
The embodiment makes a number of denotations and limits, including:
- using n to denote the number of cells in the target battery pack where an estimation of charge and/or discharge power is desired;
- using vk(t) to denote the cell voltage for cell number k in the pack, which has operational design limits so that k:l<k<_n;
umiri < vk(t) C ~,'maxmust be enforced for all - -- using zk(t) to denote the state-of-charge for cell number k in the pack, which has operational design limits < z t < zmax zmin -_ ~( ) that must be enforced for all ~ ~ 1 < k < n;
- using pk(t) to denote the cell power, which has a operational design limits so that p''''~~n ~ pk(t) ~ ~m~' must be enforced for all ~ ~ 1 ~ 1~ < n; and - using ik(t) to denote the cell current, which has a operational design limits so that 2min ~ 2x (t) C Zmax must be enforced for all ~ ~ 1 < k < n.
Modifications can be made in alternate embodiments. For example, any particular limit may be removed if desired by replacing its value by ~°°, as appropriate. As an another example, 11m1tS SLlCh aS Vmaxr Vminr Zmaxr Zminr lmaxr lmin, pmaxr pmin may furthermore be functions of temperature and other factors pertaining to the present battery pack operating condition. In one embodiment, it is assumed that the discharge current and power have positive sign and the charge current and power have negative sign. Those skilled in the art will recognize that other sign conventions may be used, and that the description of the present invention can be adapted to these conventions in a forthright manner.
In one embodiment, the model used for predicting charge assumes that the battery pack comprises ns cell modules connected in series, where each cell module comprises np individual cells connected in parallel and ns ~ 1, n~ > 1.
Other configurations are possible and are accommodated by slight modifications to the method as described.
Fig. 2 is a schematic block diagram showing the sensor components of an embodiment of the present invention.
Battery 40 is connected to load circuit 48. For example, load circuit 48 could be a motor in an Electric Vehicle (EV) or Hybrid Electric Vehicle (HEV). In some embodiments circuit 48 is a circuit that provides power and/or draws power. Measurements of battery and individual cell voltage are made with voltmeters) 44. Measurements of battery current are made with ammeter 42. Battery and individual cell temperatures are measured by temperature sensors) 46.
Voltage, current and temperature measurements are processed with arithmetic circuit 50. Arithmetic circuit (estimator means) 50, takes in the measurements from the sensor components and perform the calculation methods of the present invention for power estimation. In some embodiment, temperature is not needed in the calculation methods.
1. Calculation Based in State-of-Charge (SOC) Limits As shown in steps 10 and 20 of Figs 1A and 1B, embodiments of the present invention calculate the maximum charge/discharge current values using SOC limits. Various embodiments also have the explicit inclusion of a time horizon Ot in the calculation. The SOC limits are included as follows. First, for a constant current ik, the SOC
recurrent relationship is described as:
zk(t -E- Ot) = zk (t) - (~hOt~C)zk~ ( q ) Where zk ( t) is the present SOC for cell k, zk ( t+0t) is the predicted SOC ~t seconds into the future, C is the cell capacity in ampere-seconds, and ~i is the Coulombic efficiency factor at current level ik. Here, for simplicity of presentation, it is assumed that ~i =1 =or aiscnarge currents and ~~ - ~ < 1 for charge currents .
If there are design limits on SOC such that zmin < zk(t) < zmax for all cells in the pack, then current ik can be computed such that these limits are not exceeded.
Simple algebra gives limits based on the SOC of each cell:
2misa s k - zk (t~ - zcnin chg,soc zk (t~ - zmax Zmin,k - ~Qt~G' The pack maximum absolute currents-based only on cell SOC-are then dis,soc ~ dis,socl Zmax - ~k I1 2m~ Jk chg,soc chg,soc Zcnin - Illk 7C CZmin,k ) .
This method assumes that there is a valid SOC estimate available for every cell in the pack. If there is not the case, then an approximate remedy would be to calculate dis,soc ,~, z(t~ - zmin 2m~ ~ Ot/Cr chg,soc ," z~t) - zmax Zmin '" ~Ot~Cr where z(t) is the pack SOC.
In one embodiment of the present invention, the power predictive method can take into account more information than simply the cell SOC. For example, a Kalman filter can be used as a method to estimate all the cell SOCs in a pack.
Besides giving the SOC, Kalman filtering yields estimates of the uncertainty of the SOC estimate itself. A method of using Kalman filter to estimate SOC is described in commonly assigned U.S. Patent No. 6,534,954, hereby incorporated by reference.
Let the uncertainty have Gaussian distribution with standard deviation, as estimated by the Kalman filter, be denoted as oz. Then, the method yields a95.5o confidence that the true SOC is within the estimate ~2 aZ and a 99.70 confidence that the true SOC is within the estimate ~3 aZ.
This information can be incorporated into the estimate of~ maximum current based on SOC to have very high confidence that SOC design limits will not be violated.
This is done as (assuming a 3 aZ confidence interval):
dis,soc (zk (t) - 3~z) - zmin 2max,k Ot/rr chg,soc (zk(t) + 3Qx) - zmax Zmin,k - ~~t/C ' 2. Calculation Based on Voltage Limits Besides taking SOC limits into account, embodiments of the present invention correct a limitation in the prior art HPPC method for applying voltage limits (steps 12 and 22 of Figs. 1A and 1B). In the HPPC method, if the cell model of equation (1) is assumed, and that Rc''g and Rdls are the cell's Ohmic resistances, then equation (2) and equation (3) predict the instantaneously available current, not the constant value of current that is available for the next Ot seconds. If cases where ~t is large, the result of the calculation poses a safety or battery-health issue, as the cells may become over/under charged.
To overcome this problem, an embodiment of the present invention uses the following cell model:
vk(t + Ot) = OCV(zk(t + Ot)) - R x ik(t), (7 ) This modifies the previous cell model in equation (1).
Note that this model cannot be directly solved in closed form for the maximum current ik since zk ( t+~t) is itself a function of current (cf.(4)) and OCV(~) is a nonlinear relationship. Note that other cell models can be used as well.
Two method embodiments are directed to solving (7) for the maximum absolute value of ik(t).
2.lMethod I: Taylor-series expansion The first method uses a Taylor-series expansion to linearize the equation, so that at approximate value of i can be solved. It is assumed that OCV(~) is differentiable at point zk(t), which gives the result OCV(zk(t + ~1t)) - OCV Czk(t) - ik~C t) - OCV(zk(t)) - ik ~i~t 80CV(z) C ~ 8z I _ t z-xk ( ) + Rl(zk(t), Zk(~i~t)I C) where the first-order residual Rl(zk(t)' Zk(~i0t)~C)~Ilzk(t)II -~
as 2k(~JaOt)~C ~ ~ in 0. Note that the change in SOC over Ot seconds is generally small, so the following approximation may be applied OCV (zk (t + 4t)) = OCV (zk (t)) - ik ~C t a0 ~z (z) x=z t k( ) -Rxik which gives dis,volt OCV(zk(t)) - vm;n 2max,k ~t 80CV x + Rdis C 8z xk (t) ( g ) chg,volt ~CV (zk (t)) - ymax 2min,k - nOt 80CV z ~ cbg C 8x xk~t) + R ( 9 ) In one embodiment, both the function OCV(z) and its derivative aOCV(z)/ 8z might be computed from some known mathematical relationship for OCV(z), (e. g., Nernst's equation) using either analytic or numeric methods, or by a table lookup of empirical data. This quantity is positive for most battery electrochemistries over the entire SOC
range, so the values computed by (8) and (9) are smaller in magnitude than those from ( 2 ) and ( 3 ) for the same values of Rdis and Ring.
The HPPC procedure compensates for its inaccuracy by using modified values of Rdls and R'''g, determined experimentally, that approximate the denominator terms in ( 8 ) and ( 9 ) . This can not be accurate over the entire SOC
range, however, as 80CV(z)/ az is not constant, particularly near extreme values of z.
Discharge and charge currents with all limits enforced are computed as (steps 16 and 26 of Figs. 1A and 1B) dis r dis,soc dis,voltl Zmax - min ~2rny ~lIl Zm~ k ~ IIkIl 2max,k (10) im n - Tna.X C2rnin ~ Ink x im n ~°c ~ IIIkaJC 2m n' k ltJ
(11) and power may be calculated using the sum of all cell powers. These are equal to the product of the maximum allowed current and the predicted future voltage.
'~a chg ~ c6g pmin - nP Zminvk ~t + Ot k=1 na r l N 7LP ~ Zm a ~~iV' zk ~t~ - Zm n~iOt~C~
k= \1 - Rchg X 2m n _ i (12) Tl.a Pmtax - nP ~ Zm'axvk~t -~- ~t~.
k=1 .
na 1 7ZP ~ Zrdc,~ax ~G'v C,Zk(t) - 2m~~Ot~~i~
k= /1 - Rdis X idia max (13) Maximum and minimum cell and pack power limits may also be imposed in this calculation. Note that in all equations, 0CV ( Z) , C, Vmaxr Vminr Zmaxr Zminr lmaxr lmin, Rchgr and Rdls may be functions of temperature and other factors pertaining to the present battery pack operating conditions.
2.2 Method II: Using a Comprehensive Cell Model The method of solving (7) presented in the previous section requires less computational intensity. A second method embodiment of the present invention may be used when more computational power is available. This second method assumes a more precise mathematical model of cell dynamics, which might be in a discrete-time state-space form such as the coupled pair of equations xk[m + 1] - f (~k[m], uk[m]) ( 14 ) vk [m] - 9 ~xk [m]' uk [m]~' ( 15 ) where m is the discrete time sample index; the vector function of time xk[m] is called the estate" of the system, uk[m] is the input to the system, which includes cell current ik[m] as a component, and might also include temperature, resistance, capacity and so forth, and f(~) and g(') are functions chosen to model cell dynamics.
Alternative model forms, including continuous-time state-space forms, differential and difference equations might also be used. It is assumed that there is a method to compute SOC given the model that is implemented.
For convenience of presentation, it is assumed that the cell model is in a discrete-time state-space form. Also assume that 0t seconds may be represented in discrete time as T sample intervals. Then, this model can be used to predict cell voltage 0t seconds into the future by vk(m-t-T~ =g~~k(rn-I-T~,ZGk(7TL+T)), where xk[m+T] may be found by simulating (14) for T time samples. It is assumed that the input remains constant from time index m to m+T, so if temperature change (for example) over this interval is significant, it must be included as part of the dynamics modeled by ( 14 ) and not as a part of the measured input uk[m] .
The method then uses a bisection search algorithm to dis,volt chg,volt find 2max,x and ~min,k by looking for the ik (as a member of the uk vector) that causes equality in vmin - 9 ~xk(m + Z'~ ~ ux(m -I- ~')) ~ or 0 - g ~xk(m + T~, ~k(m -~ T~~ - v~~l~
(16) dis,volt to find 2max,k ~ and by looking for the ik that causes equality in vma,,~ - g~xk(m+T~,uk(m+T~~, or 0 _ g ~xk(m + T), uk(m + T~) - vmax (17) chg,volt to find Zmin,k _ A special case is when the state equation (14) is linear-that is, when xk (m + 1~ = Axk (m~ + Buk (m~
where A and B are constant matrices. The model presented in Section 3, entitled "An Example Cell Model," is an example where this is the case. Then, for input uk constant time m to m+T, leading to xk m + T = AT xk (m~ +~ ~ AT 1 ' B uk.
~=o Most of these terms may be pre-computed without knowledge of uk in order to speed calculation using the bisection algorithm.
dis,soc chg,soc Once the SOC-based current limits 2max,k and zmin,k are computed using (5) and (6), and the voltage-based current dis,volt chg,volt limits Zmax,k and 2min,k are computed using ( 16 ) and ( 17 ) , overall current limits may be computed using ( 10 ) and ( 11 ) (steps 16 and 26 of Figs. 1A and 1B). Power is then computed as na chg ~ chg Pmin - nP 2minvk ~t -~- Ot) k=1 ns chg ( - nP ~ Zmingk~xk(m'+T~~uk)~
k=1 with uk containing 2mn as its value for current, and n' Pmax - nP ~ Zm~~vk \t + Ot k=1 n, dis nP ~ 2max9k ~~k ~7T'G -~- ~~, 2Gk~, k=1 dis with uk containing Zm~ as its value for current.
1.2.1 Bisection search To solve ( 16 ) and ( 17 ) , a method to solve for a root of a nonlinear equation is required. In one embodiment the bisection search algorithm is used for this requirement.
The bisection search algorithm looks for a root of f(x) ( i . e. , a value of x such that f (x) =0) where it is known a priori that the root lies between values x1 < root < x2.
One way of knowing that a root lies in this interval is that the sign of f(xl) is different from the sign of f(xz), Each iteration of the bisection algorithm evaluates the function at the midpoint Xmid = (xi + xz) /2. Based on the sign of the evaluation, either x1 or xz is replaced by Xmid to retain different signs on of f(xl) and f(xz) . It is evident that the uncertainty in the location of the root is halved by this algorithmic step. The bisection algorithm repeats this iteration until the interval between xl and xz, and hence the resolution of the root of f(x) is as small as desired. If E is the desired root resolution, then the algorithm will require at most ~logz ( ~ xl - xz ~ /E) iterations. The bisection method is listed in Listing 1.
1.2.2 Finding maximum/minimum current To determine maximum discharge and charge current for any particular cell, bisection is performed on (16) and (17). Bisection is incorporated in the overall algorithm ad follows. First, three stirnulations are performed to determine cell voltages Ot seconds into the future for cell current ik = 0, ik = lminr and ik = lmax~ If cell voltages are predicted to be between vmi" and Vmax for the maximum discharge rates, then these maximum rates may be used. If the cell voltages, even during rest, are outside of bounds, then set the maximum rates to zero.. Otherwise, the true maximum rate may be found by bisecting between rate equal to zero and its maximum value. Bisection is performed between current limits (lminr 0) or (0, lmax) .
Algorithm 1 The bisection algorithm.
begin{bisection algorithm}
set x1 = first search limit set x2 = second search limit set a = desired resolution on current in bisection output set JMAX = mraximum number of bisection iterations J l iOg2 ( x2 Exl ) let func(~) be function in which to find root set the search interval dx = xa - x1 if ( ( func(xl)>_0) { // constraint: func(xl)<0 dx = -dx x1 = x2 } // root is now between (xl,xa + dx), and func(xl)<0 forj=ltoJMAX{
dx = 0.5 x dx amid = x1 + dx If ( fUnC(xmid) < 0) { x1 = amid }
if (~dx~< e) { return(x~ +0.5 x dx) }
} II have bisected too many times, so return best guess return(xl + 0.5 x dx) end{bisection algorithm}
2. An exam 1e Cell Model An example cell model for the present invention power estimation methods is presented herein, with illustrations given to show the performance of the two methods compared to the prior art PNGV HPPC method. The cell model is a discrete-time state-space model of the form of (14) and (15) that applies to battery cells. The model, named ~~Enhanced Self-Correcting Cell Model," is further described in the article ~~Advances in EKFLiPB SOC Estimation," by the inventor, published in CD-ROM and presented in Proc. 20th Electric Vehicle Symposium (EVS20) in Long Beach CA, (November 2003) and is hereby fully incorporated by reference. It is understood this model is an example model only and that a variety of suitable alternate models can be used.
The ~~Enhanced Self-Correcting Cell Model" includes effects due to open-circuit-voltage, internal resistance, voltage time constants, and hysteresis. For the purpose of example, the 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.
State-of-charge is captured by one state of the model.
This equation is zk(m+ 1~ = z~(m~ - OaOT~C)i,~(m~, where 0 T represents the inter-sample period (in seconds), and C represents the cell capacity (in ampere-seconds).
The time-constants of the cell voltage response are captured by several filter states. If there is let to be of time constants, then fk(m+ 1) = l4ffx(m) + Blik(~).
E 1E8"~ x"p The matrix f may be a diagonal matrix with real-valued entries. If so, the system is stable if all B E lLf~f "1 entries have magnitude less than one. The vector f may simply be set to of ~~l"s. The value of of and the entries in the Af matrix are chosen as part of the system identification procedure to best fit the model parameters to measured cell data.
The hysteresis level is captured by a single state hk (m + 1J - exp ~iZk ~ ~~~T ~~ hk (m]
r~,ix ('mJ'Y~T
+ 1 - exp G' ~~~ S~Uk(m])~
C
C
where Y is the hysteresis rate constant, again found by system identification.
The overall model state is xk (mJ = ( fk ('n-r,J~ hk (rrs) z~ (mJ , where the symbol () is the matrix/vector transpose operator. The state equation for the model is formed by combining all of the individual equations, above. Note, that at each time step, the state equation is linear in the input u~(mJ = ~ ik(m] 1, which speeds the prediction operation.
The output equation that combines the state values to predict cell voltage is vk(mJ = OCV(z~(m~) + Gfk(mJ - Rik(mJ + Mhk(m], where G E ILBIxn~ is a vector of constants that blend the time-constant states together in the output, R is the cell resistance (different values may be used for dis/charge), and M is the maximum hysteresis level.
The open-circuit-voltage as a function of state-of charge for example Lithium Ion Polymer Battery (LiPB) cells is plotted in Fig.3. This is an empirical relationship found by cell testing. First, the cell was fully charged (constant current to 4.2V, constant voltage to 200mA). Then, the cell was discharged at the C/25 rate until fully discharged (3.0V). The cell was then charged at the C/25 rate until the voltage was 4.2V. The low rates were used to minimize the dynamics excited in the cells. The cell voltage as a function of state of charge under discharge and under charge were averaged to compute the OCV. This has the effect of eliminating to the greatest extent possible the presence of hysteresis and ohmic resistance in the final function. For the purpose of computations involving OCV, the final curve was digitized at 200 points and stored in a table. Linear interpolation is used to look up values in the table.
The partial derivative of OCV with respect to SOC for these example cells is plotted in Fig. 4. This relationship was computed by first taking finite differences between points in the OCV plot in Fig. 3 and dividing by the distance between points (i.e., Euler's approximation to a derivative). The resulting data is too noisy to be of practical use, as shown in the gray line of Fig. 4. It was filtered using a zero-phase low-pass filter, resulting in the black line of Fig. 4, which may be used in the power calculation. This relationship is also digitized at 200 points, and linear interpolation into the table of values is used when computations requiring this function are performed.
Other parameters are fit to the cell model. In particular, the model employs four low-pass filter states (nf = 4), a nominal capacity of 7.5 Ah, and an inter-sample interval of 0T = 1s. There is very close agreement between the cell model voltage prediction and the cell true voltage.
This is illustrated in Fig. 5, which is a plot showing the voltage prediction using the cell model of the present invention. For this figure, the cell test was a sequence of sixteen UDDS cycles, performed at room temperature, separated by discharge pulses and five-minute rests, and spread over the 90o to 10o SOC range. The difference between true cell terminal voltage and estimated cell terminal voltage is very small (a root-mean-squared (RMS) voltage estimation error of less than 5mV). To better illustrate the model' s fidelity, refer to the zoom on one UDDS cycle in the 50% SOC region, shown in Fig. 6. The SOC
as a function of time is plotted in Fig. 7, which is a SOC
trace for cell test. The graph shows that SOC increases by about 5% during each UDDS cycle, but is brought down about loo during each discharge between cycles. The entire operating range for these cells (10% SOC to 90o SOC, delineated on the figure as the region between the thin dashed lines) is excited during the cell test.
3. Comparing Maximum Power Calculations The PNGV HPPC power estimation method gives a result that is a function of only SOC. Therefore, it is possible to graph available power versus SOC to summarize the algorithm calculations. The first method proposed (Method I: Taylor Series Expansion Methods) in this patent disclosure is also possible to display in this way.
Estimated power is only a function of SOC, aOCV/8z (also a function of SOC), and static limits on maximum current and power. The second method (Method II: the Comprehensive Cell Model Method), however, dynamically depends on all states of the system. Two systems at the same state of charge, but with different voltage time-constant state values or hysteresis state levels will have different amounts of power available. To compare power computed by the three methods, dynamic tests must be conducted.
For the following result, a pack of LiPB cells is assumed to be with ns - 40 and np - 1. The data to fit models was collected from prototype hand-made cells jointly developed by LG Chem (Daejeon, Korea) and Compact Power Inc.
(Monument, Colorado). Limits for the power calculations are listed in Table 1. Each cell has a nominal capacity of 7.5 Ah, and Ot was ten seconds for both charge and discharge.
Table 1: Parameters for power calculation examples.
Parameter Minimum Maximum vk(t) 3.0V 4.35V
ik (t) -200A 200A
z~(t) 0.1 0.9 pk (t) -o0 00 First, the PNGV HPPC method and Method I of the present invention are compared in Fig. 8, which is a plot comparing static maximum power calculations as functions of SOC for the PNGV HPPC method and Method I of the present invention.
The black curves correspond to charge power, and the gray curves correspond to discharge power. Note that the absolute value of power is plotted to avoid confusion due to sign conventions. Considering first the calculations of charge power, it is evident that the PNGV HPPC method produces similar values to Method I in the mid-SOC range .
The slight differences are due to the fact that the 10-second R~hg value used for the PNGV method and the derivative-modified R~hg for Method I are not identical.
Outside the mid-SOC range, the graph shows that Method I
ramps power down in the neighborhood of Zmax to avoid over-charging the cell, whereas the PNGV method has no such limits. At very low SOCs, the PNGV method over-predicts how much power is available since there are no current limits applied to the calculation. The Method I estimate is automatically lower due to the large derivative in the denominator of the calculation. This causes an anomaly near zero SOC where the method under-predicts the available charge power. However, since the cell will not be operated in this range, this is not a concern.
Considering now the discharge power curves, the comparison shows that Method I imposes limits on discharge power to ensure that the cell is not under-charged, whereas the PNGV method does not. In the SOC range from about 150 to 35%, the two methods predict similar powers. For SOC
above about 350, the power predicted by Method I saturates because the maximum discharge current limit of 200A has been reached. The PNGV method does not consider this limit.
At SOC around 99 o the graph again shows an anomaly in the Method I calculation where power is under-estimated due to the large derivative term. This apparent glitch is not a problem since the cell will not be operated in this range.
Figs. 9 through 13 show how the two main voltage-limit based methods of power estimation of the present invention (Method I and Method II) compare to the prior art PNGV
method in the dynamic cell tests shown in Fig. 5. Fig. 9 is a plot showing that discharge power capability estimates for cell cycle test comprising sixteen UDDS cycles over an SOC range of 90% down to 100. Fig. 10 is zoomed-in plot of Fig. 9, showing about one UDDS cycle. Fig. 11 is a plot showing charging power capability estimates for cell cycle test comprising sixteen UDDS cycles over an SOC range of 90o down to 100. Fig. 12 is zoomed-in plot of Fig. 11, showing about one UDDS cycle. Again, the absolute value of power is plotted.
In the discussion that follows, the results of Method II
are considered to be the "true" capability of the cell.
This assumption is justified by the fidelity of the cell model's voltage estimates, as supported by the data in Fig.
6. Fig. 9 shows that the three methods produce similar estimates. In particular, Method I and Method II appear to be nearly identical when viewed at this scale. At high SOCs, the PNGV HPPC method predicts higher power than is actually available (by as much as 9.8s), and at low SOCs, the PNGV
HPPC method under-predicts the available power. Only the methods of the present invention include SOC bounds, which explain why their predictions are so different from the PNGV HPPC estimates at low SOC. If the vehicle controller were to discharge at the rates predicted by the PNGV HPPC
method, the cell would be over-discharged in some cases (lowering its lifetime), and under-utilized in other cases.
Fig. 10 zooms in on Fig. 9(same region shown as in Fig. 6) to show greater detail. In this region, the three methods produce nearly identical predictions. A notable feature of Method II, however, is that it takes into account the entire dynamics of the cell when making a prediction.
Therefore, the strong discharges at around time 237 and 267 minutes draw the cell voltage down, and allows less discharge power than the other two methods which only consider SOC when making their estimate.
The three methods are also compared with respect to charge power, shown in Fig. 11. At this scale, the estimates appear nearly identical. Again, the PNGV HPPC
method does not consider SOC limits, so over-predicts charge power at high SOCs. It also over-predicts power at low SOCs as it ignores the increase to charge resistance at low SOC. A zoom of this plot is shown in Fig. 12, which accentuates the differences between the predictions. Here, it can be seen that the strong discharges at around time 237 and 267 minutes allow for greater charging power, as the voltage will not quickly change.
Industrial Applicability While the methods described herein, and the apparatus for carrying these methods into effect, constitute preferred embodiments of the present invention, it should be recognized that changes may be made therein without departing from the spirit or scope of the present invention, which is defined in the appended claims. For example, the steps 10, 12, 14 disclosed in Fig. 1A can be executed in different orders or used in different combinations and steps 20, 22, 24 disclosed in Fig. 1B can be executed in different orders or used in different combinations. Also, various cell models can be substituted for the purpose of estimating the maximum absolute charge/discharge power of a battery/battery cell.
A method and apparatus for calculation of power capability of battery packs using advanced cell model predictive techniques has been described in conjunction with one or more specific embodiments. This invention is defined by the following claims and their full scope ~of equivalents.
One embodiment of the present invention estimates the maximum absolute charge and/or discharge power of a battery pack. The battery pack may be, for example, a battery pack used in a hybrid electric vehicle or an electric vehicle.
The embodiment makes a number of denotations and limits, including:
- using n to denote the number of cells in the target battery pack where an estimation of charge and/or discharge power is desired;
- using vk(t) to denote the cell voltage for cell number k in the pack, which has operational design limits so that k:l<k<_n;
umiri < vk(t) C ~,'maxmust be enforced for all - -- using zk(t) to denote the state-of-charge for cell number k in the pack, which has operational design limits < z t < zmax zmin -_ ~( ) that must be enforced for all ~ ~ 1 < k < n;
- using pk(t) to denote the cell power, which has a operational design limits so that p''''~~n ~ pk(t) ~ ~m~' must be enforced for all ~ ~ 1 ~ 1~ < n; and - using ik(t) to denote the cell current, which has a operational design limits so that 2min ~ 2x (t) C Zmax must be enforced for all ~ ~ 1 < k < n.
Modifications can be made in alternate embodiments. For example, any particular limit may be removed if desired by replacing its value by ~°°, as appropriate. As an another example, 11m1tS SLlCh aS Vmaxr Vminr Zmaxr Zminr lmaxr lmin, pmaxr pmin may furthermore be functions of temperature and other factors pertaining to the present battery pack operating condition. In one embodiment, it is assumed that the discharge current and power have positive sign and the charge current and power have negative sign. Those skilled in the art will recognize that other sign conventions may be used, and that the description of the present invention can be adapted to these conventions in a forthright manner.
In one embodiment, the model used for predicting charge assumes that the battery pack comprises ns cell modules connected in series, where each cell module comprises np individual cells connected in parallel and ns ~ 1, n~ > 1.
Other configurations are possible and are accommodated by slight modifications to the method as described.
Fig. 2 is a schematic block diagram showing the sensor components of an embodiment of the present invention.
Battery 40 is connected to load circuit 48. For example, load circuit 48 could be a motor in an Electric Vehicle (EV) or Hybrid Electric Vehicle (HEV). In some embodiments circuit 48 is a circuit that provides power and/or draws power. Measurements of battery and individual cell voltage are made with voltmeters) 44. Measurements of battery current are made with ammeter 42. Battery and individual cell temperatures are measured by temperature sensors) 46.
Voltage, current and temperature measurements are processed with arithmetic circuit 50. Arithmetic circuit (estimator means) 50, takes in the measurements from the sensor components and perform the calculation methods of the present invention for power estimation. In some embodiment, temperature is not needed in the calculation methods.
1. Calculation Based in State-of-Charge (SOC) Limits As shown in steps 10 and 20 of Figs 1A and 1B, embodiments of the present invention calculate the maximum charge/discharge current values using SOC limits. Various embodiments also have the explicit inclusion of a time horizon Ot in the calculation. The SOC limits are included as follows. First, for a constant current ik, the SOC
recurrent relationship is described as:
zk(t -E- Ot) = zk (t) - (~hOt~C)zk~ ( q ) Where zk ( t) is the present SOC for cell k, zk ( t+0t) is the predicted SOC ~t seconds into the future, C is the cell capacity in ampere-seconds, and ~i is the Coulombic efficiency factor at current level ik. Here, for simplicity of presentation, it is assumed that ~i =1 =or aiscnarge currents and ~~ - ~ < 1 for charge currents .
If there are design limits on SOC such that zmin < zk(t) < zmax for all cells in the pack, then current ik can be computed such that these limits are not exceeded.
Simple algebra gives limits based on the SOC of each cell:
2misa s k - zk (t~ - zcnin chg,soc zk (t~ - zmax Zmin,k - ~Qt~G' The pack maximum absolute currents-based only on cell SOC-are then dis,soc ~ dis,socl Zmax - ~k I1 2m~ Jk chg,soc chg,soc Zcnin - Illk 7C CZmin,k ) .
This method assumes that there is a valid SOC estimate available for every cell in the pack. If there is not the case, then an approximate remedy would be to calculate dis,soc ,~, z(t~ - zmin 2m~ ~ Ot/Cr chg,soc ," z~t) - zmax Zmin '" ~Ot~Cr where z(t) is the pack SOC.
In one embodiment of the present invention, the power predictive method can take into account more information than simply the cell SOC. For example, a Kalman filter can be used as a method to estimate all the cell SOCs in a pack.
Besides giving the SOC, Kalman filtering yields estimates of the uncertainty of the SOC estimate itself. A method of using Kalman filter to estimate SOC is described in commonly assigned U.S. Patent No. 6,534,954, hereby incorporated by reference.
Let the uncertainty have Gaussian distribution with standard deviation, as estimated by the Kalman filter, be denoted as oz. Then, the method yields a95.5o confidence that the true SOC is within the estimate ~2 aZ and a 99.70 confidence that the true SOC is within the estimate ~3 aZ.
This information can be incorporated into the estimate of~ maximum current based on SOC to have very high confidence that SOC design limits will not be violated.
This is done as (assuming a 3 aZ confidence interval):
dis,soc (zk (t) - 3~z) - zmin 2max,k Ot/rr chg,soc (zk(t) + 3Qx) - zmax Zmin,k - ~~t/C ' 2. Calculation Based on Voltage Limits Besides taking SOC limits into account, embodiments of the present invention correct a limitation in the prior art HPPC method for applying voltage limits (steps 12 and 22 of Figs. 1A and 1B). In the HPPC method, if the cell model of equation (1) is assumed, and that Rc''g and Rdls are the cell's Ohmic resistances, then equation (2) and equation (3) predict the instantaneously available current, not the constant value of current that is available for the next Ot seconds. If cases where ~t is large, the result of the calculation poses a safety or battery-health issue, as the cells may become over/under charged.
To overcome this problem, an embodiment of the present invention uses the following cell model:
vk(t + Ot) = OCV(zk(t + Ot)) - R x ik(t), (7 ) This modifies the previous cell model in equation (1).
Note that this model cannot be directly solved in closed form for the maximum current ik since zk ( t+~t) is itself a function of current (cf.(4)) and OCV(~) is a nonlinear relationship. Note that other cell models can be used as well.
Two method embodiments are directed to solving (7) for the maximum absolute value of ik(t).
2.lMethod I: Taylor-series expansion The first method uses a Taylor-series expansion to linearize the equation, so that at approximate value of i can be solved. It is assumed that OCV(~) is differentiable at point zk(t), which gives the result OCV(zk(t + ~1t)) - OCV Czk(t) - ik~C t) - OCV(zk(t)) - ik ~i~t 80CV(z) C ~ 8z I _ t z-xk ( ) + Rl(zk(t), Zk(~i~t)I C) where the first-order residual Rl(zk(t)' Zk(~i0t)~C)~Ilzk(t)II -~
as 2k(~JaOt)~C ~ ~ in 0. Note that the change in SOC over Ot seconds is generally small, so the following approximation may be applied OCV (zk (t + 4t)) = OCV (zk (t)) - ik ~C t a0 ~z (z) x=z t k( ) -Rxik which gives dis,volt OCV(zk(t)) - vm;n 2max,k ~t 80CV x + Rdis C 8z xk (t) ( g ) chg,volt ~CV (zk (t)) - ymax 2min,k - nOt 80CV z ~ cbg C 8x xk~t) + R ( 9 ) In one embodiment, both the function OCV(z) and its derivative aOCV(z)/ 8z might be computed from some known mathematical relationship for OCV(z), (e. g., Nernst's equation) using either analytic or numeric methods, or by a table lookup of empirical data. This quantity is positive for most battery electrochemistries over the entire SOC
range, so the values computed by (8) and (9) are smaller in magnitude than those from ( 2 ) and ( 3 ) for the same values of Rdis and Ring.
The HPPC procedure compensates for its inaccuracy by using modified values of Rdls and R'''g, determined experimentally, that approximate the denominator terms in ( 8 ) and ( 9 ) . This can not be accurate over the entire SOC
range, however, as 80CV(z)/ az is not constant, particularly near extreme values of z.
Discharge and charge currents with all limits enforced are computed as (steps 16 and 26 of Figs. 1A and 1B) dis r dis,soc dis,voltl Zmax - min ~2rny ~lIl Zm~ k ~ IIkIl 2max,k (10) im n - Tna.X C2rnin ~ Ink x im n ~°c ~ IIIkaJC 2m n' k ltJ
(11) and power may be calculated using the sum of all cell powers. These are equal to the product of the maximum allowed current and the predicted future voltage.
'~a chg ~ c6g pmin - nP Zminvk ~t + Ot k=1 na r l N 7LP ~ Zm a ~~iV' zk ~t~ - Zm n~iOt~C~
k= \1 - Rchg X 2m n _ i (12) Tl.a Pmtax - nP ~ Zm'axvk~t -~- ~t~.
k=1 .
na 1 7ZP ~ Zrdc,~ax ~G'v C,Zk(t) - 2m~~Ot~~i~
k= /1 - Rdis X idia max (13) Maximum and minimum cell and pack power limits may also be imposed in this calculation. Note that in all equations, 0CV ( Z) , C, Vmaxr Vminr Zmaxr Zminr lmaxr lmin, Rchgr and Rdls may be functions of temperature and other factors pertaining to the present battery pack operating conditions.
2.2 Method II: Using a Comprehensive Cell Model The method of solving (7) presented in the previous section requires less computational intensity. A second method embodiment of the present invention may be used when more computational power is available. This second method assumes a more precise mathematical model of cell dynamics, which might be in a discrete-time state-space form such as the coupled pair of equations xk[m + 1] - f (~k[m], uk[m]) ( 14 ) vk [m] - 9 ~xk [m]' uk [m]~' ( 15 ) where m is the discrete time sample index; the vector function of time xk[m] is called the estate" of the system, uk[m] is the input to the system, which includes cell current ik[m] as a component, and might also include temperature, resistance, capacity and so forth, and f(~) and g(') are functions chosen to model cell dynamics.
Alternative model forms, including continuous-time state-space forms, differential and difference equations might also be used. It is assumed that there is a method to compute SOC given the model that is implemented.
For convenience of presentation, it is assumed that the cell model is in a discrete-time state-space form. Also assume that 0t seconds may be represented in discrete time as T sample intervals. Then, this model can be used to predict cell voltage 0t seconds into the future by vk(m-t-T~ =g~~k(rn-I-T~,ZGk(7TL+T)), where xk[m+T] may be found by simulating (14) for T time samples. It is assumed that the input remains constant from time index m to m+T, so if temperature change (for example) over this interval is significant, it must be included as part of the dynamics modeled by ( 14 ) and not as a part of the measured input uk[m] .
The method then uses a bisection search algorithm to dis,volt chg,volt find 2max,x and ~min,k by looking for the ik (as a member of the uk vector) that causes equality in vmin - 9 ~xk(m + Z'~ ~ ux(m -I- ~')) ~ or 0 - g ~xk(m + T~, ~k(m -~ T~~ - v~~l~
(16) dis,volt to find 2max,k ~ and by looking for the ik that causes equality in vma,,~ - g~xk(m+T~,uk(m+T~~, or 0 _ g ~xk(m + T), uk(m + T~) - vmax (17) chg,volt to find Zmin,k _ A special case is when the state equation (14) is linear-that is, when xk (m + 1~ = Axk (m~ + Buk (m~
where A and B are constant matrices. The model presented in Section 3, entitled "An Example Cell Model," is an example where this is the case. Then, for input uk constant time m to m+T, leading to xk m + T = AT xk (m~ +~ ~ AT 1 ' B uk.
~=o Most of these terms may be pre-computed without knowledge of uk in order to speed calculation using the bisection algorithm.
dis,soc chg,soc Once the SOC-based current limits 2max,k and zmin,k are computed using (5) and (6), and the voltage-based current dis,volt chg,volt limits Zmax,k and 2min,k are computed using ( 16 ) and ( 17 ) , overall current limits may be computed using ( 10 ) and ( 11 ) (steps 16 and 26 of Figs. 1A and 1B). Power is then computed as na chg ~ chg Pmin - nP 2minvk ~t -~- Ot) k=1 ns chg ( - nP ~ Zmingk~xk(m'+T~~uk)~
k=1 with uk containing 2mn as its value for current, and n' Pmax - nP ~ Zm~~vk \t + Ot k=1 n, dis nP ~ 2max9k ~~k ~7T'G -~- ~~, 2Gk~, k=1 dis with uk containing Zm~ as its value for current.
1.2.1 Bisection search To solve ( 16 ) and ( 17 ) , a method to solve for a root of a nonlinear equation is required. In one embodiment the bisection search algorithm is used for this requirement.
The bisection search algorithm looks for a root of f(x) ( i . e. , a value of x such that f (x) =0) where it is known a priori that the root lies between values x1 < root < x2.
One way of knowing that a root lies in this interval is that the sign of f(xl) is different from the sign of f(xz), Each iteration of the bisection algorithm evaluates the function at the midpoint Xmid = (xi + xz) /2. Based on the sign of the evaluation, either x1 or xz is replaced by Xmid to retain different signs on of f(xl) and f(xz) . It is evident that the uncertainty in the location of the root is halved by this algorithmic step. The bisection algorithm repeats this iteration until the interval between xl and xz, and hence the resolution of the root of f(x) is as small as desired. If E is the desired root resolution, then the algorithm will require at most ~logz ( ~ xl - xz ~ /E) iterations. The bisection method is listed in Listing 1.
1.2.2 Finding maximum/minimum current To determine maximum discharge and charge current for any particular cell, bisection is performed on (16) and (17). Bisection is incorporated in the overall algorithm ad follows. First, three stirnulations are performed to determine cell voltages Ot seconds into the future for cell current ik = 0, ik = lminr and ik = lmax~ If cell voltages are predicted to be between vmi" and Vmax for the maximum discharge rates, then these maximum rates may be used. If the cell voltages, even during rest, are outside of bounds, then set the maximum rates to zero.. Otherwise, the true maximum rate may be found by bisecting between rate equal to zero and its maximum value. Bisection is performed between current limits (lminr 0) or (0, lmax) .
Algorithm 1 The bisection algorithm.
begin{bisection algorithm}
set x1 = first search limit set x2 = second search limit set a = desired resolution on current in bisection output set JMAX = mraximum number of bisection iterations J l iOg2 ( x2 Exl ) let func(~) be function in which to find root set the search interval dx = xa - x1 if ( ( func(xl)>_0) { // constraint: func(xl)<0 dx = -dx x1 = x2 } // root is now between (xl,xa + dx), and func(xl)<0 forj=ltoJMAX{
dx = 0.5 x dx amid = x1 + dx If ( fUnC(xmid) < 0) { x1 = amid }
if (~dx~< e) { return(x~ +0.5 x dx) }
} II have bisected too many times, so return best guess return(xl + 0.5 x dx) end{bisection algorithm}
2. An exam 1e Cell Model An example cell model for the present invention power estimation methods is presented herein, with illustrations given to show the performance of the two methods compared to the prior art PNGV HPPC method. The cell model is a discrete-time state-space model of the form of (14) and (15) that applies to battery cells. The model, named ~~Enhanced Self-Correcting Cell Model," is further described in the article ~~Advances in EKFLiPB SOC Estimation," by the inventor, published in CD-ROM and presented in Proc. 20th Electric Vehicle Symposium (EVS20) in Long Beach CA, (November 2003) and is hereby fully incorporated by reference. It is understood this model is an example model only and that a variety of suitable alternate models can be used.
The ~~Enhanced Self-Correcting Cell Model" includes effects due to open-circuit-voltage, internal resistance, voltage time constants, and hysteresis. For the purpose of example, the 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.
State-of-charge is captured by one state of the model.
This equation is zk(m+ 1~ = z~(m~ - OaOT~C)i,~(m~, where 0 T represents the inter-sample period (in seconds), and C represents the cell capacity (in ampere-seconds).
The time-constants of the cell voltage response are captured by several filter states. If there is let to be of time constants, then fk(m+ 1) = l4ffx(m) + Blik(~).
E 1E8"~ x"p The matrix f may be a diagonal matrix with real-valued entries. If so, the system is stable if all B E lLf~f "1 entries have magnitude less than one. The vector f may simply be set to of ~~l"s. The value of of and the entries in the Af matrix are chosen as part of the system identification procedure to best fit the model parameters to measured cell data.
The hysteresis level is captured by a single state hk (m + 1J - exp ~iZk ~ ~~~T ~~ hk (m]
r~,ix ('mJ'Y~T
+ 1 - exp G' ~~~ S~Uk(m])~
C
C
where Y is the hysteresis rate constant, again found by system identification.
The overall model state is xk (mJ = ( fk ('n-r,J~ hk (rrs) z~ (mJ , where the symbol () is the matrix/vector transpose operator. The state equation for the model is formed by combining all of the individual equations, above. Note, that at each time step, the state equation is linear in the input u~(mJ = ~ ik(m] 1, which speeds the prediction operation.
The output equation that combines the state values to predict cell voltage is vk(mJ = OCV(z~(m~) + Gfk(mJ - Rik(mJ + Mhk(m], where G E ILBIxn~ is a vector of constants that blend the time-constant states together in the output, R is the cell resistance (different values may be used for dis/charge), and M is the maximum hysteresis level.
The open-circuit-voltage as a function of state-of charge for example Lithium Ion Polymer Battery (LiPB) cells is plotted in Fig.3. This is an empirical relationship found by cell testing. First, the cell was fully charged (constant current to 4.2V, constant voltage to 200mA). Then, the cell was discharged at the C/25 rate until fully discharged (3.0V). The cell was then charged at the C/25 rate until the voltage was 4.2V. The low rates were used to minimize the dynamics excited in the cells. The cell voltage as a function of state of charge under discharge and under charge were averaged to compute the OCV. This has the effect of eliminating to the greatest extent possible the presence of hysteresis and ohmic resistance in the final function. For the purpose of computations involving OCV, the final curve was digitized at 200 points and stored in a table. Linear interpolation is used to look up values in the table.
The partial derivative of OCV with respect to SOC for these example cells is plotted in Fig. 4. This relationship was computed by first taking finite differences between points in the OCV plot in Fig. 3 and dividing by the distance between points (i.e., Euler's approximation to a derivative). The resulting data is too noisy to be of practical use, as shown in the gray line of Fig. 4. It was filtered using a zero-phase low-pass filter, resulting in the black line of Fig. 4, which may be used in the power calculation. This relationship is also digitized at 200 points, and linear interpolation into the table of values is used when computations requiring this function are performed.
Other parameters are fit to the cell model. In particular, the model employs four low-pass filter states (nf = 4), a nominal capacity of 7.5 Ah, and an inter-sample interval of 0T = 1s. There is very close agreement between the cell model voltage prediction and the cell true voltage.
This is illustrated in Fig. 5, which is a plot showing the voltage prediction using the cell model of the present invention. For this figure, the cell test was a sequence of sixteen UDDS cycles, performed at room temperature, separated by discharge pulses and five-minute rests, and spread over the 90o to 10o SOC range. The difference between true cell terminal voltage and estimated cell terminal voltage is very small (a root-mean-squared (RMS) voltage estimation error of less than 5mV). To better illustrate the model' s fidelity, refer to the zoom on one UDDS cycle in the 50% SOC region, shown in Fig. 6. The SOC
as a function of time is plotted in Fig. 7, which is a SOC
trace for cell test. The graph shows that SOC increases by about 5% during each UDDS cycle, but is brought down about loo during each discharge between cycles. The entire operating range for these cells (10% SOC to 90o SOC, delineated on the figure as the region between the thin dashed lines) is excited during the cell test.
3. Comparing Maximum Power Calculations The PNGV HPPC power estimation method gives a result that is a function of only SOC. Therefore, it is possible to graph available power versus SOC to summarize the algorithm calculations. The first method proposed (Method I: Taylor Series Expansion Methods) in this patent disclosure is also possible to display in this way.
Estimated power is only a function of SOC, aOCV/8z (also a function of SOC), and static limits on maximum current and power. The second method (Method II: the Comprehensive Cell Model Method), however, dynamically depends on all states of the system. Two systems at the same state of charge, but with different voltage time-constant state values or hysteresis state levels will have different amounts of power available. To compare power computed by the three methods, dynamic tests must be conducted.
For the following result, a pack of LiPB cells is assumed to be with ns - 40 and np - 1. The data to fit models was collected from prototype hand-made cells jointly developed by LG Chem (Daejeon, Korea) and Compact Power Inc.
(Monument, Colorado). Limits for the power calculations are listed in Table 1. Each cell has a nominal capacity of 7.5 Ah, and Ot was ten seconds for both charge and discharge.
Table 1: Parameters for power calculation examples.
Parameter Minimum Maximum vk(t) 3.0V 4.35V
ik (t) -200A 200A
z~(t) 0.1 0.9 pk (t) -o0 00 First, the PNGV HPPC method and Method I of the present invention are compared in Fig. 8, which is a plot comparing static maximum power calculations as functions of SOC for the PNGV HPPC method and Method I of the present invention.
The black curves correspond to charge power, and the gray curves correspond to discharge power. Note that the absolute value of power is plotted to avoid confusion due to sign conventions. Considering first the calculations of charge power, it is evident that the PNGV HPPC method produces similar values to Method I in the mid-SOC range .
The slight differences are due to the fact that the 10-second R~hg value used for the PNGV method and the derivative-modified R~hg for Method I are not identical.
Outside the mid-SOC range, the graph shows that Method I
ramps power down in the neighborhood of Zmax to avoid over-charging the cell, whereas the PNGV method has no such limits. At very low SOCs, the PNGV method over-predicts how much power is available since there are no current limits applied to the calculation. The Method I estimate is automatically lower due to the large derivative in the denominator of the calculation. This causes an anomaly near zero SOC where the method under-predicts the available charge power. However, since the cell will not be operated in this range, this is not a concern.
Considering now the discharge power curves, the comparison shows that Method I imposes limits on discharge power to ensure that the cell is not under-charged, whereas the PNGV method does not. In the SOC range from about 150 to 35%, the two methods predict similar powers. For SOC
above about 350, the power predicted by Method I saturates because the maximum discharge current limit of 200A has been reached. The PNGV method does not consider this limit.
At SOC around 99 o the graph again shows an anomaly in the Method I calculation where power is under-estimated due to the large derivative term. This apparent glitch is not a problem since the cell will not be operated in this range.
Figs. 9 through 13 show how the two main voltage-limit based methods of power estimation of the present invention (Method I and Method II) compare to the prior art PNGV
method in the dynamic cell tests shown in Fig. 5. Fig. 9 is a plot showing that discharge power capability estimates for cell cycle test comprising sixteen UDDS cycles over an SOC range of 90% down to 100. Fig. 10 is zoomed-in plot of Fig. 9, showing about one UDDS cycle. Fig. 11 is a plot showing charging power capability estimates for cell cycle test comprising sixteen UDDS cycles over an SOC range of 90o down to 100. Fig. 12 is zoomed-in plot of Fig. 11, showing about one UDDS cycle. Again, the absolute value of power is plotted.
In the discussion that follows, the results of Method II
are considered to be the "true" capability of the cell.
This assumption is justified by the fidelity of the cell model's voltage estimates, as supported by the data in Fig.
6. Fig. 9 shows that the three methods produce similar estimates. In particular, Method I and Method II appear to be nearly identical when viewed at this scale. At high SOCs, the PNGV HPPC method predicts higher power than is actually available (by as much as 9.8s), and at low SOCs, the PNGV
HPPC method under-predicts the available power. Only the methods of the present invention include SOC bounds, which explain why their predictions are so different from the PNGV HPPC estimates at low SOC. If the vehicle controller were to discharge at the rates predicted by the PNGV HPPC
method, the cell would be over-discharged in some cases (lowering its lifetime), and under-utilized in other cases.
Fig. 10 zooms in on Fig. 9(same region shown as in Fig. 6) to show greater detail. In this region, the three methods produce nearly identical predictions. A notable feature of Method II, however, is that it takes into account the entire dynamics of the cell when making a prediction.
Therefore, the strong discharges at around time 237 and 267 minutes draw the cell voltage down, and allows less discharge power than the other two methods which only consider SOC when making their estimate.
The three methods are also compared with respect to charge power, shown in Fig. 11. At this scale, the estimates appear nearly identical. Again, the PNGV HPPC
method does not consider SOC limits, so over-predicts charge power at high SOCs. It also over-predicts power at low SOCs as it ignores the increase to charge resistance at low SOC. A zoom of this plot is shown in Fig. 12, which accentuates the differences between the predictions. Here, it can be seen that the strong discharges at around time 237 and 267 minutes allow for greater charging power, as the voltage will not quickly change.
Industrial Applicability While the methods described herein, and the apparatus for carrying these methods into effect, constitute preferred embodiments of the present invention, it should be recognized that changes may be made therein without departing from the spirit or scope of the present invention, which is defined in the appended claims. For example, the steps 10, 12, 14 disclosed in Fig. 1A can be executed in different orders or used in different combinations and steps 20, 22, 24 disclosed in Fig. 1B can be executed in different orders or used in different combinations. Also, various cell models can be substituted for the purpose of estimating the maximum absolute charge/discharge power of a battery/battery cell.
A method and apparatus for calculation of power capability of battery packs using advanced cell model predictive techniques has been described in conjunction with one or more specific embodiments. This invention is defined by the following claims and their full scope ~of equivalents.
Claims (13)
1. A method for estimating a maximum discharge power of a battery, comprising:
generating a signal indicative of a present state-of-charge of said battery, utilizing a sensor;
calculating said present state-of-charge of said battery based on said signal, utilizing an arithmetic circuit operably coupled to said sensor;
calculating a maximum discharge current of said battery utilizing said arithmetic circuit based on at least a minimum state-of-charge limit associated with said battery, said present state-of-charge of said battery, and a minimum voltage limit associated with said battery such that a future output voltage of said battery does not fall below said minimum voltage limit and a future state-of-charge of said battery does not fall below said minimum state-of-charge limit associated with said battery; and, calculating said maximum discharge power based on said maximum discharge current value, utilizing said arithmetic circuit.
generating a signal indicative of a present state-of-charge of said battery, utilizing a sensor;
calculating said present state-of-charge of said battery based on said signal, utilizing an arithmetic circuit operably coupled to said sensor;
calculating a maximum discharge current of said battery utilizing said arithmetic circuit based on at least a minimum state-of-charge limit associated with said battery, said present state-of-charge of said battery, and a minimum voltage limit associated with said battery such that a future output voltage of said battery does not fall below said minimum voltage limit and a future state-of-charge of said battery does not fall below said minimum state-of-charge limit associated with said battery; and, calculating said maximum discharge power based on said maximum discharge current value, utilizing said arithmetic circuit.
2. The method of claim 1, wherein said step of calculating said maximum discharge current is also based on a maximum current limit associated with said battery.
3. The method of claim 1, further comprising calculating said present state-of-charge of said battery using a Kalman filtering method.
4. The method of claim 1, wherein said battery is a battery pack comprising at least one cell.
5. A method for estimating a maximum discharge power of a battery, comprising:
generating a signal indicative of a present state-of-charge of said battery, utilizing a sensor;
calculating said present state-of-charge of said battery based on said signal, utilizing an arithmetic circuit operably coupled to said sensor;
calculating a maximum discharge current of said battery utilizing said arithmetic circuit based on at least a minimum state-of-charge limit associated with said battery, said present state-of-charge of said battery, a minimum voltage limit associated with said battery, and a cell model that is solved by a Taylor-series expansion of an open cell voltage of the battery to determine a future output voltage of the battery, such that the future output voltage of said battery does not fall below said minimum voltage limit and a future state-of-charge of said battery does not fall below said minimum state-of-charge limit associated with said battery; and, calculating said maximum discharge power based on said maximum discharge current value, utilizing said arithmetic circuit.
generating a signal indicative of a present state-of-charge of said battery, utilizing a sensor;
calculating said present state-of-charge of said battery based on said signal, utilizing an arithmetic circuit operably coupled to said sensor;
calculating a maximum discharge current of said battery utilizing said arithmetic circuit based on at least a minimum state-of-charge limit associated with said battery, said present state-of-charge of said battery, a minimum voltage limit associated with said battery, and a cell model that is solved by a Taylor-series expansion of an open cell voltage of the battery to determine a future output voltage of the battery, such that the future output voltage of said battery does not fall below said minimum voltage limit and a future state-of-charge of said battery does not fall below said minimum state-of-charge limit associated with said battery; and, calculating said maximum discharge power based on said maximum discharge current value, utilizing said arithmetic circuit.
6. The method of claim 5 wherein said cell model is solved by using a discrete time-state space model.
7. The method of claim 5 wherein said battery is a battery pack comprising at least one cell.
8. The method of claim 7 wherein said cell model is V k(t+.DELTA. t)=OCV(z k(t+.DELTA.t)) - R x i k (t) wherein v k(t+.DELTA.t) denotes a cell voltage for a cell k for a time period .DELTA.t units into the future, OCV(z k(t+ .DELTA.t)) denotes the open cell voltage as a function of a state-of-charge z k for cell k for a time period .DELTA.t units into the future, R is a constant that denotes an internal resistance of said cell k, and i k(t) denotes a cell current for cell k.
9. The method of claim 8 wherein said cell model is solved by using a discrete time-state space model.
10. The method of claim 9 wherein said discrete time-state space model is x k(m+1]=f (x k[m],u k[m]) v k[m] =g(x k[m],u k[m]) wherein m denotes a discrete time sample index, x k[m] denotes a vector function of time and a state of the battery, u k[m] denotes an input to the battery and includes cell current i k[m] as a component, and f(.cndot.) and g(.cndot.) are functions chosen to model cell dynamics.
11. The method of claim 10 wherein is found by looking for i k that causes equality in v min = g(x k[m+T],u k[m+T]) wherein g(x k[m+T], u k[m+T]) is utilized to determine the cell voltage for the cell k at a predetermined time in the future.
12. The method of claim 1, wherein said calculated maximum discharge power is checked to ensure that it falls within power limits of said battery.
13. A system for estimating a maximum discharge power of a battery, comprising:
a sensor configured to generate a signal indicative of a present state-of-charge of said battery; and an arithmetic circuit operably coupled to said sensor, said arithmetic circuit configured to calculate said present state-of-charge of said battery based on said signal, the arithmetic circuit further configured to calculate a maximum discharge current of said battery based on at least a minimum state-of-charge limit associated with said battery, said present state-of-charge of said battery, and a minimum voltage limit associated with said battery such that a future output voltage of said battery does not fall below said minimum voltage limit and a future state-of-charge of said battery does not fall below said minimum state-of-charge limit, said arithmetic circuit further configured to calculate said maximum discharge power based on said maximum discharge current value.
a sensor configured to generate a signal indicative of a present state-of-charge of said battery; and an arithmetic circuit operably coupled to said sensor, said arithmetic circuit configured to calculate said present state-of-charge of said battery based on said signal, the arithmetic circuit further configured to calculate a maximum discharge current of said battery based on at least a minimum state-of-charge limit associated with said battery, said present state-of-charge of said battery, and a minimum voltage limit associated with said battery such that a future output voltage of said battery does not fall below said minimum voltage limit and a future state-of-charge of said battery does not fall below said minimum state-of-charge limit, said arithmetic circuit further configured to calculate said maximum discharge power based on said maximum discharge current value.
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US52432603P | 2003-11-20 | 2003-11-20 | |
US60/524,326 | 2003-11-20 | ||
US10/811,088 | 2004-03-25 | ||
US10/811,088 US7321220B2 (en) | 2003-11-20 | 2004-03-25 | Method for calculating power capability of battery packs using advanced cell model predictive techniques |
PCT/KR2004/003001 WO2005050810A1 (en) | 2003-11-20 | 2004-11-19 | Method for calculating power capability of battery packs using advanced cell model predictive techniques |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2547012A1 CA2547012A1 (en) | 2005-06-02 |
CA2547012C true CA2547012C (en) | 2009-06-02 |
Family
ID=34595116
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002547012A Active CA2547012C (en) | 2003-11-20 | 2004-11-19 | Method for calculating power capability of battery packs using advanced cell model predictive techniques |
Country Status (10)
Country | Link |
---|---|
US (3) | US7321220B2 (en) |
EP (1) | EP1692754B1 (en) |
JP (1) | JP4722857B2 (en) |
KR (1) | KR100894021B1 (en) |
CN (1) | CN1883097B (en) |
BR (1) | BRPI0416652B8 (en) |
CA (1) | CA2547012C (en) |
RU (1) | RU2336618C2 (en) |
TW (1) | TWI281298B (en) |
WO (1) | WO2005050810A1 (en) |
Families Citing this family (195)
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 |
JP2006129588A (en) * | 2004-10-28 | 2006-05-18 | Sanyo Electric Co Ltd | Power control method of secondary battery, and power unit |
US8103485B2 (en) * | 2004-11-11 | 2012-01-24 | Lg Chem, Ltd. | State and parameter estimation for an electrochemical cell |
US7424663B2 (en) * | 2005-01-19 | 2008-09-09 | Intel Corporation | Lowering voltage for cache memory operation |
EP1842275A4 (en) * | 2005-01-25 | 2016-05-11 | Victhom Human Bionics Inc | Power supply charging method and device |
US8355965B2 (en) * | 2005-02-22 | 2013-01-15 | Sharp Kabushiki Kaisha | Battery exchange service system and charging method therefor, and portable device |
WO2006124130A1 (en) * | 2005-03-31 | 2006-11-23 | Energycs | Method and system for retrofitting a full hybrid to be a plug-in hybrid |
DE102005018434A1 (en) * | 2005-04-21 | 2006-10-26 | Continental Aktiengesellschaft | Motor vehicle with a pneumatic level control system |
US7683579B2 (en) * | 2005-06-14 | 2010-03-23 | Lg Chem, Ltd. | Method and apparatus of controlling for charge/discharge power of battery |
US7723957B2 (en) * | 2005-11-30 | 2010-05-25 | Lg Chem, Ltd. | System, method, and article of manufacture for determining an estimated battery parameter vector |
JP4532416B2 (en) * | 2006-01-12 | 2010-08-25 | 古河電気工業株式会社 | Battery discharge capability determination method, battery discharge capability determination device, and power supply system |
JP4773848B2 (en) * | 2006-03-03 | 2011-09-14 | プライムアースEvエナジー株式会社 | Secondary battery charge / discharge control system, battery control device, and program |
US8466684B2 (en) * | 2006-06-16 | 2013-06-18 | Chevron Technology Ventures Llc | Determination of battery predictive power limits |
JP4265629B2 (en) * | 2006-08-01 | 2009-05-20 | トヨタ自動車株式会社 | Secondary battery charge / discharge control device and hybrid vehicle equipped with the same |
TW200824169A (en) * | 2006-11-21 | 2008-06-01 | Benq Corp | Method for predicting remaining capacity of a battery |
US9013139B2 (en) * | 2007-03-26 | 2015-04-21 | The Gillette Company | Adaptive charger device and method |
WO2009059164A2 (en) * | 2007-10-31 | 2009-05-07 | Intrago Corporation | User-distributed shared vehicle system |
JP2009122056A (en) * | 2007-11-19 | 2009-06-04 | Denso Corp | Battery charge/discharge current detection device |
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 |
JP4513882B2 (en) * | 2008-03-21 | 2010-07-28 | トヨタ自動車株式会社 | Hybrid vehicle and control method thereof |
US8067111B2 (en) * | 2008-06-30 | 2011-11-29 | Lg Chem, Ltd. | Battery module having battery cell assembly with heat exchanger |
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 |
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 |
US8539408B1 (en) | 2008-07-29 | 2013-09-17 | Clarkson University | Method for thermal simulation |
JP4722976B2 (en) * | 2008-08-26 | 2011-07-13 | 本田技研工業株式会社 | Storage capacity controller |
US8202645B2 (en) | 2008-10-06 | 2012-06-19 | Lg Chem, Ltd. | Battery cell assembly and method for assembling the battery cell assembly |
TWI398658B (en) * | 2008-10-31 | 2013-06-11 | Evt Technology Co Ltd | An instant calculation method for battery pack information for electric vehicle and electric vehicle using the method |
US8887872B2 (en) * | 2008-11-17 | 2014-11-18 | Otis Elevator Company | Method of determining state of charge of energy storage system |
US8116998B2 (en) | 2009-01-30 | 2012-02-14 | Bae Systems Controls, Inc. | Battery health assessment estimator |
DE102009001300A1 (en) * | 2009-03-03 | 2010-09-09 | Robert Bosch Gmbh | Method and device for determining a characteristic quantity for detecting the stability of the vehicle electrical system |
US8004243B2 (en) * | 2009-04-08 | 2011-08-23 | Tesla Motors, Inc. | Battery capacity estimating method and apparatus |
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 |
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 |
US8852778B2 (en) | 2009-04-30 | 2014-10-07 | Lg Chem, Ltd. | Battery systems, battery modules, and method for cooling a battery module |
US9172118B2 (en) * | 2009-06-17 | 2015-10-27 | Gm Global Technology Operations, Llc. | Method and system for estimating battery life |
FR2947637B1 (en) * | 2009-07-01 | 2012-03-23 | Commissariat Energie Atomique | METHOD FOR CALIBRATING AN ELECTROCHEMICAL BATTERY |
US8703318B2 (en) * | 2009-07-29 | 2014-04-22 | Lg Chem, Ltd. | Battery module and method for cooling the battery module |
US8399118B2 (en) | 2009-07-29 | 2013-03-19 | Lg Chem, Ltd. | Battery module and method for cooling the battery module |
US20120150378A1 (en) * | 2009-08-21 | 2012-06-14 | Reva Electric Car Company (Pvt) Limited | Determination and Usage of Reserve Energy in Stored Energy Systems |
US8399119B2 (en) | 2009-08-28 | 2013-03-19 | Lg Chem, Ltd. | Battery module and method for cooling the battery module |
DE102009042194B4 (en) * | 2009-09-18 | 2019-01-31 | Bayerische Motoren Werke Aktiengesellschaft | Method for determining the operating range of a rechargeable electrical energy store |
DE102009042193A1 (en) | 2009-09-18 | 2011-03-31 | Bayerische Motoren Werke Aktiengesellschaft | Battery's internal resistance and/or impedance values evaluating method for e.g. hybrid vehicle, involves changing parameter values so that actual values match with evaluated values when evaluated values are deviated from actual values |
DE102009042192B4 (en) | 2009-09-18 | 2023-12-14 | Bayerische Motoren Werke Aktiengesellschaft | Method for estimating the internal resistance or impedance value of a battery |
DE102009049320A1 (en) | 2009-10-14 | 2011-04-21 | Bayerische Motoren Werke Aktiengesellschaft | Method for determining and / or predicting the high-current capacity of a battery |
DE102009049589A1 (en) * | 2009-10-16 | 2011-04-21 | Bayerische Motoren Werke Aktiengesellschaft | Method for determining and / or predicting the maximum performance of a battery |
DE102009045783A1 (en) * | 2009-10-19 | 2011-04-21 | Robert Bosch Gmbh | Precise power prediction method for battery packs |
FR2952235B1 (en) * | 2009-10-29 | 2015-01-16 | Commissariat Energie Atomique | METHOD FOR CHARGING OR DISCHARGING A BATTERY TO DETERMINE THE END OF CHARGE OR DISCHARGE BASED ON CURRENT MEASUREMENTS AND TEMPERATURE |
JP5496612B2 (en) * | 2009-11-11 | 2014-05-21 | 三洋電機株式会社 | Battery chargeable / dischargeable current calculation method, power supply device, and vehicle equipped with the same |
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 |
JP5517692B2 (en) * | 2010-03-26 | 2014-06-11 | 三菱重工業株式会社 | Battery pack and battery control system |
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 |
EP2579059B1 (en) * | 2010-06-07 | 2014-04-02 | Mitsubishi Electric Corporation | Charge status estimation apparatus |
JP2012039725A (en) * | 2010-08-05 | 2012-02-23 | Toyota Motor Corp | Charging method and charging system |
US8353315B2 (en) | 2010-08-23 | 2013-01-15 | Lg Chem, Ltd. | End cap |
US8758922B2 (en) | 2010-08-23 | 2014-06-24 | Lg Chem, Ltd. | Battery system and manifold assembly with two manifold members removably coupled together |
US8469404B2 (en) | 2010-08-23 | 2013-06-25 | Lg Chem, Ltd. | Connecting assembly |
US8920956B2 (en) | 2010-08-23 | 2014-12-30 | Lg Chem, Ltd. | Battery system and manifold assembly having a manifold member and a connecting fitting |
US9005799B2 (en) | 2010-08-25 | 2015-04-14 | Lg Chem, Ltd. | Battery module and methods for bonding cell terminals of battery cells together |
FR2964196B1 (en) * | 2010-08-25 | 2012-08-24 | Peugeot Citroen Automobiles Sa | DEVICE FOR DETERMINING THE SIZE (S) OF A SERIES CELL BATTERY |
US8662153B2 (en) | 2010-10-04 | 2014-03-04 | Lg Chem, Ltd. | Battery cell assembly, heat exchanger, and method for manufacturing the heat exchanger |
US8749201B2 (en) * | 2010-10-05 | 2014-06-10 | GM Global Technology Operations LLC | Battery pack capacity learn algorithm |
DE102010038017B4 (en) | 2010-10-06 | 2019-03-14 | FEV Europe GmbH | Method for determining the maximum possible current flow of a battery by a control unit and control unit |
CN106451600B (en) * | 2010-10-22 | 2019-06-18 | 核科学股份有限公司 | For the device and method of quickly charging battery |
US8400112B2 (en) * | 2010-11-10 | 2013-03-19 | Ford Global Technologies, Llc | Method for managing power limits for a battery |
WO2012070190A1 (en) * | 2010-11-25 | 2012-05-31 | パナソニック株式会社 | Charging control circuit, battery driven apparatus, charging apparatus, and charging method |
TWI428622B (en) | 2010-11-25 | 2014-03-01 | Ind Tech Res Inst | Method for checking and modulating battery capacity and power based on battery charging/discharging characteristics |
CN102195105B (en) * | 2011-01-14 | 2013-03-20 | 中兴通讯股份有限公司 | Method and apparatus for controlling accumulator charging |
CN102590749B (en) * | 2011-01-14 | 2016-09-28 | 朴昌浩 | A kind of design method for cell charge state prediction kernel |
WO2012098523A2 (en) * | 2011-01-19 | 2012-07-26 | Sendyne Corp. | Converging algorithm for real-time battery prediction |
CN102137536B (en) * | 2011-03-28 | 2014-06-04 | 天津英诺华微电子技术有限公司 | Adaptive control method for brightness of street lamp using solar energy for power supply |
US8288031B1 (en) | 2011-03-28 | 2012-10-16 | Lg Chem, Ltd. | Battery disconnect unit and method of assembling the battery disconnect unit |
DE102011007884A1 (en) * | 2011-04-21 | 2012-10-25 | Sb Limotive Company Ltd. | Method for determining a maximum available constant current of a battery |
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 |
FR2975501B1 (en) * | 2011-05-20 | 2013-05-31 | Renault Sas | METHOD FOR ESTIMATING THE CHARGE STATE OF AN ELECTRIC BATTERY |
US20120316810A1 (en) * | 2011-06-08 | 2012-12-13 | GM Global Technology Operations LLC | Battery limit calibration based on battery life and performance optimization |
US10234512B2 (en) * | 2011-06-11 | 2019-03-19 | Sendyne Corporation | Current-based cell modeling |
US8859119B2 (en) | 2011-06-30 | 2014-10-14 | 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 |
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 |
US9260033B2 (en) * | 2011-07-13 | 2016-02-16 | Sanyo Electric Co., Ltd. | Power supply device and vehicle including the same |
US9496544B2 (en) | 2011-07-28 | 2016-11-15 | Lg Chem. Ltd. | Battery modules having interconnect members with vibration dampening portions |
CN103033752B (en) * | 2011-09-30 | 2016-01-20 | 吴昌旭 | Battery of electric vehicle life-span prediction method and prolongation method |
CN103094630B (en) * | 2011-10-28 | 2015-04-15 | 东莞钜威新能源有限公司 | Battery management method and system thereof |
CN102419599B (en) * | 2011-10-31 | 2013-08-07 | 江苏科技大学 | Artificial fish swarm algorithm-based solar battery maximal power point tracking method |
JP2013115846A (en) * | 2011-11-25 | 2013-06-10 | Denso Corp | Guard processing device for battery pack |
US8718850B2 (en) * | 2011-11-30 | 2014-05-06 | Nec Laboratories America, Inc. | Systems and methods for using electric vehicles as mobile energy storage |
US8977510B2 (en) | 2011-12-15 | 2015-03-10 | Lg Chem, Ltd. | System and method for determining charging and discharging power levels for a battery pack |
CN103185865A (en) * | 2011-12-31 | 2013-07-03 | 陕西汽车集团有限责任公司 | Real-time estimation method of SOC (stress optical coefficient) closed loop of electric automobile lithium ion battery by EKF (extended kalman filter) |
US20130175996A1 (en) * | 2012-01-05 | 2013-07-11 | IDesign, Inc | Characterizing battery discharge under different loads |
DE102012200414A1 (en) * | 2012-01-12 | 2013-07-18 | Siemens Aktiengesellschaft | Method and device for determining a state of charge of an electrical energy store |
GB2499052A (en) * | 2012-02-01 | 2013-08-07 | Continental Automotive Systems | Calculating a power value in a vehicular application |
DE102012202077A1 (en) * | 2012-02-13 | 2013-08-14 | Robert Bosch Gmbh | Method for determining a current, battery management unit, battery and motor vehicle |
CN102590680A (en) * | 2012-02-29 | 2012-07-18 | 广东步步高电子工业有限公司 | Intelligent power supply capable of simulating characteristics of true battery |
DE102012204957A1 (en) * | 2012-03-28 | 2013-10-02 | Robert Bosch Gmbh | Method for determining a maximum available constant current of a battery, arrangement for carrying out such a method, battery in combination with such an arrangement and motor vehicle with such a battery |
CN102608542B (en) * | 2012-04-10 | 2013-12-11 | 吉林大学 | Method for estimating charge state of power cell |
US8922217B2 (en) * | 2012-05-08 | 2014-12-30 | GM Global Technology Operations LLC | Battery state-of-charge observer |
US9067598B2 (en) * | 2012-06-14 | 2015-06-30 | GM Global Technology Operations LLC | Method and apparatus for controlling a high-voltage electrical system for a multi-mode transmission |
KR101405354B1 (en) | 2012-11-13 | 2014-06-11 | 카코뉴에너지 주식회사 | Method for PV cell modeling |
CN103901344A (en) * | 2012-12-24 | 2014-07-02 | 财团法人金属工业研究发展中心 | Battery residual power estimating system and estimating method thereof |
DE102013000572A1 (en) | 2013-01-15 | 2014-07-17 | Rheinisch-Westfälische Technische Hochschule Aachen | Method for determining model parameters of electrochemical energy storage of e.g. electric vehicle, involves defining parameter record variant as new reference dataset to describe battery model and to determine maximum power of storage |
JP6071725B2 (en) | 2013-04-23 | 2017-02-01 | カルソニックカンセイ株式会社 | Driving force control device for electric vehicles |
JP5878495B2 (en) * | 2013-06-11 | 2016-03-08 | 株式会社豊田中央研究所 | Electric vehicle power supply system |
CN104298793B (en) * | 2013-07-16 | 2017-11-21 | 万向一二三股份公司 | A kind of model of power battery pack power limit is counter to push away dynamic algorithm |
CN103401534B (en) * | 2013-07-24 | 2015-12-23 | 中达电通股份有限公司 | Based on the method for the final voltage fast prediction of RC circuit discharge and recharge |
US9368841B2 (en) | 2013-08-30 | 2016-06-14 | Ford Global Technologies, Llc | Battery power capability estimation at vehicle start |
US10473723B2 (en) * | 2013-08-30 | 2019-11-12 | Ford Global Technologies, Llc | Parameter and state limiting in model based battery control |
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 |
CN104698382A (en) * | 2013-12-04 | 2015-06-10 | 东莞钜威新能源有限公司 | Method for predicting the SOC and SOH of battery pack |
FR3016091B1 (en) * | 2013-12-27 | 2016-01-08 | Renault Sas | METHOD AND SYSTEM FOR BATTERY MANAGEMENT FOR MOTOR VEHICLE |
CN103995232B (en) * | 2014-04-21 | 2017-01-04 | 中通客车控股股份有限公司 | A kind of detection method of lithium iron phosphate dynamic battery group peak value charge-discharge performance |
GB201407805D0 (en) | 2014-05-02 | 2014-06-18 | Dukosi Ltd | Battery condition determination |
US10830821B2 (en) * | 2014-05-05 | 2020-11-10 | Apple Inc. | Methods and apparatus for battery power and energy availability prediction |
US9419314B2 (en) * | 2014-05-12 | 2016-08-16 | GM Global Technology Operations LLC | Systems and methods for determining battery system power capability |
CN103995235A (en) * | 2014-06-02 | 2014-08-20 | 耿直 | Output current limiting method for storage battery modules in high-voltage storage battery pack |
US20160001672A1 (en) * | 2014-07-01 | 2016-01-07 | Ford Global Technologies, Llc | Equivalent circuit based battery current limit estimations |
US10451678B2 (en) | 2014-07-17 | 2019-10-22 | Ford Global Technologies, Llc | Battery system identification through impulse injection |
US9381823B2 (en) * | 2014-07-17 | 2016-07-05 | Ford Global Technologies, Llc | Real-time battery estimation |
JP6324248B2 (en) | 2014-07-17 | 2018-05-16 | 日立オートモティブシステムズ株式会社 | Battery state detection device, secondary battery system, battery state detection program, battery state detection method |
JP6316690B2 (en) | 2014-07-17 | 2018-04-25 | 日立オートモティブシステムズ株式会社 | Battery state detection device, secondary battery system, battery state detection program, battery state detection method |
US9321368B2 (en) | 2014-08-19 | 2016-04-26 | Fca Us Llc | Multi-step model predictive iterative techniques for battery system peak power estimation |
US10408880B2 (en) | 2014-08-19 | 2019-09-10 | Fca Us Llc | Techniques for robust battery state estimation |
EP3017993B1 (en) | 2014-11-07 | 2021-04-21 | Volvo Car Corporation | Power and current estimation for batteries |
CN104391251B (en) * | 2014-11-18 | 2017-04-26 | 郑州日产汽车有限公司 | Data acquisition method of electric vehicle battery management system |
FR3029298B1 (en) * | 2014-11-28 | 2016-12-30 | Renault Sa | AUTOMATIC METHOD OF ESTIMATING THE CHARGING STATE OF A CELL OF A BATTERY |
FR3029296B1 (en) | 2014-11-28 | 2016-12-30 | Renault Sa | AUTOMATIC METHOD OF ESTIMATING THE CHARGING STATE OF A CELL OF A BATTERY |
CN105891715A (en) * | 2014-12-12 | 2016-08-24 | 广西大学 | Lithium ion battery health state estimation method |
CN104483539B (en) * | 2015-01-08 | 2017-02-22 | 湖南大学 | Active power rapid measuring method based on Taylor expansion |
FI3045925T3 (en) | 2015-01-14 | 2023-05-29 | Corvus Energy Ltd | Method and system for iteratively determining state of charge of a battery cell |
CN104537268B (en) * | 2015-01-19 | 2018-08-21 | 重庆长安汽车股份有限公司 | A kind of battery maximum discharge power evaluation method and device |
CN104635165B (en) * | 2015-01-27 | 2017-03-29 | 合肥工业大学 | A kind of accurate evaluation method of light-electricity complementary electric power system accumulator dump energy |
WO2016157731A1 (en) * | 2015-03-27 | 2016-10-06 | パナソニックIpマネジメント株式会社 | Secondary battery status estimation device and status estimation method |
US9789784B2 (en) * | 2015-05-13 | 2017-10-17 | Ford Global Technologies, Llc | Maintaining a vehicle battery |
US10048320B2 (en) * | 2015-06-12 | 2018-08-14 | GM Global Technology Operations LLC | Systems and methods for estimating battery system power capability |
CN105277895B (en) * | 2015-09-30 | 2018-01-26 | 上海凌翼动力科技有限公司 | A kind of series battery power rating SOP On-line Estimation method and its application |
JP6787660B2 (en) * | 2015-12-10 | 2020-11-18 | ビークルエナジージャパン株式会社 | Battery control device, power system |
US10298042B2 (en) * | 2016-02-05 | 2019-05-21 | Nec Corporation | Resilient battery charging strategies to reduce battery degradation and self-discharging |
CN105699910A (en) * | 2016-04-21 | 2016-06-22 | 中国计量大学 | Method for on-line estimating residual electric quantity of lithium battery |
US9921272B2 (en) * | 2016-05-23 | 2018-03-20 | Lg Chem, Ltd. | System for determining a discharge power limit value and a charge power limit value of a battery cell |
EP3252616A1 (en) * | 2016-06-01 | 2017-12-06 | Honeywell spol s.r.o. | Estimating variables that are not directly measurable |
CN106054085B (en) * | 2016-07-11 | 2019-11-15 | 四川普力科技有限公司 | A method of based on temperature for estimating battery SOC |
CN106443459A (en) * | 2016-09-06 | 2017-02-22 | 中国第汽车股份有限公司 | Evaluation method of state of charge of vehicle lithium ion power battery |
CN106526490A (en) * | 2016-10-25 | 2017-03-22 | 宁德时代新能源科技股份有限公司 | Limit working current obtaining method and device |
DE102016224376B4 (en) * | 2016-12-07 | 2018-10-18 | Bayerische Motoren Werke Aktiengesellschaft | Method and control unit for operating a stationary memory |
CN106842038B (en) * | 2016-12-14 | 2020-08-14 | 广东恒沃动力科技有限公司 | Online estimation method for maximum discharge power of battery |
TWI597511B (en) * | 2016-12-23 | 2017-09-01 | Chen Tech Electric Mfg Co Ltd | Battery life cycle prediction model building method |
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 |
WO2018162023A2 (en) * | 2017-03-06 | 2018-09-13 | Volvo Truck Corporation | A battery state of power estimation method and a battery state monitoring system |
EP3593155B1 (en) * | 2017-03-06 | 2021-04-07 | Volvo Truck Corporation | A battery cell state of charge estimation method and a battery state monitoring system |
JP6729460B2 (en) * | 2017-03-17 | 2020-07-22 | トヨタ自動車株式会社 | In-vehicle battery charge controller |
CN107102271A (en) * | 2017-05-25 | 2017-08-29 | 宁德时代新能源科技股份有限公司 | Estimation method, device and system for peak power of battery pack |
CN107402356B (en) * | 2017-08-04 | 2020-03-20 | 南京南瑞继保电气有限公司 | EKF estimation lead-acid battery SOC method based on dynamic parameter identification |
US10237830B1 (en) | 2017-08-31 | 2019-03-19 | Google Llc | Dynamic battery power management |
KR101946784B1 (en) * | 2017-09-29 | 2019-02-12 | 한국과학기술원 | method for measuring battery entropy by dual kalman filter |
KR102515606B1 (en) * | 2017-10-31 | 2023-03-28 | 삼성에스디아이 주식회사 | Method, battery pack, and electronic device for displaying battery charge capacity |
KR102203245B1 (en) | 2017-11-01 | 2021-01-13 | 주식회사 엘지화학 | Apparatus and method for estimating SOC of battery |
US10549649B2 (en) * | 2017-11-10 | 2020-02-04 | GM Global Technology Operations LLC | Maximum current calculation and power prediction for a battery pack |
CN107861075B (en) * | 2017-12-24 | 2020-03-27 | 江西优特汽车技术有限公司 | Method for determining SOP of power battery |
DE102018200144A1 (en) | 2018-01-08 | 2019-07-11 | Robert Bosch Gmbh | Method and management system for controlling and monitoring multiple battery cells of a battery pack and battery pack |
CN108549746B (en) * | 2018-03-26 | 2021-11-26 | 浙江零跑科技股份有限公司 | Battery system power limit estimation method based on cell voltage |
CN108427079B (en) * | 2018-06-11 | 2019-06-07 | 西南交通大学 | A kind of power battery method for estimating remaining capacity |
CN109343370B (en) * | 2018-11-30 | 2021-10-01 | 北京宇航系统工程研究所 | Dynamic environment simulation system and method for space power supply controller |
KR102465889B1 (en) | 2018-12-18 | 2022-11-09 | 주식회사 엘지에너지솔루션 | Apparatus and Method for Controlling Charging of Secondary Battery Pack |
KR20200112248A (en) | 2019-03-21 | 2020-10-05 | 주식회사 엘지화학 | Apparatus and method for controlling battery bank |
JP6722954B1 (en) * | 2019-04-02 | 2020-07-15 | 東洋システム株式会社 | Battery residual value determination system |
CN110009528B (en) * | 2019-04-12 | 2021-06-01 | 杭州电子科技大学 | Parameter self-adaptive updating method based on optimal structure multi-dimensional Taylor network |
US11485239B2 (en) * | 2019-06-03 | 2022-11-01 | GM Global Technology Operations LLC | Power prediction for reconfigurable series-connected battery with mixed battery chemistry |
DE102019214407A1 (en) | 2019-09-20 | 2021-03-25 | Robert Bosch Gmbh | Method for determining a first voltage characteristic curve of a first electrical energy storage unit |
US11515587B2 (en) * | 2019-10-10 | 2022-11-29 | Robert Bosch Gmbh | Physics-based control of battery temperature |
CN111025172B (en) * | 2019-12-31 | 2022-03-01 | 国联汽车动力电池研究院有限责任公司 | Method for realizing rapid measurement of maximum allowable power of charging and discharging of lithium ion battery |
US11313912B2 (en) * | 2020-01-28 | 2022-04-26 | Karma Automotive Llc | Battery power limits estimation based on RC model |
US11454673B2 (en) * | 2020-02-12 | 2022-09-27 | Karma Automotive Llc | Battery current limits estimation based on RC model |
CN111289927A (en) * | 2020-03-02 | 2020-06-16 | 大陆汽车电子(长春)有限公司 | Starting signal simulation device, test method and system of intelligent battery sensor |
WO2021254620A1 (en) * | 2020-06-18 | 2021-12-23 | Volvo Truck Corporation | A method for predicting state-of-power of a multi-battery electric energy storage system |
US11522440B2 (en) * | 2020-07-29 | 2022-12-06 | Cirrus Logic, Inc. | Use of shared feedback among two or more reactive schemes |
CN112763918A (en) * | 2020-12-23 | 2021-05-07 | 重庆金康动力新能源有限公司 | Electric vehicle power battery discharge power control method |
KR20220100331A (en) * | 2021-01-08 | 2022-07-15 | 주식회사 엘지에너지솔루션 | Method for predicting maximum discharge current and battery system using the same |
CN113109716B (en) * | 2021-04-06 | 2022-08-23 | 江苏大学 | Lithium battery SOP estimation method based on electrochemical model |
WO2022241301A1 (en) * | 2021-05-14 | 2022-11-17 | Carnegie Mellon University | System and method for application-dependent selection of batteries with differentiable programming |
FR3123156A1 (en) * | 2021-05-18 | 2022-11-25 | Psa Automobiles Sa | BATTERY MANAGEMENT SYSTEM COMPRISING MEANS FOR DETECTING LIMITING CELLS, VEHICLE AND METHOD BASED ON SUCH A SYSTEM |
CN113466697B (en) * | 2021-06-10 | 2024-02-27 | 深圳拓邦股份有限公司 | SOC estimation method for battery, computer terminal and storage medium |
KR20230054191A (en) * | 2021-10-15 | 2023-04-24 | 주식회사 엘지에너지솔루션 | Battery control system for guiding soc level and method of the same |
CN114643892A (en) * | 2022-04-11 | 2022-06-21 | 广州万城万充新能源科技有限公司 | Electric vehicle charging power prediction system based on multi-mode data perception |
CN116070466B (en) * | 2023-03-08 | 2023-06-13 | 上海泰矽微电子有限公司 | Optimal SOC simulation value searching method under battery cut-off voltage |
Family Cites Families (94)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4390841A (en) * | 1980-10-14 | 1983-06-28 | Purdue Research Foundation | Monitoring apparatus and method for battery power supply |
DE69025868T2 (en) * | 1989-12-11 | 1996-09-05 | Canon Kk | charger |
CN1017561B (en) * | 1990-10-12 | 1992-07-22 | 湖南轻工研究所 | Method for predicting intermittent discharge capacity of dry battery |
JP3209457B2 (en) * | 1992-12-11 | 2001-09-17 | 本田技研工業株式会社 | Battery remaining capacity detection method |
US5825155A (en) * | 1993-08-09 | 1998-10-20 | Kabushiki Kaisha Toshiba | Battery set structure and charge/ discharge control apparatus for lithium-ion battery |
US5714866A (en) * | 1994-09-08 | 1998-02-03 | National Semiconductor Corporation | Method and apparatus for fast battery charging using neural network fuzzy logic based control |
US5578915A (en) * | 1994-09-26 | 1996-11-26 | General Motors Corporation | Dynamic battery state-of-charge and capacity determination |
US5606242A (en) * | 1994-10-04 | 1997-02-25 | Duracell, Inc. | Smart battery algorithm for reporting battery parameters to an external device |
US5633573A (en) * | 1994-11-10 | 1997-05-27 | Duracell, Inc. | Battery pack having a processor controlled battery operating system |
JPH08140209A (en) * | 1994-11-11 | 1996-05-31 | Fuji Heavy Ind Ltd | Battery managing system for electric motor vehicle |
TW269727B (en) * | 1995-04-03 | 1996-02-01 | Electrosource Inc | Battery management system |
FR2740554A1 (en) | 1995-10-31 | 1997-04-30 | Philips Electronique Lab | SYSTEM FOR MONITORING THE DISCHARGE PHASE OF THE CHARGING-DISCHARGE CYCLES OF A RECHARGEABLE BATTERY, AND HOST DEVICE PROVIDED WITH AN INTELLIGENT BATTERY |
JP3520886B2 (en) | 1996-03-08 | 2004-04-19 | サンケン電気株式会社 | Rechargeable battery status determination method |
US5694335A (en) * | 1996-03-12 | 1997-12-02 | Hollenberg; Dennis D. | Secure personal applications network |
JP3385845B2 (en) * | 1996-05-22 | 2003-03-10 | 日産自動車株式会社 | Electric vehicle power control device |
US6064180A (en) * | 1996-10-29 | 2000-05-16 | General Motors Corporation | Method and apparatus for determining battery state-of-charge using neural network architecture |
US5739670A (en) * | 1996-10-31 | 1998-04-14 | General Motors Corporation | Method for diagnosing battery condition |
JP2000504477A (en) * | 1996-11-21 | 2000-04-11 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Battery management system and battery simulator |
JPH113505A (en) * | 1997-06-11 | 1999-01-06 | Sony Corp | Magnetic head and its manufacture |
JPH1123676A (en) * | 1997-06-30 | 1999-01-29 | Sony Corp | Method and device for measuring charge characteristics of secondary cell |
JPH1132442A (en) | 1997-07-10 | 1999-02-02 | Matsushita Electric Ind Co Ltd | Method and apparatus for estimating residual capacity of storage battery |
JPH1138105A (en) * | 1997-07-15 | 1999-02-12 | Toyota Autom Loom Works Ltd | Method for calculating residual capacity of battery and method for outputting alarm to insufficient residual capacity |
US6285163B1 (en) * | 1998-05-28 | 2001-09-04 | Toyota Jidosha Kabushiki Kaisha | Means for estimating charged state of battery and method for estimating degraded state of battery |
US6018227A (en) * | 1998-06-22 | 2000-01-25 | Stryker Corporation | Battery charger especially useful with sterilizable, rechargeable battery packs |
KR100271094B1 (en) * | 1998-08-21 | 2000-11-01 | 김덕중 | Charge contoller |
US6353815B1 (en) * | 1998-11-04 | 2002-03-05 | The United States Of America As Represented By The United States Department Of Energy | Statistically qualified neuro-analytic failure detection method and system |
US6232744B1 (en) * | 1999-02-24 | 2001-05-15 | Denso Corporation | Method of controlling battery condition of self-generation electric vehicle |
EP1206826B1 (en) | 1999-05-05 | 2009-02-25 | Midtronics, Inc. | Energy management system for automotive vehicle |
JP2001095160A (en) * | 1999-09-17 | 2001-04-06 | Matsushita Electric Ind Co Ltd | Method of detecting abnormal battery cell |
DE19959019A1 (en) | 1999-12-08 | 2001-06-13 | Bosch Gmbh Robert | Method for status detection of an energy store |
DE19960761C1 (en) * | 1999-12-16 | 2001-05-23 | Daimler Chrysler Ag | Battery residual charge monitoring method uses difference between minimum current and current at intersection point between current/voltage curve and minimum voltage threshold of battery |
DE10021161A1 (en) * | 2000-04-29 | 2001-10-31 | Vb Autobatterie Gmbh | Method for determining the state of charge and the load capacity of an electric accumulator |
TW535308B (en) * | 2000-05-23 | 2003-06-01 | Canon Kk | Detecting method for detecting internal state of a rechargeable battery, detecting device for practicing said detecting method, and instrument provided with said |
US6583606B2 (en) * | 2000-05-29 | 2003-06-24 | Matsushita Electric Industrial Co., Ltd. | Method for charging determining life of and detecting state of charge of battery |
JP2002145543A (en) | 2000-11-09 | 2002-05-22 | Mitsubishi Electric Corp | Control device of elevator |
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 |
US6407532B1 (en) * | 2000-12-29 | 2002-06-18 | Nokia Mobile Phones, Ltd. | Method and apparatus for measuring battery charge and discharge current |
JP2002228730A (en) | 2001-02-06 | 2002-08-14 | Shikoku Electric Power Co Inc | Residual electric energy estimating device for secondary battery |
DE10106508A1 (en) * | 2001-02-13 | 2002-08-29 | Bosch Gmbh Robert | Method and arrangement for determining the performance of a battery |
DE10106505A1 (en) * | 2001-02-13 | 2002-08-29 | Bosch Gmbh Robert | Method and device for condition detection of technical systems such as energy storage |
JP4292721B2 (en) * | 2001-02-14 | 2009-07-08 | 株式会社日本自動車部品総合研究所 | Battery state control method for hybrid vehicle |
DE10207659B4 (en) * | 2001-02-23 | 2006-09-28 | Yazaki Corp. | A method and apparatus for estimating a terminal voltage of a battery, method and apparatus for calculating an open circuit voltage of a battery, and method and apparatus for calculating battery capacity |
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 |
CN1387279A (en) * | 2001-05-21 | 2002-12-25 | 孕龙科技股份有限公司 | Charge-discharge energy indicator and measuring method |
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 |
EP1417503B1 (en) * | 2001-06-29 | 2010-02-17 | Robert Bosch Gmbh | Methods for determining the charge state and/or the power capacity of a charge store |
US20030015993A1 (en) * | 2001-07-17 | 2003-01-23 | Sudhan Misra | Battery charging system with electronic logbook |
US7072871B1 (en) * | 2001-08-22 | 2006-07-04 | Cadex Electronics Inc. | Fuzzy logic method and apparatus for battery state of health determination |
JP3672248B2 (en) * | 2001-09-19 | 2005-07-20 | インターナショナル・ビジネス・マシーンズ・コーポレーション | Electrical apparatus, computer apparatus, intelligent battery, battery diagnosis method, battery state display method, and program |
US6727708B1 (en) * | 2001-12-06 | 2004-04-27 | Johnson Controls Technology Company | Battery monitoring system |
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 |
US20030184307A1 (en) * | 2002-02-19 | 2003-10-02 | Kozlowski James D. | Model-based predictive diagnostic tool for primary and secondary batteries |
JP4038788B2 (en) | 2002-02-22 | 2008-01-30 | アクソンデータマシン株式会社 | Battery remaining capacity determination method and apparatus |
JP2003257501A (en) * | 2002-02-27 | 2003-09-12 | Suzuki Motor Corp | Secondary battery residual capacity meter |
AU2003224494A1 (en) * | 2002-03-19 | 2003-09-29 | Mapper Lithography Ip B.V. | Direct write lithography system |
JP2004031014A (en) * | 2002-06-24 | 2004-01-29 | Nissan Motor Co Ltd | Maximum charge-discharge power calculation method and device for battery pack including parallel connection battery |
KR100471233B1 (en) * | 2002-06-26 | 2005-03-10 | 현대자동차주식회사 | Method of generating maximum charge current and maximum discharge current for battery in a hybrid electric vehicle |
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 |
JP3539424B2 (en) * | 2002-07-24 | 2004-07-07 | 日産自動車株式会社 | Electric vehicle control 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 |
US6832171B2 (en) * | 2002-12-29 | 2004-12-14 | Texas Instruments Incorporated | Circuit and method for determining battery impedance increase with aging |
US6892148B2 (en) * | 2002-12-29 | 2005-05-10 | Texas Instruments Incorporated | Circuit and method for measurement of battery capacity fade |
US7518339B2 (en) * | 2003-01-30 | 2009-04-14 | 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 |
US7199557B2 (en) * | 2003-07-01 | 2007-04-03 | Eaton Power Quality Company | Apparatus, methods and computer program products for estimation of battery reserve life using adaptively modified state of health indicator-based reserve life models |
DE10335930B4 (en) * | 2003-08-06 | 2007-08-16 | Vb Autobatterie Gmbh & Co. Kgaa | Method for determining the state of an electrochemical storage battery |
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 |
TWI251951B (en) * | 2003-09-19 | 2006-03-21 | Lg Chemical Ltd | Nonaqueous lithium secondary battery with cyclability and/or high temperature safety improved |
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 |
US20050127874A1 (en) * | 2003-12-12 | 2005-06-16 | Myoungho Lim | Method and apparatus for multiple battery cell management |
KR100651573B1 (en) * | 2003-12-18 | 2006-11-29 | 주식회사 엘지화학 | Apparatus and method for testing state of charge in battery using neural network |
JP4583765B2 (en) * | 2004-01-14 | 2010-11-17 | 富士重工業株式会社 | Remaining capacity calculation device for power storage device |
DE102004005478B4 (en) * | 2004-02-04 | 2010-01-21 | Vb Autobatterie Gmbh | Method for determining parameters for electrical states of a storage battery and monitoring device for this purpose |
US7126312B2 (en) * | 2004-07-28 | 2006-10-24 | Enerdel, Inc. | Method and apparatus for balancing multi-cell lithium battery systems |
US7233128B2 (en) * | 2004-07-30 | 2007-06-19 | Ford Global Technologies, Llc | Calculation of state of charge offset using a closed integral method |
US7525285B2 (en) * | 2004-11-11 | 2009-04-28 | Lg Chem, Ltd. | Method and system for cell equalization using state of charge |
US8103485B2 (en) * | 2004-11-11 | 2012-01-24 | Lg Chem, Ltd. | State and parameter estimation for an electrochemical cell |
US7593821B2 (en) * | 2004-11-23 | 2009-09-22 | Lg Chem, Ltd. | Method and system for joint battery state and parameter estimation |
US7315789B2 (en) * | 2004-11-23 | 2008-01-01 | Lg Chem, Ltd. | Method and system for battery parameter estimation |
US7424663B2 (en) * | 2005-01-19 | 2008-09-09 | Intel Corporation | Lowering voltage for cache memory operation |
US7197487B2 (en) * | 2005-03-16 | 2007-03-27 | Lg Chem, Ltd. | Apparatus and method for estimating battery state of charge |
US7589532B2 (en) * | 2005-08-23 | 2009-09-15 | Lg Chem, Ltd. | System and method for estimating a state vector associated with a battery |
US7446504B2 (en) * | 2005-11-10 | 2008-11-04 | Lg Chem, Ltd. | System, method, and article of manufacture for determining an estimated battery state vector |
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 |
US7521895B2 (en) * | 2006-03-02 | 2009-04-21 | Lg Chem, Ltd. | System and method for determining both an estimated battery state vector and an estimated battery parameter vector |
-
2004
- 2004-03-25 US US10/811,088 patent/US7321220B2/en active Active
- 2004-11-18 TW TW093135430A patent/TWI281298B/en active
- 2004-11-19 WO PCT/KR2004/003001 patent/WO2005050810A1/en active Application Filing
- 2004-11-19 RU RU2006121548/09A patent/RU2336618C2/en active
- 2004-11-19 KR KR1020067009587A patent/KR100894021B1/en active IP Right Grant
- 2004-11-19 BR BRPI0416652A patent/BRPI0416652B8/en active IP Right Grant
- 2004-11-19 CN CN2004800344009A patent/CN1883097B/en active Active
- 2004-11-19 JP JP2006541030A patent/JP4722857B2/en active Active
- 2004-11-19 CA CA002547012A patent/CA2547012C/en active Active
- 2004-11-19 EP EP04800108.5A patent/EP1692754B1/en active Active
-
2007
- 2007-12-21 US US11/963,307 patent/US7656122B2/en active Active
-
2009
- 2009-12-30 US US12/649,844 patent/US7969120B2/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
EP1692754B1 (en) | 2018-01-10 |
US7969120B2 (en) | 2011-06-28 |
US7321220B2 (en) | 2008-01-22 |
JP4722857B2 (en) | 2011-07-13 |
WO2005050810A1 (en) | 2005-06-02 |
KR100894021B1 (en) | 2009-04-17 |
RU2336618C2 (en) | 2008-10-20 |
EP1692754A4 (en) | 2016-03-16 |
BRPI0416652B1 (en) | 2017-11-21 |
TWI281298B (en) | 2007-05-11 |
CN1883097B (en) | 2011-11-16 |
CA2547012A1 (en) | 2005-06-02 |
US20100174500A1 (en) | 2010-07-08 |
RU2006121548A (en) | 2008-01-20 |
EP1692754A1 (en) | 2006-08-23 |
CN1883097A (en) | 2006-12-20 |
US7656122B2 (en) | 2010-02-02 |
KR20060107535A (en) | 2006-10-13 |
US20080094035A1 (en) | 2008-04-24 |
US20050110498A1 (en) | 2005-05-26 |
BRPI0416652A (en) | 2007-01-16 |
BRPI0416652B8 (en) | 2023-01-17 |
JP2007517190A (en) | 2007-06-28 |
TW200531399A (en) | 2005-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2547012C (en) | Method for calculating power capability of battery packs using advanced cell model predictive techniques | |
Plett | High-performance battery-pack power estimation using a dynamic cell model | |
Farmann et al. | Comparative study of reduced order equivalent circuit models for on-board state-of-available-power prediction of lithium-ion batteries in electric vehicles | |
Farmann et al. | A comprehensive review of on-board State-of-Available-Power prediction techniques for lithium-ion batteries in electric vehicles | |
Dong et al. | Kalman filter for onboard state of charge estimation and peak power capability analysis of lithium-ion batteries | |
Sun et al. | Model-based dynamic multi-parameter method for peak power estimation of lithium–ion batteries | |
Sun et al. | Adaptive unscented Kalman filtering for state of charge estimation of a lithium-ion battery for electric vehicles | |
Roscher et al. | Reliable state estimation of multicell lithium-ion battery systems | |
Kwak et al. | Parameter identification and SOC estimation of a battery under the hysteresis effect | |
Zhang et al. | Fractional calculus based modeling of open circuit voltage of lithium-ion batteries for electric vehicles | |
CN106772081B (en) | Battery limit charging and discharging current estimation method based on extended equivalent circuit model | |
CN107748336A (en) | The state-of-charge On-line Estimation method and system of lithium ion battery | |
CN109683101A (en) | A method of battery remaining power is obtained based on SOC-OCV curve | |
Snoussi et al. | Unknown input observer design for lithium-ion batteries SOC estimation based on a differential-algebraic model | |
EP4096962A1 (en) | Battery power limits estimation based on rc model | |
Mu et al. | Modeling, evaluation, and state estimation for batteries | |
Biswas et al. | Simultaneous state and parameter estimation of li-ion battery with one state hysteresis model using augmented unscented kalman filter | |
Speltino et al. | Experimental identification and validation of an electrochemical model of a lithium-ion battery | |
Aung et al. | State-of-charge estimation using particle swarm optimization with inverse barrier constraint in a nanosatellite | |
Santos et al. | Lead acid battery SoC estimation based on extended Kalman Filter method considering different temperature conditions | |
Naik et al. | State of Charge Estimation of Lithium-ion Batteries for Electric Vehicle. | |
Mitra et al. | Ukf based estimation of soc and core temperature of a lithium ion cell using an electrical cell model | |
Rahighi et al. | Developing a battery monitoring system software in matlab simulink environment using kalman filter | |
Rabab et al. | Equivalent Circuit Model For Sodium-Ion Batteries With Physical-Based Representations Of Their Non-Linearities | |
Chen et al. | Improving the accuracy of voltage estimation in the low charge state range at low temperature: An equivalent circuit model considering the influence of temperature on solid phase diffusion process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request |